| 1 |
# ==================== |
| 2 |
# | Model parameters | |
| 3 |
# ==================== |
| 4 |
# |
| 5 |
# Continuous equation parameters |
| 6 |
&PARM01 |
| 7 |
#- mean theta* from H&S, isothermal (T=200.K) above p=150.mb: |
| 8 |
tRef= 295.1, 295.2, 295.4, 295.6, 296.0, 296.6, 297.3, 298.2, 299.3, 300.6, |
| 9 |
302.4, 305.9, 311.8, 321.8, 338.3, 359.5, 382.8, 408.6, 437.5, 470.6, |
| 10 |
509.2, 555.8, 614.9, 696.3, 827.2, 1153.6, |
| 11 |
sRef=26*0., |
| 12 |
viscAr=0.E1, |
| 13 |
viscAh=0.E6, |
| 14 |
viscA4=0.E17, |
| 15 |
no_slip_sides=.FALSE., |
| 16 |
no_slip_bottom=.FALSE., |
| 17 |
diffKhT=0.E3, |
| 18 |
diffKrT=0., |
| 19 |
diffK4T=0.E17, |
| 20 |
diffKrS=0.E2, |
| 21 |
diffKhS=0.E3, |
| 22 |
diffK4S=0.E17, |
| 23 |
buoyancyRelation='ATMOSPHERIC', |
| 24 |
eosType='IDEALG', |
| 25 |
atm_Cp=1004.64, |
| 26 |
atm_Rq=0.6078, |
| 27 |
integr_GeoPot=2, |
| 28 |
selectFindRoSurf=1, |
| 29 |
gravity=9.80, |
| 30 |
implicitFreeSurface=.TRUE., |
| 31 |
exactConserv=.TRUE., |
| 32 |
select_rStar=2, |
| 33 |
nonlinFreeSurf=4, |
| 34 |
hFacInf=0.2, |
| 35 |
hFacSup=2.0, |
| 36 |
hFacMin=1.0, |
| 37 |
uniformLin_PhiSurf=.FALSE., |
| 38 |
vectorInvariantMomentum=.TRUE., |
| 39 |
useAbsVorticity=.TRUE., |
| 40 |
selectVortScheme=3, |
| 41 |
selectKEscheme=3, |
| 42 |
addFrictionHeating=.TRUE., |
| 43 |
#tempAdvScheme=77, |
| 44 |
saltAdvScheme=77, |
| 45 |
staggerTimeStep=.TRUE., |
| 46 |
readBinaryPrec=64, |
| 47 |
writeBinaryPrec=64, |
| 48 |
#globalFiles=.TRUE., |
| 49 |
useSingleCpuIO=.TRUE., |
| 50 |
plotLevel=0, |
| 51 |
& |
| 52 |
|
| 53 |
# Elliptic solver parameters |
| 54 |
&PARM02 |
| 55 |
cg2dMaxIters=200, |
| 56 |
cg2dTargetResWunit=1.E-17, |
| 57 |
& |
| 58 |
|
| 59 |
# Time stepping parameters |
| 60 |
&PARM03 |
| 61 |
nIter0=81000, |
| 62 |
nTimeSteps=10, |
| 63 |
#endTime=311040000., |
| 64 |
deltaT=384., |
| 65 |
abEps=0.1, |
| 66 |
forcing_In_AB=.FALSE., |
| 67 |
cAdjFreq=0., |
| 68 |
#- long run (> 1.yr): |
| 69 |
chkptFreq =15552000., |
| 70 |
pChkptFreq=31104000., |
| 71 |
dumpFreq =15552000., |
| 72 |
monitorFreq= 864000., |
| 73 |
#- short run (~ months): |
| 74 |
#chkptFreq = 864000., |
| 75 |
#pChkptFreq=2592000., |
| 76 |
#dumpFreq = 432000., |
| 77 |
#monitorFreq= 43200., |
| 78 |
#- short test: |
| 79 |
#dumpFreq =1., |
| 80 |
monitorFreq=1., |
| 81 |
& |
| 82 |
|
| 83 |
# Gridding parameters |
| 84 |
&PARM04 |
| 85 |
usingCurvilinearGrid=.TRUE., |
| 86 |
horizGridFile='dxC1_dXYa', |
| 87 |
radius_fromHorizGrid=6370.E3, |
| 88 |
delR= 1500., 2122., 3000., 4242., 6000., 7800., 9100., 9400., 9280., 8400., |
| 89 |
6980., 5800., 4820., 4004., 3328., 2766., 2298., 1910., 1586., 1318., |
| 90 |
1096., 910., 756., 628., 522., 434., |
| 91 |
& |
| 92 |
|
| 93 |
# Input datasets |
| 94 |
&PARM05 |
| 95 |
hydrogThetaFile='ini_theta_26l.bin', |
| 96 |
hydrogSaltFile ='ini_specQ_26l.bin', |
| 97 |
& |