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

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

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


Revision 1.1 - (hide annotations) (download)
Sun Jul 5 01:54:11 2015 UTC (8 years, 10 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65n, checkpoint65o, HEAD
File MIME type: text/plain
Add kpp to verification.

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

  ViewVC Help
Powered by ViewVC 1.1.22