Annotation of /MITgcm_contrib/high_res_cube/input/data

# ====================
# | 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.E-3,
 no_slip_sides=.TRUE.,
 no_slip_bottom=.TRUE.,
 diffKhT=400.,
# diffK4T=1.E10,
 diffKrT=3.E-5,
# diffK4S=1.E10,
 diffKhS=400.,
 diffKrS=3.E-5,
 tAlpha=2.E-4,
 sBeta =7.4E-4,
 gravity=9.81,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 eosType='JMD95Z',
 readBinaryPrec=32,
 writeBinaryPrec=32,
 useSingleCPUio=.TRUE.
 hFacMin=.3,
 hFacMinDr=50.,
 implicitDiffusion=.true.,
 implicitViscosity=.true.,
 allowFreezing=.true.,
# viscA4 = 1.e10,
 viscAh = 400
# cosPower = 3.,
 cosPower = 2.
 bottomDragQuadratic = 0.002,
 tempAdvScheme=4
 saltAdvScheme=4
 multiDimAdvection=.FALSE.,
 ivdc_kappa=10.,
 vectorInvariantMomentum=.TRUE.,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
# debugMode=.TRUE.
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=600,
 cg2dTargetResidual=1.E-12,
 &

# Time stepping parameters
 &PARM03
 niter0=0,
 nTimeSteps=20,
 abEps=0.1,
 deltaT=20.,
 cAdjFreq=0.,
 pChkptFreq=864000.0,
 chkptFreq=0.0,
 dumpFreq=20.0,
# tauThetaClimRelax = 5184000.,
# tauSaltClimRelax = 62208000.,
 monitorFreq=20.,
 &

# 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      ='BATHY_18360x85_filled.bin',
 hydrogThetaFile='LEVT01_JAN_18360x85x50.bin',
 hydrogSaltFile ='LEVS01_JAN_18360x85x50.bin',
#zonalWindFile  ='trenberth_taux.bin',
#meridWindFile  ='trenberth_tauy.bin',
#thetaClimFile  ='lev_surfT_cs_12m.bin',
#saltClimFile   ='lev_surfS_cs_12m.bin',
#surfQFile      ='shiQnet_cs32.bin',
#EmPmRFile      ='shiEmPR_cs32.bin',
 &
1	dimitri	1.1	# ====================
2			# \| Model parameters \|
3			# ====================
4			#
5			# Continuous equation parameters
6			&PARM01
7	dimitri	1.4	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=1.E-3,
11			no_slip_sides=.TRUE.,
12			no_slip_bottom=.TRUE.,
13			diffKhT=400.,
14			# diffK4T=1.E10,
15	dimitri	1.1	diffKrT=3.E-5,
16	dimitri	1.4	# diffK4S=1.E10,
17			diffKhS=400.,
18	dimitri	1.1	diffKrS=3.E-5,
19	dimitri	1.4	tAlpha=2.E-4,
20			sBeta =7.4E-4,
21	dimitri	1.1	gravity=9.81,
22	dimitri	1.4	rigidLid=.FALSE.,
23			implicitFreeSurface=.TRUE.,
24	dimitri	1.1	eosType='JMD95Z',
25	dimitri	1.4	readBinaryPrec=32,
26			writeBinaryPrec=32,
27			useSingleCPUio=.TRUE.
28			hFacMin=.3,
29			hFacMinDr=50.,
30			implicitDiffusion=.true.,
31			implicitViscosity=.true.,
32			allowFreezing=.true.,
33			# viscA4 = 1.e10,
34			viscAh = 400
35			# cosPower = 3.,
36			cosPower = 2.
37			bottomDragQuadratic = 0.002,
38			tempAdvScheme=4
39			saltAdvScheme=4
40			multiDimAdvection=.FALSE.,
41			ivdc_kappa=10.,
42	dimitri	1.1	vectorInvariantMomentum=.TRUE.,
43	dimitri	1.4	rigidLid=.FALSE.,
44	dimitri	1.1	implicitFreeSurface=.TRUE.,
45			exactConserv=.TRUE.,
46	dimitri	1.4	# debugMode=.TRUE.
47	dimitri	1.1	&
48
49			# Elliptic solver parameters
50			&PARM02
51	dimitri	1.4	cg2dMaxIters=600,
52			cg2dTargetResidual=1.E-12,
53	dimitri	1.1	&
54
55			# Time stepping parameters
56			&PARM03
57	dimitri	1.4	niter0=0,
58			nTimeSteps=20,
59			abEps=0.1,
60			deltaT=20.,
61			cAdjFreq=0.,
62			pChkptFreq=864000.0,
63			chkptFreq=0.0,
64			dumpFreq=20.0,
65			# tauThetaClimRelax = 5184000.,
66			# tauSaltClimRelax = 62208000.,
67			monitorFreq=20.,
68	dimitri	1.1	&
69
70			# Gridding parameters
71			&PARM04
72			usingCartesianGrid=.FALSE.,
73			usingSphericalPolarGrid=.FALSE.,
74			usingCurvilinearGrid=.TRUE.,
75	dimitri	1.4	delR = 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
76			10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04 , 19.82, 24.85,
77			31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18,
78			93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,
79			139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,
80			341.50,364.50,387.50,410.50,433.50,456.50,
81	dimitri	1.1	&
82
83			# Input datasets
84			&PARM05
85	dimitri	1.4	bathyFile ='BATHY_18360x85_filled.bin',
86			hydrogThetaFile='LEVT01_JAN_18360x85x50.bin',
87			hydrogSaltFile ='LEVS01_JAN_18360x85x50.bin',
88			#zonalWindFile ='trenberth_taux.bin',
89			#meridWindFile ='trenberth_tauy.bin',
90			#thetaClimFile ='lev_surfT_cs_12m.bin',
91			#saltClimFile ='lev_surfS_cs_12m.bin',
92			#surfQFile ='shiQnet_cs32.bin',
93			#EmPmRFile ='shiEmPR_cs32.bin',
94	dimitri	1.1	&