4 |
# |
# |
5 |
# Continuous equation parameters |
# Continuous equation parameters |
6 |
&PARM01 |
&PARM01 |
7 |
viscAr=5.E-3, |
viscAr=5.E-5, |
8 |
# This is the minimum viscosity used total=viscah+viscleith |
# This is the background viscosity used total=viscah+viscleith |
9 |
# Often U*deltaX/2 for grid noise (Griffies and Hallberg,2000) |
# viscAh=1.E2, |
10 |
viscAh=3.E1, |
# 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 |
# This is the standard Leith Viscosity, prop to grad(curl v_h) |
# This is the standard Leith Viscosity, prop to grad(curl v_h) |
26 |
viscC2leith=5., |
# viscC2leith=0.2, |
27 |
|
viscC4leith=0.2, |
28 |
# This is an additional Leith Viscosity, prop to grad(div.v_h) |
# 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) |
# It traps instabilities that have no signal in grad(curl v_h) |
30 |
# viscC2leithD=5., |
# viscC2leithD=0.2, |
31 |
# These are the minimum viscosities used with Smagorinsky |
viscC4leithD=0.2, |
|
# total=viscastrain+viscsmag, etc. |
|
|
# viscAstrain=3.E1, |
|
|
# viscAtension=3.E1, |
|
32 |
# This is harmonic Smagorinsky Coefficient |
# This is harmonic Smagorinsky Coefficient |
33 |
# 0.15=value from 3-d turbulence theory; |
# 0.2-0.9=value from 3-d turbulence theory; Smagorinsky (1993) |
34 |
# 3-4=stable oceanic values (Griffies and Hallberg, 2000) |
# 2.2-4=stable oceanic values (Griffies and Hallberg, 2000) |
35 |
# viscC2Smag=4., |
# 0.7 c2smag is approximately the same as 0.12 c2leith+c2leithd |
36 |
|
# 1.0 c2smag is approximately the same as 0.25 c2leith+c2leithd |
37 |
|
# viscC2Smag=1.0, |
38 |
|
viscC4Smag=1.0, |
39 |
# This is the max Viscosity used, be it Smag or Leith, and it scales |
# This is the max Viscosity used, be it Smag or Leith, and it scales |
40 |
# with grid size and timestep |
# with grid size and timestep |
41 |
# Often 0.5 for stability (Griffies and Hallberg,2000) |
# viscAhGridMax=1.0, |
42 |
viscAhGridMax=0.5, |
viscA4GridMax=1.0, |
43 |
|
# This uses the full grad(vort) and grad(div) |
44 |
useFullLeith=.TRUE., |
useFullLeith=.TRUE., |
45 |
|
# This uses the Strain-Tension form for calculating viscous terms |
46 |
|
# useStrainTensionVisc=.TRUE., |
47 |
no_slip_sides=.FALSE., |
no_slip_sides=.FALSE., |
48 |
no_slip_bottom=.FALSE., |
no_slip_bottom=.FALSE., |
49 |
ivdc_kappa=10., |
ivdc_kappa=10., |
82 |
# Time stepping parameters |
# Time stepping parameters |
83 |
&PARM03 |
&PARM03 |
84 |
niter0=36, |
niter0=36, |
85 |
endTime=57600., |
#endTime=86400., |
86 |
|
#- for testing, only run for 12 iterations: |
87 |
|
nTimeSteps=12, |
88 |
deltaT=1200., |
deltaT=1200., |
89 |
abEps=0.01, |
abEps=0.01, |
90 |
pchkptFreq=7200., |
#- for testing, reduce amount of output by commenting out output frequencies: |
91 |
taveFreq=7200., |
# pchkptFreq=43200., |
92 |
dumpFreq=7200., |
#taveFreq=7200., |
93 |
diagFreq=7200., |
#dumpFreq=7200., |
94 |
|
#diagFreq=7200., |
95 |
monitorFreq=1200., |
monitorFreq=1200., |
96 |
cadjfreq=0., |
cadjfreq=0., |
97 |
/ |
/ |