/[MITgcm]/MITgcm/verification/global_ocean.90x40x15/results/output_oadm.txt
ViewVC logotype

Contents of /MITgcm/verification/global_ocean.90x40x15/results/output_oadm.txt

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.8 - (show annotations) (download)
Thu Oct 6 19:50:14 2016 UTC (7 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, HEAD
Changes since 1.7: +258 -253 lines
File MIME type: text/plain
- update results after adding parenthesis in main diagonal calculation
  of implicit vertical diffusion matrix (impldiff.F & gad_implicit_r.F).

1 (PID.TID 0000.0001)
2 (PID.TID 0000.0001) // ======================================================
3 (PID.TID 0000.0001) // MITgcm UV
4 (PID.TID 0000.0001) // =========
5 (PID.TID 0000.0001) // ======================================================
6 (PID.TID 0000.0001) // execution environment starting up...
7 (PID.TID 0000.0001)
8 (PID.TID 0000.0001) // MITgcmUV version: checkpoint65z
9 (PID.TID 0000.0001) // Build user: jmc
10 (PID.TID 0000.0001) // Build host: baudelaire
11 (PID.TID 0000.0001) // Build date: Wed Oct 5 21:36:20 EDT 2016
12 (PID.TID 0000.0001)
13 (PID.TID 0000.0001) // =======================================================
14 (PID.TID 0000.0001) // Execution Environment parameter file "eedata"
15 (PID.TID 0000.0001) // =======================================================
16 (PID.TID 0000.0001) ># Example "eedata" file
17 (PID.TID 0000.0001) ># Lines beginning "#" are comments
18 (PID.TID 0000.0001) ># nTx - No. threads per process in X
19 (PID.TID 0000.0001) ># nTy - No. threads per process in Y
20 (PID.TID 0000.0001) > &EEPARMS
21 (PID.TID 0000.0001) > nTx=1,
22 (PID.TID 0000.0001) > nTy=1,
23 (PID.TID 0000.0001) > /
24 (PID.TID 0000.0001) ># Note: Some systems use / as the
25 (PID.TID 0000.0001) ># namelist terminator. Other systems
26 (PID.TID 0000.0001) ># use a & character (as shown here).
27 (PID.TID 0000.0001)
28 (PID.TID 0000.0001) // =======================================================
29 (PID.TID 0000.0001) // Computational Grid Specification ( see files "SIZE.h" )
30 (PID.TID 0000.0001) // ( and "eedata" )
31 (PID.TID 0000.0001) // =======================================================
32 (PID.TID 0000.0001) nPx = 1 ; /* No. processes in X */
33 (PID.TID 0000.0001) nPy = 1 ; /* No. processes in Y */
34 (PID.TID 0000.0001) nSx = 2 ; /* No. tiles in X per process */
35 (PID.TID 0000.0001) nSy = 2 ; /* No. tiles in Y per process */
36 (PID.TID 0000.0001) sNx = 45 ; /* Tile size in X */
37 (PID.TID 0000.0001) sNy = 20 ; /* Tile size in Y */
38 (PID.TID 0000.0001) OLx = 3 ; /* Tile overlap distance in X */
39 (PID.TID 0000.0001) OLy = 3 ; /* Tile overlap distance in Y */
40 (PID.TID 0000.0001) nTx = 1 ; /* No. threads in X per process */
41 (PID.TID 0000.0001) nTy = 1 ; /* No. threads in Y per process */
42 (PID.TID 0000.0001) Nr = 15 ; /* No. levels in the vertical */
43 (PID.TID 0000.0001) Nx = 90 ; /* Total domain size in X ( = nPx*nSx*sNx ) */
44 (PID.TID 0000.0001) Ny = 40 ; /* Total domain size in Y ( = nPy*nSy*sNy ) */
45 (PID.TID 0000.0001) nTiles = 4 ; /* Total no. tiles per process ( = nSx*nSy ) */
46 (PID.TID 0000.0001) nProcs = 1 ; /* Total no. processes ( = nPx*nPy ) */
47 (PID.TID 0000.0001) nThreads = 1 ; /* Total no. threads per process ( = nTx*nTy ) */
48 (PID.TID 0000.0001) usingMPI = F ; /* Flag used to control whether MPI is in use */
49 (PID.TID 0000.0001) /* note: To execute a program with MPI calls */
50 (PID.TID 0000.0001) /* it must be launched appropriately e.g */
51 (PID.TID 0000.0001) /* "mpirun -np 64 ......" */
52 (PID.TID 0000.0001) useCoupler= F ;/* Flag used to control communications with */
53 (PID.TID 0000.0001) /* other model components, through a coupler */
54 (PID.TID 0000.0001) debugMode = F ; /* print debug msg. (sequence of S/R calls) */
55 (PID.TID 0000.0001) printMapIncludesZeros= F ; /* print zeros in Std.Output maps */
56 (PID.TID 0000.0001) maxLengthPrt1D= 65 /* maxLength of 1D array printed to StdOut */
57 (PID.TID 0000.0001)
58 (PID.TID 0000.0001) // ======================================================
59 (PID.TID 0000.0001) // Mapping of tiles to threads
60 (PID.TID 0000.0001) // ======================================================
61 (PID.TID 0000.0001) // -o- Thread 1, tiles ( 1: 2, 1: 2)
62 (PID.TID 0000.0001)
63 (PID.TID 0000.0001) // ======================================================
64 (PID.TID 0000.0001) // Tile <-> Tile connectvity table
65 (PID.TID 0000.0001) // ======================================================
66 (PID.TID 0000.0001) // Tile number: 000001 (process no. = 000000)
67 (PID.TID 0000.0001) // WEST: Tile = 000002, Process = 000000, Comm = put
68 (PID.TID 0000.0001) // bi = 000002, bj = 000001
69 (PID.TID 0000.0001) // EAST: Tile = 000002, Process = 000000, Comm = put
70 (PID.TID 0000.0001) // bi = 000002, bj = 000001
71 (PID.TID 0000.0001) // SOUTH: Tile = 000003, Process = 000000, Comm = put
72 (PID.TID 0000.0001) // bi = 000001, bj = 000002
73 (PID.TID 0000.0001) // NORTH: Tile = 000003, Process = 000000, Comm = put
74 (PID.TID 0000.0001) // bi = 000001, bj = 000002
75 (PID.TID 0000.0001) // Tile number: 000002 (process no. = 000000)
76 (PID.TID 0000.0001) // WEST: Tile = 000001, Process = 000000, Comm = put
77 (PID.TID 0000.0001) // bi = 000001, bj = 000001
78 (PID.TID 0000.0001) // EAST: Tile = 000001, Process = 000000, Comm = put
79 (PID.TID 0000.0001) // bi = 000001, bj = 000001
80 (PID.TID 0000.0001) // SOUTH: Tile = 000004, Process = 000000, Comm = put
81 (PID.TID 0000.0001) // bi = 000002, bj = 000002
82 (PID.TID 0000.0001) // NORTH: Tile = 000004, Process = 000000, Comm = put
83 (PID.TID 0000.0001) // bi = 000002, bj = 000002
84 (PID.TID 0000.0001) // Tile number: 000003 (process no. = 000000)
85 (PID.TID 0000.0001) // WEST: Tile = 000004, Process = 000000, Comm = put
86 (PID.TID 0000.0001) // bi = 000002, bj = 000002
87 (PID.TID 0000.0001) // EAST: Tile = 000004, Process = 000000, Comm = put
88 (PID.TID 0000.0001) // bi = 000002, bj = 000002
89 (PID.TID 0000.0001) // SOUTH: Tile = 000001, Process = 000000, Comm = put
90 (PID.TID 0000.0001) // bi = 000001, bj = 000001
91 (PID.TID 0000.0001) // NORTH: Tile = 000001, Process = 000000, Comm = put
92 (PID.TID 0000.0001) // bi = 000001, bj = 000001
93 (PID.TID 0000.0001) // Tile number: 000004 (process no. = 000000)
94 (PID.TID 0000.0001) // WEST: Tile = 000003, Process = 000000, Comm = put
95 (PID.TID 0000.0001) // bi = 000001, bj = 000002
96 (PID.TID 0000.0001) // EAST: Tile = 000003, Process = 000000, Comm = put
97 (PID.TID 0000.0001) // bi = 000001, bj = 000002
98 (PID.TID 0000.0001) // SOUTH: Tile = 000002, Process = 000000, Comm = put
99 (PID.TID 0000.0001) // bi = 000002, bj = 000001
100 (PID.TID 0000.0001) // NORTH: Tile = 000002, Process = 000000, Comm = put
101 (PID.TID 0000.0001) // bi = 000002, bj = 000001
102 (PID.TID 0000.0001)
103 (PID.TID 0000.0001) INI_PARMS: opening model parameter file "data"
104 (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data
105 (PID.TID 0000.0001) // =======================================================
106 (PID.TID 0000.0001) // Parameter file "data"
107 (PID.TID 0000.0001) // =======================================================
108 (PID.TID 0000.0001) ># ====================
109 (PID.TID 0000.0001) ># | Model parameters |
110 (PID.TID 0000.0001) ># ====================
111 (PID.TID 0000.0001) >#
112 (PID.TID 0000.0001) ># Continuous equation parameters
113 (PID.TID 0000.0001) > &PARM01
114 (PID.TID 0000.0001) > tRef = 15*20.,
115 (PID.TID 0000.0001) > sRef = 15*35.,
116 (PID.TID 0000.0001) > viscAr=1.E-3,
117 (PID.TID 0000.0001) > diffKhT=0.,
118 (PID.TID 0000.0001) >#- diffKrT unused when compiled with ALLOW_3D_DIFFKR
119 (PID.TID 0000.0001) >#diffKrT=3.E-5,
120 (PID.TID 0000.0001) > diffKhS=0.,
121 (PID.TID 0000.0001) > diffKrS=3.E-5,
122 (PID.TID 0000.0001) >#ph(
123 (PID.TID 0000.0001) >#ph Many of the following params. are just for AD testing;
124 (PID.TID 0000.0001) >#ph not necesarily meaningful in a 4x4 deg. setting.
125 (PID.TID 0000.0001) > viscC4Leith=1.5,
126 (PID.TID 0000.0001) > viscC4Leithd=1.5,
127 (PID.TID 0000.0001) > viscA4GridMax=0.5,
128 (PID.TID 0000.0001) > useAreaViscLength=.TRUE.,
129 (PID.TID 0000.0001) >### viscAh=5.E5,
130 (PID.TID 0000.0001) >#ph)
131 (PID.TID 0000.0001) > bottomDragQuadratic = 0.0021,
132 (PID.TID 0000.0001) > rhoConst=1035.,
133 (PID.TID 0000.0001) > gravity=9.81,
134 (PID.TID 0000.0001) > eosType = 'JMD95Z',
135 (PID.TID 0000.0001) > ivdc_kappa=10.,
136 (PID.TID 0000.0001) > implicitDiffusion=.TRUE.,
137 (PID.TID 0000.0001) > allowFreezing=.TRUE.,
138 (PID.TID 0000.0001) > useRealFreshWaterFlux=.TRUE.,
139 (PID.TID 0000.0001) > useCDscheme=.TRUE.,
140 (PID.TID 0000.0001) > exactConserv=.TRUE.,
141 (PID.TID 0000.0001) ># turn on looped cells
142 (PID.TID 0000.0001) > hFacMin=.05,
143 (PID.TID 0000.0001) > hFacMinDr=50.,
144 (PID.TID 0000.0001) >#ph(
145 (PID.TID 0000.0001) > vectorInvariantMomentum=.TRUE.,
146 (PID.TID 0000.0001) > hFacInf=0.2,
147 (PID.TID 0000.0001) > hFacSup=2.0,
148 (PID.TID 0000.0001) > nonlinFreeSurf=2,
149 (PID.TID 0000.0001) > staggerTimeStep=.TRUE.,
150 (PID.TID 0000.0001) >### multiDimAdvection=.FALSE.,
151 (PID.TID 0000.0001) > tempAdvScheme=30,
152 (PID.TID 0000.0001) > saltAdvScheme=30,
153 (PID.TID 0000.0001) > tempVertAdvScheme=20,
154 (PID.TID 0000.0001) > saltVertAdvScheme=20,
155 (PID.TID 0000.0001) > tempImplVertAdv=.TRUE.,
156 (PID.TID 0000.0001) > saltImplVertAdv=.TRUE.,
157 (PID.TID 0000.0001) >#ph)
158 (PID.TID 0000.0001) ># set precision of data files
159 (PID.TID 0000.0001) > readBinaryPrec=32,
160 (PID.TID 0000.0001) >#debugLevel=3,
161 (PID.TID 0000.0001) > /
162 (PID.TID 0000.0001) >
163 (PID.TID 0000.0001) ># Elliptic solver parameters
164 (PID.TID 0000.0001) > &PARM02
165 (PID.TID 0000.0001) > cg2dMaxIters=500,
166 (PID.TID 0000.0001) > cg2dTargetResidual=1.E-13,
167 (PID.TID 0000.0001) > /
168 (PID.TID 0000.0001) >
169 (PID.TID 0000.0001) ># Time stepping parameters
170 (PID.TID 0000.0001) > &PARM03
171 (PID.TID 0000.0001) > nIter0 = 0,
172 (PID.TID 0000.0001) > nTimeSteps = 10,
173 (PID.TID 0000.0001) ># 100 years of integration will yield a reasonable flow field
174 (PID.TID 0000.0001) ># startTime = 0.,
175 (PID.TID 0000.0001) ># endTime = 3110400000.,
176 (PID.TID 0000.0001) > deltaTmom = 1200.0,
177 (PID.TID 0000.0001) > tauCD = 321428.,
178 (PID.TID 0000.0001) > deltaTtracer= 43200.0,
179 (PID.TID 0000.0001) > deltaTClock = 43200.0,
180 (PID.TID 0000.0001) >#deltaTfreesurf= 43200.,
181 (PID.TID 0000.0001) > abEps = 0.1,
182 (PID.TID 0000.0001) > pChkptFreq= 311040000.,
183 (PID.TID 0000.0001) > dumpFreq= 311040000.,
184 (PID.TID 0000.0001) > dumpFreq= 864000.,
185 (PID.TID 0000.0001) > adjDumpFreq= 432000.,
186 (PID.TID 0000.0001) > adjMonitorFreq=0.,
187 (PID.TID 0000.0001) ># taveFreq= 311040000.,
188 (PID.TID 0000.0001) ># taveFreq= 864000.,
189 (PID.TID 0000.0001) >#monitorFreq=31104000.,
190 (PID.TID 0000.0001) > monitorFreq=0.,
191 (PID.TID 0000.0001) ># 2 months restoring timescale for temperature
192 (PID.TID 0000.0001) > tauThetaClimRelax = 5184000.0,
193 (PID.TID 0000.0001) ># 6 months restoring timescale for salinity
194 (PID.TID 0000.0001) > tauSaltClimRelax = 15552000.0,
195 (PID.TID 0000.0001) > periodicExternalForcing=.TRUE.,
196 (PID.TID 0000.0001) > externForcingPeriod=2592000.,
197 (PID.TID 0000.0001) > externForcingCycle=31104000.,
198 (PID.TID 0000.0001) > /
199 (PID.TID 0000.0001) >
200 (PID.TID 0000.0001) ># Gridding parameters
201 (PID.TID 0000.0001) > &PARM04
202 (PID.TID 0000.0001) > usingSphericalPolarGrid=.TRUE.,
203 (PID.TID 0000.0001) > delR= 50., 70., 100., 140., 190.,
204 (PID.TID 0000.0001) > 240., 290., 340., 390., 440.,
205 (PID.TID 0000.0001) > 490., 540., 590., 640., 690.,
206 (PID.TID 0000.0001) > ygOrigin=-80.,
207 (PID.TID 0000.0001) > dySpacing=4.,
208 (PID.TID 0000.0001) > dxSpacing=4.,
209 (PID.TID 0000.0001) > /
210 (PID.TID 0000.0001) >
211 (PID.TID 0000.0001) ># Input datasets
212 (PID.TID 0000.0001) > &PARM05
213 (PID.TID 0000.0001) > bathyFile= 'bathymetry.bin',
214 (PID.TID 0000.0001) > hydrogThetaFile='lev_t.bin',
215 (PID.TID 0000.0001) > hydrogSaltFile= 'lev_s.bin',
216 (PID.TID 0000.0001) > zonalWindFile= 'trenberth_taux.bin',
217 (PID.TID 0000.0001) > meridWindFile= 'trenberth_tauy.bin',
218 (PID.TID 0000.0001) > thetaClimFile= 'lev_sst.bin',
219 (PID.TID 0000.0001) > saltClimFile= 'lev_sss.bin',
220 (PID.TID 0000.0001) > surfQFile= 'ncep_qnet.bin',
221 (PID.TID 0000.0001) ># fresh water flux is turned off, uncomment next line to turn on
222 (PID.TID 0000.0001) ># (not recommened together with surface salinity restoring)
223 (PID.TID 0000.0001) ># EmPmRFile= 'ncep_emp.bin',
224 (PID.TID 0000.0001) > /
225 (PID.TID 0000.0001)
226 (PID.TID 0000.0001) INI_PARMS ; starts to read PARM01
227 (PID.TID 0000.0001) INI_PARMS ; read PARM01 : OK
228 (PID.TID 0000.0001) S/R INI_PARMS: No request for barotropic solver
229 (PID.TID 0000.0001) S/R INI_PARMS: => Use implicitFreeSurface as default
230 (PID.TID 0000.0001) INI_PARMS ; starts to read PARM02
231 (PID.TID 0000.0001) INI_PARMS ; read PARM02 : OK
232 (PID.TID 0000.0001) INI_PARMS ; starts to read PARM03
233 (PID.TID 0000.0001) INI_PARMS ; read PARM03 : OK
234 (PID.TID 0000.0001) INI_PARMS ; starts to read PARM04
235 (PID.TID 0000.0001) INI_PARMS ; read PARM04 : OK
236 (PID.TID 0000.0001) INI_PARMS ; starts to read PARM05
237 (PID.TID 0000.0001) INI_PARMS ; read PARM05 : OK
238 (PID.TID 0000.0001) INI_PARMS: finished reading file "data"
239 (PID.TID 0000.0001) PACKAGES_BOOT: opening data.pkg
240 (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.pkg
241 (PID.TID 0000.0001) // =======================================================
242 (PID.TID 0000.0001) // Parameter file "data.pkg"
243 (PID.TID 0000.0001) // =======================================================
244 (PID.TID 0000.0001) >#
245 (PID.TID 0000.0001) ># ********
246 (PID.TID 0000.0001) ># Packages
247 (PID.TID 0000.0001) ># ********
248 (PID.TID 0000.0001) > &PACKAGES
249 (PID.TID 0000.0001) > useGMRedi = .TRUE.,
250 (PID.TID 0000.0001) > useKPP = .FALSE.,
251 (PID.TID 0000.0001) > useGrdchk = .TRUE.,
252 (PID.TID 0000.0001) ># useMNC = .TRUE.,
253 (PID.TID 0000.0001) > /
254 (PID.TID 0000.0001)
255 (PID.TID 0000.0001) PACKAGES_BOOT: finished reading data.pkg
256 (PID.TID 0000.0001) PACKAGES_BOOT: On/Off package Summary
257 -------- pkgs with a standard "usePKG" On/Off switch in "data.pkg": --------
258 pkg/gmredi compiled and used ( useGMRedi = T )
259 pkg/autodiff compiled and used ( useAUTODIFF = T )
260 pkg/grdchk compiled and used ( useGrdchk = T )
261 pkg/ctrl compiled and used ( useCTRL = T )
262 -------- pkgs without standard "usePKG" On/Off switch in "data.pkg": --------
263 pkg/generic_advdiff compiled and used ( useGAD = T )
264 pkg/mom_common compiled and used ( momStepping = T )
265 pkg/mom_vecinv compiled and used ( +vectorInvariantMomentum = T )
266 pkg/mom_fluxform compiled but not used ( & not vectorInvariantMom = F )
267 pkg/cd_code compiled and used ( useCDscheme = T )
268 pkg/debug compiled but not used ( debugMode = F )
269 pkg/rw compiled and used
270 pkg/mdsio compiled and used
271 pkg/autodiff compiled and used
272 pkg/openad compiled and used
273 pkg/cost compiled and used
274 (PID.TID 0000.0001) PACKAGES_BOOT: End of package Summary
275 (PID.TID 0000.0001)
276 (PID.TID 0000.0001) GM_READPARMS: opening data.gmredi
277 (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.gmredi
278 (PID.TID 0000.0001) // =======================================================
279 (PID.TID 0000.0001) // Parameter file "data.gmredi"
280 (PID.TID 0000.0001) // =======================================================
281 (PID.TID 0000.0001) ># GM+Redi package parameters:
282 (PID.TID 0000.0001) ># GM_Small_Number :: epsilon used in computing the slope
283 (PID.TID 0000.0001) ># GM_slopeSqCutoff :: slope^2 cut-off value
284 (PID.TID 0000.0001) >
285 (PID.TID 0000.0001) >#-from MOM :
286 (PID.TID 0000.0001) ># GM_background_K: G & Mc.W diffusion coefficient
287 (PID.TID 0000.0001) ># GM_maxSlope : max slope of isopycnals
288 (PID.TID 0000.0001) ># GM_Scrit : transition for scaling diffusion coefficient
289 (PID.TID 0000.0001) ># GM_Sd : half width scaling for diffusion coefficient
290 (PID.TID 0000.0001) ># GM_taper_scheme: slope clipping or one of the tapering schemes
291 (PID.TID 0000.0001) ># GM_Kmin_horiz : horizontal diffusion minimum value
292 (PID.TID 0000.0001) >
293 (PID.TID 0000.0001) >#-Option parameters (needs to "define" options in GMREDI_OPTIONS.h")
294 (PID.TID 0000.0001) ># GM_isopycK : isopycnal diffusion coefficient (default=GM_background_K)
295 (PID.TID 0000.0001) ># GM_AdvForm : turn on GM Advective form (default=Skew flux form)
296 (PID.TID 0000.0001) >
297 (PID.TID 0000.0001) > &GM_PARM01
298 (PID.TID 0000.0001) > GM_Small_Number = 1.D-20,
299 (PID.TID 0000.0001) > GM_slopeSqCutoff = 1.D+08,
300 (PID.TID 0000.0001) > GM_AdvForm = .TRUE.,
301 (PID.TID 0000.0001) > GM_isopycK = 1.1D+3,
302 (PID.TID 0000.0001) > GM_background_K = 0.9D+3,
303 (PID.TID 0000.0001) > GM_taper_scheme = 'dm95',
304 (PID.TID 0000.0001) > GM_maxSlope = 1.D-2,
305 (PID.TID 0000.0001) > GM_Kmin_horiz = 50.,
306 (PID.TID 0000.0001) > GM_Scrit = 4.D-3,
307 (PID.TID 0000.0001) > GM_Sd = 1.D-3,
308 (PID.TID 0000.0001) ># GM_Visbeck_alpha = 1.5D-2,
309 (PID.TID 0000.0001) > GM_Visbeck_alpha = 0.D0,
310 (PID.TID 0000.0001) > GM_Visbeck_length = 2.D+5,
311 (PID.TID 0000.0001) > GM_Visbeck_depth = 1.D+3,
312 (PID.TID 0000.0001) > GM_Visbeck_maxval_K= 2.5D+3,
313 (PID.TID 0000.0001) > /
314 (PID.TID 0000.0001) >
315 (PID.TID 0000.0001) >
316 (PID.TID 0000.0001)
317 (PID.TID 0000.0001) GM_READPARMS: finished reading data.gmredi
318 (PID.TID 0000.0001) AUTODIFF_READPARMS: opening data.autodiff
319 (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.autodiff
320 (PID.TID 0000.0001) // =======================================================
321 (PID.TID 0000.0001) // Parameter file "data.autodiff"
322 (PID.TID 0000.0001) // =======================================================
323 (PID.TID 0000.0001) ># =========================
324 (PID.TID 0000.0001) ># pkg AUTODIFF parameters :
325 (PID.TID 0000.0001) ># =========================
326 (PID.TID 0000.0001) ># inAdExact :: get an exact adjoint (no approximation) (def=.True.)
327 (PID.TID 0000.0001) >#
328 (PID.TID 0000.0001) > &AUTODIFF_PARM01
329 (PID.TID 0000.0001) ># inAdExact = .FALSE.,
330 (PID.TID 0000.0001) > /
331 (PID.TID 0000.0001)
332 (PID.TID 0000.0001) AUTODIFF_READPARMS: finished reading data.autodiff
333 (PID.TID 0000.0001) // ===================================
334 (PID.TID 0000.0001) // AUTODIFF parameters :
335 (PID.TID 0000.0001) // ===================================
336 (PID.TID 0000.0001) inAdExact = /* get an exact adjoint (no approximation) */
337 (PID.TID 0000.0001) T
338 (PID.TID 0000.0001) ;
339 (PID.TID 0000.0001) useKPPinAdMode = /* use KPP in adjoint mode */
340 (PID.TID 0000.0001) F
341 (PID.TID 0000.0001) ;
342 (PID.TID 0000.0001) useGMRediInAdMode = /* use GMRedi in adjoint mode */
343 (PID.TID 0000.0001) T
344 (PID.TID 0000.0001) ;
345 (PID.TID 0000.0001) useSEAICEinAdMode = /* use SEAICE in adjoint mode */
346 (PID.TID 0000.0001) F
347 (PID.TID 0000.0001) ;
348 (PID.TID 0000.0001) useGGL90inAdMode = /* use GGL90 in adjoint mode */
349 (PID.TID 0000.0001) F
350 (PID.TID 0000.0001) ;
351 (PID.TID 0000.0001) useSALT_PLUMEinAdMode = /* use SALT_PLUME in adjoint mode */
352 (PID.TID 0000.0001) F
353 (PID.TID 0000.0001) ;
354 (PID.TID 0000.0001) dumpAdVarExch = /* control adexch before dumpinp */
355 (PID.TID 0000.0001) 2
356 (PID.TID 0000.0001) ;
357 (PID.TID 0000.0001) mon_AdVarExch = /* control adexch before monitor */
358 (PID.TID 0000.0001) 2
359 (PID.TID 0000.0001) ;
360 (PID.TID 0000.0001) viscFacInAd = /* viscosity factor for adjoint */
361 (PID.TID 0000.0001) 1.000000000000000E+00
362 (PID.TID 0000.0001) ;
363 (PID.TID 0000.0001)
364 (PID.TID 0000.0001) OPTIM_READPARMS: opening data.optim
365 (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.optim
366 (PID.TID 0000.0001) // =======================================================
367 (PID.TID 0000.0001) // Parameter file "data.optim"
368 (PID.TID 0000.0001) // =======================================================
369 (PID.TID 0000.0001) >#
370 (PID.TID 0000.0001) ># ********************************
371 (PID.TID 0000.0001) ># Off-line optimization parameters
372 (PID.TID 0000.0001) ># ********************************
373 (PID.TID 0000.0001) > &OPTIM
374 (PID.TID 0000.0001) > optimcycle=0,
375 (PID.TID 0000.0001) > /
376 (PID.TID 0000.0001)
377 (PID.TID 0000.0001) OPTIM_READPARMS: finished reading data.optim
378 (PID.TID 0000.0001) CTRL_READPARMS: opening data.ctrl
379 (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.ctrl
380 (PID.TID 0000.0001) // =======================================================
381 (PID.TID 0000.0001) // Parameter file "data.ctrl"
382 (PID.TID 0000.0001) // =======================================================
383 (PID.TID 0000.0001) >#
384 (PID.TID 0000.0001) >#
385 (PID.TID 0000.0001) ># *********************
386 (PID.TID 0000.0001) ># ECCO controlvariables
387 (PID.TID 0000.0001) ># *********************
388 (PID.TID 0000.0001) > &CTRL_NML
389 (PID.TID 0000.0001) > ctrlUseGen=.FALSE.,
390 (PID.TID 0000.0001) ># doSinglePrecTapelev= .FALSE.,
391 (PID.TID 0000.0001) > xx_theta_file = 'xx_theta',
392 (PID.TID 0000.0001) > xx_salt_file = 'xx_salt',
393 (PID.TID 0000.0001) > xx_tr1_file = 'xx_tr1',
394 (PID.TID 0000.0001) > xx_hflux_file = 'xx_hflux',
395 (PID.TID 0000.0001) > xx_sflux_file = 'xx_sflux',
396 (PID.TID 0000.0001) > xx_tauu_file = 'xx_tauu',
397 (PID.TID 0000.0001) > xx_tauv_file = 'xx_tauv',
398 (PID.TID 0000.0001) > xx_diffkr_file = 'xx_diffkr',
399 (PID.TID 0000.0001) > xx_kapgm_file = 'xx_kapgm',
400 (PID.TID 0000.0001) > /
401 (PID.TID 0000.0001) >#
402 (PID.TID 0000.0001) ># *********************
403 (PID.TID 0000.0001) ># names for ctrl_pack/unpack
404 (PID.TID 0000.0001) ># *********************
405 (PID.TID 0000.0001) > &CTRL_PACKNAMES
406 (PID.TID 0000.0001) > /
407 (PID.TID 0000.0001) >
408 (PID.TID 0000.0001)
409 (PID.TID 0000.0001) CTRL_READPARMS: finished reading data.ctrl
410 (PID.TID 0000.0001) COST_READPARMS: opening data.cost
411 (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.cost
412 (PID.TID 0000.0001) // =======================================================
413 (PID.TID 0000.0001) // Parameter file "data.cost"
414 (PID.TID 0000.0001) // =======================================================
415 (PID.TID 0000.0001) >#
416 (PID.TID 0000.0001) >#
417 (PID.TID 0000.0001) ># ******************
418 (PID.TID 0000.0001) ># ECCO cost function
419 (PID.TID 0000.0001) ># ******************
420 (PID.TID 0000.0001) > &COST_NML
421 (PID.TID 0000.0001) >#
422 (PID.TID 0000.0001) > mult_tracer = 1.,
423 (PID.TID 0000.0001) > mult_test = 1.,
424 (PID.TID 0000.0001) > mult_atl = 1.,
425 (PID.TID 0000.0001) > /
426 (PID.TID 0000.0001)
427 (PID.TID 0000.0001) COST_READPARMS: finished reading data.cost
428 (PID.TID 0000.0001) GRDCHK_READPARMS: opening data.grdchk
429 (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.grdchk
430 (PID.TID 0000.0001) // =======================================================
431 (PID.TID 0000.0001) // Parameter file "data.grdchk"
432 (PID.TID 0000.0001) // =======================================================
433 (PID.TID 0000.0001) >
434 (PID.TID 0000.0001) ># *******************
435 (PID.TID 0000.0001) ># ECCO gradient check
436 (PID.TID 0000.0001) ># *******************
437 (PID.TID 0000.0001) > &GRDCHK_NML
438 (PID.TID 0000.0001) > grdchk_eps = 1.d-2,
439 (PID.TID 0000.0001) >### iGloPos = 76,
440 (PID.TID 0000.0001) >### jGloPos = 27,
441 (PID.TID 0000.0001) > iGloPos = 31,
442 (PID.TID 0000.0001) > jGloPos = 7,
443 (PID.TID 0000.0001) > kGloPos = 1,
444 (PID.TID 0000.0001) > iGloTile = 2,
445 (PID.TID 0000.0001) > jgloTile = 2,
446 (PID.TID 0000.0001) >### nbeg = 1,
447 (PID.TID 0000.0001) > nstep = 1,
448 (PID.TID 0000.0001) > nend = 3,
449 (PID.TID 0000.0001) > grdchkvarindex = 1,
450 (PID.TID 0000.0001) > /
451 (PID.TID 0000.0001)
452 (PID.TID 0000.0001) GRDCHK_READPARMS: finished reading data.grdchk
453 (PID.TID 0000.0001)
454 (PID.TID 0000.0001) // =======================================================
455 (PID.TID 0000.0001) // Gradient check configuration >>> START <<<
456 (PID.TID 0000.0001) // =======================================================
457 (PID.TID 0000.0001)
458 (PID.TID 0000.0001) grdchkvarindex : 1
459 (PID.TID 0000.0001) eps: 0.100E-01
460 (PID.TID 0000.0001) First location: 0
461 (PID.TID 0000.0001) Last location: 3
462 (PID.TID 0000.0001) Increment: 1
463 (PID.TID 0000.0001) grdchkWhichProc: 0
464 (PID.TID 0000.0001) iLocTile = 2 , jLocTile = 2
465 (PID.TID 0000.0001)
466 (PID.TID 0000.0001) // =======================================================
467 (PID.TID 0000.0001) // Gradient check configuration >>> END <<<
468 (PID.TID 0000.0001) // =======================================================
469 (PID.TID 0000.0001)
470 (PID.TID 0000.0001) SET_PARMS: done
471 (PID.TID 0000.0001) Enter INI_VERTICAL_GRID: setInterFDr= T ; setCenterDr= F
472 (PID.TID 0000.0001) GAD_INIT_FIXED: GAD_OlMinSize= 2 0 1
473 (PID.TID 0000.0001)
474 (PID.TID 0000.0001) // ===================================
475 (PID.TID 0000.0001) // GAD parameters :
476 (PID.TID 0000.0001) // ===================================
477 (PID.TID 0000.0001) tempAdvScheme = /* Temp. Horiz.Advection scheme selector */
478 (PID.TID 0000.0001) 30
479 (PID.TID 0000.0001) ;
480 (PID.TID 0000.0001) tempVertAdvScheme = /* Temp. Vert. Advection scheme selector */
481 (PID.TID 0000.0001) 20
482 (PID.TID 0000.0001) ;
483 (PID.TID 0000.0001) tempMultiDimAdvec = /* use Muti-Dim Advec method for Temp */
484 (PID.TID 0000.0001) T
485 (PID.TID 0000.0001) ;
486 (PID.TID 0000.0001) tempSOM_Advection = /* use 2nd Order Moment Advection for Temp */
487 (PID.TID 0000.0001) F
488 (PID.TID 0000.0001) ;
489 (PID.TID 0000.0001) AdamsBashforthGt = /* apply Adams-Bashforth extrapolation on Gt */
490 (PID.TID 0000.0001) F
491 (PID.TID 0000.0001) ;
492 (PID.TID 0000.0001) AdamsBashforth_T = /* apply Adams-Bashforth extrapolation on Temp */
493 (PID.TID 0000.0001) F
494 (PID.TID 0000.0001) ;
495 (PID.TID 0000.0001) saltAdvScheme = /* Salt. Horiz.advection scheme selector */
496 (PID.TID 0000.0001) 30
497 (PID.TID 0000.0001) ;
498 (PID.TID 0000.0001) saltVertAdvScheme = /* Salt. Vert. Advection scheme selector */
499 (PID.TID 0000.0001) 20
500 (PID.TID 0000.0001) ;
501 (PID.TID 0000.0001) saltMultiDimAdvec = /* use Muti-Dim Advec method for Salt */
502 (PID.TID 0000.0001) T
503 (PID.TID 0000.0001) ;
504 (PID.TID 0000.0001) saltSOM_Advection = /* use 2nd Order Moment Advection for Salt */
505 (PID.TID 0000.0001) F
506 (PID.TID 0000.0001) ;
507 (PID.TID 0000.0001) AdamsBashforthGs = /* apply Adams-Bashforth extrapolation on Gs */
508 (PID.TID 0000.0001) F
509 (PID.TID 0000.0001) ;
510 (PID.TID 0000.0001) AdamsBashforth_S = /* apply Adams-Bashforth extrapolation on Salt */
511 (PID.TID 0000.0001) F
512 (PID.TID 0000.0001) ;
513 (PID.TID 0000.0001) // ===================================
514 (PID.TID 0000.0001) ctrl-wet 1: nvarlength = 96912
515 (PID.TID 0000.0001) ctrl-wet 2: surface wet C = 669
516 (PID.TID 0000.0001) ctrl-wet 3: surface wet W = 640
517 (PID.TID 0000.0001) ctrl-wet 4: surface wet S = 610
518 (PID.TID 0000.0001) ctrl-wet 4a:surface wet V = 0
519 (PID.TID 0000.0001) ctrl-wet 5: 3D wet points = 8651
520 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 1 1
521 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 2 1
522 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 3 1
523 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 4 1
524 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 5 1
525 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 6 1
526 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 7 0
527 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 8 0
528 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 9 0
529 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 10 0
530 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 11 0
531 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 12 0
532 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 13 0
533 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 14 0
534 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 15 1
535 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 16 0
536 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 17 0
537 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 18 0
538 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 19 0
539 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 20 0
540 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 21 0
541 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 22 0
542 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 23 0
543 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 24 0
544 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 25 0
545 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 26 0
546 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 27 0
547 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 28 0
548 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 29 0
549 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 30 0
550 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 31 0
551 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 32 0
552 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 33 0
553 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 34 0
554 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 35 0
555 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 36 0
556 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 37 0
557 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 38 0
558 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 39 0
559 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 40 0
560 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 41 0
561 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 42 0
562 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 43 0
563 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 44 0
564 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 45 0
565 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 46 0
566 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 47 0
567 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 48 0
568 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 49 0
569 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 50 0
570 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 51 0
571 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 52 0
572 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 53 0
573 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 54 0
574 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 55 0
575 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 56 0
576 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 57 0
577 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 58 0
578 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 59 0
579 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 60 0
580 (PID.TID 0000.0001) ctrl-wet 7: flux 19890
581 (PID.TID 0000.0001) ctrl-wet 8: atmos 17302
582 (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
583 (PID.TID 0000.0001) ctrl-wet 13: global nvarlength for Nr = 15 96912
584 (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
585 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 1 2315 2149 2206 0
586 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 2 2315 2149 2206 0
587 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 3 2254 2102 2146 0
588 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 4 2215 2058 2104 0
589 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 5 2178 2027 2070 0
590 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 6 2142 1987 2029 0
591 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 7 2114 1959 2004 0
592 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 8 2076 1918 1959 0
593 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 9 2048 1887 1925 0
594 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 10 1999 1831 1869 0
595 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 11 1948 1771 1808 0
596 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 12 1850 1653 1705 0
597 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 13 1655 1404 1458 0
598 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 14 1372 1118 1164 0
599 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 15 828 623 671 0
600 (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
601 (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
602 (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
603 (PID.TID 0000.0001) ctrl_init: no. of control variables: 7
604 (PID.TID 0000.0001) ctrl_init: control vector length: 96912
605 (PID.TID 0000.0001)
606 (PID.TID 0000.0001) // =======================================================
607 (PID.TID 0000.0001) // control vector configuration >>> START <<<
608 (PID.TID 0000.0001) // =======================================================
609 (PID.TID 0000.0001)
610 (PID.TID 0000.0001) Total number of ocean points per tile:
611 (PID.TID 0000.0001) --------------------------------------
612 (PID.TID 0000.0001) snx*sny*nr = 13500
613 (PID.TID 0000.0001)
614 (PID.TID 0000.0001) Number of ocean points per tile:
615 (PID.TID 0000.0001) --------------------------------
616 (PID.TID 0000.0001) bi,bj,#(c/s/w): 0001 0001 008651 007584 008016
617 (PID.TID 0000.0001) bi,bj,#(c/s/w): 0002 0001 010278 009360 009834
618 (PID.TID 0000.0001) bi,bj,#(c/s/w): 0001 0002 003504 003232 003049
619 (PID.TID 0000.0001) bi,bj,#(c/s/w): 0002 0002 006876 006460 006425
620 (PID.TID 0000.0001)
621 (PID.TID 0000.0001) Initial state temperature contribution:
622 (PID.TID 0000.0001) Control variable index: 0101
623 (PID.TID 0000.0001)
624 (PID.TID 0000.0001) Initial state salinity contribution:
625 (PID.TID 0000.0001) Control variable index: 0102
626 (PID.TID 0000.0001)
627 (PID.TID 0000.0001) // =======================================================
628 (PID.TID 0000.0001) // control vector configuration >>> END <<<
629 (PID.TID 0000.0001) // =======================================================
630 (PID.TID 0000.0001)
631 (PID.TID 0000.0001) INI_CG2D: CG2D normalisation factor = 6.5682677425711703E-05
632 (PID.TID 0000.0001)
633 (PID.TID 0000.0001) // =======================================================
634 (PID.TID 0000.0001) // Model configuration
635 (PID.TID 0000.0001) // =======================================================
636 (PID.TID 0000.0001) //
637 (PID.TID 0000.0001) // "Physical" paramters ( PARM01 in namelist )
638 (PID.TID 0000.0001) //
639 (PID.TID 0000.0001) buoyancyRelation = /* Type of relation to get Buoyancy */
640 (PID.TID 0000.0001) 'OCEANIC'
641 (PID.TID 0000.0001) ;
642 (PID.TID 0000.0001) fluidIsAir = /* fluid major constituent is Air */
643 (PID.TID 0000.0001) F
644 (PID.TID 0000.0001) ;
645 (PID.TID 0000.0001) fluidIsWater = /* fluid major constituent is Water */
646 (PID.TID 0000.0001) T
647 (PID.TID 0000.0001) ;
648 (PID.TID 0000.0001) usingPCoords = /* use p (or p*) vertical coordinate */
649 (PID.TID 0000.0001) F
650 (PID.TID 0000.0001) ;
651 (PID.TID 0000.0001) usingZCoords = /* use z (or z*) vertical coordinate */
652 (PID.TID 0000.0001) T
653 (PID.TID 0000.0001) ;
654 (PID.TID 0000.0001) tRef = /* Reference temperature profile ( oC or K ) */
655 (PID.TID 0000.0001) 15 @ 2.000000000000000E+01 /* K = 1: 15 */
656 (PID.TID 0000.0001) ;
657 (PID.TID 0000.0001) sRef = /* Reference salinity profile ( psu ) */
658 (PID.TID 0000.0001) 15 @ 3.500000000000000E+01 /* K = 1: 15 */
659 (PID.TID 0000.0001) ;
660 (PID.TID 0000.0001) useStrainTensionVisc= /* Use StrainTension Form of Viscous Operator */
661 (PID.TID 0000.0001) F
662 (PID.TID 0000.0001) ;
663 (PID.TID 0000.0001) useVariableVisc = /* Use variable horizontal viscosity */
664 (PID.TID 0000.0001) T
665 (PID.TID 0000.0001) ;
666 (PID.TID 0000.0001) useHarmonicVisc = /* Use harmonic horizontal viscosity */
667 (PID.TID 0000.0001) F
668 (PID.TID 0000.0001) ;
669 (PID.TID 0000.0001) useBiharmonicVisc= /* Use biharmonic horiz. viscosity */
670 (PID.TID 0000.0001) T
671 (PID.TID 0000.0001) ;
672 (PID.TID 0000.0001) useSmag3D = /* Use isotropic 3-D Smagorinsky viscosity */
673 (PID.TID 0000.0001) F
674 (PID.TID 0000.0001) ;
675 (PID.TID 0000.0001) viscAh = /* Lateral harmonic viscosity ( m^2/s ) */
676 (PID.TID 0000.0001) 0.000000000000000E+00
677 (PID.TID 0000.0001) ;
678 (PID.TID 0000.0001) viscAhMax = /* Maximum lateral harmonic viscosity ( m^2/s ) */
679 (PID.TID 0000.0001) 1.000000000000000E+21
680 (PID.TID 0000.0001) ;
681 (PID.TID 0000.0001) viscAhGrid = /* Grid dependent lateral harmonic viscosity ( non-dim. ) */
682 (PID.TID 0000.0001) 0.000000000000000E+00
683 (PID.TID 0000.0001) ;
684 (PID.TID 0000.0001) useFullLeith = /* Use Full Form of Leith Viscosity on/off flag*/
685 (PID.TID 0000.0001) F
686 (PID.TID 0000.0001) ;
687 (PID.TID 0000.0001) useAreaViscLength = /* Use area for visc length instead of geom. mean*/
688 (PID.TID 0000.0001) T
689 (PID.TID 0000.0001) ;
690 (PID.TID 0000.0001) viscC2leith = /* Leith harmonic visc. factor (on grad(vort),non-dim.) */
691 (PID.TID 0000.0001) 0.000000000000000E+00
692 (PID.TID 0000.0001) ;
693 (PID.TID 0000.0001) viscC2leithD = /* Leith harmonic viscosity factor (on grad(div),non-dim.)*/
694 (PID.TID 0000.0001) 0.000000000000000E+00
695 (PID.TID 0000.0001) ;
696 (PID.TID 0000.0001) viscC2smag = /* Smagorinsky harmonic viscosity factor (non-dim.) */
697 (PID.TID 0000.0001) 0.000000000000000E+00
698 (PID.TID 0000.0001) ;
699 (PID.TID 0000.0001) viscA4 = /* Lateral biharmonic viscosity ( m^4/s ) */
700 (PID.TID 0000.0001) 0.000000000000000E+00
701 (PID.TID 0000.0001) ;
702 (PID.TID 0000.0001) viscA4Max = /* Maximum biharmonic viscosity ( m^2/s ) */
703 (PID.TID 0000.0001) 1.000000000000000E+21
704 (PID.TID 0000.0001) ;
705 (PID.TID 0000.0001) viscA4Grid = /* Grid dependent biharmonic viscosity ( non-dim. ) */
706 (PID.TID 0000.0001) 0.000000000000000E+00
707 (PID.TID 0000.0001) ;
708 (PID.TID 0000.0001) viscC4leith = /* Leith biharm viscosity factor (on grad(vort), non-dim.)*/
709 (PID.TID 0000.0001) 1.500000000000000E+00
710 (PID.TID 0000.0001) ;
711 (PID.TID 0000.0001) viscC4leithD = /* Leith biharm viscosity factor (on grad(div), non-dim.) */
712 (PID.TID 0000.0001) 1.500000000000000E+00
713 (PID.TID 0000.0001) ;
714 (PID.TID 0000.0001) viscC4Smag = /* Smagorinsky biharm viscosity factor (non-dim) */
715 (PID.TID 0000.0001) 0.000000000000000E+00
716 (PID.TID 0000.0001) ;
717 (PID.TID 0000.0001) no_slip_sides = /* Viscous BCs: No-slip sides */
718 (PID.TID 0000.0001) T
719 (PID.TID 0000.0001) ;
720 (PID.TID 0000.0001) sideDragFactor = /* side-drag scaling factor (non-dim) */
721 (PID.TID 0000.0001) 2.000000000000000E+00
722 (PID.TID 0000.0001) ;
723 (PID.TID 0000.0001) viscArNr = /* vertical profile of vertical viscosity ( m^2/s )*/
724 (PID.TID 0000.0001) 15 @ 1.000000000000000E-03 /* K = 1: 15 */
725 (PID.TID 0000.0001) ;
726 (PID.TID 0000.0001) no_slip_bottom = /* Viscous BCs: No-slip bottom */
727 (PID.TID 0000.0001) T
728 (PID.TID 0000.0001) ;
729 (PID.TID 0000.0001) bottomVisc_pCell = /* Partial-cell in bottom Visc. BC */
730 (PID.TID 0000.0001) F
731 (PID.TID 0000.0001) ;
732 (PID.TID 0000.0001) bottomDragLinear = /* linear bottom-drag coefficient ( m/s ) */
733 (PID.TID 0000.0001) 0.000000000000000E+00
734 (PID.TID 0000.0001) ;
735 (PID.TID 0000.0001) bottomDragQuadratic = /* quadratic bottom-drag coefficient (-) */
736 (PID.TID 0000.0001) 2.100000000000000E-03
737 (PID.TID 0000.0001) ;
738 (PID.TID 0000.0001) selectBotDragQuadr = /* select quadratic bottom drag options */
739 (PID.TID 0000.0001) 0
740 (PID.TID 0000.0001) ;
741 (PID.TID 0000.0001) diffKhT = /* Laplacian diffusion of heat laterally ( m^2/s ) */
742 (PID.TID 0000.0001) 0.000000000000000E+00
743 (PID.TID 0000.0001) ;
744 (PID.TID 0000.0001) diffK4T = /* Biharmonic diffusion of heat laterally ( m^4/s ) */
745 (PID.TID 0000.0001) 0.000000000000000E+00
746 (PID.TID 0000.0001) ;
747 (PID.TID 0000.0001) diffKhS = /* Laplacian diffusion of salt laterally ( m^2/s ) */
748 (PID.TID 0000.0001) 0.000000000000000E+00
749 (PID.TID 0000.0001) ;
750 (PID.TID 0000.0001) diffK4S = /* Biharmonic diffusion of salt laterally ( m^4/s ) */
751 (PID.TID 0000.0001) 0.000000000000000E+00
752 (PID.TID 0000.0001) ;
753 (PID.TID 0000.0001) diffKrNrT = /* vertical profile of vertical diffusion of Temp ( m^2/s )*/
754 (PID.TID 0000.0001) 15 @ 0.000000000000000E+00 /* K = 1: 15 */
755 (PID.TID 0000.0001) ;
756 (PID.TID 0000.0001) diffKrNrS = /* vertical profile of vertical diffusion of Salt ( m^2/s )*/
757 (PID.TID 0000.0001) 15 @ 3.000000000000000E-05 /* K = 1: 15 */
758 (PID.TID 0000.0001) ;
759 (PID.TID 0000.0001) diffKrBL79surf = /* Surface diffusion for Bryan and Lewis 79 ( m^2/s ) */
760 (PID.TID 0000.0001) 0.000000000000000E+00
761 (PID.TID 0000.0001) ;
762 (PID.TID 0000.0001) diffKrBL79deep = /* Deep diffusion for Bryan and Lewis 1979 ( m^2/s ) */
763 (PID.TID 0000.0001) 0.000000000000000E+00
764 (PID.TID 0000.0001) ;
765 (PID.TID 0000.0001) diffKrBL79scl = /* Depth scale for Bryan and Lewis 1979 ( m ) */
766 (PID.TID 0000.0001) 2.000000000000000E+02
767 (PID.TID 0000.0001) ;
768 (PID.TID 0000.0001) diffKrBL79Ho = /* Turning depth for Bryan and Lewis 1979 ( m ) */
769 (PID.TID 0000.0001) -2.000000000000000E+03
770 (PID.TID 0000.0001) ;
771 (PID.TID 0000.0001) ivdc_kappa = /* Implicit Vertical Diffusivity for Convection ( m^2/s) */
772 (PID.TID 0000.0001) 1.000000000000000E+01
773 (PID.TID 0000.0001) ;
774 (PID.TID 0000.0001) hMixCriteria= /* Criteria for mixed-layer diagnostic */
775 (PID.TID 0000.0001) -8.000000000000000E-01
776 (PID.TID 0000.0001) ;
777 (PID.TID 0000.0001) dRhoSmall = /* Parameter for mixed-layer diagnostic */
778 (PID.TID 0000.0001) 1.000000000000000E-06
779 (PID.TID 0000.0001) ;
780 (PID.TID 0000.0001) hMixSmooth= /* Smoothing parameter for mixed-layer diagnostic */
781 (PID.TID 0000.0001) 0.000000000000000E+00
782 (PID.TID 0000.0001) ;
783 (PID.TID 0000.0001) eosType = /* Type of Equation of State */
784 (PID.TID 0000.0001) 'JMD95Z'
785 (PID.TID 0000.0001) ;
786 (PID.TID 0000.0001) selectP_inEOS_Zc = /* select pressure to use in EOS (0,1,2,3) */
787 (PID.TID 0000.0001) 0
788 (PID.TID 0000.0001) 0= -g*rhoConst*z ; 1= pRef (from tRef,sRef); 2= Hyd P ; 3= Hyd+NH P
789 (PID.TID 0000.0001) ;
790 (PID.TID 0000.0001) HeatCapacity_Cp = /* Specific heat capacity ( J/kg/K ) */
791 (PID.TID 0000.0001) 3.994000000000000E+03
792 (PID.TID 0000.0001) ;
793 (PID.TID 0000.0001) celsius2K = /* 0 degree Celsius converted to Kelvin ( K ) */
794 (PID.TID 0000.0001) 2.731500000000000E+02
795 (PID.TID 0000.0001) ;
796 (PID.TID 0000.0001) rhoConst = /* Reference density (Boussinesq) ( kg/m^3 ) */
797 (PID.TID 0000.0001) 1.035000000000000E+03
798 (PID.TID 0000.0001) ;
799 (PID.TID 0000.0001) rhoFacC = /* normalized Reference density @ cell-Center (-) */
800 (PID.TID 0000.0001) 15 @ 1.000000000000000E+00 /* K = 1: 15 */
801 (PID.TID 0000.0001) ;
802 (PID.TID 0000.0001) rhoFacF = /* normalized Reference density @ W-Interface (-) */
803 (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
804 (PID.TID 0000.0001) ;
805 (PID.TID 0000.0001) rhoConstFresh = /* Fresh-water reference density ( kg/m^3 ) */
806 (PID.TID 0000.0001) 1.035000000000000E+03
807 (PID.TID 0000.0001) ;
808 (PID.TID 0000.0001) gravity = /* Gravitational acceleration ( m/s^2 ) */
809 (PID.TID 0000.0001) 9.810000000000000E+00
810 (PID.TID 0000.0001) ;
811 (PID.TID 0000.0001) gBaro = /* Barotropic gravity ( m/s^2 ) */
812 (PID.TID 0000.0001) 9.810000000000000E+00
813 (PID.TID 0000.0001) ;
814 (PID.TID 0000.0001) gravFacC = /* gravity factor (vs surf.) @ cell-Center (-) */
815 (PID.TID 0000.0001) 15 @ 1.000000000000000E+00 /* K = 1: 15 */
816 (PID.TID 0000.0001) ;
817 (PID.TID 0000.0001) gravFacF = /* gravity factor (vs surf.) @ W-Interface (-) */
818 (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
819 (PID.TID 0000.0001) ;
820 (PID.TID 0000.0001) rotationPeriod = /* Rotation Period ( s ) */
821 (PID.TID 0000.0001) 8.616400000000000E+04
822 (PID.TID 0000.0001) ;
823 (PID.TID 0000.0001) omega = /* Angular velocity ( rad/s ) */
824 (PID.TID 0000.0001) 7.292123516990375E-05
825 (PID.TID 0000.0001) ;
826 (PID.TID 0000.0001) f0 = /* Reference coriolis parameter ( 1/s ) */
827 (PID.TID 0000.0001) 1.000000000000000E-04
828 (PID.TID 0000.0001) ;
829 (PID.TID 0000.0001) beta = /* Beta ( 1/(m.s) ) */
830 (PID.TID 0000.0001) 9.999999999999999E-12
831 (PID.TID 0000.0001) ;
832 (PID.TID 0000.0001) fPrime = /* Second coriolis parameter ( 1/s ) */
833 (PID.TID 0000.0001) 0.000000000000000E+00
834 (PID.TID 0000.0001) ;
835 (PID.TID 0000.0001) rigidLid = /* Rigid lid on/off flag */
836 (PID.TID 0000.0001) F
837 (PID.TID 0000.0001) ;
838 (PID.TID 0000.0001) implicitFreeSurface = /* Implicit free surface on/off flag */
839 (PID.TID 0000.0001) T
840 (PID.TID 0000.0001) ;
841 (PID.TID 0000.0001) freeSurfFac = /* Implicit free surface factor */
842 (PID.TID 0000.0001) 1.000000000000000E+00
843 (PID.TID 0000.0001) ;
844 (PID.TID 0000.0001) implicSurfPress = /* Surface Pressure implicit factor (0-1)*/
845 (PID.TID 0000.0001) 1.000000000000000E+00
846 (PID.TID 0000.0001) ;
847 (PID.TID 0000.0001) implicDiv2Dflow = /* Barot. Flow Div. implicit factor (0-1)*/
848 (PID.TID 0000.0001) 1.000000000000000E+00
849 (PID.TID 0000.0001) ;
850 (PID.TID 0000.0001) uniformLin_PhiSurf = /* use uniform Bo_surf on/off flag*/
851 (PID.TID 0000.0001) T
852 (PID.TID 0000.0001) ;
853 (PID.TID 0000.0001) uniformFreeSurfLev = /* free-surface level-index is uniform */
854 (PID.TID 0000.0001) T
855 (PID.TID 0000.0001) ;
856 (PID.TID 0000.0001) hFacMin = /* minimum partial cell factor (hFac) */
857 (PID.TID 0000.0001) 5.000000000000000E-02
858 (PID.TID 0000.0001) ;
859 (PID.TID 0000.0001) hFacMinDr = /* minimum partial cell thickness ( m) */
860 (PID.TID 0000.0001) 5.000000000000000E+01
861 (PID.TID 0000.0001) ;
862 (PID.TID 0000.0001) exactConserv = /* Exact Volume Conservation on/off flag*/
863 (PID.TID 0000.0001) T
864 (PID.TID 0000.0001) ;
865 (PID.TID 0000.0001) linFSConserveTr = /* Tracer correction for Lin Free Surface on/off flag*/
866 (PID.TID 0000.0001) F
867 (PID.TID 0000.0001) ;
868 (PID.TID 0000.0001) nonlinFreeSurf = /* Non-linear Free Surf. options (-1,0,1,2,3)*/
869 (PID.TID 0000.0001) 2
870 (PID.TID 0000.0001) -1,0= Off ; 1,2,3= On, 2=+rescale gU,gV, 3=+update cg2d solv.
871 (PID.TID 0000.0001) ;
872 (PID.TID 0000.0001) hFacInf = /* lower threshold for hFac (nonlinFreeSurf only)*/
873 (PID.TID 0000.0001) 2.000000000000000E-01
874 (PID.TID 0000.0001) ;
875 (PID.TID 0000.0001) hFacSup = /* upper threshold for hFac (nonlinFreeSurf only)*/
876 (PID.TID 0000.0001) 2.000000000000000E+00
877 (PID.TID 0000.0001) ;
878 (PID.TID 0000.0001) select_rStar = /* r* Vertical coord. options (=0 r coord.; >0 uses r*)*/
879 (PID.TID 0000.0001) 0
880 (PID.TID 0000.0001) ;
881 (PID.TID 0000.0001) useRealFreshWaterFlux = /* Real Fresh Water Flux on/off flag*/
882 (PID.TID 0000.0001) T
883 (PID.TID 0000.0001) ;
884 (PID.TID 0000.0001) temp_EvPrRn = /* Temp. of Evap/Prec/R (UNSET=use local T)(oC)*/
885 (PID.TID 0000.0001) 1.234567000000000E+05
886 (PID.TID 0000.0001) ;
887 (PID.TID 0000.0001) salt_EvPrRn = /* Salin. of Evap/Prec/R (UNSET=use local S)(psu)*/
888 (PID.TID 0000.0001) 0.000000000000000E+00
889 (PID.TID 0000.0001) ;
890 (PID.TID 0000.0001) selectAddFluid = /* option for mass source/sink of fluid (=0: off) */
891 (PID.TID 0000.0001) 0
892 (PID.TID 0000.0001) ;
893 (PID.TID 0000.0001) temp_addMass = /* Temp. of addMass array (UNSET=use local T)(oC)*/
894 (PID.TID 0000.0001) 1.234567000000000E+05
895 (PID.TID 0000.0001) ;
896 (PID.TID 0000.0001) salt_addMass = /* Salin. of addMass array (UNSET=use local S)(psu)*/
897 (PID.TID 0000.0001) 0.000000000000000E+00
898 (PID.TID 0000.0001) ;
899 (PID.TID 0000.0001) use3Dsolver = /* use 3-D pressure solver on/off flag */
900 (PID.TID 0000.0001) F
901 (PID.TID 0000.0001) ;
902 (PID.TID 0000.0001) nonHydrostatic = /* Non-Hydrostatic on/off flag */
903 (PID.TID 0000.0001) F
904 (PID.TID 0000.0001) ;
905 (PID.TID 0000.0001) nh_Am2 = /* Non-Hydrostatic terms scaling factor */
906 (PID.TID 0000.0001) 1.000000000000000E+00
907 (PID.TID 0000.0001) ;
908 (PID.TID 0000.0001) implicitNHPress = /* Non-Hyd Pressure implicit factor (0-1)*/
909 (PID.TID 0000.0001) 1.000000000000000E+00
910 (PID.TID 0000.0001) ;
911 (PID.TID 0000.0001) selectNHfreeSurf = /* Non-Hyd (free-)Surface option */
912 (PID.TID 0000.0001) 0
913 (PID.TID 0000.0001) ;
914 (PID.TID 0000.0001) quasiHydrostatic = /* Quasi-Hydrostatic on/off flag */
915 (PID.TID 0000.0001) F
916 (PID.TID 0000.0001) ;
917 (PID.TID 0000.0001) calc_wVelocity = /* vertical velocity calculation on/off flag */
918 (PID.TID 0000.0001) T
919 (PID.TID 0000.0001) ;
920 (PID.TID 0000.0001) momStepping = /* Momentum equation on/off flag */
921 (PID.TID 0000.0001) T
922 (PID.TID 0000.0001) ;
923 (PID.TID 0000.0001) vectorInvariantMomentum= /* Vector-Invariant Momentum on/off */
924 (PID.TID 0000.0001) T
925 (PID.TID 0000.0001) ;
926 (PID.TID 0000.0001) momAdvection = /* Momentum advection on/off flag */
927 (PID.TID 0000.0001) T
928 (PID.TID 0000.0001) ;
929 (PID.TID 0000.0001) momViscosity = /* Momentum viscosity on/off flag */
930 (PID.TID 0000.0001) T
931 (PID.TID 0000.0001) ;
932 (PID.TID 0000.0001) momImplVertAdv= /* Momentum implicit vert. advection on/off*/
933 (PID.TID 0000.0001) F
934 (PID.TID 0000.0001) ;
935 (PID.TID 0000.0001) implicitViscosity = /* Implicit viscosity on/off flag */
936 (PID.TID 0000.0001) F
937 (PID.TID 0000.0001) ;
938 (PID.TID 0000.0001) implBottomFriction= /* Implicit bottom friction on/off flag */
939 (PID.TID 0000.0001) F
940 (PID.TID 0000.0001) ;
941 (PID.TID 0000.0001) metricTerms = /* metric-Terms on/off flag */
942 (PID.TID 0000.0001) T
943 (PID.TID 0000.0001) ;
944 (PID.TID 0000.0001) useNHMTerms = /* Non-Hydrostatic Metric-Terms on/off */
945 (PID.TID 0000.0001) F
946 (PID.TID 0000.0001) ;
947 (PID.TID 0000.0001) selectCoriMap = /* Coriolis Map options (0,1,2,3)*/
948 (PID.TID 0000.0001) 2
949 (PID.TID 0000.0001) 0= f-Plane ; 1= Beta-Plane ; 2= Spherical ; 3= read from file
950 (PID.TID 0000.0001) ;
951 (PID.TID 0000.0001) use3dCoriolis = /* 3-D Coriolis on/off flag */
952 (PID.TID 0000.0001) F
953 (PID.TID 0000.0001) ;
954 (PID.TID 0000.0001) useCoriolis = /* Coriolis on/off flag */
955 (PID.TID 0000.0001) T
956 (PID.TID 0000.0001) ;
957 (PID.TID 0000.0001) useCDscheme = /* CD scheme on/off flag */
958 (PID.TID 0000.0001) T
959 (PID.TID 0000.0001) ;
960 (PID.TID 0000.0001) useEnergyConservingCoriolis= /* Flx-Form Coriolis scheme flag */
961 (PID.TID 0000.0001) F
962 (PID.TID 0000.0001) ;
963 (PID.TID 0000.0001) useJamartWetPoints= /* Coriolis WetPoints method flag */
964 (PID.TID 0000.0001) F
965 (PID.TID 0000.0001) ;
966 (PID.TID 0000.0001) useJamartMomAdv= /* V.I Non-linear terms Jamart flag */
967 (PID.TID 0000.0001) F
968 (PID.TID 0000.0001) ;
969 (PID.TID 0000.0001) useAbsVorticity= /* V.I Works with f+zeta in Coriolis */
970 (PID.TID 0000.0001) F
971 (PID.TID 0000.0001) ;
972 (PID.TID 0000.0001) selectVortScheme= /* V.I Scheme selector for Vorticity-Term */
973 (PID.TID 0000.0001) 1
974 (PID.TID 0000.0001) = 0 : enstrophy (Shallow-Water Eq.) conserving scheme by Sadourny, JAS 75
975 (PID.TID 0000.0001) = 1 : same as 0 with modified hFac
976 (PID.TID 0000.0001) = 2 : energy conserving scheme (used by Sadourny in JAS 75 paper)
977 (PID.TID 0000.0001) = 3 : energy (general) and enstrophy (2D, nonDiv.) conserving scheme
978 (PID.TID 0000.0001) from Sadourny (Burridge & Haseler, ECMWF Rep.4, 1977)
979 (PID.TID 0000.0001) ;
980 (PID.TID 0000.0001) upwindVorticity= /* V.I Upwind bias vorticity flag */
981 (PID.TID 0000.0001) F
982 (PID.TID 0000.0001) ;
983 (PID.TID 0000.0001) highOrderVorticity= /* V.I High order vort. advect. flag */
984 (PID.TID 0000.0001) F
985 (PID.TID 0000.0001) ;
986 (PID.TID 0000.0001) upwindShear= /* V.I Upwind vertical Shear advection flag */
987 (PID.TID 0000.0001) F
988 (PID.TID 0000.0001) ;
989 (PID.TID 0000.0001) selectKEscheme= /* V.I Kinetic Energy scheme selector */
990 (PID.TID 0000.0001) 0
991 (PID.TID 0000.0001) ;
992 (PID.TID 0000.0001) momForcing = /* Momentum forcing on/off flag */
993 (PID.TID 0000.0001) T
994 (PID.TID 0000.0001) ;
995 (PID.TID 0000.0001) momPressureForcing = /* Momentum pressure term on/off flag */
996 (PID.TID 0000.0001) T
997 (PID.TID 0000.0001) ;
998 (PID.TID 0000.0001) implicitIntGravWave= /* Implicit Internal Gravity Wave flag */
999 (PID.TID 0000.0001) F
1000 (PID.TID 0000.0001) ;
1001 (PID.TID 0000.0001) staggerTimeStep = /* Stagger time stepping on/off flag */
1002 (PID.TID 0000.0001) T
1003 (PID.TID 0000.0001) ;
1004 (PID.TID 0000.0001) doResetHFactors = /* reset thickness factors @ each time-step */
1005 (PID.TID 0000.0001) T
1006 (PID.TID 0000.0001) ;
1007 (PID.TID 0000.0001) multiDimAdvection = /* enable/disable Multi-Dim Advection */
1008 (PID.TID 0000.0001) T
1009 (PID.TID 0000.0001) ;
1010 (PID.TID 0000.0001) useMultiDimAdvec = /* Multi-Dim Advection is/is-not used */
1011 (PID.TID 0000.0001) T
1012 (PID.TID 0000.0001) ;
1013 (PID.TID 0000.0001) implicitDiffusion = /* Implicit Diffusion on/off flag */
1014 (PID.TID 0000.0001) T
1015 (PID.TID 0000.0001) ;
1016 (PID.TID 0000.0001) tempStepping = /* Temperature equation on/off flag */
1017 (PID.TID 0000.0001) T
1018 (PID.TID 0000.0001) ;
1019 (PID.TID 0000.0001) tempAdvection = /* Temperature advection on/off flag */
1020 (PID.TID 0000.0001) T
1021 (PID.TID 0000.0001) ;
1022 (PID.TID 0000.0001) tempImplVertAdv = /* Temp. implicit vert. advection on/off */
1023 (PID.TID 0000.0001) T
1024 (PID.TID 0000.0001) ;
1025 (PID.TID 0000.0001) tempForcing = /* Temperature forcing on/off flag */
1026 (PID.TID 0000.0001) T
1027 (PID.TID 0000.0001) ;
1028 (PID.TID 0000.0001) doThetaClimRelax = /* apply SST relaxation on/off flag */
1029 (PID.TID 0000.0001) T
1030 (PID.TID 0000.0001) ;
1031 (PID.TID 0000.0001) tempIsActiveTr = /* Temp. is a dynamically Active Tracer */
1032 (PID.TID 0000.0001) T
1033 (PID.TID 0000.0001) ;
1034 (PID.TID 0000.0001) saltStepping = /* Salinity equation on/off flag */
1035 (PID.TID 0000.0001) T
1036 (PID.TID 0000.0001) ;
1037 (PID.TID 0000.0001) saltAdvection = /* Salinity advection on/off flag */
1038 (PID.TID 0000.0001) T
1039 (PID.TID 0000.0001) ;
1040 (PID.TID 0000.0001) saltImplVertAdv = /* Sali. implicit vert. advection on/off */
1041 (PID.TID 0000.0001) T
1042 (PID.TID 0000.0001) ;
1043 (PID.TID 0000.0001) saltForcing = /* Salinity forcing on/off flag */
1044 (PID.TID 0000.0001) T
1045 (PID.TID 0000.0001) ;
1046 (PID.TID 0000.0001) doSaltClimRelax = /* apply SSS relaxation on/off flag */
1047 (PID.TID 0000.0001) T
1048 (PID.TID 0000.0001) ;
1049 (PID.TID 0000.0001) saltIsActiveTr = /* Salt is a dynamically Active Tracer */
1050 (PID.TID 0000.0001) T
1051 (PID.TID 0000.0001) ;
1052 (PID.TID 0000.0001) readBinaryPrec = /* Precision used for reading binary files */
1053 (PID.TID 0000.0001) 32
1054 (PID.TID 0000.0001) ;
1055 (PID.TID 0000.0001) writeBinaryPrec = /* Precision used for writing binary files */
1056 (PID.TID 0000.0001) 32
1057 (PID.TID 0000.0001) ;
1058 (PID.TID 0000.0001) globalFiles = /* write "global" (=not per tile) files */
1059 (PID.TID 0000.0001) F
1060 (PID.TID 0000.0001) ;
1061 (PID.TID 0000.0001) useSingleCpuIO = /* only master MPI process does I/O */
1062 (PID.TID 0000.0001) F
1063 (PID.TID 0000.0001) ;
1064 (PID.TID 0000.0001) useSingleCpuInput = /* only master process reads input */
1065 (PID.TID 0000.0001) F
1066 (PID.TID 0000.0001) ;
1067 (PID.TID 0000.0001) /* debLev[*] : level of debug & auxiliary message printing */
1068 (PID.TID 0000.0001) debLevZero = 0 ; /* level of disabled aux. msg printing */
1069 (PID.TID 0000.0001) debLevA = 1 ; /* level of minimum aux. msg printing */
1070 (PID.TID 0000.0001) debLevB = 2 ; /* level of low aux. print (report read-file opening)*/
1071 (PID.TID 0000.0001) debLevC = 3 ; /* level of moderate debug prt (most pkgs debug msg) */
1072 (PID.TID 0000.0001) debLevD = 4 ; /* level of enhanced debug prt (add DEBUG_STATS prt) */
1073 (PID.TID 0000.0001) debLevE = 5 ; /* level of extensive debug printing */
1074 (PID.TID 0000.0001) debugLevel = /* select debug printing level */
1075 (PID.TID 0000.0001) 1
1076 (PID.TID 0000.0001) ;
1077 (PID.TID 0000.0001) //
1078 (PID.TID 0000.0001) // Elliptic solver(s) paramters ( PARM02 in namelist )
1079 (PID.TID 0000.0001) //
1080 (PID.TID 0000.0001) cg2dMaxIters = /* Upper limit on 2d con. grad iterations */
1081 (PID.TID 0000.0001) 500
1082 (PID.TID 0000.0001) ;
1083 (PID.TID 0000.0001) cg2dChkResFreq = /* 2d con. grad convergence test frequency */
1084 (PID.TID 0000.0001) 1
1085 (PID.TID 0000.0001) ;
1086 (PID.TID 0000.0001) cg2dUseMinResSol= /* use cg2d last-iter(=0) / min-resid.(=1) solution */
1087 (PID.TID 0000.0001) 0
1088 (PID.TID 0000.0001) ;
1089 (PID.TID 0000.0001) cg2dTargetResidual = /* 2d con. grad target residual */
1090 (PID.TID 0000.0001) 1.000000000000000E-13
1091 (PID.TID 0000.0001) ;
1092 (PID.TID 0000.0001) cg2dTargetResWunit = /* CG2d target residual [W units] */
1093 (PID.TID 0000.0001) -1.000000000000000E+00
1094 (PID.TID 0000.0001) ;
1095 (PID.TID 0000.0001) cg2dPreCondFreq = /* Freq. for updating cg2d preconditioner */
1096 (PID.TID 0000.0001) 1
1097 (PID.TID 0000.0001) ;
1098 (PID.TID 0000.0001) useSRCGSolver = /* use single reduction CG solver(s) */
1099 (PID.TID 0000.0001) F
1100 (PID.TID 0000.0001) ;
1101 (PID.TID 0000.0001) printResidualFreq = /* Freq. for printing CG residual */
1102 (PID.TID 0000.0001) 0
1103 (PID.TID 0000.0001) ;
1104 (PID.TID 0000.0001) //
1105 (PID.TID 0000.0001) // Time stepping paramters ( PARM03 in namelist )
1106 (PID.TID 0000.0001) //
1107 (PID.TID 0000.0001) deltaTMom = /* Momentum equation timestep ( s ) */
1108 (PID.TID 0000.0001) 1.200000000000000E+03
1109 (PID.TID 0000.0001) ;
1110 (PID.TID 0000.0001) deltaTFreeSurf = /* FreeSurface equation timestep ( s ) */
1111 (PID.TID 0000.0001) 1.200000000000000E+03
1112 (PID.TID 0000.0001) ;
1113 (PID.TID 0000.0001) dTtracerLev = /* Tracer equation timestep ( s ) */
1114 (PID.TID 0000.0001) 15 @ 4.320000000000000E+04 /* K = 1: 15 */
1115 (PID.TID 0000.0001) ;
1116 (PID.TID 0000.0001) deltaTClock = /* Model clock timestep ( s ) */
1117 (PID.TID 0000.0001) 4.320000000000000E+04
1118 (PID.TID 0000.0001) ;
1119 (PID.TID 0000.0001) cAdjFreq = /* Convective adjustment interval ( s ) */
1120 (PID.TID 0000.0001) 0.000000000000000E+00
1121 (PID.TID 0000.0001) ;
1122 (PID.TID 0000.0001) momForcingOutAB = /* =1: take Momentum Forcing out of Adams-Bash. stepping */
1123 (PID.TID 0000.0001) 0
1124 (PID.TID 0000.0001) ;
1125 (PID.TID 0000.0001) tracForcingOutAB = /* =1: take T,S,pTr Forcing out of Adams-Bash. stepping */
1126 (PID.TID 0000.0001) 0
1127 (PID.TID 0000.0001) ;
1128 (PID.TID 0000.0001) momDissip_In_AB = /* put Dissipation Tendency in Adams-Bash. stepping */
1129 (PID.TID 0000.0001) T
1130 (PID.TID 0000.0001) ;
1131 (PID.TID 0000.0001) doAB_onGtGs = /* apply AB on Tendencies (rather than on T,S)*/
1132 (PID.TID 0000.0001) T
1133 (PID.TID 0000.0001) ;
1134 (PID.TID 0000.0001) abEps = /* Adams-Bashforth-2 stabilizing weight */
1135 (PID.TID 0000.0001) 1.000000000000000E-01
1136 (PID.TID 0000.0001) ;
1137 (PID.TID 0000.0001) applyExchUV_early = /* Apply EXCH to U,V earlier in time-step */
1138 (PID.TID 0000.0001) F
1139 (PID.TID 0000.0001) ;
1140 (PID.TID 0000.0001) tauCD = /* CD coupling time-scale ( s ) */
1141 (PID.TID 0000.0001) 3.214280000000000E+05
1142 (PID.TID 0000.0001) ;
1143 (PID.TID 0000.0001) rCD = /* Normalised CD coupling parameter */
1144 (PID.TID 0000.0001) 9.962666600296178E-01
1145 (PID.TID 0000.0001) ;
1146 (PID.TID 0000.0001) epsAB_CD = /* AB-2 stabilizing weight for CD-scheme*/
1147 (PID.TID 0000.0001) 1.000000000000000E-01
1148 (PID.TID 0000.0001) ;
1149 (PID.TID 0000.0001) pickupStrictlyMatch= /* stop if pickup do not strictly match */
1150 (PID.TID 0000.0001) T
1151 (PID.TID 0000.0001) ;
1152 (PID.TID 0000.0001) nIter0 = /* Run starting timestep number */
1153 (PID.TID 0000.0001) 0
1154 (PID.TID 0000.0001) ;
1155 (PID.TID 0000.0001) nTimeSteps = /* Number of timesteps */
1156 (PID.TID 0000.0001) 10
1157 (PID.TID 0000.0001) ;
1158 (PID.TID 0000.0001) nEndIter = /* Run ending timestep number */
1159 (PID.TID 0000.0001) 10
1160 (PID.TID 0000.0001) ;
1161 (PID.TID 0000.0001) baseTime = /* Model base time ( s ) */
1162 (PID.TID 0000.0001) 0.000000000000000E+00
1163 (PID.TID 0000.0001) ;
1164 (PID.TID 0000.0001) startTime = /* Run start time ( s ) */
1165 (PID.TID 0000.0001) 0.000000000000000E+00
1166 (PID.TID 0000.0001) ;
1167 (PID.TID 0000.0001) endTime = /* Integration ending time ( s ) */
1168 (PID.TID 0000.0001) 4.320000000000000E+05
1169 (PID.TID 0000.0001) ;
1170 (PID.TID 0000.0001) pChkPtFreq = /* Permanent restart/pickup file interval ( s ) */
1171 (PID.TID 0000.0001) 3.110400000000000E+08
1172 (PID.TID 0000.0001) ;
1173 (PID.TID 0000.0001) chkPtFreq = /* Rolling restart/pickup file interval ( s ) */
1174 (PID.TID 0000.0001) 0.000000000000000E+00
1175 (PID.TID 0000.0001) ;
1176 (PID.TID 0000.0001) pickup_write_mdsio = /* Model IO flag. */
1177 (PID.TID 0000.0001) T
1178 (PID.TID 0000.0001) ;
1179 (PID.TID 0000.0001) pickup_read_mdsio = /* Model IO flag. */
1180 (PID.TID 0000.0001) T
1181 (PID.TID 0000.0001) ;
1182 (PID.TID 0000.0001) pickup_write_immed = /* Model IO flag. */
1183 (PID.TID 0000.0001) F
1184 (PID.TID 0000.0001) ;
1185 (PID.TID 0000.0001) writePickupAtEnd = /* Model IO flag. */
1186 (PID.TID 0000.0001) T
1187 (PID.TID 0000.0001) ;
1188 (PID.TID 0000.0001) dumpFreq = /* Model state write out interval ( s ). */
1189 (PID.TID 0000.0001) 8.640000000000000E+05
1190 (PID.TID 0000.0001) ;
1191 (PID.TID 0000.0001) dumpInitAndLast= /* write out Initial & Last iter. model state */
1192 (PID.TID 0000.0001) T
1193 (PID.TID 0000.0001) ;
1194 (PID.TID 0000.0001) snapshot_mdsio = /* Model IO flag. */
1195 (PID.TID 0000.0001) T
1196 (PID.TID 0000.0001) ;
1197 (PID.TID 0000.0001) monitorFreq = /* Monitor output interval ( s ). */
1198 (PID.TID 0000.0001) 0.000000000000000E+00
1199 (PID.TID 0000.0001) ;
1200 (PID.TID 0000.0001) monitorSelect = /* select group of variables to monitor */
1201 (PID.TID 0000.0001) 3
1202 (PID.TID 0000.0001) ;
1203 (PID.TID 0000.0001) monitor_stdio = /* Model IO flag. */
1204 (PID.TID 0000.0001) T
1205 (PID.TID 0000.0001) ;
1206 (PID.TID 0000.0001) externForcingPeriod = /* forcing period (s) */
1207 (PID.TID 0000.0001) 2.592000000000000E+06
1208 (PID.TID 0000.0001) ;
1209 (PID.TID 0000.0001) externForcingCycle = /* period of the cyle (s). */
1210 (PID.TID 0000.0001) 3.110400000000000E+07
1211 (PID.TID 0000.0001) ;
1212 (PID.TID 0000.0001) tauThetaClimRelax = /* relaxation time scale (s) */
1213 (PID.TID 0000.0001) 5.184000000000000E+06
1214 (PID.TID 0000.0001) ;
1215 (PID.TID 0000.0001) tauSaltClimRelax = /* relaxation time scale (s) */
1216 (PID.TID 0000.0001) 1.555200000000000E+07
1217 (PID.TID 0000.0001) ;
1218 (PID.TID 0000.0001) latBandClimRelax = /* max. Lat. where relaxation */
1219 (PID.TID 0000.0001) 1.800000000000000E+02
1220 (PID.TID 0000.0001) ;
1221 (PID.TID 0000.0001) //
1222 (PID.TID 0000.0001) // Gridding paramters ( PARM04 in namelist )
1223 (PID.TID 0000.0001) //
1224 (PID.TID 0000.0001) usingCartesianGrid = /* Cartesian coordinates flag ( True/False ) */
1225 (PID.TID 0000.0001) F
1226 (PID.TID 0000.0001) ;
1227 (PID.TID 0000.0001) usingCylindricalGrid = /* Cylindrical coordinates flag ( True/False ) */
1228 (PID.TID 0000.0001) F
1229 (PID.TID 0000.0001) ;
1230 (PID.TID 0000.0001) usingSphericalPolarGrid = /* Spherical coordinates flag ( True/False ) */
1231 (PID.TID 0000.0001) T
1232 (PID.TID 0000.0001) ;
1233 (PID.TID 0000.0001) usingCurvilinearGrid = /* Curvilinear coordinates flag ( True/False ) */
1234 (PID.TID 0000.0001) F
1235 (PID.TID 0000.0001) ;
1236 (PID.TID 0000.0001) selectSigmaCoord = /* Hybrid-Sigma Vert. Coordinate option */
1237 (PID.TID 0000.0001) 0
1238 (PID.TID 0000.0001) ;
1239 (PID.TID 0000.0001) rSigmaBnd = /* r/sigma transition ( units of r == m ) */
1240 (PID.TID 0000.0001) 1.234567000000000E+05
1241 (PID.TID 0000.0001) ;
1242 (PID.TID 0000.0001) rkSign = /* index orientation relative to vertical coordinate */
1243 (PID.TID 0000.0001) -1.000000000000000E+00
1244 (PID.TID 0000.0001) ;
1245 (PID.TID 0000.0001) gravitySign = /* gravity orientation relative to vertical coordinate */
1246 (PID.TID 0000.0001) -1.000000000000000E+00
1247 (PID.TID 0000.0001) ;
1248 (PID.TID 0000.0001) seaLev_Z = /* reference height of sea-level [m] */
1249 (PID.TID 0000.0001) 0.000000000000000E+00
1250 (PID.TID 0000.0001) ;
1251 (PID.TID 0000.0001) top_Pres = /* reference pressure at the top [Pa] */
1252 (PID.TID 0000.0001) 0.000000000000000E+00
1253 (PID.TID 0000.0001) ;
1254 (PID.TID 0000.0001) mass2rUnit = /* convert mass per unit area [kg/m2] to r-units [m] */
1255 (PID.TID 0000.0001) 9.661835748792270E-04
1256 (PID.TID 0000.0001) ;
1257 (PID.TID 0000.0001) rUnit2mass = /* convert r-units [m] to mass per unit area [kg/m2] */
1258 (PID.TID 0000.0001) 1.035000000000000E+03
1259 (PID.TID 0000.0001) ;
1260 (PID.TID 0000.0001) drC = /* C spacing ( units of r ) */
1261 (PID.TID 0000.0001) 2.500000000000000E+01, /* K = 1 */
1262 (PID.TID 0000.0001) 6.000000000000000E+01, /* K = 2 */
1263 (PID.TID 0000.0001) 8.500000000000000E+01, /* K = 3 */
1264 (PID.TID 0000.0001) 1.200000000000000E+02, /* K = 4 */
1265 (PID.TID 0000.0001) 1.650000000000000E+02, /* K = 5 */
1266 (PID.TID 0000.0001) 2.150000000000000E+02, /* K = 6 */
1267 (PID.TID 0000.0001) 2.650000000000000E+02, /* K = 7 */
1268 (PID.TID 0000.0001) 3.150000000000000E+02, /* K = 8 */
1269 (PID.TID 0000.0001) 3.650000000000000E+02, /* K = 9 */
1270 (PID.TID 0000.0001) 4.150000000000000E+02, /* K = 10 */
1271 (PID.TID 0000.0001) 4.650000000000000E+02, /* K = 11 */
1272 (PID.TID 0000.0001) 5.150000000000000E+02, /* K = 12 */
1273 (PID.TID 0000.0001) 5.650000000000000E+02, /* K = 13 */
1274 (PID.TID 0000.0001) 6.150000000000000E+02, /* K = 14 */
1275 (PID.TID 0000.0001) 6.650000000000000E+02, /* K = 15 */
1276 (PID.TID 0000.0001) 3.450000000000000E+02 /* K = 16 */
1277 (PID.TID 0000.0001) ;
1278 (PID.TID 0000.0001) drF = /* W spacing ( units of r ) */
1279 (PID.TID 0000.0001) 5.000000000000000E+01, /* K = 1 */
1280 (PID.TID 0000.0001) 7.000000000000000E+01, /* K = 2 */
1281 (PID.TID 0000.0001) 1.000000000000000E+02, /* K = 3 */
1282 (PID.TID 0000.0001) 1.400000000000000E+02, /* K = 4 */
1283 (PID.TID 0000.0001) 1.900000000000000E+02, /* K = 5 */
1284 (PID.TID 0000.0001) 2.400000000000000E+02, /* K = 6 */
1285 (PID.TID 0000.0001) 2.900000000000000E+02, /* K = 7 */
1286 (PID.TID 0000.0001) 3.400000000000000E+02, /* K = 8 */
1287 (PID.TID 0000.0001) 3.900000000000000E+02, /* K = 9 */
1288 (PID.TID 0000.0001) 4.400000000000000E+02, /* K = 10 */
1289 (PID.TID 0000.0001) 4.900000000000000E+02, /* K = 11 */
1290 (PID.TID 0000.0001) 5.400000000000000E+02, /* K = 12 */
1291 (PID.TID 0000.0001) 5.900000000000000E+02, /* K = 13 */
1292 (PID.TID 0000.0001) 6.400000000000000E+02, /* K = 14 */
1293 (PID.TID 0000.0001) 6.900000000000000E+02 /* K = 15 */
1294 (PID.TID 0000.0001) ;
1295 (PID.TID 0000.0001) delX = /* U spacing ( m - cartesian, degrees - spherical ) */
1296 (PID.TID 0000.0001) 90 @ 4.000000000000000E+00 /* I = 1: 90 */
1297 (PID.TID 0000.0001) ;
1298 (PID.TID 0000.0001) delY = /* V spacing ( m - cartesian, degrees - spherical ) */
1299 (PID.TID 0000.0001) 40 @ 4.000000000000000E+00 /* J = 1: 40 */
1300 (PID.TID 0000.0001) ;
1301 (PID.TID 0000.0001) xgOrigin = /* X-axis origin of West edge (cartesian: m, lat-lon: deg) */
1302 (PID.TID 0000.0001) 0.000000000000000E+00
1303 (PID.TID 0000.0001) ;
1304 (PID.TID 0000.0001) ygOrigin = /* Y-axis origin of South edge (cartesian: m, lat-lon: deg) */
1305 (PID.TID 0000.0001) -8.000000000000000E+01
1306 (PID.TID 0000.0001) ;
1307 (PID.TID 0000.0001) rSphere = /* Radius ( ignored - cartesian, m - spherical ) */
1308 (PID.TID 0000.0001) 6.370000000000000E+06
1309 (PID.TID 0000.0001) ;
1310 (PID.TID 0000.0001) deepAtmosphere = /* Deep/Shallow Atmosphere flag (True/False) */
1311 (PID.TID 0000.0001) F
1312 (PID.TID 0000.0001) ;
1313 (PID.TID 0000.0001) xC = /* xC(:,1,:,1) : P-point X coord ( deg. or m if cartesian) */
1314 (PID.TID 0000.0001) 2.000000000000000E+00, /* I = 1 */
1315 (PID.TID 0000.0001) 6.000000000000000E+00, /* I = 2 */
1316 (PID.TID 0000.0001) 1.000000000000000E+01, /* I = 3 */
1317 (PID.TID 0000.0001) . . .
1318 (PID.TID 0000.0001) 8.200000000000000E+01, /* I = 21 */
1319 (PID.TID 0000.0001) 8.600000000000000E+01, /* I = 22 */
1320 (PID.TID 0000.0001) 9.000000000000000E+01, /* I = 23 */
1321 (PID.TID 0000.0001) 9.400000000000000E+01, /* I = 24 */
1322 (PID.TID 0000.0001) 9.800000000000000E+01, /* I = 25 */
1323 (PID.TID 0000.0001) 1.020000000000000E+02, /* I = 26 */
1324 (PID.TID 0000.0001) . . .
1325 (PID.TID 0000.0001) 1.700000000000000E+02, /* I = 43 */
1326 (PID.TID 0000.0001) 1.740000000000000E+02, /* I = 44 */
1327 (PID.TID 0000.0001) 1.780000000000000E+02, /* I = 45 */
1328 (PID.TID 0000.0001) 1.820000000000000E+02, /* I = 46 */
1329 (PID.TID 0000.0001) 1.860000000000000E+02, /* I = 47 */
1330 (PID.TID 0000.0001) 1.900000000000000E+02, /* I = 48 */
1331 (PID.TID 0000.0001) . . .
1332 (PID.TID 0000.0001) 2.580000000000000E+02, /* I = 65 */
1333 (PID.TID 0000.0001) 2.620000000000000E+02, /* I = 66 */
1334 (PID.TID 0000.0001) 2.660000000000000E+02, /* I = 67 */
1335 (PID.TID 0000.0001) 2.700000000000000E+02, /* I = 68 */
1336 (PID.TID 0000.0001) 2.740000000000000E+02, /* I = 69 */
1337 (PID.TID 0000.0001) 2.780000000000000E+02, /* I = 70 */
1338 (PID.TID 0000.0001) . . .
1339 (PID.TID 0000.0001) 3.500000000000000E+02, /* I = 88 */
1340 (PID.TID 0000.0001) 3.540000000000000E+02, /* I = 89 */
1341 (PID.TID 0000.0001) 3.580000000000000E+02 /* I = 90 */
1342 (PID.TID 0000.0001) ;
1343 (PID.TID 0000.0001) yC = /* yC(1,:,1,:) : P-point Y coord ( deg. or m if cartesian) */
1344 (PID.TID 0000.0001) -7.800000000000000E+01, /* J = 1 */
1345 (PID.TID 0000.0001) -7.400000000000000E+01, /* J = 2 */
1346 (PID.TID 0000.0001) -7.000000000000000E+01, /* J = 3 */
1347 (PID.TID 0000.0001) -6.600000000000000E+01, /* J = 4 */
1348 (PID.TID 0000.0001) -6.200000000000000E+01, /* J = 5 */
1349 (PID.TID 0000.0001) -5.800000000000000E+01, /* J = 6 */
1350 (PID.TID 0000.0001) -5.400000000000000E+01, /* J = 7 */
1351 (PID.TID 0000.0001) -5.000000000000000E+01, /* J = 8 */
1352 (PID.TID 0000.0001) -4.600000000000000E+01, /* J = 9 */
1353 (PID.TID 0000.0001) -4.200000000000000E+01, /* J = 10 */
1354 (PID.TID 0000.0001) -3.800000000000000E+01, /* J = 11 */
1355 (PID.TID 0000.0001) -3.400000000000000E+01, /* J = 12 */
1356 (PID.TID 0000.0001) -3.000000000000000E+01, /* J = 13 */
1357 (PID.TID 0000.0001) -2.600000000000000E+01, /* J = 14 */
1358 (PID.TID 0000.0001) -2.200000000000000E+01, /* J = 15 */
1359 (PID.TID 0000.0001) -1.800000000000000E+01, /* J = 16 */
1360 (PID.TID 0000.0001) -1.400000000000000E+01, /* J = 17 */
1361 (PID.TID 0000.0001) -1.000000000000000E+01, /* J = 18 */
1362 (PID.TID 0000.0001) -6.000000000000000E+00, /* J = 19 */
1363 (PID.TID 0000.0001) -2.000000000000000E+00, /* J = 20 */
1364 (PID.TID 0000.0001) 2.000000000000000E+00, /* J = 21 */
1365 (PID.TID 0000.0001) 6.000000000000000E+00, /* J = 22 */
1366 (PID.TID 0000.0001) 1.000000000000000E+01, /* J = 23 */
1367 (PID.TID 0000.0001) 1.400000000000000E+01, /* J = 24 */
1368 (PID.TID 0000.0001) 1.800000000000000E+01, /* J = 25 */
1369 (PID.TID 0000.0001) 2.200000000000000E+01, /* J = 26 */
1370 (PID.TID 0000.0001) 2.600000000000000E+01, /* J = 27 */
1371 (PID.TID 0000.0001) 3.000000000000000E+01, /* J = 28 */
1372 (PID.TID 0000.0001) 3.400000000000000E+01, /* J = 29 */
1373 (PID.TID 0000.0001) 3.800000000000000E+01, /* J = 30 */
1374 (PID.TID 0000.0001) 4.200000000000000E+01, /* J = 31 */
1375 (PID.TID 0000.0001) 4.600000000000000E+01, /* J = 32 */
1376 (PID.TID 0000.0001) 5.000000000000000E+01, /* J = 33 */
1377 (PID.TID 0000.0001) 5.400000000000000E+01, /* J = 34 */
1378 (PID.TID 0000.0001) 5.800000000000000E+01, /* J = 35 */
1379 (PID.TID 0000.0001) 6.200000000000000E+01, /* J = 36 */
1380 (PID.TID 0000.0001) 6.600000000000000E+01, /* J = 37 */
1381 (PID.TID 0000.0001) 7.000000000000000E+01, /* J = 38 */
1382 (PID.TID 0000.0001) 7.400000000000000E+01, /* J = 39 */
1383 (PID.TID 0000.0001) 7.800000000000000E+01 /* J = 40 */
1384 (PID.TID 0000.0001) ;
1385 (PID.TID 0000.0001) rcoord = /* P-point R coordinate ( units of r ) */
1386 (PID.TID 0000.0001) -2.500000000000000E+01, /* K = 1 */
1387 (PID.TID 0000.0001) -8.500000000000000E+01, /* K = 2 */
1388 (PID.TID 0000.0001) -1.700000000000000E+02, /* K = 3 */
1389 (PID.TID 0000.0001) -2.900000000000000E+02, /* K = 4 */
1390 (PID.TID 0000.0001) -4.550000000000000E+02, /* K = 5 */
1391 (PID.TID 0000.0001) -6.700000000000000E+02, /* K = 6 */
1392 (PID.TID 0000.0001) -9.350000000000000E+02, /* K = 7 */
1393 (PID.TID 0000.0001) -1.250000000000000E+03, /* K = 8 */
1394 (PID.TID 0000.0001) -1.615000000000000E+03, /* K = 9 */
1395 (PID.TID 0000.0001) -2.030000000000000E+03, /* K = 10 */
1396 (PID.TID 0000.0001) -2.495000000000000E+03, /* K = 11 */
1397 (PID.TID 0000.0001) -3.010000000000000E+03, /* K = 12 */
1398 (PID.TID 0000.0001) -3.575000000000000E+03, /* K = 13 */
1399 (PID.TID 0000.0001) -4.190000000000000E+03, /* K = 14 */
1400 (PID.TID 0000.0001) -4.855000000000000E+03 /* K = 15 */
1401 (PID.TID 0000.0001) ;
1402 (PID.TID 0000.0001) rF = /* W-Interf. R coordinate ( units of r ) */
1403 (PID.TID 0000.0001) 0.000000000000000E+00, /* K = 1 */
1404 (PID.TID 0000.0001) -5.000000000000000E+01, /* K = 2 */
1405 (PID.TID 0000.0001) -1.200000000000000E+02, /* K = 3 */
1406 (PID.TID 0000.0001) -2.200000000000000E+02, /* K = 4 */
1407 (PID.TID 0000.0001) -3.600000000000000E+02, /* K = 5 */
1408 (PID.TID 0000.0001) -5.500000000000000E+02, /* K = 6 */
1409 (PID.TID 0000.0001) -7.900000000000000E+02, /* K = 7 */
1410 (PID.TID 0000.0001) -1.080000000000000E+03, /* K = 8 */
1411 (PID.TID 0000.0001) -1.420000000000000E+03, /* K = 9 */
1412 (PID.TID 0000.0001) -1.810000000000000E+03, /* K = 10 */
1413 (PID.TID 0000.0001) -2.250000000000000E+03, /* K = 11 */
1414 (PID.TID 0000.0001) -2.740000000000000E+03, /* K = 12 */
1415 (PID.TID 0000.0001) -3.280000000000000E+03, /* K = 13 */
1416 (PID.TID 0000.0001) -3.870000000000000E+03, /* K = 14 */
1417 (PID.TID 0000.0001) -4.510000000000000E+03, /* K = 15 */
1418 (PID.TID 0000.0001) -5.200000000000000E+03 /* K = 16 */
1419 (PID.TID 0000.0001) ;
1420 (PID.TID 0000.0001) deepFacC = /* deep-model grid factor @ cell-Center (-) */
1421 (PID.TID 0000.0001) 15 @ 1.000000000000000E+00 /* K = 1: 15 */
1422 (PID.TID 0000.0001) ;
1423 (PID.TID 0000.0001) deepFacF = /* deep-model grid factor @ W-Interface (-) */
1424 (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
1425 (PID.TID 0000.0001) ;
1426 (PID.TID 0000.0001) rVel2wUnit = /* convert units: rVel -> wSpeed (=1 if z-coord)*/
1427 (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
1428 (PID.TID 0000.0001) ;
1429 (PID.TID 0000.0001) wUnit2rVel = /* convert units: wSpeed -> rVel (=1 if z-coord)*/
1430 (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
1431 (PID.TID 0000.0001) ;
1432 (PID.TID 0000.0001) dBdrRef = /* Vertical grad. of reference buoyancy [(m/s/r)^2] */
1433 (PID.TID 0000.0001) 15 @ 0.000000000000000E+00 /* K = 1: 15 */
1434 (PID.TID 0000.0001) ;
1435 (PID.TID 0000.0001) rotateGrid = /* use rotated grid ( True/False ) */
1436 (PID.TID 0000.0001) F
1437 (PID.TID 0000.0001) ;
1438 (PID.TID 0000.0001) phiEuler = /* Euler angle, rotation about original z-coordinate [rad] */
1439 (PID.TID 0000.0001) 0.000000000000000E+00
1440 (PID.TID 0000.0001) ;
1441 (PID.TID 0000.0001) thetaEuler = /* Euler angle, rotation about new x-coordinate [rad] */
1442 (PID.TID 0000.0001) 0.000000000000000E+00
1443 (PID.TID 0000.0001) ;
1444 (PID.TID 0000.0001) psiEuler = /* Euler angle, rotation about new z-coordinate [rad] */
1445 (PID.TID 0000.0001) 0.000000000000000E+00
1446 (PID.TID 0000.0001) ;
1447 (PID.TID 0000.0001) dxF = /* dxF(:,1,:,1) ( units: m ) */
1448 (PID.TID 0000.0001) 90 @ 9.246038586187513E+04 /* I = 1: 90 */
1449 (PID.TID 0000.0001) ;
1450 (PID.TID 0000.0001) dxF = /* dxF(1,:,1,:) ( units: m ) */
1451 (PID.TID 0000.0001) 9.246038586187513E+04, /* J = 1 */
1452 (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 2 */
1453 (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 3 */
1454 (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 4 */
1455 (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 5 */
1456 (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 6 */
1457 (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 7 */
1458 (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 8 */
1459 (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 9 */
1460 (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 10 */
1461 (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 11 */
1462 (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 12 */
1463 (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 13 */
1464 (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 14 */
1465 (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 15 */
1466 (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 16 */
1467 (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 17 */
1468 (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 18 */
1469 (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 19 */
1470 (PID.TID 0000.0001) 2 @ 4.444389881567502E+05, /* J = 20: 21 */
1471 (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 22 */
1472 (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 23 */
1473 (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 24 */
1474 (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 25 */
1475 (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 26 */
1476 (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 27 */
1477 (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 28 */
1478 (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 29 */
1479 (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 30 */
1480 (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 31 */
1481 (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 32 */
1482 (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 33 */
1483 (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 34 */
1484 (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 35 */
1485 (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 36 */
1486 (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 37 */
1487 (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 38 */
1488 (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 39 */
1489 (PID.TID 0000.0001) 9.246038586187513E+04 /* J = 40 */
1490 (PID.TID 0000.0001) ;
1491 (PID.TID 0000.0001) dyF = /* dyF(:,1,:,1) ( units: m ) */
1492 (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1493 (PID.TID 0000.0001) ;
1494 (PID.TID 0000.0001) dyF = /* dyF(1,:,1,:) ( units: m ) */
1495 (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1496 (PID.TID 0000.0001) ;
1497 (PID.TID 0000.0001) dxG = /* dxG(:,1,:,1) ( units: m ) */
1498 (PID.TID 0000.0001) 90 @ 7.722306258078101E+04 /* I = 1: 90 */
1499 (PID.TID 0000.0001) ;
1500 (PID.TID 0000.0001) dxG = /* dxG(1,:,1,:) ( units: m ) */
1501 (PID.TID 0000.0001) 7.722306258078101E+04, /* J = 1 */
1502 (PID.TID 0000.0001) 1.075850604052271E+05, /* J = 2 */
1503 (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 3 */
1504 (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 4 */
1505 (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 5 */
1506 (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 6 */
1507 (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 7 */
1508 (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 8 */
1509 (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 9 */
1510 (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 10 */
1511 (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 11 */
1512 (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 12 */
1513 (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 13 */
1514 (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 14 */
1515 (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 15 */
1516 (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 16 */
1517 (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 17 */
1518 (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 18 */
1519 (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 19 */
1520 (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 20 */
1521 (PID.TID 0000.0001) 4.447098934081552E+05, /* J = 21 */
1522 (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 22 */
1523 (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 23 */
1524 (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 24 */
1525 (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 25 */
1526 (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 26 */
1527 (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 27 */
1528 (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 28 */
1529 (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 29 */
1530 (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 30 */
1531 (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 31 */
1532 (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 32 */
1533 (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 33 */
1534 (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 34 */
1535 (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 35 */
1536 (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 36 */
1537 (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 37 */
1538 (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 38 */
1539 (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 39 */
1540 (PID.TID 0000.0001) 1.075850604052271E+05 /* J = 40 */
1541 (PID.TID 0000.0001) ;
1542 (PID.TID 0000.0001) dyG = /* dyG(:,1,:,1) ( units: m ) */
1543 (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1544 (PID.TID 0000.0001) ;
1545 (PID.TID 0000.0001) dyG = /* dyG(1,:,1,:) ( units: m ) */
1546 (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1547 (PID.TID 0000.0001) ;
1548 (PID.TID 0000.0001) dxC = /* dxC(:,1,:,1) ( units: m ) */
1549 (PID.TID 0000.0001) 90 @ 9.246038586187513E+04 /* I = 1: 90 */
1550 (PID.TID 0000.0001) ;
1551 (PID.TID 0000.0001) dxC = /* dxC(1,:,1,:) ( units: m ) */
1552 (PID.TID 0000.0001) 9.246038586187513E+04, /* J = 1 */
1553 (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 2 */
1554 (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 3 */
1555 (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 4 */
1556 (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 5 */
1557 (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 6 */
1558 (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 7 */
1559 (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 8 */
1560 (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 9 */
1561 (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 10 */
1562 (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 11 */
1563 (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 12 */
1564 (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 13 */
1565 (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 14 */
1566 (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 15 */
1567 (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 16 */
1568 (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 17 */
1569 (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 18 */
1570 (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 19 */
1571 (PID.TID 0000.0001) 2 @ 4.444389881567502E+05, /* J = 20: 21 */
1572 (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 22 */
1573 (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 23 */
1574 (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 24 */
1575 (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 25 */
1576 (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 26 */
1577 (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 27 */
1578 (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 28 */
1579 (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 29 */
1580 (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 30 */
1581 (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 31 */
1582 (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 32 */
1583 (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 33 */
1584 (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 34 */
1585 (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 35 */
1586 (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 36 */
1587 (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 37 */
1588 (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 38 */
1589 (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 39 */
1590 (PID.TID 0000.0001) 9.246038586187513E+04 /* J = 40 */
1591 (PID.TID 0000.0001) ;
1592 (PID.TID 0000.0001) dyC = /* dyC(:,1,:,1) ( units: m ) */
1593 (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1594 (PID.TID 0000.0001) ;
1595 (PID.TID 0000.0001) dyC = /* dyC(1,:,1,:) ( units: m ) */
1596 (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1597 (PID.TID 0000.0001) ;
1598 (PID.TID 0000.0001) dxV = /* dxV(:,1,:,1) ( units: m ) */
1599 (PID.TID 0000.0001) 90 @ 7.722306258078101E+04 /* I = 1: 90 */
1600 (PID.TID 0000.0001) ;
1601 (PID.TID 0000.0001) dxV = /* dxV(1,:,1,:) ( units: m ) */
1602 (PID.TID 0000.0001) 7.722306258078101E+04, /* J = 1 */
1603 (PID.TID 0000.0001) 1.075850604052271E+05, /* J = 2 */
1604 (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 3 */
1605 (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 4 */
1606 (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 5 */
1607 (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 6 */
1608 (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 7 */
1609 (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 8 */
1610 (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 9 */
1611 (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 10 */
1612 (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 11 */
1613 (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 12 */
1614 (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 13 */
1615 (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 14 */
1616 (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 15 */
1617 (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 16 */
1618 (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 17 */
1619 (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 18 */
1620 (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 19 */
1621 (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 20 */
1622 (PID.TID 0000.0001) 4.447098934081552E+05, /* J = 21 */
1623 (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 22 */
1624 (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 23 */
1625 (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 24 */
1626 (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 25 */
1627 (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 26 */
1628 (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 27 */
1629 (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 28 */
1630 (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 29 */
1631 (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 30 */
1632 (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 31 */
1633 (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 32 */
1634 (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 33 */
1635 (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 34 */
1636 (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 35 */
1637 (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 36 */
1638 (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 37 */
1639 (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 38 */
1640 (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 39 */
1641 (PID.TID 0000.0001) 1.075850604052271E+05 /* J = 40 */
1642 (PID.TID 0000.0001) ;
1643 (PID.TID 0000.0001) dyU = /* dyU(:,1,:,1) ( units: m ) */
1644 (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1645 (PID.TID 0000.0001) ;
1646 (PID.TID 0000.0001) dyU = /* dyU(1,:,1,:) ( units: m ) */
1647 (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1648 (PID.TID 0000.0001) ;
1649 (PID.TID 0000.0001) rA = /* rA (:,1,:,1) ( units: m^2 ) */
1650 (PID.TID 0000.0001) 90 @ 4.110969866729047E+10 /* I = 1: 90 */
1651 (PID.TID 0000.0001) ;
1652 (PID.TID 0000.0001) rA = /* rA (1,:,1,:) ( units: m^2 ) */
1653 (PID.TID 0000.0001) 4.110969866729047E+10, /* J = 1 */
1654 (PID.TID 0000.0001) 5.450087291636666E+10, /* J = 2 */
1655 (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 3 */
1656 (PID.TID 0000.0001) 8.042270623659996E+10, /* J = 4 */
1657 (PID.TID 0000.0001) 9.282707674147525E+10, /* J = 5 */
1658 (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 6 */
1659 (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 7 */
1660 (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 8 */
1661 (PID.TID 0000.0001) 1.373525277677230E+11, /* J = 9 */
1662 (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 10 */
1663 (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 11 */
1664 (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 12 */
1665 (PID.TID 0000.0001) 1.712363708253574E+11, /* J = 13 */
1666 (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 14 */
1667 (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 15 */
1668 (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 16 */
1669 (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 17 */
1670 (PID.TID 0000.0001) 1.947228162702493E+11, /* J = 18 */
1671 (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 19 */
1672 (PID.TID 0000.0001) 2 @ 1.976062798008856E+11, /* J = 20: 21 */
1673 (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 22 */
1674 (PID.TID 0000.0001) 1.947228162702493E+11, /* J = 23 */
1675 (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 24 */
1676 (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 25 */
1677 (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 26 */
1678 (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 27 */
1679 (PID.TID 0000.0001) 1.712363708253574E+11, /* J = 28 */
1680 (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 29 */
1681 (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 30 */
1682 (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 31 */
1683 (PID.TID 0000.0001) 1.373525277677230E+11, /* J = 32 */
1684 (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 33 */
1685 (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 34 */
1686 (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 35 */
1687 (PID.TID 0000.0001) 9.282707674147525E+10, /* J = 36 */
1688 (PID.TID 0000.0001) 8.042270623659996E+10, /* J = 37 */
1689 (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 38 */
1690 (PID.TID 0000.0001) 5.450087291636666E+10, /* J = 39 */
1691 (PID.TID 0000.0001) 4.110969866729047E+10 /* J = 40 */
1692 (PID.TID 0000.0001) ;
1693 (PID.TID 0000.0001) rAw = /* rAw(:,1,:,1) ( units: m^2 ) */
1694 (PID.TID 0000.0001) 90 @ 4.110969866729047E+10 /* I = 1: 90 */
1695 (PID.TID 0000.0001) ;
1696 (PID.TID 0000.0001) rAw = /* rAw(1,:,1,:) ( units: m^2 ) */
1697 (PID.TID 0000.0001) 4.110969866729047E+10, /* J = 1 */
1698 (PID.TID 0000.0001) 5.450087291636666E+10, /* J = 2 */
1699 (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 3 */
1700 (PID.TID 0000.0001) 8.042270623659996E+10, /* J = 4 */
1701 (PID.TID 0000.0001) 9.282707674147525E+10, /* J = 5 */
1702 (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 6 */
1703 (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 7 */
1704 (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 8 */
1705 (PID.TID 0000.0001) 1.373525277677230E+11, /* J = 9 */
1706 (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 10 */
1707 (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 11 */
1708 (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 12 */
1709 (PID.TID 0000.0001) 1.712363708253574E+11, /* J = 13 */
1710 (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 14 */
1711 (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 15 */
1712 (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 16 */
1713 (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 17 */
1714 (PID.TID 0000.0001) 1.947228162702493E+11, /* J = 18 */
1715 (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 19 */
1716 (PID.TID 0000.0001) 2 @ 1.976062798008856E+11, /* J = 20: 21 */
1717 (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 22 */
1718 (PID.TID 0000.0001) 1.947228162702493E+11, /* J = 23 */
1719 (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 24 */
1720 (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 25 */
1721 (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 26 */
1722 (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 27 */
1723 (PID.TID 0000.0001) 1.712363708253574E+11, /* J = 28 */
1724 (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 29 */
1725 (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 30 */
1726 (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 31 */
1727 (PID.TID 0000.0001) 1.373525277677230E+11, /* J = 32 */
1728 (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 33 */
1729 (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 34 */
1730 (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 35 */
1731 (PID.TID 0000.0001) 9.282707674147525E+10, /* J = 36 */
1732 (PID.TID 0000.0001) 8.042270623659996E+10, /* J = 37 */
1733 (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 38 */
1734 (PID.TID 0000.0001) 5.450087291636666E+10, /* J = 39 */
1735 (PID.TID 0000.0001) 4.110969866729047E+10 /* J = 40 */
1736 (PID.TID 0000.0001) ;
1737 (PID.TID 0000.0001) rAs = /* rAs(:,1,:,1) ( units: m^2 ) */
1738 (PID.TID 0000.0001) 90 @ 3.433488626798251E+10 /* I = 1: 90 */
1739 (PID.TID 0000.0001) ;
1740 (PID.TID 0000.0001) rAs = /* rAs(1,:,1,:) ( units: m^2 ) */
1741 (PID.TID 0000.0001) 3.433488626798251E+10, /* J = 1 */
1742 (PID.TID 0000.0001) 4.783442523123625E+10, /* J = 2 */
1743 (PID.TID 0000.0001) 6.110091968306415E+10, /* J = 3 */
1744 (PID.TID 0000.0001) 7.406973659603818E+10, /* J = 4 */
1745 (PID.TID 0000.0001) 8.667769319778081E+10, /* J = 5 */
1746 (PID.TID 0000.0001) 9.886336479107465E+10, /* J = 6 */
1747 (PID.TID 0000.0001) 1.105673840088173E+11, /* J = 7 */
1748 (PID.TID 0000.0001) 1.217327300458638E+11, /* J = 8 */
1749 (PID.TID 0000.0001) 1.323050064586578E+11, /* J = 9 */
1750 (PID.TID 0000.0001) 1.422327061792377E+11, /* J = 10 */
1751 (PID.TID 0000.0001) 1.514674624524945E+11, /* J = 11 */
1752 (PID.TID 0000.0001) 1.599642844741385E+11, /* J = 12 */
1753 (PID.TID 0000.0001) 1.676817765813788E+11, /* J = 13 */
1754 (PID.TID 0000.0001) 1.745823399284268E+11, /* J = 14 */
1755 (PID.TID 0000.0001) 1.806323556642997E+11, /* J = 15 */
1756 (PID.TID 0000.0001) 1.858023487204767E+11, /* J = 16 */
1757 (PID.TID 0000.0001) 1.900671314104744E+11, /* J = 17 */
1758 (PID.TID 0000.0001) 1.934059261417216E+11, /* J = 18 */
1759 (PID.TID 0000.0001) 1.958024666419019E+11, /* J = 19 */
1760 (PID.TID 0000.0001) 1.972450772065981E+11, /* J = 20 */
1761 (PID.TID 0000.0001) 1.977267295821496E+11, /* J = 21 */
1762 (PID.TID 0000.0001) 1.972450772065981E+11, /* J = 22 */
1763 (PID.TID 0000.0001) 1.958024666419019E+11, /* J = 23 */
1764 (PID.TID 0000.0001) 1.934059261417216E+11, /* J = 24 */
1765 (PID.TID 0000.0001) 1.900671314104744E+11, /* J = 25 */
1766 (PID.TID 0000.0001) 1.858023487204767E+11, /* J = 26 */
1767 (PID.TID 0000.0001) 1.806323556642997E+11, /* J = 27 */
1768 (PID.TID 0000.0001) 1.745823399284268E+11, /* J = 28 */
1769 (PID.TID 0000.0001) 1.676817765813788E+11, /* J = 29 */
1770 (PID.TID 0000.0001) 1.599642844741385E+11, /* J = 30 */
1771 (PID.TID 0000.0001) 1.514674624524945E+11, /* J = 31 */
1772 (PID.TID 0000.0001) 1.422327061792377E+11, /* J = 32 */
1773 (PID.TID 0000.0001) 1.323050064586578E+11, /* J = 33 */
1774 (PID.TID 0000.0001) 1.217327300458638E+11, /* J = 34 */
1775 (PID.TID 0000.0001) 1.105673840088173E+11, /* J = 35 */
1776 (PID.TID 0000.0001) 9.886336479107465E+10, /* J = 36 */
1777 (PID.TID 0000.0001) 8.667769319778081E+10, /* J = 37 */
1778 (PID.TID 0000.0001) 7.406973659603818E+10, /* J = 38 */
1779 (PID.TID 0000.0001) 6.110091968306415E+10, /* J = 39 */
1780 (PID.TID 0000.0001) 4.783442523123625E+10 /* J = 40 */
1781 (PID.TID 0000.0001) ;
1782 (PID.TID 0000.0001) globalArea = /* Integrated horizontal Area (m^2) */
1783 (PID.TID 0000.0001) 3.450614146649749E+14
1784 (PID.TID 0000.0001) ;
1785 (PID.TID 0000.0001) // =======================================================
1786 (PID.TID 0000.0001) // End of Model config. summary
1787 (PID.TID 0000.0001) // =======================================================
1788 (PID.TID 0000.0001)
1789 (PID.TID 0000.0001) == Packages configuration : Check & print summary ==
1790 (PID.TID 0000.0001)
1791 (PID.TID 0000.0001) GMREDI_CHECK: #define GMREDI
1792 (PID.TID 0000.0001) GM_AdvForm = /* if FALSE => use SkewFlux Form */
1793 (PID.TID 0000.0001) T
1794 (PID.TID 0000.0001) ;
1795 (PID.TID 0000.0001) GM_InMomAsStress = /* if TRUE => apply as Eddy Stress */
1796 (PID.TID 0000.0001) F
1797 OAD: TIMING: stamp 0: 1475718770.202285
1798 OAD: TIMING: stamp 1: 1475718790.930249
1799 OAD: TIMING: delta stamps 1-0: 20.727964
1800 OAD: TIMING: stamp 2: 1475718837.828353
1801 OAD: TIMING: delta stamps 2-1: 46.898104
1802 OAD: TIMING: delta stamps 2-0: 67.626068
1803 OAD: TIMING: ratio stamps (2-1)/(1-0): 4.689810e+07/2.072796e+07=2.262552e+00
1804 (PID.TID 0000.0001) ;
1805 (PID.TID 0000.0001) GM_AdvSeparate = /* Calc Bolus & Euler Adv. separately */
1806 (PID.TID 0000.0001) F
1807 (PID.TID 0000.0001) ;
1808 (PID.TID 0000.0001) GM_ExtraDiag = /* Tensor Extra Diag (line 1&2) non 0 */
1809 (PID.TID 0000.0001) T
1810 (PID.TID 0000.0001) ;
1811 (PID.TID 0000.0001) GM_isopycK = /* Background Isopyc. Diffusivity [m^2/s] */
1812 (PID.TID 0000.0001) 1.100000000000000E+03
1813 (PID.TID 0000.0001) ;
1814 (PID.TID 0000.0001) GM_skewflx*K = /* Background GM_SkewFlx Diffusivity [m^2/s] */
1815 (PID.TID 0000.0001) 0.000000000000000E+00
1816 (PID.TID 0000.0001) ;
1817 (PID.TID 0000.0001) GM_advec*K = /* Backg. GM-Advec(=Bolus) Diffusivity [m^2/s]*/
1818 (PID.TID 0000.0001) 9.000000000000000E+02
1819 (PID.TID 0000.0001) ;
1820 (PID.TID 0000.0001) GM_Kmin_horiz = /* Minimum Horizontal Diffusivity [m^2/s] */
1821 (PID.TID 0000.0001) 5.000000000000000E+01
1822 (PID.TID 0000.0001) ;
1823 (PID.TID 0000.0001) GM_Visbeck_alpha = /* Visbeck alpha coeff. [-] */
1824 (PID.TID 0000.0001) 0.000000000000000E+00
1825 (PID.TID 0000.0001) ;
1826 (PID.TID 0000.0001) GM_Small_Number = /* epsilon used in slope calc */
1827 (PID.TID 0000.0001) 9.999999999999999E-21
1828 (PID.TID 0000.0001) ;
1829 (PID.TID 0000.0001) GM_slopeSqCutoff = /* Slope^2 cut-off value */
1830 (PID.TID 0000.0001) 1.000000000000000E+08
1831 (PID.TID 0000.0001) ;
1832 (PID.TID 0000.0001) GM_taper_scheme = /* Type of Tapering/Clipping scheme */
1833 (PID.TID 0000.0001) 'dm95 '
1834 (PID.TID 0000.0001) ;
1835 (PID.TID 0000.0001) GM_maxSlope = /* Maximum Slope (Tapering/Clipping) */
1836 (PID.TID 0000.0001) 1.000000000000000E-02
1837 (PID.TID 0000.0001) ;
1838 (PID.TID 0000.0001) GM_facTrL2dz = /* Minimum Trans.Layer Thick. (factor of dz) */
1839 (PID.TID 0000.0001) 1.000000000000000E+00
1840 (PID.TID 0000.0001) ;
1841 (PID.TID 0000.0001) GM_facTrL2ML = /* Max.Trans.Layer Thick. (factor of MxL Depth)*/
1842 (PID.TID 0000.0001) 5.000000000000000E+00
1843 (PID.TID 0000.0001) ;
1844 (PID.TID 0000.0001) GM_maxTransLay = /* Maximum Transition Layer Thickness [m] */
1845 (PID.TID 0000.0001) 5.000000000000000E+02
1846 (PID.TID 0000.0001) ;
1847 (PID.TID 0000.0001) GM_UseBVP = /* if TRUE => use bvp a la Ferrari et al. (2010) */
1848 (PID.TID 0000.0001) F
1849 (PID.TID 0000.0001) ;
1850 (PID.TID 0000.0001) GM_BVP_ModeNumber = /* Vertical mode number for BVP wave speed */
1851 (PID.TID 0000.0001) 1
1852 (PID.TID 0000.0001) ;
1853 (PID.TID 0000.0001) GM_BVP_cMin = /* Minimum wave speed for BVP [m/s] */
1854 (PID.TID 0000.0001) 1.000000000000000E-01
1855 (PID.TID 0000.0001) ;
1856 (PID.TID 0000.0001) GM_useSubMeso = /* if TRUE => use Sub-Meso param. (B.Fox-Kemper) */
1857 (PID.TID 0000.0001) F
1858 (PID.TID 0000.0001) ;
1859 (PID.TID 0000.0001) subMeso_Ceff = /* efficiency coeff. of Mixed-Layer Eddies [-] */
1860 (PID.TID 0000.0001) 7.000000000000001E-02
1861 (PID.TID 0000.0001) ;
1862 (PID.TID 0000.0001) subMeso_invTau = /* inverse of Sub-Meso mixing time-scale [/s] */
1863 (PID.TID 0000.0001) 2.000000000000000E-06
1864 (PID.TID 0000.0001) ;
1865 (PID.TID 0000.0001) subMeso_LfMin = /* minimum length-scale "Lf" [m] */
1866 (PID.TID 0000.0001) 1.000000000000000E+03
1867 (PID.TID 0000.0001) ;
1868 (PID.TID 0000.0001) subMeso_Lmax = /* maximum grid-scale length [m] */
1869 (PID.TID 0000.0001) 1.100000000000000E+05
1870 (PID.TID 0000.0001) ;
1871 (PID.TID 0000.0001) CTRL_CHECK: ctrl package
1872 (PID.TID 0000.0001) COST_CHECK: cost package
1873 (PID.TID 0000.0001) GRDCHK_CHECK: grdchk package
1874 (PID.TID 0000.0001) GAD_CHECK: #define ALLOW_GENERIC_ADVDIFF
1875 (PID.TID 0000.0001) // =======================================================
1876 (PID.TID 0000.0001) // Check Model config. (CONFIG_CHECK):
1877 (PID.TID 0000.0001) // CONFIG_CHECK : Normal End
1878 (PID.TID 0000.0001) // =======================================================
1879 (PID.TID 0000.0001)
1880 OAD: IT+ 1048576
1881 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
1882 (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
1883 (PID.TID 0000.0001)
1884 (PID.TID 0000.0001) // =======================================================
1885 (PID.TID 0000.0001) // Model current state
1886 (PID.TID 0000.0001) // =======================================================
1887 (PID.TID 0000.0001)
1888 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
1889 (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
1890 cg2d: Sum(rhs),rhsMax = -2.77555756156289E-15 2.53681867617525E+00
1891 cg2d: Sum(rhs),rhsMax = 5.32907051820075E-15 4.93238674419718E+00
1892 cg2d: Sum(rhs),rhsMax = 2.22044604925031E-16 6.28229902439776E+00
1893 cg2d: Sum(rhs),rhsMax = -2.04281036531029E-14 6.78919330551124E+00
1894 cg2d: Sum(rhs),rhsMax = -1.06581410364015E-14 6.88673552302266E+00
1895 cg2d: Sum(rhs),rhsMax = 1.15463194561016E-14 6.84253554757674E+00
1896 cg2d: Sum(rhs),rhsMax = -1.77635683940025E-14 7.17243286106845E+00
1897 cg2d: Sum(rhs),rhsMax = 6.21724893790088E-14 7.77283039421207E+00
1898 cg2d: Sum(rhs),rhsMax = 4.26325641456060E-14 8.10136142073580E+00
1899 OAD: DT+ 1048576
1900 OAD: IT+ 17825792
1901 OAD: IT+ 34603008
1902 OAD: DT+ 17825792
1903 OAD: IT+ 51380224
1904 OAD: IT+ 68157440
1905 cg2d: Sum(rhs),rhsMax = 3.19744231092045E-14 8.18927452945244E+00
1906 OAD: DT+ 34603008
1907 OAD: IT+ 84934656
1908 OAD: IT+ 101711872
1909 (PID.TID 0000.0001) %CHECKPOINT 10 ckptA
1910 early fc = 0.000000000000000D+00
1911 --> objf_atl(bi,bj) = 0.000000000000000D+00
1912 --> objf_atl(bi,bj) = 0.000000000000000D+00
1913 --> objf_atl(bi,bj) = 0.000000000000000D+00
1914 --> objf_atl(bi,bj) = -0.450065077785263D-01
1915 local fc = -0.450065077785263D-01
1916 global fc = -0.450065077785263D-01
1917 cg2d: Sum(rhs),rhsMax = 1.23858795651938E-15 8.04945967787813E-09
1918 cg2d: Sum(rhs),rhsMax = -2.66453525910038E-14 8.10136142073580E+00
1919 cg2d: Sum(rhs),rhsMax = 3.52791255410490E-15 9.21369200784079E-09
1920 cg2d: Sum(rhs),rhsMax = 3.19744231092045E-14 7.77283039421207E+00
1921 cg2d: Sum(rhs),rhsMax = 6.03553665379231E-15 8.74327613604373E-09
1922 cg2d: Sum(rhs),rhsMax = -2.66453525910038E-14 7.17243286106845E+00
1923 cg2d: Sum(rhs),rhsMax = 3.42260941810224E-15 5.49414350830240E-09
1924 cg2d: Sum(rhs),rhsMax = 6.03961325396085E-14 6.84253554757673E+00
1925 cg2d: Sum(rhs),rhsMax = 4.62824223390612E-15 3.68281353912157E-09
1926 cg2d: Sum(rhs),rhsMax = 3.64153152077051E-14 6.88673552302266E+00
1927 cg2d: Sum(rhs),rhsMax = -8.77076189453874E-15 5.90272604043300E-09
1928 cg2d: Sum(rhs),rhsMax = -8.43769498715119E-15 6.78919330551122E+00
1929 cg2d: Sum(rhs),rhsMax = -2.35922392732846E-15 7.26634540912282E-09
1930 cg2d: Sum(rhs),rhsMax = -6.77236045021345E-15 6.28229902439776E+00
1931 cg2d: Sum(rhs),rhsMax = -1.04916075827077E-14 8.24007965726769E-09
1932 cg2d: Sum(rhs),rhsMax = 5.71764857681956E-15 4.93238674419719E+00
1933 cg2d: Sum(rhs),rhsMax = -6.77236045021345E-15 9.24289205824485E-09
1934 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
1935 (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
1936 cg2d: Sum(rhs),rhsMax = -4.44089209850063E-16 2.53681867617526E+00
1937 cg2d: Sum(rhs),rhsMax = 1.44328993201270E-15 1.10566559735434E-08
1938 ph-pack: packing ecco_cost
1939 ph-pack: packing ecco_ctrl
1940 (PID.TID 0000.0001) // =======================================================
1941 (PID.TID 0000.0001) // Gradient-check starts (grdchk_main)
1942 (PID.TID 0000.0001) // =======================================================
1943 (PID.TID 0000.0001) grdchk reference fc: fcref = -4.50065077785263E-02
1944 grad-res -------------------------------
1945 grad-res proc # i j k bi bj iobc fc ref fc + eps fc - eps
1946 grad-res proc # i j k bi bj iobc adj grad fd grad 1 - fd/adj
1947 grad-res exact position met:
1948 grad-res 0 258 31 7 1 2 2
1949 (PID.TID 0000.0001) ====== Starts gradient-check number 1 (=ichknum) =======
1950 ph-test icomp, ncvarcomp, ichknum 258 29309 1
1951 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 0 1
1952 ph-grd -->hit<-- 31 7 1 1
1953 (PID.TID 0000.0001) grdchk pos: i,j,k= 31 7 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1
1954 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
1955 (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
1956 (PID.TID 0000.0001)
1957 (PID.TID 0000.0001) // =======================================================
1958 (PID.TID 0000.0001) // Model current state
1959 (PID.TID 0000.0001) // =======================================================
1960 (PID.TID 0000.0001)
1961 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
1962 (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
1963 cg2d: Sum(rhs),rhsMax = -3.66373598126302E-15 2.53681867617525E+00
1964 cg2d: Sum(rhs),rhsMax = -4.66293670342566E-15 4.93238674419718E+00
1965 cg2d: Sum(rhs),rhsMax = -1.90958360235527E-14 6.28229902439776E+00
1966 cg2d: Sum(rhs),rhsMax = -4.88498130835069E-15 6.78919330551119E+00
1967 cg2d: Sum(rhs),rhsMax = 2.30926389122033E-14 6.88673552302264E+00
1968 cg2d: Sum(rhs),rhsMax = 6.30606677987089E-14 6.84253554757676E+00
1969 cg2d: Sum(rhs),rhsMax = -7.46069872548105E-14 7.17243286106857E+00
1970 cg2d: Sum(rhs),rhsMax = 8.70414851306123E-14 7.77283039421237E+00
1971 cg2d: Sum(rhs),rhsMax = -5.32907051820075E-15 8.10136142073602E+00
1972 cg2d: Sum(rhs),rhsMax = 3.01980662698043E-14 8.18927452945243E+00
1973 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
1974 early fc = 0.000000000000000D+00
1975 --> objf_atl(bi,bj) = 0.000000000000000D+00
1976 --> objf_atl(bi,bj) = 0.000000000000000D+00
1977 --> objf_atl(bi,bj) = 0.000000000000000D+00
1978 --> objf_atl(bi,bj) = -0.450066523719723D-01
1979 local fc = -0.450066523719723D-01
1980 global fc = -0.450066523719723D-01
1981 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = -4.50066523719723E-02
1982 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
1983 (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
1984 (PID.TID 0000.0001)
1985 (PID.TID 0000.0001) // =======================================================
1986 (PID.TID 0000.0001) // Model current state
1987 (PID.TID 0000.0001) // =======================================================
1988 (PID.TID 0000.0001)
1989 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
1990 (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
1991 cg2d: Sum(rhs),rhsMax = -3.66373598126302E-15 2.53681867617525E+00
1992 cg2d: Sum(rhs),rhsMax = -2.05391259555654E-15 4.93238674419718E+00
1993 cg2d: Sum(rhs),rhsMax = 1.99840144432528E-15 6.28229902439777E+00
1994 cg2d: Sum(rhs),rhsMax = 3.10862446895044E-15 6.78919330551119E+00
1995 cg2d: Sum(rhs),rhsMax = 1.42108547152020E-14 6.88673552302310E+00
1996 cg2d: Sum(rhs),rhsMax = 1.42108547152020E-14 6.84253554757712E+00
1997 cg2d: Sum(rhs),rhsMax = 7.28306304154103E-14 7.17243286106881E+00
1998 cg2d: Sum(rhs),rhsMax = 7.10542735760100E-14 7.77283039421259E+00
1999 cg2d: Sum(rhs),rhsMax = 6.03961325396085E-14 8.10136142073610E+00
2000 cg2d: Sum(rhs),rhsMax = 3.55271367880050E-14 8.18927452945214E+00
2001 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2002 early fc = 0.000000000000000D+00
2003 --> objf_atl(bi,bj) = 0.000000000000000D+00
2004 --> objf_atl(bi,bj) = 0.000000000000000D+00
2005 --> objf_atl(bi,bj) = 0.000000000000000D+00
2006 --> objf_atl(bi,bj) = -0.450063631388652D-01
2007 local fc = -0.450063631388652D-01
2008 global fc = -0.450063631388652D-01
2009 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = -4.50063631388652E-02
2010 grad-res -------------------------------
2011 grad-res 0 1 31 7 1 2 2 1 -4.50065077785E-02 -4.50066523720E-02 -4.50063631389E-02
2012 grad-res 0 1 1 258 0 2 2 1 -1.44600218759E-05 -1.44616553583E-05 -1.12965418225E-04
2013 (PID.TID 0000.0001) ADM ref_cost_function = -4.50065077785263E-02
2014 (PID.TID 0000.0001) ADM adjoint_gradient = -1.44600218758839E-05
2015 (PID.TID 0000.0001) ADM finite-diff_grad = -1.44616553583027E-05
2016 (PID.TID 0000.0001) ====== End of gradient-check number 1 (ierr= 0) =======
2017 (PID.TID 0000.0001) ====== Starts gradient-check number 2 (=ichknum) =======
2018 ph-test icomp, ncvarcomp, ichknum 259 29309 2
2019 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 258 2
2020 ph-grd -->hit<-- 32 7 1 1
2021 (PID.TID 0000.0001) grdchk pos: i,j,k= 32 7 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1
2022 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2023 (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2024 (PID.TID 0000.0001)
2025 (PID.TID 0000.0001) // =======================================================
2026 (PID.TID 0000.0001) // Model current state
2027 (PID.TID 0000.0001) // =======================================================
2028 (PID.TID 0000.0001)
2029 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2030 (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2031 cg2d: Sum(rhs),rhsMax = -2.77555756156289E-15 2.53681867617525E+00
2032 cg2d: Sum(rhs),rhsMax = -4.21884749357559E-15 4.93238674419718E+00
2033 cg2d: Sum(rhs),rhsMax = -2.22044604925031E-16 6.28229902439776E+00
2034 cg2d: Sum(rhs),rhsMax = 3.64153152077051E-14 6.78919330551124E+00
2035 cg2d: Sum(rhs),rhsMax = -1.06581410364015E-14 6.88673552302285E+00
2036 cg2d: Sum(rhs),rhsMax = 7.10542735760100E-15 6.84253554757690E+00
2037 cg2d: Sum(rhs),rhsMax = -6.92779167366098E-14 7.17243286106887E+00
2038 cg2d: Sum(rhs),rhsMax = 6.03961325396085E-14 7.77283039421257E+00
2039 cg2d: Sum(rhs),rhsMax = -3.19744231092045E-14 8.10136142073614E+00
2040 cg2d: Sum(rhs),rhsMax = -3.37507799486048E-14 8.18927452945254E+00
2041 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2042 early fc = 0.000000000000000D+00
2043 --> objf_atl(bi,bj) = 0.000000000000000D+00
2044 --> objf_atl(bi,bj) = 0.000000000000000D+00
2045 --> objf_atl(bi,bj) = 0.000000000000000D+00
2046 --> objf_atl(bi,bj) = -0.450066425813319D-01
2047 local fc = -0.450066425813319D-01
2048 global fc = -0.450066425813319D-01
2049 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = -4.50066425813319E-02
2050 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2051 (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2052 (PID.TID 0000.0001)
2053 (PID.TID 0000.0001) // =======================================================
2054 (PID.TID 0000.0001) // Model current state
2055 (PID.TID 0000.0001) // =======================================================
2056 (PID.TID 0000.0001)
2057 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2058 (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2059 cg2d: Sum(rhs),rhsMax = -3.66373598126302E-15 2.53681867617525E+00
2060 cg2d: Sum(rhs),rhsMax = 7.60502771868232E-15 4.93238674419718E+00
2061 cg2d: Sum(rhs),rhsMax = -1.31006316905768E-14 6.28229902439776E+00
2062 cg2d: Sum(rhs),rhsMax = -5.32907051820075E-15 6.78919330551120E+00
2063 cg2d: Sum(rhs),rhsMax = 3.73034936274053E-14 6.88673552302311E+00
2064 cg2d: Sum(rhs),rhsMax = 6.57252030578093E-14 6.84253554757708E+00
2065 cg2d: Sum(rhs),rhsMax = -6.57252030578093E-14 7.17243286106846E+00
2066 cg2d: Sum(rhs),rhsMax = 0.00000000000000E+00 7.77283039421238E+00
2067 cg2d: Sum(rhs),rhsMax = 1.06581410364015E-13 8.10136142073615E+00
2068 cg2d: Sum(rhs),rhsMax = -6.21724893790088E-14 8.18927452945265E+00
2069 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2070 early fc = 0.000000000000000D+00
2071 --> objf_atl(bi,bj) = 0.000000000000000D+00
2072 --> objf_atl(bi,bj) = 0.000000000000000D+00
2073 --> objf_atl(bi,bj) = 0.000000000000000D+00
2074 --> objf_atl(bi,bj) = -0.450063729255709D-01
2075 local fc = -0.450063729255709D-01
2076 global fc = -0.450063729255709D-01
2077 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = -4.50063729255709E-02
2078 grad-res -------------------------------
2079 grad-res 0 2 32 7 1 2 2 1 -4.50065077785E-02 -4.50066425813E-02 -4.50063729256E-02
2080 grad-res 0 2 2 259 0 2 2 1 -1.34827845971E-05 -1.34827880468E-05 -2.55859901044E-07
2081 (PID.TID 0000.0001) ADM ref_cost_function = -4.50065077785263E-02
2082 (PID.TID 0000.0001) ADM adjoint_gradient = -1.34827845970858E-05
2083 (PID.TID 0000.0001) ADM finite-diff_grad = -1.34827880467897E-05
2084 (PID.TID 0000.0001) ====== End of gradient-check number 2 (ierr= 0) =======
2085 (PID.TID 0000.0001) ====== Starts gradient-check number 3 (=ichknum) =======
2086 ph-test icomp, ncvarcomp, ichknum 260 29309 3
2087 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 259 3
2088 ph-grd -->hit<-- 33 7 1 1
2089 (PID.TID 0000.0001) grdchk pos: i,j,k= 33 7 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1
2090 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2091 (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2092 (PID.TID 0000.0001)
2093 (PID.TID 0000.0001) // =======================================================
2094 (PID.TID 0000.0001) // Model current state
2095 (PID.TID 0000.0001) // =======================================================
2096 (PID.TID 0000.0001)
2097 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2098 (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2099 cg2d: Sum(rhs),rhsMax = -2.77555756156289E-15 2.53681867617525E+00
2100 cg2d: Sum(rhs),rhsMax = 3.10862446895044E-15 4.93238674419718E+00
2101 cg2d: Sum(rhs),rhsMax = 6.66133814775094E-15 6.28229902439776E+00
2102 cg2d: Sum(rhs),rhsMax = -2.48689957516035E-14 6.78919330551120E+00
2103 cg2d: Sum(rhs),rhsMax = -2.48689957516035E-14 6.88673552302264E+00
2104 cg2d: Sum(rhs),rhsMax = 3.37507799486048E-14 6.84253554757685E+00
2105 cg2d: Sum(rhs),rhsMax = 2.48689957516035E-14 7.17243286106846E+00
2106 cg2d: Sum(rhs),rhsMax = -7.28306304154103E-14 7.77283039421234E+00
2107 cg2d: Sum(rhs),rhsMax = 5.50670620214078E-14 8.10136142073606E+00
2108 cg2d: Sum(rhs),rhsMax = 4.97379915032070E-14 8.18927452945247E+00
2109 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2110 early fc = 0.000000000000000D+00
2111 --> objf_atl(bi,bj) = 0.000000000000000D+00
2112 --> objf_atl(bi,bj) = 0.000000000000000D+00
2113 --> objf_atl(bi,bj) = 0.000000000000000D+00
2114 --> objf_atl(bi,bj) = -0.450066011319008D-01
2115 local fc = -0.450066011319008D-01
2116 global fc = -0.450066011319008D-01
2117 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = -4.50066011319008E-02
2118 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2119 (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2120 (PID.TID 0000.0001)
2121 (PID.TID 0000.0001) // =======================================================
2122 (PID.TID 0000.0001) // Model current state
2123 (PID.TID 0000.0001) // =======================================================
2124 (PID.TID 0000.0001)
2125 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2126 (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2127 cg2d: Sum(rhs),rhsMax = -3.66373598126302E-15 2.53681867617525E+00
2128 cg2d: Sum(rhs),rhsMax = -7.71605002114484E-15 4.93238674419718E+00
2129 cg2d: Sum(rhs),rhsMax = 1.36557432028894E-14 6.28229902439776E+00
2130 cg2d: Sum(rhs),rhsMax = -1.28785870856518E-14 6.78919330551120E+00
2131 cg2d: Sum(rhs),rhsMax = 1.86517468137026E-14 6.88673552302264E+00
2132 cg2d: Sum(rhs),rhsMax = 1.77635683940025E-15 6.84253554757676E+00
2133 cg2d: Sum(rhs),rhsMax = 4.97379915032070E-14 7.17243286106837E+00
2134 cg2d: Sum(rhs),rhsMax = -2.48689957516035E-14 7.77283039421237E+00
2135 cg2d: Sum(rhs),rhsMax = -6.57252030578093E-14 8.10136142073610E+00
2136 cg2d: Sum(rhs),rhsMax = -1.17239551400417E-13 8.18927452945241E+00
2137 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2138 early fc = 0.000000000000000D+00
2139 --> objf_atl(bi,bj) = 0.000000000000000D+00
2140 --> objf_atl(bi,bj) = 0.000000000000000D+00
2141 --> objf_atl(bi,bj) = 0.000000000000000D+00
2142 --> objf_atl(bi,bj) = -0.450064147383346D-01
2143 local fc = -0.450064147383346D-01
2144 global fc = -0.450064147383346D-01
2145 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = -4.50064147383346E-02
2146 grad-res -------------------------------
2147 grad-res 0 3 33 7 1 2 2 1 -4.50065077785E-02 -4.50066011319E-02 -4.50064147383E-02
2148 grad-res 0 3 3 260 0 2 2 1 -9.32005499887E-06 -9.31967830894E-06 4.04171356050E-05
2149 (PID.TID 0000.0001) ADM ref_cost_function = -4.50065077785263E-02
2150 (PID.TID 0000.0001) ADM adjoint_gradient = -9.32005499887116E-06
2151 (PID.TID 0000.0001) ADM finite-diff_grad = -9.31967830894442E-06
2152 (PID.TID 0000.0001) ====== End of gradient-check number 3 (ierr= 0) =======
2153 (PID.TID 0000.0001) ====== Starts gradient-check number 4 (=ichknum) =======
2154 ph-test icomp, ncvarcomp, ichknum 261 29309 4
2155 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 260 4
2156 ph-grd -->hit<-- 34 7 1 1
2157 (PID.TID 0000.0001) grdchk pos: i,j,k= 34 7 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1
2158 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2159 (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2160 (PID.TID 0000.0001)
2161 (PID.TID 0000.0001) // =======================================================
2162 (PID.TID 0000.0001) // Model current state
2163 (PID.TID 0000.0001) // =======================================================
2164 (PID.TID 0000.0001)
2165 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2166 (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2167 cg2d: Sum(rhs),rhsMax = -2.77555756156289E-15 2.53681867617525E+00
2168 cg2d: Sum(rhs),rhsMax = -5.49560397189452E-15 4.93238674419718E+00
2169 cg2d: Sum(rhs),rhsMax = 4.44089209850063E-16 6.28229902439777E+00
2170 cg2d: Sum(rhs),rhsMax = 1.33226762955019E-15 6.78919330551120E+00
2171 cg2d: Sum(rhs),rhsMax = -4.70734562441066E-14 6.88673552302311E+00
2172 cg2d: Sum(rhs),rhsMax = -3.46389583683049E-14 6.84253554757713E+00
2173 cg2d: Sum(rhs),rhsMax = 1.06581410364015E-14 7.17243286106843E+00
2174 cg2d: Sum(rhs),rhsMax = 6.92779167366098E-14 7.77283039421208E+00
2175 cg2d: Sum(rhs),rhsMax = -1.77635683940025E-15 8.10136142073581E+00
2176 cg2d: Sum(rhs),rhsMax = -7.46069872548105E-14 8.18927452945245E+00
2177 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2178 early fc = 0.000000000000000D+00
2179 --> objf_atl(bi,bj) = 0.000000000000000D+00
2180 --> objf_atl(bi,bj) = 0.000000000000000D+00
2181 --> objf_atl(bi,bj) = 0.000000000000000D+00
2182 --> objf_atl(bi,bj) = -0.450066465031474D-01
2183 local fc = -0.450066465031474D-01
2184 global fc = -0.450066465031474D-01
2185 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = -4.50066465031474E-02
2186 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2187 (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2188 (PID.TID 0000.0001)
2189 (PID.TID 0000.0001) // =======================================================
2190 (PID.TID 0000.0001) // Model current state
2191 (PID.TID 0000.0001) // =======================================================
2192 (PID.TID 0000.0001)
2193 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2194 (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2195 cg2d: Sum(rhs),rhsMax = -3.66373598126302E-15 2.53681867617525E+00
2196 cg2d: Sum(rhs),rhsMax = 3.88578058618805E-15 4.93238674419718E+00
2197 cg2d: Sum(rhs),rhsMax = 1.11022302462516E-15 6.28229902439777E+00
2198 cg2d: Sum(rhs),rhsMax = -1.33226762955019E-14 6.78919330551120E+00
2199 cg2d: Sum(rhs),rhsMax = -3.73034936274053E-14 6.88673552302311E+00
2200 cg2d: Sum(rhs),rhsMax = -3.55271367880050E-15 6.84253554757709E+00
2201 cg2d: Sum(rhs),rhsMax = -3.73034936274053E-14 7.17243286106841E+00
2202 cg2d: Sum(rhs),rhsMax = 4.44089209850063E-14 7.77283039421204E+00
2203 cg2d: Sum(rhs),rhsMax = -4.97379915032070E-14 8.10136142073584E+00
2204 cg2d: Sum(rhs),rhsMax = 1.01252339845814E-13 8.18927452945255E+00
2205 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2206 early fc = 0.000000000000000D+00
2207 --> objf_atl(bi,bj) = 0.000000000000000D+00
2208 --> objf_atl(bi,bj) = 0.000000000000000D+00
2209 --> objf_atl(bi,bj) = 0.000000000000000D+00
2210 --> objf_atl(bi,bj) = -0.450063523401348D-01
2211 local fc = -0.450063523401348D-01
2212 global fc = -0.450063523401348D-01
2213 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = -4.50063523401348E-02
2214 grad-res -------------------------------
2215 grad-res 0 4 34 7 1 2 2 1 -4.50065077785E-02 -4.50066465031E-02 -4.50063523401E-02
2216 grad-res 0 4 4 261 0 2 2 1 -1.66843992494E-05 -1.47081506297E-05 1.18448892893E-01
2217 (PID.TID 0000.0001) ADM ref_cost_function = -4.50065077785263E-02
2218 (PID.TID 0000.0001) ADM adjoint_gradient = -1.66843992494025E-05
2219 (PID.TID 0000.0001) ADM finite-diff_grad = -1.47081506297192E-05
2220 (PID.TID 0000.0001) ====== End of gradient-check number 4 (ierr= 0) =======
2221 (PID.TID 0000.0001)
2222 (PID.TID 0000.0001) // =======================================================
2223 (PID.TID 0000.0001) // Gradient check results >>> START <<<
2224 (PID.TID 0000.0001) // =======================================================
2225 (PID.TID 0000.0001)
2226 (PID.TID 0000.0001) EPS = 1.000000E-02
2227 (PID.TID 0000.0001)
2228 (PID.TID 0000.0001) grdchk output h.p: Id Itile Jtile LAYER bi bj X(Id) X(Id)+/-EPS
2229 (PID.TID 0000.0001) grdchk output h.c: Id FC FC1 FC2
2230 (PID.TID 0000.0001) grdchk output h.g: Id FC1-FC2/(2*EPS) ADJ GRAD(FC) 1-FDGRD/ADGRD
2231 (PID.TID 0000.0001)
2232 (PID.TID 0000.0001) grdchk output (p): 1 31 7 1 2 2 0.000000000E+00 -1.000000000E-02
2233 (PID.TID 0000.0001) grdchk output (c): 1 -4.5006507778526E-02 -4.5006652371972E-02 -4.5006363138865E-02
2234 (PID.TID 0000.0001) grdchk output (g): 1 -1.4461655358303E-05 -1.4460021875884E-05 -1.1296541822547E-04
2235 (PID.TID 0000.0001)
2236 (PID.TID 0000.0001) grdchk output (p): 2 32 7 1 2 2 0.000000000E+00 -1.000000000E-02
2237 (PID.TID 0000.0001) grdchk output (c): 2 -4.5006507778526E-02 -4.5006642581332E-02 -4.5006372925571E-02
2238 (PID.TID 0000.0001) grdchk output (g): 2 -1.3482788046790E-05 -1.3482784597086E-05 -2.5585990104382E-07
2239 (PID.TID 0000.0001)
2240 (PID.TID 0000.0001) grdchk output (p): 3 33 7 1 2 2 0.000000000E+00 -1.000000000E-02
2241 (PID.TID 0000.0001) grdchk output (c): 3 -4.5006507778526E-02 -4.5006601131901E-02 -4.5006414738335E-02
2242 (PID.TID 0000.0001) grdchk output (g): 3 -9.3196783089444E-06 -9.3200549988712E-06 4.0417135604986E-05
2243 (PID.TID 0000.0001)
2244 (PID.TID 0000.0001) grdchk output (p): 4 34 7 1 2 2 0.000000000E+00 -1.000000000E-02
2245 (PID.TID 0000.0001) grdchk output (c): 4 -4.5006507778526E-02 -4.5006646503147E-02 -4.5006352340135E-02
2246 (PID.TID 0000.0001) grdchk output (g): 4 -1.4708150629719E-05 -1.6684399249403E-05 1.1844889289341E-01
2247 (PID.TID 0000.0001)
2248 (PID.TID 0000.0001) grdchk summary : RMS of 4 ratios = 5.9224476828567E-02
2249 (PID.TID 0000.0001)
2250 (PID.TID 0000.0001) // =======================================================
2251 (PID.TID 0000.0001) // Gradient check results >>> END <<<
2252 (PID.TID 0000.0001) // =======================================================
2253 (PID.TID 0000.0001)
2254 (PID.TID 0000.0001) Seconds in section "ALL [THE_MODEL_MAIN]":
2255 (PID.TID 0000.0001) User time: 179.97964467667043
2256 (PID.TID 0000.0001) System time: 2.6555959191173315
2257 (PID.TID 0000.0001) Wall clock time: 183.15270400047302
2258 (PID.TID 0000.0001) No. starts: 1
2259 (PID.TID 0000.0001) No. stops: 1
2260 (PID.TID 0000.0001) Seconds in section "INITIALISE_FIXED [THE_MODEL_MAIN]":
2261 (PID.TID 0000.0001) User time: 0.20396900735795498
2262 (PID.TID 0000.0001) System time: 1.99970006942749023E-002
2263 (PID.TID 0000.0001) Wall clock time: 0.23837614059448242
2264 (PID.TID 0000.0001) No. starts: 1
2265 (PID.TID 0000.0001) No. stops: 1
2266 (PID.TID 0000.0001) Seconds in section "THE_MAIN_LOOP (F) [THE_MODEL_MAIN]":
2267 (PID.TID 0000.0001) User time: 18.144242182374001
2268 (PID.TID 0000.0001) System time: 2.4916210379451513
2269 (PID.TID 0000.0001) Wall clock time: 20.727953910827637
2270 (PID.TID 0000.0001) No. starts: 1
2271 (PID.TID 0000.0001) No. stops: 1
2272 (PID.TID 0000.0001) Seconds in section "INITIALISE_VARIA [THE_MAIN_LOOP]":
2273 (PID.TID 0000.0001) User time: 1.1838391572237015
2274 (PID.TID 0000.0001) System time: 4.89937271922826767E-002
2275 (PID.TID 0000.0001) Wall clock time: 1.2693917751312256
2276 (PID.TID 0000.0001) No. starts: 9
2277 (PID.TID 0000.0001) No. stops: 9
2278 (PID.TID 0000.0001) Seconds in section "MAIN LOOP [THE_MAIN_LOOP]":
2279 (PID.TID 0000.0001) User time: 131.62598377466202
2280 (PID.TID 0000.0001) System time: 2.4566265344619751
2281 (PID.TID 0000.0001) Wall clock time: 134.42286586761475
2282 (PID.TID 0000.0001) No. starts: 9
2283 (PID.TID 0000.0001) No. stops: 9
2284 (PID.TID 0000.0001) Seconds in section "MAIN_DO_LOOP [THE_MAIN_LOOP]":
2285 (PID.TID 0000.0001) User time: 131.62499195337296
2286 (PID.TID 0000.0001) System time: 2.4556275606155396
2287 (PID.TID 0000.0001) Wall clock time: 134.42157053947449
2288 (PID.TID 0000.0001) No. starts: 9
2289 (PID.TID 0000.0001) No. stops: 9
2290 (PID.TID 0000.0001) Seconds in section "UPDATE_SURF_DR [FORWARD_STEP]":
2291 (PID.TID 0000.0001) User time: 0.12690782546997070
2292 (PID.TID 0000.0001) System time: 0.0000000000000000
2293 (PID.TID 0000.0001) Wall clock time: 0.13155961036682129
2294 (PID.TID 0000.0001) No. starts: 198
2295 (PID.TID 0000.0001) No. stops: 198
2296 (PID.TID 0000.0001) Seconds in section "LOAD_FIELDS_DRIVER [FORWARD_STEP]":
2297 (PID.TID 0000.0001) User time: 0.11702477931976318
2298 (PID.TID 0000.0001) System time: 5.99916279315948486E-003
2299 (PID.TID 0000.0001) Wall clock time: 0.12577748298645020
2300 (PID.TID 0000.0001) No. starts: 99
2301 (PID.TID 0000.0001) No. stops: 99
2302 (PID.TID 0000.0001) Seconds in section "EXTERNAL_FLDS_LOAD [LOAD_FLDS_DRIVER]":
2303 (PID.TID 0000.0001) User time: 0.11501777172088623
2304 (PID.TID 0000.0001) System time: 5.99916279315948486E-003
2305 (PID.TID 0000.0001) Wall clock time: 0.12379670143127441
2306 (PID.TID 0000.0001) No. starts: 99
2307 (PID.TID 0000.0001) No. stops: 99
2308 (PID.TID 0000.0001) Seconds in section "DO_ATMOSPHERIC_PHYS [FORWARD_STEP]":
2309 (PID.TID 0000.0001) User time: 0.11893993616104126
2310 (PID.TID 0000.0001) System time: 0.0000000000000000
2311 (PID.TID 0000.0001) Wall clock time: 0.12167930603027344
2312 (PID.TID 0000.0001) No. starts: 99
2313 (PID.TID 0000.0001) No. stops: 99
2314 (PID.TID 0000.0001) Seconds in section "DO_OCEANIC_PHYS [FORWARD_STEP]":
2315 (PID.TID 0000.0001) User time: 22.602551639080048
2316 (PID.TID 0000.0001) System time: 0.13897990435361862
2317 (PID.TID 0000.0001) Wall clock time: 22.789039850234985
2318 (PID.TID 0000.0001) No. starts: 99
2319 (PID.TID 0000.0001) No. stops: 99
2320 (PID.TID 0000.0001) Seconds in section "DYNAMICS [FORWARD_STEP]":
2321 (PID.TID 0000.0001) User time: 72.285077214241028
2322 (PID.TID 0000.0001) System time: 0.87586705386638641
2323 (PID.TID 0000.0001) Wall clock time: 73.354413986206055
2324 (PID.TID 0000.0001) No. starts: 99
2325 (PID.TID 0000.0001) No. stops: 99
2326 (PID.TID 0000.0001) Seconds in section "SOLVE_FOR_PRESSURE [FORWARD_STEP]":
2327 (PID.TID 0000.0001) User time: 5.0732443332672119
2328 (PID.TID 0000.0001) System time: 7.99794495105743408E-003
2329 (PID.TID 0000.0001) Wall clock time: 5.0852742195129395
2330 (PID.TID 0000.0001) No. starts: 99
2331 (PID.TID 0000.0001) No. stops: 99
2332 (PID.TID 0000.0001) Seconds in section "MOM_CORRECTION_STEP [FORWARD_STEP]":
2333 (PID.TID 0000.0001) User time: 0.78988099098205566
2334 (PID.TID 0000.0001) System time: 1.99997425079345703E-003
2335 (PID.TID 0000.0001) Wall clock time: 0.80133104324340820
2336 (PID.TID 0000.0001) No. starts: 99
2337 (PID.TID 0000.0001) No. stops: 99
2338 (PID.TID 0000.0001) Seconds in section "INTEGR_CONTINUITY [FORWARD_STEP]":
2339 (PID.TID 0000.0001) User time: 2.1846430301666260
2340 (PID.TID 0000.0001) System time: 4.99904155731201172E-003
2341 (PID.TID 0000.0001) Wall clock time: 2.1873135566711426
2342 (PID.TID 0000.0001) No. starts: 99
2343 (PID.TID 0000.0001) No. stops: 99
2344 (PID.TID 0000.0001) Seconds in section "CALC_SURF_DR [FORWARD_STEP]":
2345 (PID.TID 0000.0001) User time: 0.15996372699737549
2346 (PID.TID 0000.0001) System time: 0.0000000000000000
2347 (PID.TID 0000.0001) Wall clock time: 0.16621732711791992
2348 (PID.TID 0000.0001) No. starts: 99
2349 (PID.TID 0000.0001) No. stops: 99
2350 (PID.TID 0000.0001) Seconds in section "BLOCKING_EXCHANGES [FORWARD_STEP]":
2351 (PID.TID 0000.0001) User time: 1.1338067054748535
2352 (PID.TID 0000.0001) System time: 1.99997425079345703E-003
2353 (PID.TID 0000.0001) Wall clock time: 1.1271023750305176
2354 (PID.TID 0000.0001) No. starts: 198
2355 (PID.TID 0000.0001) No. stops: 198
2356 (PID.TID 0000.0001) Seconds in section "THERMODYNAMICS [FORWARD_STEP]":
2357 (PID.TID 0000.0001) User time: 46.638914704322815
2358 (PID.TID 0000.0001) System time: 1.1318282708525658
2359 (PID.TID 0000.0001) Wall clock time: 47.902785301208496
2360 (PID.TID 0000.0001) No. starts: 99
2361 (PID.TID 0000.0001) No. stops: 99
2362 (PID.TID 0000.0001) Seconds in section "TRC_CORRECTION_STEP [FORWARD_STEP]":
2363 (PID.TID 0000.0001) User time: 0.0000000000000000
2364 (PID.TID 0000.0001) System time: 0.0000000000000000
2365 (PID.TID 0000.0001) Wall clock time: 1.09863281250000000E-003
2366 (PID.TID 0000.0001) No. starts: 99
2367 (PID.TID 0000.0001) No. stops: 99
2368 (PID.TID 0000.0001) Seconds in section "COST_TILE [FORWARD_STEP]":
2369 (PID.TID 0000.0001) User time: 1.7087935209274292
2370 (PID.TID 0000.0001) System time: 0.0000000000000000
2371 (PID.TID 0000.0001) Wall clock time: 1.7157392501831055
2372 (PID.TID 0000.0001) No. starts: 99
2373 (PID.TID 0000.0001) No. stops: 99
2374 (PID.TID 0000.0001) Seconds in section "DO_THE_MODEL_IO [FORWARD_STEP]":
2375 (PID.TID 0000.0001) User time: 5.79701662063598633E-002
2376 (PID.TID 0000.0001) System time: 4.99892234802246094E-003
2377 (PID.TID 0000.0001) Wall clock time: 6.52885437011718750E-002
2378 (PID.TID 0000.0001) No. starts: 99
2379 (PID.TID 0000.0001) No. stops: 99
2380 (PID.TID 0000.0001) Seconds in section "DO_WRITE_PICKUP [FORWARD_STEP]":
2381 (PID.TID 0000.0001) User time: 3.80020141601562500E-002
2382 (PID.TID 0000.0001) System time: 4.99892234802246094E-003
2383 (PID.TID 0000.0001) Wall clock time: 4.36275005340576172E-002
2384 (PID.TID 0000.0001) No. starts: 99
2385 (PID.TID 0000.0001) No. stops: 99
2386 (PID.TID 0000.0001) Seconds in section "COST_FINAL [ADJOINT SPIN-DOWN]":
2387 (PID.TID 0000.0001) User time: 9.91821289062500000E-004
2388 (PID.TID 0000.0001) System time: 9.98973846435546875E-004
2389 (PID.TID 0000.0001) Wall clock time: 1.02782249450683594E-003
2390 (PID.TID 0000.0001) No. starts: 9
2391 (PID.TID 0000.0001) No. stops: 9
2392 (PID.TID 0000.0001) Seconds in section "THE_MAIN_LOOP (A) [THE_MODEL_MAIN]":
2393 (PID.TID 0000.0001) User time: 46.664903640747070
2394 (PID.TID 0000.0001) System time: 0.11098384857177734
2395 (PID.TID 0000.0001) Wall clock time: 46.900297880172729
2396 (PID.TID 0000.0001) No. starts: 1
2397 (PID.TID 0000.0001) No. stops: 1
2398 (PID.TID 0000.0001) Seconds in section "I/O (WRITE) [ADJOINT LOOP]":
2399 (PID.TID 0000.0001) User time: 3.79943847656250000E-002
2400 (PID.TID 0000.0001) System time: 3.99899482727050781E-003
2401 (PID.TID 0000.0001) Wall clock time: 4.17668819427490234E-002
2402 (PID.TID 0000.0001) No. starts: 1
2403 (PID.TID 0000.0001) No. stops: 1
2404 (PID.TID 0000.0001) Seconds in section "CTRL_PACK [THE_MODEL_MAIN]":
2405 (PID.TID 0000.0001) User time: 7.89871215820312500E-002
2406 (PID.TID 0000.0001) System time: 2.99906730651855469E-003
2407 (PID.TID 0000.0001) Wall clock time: 8.11889171600341797E-002
2408 (PID.TID 0000.0001) No. starts: 1
2409 (PID.TID 0000.0001) No. stops: 1
2410 (PID.TID 0000.0001) Seconds in section "CTRL_PACK [THE_MODEL_MAIN]":
2411 (PID.TID 0000.0001) User time: 7.39898681640625000E-002
2412 (PID.TID 0000.0001) System time: 5.99908828735351563E-003
2413 (PID.TID 0000.0001) Wall clock time: 7.99939632415771484E-002
2414 (PID.TID 0000.0001) No. starts: 1
2415 (PID.TID 0000.0001) No. stops: 1
2416 (PID.TID 0000.0001) Seconds in section "GRDCHK_MAIN [THE_MODEL_MAIN]":
2417 (PID.TID 0000.0001) User time: 114.81355285644531
2418 (PID.TID 0000.0001) System time: 2.39958763122558594E-002
2419 (PID.TID 0000.0001) Wall clock time: 115.12481188774109
2420 (PID.TID 0000.0001) No. starts: 1
2421 (PID.TID 0000.0001) No. stops: 1
2422 (PID.TID 0000.0001) // ======================================================
2423 (PID.TID 0000.0001) // Tile <-> Tile communication statistics
2424 (PID.TID 0000.0001) // ======================================================
2425 (PID.TID 0000.0001) // o Tile number: 000001
2426 (PID.TID 0000.0001) // No. X exchanges = 0
2427 (PID.TID 0000.0001) // Max. X spins = 0
2428 (PID.TID 0000.0001) // Min. X spins = 1000000000
2429 (PID.TID 0000.0001) // Total. X spins = 0
2430 (PID.TID 0000.0001) // Avg. X spins = 0.00E+00
2431 (PID.TID 0000.0001) // No. Y exchanges = 0
2432 (PID.TID 0000.0001) // Max. Y spins = 0
2433 (PID.TID 0000.0001) // Min. Y spins = 1000000000
2434 (PID.TID 0000.0001) // Total. Y spins = 0
2435 (PID.TID 0000.0001) // Avg. Y spins = 0.00E+00
2436 (PID.TID 0000.0001) // o Tile number: 000002
2437 (PID.TID 0000.0001) // No. X exchanges = 0
2438 (PID.TID 0000.0001) // Max. X spins = 0
2439 (PID.TID 0000.0001) // Min. X spins = 1000000000
2440 (PID.TID 0000.0001) // Total. X spins = 0
2441 (PID.TID 0000.0001) // Avg. X spins = 0.00E+00
2442 (PID.TID 0000.0001) // No. Y exchanges = 0
2443 (PID.TID 0000.0001) // Max. Y spins = 0
2444 (PID.TID 0000.0001) // Min. Y spins = 1000000000
2445 (PID.TID 0000.0001) // Total. Y spins = 0
2446 (PID.TID 0000.0001) // Avg. Y spins = 0.00E+00
2447 (PID.TID 0000.0001) // o Tile number: 000003
2448 (PID.TID 0000.0001) // No. X exchanges = 0
2449 (PID.TID 0000.0001) // Max. X spins = 0
2450 (PID.TID 0000.0001) // Min. X spins = 1000000000
2451 (PID.TID 0000.0001) // Total. X spins = 0
2452 (PID.TID 0000.0001) // Avg. X spins = 0.00E+00
2453 (PID.TID 0000.0001) // No. Y exchanges = 0
2454 (PID.TID 0000.0001) // Max. Y spins = 0
2455 (PID.TID 0000.0001) // Min. Y spins = 1000000000
2456 (PID.TID 0000.0001) // Total. Y spins = 0
2457 (PID.TID 0000.0001) // Avg. Y spins = 0.00E+00
2458 (PID.TID 0000.0001) // o Tile number: 000004
2459 (PID.TID 0000.0001) // No. X exchanges = 0
2460 (PID.TID 0000.0001) // Max. X spins = 0
2461 (PID.TID 0000.0001) // Min. X spins = 1000000000
2462 (PID.TID 0000.0001) // Total. X spins = 0
2463 (PID.TID 0000.0001) // Avg. X spins = 0.00E+00
2464 (PID.TID 0000.0001) // No. Y exchanges = 0
2465 (PID.TID 0000.0001) // Max. Y spins = 0
2466 (PID.TID 0000.0001) // Min. Y spins = 1000000000
2467 (PID.TID 0000.0001) // Total. Y spins = 0
2468 (PID.TID 0000.0001) // Avg. Y spins = 0.00E+00
2469 (PID.TID 0000.0001) // o Thread number: 000001
2470 (PID.TID 0000.0001) // No. barriers = 17872
2471 (PID.TID 0000.0001) // Max. barrier spins = 1
2472 (PID.TID 0000.0001) // Min. barrier spins = 1
2473 (PID.TID 0000.0001) // Total barrier spins = 17872
2474 (PID.TID 0000.0001) // Avg. barrier spins = 1.00E+00
2475 PROGRAM MAIN: Execution ended Normally

  ViewVC Help
Powered by ViewVC 1.1.22