MLAdjust/input.0.leithD/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.,
 ntimesteps=12,
 deltaT=1200.,
 abEps=0.01,
 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-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	ntimesteps=12,
69	deltaT=1200.,
70	abEps=0.01,
71	pchkptFreq=7200.,
72	taveFreq=7200.,
73	dumpFreq=7200.,
74	diagFreq=7200.,
75	monitorFreq=1200.,
76	cadjfreq=0.,
77	/
78
79	# Gridding parameters
80	&PARM04
81	usingCartesianGrid=.TRUE.,
82	usingSphericalPolarGrid=.FALSE.,
83	delR= 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	5d0, 5d0, 5d0, 5d0, 5d0,
91	dxspacing=1d3,
92	dyspacing=1d3,
93	/
94
95	# Input datasets
96	&PARM05
97	bathyFile='topo_sl.bin',
98	hydrogThetaFile='thetaInitial.bin',
99	hydrogSaltFile='spiceInitial.bin',
100	/