[MITgcm-devel] strange recomputation in global_ocean.cs32x15

[Jean-Michel Campin]

Hi Patrick,

In experiment global_ocean.cs32x15, I tried to generate the adjoint using the new
calling sequence in seaice+thsice (commenting out #define OLD_THSICE_CALL_SEQUENCE
 in global_ocean.cs32x15/code_ad/CPP_OPTIONS.h )
but TAF gives me 2 strange recomputation warnings like this one:

> TAF RECOMPUTATION LOOP WARNING DOLOOP_STMT ad_input_code.f:205824 in thsice_do_advect
>  extensive recomputations are required.
>  The responsible variables are : _comlev1_bibj_icefrc_123,_comlev1_bibj_iceheight_111,_comlev1_bibj_iceheight_116,_comlev1_bibj_iceheight_121,_comlev1_bibj_icemask_117,_comlev1_bibj_icevol_113,_comlev1_bibj_qice2_118,_comlev1_bibj_snowheight_112,_comlev1_bibj_snowheight_122,_comlev1_bibj_utrice_114,_comlev1_bibj_utrice_119,_comlev1_bibj_vtrice_115,_comlev1_bibj_vtrice_120,_comlev1_thsice_4_af_11,_comlev1_thsice_4_af_6,_comlev1_thsice_4_icefld_12,_comlev1_thsice_4_icefld_4,_comlev1_thsice_4_icefld_7,_comlev1_thsice_4_icefld_9,_comlev1_thsice_4_icevol_10,_comlev1_thsice_4_icevol_13,_comlev1_thsice_4_icevol_5,_comlev1_thsice_4_icevol_8,bi,bj,hicemin,hocemxl,iceheight,icemask,icemaskmin,ikey_dynamics,imax,imin,jmax,jmin,mybxhi,mybxlo,mybyhi,mybylo,myiter,mythid,ntx,nty,qice1,qice2,qsnow,rhoi,rhos,saltice,snowheight,thsice_deltat,thsice_diffk,thsiceadvscheme,usecubedsphereexchange,usediagnostics,useseaice,ygrid_rl15,ygrid_rl7,ygrid_rs18,ygrid_rs22,ygrid_rs3,ygrid_rs32,ygrid_rs44,ygrid_rs7,yseaice_dynvars_14,yseaice_dynvars_15,yw2_exch2_topo_i21,yw2_exch2_topo_i22,yw2_exch2_topo_i23,yw2_exch2_topo_i24,yw2_exch2_topo_i8,yw2_map_tile2proc2

I don't know what to do with this (but would be good to be able to 
remove the OLD_THSICE_CALL_SEQUENCE code, and to use the new calling sequence
code). If you find some time to look at this.

Cheers,
Jean-Michel