Contents of /MITgcm_contrib/high_res_cube/input/data_hr

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef               = 3*23.,3*22.,21.,2*20.,19.,2*18.,17.,2*16.,15.,14.,13.,
                      12.,11.,2*9.,8.,7.,2*6.,2*5.,3*4.,3*3.,4*2.,12*1.,
 sRef               = 50*34.5,
 viscAr= 6.0826e-04,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.TRUE.,
 diffKhT=0.,
 diffK4T=0.,
 diffKrT=0.,
 diffKhS=0.,
 diffK4S=0.,
 diffKrS=0.,
# Bryan and Lewis 1979 vertical diffusivity parameters
# with latitudinal dependence
 diffKrBL79surf      = 3.0e-5,
 diffKrBL79deep      = 1.3e-4,
 diffKrBL79scl       = 200.,
 diffKrBL79Ho        = -2000.,
 BL79LatVary         = 30.,
 diffKrBLEQsurf      = 8.0e-6,
 diffKrBLEQdeep      = 1.3e-4,
 diffKrBLEQscl       = 200.,
 diffKrBLEQHo        = -2000.,
#
 tAlpha=2.E-4,
 sBeta =7.4E-4,
 gravity=9.81,
 eosType='JMD95Z',
 readBinaryPrec=32,
 writeBinaryPrec=32,
 useSingleCPUio=.TRUE.,
 hFacMinDr=50.,
 hFacMin=0.3,
 hFacInf=0.1,
 hFacSup=5.,
# z* - begin
    select_rStar=2,
    nonlinFreeSurf=4,
# * - end
 implicitDiffusion=.true.,
 implicitViscosity=.true.,
 allowFreezing=.false.,
 viscAh=0.,
 viscAhGrid=0.,
# Following change took place starting in checkpoint57v_post
# Redefine viscC?Leith*.  Change the coefficients of Leith viscosity 
# to be proportional to enstrophy dissipation length.  Thus,
#  viscC?leith?_new=pi*(viscC?leith?_old)^(1/3)
# This is now consistent with Smagorinsky and with the manual.
# viscC4Leith=1.0,
# viscC4Leithd=1.0,
# viscC4Leith=3.1416,
# viscC4Leithd=3.1416,
 viscC4Leith=1.5,
 viscC4Leithd=1.5,
 viscA4GridMax=0.5,
##########################################################
 viscC2Leith=0.,
 viscA4=0.,
 useAreaViscLength=.TRUE.,
 sideDragFactor=0.,
 SadournyCoriolis=.FALSE.,
 useAbsVorticity=.FALSE.,
 highOrderVorticity  = .TRUE.,
 bottomDragQuadratic = 0.002,
 tempAdvScheme=33,
 saltAdvScheme=33,
 StaggerTimeStep=.true.,
 multiDimAdvection=.true.,
 vectorInvariantMomentum=.TRUE.,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 debuglevel=-1,
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=100,
 cg2dTargetResidual=1.E-5,
 &

# Time stepping parameters
 &PARM03
 niter0=216,
# endtime for December 31, 2002
 endtime=347155200.,
 abEps=0.1,
 deltaT=1200,
 cAdjFreq=0.,
 pChkptFreq=31536000.,
 chkptFreq=2592000.,
# blending: 12/24/06
  tauSaltClimRelax = 8722203.,
# tauSaltClimRelax = 3842380.,
 monitorFreq=86400.,
 forcing_In_AB = .FALSE.,
 dumpInitAndLast=.FALSE.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.FALSE.,
 usingSphericalPolarGrid=.FALSE.,
 usingCurvilinearGrid=.TRUE.,
 delR   = 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
 10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04 , 19.82, 24.85, 
 31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 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.50,249.50,272.50,295.50,318.50,
 341.50,364.50,387.50,410.50,433.50,456.50,
 &

# Input datasets
 &PARM05
 bathyFile      ='ETOPO2_510x6x510_ver5e_filled.bin',
 hydrogThetaFile='WGHC_T_510x6x510.bin',
 hydrogSaltFile ='WGHC_S_510x6x510.bin',
 &
1	# ====================
2	# \| Model parameters \|
3	# ====================
4	#
5	# Continuous equation parameters
6	&PARM01
7	tRef = 323.,322.,21.,220.,19.,218.,17.,2*16.,15.,14.,13.,
8	12.,11.,29.,8.,7.,26.,25.,34.,33.,42.,12*1.,
9	sRef = 50*34.5,
10	viscAr= 6.0826e-04,
11	no_slip_sides=.FALSE.,
12	no_slip_bottom=.TRUE.,
13	diffKhT=0.,
14	diffK4T=0.,
15	diffKrT=0.,
16	diffKhS=0.,
17	diffK4S=0.,
18	diffKrS=0.,
19	# Bryan and Lewis 1979 vertical diffusivity parameters
20	# with latitudinal dependence
21	diffKrBL79surf = 3.0e-5,
22	diffKrBL79deep = 1.3e-4,
23	diffKrBL79scl = 200.,
24	diffKrBL79Ho = -2000.,
25	BL79LatVary = 30.,
26	diffKrBLEQsurf = 8.0e-6,
27	diffKrBLEQdeep = 1.3e-4,
28	diffKrBLEQscl = 200.,
29	diffKrBLEQHo = -2000.,
30	#
31	tAlpha=2.E-4,
32	sBeta =7.4E-4,
33	gravity=9.81,
34	eosType='JMD95Z',
35	readBinaryPrec=32,
36	writeBinaryPrec=32,
37	useSingleCPUio=.TRUE.,
38	hFacMinDr=50.,
39	hFacMin=0.3,
40	hFacInf=0.1,
41	hFacSup=5.,
42	# z* - begin
43	select_rStar=2,
44	nonlinFreeSurf=4,
45	# * - end
46	implicitDiffusion=.true.,
47	implicitViscosity=.true.,
48	allowFreezing=.false.,
49	viscAh=0.,
50	viscAhGrid=0.,
51	# Following change took place starting in checkpoint57v_post
52	# Redefine viscC?Leith*. Change the coefficients of Leith viscosity
53	# to be proportional to enstrophy dissipation length. Thus,
54	# viscC?leith?_new=pi*(viscC?leith?_old)^(1/3)
55	# This is now consistent with Smagorinsky and with the manual.
56	# viscC4Leith=1.0,
57	# viscC4Leithd=1.0,
58	# viscC4Leith=3.1416,
59	# viscC4Leithd=3.1416,
60	viscC4Leith=1.5,
61	viscC4Leithd=1.5,
62	viscA4GridMax=0.5,
63	##########################################################
64	viscC2Leith=0.,
65	viscA4=0.,
66	useAreaViscLength=.TRUE.,
67	sideDragFactor=0.,
68	SadournyCoriolis=.FALSE.,
69	useAbsVorticity=.FALSE.,
70	highOrderVorticity = .TRUE.,
71	bottomDragQuadratic = 0.002,
72	tempAdvScheme=33,
73	saltAdvScheme=33,
74	StaggerTimeStep=.true.,
75	multiDimAdvection=.true.,
76	vectorInvariantMomentum=.TRUE.,
77	rigidLid=.FALSE.,
78	implicitFreeSurface=.TRUE.,
79	exactConserv=.TRUE.,
80	debuglevel=-1,
81	&
82
83	# Elliptic solver parameters
84	&PARM02
85	cg2dMaxIters=100,
86	cg2dTargetResidual=1.E-5,
87	&
88
89	# Time stepping parameters
90	&PARM03
91	niter0=216,
92	# endtime for December 31, 2002
93	endtime=347155200.,
94	abEps=0.1,
95	deltaT=1200,
96	cAdjFreq=0.,
97	pChkptFreq=31536000.,
98	chkptFreq=2592000.,
99	# blending: 12/24/06
100	tauSaltClimRelax = 8722203.,
101	# tauSaltClimRelax = 3842380.,
102	monitorFreq=86400.,
103	forcing_In_AB = .FALSE.,
104	dumpInitAndLast=.FALSE.,
105	&
106
107	# Gridding parameters
108	&PARM04
109	usingCartesianGrid=.FALSE.,
110	usingSphericalPolarGrid=.FALSE.,
111	usingCurvilinearGrid=.TRUE.,
112	delR = 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
113	10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04 , 19.82, 24.85,
114	31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18,
115	93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,
116	139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,
117	341.50,364.50,387.50,410.50,433.50,456.50,
118	&
119
120	# Input datasets
121	&PARM05
122	bathyFile ='ETOPO2_510x6x510_ver5e_filled.bin',
123	hydrogThetaFile='WGHC_T_510x6x510.bin',
124	hydrogSaltFile ='WGHC_S_510x6x510.bin',
125	&