Contents of /MITgcm_contrib/heimbach/SO4x2/input_ad_nodiva/data


# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 viscAz=1.E-3,
 viscAh=5.E5,
 viscA4=3.E4,
 diffKhT=1.E3,
 diffKzT=3.E-5,
 diffKhS=1.E3,
 diffKzS=3.E-5,
 beta=1.E-11,
 tAlpha=2.E-4,
 sBeta =7.4E-4,
 rhonil=1035.,
 gravity=9.81,
 eosType = 'JMD95Z',
 implicitDiffusion=.true.,
 implicitViscosity=.true.,
 implicitFreeSurface=.true.,
 allowFreezing=.TRUE.,
 useRealFreshWaterFlux=.TRUE.,
 useCDscheme=.FALSE.,
 useSingleCpuIO=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 staggerTimeStep=.TRUE.,
 multiDimAdvection=.FALSE.,
 tempAdvScheme=30,
 saltAdvScheme=30,
# useNHMTerms=.TRUE.,
# turn on lopped cells
 hFacMin=0.3,
 hFacMindz=50.,
# set precision of data files
 readBinaryPrec=32,
 writeBinaryPrec=32,
 globalfiles=.FALSE.,
 debugLevel=2,
 /

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

# Time stepping parameters
 &PARM03
 nIter0 =     0,
 nTimeSteps = 144,
# nTimeSteps = 17543,
 deltaTmom   = 3600.0,
 deltaTtracer= 3600.0,
 deltaTClock = 3600.0,
# tauCD       = 172800.,
# if you are using a version later than checkpoint45d on the main branch
# you can uncomment the following line and increase the time step 
# deltaTtracer and deltaTClock to 172800.0 as well to speed up the 
# asynchronous time stepping
# deltaTfreesurf = 172800.0,
 dumpInitAndLast=.TRUE.
 abEps = 0.1,
 pChkptFreq=  10454400.,
 dumpFreq=    0.,
 taveFreq=    0.,
# diagFreq=    2629830.,
 adjDumpFreq=  2613600.,
 monitorFreq= 2613600.,
# monitorFreq=1.,
# 2 months restoring timescale for temperature
# tauThetaClimRelax =  5184000.0,
# 6 months restoring timescale for salinity
# tauSaltClimRelax = 15552000.0,
 /

# Gridding parameters
 &PARM04
 usingCartesianGrid=.FALSE.,
 usingSphericalPolarGrid=.TRUE.,
 delX=90*4.0,
 delY=30*2.0,
 delZ= 10.0, 10.0, 15.0, 20.0, 20.0, 25.0, 35.0, 50.0, 
       75.0,100.0,150.0,200.0,275.0,350.0,415.0,450.0,
      500.0,500.0,500.0,500.0,500.0,500.0,500.0
 phiMin=-80.,
 thetamin=0.,
 /

# Input datasets
 &PARM05
 bathyFile='bathy_SO_4x2_smooth.bin',
 uVelInitFile='Uini_SO_4x2.bin',
 vVelInitFile='Vini_SO_4x2.bin',
 hydrogThetaFile='Tini_SO_4x2.bin',
 hydrogSaltFile='Sini_SO_4x2.bin',
 /
1
2	# ====================
3	# \| Model parameters \|
4	# ====================
5	#
6	# Continuous equation parameters
7	&PARM01
8	viscAz=1.E-3,
9	viscAh=5.E5,
10	viscA4=3.E4,
11	diffKhT=1.E3,
12	diffKzT=3.E-5,
13	diffKhS=1.E3,
14	diffKzS=3.E-5,
15	beta=1.E-11,
16	tAlpha=2.E-4,
17	sBeta =7.4E-4,
18	rhonil=1035.,
19	gravity=9.81,
20	eosType = 'JMD95Z',
21	implicitDiffusion=.true.,
22	implicitViscosity=.true.,
23	implicitFreeSurface=.true.,
24	allowFreezing=.TRUE.,
25	useRealFreshWaterFlux=.TRUE.,
26	useCDscheme=.FALSE.,
27	useSingleCpuIO=.TRUE.,
28	vectorInvariantMomentum=.TRUE.,
29	staggerTimeStep=.TRUE.,
30	multiDimAdvection=.FALSE.,
31	tempAdvScheme=30,
32	saltAdvScheme=30,
33	# useNHMTerms=.TRUE.,
34	# turn on lopped cells
35	hFacMin=0.3,
36	hFacMindz=50.,
37	# set precision of data files
38	readBinaryPrec=32,
39	writeBinaryPrec=32,
40	globalfiles=.FALSE.,
41	debugLevel=2,
42	/
43
44	# Elliptic solver parameters
45	&PARM02
46	cg2dMaxIters=100,
47	cg2dTargetResidual=1.E-13,
48	/
49
50	# Time stepping parameters
51	&PARM03
52	nIter0 = 0,
53	nTimeSteps = 144,
54	# nTimeSteps = 17543,
55	deltaTmom = 3600.0,
56	deltaTtracer= 3600.0,
57	deltaTClock = 3600.0,
58	# tauCD = 172800.,
59	# if you are using a version later than checkpoint45d on the main branch
60	# you can uncomment the following line and increase the time step
61	# deltaTtracer and deltaTClock to 172800.0 as well to speed up the
62	# asynchronous time stepping
63	# deltaTfreesurf = 172800.0,
64	dumpInitAndLast=.TRUE.
65	abEps = 0.1,
66	pChkptFreq= 10454400.,
67	dumpFreq= 0.,
68	taveFreq= 0.,
69	# diagFreq= 2629830.,
70	adjDumpFreq= 2613600.,
71	monitorFreq= 2613600.,
72	# monitorFreq=1.,
73	# 2 months restoring timescale for temperature
74	# tauThetaClimRelax = 5184000.0,
75	# 6 months restoring timescale for salinity
76	# tauSaltClimRelax = 15552000.0,
77	/
78
79	# Gridding parameters
80	&PARM04
81	usingCartesianGrid=.FALSE.,
82	usingSphericalPolarGrid=.TRUE.,
83	delX=90*4.0,
84	delY=30*2.0,
85	delZ= 10.0, 10.0, 15.0, 20.0, 20.0, 25.0, 35.0, 50.0,
86	75.0,100.0,150.0,200.0,275.0,350.0,415.0,450.0,
87	500.0,500.0,500.0,500.0,500.0,500.0,500.0
88	phiMin=-80.,
89	thetamin=0.,
90	/
91
92	# Input datasets
93	&PARM05
94	bathyFile='bathy_SO_4x2_smooth.bin',
95	uVelInitFile='Uini_SO_4x2.bin',
96	vVelInitFile='Vini_SO_4x2.bin',
97	hydrogThetaFile='Tini_SO_4x2.bin',
98	hydrogSaltFile='Sini_SO_4x2.bin',
99	/