[MITgcm-support] set initial flow field

Martin Losch martin.losch at awi.de
Tue Aug 30 04:19:12 EDT 2016


Hi Xumingquan

the parameter names are uVelInitFile and vVelInitfile, w will be computed from div(u) = 0. Be aware that specifying initial velocities with inconsistent hydrography can lead to strong initial shocks and may require very short time steps (initially).
Here are some verification experiments that use these flags

verification> grep -i uvelinitfile */input*/data | grep -v \#
advect_xz/input.nlfs/data: uVelInitFile='Udiv.bin',
advect_xz/input.pqm/data: uVelInitFile='Uvel.bin',
advect_xz/input/data: uVelInitFile='Uvel.bin',
fizhi-cs-32x32x40/input/data: uVelInitFile='uVel_cs2-32x40.d19930110z06.bin',
fizhi-cs-aqualev20/input/data: uVelInitFile='U.input.20Lev',
front_relax/input.bvp/data: uVelInitFile   ='u_ini.bin',
front_relax/input.mxl/data: uVelInitFile   ='u_ini.bin',
lab_sea/input.hb87/data: uVelInitFile    = 'U.0000008760',
obcs_ctrl/input_ad/data: uVelInitFile=   'Uini.bin',
offline_exf_seaice/input.dyn_jfnk/data: uVelInitFile    = 'uVel_3c0.bin',
offline_exf_seaice/input.dyn_lsr/data: uVelInitFile    = 'uVel_3c0.bin',
offline_exf_seaice/input.thermo/data: uVelInitFile    = 'const+20.bin',
offline_exf_seaice/input.thsice/data: uVelInitFile    = 'const+20.bin',
offline_exf_seaice/input/data: uVelInitFile    = 'uVel_3c0.bin',
offline_exf_seaice/input_ad.thsice/data: uVelInitFile    = 'const+20.bin',
offline_exf_seaice/input_ad/data: uVelInitFile    = 'const+20.bin',
seaice_itd/input/data: uVelInitFile    = 'uVel_3c0.bin',
so_box_biogeo/input/data: uVelInitFile  = 'U_ini.bin',
tutorial_deep_convection/input.smag3d/data: uVelInitFile =  'U.120mn.bin',
tutorial_deep_convection/input/data: uVelInitFile =  'U.120mn.bin',


Martin

> On 29 Aug 2016, at 17:53, elisabei <elisabei at aliyun.com> wrote:
> 
> Dear:
>       The initial salinity field and initial temperature field can been loaded by setting the parameter "hydrogThetaFile"&"hydrogSaltFile" in data.So the  initial flow field(u,v w) can been loaded?How can I do?Have some examples in code of MITgcm?
> 
> Best
> 
> xumingquan
> 2016.8.29
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list