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=0.2,
# viscC4leith=0.2,
# 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=0.2,
# viscC4leithD=0.2,
# 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)
# 0.7 c2smag is approximately the same as 0.12 c2leith+c2leithd
# 1.0 c2smag is approximately the same as 0.25 c2leith+c2leithd
 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.,
 monitorFreq=1200.,
 cadjfreq=0.,
 /

# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 usingSphericalPolarGrid=.FALSE.,
  delR= 5d0, 5d0, 5d0, 5d0, 5d0,
    5d0, 5d0, 5d0, 5d0, 5d0,
    5d0, 5d0, 5d0, 5d0, 5d0,
    5d0, 5d0, 5d0, 5d0, 5d0,
    5d0, 5d0, 5d0, 5d0, 5d0,
    5d0, 5d0, 5d0, 5d0, 5d0,
    5d0, 5d0, 5d0, 5d0, 5d0,
    5d0, 5d0, 5d0, 5d0, 5d0,
 dxspacing=1d3,
 dyspacing=1d3,
 /

# 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.3	viscC2leith=0.2,
27			# viscC4leith=0.2,
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.3	viscC2leithD=0.2,
31			# viscC4leithD=0.2,
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			# 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	baylor	1.1	# This is the max Viscosity used, be it Smag or Leith, and it scales
40			# with grid size and timestep
41	baylor	1.3	viscAhGridMax=1.0,
42			# viscA4GridMax=1.0,
43			# This uses the full grad(vort) and grad(div)
44	baylor	1.1	useFullLeith=.TRUE.,
45	baylor	1.3	# This uses the Strain-Tension form for calculating viscous terms
46			useStrainTensionVisc=.FALSE.,
47	baylor	1.1	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	baylor	1.3	#endTime=86400.,
86			#- for testing, only run for 12 iterations:
87			nTimeSteps=12,
88	baylor	1.1	deltaT=1200.,
89			abEps=0.01,
90	baylor	1.3	#- for testing, reduce amount of output by commenting out output frequencies:
91			#pchkptFreq=43200.,
92			#taveFreq=7200.,
93			#dumpFreq=7200.,
94			#diagFreq=7200.,
95	baylor	1.1	monitorFreq=1200.,
96			cadjfreq=0.,
97			/
98
99			# Gridding parameters
100			&PARM04
101			usingCartesianGrid=.TRUE.,
102			usingSphericalPolarGrid=.FALSE.,
103			delR= 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			5d0, 5d0, 5d0, 5d0, 5d0,
111			dxspacing=1d3,
112			dyspacing=1d3,
113			/
114
115			# Input datasets
116			&PARM05
117			bathyFile='topo_sl.bin',
118			hydrogThetaFile='thetaInitial.bin',
119			hydrogSaltFile='spiceInitial.bin',
120			/