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