Too small M2 tide amplitude result in estuary

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
cheng_dada
Posts: 10
Joined: Fri Feb 26, 2021 2:20 pm
Location: Shantou University

Too small M2 tide amplitude result in estuary

#1 Unread post by cheng_dada »

Hi.
I encountered too small and abnormal M2 tide amplitude result compared to observation data.
The tidal forcing is from the TPXO 7.2.
This reslut in the weak surface mixing.(or maybe tide is not the only factor affect mixing,but my m2 tide amplitude model result is too sall)
I've tried different 2D momentum boundry condition, the results didn't get any better. And I've also tried several advection and mixing (MY-25 and LMD) schemes, but the results didn't get any better.
Thanks for your tips.
HPT{FONA4Q]I{74}8{}8UDP.png
Fig 1. vertical salt field in estuary, simulated
In fact, the M2 tidal amplitude in the tidal forcing file or the observation data is 0.6
34H(6FR)Q{EO6QW0`C%1PQ9.png
Fig 2.the M2 tidal amplitude in the tidal forcing file
but the result of M2 tidal amplitude of the model is 0.08
-------------------------------------------------------------------------------------------------------------------
Here is my CPP option

Code: Select all

/***************** Model dynamics */
#define SOLVE3D
#define UV_ADV
#define UV_COR
#define UV_QDRAG
#define DJ_GRADPS

/************ Grid configuration */
#define MASKING
#define SPHERICAL
#define CURVGRID


/**** Mixing ********************/
#define TS_DIF2
#define MIX_GEO_TS
#define UV_VIS2
#define MIX_S_UV
/*#define MY25_MIXING*/
#define DIFF_GRID
#define VISC_GRID
#define TS_SMAGORINSKY
#define UV_SMAGORINSKY

#define LMD_MIXING

# ifdef MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# define RI_SPLINES
# endif

# ifdef LMD_MIXING
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_SKPP
# define LMD_BKPP
# define LMD_NONLOCAL
# define LMD_DDMIX
# define RI_SPLINES
# endif

/**** Equation of State *********/
#define NONLIN_EOS
#define SPLINES_VDIFF
#define SPLINES_VVISC
#define SALINITY

/*************** Surface Forcing */
# define SOLAR_SOURCE
# define QCORRECTION 
# define DIURNAL_SRFLUX 

/************** Other options ***/
#define RADIATION_2D
#define AVERAGES
#define ANA_BTFLUX
#define ANA_BSFLUX
#define LIMIT_BSTRESS
/*#define LIMIT_VVISC
#define LIMIT_VDIFF*/

# define WTYPE_GRID
# define ANA_WTYPE


/*#define ANA_SPONGE*/
/*#define ANA_SMFLUX*/
/*#define ANA_SMFLUX*/
/*#define ANA_STFLUX*/
/*#define ANA_SSFLUX*/
/*#define ANA_SRFLUX*/

/********  Tide  ************/

# define RAMP_TIDES
# define SSH_TIDES
# define UV_TIDES
# define ADD_M2OBC
# define ADD_FSOBC

/* nudging the time in *.in shall be set */
#undef  OUT_DOUBLE
#define RST_SINGLE
-------------------------------------------------------------------------------------------------------------------
Here is the log

Code: Select all

--------------------------------------------------------------------------------
 Model Input Parameters:  ROMS/TOMS version 4.1  
                          Tuesday - September 27, 2022 -  3:45:40 AM
--------------------------------------------------------------------------------

 Hanjiang River Plume

 Operating system  : Linux
 CPU/hardware      : x86_64
 Compiler system   : gfortran
 Compiler command  : /usr/local/mpich-3.4.3/bin/mpif90
 Compiler flags    : -frepack-arrays -fallow-argument-mismatch        -O3 -ffast-math -
 OCN Communicator  : 1140850688,  PET size = 72

 Input Script      : upwelling.in

 SVN Root URL      : 
 SVN Revision      : 

 Local Root        : /home/21rqqiao/roms/roms
 Header Dir        : /home/21rqqiao/roms/roms
 Header file       : upwelling.h
 Analytical Dir    : /home/21rqqiao/roms/roms/User/Functionals

 Resolution, Grid 01: 678x403x30,  Parallel Nodes: 72,  Tiling: 8x9


 Physical Parameters, Grid: 01
 =============================

     147808  ntimes            Number of timesteps for 3-D equations.
     75.000  dt                Timestep size (s) for 3-D equations.
         30  ndtfast           Number of timesteps for 2-D equations between
                                 each 3D timestep.
          1  ERstr             Starting ensemble/perturbation run number.
          1  ERend             Ending ensemble/perturbation run number.
          0  nrrec             Number of restart records to read from disk.
          T  LcycleRST         Switch to recycle time-records in restart file.
       1192  nRST              Number of timesteps between the writing of data
                                 into restart fields.
          1  ninfo             Number of timesteps between print of information
                                 to standard output.
          T  ldefout           Switch to create a new output NetCDF file(s).
         48  nHIS              Number of timesteps between the writing fields
                                 into history file.
     100000  ndefHIS           Number of timesteps between creation of new
                                 history files.
          0  nQCK              Number of timesteps between the writing fields
                                 into quicksave file.
          0  ntsAVG            Starting timestep for the accumulation of output
                                 time-averaged data.
       1192  nAVG              Number of timesteps between the writing of
                                 time-averaged data into averages file.
          0  ntsDIA            Starting timestep for the accumulation of output
                                 time-averaged diagnostics data.
          0  nDIA              Number of timesteps between the writing of
                                 time-averaged data into diagnostics file.
 2.0000E+00  nl_tnu2(01)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 01: temp
 2.0000E+00  nl_tnu2(02)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 02: salt
 2.5000E+01  nl_visc2          NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for momentum.
          F  LuvSponge         Turning OFF sponge on horizontal momentum.
          F  LtracerSponge(01) Turning OFF sponge on tracer 01: temp
          F  LtracerSponge(02) Turning OFF sponge on tracer 02: salt
 5.0000E-06  Akt_bak(01)       Background vertical mixing coefficient (m2/s)
                                 for tracer 01: temp
 5.0000E-06  Akt_bak(02)       Background vertical mixing coefficient (m2/s)
                                 for tracer 02: salt
 5.0000E-06  Akv_bak           Background vertical mixing coefficient (m2/s)
                                 for momentum.
 3.0000E-04  rdrg              Linear bottom drag coefficient (m/s).
 2.5000E-03  rdrg2             Quadratic bottom drag coefficient.
 2.0000E-02  Zob               Bottom roughness (m).
          1  lmd_Jwt           Jerlov water type.
          2  Vtransform        S-coordinate transformation equation.
          4  Vstretching       S-coordinate stretching function.
 3.0000E+00  theta_s           S-coordinate surface control parameter.
 4.0000E-01  theta_b           S-coordinate bottom  control parameter.
      3.000  Tcline            S-coordinate surface/bottom layer width (m) used
                                 in vertical coordinate stretching.
   1025.000  rho0              Mean density (kg/m3) for Boussinesq approximation.
  44315.000  dstart            Time-stamp assigned to model initialization (days).
      0.000  tide_start        Reference time origin for tidal forcing (days).
19000101.00  time_ref          Reference time for units attribute (yyyymmdd.dd)
 1.0000E+00  Tnudg(01)         Nudging/relaxation time scale (days)
                                 for tracer 01: temp
 1.0000E+00  Tnudg(02)         Nudging/relaxation time scale (days)
                                 for tracer 02: salt
 3.0000E+01  Znudg             Nudging/relaxation time scale (days)
                                 for free-surface.
 3.0000E+01  M2nudg            Nudging/relaxation time scale (days)
                                 for 2D momentum.
 3.0000E+01  M3nudg            Nudging/relaxation time scale (days)
                                 for 3D momentum.
 1.0000E+01  obcfac            Factor between passive and active
                                 open boundary conditions.
          F  VolCons(1)        NLM western  edge boundary volume conservation.
          F  VolCons(2)        NLM southern edge boundary volume conservation.
          F  VolCons(3)        NLM eastern  edge boundary volume conservation.
          F  VolCons(4)        NLM northern edge boundary volume conservation.
     10.000  T0                Background potential temperature (C) constant.
     34.000  S0                Background salinity (PSU) constant.
      1.000  gamma2            Slipperiness variable: free-slip (1.0) or 
                                                      no-slip (-1.0).
          T  LuvSrc            Turning ON  momentum point Sources/Sinks.
          F  LwSrc             Turning OFF volume influx point Sources/Sinks.
          T  LtracerSrc(01)    Turning ON  point Sources/Sinks on tracer 01: temp
          T  LtracerSrc(02)    Turning ON  point Sources/Sinks on tracer 02: salt
          F  LsshCLM           Turning OFF processing of SSH climatology.
          F  Lm2CLM            Turning OFF processing of 2D momentum climatology.
          T  Lm3CLM            Turning ON  processing of 3D momentum climatology.
          T  LtracerCLM(01)    Turning ON  processing of climatology tracer 01: temp
          T  LtracerCLM(02)    Turning ON  processing of climatology tracer 02: salt
          F  LnudgeM2CLM       Turning OFF nudging of 2D momentum climatology.
          T  LnudgeM3CLM       Turning ON  nudging of 3D momentum climatology.
          T  LnudgeTCLM(01)    Turning ON  nudging of climatology tracer 01: temp
          T  LnudgeTCLM(02)    Turning ON  nudging of climatology tracer 02: salt

          T  Hout(idFsur)      Write out free-surface.
          T  Hout(idUvel)      Write out 3D U-momentum component.
          T  Hout(idVvel)      Write out 3D V-momentum component.
          T  Hout(idTvar)      Write out tracer 01: temp
          T  Hout(idTvar)      Write out tracer 02: salt
          T  Hout(idVvis)      Write out vertical viscosity: AKv.
          T  Hout(idSdif)      Write out vertical diffusion: AKt(isalt).

          T  Aout(idFsur)      Write out averaged free-surface.
          T  Aout(idUvel)      Write out averaged 3D U-momentum component.
          T  Aout(idVvel)      Write out averaged 3D V-momentum component.
          T  Aout(idTvar)      Write out averaged tracer 01: temp
          T  Aout(idTvar)      Write out averaged tracer 02: salt
          T  Aout(idDano)      Write out averaged density anomaly.
          T  Aout(idVvis)      Write out averaged vertical viscosity: AKv.
          T  Aout(idSdif)      Write out averaged vertical diffusion: AKt(isalt).
          T  Aout(idHsbl)      Write out averaged depth of surface boundary layer.
          T  Aout(idHbbl)      Write out averaged depth of bottom boundary layer.

          1  inp_lib           Using standard NetCDF library for input files.
          1  out_lib           Using standard NetCDF library for output files.

  ROMS I/O variables Metadata File:  ROMS/External/varinfo.yaml

  TIDE_DATE - Checking tidal reference date for zero phase: 
                 'zero_phase_date' variable not found in: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_tide.nc
                 given tide_start =        0.0000 (days)
                 tide_DateNumber  =   693962.0000 (1900-01-01 00:00:00)

 Tile partition information for Grid 01:  678x403x30  tiling: 8x9

    
 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     194
               HaloSizeJ(1) =     114
                TileSide(1) =      91
                TileSize(1) =    4641


 Tracer Advection Scheme: NLM
 ========================

 Variable               Grid  Horizontal         Vertical
 ---------              ----  ------------       ------------

 temp                     1   Centered4          Centered4

 salt                     1   Centered4          Centered4

 Akima4     Fourth-order Akima advection
 Centered2  Second-order centered differences advection
 Centered4  Fourth-order centered differences advection
 HSIMT      Third High-order Spatial Inteporlation at Middle Time Advection with TVD limiter
 MPDATA     Multidimensional Positive Definite Advection Algorithm, recursive method
 Splines    Conservative Parabolic Splines Reconstruction Advection (only vertical; not recommended)
 Split_U3   Split third-order Upstream Advection
 Upstream3  Third-order Upstream-biased Advection (only horizontal)


 Lateral Boundary Conditions: NLM
 ============================

 Variable               Grid  West Edge    South Edge   East Edge    North Edge
 ---------              ----  ----------   ----------   ----------   ----------

 zeta                     1   Closed       Chapman Imp  Chapman Imp  Closed

 ubar                     1   Closed       Flather      Flather      Closed

 vbar                     1   Closed       Flather      Flather      Closed

 u                        1   Closed       Rad + Nud    Rad + Nud    Closed

 v                        1   Closed       Rad + Nud    Rad + Nud    Closed

 temp                     1   Closed       Rad + Nud    Rad + Nud    Closed

 salt                     1   Closed       Rad + Nud    Rad + Nud    Closed

 Activated C-preprocessing Options:

 RIVERPLUME               Hanjiang River Plume
 ADD_FSOBC                Adding tidal elevation to processed OBC data
 ADD_M2OBC                Adding tidal currents to processed OBC data
 ANA_BSFLUX               Analytical kinematic bottom salinity flux
 ANA_BTFLUX               Analytical kinematic bottom temperature flux
 ANA_WTYPE                Analytical spatially varying Jerlov water type index
 ASSUMED_SHAPE            Using assumed-shape arrays
 AVERAGES                 Writing out time-averaged nonlinear model fields
 BOUNDARY_ALLREDUCE       Using mpi_allreduce in mp_boundary routine
 COLLECT_ALLREDUCE        Using mpi_allreduce in mp_collect routine
 CURVGRID                 Orthogonal curvilinear grid
 DIAGNOSTICS_UV           Computing and writing momentum diagnostic terms
 DIFF_3DCOEF              Horizontal, time-dependent 3D diffusion coefficient
 DIFF_GRID                Horizontal diffusion coefficient scaled by grid size
 DIURNAL_SRFLUX           Modulate shortwave radiation by the local diurnal cycle
 DJ_GRADPS                Parabolic Splines density Jacobian (Shchepetkin, 2002)
 DOUBLE_PRECISION         Double precision arithmetic numerical kernel.
 LIMIT_BSTRESS            Limit bottom stress to maintain bottom velocity direction
 LMD_BKPP                 KPP bottom boundary layer mixing
 LMD_CONVEC               LMD convective mixing due to shear instability
 LMD_DDMIX                LMD double-diffusive mixing
 LMD_MIXING               Large/McWilliams/Doney interior mixing
 LMD_NONLOCAL             LMD convective nonlocal transport
 LMD_RIMIX                LMD diffusivity due to shear instability
 LMD_SKPP                 KPP surface boundary layer mixing
 MASKING                  Land/Sea masking
 MIX_GEO_TS               Mixing of tracers along geopotential surfaces
 MIX_S_UV                 Mixing of momentum along constant S-surfaces
 MPI                      MPI distributed-memory configuration
 NONLINEAR                Nonlinear Model
 NONLIN_EOS               Nonlinear Equation of State for seawater
 POWER_LAW                Power-law shape time-averaging barotropic filter
 PROFILE                  Time profiling activated
 QCORRECTION              Surface net heat flux correction
 RADIATION_2D             Use tangential phase speed in radiation conditions
 RAMP_TIDES               Ramping tidal forcing for one day
 REDUCE_ALLREDUCE         Using mpi_allreduce in mp_reduce routine
 RI_SPLINES               Parabolic Spline Reconstruction for Richardson Number
 RST_SINGLE               Single precision fields in restart NetCDF file
 SALINITY                 Using salinity
 SOLAR_SOURCE             Solar Radiation Source Term
 SOLVE3D                  Solving 3D Primitive Equations
 SPHERICAL                Spherical grid configuration
 SPLINES_VDIFF            Parabolic Spline Reconstruction for Vertical Diffusion
 SPLINES_VVISC            Parabolic Spline Reconstruction for Vertical Viscosity
 SSH_TIDES                Adding tidal elevation to SSH climatology
 TS_DIF2                  Harmonic mixing of tracers
 TS_SMAGORINSKY           Smagorinksy-like time-dependent diffusion coefficients
 UV_ADV                   Advection of momentum
 UV_COR                   Coriolis term
 UV_U3HADVECTION          Third-order upstream horizontal advection of 3D momentum
 UV_C4VADVECTION          Fourth-order centered vertical advection of momentum
 UV_QDRAG                 Quadratic bottom stress
 UV_TIDES                 Add tidal currents to 2D momentum climatologies
 UV_VIS2                  Harmonic mixing of momentum
 UV_SMAGORINSKY           Smagorinksy-like time-dependent viscosity coefficients
 VAR_RHO_2D               Variable density barotropic mode
 VISC_GRID                Horizontal viscosity coefficient scaled by grid size
 VISC_3DCOEF              Horizontal, time-dependent 3D viscosity coefficient

 CHECKADJ - use caution when activating: TS_SMAGORINSKY
            REASON: stability problems, WARNING.

 CHECKADJ - use caution when activating: UV_SMAGORINSKY
            REASON: stability problems, WARNING.

 Process Information:

 Node #    0 (pid= 1786839) is active.
 Node #    1 (pid= 1786840) is active.
 Node #    2 (pid= 1786841) is active.
 Node #    3 (pid= 1786842) is active.
 Node #    4 (pid= 1786843) is active.
 Node #    5 (pid= 1786844) is active.
 Node #    6 (pid= 1786845) is active.
 Node #    7 (pid= 1786846) is active.
 Node #    8 (pid= 1786847) is active.
 Node #    9 (pid= 1786848) is active.
 Node #   10 (pid= 1786849) is active.
 Node #   11 (pid= 1786850) is active.
 Node #   12 (pid= 1786851) is active.
 Node #   13 (pid= 1786852) is active.
 Node #   14 (pid= 1786853) is active.
 Node #   15 (pid= 1786854) is active.
 Node #   16 (pid= 1786855) is active.
 Node #   17 (pid= 1786856) is active.
 Node #   18 (pid= 1786857) is active.
 Node #   19 (pid= 1786858) is active.
 Node #   20 (pid= 1786859) is active.
 Node #   21 (pid= 1786860) is active.
 Node #   22 (pid= 1786861) is active.
 Node #   23 (pid= 1786862) is active.
 Node #   24 (pid= 1786863) is active.
 Node #   25 (pid= 1786864) is active.
 Node #   26 (pid= 1786865) is active.
 Node #   27 (pid= 1786866) is active.
 Node #   28 (pid= 1786867) is active.
 Node #   29 (pid= 1786868) is active.
 Node #   30 (pid= 1786869) is active.
 Node #   31 (pid= 1786870) is active.
 Node #   32 (pid= 1786871) is active.
 Node #   33 (pid= 1786872) is active.
 Node #   34 (pid= 1786873) is active.
 Node #   35 (pid= 1786874) is active.
 Node #   36 (pid= 1786875) is active.
 Node #   37 (pid= 1786876) is active.
 Node #   38 (pid= 1786877) is active.
 Node #   39 (pid= 1786878) is active.
 Node #   40 (pid= 1786879) is active.
 Node #   41 (pid= 1786880) is active.
 Node #   42 (pid= 1786881) is active.
 Node #   43 (pid= 1786882) is active.
 Node #   44 (pid= 1786883) is active.
 Node #   45 (pid= 1786884) is active.
 Node #   46 (pid= 1786885) is active.
 Node #   47 (pid= 1786886) is active.
 Node #   48 (pid= 1786887) is active.
 Node #   49 (pid= 1786888) is active.
 Node #   50 (pid= 1786889) is active.
 Node #   51 (pid= 1786890) is active.
 Node #   53 (pid= 1786892) is active.
 Node #   54 (pid= 1786893) is active.
 Node #   55 (pid= 1786894) is active.
 Node #   56 (pid= 1786895) is active.
 Node #   57 (pid= 1786896) is active.
 Node #   58 (pid= 1786897) is active.
 Node #   59 (pid= 1786898) is active.
 Node #   60 (pid= 1786899) is active.
 Node #   61 (pid= 1786900) is active.
 Node #   62 (pid= 1786901) is active.
 Node #   63 (pid= 1786902) is active.
 Node #   64 (pid= 1786903) is active.
 Node #   65 (pid= 1786904) is active.
 Node #   66 (pid= 1786905) is active.
 Node #   67 (pid= 1786906) is active.
 Node #   68 (pid= 1786907) is active.
 Node #   69 (pid= 1786908) is active.
 Node #   70 (pid= 1786909) is active.
 Node #   71 (pid= 1786910) is active.
 Node #   52 (pid= 1786891) is active.

 Point Sources/Sinks grid-cell flag locations counter:

    Grid: 1, file: /home/21rqqiao/roms/roms/Data/ROMS/Initial/Hanjiang.nc
                   LuvSrc = T  u-face = 9
                   LuvSrc = T  v-face = 0
                   LwSrc  = F  w-face = 0

 INITIAL: Configuring and initializing forward nonlinear model ...
 *******

  GET_GRID_NF90    - angle between XI-axis and EAST: angler
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min = -1.79942962E-01 Max =  1.56305444E+00)
  GET_GRID_NF90    - bathymetry at RHO-points: h
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  3.00000000E+00 Max =  3.00000000E+03)
  GET_GRID_NF90    - Coriolis parameter at RHO-points: f
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  4.20248309E-05 Max =  8.10971003E-05)
  GET_GRID_NF90    - reciprocal XI-grid spacing: pm
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  2.55322069E-05 Max =  2.32416345E-03)
  GET_GRID_NF90    - reciprocal ETA-grid spacing: pn
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.27345932E-04 Max =  2.65526570E-03)
  GET_GRID_NF90    - XI-derivative of inverse metric factor pn: dndx
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min = -8.99112406E+01 Max =  1.00381802E+02)
  GET_GRID_NF90    - ETA-derivative of inverse metric factor pm: dmde
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min = -1.71482438E+02 Max =  1.49350879E-01)
  GET_GRID_NF90    - x-location of RHO-points: x-rho
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  0.00000000E+00 Max =  3.37136568E+06)
  GET_GRID_NF90    - x-location of U-points: x_u
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.39926436E+04 Max =  3.35458238E+06)
  GET_GRID_NF90    - x-location of V-points: x_v
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  0.00000000E+00 Max =  3.36829988E+06)
  GET_GRID_NF90    - x-location of PSI-points: x_psi
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.40338257E+04 Max =  3.35154083E+06)
  GET_GRID_NF90    - y-location of RHO-points: y_rho
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  0.00000000E+00 Max =  5.90092190E+05)
  GET_GRID_NF90    - y-location of U-points: y_u
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  0.00000000E+00 Max =  5.85734762E+05)
  GET_GRID_NF90    - y-location of V-points: y_v
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.25213111E+03 Max =  5.86198109E+05)
  GET_GRID_NF90    - y-location of PSI-points: y-psi
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.25226748E+03 Max =  5.81863066E+05)
  GET_GRID_NF90    - longitude of RHO-points: lon_rho
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.03375440E+02 Max =  1.25442520E+02)
  GET_GRID_NF90    - longitude of U-points: lon_u
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.03560510E+02 Max =  1.25440540E+02)
  GET_GRID_NF90    - longitude of V-points: lon_v
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.03379040E+02 Max =  1.25425970E+02)
  GET_GRID_NF90    - longitude of PSI-points: lon_psi
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.03563955E+02 Max =  1.25423987E+02)
  GET_GRID_NF90    - latitude of RHO-points lat_rho
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.67945930E+01 Max =  3.38887820E+01)
  GET_GRID_NF90    - latitude of U-points: lat_u
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.67949770E+01 Max =  3.37776815E+01)
  GET_GRID_NF90    - latitude of V-points: lat_v
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.68115125E+01 Max =  3.38883325E+01)
  GET_GRID_NF90    - latitude of PSI-points lat_psi
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  1.68118850E+01 Max =  3.37772025E+01)
  GET_GRID_NF90    - mask on RHO-points: mask_rho
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  0.00000000E+00 Max =  1.00000000E+00)
  GET_GRID_NF90    - mask on U-points: mask_u
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  0.00000000E+00 Max =  1.00000000E+00)
  GET_GRID_NF90    - mask on V-points: mask_v
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  0.00000000E+00 Max =  1.00000000E+00)
  GET_GRID_NF90    - mask on PSI-points: mask_psi
                      (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/ST_grid.nc)
                      (Min =  0.00000000E+00 Max =  1.00000000E+00)

 Vertical S-coordinate System, Grid 01:

 level   S-coord     Cs-curve   Z   at hmin       at hc    half way     at hmax

    30   0.0000000   0.0000000        0.000       0.000       0.000       0.000
    29  -0.0333333  -0.0006695       -0.051      -0.051      -1.103      -2.106
    28  -0.0666667  -0.0026838       -0.104      -0.104      -4.221      -8.243
    27  -0.1000000  -0.0060605       -0.159      -0.159      -9.381     -18.463
    26  -0.1333333  -0.0108285       -0.216      -0.216     -16.626     -32.853
    25  -0.1666667  -0.0170290       -0.276      -0.276     -26.017     -51.535
    24  -0.2000000  -0.0247149       -0.337      -0.337     -37.634     -74.670
    23  -0.2333333  -0.0339514       -0.401      -0.401     -51.575    -102.452
    22  -0.2666667  -0.0448161       -0.467      -0.467     -67.956    -135.113
    21  -0.3000000  -0.0573991       -0.536      -0.536     -86.911    -172.924
    20  -0.3333333  -0.0718033       -0.608      -0.608    -108.596    -216.194
    19  -0.3666667  -0.0881450       -0.682      -0.682    -133.184    -265.270
    18  -0.4000000  -0.1065532       -0.760      -0.760    -160.868    -320.539
    17  -0.4333333  -0.1271703       -0.841      -0.841    -191.863    -382.429
    16  -0.4666667  -0.1501521       -0.925      -0.925    -226.401    -451.405
    15  -0.5000000  -0.1756670       -1.014      -1.014    -264.735    -527.973
    14  -0.5333333  -0.2038961       -1.106      -1.106    -307.136    -612.676
    13  -0.5666667  -0.2350323       -1.203      -1.203    -353.894    -706.091
    12  -0.6000000  -0.2692794       -1.304      -1.304    -405.313    -808.829
    11  -0.6333333  -0.3068504       -1.410      -1.410    -461.713    -921.530
    10  -0.6666667  -0.3479655       -1.522      -1.522    -523.424   -1044.852
     9  -0.7000000  -0.3928500       -1.639      -1.639    -590.784   -1179.470
     8  -0.7333333  -0.4417302       -1.763      -1.763    -664.131   -1326.065
     7  -0.7666667  -0.4948299       -1.892      -1.892    -743.801   -1485.305
     6  -0.8000000  -0.5523651       -2.029      -2.029    -830.118   -1657.837
     5  -0.8333333  -0.6145374       -2.172      -2.172    -923.383   -1844.268
     4  -0.8666667  -0.6815273       -2.322      -2.322   -1023.868   -2045.137
     3  -0.9000000  -0.7534850       -2.480      -2.480   -1131.796   -2260.894
     2  -0.9333333  -0.8305209       -2.646      -2.646   -1247.335   -2491.871
     1  -0.9666667  -0.9126940       -2.819      -2.819   -1370.572   -2738.244
     0  -1.0000000  -1.0000000       -3.000      -3.000   -1501.500   -3000.000

 Time Splitting Weights for Grid 01:    ndtfast =  30    nfast =  42
 ==================================

    Primary            Secondary            Accumulated to Current Step

  1-0.0008094437383769 0.0333333333333333-0.0008094437383769 0.0333333333333333
  2-0.0014053566728197 0.0333603147912792-0.0022148004111966 0.0666936481246126
  3-0.0017877524645903 0.0334071600137066-0.0040025528757869 0.1001008081383191
  4-0.0019566842408176 0.0334667517625262-0.0059592371166046 0.1335675599008453
  5-0.0019122901320372 0.0335319745705535-0.0078715272486418 0.1670995344713988
  6-0.0016548570247459 0.0335957175749547-0.0095263842733877 0.2006952520463536
  7-0.0011849025289723 0.0336508794757796-0.0107112868023600 0.2343461315221331
  8-0.0005032751608631 0.0336903762267453-0.0112145619632232 0.2680365077488784
  9 0.0003887272597151 0.0337071520654408-0.0108258347035081 0.3017436598143192
 10 0.0014892209965583 0.0336941944901169-0.0093366137069498 0.3354378543044362
 11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
 12 0.0043042707117221 0.0335513677379153-0.0022367614257357 0.4026337758325831
 13 0.0060106451121704 0.0334078920475245 0.0037738836864347 0.4360416678801076
 14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905598
 15 0.0099910761708920 0.0329439123123590 0.0216737068001212 0.5021931174029188
 16 0.0122483446563884 0.0326108764399960 0.0339220514565096 0.5348039938429148
 17 0.0146692120341107 0.0322025982847830 0.0485912634906203 0.5670065921276978
 18 0.0172400033810439 0.0317136245503127 0.0658312668716642 0.5987202166780105
 19 0.0199444086685725 0.0311389577709445 0.0857756755402367 0.6298591744489550
 20 0.0227631639997064 0.0304741441486588 0.1085388395399431 0.6603333185976138
 21 0.0256737146312910 0.0297153720153352 0.1342125541712341 0.6900486906129490
 22 0.0286498597812016 0.0288595815276255 0.1628624139524357 0.7189082721405746
 23 0.0316613792205220 0.0279045862015855 0.1945237931729577 0.7468128583421600
 24 0.0346736416507075 0.0268492068942347 0.2291974348236651 0.7736620652363948
 25 0.0376471948657328 0.0256934188392112 0.2668446296893979 0.7993554840756060
 26 0.0405373376992232 0.0244385123436867 0.3073819673886211 0.8237939964192927
 27 0.0432936737565711 0.0230872677537126 0.3506756411451922 0.8468812641730054
 28 0.0458596469320356 0.0216441452951603 0.3965352880772278 0.8685254094681656
 29 0.0481720587108284 0.0201154903974257 0.4447073467880562 0.8886408998655914
 30 0.0501605672561820 0.0185097551070648 0.4948679140442381 0.9071506549726561
 31 0.0517471682814030 0.0168377361985254 0.5466150823256412 0.9239883911711815
 32 0.0528456577069106 0.0151128305891453 0.5994607400325517 0.9391012217603267
 33 0.0533610761022577 0.0133513086655816 0.6528218161348094 0.9524525304259084
 34 0.0531891349131379 0.0115726061288397 0.7060109510479473 0.9640251365547481
 35 0.0522156244733761 0.0097996349650684 0.7582265755213234 0.9738247715198165
 36 0.0503158038019031 0.0080591141492892 0.8085423793232265 0.9818838856691057
 37 0.0473537721847153 0.0063819206892258 0.8558961515079418 0.9882658063583315
 38 0.0431818225418188 0.0048034616164019 0.8990779740497606 0.9930692679747335
 39 0.0376397765791564 0.0033640675316746 0.9367177506289170 0.9964333355064081
 40 0.0305543017255206 0.0021094083123694 0.9672720523544377 0.9985427438187775
 41 0.0217382098544505 0.0010909315881854 0.9890102622088881 0.9996336754069629
 42 0.0109897377911119 0.0003663245930371 1.0000000000000000 1.0000000000000000

 ndtfast, nfast =   30  42   nfast/ndtfast =  1.40000

 Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):

    1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.18933

 Metrics information for Grid 01:
 ===============================

 Minimum X-grid spacing, DXmin =  4.30262339E-01 km    Water points =  4.68777884E-01 km
 Maximum X-grid spacing, DXmax =  3.91662188E+01 km    Water points =  3.35665985E+01 km
 Minimum Y-grid spacing, DYmin =  3.76610144E-01 km    Water points =  3.76610144E-01 km
 Maximum Y-grid spacing, DYmax =  7.85262616E+00 km    Water points =  7.72786401E+00 km
 Minimum Z-grid spacing, DZmin =  5.10042627E-02 m     Water points =  5.10042627E-02 m
 Maximum Z-grid spacing, DZmax =  2.61756270E+02 m     Water points =  2.61756270E+02 m

 Minimum barotropic Courant Number =  3.42602393E-03
 Maximum barotropic Courant Number =  7.35849651E-01
 Maximum Coriolis   Courant Number =  6.08228252E-03

 Horizontal mixing scaled by grid area squared root, MAXVAL(grdscl) =  1.25607878E+01 km  (Water points)

 Minimum horizontal diffusion coefficient =  3.91374421E+00 m2/s
 Maximum horizontal diffusion coefficient =  1.23388747E+02 m2/s

 Minimum horizontal viscosity coefficient =  3.91374421E+00 m2/s
 Maximum horizontal viscosity coefficient =  1.23388747E+02 m2/s

   GET_NUDGCOEF_NF90 - 3D momentum inverse nudging coefficients: M3_NudgeCoef
                     (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/nudgcoef_b.nc)
                     (Min =  0.00000000E+00 Max =  4.62962963E-07)
   GET_NUDGCOEF_NF90 - potential temperature, inverse nudging coefficients: temp_NudgeCoef
                     (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/nudgcoef_b.nc)
                     (Min =  0.00000000E+00 Max =  4.62962963E-07)
   GET_NUDGCOEF_NF90 - salinity, inverse nudging coefficients: salt_NudgeCoef
                     (Grid = 01, File: /home/21rqqiao/roms/roms/Data/ROMS/Initial/nudgcoef_b.nc)
                     (Min =  0.00000000E+00 Max =  4.62962963E-07)

  GET_STATE_NF90   - NLM: state initial conditions,                       2021-05-01 00:00:00.00
                      (Grid 01, t = 44315.0000, File: Yuedong_ini_ECCO_Y2021M5.nc, Rec=0001, Index=1)
                   - free-surface
                      (Min =  2.50249269E-01 Max =  1.12928637E+00)
                   - vertically integrated u-momentum component
                      (Min = -6.30954225E-02 Max =  4.86638912E-02)
                   - vertically integrated v-momentum component
                      (Min = -4.20696065E-02 Max =  7.15219872E-02)
                   - u-momentum component
                      (Min = -3.57840844E-01 Max =  5.21045566E-01)
                   - v-momentum component
                      (Min = -8.55738064E-01 Max =  1.08812017E+00)
                   - potential temperature
                      (Min =  2.17789190E+00 Max =  2.89318040E+01)
                   - salinity
                      (Min =  3.02262479E+01 Max =  3.46201475E+01)
  GET_NGFLD_NF90   - tidal period
                      (Grid = 01, Min =  4.30819210E+04 Max =  2.38071599E+06)
  GET_2DFLD_NF90   - tidal elevation amplitude
                      (Grid = 01, Min =  5.56285643E-05 Max =  1.82574356E+01)
  GET_2DFLD_NF90   - tidal elevation phase angle
                      (Grid = 01, Min =  4.96871998E-06 Max =  6.28317701E+00)
  GET_2DFLD_NF90   - tidal current inclination angle
                      (Grid = 01, Min =  7.31209536E-07 Max =  6.28318107E+00)
  GET_2DFLD_NF90   - tidal current phase angle
                      (Grid = 01, Min =  5.76077037E-06 Max =  6.28318179E+00)
  GET_2DFLD_NF90   - maximum tidal current, ellipse major axis
                      (Grid = 01, Min =  2.63947613E-06 Max =  1.96808252E+00)
  GET_2DFLD_NF90   - minimum tidal current, ellipse minor axis
                      (Grid = 01, Min = -6.30073752E-01 Max =  5.05945978E-01)
  GET_NGFLD_NF90   - river runoff  XI-positions at RHO-points
                      (Grid = 01, Min =  5.30000000E+01 Max =  6.70000000E+02)
  GET_NGFLD_NF90   - river runoff ETA-positions at RHO-points
                      (Grid = 01, Min =  3.77000000E+02 Max =  4.04000000E+02)
  GET_NGFLD_NF90   - river runoff direction
                      (Grid = 01, Min =  0.00000000E+00 Max =  0.00000000E+00)
  GET_NGFLD_NF90   - river runoff mass transport vertical profile
                      (Grid = 01, Min =  3.33333333E-02 Max =  3.33333333E-02)
  GET_NGFLD_NF90   - river runoff mass transport,                         2021-02-19 00:00:00.00
                      (Grid= 01, Rec=50, Index=2, File: Hanjiang.nc)
                      (Tmin=      44195.0000 Tmax=      44638.0000)   t =      44244.0000
                      (Min =  1.24666670E+02 Max =  1.12000000E+04)
  GET_NGFLD_NF90   - river runoff potential temperature,                  2021-02-19 00:00:00.00
                      (Grid= 01, Rec=50, Index=2, File: Hanjiang.nc)
                      (Tmin=      44195.0000 Tmax=      44638.0000)   t =      44244.0000
                      (Min =  1.65484444E+01 Max =  1.65484444E+01)
  GET_NGFLD_NF90   - river runoff salinity,                               2021-02-19 00:00:00.00
                      (Grid= 01, Rec=50, Index=2, File: Hanjiang.nc)
                      (Tmin=      44195.0000 Tmax=      44638.0000)   t =      44244.0000
                      (Min =  0.00000000E+00 Max =  0.00000000E+00)
  GET_2DFLD_NF90   - surface u-momentum stress,                           2021-04-30 18:00:00.00
                      (Grid=01, Rec=4, Index=1, File: ST_frc_new1.nc)
                      (Tmin=      44314.0000 Tmax=      44376.7500)   t =      44314.7500
                      (Min = -1.26594052E-04 Max =  2.90227340E-05)   regrid = F
  GET_2DFLD_NF90   - surface v-momentum stress,                           2021-04-30 18:00:00.00
                      (Grid=01, Rec=4, Index=1, File: ST_frc_new1.nc)
                      (Tmin=      44314.0000 Tmax=      44376.7500)   t =      44314.7500
                      (Min = -1.95739635E-04 Max =  7.85255581E-05)   regrid = F
  GET_2DFLD_NF90   - net solar shortwave radiation flux,                  1900-04-17 12:45:04.00
                      (Grid=01, Rec=4, Index=1, File: ST_frc_new1.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        106.5313
                      (Min =  4.72144785E-05 Max =  6.62287306E-05)   regrid = F
  GET_2DFLD_NF90   - surface net heat flux,                               1900-04-17 12:45:04.00
                      (Grid=01, Rec=4, Index=1, File: ST_frc_new1.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        106.5313
                      (Min = -1.43536050E-05 Max =  3.33501383E-05)   regrid = F
  GET_2DFLD_NF90   - sea surface temperature climatology,                 1900-04-17 12:45:04.00
                      (Grid=01, Rec=4, Index=1, File: ST_frc_new1.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        106.5313
                      (Min =  9.10685860E+00 Max =  2.74500008E+01)   regrid = F
  GET_2DFLD_NF90   - surface net heat flux sensitivity to SST,            1900-04-17 12:45:04.00
                      (Grid=01, Rec=4, Index=1, File: ST_frc_new1.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        106.5313
                      (Min = -1.05751368E-05 Max = -5.28289999E-06)   regrid = F
  GET_2DFLD_NF90   - data surface net freshwater flux, (E-P)/rhow,        1900-04-17 12:45:04.00
                      (Grid=01, Rec=4, Index=1, File: ST_frc_new1.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        106.5313
                      (Min = -8.99064706E-08 Max =  3.35863642E-08)   regrid = F
  GET_NGFLD_NF90   - free-surface eastern boundary condition,             2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min =  2.50249269E-01 Max =  3.33052880E-01)
  GET_NGFLD_NF90   - free-surface southern boundary condition,            2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min =  3.27961978E-01 Max =  1.12350350E+00)
  GET_NGFLD_NF90   - 2D u-momentum eastern boundary condition,            2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min = -4.28626576E-04 Max =  5.62721167E-04)
  GET_NGFLD_NF90   - 2D v-momentum eastern boundary condition,            2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min = -2.77192669E-04 Max =  3.14600176E-04)
  GET_NGFLD_NF90   - 2D u-momentum southern boundary condition,           2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min = -5.32939421E-02 Max =  2.56420099E-02)
  GET_NGFLD_NF90   - 2D v-momentum southern boundary condition,           2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min = -3.99476294E-02 Max =  7.15219872E-02)
  GET_NGFLD_NF90   - 3D u-momentum eastern boundary condition,            2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min = -4.78231996E-02 Max =  3.79219202E-02)
  GET_NGFLD_NF90   - 3D v-momentum eastern boundary condition,            2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min = -3.30947583E-02 Max =  4.26176135E-02)
  GET_NGFLD_NF90   - 3D u-momentum southern boundary condition,           2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min = -2.71985247E-01 Max =  5.21045566E-01)
  GET_NGFLD_NF90   - 3D v-momentum southern boundary condition,           2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min = -8.55738064E-01 Max =  1.08812017E+00)
  GET_NGFLD_NF90   - temperature eastern boundary condition,              2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min =  9.48885098E+00 Max =  1.54457965E+01)
  GET_NGFLD_NF90   - salinity eastern boundary condition,                 2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min =  3.10616224E+01 Max =  3.35300721E+01)
  GET_NGFLD_NF90   - temperature southern boundary condition,             2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min =  2.17789190E+00 Max =  2.89318040E+01)
  GET_NGFLD_NF90   - salinity southern boundary condition,                2021-04-30 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44314.0000
                      (Min =  3.16789866E+01 Max =  3.46201475E+01)
  GET_3DFLD_NF90   - u-momentum component climatology,                    1900-04-17 12:45:04.00
                      (Grid=01, Rec=4, Index=1, File: ST_clm.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        106.5313
                      (Min = -4.03473266E-01 Max =  6.87254525E-01)
  GET_3DFLD_NF90   - v-momentum component climatology,                    1900-04-17 12:45:04.00
                      (Grid=01, Rec=4, Index=1, File: ST_clm.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        106.5313
                      (Min = -5.23833438E-01 Max =  5.24011684E-01)
  GET_3DFLD_NF90   - potential temperature climatology,                   1900-04-17 12:45:04.00
                      (Grid=01, Rec=4, Index=1, File: ST_clm.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        106.5313
                      (Min =  0.00000000E+00 Max =  2.92570029E+01)
  GET_3DFLD_NF90   - salinity climatology,                                1900-04-17 12:45:04.00
                      (Grid=01, Rec=4, Index=1, File: ST_clm.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        106.5313
                      (Min =  3.11222754E+01 Max =  3.50000000E+01)

 Basin information for Grid 01:

 Maximum grid stiffness ratios:  rx0 =   3.000000E-01 (Beckmann and Haidvogel)
                                 rx1 =   8.161356E+00 (Haney)

 Initial domain volumes:  TotVolume =  3.0614181190E+14 m3
                         MinCellVol =  1.2314025737E+04 m3
                         MaxCellVol =  1.2757833010E+10 m3
                            Max/Min =  1.0360407947E+06
 NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000000000001 - 000000147808)

  GET_NGFLD_NF90   - river runoff mass transport,                         2022-03-20 00:00:00.00
                      (Grid= 01, Rec=51, Index=1, File: Hanjiang.nc)
                      (Tmin=      44195.0000 Tmax=      44638.0000)   t =      44638.0000
                      (Min =  1.22666667E+02 Max =  1.12000000E+04)
  GET_NGFLD_NF90   - river runoff potential temperature,                  2022-03-20 00:00:00.00
                      (Grid= 01, Rec=51, Index=1, File: Hanjiang.nc)
                      (Tmin=      44195.0000 Tmax=      44638.0000)   t =      44638.0000
                      (Min =  1.69004134E+01 Max =  1.69004134E+01)
  GET_NGFLD_NF90   - river runoff salinity,                               2022-03-20 00:00:00.00
                      (Grid= 01, Rec=51, Index=1, File: Hanjiang.nc)
                      (Tmin=      44195.0000 Tmax=      44638.0000)   t =      44638.0000
                      (Min =  0.00000000E+00 Max =  0.00000000E+00)
  GET_2DFLD_NF90   - surface u-momentum stress,                           2021-05-01 00:00:00.00
                      (Grid=01, Rec=5, Index=2, File: ST_frc_new1.nc)
                      (Tmin=      44314.0000 Tmax=      44376.7500)   t =      44315.0000
                      (Min = -1.78876892E-04 Max =  9.21818476E-05)   regrid = F
  GET_2DFLD_NF90   - surface v-momentum stress,                           2021-05-01 00:00:00.00
                      (Grid=01, Rec=5, Index=2, File: ST_frc_new1.nc)
                      (Tmin=      44314.0000 Tmax=      44376.7500)   t =      44315.0000
                      (Min = -1.38473064E-04 Max =  2.80410051E-05)   regrid = F
  GET_2DFLD_NF90   - net solar shortwave radiation flux,                  1900-05-17 23:15:04.00
                      (Grid=01, Rec=5, Index=2, File: ST_frc_new1.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        136.9688
                      (Min =  5.30203717E-05 Max =  6.88657458E-05)   regrid = F
  GET_2DFLD_NF90   - surface net heat flux,                               1900-05-17 23:15:04.00
                      (Grid=01, Rec=5, Index=2, File: ST_frc_new1.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        136.9688
                      (Min =  9.82776273E-06 Max =  3.86719500E-05)   regrid = F
  GET_2DFLD_NF90   - sea surface temperature climatology,                 1900-05-17 23:15:04.00
                      (Grid=01, Rec=5, Index=2, File: ST_frc_new1.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        136.9688
                      (Min =  1.35000000E+01 Max =  2.92500000E+01)   regrid = F
  GET_2DFLD_NF90   - surface net heat flux sensitivity to SST,            1900-05-17 23:15:04.00
                      (Grid=01, Rec=5, Index=2, File: ST_frc_new1.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        136.9688
                      (Min = -1.02598503E-05 Max = -5.54481490E-06)   regrid = F
  GET_2DFLD_NF90   - data surface net freshwater flux, (E-P)/rhow,        1900-05-17 23:15:04.00
                      (Grid=01, Rec=5, Index=2, File: ST_frc_new1.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        136.9688
                      (Min = -6.94307630E-08 Max =  2.44499032E-08)   regrid = F
  GET_NGFLD_NF90   - free-surface eastern boundary condition,             2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min =  2.50249269E-01 Max =  3.33052880E-01)
  GET_NGFLD_NF90   - free-surface southern boundary condition,            2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min =  3.27961978E-01 Max =  1.12350350E+00)
  GET_NGFLD_NF90   - 2D u-momentum eastern boundary condition,            2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min = -4.28626576E-04 Max =  5.62721167E-04)
  GET_NGFLD_NF90   - 2D v-momentum eastern boundary condition,            2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min = -2.77192669E-04 Max =  3.14600176E-04)
  GET_NGFLD_NF90   - 2D u-momentum southern boundary condition,           2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min = -5.32939421E-02 Max =  2.56420099E-02)
  GET_NGFLD_NF90   - 2D v-momentum southern boundary condition,           2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min = -3.99476294E-02 Max =  7.15219872E-02)
  GET_NGFLD_NF90   - 3D u-momentum eastern boundary condition,            2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min = -4.78231996E-02 Max =  3.79219202E-02)
  GET_NGFLD_NF90   - 3D v-momentum eastern boundary condition,            2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min = -3.30947583E-02 Max =  4.26176135E-02)
  GET_NGFLD_NF90   - 3D u-momentum southern boundary condition,           2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min = -2.71985247E-01 Max =  5.21045566E-01)
  GET_NGFLD_NF90   - 3D v-momentum southern boundary condition,           2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min = -8.55738064E-01 Max =  1.08812017E+00)
  GET_NGFLD_NF90   - temperature eastern boundary condition,              2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min =  9.48885098E+00 Max =  1.54457965E+01)
  GET_NGFLD_NF90   - salinity eastern boundary condition,                 2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min =  3.10616224E+01 Max =  3.35300721E+01)
  GET_NGFLD_NF90   - temperature southern boundary condition,             2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min =  2.17789190E+00 Max =  2.89318040E+01)
  GET_NGFLD_NF90   - salinity southern boundary condition,                2021-05-03 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
                      (Tmin=      44314.0000 Tmax=      44347.0000)   t =      44317.0000
                      (Min =  3.16789866E+01 Max =  3.46201475E+01)
  GET_3DFLD_NF90   - u-momentum component climatology,                    1900-05-17 23:15:04.00
                      (Grid=01, Rec=5, Index=2, File: ST_clm.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        136.9688
                      (Min = -2.87496338E-01 Max =  8.53123859E-01)
  GET_3DFLD_NF90   - v-momentum component climatology,                    1900-05-17 23:15:04.00
                      (Grid=01, Rec=5, Index=2, File: ST_clm.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        136.9688
                      (Min = -4.87102904E-01 Max =  4.78485647E-01)
  GET_3DFLD_NF90   - potential temperature climatology,                   1900-05-17 23:15:04.00
                      (Grid=01, Rec=5, Index=2, File: ST_clm.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        136.9688
                      (Min =  0.00000000E+00 Max =  3.10250610E+01)
  GET_3DFLD_NF90   - salinity climatology,                                1900-05-17 23:15:04.00
                      (Grid=01, Rec=5, Index=2, File: ST_clm.nc)
                      (Tmin=         15.2188 Tmax=        350.0313)   t =        136.9688
                      (Min =  2.97641872E+01 Max =  3.50000000E+01)

Post Reply