MLAdjust/input.0.leith/data

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef(1)=29.,
 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=.TRUE.,
# This uses the Strain-Tension form for calculating viscous terms
 useStrainTensionVisc=.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.,
#- 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	# ====================
2	# \| Model parameters \|
3	# ====================
4	#
5	# Continuous equation parameters
6	&PARM01
7	tRef(1)=29.,
8	viscAr=5.E-5,
9	# This is the background viscosity used total=viscah+viscleith
10	# viscAh=1.E2,
11	# viscA4=1.E2,
12	# This is the background viscosity used on vorticity points
13	# viscAhZ=1.E2,
14	# viscA4Z=1.E2,
15	# This is the background viscosity used on divergence points
16	# viscAhD=1.E2,
17	# viscA4D=1.E2,
18	# viscAhd and viscahz are set to viscah if they are unset initially
19	#
20	# This is a maximum gridscale Reynolds number limiter
21	# ViscAhRemax=10,
22	# ViscA4Remax=10,
23	# This is a minimum autoscaling viscosity
24	# ViscAhGridMin=.01,
25	# ViscA4GridMin=.01,
26	# This is the standard Leith Viscosity, prop to grad(curl v_h)
27	viscC2leith=1.85,
28	# viscC4leith=1.85,
29	# This is an additional Leith Viscosity, prop to grad(div.v_h)
30	# It traps instabilities that have no signal in grad(curl v_h)
31	viscC2leithD=1.85,
32	# viscC4leithD=1.85,
33	# This is harmonic Smagorinsky Coefficient
34	# 0.2-0.9=value from 3-d turbulence theory; Smagorinsky (1993)
35	# 2.2-4=stable oceanic values (Griffies and Hallberg, 2000)
36	# 1.0 smag is approximately the same as 1.85 leith+leithd
37	# for this grid
38	viscC2Smag=1.0,
39	# viscC4Smag=1.0,
40	# This is the max Viscosity used, be it Smag or Leith, and it scales
41	# with grid size and timestep
42	viscAhGridMax=1.0,
43	# viscA4GridMax=1.0,
44	# This uses the full grad(vort) and grad(div)
45	useFullLeith=.TRUE.,
46	# This uses the Strain-Tension form for calculating viscous terms
47	useStrainTensionVisc=.TRUE.,
48	no_slip_sides=.FALSE.,
49	no_slip_bottom=.FALSE.,
50	ivdc_kappa=10.,
51	implicitDiffusion=.TRUE.,
52	implicitViscosity=.TRUE.,
53	f0=7.29e-5,
54	beta=0.E-11,
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	saltStepping=.TRUE.,
61	gravity=9.81,
62	rhonil=1035.,
63	implicitFreeSurface=.TRUE.,
64	usejamartwetpoints=.TRUE.,
65	rigidLid=.FALSE.,
66	staggerTimeStep=.TRUE.,
67	readBinaryPrec=64,
68	writeBinaryPrec=64,
69	tempAdvScheme=33,
70	saltAdvScheme=33,
71	nonHydrostatic=.FALSE.,
72	vectorInvariantMomentum=.TRUE.,
73	/
74
75	# Elliptic solver parameters
76	&PARM02
77	cg2dMaxIters=1000,
78	cg2dTargetResidual=1.E-11,
79	cg3dMaxIters=40,
80	cg3dTargetResidual=1.E-9,
81	/
82
83	# Time stepping parameters
84	&PARM03
85	niter0=0,
86	#endTime=86400.,
87	#- for testing, only run for 12 iterations:
88	nTimeSteps=12,
89	deltaT=1200.,
90	abEps=0.01,
91	#- for testing, reduce amount of output by commenting out output frequencies:
92	#pchkptFreq=43200.,
93	#taveFreq=7200.,
94	#dumpFreq=7200.,
95	#diagFreq=7200.,
96	monitorFreq=1200.,
97	cadjfreq=0.,
98	/
99
100	# Gridding parameters
101	&PARM04
102	usingCartesianGrid=.TRUE.,
103	usingSphericalPolarGrid=.FALSE.,
104	delR= 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	5d0, 5d0, 5d0, 5d0, 5d0,
112	dxspacing=1d3,
113	dyspacing=1d3,
114	/
115
116	# Input datasets
117	&PARM05
118	bathyFile='topo_sl.bin',
119	hydrogThetaFile='thetaInitial.bin',
120	hydrogSaltFile='spiceInitial.bin',
121	/