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