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	# ====================
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	&