BLOWUP. REASON: MaxSpeed = 2.1E+01

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
Zeng
Posts: 26
Joined: Mon Sep 19, 2022 1:06 pm
Location: Sun Yat-Sen University

BLOWUP. REASON: MaxSpeed = 2.1E+01

#1 Unread post by Zeng »

Hello everyone! I am using ROMS to simulate the diffusion of tracers over several days. But I encountered the BLOW UP REASON: MaxSpeed = 2.1E+01 at early steps(71).

The maximum water depth of research area is 142m, the resolution of grid is 0.3km. I checked the rst.nc,I found the max speed in V(-22m/s), and the maxspeed occurs at the sea surface, close to the boundary and near the shore.(In fact, the maxspeed of both the seabed and the sea surface exceeds 20, but it is the highest on the surface.)
微信截图_20230627200259.png
I have tried to find some solutions from forum. Such lke reduce DT to 10 or 5, define WET_DRY in the .h file, set the minimun depth to 2m. But all these are unuseful.The Courant number is small enough I think.

Code: Select all

Minimum barotropic Courant Number =  7.07604957E-03
 Maximum barotropic Courant Number =  5.34537607E-02
I created grid with Gridbuilder, and the values of rx1 and rx0 also meet the requirements.

Code: Select all

 Maximum grid stiffness ratios:  rx0 =   2.000000E-01 (Beckmann and Haidvogel)
                                 rx1 =   5.672658E+00 (Haney)
I have also export DEBUG=ON in build_rom.sh, unfortunately, I don't know how to degug it base on my poor ROMS konwledge.. The last solution I will take later is making the smoother bathymetry. I am confused, because I have already smooth bathymetry with Gridbuilder, need I use another methods to smooth bathymetry? if so, what kind of smooth tool should I use? Thank you for you attention, I would like to express my sincere gratitude for any kindly reply.
Last edited by Zeng on Wed Jun 28, 2023 1:20 am, edited 2 times in total.

Zeng
Posts: 26
Joined: Mon Sep 19, 2022 1:06 pm
Location: Sun Yat-Sen University

Re: BLOWUP. REASON: MaxSpeed = 2.1E+01

#2 Unread post by Zeng »

The log as shown below:

Code: Select all

(base) enzo@enzo:~/roms/project/yangjiang_0.3km$ ./romsO < roms_riverplume1_uv.in 
--------------------------------------------------------------------------------
 Model Input Parameters:  ROMS/TOMS version 4.1  
                          Tuesday - June 27, 2023 -  9:28:47 PM
--------------------------------------------------------------------------------

 River Plume Test 1, flow across u-face grid cell

 Operating system  : Linux
 CPU/hardware      : x86_64
 Compiler system   : ifort
 Compiler command  : /opt/intel/oneapi/compiler/2022.0.1/linux/bin/intel64/ifort
 Compiler flags    : -fp-model precise -heap-arrays -ip -O3 -traceback -check uninit -qopenmp -fpp 
 SVN Root URL      : https://www.myroms.org/svn/src/trunk
 SVN Revision      : 1174

 Local Root        : /home/enzo/roms/trunk
 Header Dir        : /home/enzo/roms/project/yangjiang_0.3km
 Header file       : riverplume1.h
 Analytical Dir    : /home/enzo/roms/project/yangjiang_0.3km

 Resolution, Grid 01: 1098x748x15,  Parallel Threads: 16,  Tiling: 4x4


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

      51840  ntimes            Number of timesteps for 3-D equations.
     10.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.
        720  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).
        720  nHIS              Number of timesteps between the writing fields
                                 into history file.
          0  nQCK              Number of timesteps between the writing fields
                                 into quicksave file.
          1  ntsAVG            Starting timestep for the accumulation of output
                                 time-averaged data.
        720  nAVG              Number of timesteps between the writing of
                                 time-averaged data into averages file.
 5.0000E+00  nl_tnu2(01)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 01: temp
 5.0000E+00  nl_tnu2(02)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 02: salt
 5.0000E+00  nl_tnu2(03)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 03: dye_01
 1.5000E+00  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
          F  LtracerSponge(03) Turning OFF sponge on tracer 03: dye_01
 1.0000E-06  Akt_bak(01)       Background vertical mixing coefficient (m2/s)
                                 for tracer 01: temp
 1.0000E-06  Akt_bak(02)       Background vertical mixing coefficient (m2/s)
                                 for tracer 02: salt
 1.0000E-06  Akt_bak(03)       Background vertical mixing coefficient (m2/s)
                                 for tracer 03: dye_01
 1.0000E-05  Akv_bak           Background vertical mixing coefficient (m2/s)
                                 for momentum.
 5.0000E-06  Akk_bak           Background vertical mixing coefficient (m2/s)
                                 for turbulent energy.
 5.0000E-06  Akp_bak           Background vertical mixing coefficient (m2/s)
                                 for turbulent generic statistical field.
     -1.000  gls_p             GLS stability exponent.
      0.500  gls_m             GLS turbulent kinetic energy exponent.
     -1.000  gls_n             GLS turbulent length scale exponent.
 7.6000E-06  gls_Kmin          GLS minimum value of turbulent kinetic energy.
 1.0000E-12  gls_Pmin          GLS minimum value of dissipation.
 5.4770E-01  gls_cmu0          GLS stability coefficient.
 5.5500E-01  gls_c1            GLS shear production coefficient.
 8.3300E-01  gls_c2            GLS dissipation coefficient.
-6.0000E-01  gls_c3m           GLS stable buoyancy production coefficient.
 1.0000E+00  gls_c3p           GLS unstable buoyancy production coefficient.
 2.0000E+00  gls_sigk          GLS constant Schmidt number for TKE.
 2.0000E+00  gls_sigp          GLS constant Schmidt number for PSI.
   1400.000  charnok_alpha     Charnok factor for Zos calculation.
      0.500  zos_hsig_alpha    Factor for Zos calculation using Hsig(Awave).
      0.250  sz_alpha          Factor for Wave dissipation surface tke flux .
    100.000  crgban_cw         Factor for Craig/Banner surface tke flux.
 3.0000E-04  rdrg              Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2             Quadratic bottom drag coefficient.
 2.0000E-02  Zob               Bottom roughness (m).
 2.0000E-02  Zos               Surface roughness (m).
 2.0000E+00  blk_ZQ            Height (m) of surface air humidity measurement.
 2.0000E+00  blk_ZT            Height (m) of surface air temperature measurement.
 2.0000E+00  blk_ZW            Height (m) of surface winds measurement.
 5.0000E-01  Dcrit             Minimum depth for wetting and drying (m).
          1  lmd_Jwt           Jerlov water type.
          2  Vtransform        S-coordinate transformation equation.
          4  Vstretching       S-coordinate stretching function.
 8.0000E+00  theta_s           S-coordinate surface control parameter.
 4.0000E+00  theta_b           S-coordinate bottom  control parameter.
     20.000  Tcline            S-coordinate surface/bottom layer width (m) used
                                 in vertical coordinate stretching.
   1025.000  rho0              Mean density (kg/m3) for Boussinesq approximation.
      0.000  dstart            Time-stamp assigned to model initialization (days).
       0.00  time_ref          Reference time for units attribute (yyyymmdd.dd)
 3.6000E+02  Tnudg(01)         Nudging/relaxation time scale (days)
                                 for tracer 01: temp
 3.6000E+02  Tnudg(02)         Nudging/relaxation time scale (days)
                                 for tracer 02: salt
 3.6000E+02  Tnudg(03)         Nudging/relaxation time scale (days)
                                 for tracer 03: dye_01
 0.0000E+00  Znudg             Nudging/relaxation time scale (days)
                                 for free-surface.
 0.0000E+00  M2nudg            Nudging/relaxation time scale (days)
                                 for 2D momentum.
 0.0000E+00  M3nudg            Nudging/relaxation time scale (days)
                                 for 3D momentum.
 1.0000E+00  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.
     25.000  T0                Background potential temperature (C) constant.
     35.000  S0                Background salinity (PSU) constant.
      1.000  gamma2            Slipperiness variable: free-slip (1.0) or 
                                                      no-slip (-1.0).
          F  LuvSrc            Turning OFF momentum point Sources/Sinks.
          F  LwSrc             Turning OFF volume influx point Sources/Sinks.
          F  LtracerSrc(01)    Turning OFF point Sources/Sinks on tracer 01: temp
          F  LtracerSrc(02)    Turning OFF point Sources/Sinks on tracer 02: salt
          F  LtracerSrc(03)    Turning OFF  point Sources/Sinks on tracer 03: dye_01
          F  LsshCLM           Turning OFF processing of SSH climatology.
          T  Lm2CLM            Turning ON  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
          T  LtracerCLM(03)    Turning ON  processing of climatology tracer 03: dye_01
          F  LnudgeM2CLM       Turning OFF nudging of 2D momentum climatology.
          F  LnudgeM3CLM       Turning OFF nudging of 3D momentum climatology.
          F  LnudgeTCLM(01)    Turning OFF nudging of climatology tracer 01: temp
          F  LnudgeTCLM(02)    Turning OFF nudging of climatology tracer 02: salt
          F  LnudgeTCLM(03)    Turning OFF nudging of climatology tracer 03: dye_01

          T  Hout(idFsur)      Write out free-surface.
          T  Hout(idUbar)      Write out 2D U-momentum component.
          T  Hout(idVbar)      Write out 2D V-momentum component.
          T  Hout(idUvel)      Write out 3D U-momentum component.
          T  Hout(idVvel)      Write out 3D V-momentum component.
          T  Hout(idWvel)      Write out W-momentum component.
          T  Hout(idOvel)      Write out omega vertical velocity.
          T  Hout(idTvar)      Write out tracer 01: temp
          T  Hout(idTvar)      Write out tracer 02: salt
          T  Hout(inert)       Write out inert passive tracer 01: dye_01
          T  Hout(idDano)      Write out density anomaly.
          T  Hout(idVvis)      Write out vertical viscosity: AKv.
          T  Hout(idTdif)      Write out vertical diffusion: AKt(itemp).
          T  Hout(idSdif)      Write out vertical diffusion: AKt(isalt).
          T  Hout(idMtke)      Write out turbulent kinetic energy.
          T  Hout(idMtls)      Write out turbulent generic length-scale.

          T  Aout(idFsur)      Write out averaged free-surface.
          T  Aout(idUbar)      Write out averaged 2D U-momentum component.
          T  Aout(idVbar)      Write out averaged 2D V-momentum component.
          T  Aout(idUvel)      Write out averaged 3D U-momentum component.
          T  Aout(idVvel)      Write out averaged 3D V-momentum component.
          T  Aout(idWvel)      Write out averaged W-momentum component.
          T  Aout(idOvel)      Write out averaged omega vertical velocity.
          T  Aout(idTvar)      Write out averaged tracer 01: temp
          T  Aout(idTvar)      Write out averaged tracer 02: salt
          T  Aout(inert)       Write out  averaged inert passive tracer 01: dye_01

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

 Output/Input Files:

               Output Restart File:  roms_rst.nc
               Output History File:  roms_yangjiang_his.nc
              Output Averages File:  roms_yangjiang_avg.nc
                   Input Grid File:  /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc
      Input Nonlinear Initial File:  /home/enzo/roms/project/yangjiang_0.3km/data/ic/HYCOM_GLBy0.08_2022_182_ic_YANGJIANG.nc
             Input Forcing File 01:  /home/enzo/roms/project/yangjiang_0.3km/data/Forcings/MERRA_albedo_daily_2022_07.nc
             Input Forcing File 02:  /home/enzo/roms/project/yangjiang_0.3km/data/Forcings/MERRA_cloud_3hours_2022_07.nc
             Input Forcing File 03:  /home/enzo/roms/project/yangjiang_0.3km/data/Forcings/MERRA_lwrad_down_3hours_2022_07.nc
             Input Forcing File 04:  /home/enzo/roms/project/yangjiang_0.3km/data/Forcings/MERRA_Pair_3hours_2022_07.nc
             Input Forcing File 05:  /home/enzo/roms/project/yangjiang_0.3km/data/Forcings/MERRA_Qair_3hours_2022_07.nc
             Input Forcing File 06:  /home/enzo/roms/project/yangjiang_0.3km/data/Forcings/MERRA_rain_3hours_2022_07.nc
             Input Forcing File 07:  /home/enzo/roms/project/yangjiang_0.3km/data/Forcings/MERRA_snow_3hours_2022_07.nc
             Input Forcing File 08:  /home/enzo/roms/project/yangjiang_0.3km/data/Forcings/MERRA_swrad_3hours_2022_07.nc
             Input Forcing File 09:  /home/enzo/roms/project/yangjiang_0.3km/data/Forcings/MERRA_Tair_3hours_2022_07.nc
             Input Forcing File 10:  /home/enzo/roms/project/yangjiang_0.3km/data/Forcings/MERRA_Uwind_3hours_2022_07.nc
             Input Forcing File 11:  /home/enzo/roms/project/yangjiang_0.3km/data/Forcings/MERRA_Vwind_3hours_2022_07.nc
         Input Climatology File 01:  /home/enzo/roms/project/yangjiang_0.3km/data/clm/HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc
    Input Lateral Boundary File 01:  /home/enzo/roms/project/yangjiang_0.3km/data/bdry/HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc
  ROMS I/O variables Metadata File:  /home/enzo/roms/project/yangjiang_0.3km/varinfo.yaml

 Tile partition information for Grid 01:  1098x748x15  tiling: 4x4

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1      274        1      187   768570
        1      275      549        1      187   771375
        2      550      824        1      187   771375
        3      825     1098        1      187   768570
        4        1      274      188      374   768570
        5      275      549      188      374   771375
        6      550      824      188      374   771375
        7      825     1098      188      374   768570
        8        1      274      375      561   768570
        9      275      549      375      561   771375
       10      550      824      375      561   771375
       11      825     1098      375      561   768570
       12        1      274      562      748   768570
       13      275      549      562      748   771375
       14      550      824      562      748   771375
       15      825     1098      562      748   768570

 Tile minimum and maximum fractional coordinates for Grid 01:
   (interior points only)

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50  1100.50     0.50   750.50  RHO-points
        1    -0.50  1100.50     0.50   750.50  RHO-points
        2    -0.50  1100.50     0.50   750.50  RHO-points
        3    -0.50  1099.50     0.50   750.50  RHO-points
        4     0.50  1100.50    -0.50   750.50  RHO-points
        5    -0.50  1100.50    -0.50   750.50  RHO-points
        6    -0.50  1100.50    -0.50   750.50  RHO-points
        7    -0.50  1099.50    -0.50   750.50  RHO-points
        8     0.50  1100.50    -0.50   750.50  RHO-points
        9    -0.50  1100.50    -0.50   750.50  RHO-points
       10    -0.50  1100.50    -0.50   750.50  RHO-points
       11    -0.50  1099.50    -0.50   750.50  RHO-points
       12     0.50  1100.50    -0.50   749.50  RHO-points
       13    -0.50  1100.50    -0.50   749.50  RHO-points
       14    -0.50  1100.50    -0.50   749.50  RHO-points
       15    -0.50  1099.50    -0.50   749.50  RHO-points

        0     0.00  1100.50     0.50   750.50    U-points
        1    -0.50  1100.50     0.50   750.50    U-points
        2    -0.50  1100.50     0.50   750.50    U-points
        3    -0.50  1099.00     0.50   750.50    U-points
        4     0.00  1100.50    -0.50   750.50    U-points
        5    -0.50  1100.50    -0.50   750.50    U-points
        6    -0.50  1100.50    -0.50   750.50    U-points
        7    -0.50  1099.00    -0.50   750.50    U-points
        8     0.00  1100.50    -0.50   750.50    U-points
        9    -0.50  1100.50    -0.50   750.50    U-points
       10    -0.50  1100.50    -0.50   750.50    U-points
       11    -0.50  1099.00    -0.50   750.50    U-points
       12     0.00  1100.50    -0.50   749.50    U-points
       13    -0.50  1100.50    -0.50   749.50    U-points
       14    -0.50  1100.50    -0.50   749.50    U-points
       15    -0.50  1099.00    -0.50   749.50    U-points

        0     0.50  1100.50     0.00   750.50    V-points
        1    -0.50  1100.50     0.00   750.50    V-points
        2    -0.50  1100.50     0.00   750.50    V-points
        3    -0.50  1099.50     0.00   750.50    V-points
        4     0.50  1100.50    -0.50   750.50    V-points
        5    -0.50  1100.50    -0.50   750.50    V-points
        6    -0.50  1100.50    -0.50   750.50    V-points
        7    -0.50  1099.50    -0.50   750.50    V-points
        8     0.50  1100.50    -0.50   750.50    V-points
        9    -0.50  1100.50    -0.50   750.50    V-points
       10    -0.50  1100.50    -0.50   750.50    V-points
       11    -0.50  1099.50    -0.50   750.50    V-points
       12     0.50  1100.50    -0.50   749.00    V-points
       13    -0.50  1100.50    -0.50   749.00    V-points
       14    -0.50  1100.50    -0.50   749.00    V-points
       15    -0.50  1099.50    -0.50   749.00    V-points

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

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

 temp                     1   Akima4             Akima4

 salt                     1   MPDATA             MPDATA

 dye_01                   1   HSIMT              HSIMT

 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   Chapman Exp  Chapman Exp  Chapman Exp  Chapman Exp

 ubar                     1   Shchepetkin  Shchepetkin  Shchepetkin  Shchepetkin

 vbar                     1   Shchepetkin  Shchepetkin  Shchepetkin  Shchepetkin

 u                        1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 v                        1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 temp                     1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 salt                     1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 dye_01                   1   Clamped      Clamped      Radiation    Clamped

 tke                      1   Radiation    Radiation    Radiation    Radiation

 Activated C-preprocessing Options:

 RIVERPLUME1              River Plume Test 1, flow across u-face grid cell
 ANA_BPFLUX               Analytical bottom passive tracers fluxes
 ANA_BSFLUX               Analytical kinematic bottom salinity flux
 ANA_BTFLUX               Analytical kinematic bottom temperature flux
 ANA_SPFLUX               Analytical surface passive tracer fluxes
 ASSUMED_SHAPE            Using assumed-shape arrays
 AVERAGES                 Writing out time-averaged nonlinear model fields
 BULK_FLUXES              Surface bulk fluxes parameterization
 CURVGRID                 Orthogonal curvilinear grid
 DJ_GRADPS                Parabolic Splines density Jacobian (Shchepetkin, 2002)
 DOUBLE_PRECISION         Double precision arithmetic numerical kernel.
 EMINUSP                  Compute Salt Flux using E-P
 GLS_MIXING               Generic Length-Scale turbulence closure
 LIMIT_BSTRESS            Limit bottom stress to maintain bottom velocity direction
 KANTHA_CLAYSON           Kantha and Clayson stability function formulation
 LONGWAVE_OUT             Compute outgoing longwave radiation internally
 MASKING                  Land/Sea masking
 MIX_GEO_TS               Mixing of tracers along geopotential surfaces
 MIX_S_TS                 Mixing of tracers along constant S-surfaces
 MIX_GEO_UV               Mixing of momentum along geopotential surfaces
 MIX_S_UV                 Mixing of momentum along constant S-surfaces
 NONLINEAR                Nonlinear Model
 NONLIN_EOS               Nonlinear Equation of State for seawater
 N2S2_HORAVG              Horizontal smoothing of buoyancy and shear
 OUT_DOUBLE               Double precision output fields in NetCDF files
 _OPENMP                  OpenMP parallel shared-memory directives
 PERFECT_RESTART          Processing perfect restart variables
 POWER_LAW                Power-law shape time-averaging barotropic filter
 PROFILE                  Time profiling activated
 K_GSCHEME                Third-order upstream advection of TKE fields
 RADIATION_2D             Use tangential phase speed in radiation conditions
 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
 SPLINES_VDIFF            Parabolic Spline Reconstruction for Vertical Diffusion
 SPLINES_VVISC            Parabolic Spline Reconstruction for Vertical Viscosity
 T_PASSIVE                Advecting and diffusing inert passive tracer
 TS_MPDATA_LIMIT          Further limiter in upwind corrector fluxes for stability
 TS_DIF2                  Harmonic mixing of tracers
 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_LDRAG                 Linear bottom stress
 UV_VIS2                  Harmonic mixing of momentum
 VAR_RHO_2D               Variable density barotropic mode
 WET_DRY                  Wetting and drying activated

 Process Information:

 Thread #    0 (pid=   62634) is active.
 Thread #    5 (pid=   62634) is active.
 Thread #    8 (pid=   62634) is active.
 Thread #   15 (pid=   62634) is active.
 Thread #    7 (pid=   62634) is active.
 Thread #    6 (pid=   62634) is active.
 Thread #   14 (pid=   62634) is active.
 Thread #   11 (pid=   62634) is active.
 Thread #   10 (pid=   62634) is active.
 Thread #    2 (pid=   62634) is active.
 Thread #    4 (pid=   62634) is active.
 Thread #    3 (pid=   62634) is active.
 Thread #   12 (pid=   62634) is active.
 Thread #   13 (pid=   62634) is active.
 Thread #    1 (pid=   62634) is active.
 Thread #    9 (pid=   62634) is active.

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

  GET_GRID_NF90    - angle between XI-axis and EAST: angler
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min = -2.89273051E-10 Max =  2.88658224E-10)
  GET_GRID_NF90    - bathymetry at RHO-points: h
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  2.00000000E+00 Max =  1.14589956E+02)
  GET_GRID_NF90    - Coriolis parameter at RHO-points: f
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  5.05036781E-05 Max =  5.51574581E-05)
  GET_GRID_NF90    - reciprocal XI-grid spacing: pm
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  3.37580182E-03 Max =  3.42103594E-03)
  GET_GRID_NF90    - reciprocal ETA-grid spacing: pn
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  3.38823425E-03 Max =  3.47951942E-03)
  GET_GRID_NF90    - XI-derivative of inverse metric factor pn: dndx
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min = -2.35068885E-05 Max =  2.35176815E-05)
  GET_GRID_NF90    - ETA-derivative of inverse metric factor pm: dmde
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min = -5.40388396E-03 Max =  0.00000000E+00)
  GET_GRID_NF90    - longitude of RHO-points: lon_rho
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  1.10541779E+02 Max =  1.13662488E+02)
  GET_GRID_NF90    - longitude of U-points: lon_u
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  1.10543199E+02 Max =  1.13661068E+02)
  GET_GRID_NF90    - longitude of V-points: lon_v
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  1.10541779E+02 Max =  1.13662488E+02)
  GET_GRID_NF90    - longitude of PSI-points: lon_psi
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  1.10543199E+02 Max =  1.13661068E+02)
  GET_GRID_NF90    - latitude of RHO-points lat_rho
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  2.02608738E+01 Max =  2.22225578E+01)
  GET_GRID_NF90    - latitude of U-points: lat_u
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  2.02608738E+01 Max =  2.22225578E+01)
  GET_GRID_NF90    - latitude of V-points: lat_v
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  2.02622009E+01 Max =  2.22212656E+01)
  GET_GRID_NF90    - latitude of PSI-points lat_psi
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  2.02622009E+01 Max =  2.22212656E+01)
  GET_GRID_NF90    - mask on RHO-points: mask_rho
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  0.00000000E+00 Max =  1.00000000E+00)
  GET_GRID_NF90    - mask on U-points: mask_u
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  0.00000000E+00 Max =  1.00000000E+00)
  GET_GRID_NF90    - mask on V-points: mask_v
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.nc)
                      (Min =  0.00000000E+00 Max =  1.00000000E+00)
  GET_GRID_NF90    - mask on PSI-points: mask_psi
                      (Grid = 01, File: /home/enzo/roms/project/yangjiang_0.3km/yangjiang_0.3km.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

    15   0.0000000   0.0000000        0.000       0.000       0.000       0.000
    14  -0.0666667  -0.0003983       -0.121      -0.671      -1.010      -1.174
    13  -0.1333333  -0.0017081       -0.243      -1.350      -2.060      -2.437
    12  -0.2000000  -0.0043062       -0.364      -2.043      -3.165      -3.826
    11  -0.2666667  -0.0089353       -0.486      -2.756      -4.359      -5.413
    10  -0.3333333  -0.0169029       -0.609      -3.502      -5.697      -7.325
     9  -0.4000000  -0.0304126       -0.733      -4.304      -7.276      -9.778
     8  -0.4666667  -0.0530642       -0.858      -5.197      -9.252     -13.123
     7  -0.5333333  -0.0905094       -0.986      -6.238     -11.870     -17.912
     6  -0.6000000  -0.1510347       -1.118      -7.510     -15.490     -24.952
     5  -0.6666667  -0.2452550       -1.257      -9.119     -20.572     -35.280
     4  -0.7333333  -0.3828872       -1.403     -11.162     -27.539     -49.842
     3  -0.8000000  -0.5634305       -1.557     -13.634     -36.368     -68.592
     2  -0.8666667  -0.7610631       -1.714     -16.277     -45.939     -89.008
     1  -0.9333333  -0.9210683       -1.864     -18.544     -53.876    -105.754
     0  -1.0000000  -1.0000000       -2.000     -20.000     -58.295    -114.590

 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.0022367614257356 0.4026337758325831
 13 0.0060106451121704 0.0334078920475245 0.0037738836864348 0.4360416678801076
 14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905598
 15 0.0099910761708920 0.0329439123123590 0.0216737068001213 0.5021931174029188
 16 0.0122483446563884 0.0326108764399960 0.0339220514565097 0.5348039938429148
 17 0.0146692120341107 0.0322025982847830 0.0485912634906204 0.5670065921276978
 18 0.0172400033810439 0.0317136245503127 0.0658312668716643 0.5987202166780105
 19 0.0199444086685725 0.0311389577709445 0.0857756755402368 0.6298591744489550
 20 0.0227631639997064 0.0304741441486588 0.1085388395399432 0.6603333185976138
 21 0.0256737146312911 0.0297153720153352 0.1342125541712342 0.6900486906129490
 22 0.0286498597812016 0.0288595815276255 0.1628624139524359 0.7189082721405746
 23 0.0316613792205220 0.0279045862015855 0.1945237931729578 0.7468128583421600
 24 0.0346736416507075 0.0268492068942347 0.2291974348236653 0.7736620652363948
 25 0.0376471948657328 0.0256934188392112 0.2668446296893981 0.7993554840756060
 26 0.0405373376992232 0.0244385123436867 0.3073819673886213 0.8237939964192927
 27 0.0432936737565710 0.0230872677537126 0.3506756411451924 0.8468812641730054
 28 0.0458596469320356 0.0216441452951603 0.3965352880772280 0.8685254094681656
 29 0.0481720587108284 0.0201154903974257 0.4447073467880565 0.8886408998655914
 30 0.0501605672561820 0.0185097551070648 0.4948679140442384 0.9071506549726561
 31 0.0517471682814030 0.0168377361985254 0.5466150823256415 0.9239883911711815
 32 0.0528456577069106 0.0151128305891453 0.5994607400325521 0.9391012217603267
 33 0.0533610761022577 0.0133513086655816 0.6528218161348098 0.9524525304259084
 34 0.0531891349131379 0.0115726061288397 0.7060109510479478 0.9640251365547481
 35 0.0522156244733761 0.0097996349650684 0.7582265755213239 0.9738247715198165
 36 0.0503158038019030 0.0080591141492892 0.8085423793232269 0.9818838856691057
 37 0.0473537721847153 0.0063819206892258 0.8558961515079423 0.9882658063583315
 38 0.0431818225418188 0.0048034616164019 0.8990779740497611 0.9930692679747334
 39 0.0376397765791564 0.0033640675316746 0.9367177506289175 0.9964333355064080
 40 0.0305543017255206 0.0021094083123694 0.9672720523544381 0.9985427438187774
 41 0.0217382098544504 0.0010909315881854 0.9890102622088885 0.9996336754069628
 42 0.0109897377911118 0.0003663245930371 1.0000000000000004 0.9999999999999999

 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 =  2.92309119E-01 km    Water points =  2.92309119E-01 km
 Maximum X-grid spacing, DXmax =  2.96225920E-01 km    Water points =  2.96225920E-01 km
 Minimum Y-grid spacing, DYmin =  2.87396011E-01 km    Water points =  2.87396011E-01 km
 Maximum Y-grid spacing, DYmax =  2.95138980E-01 km    Water points =  2.95138980E-01 km
 Minimum Z-grid spacing, DZmin =  1.21284539E-01 m     Water points =  1.21284539E-01 m
 Maximum Z-grid spacing, DZmax =  2.04166190E+01 m     Water points =  2.04166190E+01 m

 Minimum barotropic Courant Number =  7.07604957E-03
 Maximum barotropic Courant Number =  5.34537607E-02
 Maximum Coriolis   Courant Number =  5.51574581E-04


  GET_STATE_NF90   - NLM: state initial conditions,                       2022-07-01 00:00:00.00
                      (Grid 01, t = 738336.0000, File: HYCOM_GLBy0.08_2022_182_ic_YANGJIANG.nc, Rec=0001, Index=1)
                   - free-surface
                      (Min =  5.52000046E-01 Max =  1.32399814E+00)
                   - vertically integrated u-momentum component
                      (Min = -1.49531584E+00 Max = -4.81810697E-03)
                   - vertically integrated v-momentum component
                      (Min = -1.48637830E+00 Max =  3.42725133E-01)
                   - u-momentum component
                      (Min = -1.51859777E+00 Max =  2.22870557E-01)
                   - v-momentum component
                      (Min = -1.08283354E+00 Max =  8.55180437E-01)
                   - potential temperature
                      (Min =  1.84589802E+01 Max =  3.08091283E+01)
                   - salinity
                      (Min =  6.02299976E+00 Max =  3.45300093E+01)
                   - dye concentration, type 01
                      (Min =  0.00000000E+00 Max =  1.00000000E+04)
  GET_2DFLD_NF90   - surface u-wind component,                            2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: MERRA_Uwind_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.0000
                      (Min = -1.29363509E+01 Max = -2.07263070E+00)   regrid = T
  GET_2DFLD_NF90   - surface v-wind component,                            2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: MERRA_Vwind_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.0000
                      (Min = -1.43423418E+01 Max =  2.25609907E+01)   regrid = T
  GET_2DFLD_NF90   - surface air pressure,                                2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: MERRA_Pair_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.0000
                      (Min =  9.93625663E+04 Max =  1.00158278E+05)   regrid = T
  GET_2DFLD_NF90   - net solar shortwave radiation flux,                  2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: MERRA_swrad_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.0000
                      (Min =  3.62260939E-05 Max =  1.23798589E-04)   regrid = T
  GET_2DFLD_NF90   - downward longwave radiation flux,                    2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: MERRA_lwrad_down_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.0000
                      (Min =  1.00654447E-04 Max =  1.11343561E-04)   regrid = T
  GET_2DFLD_NF90   - surface air temperature,                             2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: MERRA_Tair_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.0000
                      (Min =  2.66196459E+01 Max =  2.93892408E+01)   regrid = T
  GET_2DFLD_NF90   - surface air relative humidity,                       2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: MERRA_Qair_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.0000
                      (Min =  1.97108250E-04 Max =  2.27929871E-04)   regrid = T
  GET_2DFLD_NF90   - rain fall rate,                                      2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: MERRA_rain_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.0000
                      (Min =  3.47117136E-05 Max =  3.04730625E-04)   regrid = T
  GET_NGFLD_NF90   - free-surface western boundary condition,             2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  0.00000000E+00 Max =  1.17714443E+00)
  GET_NGFLD_NF90   - free-surface eastern boundary condition,             2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  0.00000000E+00 Max =  1.20872233E+00)
  GET_NGFLD_NF90   - free-surface southern boundary condition,            2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  5.64784915E-01 Max =  9.96285977E-01)
  GET_NGFLD_NF90   - free-surface northern boundary condition,            2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  0.00000000E+00 Max =  1.17177368E+00)
  GET_NGFLD_NF90   - 2D u-momentum western boundary condition,            2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -6.16859707E-01 Max =  0.00000000E+00)
  GET_NGFLD_NF90   - 2D v-momentum western boundary condition,            2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -8.60267046E-01 Max =  0.00000000E+00)
  GET_NGFLD_NF90   - 2D u-momentum eastern boundary condition,            2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -8.49212362E-01 Max =  0.00000000E+00)
  GET_NGFLD_NF90   - 2D v-momentum eastern boundary condition,            2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -3.49940805E-01 Max =  1.43088329E-01)
  GET_NGFLD_NF90   - 2D u-momentum southern boundary condition,           2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -6.74444344E-01 Max = -6.04237111E-02)
  GET_NGFLD_NF90   - 2D v-momentum southern boundary condition,           2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -6.57270570E-01 Max =  1.23478427E-01)
  GET_NGFLD_NF90   - 2D u-momentum northern boundary condition,           2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -3.77390748E-01 Max =  0.00000000E+00)
  GET_NGFLD_NF90   - 2D v-momentum northern boundary condition,           2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -2.04799431E-01 Max =  1.46082149E-01)
  GET_NGFLD_NF90   - 3D u-momentum western boundary condition,            2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -6.77564960E-01 Max =  1.79639917E-02)
  GET_NGFLD_NF90   - 3D v-momentum western boundary condition,            2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -5.99894119E-01 Max =  0.00000000E+00)
  GET_NGFLD_NF90   - 3D u-momentum eastern boundary condition,            2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -1.11644364E+00 Max =  2.15449776E-02)
  GET_NGFLD_NF90   - 3D v-momentum eastern boundary condition,            2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -3.30794819E-01 Max =  8.55180437E-01)
  GET_NGFLD_NF90   - 3D u-momentum southern boundary condition,           2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -9.47471585E-01 Max =  2.00942307E-01)
  GET_NGFLD_NF90   - 3D v-momentum southern boundary condition,           2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -1.08283354E+00 Max =  8.55180437E-01)
  GET_NGFLD_NF90   - 3D u-momentum northern boundary condition,           2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -2.54918646E-01 Max =  0.00000000E+00)
  GET_NGFLD_NF90   - 3D v-momentum northern boundary condition,           2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -1.31757613E-01 Max =  1.24353775E-01)
  GET_NGFLD_NF90   - temperature western boundary condition,              2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  0.00000000E+00 Max =  3.02530003E+01)
  GET_NGFLD_NF90   - salinity western boundary condition,                 2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  0.00000000E+00 Max =  3.36716298E+01)
  GET_NGFLD_NF90   - dye concentration western boundary condition, type 01****-12-02 00:00:00.00
                      (Grid= 01, Rec=59, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin= *************** Tmax= ***************)   t = ***************
                      (Min =  0.00000000E+00 Max =  0.00000000E+00)
  GET_NGFLD_NF90   - temperature eastern boundary condition,              2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  0.00000000E+00 Max =  2.91520927E+01)
  GET_NGFLD_NF90   - salinity eastern boundary condition,                 2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  0.00000000E+00 Max =  3.45300093E+01)
  GET_NGFLD_NF90   - temperature southern boundary condition,             2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  1.86646123E+01 Max =  2.95480420E+01)
  GET_NGFLD_NF90   - salinity southern boundary condition,                2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  3.36658164E+01 Max =  3.45142197E+01)
  GET_NGFLD_NF90   - dye concentration southern boundary condition, type 0****-12-02 00:00:00.00
                      (Grid= 01, Rec=59, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin= *************** Tmax= ***************)   t = ***************
                      (Min =  0.00000000E+00 Max =  0.00000000E+00)
  GET_NGFLD_NF90   - temperature northern boundary condition,             2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  0.00000000E+00 Max =  2.92455689E+01)
  GET_NGFLD_NF90   - salinity northern boundary condition,                2022-07-01 00:00:00.00
                      (Grid= 01, Rec=1, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  0.00000000E+00 Max =  2.13523231E+01)
  GET_NGFLD_NF90   - dye concentration northern boundary condition, type 0****-12-02 00:00:00.00
                      (Grid= 01, Rec=59, Index=2, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin= *************** Tmax= ***************)   t = ***************
                      (Min =  0.00000000E+00 Max =  0.00000000E+00)
  GET_2DFLD_NF90   - vertically integrated u-momentum climatology,        2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -1.49531584E+00 Max = -4.81810697E-03)   regrid = F
  GET_2DFLD_NF90   - vertically integrated v-momentum climatology,        2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -1.48637830E+00 Max =  3.42725133E-01)   regrid = F
  GET_3DFLD_NF90   - u-momentum component climatology,                    2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -1.51859777E+00 Max =  2.22870557E-01)
  GET_3DFLD_NF90   - v-momentum component climatology,                    2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min = -1.08283354E+00 Max =  8.55180437E-01)
  GET_3DFLD_NF90   - potential temperature climatology,                   2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  1.84589802E+01 Max =  3.08091283E+01)
  GET_3DFLD_NF90   - salinity climatology,                                2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  6.02299976E+00 Max =  3.45300093E+01)
  GET_3DFLD_NF90   - dye concentration, type 01 climatology,              2022-07-01 00:00:00.00
                      (Grid=01, Rec=1, Index=1, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738336.0000
                      (Min =  0.00000000E+00 Max =  0.00000000E+00)

 Basin information for Grid 01:

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

 Initial domain volumes:  TotVolume =  2.1803736787E+12 m3
                         MinCellVol =  1.0189297887E+04 m3
                         MaxCellVol =  1.7681775255E+06 m3
                            Max/Min =  1.7353281306E+02



Zeng
Posts: 26
Joined: Mon Sep 19, 2022 1:06 pm
Location: Sun Yat-Sen University

Re: BLOWUP. REASON: MaxSpeed = 2.1E+01

#3 Unread post by Zeng »

Code: Select all

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000000000001 - 000000051840)

  GET_2DFLD_NF90   - surface u-wind component,                            2022-07-01 03:00:00.00
                      (Grid=01, Rec=2, Index=2, File: MERRA_Uwind_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.1250
                      (Min = -1.33804505E+01 Max = -8.22904751E-01)   regrid = T
  GET_2DFLD_NF90   - surface v-wind component,                            2022-07-01 03:00:00.00
                      (Grid=01, Rec=2, Index=2, File: MERRA_Vwind_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.1250
                      (Min = -1.28794699E+01 Max =  2.12873147E+01)   regrid = T
  GET_2DFLD_NF90   - surface air pressure,                                2022-07-01 03:00:00.00
                      (Grid=01, Rec=2, Index=2, File: MERRA_Pair_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.1250
                      (Min =  9.91153245E+04 Max =  9.99855503E+04)   regrid = T
  GET_2DFLD_NF90   - net solar shortwave radiation flux,                  2022-07-01 03:00:00.00
                      (Grid=01, Rec=2, Index=2, File: MERRA_swrad_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.1250
                      (Min =  4.77848364E-05 Max =  1.62826668E-04)   regrid = T
  GET_2DFLD_NF90   - downward longwave radiation flux,                    2022-07-01 03:00:00.00
                      (Grid=01, Rec=2, Index=2, File: MERRA_lwrad_down_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.1250
                      (Min =  1.02289515E-04 Max =  1.11155555E-04)   regrid = T
  GET_2DFLD_NF90   - surface air temperature,                             2022-07-01 03:00:00.00
                      (Grid=01, Rec=2, Index=2, File: MERRA_Tair_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.1250
                      (Min =  2.81791995E+01 Max =  3.00772119E+01)   regrid = T
  GET_2DFLD_NF90   - surface air relative humidity,                       2022-07-01 03:00:00.00
                      (Grid=01, Rec=2, Index=2, File: MERRA_Qair_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.1250
                      (Min =  1.99293913E-04 Max =  2.28549817E-04)   regrid = T
  GET_2DFLD_NF90   - rain fall rate,                                      2022-07-01 03:00:00.00
                      (Grid=01, Rec=2, Index=2, File: MERRA_rain_3hours_2022_07.nc)
                      (Tmin=     738336.0000 Tmax=     738366.8750)   t =     738336.1250
                      (Min =  6.06189075E-05 Max =  1.00711971E-03)   regrid = T
  GET_NGFLD_NF90   - free-surface western boundary condition,             2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -1.80538757E-01 Max =  4.81369376E-01)
  GET_NGFLD_NF90   - free-surface eastern boundary condition,             2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  0.00000000E+00 Max =  8.80197064E-01)
  GET_NGFLD_NF90   - free-surface southern boundary condition,            2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  1.85417025E-01 Max =  5.27890884E-01)
  GET_NGFLD_NF90   - free-surface northern boundary condition,            2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  0.00000000E+00 Max =  8.80197064E-01)
  GET_NGFLD_NF90   - 2D u-momentum western boundary condition,            2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  0.00000000E+00 Max =  7.20607964E-01)
  GET_NGFLD_NF90   - 2D v-momentum western boundary condition,            2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -1.06253458E-01 Max =  2.37847358E-01)
  GET_NGFLD_NF90   - 2D u-momentum eastern boundary condition,            2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -1.14490642E-01 Max =  2.97991811E-01)
  GET_NGFLD_NF90   - 2D v-momentum eastern boundary condition,            2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -6.18999842E-02 Max =  2.08595726E-01)
  GET_NGFLD_NF90   - 2D u-momentum southern boundary condition,           2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  5.64850477E-02 Max =  4.86453052E-01)
  GET_NGFLD_NF90   - 2D v-momentum southern boundary condition,           2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -8.66925456E-01 Max =  2.81053999E-01)
  GET_NGFLD_NF90   - 2D u-momentum northern boundary condition,           2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -4.17115715E-02 Max =  1.01922323E-02)
  GET_NGFLD_NF90   - 2D v-momentum northern boundary condition,           2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -4.75553479E-03 Max =  3.18649438E-02)
  GET_NGFLD_NF90   - 3D u-momentum western boundary condition,            2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  0.00000000E+00 Max =  6.99432727E-01)
  GET_NGFLD_NF90   - 3D v-momentum western boundary condition,            2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -1.39923209E-01 Max =  2.86042405E-01)
  GET_NGFLD_NF90   - 3D u-momentum eastern boundary condition,            2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -1.20662971E-01 Max =  6.70910248E-01)
  GET_NGFLD_NF90   - 3D v-momentum eastern boundary condition,            2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -4.16973931E-01 Max =  3.25142065E-01)
  GET_NGFLD_NF90   - 3D u-momentum southern boundary condition,           2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -5.93218526E-01 Max =  8.50149647E-01)
  GET_NGFLD_NF90   - 3D v-momentum southern boundary condition,           2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -1.03998822E+00 Max =  5.53817305E-01)
  GET_NGFLD_NF90   - 3D u-momentum northern boundary condition,           2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -3.71075261E-02 Max =  1.69753927E-02)
  GET_NGFLD_NF90   - 3D v-momentum northern boundary condition,           2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -5.38510554E-02 Max =  4.28894102E-02)
  GET_NGFLD_NF90   - temperature western boundary condition,              2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  0.00000000E+00 Max =  2.84980864E+01)
  GET_NGFLD_NF90   - salinity western boundary condition,                 2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  0.00000000E+00 Max =  3.35950088E+01)
  GET_NGFLD_NF90   - temperature eastern boundary condition,              2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  0.00000000E+00 Max =  2.81083902E+01)
  GET_NGFLD_NF90   - salinity eastern boundary condition,                 2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  0.00000000E+00 Max =  3.45448294E+01)
  GET_NGFLD_NF90   - temperature southern boundary condition,             2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  1.83822147E+01 Max =  2.79347030E+01)
  GET_NGFLD_NF90   - salinity southern boundary condition,                2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  3.35886191E+01 Max =  3.45448294E+01)
  GET_NGFLD_NF90   - temperature northern boundary condition,             2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  0.00000000E+00 Max =  2.81641827E+01)
  GET_NGFLD_NF90   - salinity northern boundary condition,                2022-07-02 00:00:00.00
                      (Grid= 01, Rec=2, Index=1, File: HYCOM_GLBy0.08_2022_182_188_bdry_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  0.00000000E+00 Max =  2.40835101E+01)
  GET_2DFLD_NF90   - vertically integrated u-momentum climatology,        2022-07-02 00:00:00.00
                      (Grid=01, Rec=2, Index=2, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -3.48645935E-01 Max =  7.20607964E-01)   regrid = F
  GET_2DFLD_NF90   - vertically integrated v-momentum climatology,        2022-07-02 00:00:00.00
                      (Grid=01, Rec=2, Index=2, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -8.68093244E-01 Max =  3.54839193E-01)   regrid = F
  GET_3DFLD_NF90   - u-momentum component climatology,                    2022-07-02 00:00:00.00
                      (Grid=01, Rec=2, Index=2, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -5.93218526E-01 Max =  9.05053627E-01)
  GET_3DFLD_NF90   - v-momentum component climatology,                    2022-07-02 00:00:00.00
                      (Grid=01, Rec=2, Index=2, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min = -1.03998822E+00 Max =  6.28454264E-01)
  GET_3DFLD_NF90   - potential temperature climatology,                   2022-07-02 00:00:00.00
                      (Grid=01, Rec=2, Index=2, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  1.83792989E+01 Max =  2.94984867E+01)
  GET_3DFLD_NF90   - salinity climatology,                                2022-07-02 00:00:00.00
                      (Grid=01, Rec=2, Index=2, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  8.65799904E+00 Max =  3.45448294E+01)
  GET_3DFLD_NF90   - dye concentration, type 01 climatology,              2022-07-02 00:00:00.00
                      (Grid=01, Rec=2, Index=2, File: HYCOM_GLBy0.08_2022_182_188_clim_YANGJIANG.nc)
                      (Tmin=     738336.0000 Tmax=     738342.0000)   t =     738337.0000
                      (Min =  0.00000000E+00 Max =  0.00000000E+00)

 TIME-STEP YYYY-MM-DD hh:mm:ss.ss  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME
                     C => (i,j,k)       Cu            Cv            Cw         Max Speed

         0 2022-07-01 00:00:00.00  1.529804E-01  2.692071E+02  2.693600E+02  2.259602E+12
                    (0823,463,15)  4.668495E-02  3.119485E-02  0.000000E+00  1.643748E+00
  DEF_HIS_NF90     - creating history file,            Grid 01: roms_yangjiang_his.nc
  WRT_HIS_NF90     - wrote history     fields (Index=1,1) in record = 1
  DEF_AVG_NF90     - creating average file,            Grid 01: roms_yangjiang_avg.nc
  DEF_RST_NF90     - creating restart file,            Grid 01: roms_rst.nc
         1 2022-07-01 00:00:10.00  1.741102E-01  2.724537E+02  2.726278E+02  2.232229E+12
                    (0761,526,01)  2.712961E-02  7.657837E-04  1.286470E-01  3.628447E+00
         2 2022-07-01 00:00:20.00  2.092062E-01  2.724539E+02  2.726632E+02  2.232231E+12
                    (0821,509,15)  1.336461E-02  7.444934E-03  8.878534E-01  5.684815E+00
         3 2022-07-01 00:00:30.00  2.556127E-01  2.724543E+02  2.727099E+02  2.232236E+12
                    (0821,509,15)  4.178283E-03  1.076565E-02  7.707206E-01  8.151705E+00
         4 2022-07-01 00:00:40.00  3.117414E-01  2.724547E+02  2.727664E+02  2.232243E+12
                    (0745,511,15)  1.907655E-02  4.060388E-03  6.185898E-01  1.041183E+01
         5 2022-07-01 00:00:50.00  3.732515E-01  2.724549E+02  2.728282E+02  2.232252E+12
                    (0029,257,15)  5.974186E-03  2.146107E-02  6.202628E-01  1.240293E+01
         6 2022-07-01 00:01:00.00  4.356067E-01  2.724553E+02  2.728909E+02  2.232263E+12
                    (0029,257,15)  7.799759E-03  3.442786E-02  6.546900E-01  1.405953E+01
         7 2022-07-01 00:01:10.00  4.951391E-01  2.724556E+02  2.729507E+02  2.232277E+12
                    (0030,257,15)  2.125948E-02  5.779642E-02  6.609452E-01  1.531732E+01
         8 2022-07-01 00:01:20.00  5.489688E-01  2.724559E+02  2.730048E+02  2.232293E+12
                    (0026,259,15)  2.713743E-02  1.782318E-01  6.887762E-01  1.610815E+01
         9 2022-07-01 00:01:30.00  5.958170E-01  2.724562E+02  2.730520E+02  2.232311E+12
                    (0026,259,15)  3.517906E-02  2.251958E-01  9.179379E-01  1.644753E+01
        10 2022-07-01 00:01:40.00  6.351046E-01  2.724564E+02  2.730916E+02  2.232331E+12
                    (0026,259,15)  4.254667E-02  2.589368E-01  1.144329E+00  1.638215E+01
        11 2022-07-01 00:01:50.00  6.679165E-01  2.724567E+02  2.731246E+02  2.232353E+12
                    (0030,258,15)  1.013916E-01  4.401642E-01  1.231549E+00  1.706742E+01
        12 2022-07-01 00:02:00.00  6.947804E-01  2.724570E+02  2.731518E+02  2.232378E+12
                    (0030,258,15)  1.022718E-01  4.431668E-01  1.489792E+00  1.774825E+01
        13 2022-07-01 00:02:10.00  7.165123E-01  2.724573E+02  2.731738E+02  2.232404E+12
                    (0030,258,15)  9.821906E-02  4.263160E-01  1.633539E+00  1.768254E+01
        14 2022-07-01 00:02:20.00  7.338034E-01  2.724576E+02  2.731914E+02  2.232433E+12
                    (0030,258,15)  9.116294E-02  3.966399E-01  1.698958E+00  1.765744E+01
        15 2022-07-01 00:02:30.00  7.471761E-01  2.724579E+02  2.732051E+02  2.232464E+12
                    (0030,258,15)  8.234573E-02  3.602411E-01  1.713886E+00  1.720007E+01
        16 2022-07-01 00:02:40.00  7.571425E-01  2.724582E+02  2.732153E+02  2.232498E+12
                    (0030,258,15)  7.311570E-02  3.225189E-01  1.692012E+00  1.586664E+01
        17 2022-07-01 00:02:50.00  7.644044E-01  2.724585E+02  2.732229E+02  2.232533E+12
                    (0030,258,15)  6.448894E-02  2.861387E-01  1.650649E+00  1.549316E+01
        18 2022-07-01 00:03:00.00  7.695183E-01  2.724588E+02  2.732283E+02  2.232571E+12
                    (0030,258,15)  5.657403E-02  2.516277E-01  1.604237E+00  1.494836E+01
        19 2022-07-01 00:03:10.00  7.730953E-01  2.724592E+02  2.732323E+02  2.232610E+12
                    (0030,258,15)  4.957922E-02  2.195344E-01  1.557228E+00  1.343501E+01
        20 2022-07-01 00:03:20.00  7.755271E-01  2.724595E+02  2.732351E+02  2.232652E+12
                    (0025,226,15)  1.546987E-01  3.280008E-01  1.405925E+00  1.238437E+01
        21 2022-07-01 00:03:30.00  7.780596E-01  2.724599E+02  2.732380E+02  2.232696E+12
                    (0025,226,15)  1.517887E-01  3.018653E-01  1.451394E+00  1.119239E+01
        22 2022-07-01 00:03:40.00  7.807640E-01  2.724604E+02  2.732411E+02  2.232743E+12
                    (0025,226,15)  1.493734E-01  2.798796E-01  1.419853E+00  1.035727E+01
        23 2022-07-01 00:03:50.00  7.843435E-01  2.724608E+02  2.732452E+02  2.232791E+12
                    (0017,259,15)  4.063307E-02  1.972063E-01  1.593838E+00  9.581826E+00
        24 2022-07-01 00:04:00.00  7.897136E-01  2.724614E+02  2.732511E+02  2.232841E+12
                    (0017,259,15)  4.227071E-02  2.042227E-01  1.698398E+00  8.873838E+00
        25 2022-07-01 00:04:10.00  7.960992E-01  2.724620E+02  2.732581E+02  2.232894E+12
                    (0017,259,15)  4.399571E-02  2.122367E-01  1.794164E+00  8.242543E+00
        26 2022-07-01 00:04:20.00  8.031662E-01  2.724628E+02  2.732659E+02  2.232949E+12
                    (0017,259,15)  4.587711E-02  2.206060E-01  1.894746E+00  8.117614E+00
        27 2022-07-01 00:04:30.00  8.116877E-01  2.724636E+02  2.732753E+02  2.233006E+12
                    (0017,259,15)  4.797769E-02  2.291749E-01  2.011838E+00  8.398902E+00
        28 2022-07-01 00:04:40.00  8.217641E-01  2.724646E+02  2.732863E+02  2.233065E+12
                    (0017,259,15)  5.023249E-02  2.375838E-01  2.156956E+00  8.766627E+00
        29 2022-07-01 00:04:50.00  8.325019E-01  2.724656E+02  2.732981E+02  2.233127E+12
                    (0017,259,15)  5.256898E-02  2.456998E-01  2.336366E+00  9.128858E+00
        30 2022-07-01 00:05:00.00  8.434590E-01  2.724668E+02  2.733103E+02  2.233190E+12
                    (0017,259,15)  5.490788E-02  2.532559E-01  2.554436E+00  9.476287E+00
        31 2022-07-01 00:05:10.00  8.551569E-01  2.724681E+02  2.733232E+02  2.233256E+12
                    (0017,259,15)  5.723030E-02  2.600741E-01  2.809879E+00  9.797315E+00
        32 2022-07-01 00:05:20.00  8.679532E-01  2.724695E+02  2.733374E+02  2.233324E+12
                    (0017,259,15)  5.955923E-02  2.659011E-01  3.096651E+00  1.008078E+01
        33 2022-07-01 00:05:30.00  8.817100E-01  2.724709E+02  2.733526E+02  2.233394E+12
                    (0017,259,15)  6.195360E-02  2.705415E-01  3.401821E+00  1.031948E+01
        34 2022-07-01 00:05:40.00  8.960398E-01  2.724725E+02  2.733685E+02  2.233466E+12
                    (0017,259,15)  6.446446E-02  2.737924E-01  3.706295E+00  1.051170E+01
        35 2022-07-01 00:05:50.00  9.106499E-01  2.724742E+02  2.733848E+02  2.233541E+12
                    (0004,314,15)  5.550024E-05  6.221087E-03  4.607031E+00  1.076142E+01
        36 2022-07-01 00:06:00.00  9.255326E-01  2.724759E+02  2.734015E+02  2.233617E+12
                    (0017,259,15)  6.991315E-02  2.756163E-01  4.209913E+00  1.100777E+01
        37 2022-07-01 00:06:10.00  9.408221E-01  2.724778E+02  2.734186E+02  2.233696E+12
                    (0017,259,15)  7.280943E-02  2.741114E-01  4.351823E+00  1.121645E+01
        38 2022-07-01 00:06:20.00  9.566272E-01  2.724796E+02  2.734363E+02  2.233777E+12
                    (0037,385,15)  0.000000E+00  0.000000E+00  4.967595E+00  1.139552E+01
        39 2022-07-01 00:06:30.00  9.729726E-01  2.724816E+02  2.734546E+02  2.233860E+12
                    (0036,386,15)  1.806783E-04  5.987316E-04  5.557924E+00  1.155213E+01
        40 2022-07-01 00:06:40.00  9.898123E-01  2.724836E+02  2.734734E+02  2.233945E+12
                    (0036,385,15)  3.040020E-05  8.747489E-03  6.003715E+00  1.169009E+01
        41 2022-07-01 00:06:50.00  1.007059E+00  2.724857E+02  2.734928E+02  2.234033E+12
                    (0015,360,15)  6.174161E-02  2.607138E-03  5.401228E+00  1.192173E+01
        42 2022-07-01 00:07:00.00  1.024619E+00  2.724879E+02  2.735125E+02  2.234122E+12
                    (0037,383,15)  1.506375E-02  1.167762E-01  5.199165E+00  1.224020E+01
        43 2022-07-01 00:07:10.00  1.042424E+00  2.724901E+02  2.735325E+02  2.234214E+12
                    (0018,360,15)  3.648987E-03  1.922862E-02  5.984663E+00  1.254478E+01
        44 2022-07-01 00:07:20.00  1.060459E+00  2.724924E+02  2.735529E+02  2.234308E+12
                    (0016,359,15)  1.350745E-03  3.373350E-03  7.162504E+00  1.283524E+01
        45 2022-07-01 00:07:30.00  1.078710E+00  2.724948E+02  2.735735E+02  2.234404E+12
                    (0018,222,15)  5.247130E-02  1.741643E-01  7.893377E+00  1.313766E+01
        46 2022-07-01 00:07:40.00  1.097184E+00  2.724972E+02  2.735944E+02  2.234502E+12
                    (0253,430,15)  9.557020E-02  1.780536E-01  8.442327E+00  1.351012E+01
        47 2022-07-01 00:07:50.00  1.115915E+00  2.724997E+02  2.736156E+02  2.234602E+12
                    (0016,222,10)  1.616510E-02  2.159180E-01  8.633033E+00  1.407749E+01
        48 2022-07-01 00:08:00.00  1.134891E+00  2.725023E+02  2.736372E+02  2.234704E+12
                    (0094,423,15)  9.320101E-02  3.568212E-01  9.272713E+00  1.472533E+01
        49 2022-07-01 00:08:10.00  1.154084E+00  2.725049E+02  2.736590E+02  2.234809E+12
                    (0016,223,15)  9.813783E-02  3.834881E-01  1.093703E+01  1.539604E+01
        50 2022-07-01 00:08:20.00  1.173480E+00  2.725076E+02  2.736811E+02  2.234916E+12
                    (0013,224,15)  7.015166E-03  1.147940E-01  1.207139E+01  1.566378E+01
        51 2022-07-01 00:08:30.00  1.193076E+00  2.725104E+02  2.737034E+02  2.235024E+12
                    (0019,222,15)  1.317294E-01  4.226298E-01  1.314022E+01  1.568593E+01
        52 2022-07-01 00:08:40.00  1.212765E+00  2.725132E+02  2.737260E+02  2.235135E+12
                    (0096,423,15)  0.000000E+00  8.441780E-02  1.662590E+01  1.599557E+01
        53 2022-07-01 00:08:50.00  1.232608E+00  2.725161E+02  2.737487E+02  2.235248E+12
                    (0095,423,15)  0.000000E+00  1.205388E-01  1.671981E+01  1.648972E+01
        54 2022-07-01 00:09:00.00  1.252540E+00  2.725191E+02  2.737716E+02  2.235364E+12
                    (0018,222,15)  5.572647E-02  1.711980E-01  1.394650E+01  1.723011E+01
        55 2022-07-01 00:09:10.00  1.272582E+00  2.725221E+02  2.737947E+02  2.235481E+12
                    (0094,423,15)  6.317270E-02  1.915878E-01  1.703175E+01  1.760668E+01
        56 2022-07-01 00:09:20.00  1.292743E+00  2.725253E+02  2.738180E+02  2.235600E+12
                    (0094,423,15)  1.604990E-02  1.988026E-02  1.889611E+01  1.803640E+01
        57 2022-07-01 00:09:30.00  1.313068E+00  2.725285E+02  2.738415E+02  2.235722E+12
                    (0021,222,15)  1.652673E-01  3.108906E-01  1.240386E+01  1.870935E+01
        58 2022-07-01 00:09:40.00  1.333405E+00  2.725317E+02  2.738651E+02  2.235845E+12
                    (0091,424,15)  4.133576E-02  0.000000E+00  1.270326E+01  1.951319E+01
        59 2022-07-01 00:09:50.00  1.353876E+00  2.725351E+02  2.738889E+02  2.235971E+12
                    (0021,221,15)  1.420878E-01  1.864156E-01  1.682867E+01  1.843033E+01
        60 2022-07-01 00:10:00.00  1.374432E+00  2.725385E+02  2.739129E+02  2.236099E+12
                    (0022,222,15)  9.338590E-03  4.706796E-04  1.740711E+01  1.890892E+01
        61 2022-07-01 00:10:10.00  1.395013E+00  2.725420E+02  2.739370E+02  2.236229E+12
                    (0252,431,15)  2.691726E-01  3.350836E-01  1.470122E+01  1.726274E+01
        62 2022-07-01 00:10:20.00  1.415690E+00  2.725456E+02  2.739613E+02  2.236361E+12
                    (0018,222,15)  1.818468E-03  5.261659E-02  1.685289E+01  1.804296E+01
        63 2022-07-01 00:10:30.00  1.436418E+00  2.725493E+02  2.739857E+02  2.236495E+12
                    (0253,430,15)  2.671322E-01  3.846515E-01  1.538824E+01  1.752374E+01
        64 2022-07-01 00:10:40.00  1.457242E+00  2.725530E+02  2.740102E+02  2.236631E+12
                    (0253,430,15)  2.601757E-01  3.802022E-01  2.048496E+01  1.921298E+01
        65 2022-07-01 00:10:50.00  1.478147E+00  2.725568E+02  2.740349E+02  2.236770E+12
                    (0252,431,15)  3.192724E-01  3.781810E-01  2.060153E+01  1.842632E+01
        66 2022-07-01 00:11:00.00  1.499156E+00  2.725607E+02  2.740598E+02  2.236910E+12
                    (0254,429,15)  1.103315E-01  1.795927E-01  2.068201E+01  1.878401E+01
        67 2022-07-01 00:11:10.00  1.520237E+00  2.725646E+02  2.740849E+02  2.237053E+12
                    (0255,429,15)  8.168551E-02  1.647991E-01  1.839977E+01  1.909014E+01
        68 2022-07-01 00:11:20.00  1.541395E+00  2.725687E+02  2.741101E+02  2.237197E+12
                    (0253,430,15)  3.195095E-01  4.536736E-01  1.882806E+01  1.919326E+01
        69 2022-07-01 00:11:30.00  1.562617E+00  2.725728E+02  2.741354E+02  2.237344E+12
                    (0253,430,15)  1.270807E-01  1.329343E-01  2.396751E+01  1.941464E+01
        70 2022-07-01 00:11:40.00  1.583924E+00  2.725769E+02  2.741609E+02  2.237492E+12
                    (0260,430,15)  1.952529E-01  4.283824E-01  1.283693E+01  1.990810E+01
        71 2022-07-01 00:11:50.00  1.605323E+00  2.725812E+02  2.741865E+02  2.237643E+12
                    (0016,313,15)  1.863008E-02  6.176535E-01  1.817253E+01  2.083678E+01
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 321      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 298      Source: ROMS/Drivers/nl_roms.h, ROMS_run

 Blowing-up: Saving latest model state into  RESTART file
     REASON: MaxSpeed =  2.1E+01

  WRT_RST_NF90     - wrote re-start    fields (Index=1,2) in record = 1

Elapsed wall CPU time for each process (seconds):

 Thread #    0 CPU:     414.027
 Thread #    1 CPU:     414.027
 Thread #    3 CPU:     414.027
 Thread #    8 CPU:     414.027
 Thread #    9 CPU:     414.027
 Thread #   10 CPU:     414.027
 Thread #   12 CPU:     414.027
 Thread #   14 CPU:     414.027
 Thread #   13 CPU:     414.027
 Thread #    2 CPU:     414.041
 Thread #   11 CPU:     414.042
 Thread #   15 CPU:     414.042
 Thread #    6 CPU:     414.042
 Thread #    7 CPU:     414.042
 Thread #    4 CPU:     414.043
 Thread #    5 CPU:     414.043
 Total:                6624.538
 Average:               414.034

 Nonlinear model elapsed CPU time profile, Grid: 01

  Allocation and array initialization ..............       392.368  ( 5.9229 %)
  Ocean state initialization .......................        18.456  ( 0.2786 %)
  Reading of input data ............................         5.016  ( 0.0757 %)
  Processing of input data .........................        61.348  ( 0.9261 %)
  Processing of output time averaged data ..........        82.149  ( 1.2401 %)
  Computation of vertical boundary conditions ......         1.525  ( 0.0230 %)
  Computation of global information integrals ......        23.795  ( 0.3592 %)
  Writing of output data ...........................         2.120  ( 0.0320 %)
  Model 2D kernel ..................................      3392.867  (51.2167 %)
  2D/3D coupling, vertical metrics .................        43.469  ( 0.6562 %)
  Omega vertical velocity ..........................        30.227  ( 0.4563 %)
  Equation of state for seawater ...................        42.501  ( 0.6416 %)
  Atmosphere-Ocean bulk flux parameterization ......        49.126  ( 0.7416 %)
  GLS vertical mixing parameterization .............       608.398  ( 9.1840 %)
  3D equations right-side terms ....................       252.980  ( 3.8188 %)
  3D equations predictor step ......................       324.828  ( 4.9034 %)
  Pressure gradient ................................       135.643  ( 2.0476 %)
  Harmonic mixing of tracers, S-surfaces ...........       197.644  ( 2.9835 %)
  Harmonic stress tensor, S-surfaces ...............       150.079  ( 2.2655 %)
  Corrector time-step for 3D momentum ..............        68.419  ( 1.0328 %)
  Corrector time-step for tracers ..................       580.952  ( 8.7697 %)
                                              Total:      6463.908   97.5752 %

  Unique kernel(s) regions profiled ................      6463.908   97.5752 %
  Residual, non-profiled code ......................       160.630    2.4248 %


 All percentages are with respect to total time =         6624.538


>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

 Dynamic and Automatic memory (MB) usage for Grid 01:  1098x748x15  tiling: 4x4

     tile          Dynamic        Automatic            USAGE

        0          7522.36           211.20          7733.56
        1             0.00           211.20           211.20
        2             0.00           211.20           211.20
        3             0.00           211.20           211.20
        4             0.00           211.20           211.20
        5             0.00           211.20           211.20
        6             0.00           211.20           211.20
        7             0.00           211.20           211.20
        8             0.00           211.20           211.20
        9             0.00           211.20           211.20
       10             0.00           211.20           211.20
       11             0.00           211.20           211.20
       12             0.00           211.20           211.20
       13             0.00           211.20           211.20
       14             0.00           211.20           211.20
       15             0.00           211.20           211.20

    TOTAL          7522.36          3379.20         10901.56

<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

 ROMS/TOMS - Output NetCDF summary for Grid 01:
             number of time records written in HISTORY file = 1
             number of time records written in RESTART file = 1

 Analytical header files used:

     ROMS/Functionals/ana_btflux.h
     ROMS/Functionals/ana_stflux.h

 MAIN: Abnormal termination: BLOWUP.
 REASON: MaxSpeed =  2.1E+01

Zeng
Posts: 26
Joined: Mon Sep 19, 2022 1:06 pm
Location: Sun Yat-Sen University

Re: BLOWUP. REASON: MaxSpeed = 2.1E+01

#4 Unread post by Zeng »

Hello, all guys, I checked my forcing files, I found the high values in Uwind and Vwind. So I guess the time I started simulation was during a typhoon, and my simulation domain is quete small which may leads the blowup. Now I would change the start time and try it again, I'm not sure if this help, or could you give me some other suggestions? Thank you so much! :mrgreen:
u.png
v.png

User avatar
kate
Posts: 4091
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: BLOWUP. REASON: MaxSpeed = 2.1E+01

#5 Unread post by kate »

Maybe you need an even shorter timestep during a typhoon? The model does want a shorter timestep than you might think, especially with lots of vertical levels.

Zeng
Posts: 26
Joined: Mon Sep 19, 2022 1:06 pm
Location: Sun Yat-Sen University

Re: BLOWUP. REASON: MaxSpeed = 2.1E+01

#6 Unread post by Zeng »

Dear Kate,
Thank you for your reply, besh wishes to you! :D

My grid has 15 layers in vertical direction. Actually, I don't want to simulate the situation during a typhoon, it's accidental. So I changed my start time to a calm day which is helpful, and now, with a DT of 20s, the mode has run successfully for several hours.

Kate, May I ask you another question? If I want to use bathy_smoother in Pyroms, but I only find some functions in the directory. From my understanding, we should run an executable file to smooth the bethymetry. Or, my understanding is incorrect? How can I use the bathy_smoother? Thank you again!

User avatar
kate
Posts: 4091
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: BLOWUP. REASON: MaxSpeed = 2.1E+01

#7 Unread post by kate »

I'm afraid I'm not the expert on bathy_smoother, but my understanding is that there's more than one way to smooth bathymetry in there. I use an external program that filters bathymetry in a way to target only the roughest patches first, then I finish up with one pass of a Shapiro filter. There's no reason you couldn't do that in Python instead. There's also lpsolve, which solves the filtering in a different way, but can lead to bathymetries that ROMS doesn't like.

Or you can follow me in using MOM6 which doesn't need a bathymetry smoother.

Zeng
Posts: 26
Joined: Mon Sep 19, 2022 1:06 pm
Location: Sun Yat-Sen University

Re: BLOWUP. REASON: MaxSpeed = 2.1E+01

#8 Unread post by Zeng »

Thank you, Kate, more information I have got. My simulation is still running which indicates that there are no obvious issues with bathymetry. So I'd like to keep stationary for now.

zduvims
Posts: 5
Joined: Fri Feb 10, 2023 4:06 pm
Location: Virginia Institute of Marine Science

Re: BLOWUP. REASON: MaxSpeed = 2.1E+01

#9 Unread post by zduvims »

Zeng wrote: Wed Jun 28, 2023 8:48 am Hello, all guys, I checked my forcing files, I found the high values in Uwind and Vwind. So I guess the time I started simulation was during a typhoon, and my simulation domain is quete small which may leads the blowup. Now I would change the start time and try it again, I'm not sure if this help, or could you give me some other suggestions? Thank you so much! :mrgreen:
u.pngv.png
Hi Junyang, I would like to know what software you used to check the wind fields. Thanks!

Zeng
Posts: 26
Joined: Mon Sep 19, 2022 1:06 pm
Location: Sun Yat-Sen University

Re: BLOWUP. REASON: MaxSpeed = 2.1E+01

#10 Unread post by Zeng »

Hello, zduvims, sorry for my late reply. It's Panoply, the soft can quick read your .nc file directly.

Post Reply