MLAdjust/input.0.smag/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
 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=.TRUE.,
# 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=43200.,
#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	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.0 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=.FALSE.,
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	jmc	1.7	&
72	baylor	1.1
73			# Elliptic solver parameters
74			&PARM02
75			cg2dMaxIters=1000,
76			cg2dTargetResidual=1.E-11,
77			cg3dMaxIters=40,
78			cg3dTargetResidual=1.E-9,
79	jmc	1.7	&
80	baylor	1.1
81			# Time stepping parameters
82			&PARM03
83	jmc	1.7	nIter0=0,
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	jmc	1.8	monitorSelect=2,
95	baylor	1.1	monitorFreq=1200.,
96			cadjfreq=0.,
97	jmc	1.7	&
98	baylor	1.1
99			# Gridding parameters
100			&PARM04
101			usingCartesianGrid=.TRUE.,
102	jmc	1.7	delR= 40*5.0,
103			dxSpacing=1.E3,
104			dySpacing=1.E3,
105			&
106	baylor	1.1
107			# Input datasets
108			&PARM05
109			bathyFile='topo_sl.bin',
110			hydrogThetaFile='thetaInitial.bin',
111			hydrogSaltFile='spiceInitial.bin',
112	jmc	1.7	&