MLAdjust/input.0.smag/data

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 viscAr=5.E-3,
# This is the minimum viscosity used total=viscah+viscleith
# Often U*deltaX/2 for grid noise (Griffies and Hallberg,2000)
# viscAh=3.E1,
# This is the standard Leith Viscosity, prop to grad(curl v_h)
# viscC2leith=5.,
# 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=5.,
# These are the minimum viscosities used with Smagorinsky
# total=viscastrain+viscsmag, etc.
 viscAstrain=3.E1,
 viscAtension=3.E1,
# This is harmonic Smagorinsky Coefficient
# 0.15=value from 3-d turbulence theory; 
# 3-4=stable oceanic values (Griffies and Hallberg, 2000)
 viscC2Smag=4.,
# This is the max Viscosity used, be it Smag or Leith, and it scales 
# with grid size and timestep
# Often 0.5 for stability (Griffies and Hallberg,2000)
 viscAhGridMax=0.5,
 useFullLeith=.TRUE.,
 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.,
 deltaT=1200.,
 abEps=0.01,
 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			viscAr=5.E-3,
8			# This is the minimum viscosity used total=viscah+viscleith
9			# Often U*deltaX/2 for grid noise (Griffies and Hallberg,2000)
10			# viscAh=3.E1,
11			# This is the standard Leith Viscosity, prop to grad(curl v_h)
12			# viscC2leith=5.,
13			# This is an additional Leith Viscosity, prop to grad(div.v_h)
14			# It traps instabilities that have no signal in grad(curl v_h)
15			# viscC2leithD=5.,
16			# These are the minimum viscosities used with Smagorinsky
17			# total=viscastrain+viscsmag, etc.
18			viscAstrain=3.E1,
19			viscAtension=3.E1,
20			# This is harmonic Smagorinsky Coefficient
21			# 0.15=value from 3-d turbulence theory;
22			# 3-4=stable oceanic values (Griffies and Hallberg, 2000)
23			viscC2Smag=4.,
24			# This is the max Viscosity used, be it Smag or Leith, and it scales
25			# with grid size and timestep
26			# Often 0.5 for stability (Griffies and Hallberg,2000)
27			viscAhGridMax=0.5,
28			useFullLeith=.TRUE.,
29			no_slip_sides=.FALSE.,
30			no_slip_bottom=.FALSE.,
31			ivdc_kappa=10.,
32			implicitDiffusion=.TRUE.,
33			implicitViscosity=.TRUE.,
34			f0=7.29e-5,
35			beta=0.E-11,
36			eosType='LINEAR',
37			# Using Temperature as only Active Tracer
38			tAlpha=2.E-4,
39			# Using Salinity as Spice
40			sBeta =0.E-4,
41			saltStepping=.TRUE.,
42			gravity=9.81,
43			rhonil=1035.,
44			implicitFreeSurface=.TRUE.,
45			usejamartwetpoints=.TRUE.,
46			rigidLid=.FALSE.,
47			staggerTimeStep=.TRUE.,
48			readBinaryPrec=64,
49			writeBinaryPrec=64,
50			tempAdvScheme=33,
51			saltAdvScheme=33,
52			nonHydrostatic=.FALSE.,
53			vectorInvariantMomentum=.TRUE.,
54			/
55
56			# Elliptic solver parameters
57			&PARM02
58			cg2dMaxIters=1000,
59			cg2dTargetResidual=1.E-11,
60			cg3dMaxIters=40,
61			cg3dTargetResidual=1.E-9,
62			/
63
64			# Time stepping parameters
65			&PARM03
66			niter0=0,
67			endTime=86400.,
68			deltaT=1200.,
69			abEps=0.01,
70			pchkptFreq=43200.,
71			taveFreq=7200.,
72			dumpFreq=7200.,
73			diagFreq=7200.,
74			monitorFreq=1200.,
75			cadjfreq=0.,
76			/
77
78			# Gridding parameters
79			&PARM04
80			usingCartesianGrid=.TRUE.,
81			usingSphericalPolarGrid=.FALSE.,
82			delR= 5d0, 5d0, 5d0, 5d0, 5d0,
83			5d0, 5d0, 5d0, 5d0, 5d0,
84			5d0, 5d0, 5d0, 5d0, 5d0,
85			5d0, 5d0, 5d0, 5d0, 5d0,
86			5d0, 5d0, 5d0, 5d0, 5d0,
87			5d0, 5d0, 5d0, 5d0, 5d0,
88			5d0, 5d0, 5d0, 5d0, 5d0,
89			5d0, 5d0, 5d0, 5d0, 5d0,
90			dxspacing=1d3,
91			dyspacing=1d3,
92			/
93
94			# Input datasets
95			&PARM05
96			bathyFile='topo_sl.bin',
97			hydrogThetaFile='thetaInitial.bin',
98			hydrogSaltFile='spiceInitial.bin',
99			/