[MITgcm-devel] 2nd set of calls in initialise_varia.F

[Jean-Michel Campin]

Hi Patrick,

Going to check-in the changes in initialise_varia.F.
I agree on the "remaining issue". 
Will try to look for some explanation (but not sure I will find something)

Cheers,
Jean-Michel

On Mon, Jun 25, 2012 at 03:29:11PM -0400, Patrick Heimbach wrote:
> Hi Jean-Michel,
> 
> I've test all three cases, with some small changes all work
> (I tested all with removing suggested calls in initialise_varia.F,
> so seems fine to remove those).
> 
> Remaining issue is that we can't get same results when running
> bottom_ctrl with #define NONLIN_FRSURF (and keeping nonlinFreeSurf=0)
> or instead with #undef NONLIN_FRSURF
> 
> Cheers
> -Patrick
> 
> On Jun 17, 2012, at 10:38 AM, Jean-Michel Campin wrote:
> 
> > Hi Patrick,
> > 
> > It looks like I can, in initialise_varia.F, remove the calls to:
> > INI_HFAC , INI_MASKS_ETC , INI_LINEAR_PHISURF (all 3 are within ifdef ALLOW_DEPTH_CONTROL)
> > and also INI_CG2D (and INI_CG3D) (these 2 are outside ifdef ALLOW_DEPTH_CONTROL)
> > and the results remain identical with or without these 4(+1) calls (and same
> > number of recomputation warnings).
> > 
> > I tried this in 3 Adjoint set-up:
> > 1) standard bottom_ctrl
> > 2) bottom_ctrl with #define NONLIN_FRSURF (and keeping nonlinFreeSurf=0)
> > 3) tutorial_tracer_adjsens but compiled with ALLOW_CG2D_NSA
> >    and run with nonlinFreeSurf=3 (and keeping select_rStar=2)
> > 
> > I am tempting to remove these 4(+1) calls in initialise_varia.F
> > what is your opinion on this ?
> > 
> > Cheers,
> > Jean-Michel
> > 
> > 
> > _______________________________________________
> > MITgcm-devel mailing list
> > MITgcm-devel at mitgcm.org
> > http://mitgcm.org/mailman/listinfo/mitgcm-devel
> 
> ---
> Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
> MIT | EAPS 54-1420 | 77 Massachusetts Ave | Cambridge MA 02139 USA
> FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach
> 
> 



> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel