atm_gray/inp_cs/data.mk_ini

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef=295.2, 295.5, 295.9, 296.3, 296.7, 297.1, 297.6, 298.1, 298.7, 299.3,
      300.0, 300.7, 301.9, 304.1, 308.0, 312.0, 316.0, 320.,   324., 329.,
      338.,  339., 362.3, 419.2, 573.8,
 sRef=0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 15*0.,
#sRef=25*0.,
 viscAr=0.E1,
 viscAh=0.E6,
 viscA4=0.E17,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.FALSE.,
 diffKhT=0.E3,
 diffKrT=0.,
 diffK4T=0.E17,
 diffKrS=0.E2,
 diffKhS=0.E3,
 diffK4S=0.E17,
 buoyancyRelation='ATMOSPHERIC',
 eosType='IDEALG',
#atm_Rq=0.6078,
 Integr_GeoPot=2,
 selectFindRoSurf=1,
 rotationPeriod=86400.,
 gravity=9.81,
 rhonil=1.0,
 rhoConst=1.0,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 select_rStar=2,
 nonlinFreeSurf=4,
 hFacInf=0.2,
 hFacSup=2.0,
 hFacMin=1.0,
 uniformLin_PhiSurf=.FALSE.,
#vectorInvariantMomentum=.TRUE., 
#tempAdvScheme=77,
 saltAdvScheme=77,
 staggerTimeStep=.TRUE.,
 readBinaryPrec=64,
 writeBinaryPrec=64,
 globalFiles=.TRUE.,
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=200,
 cg2dTargetResWunit=1.E-17,
 &

# Time stepping parameters
 &PARM03
 nIter0=0,
 nTimeSteps=0,
#nTimeSteps=2880,
#endTime=31104000.,
#endTime=5184000.,
#endTime=864000.,
 deltaT=300.,
 abEps=0.1,
 cAdjFreq=0.,
#chkptFreq =43200.,
#pChkptFreq=54000.,
#pChkptFreq=1200000.0,
 dumpFreq=2592000.0,
 monitorFreq=86400.,
#monitorFreq=3600.,
#dumpFreq   =1.,
#monitorFreq=1.,
 &

# Gridding parameters
 &PARM04
 usingCurvilinearGrid=.TRUE.,
 horizGridFile='dxC1_dXYa',
 radius_fromHorizGrid=6370.E3,
 delR=25*40.E2,
 Ro_SeaLevel=1.E5,
#rSphere=6370.E3,
 &

# Input datasets
 &PARM05
#hydrogThetaFile='ini_theta.bin',
 &
1	jmc	1.1	# ====================
2			# \| Model parameters \|
3			# ====================
4			#
5			# Continuous equation parameters
6			&PARM01
7			tRef=295.2, 295.5, 295.9, 296.3, 296.7, 297.1, 297.6, 298.1, 298.7, 299.3,
8			300.0, 300.7, 301.9, 304.1, 308.0, 312.0, 316.0, 320., 324., 329.,
9			338., 339., 362.3, 419.2, 573.8,
10			sRef=0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 15*0.,
11			#sRef=25*0.,
12			viscAr=0.E1,
13			viscAh=0.E6,
14			viscA4=0.E17,
15			no_slip_sides=.FALSE.,
16			no_slip_bottom=.FALSE.,
17			diffKhT=0.E3,
18			diffKrT=0.,
19			diffK4T=0.E17,
20			diffKrS=0.E2,
21			diffKhS=0.E3,
22			diffK4S=0.E17,
23			buoyancyRelation='ATMOSPHERIC',
24			eosType='IDEALG',
25			#atm_Rq=0.6078,
26			Integr_GeoPot=2,
27			selectFindRoSurf=1,
28			rotationPeriod=86400.,
29			gravity=9.81,
30			rhonil=1.0,
31			rhoConst=1.0,
32			rigidLid=.FALSE.,
33			implicitFreeSurface=.TRUE.,
34			exactConserv=.TRUE.,
35			select_rStar=2,
36			nonlinFreeSurf=4,
37			hFacInf=0.2,
38			hFacSup=2.0,
39			hFacMin=1.0,
40			uniformLin_PhiSurf=.FALSE.,
41			#vectorInvariantMomentum=.TRUE.,
42			#tempAdvScheme=77,
43			saltAdvScheme=77,
44			staggerTimeStep=.TRUE.,
45			readBinaryPrec=64,
46			writeBinaryPrec=64,
47			globalFiles=.TRUE.,
48			&
49
50			# Elliptic solver parameters
51			&PARM02
52			cg2dMaxIters=200,
53			cg2dTargetResWunit=1.E-17,
54			&
55
56			# Time stepping parameters
57			&PARM03
58			nIter0=0,
59			nTimeSteps=0,
60			#nTimeSteps=2880,
61			#endTime=31104000.,
62			#endTime=5184000.,
63			#endTime=864000.,
64			deltaT=300.,
65			abEps=0.1,
66			cAdjFreq=0.,
67			#chkptFreq =43200.,
68			#pChkptFreq=54000.,
69			#pChkptFreq=1200000.0,
70			dumpFreq=2592000.0,
71			monitorFreq=86400.,
72			#monitorFreq=3600.,
73			#dumpFreq =1.,
74			#monitorFreq=1.,
75			&
76
77			# Gridding parameters
78			&PARM04
79			usingCurvilinearGrid=.TRUE.,
80			horizGridFile='dxC1_dXYa',
81			radius_fromHorizGrid=6370.E3,
82			delR=25*40.E2,
83			Ro_SeaLevel=1.E5,
84			#rSphere=6370.E3,
85			&
86
87			# Input datasets
88			&PARM05
89			#hydrogThetaFile='ini_theta.bin',
90			&