1 |
# ==================== |
2 |
# | Model parameters | |
3 |
# ==================== |
4 |
# |
5 |
# Continuous equation parameters |
6 |
&PARM01 |
7 |
Tref = -0.194175, -0.582524, -0.970874, -1.35922, -1.74757, |
8 |
-2.13592, -2.52427, -2.91262, -3.30097, -3.68932, |
9 |
-4.07767, -4.46602, -4.85437, -5.24272, -5.63107, |
10 |
-6.01942, -6.40777, -6.79612, -7.18447, -7.57282, |
11 |
-7.96117, -8.34951, -8.73786, -9.12621, -9.51456, |
12 |
sRef= 25*1., |
13 |
viscAz=5.E-1, |
14 |
viscAh=0.0E0, |
15 |
no_slip_sides=.TRUE., |
16 |
no_slip_bottom=.TRUE., |
17 |
bottomDragQuadratic=2.E-3, |
18 |
viscA4=0.E12, |
19 |
viscC4Leith=0.1, |
20 |
diffKhT=0.0E0 |
21 |
diffKzT=0.E-2, |
22 |
diffKhS=0.E0, |
23 |
diffKzS=0.E-2, |
24 |
f0=1.0e-4, |
25 |
beta=0.E-11, |
26 |
tAlpha=2.E-4, |
27 |
sBeta =0.E-4, |
28 |
gravity=9.81, |
29 |
rhonil=1030.0, |
30 |
rigidLid=.FALSE., |
31 |
implicitFreeSurface=.TRUE., |
32 |
eosType='LINEAR', |
33 |
hFacMin=0.2, |
34 |
nonHydrostatic=.FALSE., |
35 |
multiDimAdvection=.TRUE., |
36 |
tempAdvScheme=77, |
37 |
saltAdvScheme=77, |
38 |
tracerAdvScheme=77, |
39 |
readBinaryPrec=64, |
40 |
writeBinaryPrec=64, |
41 |
writeStatePrec=64, |
42 |
globalFiles=.TRUE., |
43 |
staggerTimeStep=.TRUE., |
44 |
& |
45 |
|
46 |
# Elliptic solver parameters |
47 |
&PARM02 |
48 |
cg2dMaxIters=1000, |
49 |
cg2dTargetResidual=1.E-11, |
50 |
cg3dMaxIters=10, |
51 |
cg3dTargetResidual=1.E-13, |
52 |
& |
53 |
|
54 |
# Time stepping parameters |
55 |
&PARM03 |
56 |
niter0=0, |
57 |
#nTimeSteps=12000, |
58 |
nTimeSteps=20, |
59 |
deltaT=300., |
60 |
abEps=0.1, |
61 |
pChkptFreq=0.0, |
62 |
chkptFreq=100000.0, |
63 |
dumpFreq=30000., |
64 |
#cAdjFreq=30., |
65 |
#monitorFreq=6000., |
66 |
monitorFreq=300., |
67 |
& |
68 |
|
69 |
# Gridding parameters |
70 |
&PARM04 |
71 |
usingCartesianGrid=.TRUE., |
72 |
usingSphericalPolarGrid=.FALSE., |
73 |
# delXfile='delXvar' |
74 |
delX=200*10000.0, |
75 |
delYfile='delYvar', |
76 |
delZ = 25*144.0, |
77 |
& |
78 |
|
79 |
# Input datasets |
80 |
&PARM05 |
81 |
hydrogThetaFile='T.init', |
82 |
bathyFile='topog.slope', |
83 |
& |