1 |
#:::==================== |
2 |
# | Model parameters | |
3 |
# ==================== |
4 |
# |
5 |
# Continuous equation parameters |
6 |
&PARM01 |
7 |
# momstepping=.FALSE., |
8 |
# tempstepping=.FALSE., |
9 |
# saltstepping=.FALSE., |
10 |
Tref = 50*-1.9, |
11 |
Sref = 50*34.4, |
12 |
viscAr=1.E-3, |
13 |
#viscAh=600.0, |
14 |
viscAhGrid=0.2, |
15 |
viscA4Grid=0.02, |
16 |
no_slip_sides=.FALSE., |
17 |
no_slip_bottom=.FALSE., |
18 |
diffKhT= 100.0, |
19 |
#diffKhT= 600.0, |
20 |
diffKrT=5.E-5, |
21 |
diffKhS= 100.0, |
22 |
#diffKhS= 600.0, |
23 |
diffKrS=5.E-5, |
24 |
bottomDragQuadratic=2.5E-3, |
25 |
staggerTimestep = .true., |
26 |
tempAdvScheme=30, |
27 |
saltAdvScheme=30, |
28 |
vectorInvariantMomentum = .true., |
29 |
### momImplVertAdv=.true., |
30 |
tempImplVertAdv=.true., |
31 |
saltImplVertAdv=.true., |
32 |
eosType='JMD95Z', |
33 |
#tAlpha = 2.E-4, |
34 |
#sBeta = 7.4E-4, |
35 |
HeatCapacity_cp = 3974.0, |
36 |
rhoConst=1000., |
37 |
rhoNil=1000., |
38 |
gravity=9.81, |
39 |
gBaro=9.81, |
40 |
implicitDiffusion = .true., |
41 |
implicitViscosity = .true., |
42 |
rigidLid=.FALSE., |
43 |
implicitFreeSurface=.TRUE., |
44 |
useRealFreshWaterFlux = .true., |
45 |
exactConserv=.TRUE., |
46 |
# start nlfs with rstar |
47 |
#select_rStar=2, |
48 |
nonlinFreeSurf=4, |
49 |
hFacInf=0.2, |
50 |
hFacSup=2.0, |
51 |
# end nlfs with rstar |
52 |
hFacMin=0.2, |
53 |
#integr_GeoPot = 1, |
54 |
#ivdc_kappa = 0.005., |
55 |
#implicitDiffusion=.true., |
56 |
#implicitViscosity=.true., |
57 |
#useCDscheme=.true., |
58 |
#nonHydrostatic=.true., |
59 |
#allowFreezing = .true., |
60 |
readBinaryPrec=64, |
61 |
#writeBinaryPrec=64, |
62 |
#debuglevel = 1, |
63 |
selectCoriMap = 0, |
64 |
f0 = 0.0, |
65 |
useSingleCpuIO=.TRUE., |
66 |
globalFiles=.TRUE., |
67 |
useCoriolis=.TRUE., |
68 |
& |
69 |
|
70 |
# Elliptic solver parameters |
71 |
&PARM02 |
72 |
cg2dMaxIters=300, |
73 |
cg2dTargetResidual=1.E-11, |
74 |
#cg3dMaxIters=40, |
75 |
#cg3dTargetResidual=1.E-13, |
76 |
& |
77 |
|
78 |
#Time stepping parameters |
79 |
&PARM03 |
80 |
# niter0=0, |
81 |
# nTimeSteps=25920000., |
82 |
nTimeSteps=12, |
83 |
startTime=2592000., |
84 |
# startTime = 0., |
85 |
#startTime=62208000., |
86 |
#startTime=15552000., |
87 |
#endTime=31104000., |
88 |
deltaT=300.0, |
89 |
### deltaT=1200.0, |
90 |
forcing_In_AB = .false., |
91 |
#rCD = 0.9896, |
92 |
abEps=0.1, |
93 |
cAdjFreq = -1., |
94 |
pChkptFreq=2592000, |
95 |
chkptFreq=2592000.0, |
96 |
# pChkptFreq=600.0, |
97 |
# chkptFreq=600.0, |
98 |
dumpFreq=15552000, |
99 |
monitorFreq=15552000., |
100 |
monitorSelect=1, |
101 |
monitorFreq=1., |
102 |
& |
103 |
|
104 |
# Gridding parameters |
105 |
&PARM04 |
106 |
usingCartesianGrid=.FALSE., |
107 |
usingSphericalPolarGrid=.TRUE., |
108 |
#delR = 10., 10., 10., 10., 10., 10., 10., 10.01, 10.03, 10.11, |
109 |
#10.32, 10.8, 11.76, 13.42, 16.04, 19.82, 24.85, 31.1, 38.42, 46.5, |
110 |
#55., 63.5, 71.58, 78.9, 85.15, 90.18, 93.96, 96.58, 98.25, 99.25, |
111 |
#100.01, 101.33, 104.56, 111.33, 122.83, 139.09, 158.94, 180.83, 203.55, 226.5, |
112 |
#249.5, 272.5, 295.5, 318.5, 341.5, 364.5, 387.5, 410.5, 433.5, 456.5, |
113 |
delR=130*10., |
114 |
delX=3*.125, |
115 |
delY=400*.0078125, |
116 |
# delYFile='dlat.bin', |
117 |
# delY=30*10e3, |
118 |
xgOrigin = -105.5, |
119 |
ygOrigin = -75.4457, |
120 |
& |
121 |
|
122 |
# Input datasets |
123 |
&PARM05 |
124 |
bathyFile='bathy101.box', |
125 |
hydrogThetaFile='theta.init', |
126 |
hydrogSaltFile='salt.init', |
127 |
pSurfInitFile='etainit.round.bin', |
128 |
& |