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

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

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

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