MLAdjust/input.AhStTn/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. smag is approximately the same as 1.85 leith+leithd for this grid spacing
 viscC2smag = 1.,
#viscC4smag = 1.,
#- This is the max Viscosity used, be it Smag or Leith, and it scales
#  with grid size and timestep
 viscAhGridMax = 1.,
#viscA4GridMax = 1.,
#- This uses the full grad(vort) and grad(div)
 useFullLeith=.TRUE.,
#- This uses the Strain-Tension form for calculating viscous terms
 useStrainTensionVisc=.TRUE.,
#no_slip_sides=.FALSE.,
 no_slip_bottom=.FALSE.,
 vectorInvariantMomentum=.TRUE.,
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 f0=7.29E-5,
 beta=0.E-11,
 useJamartWetPoints=.TRUE.,
 eosType='LINEAR',
# Using Temperature as only Active Tracer
 tAlpha=2.E-4,
# Using Salinity as Spice
 sBeta =0.E-4,
 gravity=9.81,
 rhoNil=1035.,
 implicitFreeSurface=.TRUE.,
 staggerTimeStep=.TRUE.,
 tempAdvScheme=33,
 saltAdvScheme=33,
 readBinaryPrec=64,
 writeBinaryPrec=64,
 &

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

# Time stepping parameters
 &PARM03
 nIter0=0,
#endTime=864000.,
#- 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.,
 &

# 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	jmc	1.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. smag is approximately the same as 1.85 leith+leithd for this grid spacing
36			viscC2smag = 1.,
37			#viscC4smag = 1.,
38			#- This is the max Viscosity used, be it Smag or Leith, and it scales
39			# with grid size and timestep
40			viscAhGridMax = 1.,
41			#viscA4GridMax = 1.,
42			#- This uses the full grad(vort) and grad(div)
43			useFullLeith=.TRUE.,
44			#- This uses the Strain-Tension form for calculating viscous terms
45			useStrainTensionVisc=.TRUE.,
46			#no_slip_sides=.FALSE.,
47			no_slip_bottom=.FALSE.,
48			vectorInvariantMomentum=.TRUE.,
49			ivdc_kappa=10.,
50			implicitDiffusion=.TRUE.,
51			implicitViscosity=.TRUE.,
52			f0=7.29E-5,
53			beta=0.E-11,
54			useJamartWetPoints=.TRUE.,
55			eosType='LINEAR',
56			# Using Temperature as only Active Tracer
57			tAlpha=2.E-4,
58			# Using Salinity as Spice
59			sBeta =0.E-4,
60			gravity=9.81,
61			rhoNil=1035.,
62			implicitFreeSurface=.TRUE.,
63			staggerTimeStep=.TRUE.,
64			tempAdvScheme=33,
65			saltAdvScheme=33,
66			readBinaryPrec=64,
67			writeBinaryPrec=64,
68			&
69
70			# Elliptic solver parameters
71			&PARM02
72			cg2dMaxIters=1000,
73			cg2dTargetResidual=1.E-11,
74			cg3dMaxIters=40,
75			cg3dTargetResidual=1.E-9,
76			&
77
78			# Time stepping parameters
79			&PARM03
80			nIter0=0,
81			#endTime=864000.,
82			#- for testing, only run for 12 iterations:
83			nTimeSteps=12,
84			deltaT=1200.,
85			abEps=0.01,
86			#- for testing, reduce amount of output by commenting out output frequencies:
87			# pchkptFreq=43200.,
88			#taveFreq=7200.,
89			#dumpFreq=7200.,
90			#diagFreq=7200.,
91			monitorSelect=2,
92			monitorFreq=1200.,
93			&
94
95			# Gridding parameters
96			&PARM04
97			usingCartesianGrid=.TRUE.,
98			delR= 40*5.0,
99			dxSpacing=1.E3,
100			dySpacing=1.E3,
101			&
102
103			# Input datasets
104			&PARM05
105			bathyFile='topo_sl.bin',
106			hydrogThetaFile='thetaInitial.bin',
107			hydrogSaltFile='spiceInitial.bin',
108			&