Blow up problem in ROMS 3.9

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
ebpbenke
Posts: 2
Joined: Tue Dec 03, 2019 7:56 pm
Location: EBP Deutschland GmbH

Blow up problem in ROMS 3.9

#1 Unread post by ebpbenke »

Hello,

since the update from ROMS version 3.8 to 3.9 we get an BLOWUP error message of the following kind:

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

In the past we have used version 3.6, 3.7 and 3.8 for different regions. Our data were not changed.
Since the update to version 3.9 we get for all regions an BLOWUP error.

It follow parts of the output of a model run:

Code: Select all

--------------------------------------------------------------------------------
 Model Input Parameters:  ROMS/TOMS version 3.9  
                          Wednesday - May 13, 2020 -  7:21:44 AM
--------------------------------------------------------------------------------
 Operating system  : Linux
 CPU/hardware      : x86_64
 Compiler system   : ifort
 Compiler command  : /opt/nec/ve/mpi/2.6.0/bin64/mpinfort
 Compiler flags    : -fp-model precise -heap-arrays -ip -O3 -traceback -check uninit -sh 
 OCN Communicator  : 91,  PET size = 55
 Resolution, Grid 01: 414x176x32,  Parallel Nodes: 55,  Tiling: 11x5
 Tile minimum and maximum fractional coordinates for Grid 01:
   (interior points only)
     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50    36.50     0.50    34.50  RHO-points
        1    36.50    74.50     0.50    34.50  RHO-points
        2    74.50   112.50     0.50    34.50  RHO-points
        3   112.50   150.50     0.50    34.50  RHO-points
        4   150.50   188.50     0.50    34.50  RHO-points
        5   188.50   226.50     0.50    34.50  RHO-points
        6   226.50   264.50     0.50    34.50  RHO-points
        7   264.50   302.50     0.50    34.50  RHO-points
        8   302.50   340.50     0.50    34.50  RHO-points
        9   340.50   378.50     0.50    34.50  RHO-points
       10   378.50   414.50     0.50    34.50  RHO-points
       11     0.50    36.50    34.50    70.50  RHO-points
       12    36.50    74.50    34.50    70.50  RHO-points
       13    74.50   112.50    34.50    70.50  RHO-points
       14   112.50   150.50    34.50    70.50  RHO-points
       15   150.50   188.50    34.50    70.50  RHO-points
       16   188.50   226.50    34.50    70.50  RHO-points
       17   226.50   264.50    34.50    70.50  RHO-points
       18   264.50   302.50    34.50    70.50  RHO-points
       19   302.50   340.50    34.50    70.50  RHO-points
       20   340.50   378.50    34.50    70.50  RHO-points
       21   378.50   414.50    34.50    70.50  RHO-points
       22     0.50    36.50    70.50   106.50  RHO-points
       23    36.50    74.50    70.50   106.50  RHO-points
       24    74.50   112.50    70.50   106.50  RHO-points
       25   112.50   150.50    70.50   106.50  RHO-points
       26   150.50   188.50    70.50   106.50  RHO-points
       27   188.50   226.50    70.50   106.50  RHO-points
       28   226.50   264.50    70.50   106.50  RHO-points
       29   264.50   302.50    70.50   106.50  RHO-points
       30   302.50   340.50    70.50   106.50  RHO-points
       31   340.50   378.50    70.50   106.50  RHO-points
       32   378.50   414.50    70.50   106.50  RHO-points
       33     0.50    36.50   106.50   142.50  RHO-points
       34    36.50    74.50   106.50   142.50  RHO-points
       35    74.50   112.50   106.50   142.50  RHO-points
       36   112.50   150.50   106.50   142.50  RHO-points
       37   150.50   188.50   106.50   142.50  RHO-points
       38   188.50   226.50   106.50   142.50  RHO-points
       39   226.50   264.50   106.50   142.50  RHO-points
       40   264.50   302.50   106.50   142.50  RHO-points
       41   302.50   340.50   106.50   142.50  RHO-points
       42   340.50   378.50   106.50   142.50  RHO-points
       43   378.50   414.50   106.50   142.50  RHO-points
       44     0.50    36.50   142.50   176.50  RHO-points
       45    36.50    74.50   142.50   176.50  RHO-points
       46    74.50   112.50   142.50   176.50  RHO-points
       47   112.50   150.50   142.50   176.50  RHO-points
       48   150.50   188.50   142.50   176.50  RHO-points
       49   188.50   226.50   142.50   176.50  RHO-points
       50   226.50   264.50   142.50   176.50  RHO-points
       51   264.50   302.50   142.50   176.50  RHO-points
       52   302.50   340.50   142.50   176.50  RHO-points
       53   340.50   378.50   142.50   176.50  RHO-points
       54   378.50   414.50   142.50   176.50  RHO-points

        0     1.00    36.50     0.50    34.50    U-points
        1    36.50    74.50     0.50    34.50    U-points
        2    74.50   112.50     0.50    34.50    U-points
        3   112.50   150.50     0.50    34.50    U-points
        4   150.50   188.50     0.50    34.50    U-points
        5   188.50   226.50     0.50    34.50    U-points
        6   226.50   264.50     0.50    34.50    U-points
        7   264.50   302.50     0.50    34.50    U-points
        8   302.50   340.50     0.50    34.50    U-points
        9   340.50   378.50     0.50    34.50    U-points
       10   378.50   414.00     0.50    34.50    U-points
       11     1.00    36.50    34.50    70.50    U-points
       12    36.50    74.50    34.50    70.50    U-points
       13    74.50   112.50    34.50    70.50    U-points
       14   112.50   150.50    34.50    70.50    U-points
       15   150.50   188.50    34.50    70.50    U-points
       16   188.50   226.50    34.50    70.50    U-points
       17   226.50   264.50    34.50    70.50    U-points
       18   264.50   302.50    34.50    70.50    U-points
       19   302.50   340.50    34.50    70.50    U-points
       20   340.50   378.50    34.50    70.50    U-points
       21   378.50   414.00    34.50    70.50    U-points
       22     1.00    36.50    70.50   106.50    U-points
       23    36.50    74.50    70.50   106.50    U-points
       24    74.50   112.50    70.50   106.50    U-points
       25   112.50   150.50    70.50   106.50    U-points
       26   150.50   188.50    70.50   106.50    U-points
       27   188.50   226.50    70.50   106.50    U-points
       28   226.50   264.50    70.50   106.50    U-points
       29   264.50   302.50    70.50   106.50    U-points
       30   302.50   340.50    70.50   106.50    U-points
       31   340.50   378.50    70.50   106.50    U-points
       32   378.50   414.00    70.50   106.50    U-points
       33     1.00    36.50   106.50   142.50    U-points
       34    36.50    74.50   106.50   142.50    U-points
       35    74.50   112.50   106.50   142.50    U-points
       36   112.50   150.50   106.50   142.50    U-points
       37   150.50   188.50   106.50   142.50    U-points
       38   188.50   226.50   106.50   142.50    U-points
       39   226.50   264.50   106.50   142.50    U-points
       40   264.50   302.50   106.50   142.50    U-points
       41   302.50   340.50   106.50   142.50    U-points
       42   340.50   378.50   106.50   142.50    U-points
       43   378.50   414.00   106.50   142.50    U-points
       44     1.00    36.50   142.50   176.50    U-points
       45    36.50    74.50   142.50   176.50    U-points
       46    74.50   112.50   142.50   176.50    U-points
       47   112.50   150.50   142.50   176.50    U-points
       48   150.50   188.50   142.50   176.50    U-points
       49   188.50   226.50   142.50   176.50    U-points
       50   226.50   264.50   142.50   176.50    U-points
       51   264.50   302.50   142.50   176.50    U-points
       52   302.50   340.50   142.50   176.50    U-points
       53   340.50   378.50   142.50   176.50    U-points
       54   378.50   414.00   142.50   176.50    U-points

        0     0.50    36.50     1.00    34.50    V-points
        1    36.50    74.50     1.00    34.50    V-points
        2    74.50   112.50     1.00    34.50    V-points
        3   112.50   150.50     1.00    34.50    V-points
        4   150.50   188.50     1.00    34.50    V-points
        5   188.50   226.50     1.00    34.50    V-points
        6   226.50   264.50     1.00    34.50    V-points
        7   264.50   302.50     1.00    34.50    V-points
        8   302.50   340.50     1.00    34.50    V-points
        9   340.50   378.50     1.00    34.50    V-points
       10   378.50   414.50     1.00    34.50    V-points
       11     0.50    36.50    34.50    70.50    V-points
       12    36.50    74.50    34.50    70.50    V-points
       13    74.50   112.50    34.50    70.50    V-points
       14   112.50   150.50    34.50    70.50    V-points
       15   150.50   188.50    34.50    70.50    V-points
       16   188.50   226.50    34.50    70.50    V-points
       17   226.50   264.50    34.50    70.50    V-points
       18   264.50   302.50    34.50    70.50    V-points
       19   302.50   340.50    34.50    70.50    V-points
       20   340.50   378.50    34.50    70.50    V-points
       21   378.50   414.50    34.50    70.50    V-points
       22     0.50    36.50    70.50   106.50    V-points
       23    36.50    74.50    70.50   106.50    V-points
       24    74.50   112.50    70.50   106.50    V-points
       25   112.50   150.50    70.50   106.50    V-points
       26   150.50   188.50    70.50   106.50    V-points
       27   188.50   226.50    70.50   106.50    V-points
       28   226.50   264.50    70.50   106.50    V-points
       29   264.50   302.50    70.50   106.50    V-points
       30   302.50   340.50    70.50   106.50    V-points
       31   340.50   378.50    70.50   106.50    V-points
       32   378.50   414.50    70.50   106.50    V-points
       33     0.50    36.50   106.50   142.50    V-points
       34    36.50    74.50   106.50   142.50    V-points
       35    74.50   112.50   106.50   142.50    V-points
       36   112.50   150.50   106.50   142.50    V-points
       37   150.50   188.50   106.50   142.50    V-points
       38   188.50   226.50   106.50   142.50    V-points
       39   226.50   264.50   106.50   142.50    V-points
       40   264.50   302.50   106.50   142.50    V-points
       41   302.50   340.50   106.50   142.50    V-points
       42   340.50   378.50   106.50   142.50    V-points
       43   378.50   414.50   106.50   142.50    V-points
       44     0.50    36.50   142.50   176.00    V-points
       45    36.50    74.50   142.50   176.00    V-points
       46    74.50   112.50   142.50   176.00    V-points
       47   112.50   150.50   142.50   176.00    V-points
       48   150.50   188.50   142.50   176.00    V-points
       49   188.50   226.50   142.50   176.00    V-points
       50   226.50   264.50   142.50   176.00    V-points
       51   264.50   302.50   142.50   176.00    V-points
       52   302.50   340.50   142.50   176.00    V-points
       53   340.50   378.50   142.50   176.00    V-points
       54   378.50   414.50   142.50   176.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     100
               HaloSizeJ(1) =      96
                TileSide(1) =      44
                TileSize(1) =    1848


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

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

 temp                     1   varinfo.dat        ROMS/Utilit

 salt                     1   NDS/boundar        y/distribut
 Lateral Boundary Conditions: NLM
 ============================

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

 zeta                     1   Chapman Imp  Chapman Imp  Chapman Imp  Chapman Imp

 ubar                     1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 vbar                     1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

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

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

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

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

 tke                      1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

   GET_GRID    - bathymetry at RHO-points: h
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  1.00000000E+01 Max =  2.50000000E+03)
    GET_GRID    - Coriolis parameter at RHO-points: f
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  2.59725060E-05 Max =  3.58616573E-05)
    GET_GRID    - reciprocal XI-grid spacing: pm
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  3.96040493E-04 Max =  4.02057741E-04)
    GET_GRID    - reciprocal ETA-grid spacing: pn
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  3.97354603E-04 Max =  4.03391817E-04)
    GET_GRID    - x-location of RHO-points: x-rho
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  4.60475298E+06 Max =  5.66965824E+06)
    GET_GRID    - y-location of RHO-points: y_rho
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  1.14715526E+06 Max =  1.59984170E+06)
    GET_GRID    - x-location of PSI-points: x_psi
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  4.60603600E+06 Max =  5.66837522E+06)
    GET_GRID    - y-location of PSI-points: y-psi
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  1.14843404E+06 Max =  1.59856292E+06)
    GET_GRID    - x-location of U-points: x_u
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  4.60603600E+06 Max =  5.66837522E+06)
    GET_GRID    - y-location of U-points: y_u
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  1.14715526E+06 Max =  1.59984170E+06)
    GET_GRID    - x-location of V-points: x_v
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  4.60475298E+06 Max =  5.66965824E+06)
    GET_GRID    - y-location of V-points: y_v
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  1.14843404E+06 Max =  1.59856292E+06)
    GET_GRID    - latitude of RHO-points lat_rho
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  1.02613166E+01 Max =  1.42388533E+01)
    GET_GRID    - longitude of RHO-points: lon_rho
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  4.14115385E+01 Max =  5.09884615E+01)
    GET_GRID    - latitude of PSI-points lat_psi
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  1.02726327E+01 Max =  1.42277060E+01)
    GET_GRID    - longitude of PSI-points: lon_psi
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  4.14230769E+01 Max =  5.09769231E+01)
    GET_GRID    - latitude of U-points: lat_u
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  1.02613166E+01 Max =  1.42388533E+01)
    GET_GRID    - longitude of U-points: lon_u
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  4.14230769E+01 Max =  5.09769231E+01)
    GET_GRID    - latitude of V-points: lat_v
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  1.02726327E+01 Max =  1.42277060E+01)
    GET_GRID    - longitude of V-points: lon_v
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  4.14115385E+01 Max =  5.09884615E+01)
    GET_GRID    - mask on RHO-points: mask_rho
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on U-points: mask_u
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on V-points: mask_v
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - mask on PSI-points: mask_psi
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_GRID    - angle between XI-axis and EAST: angler
                   (Grid = 01, File: ../GRID_SETUP/grid.nc)
                   (Min = -1.46997540E-11 Max =  1.60250905E-11)

 Vertical S-coordinate System, Grid 01:

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

    32   0.0000000  -0.0000000       -0.000      -0.000      -0.000      -0.000
    31  -0.0312500  -0.0022549       -0.286      -1.675      -5.515      -8.425
    30  -0.0625000  -0.0048948       -0.573      -3.370     -11.478     -17.776
    29  -0.0937500  -0.0080749       -0.860      -5.091     -18.069     -28.425
    28  -0.1250000  -0.0119914       -1.147      -6.850     -25.516     -40.845
    27  -0.1562500  -0.0168927       -1.436      -8.657     -34.108     -55.632
    26  -0.1875000  -0.0230908       -1.726     -10.530     -44.207     -73.536
    25  -0.2187500  -0.0309708       -2.017     -12.486     -56.260     -95.483
    24  -0.2500000  -0.0409944       -2.310     -14.550     -70.806    -122.583
    23  -0.2812500  -0.0536924       -2.606     -16.747     -88.460    -156.112
    22  -0.3125000  -0.0696342       -2.904     -19.107    -109.885    -197.438
    21  -0.3437500  -0.0893672       -3.206     -21.656    -135.717    -247.878
    20  -0.3750000  -0.1133185       -3.512     -24.416    -166.452    -308.458
    19  -0.4062500  -0.1416675       -3.822     -27.396    -202.298    -379.609
    18  -0.4375000  -0.1742175       -4.136     -30.586    -243.028    -460.859
    17  -0.4687500  -0.2103218       -4.453     -33.954    -287.889    -550.653
    16  -0.5000000  -0.2489214       -4.772     -37.446    -335.651    -646.446
    15  -0.5312500  -0.2887177       -5.092     -40.998    -384.804    -745.115
    14  -0.5625000  -0.3284417       -5.412     -44.547    -433.873    -843.610
    13  -0.5937500  -0.3671228       -5.731     -48.044    -481.729    -939.598
    12  -0.6250000  -0.4042644       -6.049     -51.463    -527.796   -1031.885
    11  -0.6562500  -0.4398848       -6.366     -54.807    -572.095   -1120.516
    10  -0.6875000  -0.4744488       -6.681     -58.097    -615.166   -1206.608
     9  -0.7187500  -0.5087464       -6.997     -61.375    -657.927   -1292.059
     8  -0.7500000  -0.5437741       -7.313     -64.689    -701.537   -1379.265
     7  -0.7812500  -0.5806477       -7.630     -68.095    -747.293   -1470.908
     6  -0.8125000  -0.6205544       -7.951     -71.653    -796.574   -1569.842
     5  -0.8437500  -0.6647399       -8.275     -75.424    -850.828   -1679.062
     4  -0.8750000  -0.7145184       -8.604     -79.476    -911.584   -1801.727
     3  -0.9062500  -0.7712996       -8.940     -83.877    -980.480   -1941.225
     2  -0.9375000  -0.8366254       -9.283     -88.706   -1059.308   -2101.263
     1  -0.9687500  -0.9122123       -9.636     -94.048   -1150.063   -2285.967
     0  -1.0000000  -1.0000000      -10.000    -100.000   -1255.000   -2500.000

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

    Primary            Secondary            Accumulated to Current Step

  1-0.0006038999279304 0.0250000000000000-0.0006038999279304 0.0250000000000000
  2-0.0011053727845635 0.0250150974981983-0.0017092727124939 0.0500150974981983
  3-0.0015044206959555 0.0250427318178123-0.0032136934084494 0.0750578293160106
  4-0.0018010519301026 0.0250803423352112-0.0050147453385520 0.1001381716512218
  5-0.0019952879837948 0.0251253686334638-0.0070100333223467 0.1252635402846856
  6-0.0020871735042116 0.0251752508330587-0.0090972068265583 0.1504387911177443
  7-0.0020767890452586 0.0252274301706640-0.0111739958718169 0.1756662212884082
  8-0.0019642666586465 0.0252793498967954-0.0131382625304634 0.2009455711852036
  9-0.0017498083197106 0.0253284565632616-0.0148880708501740 0.2262740277484652
 10-0.0014337071879729 0.0253722017712543-0.0163217780381469 0.2516462295197196
 11-0.0010163717024450 0.0254080444509537-0.0173381497405919 0.2770542739706733
 12-0.0004983525116731 0.0254334537435148-0.0178365022522649 0.3024877277141881
 13 0.0001196277614762 0.0254459125563066-0.0177168744907887 0.3279336402704947
 14 0.0008366419202876 0.0254429218622697-0.0168802325705011 0.3533765621327644
 15 0.0016515227599278 0.0254220058142625-0.0152287098105733 0.3787985679470269
 16 0.0025628247984349 0.0253807177452643-0.0126658850121384 0.4041792856922913
 17 0.0035687831729657 0.0253166471253035-0.0090971018391727 0.4294959328175947
 18 0.0046672697013015 0.0252274275459793-0.0044298321378712 0.4547233603635741
 19 0.0058557461086128 0.0251107458034468 0.0014259139707416 0.4798341061670209
 20 0.0071312144194823 0.0249643521507315 0.0085571283902239 0.5047984583177524
 21 0.0084901645151857 0.0247860717902444 0.0170472929054096 0.5295845301079968
 22 0.0099285188562319 0.0245738176773648 0.0269758117616416 0.5541583477853615
 23 0.0114415743701613 0.0243256047059590 0.0384173861318028 0.5784839524913206
 24 0.0130239415046020 0.0240395653467049 0.0514413276364049 0.6025235178380255
 25 0.0146694804455855 0.0237139668090899 0.0661108080819904 0.6262374846471154
 26 0.0163712345011199 0.0233472297979502 0.0824820425831103 0.6495847144450656
 27 0.0181213606500220 0.0229379489354222 0.1006034032331323 0.6725226633804878
 28 0.0199110572560077 0.0224849149191717 0.1205144604891400 0.6950075782996595
 29 0.0217304889470413 0.0219871384877715 0.1422449494361813 0.7169947167874310
 30 0.0235687086599425 0.0214438762640955 0.1658136580961238 0.7384385930515265
 31 0.0254135768502523 0.0208546585475969 0.1912272349463761 0.7592932515991234
 32 0.0272516778673573 0.0202193191263406 0.2184789128137333 0.7795125707254640
 33 0.0290682334948725 0.0195380271796567 0.2475471463086058 0.7990505979051207
 34 0.0308470136562825 0.0188113213422849 0.2783941599648883 0.8178619192474056
 35 0.0325702442858409 0.0180401460008778 0.3109644042507291 0.8359020652482834
 36 0.0342185123647286 0.0172258898937318 0.3451829166154578 0.8531279551420151
 37 0.0357706681224704 0.0163704270846136 0.3809535847379282 0.8694983822266287
 38 0.0372037244036097 0.0154761603815518 0.4181573091415379 0.8849745426081804
 39 0.0384927531996420 0.0145460672714616 0.4566500623411799 0.8995206098796420
 40 0.0396107793462069 0.0135837484414705 0.4962608416873868 0.9131043583211125
 41 0.0405286713855386 0.0125934789578153 0.5367895130729254 0.9256978372789278
 42 0.0412150295941741 0.0115802621731769 0.5780045426670996 0.9372780994521047
 43 0.0416360711759210 0.0105498864333225 0.6196406138430206 0.9478279858854273
 44 0.0417555126200830 0.0095089846539245 0.6613961264631036 0.9573369705393517
 45 0.0415344492249439 0.0084650968384224 0.7029305756880475 0.9658020673777742
 46 0.0409312317865102 0.0074267356077988 0.7438618074745578 0.9732288029855730
 47 0.0399013404525124 0.0064034548131361 0.7837631479270702 0.9796322577987091
 48 0.0383972557416642 0.0054059213018233 0.8221604036687344 0.9850381791005324
 49 0.0363683267281804 0.0044459899082817 0.8585287303969148 0.9894841690088140
 50 0.0337606363915535 0.0035367817400771 0.8922893667884684 0.9930209507488911
 51 0.0305168641315885 0.0026927658302883 0.9228062309200569 0.9957137165791794
 52 0.0265761454486959 0.0019298442269986 0.9493823763687528 0.9976435608061780
 53 0.0218739287894439 0.0012654405907812 0.9712563051581966 0.9989090013969592
 54 0.0163418295573680 0.0007185923710451 0.9875981347155646 0.9996275937680043
 55 0.0099074812890402 0.0003100466321109 0.9975056160046047 0.9999376404001151
 56 0.0024943839953957 0.0000623595998849 1.0000000000000004 1.0000000000000000

 ndtfast, nfast =   40  56   nfast/ndtfast =  1.40000

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

    1.000000000000 1.039265571007 0.519632785503 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.21300

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

 Minimum X-grid spacing, DXmin =  2.48720494E+00 km    Water points =  2.48720494E+00 km
 Maximum X-grid spacing, DXmax =  2.52499433E+00 km    Water points =  2.52390409E+00 km
 Minimum Y-grid spacing, DYmin =  2.47897939E+00 km    Water points =  2.47897939E+00 km
 Maximum Y-grid spacing, DYmax =  2.51664381E+00 km    Water points =  2.51555717E+00 km
 Minimum Z-grid spacing, DZmin =  2.86140806E-01 m     Water points =  2.86140806E-01 m
 Maximum Z-grid spacing, DZmax =  2.14032856E+02 m     Water points =  2.14032856E+02 m

 Minimum barotropic Courant Number =  3.04615648E-02
 Maximum barotropic Courant Number =  4.81027948E-01
 Maximum Coriolis   Courant Number =  7.74611797E-03


 NLM: GET_STATE - Reading state initial conditions,                       2020-05-10 00:00:00.00
                   (Grid 01, t = 18980.0000, File: restart.nc, Rec=0001, Index=1)
                - free-surface
                   (Min =  1.08800694E-01 Max =  4.32979852E-01)
                - vertically integrated u-momentum component
                   (Min = -5.02590299E-01 Max =  2.85736531E-01)
                - vertically integrated v-momentum component
                   (Min = -3.90079319E-01 Max =  5.59199214E-01)
                - u-momentum component
                   (Min = -7.22452104E-01 Max =  8.92508447E-01)
                - v-momentum component
                   (Min = -1.01424241E+00 Max =  1.05941510E+00)
                - potential temperature
                   (Min =  2.84011436E+00 Max =  3.25317841E+01)
                - salinity
                   (Min =  3.48500099E+01 Max =  4.05849228E+01)
    GET_2DFLD   - surface u-wind component,                        2020-05-10 00:00:00.00
                   (Grid=01, Rec=0000001, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min = -6.64050770E+00 Max =  2.45534849E+00)      regrid = F
    GET_2DFLD   - surface v-wind component,                        2020-05-10 00:00:00.00
                   (Grid=01, Rec=0000001, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min = -3.60290575E+00 Max =  8.09830570E+00)      regrid = F
    GET_2DFLD   - surface air pressure,                            2020-05-10 00:00:00.00
                   (Grid=01, Rec=0000001, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  1.00588306E+03 Max =  1.01901312E+03)      regrid = F
    GET_2DFLD   - cloud fraction,                                  2020-05-10 00:00:00.00
                   (Grid=01, Rec=0000001, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  0.00000000E+00 Max =  9.99894500E-01)      regrid = F
    GET_2DFLD   - solar shortwave radiation flux,                  2020-05-10 00:00:00.00
                   (Grid=01, Rec=0000001, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)      regrid = F
    GET_2DFLD   - surface air temperature,                         2020-05-10 00:00:00.00
                   (Grid=01, Rec=0000001, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  1.07023497E+01 Max =  3.41326981E+01)      regrid = F
    GET_2DFLD   - surface air relative humidity,                   2020-05-10 00:00:00.00
                   (Grid=01, Rec=0000001, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  1.27292128E-01 Max =  9.57469864E-01)      regrid = F
    GET_2DFLD   - rain fall rate,                                  2020-05-10 00:00:00.00
                   (Grid=01, Rec=0000001, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  0.00000000E+00 Max =  6.76024100E-03)      regrid = F
    GET_NGFLD   - 2D u-momentum western boundary condition,        2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 2D v-momentum western boundary condition,        2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,        2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min = -8.86717141E-02 Max =  1.02867723E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,        2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min = -5.41214459E-02 Max =  1.39316723E-01)
    GET_NGFLD   - 2D u-momentum southern boundary condition,       2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 2D v-momentum southern boundary condition,       2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 2D u-momentum northern boundary condition,       2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min = -3.05563003E-01 Max =  9.68733653E-02)
    GET_NGFLD   - 2D v-momentum northern boundary condition,       2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min = -2.84521997E-01 Max =  3.47158015E-01)
    GET_NGFLD   - 3D u-momentum western boundary condition,        2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D v-momentum western boundary condition,        2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,        2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min = -5.81197202E-01 Max =  5.46317399E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,        2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min = -1.15826383E-01 Max =  8.33109260E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition,       2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D v-momentum southern boundary condition,       2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D u-momentum northern boundary condition,       2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min = -4.14090186E-01 Max =  4.77441907E-01)
    GET_NGFLD   - 3D v-momentum northern boundary condition,       2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min = -4.37167287E-01 Max =  6.27510071E-01)
    GET_NGFLD   - temperature western boundary condition,          2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  1.91904068E+01 Max =  3.05585861E+01)
    GET_NGFLD   - salinity western boundary condition,             2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  3.61785202E+01 Max =  3.73029060E+01)
    GET_NGFLD   - temperature eastern boundary condition,          2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  3.03624201E+00 Max =  3.04582500E+01)
    GET_NGFLD   - salinity eastern boundary condition,             2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  3.48864555E+01 Max =  3.66533089E+01)
    GET_NGFLD   - temperature southern boundary condition,         2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  2.72091293E+01 Max =  3.14871426E+01)
    GET_NGFLD   - salinity southern boundary condition,            2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  3.62095985E+01 Max =  3.64679527E+01)
    GET_NGFLD   - temperature northern boundary condition,         2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  3.00305533E+00 Max =  3.14966793E+01)
    GET_NGFLD   - salinity northern boundary condition,            2020-05-10 00:00:00.00
                   (Grid= 01, Rec=0000001, Index=2, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0000
                   (Min =  3.48864555E+01 Max =  4.05849228E+01)

 Basin information for Grid 01:

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

 Initial domain volumes:  TotVolume =  3.3897067099E+14 m3
                         MinCellVol =  1.8371097970E+06 m3
                         MaxCellVol =  1.3409912124E+09 m3
                            Max/Min =  7.2994614401E+02


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

    GET_2DFLD   - surface u-wind component,                        2020-05-10 01:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0417
                   (Min = -6.97507191E+00 Max =  3.25961447E+00)      regrid = F
    GET_2DFLD   - surface v-wind component,                        2020-05-10 01:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0417
                   (Min = -4.14099026E+00 Max =  8.80652809E+00)      regrid = F
    GET_2DFLD   - surface air pressure,                            2020-05-10 01:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0417
                   (Min =  1.00652405E+03 Max =  1.01928174E+03)      regrid = F
    GET_2DFLD   - cloud fraction,                                  2020-05-10 01:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0417
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)      regrid = F
    GET_2DFLD   - solar shortwave radiation flux,                  2020-05-10 00:30:00.00
                   (Grid=01, Rec=0000002, Index=2, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0208
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)      regrid = F
    GET_2DFLD   - surface air temperature,                         2020-05-10 01:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0417
                   (Min =  1.06005125E+01 Max =  3.34904709E+01)      regrid = F
    GET_2DFLD   - surface air relative humidity,                   2020-05-10 01:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0417
                   (Min =  1.59045172E-01 Max =  9.56996460E-01)      regrid = F
    GET_2DFLD   - rain fall rate,                                  2020-05-10 01:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0417
                   (Min = -1.66970437E-08 Max =  9.31945469E-05)      regrid = F
    GET_NGFLD   - 2D u-momentum western boundary condition,        2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 2D v-momentum western boundary condition,        2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,        2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min = -9.48527753E-02 Max =  9.83680189E-02)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,        2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min = -4.94436584E-02 Max =  1.42280862E-01)
    GET_NGFLD   - 2D u-momentum southern boundary condition,       2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 2D v-momentum southern boundary condition,       2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 2D u-momentum northern boundary condition,       2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min = -2.75550008E-01 Max =  9.55157056E-02)
    GET_NGFLD   - 2D v-momentum northern boundary condition,       2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min = -2.83331871E-01 Max =  2.58206248E-01)
    GET_NGFLD   - 3D u-momentum western boundary condition,        2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D v-momentum western boundary condition,        2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,        2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min = -5.62337637E-01 Max =  6.07338905E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,        2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min = -1.24728106E-01 Max =  7.97718227E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition,       2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D v-momentum southern boundary condition,       2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 3D u-momentum northern boundary condition,       2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min = -3.84849995E-01 Max =  6.07338905E-01)
    GET_NGFLD   - 3D v-momentum northern boundary condition,       2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min = -4.33289081E-01 Max =  7.15250075E-01)
    GET_NGFLD   - temperature western boundary condition,          2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  1.91842060E+01 Max =  3.07449093E+01)
    GET_NGFLD   - salinity western boundary condition,             2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  3.61758232E+01 Max =  3.72740593E+01)
    GET_NGFLD   - temperature eastern boundary condition,          2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  3.03468204E+00 Max =  3.04570637E+01)
    GET_NGFLD   - salinity eastern boundary condition,             2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  3.48864365E+01 Max =  3.66421242E+01)
    GET_NGFLD   - temperature southern boundary condition,         2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  2.74081535E+01 Max =  3.15430851E+01)
    GET_NGFLD   - salinity southern boundary condition,            2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  3.62207870E+01 Max =  3.64754257E+01)
    GET_NGFLD   - temperature northern boundary condition,         2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  3.00219035E+00 Max =  3.14611416E+01)
    GET_NGFLD   - salinity northern boundary condition,            2020-05-11 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18981.0000
                   (Min =  3.48864365E+01 Max =  4.05836525E+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

         0 2020-05-10 00:00:00.00  1.026137E-02  8.836775E+03  8.836786E+03  3.428660E+14
                     (040,166,04)  5.293244E-02  8.331740E-02  0.000000E+00  1.157315E+00
      DEF_HIS     - creating  history      file, Grid 01: history.nc
      WRT_HIS     - wrote history     fields (Index=1,1) in record = 0000001
         1 2020-05-10 00:03:36.00  1.024316E-02  8.836796E+03  8.836806E+03  3.428660E+14
                     (088,108,13)  3.099076E-02  3.725019E-02  3.367014E-01  1.159092E+00
         2 2020-05-10 00:07:12.00  1.024466E-02  8.836810E+03  8.836820E+03  3.428660E+14
                     (088,108,32)  5.618793E-02  5.638364E-02  1.297555E+00  1.157689E+00
         3 2020-05-10 00:10:48.00  1.024248E-02  8.836805E+03  8.836816E+03  3.428660E+14
                     (091,107,32)  2.408069E-02  3.946335E-03  1.302169E+00  1.166491E+00
         4 2020-05-10 00:14:24.00  1.024207E-02  8.836790E+03  8.836800E+03  3.428660E+14
                     (089,111,32)  2.423850E-03  3.212979E-04  1.061929E+00  1.182189E+00
         5 2020-05-10 00:18:00.00  1.024592E-02  8.836772E+03  8.836782E+03  3.428660E+14
                     (093,108,32)  1.513243E-02  3.594201E-02  8.981401E-01  1.197115E+00
         6 2020-05-10 00:21:36.00  1.025188E-02  8.836758E+03  8.836768E+03  3.428660E+14
                     (094,110,32)  5.811438E-03  2.282464E-02  1.434105E+00  1.180565E+00
         7 2020-05-10 00:25:12.00  1.025544E-02  8.836743E+03  8.836753E+03  3.428660E+14
                     (095,110,32)  4.593401E-03  1.736608E-02  2.469627E+00  1.191403E+00
         8 2020-05-10 00:28:48.00  1.026380E-02  8.836725E+03  8.836735E+03  3.428660E+14
                     (089,115,32)  1.700877E-03  9.763824E-04  2.818668E+00  1.205004E+00
    GET_2DFLD   - solar shortwave radiation flux,                  2020-05-10 01:30:00.00
                   (Grid=01, Rec=0000003, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0625
                   (Min =  5.08788546E-08 Max =  9.69528897E-06)      regrid = F
         9 2020-05-10 00:32:24.00  1.028225E-02  8.836708E+03  8.836719E+03  3.428660E+14
                     (099,109,32)  7.091034E-03  1.467469E-03  1.966619E+00  1.228202E+00
        10 2020-05-10 00:36:00.00  1.030198E-02  8.836705E+03  8.836715E+03  3.428660E+14
                     (089,115,32)  2.832330E-03  8.301570E-04  1.755426E+00  1.222235E+00
        11 2020-05-10 00:39:36.00  1.032081E-02  8.836721E+03  8.836732E+03  3.428660E+14
                     (099,109,32)  3.108256E-03  1.094453E-02  2.024630E+00  1.213574E+00
        12 2020-05-10 00:43:12.00  1.034507E-02  8.836746E+03  8.836756E+03  3.428660E+14
                     (030,151,32)  6.851071E-03  5.971151E-03  1.208092E+00  1.271596E+00
        13 2020-05-10 00:46:48.00  1.037875E-02  8.836759E+03  8.836769E+03  3.428660E+14
                     (104,107,32)  6.970671E-03  9.777058E-03  1.232426E+00  1.348566E+00
        14 2020-05-10 00:50:24.00  1.042582E-02  8.836753E+03  8.836763E+03  3.428660E+14
                     (411,152,11)  1.211805E-03  1.206539E-03  1.345674E+00  1.443083E+00
        15 2020-05-10 00:54:00.00  1.047916E-02  8.836733E+03  8.836743E+03  3.428660E+14
                     (411,152,11)  5.017455E-04  9.976604E-04  1.535400E+00  1.592034E+00
        16 2020-05-10 00:57:36.00  1.053735E-02  8.836707E+03  8.836718E+03  3.428660E+14
                     (411,152,11)  5.331820E-04  8.157884E-04  1.726126E+00  1.810891E+00
    GET_2DFLD   - surface u-wind component,                        2020-05-10 02:00:00.00
                   (Grid=01, Rec=0000003, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0833
                   (Min = -7.51540041E+00 Max =  3.59421420E+00)      regrid = F
    GET_2DFLD   - surface v-wind component,                        2020-05-10 02:00:00.00
                   (Grid=01, Rec=0000003, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0833
                   (Min = -4.45536852E+00 Max =  8.94251823E+00)      regrid = F
    GET_2DFLD   - surface air pressure,                            2020-05-10 02:00:00.00
                   (Grid=01, Rec=0000003, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0833
                   (Min =  1.00714557E+03 Max =  1.01988403E+03)      regrid = F
    GET_2DFLD   - cloud fraction,                                  2020-05-10 02:00:00.00
                   (Grid=01, Rec=0000003, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0833
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)      regrid = F
    GET_2DFLD   - surface air temperature,                         2020-05-10 02:00:00.00
                   (Grid=01, Rec=0000003, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0833
                   (Min =  1.03561621E+01 Max =  3.30257339E+01)      regrid = F
    GET_2DFLD   - surface air relative humidity,                   2020-05-10 02:00:00.00
                   (Grid=01, Rec=0000003, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0833
                   (Min =  1.82836895E-01 Max =  9.63333054E-01)      regrid = F
    GET_2DFLD   - rain fall rate,                                  2020-05-10 02:00:00.00
                   (Grid=01, Rec=0000003, Index=1, File: forcing.nc)
                   (Tmin=      18980.0000 Tmax=      18986.0000)      t =      18980.0833
                   (Min = -6.47586589E-08 Max =  9.31609102E-05)      regrid = F
        17 2020-05-10 01:01:12.00  1.060384E-02  8.836683E+03  8.836693E+03  3.428660E+14
                     (411,152,11)  1.562163E-03  6.246929E-04  1.913323E+00  2.069768E+00
        18 2020-05-10 01:04:48.00  1.067727E-02  8.836662E+03  8.836672E+03  3.428660E+14
                     (411,152,11)  1.528056E-04  1.303258E-04  2.090965E+00  2.374164E+00
        19 2020-05-10 01:08:24.00  1.077168E-02  8.836641E+03  8.836652E+03  3.428660E+14
                     (411,152,11)  6.730647E-03  7.023192E-05  2.252413E+00  2.716213E+00
        20 2020-05-10 01:12:00.00  1.102312E-02  8.836619E+03  8.836630E+03  3.428660E+14
                     (411,152,11)  7.819759E-03  1.367692E-03  2.400504E+00  5.703904E+00
        21 2020-05-10 01:15:36.00  1.410363E-02  8.836600E+03  8.836614E+03  3.428660E+14
                     (409,153,12)  4.705692E-04  1.525085E+00  2.390743E+00  2.505813E+01
 Found Error: 01   Line: 309      Source: ROMS/Nonlinear/main3d.F
 Found Error: 01   Line: 299      Source: ROMS/Drivers/nl_ocean.h

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

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

 All percentages are with respect to total time =          151.782

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

 Analytical header files used:

     ROMS/Functionals/ana_btflux.h

 MAIN: Abnormal termination: BLOWUP.
 REASON: MaxSpeed =  2.5E+01
At the moment we have no idea what is going wrong here.
Many thanks in advance.

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

Re: Blow up problem in ROMS 3.9

#2 Unread post by jcwarner »

a new feature in the latest ROMS is the ability to choose different advection schemes for each tracer. But you need to be careful to use the new ocean.in format.
from your output is see:

Tracer Advection Scheme: NLM
========================
Variable Grid Horizontal Vertical
--------- ---- ------------ ------------
temp 1 varinfo.dat ROMS/Utilit
salt 1 NDS/boundar y/distribut

Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Chapman Imp Chapman Imp Chapman Imp Chapman Imp
ubar 1 Rad + Nud Rad + Nud Rad + Nud Rad + Nud
....


Not sure if this is just a write screen buffer purge issue, or perhaps you dont have the correct Tracer advections chosen. try something like this:

! The user has the option of specifying the full Keyword or the first
! two letters, regardless if using uppercase or lowercase. If nested
! grids, specify values for each grid (see glossary below).
Hadvection == HSIMT \ ! temperature
HSIMT ! salinity

Vadvection == HSIMT \ ! temperature
HSIMT ! salinity

ebpbenke
Posts: 2
Joined: Tue Dec 03, 2019 7:56 pm
Location: EBP Deutschland GmbH

Re: Blow up problem in ROMS 3.9

#3 Unread post by ebpbenke »

Hello Mr. Warner,

that was the right advice. Many thanks to you.

I have added the advection configurations in my file roms.in as follows:

Code: Select all

   Hadvection == U3       \             ! temperature
                 U3                             ! salinity
   Vadvection == C4       \             ! temperature
                 C4                             ! salinity
! Adjoint-based algorithms can have different horizontal and schemes
! for active and inert tracers.
ad_Hadvection == U3       \           ! temperature
                 U3                             ! salinity

ad_Vadvection == C4       \           ! temperature
                 C4                             ! salinity
In my previuos file roms.in I haven't this configuration entries.

As second, I have updated the file varinfo.dat

With this changes the ROMS version 3.9 is running successfully.

A little bit strange is the logging output (special characters) of active processes:
Image
Attachments
roms_3_9_process_info.png

Post Reply