MLAdjust/input.0.leithD/data

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 viscAr=5.E-5,
# This is the background viscosity used total=viscah+viscleith
# viscAh=1.E2,
# viscA4=1.E2,
# This is the background viscosity used on vorticity points
# viscAhZ=1.E2,
# viscA4Z=1.E2,
# This is the background viscosity used on divergence points
# viscAhD=1.E2,
# viscA4D=1.E2,
# viscAhd and viscahz are set to viscah if they are unset initially
#
# This is a maximum gridscale Reynolds number limiter
# ViscAhRemax=10,
# ViscA4Remax=10,
# This is a minimum autoscaling viscosity
# ViscAhGridMin=.01,
# ViscA4GridMin=.01,
# This is the standard Leith Viscosity, prop to grad(curl v_h)
 viscC2leith=1.85,
# viscC4leith=1.85,
# This is an additional Leith Viscosity, prop to grad(div.v_h)
# It traps instabilities that have no signal in grad(curl v_h)
 viscC2leithD=1.85,
# viscC4leithD=1.85,
# This is harmonic Smagorinsky Coefficient
# 0.2-0.9=value from 3-d turbulence theory; Smagorinsky (1993)
# 2.2-4=stable oceanic values (Griffies and Hallberg, 2000)
# 1.0 smag is approximately the same as 1.85 leith+leithd
# for this grid
 viscC2Smag=1.0, 
# viscC4Smag=1.0,
# This is the max Viscosity used, be it Smag or Leith, and it scales 
# with grid size and timestep
 viscAhGridMax=1.0,
# viscA4GridMax=1.0,
# This uses the full grad(vort) and grad(div)
 useFullLeith=.FALSE.,
# This uses the Strain-Tension form for calculating viscous terms
 useStrainTensionVisc=.FALSE.,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.FALSE.,
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 f0=7.29e-5,
 beta=0.E-11,
 eosType='LINEAR',
# Using Temperature as only Active Tracer
 tAlpha=2.E-4,
# Using Salinity as Spice
 sBeta =0.E-4,
 saltStepping=.TRUE.,
 gravity=9.81,
 rhonil=1035.,
 implicitFreeSurface=.TRUE.,
 usejamartwetpoints=.TRUE.,
 rigidLid=.FALSE.,
 staggerTimeStep=.TRUE.,
 readBinaryPrec=64,
 writeBinaryPrec=64,
 tempAdvScheme=33,
 saltAdvScheme=33,
 nonHydrostatic=.FALSE.,
 vectorInvariantMomentum=.TRUE.,
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=1000,
 cg2dTargetResidual=1.E-11,
 cg3dMaxIters=40,
 cg3dTargetResidual=1.E-9,
 &

# Time stepping parameters
 &PARM03
 nIter0=0,
#endTime=86400.,
#- for testing, only run for 12 iterations:
 nTimeSteps=12,
 deltaT=1200.,
 abEps=0.01,
#- for testing, reduce amount of output by commenting out output frequencies:
# pchkptFreq=7200.,
# taveFreq=7200.,
# dumpFreq=7200.,
# diagFreq=7200.,
 monitorSelect=2,
 monitorFreq=1200.,
 cadjfreq=0.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 delR= 40*5.0,
 dxSpacing=1.E3,
 dySpacing=1.E3,
 &

# Input datasets
 &PARM05
 bathyFile='topo_sl.bin',
 hydrogThetaFile='thetaInitial.bin',
 hydrogSaltFile='spiceInitial.bin',
 &
1	# ====================
2	# \| Model parameters \|
3	# ====================
4	#
5	# Continuous equation parameters
6	&PARM01
7	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	# This is the standard Leith Viscosity, prop to grad(curl v_h)
26	viscC2leith=1.85,
27	# viscC4leith=1.85,
28	# 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	viscC2leithD=1.85,
31	# viscC4leithD=1.85,
32	# This is harmonic Smagorinsky Coefficient
33	# 0.2-0.9=value from 3-d turbulence theory; Smagorinsky (1993)
34	# 2.2-4=stable oceanic values (Griffies and Hallberg, 2000)
35	# 1.0 smag is approximately the same as 1.85 leith+leithd
36	# for this grid
37	viscC2Smag=1.0,
38	# viscC4Smag=1.0,
39	# This is the max Viscosity used, be it Smag or Leith, and it scales
40	# with grid size and timestep
41	viscAhGridMax=1.0,
42	# viscA4GridMax=1.0,
43	# This uses the full grad(vort) and grad(div)
44	useFullLeith=.FALSE.,
45	# This uses the Strain-Tension form for calculating viscous terms
46	useStrainTensionVisc=.FALSE.,
47	no_slip_sides=.FALSE.,
48	no_slip_bottom=.FALSE.,
49	ivdc_kappa=10.,
50	implicitDiffusion=.TRUE.,
51	implicitViscosity=.TRUE.,
52	f0=7.29e-5,
53	beta=0.E-11,
54	eosType='LINEAR',
55	# Using Temperature as only Active Tracer
56	tAlpha=2.E-4,
57	# Using Salinity as Spice
58	sBeta =0.E-4,
59	saltStepping=.TRUE.,
60	gravity=9.81,
61	rhonil=1035.,
62	implicitFreeSurface=.TRUE.,
63	usejamartwetpoints=.TRUE.,
64	rigidLid=.FALSE.,
65	staggerTimeStep=.TRUE.,
66	readBinaryPrec=64,
67	writeBinaryPrec=64,
68	tempAdvScheme=33,
69	saltAdvScheme=33,
70	nonHydrostatic=.FALSE.,
71	vectorInvariantMomentum=.TRUE.,
72	&
73
74	# Elliptic solver parameters
75	&PARM02
76	cg2dMaxIters=1000,
77	cg2dTargetResidual=1.E-11,
78	cg3dMaxIters=40,
79	cg3dTargetResidual=1.E-9,
80	&
81
82	# Time stepping parameters
83	&PARM03
84	nIter0=0,
85	#endTime=86400.,
86	#- for testing, only run for 12 iterations:
87	nTimeSteps=12,
88	deltaT=1200.,
89	abEps=0.01,
90	#- for testing, reduce amount of output by commenting out output frequencies:
91	# pchkptFreq=7200.,
92	# taveFreq=7200.,
93	# dumpFreq=7200.,
94	# diagFreq=7200.,
95	monitorSelect=2,
96	monitorFreq=1200.,
97	cadjfreq=0.,
98	&
99
100	# Gridding parameters
101	&PARM04
102	usingCartesianGrid=.TRUE.,
103	delR= 40*5.0,
104	dxSpacing=1.E3,
105	dySpacing=1.E3,
106	&
107
108	# Input datasets
109	&PARM05
110	bathyFile='topo_sl.bin',
111	hydrogThetaFile='thetaInitial.bin',
112	hydrogSaltFile='spiceInitial.bin',
113	&