1 |
baylor |
1.1 |
# ==================== |
2 |
|
|
# | Model parameters | |
3 |
|
|
# ==================== |
4 |
|
|
# |
5 |
|
|
# Continuous equation parameters |
6 |
|
|
&PARM01 |
7 |
baylor |
1.3 |
viscAr=5.E-5, |
8 |
|
|
# This is the background viscosity used total=viscah+viscleith |
9 |
|
|
# viscAh=1.E2, |
10 |
|
|
# viscA4=1.E2, |
11 |
|
|
# This is the background viscosity used on vorticity points |
12 |
|
|
# viscAhZ=1.E2, |
13 |
|
|
# viscA4Z=1.E2, |
14 |
|
|
# This is the background viscosity used on divergence points |
15 |
|
|
# viscAhD=1.E2, |
16 |
|
|
# viscA4D=1.E2, |
17 |
|
|
# viscAhd and viscahz are set to viscah if they are unset initially |
18 |
|
|
# |
19 |
|
|
# This is a maximum gridscale Reynolds number limiter |
20 |
|
|
# ViscAhRemax=10, |
21 |
|
|
# ViscA4Remax=10, |
22 |
|
|
# This is a minimum autoscaling viscosity |
23 |
|
|
# ViscAhGridMin=.01, |
24 |
|
|
# ViscA4GridMin=.01, |
25 |
baylor |
1.1 |
# This is the standard Leith Viscosity, prop to grad(curl v_h) |
26 |
baylor |
1.4 |
# viscC2leith=1.85, |
27 |
|
|
viscC4leith=1.85, |
28 |
baylor |
1.1 |
# This is an additional Leith Viscosity, prop to grad(div.v_h) |
29 |
|
|
# It traps instabilities that have no signal in grad(curl v_h) |
30 |
baylor |
1.4 |
# viscC2leithD=1.85, |
31 |
|
|
viscC4leithD=1.85, |
32 |
baylor |
1.1 |
# This is harmonic Smagorinsky Coefficient |
33 |
baylor |
1.3 |
# 0.2-0.9=value from 3-d turbulence theory; Smagorinsky (1993) |
34 |
|
|
# 2.2-4=stable oceanic values (Griffies and Hallberg, 2000) |
35 |
baylor |
1.4 |
# 1.85 smag is approximately the same as 1.85 leith+leithd |
36 |
baylor |
1.3 |
# viscC2Smag=1.0, |
37 |
|
|
viscC4Smag=1.0, |
38 |
baylor |
1.1 |
# This is the max Viscosity used, be it Smag or Leith, and it scales |
39 |
|
|
# with grid size and timestep |
40 |
baylor |
1.3 |
# viscAhGridMax=1.0, |
41 |
|
|
viscA4GridMax=1.0, |
42 |
|
|
# This uses the full grad(vort) and grad(div) |
43 |
baylor |
1.1 |
useFullLeith=.TRUE., |
44 |
baylor |
1.3 |
# This uses the Strain-Tension form for calculating viscous terms |
45 |
|
|
# useStrainTensionVisc=.TRUE., |
46 |
baylor |
1.1 |
no_slip_sides=.FALSE., |
47 |
|
|
no_slip_bottom=.FALSE., |
48 |
|
|
ivdc_kappa=10., |
49 |
|
|
implicitDiffusion=.TRUE., |
50 |
|
|
implicitViscosity=.TRUE., |
51 |
|
|
f0=7.29e-5, |
52 |
|
|
beta=0.E-11, |
53 |
|
|
eosType='LINEAR', |
54 |
|
|
# Using Temperature as only Active Tracer |
55 |
|
|
tAlpha=2.E-4, |
56 |
|
|
# Using Salinity as Spice |
57 |
|
|
sBeta =0.E-4, |
58 |
|
|
saltStepping=.TRUE., |
59 |
|
|
gravity=9.81, |
60 |
|
|
rhonil=1035., |
61 |
|
|
implicitFreeSurface=.TRUE., |
62 |
|
|
usejamartwetpoints=.TRUE., |
63 |
|
|
rigidLid=.FALSE., |
64 |
|
|
staggerTimeStep=.TRUE., |
65 |
|
|
readBinaryPrec=64, |
66 |
|
|
writeBinaryPrec=64, |
67 |
|
|
tempAdvScheme=33, |
68 |
|
|
saltAdvScheme=33, |
69 |
|
|
nonHydrostatic=.FALSE., |
70 |
|
|
vectorInvariantMomentum=.TRUE., |
71 |
|
|
/ |
72 |
|
|
|
73 |
|
|
# Elliptic solver parameters |
74 |
|
|
&PARM02 |
75 |
|
|
cg2dMaxIters=1000, |
76 |
|
|
cg2dTargetResidual=1.E-11, |
77 |
|
|
cg3dMaxIters=40, |
78 |
|
|
cg3dTargetResidual=1.E-9, |
79 |
|
|
/ |
80 |
|
|
|
81 |
|
|
# Time stepping parameters |
82 |
|
|
&PARM03 |
83 |
|
|
niter0=36, |
84 |
baylor |
1.3 |
#endTime=86400., |
85 |
|
|
#- for testing, only run for 12 iterations: |
86 |
|
|
nTimeSteps=12, |
87 |
baylor |
1.1 |
deltaT=1200., |
88 |
|
|
abEps=0.01, |
89 |
baylor |
1.3 |
#- for testing, reduce amount of output by commenting out output frequencies: |
90 |
|
|
# pchkptFreq=43200., |
91 |
|
|
#taveFreq=7200., |
92 |
|
|
#dumpFreq=7200., |
93 |
|
|
#diagFreq=7200., |
94 |
baylor |
1.1 |
monitorFreq=1200., |
95 |
|
|
cadjfreq=0., |
96 |
|
|
/ |
97 |
|
|
|
98 |
|
|
# Gridding parameters |
99 |
|
|
&PARM04 |
100 |
|
|
usingCartesianGrid=.TRUE., |
101 |
|
|
usingSphericalPolarGrid=.FALSE., |
102 |
|
|
delR= 5d0, 5d0, 5d0, 5d0, 5d0, |
103 |
|
|
5d0, 5d0, 5d0, 5d0, 5d0, |
104 |
|
|
5d0, 5d0, 5d0, 5d0, 5d0, |
105 |
|
|
5d0, 5d0, 5d0, 5d0, 5d0, |
106 |
|
|
5d0, 5d0, 5d0, 5d0, 5d0, |
107 |
|
|
5d0, 5d0, 5d0, 5d0, 5d0, |
108 |
|
|
5d0, 5d0, 5d0, 5d0, 5d0, |
109 |
|
|
5d0, 5d0, 5d0, 5d0, 5d0, |
110 |
|
|
dxspacing=1d3, |
111 |
|
|
dyspacing=1d3, |
112 |
|
|
/ |
113 |
|
|
|
114 |
|
|
# Input datasets |
115 |
|
|
&PARM05 |
116 |
|
|
bathyFile='topo_sl.bin', |
117 |
|
|
hydrogThetaFile='thetaInitial.bin', |
118 |
|
|
hydrogSaltFile='spiceInitial.bin', |
119 |
|
|
/ |