facing Error:ROMS/TOMS - Blows up ......... exit_flag: 1

[somanath]

Hi,
I am quite new to ROMS and from computer science domain.
I want to do some bench marking on our server (Linux OS,256Gb Memory,40 cores) for that I am trying to run it .

but while running I am getting Error like

ROMS/TOMS - Blows up ................ exit_flag: 1


MAIN: Abnormal termination: BLOWUP.


the input .in file configuration is like
Lm == 599 ! Number of I-direction INTERIOR RHO-points
Mm == 881 ! Number of J-direction INTERIOR RHO-points
N == 40 ! Number of vertical levels
and
NtileI == 4 ! I-direction partition
NtileJ == 10 ! J-direction partition

the .nc file grid dimension is
xi_rho = 601 ;
xi_u = 600 ;
xi_v = 601 ;
eta_rho = 883 ;
eta_u = 883 ;
eta_v = 882 ;
s_rho = 40 ;
s_w = 41 ;

The OutPut/Error I am getting is (mpiexec.hydra -np 40 ./oceanM ocean_india.in )



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

Vertical S-coordinate System, Grid 01:

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

40 0.0000000 0.0000000 0.000 0.000 0.000 0.000
39 -0.0250000 -0.0000309 0.000 0.000 0.000 0.000
38 -0.0500000 -0.0001247 0.000 0.000 0.000 0.000
37 -0.0750000 -0.0002842 0.000 0.000 0.000 0.000
36 -0.1000000 -0.0005144 0.000 0.000 0.000 0.000
35 -0.1250000 -0.0008222 0.000 0.000 0.000 0.000
34 -0.1500000 -0.0012170 0.000 0.000 0.000 0.000
33 -0.1750000 -0.0017111 0.000 0.000 0.000 0.000
32 -0.2000000 -0.0023196 0.000 0.000 0.000 0.000
31 -0.2250000 -0.0030610 0.000 0.000 0.000 0.000
30 -0.2500000 -0.0039582 0.000 0.000 0.000 0.000
29 -0.2750000 -0.0050385 0.000 0.000 0.000 0.000
28 -0.3000000 -0.0063350 0.000 0.000 0.000 0.000
27 -0.3250000 -0.0078873 0.000 0.000 0.000 0.000
26 -0.3500000 -0.0097429 0.000 0.000 0.000 0.000
25 -0.3750000 -0.0119584 0.000 0.000 0.000 0.000
24 -0.4000000 -0.0146013 0.000 0.000 0.000 0.000
23 -0.4250000 -0.0177522 0.000 0.000 0.000 0.000
22 -0.4500000 -0.0215067 0.000 0.000 0.000 0.000
21 -0.4750000 -0.0259789 0.000 0.000 0.000 0.000
20 -0.5000000 -0.0313040 0.000 0.000 0.000 0.000
19 -0.5250000 -0.0376427 0.000 0.000 0.000 0.000
18 -0.5500000 -0.0451859 0.000 0.000 0.000 0.000
17 -0.5750000 -0.0541596 0.000 0.000 0.000 0.000
16 -0.6000000 -0.0648314 0.000 0.000 0.000 0.000
15 -0.6250000 -0.0775180 0.000 0.000 0.000 0.000
14 -0.6500000 -0.0925932 0.000 0.000 0.000 0.000
13 -0.6750000 -0.1104976 0.000 0.000 0.000 0.000
12 -0.7000000 -0.1317496 0.000 0.000 0.000 0.000
11 -0.7250000 -0.1569575 0.000 0.000 0.000 0.000
10 -0.7500000 -0.1868327 0.000 0.000 0.000 0.000
9 -0.7750000 -0.2222045 0.000 0.000 0.000 0.000
8 -0.8000000 -0.2640351 0.000 0.000 0.000 0.000
7 -0.8250000 -0.3134348 0.000 0.000 0.000 0.000
6 -0.8500000 -0.3716761 0.000 0.000 0.000 0.000
5 -0.8750000 -0.4402058 0.000 0.000 0.000 0.000
4 -0.9000000 -0.5206513 0.000 0.000 0.000 0.000
3 -0.9250000 -0.6148197 0.000 0.000 0.000 0.000
2 -0.9500000 -0.7246838 0.000 0.000 0.000 0.000
1 -0.9750000 -0.8523498 0.000 0.000 0.000 0.000
0 -1.0000000 -1.0000000 0.000 0.000 0.000 0.000

Time Splitting Weights for Grid 01: ndtfast = 5 nfast = 7
==================================

Primary Secondary Accumulated to Current Step

1 0.0119306869089341 0.2000000000000000 0.0119306869089341 0.2000000000000000
2 0.0475049776008205 0.1976138626182131 0.0594356645097546 0.3976138626182131
3 0.1047629417620785 0.1881128670980491 0.1641986062718331 0.5857267297162622
4 0.1760826281731955 0.1671602787456334 0.3402812344450286 0.7528870084618955
5 0.2416469636634478 0.1319437531109943 0.5819281981084765 0.8848307615728898
6 0.2602974116474964 0.0836143603783047 0.8422256097559729 0.9684451219511945
7 0.1577743902440273 0.0315548780488055 1.0000000000000002 0.9999999999999999

ndtfast, nfast = 5 7 nfast/ndtfast = 1.40000

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

1.000000000000 1.084198606272 0.542099303136 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.00000

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

Minimum X-grid spacing, DXmin = 1.00000000E+00 km
Maximum X-grid spacing, DXmax = 1.00000000E+00 km
Minimum Y-grid spacing, DYmin = 1.00000000E+00 km
Maximum Y-grid spacing, DYmax = 1.00000000E+00 km
Minimum Z-grid spacing, DZmin = 0.00000000E+00 m
Maximum Z-grid spacing, DZmax = 0.00000000E+00 m

Minimum barotropic Courant Number = 0.00000000E+00
Maximum barotropic Courant Number = 0.00000000E+00
Maximum Coriolis Courant Number = 2.47800000E-03

GET_NGFLD - free-surface western boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -4.61830586E-01 Max = 0.00000000E+00)
GET_NGFLD - free-surface southern boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -2.39182904E-01 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum western boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -1.99090600E-01 Max = 6.18652590E-02)
GET_NGFLD - 2D v-momentum western boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -4.85758297E-02 Max = 3.84716615E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -1.44191995E-01 Max = 3.66550982E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -4.45211306E-02 Max = 3.96983512E-02)
GET_NGFLD - 3D u-momentum western boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -2.52713412E-01 Max = 2.72067100E-01)
GET_NGFLD - 3D v-momentum western boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -1.51217207E-01 Max = 1.91503897E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -2.91490614E-01 Max = 5.03728628E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -9.51866090E-01 Max = 4.69759792E-01)
GET_NGFLD - temperature western boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = 0.00000000E+00 Max = 3.02323895E+01)
GET_NGFLD - salinity western boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = 0.00000000E+00 Max = 3.65302200E+01)
GET_NGFLD - temperature southern boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = 0.00000000E+00 Max = 3.17947292E+01)
GET_NGFLD - salinity southern boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = 0.00000000E+00 Max = 3.63209915E+01)

Basin information for Grid 01:

Maximum grid stiffness ratios: rx0 = NaN (Beckmann and Haidvogel)
rx1 = NaN (Haney)

Initial basin volumes: TotVolume = 0.0000000000E+00 m3
MinVolume = 0.0000000000E+00 m3
MaxVolume = 0.0000000000E+00 m3
Max/Min = 0.0000000000E+00

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00028800)

GET_NGFLD - free-surface western boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -4.62206513E-01 Max = 0.00000000E+00)
GET_NGFLD - free-surface southern boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -2.40414307E-01 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum western boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -1.70969397E-01 Max = 6.63159266E-02)
GET_NGFLD - 2D v-momentum western boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -4.88123894E-02 Max = 3.76734100E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -2.32917905E-01 Max = 3.50951888E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -4.49662507E-02 Max = 1.17154799E-01)
GET_NGFLD - 3D u-momentum western boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -2.05577299E-01 Max = 2.85452396E-01)
GET_NGFLD - 3D v-momentum western boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -1.71727493E-01 Max = 2.02047005E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -4.01514113E-01 Max = 5.12598574E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = -7.88119018E-01 Max = 4.51970607E-01)
GET_NGFLD - temperature western boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = 0.00000000E+00 Max = 3.02420902E+01)
GET_NGFLD - salinity western boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = 0.00000000E+00 Max = 3.65283089E+01)
GET_NGFLD - temperature southern boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = 0.00000000E+00 Max = 3.18051491E+01)
GET_NGFLD - salinity southern boundary condition, t = 9136 12:00:00
(Rec=0000371, Index=1, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = 0.00000000E+00 Max = 3.63229103E+01)

STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed

0 9135 12:00:00 NaN NaN NaN 0.000000E+00
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 2.337052E-04

Elapsed CPU time (seconds):

Node # 0 CPU: 3.243
Node # 2 CPU: 3.103
Node # 4 CPU: 3.095
Node # 8 CPU: 3.111
Node # 10 CPU: 3.103
Node # 13 CPU: 3.087
Node # 14 CPU: 3.087
Node # 15 CPU: 3.087
Node # 17 CPU: 3.095
Node # 18 CPU: 3.088
Node # 23 CPU: 3.086
Node # 33 CPU: 3.098
Node # 36 CPU: 3.098
Node # 21 CPU: 3.087
Node # 25 CPU: 3.086
Node # 30 CPU: 3.086
Node # 35 CPU: 3.086
Node # 37 CPU: 3.087
Node # 19 CPU: 3.087
Node # 26 CPU: 3.085
Node # 28 CPU: 3.087
Node # 29 CPU: 3.086
Node # 32 CPU: 3.099
Node # 34 CPU: 3.099
Node # 38 CPU: 3.086
Node # 39 CPU: 3.083
Node # 27 CPU: 3.086
Node # 31 CPU: 3.087
Node # 1 CPU: 3.089
Node # 3 CPU: 3.088
Node # 9 CPU: 3.086
Node # 12 CPU: 3.087
Node # 16 CPU: 3.109
Node # 20 CPU: 3.087
Node # 24 CPU: 3.087
Node # 5 CPU: 3.088
Node # 7 CPU: 3.088
Node # 6 CPU: 3.088
Node # 11 CPU: 3.088
Node # 22 CPU: 3.086
Total: 123.768

Nonlinear model elapsed time profile, Grid: 01

Allocation and array initialization .............. 30.878 (24.9481 %)
Ocean state initialization ....................... 90.291 (72.9520 %)
Reading of input data ............................ 82.147 (66.3716 %)
Processing of input data ......................... 0.132 ( 0.1064 %)
Computation of global information integrals ...... 0.412 ( 0.3331 %)
2D/3D coupling, vertical metrics ................. 2.718 ( 2.1963 %)
Omega vertical velocity .......................... 0.313 ( 0.2526 %)
Equation of state for seawater ................... 0.946 ( 0.7647 %)
Total: 207.837 167.9249

Nonlinear model message Passage profile, Grid: 01

Message Passage: 2D halo exchanges ............... 2.960 ( 2.3917 %)
Message Passage: 3D halo exchanges ............... 0.918 ( 0.7416 %)
Message Passage: 4D halo exchanges ............... 0.124 ( 0.1000 %)
Message Passage: data broadcast .................. 79.873 (64.5346 %)
Message Passage: data reduction .................. 0.182 ( 0.1473 %)
Total: 84.057 67.9153

All percentages are with respect to total time = 123.768


ROMS/TOMS - Output NetCDF summary for Grid 01:

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_grid.h
ROMS/Functionals/ana_initial.h
ROMS/Functionals/ana_smflux.h
ROMS/Functionals/ana_stflux.h

ROMS/TOMS - Blows up ................ exit_flag: 1


MAIN: Abnormal termination: BLOWUP.


Please help me to solve this issue,
thanks in advance

[kate]

You have:

Code: Select all

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

40 0.0000000 0.0000000 0.000 0.000 0.000 0.000
39 -0.0250000 -0.0000309 0.000 0.000 0.000 0.000
38 -0.0500000 -0.0001247 0.000 0.000 0.000 0.000
37 -0.0750000 -0.0002842 0.000 0.000 0.000 0.000
36 -0.1000000 -0.0005144 0.000 0.000 0.000 0.000
35 -0.1250000 -0.0008222 0.000 0.000 0.000 0.000
34 -0.1500000 -0.0012170 0.000 0.000 0.000 0.000

You cannot have zero depth. How are you providing the depth?

[somanath]

Hi Kate,
Thanks for you quick reply...

I tried to figure out the issue but couldn't ,So I compiled a different version(3.7) but I was stuck with some other error like
GET_VARCOORDS - Cannot find "coordinates" attribute for variable: U10
in file: /home/somanath/INCUSE/input_files/U10.nc

This attribute is needed to interpolate input data
to model grid. Following CF compliance, we need:

float my_var(time, lat, lon) ;
my_var:long_name = "my variable long name" ;
my_var:units = "my variable units" ;
my_var:coordinates = "lon lat" ;
my_var:time = "my_var_time" ;


GET_2DFLD - error while reading variable: U10 at TIME index = 19

Elapsed CPU time (seconds):

Node # 0 CPU: 12.928

ROMS/TOMS - Output NetCDF summary for Grid 01:

ROMS/TOMS - Input error ............. exit_flag: 2


The U10.nc file is something like
netcdf U10 {
dimensions:
lon = 401 ;
lat = 281 ;
T_00Z = UNLIMITED ; // (121 currently)
variables:
double lon(lon) ;

The value of U10 in varinfo.dat

'U10' ! Input
'surface u-wind component'
'meter second-1' ! [m/s]
'u-wind, scalar, series'
'T_00Z'
'idUair'
'r2dvar'
1.0d0

The Dimension in cean_avg_9136.nc is

netcdf ocean_avg_9136 {
dimensions:
xi_rho = 601 ;
xi_u = 600 ;
xi_v = 601 ;
xi_psi = 600 ;
eta_rho = 883 ;
eta_u = 883 ;
eta_v = 882 ;
eta_psi = 882 ;
s_rho = 40 ;
s_w = 41 ;
tracer = 14 ;
boundary = 4 ;
ocean_time = UNLIMITED ; // (1 currently)

[kate]

I can't get to my winds files at the moment, but the ncdump of the U10 variable should look like:

Code: Select all

float U10(time, lat, lon) ;
U10:long_name = "my variable long name" ;
U10:units = "my variable units" ;
U10:coordinates = "lon lat" ;
U10:time = "my_var_time" ;

Can you show us the full "ncdump -h" of that file? I assume U10 is on some regular lat,lon grid and you want ROMS to interpolate it to your grid?

[somanath]

Hi kate,
Thanks a lot, It worked what I did I just added the value in variables part of U10.nc file as you mentioned.
But Now I am stuck with another error
INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

GET_NGFLD - O2 western boundary condition, t = 9135 12:00:00
(Rec=0000370, Index=2, File: roms_bry_1x48_2014_2015.nc)
(Tmin= 8766.5000 Tmax= 9495.5000)
(Min = 0.00000000E+00 Max = 2.08586807E+02)

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

Elapsed CPU time (seconds):


INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.
Node # 0 CPU: 15.129

ROMS/TOMS - Output NetCDF summary for Grid 01:

ROMS/TOMS - Input error ............. exit_flag: 2

Node # 3 CPU: 15.118
Node # 4 CPU: 15.118
Node # 5 CPU: 15.118
Node # 6 CPU: 15.117
Node # 8 CPU: 15.117
Node # 9 CPU: 15.117
Node # 10 CPU: 15.117
Node # 11 CPU: 15.117

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.
Node # 12 CPU: 15.117

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.
Node # 13 CPU: 15.117

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.
Node # 14 CPU: 15.117

INQUIRE - unable to find requested variable: NO3_east
in file:
file name is blank, cannot be determined.
Node # 15 CPU: 15.117

ERROR: Abnormal termination: NetCDF INPUT.
REASON: No error
Node # 1 CPU: 15.118
Node # 2 CPU: 15.118
Node # 7 CPU: 15.118

Please check this out

[kate]

The only thing that needs nitrate boundary conditions is one of the biology submodules. Did you turn one of them on?

[somanath]

yes, I am using the biology module bio_Fennel.in.

[wilkin]

Basin information for Grid 01:

Maximum grid stiffness ratios: rx0 = NaN (Beckmann and Haidvogel)
rx1 = NaN (Haney)

Initial basin volumes: TotVolume = 0.0000000000E+00 m3
MinVolume = 0.0000000000E+00 m3
MaxVolume = 0.0000000000E+00 m3
Max/Min = 0.0000000000E+00

Your ocean depth is zero everywhere.

ROMS reports that you are using an analytical grid from the ana_grid.h functional in your main source trunk:

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_grid.h
ROMS/Functionals/ana_initial.h
ROMS/Functionals/ana_smflux.h
ROMS/Functionals/ana_stflux.h

So the modifications you made to ana_grid.h have not set a finite positive bottom depth. You won't get far until you fix that.

[kate]

What boundary conditions are you providing for the ecosystem components? How are its LBC options set? Some values require external boundary fields. Do you have some you think it should be reading or do you not have them?

[rayltonbendo]

I am a new ROMS user.
I'm having problems with my case study, could someone help me?

--------------------------------------------------------------------------------
Model Input Parameters: ROMS/TOMS version 3.9
Tuesday - May 31, 2022 - 4:30:50 PM
--------------------------------------------------------------------------------

IBIS

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/gfortran
Compiler flags : -frepack-arrays -O3 -ffast-math -

SVN Root URL : https://www.myroms.org/svn/src/tags/roms-3.9
SVN Revision : 1085

Local Root : /home/joana/roms-3.9
Header Dir : /home/joana/Run1
Header file : ibis.h
Analytical Dir : /home/joana/Run1

Resolution, Grid 01: 254x110x30, Parallel Threads: 1, Tiling: 1x16


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

104832 ntimes Number of timesteps for 3-D equations.
1.000 dt Timestep size (s) for 3-D equations.
1 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.
288 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
1 nSTA Number of timesteps between the writing of data
the stations file.
T ldefout Switch to create a new output NetCDF file(s).
4320 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.
576 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
1 ntsDIA Starting timestep for the accumulation of output
time-averaged diagnostics data.
72 nDIA Number of timesteps between the writing of
time-averaged data into diagnostics file.
0.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
0.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
0.0000E+00 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
T LuvSponge Turning ON sponge on horizontal momentum.
T LtracerSponge(01) Turning ON sponge on tracer 01: temp
T LtracerSponge(02) Turning ON 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-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
5.0000E-06 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
5.0000E-06 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
1 lmd_Jwt Jerlov water type.
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
3.0000E+00 theta_s S-coordinate surface control parameter.
2.0000E+00 theta_b S-coordinate bottom control parameter.
50.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
18262.000 dstart Time-stamp assigned to model initialization (days).
19500101.00 time_ref Reference time for units attribute (yyyymmdd.dd)
3.6000E+02 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
3.6000E+02 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
3.6000E+02 Znudg Nudging/relaxation time scale (days)
for free-surface.
3.6000E+02 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
3.6000E+02 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
1.2000E+02 obcfac Factor between passive and active
open boundary conditions.
T VolCons(1) NLM western edge boundary volume conservation.
T VolCons(2) NLM southern edge boundary volume conservation.
T VolCons(3) NLM eastern edge boundary volume conservation.
T VolCons(4) NLM northern edge boundary volume conservation.
14.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
T LtracerSrc(01) Turning ON point Sources/Sinks on tracer 01: temp
T LtracerSrc(02) Turning ON point Sources/Sinks on tracer 02: salt
F LsshCLM Turning OFF processing of SSH climatology.
F Lm2CLM Turning OFF processing of 2D momentum climatology.
F Lm3CLM Turning OFF processing of 3D momentum climatology.
F LtracerCLM(01) Turning OFF processing of climatology tracer 01: temp
F LtracerCLM(02) Turning OFF processing of climatology tracer 02: salt
F LnudgeM2CLM Turning OFF nudging of 2D momentum climatology.
F LnudgeM3CLM Turning OFF nudging of 3D momentum climatology.
F LnudgeTCLM(01) Turning OFF nudging of climatology tracer 01: temp
F LnudgeTCLM(02) Turning OFF nudging of climatology tracer 02: salt

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 Aout(idFsur) Write out averaged free-surface.
T Aout(idUbar) Write out averaged 2D U-momentum component.
T Aout(idVbar) Write out averaged 2D V-momentum component.
T Aout(idUvel) Write out averaged 3D U-momentum component.
T Aout(idVvel) Write out averaged 3D V-momentum component.
T Aout(idWvel) Write out averaged W-momentum component.
T Aout(idOvel) Write out averaged omega vertical velocity.
T Aout(idTvar) Write out averaged tracer 01: temp
T Aout(idTvar) Write out averaged tracer 02: salt
T Aout(idTsur) Write out averaged surface net heat flux.
T Aout(idTsur) Write out averaged surface net salt flux.
T Aout(idSrad) Write out averaged shortwave radiation flux.


Output/Input Files:

Output Restart File: roms_rst_2000_delete.nc
Output History File: roms_his_2000_delete.nc
Output Averages File: roms_avg_2000_delete.nc
Output Diagnostics File: roms_dia_2000_delete.nc
Output Stations File: roms_sta_2000_delete.nc
Input Grid File: /home/joana/Raylton/run1/Grid/run1_grd1.nc
Input Nonlinear Initial File: /home/joana/Raylton/run1/CC_CI/cmems/initialization_conditions/ROMS_ini_run1.nc
Input Sources/Sinks File: /home/joana/Raylton/run1/CC_CI/rivers/ROMS_river_run1.nc
Input Forcing File 01: /home/joana/Raylton/run1/CC_CI/rivers/ROMS_river_run1.nc
Input Forcing File 02: /home/joana/Raylton/run1/CC_CI/era5/ecna124v_sms_era5_2000.nc
Input Forcing File 03: /home/joana/Raylton/run1/CC_CI/era5/ecna124v_shflux_era5_2000.nc
Input Forcing File 04: /home/joana/Raylton/run1/CC_CI/era5/ecna124v_swflux_era5_2000.nc
Input Forcing File 05: /home/joana/Raylton/run1/CC_CI/era5/ecna124v_swrad_era5_2000.nc
Input Lateral Boundary File 01: /home/joana/Raylton/run1/CC_CI/cmems/boundary_conditions/ROMS_BC_2000.nc
Station positions File: stations_ibis.in
ROMS I/O variables Metadata File: /home/joana/roms-3.9/ROMS/External/varinfo.dat

Tile partition information for Grid 01: 254x110x30 tiling: 1x16

tile Istr Iend Jstr Jend Npts

0 1 254 1 6 45720
1 1 254 7 13 53340
2 1 254 14 20 53340
3 1 254 21 27 53340
4 1 254 28 34 53340
5 1 254 35 41 53340
6 1 254 42 48 53340
7 1 254 49 55 53340
8 1 254 56 62 53340
9 1 254 63 69 53340
10 1 254 70 76 53340
11 1 254 77 83 53340
12 1 254 84 90 53340
13 1 254 91 97 53340
14 1 254 98 104 53340
15 1 254 105 110 45720

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 255.50 0.50 112.50 RHO-points
1 0.50 255.50 -0.50 112.50 RHO-points
2 0.50 255.50 -0.50 112.50 RHO-points
3 0.50 255.50 -0.50 112.50 RHO-points
4 0.50 255.50 -0.50 112.50 RHO-points
5 0.50 255.50 -0.50 112.50 RHO-points
6 0.50 255.50 -0.50 112.50 RHO-points
7 0.50 255.50 -0.50 112.50 RHO-points
8 0.50 255.50 -0.50 112.50 RHO-points
9 0.50 255.50 -0.50 112.50 RHO-points
10 0.50 255.50 -0.50 112.50 RHO-points
11 0.50 255.50 -0.50 112.50 RHO-points
12 0.50 255.50 -0.50 112.50 RHO-points
13 0.50 255.50 -0.50 112.50 RHO-points
14 0.50 255.50 -0.50 112.50 RHO-points
15 0.50 255.50 -0.50 111.50 RHO-points

0 0.00 255.00 0.50 112.50 U-points
1 0.00 255.00 -0.50 112.50 U-points
2 0.00 255.00 -0.50 112.50 U-points
3 0.00 255.00 -0.50 112.50 U-points
4 0.00 255.00 -0.50 112.50 U-points
5 0.00 255.00 -0.50 112.50 U-points
6 0.00 255.00 -0.50 112.50 U-points
7 0.00 255.00 -0.50 112.50 U-points
8 0.00 255.00 -0.50 112.50 U-points
9 0.00 255.00 -0.50 112.50 U-points
10 0.00 255.00 -0.50 112.50 U-points
11 0.00 255.00 -0.50 112.50 U-points
12 0.00 255.00 -0.50 112.50 U-points
13 0.00 255.00 -0.50 112.50 U-points
14 0.00 255.00 -0.50 112.50 U-points
15 0.00 255.00 -0.50 111.50 U-points

0 0.50 255.50 0.00 112.50 V-points
1 0.50 255.50 -0.50 112.50 V-points
2 0.50 255.50 -0.50 112.50 V-points
3 0.50 255.50 -0.50 112.50 V-points
4 0.50 255.50 -0.50 112.50 V-points
5 0.50 255.50 -0.50 112.50 V-points
6 0.50 255.50 -0.50 112.50 V-points
7 0.50 255.50 -0.50 112.50 V-points
8 0.50 255.50 -0.50 112.50 V-points
9 0.50 255.50 -0.50 112.50 V-points
10 0.50 255.50 -0.50 112.50 V-points
11 0.50 255.50 -0.50 112.50 V-points
12 0.50 255.50 -0.50 112.50 V-points
13 0.50 255.50 -0.50 112.50 V-points
14 0.50 255.50 -0.50 112.50 V-points
15 0.50 255.50 -0.50 111.00 V-points


Stations Parameters, Grid: 01
=============================

3 Nstation Number of stations to write out into stations file.
T Sout(idFsur) Write out free-surface.
T Sout(idUbar) Write out 2D U-momentum component.
T Sout(idVbar) Write out 2D V-momentum component.
T Sout(idu2dE) Write out 2D U-eastward at RHO-points.
T Sout(idv2dN) Write out 2D V-northward at RHO-points.
T Sout(idUvel) Write out 3D U-momentum component.
T Sout(idVvel) Write out 3D V-momentum component.
T Sout(idu3dE) Write out 3D U-eastward at RHO-points.
T Sout(idv3dN) Write out 3D V-northward at RHO-points.
T Sout(idWvel) Write out W-momentum component.
T Sout(idOvel) Write out omega vertical velocity.
T Sout(idTvar) Write out tracer 01: temp
T Sout(idTvar) Write out tracer 02: salt
T Sout(idUsms) Write out surface U-momentum stress.
T Sout(idVsms) Write out surface V-momentum stress.
T Sout(idUbms) Write out bottom U-momentum stress.
T Sout(idVbms) Write out bottom V-momentum stress.
T Sout(idTsur) Write out surface net heat flux.
T Sout(idTsur) Write out surface net salt flux.
T Sout(idSrad) Write out shortwave radiation flux.
T Sout(idDano) Write out density anomaly.
T Sout(idVvis) Write out vertical viscosity coefficient.
T Sout(idTdif) Write out vertical T-diffusion coefficient.
T Sout(idSdif) Write out vertical S-diffusion coefficient.
T Sout(idMtke) Write out turbulent kinetic energy.
T Sout(idMtls) Write out turbulent generic length-scale.

Flag and positions for station 0001: 1 -9.3500 43.8333
Flag and positions for station 0002: 1 -9.6500 42.2167
Flag and positions for station 0003: 2 50.0000 225.0000

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 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 Gradient Gradient Gradient Gradient

Activated C-preprocessing Options:

IBIS IBIS
ANA_BSFLUX Analytical kinematic bottom salinity flux
ANA_BTFLUX Analytical kinematic bottom temperature flux
ASSUMED_SHAPE Using assumed-shape arrays
AVERAGES Writing out time-averaged nonlinear model fields
DIAGNOSTICS_TS Computing and writing tracer diagnostic terms
DIAGNOSTICS_UV Computing and writing momentum diagnostic terms
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002)
DOUBLE_PRECISION Double precision arithmetic numerical kernel.
KANTHA_CLAYSON Kantha and Clayson stability function formulation
MASKING Land/Sea masking
MIX_GEO_TS Mixing of tracers along geopotential surfaces
MIX_GEO_UV Mixing of momentum along geopotential surfaces
MY25_MIXING Mellor/Yamada Level-2.5 mixing closure
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
!RST_SINGLE Double precision fields in restart NetCDF file
SALINITY Using salinity
SOLAR_SOURCE Solar Radiation Source Term
SOLVE3D Solving 3D Primitive Equations
STATIONS Writing out station data
TS_MPDATA_LIMIT Further limiter in upwind corrector fluxes for stability
TS_DIF2 Harmonic mixing of tracers
UV_ADV Advection of momentum
UV_COR Coriolis term
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum
UV_C4VADVECTION Fourth-order centered vertical advection of momentum
UV_QDRAG Quadratic bottom stress
UV_VIS2 Harmonic mixing of momentum
VAR_RHO_2D Variable density barotropic mode

Process Information:

Thread # 0 (pid= 1757994) is active.

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

GET_GRID - angle between XI-axis and EAST: angler
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.60268702E+00 Max = -1.53890563E+00)
GET_GRID - bathymetry at RHO-points: h
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.70095704E+03 Max = 5.27800395E+03)
GET_GRID - Coriolis parameter at RHO-points: f
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 8.40525343E-05 Max = 1.04707691E-04)
GET_GRID - reciprocal XI-grid spacing: pm
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 2.14445144E-04 Max = 2.14445144E-04)
GET_GRID - reciprocal ETA-grid spacing: pn
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 2.39808472E-04 Max = 2.81540215E-04)
GET_GRID - longitude of RHO-points: lon_rho
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.21700238E+01 Max = -7.07652660E+00)
GET_GRID - longitude of U-points: lon_u
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.21700238E+01 Max = -7.07652660E+00)
GET_GRID - longitude of V-points: lon_v
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.21470801E+01 Max = -7.09947028E+00)
GET_GRID - longitude of PSI-points: lon_psi
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.21470801E+01 Max = -7.09947028E+00)
GET_GRID - latitude of RHO-points lat_rho
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 3.51930988E+01 Max = 4.58865384E+01)
GET_GRID - latitude of U-points: lat_u
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 3.52140691E+01 Max = 4.58655748E+01)
GET_GRID - latitude of V-points: lat_v
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 3.51930988E+01 Max = 4.58865384E+01)
GET_GRID - latitude of PSI-points lat_psi
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 3.52140691E+01 Max = 4.58655748E+01)
GET_GRID - mask on RHO-points: mask_rho
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on U-points: mask_u
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on V-points: mask_v
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on PSI-points: mask_psi
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - horizontal viscosity sponge factor: visc_factor
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 1.00000000E+00 Max = 3.00000000E+00)
GET_GRID - horizontal diffusivity sponge factor: diff_factor
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 1.00000000E+00 Max = 3.00000000E+00)

Vertical S-coordinate System, Grid 01:

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

30 0.0000000 0.0000000 0.000 0.000 0.000 0.000
29 -0.0333333 -0.0012758 0.519 -0.865 -3.841 -8.321
28 -0.0666667 -0.0051074 5.516 -1.794 -12.129 -30.006
27 -0.1000000 -0.0115076 15.015 -2.788 -24.886 -65.120
26 -0.1333333 -0.0204966 29.051 -3.846 -42.147 -113.770
25 -0.1666667 -0.0321017 47.672 -4.969 -63.960 -176.098
24 -0.2000000 -0.0463550 70.933 -6.159 -90.380 -252.272
23 -0.2333333 -0.0632923 98.898 -7.416 -121.471 -342.479
22 -0.2666667 -0.0829499 131.630 -8.740 -157.294 -446.909
21 -0.3000000 -0.1053620 169.190 -10.134 -197.910 -565.742
20 -0.3333333 -0.1305570 211.626 -11.597 -243.367 -699.124
19 -0.3666667 -0.1585533 258.971 -13.130 -293.699 -847.153
18 -0.4000000 -0.1893544 311.232 -14.734 -348.911 -1009.847
17 -0.4333333 -0.2229439 368.380 -16.407 -408.974 -1187.119
16 -0.4666667 -0.2592789 430.339 -18.149 -473.814 -1378.747
15 -0.5000000 -0.2982845 496.978 -19.957 -543.300 -1584.338
14 -0.5333333 -0.3398462 568.096 -21.829 -617.234 -1803.293
13 -0.5666667 -0.3838043 643.415 -23.762 -695.337 -2034.778
12 -0.6000000 -0.4299468 722.561 -25.749 -777.241 -2277.684
11 -0.6333333 -0.4780043 805.063 -27.783 -862.477 -2530.602
10 -0.6666667 -0.5276456 890.341 -29.858 -950.469 -2791.802
9 -0.7000000 -0.5784757 977.702 -31.962 -1040.528 -3059.216
8 -0.7333333 -0.6300360 1066.343 -34.084 -1131.859 -3330.449
7 -0.7666667 -0.6818081 1155.355 -36.212 -1223.557 -3602.789
6 -0.8000000 -0.7332223 1243.740 -38.331 -1314.633 -3873.258
5 -0.8333333 -0.7836698 1330.430 -40.425 -1404.028 -4138.672
4 -0.8666667 -0.8325209 1414.323 -42.480 -1490.644 -4395.740
3 -0.9000000 -0.8791471 1494.317 -44.479 -1573.389 -4641.175
2 -0.9333333 -0.9229486 1569.361 -46.407 -1651.220 -4871.841
1 -0.9666667 -0.9633838 1638.505 -48.251 -1723.194 -5084.906
0 -1.0000000 -1.0000000 1700.957 -50.000 -1788.523 -5278.004

Time Splitting Weights for Grid 01: ndtfast = 1 nfast = 1
==================================

Primary Secondary Accumulated to Current Step

1 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000

ndtfast, nfast = 1 1 nfast/ndtfast = 1.00000

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

1.000000000000 1.000000000000 0.500000000000 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.00000

WARNING: unstable weights, reduce parameter Fgamma in mod_scalars.F


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

Minimum X-grid spacing, DXmin = 4.66319722E+00 km Water points = 4.66319722E+00 km
Maximum X-grid spacing, DXmax = 4.66319723E+00 km Water points = 4.66319723E+00 km
Minimum Y-grid spacing, DYmin = 3.55189044E+00 km Water points = 3.55189044E+00 km
Maximum Y-grid spacing, DYmax = 4.16999445E+00 km Water points = 4.16999445E+00 km
Minimum Z-grid spacing, DZmin = -1.31577343E+04 m Water points = 6.68401986E-01 m
Maximum Z-grid spacing, DZmax = 2.28713544E+04 m Water points = 2.72340269E+02 m

Minimum barotropic Courant Number = 5.93191142E-03
Maximum barotropic Courant Number = 7.74959916E-02
Maximum Coriolis Courant Number = 1.04707691E-04


NLM: GET_STATE - Reading state initial conditions, 2000-01-01 12:00:00.00
(Grid 01, t = 18262.5000, File: ROMS_ini_run1.nc, Rec=0001, Index=1)
- free-surface
(Min = -4.32987267E-01 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = -6.47848268E-02 Max = 1.09843572E-01)
- vertically integrated v-momentum component
(Min = -1.75664142E-01 Max = 1.38519824E-01)
- u-momentum component
(Min = -1.13222317E-01 Max = 1.28937116E-01)
- v-momentum component
(Min = -1.19600629E-01 Max = 1.26414030E-01)
- potential temperature
(Min = 0.00000000E+00 Max = 1.72034197E+01)
- salinity
(Min = 0.00000000E+00 Max = 3.64725254E+01)
GET_NGFLD - river runoff XI-positions at RHO-points
(Grid = 01, Min = 6.40000000E+01 Max = 1.04000000E+02)
GET_NGFLD - river runoff ETA-positions at RHO-points
(Grid = 01, Min = 9.60000000E+01 Max = 2.13000000E+02)
GET_NGFLD - river runoff direction
(Grid = 01, Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_NGFLD - river runoff mass transport vertical profile
(Grid = 01, Min = 0.00000000E+00 Max = 2.00000000E-01)
GET_NGFLD - river runoff potential temperature, 2000-01-01 00:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_river_run1.nc)
(Tmin= 18262.0000 Tmax= 73042.0000) t = 18262.0000
(Min = 1.47940000E+01 Max = 1.99500000E+01)
GET_NGFLD - river runoff salinity, 2000-01-01 00:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_river_run1.nc)
(Tmin= 18262.0000 Tmax= 73042.0000) t = 18262.0000
(Min = 1.57000000E+01 Max = 3.55000000E+01)
GET_2DFLD - surface u-momentum stress, 2000-01-01 11:00:00.00
(Grid=01, Rec=12, Index=1, File: ecna124v_sms_era5_2000.nc)
(Tmin= 18262.0000 Tmax= 18627.0000) t = 18262.4583
(Min = -2.62452790E-04 Max = 3.97277247E-05) regrid = T
GET_2DFLD - surface v-momentum stress, 2000-01-01 11:00:00.00
(Grid=01, Rec=12, Index=1, File: ecna124v_sms_era5_2000.nc)
(Tmin= 18262.0000 Tmax= 18627.0000) t = 18262.4583
(Min = -1.18589529E-04 Max = 1.04898793E-04) regrid = T
GET_2DFLD - solar shortwave radiation flux, 2000-01-01 11:30:00.00
(Grid=01, Rec=13, Index=1, File: ecna124v_swrad_era5_2000.nc)
(Tmin= 18261.9792 Tmax= 18626.9792) t = 18262.4792
(Min = 1.73880850E-05 Max = 1.16621777E-04) regrid = T
GET_2DFLD - surface net heat flux, 2000-01-01 11:30:00.00
(Grid=01, Rec=13, Index=1, File: ecna124v_shflux_era5_2000.nc)
(Tmin= 18261.9792 Tmax= 18626.9792) t = 18262.4792
(Min = 7.11119790E-06 Max = 5.67094035E-05) regrid = T
GET_NGFLD - 2D u-momentum western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -6.53617531E-02 Max = 9.55667123E-02)
GET_NGFLD - 2D v-momentum western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -5.20291477E-02 Max = 6.95130751E-02)
GET_NGFLD - 2D u-momentum eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -1.12437271E-01 Max = 7.48302564E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -6.42233565E-02 Max = 2.69341581E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -8.88442695E-02 Max = 1.01377249E-01)
GET_NGFLD - 2D v-momentum southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -1.06768280E-01 Max = 1.10900223E-01)
GET_NGFLD - 2D u-momentum northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -6.46682158E-02 Max = 9.21968892E-02)
GET_NGFLD - 2D v-momentum northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -1.75664142E-01 Max = 9.94500220E-02)
GET_NGFLD - 3D u-momentum western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -1.68131411E-01 Max = 1.81981206E-01)
GET_NGFLD - 3D v-momentum western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -1.47421896E-01 Max = 1.38651356E-01)
GET_NGFLD - 3D u-momentum eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -2.77225375E-01 Max = 3.22363853E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -2.80117363E-01 Max = 8.99794772E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -2.27593422E-01 Max = 2.02705994E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -2.47028366E-01 Max = 2.25783244E-01)
GET_NGFLD - 3D u-momentum northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -2.34356746E-01 Max = 1.43075764E-01)
GET_NGFLD - 3D v-momentum northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -4.88355607E-01 Max = 2.02198505E-01)
GET_NGFLD - temperature western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 1.98817396E+00 Max = 1.30999126E+01)
GET_NGFLD - salinity western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 3.48995514E+01 Max = 3.60373383E+01)
GET_NGFLD - temperature eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 2.04223514E+00 Max = 1.81060047E+01)
GET_NGFLD - salinity eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 3.49053917E+01 Max = 3.64541130E+01)
GET_NGFLD - temperature southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 1.96197915E+00 Max = 1.77900620E+01)
GET_NGFLD - salinity southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 3.48929863E+01 Max = 3.63163452E+01)
GET_NGFLD - temperature northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 0.00000000E+00 Max = 1.72163239E+01)
GET_NGFLD - salinity northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 0.00000000E+00 Max = 3.66767006E+01)

Basin information for Grid 01:

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

Initial domain volumes: TotVolume = 1.2742752400E+15 m3
MinCellVol = 1.2173304024E+07 m3
MaxCellVol = 5.0240679525E+09 m3
Max/Min = 4.1271194269E+02


NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000000043201 - 000000148032)

GET_NGFLD - river runoff potential temperature, 2000-01-31 00:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_river_run1.nc)
(Tmin= 18262.0000 Tmax= 73042.0000) t = 18292.0000
(Min = 1.47940000E+01 Max = 1.99500000E+01)
GET_NGFLD - river runoff salinity, 2000-01-31 00:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_river_run1.nc)
(Tmin= 18262.0000 Tmax= 73042.0000) t = 18292.0000
(Min = 1.57000000E+01 Max = 3.55000000E+01)
GET_2DFLD - surface u-momentum stress, 2000-01-01 12:00:00.00
(Grid=01, Rec=13, Index=2, File: ecna124v_sms_era5_2000.nc)
(Tmin= 18262.0000 Tmax= 18627.0000) t = 18262.5000
(Min = -2.65224179E-04 Max = 7.47214743E-05) regrid = T
GET_2DFLD - surface v-momentum stress, 2000-01-01 12:00:00.00
(Grid=01, Rec=13, Index=2, File: ecna124v_sms_era5_2000.nc)
(Tmin= 18262.0000 Tmax= 18627.0000) t = 18262.5000
(Min = -1.25773325E-04 Max = 1.30807716E-04) regrid = T
GET_2DFLD - solar shortwave radiation flux, 2000-01-01 12:30:00.00
(Grid=01, Rec=14, Index=2, File: ecna124v_swrad_era5_2000.nc)
(Tmin= 18261.9792 Tmax= 18626.9792) t = 18262.5208
(Min = 2.46481305E-05 Max = 1.24405759E-04) regrid = T
GET_2DFLD - surface net heat flux, 2000-01-01 12:30:00.00
(Grid=01, Rec=14, Index=2, File: ecna124v_shflux_era5_2000.nc)
(Tmin= 18261.9792 Tmax= 18626.9792) t = 18262.5208
(Min = 1.23045191E-05 Max = 7.02479404E-05) regrid = T
GET_NGFLD - 2D u-momentum western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -6.49111196E-02 Max = 9.59445834E-02)
GET_NGFLD - 2D v-momentum western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -4.95847687E-02 Max = 7.01998323E-02)
GET_NGFLD - 2D u-momentum eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -1.13728225E-01 Max = 6.92293793E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -6.42820746E-02 Max = 2.59706806E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -9.05174986E-02 Max = 9.94924977E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -1.05292454E-01 Max = 1.07233137E-01)
GET_NGFLD - 2D u-momentum northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -6.23283498E-02 Max = 7.24404529E-02)
GET_NGFLD - 2D v-momentum northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -2.05803722E-01 Max = 7.54545256E-02)
GET_NGFLD - 3D u-momentum western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -1.68550178E-01 Max = 1.80611610E-01)
GET_NGFLD - 3D v-momentum western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -1.41561508E-01 Max = 1.36841610E-01)
GET_NGFLD - 3D u-momentum eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -2.47794464E-01 Max = 3.19678664E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -3.25482637E-01 Max = 1.16334975E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -2.30487600E-01 Max = 2.05312297E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -2.35349774E-01 Max = 2.29068324E-01)
GET_NGFLD - 3D u-momentum northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -2.58162826E-01 Max = 1.40312970E-01)
GET_NGFLD - 3D v-momentum northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -5.30910313E-01 Max = 1.95586026E-01)
GET_NGFLD - temperature western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 1.98840737E+00 Max = 1.30311937E+01)
GET_NGFLD - salinity western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 3.48996239E+01 Max = 3.60337982E+01)
GET_NGFLD - temperature eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 2.04223323E+00 Max = 1.79368706E+01)
GET_NGFLD - salinity eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 3.49055405E+01 Max = 3.64442520E+01)
GET_NGFLD - temperature southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 1.96234500E+00 Max = 1.74535332E+01)
GET_NGFLD - salinity southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 3.48929443E+01 Max = 3.63216476E+01)
GET_NGFLD - temperature northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 0.00000000E+00 Max = 1.72353230E+01)
GET_NGFLD - salinity northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 0.00000000E+00 Max = 3.66278343E+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

43200 2000-01-01 12:00:00.00 NaN 2.088449E+04 NaN 1.237148E+15
(229,110,05) 4.499688E-05 1.104512E-04 0.000000E+00 4.732726E-01
Found Error: 01 Line: 314 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 302 Source: ROMS/Drivers/nl_ocean.h, ROMS_run

Elapsed wall CPU time for each process (seconds):

Thread # 0 CPU: 0.492
Total: 0.492

Nonlinear model elapsed CPU time profile, Grid: 01

Allocation and array initialization .............. 0.318 (64.7097 %)
Ocean state initialization ....................... 0.015 ( 2.9571 %)
Reading of input data ............................ 0.085 (17.1969 %)
Processing of input data ......................... 0.000 ( 0.0779 %)
Computation of global information integrals ...... 0.004 ( 0.8234 %)
2D/3D coupling, vertical metrics ................. 0.015 ( 2.9759 %)
Omega vertical velocity .......................... 0.003 ( 0.5809 %)
Equation of state for seawater ................... 0.021 ( 4.1928 %)
Total: 0.460 93.5147 %

Unique kernel(s) regions profiled ................ 0.460 93.5147 %
Residual, non-profiled code ...................... 0.032 6.4853 %


All percentages are with respect to total time = 0.492


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

Dynamic and Automatic memory (MB) usage for Grid 01: 254x110x30 tiling: 1x16

tile Dynamic Automatic USAGE

0 1134.84 14.22 1149.06
1 0.00 14.22 14.22
2 0.00 14.22 14.22
3 0.00 14.22 14.22
4 0.00 14.22 14.22
5 0.00 14.22 14.22
6 0.00 14.22 14.22
7 0.00 14.22 14.22
8 0.00 14.22 14.22
9 0.00 14.22 14.22
10 0.00 14.22 14.22
11 0.00 14.22 14.22
12 0.00 14.22 14.22
13 0.00 14.22 14.22
14 0.00 14.22 14.22
15 0.00 14.22 14.22

TOTAL 1134.84 227.54 1362.38

<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
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
Found Error: 01 Line: 361 Source: ROMS/Utility/close_io.F, close_out

ROMS/TOMS - Blows up ................ exit_flag: 1

[wilkin]

Looks like your initial conditions have a NaN in them.

STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed
0 9135 12:00:00 NaN NaN NaN 0.000000E+00

You have negative DZ in the land

Minimum Z-grid spacing, DZmin = -1.31577343E+04 m Water points = 6.68401986E-01 m

but no #define WET_DRY so you might have been unlucky with a divide by zero where h+zeta = 0.
Set the water depth to some finite positive value in the land, say h = 5, to prevent zero divisor.

Make sure there no NaNs in your initial velocity. Don't use _FillValue in the initial conditions file. ROMS wants finite values everywhere, including in the land.

But here's the thing:

Time Splitting Weights for Grid 01: ndtfast = 1 nfast = 1
==================================
Primary Secondary Accumulated to Current Step
1 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
ndtfast, nfast = 1 1 nfast/ndtfast = 1.00000
Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):
1.000000000000 1.000000000000 0.500000000000 1.000000000000 1.000000000000
Power filter parameters, Fgamma, gamma = 0.28400 0.00000
WARNING: unstable weights, reduce parameter Fgamma in mod_scalars.F

In roms.in, DT sets the baroclinic time step and this is divided by NDTFAST to get the barotropic time step. We recommend NDTFAST =~ 30 for typical realistic applications.
NDTFAST = 1 is wrong.

[rayltonbendo]

Tank's for your attention!

I have tried solving that problem and rerun the simulation, but it shows another error:

--------------------------------------------------------------------------------
Model Input Parameters: ROMS/TOMS version 3.9
Thursday - June 2, 2022 - 12:07:36 PM
--------------------------------------------------------------------------------

IBIS

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/gfortran
Compiler flags : -frepack-arrays -O3 -ffast-math -

SVN Root URL : https://www.myroms.org/svn/src/tags/roms-3.9
SVN Revision : 1085

Local Root : /home/joana/roms-3.9
Header Dir : /home/joana/Run1
Header file : ibis.h
Analytical Dir : /home/joana/Run1

Resolution, Grid 01: 254x110x30, Parallel Threads: 1, Tiling: 1x16


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

104832 ntimes Number of timesteps for 3-D equations.
300.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.
288 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
1 nSTA Number of timesteps between the writing of data
the stations file.
T ldefout Switch to create a new output NetCDF file(s).
4320 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.
576 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
1 ntsDIA Starting timestep for the accumulation of output
time-averaged diagnostics data.
72 nDIA Number of timesteps between the writing of
time-averaged data into diagnostics file.
0.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
0.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
0.0000E+00 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
T LuvSponge Turning ON sponge on horizontal momentum.
T LtracerSponge(01) Turning ON sponge on tracer 01: temp
T LtracerSponge(02) Turning ON 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-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
5.0000E-06 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
5.0000E-06 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
1 lmd_Jwt Jerlov water type.
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
3.0000E+00 theta_s S-coordinate surface control parameter.
2.0000E+00 theta_b S-coordinate bottom control parameter.
50.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
18262.000 dstart Time-stamp assigned to model initialization (days).
19500101.00 time_ref Reference time for units attribute (yyyymmdd.dd)
3.6000E+02 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
3.6000E+02 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
3.6000E+02 Znudg Nudging/relaxation time scale (days)
for free-surface.
3.6000E+02 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
3.6000E+02 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
1.2000E+02 obcfac Factor between passive and active
open boundary conditions.
T VolCons(1) NLM western edge boundary volume conservation.
T VolCons(2) NLM southern edge boundary volume conservation.
T VolCons(3) NLM eastern edge boundary volume conservation.
T VolCons(4) NLM northern edge boundary volume conservation.
14.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
T LtracerSrc(01) Turning ON point Sources/Sinks on tracer 01: temp
T LtracerSrc(02) Turning ON point Sources/Sinks on tracer 02: salt
F LsshCLM Turning OFF processing of SSH climatology.
F Lm2CLM Turning OFF processing of 2D momentum climatology.
F Lm3CLM Turning OFF processing of 3D momentum climatology.
F LtracerCLM(01) Turning OFF processing of climatology tracer 01: temp
F LtracerCLM(02) Turning OFF processing of climatology tracer 02: salt
F LnudgeM2CLM Turning OFF nudging of 2D momentum climatology.
F LnudgeM3CLM Turning OFF nudging of 3D momentum climatology.
F LnudgeTCLM(01) Turning OFF nudging of climatology tracer 01: temp
F LnudgeTCLM(02) Turning OFF nudging of climatology tracer 02: salt

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 Aout(idFsur) Write out averaged free-surface.
T Aout(idUbar) Write out averaged 2D U-momentum component.
T Aout(idVbar) Write out averaged 2D V-momentum component.
T Aout(idUvel) Write out averaged 3D U-momentum component.
T Aout(idVvel) Write out averaged 3D V-momentum component.
T Aout(idWvel) Write out averaged W-momentum component.
T Aout(idOvel) Write out averaged omega vertical velocity.
T Aout(idTvar) Write out averaged tracer 01: temp
T Aout(idTvar) Write out averaged tracer 02: salt
T Aout(idTsur) Write out averaged surface net heat flux.
T Aout(idTsur) Write out averaged surface net salt flux.
T Aout(idSrad) Write out averaged shortwave radiation flux.


Output/Input Files:

Output Restart File: roms_rst_2000.nc
Output History File: roms_his_2000.nc
Output Averages File: roms_avg_2000.nc
Output Diagnostics File: roms_dia_2000.nc
Output Stations File: roms_sta_2000.nc
Input Grid File: /home/joana/Raylton/run1/Grid/run1_grd1.nc
Input Nonlinear Initial File: /home/joana/Raylton/run1/CC_CI/cmems/initialization_conditions/ROMS_ini_run1.nc
Input Sources/Sinks File: /home/joana/Raylton/run1/CC_CI/rivers/ROMS_river_run1.nc
Input Forcing File 01: /home/joana/Raylton/run1/CC_CI/rivers/ROMS_river_run1.nc
Input Forcing File 02: /home/joana/Raylton/run1/CC_CI/era5/ecna124v_sms_era5_2000.nc
Input Forcing File 03: /home/joana/Raylton/run1/CC_CI/era5/ecna124v_shflux_era5_2000.nc
Input Forcing File 04: /home/joana/Raylton/run1/CC_CI/era5/ecna124v_swflux_era5_2000.nc
Input Forcing File 05: /home/joana/Raylton/run1/CC_CI/era5/ecna124v_swrad_era5_2000.nc
Input Lateral Boundary File 01: /home/joana/Raylton/run1/CC_CI/cmems/boundary_conditions/ROMS_BC_2000.nc
Station positions File: stations_ibis.in
ROMS I/O variables Metadata File: /home/joana/roms-3.9/ROMS/External/varinfo.dat

Tile partition information for Grid 01: 254x110x30 tiling: 1x16

tile Istr Iend Jstr Jend Npts

0 1 254 1 6 45720
1 1 254 7 13 53340
2 1 254 14 20 53340
3 1 254 21 27 53340
4 1 254 28 34 53340
5 1 254 35 41 53340
6 1 254 42 48 53340
7 1 254 49 55 53340
8 1 254 56 62 53340
9 1 254 63 69 53340
10 1 254 70 76 53340
11 1 254 77 83 53340
12 1 254 84 90 53340
13 1 254 91 97 53340
14 1 254 98 104 53340
15 1 254 105 110 45720

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 255.50 0.50 112.50 RHO-points
1 0.50 255.50 -0.50 112.50 RHO-points
2 0.50 255.50 -0.50 112.50 RHO-points
3 0.50 255.50 -0.50 112.50 RHO-points
4 0.50 255.50 -0.50 112.50 RHO-points
5 0.50 255.50 -0.50 112.50 RHO-points
6 0.50 255.50 -0.50 112.50 RHO-points
7 0.50 255.50 -0.50 112.50 RHO-points
8 0.50 255.50 -0.50 112.50 RHO-points
9 0.50 255.50 -0.50 112.50 RHO-points
10 0.50 255.50 -0.50 112.50 RHO-points
11 0.50 255.50 -0.50 112.50 RHO-points
12 0.50 255.50 -0.50 112.50 RHO-points
13 0.50 255.50 -0.50 112.50 RHO-points
14 0.50 255.50 -0.50 112.50 RHO-points
15 0.50 255.50 -0.50 111.50 RHO-points

0 0.00 255.00 0.50 112.50 U-points
1 0.00 255.00 -0.50 112.50 U-points
2 0.00 255.00 -0.50 112.50 U-points
3 0.00 255.00 -0.50 112.50 U-points
4 0.00 255.00 -0.50 112.50 U-points
5 0.00 255.00 -0.50 112.50 U-points
6 0.00 255.00 -0.50 112.50 U-points
7 0.00 255.00 -0.50 112.50 U-points
8 0.00 255.00 -0.50 112.50 U-points
9 0.00 255.00 -0.50 112.50 U-points
10 0.00 255.00 -0.50 112.50 U-points
11 0.00 255.00 -0.50 112.50 U-points
12 0.00 255.00 -0.50 112.50 U-points
13 0.00 255.00 -0.50 112.50 U-points
14 0.00 255.00 -0.50 112.50 U-points
15 0.00 255.00 -0.50 111.50 U-points

0 0.50 255.50 0.00 112.50 V-points
1 0.50 255.50 -0.50 112.50 V-points
2 0.50 255.50 -0.50 112.50 V-points
3 0.50 255.50 -0.50 112.50 V-points
4 0.50 255.50 -0.50 112.50 V-points
5 0.50 255.50 -0.50 112.50 V-points
6 0.50 255.50 -0.50 112.50 V-points
7 0.50 255.50 -0.50 112.50 V-points
8 0.50 255.50 -0.50 112.50 V-points
9 0.50 255.50 -0.50 112.50 V-points
10 0.50 255.50 -0.50 112.50 V-points
11 0.50 255.50 -0.50 112.50 V-points
12 0.50 255.50 -0.50 112.50 V-points
13 0.50 255.50 -0.50 112.50 V-points
14 0.50 255.50 -0.50 112.50 V-points
15 0.50 255.50 -0.50 111.00 V-points


Stations Parameters, Grid: 01
=============================

3 Nstation Number of stations to write out into stations file.
T Sout(idFsur) Write out free-surface.
T Sout(idUbar) Write out 2D U-momentum component.
T Sout(idVbar) Write out 2D V-momentum component.
T Sout(idu2dE) Write out 2D U-eastward at RHO-points.
T Sout(idv2dN) Write out 2D V-northward at RHO-points.
T Sout(idUvel) Write out 3D U-momentum component.
T Sout(idVvel) Write out 3D V-momentum component.
T Sout(idu3dE) Write out 3D U-eastward at RHO-points.
T Sout(idv3dN) Write out 3D V-northward at RHO-points.
T Sout(idWvel) Write out W-momentum component.
T Sout(idOvel) Write out omega vertical velocity.
T Sout(idTvar) Write out tracer 01: temp
T Sout(idTvar) Write out tracer 02: salt
T Sout(idUsms) Write out surface U-momentum stress.
T Sout(idVsms) Write out surface V-momentum stress.
T Sout(idUbms) Write out bottom U-momentum stress.
T Sout(idVbms) Write out bottom V-momentum stress.
T Sout(idTsur) Write out surface net heat flux.
T Sout(idTsur) Write out surface net salt flux.
T Sout(idSrad) Write out shortwave radiation flux.
T Sout(idDano) Write out density anomaly.
T Sout(idVvis) Write out vertical viscosity coefficient.
T Sout(idTdif) Write out vertical T-diffusion coefficient.
T Sout(idSdif) Write out vertical S-diffusion coefficient.
T Sout(idMtke) Write out turbulent kinetic energy.
T Sout(idMtls) Write out turbulent generic length-scale.

Flag and positions for station 0001: 1 -9.3500 43.8333
Flag and positions for station 0002: 1 -9.6500 42.2167
Flag and positions for station 0003: 2 50.0000 225.0000

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 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 Gradient Gradient Gradient Gradient

Activated C-preprocessing Options:

IBIS IBIS
ANA_BSFLUX Analytical kinematic bottom salinity flux
ANA_BTFLUX Analytical kinematic bottom temperature flux
ASSUMED_SHAPE Using assumed-shape arrays
AVERAGES Writing out time-averaged nonlinear model fields
DIAGNOSTICS_TS Computing and writing tracer diagnostic terms
DIAGNOSTICS_UV Computing and writing momentum diagnostic terms
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002)
DOUBLE_PRECISION Double precision arithmetic numerical kernel.
KANTHA_CLAYSON Kantha and Clayson stability function formulation
MASKING Land/Sea masking
MIX_GEO_TS Mixing of tracers along geopotential surfaces
MIX_GEO_UV Mixing of momentum along geopotential surfaces
MY25_MIXING Mellor/Yamada Level-2.5 mixing closure
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
!RST_SINGLE Double precision fields in restart NetCDF file
SALINITY Using salinity
SOLAR_SOURCE Solar Radiation Source Term
SOLVE3D Solving 3D Primitive Equations
STATIONS Writing out station data
TS_MPDATA_LIMIT Further limiter in upwind corrector fluxes for stability
TS_DIF2 Harmonic mixing of tracers
UV_ADV Advection of momentum
UV_COR Coriolis term
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum
UV_C4VADVECTION Fourth-order centered vertical advection of momentum
UV_QDRAG Quadratic bottom stress
UV_VIS2 Harmonic mixing of momentum
VAR_RHO_2D Variable density barotropic mode

Process Information:

Thread # 0 (pid= 1843949) is active.

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

GET_GRID - angle between XI-axis and EAST: angler
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.60268702E+00 Max = -1.53890563E+00)
GET_GRID - bathymetry at RHO-points: h
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.70095704E+03 Max = 5.27800395E+03)
GET_GRID - Coriolis parameter at RHO-points: f
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 8.40525343E-05 Max = 1.04707691E-04)
GET_GRID - reciprocal XI-grid spacing: pm
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 2.14445144E-04 Max = 2.14445144E-04)
GET_GRID - reciprocal ETA-grid spacing: pn
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 2.39808472E-04 Max = 2.81540215E-04)
GET_GRID - longitude of RHO-points: lon_rho
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.21700238E+01 Max = -7.07652660E+00)
GET_GRID - longitude of U-points: lon_u
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.21700238E+01 Max = -7.07652660E+00)
GET_GRID - longitude of V-points: lon_v
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.21470801E+01 Max = -7.09947028E+00)
GET_GRID - longitude of PSI-points: lon_psi
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = -1.21470801E+01 Max = -7.09947028E+00)
GET_GRID - latitude of RHO-points lat_rho
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 3.51930988E+01 Max = 4.58865384E+01)
GET_GRID - latitude of U-points: lat_u
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 3.52140691E+01 Max = 4.58655748E+01)
GET_GRID - latitude of V-points: lat_v
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 3.51930988E+01 Max = 4.58865384E+01)
GET_GRID - latitude of PSI-points lat_psi
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 3.52140691E+01 Max = 4.58655748E+01)
GET_GRID - mask on RHO-points: mask_rho
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on U-points: mask_u
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on V-points: mask_v
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on PSI-points: mask_psi
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - horizontal viscosity sponge factor: visc_factor
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 1.00000000E+00 Max = 3.00000000E+00)
GET_GRID - horizontal diffusivity sponge factor: diff_factor
(Grid = 01, File: /home/joana/Raylton/run1/Grid/run1_grd1.nc)
(Min = 1.00000000E+00 Max = 3.00000000E+00)

Vertical S-coordinate System, Grid 01:

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

30 0.0000000 0.0000000 0.000 0.000 0.000 0.000
29 -0.0333333 -0.0012758 0.519 -0.865 -3.841 -8.321
28 -0.0666667 -0.0051074 5.516 -1.794 -12.129 -30.006
27 -0.1000000 -0.0115076 15.015 -2.788 -24.886 -65.120
26 -0.1333333 -0.0204966 29.051 -3.846 -42.147 -113.770
25 -0.1666667 -0.0321017 47.672 -4.969 -63.960 -176.098
24 -0.2000000 -0.0463550 70.933 -6.159 -90.380 -252.272
23 -0.2333333 -0.0632923 98.898 -7.416 -121.471 -342.479
22 -0.2666667 -0.0829499 131.630 -8.740 -157.294 -446.909
21 -0.3000000 -0.1053620 169.190 -10.134 -197.910 -565.742
20 -0.3333333 -0.1305570 211.626 -11.597 -243.367 -699.124
19 -0.3666667 -0.1585533 258.971 -13.130 -293.699 -847.153
18 -0.4000000 -0.1893544 311.232 -14.734 -348.911 -1009.847
17 -0.4333333 -0.2229439 368.380 -16.407 -408.974 -1187.119
16 -0.4666667 -0.2592789 430.339 -18.149 -473.814 -1378.747
15 -0.5000000 -0.2982845 496.978 -19.957 -543.300 -1584.338
14 -0.5333333 -0.3398462 568.096 -21.829 -617.234 -1803.293
13 -0.5666667 -0.3838043 643.415 -23.762 -695.337 -2034.778
12 -0.6000000 -0.4299468 722.561 -25.749 -777.241 -2277.684
11 -0.6333333 -0.4780043 805.063 -27.783 -862.477 -2530.602
10 -0.6666667 -0.5276456 890.341 -29.858 -950.469 -2791.802
9 -0.7000000 -0.5784757 977.702 -31.962 -1040.528 -3059.216
8 -0.7333333 -0.6300360 1066.343 -34.084 -1131.859 -3330.449
7 -0.7666667 -0.6818081 1155.355 -36.212 -1223.557 -3602.789
6 -0.8000000 -0.7332223 1243.740 -38.331 -1314.633 -3873.258
5 -0.8333333 -0.7836698 1330.430 -40.425 -1404.028 -4138.672
4 -0.8666667 -0.8325209 1414.323 -42.480 -1490.644 -4395.740
3 -0.9000000 -0.8791471 1494.317 -44.479 -1573.389 -4641.175
2 -0.9333333 -0.9229486 1569.361 -46.407 -1651.220 -4871.841
1 -0.9666667 -0.9633838 1638.505 -48.251 -1723.194 -5084.906
0 -1.0000000 -1.0000000 1700.957 -50.000 -1788.523 -5278.004

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

Primary Secondary Accumulated to Current Step

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

ndtfast, nfast = 30 42 nfast/ndtfast = 1.40000

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

1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.18933

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

Minimum X-grid spacing, DXmin = 4.66319722E+00 km Water points = 4.66319722E+00 km
Maximum X-grid spacing, DXmax = 4.66319723E+00 km Water points = 4.66319723E+00 km
Minimum Y-grid spacing, DYmin = 3.55189044E+00 km Water points = 3.55189044E+00 km
Maximum Y-grid spacing, DYmax = 4.16999445E+00 km Water points = 4.16999445E+00 km
Minimum Z-grid spacing, DZmin = -1.31577343E+04 m Water points = 6.68401986E-01 m
Maximum Z-grid spacing, DZmax = 2.28713544E+04 m Water points = 2.72340269E+02 m

Minimum barotropic Courant Number = 5.93191142E-02
Maximum barotropic Courant Number = 7.74959916E-01
Maximum Coriolis Courant Number = 3.14123073E-02


NLM: GET_STATE - Reading state initial conditions, 2000-01-01 12:00:00.00
(Grid 01, t = 18262.5000, File: ROMS_ini_run1.nc, Rec=0001, Index=1)
- free-surface
(Min = -4.32987267E-01 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = -1.12396738E-01 Max = 1.10821879E-01)
- vertically integrated v-momentum component
(Min = -2.19604691E-01 Max = 1.38281519E-01)
- u-momentum component
(Min = -3.27469798E-01 Max = 3.22340740E-01)
- v-momentum component
(Min = -4.73506684E-01 Max = 2.94499485E-01)
- potential temperature
(Min = 0.00000000E+00 Max = 1.85206014E+01)
- salinity
(Min = 0.00000000E+00 Max = 3.68152505E+01)
GET_NGFLD - river runoff XI-positions at RHO-points
(Grid = 01, Min = 6.40000000E+01 Max = 1.04000000E+02)
GET_NGFLD - river runoff ETA-positions at RHO-points
(Grid = 01, Min = 9.60000000E+01 Max = 2.13000000E+02)
GET_NGFLD - river runoff direction
(Grid = 01, Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_NGFLD - river runoff mass transport vertical profile
(Grid = 01, Min = 0.00000000E+00 Max = 2.00000000E-01)
GET_NGFLD - river runoff potential temperature, 2000-01-01 00:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_river_run1.nc)
(Tmin= 18262.0000 Tmax= 73042.0000) t = 18262.0000
(Min = 1.47940000E+01 Max = 1.99500000E+01)
GET_NGFLD - river runoff salinity, 2000-01-01 00:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_river_run1.nc)
(Tmin= 18262.0000 Tmax= 73042.0000) t = 18262.0000
(Min = 1.57000000E+01 Max = 3.55000000E+01)
GET_2DFLD - surface u-momentum stress, 2000-01-01 11:00:00.00
(Grid=01, Rec=12, Index=1, File: ecna124v_sms_era5_2000.nc)
(Tmin= 18262.0000 Tmax= 18627.0000) t = 18262.4583
(Min = -2.62452790E-04 Max = 3.97277247E-05) regrid = T
GET_2DFLD - surface v-momentum stress, 2000-01-01 11:00:00.00
(Grid=01, Rec=12, Index=1, File: ecna124v_sms_era5_2000.nc)
(Tmin= 18262.0000 Tmax= 18627.0000) t = 18262.4583
(Min = -1.18589529E-04 Max = 1.04898793E-04) regrid = T
GET_2DFLD - solar shortwave radiation flux, 2000-01-01 11:30:00.00
(Grid=01, Rec=13, Index=1, File: ecna124v_swrad_era5_2000.nc)
(Tmin= 18261.9792 Tmax= 18626.9792) t = 18262.4792
(Min = 1.73880850E-05 Max = 1.16621777E-04) regrid = T
GET_2DFLD - surface net heat flux, 2000-01-01 11:30:00.00
(Grid=01, Rec=13, Index=1, File: ecna124v_shflux_era5_2000.nc)
(Tmin= 18261.9792 Tmax= 18626.9792) t = 18262.4792
(Min = 7.11119790E-06 Max = 5.67094035E-05) regrid = T
GET_NGFLD - 2D u-momentum western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -6.53617531E-02 Max = 9.55667123E-02)
GET_NGFLD - 2D v-momentum western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -5.20291477E-02 Max = 6.95130751E-02)
GET_NGFLD - 2D u-momentum eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -1.12437271E-01 Max = 7.48302564E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -6.42233565E-02 Max = 2.69341581E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -8.88442695E-02 Max = 1.01377249E-01)
GET_NGFLD - 2D v-momentum southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -1.06768280E-01 Max = 1.10900223E-01)
GET_NGFLD - 2D u-momentum northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -6.46682158E-02 Max = 9.21968892E-02)
GET_NGFLD - 2D v-momentum northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -1.75664142E-01 Max = 9.94500220E-02)
GET_NGFLD - 3D u-momentum western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -1.68131411E-01 Max = 1.81981206E-01)
GET_NGFLD - 3D v-momentum western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -1.47421896E-01 Max = 1.38651356E-01)
GET_NGFLD - 3D u-momentum eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -2.77225375E-01 Max = 3.22363853E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -2.80117363E-01 Max = 8.99794772E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -2.27593422E-01 Max = 2.02705994E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -2.47028366E-01 Max = 2.25783244E-01)
GET_NGFLD - 3D u-momentum northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -2.34356746E-01 Max = 1.43075764E-01)
GET_NGFLD - 3D v-momentum northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = -4.88355607E-01 Max = 2.02198505E-01)
GET_NGFLD - temperature western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 1.98817396E+00 Max = 1.30999126E+01)
GET_NGFLD - salinity western boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 3.48995514E+01 Max = 3.60373383E+01)
GET_NGFLD - temperature eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 2.04223514E+00 Max = 1.81060047E+01)
GET_NGFLD - salinity eastern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 3.49053917E+01 Max = 3.64541130E+01)
GET_NGFLD - temperature southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 1.96197915E+00 Max = 1.77900620E+01)
GET_NGFLD - salinity southern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 3.48929863E+01 Max = 3.63163452E+01)
GET_NGFLD - temperature northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 0.00000000E+00 Max = 1.72163239E+01)
GET_NGFLD - salinity northern boundary condition, 2000-01-01 12:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18262.5000
(Min = 0.00000000E+00 Max = 3.66767006E+01)

Basin information for Grid 01:

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

Initial domain volumes: TotVolume = 1.2742752400E+15 m3
MinCellVol = 1.2173304024E+07 m3
MaxCellVol = 5.0240679525E+09 m3
Max/Min = 4.1271194269E+02


NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000000000145 - 000000104976)

GET_NGFLD - river runoff potential temperature, 2000-01-31 00:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_river_run1.nc)
(Tmin= 18262.0000 Tmax= 73042.0000) t = 18292.0000
(Min = 1.47940000E+01 Max = 1.99500000E+01)
GET_NGFLD - river runoff salinity, 2000-01-31 00:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_river_run1.nc)
(Tmin= 18262.0000 Tmax= 73042.0000) t = 18292.0000
(Min = 1.57000000E+01 Max = 3.55000000E+01)
GET_2DFLD - surface u-momentum stress, 2000-01-01 12:00:00.00
(Grid=01, Rec=13, Index=2, File: ecna124v_sms_era5_2000.nc)
(Tmin= 18262.0000 Tmax= 18627.0000) t = 18262.5000
(Min = -2.65224179E-04 Max = 7.47214743E-05) regrid = T
GET_2DFLD - surface v-momentum stress, 2000-01-01 12:00:00.00
(Grid=01, Rec=13, Index=2, File: ecna124v_sms_era5_2000.nc)
(Tmin= 18262.0000 Tmax= 18627.0000) t = 18262.5000
(Min = -1.25773325E-04 Max = 1.30807716E-04) regrid = T
GET_2DFLD - solar shortwave radiation flux, 2000-01-01 12:30:00.00
(Grid=01, Rec=14, Index=2, File: ecna124v_swrad_era5_2000.nc)
(Tmin= 18261.9792 Tmax= 18626.9792) t = 18262.5208
(Min = 2.46481305E-05 Max = 1.24405759E-04) regrid = T
GET_2DFLD - surface net heat flux, 2000-01-01 12:30:00.00
(Grid=01, Rec=14, Index=2, File: ecna124v_shflux_era5_2000.nc)
(Tmin= 18261.9792 Tmax= 18626.9792) t = 18262.5208
(Min = 1.23045191E-05 Max = 7.02479404E-05) regrid = T
GET_NGFLD - 2D u-momentum western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -6.49111196E-02 Max = 9.59445834E-02)
GET_NGFLD - 2D v-momentum western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -4.95847687E-02 Max = 7.01998323E-02)
GET_NGFLD - 2D u-momentum eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -1.13728225E-01 Max = 6.92293793E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -6.42820746E-02 Max = 2.59706806E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -9.05174986E-02 Max = 9.94924977E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -1.05292454E-01 Max = 1.07233137E-01)
GET_NGFLD - 2D u-momentum northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -6.23283498E-02 Max = 7.24404529E-02)
GET_NGFLD - 2D v-momentum northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -2.05803722E-01 Max = 7.54545256E-02)
GET_NGFLD - 3D u-momentum western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -1.68550178E-01 Max = 1.80611610E-01)
GET_NGFLD - 3D v-momentum western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -1.41561508E-01 Max = 1.36841610E-01)
GET_NGFLD - 3D u-momentum eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -2.47794464E-01 Max = 3.19678664E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -3.25482637E-01 Max = 1.16334975E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -2.30487600E-01 Max = 2.05312297E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -2.35349774E-01 Max = 2.29068324E-01)
GET_NGFLD - 3D u-momentum northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -2.58162826E-01 Max = 1.40312970E-01)
GET_NGFLD - 3D v-momentum northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = -5.30910313E-01 Max = 1.95586026E-01)
GET_NGFLD - temperature western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 1.98840737E+00 Max = 1.30311937E+01)
GET_NGFLD - salinity western boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 3.48996239E+01 Max = 3.60337982E+01)
GET_NGFLD - temperature eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 2.04223323E+00 Max = 1.79368706E+01)
GET_NGFLD - salinity eastern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 3.49055405E+01 Max = 3.64442520E+01)
GET_NGFLD - temperature southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 1.96234500E+00 Max = 1.74535332E+01)
GET_NGFLD - salinity southern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 3.48929443E+01 Max = 3.63216476E+01)
GET_NGFLD - temperature northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 0.00000000E+00 Max = 1.72353230E+01)
GET_NGFLD - salinity northern boundary condition, 2000-01-02 12:00:00.00
(Grid= 01, Rec=2, Index=1, File: ROMS_BC_2000.nc)
(Tmin= 18262.5000 Tmax= 18628.5000) t = 18263.5000
(Min = 0.00000000E+00 Max = 3.66278343E+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

144 2000-01-01 12:00:00.00 2.180850E-03 2.088448E+04 2.088449E+04 1.237148E+15
(229,110,05) 1.362394E-02 3.177991E-02 0.000000E+00 4.845332E-01
DEF_HIS - creating history file, Grid 01: roms_his_2000.nc
DEF_AVG - creating average file, Grid 01: roms_avg_2000.nc
DEF_DIAGS - creating diagnostics file, Grid 01: roms_dia_2000.nc
DEF_STATION - creating stations file, Grid 01: roms_sta_2000.nc
GET_2DFLD - surface u-momentum stress, 2000-01-01 13:00:00.00
(Grid=01, Rec=14, Index=1, File: ecna124v_sms_era5_2000.nc)
(Tmin= 18262.0000 Tmax= 18627.0000) t = 18262.5417
(Min = -2.66432694E-04 Max = 1.44876284E-04) regrid = T
GET_2DFLD - surface v-momentum stress, 2000-01-01 13:00:00.00
(Grid=01, Rec=14, Index=1, File: ecna124v_sms_era5_2000.nc)
(Tmin= 18262.0000 Tmax= 18627.0000) t = 18262.5417
(Min = -1.16046625E-04 Max = 1.35441116E-04) regrid = T
145 2000-01-01 12:05:00.00 NaN NaN NaN NaN
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 -1.000000E+20
Found Error: 01 Line: 314 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 302 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,2) in record = 1

Elapsed wall CPU time for each process (seconds):

Thread # 0 CPU: 1.448
Total: 1.448

Nonlinear model elapsed CPU time profile, Grid: 01

Allocation and array initialization .............. 0.361 (24.9275 %)
Ocean state initialization ....................... 0.015 ( 1.0669 %)
Reading of input data ............................ 0.092 ( 6.3265 %)
Processing of input data ......................... 0.001 ( 0.0658 %)
Processing of output time averaged data .......... 0.132 ( 9.0899 %)
Computation of vertical boundary conditions ...... 0.000 ( 0.0188 %)
Computation of global information integrals ...... 0.008 ( 0.5601 %)
Writing of output data ........................... 0.027 ( 1.8775 %)
Model 2D kernel .................................. 0.243 (16.7837 %)
2D/3D coupling, vertical metrics ................. 0.022 ( 1.5247 %)
Omega vertical velocity .......................... 0.008 ( 0.5759 %)
Equation of state for seawater ................... 0.031 ( 2.1364 %)
My2.5 vertical mixing parameterization ........... 0.054 ( 3.7357 %)
3D equations right-side terms .................... 0.049 ( 3.3714 %)
3D equations predictor step ...................... 0.087 ( 5.9985 %)
Pressure gradient ................................ 0.017 ( 1.1805 %)
Harmonic mixing of tracers, geopotentials ........ 0.033 ( 2.2930 %)
Harmonic stress tensor, geopotentials ............ 0.048 ( 3.3110 %)
Corrector time-step for 3D momentum .............. 0.090 ( 6.1929 %)
Corrector time-step for tracers .................. 0.052 ( 3.5772 %)
Total: 1.370 94.6140 %

Unique kernel(s) regions profiled ................ 1.370 94.6140 %
Residual, non-profiled code ...................... 0.078 5.3860 %


All percentages are with respect to total time = 1.448


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

Dynamic and Automatic memory (MB) usage for Grid 01: 254x110x30 tiling: 1x16

tile Dynamic Automatic USAGE

0 1134.84 14.22 1149.06
1 0.00 14.22 14.22
2 0.00 14.22 14.22
3 0.00 14.22 14.22
4 0.00 14.22 14.22
5 0.00 14.22 14.22
6 0.00 14.22 14.22
7 0.00 14.22 14.22
8 0.00 14.22 14.22
9 0.00 14.22 14.22
10 0.00 14.22 14.22
11 0.00 14.22 14.22
12 0.00 14.22 14.22
13 0.00 14.22 14.22
14 0.00 14.22 14.22
15 0.00 14.22 14.22

TOTAL 1134.84 227.54 1362.38

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

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

Analytical header files used:

ROMS/Functionals/ana_btflux.h

MAIN: Abnormal termination: BLOWUP.
REASON: KEchar = NaN, PEchar = NaN

[wilkin]

You still have negative DZmin, so you have negative bottom depth h somewhere - probably in the land but that can still cause problems. Make h > 0 in the land.

You also have:

GET_NGFLD - river runoff mass transport vertical profile
(Grid = 01, Min = 0.00000000E+00 Max = 2.00000000E-01)
...
GET_NGFLD - river runoff salinity, 2000-01-01 00:00:00.00
(Grid= 01, Rec=1, Index=2, File: ROMS_river_run1.nc)
(Tmin= 18262.0000 Tmax= 73042.0000) t = 18262.0000
(Min = 1.57000000E+01 Max = 3.55000000E+01)

river runoff mass transport vertical profile is variable river_Vshape in your rivers file. It is the relative thickness of the layers into which the river flows. You have a 0, which is OK, but do the numbers sum to 1? If not, the volume flux is incorrectly distributed.

I also see you have rivers with salinity of 35.5. Did you mean to do that?

[rayltonbendo]

Hi Kate
Thank's for your attention!
I'm having the following problem:
Could you help me?


GET_2DFLD - solar shortwave radiation flux, 1999-01-02 17:30:00.00
(Grid=01, Rec=43, Index=1, File: ecna124v_swrad_era5_1999.nc)
(Tmin= 17896.9792 Tmax= 18260.9792) t = 17898.7292
(Min = 1.58338110E-20 Max = 6.11142803E-06) regrid = T
GET_2DFLD - surface net heat flux, 1999-01-02 17:30:00.00
(Grid=01, Rec=43, Index=1, File: ecna124v_shflux_era5_1999.nc)
(Tmin= 17896.9792 Tmax= 18260.9792) t = 17898.7292
(Min = -6.53167318E-05 Max = 3.27726450E-06) regrid = T
343 1999-01-02 16:35:00.00 8.494296E-02 2.040997E+04 2.041005E+04 1.266520E+15
(177,014,28) 2.822723E-02 5.655456E-01 3.367617E+00 3.596474E+01
Found Error: 01 Line: 314 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 302 Source: ROMS/Drivers/nl_ocean.h, ROMS_run

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

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

Elapsed wall CPU time for each process (seconds):

Thread # 0 CPU: 258.232
Total: 258.232

Nonlinear model elapsed CPU time profile, Grid: 01

Allocation and array initialization .............. 0.311 ( 0.1206 %)
Ocean state initialization ....................... 0.015 ( 0.0058 %)
Reading of input data ............................ 0.451 ( 0.1746 %)
Processing of input data ......................... 0.198 ( 0.0768 %)
Processing of output time averaged data .......... 25.063 ( 9.7057 %)
Computation of vertical boundary conditions ...... 0.108 ( 0.0417 %)
Computation of global information integrals ...... 1.428 ( 0.5531 %)
Writing of output data ........................... 1.851 ( 0.7166 %)
Model 2D kernel .................................. 84.687 (32.7949 %)
2D/3D coupling, vertical metrics ................. 2.136 ( 0.8271 %)
Omega vertical velocity .......................... 1.923 ( 0.7448 %)
Equation of state for seawater ................... 3.612 ( 1.3987 %)
My2.5 vertical mixing parameterization ........... 17.629 ( 6.8267 %)
3D equations right-side terms .................... 10.060 ( 3.8958 %)
3D equations predictor step ...................... 31.924 (12.3627 %)
Pressure gradient ................................ 4.528 ( 1.7533 %)
Harmonic mixing of tracers, geopotentials ........ 8.990 ( 3.4813 %)
Harmonic stress tensor, geopotentials ............ 14.581 ( 5.6464 %)
Corrector time-step for 3D momentum .............. 29.176 (11.2984 %)
Corrector time-step for tracers .................. 17.889 ( 6.9275 %)
Total: 256.559 99.3525 %

Unique kernel(s) regions profiled ................ 256.559 99.3525 %
Residual, non-profiled code ...................... 1.672 0.6475 %


All percentages are with respect to total time = 258.232


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

Dynamic and Automatic memory (MB) usage for Grid 01: 254x110x30 tiling: 1x16

tile Dynamic Automatic USAGE

0 1134.84 14.22 1149.06
1 0.00 14.22 14.22
2 0.00 14.22 14.22
3 0.00 14.22 14.22
4 0.00 14.22 14.22
5 0.00 14.22 14.22
6 0.00 14.22 14.22
7 0.00 14.22 14.22
8 0.00 14.22 14.22
9 0.00 14.22 14.22
10 0.00 14.22 14.22
11 0.00 14.22 14.22
12 0.00 14.22 14.22
13 0.00 14.22 14.22
14 0.00 14.22 14.22
15 0.00 14.22 14.22

TOTAL 1134.84 227.54 1362.38

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

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 = 2
number of time records written in STATION file = 343

Analytical header files used:

ROMS/Functionals/ana_btflux.h

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