MAIN: Abnormal termination: BLOWUP.

Message

Lee2015 · #1 Unread post by **Lee2015** » Wed Jun 08, 2016 2:59 am

I got a BLOWUP error, as follows. Do you guys have any ideas about it? Thanks alot!

Model Input Parameters: ROMS/TOMS version 3.7
Wednesday - June 8, 2016 - 6:37:03 PM
-----------------------------------------------------------------------------

Mexico

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -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 : 783M

Local Root : /home/lee/ROMS/src
Header Dir : /home/lee/ROMS/Projects/mexico
Header file : mexico.h
Analytical Dir: /home/lee/ROMS/Projects/mexico

Resolution, Grid 01: 0179x0134x015, Parallel Nodes: 8, Tiling: 002x004

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

8784 ntimes Number of timesteps for 3-D equations.
600.000 dt Timestep size (s) for 3-D equations.
20 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.
F LcycleRST Switch to recycle time-records in restart file.
144000 nRST Number of timesteps between the writing of data
into restart fields.
6 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
144000 nHIS Number of timesteps between the writing fields
into history file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
144000 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
5.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
5.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
1.0000E+02 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.
2.6500E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
1 lmd_Jwt Jerlov water type.
2 Vtransform S-coordinate transformation equation.
2 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.
100.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).
20130801.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.
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.
10.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
-1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
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.
T LtracerCLM(01) Turning ON processing of climatology tracer 01: temp
T LtracerCLM(02) Turning ON 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 Hout(idUsms) Write out surface U-momentum stress.
T Hout(idVsms) Write out surface V-momentum stress.
T Hout(idTsur) Write out surface net heat flux.
T Hout(idTsur) Write out surface net salt flux.
T Hout(idSrad) Write out shortwave radiation flux.
T Hout(idDano) Write out density anomaly.
T Hout(idHsbl) Write out depth of surface boundary layer.

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

Output/Input Files:

Output Restart File: Out/mexico_rst.nc
Output Averages File: Out/mexico_avg.nc
Input Grid File: In/t03031_smth9_grid.nc
Input Nonlinear Initial File: In/damee4_levfeb_b.nc
Input Forcing File 01: In/gom_shflux_era.nc
Input Forcing File 02: In/gom_sms_era.nc
Input Forcing File 03: In/gom_swflux_era.nc
Input Forcing File 04: In/gom_swrad_era.nc
Input Forcing File 05: In/gom_wind_era.nc
Input Climatology File: In/damee4_Lclm_b.nc
Input Boundary File: In/t0604_bry.nc

Tile partition information for Grid 01: 0179x0134x0015 tiling: 002x004

tile Istr Iend Jstr Jend Npts

0 1 90 1 33 44550
1 91 179 1 33 44055
2 1 90 34 67 45900
3 91 179 34 67 45390
4 1 90 68 101 45900
5 91 179 68 101 45390
6 1 90 102 134 44550
7 91 179 102 134 44055

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 90.50 0.50 33.50 RHO-points
1 90.50 179.50 0.50 33.50 RHO-points
2 0.50 90.50 33.50 67.50 RHO-points
3 90.50 179.50 33.50 67.50 RHO-points
4 0.50 90.50 67.50 101.50 RHO-points
5 90.50 179.50 67.50 101.50 RHO-points
6 0.50 90.50 101.50 134.50 RHO-points
7 90.50 179.50 101.50 134.50 RHO-points

0 1.00 90.50 0.50 33.50 U-points
1 90.50 179.00 0.50 33.50 U-points
2 1.00 90.50 33.50 67.50 U-points
3 90.50 179.00 33.50 67.50 U-points
4 1.00 90.50 67.50 101.50 U-points
5 90.50 179.00 67.50 101.50 U-points
6 1.00 90.50 101.50 134.50 U-points
7 90.50 179.00 101.50 134.50 U-points

0 0.50 90.50 1.00 33.50 V-points
1 90.50 179.50 1.00 33.50 V-points
2 0.50 90.50 33.50 67.50 V-points
3 90.50 179.50 33.50 67.50 V-points
4 0.50 90.50 67.50 101.50 V-points
5 90.50 179.50 67.50 101.50 V-points
6 0.50 90.50 101.50 134.00 V-points
7 90.50 179.50 101.50 134.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 306
HaloSizeJ(1) = 144
TileSide(1) = 96
TileSize(1) = 4032

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

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

zeta 1 Closed Gradient Gradient Closed

ubar 1 Closed Clamped Clamped Closed

vbar 1 Closed Clamped Clamped Closed

u 1 Closed Clamped Clamped Closed

v 1 Closed Clamped Clamped Closed

temp 1 Closed Clamped Clamped Closed

salt 1 Closed Clamped Clamped Closed

Activated C-preprocessing Options:

MEXICO Mexico
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
CURVGRID Orthogonal curvilinear grid.
DIFF_GRID Horizontal diffusion coefficient scaled by grid size.
DIURNAL_SRFLUX Modulate shortwave radiation by the local diurnal cycle.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
LMD_CONVEC LMD convective mixing due to shear instability.
LMD_MIXING Large/McWilliams/Doney interior mixing.
LMD_NONLOCAL LMD convective nonlocal transport.
LMD_RIMIX LMD diffusivity due to shear instability.
LMD_SKPP KPP surface boundary layer mixing.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
PERFECT_RESTART Processing perfect restart variables.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
RI_SPLINES Parabolic Spline Reconstruction for Richardson Number.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
SPLINES_VDIFF Parabolic Spline Reconstruction for Vertical Diffusion.
SPLINES_VVISC Parabolic Spline Reconstruction for Vertical Viscosity.
TS_MPDATA Recursive flux corrected MPDATA 3D 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_QDRAG Quadratic bottom stress.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
VISC_GRID Horizontal viscosity coefficient scaled by grid size.

Process Information:

Node # 0 (pid= 12099) is active.
Node # 4 (pid= 12103) is active.
Node # 6 (pid= 12105) is active.
Node # 2 (pid= 12101) is active.
Node # 3 (pid= 12102) is active.

INITIAL: Configuring and initializing forward nonlinear model ...
*******
Node # 1 (pid= 12100) is active.
Node # 5 (pid= 12104) is active.
Node # 7 (pid= 12106) is active.

Vertical S-coordinate System, Grid 01:

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

15 0.0000000 0.0000000 -0.000 0.000 0.000 0.000
14 -0.0666667 -0.0004989 2.218 -3.358 -7.757 -9.257
13 -0.1333333 -0.0032168 4.418 -6.828 -21.407 -30.564
12 -0.2000000 -0.0100379 6.583 -10.502 -45.951 -74.151
11 -0.2666667 -0.0229247 8.698 -14.480 -86.600 -150.675
10 -0.3333333 -0.0439342 10.745 -18.863 -148.816 -271.306
9 -0.4000000 -0.0752335 12.706 -23.762 -238.353 -447.812
8 -0.4666667 -0.1191055 14.563 -29.289 -361.273 -692.587
7 -0.5333333 -0.1779212 16.295 -35.563 -523.870 -1018.508
6 -0.6000000 -0.2540303 17.883 -42.702 -732.384 -1438.333
5 -0.6666667 -0.3494770 19.310 -50.807 -992.242 -1963.162
4 -0.7333333 -0.4653579 20.566 -59.935 -1306.355 -2598.950
3 -0.8000000 -0.6004851 21.663 -70.024 -1671.571 -3339.247
2 -0.8666667 -0.7487245 22.650 -80.770 -2071.601 -4150.744
1 -0.9333333 -0.8938704 23.662 -91.360 -2463.418 -4945.443
0 -1.0000000 -1.0000000 25.000 -100.000 -2751.640 -5528.281

Time Splitting Weights for Grid 01: ndtfast = 20 nfast = 29
==================================

Primary Secondary Accumulated to Current Step

1-0.0009651193358779 0.0500000000000000-0.0009651193358779 0.0500000000000000
2-0.0013488780126037 0.0500482559667939-0.0023139973484816 0.1000482559667939
3-0.0011514592651644 0.0501156998674241-0.0034654566136460 0.1501639558342180
4-0.0003735756740661 0.0501732728306823-0.0038390322877122 0.2003372286649003
5 0.0009829200513762 0.0501919516143856-0.0028561122363360 0.2505291802792859
6 0.0029141799764308 0.0501428056118168 0.0000580677400949 0.3006719858911027
7 0.0054132615310267 0.0499970966129953 0.0054713292711215 0.3506690825040980
8 0.0084687837865133 0.0497264335364439 0.0139401130576348 0.4003955160405419
9 0.0120633394191050 0.0493029943471183 0.0260034524767398 0.4496985103876601
10 0.0161716623600090 0.0486998273761630 0.0421751148367488 0.4983983377638231
11 0.0207585511322367 0.0478912442581626 0.0629336659689855 0.5462895820219857
12 0.0257765478740990 0.0468533167015507 0.0887102138430846 0.5931428987235364
13 0.0311633730493853 0.0455644893078458 0.1198735868924699 0.6387073880313822
14 0.0368391158442262 0.0440063206553765 0.1567127027366961 0.6827137086867586
15 0.0427031802506397 0.0421643648631652 0.1994158829873358 0.7248780735499238
16 0.0486309868367616 0.0400292058506332 0.2480468698240974 0.7649072794005570
17 0.0544704302037592 0.0375976565087951 0.3025173000278565 0.8025049359093521
18 0.0600380921294286 0.0348741349986072 0.3625553921572851 0.8373790709079593
19 0.0651152103984763 0.0318722303921357 0.4276706025557614 0.8692513013000950
20 0.0694434033194839 0.0286164698722119 0.4971140058752453 0.8978677711723069
21 0.0727201499285570 0.0251442997062377 0.5698341558038023 0.9230120708785446
22 0.0745940258796570 0.0215082922098099 0.6444281816834592 0.9445203630883545
23 0.0746596950216179 0.0177785909158270 0.7190878767050771 0.9622989540041815
24 0.0724526566618460 0.0140456061647461 0.7915405333669231 0.9763445601689277
25 0.0674437485167025 0.0104229733316538 0.8589842818836255 0.9867675335005816
26 0.0590334053485720 0.0070507859058187 0.9180176872321975 0.9938183194064003
27 0.0465456732896125 0.0040991156383901 0.9645633605218099 0.9979174350447905
28 0.0292219798521905 0.0017718319739095 0.9937853403740003 0.9996892670187000
29 0.0062146596259993 0.0003107329813000 0.9999999999999997 1.0000000000000000

ndtfast, nfast = 20 29 nfast/ndtfast = 1.45000

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

1.000000000000 1.060707743385 0.530353871693 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.14200

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

Minimum X-grid spacing, DXmin = 9.53173644E+00 km
Maximum X-grid spacing, DXmax = 1.06053729E+01 km
Minimum Y-grid spacing, DYmin = 1.11200424E+01 km
Maximum Y-grid spacing, DYmax = 1.11200424E+01 km
Minimum Z-grid spacing, DZmin = -2.21806494E+00 m
Maximum Z-grid spacing, DZmax = 8.11496867E+02 m

Minimum barotropic Courant Number = 5.49001381E-02
Maximum barotropic Courant Number = 9.10581682E-01
Maximum Coriolis Courant Number = 4.45309430E-02

Horizontal mixing scaled by grid size, GRDMAX = 1.08596592E+01 km

NLM: GET_STATE - Read state initial conditions, t = 1 00:00:00
(Grid 01, File: damee4_levfeb_b.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 = 0.00000000E+00 Max = 3.00944372E+01)
- salinity
(Min = 0.00000000E+00 Max = 3.68434650E+01)
GET_2DFLD - surface u-momentum stress, t = 0 22:30:00
(Rec=0000008, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -5.84637914E-04 Max = 1.22592792E-04)
GET_2DFLD - surface v-momentum stress, t = 0 22:30:00
(Rec=0000008, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -2.31960137E-04 Max = 2.06426426E-04)
GET_2DFLD - solar shortwave radiation flux, t = 0 22:30:00
(Rec=0000008, Index=1, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = 1.51806455E-05 Max = 1.16572371E-04)
GET_2DFLD - surface net heat flux, t = 0 22:30:00
(Rec=0000008, Index=1, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -9.66121467E-06 Max = 8.56924329E-05)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -7.70141743E-03 Max = 2.82108545E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.10545473E-02 Max = 5.74858561E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.16785280E-01 Max = 2.00956576E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.10545473E-02 Max = 7.94643536E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.47733809E-02 Max = 3.51793170E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.24065101E-01 Max = 8.88886154E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.89331788E-01 Max = 5.45626506E-02)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.17090203E-01 Max = 2.46778369E-01)
GET_NGFLD - temperature eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.09180876E-13 Max = 3.00770473E+01)
GET_NGFLD - salinity eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.34068986E-13 Max = 3.65945892E+01)
GET_NGFLD - temperature southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.54044008E-12 Max = 2.94915771E+01)
GET_NGFLD - salinity southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.91542485E-12 Max = 3.66457024E+01)
GET_3DFLD - potential temperature climatology, t = 0 00:00:00
(Rec=0000001, Index=1, File: damee4_Lclm_b.nc)
(Tmin= 0.0000 Tmax= 334.0000)
(Min = 0.00000000E+00 Max = 2.74143713E+01)
GET_3DFLD - salinity climatology, t = 0 00:00:00
(Rec=0000001, Index=1, File: damee4_Lclm_b.nc)
(Tmin= 0.0000 Tmax= 334.0000)
(Min = 0.00000000E+00 Max = 3.67458003E+01)

Basin information for Grid 01:

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

Initial basin volumes: TotVolume = 3.4615820664E+15 m3
MinVolume = 1.1843849025E+08 m3
MaxVolume = 9.5648797369E+10 m3
Max/Min = 8.0758203831E+02

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000145 - 00008784)

GET_2DFLD - surface u-momentum stress, t = 1 01:30:00
(Rec=0000009, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -3.85218911E-04 Max = 1.04836739E-04)
GET_2DFLD - surface v-momentum stress, t = 1 01:30:00
(Rec=0000009, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -1.19083844E-04 Max = 2.21327180E-04)
GET_2DFLD - solar shortwave radiation flux, t = 1 01:30:00
(Rec=0000009, Index=2, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -4.22234950E-20 Max = 6.54829338E-06)
GET_2DFLD - surface net heat flux, t = 1 01:30:00
(Rec=0000009, Index=2, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -6.49440989E-05 Max = -4.62032956E-06)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.88156872E-03 Max = 2.52899200E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.15718057E-02 Max = 6.22464232E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.28449917E-01 Max = 1.85216349E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.15718057E-02 Max = 9.02452767E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.58911180E-02 Max = 3.12854528E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.16634078E-01 Max = 7.79765025E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.83034706E-01 Max = 3.61946411E-02)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.13600455E-01 Max = 2.92638898E-01)
GET_NGFLD - temperature eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.08041814E-13 Max = 2.99027557E+01)
GET_NGFLD - salinity eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.34010629E-13 Max = 3.65638466E+01)
GET_NGFLD - temperature southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.51450098E-12 Max = 2.96635704E+01)
GET_NGFLD - salinity southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: t0604_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.91184872E-12 Max = 3.66002235E+01)
GET_3DFLD - potential temperature climatology, t = 31 00:00:00
(Rec=0000002, Index=2, File: damee4_Lclm_b.nc)
(Tmin= 0.0000 Tmax= 334.0000)
(Min = 0.00000000E+00 Max = 2.72762837E+01)
GET_3DFLD - salinity climatology, t = 31 00:00:00
(Rec=0000002, Index=2, File: damee4_Lclm_b.nc)
(Tmin= 0.0000 Tmax= 334.0000)
(Min = 0.00000000E+00 Max = 3.67549533E+01)

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

144 1 00:00:00 1.101055E-05 1.555875E+04 1.555875E+04 3.440743E+15
(179,065,13) 1.035187E-02 0.000000E+00 0.000000E+00 2.478671E-01
150 1 01:00:00 2.369821E-02 1.555558E+04 1.555560E+04 3.440761E+15
(108,010,15) 2.492003E-02 4.008041E-02 2.527212E+00 3.866112E+00
GET_2DFLD - surface u-momentum stress, t = 1 04:30:00
(Rec=0000010, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -1.63608568E-04 Max = 3.22725184E-05)
GET_2DFLD - surface v-momentum stress, t = 1 04:30:00
(Rec=0000010, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -7.96744387E-05 Max = 2.31890816E-04)
GET_2DFLD - solar shortwave radiation flux, t = 1 04:30:00
(Rec=0000010, Index=1, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -9.37076875E-09 Max = 9.37076875E-09)
GET_2DFLD - surface net heat flux, t = 1 04:30:00
(Rec=0000010, Index=1, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -6.52646996E-05 Max = -4.33468254E-06)
156 1 02:00:00 8.555982E-02 1.555581E+04 1.555589E+04 3.440779E+15
(045,007,15) 1.897565E-03 1.768342E-02 1.276331E+00 5.946715E+00
162 1 03:00:00 2.754587E-01 1.555611E+04 1.555638E+04 3.440796E+15
(045,040,08) 1.405551E-02 1.228169E-02 2.865928E+00 1.217377E+01
168 1 04:00:00 NaN NaN NaN NaN
(090,000,01) NaN NaN NaN 0.000000E+00

Elapsed CPU time (seconds):

Node # 0 CPU: 48.228
Node # 4 CPU: 49.788
Node # 1 CPU: 49.612
Node # 2 CPU: 49.964
Node # 5 CPU: 49.464
Node # 6 CPU: 49.516
Node # 3 CPU: 49.580
Node # 7 CPU: 49.664
Total: 395.816

Nonlinear model elapsed time profile, Grid: 01

Allocation and array initialization .............. 0.088 ( 0.0222 %)
Ocean state initialization ....................... 0.732 ( 0.1849 %)
Reading of input data ............................ 21.804 ( 5.5086 %)
Processing of input data ......................... 13.304 ( 3.3612 %)
Processing of output time averaged data .......... 0.420 ( 0.1061 %)
Computation of vertical boundary conditions ...... 1.424 ( 0.3598 %)
Computation of global information integrals ...... 1.048 ( 0.2648 %)
Model 2D kernel .................................. 295.088 (74.5518 %)
2D/3D coupling, vertical metrics ................. 9.388 ( 2.3718 %)
Omega vertical velocity .......................... 3.288 ( 0.8307 %)
Equation of state for seawater ................... 7.680 ( 1.9403 %)
KPP vertical mixing parameterization ............. 6.564 ( 1.6583 %)
3D equations right-side terms .................... 0.376 ( 0.0950 %)
3D equations predictor step ...................... 3.280 ( 0.8287 %)
Pressure gradient ................................ 0.312 ( 0.0788 %)
Harmonic mixing of tracers, geopotentials ........ 0.384 ( 0.0970 %)
Harmonic stress tensor, S-surfaces ............... 0.168 ( 0.0424 %)
Corrector time-step for 3D momentum .............. 4.528 ( 1.1440 %)
Corrector time-step for tracers .................. 10.372 ( 2.6204 %)
Total: 380.248 96.0669

Nonlinear model message Passage profile, Grid: 01

Message Passage: 2D halo exchanges ............... 311.556 (78.7123 %)
Message Passage: 3D halo exchanges ............... 21.480 ( 5.4268 %)
Message Passage: 4D halo exchanges ............... 9.604 ( 2.4264 %)
Message Passage: data broadcast .................. 18.972 ( 4.7931 %)
Message Passage: data reduction .................. 3.028 ( 0.7650 %)
Message Passage: data scattering.................. 11.032 ( 2.7872 %)
Total: 375.672 94.9108

All percentages are with respect to total time = 395.816

ROMS/TOMS - Output NetCDF summary for Grid 01:

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_srflux.h
ROMS/Functionals/ana_stflux.h

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

MAIN: Abnormal termination: BLOWUP.

wilkin · #2 Unread post by **wilkin** » Wed Jun 08, 2016 11:18 am

Check your grid. It looks like there is a point somewhere where the water depth is 25 m above 0:

Vertical S-coordinate System, Grid 01:

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

15 0.0000000 0.0000000 -0.000 0.000 0.000 0.000
14 -0.0666667 -0.0004989 2.218 -3.358 -7.757 -9.257
13 -0.1333333 -0.0032168 4.418 -6.828 -21.407 -30.564
...
1 -0.9333333 -0.8938704 23.662 -91.360 -2463.418 -4945.443
0 -1.0000000 -1.0000000 25.000 -100.000 -2751.640 -5528.281

This might explain why you have a minimum dz that is actually negative.

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

Minimum X-grid spacing, DXmin = 9.53173644E+00 km
Maximum X-grid spacing, DXmax = 1.06053729E+01 km
Minimum Y-grid spacing, DYmin = 1.11200424E+01 km
Maximum Y-grid spacing, DYmax = 1.11200424E+01 km
Minimum Z-grid spacing, DZmin = -2.21806494E+00 m
Maximum Z-grid spacing, DZmax = 8.11496867E+02 m

Lee2015 · #3 Unread post by **Lee2015** » Sat Jun 11, 2016 3:41 am

Sorry for late reply. Thanks! I will check and fix it.

Lee2015 · #4 Unread post by **Lee2015** » Sat Jun 11, 2016 10:15 am

wilkin wrote:Check your grid. It looks like there is a point somewhere where the water depth is 25 m above 0:

Vertical S-coordinate System, Grid 01:

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

15 0.0000000 0.0000000 -0.000 0.000 0.000 0.000
14 -0.0666667 -0.0004989 2.218 -3.358 -7.757 -9.257
13 -0.1333333 -0.0032168 4.418 -6.828 -21.407 -30.564
...
1 -0.9333333 -0.8938704 23.662 -91.360 -2463.418 -4945.443
0 -1.0000000 -1.0000000 25.000 -100.000 -2751.640 -5528.281
This might explain why you have a minimum dz that is actually negative.

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

Minimum X-grid spacing, DXmin = 9.53173644E+00 km
Maximum X-grid spacing, DXmax = 1.06053729E+01 km
Minimum Y-grid spacing, DYmin = 1.11200424E+01 km
Maximum Y-grid spacing, DYmax = 1.11200424E+01 km
Minimum Z-grid spacing, DZmin = -2.21806494E+00 m
Maximum Z-grid spacing, DZmax = 8.11496867E+02 m

Hi wilkin,

I have modify 25 to 0, but the same blowup error still exist.

Model Input Parameters: ROMS/TOMS version 3.7
Saturday - June 11, 2016 - 6:16:21 PM
-----------------------------------------------------------------------------

Mexico

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -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 : 783M

Local Root : /home/lee/ROMS/src
Header Dir : /home/lee/ROMS/Projects/mexico
Header file : mexico.h
Analytical Dir: /home/lee/ROMS/Projects/mexico

Resolution, Grid 01: 0179x0134x015, Parallel Nodes: 8, Tiling: 002x004

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

8784 ntimes Number of timesteps for 3-D equations.
600.000 dt Timestep size (s) for 3-D equations.
20 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.
F LcycleRST Switch to recycle time-records in restart file.
144000 nRST Number of timesteps between the writing of data
into restart fields.
6 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
144000 nHIS Number of timesteps between the writing fields
into history file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
144000 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
5.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
5.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
1.0000E+02 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.
2.6500E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
1 lmd_Jwt Jerlov water type.
2 Vtransform S-coordinate transformation equation.
2 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.
100.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).
20130801.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.
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.
10.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
-1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
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.
T LtracerCLM(01) Turning ON processing of climatology tracer 01: temp
T LtracerCLM(02) Turning ON 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 Hout(idUsms) Write out surface U-momentum stress.
T Hout(idVsms) Write out surface V-momentum stress.
T Hout(idTsur) Write out surface net heat flux.
T Hout(idTsur) Write out surface net salt flux.
T Hout(idSrad) Write out shortwave radiation flux.
T Hout(idDano) Write out density anomaly.
T Hout(idHsbl) Write out depth of surface boundary layer.

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

Output/Input Files:

Output Restart File: Out/mexico_rst.nc
Output Averages File: Out/mexico_avg.nc
Input Grid File: In/t03031_smth9_grid.nc
Input Nonlinear Initial File: In/gom_ini.nc
Input Forcing File 01: In/gom_shflux_era.nc
Input Forcing File 02: In/gom_sms_era.nc
Input Forcing File 03: In/gom_swflux_era.nc
Input Forcing File 04: In/gom_swrad_era.nc
Input Forcing File 05: In/gom_wind_era.nc
Input Climatology File: In/gom_clm.nc
Input Boundary File: In/gom_bry.nc

Tile partition information for Grid 01: 0179x0134x0015 tiling: 002x004

tile Istr Iend Jstr Jend Npts

0 1 90 1 33 44550
1 91 179 1 33 44055
2 1 90 34 67 45900
3 91 179 34 67 45390
4 1 90 68 101 45900
5 91 179 68 101 45390
6 1 90 102 134 44550
7 91 179 102 134 44055

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 90.50 0.50 33.50 RHO-points
1 90.50 179.50 0.50 33.50 RHO-points
2 0.50 90.50 33.50 67.50 RHO-points
3 90.50 179.50 33.50 67.50 RHO-points
4 0.50 90.50 67.50 101.50 RHO-points
5 90.50 179.50 67.50 101.50 RHO-points
6 0.50 90.50 101.50 134.50 RHO-points
7 90.50 179.50 101.50 134.50 RHO-points

0 1.00 90.50 0.50 33.50 U-points
1 90.50 179.00 0.50 33.50 U-points
2 1.00 90.50 33.50 67.50 U-points
3 90.50 179.00 33.50 67.50 U-points
4 1.00 90.50 67.50 101.50 U-points
5 90.50 179.00 67.50 101.50 U-points
6 1.00 90.50 101.50 134.50 U-points
7 90.50 179.00 101.50 134.50 U-points

0 0.50 90.50 1.00 33.50 V-points
1 90.50 179.50 1.00 33.50 V-points
2 0.50 90.50 33.50 67.50 V-points
3 90.50 179.50 33.50 67.50 V-points
4 0.50 90.50 67.50 101.50 V-points
5 90.50 179.50 67.50 101.50 V-points
6 0.50 90.50 101.50 134.00 V-points
7 90.50 179.50 101.50 134.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 306
HaloSizeJ(1) = 144
TileSide(1) = 96
TileSize(1) = 4032

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

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

zeta 1 Closed Gradient Gradient Closed

ubar 1 Closed Clamped Clamped Closed

vbar 1 Closed Clamped Clamped Closed

u 1 Closed Clamped Clamped Closed

v 1 Closed Clamped Clamped Closed

temp 1 Closed Clamped Clamped Closed

salt 1 Closed Clamped Clamped Closed

Activated C-preprocessing Options:

MEXICO Mexico
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
CURVGRID Orthogonal curvilinear grid.
DIFF_GRID Horizontal diffusion coefficient scaled by grid size.
DIURNAL_SRFLUX Modulate shortwave radiation by the local diurnal cycle.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
LMD_CONVEC LMD convective mixing due to shear instability.
LMD_MIXING Large/McWilliams/Doney interior mixing.
LMD_NONLOCAL LMD convective nonlocal transport.
LMD_RIMIX LMD diffusivity due to shear instability.
LMD_SKPP KPP surface boundary layer mixing.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
PERFECT_RESTART Processing perfect restart variables.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
RI_SPLINES Parabolic Spline Reconstruction for Richardson Number.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
SPLINES_VDIFF Parabolic Spline Reconstruction for Vertical Diffusion.
SPLINES_VVISC Parabolic Spline Reconstruction for Vertical Viscosity.
TS_MPDATA Recursive flux corrected MPDATA 3D 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_QDRAG Quadratic bottom stress.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
VISC_GRID Horizontal viscosity coefficient scaled by grid size.

Process Information:

Node # 0 (pid= 2254) is active.
Node # 4 (pid= 2258) is active.
Node # 5 (pid= 2259) is active.

Node # 2 (pid= 2256) is active.
INITIAL: Configuring and initializing forward nonlinear model ...
*******
Node # 1 (pid= 2255) is active.
Node # 3 (pid= 2257) is active.
Node # 6 (pid= 2260) is active.
Node # 7 (pid= 2261) is active.

Vertical S-coordinate System, Grid 01:

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

15 0.0000000 0.0000000 0.000 0.000 0.000 0.000
14 -0.0666667 -0.0004989 -0.000 -3.358 -7.765 -9.257
13 -0.1333333 -0.0032168 -0.000 -6.828 -21.449 -30.564
12 -0.2000000 -0.0100379 -0.000 -10.502 -46.079 -74.151
11 -0.2666667 -0.0229247 -0.000 -14.480 -86.890 -150.675
10 -0.3333333 -0.0439342 -0.000 -18.863 -149.370 -271.306
9 -0.4000000 -0.0752335 -0.000 -23.762 -239.299 -447.812
8 -0.4666667 -0.1191055 -0.000 -29.289 -362.767 -692.587
7 -0.5333333 -0.1779212 -0.000 -35.563 -526.099 -1018.508
6 -0.6000000 -0.2540303 -0.000 -42.702 -735.565 -1438.333
5 -0.6666667 -0.3494770 -0.000 -50.807 -996.615 -1963.162
4 -0.7333333 -0.4653579 -0.000 -59.935 -1312.176 -2598.950
3 -0.8000000 -0.6004851 -0.000 -70.024 -1679.080 -3339.247
2 -0.8666667 -0.7487245 -0.000 -80.770 -2080.962 -4150.744
1 -0.9333333 -0.8938704 -0.000 -91.360 -2474.592 -4945.443
0 -1.0000000 -1.0000000 -0.000 -100.000 -2764.140 -5528.281

Time Splitting Weights for Grid 01: ndtfast = 20 nfast = 29
==================================

Primary Secondary Accumulated to Current Step

1-0.0009651193358779 0.0500000000000000-0.0009651193358779 0.0500000000000000
2-0.0013488780126037 0.0500482559667939-0.0023139973484816 0.1000482559667939
3-0.0011514592651644 0.0501156998674241-0.0034654566136460 0.1501639558342180
4-0.0003735756740661 0.0501732728306823-0.0038390322877122 0.2003372286649003
5 0.0009829200513762 0.0501919516143856-0.0028561122363360 0.2505291802792859
6 0.0029141799764308 0.0501428056118168 0.0000580677400949 0.3006719858911027
7 0.0054132615310267 0.0499970966129953 0.0054713292711215 0.3506690825040980
8 0.0084687837865133 0.0497264335364439 0.0139401130576348 0.4003955160405419
9 0.0120633394191050 0.0493029943471183 0.0260034524767398 0.4496985103876601
10 0.0161716623600090 0.0486998273761630 0.0421751148367488 0.4983983377638231
11 0.0207585511322367 0.0478912442581626 0.0629336659689855 0.5462895820219857
12 0.0257765478740990 0.0468533167015507 0.0887102138430846 0.5931428987235364
13 0.0311633730493853 0.0455644893078458 0.1198735868924699 0.6387073880313822
14 0.0368391158442262 0.0440063206553765 0.1567127027366961 0.6827137086867586
15 0.0427031802506397 0.0421643648631652 0.1994158829873358 0.7248780735499238
16 0.0486309868367616 0.0400292058506332 0.2480468698240974 0.7649072794005570
17 0.0544704302037592 0.0375976565087951 0.3025173000278565 0.8025049359093521
18 0.0600380921294286 0.0348741349986072 0.3625553921572851 0.8373790709079593
19 0.0651152103984763 0.0318722303921357 0.4276706025557614 0.8692513013000950
20 0.0694434033194839 0.0286164698722119 0.4971140058752453 0.8978677711723069
21 0.0727201499285570 0.0251442997062377 0.5698341558038023 0.9230120708785446
22 0.0745940258796570 0.0215082922098099 0.6444281816834592 0.9445203630883545
23 0.0746596950216179 0.0177785909158270 0.7190878767050771 0.9622989540041815
24 0.0724526566618460 0.0140456061647461 0.7915405333669231 0.9763445601689277
25 0.0674437485167025 0.0104229733316538 0.8589842818836255 0.9867675335005816
26 0.0590334053485720 0.0070507859058187 0.9180176872321975 0.9938183194064003
27 0.0465456732896125 0.0040991156383901 0.9645633605218099 0.9979174350447905
28 0.0292219798521905 0.0017718319739095 0.9937853403740003 0.9996892670187000
29 0.0062146596259993 0.0003107329813000 0.9999999999999997 1.0000000000000000

ndtfast, nfast = 20 29 nfast/ndtfast = 1.45000

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

1.000000000000 1.060707743385 0.530353871693 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.14200

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

Minimum X-grid spacing, DXmin = 9.53173644E+00 km
Maximum X-grid spacing, DXmax = 1.06053729E+01 km
Minimum Y-grid spacing, DYmin = 1.11200424E+01 km
Maximum Y-grid spacing, DYmax = 1.11200424E+01 km
Minimum Z-grid spacing, DZmin = 0.00000000E+00 m
Maximum Z-grid spacing, DZmax = 8.11496867E+02 m

Minimum barotropic Courant Number = 5.49001381E-02
Maximum barotropic Courant Number = 9.10581682E-01
Maximum Coriolis Courant Number = 4.45309430E-02

Horizontal mixing scaled by grid size, GRDMAX = 1.08596592E+01 km

NLM: GET_STATE - Read state initial conditions, t = 10 00:00:00
(Grid 01, File: gom_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 = 0.00000000E+00 Max = 3.00944372E+01)
- salinity
(Min = 0.00000000E+00 Max = 3.68434650E+01)
GET_2DFLD - surface u-momentum stress, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -4.46521160E-04 Max = 7.88658003E-05)
GET_2DFLD - surface v-momentum stress, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -2.43290482E-04 Max = 2.63755963E-04)
GET_2DFLD - solar shortwave radiation flux, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = 2.46170107E-05 Max = 9.84193118E-05)
GET_2DFLD - surface net heat flux, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -3.01548208E-05 Max = 4.19633997E-05)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -7.70141743E-03 Max = 2.82108545E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.10545473E-02 Max = 5.74858561E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.16785280E-01 Max = 2.00956576E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.10545473E-02 Max = 7.94643536E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.47733809E-02 Max = 3.51793170E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.24065101E-01 Max = 8.88886154E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.89331788E-01 Max = 5.45626506E-02)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.17090203E-01 Max = 2.46778369E-01)
GET_NGFLD - temperature eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.09180876E-13 Max = 3.00770473E+01)
GET_NGFLD - salinity eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.34068986E-13 Max = 3.65945892E+01)
GET_NGFLD - temperature southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.54044008E-12 Max = 2.94915771E+01)
GET_NGFLD - salinity southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.91542485E-12 Max = 3.66457024E+01)
GET_3DFLD - potential temperature climatology, t = 0 00:00:00
(Rec=0000001, Index=1, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 59.0000)
(Min = 0.00000000E+00 Max = 2.74143713E+01)
GET_3DFLD - salinity climatology, t = 0 00:00:00
(Rec=0000001, Index=1, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 59.0000)
(Min = 0.00000000E+00 Max = 3.67458003E+01)

Basin information for Grid 01:

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

Initial basin volumes: TotVolume = 3.4615820664E+15 m3
MinVolume = 1.1843849025E+08 m3
MaxVolume = 9.5648797369E+10 m3
Max/Min = 8.0758203831E+02

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00001441 - 00008784)

GET_2DFLD - surface u-momentum stress, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -2.81609960E-04 Max = 1.96982535E-05)
GET_2DFLD - surface v-momentum stress, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -8.42378340E-05 Max = 2.03187849E-04)
GET_2DFLD - solar shortwave radiation flux, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -4.22234950E-20 Max = 4.63853085E-06)
GET_2DFLD - surface net heat flux, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -7.62745490E-05 Max = -4.25195064E-06)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.88156872E-03 Max = 2.52899200E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.15718057E-02 Max = 6.22464232E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.28449917E-01 Max = 1.85216349E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.15718057E-02 Max = 9.02452767E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.58911180E-02 Max = 3.12854528E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.16634078E-01 Max = 7.79765025E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.83034706E-01 Max = 3.61946411E-02)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.13600455E-01 Max = 2.92638898E-01)
GET_NGFLD - temperature eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.08041814E-13 Max = 2.99027557E+01)
GET_NGFLD - salinity eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.34010629E-13 Max = 3.65638466E+01)
GET_NGFLD - temperature southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.51450098E-12 Max = 2.96635704E+01)
GET_NGFLD - salinity southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.91184872E-12 Max = 3.66002235E+01)
GET_3DFLD - potential temperature climatology, t = 31 00:00:00
(Rec=0000002, Index=2, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 59.0000)
(Min = 0.00000000E+00 Max = 2.72762837E+01)
GET_3DFLD - salinity climatology, t = 31 00:00:00
(Rec=0000002, Index=2, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 59.0000)
(Min = 0.00000000E+00 Max = 3.67549533E+01)

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

1440 10 00:00:00 1.048932E-05 1.546437E+04 1.546437E+04 3.461582E+15
(090,000,01) 0.000000E+00 0.000000E+00 NaN 2.398735E-01
1446 10 01:00:00 NaN NaN NaN NaN
(090,000,01) NaN NaN NaN 0.000000E+00

Elapsed CPU time (seconds):

Node # 0 CPU: 15.040
Node # 4 CPU: 15.072
Node # 5 CPU: 14.960
Node # 2 CPU: 15.084
Node # 1 CPU: 14.944
Node # 3 CPU: 15.076
Node # 6 CPU: 15.076
Node # 7 CPU: 15.000
Total: 120.252

Nonlinear model elapsed time profile, Grid: 01

Allocation and array initialization .............. 0.072 ( 0.0599 %)
Ocean state initialization ....................... 0.744 ( 0.6187 %)
Reading of input data ............................ 11.528 ( 9.5865 %)
Processing of input data ......................... 4.212 ( 3.5026 %)
Processing of output time averaged data .......... 0.080 ( 0.0665 %)
Computation of vertical boundary conditions ...... 0.392 ( 0.3260 %)
Computation of global information integrals ...... 0.296 ( 0.2461 %)
Model 2D kernel .................................. 75.340 (62.6518 %)
2D/3D coupling, vertical metrics ................. 3.624 ( 3.0137 %)
Omega vertical velocity .......................... 1.108 ( 0.9214 %)
Equation of state for seawater ................... 2.960 ( 2.4615 %)
KPP vertical mixing parameterization ............. 1.992 ( 1.6565 %)
3D equations right-side terms .................... 0.088 ( 0.0732 %)
3D equations predictor step ...................... 0.860 ( 0.7152 %)
Pressure gradient ................................ 0.044 ( 0.0366 %)
Harmonic mixing of tracers, geopotentials ........ 0.080 ( 0.0665 %)
Harmonic stress tensor, S-surfaces ............... 0.052 ( 0.0432 %)
Corrector time-step for 3D momentum .............. 1.376 ( 1.1443 %)
Corrector time-step for tracers .................. 2.732 ( 2.2719 %)
Total: 107.580 89.4621

Nonlinear model message Passage profile, Grid: 01

Message Passage: 2D halo exchanges ............... 82.312 (68.4496 %)
Message Passage: 3D halo exchanges ............... 7.972 ( 6.6294 %)
Message Passage: 4D halo exchanges ............... 3.124 ( 2.5979 %)
Message Passage: data broadcast .................. 7.316 ( 6.0839 %)
Message Passage: data reduction .................. 1.028 ( 0.8549 %)
Message Passage: data scattering.................. 13.524 (11.2464 %)
Total: 115.276 95.8620

All percentages are with respect to total time = 120.252

ROMS/TOMS - Output NetCDF summary for Grid 01:

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_srflux.h
ROMS/Functionals/ana_stflux.h

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

MAIN: Abnormal termination: BLOWUP.

wilkin · #5 Unread post by **wilkin** » Sat Jun 11, 2016 1:08 pm

Why we cannot have zero bathymetry:
https://www.myroms.org/wiki/Frequently_ ... .3D0.29.3F

lalvarez · #6 Unread post by **lalvarez** » Sat Jun 11, 2016 1:30 pm

Hello

I guess you may check your rx0 too, it is a bit greater than 0.4 and as suggested rx0<0.4

Have a nice day

Lee2015 · #7 Unread post by **Lee2015** » Sat Jun 11, 2016 2:02 pm

wilkin and lalvarez,

Thanks for your reply!

Lee2015 · #8 Unread post by **Lee2015** » Sat Jun 11, 2016 3:23 pm

I have modify the zero value and rx0 value, but I still got the blowup error.

Model Input Parameters: ROMS/TOMS version 3.7
Saturday - June 11, 2016 - 11:10:47 PM
-----------------------------------------------------------------------------

Mexico

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -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 : 783M

Local Root : /home/lee/ROMS/src
Header Dir : /home/lee/ROMS/Projects/mexico
Header file : mexico.h
Analytical Dir: /home/lee/ROMS/Projects/mexico

Resolution, Grid 01: 0179x0134x015, Parallel Nodes: 8, Tiling: 002x004

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

8784 ntimes Number of timesteps for 3-D equations.
600.000 dt Timestep size (s) for 3-D equations.
20 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.
F LcycleRST Switch to recycle time-records in restart file.
144000 nRST Number of timesteps between the writing of data
into restart fields.
6 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
144000 nHIS Number of timesteps between the writing fields
into history file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
144000 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
5.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
5.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
1.0000E+02 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.
2.6500E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
1 lmd_Jwt Jerlov water type.
2 Vtransform S-coordinate transformation equation.
2 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.
100.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).
20130801.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.
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.
10.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
-1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
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.
T LtracerCLM(01) Turning ON processing of climatology tracer 01: temp
T LtracerCLM(02) Turning ON 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 Hout(idUsms) Write out surface U-momentum stress.
T Hout(idVsms) Write out surface V-momentum stress.
T Hout(idTsur) Write out surface net heat flux.
T Hout(idTsur) Write out surface net salt flux.
T Hout(idSrad) Write out shortwave radiation flux.
T Hout(idDano) Write out density anomaly.
T Hout(idHsbl) Write out depth of surface boundary layer.

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

Output/Input Files:

Output Restart File: Out/mexico_rst.nc
Output Averages File: Out/mexico_avg.nc
Input Grid File: In/t03031_smth9_grid.nc
Input Nonlinear Initial File: In/gom_ini.nc
Input Forcing File 01: In/gom_shflux_era.nc
Input Forcing File 02: In/gom_sms_era.nc
Input Forcing File 03: In/gom_swflux_era.nc
Input Forcing File 04: In/gom_swrad_era.nc
Input Forcing File 05: In/gom_wind_era.nc
Input Climatology File: In/gom_clm.nc
Input Boundary File: In/gom_bry.nc

Tile partition information for Grid 01: 0179x0134x0015 tiling: 002x004

tile Istr Iend Jstr Jend Npts

0 1 90 1 33 44550
1 91 179 1 33 44055
2 1 90 34 67 45900
3 91 179 34 67 45390
4 1 90 68 101 45900
5 91 179 68 101 45390
6 1 90 102 134 44550
7 91 179 102 134 44055

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 90.50 0.50 33.50 RHO-points
1 90.50 179.50 0.50 33.50 RHO-points
2 0.50 90.50 33.50 67.50 RHO-points
3 90.50 179.50 33.50 67.50 RHO-points
4 0.50 90.50 67.50 101.50 RHO-points
5 90.50 179.50 67.50 101.50 RHO-points
6 0.50 90.50 101.50 134.50 RHO-points
7 90.50 179.50 101.50 134.50 RHO-points

0 1.00 90.50 0.50 33.50 U-points
1 90.50 179.00 0.50 33.50 U-points
2 1.00 90.50 33.50 67.50 U-points
3 90.50 179.00 33.50 67.50 U-points
4 1.00 90.50 67.50 101.50 U-points
5 90.50 179.00 67.50 101.50 U-points
6 1.00 90.50 101.50 134.50 U-points
7 90.50 179.00 101.50 134.50 U-points

0 0.50 90.50 1.00 33.50 V-points
1 90.50 179.50 1.00 33.50 V-points
2 0.50 90.50 33.50 67.50 V-points
3 90.50 179.50 33.50 67.50 V-points
4 0.50 90.50 67.50 101.50 V-points
5 90.50 179.50 67.50 101.50 V-points
6 0.50 90.50 101.50 134.00 V-points
7 90.50 179.50 101.50 134.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 306
HaloSizeJ(1) = 144
TileSide(1) = 96
TileSize(1) = 4032

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

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

zeta 1 Closed Gradient Gradient Closed

ubar 1 Closed Clamped Clamped Closed

vbar 1 Closed Clamped Clamped Closed

u 1 Closed Clamped Clamped Closed

v 1 Closed Clamped Clamped Closed

temp 1 Closed Clamped Clamped Closed

salt 1 Closed Clamped Clamped Closed

Activated C-preprocessing Options:

MEXICO Mexico
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
CURVGRID Orthogonal curvilinear grid.
DIFF_GRID Horizontal diffusion coefficient scaled by grid size.
DIURNAL_SRFLUX Modulate shortwave radiation by the local diurnal cycle.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
LMD_CONVEC LMD convective mixing due to shear instability.
LMD_MIXING Large/McWilliams/Doney interior mixing.
LMD_NONLOCAL LMD convective nonlocal transport.
LMD_RIMIX LMD diffusivity due to shear instability.
LMD_SKPP KPP surface boundary layer mixing.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
PERFECT_RESTART Processing perfect restart variables.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
RI_SPLINES Parabolic Spline Reconstruction for Richardson Number.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
SPLINES_VDIFF Parabolic Spline Reconstruction for Vertical Diffusion.
SPLINES_VVISC Parabolic Spline Reconstruction for Vertical Viscosity.
TS_MPDATA Recursive flux corrected MPDATA 3D 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_QDRAG Quadratic bottom stress.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
VISC_GRID Horizontal viscosity coefficient scaled by grid size.

Process Information:

Node # 0 (pid= 11224) is active.
Node # 4 (pid= 11228) is active.
Node # 1 (pid= 11225) is active.
Node # 5 (pid= 11229) is active.
Node # 6 (pid= 11230) is active.
Node # 7 (pid= 11231) is active.

Node # 2 (pid= 11226) is active.
INITIAL: Configuring and initializing forward nonlinear model ...
Node # 3 (pid= 11227) is active.
*******

Vertical S-coordinate System, Grid 01:

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

15 0.0000000 0.0000000 0.000 0.000 0.000 0.000
14 -0.0666667 -0.0004989 -0.595 -3.358 -7.748 -9.222
13 -0.1333333 -0.0032168 -1.192 -6.828 -21.351 -30.343
12 -0.2000000 -0.0100379 -1.792 -10.502 -45.780 -73.465
11 -0.2666667 -0.0229247 -2.398 -14.480 -86.215 -149.112
10 -0.3333333 -0.0439342 -3.011 -18.863 -148.082 -268.314
9 -0.4000000 -0.0752335 -3.632 -23.762 -237.099 -442.692
8 -0.4666667 -0.1191055 -4.265 -29.289 -359.291 -684.486
7 -0.5333333 -0.1779212 -4.911 -35.563 -520.912 -1006.409
6 -0.6000000 -0.2540303 -5.572 -42.702 -728.164 -1421.063
5 -0.6666667 -0.3494770 -6.250 -50.807 -986.440 -1939.406
4 -0.7333333 -0.4653579 -6.945 -59.935 -1298.633 -2567.321
3 -0.8000000 -0.6004851 -7.658 -70.024 -1661.609 -3298.436
2 -0.8666667 -0.7487245 -8.382 -80.770 -2059.183 -4099.861
1 -0.9333333 -0.8938704 -9.103 -91.360 -2448.594 -4884.698
0 -1.0000000 -1.0000000 -9.790 -100.000 -2735.057 -5460.325

Time Splitting Weights for Grid 01: ndtfast = 20 nfast = 29
==================================

Primary Secondary Accumulated to Current Step

1-0.0009651193358779 0.0500000000000000-0.0009651193358779 0.0500000000000000
2-0.0013488780126037 0.0500482559667939-0.0023139973484816 0.1000482559667939
3-0.0011514592651644 0.0501156998674241-0.0034654566136460 0.1501639558342180
4-0.0003735756740661 0.0501732728306823-0.0038390322877122 0.2003372286649003
5 0.0009829200513762 0.0501919516143856-0.0028561122363360 0.2505291802792859
6 0.0029141799764308 0.0501428056118168 0.0000580677400949 0.3006719858911027
7 0.0054132615310267 0.0499970966129953 0.0054713292711215 0.3506690825040980
8 0.0084687837865133 0.0497264335364439 0.0139401130576348 0.4003955160405419
9 0.0120633394191050 0.0493029943471183 0.0260034524767398 0.4496985103876601
10 0.0161716623600090 0.0486998273761630 0.0421751148367488 0.4983983377638231
11 0.0207585511322367 0.0478912442581626 0.0629336659689855 0.5462895820219857
12 0.0257765478740990 0.0468533167015507 0.0887102138430846 0.5931428987235364
13 0.0311633730493853 0.0455644893078458 0.1198735868924699 0.6387073880313822
14 0.0368391158442262 0.0440063206553765 0.1567127027366961 0.6827137086867586
15 0.0427031802506397 0.0421643648631652 0.1994158829873358 0.7248780735499238
16 0.0486309868367616 0.0400292058506332 0.2480468698240974 0.7649072794005570
17 0.0544704302037592 0.0375976565087951 0.3025173000278565 0.8025049359093521
18 0.0600380921294286 0.0348741349986072 0.3625553921572851 0.8373790709079593
19 0.0651152103984763 0.0318722303921357 0.4276706025557614 0.8692513013000950
20 0.0694434033194839 0.0286164698722119 0.4971140058752453 0.8978677711723069
21 0.0727201499285570 0.0251442997062377 0.5698341558038023 0.9230120708785446
22 0.0745940258796570 0.0215082922098099 0.6444281816834592 0.9445203630883545
23 0.0746596950216179 0.0177785909158270 0.7190878767050771 0.9622989540041815
24 0.0724526566618460 0.0140456061647461 0.7915405333669231 0.9763445601689277
25 0.0674437485167025 0.0104229733316538 0.8589842818836255 0.9867675335005816
26 0.0590334053485720 0.0070507859058187 0.9180176872321975 0.9938183194064003
27 0.0465456732896125 0.0040991156383901 0.9645633605218099 0.9979174350447905
28 0.0292219798521905 0.0017718319739095 0.9937853403740003 0.9996892670187000
29 0.0062146596259993 0.0003107329813000 0.9999999999999997 1.0000000000000000

ndtfast, nfast = 20 29 nfast/ndtfast = 1.45000

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

1.000000000000 1.060707743385 0.530353871693 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.14200

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

Minimum X-grid spacing, DXmin = 9.53173644E+00 km
Maximum X-grid spacing, DXmax = 1.06053729E+01 km
Minimum Y-grid spacing, DYmin = 1.11200424E+01 km
Maximum Y-grid spacing, DYmax = 1.11200424E+01 km
Minimum Z-grid spacing, DZmin = 5.94894030E-01 m
Maximum Z-grid spacing, DZmax = 8.01424881E+02 m

Minimum barotropic Courant Number = 4.34046573E-02
Maximum barotropic Courant Number = 9.04967747E-01
Maximum Coriolis Courant Number = 4.45309430E-02

Horizontal mixing scaled by grid size, GRDMAX = 1.08596592E+01 km

NLM: GET_STATE - Read state initial conditions, t = 10 00:00:00
(Grid 01, File: gom_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 = 0.00000000E+00 Max = 3.00944372E+01)
- salinity
(Min = 0.00000000E+00 Max = 3.68434650E+01)
GET_2DFLD - surface u-momentum stress, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -4.46521160E-04 Max = 7.88658003E-05)
GET_2DFLD - surface v-momentum stress, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -2.43290482E-04 Max = 2.63755963E-04)
GET_2DFLD - solar shortwave radiation flux, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = 2.46170107E-05 Max = 9.84193118E-05)
GET_2DFLD - surface net heat flux, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -3.01548208E-05 Max = 4.19633997E-05)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -8.75419937E-03 Max = 2.78583258E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.94059033E-02 Max = 5.33198379E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.12739071E-01 Max = 1.90938376E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.94059033E-02 Max = 7.94746801E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.47942649E-02 Max = 3.51850867E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.24114767E-01 Max = 8.88591483E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.89335334E-01 Max = 5.40571101E-02)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.17034659E-01 Max = 2.46785000E-01)
GET_NGFLD - temperature eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.17184945E-13 Max = 3.00770473E+01)
GET_NGFLD - salinity eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.33525747E-13 Max = 3.66118050E+01)
GET_NGFLD - temperature southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.54043206E-12 Max = 2.94916019E+01)
GET_NGFLD - salinity southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.91539840E-12 Max = 3.66452904E+01)
GET_3DFLD - potential temperature climatology, t = 0 00:00:00
(Rec=0000001, Index=1, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = 0.00000000E+00 Max = 3.00944372E+01)
GET_3DFLD - salinity climatology, t = 0 00:00:00
(Rec=0000001, Index=1, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = 0.00000000E+00 Max = 3.68434650E+01)

Basin information for Grid 01:

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

Initial basin volumes: TotVolume = 3.4520265032E+15 m3
MinVolume = 7.9353319836E+07 m3
MaxVolume = 9.4461641422E+10 m3
Max/Min = 1.1903930625E+03

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00001441 - 00008784)

GET_2DFLD - surface u-momentum stress, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -2.81609960E-04 Max = 1.96982535E-05)
GET_2DFLD - surface v-momentum stress, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -8.42378340E-05 Max = 2.03187849E-04)
GET_2DFLD - solar shortwave radiation flux, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -4.22234950E-20 Max = 4.63853085E-06)
GET_2DFLD - surface net heat flux, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -7.62745490E-05 Max = -4.25195064E-06)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.82287770E-03 Max = 2.48840362E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.02356186E-02 Max = 5.92395216E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.24096043E-01 Max = 1.82952601E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.02356186E-02 Max = 9.02760997E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.59372259E-02 Max = 3.12866628E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.16848841E-01 Max = 7.79633820E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.83063942E-01 Max = 3.64063419E-02)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.13282099E-01 Max = 2.92650521E-01)
GET_NGFLD - temperature eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.19303584E-13 Max = 2.99027519E+01)
GET_NGFLD - salinity eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.33347667E-13 Max = 3.65748367E+01)
GET_NGFLD - temperature southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.51449003E-12 Max = 2.96636028E+01)
GET_NGFLD - salinity southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.91220173E-12 Max = 3.66001740E+01)
GET_3DFLD - potential temperature climatology, t = 31 00:00:00
(Rec=0000002, Index=2, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = 0.00000000E+00 Max = 2.99086880E+01)
GET_3DFLD - salinity climatology, t = 31 00:00:00
(Rec=0000002, Index=2, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = 0.00000000E+00 Max = 3.66465970E+01)

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

1440 10 00:00:00 1.046594E-05 1.529191E+04 1.529191E+04 3.470239E+15
(179,065,13) 1.001929E-02 0.000000E+00 0.000000E+00 2.399039E-01
1446 10 01:00:00 2.144616E-02 1.528879E+04 1.528881E+04 3.470259E+15
(108,010,15) 2.315671E-02 4.110335E-02 3.154906E+00 3.922275E+00
GET_2DFLD - surface u-momentum stress, t = 10 04:30:00
(Rec=0000082, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -2.76960764E-04 Max = 1.72726149E-05)
GET_2DFLD - surface v-momentum stress, t = 10 04:30:00
(Rec=0000082, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -7.20951979E-05 Max = 1.92146024E-04)
GET_2DFLD - solar shortwave radiation flux, t = 10 04:30:00
(Rec=0000082, Index=1, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -9.37076875E-09 Max = 9.37076875E-09)
GET_2DFLD - surface net heat flux, t = 10 04:30:00
(Rec=0000082, Index=1, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -8.50984428E-05 Max = -2.81554360E-06)
1452 10 02:00:00 7.562188E-02 1.528910E+04 1.528917E+04 3.470278E+15
(045,007,15) 4.073904E-03 1.695849E-02 1.691881E+00 6.039989E+00
1458 10 03:00:00 2.296695E-01 1.528951E+04 1.528974E+04 3.470297E+15
(045,047,08) 1.624805E-02 7.647977E-03 2.633696E+00 1.074933E+01
1464 10 04:00:00 NaN NaN NaN NaN
(090,000,01) NaN NaN NaN 0.000000E+00

Elapsed CPU time (seconds):

Node # 0 CPU: 45.568
Node # 1 CPU: 45.792
Node # 2 CPU: 45.720
Node # 4 CPU: 45.772
Node # 5 CPU: 45.880
Node # 3 CPU: 45.724
Node # 6 CPU: 45.952
Node # 7 CPU: 45.716
Total: 366.124

Nonlinear model elapsed time profile, Grid: 01

Allocation and array initialization .............. 0.096 ( 0.0262 %)
Ocean state initialization ....................... 0.928 ( 0.2535 %)
Reading of input data ............................ 18.336 ( 5.0081 %)
Processing of input data ......................... 12.888 ( 3.5201 %)
Processing of output time averaged data .......... 0.360 ( 0.0983 %)
Computation of vertical boundary conditions ...... 1.236 ( 0.3376 %)
Computation of global information integrals ...... 0.972 ( 0.2655 %)
Model 2D kernel .................................. 263.960 (72.0958 %)
2D/3D coupling, vertical metrics ................. 9.832 ( 2.6854 %)
Omega vertical velocity .......................... 3.504 ( 0.9571 %)
Equation of state for seawater ................... 8.100 ( 2.2124 %)
KPP vertical mixing parameterization ............. 6.488 ( 1.7721 %)
3D equations right-side terms .................... 0.356 ( 0.0972 %)
3D equations predictor step ...................... 3.292 ( 0.8991 %)
Pressure gradient ................................ 0.280 ( 0.0765 %)
Harmonic mixing of tracers, geopotentials ........ 0.336 ( 0.0918 %)
Harmonic stress tensor, S-surfaces ............... 0.160 ( 0.0437 %)
Corrector time-step for 3D momentum .............. 5.592 ( 1.5274 %)
Corrector time-step for tracers .................. 10.840 ( 2.9607 %)
Total: 347.556 94.9285

Nonlinear model message Passage profile, Grid: 01

Message Passage: 2D halo exchanges ............... 282.632 (77.1957 %)
Message Passage: 3D halo exchanges ............... 22.860 ( 6.2438 %)
Message Passage: 4D halo exchanges ............... 10.352 ( 2.8275 %)
Message Passage: data broadcast .................. 13.520 ( 3.6927 %)
Message Passage: data reduction .................. 2.464 ( 0.6730 %)
Message Passage: data scattering.................. 15.644 ( 4.2729 %)
Total: 347.472 94.9056

All percentages are with respect to total time = 366.124

ROMS/TOMS - Output NetCDF summary for Grid 01:

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_srflux.h
ROMS/Functionals/ana_stflux.h

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

MAIN: Abnormal termination: BLOWUP.

susonic · #9 Unread post by **susonic** » Sun Jun 12, 2016 2:47 am

Hi,
You are not allowed to have zero value in your depth.
Check out more information on the below website.

https://myroms.org/forum/viewtopic.php? ... w=previous
viewtopic.php?f=23&t=2546

kate · #10 Unread post by **kate** » Sun Jun 12, 2016 6:34 am

He's gotten beyond the zero depth.

1458 10 03:00:00 2.296695E-01 1.528951E+04 1.528974E+04 3.470297E+15
(045,047,08) 1.624805E-02 7.647977E-03 2.633696E+00 1.074933E+01
1464 10 04:00:00 NaN NaN NaN NaN
(090,000,01) NaN NaN NaN 0.000000E+00

I would set NINFO to 1 so that you can watch it blow up, perhaps get a restart record before you have NaNs everywhere. Then from the wiki FAQ:

ROMS will check for extreme values of velocity and density at run time, including looking for NaN and Inf. If found, it will save a record to the restart file and die. If this happens early in the run, you might just try a shorter timestep. If it happens after many days/years of model time, you need to dig more deeply. Find out which field went bad and where in the domain it happened. The restart record can be viewed to look for this trouble. I've had runs go bad at the bottom with a bottom stress instability, at the top with a surface forcing problem, at the boundary with boundary condition issues. I've also had it go bad in the middle of the water column in no place special due to some other instability - it goes from fine to terrible in just a few timesteps. This last can often be run past with a shorter timestep, while the other problems call for more thoughtful responses.

susonic · #11 Unread post by **susonic** » Sun Jun 12, 2016 12:03 pm

Oh! I saw wrong log file(I should have read third one). Sorry for the wrong post.

Yes, like Kate mentioned, I also recommend you to try shorter time step.

I also have a question about your boundary condition.

As far as I know, the conventional boundary condition for the sea surface and barotropic velocity are

Chapman and Flather, respectively. Would you try that one also?

-JH

wilkin · #12 Unread post by **wilkin** » Sun Jun 12, 2016 2:33 pm

You need to figure out where ROMS is blowing up in order to begin taking rational remedial steps. But you are not saving any output that would help you do that.

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

8784 ntimes Number of timesteps for 3-D equations.
600.000 dt Timestep size (s) for 3-D equations.
20 ndtfast Number of timesteps for 2-D equations between
each 3D timestep.
...
F LcycleRST Switch to recycle time-records in restart file.
144000 nRST Number of timesteps between the writing of data into restart fields.
6 ninfo Number of timesteps between print of information to standard output.
T ldefout Switch to create a new output NetCDF file(s).
144000 nHIS Number of timesteps between the writing fields into history file.
...

Kate has already suggested that you decrease NINFO to 1 instead of 6 so that you can monitor the progression to blow-up. You should do that. The information that is sent to stdout every NINFO steps includes the i,j,k coordinate of the max speed so at the moment of blow-up this gives a clue where to look at your output to see what conditions preceded the blow-up.

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

1440 10 00:00:00 1.046594E-05 1.529191E+04 1.529191E+04 3.470239E+15
(179,065,13) 1.001929E-02 0.000000E+00 0.000000E+00 2.399039E-01
1446 10 01:00:00 2.144616E-02 1.528879E+04 1.528881E+04 3.470259E+15
(108,010,15) 2.315671E-02 4.110335E-02 3.154906E+00 3.922275E+00
^^^^^^^^

But to look at what ROMS had computed so far you need to save some output. You have st nRST, nHIS and nAVG all to 144000, so you won't get any output for that many time steps, which is 1000 days. (It will be hard to do any science with so little output, too).

Until you get these issues sorted out, I suggest you save restart output very frequently, say every nRST = 6 time steps so you can see some evolution of the solution. If you set LcycleRST to TRUE then ROMS will only keep the last 2 records so the file is not going to grow absurdly large in size. At blow-up ROMS will also write the model state to a 3rd record in the restart file. It is not doing this at present because the file has not been created yet. It needs to exist for the blow-up save action to complete.

I suggest you think about your boundary condition choices, and consult with ROMS users in your immediate group. You have everything clamped (except zeta) which is going to generate all manner of wave reflections at the boundaries.

kate · #13 Unread post by **kate** » Sun Jun 12, 2016 3:50 pm

John - this is good stuff. Feel free to pad out the FAQ section of the wiki with these gems.

Is it true that the restart dump won't happen before any records are written? I know you get it for restart runs before new records are written (and it goes into record 1 instead of 3).

Lee2015 · #14 Unread post by **Lee2015** » Sun Jun 12, 2016 4:13 pm

Thank you all！I will check them...

Lee2015 · #15 Unread post by **Lee2015** » Mon Jun 13, 2016 3:00 pm

Based on your suggestions above, I got the new output, include the re-start file (I just post the download link http://www.filedropper.com/mexicorst, since I cannot post attachments here). I want to know what information the re-start file told us? the log file is as follows.

many Thanks!

Model Input Parameters: ROMS/TOMS version 3.7
Monday - June 13, 2016 - 9:26:08 PM
-----------------------------------------------------------------------------

Mexico

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -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 : 783M

Local Root : /home/lee/ROMS/src
Header Dir : /home/lee/ROMS/Projects/mexico
Header file : mexico.h
Analytical Dir: /home/lee/ROMS/Projects/mexico

Resolution, Grid 01: 0179x0134x015, Parallel Nodes: 8, Tiling: 004x002

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

21960 ntimes Number of timesteps for 3-D equations.
240.000 dt Timestep size (s) for 3-D equations.
20 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.
1440 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).
1440 nHIS Number of timesteps between the writing fields
into history file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
1440 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
5.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
5.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
1.0000E+02 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.
2.6500E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
1 lmd_Jwt Jerlov water type.
2 Vtransform S-coordinate transformation equation.
2 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.
100.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).
20130801.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.
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.
10.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
-1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
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.
T LtracerCLM(01) Turning ON processing of climatology tracer 01: temp
T LtracerCLM(02) Turning ON 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 Hout(idUsms) Write out surface U-momentum stress.
T Hout(idVsms) Write out surface V-momentum stress.
T Hout(idTsur) Write out surface net heat flux.
T Hout(idTsur) Write out surface net salt flux.
T Hout(idSrad) Write out shortwave radiation flux.
T Hout(idDano) Write out density anomaly.
T Hout(idHsbl) Write out depth of surface boundary layer.

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

Output/Input Files:

Output Restart File: Out/mexico_rst.nc
Output History File: Out/mexico_his.nc
Output Averages File: Out/mexico_avg.nc
Input Grid File: In/gom_grd.nc
Input Nonlinear Initial File: In/gom_ini.nc
Input Forcing File 01: In/gom_shflux_era.nc
Input Forcing File 02: In/gom_sms_era.nc
Input Forcing File 03: In/gom_swflux_era.nc
Input Forcing File 04: In/gom_swrad_era.nc
Input Forcing File 05: In/gom_wind_era.nc
Input Climatology File: In/gom_clm.nc
Input Boundary File: In/gom_bry.nc

Tile partition information for Grid 01: 0179x0134x0015 tiling: 004x002

tile Istr Iend Jstr Jend Npts

0 1 45 1 67 45225
1 46 90 1 67 45225
2 91 135 1 67 45225
3 136 179 1 67 44220
4 1 45 68 134 45225
5 46 90 68 134 45225
6 91 135 68 134 45225
7 136 179 68 134 44220

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 45.50 0.50 67.50 RHO-points
1 45.50 90.50 0.50 67.50 RHO-points
2 90.50 135.50 0.50 67.50 RHO-points
3 135.50 179.50 0.50 67.50 RHO-points
4 0.50 45.50 67.50 134.50 RHO-points
5 45.50 90.50 67.50 134.50 RHO-points
6 90.50 135.50 67.50 134.50 RHO-points
7 135.50 179.50 67.50 134.50 RHO-points

0 1.00 45.50 0.50 67.50 U-points
1 45.50 90.50 0.50 67.50 U-points
2 90.50 135.50 0.50 67.50 U-points
3 135.50 179.00 0.50 67.50 U-points
4 1.00 45.50 67.50 134.50 U-points
5 45.50 90.50 67.50 134.50 U-points
6 90.50 135.50 67.50 134.50 U-points
7 135.50 179.00 67.50 134.50 U-points

0 0.50 45.50 1.00 67.50 V-points
1 45.50 90.50 1.00 67.50 V-points
2 90.50 135.50 1.00 67.50 V-points
3 135.50 179.50 1.00 67.50 V-points
4 0.50 45.50 67.50 134.00 V-points
5 45.50 90.50 67.50 134.00 V-points
6 90.50 135.50 67.50 134.00 V-points
7 135.50 179.50 67.50 134.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 177
HaloSizeJ(1) = 240
TileSide(1) = 74
TileSize(1) = 3922

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

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

zeta 1 Closed Gradient Gradient Closed

ubar 1 Closed Clamped Clamped Closed

vbar 1 Closed Clamped Clamped Closed

u 1 Closed Clamped Clamped Closed

v 1 Closed Clamped Clamped Closed

temp 1 Closed Clamped Clamped Closed

salt 1 Closed Clamped Clamped Closed

Activated C-preprocessing Options:

MEXICO Mexico
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
CURVGRID Orthogonal curvilinear grid.
DIFF_GRID Horizontal diffusion coefficient scaled by grid size.
DIURNAL_SRFLUX Modulate shortwave radiation by the local diurnal cycle.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
LMD_CONVEC LMD convective mixing due to shear instability.
LMD_MIXING Large/McWilliams/Doney interior mixing.
LMD_NONLOCAL LMD convective nonlocal transport.
LMD_RIMIX LMD diffusivity due to shear instability.
LMD_SKPP KPP surface boundary layer mixing.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
PERFECT_RESTART Processing perfect restart variables.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
RI_SPLINES Parabolic Spline Reconstruction for Richardson Number.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
SPLINES_VDIFF Parabolic Spline Reconstruction for Vertical Diffusion.
SPLINES_VVISC Parabolic Spline Reconstruction for Vertical Viscosity.
TS_MPDATA Recursive flux corrected MPDATA 3D 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_QDRAG Quadratic bottom stress.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
VISC_GRID Horizontal viscosity coefficient scaled by grid size.

Process Information:

Node # 0 (pid= 13856) is active.
Node # 2 (pid= 13858) is active.
Node # 3 (pid= 13859) is active.
Node # 1 (pid= 13857) is active.

INITIAL: Configuring and initializing forward nonlinear model ...
*******
Node # 4 (pid= 13860) is active.
Node # 6 (pid= 13862) is active.
Node # 7 (pid= 13863) is active.
Node # 5 (pid= 13861) is active.

Vertical S-coordinate System, Grid 01:

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

15 0.0000000 0.0000000 0.000 0.000 0.000 0.000
14 -0.0666667 -0.0004989 -0.595 -3.358 -7.748 -9.222
13 -0.1333333 -0.0032168 -1.192 -6.828 -21.351 -30.343
12 -0.2000000 -0.0100379 -1.792 -10.502 -45.780 -73.465
11 -0.2666667 -0.0229247 -2.398 -14.480 -86.215 -149.112
10 -0.3333333 -0.0439342 -3.011 -18.863 -148.082 -268.314
9 -0.4000000 -0.0752335 -3.632 -23.762 -237.099 -442.692
8 -0.4666667 -0.1191055 -4.265 -29.289 -359.291 -684.486
7 -0.5333333 -0.1779212 -4.911 -35.563 -520.912 -1006.409
6 -0.6000000 -0.2540303 -5.572 -42.702 -728.164 -1421.063
5 -0.6666667 -0.3494770 -6.250 -50.807 -986.440 -1939.406
4 -0.7333333 -0.4653579 -6.945 -59.935 -1298.633 -2567.321
3 -0.8000000 -0.6004851 -7.658 -70.024 -1661.609 -3298.436
2 -0.8666667 -0.7487245 -8.382 -80.770 -2059.183 -4099.861
1 -0.9333333 -0.8938704 -9.103 -91.360 -2448.594 -4884.698
0 -1.0000000 -1.0000000 -9.790 -100.000 -2735.057 -5460.325

Time Splitting Weights for Grid 01: ndtfast = 20 nfast = 29
==================================

Primary Secondary Accumulated to Current Step

1-0.0009651193358779 0.0500000000000000-0.0009651193358779 0.0500000000000000
2-0.0013488780126037 0.0500482559667939-0.0023139973484816 0.1000482559667939
3-0.0011514592651644 0.0501156998674241-0.0034654566136460 0.1501639558342180
4-0.0003735756740661 0.0501732728306823-0.0038390322877122 0.2003372286649003
5 0.0009829200513762 0.0501919516143856-0.0028561122363360 0.2505291802792859
6 0.0029141799764308 0.0501428056118168 0.0000580677400949 0.3006719858911027
7 0.0054132615310267 0.0499970966129953 0.0054713292711215 0.3506690825040980
8 0.0084687837865133 0.0497264335364439 0.0139401130576348 0.4003955160405419
9 0.0120633394191050 0.0493029943471183 0.0260034524767398 0.4496985103876601
10 0.0161716623600090 0.0486998273761630 0.0421751148367488 0.4983983377638231
11 0.0207585511322367 0.0478912442581626 0.0629336659689855 0.5462895820219857
12 0.0257765478740990 0.0468533167015507 0.0887102138430846 0.5931428987235364
13 0.0311633730493853 0.0455644893078458 0.1198735868924699 0.6387073880313822
14 0.0368391158442262 0.0440063206553765 0.1567127027366961 0.6827137086867586
15 0.0427031802506397 0.0421643648631652 0.1994158829873358 0.7248780735499238
16 0.0486309868367616 0.0400292058506332 0.2480468698240974 0.7649072794005570
17 0.0544704302037592 0.0375976565087951 0.3025173000278565 0.8025049359093521
18 0.0600380921294286 0.0348741349986072 0.3625553921572851 0.8373790709079593
19 0.0651152103984763 0.0318722303921357 0.4276706025557614 0.8692513013000950
20 0.0694434033194839 0.0286164698722119 0.4971140058752453 0.8978677711723069
21 0.0727201499285570 0.0251442997062377 0.5698341558038023 0.9230120708785446
22 0.0745940258796570 0.0215082922098099 0.6444281816834592 0.9445203630883545
23 0.0746596950216179 0.0177785909158270 0.7190878767050771 0.9622989540041815
24 0.0724526566618460 0.0140456061647461 0.7915405333669231 0.9763445601689277
25 0.0674437485167025 0.0104229733316538 0.8589842818836255 0.9867675335005816
26 0.0590334053485720 0.0070507859058187 0.9180176872321975 0.9938183194064003
27 0.0465456732896125 0.0040991156383901 0.9645633605218099 0.9979174350447905
28 0.0292219798521905 0.0017718319739095 0.9937853403740003 0.9996892670187000
29 0.0062146596259993 0.0003107329813000 0.9999999999999997 1.0000000000000000

ndtfast, nfast = 20 29 nfast/ndtfast = 1.45000

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

1.000000000000 1.060707743385 0.530353871693 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.14200

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

Minimum X-grid spacing, DXmin = 9.53173644E+00 km
Maximum X-grid spacing, DXmax = 1.06053729E+01 km
Minimum Y-grid spacing, DYmin = 1.11200424E+01 km
Maximum Y-grid spacing, DYmax = 1.11200424E+01 km
Minimum Z-grid spacing, DZmin = 5.94894030E-01 m
Maximum Z-grid spacing, DZmax = 8.01424881E+02 m

Minimum barotropic Courant Number = 1.73618629E-02
Maximum barotropic Courant Number = 3.61987099E-01
Maximum Coriolis Courant Number = 1.78123772E-02

Horizontal mixing scaled by grid size, GRDMAX = 1.08596592E+01 km

NLM: GET_STATE - Read state initial conditions, t = 10 00:00:00
(Grid 01, File: gom_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 = 0.00000000E+00 Max = 3.00944372E+01)
- salinity
(Min = 0.00000000E+00 Max = 3.68434650E+01)
GET_2DFLD - surface u-momentum stress, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -4.44913672E-04 Max = 7.82397046E-05)
GET_2DFLD - surface v-momentum stress, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -2.40179656E-04 Max = 2.63159577E-04)
GET_2DFLD - solar shortwave radiation flux, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = 2.48795581E-05 Max = 9.84199738E-05)
GET_2DFLD - surface net heat flux, t = 9 22:30:00
(Rec=0000080, Index=1, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -3.01287096E-05 Max = 4.17307627E-05)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -8.51026922E-03 Max = 2.78584391E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.92203652E-02 Max = 5.33318147E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.12915687E-01 Max = 1.86779350E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.92203652E-02 Max = 7.93284103E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.46363117E-02 Max = 3.48533183E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.23634353E-01 Max = 8.88170972E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.88949215E-01 Max = 5.43890595E-02)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.17089726E-01 Max = 2.46610597E-01)
GET_NGFLD - temperature eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.14433153E-13 Max = 3.00768490E+01)
GET_NGFLD - salinity eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.35830842E-13 Max = 3.66052551E+01)
GET_NGFLD - temperature southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.54100701E-12 Max = 2.94919109E+01)
GET_NGFLD - salinity southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.90932795E-12 Max = 3.66514893E+01)
GET_3DFLD - potential temperature climatology, t = 0 00:00:00
(Rec=0000001, Index=1, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = 0.00000000E+00 Max = 3.00944372E+01)
GET_3DFLD - salinity climatology, t = 0 00:00:00
(Rec=0000001, Index=1, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = 0.00000000E+00 Max = 3.68434650E+01)

Basin information for Grid 01:

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

Initial basin volumes: TotVolume = 3.4520265032E+15 m3
MinVolume = 7.9353319836E+07 m3
MaxVolume = 9.4461641422E+10 m3
Max/Min = 1.1903930625E+03

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00003601 - 00021960)

GET_2DFLD - surface u-momentum stress, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -2.77185324E-04 Max = 1.97672739E-05)
GET_2DFLD - surface v-momentum stress, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -8.34200513E-05 Max = 2.03272235E-04)
GET_2DFLD - solar shortwave radiation flux, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -4.22234950E-20 Max = 4.63856167E-06)
GET_2DFLD - surface net heat flux, t = 10 01:30:00
(Rec=0000081, Index=2, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -7.59643827E-05 Max = -4.42213781E-06)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.69327252E-03 Max = 2.47487098E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.97281577E-02 Max = 5.92264719E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.23686261E-01 Max = 1.74726173E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.97281577E-02 Max = 9.00412723E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.61137376E-02 Max = 3.11915249E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.14210762E-01 Max = 7.79447556E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -3.80391806E-01 Max = 3.62655930E-02)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.13598138E-01 Max = 2.92345315E-01)
GET_NGFLD - temperature eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -2.13158240E-13 Max = 2.99025478E+01)
GET_NGFLD - salinity eastern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -4.35715293E-13 Max = 3.65579834E+01)
GET_NGFLD - temperature southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.51527662E-12 Max = 2.96640129E+01)
GET_NGFLD - salinity southern boundary condition, t = 31 00:00:00
(Rec=0000002, Index=1, File: gom_bry.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = -1.89893479E-12 Max = 3.66007004E+01)
GET_3DFLD - potential temperature climatology, t = 31 00:00:00
(Rec=0000002, Index=2, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = 0.00000000E+00 Max = 2.99086880E+01)
GET_3DFLD - salinity climatology, t = 31 00:00:00
(Rec=0000002, Index=2, File: gom_clm.nc)
(Tmin= 0.0000 Tmax= 61.0000)
(Min = 0.00000000E+00 Max = 3.66465970E+01)

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

3600 10 00:00:00 3.984212E-06 1.529191E+04 1.529191E+04 3.470239E+15
(179,065,14) 3.968520E-03 0.000000E+00 0.000000E+00 2.375576E-01
DEF_HIS - creating history file, Grid 01: Out/mexico_his.nc
DEF_AVG - creating average file, Grid 01: Out/mexico_avg.nc
3601 10 00:04:00 2.331361E-04 1.529169E+04 1.529169E+04 3.470241E+15
(179,065,10) 4.384972E-03 1.131634E-04 1.365077E-02 4.166317E-01
3602 10 00:08:00 7.249875E-04 1.529124E+04 1.529124E+04 3.470242E+15
(122,029,15) 2.481900E-03 3.597630E-03 3.793006E-01 7.185987E-01
3603 10 00:12:00 1.362833E-03 1.529075E+04 1.529075E+04 3.470244E+15
(110,003,15) 5.058651E-03 3.855910E-04 4.159279E-01 9.214140E-01
3604 10 00:16:00 2.136378E-03 1.529031E+04 1.529031E+04 3.470245E+15
(109,004,15) 2.499647E-03 2.152558E-04 6.768207E-01 1.130340E+00
3605 10 00:20:00 3.076398E-03 1.528993E+04 1.528993E+04 3.470246E+15
(108,001,15) 9.948771E-03 0.000000E+00 7.567061E-01 1.309797E+00
3606 10 00:24:00 4.194239E-03 1.528961E+04 1.528961E+04 3.470248E+15
(108,001,15) 6.484163E-03 0.000000E+00 9.769017E-01 1.597078E+00
3607 10 00:28:00 5.478975E-03 1.528934E+04 1.528935E+04 3.470249E+15
(019,032,15) 5.986766E-03 5.037206E-03 8.355048E-01 1.941845E+00
3608 10 00:32:00 6.931891E-03 1.528912E+04 1.528913E+04 3.470251E+15
(109,009,15) 5.106881E-02 1.700630E-02 7.999557E-01 2.559568E+00
3609 10 00:36:00 8.561587E-03 1.528893E+04 1.528894E+04 3.470252E+15
(108,009,15) 2.974816E-02 8.403897E-03 8.575925E-01 3.117455E+00
3610 10 00:40:00 1.035894E-02 1.528878E+04 1.528879E+04 3.470254E+15
(108,009,15) 3.498620E-02 1.543519E-02 1.105194E+00 3.367350E+00
3611 10 00:44:00 1.229890E-02 1.528868E+04 1.528870E+04 3.470255E+15
(108,009,15) 3.297677E-02 2.318000E-02 1.067438E+00 3.188883E+00
3612 10 00:48:00 1.436901E-02 1.528864E+04 1.528866E+04 3.470256E+15
(108,010,15) 1.774792E-02 2.876980E-02 1.306272E+00 3.313093E+00
3613 10 00:52:00 1.657581E-02 1.528866E+04 1.528868E+04 3.470258E+15
(108,010,15) 1.313418E-02 2.987557E-02 1.458393E+00 3.537982E+00
3614 10 00:56:00 1.892239E-02 1.528872E+04 1.528874E+04 3.470259E+15
(108,010,15) 3.300062E-03 2.594292E-02 1.309494E+00 3.749671E+00
3615 10 01:00:00 2.140175E-02 1.528881E+04 1.528883E+04 3.470261E+15
(108,010,15) 8.966180E-03 1.824740E-02 8.407136E-01 3.936472E+00
3616 10 01:04:00 2.400817E-02 1.528891E+04 1.528893E+04 3.470262E+15
(109,012,15) 0.000000E+00 2.740508E-02 7.131155E-01 4.097715E+00
3617 10 01:08:00 2.675458E-02 1.528898E+04 1.528901E+04 3.470263E+15
(109,012,15) 0.000000E+00 2.966750E-02 8.109075E-01 4.254162E+00
3618 10 01:12:00 2.963748E-02 1.528901E+04 1.528904E+04 3.470265E+15
(108,010,15) 2.505680E-02 2.420295E-03 1.065913E+00 4.427669E+00
3619 10 01:16:00 3.265485E-02 1.528900E+04 1.528903E+04 3.470266E+15
(108,010,15) 1.955225E-02 2.275597E-03 1.132337E+00 4.622279E+00
3620 10 01:20:00 3.579026E-02 1.528895E+04 1.528899E+04 3.470268E+15
(108,010,15) 1.218631E-02 1.108611E-03 8.360314E-01 4.822438E+00
3621 10 01:24:00 3.905354E-02 1.528890E+04 1.528894E+04 3.470269E+15
(109,012,15) 0.000000E+00 2.618865E-02 7.640983E-01 5.014980E+00
3622 10 01:28:00 4.245375E-02 1.528886E+04 1.528890E+04 3.470270E+15
(109,010,15) 3.573767E-03 2.770917E-02 9.416745E-01 5.194289E+00
GET_2DFLD - surface u-momentum stress, t = 10 04:30:00
(Rec=0000082, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -2.73595886E-04 Max = 1.66984640E-05)
GET_2DFLD - surface v-momentum stress, t = 10 04:30:00
(Rec=0000082, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -7.22489153E-05 Max = 1.92062339E-04)
GET_2DFLD - solar shortwave radiation flux, t = 10 04:30:00
(Rec=0000082, Index=1, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -9.37083168E-09 Max = 9.37083168E-09)
GET_2DFLD - surface net heat flux, t = 10 04:30:00
(Rec=0000082, Index=1, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 91.9375)
(Min = -8.47400282E-05 Max = -3.07336624E-06)
3623 10 01:32:00 4.599358E-02 1.528884E+04 1.528889E+04 3.470272E+15
(109,010,15) 4.617939E-03 2.224240E-02 8.550595E-01 5.355370E+00
3624 10 01:36:00 4.968190E-02 1.528884E+04 1.528889E+04 3.470273E+15
(121,040,15) 9.589057E-03 5.186521E-03 8.287023E-01 5.501056E+00
3625 10 01:40:00 5.352390E-02 1.528885E+04 1.528890E+04 3.470275E+15
(044,008,15) 3.196751E-03 1.062091E-02 7.694056E-01 5.619483E+00
3626 10 01:44:00 5.753623E-02 1.528887E+04 1.528893E+04 3.470276E+15
(044,008,15) 2.394450E-03 1.030454E-02 7.647498E-01 5.714415E+00
3627 10 01:48:00 6.172572E-02 1.528891E+04 1.528897E+04 3.470277E+15
(109,010,15) 8.085987E-03 5.163271E-03 8.453291E-01 5.801720E+00
3628 10 01:52:00 6.609474E-02 1.528897E+04 1.528904E+04 3.470279E+15
(049,008,15) 2.704955E-03 1.112156E-02 7.692953E-01 5.885102E+00
3629 10 01:56:00 7.066218E-02 1.528905E+04 1.528912E+04 3.470280E+15
(121,040,15) 1.101593E-02 3.912608E-03 7.780297E-01 5.964541E+00
3630 10 02:00:00 7.547039E-02 1.528913E+04 1.528921E+04 3.470282E+15
(121,040,15) 8.869402E-03 9.083749E-04 7.301998E-01 6.149596E+00
3631 10 02:04:00 8.055739E-02 1.528922E+04 1.528930E+04 3.470283E+15
(138,121,15) 7.048984E-03 8.731135E-03 6.906800E-01 6.327256E+00
3632 10 02:08:00 8.595955E-02 1.528929E+04 1.528937E+04 3.470284E+15
(136,123,15) 1.027427E-02 1.786717E-02 9.204820E-01 6.488789E+00
3633 10 02:12:00 9.172380E-02 1.528933E+04 1.528942E+04 3.470286E+15
(136,123,15) 6.703218E-03 1.772417E-02 9.423716E-01 6.634722E+00
3634 10 02:16:00 9.791744E-02 1.528935E+04 1.528945E+04 3.470287E+15
(140,121,15) 2.649422E-02 9.133879E-03 7.681593E-01 6.766490E+00
3635 10 02:20:00 1.046267E-01 1.528934E+04 1.528944E+04 3.470289E+15
(140,121,15) 2.162363E-02 7.775663E-03 8.224076E-01 6.886544E+00
3636 10 02:24:00 1.119589E-01 1.528931E+04 1.528943E+04 3.470290E+15
(109,011,15) 0.000000E+00 2.848855E-02 8.718303E-01 6.999169E+00
3637 10 02:28:00 1.200489E-01 1.528929E+04 1.528941E+04 3.470291E+15
(121,040,15) 1.426204E-02 2.905948E-03 9.133727E-01 7.110603E+00
3638 10 02:32:00 1.290457E-01 1.528927E+04 1.528940E+04 3.470293E+15
(121,040,15) 1.114994E-02 1.637831E-04 8.640343E-01 7.202179E+00
3639 10 02:36:00 1.391186E-01 1.528928E+04 1.528942E+04 3.470294E+15
(047,007,15) 6.832116E-03 1.572893E-02 7.884898E-01 7.265214E+00
3640 10 02:40:00 1.504336E-01 1.528931E+04 1.528946E+04 3.470296E+15
(047,007,15) 7.002819E-03 1.662025E-02 9.429529E-01 7.295220E+00
3641 10 02:44:00 1.631245E-01 1.528936E+04 1.528952E+04 3.470297E+15
(047,007,15) 6.181791E-03 1.535012E-02 9.674423E-01 7.445450E+00
3642 10 02:48:00 1.773106E-01 1.528942E+04 1.528959E+04 3.470298E+15
(121,040,15) 2.560130E-02 3.437448E-03 1.139805E+00 8.617138E+00
3643 10 02:52:00 1.931022E-01 1.528947E+04 1.528966E+04 3.470300E+15
(121,040,15) 2.663875E-02 1.444769E-03 1.395839E+00 9.867375E+00
3644 10 02:56:00 2.107122E-01 1.528952E+04 1.528973E+04 3.470301E+15
(121,040,15) 2.148073E-02 1.770466E-03 1.332578E+00 1.119750E+01
3645 10 03:00:00 2.303688E-01 1.528956E+04 1.528979E+04 3.470303E+15
(045,040,09) 6.122850E-03 2.276929E-03 1.237030E+00 1.257177E+01
3646 10 03:04:00 2.523248E-01 1.528959E+04 1.528984E+04 3.470304E+15
(045,040,10) 4.889920E-03 3.377145E-03 1.284976E+00 1.394920E+01
3647 10 03:08:00 2.769665E-01 1.528962E+04 1.528990E+04 3.470305E+15
(045,040,10) 5.193763E-03 4.910245E-03 1.310795E+00 1.521999E+01
3648 10 03:12:00 3.048203E-01 1.528964E+04 1.528994E+04 3.470307E+15
(045,040,11) 3.668202E-03 6.175925E-03 1.300097E+00 1.622125E+01
3649 10 03:16:00 3.364472E-01 1.528965E+04 1.528998E+04 3.470308E+15
(121,040,15) 1.733358E-02 1.677701E-02 1.431345E+00 1.671645E+01
3650 10 03:20:00 3.724959E-01 1.528964E+04 1.529001E+04 3.470310E+15
(121,040,15) 1.845981E-02 1.336600E-02 1.611351E+00 1.643193E+01
3651 10 03:24:00 4.135850E-01 1.528961E+04 1.529002E+04 3.470311E+15
(121,040,15) 1.678576E-02 7.655297E-03 1.511983E+00 1.549626E+01
3652 10 03:28:00 4.603240E-01 1.528956E+04 1.529002E+04 3.470313E+15
(014,074,15) 1.539766E-02 2.014759E-02 1.333331E+00 1.462018E+01
3653 10 03:32:00 5.132328E-01 1.528949E+04 1.529001E+04 3.470314E+15
(014,074,15) 1.244568E-02 1.654827E-02 1.264256E+00 1.411822E+01
3654 10 03:36:00 5.724773E-01 1.528941E+04 1.528999E+04 3.470315E+15
(047,042,10) 5.101470E-03 7.731331E-03 1.211294E+00 1.338915E+01
3655 10 03:40:00 6.379735E-01 1.528933E+04 1.528997E+04 3.470317E+15
(047,042,10) 5.593948E-03 7.345588E-03 1.172573E+00 1.406011E+01
3656 10 03:44:00 7.092232E-01 1.528926E+04 1.528997E+04 3.470318E+15
(121,040,15) 2.070829E-02 8.687044E-03 1.232906E+00 1.444265E+01
3657 10 03:48:00 7.853810E-01 1.528920E+04 1.528999E+04 3.470320E+15
(121,040,15) 2.144827E-02 4.054632E-03 1.379919E+00 1.445850E+01
3658 10 03:52:00 8.654652E-01 1.528917E+04 1.529003E+04 3.470321E+15
(121,040,15) 1.732117E-02 1.218217E-03 1.213976E+00 1.405132E+01
3659 10 03:56:00 9.485126E-01 1.528913E+04 1.529008E+04 3.470322E+15
(125,101,08) 2.542501E-03 6.378816E-03 1.125459E+00 1.306192E+01
3660 10 04:00:00 1.033557E+00 1.528909E+04 1.529013E+04 3.470324E+15
(125,101,09) 2.423781E-03 9.923938E-03 1.165756E+00 1.368127E+01
3661 10 04:04:00 1.119441E+00 1.528904E+04 1.529016E+04 3.470325E+15
(125,102,08) 3.118956E-03 4.393947E-02 1.169328E+00 1.411867E+01
3662 10 04:08:00 1.205361E+00 1.528897E+04 1.529018E+04 3.470327E+15
(125,102,09) 2.737442E-03 7.634651E-02 1.229841E+00 1.429991E+01
3663 10 04:12:00 1.290547E+00 1.528889E+04 1.529018E+04 3.470328E+15
(125,102,10) 1.397418E-03 1.338453E-01 1.256176E+00 1.427759E+01
3664 10 04:16:00 1.375694E+00 1.528880E+04 1.529018E+04 3.470329E+15
(125,102,11) 1.189798E-03 2.338407E-01 1.248262E+00 1.479935E+01
3665 10 04:20:00 1.467634E+00 1.528871E+04 1.529018E+04 3.470331E+15
(125,101,10) 2.241018E-03 4.462208E-01 1.394132E+00 3.398482E+01

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: 112.768
Node # 4 CPU: 113.556
Node # 6 CPU: 113.372
Node # 7 CPU: 113.460
Node # 1 CPU: 112.928
Node # 2 CPU: 112.872
Node # 3 CPU: 111.844
Node # 5 CPU: 111.656
Total: 902.456

Nonlinear model elapsed time profile, Grid: 01

Allocation and array initialization .............. 0.080 ( 0.0089 %)
Ocean state initialization ....................... 1.076 ( 0.1192 %)
Reading of input data ............................ 12.140 ( 1.3452 %)
Processing of input data ......................... 31.064 ( 3.4422 %)
Processing of output time averaged data .......... 1.112 ( 0.1232 %)
Computation of vertical boundary conditions ...... 3.408 ( 0.3776 %)
Computation of global information integrals ...... 11.672 ( 1.2934 %)
Writing of output data ........................... 5.756 ( 0.6378 %)
Model 2D kernel .................................. 709.952 (78.6689 %)
2D/3D coupling, vertical metrics ................. 17.968 ( 1.9910 %)
Omega vertical velocity .......................... 7.712 ( 0.8546 %)
Equation of state for seawater ................... 13.956 ( 1.5464 %)
KPP vertical mixing parameterization ............. 14.924 ( 1.6537 %)
3D equations right-side terms .................... 1.308 ( 0.1449 %)
3D equations predictor step ...................... 6.840 ( 0.7579 %)
Pressure gradient ................................ 0.984 ( 0.1090 %)
Harmonic mixing of tracers, geopotentials ........ 1.040 ( 0.1152 %)
Harmonic stress tensor, S-surfaces ............... 0.476 ( 0.0527 %)
Corrector time-step for 3D momentum .............. 14.020 ( 1.5535 %)
Corrector time-step for tracers .................. 22.048 ( 2.4431 %)
Total: 877.536 97.2386

Nonlinear model message Passage profile, Grid: 01

Message Passage: 2D halo exchanges ............... 739.876 (81.9847 %)
Message Passage: 3D halo exchanges ............... 45.656 ( 5.0591 %)
Message Passage: 4D halo exchanges ............... 13.224 ( 1.4653 %)
Message Passage: data broadcast .................. 15.076 ( 1.6706 %)
Message Passage: data reduction .................. 12.036 ( 1.3337 %)
Message Passage: data gathering .................. 3.468 ( 0.3843 %)
Message Passage: data scattering.................. 13.456 ( 1.4910 %)
Total: 842.792 93.3887

All percentages are with respect to total time = 902.456

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

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_srflux.h
ROMS/Functionals/ana_stflux.h

ROMS/TOMS: DONE... Monday - June 13, 2016 - 9:29:56 PM

kate · #16 Unread post by **kate** » Mon Jun 13, 2016 3:16 pm

I can't speak for the rest, but I don't want to look at your restart file. You need to start developing tools for looking at ROMS output, including restart files. The simplest way to start is with ncview, which will show you the fields in any netcdf file (or group of files). That is often good enough for this sort of blow-up.

Beyond that, you need to pick a language, whether Matlab, Ferret, Python, Fortran, NCL, whatever. It has to have a means of reading netcdf files and of creating plots. There are ROMS tools for all of the above, some easier to find than others. Anyway, pick one and start using it to see what the model output looks like.

Lee2015 · #17 Unread post by **Lee2015** » Tue Jun 14, 2016 1:50 am

kate wrote:I can't speak for the rest, but I don't want to look at your restart file. You need to start developing tools for looking at ROMS output, including restart files. The simplest way to start is with ncview, which will show you the fields in any netcdf file (or group of files). That is often good enough for this sort of blow-up.

Beyond that, you need to pick a language, whether Matlab, Ferret, Python, Fortran, NCL, whatever. It has to have a means of reading netcdf files and of creating plots. There are ROMS tools for all of the above, some easier to find than others. Anyway, pick one and start using it to see what the model output looks like.

kate, thank you! I have the tools to read netcdf, I can open the restart file. But I am a ROMS newbie, I don't know how to check the reason of the blow-up. Can you give me some suggestions about it?

lindandan · Thu Jun 16, 2016 9:09 am

kate wrote:I can't speak for the rest, but I don't want to look at your restart file. You need to start developing tools for looking at ROMS output, including restart files. The simplest way to start is with ncview, which will show you the fields in any netcdf file (or group of files). That is often good enough for this sort of blow-up.

Beyond that, you need to pick a language, whether Matlab, Ferret, Python, Fortran, NCL, whatever. It has to have a means of reading netcdf files and of creating plots. There are ROMS tools for all of the above, some easier to find than others. Anyway, pick one and start using it to see what the model output looks like.

Hello! I have a problem when run a spinup with climatology data.
I have run 855 days without error but blowup suddenly. I see the restart file,it has 11 variables 'zeta' 'u' 'v' 'ubar' 'vbar' 'temp' 'salt' 'rho' 'akv' 'akt' 'aks',and I found that there is a abnormal value at right boundary with coordinate(245,276). so I correct the value with its near value. But I just have corrected 9 variables without 'Akt''Akv''Aks',because I don't know how to correct them. And after I correct 9 variable of the restart file. I try to run with restart mode, but failed. What should I do ?
-----------------------
Next is the output
Model Input Parameters: ROMS/TOMS version 3.7
Sunday - June 12, 2016 - 12:33:04 PM
-----------------------------------------------------------------------------
Region Ocena System, 1/12 degree resolution

Operating system : Linux
CPU/hardware : x86_64
Compiler system : ifort
Compiler command : /usr/local/mpi/bin/mpif90
Compiler flags : -heap-arrays -fp-model precise -I/usr/local/mpi/include -ip -O3 -free -free -free

Input Script :
SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision : 783M

Local Root : /vol/home/student/lshw/roms/trunk
Header Dir : /vol/home/student/lshw/roms/MYPROJECT/SSC12
Header file : ssc12.h
Analytical Dir: /vol/home/student/lshw/roms/MYPROJECT/Functionals

Resolution, Grid 01: 0275x0275x032, Parallel Nodes: 16, Tiling: 004x004

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

424800 ntimes Number of timesteps for 3-D equations.
360.000 dt Timestep size (s) for 3-D equations.
72 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.
1 nrrec Number of restart records to read from disk.
T LcycleRST Switch to recycle time-records in restart file.
720 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
F ldefout Switch to create a new output NetCDF file(s).
1200 nHIS Number of timesteps between the writing fields
into history file.
14400 ndefHIS Number of timesteps between creation of new
history files.
8640000 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
864000 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
864000 ndefAVG Number of timesteps between creation of new
time-averaged file.
2.5000E+01 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
2.5000E+01 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
5.0000E+01 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
T LuvSponge Turning ON sponge on horizontal momentum.
T LtracerSponge(01) Turning ON sponge on tracer 01: temp
T LtracerSponge(02) Turning ON sponge on tracer 02: salt
1.0000E-05 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
1.0000E-05 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
1.0000E-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
5.0000E-06 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
5.0000E-06 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
2.5000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
6.0000E+00 theta_s S-coordinate surface control parameter.
5.0000E-01 theta_b S-coordinate bottom control parameter.
15.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).
-1.00 time_ref Reference time for units attribute (yyyymmdd.dd)
3.6000E+02 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
3.6000E+02 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
3.6000E+02 Znudg Nudging/relaxation time scale (days)
for free-surface.
3.6000E+02 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
3.6000E+02 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
1.2000E+02 obcfac Factor between passive and active
open boundary conditions.
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.
5.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
-1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
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.
T LtracerCLM(01) Turning ON processing of climatology tracer 01: temp
T LtracerCLM(02) Turning ON processing of climatology tracer 02: salt
F LnudgeM2CLM Turning OFF nudging of 2D momentum climatology.
F LnudgeM3CLM Turning OFF nudging of 3D momentum climatology.
T LnudgeTCLM(01) Turning ON nudging of climatology tracer 01: temp
T LnudgeTCLM(02) Turning ON 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 Hout(idUsms) Write out surface U-momentum stress.
T Hout(idVsms) Write out surface V-momentum stress.
T Hout(idTsur) Write out surface net heat flux.
T Hout(idTsur) Write out surface net salt flux.
T Hout(idSrad) Write out shortwave radiation flux.
T Hout(idDano) Write out density anomaly.
T Hout(idVvis) Write out vertical viscosity: AKv.
T Hout(idTdif) Write out vertical diffusion: AKt(itemp).
T Hout(idSdif) Write out vertical diffusion: AKt(isalt).
T Hout(idMtke) Write out turbulent kinetic energy.
T Hout(idMtls) Write out turbulent generic length-scale.

Output/Input Files:

Output Restart File: ssc12_rst.nc
Prefix for History Files: ssc12_his
Node # 9 (pid= 10133) is active.
Prefix for Averages Files: ssc12_avg
Node # 12 (pid= 10136) is active.
Input Grid File: ../Data/ssc12_grd.nc
Node # 13 (pid= 10137) is active.
Input Nonlinear Initial File: ssc12_rst.nc
Node # 8 (pid= 10132) is active.
Input Forcing File 01: ../Data/ssc12_frc.nc
Node # 10 (pid= 10134) is active.
Input Climatology File: ../Data/ssc12_clm.nc
Node # 11 (pid= 10135) is active.
Input Boundary File: ../Data/ssc12_bry.nc
Node # 14 (pid= 10138) is active.
Input/Output USER File: MyFile.dat
Node # 15 (pid= 10139) is active.

Generic User Parameters:

2.0000E+01 user(01) User parameter 01.

Tile partition information for Grid 01: 0275x0275x0032 tiling: 004x004

tile Istr Iend Jstr Jend Npts

0 1 69 1 69 152352
1 70 138 1 69 152352
2 139 207 1 69 152352
3 208 275 1 69 150144
4 1 69 70 138 152352
5 70 138 70 138 152352
6 139 207 70 138 152352
7 208 275 70 138 150144
8 1 69 139 207 152352
9 70 138 139 207 152352
10 139 207 139 207 152352
11 208 275 139 207 150144
12 1 69 208 275 150144
13 70 138 208 275 150144
14 139 207 208 275 150144
15 208 275 208 275 147968

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 69.50 0.50 69.50 RHO-points
1 69.50 138.50 0.50 69.50 RHO-points
2 138.50 207.50 0.50 69.50 RHO-points
3 207.50 275.50 0.50 69.50 RHO-points
4 0.50 69.50 69.50 138.50 RHO-points
5 69.50 138.50 69.50 138.50 RHO-points
6 138.50 207.50 69.50 138.50 RHO-points
7 207.50 275.50 69.50 138.50 RHO-points
8 0.50 69.50 138.50 207.50 RHO-points
9 69.50 138.50 138.50 207.50 RHO-points
10 138.50 207.50 138.50 207.50 RHO-points
11 207.50 275.50 138.50 207.50 RHO-points
12 0.50 69.50 207.50 275.50 RHO-points
13 69.50 138.50 207.50 275.50 RHO-points
14 138.50 207.50 207.50 275.50 RHO-points
15 207.50 275.50 207.50 275.50 RHO-points

0 1.00 69.50 0.50 69.50 U-points
1 69.50 138.50 0.50 69.50 U-points
2 138.50 207.50 0.50 69.50 U-points
3 207.50 275.00 0.50 69.50 U-points
4 1.00 69.50 69.50 138.50 U-points
5 69.50 138.50 69.50 138.50 U-points
6 138.50 207.50 69.50 138.50 U-points
7 207.50 275.00 69.50 138.50 U-points
8 1.00 69.50 138.50 207.50 U-points
9 69.50 138.50 138.50 207.50 U-points
10 138.50 207.50 138.50 207.50 U-points
11 207.50 275.00 138.50 207.50 U-points
12 1.00 69.50 207.50 275.50 U-points
13 69.50 138.50 207.50 275.50 U-points
14 138.50 207.50 207.50 275.50 U-points
15 207.50 275.00 207.50 275.50 U-points

0 0.50 69.50 1.00 69.50 V-points
1 69.50 138.50 1.00 69.50 V-points
2 138.50 207.50 1.00 69.50 V-points
3 207.50 275.50 1.00 69.50 V-points
4 0.50 69.50 69.50 138.50 V-points
5 69.50 138.50 69.50 138.50 V-points
6 138.50 207.50 69.50 138.50 V-points
7 207.50 275.50 69.50 138.50 V-points
8 0.50 69.50 138.50 207.50 V-points
9 69.50 138.50 138.50 207.50 V-points
10 138.50 207.50 138.50 207.50 V-points
11 207.50 275.50 138.50 207.50 V-points
12 0.50 69.50 207.50 275.00 V-points
13 69.50 138.50 207.50 275.00 V-points
14 138.50 207.50 207.50 275.00 V-points
15 207.50 275.50 207.50 275.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 162
HaloSizeJ(1) = 162
TileSide(1) = 75
TileSize(1) = 5625

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

v 1 Radiation Radiation Radiation Radiation

temp 1 Clamped Clamped Clamped Clamped

salt 1 Clamped Clamped Clamped Clamped

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

Activated C-preprocessing Options:

SSC12 South Sea of China Internal tidal System, 1/12 degree resolution
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_NUDGCOEF Analytical spatially varying nudging time-scales.
ANA_SPONGE Analytical enhanced viscosity/diffusion sponge.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
CURVGRID Orthogonal curvilinear grid.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_GEO_UV Mixing of momentum along geopotential surfaces.
MPI MPI distributed-memory configuration.
MY25_MIXING Mellor/Yamada Level-2.5 mixing closure.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
QCORRECTION Surface net heat flux correction.
K_GSCHEME Third-order upstream advection of TKE fields.
RADIATION_2D Use tangential phase speed in radiation conditions.
RI_SPLINES Parabolic Spline Reconstruction for Richardson Number.
!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_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_SVADVECTION Parabolic splines 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_QDRAG Quadratic bottom stress.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.

Process Information:

Node # 0 (pid= 10246) is active.
Node # 1 (pid= 10247) is active.
Node # 5 (pid= 10251) is active.
Node # 3 (pid= 10249) is active.
Node # 7 (pid= 10253) is active.
Node # 2 (pid= 10248) is active.
Node # 4 (pid= 10250) is active.
Node # 6 (pid= 10252) 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

32 0.0000000 0.0000000 0.000 0.000 0.000 0.000
31 -0.0312500 -0.0001116 -0.188 -0.235 -0.800 -1.136
30 -0.0625000 -0.0004504 -0.377 -0.472 -2.279 -3.630
29 -0.0937500 -0.0010281 -0.567 -0.711 -4.473 -7.555
28 -0.1250000 -0.0018650 -0.757 -0.951 -7.442 -13.031
27 -0.1562500 -0.0029903 -0.949 -1.194 -11.272 -20.233
26 -0.1875000 -0.0044433 -1.143 -1.440 -16.083 -29.396
25 -0.2187500 -0.0062746 -1.338 -1.688 -22.025 -40.823
24 -0.2500000 -0.0085478 -1.534 -1.939 -29.287 -54.895
23 -0.2812500 -0.0113419 -1.733 -2.194 -38.107 -72.083
22 -0.3125000 -0.0147533 -1.934 -2.454 -48.773 -92.966
21 -0.3437500 -0.0188997 -2.138 -2.720 -61.637 -118.248
20 -0.3750000 -0.0239231 -2.346 -2.992 -77.122 -148.778
19 -0.4062500 -0.0299948 -2.557 -3.272 -95.741 -185.581
18 -0.4375000 -0.0373204 -2.774 -3.561 -118.110 -229.889
17 -0.4687500 -0.0461459 -2.997 -3.862 -144.963 -283.172
16 -0.5000000 -0.0567646 -3.227 -4.176 -177.177 -347.187
15 -0.5312500 -0.0695249 -3.466 -4.506 -215.794 -424.018
14 -0.5625000 -0.0848395 -3.714 -4.855 -262.047 -516.135
13 -0.5937500 -0.1031943 -3.975 -5.227 -317.391 -626.446
12 -0.6250000 -0.1251587 -4.251 -5.626 -383.526 -758.359
11 -0.6562500 -0.1513952 -4.543 -6.057 -462.434 -915.837
10 -0.6875000 -0.1826680 -4.856 -6.526 -556.399 -1103.456
9 -0.7187500 -0.2198500 -5.192 -7.039 -668.032 -1326.437
8 -0.7500000 -0.2639239 -5.556 -7.604 -800.270 -1590.663
7 -0.7812500 -0.3159778 -5.951 -8.229 -956.366 -1902.645
6 -0.8125000 -0.3771881 -6.384 -8.923 -1139.838 -2269.423
5 -0.8437500 -0.4487876 -6.858 -9.694 -1354.372 -2698.375
4 -0.8750000 -0.5320097 -7.378 -10.553 -1603.654 -3196.882
3 -0.9062500 -0.6280041 -7.950 -11.507 -1891.123 -3771.824
2 -0.9375000 -0.7377147 -8.576 -12.564 -2219.601 -4428.850
1 -0.9687500 -0.8617150 -9.259 -13.728 -2590.801 -5171.391
0 -1.0000000 -1.0000000 -10.000 -15.000 -3004.710 -5999.420

Time Splitting Weights for Grid 01: ndtfast = 72 nfast = 98
==================================

Primary Secondary Accumulated to Current Step

1-0.0002697221589845 0.0138888888888889-0.0002697221589845 0.0138888888888889
2-0.0005179738466200 0.0138926350299859-0.0007876960056045 0.0277815239188748
3-0.0007447551074202 0.0138998291111890-0.0015324511130247 0.0416813530300638
4-0.0009500661144945 0.0139101729321254-0.0024825172275192 0.0555915259621891
5-0.0011339073179268 0.0139233682948267-0.0036164245454460 0.0695148942570158
6-0.0012962796525070 0.0139391170075756-0.0049127041979530 0.0834540112645914
7-0.0014371848048139 0.0139571208916382-0.0063498890027669 0.0974111321562297
8-0.0015566255396494 0.0139770817917051-0.0079065145424162 0.1113882139479348
9-0.0016546060858249 0.0139987015908669-0.0095611206282412 0.1253869155388017
10-0.0017311325812999 0.0140216822309478-0.0112922532095411 0.1394085977697495
11-0.0017862135776710 0.0140457257390214-0.0130784667872121 0.1534543235087709
12-0.0018198606040138 0.0140705342609335-0.0148983273912259 0.1675248577697044
13-0.0018320887900766 0.0140958101026559-0.0167304161813025 0.1816206678723603
14-0.0018229175488246 0.0141212557802959-0.0185533337301271 0.1957419236526562
15-0.0017923713183374 0.0141465740795851-0.0203457050484645 0.2098884977322413
16-0.0017404803630569 0.0141714681256731-0.0220861854115213 0.2240599658579144
17-0.0016672816343881 0.0141956414640489-0.0237534670459094 0.2382556073219634
18-0.0015728196906507 0.0142187981534154-0.0253262867365601 0.2524744054753787
19-0.0014571476763831 0.0142406428713411-0.0267834344129432 0.2667150483467199
20-0.0013203283609977 0.0142608810335131-0.0281037627739409 0.2809759293802330
21-0.0011624352367880 0.0142792189274159-0.0292661980107289 0.2952551483076488
22-0.0009835536762878 0.0142953638612601-0.0302497516870167 0.3095505121689089
23-0.0007837821489820 0.0143090243289864-0.0310335338359988 0.3238595364978953
24-0.0005632334973686 0.0143199101921667-0.0315967673333673 0.3381794466900619
25-0.0003220362723731 0.0143277328796301-0.0319188036057405 0.3525071795696921
26-0.0000603361281147 0.0143322056056353-0.0319791397338552 0.3668393851753274
27 0.0002217027239766 0.0143330436074147-0.0317574370098786 0.3811724287827420
28 0.0005238960016901 0.0143299644029150-0.0312335410081885 0.3955023931856571
29 0.0008460377792135 0.0143226880695582-0.0303875032289750 0.4098250812552152
30 0.0011878988549600 0.0143109375448469-0.0291996043740150 0.4241360188000621
31 0.0015492250600437 0.0142944389496391-0.0276503793139713 0.4384304577497012
32 0.0019297355074029 0.0142729219349163-0.0257206438065684 0.4527033796846174
33 0.0023291207815727 0.0142461200528690-0.0233915230249958 0.4669494997374865
34 0.0027470410691043 0.0142137711531249-0.0206444819558915 0.4811632708906114
35 0.0031831242296341 0.0141756178049429-0.0174613577262575 0.4953388886955543
36 0.0036369638075999 0.0141314077461980-0.0138243939186576 0.5094702964417523
37 0.0041081169846063 0.0140808943599814-0.0097162769340513 0.5235511908017337
38 0.0045961024724372 0.0140238371796396-0.0051201744616140 0.5375750279813734
39 0.0051003983467181 0.0139600024230780-0.0000197761148959 0.5515350304044513
40 0.0056204398212250 0.0138891635571513 0.0056006637063291 0.5654241939616026
41 0.0061556169628429 0.0138111018929677 0.0117562806691720 0.5792352958545702
42 0.0067052723471717 0.0137256072129282 0.0184615530163437 0.5929609030674984
43 0.0072686986547809 0.0136324784303286 0.0257302516711247 0.6065933814978269
44 0.0078451362081123 0.0135315242823455 0.0335753878792370 0.6201249057801724
45 0.0084337704490309 0.0134225640572328 0.0420091583282679 0.6335474698374053
46 0.0090337293570240 0.0133054283565518 0.0510428876852919 0.6468528981939571
47 0.0096440808080493 0.0131799598932598 0.0606869684933412 0.6600328580872169
48 0.0102638298740300 0.0130460143264814 0.0709507983673713 0.6730788724136982
49 0.0108919160629993 0.0129034611337865 0.0818427144303705 0.6859823335474847
50 0.0115272104998922 0.0127521845218004 0.0933699249302627 0.6987345180692851
51 0.0121685130479867 0.0125920843759686 0.1055384379782495 0.7113266024452537
52 0.0128145493709923 0.0124230772503021 0.1183529873492417 0.7237496796955558
53 0.0134639679357865 0.0122450973979272 0.1318169552850283 0.7359947770934830
54 0.0141153369558014 0.0120580978432635 0.1459322922408297 0.7480528749367465
55 0.0147671412750559 0.0118620514966551 0.1606994335158856 0.7599149264334017
56 0.0154177791928386 0.0116569523122794 0.1761172127087242 0.7715718787456811
57 0.0160655592290373 0.0114428164901566 0.1921827719377615 0.7830146952358378
58 0.0167086968301175 0.0112196837230866 0.2088914687678790 0.7942343789589245
59 0.0173453110157492 0.0109876184893350 0.2262367797836282 0.8052219974482595
60 0.0179734209660816 0.0107467113918941 0.2442102007497098 0.8159687088401536
61 0.0185909425496662 0.0104970805451429 0.2628011432993760 0.8264657893852965
62 0.0191956847920284 0.0102388730097309 0.2819968280914044 0.8367046623950274
63 0.0197853462848869 0.0099722662765083 0.3017821743762913 0.8466769286715357
64 0.0203575115360218 0.0096974698003293 0.3221396859123131 0.8563743984718649
65 0.0209096472597907 0.0094147265845512 0.3430493331721038 0.8657891250564161
66 0.0214390986082931 0.0091243148170541 0.3644884317803970 0.8749134398734703
67 0.0219430853431833 0.0088265495586056 0.3864315171235803 0.8837399894320759
68 0.0224186979481308 0.0085217844843947 0.4088502150717111 0.8922617739164705
69 0.0228628936819303 0.0082104136795596 0.4317131087536414 0.9004721875960301
70 0.0232724925722585 0.0078928734895328 0.4549856013258999 0.9083650610855628
71 0.0236441733500801 0.0075696444260292 0.4786297746759800 0.9159347055115920
72 0.0239744693247019 0.0072412531295003 0.5026042440006819 0.9231759586410923
73 0.0242597641994746 0.0069082743888794 0.5268640082001566 0.9300842330299717
74 0.0244962878281439 0.0065713332194423 0.5513602960283005 0.9366555662494139
75 0.0246801119118487 0.0062311069996069 0.5760404079401492 0.9428866732490209
76 0.0248071456367686 0.0058883276674979 0.6008475535769178 0.9487750009165188
77 0.0248731312524188 0.0055437839780984 0.6257206848293365 0.9543187848946172
78 0.0248736395905937 0.0051983238218148 0.6505943244199301 0.9595171087164319
79 0.0248040655249589 0.0048528566052787 0.6753983899448891 0.9643699653217107
80 0.0246596233712912 0.0045083556952099 0.7000580133161803 0.9688783210169206
81 0.0244353422283667 0.0041658609261642 0.7244933555445470 0.9730441819430847
82 0.0241260612594974 0.0038264811729924 0.7486194168040444 0.9768706631160772
83 0.0237264249147164 0.0034913969888327 0.7723458417187609 0.9803620601049099
84 0.0232308780936106 0.0031618633094617 0.7955767198123714 0.9835239234143716
85 0.0226336612488020 0.0028392122248282 0.8182103810611735 0.9863631356391998
86 0.0219288054300778 0.0025248558185948 0.8401391864912513 0.9888879914577946
87 0.0211101272691677 0.0022202890765104 0.8612493137604190 0.9911082805343050
88 0.0201712239051703 0.0019270928644386 0.8814205376655893 0.9930353733987436
89 0.0191054678506273 0.0016469369768668 0.9005260055162165 0.9946823103756104
90 0.0179060017982467 0.0013815832567192 0.9184320073144632 0.9960638936323296
91 0.0165657333682728 0.0011328887872991 0.9349977406827360 0.9971967824196287
92 0.0150773297965061 0.0009028091571842 0.9500750704792420 0.9980995915768129
93 0.0134332125629703 0.0006934017988994 0.9635082830422124 0.9987929933757124
94 0.0116255519612283 0.0005068294021915 0.9751338350034406 0.9992998227779039
95 0.0096462616083458 0.0003453634027300 0.9847800966117864 0.9996451861806338
96 0.0074869928955037 0.0002113875470585 0.9922670895072901 0.9998565737276923
97 0.0051391293792585 0.0001074015346210 0.9974062188865486 0.9999639752623133
98 0.0025937811134509 0.0000360247376868 0.9999999999999996 1.0000000000000000

ndtfast, nfast = 72 98 nfast/ndtfast = 1.36111

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

1.000000000000 1.025151749473 0.512575874736 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.24456

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

Minimum X-grid spacing, DXmin = 8.31956155E+00 km
Maximum X-grid spacing, DXmax = 9.23851059E+00 km
Minimum Y-grid spacing, DYmin = 8.32752431E+00 km
Maximum Y-grid spacing, DYmax = 9.23803707E+00 km
Minimum Z-grid spacing, DZmin = 1.87946450E-01 m
Maximum Z-grid spacing, DZmax = 8.28028490E+02 m

Minimum barotropic Courant Number = 7.65400315E-03
Maximum barotropic Courant Number = 2.03630778E-01
Maximum Coriolis Courant Number = 2.30020631E-02

NLM: GET_STATE - Read state initial conditions, t = 834 00:00:00
(Grid 01, File: ssc12_rst.nc, Rec=0001, Index=1)
- free-surface
(Min = -1.68187836E+00 Max = 7.37968147E-01)
- vertically integrated u-momentum component
(Min = -1.18573319E+00 Max = 1.04792744E+00)
- vertically integrated v-momentum component
(Min = -7.39885390E-01 Max = 1.14603496E+00)
- u-momentum component
(Min = -2.14761646E+00 Max = 2.94321512E+00)
- v-momentum component
(Min = -3.85626957E+00 Max = 1.83796010E+00)
- potential temperature
(Min = -1.72992955E+00 Max = 3.28646731E+01)
- salinity
(Min = 3.24795329E+01 Max = 3.50028380E+01)
- vertical viscosity coefficient
(Min = 0.00000000E+00 Max = 3.20311947E+02)
- temperature vertical diffusion coefficient
(Min = 0.00000000E+00 Max = 4.33677575E+02)
- salinity vertical diffusion coefficient
(Min = 0.00000000E+00 Max = 4.33677575E+02)
GET_2DFLD - surface u-momentum stress, t = 105 00:00:00
(Rec=0000004, Index=1, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -6.20676099E-05 Max = -1.88522345E-06)
GET_2DFLD - surface v-momentum stress, t = 105 00:00:00
(Rec=0000004, Index=1, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -9.36735938E-05 Max = 2.28887546E-05)
GET_2DFLD - solar shortwave radiation flux, t = 105 00:00:00
(Rec=0000004, Index=1, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 4.67014091E-05 Max = 7.53452299E-05)
GET_2DFLD - surface net heat flux, t = 105 00:00:00
(Rec=0000004, Index=1, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.61152481E-05 Max = 3.52080195E-05)
GET_2DFLD - sea surface temperature climatology, t = 105 00:00:00
(Rec=0000004, Index=1, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.78500004E+01 Max = 3.01499996E+01)
GET_2DFLD - surface net heat flux sensitivity to SST, t = 105 00:00:00
(Rec=0000004, Index=1, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.08742870E-05 Max = -5.57647613E-06)
GET_2DFLD - surface net freswater flux, (E-P), t = 105 00:00:00
(Rec=0000004, Index=1, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -5.87691981E-08 Max = 5.38352896E-08)
GET_NGFLD - free-surface western boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface eastern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface southern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface northern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum western boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum western boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum northern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - temperature western boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.58318077E+01 Max = 2.91239976E+01)
GET_NGFLD - salinity western boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.31880350E+01 Max = 3.39559378E+01)
GET_NGFLD - temperature eastern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.14744140E+00 Max = 2.87537389E+01)
GET_NGFLD - salinity eastern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.39671503E+01 Max = 3.49988728E+01)
GET_NGFLD - temperature southern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.19736436E+00 Max = 2.94572413E+01)
GET_NGFLD - salinity southern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.23207594E+01 Max = 3.50003110E+01)
GET_NGFLD - temperature northern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.08338799E+00 Max = 2.33881566E+01)
GET_NGFLD - salinity northern boundary condition, t = 105 00:00:00
(Rec=0000004, Index=2, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.27818864E+01 Max = 3.48521938E+01)
GET_3DFLD - potential temperature climatology, t = 105 00:00:00
(Rec=0000004, Index=1, File: ssc12_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.12952750E+00 Max = 2.93347012E+01)
GET_3DFLD - salinity climatology, t = 105 00:00:00
(Rec=0000004, Index=1, File: ssc12_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.23027001E+01 Max = 3.50098615E+01)

Basin information for Grid 01:

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

Initial basin volumes: TotVolume = 8.9740602851E+15 m3
MinVolume = 1.3390901769E+07 m3
MaxVolume = 7.0378199510E+10 m3
Max/Min = 5.2556728981E+03

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00200161 - 00424800)

GET_2DFLD - surface u-momentum stress, t = 135 00:00:00
(Rec=0000005, Index=2, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.82850855E-05 Max = 1.04670723E-05)
GET_2DFLD - surface v-momentum stress, t = 135 00:00:00
(Rec=0000005, Index=2, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -4.57473183E-05 Max = 2.86107910E-05)
GET_2DFLD - solar shortwave radiation flux, t = 135 00:00:00
(Rec=0000005, Index=2, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 5.37974143E-05 Max = 7.05541445E-05)
GET_2DFLD - surface net heat flux, t = 135 00:00:00
(Rec=0000005, Index=2, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 5.89348978E-06 Max = 3.80686309E-05)
GET_2DFLD - sea surface temperature climatology, t = 135 00:00:00
(Rec=0000005, Index=2, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.20245279E+01 Max = 3.10499992E+01)
GET_2DFLD - surface net heat flux sensitivity to SST, t = 135 00:00:00
(Rec=0000005, Index=2, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.03150482E-05 Max = -5.28865044E-06)
GET_2DFLD - surface net freswater flux, (E-P), t = 135 00:00:00
(Rec=0000005, Index=2, File: ssc12_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -8.59949345E-08 Max = 3.07742399E-08)
GET_NGFLD - free-surface western boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface eastern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface southern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface northern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum western boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum western boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum northern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - temperature western boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.42429398E+01 Max = 3.00986811E+01)
GET_NGFLD - salinity western boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.28348526E+01 Max = 3.40548760E+01)
GET_NGFLD - temperature eastern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.14744140E+00 Max = 2.92940963E+01)
GET_NGFLD - salinity eastern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.36917025E+01 Max = 3.50467227E+01)
GET_NGFLD - temperature southern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.19736436E+00 Max = 3.01042090E+01)
GET_NGFLD - salinity southern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.28487554E+01 Max = 3.50486625E+01)
GET_NGFLD - temperature northern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.08338799E+00 Max = 2.51038226E+01)
GET_NGFLD - salinity northern boundary condition, t = 135 00:00:00
(Rec=0000005, Index=1, File: ssc12_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.37117703E+01 Max = 3.48547595E+01)
GET_3DFLD - potential temperature climatology, t = 135 00:00:00
(Rec=0000005, Index=2, File: ssc12_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.12952750E+00 Max = 3.01322845E+01)
GET_3DFLD - salinity climatology, t = 135 00:00:00
(Rec=0000005, Index=2, File: ssc12_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.27405017E+01 Max = 3.50487220E+01)

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

200160 834 00:00:00 2.709243E-02 1.908407E+04 1.908410E+04 9.009309E+15
(275,250,27) 1.116631E-01 1.106118E-01 0.000000E+00 3.738429E+00
DEF_HIS - inquiring history file, Grid 01: ssc12_his_0014.nc
DEF_AVG - inquiring average file, Grid 01: ssc12_avg_0000.nc
200161 834 00:06:00 2.709387E-02 1.908407E+04 1.908410E+04 9.009309E+15
(274,250,26) 8.662358E-02 9.662954E-02 2.610078E-01 3.738679E+00

(omit)*********(omit)*****

205303 855 10:18:00 3.526154E-02 1.908405E+04 1.908408E+04 9.009412E+15
(275,253,12) 2.359342E-02 4.247913E-02 5.696632E-01 2.826756E+00
205304 855 10:24:00 3.527035E-02 1.908405E+04 1.908408E+04 9.009412E+15
(275,253,12) 2.475272E-02 4.158159E-02 5.812245E-01 2.823614E+00
205305 855 10:30:00 3.527900E-02 1.908405E+04 1.908408E+04 9.009412E+15
(275,253,12) 2.591119E-02 4.067038E-02 5.919800E-01 2.820217E+00
205306 855 10:36:00 3.528818E-02 1.908405E+04 1.908408E+04 9.009412E+15
(275,253,12) 2.706134E-02 3.970224E-02 6.024620E-01 2.816664E+00
205307 855 10:42:00 3.529777E-02 1.908405E+04 1.908408E+04 9.009411E+15
(275,253,12) 2.819525E-02 3.865272E-02 6.125581E-01 2.817170E+00
205308 855 10:48:00 3.530754E-02 1.908405E+04 1.908408E+04 9.009411E+15
(275,253,12) 2.927960E-02 3.753774E-02 6.226051E-01 2.826378E+00
205309 855 10:54:00 3.531776E-02 1.908405E+04 1.908408E+04 9.009411E+15
(275,253,12) 2.995494E-02 3.636968E-02 6.325719E-01 2.835087E+00
205310 855 11:00:00 3.532825E-02 1.908405E+04 1.908408E+04 9.009411E+15
(275,253,12) 3.099315E-02 3.516521E-02 6.422409E-01 2.843332E+00
205311 855 11:06:00 3.533900E-02 1.908405E+04 1.908408E+04 9.009410E+15
(275,253,12) 3.203498E-02 3.392989E-02 6.524234E-01 2.851146E+00
205312 855 11:12:00 3.534990E-02 1.908405E+04 1.908409E+04 9.009410E+15
(275,253,12) 3.305723E-02 3.272617E-02 6.627129E-01 2.858856E+00
205313 855 11:18:00 3.536101E-02 1.908405E+04 1.908409E+04 9.009410E+15
(275,253,12) 3.405883E-02 3.158176E-02 6.730664E-01 2.866121E+00
205314 855 11:24:00 3.537232E-02 1.908405E+04 1.908409E+04 9.009410E+15
(275,253,12) 3.503133E-02 3.045003E-02 6.835305E-01 2.873413E+00
205315 855 11:30:00 3.538395E-02 1.908405E+04 1.908409E+04 9.009409E+15
(275,253,12) 3.595725E-02 2.926428E-02 6.942394E-01 2.880557E+00
205316 855 11:36:00 3.539564E-02 1.908405E+04 1.908409E+04 9.009409E+15
(275,253,12) 3.683159E-02 2.799544E-02 7.049749E-01 2.889096E+00
205317 855 11:42:00 3.540743E-02 1.908405E+04 1.908409E+04 9.009409E+15
(275,253,12) 3.763728E-02 2.665123E-02 7.154729E-01 2.915840E+00
205318 855 11:48:00 3.541921E-02 1.908405E+04 1.908409E+04 9.009409E+15
(275,253,12) 3.837849E-02 2.526051E-02 7.255248E-01 2.942315E+00
205319 855 11:54:00 3.543090E-02 1.908405E+04 1.908409E+04 9.009408E+15
(275,253,12) 3.903485E-02 2.383714E-02 7.352442E-01 2.968751E+00
205320 855 12:00:00 3.544240E-02 1.908405E+04 1.908409E+04 9.009408E+15
(275,253,12) 3.959201E-02 2.245813E-02 7.445034E-01 2.994961E+00
205321 855 12:06:00 3.545381E-02 1.908405E+04 1.908409E+04 9.009408E+15
(275,253,12) 4.006172E-02 2.117694E-02 7.530285E-01 3.020709E+00
205322 855 12:12:00 3.546518E-02 1.908405E+04 1.908409E+04 9.009408E+15
(275,253,12) 4.041660E-02 1.997603E-02 7.608261E-01 3.046144E+00
205323 855 12:18:00 3.547637E-02 1.908405E+04 1.908409E+04 9.009408E+15
(275,253,13) 4.348721E-02 1.939974E-02 7.636952E-01 3.071278E+00
205324 855 12:24:00 3.548745E-02 1.908405E+04 1.908409E+04 9.009407E+15
(275,253,13) 4.360021E-02 1.811768E-02 7.705619E-01 3.095744E+00
205325 855 12:30:00 3.549807E-02 1.908405E+04 1.908409E+04 9.009407E+15
(275,253,13) 4.368998E-02 1.697693E-02 7.760241E-01 3.119359E+00
205326 855 12:36:00 3.550857E-02 1.908405E+04 1.908409E+04 9.009407E+15
(275,253,13) 4.373112E-02 1.595639E-02 7.796340E-01 3.142272E+00
205327 855 12:42:00 3.551871E-02 1.908405E+04 1.908409E+04 9.009407E+15
(275,253,13) 4.376087E-02 1.508760E-02 7.820757E-01 3.164540E+00
205328 855 12:48:00 3.552841E-02 1.908405E+04 1.908409E+04 9.009406E+15
(275,253,13) 4.377394E-02 1.438975E-02 7.830106E-01 3.186193E+00
205329 855 12:54:00 3.553786E-02 1.908405E+04 1.908409E+04 9.009406E+15
(275,253,13) 4.373087E-02 1.380425E-02 7.824235E-01 3.207101E+00
205330 855 13:00:00 3.554692E-02 1.908405E+04 1.908409E+04 9.009406E+15
(275,253,13) 4.363620E-02 1.326950E-02 7.805699E-01 3.227246E+00
205331 855 13:06:00 3.555572E-02 1.908405E+04 1.908409E+04 9.009406E+15
(275,253,13) 4.347621E-02 1.279885E-02 7.771226E-01 3.246639E+00
205332 855 13:12:00 3.556413E-02 1.908405E+04 1.908409E+04 9.009405E+15
(275,251,11) 8.687961E-03 5.868544E-02 7.606430E-01 3.265299E+00
205333 855 13:18:00 3.557216E-02 1.908405E+04 1.908409E+04 9.009405E+15
(275,251,11) 7.655248E-03 5.796280E-02 7.986381E-01 3.283115E+00
205334 855 13:24:00 3.557991E-02 1.908405E+04 1.908409E+04 9.009405E+15
(275,251,12) 7.625260E-03 5.806786E-02 8.368371E-01 3.299912E+00
205335 855 13:30:00 3.558742E-02 1.908405E+04 1.908409E+04 9.009405E+15
(275,251,12) 6.596264E-03 5.734744E-02 8.772965E-01 3.315893E+00
205336 855 13:36:00 3.559463E-02 1.908405E+04 1.908409E+04 9.009405E+15
(275,251,12) 5.537195E-03 5.650501E-02 9.163505E-01 3.331017E+00
205337 855 13:42:00 3.560156E-02 1.908405E+04 1.908409E+04 9.009404E+15
(275,251,12) 4.468625E-03 5.553002E-02 9.529627E-01 3.345418E+00
205338 855 13:48:00 3.560823E-02 1.908405E+04 1.908409E+04 9.009404E+15
(275,251,13) 3.995803E-03 5.505565E-02 9.875739E-01 3.360165E+00
205339 855 13:54:00 3.561458E-02 1.908405E+04 1.908409E+04 9.009404E+15
(275,251,13) 3.121549E-03 5.403858E-02 1.024790E+00 3.375139E+00
205340 855 14:00:00 3.562080E-02 1.908405E+04 1.908409E+04 9.009404E+15
(275,251,13) 2.427986E-03 5.308032E-02 1.060749E+00 3.389437E+00
205341 855 14:06:00 3.562679E-02 1.908405E+04 1.908409E+04 9.009403E+15
(275,251,13) 1.558878E-03 5.220892E-02 1.093671E+00 3.403209E+00
205342 855 14:12:00 3.563262E-02 1.908405E+04 1.908409E+04 9.009403E+15
(275,251,14) 9.068599E-04 5.085986E-02 1.123449E+00 3.417124E+00
205343 855 14:18:00 3.563828E-02 1.908405E+04 1.908409E+04 9.009403E+15
(275,251,14) 6.136644E-04 4.993357E-02 1.155907E+00 3.431265E+00
205344 855 14:24:00 3.564350E-02 1.908405E+04 1.908409E+04 9.009403E+15
(275,251,14) 8.381187E-05 4.897187E-02 1.184329E+00 3.445476E+00
205345 855 14:30:00 3.564900E-02 1.908405E+04 1.908409E+04 9.009402E+15
(275,251,14) 1.176028E-03 4.829006E-02 1.205966E+00 3.460747E+00
205346 855 14:36:00 3.565424E-02 1.908405E+04 1.908409E+04 9.009402E+15
(275,251,15) 3.141366E-03 4.644630E-02 1.233495E+00 3.475221E+00
205347 855 14:42:00 3.565942E-02 1.908406E+04 1.908409E+04 9.009402E+15
(275,251,15) 2.754590E-03 4.529407E-02 1.254921E+00 3.490049E+00
205348 855 14:48:00 3.566471E-02 1.908406E+04 1.908409E+04 9.009402E+15
(275,251,15) 3.662880E-03 4.342926E-02 1.270466E+00 3.507288E+00
205349 855 14:54:00 3.567121E-02 1.908405E+04 1.908409E+04 9.009401E+15
(275,251,16) 5.686473E-03 4.105372E-02 1.295611E+00 3.530932E+00
205350 855 15:00:00 3.568281E-02 1.908405E+04 1.908409E+04 9.009401E+15
(275,251,16) 9.154510E-03 3.918325E-02 1.337819E+00 3.707565E+00
205351 855 15:06:00 3.572494E-02 1.908405E+04 1.908409E+04 9.009400E+15
(275,251,17) 2.925438E-02 3.822207E-02 1.486118E+00 4.974890E+00
205352 855 15:12:00 3.620864E-02 2.626193E+04 2.626196E+04 9.009397E+15
(275,251,17) 1.398685E-01 3.509188E-02 2.106421E+00 1.220272E+01

Blowing-up: Saving latest model state into RESTART file

WRT_RST - wrote re-start fields (Index=1,1) into time record = 0000003

Elapsed CPU time (seconds):

Node # 0 CPU: 1878.147
Node # 3 CPU: 1895.773
Node # 1 CPU: 1895.909
Node # 7 CPU: 1896.513
Node # 2 CPU: 1896.583
Node # 4 CPU: 1896.520
Node # 5 CPU: 1896.442
Node # 6 CPU: 1896.219
Node # 9 CPU: 1896.308
Node # 10 CPU: 1896.471
Node # 12 CPU: 1896.603
Node # 13 CPU: 1896.665
Node # 14 CPU: 1896.622
Total: 30322.101
Node # 15 CPU: 1896.529

Nonlinear model elapsed time profile, Grid: 01

Allocation and array initialization .............. 1.099 ( 0.0036 %)
Node # 8 CPU: 1894.255
Ocean state initialization ....................... 0.299 ( 0.0010 %)
Reading of input data ............................ 18.717 ( 0.0617 %)
Processing of input data ......................... 531.208 ( 1.7519 %)
Processing of output time averaged data .......... 0.400 ( 0.0013 %)
Computation of vertical boundary conditions ...... 19.930 ( 0.0657 %)
Computation of global information integrals ...... 371.286 ( 1.2245 %)
Writing of output data ........................... 336.025 ( 1.1082 %)
Model 2D kernel .................................. 14348.311 (47.3196 %)
2D/3D coupling, vertical metrics ................. 788.128 ( 2.5992 %)
Omega vertical velocity .......................... 493.266 ( 1.6268 %)
Equation of state for seawater ................... 660.155 ( 2.1771 %)
My2.5 vertical mixing parameterization ........... 3560.626 (11.7427 %)
Node # 11 CPU: 1896.546
3D equations right-side terms .................... 760.302 ( 2.5074 %)
3D equations predictor step ...................... 1985.028 ( 6.5465 %)
Pressure gradient ................................ 890.511 ( 2.9368 %)
Harmonic mixing of tracers, geopotentials ........ 856.133 ( 2.8235 %)
Harmonic stress tensor, geopotentials ............ 1612.028 ( 5.3163 %)
Corrector time-step for 3D momentum .............. 1045.567 ( 3.4482 %)
Corrector time-step for tracers .................. 1509.167 ( 4.9771 %)
Total: 29788.185 98.2392

Nonlinear model message Passage profile, Grid: 01

Message Passage: 2D halo exchanges ............... 1790.132 ( 5.9037 %)
Message Passage: 3D halo exchanges ............... 817.576 ( 2.6963 %)
Message Passage: 4D halo exchanges ............... 200.821 ( 0.6623 %)
Message Passage: data broadcast .................. 329.280 ( 1.0859 %)
Message Passage: data reduction .................. 21.304 ( 0.0703 %)
Message Passage: data gathering .................. 34.376 ( 0.1134 %)
Message Passage: data scattering.................. 11.133 ( 0.0367 %)
Total: 3204.622 10.5686

All percentages are with respect to total time = 30322.101

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

Analytical header files used:

/vol/home/student/lshw/roms/MYPROJECT/Functionals/ana_btflux.h
/vol/home/student/lshw/roms/MYPROJECT/Functionals/ana_sponge.h
ROMS/Functionals/ana_nudgcoef.h

ROMS/TOMS: DONE... Sunday - June 12, 2016 - 1:04:46 PM

kate · #19 Unread post by **kate** » Thu Jun 16, 2016 3:54 pm

No need to double post.

Ocean Modeling Discussion

MAIN: Abnormal termination: BLOWUP.

MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.

Re: MAIN: Abnormal termination: BLOWUP.