MLAdjust/input.1.leith/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.85 smag is approximately the same as 1.85 leith+leithd
# 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=36,
#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	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.85 smag is approximately the same as 1.85 leith+leithd
36	# viscC2Smag=1.0,
37	viscC4Smag=1.0,
38	# This is the max Viscosity used, be it Smag or Leith, and it scales
39	# with grid size and timestep
40	# viscAhGridMax=1.0,
41	viscA4GridMax=1.0,
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	ivdc_kappa=10.,
49	implicitDiffusion=.TRUE.,
50	implicitViscosity=.TRUE.,
51	f0=7.29e-5,
52	beta=0.E-11,
53	eosType='LINEAR',
54	# Using Temperature as only Active Tracer
55	tAlpha=2.E-4,
56	# Using Salinity as Spice
57	sBeta =0.E-4,
58	saltStepping=.TRUE.,
59	gravity=9.81,
60	rhonil=1035.,
61	implicitFreeSurface=.TRUE.,
62	usejamartwetpoints=.TRUE.,
63	rigidLid=.FALSE.,
64	staggerTimeStep=.TRUE.,
65	readBinaryPrec=64,
66	writeBinaryPrec=64,
67	tempAdvScheme=33,
68	saltAdvScheme=33,
69	nonHydrostatic=.FALSE.,
70	vectorInvariantMomentum=.TRUE.,
71	/
72
73	# Elliptic solver parameters
74	&PARM02
75	cg2dMaxIters=1000,
76	cg2dTargetResidual=1.E-11,
77	cg3dMaxIters=40,
78	cg3dTargetResidual=1.E-9,
79	/
80
81	# Time stepping parameters
82	&PARM03
83	niter0=36,
84	#endTime=86400.,
85	#- for testing, only run for 12 iterations:
86	nTimeSteps=12,
87	deltaT=1200.,
88	abEps=0.01,
89	#- for testing, reduce amount of output by commenting out output frequencies:
90	# pchkptFreq=43200.,
91	#taveFreq=7200.,
92	#dumpFreq=7200.,
93	#diagFreq=7200.,
94	monitorFreq=1200.,
95	cadjfreq=0.,
96	/
97
98	# Gridding parameters
99	&PARM04
100	usingCartesianGrid=.TRUE.,
101	usingSphericalPolarGrid=.FALSE.,
102	delR= 5d0, 5d0, 5d0, 5d0, 5d0,
103	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	dxspacing=1d3,
111	dyspacing=1d3,
112	/
113
114	# Input datasets
115	&PARM05
116	bathyFile='topo_sl.bin',
117	hydrogThetaFile='thetaInitial.bin',
118	hydrogSaltFile='spiceInitial.bin',
119	/