An error occurred when running ROMS in COAWST

Discussion of how to use ROMS on different regional and basin scale applications.

Moderators: arango, robertson

Post Reply
Message
Author
s193264
Posts: 17
Joined: Wed Apr 06, 2022 5:59 pm
Location: College

An error occurred when running ROMS in COAWST

#1 Unread post by s193264 »

Hi everyone

I'm trying to run ROMS alone in the COAWST model.
The header file (.h) is set to #define ROMS, and the others are set to #undef.

coast.bash was able to compile successfully, and I prepared various necessary files and tried to execute the calculation.
Then an error occurred.
Below are my header files (.h), input files (.in) and output results (.out).

<.h file>

Code: Select all

/*
** svn $Id: sandy.h 25 2007-04-09 23:43:58Z jcwarner $
*******************************************************************************
** Copyright (c) 2002-2007 The ROMS/TOMS Group                               **
**   Licensed under a MIT/X style license                                    **
**   See License_ROMS.txt                                                    **
*******************************************************************************
**
** Options for Sandy Test.
**
** Application flag:   SANDY
*/

#define ROMS_MODEL
#undef NESTING
#undef WRF_MODEL
#undef SWAN_MODEL
#undef WW3_MODEL
#undef MCT_LIB
#undef MCT_INTERP_OC2AT
#undef MCT_INTERP_WV2AT
#undef MCT_INTERP_OC2WV

#if defined WRF_MODEL && (defined SWAN_MODEL || defined WW3_MODEL)
# define DRAGLIM_DAVIS
# define COARE_TAYLOR_YELLAND
#endif

#ifdef ROMS_MODEL
/* Physics + numerics */
# if defined WW3_MODEL || defined SWAN_MODEL
#  define WEC_VF
#  define WDISS_WAVEMOD
#  define UV_KIRBY
# endif
# define UV_ADV
# define UV_COR
# define UV_VIS2
# define MIX_S_UV
# undef  TS_FIXED

# undef SSW_BBL
# ifdef SSW_BBL
#  define SSW_CALC_ZNOT
#  define ANA_SEDIMENT
# else
#  define UV_LOGDRAG
# endif
# if !defined SWAN_MODEL && defined SSW_BBL
#  define ANA_WWAVE
# endif

# define DJ_GRADPS
# define TS_DIF2
# define MIX_GEO_TS
# define CURVGRID

# define SALINITY
# define SOLVE3D
# define SPLINES_VDIFF
# define SPLINES_VVISC
# define AVERAGES
# define NONLIN_EOS

/* Grid and Initial */
# define MASKING

/* Forcing */
# ifdef WRF_MODEL
#  undef  BULK_FLUXES
#  define ATM2OCN_FLUXES
#  undef  ANA_SSFLUX
#  undef  LONGWAVE_OUT
# else
#  define BULK_FLUXES
# endif
# define ANA_NUDGCOEF
# define ATM_PRESS
# define ANA_BTFLUX
# define ANA_BSFLUX
# define ANA_BPFLUX
# define ANA_SPFLUX
# define ANA_SRFLUX
# define EMINUSP
# define SOLAR_SOURCE

/* Turbulence closure */
# define GLS_MIXING
# undef  MY25_MIXING
# define LIMIT_VDIFF
# define LIMIT_VVISC

# if defined GLS_MIXING || defined MY25_MIXING
#  define KANTHA_CLAYSON
#  define N2S2_HORAVG
#  define RI_SPLINES
#  define CRAIG_BANNER
#  define CHARNOK
# endif

# define SSH_TIDES
# define UV_TIDES
# define RAMP_TIDES
# define ANA_FSOBC
# define ANA_M2OBC

/* Output */
# undef DIAGNOSTICS_UV
# undef DIAGNOSTICS_TS
#endif
<.out file>

Code: Select all

--------------------------------------------------------------------------------
 Model Input Parameters:  ROMS/TOMS version 3.9  
                          Friday - July 29, 2022 - 11:09:47 PM
--------------------------------------------------------------------------------

 Hurricane Katrina

 Operating system  : Linux
 CPU/hardware      : x86_64
 Compiler system   : gfortran
 Compiler command  : /home/concretelab/Build_COAWST/LIBRARIES/mpich/bin/mpif90
 Compiler flags    : -frepack-arrays -O3 -ffast-math -O3 -ftree-vectorize -ftree-loop-linear -funroll-loops -w -
 OCN Communicator  : 1140850688,  PET size = 1

 Input Script      : 

 SVN Root URL     : 
 SVN Revision     : 

 Local Root       : /home/concretelab/Build_COAWST/ROMS_COAWST
 Header Dir       : /home/concretelab/Build_COAWST/ROMS_COAWST/Projects/Katrina
 Header file      : katrina.h
 Analytical Dir   : /home/concretelab/Build_COAWST/ROMS_COAWST/Projects/Katrina

 Resolution, Grid 01: 98x98x15,  Parallel Nodes: 1,  Tiling: 1x1


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

      10800  ntimes            Number of timesteps for 3-D equations.
     40.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.
         60  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).
         60  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.
        432  nAVG              Number of timesteps between the writing of
                                 time-averaged data into averages file.
 2.0000E-01  nl_tnu2(01)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 01: temp
 2.0000E-01  nl_tnu2(02)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 02: salt
 1.0000E-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
 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-03  Akt_limit(01)     Vertical diffusion upper threshold (m2/s)
                                 for tracer 01: temp
 1.0000E-03  Akt_limit(02)     Vertical diffusion upper threshold (m2/s)
                                 for tracer 02: salt
 1.0000E-05  Akv_bak           Background vertical mixing coefficient (m2/s)
                                 for momentum.
 1.0000E-03  Akv_limit         Vertical viscosity upper threshold (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.
      3.000  gls_p             GLS stability exponent.
      1.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.
 1.4400E+00  gls_c1            GLS shear production coefficient.
 1.9200E+00  gls_c2            GLS dissipation coefficient.
-4.0000E-01  gls_c3m           GLS stable buoyancy production coefficient.
 1.0000E+00  gls_c3p           GLS unstable buoyancy production coefficient.
 1.0000E+00  gls_sigk          GLS constant Schmidt number for TKE.
 1.3000E+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.
      0.000  wec_alpha         WEC factor for roller/breaking energy distribution.
 3.0000E-04  rdrg              Linear bottom drag coefficient (m/s).
 2.5000E-02  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.
 1.0000E+01  blk_ZW            Height (m) of surface winds measurement.
          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.
  38591.000  dstart            Time-stamp assigned to model initialization (days).
  38591.000  tide_start        Reference time origin for tidal forcing (days).
18581117.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
 0.0000E+00  Znudg             Nudging/relaxation time scale (days)
                                 for free-surface.
 0.0000E+00  M2nudg            Nudging/relaxation time scale (days)
                                 for 2D momentum.
 1.0000E+00  M3nudg            Nudging/relaxation time scale (days)
                                 for 3D momentum.
 0.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.
     10.000  T0                Background potential temperature (C) constant.
     30.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  LsshCLM           Turning OFF processing of SSH climatology.
          F  Lm2CLM            Turning OFF processing of 2D momentum climatology.
          F  Lm3CLM            Turning OFF 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.
          F  LnudgeM3CLM       Turning OFF 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(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(idUsms)      Write out surface U-momentum stress.
          T  Hout(idVsms)      Write out surface V-momentum stress.
          T  Hout(idUbms)      Write out bottom U-momentum stress.
          T  Hout(idVbms)      Write out bottom V-momentum stress.
          T  Hout(idPair)      Write out surface air pressure.
          T  Hout(idTair)      Write out surface air temperature.
          T  Hout(idUair)      Write out surface U-wind component.
          T  Hout(idVair)      Write out surface V-wind component.
          T  Hout(idUairE)     Write out surface Eastward U-wind component.
          T  Hout(idVairN)     Write out surface Northward V-wind component.
          T  Hout(idTsur)      Write out surface net heat flux.
          T  Hout(idTsur)      Write out surface net salt flux.
          T  Hout(idSrad)      Write out shortwave radiation flux.
          T  Hout(idLrad)      Write out longwave radiation flux.
          T  Hout(idLhea)      Write out latent heat flux.
          T  Hout(idShea)      Write out sensible heat flux.
          T  Hout(idEmPf)      Write out E-P flux.
          T  Hout(idevap)      Write out evaporation rate.
          T  Hout(idrain)      Write out rain rate.
          T  Hout(idVvis)      Write out vertical viscosity: AKv.
          T  Hout(idMtke)      Write out turbulent kinetic energy.
          T  Hout(idMtls)      Write out turbulent generic length-scale.

          T  Aout(idSdif)      Write out averaged vertical diffusion: AKt(isalt).

 Output/Input Files:

               Output Restart File:  Katrina_ocean_rst.nc
               Output History File:  Katrina_ocean_his.nc
              Output Averages File:  Katrina_ocean_avg.nc
                   Input Grid File:  Projects/Katrina/roms_grid_Katrina.nc
      Input Nonlinear Initial File:  Projects/Katrina/Katrina_ini.nc
                Tidal Forcing File:  Projects/Katrina/tide_forc_Katrina.nc
             Input Forcing File 01:  Projects/Katrina/romsforc_NARR_Katrina2005.nc
         Input Climatology File 01:  Projects/Katrina/Katrina_clim.nc
    Input Lateral Boundary File 01:  Projects/Katrina/Katrina_bry.nc
  ROMS I/O variables Metadata File:  ROMS/External/varinfo.dat

 Tile partition information for Grid 01:  98x98x15  tiling: 1x1

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1       98        1       98   144060

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

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50    98.50     0.50    98.50  RHO-points

        0     1.00    98.00     0.50    98.50    U-points

        0     0.50    98.50     1.00    98.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     218
               HaloSizeJ(1) =     218
                TileSide(1) =     103
                TileSize(1) =   10609


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

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

 temp                     1   Upstream3          Centered4

 salt                     1   Upstream3          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  Chapman Imp

 ubar                     1   Closed       Flather      Flather      Flather

 vbar                     1   Closed       Flather      Flather      Flather

 u                        1   Closed       Nested       Nested       Nested

 v                        1   Closed       Nested       Nested       Nested

 temp                     1   Gradient     Nested       Nested       Gradient

 salt                     1   Gradient     Nested       Nested       Gradient

 tke                      1   Closed       Gradient     Gradient     Gradient

 Activated C-preprocessing Options:

 KATRINA                  Hurricane Katrina
 ANA_BSFLUX               Analytical kinematic bottom salinity flux
 ANA_BTFLUX               Analytical kinematic bottom temperature flux
 ANA_FSOBC                Analytical free-surface boundary conditions
 ANA_M2OBC                Analytical 2D momentum boundary conditions
 ANA_NUDGCOEF             Analytical spatially varying nudging time-scales
 ANA_SRFLUX               Analytical kinematic shortwave radiation flux
 ASSUMED_SHAPE            Using assumed-shape arrays
 ATM_PRESS                Impose atmospheric pressure onto sea surface
 AVERAGES                 Writing out time-averaged nonlinear model fields
 !BOUNDARY_ALLGATHER      Using mpi_allreduce in mp_boundary routine
 BULK_FLUXES              Surface bulk fluxes parameterization
 !COLLECT_ALL...          Using mpi_isend/mpi_recv in mp_collect routine
 CHARNOK                  Charnok surface roughness from wind stress
 CRAIG_BANNER             Craig and Banner wave breaking surface flux
 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
 KANTHA_CLAYSON           Kantha and Clayson stability function formulation
 LIMIT_VDIFF              Impose an upper limit on vertical diffusion coefficient
 LIMIT_VVISC              Impose an upper limit on vertical viscosity coefficient
 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
 N2S2_HORAVG              Horizontal smoothing of buoyancy and shear
 POWER_LAW                Power-law shape time-averaging barotropic filter
 PROFILE                  Time profiling activated
 K_GSCHEME                Third-order upstream advection of TKE fields
 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              Double 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
 SSH_TIDES                Add tidal elevation to SSH climatology
 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_LOGDRAG               Logarithmic bottom stress
 UV_TIDES                 Add tidal currents to 2D momentum climatologies
 UV_VIS2                  Harmonic mixing of momentum
 VAR_RHO_2D               Variable density barotropic mode

 Process Information:

 Node #    0 (pid=  105506) is active.

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

    GET_GRID    - angle between XI-axis and EAST: angler
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -2.10128491E-12 Max =  2.88440987E-12)
    GET_GRID    - bathymetry at RHO-points: h
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -2.15744419E+02 Max =  3.33282292E+03)
    GET_GRID    - Coriolis parameter at RHO-points: f
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  6.41047590E-05 Max =  7.91128633E-05)
    GET_GRID    - reciprocal XI-grid spacing: pm
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  1.40936455E-04 Max =  1.50689689E-04)
    GET_GRID    - reciprocal ETA-grid spacing: pn
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  1.31388388E-04 Max =  1.31388389E-04)
    GET_GRID    - XI-derivative of inverse metric factor pn: dndx
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -5.92140623E-07 Max =  5.92140623E-07)
    GET_GRID    - ETA-derivative of inverse metric factor pm: dmde
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -5.10947184E+00 Max =  0.00000000E+00)
    GET_GRID    - longitude of RHO-points: lon_rho
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30376282E+01 Max = -8.60048832E+01)
    GET_GRID    - longitude of U-points: lon_u
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30021093E+01 Max = -8.60404021E+01)
    GET_GRID    - longitude of V-points: lon_v
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30376282E+01 Max = -8.60048832E+01)
    GET_GRID    - longitude of PSI-points: lon_psi
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30021093E+01 Max = -8.60404021E+01)
    GET_GRID    - latitude of RHO-points lat_rho
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.60755171E+01 Max =  3.28514916E+01)
    GET_GRID    - latitude of U-points: lat_u
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.60755171E+01 Max =  3.28514916E+01)
    GET_GRID    - latitude of V-points: lat_v
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.61097392E+01 Max =  3.28172695E+01)
    GET_GRID    - latitude of PSI-points lat_psi
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.61097392E+01 Max =  3.28172695E+01)
    GET_GRID    - mask on RHO-points: mask_rho
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on U-points: mask_u
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on V-points: mask_v
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on PSI-points: mask_psi
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.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       -1.375      -0.671      -1.929      -2.645
    13  -0.1333333  -0.0017081       -2.533      -1.350      -5.261      -8.310
    12  -0.2000000  -0.0043062       -3.385      -2.043     -10.576     -18.242
    11  -0.2666667  -0.0089353       -3.754      -2.756     -19.015     -34.904
    10  -0.3333333  -0.0169029       -3.329      -3.502     -32.592     -62.625
     9  -0.4000000  -0.0304126       -1.586      -4.304     -54.697    -108.707
     8  -0.4666667  -0.0530642        2.331      -5.197     -90.870    -185.076
     7  -0.5333333  -0.0905094        9.766      -6.238    -149.807    -310.455
     6  -0.6000000  -0.1510347       22.688      -7.510    -244.259    -512.298
     5  -0.6666667  -0.2452550       43.623      -9.119    -390.561    -825.770
     4  -0.7333333  -0.3828872       74.881     -11.162    -603.665   -1283.062
     3  -0.8000000  -0.5634305      116.342     -13.634    -882.800   -1882.517
     2  -0.8666667  -0.7610631      161.867     -16.277   -1188.232   -2538.588
     1  -0.9333333  -0.9210683      198.445     -18.544   -1435.763   -3070.001
     0  -1.0000000  -1.0000000      215.744     -20.000   -1558.539   -3332.823

 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 =  6.63615411E+00 km    Water points =  6.79507250E+00 km
 Maximum X-grid spacing, DXmax =  7.09539627E+00 km    Water points =  7.09539627E+00 km
 Minimum Y-grid spacing, DYmin =  7.61102264E+00 km    Water points =  7.61102264E+00 km
 Maximum Y-grid spacing, DYmax =  7.61102265E+00 km    Water points =  7.61102265E+00 km
 Minimum Z-grid spacing, DZmin = -8.41525285E+01 m     Water points = -1.38040192E-01 m
 Maximum Z-grid spacing, DZmax =  4.73552371E+01 m     Water points =  6.56070675E+02 m

 Minimum barotropic Courant Number =  7.42536615E-07
 Maximum barotropic Courant Number =  4.64534958E-02
 Maximum Coriolis   Courant Number =  2.97491249E-03


 NLM: GET_STATE - Reading state initial conditions,                       2005-08-27 00:00:00.00
                   (Grid 01, t = 53609.0000, File: Katrina_ini.nc, Rec=0001, Index=1)
                - free-surface
                   (Min = -3.54594313E-01 Max =  4.45445436E-01)
                - vertically integrated u-momentum component
                   (Min = -1.41692584E-01 Max =  2.49862346E-01)
                - vertically integrated v-momentum component
                   (Min = -1.73684799E-01 Max =  1.14896163E-01)
                - u-momentum component
                   (Min = -9.26025459E-01 Max =  1.03917037E+00)
                - v-momentum component
                   (Min = -9.95687677E-01 Max =  5.99752986E-01)
                - potential temperature
                   (Min =  4.15533894E+00 Max =  3.54882013E+01)
                - salinity
                   (Min = -1.15116346E+01 Max =  3.68338484E+01)
    GET_NGFLD   - tidal period
                   (Grid = 01, Min =  1.49047215E+04 Max =  9.67260840E+04)
    GET_2DFLD   - tidal elevation amplitude
                   (Grid = 01, Min =  0.00000000E+00 Max =  2.09046983E-01)
    GET_2DFLD   - tidal elevation phase angle
                   (Grid = 01, Min = -1.29687416E+01 Max =  4.82525186E+00)
    GET_2DFLD   - tidal current inclination angle
                   (Grid = 01, Min =  0.00000000E+00 Max =  6.28316919E+00)
    GET_2DFLD   - tidal current phase angle
                   (Grid = 01, Min =  0.00000000E+00 Max =  6.28306597E+00)
    GET_2DFLD   - maximum tidal current, ellipse major axis
                   (Grid = 01, Min =  0.00000000E+00 Max =  3.46792927E-01)
    GET_2DFLD   - minimum tidal current, ellipse minor axis
                   (Grid = 01, Min = -8.07614703E-02 Max =  4.40275504E-02)
    GET_2DFLD   - surface u-wind component,                        2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -1.07696928E+01 Max = -8.92840971E-01)      regrid = T
    GET_2DFLD   - surface v-wind component,                        2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -8.16268255E+00 Max =  3.43387192E+00)      regrid = T
    GET_2DFLD   - surface air pressure,                            2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  1.00612364E+03 Max =  1.01279449E+03)      regrid = T
    GET_2DFLD   - net longwave radiation flux,                     2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -1.29572949E-05 Max = -6.30481076E-07)      regrid = T
    GET_2DFLD   - surface air temperature,                         2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  2.74625011E+01 Max =  3.52368334E+01)      regrid = T
    GET_2DFLD   - surface air relative humidity,                   2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  4.85722983E-01 Max =  8.13929408E-01)      regrid = T
    GET_2DFLD   - rain fall rate,                                  2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  3.83558273E-07 Max =  3.75354200E-05)      regrid = T
    GET_NGFLD   - 3D u-momentum eastern boundary condition,        2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -3.29027653E-01 Max =  2.30727926E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,        2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -8.52026403E-01 Max =  1.83875933E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -9.26025450E-01 Max =  6.07419252E-01)
    GET_NGFLD   - 3D v-momentum southern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -9.95687664E-01 Max =  5.74938476E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D v-momentum northern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - temperature eastern boundary condition,          2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  3.17894020E+01)
    GET_NGFLD   - salinity eastern boundary condition,             2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  3.63082657E+01)
    GET_NGFLD   - temperature southern boundary condition,         2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  4.22060776E+00 Max =  3.13313637E+01)
    GET_NGFLD   - salinity southern boundary condition,            2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  3.48331337E+01 Max =  3.68266792E+01)
    GET_3DFLD   - potential temperature climatology,               2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  4.15533894E+00 Max =  3.54882013E+01)
    GET_3DFLD   - salinity climatology,                            2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -1.15116346E+01 Max =  3.68338484E+01)

 Basin information for Grid 01:

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

 Initial domain volumes:  TotVolume =  3.7806098739E+14 m3
                         MinCellVol = -7.1542085900E+06 m3
                         MaxCellVol =  3.4742261021E+10 m3
                            Max/Min = -4.8561990587E+03


 NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000032438881 - 000032449680)

    GET_2DFLD   - surface u-wind component,                        2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min = -1.20530137E+01 Max =  5.51073609E-01)      regrid = T
    GET_2DFLD   - surface v-wind component,                        2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min = -7.55569400E+00 Max =  3.14627136E+00)      regrid = T
    GET_2DFLD   - surface air pressure,                            2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  1.00585813E+03 Max =  1.01232250E+03)      regrid = T
    GET_2DFLD   - net longwave radiation flux,                     2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min = -1.34115722E-05 Max = -1.42211893E-06)      regrid = T
    GET_2DFLD   - surface air temperature,                         2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  2.38640230E+01 Max =  3.15746287E+01)      regrid = T
    GET_2DFLD   - surface air relative humidity,                   2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  6.08187675E-01 Max =  8.83000289E-01)      regrid = T
    GET_2DFLD   - rain fall rate,                                  2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  2.56529711E-06 Max =  3.64055637E-04)      regrid = T
    GET_NGFLD   - 3D u-momentum eastern boundary condition,        2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -1.06390405E+00 Max =  2.70750403E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,        2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -4.67574030E-01 Max =  3.27160418E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -1.06390405E+00 Max =  5.65890431E-01)
    GET_NGFLD   - 3D v-momentum southern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -9.68096316E-01 Max =  6.13836169E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D v-momentum northern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - temperature eastern boundary condition,          2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  3.17175808E+01)
    GET_NGFLD   - salinity eastern boundary condition,             2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  3.63111801E+01)
    GET_NGFLD   - temperature southern boundary condition,         2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  4.21079731E+00 Max =  3.10041027E+01)
    GET_NGFLD   - salinity southern boundary condition,            2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  3.48317986E+01 Max =  3.66650505E+01)
    GET_3DFLD   - potential temperature climatology,               2005-08-28 00:00:00.00
                   (Grid=01, Rec=2, Index=2, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  4.19558276E+00 Max =  3.47792665E+01)
    GET_3DFLD   - salinity climatology,                            2005-08-28 00:00:00.00
                   (Grid=01, Rec=2, Index=2, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -1.77656485E+01 Max =  3.68394868E+01)

 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

  32438880 2005-08-27 00:00:00.00           NaN           NaN           NaN           NaN
                       (54,22,12)  5.507628E-03  1.540193E-03  0.000000E+00  1.049467E+00
 Found Error: 01   Line: 341      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 331      Source: ROMS/Drivers/nl_ocean.h, ROMS_run

Elapsed wall CPU time for each process (seconds):

 Node   #    0 CPU:       0.184
 Total:                   0.184

 Nonlinear model elapsed CPU time profile, Grid: 01

  Allocation and array initialization ..............         0.031  (16.7856 %)
  Ocean state initialization .......................         0.003  ( 1.7338 %)
  Reading of input data ............................         0.098  (53.5327 %)
  Processing of input data .........................         0.001  ( 0.4322 %)
  Computation of global information integrals ......         0.001  ( 0.6048 %)
  2D/3D coupling, vertical metrics .................         0.003  ( 1.4888 %)
  Omega vertical velocity ..........................         0.001  ( 0.2842 %)
  Equation of state for seawater ...................         0.006  ( 3.5079 %)
                                              Total:         0.144   78.3701 %

  Unique kernel(s) regions profiled ................         0.144   78.3701 %
  Residual, non-profiled code ......................         0.040   21.6299 %


 All percentages are with respect to total time =            0.184


 MPI communications profile, Grid: 01

  Message Passage: 2D halo exchanges ...............         0.000  ( 0.0240 %)
  Message Passage: 3D halo exchanges ...............         0.000  ( 0.0120 %)
  Message Passage: 4D halo exchanges ...............         0.000  ( 0.0016 %)
  Message Passage: data broadcast ..................         0.002  ( 0.9532 %)
  Message Passage: data reduction ..................         0.000  ( 0.0343 %)
  Message Passage: data scattering..................         0.007  ( 3.5776 %)
  Message Passage: point data gathering ............         0.000  ( 0.0011 %)
  Message Passage: synchronization barrier .........         0.000  ( 0.0022 %)
                                              Total:         0.008    4.6059 %

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

 Dynamic and Automatic memory (MB) usage for Grid 01:  98x98x15  tiling: 1x1

     tile          Dynamic        Automatic            USAGE      MPI-Buffers

        0            79.38             9.74            89.13             1.20

    TOTAL            79.38             9.74            89.13             1.20

<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
 Found Error: 01   Line: 55       Source: ROMS/Utility/close_io.F, close_inp.F

 ROMS/TOMS - Output NetCDF summary for Grid 01:

 Analytical header files used:

     ROMS/Functionals/ana_btflux.h
     ROMS/Functionals/ana_fsobc.h
     ROMS/Functionals/ana_m2obc.h
     ROMS/Functionals/ana_nudgcoef.h
     ROMS/Functionals/ana_srflux.h
 Found Error: 01   Line: 437      Source: ROMS/Utility/close_io.F, close_out

 ROMS/TOMS - Blows up ................ exit_flag: 1
Any advice is fine.
Please help me.
Thank you

jcwarner
Posts: 1098
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: An error occurred when running ROMS in COAWST

#2 Unread post by jcwarner »

GET_GRID - bathymetry at RHO-points: h
(Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
(Min = -2.15744419E+02 Max = 3.33282292E+03)

you need to activate wet dry if you have negative h values.
Suggest you take a look at your grid and decide if you are really resolving wet/dry processes.
if not, then maybe set a minimum depth of 2 m or something like that.
-j

s193264
Posts: 17
Joined: Wed Apr 06, 2022 5:59 pm
Location: College

Re: An error occurred when running ROMS in COAWST

#3 Unread post by s193264 »

jcwarner wrote: Fri Jul 29, 2022 2:58 pm GET_GRID - bathymetry at RHO-points: h
(Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
(Min = -2.15744419E+02 Max = 3.33282292E+03)

you need to activate wet dry if you have negative h values.
Suggest you take a look at your grid and decide if you are really resolving wet/dry processes.
if not, then maybe set a minimum depth of 2 m or something like that.
-j
Dear jcwarner

Thank you for your reply.
As per your advice, I changed the minimum value to 1m and ran it again.
Then the following error occurred.

Code: Select all

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
  43251840 2005-08-27 00:00:00.00  1.244784E-02  1.152586E+04  1.152587E+04  3.786342E+14
                       (54,22,12)  4.130721E-03  1.155144E-03  0.000000E+00  1.049467E+00
      DEF_HIS     - creating  history      file, Grid 01: Katrina_ocean_his.nc
      WRT_HIS     - wrote history     fields (Index=1,1) in record = 1
      DEF_AVG     - creating  average      file, Grid 01: Katrina_ocean_avg.nc
  43251841 2005-08-27 00:00:30.00  1.235159E-02  1.152582E+04  1.152583E+04  3.786331E+14
                       (58,55,11)  2.053789E-04  6.708336E-05  7.951317E-02  5.845875E+00
  43251842 2005-08-27 00:01:00.00           NaN           NaN           NaN           NaN
                       (92,18,15)  4.451173E-04      Infinity  9.036695E-03  0.000000E+00
 Found Error: 01   Line: 341      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 331      Source: ROMS/Drivers/nl_ocean.h, ROMS_run
 Blowing-up: Saving latest model state into  RESTART file
     REASON: KEchar =      NaN, PEchar =      NaN
      WRT_RST     - wrote re-start    fields (Index=1,1) in record = 1
I thought that this was because the time step was large,
so I tried setting dt=20 and dtfast=18, but this time I got a different error.

Code: Select all

 Process Information:

 Node #    0 (pid=  129799) is active.

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

    GET_GRID    - angle between XI-axis and EAST: angler
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -2.04122474E-12 Max =  2.76144271E-12)
    GET_GRID    - bathymetry at RHO-points: h
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  5.95202990E-03 Max =  3.35037820E+03)
    GET_GRID    - Coriolis parameter at RHO-points: f
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  6.41047132E-05 Max =  7.91128070E-05)
    GET_GRID    - reciprocal XI-grid spacing: pm
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  1.40936443E-04 Max =  1.50689667E-04)
    GET_GRID    - reciprocal ETA-grid spacing: pn
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  1.31369178E-04 Max =  1.31413711E-04)
    GET_GRID    - XI-derivative of inverse metric factor pn: dndx
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -5.92198376E-07 Max =  5.92198376E-07)
    GET_GRID    - ETA-derivative of inverse metric factor pm: dmde
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -5.10897846E+00 Max =  0.00000000E+00)
    GET_GRID    - longitude of RHO-points: lon_rho
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30376282E+01 Max = -8.60048832E+01)
    GET_GRID    - longitude of U-points: lon_u
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30021093E+01 Max = -8.60404021E+01)
    GET_GRID    - longitude of V-points: lon_v
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30376282E+01 Max = -8.60048832E+01)
    GET_GRID    - longitude of PSI-points: lon_psi
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30021093E+01 Max = -8.60404021E+01)
    GET_GRID    - latitude of RHO-points lat_rho
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.60754971E+01 Max =  3.28514653E+01)
    GET_GRID    - latitude of U-points: lat_u
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.60754971E+01 Max =  3.28514653E+01)
    GET_GRID    - latitude of V-points: lat_v
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.61097280E+01 Max =  3.28172547E+01)
    GET_GRID    - latitude of PSI-points lat_psi
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.61097280E+01 Max =  3.28172547E+01)
    GET_GRID    - mask on RHO-points: mask_rho
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on U-points: mask_u
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on V-points: mask_v
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on PSI-points: mask_psi
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.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.000      -0.671      -1.977      -2.652
    13  -0.1333333  -0.0017081       -0.001      -1.350      -5.463      -8.340
    12  -0.2000000  -0.0043062       -0.001      -2.043     -11.081     -18.318
    11  -0.2666667  -0.0089353       -0.002      -2.756     -20.062     -35.061
    10  -0.3333333  -0.0169029       -0.002      -3.502     -34.570     -62.922
     9  -0.4000000  -0.0304126       -0.002      -4.304     -58.251    -109.242
     8  -0.4666667  -0.0530642       -0.003      -5.197     -97.067    -186.008
     7  -0.5333333  -0.0905094       -0.003      -6.238    -160.373    -312.044
     6  -0.6000000  -0.1510347       -0.004      -7.510    -261.885    -514.949
     5  -0.6666667  -0.2452550       -0.004      -9.119    -419.178    -830.075
     4  -0.7333333  -0.3828872       -0.004     -11.162    -648.336   -1289.784
     3  -0.8000000  -0.5634305       -0.005     -13.634    -948.530   -1892.409
     2  -0.8666667  -0.7610631       -0.005     -16.277   -1277.014   -2551.949
     1  -0.9333333  -0.9210683       -0.006     -18.544   -1543.209   -3086.171
     0  -1.0000000  -1.0000000       -0.006     -20.000   -1675.192   -3350.378

 Time Splitting Weights for Grid 01:    ndtfast =  18    nfast =  26
 ==================================

    Primary            Secondary            Accumulated to Current Step

  1-0.0009217619993218 0.0555555555555555-0.0009217619993218 0.0555555555555555
  2-0.0010937935934335 0.0556067645555179-0.0020155555927553 0.1111623201110734
  3-0.0005164538239966 0.0556675308662642-0.0025320094167519 0.1668298509773376
  4 0.0008088610358613 0.0556962227453751-0.0017231483808905 0.2225260737227127
  5 0.0028785206760085 0.0556512860211606 0.0011553722951179 0.2781773597438733
  6 0.0056849852215555 0.0554913682058268 0.0068403575166734 0.3336687279497000
  7 0.0092146509828875 0.0551755356935181 0.0160550084995610 0.3888442636432182
  8 0.0134452174834803 0.0546636106389133 0.0295002259830412 0.4435078742821315
  9 0.0183425757655015 0.0539166541120533 0.0478428017485428 0.4974245283941847
 10 0.0238572179731971 0.0528976221250809 0.0717000197217398 0.5503221505192657
 11 0.0299201682140617 0.0515722211265700 0.1016201879358016 0.6018943716458357
 12 0.0364384346977950 0.0499099895591221 0.1380586226335966 0.6518043612049578
 13 0.0432899831530415 0.0478856320759113 0.1813486057866380 0.6996899932808691
 14 0.0503182315219162 0.0454806330118534 0.2316668373085543 0.7451706262927226
 15 0.0573260659323152 0.0426851757050803 0.2889929032408694 0.7878558019978029
 16 0.0640693779480095 0.0395003942643961 0.3530622811888789 0.8273561962621990
 17 0.0702501230965260 0.0359409843783956 0.4233124042854050 0.8632971806405947
 18 0.0755089006748113 0.0320381997619219 0.4988213049602163 0.8953353804025166
 19 0.0794170548326810 0.0278432608355435 0.5782383597928972 0.9231786412380601
 20 0.0814682969340545 0.0234312022337279 0.6597066567269516 0.9466098434717881
 21 0.0810698491959736 0.0189051857373916 0.7407765059229252 0.9655150292091796
 22 0.0775331096054062 0.0144013052265041 0.8183096155283314 0.9799163344356838
 23 0.0700638381138352 0.0100939102484260 0.8883734536421666 0.9900102446841098
 24 0.0577518641096317 0.0062014747976574 0.9461253177517983 0.9962117194817672
 25 0.0395603151682136 0.0029930379026779 0.9856856329200119 0.9992047573844451
 26 0.0143143670799881 0.0007952426155549 0.9999999999999999 1.0000000000000000

 ndtfast, nfast =   18  26   nfast/ndtfast =  1.44444

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

    1.000000000000 1.064325645941 0.532162822970 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.12622

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

 Minimum X-grid spacing, DXmin =  6.63615510E+00 km    Water points =  6.79507250E+00 km
 Maximum X-grid spacing, DXmax =  7.09539688E+00 km    Water points =  7.09539688E+00 km
 Minimum Y-grid spacing, DYmin =  7.60955604E+00 km    Water points =  7.61102264E+00 km
 Maximum Y-grid spacing, DYmax =  7.61213562E+00 km    Water points =  7.61213562E+00 km
 Minimum Z-grid spacing, DZmin =  3.96684645E-04 m     Water points =  3.96684645E-04 m
 Maximum Z-grid spacing, DZmax =  6.59540089E+02 m     Water points =  6.59540089E+02 m

 Minimum barotropic Courant Number =  5.26636174E-05
 Maximum barotropic Courant Number =  3.88104259E-02
 Maximum Coriolis   Courant Number =  1.48745624E-03


 NLM: GET_STATE - Reading state initial conditions,                       2005-08-27 00:00:00.00
                   (Grid 01, t = 53609.0000, File: Katrina_ini.nc, Rec=0001, Index=1)
                - free-surface
                   (Min = -3.54594313E-01 Max =  4.45445436E-01)
                - vertically integrated u-momentum component
                   (Min = -1.41692584E-01 Max =  2.49862346E-01)
                - vertically integrated v-momentum component
                   (Min = -1.73684799E-01 Max =  1.14896163E-01)
                - u-momentum component
                   (Min = -9.26025459E-01 Max =  1.03917037E+00)
                - v-momentum component
                   (Min = -9.95687677E-01 Max =  5.99752986E-01)
                - potential temperature
                   (Min =  4.15533894E+00 Max =  3.54882013E+01)
                - salinity
                   (Min = -1.15116346E+01 Max =  3.68338484E+01)
    GET_NGFLD   - tidal period
                   (Grid = 01, Min =  1.49047215E+04 Max =  9.67260840E+04)
    GET_2DFLD   - tidal elevation amplitude
                   (Grid = 01, Min =  0.00000000E+00 Max =  2.09046983E-01)
    GET_2DFLD   - tidal elevation phase angle
                   (Grid = 01, Min = -1.29687416E+01 Max =  4.82525186E+00)
    GET_2DFLD   - tidal current inclination angle
                   (Grid = 01, Min =  0.00000000E+00 Max =  6.28316919E+00)
    GET_2DFLD   - tidal current phase angle
                   (Grid = 01, Min =  0.00000000E+00 Max =  6.28306597E+00)
    GET_2DFLD   - maximum tidal current, ellipse major axis
                   (Grid = 01, Min =  0.00000000E+00 Max =  3.46792927E-01)
    GET_2DFLD   - minimum tidal current, ellipse minor axis
                   (Grid = 01, Min = -8.07614703E-02 Max =  4.40275504E-02)
    GET_2DFLD   - surface u-wind component,                        2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -1.07696928E+01 Max = -8.92840971E-01)      regrid = T
    GET_2DFLD   - surface v-wind component,                        2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -8.16269459E+00 Max =  3.43376645E+00)      regrid = T
    GET_2DFLD   - surface air pressure,                            2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  1.00612362E+03 Max =  1.01279447E+03)      regrid = T
    GET_2DFLD   - net longwave radiation flux,                     2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -1.29572887E-05 Max = -6.30481076E-07)      regrid = T
    GET_2DFLD   - surface air temperature,                         2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  2.74624980E+01 Max =  3.52368103E+01)      regrid = T
    GET_2DFLD   - surface air relative humidity,                   2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  4.85723022E-01 Max =  8.13929408E-01)      regrid = T
    GET_2DFLD   - rain fall rate,                                  2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  3.83558273E-07 Max =  3.75390348E-05)      regrid = T
    GET_NGFLD   - 3D u-momentum eastern boundary condition,        2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -3.29027653E-01 Max =  2.30727926E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,        2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -8.52026403E-01 Max =  1.83875933E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -9.26025450E-01 Max =  6.07419252E-01)
    GET_NGFLD   - 3D v-momentum southern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -9.95687664E-01 Max =  5.74938476E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D v-momentum northern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - temperature eastern boundary condition,          2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  3.17894020E+01)
    GET_NGFLD   - salinity eastern boundary condition,             2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  3.63082657E+01)
    GET_NGFLD   - temperature southern boundary condition,         2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  4.22060776E+00 Max =  3.13313637E+01)
    GET_NGFLD   - salinity southern boundary condition,            2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  3.48331337E+01 Max =  3.68266792E+01)
    GET_3DFLD   - potential temperature climatology,               2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  4.15533894E+00 Max =  3.54882013E+01)
    GET_3DFLD   - salinity climatology,                            2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -1.15116346E+01 Max =  3.68338484E+01)

 Basin information for Grid 01:

 Maximum grid stiffness ratios:  rx0 =   9.932805E-01 (Beckmann and Haidvogel)
                                 rx1 =   2.837318E+01 (Haney)

 Initial domain volumes:  TotVolume =  3.7844634699E+14 m3
                         MinCellVol =  2.0730148611E+04 m3
                         MaxCellVol =  3.4788293472E+10 m3
                            Max/Min =  1.6781497385E+06


 NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000064877761 - 000064888560)

    GET_2DFLD   - surface u-wind component,                        2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min = -1.20530137E+01 Max =  5.51073609E-01)      regrid = T
    GET_2DFLD   - surface v-wind component,                        2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min = -7.55569297E+00 Max =  3.14627136E+00)      regrid = T
    GET_2DFLD   - surface air pressure,                            2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  1.00585810E+03 Max =  1.01232248E+03)      regrid = T
    GET_2DFLD   - net longwave radiation flux,                     2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min = -1.34115877E-05 Max = -1.42211893E-06)      regrid = T
    GET_2DFLD   - surface air temperature,                         2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  2.38640228E+01 Max =  3.15746287E+01)      regrid = T
    GET_2DFLD   - surface air relative humidity,                   2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  6.08187971E-01 Max =  8.83000289E-01)      regrid = T
    GET_2DFLD   - rain fall rate,                                  2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  2.56529711E-06 Max =  3.64055637E-04)      regrid = T
    GET_NGFLD   - 3D u-momentum eastern boundary condition,        2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -1.06390405E+00 Max =  2.70750403E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,        2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -4.67574030E-01 Max =  3.27160418E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -1.06390405E+00 Max =  5.65890431E-01)
    GET_NGFLD   - 3D v-momentum southern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -9.68096316E-01 Max =  6.13836169E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D v-momentum northern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - temperature eastern boundary condition,          2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  3.17175808E+01)
    GET_NGFLD   - salinity eastern boundary condition,             2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  3.63111801E+01)
    GET_NGFLD   - temperature southern boundary condition,         2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  4.21079731E+00 Max =  3.10041027E+01)
    GET_NGFLD   - salinity southern boundary condition,            2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  3.48317986E+01 Max =  3.66650505E+01)
    GET_3DFLD   - potential temperature climatology,               2005-08-28 00:00:00.00
                   (Grid=01, Rec=2, Index=2, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  4.19558276E+00 Max =  3.47792665E+01)
    GET_3DFLD   - salinity climatology,                            2005-08-28 00:00:00.00
                   (Grid=01, Rec=2, Index=2, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -1.77656485E+01 Max =  3.68394868E+01)

 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

  64877760 2005-08-27 00:00:00.00  1.244784E-02  1.152586E+04  1.152587E+04  3.786342E+14
                       (54,22,12)  2.753814E-03  7.700963E-04  0.000000E+00  1.049467E+00
      DEF_HIS     - creating  history      file, Grid 01: Katrina_ocean_his.nc
      WRT_HIS     - wrote history     fields (Index=1,1) in record = 1
      DEF_AVG     - creating  average      file, Grid 01: Katrina_ocean_avg.nc
  64877761 2005-08-27 00:00:20.00  1.237147E-02  1.152583E+04  1.152584E+04  3.786334E+14
                       (58,55,11)  1.284955E-04  4.551020E-05  5.362781E-02  4.036428E+00
  64877762 2005-08-27 00:00:40.00  1.230128E-01  1.152581E+04  1.152593E+04  3.786328E+14
                       (95,64,01)  9.269159E+00  0.000000E+00  1.351492E-02  3.890935E+03
 DIAG speed trouble    3890.9349835834601     
 DIAG speed ijk    5431.7885291991397               96          64           1           1
 Found Error: 01   Line: 341      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 331      Source: ROMS/Drivers/nl_ocean.h, ROMS_run

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

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

Elapsed wall CPU time for each process (seconds):

 Node   #    0 CPU:       0.451
 Total:                   0.451

 Nonlinear model elapsed CPU time profile, Grid: 01

  Allocation and array initialization ..............         0.030  ( 6.6131 %)
  Ocean state initialization .......................         0.004  ( 0.8631 %)
  Reading of input data ............................         0.071  (15.8191 %)
  Processing of input data .........................         0.002  ( 0.5284 %)
  Processing of output time averaged data ..........         0.001  ( 0.1359 %)
  Computation of vertical boundary conditions ......         0.000  ( 0.0743 %)
  Computation of global information integrals ......         0.004  ( 0.8595 %)
  Writing of output data ...........................         0.026  ( 5.7028 %)
  Model 2D kernel ..................................         0.047  (10.4547 %)
  Tidal forcing ....................................         0.001  ( 0.2089 %)
  2D/3D coupling, vertical metrics .................         0.005  ( 1.2165 %)
  Omega vertical velocity ..........................         0.003  ( 0.6947 %)
  Equation of state for seawater ...................         0.013  ( 2.9513 %)
  Atmosphere-Ocean bulk flux parameterization ......         0.021  ( 4.7455 %)
  GLS vertical mixing parameterization .............         0.114  (25.2444 %)
  3D equations right-side terms ....................         0.008  ( 1.8558 %)
  3D equations predictor step ......................         0.021  ( 4.6599 %)
  Pressure gradient ................................         0.008  ( 1.6738 %)
  Harmonic mixing of tracers, geopotentials ........         0.008  ( 1.7935 %)
  Harmonic stress tensor, S-surfaces ...............         0.004  ( 0.8223 %)
  Corrector time-step for 3D momentum ..............         0.009  ( 1.9268 %)
  Corrector time-step for tracers ..................         0.013  ( 2.8220 %)
                                              Total:         0.413   91.6661 %

  Unique kernel(s) regions profiled ................         0.413   91.6661 %
  Residual, non-profiled code ......................         0.038    8.3339 %


 All percentages are with respect to total time =            0.451


 MPI communications profile, Grid: 01

  Message Passage: 2D halo exchanges ...............         0.000  ( 0.0889 %)
  Message Passage: 3D halo exchanges ...............         0.000  ( 0.0186 %)
  Message Passage: 4D halo exchanges ...............         0.000  ( 0.0055 %)
  Message Passage: data broadcast ..................         0.002  ( 0.4952 %)
  Message Passage: data reduction ..................         0.000  ( 0.0162 %)
  Message Passage: data gathering ..................         0.009  ( 1.9188 %)
  Message Passage: data scattering..................         0.007  ( 1.4795 %)
  Message Passage: boundary data gathering .........         0.000  ( 0.0047 %)
  Message Passage: point data gathering ............         0.000  ( 0.0002 %)
  Message Passage: synchronization barrier .........         0.000  ( 0.0007 %)
                                              Total:         0.018    4.0283 %

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

 Dynamic and Automatic memory (MB) usage for Grid 01:  98x98x15  tiling: 1x1

     tile          Dynamic        Automatic            USAGE      MPI-Buffers

        0            79.38             9.74            89.13             1.36

    TOTAL            79.38             9.74            89.13             1.36

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

 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_fsobc.h
     ROMS/Functionals/ana_m2obc.h
     ROMS/Functionals/ana_nudgcoef.h
     ROMS/Functionals/ana_srflux.h

 MAIN: Abnormal termination: BLOWUP.
 REASON: MaxSpeed =  3.9E+03
What should I do?

Input/Output

Code: Select all

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

      10800  ntimes            Number of timesteps for 3-D equations.
     20.000  dt                Timestep size (s) for 3-D equations.
         18  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.
         60  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).
         60  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.
        432  nAVG              Number of timesteps between the writing of
                                 time-averaged data into averages file.
 2.0000E-01  nl_tnu2(01)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 01: temp
 2.0000E-01  nl_tnu2(02)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 02: salt
 1.0000E-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
 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-03  Akt_limit(01)     Vertical diffusion upper threshold (m2/s)
                                 for tracer 01: temp
 1.0000E-03  Akt_limit(02)     Vertical diffusion upper threshold (m2/s)
                                 for tracer 02: salt
 1.0000E-05  Akv_bak           Background vertical mixing coefficient (m2/s)
                                 for momentum.
 1.0000E-03  Akv_limit         Vertical viscosity upper threshold (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.
      3.000  gls_p             GLS stability exponent.
      1.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.
 1.4400E+00  gls_c1            GLS shear production coefficient.
 1.9200E+00  gls_c2            GLS dissipation coefficient.
-4.0000E-01  gls_c3m           GLS stable buoyancy production coefficient.
 1.0000E+00  gls_c3p           GLS unstable buoyancy production coefficient.
 1.0000E+00  gls_sigk          GLS constant Schmidt number for TKE.
 1.3000E+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.
      0.000  wec_alpha         WEC factor for roller/breaking energy distribution.
 3.0000E-04  rdrg              Linear bottom drag coefficient (m/s).
 2.5000E-02  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.
 1.0000E+01  blk_ZW            Height (m) of surface winds measurement.
          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.
  38591.000  dstart            Time-stamp assigned to model initialization (days).
  38591.000  tide_start        Reference time origin for tidal forcing (days).
18581117.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
 0.0000E+00  Znudg             Nudging/relaxation time scale (days)
                                 for free-surface.
 0.0000E+00  M2nudg            Nudging/relaxation time scale (days)
                                 for 2D momentum.
 1.0000E+00  M3nudg            Nudging/relaxation time scale (days)
                                 for 3D momentum.
 0.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.
     10.000  T0                Background potential temperature (C) constant.
     30.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  LsshCLM           Turning OFF processing of SSH climatology.
          F  Lm2CLM            Turning OFF processing of 2D momentum climatology.
          F  Lm3CLM            Turning OFF 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.
          F  LnudgeM3CLM       Turning OFF 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(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(idUsms)      Write out surface U-momentum stress.
          T  Hout(idVsms)      Write out surface V-momentum stress.
          T  Hout(idUbms)      Write out bottom U-momentum stress.
          T  Hout(idVbms)      Write out bottom V-momentum stress.
          T  Hout(idPair)      Write out surface air pressure.
          T  Hout(idTair)      Write out surface air temperature.
          T  Hout(idUair)      Write out surface U-wind component.
          T  Hout(idVair)      Write out surface V-wind component.
          T  Hout(idUairE)     Write out surface Eastward U-wind component.
          T  Hout(idVairN)     Write out surface Northward V-wind component.
          T  Hout(idTsur)      Write out surface net heat flux.
          T  Hout(idTsur)      Write out surface net salt flux.
          T  Hout(idSrad)      Write out shortwave radiation flux.
          T  Hout(idLrad)      Write out longwave radiation flux.
          T  Hout(idLhea)      Write out latent heat flux.
          T  Hout(idShea)      Write out sensible heat flux.
          T  Hout(idEmPf)      Write out E-P flux.
          T  Hout(idevap)      Write out evaporation rate.
          T  Hout(idrain)      Write out rain rate.
          T  Hout(idVvis)      Write out vertical viscosity: AKv.
          T  Hout(idMtke)      Write out turbulent kinetic energy.
          T  Hout(idMtls)      Write out turbulent generic length-scale.

          T  Aout(idSdif)      Write out averaged vertical diffusion: AKt(isalt).

 Output/Input Files:

               Output Restart File:  Katrina_ocean_rst.nc
               Output History File:  Katrina_ocean_his.nc
              Output Averages File:  Katrina_ocean_avg.nc
                   Input Grid File:  Projects/Katrina/roms_grid_Katrina.nc
      Input Nonlinear Initial File:  Projects/Katrina/Katrina_ini.nc
                Tidal Forcing File:  Projects/Katrina/tide_forc_Katrina.nc
             Input Forcing File 01:  Projects/Katrina/romsforc_NARR_Katrina2005.nc
         Input Climatology File 01:  Projects/Katrina/Katrina_clim.nc
    Input Lateral Boundary File 01:  Projects/Katrina/Katrina_bry.nc
  ROMS I/O variables Metadata File:  ROMS/External/varinfo.dat

 Tile partition information for Grid 01:  98x98x15  tiling: 1x1

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1       98        1       98   144060

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

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50    98.50     0.50    98.50  RHO-points

        0     1.00    98.00     0.50    98.50    U-points

        0     0.50    98.50     1.00    98.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     218
               HaloSizeJ(1) =     218
                TileSide(1) =     103
                TileSize(1) =   10609


 Tracer Advection Scheme: NLM

s193264
Posts: 17
Joined: Wed Apr 06, 2022 5:59 pm
Location: College

Re: An error occurred when running ROMS in COAWST

#4 Unread post by s193264 »

I have emailed you the climatology, boundary conditions, and initial condition files that I used.
I created these using the seapy utility.
Are there any strange parts?

jcwarner
Posts: 1098
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: An error occurred when running ROMS in COAWST

#5 Unread post by jcwarner »

I see still:

GET_GRID - bathymetry at RHO-points: h
(Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
(Min = 5.95202990E-03 Max = 3.35037820E+03)

that is a min of 0.005 m . I dont think you need a min of 5 mm.
try again. in the grid file to set h to have a minimum depth. maybe something like

h=ncread('roms_grid_Katrina.nc','h');
h(h<3)=3;
ncwrite('roms_grid_Katrina.nc','h',h);

if it blows up again, especially if it tells you the i,j indices. then you need to look at the rst file and see what is going on at that location.
-j

s193264
Posts: 17
Joined: Wed Apr 06, 2022 5:59 pm
Location: College

Re: An error occurred when running ROMS in COAWST

#6 Unread post by s193264 »

jcwarner wrote: Sat Jul 30, 2022 1:52 pm I see still:

GET_GRID - bathymetry at RHO-points: h
(Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
(Min = 5.95202990E-03 Max = 3.35037820E+03)

that is a min of 0.005 m . I dont think you need a min of 5 mm.
try again. in the grid file to set h to have a minimum depth. maybe something like

h=ncread('roms_grid_Katrina.nc','h');
h(h<3)=3;
ncwrite('roms_grid_Katrina.nc','h',h);

if it blows up again, especially if it tells you the i,j indices. then you need to look at the rst file and see what is going on at that location.
-j

Dear jcwarner

Thank you for your advice.
I used MATLAB and changed the minimum value of h to 3 and ran it again.
Based on my output, the problem seems to have been solved by your advice.

But then I get another error.
I checked the rst file, but I'm not sure what's going on.
I don't know how to send the rst file, so I will send it to your email (jcwarner@usgs.gov).
If you have time, I would appreciate it if you could check it.
I'm sorry, but I'd appreciate it if you could let me know.

<roms.out>

Code: Select all

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

      10800  ntimes            Number of timesteps for 3-D equations.
     30.000  dt                Timestep size (s) for 3-D equations.
         28  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.
         60  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).
         60  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.
        432  nAVG              Number of timesteps between the writing of
                                 time-averaged data into averages file.
 2.0000E-01  nl_tnu2(01)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 01: temp
 2.0000E-01  nl_tnu2(02)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 02: salt
 1.0000E-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
 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-03  Akt_limit(01)     Vertical diffusion upper threshold (m2/s)
                                 for tracer 01: temp
 1.0000E-03  Akt_limit(02)     Vertical diffusion upper threshold (m2/s)
                                 for tracer 02: salt
 1.0000E-05  Akv_bak           Background vertical mixing coefficient (m2/s)
                                 for momentum.
 1.0000E-03  Akv_limit         Vertical viscosity upper threshold (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.
      3.000  gls_p             GLS stability exponent.
      1.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.
 1.4400E+00  gls_c1            GLS shear production coefficient.
 1.9200E+00  gls_c2            GLS dissipation coefficient.
-4.0000E-01  gls_c3m           GLS stable buoyancy production coefficient.
 1.0000E+00  gls_c3p           GLS unstable buoyancy production coefficient.
 1.0000E+00  gls_sigk          GLS constant Schmidt number for TKE.
 1.3000E+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.
      0.000  wec_alpha         WEC factor for roller/breaking energy distribution.
 3.0000E-04  rdrg              Linear bottom drag coefficient (m/s).
 2.5000E-02  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.
 1.0000E+01  blk_ZW            Height (m) of surface winds measurement.
          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.
  38591.000  dstart            Time-stamp assigned to model initialization (days).
  38591.000  tide_start        Reference time origin for tidal forcing (days).
18581117.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
 0.0000E+00  Znudg             Nudging/relaxation time scale (days)
                                 for free-surface.
 0.0000E+00  M2nudg            Nudging/relaxation time scale (days)
                                 for 2D momentum.
 1.0000E+00  M3nudg            Nudging/relaxation time scale (days)
                                 for 3D momentum.
 0.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.
     10.000  T0                Background potential temperature (C) constant.
     30.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  LsshCLM           Turning OFF processing of SSH climatology.
          F  Lm2CLM            Turning OFF processing of 2D momentum climatology.
          F  Lm3CLM            Turning OFF 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.
          F  LnudgeM3CLM       Turning OFF 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(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(idUsms)      Write out surface U-momentum stress.
          T  Hout(idVsms)      Write out surface V-momentum stress.
          T  Hout(idUbms)      Write out bottom U-momentum stress.
          T  Hout(idVbms)      Write out bottom V-momentum stress.
          T  Hout(idPair)      Write out surface air pressure.
          T  Hout(idTair)      Write out surface air temperature.
          T  Hout(idUair)      Write out surface U-wind component.
          T  Hout(idVair)      Write out surface V-wind component.
          T  Hout(idUairE)     Write out surface Eastward U-wind component.
          T  Hout(idVairN)     Write out surface Northward V-wind component.
          T  Hout(idTsur)      Write out surface net heat flux.
          T  Hout(idTsur)      Write out surface net salt flux.
          T  Hout(idSrad)      Write out shortwave radiation flux.
          T  Hout(idLrad)      Write out longwave radiation flux.
          T  Hout(idLhea)      Write out latent heat flux.
          T  Hout(idShea)      Write out sensible heat flux.
          T  Hout(idEmPf)      Write out E-P flux.
          T  Hout(idevap)      Write out evaporation rate.
          T  Hout(idrain)      Write out rain rate.
          T  Hout(idVvis)      Write out vertical viscosity: AKv.
          T  Hout(idMtke)      Write out turbulent kinetic energy.
          T  Hout(idMtls)      Write out turbulent generic length-scale.

          T  Aout(idSdif)      Write out averaged vertical diffusion: AKt(isalt).

 Output/Input Files:

               Output Restart File:  Katrina_ocean_rst.nc
               Output History File:  Katrina_ocean_his.nc
              Output Averages File:  Katrina_ocean_avg.nc
                   Input Grid File:  Projects/Katrina/roms_grid_Katrina.nc
      Input Nonlinear Initial File:  Projects/Katrina/Katrina_ini.nc
                Tidal Forcing File:  Projects/Katrina/tide_forc_Katrina.nc
             Input Forcing File 01:  Projects/Katrina/romsforc_NARR_Katrina2005.nc
         Input Climatology File 01:  Projects/Katrina/Katrina_clim.nc
    Input Lateral Boundary File 01:  Projects/Katrina/Katrina_bry.nc
  ROMS I/O variables Metadata File:  ROMS/External/varinfo.dat

 Tile partition information for Grid 01:  98x98x15  tiling: 1x1

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1       98        1       98   144060

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

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50    98.50     0.50    98.50  RHO-points

        0     1.00    98.00     0.50    98.50    U-points

        0     0.50    98.50     1.00    98.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     218
               HaloSizeJ(1) =     218
                TileSide(1) =     103
                TileSize(1) =   10609


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

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

 temp                     1   Upstream3          Centered4

 salt                     1   Upstream3          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  Chapman Imp

 ubar                     1   Closed       Flather      Flather      Flather

 vbar                     1   Closed       Flather      Flather      Flather

 u                        1   Closed       Nested       Nested       Nested

 v                        1   Closed       Nested       Nested       Nested

 temp                     1   Gradient     Nested       Nested       Gradient

 salt                     1   Gradient     Nested       Nested       Gradient

 tke                      1   Closed       Gradient     Gradient     Gradient

 Activated C-preprocessing Options:

 KATRINA                  Hurricane Katrina
 ANA_BSFLUX               Analytical kinematic bottom salinity flux
 ANA_BTFLUX               Analytical kinematic bottom temperature flux
 ANA_FSOBC                Analytical free-surface boundary conditions
 ANA_M2OBC                Analytical 2D momentum boundary conditions
 ANA_NUDGCOEF             Analytical spatially varying nudging time-scales
 ANA_SRFLUX               Analytical kinematic shortwave radiation flux
 ASSUMED_SHAPE            Using assumed-shape arrays
 ATM_PRESS                Impose atmospheric pressure onto sea surface
 AVERAGES                 Writing out time-averaged nonlinear model fields
 !BOUNDARY_ALLGATHER      Using mpi_allreduce in mp_boundary routine
 BULK_FLUXES              Surface bulk fluxes parameterization
 !COLLECT_ALL...          Using mpi_isend/mpi_recv in mp_collect routine
 CHARNOK                  Charnok surface roughness from wind stress
 CRAIG_BANNER             Craig and Banner wave breaking surface flux
 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
 KANTHA_CLAYSON           Kantha and Clayson stability function formulation
 LIMIT_VDIFF              Impose an upper limit on vertical diffusion coefficient
 LIMIT_VVISC              Impose an upper limit on vertical viscosity coefficient
 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
 N2S2_HORAVG              Horizontal smoothing of buoyancy and shear
 POWER_LAW                Power-law shape time-averaging barotropic filter
 PROFILE                  Time profiling activated
 K_GSCHEME                Third-order upstream advection of TKE fields
 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              Double 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
 SSH_TIDES                Add tidal elevation to SSH climatology
 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_LOGDRAG               Logarithmic bottom stress
 UV_TIDES                 Add tidal currents to 2D momentum climatologies
 UV_VIS2                  Harmonic mixing of momentum
 VAR_RHO_2D               Variable density barotropic mode

 Process Information:

 Node #    0 (pid=  137729) is active.

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

    GET_GRID    - angle between XI-axis and EAST: angler
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -2.04122474E-12 Max =  2.76144271E-12)
    GET_GRID    - bathymetry at RHO-points: h
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  3.00000000E+00 Max =  3.35037820E+03)
    GET_GRID    - Coriolis parameter at RHO-points: f
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  6.41047132E-05 Max =  7.91128070E-05)
    GET_GRID    - reciprocal XI-grid spacing: pm
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  1.40936443E-04 Max =  1.50689667E-04)
    GET_GRID    - reciprocal ETA-grid spacing: pn
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  1.31369178E-04 Max =  1.31413711E-04)
    GET_GRID    - XI-derivative of inverse metric factor pn: dndx
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -5.92198376E-07 Max =  5.92198376E-07)
    GET_GRID    - ETA-derivative of inverse metric factor pm: dmde
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -5.10897846E+00 Max =  0.00000000E+00)
    GET_GRID    - longitude of RHO-points: lon_rho
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30376282E+01 Max = -8.60048832E+01)
    GET_GRID    - longitude of U-points: lon_u
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30021093E+01 Max = -8.60404021E+01)
    GET_GRID    - longitude of V-points: lon_v
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30376282E+01 Max = -8.60048832E+01)
    GET_GRID    - longitude of PSI-points: lon_psi
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min = -9.30021093E+01 Max = -8.60404021E+01)
    GET_GRID    - latitude of RHO-points lat_rho
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.60754971E+01 Max =  3.28514653E+01)
    GET_GRID    - latitude of U-points: lat_u
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.60754971E+01 Max =  3.28514653E+01)
    GET_GRID    - latitude of V-points: lat_v
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.61097280E+01 Max =  3.28172547E+01)
    GET_GRID    - latitude of PSI-points lat_psi
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  2.61097280E+01 Max =  3.28172547E+01)
    GET_GRID    - mask on RHO-points: mask_rho
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on U-points: mask_u
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on V-points: mask_v
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on PSI-points: mask_psi
                   (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.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.174      -0.671      -1.978      -2.652
    13  -0.1333333  -0.0017081       -0.348      -1.350      -5.465      -8.340
    12  -0.2000000  -0.0043062       -0.523      -2.043     -11.088     -18.318
    11  -0.2666667  -0.0089353       -0.699      -2.756     -20.076     -35.061
    10  -0.3333333  -0.0169029       -0.876      -3.502     -34.595     -62.922
     9  -0.4000000  -0.0304126       -1.055      -4.304     -58.297    -109.242
     8  -0.4666667  -0.0530642       -1.238      -5.197     -97.147    -186.008
     7  -0.5333333  -0.0905094       -1.427      -6.238    -160.508    -312.044
     6  -0.6000000  -0.1510347       -1.624      -7.510    -262.112    -514.949
     5  -0.6666667  -0.2452550       -1.835      -9.119    -419.545    -830.075
     4  -0.7333333  -0.3828872       -2.063     -11.162    -648.909   -1289.784
     3  -0.8000000  -0.5634305       -2.307     -13.634    -949.373   -1892.409
     2  -0.8666667  -0.7610631       -2.559     -16.277   -1278.153   -2551.949
     1  -0.9333333  -0.9210683       -2.795     -18.544   -1544.588   -3086.171
     0  -1.0000000  -1.0000000       -3.000     -20.000   -1676.689   -3350.378

 Time Splitting Weights for Grid 01:    ndtfast =  28    nfast =  40
 ==================================

    Primary            Secondary            Accumulated to Current Step

  1-0.0008549076326633 0.0357142857142857-0.0008549076326633 0.0357142857142857
  2-0.0014558918889834 0.0357448181297380-0.0023107995216466 0.0714591038440237
  3-0.0018029740646893 0.0357968142686302-0.0041137735863359 0.1072559181126539
  4-0.0018962369765050 0.0358612061995120-0.0060100105628409 0.1431171243121659
  5-0.0017358959479127 0.0359289289486729-0.0077459065107536 0.1790460532608388
  6-0.0013223981892219 0.0359909252325269-0.0090683046999755 0.2150369784933657
  7-0.0006565505719430 0.0360381537392848-0.0097248552719185 0.2510751322326505
  8 0.0002603242025330 0.0360616019739971-0.0094645310693855 0.2871367342066476
  9 0.0014262030399494 0.0360523046810495-0.0080383280294362 0.3231890388876971
 10 0.0028381494958941 0.0360013688581941-0.0052001785335421 0.3591904077458913
 11 0.0044920724242039 0.0359000063761979-0.0007081061093382 0.3950904141220892
 12 0.0063824562310637 0.0357395752181906 0.0056743501217256 0.4308299893402799
 13 0.0085020627348001 0.0355116303527955 0.0141764128565257 0.4663416196930754
 14 0.0108416046313699 0.0352079852551241 0.0250180174878955 0.5015496049481996
 15 0.0133893905655432 0.0348207850897180 0.0384074080534388 0.5363703900379175
 16 0.0161309418077813 0.0343425925695200 0.0545383498612201 0.5707129826074375
 17 0.0190485805368085 0.0337664875049564 0.0735869303980286 0.6044794701123939
 18 0.0221209897278798 0.0330861810572133 0.0957079201259084 0.6375656511696072
 19 0.0253227446467419 0.0322961457097890 0.1210306647726504 0.6698617968793962
 20 0.0286238159492901 0.0313917619724054 0.1496544807219404 0.7012535588518015
 21 0.0319890443869186 0.0303694828313593 0.1816435251088590 0.7316230416831608
 22 0.0353775871175665 0.0292270169603979 0.2170211122264255 0.7608500586435587
 23 0.0387423356224582 0.0279635317061991 0.2557634478488837 0.7888135903497578
 24 0.0420293052285376 0.0265798768625399 0.2977927530774214 0.8153934672122978
 25 0.0451769962365983 0.0250788302472350 0.3429697493140197 0.8404722974595328
 26 0.0481157266551073 0.0234653660959279 0.3910854759691270 0.8639376635554606
 27 0.0507669365397242 0.0217469472868169 0.4418524125088512 0.8856846108422776
 28 0.0530424639385140 0.0199338424103982 0.4948948764473652 0.9056184532526758
 29 0.0548437924428562 0.0180394686983084 0.5497386688902214 0.9236579219509842
 30 0.0560612703440466 0.0160807618253492 0.6057999392342680 0.9397386837763334
 31 0.0565733013955951 0.0140785735987762 0.6623732406298630 0.9538172573751096
 32 0.0562455071812177 0.0120580985489335 0.7186187478110807 0.9658753559240431
 33 0.0549298610885234 0.0100493304353186 0.7735486088996041 0.9759246863593616
 34 0.0524637938883948 0.0080875496821570 0.8260124027879989 0.9840122360415186
 35 0.0486692709200651 0.0062138427575715 0.8746816737080639 0.9902260787990901
 36 0.0433518408818877 0.0044756545104263 0.9180335145899516 0.9947017333095164
 37 0.0362996562278023 0.0029273744789303 0.9543331708177540 0.9976291077884467
 38 0.0272824651694943 0.0016309581850802 0.9816156359872483 0.9992600659735269
 39 0.0160505752842493 0.0006565844290269 0.9976662112714976 0.9999166504025537
 40 0.0023337887285027 0.0000833495974465 1.0000000000000002 1.0000000000000002

 ndtfast, nfast =   28  40   nfast/ndtfast =  1.42857

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

    1.000000000000 1.049755056052 0.524877528026 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.18257

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

 Minimum X-grid spacing, DXmin =  6.63615510E+00 km    Water points =  6.79507250E+00 km
 Maximum X-grid spacing, DXmax =  7.09539688E+00 km    Water points =  7.09539688E+00 km
 Minimum Y-grid spacing, DYmin =  7.60955604E+00 km    Water points =  7.61102264E+00 km
 Maximum Y-grid spacing, DYmax =  7.61213562E+00 km    Water points =  7.61213562E+00 km
 Minimum Z-grid spacing, DZmin =  1.74068899E-01 m     Water points =  1.74068899E-01 m
 Maximum Z-grid spacing, DZmax =  6.59540089E+02 m     Water points =  6.59540089E+02 m

 Minimum barotropic Courant Number =  1.13632420E-03
 Maximum barotropic Courant Number =  3.74243392E-02
 Maximum Coriolis   Courant Number =  2.23118437E-03


 NLM: GET_STATE - Reading state initial conditions,                       2005-08-27 00:00:00.00
                   (Grid 01, t = 53609.0000, File: Katrina_ini.nc, Rec=0001, Index=1)
                - free-surface
                   (Min = -3.54594313E-01 Max =  4.45445436E-01)
                - vertically integrated u-momentum component
                   (Min = -1.41080142E-01 Max =  2.52036088E-01)
                - vertically integrated v-momentum component
                   (Min = -1.71133164E-01 Max =  1.16655195E-01)
                - u-momentum component
                   (Min = -9.26081898E-01 Max =  1.03925475E+00)
                - v-momentum component
                   (Min = -9.95697845E-01 Max =  5.99459860E-01)
                - potential temperature
                   (Min =  4.15573828E+00 Max =  3.54843508E+01)
                - salinity
                   (Min = -1.15371107E+01 Max =  3.68318607E+01)
    GET_NGFLD   - tidal period
                   (Grid = 01, Min =  1.49047215E+04 Max =  9.67260840E+04)
    GET_2DFLD   - tidal elevation amplitude
                   (Grid = 01, Min =  0.00000000E+00 Max =  2.09046983E-01)
    GET_2DFLD   - tidal elevation phase angle
                   (Grid = 01, Min = -1.29687416E+01 Max =  4.82525186E+00)
    GET_2DFLD   - tidal current inclination angle
                   (Grid = 01, Min =  0.00000000E+00 Max =  6.28316919E+00)
    GET_2DFLD   - tidal current phase angle
                   (Grid = 01, Min =  0.00000000E+00 Max =  6.28306597E+00)
    GET_2DFLD   - maximum tidal current, ellipse major axis
                   (Grid = 01, Min =  0.00000000E+00 Max =  3.46792927E-01)
    GET_2DFLD   - minimum tidal current, ellipse minor axis
                   (Grid = 01, Min = -8.07614703E-02 Max =  4.40275504E-02)
    GET_2DFLD   - surface u-wind component,                        2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -1.07696928E+01 Max = -8.92840971E-01)      regrid = T
    GET_2DFLD   - surface v-wind component,                        2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -8.16269459E+00 Max =  3.43376645E+00)      regrid = T
    GET_2DFLD   - surface air pressure,                            2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  1.00612362E+03 Max =  1.01279447E+03)      regrid = T
    GET_2DFLD   - net longwave radiation flux,                     2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -1.29572887E-05 Max = -6.30481076E-07)      regrid = T
    GET_2DFLD   - surface air temperature,                         2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  2.74624980E+01 Max =  3.52368103E+01)      regrid = T
    GET_2DFLD   - surface air relative humidity,                   2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  4.85723022E-01 Max =  8.13929408E-01)      regrid = T
    GET_2DFLD   - rain fall rate,                                  2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  3.83558273E-07 Max =  3.75390348E-05)      regrid = T
    GET_NGFLD   - 3D u-momentum eastern boundary condition,        2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -3.29044104E-01 Max =  2.31331423E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,        2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -8.52144778E-01 Max =  1.83875933E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -9.26081896E-01 Max =  6.07348561E-01)
    GET_NGFLD   - 3D v-momentum southern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -9.95697856E-01 Max =  5.74903667E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D v-momentum northern boundary condition,       2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - temperature eastern boundary condition,          2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  3.17898636E+01)
    GET_NGFLD   - salinity eastern boundary condition,             2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  0.00000000E+00 Max =  3.63082619E+01)
    GET_NGFLD   - temperature southern boundary condition,         2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  4.21922779E+00 Max =  3.13313484E+01)
    GET_NGFLD   - salinity southern boundary condition,            2005-08-27 00:00:00.00
                   (Grid= 01, Rec=1, Index=2, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  3.48222885E+01 Max =  3.68282700E+01)
    GET_3DFLD   - potential temperature climatology,               2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min =  4.15573828E+00 Max =  3.54843508E+01)
    GET_3DFLD   - salinity climatology,                            2005-08-27 00:00:00.00
                   (Grid=01, Rec=1, Index=1, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.0000
                   (Min = -1.15371107E+01 Max =  3.68318607E+01)

 Basin information for Grid 01:

 Maximum grid stiffness ratios:  rx0 =   7.239834E-01 (Beckmann and Haidvogel)
                                 rx1 =   1.937196E+01 (Haney)

 Initial domain volumes:  TotVolume =  3.7847305688E+14 m3
                         MinCellVol =  9.0023996777E+06 m3
                         MaxCellVol =  3.4788293472E+10 m3
                            Max/Min =  3.8643355902E+03


 NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000043251841 - 000043262640)

    GET_2DFLD   - surface u-wind component,                        2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min = -1.20530137E+01 Max =  5.51073609E-01)      regrid = T
    GET_2DFLD   - surface v-wind component,                        2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min = -7.55569297E+00 Max =  3.14627136E+00)      regrid = T
    GET_2DFLD   - surface air pressure,                            2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  1.00585810E+03 Max =  1.01232248E+03)      regrid = T
    GET_2DFLD   - net longwave radiation flux,                     2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min = -1.34115877E-05 Max = -1.42211893E-06)      regrid = T
    GET_2DFLD   - surface air temperature,                         2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  2.38640228E+01 Max =  3.15746287E+01)      regrid = T
    GET_2DFLD   - surface air relative humidity,                   2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  6.08187971E-01 Max =  8.83000289E-01)      regrid = T
    GET_2DFLD   - rain fall rate,                                  2005-08-27 03:00:00.00
                   (Grid=01, Rec=2, Index=2, File: romsforc_NARR_Katrina2005.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53609.1250
                   (Min =  2.56529711E-06 Max =  3.64055637E-04)      regrid = T
    GET_NGFLD   - 3D u-momentum eastern boundary condition,        2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -1.06418061E+00 Max =  2.73843646E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,        2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -4.67660874E-01 Max =  3.26894641E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -1.06418061E+00 Max =  5.55075884E-01)
    GET_NGFLD   - 3D v-momentum southern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -9.68237638E-01 Max =  6.13834620E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D v-momentum northern boundary condition,       2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - temperature eastern boundary condition,          2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  3.17182083E+01)
    GET_NGFLD   - salinity eastern boundary condition,             2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  0.00000000E+00 Max =  3.63111801E+01)
    GET_NGFLD   - temperature southern boundary condition,         2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  4.21179533E+00 Max =  3.10040932E+01)
    GET_NGFLD   - salinity southern boundary condition,            2005-08-28 00:00:00.00
                   (Grid= 01, Rec=2, Index=1, File: Katrina_bry.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  3.48143768E+01 Max =  3.66702194E+01)
    GET_3DFLD   - potential temperature climatology,               2005-08-28 00:00:00.00
                   (Grid=01, Rec=2, Index=2, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min =  4.20409152E+00 Max =  3.47801044E+01)
    GET_3DFLD   - salinity climatology,                            2005-08-28 00:00:00.00
                   (Grid=01, Rec=2, Index=2, File: Katrina_clim.nc)
                   (Tmin=      53609.0000 Tmax=      53614.0000)      t =      53610.0000
                   (Min = -1.78632816E+01 Max =  3.68393278E+01)

 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

  43251840 2005-08-27 00:00:00.00  1.242553E-02  1.151213E+04  1.151214E+04  3.790855E+14
                       (54,22,12)  4.132545E-03  1.152803E-03  0.000000E+00  1.049187E+00
      DEF_HIS     - creating  history      file, Grid 01: Katrina_ocean_his.nc
      WRT_HIS     - wrote history     fields (Index=1,1) in record = 1
      DEF_AVG     - creating  average      file, Grid 01: Katrina_ocean_avg.nc
  43251841 2005-08-27 00:00:30.00  1.232743E-02  1.151209E+04  1.151210E+04  3.790844E+14
                       (51,43,01)  1.401545E-04  6.384282E-05  9.199845E-03  1.049271E+00
  43251842 2005-08-27 00:01:00.00  1.226325E-02  1.151206E+04  1.151208E+04  3.790835E+14
                       (01,01,15)  8.649880E-04  1.029863E-03  5.399073E-01  1.050004E+00
  43251843 2005-08-27 00:01:30.00  1.219683E-02  1.151204E+04  1.151206E+04  3.790828E+14
                       (02,02,15)  1.427515E-03  9.590337E-04  3.432228E-01  1.051444E+00
  43251844 2005-08-27 00:02:00.00  1.213459E-02  1.151203E+04  1.151204E+04  3.790821E+14
                       (02,02,15)  1.351349E-03  8.934342E-04  3.984512E-01  1.053346E+00
  43251845 2005-08-27 00:02:30.00  1.208108E-02  1.151201E+04  1.151202E+04  3.790813E+14
                       (02,03,15)  1.015484E-03  7.323119E-04  2.872467E-01  2.605005E+00
  43251846 2005-08-27 00:03:00.00  1.225973E-02  1.151199E+04  1.151200E+04  3.790805E+14
                       (03,03,15)  1.196924E-03  8.446957E-04  3.284440E-01  9.310640E+01
 DIAG speed trouble    93.106400122773167     
 DIAG speed ijk    131.67231528779021               96          64           1           1
 Found Error: 01   Line: 341      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 331      Source: ROMS/Drivers/nl_ocean.h, ROMS_run

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

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

Elapsed wall CPU time for each process (seconds):

 Node   #    0 CPU:       1.096
 Total:                   1.096

 Nonlinear model elapsed CPU time profile, Grid: 01

  Allocation and array initialization ..............         0.029  ( 2.6845 %)
  Ocean state initialization .......................         0.003  ( 0.2821 %)
  Reading of input data ............................         0.070  ( 6.3672 %)
  Processing of input data .........................         0.005  ( 0.4667 %)
  Processing of output time averaged data ..........         0.001  ( 0.1348 %)
  Computation of vertical boundary conditions ......         0.001  ( 0.0842 %)
  Computation of global information integrals ......         0.009  ( 0.8505 %)
  Writing of output data ...........................         0.033  ( 3.0298 %)
  Model 2D kernel ..................................         0.224  (20.4799 %)
  Tidal forcing ....................................         0.003  ( 0.2485 %)
  2D/3D coupling, vertical metrics .................         0.010  ( 0.9332 %)
  Omega vertical velocity ..........................         0.009  ( 0.7986 %)
  Equation of state for seawater ...................         0.027  ( 2.4262 %)
  Atmosphere-Ocean bulk flux parameterization ......         0.065  ( 5.9534 %)
  GLS vertical mixing parameterization .............         0.349  (31.8597 %)
  3D equations right-side terms ....................         0.025  ( 2.2719 %)
  3D equations predictor step ......................         0.064  ( 5.8738 %)
  Pressure gradient ................................         0.022  ( 2.0415 %)
  Harmonic mixing of tracers, geopotentials ........         0.025  ( 2.2475 %)
  Harmonic stress tensor, S-surfaces ...............         0.011  ( 1.0311 %)
  Corrector time-step for 3D momentum ..............         0.026  ( 2.4149 %)
  Corrector time-step for tracers ..................         0.039  ( 3.5814 %)
                                              Total:         1.053   96.0613 %

  Unique kernel(s) regions profiled ................         1.053   96.0613 %
  Residual, non-profiled code ......................         0.043    3.9387 %


 All percentages are with respect to total time =            1.096


 MPI communications profile, Grid: 01

  Message Passage: 2D halo exchanges ...............         0.002  ( 0.1480 %)
  Message Passage: 3D halo exchanges ...............         0.000  ( 0.0195 %)
  Message Passage: 4D halo exchanges ...............         0.000  ( 0.0053 %)
  Message Passage: data broadcast ..................         0.002  ( 0.1903 %)
  Message Passage: data reduction ..................         0.000  ( 0.0110 %)
  Message Passage: data gathering ..................         0.008  ( 0.7669 %)
  Message Passage: data scattering..................         0.007  ( 0.5963 %)
  Message Passage: boundary data gathering .........         0.000  ( 0.0057 %)
  Message Passage: point data gathering ............         0.000  ( 0.0001 %)
  Message Passage: synchronization barrier .........         0.000  ( 0.0002 %)
                                              Total:         0.019    1.7432 %

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

 Dynamic and Automatic memory (MB) usage for Grid 01:  98x98x15  tiling: 1x1

     tile          Dynamic        Automatic            USAGE      MPI-Buffers

        0            79.38             9.74            89.13             1.36

    TOTAL            79.38             9.74            89.13             1.36

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

 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_fsobc.h
     ROMS/Functionals/ana_m2obc.h
     ROMS/Functionals/ana_nudgcoef.h
     ROMS/Functionals/ana_srflux.h

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


User avatar
arango
Site Admin
Posts: 1274
Joined: Wed Feb 26, 2003 4:41 pm
Location: DMCS, Rutgers University
Contact:

Re: An error occurred when running ROMS in COAWST

#7 Unread post by arango »

Gosh, you need to pay attention to what ROMS tells you in the standard output. Where in the ocean can we have velocities of 93 m/s, and what does that mean? Your application is blowing-up because of a violation of the CFL condition. If you don't know what a CFL condition is and how to remediate it, I will recommend reading the literature or searching in your browser. We don't have time to guide every step that you take. It is advantageous for numerical modeling users to take basic Geophysical Fluid Dynamics courses during their studies.

s193264
Posts: 17
Joined: Wed Apr 06, 2022 5:59 pm
Location: College

Re: An error occurred when running ROMS in COAWST

#8 Unread post by s193264 »

arango wrote: Sat Jul 30, 2022 2:44 pm Gosh, you need to pay attention to what ROMS tells you in the standard output. Where in the ocean can we have velocities of 93 m/s, and what does that mean? Your application is blowing-up because of a violation of the CFL condition. If you don't know what a CFL condition is and how to remediate it, I will recommend reading the literature or searching in your browser. We don't have time to guide every step that you take. It is advantageous for numerical modeling users to take basic Geophysical Fluid Dynamics courses during their studies.
Hi arango.
Thanks for your advice.
First of all, I am sorry that you are hearing my problems with basic numerical calculations due to my lack of learning.
However, since I originally did not major in marine systems, I think that my knowledge of Geophysical Fluid Dynamics is lacking.
In the future, I plan to study these fields in addition to operating the ROMS model.

Now, as you said, I checked my browser for violations of CFL conditions and how to fix them.
According to it, it was found that it is necessary to reduce the time step in order to solve the CFL condition violation.
However, in my calculation conditions, dt is 30, which seems to be sufficiently small.
Is my perception wrong? ?
Or is there a way to resolve CFL condition violations that I haven't looked into yet?

jcwarner
Posts: 1098
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: An error occurred when running ROMS in COAWST

#9 Unread post by jcwarner »

DIAG speed trouble 93.106400122773167
DIAG speed ijk 131.67231528779021 96 64 1 1
Found Error: 01 Line: 341 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 331 Source: ROMS/Drivers/nl_ocean.h, ROMS_run
Blowing-up: Saving latest model state into RESTART file
REASON: MaxSpeed = 9.3E+01
WRT_RST - wrote re-start fields (Index=1,1) in record = 1

please ncview the rst file, look near i = 96, j = 64 and see what it looks like.

s193264
Posts: 17
Joined: Wed Apr 06, 2022 5:59 pm
Location: College

Re: An error occurred when running ROMS in COAWST

#10 Unread post by s193264 »

Dear jcwarner, arango

Thanks to Arango and jcwarner, I was able to do it successfully.
The cause of the error was a violation of the CFL condition due to the small timestep.
Thank you guys.

User avatar
wilkin
Posts: 756
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: An error occurred when running ROMS in COAWST

#11 Unread post by wilkin »

I notice you have salinities in the initial conditions (and boundary conditions) of minus 11

Code: Select all

NLM: GET_STATE - Reading state initial conditions,                       2005-08-27 00:00:00.00
                   (Grid 01, t = 53609.0000, File: Katrina_ini.nc, Rec=0001, Index=1)
 ...
                - salinity
                   (Min = -1.15371107E+01 Max =  3.68318607E+01)
That is not going to go well with the nonlinear equation of state.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

Post Reply