grid-tide blowup

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
AMIRHOSSEIN
Posts: 43
Joined: Tue Feb 11, 2020 6:03 pm
Location: IUST

grid-tide blowup

#1 Unread post by AMIRHOSSEIN »

Hi,

The first question, which is a bit basic, what are the priorities for checking the output in .log file?

The second question is that the grid file has 30 vertical layers, but by changing the boundary and time, I could not prevent the model error(in 30 vertical layers), but considering 16 layers (the grid file is the same 30), the model was executed correctly for a year and I do not know where the problem comes from?!

The third question is that How can the following errors occur?
MaxSpeed & MaxDensity

Thanks a lot
Attachments
roms-a-year.log
(55.18 MiB) Downloaded 257 times
roms.log
(3.66 MiB) Downloaded 142 times

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

Re: grid-tide blowup

#2 Unread post by kate »

So the same configuration runs with 16 vertical levels and blows up with 30? This is not surprising. When you add more vertical levels, the vertical advection requires a shorter timestep. Your run ran long enough that it seems fine except for the timestep limitation. Try a smaller DT. Your vertical CFL number got too big here:

Code: Select all

     20689 0001-01-08 04:24:30.00  9.228638E-02  1.185957E+04  1.185966E+04  1.811791E+14
                     (292,059,30)  2.852839E-02  1.282760E-02  1.111797E+01  2.098128E+01
(1.111797E+01) Remember, this needs to stay less than one.

AMIRHOSSEIN
Posts: 43
Joined: Tue Feb 11, 2020 6:03 pm
Location: IUST

Re: grid-tide blowup

#3 Unread post by AMIRHOSSEIN »

Hi kate,
thanks a lot for your response :D
I'll do it and announce the result.

AMIRHOSSEIN
Posts: 43
Joined: Tue Feb 11, 2020 6:03 pm
Location: IUST

Re: grid-tide blowup

#4 Unread post by AMIRHOSSEIN »

I face this error when I use a shorter time step with 30 vertical levels:

--------------------------------------------------------------------------------
Model Input Parameters: ROMS/TOMS version 3.9
Saturday - May 15, 2021 - 1:30:09 AM
--------------------------------------------------------------------------------

Wind-Driven Upwelling/Downwelling over a Periodic Channel

Operating system : Linux
CPU/hardware : x86_64
Compiler system : ifort
Compiler command : /home/siadat/LIBRARIES/mvapich/bin/mpif90
Compiler flags : -fp-model precise -heap-arrays -ip -O3 -traceback -check uninit
OCN Communicator : 1140850688, PET size = 36

Input Script : roms_upwelling.in

GIT Root URL : https://www.myroms.org/git/src
GIT Revision : 766ca53c47d4e26e0329b305ac9cd9caaab44da6
SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision : 1061

Local Root : /home/siadat/AMIR/src
Header Dir : /home/siadat/AMIR/grid-tide-sazman/only-tide-run
Header file : upwelling.h
Analytical Dir : /home/siadat/AMIR/grid-tide-sazman/only-tide-run

Resolution, Grid 01: 296x197x30, Parallel Nodes: 36, Tiling: 9x4


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

2102400 ntimes Number of timesteps for 3-D equations.
15.000 dt Timestep size (s) for 3-D equations.
3 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.
T ldefout Switch to create a new output NetCDF file(s).
72 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.
72 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
5.0000E+00 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
F LuvSponge Turning OFF sponge on horizontal momentum.
F LtracerSponge(01) Turning OFF sponge on tracer 01: temp
F LtracerSponge(02) Turning OFF sponge on tracer 02: salt
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.
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.0000E-01 Dcrit Minimum depth for wetting and drying (m).
2 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
6.0000E+00 theta_s S-coordinate surface control parameter.
0.0000E+00 theta_b S-coordinate bottom control parameter.
10.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
0.500 dstart Time-stamp assigned to model initialization (days).
*********** tide_start Reference time origin for tidal forcing (days).
20121028.00 time_ref Reference time for units attribute (yyyymmdd.dd)
0.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
0.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
0.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
0.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
0.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
2.0000E+00 obcfac Factor between passive and active
open boundary conditions.
F VolCons(1) NLM western edge boundary volume conservation.
F VolCons(2) NLM southern edge boundary volume conservation.
F VolCons(3) NLM eastern edge boundary volume conservation.
F VolCons(4) NLM northern edge boundary volume conservation.
14.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1027.000 R0 Background density (kg/m3) used in linear Equation
of State.
1.7000E-04 Tcoef Thermal expansion coefficient (1/Celsius).
0.0000E+00 Scoef Saline contraction coefficient (1/PSU).
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
F LtracerSrc(01) Turning OFF point Sources/Sinks on tracer 01: temp
F LtracerSrc(02) Turning OFF point Sources/Sinks on tracer 02: salt
F LsshCLM Turning OFF processing of SSH climatology.
F Lm2CLM Turning OFF processing of 2D momentum climatology.
F Lm3CLM Turning OFF processing of 3D momentum climatology.
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 Dout(M2rate) Write out 2D momentum acceleration.
T Dout(M2pgrd) Write out 2D momentum pressure gradient.
T Dout(M2fcor) Write out 2D momentum Coriolis force.
T Dout(M2hadv) Write out 2D momentum horizontal advection.
T Dout(M2xadv) Write out 2D momentum horizontal X-advection.
T Dout(M2yadv) Write out 2D momentum horizontal Y-advection.
T Dout(M2hvis) Write out 2D momentum horizontal viscosity.
T Dout(M2xvis) Write out 2D momentum horizontal X-viscosity.
T Dout(M2yvis) Write out 2D momentum horizontal Y-viscosity.
T Dout(M2sstr) Write out 2D momentum surface stress.
T Dout(M2bstr) Write out 2D momentum bottom stress.

T Dout(M3rate) Write out 3D momentum acceleration.
T Dout(M3pgrd) Write out 3D momentum pressure gradient.
T Dout(M3fcor) Write out 3D momentum Coriolis force.
T Dout(M3hadv) Write out 3D momentum horizontal advection.
T Dout(M3xadv) Write out 3D momentum horizontal X-advection.
T Dout(M3yadv) Write out 3D momentum horizontal Y-advection.
T Dout(M3vadv) Write out 3D momentum vertical advection.
T Dout(M3hvis) Write out 3D momentum horizontal viscosity.
T Dout(M3xvis) Write out 3D momentum horizontal X-viscosity.
T Dout(M3yvis) Write out 3D momentum horizontal Y-viscosity.
T Dout(M3vvis) Write out 3D momentum vertical viscosity.

T Dout(iTrate) Write out rate of change of tracer 01: temp
T Dout(iTrate) Write out rate of change of tracer 02: salt
T Dout(iThadv) Write out horizontal advection, tracer 01: temp
T Dout(iThadv) Write out horizontal advection, tracer 02: salt
T Dout(iTxadv) Write out horizontal X-advection, tracer 01: temp
T Dout(iTxadv) Write out horizontal X-advection, tracer 02: salt
T Dout(iTyadv) Write out horizontal Y-advection, tracer 01: temp
T Dout(iTyadv) Write out horizontal Y-advection, tracer 02: salt
T Dout(iTvadv) Write out vertical advection, tracer 01: temp
T Dout(iTvadv) Write out vertical advection, tracer 02: salt
T Dout(iThdif) Write out horizontal diffusion, tracer 01: temp
T Dout(iThdif) Write out horizontal diffusion, tracer 02: salt
T Dout(iTxdif) Write out horizontal X-diffusion, tracer 01: temp
T Dout(iTxdif) Write out horizontal X-diffusion, tracer 02: salt
T Dout(iTydif) Write out horizontal Y-diffusion , tracer 01: temp
T Dout(iTydif) Write out horizontal Y-diffusion , tracer 02: salt
T Dout(iTvdif) Write out vertical diffusion, tracer 01: temp
T Dout(iTvdif) Write out vertical diffusion, tracer 02: salt

Output/Input Files:

Output Restart File: roms_rst.nc
Output History File: roms_his.nc
Output Averages File: roms_avg.nc
Output Diagnostics File: roms_dia.nc
Input Grid File: Persian_roms_grid.nc
Tidal Forcing File: tide_forc_Persian.nc
ROMS I/O variables Metadata File: /home/siadat/AMIR/grid-tide-sazman/only-tide-run/varinfo.dat

Tile partition information for Grid 01: 296x197x30 tiling: 9x4

tile Istr Iend Jstr Jend Npts

0 1 33 1 49 48510
1 34 66 1 49 48510
2 67 99 1 49 48510
3 100 132 1 49 48510
4 133 165 1 49 48510
5 166 198 1 49 48510
6 199 231 1 49 48510
7 232 264 1 49 48510
8 265 296 1 49 47040
9 1 33 50 99 49500
10 34 66 50 99 49500
11 67 99 50 99 49500
12 100 132 50 99 49500
13 133 165 50 99 49500
14 166 198 50 99 49500
15 199 231 50 99 49500
16 232 264 50 99 49500
17 265 296 50 99 48000
18 1 33 100 149 49500
19 34 66 100 149 49500
20 67 99 100 149 49500
21 100 132 100 149 49500
22 133 165 100 149 49500
23 166 198 100 149 49500
24 199 231 100 149 49500
25 232 264 100 149 49500
26 265 296 100 149 48000
27 1 33 150 197 47520
28 34 66 150 197 47520
29 67 99 150 197 47520
30 100 132 150 197 47520
31 133 165 150 197 47520
32 166 198 150 197 47520
33 199 231 150 197 47520
34 232 264 150 197 47520
35 265 296 150 197 46080

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 33.50 0.50 49.50 RHO-points
1 33.50 66.50 0.50 49.50 RHO-points
2 66.50 99.50 0.50 49.50 RHO-points
3 99.50 132.50 0.50 49.50 RHO-points
4 132.50 165.50 0.50 49.50 RHO-points
5 165.50 198.50 0.50 49.50 RHO-points
6 198.50 231.50 0.50 49.50 RHO-points
7 231.50 264.50 0.50 49.50 RHO-points
8 264.50 296.50 0.50 49.50 RHO-points
9 0.50 33.50 49.50 99.50 RHO-points
10 33.50 66.50 49.50 99.50 RHO-points
11 66.50 99.50 49.50 99.50 RHO-points
12 99.50 132.50 49.50 99.50 RHO-points
13 132.50 165.50 49.50 99.50 RHO-points
14 165.50 198.50 49.50 99.50 RHO-points
15 198.50 231.50 49.50 99.50 RHO-points
16 231.50 264.50 49.50 99.50 RHO-points
17 264.50 296.50 49.50 99.50 RHO-points
18 0.50 33.50 99.50 149.50 RHO-points
19 33.50 66.50 99.50 149.50 RHO-points
20 66.50 99.50 99.50 149.50 RHO-points
21 99.50 132.50 99.50 149.50 RHO-points
22 132.50 165.50 99.50 149.50 RHO-points
23 165.50 198.50 99.50 149.50 RHO-points
24 198.50 231.50 99.50 149.50 RHO-points
25 231.50 264.50 99.50 149.50 RHO-points
26 264.50 296.50 99.50 149.50 RHO-points
27 0.50 33.50 149.50 197.50 RHO-points
28 33.50 66.50 149.50 197.50 RHO-points
29 66.50 99.50 149.50 197.50 RHO-points
30 99.50 132.50 149.50 197.50 RHO-points
31 132.50 165.50 149.50 197.50 RHO-points
32 165.50 198.50 149.50 197.50 RHO-points
33 198.50 231.50 149.50 197.50 RHO-points
34 231.50 264.50 149.50 197.50 RHO-points
35 264.50 296.50 149.50 197.50 RHO-points

0 1.00 33.50 0.50 49.50 U-points
1 33.50 66.50 0.50 49.50 U-points
2 66.50 99.50 0.50 49.50 U-points
3 99.50 132.50 0.50 49.50 U-points
4 132.50 165.50 0.50 49.50 U-points
5 165.50 198.50 0.50 49.50 U-points
6 198.50 231.50 0.50 49.50 U-points
7 231.50 264.50 0.50 49.50 U-points
8 264.50 296.00 0.50 49.50 U-points
9 1.00 33.50 49.50 99.50 U-points
10 33.50 66.50 49.50 99.50 U-points
11 66.50 99.50 49.50 99.50 U-points
12 99.50 132.50 49.50 99.50 U-points
13 132.50 165.50 49.50 99.50 U-points
14 165.50 198.50 49.50 99.50 U-points
15 198.50 231.50 49.50 99.50 U-points
16 231.50 264.50 49.50 99.50 U-points
17 264.50 296.00 49.50 99.50 U-points
18 1.00 33.50 99.50 149.50 U-points
19 33.50 66.50 99.50 149.50 U-points
20 66.50 99.50 99.50 149.50 U-points
21 99.50 132.50 99.50 149.50 U-points
22 132.50 165.50 99.50 149.50 U-points
23 165.50 198.50 99.50 149.50 U-points
24 198.50 231.50 99.50 149.50 U-points
25 231.50 264.50 99.50 149.50 U-points
26 264.50 296.00 99.50 149.50 U-points
27 1.00 33.50 149.50 197.50 U-points
28 33.50 66.50 149.50 197.50 U-points
29 66.50 99.50 149.50 197.50 U-points
30 99.50 132.50 149.50 197.50 U-points
31 132.50 165.50 149.50 197.50 U-points
32 165.50 198.50 149.50 197.50 U-points
33 198.50 231.50 149.50 197.50 U-points
34 231.50 264.50 149.50 197.50 U-points
35 264.50 296.00 149.50 197.50 U-points

0 0.50 33.50 1.00 49.50 V-points
1 33.50 66.50 1.00 49.50 V-points
2 66.50 99.50 1.00 49.50 V-points
3 99.50 132.50 1.00 49.50 V-points
4 132.50 165.50 1.00 49.50 V-points
5 165.50 198.50 1.00 49.50 V-points
6 198.50 231.50 1.00 49.50 V-points
7 231.50 264.50 1.00 49.50 V-points
8 264.50 296.50 1.00 49.50 V-points
9 0.50 33.50 49.50 99.50 V-points
10 33.50 66.50 49.50 99.50 V-points
11 66.50 99.50 49.50 99.50 V-points
12 99.50 132.50 49.50 99.50 V-points
13 132.50 165.50 49.50 99.50 V-points
14 165.50 198.50 49.50 99.50 V-points
15 198.50 231.50 49.50 99.50 V-points
16 231.50 264.50 49.50 99.50 V-points
17 264.50 296.50 49.50 99.50 V-points
18 0.50 33.50 99.50 149.50 V-points
19 33.50 66.50 99.50 149.50 V-points
20 66.50 99.50 99.50 149.50 V-points
21 99.50 132.50 99.50 149.50 V-points
22 132.50 165.50 99.50 149.50 V-points
23 165.50 198.50 99.50 149.50 V-points
24 198.50 231.50 99.50 149.50 V-points
25 231.50 264.50 99.50 149.50 V-points
26 264.50 296.50 99.50 149.50 V-points
27 0.50 33.50 149.50 197.00 V-points
28 33.50 66.50 149.50 197.00 V-points
29 66.50 99.50 149.50 197.00 V-points
30 99.50 132.50 149.50 197.00 V-points
31 132.50 165.50 149.50 197.00 V-points
32 165.50 198.50 149.50 197.00 V-points
33 198.50 231.50 149.50 197.00 V-points
34 231.50 264.50 149.50 197.00 V-points
35 264.50 296.50 149.50 197.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 90
HaloSizeJ(1) = 124
TileSide(1) = 56
TileSize(1) = 2184


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

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

temp 1 Upstream3 Centered4

salt 1 Upstream3 Centered4

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


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

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

zeta 1 Closed Radiation Radiation Closed

ubar 1 Closed Radiation Radiation Closed

vbar 1 Closed Radiation Radiation Closed

u 1 Closed Radiation Radiation Closed

v 1 Closed Radiation Radiation Closed

temp 1 Closed Radiation Radiation Closed

salt 1 Closed Radiation Radiation Closed

Activated C-preprocessing Options:

UPWELLING Wind-Driven Upwelling/Downwelling over a Periodic Channel
ANA_BSFLUX Analytical kinematic bottom salinity flux
ANA_BTFLUX Analytical kinematic bottom temperature flux
ANA_FSOBC Analytical free-surface boundary conditions
ANA_INITIAL Analytical initial conditions
ANA_M2OBC Analytical 2D momentum boundary conditions
ANA_SMFLUX Analytical kinematic surface momentum flux
ANA_SSFLUX Analytical kinematic surface salinity flux
ANA_STFLUX Analytical kinematic surface temperature flux
ASSUMED_SHAPE Using assumed-shape arrays
AVERAGES Writing out time-averaged nonlinear model fields
!BOUNDARY_ALLGATHER Using mpi_allreduce in mp_boundary routine
!COLLECT_ALL... Using mpi_isend/mpi_recv in mp_collect routine
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.
LIMIT_BSTRESS Limit bottom stress to maintain bottom velocity direction
MASKING Land/Sea masking
MIX_S_TS Mixing of tracers along constant S-surfaces
MIX_S_UV Mixing of momentum along constant S-surfaces
MPI MPI distributed-memory configuration
NONLINEAR Nonlinear Model
!NONLIN_EOS Linear Equation of State for seawater
POWER_LAW Power-law shape time-averaging barotropic filter
PROFILE Time profiling activated
RAMP_TIDES Ramping tidal forcing for one day
REDUCE_ALLREDUCE Using mpi_allreduce in mp_reduce routine
!RST_SINGLE Double precision fields in restart NetCDF file
SALINITY Using salinity
SOLVE3D Solving 3D Primitive Equations
SSH_TIDES Add tidal elevation to SSH climatology
TS_DIF2 Harmonic mixing of tracers
UV_ADV Advection of momentum
UV_COR Coriolis term
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum
UV_C4VADVECTION Fourth-order centered vertical advection of momentum
UV_LDRAG Linear bottom stress
UV_TIDES Add tidal currents to 2D momentum climatologies
UV_VIS2 Harmonic mixing of momentum
VAR_RHO_2D Variable density barotropic mode
WET_DRY Wetting and drying activated

Process Information:

Node # 0 (pid= 15974) is active.
Node # 1 (pid= 15975) is active.
Node # 2 (pid= 15976) is active.
Node # 3 (pid= 15977) is active.
Node # 4 (pid= 15978) is active.
Node # 5 (pid= 15979) is active.
Node # 6 (pid= 15980) is active.
Node # 7 (pid= 15981) is active.
Node # 10 (pid= 15984) is active.
Node # 11 (pid= 15985) is active.
Node # 12 (pid= 15986) is active.
Node # 13 (pid= 15987) is active.
Node # 14 (pid= 15988) is active.
Node # 15 (pid= 15989) is active.
Node # 16 (pid= 15990) is active.
Node # 17 (pid= 15991) is active.
Node # 18 (pid= 15992) is active.
Node # 19 (pid= 15993) is active.
Node # 20 (pid= 15994) is active.
Node # 21 (pid= 15995) is active.
Node # 22 (pid= 15996) is active.
Node # 23 (pid= 15997) is active.
Node # 24 (pid= 15998) is active.
Node # 25 (pid= 15999) is active.
Node # 27 (pid= 16001) is active.
Node # 28 (pid= 16002) is active.
Node # 29 (pid= 16003) is active.
Node # 30 (pid= 16004) is active.
Node # 31 (pid= 16005) is active.
Node # 32 (pid= 16006) is active.
Node # 33 (pid= 16007) is active.
Node # 34 (pid= 16008) is active.
Node # 35 (pid= 16009) is active.
Node # 26 (pid= 16000) is active.
Node # 8 (pid= 15982) is active.
Node # 9 (pid= 15983) is active.

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

GET_GRID - bathymetry at RHO-points: h
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 1.86978871E+00 Max = 3.47617100E+03)
GET_GRID - Coriolis parameter at RHO-points: f
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 5.69842278E-05 Max = 7.40193949E-05)
GET_GRID - reciprocal XI-grid spacing: pm
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 2.18316466E-04 Max = 2.33234175E-04)
GET_GRID - reciprocal ETA-grid spacing: pn
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 2.37580994E-04 Max = 2.37580994E-04)
GET_GRID - x-location of RHO-points: x-rho
(Grid = 01, File: Persian_roms_grid.nc)
(Min = -7.23213034E+05 Max = 7.23213034E+05)
GET_GRID - y-location of RHO-points: y_rho
(Grid = 01, File: Persian_roms_grid.nc)
(Min = -2.74066636E+06 Max = -1.84531373E+06)
GET_GRID - x-location of PSI-points: x_psi
(Grid = 01, File: Persian_roms_grid.nc)
(Min = -7.20591582E+05 Max = 7.20591582E+05)
GET_GRID - y-location of PSI-points: y-psi
(Grid = 01, File: Persian_roms_grid.nc)
(Min = -2.73842625E+06 Max = -1.84763986E+06)
GET_GRID - x-location of U-points: x_u
(Grid = 01, File: Persian_roms_grid.nc)
(Min = -7.20783957E+05 Max = 7.20783957E+05)
GET_GRID - y-location of U-points: y_u
(Grid = 01, File: Persian_roms_grid.nc)
(Min = -2.74066636E+06 Max = -1.84550147E+06)
GET_GRID - x-location of V-points: x_v
(Grid = 01, File: Persian_roms_grid.nc)
(Min = -7.23020011E+05 Max = 7.23020011E+05)
GET_GRID - y-location of V-points: y_v
(Grid = 01, File: Persian_roms_grid.nc)
(Min = -2.73842625E+06 Max = -1.84745207E+06)
GET_GRID - latitude of RHO-points lat_rho
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 2.30000000E+01 Max = 3.05000000E+01)
GET_GRID - longitude of RHO-points: lon_rho
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 4.77000000E+01 Max = 6.10000000E+01)
GET_GRID - latitude of PSI-points lat_psi
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 2.30189394E+01 Max = 3.04810606E+01)
GET_GRID - longitude of PSI-points: lon_psi
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 4.77223906E+01 Max = 6.09776094E+01)
GET_GRID - latitude of U-points: lat_u
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 2.30000000E+01 Max = 3.05000000E+01)
GET_GRID - longitude of U-points: lon_u
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 4.77223906E+01 Max = 6.09776094E+01)
GET_GRID - latitude of V-points: lat_v
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 2.30189394E+01 Max = 3.04810606E+01)
GET_GRID - longitude of V-points: lon_v
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 4.77000000E+01 Max = 6.10000000E+01)
GET_GRID - mask on RHO-points: mask_rho
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on U-points: mask_u
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on V-points: mask_v
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on PSI-points: mask_psi
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 0.00000000E+00 Max = 2.00000000E+00)
GET_GRID - angle between XI-axis and EAST: angler
(Grid = 01, File: Persian_roms_grid.nc)
(Min = 0.00000000E+00 Max = 0.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.0009981 -0.053 -0.172 -2.057 -3.792
28 -0.0666667 -0.0020363 -0.106 -0.344 -4.184 -7.723
27 -0.1000000 -0.0031562 -0.158 -0.516 -6.452 -11.937
26 -0.1333333 -0.0044028 -0.211 -0.689 -8.939 -16.591
25 -0.1666667 -0.0058261 -0.264 -0.862 -11.731 -21.856
24 -0.2000000 -0.0074832 -0.317 -1.037 -14.928 -27.933
23 -0.2333333 -0.0094406 -0.370 -1.214 -18.644 -35.050
22 -0.2666667 -0.0117770 -0.424 -1.392 -23.015 -43.480
21 -0.3000000 -0.0145859 -0.477 -1.573 -28.203 -53.549
20 -0.3333333 -0.0179803 -0.530 -1.757 -34.404 -65.647
19 -0.3666667 -0.0220963 -0.584 -1.944 -41.852 -80.246
18 -0.4000000 -0.0270990 -0.638 -2.135 -50.833 -97.919
17 -0.4333333 -0.0331894 -0.692 -2.333 -61.696 -119.362
16 -0.4666667 -0.0406117 -0.747 -2.536 -74.861 -145.422
15 -0.5000000 -0.0496640 -0.802 -2.748 -90.844 -177.131
14 -0.5333333 -0.0607094 -0.858 -2.970 -110.274 -215.749
13 -0.5666667 -0.0741913 -0.914 -3.204 -133.917 -262.812
12 -0.6000000 -0.0906508 -0.972 -3.453 -162.708 -320.196
11 -0.6333333 -0.1107484 -1.030 -3.720 -197.790 -390.191
10 -0.6666667 -0.1352907 -1.090 -4.010 -240.557 -475.592
9 -0.7000000 -0.1652627 -1.151 -4.326 -292.712 -579.814
8 -0.7333333 -0.2018673 -1.215 -4.676 -356.336 -707.025
7 -0.7666667 -0.2465736 -1.280 -5.066 -433.968 -862.318
6 -0.8000000 -0.3011757 -1.349 -5.506 -528.710 -1051.912
5 -0.8333333 -0.3678650 -1.421 -6.006 -644.353 -1283.403
4 -0.8666667 -0.4493181 -1.498 -6.580 -785.523 -1566.068
3 -0.9000000 -0.5488038 -1.579 -7.244 -957.873 -1911.238
2 -0.9333333 -0.6703150 -1.668 -8.018 -1168.307 -2332.752
1 -0.9666667 -0.8187283 -1.764 -8.927 -1425.256 -2847.515
0 -1.0000000 -1.0000000 -1.870 -10.000 -1739.020 -3476.171

Time Splitting Weights for Grid 01: ndtfast = 3 nfast = 4
==================================

Primary Secondary Accumulated to Current Step

1 0.0558333333428837 0.3333333333332940 0.0558333333428837 0.3333333333332940
2 0.2149999999995108 0.3147222222190016 0.2708333333423945 0.6480555555522955
3 0.4024999999719729 0.2430555555525065 0.6733333333143674 0.8911111111048020
4 0.3266666666856327 0.1088888888951980 1.0000000000000000 1.0000000000000000

ndtfast, nfast = 3 4 nfast/ndtfast = 1.33333

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

1.000000000000 1.085000000007 0.542500000003 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.00000

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

Minimum X-grid spacing, DXmin = 4.28753633E+00 km Water points = 4.29917804E+00 km
Maximum X-grid spacing, DXmax = 4.58050654E+00 km Water points = 4.58050654E+00 km
Minimum Y-grid spacing, DYmin = 4.20909091E+00 km Water points = 4.20909091E+00 km
Maximum Y-grid spacing, DYmax = 4.20909091E+00 km Water points = 4.20909091E+00 km
Minimum Z-grid spacing, DZmin = 5.28023281E-02 m Water points = 5.58882656E-02 m
Maximum Z-grid spacing, DZmax = 6.28656359E+02 m Water points = 6.28656359E+02 m

Minimum barotropic Courant Number = 7.24874895E-03
Maximum barotropic Courant Number = 2.98941489E-01
Maximum Coriolis Courant Number = 1.10155619E-03

ANA_INITIAL - vertically integrated u-momentum component: ubar
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
ANA_INITIAL - vertically integrated v-momentum component: vbar
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
ANA_INITIAL - free-surface: zeta
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
ANA_INITIAL - u-momentum component: u
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
ANA_INITIAL - v-momentum component: v
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
ANA_INITIAL - potential temperature: temp
(Grid = 01, Min = 1.40000000E+01 Max = 2.19957770E+01)
ANA_INITIAL - salinity: salt
(Grid = 01, Min = 3.50000000E+01 Max = 3.50000000E+01)
GET_NGFLD - tidal period
(Grid = 01, Min = 2.19720215E+04 Max = 2.38071489E+06)
GET_2DFLD - tidal elevation amplitude
(Grid = 01, Min = -1.23526398E-14 Max = 8.75923031E-01)
GET_2DFLD - tidal elevation phase angle
(Grid = 01, Min = -5.49838079E+00 Max = 8.47757976E+00)
GET_2DFLD - tidal current inclination angle
(Grid = 01, Min = 0.00000000E+00 Max = 6.28317385E+00)
GET_2DFLD - tidal current phase angle
(Grid = 01, Min = 0.00000000E+00 Max = 6.28316977E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Grid = 01, Min = 0.00000000E+00 Max = 6.48205227E-01)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Grid = 01, Min = -1.94831549E-01 Max = 2.54299358E-01)

Basin information for Grid 01:

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

Initial domain volumes: TotVolume = 1.7663437614E+14 m3
MinCellVol = 1.0197991434E+06 m3
MaxCellVol = 1.1942855421E+10 m3
Max/Min = 1.1710987893E+04


7360 2012-10-29 18:40:00.00 3.382645E-01 1.202837E+04 1.202871E+04 1.781154E+14
(231,058,30) 5.546413E-03 3.282268E-03 4.716395E+01 1.102684E+01
7361 2012-10-29 18:40:15.00 3.248345E-01 1.202846E+04 1.202879E+04 1.781171E+14
(231,058,30) 4.944411E-03 2.938754E-03 5.116505E+01 1.101317E+01
7362 2012-10-29 18:40:30.00 3.311236E-01 1.202853E+04 1.202886E+04 1.781190E+14
(226,061,30) 4.219131E-03 1.991485E-03 4.882558E+01 1.098131E+01
7363 2012-10-29 18:40:45.00 3.514411E-01 1.202864E+04 1.202899E+04 1.781216E+14
(226,061,30) 5.287721E-03 1.984694E-03 4.620421E+01 1.093203E+01
7364 2012-10-29 18:41:00.00 3.661830E-01 1.202879E+04 1.202916E+04 1.781246E+14
(222,030,30) 8.704506E-04 1.984459E-03 4.140152E+01 1.086639E+01
7365 2012-10-29 18:41:15.00 3.613346E-01 1.202898E+04 1.202934E+04 1.781277E+14
(217,031,30) 1.188826E-04 3.934099E-03 3.474698E+01 1.078582E+01
7366 2012-10-29 18:41:30.00 3.422162E-01 1.202917E+04 1.202951E+04 1.781308E+14
(290,055,30) 6.375751E-03 2.577218E-03 4.430922E+01 1.069207E+01
7367 2012-10-29 18:41:45.00 3.278931E-01 1.202934E+04 1.202967E+04 1.781341E+14
(290,055,30) 5.785064E-03 1.180376E-03 4.433389E+01 1.058723E+01
7368 2012-10-29 18:42:00.00 3.316144E-01 1.202947E+04 1.202980E+04 1.781374E+14
(219,061,30) 3.756705E-03 2.511664E-03 4.667482E+01 1.047369E+01
7369 2012-10-29 18:42:15.00 3.483822E-01 1.202954E+04 1.202988E+04 1.781404E+14
(219,061,30) 4.099640E-03 1.219892E-03 4.800742E+01 1.035414E+01
7370 2012-10-29 18:42:30.00 3.615934E-01 1.202957E+04 1.202993E+04 1.781428E+14
(226,060,30) 4.365490E-03 1.548314E-04 3.778717E+01 1.023151E+01
7371 2012-10-29 18:42:45.00 3.597649E-01 1.202964E+04 1.203000E+04 1.781451E+14
(208,042,30) 4.137802E-04 1.567138E-03 3.598323E+01 1.010889E+01
Found Error: 01 Line: 315 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 303 Source: ROMS/Drivers/nl_ocean.h, ROMS_run

Blowing-up: Saving latest model state into RESTART file
REASON: MaxDensity = 9.2E+02

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

Elapsed wall CPU time for each process (seconds):

Node # 0 CPU: 2011.695
Node # 1 CPU: 2104.411
Node # 2 CPU: 2085.632
Node # 3 CPU: 2105.727
Node # 4 CPU: 2105.241
Node # 5 CPU: 2044.155
Node # 6 CPU: 2103.442
Node # 7 CPU: 2098.814
Node # 8 CPU: 2106.134
Node # 9 CPU: 2105.983
Node # 10 CPU: 2104.696
Node # 11 CPU: 2104.578
Node # 12 CPU: 2105.703
Node # 13 CPU: 2089.597
Node # 14 CPU: 2106.216
Node # 15 CPU: 2104.822
Node # 16 CPU: 2105.390
Node # 17 CPU: 2105.465
Node # 18 CPU: 2081.495
Node # 19 CPU: 2105.477
Node # 20 CPU: 2087.576
Node # 21 CPU: 2105.457
Node # 22 CPU: 2105.259
Node # 23 CPU: 2086.149
Node # 24 CPU: 2071.530
Node # 25 CPU: 2096.574
Node # 26 CPU: 2106.401
Node # 27 CPU: 2106.110
Node # 28 CPU: 2105.914
Node # 29 CPU: 2104.901
Node # 30 CPU: 2105.715
Node # 31 CPU: 2104.445
Node # 32 CPU: 2105.697
Node # 33 CPU: 2105.275
Node # 34 CPU: 2105.748
Node # 35 CPU: 2105.149
Total: 75492.572
Average: 2097.016
Minimum: 2011.695
Maximum: 2106.401

Nonlinear model elapsed CPU time profile, Grid: 01

Allocation and array initialization .............. 5.881 ( 0.0078 %)
Ocean state initialization ....................... 0.743 ( 0.0010 %)
Reading of input data ............................ 14.296 ( 0.0189 %)
Processing of input data ......................... 281.244 ( 0.3725 %)
Processing of output time averaged data .......... 10469.466 (13.8682 %)
Computation of vertical boundary conditions ...... 57.113 ( 0.0757 %)
Computation of global information integrals ...... 770.025 ( 1.0200 %)
Writing of output data ........................... 17555.355 (23.2544 %)
Model 2D kernel .................................. 7408.343 ( 9.8133 %)
Tidal forcing .................................... 1492.521 ( 1.9770 %)
2D/3D coupling, vertical metrics ................. 1654.458 ( 2.1916 %)
Omega vertical velocity .......................... 984.658 ( 1.3043 %)
Equation of state for seawater ................... 831.253 ( 1.1011 %)
3D equations right-side terms .................... 3756.312 ( 4.9757 %)
3D equations predictor step ...................... 9434.468 (12.4972 %)
Pressure gradient ................................ 1416.116 ( 1.8758 %)
Harmonic mixing of tracers, S-surfaces ........... 1105.157 ( 1.4639 %)
Harmonic stress tensor, S-surfaces ............... 1010.786 ( 1.3389 %)
Corrector time-step for 3D momentum .............. 9252.722 (12.2565 %)
Corrector time-step for tracers .................. 7012.643 ( 9.2892 %)
Total: 74513.559 98.7032 %

Unique kernel(s) regions profiled ................ 74513.559 98.7032 %
Residual, non-profiled code ...................... 979.012 1.2968 %


All percentages are with respect to total time = 75492.572


MPI communications profile, Grid: 01

Message Passage: 2D halo exchanges ............... 3358.061 ( 4.4482 %)
Message Passage: 3D halo exchanges ............... 2212.833 ( 2.9312 %)
Message Passage: 4D halo exchanges ............... 1446.070 ( 1.9155 %)
Message Passage: data broadcast .................. 4277.199 ( 5.6657 %)
Message Passage: data reduction .................. 131.596 ( 0.1743 %)
Message Passage: data gathering .................. 5572.557 ( 7.3816 %)
Message Passage: data scattering.................. 8.030 ( 0.0106 %)
Message Passage: point data gathering ............ 0.002 ( 0.0000 %)
Message Passage: synchronization barrier ......... 0.544 ( 0.0007 %)
Total: 17006.890 22.5279 %

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

Dynamic and Automatic memory (MB) usage for Grid 01: 296x197x30 tiling: 9x4

tile Dynamic Automatic USAGE MPI-Buffers

0 69.04 29.41 98.45 19.50
1 70.93 29.41 100.34 19.50
2 70.93 29.41 100.34 19.50
3 70.93 29.41 100.34 19.50
4 70.93 29.41 100.34 19.50
5 70.93 29.41 100.34 19.50
6 70.93 29.41 100.34 19.50
7 70.93 29.41 100.34 19.50
8 69.04 29.41 98.45 19.50
9 71.66 29.41 101.07 19.50
10 73.62 29.41 103.03 19.50
11 73.62 29.41 103.03 19.50
12 73.62 29.41 103.03 19.50
13 73.62 29.41 103.03 19.50
14 73.62 29.41 103.03 19.50
15 73.62 29.41 103.03 19.50
16 73.62 29.41 103.03 19.50
17 71.66 29.41 101.07 19.50
18 71.66 29.41 101.07 19.50
19 73.62 29.41 103.03 19.50
20 73.62 29.41 103.03 19.50
21 73.62 29.41 103.03 19.50
22 73.62 29.41 103.03 19.50
23 73.62 29.41 103.03 19.50
24 73.62 29.41 103.03 19.50
25 73.62 29.41 103.03 19.50
26 71.66 29.41 101.07 19.50
27 67.73 29.41 97.15 19.50
28 69.59 29.41 99.00 19.50
29 69.59 29.41 99.00 19.50
30 69.59 29.41 99.00 19.50
31 69.59 29.41 99.00 19.50
32 69.59 29.41 99.00 19.50
33 69.59 29.41 99.00 19.50
34 69.59 29.41 99.00 19.50
35 67.73 29.41 97.15 19.50

TOTAL 2574.46 1058.90 3633.35 701.96

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

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

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_fsobc.h
ROMS/Functionals/ana_initial.h
ROMS/Functionals/ana_m2obc.h
ROMS/Functionals/ana_smflux.h
ROMS/Functionals/ana_stflux.h

MAIN: Abnormal termination: BLOWUP.
REASON: MaxDensity = 9.2E+02

Post Reply