/[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.5 by jmc,
Tue Dec  5 20:55:00 2006 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:  checkpoint58n_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 Aug 23 11:34:15 EDT 2006
+ (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 21
  (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 31
+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 90
+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 127
+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 146
+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 181
+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 202
+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 216
+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 239
+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 257
+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)     ;
-Line 273
+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) SET_PARMS: done
  (PID.TID 0000.0001) Enter INI_VERTICAL_GRID: setInterFDr=    T ; setCenterDr=    F
- (PID.TID 0000.0001) %MON XC_max                       =   1.5250000000000E+06
+ (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 352
+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)  MDSREADFIELD: opening global file: topog.box
+ (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 369
+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 384
+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 399
+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 414
+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 441
+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) AdamsBashforthGt =   /* apply Adams-Bashforth extrapolation on Gt */
+ (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) 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 */
-Line 456
+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) AdamsBashforthGs =   /* apply Adams-Bashforth extrapolation on Gs */
+ (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) 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=  /* fluid 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 503
+Line 589
  (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 554
+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 584
+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)                 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)     ;
-Line 636
+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) uniformLin_PhiSurf =   /* use uniform Bo_surf on/off flag*/
+ (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)                   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) 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) 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 */
-Line 685
+Line 812
  (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)     ;
-Line 700
+Line 833
  (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 */
-Line 712
+Line 845
  (PID.TID 0000.0001) useNHMTerms = /* Non-Hydrostatic Metric-Terms on/off */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) useConstantF =  /* use Constant f0 Coriolis flag */
+ (PID.TID 0000.0001) selectCoriMap = /* Coriolis Map options (0,1,2,3)*/
- (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)                       1
- (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001)     0= f-Plane ; 1= Beta-Plane ; 2= Spherical ; 3= read from file
- (PID.TID 0000.0001) useBetaPlaneF = /* use Beta-Plane Coriolis flag */
- (PID.TID 0000.0001)                   T
- (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) useSphereF  =   /* use Spherical Coriolis flag */
- (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) use3dCoriolis = /* 3-D Coriolis on/off flag */
  (PID.TID 0000.0001)                   F
-Line 730
+Line 858
  (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 772
+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 781
+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) saltIsActiveTr = /* Salt  is a dynamically Active Tracer */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001)  readBinaryPrec = /* Precision used for reading binary files */
  (PID.TID 0000.0001)                      32
  (PID.TID 0000.0001)     ;
-Line 806
+Line 948
  (PID.TID 0000.0001)                      64
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001)  globalFiles = /* write "global" (=not per tile) files */
- (PID.TID 0000.0001)                   T
+ (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)  debugMode  =  /* Debug Mode on/off flag */
+ (PID.TID 0000.0001) /* debLev[*]  : level of debug & auxiliary message printing */
- (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001) debLevZero =  0 ; /* level of disabled aux. msg printing */
- (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001)    debLevA =  1 ; /* level of minimum  aux. msg printing */
- (PID.TID 0000.0001)    debLevA  =  /* 1rst level of debugging */
+ (PID.TID 0000.0001)    debLevB =  2 ; /* level of low aux. print (report read-file opening)*/
- (PID.TID 0000.0001)                       1
+ (PID.TID 0000.0001)    debLevC =  3 ; /* level of moderate debug prt (most pkgs debug msg) */
- (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001)    debLevD =  4 ; /* level of enhanced debug prt (add DEBUG_STATS prt) */
- (PID.TID 0000.0001)    debLevB  =  /* 2nd  level of debugging */
+ (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)  debugLevel =  /* select debugging level */
- (PID.TID 0000.0001)                       1
- (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) //
  (PID.TID 0000.0001) // Elliptic solver(s) paramters ( PARM02 in namelist )
  (PID.TID 0000.0001) //
-Line 841
+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)                  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) deltaTmom =   /* 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) deltaTClock  =   /* Model clock timestep ( s ) */
-Line 880
+Line 1020
  (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 904
+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 916
+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 932
+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 )
-Line 949
+Line 1107
  (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) gravitySign = /* gravity 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) mass2rUnit = /* convert mass per unit area [kg/m2] to r-units [m] */
- (PID.TID 0000.0001)                 1.000000000000000E+00
+ (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 967
+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)                 6.370000000000000E+06
-Line 985
+Line 1152
  (PID.TID 0000.0001) deepAtmosphere = /* Deep/Shallow Atmosphere flag (True/False) */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) xcoord =   /* P-point X coord ( m - cartesian, degrees - spherical ) */
+ (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 1016
+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 ) */
-Line 1064
+Line 1233
  (PID.TID 0000.0001) deepFacF = /* deep-model grid factor @ W-Interface (-) */
  (PID.TID 0000.0001)     2 @  1.000000000000000E+00              /* K =  1:  2 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dBdrRef = /* Vertical gradient of reference boyancy [(m/s/r)^2)] */
+ (PID.TID 0000.0001) rVel2wUnit = /* convert units: rVel -> wSpeed (=1 if z-coord)*/
+ (PID.TID 0000.0001)     2 @  1.000000000000000E+00              /* K =  1:  2 */
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) wUnit2rVel = /* convert units: wSpeed -> rVel (=1 if z-coord)*/
+ (PID.TID 0000.0001)     2 @  1.000000000000000E+00              /* K =  1:  2 */
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) dBdrRef = /* Vertical grad. of reference buoyancy [(m/s/r)^2] */
  (PID.TID 0000.0001)                 0.000000000000000E+00       /* K =  1 */
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) rotateGrid = /* use rotated grid ( True/False ) */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) phiEuler = /* Euler angle, rotation about original z-coordinate [rad] */
+ (PID.TID 0000.0001)                 0.000000000000000E+00
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) thetaEuler = /* Euler angle, rotation about new x-coordinate [rad] */
+ (PID.TID 0000.0001)                 0.000000000000000E+00
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) psiEuler = /* Euler angle, rotation about new z-coordinate [rad] */
+ (PID.TID 0000.0001)                 0.000000000000000E+00
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) dxF =  /* dxF(:,1,:,1) ( units: m ) */
- (PID.TID 0000.0001)    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) 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              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
  (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              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
  (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              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
  (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              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
  (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              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
  (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              /* I =  1: 32 */
  (PID.TID 0000.0001)     ;
  (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              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) dxC =  /* dxC(:,1,:,1) ( units: m ) */
- (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) 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              /* J =  1: 32 */
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) dyC =  /* dyC(:,1,:,1) ( units: m ) */
- (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) dyC =  /* dyC(1,:,1,:) ( units: m ) */
- (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) dxV =  /* dxV(:,1,:,1) ( units: m ) */
- (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) dxV =  /* dxV(1,:,1,:) ( units: m ) */
- (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) dyU =  /* dyU(:,1,:,1) ( units: m ) */
- (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) dyU =  /* dyU(1,:,1,:) ( units: m ) */
- (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) rA  =  /* rA (:,1,:,1) ( units: m^2 ) */
- (PID.TID 0000.0001)    31 @  2.500000000000000E+09              /* I =  1: 31 */
+ (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)    31 @  2.500000000000000E+09              /* J =  1: 31 */
+ (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)    31 @  2.500000000000000E+09              /* I =  1: 31 */
+ (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)    31 @  2.500000000000000E+09              /* J =  1: 31 */
+ (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)    31 @  2.500000000000000E+09              /* I =  1: 31 */
+ (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)    31 @  2.500000000000000E+09              /* J =  1: 31 */
+ (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
-Line 1140
+Line 1327
  (PID.TID 0000.0001) // End of Model config. summary
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001)
- (PID.TID 0000.0001)  MDSREADFIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001) == Packages configuration : Check & print summary ==
- (PID.TID 0000.0001)  MDSREADFIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)
- (PID.TID 0000.0001)  MDSREADFIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001) GAD_CHECK: #define ALLOW_GENERIC_ADVDIFF
- (PID.TID 0000.0001)  MDSREADFIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001) // =======================================================
- (PID.TID 0000.0001)  MDSREADFIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001) // CONFIG_CHECK : Normal End
- (PID.TID 0000.0001)  MDSREADFIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001) // =======================================================
- (PID.TID 0000.0001)  MDSREADFIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)
- (PID.TID 0000.0001)  MDSREADFIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)  MDS_READ_META: opening file: pickup.0000200000.meta
- (PID.TID 0000.0001)  MDSREADFIELD: opening global file: pickup.0000200000.data
+ (PID.TID 0000.0001)  nRecords =  11 ; filePrec =  64 ; fileIter =    200000
- (PID.TID 0000.0001)  MDSREADFIELD: opening global file: trinit1
+ (PID.TID 0000.0001)     nDims =   2 , dims:
- (PID.TID 0000.0001)  MDSREADFIELD: opening global file: windx.cos_y
+ (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 1162
+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 1190
+Line 1392
  (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 dynstat_sst_max              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_min              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_mean             =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sst_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_max              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_min              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_mean             =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_max          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_min          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_mean         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_sd           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_del2         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_max           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_min           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_mean          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_sd            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_del2          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_max         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_min         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_mean        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_sd          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_del2        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fu_max            =   4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_min            =  -4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_mean           =  -8.5379215011820E-03
- (PID.TID 0000.0001) %MON extforcing_fu_sd             =   3.5065777661936E-01
- (PID.TID 0000.0001) %MON extforcing_fu_del2           =   8.8756883092094E-04
- (PID.TID 0000.0001) %MON extforcing_fv_max            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_min            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_mean           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_sd             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_del2           =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347586549727E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079091217551E-02
  (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.0006153255701E-09
-Line 1236
+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 1244
+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.18737609145978E-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.80030121840307E-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 1255
+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.3474823364103E-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.0969823047232E-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.1260621643150E-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.4764108425006E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.4764106171206E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -4.1070490460843E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -4.1070489119738E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -3.6429116802651E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   5.6690808732883E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.4198835094896E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.4198834566473E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0842430221123E-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 1283
+Line 1451
  (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 dynstat_sst_max              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_min              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_mean             =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sst_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_max              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_min              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_mean             =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_max          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_min          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_mean         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_sd           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_del2         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_max           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_min           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_mean          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_sd            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_del2          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_max         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_min         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_mean        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_sd          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_del2        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fu_max            =   4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_min            =  -4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_mean           =  -8.5379215011820E-03
- (PID.TID 0000.0001) %MON extforcing_fu_sd             =   3.5065777661936E-01
- (PID.TID 0000.0001) %MON extforcing_fu_del2           =   8.8756883092094E-04
- (PID.TID 0000.0001) %MON extforcing_fv_max            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_min            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_mean           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_sd             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_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.2856392368674E-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 1329
+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         =  -7.2858233592530E-20
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.1338161746577E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -3.6429116796265E-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.29182940352864E-05  7.14831409113866E-01
+  cg2d: Sum(rhs),rhsMax =   1.29182940431516E-05  7.14831409113866E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.17016481019693E-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.01765683373663E-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 1350
+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.3472844558261E-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.2511462253914E-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.5599180900709E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.3522319386033E-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.4668357335031E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.4668355248868E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.2037681229413E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.2037683585659E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   5.4043510948307E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -2.6053683253500E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   7.7071975512831E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   7.7071974717058E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.4897450760562E-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 1378
+Line 1511
  (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 dynstat_sst_max              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_min              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_mean             =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sst_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_max              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_min              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_mean             =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_max          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_min          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_mean         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_sd           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_del2         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_max           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_min           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_mean          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_sd            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_del2          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_max         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_min         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_mean        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_sd          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_del2        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fu_max            =   4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_min            =  -4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_mean           =  -8.5379215011820E-03
- (PID.TID 0000.0001) %MON extforcing_fu_sd             =   3.5065777661936E-01
- (PID.TID 0000.0001) %MON extforcing_fu_del2           =   8.8756883092094E-04
- (PID.TID 0000.0001) %MON extforcing_fv_max            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_min            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_mean           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_sd             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_del2           =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347585314096E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093466561E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.5630144983530E-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.7522710928396E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004696744475E-02
-Line 1424
+Line 1522
  (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.0808702170583E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -5.2107366506999E-20
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   5.4043510852915E-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.29179873122401E-05  7.14831436211394E-01
+  cg2d: Sum(rhs),rhsMax =   1.29179873209484E-05  7.14831436211394E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.16678717243530E-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.60155251655069E-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.2594517347064E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594517347063E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3475494279147E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3475494280297E-06
  (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.1735243795035E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   5.1237305261427E-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.2601410072550E-10
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -1.2601408525219E-10
  (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.8420639745891E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   3.8420640602708E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.6656310558319E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.6656308397651E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -9.0801484598699E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -4.8013098832578E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.2298095642515E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.2298093646288E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.1228957716442E-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 1473
+Line 1571
  (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 dynstat_sst_max              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_min              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_mean             =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sst_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_max              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_min              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_mean             =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_max          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_min          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_mean         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_sd           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_del2         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_max           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_min           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_mean          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_sd            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_del2          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_max         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_min         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_mean        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_sd          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_del2        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fu_max            =   4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_min            =  -4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_mean           =  -8.5379215011820E-03
- (PID.TID 0000.0001) %MON extforcing_fu_sd             =   3.5065777661936E-01
- (PID.TID 0000.0001) %MON extforcing_fu_del2           =   8.8756883092094E-04
- (PID.TID 0000.0001) %MON extforcing_fv_max            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_min            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_mean           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_sd             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_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.0736511796713E-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.7522696832281E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004696439068E-02
-Line 1519
+Line 1582
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759247736E-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.8160296919740E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -9.6026197665155E-20
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -9.0801484598699E-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.29183979198783E-05  7.14831406464009E-01
+  cg2d: Sum(rhs),rhsMax =   1.29183979185044E-05  7.14831406464009E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.00714383795439E-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.27800113850331E-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) // =======================================================
-Line 1540
+Line 1603
  (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.2594517219516E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472375531510E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472375532881E-06
- (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574660860786E-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.1735249377186E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -2.4376814494190E-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.5599188258601E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -1.9581025602098E-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.9360995553434E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.9361001253128E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -4.0912349000573E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -4.0912351571023E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -4.9279197671656E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -9.9820383471874E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.9178890259297E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.9178891520071E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.9047004703512E-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 1568
+Line 1631
  (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 dynstat_sst_max              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_min              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_mean             =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sst_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_max              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_min              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_mean             =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_max          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_min          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_mean         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_sd           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_del2         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_max           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_min           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_mean          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_sd            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_del2          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_max         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_min         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_mean        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_sd          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_del2        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fu_max            =   4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_min            =  -4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_mean           =  -8.5379215011820E-03
- (PID.TID 0000.0001) %MON extforcing_fu_sd             =   3.5065777661936E-01
- (PID.TID 0000.0001) %MON extforcing_fu_del2           =   8.8756883092094E-04
- (PID.TID 0000.0001) %MON extforcing_fv_max            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_min            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_mean           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_sd             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_del2           =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347584056695E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093064320E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   3.2729879200459E-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.7522688592639E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004696598940E-02
- (PID.TID 0000.0001) %MON ke_mean                      =   9.4023111519788E-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.4130473004107E-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         =  -9.8558395744748E-20
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -1.9964076694375E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -4.9279197872374E-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.29179148317393E-05  7.14831431051197E-01
+  cg2d: Sum(rhs),rhsMax =   1.29179148262715E-05  7.14831431051196E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.17107070057593E-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.97646073430777E-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.1785902732705E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594507560530E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594507560532E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3475311734807E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3475311736921E-06
  (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.0868961106712E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735243747491E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -3.7239070551786E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -3.7239074579678E-10
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548180547E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408278058110E-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.1283018928394E-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.2239052616060E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.2239063447341E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.9653511904180E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.9653501883149E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   7.6327672767378E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   1.7386250129485E-20
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   6.7069706964106E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   6.7069706534392E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   2.9393667424149E-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 1663
+Line 1691
  (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 dynstat_sst_max              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_min              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_mean             =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sst_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_max              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_min              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_mean             =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_max          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_min          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_mean         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_sd           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_del2         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_max           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_min           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_mean          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_sd            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_del2          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_max         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_min         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_mean        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_sd          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_del2        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fu_max            =   4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_min            =  -4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_mean           =  -8.5379215011820E-03
- (PID.TID 0000.0001) %MON extforcing_fu_sd             =   3.5065777661936E-01
- (PID.TID 0000.0001) %MON extforcing_fu_del2           =   8.8756883092094E-04
- (PID.TID 0000.0001) %MON extforcing_fv_max            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_min            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_mean           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_sd             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_del2           =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347584442685E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093939875E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093939876E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   1.7791242092848E-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.7522719925861E-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.4023105216738E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.1062500000000E+16
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759598706E-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         =   1.5265534574029E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   3.4772500258970E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   7.6327672870146E-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.29183692763792E-05  7.14831428086830E-01
+  cg2d: Sum(rhs),rhsMax =   1.29183692794133E-05  7.14831428086836E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.28565909704110E-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.56127649303186E-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) // =======================================================
-Line 1730
+Line 1723
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0001200000000E+09
  (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785900288036E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594511076654E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3470830988358E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3470830991681E-06
  (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.0868966158535E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735239645687E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   3.2506628225280E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   3.2506624141082E-10
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548296738E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408277509125E-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.8431378272519E-11
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   9.8431432383041E-11
  (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.9895875416696E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.9895877633244E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.8704261407256E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.8704262655228E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.3276361734547E-20
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   9.1993590648481E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.2304091185756E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.2304090930841E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.6659607108016E-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 1758
+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 dynstat_sst_max              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_min              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_mean             =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sst_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_max              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_min              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_mean             =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_max          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_min          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_mean         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_sd           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_del2         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_max           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_min           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_mean          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_sd            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_del2          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_max         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_min         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_mean        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_sd          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_del2        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fu_max            =   4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_min            =  -4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_mean           =  -8.5379215011820E-03
- (PID.TID 0000.0001) %MON extforcing_fu_sd             =   3.5065777661936E-01
- (PID.TID 0000.0001) %MON extforcing_fu_del2           =   8.8756883092094E-04
- (PID.TID 0000.0001) %MON extforcing_fv_max            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_min            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_mean           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_sd             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_del2           =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347586463414E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079095296994E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.3916700333357E-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.7522689154176E-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.4023106862759E-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.4130473019324E-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         =  -2.6552723590489E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.8398718129696E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -1.3276361795244E-19
+ (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.29176754176067E-05  7.14831452448172E-01
+  cg2d: Sum(rhs),rhsMax =   1.29176754160021E-05  7.14831452448173E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.00898122153472E-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.56122073215609E-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) // =======================================================
-Line 1825
+Line 1783
  (PID.TID 0000.0001) %MON time_secondsf                =   4.0001400000000E+09
  (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785900384355E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594517423048E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3473803444489E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3473803447387E-06
  (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.0868959295820E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868959295821E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735246434860E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735246434859E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   1.6257393654592E-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.4408278720511E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326095407372E-05
  (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.1358487064641E-09
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -2.1358486685855E-09
  (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.9678630866110E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.9678631760180E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.2316897902638E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.2316893395036E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   2.4088764448761E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -6.9450454176062E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.2289378318267E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.2289378984835E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.3867028240911E-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 1853
+Line 1811
  (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 dynstat_sst_max              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_min              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_mean             =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sst_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_max              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_min              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_mean             =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_max          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_min          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_mean         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_sd           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_del2         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_max           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_min           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_mean          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_sd            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_del2          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_max         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_min         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_mean        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_sd          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_del2        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fu_max            =   4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_min            =  -4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_mean           =  -8.5379215011820E-03
- (PID.TID 0000.0001) %MON extforcing_fu_sd             =   3.5065777661936E-01
- (PID.TID 0000.0001) %MON extforcing_fu_del2           =   8.8756883092094E-04
- (PID.TID 0000.0001) %MON extforcing_fv_max            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_min            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_mean           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_sd             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_del2           =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347583718328E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079092310019E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.5853518322110E-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.7522696504045E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004696211241E-02
-Line 1899
+Line 1822
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.3833759310686E-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         =   4.8177527253237E-20
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -1.3890090835212E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   2.4088763626619E-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.29181366819893E-05  7.14831382145344E-01
+  cg2d: Sum(rhs),rhsMax =   1.29181366901321E-05  7.14831382145344E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.06514671227263E-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.65957983915145E-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.1785900968546E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   2.1785900968545E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594503936919E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594503936921E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472411337467E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472411339237E-06
  (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.0868963218162E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868963218163E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735245595629E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -4.2115283125600E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -4.2115275950292E-10
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065549284519E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408284752086E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326096277653E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773800541E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599182940112E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -9.2460704000194E-11
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -9.2460769736842E-11
  (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.6578598152846E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.6578603666276E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.4862337270752E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.4862335808575E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.0724021595413E-20
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   7.6360558554278E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   7.3628383901427E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   7.3628384678627E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.1854211698628E-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 1948
+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 dynstat_sst_max              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_min              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_mean             =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sst_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_max              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_min              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_mean             =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_max          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_min          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_mean         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_sd           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_del2         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_max           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_min           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_mean          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_sd            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_del2          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_max         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_min         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_mean        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_sd          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_del2        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fu_max            =   4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_min            =  -4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_mean           =  -8.5379215011820E-03
- (PID.TID 0000.0001) %MON extforcing_fu_sd             =   3.5065777661936E-01
- (PID.TID 0000.0001) %MON extforcing_fu_del2           =   8.8756883092094E-04
- (PID.TID 0000.0001) %MON extforcing_fv_max            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_min            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_mean           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_sd             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_del2           =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347585287265E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079095202163E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.1262878522277E-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.7522705426052E-06
+ (PID.TID 0000.0001) %MON pe_b_mean                    =   5.7522705426051E-06
  (PID.TID 0000.0001) %MON ke_max                       =   1.4004697375029E-02
- (PID.TID 0000.0001) %MON ke_mean                      =   9.4023110830423E-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.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.4130473004553E-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         =  -2.1448043230970E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.5272111710856E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =  -1.0724021615485E-19
+ (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.29179200117433E-05  7.14831450919334E-01
+  cg2d: Sum(rhs),rhsMax =   1.29179200107910E-05  7.14831450919334E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    1.11079791354847E-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.08082503103431E-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.1785900450429E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594527421719E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594527421720E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3476472944075E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3476472946893E-06
- (PID.TID 0000.0001) %MON dynstat_eta_sd               =   7.6574670300973E-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.0868964909379E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868964909380E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735244831868E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   3.1056225620952E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   3.1056227065058E-10
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065548419668E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408284028509E-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.5599182338730E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.5253257552387E-10
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =   1.5253261266447E-10
  (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.6066183432995E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   2.6066179811896E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.9599356590770E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -3.9599354909733E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   5.1877166718076E-21
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   9.8743377316554E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.7292355947767E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   8.7292357087263E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   3.9402321240493E-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 2043
+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 dynstat_sst_max              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_min              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_mean             =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sst_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_max              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_min              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_mean             =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_max          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_min          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_mean         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_sd           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_del2         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_max           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_min           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_mean          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_sd            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_del2          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_max         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_min         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_mean        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_sd          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_del2        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fu_max            =   4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_min            =  -4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_mean           =  -8.5379215011820E-03
- (PID.TID 0000.0001) %MON extforcing_fu_sd             =   3.5065777661936E-01
- (PID.TID 0000.0001) %MON extforcing_fu_del2           =   8.8756883092094E-04
- (PID.TID 0000.0001) %MON extforcing_fv_max            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_min            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_mean           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_sd             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_del2           =   0.0000000000000E+00
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347585963752E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079093571044E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   3.1679485272616E-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.7522702775595E-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.4023105744071E-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.9603605531574E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =   1.0735859124867E-04
- (PID.TID 0000.0001) %MON vort_a_sd                    =   4.4130473014548E-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         =   1.0375433343615E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   1.9748675463311E-19
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   5.1877166718076E-20
+ (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.29185494260192E-05  7.14831404157640E-01
+  cg2d: Sum(rhs),rhsMax =   1.29185494190942E-05  7.14831404157643E-01
- (PID.TID 0000.0001)                    cg2d_init_res =    9.57340498015102E-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.79891798066715E-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.1785899418361E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594507089043E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -9.2594507089044E-02
- (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472039449101E-06
+ (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -8.3472039451370E-06
  (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.0868961191918E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   7.0868961191919E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735249175663E-02
+ (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -5.1735249175662E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   2.3555323438000E-10
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =   2.3555353362573E-10
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   2.3065549206510E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   6.4408284343721E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   8.2326094458559E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.6769773535244E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -5.5599188270228E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -2.7314099682717E-09
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -2.7314097624269E-09
- (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608826082411E-02
+ (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   3.7608826082412E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   1.6259700240885E-03
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   4.3373808400809E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   3.1423902139068E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   3.1423902697861E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.8826322350651E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -2.8826323840767E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   1.6364327484810E-20
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.5300064542368E-21
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.4185795126640E-11
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   9.4185792767575E-11
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.2311135596787E-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 2138
+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 dynstat_sst_max              =   2.0000000000000E+01
+ (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347584476768E-02
- (PID.TID 0000.0001) %MON dynstat_sst_min              =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_mean             =   2.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sst_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sst_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_max              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_min              =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_mean             =   1.0000000000000E+01
- (PID.TID 0000.0001) %MON dynstat_sss_sd               =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON dynstat_sss_del2             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_max          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_min          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_mean         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_sd           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qnet_del2         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_max           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_min           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_mean          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_sd            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_qsw_del2          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_max         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_min         =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_mean        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_sd          =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_empmr_del2        =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fu_max            =   4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_min            =  -4.9931475520134E-01
- (PID.TID 0000.0001) %MON extforcing_fu_mean           =  -8.5379215011820E-03
- (PID.TID 0000.0001) %MON extforcing_fu_sd             =   3.5065777661936E-01
- (PID.TID 0000.0001) %MON extforcing_fu_del2           =   8.8756883092094E-04
- (PID.TID 0000.0001) %MON extforcing_fv_max            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_min            =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_mean           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_sd             =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON extforcing_fv_del2           =   0.0000000000000E+00
- (PID.TID 0000.0001) %MON advcfl_uvel_max              =   2.8347584476767E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   6.7079094140976E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   2.5139121711254E-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.7522690130739E-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.4023111433233E-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.4130473016938E-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         =   3.2728654970101E-19
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =  -3.0600129084735E-20
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.6364327485050E-19
+ (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.378941991
+ (PID.TID 0000.0001)           User time:  0.14999999999999999
- (PID.TID 0000.0001)         System time:  0.021996999
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.566043139
+ (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.0539919997
+ (PID.TID 0000.0001)           User time:  2.00000000000000004E-002
- (PID.TID 0000.0001)         System time:  0.00400000007
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.101105928
+ (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.318950992
+ (PID.TID 0000.0001)           User time:  0.13000000000000000
- (PID.TID 0000.0001)         System time:  0.0149970003
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.425887823
+ (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.00699899718
+ (PID.TID 0000.0001)           User time:  9.99999999999999847E-003
- (PID.TID 0000.0001)         System time:  0.00199900009
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.0115380287
+ (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 "MAIN LOOP           [THE_MAIN_LOOP]":
- (PID.TID 0000.0001)           User time:  0.302953988
+ (PID.TID 0000.0001)           User time:  0.12000000000000000
- (PID.TID 0000.0001)         System time:  0.0109980004
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.377712965
+ (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.302953988
+ (PID.TID 0000.0001)           User time:  0.11999999999999997
- (PID.TID 0000.0001)         System time:  0.00999800116
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.377342701
+ (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 "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.00101208687
+ (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.
+ (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.000341892242
+ (PID.TID 0000.0001)     Wall clock time:  9.84668731689453125E-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_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.000334978104
+ (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.00699800253
+ (PID.TID 0000.0001)           User time:   0.0000000000000000
- (PID.TID 0000.0001)         System time:  0.00100000016
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.00552582741
+ (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.052993983
+ (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.0564746857
+ (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.0579909682
+ (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.0564823151
+ (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.0849860311
+ (PID.TID 0000.0001)           User time:  4.99999999999999889E-002
- (PID.TID 0000.0001)         System time:  0.
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.0847268105
+ (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 "MOM_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.00559735298
+ (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 "TRC_CORRECTION_STEP [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.00299900025
+ (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.00264954567
+ (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.
+ (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.00196218491
+ (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 "MONITOR             [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.071989052
+ (PID.TID 0000.0001)           User time:  1.00000000000000089E-002
- (PID.TID 0000.0001)         System time:  0.00199900009
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.0741562843
+ (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.015997991
+ (PID.TID 0000.0001)           User time:   0.0000000000000000
- (PID.TID 0000.0001)         System time:  0.00499900058
+ (PID.TID 0000.0001)         System time:   0.0000000000000000
- (PID.TID 0000.0001)     Wall clock time:  0.0834841728
+ (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 "DO_WRITE_PICKUP     [FORWARD_STEP]":
- (PID.TID 0000.0001)           User time:  0.
+ (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.000344276428
+ (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 "DO_WRITE_PICKUP        [THE_MODEL_MAIN]":
- (PID.TID 0000.0001)           User time:  0.00599899888
- (PID.TID 0000.0001)         System time:  0.00299999863
- (PID.TID 0000.0001)     Wall clock time:  0.0388941765
- (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 2329
+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 =           8136
+ (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 =           8136
+ (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.5
 


changed lines


 
Added in v.1.6
 Legend:



Removed from v.1.5
 


changed lines


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

	ViewVC Help
Powered by ViewVC 1.1.22