Contents of /MITgcm_contrib/ksnow/press_release/input/data

#:::====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
# momViscosity=.FALSE.,
# momForcing=.FALSE.,
# momAdvection=.False.,
# momstepping=.FALSE.,
 tempstepping=.FALSE.,
 saltstepping=.FALSE.,
 Tref = 50*-1.9,
 Sref = 50*34.4,
 viscAr=2.E-3,
#viscAh=600.0,
 viscAhGrid=0.4,
# viscA4Grid=0.02,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.TRUE.,
# bottomDragQuadratic=2.5E-3,
 diffKhT= 100.0,
#diffKhT= 600.0,
 diffKrT=5.E-5,
 diffKhS= 100.0,
#diffKhS= 600.0,
 diffKrS=5.E-5,
 bottomDragQuadratic=0.,
 staggerTimestep = .true.,
 tempAdvScheme=30,
 saltAdvScheme=30,
 vectorInvariantMomentum = .true.,
### momImplVertAdv=.true.,
 tempImplVertAdv=.true.,
 saltImplVertAdv=.true.,
 eosType='JMD95Z',
#tAlpha = 2.E-4,
#sBeta  = 7.4E-4,
 HeatCapacity_cp = 3974.0,
 rhoConst=1030.,
 rhoNil=1030.,
 gravity=9.81,
 gBaro=9.81,
 implicitDiffusion = .true.,
 implicitViscosity = .true.,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 useRealFreshWaterFlux = .true.,
 exactConserv=.TRUE.,
# start nlfs with rstar
#select_rStar=2,
 nonlinFreeSurf=4,
 hFacInf=0.2,
 hFacSup=2.0,
# end nlfs with rstar
 hFacMin=0.2,
#integr_GeoPot = 1,
#ivdc_kappa = 0.005.,
#useCDscheme=.true.,
#allowFreezing = .true.,
 readBinaryPrec=64,
#writeBinaryPrec=64,
 debuglevel = 0,
 f0 = 0.0,
 useSingleCpuIO=.TRUE.,
 globalFiles=.TRUE.,
 bottomVisc_pCell=.TRUE.,
 selectCoriMap=0,
 f0=-1.e-4
 selectImplicitDrag=2
 selectBotDragQuadr=2
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=600,
 cg2dTargetResidual=1.E-12,
 cg2dMinColumnEps=1,
 pReleaseVisc=8,
 printResidualFreq=200,
 &

 &PARM03
 niter0=0,
# nTimeSteps=8640,
#startTime=86400.,
 endTime=93312000.,
 deltaT=100,
 forcing_In_AB=.false.,
 abEps=0.1,
 cAdjFreq = -1.,
 pChkptFreq=432000.0,
 chkptFreq=432000.,
 dumpFreq=216000
 monitorSelect=1,
 monitorFreq=43200.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 usingSphericalPolarGrid=.FALSE.,
#delR = 10., 10., 10., 10., 10., 10., 10., 10.01, 10.03, 10.11,
#10.32, 10.8, 11.76, 13.42, 16.04, 19.82, 24.85, 31.1, 38.42, 46.5,
#55., 63.5, 71.58, 78.9, 85.15, 90.18, 93.96, 96.58, 98.25, 99.25,
#100.01, 101.33, 104.56, 111.33, 122.83, 139.09, 158.94, 180.83, 203.55, 226.5,
#249.5, 272.5, 295.5, 318.5, 341.5, 364.5, 387.5, 410.5, 433.5, 456.5,
 delR=105*20.,
 delX=20*1600,
 delY=200*1600,
# delYFile='dlat.bin',
# delY=30*10e3,
# xgOrigin = -105.5,
# ygOrigin =  -75.4457,
 &

# Input datasets
 &PARM05
 bathyFile='bathyS1.box',
 hydrogThetaFile='theta.init',  
 hydrogSaltFile='salt.init',
 pSurfInitFile='etainit.bin',
 &
1	#:::====================
2	# \| Model parameters \|
3	# ====================
4	#
5	# Continuous equation parameters
6	&PARM01
7	# momViscosity=.FALSE.,
8	# momForcing=.FALSE.,
9	# momAdvection=.False.,
10	# momstepping=.FALSE.,
11	tempstepping=.FALSE.,
12	saltstepping=.FALSE.,
13	Tref = 50*-1.9,
14	Sref = 50*34.4,
15	viscAr=2.E-3,
16	#viscAh=600.0,
17	viscAhGrid=0.4,
18	# viscA4Grid=0.02,
19	no_slip_sides=.FALSE.,
20	no_slip_bottom=.TRUE.,
21	# bottomDragQuadratic=2.5E-3,
22	diffKhT= 100.0,
23	#diffKhT= 600.0,
24	diffKrT=5.E-5,
25	diffKhS= 100.0,
26	#diffKhS= 600.0,
27	diffKrS=5.E-5,
28	bottomDragQuadratic=0.,
29	staggerTimestep = .true.,
30	tempAdvScheme=30,
31	saltAdvScheme=30,
32	vectorInvariantMomentum = .true.,
33	### momImplVertAdv=.true.,
34	tempImplVertAdv=.true.,
35	saltImplVertAdv=.true.,
36	eosType='JMD95Z',
37	#tAlpha = 2.E-4,
38	#sBeta = 7.4E-4,
39	HeatCapacity_cp = 3974.0,
40	rhoConst=1030.,
41	rhoNil=1030.,
42	gravity=9.81,
43	gBaro=9.81,
44	implicitDiffusion = .true.,
45	implicitViscosity = .true.,
46	rigidLid=.FALSE.,
47	implicitFreeSurface=.TRUE.,
48	useRealFreshWaterFlux = .true.,
49	exactConserv=.TRUE.,
50	# start nlfs with rstar
51	#select_rStar=2,
52	nonlinFreeSurf=4,
53	hFacInf=0.2,
54	hFacSup=2.0,
55	# end nlfs with rstar
56	hFacMin=0.2,
57	#integr_GeoPot = 1,
58	#ivdc_kappa = 0.005.,
59	#useCDscheme=.true.,
60	#allowFreezing = .true.,
61	readBinaryPrec=64,
62	#writeBinaryPrec=64,
63	debuglevel = 0,
64	f0 = 0.0,
65	useSingleCpuIO=.TRUE.,
66	globalFiles=.TRUE.,
67	bottomVisc_pCell=.TRUE.,
68	selectCoriMap=0,
69	f0=-1.e-4
70	selectImplicitDrag=2
71	selectBotDragQuadr=2
72	&
73
74	# Elliptic solver parameters
75	&PARM02
76	cg2dMaxIters=600,
77	cg2dTargetResidual=1.E-12,
78	cg2dMinColumnEps=1,
79	pReleaseVisc=8,
80	printResidualFreq=200,
81	&
82
83	&PARM03
84	niter0=0,
85	# nTimeSteps=8640,
86	#startTime=86400.,
87	endTime=93312000.,
88	deltaT=100,
89	forcing_In_AB=.false.,
90	abEps=0.1,
91	cAdjFreq = -1.,
92	pChkptFreq=432000.0,
93	chkptFreq=432000.,
94	dumpFreq=216000
95	monitorSelect=1,
96	monitorFreq=43200.,
97	&
98
99	# Gridding parameters
100	&PARM04
101	usingCartesianGrid=.TRUE.,
102	usingSphericalPolarGrid=.FALSE.,
103	#delR = 10., 10., 10., 10., 10., 10., 10., 10.01, 10.03, 10.11,
104	#10.32, 10.8, 11.76, 13.42, 16.04, 19.82, 24.85, 31.1, 38.42, 46.5,
105	#55., 63.5, 71.58, 78.9, 85.15, 90.18, 93.96, 96.58, 98.25, 99.25,
106	#100.01, 101.33, 104.56, 111.33, 122.83, 139.09, 158.94, 180.83, 203.55, 226.5,
107	#249.5, 272.5, 295.5, 318.5, 341.5, 364.5, 387.5, 410.5, 433.5, 456.5,
108	delR=105*20.,
109	delX=20*1600,
110	delY=200*1600,
111	# delYFile='dlat.bin',
112	# delY=30*10e3,
113	# xgOrigin = -105.5,
114	# ygOrigin = -75.4457,
115	&
116
117	# Input datasets
118	&PARM05
119	bathyFile='bathyS1.box',
120	hydrogThetaFile='theta.init',
121	hydrogSaltFile='salt.init',
122	pSurfInitFile='etainit.bin',
123	&