/[MITgcm]/MITgcm/verification/global_ocean.cs32x15/results/output.thsice.txt

Diff of /MITgcm/verification/global_ocean.cs32x15/results/output.thsice.txt

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

-revision 1.13 by jmc,
Wed Nov  3 18:23:31 2010 UTC
+revision 1.14 by jmc,
Mon Nov  7 21:59:40 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:  checkpoint62m
+ (PID.TID 0000.0001) // MITgcmUV version:  checkpoint63e
  (PID.TID 0000.0001) // Build user:        jmc
  (PID.TID 0000.0001) // Build host:        baudelaire
- (PID.TID 0000.0001) // Build date:        Wed Nov  3 03:25:03 EDT 2010
+ (PID.TID 0000.0001) // Build date:        Mon Nov  7 15:50:35 EST 2011
  (PID.TID 0000.0001)
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Execution Environment parameter file "eedata"
 Line 52
  (PID.TID 0000.0001)                   /*  "mpirun -np 64 ......"                    */
  (PID.TID 0000.0001) useCoupler=    F ;/* Flag used to control communications with   */
  (PID.TID 0000.0001)                   /*  other model components, through a coupler */
+ (PID.TID 0000.0001) debugMode =    F ; /* print debug msg. (sequence of S/R calls)  */
  (PID.TID 0000.0001) printMapIncludesZeros=    F ; /* print zeros in Std.Output maps */
  (PID.TID 0000.0001) maxLengthPrt1D=   65 /* maxLength of 1D array printed to StdOut */
  (PID.TID 0000.0001)
-Line 69
+Line 70
  (PID.TID 0000.0001)  write to log-file: w2_tile_topology.0000.log
  (PID.TID 0000.0001) =====       setting W2 TOPOLOGY: Done
  (PID.TID 0000.0001)
- (PID.TID 0000.0001) DEBUG_MSG: ENTERED S/R THE_MODEL_MAIN
+ (PID.TID 0000.0001)  INI_PARMS: opening model parameter file "data"
- (PID.TID 0000.0001) DEBUG_MSG: CALLING S/R INITIALISE_FIXED
+ (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 165
+Line 166
  (PID.TID 0000.0001) > saltClimFile   ='lev_surfS_cs_12m.bin',
  (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 210
+Line 212
  (PID.TID 0000.0001) >
  (PID.TID 0000.0001) > &GM_PARM01
  (PID.TID 0000.0001) >  GM_background_K    = 800.,
+ (PID.TID 0000.0001) ># GM_AdvForm         = .TRUE.,
+ (PID.TID 0000.0001) ># GM_background_K    = 10.,
+ (PID.TID 0000.0001) ># GM_Visbeck_alpha   = 1.e-2,
+ (PID.TID 0000.0001) >  GM_Visbeck_length  = 2.e+5,
+ (PID.TID 0000.0001) >  GM_Visbeck_depth   = 1.e+3,
+ (PID.TID 0000.0001) >  GM_Visbeck_maxval_K= 2.5e+3,
  (PID.TID 0000.0001) >  GM_taper_scheme    = 'gkw91',
+ (PID.TID 0000.0001) ># GM_Small_Number    = 1.e-18,
  (PID.TID 0000.0001) >  GM_maxSlope        = 1.e-2,
  (PID.TID 0000.0001) >  GM_Kmin_horiz      = 50.,
  (PID.TID 0000.0001) >  GM_Scrit           = 4.e-3,
  (PID.TID 0000.0001) >  GM_Sd              = 1.e-3,
- (PID.TID 0000.0001) ># GM_Visbeck_alpha   = 0.,
- (PID.TID 0000.0001) ># GM_Visbeck_length  = 2.e+5,
- (PID.TID 0000.0001) ># GM_Visbeck_depth   = 1.e+3,
- (PID.TID 0000.0001) ># GM_Visbeck_maxval_K= 2.5e+3,
  (PID.TID 0000.0001) > /
  (PID.TID 0000.0001) >
- (PID.TID 0000.0001) >
  (PID.TID 0000.0001)
  (PID.TID 0000.0001)  GM_READPARMS: finished reading data.gmredi
  (PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.blk
-Line 399
+Line 403
  (PID.TID 0000.0001) // Parameter file "data.diagnostics"
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) ># Diagnostic Package Choices
- (PID.TID 0000.0001) >#-----------------
+ (PID.TID 0000.0001) >#--------------------
- (PID.TID 0000.0001) ># for each output-stream:
+ (PID.TID 0000.0001) >#  dumpAtLast (logical): always write output at the end of simulation (default=F)
- (PID.TID 0000.0001) >#  filename(n) : prefix of the output file name (only 8.c long) for outp.stream n
+ (PID.TID 0000.0001) >#  diag_mnc   (logical): write to NetCDF files (default=useMNC)
+ (PID.TID 0000.0001) >#--for each output-stream:
+ (PID.TID 0000.0001) >#  fileName(n) : prefix of the output file name (max 80c long) for outp.stream n
  (PID.TID 0000.0001) >#  frequency(n):< 0 : write snap-shot output every |frequency| seconds
  (PID.TID 0000.0001) >#               > 0 : write time-average output every frequency seconds
  (PID.TID 0000.0001) >#  timePhase(n)     : write at time = timePhase + multiple of |frequency|
- (PID.TID 0000.0001) >#  averagingFreq(n) : frequency (in s) for periodic averaging interval
+ (PID.TID 0000.0001) >#    averagingFreq  : frequency (in s) for periodic averaging interval
- (PID.TID 0000.0001) >#  averagingPhase(n): phase     (in s) for periodic averaging interval
+ (PID.TID 0000.0001) >#    averagingPhase : phase     (in s) for periodic averaging interval
- (PID.TID 0000.0001) >#  repeatCycle(n)   : number of averaging intervals in 1 cycle
+ (PID.TID 0000.0001) >#    repeatCycle    : number of averaging intervals in 1 cycle
  (PID.TID 0000.0001) >#  levels(:,n) : list of levels to write to file (Notes: declared as REAL)
- (PID.TID 0000.0001) >#                 when this entry is missing, select all common levels of this list
+ (PID.TID 0000.0001) >#                when this entry is missing, select all common levels of this list
- (PID.TID 0000.0001) >#  fields(:,n) : list of diagnostics fields (8.c) (see "available_diagnostics.log"
+ (PID.TID 0000.0001) >#  fields(:,n) : list of selected diagnostics fields (8.c) in outp.stream n
- (PID.TID 0000.0001) >#                 file for the list of all available diag. in this particular config)
+ (PID.TID 0000.0001) >#                (see "available_diagnostics.log" file for the full list of diags)
- (PID.TID 0000.0001) >#-----------------
+ (PID.TID 0000.0001) >#  missing_value(n) : missing value for real-type fields in output file "n"
+ (PID.TID 0000.0001) >#  fileFlags(n)     : specific code (8c string) for output file "n"
+ (PID.TID 0000.0001) >#--------------------
  (PID.TID 0000.0001) > &DIAGNOSTICS_LIST
  (PID.TID 0000.0001) ># diag_mnc     = .FALSE.,
  (PID.TID 0000.0001) ># dumpAtLast   = .TRUE.,
- (PID.TID 0000.0001) >   fields(1,1) = 'ETAN    ','ETANSQ  ','DETADT2 ','PHIBOT  ','PHIBOTSQ',
+ (PID.TID 0000.0001) >#--
- (PID.TID 0000.0001) >                 'oceTAUX ','oceTAUY ','TRELAX  ','SRELAX  ',
+ (PID.TID 0000.0001) >  fields(1:17,1) = 'ETAN    ','ETANSQ  ','DETADT2 ','PHIBOT  ','PHIBOTSQ',
- (PID.TID 0000.0001) >                 'surForcT','surForcS','TFLUX   ','SFLUX   ',
+ (PID.TID 0000.0001) >                   'oceTAUX ','oceTAUY ','TRELAX  ','SRELAX  ',
- (PID.TID 0000.0001) >                 'oceQsw  ','oceQnet ','oceSflux','oceFWflx',
+ (PID.TID 0000.0001) >                   'surForcT','surForcS','TFLUX   ','SFLUX   ',
- (PID.TID 0000.0001) >   levels(1,1) = 1.,
+ (PID.TID 0000.0001) >                   'oceQsw  ','oceQnet ','oceSflux','oceFWflx',
- (PID.TID 0000.0001) >   filename(1) = 'surfDiag',
+ (PID.TID 0000.0001) >   fileName(1) = 'surfDiag',
  (PID.TID 0000.0001) >  frequency(1) = 1555200000.,
- (PID.TID 0000.0001) >   fields(1,2) = 'UVEL    ','VVEL    ','WVEL    ','PHIHYD  ',
+ (PID.TID 0000.0001) >
- (PID.TID 0000.0001) >                 'VVELMASS','UVELMASS','WVELSQ  ',
+ (PID.TID 0000.0001) >  fields(1:15,2) = 'UVEL    ','VVEL    ','WVEL    ','PHIHYD  ',
- (PID.TID 0000.0001) >                 'THETA   ','UTHMASS ','VTHMASS ','WTHMASS ',
+ (PID.TID 0000.0001) >                   'VVELMASS','UVELMASS','WVELSQ  ',
- (PID.TID 0000.0001) >                 'SALT    ','USLTMASS','VSLTMASS','WSLTMASS',
+ (PID.TID 0000.0001) >                   'THETA   ','UTHMASS ','VTHMASS ','WTHMASS ',
+ (PID.TID 0000.0001) >                   'SALT    ','USLTMASS','VSLTMASS','WSLTMASS',
  (PID.TID 0000.0001) ># do not specify levels => all levels are selected
- (PID.TID 0000.0001) >   filename(2) = 'dynDiag',
+ (PID.TID 0000.0001) >   fileName(2) = 'dynDiag',
  (PID.TID 0000.0001) >  frequency(2) = 1555200000.,
- (PID.TID 0000.0001) >   fields(1,3) = 'DRHODR  ','RHOAnoma','CONVADJ ',
+ (PID.TID 0000.0001) >
- (PID.TID 0000.0001) >                 'GM_PsiX ','GM_PsiY ',
+ (PID.TID 0000.0001) >  fields(1:12,3) = 'DRHODR  ','RHOAnoma','CONVADJ ',
- (PID.TID 0000.0001) >                 'GM_Kwx  ','GM_Kwy  ','GM_Kwz  ',
+ (PID.TID 0000.0001) >                   'GM_PsiX ','GM_PsiY ',
- (PID.TID 0000.0001) >                 'GM_Kux  ','GM_Kvy  ',
+ (PID.TID 0000.0001) >                   'GM_Kwx  ','GM_Kwy  ','GM_Kwz  ',
- (PID.TID 0000.0001) >                 'GM_Kuz  ','GM_Kvz  ',
+ (PID.TID 0000.0001) >                   'GM_Kux  ','GM_Kvy  ',
+ (PID.TID 0000.0001) >                   'GM_Kuz  ','GM_Kvz  ',
  (PID.TID 0000.0001) >#- disable this output list by commenting out the file name
- (PID.TID 0000.0001) >#  filename(3) = 'oceDiag',
+ (PID.TID 0000.0001) >#  fileName(3) = 'oceDiag',
  (PID.TID 0000.0001) >  frequency(3) = 1555200000.,
- (PID.TID 0000.0001) >   fields(1,4) = 'ADVx_TH ','ADVy_TH ','ADVr_TH ',
+ (PID.TID 0000.0001) >
- (PID.TID 0000.0001) >                 'DFxE_TH ','DFyE_TH ','DFrE_TH ',
+ (PID.TID 0000.0001) >  fields(1:7,4)  = 'ADVx_TH ','ADVy_TH ','ADVr_TH ',
- (PID.TID 0000.0001) >                 'DFrI_TH ',
+ (PID.TID 0000.0001) >                   'DIFx_TH ','DIFy_TH ','DFrE_TH ',
- (PID.TID 0000.0001) >#                'ADVx_SLT',
+ (PID.TID 0000.0001) >                   'DFrI_TH ',
- (PID.TID 0000.0001) >#  filename(4) = 'flxDiag',
+ (PID.TID 0000.0001) >#                  'ADVx_SLT',
+ (PID.TID 0000.0001) >#  fileName(4) = 'flxDiag',
  (PID.TID 0000.0001) >  frequency(4) = 1728000.,
- (PID.TID 0000.0001) >   fields(1,5) = 'SI_Fract','SI_Thick','SI_SnowH',
+ (PID.TID 0000.0001) >
- (PID.TID 0000.0001) >                 'SI_Tsrf ','SI_Tice1','SI_Tice2',
+ (PID.TID 0000.0001) >  fields(1:16,5) = 'SI_Fract','SI_Thick','SI_SnowH',
- (PID.TID 0000.0001) >                 'SI_Qice1','SI_Qice2','SIsnwAge',
+ (PID.TID 0000.0001) >                   'SI_Tsrf ','SI_Tice1','SI_Tice2',
- (PID.TID 0000.0001) >                 'SIsnwPrc','SIalbedo',
+ (PID.TID 0000.0001) >                   'SI_Qice1','SI_Qice2','SIsnwAge',
- (PID.TID 0000.0001) >                 'SIflx2oc','SIfrw2oc','SIsaltFx',
+ (PID.TID 0000.0001) >                   'SIsnwPrc','SIalbedo',
- (PID.TID 0000.0001) >                 'SIflxAtm','SIfrwAtm',
+ (PID.TID 0000.0001) >                   'SIflx2oc','SIfrw2oc','SIsaltFx',
- (PID.TID 0000.0001) >#                'SItOcMxL','SIsOcMxL',
+ (PID.TID 0000.0001) >                   'SIflxAtm','SIfrwAtm',
- (PID.TID 0000.0001) >   filename(5) = 'thSIceDiag',
+ (PID.TID 0000.0001) >#                  'SItOcMxL','SIsOcMxL',
+ (PID.TID 0000.0001) >   fileName(5) = 'thSIceDiag',
  (PID.TID 0000.0001) >  frequency(5) = 1555200000.,
  (PID.TID 0000.0001) > averagingFreq(5) = 2592000.,
  (PID.TID 0000.0001) >   repeatCycle(5) = 12,
  (PID.TID 0000.0001) >#--- to check conservation of Volume, Heat & Salt (with z*+realFW+AB_2):
- (PID.TID 0000.0001) >   fields(1,6) = 'ETAN    ','oceFWflx',
+ (PID.TID 0000.0001) >  fields(1:2,6)  = 'ETAN    ','oceFWflx',
- (PID.TID 0000.0001) >   filename(6) = 'surfInst',
+ (PID.TID 0000.0001) >   fileName(6) = 'surfInst',
  (PID.TID 0000.0001) >  frequency(6) = -777600.,
  (PID.TID 0000.0001) >  timePhase(6) = 0.,
- (PID.TID 0000.0001) >   fields(1,7) = 'ETAN    ',
+ (PID.TID 0000.0001) >  fields(1,7)  = 'ETAN    ',
- (PID.TID 0000.0001) >   filename(7) = 'etaInst',
+ (PID.TID 0000.0001) >   fileName(7) = 'etaInst',
  (PID.TID 0000.0001) >  frequency(7) = -777600.,
  (PID.TID 0000.0001) >  timePhase(7) = 86400.,
- (PID.TID 0000.0001) >   fields(1,8) = 'THETA   ','SALT    ',
+ (PID.TID 0000.0001) >  fields(1:2,8)  = 'THETA   ','SALT    ',
- (PID.TID 0000.0001) >   filename(8) = 'dynInst',
+ (PID.TID 0000.0001) >   fileName(8) = 'dynInst',
  (PID.TID 0000.0001) >  frequency(8) = -777600.,
  (PID.TID 0000.0001) >  timePhase(8) = 0.,
  (PID.TID 0000.0001) > /
  (PID.TID 0000.0001) >
+ (PID.TID 0000.0001) >#--------------------
  (PID.TID 0000.0001) ># Parameter for Diagnostics of per level statistics:
- (PID.TID 0000.0001) >#-----------------
+ (PID.TID 0000.0001) >#--------------------
- (PID.TID 0000.0001) ># for each output-stream:
+ (PID.TID 0000.0001) >#  diagSt_mnc (logical): write stat-diags to NetCDF files (default=diag_mnc)
- (PID.TID 0000.0001) >#  stat_fname(n) : prefix of the output file name (only 8.c long) for outp.stream n
+ (PID.TID 0000.0001) >#  diagSt_regMaskFile : file containing the region-mask to read-in
+ (PID.TID 0000.0001) >#  nSetRegMskFile   : number of region-mask sets within the region-mask file
+ (PID.TID 0000.0001) >#  set_regMask(i)   : region-mask set-index that identifies the region "i"
+ (PID.TID 0000.0001) >#  val_regMask(i)   : region "i" identifier value in the region mask
+ (PID.TID 0000.0001) >#--for each output-stream:
+ (PID.TID 0000.0001) >#  stat_fName(n) : prefix of the output file name (max 80c long) for outp.stream n
  (PID.TID 0000.0001) >#  stat_freq(n):< 0 : write snap-shot output every |stat_freq| seconds
  (PID.TID 0000.0001) >#               > 0 : write time-average output every stat_freq seconds
  (PID.TID 0000.0001) >#  stat_phase(n)    : write at time = stat_phase + multiple of |stat_freq|
  (PID.TID 0000.0001) >#  stat_region(:,n) : list of "regions" (default: 1 region only=global)
- (PID.TID 0000.0001) >#  stat_fields(:,n) : list of diagnostics fields (8.c) (see "available_diagnostics.log"
+ (PID.TID 0000.0001) >#  stat_fields(:,n) : list of selected diagnostics fields (8.c) in outp.stream n
- (PID.TID 0000.0001) >#                 file for the list of all available diag. in this particular config)
+ (PID.TID 0000.0001) >#                (see "available_diagnostics.log" file for the full list of diags)
- (PID.TID 0000.0001) >#-----------------
+ (PID.TID 0000.0001) >#--------------------
  (PID.TID 0000.0001) > &DIAG_STATIS_PARMS
  (PID.TID 0000.0001) >#- regional mask: 3 lat. band: 1 : y <= -24 ; 2 : -24<y<24 ; 3 : 24 <= y
  (PID.TID 0000.0001) > diagSt_regMaskFile='regMask_lat24.bin',
- (PID.TID 0000.0001) > nSetRegMskFile=1,
+ (PID.TID 0000.0001) > nSetRegMskFile = 1,
- (PID.TID 0000.0001) > set_regMask(1)= 1,  1,  1,
+ (PID.TID 0000.0001) > set_regMask(1:3) = 1,  1,  1,
- (PID.TID 0000.0001) > val_regMask(1)= 1., 2., 3.,
+ (PID.TID 0000.0001) > val_regMask(1:3) = 1., 2., 3.,
  (PID.TID 0000.0001) >#---
- (PID.TID 0000.0001) > stat_fields(1,1)= 'ETAN    ','ETANSQ  ','DETADT2 ',
+ (PID.TID 0000.0001) > stat_fields(1:6,1)  = 'ETAN    ','ETANSQ  ','DETADT2 ',
- (PID.TID 0000.0001) >#                  'UVEL    ','VVEL    ','WVEL    ',
+ (PID.TID 0000.0001) >#                      'UVEL    ','VVEL    ','WVEL    ',
- (PID.TID 0000.0001) >                   'THETA   ','SALT    ','CONVADJ ',
+ (PID.TID 0000.0001) >                       'THETA   ','SALT    ','CONVADJ ',
- (PID.TID 0000.0001) >    stat_fname(1)= 'dynStDiag',
+ (PID.TID 0000.0001) >  stat_fName(1) = 'dynStDiag',
- (PID.TID 0000.0001) >     stat_freq(1)= 864000.,
+ (PID.TID 0000.0001) >   stat_freq(1) = 864000.,
  (PID.TID 0000.0001) >#--- to check conservation of Volume, Heat & Salt (with z*+realFW+AB_2):
- (PID.TID 0000.0001) > stat_fields(1,2)= 'oceTAUX ','oceTAUY ','TRELAX  ','SRELAX  ',
+ (PID.TID 0000.0001) > stat_fields(1:11,2) = 'oceTAUX ','oceTAUY ','TRELAX  ','SRELAX  ',
- (PID.TID 0000.0001) >                   'surForcT','surForcS','TFLUX   ','SFLUX   ',
+ (PID.TID 0000.0001) >                       'surForcT','surForcS','TFLUX   ','SFLUX   ',
- (PID.TID 0000.0001) >                   'oceQnet ','oceSflux','oceFWflx',
+ (PID.TID 0000.0001) >                       'oceQnet ','oceSflux','oceFWflx',
- (PID.TID 0000.0001) >    stat_fname(2)= 'surfStDiag',
+ (PID.TID 0000.0001) >    stat_fName(2) = 'surfStDiag',
- (PID.TID 0000.0001) >     stat_freq(2)= 777600.,
+ (PID.TID 0000.0001) >     stat_freq(2) = 777600.,
- (PID.TID 0000.0001) >    stat_phase(2)= 86400.,
+ (PID.TID 0000.0001) >    stat_phase(2) = 86400.,
- (PID.TID 0000.0001) > stat_fields(1,3)= 'ETAN    ','THETA   ','SALT    ',
+ (PID.TID 0000.0001) > stat_fields(1:3,3)  = 'ETAN    ','THETA   ','SALT    ',
- (PID.TID 0000.0001) >    stat_fname(3)= 'instStDiag',
+ (PID.TID 0000.0001) >    stat_fName(3) = 'instStDiag',
- (PID.TID 0000.0001) >     stat_freq(3)= -777600.,
+ (PID.TID 0000.0001) >     stat_freq(3) = -777600.,
- (PID.TID 0000.0001) >    stat_phase(3)= 86400.,
+ (PID.TID 0000.0001) >    stat_phase(3) = 86400.,
  (PID.TID 0000.0001) >#---
- (PID.TID 0000.0001) > stat_fields(1,5)= 'SI_Fract','SI_Thick','SI_SnowH',
+ (PID.TID 0000.0001) > stat_fields(1:16,5) = 'SI_Fract','SI_Thick','SI_SnowH',
- (PID.TID 0000.0001) >                   'SI_Tsrf ','SI_Tice1','SI_Tice2',
+ (PID.TID 0000.0001) >                       'SI_Tsrf ','SI_Tice1','SI_Tice2',
- (PID.TID 0000.0001) >                   'SI_Qice1','SI_Qice2',
+ (PID.TID 0000.0001) >                       'SI_Qice1','SI_Qice2',
- (PID.TID 0000.0001) >                   'SIsnwPrc','SIalbedo','SIsnwAge',
+ (PID.TID 0000.0001) >                       'SIsnwPrc','SIalbedo','SIsnwAge',
- (PID.TID 0000.0001) >                   'SIflx2oc','SIfrw2oc','SIsaltFx',
+ (PID.TID 0000.0001) >                       'SIflx2oc','SIfrw2oc','SIsaltFx',
- (PID.TID 0000.0001) >                   'SIflxAtm','SIfrwAtm',
+ (PID.TID 0000.0001) >                       'SIflxAtm','SIfrwAtm',
- (PID.TID 0000.0001) >#                  'SItOcMxL','SIsOcMxL',
+ (PID.TID 0000.0001) >#                      'SItOcMxL','SIsOcMxL',
- (PID.TID 0000.0001) > stat_region(1,5)= 1, 3, 0,
+ (PID.TID 0000.0001) > stat_region(1:3,5)  = 1, 3, 0,
- (PID.TID 0000.0001) >    stat_fname(5)= 'thSIceStDiag',
+ (PID.TID 0000.0001) >  stat_fName(5) = 'thSIceStDiag',
- (PID.TID 0000.0001) >     stat_freq(5)= 864000.,
+ (PID.TID 0000.0001) >   stat_freq(5) = 864000.,
  (PID.TID 0000.0001) > /
  (PID.TID 0000.0001) >
  (PID.TID 0000.0001)
-Line 522
+Line 540
  (PID.TID 0000.0001) S/R DIAGNOSTICS_READPARMS, read namelist "diagnostics_list": OK
  (PID.TID 0000.0001) S/R DIAGNOSTICS_READPARMS, read namelist "DIAG_STATIS_PARMS": start
  (PID.TID 0000.0001) S/R DIAGNOSTICS_READPARMS, read namelist "DIAG_STATIS_PARMS": OK
+ (PID.TID 0000.0001)  DIAGNOSTICS_READPARMS: global parameter summary:
+ (PID.TID 0000.0001)  dumpAtLast = /* always write time-ave diags at the end */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001)  diag_mnc =   /* write NetCDF output files */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001)  useMissingValue = /* put MissingValue where mask = 0 */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001)  diagCG_maxIters = /* max number of iters in diag_cg2d */
+ (PID.TID 0000.0001)                     200
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001)  diagCG_resTarget = /* residual target for diag_cg2d */
+ (PID.TID 0000.0001)                 1.000000000000000E-07
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) -----------------------------------------------------
  (PID.TID 0000.0001)  DIAGNOSTICS_READPARMS: active diagnostics summary:
  (PID.TID 0000.0001) -----------------------------------------------------
-Line 529
+Line 563
  (PID.TID 0000.0001) Output Frequency: 1555200000.000000 ; Phase:           0.000000
  (PID.TID 0000.0001)  Averaging Freq.: 1555200000.000000 , Phase:           0.000000 , Cycle:   1
  (PID.TID 0000.0001)  missing value:  1.234567000000E+05 ; for integers:   123456789
- (PID.TID 0000.0001)  Levels:       1.
+ (PID.TID 0000.0001)  Levels:    will be set later
  (PID.TID 0000.0001)  Fields:    ETAN     ETANSQ   DETADT2  PHIBOT   PHIBOTSQ oceTAUX  oceTAUY  TRELAX   SRELAX   surForcT
  (PID.TID 0000.0001)  Fields:    surForcS TFLUX    SFLUX    oceQsw   oceQnet  oceSflux oceFWflx
  (PID.TID 0000.0001) Creating Output Stream: dynDiag
-Line 686
+Line 720
  (PID.TID 0000.0001) %MON AngleSN_sd                   =   6.2496317138039E-01
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: bathy_Hmin50.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.200000000000000E+03
  (PID.TID 0000.0001) // CMAX =         -5.000000000000000E+01
  (PID.TID 0000.0001) // CINT =          1.907407407407407E+02
-Line 728
+Line 762
  (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 =         -7.105427357601002E-15
  (PID.TID 0000.0001) // CMAX =         -7.105427357601002E-15
  (PID.TID 0000.0001) // CINT =          0.000000000000000E+00
-Line 743
+Line 777
  (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
-Line 758
+Line 792
  (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
-Line 773
+Line 807
  (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
-Line 787
+Line 821
  (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 830
+Line 865
  (PID.TID 0000.0001) // ===================================
  (PID.TID 0000.0001) ------------------------------------------------------------
  (PID.TID 0000.0001) DIAGNOSTICS_SET_LEVELS: done
- (PID.TID 0000.0001)  Total Nb of available Diagnostics: ndiagt=   227
+ (PID.TID 0000.0001)  Total Nb of available Diagnostics: ndiagt=   240
  (PID.TID 0000.0001)  write list of available Diagnostics to file: available_diagnostics.log
  (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    23 ETAN
  (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    24 ETANSQ
  (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    25 DETADT2
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    69 PHIBOT
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    71 PHIBOT
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    70 PHIBOTSQ
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    72 PHIBOTSQ
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    75 oceTAUX
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    78 oceTAUX
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    76 oceTAUY
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    79 oceTAUY
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    84 TRELAX
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    87 TRELAX
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    85 SRELAX
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    88 SRELAX
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    86 surForcT
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    89 surForcT
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    87 surForcS
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    90 surForcS
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    88 TFLUX
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    91 TFLUX
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    89 SFLUX
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    92 SFLUX
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    82 oceQsw
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    85 oceQsw
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    81 oceQnet
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    84 oceQnet
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    80 oceSflux
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    83 oceSflux
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    79 oceFWflx
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    82 oceFWflx
  (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    30 UVEL
  (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    31 VVEL
  (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    32 WVEL
- (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    67 PHIHYD
+ (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    69 PHIHYD
  (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    46 VVELMASS
  (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    45 UVELMASS
  (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    38 WVELSQ
  (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    26 THETA
- (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    48 UTHMASS
+ (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    50 UTHMASS
- (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    49 VTHMASS
+ (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    51 VTHMASS
- (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    50 WTHMASS
+ (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    52 WTHMASS
  (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    27 SALT
- (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    51 USLTMASS
+ (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    53 USLTMASS
- (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    52 VSLTMASS
+ (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    54 VSLTMASS
- (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    53 WSLTMASS
+ (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    55 WSLTMASS
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   194 SI_Fract
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   207 SI_Fract
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   195 SI_Thick
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   208 SI_Thick
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   194 SI_Fract has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter-mate #   207 SI_Fract is already set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   196 SI_SnowH
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   209 SI_SnowH
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   194 SI_Fract has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter-mate #   207 SI_Fract is already set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   197 SI_Tsrf
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   210 SI_Tsrf
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   194 SI_Fract has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter-mate #   207 SI_Fract is already set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   198 SI_Tice1
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   211 SI_Tice1
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   194 SI_Fract has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter-mate #   207 SI_Fract is already set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   199 SI_Tice2
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   212 SI_Tice2
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   194 SI_Fract has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter-mate #   207 SI_Fract is already set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   200 SI_Qice1
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   213 SI_Qice1
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   195 SI_Thick has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter-mate #   208 SI_Thick is already set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   201 SI_Qice2
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   214 SI_Qice2
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   195 SI_Thick has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter-mate #   208 SI_Thick is already set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   202 SIsnwAge
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   215 SIsnwAge
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   205 SIsnwPrc
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   218 SIsnwPrc
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Count.Diag #   203 SI_FrcFx
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Mate Diag. #   216 SI_FrcFx
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   204 SIalbedo
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   217 SIalbedo
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   203 SI_FrcFx has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter-mate #   216 SI_FrcFx is already set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   208 SIflx2oc
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   221 SIflx2oc
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   209 SIfrw2oc
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   222 SIfrw2oc
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   210 SIsaltFx
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   223 SIsaltFx
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   206 SIflxAtm
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   219 SIflxAtm
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   207 SIfrwAtm
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x 12 Levels for Diagnostic #   220 SIfrwAtm
  (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    23 ETAN
- (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    79 oceFWflx
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    82 oceFWflx
  (PID.TID 0000.0001) SETDIAG: Allocate  1 x  1 Levels for Diagnostic #    23 ETAN
  (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    26 THETA
  (PID.TID 0000.0001) SETDIAG: Allocate 15 x  1 Levels for Diagnostic #    27 SALT
  (PID.TID 0000.0001)   space allocated for all diagnostics:     479 levels
- (PID.TID 0000.0001)   set mate pointer for diag #    75  oceTAUX  , Parms: UU      U1 , mate:    76
+ (PID.TID 0000.0001)   set mate pointer for diag #    78  oceTAUX  , Parms: UU      U1 , mate:    79
- (PID.TID 0000.0001)   set mate pointer for diag #    76  oceTAUY  , Parms: VV      U1 , mate:    75
+ (PID.TID 0000.0001)   set mate pointer for diag #    79  oceTAUY  , Parms: VV      U1 , mate:    78
  (PID.TID 0000.0001)   set mate pointer for diag #    30  UVEL     , Parms: UUR     MR , mate:    31
  (PID.TID 0000.0001)   set mate pointer for diag #    31  VVEL     , Parms: VVR     MR , mate:    30
  (PID.TID 0000.0001)   set mate pointer for diag #    46  VVELMASS , Parms: VVr     MR , mate:    45
  (PID.TID 0000.0001)   set mate pointer for diag #    45  UVELMASS , Parms: UUr     MR , mate:    46
- (PID.TID 0000.0001)   set mate pointer for diag #    48  UTHMASS  , Parms: UUr     MR , mate:    49
+ (PID.TID 0000.0001)   set mate pointer for diag #    50  UTHMASS  , Parms: UUr     MR , mate:    51
- (PID.TID 0000.0001)   set mate pointer for diag #    49  VTHMASS  , Parms: VVr     MR , mate:    48
+ (PID.TID 0000.0001)   set mate pointer for diag #    51  VTHMASS  , Parms: VVr     MR , mate:    50
- (PID.TID 0000.0001)   set mate pointer for diag #    51  USLTMASS , Parms: UUr     MR , mate:    52
+ (PID.TID 0000.0001)   set mate pointer for diag #    53  USLTMASS , Parms: UUr     MR , mate:    54
- (PID.TID 0000.0001)   set mate pointer for diag #    52  VSLTMASS , Parms: VVr     MR , mate:    51
+ (PID.TID 0000.0001)   set mate pointer for diag #    54  VSLTMASS , Parms: VVr     MR , mate:    53
+ (PID.TID 0000.0001) DIAGNOSTICS_SET_POINTERS: Set levels for Outp.Stream: surfDiag
+ (PID.TID 0000.0001)  Levels:       1.
  (PID.TID 0000.0001) DIAGNOSTICS_SET_POINTERS: Set levels for Outp.Stream: dynDiag
  (PID.TID 0000.0001)  Levels:       1.   2.   3.   4.   5.   6.   7.   8.   9.  10.  11.  12.  13.  14.  15.
  (PID.TID 0000.0001) DIAGNOSTICS_SET_POINTERS: Set levels for Outp.Stream: thSIceDiag
-Line 928
+Line 965
  (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    25 DETADT2
  (PID.TID 0000.0001) SETDIAG: Allocate 15 Levels for Stats-Diag #    26 THETA
  (PID.TID 0000.0001) SETDIAG: Allocate 15 Levels for Stats-Diag #    27 SALT
- (PID.TID 0000.0001) SETDIAG: Allocate 15 Levels for Stats-Diag #    74 CONVADJ
+ (PID.TID 0000.0001) SETDIAG: Allocate 15 Levels for Stats-Diag #    77 CONVADJ
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    75 oceTAUX
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    78 oceTAUX
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    76 oceTAUY
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    79 oceTAUY
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    84 TRELAX
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    87 TRELAX
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    85 SRELAX
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    88 SRELAX
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    86 surForcT
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    89 surForcT
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    87 surForcS
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    90 surForcS
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    88 TFLUX
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    91 TFLUX
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    89 SFLUX
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    92 SFLUX
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    81 oceQnet
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    84 oceQnet
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    80 oceSflux
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    83 oceSflux
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    79 oceFWflx
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    82 oceFWflx
  (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #    23 ETAN
  (PID.TID 0000.0001) SETDIAG: Allocate 15 Levels for Stats-Diag #    26 THETA
  (PID.TID 0000.0001) SETDIAG: Allocate 15 Levels for Stats-Diag #    27 SALT
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   194 SI_Fract
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   207 SI_Fract
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   195 SI_Thick
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   208 SI_Thick
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   194 SI_Fract has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   207 SI_Fract has already been set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   196 SI_SnowH
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   209 SI_SnowH
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   194 SI_Fract has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   207 SI_Fract has already been set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   197 SI_Tsrf
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   210 SI_Tsrf
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   194 SI_Fract has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   207 SI_Fract has already been set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   198 SI_Tice1
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   211 SI_Tice1
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   194 SI_Fract has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   207 SI_Fract has already been set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   199 SI_Tice2
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   212 SI_Tice2
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   194 SI_Fract has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   207 SI_Fract has already been set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   200 SI_Qice1
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   213 SI_Qice1
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   195 SI_Thick has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   208 SI_Thick has already been set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   201 SI_Qice2
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   214 SI_Qice2
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   195 SI_Thick has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   208 SI_Thick has already been set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   205 SIsnwPrc
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   218 SIsnwPrc
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Count.Diag #   203 SI_FrcFx
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Count.Diag #   216 SI_FrcFx
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   204 SIalbedo
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   217 SIalbedo
- (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   203 SI_FrcFx has already been set
+ (PID.TID 0000.0001) - NOTE - SETDIAG: Counter Diagnostic #   216 SI_FrcFx has already been set
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   202 SIsnwAge
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   215 SIsnwAge
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   208 SIflx2oc
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   221 SIflx2oc
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   209 SIfrw2oc
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   222 SIfrw2oc
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   210 SIsaltFx
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   223 SIsaltFx
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   206 SIflxAtm
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   219 SIflxAtm
- (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   207 SIfrwAtm
+ (PID.TID 0000.0001) SETDIAG: Allocate  1 Levels for Stats-Diag #   220 SIfrwAtm
  (PID.TID 0000.0001)   space allocated for all stats-diags:     107 levels
  (PID.TID 0000.0001) DIAGSTATS_SET_POINTERS: done
  (PID.TID 0000.0001) ------------------------------------------------------------
-Line 1164
+Line 1201
  (PID.TID 0000.0001) beta =   /* Beta ( 1/(m.s) ) */
  (PID.TID 0000.0001)                 9.999999999999999E-12
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) fPrime =   /* Second coriolis parameter ( 1/s ) */
+ (PID.TID 0000.0001)                 0.000000000000000E+00
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) rigidLid =   /* Rigid lid on/off flag */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
-Line 1194
+Line 1234
  (PID.TID 0000.0001) hFacMinDr = /* minimum partial cell thickness ( m) */
  (PID.TID 0000.0001)                 1.000000000000000E-01
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) nonlinFreeSurf =  /* Non-linear Free Surf. options (-1,0,1,2,3)*/
+ (PID.TID 0000.0001) nonlinFreeSurf = /* Non-linear Free Surf. options (-1,0,1,2,3)*/
  (PID.TID 0000.0001)                       4
- (PID.TID 0000.0001)     ;
  (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)     ;
-Line 1267
+Line 1307
  (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)                       2
- (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)                   F
- (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) useSphereF  =   /* use Spherical Coriolis flag */
- (PID.TID 0000.0001)                   T
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) use3dCoriolis = /* 3-D Coriolis on/off flag */
  (PID.TID 0000.0001)                   F
-Line 1285
+Line 1320
  (PID.TID 0000.0001) useCDscheme =  /* CD scheme on/off flag */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) useEnergyConservingCoriolis= /* Flx-Form Coriolis scheme flag */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) useJamartWetPoints= /* Coriolis WetPoints method flag */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) useJamartMomAdv= /* V.I. Non-linear terms Jamart flag */
+ (PID.TID 0000.0001) useJamartMomAdv= /* V.I Non-linear terms Jamart flag */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) useAbsVorticity= /* Work with f+zeta in Coriolis */
+ (PID.TID 0000.0001) useAbsVorticity= /* V.I Works with f+zeta in Coriolis */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) selectVortScheme= /* Scheme selector for Vorticity-Term */
+ (PID.TID 0000.0001) selectVortScheme= /* V.I Scheme selector for Vorticity-Term */
  (PID.TID 0000.0001)                       1
  (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
-Line 1302
+Line 1340
  (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= /* Upwind bias vorticity flag */
+ (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 1377
+Line 1415
  (PID.TID 0000.0001)  useSingleCpuIO = /* only master MPI process does I/O */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001)  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 1410
+Line 1446
  (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) //
-Line 1446
+Line 1485
  (PID.TID 0000.0001) pickupStrictlyMatch= /* stop if pickup do not strictly match */
  (PID.TID 0000.0001)                   T
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) nIter0 =   /* Run starting timestep number  */
+ (PID.TID 0000.0001) nIter0   =   /* Run starting timestep number */
  (PID.TID 0000.0001)                   36000
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) nTimeSteps =   /* Number of timesteps */
+ (PID.TID 0000.0001) nTimeSteps = /* Number of timesteps */
  (PID.TID 0000.0001)                      20
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) baseTime =   /* Model base time ( s ). */
+ (PID.TID 0000.0001) nEndIter =   /* Run ending timestep number */
+ (PID.TID 0000.0001)                   36020
+ (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)                 3.110400000000000E+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)                 3.112128000000000E+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)                 3.110400000000000E+07
  (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 1603
+Line 1645
  (PID.TID 0000.0001) rSphere =  /* Radius ( ignored - cartesian, m - spherical ) */
  (PID.TID 0000.0001)                 6.370000000000000E+06
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) radius_fromHorizGrid = /* sphere Radius of input horiz. grid */
+ (PID.TID 0000.0001)                 6.370000000000000E+06
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) deepAtmosphere = /* Deep/Shallow Atmosphere flag (True/False) */
  (PID.TID 0000.0001)                   F
  (PID.TID 0000.0001)     ;
-Line 1701
+Line 1746
  (PID.TID 0000.0001) wUnit2rVel = /* convert units: wSpeed -> rVel (=1 if z-coord)*/
  (PID.TID 0000.0001)    16 @  1.000000000000000E+00              /* K =  1: 16 */
  (PID.TID 0000.0001)     ;
- (PID.TID 0000.0001) dBdrRef = /* Vertical gradient of reference boyancy [(m/s/r)^2)] */
+ (PID.TID 0000.0001) dBdrRef = /* Vertical grad. of reference buoyancy [(m/s/r)^2] */
  (PID.TID 0000.0001)    15 @  0.000000000000000E+00              /* K =  1: 15 */
  (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) rotateGrid = /* use rotated grid ( True/False ) */
-Line 2253
+Line 2298
  (PID.TID 0000.0001)
  (PID.TID 0000.0001) == Packages configuration : Check & print summary ==
  (PID.TID 0000.0001)
- (PID.TID 0000.0001) GAD_CHECK: #define ALLOW_GENERIC_ADVDIFF
  (PID.TID 0000.0001) GMREDI_CHECK: #define GMREDI
  (PID.TID 0000.0001) GM_AdvForm =     /* if FALSE => use SkewFlux Form */
  (PID.TID 0000.0001)                   F
-Line 2303
+Line 2347
  (PID.TID 0000.0001) GM_maxTransLay = /* Maximum Transition Layer Thickness [m] */
  (PID.TID 0000.0001)                 5.000000000000000E+02
  (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) GM_UseBVP = /* if TRUE => use bvp a la Ferrari et al. (2010) */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) GM_BVP_ModeNumber = /* Vertical mode number for BVP wave speed */
+ (PID.TID 0000.0001)                       1
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) GM_BVP_cMin = /* Minimum wave speed for BVP [m/s] */
+ (PID.TID 0000.0001)                 1.000000000000000E-01
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) GM_useSubMeso = /* if TRUE => use Sub-Meso param. (B.Fox-Kemper) */
+ (PID.TID 0000.0001)                   F
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) subMeso_Ceff = /* efficiency coeff. of Mixed-Layer Eddies [-] */
+ (PID.TID 0000.0001)                 7.000000000000001E-02
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) subMeso_invTau = /* inverse of Sub-Meso mixing time-scale [/s] */
+ (PID.TID 0000.0001)                 2.000000000000000E-06
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) subMeso_LfMin = /* minimum length-scale "Lf" [m] */
+ (PID.TID 0000.0001)                 1.000000000000000E+03
+ (PID.TID 0000.0001)     ;
+ (PID.TID 0000.0001) subMeso_Lmax = /* maximum grid-scale length [m] */
+ (PID.TID 0000.0001)                 1.100000000000000E+05
+ (PID.TID 0000.0001)     ;
  (PID.TID 0000.0001) THSICE_CHECK: #define THSICE
+ (PID.TID 0000.0001) GAD_CHECK: #define ALLOW_GENERIC_ADVDIFF
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // CONFIG_CHECK : Normal End
  (PID.TID 0000.0001) // =======================================================
-Line 2356 
 listId=    1 ; file name: surfDiag
+Line 2425 
 listId=    1 ; file name: surfDiag
 |ETAN    |      1 |      0 |   1 |       0 |
 |ETANSQ  |      2 |      0 |   1 |       0 |
 |DETADT2 |      3 |      0 |   1 |       0 |
 |PHIBOT  |      4 |      0 |   1 |       0 |
 |PHIBOTSQ|      5 |      0 |   1 |       0 |
 |oceTAUX |      6 |      7 |   1 |       0 |       0 |
 |oceTAUY |      7 |      6 |   1 |       0 |       0 |
 |TRELAX  |      8 |      0 |   1 |       0 |
 |SRELAX  |      9 |      0 |   1 |       0 |
 |surForcT|     10 |      0 |   1 |       0 |
 |surForcS|     11 |      0 |   1 |       0 |
 |TFLUX   |     12 |      0 |   1 |       0 |
 |SFLUX   |     13 |      0 |   1 |       0 |
 |oceQsw  |     14 |      0 |   1 |       0 |
 |oceQnet |     15 |      0 |   1 |       0 |
 |oceSflux|     16 |      0 |   1 |       0 |
 |oceFWflx|     17 |      0 |   1 |       0 |
  ------------------------------------------------------------------------
  listId=    2 ; file name: dynDiag
   nFlds, nActive,       freq     &     phase        , nLev
-Line 2379 
 listId=    2 ; file name: dynDiag
+Line 2448 
 listId=    2 ; file name: dynDiag
 |UVEL    |     18 |     33 |  15 |       0 |       0 |
 |VVEL    |     33 |     18 |  15 |       0 |       0 |
 |WVEL    |     48 |      0 |  15 |       0 |
 |PHIHYD  |     63 |      0 |  15 |       0 |
 |VVELMASS|     78 |     93 |  15 |       0 |       0 |
 |UVELMASS|     93 |     78 |  15 |       0 |       0 |
 |WVELSQ  |    108 |      0 |  15 |       0 |
 |THETA   |    123 |      0 |  15 |       0 |
 |UTHMASS |    138 |    153 |  15 |       0 |       0 |
 |VTHMASS |    153 |    138 |  15 |       0 |       0 |
 |WTHMASS |    168 |      0 |  15 |       0 |
 |SALT    |    183 |      0 |  15 |       0 |
 |USLTMASS|    198 |    213 |  15 |       0 |       0 |
 |VSLTMASS|    213 |    198 |  15 |       0 |       0 |
 |WSLTMASS|    228 |      0 |  15 |       0 |
  ------------------------------------------------------------------------
  listId=    3 ; file name: thSIceDiag
   nFlds, nActive,       freq     &     phase        , nLev
  |   17  |1555200000.000000         0.000000 |   1
   levels:   1
   diag# | name   |   ipt  |  iMate | kLev|   count |   mate.C|
 |SI_Fract|    243 |      0 |   1 |      0(x 12) |
 |SI_Thick|    255 |    243 |   1 |      0(x 12) |
 |SI_SnowH|    267 |    243 |   1 |      0(x 12) |
 |SI_Tsrf |    279 |    243 |   1 |      0(x 12) |
 |SI_Tice1|    291 |    243 |   1 |      0(x 12) |
 |SI_Tice2|    303 |    243 |   1 |      0(x 12) |
 |SI_Qice1|    315 |    255 |   1 |      0(x 12) |
 |SI_Qice2|    327 |    255 |   1 |      0(x 12) |
 |SIsnwAge|    339 |      0 |   1 |      0(x 12) |
 |SIsnwPrc|    351 |    363 |   1 |      0(x 12) |
 |SIalbedo|    375 |    363 |   1 |      0(x 12) |
 |SIflx2oc|    387 |      0 |   1 |      0(x 12) |
 |SIfrw2oc|    399 |      0 |   1 |      0(x 12) |
 |SIsaltFx|    411 |      0 |   1 |      0(x 12) |
 |SIflxAtm|    423 |      0 |   1 |      0(x 12) |
 |SIfrwAtm|    435 |      0 |   1 |      0(x 12) |
 |SI_FrcFx|    363 |      0 |   1 |      0(x 12) |
  ------------------------------------------------------------------------
  listId=    4 ; file name: surfInst
   nFlds, nActive,       freq     &     phase        , nLev
-Line 2421 
 listId=    4 ; file name: surfInst
+Line 2490 
 listId=    4 ; file name: surfInst
   levels:   1
   diag# | name   |   ipt  |  iMate | kLev|   count |   mate.C|
 |ETAN    |    447 |      0 |   1 |       0 |
 |oceFWflx|    448 |      0 |   1 |       0 |
  ------------------------------------------------------------------------
  listId=    5 ; file name: etaInst
   nFlds, nActive,       freq     &     phase        , nLev
-Line 2450 
 listId=   1 ; file name: dynStDiag
+Line 2519 
 listId=   1 ; file name: dynStDiag
 |DETADT2 |      3 |      0 | 0.00000E+00 |
 |THETA   |      4 |      0 | 0.00000E+00 |
 |SALT    |     19 |      0 | 0.00000E+00 |
 |CONVADJ |     34 |      0 | 0.00000E+00 |
  ------------------------------------------------------------------------
  listId=   2 ; file name: surfStDiag
   nFlds, nActive,       freq     &     phase        |
  |   11  |    777600.000000     86400.000000 |
   Regions:   0
   diag# | name   |   ipt  |  iMate |    Volume   |   mate-Vol. |
 |oceTAUX |     49 |      0 | 0.00000E+00 |
 |oceTAUY |     50 |      0 | 0.00000E+00 |
 |TRELAX  |     51 |      0 | 0.00000E+00 |
 |SRELAX  |     52 |      0 | 0.00000E+00 |
 |surForcT|     53 |      0 | 0.00000E+00 |
 |surForcS|     54 |      0 | 0.00000E+00 |
 |TFLUX   |     55 |      0 | 0.00000E+00 |
 |SFLUX   |     56 |      0 | 0.00000E+00 |
 |oceQnet |     57 |      0 | 0.00000E+00 |
 |oceSflux|     58 |      0 | 0.00000E+00 |
 |oceFWflx|     59 |      0 | 0.00000E+00 |
  ------------------------------------------------------------------------
  listId=   3 ; file name: instStDiag
   nFlds, nActive,       freq     &     phase        |
-Line 2483 
 listId=   4 ; file name: thSIceStDiag
+Line 2552 
 listId=   4 ; file name: thSIceStDiag
  |   17  |    864000.000000         0.000000 |
   Regions:   0  1  3
   diag# | name   |   ipt  |  iMate |    Volume   |   mate-Vol. |
 |SI_Fract|     91 |      0 | 0.00000E+00 |
 |SI_Thick|     92 |     91 | 0.00000E+00 | 0.00000E+00 |
 |SI_SnowH|     93 |     91 | 0.00000E+00 | 0.00000E+00 |
 |SI_Tsrf |     94 |     91 | 0.00000E+00 | 0.00000E+00 |
 |SI_Tice1|     95 |     91 | 0.00000E+00 | 0.00000E+00 |
 |SI_Tice2|     96 |     91 | 0.00000E+00 | 0.00000E+00 |
 |SI_Qice1|     97 |     92 | 0.00000E+00 | 0.00000E+00 |
 |SI_Qice2|     98 |     92 | 0.00000E+00 | 0.00000E+00 |
 |SIsnwPrc|     99 |    100 | 0.00000E+00 | 0.00000E+00 |
 |SIalbedo|    101 |    100 | 0.00000E+00 | 0.00000E+00 |
 |SIsnwAge|    102 |      0 | 0.00000E+00 |
 |SIflx2oc|    103 |      0 | 0.00000E+00 |
 |SIfrw2oc|    104 |      0 | 0.00000E+00 |
 |SIsaltFx|    105 |      0 | 0.00000E+00 |
 |SIflxAtm|    106 |      0 | 0.00000E+00 |
 |SIfrwAtm|    107 |      0 | 0.00000E+00 |
 |SI_FrcFx|    100 |      0 | 0.00000E+00 |
  ------------------------------------------------------------------------
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: pickup_ic.0000036000
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: pickup_ic.0000036000
-Line 2635 
 listId=   4 ; file name: thSIceStDiag
+Line 2704 
 listId=   4 ; file name: thSIceStDiag
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
- S/R BULKF_FIELDS_LOAD: Reading new data:   12    1     36000  3.110400000000E+09
+  BULKF_FIELDS_LOAD,     36000 : iP,iLd,i0,i1=   12    0   12    1 ; Wght=  0.5000000000  0.5000000000
+  BULKF_FIELDS_LOAD, it=     36000 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_tair_cs.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_tair_cs.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_qair_cs.bin
-Line 2648 
 S/R BULKF_FIELDS_LOAD: Reading new data:
+Line 2718 
 S/R BULKF_FIELDS_LOAD: Reading new data:
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_windspeed_cs.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_pr_scal_cs.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_pr_scal_cs.bin
- S/R EXTERNAL_FIELDS_LOAD: Reading new data:   12    1     36000  3.110400000000E+09
+  EXTERNAL_FIELDS_LOAD,     36000 : iP,iLd,i0,i1=   12    0   12    1 ; Wght=  0.5000000000  0.5000000000
+  EXTERNAL_FIELDS_LOAD, it=     36000 : Reading new data, i0,i1=   12    1 (prev=   12    0 )
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: trenberth_taux.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: trenberth_taux.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: trenberth_tauy.bin
-Line 2663 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 2734 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
   cg2d: Sum(rhs),rhsMax =   3.91123171841600E+02  2.29216909151821E+02
  (PID.TID 0000.0001)                    cg2d_init_res =    1.95242193209935E+01
  (PID.TID 0000.0001)                       cg2d_iters =    72
- (PID.TID 0000.0001)                         cg2d_res =    5.46408886398732E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.46408886396900E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36001
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1104864000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.3388365064227E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.3388365064228E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1986231184597E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8173432388293E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.7027393996681E+00
-Line 2686 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 2757 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.6363744991273E-05
  (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.3985212379618E-05
  (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7706298341014E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   1.7416098782924E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =   1.7416098782922E-09
  (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6302076880399E-06
  (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0408133923544E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   3.0143977320706E+01
-Line 2738 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 2809 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259583286529E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787190486060E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817278636374E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   3.7878597719727E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   3.7878597719771E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   9.2604174719932E-10
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   9.2604174790098E-10
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 2786 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 2857 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
   Computing Diagnostic #     23  ETAN         Counter:       1   Parms: SM      M1
-  Computing Diagnostic #     79  oceFWflx     Counter:       1   Parms: SM      U1
+  Computing Diagnostic #     82  oceFWflx     Counter:       1   Parms: SM      U1
   Computing Diagnostic #     26  THETA        Counter:       1   Parms: SMR     MR
   Computing Diagnostic #     27  SALT         Counter:       1   Parms: SMR     MR
-  Compute Stats, Diag. #     75  oceTAUX   vol(   0 ): 3.457E+14  Parms: UU      U1
+  Compute Stats, Diag. #     78  oceTAUX   vol(   0 ): 3.457E+14  Parms: UU      U1
-  Compute Stats, Diag. #     76  oceTAUY   vol(   0 ): 3.478E+14  Parms: VV      U1
+  Compute Stats, Diag. #     79  oceTAUY   vol(   0 ): 3.478E+14  Parms: VV      U1
-  Compute Stats, Diag. #     84  TRELAX    vol(   0 ): 3.639E+14  Parms: SM      U1
+  Compute Stats, Diag. #     87  TRELAX    vol(   0 ): 3.639E+14  Parms: SM      U1
-  Compute Stats, Diag. #     85  SRELAX    vol(   0 ): 3.639E+14  Parms: SM      U1
+  Compute Stats, Diag. #     88  SRELAX    vol(   0 ): 3.639E+14  Parms: SM      U1
-  Compute Stats, Diag. #     86  surForcT  vol(   0 ): 3.639E+14  Parms: SM      U1
+  Compute Stats, Diag. #     89  surForcT  vol(   0 ): 3.639E+14  Parms: SM      U1
-  Compute Stats, Diag. #     87  surForcS  vol(   0 ): 3.639E+14  Parms: SM      U1
+  Compute Stats, Diag. #     90  surForcS  vol(   0 ): 3.639E+14  Parms: SM      U1
-  Compute Stats, Diag. #     88  TFLUX     vol(   0 ): 3.639E+14  Parms: SM      U1
+  Compute Stats, Diag. #     91  TFLUX     vol(   0 ): 3.639E+14  Parms: SM      U1
-  Compute Stats, Diag. #     89  SFLUX     vol(   0 ): 3.639E+14  Parms: SM      U1
+  Compute Stats, Diag. #     92  SFLUX     vol(   0 ): 3.639E+14  Parms: SM      U1
-  Compute Stats, Diag. #     81  oceQnet   vol(   0 ): 3.639E+14  Parms: SM      U1
+  Compute Stats, Diag. #     84  oceQnet   vol(   0 ): 3.639E+14  Parms: SM      U1
-  Compute Stats, Diag. #     80  oceSflux  vol(   0 ): 3.639E+14  Parms: SM      U1
+  Compute Stats, Diag. #     83  oceSflux  vol(   0 ): 3.639E+14  Parms: SM      U1
-  Compute Stats, Diag. #     79  oceFWflx  vol(   0 ): 3.639E+14  Parms: SM      U1
+  Compute Stats, Diag. #     82  oceFWflx  vol(   0 ): 3.639E+14  Parms: SM      U1
-  cg2d: Sum(rhs),rhsMax =   3.91131921824871E+02  2.29255287241867E+02
+  BULKF_FIELDS_LOAD,     36001 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.4666666667  0.5333333333
- (PID.TID 0000.0001)                    cg2d_init_res =    2.74856936053468E+00
+  EXTERNAL_FIELDS_LOAD,     36001 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.4666666667  0.5333333333
+  cg2d: Sum(rhs),rhsMax =   3.91131921824870E+02  2.29255287241867E+02
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.74856936053251E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.38174366779394E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.38174366778802E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36002
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1105728000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.3253161004735E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.3253161004741E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1985096656221E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8174733810956E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.7022747292739E+00
-Line 2826 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 2899 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.6319623354942E-05
  (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.3878461487957E-05
  (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7741686258955E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2005582791155E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2005582791150E-09
  (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6256908530211E-06
  (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0321202865293E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   3.0124603885919E+01
-Line 2878 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 2951 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259583113972E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787196681009E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817293504035E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   8.5120869160650E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   8.5120869160542E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.8910280778552E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.8910280776704E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 2929 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3002 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
   Compute Stats, Diag. #     23  ETAN      vol(   0 ): 3.639E+14  Parms: SM      M1
   Compute Stats, Diag. #     26  THETA     vol(   0 ): 1.340E+18  Parms: SMR     MR
   Compute Stats, Diag. #     27  SALT      vol(   0 ): 1.340E+18  Parms: SMR     MR
+  BULKF_FIELDS_LOAD,     36002 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.4333333333  0.5666666667
+  EXTERNAL_FIELDS_LOAD,     36002 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.4333333333  0.5666666667
   cg2d: Sum(rhs),rhsMax =   3.91306457345962E+02  2.29299174390983E+02
- (PID.TID 0000.0001)                    cg2d_init_res =    6.05881510161922E+00
+ (PID.TID 0000.0001)                    cg2d_init_res =    6.05881510161933E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.26814301672221E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.26814301673426E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36003
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1106592000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.3125883441893E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.3125883441881E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1983414707440E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8200693230051E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.7027138668921E+00
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1151932971679E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1151932971680E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7366837431578E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.7047805389613E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2788824696051E-04
-Line 2955 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3030 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.6276318679090E-05
  (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.3856395721696E-05
  (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7771419054994E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -6.4725927528087E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -6.4725927528003E-10
  (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6253268016134E-06
  (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0333590472275E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   3.0104705418602E+01
-Line 3007 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3082 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582970957E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787196249296E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817294256821E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.0416070374144E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.0416070374209E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.4810448161579E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.4810448162180E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 3054 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3129 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
+  BULKF_FIELDS_LOAD,     36003 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.4000000000  0.6000000000
+  EXTERNAL_FIELDS_LOAD,     36003 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.4000000000  0.6000000000
   cg2d: Sum(rhs),rhsMax =   3.91399720621000E+02  2.29353133798250E+02
- (PID.TID 0000.0001)                    cg2d_init_res =    4.24480993413539E+00
+ (PID.TID 0000.0001)                    cg2d_init_res =    4.24480993414271E+00
  (PID.TID 0000.0001)                       cg2d_iters =    68
- (PID.TID 0000.0001)                         cg2d_res =    5.39090247348569E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.39090247345064E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36004
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1107456000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.3006737311701E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.3006737311703E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1981545912127E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8214564678268E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.7027252387751E+00
-Line 3080 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3157 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.6233996433000E-05
  (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.3910409654922E-05
  (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7794292282195E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -8.6206715887120E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -8.6206715887194E-10
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6243834885156E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6243834885155E-06
  (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0332213936820E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   3.0084398517367E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.3484357320152E+00
-Line 3132 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3209 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582848045E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787197131553E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817297555102E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.5113978542746E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.5113978542828E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.6693583634293E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.6693583635308E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 3179 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3256 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
+  BULKF_FIELDS_LOAD,     36004 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.3666666667  0.6333333333
+  EXTERNAL_FIELDS_LOAD,     36004 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.3666666667  0.6333333333
   cg2d: Sum(rhs),rhsMax =   3.91464684334993E+02  2.29415843069631E+02
- (PID.TID 0000.0001)                    cg2d_init_res =    3.32097632146683E+00
+ (PID.TID 0000.0001)                    cg2d_init_res =    3.32097632147858E+00
  (PID.TID 0000.0001)                       cg2d_iters =    67
- (PID.TID 0000.0001)                         cg2d_res =    5.58497136618731E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.58497136615296E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36005
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1108320000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2890681292899E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2890681292887E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1979604066709E+01
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1979604066710E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8224227010996E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.7024609847232E+00
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1226410257466E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1226410257467E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7422409787436E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6976679326141E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2669522207994E-04
-Line 3203 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3282 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.2379919086165E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3480489249256E-02
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.6192949048934E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.4614302458593E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.4614302458597E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7810860499615E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7810860499616E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -9.8426246410208E-10
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -9.8426246410086E-10
  (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6236585157975E-06
  (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0333789348411E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   3.0063765078410E+01
-Line 3245 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3324 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.3568352349763E-04
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.6510979690671E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4948231906690E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8070126308180E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8070126308182E-02
  (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4157333978968E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9317608591373E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.6375236009244E-02
-Line 3257 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3336 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582744004E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787197416595E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817299350069E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.7332195251879E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.7332195251893E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.7649090156800E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.7649090157239E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 3304 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3383 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
-  cg2d: Sum(rhs),rhsMax =   3.91499407618689E+02  2.29487198217044E+02
+  BULKF_FIELDS_LOAD,     36005 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.3333333333  0.6666666667
- (PID.TID 0000.0001)                    cg2d_init_res =    2.92897091813487E+00
+  EXTERNAL_FIELDS_LOAD,     36005 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.3333333333  0.6666666667
+  cg2d: Sum(rhs),rhsMax =   3.91499407618690E+02  2.29487198217044E+02
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.92897091812492E+00
  (PID.TID 0000.0001)                       cg2d_iters =    67
- (PID.TID 0000.0001)                         cg2d_res =    5.52939183479604E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.52939183475855E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36006
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1109184000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2774338669258E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2774338669257E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1977591659910E+01
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1977591659909E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8229391554397E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.7020417777615E+00
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1247541957359E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1247541957358E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7440136310914E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6939362810533E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2611418104018E-04
-Line 3328 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3409 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.2330017164161E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3481279915117E-02
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.6153260274843E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.5348942330834E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.5348942330838E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7821510534710E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7821510534711E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.0976019212687E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.0976019212710E-09
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6230663056738E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6230663056739E-06
  (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0336638819912E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   3.0042876824002E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.3414282692797E+00
-Line 3355 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3436 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_qsw_del2             =   6.7110784977734E-01
  (PID.TID 0000.0001) %MON forcing_empmr_max            =   1.7900213910895E-03
  (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.6484175422527E-04
- (PID.TID 0000.0001) %MON forcing_empmr_mean           =   6.1866926168019E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_mean           =   6.1866926168018E-07
  (PID.TID 0000.0001) %MON forcing_empmr_sd             =   6.2707843278933E-05
  (PID.TID 0000.0001) %MON forcing_empmr_del2           =   1.5687472978895E-06
  (PID.TID 0000.0001) %MON forcing_fu_max               =   2.4736588936290E-01
-Line 3370 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3451 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.3667008926707E-04
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.6405142973951E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.5023097808434E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8104849441471E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8104849441473E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4189491076931E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4189491076933E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9315100138425E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.6258217139798E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6195939414723E-04
-Line 3382 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3463 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582657991E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787197542056E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817300431838E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.8574529478509E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.8574529478393E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.8270703060942E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.8270703059336E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 3429 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3510 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
-  cg2d: Sum(rhs),rhsMax =   3.91509297222492E+02  2.29566698430078E+02
+  BULKF_FIELDS_LOAD,     36006 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.3000000000  0.7000000000
- (PID.TID 0000.0001)                    cg2d_init_res =    2.77965884895296E+00
+  EXTERNAL_FIELDS_LOAD,     36006 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.3000000000  0.7000000000
+  cg2d: Sum(rhs),rhsMax =   3.91509297222492E+02  2.29566698430077E+02
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.77965884895224E+00
  (PID.TID 0000.0001)                       cg2d_iters =    68
- (PID.TID 0000.0001)                         cg2d_res =    5.36088828903311E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.36088828904090E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36007
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1110048000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2655174601428E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2655174601438E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1975499478975E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8230862477727E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.7015447107648E+00
  (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1263285270775E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7451052652463E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6900862111964E-01
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2557746617453E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2557746617452E-04
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   1.2980438384365E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8683108917509E-05
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8683108917508E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.9516646725869E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -2.1116269833282E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.2266492485912E-04
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.2266492485911E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3481722312011E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.6114961343228E-05
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.6114961343229E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.6107233909637E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.6107233909641E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7826799941972E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7826799941973E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.1668819704194E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.1668819704190E-09
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6225102263581E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6225102263582E-06
  (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0337458717576E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   3.0021802048657E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.3488689733355E+00
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.5032138058379E+00
  (PID.TID 0000.0001) %MON dynstat_theta_sd             =   4.7493469843194E+00
- (PID.TID 0000.0001) %MON dynstat_theta_del2           =   2.9518636745365E-03
+ (PID.TID 0000.0001) %MON dynstat_theta_del2           =   2.9518636745366E-03
  (PID.TID 0000.0001) %MON dynstat_salt_max             =   4.0831024834176E+01
  (PID.TID 0000.0001) %MON dynstat_salt_min             =   1.8034317400091E+01
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4756113372676E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   5.1183253054356E-01
- (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.3912265362257E-04
+ (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.3912265362258E-04
  (PID.TID 0000.0001) %MON forcing_qnet_max             =   7.9038456668651E+02
  (PID.TID 0000.0001) %MON forcing_qnet_min             =  -3.6658224186914E+02
  (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -5.8071040373559E+00
  (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.2059615685668E+02
- (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0630450047630E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0630450047629E+00
  (PID.TID 0000.0001) %MON forcing_qsw_max              =  -0.0000000000000E+00
  (PID.TID 0000.0001) %MON forcing_qsw_min              =  -3.4073236530849E+02
  (PID.TID 0000.0001) %MON forcing_qsw_mean             =  -1.8655107526619E+02
-Line 3480 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3563 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_qsw_del2             =   6.6999148226222E-01
  (PID.TID 0000.0001) %MON forcing_empmr_max            =   1.0846852335446E-03
  (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.5175627072566E-04
- (PID.TID 0000.0001) %MON forcing_empmr_mean           =   1.7620435717364E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_mean           =   1.7620435717365E-07
  (PID.TID 0000.0001) %MON forcing_empmr_sd             =   5.5235346001753E-05
  (PID.TID 0000.0001) %MON forcing_empmr_del2           =   1.1727819337975E-06
  (PID.TID 0000.0001) %MON forcing_fu_max               =   2.4731784009269E-01
-Line 3494 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3577 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_sd                =   7.5544817904543E-02
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.3785256950118E-04
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.6295947689734E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.5070313531068E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.5070313531069E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8122094905147E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8122094905149E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4208557283339E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4208557283342E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9308648057314E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.6138673439951E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6194057222906E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.3395905563296E+18
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.2381845512428E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.3862365782034E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.3862365782033E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =  -2.0549865324846E-05
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582585218E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787197428633E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817300951655E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.9593287171036E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.9593287170884E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.8724823264992E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.8724823264087E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 3554 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3637 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
+  BULKF_FIELDS_LOAD,     36007 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.2666666667  0.7333333333
+  EXTERNAL_FIELDS_LOAD,     36007 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.2666666667  0.7333333333
   cg2d: Sum(rhs),rhsMax =   3.91503768623070E+02  2.29653736900058E+02
- (PID.TID 0000.0001)                    cg2d_init_res =    2.70203148484498E+00
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.70203148482061E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.11577335738897E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.11577335743439E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36008
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1110912000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2531873118639E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2531873118654E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1973319578344E+01
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1973319578343E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8230040185345E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.7010267589623E+00
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1276522122115E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1276522122114E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7455221937551E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6861149336031E-01
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2508413735098E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2508413735097E-04
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   1.2978388120682E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8658040502310E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.9521111920796E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -2.1122941163268E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.2193417761524E-04
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.2193417761523E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3481775601670E-02
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.6078030979900E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.6876977669111E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.6876977669117E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7827057529016E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7827057529017E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.1946867500248E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.1946867500240E-09
  (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6220873878842E-06
  (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0338023837694E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   3.0000605634138E+01
-Line 3592 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3677 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_salt_min             =   1.8034875180377E+01
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4756114166606E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   5.1182559672501E-01
- (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.3945736322431E-04
+ (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.3945736322432E-04
- (PID.TID 0000.0001) %MON forcing_qnet_max             =   7.9361545319753E+02
+ (PID.TID 0000.0001) %MON forcing_qnet_max             =   7.9361545319752E+02
  (PID.TID 0000.0001) %MON forcing_qnet_min             =  -2.3631242929674E+02
- (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -5.8817586704126E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -5.8817586704128E+00
  (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.2044935518654E+02
  (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0527159003843E+00
  (PID.TID 0000.0001) %MON forcing_qsw_max              =  -0.0000000000000E+00
-Line 3605 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3690 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_qsw_del2             =   6.6740856327453E-01
  (PID.TID 0000.0001) %MON forcing_empmr_max            =   7.7364662397364E-04
  (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.3749339523291E-04
- (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -9.8503774950141E-08
+ (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -9.8503774950161E-08
  (PID.TID 0000.0001) %MON forcing_empmr_sd             =   5.2342343064608E-05
  (PID.TID 0000.0001) %MON forcing_empmr_del2           =   1.0366783585396E-06
  (PID.TID 0000.0001) %MON forcing_fu_max               =   2.4726979082249E-01
-Line 3619 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3704 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_sd                =   7.5645935346452E-02
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.3922805906111E-04
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.6183314724195E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.5089776686688E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.5089776686689E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8122934736114E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8122934736117E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4209496264105E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4209496264109E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9300326567700E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.6016599334840E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6191627730202E-04
-Line 3632 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3717 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582519369E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787197077809E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817301200150E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   8.0187474751627E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   8.0187474751770E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9004918473538E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9004918476018E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 3679 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3764 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
-  cg2d: Sum(rhs),rhsMax =   3.91486139871320E+02  2.29746963060133E+02
+  BULKF_FIELDS_LOAD,     36008 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.2333333333  0.7666666667
- (PID.TID 0000.0001)                    cg2d_init_res =    2.67156784229563E+00
+  EXTERNAL_FIELDS_LOAD,     36008 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.2333333333  0.7666666667
+  cg2d: Sum(rhs),rhsMax =   3.91486139871320E+02  2.29746963060132E+02
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.67156784229858E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.44418240422538E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.44418240421608E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36009
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1111776000000E+09
  (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2403424935864E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1971042154524E+01
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1971042154525E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8227418185239E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.7005245530022E+00
  (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1288524998585E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7452811878665E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6820187944108E-01
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2461895251606E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2461895251605E-04
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   1.2976354693890E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8633532082744E-05
- (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.9519195365314E-01
+ (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.9519195365315E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -2.1128980936852E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.2116387204702E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3481408306787E-02
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.6042541577189E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.7646784032062E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.7646784032063E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7822554338372E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7822554338373E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2228266221177E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2228266221176E-09
  (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6218065611361E-06
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0338994360841E-08
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0338994360840E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9979348113421E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.3852753974119E+00
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.5033482231088E+00
-Line 3720 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3807 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.3967010093795E-04
  (PID.TID 0000.0001) %MON forcing_qnet_max             =   7.9679638952914E+02
  (PID.TID 0000.0001) %MON forcing_qnet_min             =  -2.3593040361628E+02
- (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -5.9424603634831E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -5.9424603634833E+00
  (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.2037409410669E+02
  (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0482997269942E+00
  (PID.TID 0000.0001) %MON forcing_qsw_max              =  -0.0000000000000E+00
-Line 3730 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3817 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_qsw_del2             =   6.6522640536388E-01
  (PID.TID 0000.0001) %MON forcing_empmr_max            =   5.1402201507388E-04
  (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.0316831189022E-04
- (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -3.1409376269531E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -3.1409376269534E-07
  (PID.TID 0000.0001) %MON forcing_empmr_sd             =   4.9947911782859E-05
- (PID.TID 0000.0001) %MON forcing_empmr_del2           =   9.5635737137828E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_del2           =   9.5635737137827E-07
  (PID.TID 0000.0001) %MON forcing_fu_max               =   2.4722174155228E-01
  (PID.TID 0000.0001) %MON forcing_fu_min               =  -2.2809828301201E-01
  (PID.TID 0000.0001) %MON forcing_fu_mean              =  -4.2000491789953E-03
-Line 3744 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3831 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_sd                =   7.5757214244312E-02
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.4079325037298E-04
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.6067140446864E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.5081422691354E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.5081422691355E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8108252635297E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8108252635300E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4193288440915E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4193288440918E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9291659710253E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.5891962052562E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6188692597221E-04
-Line 3757 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3844 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582455302E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787196613584E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817301297664E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   8.0356881165319E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   8.0356881165393E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9156129485757E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9156129486127E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 3804 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3891 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
-  cg2d: Sum(rhs),rhsMax =   3.91457886686040E+02  2.29845851768187E+02
+  BULKF_FIELDS_LOAD,     36009 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.2000000000  0.8000000000
- (PID.TID 0000.0001)                    cg2d_init_res =    2.66058318162855E+00
+  EXTERNAL_FIELDS_LOAD,     36009 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.2000000000  0.8000000000
+  cg2d: Sum(rhs),rhsMax =   3.91457886686041E+02  2.29845851768189E+02
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.66058318164582E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.50267762080705E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.50267762090410E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36010
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1112640000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2270180666403E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2270180666398E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1968648953259E+01
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1968648953258E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8223215967416E-01
- (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.7000282875400E+00
+ (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.7000282875399E+00
  (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1298958273414E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7444086311143E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6777942682888E-01
-Line 3828 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3917 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.2040491683013E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3480602028705E-02
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.6008541238993E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.8405828655294E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.8405828655300E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7813644398365E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7813644398366E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2451958041365E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2451958041335E-09
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6215060670555E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6215060670554E-06
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0336083140411E-08
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0336083140410E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9958085669613E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.4034792315777E+00
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.5034163342807E+00
-Line 3842 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3931 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_salt_min             =   1.8036015600309E+01
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4756115375757E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   5.1181206760147E-01
- (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.3999183516064E-04
+ (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.3999183516063E-04
  (PID.TID 0000.0001) %MON forcing_qnet_max             =   7.9995369355762E+02
  (PID.TID 0000.0001) %MON forcing_qnet_min             =  -2.3555036982624E+02
- (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -5.9980072122519E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -5.9980072122523E+00
  (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.2037486174731E+02
  (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0515774169975E+00
  (PID.TID 0000.0001) %MON forcing_qsw_max              =  -0.0000000000000E+00
-Line 3855 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3944 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_qsw_del2             =   6.6429930669879E-01
  (PID.TID 0000.0001) %MON forcing_empmr_max            =   3.9430549916432E-04
  (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.0800471653405E-04
- (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -5.0339067668266E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -5.0339067668273E-07
  (PID.TID 0000.0001) %MON forcing_empmr_sd             =   4.9099688882812E-05
  (PID.TID 0000.0001) %MON forcing_empmr_del2           =   9.3433498831442E-07
  (PID.TID 0000.0001) %MON forcing_fu_max               =   2.4734156444157E-01
-Line 3869 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 3958 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_sd                =   7.5878609891640E-02
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.4254447090673E-04
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.5947324872674E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.5045355374288E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.5045355374289E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8079202868631E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8079202868636E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4161211752725E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4161211752730E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9282423423910E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.5764730831735E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6185301693979E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.3395905635104E+18
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.2529636027581E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.3965864716742E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.3965864716741E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =  -2.0549865324846E-05
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582389684E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787196075333E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817301306510E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   8.0271970263851E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   8.0271970263847E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9360707199906E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9360707198271E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 3930 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4019 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
   Computing Diagnostic #     23  ETAN         Counter:       1   Parms: SM      M1
-  Computing Diagnostic #     79  oceFWflx     Counter:       1   Parms: SM      U1
+  Computing Diagnostic #     82  oceFWflx     Counter:       1   Parms: SM      U1
   Computing Diagnostic #     26  THETA        Counter:       1   Parms: SMR     MR
   Computing Diagnostic #     27  SALT         Counter:       1   Parms: SMR     MR
   Compute Stats, Diag. #     23  ETAN      vol(   0 ): 3.639E+15  Parms: SM      M1
-Line 3938 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4027 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
   Compute Stats, Diag. #     25  DETADT2   vol(   0 ): 3.639E+15  Parms: SM      M1
   Compute Stats, Diag. #     26  THETA     vol(   0 ): 1.340E+19  Parms: SMR     MR
   Compute Stats, Diag. #     27  SALT      vol(   0 ): 1.340E+19  Parms: SMR     MR
-  Compute Stats, Diag. #     74  CONVADJ   vol(   0 ): 1.315E+19  Parms: SMR     LR
+  Compute Stats, Diag. #     77  CONVADJ   vol(   0 ): 1.315E+19  Parms: SMR     LR
-  Compute Stats, Diag. #     75  oceTAUX   vol(   0 ): 3.111E+15  Parms: UU      U1
+  Compute Stats, Diag. #     78  oceTAUX   vol(   0 ): 3.111E+15  Parms: UU      U1
-  Compute Stats, Diag. #     76  oceTAUY   vol(   0 ): 3.130E+15  Parms: VV      U1
+  Compute Stats, Diag. #     79  oceTAUY   vol(   0 ): 3.130E+15  Parms: VV      U1
-  Compute Stats, Diag. #     84  TRELAX    vol(   0 ): 3.275E+15  Parms: SM      U1
+  Compute Stats, Diag. #     87  TRELAX    vol(   0 ): 3.275E+15  Parms: SM      U1
-  Compute Stats, Diag. #     85  SRELAX    vol(   0 ): 3.275E+15  Parms: SM      U1
+  Compute Stats, Diag. #     88  SRELAX    vol(   0 ): 3.275E+15  Parms: SM      U1
-  Compute Stats, Diag. #     86  surForcT  vol(   0 ): 3.275E+15  Parms: SM      U1
+  Compute Stats, Diag. #     89  surForcT  vol(   0 ): 3.275E+15  Parms: SM      U1
-  Compute Stats, Diag. #     87  surForcS  vol(   0 ): 3.275E+15  Parms: SM      U1
+  Compute Stats, Diag. #     90  surForcS  vol(   0 ): 3.275E+15  Parms: SM      U1
-  Compute Stats, Diag. #     88  TFLUX     vol(   0 ): 3.275E+15  Parms: SM      U1
+  Compute Stats, Diag. #     91  TFLUX     vol(   0 ): 3.275E+15  Parms: SM      U1
-  Compute Stats, Diag. #     89  SFLUX     vol(   0 ): 3.275E+15  Parms: SM      U1
+  Compute Stats, Diag. #     92  SFLUX     vol(   0 ): 3.275E+15  Parms: SM      U1
-  Compute Stats, Diag. #     81  oceQnet   vol(   0 ): 3.275E+15  Parms: SM      U1
+  Compute Stats, Diag. #     84  oceQnet   vol(   0 ): 3.275E+15  Parms: SM      U1
-  Compute Stats, Diag. #     80  oceSflux  vol(   0 ): 3.275E+15  Parms: SM      U1
+  Compute Stats, Diag. #     83  oceSflux  vol(   0 ): 3.275E+15  Parms: SM      U1
-  Compute Stats, Diag. #     79  oceFWflx  vol(   0 ): 3.275E+15  Parms: SM      U1
+  Compute Stats, Diag. #     82  oceFWflx  vol(   0 ): 3.275E+15  Parms: SM      U1
-  Compute Stats, Diag. #    194  SI_Fract  vol(   0 ): 3.639E+15  Parms: SM P    M1
+  Compute Stats, Diag. #    207  SI_Fract  vol(   0 ): 3.639E+15  Parms: SM P    M1
-  Compute Stats, Diag. #    194  SI_Fract  vol(   1 ): 1.289E+15  Parms: SM P    M1
+  Compute Stats, Diag. #    207  SI_Fract  vol(   1 ): 1.289E+15  Parms: SM P    M1
-  Compute Stats, Diag. #    194  SI_Fract  vol(   3 ): 7.935E+14  Parms: SM P    M1
+  Compute Stats, Diag. #    207  SI_Fract  vol(   3 ): 7.935E+14  Parms: SM P    M1
-  Compute Stats, Diag. #    195  SI_Thick  vol(   0 ): 2.483E+14  Parms: SM PC   M1
+  Compute Stats, Diag. #    208  SI_Thick  vol(   0 ): 2.483E+14  Parms: SM PC   M1
-     use Counter Mate  #    194  SI_Fract  vol(   0 ): 3.639E+15 integral 2.483E+14
+     use Counter Mate  #    207  SI_Fract  vol(   0 ): 3.639E+15 integral 2.483E+14
-  Compute Stats, Diag. #    195  SI_Thick  vol(   1 ): 1.189E+14  Parms: SM PC   M1
+  Compute Stats, Diag. #    208  SI_Thick  vol(   1 ): 1.189E+14  Parms: SM PC   M1
-     use Counter Mate  #    194  SI_Fract  vol(   1 ): 1.289E+15 integral 1.189E+14
+     use Counter Mate  #    207  SI_Fract  vol(   1 ): 1.289E+15 integral 1.189E+14
-  Compute Stats, Diag. #    195  SI_Thick  vol(   3 ): 1.294E+14  Parms: SM PC   M1
+  Compute Stats, Diag. #    208  SI_Thick  vol(   3 ): 1.294E+14  Parms: SM PC   M1
-     use Counter Mate  #    194  SI_Fract  vol(   3 ): 7.935E+14 integral 1.294E+14
+     use Counter Mate  #    207  SI_Fract  vol(   3 ): 7.935E+14 integral 1.294E+14
-  Compute Stats, Diag. #    196  SI_SnowH  vol(   0 ): 2.483E+14  Parms: SM PC   M1
+  Compute Stats, Diag. #    209  SI_SnowH  vol(   0 ): 2.483E+14  Parms: SM PC   M1
-     use Counter Mate  #    194  SI_Fract  vol(   0 ): 3.639E+15 integral 2.483E+14
+     use Counter Mate  #    207  SI_Fract  vol(   0 ): 3.639E+15 integral 2.483E+14
-  Compute Stats, Diag. #    196  SI_SnowH  vol(   1 ): 1.189E+14  Parms: SM PC   M1
+  Compute Stats, Diag. #    209  SI_SnowH  vol(   1 ): 1.189E+14  Parms: SM PC   M1
-     use Counter Mate  #    194  SI_Fract  vol(   1 ): 1.289E+15 integral 1.189E+14
+     use Counter Mate  #    207  SI_Fract  vol(   1 ): 1.289E+15 integral 1.189E+14
-  Compute Stats, Diag. #    196  SI_SnowH  vol(   3 ): 1.294E+14  Parms: SM PC   M1
+  Compute Stats, Diag. #    209  SI_SnowH  vol(   3 ): 1.294E+14  Parms: SM PC   M1
-     use Counter Mate  #    194  SI_Fract  vol(   3 ): 7.935E+14 integral 1.294E+14
+     use Counter Mate  #    207  SI_Fract  vol(   3 ): 7.935E+14 integral 1.294E+14
-  Compute Stats, Diag. #    197  SI_Tsrf   vol(   0 ): 2.483E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    210  SI_Tsrf   vol(   0 ): 2.483E+14  Parms: SM  C   M1
-     use Counter Mate  #    194  SI_Fract  vol(   0 ): 3.639E+15 integral 2.483E+14
+     use Counter Mate  #    207  SI_Fract  vol(   0 ): 3.639E+15 integral 2.483E+14
-  Compute Stats, Diag. #    197  SI_Tsrf   vol(   1 ): 1.189E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    210  SI_Tsrf   vol(   1 ): 1.189E+14  Parms: SM  C   M1
-     use Counter Mate  #    194  SI_Fract  vol(   1 ): 1.289E+15 integral 1.189E+14
+     use Counter Mate  #    207  SI_Fract  vol(   1 ): 1.289E+15 integral 1.189E+14
-  Compute Stats, Diag. #    197  SI_Tsrf   vol(   3 ): 1.294E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    210  SI_Tsrf   vol(   3 ): 1.294E+14  Parms: SM  C   M1
-     use Counter Mate  #    194  SI_Fract  vol(   3 ): 7.935E+14 integral 1.294E+14
+     use Counter Mate  #    207  SI_Fract  vol(   3 ): 7.935E+14 integral 1.294E+14
-  Compute Stats, Diag. #    198  SI_Tice1  vol(   0 ): 2.483E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    211  SI_Tice1  vol(   0 ): 2.483E+14  Parms: SM  C   M1
-     use Counter Mate  #    194  SI_Fract  vol(   0 ): 3.639E+15 integral 2.483E+14
+     use Counter Mate  #    207  SI_Fract  vol(   0 ): 3.639E+15 integral 2.483E+14
-  Compute Stats, Diag. #    198  SI_Tice1  vol(   1 ): 1.189E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    211  SI_Tice1  vol(   1 ): 1.189E+14  Parms: SM  C   M1
-     use Counter Mate  #    194  SI_Fract  vol(   1 ): 1.289E+15 integral 1.189E+14
+     use Counter Mate  #    207  SI_Fract  vol(   1 ): 1.289E+15 integral 1.189E+14
-  Compute Stats, Diag. #    198  SI_Tice1  vol(   3 ): 1.294E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    211  SI_Tice1  vol(   3 ): 1.294E+14  Parms: SM  C   M1
-     use Counter Mate  #    194  SI_Fract  vol(   3 ): 7.935E+14 integral 1.294E+14
+     use Counter Mate  #    207  SI_Fract  vol(   3 ): 7.935E+14 integral 1.294E+14
-  Compute Stats, Diag. #    199  SI_Tice2  vol(   0 ): 2.483E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    212  SI_Tice2  vol(   0 ): 2.483E+14  Parms: SM  C   M1
-     use Counter Mate  #    194  SI_Fract  vol(   0 ): 3.639E+15 integral 2.483E+14
+     use Counter Mate  #    207  SI_Fract  vol(   0 ): 3.639E+15 integral 2.483E+14
-  Compute Stats, Diag. #    199  SI_Tice2  vol(   1 ): 1.189E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    212  SI_Tice2  vol(   1 ): 1.189E+14  Parms: SM  C   M1
-     use Counter Mate  #    194  SI_Fract  vol(   1 ): 1.289E+15 integral 1.189E+14
+     use Counter Mate  #    207  SI_Fract  vol(   1 ): 1.289E+15 integral 1.189E+14
-  Compute Stats, Diag. #    199  SI_Tice2  vol(   3 ): 1.294E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    212  SI_Tice2  vol(   3 ): 1.294E+14  Parms: SM  C   M1
-     use Counter Mate  #    194  SI_Fract  vol(   3 ): 7.935E+14 integral 1.294E+14
+     use Counter Mate  #    207  SI_Fract  vol(   3 ): 7.935E+14 integral 1.294E+14
-  Compute Stats, Diag. #    200  SI_Qice1  vol(   0 ): 1.136E+15  Parms: SM  C   M1
+  Compute Stats, Diag. #    213  SI_Qice1  vol(   0 ): 1.136E+15  Parms: SM  C   M1
-     use Counter Mate  #    195  SI_Thick  vol(   0 ): 2.483E+14 integral 1.136E+15
+     use Counter Mate  #    208  SI_Thick  vol(   0 ): 2.483E+14 integral 1.136E+15
-  Compute Stats, Diag. #    200  SI_Qice1  vol(   1 ): 6.824E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    213  SI_Qice1  vol(   1 ): 6.824E+14  Parms: SM  C   M1
-     use Counter Mate  #    195  SI_Thick  vol(   1 ): 1.189E+14 integral 6.824E+14
+     use Counter Mate  #    208  SI_Thick  vol(   1 ): 1.189E+14 integral 6.824E+14
-  Compute Stats, Diag. #    200  SI_Qice1  vol(   3 ): 4.532E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    213  SI_Qice1  vol(   3 ): 4.532E+14  Parms: SM  C   M1
-     use Counter Mate  #    195  SI_Thick  vol(   3 ): 1.294E+14 integral 4.532E+14
+     use Counter Mate  #    208  SI_Thick  vol(   3 ): 1.294E+14 integral 4.532E+14
-  Compute Stats, Diag. #    201  SI_Qice2  vol(   0 ): 1.136E+15  Parms: SM  C   M1
+  Compute Stats, Diag. #    214  SI_Qice2  vol(   0 ): 1.136E+15  Parms: SM  C   M1
-     use Counter Mate  #    195  SI_Thick  vol(   0 ): 2.483E+14 integral 1.136E+15
+     use Counter Mate  #    208  SI_Thick  vol(   0 ): 2.483E+14 integral 1.136E+15
-  Compute Stats, Diag. #    201  SI_Qice2  vol(   1 ): 6.824E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    214  SI_Qice2  vol(   1 ): 6.824E+14  Parms: SM  C   M1
-     use Counter Mate  #    195  SI_Thick  vol(   1 ): 1.189E+14 integral 6.824E+14
+     use Counter Mate  #    208  SI_Thick  vol(   1 ): 1.189E+14 integral 6.824E+14
-  Compute Stats, Diag. #    201  SI_Qice2  vol(   3 ): 4.532E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    214  SI_Qice2  vol(   3 ): 4.532E+14  Parms: SM  C   M1
-     use Counter Mate  #    195  SI_Thick  vol(   3 ): 1.294E+14 integral 4.532E+14
+     use Counter Mate  #    208  SI_Thick  vol(   3 ): 1.294E+14 integral 4.532E+14
-  Compute Stats, Diag. #    205  SIsnwPrc  vol(   0 ): 2.483E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    218  SIsnwPrc  vol(   0 ): 2.483E+14  Parms: SM  C   M1
-     use Counter Mate  #    203  SI_FrcFx  vol(   0 ): 3.639E+15 integral 2.483E+14
+     use Counter Mate  #    216  SI_FrcFx  vol(   0 ): 3.639E+15 integral 2.483E+14
-  Compute Stats, Diag. #    205  SIsnwPrc  vol(   1 ): 1.189E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    218  SIsnwPrc  vol(   1 ): 1.189E+14  Parms: SM  C   M1
-     use Counter Mate  #    203  SI_FrcFx  vol(   1 ): 1.289E+15 integral 1.189E+14
+     use Counter Mate  #    216  SI_FrcFx  vol(   1 ): 1.289E+15 integral 1.189E+14
-  Compute Stats, Diag. #    205  SIsnwPrc  vol(   3 ): 1.294E+14  Parms: SM  C   M1
+  Compute Stats, Diag. #    218  SIsnwPrc  vol(   3 ): 1.294E+14  Parms: SM  C   M1
-     use Counter Mate  #    203  SI_FrcFx  vol(   3 ): 7.935E+14 integral 1.294E+14
+     use Counter Mate  #    216  SI_FrcFx  vol(   3 ): 7.935E+14 integral 1.294E+14
-  Compute Stats, Diag. #    204  SIalbedo  vol(   0 ): 2.483E+14  Parms: SM PC   M1
+  Compute Stats, Diag. #    217  SIalbedo  vol(   0 ): 2.483E+14  Parms: SM PC   M1
-     use Counter Mate  #    203  SI_FrcFx  vol(   0 ): 3.639E+15 integral 2.483E+14
+     use Counter Mate  #    216  SI_FrcFx  vol(   0 ): 3.639E+15 integral 2.483E+14
-  Compute Stats, Diag. #    204  SIalbedo  vol(   1 ): 1.189E+14  Parms: SM PC   M1
+  Compute Stats, Diag. #    217  SIalbedo  vol(   1 ): 1.189E+14  Parms: SM PC   M1
-     use Counter Mate  #    203  SI_FrcFx  vol(   1 ): 1.289E+15 integral 1.189E+14
+     use Counter Mate  #    216  SI_FrcFx  vol(   1 ): 1.289E+15 integral 1.189E+14
-  Compute Stats, Diag. #    204  SIalbedo  vol(   3 ): 1.294E+14  Parms: SM PC   M1
+  Compute Stats, Diag. #    217  SIalbedo  vol(   3 ): 1.294E+14  Parms: SM PC   M1
-     use Counter Mate  #    203  SI_FrcFx  vol(   3 ): 7.935E+14 integral 1.294E+14
+     use Counter Mate  #    216  SI_FrcFx  vol(   3 ): 7.935E+14 integral 1.294E+14
-  Compute Stats, Diag. #    202  SIsnwAge  vol(   0 ): 3.639E+15  Parms: SM P    M1
+  Compute Stats, Diag. #    215  SIsnwAge  vol(   0 ): 3.639E+15  Parms: SM P    M1
-  Compute Stats, Diag. #    202  SIsnwAge  vol(   1 ): 1.289E+15  Parms: SM P    M1
+  Compute Stats, Diag. #    215  SIsnwAge  vol(   1 ): 1.289E+15  Parms: SM P    M1
-  Compute Stats, Diag. #    202  SIsnwAge  vol(   3 ): 7.935E+14  Parms: SM P    M1
+  Compute Stats, Diag. #    215  SIsnwAge  vol(   3 ): 7.935E+14  Parms: SM P    M1
-  Compute Stats, Diag. #    208  SIflx2oc  vol(   0 ): 3.639E+15  Parms: SM      M1
+  Compute Stats, Diag. #    221  SIflx2oc  vol(   0 ): 3.639E+15  Parms: SM      M1
-  Compute Stats, Diag. #    208  SIflx2oc  vol(   1 ): 1.289E+15  Parms: SM      M1
+  Compute Stats, Diag. #    221  SIflx2oc  vol(   1 ): 1.289E+15  Parms: SM      M1
-  Compute Stats, Diag. #    208  SIflx2oc  vol(   3 ): 7.935E+14  Parms: SM      M1
+  Compute Stats, Diag. #    221  SIflx2oc  vol(   3 ): 7.935E+14  Parms: SM      M1
-  Compute Stats, Diag. #    209  SIfrw2oc  vol(   0 ): 3.639E+15  Parms: SM      M1
+  Compute Stats, Diag. #    222  SIfrw2oc  vol(   0 ): 3.639E+15  Parms: SM      M1
-  Compute Stats, Diag. #    209  SIfrw2oc  vol(   1 ): 1.289E+15  Parms: SM      M1
+  Compute Stats, Diag. #    222  SIfrw2oc  vol(   1 ): 1.289E+15  Parms: SM      M1
-  Compute Stats, Diag. #    209  SIfrw2oc  vol(   3 ): 7.935E+14  Parms: SM      M1
+  Compute Stats, Diag. #    222  SIfrw2oc  vol(   3 ): 7.935E+14  Parms: SM      M1
-  Compute Stats, Diag. #    210  SIsaltFx  vol(   0 ): 3.639E+15  Parms: SM      M1
+  Compute Stats, Diag. #    223  SIsaltFx  vol(   0 ): 3.639E+15  Parms: SM      M1
-  Compute Stats, Diag. #    210  SIsaltFx  vol(   1 ): 1.289E+15  Parms: SM      M1
+  Compute Stats, Diag. #    223  SIsaltFx  vol(   1 ): 1.289E+15  Parms: SM      M1
-  Compute Stats, Diag. #    210  SIsaltFx  vol(   3 ): 7.935E+14  Parms: SM      M1
+  Compute Stats, Diag. #    223  SIsaltFx  vol(   3 ): 7.935E+14  Parms: SM      M1
-  Compute Stats, Diag. #    206  SIflxAtm  vol(   0 ): 3.639E+15  Parms: SM      M1
+  Compute Stats, Diag. #    219  SIflxAtm  vol(   0 ): 3.639E+15  Parms: SM      M1
-  Compute Stats, Diag. #    206  SIflxAtm  vol(   1 ): 1.289E+15  Parms: SM      M1
+  Compute Stats, Diag. #    219  SIflxAtm  vol(   1 ): 1.289E+15  Parms: SM      M1
-  Compute Stats, Diag. #    206  SIflxAtm  vol(   3 ): 7.935E+14  Parms: SM      M1
+  Compute Stats, Diag. #    219  SIflxAtm  vol(   3 ): 7.935E+14  Parms: SM      M1
-  Compute Stats, Diag. #    207  SIfrwAtm  vol(   0 ): 3.639E+15  Parms: SM      M1
+  Compute Stats, Diag. #    220  SIfrwAtm  vol(   0 ): 3.639E+15  Parms: SM      M1
-  Compute Stats, Diag. #    207  SIfrwAtm  vol(   1 ): 1.289E+15  Parms: SM      M1
+  Compute Stats, Diag. #    220  SIfrwAtm  vol(   1 ): 1.289E+15  Parms: SM      M1
-  Compute Stats, Diag. #    207  SIfrwAtm  vol(   3 ): 7.935E+14  Parms: SM      M1
+  Compute Stats, Diag. #    220  SIfrwAtm  vol(   3 ): 7.935E+14  Parms: SM      M1
+  BULKF_FIELDS_LOAD,     36010 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.1666666667  0.8333333333
+  EXTERNAL_FIELDS_LOAD,     36010 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.1666666667  0.8333333333
   cg2d: Sum(rhs),rhsMax =   3.91430150858270E+02  2.29949386798939E+02
- (PID.TID 0000.0001)                    cg2d_init_res =    2.65469956429210E+00
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.65469956429600E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.65998377899878E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.65998377904557E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36011
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1113504000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2134361015207E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.2134361015204E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1966126084823E+01
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1966126084824E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8219090698402E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.6995361430003E+00
  (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1308339919130E-02
-Line 4049 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4140 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.1969329645852E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3479353123246E-02
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.5976065281887E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.9147430895795E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.9147430895805E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7800718769998E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7800718770000E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2146281984222E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2146281984232E-09
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6213071088479E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6213071088481E-06
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0332712330234E-08
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0332712330237E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9936870413488E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.4216752708134E+00
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.5034854158603E+00
  (PID.TID 0000.0001) %MON dynstat_theta_sd             =   4.7497764225123E+00
- (PID.TID 0000.0001) %MON dynstat_theta_del2           =   2.9662215685234E-03
+ (PID.TID 0000.0001) %MON dynstat_theta_del2           =   2.9662215685233E-03
  (PID.TID 0000.0001) %MON dynstat_salt_max             =   4.0836863540121E+01
  (PID.TID 0000.0001) %MON dynstat_salt_min             =   1.8036636154289E+01
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4756115946852E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   5.1180506899605E-01
- (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4043375234564E-04
+ (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4043375234563E-04
  (PID.TID 0000.0001) %MON forcing_qnet_max             =   8.0308307135363E+02
  (PID.TID 0000.0001) %MON forcing_qnet_min             =  -2.3517228105122E+02
- (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.1221141630004E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.1221141630012E+00
  (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.2034715145913E+02
  (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0552800356491E+00
  (PID.TID 0000.0001) %MON forcing_qsw_max              =  -0.0000000000000E+00
-Line 4076 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4167 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_qsw_del2             =   6.6139123592579E-01
  (PID.TID 0000.0001) %MON forcing_empmr_max            =   7.5874665778011E-04
  (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.1398312633061E-04
- (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -4.9417285066546E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -4.9417285066562E-07
  (PID.TID 0000.0001) %MON forcing_empmr_sd             =   4.9699306820229E-05
  (PID.TID 0000.0001) %MON forcing_empmr_del2           =   1.0426714424877E-06
  (PID.TID 0000.0001) %MON forcing_fu_max               =   2.4760593894036E-01
-Line 4090 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4181 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_sd                =   7.6010073816084E-02
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.4447772305236E-04
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.5823789092782E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4982102173211E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4982102173212E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8037060442559E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.8037060442566E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4114675933775E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4114675933783E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9272640480159E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.5634894025277E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6181510726792E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.3395905788017E+18
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.2580641306409E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.3999207920377E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.3999207920376E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =  -2.0549865324846E-05
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582321038E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787195461354E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817301277590E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.9380565196205E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.9380565196105E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9276035065014E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9276035065263E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 4154 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4245 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
   Compute Stats, Diag. #     23  ETAN      vol(   0 ): 3.639E+14  Parms: SM      M1
   Compute Stats, Diag. #     26  THETA     vol(   0 ): 1.340E+18  Parms: SMR     MR
   Compute Stats, Diag. #     27  SALT      vol(   0 ): 1.340E+18  Parms: SMR     MR
-  cg2d: Sum(rhs),rhsMax =   3.91394363526222E+02  2.30056289531097E+02
+  BULKF_FIELDS_LOAD,     36011 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.1333333333  0.8666666667
- (PID.TID 0000.0001)                    cg2d_init_res =    2.65157687489623E+00
+  EXTERNAL_FIELDS_LOAD,     36011 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.1333333333  0.8666666667
+  cg2d: Sum(rhs),rhsMax =   3.91394363526221E+02  2.30056289531097E+02
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.65157687492423E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.67356832308578E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.67356832311039E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36012
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1114368000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1998491745812E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1998491745818E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1963458820621E+01
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1963458820619E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8213767894510E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.6990481496258E+00
  (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1317227171078E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7409176188108E-01
+ (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7409176188109E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6689498520935E-01
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2311896001014E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2311896001013E-04
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   1.2970908625051E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8567847521478E-05
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8567847521477E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.9475928032105E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -2.1143657391300E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.1904277143128E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3477673202936E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.5945146016758E-05
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.5945146016759E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.9861083579740E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   8.9861083579752E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7784243895153E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7784243895155E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2395274180234E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2395274180248E-09
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6212218443041E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6212218443044E-06
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0328280881357E-08
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0328280881362E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9915750600271E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.4398629181247E+00
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.5035551049237E+00
-Line 4192 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4285 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_salt_min             =   1.8037366081844E+01
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4756116432199E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   5.1179759251744E-01
- (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4093079586567E-04
+ (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4093079586562E-04
  (PID.TID 0000.0001) %MON forcing_qnet_max             =   8.0618610184580E+02
  (PID.TID 0000.0001) %MON forcing_qnet_min             =  -2.3479608753032E+02
- (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.1971365841535E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.1971365841552E+00
  (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.2039635648387E+02
  (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0682566711480E+00
  (PID.TID 0000.0001) %MON forcing_qsw_max              =  -0.0000000000000E+00
-Line 4205 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4298 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_qsw_del2             =   6.5859164944186E-01
  (PID.TID 0000.0001) %MON forcing_empmr_max            =   7.4979533925785E-04
  (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.1808628194163E-04
- (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -6.3762754954440E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -6.3762754954479E-07
- (PID.TID 0000.0001) %MON forcing_empmr_sd             =   4.9425833552214E-05
+ (PID.TID 0000.0001) %MON forcing_empmr_sd             =   4.9425833552215E-05
  (PID.TID 0000.0001) %MON forcing_empmr_del2           =   1.0347184003415E-06
  (PID.TID 0000.0001) %MON forcing_fu_max               =   2.4787031343914E-01
  (PID.TID 0000.0001) %MON forcing_fu_min               =  -2.3765894517478E-01
-Line 4219 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4312 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_sd                =   7.6151553873676E-02
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.4658872547121E-04
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.5696480452679E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4892826481489E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4892826481490E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7983346133631E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7983346133638E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4055361383928E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.4055361383936E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9262743252327E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.5502464042277E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6177379298924E-04
-Line 4232 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4325 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582249562E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787194883184E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817301417921E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.8714156026242E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.8714156026132E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9411253968855E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9411253968684E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 4279 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4372 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
-  cg2d: Sum(rhs),rhsMax =   3.91348280743500E+02  2.30165112012775E+02
+  BULKF_FIELDS_LOAD,     36012 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.1000000000  0.9000000000
- (PID.TID 0000.0001)                    cg2d_init_res =    2.62880540543944E+00
+  EXTERNAL_FIELDS_LOAD,     36012 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.1000000000  0.9000000000
+  cg2d: Sum(rhs),rhsMax =   3.91348280743500E+02  2.30165112012776E+02
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.62880540546608E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.60159325000787E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.60159325001877E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36013
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1115232000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1867492710061E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1867492710067E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1960638945423E+01
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1960638945424E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8206913803937E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.6985631424223E+00
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1326332475376E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1326332475375E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7383923997945E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6643270121657E-01
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2249828403734E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2249828403735E-04
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   1.2969437724762E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8549735706081E-05
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8549735706080E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.9449490032690E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -2.1147499747559E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.1845205516326E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3475589763592E-02
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.5915816178033E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.0537526277640E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.0537526277643E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7764753324981E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7764753324982E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2710409310934E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2710409310936E-09
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6212245335611E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6212245335613E-06
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0322289410234E-08
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0322289410239E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9894770906997E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.4580265526266E+00
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.5036251061250E+00
-Line 4317 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4412 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_salt_min             =   1.8038470005667E+01
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4756116802659E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   5.1178903999054E-01
- (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4207156649405E-04
+ (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4207156649406E-04
  (PID.TID 0000.0001) %MON forcing_qnet_max             =   8.0926481943961E+02
  (PID.TID 0000.0001) %MON forcing_qnet_min             =  -2.3442174029807E+02
- (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.2280625490105E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.2280625490118E+00
  (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.2045340658252E+02
- (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0657631772278E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0657631772275E+00
  (PID.TID 0000.0001) %MON forcing_qsw_max              =  -0.0000000000000E+00
  (PID.TID 0000.0001) %MON forcing_qsw_min              =  -3.3950810942532E+02
  (PID.TID 0000.0001) %MON forcing_qsw_mean             =  -1.8661640031248E+02
-Line 4330 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4425 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_qsw_del2             =   6.5615450436449E-01
  (PID.TID 0000.0001) %MON forcing_empmr_max            =   4.0677668889919E-04
  (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.2063164610717E-04
- (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -8.2106293319219E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -8.2106293319238E-07
  (PID.TID 0000.0001) %MON forcing_empmr_sd             =   4.8530127302672E-05
- (PID.TID 0000.0001) %MON forcing_empmr_del2           =   9.3974825927882E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_del2           =   9.3974825927881E-07
  (PID.TID 0000.0001) %MON forcing_fu_max               =   2.4813468793793E-01
  (PID.TID 0000.0001) %MON forcing_fu_min               =  -2.4114197747619E-01
  (PID.TID 0000.0001) %MON forcing_fu_mean              =  -4.0261334360615E-03
-Line 4345 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4440 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.4887295503804E-04
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.5565367943406E-02
  (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4777586966682E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7919799519936E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7919799519943E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.3985189134404E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.3985189134411E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9252927369980E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.5367472025258E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6172970201004E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.3395906131820E+18
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.2683575705842E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.4061523610238E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.4061523610237E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =  -2.0549865324846E-05
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582175931E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787194241417E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817301510240E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.7982344030099E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.7982344030302E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9600821613348E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9600821613194E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 4404 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4499 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
-  cg2d: Sum(rhs),rhsMax =   3.91298853601252E+02  2.30274254068618E+02
+  BULKF_FIELDS_LOAD,     36013 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.0666666667  0.9333333333
- (PID.TID 0000.0001)                    cg2d_init_res =    2.60691017136006E+00
+  EXTERNAL_FIELDS_LOAD,     36013 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.0666666667  0.9333333333
+  cg2d: Sum(rhs),rhsMax =   3.91298853601251E+02  2.30274254068619E+02
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.60691017136230E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.55449881168837E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.55449881164770E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36014
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1116096000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1737406658910E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1737406658917E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1957664827300E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8199562292381E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.6980787411116E+00
  (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1333899740921E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7354201187243E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6595693345230E-01
- (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2179966155015E-04
+ (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.2179966155014E-04
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   1.2968171049948E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8534012499149E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.9417552573826E-01
-Line 4428 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4525 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.1790641177427E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3473142215923E-02
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.5888122688380E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.1166675438452E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.1166675438450E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7742754454093E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7742754454094E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2779393595526E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2779393595514E-09
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6211943184065E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6211943184067E-06
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0313154771110E-08
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0313154771115E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9876343281487E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.4761606544430E+00
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.5036957841209E+00
-Line 4445 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4542 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4238605565147E-04
  (PID.TID 0000.0001) %MON forcing_qnet_max             =   8.1232886498818E+02
  (PID.TID 0000.0001) %MON forcing_qnet_min             =  -2.3404919253651E+02
- (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.3145784638051E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.3145784638054E+00
  (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.2058917960438E+02
  (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0786348712968E+00
  (PID.TID 0000.0001) %MON forcing_qsw_max              =  -0.0000000000000E+00
-Line 4469 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4566 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_sd                =   7.6464336061349E-02
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.5132568855654E-04
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.5430431180769E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4638376082911E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4638376082910E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7848074899381E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7848074899386E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.3905985873335E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.3905985873339E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9243178654274E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.5229955589623E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6168347956035E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.3395906381231E+18
- (PID.TID 0000.0001) %MON vort_r_min                   =  -1.2734979596035E-06
+ (PID.TID 0000.0001) %MON vort_r_min                   =  -1.2734979596036E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.4089723554284E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.4089723554283E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =  -2.0549865324846E-05
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582101505E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787193509048E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817301427182E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.7079818895513E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.7079818895520E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9614681891076E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9614681890864E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 4516 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4613 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON thSI_Tic1_ave_N              =  -1.4138017736090E+01
  (PID.TID 0000.0001) %MON thSI_Tic1_min_S              =  -1.4660760778607E+01
  (PID.TID 0000.0001) %MON thSI_Tic1_min_N              =  -2.1501493250135E+01
- (PID.TID 0000.0001) %MON thSI_Tic1_max_S              =  -1.6093929460901E-01
+ (PID.TID 0000.0001) %MON thSI_Tic1_max_S              =  -1.6093929460902E-01
  (PID.TID 0000.0001) %MON thSI_Tic1_max_N              =  -8.6662576935006E-01
  (PID.TID 0000.0001) %MON thSI_Tic2_ave_G              =  -3.9102525212702E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_ave_S              =  -3.0534402067594E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_ave_N              =  -5.1608571068800E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_min_S              =  -8.3951639892979E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_min_N              =  -7.8438840895258E+00
- (PID.TID 0000.0001) %MON thSI_Tic2_max_S              =  -9.6811706984241E-01
+ (PID.TID 0000.0001) %MON thSI_Tic2_max_S              =  -9.6811706984242E-01
- (PID.TID 0000.0001) %MON thSI_Tic2_max_N              =  -1.3623327669032E+00
+ (PID.TID 0000.0001) %MON thSI_Tic2_max_N              =  -1.3623327669031E+00
  (PID.TID 0000.0001) %MON thSI_TotEnerg_G              =  -3.8079763240338E+22
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
+  BULKF_FIELDS_LOAD,     36014 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.0333333333  0.9666666667
+  EXTERNAL_FIELDS_LOAD,     36014 : iP,iLd,i0,i1=   12    1   12    1 ; Wght=  0.0333333333  0.9666666667
   cg2d: Sum(rhs),rhsMax =   3.91244729458466E+02  2.30381959125033E+02
- (PID.TID 0000.0001)                    cg2d_init_res =    2.58048674627573E+00
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.58048674628469E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.49118874123129E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.49118874126113E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36015
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1116960000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1611388842789E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1611388842791E-01
  (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1954543683474E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8191512175723E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.6975911317112E+00
-Line 4553 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4652 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.1737971311213E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3470381798103E-02
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.5862052890844E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.1739669251881E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.1739669251871E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7718879980210E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7718879980212E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2817932423782E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2817932423771E-09
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6212282543157E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6212282543160E-06
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0302273291719E-08
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0302273291725E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9861539224238E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.4942552966538E+00
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.5037669637992E+00
  (PID.TID 0000.0001) %MON dynstat_theta_sd             =   4.7503177721078E+00
- (PID.TID 0000.0001) %MON dynstat_theta_del2           =   2.9804132318125E-03
+ (PID.TID 0000.0001) %MON dynstat_theta_del2           =   2.9804132318124E-03
  (PID.TID 0000.0001) %MON dynstat_salt_max             =   4.0840584655423E+01
  (PID.TID 0000.0001) %MON dynstat_salt_min             =   1.8093809723115E+01
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4756117465206E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   5.1177485228369E-01
- (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4266011393589E-04
+ (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4266011393591E-04
  (PID.TID 0000.0001) %MON forcing_qnet_max             =   8.1538008368983E+02
  (PID.TID 0000.0001) %MON forcing_qnet_min             =  -2.3367839909710E+02
- (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.3824683634678E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.3824683634661E+00
  (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.2067440498197E+02
  (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0862070474926E+00
  (PID.TID 0000.0001) %MON forcing_qsw_max              =  -0.0000000000000E+00
-Line 4580 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4679 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_qsw_del2             =   6.5229829940897E-01
  (PID.TID 0000.0001) %MON forcing_empmr_max            =   2.6780416732051E-04
  (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.2195583800067E-04
- (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -9.6433689129306E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -9.6433689129253E-07
  (PID.TID 0000.0001) %MON forcing_empmr_sd             =   4.8116623436947E-05
- (PID.TID 0000.0001) %MON forcing_empmr_del2           =   9.0292527184252E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_del2           =   9.0292527184253E-07
  (PID.TID 0000.0001) %MON forcing_fu_max               =   2.4866343693550E-01
  (PID.TID 0000.0001) %MON forcing_fu_min               =  -2.4810804207900E-01
  (PID.TID 0000.0001) %MON forcing_fu_mean              =  -3.9391755645947E-03
-Line 4593 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4692 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_mean              =  -1.3739459161841E-02
  (PID.TID 0000.0001) %MON forcing_fv_sd                =   7.6635516474945E-02
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.5394204351914E-04
- (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.5291652088080E-02
+ (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.5291652088081E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4476919253969E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4476919253968E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7770235105288E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7770235105293E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.3820029943965E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.3820029943971E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9233537690169E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.5089949808663E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6163577938540E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.3395906648743E+18
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.2786015054104E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.4115440661905E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.4115440661904E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =  -2.0549865324846E-05
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259582029271E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787192758146E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817301331894E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.6049992196079E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.6049992196000E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9630823442856E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9630823443357E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 4642 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4741 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON thSI_Tic1_min_S              =  -1.4529752898170E+01
  (PID.TID 0000.0001) %MON thSI_Tic1_min_N              =  -2.1566733407177E+01
  (PID.TID 0000.0001) %MON thSI_Tic1_max_S              =  -1.5593667490538E-01
- (PID.TID 0000.0001) %MON thSI_Tic1_max_N              =  -8.9121204477310E-01
+ (PID.TID 0000.0001) %MON thSI_Tic1_max_N              =  -8.9121204477309E-01
  (PID.TID 0000.0001) %MON thSI_Tic2_ave_G              =  -3.9203009705686E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_ave_S              =  -3.0494517294379E+00
- (PID.TID 0000.0001) %MON thSI_Tic2_ave_N              =  -5.1882722150832E+00
+ (PID.TID 0000.0001) %MON thSI_Tic2_ave_N              =  -5.1882722150831E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_min_S              =  -8.3614740129122E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_min_N              =  -7.8697514121839E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_max_S              =  -9.9312911587937E-01
-Line 4654 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4753 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
- S/R BULKF_FIELDS_LOAD: Reading new data:    1    2     36015  3.111696000000E+09
+  BULKF_FIELDS_LOAD,     36015 : iP,iLd,i0,i1=   12    1    1    2 ; Wght=  1.0000000000  0.0000000000
+  BULKF_FIELDS_LOAD, it=     36015 : Reading new data, i0,i1=    1    2 (prev=   12    1 )
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_tair_cs.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_tair_cs.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_qair_cs.bin
-Line 4667 
 S/R BULKF_FIELDS_LOAD: Reading new data:
+Line 4767 
 S/R BULKF_FIELDS_LOAD: Reading new data:
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_windspeed_cs.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_pr_scal_cs.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: ncep_pr_scal_cs.bin
- S/R EXTERNAL_FIELDS_LOAD: Reading new data:    1    2     36015  3.111696000000E+09
+  EXTERNAL_FIELDS_LOAD,     36015 : iP,iLd,i0,i1=   12    1    1    2 ; Wght=  1.0000000000  0.0000000000
+  EXTERNAL_FIELDS_LOAD, it=     36015 : Reading new data, i0,i1=    1    2 (prev=   12    1 )
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: trenberth_taux.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: trenberth_taux.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: trenberth_tauy.bin
-Line 4676 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4777 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: lev_surfT_cs_12m.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: lev_surfS_cs_12m.bin
  (PID.TID 0000.0001)  MDS_READ_FIELD: opening global file: lev_surfS_cs_12m.bin
-  cg2d: Sum(rhs),rhsMax =   3.91185325996692E+02  2.30486310409112E+02
+  cg2d: Sum(rhs),rhsMax =   3.91185325996692E+02  2.30486310409111E+02
- (PID.TID 0000.0001)                    cg2d_init_res =    2.53105133743304E+00
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.53105133744283E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.33592759075956E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.33592759064611E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36016
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1117824000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1488498299476E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1488498299466E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1951284945205E+01
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1951284945204E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8182676843235E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.6970909606190E+00
  (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1348455677511E-02
-Line 4697 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4798 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8510967564953E-05
  (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.9338801491522E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -2.1156197882242E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.1684625488476E-04
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.1684625488475E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3467369425564E-02
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.5837601101632E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.2252842315457E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.2252842315442E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7694033076195E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7694033076198E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2898062539275E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2898062539274E-09
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6212226293989E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6212226293993E-06
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0288248466499E-08
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0288248466507E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9846689842694E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.5122955476190E+00
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.5038383551257E+00
  (PID.TID 0000.0001) %MON dynstat_theta_sd             =   4.7504717729839E+00
- (PID.TID 0000.0001) %MON dynstat_theta_del2           =   2.9850292200184E-03
+ (PID.TID 0000.0001) %MON dynstat_theta_del2           =   2.9850292200182E-03
  (PID.TID 0000.0001) %MON dynstat_salt_max             =   4.0841190697168E+01
  (PID.TID 0000.0001) %MON dynstat_salt_min             =   1.8095233471110E+01
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4756117723137E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   5.1176820749556E-01
- (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4334743837496E-04
+ (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4334743837494E-04
  (PID.TID 0000.0001) %MON forcing_qnet_max             =   8.1841713990377E+02
  (PID.TID 0000.0001) %MON forcing_qnet_min             =  -2.3330931537991E+02
- (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.4213345363429E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.4213345363396E+00
- (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.2080732031986E+02
+ (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.2080732031987E+02
- (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0946546910177E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0946546910178E+00
  (PID.TID 0000.0001) %MON forcing_qsw_max              =  -0.0000000000000E+00
  (PID.TID 0000.0001) %MON forcing_qsw_min              =  -3.3889598148373E+02
  (PID.TID 0000.0001) %MON forcing_qsw_mean             =  -1.8664899329531E+02
  (PID.TID 0000.0001) %MON forcing_qsw_sd               =   7.8879491330213E+01
  (PID.TID 0000.0001) %MON forcing_qsw_del2             =   6.4977660420779E-01
  (PID.TID 0000.0001) %MON forcing_empmr_max            =   2.3999387230276E-04
- (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.2098824636312E-04
+ (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.2098824636313E-04
- (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -1.0583992043643E-06
+ (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -1.0583992043635E-06
  (PID.TID 0000.0001) %MON forcing_empmr_sd             =   4.7942961153033E-05
- (PID.TID 0000.0001) %MON forcing_empmr_del2           =   8.8381054011945E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_del2           =   8.8381054011946E-07
  (PID.TID 0000.0001) %MON forcing_fu_max               =   2.4892781143428E-01
  (PID.TID 0000.0001) %MON forcing_fu_min               =  -2.5159107438040E-01
  (PID.TID 0000.0001) %MON forcing_fu_mean              =  -3.8956966288612E-03
-Line 4741 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4842 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_sd                =   7.6816469815204E-02
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.5671701728584E-04
  (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.5149010661783E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4295111209769E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4295111209766E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7689224822011E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7689224822018E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.3730573652917E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.3730573652925E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9223736688017E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.4947482410815E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6158724773682E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.3395906941676E+18
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.2836458756132E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.4138278528544E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.4138278528543E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =  -2.0549865324846E-05
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259581961757E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787191973045E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817301136683E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.5167045854505E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.5167045854573E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9657032865232E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9657032865468E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 4789 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4890 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON thSI_Tic1_min_S              =  -1.4400240726241E+01
  (PID.TID 0000.0001) %MON thSI_Tic1_min_N              =  -2.1666976584978E+01
  (PID.TID 0000.0001) %MON thSI_Tic1_max_S              =  -1.5122458804757E-01
- (PID.TID 0000.0001) %MON thSI_Tic1_max_N              =  -9.1696099908928E-01
+ (PID.TID 0000.0001) %MON thSI_Tic1_max_N              =  -9.1696099908927E-01
  (PID.TID 0000.0001) %MON thSI_Tic2_ave_G              =  -3.9301348142310E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_ave_S              =  -3.0453522899102E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_ave_N              =  -5.2152772665678E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_min_S              =  -8.3266000794908E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_min_N              =  -7.8948015172622E+00
- (PID.TID 0000.0001) %MON thSI_Tic2_max_S              =  -9.6594499550928E-01
+ (PID.TID 0000.0001) %MON thSI_Tic2_max_S              =  -9.6594499550929E-01
  (PID.TID 0000.0001) %MON thSI_Tic2_max_N              =  -1.3897142727310E+00
  (PID.TID 0000.0001) %MON thSI_TotEnerg_G              =  -3.8054401941710E+22
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
+  BULKF_FIELDS_LOAD,     36016 : iP,iLd,i0,i1=    1    2    1    2 ; Wght=  0.9666666667  0.0333333333
+  EXTERNAL_FIELDS_LOAD,     36016 : iP,iLd,i0,i1=    1    2    1    2 ; Wght=  0.9666666667  0.0333333333
   cg2d: Sum(rhs),rhsMax =   3.91137986212888E+02  2.30586335088459E+02
- (PID.TID 0000.0001)                    cg2d_init_res =    2.50737453752571E+00
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.50737453752080E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.29878022138561E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.29878022130434E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36017
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1118688000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1382313292827E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1382313292837E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1948141799486E+01
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1948141799485E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8175635793482E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.6966340017210E+00
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1356280157183E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1356280157182E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7257288994703E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6445237182875E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.1939423111846E-04
-Line 4824 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4927 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -2.1157503819145E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.1630896156388E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3463815777017E-02
- (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.5814557660633E-05
+ (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.5814557660634E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.2718497826848E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.2718497826835E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7668729590302E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7668729590304E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2450922336620E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.2450922336658E-09
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6211465608393E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6211465608398E-06
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0272142175052E-08
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0272142175063E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9831937878303E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.5303372884598E+00
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.5039114366013E+00
  (PID.TID 0000.0001) %MON dynstat_theta_sd             =   4.7506381250515E+00
- (PID.TID 0000.0001) %MON dynstat_theta_del2           =   2.9893825365562E-03
+ (PID.TID 0000.0001) %MON dynstat_theta_del2           =   2.9893825365557E-03
  (PID.TID 0000.0001) %MON dynstat_salt_max             =   4.0841677349579E+01
  (PID.TID 0000.0001) %MON dynstat_salt_min             =   1.8096651669656E+01
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4756118155337E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   5.1176105627020E-01
- (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4399154994265E-04
+ (PID.TID 0000.0001) %MON dynstat_salt_del2            =   8.4399154994261E-04
- (PID.TID 0000.0001) %MON forcing_qnet_max             =   8.1145496619929E+02
+ (PID.TID 0000.0001) %MON forcing_qnet_max             =   8.1145496619928E+02
  (PID.TID 0000.0001) %MON forcing_qnet_min             =  -2.3176764085007E+02
- (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.6472969211290E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_mean            =  -6.6472969211299E+00
- (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.1870823731159E+02
+ (PID.TID 0000.0001) %MON forcing_qnet_sd              =   1.1870823731158E+02
- (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0846592320291E+00
+ (PID.TID 0000.0001) %MON forcing_qnet_del2            =   1.0846592320292E+00
  (PID.TID 0000.0001) %MON forcing_qsw_max              =  -0.0000000000000E+00
  (PID.TID 0000.0001) %MON forcing_qsw_min              =  -3.3787837809554E+02
  (PID.TID 0000.0001) %MON forcing_qsw_mean             =  -1.8661941074957E+02
-Line 4852 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4955 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_qsw_del2             =   6.4439992217009E-01
  (PID.TID 0000.0001) %MON forcing_empmr_max            =   4.6164105776294E-04
  (PID.TID 0000.0001) %MON forcing_empmr_min            =  -3.1726008424470E-04
- (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -8.4345908496529E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_mean           =  -8.4345908496294E-07
- (PID.TID 0000.0001) %MON forcing_empmr_sd             =   4.8098049809882E-05
+ (PID.TID 0000.0001) %MON forcing_empmr_sd             =   4.8098049809889E-05
- (PID.TID 0000.0001) %MON forcing_empmr_del2           =   8.8858766218931E-07
+ (PID.TID 0000.0001) %MON forcing_empmr_del2           =   8.8858766219010E-07
  (PID.TID 0000.0001) %MON forcing_fu_max               =   2.4586080052368E-01
  (PID.TID 0000.0001) %MON forcing_fu_min               =  -2.4690104382675E-01
  (PID.TID 0000.0001) %MON forcing_fu_mean              =  -3.9278253198990E-03
-Line 4865 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 4968 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON forcing_fv_mean              =  -1.3741594714266E-02
  (PID.TID 0000.0001) %MON forcing_fv_sd                =   7.6550053729987E-02
  (PID.TID 0000.0001) %MON forcing_fv_del2              =   2.5418174123344E-04
- (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.5003708964773E-02
+ (PID.TID 0000.0001) %MON advcfl_uvel_max              =   7.5003708964774E-02
- (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4094575637206E-02
+ (PID.TID 0000.0001) %MON advcfl_vvel_max              =   8.4094575637203E-02
- (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7606725909514E-02
+ (PID.TID 0000.0001) %MON advcfl_wvel_max              =   5.7606725909521E-02
- (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.3639473937606E-02
+ (PID.TID 0000.0001) %MON advcfl_W_hf_max              =   6.3639473937613E-02
  (PID.TID 0000.0001) %MON pe_b_mean                    =  -1.9215925032865E-03
  (PID.TID 0000.0001) %MON ke_max                       =   3.4804936315104E-02
  (PID.TID 0000.0001) %MON ke_mean                      =   1.6153531828135E-04
  (PID.TID 0000.0001) %MON ke_vol                       =   1.3395907263182E+18
  (PID.TID 0000.0001) %MON vort_r_min                   =  -1.2884370648247E-06
- (PID.TID 0000.0001) %MON vort_r_max                   =   1.4155833308234E-06
+ (PID.TID 0000.0001) %MON vort_r_max                   =   1.4155833308232E-06
  (PID.TID 0000.0001) %MON vort_a_mean                  =  -2.0549865324846E-05
  (PID.TID 0000.0001) %MON vort_a_sd                    =   7.5259581905747E-05
  (PID.TID 0000.0001) %MON vort_p_mean                  =  -2.4787191144727E-05
  (PID.TID 0000.0001) %MON vort_p_sd                    =   1.2817300813596E-04
- (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.3816491965627E-08
+ (PID.TID 0000.0001) %MON surfExpan_theta_mean         =   7.3816491965477E-08
- (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9539701124791E-08
+ (PID.TID 0000.0001) %MON surfExpan_salt_mean          =   1.9539701123912E-08
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
-Line 4921 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 5024 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON thSI_Tic2_min_S              =  -8.2917468067431E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_min_N              =  -7.9190492057570E+00
  (PID.TID 0000.0001) %MON thSI_Tic2_max_S              =  -9.9363404317210E-01
- (PID.TID 0000.0001) %MON thSI_Tic2_max_N              =  -1.4028195218075E+00
+ (PID.TID 0000.0001) %MON thSI_Tic2_max_N              =  -1.4028195218074E+00
  (PID.TID 0000.0001) %MON thSI_TotEnerg_G              =  -3.8043484385276E+22
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // End MONITOR Therm.SeaIce statistics
  (PID.TID 0000.0001) // =======================================================
-  cg2d: Sum(rhs),rhsMax =   3.91105770461027E+02  2.30680090912679E+02
+  BULKF_FIELDS_LOAD,     36017 : iP,iLd,i0,i1=    1    2    1    2 ; Wght=  0.9333333333  0.0666666667
- (PID.TID 0000.0001)                    cg2d_init_res =    2.50818125400228E+00
+  EXTERNAL_FIELDS_LOAD,     36017 : iP,iLd,i0,i1=    1    2    1    2 ; Wght=  0.9333333333  0.0666666667
+  cg2d: Sum(rhs),rhsMax =   3.91105770461028E+02  2.30680090912679E+02
+ (PID.TID 0000.0001)                    cg2d_init_res =    2.50818125399965E+00
  (PID.TID 0000.0001)                       cg2d_iters =    69
- (PID.TID 0000.0001)                         cg2d_res =    5.19308520689081E-07
+ (PID.TID 0000.0001)                         cg2d_res =    5.19308520688702E-07
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
  (PID.TID 0000.0001) // =======================================================
  (PID.TID 0000.0001) %MON time_tsnumber                =                 36018
  (PID.TID 0000.0001) %MON time_secondsf                =   3.1119552000000E+09
- (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1281776687796E-01
+ (PID.TID 0000.0001) %MON dynstat_eta_max              =   7.1281776687782E-01
- (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1945177843112E+01
+ (PID.TID 0000.0001) %MON dynstat_eta_min              =  -1.1945177843113E+01
  (PID.TID 0000.0001) %MON dynstat_eta_mean             =  -5.8170844206080E-01
  (PID.TID 0000.0001) %MON dynstat_eta_sd               =   1.6962216825110E+00
- (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1363887962690E-02
+ (PID.TID 0000.0001) %MON dynstat_eta_del2             =   2.1363887962689E-02
  (PID.TID 0000.0001) %MON dynstat_uvel_max             =   1.7265507262978E-01
- (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6392889537648E-01
+ (PID.TID 0000.0001) %MON dynstat_uvel_min             =  -2.6392889537649E-01
  (PID.TID 0000.0001) %MON dynstat_uvel_mean            =  -5.1855701120018E-04
  (PID.TID 0000.0001) %MON dynstat_uvel_sd              =   1.2965291328522E-02
- (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8496442072577E-05
+ (PID.TID 0000.0001) %MON dynstat_uvel_del2            =   5.8496442072576E-05
- (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.9242876818253E-01
+ (PID.TID 0000.0001) %MON dynstat_vvel_max             =   1.9242876818252E-01
  (PID.TID 0000.0001) %MON dynstat_vvel_min             =  -2.1157974079155E-01
- (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.1574615123888E-04
+ (PID.TID 0000.0001) %MON dynstat_vvel_mean            =  -4.1574615123887E-04
  (PID.TID 0000.0001) %MON dynstat_vvel_sd              =   1.3459880644489E-02
  (PID.TID 0000.0001) %MON dynstat_vvel_del2            =   5.5792553791920E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.3126512324643E-05
+ (PID.TID 0000.0001) %MON dynstat_wvel_max             =   9.3126512324644E-05
- (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7643750065875E-04
+ (PID.TID 0000.0001) %MON dynstat_wvel_min             =  -1.7643750065878E-04
- (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.1941350752638E-09
+ (PID.TID 0000.0001) %MON dynstat_wvel_mean            =  -1.1941350752572E-09
- (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6209715638186E-06
+ (PID.TID 0000.0001) %MON dynstat_wvel_sd              =   4.6209715638190E-06
- (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0251749048513E-08
+ (PID.TID 0000.0001) %MON dynstat_wvel_del2            =   4.0251749048524E-08
  (PID.TID 0000.0001) %MON dynstat_theta_max            =   2.9817319827279E+01
  (PID.TID 0000.0001) %MON dynstat_theta_min            =  -3.5483720361727E+00
  (PID.TID 0000.0001) %MON dynstat_theta_mean           =   3.5039864237419E+00
-Line 4964 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
+Line 5069 
 S/R EXTERNAL_FIELDS_LOAD: Reading new da
  (PID.TID 0000.0001) %MON dynstat_salt_min             =   1.8098068269031E+01
  (PID.TID 0000.0001) %MON dynstat_salt_mean            =   3.4756118806044E+01
  (PID.TID 0000.0001) %MON dynstat_salt_sd              =   5.1175276042748E-01
 </