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

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

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


Revision 1.9 - (hide annotations) (download)
Wed Jun 26 22:39:18 2013 UTC (10 years, 11 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint65, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, HEAD
Changes since 1.8: +439 -1261 lines
File MIME type: text/plain
reduce length of the test to 20.iters (same as fwd test)

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.9 (PID.TID 0000.0001) // MITgcmUV version: checkpoint64i
9 jmc 1.4 (PID.TID 0000.0001) // Build user: jmc
10 jmc 1.6 (PID.TID 0000.0001) // Build host: baudelaire
11 jmc 1.9 (PID.TID 0000.0001) // Build date: Thu Jun 20 17:46:13 EDT 2013
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 jmc 1.6 (PID.TID 0000.0001) > /
24 jmc 1.1 (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.5 (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.5 (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.7 (PID.TID 0000.0001) debugMode = F ; /* print debug msg. (sequence of S/R calls) */
55 jmc 1.5 (PID.TID 0000.0001) printMapIncludesZeros= F ; /* print zeros in Std.Output maps */
56     (PID.TID 0000.0001) maxLengthPrt1D= 65 /* maxLength of 1D array printed to StdOut */
57 jmc 1.1 (PID.TID 0000.0001)
58     (PID.TID 0000.0001) // ======================================================
59     (PID.TID 0000.0001) // Mapping of tiles to threads
60     (PID.TID 0000.0001) // ======================================================
61     (PID.TID 0000.0001) // -o- Thread 1, tiles ( 1: 1, 1: 1)
62     (PID.TID 0000.0001)
63     (PID.TID 0000.0001) // ======================================================
64     (PID.TID 0000.0001) // Tile <-> Tile connectvity table
65     (PID.TID 0000.0001) // ======================================================
66 jmc 1.9 (PID.TID 0000.0001) // Tile number: 000001 (process no. = 000000)
67     (PID.TID 0000.0001) // WEST: Tile = 000001, Process = 000000, Comm = put
68 jmc 1.1 (PID.TID 0000.0001) // bi = 000001, bj = 000001
69 jmc 1.9 (PID.TID 0000.0001) // EAST: Tile = 000001, Process = 000000, Comm = put
70 jmc 1.1 (PID.TID 0000.0001) // bi = 000001, bj = 000001
71 jmc 1.9 (PID.TID 0000.0001) // SOUTH: Tile = 000001, Process = 000000, Comm = put
72 jmc 1.1 (PID.TID 0000.0001) // bi = 000001, bj = 000001
73 jmc 1.9 (PID.TID 0000.0001) // NORTH: Tile = 000001, Process = 000000, Comm = put
74 jmc 1.1 (PID.TID 0000.0001) // bi = 000001, bj = 000001
75     (PID.TID 0000.0001)
76 jmc 1.7 (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.7 (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 jmc 1.7 (PID.TID 0000.0001) > diffKhT=0.,
92 jmc 1.1 (PID.TID 0000.0001) > diffKrT=3.E-5,
93 jmc 1.7 (PID.TID 0000.0001) > diffKhS=0.,
94 jmc 1.1 (PID.TID 0000.0001) > diffKrS=3.E-5,
95 jmc 1.7 (PID.TID 0000.0001) > rhoConst=1035.,
96     (PID.TID 0000.0001) > rhoConstFresh=1000.,
97 jmc 1.1 (PID.TID 0000.0001) > eosType = 'POLY3',
98 jmc 1.7 (PID.TID 0000.0001) > ivdc_kappa=10.,
99 jmc 1.1 (PID.TID 0000.0001) > implicitDiffusion=.TRUE.,
100 jmc 1.7 (PID.TID 0000.0001) > allowFreezing=.TRUE.,
101     (PID.TID 0000.0001) >#exactConserv=.TRUE.,
102     (PID.TID 0000.0001) >#- runoff missing in SSH forcing => comment out realFreshWaterFlux
103     (PID.TID 0000.0001) >#useRealFreshWaterFlux=.TRUE.,
104 jmc 1.1 (PID.TID 0000.0001) > useCDscheme=.TRUE.,
105     (PID.TID 0000.0001) ># turn on looped cells
106     (PID.TID 0000.0001) > hFacMin=.05,
107     (PID.TID 0000.0001) > hFacMindr=50.,
108     (PID.TID 0000.0001) ># set precision of data files
109     (PID.TID 0000.0001) > readBinaryPrec=32,
110 jmc 1.6 (PID.TID 0000.0001) > /
111 jmc 1.1 (PID.TID 0000.0001) >
112     (PID.TID 0000.0001) ># Elliptic solver parameters
113     (PID.TID 0000.0001) > &PARM02
114     (PID.TID 0000.0001) > cg2dMaxIters=500,
115     (PID.TID 0000.0001) > cg2dTargetResidual=1.E-13,
116 jmc 1.6 (PID.TID 0000.0001) > /
117 jmc 1.1 (PID.TID 0000.0001) >
118     (PID.TID 0000.0001) ># Time stepping parameters
119     (PID.TID 0000.0001) > &PARM03
120     (PID.TID 0000.0001) > nIter0 = 0,
121 jmc 1.9 (PID.TID 0000.0001) > nTimeSteps = 20,
122 jmc 1.1 (PID.TID 0000.0001) ># 100 years of integration will yield a reasonable flow field
123     (PID.TID 0000.0001) ># startTime = 0.,
124     (PID.TID 0000.0001) ># endTime = 3110400000.,
125 jmc 1.7 (PID.TID 0000.0001) > deltaTmom = 1200.,
126 jmc 1.1 (PID.TID 0000.0001) > tauCD = 321428.,
127 jmc 1.7 (PID.TID 0000.0001) > deltaTtracer = 43200.,
128     (PID.TID 0000.0001) > deltaTClock = 43200.,
129     (PID.TID 0000.0001) > deltaTfreesurf=43200.,
130 jmc 1.1 (PID.TID 0000.0001) > abEps = 0.1,
131     (PID.TID 0000.0001) > pChkptFreq= 311040000.,
132     (PID.TID 0000.0001) > dumpFreq= 311040000.,
133 jmc 1.7 (PID.TID 0000.0001) > taveFreq= 311040000.,
134 jmc 1.1 (PID.TID 0000.0001) > dumpFreq= 864000.,
135     (PID.TID 0000.0001) > taveFreq= 864000.,
136     (PID.TID 0000.0001) >#monitorFreq=31104000.,
137 jmc 1.7 (PID.TID 0000.0001) > monitorFreq=2592000.,
138 jmc 1.1 (PID.TID 0000.0001) ># 2 months restoring timescale for temperature
139 jmc 1.7 (PID.TID 0000.0001) >#tauThetaClimRelax = 5184000.0,
140 jmc 1.1 (PID.TID 0000.0001) ># 6 months restoring timescale for salinity
141     (PID.TID 0000.0001) > tauSaltClimRelax = 15552000.0,
142 jmc 1.7 (PID.TID 0000.0001) >#periodicExternalForcing=.TRUE.,
143     (PID.TID 0000.0001) >#externForcingPeriod=2592000.,
144     (PID.TID 0000.0001) >#externForcingCycle=31104000.,
145 jmc 1.6 (PID.TID 0000.0001) > /
146 jmc 1.1 (PID.TID 0000.0001) >
147     (PID.TID 0000.0001) ># Gridding parameters
148     (PID.TID 0000.0001) > &PARM04
149     (PID.TID 0000.0001) > usingSphericalPolarGrid=.TRUE.,
150 jmc 1.7 (PID.TID 0000.0001) > delR= 50., 70., 100., 140., 190.,
151     (PID.TID 0000.0001) > 240., 290., 340., 390., 440.,
152     (PID.TID 0000.0001) > 490., 540., 590., 640., 690.,
153 jmc 1.5 (PID.TID 0000.0001) > ygOrigin=-80.,
154 jmc 1.1 (PID.TID 0000.0001) > dySpacing=4.,
155     (PID.TID 0000.0001) > dxSpacing=4.,
156 jmc 1.6 (PID.TID 0000.0001) > /
157 jmc 1.1 (PID.TID 0000.0001) >
158     (PID.TID 0000.0001) ># Input datasets
159     (PID.TID 0000.0001) > &PARM05
160     (PID.TID 0000.0001) > bathyFile= 'Topog.bin',
161     (PID.TID 0000.0001) > hydrogThetaFile='Lev_clim_theta.bin',
162     (PID.TID 0000.0001) > hydrogSaltFile= 'Lev_clim_salt.bin',
163     (PID.TID 0000.0001) > zonalWindFile= 'Taux.bin',
164     (PID.TID 0000.0001) > meridWindFile= 'Tauy.bin',
165 jmc 1.7 (PID.TID 0000.0001) >#thetaClimFile= 'SST.bin',
166 jmc 1.1 (PID.TID 0000.0001) > saltClimFile= 'SSS.bin',
167 jmc 1.6 (PID.TID 0000.0001) > /
168 jmc 1.1 (PID.TID 0000.0001)
169 jmc 1.7 (PID.TID 0000.0001) INI_PARMS ; starts to read PARM01
170     (PID.TID 0000.0001) INI_PARMS ; read PARM01 : OK
171 jmc 1.1 (PID.TID 0000.0001) S/R INI_PARMS: No request for barotropic solver
172     (PID.TID 0000.0001) S/R INI_PARMS: => Use implicitFreeSurface as default
173 jmc 1.7 (PID.TID 0000.0001) INI_PARMS ; starts to read PARM02
174     (PID.TID 0000.0001) INI_PARMS ; read PARM02 : OK
175     (PID.TID 0000.0001) INI_PARMS ; starts to read PARM03
176     (PID.TID 0000.0001) INI_PARMS ; read PARM03 : OK
177     (PID.TID 0000.0001) INI_PARMS ; starts to read PARM04
178     (PID.TID 0000.0001) INI_PARMS ; read PARM04 : OK
179     (PID.TID 0000.0001) INI_PARMS ; starts to read PARM05
180     (PID.TID 0000.0001) INI_PARMS ; read PARM05 : OK
181     (PID.TID 0000.0001) INI_PARMS: finished reading file "data"
182 jmc 1.1 (PID.TID 0000.0001) PACKAGES_BOOT: opening data.pkg
183     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.pkg
184     (PID.TID 0000.0001) // =======================================================
185     (PID.TID 0000.0001) // Parameter file "data.pkg"
186     (PID.TID 0000.0001) // =======================================================
187     (PID.TID 0000.0001) ># Packages
188     (PID.TID 0000.0001) > &PACKAGES
189     (PID.TID 0000.0001) > useGMRedi=.TRUE.,
190     (PID.TID 0000.0001) > useEBM=.TRUE.,
191     (PID.TID 0000.0001) > useGrdchk=.TRUE.,
192 jmc 1.6 (PID.TID 0000.0001) > /
193 jmc 1.1 (PID.TID 0000.0001) >
194     (PID.TID 0000.0001)
195     (PID.TID 0000.0001) PACKAGES_BOOT: finished reading data.pkg
196     (PID.TID 0000.0001) GM_READPARMS: opening data.gmredi
197     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.gmredi
198     (PID.TID 0000.0001) // =======================================================
199     (PID.TID 0000.0001) // Parameter file "data.gmredi"
200     (PID.TID 0000.0001) // =======================================================
201     (PID.TID 0000.0001) ># GM+Redi package parameters:
202     (PID.TID 0000.0001) >
203     (PID.TID 0000.0001) >#-from MOM :
204     (PID.TID 0000.0001) ># GM_background_K: G & Mc.W diffusion coefficient
205     (PID.TID 0000.0001) ># GM_maxSlope : max slope of isopycnals
206     (PID.TID 0000.0001) ># GM_Scrit : transition for scaling diffusion coefficient
207     (PID.TID 0000.0001) ># GM_Sd : half width scaling for diffusion coefficient
208     (PID.TID 0000.0001) ># GM_taper_scheme: slope clipping or one of the tapering schemes
209     (PID.TID 0000.0001) ># GM_Kmin_horiz : horizontal diffusion minimum value
210     (PID.TID 0000.0001) >
211     (PID.TID 0000.0001) >#-Option parameters (needs to "define" options in GMREDI_OPTIONS.h")
212     (PID.TID 0000.0001) ># GM_isopycK : isopycnal diffusion coefficient (default=GM_background_K)
213     (PID.TID 0000.0001) ># GM_AdvForm : turn on GM Advective form (default=Skew flux form)
214     (PID.TID 0000.0001) >
215     (PID.TID 0000.0001) > &GM_PARM01
216     (PID.TID 0000.0001) > GM_background_K = 1.e+3,
217     (PID.TID 0000.0001) > GM_taper_scheme = 'gkw91',
218     (PID.TID 0000.0001) > GM_maxSlope = 1.e-2,
219     (PID.TID 0000.0001) > GM_Kmin_horiz = 50.,
220     (PID.TID 0000.0001) > GM_Scrit = 4.e-3,
221     (PID.TID 0000.0001) > GM_Sd = 1.e-3,
222     (PID.TID 0000.0001) ># GM_Visbeck_alpha = 0.,
223     (PID.TID 0000.0001) ># GM_Visbeck_length = 2.e+5,
224     (PID.TID 0000.0001) ># GM_Visbeck_depth = 1.e+3,
225     (PID.TID 0000.0001) ># GM_Visbeck_maxval_K= 2.5e+3,
226 jmc 1.6 (PID.TID 0000.0001) > /
227 jmc 1.1 (PID.TID 0000.0001) >
228     (PID.TID 0000.0001) >
229     (PID.TID 0000.0001)
230     (PID.TID 0000.0001) GM_READPARMS: finished reading data.gmredi
231 jmc 1.9 (PID.TID 0000.0001) EBM_READPARMS: opening data.ebm
232     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.ebm
233     (PID.TID 0000.0001) // =======================================================
234     (PID.TID 0000.0001) // Parameter file "data.ebm"
235     (PID.TID 0000.0001) // =======================================================
236     (PID.TID 0000.0001) ># Packages
237     (PID.TID 0000.0001) > &EBM_PARM01
238     (PID.TID 0000.0001) ># convert from mm/year to m/sec
239     (PID.TID 0000.0001) > scale_runoff = 3.170979198E-11,
240     (PID.TID 0000.0001) > tauThetaZonRelax = 5184000.0,
241     (PID.TID 0000.0001) > RunoffFile = 'Runoff.bin',
242     (PID.TID 0000.0001) > /
243     (PID.TID 0000.0001)
244     (PID.TID 0000.0001) EBM_READPARMS: finished reading data.ebm
245     (PID.TID 0000.0001) EBM_READPARMS: set tauThetaClimRelax to tauThetaZonRelax value
246     (PID.TID 0000.0001) AUTODIFF_READPARMS: opening data.autodiff
247     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.autodiff
248     (PID.TID 0000.0001) // =======================================================
249     (PID.TID 0000.0001) // Parameter file "data.autodiff"
250     (PID.TID 0000.0001) // =======================================================
251     (PID.TID 0000.0001) ># =========================
252     (PID.TID 0000.0001) ># pkg AUTODIFF parameters :
253     (PID.TID 0000.0001) ># =========================
254     (PID.TID 0000.0001) ># inAdExact :: get an exact adjoint (no approximation) (def=.True.)
255     (PID.TID 0000.0001) >#
256     (PID.TID 0000.0001) > &AUTODIFF_PARM01
257     (PID.TID 0000.0001) ># inAdExact = .FALSE.,
258     (PID.TID 0000.0001) > /
259     (PID.TID 0000.0001)
260     (PID.TID 0000.0001) AUTODIFF_READPARMS: finished reading data.autodiff
261     (PID.TID 0000.0001) // ===================================
262     (PID.TID 0000.0001) // AUTODIFF parameters :
263     (PID.TID 0000.0001) // ===================================
264     (PID.TID 0000.0001) inAdExact = /* get an exact adjoint (no approximation) */
265     (PID.TID 0000.0001) T
266     (PID.TID 0000.0001) ;
267     (PID.TID 0000.0001) useKPPinAdMode = /* use KPP in adjoint mode */
268     (PID.TID 0000.0001) F
269     (PID.TID 0000.0001) ;
270     (PID.TID 0000.0001) useGMRediInAdMode = /* use GMRedi in adjoint mode */
271     (PID.TID 0000.0001) T
272     (PID.TID 0000.0001) ;
273     (PID.TID 0000.0001) useSEAICEinAdMode = /* use SEAICE in adjoint mode */
274     (PID.TID 0000.0001) F
275     (PID.TID 0000.0001) ;
276     (PID.TID 0000.0001) useGGL90inAdMode = /* use GGL90 in adjoint mode */
277     (PID.TID 0000.0001) F
278     (PID.TID 0000.0001) ;
279     (PID.TID 0000.0001) useSALT_PLUMEinAdMode = /* use SALT_PLUME in adjoint mode */
280     (PID.TID 0000.0001) F
281     (PID.TID 0000.0001) ;
282     (PID.TID 0000.0001) dumpAdVarExch = /* control adexch before dumpinp */
283     (PID.TID 0000.0001) 2
284     (PID.TID 0000.0001) ;
285     (PID.TID 0000.0001) mon_AdVarExch = /* control adexch before monitor */
286     (PID.TID 0000.0001) 2
287     (PID.TID 0000.0001) ;
288     (PID.TID 0000.0001)
289 jmc 1.1 (PID.TID 0000.0001) OPTIM_READPARMS: opening data.optim
290     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.optim
291     (PID.TID 0000.0001) // =======================================================
292     (PID.TID 0000.0001) // Parameter file "data.optim"
293     (PID.TID 0000.0001) // =======================================================
294     (PID.TID 0000.0001) >#
295     (PID.TID 0000.0001) ># ********************************
296     (PID.TID 0000.0001) ># Off-line optimization parameters
297     (PID.TID 0000.0001) ># ********************************
298     (PID.TID 0000.0001) > &OPTIM
299     (PID.TID 0000.0001) > optimcycle=0,
300 jmc 1.6 (PID.TID 0000.0001) > /
301 jmc 1.1 (PID.TID 0000.0001)
302     (PID.TID 0000.0001) OPTIM_READPARMS: finished reading data.optim
303     (PID.TID 0000.0001) CTRL_READPARMS: opening data.ctrl
304     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.ctrl
305     (PID.TID 0000.0001) // =======================================================
306     (PID.TID 0000.0001) // Parameter file "data.ctrl"
307     (PID.TID 0000.0001) // =======================================================
308     (PID.TID 0000.0001) >#
309     (PID.TID 0000.0001) >#
310     (PID.TID 0000.0001) ># *********************
311     (PID.TID 0000.0001) ># ECCO controlvariables
312     (PID.TID 0000.0001) ># *********************
313     (PID.TID 0000.0001) > &CTRL_NML
314     (PID.TID 0000.0001) > xx_theta_file = 'xx_theta',
315     (PID.TID 0000.0001) > xx_salt_file = 'xx_salt',
316     (PID.TID 0000.0001) > xx_tr1_file = 'xx_tr1',
317     (PID.TID 0000.0001) > xx_hflux_file = 'xx_hflux',
318     (PID.TID 0000.0001) > xx_sflux_file = 'xx_sflux',
319     (PID.TID 0000.0001) > xx_tauu_file = 'xx_tauu',
320     (PID.TID 0000.0001) > xx_tauv_file = 'xx_tauv',
321     (PID.TID 0000.0001) > xx_diffkr_file = 'xx_diffkr',
322     (PID.TID 0000.0001) > xx_kapgm_file = 'xx_kapgm',
323 jmc 1.6 (PID.TID 0000.0001) > /
324 jmc 1.1 (PID.TID 0000.0001) >#
325     (PID.TID 0000.0001) ># *********************
326     (PID.TID 0000.0001) ># names for ctrl_pack/unpack
327     (PID.TID 0000.0001) ># *********************
328     (PID.TID 0000.0001) > &CTRL_PACKNAMES
329 jmc 1.6 (PID.TID 0000.0001) > /
330 jmc 1.1 (PID.TID 0000.0001) >
331     (PID.TID 0000.0001)
332     (PID.TID 0000.0001) CTRL_READPARMS: finished reading data.ctrl
333     (PID.TID 0000.0001) COST_READPARMS: opening data.cost
334     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.cost
335     (PID.TID 0000.0001) // =======================================================
336     (PID.TID 0000.0001) // Parameter file "data.cost"
337     (PID.TID 0000.0001) // =======================================================
338     (PID.TID 0000.0001) >#
339     (PID.TID 0000.0001) >#
340     (PID.TID 0000.0001) ># ******************
341     (PID.TID 0000.0001) ># ECCO cost function
342     (PID.TID 0000.0001) ># ******************
343     (PID.TID 0000.0001) > &COST_NML
344     (PID.TID 0000.0001) >#
345     (PID.TID 0000.0001) > mult_tracer = 1.,
346     (PID.TID 0000.0001) > mult_test = 1.,
347     (PID.TID 0000.0001) > mult_atl = 1.,
348 jmc 1.6 (PID.TID 0000.0001) > /
349 jmc 1.1 (PID.TID 0000.0001)
350     (PID.TID 0000.0001) COST_READPARMS: finished reading data.cost
351     (PID.TID 0000.0001) GRDCHK_READPARMS: opening data.grdchk
352     (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.grdchk
353     (PID.TID 0000.0001) // =======================================================
354     (PID.TID 0000.0001) // Parameter file "data.grdchk"
355     (PID.TID 0000.0001) // =======================================================
356     (PID.TID 0000.0001) >
357     (PID.TID 0000.0001) ># *******************
358     (PID.TID 0000.0001) ># ECCO gradient check
359     (PID.TID 0000.0001) ># *******************
360     (PID.TID 0000.0001) > &GRDCHK_NML
361     (PID.TID 0000.0001) > grdchk_eps = 1.d-2,
362     (PID.TID 0000.0001) > nbeg = 1,
363     (PID.TID 0000.0001) > nstep = 1,
364     (PID.TID 0000.0001) > nend = 2,
365     (PID.TID 0000.0001) > grdchkvarindex = 1,
366 jmc 1.6 (PID.TID 0000.0001) > /
367 jmc 1.1 (PID.TID 0000.0001)
368     (PID.TID 0000.0001) GRDCHK_READPARMS: finished reading data.grdchk
369     (PID.TID 0000.0001)
370     (PID.TID 0000.0001) // =======================================================
371     (PID.TID 0000.0001) // Gradient check configuration >>> START <<<
372     (PID.TID 0000.0001) // =======================================================
373     (PID.TID 0000.0001)
374     (PID.TID 0000.0001) eps: 0.100E-01
375     (PID.TID 0000.0001) First location: 1
376     (PID.TID 0000.0001) Last location: 2
377     (PID.TID 0000.0001) Increment: 1
378 jmc 1.7 (PID.TID 0000.0001) grdchkWhichProc: 0
379     (PID.TID 0000.0001) iLocTile = 1 , jLocTile = 1
380 jmc 1.1 (PID.TID 0000.0001)
381     (PID.TID 0000.0001) // =======================================================
382     (PID.TID 0000.0001) // Gradient check configuration >>> END <<<
383     (PID.TID 0000.0001) // =======================================================
384     (PID.TID 0000.0001)
385     (PID.TID 0000.0001) SET_PARMS: done
386     (PID.TID 0000.0001) Enter INI_VERTICAL_GRID: setInterFDr= T ; setCenterDr= F
387 jmc 1.7 (PID.TID 0000.0001) %MON XC_max = 3.5800000000000E+02
388     (PID.TID 0000.0001) %MON XC_min = 2.0000000000000E+00
389     (PID.TID 0000.0001) %MON XC_mean = 1.8000000000000E+02
390     (PID.TID 0000.0001) %MON XC_sd = 1.0391663325314E+02
391     (PID.TID 0000.0001) %MON XG_max = 3.5600000000000E+02
392     (PID.TID 0000.0001) %MON XG_min = 0.0000000000000E+00
393     (PID.TID 0000.0001) %MON XG_mean = 1.7800000000000E+02
394     (PID.TID 0000.0001) %MON XG_sd = 1.0391663325314E+02
395     (PID.TID 0000.0001) %MON DXC_max = 4.4443898815675E+05
396     (PID.TID 0000.0001) %MON DXC_min = 9.2460385861875E+04
397     (PID.TID 0000.0001) %MON DXC_mean = 3.1372497618153E+05
398     (PID.TID 0000.0001) %MON DXC_sd = 1.1216447457560E+05
399     (PID.TID 0000.0001) %MON DXF_max = 4.4443898815675E+05
400     (PID.TID 0000.0001) %MON DXF_min = 9.2460385861875E+04
401     (PID.TID 0000.0001) %MON DXF_mean = 3.1372497618153E+05
402     (PID.TID 0000.0001) %MON DXF_sd = 1.1216447457560E+05
403     (PID.TID 0000.0001) %MON DXG_max = 4.4470989340816E+05
404     (PID.TID 0000.0001) %MON DXG_min = 7.7223062580781E+04
405     (PID.TID 0000.0001) %MON DXG_mean = 3.1353386340260E+05
406     (PID.TID 0000.0001) %MON DXG_sd = 1.1256651772502E+05
407     (PID.TID 0000.0001) %MON DXV_max = 4.4470989340816E+05
408     (PID.TID 0000.0001) %MON DXV_min = 7.7223062580781E+04
409     (PID.TID 0000.0001) %MON DXV_mean = 3.1353386340260E+05
410     (PID.TID 0000.0001) %MON DXV_sd = 1.1256651772502E+05
411     (PID.TID 0000.0001) %MON YC_max = 7.8000000000000E+01
412     (PID.TID 0000.0001) %MON YC_min = -7.8000000000000E+01
413     (PID.TID 0000.0001) %MON YC_mean = 0.0000000000000E+00
414     (PID.TID 0000.0001) %MON YC_sd = 4.6173585522461E+01
415     (PID.TID 0000.0001) %MON YG_max = 7.6000000000000E+01
416     (PID.TID 0000.0001) %MON YG_min = -8.0000000000000E+01
417     (PID.TID 0000.0001) %MON YG_mean = -2.0000000000000E+00
418     (PID.TID 0000.0001) %MON YG_sd = 4.6173585522461E+01
419     (PID.TID 0000.0001) %MON DYC_max = 4.4470989340816E+05
420     (PID.TID 0000.0001) %MON DYC_min = 4.4470989340816E+05
421     (PID.TID 0000.0001) %MON DYC_mean = 4.4470989340814E+05
422     (PID.TID 0000.0001) %MON DYC_sd = 1.9848812371492E-08
423     (PID.TID 0000.0001) %MON DYF_max = 4.4470989340816E+05
424     (PID.TID 0000.0001) %MON DYF_min = 4.4470989340816E+05
425     (PID.TID 0000.0001) %MON DYF_mean = 4.4470989340814E+05
426     (PID.TID 0000.0001) %MON DYF_sd = 1.9848812371492E-08
427     (PID.TID 0000.0001) %MON DYG_max = 4.4470989340816E+05
428     (PID.TID 0000.0001) %MON DYG_min = 4.4470989340816E+05
429     (PID.TID 0000.0001) %MON DYG_mean = 4.4470989340814E+05
430     (PID.TID 0000.0001) %MON DYG_sd = 1.9848812371492E-08
431     (PID.TID 0000.0001) %MON DYU_max = 4.4470989340816E+05
432     (PID.TID 0000.0001) %MON DYU_min = 4.4470989340816E+05
433     (PID.TID 0000.0001) %MON DYU_mean = 4.4470989340814E+05
434     (PID.TID 0000.0001) %MON DYU_sd = 1.9848812371492E-08
435     (PID.TID 0000.0001) %MON RA_max = 1.9760627980089E+11
436     (PID.TID 0000.0001) %MON RA_min = 4.1109698667290E+10
437     (PID.TID 0000.0001) %MON RA_mean = 1.3948826965197E+11
438     (PID.TID 0000.0001) %MON RA_sd = 4.9870522472902E+10
439     (PID.TID 0000.0001) %MON RAW_max = 1.9760627980089E+11
440     (PID.TID 0000.0001) %MON RAW_min = 4.1109698667290E+10
441     (PID.TID 0000.0001) %MON RAW_mean = 1.3948826965197E+11
442     (PID.TID 0000.0001) %MON RAW_sd = 4.9870522472902E+10
443     (PID.TID 0000.0001) %MON RAS_max = 1.9772672958215E+11
444     (PID.TID 0000.0001) %MON RAS_min = 3.4334886267983E+10
445     (PID.TID 0000.0001) %MON RAS_mean = 1.3940329716694E+11
446     (PID.TID 0000.0001) %MON RAS_sd = 5.0049278732354E+10
447     (PID.TID 0000.0001) %MON RAZ_max = 1.9772672958215E+11
448     (PID.TID 0000.0001) %MON RAZ_min = 3.4334886267983E+10
449     (PID.TID 0000.0001) %MON RAZ_mean = 1.3940329716694E+11
450     (PID.TID 0000.0001) %MON RAZ_sd = 5.0049278732354E+10
451     (PID.TID 0000.0001) %MON AngleCS_max = 1.0000000000000E+00
452     (PID.TID 0000.0001) %MON AngleCS_min = 1.0000000000000E+00
453     (PID.TID 0000.0001) %MON AngleCS_mean = 1.0000000000000E+00
454     (PID.TID 0000.0001) %MON AngleCS_sd = 0.0000000000000E+00
455     (PID.TID 0000.0001) %MON AngleSN_max = 0.0000000000000E+00
456     (PID.TID 0000.0001) %MON AngleSN_min = 0.0000000000000E+00
457     (PID.TID 0000.0001) %MON AngleSN_mean = 0.0000000000000E+00
458     (PID.TID 0000.0001) %MON AngleSN_sd = 0.0000000000000E+00
459     (PID.TID 0000.0001) GAD_INIT_FIXED: GAD_OlMinSize= 1 0 1
460 jmc 1.1 (PID.TID 0000.0001)
461     (PID.TID 0000.0001) // ===================================
462     (PID.TID 0000.0001) // GAD parameters :
463     (PID.TID 0000.0001) // ===================================
464     (PID.TID 0000.0001) tempAdvScheme = /* Temp. Horiz.Advection scheme selector */
465     (PID.TID 0000.0001) 2
466     (PID.TID 0000.0001) ;
467     (PID.TID 0000.0001) tempVertAdvScheme = /* Temp. Vert. Advection scheme selector */
468     (PID.TID 0000.0001) 2
469     (PID.TID 0000.0001) ;
470     (PID.TID 0000.0001) tempMultiDimAdvec = /* use Muti-Dim Advec method for Temp */
471     (PID.TID 0000.0001) F
472     (PID.TID 0000.0001) ;
473     (PID.TID 0000.0001) tempSOM_Advection = /* use 2nd Order Moment Advection for Temp */
474     (PID.TID 0000.0001) F
475     (PID.TID 0000.0001) ;
476     (PID.TID 0000.0001) AdamsBashforthGt = /* apply Adams-Bashforth extrapolation on Gt */
477     (PID.TID 0000.0001) T
478     (PID.TID 0000.0001) ;
479     (PID.TID 0000.0001) AdamsBashforth_T = /* apply Adams-Bashforth extrapolation on Temp */
480     (PID.TID 0000.0001) F
481     (PID.TID 0000.0001) ;
482     (PID.TID 0000.0001) saltAdvScheme = /* Salt. Horiz.advection scheme selector */
483     (PID.TID 0000.0001) 2
484     (PID.TID 0000.0001) ;
485     (PID.TID 0000.0001) saltVertAdvScheme = /* Salt. Vert. Advection scheme selector */
486     (PID.TID 0000.0001) 2
487     (PID.TID 0000.0001) ;
488     (PID.TID 0000.0001) saltMultiDimAdvec = /* use Muti-Dim Advec method for Salt */
489     (PID.TID 0000.0001) F
490     (PID.TID 0000.0001) ;
491     (PID.TID 0000.0001) saltSOM_Advection = /* use 2nd Order Moment Advection for Salt */
492     (PID.TID 0000.0001) F
493     (PID.TID 0000.0001) ;
494     (PID.TID 0000.0001) AdamsBashforthGs = /* apply Adams-Bashforth extrapolation on Gs */
495     (PID.TID 0000.0001) T
496     (PID.TID 0000.0001) ;
497     (PID.TID 0000.0001) AdamsBashforth_S = /* apply Adams-Bashforth extrapolation on Salt */
498     (PID.TID 0000.0001) F
499     (PID.TID 0000.0001) ;
500     (PID.TID 0000.0001) // ===================================
501     (PID.TID 0000.0001) ctrl-wet 1: nvarlength = 91470
502     (PID.TID 0000.0001) ctrl-wet 2: surface wet C = 2158
503     (PID.TID 0000.0001) ctrl-wet 3: surface wet W = 2052
504     (PID.TID 0000.0001) ctrl-wet 4: surface wet S = 2004
505     (PID.TID 0000.0001) ctrl-wet 4a:surface wet V = 0
506     (PID.TID 0000.0001) ctrl-wet 5: 3D wet points = 29138
507     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 1 1
508     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 2 1
509     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 3 0
510     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 4 0
511     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 5 1
512     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 6 1
513     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 7 0
514     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 8 0
515     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 9 0
516     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 10 0
517     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 11 0
518     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 12 0
519     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 13 0
520     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 14 0
521     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 15 1
522     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 16 0
523     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 17 0
524     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 18 0
525     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 19 0
526     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 20 0
527     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 21 0
528     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 22 0
529     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 23 0
530     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 24 0
531     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 25 0
532     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 26 0
533     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 27 0
534     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 28 0
535     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 29 0
536     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 30 0
537     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 31 0
538     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 32 0
539     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 33 0
540     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 34 0
541     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 35 0
542     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 36 0
543     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 37 0
544     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 38 0
545     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 39 0
546     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 40 0
547     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 41 0
548     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 42 0
549     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 43 0
550     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 44 0
551     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 45 0
552     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 46 0
553     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 47 0
554     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 48 0
555     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 49 0
556     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 50 0
557     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 51 0
558     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 52 0
559     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 53 0
560     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 54 0
561     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 55 0
562     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 56 0
563     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 57 0
564     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 58 0
565     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 59 0
566     (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 60 0
567     (PID.TID 0000.0001) ctrl-wet 7: flux 62332
568     (PID.TID 0000.0001) ctrl-wet 8: atmos 58276
569     (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
570     (PID.TID 0000.0001) ctrl-wet 13: global nvarlength for Nr = 15 91470
571     (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
572     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 1 2158 2004 2052 0
573     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 2 2155 2003 2049 0
574     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 3 2151 2000 2047 0
575     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 4 2143 1993 2040 0
576     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 5 2141 1993 2038 0
577     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 6 2139 1991 2036 0
578     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 7 2126 1977 2022 0
579     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 8 2108 1961 2003 0
580     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 9 2101 1954 1995 0
581     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 10 2086 1938 1981 0
582     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 11 2054 1903 1946 0
583     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 12 1998 1841 1879 0
584     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 13 1844 1650 1702 0
585     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 14 1180 951 996 0
586     (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 15 754 556 602 0
587     (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
588     (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
589     (PID.TID 0000.0001) ctrl-wet -------------------------------------------------
590     (PID.TID 0000.0001) ctrl_init: no. of control variables: 5
591     (PID.TID 0000.0001) ctrl_init: control vector length: 91470
592 jmc 1.7 (PID.TID 0000.0001) %MON fCori_max = 1.4265546244797E-04
593     (PID.TID 0000.0001) %MON fCori_min = -1.4265546244797E-04
594     (PID.TID 0000.0001) %MON fCori_mean = -1.2347858075529E-21
595     (PID.TID 0000.0001) %MON fCori_sd = 9.6599226301547E-05
596     (PID.TID 0000.0001) %MON fCoriG_max = 1.4151032568025E-04
597     (PID.TID 0000.0001) %MON fCoriG_min = -1.4362679550910E-04
598     (PID.TID 0000.0001) %MON fCoriG_mean = -3.5906698877274E-06
599     (PID.TID 0000.0001) %MON fCoriG_sd = 9.6548915696244E-05
600     (PID.TID 0000.0001) %MON fCoriCos_max = 1.4575362704741E-04
601     (PID.TID 0000.0001) %MON fCoriCos_min = 3.0322354601388E-05
602     (PID.TID 0000.0001) %MON fCoriCos_mean = 1.0288600773632E-04
603     (PID.TID 0000.0001) %MON fCoriCos_sd = 3.6784304327266E-05
604 jmc 1.1 (PID.TID 0000.0001) INI_CG2D: CG2D normalisation factor = 9.0544177972206531E-05
605     (PID.TID 0000.0001)
606     (PID.TID 0000.0001) // =======================================================
607     (PID.TID 0000.0001) // Model configuration
608     (PID.TID 0000.0001) // =======================================================
609     (PID.TID 0000.0001) //
610     (PID.TID 0000.0001) // "Physical" paramters ( PARM01 in namelist )
611     (PID.TID 0000.0001) //
612     (PID.TID 0000.0001) buoyancyRelation = /* Type of relation to get Buoyancy */
613     (PID.TID 0000.0001) 'OCEANIC'
614     (PID.TID 0000.0001) ;
615 jmc 1.6 (PID.TID 0000.0001) fluidIsAir = /* fluid major constituent is Air */
616 jmc 1.1 (PID.TID 0000.0001) F
617     (PID.TID 0000.0001) ;
618 jmc 1.6 (PID.TID 0000.0001) fluidIsWater = /* fluid major constituent is Water */
619 jmc 1.1 (PID.TID 0000.0001) T
620     (PID.TID 0000.0001) ;
621     (PID.TID 0000.0001) usingPCoords = /* use p (or p*) vertical coordinate */
622     (PID.TID 0000.0001) F
623     (PID.TID 0000.0001) ;
624     (PID.TID 0000.0001) usingZCoords = /* use z (or z*) vertical coordinate */
625     (PID.TID 0000.0001) T
626     (PID.TID 0000.0001) ;
627     (PID.TID 0000.0001) tRef = /* Reference temperature profile ( oC or K ) */
628     (PID.TID 0000.0001) 15 @ 2.000000000000000E+01 /* K = 1: 15 */
629     (PID.TID 0000.0001) ;
630     (PID.TID 0000.0001) sRef = /* Reference salinity profile ( psu ) */
631     (PID.TID 0000.0001) 15 @ 3.500000000000000E+01 /* K = 1: 15 */
632     (PID.TID 0000.0001) ;
633 jmc 1.6 (PID.TID 0000.0001) viscAh = /* Lateral eddy viscosity ( m^2/s ) */
634 jmc 1.1 (PID.TID 0000.0001) 5.000000000000000E+05
635     (PID.TID 0000.0001) ;
636 jmc 1.6 (PID.TID 0000.0001) viscAhMax = /* Maximum lateral eddy viscosity ( m^2/s ) */
637 jmc 1.1 (PID.TID 0000.0001) 1.000000000000000E+21
638     (PID.TID 0000.0001) ;
639 jmc 1.6 (PID.TID 0000.0001) viscAhGrid = /* Grid dependent lateral eddy viscosity ( non-dim. ) */
640 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
641     (PID.TID 0000.0001) ;
642 jmc 1.6 (PID.TID 0000.0001) useFullLeith = /* Use Full Form of Leith Viscosity on/off flag*/
643 jmc 1.1 (PID.TID 0000.0001) F
644     (PID.TID 0000.0001) ;
645 jmc 1.6 (PID.TID 0000.0001) useStrainTensionVisc= /* Use StrainTension Form of Viscous Operator flag*/
646 jmc 1.1 (PID.TID 0000.0001) F
647     (PID.TID 0000.0001) ;
648 jmc 1.6 (PID.TID 0000.0001) useAreaViscLength = /* Use area for visc length instead of geom. mean*/
649 jmc 1.1 (PID.TID 0000.0001) F
650     (PID.TID 0000.0001) ;
651     (PID.TID 0000.0001) viscC2leith = /* Leith harmonic visc. factor (on grad(vort),non-dim.) */
652     (PID.TID 0000.0001) 0.000000000000000E+00
653     (PID.TID 0000.0001) ;
654 jmc 1.6 (PID.TID 0000.0001) viscC2leithD = /* Leith harmonic viscosity factor (on grad(div),non-dim.)*/
655 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
656     (PID.TID 0000.0001) ;
657 jmc 1.6 (PID.TID 0000.0001) viscC2smag = /* Smagorinsky harmonic viscosity factor (non-dim.) */
658 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
659     (PID.TID 0000.0001) ;
660 jmc 1.6 (PID.TID 0000.0001) viscA4 = /* Lateral biharmonic viscosity ( m^4/s ) */
661 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
662     (PID.TID 0000.0001) ;
663 jmc 1.6 (PID.TID 0000.0001) viscA4Max = /* Maximum biharmonic viscosity ( m^2/s ) */
664 jmc 1.1 (PID.TID 0000.0001) 1.000000000000000E+21
665     (PID.TID 0000.0001) ;
666 jmc 1.6 (PID.TID 0000.0001) viscA4Grid = /* Grid dependent biharmonic viscosity ( non-dim. ) */
667 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
668     (PID.TID 0000.0001) ;
669 jmc 1.6 (PID.TID 0000.0001) viscC4leith = /* Leith biharm viscosity factor (on grad(vort), non-dim.)*/
670 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
671     (PID.TID 0000.0001) ;
672     (PID.TID 0000.0001) viscC4leithD = /* Leith biharm viscosity factor (on grad(div), non-dim.) */
673     (PID.TID 0000.0001) 0.000000000000000E+00
674     (PID.TID 0000.0001) ;
675     (PID.TID 0000.0001) viscC4Smag = /* Smagorinsky biharm viscosity factor (non-dim) */
676     (PID.TID 0000.0001) 0.000000000000000E+00
677     (PID.TID 0000.0001) ;
678     (PID.TID 0000.0001) no_slip_sides = /* Viscous BCs: No-slip sides */
679     (PID.TID 0000.0001) T
680     (PID.TID 0000.0001) ;
681     (PID.TID 0000.0001) sideDragFactor = /* side-drag scaling factor (non-dim) */
682     (PID.TID 0000.0001) 2.000000000000000E+00
683     (PID.TID 0000.0001) ;
684 jmc 1.6 (PID.TID 0000.0001) viscArNr = /* vertical profile of vertical viscosity ( m^2/s )*/
685     (PID.TID 0000.0001) 15 @ 1.000000000000000E-03 /* K = 1: 15 */
686 jmc 1.1 (PID.TID 0000.0001) ;
687     (PID.TID 0000.0001) no_slip_bottom = /* Viscous BCs: No-slip bottom */
688     (PID.TID 0000.0001) T
689     (PID.TID 0000.0001) ;
690 heimbach 1.3 (PID.TID 0000.0001) bottomDragLinear = /* linear bottom-drag coefficient ( m/s ) */
691 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
692     (PID.TID 0000.0001) ;
693 jmc 1.5 (PID.TID 0000.0001) bottomDragQuadratic = /* quadratic bottom-drag coefficient (-) */
694 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
695     (PID.TID 0000.0001) ;
696     (PID.TID 0000.0001) diffKhT = /* Laplacian diffusion of heat laterally ( m^2/s ) */
697     (PID.TID 0000.0001) 0.000000000000000E+00
698     (PID.TID 0000.0001) ;
699 heimbach 1.3 (PID.TID 0000.0001) diffK4T = /* Biharmonic diffusion of heat laterally ( m^4/s ) */
700 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
701     (PID.TID 0000.0001) ;
702     (PID.TID 0000.0001) diffKhS = /* Laplacian diffusion of salt laterally ( m^2/s ) */
703     (PID.TID 0000.0001) 0.000000000000000E+00
704     (PID.TID 0000.0001) ;
705 heimbach 1.3 (PID.TID 0000.0001) diffK4S = /* Biharmonic diffusion of salt laterally ( m^4/s ) */
706 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
707     (PID.TID 0000.0001) ;
708     (PID.TID 0000.0001) diffKrNrT = /* vertical profile of vertical diffusion of Temp ( m^2/s )*/
709     (PID.TID 0000.0001) 15 @ 3.000000000000000E-05 /* K = 1: 15 */
710     (PID.TID 0000.0001) ;
711     (PID.TID 0000.0001) diffKrNrS = /* vertical profile of vertical diffusion of Salt ( m^2/s )*/
712     (PID.TID 0000.0001) 15 @ 3.000000000000000E-05 /* K = 1: 15 */
713     (PID.TID 0000.0001) ;
714 heimbach 1.3 (PID.TID 0000.0001) diffKrBL79surf = /* Surface diffusion for Bryan and Lewis 79 ( m^2/s ) */
715 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
716     (PID.TID 0000.0001) ;
717 heimbach 1.3 (PID.TID 0000.0001) diffKrBL79deep = /* Deep diffusion for Bryan and Lewis 1979 ( m^2/s ) */
718 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
719     (PID.TID 0000.0001) ;
720 heimbach 1.3 (PID.TID 0000.0001) diffKrBL79scl = /* Depth scale for Bryan and Lewis 1979 ( m ) */
721 jmc 1.1 (PID.TID 0000.0001) 2.000000000000000E+02
722     (PID.TID 0000.0001) ;
723 heimbach 1.3 (PID.TID 0000.0001) diffKrBL79Ho = /* Turning depth for Bryan and Lewis 1979 ( m ) */
724 jmc 1.1 (PID.TID 0000.0001) -2.000000000000000E+03
725     (PID.TID 0000.0001) ;
726 jmc 1.6 (PID.TID 0000.0001) ivdc_kappa = /* Implicit Vertical Diffusivity for Convection ( m^2/s) */
727 jmc 1.7 (PID.TID 0000.0001) 1.000000000000000E+01
728 jmc 1.1 (PID.TID 0000.0001) ;
729     (PID.TID 0000.0001) hMixCriteria= /* Criteria for mixed-layer diagnostic */
730     (PID.TID 0000.0001) -8.000000000000000E-01
731     (PID.TID 0000.0001) ;
732 jmc 1.6 (PID.TID 0000.0001) dRhoSmall = /* Parameter for mixed-layer diagnostic */
733 jmc 1.5 (PID.TID 0000.0001) 1.000000000000000E-06
734     (PID.TID 0000.0001) ;
735     (PID.TID 0000.0001) hMixSmooth= /* Smoothing parameter for mixed-layer diagnostic */
736     (PID.TID 0000.0001) 0.000000000000000E+00
737     (PID.TID 0000.0001) ;
738 jmc 1.1 (PID.TID 0000.0001) eosType = /* Type of Equation of State */
739     (PID.TID 0000.0001) 'POLY3 '
740     (PID.TID 0000.0001) ;
741     (PID.TID 0000.0001) 1 13.497 32.600 24.548 -0.202 0.773 -0.005 -0.002 0.000 0.000 0.000 0.000 0.000
742     (PID.TID 0000.0001) 2 13.488 32.600 24.819 -0.204 0.773 -0.005 -0.002 0.000 0.000 0.000 0.000 0.000
743     (PID.TID 0000.0001) 3 13.476 32.600 25.201 -0.205 0.772 -0.005 -0.002 0.000 0.000 0.000 0.000 0.000
744 jmc 1.6 (PID.TID 0000.0001) 4 8.469 35.150 28.642 -0.163 0.784 -0.005 -0.002 -0.000 0.000 0.000 0.000 0.000
745     (PID.TID 0000.0001) 5 8.451 35.150 29.393 -0.167 0.782 -0.005 -0.002 -0.000 0.000 0.000 0.000 0.000
746 jmc 1.1 (PID.TID 0000.0001) 6 5.940 34.900 30.563 -0.144 0.786 -0.005 -0.002 -0.001 0.000 0.000 0.000 0.000
747 jmc 1.6 (PID.TID 0000.0001) 7 4.425 34.900 31.985 -0.134 0.788 -0.005 -0.003 -0.000 0.000 0.000 0.000 0.000
748     (PID.TID 0000.0001) 8 2.911 34.750 33.505 -0.125 0.789 -0.006 -0.003 -0.000 0.000 0.000 0.000 0.000
749 jmc 1.1 (PID.TID 0000.0001) 9 1.893 34.800 35.341 -0.123 0.789 -0.006 -0.003 0.000 0.000 0.000 0.000 0.000
750     (PID.TID 0000.0001) 10 1.367 34.800 37.299 -0.129 0.787 -0.006 -0.003 0.000 0.000 0.000 0.000 0.000
751     (PID.TID 0000.0001) 11 1.327 34.800 39.405 -0.141 0.783 -0.005 -0.003 0.000 0.000 0.000 0.000 0.000
752     (PID.TID 0000.0001) 12 0.791 34.800 41.788 -0.149 0.780 -0.005 -0.003 0.000 0.000 0.000 0.000 0.000
753     (PID.TID 0000.0001) 13 0.736 34.750 44.259 -0.164 0.776 -0.005 -0.002 0.000 0.000 0.000 0.000 0.000
754     (PID.TID 0000.0001) 14 0.186 34.750 47.044 -0.175 0.773 -0.005 -0.002 0.000 0.000 0.000 0.000 0.000
755     (PID.TID 0000.0001) 15 0.595 34.750 49.835 -0.196 0.767 -0.004 -0.002 0.000 0.000 0.000 0.000 0.000
756 jmc 1.9 (PID.TID 0000.0001) celsius2K = /* 0 degree Celsius converted to Kelvin ( K ) */
757     (PID.TID 0000.0001) 2.731500000000000E+02
758 jmc 1.1 (PID.TID 0000.0001) ;
759 jmc 1.9 (PID.TID 0000.0001) rhoConst = /* Reference density (Boussinesq) ( kg/m^3 ) */
760 jmc 1.1 (PID.TID 0000.0001) 1.035000000000000E+03
761     (PID.TID 0000.0001) ;
762     (PID.TID 0000.0001) rhoFacC = /* normalized Reference density @ cell-Center (-) */
763     (PID.TID 0000.0001) 15 @ 1.000000000000000E+00 /* K = 1: 15 */
764     (PID.TID 0000.0001) ;
765     (PID.TID 0000.0001) rhoFacF = /* normalized Reference density @ W-Interface (-) */
766     (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
767     (PID.TID 0000.0001) ;
768 jmc 1.9 (PID.TID 0000.0001) rhoConstFresh = /* Fresh-water reference density ( kg/m^3 ) */
769 jmc 1.7 (PID.TID 0000.0001) 1.000000000000000E+03
770 jmc 1.1 (PID.TID 0000.0001) ;
771     (PID.TID 0000.0001) gravity = /* Gravitational acceleration ( m/s^2 ) */
772     (PID.TID 0000.0001) 9.810000000000000E+00
773     (PID.TID 0000.0001) ;
774     (PID.TID 0000.0001) gBaro = /* Barotropic gravity ( m/s^2 ) */
775     (PID.TID 0000.0001) 9.810000000000000E+00
776     (PID.TID 0000.0001) ;
777     (PID.TID 0000.0001) rotationPeriod = /* Rotation Period ( s ) */
778 jmc 1.7 (PID.TID 0000.0001) 8.616400000000000E+04
779 jmc 1.1 (PID.TID 0000.0001) ;
780     (PID.TID 0000.0001) omega = /* Angular velocity ( rad/s ) */
781 jmc 1.7 (PID.TID 0000.0001) 7.292123516990375E-05
782 jmc 1.1 (PID.TID 0000.0001) ;
783     (PID.TID 0000.0001) f0 = /* Reference coriolis parameter ( 1/s ) */
784     (PID.TID 0000.0001) 1.000000000000000E-04
785     (PID.TID 0000.0001) ;
786     (PID.TID 0000.0001) beta = /* Beta ( 1/(m.s) ) */
787     (PID.TID 0000.0001) 9.999999999999999E-12
788     (PID.TID 0000.0001) ;
789 jmc 1.7 (PID.TID 0000.0001) fPrime = /* Second coriolis parameter ( 1/s ) */
790     (PID.TID 0000.0001) 0.000000000000000E+00
791     (PID.TID 0000.0001) ;
792 jmc 1.6 (PID.TID 0000.0001) rigidLid = /* Rigid lid on/off flag */
793     (PID.TID 0000.0001) F
794 jmc 1.1 (PID.TID 0000.0001) ;
795     (PID.TID 0000.0001) implicitFreeSurface = /* Implicit free surface on/off flag */
796     (PID.TID 0000.0001) T
797     (PID.TID 0000.0001) ;
798 jmc 1.6 (PID.TID 0000.0001) freeSurfFac = /* Implicit free surface factor */
799     (PID.TID 0000.0001) 1.000000000000000E+00
800 jmc 1.1 (PID.TID 0000.0001) ;
801 jmc 1.6 (PID.TID 0000.0001) implicSurfPress = /* Surface Pressure implicit factor (0-1)*/
802 jmc 1.1 (PID.TID 0000.0001) 1.000000000000000E+00
803     (PID.TID 0000.0001) ;
804 jmc 1.6 (PID.TID 0000.0001) implicDiv2Dflow = /* Barot. Flow Div. implicit factor (0-1)*/
805 jmc 1.1 (PID.TID 0000.0001) 1.000000000000000E+00
806     (PID.TID 0000.0001) ;
807 jmc 1.9 (PID.TID 0000.0001) uniformLin_PhiSurf = /* use uniform Bo_surf on/off flag*/
808     (PID.TID 0000.0001) T
809 jmc 1.1 (PID.TID 0000.0001) ;
810 jmc 1.9 (PID.TID 0000.0001) uniformFreeSurfLev = /* free-surface level-index is uniform */
811 jmc 1.1 (PID.TID 0000.0001) T
812     (PID.TID 0000.0001) ;
813 jmc 1.6 (PID.TID 0000.0001) hFacMin = /* minimum partial cell factor (hFac) */
814     (PID.TID 0000.0001) 5.000000000000000E-02
815     (PID.TID 0000.0001) ;
816     (PID.TID 0000.0001) hFacMinDr = /* minimum partial cell thickness ( m) */
817     (PID.TID 0000.0001) 5.000000000000000E-02
818     (PID.TID 0000.0001) ;
819 jmc 1.9 (PID.TID 0000.0001) exactConserv = /* Exact Volume Conservation on/off flag*/
820     (PID.TID 0000.0001) F
821     (PID.TID 0000.0001) ;
822     (PID.TID 0000.0001) linFSConserveTr = /* Tracer correction for Lin Free Surface on/off flag*/
823     (PID.TID 0000.0001) F
824     (PID.TID 0000.0001) ;
825 jmc 1.7 (PID.TID 0000.0001) nonlinFreeSurf = /* Non-linear Free Surf. options (-1,0,1,2,3)*/
826 jmc 1.1 (PID.TID 0000.0001) 0
827 jmc 1.7 (PID.TID 0000.0001) -1,0= Off ; 1,2,3= On, 2=+rescale gU,gV, 3=+update cg2d solv.
828 jmc 1.1 (PID.TID 0000.0001) ;
829     (PID.TID 0000.0001) hFacInf = /* lower threshold for hFac (nonlinFreeSurf only)*/
830     (PID.TID 0000.0001) 2.000000000000000E-01
831     (PID.TID 0000.0001) ;
832     (PID.TID 0000.0001) hFacSup = /* upper threshold for hFac (nonlinFreeSurf only)*/
833     (PID.TID 0000.0001) 2.000000000000000E+00
834     (PID.TID 0000.0001) ;
835 jmc 1.6 (PID.TID 0000.0001) select_rStar = /* r* Vertical coord. options (=0 r coord.; >0 uses r*)*/
836 jmc 1.4 (PID.TID 0000.0001) 0
837     (PID.TID 0000.0001) ;
838     (PID.TID 0000.0001) useRealFreshWaterFlux = /* Real Fresh Water Flux on/off flag*/
839 jmc 1.7 (PID.TID 0000.0001) F
840 jmc 1.1 (PID.TID 0000.0001) ;
841 jmc 1.4 (PID.TID 0000.0001) temp_EvPrRn = /* Temp. of Evap/Prec/R (UNSET=use local T)(oC)*/
842     (PID.TID 0000.0001) 1.234567000000000E+05
843     (PID.TID 0000.0001) ;
844 jmc 1.6 (PID.TID 0000.0001) salt_EvPrRn = /* Salin. of Evap/Prec/R (UNSET=use local S)(psu)*/
845 jmc 1.4 (PID.TID 0000.0001) 0.000000000000000E+00
846     (PID.TID 0000.0001) ;
847 jmc 1.9 (PID.TID 0000.0001) selectAddFluid = /* option for mass source/sink of fluid (=0: off) */
848     (PID.TID 0000.0001) 0
849     (PID.TID 0000.0001) ;
850 jmc 1.6 (PID.TID 0000.0001) temp_addMass = /* Temp. of addMass array (UNSET=use local T)(oC)*/
851     (PID.TID 0000.0001) 1.234567000000000E+05
852     (PID.TID 0000.0001) ;
853     (PID.TID 0000.0001) salt_addMass = /* Salin. of addMass array (UNSET=use local S)(psu)*/
854     (PID.TID 0000.0001) 0.000000000000000E+00
855     (PID.TID 0000.0001) ;
856     (PID.TID 0000.0001) convertFW2Salt = /* convert F.W. Flux to Salt Flux (-1=use local S)(psu)*/
857 jmc 1.7 (PID.TID 0000.0001) 3.500000000000000E+01
858 jmc 1.1 (PID.TID 0000.0001) ;
859     (PID.TID 0000.0001) use3Dsolver = /* use 3-D pressure solver on/off flag */
860     (PID.TID 0000.0001) F
861     (PID.TID 0000.0001) ;
862     (PID.TID 0000.0001) nonHydrostatic = /* Non-Hydrostatic on/off flag */
863     (PID.TID 0000.0001) F
864     (PID.TID 0000.0001) ;
865     (PID.TID 0000.0001) nh_Am2 = /* Non-Hydrostatic terms scaling factor */
866     (PID.TID 0000.0001) 1.000000000000000E+00
867     (PID.TID 0000.0001) ;
868 jmc 1.6 (PID.TID 0000.0001) implicitNHPress = /* Non-Hyd Pressure implicit factor (0-1)*/
869     (PID.TID 0000.0001) 1.000000000000000E+00
870     (PID.TID 0000.0001) ;
871     (PID.TID 0000.0001) selectNHfreeSurf = /* Non-Hyd (free-)Surface option */
872     (PID.TID 0000.0001) 0
873     (PID.TID 0000.0001) ;
874 jmc 1.1 (PID.TID 0000.0001) quasiHydrostatic = /* Quasi-Hydrostatic on/off flag */
875     (PID.TID 0000.0001) F
876     (PID.TID 0000.0001) ;
877 jmc 1.9 (PID.TID 0000.0001) calc_wVelocity = /* vertical velocity calculation on/off flag */
878     (PID.TID 0000.0001) T
879     (PID.TID 0000.0001) ;
880 jmc 1.1 (PID.TID 0000.0001) momStepping = /* Momentum equation on/off flag */
881     (PID.TID 0000.0001) T
882     (PID.TID 0000.0001) ;
883     (PID.TID 0000.0001) vectorInvariantMomentum= /* Vector-Invariant Momentum on/off */
884     (PID.TID 0000.0001) F
885     (PID.TID 0000.0001) ;
886     (PID.TID 0000.0001) momAdvection = /* Momentum advection on/off flag */
887     (PID.TID 0000.0001) T
888     (PID.TID 0000.0001) ;
889     (PID.TID 0000.0001) momViscosity = /* Momentum viscosity on/off flag */
890     (PID.TID 0000.0001) T
891     (PID.TID 0000.0001) ;
892 jmc 1.6 (PID.TID 0000.0001) momImplVertAdv= /* Momentum implicit vert. advection on/off*/
893 jmc 1.1 (PID.TID 0000.0001) F
894     (PID.TID 0000.0001) ;
895     (PID.TID 0000.0001) implicitViscosity = /* Implicit viscosity on/off flag */
896     (PID.TID 0000.0001) F
897     (PID.TID 0000.0001) ;
898     (PID.TID 0000.0001) metricTerms = /* metric-Terms on/off flag */
899     (PID.TID 0000.0001) T
900     (PID.TID 0000.0001) ;
901     (PID.TID 0000.0001) useNHMTerms = /* Non-Hydrostatic Metric-Terms on/off */
902     (PID.TID 0000.0001) F
903     (PID.TID 0000.0001) ;
904 jmc 1.7 (PID.TID 0000.0001) selectCoriMap = /* Coriolis Map options (0,1,2,3)*/
905     (PID.TID 0000.0001) 2
906     (PID.TID 0000.0001) 0= f-Plane ; 1= Beta-Plane ; 2= Spherical ; 3= read from file
907 jmc 1.1 (PID.TID 0000.0001) ;
908     (PID.TID 0000.0001) use3dCoriolis = /* 3-D Coriolis on/off flag */
909     (PID.TID 0000.0001) F
910     (PID.TID 0000.0001) ;
911     (PID.TID 0000.0001) useCoriolis = /* Coriolis on/off flag */
912     (PID.TID 0000.0001) T
913     (PID.TID 0000.0001) ;
914     (PID.TID 0000.0001) useCDscheme = /* CD scheme on/off flag */
915     (PID.TID 0000.0001) T
916     (PID.TID 0000.0001) ;
917 jmc 1.7 (PID.TID 0000.0001) useEnergyConservingCoriolis= /* Flx-Form Coriolis scheme flag */
918     (PID.TID 0000.0001) F
919     (PID.TID 0000.0001) ;
920 jmc 1.1 (PID.TID 0000.0001) useJamartWetPoints= /* Coriolis WetPoints method flag */
921     (PID.TID 0000.0001) F
922     (PID.TID 0000.0001) ;
923 jmc 1.7 (PID.TID 0000.0001) useJamartMomAdv= /* V.I Non-linear terms Jamart flag */
924 jmc 1.1 (PID.TID 0000.0001) F
925     (PID.TID 0000.0001) ;
926 jmc 1.7 (PID.TID 0000.0001) useAbsVorticity= /* V.I Works with f+zeta in Coriolis */
927 jmc 1.1 (PID.TID 0000.0001) F
928     (PID.TID 0000.0001) ;
929 jmc 1.7 (PID.TID 0000.0001) selectVortScheme= /* V.I Scheme selector for Vorticity-Term */
930 jmc 1.4 (PID.TID 0000.0001) 123456789
931     (PID.TID 0000.0001) = 0 : enstrophy (Shallow-Water Eq.) conserving scheme by Sadourny, JAS 75
932     (PID.TID 0000.0001) = 1 : same as 0 with modified hFac
933     (PID.TID 0000.0001) = 2 : energy conserving scheme (used by Sadourny in JAS 75 paper)
934     (PID.TID 0000.0001) = 3 : energy (general) and enstrophy (2D, nonDiv.) conserving scheme
935     (PID.TID 0000.0001) from Sadourny (Burridge & Haseler, ECMWF Rep.4, 1977)
936     (PID.TID 0000.0001) ;
937 jmc 1.7 (PID.TID 0000.0001) upwindVorticity= /* V.I Upwind bias vorticity flag */
938 jmc 1.1 (PID.TID 0000.0001) F
939     (PID.TID 0000.0001) ;
940 jmc 1.7 (PID.TID 0000.0001) highOrderVorticity= /* V.I High order vort. advect. flag */
941 jmc 1.1 (PID.TID 0000.0001) F
942     (PID.TID 0000.0001) ;
943 jmc 1.7 (PID.TID 0000.0001) upwindShear= /* V.I Upwind vertical Shear advection flag */
944 jmc 1.1 (PID.TID 0000.0001) F
945     (PID.TID 0000.0001) ;
946 jmc 1.7 (PID.TID 0000.0001) selectKEscheme= /* V.I Kinetic Energy scheme selector */
947 jmc 1.1 (PID.TID 0000.0001) 0
948     (PID.TID 0000.0001) ;
949     (PID.TID 0000.0001) momForcing = /* Momentum forcing on/off flag */
950     (PID.TID 0000.0001) T
951     (PID.TID 0000.0001) ;
952     (PID.TID 0000.0001) momPressureForcing = /* Momentum pressure term on/off flag */
953     (PID.TID 0000.0001) T
954     (PID.TID 0000.0001) ;
955     (PID.TID 0000.0001) implicitIntGravWave= /* Implicit Internal Gravity Wave flag */
956     (PID.TID 0000.0001) F
957     (PID.TID 0000.0001) ;
958 jmc 1.9 (PID.TID 0000.0001) staggerTimeStep = /* Stagger time stepping on/off flag */
959     (PID.TID 0000.0001) F
960     (PID.TID 0000.0001) ;
961     (PID.TID 0000.0001) doResetHFactors = /* reset thickness factors @ each time-step */
962 jmc 1.1 (PID.TID 0000.0001) F
963     (PID.TID 0000.0001) ;
964 jmc 1.9 (PID.TID 0000.0001) multiDimAdvection = /* enable/disable Multi-Dim Advection */
965 jmc 1.1 (PID.TID 0000.0001) T
966     (PID.TID 0000.0001) ;
967     (PID.TID 0000.0001) useMultiDimAdvec = /* Multi-Dim Advection is/is-not used */
968     (PID.TID 0000.0001) F
969     (PID.TID 0000.0001) ;
970 jmc 1.6 (PID.TID 0000.0001) implicitDiffusion = /* Implicit Diffusion on/off flag */
971 jmc 1.1 (PID.TID 0000.0001) T
972     (PID.TID 0000.0001) ;
973     (PID.TID 0000.0001) tempStepping = /* Temperature equation on/off flag */
974     (PID.TID 0000.0001) T
975     (PID.TID 0000.0001) ;
976 jmc 1.9 (PID.TID 0000.0001) tempAdvection = /* Temperature advection on/off flag */
977 jmc 1.1 (PID.TID 0000.0001) T
978     (PID.TID 0000.0001) ;
979 jmc 1.6 (PID.TID 0000.0001) tempImplVertAdv = /* Temp. implicit vert. advection on/off */
980 jmc 1.1 (PID.TID 0000.0001) F
981     (PID.TID 0000.0001) ;
982     (PID.TID 0000.0001) tempForcing = /* Temperature forcing on/off flag */
983     (PID.TID 0000.0001) T
984     (PID.TID 0000.0001) ;
985 jmc 1.9 (PID.TID 0000.0001) doThetaClimRelax = /* apply SST relaxation on/off flag */
986     (PID.TID 0000.0001) T
987     (PID.TID 0000.0001) ;
988 jmc 1.4 (PID.TID 0000.0001) tempIsActiveTr = /* Temp. is a dynamically Active Tracer */
989     (PID.TID 0000.0001) T
990     (PID.TID 0000.0001) ;
991 jmc 1.1 (PID.TID 0000.0001) saltStepping = /* Salinity equation on/off flag */
992     (PID.TID 0000.0001) T
993     (PID.TID 0000.0001) ;
994 jmc 1.9 (PID.TID 0000.0001) saltAdvection = /* Salinity advection on/off flag */
995 jmc 1.1 (PID.TID 0000.0001) T
996     (PID.TID 0000.0001) ;
997 jmc 1.6 (PID.TID 0000.0001) saltImplVertAdv = /* Sali. implicit vert. advection on/off */
998 jmc 1.1 (PID.TID 0000.0001) F
999     (PID.TID 0000.0001) ;
1000     (PID.TID 0000.0001) saltForcing = /* Salinity forcing on/off flag */
1001     (PID.TID 0000.0001) T
1002     (PID.TID 0000.0001) ;
1003 jmc 1.9 (PID.TID 0000.0001) doSaltClimRelax = /* apply SSS relaxation on/off flag */
1004     (PID.TID 0000.0001) T
1005     (PID.TID 0000.0001) ;
1006 jmc 1.4 (PID.TID 0000.0001) saltIsActiveTr = /* Salt is a dynamically Active Tracer */
1007     (PID.TID 0000.0001) T
1008     (PID.TID 0000.0001) ;
1009 jmc 1.1 (PID.TID 0000.0001) readBinaryPrec = /* Precision used for reading binary files */
1010     (PID.TID 0000.0001) 32
1011     (PID.TID 0000.0001) ;
1012     (PID.TID 0000.0001) writeBinaryPrec = /* Precision used for writing binary files */
1013     (PID.TID 0000.0001) 32
1014     (PID.TID 0000.0001) ;
1015     (PID.TID 0000.0001) globalFiles = /* write "global" (=not per tile) files */
1016     (PID.TID 0000.0001) F
1017     (PID.TID 0000.0001) ;
1018     (PID.TID 0000.0001) useSingleCpuIO = /* only master MPI process does I/O */
1019     (PID.TID 0000.0001) F
1020     (PID.TID 0000.0001) ;
1021 jmc 1.7 (PID.TID 0000.0001) /* debLev[*] : level of debug & auxiliary message printing */
1022     (PID.TID 0000.0001) debLevZero = 0 ; /* level of disabled aux. msg printing */
1023     (PID.TID 0000.0001) debLevA = 1 ; /* level of minimum aux. msg printing */
1024     (PID.TID 0000.0001) debLevB = 2 ; /* level of low aux. print (report read-file opening)*/
1025     (PID.TID 0000.0001) debLevC = 3 ; /* level of moderate debug prt (most pkgs debug msg) */
1026     (PID.TID 0000.0001) debLevD = 4 ; /* level of enhanced debug prt (add DEBUG_STATS prt) */
1027     (PID.TID 0000.0001) debLevE = 5 ; /* level of extensive debug printing */
1028     (PID.TID 0000.0001) debugLevel = /* select debug printing level */
1029 jmc 1.1 (PID.TID 0000.0001) 1
1030     (PID.TID 0000.0001) ;
1031     (PID.TID 0000.0001) //
1032     (PID.TID 0000.0001) // Elliptic solver(s) paramters ( PARM02 in namelist )
1033     (PID.TID 0000.0001) //
1034     (PID.TID 0000.0001) cg2dMaxIters = /* Upper limit on 2d con. grad iterations */
1035     (PID.TID 0000.0001) 500
1036     (PID.TID 0000.0001) ;
1037     (PID.TID 0000.0001) cg2dChkResFreq = /* 2d con. grad convergence test frequency */
1038     (PID.TID 0000.0001) 1
1039     (PID.TID 0000.0001) ;
1040 jmc 1.9 (PID.TID 0000.0001) cg2dUseMinResSol= /* use cg2d last-iter(=0) / min-resid.(=1) solution */
1041     (PID.TID 0000.0001) 0
1042     (PID.TID 0000.0001) ;
1043 jmc 1.1 (PID.TID 0000.0001) cg2dTargetResidual = /* 2d con. grad target residual */
1044     (PID.TID 0000.0001) 1.000000000000000E-13
1045     (PID.TID 0000.0001) ;
1046     (PID.TID 0000.0001) cg2dTargetResWunit = /* CG2d target residual [W units] */
1047     (PID.TID 0000.0001) -1.000000000000000E+00
1048     (PID.TID 0000.0001) ;
1049     (PID.TID 0000.0001) cg2dPreCondFreq = /* Freq. for updating cg2d preconditioner */
1050     (PID.TID 0000.0001) 1
1051     (PID.TID 0000.0001) ;
1052 jmc 1.6 (PID.TID 0000.0001) useSRCGSolver = /* use single reduction CG solver(s) */
1053     (PID.TID 0000.0001) F
1054     (PID.TID 0000.0001) ;
1055 jmc 1.7 (PID.TID 0000.0001) printResidualFreq = /* Freq. for printing CG residual */
1056     (PID.TID 0000.0001) 0
1057     (PID.TID 0000.0001) ;
1058 jmc 1.1 (PID.TID 0000.0001) //
1059     (PID.TID 0000.0001) // Time stepping paramters ( PARM03 in namelist )
1060     (PID.TID 0000.0001) //
1061 jmc 1.9 (PID.TID 0000.0001) deltaTMom = /* Momentum equation timestep ( s ) */
1062 jmc 1.1 (PID.TID 0000.0001) 1.200000000000000E+03
1063     (PID.TID 0000.0001) ;
1064 jmc 1.9 (PID.TID 0000.0001) deltaTFreeSurf = /* FreeSurface equation timestep ( s ) */
1065 jmc 1.7 (PID.TID 0000.0001) 4.320000000000000E+04
1066 jmc 1.1 (PID.TID 0000.0001) ;
1067 jmc 1.5 (PID.TID 0000.0001) dTtracerLev = /* Tracer equation timestep ( s ) */
1068 jmc 1.1 (PID.TID 0000.0001) 15 @ 4.320000000000000E+04 /* K = 1: 15 */
1069     (PID.TID 0000.0001) ;
1070     (PID.TID 0000.0001) deltaTClock = /* Model clock timestep ( s ) */
1071     (PID.TID 0000.0001) 4.320000000000000E+04
1072     (PID.TID 0000.0001) ;
1073     (PID.TID 0000.0001) cAdjFreq = /* Convective adjustment interval ( s ) */
1074     (PID.TID 0000.0001) 0.000000000000000E+00
1075     (PID.TID 0000.0001) ;
1076     (PID.TID 0000.0001) momForcingOutAB = /* =1: take Momentum Forcing out of Adams-Bash. stepping */
1077     (PID.TID 0000.0001) 0
1078     (PID.TID 0000.0001) ;
1079     (PID.TID 0000.0001) tracForcingOutAB = /* =1: take T,S,pTr Forcing out of Adams-Bash. stepping */
1080     (PID.TID 0000.0001) 0
1081     (PID.TID 0000.0001) ;
1082     (PID.TID 0000.0001) momDissip_In_AB = /* put Dissipation Tendency in Adams-Bash. stepping */
1083     (PID.TID 0000.0001) T
1084     (PID.TID 0000.0001) ;
1085     (PID.TID 0000.0001) doAB_onGtGs = /* apply AB on Tendencies (rather than on T,S)*/
1086     (PID.TID 0000.0001) T
1087     (PID.TID 0000.0001) ;
1088     (PID.TID 0000.0001) abEps = /* Adams-Bashforth-2 stabilizing weight */
1089     (PID.TID 0000.0001) 1.000000000000000E-01
1090     (PID.TID 0000.0001) ;
1091     (PID.TID 0000.0001) tauCD = /* CD coupling time-scale ( s ) */
1092     (PID.TID 0000.0001) 3.214280000000000E+05
1093     (PID.TID 0000.0001) ;
1094     (PID.TID 0000.0001) rCD = /* Normalised CD coupling parameter */
1095     (PID.TID 0000.0001) 9.962666600296178E-01
1096     (PID.TID 0000.0001) ;
1097 jmc 1.6 (PID.TID 0000.0001) epsAB_CD = /* AB-2 stabilizing weight for CD-scheme*/
1098     (PID.TID 0000.0001) 1.000000000000000E-01
1099     (PID.TID 0000.0001) ;
1100 heimbach 1.3 (PID.TID 0000.0001) pickupStrictlyMatch= /* stop if pickup do not strictly match */
1101     (PID.TID 0000.0001) T
1102     (PID.TID 0000.0001) ;
1103 jmc 1.7 (PID.TID 0000.0001) nIter0 = /* Run starting timestep number */
1104 heimbach 1.3 (PID.TID 0000.0001) 0
1105     (PID.TID 0000.0001) ;
1106 jmc 1.7 (PID.TID 0000.0001) nTimeSteps = /* Number of timesteps */
1107 jmc 1.9 (PID.TID 0000.0001) 20
1108 jmc 1.7 (PID.TID 0000.0001) ;
1109     (PID.TID 0000.0001) nEndIter = /* Run ending timestep number */
1110 jmc 1.9 (PID.TID 0000.0001) 20
1111 heimbach 1.3 (PID.TID 0000.0001) ;
1112 jmc 1.7 (PID.TID 0000.0001) baseTime = /* Model base time ( s ) */
1113 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
1114     (PID.TID 0000.0001) ;
1115 jmc 1.7 (PID.TID 0000.0001) startTime = /* Run start time ( s ) */
1116 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
1117     (PID.TID 0000.0001) ;
1118 jmc 1.7 (PID.TID 0000.0001) endTime = /* Integration ending time ( s ) */
1119 jmc 1.9 (PID.TID 0000.0001) 8.640000000000000E+05
1120 jmc 1.1 (PID.TID 0000.0001) ;
1121 jmc 1.7 (PID.TID 0000.0001) pChkPtFreq = /* Permanent restart/pickup file interval ( s ) */
1122 jmc 1.1 (PID.TID 0000.0001) 3.110400000000000E+08
1123     (PID.TID 0000.0001) ;
1124 jmc 1.7 (PID.TID 0000.0001) chkPtFreq = /* Rolling restart/pickup file interval ( s ) */
1125 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
1126     (PID.TID 0000.0001) ;
1127     (PID.TID 0000.0001) pickup_write_mdsio = /* Model IO flag. */
1128     (PID.TID 0000.0001) T
1129     (PID.TID 0000.0001) ;
1130     (PID.TID 0000.0001) pickup_read_mdsio = /* Model IO flag. */
1131     (PID.TID 0000.0001) T
1132     (PID.TID 0000.0001) ;
1133     (PID.TID 0000.0001) pickup_write_immed = /* Model IO flag. */
1134     (PID.TID 0000.0001) F
1135     (PID.TID 0000.0001) ;
1136 jmc 1.4 (PID.TID 0000.0001) writePickupAtEnd = /* Model IO flag. */
1137     (PID.TID 0000.0001) T
1138     (PID.TID 0000.0001) ;
1139 jmc 1.1 (PID.TID 0000.0001) dumpFreq = /* Model state write out interval ( s ). */
1140     (PID.TID 0000.0001) 8.640000000000000E+05
1141     (PID.TID 0000.0001) ;
1142     (PID.TID 0000.0001) dumpInitAndLast= /* write out Initial & Last iter. model state */
1143     (PID.TID 0000.0001) T
1144     (PID.TID 0000.0001) ;
1145     (PID.TID 0000.0001) snapshot_mdsio = /* Model IO flag. */
1146     (PID.TID 0000.0001) T
1147     (PID.TID 0000.0001) ;
1148     (PID.TID 0000.0001) monitorFreq = /* Monitor output interval ( s ). */
1149 jmc 1.7 (PID.TID 0000.0001) 2.592000000000000E+06
1150 jmc 1.1 (PID.TID 0000.0001) ;
1151 heimbach 1.3 (PID.TID 0000.0001) monitorSelect = /* select group of variables to monitor */
1152     (PID.TID 0000.0001) 3
1153     (PID.TID 0000.0001) ;
1154 jmc 1.1 (PID.TID 0000.0001) monitor_stdio = /* Model IO flag. */
1155     (PID.TID 0000.0001) T
1156     (PID.TID 0000.0001) ;
1157     (PID.TID 0000.0001) externForcingPeriod = /* forcing period (s) */
1158 jmc 1.7 (PID.TID 0000.0001) 0.000000000000000E+00
1159 jmc 1.1 (PID.TID 0000.0001) ;
1160     (PID.TID 0000.0001) externForcingCycle = /* period of the cyle (s). */
1161 jmc 1.7 (PID.TID 0000.0001) 0.000000000000000E+00
1162 jmc 1.1 (PID.TID 0000.0001) ;
1163     (PID.TID 0000.0001) tauThetaClimRelax = /* relaxation time scale (s) */
1164 jmc 1.8 (PID.TID 0000.0001) 5.184000000000000E+06
1165 jmc 1.1 (PID.TID 0000.0001) ;
1166     (PID.TID 0000.0001) tauSaltClimRelax = /* relaxation time scale (s) */
1167     (PID.TID 0000.0001) 1.555200000000000E+07
1168     (PID.TID 0000.0001) ;
1169     (PID.TID 0000.0001) latBandClimRelax = /* max. Lat. where relaxation */
1170     (PID.TID 0000.0001) 1.800000000000000E+02
1171     (PID.TID 0000.0001) ;
1172     (PID.TID 0000.0001) //
1173     (PID.TID 0000.0001) // Gridding paramters ( PARM04 in namelist )
1174     (PID.TID 0000.0001) //
1175     (PID.TID 0000.0001) usingCartesianGrid = /* Cartesian coordinates flag ( True/False ) */
1176     (PID.TID 0000.0001) F
1177     (PID.TID 0000.0001) ;
1178     (PID.TID 0000.0001) usingCylindricalGrid = /* Cylindrical coordinates flag ( True/False ) */
1179     (PID.TID 0000.0001) F
1180     (PID.TID 0000.0001) ;
1181     (PID.TID 0000.0001) usingSphericalPolarGrid = /* Spherical coordinates flag ( True/False ) */
1182     (PID.TID 0000.0001) T
1183     (PID.TID 0000.0001) ;
1184     (PID.TID 0000.0001) usingCurvilinearGrid = /* Curvilinear coordinates flag ( True/False ) */
1185     (PID.TID 0000.0001) F
1186     (PID.TID 0000.0001) ;
1187 jmc 1.6 (PID.TID 0000.0001) selectSigmaCoord = /* Hybrid-Sigma Vert. Coordinate option */
1188     (PID.TID 0000.0001) 0
1189     (PID.TID 0000.0001) ;
1190 heimbach 1.3 (PID.TID 0000.0001) Ro_SeaLevel = /* r(1) ( units of r == m ) */
1191 jmc 1.1 (PID.TID 0000.0001) 0.000000000000000E+00
1192     (PID.TID 0000.0001) ;
1193 jmc 1.6 (PID.TID 0000.0001) rSigmaBnd = /* r/sigma transition ( units of r == m ) */
1194     (PID.TID 0000.0001) 1.234567000000000E+05
1195     (PID.TID 0000.0001) ;
1196 jmc 1.1 (PID.TID 0000.0001) rkSign = /* index orientation relative to vertical coordinate */
1197     (PID.TID 0000.0001) -1.000000000000000E+00
1198     (PID.TID 0000.0001) ;
1199     (PID.TID 0000.0001) gravitySign = /* gravity orientation relative to vertical coordinate */
1200     (PID.TID 0000.0001) -1.000000000000000E+00
1201     (PID.TID 0000.0001) ;
1202     (PID.TID 0000.0001) mass2rUnit = /* convert mass per unit area [kg/m2] to r-units [m] */
1203     (PID.TID 0000.0001) 9.661835748792270E-04
1204     (PID.TID 0000.0001) ;
1205     (PID.TID 0000.0001) rUnit2mass = /* convert r-units [m] to mass per unit area [kg/m2] */
1206     (PID.TID 0000.0001) 1.035000000000000E+03
1207     (PID.TID 0000.0001) ;
1208     (PID.TID 0000.0001) drC = /* C spacing ( units of r ) */
1209     (PID.TID 0000.0001) 2.500000000000000E+01, /* K = 1 */
1210     (PID.TID 0000.0001) 6.000000000000000E+01, /* K = 2 */
1211     (PID.TID 0000.0001) 8.500000000000000E+01, /* K = 3 */
1212     (PID.TID 0000.0001) 1.200000000000000E+02, /* K = 4 */
1213     (PID.TID 0000.0001) 1.650000000000000E+02, /* K = 5 */
1214     (PID.TID 0000.0001) 2.150000000000000E+02, /* K = 6 */
1215     (PID.TID 0000.0001) 2.650000000000000E+02, /* K = 7 */
1216     (PID.TID 0000.0001) 3.150000000000000E+02, /* K = 8 */
1217     (PID.TID 0000.0001) 3.650000000000000E+02, /* K = 9 */
1218     (PID.TID 0000.0001) 4.150000000000000E+02, /* K = 10 */
1219     (PID.TID 0000.0001) 4.650000000000000E+02, /* K = 11 */
1220     (PID.TID 0000.0001) 5.150000000000000E+02, /* K = 12 */
1221     (PID.TID 0000.0001) 5.650000000000000E+02, /* K = 13 */
1222     (PID.TID 0000.0001) 6.150000000000000E+02, /* K = 14 */
1223     (PID.TID 0000.0001) 6.650000000000000E+02 /* K = 15 */
1224     (PID.TID 0000.0001) ;
1225     (PID.TID 0000.0001) drF = /* W spacing ( units of r ) */
1226     (PID.TID 0000.0001) 5.000000000000000E+01, /* K = 1 */
1227     (PID.TID 0000.0001) 7.000000000000000E+01, /* K = 2 */
1228     (PID.TID 0000.0001) 1.000000000000000E+02, /* K = 3 */
1229     (PID.TID 0000.0001) 1.400000000000000E+02, /* K = 4 */
1230     (PID.TID 0000.0001) 1.900000000000000E+02, /* K = 5 */
1231     (PID.TID 0000.0001) 2.400000000000000E+02, /* K = 6 */
1232     (PID.TID 0000.0001) 2.900000000000000E+02, /* K = 7 */
1233     (PID.TID 0000.0001) 3.400000000000000E+02, /* K = 8 */
1234     (PID.TID 0000.0001) 3.900000000000000E+02, /* K = 9 */
1235     (PID.TID 0000.0001) 4.400000000000000E+02, /* K = 10 */
1236     (PID.TID 0000.0001) 4.900000000000000E+02, /* K = 11 */
1237     (PID.TID 0000.0001) 5.400000000000000E+02, /* K = 12 */
1238     (PID.TID 0000.0001) 5.900000000000000E+02, /* K = 13 */
1239     (PID.TID 0000.0001) 6.400000000000000E+02, /* K = 14 */
1240     (PID.TID 0000.0001) 6.900000000000000E+02 /* K = 15 */
1241     (PID.TID 0000.0001) ;
1242 jmc 1.5 (PID.TID 0000.0001) delX = /* U spacing ( m - cartesian, degrees - spherical ) */
1243 jmc 1.1 (PID.TID 0000.0001) 90 @ 4.000000000000000E+00 /* I = 1: 90 */
1244     (PID.TID 0000.0001) ;
1245 jmc 1.5 (PID.TID 0000.0001) delY = /* V spacing ( m - cartesian, degrees - spherical ) */
1246 jmc 1.1 (PID.TID 0000.0001) 40 @ 4.000000000000000E+00 /* J = 1: 40 */
1247     (PID.TID 0000.0001) ;
1248 jmc 1.9 (PID.TID 0000.0001) xgOrigin = /* X-axis origin of West edge (cartesian: m, lat-lon: deg) */
1249 jmc 1.5 (PID.TID 0000.0001) 0.000000000000000E+00
1250     (PID.TID 0000.0001) ;
1251 jmc 1.9 (PID.TID 0000.0001) ygOrigin = /* Y-axis origin of South edge (cartesian: m, lat-lon: deg) */
1252 jmc 1.1 (PID.TID 0000.0001) -8.000000000000000E+01
1253     (PID.TID 0000.0001) ;
1254     (PID.TID 0000.0001) rSphere = /* Radius ( ignored - cartesian, m - spherical ) */
1255     (PID.TID 0000.0001) 6.370000000000000E+06
1256     (PID.TID 0000.0001) ;
1257     (PID.TID 0000.0001) deepAtmosphere = /* Deep/Shallow Atmosphere flag (True/False) */
1258     (PID.TID 0000.0001) F
1259     (PID.TID 0000.0001) ;
1260 jmc 1.5 (PID.TID 0000.0001) xC = /* xC(:,1,:,1) : P-point X coord ( deg. or m if cartesian) */
1261 jmc 1.1 (PID.TID 0000.0001) 2.000000000000000E+00, /* I = 1 */
1262     (PID.TID 0000.0001) 6.000000000000000E+00, /* I = 2 */
1263     (PID.TID 0000.0001) 1.000000000000000E+01, /* I = 3 */
1264 jmc 1.5 (PID.TID 0000.0001) . . .
1265 jmc 1.1 (PID.TID 0000.0001) 8.200000000000000E+01, /* I = 21 */
1266     (PID.TID 0000.0001) 8.600000000000000E+01, /* I = 22 */
1267     (PID.TID 0000.0001) 9.000000000000000E+01, /* I = 23 */
1268     (PID.TID 0000.0001) 9.400000000000000E+01, /* I = 24 */
1269     (PID.TID 0000.0001) 9.800000000000000E+01, /* I = 25 */
1270     (PID.TID 0000.0001) 1.020000000000000E+02, /* I = 26 */
1271 jmc 1.5 (PID.TID 0000.0001) . . .
1272 jmc 1.1 (PID.TID 0000.0001) 1.700000000000000E+02, /* I = 43 */
1273     (PID.TID 0000.0001) 1.740000000000000E+02, /* I = 44 */
1274     (PID.TID 0000.0001) 1.780000000000000E+02, /* I = 45 */
1275     (PID.TID 0000.0001) 1.820000000000000E+02, /* I = 46 */
1276     (PID.TID 0000.0001) 1.860000000000000E+02, /* I = 47 */
1277     (PID.TID 0000.0001) 1.900000000000000E+02, /* I = 48 */
1278 jmc 1.5 (PID.TID 0000.0001) . . .
1279 jmc 1.1 (PID.TID 0000.0001) 2.580000000000000E+02, /* I = 65 */
1280     (PID.TID 0000.0001) 2.620000000000000E+02, /* I = 66 */
1281     (PID.TID 0000.0001) 2.660000000000000E+02, /* I = 67 */
1282     (PID.TID 0000.0001) 2.700000000000000E+02, /* I = 68 */
1283     (PID.TID 0000.0001) 2.740000000000000E+02, /* I = 69 */
1284     (PID.TID 0000.0001) 2.780000000000000E+02, /* I = 70 */
1285 jmc 1.5 (PID.TID 0000.0001) . . .
1286 jmc 1.1 (PID.TID 0000.0001) 3.500000000000000E+02, /* I = 88 */
1287     (PID.TID 0000.0001) 3.540000000000000E+02, /* I = 89 */
1288     (PID.TID 0000.0001) 3.580000000000000E+02 /* I = 90 */
1289     (PID.TID 0000.0001) ;
1290 jmc 1.5 (PID.TID 0000.0001) yC = /* yC(1,:,1,:) : P-point Y coord ( deg. or m if cartesian) */
1291 jmc 1.1 (PID.TID 0000.0001) -7.800000000000000E+01, /* J = 1 */
1292     (PID.TID 0000.0001) -7.400000000000000E+01, /* J = 2 */
1293     (PID.TID 0000.0001) -7.000000000000000E+01, /* J = 3 */
1294     (PID.TID 0000.0001) -6.600000000000000E+01, /* J = 4 */
1295     (PID.TID 0000.0001) -6.200000000000000E+01, /* J = 5 */
1296     (PID.TID 0000.0001) -5.800000000000000E+01, /* J = 6 */
1297     (PID.TID 0000.0001) -5.400000000000000E+01, /* J = 7 */
1298     (PID.TID 0000.0001) -5.000000000000000E+01, /* J = 8 */
1299     (PID.TID 0000.0001) -4.600000000000000E+01, /* J = 9 */
1300     (PID.TID 0000.0001) -4.200000000000000E+01, /* J = 10 */
1301     (PID.TID 0000.0001) -3.800000000000000E+01, /* J = 11 */
1302     (PID.TID 0000.0001) -3.400000000000000E+01, /* J = 12 */
1303     (PID.TID 0000.0001) -3.000000000000000E+01, /* J = 13 */
1304     (PID.TID 0000.0001) -2.600000000000000E+01, /* J = 14 */
1305     (PID.TID 0000.0001) -2.200000000000000E+01, /* J = 15 */
1306     (PID.TID 0000.0001) -1.800000000000000E+01, /* J = 16 */
1307     (PID.TID 0000.0001) -1.400000000000000E+01, /* J = 17 */
1308     (PID.TID 0000.0001) -1.000000000000000E+01, /* J = 18 */
1309     (PID.TID 0000.0001) -6.000000000000000E+00, /* J = 19 */
1310     (PID.TID 0000.0001) -2.000000000000000E+00, /* J = 20 */
1311     (PID.TID 0000.0001) 2.000000000000000E+00, /* J = 21 */
1312     (PID.TID 0000.0001) 6.000000000000000E+00, /* J = 22 */
1313     (PID.TID 0000.0001) 1.000000000000000E+01, /* J = 23 */
1314     (PID.TID 0000.0001) 1.400000000000000E+01, /* J = 24 */
1315     (PID.TID 0000.0001) 1.800000000000000E+01, /* J = 25 */
1316     (PID.TID 0000.0001) 2.200000000000000E+01, /* J = 26 */
1317     (PID.TID 0000.0001) 2.600000000000000E+01, /* J = 27 */
1318     (PID.TID 0000.0001) 3.000000000000000E+01, /* J = 28 */
1319     (PID.TID 0000.0001) 3.400000000000000E+01, /* J = 29 */
1320     (PID.TID 0000.0001) 3.800000000000000E+01, /* J = 30 */
1321     (PID.TID 0000.0001) 4.200000000000000E+01, /* J = 31 */
1322     (PID.TID 0000.0001) 4.600000000000000E+01, /* J = 32 */
1323     (PID.TID 0000.0001) 5.000000000000000E+01, /* J = 33 */
1324     (PID.TID 0000.0001) 5.400000000000000E+01, /* J = 34 */
1325     (PID.TID 0000.0001) 5.800000000000000E+01, /* J = 35 */
1326     (PID.TID 0000.0001) 6.200000000000000E+01, /* J = 36 */
1327     (PID.TID 0000.0001) 6.600000000000000E+01, /* J = 37 */
1328     (PID.TID 0000.0001) 7.000000000000000E+01, /* J = 38 */
1329     (PID.TID 0000.0001) 7.400000000000000E+01, /* J = 39 */
1330     (PID.TID 0000.0001) 7.800000000000000E+01 /* J = 40 */
1331     (PID.TID 0000.0001) ;
1332 jmc 1.5 (PID.TID 0000.0001) rcoord = /* P-point R coordinate ( units of r ) */
1333 jmc 1.1 (PID.TID 0000.0001) -2.500000000000000E+01, /* K = 1 */
1334     (PID.TID 0000.0001) -8.500000000000000E+01, /* K = 2 */
1335     (PID.TID 0000.0001) -1.700000000000000E+02, /* K = 3 */
1336     (PID.TID 0000.0001) -2.900000000000000E+02, /* K = 4 */
1337     (PID.TID 0000.0001) -4.550000000000000E+02, /* K = 5 */
1338     (PID.TID 0000.0001) -6.700000000000000E+02, /* K = 6 */
1339     (PID.TID 0000.0001) -9.350000000000000E+02, /* K = 7 */
1340     (PID.TID 0000.0001) -1.250000000000000E+03, /* K = 8 */
1341     (PID.TID 0000.0001) -1.615000000000000E+03, /* K = 9 */
1342     (PID.TID 0000.0001) -2.030000000000000E+03, /* K = 10 */
1343     (PID.TID 0000.0001) -2.495000000000000E+03, /* K = 11 */
1344     (PID.TID 0000.0001) -3.010000000000000E+03, /* K = 12 */
1345     (PID.TID 0000.0001) -3.575000000000000E+03, /* K = 13 */
1346     (PID.TID 0000.0001) -4.190000000000000E+03, /* K = 14 */
1347     (PID.TID 0000.0001) -4.855000000000000E+03 /* K = 15 */
1348     (PID.TID 0000.0001) ;
1349     (PID.TID 0000.0001) rF = /* W-Interf. R coordinate ( units of r ) */
1350     (PID.TID 0000.0001) 0.000000000000000E+00, /* K = 1 */
1351     (PID.TID 0000.0001) -5.000000000000000E+01, /* K = 2 */
1352     (PID.TID 0000.0001) -1.200000000000000E+02, /* K = 3 */
1353     (PID.TID 0000.0001) -2.200000000000000E+02, /* K = 4 */
1354     (PID.TID 0000.0001) -3.600000000000000E+02, /* K = 5 */
1355     (PID.TID 0000.0001) -5.500000000000000E+02, /* K = 6 */
1356     (PID.TID 0000.0001) -7.900000000000000E+02, /* K = 7 */
1357     (PID.TID 0000.0001) -1.080000000000000E+03, /* K = 8 */
1358     (PID.TID 0000.0001) -1.420000000000000E+03, /* K = 9 */
1359     (PID.TID 0000.0001) -1.810000000000000E+03, /* K = 10 */
1360     (PID.TID 0000.0001) -2.250000000000000E+03, /* K = 11 */
1361     (PID.TID 0000.0001) -2.740000000000000E+03, /* K = 12 */
1362     (PID.TID 0000.0001) -3.280000000000000E+03, /* K = 13 */
1363     (PID.TID 0000.0001) -3.870000000000000E+03, /* K = 14 */
1364     (PID.TID 0000.0001) -4.510000000000000E+03, /* K = 15 */
1365     (PID.TID 0000.0001) -5.200000000000000E+03 /* K = 16 */
1366     (PID.TID 0000.0001) ;
1367     (PID.TID 0000.0001) deepFacC = /* deep-model grid factor @ cell-Center (-) */
1368     (PID.TID 0000.0001) 15 @ 1.000000000000000E+00 /* K = 1: 15 */
1369     (PID.TID 0000.0001) ;
1370     (PID.TID 0000.0001) deepFacF = /* deep-model grid factor @ W-Interface (-) */
1371     (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
1372     (PID.TID 0000.0001) ;
1373     (PID.TID 0000.0001) rVel2wUnit = /* convert units: rVel -> wSpeed (=1 if z-coord)*/
1374     (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
1375     (PID.TID 0000.0001) ;
1376     (PID.TID 0000.0001) wUnit2rVel = /* convert units: wSpeed -> rVel (=1 if z-coord)*/
1377     (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */
1378     (PID.TID 0000.0001) ;
1379 jmc 1.7 (PID.TID 0000.0001) dBdrRef = /* Vertical grad. of reference buoyancy [(m/s/r)^2] */
1380 jmc 1.1 (PID.TID 0000.0001) 15 @ 0.000000000000000E+00 /* K = 1: 15 */
1381     (PID.TID 0000.0001) ;
1382 jmc 1.4 (PID.TID 0000.0001) rotateGrid = /* use rotated grid ( True/False ) */
1383     (PID.TID 0000.0001) F
1384     (PID.TID 0000.0001) ;
1385     (PID.TID 0000.0001) phiEuler = /* Euler angle, rotation about original z-coordinate [rad] */
1386     (PID.TID 0000.0001) 0.000000000000000E+00
1387     (PID.TID 0000.0001) ;
1388     (PID.TID 0000.0001) thetaEuler = /* Euler angle, rotation about new x-coordinate [rad] */
1389     (PID.TID 0000.0001) 0.000000000000000E+00
1390     (PID.TID 0000.0001) ;
1391     (PID.TID 0000.0001) psiEuler = /* Euler angle, rotation about new z-coordinate [rad] */
1392     (PID.TID 0000.0001) 0.000000000000000E+00
1393     (PID.TID 0000.0001) ;
1394 jmc 1.1 (PID.TID 0000.0001) dxF = /* dxF(:,1,:,1) ( units: m ) */
1395     (PID.TID 0000.0001) 90 @ 9.246038586187513E+04 /* I = 1: 90 */
1396     (PID.TID 0000.0001) ;
1397     (PID.TID 0000.0001) dxF = /* dxF(1,:,1,:) ( units: m ) */
1398     (PID.TID 0000.0001) 9.246038586187513E+04, /* J = 1 */
1399     (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 2 */
1400     (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 3 */
1401     (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 4 */
1402     (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 5 */
1403     (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 6 */
1404     (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 7 */
1405     (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 8 */
1406     (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 9 */
1407     (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 10 */
1408     (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 11 */
1409     (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 12 */
1410     (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 13 */
1411     (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 14 */
1412     (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 15 */
1413     (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 16 */
1414     (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 17 */
1415     (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 18 */
1416     (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 19 */
1417     (PID.TID 0000.0001) 2 @ 4.444389881567502E+05, /* J = 20: 21 */
1418     (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 22 */
1419     (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 23 */
1420     (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 24 */
1421     (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 25 */
1422     (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 26 */
1423     (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 27 */
1424     (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 28 */
1425     (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 29 */
1426     (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 30 */
1427     (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 31 */
1428     (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 32 */
1429     (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 33 */
1430     (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 34 */
1431     (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 35 */
1432     (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 36 */
1433     (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 37 */
1434     (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 38 */
1435     (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 39 */
1436     (PID.TID 0000.0001) 9.246038586187513E+04 /* J = 40 */
1437     (PID.TID 0000.0001) ;
1438     (PID.TID 0000.0001) dyF = /* dyF(:,1,:,1) ( units: m ) */
1439     (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1440     (PID.TID 0000.0001) ;
1441     (PID.TID 0000.0001) dyF = /* dyF(1,:,1,:) ( units: m ) */
1442     (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1443     (PID.TID 0000.0001) ;
1444     (PID.TID 0000.0001) dxG = /* dxG(:,1,:,1) ( units: m ) */
1445     (PID.TID 0000.0001) 90 @ 7.722306258078101E+04 /* I = 1: 90 */
1446     (PID.TID 0000.0001) ;
1447     (PID.TID 0000.0001) dxG = /* dxG(1,:,1,:) ( units: m ) */
1448     (PID.TID 0000.0001) 7.722306258078101E+04, /* J = 1 */
1449     (PID.TID 0000.0001) 1.075850604052271E+05, /* J = 2 */
1450     (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 3 */
1451     (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 4 */
1452     (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 5 */
1453     (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 6 */
1454     (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 7 */
1455     (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 8 */
1456     (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 9 */
1457     (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 10 */
1458     (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 11 */
1459     (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 12 */
1460     (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 13 */
1461     (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 14 */
1462     (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 15 */
1463     (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 16 */
1464     (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 17 */
1465     (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 18 */
1466     (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 19 */
1467     (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 20 */
1468     (PID.TID 0000.0001) 4.447098934081552E+05, /* J = 21 */
1469     (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 22 */
1470     (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 23 */
1471     (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 24 */
1472     (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 25 */
1473     (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 26 */
1474     (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 27 */
1475     (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 28 */
1476     (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 29 */
1477     (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 30 */
1478     (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 31 */
1479     (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 32 */
1480     (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 33 */
1481     (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 34 */
1482     (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 35 */
1483     (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 36 */
1484     (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 37 */
1485     (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 38 */
1486     (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 39 */
1487     (PID.TID 0000.0001) 1.075850604052271E+05 /* J = 40 */
1488     (PID.TID 0000.0001) ;
1489     (PID.TID 0000.0001) dyG = /* dyG(:,1,:,1) ( units: m ) */
1490     (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1491     (PID.TID 0000.0001) ;
1492     (PID.TID 0000.0001) dyG = /* dyG(1,:,1,:) ( units: m ) */
1493     (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1494     (PID.TID 0000.0001) ;
1495     (PID.TID 0000.0001) dxC = /* dxC(:,1,:,1) ( units: m ) */
1496     (PID.TID 0000.0001) 90 @ 9.246038586187513E+04 /* I = 1: 90 */
1497     (PID.TID 0000.0001) ;
1498     (PID.TID 0000.0001) dxC = /* dxC(1,:,1,:) ( units: m ) */
1499     (PID.TID 0000.0001) 9.246038586187513E+04, /* J = 1 */
1500     (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 2 */
1501     (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 3 */
1502     (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 4 */
1503     (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 5 */
1504     (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 6 */
1505     (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 7 */
1506     (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 8 */
1507     (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 9 */
1508     (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 10 */
1509     (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 11 */
1510     (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 12 */
1511     (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 13 */
1512     (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 14 */
1513     (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 15 */
1514     (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 16 */
1515     (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 17 */
1516     (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 18 */
1517     (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 19 */
1518     (PID.TID 0000.0001) 2 @ 4.444389881567502E+05, /* J = 20: 21 */
1519     (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 22 */
1520     (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 23 */
1521     (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 24 */
1522     (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 25 */
1523     (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 26 */
1524     (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 27 */
1525     (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 28 */
1526     (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 29 */
1527     (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 30 */
1528     (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 31 */
1529     (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 32 */
1530     (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 33 */
1531     (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 34 */
1532     (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 35 */
1533     (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 36 */
1534     (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 37 */
1535     (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 38 */
1536     (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 39 */
1537     (PID.TID 0000.0001) 9.246038586187513E+04 /* J = 40 */
1538     (PID.TID 0000.0001) ;
1539     (PID.TID 0000.0001) dyC = /* dyC(:,1,:,1) ( units: m ) */
1540     (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1541     (PID.TID 0000.0001) ;
1542     (PID.TID 0000.0001) dyC = /* dyC(1,:,1,:) ( units: m ) */
1543     (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1544     (PID.TID 0000.0001) ;
1545     (PID.TID 0000.0001) dxV = /* dxV(:,1,:,1) ( units: m ) */
1546     (PID.TID 0000.0001) 90 @ 7.722306258078101E+04 /* I = 1: 90 */
1547     (PID.TID 0000.0001) ;
1548     (PID.TID 0000.0001) dxV = /* dxV(1,:,1,:) ( units: m ) */
1549     (PID.TID 0000.0001) 7.722306258078101E+04, /* J = 1 */
1550     (PID.TID 0000.0001) 1.075850604052271E+05, /* J = 2 */
1551     (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 3 */
1552     (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 4 */
1553     (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 5 */
1554     (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 6 */
1555     (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 7 */
1556     (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 8 */
1557     (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 9 */
1558     (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 10 */
1559     (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 11 */
1560     (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 12 */
1561     (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 13 */
1562     (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 14 */
1563     (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 15 */
1564     (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 16 */
1565     (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 17 */
1566     (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 18 */
1567     (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 19 */
1568     (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 20 */
1569     (PID.TID 0000.0001) 4.447098934081552E+05, /* J = 21 */
1570     (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 22 */
1571     (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 23 */
1572     (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 24 */
1573     (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 25 */
1574     (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 26 */
1575     (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 27 */
1576     (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 28 */
1577     (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 29 */
1578     (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 30 */
1579     (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 31 */
1580     (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 32 */
1581     (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 33 */
1582     (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 34 */
1583     (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 35 */
1584     (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 36 */
1585     (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 37 */
1586     (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 38 */
1587     (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 39 */
1588     (PID.TID 0000.0001) 1.075850604052271E+05 /* J = 40 */
1589     (PID.TID 0000.0001) ;
1590     (PID.TID 0000.0001) dyU = /* dyU(:,1,:,1) ( units: m ) */
1591     (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */
1592     (PID.TID 0000.0001) ;
1593     (PID.TID 0000.0001) dyU = /* dyU(1,:,1,:) ( units: m ) */
1594     (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */
1595     (PID.TID 0000.0001) ;
1596     (PID.TID 0000.0001) rA = /* rA (:,1,:,1) ( units: m^2 ) */
1597     (PID.TID 0000.0001) 90 @ 4.110969866729047E+10 /* I = 1: 90 */
1598     (PID.TID 0000.0001) ;
1599     (PID.TID 0000.0001) rA = /* rA (1,:,1,:) ( units: m^2 ) */
1600     (PID.TID 0000.0001) 4.110969866729047E+10, /* J = 1 */
1601     (PID.TID 0000.0001) 5.450087291636665E+10, /* J = 2 */
1602     (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 3 */
1603     (PID.TID 0000.0001) 8.042270623659998E+10, /* J = 4 */
1604     (PID.TID 0000.0001) 9.282707674147523E+10, /* J = 5 */
1605     (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 6 */
1606     (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 7 */
1607     (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 8 */
1608     (PID.TID 0000.0001) 1.373525277677229E+11, /* J = 9 */
1609     (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 10 */
1610     (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 11 */
1611     (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 12 */
1612     (PID.TID 0000.0001) 1.712363708253573E+11, /* J = 13 */
1613     (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 14 */
1614     (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 15 */
1615     (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 16 */
1616     (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 17 */
1617     (PID.TID 0000.0001) 1.947228162702492E+11, /* J = 18 */
1618     (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 19 */
1619     (PID.TID 0000.0001) 2 @ 1.976062798008856E+11, /* J = 20: 21 */
1620     (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 22 */
1621     (PID.TID 0000.0001) 1.947228162702492E+11, /* J = 23 */
1622     (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 24 */
1623     (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 25 */
1624     (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 26 */
1625     (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 27 */
1626     (PID.TID 0000.0001) 1.712363708253573E+11, /* J = 28 */
1627     (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 29 */
1628     (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 30 */
1629     (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 31 */
1630     (PID.TID 0000.0001) 1.373525277677229E+11, /* J = 32 */
1631     (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 33 */
1632     (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 34 */
1633     (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 35 */
1634     (PID.TID 0000.0001) 9.282707674147523E+10, /* J = 36 */
1635     (PID.TID 0000.0001) 8.042270623659998E+10, /* J = 37 */
1636     (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 38 */
1637     (PID.TID 0000.0001) 5.450087291636665E+10, /* J = 39 */
1638     (PID.TID 0000.0001) 4.110969866729047E+10 /* J = 40 */
1639     (PID.TID 0000.0001) ;
1640     (PID.TID 0000.0001) rAw = /* rAw(:,1,:,1) ( units: m^2 ) */
1641     (PID.TID 0000.0001) 90 @ 4.110969866729047E+10 /* I = 1: 90 */
1642     (PID.TID 0000.0001) ;
1643     (PID.TID 0000.0001) rAw = /* rAw(1,:,1,:) ( units: m^2 ) */
1644     (PID.TID 0000.0001) 4.110969866729047E+10, /* J = 1 */
1645     (PID.TID 0000.0001) 5.450087291636665E+10, /* J = 2 */
1646     (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 3 */
1647     (PID.TID 0000.0001) 8.042270623659998E+10, /* J = 4 */
1648     (PID.TID 0000.0001) 9.282707674147523E+10, /* J = 5 */
1649     (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 6 */
1650     (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 7 */
1651     (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 8 */
1652     (PID.TID 0000.0001) 1.373525277677229E+11, /* J = 9 */
1653     (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 10 */
1654     (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 11 */
1655     (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 12 */
1656     (PID.TID 0000.0001) 1.712363708253573E+11, /* J = 13 */
1657     (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 14 */
1658     (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 15 */
1659     (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 16 */
1660     (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 17 */
1661     (PID.TID 0000.0001) 1.947228162702492E+11, /* J = 18 */
1662     (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 19 */
1663     (PID.TID 0000.0001) 2 @ 1.976062798008856E+11, /* J = 20: 21 */
1664     (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 22 */
1665     (PID.TID 0000.0001) 1.947228162702492E+11, /* J = 23 */
1666     (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 24 */
1667     (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 25 */
1668     (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 26 */
1669     (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 27 */
1670     (PID.TID 0000.0001) 1.712363708253573E+11, /* J = 28 */
1671     (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 29 */
1672     (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 30 */
1673     (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 31 */
1674     (PID.TID 0000.0001) 1.373525277677229E+11, /* J = 32 */
1675     (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 33 */
1676     (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 34 */
1677     (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 35 */
1678     (PID.TID 0000.0001) 9.282707674147523E+10, /* J = 36 */
1679     (PID.TID 0000.0001) 8.042270623659998E+10, /* J = 37 */
1680     (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 38 */
1681     (PID.TID 0000.0001) 5.450087291636665E+10, /* J = 39 */
1682     (PID.TID 0000.0001) 4.110969866729047E+10 /* J = 40 */
1683     (PID.TID 0000.0001) ;
1684     (PID.TID 0000.0001) rAs = /* rAs(:,1,:,1) ( units: m^2 ) */
1685     (PID.TID 0000.0001) 90 @ 3.433488626798250E+10 /* I = 1: 90 */
1686     (PID.TID 0000.0001) ;
1687     (PID.TID 0000.0001) rAs = /* rAs(1,:,1,:) ( units: m^2 ) */
1688     (PID.TID 0000.0001) 3.433488626798250E+10, /* J = 1 */
1689     (PID.TID 0000.0001) 4.783442523123625E+10, /* J = 2 */
1690     (PID.TID 0000.0001) 6.110091968306414E+10, /* J = 3 */
1691     (PID.TID 0000.0001) 7.406973659603818E+10, /* J = 4 */
1692     (PID.TID 0000.0001) 8.667769319778079E+10, /* J = 5 */
1693     (PID.TID 0000.0001) 9.886336479107463E+10, /* J = 6 */
1694     (PID.TID 0000.0001) 1.105673840088173E+11, /* J = 7 */
1695     (PID.TID 0000.0001) 1.217327300458638E+11, /* J = 8 */
1696     (PID.TID 0000.0001) 1.323050064586578E+11, /* J = 9 */
1697     (PID.TID 0000.0001) 1.422327061792377E+11, /* J = 10 */
1698     (PID.TID 0000.0001) 1.514674624524945E+11, /* J = 11 */
1699     (PID.TID 0000.0001) 1.599642844741385E+11, /* J = 12 */
1700     (PID.TID 0000.0001) 1.676817765813788E+11, /* J = 13 */
1701     (PID.TID 0000.0001) 1.745823399284268E+11, /* J = 14 */
1702     (PID.TID 0000.0001) 1.806323556642996E+11, /* J = 15 */
1703     (PID.TID 0000.0001) 1.858023487204767E+11, /* J = 16 */
1704     (PID.TID 0000.0001) 1.900671314104743E+11, /* J = 17 */
1705     (PID.TID 0000.0001) 1.934059261417216E+11, /* J = 18 */
1706     (PID.TID 0000.0001) 1.958024666419019E+11, /* J = 19 */
1707     (PID.TID 0000.0001) 1.972450772065981E+11, /* J = 20 */
1708     (PID.TID 0000.0001) 1.977267295821495E+11, /* J = 21 */
1709     (PID.TID 0000.0001) 1.972450772065981E+11, /* J = 22 */
1710     (PID.TID 0000.0001) 1.958024666419019E+11, /* J = 23 */
1711     (PID.TID 0000.0001) 1.934059261417216E+11, /* J = 24 */
1712     (PID.TID 0000.0001) 1.900671314104743E+11, /* J = 25 */
1713     (PID.TID 0000.0001) 1.858023487204767E+11, /* J = 26 */
1714     (PID.TID 0000.0001) 1.806323556642996E+11, /* J = 27 */
1715     (PID.TID 0000.0001) 1.745823399284268E+11, /* J = 28 */
1716     (PID.TID 0000.0001) 1.676817765813788E+11, /* J = 29 */
1717     (PID.TID 0000.0001) 1.599642844741385E+11, /* J = 30 */
1718     (PID.TID 0000.0001) 1.514674624524945E+11, /* J = 31 */
1719     (PID.TID 0000.0001) 1.422327061792377E+11, /* J = 32 */
1720     (PID.TID 0000.0001) 1.323050064586578E+11, /* J = 33 */
1721     (PID.TID 0000.0001) 1.217327300458638E+11, /* J = 34 */
1722     (PID.TID 0000.0001) 1.105673840088173E+11, /* J = 35 */
1723     (PID.TID 0000.0001) 9.886336479107463E+10, /* J = 36 */
1724     (PID.TID 0000.0001) 8.667769319778079E+10, /* J = 37 */
1725     (PID.TID 0000.0001) 7.406973659603818E+10, /* J = 38 */
1726     (PID.TID 0000.0001) 6.110091968306414E+10, /* J = 39 */
1727     (PID.TID 0000.0001) 4.783442523123625E+10 /* J = 40 */
1728     (PID.TID 0000.0001) ;
1729     (PID.TID 0000.0001) globalArea = /* Integrated horizontal Area (m^2) */
1730     (PID.TID 0000.0001) 3.335490855875365E+14
1731     (PID.TID 0000.0001) ;
1732     (PID.TID 0000.0001) // =======================================================
1733     (PID.TID 0000.0001) // End of Model config. summary
1734     (PID.TID 0000.0001) // =======================================================
1735     (PID.TID 0000.0001)
1736 jmc 1.4 (PID.TID 0000.0001) == Packages configuration : Check & print summary ==
1737     (PID.TID 0000.0001)
1738     (PID.TID 0000.0001) GMREDI_CHECK: #define GMREDI
1739     (PID.TID 0000.0001) GM_AdvForm = /* if FALSE => use SkewFlux Form */
1740     (PID.TID 0000.0001) F
1741     (PID.TID 0000.0001) ;
1742     (PID.TID 0000.0001) GM_InMomAsStress = /* if TRUE => apply as Eddy Stress */
1743     (PID.TID 0000.0001) F
1744     (PID.TID 0000.0001) ;
1745     (PID.TID 0000.0001) GM_AdvSeparate = /* Calc Bolus & Euler Adv. separately */
1746     (PID.TID 0000.0001) F
1747     (PID.TID 0000.0001) ;
1748     (PID.TID 0000.0001) GM_ExtraDiag = /* Tensor Extra Diag (line 1&2) non 0 */
1749     (PID.TID 0000.0001) F
1750     (PID.TID 0000.0001) ;
1751     (PID.TID 0000.0001) GM_isopycK = /* Background Isopyc. Diffusivity [m^2/s] */
1752     (PID.TID 0000.0001) 1.000000000000000E+03
1753     (PID.TID 0000.0001) ;
1754     (PID.TID 0000.0001) GM_skewflx*K = /* Background GM_SkewFlx Diffusivity [m^2/s] */
1755     (PID.TID 0000.0001) 1.000000000000000E+03
1756     (PID.TID 0000.0001) ;
1757     (PID.TID 0000.0001) GM_advec*K = /* Backg. GM-Advec(=Bolus) Diffusivity [m^2/s]*/
1758     (PID.TID 0000.0001) 0.000000000000000E+00
1759     (PID.TID 0000.0001) ;
1760     (PID.TID 0000.0001) GM_Kmin_horiz = /* Minimum Horizontal Diffusivity [m^2/s] */
1761 jmc 1.6 (PID.TID 0000.0001) 5.000000000000000E+01
1762 jmc 1.4 (PID.TID 0000.0001) ;
1763     (PID.TID 0000.0001) GM_Visbeck_alpha = /* Visbeck alpha coeff. [-] */
1764     (PID.TID 0000.0001) 0.000000000000000E+00
1765     (PID.TID 0000.0001) ;
1766     (PID.TID 0000.0001) GM_Small_Number = /* epsilon used in slope calc */
1767     (PID.TID 0000.0001) 9.999999999999999E-21
1768     (PID.TID 0000.0001) ;
1769     (PID.TID 0000.0001) GM_slopeSqCutoff = /* Slope^2 cut-off value */
1770     (PID.TID 0000.0001) 1.000000000000000E+48
1771     (PID.TID 0000.0001) ;
1772     (PID.TID 0000.0001) GM_taper_scheme = /* Type of Tapering/Clipping scheme */
1773     (PID.TID 0000.0001) 'gkw91 '
1774     (PID.TID 0000.0001) ;
1775     (PID.TID 0000.0001) GM_maxSlope = /* Maximum Slope (Tapering/Clipping) */
1776     (PID.TID 0000.0001) 1.000000000000000E-02
1777     (PID.TID 0000.0001) ;
1778     (PID.TID 0000.0001) GM_facTrL2dz = /* Minimum Trans.Layer Thick. (factor of dz) */
1779     (PID.TID 0000.0001) 1.000000000000000E+00
1780     (PID.TID 0000.0001) ;
1781     (PID.TID 0000.0001) GM_facTrL2ML = /* Max.Trans.Layer Thick. (factor of MxL Depth)*/
1782     (PID.TID 0000.0001) 5.000000000000000E+00
1783     (PID.TID 0000.0001) ;
1784     (PID.TID 0000.0001) GM_maxTransLay = /* Maximum Transition Layer Thickness [m] */
1785     (PID.TID 0000.0001) 5.000000000000000E+02
1786     (PID.TID 0000.0001) ;
1787 jmc 1.7 (PID.TID 0000.0001) GM_UseBVP = /* if TRUE => use bvp a la Ferrari et al. (2010) */
1788     (PID.TID 0000.0001) F
1789     (PID.TID 0000.0001) ;
1790     (PID.TID 0000.0001) GM_BVP_ModeNumber = /* Vertical mode number for BVP wave speed */
1791     (PID.TID 0000.0001) 1
1792     (PID.TID 0000.0001) ;
1793     (PID.TID 0000.0001) GM_BVP_cMin = /* Minimum wave speed for BVP [m/s] */
1794     (PID.TID 0000.0001) 1.000000000000000E-01
1795     (PID.TID 0000.0001) ;
1796     (PID.TID 0000.0001) GM_useSubMeso = /* if TRUE => use Sub-Meso param. (B.Fox-Kemper) */
1797     (PID.TID 0000.0001) F
1798     (PID.TID 0000.0001) ;
1799     (PID.TID 0000.0001) subMeso_Ceff = /* efficiency coeff. of Mixed-Layer Eddies [-] */
1800     (PID.TID 0000.0001) 7.000000000000001E-02
1801     (PID.TID 0000.0001) ;
1802     (PID.TID 0000.0001) subMeso_invTau = /* inverse of Sub-Meso mixing time-scale [/s] */
1803     (PID.TID 0000.0001) 2.000000000000000E-06
1804     (PID.TID 0000.0001) ;
1805     (PID.TID 0000.0001) subMeso_LfMin = /* minimum length-scale "Lf" [m] */
1806     (PID.TID 0000.0001) 1.000000000000000E+03
1807     (PID.TID 0000.0001) ;
1808     (PID.TID 0000.0001) subMeso_Lmax = /* maximum grid-scale length [m] */
1809     (PID.TID 0000.0001) 1.100000000000000E+05
1810     (PID.TID 0000.0001) ;
1811 jmc 1.4 (PID.TID 0000.0001) CTRL_CHECK: ctrl package
1812     (PID.TID 0000.0001) COST_CHECK: cost package
1813 jmc 1.6 (PID.TID 0000.0001) GRDCHK_CHECK: grdchk package
1814 jmc 1.7 (PID.TID 0000.0001) GAD_CHECK: #define ALLOW_GENERIC_ADVDIFF
1815 jmc 1.4 (PID.TID 0000.0001) // =======================================================
1816 jmc 1.9 (PID.TID 0000.0001) // Check Model config. (CONFIG_CHECK):
1817 jmc 1.4 (PID.TID 0000.0001) // CONFIG_CHECK : Normal End
1818     (PID.TID 0000.0001) // =======================================================
1819     (PID.TID 0000.0001)
1820 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
1821     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
1822     (PID.TID 0000.0001)
1823 jmc 1.7 (PID.TID 0000.0001) // =======================================================
1824     (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
1825     (PID.TID 0000.0001) // =======================================================
1826     (PID.TID 0000.0001) %MON time_tsnumber = 0
1827     (PID.TID 0000.0001) %MON time_secondsf = 0.0000000000000E+00
1828     (PID.TID 0000.0001) %MON dynstat_eta_max = 0.0000000000000E+00
1829     (PID.TID 0000.0001) %MON dynstat_eta_min = 0.0000000000000E+00
1830     (PID.TID 0000.0001) %MON dynstat_eta_mean = 0.0000000000000E+00
1831     (PID.TID 0000.0001) %MON dynstat_eta_sd = 0.0000000000000E+00
1832     (PID.TID 0000.0001) %MON dynstat_eta_del2 = 0.0000000000000E+00
1833     (PID.TID 0000.0001) %MON dynstat_uvel_max = 0.0000000000000E+00
1834     (PID.TID 0000.0001) %MON dynstat_uvel_min = 0.0000000000000E+00
1835     (PID.TID 0000.0001) %MON dynstat_uvel_mean = 0.0000000000000E+00
1836     (PID.TID 0000.0001) %MON dynstat_uvel_sd = 0.0000000000000E+00
1837     (PID.TID 0000.0001) %MON dynstat_uvel_del2 = 0.0000000000000E+00
1838     (PID.TID 0000.0001) %MON dynstat_vvel_max = 0.0000000000000E+00
1839     (PID.TID 0000.0001) %MON dynstat_vvel_min = 0.0000000000000E+00
1840     (PID.TID 0000.0001) %MON dynstat_vvel_mean = 0.0000000000000E+00
1841     (PID.TID 0000.0001) %MON dynstat_vvel_sd = 0.0000000000000E+00
1842     (PID.TID 0000.0001) %MON dynstat_vvel_del2 = 0.0000000000000E+00
1843     (PID.TID 0000.0001) %MON dynstat_wvel_max = -0.0000000000000E+00
1844     (PID.TID 0000.0001) %MON dynstat_wvel_min = -0.0000000000000E+00
1845     (PID.TID 0000.0001) %MON dynstat_wvel_mean = 0.0000000000000E+00
1846     (PID.TID 0000.0001) %MON dynstat_wvel_sd = 0.0000000000000E+00
1847     (PID.TID 0000.0001) %MON dynstat_wvel_del2 = 0.0000000000000E+00
1848     (PID.TID 0000.0001) %MON dynstat_theta_max = 2.9329580307007E+01
1849     (PID.TID 0000.0001) %MON dynstat_theta_min = -1.8507585525513E+00
1850     (PID.TID 0000.0001) %MON dynstat_theta_mean = 3.4385862740244E+00
1851     (PID.TID 0000.0001) %MON dynstat_theta_sd = 4.1126810597647E+00
1852     (PID.TID 0000.0001) %MON dynstat_theta_del2 = 5.2167403232002E-03
1853     (PID.TID 0000.0001) %MON dynstat_salt_max = 3.7439460754395E+01
1854     (PID.TID 0000.0001) %MON dynstat_salt_min = 3.1028352737427E+01
1855     (PID.TID 0000.0001) %MON dynstat_salt_mean = 3.4715046474649E+01
1856     (PID.TID 0000.0001) %MON dynstat_salt_sd = 2.7634376761929E-01
1857     (PID.TID 0000.0001) %MON dynstat_salt_del2 = 8.6201567836572E-04
1858     (PID.TID 0000.0001) %MON dynstat_sst_max = 2.9329580307007E+01
1859     (PID.TID 0000.0001) %MON dynstat_sst_min = -1.6685600280762E+00
1860     (PID.TID 0000.0001) %MON dynstat_sst_mean = 1.8594263769391E+01
1861     (PID.TID 0000.0001) %MON dynstat_sst_sd = 9.1650029901491E+00
1862     (PID.TID 0000.0001) %MON dynstat_sst_del2 = 2.5285957321204E-02
1863     (PID.TID 0000.0001) %MON dynstat_sss_max = 3.7439460754395E+01
1864     (PID.TID 0000.0001) %MON dynstat_sss_min = 3.1028352737427E+01
1865     (PID.TID 0000.0001) %MON dynstat_sss_mean = 3.4857417798106E+01
1866     (PID.TID 0000.0001) %MON dynstat_sss_sd = 9.3856650790078E-01
1867     (PID.TID 0000.0001) %MON dynstat_sss_del2 = 7.1164020534050E-03
1868 jmc 1.8 (PID.TID 0000.0001) %MON forcing_qnet_max = 0.0000000000000E+00
1869     (PID.TID 0000.0001) %MON forcing_qnet_min = 0.0000000000000E+00
1870     (PID.TID 0000.0001) %MON forcing_qnet_mean = 0.0000000000000E+00
1871     (PID.TID 0000.0001) %MON forcing_qnet_sd = 0.0000000000000E+00
1872     (PID.TID 0000.0001) %MON forcing_qnet_del2 = 0.0000000000000E+00
1873 jmc 1.7 (PID.TID 0000.0001) %MON forcing_qsw_max = 0.0000000000000E+00
1874     (PID.TID 0000.0001) %MON forcing_qsw_min = 0.0000000000000E+00
1875     (PID.TID 0000.0001) %MON forcing_qsw_mean = 0.0000000000000E+00
1876     (PID.TID 0000.0001) %MON forcing_qsw_sd = 0.0000000000000E+00
1877     (PID.TID 0000.0001) %MON forcing_qsw_del2 = 0.0000000000000E+00
1878 jmc 1.8 (PID.TID 0000.0001) %MON forcing_empmr_max = 0.0000000000000E+00
1879     (PID.TID 0000.0001) %MON forcing_empmr_min = 0.0000000000000E+00
1880     (PID.TID 0000.0001) %MON forcing_empmr_mean = 0.0000000000000E+00
1881     (PID.TID 0000.0001) %MON forcing_empmr_sd = 0.0000000000000E+00
1882     (PID.TID 0000.0001) %MON forcing_empmr_del2 = 0.0000000000000E+00
1883 jmc 1.7 (PID.TID 0000.0001) %MON forcing_fu_max = 1.6004334390163E-01
1884     (PID.TID 0000.0001) %MON forcing_fu_min = -1.0329221189022E-01
1885     (PID.TID 0000.0001) %MON forcing_fu_mean = 4.1280488583042E-03
1886     (PID.TID 0000.0001) %MON forcing_fu_sd = 5.5800359441598E-02
1887     (PID.TID 0000.0001) %MON forcing_fu_del2 = 3.2370368543277E-04
1888     (PID.TID 0000.0001) %MON forcing_fv_max = 8.4633164107800E-02
1889     (PID.TID 0000.0001) %MON forcing_fv_min = -7.0490099489689E-02
1890     (PID.TID 0000.0001) %MON forcing_fv_mean = 2.8803575798275E-03
1891     (PID.TID 0000.0001) %MON forcing_fv_sd = 2.2193297771215E-02
1892     (PID.TID 0000.0001) %MON forcing_fv_del2 = 2.1170835219132E-04
1893     (PID.TID 0000.0001) %MON advcfl_uvel_max = 0.0000000000000E+00
1894     (PID.TID 0000.0001) %MON advcfl_vvel_max = 0.0000000000000E+00
1895     (PID.TID 0000.0001) %MON advcfl_wvel_max = 0.0000000000000E+00
1896     (PID.TID 0000.0001) %MON advcfl_W_hf_max = 0.0000000000000E+00
1897     (PID.TID 0000.0001) %MON pe_b_mean = 0.0000000000000E+00
1898     (PID.TID 0000.0001) %MON ke_max = 0.0000000000000E+00
1899     (PID.TID 0000.0001) %MON ke_mean = 0.0000000000000E+00
1900     (PID.TID 0000.0001) %MON ke_vol = 1.4101703507226E+18
1901     (PID.TID 0000.0001) %MON vort_r_min = 0.0000000000000E+00
1902     (PID.TID 0000.0001) %MON vort_r_max = 0.0000000000000E+00
1903     (PID.TID 0000.0001) %MON vort_a_mean = -1.9976608796703E-05
1904     (PID.TID 0000.0001) %MON vort_a_sd = 7.3718720794771E-05
1905     (PID.TID 0000.0001) %MON vort_p_mean = -2.3730263022062E-05
1906     (PID.TID 0000.0001) %MON vort_p_sd = 1.0369776513329E-04
1907     (PID.TID 0000.0001) %MON surfExpan_theta_mean = 0.0000000000000E+00
1908     (PID.TID 0000.0001) %MON surfExpan_salt_mean = 0.0000000000000E+00
1909     (PID.TID 0000.0001) // =======================================================
1910     (PID.TID 0000.0001) // End MONITOR dynamic field statistics
1911     (PID.TID 0000.0001) // =======================================================
1912 jmc 1.6 SH 5.2213892178558012 23.296536408581801
1913     NH 10.477982531836687 25.564398561832434
1914     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
1915 jmc 1.7 cg2d: Sum(rhs),rhsMax = -9.40914013369820E-15 7.56835956534411E+00
1916     SH 5.2201194813001734 23.297498110483730
1917     NH 10.466538306882399 25.565240031177893
1918 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.41900380334903E-14 8.72197575655387E+00
1919 jmc 1.7 SH 5.2213217366822278 23.295110541783401
1920     NH 10.461108876579829 25.564424809914442
1921 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.12236608895699E-13 8.58491923240257E+00
1922 jmc 1.7 SH 5.2239758527201161 23.292330991848473
1923     NH 10.454544438170274 25.562857148937155
1924 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.95340271735844E-13 8.57300301452757E+00
1925 jmc 1.7 SH 5.2277890242648937 23.291025332157801
1926     NH 10.447075395367449 25.560581831928403
1927 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.41667927389122E-13 8.54199816412785E+00
1928 jmc 1.7 SH 5.2323630972755382 23.291656354054112
1929     NH 10.438753008306207 25.558131387204696
1930 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.55049323275125E-13 8.39983291285462E+00
1931 jmc 1.7 SH 5.2374226153361292 23.293835175661229
1932     NH 10.429688799623705 25.556353938299310
1933 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.03798652900866E-13 8.14938897866034E+00
1934 jmc 1.7 SH 5.2428098519675359 23.296876966114610
1935     NH 10.420249419348920 25.555836778851869
1936 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.36120020705266E-13 7.83063879768698E+00
1937 jmc 1.7 SH 5.2484301144941696 23.300155030698079
1938     NH 10.410965029195040 25.556641789184368
1939 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.79269401942872E-13 7.50867967598465E+00
1940 jmc 1.7 SH 5.2542041359492941 23.303266303870714
1941     NH 10.402332418694868 25.558401819850985
1942 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.47919479285019E-13 7.13635421945479E+00
1943 jmc 1.7 SH 5.2600399099036395 23.306034378293475
1944     NH 10.394646890277841 25.560572845205797
1945 jmc 1.9 cg2d: Sum(rhs),rhsMax = 5.92303983637521E-13 6.72826174090860E+00
1946 jmc 1.7 SH 5.2658237553761751 23.308436294575813
1947     NH 10.387931796005716 25.562670281538317
1948 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.29854199510038E-13 6.29061611247720E+00
1949 jmc 1.7 SH 5.2714235394728348 23.310527091229915
1950     NH 10.381968237103536 25.564392021505057
1951 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.86194679066654E-13 5.82354057585758E+00
1952 jmc 1.7 SH 5.2766969400194057 23.312391248058475
1953     NH 10.376388653533922 25.565627357111339
1954 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.76625944373166E-13 5.32772671031802E+00
1955 jmc 1.7 SH 5.2815003385481987 23.314118302141765
1956     NH 10.370786591431056 25.566397306321107
1957 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.28806426530548E-13 4.81328439449781E+00
1958 jmc 1.7 SH 5.2856965016658819 23.315790314925948
1959     NH 10.364805610096198 25.566780027675406
1960 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.45976588991925E-13 4.27688077868645E+00
1961 jmc 1.7 SH 5.2891603915413157 23.317473118939176
1962     NH 10.358190016883384 25.566855652924460
1963 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.47159800906383E-13 4.08273968506819E+00
1964 jmc 1.7 SH 5.2917830540431510 23.319209962725836
1965     NH 10.350797384894520 25.566681379427052
1966 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.53821139054133E-13 4.69596425100999E+00
1967 jmc 1.7 SH 5.2934738052695707 23.321019314889838
1968 jmc 1.9 NH 10.342582254673175 25.566293875024940
1969     cg2d: Sum(rhs),rhsMax = 2.99919811208582E-13 5.24143688984563E+00
1970 jmc 1.7 SH 5.2941611378271887 23.322897857983492
1971     NH 10.333566201914152 25.565715210821566
1972 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.07545655608976E-13 5.71014649701814E+00
1973 jmc 1.1 (PID.TID 0000.0001) // Time-average data written, t-step 20
1974     (PID.TID 0000.0001)
1975 jmc 1.9 (PID.TID 0000.0001) %CHECKPOINT 20 ckptA
1976     ph-1 in thsice_cost_final
1977     ph-3 in thsice_cost_final
1978 jmc 1.6 --> objf_test(bi,bj) = 0.000000000000000D+00
1979     --> objf_tracer(bi,bj) = 0.000000000000000D+00
1980 jmc 1.9 --> objf_atl(bi,bj) = -0.371967820065626D-01
1981     local fc = -0.371967820065626D-01
1982     global fc = -0.371967820065626D-01
1983 jmc 1.6 SH 5.2213892178558012 23.296536408581801
1984     NH 10.477982531836687 25.564398561832434
1985     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
1986 jmc 1.7 cg2d: Sum(rhs),rhsMax = -9.40914013369820E-15 7.56835956534411E+00
1987     SH 5.2201194813001734 23.297498110483730
1988     NH 10.466538306882399 25.565240031177893
1989 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.41900380334903E-14 8.72197575655387E+00
1990 jmc 1.7 SH 5.2213217366822278 23.295110541783401
1991     NH 10.461108876579829 25.564424809914442
1992 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.12236608895699E-13 8.58491923240257E+00
1993 jmc 1.7 SH 5.2239758527201161 23.292330991848473
1994     NH 10.454544438170274 25.562857148937155
1995 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.95340271735844E-13 8.57300301452757E+00
1996 jmc 1.7 SH 5.2277890242648937 23.291025332157801
1997     NH 10.447075395367449 25.560581831928403
1998 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.41667927389122E-13 8.54199816412785E+00
1999 jmc 1.7 SH 5.2323630972755382 23.291656354054112
2000     NH 10.438753008306207 25.558131387204696
2001 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.55049323275125E-13 8.39983291285462E+00
2002 jmc 1.7 SH 5.2374226153361292 23.293835175661229
2003     NH 10.429688799623705 25.556353938299310
2004 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.03798652900866E-13 8.14938897866034E+00
2005 jmc 1.7 SH 5.2428098519675359 23.296876966114610
2006     NH 10.420249419348920 25.555836778851869
2007 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.36120020705266E-13 7.83063879768698E+00
2008 jmc 1.7 SH 5.2484301144941696 23.300155030698079
2009     NH 10.410965029195040 25.556641789184368
2010 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.79269401942872E-13 7.50867967598465E+00
2011 jmc 1.7 SH 5.2542041359492941 23.303266303870714
2012     NH 10.402332418694868 25.558401819850985
2013 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.47919479285019E-13 7.13635421945479E+00
2014 jmc 1.7 SH 5.2600399099036395 23.306034378293475
2015     NH 10.394646890277841 25.560572845205797
2016 jmc 1.9 cg2d: Sum(rhs),rhsMax = 5.92303983637521E-13 6.72826174090860E+00
2017 jmc 1.7 SH 5.2658237553761751 23.308436294575813
2018     NH 10.387931796005716 25.562670281538317
2019 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.29854199510038E-13 6.29061611247720E+00
2020 jmc 1.7 SH 5.2714235394728348 23.310527091229915
2021     NH 10.381968237103536 25.564392021505057
2022 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.86194679066654E-13 5.82354057585758E+00
2023 jmc 1.7 SH 5.2766969400194057 23.312391248058475
2024     NH 10.376388653533922 25.565627357111339
2025 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.76625944373166E-13 5.32772671031802E+00
2026 jmc 1.7 SH 5.2815003385481987 23.314118302141765
2027     NH 10.370786591431056 25.566397306321107
2028 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.28806426530548E-13 4.81328439449781E+00
2029 jmc 1.7 SH 5.2856965016658819 23.315790314925948
2030     NH 10.364805610096198 25.566780027675406
2031 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.45976588991925E-13 4.27688077868645E+00
2032 jmc 1.7 SH 5.2891603915413157 23.317473118939176
2033     NH 10.358190016883384 25.566855652924460
2034 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.47159800906383E-13 4.08273968506819E+00
2035 jmc 1.7 SH 5.2917830540431510 23.319209962725836
2036     NH 10.350797384894520 25.566681379427052
2037 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.53821139054133E-13 4.69596425100999E+00
2038 jmc 1.7 SH 5.2934738052695707 23.321019314889838
2039 jmc 1.9 NH 10.342582254673175 25.566293875024940
2040     cg2d: Sum(rhs),rhsMax = 2.99919811208582E-13 5.24143688984563E+00
2041 jmc 1.7 SH 5.2941611378271887 23.322897857983492
2042     NH 10.333566201914152 25.565715210821566
2043 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.07545655608976E-13 5.71014649701814E+00
2044 jmc 1.6 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2045 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.99919811208582E-13 5.24143688984563E+00
2046     cg2d: Sum(rhs),rhsMax = 3.07656677911439E-13 5.71014649701814E+00
2047 jmc 1.6 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2048     (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2049 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.77635683940025E-15 1.01287717846603E-08
2050     cg2d: Sum(rhs),rhsMax = -7.69234394064344E-14 1.09879179527660E-08
2051     cg2d: Sum(rhs),rhsMax = 2.86194679066654E-13 5.82354057585758E+00
2052     cg2d: Sum(rhs),rhsMax = 2.81170919880225E-13 5.32772671031802E+00
2053     cg2d: Sum(rhs),rhsMax = 3.40345807092746E-13 4.81328439449783E+00
2054     cg2d: Sum(rhs),rhsMax = 4.01782773717940E-13 4.27688077868644E+00
2055     cg2d: Sum(rhs),rhsMax = 3.81292220019702E-13 4.08273968506817E+00
2056     cg2d: Sum(rhs),rhsMax = 3.65686647629815E-13 4.69596425100994E+00
2057     cg2d: Sum(rhs),rhsMax = -2.55653137548606E-13 1.15636713703768E-08
2058     cg2d: Sum(rhs),rhsMax = -4.60183974260175E-13 1.26536107829254E-08
2059     cg2d: Sum(rhs),rhsMax = -6.75645303593875E-13 1.39162850775954E-08
2060     cg2d: Sum(rhs),rhsMax = -9.19134560128931E-13 1.53380614290839E-08
2061     cg2d: Sum(rhs),rhsMax = -1.17611129168971E-12 1.69103114062478E-08
2062     cg2d: Sum(rhs),rhsMax = -1.46379436349875E-12 1.86169355776582E-08
2063     cg2d: Sum(rhs),rhsMax = 3.03798652900866E-13 8.14938897866034E+00
2064     cg2d: Sum(rhs),rhsMax = 3.34871019802563E-13 7.83063879768698E+00
2065     cg2d: Sum(rhs),rhsMax = 4.69978223005540E-13 7.50867967598465E+00
2066     cg2d: Sum(rhs),rhsMax = 4.42042236148410E-13 7.13635421945480E+00
2067     cg2d: Sum(rhs),rhsMax = 6.15209272414319E-13 6.72826174090860E+00
2068     cg2d: Sum(rhs),rhsMax = 3.95759813809349E-13 6.29061611247718E+00
2069     cg2d: Sum(rhs),rhsMax = -1.63189253776785E-12 2.04399954066654E-08
2070     cg2d: Sum(rhs),rhsMax = -1.69365216295958E-12 2.23569887044410E-08
2071     cg2d: Sum(rhs),rhsMax = -1.80626173740261E-12 2.43406113494753E-08
2072     cg2d: Sum(rhs),rhsMax = -1.89170554221185E-12 2.63595082858861E-08
2073     cg2d: Sum(rhs),rhsMax = -1.86102522281573E-12 2.83800226466908E-08
2074     cg2d: Sum(rhs),rhsMax = -1.74945474773480E-12 3.03683331156591E-08
2075 jmc 1.6 (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2076 jmc 1.7 cg2d: Sum(rhs),rhsMax = -9.40914013369820E-15 7.56835956534411E+00
2077 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.41900380334903E-14 8.72197575655387E+00
2078     cg2d: Sum(rhs),rhsMax = 1.04846686888038E-13 8.58491923240257E+00
2079     cg2d: Sum(rhs),rhsMax = 1.93577792684252E-13 8.57300301452757E+00
2080     cg2d: Sum(rhs),rhsMax = 1.48010076417293E-13 8.54199816412785E+00
2081     cg2d: Sum(rhs),rhsMax = 3.96654931122953E-13 8.39983291285461E+00
2082     cg2d: Sum(rhs),rhsMax = -1.60096241819119E-12 3.22921328128692E-08
2083     cg2d: Sum(rhs),rhsMax = -1.44032962640095E-12 3.41215550713648E-08
2084     cg2d: Sum(rhs),rhsMax = -1.32341707037575E-12 3.58294707242427E-08
2085     cg2d: Sum(rhs),rhsMax = -1.25662501238022E-12 3.73920335700041E-08
2086     cg2d: Sum(rhs),rhsMax = -1.21928248747460E-12 3.87900997765171E-08
2087     cg2d: Sum(rhs),rhsMax = -1.22127048057807E-12 4.00072765226610E-08
2088 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2089     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2090     (PID.TID 0000.0001)
2091     ph-pack: packing ecco_cost
2092     ph-pack: packing ecco_ctrl
2093 jmc 1.9 (PID.TID 0000.0001) // =======================================================
2094     (PID.TID 0000.0001) // Gradient-check starts (grdchk_main)
2095     (PID.TID 0000.0001) // =======================================================
2096     (PID.TID 0000.0001) grdchk reference fc: fcref = -3.71967820065626E-02
2097 jmc 1.1 grad-res -------------------------------
2098 jmc 1.9 grad-res proc # i j k bi bj iobc fc ref fc + eps fc - eps
2099     grad-res proc # i j k bi bj iobc adj grad fd grad 1 - fd/adj
2100     (PID.TID 0000.0001) ====== Starts gradient-check number 1 (=ichknum) =======
2101 jmc 1.7 ph-test icomp, ncvarcomp, ichknum 1 29138 1
2102     ph-grd _loc: bi, bj, icomptest, ichknum 1 1 0 1
2103 jmc 1.6 ph-grd -->hit<-- 44 2 1 1
2104 jmc 1.9 (PID.TID 0000.0001) grdchk pos: i,j,k= 44 2 1 ; bi,bj= 1 1 ; iobc= 1 ; rec= 1
2105 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2106     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2107     (PID.TID 0000.0001)
2108     (PID.TID 0000.0001) // =======================================================
2109     (PID.TID 0000.0001) // Model current state
2110     (PID.TID 0000.0001) // =======================================================
2111     (PID.TID 0000.0001)
2112 jmc 1.6 SH 5.2214032038697837 23.296536408581801
2113     NH 10.477982531836687 25.564398561832434
2114     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2115 jmc 1.7 cg2d: Sum(rhs),rhsMax = -9.40914013369820E-15 7.56835956534411E+00
2116     SH 5.2201334211142694 23.297498128510995
2117     NH 10.466538306882399 25.565240031177893
2118 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.64451785522601E-15 8.72197575654855E+00
2119 jmc 1.7 SH 5.2213357119372290 23.295110577687069
2120     NH 10.461108876580283 25.564424809919160
2121 jmc 1.9 cg2d: Sum(rhs),rhsMax = 9.06288932789323E-14 8.58491923233593E+00
2122 jmc 1.7 SH 5.2239898826075546 23.292331045693061
2123     NH 10.454544438174480 25.562857148972114
2124 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.16996559609939E-13 8.57300301424113E+00
2125 jmc 1.7 SH 5.2278030705794549 23.291025403982246
2126     NH 10.447075395385923 25.560581832030721
2127 jmc 1.9 cg2d: Sum(rhs),rhsMax = 9.99027249815043E-14 8.54199816328915E+00
2128 jmc 1.7 SH 5.2323771053971768 23.291656443776731
2129     NH 10.438753008359129 25.558131387418598
2130 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.45390382960886E-13 8.39983291096013E+00
2131 jmc 1.7 SH 5.2374365523609754 23.293835283113992
2132     NH 10.429688799736539 25.556353938648101
2133 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.25781068788444E-13 8.14938897519342E+00
2134 jmc 1.7 SH 5.2428237119067944 23.296877091127328
2135     NH 10.420249419540937 25.555836779317190
2136 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.74804354219549E-13 7.83063879238954E+00
2137 jmc 1.7 SH 5.2484439044173996 23.300155173134044
2138     NH 10.410965029465672 25.556641789710397
2139 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.34437208429728E-13 7.50867966909402E+00
2140 jmc 1.7 SH 5.2542178634641346 23.303266463613056
2141     NH 10.402332419017455 25.558401820371557
2142 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.83987849797529E-13 7.13635421172750E+00
2143 jmc 1.7 SH 5.2600535783429905 23.306034555220720
2144     NH 10.394646890607135 25.560572845673278
2145 jmc 1.9 cg2d: Sum(rhs),rhsMax = 6.25076379545675E-13 6.72826173338455E+00
2146 jmc 1.7 SH 5.2658373644449625 23.308436488551081
2147     NH 10.387931796293969 25.562670281937017
2148 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.21132009872827E-13 6.29061610611811E+00
2149 jmc 1.7 SH 5.2714370875355030 23.310527302104560
2150     NH 10.381968237317153 25.564392021846800
2151 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.49009146635615E-13 5.82354057125141E+00
2152 jmc 1.7 SH 5.2767104258373934 23.312391475681295
2153     NH 10.376388653662389 25.565627357421590
2154 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.11657080750882E-13 5.32772670756981E+00
2155 jmc 1.7 SH 5.2815137620302721 23.314118546367354
2156 jmc 1.9 NH 10.370786591485682 25.566397306626698
2157     cg2d: Sum(rhs),rhsMax = 3.16899284591443E-13 4.81328439334297E+00
2158 jmc 1.7 SH 5.2857098639553328 23.315790575619417
2159     NH 10.364805610102223 25.566780027998220
2160 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.66043870162025E-13 4.27688077870375E+00
2161 jmc 1.7 SH 5.2891736947930212 23.317473395977515
2162     NH 10.358190016870481 25.566855653280697
2163 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.64673569119844E-13 4.08273968597044E+00
2164 jmc 1.7 SH 5.2917963010744984 23.319210255997064
2165     NH 10.350797384889006 25.566681379829447
2166 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.58567342484406E-13 4.69596425233109E+00
2167 jmc 1.7 SH 5.2934869991663049 23.321019624290443
2168 jmc 1.9 NH 10.342582254694207 25.566293875484519
2169     cg2d: Sum(rhs),rhsMax = 3.22568360910935E-13 5.24143689138270E+00
2170 jmc 1.7 SH 5.2941742815707471 23.322898183415759
2171     NH 10.333566201972985 25.565715211348220
2172 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.92238455656957E-13 5.71014649850039E+00
2173 jmc 1.6 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2174 jmc 1.9 ph-1 in thsice_cost_final
2175     ph-3 in thsice_cost_final
2176 jmc 1.6 --> objf_test(bi,bj) = 0.000000000000000D+00
2177     --> objf_tracer(bi,bj) = 0.000000000000000D+00
2178 jmc 1.9 --> objf_atl(bi,bj) = -0.371967821606290D-01
2179     local fc = -0.371967821606290D-01
2180     global fc = -0.371967821606290D-01
2181     (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = -3.71967821606290E-02
2182 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2183     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2184     (PID.TID 0000.0001)
2185     (PID.TID 0000.0001) // =======================================================
2186     (PID.TID 0000.0001) // Model current state
2187     (PID.TID 0000.0001) // =======================================================
2188     (PID.TID 0000.0001)
2189 jmc 1.6 SH 5.2213752318418187 23.296536408581801
2190     NH 10.477982531836687 25.564398561832434
2191     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2192 jmc 1.7 cg2d: Sum(rhs),rhsMax = -9.40914013369820E-15 7.56835956534411E+00
2193     SH 5.2201055417050384 23.297498092456351
2194     NH 10.466538306882399 25.565240031177893
2195 jmc 1.9 cg2d: Sum(rhs),rhsMax = -9.14546216534973E-15 8.72197575655921E+00
2196 jmc 1.7 SH 5.2213077618590091 23.295110505880132
2197     NH 10.461108876579431 25.564424809909838
2198 jmc 1.9 cg2d: Sum(rhs),rhsMax = 9.89555659636210E-14 8.58491923246909E+00
2199 jmc 1.7 SH 5.2239618235048511 23.292330938005250
2200     NH 10.454544438166067 25.562857148902197
2201 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.73208669629332E-13 8.57300301481357E+00
2202 jmc 1.7 SH 5.2277749788662504 23.291025260335857
2203     NH 10.447075395348975 25.560581831826198
2204 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.40477907084602E-13 8.54199816496515E+00
2205 jmc 1.7 SH 5.2323490903098673 23.291656264335529
2206     NH 10.438753008253457 25.558131386991079
2207 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.51343953930439E-13 8.39983291474608E+00
2208 jmc 1.7 SH 5.2374086797000814 23.293835068213980
2209     NH 10.429688799511098 25.556353937950860
2210 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.16656423304806E-13 8.14938898212175E+00
2211 jmc 1.7 SH 5.2427959936439947 23.296876841109395
2212     NH 10.420249419157130 25.555836778387174
2213 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.04655475794158E-13 7.83063880297591E+00
2214 jmc 1.7 SH 5.2484163264093127 23.300154888271607
2215     NH 10.410965028924863 25.556641788659135
2216 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.41667535877599E-13 7.50867968286433E+00
2217 jmc 1.7 SH 5.2541904104915034 23.303266144140309
2218     NH 10.402332418372737 25.558401819331152
2219 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.82246187427648E-13 7.13635422716987E+00
2220 jmc 1.7 SH 5.2600262437351830 23.306034201381010
2221     NH 10.394646889949001 25.560572844738886
2222 jmc 1.9 cg2d: Sum(rhs),rhsMax = 5.82763004519649E-13 6.72826174842132E+00
2223 jmc 1.7 SH 5.2658101487857607 23.308436100618223
2224     NH 10.387931795717861 25.562670281140015
2225 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.07635861229727E-13 6.29061611882678E+00
2226 jmc 1.7 SH 5.2714099940886285 23.310526880376244
2227     NH 10.381968236890202 25.564392021163712
2228 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.09655079355764E-13 5.82354058045789E+00
2229 jmc 1.7 SH 5.2766834570720675 23.312391020460154
2230 jmc 1.9 NH 10.376388653405570 25.565627356801258
2231     cg2d: Sum(rhs),rhsMax = 3.50788842418126E-13 5.32772671306349E+00
2232 jmc 1.7 SH 5.2814869181208906 23.314118057944484
2233     NH 10.370786591376486 25.566397306015574
2234 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.24368873375158E-13 4.81328439565244E+00
2235 jmc 1.7 SH 5.2856831426076951 23.315790054264824
2236     NH 10.364805610090059 25.566780027352706
2237 jmc 1.9 cg2d: Sum(rhs),rhsMax = 5.19966014689288E-13 4.27688077867098E+00
2238 jmc 1.7 SH 5.2891470916899834 23.317472841937274
2239     NH 10.358190016896174 25.566855652568336
2240 jmc 1.9 cg2d: Sum(rhs),rhsMax = 5.28826982204578E-13 4.08273968416910E+00
2241 jmc 1.7 SH 5.2917698105745785 23.319209669495365
2242     NH 10.350797384899920 25.566681379024772
2243 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.09353106967103E-13 4.69596424969276E+00
2244 jmc 1.7 SH 5.2934606150915329 23.321019005534538
2245     NH 10.342582254652029 25.566293874565417
2246 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.46753745109163E-13 5.24143688831332E+00
2247 jmc 1.7 SH 5.2941479979523365 23.322897532601303
2248     NH 10.333566201855206 25.565715210295195
2249 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.75068032187897E-13 5.71014649554047E+00
2250 jmc 1.6 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2251 jmc 1.9 ph-1 in thsice_cost_final
2252     ph-3 in thsice_cost_final
2253 jmc 1.6 --> objf_test(bi,bj) = 0.000000000000000D+00
2254     --> objf_tracer(bi,bj) = 0.000000000000000D+00
2255 jmc 1.9 --> objf_atl(bi,bj) = -0.371967818527124D-01
2256     local fc = -0.371967818527124D-01
2257     global fc = -0.371967818527124D-01
2258     (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = -3.71967818527124E-02
2259 jmc 1.1 grad-res -------------------------------
2260 jmc 1.9 grad-res 0 1 44 2 1 1 1 1 -3.71967820066E-02 -3.71967821606E-02 -3.71967818527E-02
2261     grad-res 0 1 1 1 0 1 1 1 -1.53957497164E-08 -1.53958287091E-08 -5.13081585751E-06
2262     (PID.TID 0000.0001) ADM ref_cost_function = -3.71967820065626E-02
2263     (PID.TID 0000.0001) ADM adjoint_gradient = -1.53957497163737E-08
2264     (PID.TID 0000.0001) ADM finite-diff_grad = -1.53958287091305E-08
2265     (PID.TID 0000.0001) ====== End of gradient-check number 1 (ierr= 0) =======
2266     (PID.TID 0000.0001) ====== Starts gradient-check number 2 (=ichknum) =======
2267 jmc 1.7 ph-test icomp, ncvarcomp, ichknum 2 29138 2
2268     ph-grd _loc: bi, bj, icomptest, ichknum 1 1 1 2
2269 jmc 1.6 ph-grd -->hit<-- 45 2 1 1
2270 jmc 1.9 (PID.TID 0000.0001) grdchk pos: i,j,k= 45 2 1 ; bi,bj= 1 1 ; iobc= 1 ; rec= 1
2271 jmc 1.1 (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.6 SH 5.2214032038697837 23.296536408581801
2279     NH 10.477982531836687 25.564398561832434
2280     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2281 jmc 1.7 cg2d: Sum(rhs),rhsMax = -9.40914013369820E-15 7.56835956534411E+00
2282     SH 5.2201334270583857 23.297498128510995
2283     NH 10.466538306882399 25.565240031177893
2284 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.16528184271425E-14 8.72197575654960E+00
2285 jmc 1.7 SH 5.2213357268121854 23.295110577694288
2286     NH 10.461108876580226 25.564424809918648
2287 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.24230487008603E-13 8.58491923234548E+00
2288 jmc 1.7 SH 5.2239899090945414 23.292331045714661
2289     NH 10.454544438173912 25.562857148969101
2290 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.67837965747708E-13 8.57300301428348E+00
2291 jmc 1.7 SH 5.2278031079794687 23.291025404036702
2292     NH 10.447075395383706 25.560581832021569
2293 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.43406120312051E-13 8.54199816341439E+00
2294 jmc 1.7 SH 5.2323771548610694 23.291656443887064
2295     NH 10.438753008352990 25.558131387399897
2296 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.97605559587788E-13 8.39983291123391E+00
2297 jmc 1.7 SH 5.2374366157289387 23.293835283304702
2298     NH 10.429688799724431 25.556353938620191
2299 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.72656897060131E-13 8.14938897566263E+00
2300 jmc 1.7 SH 5.2428237903615695 23.296877091418480
2301     NH 10.420249419522293 25.555836779285016
2302 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.74061892571831E-13 7.83063879303954E+00
2303 jmc 1.7 SH 5.2484439982684989 23.300155173541270
2304     NH 10.410965029442991 25.556641789681009
2305 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.81559236931162E-13 7.50867966982532E+00
2306 jmc 1.7 SH 5.2542179726058862 23.303266464150965
2307     NH 10.402332418996139 25.558401820350298
2308 jmc 1.9 cg2d: Sum(rhs),rhsMax = 5.26384491550402E-13 7.13635421237813E+00
2309 jmc 1.7 SH 5.2600537024886762 23.306034555906137
2310     NH 10.394646890593151 25.560572845661966
2311 jmc 1.9 cg2d: Sum(rhs),rhsMax = 5.52648204976691E-13 6.72826173378720E+00
2312 jmc 1.7 SH 5.2658375031680862 23.308436489402823
2313     NH 10.387931796291639 25.562670281933833
2314 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.77260722661532E-13 6.29061610616139E+00
2315 jmc 1.7 SH 5.2714372402481331 23.310527303142294
2316     NH 10.381968237327385 25.564392021848789
2317 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.18946258293096E-13 5.82354057091309E+00
2318 jmc 1.7 SH 5.2767105917858430 23.312391476923892
2319     NH 10.376388653682454 25.565627357426592
2320 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.64649413495022E-13 5.32772670691110E+00
2321 jmc 1.7 SH 5.2815139403248281 23.314118547832152
2322     NH 10.370786591511035 25.566397306633576
2323 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.79966891284056E-13 4.81328439247409E+00
2324 jmc 1.7 SH 5.2857100536156167 23.315790577322105
2325     NH 10.364805610127746 25.566780028006690
2326 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.49064396779164E-13 4.27688077774310E+00
2327 jmc 1.7 SH 5.2891738947913041 23.317473397931735
2328     NH 10.358190016891854 25.566855653290588
2329 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.27838320327112E-13 4.08273968521706E+00
2330 jmc 1.7 SH 5.2917965103699203 23.319210258214753
2331     NH 10.350797384903501 25.566681379840020
2332 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.77715350716079E-13 4.69596425165899E+00
2333 jmc 1.7 SH 5.2934872167308527 23.321019626781776
2334     NH 10.342582254700403 25.566293875494353
2335 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.89112353449406E-13 5.24143689084384E+00
2336 jmc 1.7 SH 5.2941745064091492 23.322898186189491
2337     NH 10.333566201970655 25.565715211355837
2338 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.17947057570933E-13 5.71014649814151E+00
2339 jmc 1.6 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2340 jmc 1.9 ph-1 in thsice_cost_final
2341     ph-3 in thsice_cost_final
2342 jmc 1.6 --> objf_test(bi,bj) = 0.000000000000000D+00
2343     --> objf_tracer(bi,bj) = 0.000000000000000D+00
2344 jmc 1.9 --> objf_atl(bi,bj) = -0.371967821434007D-01
2345     local fc = -0.371967821434007D-01
2346     global fc = -0.371967821434007D-01
2347     (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = -3.71967821434007E-02
2348 jmc 1.1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
2349     (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
2350     (PID.TID 0000.0001)
2351     (PID.TID 0000.0001) // =======================================================
2352     (PID.TID 0000.0001) // Model current state
2353     (PID.TID 0000.0001) // =======================================================
2354     (PID.TID 0000.0001)
2355 jmc 1.6 SH 5.2213752318418187 23.296536408581801
2356     NH 10.477982531836687 25.564398561832434
2357     (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
2358 jmc 1.7 cg2d: Sum(rhs),rhsMax = -9.40914013369820E-15 7.56835956534411E+00
2359     SH 5.2201055357530208 23.297498092456351
2360     NH 10.466538306882399 25.565240031177893
2361 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.11001224773122E-14 8.72197575655813E+00
2362 jmc 1.7 SH 5.2213077469803579 23.295110505872913
2363     NH 10.461108876579488 25.564424809910292
2364 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.28824034772990E-13 8.58491923245957E+00
2365 jmc 1.7 SH 5.2239617970200811 23.292330937983650
2366 jmc 1.9 NH 10.454544438166636 25.562857148905209
2367     cg2d: Sum(rhs),rhsMax = 1.69201458399826E-13 8.57300301477123E+00
2368 jmc 1.7 SH 5.2277749414750474 23.291025260281288
2369     NH 10.447075395351249 25.560581831835407
2370 jmc 1.9 cg2d: Sum(rhs),rhsMax = 1.47128836891497E-13 8.54199816484020E+00
2371 jmc 1.7 SH 5.2323490408621751 23.291656264225253
2372     NH 10.438753008259539 25.558131387009780
2373 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.53439499889419E-13 8.39983291447280E+00
2374 jmc 1.7 SH 5.2374086163554239 23.293835068023270
2375     NH 10.429688799523205 25.556353937978770
2376 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.79144762860284E-13 8.14938898165323E+00
2377 jmc 1.7 SH 5.2427959152182666 23.296876840818243
2378     NH 10.420249419175775 25.555836778419348
2379 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.89320520266523E-13 7.83063880232717E+00
2380 jmc 1.7 SH 5.2484162325908414 23.300154887864437
2381     NH 10.410965028947487 25.556641788688466
2382 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.54997151067005E-13 7.50867968213432E+00
2383 jmc 1.7 SH 5.2541903013832325 23.303266143602514
2384     NH 10.402332418393996 25.558401819352241
2385 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.30641633464290E-13 7.13635422652038E+00
2386 jmc 1.7 SH 5.2600261196207043 23.306034200695763
2387     NH 10.394646889962928 25.560572844750197
2388 jmc 1.9 cg2d: Sum(rhs),rhsMax = 5.91811322170344E-13 6.72826174801940E+00
2389 jmc 1.7 SH 5.2658100100883871 23.308436099766652
2390     NH 10.387931795720192 25.562670281143255
2391 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.36751460050522E-13 6.29061611878358E+00
2392 jmc 1.7 SH 5.2714098413931652 23.310526879338795
2393     NH 10.381968236880027 25.564392021161666
2394 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.51763887515466E-13 5.82354058079576E+00
2395 jmc 1.7 SH 5.2766832911291885 23.312391019217841
2396     NH 10.376388653385504 25.565627356796313
2397 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.47232789796215E-13 5.32772671372097E+00
2398 jmc 1.7 SH 5.2814867398175238 23.314118056479856
2399     NH 10.370786591351134 25.566397306008753
2400 jmc 1.9 cg2d: Sum(rhs),rhsMax = 3.26780269510607E-13 4.81328439651968E+00
2401 jmc 1.7 SH 5.2856829529216611 23.315790052562363
2402     NH 10.364805610064536 25.566780027344180
2403 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.56717996755174E-13 4.27688077962971E+00
2404     SH 5.2891468916468511 23.317472839983111
2405 jmc 1.7 NH 10.358190016874801 25.566855652558388
2406 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.40279757096818E-13 4.08273968492202E+00
2407 jmc 1.7 SH 5.2917696012132183 23.319209667277732
2408 jmc 1.9 NH 10.350797384885482 25.566681379014256
2409     cg2d: Sum(rhs),rhsMax = 3.33642835581571E-13 4.69596425036451E+00
2410 jmc 1.7 SH 5.2934603974383663 23.321019003043148
2411     NH 10.342582254645833 25.566293874555640
2412 jmc 1.9 cg2d: Sum(rhs),rhsMax = 4.16333634234434E-13 5.24143688885200E+00
2413 jmc 1.7 SH 5.2941477730013276 23.322897529827458
2414     NH 10.333566201857593 25.565715210287578
2415 jmc 1.9 cg2d: Sum(rhs),rhsMax = 2.30045149596236E-13 5.71014649589941E+00
2416 jmc 1.6 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
2417 jmc 1.9 ph-1 in thsice_cost_final
2418     ph-3 in thsice_cost_final
2419 jmc 1.6 --> objf_test(bi,bj) = 0.000000000000000D+00
2420     --> objf_tracer(bi,bj) = 0.000000000000000D+00
2421 jmc 1.9 --> objf_atl(bi,bj) = -0.371967818699053D-01
2422     local fc = -0.371967818699053D-01
2423     global fc = -0.371967818699053D-01
2424     (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = -3.71967818699053E-02
2425 jmc 1.1 grad-res -------------------------------
2426 jmc 1.9 grad-res 0 2 45 2 1 1 1 1 -3.71967820066E-02 -3.71967821434E-02 -3.71967818699E-02
2427     grad-res 0 2 2 2 0 1 1 1 -1.36746649794E-08 -1.36747717316E-08 -7.80657355226E-06
2428     (PID.TID 0000.0001) ADM ref_cost_function = -3.71967820065626E-02
2429     (PID.TID 0000.0001) ADM adjoint_gradient = -1.36746649793641E-08
2430     (PID.TID 0000.0001) ADM finite-diff_grad = -1.36747717316421E-08
2431     (PID.TID 0000.0001) ====== End of gradient-check number 2 (ierr= 0) =======
2432 jmc 1.1 (PID.TID 0000.0001)
2433     (PID.TID 0000.0001) // =======================================================
2434     (PID.TID 0000.0001) // Gradient check results >>> START <<<
2435     (PID.TID 0000.0001) // =======================================================
2436     (PID.TID 0000.0001)
2437 jmc 1.6 (PID.TID 0000.0001) EPS = 1.000000E-02
2438 jmc 1.1 (PID.TID 0000.0001)
2439 jmc 1.7 (PID.TID 0000.0001) grdchk output h.p: Id Itile Jtile LAYER bi bj X(Id) X(Id)+/-EPS
2440 jmc 1.6 (PID.TID 0000.0001) grdchk output h.c: Id FC FC1 FC2
2441     (PID.TID 0000.0001) grdchk output h.g: Id FC1-FC2/(2*EPS) ADJ GRAD(FC) 1-FDGRD/ADGRD
2442     (PID.TID 0000.0001)
2443 jmc 1.7 (PID.TID 0000.0001) grdchk output (p): 1 44 2 1 1 1 0.000000000E+00 -1.000000000E-02
2444 jmc 1.9 (PID.TID 0000.0001) grdchk output (c): 1 -3.7196782006563E-02 -3.7196782160629E-02 -3.7196781852712E-02
2445     (PID.TID 0000.0001) grdchk output (g): 1 -1.5395828709130E-08 -1.5395749716374E-08 -5.1308158575125E-06
2446 jmc 1.7 (PID.TID 0000.0001)
2447     (PID.TID 0000.0001) grdchk output (p): 2 45 2 1 1 1 0.000000000E+00 -1.000000000E-02
2448 jmc 1.9 (PID.TID 0000.0001) grdchk output (c): 2 -3.7196782006563E-02 -3.7196782143401E-02 -3.7196781869905E-02
2449     (PID.TID 0000.0001) grdchk output (g): 2 -1.3674771731642E-08 -1.3674664979364E-08 -7.8065735522603E-06
2450     (PID.TID 0000.0001)
2451     (PID.TID 0000.0001) grdchk summary : RMS of 2 ratios = 6.6055984585256E-06
2452 jmc 1.1 (PID.TID 0000.0001)
2453     (PID.TID 0000.0001) // =======================================================
2454     (PID.TID 0000.0001) // Gradient check results >>> END <<<
2455     (PID.TID 0000.0001) // =======================================================
2456     (PID.TID 0000.0001)
2457     (PID.TID 0000.0001) Seconds in section "ALL [THE_MODEL_MAIN]":
2458 jmc 1.9 (PID.TID 0000.0001) User time: 159.14000000000001
2459     (PID.TID 0000.0001) System time: 0.48999999999999999
2460     (PID.TID 0000.0001) Wall clock time: 161.51889920234680
2461 jmc 1.6 (PID.TID 0000.0001) No. starts: 1
2462     (PID.TID 0000.0001) No. stops: 1
2463 jmc 1.1 (PID.TID 0000.0001) Seconds in section "INITIALISE_FIXED [THE_MODEL_MAIN]":
2464 jmc 1.9 (PID.TID 0000.0001) User time: 0.22000000000000000
2465 jmc 1.8 (PID.TID 0000.0001) System time: 2.00000000000000004E-002
2466 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 0.56264901161193848
2467 jmc 1.6 (PID.TID 0000.0001) No. starts: 1
2468     (PID.TID 0000.0001) No. stops: 1
2469 jmc 1.1 (PID.TID 0000.0001) Seconds in section "ADTHE_MAIN_LOOP [ADJOINT RUN]":
2470 jmc 1.9 (PID.TID 0000.0001) User time: 94.980000000000004
2471     (PID.TID 0000.0001) System time: 0.45000000000000001
2472     (PID.TID 0000.0001) Wall clock time: 96.793227910995483
2473 jmc 1.6 (PID.TID 0000.0001) No. starts: 1
2474     (PID.TID 0000.0001) No. stops: 1
2475 jmc 1.9 (PID.TID 0000.0001) Seconds in section "FORWARD_STEP [MAIN_DO_LOOP]":
2476     (PID.TID 0000.0001) User time: 95.480000000000103
2477     (PID.TID 0000.0001) System time: 7.00000000000000622E-002
2478     (PID.TID 0000.0001) Wall clock time: 95.983871936798096
2479     (PID.TID 0000.0001) No. starts: 120
2480     (PID.TID 0000.0001) No. stops: 120
2481 jmc 1.1 (PID.TID 0000.0001) Seconds in section "LOAD_FIELDS_DRIVER [FORWARD_STEP]":
2482 jmc 1.8 (PID.TID 0000.0001) User time: 9.99999999999090505E-003
2483 jmc 1.6 (PID.TID 0000.0001) System time: 0.0000000000000000
2484 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 3.36146354675292969E-003
2485     (PID.TID 0000.0001) No. starts: 120
2486     (PID.TID 0000.0001) No. stops: 120
2487 jmc 1.8 (PID.TID 0000.0001) Seconds in section "EXTERNAL_FLDS_LOAD [LOAD_FLDS_DRIVER]":
2488 jmc 1.9 (PID.TID 0000.0001) User time: 9.99999999999090505E-003
2489 jmc 1.8 (PID.TID 0000.0001) System time: 0.0000000000000000
2490 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 1.33252143859863281E-003
2491     (PID.TID 0000.0001) No. starts: 140
2492     (PID.TID 0000.0001) No. stops: 140
2493 jmc 1.1 (PID.TID 0000.0001) Seconds in section "EBM [FORWARD_STEP]":
2494 jmc 1.9 (PID.TID 0000.0001) User time: 0.28999999999999204
2495 jmc 1.6 (PID.TID 0000.0001) System time: 0.0000000000000000
2496 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 0.28806734085083008
2497     (PID.TID 0000.0001) No. starts: 120
2498     (PID.TID 0000.0001) No. stops: 120
2499 jmc 1.1 (PID.TID 0000.0001) Seconds in section "DO_ATMOSPHERIC_PHYS [FORWARD_STEP]":
2500 jmc 1.9 (PID.TID 0000.0001) User time: 0.0000000000000000
2501 jmc 1.6 (PID.TID 0000.0001) System time: 0.0000000000000000
2502 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 1.10316276550292969E-003
2503     (PID.TID 0000.0001) No. starts: 120
2504     (PID.TID 0000.0001) No. stops: 120
2505 jmc 1.1 (PID.TID 0000.0001) Seconds in section "DO_OCEANIC_PHYS [FORWARD_STEP]":
2506 jmc 1.9 (PID.TID 0000.0001) User time: 12.789999999999964
2507     (PID.TID 0000.0001) System time: 0.0000000000000000
2508     (PID.TID 0000.0001) Wall clock time: 12.798147439956665
2509     (PID.TID 0000.0001) No. starts: 120
2510     (PID.TID 0000.0001) No. stops: 120
2511 jmc 1.1 (PID.TID 0000.0001) Seconds in section "THERMODYNAMICS [FORWARD_STEP]":
2512 jmc 1.9 (PID.TID 0000.0001) User time: 24.959999999999980
2513     (PID.TID 0000.0001) System time: 1.00000000000000089E-002
2514     (PID.TID 0000.0001) Wall clock time: 25.012357234954834
2515     (PID.TID 0000.0001) No. starts: 120
2516     (PID.TID 0000.0001) No. stops: 120
2517 jmc 1.1 (PID.TID 0000.0001) Seconds in section "DYNAMICS [FORWARD_STEP]":
2518 jmc 1.9 (PID.TID 0000.0001) User time: 35.149999999999977
2519     (PID.TID 0000.0001) System time: 0.0000000000000000
2520     (PID.TID 0000.0001) Wall clock time: 35.247949361801147
2521     (PID.TID 0000.0001) No. starts: 120
2522     (PID.TID 0000.0001) No. stops: 120
2523 jmc 1.1 (PID.TID 0000.0001) Seconds in section "SOLVE_FOR_PRESSURE [FORWARD_STEP]":
2524 jmc 1.9 (PID.TID 0000.0001) User time: 16.409999999999997
2525     (PID.TID 0000.0001) System time: 0.0000000000000000
2526     (PID.TID 0000.0001) Wall clock time: 16.457117080688477
2527     (PID.TID 0000.0001) No. starts: 120
2528     (PID.TID 0000.0001) No. stops: 120
2529 jmc 1.1 (PID.TID 0000.0001) Seconds in section "MOM_CORRECTION_STEP [FORWARD_STEP]":
2530 jmc 1.9 (PID.TID 0000.0001) User time: 0.71000000000000796
2531 jmc 1.8 (PID.TID 0000.0001) System time: 0.0000000000000000
2532 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 0.69007325172424316
2533     (PID.TID 0000.0001) No. starts: 120
2534     (PID.TID 0000.0001) No. stops: 120
2535     (PID.TID 0000.0001) Seconds in section "INTEGR_CONTINUITY [FORWARD_STEP]":
2536     (PID.TID 0000.0001) User time: 1.1999999999999886
2537     (PID.TID 0000.0001) System time: 0.0000000000000000
2538     (PID.TID 0000.0001) Wall clock time: 1.2197821140289307
2539     (PID.TID 0000.0001) No. starts: 120
2540     (PID.TID 0000.0001) No. stops: 120
2541 jmc 1.1 (PID.TID 0000.0001) Seconds in section "TRC_CORRECTION_STEP [FORWARD_STEP]":
2542 jmc 1.9 (PID.TID 0000.0001) User time: 0.31999999999990791
2543 jmc 1.6 (PID.TID 0000.0001) System time: 0.0000000000000000
2544 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 0.35293221473693848
2545     (PID.TID 0000.0001) No. starts: 120
2546     (PID.TID 0000.0001) No. stops: 120
2547 jmc 1.1 (PID.TID 0000.0001) Seconds in section "BLOCKING_EXCHANGES [FORWARD_STEP]":
2548 jmc 1.9 (PID.TID 0000.0001) User time: 0.49000000000009436
2549 jmc 1.8 (PID.TID 0000.0001) System time: 0.0000000000000000
2550 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 0.47250080108642578
2551     (PID.TID 0000.0001) No. starts: 120
2552     (PID.TID 0000.0001) No. stops: 120
2553 jmc 1.1 (PID.TID 0000.0001) Seconds in section "DO_STATEVARS_TAVE [FORWARD_STEP]":
2554 jmc 1.9 (PID.TID 0000.0001) User time: 0.95999999999997954
2555 jmc 1.8 (PID.TID 0000.0001) System time: 0.0000000000000000
2556 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 0.98884296417236328
2557     (PID.TID 0000.0001) No. starts: 120
2558     (PID.TID 0000.0001) No. stops: 120
2559 jmc 1.7 (PID.TID 0000.0001) Seconds in section "MONITOR [FORWARD_STEP]":
2560 jmc 1.9 (PID.TID 0000.0001) User time: 0.0000000000000000
2561 jmc 1.6 (PID.TID 0000.0001) System time: 0.0000000000000000
2562 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 2.21967697143554688E-004
2563     (PID.TID 0000.0001) No. starts: 20
2564     (PID.TID 0000.0001) No. stops: 20
2565 jmc 1.1 (PID.TID 0000.0001) Seconds in section "COST_TILE [FORWARD_STEP]":
2566 jmc 1.9 (PID.TID 0000.0001) User time: 1.8700000000000045
2567 jmc 1.8 (PID.TID 0000.0001) System time: 0.0000000000000000
2568 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 1.8560354709625244
2569     (PID.TID 0000.0001) No. starts: 120
2570     (PID.TID 0000.0001) No. stops: 120
2571 jmc 1.1 (PID.TID 0000.0001) Seconds in section "DO_THE_MODEL_IO [FORWARD_STEP]":
2572 jmc 1.9 (PID.TID 0000.0001) User time: 0.24000000000000909
2573     (PID.TID 0000.0001) System time: 1.00000000000000089E-002
2574     (PID.TID 0000.0001) Wall clock time: 0.36638045310974121
2575     (PID.TID 0000.0001) No. starts: 120
2576     (PID.TID 0000.0001) No. stops: 120
2577     (PID.TID 0000.0001) Seconds in section "DO_WRITE_PICKUP [FORWARD_STEP]":
2578     (PID.TID 0000.0001) User time: 6.00000000000022737E-002
2579 jmc 1.8 (PID.TID 0000.0001) System time: 5.00000000000000444E-002
2580 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 0.20424556732177734
2581     (PID.TID 0000.0001) No. starts: 120
2582     (PID.TID 0000.0001) No. stops: 120
2583 jmc 1.1 (PID.TID 0000.0001) Seconds in section "CTRL_PACK [THE_MODEL_MAIN]":
2584 jmc 1.9 (PID.TID 0000.0001) User time: 6.00000000000022737E-002
2585     (PID.TID 0000.0001) System time: 0.0000000000000000
2586     (PID.TID 0000.0001) Wall clock time: 8.15470218658447266E-002
2587 jmc 1.6 (PID.TID 0000.0001) No. starts: 1
2588     (PID.TID 0000.0001) No. stops: 1
2589 jmc 1.1 (PID.TID 0000.0001) Seconds in section "CTRL_PACK [THE_MODEL_MAIN]":
2590 jmc 1.7 (PID.TID 0000.0001) User time: 6.00000000000022737E-002
2591 jmc 1.9 (PID.TID 0000.0001) System time: 9.99999999999995337E-003
2592     (PID.TID 0000.0001) Wall clock time: 5.83710670471191406E-002
2593 jmc 1.6 (PID.TID 0000.0001) No. starts: 1
2594     (PID.TID 0000.0001) No. stops: 1
2595 jmc 1.1 (PID.TID 0000.0001) Seconds in section "GRDCHK_MAIN [THE_MODEL_MAIN]":
2596 jmc 1.9 (PID.TID 0000.0001) User time: 63.820000000000007
2597     (PID.TID 0000.0001) System time: 1.00000000000000089E-002
2598     (PID.TID 0000.0001) Wall clock time: 64.022987842559814
2599 jmc 1.6 (PID.TID 0000.0001) No. starts: 1
2600     (PID.TID 0000.0001) No. stops: 1
2601 jmc 1.1 (PID.TID 0000.0001) Seconds in section "INITIALISE_VARIA [THE_MAIN_LOOP]":
2602 jmc 1.9 (PID.TID 0000.0001) User time: 0.31999999999999318
2603 jmc 1.6 (PID.TID 0000.0001) System time: 1.00000000000000089E-002
2604 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 0.31936240196228027
2605 jmc 1.6 (PID.TID 0000.0001) No. starts: 4
2606     (PID.TID 0000.0001) No. stops: 4
2607 jmc 1.1 (PID.TID 0000.0001) Seconds in section "MAIN LOOP [THE_MAIN_LOOP]":
2608 jmc 1.9 (PID.TID 0000.0001) User time: 63.429999999999978
2609     (PID.TID 0000.0001) System time: 0.0000000000000000
2610     (PID.TID 0000.0001) Wall clock time: 63.580344915390015
2611 jmc 1.6 (PID.TID 0000.0001) No. starts: 4
2612     (PID.TID 0000.0001) No. stops: 4
2613 jmc 1.9 (PID.TID 0000.0001) Seconds in section "MAIN_DO_LOOP [THE_MAIN_LOOP]":
2614     (PID.TID 0000.0001) User time: 62.530000000000058
2615     (PID.TID 0000.0001) System time: 0.0000000000000000
2616     (PID.TID 0000.0001) Wall clock time: 62.677838802337646
2617     (PID.TID 0000.0001) No. starts: 80
2618     (PID.TID 0000.0001) No. stops: 80
2619 jmc 1.1 (PID.TID 0000.0001) Seconds in section "COST_FINAL [ADJOINT SPIN-DOWN]":
2620 jmc 1.6 (PID.TID 0000.0001) User time: 0.0000000000000000
2621     (PID.TID 0000.0001) System time: 0.0000000000000000
2622 jmc 1.9 (PID.TID 0000.0001) Wall clock time: 3.49283218383789063E-004
2623 jmc 1.6 (PID.TID 0000.0001) No. starts: 4
2624     (PID.TID 0000.0001) No. stops: 4
2625 jmc 1.1 (PID.TID 0000.0001) // ======================================================
2626     (PID.TID 0000.0001) // Tile <-> Tile communication statistics
2627     (PID.TID 0000.0001) // ======================================================
2628     (PID.TID 0000.0001) // o Tile number: 000001
2629     (PID.TID 0000.0001) // No. X exchanges = 0
2630     (PID.TID 0000.0001) // Max. X spins = 0
2631     (PID.TID 0000.0001) // Min. X spins = 1000000000
2632     (PID.TID 0000.0001) // Total. X spins = 0
2633     (PID.TID 0000.0001) // Avg. X spins = 0.00E+00
2634     (PID.TID 0000.0001) // No. Y exchanges = 0
2635     (PID.TID 0000.0001) // Max. Y spins = 0
2636     (PID.TID 0000.0001) // Min. Y spins = 1000000000
2637     (PID.TID 0000.0001) // Total. Y spins = 0
2638     (PID.TID 0000.0001) // Avg. Y spins = 0.00E+00
2639     (PID.TID 0000.0001) // o Thread number: 000001
2640 jmc 1.9 (PID.TID 0000.0001) // No. barriers = 117660
2641 jmc 1.1 (PID.TID 0000.0001) // Max. barrier spins = 1
2642     (PID.TID 0000.0001) // Min. barrier spins = 1
2643 jmc 1.9 (PID.TID 0000.0001) // Total barrier spins = 117660
2644 jmc 1.1 (PID.TID 0000.0001) // Avg. barrier spins = 1.00E+00
2645 jmc 1.7 PROGRAM MAIN: Execution ended Normally

  ViewVC Help
Powered by ViewVC 1.1.22