1 |
# ==================== |
2 |
# | Model parameters | |
3 |
# ==================== |
4 |
# |
5 |
# Continuous equation parameters |
6 |
&PARM01 |
7 |
tRef = 3*23.,3*22.,21.,2*20.,19.,2*18.,17.,2*16.,15.,14.,13., |
8 |
12.,11.,2*9.,8.,7.,2*6.,2*5.,3*4.,3*3.,4*2.,12*1., |
9 |
sRef = 50*34.5, |
10 |
no_slip_sides = .FALSE., |
11 |
no_slip_bottom = .TRUE., |
12 |
# |
13 |
## viscAh = 1.E3, |
14 |
## viscAr = 1.E-4, |
15 |
## diffKhT=1.E3, |
16 |
## diffKhS=1.E3, |
17 |
## diffK4T=0., |
18 |
## diffK4S=0., |
19 |
## diffKrT=2.E-5, |
20 |
## diffKrS=2.E-5, |
21 |
# |
22 |
viscAr=1.e-4, |
23 |
### viscAr=1.e-4, |
24 |
viscAh=0.5e4, |
25 |
### viscAh=1.0e4, |
26 |
### viscA4=1.E10 |
27 |
diffKhT=1.E2, |
28 |
diffKhS=1.E2, |
29 |
# diffK4T=1.E10, |
30 |
# diffK4S=1.E10, |
31 |
diffKrT=1.e-5, |
32 |
diffKrS=1.e-5, |
33 |
# |
34 |
### diffKrBL79surf=0.1E-4, |
35 |
### diffKrBL79deep=1.0E-4, |
36 |
bottomDragQuadratic = 0.0001, |
37 |
### bottomDragQuadratic = 0.001, |
38 |
### ivdc_kappa=100., |
39 |
implicitDiffusion=.TRUE., |
40 |
implicitViscosity=.TRUE., |
41 |
gravity=9.81, |
42 |
rhonil=1027., |
43 |
rhoConstFresh=1000., |
44 |
eosType='JMD95Z', |
45 |
implicitFreeSurface=.TRUE., |
46 |
exactConserv=.TRUE., |
47 |
useSingleCpuIO=.TRUE., |
48 |
tempAdvScheme=30, |
49 |
saltAdvScheme=30, |
50 |
#ph( |
51 |
staggerTimeStep=.TRUE., |
52 |
vectorInvariantMomentum=.FALSE., |
53 |
multiDimAdvection=.FALSE. |
54 |
useCDscheme=.TRUE., |
55 |
useJamartWetPoints=.TRUE., |
56 |
#ph) |
57 |
### select_rStar=2, |
58 |
### hFacInf=0.2, |
59 |
### hFacSup=2.0, |
60 |
### nonlinFreeSurf=1, |
61 |
useRealFreshWaterFlux=.FALSE., |
62 |
allowFreezing=.TRUE., |
63 |
#ph( |
64 |
# hFacMin=.3, |
65 |
# hFacMinDr=50., |
66 |
#ph) |
67 |
readBinaryPrec=64, |
68 |
writeBinaryPrec=64, |
69 |
debugLevel=1, |
70 |
/ |
71 |
|
72 |
# Elliptic solver parameters |
73 |
&PARM02 |
74 |
cg2dMaxIters=500, |
75 |
cg2dTargetResWunit=1.E-15, |
76 |
/ |
77 |
|
78 |
# Time stepping parameters |
79 |
&PARM03 |
80 |
nIter0=1, |
81 |
# 1 year at 600sec |
82 |
### nTimeSteps=52560, |
83 |
# 1 year at 300sec |
84 |
### nTimeSteps=105120, |
85 |
### nTimeSteps=8640, |
86 |
# 30.5 days at 60sec |
87 |
endTime=2635200., |
88 |
### endTime=172800., |
89 |
### endTime=31536000., |
90 |
# |
91 |
deltaT = 3600., |
92 |
#deltaTmom =1200., |
93 |
#deltaTtracer=86400., |
94 |
#deltaTfreesurf=86400., |
95 |
#deltaTClock =86400., |
96 |
abEps = 0.1, |
97 |
tauCD=172800.0, |
98 |
cAdjFreq=1., |
99 |
forcing_In_AB=.FALSE., |
100 |
momDissip_In_AB=.FALSE., |
101 |
pChkptFreq =2635200.0, |
102 |
chkptFreq =2635200.0, |
103 |
# taveFreq =2635200.0, |
104 |
dumpFreq = 864000.0, |
105 |
monitorFreq = 43200.0, |
106 |
#latBandClimRelax=50., |
107 |
### tauThetaClimRelax= 2592000., |
108 |
### tauSaltClimRelax = 2592000., |
109 |
pickupStrictlyMatch=.FALSE., |
110 |
/ |
111 |
|
112 |
# Gridding parameters |
113 |
&PARM04 |
114 |
usingCurvilinearGrid=.TRUE., |
115 |
###horizGridFile='grid_cs32', |
116 |
delR = |
117 |
10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01, |
118 |
10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04, 19.82, 24.85, |
119 |
31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18, |
120 |
93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83, |
121 |
139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50, |
122 |
341.50,364.50,387.50,410.50,433.50,456.50, |
123 |
/ |
124 |
|
125 |
# Input datasets |
126 |
&PARM05 |
127 |
### bathyFile ='llc96x96x288_bathy_mindepth_noisolated.bin', |
128 |
bathyFile ='llc96x96x288_bathy_shap_filt_44.bin', |
129 |
# hydrogThetaFile='DiagTheta_diffus_shap_44.0000000144.bin', |
130 |
# hydrogSaltFile ='DiagSalt_diffus_shap_44.0000000144.bin', |
131 |
hydrogThetaFile='DiagTheta_diffus_nopart.0000000072.bin', |
132 |
hydrogSaltFile ='DiagSalt_diffus_nopart.0000000072.bin', |
133 |
# |
134 |
/ |