Very early blowing up problem

[lcbernardo]

Hi everyone,

I am quite new to ROMS and would like to know how to use it properly. So far, I have successfully run some of the test cases, and also tried running simple coastal applications for some of my study sites. For these cases, I did not really need to have 3D boundary and initial conditions, and could manage with some simple tidal forcing.

Lately however, I am trying to perform some regional simulations, but I am having a difficult time getting them to run, and my early attempts have all led to the model blowing up after a few time-steps. I am trying to run a regional simulation for areas south of Japan, with which I hope to eventually use to try the nesting features.

By modifying some of the available ROMS matlab scripts, I generated some test initial and boundary forcing files which I derived from a non-ROMS coarse-resolution model developed by a colleague. I generated initial and boundary conditions for temperature, salt, u, v, ubar, vbar, and zeta, but since the model I am running is regional in scale, I have not added tidal forcing, which was also the practice followed by my colleague. However, the runs I try keep blowing up after only a few time steps, and I am really confused regarding what the problem is. I don't think that the inputs I have provided are too unrealistic to lead to blowing up right away, but more likely it is due to my inexperience. I would like to ask for help regarding how to properly debug the run.

The following are the options from my cppdefs.h file:

#define UV_VIS2
#define MIX_GEO_UV
#define MIX_GEO_TS
#undef ANA_GRID
#define MASKING
#define UV_ADV
#define UV_COR
#define TS_DIF2
#define TS_U3HADVECTION
#define TS_SADVECTION

#define NONLIN_EOS
#define SALINITY
#define DJ_GRADPS
#define SOLVE3D
#define SPLINES
#define WET_DRY

/*** Option for Boundary condition ***/

#define RADIATION_2D

#define SOLAR_SOURCE

#define BULK_FLUXES
#ifdef BULK_FLUXES
#define LONGWAVE
#define EMINUSP
#else
#define ANA_SMFLUX
#define ANA_STFLUX
#endif

#define UV_LOGDRAG
#define SSW_LOGINT
#define LIMIT_BSTRESS

#ifdef SOLVE3D
#define GLS_MIXING
#endif

#if defined GLS_MIXING
#define KANTHA_CLAYSON
#undef CANUTO_A
#define N2S2_HORAVG
#endif

#define ANA_BPFLUX
#define ANA_BTFLUX
#define ANA_BSFLUX
#define ANA_SPFLUX



When I try to run the model, I then get the following output, which includes the information on my run configuration:



Model Input Parameters: ROMS/TOMS version 3.7
Wednesday - January 7, 2015 - 9:33:37 PM
-----------------------------------------------------------------------------

New Regional Case

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/apps/mvapich2/1.5.1/gcc/bin/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none

Input Script :

SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision : exported

Local Root : /home0/usr8/13D51333/ROMS/trunk
Header Dir : /home/usr8/13D51333/ROMS/Projects/lor_regional_test
Header file : regional_test.h
Analytical Dir: /home/usr8/13D51333/ROMS/Projects/lor_regional_test

Resolution, Grid 01: 0308x0278x030, Parallel Nodes: 8, Tiling: 004x002


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

172800 ntimes Number of timesteps for 3-D equations.
0.500 dt Timestep size (s) for 3-D equations.
80 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.
600 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).
360 nHIS Number of timesteps between the writing fields
into history 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
1.0000E-01 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
F LuvSponge Turning OFF sponge on horizontal momentum.
F LtracerSponge(01) Turning OFF sponge on tracer 01: temp
F LtracerSponge(02) Turning OFF sponge on tracer 02: salt
5.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
5.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
5.0000E-06 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.000 gls_p GLS stability exponent.
1.500 gls_m GLS turbulent kinetic energy exponent.
-1.000 gls_n GLS turbulent length scale exponent.
7.6000E-06 gls_Kmin GLS minimum value of turbulent kinetic energy.
1.0000E-12 gls_Pmin GLS minimum value of dissipation.
5.4770E-01 gls_cmu0 GLS stability coefficient.
1.4400E+00 gls_c1 GLS shear production coefficient.
1.9200E+00 gls_c2 GLS dissipation coefficient.
-4.0000E-01 gls_c3m GLS stable buoyancy production coefficient.
1.0000E+00 gls_c3p GLS unstable buoyancy production coefficient.
1.0000E+00 gls_sigk GLS constant Schmidt number for TKE.
1.3000E+00 gls_sigp GLS constant Schmidt number for PSI.
1400.000 charnok_alpha Charnok factor for Zos calculation.
0.500 zos_hsig_alpha Factor for Zos calculation using Hsig(Awave).
0.250 sz_alpha Factor for Wave dissipation surface tke flux .
100.000 crgban_cw Factor for Craig/Banner surface tke flux.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
2.5000E-02 rdrg2 Quadratic bottom drag coefficient.
3.0000E-02 Zob Bottom roughness (m).
5.0000E-01 Zos Surface roughness (m).
1.0000E+01 blk_ZQ Height (m) of surface air humidity measurement.
1.0000E+01 blk_ZT Height (m) of surface air temperature measurement.
1.0000E+01 blk_ZW Height (m) of surface winds measurement.
1.0000E-01 Dcrit Minimum depth for wetting and drying (m).
4 lmd_Jwt Jerlov water type.
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
7.0000E+00 theta_s S-coordinate surface control parameter.
1.0000E-01 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).
20140610.00 time_ref Reference time for units attribute (yyyymmdd.dd)
1.8000E+02 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
1.8000E+02 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
1.8000E+02 Znudg Nudging/relaxation time scale (days)
for free-surface.
1.8000E+02 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
1.8000E+02 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
1.2000E+02 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.
27.000 T0 Background potential temperature (C) constant.
34.000 S0 Background salinity (PSU) constant.
-1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
F LtracerSrc(01) Turning OFF point Sources/Sinks on tracer 01: temp
F LtracerSrc(02) Turning OFF point Sources/Sinks on tracer 02: salt
F LsshCLM Turning OFF processing of SSH climatology.
F Lm2CLM Turning OFF processing of 2D momentum climatology.
F Lm3CLM Turning OFF processing of 3D momentum climatology.
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

Output/Input Files:

Output Restart File: ocean_rst.nc
Output History File: ocean_his.nc
Input Grid File: Data/regional_grd.nc
Input Nonlinear Initial File: Data/regional_init_drlu_corr1.nc
Input Forcing File 01: Data/regional_test_frc.nc
Input Boundary File: Data/regional_bry_drlu_corr1.nc

Tile partition information for Grid 01: 0308x0278x0030 tiling: 004x002

tile Istr Iend Jstr Jend Npts

0 1 77 1 139 321090
1 78 154 1 139 321090
2 155 231 1 139 321090
3 232 308 1 139 321090
4 1 77 140 278 321090
5 78 154 140 278 321090
6 155 231 140 278 321090
7 232 308 140 278 321090

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 77.50 0.50 139.50 RHO-points
1 77.50 154.50 0.50 139.50 RHO-points
2 154.50 231.50 0.50 139.50 RHO-points
3 231.50 308.50 0.50 139.50 RHO-points
4 0.50 77.50 139.50 278.50 RHO-points
5 77.50 154.50 139.50 278.50 RHO-points
6 154.50 231.50 139.50 278.50 RHO-points
7 231.50 308.50 139.50 278.50 RHO-points

0 1.00 77.50 0.50 139.50 U-points
1 77.50 154.50 0.50 139.50 U-points
2 154.50 231.50 0.50 139.50 U-points
3 231.50 308.00 0.50 139.50 U-points
4 1.00 77.50 139.50 278.50 U-points
5 77.50 154.50 139.50 278.50 U-points
6 154.50 231.50 139.50 278.50 U-points
7 231.50 308.00 139.50 278.50 U-points

0 0.50 77.50 1.00 139.50 V-points
1 77.50 154.50 1.00 139.50 V-points
2 154.50 231.50 1.00 139.50 V-points
3 231.50 308.50 1.00 139.50 V-points
4 0.50 77.50 139.50 278.00 V-points
5 77.50 154.50 139.50 278.00 V-points
6 154.50 231.50 139.50 278.00 V-points
7 231.50 308.50 139.50 278.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 178
HaloSizeJ(1) = 302
TileSide(1) = 145
TileSize(1) = 12035


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

vbar 1 Flather Flather Flather Flather

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:

REGIONAL_TEST New Regional Case
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
BULK_FLUXES Surface bulk fluxes parameterization.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EMINUSP Compute Salt Flux using E-P.
GLS_MIXING Generic Length-Scale turbulence closure.
LIMIT_BSTRESS Limit bottom stress to maintain bottom velocity direction.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
LONGWAVE Compute net longwave radiation internally.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_GEO_UV Mixing of momentum along geopotential surfaces.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
RADIATION_2D Use tangential phase speed in radiation conditions.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
SPLINES Conservative parabolic spline reconstruction.
TS_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
UV_ADV Advection of momentum.
UV_COR Coriolis term.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_C4VADVECTION Fourth-order centered vertical advection of momentum.
UV_LOGDRAG Logarithmic bottom stress.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
WET_DRY Wetting and drying activated.

Process Information:

Node # 0 (pid= 10050) is active.
Node # 1 (pid= 10051) is active.
Node # 2 (pid= 10052) is active.
Node # 3 (pid= 10053) is active.
Node # 4 (pid= 10054) is active.
Node # 5 (pid= 10055) is active.
Node # 6 (pid= 10056) is active.
Node # 7 (pid= 10057) is active.

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

30 0.0000000 0.0000000 -0.000 0.000 0.000 0.000
29 -0.0333333 -0.0000525 6.661 -3.339 -6.444 -6.810
28 -0.0666667 -0.0002129 13.312 -6.688 -13.220 -14.315
27 -0.1000000 -0.0004899 19.951 -10.049 -20.356 -22.571
26 -0.1333333 -0.0008987 26.577 -13.423 -27.897 -31.677
25 -0.1666667 -0.0014616 33.187 -16.813 -35.913 -41.775
24 -0.2000000 -0.0022095 39.779 -20.221 -44.498 -53.065
23 -0.2333333 -0.0031830 46.348 -23.652 -53.779 -65.809
22 -0.2666667 -0.0044355 52.890 -27.110 -63.919 -80.349
21 -0.3000000 -0.0060353 59.396 -30.604 -75.128 -97.128
20 -0.3333333 -0.0080696 65.860 -34.140 -87.676 -116.706
19 -0.3666667 -0.0106496 72.268 -37.732 -101.905 -139.800
18 -0.4000000 -0.0139158 78.608 -41.392 -118.247 -167.314
17 -0.4333333 -0.0180461 84.862 -45.138 -137.250 -200.396
16 -0.4666667 -0.0232653 91.007 -48.993 -159.607 -240.493
15 -0.5000000 -0.0298569 97.014 -52.986 -186.191 -289.432
14 -0.5333333 -0.0381784 102.849 -57.151 -218.103 -349.516
13 -0.5666667 -0.0486802 108.465 -61.535 -256.730 -423.648
12 -0.6000000 -0.0619292 113.807 -66.193 -303.818 -515.477
11 -0.6333333 -0.0786377 118.803 -71.197 -361.561 -629.597
10 -0.6666667 -0.0997005 123.363 -76.637 -432.716 -771.769
9 -0.7000000 -0.1262391 127.376 -82.624 -520.735 -949.219
8 -0.7333333 -0.1596562 130.701 -89.299 -629.940 -1170.985
7 -0.7666667 -0.2017025 133.163 -96.837 -765.722 -1448.346
6 -0.8000000 -0.2545552 134.544 -105.456 -934.787 -1795.330
5 -0.8333333 -0.3209111 134.576 -115.424 -1145.441 -2229.309
4 -0.8666667 -0.4040926 132.924 -127.076 -1407.916 -2771.690
3 -0.9000000 -0.5081650 129.183 -140.817 -1734.735 -3448.664
2 -0.9333333 -0.6380580 122.861 -157.139 -2141.079 -4291.991
1 -0.9666667 -0.7996795 113.365 -176.635 -2645.144 -5339.734
0 -1.0000000 -1.0000000 100.000 -200.000 -3268.400 -6636.800

Time Splitting Weights for Grid 01: ndtfast = 80 nfast = 109
==================================

Primary Secondary Accumulated to Current Step

1-0.0002288942706210 0.0125000000000000-0.0002288942706210 0.0125000000000000
2-0.0004416899831918 0.0125028611783828-0.0006705842538128 0.0250028611783828
3-0.0006383871597875 0.0125083823031727-0.0013089714136003 0.0375112434815554
4-0.0008189858862562 0.0125163621426700-0.0021279572998565 0.0500276056242254
5-0.0009834863858029 0.0125265994662482-0.0031114436856593 0.0625542050904736
6-0.0011318891220071 0.0125388930460707-0.0042433328076664 0.0750930981365444
7-0.0012641949312745 0.0125530416600958-0.0055075277389409 0.0876461397966402
8-0.0013804051847213 0.0125688440967368-0.0068879329236623 0.1002149838933770
9-0.0014805219794934 0.0125860991615458-0.0083684549031557 0.1128010830549228
10-0.0015645483595181 0.0126046056862894-0.0099330032626738 0.1254056887412122
11-0.0016324885656902 0.0126241625407834-0.0115654918283639 0.1380298512819956
12-0.0016843483154913 0.0126445686478546-0.0132498401438553 0.1506744199298502
13-0.0017201351120431 0.0126656230017982-0.0149699752558984 0.1633400429316484
14-0.0017398585825937 0.0126871246906987-0.0167098338384921 0.1760271676223471
15-0.0017435308464381 0.0127088729229812-0.0184533646849303 0.1887360405453282
16-0.0017311669122724 0.0127306670585616-0.0201845315972026 0.2014667076038899
17-0.0017027851049807 0.0127523066449650-0.0218873167021834 0.2142190142488549
18-0.0016584075218572 0.0127735914587773-0.0235457242240405 0.2269926057076322
19-0.0015980605182600 0.0127943215528005-0.0251437847423005 0.2397869272604327
20-0.0015217752227004 0.0128142973092788-0.0266655599650010 0.2526012245697115
21-0.0014295880813646 0.0128333194995625-0.0280951480463656 0.2654345440692740
22-0.0013215414320692 0.0128511893505796-0.0294166894784347 0.2782857334198535
23-0.0011976841076507 0.0128677086184804-0.0306143735860854 0.2911534420383339
24-0.0010580720687888 0.0128826796698261-0.0316724456548742 0.3040361217081600
25-0.0009027690662622 0.0128959055706859-0.0325752147211364 0.3169320272788460
26-0.0007318473326393 0.0129071901840142-0.0333070620537757 0.3298392174628602
27-0.0005453883034021 0.0129163382756722-0.0338524503571778 0.3427555557385323
28-0.0003434833675031 0.0129231556294647-0.0341959337246809 0.3556787113679971
29-0.0001262346473571 0.0129274491715585-0.0343221683720380 0.3686061605395556
30 0.0001062441917348 0.0129290271046505-0.0342159241803032 0.3815351876442061
31 0.0003538271027266 0.0129276990522538-0.0338620970775767 0.3944628866964598
32 0.0006163747885352 0.0129232762134697-0.0332457222890414 0.4073861629099296
33 0.0008937338043156 0.0129155715286130-0.0323519884847259 0.4203017344385426
34 0.0011857356303006 0.0129043998560591-0.0311662528544253 0.4332061342946016
35 0.0014921957152090 0.0128895781606803-0.0296740571392163 0.4460957124552820
36 0.0018129124902188 0.0128709257142402-0.0278611446489975 0.4589666381695222
37 0.0021476663535069 0.0128482643081125-0.0257134782954906 0.4718149024776346
38 0.0024962186253557 0.0128214184786936-0.0232172596701349 0.4846363209563282
39 0.0028583104738258 0.0127902157458767-0.0203589491963091 0.4974265367022049
40 0.0032336618109948 0.0127544868649539-0.0171252873853142 0.5101810235671588
41 0.0036219701597632 0.0127140660923164-0.0135033172255511 0.5228950896594752
42 0.0040229094912255 0.0126687914653194-0.0094804077343255 0.5355638811247946
43 0.0044361290326093 0.0126185050966791-0.0050442787017162 0.5481823862214736
44 0.0048612520457791 0.0125630534837715-0.0001830266559371 0.5607454397052450
45 0.0052978745763075 0.0125022878331992 0.0051148479203703 0.5732477275384442
46 0.0057455641731126 0.0124360644009954 0.0108604120934830 0.5856837919394396
47 0.0062038585786618 0.0123642448488315 0.0170642706721447 0.5980480367882711
48 0.0066722643897413 0.0122866966165982 0.0237365350618861 0.6103347334048693
49 0.0071502556887934 0.0122032933117264 0.0308867907506794 0.6225380267165957
50 0.0076372726458184 0.0121139151156165 0.0385240633964979 0.6346519418322123
51 0.0081327200908446 0.0120184492075438 0.0466567834873425 0.6466703910397561
52 0.0086359660569633 0.0119167902064082 0.0552927495443058 0.6585871812461643
53 0.0091463402939309 0.0118088406306962 0.0644390898382367 0.6703960218768605
54 0.0096631327523374 0.0116945113770221 0.0741022225905740 0.6820905332538826
55 0.0101855920383408 0.0115737222176178 0.0842878146289148 0.6936642554715005
56 0.0107129238389685 0.0114464023171386 0.0950007384678833 0.7051106577886390
57 0.0112442893179843 0.0113124907691515 0.1062450277858676 0.7164231485577904
58 0.0117788034823227 0.0111719371526767 0.1180238312681904 0.7275950857104670
59 0.0123155335190887 0.0110247021091476 0.1303393647872791 0.7386197878196147
60 0.0128534971031246 0.0108707579401590 0.1431928618904036 0.7494905457597737
61 0.0133916606751427 0.0107100892263700 0.1565845225655464 0.7602006349861437
62 0.0139289376904250 0.0105426934679307 0.1705134602559714 0.7707433284540743
63 0.0144641868380886 0.0103685817468004 0.1849776470940600 0.7811119102008747
64 0.0149962102309176 0.0101877794113243 0.1999738573249776 0.7912996896121989
65 0.0155237515657621 0.0100003267834378 0.2154976088907397 0.8013000163956366
66 0.0160454942545025 0.0098062798888658 0.2315431031452422 0.8111062962845024
67 0.0165600595255809 0.0096057112106845 0.2481031626708231 0.8207120074951868
68 0.0170660044960986 0.0093987104666147 0.2651691671669217 0.8301107179618015
69 0.0175618202144801 0.0091853854104135 0.2827309873814018 0.8392961033722149
70 0.0180459296737036 0.0089658626577325 0.3007769170551054 0.8482619660299474
71 0.0185166857950970 0.0087402885368112 0.3192936028502024 0.8570022545667586
72 0.0189723693827021 0.0085088299643725 0.3382659722329045 0.8655110845311311
73 0.0194111870482028 0.0082716753470887 0.3576771592811073 0.8737827598782199
74 0.0198312691064217 0.0080290355089862 0.3775084283875290 0.8818117953872060
75 0.0202306674413824 0.0077811446451559 0.3977390958289113 0.8895929400323619
76 0.0206073533429372 0.0075282613021386 0.4183464491718485 0.8971212013345005
77 0.0209592153139631 0.0072706693853519 0.4393056644858117 0.9043918707198524
78 0.0212840568481225 0.0070086791939274 0.4605897213339341 0.9114005499137797
79 0.0215795941781907 0.0067426284833258 0.4821693155121248 0.9181431783971056
80 0.0218434539949506 0.0064728835560984 0.5040127695070754 0.9246160619532040
81 0.0220731711366529 0.0061998403811616 0.5260859406437283 0.9308159023343656
82 0.0222661862490424 0.0059239257419534 0.5483521268927706 0.9367398280763191
83 0.0224198434159520 0.0056455984138404 0.5707719703087226 0.9423854264901594
84 0.0225313877604618 0.0053653503711410 0.5933033580691844 0.9477507768613004
85 0.0225979630166251 0.0050837080241352 0.6159013210858095 0.9528344848854355
86 0.0226166090717606 0.0048012334864274 0.6385179301575701 0.9576357183718629
87 0.0225842594793116 0.0045185258730304 0.6611021896368817 0.9621542442448933
88 0.0224977389422703 0.0042362226295390 0.6835999285791520 0.9663904668744322
89 0.0223537607671699 0.0039550008927606 0.7059536893463219 0.9703454677671928
90 0.0221489242886421 0.0036755788831710 0.7281026136349640 0.9740210466503637
91 0.0218797122645414 0.0033987173295630 0.7499823258995053 0.9774197639799267
92 0.0215424882416354 0.0031252209262562 0.7715248141411407 0.9805449849061829
93 0.0211334938918622 0.0028559398232357 0.7926583080330030 0.9834009247294186
94 0.0206488463191534 0.0025917711495875 0.8133071543521564 0.9859926958790061
95 0.0200845353368239 0.0023336605705980 0.8333916896889804 0.9883263564496041
96 0.0194364207155279 0.0020826038788877 0.8528281104045082 0.9904089603284918
97 0.0187002294017812 0.0018396486199437 0.8715283398062894 0.9922486089484355
98 0.0178715527070506 0.0016058957524214 0.8893998925133400 0.9938545047008569
99 0.0169458434674082 0.0013825013435833 0.9063457359807481 0.9952370060444402
100 0.0159184131737539 0.0011706783002407 0.9222641491545019 0.9964076843446809
101 0.0147844290726025 0.0009716981355687 0.9370485782271044 0.9973793824802496
102 0.0135389112374388 0.0007868927721612 0.9505874894645432 0.9981662752524108
103 0.0121767296106375 0.0006176563816932 0.9627642190751807 0.9987839316341039
104 0.0106926010159510 0.0004654472615602 0.9734568200911318 0.9992493788956642
105 0.0090810861415623 0.0003317897488609 0.9825379062326941 0.9995811686445251
106 0.0073365864937051 0.0002182761720913 0.9898744927263992 0.9997994448166164
107 0.0054533413208499 0.0001265688409200 0.9953278340472491 0.9999260136575364
108 0.0034254245084567 0.0000584020744094 0.9987532585557057 0.9999844157319459
109 0.0012467414442944 0.0000155842680537 1.0000000000000002 0.9999999999999996

ndtfast, nfast = 80 109 nfast/ndtfast = 1.36250

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

1.000000000000 1.023119861011 0.511559930505 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.24850

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

Minimum X-grid spacing, DXmin = 1.50000000E+00 km
Maximum X-grid spacing, DXmax = 1.50000000E+00 km
Minimum Y-grid spacing, DYmin = 1.50000000E+00 km
Maximum Y-grid spacing, DYmax = 1.50000000E+00 km
Minimum Z-grid spacing, DZmin = -6.66141633E+00 m
Maximum Z-grid spacing, DZmax = 1.29706643E+03 m

Minimum barotropic Courant Number = 4.22908937E-06
Maximum barotropic Courant Number = 1.50354913E-03
Maximum Coriolis Courant Number = 3.30306627E-05


NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: regional_init_drlu_corr1.nc, Rec=0001, Index=1)
- free-surface
(Min = -1.31885156E-01 Max = 8.38945185E-01)
- vertically integrated u-momentum component
(Min = -4.61772149E-01 Max = 9.95705068E-01)
- vertically integrated v-momentum component
(Min = -3.29916164E-01 Max = 1.17622299E+00)
- u-momentum component
(Min = -1.92722530E+00 Max = 1.09652061E+00)
- v-momentum component
(Min = -4.74239769E-01 Max = 1.63825437E+00)
- potential temperature
(Min = 1.02033235E+00 Max = 3.05137677E+01)
- salinity
(Min = 3.37728616E+01 Max = 3.49289163E+01)
GET_2DFLD - surface u-wind component, t = 0 00:00:00
(Rec=0000001, Index=1, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface v-wind component, t = 0 00:00:00
(Rec=0000001, Index=1, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface air pressure, t = 0 00:00:00
(Rec=0000001, Index=1, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 1.01325000E+03 Max = 1.01325000E+03)
GET_2DFLD - cloud fraction, t = 0 00:00:00
(Rec=0000001, Index=1, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 3.00000012E-01 Max = 3.00000012E-01)
GET_2DFLD - solar shortwave radiation flux, t = 0 00:00:00
(Rec=0000001, Index=1, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface air temperature, t = 0 00:00:00
(Rec=0000001, Index=1, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 2.70000000E+01 Max = 2.70000000E+01)
GET_2DFLD - surface air relative humidity, t = 0 00:00:00
(Rec=0000001, Index=1, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 3.00000000E-01 Max = 3.00000000E-01)
GET_2DFLD - rain fall rate, t = 0 00:00:00
(Rec=0000001, Index=1, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.31885156E-01 Max = 9.00821239E-02)
GET_NGFLD - free-surface eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 4.19598401E-01 Max = 8.38945210E-01)
GET_NGFLD - free-surface southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 4.48445268E-02 Max = 8.15961599E-01)
GET_NGFLD - free-surface northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.31885156E-01 Max = 4.19598401E-01)
GET_NGFLD - 2D u-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 2.28824317E-02 Max = 1.94436431E-01)
GET_NGFLD - 2D v-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -3.11207268E-02 Max = 1.55748472E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.58161715E-01 Max = 3.44468176E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -8.86774212E-02 Max = 3.12169701E-01)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -5.30717000E-02 Max = 4.29459512E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -4.99790981E-02 Max = 3.84470344E-01)
GET_NGFLD - 2D u-momentum northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -2.82011032E-01 Max = 5.95688879E-01)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.64574408E-03 Max = 5.13584197E-01)
GET_NGFLD - 3D u-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -6.17085285E-02 Max = 3.48345220E-01)
GET_NGFLD - 3D v-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -3.22692543E-01 Max = 3.33609760E-01)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -5.79143584E-01 Max = 6.75246298E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -4.26629543E-01 Max = 6.75571918E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.11535013E-01 Max = 5.66532791E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -3.09145749E-01 Max = 1.36818326E+00)
GET_NGFLD - 3D u-momentum northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.51778960E+00 Max = 8.31422031E-01)
GET_NGFLD - 3D v-momentum northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -2.81586587E-01 Max = 7.91704059E-01)
GET_NGFLD - temperature western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 1.91614361E+01 Max = 2.70369129E+01)
GET_NGFLD - salinity western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 3.40383606E+01 Max = 3.44423103E+01)
GET_NGFLD - temperature eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 1.02483130E+00 Max = 2.92381973E+01)
GET_NGFLD - salinity eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 3.41783676E+01 Max = 3.49039421E+01)
GET_NGFLD - temperature southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 1.02198255E+00 Max = 2.92381973E+01)
GET_NGFLD - salinity southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 3.39681854E+01 Max = 3.48956261E+01)
GET_NGFLD - temperature northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 4.22208881E+00 Max = 2.79137688E+01)
GET_NGFLD - salinity northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 3.40822067E+01 Max = 3.48083572E+01)

Basin information for Grid 01:

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

Initial basin volumes: TotVolume = 2.7401513506E+14 m3
MinVolume = 3.9369930300E+03 m3
MaxVolume = 2.9183994588E+09 m3
Max/Min = 7.4127625743E+05

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

GET_2DFLD - surface u-wind component, t = 15 00:00:00
(Rec=0000002, Index=2, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface v-wind component, t = 15 00:00:00
(Rec=0000002, Index=2, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface air pressure, t = 15 00:00:00
(Rec=0000002, Index=2, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 1.01325000E+03 Max = 1.01325000E+03)
GET_2DFLD - cloud fraction, t = 15 00:00:00
(Rec=0000002, Index=2, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 3.00000012E-01 Max = 3.00000012E-01)
GET_2DFLD - solar shortwave radiation flux, t = 15 00:00:00
(Rec=0000002, Index=2, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface air temperature, t = 15 00:00:00
(Rec=0000002, Index=2, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 2.70000000E+01 Max = 2.70000000E+01)
GET_2DFLD - surface air relative humidity, t = 15 00:00:00
(Rec=0000002, Index=2, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 3.00000000E-01 Max = 3.00000000E-01)
GET_2DFLD - rain fall rate, t = 15 00:00:00
(Rec=0000002, Index=2, File: regional_test_frc.nc)
(Tmin= 0.0000 Tmax= 60.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface western boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.52175546E-01 Max = 7.68310875E-02)
GET_NGFLD - free-surface eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 4.30304825E-01 Max = 8.48308504E-01)
GET_NGFLD - free-surface southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 4.04680446E-02 Max = 8.26480865E-01)
GET_NGFLD - free-surface northern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.52175546E-01 Max = 4.30304825E-01)
GET_NGFLD - 2D u-momentum western boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 2.28352156E-02 Max = 1.97847337E-01)
GET_NGFLD - 2D v-momentum western boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -2.89780982E-02 Max = 1.83818042E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.55213594E-01 Max = 3.32948357E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.06036112E-01 Max = 3.06090325E-01)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -5.12533821E-02 Max = 3.72468010E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -4.82416712E-02 Max = 3.37830335E-01)
GET_NGFLD - 2D u-momentum northern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -2.71638721E-01 Max = 5.76307297E-01)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -8.38926528E-03 Max = 5.05541265E-01)
GET_NGFLD - 3D u-momentum western boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -7.09301829E-02 Max = 3.03893089E-01)
GET_NGFLD - 3D v-momentum western boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -3.35006624E-01 Max = 2.82401145E-01)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -5.83098829E-01 Max = 6.95655107E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -4.28113699E-01 Max = 6.74043417E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.17604315E-01 Max = 5.70286274E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -2.87386745E-01 Max = 1.37538004E+00)
GET_NGFLD - 3D u-momentum northern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -1.47611964E+00 Max = 8.97559881E-01)
GET_NGFLD - 3D v-momentum northern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = -2.46511430E-01 Max = 8.37805092E-01)
GET_NGFLD - temperature western boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 1.95922852E+01 Max = 2.71273766E+01)
GET_NGFLD - salinity western boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 3.40337448E+01 Max = 3.44209557E+01)
GET_NGFLD - temperature eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 1.02496529E+00 Max = 2.92330093E+01)
GET_NGFLD - salinity eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 3.41968346E+01 Max = 3.49070244E+01)
GET_NGFLD - temperature southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 1.02208340E+00 Max = 2.92330093E+01)
GET_NGFLD - salinity southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 3.39651375E+01 Max = 3.48979721E+01)
GET_NGFLD - temperature northern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 4.21473789E+00 Max = 2.78516216E+01)
GET_NGFLD - salinity northern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: regional_bry_drlu_corr1.nc)
(Tmin= 0.0000 Tmax= 10.0000)
(Min = 3.40762558E+01 Max = 3.48148994E+01)

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

0 0 00:00:00 1.680682E-02 1.563801E+04 1.563802E+04 2.740898E+14
(026,010,30) 2.011446E-04 5.248219E-04 0.000000E+00 1.923118E+00
DEF_HIS - creating history file, Grid 01: ocean_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
1 0 00:00:00 1.680581E-02 1.563801E+04 1.563802E+04 2.740898E+14
(145,192,01) 1.327778E-04 4.213520E-05 3.403248E-01 1.922544E+00
2 0 00:00:01 1.680549E-02 1.563801E+04 1.563802E+04 2.740898E+14
(023,072,30) 2.122700E-06 5.388221E-05 4.407313E-01 1.921968E+00
3 0 00:00:01 1.680513E-02 1.563801E+04 1.563802E+04 2.740898E+14
(023,072,30) 2.125580E-06 5.383557E-05 4.346372E-01 1.921391E+00
4 0 00:00:02 1.680472E-02 1.563801E+04 1.563802E+04 2.740898E+14
(023,072,30) 2.130057E-06 5.382015E-05 4.286286E-01 1.920812E+00
5 0 00:00:02 1.680425E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.136232E-06 5.381432E-05 4.226885E-01 1.920233E+00
6 0 00:00:03 1.680375E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.143957E-06 5.380407E-05 4.168178E-01 1.919652E+00
7 0 00:00:03 1.680322E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.153128E-06 5.378048E-05 4.110183E-01 1.919070E+00
8 0 00:00:04 1.680265E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.163555E-06 5.371991E-05 4.052909E-01 1.918487E+00
9 0 00:00:04 1.680205E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.174841E-06 5.357774E-05 3.996387E-01 1.917902E+00
10 0 00:00:05 1.680141E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.186463E-06 5.331565E-05 3.940681E-01 1.917316E+00
11 0 00:00:05 1.680075E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.198272E-06 5.296797E-05 3.885821E-01 1.916729E+00
12 0 00:00:06 1.680008E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.210937E-06 5.261222E-05 3.831653E-01 1.916141E+00
13 0 00:00:06 1.679938E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.225278E-06 5.225010E-05 3.777866E-01 1.915551E+00
14 0 00:00:07 1.679867E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.241433E-06 5.184308E-05 3.724261E-01 1.914961E+00
15 0 00:00:07 1.679794E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.259113E-06 5.139675E-05 3.670941E-01 1.914369E+00
16 0 00:00:08 1.679720E-02 1.563801E+04 1.563803E+04 2.740898E+14
(023,072,30) 2.278160E-06 5.093374E-05 3.618236E-01 1.913776E+00
17 0 00:00:08 1.679645E-02 1.563801E+04 1.563803E+04 2.740899E+14
(023,072,30) 2.298625E-06 5.046674E-05 3.566316E-01 1.913181E+00
18 0 00:00:09 1.679568E-02 1.563801E+04 1.563803E+04 2.740899E+14
(023,072,30) 2.320609E-06 5.000391E-05 3.515136E-01 1.912586E+00
19 0 00:00:09 1.679490E-02 1.563801E+04 1.563803E+04 2.740899E+14
(023,072,30) 2.344177E-06 4.955156E-05 3.464624E-01 1.911989E+00
20 0 00:00:10 1.679410E-02 1.563801E+04 1.563803E+04 2.740899E+14
(023,072,30) 2.369369E-06 4.911461E-05 3.414734E-01 1.911391E+00
21 0 00:00:10 1.679330E-02 1.563801E+04 1.563803E+04 2.740899E+14
(023,072,30) 2.396206E-06 4.869705E-05 3.365423E-01 1.910792E+00
22 0 00:00:11 1.679249E-02 1.563801E+04 1.563803E+04 2.740899E+14
(023,072,30) 2.424689E-06 4.830183E-05 3.316657E-01 1.910192E+00
23 0 00:00:11 1.679166E-02 1.563801E+04 1.563803E+04 2.740899E+14
(023,072,30) 2.454811E-06 4.792996E-05 3.268427E-01 1.909591E+00
24 0 00:00:12 1.679082E-02 1.563801E+04 1.563803E+04 2.740899E+14
(023,072,30) 2.486561E-06 4.758014E-05 3.220727E-01 1.908989E+00
25 0 00:00:12 1.678998E-02 1.563801E+04 1.563803E+04 2.740899E+14
(023,072,30) 2.519923E-06 4.724955E-05 3.173542E-01 1.908385E+00
26 0 00:00:13 1.678913E-02 1.563802E+04 1.563803E+04 2.740899E+14
(023,072,30) 2.554880E-06 4.693526E-05 3.126854E-01 1.907781E+00
27 0 00:00:13 1.678827E-02 NaN NaN 2.740899E+14
(023,072,30) 2.591410E-06 4.663499E-05 3.080646E-01 1.907175E+00

Blowing-up: Saving latest model state into RESTART file

WRT_RST - wrote re-start fields (Index=2,2) into time record = 0000001

Elapsed CPU time (seconds):

Node # 0 CPU: 51.395
Node # 1 CPU: 52.451
Node # 2 CPU: 52.451
Node # 3 CPU: 52.451
Node # 4 CPU: 52.451
Node # 5 CPU: 52.451
Node # 6 CPU: 52.451
Node # 7 CPU: 52.447

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

Analytical header files used:

ROMS/Functionals/ana_btflux.h

ROMS/TOMS: DONE... Wednesday - January 7, 2015 - 9:34:30 PM



I previously thought the problem had to do with some grossly unrealistic values generated when I created the forcing files. But I have checked for these and found none. Hence, I would greatly appreciate any advice or suggestions on how to identify the problems. Thank you very much!

Best regards,
Lawrence

[kate]

This is worrisome:

Maximum grid stiffness ratios: rx0 = 9.993076E-01 (Beckmann and Haidvogel)

Perhaps you should try some bathymetry smoothing. It really does make a difference to model stability.

[lcbernardo]

Hi Kate,

Thank you so much for pointing this out, and I'll remember it as I continue my learning. I'll attempt some bathymetry smoothing and see what difference it can make.

In relation to this, when I was simulating smaller domains using other models, I would start with a relatively "simple" configuration, such as tidally forced open boundaries and homogenous temperature and salinity profiles. Using this simplified configuration, I could identify problems related to bathymetry. For my present regional case, I immediately moved into creating 3D initial and boundary conditions, but is there a simpler configuration I can try to have a good chance of detecting bathymetry related problems early on?

Lawrence

[d.kobashi]

A few comments aside from what Kate suggested.

1. Turn off Wet/Dry and see if the model still blows up.
2. You said you turned on TS_SADVECTION, but your ROMS log says TS_C4VADVECTION. So I am confused. Don't use TS_SADVECTION.
3. Try other GLS mixing like k-omega, not k-e.
4. You max courant number is too low. You are wasting CPU resources. This will not blow up the model though.
5. Too much bathy smoothing will compromise the accuracy of the model so be careful. Turn on UV_DIAGNOSTICS and TS_DIAGNOSTICS. In the diagnostic file, you can see 3D pressure gradient terms. If the bottom values are too high, pressure gradient error may be high. you may want to reduce it. Your rx1 value looks ok, so I am not sure if that's the case, though.

Hope this helps. Good luck.

DJ Kobashi@TAMU

[lcbernardo]

Hi DJ,

Thank you for the additional comments and suggestions. So far, doing some bathymetric smoothing seems to have changed the results somewhat, but indeed, I have to be careful about overdoing it. I'll also turn on the diagnostics for my succeeding tests.

I also want to maximize the available computational resources, so I appreciate your comment about the Courant number. And since I'm working on a relatively regional scale, it does makes sense to turn off the wetting and drying for now. Thanks again!

Lawrence

[lanerolle]

I can suggest:

1. Not using the spline-based vertical advection for momentum & tracers but instead using say the default 4th order centered schemes
2. Switching off the SPLINE CPP option
3. Reducing the time step
4. Changing the baroclinic:barotropic splitting ratio (NDTFAST) - I usually use 20, 30 or 40

[lcbernardo]

Hi lanerolle,

Thank you for your suggestions. I'll try not using the spline-related schemes in one of my test configurations. As for the timestep, I was of the thinking that it was already very small given that the horizontal grid resolution is greater than 1 km. But maybe it is necessary to handle the bathymetric complexity? I had my NDTFAST orginally at 20, and tried to increase it to see if it would help solve the problem. I will decrease it to one of the values you suggested.

Lawrence

[kate]

On the timestep (both dt and ndtfast), you try things and then if it's stable, try a longer timestep. If it's unstable, try a shorter timestep. You might be right at the edge, in which it runs for thousands of steps, then blows up. Sometimes you can get through those by using a shorter timestep for a few days, then go back to what you had before. The Courant numbers are there for guidance, but of course depend on the instantaneous currants.

[lcbernardo]

Hi Kate,

Thank you for the guidance on selecting the proper timestep. In numerical models I've tried before, ocean or other, selecting the proper timestep indeed involved something similar to what you have suggested. What's new to me is the use of not one but two (barotropic and baroclinic), which I will need to read up upon.

Lawrence

[d.kobashi]

Splitting time step is a way to run a full 3D ocean model more computationally efficiently as baroclinic mode is much slower than barotropic mode. It is used for recent ocean models: structured, unstructured models alike. GETM, MOM, FVCOM etc.

By the way, what I suggested was to identify the problem you are facing, not to make the model work.
You included lots of options. Ideally, you need to check them one by one to identify what option or input data cause the problem. Fortunately, there is lots of information on this forum and many of the problems were identified by ROMS users and developers in the past. So this forum does help a lot to a newbie (I still can call myself a newbie compared to other experienced modelers.....).

When ROMS blows up, the first thing to do (at least to me) is to check restart file and see if there is any suspicious values (e.g. current speed higher than 5 m/s, water level higher than 20 m etc.). And also check input data to see if there are any NaNs etc. in the input files. Debug mode also often helps you identify problems (-g in fortran option).

Good luck

DJ Kobashi@TAMU

[lcbernardo]

Hi DJ and all,

Thank you for the explanation of the time-splitting feature. I'll try to look more into the details as I become more familiar with the model.

I'd like to report that with my recent runs, the early blowing up problem stopped occurring after doing some bathymetry smoothing. Now I'm trying to check the differences when turning some features on or off, like wetting and drying and the splines. Of course, I continue to improve the bathymetric data set in such a way that the important features are still well-reproduced in spite of the smoothing.

As for the forums, I have indeed found the discussions helpful as I go along, and the bookmark feature has really come in handy as I continue to learn. I think learning how to implement tides will come next for me.

I'm also now looking into the proper boundary conditions. And if I may throw in an additional question, for the case wherein you have boundary forcing data for all relevant fields (zeta, u, v, ubar, vbar, salt, temperature), what would be the best combination of boundary conditions to use? Following some of the other users in the forums, I have chosen a combination of Chapman and Flather, with Rad+Nudging as in the following:

LBC(isFsur) == Cha Cha Cha Cha ! free-surface
LBC(isUbar) == Fla Fla Fla Fla ! 2D U-momentum
LBC(isVbar) == Fla Fla Fla Fla ! 2D V-momentum
LBC(isUvel) == RadNud RadNud RadNud RadNud ! 3D U-momentum
LBC(isVvel) == RadNud RadNud RadNud RadNud ! 3D V-momentum
LBC(isMtke) == Gra Gra Gra Gra ! mixing TKE

LBC(isTvar) == RadNud RadNud RadNud RadNud \ ! temperature
RadNud RadNud RadNud RadNud ! salinity

However, would it make more sense to use the clamped boundary condition for all fields, given that the input data is quite "complete"?

Once again, thank you all for your kind assistance.

Lawrence

[lizhicheng]

I try to add WET_DRY cpp option, and smooth the topography, but the model early blowing up? The rx0 is high?

Basin information for Grid 01:

Maximum grid stiffness ratios: rx0 = 1.547781E+03 (Beckmann and Haidvogel)
rx1 = 1.229183E+05 (Haney)

Initial domain volumes: TotVolume = -4.2299594144E+13 m3
MinCellVol = -9.0403834918E+08 m3
MaxCellVol = 1.6578109192E+08 m3
Max/Min = -1.8337838441E-01

[lcbernardo]

Hi Lizhicheng,

Your rx0 and rx1 values are still too high I think. But I think it would be better to show your complete log file to see what other possible problems the model may have.

As for bathymetry smoothing, you can try to use GridBuilder (https://austides.com/downloads/) or other tools which allow you to specify your target rx0 and rx1 values.

Lawrence

[kate]

When doing wetting and drying, it is possible to have huge rx0 numbers in the wet-dry zone and still be OK. It is also possible to be in trouble there, so you need to really know what you are doing. I know that a "dry" cell 1 m above sea level next to a fjord can spell trouble.

There might still be code tweaks that John Warner and I have that Hernan doesn't for making wetting and drying behave better.