blowing up error

[javadi]

Dear users,
I have faced blowing up when running a model based on my domain. I have made my all grid, initial, and forcing files. In the upwelling example, I added my input files and I started running the model. The model has faced blowing up after 5-time steps without any changing in upwelling.in(dt=300). I have changed my dt to 150 s, the model ran for 12-time steps,s and the model run for 25 timesteps after decreasing dt to 75 s. Surprisingly, I have considered dt=15 and 5 s but the model have faced this problem after only 3 timesteps. It means that the dt=75 s is the best choice for me in this case. Without my forcing file, I have faced blowing up too.

Please help with your interesting comments.



--------------------------------------------------------------------------------
Model Input Parameters: ROMS/TOMS version 3.9
Tuesday - December 15, 2020 - 11:44:08 AM
--------------------------------------------------------------------------------

Wind-Driven Upwelling/Downwelling over a Periodic Channel

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/trunk
SVN Revision : 1042

Local Root : /home/jav/Roms/trunk
Header Dir : /home/jav/Test
Header file : upwelling.h
Analytical Dir : /home/javi/Test

Resolution, Grid 01: 236x405x32, Parallel Threads: 1, Tiling: 1x1


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

1440 ntimes Number of timesteps for 3-D equations.
10.000 dt Timestep size (s) for 3-D equations.
30 ndtfast Number of timesteps for 2-D equations between
each 3D timestep.
1 ERstr Starting ensemble/perturbation run number.
1 ERend Ending ensemble/perturbation run number.
0 nrrec Number of restart records to read from disk.
T LcycleRST Switch to recycle time-records in restart file.
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).
2 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
7.0000E+00 theta_s S-coordinate surface control parameter.
2.0000E+00 theta_b S-coordinate bottom control parameter.
200.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
0.000 dstart Time-stamp assigned to model initialization (days).
20180101.06 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.
0.0000E+00 obcfac Factor between passive and active
open boundary conditions.
F VolCons(1) NLM western edge boundary volume conservation.
F VolCons(2) NLM southern edge boundary volume conservation.
F VolCons(3) NLM eastern edge boundary volume conservation.
F VolCons(4) NLM northern edge boundary volume conservation.
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: croco_grd.nc
Input Nonlinear Initial File: croco_ini.nc
Input Forcing File 01: gom_cloud_era.nc
Input Forcing File 02: gom_latent_era.nc
Input Forcing File 03: gom_lwrad_era.nc
Input Forcing File 04: gom_Pair_era.nc
Input Forcing File 05: gom_Qair_era.nc
Input Forcing File 06: gom_rain_era.nc
Input Forcing File 07: gom_sensible_era.nc
Input Forcing File 08: gom_shflux_era.nc
Input Forcing File 09: gom_sms_era.nc
Input Forcing File 10: gom_swflux_era.nc
Input Forcing File 11: gom_swrad_era.nc
Input Forcing File 12: gom_Tair_era.nc
Input Forcing File 13: gom_wind_era.nc
ROMS I/O variables Metadata File: varinfo.dat

Tile partition information for Grid 01: 236x405x32 tiling: 1x1

tile Istr Iend Jstr Jend Npts

0 1 236 1 405 3058560

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

tile Xmin Xmax Ymin Ymax grid

0 -1.50 238.50 0.50 405.50 RHO-points

0 -2.00 238.00 0.50 405.50 U-points

0 -1.50 238.50 0.00 405.00 V-points

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 Periodic Closed Periodic Closed

ubar 1 Periodic Closed Periodic Closed

vbar 1 Periodic Closed Periodic Closed

u 1 Periodic Closed Periodic Closed

v 1 Periodic Closed Periodic Closed

temp 1 Periodic Closed Periodic Closed

salt 1 Periodic Closed Periodic 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_SSFLUX Analytical kinematic surface salinity flux
ANA_STFLUX Analytical kinematic surface temperature flux
ANA_VMIX Analytical vertical mixing coefficients
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.
MIX_S_TS Mixing of tracers along constant S-surfaces
MIX_S_UV Mixing of momentum along constant S-surfaces
NONLINEAR Nonlinear Model
!NONLIN_EOS Linear Equation of State for seawater
POWER_LAW Power-law shape time-averaging barotropic filter
PROFILE Time profiling activated
!RST_SINGLE Double precision fields in restart NetCDF file
SALINITY Using salinity
SOLVE3D Solving 3D Primitive Equations
SPLINES_VDIFF Parabolic Spline Reconstruction for Vertical Diffusion
SPLINES_VVISC Parabolic Spline Reconstruction for Vertical Viscosity
TS_DIF2 Harmonic mixing of tracers
UV_ADV Advection of momentum
UV_COR Coriolis term
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum
UV_C4VADVECTION Fourth-order centered vertical advection of momentum
UV_LDRAG Linear bottom stress
UV_VIS2 Harmonic mixing of momentum
VAR_RHO_2D Variable density barotropic mode

Process Information:

Thread # 0 (pid= 13456) is active.

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

GET_GRID - angle between XI-axis and EAST: angler
(Grid = 01, File: croco_grd.nc)
(Min = 2.05175782E-04 Max = 2.59455950E-04)
GET_GRID - bathymetry at RHO-points: h
(Grid = 01, File: croco_grd.nc)
(Min = 7.50000000E+01 Max = 1.21280280E+02)
GET_GRID - Coriolis parameter at RHO-points: f
(Grid = 01, File: croco_grd.nc)
(Min = 8.57239581E-05 Max = 1.08402615E-04)
GET_GRID - reciprocal XI-grid spacing: pm
(Grid = 01, File: croco_grd.nc)
(Min = 2.78060917E-04 Max = 3.36270552E-04)
GET_GRID - reciprocal ETA-grid spacing: pn
(Grid = 01, File: croco_grd.nc)
(Min = 2.78118000E-04 Max = 3.36008921E-04)
GET_GRID - x-location of RHO-points: x-rho
(Grid = 01, File: croco_grd.nc)
(Min = 0.00000000E+00 Max = 8.52331219E+05)
GET_GRID - x-location of U-points: x_u
(Grid = 01, File: croco_grd.nc)
(Min = 1.48689797E+03 Max = 8.50533052E+05)
GET_GRID - x-location of V-points: x_v
(Grid = 01, File: croco_grd.nc)
(Min = 0.00000000E+00 Max = 8.52156273E+05)
GET_GRID - x-location of PSI-points: x_psi
(Grid = 01, File: croco_grd.nc)
(Min = 1.48728387E+03 Max = 8.50358475E+05)
GET_GRID - y-location of RHO-points: y_rho
(Grid = 01, File: croco_grd.nc)
(Min = 0.00000000E+00 Max = 1.33494208E+06)
GET_GRID - y-location of U-points: y_u
(Grid = 01, File: croco_grd.nc)
(Min = 0.00000000E+00 Max = 1.33494208E+06)
GET_GRID - y-location of V-points: y_v
(Grid = 01, File: croco_grd.nc)
(Min = 1.79779806E+03 Max = 1.33345403E+06)
GET_GRID - y-location of PSI-points: y-psi
(Grid = 01, File: croco_grd.nc)
(Min = 1.79779806E+03 Max = 1.33345403E+06)
GET_GRID - longitude of RHO-points: lon_rho
(Grid = 01, File: croco_grd.nc)
(Min = 4.60000000E+01 Max = 5.54800000E+01)
GET_GRID - longitude of U-points: lon_u
(Grid = 01, File: croco_grd.nc)
(Min = 4.60200000E+01 Max = 5.54600000E+01)
GET_GRID - longitude of V-points: lon_v
(Grid = 01, File: croco_grd.nc)
(Min = 4.60000000E+01 Max = 5.54800000E+01)
GET_GRID - longitude of PSI-points: lon_psi
(Grid = 01, File: croco_grd.nc)
(Min = 4.60200000E+01 Max = 5.54600000E+01)
GET_GRID - latitude of RHO-points lat_rho
(Grid = 01, File: croco_grd.nc)
(Min = 3.60000000E+01 Max = 4.80121300E+01)
GET_GRID - latitude of U-points: lat_u
(Grid = 01, File: croco_grd.nc)
(Min = 3.60000000E+01 Max = 4.80121300E+01)
GET_GRID - latitude of V-points: lat_v
(Grid = 01, File: croco_grd.nc)
(Min = 3.60161803E+01 Max = 4.79987436E+01)
GET_GRID - latitude of PSI-points lat_psi
(Grid = 01, File: croco_grd.nc)
(Min = 3.60161803E+01 Max = 4.79987436E+01)

Vertical S-coordinate System, Grid 01:

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

32 0.0000000 -0.0000000 -0.000 -0.000 -0.000 -0.000
31 -0.0312500 -0.0005983 -1.717 -2.387 -2.077 -2.387
30 -0.0625000 -0.0017244 -3.444 -4.798 -4.170 -4.798
29 -0.0937500 -0.0036555 -5.188 -7.245 -6.290 -7.245
28 -0.1250000 -0.0068265 -6.958 -9.750 -8.450 -9.750
27 -0.1562500 -0.0119111 -8.766 -12.342 -10.671 -12.342
26 -0.1875000 -0.0199406 -10.635 -15.069 -12.988 -15.069
25 -0.2187500 -0.0324732 -12.596 -18.002 -15.450 -18.002
24 -0.2500000 -0.0518170 -14.696 -21.247 -18.133 -21.247
23 -0.2812500 -0.0812961 -17.004 -24.956 -21.142 -24.956
22 -0.3125000 -0.1254893 -19.612 -29.338 -24.627 -29.338
21 -0.3437500 -0.1902675 -22.642 -34.663 -28.777 -34.663
20 -0.3750000 -0.2822874 -26.229 -41.236 -33.807 -41.236
19 -0.4062500 -0.4074733 -30.494 -49.326 -39.909 -49.326
18 -0.4375000 -0.5682244 -35.486 -59.045 -47.159 -59.045
17 -0.4687500 -0.7600506 -41.115 -70.187 -55.414 -70.187
16 -0.5000000 -0.9698301 -47.110 -82.150 -64.248 -82.150
15 -0.5312500 -1.1781602 -53.076 -94.047 -73.036 -94.047
14 -0.5625000 -1.3655685 -58.614 -104.986 -81.147 -104.986
13 -0.5937500 -1.5187208 -63.451 -114.357 -88.152 -114.357
12 -0.6250000 -1.6327622 -67.488 -121.938 -93.894 -121.938
11 -0.6562500 -1.7095563 -70.764 -127.813 -98.432 -127.813
10 -0.6875000 -1.7542959 -73.383 -132.220 -101.934 -132.220
9 -0.7187500 -1.7726751 -75.464 -135.421 -104.586 -135.421
8 -0.7500000 -1.7693246 -77.100 -137.627 -106.535 -137.627
7 -0.7812500 -1.7472460 -78.353 -138.976 -107.879 -138.976
6 -0.8125000 -1.7077733 -79.250 -139.528 -108.661 -139.528
5 -0.8437500 -1.6507162 -79.787 -139.275 -108.875 -139.275
4 -0.8750000 -1.5745057 -79.933 -138.145 -108.470 -138.145
3 -0.9062500 -1.4762679 -79.628 -136.007 -107.354 -136.007
2 -0.9375000 -1.3518043 -78.787 -132.668 -105.391 -132.668
1 -0.9687500 -1.1954774 -77.294 -127.870 -102.398 -127.870
0 -1.0000000 -1.0000000 -75.000 -121.280 -98.140 -121.280

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 = 2.97379593E+00 km
Maximum X-grid spacing, DXmax = 3.59633426E+00 km
Minimum Y-grid spacing, DYmin = 2.97611146E+00 km
Maximum Y-grid spacing, DYmax = 3.59559612E+00 km
Minimum Z-grid spacing, DZmin = -6.59006007E+00 m
Maximum Z-grid spacing, DZmax = 1.19634782E+01 m

Minimum barotropic Courant Number = 3.55584979E-03
Maximum barotropic Courant Number = 5.03180240E-03
Maximum Coriolis Courant Number = 1.08402615E-03


NLM: GET_STATE - Reading state initial conditions, 2018-01-01 01:30:00.00
(Grid 01, t = 0.0000, File: croco_ini.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 7.05641875E-02 Max = 2.01933636E+01)
- salinity
(Min = 1.01258829E+01 Max = 1.46311805E+01)
GET_2DFLD - surface u-momentum stress, 2018-01-01 01:30:00.00
(Grid=01, Rec=1, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0000 Tmax= 30.8750) t = 0.0000
(Min = -1.22882390E-04 Max = 2.05808057E-04) regrid = T
GET_2DFLD - surface v-momentum stress, 2018-01-01 01:30:00.00
(Grid=01, Rec=1, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0000 Tmax= 30.8750) t = 0.0000
(Min = -3.57454272E-05 Max = 3.87079847E-04) regrid = T

Basin information for Grid 01:

Maximum grid stiffness ratios: rx0 = 7.632490E-02 (Beckmann and Haidvogel)
rx1 = 1.934840E+04 (Haney)

Initial domain volumes: TotVolume = 7.7739947385E+13 m3
MinCellVol = -6.8815918813E+07 m3
MaxCellVol = 1.2492719892E+08 m3
Max/Min = -1.8153822702E+00


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

GET_2DFLD - surface u-momentum stress, 2018-01-01 04:30:00.00
(Grid=01, Rec=2, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0000 Tmax= 30.8750) t = 0.1250
(Min = -1.75362703E-04 Max = 1.57320342E-04) regrid = T
GET_2DFLD - surface v-momentum stress, 2018-01-01 04:30:00.00
(Grid=01, Rec=2, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0000 Tmax= 30.8750) t = 0.1250
(Min = -3.97578254E-05 Max = 4.33428997E-04) regrid = T

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

0 2018-01-01 01:30:00.00 0.000000E+00 3.691321E+02 3.691321E+02 7.773995E+13
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: roms_his.nc
WRT_HIS - wrote history fields (Index=1,1) in record = 1
DEF_AVG - creating average file, Grid 01: roms_avg.nc
DEF_DIAGS - creating diagnostics file, Grid 01: roms_dia.nc
1 2018-01-01 01:30:10.00 2.874661E-09 3.691321E+02 3.691321E+02 7.773995E+13
(147,236,05) 9.764726E-05 1.403401E-05 -0.000000E+00 4.519250E-02
2 2018-01-01 01:30:20.00 -5.907644E-07 3.691321E+02 3.691321E+02 7.773995E+13
(147,236,06) 3.261724E-02 7.156293E-04 3.810380E-03 1.740675E+01
3 2018-01-01 01:30:30.00 -7.108906E-02 3.691321E+02 3.690610E+02 7.773995E+13
(147,236,06) 1.160866E+01 6.953663E-02 1.320967E+00 6.188637E+03
Found Error: 01 Line: 313 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 300 Source: ROMS/Drivers/nl_ocean.h

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

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

Elapsed wall CPU time for each process (seconds):

Thread # 0 CPU: 61.562
Total: 61.562

Nonlinear model elapsed CPU time profile, Grid: 01

Allocation and array initialization .............. 10.423 (16.9308 %)
Ocean state initialization ....................... 1.948 ( 3.1647 %)
Reading of input data ............................ 0.156 ( 0.2528 %)
Processing of input data ......................... 0.084 ( 0.1357 %)
Processing of output time averaged data .......... 13.829 (22.4642 %)
Computation of vertical boundary conditions ...... 0.004 ( 0.0073 %)
Computation of global information integrals ...... 0.383 ( 0.6225 %)
Writing of output data ........................... 1.327 ( 2.1562 %)
Model 2D kernel .................................. 9.668 (15.7039 %)
2D/3D coupling, vertical metrics ................. 0.449 ( 0.7291 %)
Omega vertical velocity .......................... 0.141 ( 0.2291 %)
Equation of state for seawater ................... 0.254 ( 0.4132 %)
3D equations right-side terms .................... 2.974 ( 4.8307 %)
3D equations predictor step ...................... 7.901 (12.8342 %)
Pressure gradient ................................ 1.611 ( 2.6167 %)
Harmonic mixing of tracers, S-surfaces ........... 1.014 ( 1.6467 %)
Harmonic stress tensor, S-surfaces ............... 0.813 ( 1.3212 %)
Corrector time-step for 3D momentum .............. 2.695 ( 4.3779 %)
Corrector time-step for tracers .................. 3.327 ( 5.4040 %)
Total: 59.002 95.8409 %

Unique kernel(s) regions profiled ................ 59.002 95.8409 %
Residual, non-profiled code ...................... 2.560 4.1591 %


All percentages are with respect to total time = 61.562


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

Dynamic and Automatic memory (MB) usage for Grid 01: 236x405x32 tiling: 1x1

tile Dynamic Automatic USAGE

0 3724.86 184.06 3908.92

TOTAL 3724.86 184.06 3908.92

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

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

Analytical header files used:

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

MAIN: Abnormal termination: BLOWUP.
REASON: MaxSpeed = 6.2E+03
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO
ERROR STOP

Error termination. Backtrace:
#0 0x7efd361ffdf1 in ???
#1 0x7efd36200939 in ???
#2 0x7efd36202047 in ???
#3 0x403994 in ???
#4 0x7efd35e61041 in ???
#5 0x4039ed in ???
#6 0xffffffffffffffff in ???
[root@localhost Test]#

[kate]

This is when you need to start looking at the model's netCDF output as well as just the text. Where is the model blowing up and in what way? Is it a boundary problem or not?

[hezhi]

hi kate
The surface velocity is always 1.4 m/s at boundary。I want to konw whether .h file is important for simulation or not。For example I choose TS_U3HADVECTION and
TS_A4HADVECTION or TS_A4VADVECTION and TS_C4VADVECTION.THere are different for result u and result v？I am newer。


good luck

[kate]

You don't have an open boundary. Where is the velocity always 1.4 m/s? Do you expect it to be less? Is that where the blow-up is happening or not? Actually, it reports: REASON: MaxSpeed = 6.2E+03
Rather a larger velocity. So where is that and why?

I usually use U3H advection and some fourth-order vertical advection. I use the new HSIMT for bio tracers to keep them well-behaved. My typical mode of blowing up is a vertical advection problem, mid-depth, something that can be fixed with a smaller timestep.

Your cppdefs look fine except that I would only use MIX_S_TS over a flat bottom. Try one of the rotated mixing tensor options instead.
If you want, you can try a quadratic bottom drag if you turn on the bottom drag limiter. I forget, is that an option in Hernan's code? It should be.

[javadi]

Dear kate and hezi

I try to run upwelling example and here is my upwelling.h file. I guess that the problem comes from my initial condition. Without initial condition, the model ran without any problem using ANA file. Surprisingly, when I apply my grid file, the minimum z grid and minimum volume are positive. However, after using the initial condition both parameters are negative(?).

Sincerely Javad Babagoli

[kate]

I bet you have some vertically unstable T and S profiles. If they are what you want, you should turn on a vertical mixing scheme like GLS or LMD.

[javadi]

Dear kate
Thank you so much for your response. I have another question in this case. As there is not GLS or LMD flags in my upwelling.h, I should copy this flag in cppdefs and add this flag in my upwelling.h like this:
#define LMD_MIXING
Or should I another action.
Please help me in this case

[kate]

Yes, exactly. However, you know how you have the #ifdef GLS_MIXING... section? Here's what I have for LMD:

Code: Select all

# define LMD_MIXING
# ifdef LMD_MIXING
#  define LMD_RIMIX
#  define LMD_CONVEC
#  define LMD_SKPP
#  define LI_FOX_KEMPER
#  undef LMD_BKPP
#  define LMD_NONLOCAL
#  define LMD_SHAPIRO
#  define LMD_DDMIX
#  define LIMIT_VDIFF
# endif

You might want to turn off DDMIX (double diffusion), depending. Oh, and I believe the trunk code won't have LI_FOX_KEMPER.

[javadi]

Dear Kate

Thank you so much for your cooperation. I have applied your suggestion and I have changed my upwelling.h. However, I have faced this error:

CHECKDEFS - only one vertical closure scheme is allowed.
Found Error: 05 Line: 1255 Source: ROMS/Utility/inp_par.F
Found Error: 05 Line: 112 Source: ROMS/Drivers/nl_ocean.h

upwelling.h file:

/*
** svn $Id: upwelling.h 1001 2020-01-10 22:41:16Z arango $
*******************************************************************************
** Copyright (c) 2002-2020 The ROMS/TOMS Group **
** Licensed under a MIT/X style license **
** See License_ROMS.txt **
*******************************************************************************
**
** Options for Upwelling Test.
**
** Application flag: UPWELLING
** Input script: roms_upwelling.in
*/

#define UV_ADV
#define UV_COR
#define UV_LDRAG
#define UV_VIS2
#undef MIX_GEO_UV
#define MIX_S_UV
#define SPLINES_VDIFF
#define SPLINES_VVISC
#define DJ_GRADPS
#define TS_DIF2
#undef TS_DIF4
#undef MIX_GEO_TS
#define MIX_S_TS

# define LMD_MIXING
# ifdef LMD_MIXING
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_SKPP
# define LI_FOX_KEMPER
# undef LMD_BKPP
# define LMD_NONLOCAL
# define LMD_SHAPIRO
# define LMD_DDMIX
# define LIMIT_VDIFF
# endif

#define SALINITY
#define SOLVE3D
#define AVERAGES
#define DIAGNOSTICS_TS
#define DIAGNOSTICS_UV

#undef ANA_GRID
#undef ANA_INITIAL
#undef ANA_SMFLUX
#define ANA_STFLUX
#define ANA_SSFLUX
#define ANA_BTFLUX
#define ANA_BSFLUX

#if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# define RI_SPLINES
#else
# define ANA_VMIX
#endif

#if defined BIO_FENNEL || defined ECOSIM || \
defined NPZD_POWELL || defined NEMURO
# define ANA_BIOLOGY
# define ANA_SPFLUX
# define ANA_BPFLUX
# define ANA_SRFLUX
#endif

#if defined NEMURO
# define HOLLING_GRAZING
# undef IVLEV_EXPLICIT
#endif

#ifdef BIO_FENNEL
# define CARBON
# define DENITRIFICATION
# define BIO_SEDIMENT
# define DIAGNOSTICS_BIO
#endif

#ifdef PERFECT_RESTART
# undef AVERAGES
# undef DIAGNOSTICS_BIO
# undef DIAGNOSTICS_TS
# undef DIAGNOSTICS_UV
# define OUT_DOUBLE
#endif

[wilkin]

You have defined LMD_MIXING, but then you also have this block of definitions ...

Code: Select all

#if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# define RI_SPLINES
#else
# define ANA_VMIX
#endif

The else clause activates #define ANA_VMIX, so you have a conflict, both LMD_MIXING and ANA_VMIX are defined, hence the error from CHECKDEFS.

[javadi]

Dear wikin

Thank you so much for your help.I have reformed my .h based on new conditions and I have turned ANA_VMIX off but the model is unstable too.

[wilkin]

Is this supposed to a real ocean, or something idealized? You have some odd looking ranges of initial temperature and salinity:

- potential temperature
(Min = 7.05641875E-02 Max = 2.01933636E+01)
- salinity
(Min = 1.01258829E+01 Max = 1.46311805E+01)

0.7 C to 20 C for temperature.
10 to 14 for salinity.

Are these backwards?

[javadi]

Dear wikin

Yes. It is a real area. The values for temperature and salinity are correct based on my previous experiences.

[javadi]

It can lead to instability as we do not consider ice in the model(?)The water converts to ice when the temperature is approximately zero. The observation supports this; however, in the model, I do not consider ice. If this leads to instability, how can I add ice to the model?

[kate]

Do you want to have ice in the model? There are numerous options to do so, including my branch, COAWST, Metroms, and Scott Durski's branch. Metroms (on github) has a version coupled to the community CICE model, the most complete/mature of the ice models.

On the other hand, before I had ice in the model, it was happily running with temperatures down to -20 C. No problem, except the equation of state isn't valid there.

[wilkin]

Then I'm out of suggestions.

But I add a few comments ...

The code base you are using must be extremely old because in roms_ini.nc you have a variable scrum_time, in addition to ocean_time. The SCRUM model went away a decade or more ago, so you should be careful you have not inherited some outdated code in preparing your inputs.

Also, the negative dz is a worry. You need to track that down.

And finally, the combination of Vtransform = 2 with Vstretching = 1 is an odd one. The Vtransform = 2 coordinate was never a part of SCRUM, so I worry that your code to make the initial conditions is not interpreting the vertical coordinate correctly.

Since you are starting out, why not make small incremental changes on the UPWELLING test case? So, keep Vtransform etc. the same as in the analytical case, and only change your initial conditions at first. Maybe keep the ANA surface forcing too until you know your initial conditions are not the problem.