(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:  checkpoint69j
(PID.TID 0000.0001) // Build user:        jm_c
(PID.TID 0000.0001) // Build host:        baudelaire.mit.edu
(PID.TID 0000.0001) // Build date:        Thu Dec 25 07:47:57 PM EST 2025
(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
(PID.TID 0000.0001) ># namelist terminator. Other systems
(PID.TID 0000.0001) ># use a / character (as shown here).
(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 =   20 ; /* 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= 16.0 , 15.2 , 14.5 , 13.9 , 13.3 ,
(PID.TID 0000.0001) >       12.4 , 11.3 ,  9.9 ,  8.4 ,  6.7 ,
(PID.TID 0000.0001) >        5.2 ,  3.8 ,  2.9 ,  2.3 ,  1.8 ,
(PID.TID 0000.0001) >        1.5 ,  1.1 ,  0.8 ,  0.66,  0.63,
(PID.TID 0000.0001) > sRef= 34.65, 34.75, 34.82, 34.87, 34.90,
(PID.TID 0000.0001) >       34.90, 34.86, 34.78, 34.69, 34.60,
(PID.TID 0000.0001) >       34.58, 34.62, 34.68, 34.72, 34.73,
(PID.TID 0000.0001) >       34.74, 34.73, 34.73, 34.72, 34.72,
(PID.TID 0000.0001) > viscAr=1.E-3,
(PID.TID 0000.0001) > viscAh=5.E5,
(PID.TID 0000.0001) > no_slip_sides=.FALSE.,
(PID.TID 0000.0001) > no_slip_bottom=.TRUE.,
(PID.TID 0000.0001) > diffKhT=1.E3,
(PID.TID 0000.0001) >#- diffKrT unused when compiled with ALLOW_3D_DIFFKR
(PID.TID 0000.0001) >#diffKrT=3.E-5,
(PID.TID 0000.0001) > diffKhS=1.E3,
(PID.TID 0000.0001) > diffKrS=3.E-5,
(PID.TID 0000.0001) > tAlpha=2.E-4,
(PID.TID 0000.0001) > sBeta =7.4E-4,
(PID.TID 0000.0001) > eosType='LINEAR',
(PID.TID 0000.0001) > gravity=9.81,
(PID.TID 0000.0001) > implicitDiffusion=.TRUE.,
(PID.TID 0000.0001) > rigidLid=.FALSE.,
(PID.TID 0000.0001) > implicitFreeSurface=.TRUE.,
(PID.TID 0000.0001) > exactConserv=.TRUE.,
(PID.TID 0000.0001) > useCDscheme=.TRUE.,
(PID.TID 0000.0001) >#-nonlinear free surface
(PID.TID 0000.0001) > hFacInf=0.2,
(PID.TID 0000.0001) > hFacSup=2.0,
(PID.TID 0000.0001) > nonlinFreeSurf=3,
(PID.TID 0000.0001) > select_rStar=1,
(PID.TID 0000.0001) >#-I/O settings
(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) > useNSACGSolver=.FALSE.,
(PID.TID 0000.0001) > cg2dMaxIters=200,
(PID.TID 0000.0001) > cg2dTargetResidual=1.E-17,
(PID.TID 0000.0001) >#cg2dTargetResWunit=4.6E-23,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Time stepping parameters
(PID.TID 0000.0001) > &PARM03
(PID.TID 0000.0001) ># startTime         = 49766400000.,
(PID.TID 0000.0001) > startTime         = 0.,
(PID.TID 0000.0001) ># endTime           = 49767264000.,
(PID.TID 0000.0001) > endTime           = 345600.,
(PID.TID 0000.0001) > deltaTClock       = 86400.,
(PID.TID 0000.0001) > deltaTFreeSurf    = 86400.,
(PID.TID 0000.0001) > deltaTMom         = 2400.,
(PID.TID 0000.0001) > tauCD             = 321428.,
(PID.TID 0000.0001) > cAdjFreq          = -1.,
(PID.TID 0000.0001) > abEps             = 0.1,
(PID.TID 0000.0001) > pChkptFreq        = 86400000.,
(PID.TID 0000.0001) > chkptFreq         = 86400000.,
(PID.TID 0000.0001) > dumpFreq          = 86400.,
(PID.TID 0000.0001) > monitorFreq       = 43200.,
(PID.TID 0000.0001) > adjMonitorFreq    = 43200.,
(PID.TID 0000.0001) > monitorSelect     = 4,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > tauThetaClimRelax = 2592000.,
(PID.TID 0000.0001) > tauSaltClimRelax  = 2592000.,
(PID.TID 0000.0001) >#tauTr1ClimRelax   = 31104000., <= moved to data.ptracers
(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.,  50.,  55.,  60.,
(PID.TID 0000.0001) >        65.,  70.,  80.,  95.,
(PID.TID 0000.0001) >       120., 155., 200., 260.,
(PID.TID 0000.0001) >       320., 400., 480., 570.,
(PID.TID 0000.0001) >       655., 725., 775., 815.,
(PID.TID 0000.0001) > ygOrigin=-80.,
(PID.TID 0000.0001) > delY=40*4.,
(PID.TID 0000.0001) > delX=90*4.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Input datasets
(PID.TID 0000.0001) > &PARM05
(PID.TID 0000.0001) > hydrogThetaFile='theta.bin',
(PID.TID 0000.0001) > hydrogSaltFile='salt.bin',
(PID.TID 0000.0001) > bathyFile='topog.bin',
(PID.TID 0000.0001) > zonalWindFile='windx.bin',
(PID.TID 0000.0001) > meridWindFile='windy.bin',
(PID.TID 0000.0001) > thetaClimFile='SST.bin'
(PID.TID 0000.0001) > saltClimFile='SSS.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)  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) > usePtracers  = .TRUE.,
(PID.TID 0000.0001) > useGMRedi    = .TRUE.,
(PID.TID 0000.0001) > useKPP       = .FALSE.,
(PID.TID 0000.0001) > useGrdchk    = .TRUE.,
(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/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 )
 pkg/ptracers             compiled   and   used ( usePTRACERS              = 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/monitor              compiled   and   used ( monitorFreq > 0.         = 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-20,
(PID.TID 0000.0001) >  GM_slopeSqCutoff = 1.D+08,
(PID.TID 0000.0001) >  GM_AdvForm         = .TRUE.,
(PID.TID 0000.0001) >  GM_isopycK         = 1.1D+3,
(PID.TID 0000.0001) >  GM_background_K    = 0.9D+3,
(PID.TID 0000.0001) >#  GM_background_K    = 1.D+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.,
(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)  PTRACERS_READPARMS: opening data.ptracers
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.ptracers
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.ptracers"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># This set-up makes is supposed to immitate salt. That is, PTRACER(:,:,:,1)
(PID.TID 0000.0001) ># should look exactly like S(:,:,:).
(PID.TID 0000.0001) > &PTRACERS_PARM01
(PID.TID 0000.0001) > PTRACERS_numInUse=1,
(PID.TID 0000.0001) > PTRACERS_advScheme=2,
(PID.TID 0000.0001) > PTRACERS_diffKh=1.E3,
(PID.TID 0000.0001) > PTRACERS_diffKr=3.E-5,
(PID.TID 0000.0001) > PTRACERS_initialFile=' ',
(PID.TID 0000.0001) > PTRACERS_ref=1., 19*0.,
(PID.TID 0000.0001) > tauTr1ClimRelax= 31104000.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  PTRACERS_READPARMS: finished reading data.ptracers
(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) > xx_genarr3d_weight(1) = 'ones_64b.bin',
(PID.TID 0000.0001) > xx_genarr3d_file(1)='xx_ptr1',
(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) // =======================================================
(PID.TID 0000.0001) // cost configuration  >>> START <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) lastinterval =   /* cost interval over which to average ( s ). */
(PID.TID 0000.0001)                 2.592000000000000E+06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cost_mask_file = /* file name of cost mask file */
(PID.TID 0000.0001)               ''
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // cost configuration  >>> END <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
(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-4,
(PID.TID 0000.0001) >### nbeg             = 1,
(PID.TID 0000.0001) > iGloPos          = 15,
(PID.TID 0000.0001) > jGloPos          = 5,
(PID.TID 0000.0001) > kGloPos          = 4,
(PID.TID 0000.0001) > iGloTile         = 1,
(PID.TID 0000.0001) > jgloTile         = 1,
(PID.TID 0000.0001) > nstep            = 1,
(PID.TID 0000.0001) > nend             = 4,
(PID.TID 0000.0001) > grdchkvarname    = "xx_ptr1",
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001) GRDCHK_READPARMS: finished reading data.grdchk
(PID.TID 0000.0001) SET_PARMS: done
(PID.TID 0000.0001) Enter INI_VERTICAL_GRID: setInterFDr=    T ; setCenterDr=    F
(PID.TID 0000.0001) %MON XC_max                       =   3.5800000000000E+02
(PID.TID 0000.0001) %MON XC_min                       =   2.0000000000000E+00
(PID.TID 0000.0001) %MON XC_mean                      =   1.8000000000000E+02
(PID.TID 0000.0001) %MON XC_sd                        =   1.0391663325314E+02
(PID.TID 0000.0001) %MON XG_max                       =   3.5600000000000E+02
(PID.TID 0000.0001) %MON XG_min                       =   0.0000000000000E+00
(PID.TID 0000.0001) %MON XG_mean                      =   1.7800000000000E+02
(PID.TID 0000.0001) %MON XG_sd                        =   1.0391663325314E+02
(PID.TID 0000.0001) %MON DXC_max                      =   4.4443898815675E+05
(PID.TID 0000.0001) %MON DXC_min                      =   9.2460385861875E+04
(PID.TID 0000.0001) %MON DXC_mean                     =   3.1372497618153E+05
(PID.TID 0000.0001) %MON DXC_sd                       =   1.1216447457560E+05
(PID.TID 0000.0001) %MON DXF_max                      =   4.4443898815675E+05
(PID.TID 0000.0001) %MON DXF_min                      =   9.2460385861875E+04
(PID.TID 0000.0001) %MON DXF_mean                     =   3.1372497618153E+05
(PID.TID 0000.0001) %MON DXF_sd                       =   1.1216447457560E+05
(PID.TID 0000.0001) %MON DXG_max                      =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DXG_min                      =   7.7223062580781E+04
(PID.TID 0000.0001) %MON DXG_mean                     =   3.1353386340261E+05
(PID.TID 0000.0001) %MON DXG_sd                       =   1.1256651772502E+05
(PID.TID 0000.0001) %MON DXV_max                      =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DXV_min                      =   7.7223062580781E+04
(PID.TID 0000.0001) %MON DXV_mean                     =   3.1353386340261E+05
(PID.TID 0000.0001) %MON DXV_sd                       =   1.1256651772502E+05
(PID.TID 0000.0001) %MON YC_max                       =   7.8000000000000E+01
(PID.TID 0000.0001) %MON YC_min                       =  -7.8000000000000E+01
(PID.TID 0000.0001) %MON YC_mean                      =   0.0000000000000E+00
(PID.TID 0000.0001) %MON YC_sd                        =   4.6173585522461E+01
(PID.TID 0000.0001) %MON YG_max                       =   7.6000000000000E+01
(PID.TID 0000.0001) %MON YG_min                       =  -8.0000000000000E+01
(PID.TID 0000.0001) %MON YG_mean                      =  -2.0000000000000E+00
(PID.TID 0000.0001) %MON YG_sd                        =   4.6173585522461E+01
(PID.TID 0000.0001) %MON DYC_max                      =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYC_min                      =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYC_mean                     =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYC_sd                       =   5.2968971431255E-09
(PID.TID 0000.0001) %MON DYF_max                      =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYF_min                      =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYF_mean                     =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYF_sd                       =   5.2968971431255E-09
(PID.TID 0000.0001) %MON DYG_max                      =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYG_min                      =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYG_mean                     =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYG_sd                       =   5.2968971431255E-09
(PID.TID 0000.0001) %MON DYU_max                      =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYU_min                      =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYU_mean                     =   4.4470989340816E+05
(PID.TID 0000.0001) %MON DYU_sd                       =   5.2968971431255E-09
(PID.TID 0000.0001) %MON RA_max                       =   1.9760627980089E+11
(PID.TID 0000.0001) %MON RA_min                       =   4.1109698667290E+10
(PID.TID 0000.0001) %MON RA_mean                      =   1.3948826965196E+11
(PID.TID 0000.0001) %MON RA_sd                        =   4.9870522472902E+10
(PID.TID 0000.0001) %MON RAW_max                      =   1.9760627980089E+11
(PID.TID 0000.0001) %MON RAW_min                      =   4.1109698667290E+10
(PID.TID 0000.0001) %MON RAW_mean                     =   1.3948826965196E+11
(PID.TID 0000.0001) %MON RAW_sd                       =   4.9870522472902E+10
(PID.TID 0000.0001) %MON RAS_max                      =   1.9772672958215E+11
(PID.TID 0000.0001) %MON RAS_min                      =   3.4334886267983E+10
(PID.TID 0000.0001) %MON RAS_mean                     =   1.3940329716694E+11
(PID.TID 0000.0001) %MON RAS_sd                       =   5.0049278732353E+10
(PID.TID 0000.0001) %MON RAZ_max                      =   1.9772672958215E+11
(PID.TID 0000.0001) %MON RAZ_min                      =   3.4334886267983E+10
(PID.TID 0000.0001) %MON RAZ_mean                     =   1.3940329716694E+11
(PID.TID 0000.0001) %MON RAZ_sd                       =   5.0049278732353E+10
(PID.TID 0000.0001) %MON AngleCS_max                  =   1.0000000000000E+00
(PID.TID 0000.0001) %MON AngleCS_min                  =   1.0000000000000E+00
(PID.TID 0000.0001) %MON AngleCS_mean                 =   1.0000000000000E+00
(PID.TID 0000.0001) %MON AngleCS_sd                   =   0.0000000000000E+00
(PID.TID 0000.0001) %MON AngleSN_max                  =   0.0000000000000E+00
(PID.TID 0000.0001) %MON AngleSN_min                  =   0.0000000000000E+00
(PID.TID 0000.0001) %MON AngleSN_mean                 =   0.0000000000000E+00
(PID.TID 0000.0001) %MON AngleSN_sd                   =   0.0000000000000E+00
(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_INIT_FIXED: ivar=   1 = number of CTRL variables defined
(PID.TID 0000.0001) 
(PID.TID 0000.0001) ctrl-wet 1:    nvarlength =        37994
(PID.TID 0000.0001) ctrl-wet 2: surface wet C =          616
(PID.TID 0000.0001) ctrl-wet 3: surface wet W =          586
(PID.TID 0000.0001) ctrl-wet 4: surface wet S =          557
(PID.TID 0000.0001) ctrl-wet 5: 3D wet points =        10880
(PID.TID 0000.0001) ctrl-wet 6: no recs for ivar =     1           1
(PID.TID 0000.0001) ctrl-wet 6: no recs for ivar =     2           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for ivar =     3           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for ivar =     4           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for ivar =     5           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for ivar =     6           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for ivar =     7           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for ivar =     8           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for ivar =     9           0
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl-wet 13: global nvarlength for Nr =   20       37994
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=    1        2158        2004        2052
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=    2        2155        2003        2049
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=    3        2154        2002        2048
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=    4        2151        2000        2047
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=    5        2148        1997        2044
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=    6        2143        1993        2040
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=    7        2143        1993        2040
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=    8        2141        1993        2038
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=    9        2139        1991        2036
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=   10        2134        1986        2030
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=   11        2126        1977        2022
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=   12        2108        1961        2003
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=   13        2101        1954        1995
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=   14        2086        1938        1981
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=   15        2054        1903        1946
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=   16        1998        1841        1879
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=   17        1844        1650        1702
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=   18        1180         951         996
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=   19         754         556         602
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W k=   20         277         168         186
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl_init_wet: no. of control variables:            1
(PID.TID 0000.0001) ctrl_init_wet: control vector length:           37994
(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 =    18000
(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): 001 001   10880    9660   10196
(PID.TID 0000.0001)  bi,bj,#(c/s/w): 002 001   12622   11628   12169
(PID.TID 0000.0001)  bi,bj,#(c/s/w): 001 002    4933    4598    4381
(PID.TID 0000.0001)  bi,bj,#(c/s/w): 002 002    9559    8975    8990
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  -> 3d control, genarr3d no.  1 is in use
(PID.TID 0000.0001)       file       = xx_ptr1
(PID.TID 0000.0001)       ncvartype  = Arr3D
(PID.TID 0000.0001)       index      =     1  (use this for pkg/grdchk)
(PID.TID 0000.0001)       ncvarindex =     1
(PID.TID 0000.0001)       weight     = ones_64b.bin
(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) %MON fCori_max                    =   1.4265546244797E-04
(PID.TID 0000.0001) %MON fCori_min                    =  -1.4265546244797E-04
(PID.TID 0000.0001) %MON fCori_mean                   =   4.7801269115805E-19
(PID.TID 0000.0001) %MON fCori_sd                     =   9.6599226301546E-05
(PID.TID 0000.0001) %MON fCoriG_max                   =   1.4151032568025E-04
(PID.TID 0000.0001) %MON fCoriG_min                   =  -1.4362679550910E-04
(PID.TID 0000.0001) %MON fCoriG_mean                  =  -3.5906698877273E-06
(PID.TID 0000.0001) %MON fCoriG_sd                    =   9.6548915696246E-05
(PID.TID 0000.0001) %MON fCoriCos_max                 =   1.4575362704741E-04
(PID.TID 0000.0001) %MON fCoriCos_min                 =   3.0322354601388E-05
(PID.TID 0000.0001) %MON fCoriCos_mean                =   1.0288600773633E-04
(PID.TID 0000.0001) %MON fCoriCos_sd                  =   3.6784304327266E-05
(PID.TID 0000.0001) INI_CG2D: CG2D normalisation factor =  8.8319877372200804E-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)                 1.600000000000000E+01,      /* K =  1 */
(PID.TID 0000.0001)                 1.520000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                 1.450000000000000E+01,      /* K =  3 */
(PID.TID 0000.0001)                 1.390000000000000E+01,      /* K =  4 */
(PID.TID 0000.0001)                 1.330000000000000E+01,      /* K =  5 */
(PID.TID 0000.0001)                 1.240000000000000E+01,      /* K =  6 */
(PID.TID 0000.0001)                 1.130000000000000E+01,      /* K =  7 */
(PID.TID 0000.0001)                 9.900000000000000E+00,      /* K =  8 */
(PID.TID 0000.0001)                 8.400000000000000E+00,      /* K =  9 */
(PID.TID 0000.0001)                 6.700000000000000E+00,      /* K = 10 */
(PID.TID 0000.0001)                 5.200000000000000E+00,      /* K = 11 */
(PID.TID 0000.0001)                 3.800000000000000E+00,      /* K = 12 */
(PID.TID 0000.0001)                 2.900000000000000E+00,      /* K = 13 */
(PID.TID 0000.0001)                 2.300000000000000E+00,      /* K = 14 */
(PID.TID 0000.0001)                 1.800000000000000E+00,      /* K = 15 */
(PID.TID 0000.0001)                 1.500000000000000E+00,      /* K = 16 */
(PID.TID 0000.0001)                 1.100000000000000E+00,      /* K = 17 */
(PID.TID 0000.0001)                 8.000000000000000E-01,      /* K = 18 */
(PID.TID 0000.0001)                 6.600000000000000E-01,      /* K = 19 */
(PID.TID 0000.0001)                 6.300000000000000E-01       /* K = 20 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) sRef =   /* Reference salinity profile ( g/kg ) */
(PID.TID 0000.0001)                 3.465000000000000E+01,      /* K =  1 */
(PID.TID 0000.0001)                 3.475000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                 3.482000000000000E+01,      /* K =  3 */
(PID.TID 0000.0001)                 3.487000000000000E+01,      /* K =  4 */
(PID.TID 0000.0001)     2 @  3.490000000000000E+01,             /* K =  5:  6 */
(PID.TID 0000.0001)                 3.486000000000000E+01,      /* K =  7 */
(PID.TID 0000.0001)                 3.478000000000000E+01,      /* K =  8 */
(PID.TID 0000.0001)                 3.469000000000000E+01,      /* K =  9 */
(PID.TID 0000.0001)                 3.460000000000000E+01,      /* K = 10 */
(PID.TID 0000.0001)                 3.458000000000000E+01,      /* K = 11 */
(PID.TID 0000.0001)                 3.462000000000000E+01,      /* K = 12 */
(PID.TID 0000.0001)                 3.468000000000000E+01,      /* K = 13 */
(PID.TID 0000.0001)                 3.472000000000000E+01,      /* K = 14 */
(PID.TID 0000.0001)                 3.473000000000000E+01,      /* K = 15 */
(PID.TID 0000.0001)                 3.474000000000000E+01,      /* K = 16 */
(PID.TID 0000.0001)     2 @  3.473000000000000E+01,             /* K = 17: 18 */
(PID.TID 0000.0001)     2 @  3.472000000000000E+01              /* K = 19: 20 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoRef =   /* Density vertical profile from (Ref,sRef)( kg/m^3 ) */
(PID.TID 0000.0001)                 9.998000000000000E+02,      /* K =  1 */
(PID.TID 0000.0001)                 1.000033953200000E+03,      /* K =  2 */
(PID.TID 0000.0001)                 1.000225714840000E+03,      /* K =  3 */
(PID.TID 0000.0001)                 1.000382683440000E+03,      /* K =  4 */
(PID.TID 0000.0001)                 1.000524855000000E+03,      /* K =  5 */
(PID.TID 0000.0001)                 1.000704819000000E+03,      /* K =  6 */
(PID.TID 0000.0001)                 1.000895180920000E+03,      /* K =  7 */
(PID.TID 0000.0001)                 1.001115936760000E+03,      /* K =  8 */
(PID.TID 0000.0001)                 1.001349290080000E+03,      /* K =  9 */
(PID.TID 0000.0001)                 1.001622635400000E+03,      /* K = 10 */
(PID.TID 0000.0001)                 1.001907778360000E+03,      /* K = 11 */
(PID.TID 0000.0001)                 1.002217316440000E+03,      /* K = 12 */
(PID.TID 0000.0001)                 1.002441671560000E+03,      /* K = 13 */
(PID.TID 0000.0001)                 1.002591241640000E+03,      /* K = 14 */
(PID.TID 0000.0001)                 1.002698620160000E+03,      /* K = 15 */
(PID.TID 0000.0001)                 1.002766006680000E+03,      /* K = 16 */
(PID.TID 0000.0001)                 1.002838592160000E+03,      /* K = 17 */
(PID.TID 0000.0001)                 1.002898580160000E+03,      /* K = 18 */
(PID.TID 0000.0001)                 1.002919176040000E+03,      /* K = 19 */
(PID.TID 0000.0001)                 1.002925174840000E+03       /* K = 20 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dBdrRef = /* Vertical grad. of reference buoyancy [(m/s/r)^2] */
(PID.TID 0000.0001)                 0.000000000000000E+00,      /* K =  1 */
(PID.TID 0000.0001)                 4.591080000000024E-05,      /* K =  2 */
(PID.TID 0000.0001)                 3.583920000000001E-05,      /* K =  3 */
(PID.TID 0000.0001)                 2.678556521739093E-05,      /* K =  4 */
(PID.TID 0000.0001)                 2.231971200000012E-05,      /* K =  5 */
(PID.TID 0000.0001)                 2.616000000000001E-05,      /* K =  6 */
(PID.TID 0000.0001)                 2.490432000000007E-05,      /* K =  7 */
(PID.TID 0000.0001)                 2.475483428571444E-05,      /* K =  8 */
(PID.TID 0000.0001)                 2.129910697674396E-05,      /* K =  9 */
(PID.TID 0000.0001)                 1.950584727272747E-05,      /* K = 10 */
(PID.TID 0000.0001)                 1.576232112676044E-05,      /* K = 11 */
(PID.TID 0000.0001)                 1.320511304347824E-05,      /* K = 12 */
(PID.TID 0000.0001)                 7.590910344827643E-06,      /* K = 13 */
(PID.TID 0000.0001)                 4.076599999999984E-06,      /* K = 14 */
(PID.TID 0000.0001)                 2.394531818181785E-06,      /* K = 15 */
(PID.TID 0000.0001)                 1.259417142857214E-06,      /* K = 16 */
(PID.TID 0000.0001)                 1.162785306122388E-06,      /* K = 17 */
(PID.TID 0000.0001)                 8.530434782608699E-07,      /* K = 18 */
(PID.TID 0000.0001)                 2.694480000000193E-07,      /* K = 19 */
(PID.TID 0000.0001)                 7.403773584905668E-08       /* K = 20 */
(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)                   F
(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)    20 @  1.000000000000000E-03              /* K =  1: 20 */
(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)                 1.000000000000000E+03
(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)                 1.000000000000000E+03
(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)    20 @  0.000000000000000E+00              /* K =  1: 20 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrNrS = /* vertical profile of vertical diffusion of Salt ( m^2/s )*/
(PID.TID 0000.0001)    20 @  3.000000000000000E-05              /* K =  1: 20 */
(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)                 0.000000000000000E+00
(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)               'LINEAR'
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tAlpha = /* Linear EOS thermal expansion coefficient ( 1/oC ) */
(PID.TID 0000.0001)                 2.000000000000000E-04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) sBeta  = /* Linear EOS haline contraction coefficient ( 1/(g/kg) ) */
(PID.TID 0000.0001)                 7.400000000000000E-04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoNil    = /* Reference density for Linear EOS ( kg/m^3 ) */
(PID.TID 0000.0001)                 9.998000000000000E+02
(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)                 9.998000000000000E+02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoFacC = /* normalized Reference density @ cell-Center (-) */
(PID.TID 0000.0001)    20 @  1.000000000000000E+00              /* K =  1: 20 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoFacF = /* normalized Reference density @ W-Interface (-) */
(PID.TID 0000.0001)    21 @  1.000000000000000E+00              /* K =  1: 21 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoConstFresh = /* Fresh-water reference density ( kg/m^3 ) */
(PID.TID 0000.0001)                 9.998000000000000E+02
(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)    20 @  1.000000000000000E+00              /* K =  1: 20 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) gravFacF = /* gravity factor (vs surf.) @ W-Interface (-) */
(PID.TID 0000.0001)    21 @  1.000000000000000E+00              /* K =  1: 21 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rotationPeriod =   /* Rotation Period ( s ) */
(PID.TID 0000.0001)                 8.616400000000000E+04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) omega =   /* Angular velocity ( rad/s ) */
(PID.TID 0000.0001)                 7.292123516990375E-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)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hFacMinDr = /* minimum partial cell thickness ( m) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) exactConserv = /* Update etaN from continuity Eq on/off flag */
(PID.TID 0000.0001)                   T
(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)                       3
(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)                       1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useRealFreshWaterFlux = /* Real Fresh Water Flux on/off flag*/
(PID.TID 0000.0001)                   F
(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)                 3.500000000000000E+01
(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) selectMetricTerms= /* Scheme selector for Horizontal Metric Terms */
(PID.TID 0000.0001)                       1
(PID.TID 0000.0001)    = 0 : Off (ignore Spherical/Cylindrical Metric Terms)
(PID.TID 0000.0001)    = 1 : original discretization
(PID.TID 0000.0001)    = 2 : using (Spherical) grid-spacing
(PID.TID 0000.0001)    = 3 : as 2 but gU-Metric inside Advection
(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) 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) 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) select3dCoriScheme= /* Scheme selector for 3-D Coriolis-Term */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)    = 0 : Off (ignore 3-D Coriolis Terms in Omega.Cos(Lat) )
(PID.TID 0000.0001)    = 1 : original discretization ; = 2 : using averaged Transport
(PID.TID 0000.0001)    = 3 : same as 2 with hFac in gW_Cor
(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)    = 4 : hFac weighted average (Angular Mom. conserving)
(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)                   T
(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) selectPenetratingSW = /* short wave penetration selector */
(PID.TID 0000.0001)                       1
(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)                     200
(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-17
(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)                      -1
(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)                 2.400000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deltaTFreeSurf = /* FreeSurface equation timestep ( s ) */
(PID.TID 0000.0001)                 8.640000000000000E+04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dTtracerLev =  /* Tracer equation timestep ( s ) */
(PID.TID 0000.0001)    20 @  8.640000000000000E+04              /* K =  1: 20 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deltaTClock  =   /* Model clock timestep ( s ) */
(PID.TID 0000.0001)                 8.640000000000000E+04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cAdjFreq =   /* Convective adjustment interval ( s ) */
(PID.TID 0000.0001)                 8.640000000000000E+04
(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.925333200592357E-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)                 3.456000000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) pChkPtFreq = /* Permanent restart/pickup file interval ( s ) */
(PID.TID 0000.0001)                 8.640000000000000E+07
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) chkPtFreq  = /* Rolling restart/pickup file interval ( s ) */
(PID.TID 0000.0001)                 8.640000000000000E+07
(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)                 8.640000000000000E+04
(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)                 4.320000000000000E+04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) monitorSelect = /* select group of variables to monitor */
(PID.TID 0000.0001)                       4
(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)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) externForcingCycle =   /* period of the cyle (s). */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tauThetaClimRelax =   /* relaxation time scale (s) */
(PID.TID 0000.0001)                 2.592000000000000E+06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tauSaltClimRelax =   /* relaxation time scale (s) */
(PID.TID 0000.0001)                 2.592000000000000E+06
(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)                 1.000200040008002E-03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rUnit2mass = /* convert r-units [m] to mass per unit area [kg/m2] */
(PID.TID 0000.0001)                 9.998000000000000E+02
(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)                 5.000000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                 5.250000000000000E+01,      /* K =  3 */
(PID.TID 0000.0001)                 5.750000000000000E+01,      /* K =  4 */
(PID.TID 0000.0001)                 6.250000000000000E+01,      /* K =  5 */
(PID.TID 0000.0001)                 6.750000000000000E+01,      /* K =  6 */
(PID.TID 0000.0001)                 7.500000000000000E+01,      /* K =  7 */
(PID.TID 0000.0001)                 8.750000000000000E+01,      /* K =  8 */
(PID.TID 0000.0001)                 1.075000000000000E+02,      /* K =  9 */
(PID.TID 0000.0001)                 1.375000000000000E+02,      /* K = 10 */
(PID.TID 0000.0001)                 1.775000000000000E+02,      /* K = 11 */
(PID.TID 0000.0001)                 2.300000000000000E+02,      /* K = 12 */
(PID.TID 0000.0001)                 2.900000000000000E+02,      /* K = 13 */
(PID.TID 0000.0001)                 3.600000000000000E+02,      /* K = 14 */
(PID.TID 0000.0001)                 4.400000000000000E+02,      /* K = 15 */
(PID.TID 0000.0001)                 5.250000000000000E+02,      /* K = 16 */
(PID.TID 0000.0001)                 6.125000000000000E+02,      /* K = 17 */
(PID.TID 0000.0001)                 6.900000000000000E+02,      /* K = 18 */
(PID.TID 0000.0001)                 7.500000000000000E+02,      /* K = 19 */
(PID.TID 0000.0001)                 7.950000000000000E+02,      /* K = 20 */
(PID.TID 0000.0001)                 4.075000000000000E+02       /* K = 21 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) drF =   /* W spacing ( units of r ) */
(PID.TID 0000.0001)     2 @  5.000000000000000E+01,             /* K =  1:  2 */
(PID.TID 0000.0001)                 5.500000000000000E+01,      /* K =  3 */
(PID.TID 0000.0001)                 6.000000000000000E+01,      /* K =  4 */
(PID.TID 0000.0001)                 6.500000000000000E+01,      /* K =  5 */
(PID.TID 0000.0001)                 7.000000000000000E+01,      /* K =  6 */
(PID.TID 0000.0001)                 8.000000000000000E+01,      /* K =  7 */
(PID.TID 0000.0001)                 9.500000000000000E+01,      /* K =  8 */
(PID.TID 0000.0001)                 1.200000000000000E+02,      /* K =  9 */
(PID.TID 0000.0001)                 1.550000000000000E+02,      /* K = 10 */
(PID.TID 0000.0001)                 2.000000000000000E+02,      /* K = 11 */
(PID.TID 0000.0001)                 2.600000000000000E+02,      /* K = 12 */
(PID.TID 0000.0001)                 3.200000000000000E+02,      /* K = 13 */
(PID.TID 0000.0001)                 4.000000000000000E+02,      /* K = 14 */
(PID.TID 0000.0001)                 4.800000000000000E+02,      /* K = 15 */
(PID.TID 0000.0001)                 5.700000000000000E+02,      /* K = 16 */
(PID.TID 0000.0001)                 6.550000000000000E+02,      /* K = 17 */
(PID.TID 0000.0001)                 7.250000000000000E+02,      /* K = 18 */
(PID.TID 0000.0001)                 7.750000000000000E+02,      /* K = 19 */
(PID.TID 0000.0001)                 8.150000000000000E+02       /* K = 20 */
(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)                -7.500000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                -1.275000000000000E+02,      /* K =  3 */
(PID.TID 0000.0001)                -1.850000000000000E+02,      /* K =  4 */
(PID.TID 0000.0001)                -2.475000000000000E+02,      /* K =  5 */
(PID.TID 0000.0001)                -3.150000000000000E+02,      /* K =  6 */
(PID.TID 0000.0001)                -3.900000000000000E+02,      /* K =  7 */
(PID.TID 0000.0001)                -4.775000000000000E+02,      /* K =  8 */
(PID.TID 0000.0001)                -5.850000000000000E+02,      /* K =  9 */
(PID.TID 0000.0001)                -7.225000000000000E+02,      /* K = 10 */
(PID.TID 0000.0001)                -9.000000000000000E+02,      /* K = 11 */
(PID.TID 0000.0001)                -1.130000000000000E+03,      /* K = 12 */
(PID.TID 0000.0001)                -1.420000000000000E+03,      /* K = 13 */
(PID.TID 0000.0001)                -1.780000000000000E+03,      /* K = 14 */
(PID.TID 0000.0001)                -2.220000000000000E+03,      /* K = 15 */
(PID.TID 0000.0001)                -2.745000000000000E+03,      /* K = 16 */
(PID.TID 0000.0001)                -3.357500000000000E+03,      /* K = 17 */
(PID.TID 0000.0001)                -4.047500000000000E+03,      /* K = 18 */
(PID.TID 0000.0001)                -4.797500000000000E+03,      /* K = 19 */
(PID.TID 0000.0001)                -5.592500000000000E+03       /* K = 20 */
(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.000000000000000E+02,      /* K =  3 */
(PID.TID 0000.0001)                -1.550000000000000E+02,      /* K =  4 */
(PID.TID 0000.0001)                -2.150000000000000E+02,      /* K =  5 */
(PID.TID 0000.0001)                -2.800000000000000E+02,      /* K =  6 */
(PID.TID 0000.0001)                -3.500000000000000E+02,      /* K =  7 */
(PID.TID 0000.0001)                -4.300000000000000E+02,      /* K =  8 */
(PID.TID 0000.0001)                -5.250000000000000E+02,      /* K =  9 */
(PID.TID 0000.0001)                -6.450000000000000E+02,      /* K = 10 */
(PID.TID 0000.0001)                -8.000000000000000E+02,      /* K = 11 */
(PID.TID 0000.0001)                -1.000000000000000E+03,      /* K = 12 */
(PID.TID 0000.0001)                -1.260000000000000E+03,      /* K = 13 */
(PID.TID 0000.0001)                -1.580000000000000E+03,      /* K = 14 */
(PID.TID 0000.0001)                -1.980000000000000E+03,      /* K = 15 */
(PID.TID 0000.0001)                -2.460000000000000E+03,      /* K = 16 */
(PID.TID 0000.0001)                -3.030000000000000E+03,      /* K = 17 */
(PID.TID 0000.0001)                -3.685000000000000E+03,      /* K = 18 */
(PID.TID 0000.0001)                -4.410000000000000E+03,      /* K = 19 */
(PID.TID 0000.0001)                -5.185000000000000E+03,      /* K = 20 */
(PID.TID 0000.0001)                -6.000000000000000E+03       /* K = 21 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deepFacC = /* deep-model grid factor @ cell-Center (-) */
(PID.TID 0000.0001)    20 @  1.000000000000000E+00              /* K =  1: 20 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deepFacF = /* deep-model grid factor @ W-Interface (-) */
(PID.TID 0000.0001)    21 @  1.000000000000000E+00              /* K =  1: 21 */
(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.335490855875281E+14
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rAc_3dMean = /* 3-D Averaged grid-cell Area (m^2) */
(PID.TID 0000.0001)                 1.557727218313311E+11
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) n2dWetPts = /* Number of wet surface points (-) */
(PID.TID 0000.0001)                 2.158000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) n3dWetPts = /* Number of wet grid points (-) */
(PID.TID 0000.0001)                 3.799400000000000E+04
(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)                   T
(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)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_isopycK =    /* Background Isopyc. Diffusivity [m^2/s] */
(PID.TID 0000.0001)                 1.100000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_advec*K =    /* Backg. GM-Advec(=Bolus) Diffusivity [m^2/s] */
(PID.TID 0000.0001)                 9.000000000000000E+02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_skewflx*K =  /* Background GM_SkewFlx Diffusivity [m^2/s] */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_isoFac_calcK = /* Fraction of dynamic K added to Redi tensor */
(PID.TID 0000.0001)                 1.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)                 9.999999999999999E-21
(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) GM_useGEOM = /* using GEOMETRIC */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_CHECK: #define ALLOW_PTRACERS
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) // PTRACERS parameters
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) PTRACERS_numInUse = /* number of tracers */
(PID.TID 0000.0001)                       1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_Iter0 = /* timestep number when tracers are initialized */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_startAllTrc =/* all tracers start @ startTime */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_doAB_onGpTr =/* apply AB on Tendencies (rather than on Tracers) */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_addSrelax2EmP =/* add Salt relaxation to EmP */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_dTLev =   /* Ptracer timestep ( s ) */
(PID.TID 0000.0001)    20 @  8.640000000000000E+04              /* K =  1: 20 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_monitorFreq = /* Frequency^-1 for monitor output (s) */
(PID.TID 0000.0001)                 4.320000000000000E+04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_dumpFreq = /* Frequency^-1 for snapshot output (s) */
(PID.TID 0000.0001)                 8.640000000000000E+04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_useRecords = /* all tracers in 1 file */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_snapshot_mnc = /* use MNC for snapshot output */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_pickup_write_mnc = /* use MNC for writing pickups */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_pickup_read_mnc = /* use MNC for reading pickups */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001)  -----------------------------------
(PID.TID 0000.0001)  tracer number :    1
(PID.TID 0000.0001) PTRACERS_names = /* Tracer short name */
(PID.TID 0000.0001)               ''
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_long_names = /* Tracer long name */
(PID.TID 0000.0001)               ''
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_ioLabel = /* tracer IO Label */
(PID.TID 0000.0001)               '01'
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_advScheme = /* Advection Scheme */
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_SOM_Advection = /* tracer uses SOM advection scheme */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_ImplVertAdv = /* implicit vert. advection flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_MultiDimAdv = /* tracer uses Multi-Dim advection */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_AdamsBashGtr = /* apply AB on tracer tendency */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_AdamsBash_Tr = /* apply AB on passive tracer */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_diffKh = /* Laplacian Diffusivity */
(PID.TID 0000.0001)                 1.000000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_diffK4 = /* Biharmonic Diffusivity */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_diffKrNr = /* Vertical Diffusivity */
(PID.TID 0000.0001)    20 @  3.000000000000000E-05              /* K =  1: 20 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_useGMRedi = /* apply GM-Redi */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_useDWNSLP = /* apply DOWN-SLOPE Flow */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_useKPP = /* apply KPP scheme */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_ref = /* Reference vertical profile */
(PID.TID 0000.0001)                 1.000000000000000E+00,      /* K =  1 */
(PID.TID 0000.0001)    19 @  0.000000000000000E+00              /* K =  2: 20 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) PTRACERS_EvPrRn =/* tracer conc. in Evap. & Rain */
(PID.TID 0000.0001)                 1.234567000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001)  -----------------------------------
(PID.TID 0000.0001) CTRL_CHECK:  --> Starts to check CTRL set-up
(PID.TID 0000.0001) CTRL_CHECK: pTracer #   1 in Gen-CTRL
(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) 
(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 :                      1
(PID.TID 0000.0001)   matching CTRL xx_file:       "xx_ptr1"
(PID.TID 0000.0001)   eps =                         1.000E-04
(PID.TID 0000.0001)   First location:                       0
(PID.TID 0000.0001)   Last location:                        4
(PID.TID 0000.0001)   Increment:                            1
(PID.TID 0000.0001)   grdchkWhichProc:                      0
(PID.TID 0000.0001)   iLocTile =      1 ,   jLocTile =      1
(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) 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) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) %MON time_tsnumber                =                     0
(PID.TID 0000.0001) %MON time_secondsf                =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_eta_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_eta_min              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_eta_mean             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_eta_sd               =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_eta_del2             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_uvel_max             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_uvel_min             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_uvel_mean            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_uvel_sd              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_uvel_del2            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_vvel_max             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_vvel_min             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_vvel_mean            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_vvel_sd              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_vvel_del2            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_wvel_max             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_wvel_min             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_wvel_mean            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_wvel_sd              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_wvel_del2            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9329580307007E+01
(PID.TID 0000.0001) %MON dynstat_theta_min            =  -1.8496928215027E+00
(PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.4434668483202E+00
(PID.TID 0000.0001) %MON dynstat_theta_sd             =   4.2513626348502E+00
(PID.TID 0000.0001) %MON dynstat_theta_del2           =   5.1755281830340E-03
(PID.TID 0000.0001) %MON dynstat_salt_max             =   3.7439460754395E+01
(PID.TID 0000.0001) %MON dynstat_salt_min             =   3.1028352737427E+01
(PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4718545511620E+01
(PID.TID 0000.0001) %MON dynstat_salt_sd              =   2.8229416934120E-01
(PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.2249504761060E-04
(PID.TID 0000.0001) %MON forcing_qnet_max             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_min             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_mean            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_sd              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_del2            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_min              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_mean             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_sd               =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_del2             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_max            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_min            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_mean           =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_sd             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_del2           =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_fu_max               =   1.6004334390163E-01
(PID.TID 0000.0001) %MON forcing_fu_min               =  -1.0329221189022E-01
(PID.TID 0000.0001) %MON forcing_fu_mean              =   4.1280488583042E-03
(PID.TID 0000.0001) %MON forcing_fu_sd                =   5.5800359441598E-02
(PID.TID 0000.0001) %MON forcing_fu_del2              =   3.2370368543277E-04
(PID.TID 0000.0001) %MON forcing_fv_max               =   8.4633164107800E-02
(PID.TID 0000.0001) %MON forcing_fv_min               =  -7.0490099489689E-02
(PID.TID 0000.0001) %MON forcing_fv_mean              =   2.8803575798275E-03
(PID.TID 0000.0001) %MON forcing_fv_sd                =   2.2193297771215E-02
(PID.TID 0000.0001) %MON forcing_fv_del2              =   2.1170835219132E-04
(PID.TID 0000.0001) %MON trAdv_CFL_u_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON trAdv_CFL_v_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON trAdv_CFL_w_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_uvel_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_vvel_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_wvel_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_W_hf_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON pe_b_mean                    =   0.0000000000000E+00
(PID.TID 0000.0001) %MON ke_max                       =   0.0000000000000E+00
(PID.TID 0000.0001) %MON ke_mean                      =   0.0000000000000E+00
(PID.TID 0000.0001) %MON ke_vol                       =   1.4474490094807E+18
(PID.TID 0000.0001) %MON vort_r_min                   =   0.0000000000000E+00
(PID.TID 0000.0001) %MON vort_r_max                   =   0.0000000000000E+00
(PID.TID 0000.0001) %MON vort_a_mean                  =  -1.8620194793381E-05
(PID.TID 0000.0001) %MON vort_a_sd                    =   7.3819348238724E-05
(PID.TID 0000.0001) %MON vort_p_mean                  =  -2.1817369821813E-05
(PID.TID 0000.0001) %MON vort_p_sd                    =   1.0010250762040E-04
(PID.TID 0000.0001) %MON surfExpan_theta_mean         =   0.0000000000000E+00
(PID.TID 0000.0001) %MON surfExpan_salt_mean          =   0.0000000000000E+00
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End MONITOR dynamic field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Begin MONITOR ptracer field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) %MON trcstat_ptracer01_max        =   1.0000000000000E+00
(PID.TID 0000.0001) %MON trcstat_ptracer01_min        =   0.0000000000000E+00
(PID.TID 0000.0001) %MON trcstat_ptracer01_mean       =   1.1521963240253E-02
(PID.TID 0000.0001) %MON trcstat_ptracer01_sd         =   1.0672023052516E-01
(PID.TID 0000.0001) %MON trcstat_ptracer01_del2       =   0.0000000000000E+00
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End MONITOR ptracer field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
(PID.TID 0000.0001)      cg2d_init_res =   3.25961757826532E+00
(PID.TID 0000.0001)      cg2d_iters(min,last) =      -1     186
(PID.TID 0000.0001)      cg2d_last_res =   9.82328669371588E-18
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) %MON time_tsnumber                =                     1
(PID.TID 0000.0001) %MON time_secondsf                =   8.6400000000000E+04
(PID.TID 0000.0001) %MON dynstat_eta_max              =   1.6743052600423E+00
(PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.3455863869851E+00
(PID.TID 0000.0001) %MON dynstat_eta_mean             =  -1.4053403839328E-18
(PID.TID 0000.0001) %MON dynstat_eta_sd               =   4.3408438286638E-01
(PID.TID 0000.0001) %MON dynstat_eta_del2             =   3.4693072840010E-03
(PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.6568006879409E-02
(PID.TID 0000.0001) %MON dynstat_uvel_min             =  -9.9209029407204E-02
(PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -1.0862912534280E-03
(PID.TID 0000.0001) %MON dynstat_uvel_sd              =   3.1095830800429E-03
(PID.TID 0000.0001) %MON dynstat_uvel_del2            =   3.2885831448895E-05
(PID.TID 0000.0001) %MON dynstat_vvel_max             =   9.0871363656057E-02
(PID.TID 0000.0001) %MON dynstat_vvel_min             =  -1.3005606017108E-01
(PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.2707164556280E-03
(PID.TID 0000.0001) %MON dynstat_vvel_sd              =   4.9820272816924E-03
(PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.0045244193533E-05
(PID.TID 0000.0001) %MON dynstat_wvel_max             =   4.4119853376755E-04
(PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.9804148516926E-04
(PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -9.5396253396687E-08
(PID.TID 0000.0001) %MON dynstat_wvel_sd              =   1.7868311507756E-05
(PID.TID 0000.0001) %MON dynstat_wvel_del2            =   1.9816405776776E-07
(PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9327980755348E+01
(PID.TID 0000.0001) %MON dynstat_theta_min            =  -1.8513879261731E+00
(PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.4434668483202E+00
(PID.TID 0000.0001) %MON dynstat_theta_sd             =   4.2510831025359E+00
(PID.TID 0000.0001) %MON dynstat_theta_del2           =   5.1683993482881E-03
(PID.TID 0000.0001) %MON dynstat_salt_max             =   3.7438512393600E+01
(PID.TID 0000.0001) %MON dynstat_salt_min             =   3.1031060186620E+01
(PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4718545511620E+01
(PID.TID 0000.0001) %MON dynstat_salt_sd              =   2.8219507328089E-01
(PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.2344897487274E-04
(PID.TID 0000.0001) %MON forcing_qnet_max             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_min             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_mean            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_sd              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_del2            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_min              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_mean             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_sd               =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_del2             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_max            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_min            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_mean           =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_sd             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_del2           =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_fu_max               =   1.6004334390163E-01
(PID.TID 0000.0001) %MON forcing_fu_min               =  -1.0329221189022E-01
(PID.TID 0000.0001) %MON forcing_fu_mean              =   4.1280488583042E-03
(PID.TID 0000.0001) %MON forcing_fu_sd                =   5.5800359441598E-02
(PID.TID 0000.0001) %MON forcing_fu_del2              =   3.2370368543277E-04
(PID.TID 0000.0001) %MON forcing_fv_max               =   8.4633164107800E-02
(PID.TID 0000.0001) %MON forcing_fv_min               =  -7.0490099489689E-02
(PID.TID 0000.0001) %MON forcing_fv_mean              =   2.8803575798275E-03
(PID.TID 0000.0001) %MON forcing_fv_sd                =   2.2193297771215E-02
(PID.TID 0000.0001) %MON forcing_fv_del2              =   2.1170835219132E-04
(PID.TID 0000.0001) %MON trAdv_CFL_u_max              =   2.6977069407749E-02
(PID.TID 0000.0001) %MON trAdv_CFL_v_max              =   1.1031754027236E-02
(PID.TID 0000.0001) %MON trAdv_CFL_w_max              =   2.8841667303896E-02
(PID.TID 0000.0001) %MON advcfl_uvel_max              =   3.4389055671184E-02
(PID.TID 0000.0001) %MON advcfl_vvel_max              =   2.5267806642809E-02
(PID.TID 0000.0001) %MON advcfl_wvel_max              =   1.5537798739418E-01
(PID.TID 0000.0001) %MON advcfl_W_hf_max              =   1.7367998966642E-01
(PID.TID 0000.0001) %MON pe_b_mean                    =   2.1298244853148E-04
(PID.TID 0000.0001) %MON ke_max                       =   4.9157938604515E-03
(PID.TID 0000.0001) %MON ke_mean                      =   1.6900689799910E-05
(PID.TID 0000.0001) %MON ke_vol                       =   1.4474490094807E+18
(PID.TID 0000.0001) %MON vort_r_min                   =  -3.0400601742071E-07
(PID.TID 0000.0001) %MON vort_r_max                   =   2.9340104619709E-07
(PID.TID 0000.0001) %MON vort_a_mean                  =  -1.8620194793382E-05
(PID.TID 0000.0001) %MON vort_a_sd                    =   7.3819142035417E-05
(PID.TID 0000.0001) %MON vort_p_mean                  =  -2.1817369821813E-05
(PID.TID 0000.0001) %MON vort_p_sd                    =   1.0010207052842E-04
(PID.TID 0000.0001) %MON surfExpan_theta_mean         =   6.0600014278862E-06
(PID.TID 0000.0001) %MON surfExpan_salt_mean          =   2.1441792214384E-07
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End MONITOR dynamic field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Begin MONITOR ptracer field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) %MON trcstat_ptracer01_max        =   1.0042254150255E+00
(PID.TID 0000.0001) %MON trcstat_ptracer01_min        =   4.2092357168592E-79
(PID.TID 0000.0001) %MON trcstat_ptracer01_mean       =   1.1521963240253E-02
(PID.TID 0000.0001) %MON trcstat_ptracer01_sd         =   1.0643562464203E-01
(PID.TID 0000.0001) %MON trcstat_ptracer01_del2       =   2.7841954026231E-05
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End MONITOR ptracer field statistics
(PID.TID 0000.0001) // =======================================================
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
(PID.TID 0000.0001)      cg2d_init_res =   2.66238707131540E-01
(PID.TID 0000.0001)      cg2d_iters(min,last) =      -1     180
(PID.TID 0000.0001)      cg2d_last_res =   8.76220416500392E-18
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) %MON time_tsnumber                =                     2
(PID.TID 0000.0001) %MON time_secondsf                =   1.7280000000000E+05
(PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1722938479294E+00
(PID.TID 0000.0001) %MON dynstat_eta_min              =  -2.1180052058260E+00
(PID.TID 0000.0001) %MON dynstat_eta_mean             =   5.9961189714465E-17
(PID.TID 0000.0001) %MON dynstat_eta_sd               =   6.7566276976025E-01
(PID.TID 0000.0001) %MON dynstat_eta_del2             =   3.5963390736788E-03
(PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.5349909649340E-01
(PID.TID 0000.0001) %MON dynstat_uvel_min             =  -1.9061986513844E-01
(PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -1.6400805887536E-03
(PID.TID 0000.0001) %MON dynstat_uvel_sd              =   5.7854632395256E-03
(PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.1858910099046E-05
(PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.7704435114990E-01
(PID.TID 0000.0001) %MON dynstat_vvel_min             =  -2.5435925132025E-01
(PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.4264608498679E-03
(PID.TID 0000.0001) %MON dynstat_vvel_sd              =   8.8549455675062E-03
(PID.TID 0000.0001) %MON dynstat_vvel_del2            =   7.5495195703149E-05
(PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.8364156103437E-04
(PID.TID 0000.0001) %MON dynstat_wvel_min             =  -5.9524948848794E-04
(PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -6.5561960226668E-08
(PID.TID 0000.0001) %MON dynstat_wvel_sd              =   3.4064111608528E-05
(PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.7485366218700E-07
(PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9326233278365E+01
(PID.TID 0000.0001) %MON dynstat_theta_min            =  -1.8491828881337E+00
(PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.4433949924863E+00
(PID.TID 0000.0001) %MON dynstat_theta_sd             =   4.2508207215487E+00
(PID.TID 0000.0001) %MON dynstat_theta_del2           =   5.1215119746026E-03
(PID.TID 0000.0001) %MON dynstat_salt_max             =   3.7437632997837E+01
(PID.TID 0000.0001) %MON dynstat_salt_min             =   3.1033203489564E+01
(PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4718542823795E+01
(PID.TID 0000.0001) %MON dynstat_salt_sd              =   2.8211408444131E-01
(PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.2272810197565E-04
(PID.TID 0000.0001) %MON forcing_qnet_max             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_min             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_mean            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_sd              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_del2            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_min              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_mean             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_sd               =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_del2             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_max            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_min            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_mean           =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_sd             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_del2           =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_fu_max               =   1.6004334390163E-01
(PID.TID 0000.0001) %MON forcing_fu_min               =  -1.0329221189022E-01
(PID.TID 0000.0001) %MON forcing_fu_mean              =   4.1280488583042E-03
(PID.TID 0000.0001) %MON forcing_fu_sd                =   5.5800359441598E-02
(PID.TID 0000.0001) %MON forcing_fu_del2              =   3.2370368543277E-04
(PID.TID 0000.0001) %MON forcing_fv_max               =   8.4633164107800E-02
(PID.TID 0000.0001) %MON forcing_fv_min               =  -7.0490099489689E-02
(PID.TID 0000.0001) %MON forcing_fv_mean              =   2.8803575798275E-03
(PID.TID 0000.0001) %MON forcing_fv_sd                =   2.2193297771215E-02
(PID.TID 0000.0001) %MON forcing_fv_del2              =   2.1170835219132E-04
(PID.TID 0000.0001) %MON trAdv_CFL_u_max              =   3.4396042957797E-02
(PID.TID 0000.0001) %MON trAdv_CFL_v_max              =   2.5477453623859E-02
(PID.TID 0000.0001) %MON trAdv_CFL_w_max              =   1.4418143360672E-01
(PID.TID 0000.0001) %MON advcfl_uvel_max              =   5.9828840721428E-02
(PID.TID 0000.0001) %MON advcfl_vvel_max              =   4.9417923099586E-02
(PID.TID 0000.0001) %MON advcfl_wvel_max              =   3.1131438583132E-01
(PID.TID 0000.0001) %MON advcfl_W_hf_max              =   3.4786496603512E-01
(PID.TID 0000.0001) %MON pe_b_mean                    =   5.1600685488463E-04
(PID.TID 0000.0001) %MON ke_max                       =   1.8344723279920E-02
(PID.TID 0000.0001) %MON ke_mean                      =   5.2871985504159E-05
(PID.TID 0000.0001) %MON ke_vol                       =   1.4474490094807E+18
(PID.TID 0000.0001) %MON vort_r_min                   =  -5.9016503193913E-07
(PID.TID 0000.0001) %MON vort_r_max                   =   5.6830858846864E-07
(PID.TID 0000.0001) %MON vort_a_mean                  =  -1.8620194793382E-05
(PID.TID 0000.0001) %MON vort_a_sd                    =   7.3819041903640E-05
(PID.TID 0000.0001) %MON vort_p_mean                  =  -2.1817354111469E-05
(PID.TID 0000.0001) %MON vort_p_sd                    =   1.0010698930140E-04
(PID.TID 0000.0001) %MON surfExpan_theta_mean         =   4.1467193481611E-06
(PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.8391512493633E-07
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End MONITOR dynamic field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Begin MONITOR ptracer field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) %MON trcstat_ptracer01_max        =   1.0135028593422E+00
(PID.TID 0000.0001) %MON trcstat_ptracer01_min        =  -2.3477632371764E-02
(PID.TID 0000.0001) %MON trcstat_ptracer01_mean       =   1.1521869985534E-02
(PID.TID 0000.0001) %MON trcstat_ptracer01_sd         =   1.0617001604937E-01
(PID.TID 0000.0001) %MON trcstat_ptracer01_del2       =   4.4700606541972E-05
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End MONITOR ptracer field statistics
(PID.TID 0000.0001) // =======================================================
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
(PID.TID 0000.0001)      cg2d_init_res =   3.92551697890786E-01
(PID.TID 0000.0001)      cg2d_iters(min,last) =      -1     182
(PID.TID 0000.0001)      cg2d_last_res =   7.51755496846236E-18
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) %MON time_tsnumber                =                     3
(PID.TID 0000.0001) %MON time_secondsf                =   2.5920000000000E+05
(PID.TID 0000.0001) %MON dynstat_eta_max              =   2.4140754918099E+00
(PID.TID 0000.0001) %MON dynstat_eta_min              =  -2.3612137980946E+00
(PID.TID 0000.0001) %MON dynstat_eta_mean             =  -1.7238842042909E-16
(PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6426133342661E-01
(PID.TID 0000.0001) %MON dynstat_eta_del2             =   3.3632240946691E-03
(PID.TID 0000.0001) %MON dynstat_uvel_max             =   2.2141094835870E-01
(PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.7588418338258E-01
(PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -1.6665508561476E-03
(PID.TID 0000.0001) %MON dynstat_uvel_sd              =   8.3985173092517E-03
(PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.6142553628778E-05
(PID.TID 0000.0001) %MON dynstat_vvel_max             =   2.6054430028602E-01
(PID.TID 0000.0001) %MON dynstat_vvel_min             =  -3.7485964306326E-01
(PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.1207885178596E-03
(PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.2354962037809E-02
(PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.0292211707282E-04
(PID.TID 0000.0001) %MON dynstat_wvel_max             =   1.3121798223018E-03
(PID.TID 0000.0001) %MON dynstat_wvel_min             =  -8.7994735809658E-04
(PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -3.1338158955944E-08
(PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.8185884395796E-05
(PID.TID 0000.0001) %MON dynstat_wvel_del2            =   5.1970403022149E-07
(PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9324393327545E+01
(PID.TID 0000.0001) %MON dynstat_theta_min            =  -1.8433200099543E+00
(PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.4434166596962E+00
(PID.TID 0000.0001) %MON dynstat_theta_sd             =   4.2506397887185E+00
(PID.TID 0000.0001) %MON dynstat_theta_del2           =   5.0622124453523E-03
(PID.TID 0000.0001) %MON dynstat_salt_max             =   3.7436798088433E+01
(PID.TID 0000.0001) %MON dynstat_salt_min             =   3.1035935877654E+01
(PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4718542808429E+01
(PID.TID 0000.0001) %MON dynstat_salt_sd              =   2.8204831087947E-01
(PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.2068695865780E-04
(PID.TID 0000.0001) %MON forcing_qnet_max             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_min             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_mean            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_sd              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_del2            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_min              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_mean             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_sd               =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_del2             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_max            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_min            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_mean           =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_sd             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_del2           =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_fu_max               =   1.6004334390163E-01
(PID.TID 0000.0001) %MON forcing_fu_min               =  -1.0329221189022E-01
(PID.TID 0000.0001) %MON forcing_fu_mean              =   4.1280488583042E-03
(PID.TID 0000.0001) %MON forcing_fu_sd                =   5.5800359441598E-02
(PID.TID 0000.0001) %MON forcing_fu_del2              =   3.2370368543277E-04
(PID.TID 0000.0001) %MON forcing_fv_max               =   8.4633164107800E-02
(PID.TID 0000.0001) %MON forcing_fv_min               =  -7.0490099489689E-02
(PID.TID 0000.0001) %MON forcing_fv_mean              =   2.8803575798275E-03
(PID.TID 0000.0001) %MON forcing_fv_sd                =   2.2193297771215E-02
(PID.TID 0000.0001) %MON forcing_fv_del2              =   2.1170835219132E-04
(PID.TID 0000.0001) %MON trAdv_CFL_u_max              =   5.9849338443010E-02
(PID.TID 0000.0001) %MON trAdv_CFL_v_max              =   4.9923946694845E-02
(PID.TID 0000.0001) %MON trAdv_CFL_w_max              =   2.8742625017039E-01
(PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.7134445371947E-02
(PID.TID 0000.0001) %MON advcfl_vvel_max              =   7.2829216621317E-02
(PID.TID 0000.0001) %MON advcfl_wvel_max              =   4.5837464677179E-01
(PID.TID 0000.0001) %MON advcfl_W_hf_max              =   5.1296969109274E-01
(PID.TID 0000.0001) %MON pe_b_mean                    =   6.6020569779455E-04
(PID.TID 0000.0001) %MON ke_max                       =   3.9075413857355E-02
(PID.TID 0000.0001) %MON ke_mean                      =   1.0304434708380E-04
(PID.TID 0000.0001) %MON ke_vol                       =   1.4474490094807E+18
(PID.TID 0000.0001) %MON vort_r_min                   =  -8.6228013318125E-07
(PID.TID 0000.0001) %MON vort_r_max                   =   8.3600040616677E-07
(PID.TID 0000.0001) %MON vort_a_mean                  =  -1.8620194793382E-05
(PID.TID 0000.0001) %MON vort_a_sd                    =   7.3819044706365E-05
(PID.TID 0000.0001) %MON vort_p_mean                  =  -2.1817357372569E-05
(PID.TID 0000.0001) %MON vort_p_sd                    =   1.0011062181928E-04
(PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.7664004093567E-06
(PID.TID 0000.0001) %MON surfExpan_salt_mean          =   7.4273904249391E-08
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End MONITOR dynamic field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Begin MONITOR ptracer field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) %MON trcstat_ptracer01_max        =   1.0365841835399E+00
(PID.TID 0000.0001) %MON trcstat_ptracer01_min        =  -6.3110104723254E-02
(PID.TID 0000.0001) %MON trcstat_ptracer01_mean       =   1.1521451706624E-02
(PID.TID 0000.0001) %MON trcstat_ptracer01_sd         =   1.0591631665363E-01
(PID.TID 0000.0001) %MON trcstat_ptracer01_del2       =   6.1213823969941E-05
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End MONITOR ptracer field statistics
(PID.TID 0000.0001) // =======================================================
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001)      cg2d_init_res =   4.81371497819798E-01
(PID.TID 0000.0001)      cg2d_iters(min,last) =      -1     182
(PID.TID 0000.0001)      cg2d_last_res =   8.51526785967628E-18
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) %MON time_tsnumber                =                     4
(PID.TID 0000.0001) %MON time_secondsf                =   3.4560000000000E+05
(PID.TID 0000.0001) %MON dynstat_eta_max              =   2.4965434427473E+00
(PID.TID 0000.0001) %MON dynstat_eta_min              =  -2.3353439381808E+00
(PID.TID 0000.0001) %MON dynstat_eta_mean             =   6.7456338428773E-17
(PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.7757200900409E-01
(PID.TID 0000.0001) %MON dynstat_eta_del2             =   3.1485297038805E-03
(PID.TID 0000.0001) %MON dynstat_uvel_max             =   2.7975586143468E-01
(PID.TID 0000.0001) %MON dynstat_uvel_min             =  -3.5403346233508E-01
(PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -1.4166161294004E-03
(PID.TID 0000.0001) %MON dynstat_uvel_sd              =   1.1136898417508E-02
(PID.TID 0000.0001) %MON dynstat_uvel_del2            =   1.0611319362049E-04
(PID.TID 0000.0001) %MON dynstat_vvel_max             =   3.4070558210048E-01
(PID.TID 0000.0001) %MON dynstat_vvel_min             =  -4.8895811484273E-01
(PID.TID 0000.0001) %MON dynstat_vvel_mean            =   6.4330455157150E-04
(PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.5264434258407E-02
(PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.2184250405457E-04
(PID.TID 0000.0001) %MON dynstat_wvel_max             =   1.7199202381894E-03
(PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.1447923391891E-03
(PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -9.0867344519041E-09
(PID.TID 0000.0001) %MON dynstat_wvel_sd              =   6.0062696933672E-05
(PID.TID 0000.0001) %MON dynstat_wvel_del2            =   6.3099042930147E-07
(PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9322579385039E+01
(PID.TID 0000.0001) %MON dynstat_theta_min            =  -1.8143597695966E+00
(PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.4434413012240E+00
(PID.TID 0000.0001) %MON dynstat_theta_sd             =   4.2504854943091E+00
(PID.TID 0000.0001) %MON dynstat_theta_del2           =   5.0063053064132E-03
(PID.TID 0000.0001) %MON dynstat_salt_max             =   3.7436005928498E+01
(PID.TID 0000.0001) %MON dynstat_salt_min             =   3.1038905895788E+01
(PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4718543405790E+01
(PID.TID 0000.0001) %MON dynstat_salt_sd              =   2.8199011833275E-01
(PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.1826785630806E-04
(PID.TID 0000.0001) %MON forcing_qnet_max             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_min             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_mean            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_sd              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qnet_del2            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_max              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_min              =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_mean             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_sd               =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_qsw_del2             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_max            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_min            =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_mean           =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_sd             =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_empmr_del2           =   0.0000000000000E+00
(PID.TID 0000.0001) %MON forcing_fu_max               =   1.6004334390163E-01
(PID.TID 0000.0001) %MON forcing_fu_min               =  -1.0329221189022E-01
(PID.TID 0000.0001) %MON forcing_fu_mean              =   4.1280488583042E-03
(PID.TID 0000.0001) %MON forcing_fu_sd                =   5.5800359441598E-02
(PID.TID 0000.0001) %MON forcing_fu_del2              =   3.2370368543277E-04
(PID.TID 0000.0001) %MON forcing_fv_max               =   8.4633164107800E-02
(PID.TID 0000.0001) %MON forcing_fv_min               =  -7.0490099489689E-02
(PID.TID 0000.0001) %MON forcing_fv_mean              =   2.8803575798275E-03
(PID.TID 0000.0001) %MON forcing_fv_sd                =   2.2193297771215E-02
(PID.TID 0000.0001) %MON forcing_fv_del2              =   2.1170835219132E-04
(PID.TID 0000.0001) %MON trAdv_CFL_u_max              =   7.7161954703140E-02
(PID.TID 0000.0001) %MON trAdv_CFL_v_max              =   7.3449834843828E-02
(PID.TID 0000.0001) %MON trAdv_CFL_w_max              =   4.2232817386396E-01
(PID.TID 0000.0001) %MON advcfl_uvel_max              =   8.5910483431359E-02
(PID.TID 0000.0001) %MON advcfl_vvel_max              =   9.4996719768585E-02
(PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.9435516508114E-01
(PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.6665346136485E-01
(PID.TID 0000.0001) %MON pe_b_mean                    =   6.8340276198649E-04
(PID.TID 0000.0001) %MON ke_max                       =   6.5658211688714E-02
(PID.TID 0000.0001) %MON ke_mean                      =   1.6308888218461E-04
(PID.TID 0000.0001) %MON ke_vol                       =   1.4474490094807E+18
(PID.TID 0000.0001) %MON vort_r_min                   =  -1.1142780134389E-06
(PID.TID 0000.0001) %MON vort_r_max                   =   1.0945056568835E-06
(PID.TID 0000.0001) %MON vort_a_mean                  =  -1.8620194793382E-05
(PID.TID 0000.0001) %MON vort_a_sd                    =   7.3819102578469E-05
(PID.TID 0000.0001) %MON vort_p_mean                  =  -2.1817369223887E-05
(PID.TID 0000.0001) %MON vort_p_sd                    =   1.0011199422819E-04
(PID.TID 0000.0001) %MON surfExpan_theta_mean         =   3.6636405801389E-07
(PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -1.3479900611624E-08
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End MONITOR dynamic field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Begin MONITOR ptracer field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) %MON trcstat_ptracer01_max        =   1.0640605525151E+00
(PID.TID 0000.0001) %MON trcstat_ptracer01_min        =  -1.1381665492144E-01
(PID.TID 0000.0001) %MON trcstat_ptracer01_mean       =   1.1520685994679E-02
(PID.TID 0000.0001) %MON trcstat_ptracer01_sd         =   1.0567383186993E-01
(PID.TID 0000.0001) %MON trcstat_ptracer01_del2       =   7.7833603450913E-05
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End MONITOR ptracer field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) %CHECKPOINT         4 ckptA
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471774570D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723380947D+15
(PID.TID 0000.0001)  global fc =  0.184102723380947D+15
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
 cg2d: Sum(rhs),rhsMax =   0.00000000000000E+00  0.00000000000000E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =   3.44169137633799E-15  1.47881353008629E+03
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =   2.55351295663786E-15  7.98899709494860E+03
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.02695629777827E-15  2.20410845001160E+04
(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) 
(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) 
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471774570D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723380947D+15
(PID.TID 0000.0001)  global fc =  0.184102723380947D+15
(PID.TID 0000.0001) grdchk reference fc: fcref       =  1.84102723380947E+14
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 1881   15    5    4    1    1
(PID.TID 0000.0001) ====== Starts gradient-check number   1 (=ichknum) =======
 ph-test icomp, ncvarcomp, ichknum         1881       37994           1
 ph-grd _loc: bi, bj, icomptest, ichknum            1           1           0           1
 ph-grd _loc: bi, bj, icomptest, ichknum            1           1         616           1
 ph-grd _loc: bi, bj, icomptest, ichknum            1           1        1232           1
 ph-grd _loc: bi, bj, icomptest, ichknum            1           1        1848           1
 ph-grd -->hit<--           15           5           4           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   15    5    4 ; bi,bj=   1   1 ; 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) 
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471742246D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723377715D+15
(PID.TID 0000.0001)  global fc =  0.184102723377715D+15
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  1.84102723377715E+14
(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) 
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471806894D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723384179D+15
(PID.TID 0000.0001)  global fc =  0.184102723384179D+15
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  1.84102723384179E+14
grad-res -------------------------------
 grad-res     0    1   15    5    4    1    1    1   1.84102723381E+14  1.84102723378E+14  1.84102723384E+14
 grad-res     0    1    1 1881    0    1    1    1  -3.23242721513E+07 -3.23242187500E+07  1.65204948521E-06
(PID.TID 0000.0001)  ADM  ref_cost_function      =  1.84102723380947E+14
(PID.TID 0000.0001)  ADM  adjoint_gradient       = -3.23242721512972E+07
(PID.TID 0000.0001)  ADM  finite-diff_grad       = -3.23242187500000E+07
(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         1882       37994           2
 ph-grd _loc: bi, bj, icomptest, ichknum            1           1        1881           2
 ph-grd -->hit<--           16           5           4           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   16    5    4 ; bi,bj=   1   1 ; 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) 
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471753096D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723378800D+15
(PID.TID 0000.0001)  global fc =  0.184102723378800D+15
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  1.84102723378800E+14
(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) 
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471796044D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723383094D+15
(PID.TID 0000.0001)  global fc =  0.184102723383094D+15
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  1.84102723383094E+14
grad-res -------------------------------
 grad-res     0    2   16    5    4    1    1    1   1.84102723381E+14  1.84102723379E+14  1.84102723383E+14
 grad-res     0    2    2 1882    0    1    1    1  -2.14737521330E+07 -2.14737500000E+07  9.93301508823E-08
(PID.TID 0000.0001)  ADM  ref_cost_function      =  1.84102723380947E+14
(PID.TID 0000.0001)  ADM  adjoint_gradient       = -2.14737521329910E+07
(PID.TID 0000.0001)  ADM  finite-diff_grad       = -2.14737500000000E+07
(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         1883       37994           3
 ph-grd _loc: bi, bj, icomptest, ichknum            1           1        1882           3
 ph-grd -->hit<--           17           5           4           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   17    5    4 ; bi,bj=   1   1 ; 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) 
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471730206D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723376511D+15
(PID.TID 0000.0001)  global fc =  0.184102723376511D+15
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  1.84102723376511E+14
(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) 
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471818934D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723385383D+15
(PID.TID 0000.0001)  global fc =  0.184102723385383D+15
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  1.84102723385383E+14
grad-res -------------------------------
 grad-res     0    3   17    5    4    1    1    1   1.84102723381E+14  1.84102723377E+14  1.84102723385E+14
 grad-res     0    3    3 1883    0    1    1    1  -4.43643318051E+07 -4.43642187500E+07  2.54833250102E-06
(PID.TID 0000.0001)  ADM  ref_cost_function      =  1.84102723380947E+14
(PID.TID 0000.0001)  ADM  adjoint_gradient       = -4.43643318050686E+07
(PID.TID 0000.0001)  ADM  finite-diff_grad       = -4.43642187500000E+07
(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         1884       37994           4
 ph-grd _loc: bi, bj, icomptest, ichknum            1           1        1883           4
 ph-grd -->hit<--           18           5           4           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   18    5    4 ; bi,bj=   1   1 ; 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) 
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471743619D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723377852D+15
(PID.TID 0000.0001)  global fc =  0.184102723377852D+15
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  1.84102723377852E+14
(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) 
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471805521D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723384042D+15
(PID.TID 0000.0001)  global fc =  0.184102723384042D+15
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  1.84102723384042E+14
grad-res -------------------------------
 grad-res     0    4   18    5    4    1    1    1   1.84102723381E+14  1.84102723378E+14  1.84102723384E+14
 grad-res     0    4    4 1884    0    1    1    1  -3.09510269421E+07 -3.09510937500E+07 -2.15850315688E-06
(PID.TID 0000.0001)  ADM  ref_cost_function      =  1.84102723380947E+14
(PID.TID 0000.0001)  ADM  adjoint_gradient       = -3.09510269421106E+07
(PID.TID 0000.0001)  ADM  finite-diff_grad       = -3.09510937500000E+07
(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         1885       37994           5
 ph-grd _loc: bi, bj, icomptest, ichknum            1           1        1884           5
 ph-grd -->hit<--           19           5           4           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   19    5    4 ; bi,bj=   1   1 ; 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) 
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471741881D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723377678D+15
(PID.TID 0000.0001)  global fc =  0.184102723377678D+15
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  1.84102723377678E+14
(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) 
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -3.49720252756924E-15  8.55254939901145E+00
 cg2d: Sum(rhs),rhsMax =   1.12410081243297E-14  9.18227067226688E+00
 cg2d: Sum(rhs),rhsMax =  -5.44009282066327E-15  8.71966856788159E+00
 cg2d: Sum(rhs),rhsMax =  -2.55351295663786E-15  8.12211131921598E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
(PID.TID 0000.0001)   early fc =  0.000000000000000D+00
 --> objf_tracer(bi,bj)    =  0.513236471807259D+14
 --> objf_tracer(bi,bj)    =  0.579123484504980D+14
 --> objf_tracer(bi,bj)    =  0.267040314490043D+14
 --> objf_tracer(bi,bj)    =  0.481626963039877D+14
(PID.TID 0000.0001)   local fc =  0.184102723384216D+15
(PID.TID 0000.0001)  global fc =  0.184102723384216D+15
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  1.84102723384216E+14
grad-res -------------------------------
 grad-res     0    5   19    5    4    1    1    1   1.84102723381E+14  1.84102723378E+14  1.84102723384E+14
 grad-res     0    5    5 1885    0    1    1    1  -3.26888456348E+07 -3.26889062500E+07 -1.85430783506E-06
(PID.TID 0000.0001)  ADM  ref_cost_function      =  1.84102723380947E+14
(PID.TID 0000.0001)  ADM  adjoint_gradient       = -3.26888456348174E+07
(PID.TID 0000.0001)  ADM  finite-diff_grad       = -3.26889062500000E+07
(PID.TID 0000.0001) ====== End of gradient-check number   5 (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-04 ; grdchk CTRL var/file name: "xx_ptr1"
(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    15     5     4    1    1   0.000000000E+00 -1.000000000E-04
(PID.TID 0000.0001) grdchk output (c):   1  1.8410272338095E+14  1.8410272337771E+14  1.8410272338418E+14
(PID.TID 0000.0001) grdchk output (g):   1    -3.2324218750000E+07 -3.2324272151297E+07  1.6520494852124E-06
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   2    16     5     4    1    1   0.000000000E+00 -1.000000000E-04
(PID.TID 0000.0001) grdchk output (c):   2  1.8410272338095E+14  1.8410272337880E+14  1.8410272338309E+14
(PID.TID 0000.0001) grdchk output (g):   2    -2.1473750000000E+07 -2.1473752132991E+07  9.9330150882260E-08
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   3    17     5     4    1    1   0.000000000E+00 -1.000000000E-04
(PID.TID 0000.0001) grdchk output (c):   3  1.8410272338095E+14  1.8410272337651E+14  1.8410272338538E+14
(PID.TID 0000.0001) grdchk output (g):   3    -4.4364218750000E+07 -4.4364331805069E+07  2.5483325010178E-06
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   4    18     5     4    1    1   0.000000000E+00 -1.000000000E-04
(PID.TID 0000.0001) grdchk output (c):   4  1.8410272338095E+14  1.8410272337785E+14  1.8410272338404E+14
(PID.TID 0000.0001) grdchk output (g):   4    -3.0951093750000E+07 -3.0951026942111E+07 -2.1585031568794E-06
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   5    19     5     4    1    1   0.000000000E+00 -1.000000000E-04
(PID.TID 0000.0001) grdchk output (c):   5  1.8410272338095E+14  1.8410272337768E+14  1.8410272338422E+14
(PID.TID 0000.0001) grdchk output (g):   5    -3.2688906250000E+07 -3.2688845634817E+07 -1.8543078350586E-06
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk  summary  :  RMS of    5 ratios =  1.8617586278369E-06
(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:   22.073165626032278
(PID.TID 0000.0001)         System time:  0.43623100058175623
(PID.TID 0000.0001)     Wall clock time:   22.658195972442627
(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:   6.3411998329684138E-002
(PID.TID 0000.0001)         System time:   7.9669998958706856E-003
(PID.TID 0000.0001)     Wall clock time:   7.2613954544067383E-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:   13.582355692982674
(PID.TID 0000.0001)         System time:   1.2097969651222229E-002
(PID.TID 0000.0001)     Wall clock time:   13.667284488677979
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "UPDATE_R_STAR       [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.29996734857559204
(PID.TID 0000.0001)         System time:   2.2301077842712402E-004
(PID.TID 0000.0001)     Wall clock time:  0.30214715003967285
(PID.TID 0000.0001)          No. starts:          96
(PID.TID 0000.0001)           No. stops:          96
(PID.TID 0000.0001)   Seconds in section "LOAD_FIELDS_DRIVER  [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   2.0592808723449707E-003
(PID.TID 0000.0001)         System time:   5.6974589824676514E-005
(PID.TID 0000.0001)     Wall clock time:   2.1297931671142578E-003
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "EXTERNAL_FLDS_LOAD [LOAD_FLDS_DRIVER]":
(PID.TID 0000.0001)           User time:   4.2572617530822754E-004
(PID.TID 0000.0001)         System time:   1.4029443264007568E-005
(PID.TID 0000.0001)     Wall clock time:   4.4250488281250000E-004
(PID.TID 0000.0001)          No. starts:          52
(PID.TID 0000.0001)           No. stops:          52
(PID.TID 0000.0001)   Seconds in section "CTRL_MAP_FORCING  [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   3.5311281681060791E-004
(PID.TID 0000.0001)         System time:   1.5005469322204590E-005
(PID.TID 0000.0001)     Wall clock time:   3.7217140197753906E-004
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "DO_ATMOSPHERIC_PHYS [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   1.3741359114646912E-002
(PID.TID 0000.0001)         System time:   4.8198550939559937E-004
(PID.TID 0000.0001)     Wall clock time:   1.4318466186523438E-002
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "DO_OCEANIC_PHYS     [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   2.0580911338329315
(PID.TID 0000.0001)         System time:   2.4110227823257446E-003
(PID.TID 0000.0001)     Wall clock time:   2.0733504295349121
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "THERMODYNAMICS      [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   5.0421176403760910
(PID.TID 0000.0001)         System time:   1.1450201272964478E-003
(PID.TID 0000.0001)     Wall clock time:   5.0700449943542480
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "DYNAMICS            [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   3.1989133358001709
(PID.TID 0000.0001)         System time:   9.1300904750823975E-004
(PID.TID 0000.0001)     Wall clock time:   3.2167551517486572
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "UPDATE_CG2D         [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   7.0463120937347412E-002
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:   7.0925474166870117E-002
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "SOLVE_FOR_PRESSURE  [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   1.9801445007324219
(PID.TID 0000.0001)         System time:   9.9837779998779297E-006
(PID.TID 0000.0001)     Wall clock time:   1.9884123802185059
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "MOM_CORRECTION_STEP [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.10141700506210327
(PID.TID 0000.0001)         System time:   1.6987323760986328E-005
(PID.TID 0000.0001)     Wall clock time:  0.10196590423583984
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "INTEGR_CONTINUITY   [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.19059416651725769
(PID.TID 0000.0001)         System time:   2.9802322387695312E-006
(PID.TID 0000.0001)     Wall clock time:  0.19193315505981445
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "CALC_R_STAR         [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   2.1690011024475098E-002
(PID.TID 0000.0001)         System time:   6.0200691223144531E-006
(PID.TID 0000.0001)     Wall clock time:   2.1806478500366211E-002
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "TRC_CORRECTION_STEP [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.30075025558471680
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  0.30224227905273438
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "BLOCKING_EXCHANGES  [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.12095633149147034
(PID.TID 0000.0001)         System time:   7.0035457611083984E-006
(PID.TID 0000.0001)     Wall clock time:  0.12170648574829102
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "MONITOR             [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   7.8891873359680176E-002
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:   7.9274177551269531E-002
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "COST_TILE           [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   1.5914142131805420E-003
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:   1.6038417816162109E-003
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "DO_THE_MODEL_IO     [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   6.3320755958557129E-002
(PID.TID 0000.0001)         System time:   5.7230070233345032E-003
(PID.TID 0000.0001)     Wall clock time:   6.9384813308715820E-002
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "PTRACERS_RESET      [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   3.7530064582824707E-004
(PID.TID 0000.0001)         System time:   3.0025839805603027E-006
(PID.TID 0000.0001)     Wall clock time:   3.7240982055664062E-004
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "DO_WRITE_PICKUP     [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   2.8336584568023682E-002
(PID.TID 0000.0001)         System time:   9.7899138927459717E-004
(PID.TID 0000.0001)     Wall clock time:   3.0234813690185547E-002
(PID.TID 0000.0001)          No. starts:          48
(PID.TID 0000.0001)           No. stops:          48
(PID.TID 0000.0001)   Seconds in section "GRDCHK_MAIN         [THE_MODEL_MAIN]":
(PID.TID 0000.0001)           User time:   12.900241851806641
(PID.TID 0000.0001)         System time:   1.4885991811752319E-002
(PID.TID 0000.0001)     Wall clock time:   12.984414815902710
(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.53582382202148438
(PID.TID 0000.0001)         System time:   7.9140067100524902E-003
(PID.TID 0000.0001)     Wall clock time:  0.54688954353332520
(PID.TID 0000.0001)          No. starts:          11
(PID.TID 0000.0001)           No. stops:          11
(PID.TID 0000.0001)   Seconds in section "MAIN LOOP           [THE_MAIN_LOOP]":
(PID.TID 0000.0001)           User time:   12.284673690795898
(PID.TID 0000.0001)         System time:   1.2099742889404297E-003
(PID.TID 0000.0001)     Wall clock time:   12.351166725158691
(PID.TID 0000.0001)          No. starts:          11
(PID.TID 0000.0001)           No. stops:          11
(PID.TID 0000.0001)   Seconds in section "MAIN_DO_LOOP        [THE_MAIN_LOOP]":
(PID.TID 0000.0001)           User time:   12.283740997314453
(PID.TID 0000.0001)         System time:   1.1189877986907959E-003
(PID.TID 0000.0001)     Wall clock time:   12.350131034851074
(PID.TID 0000.0001)          No. starts:          44
(PID.TID 0000.0001)           No. stops:          44
(PID.TID 0000.0001)   Seconds in section "COST_FINAL         [ADJOINT SPIN-DOWN]":
(PID.TID 0000.0001)           User time:   4.4441223144531250E-004
(PID.TID 0000.0001)         System time:   7.9989433288574219E-005
(PID.TID 0000.0001)     Wall clock time:   5.2833557128906250E-004
(PID.TID 0000.0001)          No. starts:          11
(PID.TID 0000.0001)           No. stops:          11
(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 =          66982
(PID.TID 0000.0001) //      Max. barrier spins =              1
(PID.TID 0000.0001) //      Min. barrier spins =              1
(PID.TID 0000.0001) //     Total barrier spins =          66982
(PID.TID 0000.0001) //      Avg. barrier spins =       1.00E+00
PROGRAM MAIN: Execution ended Normally
