(PID.TID 0000.0001) 
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) //                      MITgcm UV
(PID.TID 0000.0001) //                      =========
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // execution environment starting up...
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // MITgcmUV version:  checkpoint68q
(PID.TID 0000.0001) // Build user:        jmc
(PID.TID 0000.0001) // Build host:        jaures.mit.edu
(PID.TID 0000.0001) // Build date:        Tue Jul 25 03:57:55 PM EDT 2023
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Execution Environment parameter file "eedata"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># Example "eedata" file
(PID.TID 0000.0001) ># Lines beginning "#" are comments
(PID.TID 0000.0001) ># nTx - No. threads per process in X
(PID.TID 0000.0001) ># nTy - No. threads per process in Y
(PID.TID 0000.0001) > &EEPARMS
(PID.TID 0000.0001) > nTx=1,
(PID.TID 0000.0001) > nTy=1,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) ># Note: Some systems use & as the namelist terminator (as shown here).
(PID.TID 0000.0001) >#       Other systems use a / character.
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Computational Grid Specification ( see files "SIZE.h" )
(PID.TID 0000.0001) //                                  ( and "eedata"       )
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)      nPx =    1 ; /* No. processes in X */
(PID.TID 0000.0001)      nPy =    1 ; /* No. processes in Y */
(PID.TID 0000.0001)      nSx =    2 ; /* No. tiles in X per process */
(PID.TID 0000.0001)      nSy =    2 ; /* No. tiles in Y per process */
(PID.TID 0000.0001)      sNx =   45 ; /* Tile size in X */
(PID.TID 0000.0001)      sNy =   20 ; /* Tile size in Y */
(PID.TID 0000.0001)      OLx =    3 ; /* Tile overlap distance in X */
(PID.TID 0000.0001)      OLy =    3 ; /* Tile overlap distance in Y */
(PID.TID 0000.0001)      nTx =    1 ; /* No. threads in X per process */
(PID.TID 0000.0001)      nTy =    1 ; /* No. threads in Y per process */
(PID.TID 0000.0001)       Nr =   15 ; /* No. levels in the vertical   */
(PID.TID 0000.0001)       Nx =   90 ; /* Total domain size in X ( = nPx*nSx*sNx ) */
(PID.TID 0000.0001)       Ny =   40 ; /* Total domain size in Y ( = nPy*nSy*sNy ) */
(PID.TID 0000.0001)   nTiles =    4 ; /* Total no. tiles per process ( = nSx*nSy ) */
(PID.TID 0000.0001)   nProcs =    1 ; /* Total no. processes ( = nPx*nPy ) */
(PID.TID 0000.0001) nThreads =    1 ; /* Total no. threads per process ( = nTx*nTy ) */
(PID.TID 0000.0001) usingMPI =    F ; /* Flag used to control whether MPI is in use */
(PID.TID 0000.0001)                   /*  note: To execute a program with MPI calls */
(PID.TID 0000.0001)                   /*  it must be launched appropriately e.g     */
(PID.TID 0000.0001)                   /*  "mpirun -np 64 ......"                    */
(PID.TID 0000.0001) useCoupler=   F ; /* Flag used to control communications with   */
(PID.TID 0000.0001)                   /*  other model components, through a coupler */
(PID.TID 0000.0001) useNest2W_parent =    F ;/* Control 2-W Nesting comm */
(PID.TID 0000.0001) useNest2W_child  =    F ;/* Control 2-W Nesting comm */
(PID.TID 0000.0001) debugMode =    F ; /* print debug msg. (sequence of S/R calls)  */
(PID.TID 0000.0001) printMapIncludesZeros=    F ; /* print zeros in Std.Output maps */
(PID.TID 0000.0001) maxLengthPrt1D=   65 /* maxLength of 1D array printed to StdOut */
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // Mapping of tiles to threads
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // -o- Thread   1, tiles (   1:   2,   1:   2)
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // Tile <-> Tile connectvity table
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // Tile number: 000001 (process no. = 000000)
(PID.TID 0000.0001) //        WEST: Tile = 000002, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000001
(PID.TID 0000.0001) //        EAST: Tile = 000002, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000001
(PID.TID 0000.0001) //       SOUTH: Tile = 000003, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000002
(PID.TID 0000.0001) //       NORTH: Tile = 000003, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000002
(PID.TID 0000.0001) // Tile number: 000002 (process no. = 000000)
(PID.TID 0000.0001) //        WEST: Tile = 000001, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001) //        EAST: Tile = 000001, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001) //       SOUTH: Tile = 000004, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000002
(PID.TID 0000.0001) //       NORTH: Tile = 000004, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000002
(PID.TID 0000.0001) // Tile number: 000003 (process no. = 000000)
(PID.TID 0000.0001) //        WEST: Tile = 000004, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000002
(PID.TID 0000.0001) //        EAST: Tile = 000004, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000002
(PID.TID 0000.0001) //       SOUTH: Tile = 000001, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001) //       NORTH: Tile = 000001, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001) // Tile number: 000004 (process no. = 000000)
(PID.TID 0000.0001) //        WEST: Tile = 000003, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000002
(PID.TID 0000.0001) //        EAST: Tile = 000003, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000002
(PID.TID 0000.0001) //       SOUTH: Tile = 000002, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000001
(PID.TID 0000.0001) //       NORTH: Tile = 000002, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000001
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  INI_PARMS: opening model parameter file "data"
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># ====================
(PID.TID 0000.0001) ># | Model parameters |
(PID.TID 0000.0001) ># ====================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># Continuous equation parameters
(PID.TID 0000.0001) > &PARM01
(PID.TID 0000.0001) > tRef = 15*20.,
(PID.TID 0000.0001) > sRef = 15*35.,
(PID.TID 0000.0001) > viscAr=1.E-3,
(PID.TID 0000.0001) > viscAh=5.E5,
(PID.TID 0000.0001) > diffKhT=0.0,
(PID.TID 0000.0001) > diffKrT=3.E-5,
(PID.TID 0000.0001) > diffKhS=0.0,
(PID.TID 0000.0001) > diffKrS=3.E-5,
(PID.TID 0000.0001) > rhoConst=1035.,
(PID.TID 0000.0001) > rotationPeriod=86400.,
(PID.TID 0000.0001) > gravity=9.81,
(PID.TID 0000.0001) > eosType = 'JMD95Z',
(PID.TID 0000.0001) > ivdc_kappa=100.,
(PID.TID 0000.0001) > implicitDiffusion=.TRUE.,
(PID.TID 0000.0001) > allowFreezing=.TRUE.,
(PID.TID 0000.0001) > useRealFreshWaterFlux=.TRUE.,
(PID.TID 0000.0001) > useCDscheme=.TRUE.,
(PID.TID 0000.0001) ># turn on looped cells
(PID.TID 0000.0001) > hFacMin=.05,
(PID.TID 0000.0001) > hFacMindr=50.,
(PID.TID 0000.0001) ># set precision of data files
(PID.TID 0000.0001) > readBinaryPrec=32,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Elliptic solver parameters
(PID.TID 0000.0001) > &PARM02
(PID.TID 0000.0001) > cg2dMaxIters=1000,
(PID.TID 0000.0001) > cg2dTargetResidual=1.E-19,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Time stepping parameters
(PID.TID 0000.0001) > &PARM03
(PID.TID 0000.0001) > nIter0 =     0,
(PID.TID 0000.0001) > nTimeSteps = 4,
(PID.TID 0000.0001) ># 100 years of integration will yield a reasonable flow field
(PID.TID 0000.0001) ># startTime  =          0.,
(PID.TID 0000.0001) ># endTime    = 3110400000.,
(PID.TID 0000.0001) > deltaTMom = 1200.0,
(PID.TID 0000.0001) > tauCD =     321428.,
(PID.TID 0000.0001) > deltaTtracer= 43200.0,
(PID.TID 0000.0001) > deltaTClock = 43200.0,
(PID.TID 0000.0001) ># if you are using a version later than checkpoint45d on the main branch
(PID.TID 0000.0001) ># you can uncomment the following line and increase the time step
(PID.TID 0000.0001) ># deltaTtracer and deltaTClock to 172800.0 as well to speed up the
(PID.TID 0000.0001) ># asynchronous time stepping
(PID.TID 0000.0001) ># deltaTfreesurf = 172800.0,
(PID.TID 0000.0001) > abEps = 0.1,
(PID.TID 0000.0001) > pChkptFreq= 311040000.,
(PID.TID 0000.0001) > dumpFreq       = 2592000.,
(PID.TID 0000.0001) > adjDumpFreq    = 2592000.,
(PID.TID 0000.0001) > monitorFreq    = 0.,
(PID.TID 0000.0001) > adjMonitorFreq = 0.,
(PID.TID 0000.0001) ># 2 months restoring timescale for temperature
(PID.TID 0000.0001) > tauThetaClimRelax =  5184000.0,
(PID.TID 0000.0001) ># 6 months restoring timescale for salinity
(PID.TID 0000.0001) > tauSaltClimRelax = 15552000.0,
(PID.TID 0000.0001) > periodicExternalForcing=.TRUE.,
(PID.TID 0000.0001) > externForcingPeriod=2592000.,
(PID.TID 0000.0001) > externForcingCycle=31104000.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Gridding parameters
(PID.TID 0000.0001) > &PARM04
(PID.TID 0000.0001) > usingSphericalPolarGrid=.TRUE.,
(PID.TID 0000.0001) > delR= 50.,  70., 100., 140., 190.,
(PID.TID 0000.0001) >      240., 290., 340., 390., 440.,
(PID.TID 0000.0001) >      490., 540., 590., 640., 690.,
(PID.TID 0000.0001) > ygOrigin=-80.,
(PID.TID 0000.0001) > dySpacing=4.,
(PID.TID 0000.0001) > dxSpacing=4.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Input datasets
(PID.TID 0000.0001) > &PARM05
(PID.TID 0000.0001) > bathyFile=      'bathymetry.bin',
(PID.TID 0000.0001) > hydrogThetaFile='lev_t.bin',
(PID.TID 0000.0001) > hydrogSaltFile= 'lev_s.bin',
(PID.TID 0000.0001) > zonalWindFile=  'trenberth_taux.bin',
(PID.TID 0000.0001) > meridWindFile=  'trenberth_tauy.bin',
(PID.TID 0000.0001) > thetaClimFile=  'lev_sst.bin',
(PID.TID 0000.0001) > saltClimFile=   'lev_sss.bin',
(PID.TID 0000.0001) > surfQnetFile=   'ncep_qnet.bin',
(PID.TID 0000.0001) ># fresh water flux is turned off, uncomment next line to turn on
(PID.TID 0000.0001) ># (not recommened together with surface salinity restoring)
(PID.TID 0000.0001) ># EmPmRFile=      'ncep_emp.bin',
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  INI_PARMS ; starts to read PARM01
(PID.TID 0000.0001)  INI_PARMS ; read PARM01 : OK
(PID.TID 0000.0001) S/R INI_PARMS: No request for barotropic solver
(PID.TID 0000.0001) S/R INI_PARMS: => Use implicitFreeSurface as default
(PID.TID 0000.0001)  INI_PARMS ; starts to read PARM02
(PID.TID 0000.0001)  INI_PARMS ; read PARM02 : OK
(PID.TID 0000.0001)  INI_PARMS ; starts to read PARM03
(PID.TID 0000.0001)  INI_PARMS ; read PARM03 : OK
(PID.TID 0000.0001)  INI_PARMS ; starts to read PARM04
(PID.TID 0000.0001)  INI_PARMS ; read PARM04 : OK
(PID.TID 0000.0001)  INI_PARMS ; starts to read PARM05
(PID.TID 0000.0001)  INI_PARMS ; read PARM05 : OK
(PID.TID 0000.0001)  INI_PARMS: finished reading file "data"
(PID.TID 0000.0001)  PACKAGES_BOOT: opening data.pkg
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.pkg
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.pkg"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># ********
(PID.TID 0000.0001) ># Packages
(PID.TID 0000.0001) ># ********
(PID.TID 0000.0001) > &PACKAGES
(PID.TID 0000.0001) > useGMRedi    = .TRUE.,
(PID.TID 0000.0001) > useGGL90     = .FALSE.,
(PID.TID 0000.0001) > useKPP       = .FALSE.,
(PID.TID 0000.0001) > useGrdchk    = .TRUE.,
(PID.TID 0000.0001) > useMNC       = .FALSE.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  PACKAGES_BOOT: finished reading data.pkg
(PID.TID 0000.0001)  PACKAGES_BOOT: On/Off package Summary
 --------  pkgs with a standard "usePKG" On/Off switch in "data.pkg":  --------
 pkg/ggl90                compiled but not used ( useGGL90                 = F )
 pkg/kpp                  compiled but not used ( useKPP                   = F )
 pkg/gmredi               compiled   and   used ( useGMRedi                = T )
 pkg/autodiff             compiled   and   used ( useAUTODIFF              = T )
 pkg/grdchk               compiled   and   used ( useGrdchk                = T )
 pkg/ctrl                 compiled   and   used ( useCTRL                  = T )
 -------- pkgs without standard "usePKG" On/Off switch in "data.pkg":  --------
 pkg/generic_advdiff      compiled   and   used ( useGAD                   = T )
 pkg/mom_common           compiled   and   used ( momStepping              = T )
 pkg/mom_vecinv           compiled but not used ( +vectorInvariantMomentum = F )
 pkg/mom_fluxform         compiled   and   used ( & not vectorInvariantMom = T )
 pkg/cd_code              compiled   and   used ( useCDscheme              = T )
 pkg/debug                compiled but not used ( debugMode                = F )
 pkg/rw                   compiled   and   used
 pkg/mdsio                compiled   and   used
 pkg/autodiff             compiled   and   used
 pkg/cost                 compiled   and   used
(PID.TID 0000.0001)  PACKAGES_BOOT: End of package Summary
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  GM_READPARMS: opening data.gmredi
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.gmredi
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.gmredi"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># GM+Redi package parameters:
(PID.TID 0000.0001) >#     GM_Small_Number  :: epsilon used in computing the slope
(PID.TID 0000.0001) >#     GM_slopeSqCutoff :: slope^2 cut-off value
(PID.TID 0000.0001) >
(PID.TID 0000.0001) >#-from MOM :
(PID.TID 0000.0001) ># GM_background_K: 	G & Mc.W  diffusion coefficient
(PID.TID 0000.0001) ># GM_maxSlope    :	max slope of isopycnals
(PID.TID 0000.0001) ># GM_Scrit       :	transition for scaling diffusion coefficient
(PID.TID 0000.0001) ># GM_Sd          :	half width scaling for diffusion coefficient
(PID.TID 0000.0001) ># GM_taper_scheme:	slope clipping or one of the tapering schemes
(PID.TID 0000.0001) ># GM_Kmin_horiz  :	horizontal diffusion minimum value
(PID.TID 0000.0001) >
(PID.TID 0000.0001) >#-Option parameters (needs to "define" options in GMREDI_OPTIONS.h")
(PID.TID 0000.0001) ># GM_isopycK     :	isopycnal diffusion coefficient (default=GM_background_K)
(PID.TID 0000.0001) ># GM_AdvForm     :	turn on GM Advective form       (default=Skew flux form)
(PID.TID 0000.0001) >
(PID.TID 0000.0001) > &GM_PARM01
(PID.TID 0000.0001) >  GM_Small_Number  = 1.D-12,
(PID.TID 0000.0001) >  GM_slopeSqCutoff = 1.D+08,
(PID.TID 0000.0001) >  GM_AdvForm         = .FALSE.,
(PID.TID 0000.0001) >  GM_isopycK         = 1.0D+3,
(PID.TID 0000.0001) >  GM_background_K    = 1.0D+3,
(PID.TID 0000.0001) >  GM_taper_scheme    = 'dm95',
(PID.TID 0000.0001) >  GM_maxSlope        = 1.D-2,
(PID.TID 0000.0001) >  GM_Kmin_horiz      = 50.,
(PID.TID 0000.0001) >  GM_Scrit           = 4.D-3,
(PID.TID 0000.0001) >  GM_Sd              = 1.D-3,
(PID.TID 0000.0001) >#  GM_Visbeck_alpha   = 1.5D-2,
(PID.TID 0000.0001) >#  GM_Visbeck_alpha   = 0.D0,
(PID.TID 0000.0001) >#  GM_Visbeck_length  = 2.D+5,
(PID.TID 0000.0001) >#  GM_Visbeck_depth   = 1.D+3,
(PID.TID 0000.0001) >#  GM_Visbeck_maxval_K= 2.5D+3,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) >
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  GM_READPARMS: finished reading data.gmredi
(PID.TID 0000.0001) AUTODIFF_READPARMS: opening data.autodiff
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.autodiff
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.autodiff"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># =========================
(PID.TID 0000.0001) ># pkg AUTODIFF parameters :
(PID.TID 0000.0001) ># =========================
(PID.TID 0000.0001) >#  inAdExact :: get an exact adjoint (no approximation) (def=.True.)
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &AUTODIFF_PARM01
(PID.TID 0000.0001) ># inAdExact = .FALSE.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001) AUTODIFF_READPARMS: finished reading data.autodiff
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) // AUTODIFF parameters :
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) inAdExact = /* get an exact adjoint (no approximation) */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useApproxAdvectionInAdMode = /* approximate AD-advection */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cg2dFullAdjoint = /* use full hand written cg2d adjoint (no approximation) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useKPPinAdMode = /* use KPP in adjoint mode */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useGMRediInAdMode = /* use GMRedi in adjoint mode */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useSEAICEinAdMode = /* use SEAICE in adjoint mode */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useGGL90inAdMode = /* use GGL90 in adjoint mode */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useSALT_PLUMEinAdMode = /* use SALT_PLUME in adjoint mode */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dumpAdVarExch = /* control adexch before dumpinp */
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) mon_AdVarExch = /* control adexch before monitor */
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) viscFacInFw = /* viscosity factor for forward model */
(PID.TID 0000.0001)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) viscFacInAd = /* viscosity factor for adjoint */
(PID.TID 0000.0001)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) SIregFacInAd = /* sea ice factor for adjoint model */
(PID.TID 0000.0001)                 1.234567000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) SIregFacInFw = /* sea ice factor for forward model */
(PID.TID 0000.0001)                 1.234567000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) 
(PID.TID 0000.0001) OPTIM_READPARMS: opening data.optim
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.optim
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.optim"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># ********************************
(PID.TID 0000.0001) ># Off-line optimization parameters
(PID.TID 0000.0001) ># ********************************
(PID.TID 0000.0001) > &OPTIM
(PID.TID 0000.0001) > optimcycle=0,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001) OPTIM_READPARMS: finished reading data.optim
(PID.TID 0000.0001) CTRL_READPARMS: opening data.ctrl
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.ctrl
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.ctrl"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) ># ECCO controlvariables
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) > &CTRL_NML
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) ># names for ctrl_pack/unpack
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) > &CTRL_PACKNAMES
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) ># names for CTRL_GENARR, CTRL_GENTIM
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) > &CTRL_NML_GENARR
(PID.TID 0000.0001) ># the bounds are commented out to reproduce old results with
(PID.TID 0000.0001) ># the CPP-flag DISABLE_CTRL_THETA_LIMIT defined
(PID.TID 0000.0001) > xx_genarr3d_file(1)       = 'xx_theta',
(PID.TID 0000.0001) > xx_genarr3d_weight(1)     = 'ones_64b.bin',
(PID.TID 0000.0001) >#xx_genarr3d_bounds(1:5,1) = -2.,-1.9,39.,40.,5.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_genarr3d_file(2)       = 'xx_salt',
(PID.TID 0000.0001) > xx_genarr3d_weight(2)     = 'ones_64b.bin',
(PID.TID 0000.0001) >#xx_genarr3d_bounds(1:5,2) = 29.,29.5,40.5,41.,5.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_genarr2d_file(1)       = 'xx_qnet',
(PID.TID 0000.0001) > xx_genarr2d_weight(1)     = 'ones_64b.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_genarr2d_file(2)       = 'xx_empmr',
(PID.TID 0000.0001) > xx_genarr2d_weight(2)     = 'ones_64b.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_genarr2d_file(3)       = 'xx_fu',
(PID.TID 0000.0001) > xx_genarr2d_weight(3)     = 'ones_64b.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_genarr2d_file(4)       = 'xx_fv',
(PID.TID 0000.0001) > xx_genarr2d_weight(4)     = 'ones_64b.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001) CTRL_READPARMS: finished reading data.ctrl
(PID.TID 0000.0001) read-write ctrl files from current run directory
(PID.TID 0000.0001) COST_READPARMS: opening data.cost
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.cost
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.cost"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># ******************
(PID.TID 0000.0001) ># ECCO cost function
(PID.TID 0000.0001) ># ******************
(PID.TID 0000.0001) > &COST_NML
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > mult_tracer  = 1.,
(PID.TID 0000.0001) > mult_test    = 1.,
(PID.TID 0000.0001) > mult_atl     = 1.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001) COST_READPARMS: finished reading data.cost
(PID.TID 0000.0001) GRDCHK_READPARMS: opening data.grdchk
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.grdchk
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.grdchk"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># *******************
(PID.TID 0000.0001) ># ECCO gradient check
(PID.TID 0000.0001) ># *******************
(PID.TID 0000.0001) > &GRDCHK_NML
(PID.TID 0000.0001) > grdchk_eps       = 1.d-2,
(PID.TID 0000.0001) > iGloTile         = 2,
(PID.TID 0000.0001) > jGloTile         = 2,
(PID.TID 0000.0001) > iGloPos          = 26,
(PID.TID 0000.0001) >#jGloPos          = 19,
(PID.TID 0000.0001) ># This original setting would start the gradient check for temperatures
(PID.TID 0000.0001) ># that are below the freezing point in the inital conditions. Because of
(PID.TID 0000.0001) ># allowFreezing=True in data&PARM01 and the unphysical nature of this
(PID.TID 0000.0001) ># freezing parameterisation, the FD gradients depend very much on the
(PID.TID 0000.0001) ># whether they are evaluated above and below freezing or not leading
(PID.TID 0000.0001) ># to inconsistent results. Moving the starting point to jGloPos=17
(PID.TID 0000.0001) ># avoids this. This is similar to input_ad but here we use jGloTile = 2.
(PID.TID 0000.0001) > jGloPos          = 17,
(PID.TID 0000.0001) > kGloPos          = 1,
(PID.TID 0000.0001) >### nbeg             = 1,
(PID.TID 0000.0001) > nstep            = 1,
(PID.TID 0000.0001) > nend             = 7,
(PID.TID 0000.0001) > grdchkvarindex   = 201,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001) GRDCHK_READPARMS: finished reading data.grdchk
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Gradient check configuration  >>> START <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
(PID.TID 0000.0001)   grdchkvarindex :                        201
(PID.TID 0000.0001)   eps:                              0.100E-01
(PID.TID 0000.0001)   First location:                           0
(PID.TID 0000.0001)   Last location:                            7
(PID.TID 0000.0001)   Increment:                                1
(PID.TID 0000.0001)   grdchkWhichProc:                          0
(PID.TID 0000.0001)   iLocTile =       2  ,    jLocTile =       2
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Gradient check configuration  >>> END <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
(PID.TID 0000.0001) SET_PARMS: done
(PID.TID 0000.0001) Enter INI_VERTICAL_GRID: setInterFDr=    T ; setCenterDr=    F
(PID.TID 0000.0001) GAD_INIT_FIXED: GAD_OlMinSize=  1  0  1
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) // GAD parameters :
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) tempAdvScheme =   /* Temp. Horiz.Advection scheme selector */
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempVertAdvScheme =   /* Temp. Vert. Advection scheme selector */
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempMultiDimAdvec =   /* use Muti-Dim Advec method for Temp */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempSOM_Advection = /* use 2nd Order Moment Advection for Temp */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) AdamsBashforthGt = /* apply Adams-Bashforth extrapolation on Gt */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) AdamsBashforth_T = /* apply Adams-Bashforth extrapolation on Temp */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltAdvScheme =   /* Salt. Horiz.advection scheme selector */
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltVertAdvScheme =   /* Salt. Vert. Advection scheme selector */
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltMultiDimAdvec =   /* use Muti-Dim Advec method for Salt */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltSOM_Advection = /* use 2nd Order Moment Advection for Salt */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) AdamsBashforthGs = /* apply Adams-Bashforth extrapolation on Gs */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) AdamsBashforth_S = /* apply Adams-Bashforth extrapolation on Salt */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) ctrl-wet 1:    nvarlength =        67878
(PID.TID 0000.0001) ctrl-wet 2: surface wet C =          669
(PID.TID 0000.0001) ctrl-wet 3: surface wet W =          640
(PID.TID 0000.0001) ctrl-wet 4: surface wet S =          610
(PID.TID 0000.0001) ctrl-wet 4a:surface wet V =            0
(PID.TID 0000.0001) ctrl-wet 5: 3D wet points =         8651
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     1           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     2           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     3           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     4           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     5           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     6           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     7           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     8           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     9           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    10           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    11           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    12           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    13           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    14           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    15           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    16           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    17           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    18           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    19           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    20           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    21           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    22           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    23           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    24           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    25           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    26           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    27           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    28           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    29           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    30           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    31           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    32           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    33           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    34           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    35           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    36           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    37           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    38           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    39           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    40           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    41           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    42           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    43           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    44           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    45           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    46           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    47           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    48           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    49           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    50           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    51           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    52           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    53           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    54           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    55           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    56           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    57           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    58           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    59           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    60           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    61           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    62           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    63           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    64           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    65           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    66           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    67           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    68           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    69           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    70           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    71           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    72           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    73           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    74           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    75           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    76           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    77           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    78           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    79           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    80           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    81           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    82           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    83           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    84           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    85           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    86           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    87           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    88           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    89           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    90           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    91           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    92           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    93           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    94           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    95           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    96           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    97           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    98           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    99           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   100           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   101           1
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   102           1
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   103           1
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   104           1
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   105           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   106           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   107           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   108           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   109           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   110           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   111           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   112           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   113           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   114           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   115           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   116           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   117           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   118           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   119           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   120           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   121           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   122           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   123           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   124           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   125           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   126           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   127           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   128           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   129           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   130           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   131           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   132           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   133           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   134           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   135           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   136           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   137           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   138           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   139           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   140           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   141           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   142           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   143           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   144           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   145           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   146           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   147           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   148           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   149           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   150           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   151           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   152           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   153           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   154           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   155           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   156           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   157           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   158           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   159           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   160           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   161           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   162           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   163           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   164           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   165           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   166           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   167           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   168           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   169           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   170           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   171           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   172           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   173           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   174           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   175           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   176           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   177           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   178           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   179           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   180           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   181           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   182           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   183           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   184           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   185           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   186           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   187           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   188           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   189           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   190           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   191           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   192           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   193           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   194           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   195           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   196           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   197           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   198           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   199           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   200           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   201           1
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   202           1
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   203           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   204           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   205           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   206           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   207           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   208           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   209           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   210           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   211           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   212           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   213           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   214           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   215           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   216           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   217           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   218           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   219           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   220           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   221           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   222           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   223           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   224           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   225           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   226           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   227           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   228           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   229           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   230           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   231           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   232           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   233           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   234           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   235           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   236           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   237           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   238           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   239           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   240           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   241           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   242           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   243           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   244           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   245           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   246           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   247           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   248           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   249           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   250           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   251           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   252           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   253           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   254           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   255           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   256           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   257           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   258           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   259           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   260           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   261           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   262           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   263           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   264           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   265           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   266           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   267           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   268           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   269           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   270           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   271           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   272           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   273           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   274           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   275           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   276           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   277           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   278           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   279           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   280           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   281           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   282           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   283           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   284           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   285           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   286           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   287           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   288           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   289           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   290           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   291           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   292           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   293           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   294           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   295           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   296           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   297           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   298           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   299           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   300           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   301           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   302           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   303           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   304           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   305           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   306           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   307           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   308           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   309           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   310           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   311           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   312           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   313           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   314           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   315           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   316           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   317           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   318           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   319           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   320           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   321           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   322           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   323           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   324           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   325           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   326           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   327           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   328           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   329           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   330           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   331           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   332           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   333           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   334           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   335           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   336           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   337           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   338           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   339           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   340           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   341           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   342           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   343           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   344           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   345           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   346           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   347           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   348           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   349           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   350           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   351           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   352           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   353           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   354           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   355           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   356           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   357           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   358           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   359           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   360           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   361           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   362           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   363           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   364           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   365           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   366           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   367           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   368           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   369           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   370           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   371           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   372           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   373           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   374           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   375           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   376           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   377           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   378           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   379           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   380           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   381           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   382           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   383           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   384           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   385           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   386           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   387           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   388           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   389           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   390           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   391           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   392           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   393           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   394           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   395           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   396           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   397           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   398           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   399           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =   400           0
(PID.TID 0000.0001) ctrl-wet 7: flux         17302
(PID.TID 0000.0001) ctrl-wet 8: atmos        17302
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl-wet 13: global nvarlength for Nr =   15       67878
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    1        2315        2149        2206           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    2        2315        2149        2206           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    3        2254        2102        2146           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    4        2215        2058        2104           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    5        2178        2027        2070           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    6        2142        1987        2029           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    7        2114        1959        2004           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    8        2076        1918        1959           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    9        2048        1887        1925           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   10        1999        1831        1869           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   11        1948        1771        1808           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   12        1850        1653        1705           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   13        1655        1404        1458           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   14        1372        1118        1164           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   15         828         623         671           0
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl_init_wet: no. of control variables:            6
(PID.TID 0000.0001) ctrl_init_wet: control vector length:           67878
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // control vector configuration  >>> START <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  Total number of ocean points per tile:
(PID.TID 0000.0001)  --------------------------------------
(PID.TID 0000.0001)  snx*sny*nr =    13500
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  Number of ocean points per tile:
(PID.TID 0000.0001)  --------------------------------
(PID.TID 0000.0001)  bi,bj,#(c/s/w): 0001 0001 008651 007584 008016
(PID.TID 0000.0001)  bi,bj,#(c/s/w): 0002 0001 010278 009360 009834
(PID.TID 0000.0001)  bi,bj,#(c/s/w): 0001 0002 003504 003232 003049
(PID.TID 0000.0001)  bi,bj,#(c/s/w): 0002 0002 006876 006460 006425
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  Settings of generic controls:
(PID.TID 0000.0001)  -----------------------------
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  -> 2D control, genarr2d no.  1 is in use
(PID.TID 0000.0001)       file       = xx_qnet
(PID.TID 0000.0001)       weight     = ones_64b.bin
(PID.TID 0000.0001)       index      =  0101
(PID.TID 0000.0001)       ncvarindex =  0201
(PID.TID 0000.0001)  -> 2D control, genarr2d no.  2 is in use
(PID.TID 0000.0001)       file       = xx_empmr
(PID.TID 0000.0001)       weight     = ones_64b.bin
(PID.TID 0000.0001)       index      =  0102
(PID.TID 0000.0001)       ncvarindex =  0202
(PID.TID 0000.0001)  -> 2D control, genarr2d no.  3 is in use
(PID.TID 0000.0001)       file       = xx_fu
(PID.TID 0000.0001)       weight     = ones_64b.bin
(PID.TID 0000.0001)       index      =  0103
(PID.TID 0000.0001)       ncvarindex =  0203
(PID.TID 0000.0001)  -> 2D control, genarr2d no.  4 is in use
(PID.TID 0000.0001)       file       = xx_fv
(PID.TID 0000.0001)       weight     = ones_64b.bin
(PID.TID 0000.0001)       index      =  0104
(PID.TID 0000.0001)       ncvarindex =  0204
(PID.TID 0000.0001)  -> 3d control, genarr3d no.  1 is in use
(PID.TID 0000.0001)       file       = xx_theta
(PID.TID 0000.0001)       weight     = ones_64b.bin
(PID.TID 0000.0001)       index      =  0201
(PID.TID 0000.0001)       ncvarindex =  0301
(PID.TID 0000.0001)  -> 3d control, genarr3d no.  2 is in use
(PID.TID 0000.0001)       file       = xx_salt
(PID.TID 0000.0001)       weight     = ones_64b.bin
(PID.TID 0000.0001)       index      =  0202
(PID.TID 0000.0001)       ncvarindex =  0302
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // control vector configuration  >>> END <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
(PID.TID 0000.0001) INI_CG2D: CG2D normalisation factor =  6.5682677425711703E-05
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model configuration
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) //
(PID.TID 0000.0001) // "Physical" paramters ( PARM01 in namelist )
(PID.TID 0000.0001) //
(PID.TID 0000.0001) buoyancyRelation = /* Type of relation to get Buoyancy */
(PID.TID 0000.0001)               'OCEANIC'
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) fluidIsAir   =  /* fluid major constituent is Air */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) fluidIsWater =  /* fluid major constituent is Water */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingPCoords =  /* use p (or p*) vertical coordinate */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingZCoords =  /* use z (or z*) vertical coordinate */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tRef =   /* Reference temperature profile ( oC or K ) */
(PID.TID 0000.0001)    15 @  2.000000000000000E+01              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) sRef =   /* Reference salinity profile ( g/kg ) */
(PID.TID 0000.0001)    15 @  3.500000000000000E+01              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoRef =   /* Density vertical profile from (Ref,sRef)( kg/m^3 ) */
(PID.TID 0000.0001)                 1.024872626184147E+03,      /* K =  1 */
(PID.TID 0000.0001)                 1.025135462285008E+03,      /* K =  2 */
(PID.TID 0000.0001)                 1.025507198938228E+03,      /* K =  3 */
(PID.TID 0000.0001)                 1.026030780760464E+03,      /* K =  4 */
(PID.TID 0000.0001)                 1.026748377776259E+03,      /* K =  5 */
(PID.TID 0000.0001)                 1.027679406285166E+03,      /* K =  6 */
(PID.TID 0000.0001)                 1.028820735595355E+03,      /* K =  7 */
(PID.TID 0000.0001)                 1.030168558073105E+03,      /* K =  8 */
(PID.TID 0000.0001)                 1.031718419899614E+03,      /* K =  9 */
(PID.TID 0000.0001)                 1.033465256541184E+03,      /* K = 10 */
(PID.TID 0000.0001)                 1.035403432414885E+03,      /* K = 11 */
(PID.TID 0000.0001)                 1.037526784183520E+03,      /* K = 12 */
(PID.TID 0000.0001)                 1.039828667078104E+03,      /* K = 13 */
(PID.TID 0000.0001)                 1.042302003623418E+03,      /* K = 14 */
(PID.TID 0000.0001)                 1.044939334132512E+03       /* K = 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dBdrRef = /* Vertical grad. of reference buoyancy [(m/s/r)^2] */
(PID.TID 0000.0001)    15 @  0.000000000000000E+00              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useStrainTensionVisc= /* Use StrainTension Form of Viscous Operator */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useVariableVisc = /* Use variable horizontal viscosity */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useHarmonicVisc = /* Use harmonic horizontal viscosity */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useBiharmonicVisc= /* Use biharmonic horiz.  viscosity */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useSmag3D = /* Use isotropic 3-D Smagorinsky viscosity */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) viscAh  =   /* Lateral harmonic viscosity ( m^2/s ) */
(PID.TID 0000.0001)                 5.000000000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) viscA4  =   /* Lateral biharmonic viscosity ( m^4/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) no_slip_sides =  /* Viscous BCs: No-slip sides */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) sideDragFactor = /* side-drag scaling factor (non-dim) */
(PID.TID 0000.0001)                 2.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) viscArNr = /* vertical profile of vertical viscosity ( m^2/s )*/
(PID.TID 0000.0001)    15 @  1.000000000000000E-03              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) no_slip_bottom =  /* Viscous BCs: No-slip bottom */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) bottomVisc_pCell = /* Partial-cell in bottom Visc. BC */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) bottomDragLinear = /* linear bottom-drag coefficient ( m/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) bottomDragQuadratic = /* quadratic bottom-drag coefficient (-) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectBotDragQuadr = /* select quadratic bottom drag options */
(PID.TID 0000.0001)                      -1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKhT =   /* Laplacian diffusion of heat laterally ( m^2/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffK4T =   /* Biharmonic diffusion of heat laterally ( m^4/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKhS =   /* Laplacian diffusion of salt laterally ( m^2/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffK4S =   /* Biharmonic diffusion of salt laterally ( m^4/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrNrT = /* vertical profile of vertical diffusion of Temp ( m^2/s )*/
(PID.TID 0000.0001)    15 @  3.000000000000000E-05              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrNrS = /* vertical profile of vertical diffusion of Salt ( m^2/s )*/
(PID.TID 0000.0001)    15 @  3.000000000000000E-05              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrBL79surf = /* Surface diffusion for Bryan and Lewis 79 ( m^2/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrBL79deep = /* Deep diffusion for Bryan and Lewis 1979 ( m^2/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrBL79scl = /* Depth scale for Bryan and Lewis 1979 ( m ) */
(PID.TID 0000.0001)                 2.000000000000000E+02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrBL79Ho = /* Turning depth for Bryan and Lewis 1979 ( m ) */
(PID.TID 0000.0001)                -2.000000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) ivdc_kappa = /* Implicit Vertical Diffusivity for Convection ( m^2/s) */
(PID.TID 0000.0001)                 1.000000000000000E+02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hMixCriteria=  /* Criteria for mixed-layer diagnostic */
(PID.TID 0000.0001)                -8.000000000000000E-01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dRhoSmall =  /* Parameter for mixed-layer diagnostic */
(PID.TID 0000.0001)                 1.000000000000000E-06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hMixSmooth=  /* Smoothing parameter for mixed-layer diagnostic */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) eosType =  /* Type of Equation of State */
(PID.TID 0000.0001)               'JMD95Z'
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) eosRefP0 = /* Reference atmospheric pressure for EOS ( Pa ) */
(PID.TID 0000.0001)                 1.013250000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectP_inEOS_Zc = /* select pressure to use in EOS (0,1,2,3) */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     0= -g*rhoConst*z ; 1= pRef (from tRef,sRef); 2= Hyd P ; 3= Hyd+NH P
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) surf_pRef = /* Surface reference pressure ( Pa ) */
(PID.TID 0000.0001)                 1.013250000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) HeatCapacity_Cp =  /* Specific heat capacity ( J/kg/K ) */
(PID.TID 0000.0001)                 3.994000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) celsius2K = /* 0 degree Celsius converted to Kelvin ( K ) */
(PID.TID 0000.0001)                 2.731500000000000E+02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoConst  = /* Reference density (Boussinesq)  ( kg/m^3 ) */
(PID.TID 0000.0001)                 1.035000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoFacC = /* normalized Reference density @ cell-Center (-) */
(PID.TID 0000.0001)    15 @  1.000000000000000E+00              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoFacF = /* normalized Reference density @ W-Interface (-) */
(PID.TID 0000.0001)    16 @  1.000000000000000E+00              /* K =  1: 16 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoConstFresh = /* Fresh-water reference density ( kg/m^3 ) */
(PID.TID 0000.0001)                 1.035000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) gravity =   /* Gravitational acceleration ( m/s^2 ) */
(PID.TID 0000.0001)                 9.810000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) gBaro =   /* Barotropic gravity ( m/s^2 ) */
(PID.TID 0000.0001)                 9.810000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) gravFacC = /* gravity factor (vs surf.) @ cell-Center (-) */
(PID.TID 0000.0001)    15 @  1.000000000000000E+00              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) gravFacF = /* gravity factor (vs surf.) @ W-Interface (-) */
(PID.TID 0000.0001)    16 @  1.000000000000000E+00              /* K =  1: 16 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rotationPeriod =   /* Rotation Period ( s ) */
(PID.TID 0000.0001)                 8.640000000000000E+04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) omega =   /* Angular velocity ( rad/s ) */
(PID.TID 0000.0001)                 7.272205216643040E-05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) f0 =   /* Reference coriolis parameter ( 1/s ) */
(PID.TID 0000.0001)                 1.000000000000000E-04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) beta =   /* Beta ( 1/(m.s) ) */
(PID.TID 0000.0001)                 9.999999999999999E-12
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) fPrime =   /* Second coriolis parameter ( 1/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rigidLid =   /* Rigid lid on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) implicitFreeSurface =   /* Implicit free surface on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) freeSurfFac =   /* Implicit free surface factor */
(PID.TID 0000.0001)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) implicSurfPress =  /* Surface Pressure implicit factor (0-1) */
(PID.TID 0000.0001)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) implicDiv2DFlow =  /* Barot. Flow Div. implicit factor (0-1) */
(PID.TID 0000.0001)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) uniformLin_PhiSurf = /* use uniform Bo_surf on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) uniformFreeSurfLev = /* free-surface level-index is uniform */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) sIceLoadFac =  /* scale factor for sIceLoad (0-1) */
(PID.TID 0000.0001)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hFacMin =   /* minimum partial cell factor (hFac) */
(PID.TID 0000.0001)                 5.000000000000000E-02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hFacMinDr = /* minimum partial cell thickness ( m) */
(PID.TID 0000.0001)                 5.000000000000000E+01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) exactConserv =  /* Exact Volume Conservation on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) linFSConserveTr = /* Tracer correction for Lin Free Surface on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nonlinFreeSurf = /* Non-linear Free Surf. options (-1,0,1,2,3)*/
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)      -1,0= Off ; 1,2,3= On, 2=+rescale gU,gV, 3=+update cg2d solv.
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hFacInf =   /* lower threshold for hFac (nonlinFreeSurf only)*/
(PID.TID 0000.0001)                 2.000000000000000E-01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hFacSup =   /* upper threshold for hFac (nonlinFreeSurf only)*/
(PID.TID 0000.0001)                 2.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) select_rStar = /* r* Vertical coord. options (=0 r coord.; >0 uses r*)*/
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useRealFreshWaterFlux = /* Real Fresh Water Flux on/off flag*/
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) temp_EvPrRn = /* Temp. of Evap/Prec/R (UNSET=use local T)(oC)*/
(PID.TID 0000.0001)                 1.234567000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) salt_EvPrRn = /* Salin. of Evap/Prec/R (UNSET=use local S)(g/kg)*/
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectAddFluid = /* option for mass source/sink of fluid (=0: off) */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) temp_addMass = /* Temp. of addMass array (UNSET=use local T)(oC)*/
(PID.TID 0000.0001)                 1.234567000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) salt_addMass = /* Salin. of addMass array (UNSET=use local S)(g/kg)*/
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) convertFW2Salt = /* convert F.W. Flux to Salt Flux (-1=use local S)(g/kg)*/
(PID.TID 0000.0001)                -1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) use3Dsolver = /* use 3-D pressure solver on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nonHydrostatic =  /* Non-Hydrostatic on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nh_Am2 = /* Non-Hydrostatic terms scaling factor */
(PID.TID 0000.0001)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) implicitNHPress = /* Non-Hyd Pressure implicit factor (0-1)*/
(PID.TID 0000.0001)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectNHfreeSurf = /* Non-Hyd (free-)Surface option */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) quasiHydrostatic = /* Quasi-Hydrostatic on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) calc_wVelocity = /* vertical velocity calculation on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momStepping =  /* Momentum equation on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) vectorInvariantMomentum= /* Vector-Invariant Momentum on/off */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momAdvection =  /* Momentum advection on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momViscosity =  /* Momentum viscosity on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momImplVertAdv= /* Momentum implicit vert. advection on/off*/
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) implicitViscosity = /* Implicit viscosity on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectImplicitDrag= /* Implicit bot Drag options (0,1,2)*/
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     0= Expl. ; 1= Impl. on provis. Vel ; 2= Fully Impl (with surf.P)
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) metricTerms =  /* metric-Terms on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useNHMTerms = /* Non-Hydrostatic Metric-Terms on/off */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectCoriMap = /* Coriolis Map options (0,1,2,3)*/
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     0= f-Plane ; 1= Beta-Plane ; 2= Spherical ; 3= read from file
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) use3dCoriolis = /* 3-D Coriolis on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useCoriolis =  /* Coriolis on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useCDscheme =  /* CD scheme on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectCoriScheme= /* Scheme selector for Coriolis-Term */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)    = 0 : original discretization (simple averaging, no hFac)
(PID.TID 0000.0001)    = 1 : Wet-point averaging (Jamar & Ozer 1986)
(PID.TID 0000.0001)    = 2 : energy conserving scheme (no hFac weight)
(PID.TID 0000.0001)    = 3 : energy conserving scheme using Wet-point averaging
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momForcing =  /* Momentum forcing on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momTidalForcing = /* Momentum Tidal forcing on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momPressureForcing =  /* Momentum pressure term on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) implicitIntGravWave= /* Implicit Internal Gravity Wave flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) staggerTimeStep =    /* Stagger time stepping on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) doResetHFactors = /* reset thickness factors @ each time-step */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) multiDimAdvection =  /* enable/disable Multi-Dim Advection */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useMultiDimAdvec =   /* Multi-Dim Advection is/is-not used */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) implicitDiffusion = /* Implicit Diffusion on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempStepping =  /* Temperature equation on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempAdvection = /* Temperature advection on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempImplVertAdv = /* Temp. implicit vert. advection on/off */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempForcing  =  /* Temperature forcing on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) doThetaClimRelax = /* apply SST relaxation on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempIsActiveTr = /* Temp. is a dynamically Active Tracer */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltStepping =  /* Salinity equation on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltAdvection = /* Salinity advection on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltImplVertAdv = /* Sali. implicit vert. advection on/off */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltForcing  =  /* Salinity forcing on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) doSaltClimRelax = /* apply SSS relaxation on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltIsActiveTr = /* Salt  is a dynamically Active Tracer */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001)  readBinaryPrec = /* Precision used for reading binary files */
(PID.TID 0000.0001)                      32
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) writeBinaryPrec = /* Precision used for writing binary files */
(PID.TID 0000.0001)                      32
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001)  rwSuffixType =   /* select format of mds file suffix */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)    = 0 : myIter (I10.10) ;   = 1 : 100*myTime (100th sec) ;
(PID.TID 0000.0001)    = 2 : myTime (seconds);   = 3 : myTime/360 (10th of hr);
(PID.TID 0000.0001)    = 4 : myTime/3600 (hours)
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001)  globalFiles = /* write "global" (=not per tile) files */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001)  useSingleCpuIO = /* only master MPI process does I/O */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001)  useSingleCpuInput = /* only master process reads input */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) /* debLev[*]  : level of debug & auxiliary message printing */
(PID.TID 0000.0001) debLevZero =  0 ; /* level of disabled aux. msg printing */
(PID.TID 0000.0001)    debLevA =  1 ; /* level of minimum  aux. msg printing */
(PID.TID 0000.0001)    debLevB =  2 ; /* level of low aux. print (report read-file opening)*/
(PID.TID 0000.0001)    debLevC =  3 ; /* level of moderate debug prt (most pkgs debug msg) */
(PID.TID 0000.0001)    debLevD =  4 ; /* level of enhanced debug prt (add DEBUG_STATS prt) */
(PID.TID 0000.0001)    debLevE =  5 ; /* level of extensive debug printing */
(PID.TID 0000.0001) debugLevel =  /* select debug printing level */
(PID.TID 0000.0001)                       1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001)  plotLevel =  /* select PLOT_FIELD printing level */
(PID.TID 0000.0001)                       1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) //
(PID.TID 0000.0001) // Elliptic solver(s) paramters ( PARM02 in namelist )
(PID.TID 0000.0001) //
(PID.TID 0000.0001) cg2dMaxIters =   /* Upper limit on 2d con. grad iterations  */
(PID.TID 0000.0001)                    1000
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cg2dMinItersNSA =   /* Minimum number of iterations of 2d con. grad solver  */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cg2dUseMinResSol= /* use cg2d last-iter(=0) / min-resid.(=1) solution */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cg2dTargetResidual =   /* 2d con. grad target residual  */
(PID.TID 0000.0001)                 1.000000000000000E-19
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cg2dTargetResWunit =   /* CG2d target residual [W units] */
(PID.TID 0000.0001)                -1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cg2dPreCondFreq =   /* Freq. for updating cg2d preconditioner */
(PID.TID 0000.0001)                       1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useSRCGSolver =  /* use single reduction CG solver(s) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useNSACGSolver =  /* use not-self-adjoint CG solver */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) printResidualFreq = /* Freq. for printing CG residual */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) //
(PID.TID 0000.0001) // Time stepping paramters ( PARM03 in namelist )
(PID.TID 0000.0001) //
(PID.TID 0000.0001) deltaTMom =   /* Momentum equation timestep ( s ) */
(PID.TID 0000.0001)                 1.200000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deltaTFreeSurf = /* FreeSurface equation timestep ( s ) */
(PID.TID 0000.0001)                 1.200000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dTtracerLev =  /* Tracer equation timestep ( s ) */
(PID.TID 0000.0001)    15 @  4.320000000000000E+04              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deltaTClock  =   /* Model clock timestep ( s ) */
(PID.TID 0000.0001)                 4.320000000000000E+04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cAdjFreq =   /* Convective adjustment interval ( s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momForcingOutAB = /* =1: take Momentum Forcing out of Adams-Bash. stepping */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tracForcingOutAB = /* =1: take T,S,pTr Forcing out of Adams-Bash. stepping */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momDissip_In_AB = /* put Dissipation Tendency in Adams-Bash. stepping */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) doAB_onGtGs = /* apply AB on Tendencies (rather than on T,S)*/
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) abEps =   /* Adams-Bashforth-2 stabilizing weight */
(PID.TID 0000.0001)                 1.000000000000000E-01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) applyExchUV_early = /* Apply EXCH to U,V earlier in time-step */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tauCD =   /* CD coupling time-scale ( s ) */
(PID.TID 0000.0001)                 3.214280000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rCD =   /* Normalised CD coupling parameter */
(PID.TID 0000.0001)                 9.962666600296178E-01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) epsAB_CD = /* AB-2 stabilizing weight for CD-scheme*/
(PID.TID 0000.0001)                 1.000000000000000E-01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) pickupStrictlyMatch= /* stop if pickup do not strictly match */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nIter0   =   /* Run starting timestep number */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nTimeSteps = /* Number of timesteps */
(PID.TID 0000.0001)                       4
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nEndIter =   /* Run ending timestep number */
(PID.TID 0000.0001)                       4
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) baseTime =   /* Model base time ( s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) startTime =  /* Run start time ( s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) endTime  =   /* Integration ending time ( s ) */
(PID.TID 0000.0001)                 1.728000000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) pChkPtFreq = /* Permanent restart/pickup file interval ( s ) */
(PID.TID 0000.0001)                 3.110400000000000E+08
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) chkPtFreq  = /* Rolling restart/pickup file interval ( s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) pickup_write_mdsio =   /* Model IO flag. */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) pickup_read_mdsio =   /* Model IO flag. */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) writePickupAtEnd =   /* Model IO flag. */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dumpFreq =   /* Model state write out interval ( s ). */
(PID.TID 0000.0001)                 2.592000000000000E+06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dumpInitAndLast= /* write out Initial & Last iter. model state */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) snapshot_mdsio =   /* Model IO flag. */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) monitorFreq =   /* Monitor output interval ( s ). */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) monitorSelect = /* select group of variables to monitor */
(PID.TID 0000.0001)                       3
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) monitor_stdio =   /* Model IO flag. */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) externForcingPeriod =   /* forcing period (s) */
(PID.TID 0000.0001)                 2.592000000000000E+06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) externForcingCycle =   /* period of the cyle (s). */
(PID.TID 0000.0001)                 3.110400000000000E+07
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tauThetaClimRelax =   /* relaxation time scale (s) */
(PID.TID 0000.0001)                 5.184000000000000E+06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tauSaltClimRelax =   /* relaxation time scale (s) */
(PID.TID 0000.0001)                 1.555200000000000E+07
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) latBandClimRelax =   /* max. Lat. where relaxation */
(PID.TID 0000.0001)                 1.800000000000000E+02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) //
(PID.TID 0000.0001) // Gridding paramters ( PARM04 in namelist )
(PID.TID 0000.0001) //
(PID.TID 0000.0001) usingCartesianGrid = /* Cartesian coordinates flag ( True/False ) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingCylindricalGrid = /* Cylindrical coordinates flag ( True/False ) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingSphericalPolarGrid = /* Spherical coordinates flag ( True/False ) */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingCurvilinearGrid = /* Curvilinear coordinates flag ( True/False ) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useMin4hFacEdges = /* set hFacW,S as minimum of adjacent hFacC factor */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) interViscAr_pCell = /* account for partial-cell in interior vert. viscosity */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) interDiffKr_pCell = /* account for partial-cell in interior vert. diffusion */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) pCellMix_select = /* option to enhance mixing near surface & bottom */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectSigmaCoord = /* Hybrid-Sigma Vert. Coordinate option */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rSigmaBnd = /* r/sigma transition ( units of r ==  m ) */
(PID.TID 0000.0001)                 1.234567000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rkSign =   /* index orientation relative to vertical coordinate */
(PID.TID 0000.0001)                -1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) gravitySign = /* gravity orientation relative to vertical coordinate */
(PID.TID 0000.0001)                -1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) seaLev_Z =  /* reference height of sea-level [m] */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) top_Pres =  /* reference pressure at the top [Pa] */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) mass2rUnit = /* convert mass per unit area [kg/m2] to r-units [m] */
(PID.TID 0000.0001)                 9.661835748792270E-04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rUnit2mass = /* convert r-units [m] to mass per unit area [kg/m2] */
(PID.TID 0000.0001)                 1.035000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) drC =   /* C spacing ( units of r ) */
(PID.TID 0000.0001)                 2.500000000000000E+01,      /* K =  1 */
(PID.TID 0000.0001)                 6.000000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                 8.500000000000000E+01,      /* K =  3 */
(PID.TID 0000.0001)                 1.200000000000000E+02,      /* K =  4 */
(PID.TID 0000.0001)                 1.650000000000000E+02,      /* K =  5 */
(PID.TID 0000.0001)                 2.150000000000000E+02,      /* K =  6 */
(PID.TID 0000.0001)                 2.650000000000000E+02,      /* K =  7 */
(PID.TID 0000.0001)                 3.150000000000000E+02,      /* K =  8 */
(PID.TID 0000.0001)                 3.650000000000000E+02,      /* K =  9 */
(PID.TID 0000.0001)                 4.150000000000000E+02,      /* K = 10 */
(PID.TID 0000.0001)                 4.650000000000000E+02,      /* K = 11 */
(PID.TID 0000.0001)                 5.150000000000000E+02,      /* K = 12 */
(PID.TID 0000.0001)                 5.650000000000000E+02,      /* K = 13 */
(PID.TID 0000.0001)                 6.150000000000000E+02,      /* K = 14 */
(PID.TID 0000.0001)                 6.650000000000000E+02,      /* K = 15 */
(PID.TID 0000.0001)                 3.450000000000000E+02       /* K = 16 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) drF =   /* W spacing ( units of r ) */
(PID.TID 0000.0001)                 5.000000000000000E+01,      /* K =  1 */
(PID.TID 0000.0001)                 7.000000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                 1.000000000000000E+02,      /* K =  3 */
(PID.TID 0000.0001)                 1.400000000000000E+02,      /* K =  4 */
(PID.TID 0000.0001)                 1.900000000000000E+02,      /* K =  5 */
(PID.TID 0000.0001)                 2.400000000000000E+02,      /* K =  6 */
(PID.TID 0000.0001)                 2.900000000000000E+02,      /* K =  7 */
(PID.TID 0000.0001)                 3.400000000000000E+02,      /* K =  8 */
(PID.TID 0000.0001)                 3.900000000000000E+02,      /* K =  9 */
(PID.TID 0000.0001)                 4.400000000000000E+02,      /* K = 10 */
(PID.TID 0000.0001)                 4.900000000000000E+02,      /* K = 11 */
(PID.TID 0000.0001)                 5.400000000000000E+02,      /* K = 12 */
(PID.TID 0000.0001)                 5.900000000000000E+02,      /* K = 13 */
(PID.TID 0000.0001)                 6.400000000000000E+02,      /* K = 14 */
(PID.TID 0000.0001)                 6.900000000000000E+02       /* K = 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) delX = /* U spacing ( m - cartesian, degrees - spherical ) */
(PID.TID 0000.0001)    90 @  4.000000000000000E+00              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) delY = /* V spacing ( m - cartesian, degrees - spherical ) */
(PID.TID 0000.0001)    40 @  4.000000000000000E+00              /* J =  1: 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) xgOrigin = /* X-axis origin of West  edge (cartesian: m, lat-lon: deg) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) ygOrigin = /* Y-axis origin of South edge (cartesian: m, lat-lon: deg) */
(PID.TID 0000.0001)                -8.000000000000000E+01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rSphere =  /* Radius ( ignored - cartesian, m - spherical ) */
(PID.TID 0000.0001)                 6.370000000000000E+06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deepAtmosphere = /* Deep/Shallow Atmosphere flag (True/False) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) xC =  /* xC(:,1,:,1) : P-point X coord ( deg. or m if cartesian) */
(PID.TID 0000.0001)                 2.000000000000000E+00,      /* I =  1 */
(PID.TID 0000.0001)                 6.000000000000000E+00,      /* I =  2 */
(PID.TID 0000.0001)                 1.000000000000000E+01,      /* I =  3 */
(PID.TID 0000.0001)      . . .
(PID.TID 0000.0001)                 8.200000000000000E+01,      /* I = 21 */
(PID.TID 0000.0001)                 8.600000000000000E+01,      /* I = 22 */
(PID.TID 0000.0001)                 9.000000000000000E+01,      /* I = 23 */
(PID.TID 0000.0001)                 9.400000000000000E+01,      /* I = 24 */
(PID.TID 0000.0001)                 9.800000000000000E+01,      /* I = 25 */
(PID.TID 0000.0001)                 1.020000000000000E+02,      /* I = 26 */
(PID.TID 0000.0001)      . . .
(PID.TID 0000.0001)                 1.700000000000000E+02,      /* I = 43 */
(PID.TID 0000.0001)                 1.740000000000000E+02,      /* I = 44 */
(PID.TID 0000.0001)                 1.780000000000000E+02,      /* I = 45 */
(PID.TID 0000.0001)                 1.820000000000000E+02,      /* I = 46 */
(PID.TID 0000.0001)                 1.860000000000000E+02,      /* I = 47 */
(PID.TID 0000.0001)                 1.900000000000000E+02,      /* I = 48 */
(PID.TID 0000.0001)      . . .
(PID.TID 0000.0001)                 2.580000000000000E+02,      /* I = 65 */
(PID.TID 0000.0001)                 2.620000000000000E+02,      /* I = 66 */
(PID.TID 0000.0001)                 2.660000000000000E+02,      /* I = 67 */
(PID.TID 0000.0001)                 2.700000000000000E+02,      /* I = 68 */
(PID.TID 0000.0001)                 2.740000000000000E+02,      /* I = 69 */
(PID.TID 0000.0001)                 2.780000000000000E+02,      /* I = 70 */
(PID.TID 0000.0001)      . . .
(PID.TID 0000.0001)                 3.500000000000000E+02,      /* I = 88 */
(PID.TID 0000.0001)                 3.540000000000000E+02,      /* I = 89 */
(PID.TID 0000.0001)                 3.580000000000000E+02       /* I = 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) yC =  /* yC(1,:,1,:) : P-point Y coord ( deg. or m if cartesian) */
(PID.TID 0000.0001)                -7.800000000000000E+01,      /* J =  1 */
(PID.TID 0000.0001)                -7.400000000000000E+01,      /* J =  2 */
(PID.TID 0000.0001)                -7.000000000000000E+01,      /* J =  3 */
(PID.TID 0000.0001)                -6.600000000000000E+01,      /* J =  4 */
(PID.TID 0000.0001)                -6.200000000000000E+01,      /* J =  5 */
(PID.TID 0000.0001)                -5.800000000000000E+01,      /* J =  6 */
(PID.TID 0000.0001)                -5.400000000000000E+01,      /* J =  7 */
(PID.TID 0000.0001)                -5.000000000000000E+01,      /* J =  8 */
(PID.TID 0000.0001)                -4.600000000000000E+01,      /* J =  9 */
(PID.TID 0000.0001)                -4.200000000000000E+01,      /* J = 10 */
(PID.TID 0000.0001)                -3.800000000000000E+01,      /* J = 11 */
(PID.TID 0000.0001)                -3.400000000000000E+01,      /* J = 12 */
(PID.TID 0000.0001)                -3.000000000000000E+01,      /* J = 13 */
(PID.TID 0000.0001)                -2.600000000000000E+01,      /* J = 14 */
(PID.TID 0000.0001)                -2.200000000000000E+01,      /* J = 15 */
(PID.TID 0000.0001)                -1.800000000000000E+01,      /* J = 16 */
(PID.TID 0000.0001)                -1.400000000000000E+01,      /* J = 17 */
(PID.TID 0000.0001)                -1.000000000000000E+01,      /* J = 18 */
(PID.TID 0000.0001)                -6.000000000000000E+00,      /* J = 19 */
(PID.TID 0000.0001)                -2.000000000000000E+00,      /* J = 20 */
(PID.TID 0000.0001)                 2.000000000000000E+00,      /* J = 21 */
(PID.TID 0000.0001)                 6.000000000000000E+00,      /* J = 22 */
(PID.TID 0000.0001)                 1.000000000000000E+01,      /* J = 23 */
(PID.TID 0000.0001)                 1.400000000000000E+01,      /* J = 24 */
(PID.TID 0000.0001)                 1.800000000000000E+01,      /* J = 25 */
(PID.TID 0000.0001)                 2.200000000000000E+01,      /* J = 26 */
(PID.TID 0000.0001)                 2.600000000000000E+01,      /* J = 27 */
(PID.TID 0000.0001)                 3.000000000000000E+01,      /* J = 28 */
(PID.TID 0000.0001)                 3.400000000000000E+01,      /* J = 29 */
(PID.TID 0000.0001)                 3.800000000000000E+01,      /* J = 30 */
(PID.TID 0000.0001)                 4.200000000000000E+01,      /* J = 31 */
(PID.TID 0000.0001)                 4.600000000000000E+01,      /* J = 32 */
(PID.TID 0000.0001)                 5.000000000000000E+01,      /* J = 33 */
(PID.TID 0000.0001)                 5.400000000000000E+01,      /* J = 34 */
(PID.TID 0000.0001)                 5.800000000000000E+01,      /* J = 35 */
(PID.TID 0000.0001)                 6.200000000000000E+01,      /* J = 36 */
(PID.TID 0000.0001)                 6.600000000000000E+01,      /* J = 37 */
(PID.TID 0000.0001)                 7.000000000000000E+01,      /* J = 38 */
(PID.TID 0000.0001)                 7.400000000000000E+01,      /* J = 39 */
(PID.TID 0000.0001)                 7.800000000000000E+01       /* J = 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rcoord = /* P-point R coordinate (  units of r ) */
(PID.TID 0000.0001)                -2.500000000000000E+01,      /* K =  1 */
(PID.TID 0000.0001)                -8.500000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                -1.700000000000000E+02,      /* K =  3 */
(PID.TID 0000.0001)                -2.900000000000000E+02,      /* K =  4 */
(PID.TID 0000.0001)                -4.550000000000000E+02,      /* K =  5 */
(PID.TID 0000.0001)                -6.700000000000000E+02,      /* K =  6 */
(PID.TID 0000.0001)                -9.350000000000000E+02,      /* K =  7 */
(PID.TID 0000.0001)                -1.250000000000000E+03,      /* K =  8 */
(PID.TID 0000.0001)                -1.615000000000000E+03,      /* K =  9 */
(PID.TID 0000.0001)                -2.030000000000000E+03,      /* K = 10 */
(PID.TID 0000.0001)                -2.495000000000000E+03,      /* K = 11 */
(PID.TID 0000.0001)                -3.010000000000000E+03,      /* K = 12 */
(PID.TID 0000.0001)                -3.575000000000000E+03,      /* K = 13 */
(PID.TID 0000.0001)                -4.190000000000000E+03,      /* K = 14 */
(PID.TID 0000.0001)                -4.855000000000000E+03       /* K = 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rF =   /* W-Interf. R coordinate (  units of r ) */
(PID.TID 0000.0001)                 0.000000000000000E+00,      /* K =  1 */
(PID.TID 0000.0001)                -5.000000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                -1.200000000000000E+02,      /* K =  3 */
(PID.TID 0000.0001)                -2.200000000000000E+02,      /* K =  4 */
(PID.TID 0000.0001)                -3.600000000000000E+02,      /* K =  5 */
(PID.TID 0000.0001)                -5.500000000000000E+02,      /* K =  6 */
(PID.TID 0000.0001)                -7.900000000000000E+02,      /* K =  7 */
(PID.TID 0000.0001)                -1.080000000000000E+03,      /* K =  8 */
(PID.TID 0000.0001)                -1.420000000000000E+03,      /* K =  9 */
(PID.TID 0000.0001)                -1.810000000000000E+03,      /* K = 10 */
(PID.TID 0000.0001)                -2.250000000000000E+03,      /* K = 11 */
(PID.TID 0000.0001)                -2.740000000000000E+03,      /* K = 12 */
(PID.TID 0000.0001)                -3.280000000000000E+03,      /* K = 13 */
(PID.TID 0000.0001)                -3.870000000000000E+03,      /* K = 14 */
(PID.TID 0000.0001)                -4.510000000000000E+03,      /* K = 15 */
(PID.TID 0000.0001)                -5.200000000000000E+03       /* K = 16 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deepFacC = /* deep-model grid factor @ cell-Center (-) */
(PID.TID 0000.0001)    15 @  1.000000000000000E+00              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deepFacF = /* deep-model grid factor @ W-Interface (-) */
(PID.TID 0000.0001)    16 @  1.000000000000000E+00              /* K =  1: 16 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rotateGrid = /* use rotated grid ( True/False ) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) phiEuler = /* Euler angle, rotation about original z-coordinate [rad] */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) thetaEuler = /* Euler angle, rotation about new x-coordinate [rad] */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) psiEuler = /* Euler angle, rotation about new z-coordinate [rad] */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dxF =  /* dxF(:,1,:,1) ( units: m ) */
(PID.TID 0000.0001)    90 @  9.246038586187513E+04              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dxF =  /* dxF(1,:,1,:) ( units: m ) */
(PID.TID 0000.0001)                 9.246038586187513E+04,      /* J =  1 */
(PID.TID 0000.0001)                 1.225786591246834E+05,      /* J =  2 */
(PID.TID 0000.0001)                 1.520997414818001E+05,      /* J =  3 */
(PID.TID 0000.0001)                 1.808798091874300E+05,      /* J =  4 */
(PID.TID 0000.0001)                 2.087786486446736E+05,      /* J =  5 */
(PID.TID 0000.0001)                 2.356603395120765E+05,      /* J =  6 */
(PID.TID 0000.0001)                 2.613939168938713E+05,      /* J =  7 */
(PID.TID 0000.0001)                 2.858540093877838E+05,      /* J =  8 */
(PID.TID 0000.0001)                 3.089214498819034E+05,      /* J =  9 */
(PID.TID 0000.0001)                 3.304838561248741E+05,      /* J = 10 */
(PID.TID 0000.0001)                 3.504361782409254E+05,      /* J = 11 */
(PID.TID 0000.0001)                 3.686812105223082E+05,      /* J = 12 */
(PID.TID 0000.0001)                 3.851300650057323E+05,      /* J = 13 */
(PID.TID 0000.0001)                 3.997026045255871E+05,      /* J = 14 */
(PID.TID 0000.0001)                 4.123278331341585E+05,      /* J = 15 */
(PID.TID 0000.0001)                 4.229442419867491E+05,      /* J = 16 */
(PID.TID 0000.0001)                 4.315001090065870E+05,      /* J = 17 */
(PID.TID 0000.0001)                 4.379537508695838E+05,      /* J = 18 */
(PID.TID 0000.0001)                 4.422737260813012E+05,      /* J = 19 */
(PID.TID 0000.0001)     2 @  4.444389881567502E+05,             /* J = 20: 21 */
(PID.TID 0000.0001)                 4.422737260813012E+05,      /* J = 22 */
(PID.TID 0000.0001)                 4.379537508695838E+05,      /* J = 23 */
(PID.TID 0000.0001)                 4.315001090065870E+05,      /* J = 24 */
(PID.TID 0000.0001)                 4.229442419867491E+05,      /* J = 25 */
(PID.TID 0000.0001)                 4.123278331341585E+05,      /* J = 26 */
(PID.TID 0000.0001)                 3.997026045255871E+05,      /* J = 27 */
(PID.TID 0000.0001)                 3.851300650057323E+05,      /* J = 28 */
(PID.TID 0000.0001)                 3.686812105223082E+05,      /* J = 29 */
(PID.TID 0000.0001)                 3.504361782409254E+05,      /* J = 30 */
(PID.TID 0000.0001)                 3.304838561248741E+05,      /* J = 31 */
(PID.TID 0000.0001)                 3.089214498819034E+05,      /* J = 32 */
(PID.TID 0000.0001)                 2.858540093877838E+05,      /* J = 33 */
(PID.TID 0000.0001)                 2.613939168938713E+05,      /* J = 34 */
(PID.TID 0000.0001)                 2.356603395120765E+05,      /* J = 35 */
(PID.TID 0000.0001)                 2.087786486446736E+05,      /* J = 36 */
(PID.TID 0000.0001)                 1.808798091874300E+05,      /* J = 37 */
(PID.TID 0000.0001)                 1.520997414818001E+05,      /* J = 38 */
(PID.TID 0000.0001)                 1.225786591246834E+05,      /* J = 39 */
(PID.TID 0000.0001)                 9.246038586187513E+04       /* J = 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dyF =  /* dyF(:,1,:,1) ( units: m ) */
(PID.TID 0000.0001)    90 @  4.447098934081552E+05              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dyF =  /* dyF(1,:,1,:) ( units: m ) */
(PID.TID 0000.0001)    40 @  4.447098934081552E+05              /* J =  1: 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dxG =  /* dxG(:,1,:,1) ( units: m ) */
(PID.TID 0000.0001)    90 @  7.722306258078101E+04              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dxG =  /* dxG(1,:,1,:) ( units: m ) */
(PID.TID 0000.0001)                 7.722306258078101E+04,      /* J =  1 */
(PID.TID 0000.0001)                 1.075850604052271E+05,      /* J =  2 */
(PID.TID 0000.0001)                 1.374229146297914E+05,      /* J =  3 */
(PID.TID 0000.0001)                 1.665912582279823E+05,      /* J =  4 */
(PID.TID 0000.0001)                 1.949479859617814E+05,      /* J =  5 */
(PID.TID 0000.0001)                 2.223549467040777E+05,      /* J =  6 */
(PID.TID 0000.0001)                 2.486786164970726E+05,      /* J =  7 */
(PID.TID 0000.0001)                 2.737907490675810E+05,      /* J =  8 */
(PID.TID 0000.0001)                 2.975690006299821E+05,      /* J =  9 */
(PID.TID 0000.0001)                 3.198975259328452E+05,      /* J = 10 */
(PID.TID 0000.0001)                 3.406675426453503E+05,      /* J = 11 */
(PID.TID 0000.0001)                 3.597778613338690E+05,      /* J = 12 */
(PID.TID 0000.0001)                 3.771353784467131E+05,      /* J = 13 */
(PID.TID 0000.0001)                 3.926555299052806E+05,      /* J = 14 */
(PID.TID 0000.0001)                 4.062627030917454E+05,      /* J = 15 */
(PID.TID 0000.0001)                 4.178906052261313E+05,      /* J = 16 */
(PID.TID 0000.0001)                 4.274825863380723E+05,      /* J = 17 */
(PID.TID 0000.0001)                 4.349919152597734E+05,      /* J = 18 */
(PID.TID 0000.0001)                 4.403820072955634E+05,      /* J = 19 */
(PID.TID 0000.0001)                 4.436266024588540E+05,      /* J = 20 */
(PID.TID 0000.0001)                 4.447098934081552E+05,      /* J = 21 */
(PID.TID 0000.0001)                 4.436266024588540E+05,      /* J = 22 */
(PID.TID 0000.0001)                 4.403820072955634E+05,      /* J = 23 */
(PID.TID 0000.0001)                 4.349919152597734E+05,      /* J = 24 */
(PID.TID 0000.0001)                 4.274825863380723E+05,      /* J = 25 */
(PID.TID 0000.0001)                 4.178906052261313E+05,      /* J = 26 */
(PID.TID 0000.0001)                 4.062627030917454E+05,      /* J = 27 */
(PID.TID 0000.0001)                 3.926555299052806E+05,      /* J = 28 */
(PID.TID 0000.0001)                 3.771353784467131E+05,      /* J = 29 */
(PID.TID 0000.0001)                 3.597778613338690E+05,      /* J = 30 */
(PID.TID 0000.0001)                 3.406675426453503E+05,      /* J = 31 */
(PID.TID 0000.0001)                 3.198975259328452E+05,      /* J = 32 */
(PID.TID 0000.0001)                 2.975690006299821E+05,      /* J = 33 */
(PID.TID 0000.0001)                 2.737907490675810E+05,      /* J = 34 */
(PID.TID 0000.0001)                 2.486786164970726E+05,      /* J = 35 */
(PID.TID 0000.0001)                 2.223549467040777E+05,      /* J = 36 */
(PID.TID 0000.0001)                 1.949479859617814E+05,      /* J = 37 */
(PID.TID 0000.0001)                 1.665912582279823E+05,      /* J = 38 */
(PID.TID 0000.0001)                 1.374229146297914E+05,      /* J = 39 */
(PID.TID 0000.0001)                 1.075850604052271E+05       /* J = 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dyG =  /* dyG(:,1,:,1) ( units: m ) */
(PID.TID 0000.0001)    90 @  4.447098934081552E+05              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dyG =  /* dyG(1,:,1,:) ( units: m ) */
(PID.TID 0000.0001)    40 @  4.447098934081552E+05              /* J =  1: 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dxC =  /* dxC(:,1,:,1) ( units: m ) */
(PID.TID 0000.0001)    90 @  9.246038586187513E+04              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dxC =  /* dxC(1,:,1,:) ( units: m ) */
(PID.TID 0000.0001)                 9.246038586187513E+04,      /* J =  1 */
(PID.TID 0000.0001)                 1.225786591246834E+05,      /* J =  2 */
(PID.TID 0000.0001)                 1.520997414818001E+05,      /* J =  3 */
(PID.TID 0000.0001)                 1.808798091874300E+05,      /* J =  4 */
(PID.TID 0000.0001)                 2.087786486446736E+05,      /* J =  5 */
(PID.TID 0000.0001)                 2.356603395120765E+05,      /* J =  6 */
(PID.TID 0000.0001)                 2.613939168938713E+05,      /* J =  7 */
(PID.TID 0000.0001)                 2.858540093877838E+05,      /* J =  8 */
(PID.TID 0000.0001)                 3.089214498819034E+05,      /* J =  9 */
(PID.TID 0000.0001)                 3.304838561248741E+05,      /* J = 10 */
(PID.TID 0000.0001)                 3.504361782409254E+05,      /* J = 11 */
(PID.TID 0000.0001)                 3.686812105223082E+05,      /* J = 12 */
(PID.TID 0000.0001)                 3.851300650057323E+05,      /* J = 13 */
(PID.TID 0000.0001)                 3.997026045255871E+05,      /* J = 14 */
(PID.TID 0000.0001)                 4.123278331341585E+05,      /* J = 15 */
(PID.TID 0000.0001)                 4.229442419867491E+05,      /* J = 16 */
(PID.TID 0000.0001)                 4.315001090065870E+05,      /* J = 17 */
(PID.TID 0000.0001)                 4.379537508695838E+05,      /* J = 18 */
(PID.TID 0000.0001)                 4.422737260813012E+05,      /* J = 19 */
(PID.TID 0000.0001)     2 @  4.444389881567502E+05,             /* J = 20: 21 */
(PID.TID 0000.0001)                 4.422737260813012E+05,      /* J = 22 */
(PID.TID 0000.0001)                 4.379537508695838E+05,      /* J = 23 */
(PID.TID 0000.0001)                 4.315001090065870E+05,      /* J = 24 */
(PID.TID 0000.0001)                 4.229442419867491E+05,      /* J = 25 */
(PID.TID 0000.0001)                 4.123278331341585E+05,      /* J = 26 */
(PID.TID 0000.0001)                 3.997026045255871E+05,      /* J = 27 */
(PID.TID 0000.0001)                 3.851300650057323E+05,      /* J = 28 */
(PID.TID 0000.0001)                 3.686812105223082E+05,      /* J = 29 */
(PID.TID 0000.0001)                 3.504361782409254E+05,      /* J = 30 */
(PID.TID 0000.0001)                 3.304838561248741E+05,      /* J = 31 */
(PID.TID 0000.0001)                 3.089214498819034E+05,      /* J = 32 */
(PID.TID 0000.0001)                 2.858540093877838E+05,      /* J = 33 */
(PID.TID 0000.0001)                 2.613939168938713E+05,      /* J = 34 */
(PID.TID 0000.0001)                 2.356603395120765E+05,      /* J = 35 */
(PID.TID 0000.0001)                 2.087786486446736E+05,      /* J = 36 */
(PID.TID 0000.0001)                 1.808798091874300E+05,      /* J = 37 */
(PID.TID 0000.0001)                 1.520997414818001E+05,      /* J = 38 */
(PID.TID 0000.0001)                 1.225786591246834E+05,      /* J = 39 */
(PID.TID 0000.0001)                 9.246038586187513E+04       /* J = 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dyC =  /* dyC(:,1,:,1) ( units: m ) */
(PID.TID 0000.0001)    90 @  4.447098934081552E+05              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dyC =  /* dyC(1,:,1,:) ( units: m ) */
(PID.TID 0000.0001)    40 @  4.447098934081552E+05              /* J =  1: 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dxV =  /* dxV(:,1,:,1) ( units: m ) */
(PID.TID 0000.0001)    90 @  7.722306258078101E+04              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dxV =  /* dxV(1,:,1,:) ( units: m ) */
(PID.TID 0000.0001)                 7.722306258078101E+04,      /* J =  1 */
(PID.TID 0000.0001)                 1.075850604052271E+05,      /* J =  2 */
(PID.TID 0000.0001)                 1.374229146297914E+05,      /* J =  3 */
(PID.TID 0000.0001)                 1.665912582279823E+05,      /* J =  4 */
(PID.TID 0000.0001)                 1.949479859617814E+05,      /* J =  5 */
(PID.TID 0000.0001)                 2.223549467040777E+05,      /* J =  6 */
(PID.TID 0000.0001)                 2.486786164970726E+05,      /* J =  7 */
(PID.TID 0000.0001)                 2.737907490675810E+05,      /* J =  8 */
(PID.TID 0000.0001)                 2.975690006299821E+05,      /* J =  9 */
(PID.TID 0000.0001)                 3.198975259328452E+05,      /* J = 10 */
(PID.TID 0000.0001)                 3.406675426453503E+05,      /* J = 11 */
(PID.TID 0000.0001)                 3.597778613338690E+05,      /* J = 12 */
(PID.TID 0000.0001)                 3.771353784467131E+05,      /* J = 13 */
(PID.TID 0000.0001)                 3.926555299052806E+05,      /* J = 14 */
(PID.TID 0000.0001)                 4.062627030917454E+05,      /* J = 15 */
(PID.TID 0000.0001)                 4.178906052261313E+05,      /* J = 16 */
(PID.TID 0000.0001)                 4.274825863380723E+05,      /* J = 17 */
(PID.TID 0000.0001)                 4.349919152597734E+05,      /* J = 18 */
(PID.TID 0000.0001)                 4.403820072955634E+05,      /* J = 19 */
(PID.TID 0000.0001)                 4.436266024588540E+05,      /* J = 20 */
(PID.TID 0000.0001)                 4.447098934081552E+05,      /* J = 21 */
(PID.TID 0000.0001)                 4.436266024588540E+05,      /* J = 22 */
(PID.TID 0000.0001)                 4.403820072955634E+05,      /* J = 23 */
(PID.TID 0000.0001)                 4.349919152597734E+05,      /* J = 24 */
(PID.TID 0000.0001)                 4.274825863380723E+05,      /* J = 25 */
(PID.TID 0000.0001)                 4.178906052261313E+05,      /* J = 26 */
(PID.TID 0000.0001)                 4.062627030917454E+05,      /* J = 27 */
(PID.TID 0000.0001)                 3.926555299052806E+05,      /* J = 28 */
(PID.TID 0000.0001)                 3.771353784467131E+05,      /* J = 29 */
(PID.TID 0000.0001)                 3.597778613338690E+05,      /* J = 30 */
(PID.TID 0000.0001)                 3.406675426453503E+05,      /* J = 31 */
(PID.TID 0000.0001)                 3.198975259328452E+05,      /* J = 32 */
(PID.TID 0000.0001)                 2.975690006299821E+05,      /* J = 33 */
(PID.TID 0000.0001)                 2.737907490675810E+05,      /* J = 34 */
(PID.TID 0000.0001)                 2.486786164970726E+05,      /* J = 35 */
(PID.TID 0000.0001)                 2.223549467040777E+05,      /* J = 36 */
(PID.TID 0000.0001)                 1.949479859617814E+05,      /* J = 37 */
(PID.TID 0000.0001)                 1.665912582279823E+05,      /* J = 38 */
(PID.TID 0000.0001)                 1.374229146297914E+05,      /* J = 39 */
(PID.TID 0000.0001)                 1.075850604052271E+05       /* J = 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dyU =  /* dyU(:,1,:,1) ( units: m ) */
(PID.TID 0000.0001)    90 @  4.447098934081552E+05              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dyU =  /* dyU(1,:,1,:) ( units: m ) */
(PID.TID 0000.0001)    40 @  4.447098934081552E+05              /* J =  1: 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rA  =  /* rA (:,1,:,1) ( units: m^2 ) */
(PID.TID 0000.0001)    90 @  4.110969866729047E+10              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rA  =  /* rA (1,:,1,:) ( units: m^2 ) */
(PID.TID 0000.0001)                 4.110969866729047E+10,      /* J =  1 */
(PID.TID 0000.0001)                 5.450087291636665E+10,      /* J =  2 */
(PID.TID 0000.0001)                 6.762652439100235E+10,      /* J =  3 */
(PID.TID 0000.0001)                 8.042270623659998E+10,      /* J =  4 */
(PID.TID 0000.0001)                 9.282707674147523E+10,      /* J =  5 */
(PID.TID 0000.0001)                 1.047792030594103E+11,      /* J =  6 */
(PID.TID 0000.0001)                 1.162208556324093E+11,      /* J =  7 */
(PID.TID 0000.0001)                 1.270962918792468E+11,      /* J =  8 */
(PID.TID 0000.0001)                 1.373525277677229E+11,      /* J =  9 */
(PID.TID 0000.0001)                 1.469395959475426E+11,      /* J = 10 */
(PID.TID 0000.0001)                 1.558107891862220E+11,      /* J = 11 */
(PID.TID 0000.0001)                 1.639228879220326E+11,      /* J = 12 */
(PID.TID 0000.0001)                 1.712363708253573E+11,      /* J = 13 */
(PID.TID 0000.0001)                 1.777156073426415E+11,      /* J = 14 */
(PID.TID 0000.0001)                 1.833290312848623E+11,      /* J = 15 */
(PID.TID 0000.0001)                 1.880492946148330E+11,      /* J = 16 */
(PID.TID 0000.0001)                 1.918534006840893E+11,      /* J = 17 */
(PID.TID 0000.0001)                 1.947228162702492E+11,      /* J = 18 */
(PID.TID 0000.0001)                 1.966435618690094E+11,      /* J = 19 */
(PID.TID 0000.0001)     2 @  1.976062798008856E+11,             /* J = 20: 21 */
(PID.TID 0000.0001)                 1.966435618690094E+11,      /* J = 22 */
(PID.TID 0000.0001)                 1.947228162702492E+11,      /* J = 23 */
(PID.TID 0000.0001)                 1.918534006840893E+11,      /* J = 24 */
(PID.TID 0000.0001)                 1.880492946148330E+11,      /* J = 25 */
(PID.TID 0000.0001)                 1.833290312848623E+11,      /* J = 26 */
(PID.TID 0000.0001)                 1.777156073426415E+11,      /* J = 27 */
(PID.TID 0000.0001)                 1.712363708253573E+11,      /* J = 28 */
(PID.TID 0000.0001)                 1.639228879220326E+11,      /* J = 29 */
(PID.TID 0000.0001)                 1.558107891862220E+11,      /* J = 30 */
(PID.TID 0000.0001)                 1.469395959475426E+11,      /* J = 31 */
(PID.TID 0000.0001)                 1.373525277677229E+11,      /* J = 32 */
(PID.TID 0000.0001)                 1.270962918792468E+11,      /* J = 33 */
(PID.TID 0000.0001)                 1.162208556324093E+11,      /* J = 34 */
(PID.TID 0000.0001)                 1.047792030594103E+11,      /* J = 35 */
(PID.TID 0000.0001)                 9.282707674147523E+10,      /* J = 36 */
(PID.TID 0000.0001)                 8.042270623659998E+10,      /* J = 37 */
(PID.TID 0000.0001)                 6.762652439100235E+10,      /* J = 38 */
(PID.TID 0000.0001)                 5.450087291636665E+10,      /* J = 39 */
(PID.TID 0000.0001)                 4.110969866729047E+10       /* J = 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rAw =  /* rAw(:,1,:,1) ( units: m^2 ) */
(PID.TID 0000.0001)    90 @  4.110969866729047E+10              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rAw =  /* rAw(1,:,1,:) ( units: m^2 ) */
(PID.TID 0000.0001)                 4.110969866729047E+10,      /* J =  1 */
(PID.TID 0000.0001)                 5.450087291636665E+10,      /* J =  2 */
(PID.TID 0000.0001)                 6.762652439100235E+10,      /* J =  3 */
(PID.TID 0000.0001)                 8.042270623659998E+10,      /* J =  4 */
(PID.TID 0000.0001)                 9.282707674147523E+10,      /* J =  5 */
(PID.TID 0000.0001)                 1.047792030594103E+11,      /* J =  6 */
(PID.TID 0000.0001)                 1.162208556324093E+11,      /* J =  7 */
(PID.TID 0000.0001)                 1.270962918792468E+11,      /* J =  8 */
(PID.TID 0000.0001)                 1.373525277677229E+11,      /* J =  9 */
(PID.TID 0000.0001)                 1.469395959475426E+11,      /* J = 10 */
(PID.TID 0000.0001)                 1.558107891862220E+11,      /* J = 11 */
(PID.TID 0000.0001)                 1.639228879220326E+11,      /* J = 12 */
(PID.TID 0000.0001)                 1.712363708253573E+11,      /* J = 13 */
(PID.TID 0000.0001)                 1.777156073426415E+11,      /* J = 14 */
(PID.TID 0000.0001)                 1.833290312848623E+11,      /* J = 15 */
(PID.TID 0000.0001)                 1.880492946148330E+11,      /* J = 16 */
(PID.TID 0000.0001)                 1.918534006840893E+11,      /* J = 17 */
(PID.TID 0000.0001)                 1.947228162702492E+11,      /* J = 18 */
(PID.TID 0000.0001)                 1.966435618690094E+11,      /* J = 19 */
(PID.TID 0000.0001)     2 @  1.976062798008856E+11,             /* J = 20: 21 */
(PID.TID 0000.0001)                 1.966435618690094E+11,      /* J = 22 */
(PID.TID 0000.0001)                 1.947228162702492E+11,      /* J = 23 */
(PID.TID 0000.0001)                 1.918534006840893E+11,      /* J = 24 */
(PID.TID 0000.0001)                 1.880492946148330E+11,      /* J = 25 */
(PID.TID 0000.0001)                 1.833290312848623E+11,      /* J = 26 */
(PID.TID 0000.0001)                 1.777156073426415E+11,      /* J = 27 */
(PID.TID 0000.0001)                 1.712363708253573E+11,      /* J = 28 */
(PID.TID 0000.0001)                 1.639228879220326E+11,      /* J = 29 */
(PID.TID 0000.0001)                 1.558107891862220E+11,      /* J = 30 */
(PID.TID 0000.0001)                 1.469395959475426E+11,      /* J = 31 */
(PID.TID 0000.0001)                 1.373525277677229E+11,      /* J = 32 */
(PID.TID 0000.0001)                 1.270962918792468E+11,      /* J = 33 */
(PID.TID 0000.0001)                 1.162208556324093E+11,      /* J = 34 */
(PID.TID 0000.0001)                 1.047792030594103E+11,      /* J = 35 */
(PID.TID 0000.0001)                 9.282707674147523E+10,      /* J = 36 */
(PID.TID 0000.0001)                 8.042270623659998E+10,      /* J = 37 */
(PID.TID 0000.0001)                 6.762652439100235E+10,      /* J = 38 */
(PID.TID 0000.0001)                 5.450087291636665E+10,      /* J = 39 */
(PID.TID 0000.0001)                 4.110969866729047E+10       /* J = 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rAs =  /* rAs(:,1,:,1) ( units: m^2 ) */
(PID.TID 0000.0001)    90 @  3.433488626798250E+10              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rAs =  /* rAs(1,:,1,:) ( units: m^2 ) */
(PID.TID 0000.0001)                 3.433488626798250E+10,      /* J =  1 */
(PID.TID 0000.0001)                 4.783442523123625E+10,      /* J =  2 */
(PID.TID 0000.0001)                 6.110091968306414E+10,      /* J =  3 */
(PID.TID 0000.0001)                 7.406973659603818E+10,      /* J =  4 */
(PID.TID 0000.0001)                 8.667769319778079E+10,      /* J =  5 */
(PID.TID 0000.0001)                 9.886336479107463E+10,      /* J =  6 */
(PID.TID 0000.0001)                 1.105673840088173E+11,      /* J =  7 */
(PID.TID 0000.0001)                 1.217327300458638E+11,      /* J =  8 */
(PID.TID 0000.0001)                 1.323050064586578E+11,      /* J =  9 */
(PID.TID 0000.0001)                 1.422327061792377E+11,      /* J = 10 */
(PID.TID 0000.0001)                 1.514674624524945E+11,      /* J = 11 */
(PID.TID 0000.0001)                 1.599642844741385E+11,      /* J = 12 */
(PID.TID 0000.0001)                 1.676817765813788E+11,      /* J = 13 */
(PID.TID 0000.0001)                 1.745823399284268E+11,      /* J = 14 */
(PID.TID 0000.0001)                 1.806323556642996E+11,      /* J = 15 */
(PID.TID 0000.0001)                 1.858023487204767E+11,      /* J = 16 */
(PID.TID 0000.0001)                 1.900671314104743E+11,      /* J = 17 */
(PID.TID 0000.0001)                 1.934059261417216E+11,      /* J = 18 */
(PID.TID 0000.0001)                 1.958024666419019E+11,      /* J = 19 */
(PID.TID 0000.0001)                 1.972450772065981E+11,      /* J = 20 */
(PID.TID 0000.0001)                 1.977267295821495E+11,      /* J = 21 */
(PID.TID 0000.0001)                 1.972450772065981E+11,      /* J = 22 */
(PID.TID 0000.0001)                 1.958024666419019E+11,      /* J = 23 */
(PID.TID 0000.0001)                 1.934059261417216E+11,      /* J = 24 */
(PID.TID 0000.0001)                 1.900671314104743E+11,      /* J = 25 */
(PID.TID 0000.0001)                 1.858023487204767E+11,      /* J = 26 */
(PID.TID 0000.0001)                 1.806323556642996E+11,      /* J = 27 */
(PID.TID 0000.0001)                 1.745823399284268E+11,      /* J = 28 */
(PID.TID 0000.0001)                 1.676817765813788E+11,      /* J = 29 */
(PID.TID 0000.0001)                 1.599642844741385E+11,      /* J = 30 */
(PID.TID 0000.0001)                 1.514674624524945E+11,      /* J = 31 */
(PID.TID 0000.0001)                 1.422327061792377E+11,      /* J = 32 */
(PID.TID 0000.0001)                 1.323050064586578E+11,      /* J = 33 */
(PID.TID 0000.0001)                 1.217327300458638E+11,      /* J = 34 */
(PID.TID 0000.0001)                 1.105673840088173E+11,      /* J = 35 */
(PID.TID 0000.0001)                 9.886336479107463E+10,      /* J = 36 */
(PID.TID 0000.0001)                 8.667769319778079E+10,      /* J = 37 */
(PID.TID 0000.0001)                 7.406973659603818E+10,      /* J = 38 */
(PID.TID 0000.0001)                 6.110091968306414E+10,      /* J = 39 */
(PID.TID 0000.0001)                 4.783442523123625E+10       /* J = 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) globalArea = /* Integrated horizontal Area (m^2) */
(PID.TID 0000.0001)                 3.450614146649749E+14
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End of Model config. summary
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
(PID.TID 0000.0001) == Packages configuration : Check & print summary ==
(PID.TID 0000.0001) 
(PID.TID 0000.0001) GMREDI_CHECK: #define GMREDI
(PID.TID 0000.0001) GM_AdvForm =     /* if FALSE => use SkewFlux Form */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_InMomAsStress = /* if TRUE => apply as Eddy Stress */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_AdvSeparate = /* Calc Bolus & Euler Adv. separately */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_ExtraDiag =   /* Tensor Extra Diag (line 1&2) non 0 */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_isopycK =    /* Background Isopyc. Diffusivity [m^2/s] */
(PID.TID 0000.0001)                 1.000000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_skewflx*K =  /* Background GM_SkewFlx Diffusivity [m^2/s] */
(PID.TID 0000.0001)                 1.000000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_advec*K =    /* Backg. GM-Advec(=Bolus) Diffusivity [m^2/s]*/
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_Kmin_horiz = /* Minimum Horizontal Diffusivity [m^2/s] */
(PID.TID 0000.0001)                 5.000000000000000E+01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_Visbeck_alpha = /* Visbeck alpha coeff. [-] */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_Small_Number =  /* epsilon used in slope calc */
(PID.TID 0000.0001)                 1.000000000000000E-12
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_slopeSqCutoff = /* Slope^2 cut-off value */
(PID.TID 0000.0001)                 1.000000000000000E+08
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_taper_scheme =  /* Type of Tapering/Clipping scheme */
(PID.TID 0000.0001)               'dm95                                    '
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_maxSlope =  /* Maximum Slope (Tapering/Clipping) */
(PID.TID 0000.0001)                 1.000000000000000E-02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_facTrL2dz = /* Minimum Trans.Layer Thick. (factor of dz) */
(PID.TID 0000.0001)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_facTrL2ML = /* Max.Trans.Layer Thick. (factor of MxL Depth)*/
(PID.TID 0000.0001)                 5.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_maxTransLay = /* Maximum Transition Layer Thickness [m] */
(PID.TID 0000.0001)                 5.000000000000000E+02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_UseBVP = /* if TRUE => use bvp a la Ferrari et al. (2010) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_BVP_ModeNumber = /* Vertical mode number for BVP wave speed */
(PID.TID 0000.0001)                       1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_BVP_cMin = /* Minimum wave speed for BVP [m/s] */
(PID.TID 0000.0001)                 1.000000000000000E-01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_useSubMeso = /* if TRUE => use Sub-Meso param. (B.Fox-Kemper) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) subMeso_Ceff = /* efficiency coeff. of Mixed-Layer Eddies [-] */
(PID.TID 0000.0001)                 7.000000000000001E-02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) subMeso_invTau = /* inverse of Sub-Meso mixing time-scale [/s] */
(PID.TID 0000.0001)                 2.000000000000000E-06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) subMeso_LfMin = /* minimum length-scale "Lf" [m] */
(PID.TID 0000.0001)                 1.000000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) subMeso_Lmax = /* maximum grid-scale length [m] */
(PID.TID 0000.0001)                 1.100000000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_useLeithQG = /* if TRUE => add QG Leith viscosity to GMRedi tensor */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) CTRL_CHECK:  --> Starts to check CTRL set-up
(PID.TID 0000.0001) CTRL_CHECK:  <-- Ends Normally
(PID.TID 0000.0001) 
(PID.TID 0000.0001) COST_CHECK: #define ALLOW_COST
(PID.TID 0000.0001) GRDCHK_CHECK: grdchk package
(PID.TID 0000.0001) GAD_CHECK: #define ALLOW_GENERIC_ADVDIFF
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Check Model config. (CONFIG_CHECK):
(PID.TID 0000.0001) // CONFIG_CHECK : Normal End
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.15463194561016E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   8.88178419700125E-16  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -1.84297022087776E-14  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   1.33226762955019E-14  6.70003306073569E+00
(PID.TID 0000.0001) %CHECKPOINT         4 ckptA
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205003975202165D+06
(PID.TID 0000.0001)   local fc =  0.804687744759400D+06
(PID.TID 0000.0001)  global fc =  0.804687744759400D+06
 cg2d: Sum(rhs),rhsMax =   1.33226762955019E-14  6.70003306073569E+00
 cg2d: Sum(rhs),rhsMax =   1.33226762955019E-14  6.70003306073569E+00
 cg2d: Sum(rhs),rhsMax =   0.00000000000000E+00  0.00000000000000E+00
 cg2d: Sum(rhs),rhsMax =  -1.84297022087776E-14  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =  -1.84297022087776E-14  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =  -2.22044604925031E-16  7.86983807726830E-04
 cg2d: Sum(rhs),rhsMax =   8.88178419700125E-16  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =   8.88178419700125E-16  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =   1.15463194561016E-14  1.33029393421496E-03
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.15463194561016E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   1.15463194561016E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   6.21724893790088E-15  1.40191055720658E-03
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 ph-pack: packing ecco_cost
 ph-pack: packing ecco_ctrl
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Gradient-check starts (grdchk_main)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.15463194561016E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   8.88178419700125E-16  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -1.84297022087776E-14  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   1.33226762955019E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205003975202165D+06
(PID.TID 0000.0001)   local fc =  0.804687744759400D+06
(PID.TID 0000.0001)  global fc =  0.804687744759400D+06
(PID.TID 0000.0001) grdchk reference fc: fcref       =  8.04687744759400E+05
grad-res -------------------------------
 grad-res  proc    #    i    j    k   bi   bj iobc       fc ref            fc + eps           fc - eps
 grad-res  proc    #    i    j    k   bi   bj iobc      adj grad            fd grad          1 - fd/adj
 grad-res exact position met: 
 grad-res     0  526   26   17    1    2    2
(PID.TID 0000.0001) ====== Starts gradient-check number   1 (=ichknum) =======
 ph-test icomp, ncvarcomp, ichknum          526       29309           1
 ph-grd _loc: bi, bj, icomptest, ichknum            2           2           0           1
 ph-grd -->hit<--           26          17           1           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   26   17    1 ; bi,bj=   2   2 ; iobc=  1 ; rec=   1
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.13242748511766E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   3.99680288865056E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -6.10622663543836E-15  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   2.75335310107039E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205003997193242D+06
(PID.TID 0000.0001)   local fc =  0.804687766750477D+06
(PID.TID 0000.0001)  global fc =  0.804687766750477D+06
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  8.04687766750477E+05
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.13242748511766E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =  -4.21884749357559E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -3.99680288865056E-15  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   1.82076576038526E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205003953394813D+06
(PID.TID 0000.0001)   local fc =  0.804687722952049D+06
(PID.TID 0000.0001)  global fc =  0.804687722952049D+06
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  8.04687722952049E+05
grad-res -------------------------------
 grad-res     0    1   26   17    1    2    2    1   8.04687744759E+05  8.04687766750E+05  8.04687722952E+05
 grad-res     0    1    1  526    0    2    2    1   2.18992150034E+00  2.18992143637E+00  2.92152290138E-08
(PID.TID 0000.0001)  ADM  ref_cost_function      =  8.04687744759400E+05
(PID.TID 0000.0001)  ADM  adjoint_gradient       =  2.18992150034475E+00
(PID.TID 0000.0001)  ADM  finite-diff_grad       =  2.18992143636569E+00
(PID.TID 0000.0001) ====== End of gradient-check number   1 (ierr=  0) =======
(PID.TID 0000.0001) ====== Starts gradient-check number   2 (=ichknum) =======
 ph-test icomp, ncvarcomp, ichknum          527       29309           2
 ph-grd _loc: bi, bj, icomptest, ichknum            2           2         526           2
 ph-grd -->hit<--           30          17           1           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   30   17    1 ; bi,bj=   2   2 ; iobc=  1 ; rec=   1
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.15463194561016E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   6.99440505513849E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -5.88418203051333E-15  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   2.57571741713036E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205003984692549D+06
(PID.TID 0000.0001)   local fc =  0.804687754249784D+06
(PID.TID 0000.0001)  global fc =  0.804687754249784D+06
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  8.04687754249784E+05
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.16573417585641E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   3.49720252756924E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -6.99440505513849E-15  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   9.76996261670138E-15  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205003965894237D+06
(PID.TID 0000.0001)   local fc =  0.804687735451472D+06
(PID.TID 0000.0001)  global fc =  0.804687735451472D+06
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  8.04687735451472E+05
grad-res -------------------------------
 grad-res     0    2   30   17    1    2    2    1   8.04687744759E+05  8.04687754250E+05  8.04687735451E+05
 grad-res     0    2    2  527    0    2    2    1   9.39915722870E-01  9.39915596973E-01  1.33944963787E-07
(PID.TID 0000.0001)  ADM  ref_cost_function      =  8.04687744759400E+05
(PID.TID 0000.0001)  ADM  adjoint_gradient       =  9.39915722870128E-01
(PID.TID 0000.0001)  ADM  finite-diff_grad       =  9.39915596973151E-01
(PID.TID 0000.0001) ====== End of gradient-check number   2 (ierr=  0) =======
(PID.TID 0000.0001) ====== Starts gradient-check number   3 (=ichknum) =======
 ph-test icomp, ncvarcomp, ichknum          528       29309           3
 ph-grd _loc: bi, bj, icomptest, ichknum            2           2         527           3
 ph-grd -->hit<--           31          17           1           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   31   17    1 ; bi,bj=   2   2 ; iobc=  1 ; rec=   1
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.13242748511766E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   1.94289029309402E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -9.32587340685131E-15  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   1.33226762955019E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205003980663090D+06
(PID.TID 0000.0001)   local fc =  0.804687750220325D+06
(PID.TID 0000.0001)  global fc =  0.804687750220325D+06
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  8.04687750220325E+05
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.16573417585641E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -1.31006316905768E-14  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   1.95399252334028E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205003969921301D+06
(PID.TID 0000.0001)   local fc =  0.804687739478537D+06
(PID.TID 0000.0001)  global fc =  0.804687739478537D+06
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  8.04687739478537E+05
grad-res -------------------------------
 grad-res     0    3   31   17    1    2    2    1   8.04687744759E+05  8.04687750220E+05  8.04687739479E+05
 grad-res     0    3    3  528    0    2    2    1   5.37089397047E-01  5.37089438876E-01 -7.78807265345E-08
(PID.TID 0000.0001)  ADM  ref_cost_function      =  8.04687744759400E+05
(PID.TID 0000.0001)  ADM  adjoint_gradient       =  5.37089397047225E-01
(PID.TID 0000.0001)  ADM  finite-diff_grad       =  5.37089438876137E-01
(PID.TID 0000.0001) ====== End of gradient-check number   3 (ierr=  0) =======
(PID.TID 0000.0001) ====== Starts gradient-check number   4 (=ichknum) =======
 ph-test icomp, ncvarcomp, ichknum          529       29309           4
 ph-grd _loc: bi, bj, icomptest, ichknum            2           2         528           4
 ph-grd -->hit<--           32          17           1           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   32   17    1 ; bi,bj=   2   2 ; iobc=  1 ; rec=   1
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.16573417585641E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -9.88098491916389E-15  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   1.99840144432528E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205003997321381D+06
(PID.TID 0000.0001)   local fc =  0.804687766878617D+06
(PID.TID 0000.0001)  global fc =  0.804687766878617D+06
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  8.04687766878617E+05
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.16573417585641E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   9.76996261670138E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =   6.66133814775094E-16  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   2.39808173319034E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205003953270627D+06
(PID.TID 0000.0001)   local fc =  0.804687722827862D+06
(PID.TID 0000.0001)  global fc =  0.804687722827862D+06
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  8.04687722827862E+05
grad-res -------------------------------
 grad-res     0    4   32   17    1    2    2    1   8.04687744759E+05  8.04687766879E+05  8.04687722828E+05
 grad-res     0    4    4  529    0    2    2    1   2.20253900008E+00  2.20253771986E+00  5.81250154674E-07
(PID.TID 0000.0001)  ADM  ref_cost_function      =  8.04687744759400E+05
(PID.TID 0000.0001)  ADM  adjoint_gradient       =  2.20253900008494E+00
(PID.TID 0000.0001)  ADM  finite-diff_grad       =  2.20253771985881E+00
(PID.TID 0000.0001) ====== End of gradient-check number   4 (ierr=  0) =======
(PID.TID 0000.0001) ====== Starts gradient-check number   5 (=ichknum) =======
 ph-test icomp, ncvarcomp, ichknum          530       29309           5
 ph-grd _loc: bi, bj, icomptest, ichknum            2           2         529           5
 ph-grd -->hit<--           36          17           1           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   36   17    1 ; bi,bj=   2   2 ; iobc=  1 ; rec=   1
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.15463194561016E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =  -6.10622663543836E-16  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -1.25455201782643E-14  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   2.79776202205539E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205004049861352D+06
(PID.TID 0000.0001)   local fc =  0.804687819418588D+06
(PID.TID 0000.0001)  global fc =  0.804687819418588D+06
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  8.04687819418588E+05
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.15463194561016E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   1.05471187339390E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -8.88178419700125E-16  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   1.24344978758018E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089993D+06
 --> objf_test(bi,bj)        =  0.205003900726674D+06
(PID.TID 0000.0001)   local fc =  0.804687670283909D+06
(PID.TID 0000.0001)  global fc =  0.804687670283909D+06
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  8.04687670283909E+05
grad-res -------------------------------
 grad-res     0    5   36   17    1    2    2    1   8.04687744759E+05  8.04687819419E+05  8.04687670284E+05
 grad-res     0    5    5  530    0    2    2    1   7.45673438915E+00  7.45673393249E+00  6.12419567281E-08
(PID.TID 0000.0001)  ADM  ref_cost_function      =  8.04687744759400E+05
(PID.TID 0000.0001)  ADM  adjoint_gradient       =  7.45673438915081E+00
(PID.TID 0000.0001)  ADM  finite-diff_grad       =  7.45673393248580E+00
(PID.TID 0000.0001) ====== End of gradient-check number   5 (ierr=  0) =======
(PID.TID 0000.0001) ====== Starts gradient-check number   6 (=ichknum) =======
 ph-test icomp, ncvarcomp, ichknum          531       29309           6
 ph-grd _loc: bi, bj, icomptest, ichknum            2           2         530           6
 ph-grd -->hit<--           37          17           1           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   37   17    1 ; bi,bj=   2   2 ; iobc=  1 ; rec=   1
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.17683640610267E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   2.60902410786912E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -1.27675647831893E-14  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   1.33226762955019E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089995D+06
 --> objf_test(bi,bj)        =  0.205004038769940D+06
(PID.TID 0000.0001)   local fc =  0.804687808327177D+06
(PID.TID 0000.0001)  global fc =  0.804687808327177D+06
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  8.04687808327177E+05
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.15463194561016E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   6.66133814775094E-16  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -1.22124532708767E-14  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   2.57571741713036E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089991D+06
 --> objf_test(bi,bj)        =  0.205003911816880D+06
(PID.TID 0000.0001)   local fc =  0.804687681374114D+06
(PID.TID 0000.0001)  global fc =  0.804687681374114D+06
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  8.04687681374114E+05
grad-res -------------------------------
 grad-res     0    6   37   17    1    2    2    1   8.04687744759E+05  8.04687808327E+05  8.04687681374E+05
 grad-res     0    6    6  531    0    2    2    1   6.34765354129E+00  6.34765317664E+00  5.74465155445E-08
(PID.TID 0000.0001)  ADM  ref_cost_function      =  8.04687744759400E+05
(PID.TID 0000.0001)  ADM  adjoint_gradient       =  6.34765354128608E+00
(PID.TID 0000.0001)  ADM  finite-diff_grad       =  6.34765317663550E+00
(PID.TID 0000.0001) ====== End of gradient-check number   6 (ierr=  0) =======
(PID.TID 0000.0001) ====== Starts gradient-check number   7 (=ichknum) =======
 ph-test icomp, ncvarcomp, ichknum          532       29309           7
 ph-grd _loc: bi, bj, icomptest, ichknum            2           2         531           7
 ph-grd -->hit<--           38          17           1           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   38   17    1 ; bi,bj=   2   2 ; iobc=  1 ; rec=   1
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.16573417585641E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   4.27435864480685E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -8.65973959207622E-15  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   1.86517468137026E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089996D+06
 --> objf_test(bi,bj)        =  0.205004039670558D+06
(PID.TID 0000.0001)   local fc =  0.804687809227796D+06
(PID.TID 0000.0001)  global fc =  0.804687809227796D+06
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  8.04687809227796E+05
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.15463194561016E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   5.55111512312578E-16  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -2.10942374678780E-15  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   5.77315972805081E-15  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089990D+06
 --> objf_test(bi,bj)        =  0.205003910916544D+06
(PID.TID 0000.0001)   local fc =  0.804687680473777D+06
(PID.TID 0000.0001)  global fc =  0.804687680473777D+06
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  8.04687680473777E+05
grad-res -------------------------------
 grad-res     0    7   38   17    1    2    2    1   8.04687744759E+05  8.04687809228E+05  8.04687680474E+05
 grad-res     0    7    7  532    0    2    2    1   6.43770134172E+00  6.43770095776E+00  5.96433663569E-08
(PID.TID 0000.0001)  ADM  ref_cost_function      =  8.04687744759400E+05
(PID.TID 0000.0001)  ADM  adjoint_gradient       =  6.43770134172453E+00
(PID.TID 0000.0001)  ADM  finite-diff_grad       =  6.43770095775835E+00
(PID.TID 0000.0001) ====== End of gradient-check number   7 (ierr=  0) =======
(PID.TID 0000.0001) ====== Starts gradient-check number   8 (=ichknum) =======
 ph-test icomp, ncvarcomp, ichknum          533       29309           8
 ph-grd _loc: bi, bj, icomptest, ichknum            2           2         532           8
 ph-grd -->hit<--           39          17           1           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   39   17    1 ; bi,bj=   2   2 ; iobc=  1 ; rec=   1
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.16573417585641E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =   3.55271367880050E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -7.32747196252603E-15  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   2.08721928629529E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804090017D+06
 --> objf_test(bi,bj)        =  0.205004042215454D+06
(PID.TID 0000.0001)   local fc =  0.804687811772713D+06
(PID.TID 0000.0001)  global fc =  0.804687811772713D+06
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  8.04687811772713E+05
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =   1.14352971536391E-14  2.53674886388737E+00
 cg2d: Sum(rhs),rhsMax =  -1.88737914186277E-15  5.04325900480939E+00
 cg2d: Sum(rhs),rhsMax =  -2.20934381900406E-14  6.29008403763954E+00
 cg2d: Sum(rhs),rhsMax =   2.62012633811537E-14  6.70003306073569E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)        =  0.221431048007422D+06
 --> objf_test(bi,bj)        =  0.241096917459820D+06
 --> objf_test(bi,bj)        =  0.137155804089968D+06
 --> objf_test(bi,bj)        =  0.205003908381863D+06
(PID.TID 0000.0001)   local fc =  0.804687677939074D+06
(PID.TID 0000.0001)  global fc =  0.804687677939074D+06
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  8.04687677939074E+05
grad-res -------------------------------
 grad-res     0    8   39   17    1    2    2    1   8.04687744759E+05  8.04687811773E+05  8.04687677939E+05
 grad-res     0    8    8  533    0    2    2    1   6.69168724574E+00  6.69168197201E+00  7.88101598759E-07
(PID.TID 0000.0001)  ADM  ref_cost_function      =  8.04687744759400E+05
(PID.TID 0000.0001)  ADM  adjoint_gradient       =  6.69168724573575E+00
(PID.TID 0000.0001)  ADM  finite-diff_grad       =  6.69168197200634E+00
(PID.TID 0000.0001) ====== End of gradient-check number   8 (ierr=  0) =======
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Gradient check results  >>> START <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  EPS =   1.000000E-02
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output h.p:  Id Itile Jtile LAYER   bi   bj   X(Id)           X(Id)+/-EPS
(PID.TID 0000.0001) grdchk output h.c:  Id  FC                   FC1                  FC2
(PID.TID 0000.0001) grdchk output h.g:  Id     FC1-FC2/(2*EPS)      ADJ GRAD(FC)         1-FDGRD/ADGRD
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   1    26    17     1    2    2   0.000000000E+00 -1.000000000E-02
(PID.TID 0000.0001) grdchk output (c):   1  8.0468774475940E+05  8.0468776675048E+05  8.0468772295205E+05
(PID.TID 0000.0001) grdchk output (g):   1     2.1899214363657E+00  2.1899215003448E+00  2.9215229013779E-08
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   2    30    17     1    2    2   0.000000000E+00 -1.000000000E-02
(PID.TID 0000.0001) grdchk output (c):   2  8.0468774475940E+05  8.0468775424978E+05  8.0468773545147E+05
(PID.TID 0000.0001) grdchk output (g):   2     9.3991559697315E-01  9.3991572287013E-01  1.3394496378716E-07
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   3    31    17     1    2    2   0.000000000E+00 -1.000000000E-02
(PID.TID 0000.0001) grdchk output (c):   3  8.0468774475940E+05  8.0468775022033E+05  8.0468773947854E+05
(PID.TID 0000.0001) grdchk output (g):   3     5.3708943887614E-01  5.3708939704722E-01 -7.7880726534474E-08
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   4    32    17     1    2    2   0.000000000E+00 -1.000000000E-02
(PID.TID 0000.0001) grdchk output (c):   4  8.0468774475940E+05  8.0468776687862E+05  8.0468772282786E+05
(PID.TID 0000.0001) grdchk output (g):   4     2.2025377198588E+00  2.2025390000849E+00  5.8125015467425E-07
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   5    36    17     1    2    2   0.000000000E+00 -1.000000000E-02
(PID.TID 0000.0001) grdchk output (c):   5  8.0468774475940E+05  8.0468781941859E+05  8.0468767028391E+05
(PID.TID 0000.0001) grdchk output (g):   5     7.4567339324858E+00  7.4567343891508E+00  6.1241956728075E-08
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   6    37    17     1    2    2   0.000000000E+00 -1.000000000E-02
(PID.TID 0000.0001) grdchk output (c):   6  8.0468774475940E+05  8.0468780832718E+05  8.0468768137411E+05
(PID.TID 0000.0001) grdchk output (g):   6     6.3476531766355E+00  6.3476535412861E+00  5.7446515544513E-08
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   7    38    17     1    2    2   0.000000000E+00 -1.000000000E-02
(PID.TID 0000.0001) grdchk output (c):   7  8.0468774475940E+05  8.0468780922780E+05  8.0468768047378E+05
(PID.TID 0000.0001) grdchk output (g):   7     6.4377009577584E+00  6.4377013417245E+00  5.9643366356887E-08
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   8    39    17     1    2    2   0.000000000E+00 -1.000000000E-02
(PID.TID 0000.0001) grdchk output (c):   8  8.0468774475940E+05  8.0468781177271E+05  8.0468767793907E+05
(PID.TID 0000.0001) grdchk output (g):   8     6.6916819720063E+00  6.6916872457358E+00  7.8810159875875E-07
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk  summary  :  RMS of    8 ratios =  3.5256631882052E-07
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Gradient check results  >>> END <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
(PID.TID 0000.0001)   Seconds in section "ALL                    [THE_MODEL_MAIN]":
(PID.TID 0000.0001)           User time:   16.808654146967456
(PID.TID 0000.0001)         System time:  0.43568199360743165
(PID.TID 0000.0001)     Wall clock time:   17.317790031433105
(PID.TID 0000.0001)          No. starts:           1
(PID.TID 0000.0001)           No. stops:           1
(PID.TID 0000.0001)   Seconds in section "INITIALISE_FIXED       [THE_MODEL_MAIN]":
(PID.TID 0000.0001)           User time:   5.5339000653475523E-002
(PID.TID 0000.0001)         System time:   7.7159998472779989E-003
(PID.TID 0000.0001)     Wall clock time:   6.4998149871826172E-002
(PID.TID 0000.0001)          No. starts:           1
(PID.TID 0000.0001)           No. stops:           1
(PID.TID 0000.0001)   Seconds in section "FORWARD_STEP        [MAIN_DO_LOOP]":
(PID.TID 0000.0001)           User time:   12.721057251095772
(PID.TID 0000.0001)         System time:   8.7109729647636414E-003
(PID.TID 0000.0001)     Wall clock time:   12.763081073760986
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "LOAD_FIELDS_DRIVER  [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   3.6277584731578827E-002
(PID.TID 0000.0001)         System time:   1.1079423129558563E-003
(PID.TID 0000.0001)     Wall clock time:   3.7623882293701172E-002
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "EXTERNAL_FLDS_LOAD [LOAD_FLDS_DRIVER]":
(PID.TID 0000.0001)           User time:   3.7630975246429443E-002
(PID.TID 0000.0001)         System time:   1.1209733784198761E-003
(PID.TID 0000.0001)     Wall clock time:   3.9017915725708008E-002
(PID.TID 0000.0001)          No. starts:          76
(PID.TID 0000.0001)           No. stops:          76
(PID.TID 0000.0001)   Seconds in section "CTRL_MAP_FORCING  [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   4.0344148874282837E-004
(PID.TID 0000.0001)         System time:   7.9870223999023438E-006
(PID.TID 0000.0001)     Wall clock time:   3.8766860961914062E-004
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "DO_ATMOSPHERIC_PHYS [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   1.6113817691802979E-002
(PID.TID 0000.0001)         System time:   3.0805170536041260E-004
(PID.TID 0000.0001)     Wall clock time:   1.6480684280395508E-002
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "DO_OCEANIC_PHYS     [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   2.6618407368659973
(PID.TID 0000.0001)         System time:   7.4796751141548157E-004
(PID.TID 0000.0001)     Wall clock time:   2.6698219776153564
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "THERMODYNAMICS      [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   3.8026328682899475
(PID.TID 0000.0001)         System time:   2.1102651953697205E-004
(PID.TID 0000.0001)     Wall clock time:   3.8130850791931152
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "DYNAMICS            [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   4.7135322839021683
(PID.TID 0000.0001)         System time:   1.0105967521667480E-004
(PID.TID 0000.0001)     Wall clock time:   4.7255547046661377
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "SOLVE_FOR_PRESSURE  [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.80988600850105286
(PID.TID 0000.0001)         System time:   1.0820142924785614E-003
(PID.TID 0000.0001)     Wall clock time:  0.81301522254943848
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "MOM_CORRECTION_STEP [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.15380933880805969
(PID.TID 0000.0001)         System time:   3.1944364309310913E-005
(PID.TID 0000.0001)     Wall clock time:  0.15430212020874023
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "INTEGR_CONTINUITY   [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.15073117613792419
(PID.TID 0000.0001)         System time:   2.8993934392929077E-005
(PID.TID 0000.0001)     Wall clock time:  0.15116453170776367
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "TRC_CORRECTION_STEP [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   3.9342045783996582E-004
(PID.TID 0000.0001)         System time:   2.9951333999633789E-006
(PID.TID 0000.0001)     Wall clock time:   4.0388107299804688E-004
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "BLOCKING_EXCHANGES  [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.12369385361671448
(PID.TID 0000.0001)         System time:   2.3990869522094727E-005
(PID.TID 0000.0001)     Wall clock time:  0.12417626380920410
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "COST_TILE           [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.22009587287902832
(PID.TID 0000.0001)         System time:   5.0067901611328125E-006
(PID.TID 0000.0001)     Wall clock time:  0.22075891494750977
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "DO_THE_MODEL_IO     [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   8.8788270950317383E-003
(PID.TID 0000.0001)         System time:   1.9389837980270386E-003
(PID.TID 0000.0001)     Wall clock time:   1.0863542556762695E-002
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "DO_WRITE_PICKUP     [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   1.6277760267257690E-002
(PID.TID 0000.0001)         System time:   3.0409805476665497E-003
(PID.TID 0000.0001)     Wall clock time:   1.9399881362915039E-002
(PID.TID 0000.0001)          No. starts:          72
(PID.TID 0000.0001)           No. stops:          72
(PID.TID 0000.0001)   Seconds in section "CTRL_PACK           [THE_MODEL_MAIN]":
(PID.TID 0000.0001)           User time:   1.3020038604736328E-002
(PID.TID 0000.0001)         System time:   6.9737434387207031E-006
(PID.TID 0000.0001)     Wall clock time:   1.3069152832031250E-002
(PID.TID 0000.0001)          No. starts:           1
(PID.TID 0000.0001)           No. stops:           1
(PID.TID 0000.0001)   Seconds in section "CTRL_PACK     [THE_MODEL_MAIN]":
(PID.TID 0000.0001)           User time:   1.0198116302490234E-002
(PID.TID 0000.0001)         System time:   1.9150078296661377E-003
(PID.TID 0000.0001)     Wall clock time:   1.2162923812866211E-002
(PID.TID 0000.0001)          No. starts:           1
(PID.TID 0000.0001)           No. stops:           1
(PID.TID 0000.0001)   Seconds in section "GRDCHK_MAIN         [THE_MODEL_MAIN]":
(PID.TID 0000.0001)           User time:   12.532314300537109
(PID.TID 0000.0001)         System time:   3.1001001596450806E-002
(PID.TID 0000.0001)     Wall clock time:   12.603947162628174
(PID.TID 0000.0001)          No. starts:           1
(PID.TID 0000.0001)           No. stops:           1
(PID.TID 0000.0001)   Seconds in section "INITIALISE_VARIA    [THE_MAIN_LOOP]":
(PID.TID 0000.0001)           User time:  0.44882678985595703
(PID.TID 0000.0001)         System time:   1.6971945762634277E-002
(PID.TID 0000.0001)     Wall clock time:  0.46761512756347656
(PID.TID 0000.0001)          No. starts:          17
(PID.TID 0000.0001)           No. stops:          17
(PID.TID 0000.0001)   Seconds in section "MAIN LOOP           [THE_MAIN_LOOP]":
(PID.TID 0000.0001)           User time:   11.991822719573975
(PID.TID 0000.0001)         System time:   3.3449828624725342E-003
(PID.TID 0000.0001)     Wall clock time:   12.026318788528442
(PID.TID 0000.0001)          No. starts:          17
(PID.TID 0000.0001)           No. stops:          17
(PID.TID 0000.0001)   Seconds in section "MAIN_DO_LOOP        [THE_MAIN_LOOP]":
(PID.TID 0000.0001)           User time:   11.990252017974854
(PID.TID 0000.0001)         System time:   3.3420026302337646E-003
(PID.TID 0000.0001)     Wall clock time:   12.024749517440796
(PID.TID 0000.0001)          No. starts:          68
(PID.TID 0000.0001)           No. stops:          68
(PID.TID 0000.0001)   Seconds in section "COST_FINAL         [ADJOINT SPIN-DOWN]":
(PID.TID 0000.0001)           User time:   9.7799301147460938E-004
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:   9.8633766174316406E-004
(PID.TID 0000.0001)          No. starts:          17
(PID.TID 0000.0001)           No. stops:          17
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // Tile <-> Tile communication statistics
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // o Tile number: 000001
(PID.TID 0000.0001) //         No. X exchanges =              0
(PID.TID 0000.0001) //            Max. X spins =              0
(PID.TID 0000.0001) //            Min. X spins =     1000000000
(PID.TID 0000.0001) //          Total. X spins =              0
(PID.TID 0000.0001) //            Avg. X spins =       0.00E+00
(PID.TID 0000.0001) //         No. Y exchanges =              0
(PID.TID 0000.0001) //            Max. Y spins =              0
(PID.TID 0000.0001) //            Min. Y spins =     1000000000
(PID.TID 0000.0001) //          Total. Y spins =              0
(PID.TID 0000.0001) //            Avg. Y spins =       0.00E+00
(PID.TID 0000.0001) // o Tile number: 000002
(PID.TID 0000.0001) //         No. X exchanges =              0
(PID.TID 0000.0001) //            Max. X spins =              0
(PID.TID 0000.0001) //            Min. X spins =     1000000000
(PID.TID 0000.0001) //          Total. X spins =              0
(PID.TID 0000.0001) //            Avg. X spins =       0.00E+00
(PID.TID 0000.0001) //         No. Y exchanges =              0
(PID.TID 0000.0001) //            Max. Y spins =              0
(PID.TID 0000.0001) //            Min. Y spins =     1000000000
(PID.TID 0000.0001) //          Total. Y spins =              0
(PID.TID 0000.0001) //            Avg. Y spins =       0.00E+00
(PID.TID 0000.0001) // o Tile number: 000003
(PID.TID 0000.0001) //         No. X exchanges =              0
(PID.TID 0000.0001) //            Max. X spins =              0
(PID.TID 0000.0001) //            Min. X spins =     1000000000
(PID.TID 0000.0001) //          Total. X spins =              0
(PID.TID 0000.0001) //            Avg. X spins =       0.00E+00
(PID.TID 0000.0001) //         No. Y exchanges =              0
(PID.TID 0000.0001) //            Max. Y spins =              0
(PID.TID 0000.0001) //            Min. Y spins =     1000000000
(PID.TID 0000.0001) //          Total. Y spins =              0
(PID.TID 0000.0001) //            Avg. Y spins =       0.00E+00
(PID.TID 0000.0001) // o Tile number: 000004
(PID.TID 0000.0001) //         No. X exchanges =              0
(PID.TID 0000.0001) //            Max. X spins =              0
(PID.TID 0000.0001) //            Min. X spins =     1000000000
(PID.TID 0000.0001) //          Total. X spins =              0
(PID.TID 0000.0001) //            Avg. X spins =       0.00E+00
(PID.TID 0000.0001) //         No. Y exchanges =              0
(PID.TID 0000.0001) //            Max. Y spins =              0
(PID.TID 0000.0001) //            Min. Y spins =     1000000000
(PID.TID 0000.0001) //          Total. Y spins =              0
(PID.TID 0000.0001) //            Avg. Y spins =       0.00E+00
(PID.TID 0000.0001) // o Thread number: 000001
(PID.TID 0000.0001) //            No. barriers =          19798
(PID.TID 0000.0001) //      Max. barrier spins =              1
(PID.TID 0000.0001) //      Min. barrier spins =              1
(PID.TID 0000.0001) //     Total barrier spins =          19798
(PID.TID 0000.0001) //      Avg. barrier spins =       1.00E+00
PROGRAM MAIN: Execution ended Normally
