(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:  checkpoint64u
(PID.TID 0000.0001) // Build user:        jmc
(PID.TID 0000.0001) // Build host:        baudelaire
(PID.TID 0000.0001) // Build date:        Wed Mar 19 11:30:58 EDT 2014
(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 =    2 ; /* Tile overlap distance in X */
(PID.TID 0000.0001)      OLy =    2 ; /* Tile overlap distance in Y */
(PID.TID 0000.0001)      nTx =    1 ; /* No. threads in X per process */
(PID.TID 0000.0001)      nTy =    1 ; /* No. threads in Y per process */
(PID.TID 0000.0001)       Nr =   15 ; /* No. levels in the vertical   */
(PID.TID 0000.0001)       Nx =   90 ; /* Total domain size in X ( = nPx*nSx*sNx ) */
(PID.TID 0000.0001)       Ny =   40 ; /* Total domain size in Y ( = nPy*nSy*sNy ) */
(PID.TID 0000.0001)   nTiles =    4 ; /* Total no. tiles per process ( = nSx*nSy ) */
(PID.TID 0000.0001)   nProcs =    1 ; /* Total no. processes ( = nPx*nPy ) */
(PID.TID 0000.0001) nThreads =    1 ; /* Total no. threads per process ( = nTx*nTy ) */
(PID.TID 0000.0001) usingMPI =    F ; /* Flag used to control whether MPI is in use */
(PID.TID 0000.0001)                   /*  note: To execute a program with MPI calls */
(PID.TID 0000.0001)                   /*  it must be launched appropriately e.g     */
(PID.TID 0000.0001)                   /*  "mpirun -np 64 ......"                    */
(PID.TID 0000.0001) useCoupler=    F ;/* Flag used to control communications with   */
(PID.TID 0000.0001)                   /*  other model components, through a coupler */
(PID.TID 0000.0001) debugMode =    F ; /* print debug msg. (sequence of S/R calls)  */
(PID.TID 0000.0001) printMapIncludesZeros=    F ; /* print zeros in Std.Output maps */
(PID.TID 0000.0001) maxLengthPrt1D=   65 /* maxLength of 1D array printed to StdOut */
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // Mapping of tiles to threads
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // -o- Thread   1, tiles (   1:   2,   1:   2)
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // Tile <-> Tile connectvity table
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // Tile number: 000001 (process no. = 000000)
(PID.TID 0000.0001) //        WEST: Tile = 000002, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000001
(PID.TID 0000.0001) //        EAST: Tile = 000002, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000001
(PID.TID 0000.0001) //       SOUTH: Tile = 000003, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000002
(PID.TID 0000.0001) //       NORTH: Tile = 000003, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000002
(PID.TID 0000.0001) // Tile number: 000002 (process no. = 000000)
(PID.TID 0000.0001) //        WEST: Tile = 000001, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001) //        EAST: Tile = 000001, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001) //       SOUTH: Tile = 000004, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000002
(PID.TID 0000.0001) //       NORTH: Tile = 000004, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000002
(PID.TID 0000.0001) // Tile number: 000003 (process no. = 000000)
(PID.TID 0000.0001) //        WEST: Tile = 000004, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000002
(PID.TID 0000.0001) //        EAST: Tile = 000004, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000002
(PID.TID 0000.0001) //       SOUTH: Tile = 000001, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001) //       NORTH: Tile = 000001, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001) // Tile number: 000004 (process no. = 000000)
(PID.TID 0000.0001) //        WEST: Tile = 000003, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000002
(PID.TID 0000.0001) //        EAST: Tile = 000003, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000002
(PID.TID 0000.0001) //       SOUTH: Tile = 000002, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000001
(PID.TID 0000.0001) //       NORTH: Tile = 000002, Process = 000000, Comm = put
(PID.TID 0000.0001) //                bi = 000002, bj = 000001
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  INI_PARMS: opening model parameter file "data"
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># ====================
(PID.TID 0000.0001) ># | Model parameters |
(PID.TID 0000.0001) ># ====================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># Continuous equation parameters
(PID.TID 0000.0001) > &PARM01
(PID.TID 0000.0001) > tRef = 15*20.,
(PID.TID 0000.0001) > sRef = 15*35.,
(PID.TID 0000.0001) > viscAr=1.E-3,
(PID.TID 0000.0001) > viscAh=5.E5,
(PID.TID 0000.0001) > diffKhT=0.,
(PID.TID 0000.0001) > diffKrT=3.E-5,
(PID.TID 0000.0001) > diffKhS=0.,
(PID.TID 0000.0001) > diffKrS=3.E-5,
(PID.TID 0000.0001) > rhonil=1035.,
(PID.TID 0000.0001) > rhoConstFresh=1000.,
(PID.TID 0000.0001) > eosType = 'JMD95Z',
(PID.TID 0000.0001) > ivdc_kappa=100.,
(PID.TID 0000.0001) > implicitDiffusion=.TRUE.,
(PID.TID 0000.0001) > allowFreezing=.TRUE.,
(PID.TID 0000.0001) > exactConserv=.TRUE.,
(PID.TID 0000.0001) > useRealFreshWaterFlux=.TRUE.,
(PID.TID 0000.0001) > useCDscheme=.TRUE.,
(PID.TID 0000.0001) ># turn on looped cells
(PID.TID 0000.0001) > hFacMin=.05,
(PID.TID 0000.0001) > hFacMindr=50.,
(PID.TID 0000.0001) ># set precision of data files
(PID.TID 0000.0001) > readBinaryPrec=32,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Elliptic solver parameters
(PID.TID 0000.0001) > &PARM02
(PID.TID 0000.0001) > cg2dMaxIters=500,
(PID.TID 0000.0001) > cg2dTargetResidual=1.E-19,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Time stepping parameters
(PID.TID 0000.0001) > &PARM03
(PID.TID 0000.0001) > nIter0=       0,
(PID.TID 0000.0001) >#nTimeSteps = 360,
(PID.TID 0000.0001) > nTimeSteps = 10,
(PID.TID 0000.0001) > deltaTmom = 1800.,
(PID.TID 0000.0001) > tauCD =     321428.,
(PID.TID 0000.0001) > deltaTtracer= 86400.,
(PID.TID 0000.0001) > deltaTClock = 86400.,
(PID.TID 0000.0001) > deltaTfreesurf= 86400.,
(PID.TID 0000.0001) > abEps = 0.1,
(PID.TID 0000.0001) > pChkptFreq= 2592000.,
(PID.TID 0000.0001) > dumpFreq=   2592000.,
(PID.TID 0000.0001) ># taveFreq=   2592000.,
(PID.TID 0000.0001) > monitorFreq=0.,
(PID.TID 0000.0001) ># 2 months restoring timescale for temperature
(PID.TID 0000.0001) >#tauThetaClimRelax=  5184000.,
(PID.TID 0000.0001) ># 6 months restoring timescale for salinity
(PID.TID 0000.0001) > tauSaltClimRelax = 15552000.,
(PID.TID 0000.0001) > periodicExternalForcing=.TRUE.,
(PID.TID 0000.0001) > externForcingPeriod=2592000.,
(PID.TID 0000.0001) > externForcingCycle=31104000.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Gridding parameters
(PID.TID 0000.0001) > &PARM04
(PID.TID 0000.0001) > usingSphericalPolarGrid=.TRUE.,
(PID.TID 0000.0001) > delR= 50.,  70., 100., 140., 190.,
(PID.TID 0000.0001) >      240., 290., 340., 390., 440.,
(PID.TID 0000.0001) >      490., 540., 590., 640., 690.,
(PID.TID 0000.0001) > ygOrigin=-80.,
(PID.TID 0000.0001) > dySpacing=4.,
(PID.TID 0000.0001) > dxSpacing=4.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Input datasets
(PID.TID 0000.0001) > &PARM05
(PID.TID 0000.0001) > bathyFile=      'bathymetry.bin',
(PID.TID 0000.0001) > hydrogThetaFile='lev_t.bin',
(PID.TID 0000.0001) > hydrogSaltFile= 'lev_s.bin',
(PID.TID 0000.0001) > zonalWindFile=  'trenberth_taux.bin',
(PID.TID 0000.0001) > meridWindFile=  'trenberth_tauy.bin',
(PID.TID 0000.0001) > thetaClimFile=  'lev_sst.bin',
(PID.TID 0000.0001) > saltClimFile=   'lev_sss.bin',
(PID.TID 0000.0001) > surfQFile=      'ncep_qnet.bin',
(PID.TID 0000.0001) ># fresh water flux is turned off, uncomment next line to turn on
(PID.TID 0000.0001) ># (not recommened together with surface salinity restoring)
(PID.TID 0000.0001) ># EmPmRFile=      'ncep_emp.bin',
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  INI_PARMS ; starts to read PARM01
(PID.TID 0000.0001)  INI_PARMS ; read PARM01 : OK
(PID.TID 0000.0001) S/R INI_PARMS: No request for barotropic solver
(PID.TID 0000.0001) S/R INI_PARMS: => Use implicitFreeSurface as default
(PID.TID 0000.0001)  INI_PARMS ; starts to read PARM02
(PID.TID 0000.0001)  INI_PARMS ; read PARM02 : OK
(PID.TID 0000.0001)  INI_PARMS ; starts to read PARM03
(PID.TID 0000.0001)  INI_PARMS ; read PARM03 : OK
(PID.TID 0000.0001)  INI_PARMS ; starts to read PARM04
(PID.TID 0000.0001)  INI_PARMS ; read PARM04 : OK
(PID.TID 0000.0001)  INI_PARMS ; starts to read PARM05
(PID.TID 0000.0001)  INI_PARMS ; read PARM05 : OK
(PID.TID 0000.0001)  INI_PARMS: finished reading file "data"
(PID.TID 0000.0001)  PACKAGES_BOOT: opening data.pkg
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.pkg
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.pkg"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># ********
(PID.TID 0000.0001) ># Packages
(PID.TID 0000.0001) ># ********
(PID.TID 0000.0001) > &PACKAGES
(PID.TID 0000.0001) > useGMRedi    = .TRUE.,
(PID.TID 0000.0001) > 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/gmredi               compiled   and   used ( useGMRedi                = T )
 pkg/grdchk               compiled   and   used ( useGrdchk                = T )
 -------- pkgs without standard "usePKG" On/Off switch in "data.pkg":  --------
 pkg/generic_advdiff      compiled   and   used ( useGAD                   = T )
 pkg/mom_common           compiled   and   used ( momStepping              = T )
 pkg/mom_vecinv           compiled but not used ( +vectorInvariantMomentum = F )
 pkg/mom_fluxform         compiled   and   used ( & not vectorInvariantMom = T )
 pkg/cd_code              compiled   and   used ( useCDscheme              = T )
 pkg/debug                compiled but not used ( debugMode                = F )
 pkg/rw                   compiled   and   used
 pkg/mdsio                compiled   and   used
 pkg/autodiff             compiled   and   used
 pkg/openad               compiled   and   used
 pkg/cost                 compiled   and   used
 pkg/ctrl                 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_background_K    = 1.e+3,
(PID.TID 0000.0001) >  GM_taper_scheme    = 'dm95',
(PID.TID 0000.0001) >  GM_maxSlope        = 1.e-2,
(PID.TID 0000.0001) >  GM_Kmin_horiz      = 50.,
(PID.TID 0000.0001) >  GM_Scrit           = 4.e-3,
(PID.TID 0000.0001) >  GM_Sd              = 1.e-3,
(PID.TID 0000.0001) ># GM_Visbeck_alpha   = 0.,
(PID.TID 0000.0001) ># GM_Visbeck_length  = 2.e+5,
(PID.TID 0000.0001) ># GM_Visbeck_depth   = 1.e+3,
(PID.TID 0000.0001) ># GM_Visbeck_maxval_K= 2.5e+3,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) >
(PID.TID 0000.0001) 
(PID.TID 0000.0001)  GM_READPARMS: finished reading data.gmredi
(PID.TID 0000.0001) AUTODIFF_READPARMS: opening data.autodiff
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.autodiff
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.autodiff"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># =========================
(PID.TID 0000.0001) ># pkg AUTODIFF parameters :
(PID.TID 0000.0001) ># =========================
(PID.TID 0000.0001) >#  inAdExact :: get an exact adjoint (no approximation) (def=.True.)
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &AUTODIFF_PARM01
(PID.TID 0000.0001) ># inAdExact = .FALSE.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001) AUTODIFF_READPARMS: finished reading data.autodiff
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) // AUTODIFF parameters :
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) inAdExact = /* get an exact adjoint (no approximation) */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) 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) 
(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) > numiter=1,
(PID.TID 0000.0001) > nfunc=3,
(PID.TID 0000.0001) > fmin=5.74,
(PID.TID 0000.0001) > iprint=10,
(PID.TID 0000.0001) > nupdate=8,
(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) >#
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) ># ECCO controlvariables
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) > &CTRL_NML
(PID.TID 0000.0001) > xx_hfluxm_file      = 'xx_hfluxm',
(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) CTRL_READPARMS: finished reading data.ctrl
(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_temp_tut  = 1.,
(PID.TID 0000.0001) > mult_hflux_tut = 2.,
(PID.TID 0000.0001) > lastinterval= 864000.,
(PID.TID 0000.0001) >#lastinterval=31104000.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001) COST_READPARMS: finished reading data.cost
(PID.TID 0000.0001) S/R COST_WEIGHTS: Temperature weights loaded from: Err_levitus_15layer.bin
(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-1,
(PID.TID 0000.0001) ># iGloPos          = 76,
(PID.TID 0000.0001) ># jGloPos          = 27,
(PID.TID 0000.0001) ># kGloPos          = 1,
(PID.TID 0000.0001) > nbeg             = 1,
(PID.TID 0000.0001) > nstep            = 1,
(PID.TID 0000.0001) > nend             = 3,
(PID.TID 0000.0001) > grdchkvarindex   =24,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) 
(PID.TID 0000.0001) GRDCHK_READPARMS: finished reading data.grdchk
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Gradient check configuration  >>> START <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
(PID.TID 0000.0001)   eps:                              0.100E+00
(PID.TID 0000.0001)   First location:                           1
(PID.TID 0000.0001)   Last location:                            3
(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) SET_PARMS: done
(PID.TID 0000.0001) Enter INI_VERTICAL_GRID: setInterFDr=    T ; setCenterDr=    F
(PID.TID 0000.0001) GAD_INIT_FIXED: GAD_OlMinSize=  1  0  1
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) // GAD parameters :
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) tempAdvScheme =   /* Temp. Horiz.Advection scheme selector */
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempVertAdvScheme =   /* Temp. Vert. Advection scheme selector */
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempMultiDimAdvec =   /* use Muti-Dim Advec method for Temp */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempSOM_Advection = /* use 2nd Order Moment Advection for Temp */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) AdamsBashforthGt = /* apply Adams-Bashforth extrapolation on Gt */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) AdamsBashforth_T = /* apply Adams-Bashforth extrapolation on Temp */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltAdvScheme =   /* Salt. Horiz.advection scheme selector */
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltVertAdvScheme =   /* Salt. Vert. Advection scheme selector */
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltMultiDimAdvec =   /* use Muti-Dim Advec method for Salt */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) saltSOM_Advection = /* use 2nd Order Moment Advection for Salt */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) AdamsBashforthGs = /* apply Adams-Bashforth extrapolation on Gs */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) AdamsBashforth_S = /* apply Adams-Bashforth extrapolation on Salt */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) ctrl-wet 1:    nvarlength =         2315
(PID.TID 0000.0001) ctrl-wet 2: surface wet C =          669
(PID.TID 0000.0001) ctrl-wet 3: surface wet W =          640
(PID.TID 0000.0001) ctrl-wet 4: surface wet S =          610
(PID.TID 0000.0001) ctrl-wet 4a:surface wet V =            0
(PID.TID 0000.0001) ctrl-wet 5: 3D wet points =         8651
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     1           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     2           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     3           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     4           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     5           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     6           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     7           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     8           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =     9           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    10           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    11           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    12           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    13           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    14           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    15           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    16           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    17           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    18           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    19           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    20           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    21           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    22           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    23           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    24           1
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    25           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    26           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    27           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    28           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    29           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    30           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    31           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    32           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    33           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    34           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    35           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    36           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    37           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    38           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    39           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    40           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    41           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    42           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    43           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    44           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    45           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    46           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    47           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    48           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    49           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    50           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    51           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    52           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    53           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    54           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    55           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    56           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    57           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    58           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    59           0
(PID.TID 0000.0001) ctrl-wet 6: no recs for i =    60           0
(PID.TID 0000.0001) ctrl-wet 7: flux         17302
(PID.TID 0000.0001) ctrl-wet 8: atmos        17302
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl-wet 13: global nvarlength for Nr =   15        2315
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    1        2315        2149        2206           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    2        2315        2149        2206           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    3        2254        2102        2146           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    4        2215        2058        2104           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    5        2178        2027        2070           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    6        2142        1987        2029           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    7        2114        1959        2004           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    8        2076        1918        1959           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=    9        2048        1887        1925           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   10        1999        1831        1869           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   11        1948        1771        1808           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   12        1850        1653        1705           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   13        1655        1404        1458           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   14        1372        1118        1164           0
(PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k=   15         828         623         671           0
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl-wet -------------------------------------------------
(PID.TID 0000.0001) ctrl_init: no. of control variables:            1
(PID.TID 0000.0001) ctrl_init: control vector length:            2315
(PID.TID 0000.0001) INI_CG2D: CG2D normalisation factor =  6.5682677425711703E-05
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model configuration
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) //
(PID.TID 0000.0001) // "Physical" paramters ( PARM01 in namelist )
(PID.TID 0000.0001) //
(PID.TID 0000.0001) buoyancyRelation = /* Type of relation to get Buoyancy */
(PID.TID 0000.0001)               'OCEANIC'
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) fluidIsAir   =  /* fluid major constituent is Air */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) fluidIsWater =  /* fluid major constituent is Water */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingPCoords =  /* use p (or p*) vertical coordinate */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingZCoords =  /* use z (or z*) vertical coordinate */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tRef =   /* Reference temperature profile ( oC or K ) */
(PID.TID 0000.0001)    15 @  2.000000000000000E+01              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) sRef =   /* Reference salinity profile ( psu ) */
(PID.TID 0000.0001)    15 @  3.500000000000000E+01              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useStrainTensionVisc= /* Use StrainTension Form of Viscous Operator */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useVariableVisc = /* Use variable horizontal viscosity */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useHarmonicVisc = /* Use harmonic horizontal viscosity */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useBiharmonicVisc= /* Use biharmonic horiz.  viscosity */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useSmag3D = /* Use isotropic 3-D Smagorinsky viscosity */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) viscAh  =   /* Lateral harmonic viscosity ( m^2/s ) */
(PID.TID 0000.0001)                 5.000000000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) viscA4  =   /* Lateral biharmonic viscosity ( m^4/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) no_slip_sides =  /* Viscous BCs: No-slip sides */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) sideDragFactor = /* side-drag scaling factor (non-dim) */
(PID.TID 0000.0001)                 2.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) viscArNr = /* vertical profile of vertical viscosity ( m^2/s )*/
(PID.TID 0000.0001)    15 @  1.000000000000000E-03              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) no_slip_bottom =  /* Viscous BCs: No-slip bottom */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) 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) diffKhT =   /* Laplacian diffusion of heat laterally ( m^2/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffK4T =   /* Biharmonic diffusion of heat laterally ( m^4/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKhS =   /* Laplacian diffusion of salt laterally ( m^2/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffK4S =   /* Biharmonic diffusion of salt laterally ( m^4/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrNrT = /* vertical profile of vertical diffusion of Temp ( m^2/s )*/
(PID.TID 0000.0001)    15 @  3.000000000000000E-05              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrNrS = /* vertical profile of vertical diffusion of Salt ( m^2/s )*/
(PID.TID 0000.0001)    15 @  3.000000000000000E-05              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrBL79surf = /* Surface diffusion for Bryan and Lewis 79 ( m^2/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrBL79deep = /* Deep diffusion for Bryan and Lewis 1979 ( m^2/s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrBL79scl = /* Depth scale for Bryan and Lewis 1979 ( m ) */
(PID.TID 0000.0001)                 2.000000000000000E+02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) diffKrBL79Ho = /* Turning depth for Bryan and Lewis 1979 ( m ) */
(PID.TID 0000.0001)                -2.000000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) ivdc_kappa = /* Implicit Vertical Diffusivity for Convection ( m^2/s) */
(PID.TID 0000.0001)                 1.000000000000000E+02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hMixCriteria=  /* Criteria for mixed-layer diagnostic */
(PID.TID 0000.0001)                -8.000000000000000E-01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dRhoSmall =  /* Parameter for mixed-layer diagnostic */
(PID.TID 0000.0001)                 1.000000000000000E-06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hMixSmooth=  /* Smoothing parameter for mixed-layer diagnostic */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) eosType =  /* Type of Equation of State */
(PID.TID 0000.0001)               'JMD95Z'
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) celsius2K = /* 0 degree Celsius converted to Kelvin ( K ) */
(PID.TID 0000.0001)                 2.731500000000000E+02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoConst  = /* Reference density (Boussinesq)  ( kg/m^3 ) */
(PID.TID 0000.0001)                 1.035000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoFacC = /* normalized Reference density @ cell-Center (-) */
(PID.TID 0000.0001)    15 @  1.000000000000000E+00              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoFacF = /* normalized Reference density @ W-Interface (-) */
(PID.TID 0000.0001)    16 @  1.000000000000000E+00              /* K =  1: 16 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rhoConstFresh = /* Fresh-water reference density ( kg/m^3 ) */
(PID.TID 0000.0001)                 1.000000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) gravity =   /* Gravitational acceleration ( m/s^2 ) */
(PID.TID 0000.0001)                 9.810000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) gBaro =   /* Barotropic gravity ( m/s^2 ) */
(PID.TID 0000.0001)                 9.810000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) 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) hFacMin =   /* minimum partial cell factor (hFac) */
(PID.TID 0000.0001)                 5.000000000000000E-02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hFacMinDr = /* minimum partial cell thickness ( m) */
(PID.TID 0000.0001)                 5.000000000000000E+01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) exactConserv =  /* Exact Volume Conservation on/off flag*/
(PID.TID 0000.0001)                   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)                       0
(PID.TID 0000.0001)      -1,0= Off ; 1,2,3= On, 2=+rescale gU,gV, 3=+update cg2d solv.
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hFacInf =   /* lower threshold for hFac (nonlinFreeSurf only)*/
(PID.TID 0000.0001)                 2.000000000000000E-01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) hFacSup =   /* upper threshold for hFac (nonlinFreeSurf only)*/
(PID.TID 0000.0001)                 2.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) select_rStar = /* r* Vertical coord. options (=0 r coord.; >0 uses r*)*/
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useRealFreshWaterFlux = /* Real Fresh Water Flux on/off flag*/
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) temp_EvPrRn = /* Temp. of Evap/Prec/R (UNSET=use local T)(oC)*/
(PID.TID 0000.0001)                 1.234567000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) salt_EvPrRn = /* Salin. of Evap/Prec/R (UNSET=use local S)(psu)*/
(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)(psu)*/
(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)(psu)*/
(PID.TID 0000.0001)                -1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) use3Dsolver = /* use 3-D pressure solver on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nonHydrostatic =  /* Non-Hydrostatic on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nh_Am2 = /* Non-Hydrostatic terms scaling factor */
(PID.TID 0000.0001)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) implicitNHPress = /* Non-Hyd Pressure implicit factor (0-1)*/
(PID.TID 0000.0001)                 1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectNHfreeSurf = /* Non-Hyd (free-)Surface option */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) quasiHydrostatic = /* Quasi-Hydrostatic on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) calc_wVelocity = /* vertical velocity calculation on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momStepping =  /* Momentum equation on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) vectorInvariantMomentum= /* Vector-Invariant Momentum on/off */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momAdvection =  /* Momentum advection on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momViscosity =  /* Momentum viscosity on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momImplVertAdv= /* Momentum implicit vert. advection on/off*/
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) implicitViscosity = /* Implicit viscosity on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) metricTerms =  /* metric-Terms on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useNHMTerms = /* Non-Hydrostatic Metric-Terms on/off */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectCoriMap = /* Coriolis Map options (0,1,2,3)*/
(PID.TID 0000.0001)                       2
(PID.TID 0000.0001)     0= f-Plane ; 1= Beta-Plane ; 2= Spherical ; 3= read from file
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) use3dCoriolis = /* 3-D Coriolis on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useCoriolis =  /* Coriolis on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useCDscheme =  /* CD scheme on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useEnergyConservingCoriolis= /* Flx-Form Coriolis scheme flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useJamartWetPoints= /* Coriolis WetPoints method flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useJamartMomAdv= /* V.I Non-linear terms Jamart flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useAbsVorticity= /* V.I Works with f+zeta in Coriolis */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectVortScheme= /* V.I Scheme selector for Vorticity-Term */
(PID.TID 0000.0001)               123456789
(PID.TID 0000.0001)    = 0 : enstrophy (Shallow-Water Eq.) conserving scheme by Sadourny, JAS 75
(PID.TID 0000.0001)    = 1 : same as 0 with modified hFac
(PID.TID 0000.0001)    = 2 : energy conserving scheme (used by Sadourny in JAS 75 paper)
(PID.TID 0000.0001)    = 3 : energy (general) and enstrophy (2D, nonDiv.) conserving scheme
(PID.TID 0000.0001)          from Sadourny (Burridge & Haseler, ECMWF Rep.4, 1977)
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) upwindVorticity= /* V.I Upwind bias vorticity flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) highOrderVorticity= /* V.I High order vort. advect. flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) upwindShear= /* V.I Upwind vertical Shear advection flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectKEscheme= /* V.I Kinetic Energy scheme selector */
(PID.TID 0000.0001)                       0
(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) momPressureForcing =  /* Momentum pressure term on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) implicitIntGravWave= /* Implicit Internal Gravity Wave flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) staggerTimeStep =    /* Stagger time stepping on/off flag */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) doResetHFactors = /* reset thickness factors @ each time-step */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) multiDimAdvection =  /* enable/disable Multi-Dim Advection */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useMultiDimAdvec =   /* Multi-Dim Advection is/is-not used */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) implicitDiffusion = /* Implicit Diffusion on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempStepping =  /* Temperature equation on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempAdvection = /* Temperature advection on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempImplVertAdv = /* Temp. implicit vert. advection on/off */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tempForcing  =  /* Temperature forcing on/off flag */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) doThetaClimRelax = /* apply SST relaxation on/off flag */
(PID.TID 0000.0001)                   F
(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)  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) //
(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)                     500
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cg2dChkResFreq =   /* 2d con. grad convergence test frequency */
(PID.TID 0000.0001)                       1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cg2dUseMinResSol= /* use cg2d last-iter(=0) / min-resid.(=1) solution */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cg2dTargetResidual =   /* 2d con. grad target residual  */
(PID.TID 0000.0001)                 1.000000000000000E-19
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cg2dTargetResWunit =   /* CG2d target residual [W units] */
(PID.TID 0000.0001)                -1.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) cg2dPreCondFreq =   /* Freq. for updating cg2d preconditioner */
(PID.TID 0000.0001)                       1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) useSRCGSolver =  /* use single reduction CG solver(s) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) printResidualFreq = /* Freq. for printing CG residual */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) //
(PID.TID 0000.0001) // Time stepping paramters ( PARM03 in namelist )
(PID.TID 0000.0001) //
(PID.TID 0000.0001) deltaTMom =   /* Momentum equation timestep ( s ) */
(PID.TID 0000.0001)                 1.800000000000000E+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)    15 @  8.640000000000000E+04              /* K =  1: 15 */
(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)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momForcingOutAB = /* =1: take Momentum Forcing out of Adams-Bash. stepping */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tracForcingOutAB = /* =1: take T,S,pTr Forcing out of Adams-Bash. stepping */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) momDissip_In_AB = /* put Dissipation Tendency in Adams-Bash. stepping */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) doAB_onGtGs = /* apply AB on Tendencies (rather than on T,S)*/
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) abEps =   /* Adams-Bashforth-2 stabilizing weight */
(PID.TID 0000.0001)                 1.000000000000000E-01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) 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.943999900444267E-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)                      10
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nEndIter =   /* Run ending timestep number */
(PID.TID 0000.0001)                      10
(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)                 8.640000000000000E+05
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) pChkPtFreq = /* Permanent restart/pickup file interval ( s ) */
(PID.TID 0000.0001)                 2.592000000000000E+06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) chkPtFreq  = /* Rolling restart/pickup file interval ( s ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) pickup_write_mdsio =   /* Model IO flag. */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) pickup_read_mdsio =   /* Model IO flag. */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) pickup_write_immed =   /* Model IO flag. */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) writePickupAtEnd =   /* Model IO flag. */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dumpFreq =   /* Model state write out interval ( s ). */
(PID.TID 0000.0001)                 2.592000000000000E+06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dumpInitAndLast= /* write out Initial & Last iter. model state */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) snapshot_mdsio =   /* Model IO flag. */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) monitorFreq =   /* Monitor output interval ( s ). */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) monitorSelect = /* select group of variables to monitor */
(PID.TID 0000.0001)                       3
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) monitor_stdio =   /* Model IO flag. */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) externForcingPeriod =   /* forcing period (s) */
(PID.TID 0000.0001)                 2.592000000000000E+06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) externForcingCycle =   /* period of the cyle (s). */
(PID.TID 0000.0001)                 3.110400000000000E+07
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tauThetaClimRelax =   /* relaxation time scale (s) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) tauSaltClimRelax =   /* relaxation time scale (s) */
(PID.TID 0000.0001)                 1.555200000000000E+07
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) latBandClimRelax =   /* max. Lat. where relaxation */
(PID.TID 0000.0001)                 1.800000000000000E+02
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) //
(PID.TID 0000.0001) // Gridding paramters ( PARM04 in namelist )
(PID.TID 0000.0001) //
(PID.TID 0000.0001) usingCartesianGrid = /* Cartesian coordinates flag ( True/False ) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingCylindricalGrid = /* Cylindrical coordinates flag ( True/False ) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingSphericalPolarGrid = /* Spherical coordinates flag ( True/False ) */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingCurvilinearGrid = /* Curvilinear coordinates flag ( True/False ) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) selectSigmaCoord = /* Hybrid-Sigma Vert. Coordinate option */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) Ro_SeaLevel = /* r(1) ( units of r ==  m ) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(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) mass2rUnit = /* convert mass per unit area [kg/m2] to r-units [m] */
(PID.TID 0000.0001)                 9.661835748792270E-04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rUnit2mass = /* convert r-units [m] to mass per unit area [kg/m2] */
(PID.TID 0000.0001)                 1.035000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) drC =   /* C spacing ( units of r ) */
(PID.TID 0000.0001)                 2.500000000000000E+01,      /* K =  1 */
(PID.TID 0000.0001)                 6.000000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                 8.500000000000000E+01,      /* K =  3 */
(PID.TID 0000.0001)                 1.200000000000000E+02,      /* K =  4 */
(PID.TID 0000.0001)                 1.650000000000000E+02,      /* K =  5 */
(PID.TID 0000.0001)                 2.150000000000000E+02,      /* K =  6 */
(PID.TID 0000.0001)                 2.650000000000000E+02,      /* K =  7 */
(PID.TID 0000.0001)                 3.150000000000000E+02,      /* K =  8 */
(PID.TID 0000.0001)                 3.650000000000000E+02,      /* K =  9 */
(PID.TID 0000.0001)                 4.150000000000000E+02,      /* K = 10 */
(PID.TID 0000.0001)                 4.650000000000000E+02,      /* K = 11 */
(PID.TID 0000.0001)                 5.150000000000000E+02,      /* K = 12 */
(PID.TID 0000.0001)                 5.650000000000000E+02,      /* K = 13 */
(PID.TID 0000.0001)                 6.150000000000000E+02,      /* K = 14 */
(PID.TID 0000.0001)                 6.650000000000000E+02,      /* K = 15 */
(PID.TID 0000.0001)                 3.450000000000000E+02       /* K = 16 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) drF =   /* W spacing ( units of r ) */
(PID.TID 0000.0001)                 5.000000000000000E+01,      /* K =  1 */
(PID.TID 0000.0001)                 7.000000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                 1.000000000000000E+02,      /* K =  3 */
(PID.TID 0000.0001)                 1.400000000000000E+02,      /* K =  4 */
(PID.TID 0000.0001)                 1.900000000000000E+02,      /* K =  5 */
(PID.TID 0000.0001)                 2.400000000000000E+02,      /* K =  6 */
(PID.TID 0000.0001)                 2.900000000000000E+02,      /* K =  7 */
(PID.TID 0000.0001)                 3.400000000000000E+02,      /* K =  8 */
(PID.TID 0000.0001)                 3.900000000000000E+02,      /* K =  9 */
(PID.TID 0000.0001)                 4.400000000000000E+02,      /* K = 10 */
(PID.TID 0000.0001)                 4.900000000000000E+02,      /* K = 11 */
(PID.TID 0000.0001)                 5.400000000000000E+02,      /* K = 12 */
(PID.TID 0000.0001)                 5.900000000000000E+02,      /* K = 13 */
(PID.TID 0000.0001)                 6.400000000000000E+02,      /* K = 14 */
(PID.TID 0000.0001)                 6.900000000000000E+02       /* K = 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) delX = /* U spacing ( m - cartesian, degrees - spherical ) */
(PID.TID 0000.0001)    90 @  4.000000000000000E+00              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) delY = /* V spacing ( m - cartesian, degrees - spherical ) */
(PID.TID 0000.0001)    40 @  4.000000000000000E+00              /* J =  1: 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) xgOrigin = /* X-axis origin of West  edge (cartesian: m, lat-lon: deg) */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) ygOrigin = /* Y-axis origin of South edge (cartesian: m, lat-lon: deg) */
(PID.TID 0000.0001)                -8.000000000000000E+01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rSphere =  /* Radius ( ignored - cartesian, m - spherical ) */
(PID.TID 0000.0001)                 6.370000000000000E+06
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deepAtmosphere = /* Deep/Shallow Atmosphere flag (True/False) */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) xC =  /* xC(:,1,:,1) : P-point X coord ( deg. or m if cartesian) */
(PID.TID 0000.0001)                 2.000000000000000E+00,      /* I =  1 */
(PID.TID 0000.0001)                 6.000000000000000E+00,      /* I =  2 */
(PID.TID 0000.0001)                 1.000000000000000E+01,      /* I =  3 */
(PID.TID 0000.0001)      . . .
(PID.TID 0000.0001)                 8.200000000000000E+01,      /* I = 21 */
(PID.TID 0000.0001)                 8.600000000000000E+01,      /* I = 22 */
(PID.TID 0000.0001)                 9.000000000000000E+01,      /* I = 23 */
(PID.TID 0000.0001)                 9.400000000000000E+01,      /* I = 24 */
(PID.TID 0000.0001)                 9.800000000000000E+01,      /* I = 25 */
(PID.TID 0000.0001)                 1.020000000000000E+02,      /* I = 26 */
(PID.TID 0000.0001)      . . .
(PID.TID 0000.0001)                 1.700000000000000E+02,      /* I = 43 */
(PID.TID 0000.0001)                 1.740000000000000E+02,      /* I = 44 */
(PID.TID 0000.0001)                 1.780000000000000E+02,      /* I = 45 */
(PID.TID 0000.0001)                 1.820000000000000E+02,      /* I = 46 */
(PID.TID 0000.0001)                 1.860000000000000E+02,      /* I = 47 */
(PID.TID 0000.0001)                 1.900000000000000E+02,      /* I = 48 */
(PID.TID 0000.0001)      . . .
(PID.TID 0000.0001)                 2.580000000000000E+02,      /* I = 65 */
(PID.TID 0000.0001)                 2.620000000000000E+02,      /* I = 66 */
(PID.TID 0000.0001)                 2.660000000000000E+02,      /* I = 67 */
(PID.TID 0000.0001)                 2.700000000000000E+02,      /* I = 68 */
(PID.TID 0000.0001)                 2.740000000000000E+02,      /* I = 69 */
(PID.TID 0000.0001)                 2.780000000000000E+02,      /* I = 70 */
(PID.TID 0000.0001)      . . .
(PID.TID 0000.0001)                 3.500000000000000E+02,      /* I = 88 */
(PID.TID 0000.0001)                 3.540000000000000E+02,      /* I = 89 */
(PID.TID 0000.0001)                 3.580000000000000E+02       /* I = 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) yC =  /* yC(1,:,1,:) : P-point Y coord ( deg. or m if cartesian) */
(PID.TID 0000.0001)                -7.800000000000000E+01,      /* J =  1 */
(PID.TID 0000.0001)                -7.400000000000000E+01,      /* J =  2 */
(PID.TID 0000.0001)                -7.000000000000000E+01,      /* J =  3 */
(PID.TID 0000.0001)                -6.600000000000000E+01,      /* J =  4 */
(PID.TID 0000.0001)                -6.200000000000000E+01,      /* J =  5 */
(PID.TID 0000.0001)                -5.800000000000000E+01,      /* J =  6 */
(PID.TID 0000.0001)                -5.400000000000000E+01,      /* J =  7 */
(PID.TID 0000.0001)                -5.000000000000000E+01,      /* J =  8 */
(PID.TID 0000.0001)                -4.600000000000000E+01,      /* J =  9 */
(PID.TID 0000.0001)                -4.200000000000000E+01,      /* J = 10 */
(PID.TID 0000.0001)                -3.800000000000000E+01,      /* J = 11 */
(PID.TID 0000.0001)                -3.400000000000000E+01,      /* J = 12 */
(PID.TID 0000.0001)                -3.000000000000000E+01,      /* J = 13 */
(PID.TID 0000.0001)                -2.600000000000000E+01,      /* J = 14 */
(PID.TID 0000.0001)                -2.200000000000000E+01,      /* J = 15 */
(PID.TID 0000.0001)                -1.800000000000000E+01,      /* J = 16 */
(PID.TID 0000.0001)                -1.400000000000000E+01,      /* J = 17 */
(PID.TID 0000.0001)                -1.000000000000000E+01,      /* J = 18 */
(PID.TID 0000.0001)                -6.000000000000000E+00,      /* J = 19 */
(PID.TID 0000.0001)                -2.000000000000000E+00,      /* J = 20 */
(PID.TID 0000.0001)                 2.000000000000000E+00,      /* J = 21 */
(PID.TID 0000.0001)                 6.000000000000000E+00,      /* J = 22 */
(PID.TID 0000.0001)                 1.000000000000000E+01,      /* J = 23 */
(PID.TID 0000.0001)                 1.400000000000000E+01,      /* J = 24 */
(PID.TID 0000.0001)                 1.800000000000000E+01,      /* J = 25 */
(PID.TID 0000.0001)                 2.200000000000000E+01,      /* J = 26 */
(PID.TID 0000.0001)                 2.600000000000000E+01,      /* J = 27 */
(PID.TID 0000.0001)                 3.000000000000000E+01,      /* J = 28 */
(PID.TID 0000.0001)                 3.400000000000000E+01,      /* J = 29 */
(PID.TID 0000.0001)                 3.800000000000000E+01,      /* J = 30 */
(PID.TID 0000.0001)                 4.200000000000000E+01,      /* J = 31 */
(PID.TID 0000.0001)                 4.600000000000000E+01,      /* J = 32 */
(PID.TID 0000.0001)                 5.000000000000000E+01,      /* J = 33 */
(PID.TID 0000.0001)                 5.400000000000000E+01,      /* J = 34 */
(PID.TID 0000.0001)                 5.800000000000000E+01,      /* J = 35 */
(PID.TID 0000.0001)                 6.200000000000000E+01,      /* J = 36 */
(PID.TID 0000.0001)                 6.600000000000000E+01,      /* J = 37 */
(PID.TID 0000.0001)                 7.000000000000000E+01,      /* J = 38 */
(PID.TID 0000.0001)                 7.400000000000000E+01,      /* J = 39 */
(PID.TID 0000.0001)                 7.800000000000000E+01       /* J = 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rcoord = /* P-point R coordinate (  units of r ) */
(PID.TID 0000.0001)                -2.500000000000000E+01,      /* K =  1 */
(PID.TID 0000.0001)                -8.500000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                -1.700000000000000E+02,      /* K =  3 */
(PID.TID 0000.0001)                -2.900000000000000E+02,      /* K =  4 */
(PID.TID 0000.0001)                -4.550000000000000E+02,      /* K =  5 */
(PID.TID 0000.0001)                -6.700000000000000E+02,      /* K =  6 */
(PID.TID 0000.0001)                -9.350000000000000E+02,      /* K =  7 */
(PID.TID 0000.0001)                -1.250000000000000E+03,      /* K =  8 */
(PID.TID 0000.0001)                -1.615000000000000E+03,      /* K =  9 */
(PID.TID 0000.0001)                -2.030000000000000E+03,      /* K = 10 */
(PID.TID 0000.0001)                -2.495000000000000E+03,      /* K = 11 */
(PID.TID 0000.0001)                -3.010000000000000E+03,      /* K = 12 */
(PID.TID 0000.0001)                -3.575000000000000E+03,      /* K = 13 */
(PID.TID 0000.0001)                -4.190000000000000E+03,      /* K = 14 */
(PID.TID 0000.0001)                -4.855000000000000E+03       /* K = 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rF =   /* W-Interf. R coordinate (  units of r ) */
(PID.TID 0000.0001)                 0.000000000000000E+00,      /* K =  1 */
(PID.TID 0000.0001)                -5.000000000000000E+01,      /* K =  2 */
(PID.TID 0000.0001)                -1.200000000000000E+02,      /* K =  3 */
(PID.TID 0000.0001)                -2.200000000000000E+02,      /* K =  4 */
(PID.TID 0000.0001)                -3.600000000000000E+02,      /* K =  5 */
(PID.TID 0000.0001)                -5.500000000000000E+02,      /* K =  6 */
(PID.TID 0000.0001)                -7.900000000000000E+02,      /* K =  7 */
(PID.TID 0000.0001)                -1.080000000000000E+03,      /* K =  8 */
(PID.TID 0000.0001)                -1.420000000000000E+03,      /* K =  9 */
(PID.TID 0000.0001)                -1.810000000000000E+03,      /* K = 10 */
(PID.TID 0000.0001)                -2.250000000000000E+03,      /* K = 11 */
(PID.TID 0000.0001)                -2.740000000000000E+03,      /* K = 12 */
(PID.TID 0000.0001)                -3.280000000000000E+03,      /* K = 13 */
(PID.TID 0000.0001)                -3.870000000000000E+03,      /* K = 14 */
(PID.TID 0000.0001)                -4.510000000000000E+03,      /* K = 15 */
(PID.TID 0000.0001)                -5.200000000000000E+03       /* K = 16 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deepFacC = /* deep-model grid factor @ cell-Center (-) */
(PID.TID 0000.0001)    15 @  1.000000000000000E+00              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deepFacF = /* deep-model grid factor @ W-Interface (-) */
(PID.TID 0000.0001)    16 @  1.000000000000000E+00              /* K =  1: 16 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rVel2wUnit = /* convert units: rVel -> wSpeed (=1 if z-coord)*/
(PID.TID 0000.0001)    16 @  1.000000000000000E+00              /* K =  1: 16 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) wUnit2rVel = /* convert units: wSpeed -> rVel (=1 if z-coord)*/
(PID.TID 0000.0001)    16 @  1.000000000000000E+00              /* K =  1: 16 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) dBdrRef = /* Vertical grad. of reference buoyancy [(m/s/r)^2] */
(PID.TID 0000.0001)    15 @  0.000000000000000E+00              /* K =  1: 15 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) 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.450087291636666E+10,      /* J =  2 */
(PID.TID 0000.0001)                 6.762652439100235E+10,      /* J =  3 */
(PID.TID 0000.0001)                 8.042270623659996E+10,      /* J =  4 */
(PID.TID 0000.0001)                 9.282707674147525E+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.373525277677230E+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.712363708253574E+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.947228162702493E+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.947228162702493E+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.712363708253574E+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.373525277677230E+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.282707674147525E+10,      /* J = 36 */
(PID.TID 0000.0001)                 8.042270623659996E+10,      /* J = 37 */
(PID.TID 0000.0001)                 6.762652439100235E+10,      /* J = 38 */
(PID.TID 0000.0001)                 5.450087291636666E+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.450087291636666E+10,      /* J =  2 */
(PID.TID 0000.0001)                 6.762652439100235E+10,      /* J =  3 */
(PID.TID 0000.0001)                 8.042270623659996E+10,      /* J =  4 */
(PID.TID 0000.0001)                 9.282707674147525E+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.373525277677230E+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.712363708253574E+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.947228162702493E+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.947228162702493E+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.712363708253574E+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.373525277677230E+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.282707674147525E+10,      /* J = 36 */
(PID.TID 0000.0001)                 8.042270623659996E+10,      /* J = 37 */
(PID.TID 0000.0001)                 6.762652439100235E+10,      /* J = 38 */
(PID.TID 0000.0001)                 5.450087291636666E+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.433488626798251E+10              /* I =  1: 90 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) rAs =  /* rAs(1,:,1,:) ( units: m^2 ) */
(PID.TID 0000.0001)                 3.433488626798251E+10,      /* J =  1 */
(PID.TID 0000.0001)                 4.783442523123625E+10,      /* J =  2 */
(PID.TID 0000.0001)                 6.110091968306415E+10,      /* J =  3 */
(PID.TID 0000.0001)                 7.406973659603818E+10,      /* J =  4 */
(PID.TID 0000.0001)                 8.667769319778081E+10,      /* J =  5 */
(PID.TID 0000.0001)                 9.886336479107465E+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.806323556642997E+11,      /* J = 15 */
(PID.TID 0000.0001)                 1.858023487204767E+11,      /* J = 16 */
(PID.TID 0000.0001)                 1.900671314104744E+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.977267295821496E+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.900671314104744E+11,      /* J = 25 */
(PID.TID 0000.0001)                 1.858023487204767E+11,      /* J = 26 */
(PID.TID 0000.0001)                 1.806323556642997E+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.886336479107465E+10,      /* J = 36 */
(PID.TID 0000.0001)                 8.667769319778081E+10,      /* J = 37 */
OAD: TIMING: stamp 0: 1395243753.634389
OAD: TIMING: stamp 1: 1395243762.189365
OAD: TIMING: delta stamps 1-0: 8.554976
(PID.TID 0000.0001)                 7.406973659603818E+10,      /* J = 38 */
(PID.TID 0000.0001)                 6.110091968306415E+10,      /* J = 39 */
(PID.TID 0000.0001)                 4.783442523123625E+10       /* J = 40 */
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) globalArea = /* Integrated horizontal Area (m^2) */
(PID.TID 0000.0001)                 3.450614146649749E+14
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // End of Model config. summary
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
(PID.TID 0000.0001) == Packages configuration : Check & print summary ==
(PID.TID 0000.0001) 
(PID.TID 0000.0001) GMREDI_CHECK: #define GMREDI
(PID.TID 0000.0001) GM_AdvForm =     /* if FALSE => use SkewFlux Form */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_InMomAsStress = /* if TRUE => apply as Eddy Stress */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_AdvSeparate = /* Calc Bolus & Euler Adv. separately */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_ExtraDiag =   /* Tensor Extra Diag (line 1&2) non 0 */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_isopycK =    /* Background Isopyc. Diffusivity [m^2/s] */
(PID.TID 0000.0001)                 1.000000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_skewflx*K =  /* Background GM_SkewFlx Diffusivity [m^2/s] */
(PID.TID 0000.0001)                 1.000000000000000E+03
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_advec*K =    /* Backg. GM-Advec(=Bolus) Diffusivity [m^2/s]*/
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_Kmin_horiz = /* Minimum Horizontal Diffusivity [m^2/s] */
(PID.TID 0000.0001)                 5.000000000000000E+01
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_Visbeck_alpha = /* Visbeck alpha coeff. [-] */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) GM_Small_Number =  /* epsilon used in slope calc */
(PID.TID 0000.0001)                 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) CTRL_CHECK: ctrl package
(PID.TID 0000.0001) COST_CHECK: cost package
(PID.TID 0000.0001) GRDCHK_CHECK: grdchk package
(PID.TID 0000.0001) GAD_CHECK: #define ALLOW_GENERIC_ADVDIFF
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Check Model config. (CONFIG_CHECK):
(PID.TID 0000.0001) // CONFIG_CHECK : Normal End
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 OAD: IT+      1048576
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -4.66293670342566E-15  2.54952604282385E+00
 cg2d: Sum(rhs),rhsMax =   3.94129173741931E-15  2.62363057992942E+00
 cg2d: Sum(rhs),rhsMax =  -1.09912079437890E-14  2.39895957779123E+00
 cg2d: Sum(rhs),rhsMax =   8.88178419700125E-16  2.21798267529300E+00
 cg2d: Sum(rhs),rhsMax =  -5.32907051820075E-15  2.01662326576520E+00
 cg2d: Sum(rhs),rhsMax =  -5.77315972805081E-14  1.79933840571731E+00
 cg2d: Sum(rhs),rhsMax =  -3.46389583683049E-14  1.57603319715841E+00
 cg2d: Sum(rhs),rhsMax =  -1.99840144432528E-14  1.36017289258806E+00
 cg2d: Sum(rhs),rhsMax =   2.87547763377916E-14  1.18689213330771E+00
 OAD: DT+      1048576
 OAD: IT+     17825792
 OAD: IT+     34603008
 OAD: DT+     17825792
 OAD: IT+     51380224
 OAD: IT+     68157440
 cg2d: Sum(rhs),rhsMax =  -3.09752223870419E-14  1.12112051380257E+00
(PID.TID 0000.0001) %CHECKPOINT        10 ckptA
  early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.141635681714548D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.185615632957278D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.140673784968176D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.152098128542327D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
  local fc =  0.620023228182329D+01
 global fc =  0.620023228182329D+01
 cg2d: Sum(rhs),rhsMax =   0.00000000000000E+00  0.00000000000000E+00
 cg2d: Sum(rhs),rhsMax =  -4.51860771022439E-14  1.18689213330773E+00
 cg2d: Sum(rhs),rhsMax =  -2.22044604925031E-15  3.97560864970023E-08
 cg2d: Sum(rhs),rhsMax =   2.26485497023532E-14  1.36017289258806E+00
 cg2d: Sum(rhs),rhsMax =  -4.44089209850063E-16  3.41865552528313E-08
 cg2d: Sum(rhs),rhsMax =   2.26485497023532E-14  1.57603319715842E+00
 cg2d: Sum(rhs),rhsMax =   1.33226762955019E-15  4.24929829971601E-08
OAD: TIMING: stamp 2: 1395243778.892145
OAD: TIMING: delta stamps 2-1: 16.702780
OAD: TIMING: delta stamps 2-0: 25.257756
OAD: TIMING: ratio stamps (2-1)/(1-0): 1.670278e+07/8.554976e+06=1.952405e+00
 cg2d: Sum(rhs),rhsMax =  -2.26485497023532E-14  1.79933840571732E+00
 cg2d: Sum(rhs),rhsMax =  -2.22044604925031E-15  8.44536286804658E-08
 cg2d: Sum(rhs),rhsMax =   0.00000000000000E+00  2.01662326576520E+00
 cg2d: Sum(rhs),rhsMax =  -3.33066907387547E-16  1.36480754116999E-07
 cg2d: Sum(rhs),rhsMax =   1.06581410364015E-14  2.21798267529300E+00
 cg2d: Sum(rhs),rhsMax =  -2.38697950294409E-15  2.12919271877336E-07
 cg2d: Sum(rhs),rhsMax =   1.93178806284777E-14  2.39895957779122E+00
 cg2d: Sum(rhs),rhsMax =   6.52256026967279E-16  3.15770102898000E-07
 cg2d: Sum(rhs),rhsMax =  -5.71764857681956E-15  2.62363057992944E+00
 cg2d: Sum(rhs),rhsMax =   1.33226762955019E-15  4.35249443192189E-07
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -4.44089209850063E-15  2.54952604282385E+00
 cg2d: Sum(rhs),rhsMax =   3.10862446895044E-15  5.76014300791293E-07
 ph-pack: packing ecco_cost
 ph-pack: packing ecco_ctrl
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Gradient-check starts (grdchk_main)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) grdchk reference fc: fcref       =  6.20023228182329E+00
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
(PID.TID 0000.0001) ====== Starts gradient-check number   1 (=ichknum) =======
 ph-test icomp, ncvarcomp, ichknum            1        2315           1
 ph-grd _loc: bi, bj, icomptest, ichknum            1           1           0           1
 ph-grd -->hit<--           43           2           1           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   43    2    1 ; 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) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    1 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -4.66293670342566E-15  2.54952604282385E+00
 cg2d: Sum(rhs),rhsMax =   3.94129173741931E-15  2.62363057992942E+00
 cg2d: Sum(rhs),rhsMax =   1.35447209004269E-14  2.39895957779137E+00
 cg2d: Sum(rhs),rhsMax =   1.33226762955019E-14  2.21798267529369E+00
 cg2d: Sum(rhs),rhsMax =  -1.90958360235527E-14  2.01662326576458E+00
 cg2d: Sum(rhs),rhsMax =  -3.55271367880050E-14  1.79933840571470E+00
 cg2d: Sum(rhs),rhsMax =   3.55271367880050E-15  1.57603319715180E+00
 cg2d: Sum(rhs),rhsMax =   1.86517468137026E-14  1.36017289257808E+00
 cg2d: Sum(rhs),rhsMax =  -1.13242748511766E-14  1.18689213332293E+00
 cg2d: Sum(rhs),rhsMax =   2.10942374678780E-15  1.12112051381376E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
  early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.141635654680537D+01
 --> objf_hflux_tut(bi,bj) =  0.685376854919853D-08
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.185615632956738D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.140673784968281D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.152098128542393D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
  local fc =  0.620023202518704D+01
 global fc =  0.620023202518704D+01
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  6.20023202518704E+00
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    1 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -4.66293670342566E-15  2.54952604282385E+00
 cg2d: Sum(rhs),rhsMax =   3.94129173741931E-15  2.62363057992942E+00
 cg2d: Sum(rhs),rhsMax =   1.96509475358653E-14  2.39895957779109E+00
 cg2d: Sum(rhs),rhsMax =  -1.73194791841524E-14  2.21798267529292E+00
 cg2d: Sum(rhs),rhsMax =   7.10542735760100E-15  2.01662326576495E+00
 cg2d: Sum(rhs),rhsMax =  -1.77635683940025E-15  1.79933840571969E+00
 cg2d: Sum(rhs),rhsMax =   5.10702591327572E-14  1.57603319716532E+00
 cg2d: Sum(rhs),rhsMax =  -4.44089209850063E-15  1.36017289259794E+00
 cg2d: Sum(rhs),rhsMax =   4.37427871702312E-14  1.18689213328779E+00
 cg2d: Sum(rhs),rhsMax =   5.32907051820075E-15  1.12112051378184E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
  early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.141635708756645D+01
 --> objf_hflux_tut(bi,bj) =  0.685376854919853D-08
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.185615632957816D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.140673784968063D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.152098128542264D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
  local fc =  0.620023256595543D+01
 global fc =  0.620023256595543D+01
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  6.20023256595543E+00
grad-res -------------------------------
 grad-res     0    1   43    2    1    1    1    1   6.20023228182E+00  6.20023202519E+00  6.20023256596E+00
 grad-res     0    1    1    1    0    1    1    1  -2.70384203444E-06 -2.70384194145E-06  3.43918964418E-08
(PID.TID 0000.0001)  ADM  ref_cost_function      =  6.20023228182329E+00
(PID.TID 0000.0001)  ADM  adjoint_gradient       = -2.70384203444404E-06
(PID.TID 0000.0001)  ADM  finite-diff_grad       = -2.70384194145379E-06
(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            2        2315           2
 ph-grd _loc: bi, bj, icomptest, ichknum            1           1           1           2
 ph-grd -->hit<--           44           2           1           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   44    2    1 ; 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) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    1 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -4.66293670342566E-15  2.54952604282385E+00
 cg2d: Sum(rhs),rhsMax =   3.94129173741931E-15  2.62363057992942E+00
 cg2d: Sum(rhs),rhsMax =   1.11022302462516E-15  2.39895957779139E+00
 cg2d: Sum(rhs),rhsMax =  -8.88178419700125E-16  2.21798267529377E+00
 cg2d: Sum(rhs),rhsMax =  -3.19744231092045E-14  2.01662326576567E+00
 cg2d: Sum(rhs),rhsMax =   9.76996261670138E-15  1.79933840571689E+00
 cg2d: Sum(rhs),rhsMax =   1.99840144432528E-14  1.57603319715508E+00
 cg2d: Sum(rhs),rhsMax =  -1.28785870856518E-14  1.36017289258229E+00
 cg2d: Sum(rhs),rhsMax =  -2.08721928629529E-14  1.18689213332388E+00
 cg2d: Sum(rhs),rhsMax =  -3.08642000845794E-14  1.12112051381844E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
  early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.141635653982773D+01
 --> objf_hflux_tut(bi,bj) =  0.365120182432684D-08
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.185615632953011D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.140673784968275D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.152098128542403D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
  local fc =  0.620023201176702D+01
 global fc =  0.620023201176702D+01
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  6.20023201176702E+00
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    1 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -4.66293670342566E-15  2.54952604282385E+00
 cg2d: Sum(rhs),rhsMax =   3.94129173741931E-15  2.62363057992942E+00
 cg2d: Sum(rhs),rhsMax =   8.77076189453874E-15  2.39895957779107E+00
 cg2d: Sum(rhs),rhsMax =   1.42108547152020E-14  2.21798267529180E+00
 cg2d: Sum(rhs),rhsMax =  -3.55271367880050E-15  2.01662326576646E+00
 cg2d: Sum(rhs),rhsMax =   2.08721928629529E-14  1.79933840571824E+00
 cg2d: Sum(rhs),rhsMax =  -6.66133814775094E-15  1.57603319716198E+00
 cg2d: Sum(rhs),rhsMax =  -2.79776202205539E-14  1.36017289259333E+00
 cg2d: Sum(rhs),rhsMax =   2.68673971959288E-14  1.18689213328727E+00
 cg2d: Sum(rhs),rhsMax =   1.39888101102770E-14  1.12112051377675E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
  early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.141635709454111D+01
 --> objf_hflux_tut(bi,bj) =  0.365120182432684D-08
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.185615632961544D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.140673784968070D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.152098128542254D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
  local fc =  0.620023256656219D+01
 global fc =  0.620023256656219D+01
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  6.20023256656219E+00
grad-res -------------------------------
 grad-res     0    2   44    2    1    1    1    1   6.20023228182E+00  6.20023201177E+00  6.20023256656E+00
 grad-res     0    2    2    2    0    1    1    1  -2.77397605796E-06 -2.77397586235E-06  7.05170787318E-08
(PID.TID 0000.0001)  ADM  ref_cost_function      =  6.20023228182329E+00
(PID.TID 0000.0001)  ADM  adjoint_gradient       = -2.77397605795953E-06
(PID.TID 0000.0001)  ADM  finite-diff_grad       = -2.77397586234684E-06
(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            3        2315           3
 ph-grd _loc: bi, bj, icomptest, ichknum            1           1           2           3
 ph-grd -->hit<--           45           2           1           1
(PID.TID 0000.0001) grdchk pos: i,j,k=   45    2    1 ; 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) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    1 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -4.66293670342566E-15  2.54952604282385E+00
 cg2d: Sum(rhs),rhsMax =   3.94129173741931E-15  2.62363057992942E+00
 cg2d: Sum(rhs),rhsMax =   3.21964677141295E-15  2.39895957779141E+00
 cg2d: Sum(rhs),rhsMax =  -1.28785870856518E-14  2.21798267529371E+00
 cg2d: Sum(rhs),rhsMax =   2.22044604925031E-15  2.01662326576542E+00
 cg2d: Sum(rhs),rhsMax =  -3.06421554796543E-14  1.79933840571639E+00
 cg2d: Sum(rhs),rhsMax =  -3.33066907387547E-14  1.57603319715464E+00
 cg2d: Sum(rhs),rhsMax =  -1.82076576038526E-14  1.36017289258113E+00
 cg2d: Sum(rhs),rhsMax =   1.66533453693773E-14  1.18689213332621E+00
 cg2d: Sum(rhs),rhsMax =  -1.55431223447522E-15  1.12112051381833E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
  early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.141635655112324D+01
 --> objf_hflux_tut(bi,bj) =  0.478358297175945D-08
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.185615632654069D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.140673784968276D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.152098128542402D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
  local fc =  0.620023202233788D+01
 global fc =  0.620023202233788D+01
(PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus  =  6.20023202233788E+00
(PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001) 
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) 
 EXTERNAL_FIELDS_LOAD, it=         0 : Reading new data, i0,i1=   12    1 (prev=   12    1 )
(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
 cg2d: Sum(rhs),rhsMax =  -4.66293670342566E-15  2.54952604282385E+00
 cg2d: Sum(rhs),rhsMax =   3.88578058618805E-16  2.62363057992942E+00
 cg2d: Sum(rhs),rhsMax =   1.07691633388640E-14  2.39895957779110E+00
 cg2d: Sum(rhs),rhsMax =   1.28785870856518E-14  2.21798267529208E+00
 cg2d: Sum(rhs),rhsMax =  -1.33226762955019E-15  2.01662326576635E+00
 cg2d: Sum(rhs),rhsMax =  -1.50990331349021E-14  1.79933840571899E+00
 cg2d: Sum(rhs),rhsMax =  -1.73194791841524E-14  1.57603319716224E+00
 cg2d: Sum(rhs),rhsMax =  -2.35367281220533E-14  1.36017289259477E+00
 cg2d: Sum(rhs),rhsMax =   1.35447209004269E-14  1.18689213328675E+00
 cg2d: Sum(rhs),rhsMax =   4.21884749357559E-15  1.12112051377829E+00
(PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE
  early fc =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.141635708324548D+01
 --> objf_hflux_tut(bi,bj) =  0.478358297175945D-08
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.185615633260498D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.140673784968068D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
 --> objf_test(bi,bj)   =  0.000000000000000D+00
 --> objf_tracer(bi,bj) =  0.000000000000000D+00
 --> objf_atl(bi,bj)    =  0.000000000000000D+00
 --> objf_temp_tut(bi,bj)   =  0.152098128542254D+01
 --> objf_hflux_tut(bi,bj) =  0.000000000000000D+00
  local fc =  0.620023256052084D+01
 global fc =  0.620023256052084D+01
(PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus =  6.20023256052084E+00
grad-res -------------------------------
 grad-res     0    3   45    2    1    1    1    1   6.20023228182E+00  6.20023202234E+00  6.20023256052E+00
 grad-res     0    3    3    3    0    1    1    1  -2.69091500991E-06 -2.69091481542E-06  7.22781396956E-08
(PID.TID 0000.0001)  ADM  ref_cost_function      =  6.20023228182329E+00
(PID.TID 0000.0001)  ADM  adjoint_gradient       = -2.69091500991183E-06
(PID.TID 0000.0001)  ADM  finite-diff_grad       = -2.69091481541750E-06
(PID.TID 0000.0001) ====== End of gradient-check number   3 (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-01
(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    43     2     1    1    1   0.000000000E+00 -1.000000000E-01
(PID.TID 0000.0001) grdchk output (c):   1  6.2002322818233E+00  6.2002320251870E+00  6.2002325659554E+00
(PID.TID 0000.0001) grdchk output (g):   1    -2.7038419414538E-06 -2.7038420344440E-06  3.4391896441832E-08
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   2    44     2     1    1    1   0.000000000E+00 -1.000000000E-01
(PID.TID 0000.0001) grdchk output (c):   2  6.2002322818233E+00  6.2002320117670E+00  6.2002325665622E+00
(PID.TID 0000.0001) grdchk output (g):   2    -2.7739758623468E-06 -2.7739760579595E-06  7.0517078731847E-08
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk output (p):   3    45     2     1    1    1   0.000000000E+00 -1.000000000E-01
(PID.TID 0000.0001) grdchk output (c):   3  6.2002322818233E+00  6.2002320223379E+00  6.2002325605208E+00
(PID.TID 0000.0001) grdchk output (g):   3    -2.6909148154175E-06 -2.6909150099118E-06  7.2278139695570E-08
(PID.TID 0000.0001) 
(PID.TID 0000.0001) grdchk  summary  :  RMS of    3 ratios =  6.1588934102352E-08
(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:   55.410000000000004
(PID.TID 0000.0001)         System time:   1.2700000000000000
(PID.TID 0000.0001)     Wall clock time:   56.830909013748169
(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:  8.00000000000000017E-002
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  8.45730304718017578E-002
(PID.TID 0000.0001)          No. starts:           1
(PID.TID 0000.0001)           No. stops:           1
(PID.TID 0000.0001)   Seconds in section "THE_MAIN_LOOP (F)      [THE_MODEL_MAIN]":
(PID.TID 0000.0001)           User time:   7.3300000000000001
(PID.TID 0000.0001)         System time:   1.2000000000000000
(PID.TID 0000.0001)     Wall clock time:   8.5549681186676025
(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.28000000000000114
(PID.TID 0000.0001)         System time:  3.00000000000000266E-002
(PID.TID 0000.0001)     Wall clock time:  0.31980967521667480
(PID.TID 0000.0001)          No. starts:           7
(PID.TID 0000.0001)           No. stops:           7
(PID.TID 0000.0001)   Seconds in section "MAIN LOOP           [THE_MAIN_LOOP]":
(PID.TID 0000.0001)           User time:   38.419999999999995
(PID.TID 0000.0001)         System time:   1.1799999999999993
(PID.TID 0000.0001)     Wall clock time:   39.688550233840942
(PID.TID 0000.0001)          No. starts:           7
(PID.TID 0000.0001)           No. stops:           7
(PID.TID 0000.0001)   Seconds in section "MAIN_DO_LOOP        [THE_MAIN_LOOP]":
(PID.TID 0000.0001)           User time:   38.389999999999993
(PID.TID 0000.0001)         System time:   1.1799999999999993
(PID.TID 0000.0001)     Wall clock time:   39.664992332458496
(PID.TID 0000.0001)          No. starts:           7
(PID.TID 0000.0001)           No. stops:           7
(PID.TID 0000.0001)   Seconds in section "LOAD_FIELDS_DRIVER  [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  6.00000000000022737E-002
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  4.45992946624755859E-002
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "EXTERNAL_FLDS_LOAD [LOAD_FLDS_DRIVER]":
(PID.TID 0000.0001)           User time:  4.99999999999971578E-002
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  4.31044101715087891E-002
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "DO_ATMOSPHERIC_PHYS [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   0.0000000000000000
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  7.67707824707031250E-004
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "DO_OCEANIC_PHYS     [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   7.2799999999999798
(PID.TID 0000.0001)         System time:  0.13000000000000012
(PID.TID 0000.0001)     Wall clock time:   7.4238917827606201
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "THERMODYNAMICS      [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   9.5000000000000355
(PID.TID 0000.0001)         System time:  0.20000000000000018
(PID.TID 0000.0001)     Wall clock time:   9.7092530727386475
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "DYNAMICS            [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   20.120000000000033
(PID.TID 0000.0001)         System time:  0.61999999999999988
(PID.TID 0000.0001)     Wall clock time:   20.794983863830566
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "SOLVE_FOR_PRESSURE  [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:   7.0199999999999960
(PID.TID 0000.0001)         System time:  1.00000000000000089E-002
(PID.TID 0000.0001)     Wall clock time:   7.0264461040496826
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "MOM_CORRECTION_STEP [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.22000000000000597
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  0.24728155136108398
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "INTEGR_CONTINUITY   [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.73999999999998778
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  0.72832798957824707
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "TRC_CORRECTION_STEP [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.11999999999999034
(PID.TID 0000.0001)         System time:  1.00000000000000089E-002
(PID.TID 0000.0001)     Wall clock time:  0.13600325584411621
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "BLOCKING_EXCHANGES  [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.36000000000000654
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  0.35502862930297852
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "COST_TILE           [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  0.47000000000000597
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  0.49166226387023926
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "DO_THE_MODEL_IO     [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  2.00000000000031264E-002
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  2.26461887359619141E-002
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "DO_WRITE_PICKUP     [FORWARD_STEP]":
(PID.TID 0000.0001)           User time:  3.00000000000011369E-002
(PID.TID 0000.0001)         System time:  1.00000000000000089E-002
(PID.TID 0000.0001)     Wall clock time:  4.09979820251464844E-002
(PID.TID 0000.0001)          No. starts:          79
(PID.TID 0000.0001)           No. stops:          79
(PID.TID 0000.0001)   Seconds in section "COST_FINAL         [ADJOINT SPIN-DOWN]":
(PID.TID 0000.0001)           User time:  3.00000000000011369E-002
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  2.33662128448486328E-002
(PID.TID 0000.0001)          No. starts:           7
(PID.TID 0000.0001)           No. stops:           7
(PID.TID 0000.0001)   Seconds in section "THE_MAIN_LOOP (A)      [THE_MODEL_MAIN]":
(PID.TID 0000.0001)           User time:   16.609999999999999
(PID.TID 0000.0001)         System time:  5.00000000000000444E-002
(PID.TID 0000.0001)     Wall clock time:   16.703561067581177
(PID.TID 0000.0001)          No. starts:           1
(PID.TID 0000.0001)           No. stops:           1
(PID.TID 0000.0001)   Seconds in section "CTRL_PACK           [THE_MODEL_MAIN]":
(PID.TID 0000.0001)           User time:   0.0000000000000000
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  5.98597526550292969E-003
(PID.TID 0000.0001)          No. starts:           1
(PID.TID 0000.0001)           No. stops:           1
(PID.TID 0000.0001)   Seconds in section "CTRL_PACK     [THE_MODEL_MAIN]":
(PID.TID 0000.0001)           User time:   0.0000000000000000
(PID.TID 0000.0001)         System time:   0.0000000000000000
(PID.TID 0000.0001)     Wall clock time:  5.05590438842773438E-003
(PID.TID 0000.0001)          No. starts:           1
(PID.TID 0000.0001)           No. stops:           1
(PID.TID 0000.0001)   Seconds in section "GRDCHK_MAIN         [THE_MODEL_MAIN]":
(PID.TID 0000.0001)           User time:   31.390000000000004
(PID.TID 0000.0001)         System time:  2.00000000000000178E-002
(PID.TID 0000.0001)     Wall clock time:   31.476690053939819
(PID.TID 0000.0001)          No. starts:           1
(PID.TID 0000.0001)           No. stops:           1
(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 =          96128
(PID.TID 0000.0001) //      Max. barrier spins =              1
(PID.TID 0000.0001) //      Min. barrier spins =              1
(PID.TID 0000.0001) //     Total barrier spins =          96128
(PID.TID 0000.0001) //      Avg. barrier spins =       1.00E+00
PROGRAM MAIN: Execution ended Normally
