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,
 viscAr= 1.77e-5
 no_slip_sides=.TRUE.,
 no_slip_bottom=.TRUE.,
 diffKhT=0.,
 diffK4T=0.,
 diffKrT=1.51e-5
 diffKhS=0.,
 diffK4S=0.,
 diffKrS=1.51e-5
 tAlpha=2.E-4,
 sBeta =7.4E-4,
 gravity=9.81,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 eosType='JMD95Z',
 readBinaryPrec=32,
 writeBinaryPrec=32,
 useSingleCPUio=.TRUE.
 hFacMinDr=50.,
 hFacMin=0.1,
 hFacInf=0.1,
 hFacSup=5.,
# z* - begin
    select_rStar=2,
    nonlinFreeSurf=4,
# * - end
 implicitDiffusion=.true.,
 implicitViscosity=.true.,
 allowFreezing=.false.,
 viscA4     = 0.,
 viscA4Grid = 2e-3,
 viscAh     = 0.,
 bottomDragQuadratic = 0.002,
 tempAdvScheme=33
 saltAdvScheme=33
 StaggerTimeStep=.true.
 multiDimAdvection=.true.
 vectorInvariantMomentum=.TRUE.,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 debuglevel=0
 &

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

# Time stepping parameters
 &PARM03
 niter0=0,
 endTime=259200
 abEps=0.1,
 deltaT=300.,
 cAdjFreq=0.,
 pChkptFreq=259200
 tavefreq=259200
 tauSaltClimRelax = 3842380.,
 monitorFreq=1
 &

# 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_18360x85_ver5e_filled.bin',
 hydrogThetaFile='LEVT01_JAN_18360x85x50.bin',
 hydrogSaltFile ='LEVS01_JAN_18360x85x50.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.5	viscAr= 1.77e-5
11	dimitri	1.1	no_slip_sides=.TRUE.,
12			no_slip_bottom=.TRUE.,
13	dimitri	1.2	diffKhT=0.,
14	dimitri	1.3	diffK4T=0.,
15	dimitri	1.5	diffKrT=1.51e-5
16	dimitri	1.2	diffKhS=0.,
17	dimitri	1.3	diffK4S=0.,
18	dimitri	1.5	diffKrS=1.51e-5
19	dimitri	1.1	tAlpha=2.E-4,
20			sBeta =7.4E-4,
21			gravity=9.81,
22			rigidLid=.FALSE.,
23			implicitFreeSurface=.TRUE.,
24			eosType='JMD95Z',
25			readBinaryPrec=32,
26			writeBinaryPrec=32,
27			useSingleCPUio=.TRUE.
28			hFacMinDr=50.,
29	dimitri	1.4	hFacMin=0.1,
30			hFacInf=0.1,
31			hFacSup=5.,
32			# z* - begin
33			select_rStar=2,
34			nonlinFreeSurf=4,
35			# * - end
36	dimitri	1.1	implicitDiffusion=.true.,
37			implicitViscosity=.true.,
38	dimitri	1.4	allowFreezing=.false.,
39	dimitri	1.6	viscA4 = 0.,
40	dimitri	1.8	viscA4Grid = 2e-3,
41	dimitri	1.6	viscAh = 0.,
42	dimitri	1.1	bottomDragQuadratic = 0.002,
43	dimitri	1.2	tempAdvScheme=33
44			saltAdvScheme=33
45			StaggerTimeStep=.true.
46			multiDimAdvection=.true.
47	dimitri	1.1	vectorInvariantMomentum=.TRUE.,
48			rigidLid=.FALSE.,
49			implicitFreeSurface=.TRUE.,
50			exactConserv=.TRUE.,
51	dimitri	1.4	debuglevel=0
52	dimitri	1.1	&
53
54			# Elliptic solver parameters
55			&PARM02
56	dimitri	1.4	cg2dMaxIters=100,
57			cg2dTargetResidual=1.E-5,
58	dimitri	1.1	&
59
60			# Time stepping parameters
61			&PARM03
62			niter0=0,
63	dimitri	1.7	endTime=259200
64	dimitri	1.1	abEps=0.1,
65	dimitri	1.4	deltaT=300.,
66	dimitri	1.1	cAdjFreq=0.,
67	dimitri	1.4	pChkptFreq=259200
68			tavefreq=259200
69	dimitri	1.5	tauSaltClimRelax = 3842380.,
70	dimitri	1.4	monitorFreq=1
71	dimitri	1.1	&
72
73			# Gridding parameters
74			&PARM04
75			usingCartesianGrid=.FALSE.,
76			usingSphericalPolarGrid=.FALSE.,
77			usingCurvilinearGrid=.TRUE.,
78			delR = 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
79			10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04 , 19.82, 24.85,
80			31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18,
81			93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,
82			139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,
83			341.50,364.50,387.50,410.50,433.50,456.50,
84			&
85
86			# Input datasets
87			&PARM05
88	dimitri	1.5	bathyFile ='ETOPO2_18360x85_ver5e_filled.bin',
89	dimitri	1.1	hydrogThetaFile='LEVT01_JAN_18360x85x50.bin',
90			hydrogSaltFile ='LEVS01_JAN_18360x85x50.bin',
91			&