Contents of /MITgcm_contrib/atnguyen/llc_90/aste_90x150x60/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,
# try noslip with harmonic
 no_slip_sides  = .TRUE.,
 no_slip_bottom = .TRUE.,
#
# viscAr=5.e-4,
# viscA4Grid = 0.2,
 viscAr=1.e-3,
 viscAh=1.E0,
 viscAhGrid = 2.E-2,
 diffKhT=1.E1,
 diffKhS=1.E1,
 diffKrT=1.E-5,
 diffKrS=1.E-5,
 bottomDragQuadratic = 0.001,
#when using ggl90
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 gravity=9.81,
 rhonil=1027.,
 rhoConstFresh=1000.,
 eosType='JMD95Z',
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 useSingleCpuIO=.TRUE.,
 tempAdvScheme=30,
 saltAdvScheme=30,
 tempVertAdvScheme=3,
 saltVertAdvScheme=3,
 tempImplVertAdv=.TRUE.,
 saltImplVertAdv=.TRUE.,
 staggerTimeStep=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 useJamartWetPoints=.TRUE.,
 hFacMin=.2,
 hFacMinDr=50.,
 useRealFreshWaterFlux=.TRUE.,
# Gael's flags:
# balanceSaltClimRelax=.TRUE.,
 select_rStar=2,
 nonlinFreeSurf=4,
 allowFreezing=.FALSE.,
 readBinaryPrec=32,
 writeBinaryPrec=32,
 debugLevel=1,
 /

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=300,
 cg2dTargetResWunit=1.E-10,
 /

# Time stepping parameters
 &PARM03
 pickupStrictlyMatch=.FALSE.,
 nIter0=1,
# nIter0=1440,
#365 days 24*365=8760
#20years
 nTimeSteps=175200,
 deltaT = 3600.,
 abEps = 0.1,
 forcing_In_AB=.FALSE.,
 momDissip_In_AB=.FALSE.,
#when using ab3:
 doAB_onGtGs=.FALSE.,
 alph_AB=0.5,
 beta_AB=0.281105,
 pChkptFreq  =31536000.0,
# pChkptFreq  =2635200.0,
# pChkptFreq  =878400.0,
 chkptFreq   =2635200.0,
# dumpFreq    = 86400.0,
# monitorFreq = 3600.0,
 monitorFreq = 21600.0,
 /

# Gridding parameters
 &PARM04
 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_eccollcaste_90x50_min2pts.bin',
 hydrogThetaFile='T_OWPv1_M_eccollc_90x50.bin',
 hydrogSaltFile ='S_OWPv1_M_eccollc_90x50.bin',
 viscA4Dfile    ='viscA4Dfld_eccollc_90x50.bin',
 viscA4Zfile    ='viscA4Zfld_eccollc_90x50.bin',
# diffKrFile     ='DIFFKR_2_20_1_lat6070_cube81_reduced.bin',
 /
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	# try noslip with harmonic
11	no_slip_sides = .TRUE.,
12	no_slip_bottom = .TRUE.,
13	#
14	# viscAr=5.e-4,
15	# viscA4Grid = 0.2,
16	viscAr=1.e-3,
17	viscAh=1.E0,
18	viscAhGrid = 2.E-2,
19	diffKhT=1.E1,
20	diffKhS=1.E1,
21	diffKrT=1.E-5,
22	diffKrS=1.E-5,
23	bottomDragQuadratic = 0.001,
24	#when using ggl90
25	ivdc_kappa=10.,
26	implicitDiffusion=.TRUE.,
27	implicitViscosity=.TRUE.,
28	gravity=9.81,
29	rhonil=1027.,
30	rhoConstFresh=1000.,
31	eosType='JMD95Z',
32	implicitFreeSurface=.TRUE.,
33	exactConserv=.TRUE.,
34	useSingleCpuIO=.TRUE.,
35	tempAdvScheme=30,
36	saltAdvScheme=30,
37	tempVertAdvScheme=3,
38	saltVertAdvScheme=3,
39	tempImplVertAdv=.TRUE.,
40	saltImplVertAdv=.TRUE.,
41	staggerTimeStep=.TRUE.,
42	vectorInvariantMomentum=.TRUE.,
43	useJamartWetPoints=.TRUE.,
44	hFacMin=.2,
45	hFacMinDr=50.,
46	useRealFreshWaterFlux=.TRUE.,
47	# Gael's flags:
48	# balanceSaltClimRelax=.TRUE.,
49	select_rStar=2,
50	nonlinFreeSurf=4,
51	allowFreezing=.FALSE.,
52	readBinaryPrec=32,
53	writeBinaryPrec=32,
54	debugLevel=1,
55	/
56
57	# Elliptic solver parameters
58	&PARM02
59	cg2dMaxIters=300,
60	cg2dTargetResWunit=1.E-10,
61	/
62
63	# Time stepping parameters
64	&PARM03
65	pickupStrictlyMatch=.FALSE.,
66	nIter0=1,
67	# nIter0=1440,
68	#365 days 24*365=8760
69	#20years
70	nTimeSteps=175200,
71	deltaT = 3600.,
72	abEps = 0.1,
73	forcing_In_AB=.FALSE.,
74	momDissip_In_AB=.FALSE.,
75	#when using ab3:
76	doAB_onGtGs=.FALSE.,
77	alph_AB=0.5,
78	beta_AB=0.281105,
79	pChkptFreq =31536000.0,
80	# pChkptFreq =2635200.0,
81	# pChkptFreq =878400.0,
82	chkptFreq =2635200.0,
83	# dumpFreq = 86400.0,
84	# monitorFreq = 3600.0,
85	monitorFreq = 21600.0,
86	/
87
88	# Gridding parameters
89	&PARM04
90	usingCurvilinearGrid=.TRUE.,
91	delR =
92	10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
93	10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04, 19.82, 24.85,
94	31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18,
95	93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,
96	139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,
97	341.50,364.50,387.50,410.50,433.50,456.50,
98	/
99
100	# Input datasets
101	&PARM05
102	bathyFile = 'bathy_eccollcaste_90x50_min2pts.bin',
103	hydrogThetaFile='T_OWPv1_M_eccollc_90x50.bin',
104	hydrogSaltFile ='S_OWPv1_M_eccollc_90x50.bin',
105	viscA4Dfile ='viscA4Dfld_eccollc_90x50.bin',
106	viscA4Zfile ='viscA4Zfld_eccollc_90x50.bin',
107	# diffKrFile ='DIFFKR_2_20_1_lat6070_cube81_reduced.bin',
108	/