[MITgcm-support] advection scheme
Julie Deshayes
jdeshayes at whoi.edu
Mon Aug 27 15:48:58 EDT 2007
hello,
I am running a configuration of the model at high latitude (f=1e-4 /s).
my domain is 700km x 700km and the resolution is 10 km. I do not have as
many eddies as I expected, so I am running some tests to decrease the
values of viscAh and diffKhT that I am using now (see data file at the
end of my email). I hope that this will give me more eddies. Besides, I
am thinking of changing the advection scheme (currently centered 4th
order) to a non-linear flux limited scheme. Would you recommend me to do
that ? or, considering my configuration, what advection scheme would you
recommend me to use ?
thank you in advance for your help,
Julie
# Model parameters
#
# Continuous equation parameters
&PARM01
readBinaryPrec=64,
useSingleCpuIO=.TRUE.,
tRef= 5.35, 5.20, 5.05, 4.9, 4.75, 4.6, 4.45, 4.3, 4.15, 4.0,
sRef=10*10.,
rigidLid=.FALSE.,
implicitFreeSurface=.TRUE.,
hFacMin=0.2,
no_slip_sides=.TRUE.,
no_slip_bottom=.TRUE.,
bottomDragLinear=1.e-6,
eosType='LINEAR',
tAlpha=2.E-4,
sBeta =0.,
saltStepping=.FALSE.,
tempAdvScheme=4,
saltAdvScheme=4,
diffKhT=120.,
diffKzT=1.E-5,
f0=1.0e-4,
beta=0.E-11,
# useCDscheme=.FALSE.,
useNHMTerms=.FALSE.,
implicitDiffusion=.TRUE.,
ivdc_kappa=1000.0,
viscAz=2.E-4,
viscAh=200.,
&
#
# Elliptic solver parameters
&PARM02
cg2dMaxIters=1000,
cg2dTargetResidual=1.E-11,
&
#
# Time stepping parameters
&PARM03
startTime=622080000.,
endTime=933120000.,
deltaT=1500.0,
abEps=0.1,
pChkptFreq=0.0,
chkptFreq=311040000.,
dumpFreq=311040000.,
# taveFreq=86400.,
cAdjFreq=0,
&
More information about the MITgcm-support
mailing list