/[MITgcm]/MITgcm/verification/OpenAD/results/output_oadm.ggl90.txt
ViewVC logotype

Annotation of /MITgcm/verification/OpenAD/results/output_oadm.ggl90.txt

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


Revision 1.1 - (hide annotations) (download)
Thu Jul 2 15:57:20 2015 UTC (8 years, 10 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65n, checkpoint65o, HEAD
File MIME type: text/plain
Adding ggl90 preliminary ggl90 test to OpenAD test suite.
Config. with ggl90 currently diverges after 64 time steps.

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

  ViewVC Help
Powered by ViewVC 1.1.22