tutorial_tracer_adjsens/input_ad.som81/data

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef= 16.0 , 15.2 , 14.5 , 13.9 , 13.3 ,
       12.4 , 11.3 ,  9.9 ,  8.4 ,  6.7 ,
        5.2 ,  3.8 ,  2.9 ,  2.3 ,  1.8 ,
        1.5 ,  1.1 ,  0.8 ,  0.66,  0.63,
 sRef= 34.65, 34.75, 34.82, 34.87, 34.90,
       34.90, 34.86, 34.78, 34.69, 34.60,
       34.58, 34.62, 34.68, 34.72, 34.73,
       34.74, 34.73, 34.73, 34.72, 34.72,
 viscAr=1.E-3,
 viscAh=5.E5,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.TRUE.,
 diffKhT=1.E3,
 diffKrT=3.E-5,
 diffKhS=1.E3,
 diffKrS=3.E-5,
 tAlpha=2.E-4,
 sBeta =7.4E-4,
 eosType='LINEAR',
 gravity=9.81,
 implicitDiffusion=.TRUE.,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 useCDscheme=.TRUE.,
#-nonlinear free surface
 hFacInf=0.2,
 hFacSup=2.0,
 nonlinFreeSurf=3,
 select_rStar=2,
#
 staggerTimeStep=.TRUE.,
 tempAdvScheme=80,
 saltAdvScheme=81,
#-I/O settings
 readBinaryPrec=32,
 &

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

# Time stepping parameters
 &PARM03
# startTime         = 49766400000.,
 startTime         = 0.,
# endTime           = 49767264000.,
 endTime           = 345600.,
 deltaTClock       = 86400.,
 deltaTFreeSurf    = 86400.,
 deltaTmom         = 2400.,
 tauCD             = 321428.,
 cAdjFreq          = -1.,
 abEps             = 0.1,
 pChkptFreq        = 86400000.,
 chkptFreq         = 86400000.,
 dumpFreq          = 86400.,
 monitorFreq       = 43200.,
 adjMonitorFreq    = 43200.,
#
 tauThetaClimRelax = 2592000.,
 tauSaltClimRelax  = 2592000.,
#tauTr1ClimRelax   = 31104000., <= moved to data.ptracers
 &

# Gridding parameters
 &PARM04
 usingSphericalPolarGrid=.TRUE.,
 delR=  50.,  50.,  55.,  60.,
        65.,  70.,  80.,  95.,
       120., 155., 200., 260.,
       320., 400., 480., 570.,
       655., 725., 775., 815.,
 ygOrigin=-80.,
 delY=40*4.,
 delX=90*4.,
 &

# Input datasets
 &PARM05
 hydrogThetaFile='theta.bin',
 hydrogSaltFile='salt.bin',
 bathyFile='topog.bin',
 zonalWindFile='windx.bin',
 meridWindFile='windy.bin',
 thetaClimFile='SST.bin'
 saltClimFile='SSS.bin'
 &
1	heimbach	1.1	# ====================
2			# \| Model parameters \|
3			# ====================
4			#
5			# Continuous equation parameters
6			&PARM01
7			tRef= 16.0 , 15.2 , 14.5 , 13.9 , 13.3 ,
8			12.4 , 11.3 , 9.9 , 8.4 , 6.7 ,
9			5.2 , 3.8 , 2.9 , 2.3 , 1.8 ,
10			1.5 , 1.1 , 0.8 , 0.66, 0.63,
11			sRef= 34.65, 34.75, 34.82, 34.87, 34.90,
12			34.90, 34.86, 34.78, 34.69, 34.60,
13			34.58, 34.62, 34.68, 34.72, 34.73,
14			34.74, 34.73, 34.73, 34.72, 34.72,
15	jmc	1.2	viscAr=1.E-3,
16	heimbach	1.1	viscAh=5.E5,
17			no_slip_sides=.FALSE.,
18			no_slip_bottom=.TRUE.,
19			diffKhT=1.E3,
20	jmc	1.2	diffKrT=3.E-5,
21	heimbach	1.1	diffKhS=1.E3,
22	jmc	1.2	diffKrS=3.E-5,
23	heimbach	1.1	tAlpha=2.E-4,
24			sBeta =7.4E-4,
25	jmc	1.2	eosType='LINEAR',
26	heimbach	1.1	gravity=9.81,
27	jmc	1.2	implicitDiffusion=.TRUE.,
28	heimbach	1.1	rigidLid=.FALSE.,
29			implicitFreeSurface=.TRUE.,
30			exactConserv=.TRUE.,
31			useCDscheme=.TRUE.,
32	jmc	1.2	#-nonlinear free surface
33	heimbach	1.1	hFacInf=0.2,
34			hFacSup=2.0,
35			nonlinFreeSurf=3,
36			select_rStar=2,
37			#
38	jmc	1.2	staggerTimeStep=.TRUE.,
39			tempAdvScheme=80,
40	heimbach	1.1	saltAdvScheme=81,
41	jmc	1.2	#-I/O settings
42			readBinaryPrec=32,
43	heimbach	1.1	&
44
45			# Elliptic solver parameters
46			&PARM02
47			cg2dMaxIters=100,
48			cg2dTargetResidual=1.E-17,
49			&
50
51			# Time stepping parameters
52			&PARM03
53			# startTime = 49766400000.,
54			startTime = 0.,
55			# endTime = 49767264000.,
56			endTime = 345600.,
57	jmc	1.2	deltaTClock = 86400.,
58			deltaTFreeSurf = 86400.,
59			deltaTmom = 2400.,
60			tauCD = 321428.,
61			cAdjFreq = -1.,
62			abEps = 0.1,
63	heimbach	1.1	pChkptFreq = 86400000.,
64			chkptFreq = 86400000.,
65			dumpFreq = 86400.,
66			monitorFreq = 43200.,
67			adjMonitorFreq = 43200.,
68			#
69	jmc	1.2	tauThetaClimRelax = 2592000.,
70			tauSaltClimRelax = 2592000.,
71	heimbach	1.1	#tauTr1ClimRelax = 31104000., <= moved to data.ptracers
72			&
73
74			# Gridding parameters
75			&PARM04
76			usingSphericalPolarGrid=.TRUE.,
77	jmc	1.2	delR= 50., 50., 55., 60.,
78			65., 70., 80., 95.,
79			120., 155., 200., 260.,
80			320., 400., 480., 570.,
81			655., 725., 775., 815.,
82	heimbach	1.1	ygOrigin=-80.,
83			delY=40*4.,
84			delX=90*4.,
85			&
86
87			# Input datasets
88			&PARM05
89			hydrogThetaFile='theta.bin',
90			hydrogSaltFile='salt.bin',
91			bathyFile='topog.bin',
92			zonalWindFile='windx.bin',
93			meridWindFile='windy.bin',
94			thetaClimFile='SST.bin'
95			saltClimFile='SSS.bin'
96			&