[Mitgcm-support] RE: Momentum timestep

mitgcm-support at dev.mitgcm.org mitgcm-support at dev.mitgcm.org
Wed Jul 9 15:47:06 EDT 2003


Do we have to introduce the useCD flag? If so why can't it default to
true, which is closer to how it is now (with no flag)? 

Chris



> -----Original Message-----
> From: Jean-Michel Campin [mailto:jmc at gulf.mit.edu] 
> Sent: Wednesday, April 16, 2003 11:41 AM
> To: support at mitgcm.org
> Subject: Momentum timestep
> 
> 
> Hi,
> 
> I'am going to check-in some changes in the momentum timestepping,
> 1) store u*,v* in gU,V instead of in gu,vNm1 (tracers timestepping is 
>    already done in the same way). This affects several S/R
>    but has no impact on output.
> 2) move the forcing and the CD scheme from mom_fluxform & mom_vecinv 
>    to timestep.F 
> 
> The motivations for (2) are:
>  a) to allow to put the momForcing out of the Adams-Bashforth.
>   (since it seems logical to keep the AB in timestep.F, then 
> the forcing
>    has to follows).
>  b) CD scheme needs to be called after the forcing, so it has 
> to move also.
>    Anyway, uVelD & vVelD are step-forward in mom_cdscheme, 
> and seems logical 
>    to have also this call in timestep.F.
>    Doing like this allows to use the CD scheme with mom_vecinv.
>  c) this open the way for doing the AB on u,v instead of gU,gV.
>   [this improve the energy conservation and was actually used for the 
>    NLFS-conservation paper]
> 
> (2) affects the outputs (at the truncation level) since the sequence 
>   of addition to gU,gV (+Forcing , +Coriolis, + MetricTerm 
> ... ) is modified.
>   (will need to update the output.txt of the high sensitive test exp,
>    glob-ocean-pressure, lab_sea and nat_box)
> 
> Also I propose to add a running flag parameter for the CD-scheme,
>  (useCDscheme) default value =FALSE 
>  This means that when using the CD scheme, 
> "useCDscheme=.TRUE." will have 
>  to be added to the "data" file.
>  (I plan to add "WARNING" when useCDscheme=.TRUE. is missing)
> 
> Any comments ?
> 
> See you,
> 
> Jean-Michel
> 




More information about the MITgcm-support mailing list