/[MITgcm]/MITgcm/verification/OpenAD/results/output_adm.txt
ViewVC logotype

Annotation of /MITgcm/verification/OpenAD/results/output_adm.txt

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

Revision 1.5 - (hide annotations) (download)
Tue Jul 29 19:02:07 2014 UTC (9 years, 9 months ago) by jmc
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, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, HEAD
Changes since 1.4: +490 -422 lines
File MIME type: text/plain 
- use default cd_code options (#undef CD_CODE_NO_AB_MOMENTUM &
  #undef CD_CODE_NO_AB_CORIOLIS) for both AD and OAD built
- and update the corresponding output.
1 jmc 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 jmc 1.5 (PID.TID 0000.0001) // MITgcmUV version: checkpoint65a
9 jmc 1.3 (PID.TID 0000.0001) // Build user: jmc
10     (PID.TID 0000.0001) // Build host: baudelaire
11 jmc 1.5 (PID.TID 0000.0001) // Build date: Tue Jul 29 13:06:16 EDT 2014
12 jmc 1.1 (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 jmc 1.3 (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 jmc 1.1 (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 jmc 1.3 (PID.TID 0000.0001) useCoupler= F ;/* Flag used to control communications with */
53 jmc 1.1 (PID.TID 0000.0001) /* other model components, through a coupler */
54 jmc 1.3 (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 jmc 1.1 (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 jmc 1.5 (PID.TID 0000.0001) // Tile number: 000001 (process no. = 000000)
67     (PID.TID 0000.0001) // WEST: Tile = 000001, Process = 000000, Comm = put
68 jmc 1.1 (PID.TID 0000.0001) // bi = 000001, bj = 000001
69 jmc 1.5 (PID.TID 0000.0001) // EAST: Tile = 000001, Process = 000000, Comm = put
70 jmc 1.1 (PID.TID 0000.0001) // bi = 000001, bj = 000001
71 jmc 1.5 (PID.TID 0000.0001) // SOUTH: Tile = 000001, Process = 000000, Comm = put
72 jmc 1.1 (PID.TID 0000.0001) // bi = 000001, bj = 000001
73 jmc 1.5 (PID.TID 0000.0001) // NORTH: Tile = 000001, Process = 000000, Comm = put
74 jmc 1.1 (PID.TID 0000.0001) // bi = 000001, bj = 000001
75     (PID.TID 0000.0001)
76 jmc 1.3 (PID.TID 0000.0001) INI_PARMS: opening model parameter file "data"
77     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data
78 jmc 1.1 (PID.TID 0000.0001) // =======================================================
79 jmc 1.3 (PID.TID 0000.0001) // Parameter file "data"
80 jmc 1.1 (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) > rhonil=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) > allowFreezing=.TRUE.,
102     (PID.TID 0000.0001) > useRealFreshWaterFlux=.TRUE.,
103     (PID.TID 0000.0001) > useCDscheme=.TRUE.,
104     (PID.TID 0000.0001) ># turn on looped cells
105     (PID.TID 0000.0001) > hFacMin=.05,
106     (PID.TID 0000.0001) > hFacMindr=50.,
107     (PID.TID 0000.0001) ># set precision of data files
108     (PID.TID 0000.0001) > readBinaryPrec=32,
109 jmc 1.3 (PID.TID 0000.0001) > /
110 jmc 1.1 (PID.TID 0000.0001) >
111     (PID.TID 0000.0001) ># Elliptic solver parameters
112     (PID.TID 0000.0001) > &PARM02
113 heimbach 1.2 (PID.TID 0000.0001) > cg2dMaxIters=1000,
114     (PID.TID 0000.0001) > cg2dTargetResidual=1.E-19,
115 jmc 1.3 (PID.TID 0000.0001) > /
116 jmc 1.1 (PID.TID 0000.0001) >
117     (PID.TID 0000.0001) ># Time stepping parameters
118     (PID.TID 0000.0001) > &PARM03
119 heimbach 1.2 (PID.TID 0000.0001) > nIter0 = 0,
120     (PID.TID 0000.0001) > nTimeSteps = 4,
121 jmc 1.1 (PID.TID 0000.0001) ># 100 years of integration will yield a reasonable flow field
122     (PID.TID 0000.0001) ># startTime = 0.,
123     (PID.TID 0000.0001) ># endTime = 3110400000.,
124     (PID.TID 0000.0001) > deltaTmom = 1200.0,
125     (PID.TID 0000.0001) > tauCD = 321428.,
126     (PID.TID 0000.0001) > deltaTtracer= 43200.0,
127     (PID.TID 0000.0001) > deltaTClock = 43200.0,
128     (PID.TID 0000.0001) ># if you are using a version later than checkpoint45d on the main branch
129     (PID.TID 0000.0001) ># you can uncomment the following line and increase the time step
130     (PID.TID 0000.0001) ># deltaTtracer and deltaTClock to 172800.0 as well to speed up the
131     (PID.TID 0000.0001) ># asynchronous time stepping
132     (PID.TID 0000.0001) ># deltaTfreesurf = 172800.0,
133     (PID.TID 0000.0001) > abEps = 0.1,
134     (PID.TID 0000.0001) > pChkptFreq= 311040000.,
135     (PID.TID 0000.0001) > dumpFreq = 2592000.,
136     (PID.TID 0000.0001) > adjDumpFreq = 2592000.,
137     (PID.TID 0000.0001) > monitorFreq = 2592000.,
138     (PID.TID 0000.0001) > adjMonitorFreq = 2592000.,
139     (PID.TID 0000.0001) ># 2 months restoring timescale for temperature
140     (PID.TID 0000.0001) > tauThetaClimRelax = 5184000.0,
141     (PID.TID 0000.0001) ># 6 months restoring timescale for salinity
142     (PID.TID 0000.0001) > tauSaltClimRelax = 15552000.0,
143     (PID.TID 0000.0001) > periodicExternalForcing=.TRUE.,
144     (PID.TID 0000.0001) > externForcingPeriod=2592000.,
145     (PID.TID 0000.0001) > externForcingCycle=31104000.,
146 jmc 1.3 (PID.TID 0000.0001) > /
147 jmc 1.1 (PID.TID 0000.0001) >
148     (PID.TID 0000.0001) ># Gridding parameters
149     (PID.TID 0000.0001) > &PARM04
150     (PID.TID 0000.0001) > usingSphericalPolarGrid=.TRUE.,
151     (PID.TID 0000.0001) > delR= 50., 70., 100., 140., 190.,
152     (PID.TID 0000.0001) > 240., 290., 340., 390., 440.,
153     (PID.TID 0000.0001) > 490., 540., 590., 640., 690.,
154 jmc 1.3 (PID.TID 0000.0001) > ygOrigin=-80.,
155 jmc 1.1 (PID.TID 0000.0001) > dySpacing=4.,
156     (PID.TID 0000.0001) > dxSpacing=4.,
157 jmc 1.3 (PID.TID 0000.0001) > /
158 jmc 1.1 (PID.TID 0000.0001) >
159     (PID.TID 0000.0001) ># Input datasets
160     (PID.TID 0000.0001) > &PARM05
161     (PID.TID 0000.0001) > bathyFile= 'bathymetry.bin',
162     (PID.TID 0000.0001) > hydrogThetaFile='lev_t.bin',
163     (PID.TID 0000.0001) > hydrogSaltFile= 'lev_s.bin',
164     (PID.TID 0000.0001) > zonalWindFile= 'trenberth_taux.bin',
165     (PID.TID 0000.0001) > meridWindFile= 'trenberth_tauy.bin',
166     (PID.TID 0000.0001) > thetaClimFile= 'lev_sst.bin',
167     (PID.TID 0000.0001) > saltClimFile= 'lev_sss.bin',
168     (PID.TID 0000.0001) > surfQFile= 'ncep_qnet.bin',
169     (PID.TID 0000.0001) ># fresh water flux is turned off, uncomment next line to turn on
170     (PID.TID 0000.0001) ># (not recommened together with surface salinity restoring)
171     (PID.TID 0000.0001) ># EmPmRFile= 'ncep_emp.bin',
172 jmc 1.3 (PID.TID 0000.0001) > /
173 jmc 1.1 (PID.TID 0000.0001)
174 jmc 1.3 (PID.TID 0000.0001) INI_PARMS ; starts to read PARM01
175     (PID.TID 0000.0001) INI_PARMS ; read PARM01 : OK
176 jmc 1.1 (PID.TID 0000.0001) S/R INI_PARMS: No request for barotropic solver
177     (PID.TID 0000.0001) S/R INI_PARMS: => Use implicitFreeSurface as default
178 jmc 1.3 (PID.TID 0000.0001) INI_PARMS ; starts to read PARM02
179     (PID.TID 0000.0001) INI_PARMS ; read PARM02 : OK
180     (PID.TID 0000.0001) INI_PARMS ; starts to read PARM03
181     (PID.TID 0000.0001) INI_PARMS ; read PARM03 : OK
182     (PID.TID 0000.0001) INI_PARMS ; starts to read PARM04
183     (PID.TID 0000.0001) INI_PARMS ; read PARM04 : OK
184     (PID.TID 0000.0001) INI_PARMS ; starts to read PARM05
185     (PID.TID 0000.0001) INI_PARMS ; read PARM05 : OK
186     (PID.TID 0000.0001) INI_PARMS: finished reading file "data"
187 jmc 1.1 (PID.TID 0000.0001) PACKAGES_BOOT: opening data.pkg
188     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.pkg
189     (PID.TID 0000.0001) // =======================================================
190     (PID.TID 0000.0001) // Parameter file "data.pkg"
191     (PID.TID 0000.0001) // =======================================================
192     (PID.TID 0000.0001) >#
193     (PID.TID 0000.0001) ># ********
194     (PID.TID 0000.0001) ># Packages
195     (PID.TID 0000.0001) ># ********
196     (PID.TID 0000.0001) > &PACKAGES
197     (PID.TID 0000.0001) > useGMRedi = .TRUE.,
198     (PID.TID 0000.0001) > useKPP = .FALSE.,
199     (PID.TID 0000.0001) > useGrdchk = .TRUE.,
200     (PID.TID 0000.0001) > useMNC = .FALSE.,
201 jmc 1.3 (PID.TID 0000.0001) > /
202 jmc 1.1 (PID.TID 0000.0001)
203     (PID.TID 0000.0001) PACKAGES_BOOT: finished reading data.pkg
204 jmc 1.5 (PID.TID 0000.0001) PACKAGES_BOOT: On/Off package Summary
205     -------- pkgs with a standard "usePKG" On/Off switch in "data.pkg": --------
206     pkg/gmredi compiled and used ( useGMRedi = T )
207     pkg/grdchk compiled and used ( useGrdchk = T )
208     pkg/mnc compiled but not used ( useMNC = F )
209     -------- pkgs without standard "usePKG" On/Off switch in "data.pkg": --------
210     pkg/generic_advdiff compiled and used ( useGAD = T )
211     pkg/mom_common compiled and used ( momStepping = T )
212     pkg/mom_vecinv compiled but not used ( +vectorInvariantMomentum = F )
213     pkg/mom_fluxform compiled and used ( & not vectorInvariantMom = T )
214     pkg/cd_code compiled and used ( useCDscheme = T )
215     pkg/monitor compiled and used ( monitorFreq > 0. = T )
216     pkg/debug compiled but not used ( debugMode = F )
217     pkg/rw compiled and used
218     pkg/mdsio compiled and used
219     pkg/autodiff compiled and used
220     pkg/cost compiled and used
221     pkg/ctrl compiled and used
222     (PID.TID 0000.0001) PACKAGES_BOOT: End of package Summary
223     (PID.TID 0000.0001)
224 jmc 1.1 (PID.TID 0000.0001) GM_READPARMS: opening data.gmredi
225     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.gmredi
226     (PID.TID 0000.0001) // =======================================================
227     (PID.TID 0000.0001) // Parameter file "data.gmredi"
228     (PID.TID 0000.0001) // =======================================================
229     (PID.TID 0000.0001) ># GM+Redi package parameters:
230     (PID.TID 0000.0001) ># GM_Small_Number :: epsilon used in computing the slope
231     (PID.TID 0000.0001) ># GM_slopeSqCutoff :: slope^2 cut-off value
232     (PID.TID 0000.0001) >
233     (PID.TID 0000.0001) >#-from MOM :
234     (PID.TID 0000.0001) ># GM_background_K: G & Mc.W diffusion coefficient
235     (PID.TID 0000.0001) ># GM_maxSlope : max slope of isopycnals
236     (PID.TID 0000.0001) ># GM_Scrit : transition for scaling diffusion coefficient
237     (PID.TID 0000.0001) ># GM_Sd : half width scaling for diffusion coefficient
238     (PID.TID 0000.0001) ># GM_taper_scheme: slope clipping or one of the tapering schemes
239     (PID.TID 0000.0001) ># GM_Kmin_horiz : horizontal diffusion minimum value
240     (PID.TID 0000.0001) >
241     (PID.TID 0000.0001) >#-Option parameters (needs to "define" options in GMREDI_OPTIONS.h")
242     (PID.TID 0000.0001) ># GM_isopycK : isopycnal diffusion coefficient (default=GM_background_K)
243     (PID.TID 0000.0001) ># GM_AdvForm : turn on GM Advective form (default=Skew flux form)
244     (PID.TID 0000.0001) >
245     (PID.TID 0000.0001) > &GM_PARM01
246     (PID.TID 0000.0001) > GM_Small_Number = 1.D-12,
247     (PID.TID 0000.0001) > GM_slopeSqCutoff = 1.D+08,
248     (PID.TID 0000.0001) > GM_AdvForm = .FALSE.,
249     (PID.TID 0000.0001) > GM_isopycK = 1.0D+3,
250     (PID.TID 0000.0001) > GM_background_K = 1.0D+3,
251     (PID.TID 0000.0001) > GM_taper_scheme = 'dm95',
252     (PID.TID 0000.0001) > GM_maxSlope = 1.D-2,
253     (PID.TID 0000.0001) > GM_Kmin_horiz = 50.,
254     (PID.TID 0000.0001) > GM_Scrit = 4.D-3,
255     (PID.TID 0000.0001) > GM_Sd = 1.D-3,
256     (PID.TID 0000.0001) ># GM_Visbeck_alpha = 1.5D-2,
257     (PID.TID 0000.0001) ># GM_Visbeck_alpha = 0.D0,
258     (PID.TID 0000.0001) ># GM_Visbeck_length = 2.D+5,
259     (PID.TID 0000.0001) ># GM_Visbeck_depth = 1.D+3,
260     (PID.TID 0000.0001) ># GM_Visbeck_maxval_K= 2.5D+3,
261 jmc 1.3 (PID.TID 0000.0001) > /
262 jmc 1.1 (PID.TID 0000.0001) >
263     (PID.TID 0000.0001) >
264     (PID.TID 0000.0001)
265     (PID.TID 0000.0001) GM_READPARMS: finished reading data.gmredi
266 jmc 1.5 (PID.TID 0000.0001) AUTODIFF_READPARMS: opening data.autodiff
267     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.autodiff
268     (PID.TID 0000.0001) // =======================================================
269     (PID.TID 0000.0001) // Parameter file "data.autodiff"
270     (PID.TID 0000.0001) // =======================================================
271     (PID.TID 0000.0001) ># =========================
272     (PID.TID 0000.0001) ># pkg AUTODIFF parameters :
273     (PID.TID 0000.0001) ># =========================
274     (PID.TID 0000.0001) ># inAdExact :: get an exact adjoint (no approximation) (def=.True.)
275     (PID.TID 0000.0001) >#
276     (PID.TID 0000.0001) > &AUTODIFF_PARM01
277     (PID.TID 0000.0001) ># inAdExact = .FALSE.,
278     (PID.TID 0000.0001) > /
279     (PID.TID 0000.0001)
280     (PID.TID 0000.0001) AUTODIFF_READPARMS: finished reading data.autodiff
281     (PID.TID 0000.0001) // ===================================
282     (PID.TID 0000.0001) // AUTODIFF parameters :
283     (PID.TID 0000.0001) // ===================================
284     (PID.TID 0000.0001) inAdExact = /* get an exact adjoint (no approximation) */
285     (PID.TID 0000.0001) T
286     (PID.TID 0000.0001) ;
287     (PID.TID 0000.0001) useKPPinAdMode = /* use KPP in adjoint mode */
288     (PID.TID 0000.0001) F
289     (PID.TID 0000.0001) ;
290     (PID.TID 0000.0001) useGMRediInAdMode = /* use GMRedi in adjoint mode */
291     (PID.TID 0000.0001) T
292     (PID.TID 0000.0001) ;
293     (PID.TID 0000.0001) useSEAICEinAdMode = /* use SEAICE in adjoint mode */
294     (PID.TID 0000.0001) F
295     (PID.TID 0000.0001) ;
296     (PID.TID 0000.0001) useGGL90inAdMode = /* use GGL90 in adjoint mode */
297     (PID.TID 0000.0001) F
298     (PID.TID 0000.0001) ;
299     (PID.TID 0000.0001) useSALT_PLUMEinAdMode = /* use SALT_PLUME in adjoint mode */
300     (PID.TID 0000.0001) F
301     (PID.TID 0000.0001) ;
302     (PID.TID 0000.0001) dumpAdVarExch = /* control adexch before dumpinp */
303     (PID.TID 0000.0001) 2
304     (PID.TID 0000.0001) ;
305     (PID.TID 0000.0001) mon_AdVarExch = /* control adexch before monitor */
306     (PID.TID 0000.0001) 2
307     (PID.TID 0000.0001) ;
308     (PID.TID 0000.0001) viscFacInAd = /* viscosity factor for adjoint */
309     (PID.TID 0000.0001) 1.000000000000000E+00
310     (PID.TID 0000.0001) ;
311     (PID.TID 0000.0001)
312 jmc 1.1 (PID.TID 0000.0001) OPTIM_READPARMS: opening data.optim
313     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.optim
314     (PID.TID 0000.0001) // =======================================================
315     (PID.TID 0000.0001) // Parameter file "data.optim"
316     (PID.TID 0000.0001) // =======================================================
317     (PID.TID 0000.0001) >#
318     (PID.TID 0000.0001) ># ********************************
319     (PID.TID 0000.0001) ># Off-line optimization parameters
320     (PID.TID 0000.0001) ># ********************************
321     (PID.TID 0000.0001) > &OPTIM
322     (PID.TID 0000.0001) > optimcycle=0,
323 jmc 1.3 (PID.TID 0000.0001) > /
324 jmc 1.1 (PID.TID 0000.0001)
325     (PID.TID 0000.0001) OPTIM_READPARMS: finished reading data.optim
326     (PID.TID 0000.0001) CTRL_READPARMS: opening data.ctrl
327     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.ctrl
328     (PID.TID 0000.0001) // =======================================================
329     (PID.TID 0000.0001) // Parameter file "data.ctrl"
330     (PID.TID 0000.0001) // =======================================================
331     (PID.TID 0000.0001) >#
332     (PID.TID 0000.0001) >#
333     (PID.TID 0000.0001) ># *********************
334     (PID.TID 0000.0001) ># ECCO controlvariables
335     (PID.TID 0000.0001) ># *********************
336     (PID.TID 0000.0001) > &CTRL_NML
337     (PID.TID 0000.0001) > xx_theta_file = 'xx_theta',
338     (PID.TID 0000.0001) > xx_salt_file = 'xx_salt',
339     (PID.TID 0000.0001) > xx_tr1_file = 'xx_tr1',
340     (PID.TID 0000.0001) > xx_hflux_file = 'xx_hflux',
341     (PID.TID 0000.0001) > xx_sflux_file = 'xx_sflux',
342     (PID.TID 0000.0001) > xx_tauu_file = 'xx_tauu',
343     (PID.TID 0000.0001) > xx_tauv_file = 'xx_tauv',
344     (PID.TID 0000.0001) > xx_diffkr_file = 'xx_diffkr',
345     (PID.TID 0000.0001) > xx_kapgm_file = 'xx_kapgm',
346 jmc 1.3 (PID.TID 0000.0001) > /
347 jmc 1.1 (PID.TID 0000.0001) >#
348     (PID.TID 0000.0001) ># *********************
349     (PID.TID 0000.0001) ># names for ctrl_pack/unpack
350     (PID.TID 0000.0001) ># *********************
351     (PID.TID 0000.0001) > &CTRL_PACKNAMES
352 jmc 1.3 (PID.TID 0000.0001) > /
353 jmc 1.1 (PID.TID 0000.0001) >
354     (PID.TID 0000.0001)
355     (PID.TID 0000.0001) CTRL_READPARMS: finished reading data.ctrl
356 jmc 1.5 (PID.TID 0000.0001) useSmoothCorrel2DinAdMode = /* use ctrlSmoothCorrel2D in adjoint mode */
357     (PID.TID 0000.0001) F
358     (PID.TID 0000.0001) ;
359 jmc 1.1 (PID.TID 0000.0001) COST_READPARMS: opening data.cost
360     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.cost
361     (PID.TID 0000.0001) // =======================================================
362     (PID.TID 0000.0001) // Parameter file "data.cost"
363     (PID.TID 0000.0001) // =======================================================
364     (PID.TID 0000.0001) >#
365     (PID.TID 0000.0001) >#
366     (PID.TID 0000.0001) ># ******************
367     (PID.TID 0000.0001) ># ECCO cost function
368     (PID.TID 0000.0001) ># ******************
369     (PID.TID 0000.0001) > &COST_NML
370     (PID.TID 0000.0001) >#
371     (PID.TID 0000.0001) > mult_tracer = 1.,
372     (PID.TID 0000.0001) > mult_test = 1.,
373     (PID.TID 0000.0001) > mult_atl = 1.,
374 jmc 1.3 (PID.TID 0000.0001) > /
375 jmc 1.1 (PID.TID 0000.0001)
376     (PID.TID 0000.0001) COST_READPARMS: finished reading data.cost
377     (PID.TID 0000.0001) GRDCHK_READPARMS: opening data.grdchk
378     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.grdchk
379     (PID.TID 0000.0001) // =======================================================
380     (PID.TID 0000.0001) // Parameter file "data.grdchk"
381     (PID.TID 0000.0001) // =======================================================
382     (PID.TID 0000.0001) >
383     (PID.TID 0000.0001) ># *******************
384     (PID.TID 0000.0001) ># ECCO gradient check
385     (PID.TID 0000.0001) ># *******************
386     (PID.TID 0000.0001) > &GRDCHK_NML
387     (PID.TID 0000.0001) > grdchk_eps = 1.d-2,
388 heimbach 1.2 (PID.TID 0000.0001) > iGloPos = 71,
389     (PID.TID 0000.0001) > jGloPos = 39,
390 jmc 1.1 (PID.TID 0000.0001) > kGloPos = 1,
391     (PID.TID 0000.0001) >### nbeg = 1,
392     (PID.TID 0000.0001) > nstep = 1,
393 heimbach 1.2 (PID.TID 0000.0001) > nend = 7,
394 jmc 1.1 (PID.TID 0000.0001) > grdchkvarindex = 1,
395 jmc 1.3 (PID.TID 0000.0001) > /
396 jmc 1.1 (PID.TID 0000.0001)
397     (PID.TID 0000.0001) GRDCHK_READPARMS: finished reading data.grdchk
398     (PID.TID 0000.0001)
399     (PID.TID 0000.0001) // =======================================================
400     (PID.TID 0000.0001) // Gradient check configuration >>> START <<<
401     (PID.TID 0000.0001) // =======================================================
402     (PID.TID 0000.0001)
403     (PID.TID 0000.0001) eps: 0.100E-01
404     (PID.TID 0000.0001) First location: 0
405 heimbach 1.2 (PID.TID 0000.0001) Last location: 7
406 jmc 1.1 (PID.TID 0000.0001) Increment: 1
407 jmc 1.4 (PID.TID 0000.0001) grdchkWhichProc: 0
408     (PID.TID 0000.0001) iLocTile = 1 , jLocTile = 1
409 jmc 1.1 (PID.TID 0000.0001)
410     (PID.TID 0000.0001) // =======================================================
411     (PID.TID 0000.0001) // Gradient check configuration >>> END <<<
412     (PID.TID 0000.0001) // =======================================================
413     (PID.TID 0000.0001)
414     (PID.TID 0000.0001) SET_PARMS: done
415     (PID.TID 0000.0001) Enter INI_VERTICAL_GRID: setInterFDr= T ; setCenterDr= F
416     (PID.TID 0000.0001) %MON XC_max = 3.5800000000000E+02
417     (PID.TID 0000.0001) %MON XC_min = 2.0000000000000E+00
418     (PID.TID 0000.0001) %MON XC_mean = 1.8000000000000E+02
419     (PID.TID 0000.0001) %MON XC_sd = 1.0391663325314E+02
420     (PID.TID 0000.0001) %MON XG_max = 3.5600000000000E+02
421 heimbach 1.2 (PID.TID 0000.0001) %MON XG_min = 0.0000000000000E+00
422     (PID.TID 0000.0001) %MON XG_mean = 1.7800000000000E+02
423     (PID.TID 0000.0001) %MON XG_sd = 1.0391663325314E+02
424 jmc 1.1 (PID.TID 0000.0001) %MON DXC_max = 4.4443898815675E+05
425     (PID.TID 0000.0001) %MON DXC_min = 9.2460385861875E+04
426     (PID.TID 0000.0001) %MON DXC_mean = 3.1372497618153E+05
427     (PID.TID 0000.0001) %MON DXC_sd = 1.1216447457560E+05
428     (PID.TID 0000.0001) %MON DXF_max = 4.4443898815675E+05
429     (PID.TID 0000.0001) %MON DXF_min = 9.2460385861875E+04
430     (PID.TID 0000.0001) %MON DXF_mean = 3.1372497618153E+05
431     (PID.TID 0000.0001) %MON DXF_sd = 1.1216447457560E+05
432     (PID.TID 0000.0001) %MON DXG_max = 4.4470989340816E+05
433     (PID.TID 0000.0001) %MON DXG_min = 7.7223062580781E+04
434     (PID.TID 0000.0001) %MON DXG_mean = 3.1353386340260E+05
435     (PID.TID 0000.0001) %MON DXG_sd = 1.1256651772502E+05
436     (PID.TID 0000.0001) %MON DXV_max = 4.4470989340816E+05
437     (PID.TID 0000.0001) %MON DXV_min = 7.7223062580781E+04
438     (PID.TID 0000.0001) %MON DXV_mean = 3.1353386340260E+05
439     (PID.TID 0000.0001) %MON DXV_sd = 1.1256651772502E+05
440     (PID.TID 0000.0001) %MON YC_max = 7.8000000000000E+01
441     (PID.TID 0000.0001) %MON YC_min = -7.8000000000000E+01
442     (PID.TID 0000.0001) %MON YC_mean = 0.0000000000000E+00
443     (PID.TID 0000.0001) %MON YC_sd = 4.6173585522461E+01
444     (PID.TID 0000.0001) %MON YG_max = 7.6000000000000E+01
445     (PID.TID 0000.0001) %MON YG_min = -8.0000000000000E+01
446 heimbach 1.2 (PID.TID 0000.0001) %MON YG_mean = -2.0000000000000E+00
447     (PID.TID 0000.0001) %MON YG_sd = 4.6173585522461E+01
448 jmc 1.1 (PID.TID 0000.0001) %MON DYC_max = 4.4470989340816E+05
449     (PID.TID 0000.0001) %MON DYC_min = 4.4470989340816E+05
450     (PID.TID 0000.0001) %MON DYC_mean = 4.4470989340814E+05
451     (PID.TID 0000.0001) %MON DYC_sd = 1.9848812371492E-08
452     (PID.TID 0000.0001) %MON DYF_max = 4.4470989340816E+05
453     (PID.TID 0000.0001) %MON DYF_min = 4.4470989340816E+05
454     (PID.TID 0000.0001) %MON DYF_mean = 4.4470989340814E+05
455     (PID.TID 0000.0001) %MON DYF_sd = 1.9848812371492E-08
456     (PID.TID 0000.0001) %MON DYG_max = 4.4470989340816E+05
457     (PID.TID 0000.0001) %MON DYG_min = 4.4470989340816E+05
458     (PID.TID 0000.0001) %MON DYG_mean = 4.4470989340814E+05
459     (PID.TID 0000.0001) %MON DYG_sd = 1.9848812371492E-08
460     (PID.TID 0000.0001) %MON DYU_max = 4.4470989340816E+05
461     (PID.TID 0000.0001) %MON DYU_min = 4.4470989340816E+05
462     (PID.TID 0000.0001) %MON DYU_mean = 4.4470989340814E+05
463     (PID.TID 0000.0001) %MON DYU_sd = 1.9848812371492E-08
464     (PID.TID 0000.0001) %MON RA_max = 1.9760627980089E+11
465     (PID.TID 0000.0001) %MON RA_min = 4.1109698667290E+10
466     (PID.TID 0000.0001) %MON RA_mean = 1.3948826965197E+11
467     (PID.TID 0000.0001) %MON RA_sd = 4.9870522472902E+10
468     (PID.TID 0000.0001) %MON RAW_max = 1.9760627980089E+11
469     (PID.TID 0000.0001) %MON RAW_min = 4.1109698667290E+10
470     (PID.TID 0000.0001) %MON RAW_mean = 1.3948826965197E+11
471     (PID.TID 0000.0001) %MON RAW_sd = 4.9870522472902E+10
472     (PID.TID 0000.0001) %MON RAS_max = 1.9772672958215E+11
473     (PID.TID 0000.0001) %MON RAS_min = 3.4334886267983E+10
474     (PID.TID 0000.0001) %MON RAS_mean = 1.3940329716694E+11
475     (PID.TID 0000.0001) %MON RAS_sd = 5.0049278732354E+10
476     (PID.TID 0000.0001) %MON RAZ_max = 1.9772672958215E+11
477     (PID.TID 0000.0001) %MON RAZ_min = 3.4334886267983E+10
478     (PID.TID 0000.0001) %MON RAZ_mean = 1.3940329716694E+11
479     (PID.TID 0000.0001) %MON RAZ_sd = 5.0049278732354E+10
480     (PID.TID 0000.0001) %MON AngleCS_max = 1.0000000000000E+00
481     (PID.TID 0000.0001) %MON AngleCS_min = 1.0000000000000E+00
482     (PID.TID 0000.0001) %MON AngleCS_mean = 1.0000000000000E+00
483     (PID.TID 0000.0001) %MON AngleCS_sd = 0.0000000000000E+00
484     (PID.TID 0000.0001) %MON AngleSN_max = 0.0000000000000E+00
485     (PID.TID 0000.0001) %MON AngleSN_min = 0.0000000000000E+00
486     (PID.TID 0000.0001) %MON AngleSN_mean = 0.0000000000000E+00
487     (PID.TID 0000.0001) %MON AngleSN_sd = 0.0000000000000E+00
488 jmc 1.3 (PID.TID 0000.0001) GAD_INIT_FIXED: GAD_OlMinSize= 1 0 1
489 jmc 1.1 (PID.TID 0000.0001)
490     (PID.TID 0000.0001) // ===================================
491     (PID.TID 0000.0001) // GAD parameters :
492     (PID.TID 0000.0001) // ===================================
493     (PID.TID 0000.0001) tempAdvScheme = /* Temp. Horiz.Advection scheme selector */
494     (PID.TID 0000.0001) 2
495     (PID.TID 0000.0001) ;
496     (PID.TID 0000.0001) tempVertAdvScheme = /* Temp. Vert. Advection scheme selector */
497     (PID.TID 0000.0001) 2
498     (PID.TID 0000.0001) ;
499     (PID.TID 0000.0001) tempMultiDimAdvec = /* use Muti-Dim Advec method for Temp */
500     (PID.TID 0000.0001) F
501     (PID.TID 0000.0001) ;
502     (PID.TID 0000.0001) tempSOM_Advection = /* use 2nd Order Moment Advection for Temp */
503     (PID.TID 0000.0001) F
504     (PID.TID 0000.0001) ;
505     (PID.TID 0000.0001) AdamsBashforthGt = /* apply Adams-Bashforth extrapolation on Gt */
506     (PID.TID 0000.0001) T
507     (PID.TID 0000.0001) ;
508     (PID.TID 0000.0001) AdamsBashforth_T = /* apply Adams-Bashforth extrapolation on Temp */
509     (PID.TID 0000.0001) F
510     (PID.TID 0000.0001) ;
511     (PID.TID 0000.0001) saltAdvScheme = /* Salt. Horiz.advection scheme selector */
512     (PID.TID 0000.0001) 2
513     (PID.TID 0000.0001) ;
514     (PID.TID 0000.0001) saltVertAdvScheme = /* Salt. Vert. Advection scheme selector */
515     (PID.TID 0000.0001) 2
516     (PID.TID 0000.0001) ;
517     (PID.TID 0000.0001) saltMultiDimAdvec = /* use Muti-Dim Advec method for Salt */
518     (PID.TID 0000.0001) F
519     (PID.TID 0000.0001) ;
520     (PID.TID 0000.0001) saltSOM_Advection = /* use 2nd Order Moment Advection for Salt */
521     (PID.TID 0000.0001) F
522     (PID.TID 0000.0001) ;
523     (PID.TID 0000.0001) AdamsBashforthGs = /* apply Adams-Bashforth extrapolation on Gs */
524     (PID.TID 0000.0001) T
525     (PID.TID 0000.0001) ;
526     (PID.TID 0000.0001) AdamsBashforth_S = /* apply Adams-Bashforth extrapolation on Salt */
527     (PID.TID 0000.0001) F
528     (PID.TID 0000.0001) ;
529     (PID.TID 0000.0001) // ===================================
530     (PID.TID 0000.0001) ctrl-wet 1: nvarlength = 117236
531     (PID.TID 0000.0001) ctrl-wet 2: surface wet C = 2315
532     (PID.TID 0000.0001) ctrl-wet 3: surface wet W = 2206
533     (PID.TID 0000.0001) ctrl-wet 4: surface wet S = 2149
534     (PID.TID 0000.0001) ctrl-wet 4a:surface wet V = 0
535     (PID.TID 0000.0001) ctrl-wet 5: 3D wet points = 29309
536     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 1 1
537     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 2 1
538     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 3 0
539     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 4 0
540     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 5 0
541     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 6 0
542     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 7 0
543     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 8 0
544     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 9 0
545     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 10 0
546     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 11 0
547     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 12 0
548     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 13 0
549     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 14 0
550     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 15 1
551     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 16 1
552     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 17 0
553     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 18 0
554     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 19 0
555     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 20 0
556     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 21 0
557     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 22 0
558     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 23 0
559     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 24 0
560     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 25 0
561     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 26 0
562     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 27 0
563     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 28 0
564     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 29 0
565     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 30 0
566     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 31 0
567     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 32 0
568     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 33 0
569     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 34 0
570     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 35 0
571     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 36 0
572     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 37 0
573     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 38 0
574     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 39 0
575     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 40 0
576     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 41 0
577     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 42 0
578     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 43 0
579     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 44 0
580     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 45 0
581     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 46 0
582     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 47 0
583     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 48 0
584     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 49 0
585     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 50 0
586     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 51 0
587     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 52 0
588     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 53 0
589     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 54 0
590     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 55 0
591     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 56 0
592     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 57 0
593     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 58 0
594     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 59 0
595     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 60 0
596     (PID.TID 0000.0001) ctrl-wet 7: flux 58618
597     (PID.TID 0000.0001) ctrl-wet 8: atmos 58618
598     (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
599     (PID.TID 0000.0001) ctrl-wet 13: global nvarlength for Nr = 15 117236
600     (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
601     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 1 2315 2149 2206 0
602     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 2 2315 2149 2206 0
603     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 3 2254 2102 2146 0
604     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 4 2215 2058 2104 0
605     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 5 2178 2027 2070 0
606     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 6 2142 1987 2029 0
607     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 7 2114 1959 2004 0
608     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 8 2076 1918 1959 0
609     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 9 2048 1887 1925 0
610     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 10 1999 1831 1869 0
611     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 11 1948 1771 1808 0
612     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 12 1850 1653 1705 0
613     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 13 1655 1404 1458 0
614     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 14 1372 1118 1164 0
615     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 15 828 623 671 0
616     (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
617     (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
618     (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
619     (PID.TID 0000.0001) ctrl_init: no. of control variables: 4
620     (PID.TID 0000.0001) ctrl_init: control vector length: 117236
621     (PID.TID 0000.0001) %MON fCori_max = 1.4226580169407E-04
622     (PID.TID 0000.0001) %MON fCori_min = -1.4226580169407E-04
623     (PID.TID 0000.0001) %MON fCori_mean = 7.5291817533716E-23
624     (PID.TID 0000.0001) %MON fCori_sd = 9.6335367303778E-05
625     (PID.TID 0000.0001) %MON fCoriG_max = 1.4112379284621E-04
626     (PID.TID 0000.0001) %MON fCoriG_min = -1.4323448157692E-04
627 heimbach 1.2 (PID.TID 0000.0001) %MON fCoriG_mean = -3.5808620394229E-06
628     (PID.TID 0000.0001) %MON fCoriG_sd = 9.6285194120965E-05
629 jmc 1.1 (PID.TID 0000.0001) %MON fCoriCos_max = 1.4535550371427E-04
630     (PID.TID 0000.0001) %MON fCoriCos_min = 3.0239529651320E-05
631     (PID.TID 0000.0001) %MON fCoriCos_mean = 1.0260497651149E-04
632     (PID.TID 0000.0001) %MON fCoriCos_sd = 3.6683828681187E-05
633     (PID.TID 0000.0001) INI_CG2D: CG2D normalisation factor = 6.5682677425711703E-05
634     (PID.TID 0000.0001)
635     (PID.TID 0000.0001) // =======================================================
636     (PID.TID 0000.0001) // Model configuration
637     (PID.TID 0000.0001) // =======================================================
638     (PID.TID 0000.0001) //
639     (PID.TID 0000.0001) // "Physical" paramters ( PARM01 in namelist )
640     (PID.TID 0000.0001) //
641     (PID.TID 0000.0001) buoyancyRelation = /* Type of relation to get Buoyancy */
642     (PID.TID 0000.0001) 'OCEANIC'
643     (PID.TID 0000.0001) ;
644 jmc 1.3 (PID.TID 0000.0001) fluidIsAir = /* fluid major constituent is Air */
645 jmc 1.1 (PID.TID 0000.0001) F
646     (PID.TID 0000.0001) ;
647 jmc 1.3 (PID.TID 0000.0001) fluidIsWater = /* fluid major constituent is Water */
648 jmc 1.1 (PID.TID 0000.0001) T
649     (PID.TID 0000.0001) ;
650     (PID.TID 0000.0001) usingPCoords = /* use p (or p*) vertical coordinate */
651     (PID.TID 0000.0001) F
652     (PID.TID 0000.0001) ;
653     (PID.TID 0000.0001) usingZCoords = /* use z (or z*) vertical coordinate */
654     (PID.TID 0000.0001) T
655     (PID.TID 0000.0001) ;
656     (PID.TID 0000.0001) tRef = /* Reference temperature profile ( oC or K ) */
657     (PID.TID 0000.0001) 15 @ 2.000000000000000E+01 /* K = 1: 15 */
658     (PID.TID 0000.0001) ;
659     (PID.TID 0000.0001) sRef = /* Reference salinity profile ( psu ) */
660     (PID.TID 0000.0001) 15 @ 3.500000000000000E+01 /* K = 1: 15 */
661     (PID.TID 0000.0001) ;
662 jmc 1.5 (PID.TID 0000.0001) useStrainTensionVisc= /* Use StrainTension Form of Viscous Operator */
663     (PID.TID 0000.0001) F
664 jmc 1.1 (PID.TID 0000.0001) ;
665 jmc 1.5 (PID.TID 0000.0001) useVariableVisc = /* Use variable horizontal viscosity */
666     (PID.TID 0000.0001) F
667 jmc 1.1 (PID.TID 0000.0001) ;
668 jmc 1.5 (PID.TID 0000.0001) useHarmonicVisc = /* Use harmonic horizontal viscosity */
669     (PID.TID 0000.0001) T
670 jmc 1.1 (PID.TID 0000.0001) ;
671 jmc 1.5 (PID.TID 0000.0001) useBiharmonicVisc= /* Use biharmonic horiz. viscosity */
672 jmc 1.1 (PID.TID 0000.0001) F
673     (PID.TID 0000.0001) ;
674 jmc 1.5 (PID.TID 0000.0001) useSmag3D = /* Use isotropic 3-D Smagorinsky viscosity */
675 jmc 1.1 (PID.TID 0000.0001) F
676     (PID.TID 0000.0001) ;
677 jmc 1.5 (PID.TID 0000.0001) viscAh = /* Lateral harmonic viscosity ( m^2/s ) */
678     (PID.TID 0000.0001) 5.000000000000000E+05
679 jmc 1.1 (PID.TID 0000.0001) ;
680 jmc 1.3 (PID.TID 0000.0001) viscA4 = /* Lateral biharmonic viscosity ( m^4/s ) */
681 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
682     (PID.TID 0000.0001) ;
683     (PID.TID 0000.0001) no_slip_sides = /* Viscous BCs: No-slip sides */
684     (PID.TID 0000.0001) T
685     (PID.TID 0000.0001) ;
686     (PID.TID 0000.0001) sideDragFactor = /* side-drag scaling factor (non-dim) */
687     (PID.TID 0000.0001) 2.000000000000000E+00
688     (PID.TID 0000.0001) ;
689 jmc 1.3 (PID.TID 0000.0001) viscArNr = /* vertical profile of vertical viscosity ( m^2/s )*/
690     (PID.TID 0000.0001) 15 @ 1.000000000000000E-03 /* K = 1: 15 */
691 jmc 1.1 (PID.TID 0000.0001) ;
692     (PID.TID 0000.0001) no_slip_bottom = /* Viscous BCs: No-slip bottom */
693     (PID.TID 0000.0001) T
694     (PID.TID 0000.0001) ;
695 jmc 1.3 (PID.TID 0000.0001) bottomDragLinear = /* linear bottom-drag coefficient ( m/s ) */
696 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
697     (PID.TID 0000.0001) ;
698 jmc 1.3 (PID.TID 0000.0001) bottomDragQuadratic = /* quadratic bottom-drag coefficient (-) */
699 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
700     (PID.TID 0000.0001) ;
701     (PID.TID 0000.0001) diffKhT = /* Laplacian diffusion of heat laterally ( m^2/s ) */
702     (PID.TID 0000.0001) 0.000000000000000E+00
703     (PID.TID 0000.0001) ;
704 jmc 1.3 (PID.TID 0000.0001) diffK4T = /* Biharmonic diffusion of heat laterally ( m^4/s ) */
705 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
706     (PID.TID 0000.0001) ;
707     (PID.TID 0000.0001) diffKhS = /* Laplacian diffusion of salt laterally ( m^2/s ) */
708     (PID.TID 0000.0001) 0.000000000000000E+00
709     (PID.TID 0000.0001) ;
710 jmc 1.3 (PID.TID 0000.0001) diffK4S = /* Biharmonic diffusion of salt laterally ( m^4/s ) */
711 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
712     (PID.TID 0000.0001) ;
713     (PID.TID 0000.0001) diffKrNrT = /* vertical profile of vertical diffusion of Temp ( m^2/s )*/
714     (PID.TID 0000.0001) 15 @ 3.000000000000000E-05 /* K = 1: 15 */
715     (PID.TID 0000.0001) ;
716     (PID.TID 0000.0001) diffKrNrS = /* vertical profile of vertical diffusion of Salt ( m^2/s )*/
717     (PID.TID 0000.0001) 15 @ 3.000000000000000E-05 /* K = 1: 15 */
718     (PID.TID 0000.0001) ;
719 jmc 1.3 (PID.TID 0000.0001) diffKrBL79surf = /* Surface diffusion for Bryan and Lewis 79 ( m^2/s ) */
720 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
721     (PID.TID 0000.0001) ;
722 jmc 1.3 (PID.TID 0000.0001) diffKrBL79deep = /* Deep diffusion for Bryan and Lewis 1979 ( m^2/s ) */
723 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
724     (PID.TID 0000.0001) ;
725 jmc 1.3 (PID.TID 0000.0001) diffKrBL79scl = /* Depth scale for Bryan and Lewis 1979 ( m ) */
726 jmc 1.1 (PID.TID 0000.0001) 2.000000000000000E+02
727     (PID.TID 0000.0001) ;
728 jmc 1.3 (PID.TID 0000.0001) diffKrBL79Ho = /* Turning depth for Bryan and Lewis 1979 ( m ) */
729 jmc 1.1 (PID.TID 0000.0001) -2.000000000000000E+03
730     (PID.TID 0000.0001) ;
731 jmc 1.3 (PID.TID 0000.0001) ivdc_kappa = /* Implicit Vertical Diffusivity for Convection ( m^2/s) */
732 jmc 1.1 (PID.TID 0000.0001) 1.000000000000000E+02
733     (PID.TID 0000.0001) ;
734     (PID.TID 0000.0001) hMixCriteria= /* Criteria for mixed-layer diagnostic */
735     (PID.TID 0000.0001) -8.000000000000000E-01
736     (PID.TID 0000.0001) ;
737 jmc 1.3 (PID.TID 0000.0001) dRhoSmall = /* Parameter for mixed-layer diagnostic */
738     (PID.TID 0000.0001) 1.000000000000000E-06
739     (PID.TID 0000.0001) ;
740     (PID.TID 0000.0001) hMixSmooth= /* Smoothing parameter for mixed-layer diagnostic */
741     (PID.TID 0000.0001) 0.000000000000000E+00
742     (PID.TID 0000.0001) ;
743 jmc 1.1 (PID.TID 0000.0001) eosType = /* Type of Equation of State */
744     (PID.TID 0000.0001) 'JMD95Z'
745     (PID.TID 0000.0001) ;
746 jmc 1.5 (PID.TID 0000.0001) HeatCapacity_Cp = /* Specific heat capacity ( J/kg/K ) */
747     (PID.TID 0000.0001) 3.994000000000000E+03
748 jmc 1.1 (PID.TID 0000.0001) ;
749 jmc 1.5 (PID.TID 0000.0001) celsius2K = /* 0 degree Celsius converted to Kelvin ( K ) */
750     (PID.TID 0000.0001) 2.731500000000000E+02
751 jmc 1.1 (PID.TID 0000.0001) ;
752 jmc 1.5 (PID.TID 0000.0001) rhoConst = /* Reference density (Boussinesq) ( kg/m^3 ) */
753 jmc 1.1 (PID.TID 0000.0001) 1.035000000000000E+03
754     (PID.TID 0000.0001) ;
755     (PID.TID 0000.0001) rhoFacC = /* normalized Reference density @ cell-Center (-) */
756     (PID.TID 0000.0001) 15 @ 1.000000000000000E+00 /* K = 1: 15 */
757     (PID.TID 0000.0001) ;
758     (PID.TID 0000.0001) rhoFacF = /* normalized Reference density @ W-Interface (-) */
759     (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
760     (PID.TID 0000.0001) ;
761 jmc 1.5 (PID.TID 0000.0001) rhoConstFresh = /* Fresh-water reference density ( kg/m^3 ) */
762 jmc 1.1 (PID.TID 0000.0001) 1.035000000000000E+03
763     (PID.TID 0000.0001) ;
764     (PID.TID 0000.0001) gravity = /* Gravitational acceleration ( m/s^2 ) */
765     (PID.TID 0000.0001) 9.810000000000000E+00
766     (PID.TID 0000.0001) ;
767     (PID.TID 0000.0001) gBaro = /* Barotropic gravity ( m/s^2 ) */
768     (PID.TID 0000.0001) 9.810000000000000E+00
769     (PID.TID 0000.0001) ;
770     (PID.TID 0000.0001) rotationPeriod = /* Rotation Period ( s ) */
771     (PID.TID 0000.0001) 8.640000000000000E+04
772     (PID.TID 0000.0001) ;
773     (PID.TID 0000.0001) omega = /* Angular velocity ( rad/s ) */
774     (PID.TID 0000.0001) 7.272205216643040E-05
775     (PID.TID 0000.0001) ;
776     (PID.TID 0000.0001) f0 = /* Reference coriolis parameter ( 1/s ) */
777     (PID.TID 0000.0001) 1.000000000000000E-04
778     (PID.TID 0000.0001) ;
779     (PID.TID 0000.0001) beta = /* Beta ( 1/(m.s) ) */
780     (PID.TID 0000.0001) 9.999999999999999E-12
781     (PID.TID 0000.0001) ;
782 jmc 1.3 (PID.TID 0000.0001) fPrime = /* Second coriolis parameter ( 1/s ) */
783     (PID.TID 0000.0001) 0.000000000000000E+00
784     (PID.TID 0000.0001) ;
785     (PID.TID 0000.0001) rigidLid = /* Rigid lid on/off flag */
786     (PID.TID 0000.0001) F
787 jmc 1.1 (PID.TID 0000.0001) ;
788     (PID.TID 0000.0001) implicitFreeSurface = /* Implicit free surface on/off flag */
789     (PID.TID 0000.0001) T
790     (PID.TID 0000.0001) ;
791 jmc 1.3 (PID.TID 0000.0001) freeSurfFac = /* Implicit free surface factor */
792     (PID.TID 0000.0001) 1.000000000000000E+00
793 jmc 1.1 (PID.TID 0000.0001) ;
794 jmc 1.3 (PID.TID 0000.0001) implicSurfPress = /* Surface Pressure implicit factor (0-1)*/
795 jmc 1.1 (PID.TID 0000.0001) 1.000000000000000E+00
796     (PID.TID 0000.0001) ;
797 jmc 1.3 (PID.TID 0000.0001) implicDiv2Dflow = /* Barot. Flow Div. implicit factor (0-1)*/
798 jmc 1.1 (PID.TID 0000.0001) 1.000000000000000E+00
799     (PID.TID 0000.0001) ;
800 jmc 1.5 (PID.TID 0000.0001) uniformLin_PhiSurf = /* use uniform Bo_surf on/off flag*/
801     (PID.TID 0000.0001) T
802 jmc 1.1 (PID.TID 0000.0001) ;
803 jmc 1.5 (PID.TID 0000.0001) uniformFreeSurfLev = /* free-surface level-index is uniform */
804 jmc 1.1 (PID.TID 0000.0001) T
805     (PID.TID 0000.0001) ;
806 jmc 1.3 (PID.TID 0000.0001) hFacMin = /* minimum partial cell factor (hFac) */
807     (PID.TID 0000.0001) 5.000000000000000E-02
808     (PID.TID 0000.0001) ;
809     (PID.TID 0000.0001) hFacMinDr = /* minimum partial cell thickness ( m) */
810 jmc 1.5 (PID.TID 0000.0001) 5.000000000000000E+01
811     (PID.TID 0000.0001) ;
812     (PID.TID 0000.0001) exactConserv = /* Exact Volume Conservation on/off flag*/
813     (PID.TID 0000.0001) F
814     (PID.TID 0000.0001) ;
815     (PID.TID 0000.0001) linFSConserveTr = /* Tracer correction for Lin Free Surface on/off flag*/
816     (PID.TID 0000.0001) F
817 jmc 1.3 (PID.TID 0000.0001) ;
818     (PID.TID 0000.0001) nonlinFreeSurf = /* Non-linear Free Surf. options (-1,0,1,2,3)*/
819 jmc 1.1 (PID.TID 0000.0001) 0
820 jmc 1.3 (PID.TID 0000.0001) -1,0= Off ; 1,2,3= On, 2=+rescale gU,gV, 3=+update cg2d solv.
821 jmc 1.1 (PID.TID 0000.0001) ;
822     (PID.TID 0000.0001) hFacInf = /* lower threshold for hFac (nonlinFreeSurf only)*/
823     (PID.TID 0000.0001) 2.000000000000000E-01
824     (PID.TID 0000.0001) ;
825     (PID.TID 0000.0001) hFacSup = /* upper threshold for hFac (nonlinFreeSurf only)*/
826     (PID.TID 0000.0001) 2.000000000000000E+00
827     (PID.TID 0000.0001) ;
828 jmc 1.3 (PID.TID 0000.0001) select_rStar = /* r* Vertical coord. options (=0 r coord.; >0 uses r*)*/
829     (PID.TID 0000.0001) 0
830     (PID.TID 0000.0001) ;
831     (PID.TID 0000.0001) useRealFreshWaterFlux = /* Real Fresh Water Flux on/off flag*/
832 jmc 1.1 (PID.TID 0000.0001) T
833     (PID.TID 0000.0001) ;
834 jmc 1.3 (PID.TID 0000.0001) temp_EvPrRn = /* Temp. of Evap/Prec/R (UNSET=use local T)(oC)*/
835     (PID.TID 0000.0001) 1.234567000000000E+05
836     (PID.TID 0000.0001) ;
837     (PID.TID 0000.0001) salt_EvPrRn = /* Salin. of Evap/Prec/R (UNSET=use local S)(psu)*/
838     (PID.TID 0000.0001) 0.000000000000000E+00
839     (PID.TID 0000.0001) ;
840 jmc 1.5 (PID.TID 0000.0001) selectAddFluid = /* option for mass source/sink of fluid (=0: off) */
841     (PID.TID 0000.0001) 0
842     (PID.TID 0000.0001) ;
843 jmc 1.3 (PID.TID 0000.0001) temp_addMass = /* Temp. of addMass array (UNSET=use local T)(oC)*/
844     (PID.TID 0000.0001) 1.234567000000000E+05
845     (PID.TID 0000.0001) ;
846     (PID.TID 0000.0001) salt_addMass = /* Salin. of addMass array (UNSET=use local S)(psu)*/
847     (PID.TID 0000.0001) 0.000000000000000E+00
848     (PID.TID 0000.0001) ;
849     (PID.TID 0000.0001) convertFW2Salt = /* convert F.W. Flux to Salt Flux (-1=use local S)(psu)*/
850 jmc 1.1 (PID.TID 0000.0001) -1.000000000000000E+00
851     (PID.TID 0000.0001) ;
852     (PID.TID 0000.0001) use3Dsolver = /* use 3-D pressure solver on/off flag */
853     (PID.TID 0000.0001) F
854     (PID.TID 0000.0001) ;
855     (PID.TID 0000.0001) nonHydrostatic = /* Non-Hydrostatic on/off flag */
856     (PID.TID 0000.0001) F
857     (PID.TID 0000.0001) ;
858     (PID.TID 0000.0001) nh_Am2 = /* Non-Hydrostatic terms scaling factor */
859     (PID.TID 0000.0001) 1.000000000000000E+00
860     (PID.TID 0000.0001) ;
861 jmc 1.3 (PID.TID 0000.0001) implicitNHPress = /* Non-Hyd Pressure implicit factor (0-1)*/
862     (PID.TID 0000.0001) 1.000000000000000E+00
863     (PID.TID 0000.0001) ;
864     (PID.TID 0000.0001) selectNHfreeSurf = /* Non-Hyd (free-)Surface option */
865     (PID.TID 0000.0001) 0
866     (PID.TID 0000.0001) ;
867 jmc 1.1 (PID.TID 0000.0001) quasiHydrostatic = /* Quasi-Hydrostatic on/off flag */
868     (PID.TID 0000.0001) F
869     (PID.TID 0000.0001) ;
870 jmc 1.5 (PID.TID 0000.0001) calc_wVelocity = /* vertical velocity calculation on/off flag */
871     (PID.TID 0000.0001) T
872     (PID.TID 0000.0001) ;
873 jmc 1.1 (PID.TID 0000.0001) momStepping = /* Momentum equation on/off flag */
874     (PID.TID 0000.0001) T
875     (PID.TID 0000.0001) ;
876     (PID.TID 0000.0001) vectorInvariantMomentum= /* Vector-Invariant Momentum on/off */
877     (PID.TID 0000.0001) F
878     (PID.TID 0000.0001) ;
879     (PID.TID 0000.0001) momAdvection = /* Momentum advection on/off flag */
880     (PID.TID 0000.0001) T
881     (PID.TID 0000.0001) ;
882     (PID.TID 0000.0001) momViscosity = /* Momentum viscosity on/off flag */
883     (PID.TID 0000.0001) T
884     (PID.TID 0000.0001) ;
885 jmc 1.3 (PID.TID 0000.0001) momImplVertAdv= /* Momentum implicit vert. advection on/off*/
886 jmc 1.1 (PID.TID 0000.0001) F
887     (PID.TID 0000.0001) ;
888     (PID.TID 0000.0001) implicitViscosity = /* Implicit viscosity on/off flag */
889     (PID.TID 0000.0001) F
890     (PID.TID 0000.0001) ;
891     (PID.TID 0000.0001) metricTerms = /* metric-Terms on/off flag */
892     (PID.TID 0000.0001) T
893     (PID.TID 0000.0001) ;
894     (PID.TID 0000.0001) useNHMTerms = /* Non-Hydrostatic Metric-Terms on/off */
895     (PID.TID 0000.0001) F
896     (PID.TID 0000.0001) ;
897 jmc 1.3 (PID.TID 0000.0001) selectCoriMap = /* Coriolis Map options (0,1,2,3)*/
898     (PID.TID 0000.0001) 2
899     (PID.TID 0000.0001) 0= f-Plane ; 1= Beta-Plane ; 2= Spherical ; 3= read from file
900 jmc 1.1 (PID.TID 0000.0001) ;
901     (PID.TID 0000.0001) use3dCoriolis = /* 3-D Coriolis on/off flag */
902     (PID.TID 0000.0001) F
903     (PID.TID 0000.0001) ;
904     (PID.TID 0000.0001) useCoriolis = /* Coriolis on/off flag */
905     (PID.TID 0000.0001) T
906     (PID.TID 0000.0001) ;
907     (PID.TID 0000.0001) useCDscheme = /* CD scheme on/off flag */
908     (PID.TID 0000.0001) T
909     (PID.TID 0000.0001) ;
910 jmc 1.3 (PID.TID 0000.0001) useEnergyConservingCoriolis= /* Flx-Form Coriolis scheme flag */
911     (PID.TID 0000.0001) F
912     (PID.TID 0000.0001) ;
913 jmc 1.1 (PID.TID 0000.0001) useJamartWetPoints= /* Coriolis WetPoints method flag */
914     (PID.TID 0000.0001) F
915     (PID.TID 0000.0001) ;
916 jmc 1.3 (PID.TID 0000.0001) useJamartMomAdv= /* V.I Non-linear terms Jamart flag */
917 jmc 1.1 (PID.TID 0000.0001) F
918     (PID.TID 0000.0001) ;
919 jmc 1.3 (PID.TID 0000.0001) useAbsVorticity= /* V.I Works with f+zeta in Coriolis */
920 jmc 1.1 (PID.TID 0000.0001) F
921     (PID.TID 0000.0001) ;
922 jmc 1.3 (PID.TID 0000.0001) selectVortScheme= /* V.I Scheme selector for Vorticity-Term */
923     (PID.TID 0000.0001) 123456789
924     (PID.TID 0000.0001) = 0 : enstrophy (Shallow-Water Eq.) conserving scheme by Sadourny, JAS 75
925     (PID.TID 0000.0001) = 1 : same as 0 with modified hFac
926     (PID.TID 0000.0001) = 2 : energy conserving scheme (used by Sadourny in JAS 75 paper)
927     (PID.TID 0000.0001) = 3 : energy (general) and enstrophy (2D, nonDiv.) conserving scheme
928     (PID.TID 0000.0001) from Sadourny (Burridge & Haseler, ECMWF Rep.4, 1977)
929 jmc 1.1 (PID.TID 0000.0001) ;
930 jmc 1.3 (PID.TID 0000.0001) upwindVorticity= /* V.I Upwind bias vorticity flag */
931 jmc 1.1 (PID.TID 0000.0001) F
932     (PID.TID 0000.0001) ;
933 jmc 1.3 (PID.TID 0000.0001) highOrderVorticity= /* V.I High order vort. advect. flag */
934 jmc 1.1 (PID.TID 0000.0001) F
935     (PID.TID 0000.0001) ;
936 jmc 1.3 (PID.TID 0000.0001) upwindShear= /* V.I Upwind vertical Shear advection flag */
937 jmc 1.1 (PID.TID 0000.0001) F
938     (PID.TID 0000.0001) ;
939 jmc 1.3 (PID.TID 0000.0001) selectKEscheme= /* V.I Kinetic Energy scheme selector */
940 jmc 1.1 (PID.TID 0000.0001) 0
941     (PID.TID 0000.0001) ;
942     (PID.TID 0000.0001) momForcing = /* Momentum forcing on/off flag */
943     (PID.TID 0000.0001) T
944     (PID.TID 0000.0001) ;
945     (PID.TID 0000.0001) momPressureForcing = /* Momentum pressure term on/off flag */
946     (PID.TID 0000.0001) T
947     (PID.TID 0000.0001) ;
948     (PID.TID 0000.0001) implicitIntGravWave= /* Implicit Internal Gravity Wave flag */
949     (PID.TID 0000.0001) F
950     (PID.TID 0000.0001) ;
951 jmc 1.5 (PID.TID 0000.0001) staggerTimeStep = /* Stagger time stepping on/off flag */
952     (PID.TID 0000.0001) F
953     (PID.TID 0000.0001) ;
954     (PID.TID 0000.0001) doResetHFactors = /* reset thickness factors @ each time-step */
955 jmc 1.1 (PID.TID 0000.0001) F
956     (PID.TID 0000.0001) ;
957 jmc 1.5 (PID.TID 0000.0001) multiDimAdvection = /* enable/disable Multi-Dim Advection */
958 jmc 1.1 (PID.TID 0000.0001) T
959     (PID.TID 0000.0001) ;
960     (PID.TID 0000.0001) useMultiDimAdvec = /* Multi-Dim Advection is/is-not used */
961     (PID.TID 0000.0001) F
962     (PID.TID 0000.0001) ;
963 jmc 1.3 (PID.TID 0000.0001) implicitDiffusion = /* Implicit Diffusion on/off flag */
964 jmc 1.1 (PID.TID 0000.0001) T
965     (PID.TID 0000.0001) ;
966     (PID.TID 0000.0001) tempStepping = /* Temperature equation on/off flag */
967     (PID.TID 0000.0001) T
968     (PID.TID 0000.0001) ;
969 jmc 1.5 (PID.TID 0000.0001) tempAdvection = /* Temperature advection on/off flag */
970 jmc 1.1 (PID.TID 0000.0001) T
971     (PID.TID 0000.0001) ;
972 jmc 1.3 (PID.TID 0000.0001) tempImplVertAdv = /* Temp. implicit vert. advection on/off */
973 jmc 1.1 (PID.TID 0000.0001) F
974     (PID.TID 0000.0001) ;
975     (PID.TID 0000.0001) tempForcing = /* Temperature forcing on/off flag */
976     (PID.TID 0000.0001) T
977     (PID.TID 0000.0001) ;
978 jmc 1.5 (PID.TID 0000.0001) doThetaClimRelax = /* apply SST relaxation on/off flag */
979     (PID.TID 0000.0001) T
980     (PID.TID 0000.0001) ;
981 jmc 1.3 (PID.TID 0000.0001) tempIsActiveTr = /* Temp. is a dynamically Active Tracer */
982     (PID.TID 0000.0001) T
983     (PID.TID 0000.0001) ;
984 jmc 1.1 (PID.TID 0000.0001) saltStepping = /* Salinity equation on/off flag */
985     (PID.TID 0000.0001) T
986     (PID.TID 0000.0001) ;
987 jmc 1.5 (PID.TID 0000.0001) saltAdvection = /* Salinity advection on/off flag */
988 jmc 1.1 (PID.TID 0000.0001) T
989     (PID.TID 0000.0001) ;
990 jmc 1.3 (PID.TID 0000.0001) saltImplVertAdv = /* Sali. implicit vert. advection on/off */
991 jmc 1.1 (PID.TID 0000.0001) F
992     (PID.TID 0000.0001) ;
993     (PID.TID 0000.0001) saltForcing = /* Salinity forcing on/off flag */
994     (PID.TID 0000.0001) T
995     (PID.TID 0000.0001) ;
996 jmc 1.5 (PID.TID 0000.0001) doSaltClimRelax = /* apply SSS relaxation on/off flag */
997     (PID.TID 0000.0001) T
998     (PID.TID 0000.0001) ;
999 jmc 1.3 (PID.TID 0000.0001) saltIsActiveTr = /* Salt is a dynamically Active Tracer */
1000     (PID.TID 0000.0001) T
1001     (PID.TID 0000.0001) ;
1002 jmc 1.1 (PID.TID 0000.0001) readBinaryPrec = /* Precision used for reading binary files */
1003     (PID.TID 0000.0001) 32
1004     (PID.TID 0000.0001) ;
1005     (PID.TID 0000.0001) writeBinaryPrec = /* Precision used for writing binary files */
1006     (PID.TID 0000.0001) 32
1007     (PID.TID 0000.0001) ;
1008     (PID.TID 0000.0001) globalFiles = /* write "global" (=not per tile) files */
1009     (PID.TID 0000.0001) F
1010     (PID.TID 0000.0001) ;
1011     (PID.TID 0000.0001) useSingleCpuIO = /* only master MPI process does I/O */
1012     (PID.TID 0000.0001) F
1013     (PID.TID 0000.0001) ;
1014 jmc 1.5 (PID.TID 0000.0001) useSingleCpuInput = /* only master process reads input */
1015     (PID.TID 0000.0001) F
1016     (PID.TID 0000.0001) ;
1017 jmc 1.3 (PID.TID 0000.0001) /* debLev[*] : level of debug & auxiliary message printing */
1018     (PID.TID 0000.0001) debLevZero = 0 ; /* level of disabled aux. msg printing */
1019     (PID.TID 0000.0001) debLevA = 1 ; /* level of minimum aux. msg printing */
1020     (PID.TID 0000.0001) debLevB = 2 ; /* level of low aux. print (report read-file opening)*/
1021     (PID.TID 0000.0001) debLevC = 3 ; /* level of moderate debug prt (most pkgs debug msg) */
1022     (PID.TID 0000.0001) debLevD = 4 ; /* level of enhanced debug prt (add DEBUG_STATS prt) */
1023     (PID.TID 0000.0001) debLevE = 5 ; /* level of extensive debug printing */
1024     (PID.TID 0000.0001) debugLevel = /* select debug printing level */
1025 jmc 1.1 (PID.TID 0000.0001) 1
1026     (PID.TID 0000.0001) ;
1027     (PID.TID 0000.0001) //
1028     (PID.TID 0000.0001) // Elliptic solver(s) paramters ( PARM02 in namelist )
1029     (PID.TID 0000.0001) //
1030     (PID.TID 0000.0001) cg2dMaxIters = /* Upper limit on 2d con. grad iterations */
1031 heimbach 1.2 (PID.TID 0000.0001) 1000
1032 jmc 1.1 (PID.TID 0000.0001) ;
1033     (PID.TID 0000.0001) cg2dChkResFreq = /* 2d con. grad convergence test frequency */
1034     (PID.TID 0000.0001) 1
1035     (PID.TID 0000.0001) ;
1036 jmc 1.5 (PID.TID 0000.0001) cg2dUseMinResSol= /* use cg2d last-iter(=0) / min-resid.(=1) solution */
1037     (PID.TID 0000.0001) 0
1038     (PID.TID 0000.0001) ;
1039 jmc 1.1 (PID.TID 0000.0001) cg2dTargetResidual = /* 2d con. grad target residual */
1040 heimbach 1.2 (PID.TID 0000.0001) 1.000000000000000E-19
1041 jmc 1.1 (PID.TID 0000.0001) ;
1042     (PID.TID 0000.0001) cg2dTargetResWunit = /* CG2d target residual [W units] */
1043     (PID.TID 0000.0001) -1.000000000000000E+00
1044     (PID.TID 0000.0001) ;
1045     (PID.TID 0000.0001) cg2dPreCondFreq = /* Freq. for updating cg2d preconditioner */
1046     (PID.TID 0000.0001) 1
1047     (PID.TID 0000.0001) ;
1048 jmc 1.3 (PID.TID 0000.0001) useSRCGSolver = /* use single reduction CG solver(s) */
1049     (PID.TID 0000.0001) F
1050     (PID.TID 0000.0001) ;
1051     (PID.TID 0000.0001) printResidualFreq = /* Freq. for printing CG residual */
1052     (PID.TID 0000.0001) 0
1053     (PID.TID 0000.0001) ;
1054 jmc 1.1 (PID.TID 0000.0001) //
1055     (PID.TID 0000.0001) // Time stepping paramters ( PARM03 in namelist )
1056     (PID.TID 0000.0001) //
1057 jmc 1.5 (PID.TID 0000.0001) deltaTMom = /* Momentum equation timestep ( s ) */
1058 jmc 1.1 (PID.TID 0000.0001) 1.200000000000000E+03
1059     (PID.TID 0000.0001) ;
1060 jmc 1.5 (PID.TID 0000.0001) deltaTFreeSurf = /* FreeSurface equation timestep ( s ) */
1061 jmc 1.1 (PID.TID 0000.0001) 1.200000000000000E+03
1062     (PID.TID 0000.0001) ;
1063 jmc 1.3 (PID.TID 0000.0001) dTtracerLev = /* Tracer equation timestep ( s ) */
1064 jmc 1.1 (PID.TID 0000.0001) 15 @ 4.320000000000000E+04 /* K = 1: 15 */
1065     (PID.TID 0000.0001) ;
1066     (PID.TID 0000.0001) deltaTClock = /* Model clock timestep ( s ) */
1067     (PID.TID 0000.0001) 4.320000000000000E+04
1068     (PID.TID 0000.0001) ;
1069     (PID.TID 0000.0001) cAdjFreq = /* Convective adjustment interval ( s ) */
1070     (PID.TID 0000.0001) 0.000000000000000E+00
1071     (PID.TID 0000.0001) ;
1072     (PID.TID 0000.0001) momForcingOutAB = /* =1: take Momentum Forcing out of Adams-Bash. stepping */
1073     (PID.TID 0000.0001) 0
1074     (PID.TID 0000.0001) ;
1075     (PID.TID 0000.0001) tracForcingOutAB = /* =1: take T,S,pTr Forcing out of Adams-Bash. stepping */
1076     (PID.TID 0000.0001) 0
1077     (PID.TID 0000.0001) ;
1078     (PID.TID 0000.0001) momDissip_In_AB = /* put Dissipation Tendency in Adams-Bash. stepping */
1079     (PID.TID 0000.0001) T
1080     (PID.TID 0000.0001) ;
1081     (PID.TID 0000.0001) doAB_onGtGs = /* apply AB on Tendencies (rather than on T,S)*/
1082     (PID.TID 0000.0001) T
1083     (PID.TID 0000.0001) ;
1084     (PID.TID 0000.0001) abEps = /* Adams-Bashforth-2 stabilizing weight */
1085     (PID.TID 0000.0001) 1.000000000000000E-01
1086     (PID.TID 0000.0001) ;
1087     (PID.TID 0000.0001) tauCD = /* CD coupling time-scale ( s ) */
1088     (PID.TID 0000.0001) 3.214280000000000E+05
1089     (PID.TID 0000.0001) ;
1090     (PID.TID 0000.0001) rCD = /* Normalised CD coupling parameter */
1091     (PID.TID 0000.0001) 9.962666600296178E-01
1092     (PID.TID 0000.0001) ;
1093 jmc 1.3 (PID.TID 0000.0001) epsAB_CD = /* AB-2 stabilizing weight for CD-scheme*/
1094     (PID.TID 0000.0001) 1.000000000000000E-01
1095     (PID.TID 0000.0001) ;
1096 heimbach 1.2 (PID.TID 0000.0001) pickupStrictlyMatch= /* stop if pickup do not strictly match */
1097     (PID.TID 0000.0001) T
1098     (PID.TID 0000.0001) ;
1099 jmc 1.3 (PID.TID 0000.0001) nIter0 = /* Run starting timestep number */
1100 heimbach 1.2 (PID.TID 0000.0001) 0
1101     (PID.TID 0000.0001) ;
1102 jmc 1.3 (PID.TID 0000.0001) nTimeSteps = /* Number of timesteps */
1103     (PID.TID 0000.0001) 4
1104     (PID.TID 0000.0001) ;
1105     (PID.TID 0000.0001) nEndIter = /* Run ending timestep number */
1106 heimbach 1.2 (PID.TID 0000.0001) 4
1107     (PID.TID 0000.0001) ;
1108 jmc 1.3 (PID.TID 0000.0001) baseTime = /* Model base time ( s ) */
1109 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
1110     (PID.TID 0000.0001) ;
1111 jmc 1.3 (PID.TID 0000.0001) startTime = /* Run start time ( s ) */
1112 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
1113     (PID.TID 0000.0001) ;
1114 jmc 1.3 (PID.TID 0000.0001) endTime = /* Integration ending time ( s ) */
1115 heimbach 1.2 (PID.TID 0000.0001) 1.728000000000000E+05
1116 jmc 1.1 (PID.TID 0000.0001) ;
1117 jmc 1.3 (PID.TID 0000.0001) pChkPtFreq = /* Permanent restart/pickup file interval ( s ) */
1118 jmc 1.1 (PID.TID 0000.0001) 3.110400000000000E+08
1119     (PID.TID 0000.0001) ;
1120 jmc 1.3 (PID.TID 0000.0001) chkPtFreq = /* Rolling restart/pickup file interval ( s ) */
1121 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
1122     (PID.TID 0000.0001) ;
1123     (PID.TID 0000.0001) pickup_write_mdsio = /* Model IO flag. */
1124     (PID.TID 0000.0001) T
1125     (PID.TID 0000.0001) ;
1126     (PID.TID 0000.0001) pickup_read_mdsio = /* Model IO flag. */
1127     (PID.TID 0000.0001) T
1128     (PID.TID 0000.0001) ;
1129     (PID.TID 0000.0001) pickup_write_mnc = /* Model IO flag. */
1130     (PID.TID 0000.0001) F
1131     (PID.TID 0000.0001) ;
1132     (PID.TID 0000.0001) pickup_read_mnc = /* Model IO flag. */
1133     (PID.TID 0000.0001) F
1134     (PID.TID 0000.0001) ;
1135     (PID.TID 0000.0001) pickup_write_immed = /* Model IO flag. */
1136     (PID.TID 0000.0001) F
1137     (PID.TID 0000.0001) ;
1138 jmc 1.3 (PID.TID 0000.0001) writePickupAtEnd = /* Model IO flag. */
1139     (PID.TID 0000.0001) T
1140     (PID.TID 0000.0001) ;
1141 jmc 1.1 (PID.TID 0000.0001) dumpFreq = /* Model state write out interval ( s ). */
1142     (PID.TID 0000.0001) 2.592000000000000E+06
1143     (PID.TID 0000.0001) ;
1144     (PID.TID 0000.0001) dumpInitAndLast= /* write out Initial & Last iter. model state */
1145     (PID.TID 0000.0001) T
1146     (PID.TID 0000.0001) ;
1147     (PID.TID 0000.0001) snapshot_mdsio = /* Model IO flag. */
1148     (PID.TID 0000.0001) T
1149     (PID.TID 0000.0001) ;
1150     (PID.TID 0000.0001) snapshot_mnc = /* Model IO flag. */
1151     (PID.TID 0000.0001) F
1152     (PID.TID 0000.0001) ;
1153     (PID.TID 0000.0001) monitorFreq = /* Monitor output interval ( s ). */
1154     (PID.TID 0000.0001) 2.592000000000000E+06
1155     (PID.TID 0000.0001) ;
1156 jmc 1.3 (PID.TID 0000.0001) monitorSelect = /* select group of variables to monitor */
1157     (PID.TID 0000.0001) 3
1158     (PID.TID 0000.0001) ;
1159 jmc 1.1 (PID.TID 0000.0001) monitor_stdio = /* Model IO flag. */
1160     (PID.TID 0000.0001) T
1161     (PID.TID 0000.0001) ;
1162     (PID.TID 0000.0001) monitor_mnc = /* Model IO flag. */
1163     (PID.TID 0000.0001) F
1164     (PID.TID 0000.0001) ;
1165     (PID.TID 0000.0001) externForcingPeriod = /* forcing period (s) */
1166     (PID.TID 0000.0001) 2.592000000000000E+06
1167     (PID.TID 0000.0001) ;
1168     (PID.TID 0000.0001) externForcingCycle = /* period of the cyle (s). */
1169     (PID.TID 0000.0001) 3.110400000000000E+07
1170     (PID.TID 0000.0001) ;
1171     (PID.TID 0000.0001) tauThetaClimRelax = /* relaxation time scale (s) */
1172     (PID.TID 0000.0001) 5.184000000000000E+06
1173     (PID.TID 0000.0001) ;
1174     (PID.TID 0000.0001) tauSaltClimRelax = /* relaxation time scale (s) */
1175     (PID.TID 0000.0001) 1.555200000000000E+07
1176     (PID.TID 0000.0001) ;
1177     (PID.TID 0000.0001) latBandClimRelax = /* max. Lat. where relaxation */
1178     (PID.TID 0000.0001) 1.800000000000000E+02
1179     (PID.TID 0000.0001) ;
1180     (PID.TID 0000.0001) //
1181     (PID.TID 0000.0001) // Gridding paramters ( PARM04 in namelist )
1182     (PID.TID 0000.0001) //
1183     (PID.TID 0000.0001) usingCartesianGrid = /* Cartesian coordinates flag ( True/False ) */
1184     (PID.TID 0000.0001) F
1185     (PID.TID 0000.0001) ;
1186     (PID.TID 0000.0001) usingCylindricalGrid = /* Cylindrical coordinates flag ( True/False ) */
1187     (PID.TID 0000.0001) F
1188     (PID.TID 0000.0001) ;
1189     (PID.TID 0000.0001) usingSphericalPolarGrid = /* Spherical coordinates flag ( True/False ) */
1190     (PID.TID 0000.0001) T
1191     (PID.TID 0000.0001) ;
1192     (PID.TID 0000.0001) usingCurvilinearGrid = /* Curvilinear coordinates flag ( True/False ) */
1193     (PID.TID 0000.0001) F
1194     (PID.TID 0000.0001) ;
1195 jmc 1.3 (PID.TID 0000.0001) selectSigmaCoord = /* Hybrid-Sigma Vert. Coordinate option */
1196     (PID.TID 0000.0001) 0
1197     (PID.TID 0000.0001) ;
1198     (PID.TID 0000.0001) Ro_SeaLevel = /* r(1) ( units of r == m ) */
1199 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
1200     (PID.TID 0000.0001) ;
1201 jmc 1.3 (PID.TID 0000.0001) rSigmaBnd = /* r/sigma transition ( units of r == m ) */
1202     (PID.TID 0000.0001) 1.234567000000000E+05
1203     (PID.TID 0000.0001) ;
1204 jmc 1.1 (PID.TID 0000.0001) rkSign = /* index orientation relative to vertical coordinate */
1205     (PID.TID 0000.0001) -1.000000000000000E+00
1206     (PID.TID 0000.0001) ;
1207     (PID.TID 0000.0001) gravitySign = /* gravity orientation relative to vertical coordinate */
1208     (PID.TID 0000.0001) -1.000000000000000E+00
1209     (PID.TID 0000.0001) ;
1210     (PID.TID 0000.0001) mass2rUnit = /* convert mass per unit area [kg/m2] to r-units [m] */
1211     (PID.TID 0000.0001) 9.661835748792270E-04
1212     (PID.TID 0000.0001) ;
1213     (PID.TID 0000.0001) rUnit2mass = /* convert r-units [m] to mass per unit area [kg/m2] */
1214     (PID.TID 0000.0001) 1.035000000000000E+03
1215     (PID.TID 0000.0001) ;
1216     (PID.TID 0000.0001) drC = /* C spacing ( units of r ) */
1217     (PID.TID 0000.0001) 2.500000000000000E+01, /* K = 1 */
1218     (PID.TID 0000.0001) 6.000000000000000E+01, /* K = 2 */
1219     (PID.TID 0000.0001) 8.500000000000000E+01, /* K = 3 */
1220     (PID.TID 0000.0001) 1.200000000000000E+02, /* K = 4 */
1221     (PID.TID 0000.0001) 1.650000000000000E+02, /* K = 5 */
1222     (PID.TID 0000.0001) 2.150000000000000E+02, /* K = 6 */
1223     (PID.TID 0000.0001) 2.650000000000000E+02, /* K = 7 */
1224     (PID.TID 0000.0001) 3.150000000000000E+02, /* K = 8 */
1225     (PID.TID 0000.0001) 3.650000000000000E+02, /* K = 9 */
1226     (PID.TID 0000.0001) 4.150000000000000E+02, /* K = 10 */
1227     (PID.TID 0000.0001) 4.650000000000000E+02, /* K = 11 */
1228     (PID.TID 0000.0001) 5.150000000000000E+02, /* K = 12 */
1229     (PID.TID 0000.0001) 5.650000000000000E+02, /* K = 13 */
1230     (PID.TID 0000.0001) 6.150000000000000E+02, /* K = 14 */
1231 jmc 1.5 (PID.TID 0000.0001) 6.650000000000000E+02, /* K = 15 */
1232     (PID.TID 0000.0001) 3.450000000000000E+02 /* K = 16 */
1233 jmc 1.1 (PID.TID 0000.0001) ;
1234     (PID.TID 0000.0001) drF = /* W spacing ( units of r ) */
1235     (PID.TID 0000.0001) 5.000000000000000E+01, /* K = 1 */
1236     (PID.TID 0000.0001) 7.000000000000000E+01, /* K = 2 */
1237     (PID.TID 0000.0001) 1.000000000000000E+02, /* K = 3 */
1238     (PID.TID 0000.0001) 1.400000000000000E+02, /* K = 4 */
1239     (PID.TID 0000.0001) 1.900000000000000E+02, /* K = 5 */
1240     (PID.TID 0000.0001) 2.400000000000000E+02, /* K = 6 */
1241     (PID.TID 0000.0001) 2.900000000000000E+02, /* K = 7 */
1242     (PID.TID 0000.0001) 3.400000000000000E+02, /* K = 8 */
1243     (PID.TID 0000.0001) 3.900000000000000E+02, /* K = 9 */
1244     (PID.TID 0000.0001) 4.400000000000000E+02, /* K = 10 */
1245     (PID.TID 0000.0001) 4.900000000000000E+02, /* K = 11 */
1246     (PID.TID 0000.0001) 5.400000000000000E+02, /* K = 12 */
1247     (PID.TID 0000.0001) 5.900000000000000E+02, /* K = 13 */
1248     (PID.TID 0000.0001) 6.400000000000000E+02, /* K = 14 */
1249     (PID.TID 0000.0001) 6.900000000000000E+02 /* K = 15 */
1250     (PID.TID 0000.0001) ;
1251 jmc 1.3 (PID.TID 0000.0001) delX = /* U spacing ( m - cartesian, degrees - spherical ) */
1252 jmc 1.1 (PID.TID 0000.0001) 90 @ 4.000000000000000E+00 /* I = 1: 90 */
1253     (PID.TID 0000.0001) ;
1254 jmc 1.3 (PID.TID 0000.0001) delY = /* V spacing ( m - cartesian, degrees - spherical ) */
1255 jmc 1.1 (PID.TID 0000.0001) 40 @ 4.000000000000000E+00 /* J = 1: 40 */
1256     (PID.TID 0000.0001) ;
1257 jmc 1.5 (PID.TID 0000.0001) xgOrigin = /* X-axis origin of West edge (cartesian: m, lat-lon: deg) */
1258 jmc 1.3 (PID.TID 0000.0001) 0.000000000000000E+00
1259     (PID.TID 0000.0001) ;
1260 jmc 1.5 (PID.TID 0000.0001) ygOrigin = /* Y-axis origin of South edge (cartesian: m, lat-lon: deg) */
1261 jmc 1.1 (PID.TID 0000.0001) -8.000000000000000E+01
1262     (PID.TID 0000.0001) ;
1263     (PID.TID 0000.0001) rSphere = /* Radius ( ignored - cartesian, m - spherical ) */
1264     (PID.TID 0000.0001) 6.370000000000000E+06
1265     (PID.TID 0000.0001) ;
1266     (PID.TID 0000.0001) deepAtmosphere = /* Deep/Shallow Atmosphere flag (True/False) */
1267     (PID.TID 0000.0001) F
1268     (PID.TID 0000.0001) ;
1269 jmc 1.3 (PID.TID 0000.0001) xC = /* xC(:,1,:,1) : P-point X coord ( deg. or m if cartesian) */
1270 jmc 1.1 (PID.TID 0000.0001) 2.000000000000000E+00, /* I = 1 */
1271     (PID.TID 0000.0001) 6.000000000000000E+00, /* I = 2 */
1272     (PID.TID 0000.0001) 1.000000000000000E+01, /* I = 3 */
1273 jmc 1.3 (PID.TID 0000.0001) . . .
1274 jmc 1.1 (PID.TID 0000.0001) 8.200000000000000E+01, /* I = 21 */
1275     (PID.TID 0000.0001) 8.600000000000000E+01, /* I = 22 */
1276     (PID.TID 0000.0001) 9.000000000000000E+01, /* I = 23 */
1277     (PID.TID 0000.0001) 9.400000000000000E+01, /* I = 24 */
1278     (PID.TID 0000.0001) 9.800000000000000E+01, /* I = 25 */
1279     (PID.TID 0000.0001) 1.020000000000000E+02, /* I = 26 */
1280 jmc 1.3 (PID.TID 0000.0001) . . .
1281 jmc 1.1 (PID.TID 0000.0001) 1.700000000000000E+02, /* I = 43 */
1282     (PID.TID 0000.0001) 1.740000000000000E+02, /* I = 44 */
1283     (PID.TID 0000.0001) 1.780000000000000E+02, /* I = 45 */
1284     (PID.TID 0000.0001) 1.820000000000000E+02, /* I = 46 */
1285     (PID.TID 0000.0001) 1.860000000000000E+02, /* I = 47 */
1286     (PID.TID 0000.0001) 1.900000000000000E+02, /* I = 48 */
1287 jmc 1.3 (PID.TID 0000.0001) . . .
1288 jmc 1.1 (PID.TID 0000.0001) 2.580000000000000E+02, /* I = 65 */
1289     (PID.TID 0000.0001) 2.620000000000000E+02, /* I = 66 */
1290     (PID.TID 0000.0001) 2.660000000000000E+02, /* I = 67 */
1291     (PID.TID 0000.0001) 2.700000000000000E+02, /* I = 68 */
1292     (PID.TID 0000.0001) 2.740000000000000E+02, /* I = 69 */
1293     (PID.TID 0000.0001) 2.780000000000000E+02, /* I = 70 */
1294 jmc 1.3 (PID.TID 0000.0001) . . .
1295 jmc 1.1 (PID.TID 0000.0001) 3.500000000000000E+02, /* I = 88 */
1296     (PID.TID 0000.0001) 3.540000000000000E+02, /* I = 89 */
1297     (PID.TID 0000.0001) 3.580000000000000E+02 /* I = 90 */
1298     (PID.TID 0000.0001) ;
1299 jmc 1.3 (PID.TID 0000.0001) yC = /* yC(1,:,1,:) : P-point Y coord ( deg. or m if cartesian) */
1300 jmc 1.1 (PID.TID 0000.0001) -7.800000000000000E+01, /* J = 1 */
1301     (PID.TID 0000.0001) -7.400000000000000E+01, /* J = 2 */
1302     (PID.TID 0000.0001) -7.000000000000000E+01, /* J = 3 */
1303     (PID.TID 0000.0001) -6.600000000000000E+01, /* J = 4 */
1304     (PID.TID 0000.0001) -6.200000000000000E+01, /* J = 5 */
1305     (PID.TID 0000.0001) -5.800000000000000E+01, /* J = 6 */
1306     (PID.TID 0000.0001) -5.400000000000000E+01, /* J = 7 */
1307     (PID.TID 0000.0001) -5.000000000000000E+01, /* J = 8 */
1308     (PID.TID 0000.0001) -4.600000000000000E+01, /* J = 9 */
1309     (PID.TID 0000.0001) -4.200000000000000E+01, /* J = 10 */
1310     (PID.TID 0000.0001) -3.800000000000000E+01, /* J = 11 */
1311     (PID.TID 0000.0001) -3.400000000000000E+01, /* J = 12 */
1312     (PID.TID 0000.0001) -3.000000000000000E+01, /* J = 13 */
1313     (PID.TID 0000.0001) -2.600000000000000E+01, /* J = 14 */
1314     (PID.TID 0000.0001) -2.200000000000000E+01, /* J = 15 */
1315     (PID.TID 0000.0001) -1.800000000000000E+01, /* J = 16 */
1316     (PID.TID 0000.0001) -1.400000000000000E+01, /* J = 17 */
1317     (PID.TID 0000.0001) -1.000000000000000E+01, /* J = 18 */
1318     (PID.TID 0000.0001) -6.000000000000000E+00, /* J = 19 */
1319     (PID.TID 0000.0001) -2.000000000000000E+00, /* J = 20 */
1320     (PID.TID 0000.0001) 2.000000000000000E+00, /* J = 21 */
1321     (PID.TID 0000.0001) 6.000000000000000E+00, /* J = 22 */
1322     (PID.TID 0000.0001) 1.000000000000000E+01, /* J = 23 */
1323     (PID.TID 0000.0001) 1.400000000000000E+01, /* J = 24 */
1324     (PID.TID 0000.0001) 1.800000000000000E+01, /* J = 25 */
1325     (PID.TID 0000.0001) 2.200000000000000E+01, /* J = 26 */
1326     (PID.TID 0000.0001) 2.600000000000000E+01, /* J = 27 */
1327     (PID.TID 0000.0001) 3.000000000000000E+01, /* J = 28 */
1328     (PID.TID 0000.0001) 3.400000000000000E+01, /* J = 29 */
1329     (PID.TID 0000.0001) 3.800000000000000E+01, /* J = 30 */
1330     (PID.TID 0000.0001) 4.200000000000000E+01, /* J = 31 */
1331     (PID.TID 0000.0001) 4.600000000000000E+01, /* J = 32 */
1332     (PID.TID 0000.0001) 5.000000000000000E+01, /* J = 33 */
1333     (PID.TID 0000.0001) 5.400000000000000E+01, /* J = 34 */
1334     (PID.TID 0000.0001) 5.800000000000000E+01, /* J = 35 */
1335     (PID.TID 0000.0001) 6.200000000000000E+01, /* J = 36 */
1336     (PID.TID 0000.0001) 6.600000000000000E+01, /* J = 37 */
1337     (PID.TID 0000.0001) 7.000000000000000E+01, /* J = 38 */
1338     (PID.TID 0000.0001) 7.400000000000000E+01, /* J = 39 */
1339     (PID.TID 0000.0001) 7.800000000000000E+01 /* J = 40 */
1340     (PID.TID 0000.0001) ;
1341 jmc 1.3 (PID.TID 0000.0001) rcoord = /* P-point R coordinate ( units of r ) */
1342 jmc 1.1 (PID.TID 0000.0001) -2.500000000000000E+01, /* K = 1 */
1343     (PID.TID 0000.0001) -8.500000000000000E+01, /* K = 2 */
1344     (PID.TID 0000.0001) -1.700000000000000E+02, /* K = 3 */
1345     (PID.TID 0000.0001) -2.900000000000000E+02, /* K = 4 */
1346     (PID.TID 0000.0001) -4.550000000000000E+02, /* K = 5 */
1347     (PID.TID 0000.0001) -6.700000000000000E+02, /* K = 6 */
1348     (PID.TID 0000.0001) -9.350000000000000E+02, /* K = 7 */
1349     (PID.TID 0000.0001) -1.250000000000000E+03, /* K = 8 */
1350     (PID.TID 0000.0001) -1.615000000000000E+03, /* K = 9 */
1351     (PID.TID 0000.0001) -2.030000000000000E+03, /* K = 10 */
1352     (PID.TID 0000.0001) -2.495000000000000E+03, /* K = 11 */
1353     (PID.TID 0000.0001) -3.010000000000000E+03, /* K = 12 */
1354     (PID.TID 0000.0001) -3.575000000000000E+03, /* K = 13 */
1355     (PID.TID 0000.0001) -4.190000000000000E+03, /* K = 14 */
1356     (PID.TID 0000.0001) -4.855000000000000E+03 /* K = 15 */
1357     (PID.TID 0000.0001) ;
1358     (PID.TID 0000.0001) rF = /* W-Interf. R coordinate ( units of r ) */
1359     (PID.TID 0000.0001) 0.000000000000000E+00, /* K = 1 */
1360     (PID.TID 0000.0001) -5.000000000000000E+01, /* K = 2 */
1361     (PID.TID 0000.0001) -1.200000000000000E+02, /* K = 3 */
1362     (PID.TID 0000.0001) -2.200000000000000E+02, /* K = 4 */
1363     (PID.TID 0000.0001) -3.600000000000000E+02, /* K = 5 */
1364     (PID.TID 0000.0001) -5.500000000000000E+02, /* K = 6 */
1365     (PID.TID 0000.0001) -7.900000000000000E+02, /* K = 7 */
1366     (PID.TID 0000.0001) -1.080000000000000E+03, /* K = 8 */
1367     (PID.TID 0000.0001) -1.420000000000000E+03, /* K = 9 */
1368     (PID.TID 0000.0001) -1.810000000000000E+03, /* K = 10 */
1369     (PID.TID 0000.0001) -2.250000000000000E+03, /* K = 11 */
1370     (PID.TID 0000.0001) -2.740000000000000E+03, /* K = 12 */
1371     (PID.TID 0000.0001) -3.280000000000000E+03, /* K = 13 */
1372     (PID.TID 0000.0001) -3.870000000000000E+03, /* K = 14 */
1373     (PID.TID 0000.0001) -4.510000000000000E+03, /* K = 15 */
1374     (PID.TID 0000.0001) -5.200000000000000E+03 /* K = 16 */
1375     (PID.TID 0000.0001) ;
1376     (PID.TID 0000.0001) deepFacC = /* deep-model grid factor @ cell-Center (-) */
1377     (PID.TID 0000.0001) 15 @ 1.000000000000000E+00 /* K = 1: 15 */
1378     (PID.TID 0000.0001) ;
1379     (PID.TID 0000.0001) deepFacF = /* deep-model grid factor @ W-Interface (-) */
1380     (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
1381     (PID.TID 0000.0001) ;
1382     (PID.TID 0000.0001) rVel2wUnit = /* convert units: rVel -> wSpeed (=1 if z-coord)*/
1383     (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
1384     (PID.TID 0000.0001) ;
1385     (PID.TID 0000.0001) wUnit2rVel = /* convert units: wSpeed -> rVel (=1 if z-coord)*/
1386     (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
1387     (PID.TID 0000.0001) ;
1388 jmc 1.3 (PID.TID 0000.0001) dBdrRef = /* Vertical grad. of reference buoyancy [(m/s/r)^2] */
1389 jmc 1.1 (PID.TID 0000.0001) 15 @ 0.000000000000000E+00 /* K = 1: 15 */
1390     (PID.TID 0000.0001) ;
1391 jmc 1.3 (PID.TID 0000.0001) rotateGrid = /* use rotated grid ( True/False ) */
1392     (PID.TID 0000.0001) F
1393     (PID.TID 0000.0001) ;
1394     (PID.TID 0000.0001) phiEuler = /* Euler angle, rotation about original z-coordinate [rad] */
1395     (PID.TID 0000.0001) 0.000000000000000E+00
1396     (PID.TID 0000.0001) ;
1397     (PID.TID 0000.0001) thetaEuler = /* Euler angle, rotation about new x-coordinate [rad] */
1398     (PID.TID 0000.0001) 0.000000000000000E+00
1399     (PID.TID 0000.0001) ;
1400     (PID.TID 0000.0001) psiEuler = /* Euler angle, rotation about new z-coordinate [rad] */
1401     (PID.TID 0000.0001) 0.000000000000000E+00
1402     (PID.TID 0000.0001) ;
1403 jmc 1.1 (PID.TID 0000.0001) dxF = /* dxF(:,1,:,1) ( units: m ) */
1404     (PID.TID 0000.0001) 90 @ 9.246038586187513E+04 /* I = 1: 90 */
1405     (PID.TID 0000.0001) ;
1406     (PID.TID 0000.0001) dxF = /* dxF(1,:,1,:) ( units: m ) */
1407     (PID.TID 0000.0001) 9.246038586187513E+04, /* J = 1 */
1408     (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 2 */
1409     (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 3 */
1410     (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 4 */
1411     (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 5 */
1412     (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 6 */
1413     (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 7 */
1414     (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 8 */
1415     (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 9 */
1416     (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 10 */
1417     (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 11 */
1418     (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 12 */
1419     (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 13 */
1420     (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 14 */
1421     (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 15 */
1422     (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 16 */
1423     (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 17 */
1424     (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 18 */
1425     (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 19 */
1426     (PID.TID 0000.0001) 2 @ 4.444389881567502E+05, /* J = 20: 21 */
1427     (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 22 */
1428     (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 23 */
1429     (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 24 */
1430     (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 25 */
1431     (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 26 */
1432     (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 27 */
1433     (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 28 */
1434     (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 29 */
1435     (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 30 */
1436     (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 31 */
1437     (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 32 */
1438     (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 33 */
1439     (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 34 */
1440     (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 35 */
1441     (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 36 */
1442     (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 37 */
1443     (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 38 */
1444     (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 39 */
1445     (PID.TID 0000.0001) 9.246038586187513E+04 /* J = 40 */
1446     (PID.TID 0000.0001) ;
1447     (PID.TID 0000.0001) dyF = /* dyF(:,1,:,1) ( units: m ) */
1448     (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1449     (PID.TID 0000.0001) ;
1450     (PID.TID 0000.0001) dyF = /* dyF(1,:,1,:) ( units: m ) */
1451     (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1452     (PID.TID 0000.0001) ;
1453     (PID.TID 0000.0001) dxG = /* dxG(:,1,:,1) ( units: m ) */
1454     (PID.TID 0000.0001) 90 @ 7.722306258078101E+04 /* I = 1: 90 */
1455     (PID.TID 0000.0001) ;
1456     (PID.TID 0000.0001) dxG = /* dxG(1,:,1,:) ( units: m ) */
1457     (PID.TID 0000.0001) 7.722306258078101E+04, /* J = 1 */
1458     (PID.TID 0000.0001) 1.075850604052271E+05, /* J = 2 */
1459     (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 3 */
1460     (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 4 */
1461     (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 5 */
1462     (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 6 */
1463     (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 7 */
1464     (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 8 */
1465     (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 9 */
1466     (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 10 */
1467     (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 11 */
1468     (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 12 */
1469     (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 13 */
1470     (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 14 */
1471     (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 15 */
1472     (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 16 */
1473     (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 17 */
1474     (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 18 */
1475     (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 19 */
1476     (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 20 */
1477     (PID.TID 0000.0001) 4.447098934081552E+05, /* J = 21 */
1478     (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 22 */
1479     (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 23 */
1480     (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 24 */
1481     (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 25 */
1482     (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 26 */
1483     (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 27 */
1484     (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 28 */
1485     (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 29 */
1486     (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 30 */
1487     (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 31 */
1488     (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 32 */
1489     (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 33 */
1490     (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 34 */
1491     (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 35 */
1492     (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 36 */
1493     (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 37 */
1494     (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 38 */
1495     (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 39 */
1496     (PID.TID 0000.0001) 1.075850604052271E+05 /* J = 40 */
1497     (PID.TID 0000.0001) ;
1498     (PID.TID 0000.0001) dyG = /* dyG(:,1,:,1) ( units: m ) */
1499     (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1500     (PID.TID 0000.0001) ;
1501     (PID.TID 0000.0001) dyG = /* dyG(1,:,1,:) ( units: m ) */
1502     (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1503     (PID.TID 0000.0001) ;
1504     (PID.TID 0000.0001) dxC = /* dxC(:,1,:,1) ( units: m ) */
1505     (PID.TID 0000.0001) 90 @ 9.246038586187513E+04 /* I = 1: 90 */
1506     (PID.TID 0000.0001) ;
1507     (PID.TID 0000.0001) dxC = /* dxC(1,:,1,:) ( units: m ) */
1508     (PID.TID 0000.0001) 9.246038586187513E+04, /* J = 1 */
1509     (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 2 */
1510     (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 3 */
1511     (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 4 */
1512     (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 5 */
1513     (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 6 */
1514     (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 7 */
1515     (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 8 */
1516     (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 9 */
1517     (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 10 */
1518     (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 11 */
1519     (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 12 */
1520     (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 13 */
1521     (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 14 */
1522     (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 15 */
1523     (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 16 */
1524     (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 17 */
1525     (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 18 */
1526     (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 19 */
1527     (PID.TID 0000.0001) 2 @ 4.444389881567502E+05, /* J = 20: 21 */
1528     (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 22 */
1529     (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 23 */
1530     (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 24 */
1531     (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 25 */
1532     (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 26 */
1533     (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 27 */
1534     (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 28 */
1535     (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 29 */
1536     (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 30 */
1537     (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 31 */
1538     (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 32 */
1539     (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 33 */
1540     (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 34 */
1541     (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 35 */
1542     (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 36 */
1543     (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 37 */
1544     (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 38 */
1545     (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 39 */
1546     (PID.TID 0000.0001) 9.246038586187513E+04 /* J = 40 */
1547     (PID.TID 0000.0001) ;
1548     (PID.TID 0000.0001) dyC = /* dyC(:,1,:,1) ( units: m ) */
1549     (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1550     (PID.TID 0000.0001) ;
1551     (PID.TID 0000.0001) dyC = /* dyC(1,:,1,:) ( units: m ) */
1552     (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1553     (PID.TID 0000.0001) ;
1554     (PID.TID 0000.0001) dxV = /* dxV(:,1,:,1) ( units: m ) */
1555     (PID.TID 0000.0001) 90 @ 7.722306258078101E+04 /* I = 1: 90 */
1556     (PID.TID 0000.0001) ;
1557     (PID.TID 0000.0001) dxV = /* dxV(1,:,1,:) ( units: m ) */
1558     (PID.TID 0000.0001) 7.722306258078101E+04, /* J = 1 */
1559     (PID.TID 0000.0001) 1.075850604052271E+05, /* J = 2 */
1560     (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 3 */
1561     (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 4 */
1562     (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 5 */
1563     (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 6 */
1564     (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 7 */
1565     (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 8 */
1566     (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 9 */
1567     (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 10 */
1568     (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 11 */
1569     (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 12 */
1570     (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 13 */
1571     (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 14 */
1572     (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 15 */
1573     (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 16 */
1574     (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 17 */
1575     (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 18 */
1576     (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 19 */
1577     (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 20 */
1578     (PID.TID 0000.0001) 4.447098934081552E+05, /* J = 21 */
1579     (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 22 */
1580     (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 23 */
1581     (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 24 */
1582     (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 25 */
1583     (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 26 */
1584     (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 27 */
1585     (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 28 */
1586     (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 29 */
1587     (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 30 */
1588     (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 31 */
1589     (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 32 */
1590     (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 33 */
1591     (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 34 */
1592     (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 35 */
1593     (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 36 */
1594     (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 37 */
1595     (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 38 */
1596     (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 39 */
1597     (PID.TID 0000.0001) 1.075850604052271E+05 /* J = 40 */
1598     (PID.TID 0000.0001) ;
1599     (PID.TID 0000.0001) dyU = /* dyU(:,1,:,1) ( units: m ) */
1600     (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1601     (PID.TID 0000.0001) ;
1602     (PID.TID 0000.0001) dyU = /* dyU(1,:,1,:) ( units: m ) */
1603     (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1604     (PID.TID 0000.0001) ;
1605     (PID.TID 0000.0001) rA = /* rA (:,1,:,1) ( units: m^2 ) */
1606     (PID.TID 0000.0001) 90 @ 4.110969866729047E+10 /* I = 1: 90 */
1607     (PID.TID 0000.0001) ;
1608     (PID.TID 0000.0001) rA = /* rA (1,:,1,:) ( units: m^2 ) */
1609     (PID.TID 0000.0001) 4.110969866729047E+10, /* J = 1 */
1610     (PID.TID 0000.0001) 5.450087291636665E+10, /* J = 2 */
1611     (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 3 */
1612     (PID.TID 0000.0001) 8.042270623659998E+10, /* J = 4 */
1613     (PID.TID 0000.0001) 9.282707674147523E+10, /* J = 5 */
1614     (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 6 */
1615     (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 7 */
1616     (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 8 */
1617     (PID.TID 0000.0001) 1.373525277677229E+11, /* J = 9 */
1618     (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 10 */
1619     (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 11 */
1620     (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 12 */
1621     (PID.TID 0000.0001) 1.712363708253573E+11, /* J = 13 */
1622     (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 14 */
1623     (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 15 */
1624     (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 16 */
1625     (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 17 */
1626     (PID.TID 0000.0001) 1.947228162702492E+11, /* J = 18 */
1627     (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 19 */
1628     (PID.TID 0000.0001) 2 @ 1.976062798008856E+11, /* J = 20: 21 */
1629     (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 22 */
1630     (PID.TID 0000.0001) 1.947228162702492E+11, /* J = 23 */
1631     (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 24 */
1632     (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 25 */
1633     (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 26 */
1634     (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 27 */
1635     (PID.TID 0000.0001) 1.712363708253573E+11, /* J = 28 */
1636     (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 29 */
1637     (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 30 */
1638     (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 31 */
1639     (PID.TID 0000.0001) 1.373525277677229E+11, /* J = 32 */
1640     (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 33 */
1641     (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 34 */
1642     (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 35 */
1643     (PID.TID 0000.0001) 9.282707674147523E+10, /* J = 36 */
1644     (PID.TID 0000.0001) 8.042270623659998E+10, /* J = 37 */
1645     (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 38 */
1646     (PID.TID 0000.0001) 5.450087291636665E+10, /* J = 39 */
1647     (PID.TID 0000.0001) 4.110969866729047E+10 /* J = 40 */
1648     (PID.TID 0000.0001) ;
1649     (PID.TID 0000.0001) rAw = /* rAw(:,1,:,1) ( units: m^2 ) */
1650     (PID.TID 0000.0001) 90 @ 4.110969866729047E+10 /* I = 1: 90 */
1651     (PID.TID 0000.0001) ;
1652     (PID.TID 0000.0001) rAw = /* rAw(1,:,1,:) ( units: m^2 ) */
1653     (PID.TID 0000.0001) 4.110969866729047E+10, /* J = 1 */
1654     (PID.TID 0000.0001) 5.450087291636665E+10, /* J = 2 */
1655     (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 3 */
1656     (PID.TID 0000.0001) 8.042270623659998E+10, /* J = 4 */
1657     (PID.TID 0000.0001) 9.282707674147523E+10, /* J = 5 */
1658     (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 6 */
1659     (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 7 */
1660     (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 8 */
1661     (PID.TID 0000.0001) 1.373525277677229E+11, /* J = 9 */
1662     (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 10 */
1663     (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 11 */
1664     (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 12 */
1665     (PID.TID 0000.0001) 1.712363708253573E+11, /* J = 13 */
1666     (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 14 */
1667     (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 15 */
1668     (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 16 */
1669     (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 17 */
1670     (PID.TID 0000.0001) 1.947228162702492E+11, /* J = 18 */
1671     (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 19 */
1672     (PID.TID 0000.0001) 2 @ 1.976062798008856E+11, /* J = 20: 21 */
1673     (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 22 */
1674     (PID.TID 0000.0001) 1.947228162702492E+11, /* J = 23 */
1675     (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 24 */
1676     (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 25 */
1677     (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 26 */
1678     (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 27 */
1679     (PID.TID 0000.0001) 1.712363708253573E+11, /* J = 28 */
1680     (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 29 */
1681     (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 30 */
1682     (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 31 */
1683     (PID.TID 0000.0001) 1.373525277677229E+11, /* J = 32 */
1684     (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 33 */
1685     (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 34 */
1686     (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 35 */
1687     (PID.TID 0000.0001) 9.282707674147523E+10, /* J = 36 */
1688     (PID.TID 0000.0001) 8.042270623659998E+10, /* J = 37 */
1689     (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 38 */
1690     (PID.TID 0000.0001) 5.450087291636665E+10, /* J = 39 */
1691     (PID.TID 0000.0001) 4.110969866729047E+10 /* J = 40 */
1692     (PID.TID 0000.0001) ;
1693     (PID.TID 0000.0001) rAs = /* rAs(:,1,:,1) ( units: m^2 ) */
1694     (PID.TID 0000.0001) 90 @ 3.433488626798250E+10 /* I = 1: 90 */
1695     (PID.TID 0000.0001) ;
1696     (PID.TID 0000.0001) rAs = /* rAs(1,:,1,:) ( units: m^2 ) */
1697     (PID.TID 0000.0001) 3.433488626798250E+10, /* J = 1 */
1698     (PID.TID 0000.0001) 4.783442523123625E+10, /* J = 2 */
1699     (PID.TID 0000.0001) 6.110091968306414E+10, /* J = 3 */
1700     (PID.TID 0000.0001) 7.406973659603818E+10, /* J = 4 */
1701     (PID.TID 0000.0001) 8.667769319778079E+10, /* J = 5 */
1702     (PID.TID 0000.0001) 9.886336479107463E+10, /* J = 6 */
1703     (PID.TID 0000.0001) 1.105673840088173E+11, /* J = 7 */
1704     (PID.TID 0000.0001) 1.217327300458638E+11, /* J = 8 */
1705     (PID.TID 0000.0001) 1.323050064586578E+11, /* J = 9 */
1706     (PID.TID 0000.0001) 1.422327061792377E+11, /* J = 10 */
1707     (PID.TID 0000.0001) 1.514674624524945E+11, /* J = 11 */
1708     (PID.TID 0000.0001) 1.599642844741385E+11, /* J = 12 */
1709     (PID.TID 0000.0001) 1.676817765813788E+11, /* J = 13 */
1710     (PID.TID 0000.0001) 1.745823399284268E+11, /* J = 14 */
1711     (PID.TID 0000.0001) 1.806323556642996E+11, /* J = 15 */
1712     (PID.TID 0000.0001) 1.858023487204767E+11, /* J = 16 */
1713     (PID.TID 0000.0001) 1.900671314104743E+11, /* J = 17 */
1714     (PID.TID 0000.0001) 1.934059261417216E+11, /* J = 18 */
1715     (PID.TID 0000.0001) 1.958024666419019E+11, /* J = 19 */
1716     (PID.TID 0000.0001) 1.972450772065981E+11, /* J = 20 */
1717     (PID.TID 0000.0001) 1.977267295821495E+11, /* J = 21 */
1718     (PID.TID 0000.0001) 1.972450772065981E+11, /* J = 22 */
1719     (PID.TID 0000.0001) 1.958024666419019E+11, /* J = 23 */
1720     (PID.TID 0000.0001) 1.934059261417216E+11, /* J = 24 */
1721     (PID.TID 0000.0001) 1.900671314104743E+11, /* J = 25 */
1722     (PID.TID 0000.0001) 1.858023487204767E+11, /* J = 26 */
1723     (PID.TID 0000.0001) 1.806323556642996E+11, /* J = 27 */
1724     (PID.TID 0000.0001) 1.745823399284268E+11, /* J = 28 */
1725     (PID.TID 0000.0001) 1.676817765813788E+11, /* J = 29 */
1726     (PID.TID 0000.0001) 1.599642844741385E+11, /* J = 30 */
1727     (PID.TID 0000.0001) 1.514674624524945E+11, /* J = 31 */
1728     (PID.TID 0000.0001) 1.422327061792377E+11, /* J = 32 */
1729     (PID.TID 0000.0001) 1.323050064586578E+11, /* J = 33 */
1730     (PID.TID 0000.0001) 1.217327300458638E+11, /* J = 34 */
1731     (PID.TID 0000.0001) 1.105673840088173E+11, /* J = 35 */
1732     (PID.TID 0000.0001) 9.886336479107463E+10, /* J = 36 */
1733     (PID.TID 0000.0001) 8.667769319778079E+10, /* J = 37 */
1734     (PID.TID 0000.0001) 7.406973659603818E+10, /* J = 38 */
1735     (PID.TID 0000.0001) 6.110091968306414E+10, /* J = 39 */
1736     (PID.TID 0000.0001) 4.783442523123625E+10 /* J = 40 */
1737     (PID.TID 0000.0001) ;
1738     (PID.TID 0000.0001) globalArea = /* Integrated horizontal Area (m^2) */
1739     (PID.TID 0000.0001) 3.450614146649838E+14
1740     (PID.TID 0000.0001) ;
1741     (PID.TID 0000.0001) // =======================================================
1742     (PID.TID 0000.0001) // End of Model config. summary
1743     (PID.TID 0000.0001) // =======================================================
1744     (PID.TID 0000.0001)
1745 jmc 1.3 (PID.TID 0000.0001) == Packages configuration : Check & print summary ==
1746     (PID.TID 0000.0001)
1747     (PID.TID 0000.0001) GMREDI_CHECK: #define GMREDI
1748     (PID.TID 0000.0001) GM_AdvForm = /* if FALSE => use SkewFlux Form */
1749     (PID.TID 0000.0001) F
1750     (PID.TID 0000.0001) ;
1751     (PID.TID 0000.0001) GM_InMomAsStress = /* if TRUE => apply as Eddy Stress */
1752     (PID.TID 0000.0001) F
1753     (PID.TID 0000.0001) ;
1754     (PID.TID 0000.0001) GM_AdvSeparate = /* Calc Bolus & Euler Adv. separately */
1755     (PID.TID 0000.0001) F
1756     (PID.TID 0000.0001) ;
1757     (PID.TID 0000.0001) GM_ExtraDiag = /* Tensor Extra Diag (line 1&2) non 0 */
1758     (PID.TID 0000.0001) F
1759     (PID.TID 0000.0001) ;
1760     (PID.TID 0000.0001) GM_isopycK = /* Background Isopyc. Diffusivity [m^2/s] */
1761     (PID.TID 0000.0001) 1.000000000000000E+03
1762     (PID.TID 0000.0001) ;
1763     (PID.TID 0000.0001) GM_skewflx*K = /* Background GM_SkewFlx Diffusivity [m^2/s] */
1764     (PID.TID 0000.0001) 1.000000000000000E+03
1765     (PID.TID 0000.0001) ;
1766     (PID.TID 0000.0001) GM_advec*K = /* Backg. GM-Advec(=Bolus) Diffusivity [m^2/s]*/
1767     (PID.TID 0000.0001) 0.000000000000000E+00
1768     (PID.TID 0000.0001) ;
1769     (PID.TID 0000.0001) GM_Kmin_horiz = /* Minimum Horizontal Diffusivity [m^2/s] */
1770     (PID.TID 0000.0001) 5.000000000000000E+01
1771     (PID.TID 0000.0001) ;
1772     (PID.TID 0000.0001) GM_Visbeck_alpha = /* Visbeck alpha coeff. [-] */
1773     (PID.TID 0000.0001) 0.000000000000000E+00
1774     (PID.TID 0000.0001) ;
1775     (PID.TID 0000.0001) GM_Small_Number = /* epsilon used in slope calc */
1776     (PID.TID 0000.0001) 1.000000000000000E-12
1777     (PID.TID 0000.0001) ;
1778     (PID.TID 0000.0001) GM_slopeSqCutoff = /* Slope^2 cut-off value */
1779     (PID.TID 0000.0001) 1.000000000000000E+08
1780     (PID.TID 0000.0001) ;
1781     (PID.TID 0000.0001) GM_taper_scheme = /* Type of Tapering/Clipping scheme */
1782     (PID.TID 0000.0001) 'dm95 '
1783     (PID.TID 0000.0001) ;
1784     (PID.TID 0000.0001) GM_maxSlope = /* Maximum Slope (Tapering/Clipping) */
1785     (PID.TID 0000.0001) 1.000000000000000E-02
1786     (PID.TID 0000.0001) ;
1787     (PID.TID 0000.0001) GM_facTrL2dz = /* Minimum Trans.Layer Thick. (factor of dz) */
1788     (PID.TID 0000.0001) 1.000000000000000E+00
1789     (PID.TID 0000.0001) ;
1790     (PID.TID 0000.0001) GM_facTrL2ML = /* Max.Trans.Layer Thick. (factor of MxL Depth)*/
1791     (PID.TID 0000.0001) 5.000000000000000E+00
1792     (PID.TID 0000.0001) ;
1793     (PID.TID 0000.0001) GM_maxTransLay = /* Maximum Transition Layer Thickness [m] */
1794     (PID.TID 0000.0001) 5.000000000000000E+02
1795     (PID.TID 0000.0001) ;
1796     (PID.TID 0000.0001) GM_UseBVP = /* if TRUE => use bvp a la Ferrari et al. (2010) */
1797     (PID.TID 0000.0001) F
1798     (PID.TID 0000.0001) ;
1799     (PID.TID 0000.0001) GM_BVP_ModeNumber = /* Vertical mode number for BVP wave speed */
1800     (PID.TID 0000.0001) 1
1801     (PID.TID 0000.0001) ;
1802     (PID.TID 0000.0001) GM_BVP_cMin = /* Minimum wave speed for BVP [m/s] */
1803     (PID.TID 0000.0001) 1.000000000000000E-01
1804     (PID.TID 0000.0001) ;
1805 jmc 1.4 (PID.TID 0000.0001) GM_useSubMeso = /* if TRUE => use Sub-Meso param. (B.Fox-Kemper) */
1806     (PID.TID 0000.0001) F
1807     (PID.TID 0000.0001) ;
1808     (PID.TID 0000.0001) subMeso_Ceff = /* efficiency coeff. of Mixed-Layer Eddies [-] */
1809     (PID.TID 0000.0001) 7.000000000000001E-02
1810     (PID.TID 0000.0001) ;
1811     (PID.TID 0000.0001) subMeso_invTau = /* inverse of Sub-Meso mixing time-scale [/s] */
1812     (PID.TID 0000.0001) 2.000000000000000E-06
1813     (PID.TID 0000.0001) ;
1814     (PID.TID 0000.0001) subMeso_LfMin = /* minimum length-scale "Lf" [m] */
1815     (PID.TID 0000.0001) 1.000000000000000E+03
1816     (PID.TID 0000.0001) ;
1817     (PID.TID 0000.0001) subMeso_Lmax = /* maximum grid-scale length [m] */
1818     (PID.TID 0000.0001) 1.100000000000000E+05
1819     (PID.TID 0000.0001) ;
1820 jmc 1.3 (PID.TID 0000.0001) CTRL_CHECK: ctrl package
1821     (PID.TID 0000.0001) COST_CHECK: cost package
1822     (PID.TID 0000.0001) GRDCHK_CHECK: grdchk package
1823     (PID.TID 0000.0001) GAD_CHECK: #define ALLOW_GENERIC_ADVDIFF
1824 jmc 1.1 (PID.TID 0000.0001) // =======================================================
1825 jmc 1.5 (PID.TID 0000.0001) // Check Model config. (CONFIG_CHECK):
1826 jmc 1.3 (PID.TID 0000.0001) // CONFIG_CHECK : Normal End
1827 jmc 1.1 (PID.TID 0000.0001) // =======================================================
1828     (PID.TID 0000.0001)
1829     (PID.TID 0000.0001) Start initial hydrostatic pressure computation
1830     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
1831     (PID.TID 0000.0001)
1832     (PID.TID 0000.0001) // =======================================================
1833     (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
1834     (PID.TID 0000.0001) // =======================================================
1835     (PID.TID 0000.0001) %MON time_tsnumber = 0
1836     (PID.TID 0000.0001) %MON time_secondsf = 0.0000000000000E+00
1837     (PID.TID 0000.0001) %MON dynstat_eta_max = 0.0000000000000E+00
1838     (PID.TID 0000.0001) %MON dynstat_eta_min = 0.0000000000000E+00
1839     (PID.TID 0000.0001) %MON dynstat_eta_mean = 0.0000000000000E+00
1840     (PID.TID 0000.0001) %MON dynstat_eta_sd = 0.0000000000000E+00
1841     (PID.TID 0000.0001) %MON dynstat_eta_del2 = 0.0000000000000E+00
1842     (PID.TID 0000.0001) %MON dynstat_uvel_max = 0.0000000000000E+00
1843     (PID.TID 0000.0001) %MON dynstat_uvel_min = 0.0000000000000E+00
1844     (PID.TID 0000.0001) %MON dynstat_uvel_mean = 0.0000000000000E+00
1845     (PID.TID 0000.0001) %MON dynstat_uvel_sd = 0.0000000000000E+00
1846     (PID.TID 0000.0001) %MON dynstat_uvel_del2 = 0.0000000000000E+00
1847     (PID.TID 0000.0001) %MON dynstat_vvel_max = 0.0000000000000E+00
1848     (PID.TID 0000.0001) %MON dynstat_vvel_min = 0.0000000000000E+00
1849     (PID.TID 0000.0001) %MON dynstat_vvel_mean = 0.0000000000000E+00
1850     (PID.TID 0000.0001) %MON dynstat_vvel_sd = 0.0000000000000E+00
1851     (PID.TID 0000.0001) %MON dynstat_vvel_del2 = 0.0000000000000E+00
1852 jmc 1.3 (PID.TID 0000.0001) %MON dynstat_wvel_max = -0.0000000000000E+00
1853     (PID.TID 0000.0001) %MON dynstat_wvel_min = -0.0000000000000E+00
1854 jmc 1.1 (PID.TID 0000.0001) %MON dynstat_wvel_mean = 0.0000000000000E+00
1855     (PID.TID 0000.0001) %MON dynstat_wvel_sd = 0.0000000000000E+00
1856     (PID.TID 0000.0001) %MON dynstat_wvel_del2 = 0.0000000000000E+00
1857     (PID.TID 0000.0001) %MON dynstat_theta_max = 2.9733388900757E+01
1858     (PID.TID 0000.0001) %MON dynstat_theta_min = -1.9000000000000E+00
1859     (PID.TID 0000.0001) %MON dynstat_theta_mean = 3.6197074666129E+00
1860     (PID.TID 0000.0001) %MON dynstat_theta_sd = 4.4174614518568E+00
1861 jmc 1.3 (PID.TID 0000.0001) %MON dynstat_theta_del2 = 5.6229722638211E-03
1862 jmc 1.1 (PID.TID 0000.0001) %MON dynstat_salt_max = 3.7475627899170E+01
1863     (PID.TID 0000.0001) %MON dynstat_salt_min = 2.9752769470215E+01
1864     (PID.TID 0000.0001) %MON dynstat_salt_mean = 3.4718006161475E+01
1865     (PID.TID 0000.0001) %MON dynstat_salt_sd = 2.9628819587990E-01
1866 jmc 1.3 (PID.TID 0000.0001) %MON dynstat_salt_del2 = 1.0127003139695E-03
1867 jmc 1.1 (PID.TID 0000.0001) %MON dynstat_sst_max = 2.9733388900757E+01
1868     (PID.TID 0000.0001) %MON dynstat_sst_min = -1.9000000000000E+00
1869     (PID.TID 0000.0001) %MON dynstat_sst_mean = 1.8484768169451E+01
1870     (PID.TID 0000.0001) %MON dynstat_sst_sd = 9.1795873467430E+00
1871 jmc 1.3 (PID.TID 0000.0001) %MON dynstat_sst_del2 = 2.7217991810411E-02
1872 jmc 1.1 (PID.TID 0000.0001) %MON dynstat_sss_max = 3.7475627899170E+01
1873     (PID.TID 0000.0001) %MON dynstat_sss_min = 2.9752769470215E+01
1874     (PID.TID 0000.0001) %MON dynstat_sss_mean = 3.4849925275207E+01
1875     (PID.TID 0000.0001) %MON dynstat_sss_sd = 9.6344930396289E-01
1876 jmc 1.3 (PID.TID 0000.0001) %MON dynstat_sss_del2 = 8.0490524880020E-03
1877     (PID.TID 0000.0001) %MON forcing_qnet_max = 4.4161361694336E+02
1878     (PID.TID 0000.0001) %MON forcing_qnet_min = -1.8132531738281E+02
1879     (PID.TID 0000.0001) %MON forcing_qnet_mean = -1.6512526001044E+01
1880     (PID.TID 0000.0001) %MON forcing_qnet_sd = 1.0987288293460E+02
1881     (PID.TID 0000.0001) %MON forcing_qnet_del2 = 1.0245701543512E+00
1882     (PID.TID 0000.0001) %MON forcing_qsw_max = 0.0000000000000E+00
1883     (PID.TID 0000.0001) %MON forcing_qsw_min = 0.0000000000000E+00
1884     (PID.TID 0000.0001) %MON forcing_qsw_mean = 0.0000000000000E+00
1885     (PID.TID 0000.0001) %MON forcing_qsw_sd = 0.0000000000000E+00
1886     (PID.TID 0000.0001) %MON forcing_qsw_del2 = 0.0000000000000E+00
1887     (PID.TID 0000.0001) %MON forcing_empmr_max = 0.0000000000000E+00
1888     (PID.TID 0000.0001) %MON forcing_empmr_min = 0.0000000000000E+00
1889     (PID.TID 0000.0001) %MON forcing_empmr_mean = 0.0000000000000E+00
1890     (PID.TID 0000.0001) %MON forcing_empmr_sd = 0.0000000000000E+00
1891     (PID.TID 0000.0001) %MON forcing_empmr_del2 = 0.0000000000000E+00
1892     (PID.TID 0000.0001) %MON forcing_fu_max = 2.7877888083458E-01
1893     (PID.TID 0000.0001) %MON forcing_fu_min = -1.7088057100773E-01
1894     (PID.TID 0000.0001) %MON forcing_fu_mean = 1.6001735152566E-02
1895     (PID.TID 0000.0001) %MON forcing_fu_sd = 9.1062350670783E-02
1896     (PID.TID 0000.0001) %MON forcing_fu_del2 = 7.2483493833855E-04
1897     (PID.TID 0000.0001) %MON forcing_fv_max = 2.3312132060528E-01
1898     (PID.TID 0000.0001) %MON forcing_fv_min = -1.7813690006733E-01
1899     (PID.TID 0000.0001) %MON forcing_fv_mean = -5.2588977395051E-03
1900     (PID.TID 0000.0001) %MON forcing_fv_sd = 4.8248606070107E-02
1901     (PID.TID 0000.0001) %MON forcing_fv_del2 = 4.9948324922575E-04
1902 jmc 1.1 (PID.TID 0000.0001) %MON advcfl_uvel_max = 0.0000000000000E+00
1903     (PID.TID 0000.0001) %MON advcfl_vvel_max = 0.0000000000000E+00
1904     (PID.TID 0000.0001) %MON advcfl_wvel_max = 0.0000000000000E+00
1905     (PID.TID 0000.0001) %MON advcfl_W_hf_max = 0.0000000000000E+00
1906     (PID.TID 0000.0001) %MON pe_b_mean = 0.0000000000000E+00
1907     (PID.TID 0000.0001) %MON ke_max = 0.0000000000000E+00
1908     (PID.TID 0000.0001) %MON ke_mean = 0.0000000000000E+00
1909     (PID.TID 0000.0001) %MON ke_vol = 1.3226782436723E+18
1910     (PID.TID 0000.0001) %MON vort_r_min = 0.0000000000000E+00
1911     (PID.TID 0000.0001) %MON vort_r_max = 0.0000000000000E+00
1912     (PID.TID 0000.0001) %MON vort_a_mean = -2.5205769728595E-05
1913     (PID.TID 0000.0001) %MON vort_a_sd = 7.3403301351718E-05
1914     (PID.TID 0000.0001) %MON vort_p_mean = -3.1758669979758E-05
1915     (PID.TID 0000.0001) %MON vort_p_sd = 1.3126616755304E-04
1916     (PID.TID 0000.0001) %MON surfExpan_theta_mean = 0.0000000000000E+00
1917     (PID.TID 0000.0001) %MON surfExpan_salt_mean = 0.0000000000000E+00
1918     (PID.TID 0000.0001) // =======================================================
1919     (PID.TID 0000.0001) // End MONITOR dynamic field statistics
1920     (PID.TID 0000.0001) // =======================================================
1921 jmc 1.3 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 jmc 1.1 cg2d: Sum(rhs),rhsMax = 1.05956909912663E-14 2.53674886388737E+00
1924 jmc 1.5 cg2d: Sum(rhs),rhsMax = -9.29811783123569E-16 5.04325900480939E+00
1925     cg2d: Sum(rhs),rhsMax = -1.94913529760754E-14 6.29008403763971E+00
1926     cg2d: Sum(rhs),rhsMax = -5.78946612872500E-14 6.70003306073551E+00
1927 jmc 1.3 (PID.TID 0000.0001) %CHECKPOINT 4 ckptA
1928 jmc 1.5 early fc = 0.000000000000000D+00
1929     --> objf_test(bi,bj) = 0.804687744759401D+06
1930 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
1931     --> objf_atl(bi,bj) = 0.000000000000000D+00
1932 jmc 1.5 local fc = 0.804687744759401D+06
1933     global fc = 0.804687744759401D+06
1934 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 1 )
1935     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
1936 jmc 1.1 cg2d: Sum(rhs),rhsMax = 1.05956909912663E-14 2.53674886388737E+00
1937 jmc 1.5 cg2d: Sum(rhs),rhsMax = -9.29811783123569E-16 5.04325900480939E+00
1938     cg2d: Sum(rhs),rhsMax = -1.94913529760754E-14 6.29008403763971E+00
1939     cg2d: Sum(rhs),rhsMax = -5.78946612872500E-14 6.70003306073551E+00
1940 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
1941 jmc 1.5 cg2d: Sum(rhs),rhsMax = -5.78946612872500E-14 6.70003306073551E+00
1942 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
1943     (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
1944 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 0.00000000000000E+00 0.00000000000000E+00
1945 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 1 )
1946     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
1947 jmc 1.1 cg2d: Sum(rhs),rhsMax = 1.05956909912663E-14 2.53674886388737E+00
1948 jmc 1.5 cg2d: Sum(rhs),rhsMax = -9.29811783123569E-16 5.04325900480939E+00
1949     cg2d: Sum(rhs),rhsMax = -1.94913529760754E-14 6.29008403763971E+00
1950     cg2d: Sum(rhs),rhsMax = 1.73456084565093E-15 7.86983807726851E-04
1951     cg2d: Sum(rhs),rhsMax = 1.30787308066926E-14 1.33029393421495E-03
1952     cg2d: Sum(rhs),rhsMax = -1.37259995036665E-15 1.40191055720657E-03
1953 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
1954     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
1955     (PID.TID 0000.0001)
1956     (PID.TID 0000.0001) // =======================================================
1957     (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
1958     (PID.TID 0000.0001) // =======================================================
1959     (PID.TID 0000.0001) %MON ad_time_tsnumber = 0
1960     (PID.TID 0000.0001) %MON ad_time_secondsf = 0.0000000000000E+00
1961 jmc 1.5 (PID.TID 0000.0001) %MON ad_dynstat_adeta_max = 1.3550683205742E+02
1962     (PID.TID 0000.0001) %MON ad_dynstat_adeta_min = -4.3068281185577E+01
1963     (PID.TID 0000.0001) %MON ad_dynstat_adeta_mean = 3.0864811552139E-01
1964     (PID.TID 0000.0001) %MON ad_dynstat_adeta_sd = 1.6951595218077E+01
1965     (PID.TID 0000.0001) %MON ad_dynstat_adeta_del2 = 4.8323265865024E-01
1966     (PID.TID 0000.0001) %MON ad_dynstat_aduvel_max = 5.5954280681864E+01
1967     (PID.TID 0000.0001) %MON ad_dynstat_aduvel_min = -1.0780668343557E+02
1968     (PID.TID 0000.0001) %MON ad_dynstat_aduvel_mean = -3.6363581261081E+00
1969     (PID.TID 0000.0001) %MON ad_dynstat_aduvel_sd = 1.4356156918440E+01
1970     (PID.TID 0000.0001) %MON ad_dynstat_aduvel_del2 = 5.7258676444573E-02
1971     (PID.TID 0000.0001) %MON ad_dynstat_advvel_max = 1.6191225597068E+02
1972     (PID.TID 0000.0001) %MON ad_dynstat_advvel_min = -1.5508690881795E+02
1973     (PID.TID 0000.0001) %MON ad_dynstat_advvel_mean = -3.5089813411514E+00
1974     (PID.TID 0000.0001) %MON ad_dynstat_advvel_sd = 3.3429443019238E+01
1975     (PID.TID 0000.0001) %MON ad_dynstat_advvel_del2 = 1.0630213533164E-01
1976     (PID.TID 0000.0001) %MON ad_dynstat_adwvel_max = 9.6291369464428E+03
1977     (PID.TID 0000.0001) %MON ad_dynstat_adwvel_min = -8.0212598087632E+04
1978     (PID.TID 0000.0001) %MON ad_dynstat_adwvel_mean = -2.8566486397398E+02
1979     (PID.TID 0000.0001) %MON ad_dynstat_adwvel_sd = 3.1351831455068E+03
1980     (PID.TID 0000.0001) %MON ad_dynstat_adwvel_del2 = 1.7209629792256E+01
1981     (PID.TID 0000.0001) %MON ad_dynstat_adtheta_max = 8.2651719664676E+01
1982     (PID.TID 0000.0001) %MON ad_dynstat_adtheta_min = -4.8880614970974E+01
1983     (PID.TID 0000.0001) %MON ad_dynstat_adtheta_mean = 4.6180891894873E-01
1984     (PID.TID 0000.0001) %MON ad_dynstat_adtheta_sd = 4.5299733854328E+00
1985     (PID.TID 0000.0001) %MON ad_dynstat_adtheta_del2 = 2.0428647971614E-02
1986     (PID.TID 0000.0001) %MON ad_dynstat_adsalt_max = 1.7196765112889E+02
1987     (PID.TID 0000.0001) %MON ad_dynstat_adsalt_min = -1.7060821975003E+02
1988     (PID.TID 0000.0001) %MON ad_dynstat_adsalt_mean = 3.8434763110441E-02
1989     (PID.TID 0000.0001) %MON ad_dynstat_adsalt_sd = 2.8167187528791E+00
1990     (PID.TID 0000.0001) %MON ad_dynstat_adsalt_del2 = 4.7705225892906E-02
1991     (PID.TID 0000.0001) %MON ad_dynstat_adsst_max = 8.2651719664676E+01
1992     (PID.TID 0000.0001) %MON ad_dynstat_adsst_min = -3.6562234334533E+00
1993     (PID.TID 0000.0001) %MON ad_dynstat_adsst_mean = 3.5225201180820E+01
1994     (PID.TID 0000.0001) %MON ad_dynstat_adsst_sd = 1.7747800371256E+01
1995     (PID.TID 0000.0001) %MON ad_dynstat_adsst_del2 = 1.9384329002905E-01
1996     (PID.TID 0000.0001) %MON ad_dynstat_adsss_max = 3.6292301866090E+01
1997     (PID.TID 0000.0001) %MON ad_dynstat_adsss_min = -1.7060821975003E+02
1998     (PID.TID 0000.0001) %MON ad_dynstat_adsss_mean = -3.4511622170179E-01
1999     (PID.TID 0000.0001) %MON ad_dynstat_adsss_sd = 6.4015883269168E+00
2000     (PID.TID 0000.0001) %MON ad_dynstat_adsss_del2 = 3.8654377162560E-01
2001 jmc 1.3 (PID.TID 0000.0001) %MON ad_forcing_adqnet_max = 0.0000000000000E+00
2002     (PID.TID 0000.0001) %MON ad_forcing_adqnet_min = 0.0000000000000E+00
2003     (PID.TID 0000.0001) %MON ad_forcing_adqnet_mean = 0.0000000000000E+00
2004     (PID.TID 0000.0001) %MON ad_forcing_adqnet_sd = 0.0000000000000E+00
2005     (PID.TID 0000.0001) %MON ad_forcing_adqnet_del2 = 0.0000000000000E+00
2006 jmc 1.5 (PID.TID 0000.0001) %MON ad_forcing_adempmr_max = 2.5195521881915E+03
2007     (PID.TID 0000.0001) %MON ad_forcing_adempmr_min = -1.1738957059563E+04
2008     (PID.TID 0000.0001) %MON ad_forcing_adempmr_mean = -1.4265369061821E+01
2009     (PID.TID 0000.0001) %MON ad_forcing_adempmr_sd = 4.3163518939887E+02
2010     (PID.TID 0000.0001) %MON ad_forcing_adempmr_del2 = 2.3357369077085E+01
2011 jmc 1.3 (PID.TID 0000.0001) %MON ad_forcing_adfu_max = 0.0000000000000E+00
2012     (PID.TID 0000.0001) %MON ad_forcing_adfu_min = 0.0000000000000E+00
2013     (PID.TID 0000.0001) %MON ad_forcing_adfu_mean = 0.0000000000000E+00
2014     (PID.TID 0000.0001) %MON ad_forcing_adfu_sd = 0.0000000000000E+00
2015     (PID.TID 0000.0001) %MON ad_forcing_adfu_del2 = 0.0000000000000E+00
2016     (PID.TID 0000.0001) %MON ad_forcing_adfv_max = 0.0000000000000E+00
2017     (PID.TID 0000.0001) %MON ad_forcing_adfv_min = 0.0000000000000E+00
2018     (PID.TID 0000.0001) %MON ad_forcing_adfv_mean = 0.0000000000000E+00
2019     (PID.TID 0000.0001) %MON ad_forcing_adfv_sd = 0.0000000000000E+00
2020     (PID.TID 0000.0001) %MON ad_forcing_adfv_del2 = 0.0000000000000E+00
2021 jmc 1.1 (PID.TID 0000.0001) // =======================================================
2022     (PID.TID 0000.0001) // End MONITOR dynamic field statistics
2023     (PID.TID 0000.0001) // =======================================================
2024     ph-pack: packing ecco_cost
2025     ph-pack: packing ecco_ctrl
2026 jmc 1.5 (PID.TID 0000.0001) // =======================================================
2027     (PID.TID 0000.0001) // Gradient-check starts (grdchk_main)
2028     (PID.TID 0000.0001) // =======================================================
2029     (PID.TID 0000.0001) grdchk reference fc: fcref = 8.04687744759401E+05
2030 jmc 1.1 grad-res -------------------------------
2031 jmc 1.5 grad-res proc # i j k bi bj iobc fc ref fc + eps fc - eps
2032     grad-res proc # i j k bi bj iobc adj grad fd grad 1 - fd/adj
2033 heimbach 1.2 grad-res closest next position:
2034     grad-res 0 2283 71 39 1 1 1
2035 jmc 1.5 (PID.TID 0000.0001) ====== Starts gradient-check number 1 (=ichknum) =======
2036 jmc 1.4 ph-test icomp, ncvarcomp, ichknum 2283 29309 1
2037     ph-grd _loc: bi, bj, icomptest, ichknum 1 1 0 1
2038 jmc 1.3 ph-grd -->hit<-- 73 39 1 1
2039 jmc 1.5 (PID.TID 0000.0001) grdchk pos: i,j,k= 73 39 1 ; bi,bj= 1 1 ; iobc= 1 ; rec= 1
2040 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2041     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2042     (PID.TID 0000.0001)
2043     (PID.TID 0000.0001) // =======================================================
2044     (PID.TID 0000.0001) // Model current state
2045     (PID.TID 0000.0001) // =======================================================
2046     (PID.TID 0000.0001)
2047 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2048     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2049 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.06234465668820E-14 2.53674886388737E+00
2050 jmc 1.5 cg2d: Sum(rhs),rhsMax = 1.25732757538799E-14 5.04325900480939E+00
2051     cg2d: Sum(rhs),rhsMax = -1.31734900765679E-14 6.29008403763971E+00
2052     cg2d: Sum(rhs),rhsMax = -7.24940940610708E-14 6.70003306073551E+00
2053 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2054 jmc 1.5 early fc = 0.000000000000000D+00
2055     --> objf_test(bi,bj) = 0.804687708291728D+06
2056 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2057     --> objf_atl(bi,bj) = 0.000000000000000D+00
2058 jmc 1.5 local fc = 0.804687708291728D+06
2059     global fc = 0.804687708291728D+06
2060     (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687708291728E+05
2061 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2062     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2063     (PID.TID 0000.0001)
2064     (PID.TID 0000.0001) // =======================================================
2065     (PID.TID 0000.0001) // Model current state
2066     (PID.TID 0000.0001) // =======================================================
2067     (PID.TID 0000.0001)
2068 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2069     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2070 jmc 1.1 cg2d: Sum(rhs),rhsMax = 1.05956909912663E-14 2.53674886388737E+00
2071 jmc 1.5 cg2d: Sum(rhs),rhsMax = -9.29811783123569E-16 5.04325900480939E+00
2072     cg2d: Sum(rhs),rhsMax = -1.94913529760754E-14 6.29008403763971E+00
2073     cg2d: Sum(rhs),rhsMax = -5.78946612872500E-14 6.70003306073551E+00
2074 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2075 jmc 1.5 early fc = 0.000000000000000D+00
2076     --> objf_test(bi,bj) = 0.804687744759401D+06
2077 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2078     --> objf_atl(bi,bj) = 0.000000000000000D+00
2079 jmc 1.5 local fc = 0.804687744759401D+06
2080     global fc = 0.804687744759401D+06
2081     (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687744759401E+05
2082 jmc 1.1 grad-res -------------------------------
2083 jmc 1.5 grad-res 0 1 73 39 1 1 1 1 8.04687744759E+05 8.04687708292E+05 8.04687744759E+05
2084     grad-res 0 1 1 2283 0 1 1 1 -3.65575776577E+00 -1.82338361046E+00 5.01229641763E-01
2085     (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759401E+05
2086     (PID.TID 0000.0001) ADM adjoint_gradient = -3.65575776577278E+00
2087     (PID.TID 0000.0001) ADM finite-diff_grad = -1.82338361046277E+00
2088     (PID.TID 0000.0001) ====== End of gradient-check number 1 (ierr= 0) =======
2089     (PID.TID 0000.0001) ====== Starts gradient-check number 2 (=ichknum) =======
2090 jmc 1.4 ph-test icomp, ncvarcomp, ichknum 2284 29309 2
2091     ph-grd _loc: bi, bj, icomptest, ichknum 1 1 2283 2
2092 jmc 1.3 ph-grd -->hit<-- 74 39 1 1
2093 jmc 1.5 (PID.TID 0000.0001) grdchk pos: i,j,k= 74 39 1 ; bi,bj= 1 1 ; iobc= 1 ; rec= 1
2094 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2095     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2096     (PID.TID 0000.0001)
2097     (PID.TID 0000.0001) // =======================================================
2098     (PID.TID 0000.0001) // Model current state
2099     (PID.TID 0000.0001) // =======================================================
2100     (PID.TID 0000.0001)
2101 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2102     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2103 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.06234465668820E-14 2.53674886388737E+00
2104 jmc 1.5 cg2d: Sum(rhs),rhsMax = 1.29826704942104E-14 5.04325900480939E+00
2105     cg2d: Sum(rhs),rhsMax = -1.90958360235527E-14 6.29008403763971E+00
2106     cg2d: Sum(rhs),rhsMax = -6.09512440519211E-14 6.70003306073551E+00
2107 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2108 jmc 1.5 early fc = 0.000000000000000D+00
2109     --> objf_test(bi,bj) = 0.804687708500864D+06
2110 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2111     --> objf_atl(bi,bj) = 0.000000000000000D+00
2112 jmc 1.5 local fc = 0.804687708500864D+06
2113     global fc = 0.804687708500864D+06
2114     (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687708500864E+05
2115 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2116     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2117     (PID.TID 0000.0001)
2118     (PID.TID 0000.0001) // =======================================================
2119     (PID.TID 0000.0001) // Model current state
2120     (PID.TID 0000.0001) // =======================================================
2121     (PID.TID 0000.0001)
2122 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2123     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2124 jmc 1.1 cg2d: Sum(rhs),rhsMax = 1.05956909912663E-14 2.53674886388737E+00
2125 jmc 1.5 cg2d: Sum(rhs),rhsMax = -9.29811783123569E-16 5.04325900480939E+00
2126     cg2d: Sum(rhs),rhsMax = -1.94913529760754E-14 6.29008403763971E+00
2127     cg2d: Sum(rhs),rhsMax = -5.78946612872500E-14 6.70003306073551E+00
2128 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2129 jmc 1.5 early fc = 0.000000000000000D+00
2130     --> objf_test(bi,bj) = 0.804687744759401D+06
2131 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2132     --> objf_atl(bi,bj) = 0.000000000000000D+00
2133 jmc 1.5 local fc = 0.804687744759401D+06
2134     global fc = 0.804687744759401D+06
2135     (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687744759401E+05
2136 jmc 1.1 grad-res -------------------------------
2137 jmc 1.5 grad-res 0 2 74 39 1 1 1 1 8.04687744759E+05 8.04687708501E+05 8.04687744759E+05
2138     grad-res 0 2 2 2284 0 1 1 1 -3.63465932636E+00 -1.81292684865E+00 5.01211341734E-01
2139     (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759401E+05
2140     (PID.TID 0000.0001) ADM adjoint_gradient = -3.63465932636325E+00
2141     (PID.TID 0000.0001) ADM finite-diff_grad = -1.81292684865184E+00
2142     (PID.TID 0000.0001) ====== End of gradient-check number 2 (ierr= 0) =======
2143     (PID.TID 0000.0001) ====== Starts gradient-check number 3 (=ichknum) =======
2144 jmc 1.4 ph-test icomp, ncvarcomp, ichknum 2285 29309 3
2145     ph-grd _loc: bi, bj, icomptest, ichknum 1 1 2284 3
2146 jmc 1.3 ph-grd -->hit<-- 75 39 1 1
2147 jmc 1.5 (PID.TID 0000.0001) grdchk pos: i,j,k= 75 39 1 ; bi,bj= 1 1 ; iobc= 1 ; rec= 1
2148 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2149     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2150     (PID.TID 0000.0001)
2151     (PID.TID 0000.0001) // =======================================================
2152     (PID.TID 0000.0001) // Model current state
2153     (PID.TID 0000.0001) // =======================================================
2154     (PID.TID 0000.0001)
2155 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2156     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2157 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.06234465668820E-14 2.53674886388737E+00
2158 jmc 1.5 cg2d: Sum(rhs),rhsMax = 1.34059430223488E-14 5.04325900480939E+00
2159     cg2d: Sum(rhs),rhsMax = -2.10213790818869E-14 6.29008403763971E+00
2160     cg2d: Sum(rhs),rhsMax = -5.99000016254791E-14 6.70003306073551E+00
2161 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2162 jmc 1.5 early fc = 0.000000000000000D+00
2163     --> objf_test(bi,bj) = 0.804687708506845D+06
2164 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2165     --> objf_atl(bi,bj) = 0.000000000000000D+00
2166 jmc 1.5 local fc = 0.804687708506845D+06
2167     global fc = 0.804687708506845D+06
2168     (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687708506845E+05
2169 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2170     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2171     (PID.TID 0000.0001)
2172     (PID.TID 0000.0001) // =======================================================
2173     (PID.TID 0000.0001) // Model current state
2174     (PID.TID 0000.0001) // =======================================================
2175     (PID.TID 0000.0001)
2176 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2177     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2178 jmc 1.1 cg2d: Sum(rhs),rhsMax = 1.05956909912663E-14 2.53674886388737E+00
2179 jmc 1.5 cg2d: Sum(rhs),rhsMax = -9.29811783123569E-16 5.04325900480939E+00
2180     cg2d: Sum(rhs),rhsMax = -1.94913529760754E-14 6.29008403763971E+00
2181     cg2d: Sum(rhs),rhsMax = -5.78946612872500E-14 6.70003306073551E+00
2182 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2183 jmc 1.5 early fc = 0.000000000000000D+00
2184     --> objf_test(bi,bj) = 0.804687744759401D+06
2185 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2186     --> objf_atl(bi,bj) = 0.000000000000000D+00
2187 jmc 1.5 local fc = 0.804687744759401D+06
2188     global fc = 0.804687744759401D+06
2189     (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687744759401E+05
2190 jmc 1.1 grad-res -------------------------------
2191 jmc 1.5 grad-res 0 3 75 39 1 1 1 1 8.04687744759E+05 8.04687708507E+05 8.04687744759E+05
2192     grad-res 0 3 3 2285 0 1 1 1 -3.63418584966E+00 -1.81262780097E+00 5.01228644885E-01
2193     (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759401E+05
2194     (PID.TID 0000.0001) ADM adjoint_gradient = -3.63418584965922E+00
2195     (PID.TID 0000.0001) ADM finite-diff_grad = -1.81262780097313E+00
2196     (PID.TID 0000.0001) ====== End of gradient-check number 3 (ierr= 0) =======
2197     (PID.TID 0000.0001) ====== Starts gradient-check number 4 (=ichknum) =======
2198 jmc 1.4 ph-test icomp, ncvarcomp, ichknum 2286 29309 4
2199     ph-grd _loc: bi, bj, icomptest, ichknum 1 1 2285 4
2200 jmc 1.3 ph-grd -->hit<-- 76 39 1 1
2201 jmc 1.5 (PID.TID 0000.0001) grdchk pos: i,j,k= 76 39 1 ; bi,bj= 1 1 ; iobc= 1 ; rec= 1
2202 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2203     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2204     (PID.TID 0000.0001)
2205     (PID.TID 0000.0001) // =======================================================
2206     (PID.TID 0000.0001) // Model current state
2207     (PID.TID 0000.0001) // =======================================================
2208     (PID.TID 0000.0001)
2209 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2210     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2211 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.06234465668820E-14 2.53674886388737E+00
2212 jmc 1.5 cg2d: Sum(rhs),rhsMax = -8.32667268468867E-17 5.04325900480939E+00
2213     cg2d: Sum(rhs),rhsMax = -2.22183382803109E-14 6.29008403763971E+00
2214     cg2d: Sum(rhs),rhsMax = -5.94559124156291E-14 6.70003306073551E+00
2215 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2216 jmc 1.5 early fc = 0.000000000000000D+00
2217     --> objf_test(bi,bj) = 0.804687712788332D+06
2218 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2219     --> objf_atl(bi,bj) = 0.000000000000000D+00
2220 jmc 1.5 local fc = 0.804687712788332D+06
2221     global fc = 0.804687712788332D+06
2222     (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687712788332E+05
2223 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2224     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2225     (PID.TID 0000.0001)
2226     (PID.TID 0000.0001) // =======================================================
2227     (PID.TID 0000.0001) // Model current state
2228     (PID.TID 0000.0001) // =======================================================
2229     (PID.TID 0000.0001)
2230 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2231     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2232 jmc 1.1 cg2d: Sum(rhs),rhsMax = 1.05956909912663E-14 2.53674886388737E+00
2233 jmc 1.5 cg2d: Sum(rhs),rhsMax = -9.29811783123569E-16 5.04325900480939E+00
2234     cg2d: Sum(rhs),rhsMax = -1.94913529760754E-14 6.29008403763971E+00
2235     cg2d: Sum(rhs),rhsMax = -5.78946612872500E-14 6.70003306073551E+00
2236 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2237 jmc 1.5 early fc = 0.000000000000000D+00
2238     --> objf_test(bi,bj) = 0.804687744759401D+06
2239 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2240     --> objf_atl(bi,bj) = 0.000000000000000D+00
2241 jmc 1.5 local fc = 0.804687744759401D+06
2242     global fc = 0.804687744759401D+06
2243     (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687744759401E+05
2244 jmc 1.1 grad-res -------------------------------
2245 jmc 1.5 grad-res 0 4 76 39 1 1 1 1 8.04687744759E+05 8.04687712788E+05 8.04687744759E+05
2246     grad-res 0 4 4 2286 0 1 1 1 -1.65446443816E-02 -1.59855341772E+00 -9.56205970249E+01
2247     (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759401E+05
2248     (PID.TID 0000.0001) ADM adjoint_gradient = -1.65446443816105E-02
2249     (PID.TID 0000.0001) ADM finite-diff_grad = -1.59855341771618E+00
2250     (PID.TID 0000.0001) ====== End of gradient-check number 4 (ierr= 0) =======
2251     (PID.TID 0000.0001) ====== Starts gradient-check number 5 (=ichknum) =======
2252 jmc 1.4 ph-test icomp, ncvarcomp, ichknum 2287 29309 5
2253     ph-grd _loc: bi, bj, icomptest, ichknum 1 1 2286 5
2254 jmc 1.3 ph-grd -->hit<-- 85 39 1 1
2255 jmc 1.5 (PID.TID 0000.0001) grdchk pos: i,j,k= 85 39 1 ; bi,bj= 1 1 ; iobc= 1 ; rec= 1
2256 heimbach 1.2 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2257     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2258     (PID.TID 0000.0001)
2259     (PID.TID 0000.0001) // =======================================================
2260     (PID.TID 0000.0001) // Model current state
2261     (PID.TID 0000.0001) // =======================================================
2262     (PID.TID 0000.0001)
2263 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2264     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2265 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.05956909912663E-14 2.53674886388737E+00
2266 jmc 1.5 cg2d: Sum(rhs),rhsMax = 1.49602552568240E-14 5.04325900480939E+00
2267     cg2d: Sum(rhs),rhsMax = -1.53939361258182E-14 6.29008403763971E+00
2268     cg2d: Sum(rhs),rhsMax = -3.48401862915182E-14 6.70003306073551E+00
2269 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2270 jmc 1.5 early fc = 0.000000000000000D+00
2271     --> objf_test(bi,bj) = 0.804687749551412D+06
2272 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2273     --> objf_atl(bi,bj) = 0.000000000000000D+00
2274 jmc 1.5 local fc = 0.804687749551412D+06
2275     global fc = 0.804687749551412D+06
2276     (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687749551412E+05
2277 heimbach 1.2 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2278     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2279     (PID.TID 0000.0001)
2280     (PID.TID 0000.0001) // =======================================================
2281     (PID.TID 0000.0001) // Model current state
2282     (PID.TID 0000.0001) // =======================================================
2283     (PID.TID 0000.0001)
2284 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2285     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2286 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.05956909912663E-14 2.53674886388737E+00
2287 jmc 1.5 cg2d: Sum(rhs),rhsMax = 1.18377530000657E-14 5.04325900480939E+00
2288     cg2d: Sum(rhs),rhsMax = -1.20944920745103E-14 6.29008403763971E+00
2289     cg2d: Sum(rhs),rhsMax = -4.83744988510892E-14 6.70003306073551E+00
2290 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2291 jmc 1.5 early fc = 0.000000000000000D+00
2292     --> objf_test(bi,bj) = 0.804687740147925D+06
2293 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2294     --> objf_atl(bi,bj) = 0.000000000000000D+00
2295 jmc 1.5 local fc = 0.804687740147925D+06
2296     global fc = 0.804687740147925D+06
2297     (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687740147925E+05
2298 heimbach 1.2 grad-res -------------------------------
2299 jmc 1.5 grad-res 0 5 85 39 1 1 1 1 8.04687744759E+05 8.04687749551E+05 8.04687740148E+05
2300     grad-res 0 5 5 2287 0 1 1 1 4.70174487688E-01 4.70174360089E-01 2.71386260242E-07
2301     (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759401E+05
2302     (PID.TID 0000.0001) ADM adjoint_gradient = 4.70174487687900E-01
2303     (PID.TID 0000.0001) ADM finite-diff_grad = 4.70174360089004E-01
2304     (PID.TID 0000.0001) ====== End of gradient-check number 5 (ierr= 0) =======
2305     (PID.TID 0000.0001) ====== Starts gradient-check number 6 (=ichknum) =======
2306 jmc 1.4 ph-test icomp, ncvarcomp, ichknum 2288 29309 6
2307     ph-grd _loc: bi, bj, icomptest, ichknum 1 1 2287 6
2308 jmc 1.3 ph-grd -->hit<-- 86 39 1 1
2309 jmc 1.5 (PID.TID 0000.0001) grdchk pos: i,j,k= 86 39 1 ; bi,bj= 1 1 ; iobc= 1 ; rec= 1
2310 heimbach 1.2 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2311     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2312     (PID.TID 0000.0001)
2313     (PID.TID 0000.0001) // =======================================================
2314     (PID.TID 0000.0001) // Model current state
2315     (PID.TID 0000.0001) // =======================================================
2316     (PID.TID 0000.0001)
2317 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2318     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2319 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.06234465668820E-14 2.53674886388737E+00
2320 jmc 1.5 cg2d: Sum(rhs),rhsMax = 7.74380559676047E-15 5.04325900480939E+00
2321     cg2d: Sum(rhs),rhsMax = -1.41796296926344E-14 6.29008403763971E+00
2322     cg2d: Sum(rhs),rhsMax = -2.16354711923827E-14 6.70003306073551E+00
2323 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2324 jmc 1.5 early fc = 0.000000000000000D+00
2325     --> objf_test(bi,bj) = 0.804687751308914D+06
2326 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2327     --> objf_atl(bi,bj) = 0.000000000000000D+00
2328 jmc 1.5 local fc = 0.804687751308914D+06
2329     global fc = 0.804687751308914D+06
2330     (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687751308914E+05
2331 heimbach 1.2 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2332     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2333     (PID.TID 0000.0001)
2334     (PID.TID 0000.0001) // =======================================================
2335     (PID.TID 0000.0001) // Model current state
2336     (PID.TID 0000.0001) // =======================================================
2337     (PID.TID 0000.0001)
2338 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2339     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2340 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.05956909912663E-14 2.53674886388737E+00
2341 jmc 1.5 cg2d: Sum(rhs),rhsMax = 1.48769885299771E-14 5.04325900480939E+00
2342     cg2d: Sum(rhs),rhsMax = -6.53643805748061E-15 6.29008403763971E+00
2343     cg2d: Sum(rhs),rhsMax = -1.75866265994529E-14 6.70003306073551E+00
2344 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2345 jmc 1.5 early fc = 0.000000000000000D+00
2346     --> objf_test(bi,bj) = 0.804687738384548D+06
2347 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2348     --> objf_atl(bi,bj) = 0.000000000000000D+00
2349 jmc 1.5 local fc = 0.804687738384548D+06
2350     global fc = 0.804687738384548D+06
2351     (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687738384548E+05
2352 heimbach 1.2 grad-res -------------------------------
2353 jmc 1.5 grad-res 0 6 86 39 1 1 1 1 8.04687744759E+05 8.04687751309E+05 8.04687738385E+05
2354     grad-res 0 6 6 2288 0 1 1 1 6.46218052670E-01 6.46218325710E-01 -4.22520225252E-07
2355     (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759401E+05
2356     (PID.TID 0000.0001) ADM adjoint_gradient = 6.46218052669727E-01
2357     (PID.TID 0000.0001) ADM finite-diff_grad = 6.46218325709924E-01
2358     (PID.TID 0000.0001) ====== End of gradient-check number 6 (ierr= 0) =======
2359     (PID.TID 0000.0001) ====== Starts gradient-check number 7 (=ichknum) =======
2360 jmc 1.4 ph-test icomp, ncvarcomp, ichknum 2289 29309 7
2361     ph-grd _loc: bi, bj, icomptest, ichknum 1 1 2288 7
2362 jmc 1.3 ph-grd -->hit<-- 87 39 1 1
2363 jmc 1.5 (PID.TID 0000.0001) grdchk pos: i,j,k= 87 39 1 ; bi,bj= 1 1 ; iobc= 1 ; rec= 1
2364 heimbach 1.2 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2365     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2366     (PID.TID 0000.0001)
2367     (PID.TID 0000.0001) // =======================================================
2368     (PID.TID 0000.0001) // Model current state
2369     (PID.TID 0000.0001) // =======================================================
2370     (PID.TID 0000.0001)
2371 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2372     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2373 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.06512021424976E-14 2.53674886388737E+00
2374 jmc 1.5 cg2d: Sum(rhs),rhsMax = 9.28424004342787E-15 5.04325900480939E+00
2375     cg2d: Sum(rhs),rhsMax = -2.49869569479699E-14 6.29008403763971E+00
2376     cg2d: Sum(rhs),rhsMax = -2.69367861349679E-14 6.70003306073551E+00
2377 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2378 jmc 1.5 early fc = 0.000000000000000D+00
2379     --> objf_test(bi,bj) = 0.804687749828302D+06
2380 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2381     --> objf_atl(bi,bj) = 0.000000000000000D+00
2382 jmc 1.5 local fc = 0.804687749828302D+06
2383     global fc = 0.804687749828302D+06
2384     (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687749828302E+05
2385 heimbach 1.2 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2386     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2387     (PID.TID 0000.0001)
2388     (PID.TID 0000.0001) // =======================================================
2389     (PID.TID 0000.0001) // Model current state
2390     (PID.TID 0000.0001) // =======================================================
2391     (PID.TID 0000.0001)
2392 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2393     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2394 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.06512021424976E-14 2.53674886388737E+00
2395 jmc 1.5 cg2d: Sum(rhs),rhsMax = 6.73766598069392E-15 5.04325900480939E+00
2396     cg2d: Sum(rhs),rhsMax = -1.15185638804860E-14 6.29008403763971E+00
2397     cg2d: Sum(rhs),rhsMax = -6.89726054048379E-14 6.70003306073551E+00
2398 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2399 jmc 1.5 early fc = 0.000000000000000D+00
2400     --> objf_test(bi,bj) = 0.804687739866778D+06
2401 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2402     --> objf_atl(bi,bj) = 0.000000000000000D+00
2403 jmc 1.5 local fc = 0.804687739866778D+06
2404     global fc = 0.804687739866778D+06
2405     (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687739866778E+05
2406 heimbach 1.2 grad-res -------------------------------
2407 jmc 1.5 grad-res 0 7 87 39 1 1 1 1 8.04687744759E+05 8.04687749828E+05 8.04687739867E+05
2408     grad-res 0 7 7 2289 0 1 1 1 4.98076169749E-01 4.98076213989E-01 -8.88226561191E-08
2409     (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759401E+05
2410     (PID.TID 0000.0001) ADM adjoint_gradient = 4.98076169748959E-01
2411     (PID.TID 0000.0001) ADM finite-diff_grad = 4.98076213989407E-01
2412     (PID.TID 0000.0001) ====== End of gradient-check number 7 (ierr= 0) =======
2413     (PID.TID 0000.0001) ====== Starts gradient-check number 8 (=ichknum) =======
2414 jmc 1.4 ph-test icomp, ncvarcomp, ichknum 2290 29309 8
2415     ph-grd _loc: bi, bj, icomptest, ichknum 1 1 2289 8
2416 jmc 1.3 ph-grd -->hit<-- 88 39 1 1
2417 jmc 1.5 (PID.TID 0000.0001) grdchk pos: i,j,k= 88 39 1 ; bi,bj= 1 1 ; iobc= 1 ; rec= 1
2418 heimbach 1.2 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2419     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2420     (PID.TID 0000.0001)
2421     (PID.TID 0000.0001) // =======================================================
2422     (PID.TID 0000.0001) // Model current state
2423     (PID.TID 0000.0001) // =======================================================
2424     (PID.TID 0000.0001)
2425 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2426     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2427 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.06373243546898E-14 2.53674886388737E+00
2428 jmc 1.5 cg2d: Sum(rhs),rhsMax = 1.77635683940025E-14 5.04325900480939E+00
2429     cg2d: Sum(rhs),rhsMax = -3.35842464949110E-14 6.29008403763971E+00
2430     cg2d: Sum(rhs),rhsMax = -3.34871019802563E-14 6.70003306073551E+00
2431 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2432 jmc 1.5 early fc = 0.000000000000000D+00
2433     --> objf_test(bi,bj) = 0.804687750914794D+06
2434 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2435     --> objf_atl(bi,bj) = 0.000000000000000D+00
2436 jmc 1.5 local fc = 0.804687750914794D+06
2437     global fc = 0.804687750914794D+06
2438     (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687750914794E+05
2439 heimbach 1.2 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2440     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2441     (PID.TID 0000.0001)
2442     (PID.TID 0000.0001) // =======================================================
2443     (PID.TID 0000.0001) // Model current state
2444     (PID.TID 0000.0001) // =======================================================
2445     (PID.TID 0000.0001)
2446 jmc 1.3 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 )
2447     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2448 heimbach 1.2 cg2d: Sum(rhs),rhsMax = 1.06095687790742E-14 2.53674886388737E+00
2449 jmc 1.5 cg2d: Sum(rhs),rhsMax = 2.23779328401008E-14 5.04325900480939E+00
2450     cg2d: Sum(rhs),rhsMax = -1.94289029309402E-14 6.29008403763971E+00
2451     cg2d: Sum(rhs),rhsMax = -7.71743779992562E-14 6.70003306073551E+00
2452 jmc 1.3 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2453 jmc 1.5 early fc = 0.000000000000000D+00
2454     --> objf_test(bi,bj) = 0.804687738782400D+06
2455 jmc 1.3 --> objf_tracer(bi,bj) = 0.000000000000000D+00
2456     --> objf_atl(bi,bj) = 0.000000000000000D+00
2457 jmc 1.5 local fc = 0.804687738782400D+06
2458     global fc = 0.804687738782400D+06
2459     (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687738782400E+05
2460 heimbach 1.2 grad-res -------------------------------
2461 jmc 1.5 grad-res 0 8 88 39 1 1 1 1 8.04687744759E+05 8.04687750915E+05 8.04687738782E+05
2462     grad-res 0 8 8 2290 0 1 1 1 6.06619515888E-01 6.06619729660E-01 -3.52399837977E-07
2463     (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759401E+05
2464     (PID.TID 0000.0001) ADM adjoint_gradient = 6.06619515887832E-01
2465     (PID.TID 0000.0001) ADM finite-diff_grad = 6.06619729660451E-01
2466     (PID.TID 0000.0001) ====== End of gradient-check number 8 (ierr= 0) =======
2467 heimbach 1.2 (PID.TID 0000.0001)
2468     (PID.TID 0000.0001) // =======================================================
2469     (PID.TID 0000.0001) // Gradient check results >>> START <<<
2470     (PID.TID 0000.0001) // =======================================================
2471     (PID.TID 0000.0001)
2472 jmc 1.3 (PID.TID 0000.0001) EPS = 1.000000E-02
2473 heimbach 1.2 (PID.TID 0000.0001)
2474 jmc 1.4 (PID.TID 0000.0001) grdchk output h.p: Id Itile Jtile LAYER bi bj X(Id) X(Id)+/-EPS
2475 jmc 1.3 (PID.TID 0000.0001) grdchk output h.c: Id FC FC1 FC2
2476     (PID.TID 0000.0001) grdchk output h.g: Id FC1-FC2/(2*EPS) ADJ GRAD(FC) 1-FDGRD/ADGRD
2477     (PID.TID 0000.0001)
2478 jmc 1.4 (PID.TID 0000.0001) grdchk output (p): 1 73 39 1 1 1 0.000000000E+00 -1.000000000E-02
2479 jmc 1.5 (PID.TID 0000.0001) grdchk output (c): 1 8.0468774475940E+05 8.0468770829173E+05 8.0468774475940E+05
2480     (PID.TID 0000.0001) grdchk output (g): 1 -1.8233836104628E+00 -3.6557577657728E+00 5.0122964176284E-01
2481 jmc 1.3 (PID.TID 0000.0001)
2482 jmc 1.4 (PID.TID 0000.0001) grdchk output (p): 2 74 39 1 1 1 0.000000000E+00 -1.000000000E-02
2483 jmc 1.5 (PID.TID 0000.0001) grdchk output (c): 2 8.0468774475940E+05 8.0468770850086E+05 8.0468774475940E+05
2484     (PID.TID 0000.0001) grdchk output (g): 2 -1.8129268486518E+00 -3.6346593263633E+00 5.0121134173369E-01
2485 jmc 1.3 (PID.TID 0000.0001)
2486 jmc 1.4 (PID.TID 0000.0001) grdchk output (p): 3 75 39 1 1 1 0.000000000E+00 -1.000000000E-02
2487 jmc 1.5 (PID.TID 0000.0001) grdchk output (c): 3 8.0468774475940E+05 8.0468770850684E+05 8.0468774475940E+05
2488     (PID.TID 0000.0001) grdchk output (g): 3 -1.8126278009731E+00 -3.6341858496592E+00 5.0122864488532E-01
2489 jmc 1.3 (PID.TID 0000.0001)
2490 jmc 1.4 (PID.TID 0000.0001) grdchk output (p): 4 76 39 1 1 1 0.000000000E+00 -1.000000000E-02
2491 jmc 1.5 (PID.TID 0000.0001) grdchk output (c): 4 8.0468774475940E+05 8.0468771278833E+05 8.0468774475940E+05
2492     (PID.TID 0000.0001) grdchk output (g): 4 -1.5985534177162E+00 -1.6544644381610E-02 -9.5620597024919E+01
2493 jmc 1.3 (PID.TID 0000.0001)
2494 jmc 1.4 (PID.TID 0000.0001) grdchk output (p): 5 85 39 1 1 1 0.000000000E+00 -1.000000000E-02
2495 jmc 1.5 (PID.TID 0000.0001) grdchk output (c): 5 8.0468774475940E+05 8.0468774955141E+05 8.0468774014792E+05
2496     (PID.TID 0000.0001) grdchk output (g): 5 4.7017436008900E-01 4.7017448768790E-01 2.7138626024215E-07
2497 jmc 1.3 (PID.TID 0000.0001)
2498 jmc 1.4 (PID.TID 0000.0001) grdchk output (p): 6 86 39 1 1 1 0.000000000E+00 -1.000000000E-02
2499 jmc 1.5 (PID.TID 0000.0001) grdchk output (c): 6 8.0468774475940E+05 8.0468775130891E+05 8.0468773838455E+05
2500     (PID.TID 0000.0001) grdchk output (g): 6 6.4621832570992E-01 6.4621805266973E-01 -4.2252022525169E-07
2501 jmc 1.3 (PID.TID 0000.0001)
2502 jmc 1.4 (PID.TID 0000.0001) grdchk output (p): 7 87 39 1 1 1 0.000000000E+00 -1.000000000E-02
2503 jmc 1.5 (PID.TID 0000.0001) grdchk output (c): 7 8.0468774475940E+05 8.0468774982830E+05 8.0468773986678E+05
2504     (PID.TID 0000.0001) grdchk output (g): 7 4.9807621398941E-01 4.9807616974896E-01 -8.8822656119092E-08
2505 jmc 1.3 (PID.TID 0000.0001)
2506 jmc 1.4 (PID.TID 0000.0001) grdchk output (p): 8 88 39 1 1 1 0.000000000E+00 -1.000000000E-02
2507 jmc 1.5 (PID.TID 0000.0001) grdchk output (c): 8 8.0468774475940E+05 8.0468775091479E+05 8.0468773878240E+05
2508     (PID.TID 0000.0001) grdchk output (g): 8 6.0661972966045E-01 6.0661951588783E-01 -3.5239983797730E-07
2509     (PID.TID 0000.0001)
2510     (PID.TID 0000.0001) grdchk summary : RMS of 8 ratios = 3.3808379600181E+01
2511 heimbach 1.2 (PID.TID 0000.0001)
2512     (PID.TID 0000.0001) // =======================================================
2513     (PID.TID 0000.0001) // Gradient check results >>> END <<<
2514 jmc 1.1 (PID.TID 0000.0001) // =======================================================
2515     (PID.TID 0000.0001)
2516     (PID.TID 0000.0001) Seconds in section "ALL [THE_MODEL_MAIN]":
2517 jmc 1.5 (PID.TID 0000.0001) User time: 67.030000000000001
2518     (PID.TID 0000.0001) System time: 0.20999999999999999
2519     (PID.TID 0000.0001) Wall clock time: 67.419028997421265
2520 jmc 1.3 (PID.TID 0000.0001) No. starts: 1
2521     (PID.TID 0000.0001) No. stops: 1
2522 jmc 1.1 (PID.TID 0000.0001) Seconds in section "INITIALISE_FIXED [THE_MODEL_MAIN]":
2523 jmc 1.5 (PID.TID 0000.0001) User time: 0.17000000000000001
2524 jmc 1.3 (PID.TID 0000.0001) System time: 2.00000000000000004E-002
2525 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 0.19813203811645508
2526 jmc 1.3 (PID.TID 0000.0001) No. starts: 1
2527     (PID.TID 0000.0001) No. stops: 1
2528 jmc 1.1 (PID.TID 0000.0001) Seconds in section "ADTHE_MAIN_LOOP [ADJOINT RUN]":
2529 jmc 1.5 (PID.TID 0000.0001) User time: 17.829999999999998
2530     (PID.TID 0000.0001) System time: 0.15000000000000002
2531     (PID.TID 0000.0001) Wall clock time: 18.029143095016479
2532 jmc 1.3 (PID.TID 0000.0001) No. starts: 1
2533     (PID.TID 0000.0001) No. stops: 1
2534 jmc 1.5 (PID.TID 0000.0001) Seconds in section "FORWARD_STEP [MAIN_DO_LOOP]":
2535     (PID.TID 0000.0001) User time: 51.079999999999927
2536     (PID.TID 0000.0001) System time: 2.00000000000000178E-002
2537     (PID.TID 0000.0001) Wall clock time: 51.219977378845215
2538     (PID.TID 0000.0001) No. starts: 72
2539     (PID.TID 0000.0001) No. stops: 72
2540 jmc 1.1 (PID.TID 0000.0001) Seconds in section "LOAD_FIELDS_DRIVER [FORWARD_STEP]":
2541 jmc 1.5 (PID.TID 0000.0001) User time: 9.99999999999943157E-002
2542 jmc 1.4 (PID.TID 0000.0001) System time: 0.0000000000000000
2543 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 0.10257363319396973
2544 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2545     (PID.TID 0000.0001) No. stops: 72
2546 jmc 1.1 (PID.TID 0000.0001) Seconds in section "EXTERNAL_FLDS_LOAD [LOAD_FLDS_DRIVER]":
2547 jmc 1.5 (PID.TID 0000.0001) User time: 9.99999999999943157E-002
2548 jmc 1.4 (PID.TID 0000.0001) System time: 0.0000000000000000
2549 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 0.10687303543090820
2550 jmc 1.3 (PID.TID 0000.0001) No. starts: 76
2551     (PID.TID 0000.0001) No. stops: 76
2552 jmc 1.1 (PID.TID 0000.0001) Seconds in section "DO_ATMOSPHERIC_PHYS [FORWARD_STEP]":
2553 jmc 1.5 (PID.TID 0000.0001) User time: 9.99999999999090505E-003
2554 jmc 1.3 (PID.TID 0000.0001) System time: 0.0000000000000000
2555 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 6.55889511108398438E-004
2556 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2557     (PID.TID 0000.0001) No. stops: 72
2558 jmc 1.1 (PID.TID 0000.0001) Seconds in section "DO_OCEANIC_PHYS [FORWARD_STEP]":
2559 jmc 1.5 (PID.TID 0000.0001) User time: 8.1899999999999977
2560     (PID.TID 0000.0001) System time: 1.00000000000000089E-002
2561     (PID.TID 0000.0001) Wall clock time: 8.2164540290832520
2562 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2563     (PID.TID 0000.0001) No. stops: 72
2564 jmc 1.1 (PID.TID 0000.0001) Seconds in section "THERMODYNAMICS [FORWARD_STEP]":
2565 jmc 1.5 (PID.TID 0000.0001) User time: 16.469999999999970
2566 jmc 1.4 (PID.TID 0000.0001) System time: 0.0000000000000000
2567 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 16.527510404586792
2568 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2569     (PID.TID 0000.0001) No. stops: 72
2570 jmc 1.1 (PID.TID 0000.0001) Seconds in section "DYNAMICS [FORWARD_STEP]":
2571 jmc 1.5 (PID.TID 0000.0001) User time: 19.130000000000010
2572 jmc 1.4 (PID.TID 0000.0001) System time: 1.00000000000000089E-002
2573 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 19.160159349441528
2574 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2575     (PID.TID 0000.0001) No. stops: 72
2576 jmc 1.1 (PID.TID 0000.0001) Seconds in section "SOLVE_FOR_PRESSURE [FORWARD_STEP]":
2577 jmc 1.5 (PID.TID 0000.0001) User time: 4.3899999999999864
2578     (PID.TID 0000.0001) System time: 0.0000000000000000
2579     (PID.TID 0000.0001) Wall clock time: 4.3996253013610840
2580 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2581     (PID.TID 0000.0001) No. stops: 72
2582 jmc 1.1 (PID.TID 0000.0001) Seconds in section "MOM_CORRECTION_STEP [FORWARD_STEP]":
2583 jmc 1.5 (PID.TID 0000.0001) User time: 0.39999999999999147
2584 jmc 1.3 (PID.TID 0000.0001) System time: 0.0000000000000000
2585 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 0.42065262794494629
2586     (PID.TID 0000.0001) No. starts: 72
2587     (PID.TID 0000.0001) No. stops: 72
2588     (PID.TID 0000.0001) Seconds in section "INTEGR_CONTINUITY [FORWARD_STEP]":
2589     (PID.TID 0000.0001) User time: 0.72000000000004150
2590     (PID.TID 0000.0001) System time: 0.0000000000000000
2591     (PID.TID 0000.0001) Wall clock time: 0.71262645721435547
2592 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2593     (PID.TID 0000.0001) No. stops: 72
2594 jmc 1.1 (PID.TID 0000.0001) Seconds in section "TRC_CORRECTION_STEP [FORWARD_STEP]":
2595 jmc 1.5 (PID.TID 0000.0001) User time: 0.21000000000000796
2596 jmc 1.3 (PID.TID 0000.0001) System time: 0.0000000000000000
2597 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 0.21013736724853516
2598 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2599     (PID.TID 0000.0001) No. stops: 72
2600 jmc 1.1 (PID.TID 0000.0001) Seconds in section "BLOCKING_EXCHANGES [FORWARD_STEP]":
2601 jmc 1.5 (PID.TID 0000.0001) User time: 0.28999999999997783
2602 jmc 1.3 (PID.TID 0000.0001) System time: 0.0000000000000000
2603 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 0.27254319190979004
2604 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2605     (PID.TID 0000.0001) No. stops: 72
2606 jmc 1.1 (PID.TID 0000.0001) Seconds in section "MONITOR [FORWARD_STEP]":
2607 jmc 1.5 (PID.TID 0000.0001) User time: 1.00000000000051159E-002
2608 jmc 1.3 (PID.TID 0000.0001) System time: 0.0000000000000000
2609 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 7.18832015991210938E-004
2610 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2611     (PID.TID 0000.0001) No. stops: 72
2612 jmc 1.1 (PID.TID 0000.0001) Seconds in section "COST_TILE [FORWARD_STEP]":
2613 jmc 1.5 (PID.TID 0000.0001) User time: 1.0799999999999557
2614 jmc 1.3 (PID.TID 0000.0001) System time: 0.0000000000000000
2615 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 1.1138741970062256
2616 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2617     (PID.TID 0000.0001) No. stops: 72
2618 jmc 1.1 (PID.TID 0000.0001) Seconds in section "DO_THE_MODEL_IO [FORWARD_STEP]":
2619 jmc 1.5 (PID.TID 0000.0001) User time: 1.99999999999960210E-002
2620     (PID.TID 0000.0001) System time: 0.0000000000000000
2621     (PID.TID 0000.0001) Wall clock time: 2.31816768646240234E-002
2622 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2623     (PID.TID 0000.0001) No. stops: 72
2624 jmc 1.1 (PID.TID 0000.0001) Seconds in section "DO_WRITE_PICKUP [FORWARD_STEP]":
2625 jmc 1.5 (PID.TID 0000.0001) User time: 5.00000000000113687E-002
2626     (PID.TID 0000.0001) System time: 0.0000000000000000
2627     (PID.TID 0000.0001) Wall clock time: 4.88457679748535156E-002
2628 jmc 1.3 (PID.TID 0000.0001) No. starts: 72
2629     (PID.TID 0000.0001) No. stops: 72
2630 jmc 1.5 (PID.TID 0000.0001) Seconds in section "I/O (WRITE) [ADJOINT LOOP]":
2631     (PID.TID 0000.0001) User time: 4.00000000000027001E-002
2632     (PID.TID 0000.0001) System time: 0.0000000000000000
2633     (PID.TID 0000.0001) Wall clock time: 4.38201427459716797E-002
2634     (PID.TID 0000.0001) No. starts: 1
2635     (PID.TID 0000.0001) No. stops: 1
2636 jmc 1.1 (PID.TID 0000.0001) Seconds in section "CTRL_PACK [THE_MODEL_MAIN]":
2637 jmc 1.5 (PID.TID 0000.0001) User time: 5.00000000000007105E-002
2638     (PID.TID 0000.0001) System time: 9.99999999999998113E-003
2639     (PID.TID 0000.0001) Wall clock time: 5.79869747161865234E-002
2640 jmc 1.3 (PID.TID 0000.0001) No. starts: 1
2641     (PID.TID 0000.0001) No. stops: 1
2642 jmc 1.1 (PID.TID 0000.0001) Seconds in section "CTRL_PACK [THE_MODEL_MAIN]":
2643 jmc 1.5 (PID.TID 0000.0001) User time: 5.99999999999987210E-002
2644 jmc 1.3 (PID.TID 0000.0001) System time: 0.0000000000000000
2645 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 5.73489665985107422E-002
2646 jmc 1.3 (PID.TID 0000.0001) No. starts: 1
2647     (PID.TID 0000.0001) No. stops: 1
2648 jmc 1.1 (PID.TID 0000.0001) Seconds in section "GRDCHK_MAIN [THE_MODEL_MAIN]":
2649 jmc 1.5 (PID.TID 0000.0001) User time: 48.920000000000002
2650     (PID.TID 0000.0001) System time: 2.99999999999999989E-002
2651     (PID.TID 0000.0001) Wall clock time: 49.076337099075317
2652 jmc 1.3 (PID.TID 0000.0001) No. starts: 1
2653     (PID.TID 0000.0001) No. stops: 1
2654 jmc 1.1 (PID.TID 0000.0001) Seconds in section "INITIALISE_VARIA [THE_MAIN_LOOP]":
2655 jmc 1.5 (PID.TID 0000.0001) User time: 1.3499999999999943
2656     (PID.TID 0000.0001) System time: 9.99999999999998113E-003
2657     (PID.TID 0000.0001) Wall clock time: 1.3673486709594727
2658 jmc 1.3 (PID.TID 0000.0001) No. starts: 16
2659     (PID.TID 0000.0001) No. stops: 16
2660 jmc 1.1 (PID.TID 0000.0001) Seconds in section "MAIN LOOP [THE_MAIN_LOOP]":
2661 jmc 1.5 (PID.TID 0000.0001) User time: 47.329999999999998
2662     (PID.TID 0000.0001) System time: 0.0000000000000000
2663     (PID.TID 0000.0001) Wall clock time: 47.438304901123047
2664 jmc 1.3 (PID.TID 0000.0001) No. starts: 16
2665     (PID.TID 0000.0001) No. stops: 16
2666 jmc 1.5 (PID.TID 0000.0001) Seconds in section "MAIN_DO_LOOP [THE_MAIN_LOOP]":
2667     (PID.TID 0000.0001) User time: 45.339999999999961
2668     (PID.TID 0000.0001) System time: 0.0000000000000000
2669     (PID.TID 0000.0001) Wall clock time: 45.451010704040527
2670 jmc 1.3 (PID.TID 0000.0001) No. starts: 64
2671     (PID.TID 0000.0001) No. stops: 64
2672 jmc 1.1 (PID.TID 0000.0001) Seconds in section "COST_FINAL [ADJOINT SPIN-DOWN]":
2673 jmc 1.3 (PID.TID 0000.0001) User time: 0.0000000000000000
2674     (PID.TID 0000.0001) System time: 0.0000000000000000
2675 jmc 1.5 (PID.TID 0000.0001) Wall clock time: 1.77431106567382813E-003
2676 jmc 1.3 (PID.TID 0000.0001) No. starts: 16
2677     (PID.TID 0000.0001) No. stops: 16
2678 jmc 1.1 (PID.TID 0000.0001) // ======================================================
2679     (PID.TID 0000.0001) // Tile <-> Tile communication statistics
2680     (PID.TID 0000.0001) // ======================================================
2681     (PID.TID 0000.0001) // o Tile number: 000001
2682     (PID.TID 0000.0001) // No. X exchanges = 0
2683     (PID.TID 0000.0001) // Max. X spins = 0
2684     (PID.TID 0000.0001) // Min. X spins = 1000000000
2685     (PID.TID 0000.0001) // Total. X spins = 0
2686     (PID.TID 0000.0001) // Avg. X spins = 0.00E+00
2687     (PID.TID 0000.0001) // No. Y exchanges = 0
2688     (PID.TID 0000.0001) // Max. Y spins = 0
2689     (PID.TID 0000.0001) // Min. Y spins = 1000000000
2690     (PID.TID 0000.0001) // Total. Y spins = 0
2691     (PID.TID 0000.0001) // Avg. Y spins = 0.00E+00
2692     (PID.TID 0000.0001) // o Thread number: 000001
2693 jmc 1.5 (PID.TID 0000.0001) // No. barriers = 19380
2694 jmc 1.1 (PID.TID 0000.0001) // Max. barrier spins = 1
2695     (PID.TID 0000.0001) // Min. barrier spins = 1
2696 jmc 1.5 (PID.TID 0000.0001) // Total barrier spins = 19380
2697 jmc 1.1 (PID.TID 0000.0001) // Avg. barrier spins = 1.00E+00
2698 jmc 1.3 PROGRAM MAIN: Execution ended Normally
  ViewVC Help
Powered by ViewVC 1.1.22