[MITgcm-devel] Inverted echo sounders

[Matthew Mazloff]

Hi Patrick,

I want to constrain to daily averaged travel time for inverted echo  
sounders -- so I put
# define ALLOW_IESTAU_COST_CONTRIBUTION
in ECCO_CPPOPTIONS.h

with this defined cost_averagesfields calls a new routine:
cost_averages_ies.F

Here T,S, and P are used to get sound speed.  Sound speed is then  
integrated up to get daily average travel time and this is saved as  
file iestaubar

The files I had to make minor additions too are:
ctrl_dummy.h   (declared xx_iestaubar_mean_dummy)
cost_averages_bar_directives.h  (CADJ STORE xx_iestaubar_mean_dummy =  
onetape)
ecco_cost.h  (declared iestaubar, iestaubarfile)
cost_averagesinit.F  (iestaubar, xx_iestaubar_mean_dummy)
ecco_cost_init_barfiles.F (fnameiestaubar,adfnameiestaubar)
and I added cost_averages_ies.f to ecco_ad_diff.list
cost_averagesfields.F  (now calls cost_averages_ies)

The code compiles but gives major recomputations.  See below

Is there any place obvious I need to add a store?

Thanks
Matt










 From taf_ad.log:

TAF RECOMPUTATION LOOP WARNING DOLOOP_STMT ad_input_code.f:441434 in  
the_main_loop
  extensive recomputations are required.
  The responsible variables are :  
psbar 
,sstbar,wbar,xx_psbar_mean_dummy,xx_sstbar_mean_dummy,xx_wbar_mean_dummy
TAF RECOMPUTATION WARNING DOLOOP_STMT ad_input_code.f:441434 in  
the_main_loop
  triggered by      :  
apressure0,apressure1,aqh0,aqh1,atemp0,atemp1,empmr0,empmr1,etan,evap,gs,gsnm1,gt,gtnm1,gunm1,gvnm1,hflux0,hflux1,lwdown0,lwdown1,lwflux0,lwflux1,obns,obns0,obns1,obnt,obnt0,obnt1,obnu0,obnu1,obnv0,obnv1,obss,obss0,obss1,obst,obst0,obst1,obsu0,obsu1,obsv0,obsv1,obws,obws0,obws1,obwt,obwt0,obwt1,obwu0,obwu1,obwv0,obwv1,pload0,pload1,precip0,precip1,psbar,qnet0,qnet1,qsw0,qsw1,runoff0,runoff1,salt,saltflux0,saltflux1,sflux0,sflux1,siceload,sss0,sss1,sst0,sst1,sstbar,surfaceforcingtice,swdown0,swdown1,swflux0,swflux1,taux0,taux1,tauy0,tauy1,theta,totphihyd,ustress0,ustress1,uvel,uwind0,uwind1,vstress0,vstress1,vvel,vwind0,vwind1,wbar,wspeed0,wspeed1,wvel,xx_aqh0,xx_aqh1,xx_atemp0,xx_atemp1,xx_obcsn0,xx_obcsn1,xx_obcss0,xx_obcss1,xx_obcsw0,xx_obcsw1,xx_precip0,xx_precip1,xx_psbar_mean_dummy,xx_sstbar_mean_dummy,xx_swdown0,xx_swdown1,xx_uwind0,xx_uwind1,xx_vwind0,xx_vwind1,xx_wbar_mean_dummy
  still to be recomputed :  
apressure0,apressure1,aqh0,aqh1,atemp0,atemp1,empmr0,empmr1,etan,evap,gs,gsnm1,gt,gtnm1,gunm1,gvnm1,hflux0,hflux1,lwdown0,lwdown1,lwflux0,lwflux1,obns,obns0,obns1,obnt,obnt0,obnt1,obnu0,obnu1,obnv0,obnv1,obss,obss0,obss1,obst,obst0,obst1,obsu0,obsu1,obsv0,obsv1,obws,obws0,obws1,obwt,obwt0,obwt1,obwu0,obwu1,obwv0,obwv1,pload0,pload1,precip0,precip1,psbar,qnet0,qnet1,qsw0,qsw1,runoff0,runoff1,salt,saltflux0,saltflux1,sflux0,sflux1,siceload,sss0,sss1,sst0,sst1,sstbar,surfaceforcingtice,swdown0,swdown1,swflux0,swflux1,taux0,taux1,tauy0,tauy1,theta,totphihyd,ustress0,ustress1,uvel,uwind0,uwind1,vstress0,vstress1,vvel,vwind0,vwind1,wbar,wspeed0,wspeed1,wvel,xx_aqh0,xx_aqh1,xx_atemp0,xx_atemp1,xx_obcsn0,xx_obcsn1,xx_obcss0,xx_obcss1,xx_obcsw0,xx_obcsw1,xx_precip0,xx_precip1,xx_psbar_mean_dummy,xx_sstbar_mean_dummy,xx_swdown0,xx_swdown1,xx_uwind0,xx_uwind1,xx_vwind0,xx_vwind1,xx_wbar_mean_dummy
TAF RECOMPUTATION WARNING SEQU_STMT ad_input_code.f:441433 in  
the_main_loop
  triggered by      :  
apressure0,apressure1,aqh0,aqh1,atemp0,atemp1,empmr0,empmr1,etan,evap,gs,gsnm1,gt,gtnm1,gunm1,gvnm1,hflux0,hflux1,lwdown0,lwdown1,lwflux0,lwflux1,obns,obns0,obns1,obnt,obnt0,obnt1,obnu0,obnu1,obnv0,obnv1,obss,obss0,obss1,obst,obst0,obst1,obsu0,obsu1,obsv0,obsv1,obws,obws0,obws1,obwt,obwt0,obwt1,obwu0,obwu1,obwv0,obwv1,pload0,pload1,precip0,precip1,psbar,qnet0,qnet1,qsw0,qsw1,runoff0,runoff1,salt,saltflux0,saltflux1,sflux0,sflux1,siceload,sss0,sss1,sst0,sst1,sstbar,surfaceforcingtice,swdown0,swdown1,swflux0,swflux1,taux0,taux1,tauy0,tauy1,theta,totphihyd,ustress0,ustress1,uvel,uwind0,uwind1,vstress0,vstress1,vvel,vwind0,vwind1,wbar,wspeed0,wspeed1,wvel,xx_aqh0,xx_aqh1,xx_atemp0,xx_atemp1,xx_obcsn0,xx_obcsn1,xx_obcss0,xx_obcss1,xx_obcsw0,xx_obcsw1,xx_precip0,xx_precip1,xx_psbar_mean_dummy,xx_sstbar_mean_dummy,xx_swdown0,xx_swdown1,xx_uwind0,xx_uwind1,xx_vwind0,xx_vwind1,xx_wbar_mean_dummy
  still to be recomputed :  
psbar 
,sstbar,wbar,xx_psbar_mean_dummy,xx_sstbar_mean_dummy,xx_wbar_mean_dummy
TAF INFO    ad_input_code.f:441474 : tape computations inserted
TAF RECOMPUTATION LOOP WARNING DOLOOP_STMT ad_input_code.f:441434 in  
the_main_loop
  extensive recomputations are required.
  The responsible variables are :  
psbar 
,sstbar,wbar,xx_psbar_mean_dummy,xx_sstbar_mean_dummy,xx_wbar_mean_dummy
*WARNING* : default new identifier ilev_21                          
changed to ilev_22
TAF RECOMPUTATION WARNING DOLOOP_STMT ad_input_code.f:441434 in  
the_main_loop
.
.
.
.