/[MITgcm]/MITgcm/verification/matrix_example/results/output.txt

Diff of /MITgcm/verification/matrix_example/results/output.txt

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

-revision 1.4 by jmc,
Thu Oct 13 01:21:32 2005 UTC
+revision 1.6 by jmc,
Mon Nov  7 21:59:42 2011 UTC
 Line 5
  (PID.TID 0000.0001) // ======================================================
  (PID.TID 0000.0001) // execution environment starting up...
  (PID.TID 0000.0001)
- (PID.TID 0000.0001) // MITgcmUV version:  checkpoint57s_post
+ (PID.TID 0000.0001) // MITgcmUV version:  checkpoint63e
  (PID.TID 0000.0001) // Build user:        jmc
- (PID.TID 0000.0001) // Build host:        faulks
+ (PID.TID 0000.0001) // Build host:        baudelaire
- (PID.TID 0000.0001) // Build date:        Wed Sep 28 21:18:19 EDT 2005
+ (PID.TID 0000.0001) // Build date:        Mon Nov  7 16:06:29 EST 2011
  (PID.TID 0000.0001)
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Execution Environment parameter file "eedata"
 Line 18
  (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,nTy=1
+ (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) ># Note: Some systems use & as the namelist terminator (as shown here).
- (PID.TID 0000.0001) ># namelist terminator. Other systems
+ (PID.TID 0000.0001) >#       Other systems use a / character.
- (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" )
 Line 30
  (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 =    1 ; /* No. tiles in X per process */
+ (PID.TID 0000.0001)      nSx =    2 ; /* No. tiles in X per process */
- (PID.TID 0000.0001)      nSy =    1 ; /* No. tiles in Y per process */
+ (PID.TID 0000.0001)      nSy =    4 ; /* No. tiles in Y per process */
- (PID.TID 0000.0001)      sNx =   31 ; /* Tile size in X */
+ (PID.TID 0000.0001)      sNx =   16 ; /* Tile size in X */
- (PID.TID 0000.0001)      sNy =   31 ; /* Tile size in Y */
+ (PID.TID 0000.0001)      sNy =    8 ; /* Tile size in Y */
  (PID.TID 0000.0001)      OLx =    3 ; /* Tile overlap distance in X */
  (PID.TID 0000.0001)      OLy =    3 ; /* Tile overlap distance in Y */
  (PID.TID 0000.0001)      nTx =    1 ; /* No. threads in X per process */
  (PID.TID 0000.0001)      nTy =    1 ; /* No. threads in Y per process */
  (PID.TID 0000.0001)       Nr =    1 ; /* No. levels in the vertical   */
- (PID.TID 0000.0001)       nX =   31 ; /* Total domain size in X ( = nPx*nSx*sNx ) */
+ (PID.TID 0000.0001)       Nx =   32 ; /* Total domain size in X ( = nPx*nSx*sNx ) */
- (PID.TID 0000.0001)       nY =   31 ; /* Total domain size in Y ( = nPy*nSy*sNy ) */
+ (PID.TID 0000.0001)       Ny =   32 ; /* Total domain size in Y ( = nPy*nSy*sNy ) */
- (PID.TID 0000.0001)   nTiles =    1 ; /* Total no. tiles per process ( = nSx*nSy ) */
+ (PID.TID 0000.0001)   nTiles =    8 ; /* 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) 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:   1,   1:   1)
+ (PID.TID 0000.0001) // -o- Thread   1, tiles (   1:   2,   1:   4)
  (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. = 000001)
+ (PID.TID 0000.0001) //        WEST: Tile = 000002, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000001
+ (PID.TID 0000.0001) //        EAST: Tile = 000002, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000001
+ (PID.TID 0000.0001) //       SOUTH: Tile = 000007, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000004
+ (PID.TID 0000.0001) //       NORTH: Tile = 000003, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000002
+ (PID.TID 0000.0001) // Tile number: 000002 (process no. = 000001)
  (PID.TID 0000.0001) //        WEST: Tile = 000001, Process = 000001, Comm = put
  (PID.TID 0000.0001) //                bi = 000001, bj = 000001
  (PID.TID 0000.0001) //        EAST: Tile = 000001, Process = 000001, Comm = put
  (PID.TID 0000.0001) //                bi = 000001, bj = 000001
+ (PID.TID 0000.0001) //       SOUTH: Tile = 000008, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000004
+ (PID.TID 0000.0001) //       NORTH: Tile = 000004, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000002
+ (PID.TID 0000.0001) // Tile number: 000003 (process no. = 000001)
+ (PID.TID 0000.0001) //        WEST: Tile = 000004, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000002
+ (PID.TID 0000.0001) //        EAST: Tile = 000004, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000002
  (PID.TID 0000.0001) //       SOUTH: Tile = 000001, Process = 000001, Comm = put
  (PID.TID 0000.0001) //                bi = 000001, bj = 000001
+ (PID.TID 0000.0001) //       NORTH: Tile = 000005, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000003
+ (PID.TID 0000.0001) // Tile number: 000004 (process no. = 000001)
+ (PID.TID 0000.0001) //        WEST: Tile = 000003, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000002
+ (PID.TID 0000.0001) //        EAST: Tile = 000003, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000002
+ (PID.TID 0000.0001) //       SOUTH: Tile = 000002, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000001
+ (PID.TID 0000.0001) //       NORTH: Tile = 000006, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000003
+ (PID.TID 0000.0001) // Tile number: 000005 (process no. = 000001)
+ (PID.TID 0000.0001) //        WEST: Tile = 000006, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000003
+ (PID.TID 0000.0001) //        EAST: Tile = 000006, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000003
+ (PID.TID 0000.0001) //       SOUTH: Tile = 000003, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000002
+ (PID.TID 0000.0001) //       NORTH: Tile = 000007, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000004
+ (PID.TID 0000.0001) // Tile number: 000006 (process no. = 000001)
+ (PID.TID 0000.0001) //        WEST: Tile = 000005, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000003
+ (PID.TID 0000.0001) //        EAST: Tile = 000005, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000003
+ (PID.TID 0000.0001) //       SOUTH: Tile = 000004, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000002
+ (PID.TID 0000.0001) //       NORTH: Tile = 000008, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000004
+ (PID.TID 0000.0001) // Tile number: 000007 (process no. = 000001)
+ (PID.TID 0000.0001) //        WEST: Tile = 000008, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000004
+ (PID.TID 0000.0001) //        EAST: Tile = 000008, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000004
+ (PID.TID 0000.0001) //       SOUTH: Tile = 000005, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000003
  (PID.TID 0000.0001) //       NORTH: Tile = 000001, Process = 000001, Comm = put
  (PID.TID 0000.0001) //                bi = 000001, bj = 000001
+ (PID.TID 0000.0001) // Tile number: 000008 (process no. = 000001)
+ (PID.TID 0000.0001) //        WEST: Tile = 000007, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000004
+ (PID.TID 0000.0001) //        EAST: Tile = 000007, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000001, bj = 000004
+ (PID.TID 0000.0001) //       SOUTH: Tile = 000006, Process = 000001, Comm = put
+ (PID.TID 0000.0001) //                bi = 000002, bj = 000003
+ (PID.TID 0000.0001) //       NORTH: Tile = 000002, Process = 000001, 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) // Model parameter file "data"
+ (PID.TID 0000.0001) // Parameter file "data"
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) ># ====================
  (PID.TID 0000.0001) ># | Model parameters |
-Line 89
+Line 157
  (PID.TID 0000.0001) > rhonil=1035.,
  (PID.TID 0000.0001) > gravity=9.81,
  (PID.TID 0000.0001) > eosType='LINEAR',
+ (PID.TID 0000.0001) > tAlpha=2.E-4,
+ (PID.TID 0000.0001) > sBeta =0.E-4,
  (PID.TID 0000.0001) ># ivdc_kappa=50.,
  (PID.TID 0000.0001) ># implicitViscosity=.TRUE.,
  (PID.TID 0000.0001) ># implicitDiffusion=.TRUE.,
- (PID.TID 0000.0001) ># useOldFreezing=.TRUE.,
- (PID.TID 0000.0001) > useCDscheme=.FALSE.,
- (PID.TID 0000.0001) ># useNHMTerms=.TRUE.,
  (PID.TID 0000.0001) ># turn on looped cells
  (PID.TID 0000.0001) ># hFacMin=.2,
  (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) > globalFiles=.TRUE.,
- (PID.TID 0000.0001) > useSingleCpuIO=.TRUE.
  (PID.TID 0000.0001) > beta=1.E-11,
  (PID.TID 0000.0001) ># f0=1.e-4,
- (PID.TID 0000.0001) > tAlpha=2.E-4,
- (PID.TID 0000.0001) > sBeta =0.E-4,
- (PID.TID 0000.0001) ># gBaro=9.81,
  (PID.TID 0000.0001) > rigidLid=.FALSE.,
  (PID.TID 0000.0001) > implicitFreeSurface=.TRUE.,
- (PID.TID 0000.0001) ># tempAdvScheme       = 2
+ (PID.TID 0000.0001) ># set precision of data files
- (PID.TID 0000.0001) ># saltAdvScheme       = 2
+ (PID.TID 0000.0001) ># readBinaryPrec=32,
  (PID.TID 0000.0001) > writeBinaryPrec=64,
- (PID.TID 0000.0001) > multiDimAdvection   = .TRUE.
+ (PID.TID 0000.0001) > useSingleCpuIO=.TRUE.,
  (PID.TID 0000.0001) > /
  (PID.TID 0000.0001) >
  (PID.TID 0000.0001) ># Elliptic solver parameters
-Line 126
+Line 186
  (PID.TID 0000.0001) >#startTime = 0.,
  (PID.TID 0000.0001) > nIter0  = 200000,
  (PID.TID 0000.0001) > nTimeSteps = 10,
- (PID.TID 0000.0001) > deltaTmom = 20000.0,
+ (PID.TID 0000.0001) > deltaTmom   = 20000.0,
- (PID.TID 0000.0001) ># tauCD = 172800.,
+ (PID.TID 0000.0001) > deltaTtracer= 20000.0,
- (PID.TID 0000.0001) > deltaTtracer = 20000.0,
  (PID.TID 0000.0001) > deltaTClock = 20000.0,
  (PID.TID 0000.0001) > abEps = 0.1,
  (PID.TID 0000.0001) >#cAdjFreq = -1,
+ (PID.TID 0000.0001) > pickupStrictlyMatch=.FALSE.,
  (PID.TID 0000.0001) > pChkptFreq= 500000000.0,
  (PID.TID 0000.0001) > dumpFreq=   200000000.0,
- (PID.TID 0000.0001) > taveFreq=   0.,
- (PID.TID 0000.0001) ># taveDumpFreq=155520000.,
  (PID.TID 0000.0001) > monitorFreq=2.,
+ (PID.TID 0000.0001) > monitorSelect=2,
  (PID.TID 0000.0001) ># tauThetaClimRelax = 864000.0,
  (PID.TID 0000.0001) ># tauSaltClimRelax = 0.0,
  (PID.TID 0000.0001) ># periodicExternalForcing=.FALSE.,
-Line 145
+Line 204
  (PID.TID 0000.0001) ># Gridding parameters
  (PID.TID 0000.0001) > &PARM04
  (PID.TID 0000.0001) > usingCartesianGrid=.TRUE.,
- (PID.TID 0000.0001) > usingSphericalPolarGrid=.FALSE.,
+ (PID.TID 0000.0001) > delX=32*50.E3,
- (PID.TID 0000.0001) > delX=31*50E3,
+ (PID.TID 0000.0001) > delY=32*50.E3,
- (PID.TID 0000.0001) > delY=31*50E3,
  (PID.TID 0000.0001) > delR=5000.,
- (PID.TID 0000.0001) ># phiMin=0.,
+ (PID.TID 0000.0001) > ygOrigin=-50.E3,
  (PID.TID 0000.0001) > /
  (PID.TID 0000.0001) >
  (PID.TID 0000.0001) ># Input datasets
  (PID.TID 0000.0001) > &PARM05
- (PID.TID 0000.0001) > bathyFile='topog.box',
+ (PID.TID 0000.0001) > bathyFile='topo_box.bin',
  (PID.TID 0000.0001) > hydrogThetaFile=,
  (PID.TID 0000.0001) > hydrogSaltFile=,
- (PID.TID 0000.0001) > zonalWindFile='windx.cos_y',
+ (PID.TID 0000.0001) > zonalWindFile='taux_cosY.bin',
  (PID.TID 0000.0001) > meridWindFile=,
  (PID.TID 0000.0001) > thetaClimFile=,
  (PID.TID 0000.0001) > saltClimFile=,
  (PID.TID 0000.0001) > /
  (PID.TID 0000.0001)
- (PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM01
+ (PID.TID 0000.0001)  INI_PARMS ; starts to read PARM01
- (PID.TID 0000.0001) S/R INI_PARMS ; read PARM01 : OK
+ (PID.TID 0000.0001)  INI_PARMS ; read PARM01 : OK
- (PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM02
+ (PID.TID 0000.0001)  INI_PARMS ; starts to read PARM02
- (PID.TID 0000.0001) S/R INI_PARMS ; read PARM02 : OK
+ (PID.TID 0000.0001)  INI_PARMS ; read PARM02 : OK
- (PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM03
+ (PID.TID 0000.0001)  INI_PARMS ; starts to read PARM03
- (PID.TID 0000.0001) S/R INI_PARMS ; read PARM03 : OK
+ (PID.TID 0000.0001)  INI_PARMS ; read PARM03 : OK
- (PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM04
+ (PID.TID 0000.0001)  INI_PARMS ; starts to read PARM04
- (PID.TID 0000.0001) S/R INI_PARMS ; read PARM04 : OK
+ (PID.TID 0000.0001)  INI_PARMS ; read PARM04 : OK
- (PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM05
+ (PID.TID 0000.0001)  INI_PARMS ; starts to read PARM05
- (PID.TID 0000.0001) S/R INI_PARMS ; read PARM05 : OK
+ (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) // =======================================================
-Line 180
+Line 239
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) ># Packages
  (PID.TID 0000.0001) > &PACKAGES
- (PID.TID 0000.0001) > usePTRACERS=.TRUE.
+ (PID.TID 0000.0001) > usePTRACERS=.TRUE.,
- (PID.TID 0000.0001) > useMATRIX=.TRUE.
+ (PID.TID 0000.0001) > useMATRIX  =.TRUE.,
  (PID.TID 0000.0001) > /
  (PID.TID 0000.0001) >
  (PID.TID 0000.0001)
-Line 201
+Line 260
  (PID.TID 0000.0001) > PTRACERS_advScheme(1)=30,
  (PID.TID 0000.0001) > PTRACERS_diffKh(1)=5.E3,
  (PID.TID 0000.0001) > PTRACERS_diffKr(1)=5.E-5,
- (PID.TID 0000.0001) > PTRACERS_initialFile(1)='trinit1',
+ (PID.TID 0000.0001) > PTRACERS_initialFile(1)='tr1_ini.bin',
  (PID.TID 0000.0001) > /
  (PID.TID 0000.0001) >
  (PID.TID 0000.0001)
-Line 215
+Line 274
  (PID.TID 0000.0001) PTRACERS_Iter0 = /* timestep number when tracers are initialized */
  (PID.TID 0000.0001)                  200000
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) PTRACERS_addSrelax2EmP =/* add Salt relaxation to EmP */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) PTRACERS_dTLev =   /* Ptracer timestep ( s ) */
+ (PID.TID 0000.0001)                 2.000000000000000E+04       /* K =  1 */
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) PTRACERS_dumpFreq = /* Frequency^-1 for snapshot output (s) */
  (PID.TID 0000.0001)                 0.000000000000000E+00
  (PID.TID 0000.0001)     ;
-Line 238
+Line 303
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001)  -----------------------------------
  (PID.TID 0000.0001)  tracer number :    1
+ (PID.TID 0000.0001) PTRACERS_names = /* Tracer short name */
+ (PID.TID 0000.0001)               ''
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) PTRACERS_long_names = /* Tracer long name */
+ (PID.TID 0000.0001)               ''
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) PTRACERS_ioLabel = /* tracer IO Label */
+ (PID.TID 0000.0001)               '01'
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) PTRACERS_advScheme = /* Advection Scheme */
  (PID.TID 0000.0001)                      30
  (PID.TID 0000.0001)     ;
-Line 256
+Line 330
  (PID.TID 0000.0001) PTRACERS_useGMRedi = /* apply GM-Redi */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) PTRACERS_useDWNSLP = /* apply DOWN-SLOPE Flow */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) PTRACERS_useKPP = /* apply KPP scheme */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) PTRACERS_ref = /* Reference vertical profile */
+ (PID.TID 0000.0001)                 0.000000000000000E+00       /* K =  1 */
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) PTRACERS_EvPrRn =/* tracer conc. in Evap. & Rain */
+ (PID.TID 0000.0001)                 1.234567000000000E+05
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001)  -----------------------------------
  (PID.TID 0000.0001)  MATRIX_READPARMS: opening data.matrix
  (PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.matrix
-Line 266
+Line 349
  (PID.TID 0000.0001) // Parameter file "data.matrix"
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) > &MATRIX_PARM01
- (PID.TID 0000.0001) > expMatrixWriteTime=20000.0,
+ (PID.TID 0000.0001) > expMatrixWriteTime=20000.,
- (PID.TID 0000.0001) > impMatrixWriteTime=20000.0
+ (PID.TID 0000.0001) > impMatrixWriteTime=20000.,
  (PID.TID 0000.0001) > /
  (PID.TID 0000.0001)
  (PID.TID 0000.0001)  MATRIX_READPARMS: finished reading data.matrix
- (PID.TID 0000.0001) %MON XC_max                       =   1.5250000000000E+06
+ (PID.TID 0000.0001) SET_PARMS: done
+ (PID.TID 0000.0001) Enter INI_VERTICAL_GRID: setInterFDr=    T ; setCenterDr=    F
+ (PID.TID 0000.0001) %MON XC_max                       =   1.5750000000000E+06
  (PID.TID 0000.0001) %MON XC_min                       =   2.5000000000000E+04
- (PID.TID 0000.0001) %MON XC_mean                      =   7.7500000000000E+05
+ (PID.TID 0000.0001) %MON XC_mean                      =   8.0000000000000E+05
- (PID.TID 0000.0001) %MON XC_sd                        =   4.4721359549996E+05
+ (PID.TID 0000.0001) %MON XC_sd                        =   4.6165463281548E+05
- (PID.TID 0000.0001) %MON XG_max                       =   1.5000000000000E+06
+ (PID.TID 0000.0001) %MON XG_max                       =   1.5500000000000E+06
- (PID.TID 0000.0001) %MON XG_min                       =   5.0000000000000E+04
+ (PID.TID 0000.0001) %MON XG_min                       =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON XG_mean                      =   7.7500000000000E+05
- (PID.TID 0000.0001) %MON XG_sd                        =   4.3277207241996E+05
+ (PID.TID 0000.0001) %MON XG_sd                        =   4.6165463281548E+05
  (PID.TID 0000.0001) %MON DXC_max                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DXC_min                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DXC_mean                     =   5.0000000000000E+04
- (PID.TID 0000.0001) %MON DXC_sd                       =   7.2759576141834E-12
+ (PID.TID 0000.0001) %MON DXC_sd                       =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON DXF_max                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DXF_min                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DXF_mean                     =   5.0000000000000E+04
- (PID.TID 0000.0001) %MON DXF_sd                       =   7.2759576141834E-12
+ (PID.TID 0000.0001) %MON DXF_sd                       =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON DXG_max                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DXG_min                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DXG_mean                     =   5.0000000000000E+04
- (PID.TID 0000.0001) %MON DXG_sd                       =   7.2759576141834E-12
+ (PID.TID 0000.0001) %MON DXG_sd                       =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON DXV_max                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DXV_min                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DXV_mean                     =   5.0000000000000E+04
- (PID.TID 0000.0001) %MON DXV_sd                       =   7.2759576141834E-12
+ (PID.TID 0000.0001) %MON DXV_sd                       =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON YC_max                       =   1.5250000000000E+06
- (PID.TID 0000.0001) %MON YC_min                       =   2.5000000000000E+04
+ (PID.TID 0000.0001) %MON YC_min                       =  -2.5000000000000E+04
- (PID.TID 0000.0001) %MON YC_mean                      =   7.7500000000000E+05
+ (PID.TID 0000.0001) %MON YC_mean                      =   7.5000000000000E+05
- (PID.TID 0000.0001) %MON YC_sd                        =   4.4721359549996E+05
+ (PID.TID 0000.0001) %MON YC_sd                        =   4.6165463281548E+05
  (PID.TID 0000.0001) %MON YG_max                       =   1.5000000000000E+06
- (PID.TID 0000.0001) %MON YG_min                       =   5.0000000000000E+04
+ (PID.TID 0000.0001) %MON YG_min                       =  -5.0000000000000E+04
- (PID.TID 0000.0001) %MON YG_mean                      =   7.7500000000000E+05
+ (PID.TID 0000.0001) %MON YG_mean                      =   7.2500000000000E+05
- (PID.TID 0000.0001) %MON YG_sd                        =   4.3277207241996E+05
+ (PID.TID 0000.0001) %MON YG_sd                        =   4.6165463281548E+05
  (PID.TID 0000.0001) %MON DYC_max                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DYC_min                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DYC_mean                     =   5.0000000000000E+04
- (PID.TID 0000.0001) %MON DYC_sd                       =   7.2759576141834E-12
+ (PID.TID 0000.0001) %MON DYC_sd                       =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON DYF_max                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DYF_min                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DYF_mean                     =   5.0000000000000E+04
- (PID.TID 0000.0001) %MON DYF_sd                       =   7.2759576141834E-12
+ (PID.TID 0000.0001) %MON DYF_sd                       =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON DYG_max                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DYG_min                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DYG_mean                     =   5.0000000000000E+04
- (PID.TID 0000.0001) %MON DYG_sd                       =   7.2759576141834E-12
+ (PID.TID 0000.0001) %MON DYG_sd                       =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON DYU_max                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DYU_min                      =   5.0000000000000E+04
  (PID.TID 0000.0001) %MON DYU_mean                     =   5.0000000000000E+04
- (PID.TID 0000.0001) %MON DYU_sd                       =   7.2759576141834E-12
+ (PID.TID 0000.0001) %MON DYU_sd                       =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON RA_max                       =   2.5000000000000E+09
  (PID.TID 0000.0001) %MON RA_min                       =   2.5000000000000E+09
  (PID.TID 0000.0001) %MON RA_mean                      =   2.5000000000000E+09
-Line 343
+Line 428
  (PID.TID 0000.0001) %MON AngleSN_min                  =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON AngleSN_mean                 =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON AngleSN_sd                   =   0.0000000000000E+00
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: topo_box.bin
  (PID.TID 0000.0001) // =======================================================
- (PID.TID 0000.0001) // Field Model R_low (ini_masks_etc) at iteration          1
+ (PID.TID 0000.0001) // Field Model R_low (ini_masks_etc)
  (PID.TID 0000.0001) // CMIN =         -5.000000000000000E+03
  (PID.TID 0000.0001) // CMAX =         -5.000000000000000E+03
  (PID.TID 0000.0001) // CINT =          0.000000000000000E+00
  (PID.TID 0000.0001) // SYMBOLS (CMIN->CMAX): -abcdefghijklmnopqrstuvwxyz+
  (PID.TID 0000.0001) //                  0.0: .
- (PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):(  -2:  34:   1)
+ (PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):(    -2:    35:     1)
- (PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):(  34:  -2:  -1)
+ (PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):(    35:    -2:    -1)
  (PID.TID 0000.0001) // RANGE K (Lo:Hi:Step):(   1:   1:   1)
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // =======================================================
-Line 359
+Line 445
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001)
  (PID.TID 0000.0001) // =======================================================
- (PID.TID 0000.0001) // Field Model Ro_surf (ini_masks_etc) at iteration          1
+ (PID.TID 0000.0001) // Field Model Ro_surf (ini_masks_etc)
  (PID.TID 0000.0001) // CMIN =          1.000000000000000E+32
  (PID.TID 0000.0001) // CMAX =         -1.000000000000000E+32
  (PID.TID 0000.0001) // CINT =          0.000000000000000E+00
  (PID.TID 0000.0001) // SYMBOLS (CMIN->CMAX): -abcdefghijklmnopqrstuvwxyz+
  (PID.TID 0000.0001) //                  0.0: .
- (PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):(  -2:  34:   1)
+ (PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):(    -2:    35:     1)
- (PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):(  34:  -2:  -1)
+ (PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):(    35:    -2:    -1)
  (PID.TID 0000.0001) // RANGE K (Lo:Hi:Step):(   1:   1:   1)
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // =======================================================
-Line 374
+Line 460
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001)
  (PID.TID 0000.0001) // =======================================================
- (PID.TID 0000.0001) // Field hFacC at iteration          1
+ (PID.TID 0000.0001) // Field hFacC at iteration          0
  (PID.TID 0000.0001) // CMIN =          1.000000000000000E+00
  (PID.TID 0000.0001) // CMAX =          1.000000000000000E+00
  (PID.TID 0000.0001) // CINT =          0.000000000000000E+00
  (PID.TID 0000.0001) // SYMBOLS (CMIN->CMAX): -abcdefghijklmnopqrstuvwxyz+
  (PID.TID 0000.0001) //                  0.0: .
- (PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):(  -2:  34:   1)
+ (PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):(    -2:    35:     1)
- (PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):(  34:  -2:  -1)
+ (PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):(    35:    -2:    -1)
  (PID.TID 0000.0001) // RANGE K (Lo:Hi:Step):(   1:   1:   1)
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // =======================================================
-Line 389
+Line 475
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001)
  (PID.TID 0000.0001) // =======================================================
- (PID.TID 0000.0001) // Field hFacW at iteration          1
+ (PID.TID 0000.0001) // Field hFacW at iteration          0
  (PID.TID 0000.0001) // CMIN =          1.000000000000000E+00
  (PID.TID 0000.0001) // CMAX =          1.000000000000000E+00
  (PID.TID 0000.0001) // CINT =          0.000000000000000E+00
  (PID.TID 0000.0001) // SYMBOLS (CMIN->CMAX): -abcdefghijklmnopqrstuvwxyz+
  (PID.TID 0000.0001) //                  0.0: .
- (PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):(  -2:  34:   1)
+ (PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):(    -2:    35:     1)
- (PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):(  34:  -2:  -1)
+ (PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):(    35:    -2:    -1)
  (PID.TID 0000.0001) // RANGE K (Lo:Hi:Step):(   1:   1:   1)
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // =======================================================
-Line 404
+Line 490
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001)
  (PID.TID 0000.0001) // =======================================================
- (PID.TID 0000.0001) // Field hFacS at iteration          1
+ (PID.TID 0000.0001) // Field hFacS at iteration          0
  (PID.TID 0000.0001) // CMIN =          1.000000000000000E+00
  (PID.TID 0000.0001) // CMAX =          1.000000000000000E+00
  (PID.TID 0000.0001) // CINT =          0.000000000000000E+00
  (PID.TID 0000.0001) // SYMBOLS (CMIN->CMAX): -abcdefghijklmnopqrstuvwxyz+
  (PID.TID 0000.0001) //                  0.0: .
- (PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):(  -2:  34:   1)
+ (PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):(    -2:    35:     1)
- (PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):(  34:  -2:  -1)
+ (PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):(    35:    -2:    -1)
  (PID.TID 0000.0001) // RANGE K (Lo:Hi:Step):(   1:   1:   1)
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // END OF FIELD                                          =
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001)
+ (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 :
-Line 431
+Line 518
  (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) tempAdamsBashforth =   /* use Adams-Bashforth time-stepping for Temp */
+ (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)     ;
-Line 443
+Line 536
  (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) saltAdamsBashforth =   /* use Adams-Bashforth time-stepping for Salt */
+ (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) PTRACERS_INIT_FIXED: PTRACERS_SOM_Advection(   1) =     F
+ (PID.TID 0000.0001) PTRACERS_INIT_FIXED: GAD_OlMinSize=  2  0  1
  (PID.TID 0000.0001) %MON fCori_max                    =   1.1525000000000E-04
- (PID.TID 0000.0001) %MON fCori_min                    =   1.0025000000000E-04
+ (PID.TID 0000.0001) %MON fCori_min                    =   9.9750000000000E-05
- (PID.TID 0000.0001) %MON fCori_mean                   =   1.0775000000000E-04
+ (PID.TID 0000.0001) %MON fCori_mean                   =   1.0750000000000E-04
- (PID.TID 0000.0001) %MON fCori_sd                     =   4.4721359549996E-06
+ (PID.TID 0000.0001) %MON fCori_sd                     =   4.6165463281548E-06
  (PID.TID 0000.0001) %MON fCoriG_max                   =   1.1500000000000E-04
- (PID.TID 0000.0001) %MON fCoriG_min                   =   1.0000000000000E-04
+ (PID.TID 0000.0001) %MON fCoriG_min                   =   9.9500000000000E-05
- (PID.TID 0000.0001) %MON fCoriG_mean                  =   1.0750000000000E-04
+ (PID.TID 0000.0001) %MON fCoriG_mean                  =   1.0725000000000E-04
- (PID.TID 0000.0001) %MON fCoriG_sd                    =   4.4721359549996E-06
+ (PID.TID 0000.0001) %MON fCoriG_sd                    =   4.6165463281548E-06
  (PID.TID 0000.0001) %MON fCoriCos_max                 =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON fCoriCos_min                 =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON fCoriCos_mean                =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON fCoriCos_sd                  =   0.0000000000000E+00
  (PID.TID 0000.0001) INI_CG2D: CG2D normalisation factor =  2.0000000000000001E-04
  (PID.TID 0000.0001)
- (PID.TID 0000.0001) CONFIG_CHECK: OK
  (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 = OCEANIC
+ (PID.TID 0000.0001) buoyancyRelation = /* Type of relation to get Buoyancy */
- (PID.TID 0000.0001) fluidIsAir  =  /* fluid major constituent is Air */
+ (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=  /* fuild major constituent is Water */
+ (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 */
-Line 481
+Line 583
  (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 oK ) */
+ (PID.TID 0000.0001) tRef =   /* Reference temperature profile ( oC or K ) */
  (PID.TID 0000.0001)                 2.000000000000000E+01       /* K =  1 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) sRef =   /* Reference salinity profile ( ppt ) */
+ (PID.TID 0000.0001) sRef =   /* Reference salinity profile ( psu ) */
  (PID.TID 0000.0001)                 1.000000000000000E+01       /* K =  1 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) viscAh =   /* Lateral eddy viscosity ( m^2/s ) */
+ (PID.TID 0000.0001) viscAh  =   /* Lateral eddy viscosity ( m^2/s ) */
  (PID.TID 0000.0001)                 5.000000000000000E+03
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) viscAhMax =   /* Maximum lateral eddy viscosity ( m^2/s ) */
+ (PID.TID 0000.0001) viscAhMax = /* Maximum lateral eddy viscosity ( m^2/s ) */
  (PID.TID 0000.0001)                 1.000000000000000E+21
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) viscAhGrid =   /* Grid dependent lateral eddy viscosity ( non-dim. ) */
+ (PID.TID 0000.0001) viscAhGrid = /* Grid dependent lateral eddy viscosity ( non-dim. ) */
  (PID.TID 0000.0001)                 0.000000000000000E+00
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) useFullLeith =   /* Use Full Form of Leith Viscosity on/off flag*/
+ (PID.TID 0000.0001) useFullLeith = /* Use Full Form of Leith Viscosity on/off flag*/
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) useStrainTensionVisc =   /* Use StrainTension Form of Viscous Operator on/off flag*/
+ (PID.TID 0000.0001) useStrainTensionVisc= /* Use StrainTension Form of Viscous Operator flag*/
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) useAreaViscLength =   /* Use area for visc length instead of geom. mean*/
+ (PID.TID 0000.0001) useAreaViscLength = /* Use area for visc length instead of geom. mean*/
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) viscC2leith = /* Leith harmonic visc. factor (on grad(vort),non-dim.) */
  (PID.TID 0000.0001)                 0.000000000000000E+00
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) viscC2leithD = /* Leith harmonic viscosity factor (on grad(div),non-dim.) */
+ (PID.TID 0000.0001) viscC2leithD = /* Leith harmonic viscosity factor (on grad(div),non-dim.)*/
  (PID.TID 0000.0001)                 0.000000000000000E+00
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) viscC2smag =   /* Smagorinsky harmonic viscosity factor (non-dim.) */
+ (PID.TID 0000.0001) viscC2smag = /* Smagorinsky harmonic viscosity factor (non-dim.) */
  (PID.TID 0000.0001)                 0.000000000000000E+00
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) viscA4 =   /* Lateral biharmonic viscosity ( m^4/s ) */
+ (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) viscA4Max =   /* Maximum biharmonic viscosity ( m^2/s ) */
+ (PID.TID 0000.0001) viscA4Max = /* Maximum biharmonic viscosity ( m^2/s ) */
  (PID.TID 0000.0001)                 1.000000000000000E+21
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) viscA4Grid =   /* Grid dependent biharmonic viscosity ( non-dim. ) */
+ (PID.TID 0000.0001) viscA4Grid = /* Grid dependent biharmonic viscosity ( non-dim. ) */
  (PID.TID 0000.0001)                 0.000000000000000E+00
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) viscC4leith = /* Leith biharm viscosity factor (on grad(vort), non-dim.) */
+ (PID.TID 0000.0001) viscC4leith = /* Leith biharm viscosity factor (on grad(vort), non-dim.)*/
  (PID.TID 0000.0001)                 0.000000000000000E+00
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) viscC4leithD = /* Leith biharm viscosity factor (on grad(div), non-dim.) */
-Line 538
+Line 640
  (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) viscAr =   /* Vertical eddy viscosity ( units of r^2/s ) */
+ (PID.TID 0000.0001) viscArNr = /* vertical profile of vertical viscosity ( m^2/s )*/
- (PID.TID 0000.0001)                 1.000000000000000E-02
+ (PID.TID 0000.0001)                 1.000000000000000E-02       /* K =  1 */
  (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 ( 1/s ) */
+ (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 coeff. ( 1/m ) */
+ (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)                 5.000000000000000E+03
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) diffK4T =   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */
+ (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 =   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */
+ (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 )*/
-Line 568
+Line 670
  (PID.TID 0000.0001) diffKrNrS = /* vertical profile of vertical diffusion of Salt ( m^2/s )*/
  (PID.TID 0000.0001)                 0.000000000000000E+00       /* K =  1 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) diffKrBL79surf =   /* Surface diffusion for Bryan and Lewis 1979 ( m^2/s ) */
+ (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) 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) 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) diffKrBL79Ho = /* Turning depth for Bryan and Lewis 1979 ( m ) */
  (PID.TID 0000.0001)                -2.000000000000000E+03
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001)  Equation of State : eosType = LINEAR
+ (PID.TID 0000.0001) ivdc_kappa = /* Implicit Vertical Diffusivity for Convection ( m^2/s) */
- (PID.TID 0000.0001) tAlpha =   /* Linear EOS thermal expansion coefficient ( 1/degree ) */
+ (PID.TID 0000.0001)                 0.000000000000000E+00
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) hMixCriteria=  /* Criteria for mixed-layer diagnostic */
+ (PID.TID 0000.0001)                -8.000000000000000E-01
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) dRhoSmall =  /* Parameter for mixed-layer diagnostic */
+ (PID.TID 0000.0001)                 1.000000000000000E-06
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) hMixSmooth=  /* Smoothing parameter for mixed-layer diagnostic */
+ (PID.TID 0000.0001)                 0.000000000000000E+00
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) eosType =  /* Type of Equation of State */
+ (PID.TID 0000.0001)               'LINEAR'
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) tAlpha =   /* Linear EOS thermal expansion coefficient ( 1/oC ) */
  (PID.TID 0000.0001)                 2.000000000000000E-04
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) sBeta =   /* Linear EOS haline contraction coefficient ( 1/ppt ) */
+ (PID.TID 0000.0001) sBeta =   /* Linear EOS haline contraction coefficient ( 1/psu ) */
  (PID.TID 0000.0001)                 0.000000000000000E+00
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) rhonil =   /* Reference density ( kg/m^3 ) */
-Line 593
+Line 709
  (PID.TID 0000.0001) rhoConst =   /* Reference density ( 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)                 1.000000000000000E+00       /* K =  1 */
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) rhoFacF = /* normalized Reference density @ W-Interface (-) */
+ (PID.TID 0000.0001)     2 @  1.000000000000000E+00              /* K =  1:  2 */
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) rhoConstFresh =   /* Reference density ( kg/m^3 ) */
  (PID.TID 0000.0001)                 1.035000000000000E+03
  (PID.TID 0000.0001)     ;
-Line 603
+Line 725
  (PID.TID 0000.0001)                 9.810000000000000E+00
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) rotationPeriod =   /* Rotation Period ( s ) */
- (PID.TID 0000.0001)                 8.640000000000000E+04
+ (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.272205216643040E-05
+ (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
-Line 614
+Line 736
  (PID.TID 0000.0001) beta =   /* Beta ( 1/(m.s) ) */
  (PID.TID 0000.0001)                 9.999999999999999E-12
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) freeSurfFac =   /* Implicit free surface factor */
+ (PID.TID 0000.0001) fPrime =   /* Second coriolis parameter ( 1/s ) */
- (PID.TID 0000.0001)                 1.000000000000000E+00
+ (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) rigidLid =   /* Rigid lid on/off flag */
+ (PID.TID 0000.0001) freeSurfFac =   /* Implicit free surface factor */
- (PID.TID 0000.0001)                   F
+ (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) 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) implicDiv2Dflow =  /* Barot. Flow Div. implicit factor (0-1)*/
  (PID.TID 0000.0001)                 1.000000000000000E+00
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) exactConserv =   /* Exact Volume Conservation on/off flag*/
+ (PID.TID 0000.0001) exactConserv =  /* Exact Volume Conservation on/off flag*/
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) linFSConserveTr = /* Tracer correction for Lin Free Surface on/off flag*/
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) uniformLin_PhiSurf =   /* use uniform Bo_surf on/off flag*/
+ (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) nonlinFreeSurf =   /* Non-linear Free Surf. options (-1,0,1,2,3)*/
+ (PID.TID 0000.0001) hFacMin =   /* minimum partial cell factor (hFac) */
- (PID.TID 0000.0001)                       0
+ (PID.TID 0000.0001)                 1.000000000000000E+00
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) hFacMinDr = /* minimum partial cell thickness ( m) */
+ (PID.TID 0000.0001)                 1.000000000000000E+00
+ (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* Coordinate options (not yet implemented)*/
+ (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) selectAddFluid = /* option for mass source/sink of fluid (=0: off) */
+ (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)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) convertFW2Salt = /* convert F.W. Flux to Salt Flux (-1=use local S)(ppt)*/
+ (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) 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)                 3.500000000000000E+01
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) use3Dsolver = /* use 3-D pressure solver on/off flag */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) nonHydrostatic =  /* Non-Hydrostatic on/off flag */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) nh_Am2 = /* Non-Hydrostatic terms scaling factor */
  (PID.TID 0000.0001)                 1.000000000000000E+00
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) implicitNHPress = /* Non-Hyd Pressure implicit factor (0-1)*/
+ (PID.TID 0000.0001)                 1.000000000000000E+00
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) selectNHfreeSurf = /* Non-Hyd (free-)Surface option */
+ (PID.TID 0000.0001)                       0
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) quasiHydrostatic = /* Quasi-Hydrostatic on/off flag */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) 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) 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)                   F
+ (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)                       1
+ (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)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) useJamartWetPoints= /* Coriolis WetPoints method flag */
+ (PID.TID 0000.0001) useEnergyConservingCoriolis= /* Flx-Form Coriolis scheme 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) useJamartWetPoints= /* Coriolis WetPoints method flag */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) SadournyCoriolis= /* Sadourny Coriolis discr. flag */
+ (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) upwindVorticity= /* Upwind bias vorticity flag */
+ (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) useAbsVorticity= /* Work with f+zeta in Coriolis */
+ (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= /* High order interp. of vort. flag */
+ (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= /* Upwind vertical Shear advection flag */
+ (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= /* Kinetic Energy scheme selector */
+ (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 */
-Line 711
+Line 896
  (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)     ;
-Line 720
+Line 908
  (PID.TID 0000.0001) useMultiDimAdvec =   /* Multi-Dim Advection is/is-not used */
  (PID.TID 0000.0001)                   T
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) implicitDiffusion =/* Implicit Diffusion on/off flag */
+ (PID.TID 0000.0001) implicitDiffusion = /* Implicit Diffusion on/off flag */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) tempStepping =  /* Temperature equation on/off flag */
-Line 729
+Line 917
  (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) 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) 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) 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)  debugMode  =  /* Debug Mode on/off flag */
+ (PID.TID 0000.0001) saltIsActiveTr = /* Salt  is a dynamically Active Tracer */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001)    debLevA  =  /* 1rst level of debugging */
+ (PID.TID 0000.0001)  readBinaryPrec = /* Precision used for reading binary files */
- (PID.TID 0000.0001)                       1
+ (PID.TID 0000.0001)                      32
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001)    debLevB  =  /* 2nd  level of debugging */
+ (PID.TID 0000.0001) writeBinaryPrec = /* Precision used for writing binary files */
- (PID.TID 0000.0001)                       2
+ (PID.TID 0000.0001)                      64
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001)  debugLevel =  /* select debugging level */
+ (PID.TID 0000.0001)  globalFiles = /* write "global" (=not per tile) files */
- (PID.TID 0000.0001)                       1
+ (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)                   T
+ (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)                       2
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) //
  (PID.TID 0000.0001) // Elliptic solver(s) paramters ( PARM02 in namelist )
-Line 777
+Line 981
  (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) nIter0 =   /* Run starting timestep number  */
+ (PID.TID 0000.0001) deltaTmom =   /* Momentum equation timestep ( s ) */
- (PID.TID 0000.0001)                  200000
- (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) deltatTmom =   /* Momentum equation timestep ( s ) */
  (PID.TID 0000.0001)                 2.000000000000000E+04
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) deltaTfreesurf =   /* FreeSurface equation timestep ( s ) */
+ (PID.TID 0000.0001) deltaTfreesurf = /* FreeSurface equation timestep ( s ) */
  (PID.TID 0000.0001)                 2.000000000000000E+04
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dTtracerLev =   /* Tracer equation timestep ( s ) */
+ (PID.TID 0000.0001) dTtracerLev =  /* Tracer equation timestep ( s ) */
  (PID.TID 0000.0001)                 2.000000000000000E+04       /* K =  1 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) deltatTClock  =   /* Model clock timestep ( s ) */
+ (PID.TID 0000.0001) deltaTClock  =   /* Model clock timestep ( s ) */
  (PID.TID 0000.0001)                 2.000000000000000E+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) forcing_In_AB =   /* put T,S Forcing in Adams-Bash. stepping */
+ (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) baseTime =   /* Model base time ( s ). */
+ (PID.TID 0000.0001) pickupStrictlyMatch= /* stop if pickup do not strictly match */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) nIter0   =   /* Run starting timestep number */
+ (PID.TID 0000.0001)                  200000
+ (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)                  200010
+ (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) startTime =  /* Run start time ( s ) */
  (PID.TID 0000.0001)                 4.000000000000000E+09
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) endTime =   /* Integration ending time ( s ). */
+ (PID.TID 0000.0001) endTime  =   /* Integration ending time ( s ) */
  (PID.TID 0000.0001)                 4.000200000000000E+09
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) pChkPtFreq =   /* Permanent restart/checkpoint file interval ( s ). */
+ (PID.TID 0000.0001) pChkPtFreq = /* Permanent restart/pickup file interval ( s ) */
  (PID.TID 0000.0001)                 5.000000000000000E+08
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) chkPtFreq =   /* Rolling restart/checkpoint file interval ( s ). */
+ (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. */
-Line 831
+Line 1056
  (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.000000000000000E+08
  (PID.TID 0000.0001)     ;
-Line 843
+Line 1071
  (PID.TID 0000.0001) monitorFreq =   /* Monitor output interval ( s ). */
  (PID.TID 0000.0001)                 2.000000000000000E+00
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) monitorSelect = /* select group of variables to monitor */
+ (PID.TID 0000.0001)                       2
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) monitor_stdio =   /* Model IO flag. */
  (PID.TID 0000.0001)                   T
  (PID.TID 0000.0001)     ;
-Line 859
+Line 1090
  (PID.TID 0000.0001)                 0.000000000000000E+00
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) latBandClimRelax =   /* max. Lat. where relaxation */
- (PID.TID 0000.0001)                 4.650000000000000E+06
+ (PID.TID 0000.0001)                 4.800000000000000E+06
  (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) usingCartesianGrid = /* Cartesian coordinates flag ( True/False ) */
  (PID.TID 0000.0001)                   T
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) usingSphericalPolarGrid =   /* Spherical coordinates flag ( True / False ) */
+ (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)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) usingCylindricalGrid =   /* Spherical coordinates flag ( True / False ) */
+ (PID.TID 0000.0001) usingCurvilinearGrid = /* Curvilinear coordinates flag ( True/False ) */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) Ro_SeaLevel =   /* r(1) ( units of r ) */
+ (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) horiVertRatio =   /* Ratio on units : Horiz - Vertical */
+ (PID.TID 0000.0001) gravitySign = /* gravity orientation relative to vertical coordinate */
- (PID.TID 0000.0001)                 1.000000000000000E+00
+ (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+03       /* K =  1 */
-Line 888
+Line 1134
  (PID.TID 0000.0001) drF =   /* W spacing ( units of r ) */
  (PID.TID 0000.0001)                 5.000000000000000E+03       /* K =  1 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) delX =   /* U spacing ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) delX = /* U spacing ( m - cartesian, degrees - spherical ) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* I =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) delY =   /* V spacing ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) delY = /* V spacing ( m - cartesian, degrees - spherical ) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* J =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) phiMin =    /* South edge (ignored - cartesian, degrees - spherical ) */
+ (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) thetaMin =    /* West edge ( ignored - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) ygOrigin = /* Y-axis origin of South edge (cartesian: m, lat-lon: deg.) */
- (PID.TID 0000.0001)                 0.000000000000000E+00
+ (PID.TID 0000.0001)                -5.000000000000000E+04
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) rSphere =    /* Radius ( ignored - cartesian, m - spherical ) */
+ (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) xcoord =   /* P-point X coord ( m - cartesian, degrees - spherical ) */
+ (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.500000000000000E+04,      /* I =  1 */
  (PID.TID 0000.0001)                 7.500000000000000E+04,      /* I =  2 */
  (PID.TID 0000.0001)                 1.250000000000000E+05,      /* I =  3 */
-Line 934
+Line 1183
  (PID.TID 0000.0001)                 1.375000000000000E+06,      /* I = 28 */
  (PID.TID 0000.0001)                 1.425000000000000E+06,      /* I = 29 */
  (PID.TID 0000.0001)                 1.475000000000000E+06,      /* I = 30 */
- (PID.TID 0000.0001)                 1.525000000000000E+06       /* I = 31 */
+ (PID.TID 0000.0001)                 1.525000000000000E+06,      /* I = 31 */
+ (PID.TID 0000.0001)                 1.575000000000000E+06       /* I = 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) ycoord =   /* P-point Y coord ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) yC =  /* yC(1,:,1,:) : P-point Y coord ( deg. or m if cartesian) */
- (PID.TID 0000.0001)                 2.500000000000000E+04,      /* J =  1 */
+ (PID.TID 0000.0001)                -2.500000000000000E+04,      /* J =  1 */
- (PID.TID 0000.0001)                 7.500000000000000E+04,      /* J =  2 */
+ (PID.TID 0000.0001)                 2.500000000000000E+04,      /* J =  2 */
- (PID.TID 0000.0001)                 1.250000000000000E+05,      /* J =  3 */
+ (PID.TID 0000.0001)                 7.500000000000000E+04,      /* J =  3 */
- (PID.TID 0000.0001)                 1.750000000000000E+05,      /* J =  4 */
+ (PID.TID 0000.0001)                 1.250000000000000E+05,      /* J =  4 */
- (PID.TID 0000.0001)                 2.250000000000000E+05,      /* J =  5 */
+ (PID.TID 0000.0001)                 1.750000000000000E+05,      /* J =  5 */
- (PID.TID 0000.0001)                 2.750000000000000E+05,      /* J =  6 */
+ (PID.TID 0000.0001)                 2.250000000000000E+05,      /* J =  6 */
- (PID.TID 0000.0001)                 3.250000000000000E+05,      /* J =  7 */
+ (PID.TID 0000.0001)                 2.750000000000000E+05,      /* J =  7 */
- (PID.TID 0000.0001)                 3.750000000000000E+05,      /* J =  8 */
+ (PID.TID 0000.0001)                 3.250000000000000E+05,      /* J =  8 */
- (PID.TID 0000.0001)                 4.250000000000000E+05,      /* J =  9 */
+ (PID.TID 0000.0001)                 3.750000000000000E+05,      /* J =  9 */
- (PID.TID 0000.0001)                 4.750000000000000E+05,      /* J = 10 */
+ (PID.TID 0000.0001)                 4.250000000000000E+05,      /* J = 10 */
- (PID.TID 0000.0001)                 5.250000000000000E+05,      /* J = 11 */
+ (PID.TID 0000.0001)                 4.750000000000000E+05,      /* J = 11 */
- (PID.TID 0000.0001)                 5.750000000000000E+05,      /* J = 12 */
+ (PID.TID 0000.0001)                 5.250000000000000E+05,      /* J = 12 */
- (PID.TID 0000.0001)                 6.250000000000000E+05,      /* J = 13 */
+ (PID.TID 0000.0001)                 5.750000000000000E+05,      /* J = 13 */
- (PID.TID 0000.0001)                 6.750000000000000E+05,      /* J = 14 */
+ (PID.TID 0000.0001)                 6.250000000000000E+05,      /* J = 14 */
- (PID.TID 0000.0001)                 7.250000000000000E+05,      /* J = 15 */
+ (PID.TID 0000.0001)                 6.750000000000000E+05,      /* J = 15 */
- (PID.TID 0000.0001)                 7.750000000000000E+05,      /* J = 16 */
+ (PID.TID 0000.0001)                 7.250000000000000E+05,      /* J = 16 */
- (PID.TID 0000.0001)                 8.250000000000000E+05,      /* J = 17 */
+ (PID.TID 0000.0001)                 7.750000000000000E+05,      /* J = 17 */
- (PID.TID 0000.0001)                 8.750000000000000E+05,      /* J = 18 */
+ (PID.TID 0000.0001)                 8.250000000000000E+05,      /* J = 18 */
- (PID.TID 0000.0001)                 9.250000000000000E+05,      /* J = 19 */
+ (PID.TID 0000.0001)                 8.750000000000000E+05,      /* J = 19 */
- (PID.TID 0000.0001)                 9.750000000000000E+05,      /* J = 20 */
+ (PID.TID 0000.0001)                 9.250000000000000E+05,      /* J = 20 */
- (PID.TID 0000.0001)                 1.025000000000000E+06,      /* J = 21 */
+ (PID.TID 0000.0001)                 9.750000000000000E+05,      /* J = 21 */
- (PID.TID 0000.0001)                 1.075000000000000E+06,      /* J = 22 */
+ (PID.TID 0000.0001)                 1.025000000000000E+06,      /* J = 22 */
- (PID.TID 0000.0001)                 1.125000000000000E+06,      /* J = 23 */
+ (PID.TID 0000.0001)                 1.075000000000000E+06,      /* J = 23 */
- (PID.TID 0000.0001)                 1.175000000000000E+06,      /* J = 24 */
+ (PID.TID 0000.0001)                 1.125000000000000E+06,      /* J = 24 */
- (PID.TID 0000.0001)                 1.225000000000000E+06,      /* J = 25 */
+ (PID.TID 0000.0001)                 1.175000000000000E+06,      /* J = 25 */
- (PID.TID 0000.0001)                 1.275000000000000E+06,      /* J = 26 */
+ (PID.TID 0000.0001)                 1.225000000000000E+06,      /* J = 26 */
- (PID.TID 0000.0001)                 1.325000000000000E+06,      /* J = 27 */
+ (PID.TID 0000.0001)                 1.275000000000000E+06,      /* J = 27 */
- (PID.TID 0000.0001)                 1.375000000000000E+06,      /* J = 28 */
+ (PID.TID 0000.0001)                 1.325000000000000E+06,      /* J = 28 */
- (PID.TID 0000.0001)                 1.425000000000000E+06,      /* J = 29 */
+ (PID.TID 0000.0001)                 1.375000000000000E+06,      /* J = 29 */
- (PID.TID 0000.0001)                 1.475000000000000E+06,      /* J = 30 */
+ (PID.TID 0000.0001)                 1.425000000000000E+06,      /* J = 30 */
- (PID.TID 0000.0001)                 1.525000000000000E+06       /* J = 31 */
+ (PID.TID 0000.0001)                 1.475000000000000E+06,      /* J = 31 */
+ (PID.TID 0000.0001)                 1.525000000000000E+06       /* J = 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) rcoord =   /* P-point R coordinate (  units of r ) */
+ (PID.TID 0000.0001) rcoord = /* P-point R coordinate (  units of r ) */
  (PID.TID 0000.0001)                -2.500000000000000E+03       /* K =  1 */
  (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+03       /* K =  2 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dxF =  /* dxF(:,1,:,1) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) deepFacC = /* deep-model grid factor @ cell-Center (-) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* I =  1: 31 */
+ (PID.TID 0000.0001)                 1.000000000000000E+00       /* K =  1 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dxF =  /* dxF(1,:,1,:) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) deepFacF = /* deep-model grid factor @ W-Interface (-) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* J =  1: 31 */
+ (PID.TID 0000.0001)     2 @  1.000000000000000E+00              /* K =  1:  2 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dyF =  /* dyF(:,1,:,1) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) rVel2wUnit = /* convert units: rVel -> wSpeed (=1 if z-coord)*/
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* I =  1: 31 */
+ (PID.TID 0000.0001)     2 @  1.000000000000000E+00              /* K =  1:  2 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dyF =  /* dyF(1,:,1,:) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) wUnit2rVel = /* convert units: wSpeed -> rVel (=1 if z-coord)*/
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* J =  1: 31 */
+ (PID.TID 0000.0001)     2 @  1.000000000000000E+00              /* K =  1:  2 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dxG =  /* dxG(:,1,:,1) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dBdrRef = /* Vertical grad. of reference buoyancy [(m/s/r)^2] */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* I =  1: 31 */
+ (PID.TID 0000.0001)                 0.000000000000000E+00       /* K =  1 */
+ (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) dxG =  /* dxG(1,:,1,:) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) phiEuler = /* Euler angle, rotation about original z-coordinate [rad] */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* J =  1: 31 */
+ (PID.TID 0000.0001)                 0.000000000000000E+00
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dyG =  /* dyG(:,1,:,1) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) thetaEuler = /* Euler angle, rotation about new x-coordinate [rad] */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* I =  1: 31 */
+ (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) dyG =  /* dyG(1,:,1,:) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dxF =  /* dxF(:,1,:,1) ( units: m ) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* J =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dxC =  /* dxC(:,1,:,1) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dxF =  /* dxF(1,:,1,:) ( units: m ) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* I =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dxC =  /* dxC(1,:,1,:) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dyF =  /* dyF(:,1,:,1) ( units: m ) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* J =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dyC =  /* dyC(:,1,:,1) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dyF =  /* dyF(1,:,1,:) ( units: m ) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* I =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dyC =  /* dyC(1,:,1,:) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dxG =  /* dxG(:,1,:,1) ( units: m ) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* J =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dxV =  /* dxV(:,1,:,1) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dxG =  /* dxG(1,:,1,:) ( units: m ) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* I =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dxV =  /* dxV(1,:,1,:) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dyG =  /* dyG(:,1,:,1) ( units: m ) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* J =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dyU =  /* dyU(:,1,:,1) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dyG =  /* dyG(1,:,1,:) ( units: m ) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* I =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dyU =  /* dyU(1,:,1,:) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dxC =  /* dxC(:,1,:,1) ( units: m ) */
- (PID.TID 0000.0001)    31 @  5.000000000000000E+04              /* J =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) rA =  /* rA(:,1,:,1) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dxC =  /* dxC(1,:,1,:) ( units: m ) */
- (PID.TID 0000.0001)    31 @  2.500000000000000E+09              /* I =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) rA =  /* rA(1,:,1,:) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dyC =  /* dyC(:,1,:,1) ( units: m ) */
- (PID.TID 0000.0001)    31 @  2.500000000000000E+09              /* J =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) rAw =  /* rAw(:,1,:,1) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dyC =  /* dyC(1,:,1,:) ( units: m ) */
- (PID.TID 0000.0001)    31 @  2.500000000000000E+09              /* I =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) rAw =  /* rAw(1,:,1,:) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dxV =  /* dxV(:,1,:,1) ( units: m ) */
- (PID.TID 0000.0001)    31 @  2.500000000000000E+09              /* J =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) rAs =  /* rAs(:,1,:,1) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dxV =  /* dxV(1,:,1,:) ( units: m ) */
- (PID.TID 0000.0001)    31 @  2.500000000000000E+09              /* I =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) rAs =  /* rAs(1,:,1,:) ( m - cartesian, degrees - spherical ) */
+ (PID.TID 0000.0001) dyU =  /* dyU(:,1,:,1) ( units: m ) */
- (PID.TID 0000.0001)    31 @  2.500000000000000E+09              /* J =  1: 31 */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) dyU =  /* dyU(1,:,1,:) ( units: m ) */
+ (PID.TID 0000.0001)    32 @  5.000000000000000E+04              /* J =  1: 32 */
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) rA  =  /* rA (:,1,:,1) ( units: m^2 ) */
+ (PID.TID 0000.0001)    32 @  2.500000000000000E+09              /* I =  1: 32 */
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) rA  =  /* rA (1,:,1,:) ( units: m^2 ) */
+ (PID.TID 0000.0001)    32 @  2.500000000000000E+09              /* J =  1: 32 */
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) rAw =  /* rAw(:,1,:,1) ( units: m^2 ) */
+ (PID.TID 0000.0001)    32 @  2.500000000000000E+09              /* I =  1: 32 */
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) rAw =  /* rAw(1,:,1,:) ( units: m^2 ) */
+ (PID.TID 0000.0001)    32 @  2.500000000000000E+09              /* J =  1: 32 */
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) rAs =  /* rAs(:,1,:,1) ( units: m^2 ) */
+ (PID.TID 0000.0001)    32 @  2.500000000000000E+09              /* I =  1: 32 */
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) rAs =  /* rAs(1,:,1,:) ( units: m^2 ) */
+ (PID.TID 0000.0001)    32 @  2.500000000000000E+09              /* J =  1: 32 */
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) globalArea = /* Integrated horizontal Area (m^2) */
+ (PID.TID 0000.0001)                 2.212500000000000E+12
+ (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) GAD_CHECK: #define ALLOW_GENERIC_ADVDIFF
+ (PID.TID 0000.0001) // =======================================================
+ (PID.TID 0000.0001) // CONFIG_CHECK : Normal End
+ (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001)
+ (PID.TID 0000.0001)  MDS_READ_META: opening file: pickup.0000200000.meta
+ (PID.TID 0000.0001)  nRecords =  11 ; filePrec =  64 ; fileIter =    200000
+ (PID.TID 0000.0001)     nDims =   2 , dims:
+ (PID.TID 0000.0001)    1:  32   1  32
+ (PID.TID 0000.0001)    2:  32   1  32
+ (PID.TID 0000.0001)     nFlds =   0 , nFl3D =   0 , fields:
+ (PID.TID 0000.0001)   nTimRec =   0 , timeList:
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001) READ_MFLDS_CHECK: - normal end ; reset MFLDS file-name: pickup.0000200000
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: taux_cosY.bin
+ (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: tr1_ini.bin
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Model current state
  (PID.TID 0000.0001) // =======================================================
-Line 1054
+Line 1364
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0000000000000E+09
  (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785897865594E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594533267691E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472183666459E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472183666483E-06
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574659699777E-02
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   1.5101959851848E-03
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.4948430419136E-04
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868966374318E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735246373737E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   3.8479431497424E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   3.8479432529240E-10
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548220375E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408284497952E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326100839942E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769772804388E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599180849815E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -3.4790032297163E-10
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -3.4790032264254E-10
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608824106565E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259699447259E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373805962517E-04
  (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.5007691569626E-10
  (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.4463438512757E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   8.2514292094692E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   8.2514292197460E-21
  (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.4348255593201E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.6888856562121E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   5.6892816450992E-12
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   2.0000000000000E+01
-Line 1093
+Line 1403
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759337772E-06
  (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603605988726E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473007937E-06
+ (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473007916E-06
  (PID.TID 0000.0001) %MON vort_p_mean                  =   1.1765497076023E-04
  (PID.TID 0000.0001) %MON vort_p_sd                    =   5.6390875812523E-05
  (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.6502858439492E-19
-Line 1101
+Line 1411
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-  cg2d: Sum(rhs),rhsMax =   1.29178857595467E-05  7.14831396830014E-01
+ (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
- (PID.TID 0000.0001)                    cg2d_init_res =    1.18737609102005E-06
+  cg2d: Sum(rhs),rhsMax =   1.29178857548595E-05  7.14831396830014E-01
+ (PID.TID 0000.0001)                    cg2d_init_res =    1.18737609142484E-06
  (PID.TID 0000.0001)                       cg2d_iters =    12
- (PID.TID 0000.0001)                         cg2d_res =    9.80030121455239E-08
+ (PID.TID 0000.0001)                         cg2d_res =    9.80030121560800E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 1112
+Line 1423
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0000200000000E+09
  (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785897698663E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594504756348E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3474823363410E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3474823364294E-06
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574664295944E-02
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   1.5101960936734E-03
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.4948430251763E-04
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868962871165E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735249995004E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -3.0969828630414E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -3.0969824086257E-10
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548369538E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408277703333E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326092987266E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769772979212E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599181707757E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -2.1260621651301E-09
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -2.1260621657529E-09
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608824990142E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259699288049E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373806964704E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.4764107977971E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.4764106171206E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -4.1070490945131E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -4.1070489119738E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   5.9235338542589E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   5.6690808732883E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.4198834485406E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.4198834566473E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0842429779504E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   6.7575867429804E-12
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   2.0000000000000E+01
-Line 1142
+Line 1453
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347585148466E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079091916849E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   3.2856392756104E-09
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   3.2856391295791E-09
  (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522693753644E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004695949551E-02
-Line 1151
+Line 1462
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833758984423E-06
  (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603605703614E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130472993007E-06
+ (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130472993037E-06
  (PID.TID 0000.0001) %MON vort_p_mean                  =   1.1765497076023E-04
  (PID.TID 0000.0001) %MON vort_p_sd                    =   5.6390875807088E-05
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.1847067658338E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.1338161746577E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   5.9235338291692E-20
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   5.6690808732883E-20
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-  Writing explicit matrix : 1 1 1
+ (PID.TID 0000.0001) Writing explicit matrix :   1     1     1
-  Writing implicit matrix : 1 1 1
+ (PID.TID 0000.0001) Writing implicit matrix :   1     1     1
-  cg2d: Sum(rhs),rhsMax =   1.29182940463695E-05  7.14831409113866E-01
+  cg2d: Sum(rhs),rhsMax =   1.29182940431516E-05  7.14831409113866E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.17016480835591E-06
+ (PID.TID 0000.0001)                    cg2d_init_res =    1.17016481009178E-06
  (PID.TID 0000.0001)                       cg2d_iters =    17
- (PID.TID 0000.0001)                         cg2d_res =    8.01765680838308E-08
+ (PID.TID 0000.0001)                         cg2d_res =    8.01765682060520E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 1172
+Line 1483
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0000400000000E+09
  (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785901123444E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594502499767E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472844558302E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472844558446E-06
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574675727548E-02
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   1.5101960166437E-03
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.4948430210992E-04
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868963285239E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735244543210E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -3.2511475273962E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -3.2511464199561E-10
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548754140E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408278888938E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326092898174E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773366640E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599180900710E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599180900709E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.3522315324072E-10
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.3522318375365E-10
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608824943704E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259699517013E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373807238716E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.4668356589973E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.4668355248868E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.2037677466869E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.2037683585659E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -8.1737367909844E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -2.6053683253500E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   7.7071974524504E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   7.7071974717058E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.4897446889964E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   5.7006267237595E-12
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   2.0000000000000E+01
-Line 1201
+Line 1512
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347585314096E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093466560E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093466561E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.5630141973495E-09
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.5630146868527E-09
  (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522710928395E-06
+ (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522710928396E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004696744475E-02
- (PID.TID 0000.0001) %MON ke_mean                      =   9.4023106287814E-04
+ (PID.TID 0000.0001) %MON ke_mean                      =   9.4023106287815E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.1062500000000E+16
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759183815E-06
  (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603606867501E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130472998718E-06
+ (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130472998737E-06
  (PID.TID 0000.0001) %MON vort_p_mean                  =   1.1765497076023E-04
  (PID.TID 0000.0001) %MON vort_p_sd                    =   5.6390875806388E-05
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -1.6347473602041E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -5.2107366506999E-20
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -8.1737368010203E-20
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -2.6053683253500E-20
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-  Writing explicit matrix : 1 2 1
+ (PID.TID 0000.0001) Writing explicit matrix :   1     2     1
-  Writing implicit matrix : 1 2 1
+ (PID.TID 0000.0001) Writing implicit matrix :   1     2     1
-  cg2d: Sum(rhs),rhsMax =   1.29179873280191E-05  7.14831436211391E-01
+  cg2d: Sum(rhs),rhsMax =   1.29179873209484E-05  7.14831436211394E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.16678716151272E-06
+ (PID.TID 0000.0001)                    cg2d_init_res =    1.16678716605319E-06
  (PID.TID 0000.0001)                       cg2d_iters =    17
- (PID.TID 0000.0001)                         cg2d_res =    9.60155228134612E-08
+ (PID.TID 0000.0001)                         cg2d_res =    9.60155243853393E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                200003
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0000600000000E+09
  (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785899266340E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594517347060E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594517347063E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3475494278567E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3475494280297E-06
- (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574666345089E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574666345090E-02
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   1.5101960764780E-03
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.4948431008927E-04
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868961866603E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735243795034E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735243795035E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   5.1237290836000E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   5.1237305701754E-10
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548350721E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408271646423E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326089233615E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773251636E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599177397032E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -1.2601409283126E-10
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -1.2601408525219E-10
- (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608824039992E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608824039993E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259699567303E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373807697362E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   3.8420636728406E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   3.8420640602708E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.6656310223043E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.6656308397651E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -6.8994165308731E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -4.8013098832578E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.2298091440428E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.2298093646288E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.1228954971952E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   6.5454873792451E-12
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   2.0000000000000E+01
-Line 1262
+Line 1573
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347584746641E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093006543E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   3.0736509382725E-09
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   3.0736512482166E-09
  (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522696832280E-06
+ (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522696832281E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004696439068E-02
- (PID.TID 0000.0001) %MON ke_mean                      =   9.4023102105307E-04
+ (PID.TID 0000.0001) %MON ke_mean                      =   9.4023102105309E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.1062500000000E+16
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759247736E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603605628026E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603605628027E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473009715E-06
+ (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473009712E-06
  (PID.TID 0000.0001) %MON vort_p_mean                  =   1.1765497076023E-04
  (PID.TID 0000.0001) %MON vort_p_sd                    =   5.6390875806486E-05
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -1.3798833061746E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -9.6026197665155E-20
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -6.8994165308731E-20
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -4.8013098832578E-20
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-  Writing explicit matrix : 1 3 1
+ (PID.TID 0000.0001) Writing explicit matrix :   1     3     1
-  Writing implicit matrix : 1 3 1
+ (PID.TID 0000.0001) Writing implicit matrix :   1     3     1
-  cg2d: Sum(rhs),rhsMax =   1.29183979177688E-05  7.14831406463998E-01
+  cg2d: Sum(rhs),rhsMax =   1.29183979185044E-05  7.14831406464009E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.00714382982494E-06
+ (PID.TID 0000.0001)                    cg2d_init_res =    1.00714381959586E-06
  (PID.TID 0000.0001)                       cg2d_iters =    13
- (PID.TID 0000.0001)                         cg2d_res =    9.27800094038682E-08
+ (PID.TID 0000.0001)                         cg2d_res =    9.27800132466659E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                200004
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0000800000000E+09
  (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785898191342E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594517219513E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594517219516E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472375530469E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472375532881E-06
- (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574660860784E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574660860785E-02
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   1.5101960048114E-03
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.4948430256415E-04
- (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868960141737E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868960141738E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735249377185E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735249377186E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -2.4376806381101E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -2.4376801736111E-10
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065549692281E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408276963260E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326090200297E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773266080E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599188258600E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599188258601E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -1.9581025397157E-09
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -1.9581025783991E-09
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608825808309E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259699681986E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373807480635E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.9360993839800E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.9361001253128E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -4.0912341997027E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -4.0912351571023E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -9.6157740395260E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -9.9820383471874E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.9178887359750E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.9178891520071E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.9047003233345E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   6.0338578939217E-12
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   2.0000000000000E+01
-Line 1321
+Line 1632
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347584056695E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093064319E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093064320E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   3.2729873597622E-09
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   3.2729881256819E-09
  (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522688592638E-06
+ (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522688592639E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004696598940E-02
- (PID.TID 0000.0001) %MON ke_mean                      =   9.4023111519785E-04
+ (PID.TID 0000.0001) %MON ke_mean                      =   9.4023111519787E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.1062500000000E+16
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833758916884E-06
  (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603606702113E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473004104E-06
+ (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473004132E-06
  (PID.TID 0000.0001) %MON vort_p_mean                  =   1.1765497076023E-04
  (PID.TID 0000.0001) %MON vort_p_sd                    =   5.6390875807409E-05
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -1.9231548080337E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -1.9964076694375E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -9.6157740401683E-20
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -9.9820383471874E-20
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-  Writing explicit matrix : 1 4 1
+ (PID.TID 0000.0001) Writing explicit matrix :   1     4     1
-  Writing implicit matrix : 1 4 1
+ (PID.TID 0000.0001) Writing implicit matrix :   1     4     1
-  cg2d: Sum(rhs),rhsMax =   1.29179148264120E-05  7.14831431051186E-01
+  cg2d: Sum(rhs),rhsMax =   1.29179148262715E-05  7.14831431051196E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.17107066849388E-06
+ (PID.TID 0000.0001)                    cg2d_init_res =    1.17107068503552E-06
  (PID.TID 0000.0001)                       cg2d_iters =    19
- (PID.TID 0000.0001)                         cg2d_res =    6.97646058147324E-08
+ (PID.TID 0000.0001)                         cg2d_res =    6.97646068028989E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                200005
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0001000000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785902732704E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785902732705E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594507560529E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594507560532E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3475311731201E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3475311736921E-06
- (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574681716266E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574681716268E-02
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   1.5101960881379E-03
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.4948431591928E-04
- (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868961106711E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868961106712E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735243747490E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735243747491E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -3.7239079663835E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -3.7239074579678E-10
- (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548180546E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548180547E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408278058108E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326093545990E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773484969E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599188084810E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599188084811E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.1283009414405E-10
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.1283018092105E-10
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608825000703E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259699502017E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373807615593E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.2239063754678E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.2239063447341E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.9653514884412E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.9653501883149E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -2.3032311023881E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   1.7386250129485E-20
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   6.7069704365848E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   6.7069706534392E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   2.9393667121900E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.7734032491395E-12
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   2.0000000000000E+01
-Line 1381
+Line 1692
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347584442685E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093939874E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093939876E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   1.7791251003742E-09
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   1.7791250757873E-09
  (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522719925857E-06
+ (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522719925860E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004696805575E-02
- (PID.TID 0000.0001) %MON ke_mean                      =   9.4023105216735E-04
+ (PID.TID 0000.0001) %MON ke_mean                      =   9.4023105216738E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.1062500000000E+16
- (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759598705E-06
+ (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759598706E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603606322757E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603606322758E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130472999792E-06
+ (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130472999764E-06
  (PID.TID 0000.0001) %MON vort_p_mean                  =   1.1765497076023E-04
  (PID.TID 0000.0001) %MON vort_p_sd                    =   5.6390875804351E-05
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -4.6064622047762E-20
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   3.4772500258970E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -2.3032311023881E-20
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.7386250129485E-19
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-  Writing explicit matrix : 1 5 1
+ (PID.TID 0000.0001) Writing explicit matrix :   1     5     1
-  Writing implicit matrix : 1 5 1
+ (PID.TID 0000.0001) Writing implicit matrix :   1     5     1
-  cg2d: Sum(rhs),rhsMax =   1.29183692857450E-05  7.14831428086815E-01
+  cg2d: Sum(rhs),rhsMax =   1.29183692794133E-05  7.14831428086836E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.28565905300292E-06
+ (PID.TID 0000.0001)                    cg2d_init_res =    1.28565909253697E-06
  (PID.TID 0000.0001)                       cg2d_iters =    19
- (PID.TID 0000.0001)                         cg2d_res =    8.56127615452914E-08
+ (PID.TID 0000.0001)                         cg2d_res =    8.56127643088673E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                200006
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0001200000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785900288035E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785900288036E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594511076651E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594511076654E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3470830987182E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3470830991681E-06
- (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574661234563E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574661234564E-02
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   1.5101959754559E-03
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.4948430958426E-04
- (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868966158534E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868966158535E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735239645686E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735239645687E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   3.2506616721830E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   3.2506624141082E-10
- (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548296737E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548296738E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408277509124E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326093112029E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773824248E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599185452787E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599185452786E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   9.8431416728594E-11
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   9.8431432383041E-11
- (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608825382468E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608825382469E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259700071435E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373808478637E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.9895875044167E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.9895877633244E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.8704262506217E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.8704262655228E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -4.0330181288333E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   9.1993590648481E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.2304088274051E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.2304090930841E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.6659601838956E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   6.1310667548702E-12
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   2.0000000000000E+01
-Line 1440
+Line 1751
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   1.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347586463413E-02
+ (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347586463414E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079095296993E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079095296994E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.3916700035334E-09
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.3916702106595E-09
  (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522689154175E-06
+ (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522689154177E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004697420082E-02
- (PID.TID 0000.0001) %MON ke_mean                      =   9.4023106862756E-04
+ (PID.TID 0000.0001) %MON ke_mean                      =   9.4023106862760E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.1062500000000E+16
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833760157439E-06
  (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603606778246E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473019331E-06
+ (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473019335E-06
  (PID.TID 0000.0001) %MON vort_p_mean                  =   1.1765497076023E-04
  (PID.TID 0000.0001) %MON vort_p_sd                    =   5.6390875806540E-05
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -8.0660361754524E-20
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.8398718129696E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -4.0330180877262E-20
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   9.1993590648481E-20
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-  Writing explicit matrix : 1 6 1
+ (PID.TID 0000.0001) Writing explicit matrix :   1     6     1
-  Writing implicit matrix : 1 6 1
+ (PID.TID 0000.0001) Writing implicit matrix :   1     6     1
-  cg2d: Sum(rhs),rhsMax =   1.29176754107927E-05  7.14831452448160E-01
+  cg2d: Sum(rhs),rhsMax =   1.29176754160021E-05  7.14831452448173E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.00898117217147E-06
+ (PID.TID 0000.0001)                    cg2d_init_res =    1.00898122896268E-06
  (PID.TID 0000.0001)                       cg2d_iters =    12
- (PID.TID 0000.0001)                         cg2d_res =    8.56122025051976E-08
+ (PID.TID 0000.0001)                         cg2d_res =    8.56122077076605E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                200007
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0001400000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785900384354E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785900384355E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594517423044E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594517423048E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3473803442196E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3473803447387E-06
- (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574666126631E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574666126633E-02
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   1.5101960444399E-03
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.4948431850206E-04
- (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868959295819E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868959295821E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735246434857E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735246434859E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   1.6257411041474E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   1.6257412554825E-10
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548016988E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408278720510E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326095407372E-05
- (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773077504E-01
+ (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773077505E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599187702811E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599187702810E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -2.1358485087293E-09
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -2.1358486685855E-09
- (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608825053298E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608825053299E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259699481704E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373807310502E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.9678630344570E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.9678631760180E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.2316893190145E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.2316893395036E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   5.0347982405943E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -6.9450454176062E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.2289373012666E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.2289378984835E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.3867028630584E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   5.5305305484145E-12
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   2.0000000000000E+01
-Line 1501
+Line 1812
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347583718328E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079092310018E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079092310019E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.5853514552116E-09
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.5853514716029E-09
  (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522696504042E-06
+ (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522696504045E-06
- (PID.TID 0000.0001) %MON ke_max                       =   1.4004696211240E-02
+ (PID.TID 0000.0001) %MON ke_max                       =   1.4004696211241E-02
- (PID.TID 0000.0001) %MON ke_mean                      =   9.4023105043368E-04
+ (PID.TID 0000.0001) %MON ke_mean                      =   9.4023105043372E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.1062500000000E+16
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759310686E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603605780201E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603605780202E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473005755E-06
+ (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473005734E-06
  (PID.TID 0000.0001) %MON vort_p_mean                  =   1.1765497076023E-04
  (PID.TID 0000.0001) %MON vort_p_sd                    =   5.6390875803626E-05
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.0069596481189E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -1.3890090835212E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   5.0347982405943E-20
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -6.9450454176062E-20
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-  Writing explicit matrix : 1 7 1
+ (PID.TID 0000.0001) Writing explicit matrix :   1     7     1
-  Writing implicit matrix : 1 7 1
+ (PID.TID 0000.0001) Writing implicit matrix :   1     7     1
-  cg2d: Sum(rhs),rhsMax =   1.29181366816476E-05  7.14831382145328E-01
+  cg2d: Sum(rhs),rhsMax =   1.29181366901321E-05  7.14831382145344E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.06514667403524E-06
+ (PID.TID 0000.0001)                    cg2d_init_res =    1.06514671138306E-06
  (PID.TID 0000.0001)                       cg2d_iters =    17
- (PID.TID 0000.0001)                         cg2d_res =    7.65957967067599E-08
+ (PID.TID 0000.0001)                         cg2d_res =    7.65957981700821E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                200008
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0001600000000E+09
  (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785900968545E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594503936920E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594503936921E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472411332483E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472411339237E-06
- (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574672065167E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574672065171E-02
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   1.5101960291410E-03
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.4948430485202E-04
- (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868963218160E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868963218163E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735245595628E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735245595629E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -4.2115287475889E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -4.2115275950292E-10
- (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065549284518E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065549284519E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408284752083E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326096277653E-05
- (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773800540E-01
+ (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773800541E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599182940113E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599182940112E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -9.2460814655326E-11
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -9.2460769736842E-11
- (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608825868659E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608825868660E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259699978430E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373808056058E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.6578594874591E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.6578603666276E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.4862330537289E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.4862335808575E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   6.1076936853976E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   7.6360558554278E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   7.3628382175180E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   7.3628384678627E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.1854211844086E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.9668075865357E-12
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   2.0000000000000E+01
-Line 1560
+Line 1871
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   1.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347585287264E-02
+ (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347585287265E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079095202162E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079095202163E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.1262875899673E-09
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.1262882933021E-09
  (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522705426046E-06
+ (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522705426051E-06
- (PID.TID 0000.0001) %MON ke_max                       =   1.4004697375028E-02
+ (PID.TID 0000.0001) %MON ke_max                       =   1.4004697375029E-02
- (PID.TID 0000.0001) %MON ke_mean                      =   9.4023110830419E-04
+ (PID.TID 0000.0001) %MON ke_mean                      =   9.4023110830424E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.1062500000000E+16
- (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759546307E-06
+ (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759546308E-06
  (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603606346513E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473004550E-06
+ (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473004535E-06
  (PID.TID 0000.0001) %MON vort_p_mean                  =   1.1765497076023E-04
  (PID.TID 0000.0001) %MON vort_p_sd                    =   5.6390875810885E-05
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.2215387290508E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.5272111710856E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   6.1076936452539E-20
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   7.6360558554278E-20
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-  Writing explicit matrix : 1 8 1
+ (PID.TID 0000.0001) Writing explicit matrix :   1     8     1
-  Writing implicit matrix : 1 8 1
+ (PID.TID 0000.0001) Writing implicit matrix :   1     8     1
-  cg2d: Sum(rhs),rhsMax =   1.29179200165051E-05  7.14831450919323E-01
+  cg2d: Sum(rhs),rhsMax =   1.29179200107910E-05  7.14831450919334E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.11079783471147E-06
+ (PID.TID 0000.0001)                    cg2d_init_res =    1.11079790209677E-06
  (PID.TID 0000.0001)                       cg2d_iters =    18
- (PID.TID 0000.0001)                         cg2d_res =    9.08082423224432E-08
+ (PID.TID 0000.0001)                         cg2d_res =    9.08082512138358E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                200009
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0001800000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785900450428E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785900450429E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594527421711E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594527421720E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3476472939536E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3476472946893E-06
- (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574670300971E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574670300974E-02
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   1.5101961249859E-03
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.4948431760482E-04
- (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868964909377E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868964909380E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735244831866E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735244831868E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   3.1056203679474E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   3.1056227065058E-10
- (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548419667E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548419668E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408284028507E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326096081573E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773392761E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599182338729E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599182338730E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.5253266357636E-10
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.5253261266447E-10
- (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608824999184E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608824999185E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259699624082E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373808125645E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.6066178307228E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.6066179811896E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.9599348236807E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.9599354909733E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   1.1936681192639E-20
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   9.8743377316554E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.7292347605710E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.7292357087263E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.9402317998055E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   6.3133947213713E-12
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   2.0000000000000E+01
-Line 1620
+Line 1931
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   1.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347585963751E-02
+ (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347585963752E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093571043E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093571044E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   3.1679478589445E-09
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   3.1679483927786E-09
  (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522702775592E-06
+ (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522702775597E-06
- (PID.TID 0000.0001) %MON ke_max                       =   1.4004696677850E-02
+ (PID.TID 0000.0001) %MON ke_max                       =   1.4004696677851E-02
- (PID.TID 0000.0001) %MON ke_mean                      =   9.4023105744067E-04
+ (PID.TID 0000.0001) %MON ke_mean                      =   9.4023105744072E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.1062500000000E+16
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759374588E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603605531573E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603605531574E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473014550E-06
+ (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473014525E-06
  (PID.TID 0000.0001) %MON vort_p_mean                  =   1.1765497076023E-04
  (PID.TID 0000.0001) %MON vort_p_sd                    =   5.6390875805417E-05
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   2.3873362385519E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.9748675463311E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.1936681192759E-19
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   9.8743377316554E-20
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-  Writing explicit matrix : 1 9 1
+ (PID.TID 0000.0001) Writing explicit matrix :   1     9     1
-  Writing implicit matrix : 1 9 1
+ (PID.TID 0000.0001) Writing implicit matrix :   1     9     1
-  cg2d: Sum(rhs),rhsMax =   1.29185494079313E-05  7.14831404157621E-01
+  cg2d: Sum(rhs),rhsMax =   1.29185494190942E-05  7.14831404157643E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    9.57340449163656E-07
+ (PID.TID 0000.0001)                    cg2d_init_res =    9.57340490366612E-07
  (PID.TID 0000.0001)                       cg2d_iters =    10
- (PID.TID 0000.0001)                         cg2d_res =    9.79891654007067E-08
+ (PID.TID 0000.0001)                         cg2d_res =    9.79891780984708E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                200010
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0002000000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785899418360E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785899418361E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594507089039E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594507089044E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472039444729E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472039451370E-06
- (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574661884553E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574661884555E-02
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   1.5101959940328E-03
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.4948431042020E-04
- (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868961191916E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868961191919E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735249175660E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735249175662E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   2.3555354848511E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   2.3555353362573E-10
- (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065549206509E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065549206510E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408284343718E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326094458559E-05
- (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773535243E-01
+ (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773535244E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599188270227E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599188270228E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -2.7314094776633E-09
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -2.7314097624269E-09
- (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608826082410E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608826082412E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259700240884E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373808400809E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   3.1423899643123E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   3.1423902697861E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.8826323803514E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.8826323840767E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -6.4332620610712E-22
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.5300064542368E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.4185781304393E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.4185792767575E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.2311129461481E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   6.9505576960878E-12
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =   2.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   2.0000000000000E+01
-Line 1680
+Line 1991
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   1.0000000000000E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347584476766E-02
+ (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347584476768E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079094140974E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079094140976E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.5139119714499E-09
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.5139122158289E-09
  (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522690130736E-06
+ (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522690130740E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004696994316E-02
- (PID.TID 0000.0001) %MON ke_mean                      =   9.4023111433226E-04
+ (PID.TID 0000.0001) %MON ke_mean                      =   9.4023111433234E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.1062500000000E+16
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759963308E-06
  (PID.TID 0000.0001) %MON vort_r_max                   =   1.9603606247359E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473016947E-06
+ (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473016939E-06
  (PID.TID 0000.0001) %MON vort_p_mean                  =   1.1765497076023E-04
  (PID.TID 0000.0001) %MON vort_p_sd                    =   5.6390875808950E-05
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -1.2866524719480E-20
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -3.0600129084735E-20
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -6.4332623597400E-21
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -1.5300064542368E-20
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-  Writing explicit matrix : 1 10 1
+ (PID.TID 0000.0001) Writing explicit matrix :   1    10     1
-  Writing implicit matrix : 1 10 1
+ (PID.TID 0000.0001) Writing implicit matrix :   1    10     1
  (PID.TID 0000.0001) %CHECKPOINT    200010 ckptA
  (PID.TID 0000.0001)   Seconds in section "ALL                    [THE_MODEL_MAIN]":
- (PID.TID 0000.0001)           User time:  0.360944994
+ (PID.TID 0000.0001)           User time:  0.14999999999999999
- (PID.TID 0000.0001)         System time:  0.0359949993
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.589018822
+ (PID.TID 0000.0001)     Wall clock time:  0.17843103408813477
- (PID.TID 0000.0001)          No. starts: 1
+ (PID.TID 0000.0001)          No. starts:           1
- (PID.TID 0000.0001)           No. stops: 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:  0.0439930002
+ (PID.TID 0000.0001)           User time:  2.00000000000000004E-002
- (PID.TID 0000.0001)         System time:  0.0089989996
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.127004147
+ (PID.TID 0000.0001)     Wall clock time:  3.89850139617919922E-002
- (PID.TID 0000.0001)          No. starts: 1
+ (PID.TID 0000.0001)          No. starts:           1
- (PID.TID 0000.0001)           No. stops: 1
+ (PID.TID 0000.0001)           No. stops:           1
  (PID.TID 0000.0001)   Seconds in section "THE_MAIN_LOOP          [THE_MODEL_MAIN]":
- (PID.TID 0000.0001)           User time:  0.307953987
+ (PID.TID 0000.0001)           User time:  0.13000000000000000
- (PID.TID 0000.0001)         System time:  0.0209959988
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.444175959
+ (PID.TID 0000.0001)     Wall clock time:  0.13941287994384766
- (PID.TID 0000.0001)          No. starts: 1
+ (PID.TID 0000.0001)          No. starts:           1
- (PID.TID 0000.0001)           No. stops: 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.00899900123
+ (PID.TID 0000.0001)           User time:  9.99999999999999847E-003
- (PID.TID 0000.0001)         System time:  0.
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.0116119385
+ (PID.TID 0000.0001)     Wall clock time:  1.92849636077880859E-002
- (PID.TID 0000.0001)          No. starts: 1
+ (PID.TID 0000.0001)          No. starts:           1
- (PID.TID 0000.0001)           No. stops: 1
+ (PID.TID 0000.0001)           No. stops:           1
- (PID.TID 0000.0001)   Seconds in section "MONITOR             [THE_MAIN_LOOP]":
- (PID.TID 0000.0001)           User time:  0.00399899855
- (PID.TID 0000.0001)         System time:  0.
- (PID.TID 0000.0001)     Wall clock time:  0.00441098213
- (PID.TID 0000.0001)          No. starts: 1
- (PID.TID 0000.0001)           No. stops: 1
- (PID.TID 0000.0001)   Seconds in section "DO_THE_MODEL_IO     [THE_MAIN_LOOP]":
- (PID.TID 0000.0001)           User time:  0.00300000235
- (PID.TID 0000.0001)         System time:  0.00499900058
- (PID.TID 0000.0001)     Wall clock time:  0.0296480656
- (PID.TID 0000.0001)          No. starts: 1
- (PID.TID 0000.0001)           No. stops: 1
  (PID.TID 0000.0001)   Seconds in section "MAIN LOOP           [THE_MAIN_LOOP]":
- (PID.TID 0000.0001)           User time:  0.291955985
+ (PID.TID 0000.0001)           User time:  0.12000000000000000
- (PID.TID 0000.0001)         System time:  0.0159969982
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.398317814
+ (PID.TID 0000.0001)     Wall clock time:  0.12009620666503906
- (PID.TID 0000.0001)          No. starts: 1
+ (PID.TID 0000.0001)          No. starts:           1
- (PID.TID 0000.0001)           No. stops: 1
+ (PID.TID 0000.0001)           No. stops:           1
  (PID.TID 0000.0001)   Seconds in section "FORWARD_STEP        [THE_MAIN_LOOP]":
- (PID.TID 0000.0001)           User time:  0.291955985
+ (PID.TID 0000.0001)           User time:  0.11999999999999997
- (PID.TID 0000.0001)         System time:  0.0159969982
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.39794302
+ (PID.TID 0000.0001)     Wall clock time:  0.11999726295471191
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
- (PID.TID 0000.0001)   Seconds in section "EXTERNAL_FIELDS_LOAD[FORWARD_STEP]":
+ (PID.TID 0000.0001)   Seconds in section "LOAD_FIELDS_DRIVER  [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.
+ (PID.TID 0000.0001)           User time:   0.0000000000000000
- (PID.TID 0000.0001)         System time:  0.
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.000340938568
+ (PID.TID 0000.0001)     Wall clock time:  2.80857086181640625E-004
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
+ (PID.TID 0000.0001)   Seconds in section "EXTERNAL_FLDS_LOAD [LOAD_FLDS_DRIVER]":
+ (PID.TID 0000.0001)           User time:   0.0000000000000000
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
+ (PID.TID 0000.0001)     Wall clock time:  9.84668731689453125E-005
+ (PID.TID 0000.0001)          No. starts:          10
+ (PID.TID 0000.0001)           No. stops:          10
  (PID.TID 0000.0001)   Seconds in section "DO_ATMOSPHERIC_PHYS [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.
+ (PID.TID 0000.0001)           User time:   0.0000000000000000
- (PID.TID 0000.0001)         System time:  0.
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.000336885452
+ (PID.TID 0000.0001)     Wall clock time:  9.03606414794921875E-005
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
  (PID.TID 0000.0001)   Seconds in section "DO_OCEANIC_PHYS     [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.00799699873
+ (PID.TID 0000.0001)           User time:   0.0000000000000000
- (PID.TID 0000.0001)         System time:  0.
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.0054526329
+ (PID.TID 0000.0001)     Wall clock time:  2.96044349670410156E-003
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
  (PID.TID 0000.0001)   Seconds in section "THERMODYNAMICS      [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.0529929921
+ (PID.TID 0000.0001)           User time:  3.00000000000000266E-002
- (PID.TID 0000.0001)         System time:  0.00100000016
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.0566320419
+ (PID.TID 0000.0001)     Wall clock time:  3.18973064422607422E-002
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
  (PID.TID 0000.0001)   Seconds in section "DYNAMICS            [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.055990994
+ (PID.TID 0000.0001)           User time:  2.00000000000000178E-002
- (PID.TID 0000.0001)         System time:  0.
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.0537729263
+ (PID.TID 0000.0001)     Wall clock time:  2.98519134521484375E-002
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
  (PID.TID 0000.0001)   Seconds in section "SOLVE_FOR_PRESSURE  [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.105984047
+ (PID.TID 0000.0001)           User time:  4.99999999999999889E-002
- (PID.TID 0000.0001)         System time:  0.00299999863
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.108469009
+ (PID.TID 0000.0001)     Wall clock time:  2.74558067321777344E-002
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
- (PID.TID 0000.0001)   Seconds in section "UV_CORRECTION_STEP  [FORWARD_STEP]":
+ (PID.TID 0000.0001)   Seconds in section "MOM_CORRECTION_STEP [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.00499999523
+ (PID.TID 0000.0001)           User time:   0.0000000000000000
- (PID.TID 0000.0001)         System time:  0.
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.00538182259
+ (PID.TID 0000.0001)     Wall clock time:  2.80427932739257813E-003
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
- (PID.TID 0000.0001)   Seconds in section "TS_CORRECTION_STEP  [FORWARD_STEP]":
+ (PID.TID 0000.0001)   Seconds in section "TRC_CORRECTION_STEP [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.00399999321
+ (PID.TID 0000.0001)           User time:   0.0000000000000000
- (PID.TID 0000.0001)         System time:  0.
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.00273919106
+ (PID.TID 0000.0001)     Wall clock time:  1.46293640136718750E-003
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
  (PID.TID 0000.0001)   Seconds in section "BLOCKING_EXCHANGES  [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.000999003649
+ (PID.TID 0000.0001)           User time:   0.0000000000000000
- (PID.TID 0000.0001)         System time:  0.
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.00204181671
+ (PID.TID 0000.0001)     Wall clock time:  1.75905227661132813E-003
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
- (PID.TID 0000.0001)   Seconds in section "DO_STATEVARS_TAVE   [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.
- (PID.TID 0000.0001)         System time:  0.
- (PID.TID 0000.0001)     Wall clock time:  0.000355958939
- (PID.TID 0000.0001)          No. starts: 10
- (PID.TID 0000.0001)           No. stops: 10
  (PID.TID 0000.0001)   Seconds in section "MONITOR             [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.0429939628
+ (PID.TID 0000.0001)           User time:  1.00000000000000089E-002
- (PID.TID 0000.0001)         System time:  0.00100000016
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.0436413288
+ (PID.TID 0000.0001)     Wall clock time:  1.37252807617187500E-002
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
  (PID.TID 0000.0001)   Seconds in section "DO_THE_MODEL_IO     [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.0129990131
+ (PID.TID 0000.0001)           User time:   0.0000000000000000
- (PID.TID 0000.0001)         System time:  0.0109969992
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.113498926
+ (PID.TID 0000.0001)     Wall clock time:  5.23114204406738281E-003
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
- (PID.TID 0000.0001)   Seconds in section "WRITE_CHECKPOINT    [FORWARD_STEP]":
+ (PID.TID 0000.0001)   Seconds in section "DO_WRITE_PICKUP     [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.000999994576
+ (PID.TID 0000.0001)           User time:  9.99999999999998113E-003
- (PID.TID 0000.0001)         System time:  0.
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.000370264053
+ (PID.TID 0000.0001)     Wall clock time:  1.31797790527343750E-003
- (PID.TID 0000.0001)          No. starts: 10
+ (PID.TID 0000.0001)          No. starts:          10
- (PID.TID 0000.0001)           No. stops: 10
+ (PID.TID 0000.0001)           No. stops:          10
- (PID.TID 0000.0001)   Seconds in section "WRITE_CHECKPOINT       [THE_MODEL_MAIN]":
- (PID.TID 0000.0001)           User time:  0.00899800658
- (PID.TID 0000.0001)         System time:  0.00600000098
- (PID.TID 0000.0001)     Wall clock time:  0.0176517963
- (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) // ======================================================
-Line 1848
+Line 2141
  (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 Tile number: 000005
+ (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: 000006
+ (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: 000007
+ (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: 000008
+ (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 =           5650
+ (PID.TID 0000.0001) //            No. barriers =           2810
  (PID.TID 0000.0001) //      Max. barrier spins =              1
  (PID.TID 0000.0001) //      Min. barrier spins =              1
- (PID.TID 0000.0001) //     Total barrier spins =           5650
+ (PID.TID 0000.0001) //     Total barrier spins =           2810
  (PID.TID 0000.0001) //      Avg. barrier spins =       1.00E+00
+ PROGRAM MAIN: Execution ended Normally

 Legend:



Removed from v.1.4
 


changed lines


 
Added in v.1.6
 Legend:



Removed from v.1.4
 


changed lines


 
Added in v.1.6
-Removed from v.1.4
+Added in v.1.6

	ViewVC Help
Powered by ViewVC 1.1.22