global_oce_in_p/input/data.tex

\begin{verbatim}
     1  # ====================
     2  # | Model parameters |
     3  # ====================
     4  #
     5  # Continuous equation parameters
     6   &PARM01
     7   tRef=15*20.,
     8   sRef=15*35.,
     9   viscAh =3.E5,
    10   no_slip_sides=.TRUE.,
    11   viscAr =1.721611620915750e5,
    12  #viscAz =1.67E-3,
    13   no_slip_bottom=.FALSE.,
    14   diffKhT=1.E3, 
    15   diffKrT=5.154525811125000e3,
    16  #diffKzT=0.5E-4,
    17   diffKhS=1.E3,  
    18   diffKrS=5.154525811125000e3,
    19  #diffKzS=0.5E-4,
    20   implicitDiffusion=.TRUE.,
    21   ivdc_kappa=1.030905162225000e9,
    22  #ivdc_kappa=10.0,
    23   gravity=9.81,
    24   gravitySign=-1.D0,
    25   rhonil=1035.,
    26   buoyancyRelation='OCEANICP',
    27   eosType='JMD95P',
    28   rigidLid=.FALSE.,
    29   implicitFreeSurface=.TRUE.,
    30   exactConserv=.TRUE.,
    31   nonlinFreeSurf=3,
    32   hFacInf=0.2,
    33   hFacSup=2.0,
    34   useRealFreshWaterFlux=.FALSE.,
    35   readBinaryPrec=64,
    36   writeBinaryPrec=64,
    37   cosPower=0.5,
    38   &
    39  
    40  # Elliptic solver parameters
    41   &PARM02
    42   cg2dMaxIters=200,
    43   cg2dTargetResidual=1.E-9,
    44   &
    45  
    46  # Time stepping parameters
    47   &PARM03
    48   startTime =           0.,
    49   endTime   =     8640000.,
    50  #endTime   = 62208000000.,
    51   deltaTmom      =   1200.0,
    52   deltaTtracer   = 172800.0,
    53   deltaTfreesurf = 172800.0,
    54   deltaTClock    = 172800.0,
    55   pChkptFreq  = 3110400000.,
    56   dumpFreq    = 3110400000.,
    57   taveFreq    = 3110400000.,
    58   monitorFreq =   31104000.,
    59   periodicExternalForcing=.TRUE.,
    60   externForcingPeriod=2592000.,
    61   externForcingCycle=31104000.,
    62   tauThetaClimRelax=5184000.0,
    63   abEps=0.1,
    64   &
    65  
    66  # Gridding parameters
    67   &PARM04
    68   usingCartesianGrid=.FALSE.,
    69   usingSphericalPolarGrid=.TRUE.,
    70   dXspacing=4.,
    71   dYspacing=4.,
    72   Ro_SeaLevel=53023122.566084,
    73   groundAtK1=.TRUE.,
    74   delR=7103300.720021, 6570548.440790, 6041670.010249,
    75        5516436.666057, 4994602.034410, 4475903.435290,
    76        3960063.245801, 3446790.312651, 2935781.405664,
    77        2426722.705046, 1919291.315988, 1413156.804970,
    78        1008846.750166,  705919.025481,  504089.693499,
    79   ygOrigin=-80.,
    80   &
    81  
    82  # Input datasets
    83   &PARM05
    84   topoFile       ='topog.bin',
    85   pLoadFile      ='deltageopotjmd95.bin',
    86   hydrogThetaFile='lev_t.bin',
    87   hydrogSaltFile ='lev_s.bin',
    88   zonalWindFile  ='trenberth_taux.bin',
    89   meridWindFile  ='trenberth_tauy.bin',
    90   thetaClimFile  ='lev_sst.bin',
    91  #saltClimFile   ='lev_sss.bin',
    92   surfQFile      ='shi_qnet.bin',
    93   EmPmRFile      ='shi_empmr.bin',
    94   &
\end{verbatim}
1	mlosch	1.1	\begin{verbatim}
2			1 # ====================
3			2 # \| Model parameters \|
4			3 # ====================
5			4 #
6			5 # Continuous equation parameters
7			6 &PARM01
8			7 tRef=15*20.,
9			8 sRef=15*35.,
10			9 viscAh =3.E5,
11			10 no_slip_sides=.TRUE.,
12			11 viscAr =1.721611620915750e5,
13			12 #viscAz =1.67E-3,
14			13 no_slip_bottom=.FALSE.,
15			14 diffKhT=1.E3,
16			15 diffKrT=5.154525811125000e3,
17			16 #diffKzT=0.5E-4,
18			17 diffKhS=1.E3,
19			18 diffKrS=5.154525811125000e3,
20			19 #diffKzS=0.5E-4,
21			20 implicitDiffusion=.TRUE.,
22			21 ivdc_kappa=1.030905162225000e9,
23			22 #ivdc_kappa=10.0,
24			23 gravity=9.81,
25			24 gravitySign=-1.D0,
26			25 rhonil=1035.,
27			26 buoyancyRelation='OCEANICP',
28			27 eosType='JMD95P',
29			28 rigidLid=.FALSE.,
30			29 implicitFreeSurface=.TRUE.,
31			30 exactConserv=.TRUE.,
32			31 nonlinFreeSurf=3,
33			32 hFacInf=0.2,
34			33 hFacSup=2.0,
35			34 useRealFreshWaterFlux=.FALSE.,
36			35 readBinaryPrec=64,
37			36 writeBinaryPrec=64,
38			37 cosPower=0.5,
39			38 &
40			39
41			40 # Elliptic solver parameters
42			41 &PARM02
43			42 cg2dMaxIters=200,
44			43 cg2dTargetResidual=1.E-9,
45			44 &
46			45
47			46 # Time stepping parameters
48			47 &PARM03
49			48 startTime = 0.,
50			49 endTime = 8640000.,
51			50 #endTime = 62208000000.,
52			51 deltaTmom = 1200.0,
53			52 deltaTtracer = 172800.0,
54			53 deltaTfreesurf = 172800.0,
55			54 deltaTClock = 172800.0,
56			55 pChkptFreq = 3110400000.,
57			56 dumpFreq = 3110400000.,
58			57 taveFreq = 3110400000.,
59			58 monitorFreq = 31104000.,
60			59 periodicExternalForcing=.TRUE.,
61			60 externForcingPeriod=2592000.,
62			61 externForcingCycle=31104000.,
63			62 tauThetaClimRelax=5184000.0,
64			63 abEps=0.1,
65			64 &
66			65
67			66 # Gridding parameters
68			67 &PARM04
69			68 usingCartesianGrid=.FALSE.,
70			69 usingSphericalPolarGrid=.TRUE.,
71			70 dXspacing=4.,
72			71 dYspacing=4.,
73			72 Ro_SeaLevel=53023122.566084,
74			73 groundAtK1=.TRUE.,
75			74 delR=7103300.720021, 6570548.440790, 6041670.010249,
76			75 5516436.666057, 4994602.034410, 4475903.435290,
77			76 3960063.245801, 3446790.312651, 2935781.405664,
78			77 2426722.705046, 1919291.315988, 1413156.804970,
79			78 1008846.750166, 705919.025481, 504089.693499,
80	jmc	1.2	79 ygOrigin=-80.,
81	mlosch	1.1	80 &
82			81
83			82 # Input datasets
84			83 &PARM05
85			84 topoFile ='topog.bin',
86			85 pLoadFile ='deltageopotjmd95.bin',
87			86 hydrogThetaFile='lev_t.bin',
88			87 hydrogSaltFile ='lev_s.bin',
89			88 zonalWindFile ='trenberth_taux.bin',
90			89 meridWindFile ='trenberth_tauy.bin',
91			90 thetaClimFile ='lev_sst.bin',
92			91 #saltClimFile ='lev_sss.bin',
93			92 surfQFile ='shi_qnet.bin',
94			93 EmPmRFile ='shi_empmr.bin',
95			94 &
96			\end{verbatim}