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 |
#tRef=26*280., |
12 |
sRef=26*0., |
13 |
viscAr=0.E1, |
14 |
viscAh=0.E6, |
15 |
viscA4=0.E17, |
16 |
no_slip_sides=.FALSE., |
17 |
no_slip_bottom=.FALSE., |
18 |
diffKhT=0.E3, |
19 |
diffKrT=0., |
20 |
diffK4T=0.E17, |
21 |
diffKrS=0.E2, |
22 |
diffKhS=0.E3, |
23 |
diffK4S=0.E17, |
24 |
buoyancyRelation='ATMOSPHERIC', |
25 |
eosType='IDEALG', |
26 |
atm_Cp=1004.64, |
27 |
atm_Rq=0.6078, |
28 |
integr_GeoPot=2, |
29 |
selectFindRoSurf=1, |
30 |
gravity=9.80, |
31 |
rigidLid=.FALSE., |
32 |
implicitFreeSurface=.TRUE., |
33 |
exactConserv=.TRUE., |
34 |
select_rStar=2, |
35 |
nonlinFreeSurf=4, |
36 |
hFacInf=0.2, |
37 |
hFacSup=2.0, |
38 |
hFacMin=1.0, |
39 |
uniformLin_PhiSurf=.FALSE., |
40 |
vectorInvariantMomentum=.TRUE., |
41 |
useAbsVorticity=.TRUE., |
42 |
selectVortScheme=3, |
43 |
selectKEscheme=3, |
44 |
addFrictionHeating=.TRUE., |
45 |
#tempAdvScheme=77, |
46 |
saltAdvScheme=77, |
47 |
staggerTimeStep=.TRUE., |
48 |
readBinaryPrec=64, |
49 |
writeBinaryPrec=64, |
50 |
#globalFiles=.TRUE., |
51 |
#useSingleCpuIO=.TRUE., |
52 |
debugLevel=2, |
53 |
plotLevel=0, |
54 |
& |
55 |
|
56 |
# Elliptic solver parameters |
57 |
&PARM02 |
58 |
cg2dMaxIters=200, |
59 |
cg2dTargetResWunit=1.E-17, |
60 |
& |
61 |
|
62 |
# Time stepping parameters |
63 |
&PARM03 |
64 |
nIter0=0, |
65 |
# 10.yrs or 8.yrs: |
66 |
#nTimeSteps=864000, |
67 |
nTimeSteps=691200, |
68 |
#endTime=311040000., |
69 |
#endTime=5184000., |
70 |
deltaT=360., |
71 |
abEps=0.1, |
72 |
forcing_In_AB=.FALSE., |
73 |
cAdjFreq=0., |
74 |
#- long run (> 1.yr): |
75 |
#pChkptFreq =311040000., |
76 |
#chkptFreq = 77760000., |
77 |
pChkptFreq =248832000., |
78 |
chkptFreq = 62208000., |
79 |
dumpFreq = 31104000., |
80 |
monitorFreq= 2592000., |
81 |
#- short run (~ months): |
82 |
#pChkptFreq = 2592000., |
83 |
#chkptFreq = 864000., |
84 |
#dumpFreq = 432000., |
85 |
#monitorFreq= 43200., |
86 |
#- short test: |
87 |
nTimeSteps=40, |
88 |
#dumpFreq =1., |
89 |
monitorFreq=1., |
90 |
& |
91 |
|
92 |
# Gridding parameters |
93 |
&PARM04 |
94 |
usingCurvilinearGrid=.TRUE., |
95 |
horizGridFile='dxC1_dXYa', |
96 |
radius_fromHorizGrid=6370.E3, |
97 |
#delR=25*40.E2, |
98 |
delR= 1500., 2122., 3000., 4242., 6000., 7800., 9100., 9400., 9280., 8400., |
99 |
6980., 5800., 4820., 4004., 3328., 2766., 2298., 1910., 1586., 1318., |
100 |
1096., 910., 756., 628., 522., 434., |
101 |
& |
102 |
|
103 |
# Input datasets |
104 |
&PARM05 |
105 |
hydrogThetaFile='ini_theta_26l.bin', |
106 |
hydrogSaltFile ='ini_specQ_26l.bin', |
107 |
& |