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.,
 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	# ====================
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	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	/