1 |
# ==================== |
2 |
# | Model parameters | |
3 |
# ==================== |
4 |
# |
5 |
# Continuous equation parameters |
6 |
&PARM01 |
7 |
tRef=295.1, 295.3, 295.4, 295.6, 295.8, 296.0, 296.2, 296.4, 296.6, 296.8, |
8 |
297.0, 297.3, 297.5, 297.7, 298.0, 298.3, 298.5, 298.8, 299.1, 299.5, |
9 |
299.8, 300.1, 300.5, 300.9, 301.5, 302.3, 303.4, 304.9, 306.8, 309.4, |
10 |
312.9, 317.8, 324.7, 336.1, 352.6, 373.4, 401.2, 441.6, 511.0, 699.5, |
11 |
#tRef=295.2, 295.5, 295.9, 296.3, 296.7, 297.1, 297.6, 298.1, 298.7, 299.3, |
12 |
# 300.0, 300.7, 301.9, 304.1, 308.0, 315.1, 329.5, 362.3, 419.2, 573.8, |
13 |
#tRef=295.3, 296.1, 296.9, 297.9, 299.0, 300.3, 302.8, 311.0, 343.9, 470.7, |
14 |
sRef=40*0.0, |
15 |
#sRef=1.0E-3,5.0E-4,1.E-4,5*1.E-6,1.E-8,1.E-9, |
16 |
viscAr=0.E1, |
17 |
viscAh=0.E6, |
18 |
viscA4=0.E17, |
19 |
no_slip_sides=.FALSE., |
20 |
no_slip_bottom=.FALSE., |
21 |
diffKhT=0.E3, |
22 |
diffKrT=0., |
23 |
diffK4T=0.E17, |
24 |
diffKrS=0.E2, |
25 |
diffKhS=0.E3, |
26 |
diffK4S=0.E17, |
27 |
buoyancyRelation='ATMOSPHERIC', |
28 |
eosType='IDEALGAS', |
29 |
atm_Rq=0.6078, |
30 |
#Integr_GeoPot=2, |
31 |
rotationPeriod=86400., |
32 |
gravity=9.81, |
33 |
rhonil=1.0, |
34 |
rhoConst=1.0, |
35 |
rigidLid=.FALSE., |
36 |
implicitFreeSurface=.TRUE., |
37 |
exactConserv=.TRUE., |
38 |
uniformLin_PhiSurf=.FALSE., |
39 |
staggerTimeStep=.TRUE., |
40 |
nonHydrostatic=.FALSE., |
41 |
readBinaryPrec=64, |
42 |
writeBinaryPrec=64, |
43 |
hFacMin=0.5, |
44 |
nonlinFreeSurf=4, |
45 |
select_rStar=2, |
46 |
saltAdvScheme=2, |
47 |
vectorInvariantMomentum=.TRUE., |
48 |
#- here are the new things: |
49 |
useAbsVorticity=.TRUE., |
50 |
SadournyCoriolis=.TRUE., |
51 |
selectKEscheme=3, |
52 |
& |
53 |
|
54 |
# Elliptic solver parameters |
55 |
&PARM02 |
56 |
cg2dMaxIters=200, |
57 |
#cg2dTargetResidual=1.E-12, |
58 |
cg2dTargetResWunit=1.E-17, |
59 |
& |
60 |
|
61 |
# Time stepping parameters |
62 |
&PARM03 |
63 |
nIter0=0, |
64 |
nTimeSteps=6, |
65 |
#endTime=864000.0, |
66 |
deltaT=120.0, |
67 |
forcing_In_AB=.FALSE., |
68 |
abEps=0.1, |
69 |
cAdjFreq=0., |
70 |
pChkptFreq=0.0, |
71 |
ChkptFreq=0.0, |
72 |
dumpFreq=0.0, |
73 |
dumpInitAndLast=.FALSE., |
74 |
monitorFreq=120., |
75 |
taveFreq=0.0, |
76 |
tave_lastiter=1., |
77 |
# for daily testing: to save disk space |
78 |
writePickupAtEnd=.FALSE., |
79 |
& |
80 |
|
81 |
# Gridding parameters |
82 |
&PARM04 |
83 |
usingCartesianGrid=.FALSE., |
84 |
usingSphericalPolarGrid=.FALSE., |
85 |
usingCurvilinearGrid=.TRUE., |
86 |
horizGridFile='dxC1_dXYa', |
87 |
#delR=40*0.25325E4, |
88 |
#Ro_SeaLevel=1.013E5, |
89 |
delR=40*0.25E4, |
90 |
Ro_SeaLevel=1.E5, |
91 |
rSphere=6371.E3, |
92 |
& |
93 |
|
94 |
# Input datasets |
95 |
&PARM05 |
96 |
topoFile='top19232.data', |
97 |
uVelInitFile='uVel_cs2-32x40.d19930110z06.bin', |
98 |
vVelInitFile='vVel_cs2-32x40.d19930110z06.bin', |
99 |
hydrogThetaFile='TH_cs2-32x40.d19930110z06.bin', |
100 |
hydrogSaltFile='SHwet_cs32x40.d19930110z06.bin', |
101 |
& |