Blowing-up: Saving latest model state into RESTART file

[shiyf]

My model is blowing and I don't know why. This problem seems to have been discussed many times.

Reference before discussion,I changed my time step, cheaked the grid file, DxMIN,DYmin,DZmin,and boundary conditions.

However, the model is still blowing, I don't know how to solve this problem,
could you give me some useful advice.



Model Input Parameters: ROMS/TOMS version 3.7
Monday - June 26, 2017 - 2:28:18 AM
-----------------------------------------------------------------------------

SHIYF

Operating system : Linux
CPU/hardware : x86_64
Compiler system : ifort
Compiler command : /opt/sgi/mpt/mpt-2.08/bin/mpif90
Compiler flags : -heap-arrays -fp-model precise -ip -O3 -free -free -free

Input Script :

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

Local Root : /IS5000/users/shiyf/roms/trunk
Header Dir : /IS5000/users/shiyf/roms/test1
Header file : shiyf.h
Analytical Dir: /IS5000/users/shiyf/roms/test1

Resolution, Grid 01: 0199x0259x020, Parallel Nodes: 20, Tiling: 004x005


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

3144960 ntimes Number of timesteps for 3-D equations.
10.000 dt Timestep size (s) for 3-D equations.
15 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.
60 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
360 nHIS Number of timesteps between the writing fields
into history file.
8640 ndefHIS Number of timesteps between creation of new
history files.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
360 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
8640 ndefAVG Number of timesteps between creation of new
time-averaged file.
0.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
0.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
5.0000E+00 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
F LuvSponge Turning OFF sponge on horizontal momentum.
F LtracerSponge(01) Turning OFF sponge on tracer 01: temp
F LtracerSponge(02) Turning OFF sponge on tracer 02: salt
1.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
1.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
1.0000E-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
5.0000E-06 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
5.0000E-06 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
3.000 gls_p GLS stability exponent.
1.500 gls_m GLS turbulent kinetic energy exponent.
-1.000 gls_n GLS turbulent length scale exponent.
7.6000E-06 gls_Kmin GLS minimum value of turbulent kinetic energy.
1.0000E-12 gls_Pmin GLS minimum value of dissipation.
5.4770E-01 gls_cmu0 GLS stability coefficient.
1.4400E+00 gls_c1 GLS shear production coefficient.
1.9200E+00 gls_c2 GLS dissipation coefficient.
-4.0000E-01 gls_c3m GLS stable buoyancy production coefficient.
1.0000E+00 gls_c3p GLS unstable buoyancy production coefficient.
1.0000E+00 gls_sigk GLS constant Schmidt number for TKE.
1.3000E+00 gls_sigp GLS constant Schmidt number for PSI.
1400.000 charnok_alpha Charnok factor for Zos calculation.
0.500 zos_hsig_alpha Factor for Zos calculation using Hsig(Awave).
0.250 sz_alpha Factor for Wave dissipation surface tke flux .
100.000 crgban_cw Factor for Craig/Banner surface tke flux.
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).
2.0000E-02 Zos Surface roughness (m).
1 lmd_Jwt Jerlov water type.
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
7.0000E+00 theta_s S-coordinate surface control parameter.
1.0000E-01 theta_b S-coordinate bottom control parameter.
200.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
0.125 dstart Time-stamp assigned to model initialization (days).
20160101.00 time_ref Reference time for units attribute (yyyymmdd.dd)
2.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
2.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.
T VolCons(1) NLM western edge boundary volume conservation.
T VolCons(2) NLM southern edge boundary volume conservation.
T VolCons(3) NLM eastern edge boundary volume conservation.
T VolCons(4) NLM northern edge boundary volume conservation.
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.
F LtracerCLM(01) Turning OFF processing of climatology tracer 01: temp
F LtracerCLM(02) Turning OFF processing of climatology tracer 02: salt
F LnudgeM2CLM Turning OFF nudging of 2D momentum climatology.
F LnudgeM3CLM Turning OFF nudging of 3D momentum climatology.
F LnudgeTCLM(01) Turning OFF nudging of climatology tracer 01: temp
F LnudgeTCLM(02) Turning OFF nudging of climatology tracer 02: salt
T Hout(idFsur) Write out free-surface.
T Hout(idUbar) Write out 2D U-momentum component.
T Hout(idVbar) Write out 2D V-momentum component.
T Hout(idUvel) Write out 3D U-momentum component.
T Hout(idVvel) Write out 3D V-momentum component.
T Hout(idWvel) Write out W-momentum component.
T Hout(idOvel) Write out omega vertical velocity.
T Hout(idTvar) Write out tracer 01: temp
T Hout(idTvar) Write out tracer 02: salt

T Aout(idFsur) Write out averaged free-surface.
T Aout(idUbar) Write out averaged 2D U-momentum component.
T Aout(idVbar) Write out averaged 2D V-momentum component.
T Aout(idUvel) Write out averaged 3D U-momentum component.
T Aout(idVvel) Write out averaged 3D V-momentum component.
T Aout(idWvel) Write out averaged W-momentum component.
T Aout(idOvel) Write out averaged omega vertical velocity.
T Aout(idTvar) Write out averaged tracer 01: temp
T Aout(idTvar) Write out averaged tracer 02: salt

Output/Input Files:

Output Restart File: ocean_rst.nc
Prefix for History Files: ocean_his
Prefix for Averages Files: ocean_avg
Input Grid File: /IS5000/users/shiyf/roms/test1/roms_grid.nc
Input Nonlinear Initial File: /IS5000/users/shiyf/roms/test1/roms_ini_linear.nc
Input Forcing File 01: /IS5000/users/shiyf/roms/test1/gom_shflux_era.nc
Input Forcing File 02: /IS5000/users/shiyf/roms/test1/gom_sms_era.nc
Input Forcing File 03: /IS5000/users/shiyf/roms/test1/gom_SST_era.nc
Input Forcing File 04: /IS5000/users/shiyf/roms/test1/gom_swflux_era.nc
Input Forcing File 05: /IS5000/users/shiyf/roms/test1/gom_swrad_era.nc
Input Boundary File: /IS5000/users/shiyf/roms/test1/roms_bry.nc

Tile partition information for Grid 01: 0199x0259x0020 tiling: 004x005

tile Istr Iend Jstr Jend Npts

0 1 50 1 52 52000
1 51 100 1 52 52000
2 101 150 1 52 52000
3 151 199 1 52 50960
4 1 50 53 104 52000
5 51 100 53 104 52000
6 101 150 53 104 52000
7 151 199 53 104 50960
8 1 50 105 156 52000
9 51 100 105 156 52000
10 101 150 105 156 52000
11 151 199 105 156 50960
12 1 50 157 208 52000
13 51 100 157 208 52000
14 101 150 157 208 52000
15 151 199 157 208 50960
16 1 50 209 259 51000
17 51 100 209 259 51000
18 101 150 209 259 51000
19 151 199 209 259 49980

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 50.50 0.50 52.50 RHO-points
1 50.50 100.50 0.50 52.50 RHO-points
2 100.50 150.50 0.50 52.50 RHO-points
3 150.50 199.50 0.50 52.50 RHO-points
4 0.50 50.50 52.50 104.50 RHO-points
5 50.50 100.50 52.50 104.50 RHO-points
6 100.50 150.50 52.50 104.50 RHO-points
7 150.50 199.50 52.50 104.50 RHO-points
8 0.50 50.50 104.50 156.50 RHO-points
9 50.50 100.50 104.50 156.50 RHO-points
10 100.50 150.50 104.50 156.50 RHO-points
11 150.50 199.50 104.50 156.50 RHO-points
12 0.50 50.50 156.50 208.50 RHO-points
13 50.50 100.50 156.50 208.50 RHO-points
14 100.50 150.50 156.50 208.50 RHO-points
15 150.50 199.50 156.50 208.50 RHO-points
16 0.50 50.50 208.50 259.50 RHO-points
17 50.50 100.50 208.50 259.50 RHO-points
18 100.50 150.50 208.50 259.50 RHO-points
19 150.50 199.50 208.50 259.50 RHO-points

0 1.00 50.50 0.50 52.50 U-points
1 50.50 100.50 0.50 52.50 U-points
2 100.50 150.50 0.50 52.50 U-points
3 150.50 199.00 0.50 52.50 U-points
4 1.00 50.50 52.50 104.50 U-points
5 50.50 100.50 52.50 104.50 U-points
6 100.50 150.50 52.50 104.50 U-points
7 150.50 199.00 52.50 104.50 U-points
8 1.00 50.50 104.50 156.50 U-points
9 50.50 100.50 104.50 156.50 U-points
10 100.50 150.50 104.50 156.50 U-points
11 150.50 199.00 104.50 156.50 U-points
12 1.00 50.50 156.50 208.50 U-points
13 50.50 100.50 156.50 208.50 U-points
14 100.50 150.50 156.50 208.50 U-points
15 150.50 199.00 156.50 208.50 U-points
16 1.00 50.50 208.50 259.50 U-points
17 50.50 100.50 208.50 259.50 U-points
18 100.50 150.50 208.50 259.50 U-points
19 150.50 199.00 208.50 259.50 U-points

0 0.50 50.50 1.00 52.50 V-points
1 50.50 100.50 1.00 52.50 V-points
2 100.50 150.50 1.00 52.50 V-points
3 150.50 199.50 1.00 52.50 V-points
4 0.50 50.50 52.50 104.50 V-points
5 50.50 100.50 52.50 104.50 V-points
6 100.50 150.50 52.50 104.50 V-points
7 150.50 199.50 52.50 104.50 V-points
8 0.50 50.50 104.50 156.50 V-points
9 50.50 100.50 104.50 156.50 V-points
10 100.50 150.50 104.50 156.50 V-points
11 150.50 199.50 104.50 156.50 V-points
12 0.50 50.50 156.50 208.50 V-points
13 50.50 100.50 156.50 208.50 V-points
14 100.50 150.50 156.50 208.50 V-points
15 150.50 199.50 156.50 208.50 V-points
16 0.50 50.50 208.50 259.00 V-points
17 50.50 100.50 208.50 259.00 V-points
18 100.50 150.50 208.50 259.00 V-points
19 150.50 199.50 208.50 259.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 124
HaloSizeJ(1) = 128
TileSide(1) = 58
TileSize(1) = 3248


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

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

zeta 1 Closed Clamped Clamped 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

tke 1 Closed Clamped Clamped Closed

Activated C-preprocessing Options:

SHIYF SHIYF
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
CURVGRID Orthogonal curvilinear grid.
DIFF_GRID Horizontal diffusion coefficient scaled by grid size.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
GLS_MIXING Generic Length-Scale turbulence closure.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
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.
NO_WRITE_GRID Not Writing grid arrays into NetCDF ouput files.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
RADIATION_2D Use tangential phase speed in radiation conditions.
RI_SPLINES Parabolic Spline Reconstruction for Richardson Number.
RST_SINGLE Single 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_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
UV_ADV Advection of momentum.
UV_COR Coriolis term.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_C4VADVECTION Fourth-order centered vertical advection of momentum.
UV_LDRAG Linear 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= 70272) is active.
Node # 3 (pid= 70275) is active.
Node # 19 (pid= 70291) is active.
Node # 11 (pid= 70283) is active.
Node # 14 (pid= 70286) is active.
Node # 18 (pid= 70290) is active.
Node # 1 (pid= 70273) is active.
Node # 2 (pid= 70274) is active.
Node # 4 (pid= 70276) is active.
Node # 5 (pid= 70277) is active.
Node # 6 (pid= 70278) is active.
Node # 7 (pid= 70279) is active.
Node # 8 (pid= 70280) is active.
Node # 9 (pid= 70281) is active.
Node # 10 (pid= 70282) is active.
Node # 12 (pid= 70284) is active.
Node # 13 (pid= 70285) is active.
Node # 16 (pid= 70288) is active.
Node # 17 (pid= 70289) is active.
Node # 15 (pid= 70287) 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

20 0.0000000 0.0000000 0.000 0.000 0.000 0.000
19 -0.0500000 -0.0001188 -0.909 -5.012 -9.914 -10.602
18 -0.1000000 -0.0004899 -1.819 -10.049 -20.723 -23.035
17 -0.1500000 -0.0011592 -2.729 -15.116 -32.590 -37.635
16 -0.2000000 -0.0022095 -3.640 -20.221 -45.809 -55.001
15 -0.2500000 -0.0037705 -4.552 -25.377 -60.841 -76.076
14 -0.3000000 -0.0060353 -5.466 -30.604 -78.370 -102.262
13 -0.3500000 -0.0092834 -6.381 -35.928 -99.389 -135.590
12 -0.4000000 -0.0139158 -7.298 -41.392 -125.320 -178.970
11 -0.4500000 -0.0205034 -8.219 -47.050 -158.189 -236.549
10 -0.5000000 -0.0298569 -9.145 -52.986 -200.874 -314.215
9 -0.5500000 -0.0431244 -10.078 -59.312 -257.448 -420.305
8 -0.6000000 -0.0619292 -11.022 -66.193 -333.671 -566.608
7 -0.6500000 -0.0885607 -11.979 -73.856 -437.669 -769.751
6 -0.7000000 -0.1262391 -12.957 -82.624 -580.869 -1053.118
5 -0.7500000 -0.1794756 -13.963 -92.948 -779.279 -1449.472
4 -0.8000000 -0.2545552 -15.008 -105.456 -1055.202 -2004.455
3 -0.8500000 -0.3601634 -16.109 -121.016 -1439.461 -2781.143
2 -0.9000000 -0.5081650 -17.288 -140.817 -1974.160 -3865.701
1 -0.9500000 -0.7144926 -18.572 -166.449 -2715.836 -5373.835
0 -1.0000000 -1.0000000 -20.000 -200.000 -3738.494 -7456.989

Time Splitting Weights for Grid 01: ndtfast = 15 nfast = 22
==================================

Primary Secondary Accumulated to Current Step

1-0.0006740978872509 0.0666666666666667-0.0006740978872509 0.0666666666666667
2-0.0001913161153518 0.0667116065258167-0.0008654140026027 0.1333782731924834
3 0.0014472038369421 0.0667243609335069 0.0005817898343395 0.2001026341259903
4 0.0042370228855824 0.0666278806777107 0.0048188127199218 0.2667305148037010
5 0.0081665994120431 0.0663454124853385 0.0129854121319649 0.3330759272890396
6 0.0132119623624627 0.0658009725245357 0.0261973744944276 0.3988768998135752
7 0.0193298623751124 0.0649201750337048 0.0455272368695400 0.4637970748472801
8 0.0264494009361902 0.0636315175420307 0.0719766378057303 0.5274285923893107
9 0.0344621375639418 0.0618682241462847 0.1064387753696720 0.5892968165355954
10 0.0432106750211074 0.0595707483086885 0.1496494503907794 0.6488675648442840
11 0.0524757225556949 0.0566900366406147 0.2021251729464743 0.7055576014848987
12 0.0619616371700794 0.0531916551369017 0.2640868101165537 0.7587492566218004
13 0.0712804429184289 0.0490608793255631 0.3353672530349826 0.8078101359473635
14 0.0799343282324558 0.0443088497976678 0.4153015812674384 0.8521189857450313
15 0.0872966212754959 0.0389798945821708 0.5025982025429343 0.8910988803272021
16 0.0925912433249121 0.0331601198304711 0.5951894458678464 0.9242590001576731
17 0.0948706401828253 0.0269873702754769 0.6900600860506717 0.9512463704331501
18 0.0929921916151718 0.0206626609299552 0.7830522776658435 0.9719090313631054
19 0.0855930988190858 0.0144631814889438 0.8686453764849293 0.9863722128520491
20 0.0710637499186094 0.0087569749010047 0.9397091264035388 0.9951291877530538
21 0.0475195634887281 0.0040193915730974 0.9872286898922669 0.9991485793261512
22 0.0127713101077333 0.0008514206738489 1.0000000000000002 1.0000000000000000

ndtfast, nfast = 15 22 nfast/ndtfast = 1.46667

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

1.000000000000 1.070450639131 0.535225319565 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.09467

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

Minimum X-grid spacing, DXmin = 8.39240211E+00 km
Maximum X-grid spacing, DXmax = 1.07411362E+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 = 9.09306915E-01 m
Maximum Z-grid spacing, DZmax = 2.08315423E+03 m

Minimum barotropic Courant Number = 1.20930843E-03
Maximum barotropic Courant Number = 2.41567963E-02
Maximum Coriolis Courant Number = 9.56534064E-04

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

NLM: GET_STATE - Read state initial conditions, t = 0 03:00:00
(Grid 01, File: roms_ini_linear.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 = 2.84540166E+01)
- salinity
(Min = 0.00000000E+00 Max = 3.49894417E+01)
GET_2DFLD - surface u-momentum stress, t = 0 01:30:00
(Rec=0000001, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 365.9375)
(Min = -7.72763229E-04 Max = 4.49069785E-04)
GET_2DFLD - surface v-momentum stress, t = 0 01:30:00
(Rec=0000001, Index=1, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 365.9375)
(Min = -4.32275935E-04 Max = 7.17418147E-04)
GET_2DFLD - solar shortwave radiation flux, t = 0 01:30:00
(Rec=0000001, Index=1, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 365.9375)
(Min = 3.08818936E-05 Max = 1.46722456E-04)
GET_2DFLD - surface net heat flux, t = 0 01:30:00
(Rec=0000001, Index=1, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 365.9375)
(Min = -4.02489277E-05 Max = 9.00919604E-05)
GET_2DFLD - surface net freswater flux, (E-P), t = 0 01:30:00
(Rec=0000001, Index=1, File: gom_swflux_era.nc)
(Tmin= 0.0625 Tmax= 365.9375)
(Min = -2.30701820E-07 Max = 1.71600979E-07)
GET_NGFLD - free-surface eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -2.64999811E-02 Max = 1.27974999E+00)
GET_NGFLD - free-surface southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 7.58500040E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -1.36086360E-01 Max = 3.23839605E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -6.63492829E-02 Max = 1.19354643E-01)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -7.26458132E-02 Max = 4.06116545E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -6.86429963E-02 Max = 8.93137828E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -9.17234361E-01 Max = 1.03719985E+00)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -4.95005399E-01 Max = 4.40441459E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -6.32066369E-01 Max = 3.68236303E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -3.90926629E-01 Max = 4.29749966E-01)
GET_NGFLD - temperature eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -2.06360966E-01 Max = 2.80411301E+01)
GET_NGFLD - salinity eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 3.52298317E+01)
GET_NGFLD - temperature southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -1.81017101E-01 Max = 2.85190125E+01)
GET_NGFLD - salinity southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 3.56930046E+01)

Basin information for Grid 01:

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

Initial basin volumes: TotVolume = 1.1284720359E+16 m3
MinVolume = 8.5117018506E+07 m3
MaxVolume = 2.3425697327E+11 m3
Max/Min = 2.7521755036E+03

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

GET_2DFLD - surface u-momentum stress, t = 0 04:30:00
(Rec=0000002, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 365.9375)
(Min = -1.02063493E-03 Max = 8.13465144E-04)
GET_2DFLD - surface v-momentum stress, t = 0 04:30:00
(Rec=0000002, Index=2, File: gom_sms_era.nc)
(Tmin= 0.0625 Tmax= 365.9375)
(Min = -4.16719957E-04 Max = 1.03997598E-03)
GET_2DFLD - solar shortwave radiation flux, t = 0 04:30:00
(Rec=0000002, Index=2, File: gom_swrad_era.nc)
(Tmin= 0.0625 Tmax= 365.9375)
(Min = 3.71969539E-05 Max = 1.63637222E-04)
GET_2DFLD - surface net heat flux, t = 0 04:30:00
(Rec=0000002, Index=2, File: gom_shflux_era.nc)
(Tmin= 0.0625 Tmax= 365.9375)
(Min = -1.98255670E-05 Max = 1.29853056E-04)
GET_2DFLD - surface net freswater flux, (E-P), t = 0 04:30:00
(Rec=0000002, Index=2, File: gom_swflux_era.nc)
(Tmin= 0.0625 Tmax= 365.9375)
(Min = -2.94210686E-07 Max = 1.59138668E-07)
GET_NGFLD - free-surface eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 1.29525006E+00)
GET_NGFLD - free-surface southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 7.86999702E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -1.32524595E-01 Max = 3.17164749E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -6.48542270E-02 Max = 9.82362702E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -6.49622902E-02 Max = 3.82685363E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -9.62820500E-02 Max = 9.34708491E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -6.93626523E-01 Max = 1.17415941E+00)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -4.41500038E-01 Max = 3.29676598E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -5.44846058E-01 Max = 5.32745183E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -3.51513118E-01 Max = 3.22937965E-01)
GET_NGFLD - temperature eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -2.06360966E-01 Max = 2.79875050E+01)
GET_NGFLD - salinity eastern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 3.51747093E+01)
GET_NGFLD - temperature southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -1.81132615E-01 Max = 2.84168167E+01)
GET_NGFLD - salinity southern boundary condition, t = 1 00:00:00
(Rec=0000002, Index=1, File: roms_bry.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 3.57525291E+01)

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

0 0 03:00:00 9.755565E-06 2.306463E+04 2.306463E+04 1.131373E+16
(199,180,20) 5.557101E-04 0.000000E+00 0.000000E+00 7.329277E-01
DEF_HIS - creating history file, Grid 01: ocean_his_0001.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file, Grid 01: ocean_avg_0001.nc
1 0 03:00:10 1.259702E-03 2.306463E+04 2.306463E+04 1.131373E+16
(199,180,12) 5.835219E-04 1.257146E-06 5.050335E-03 1.288762E+00
2 0 03:00:20 4.658913E-03 2.306462E+04 2.306463E+04 1.131373E+16
(185,165,20) 1.745415E-05 2.990607E-05 8.671431E-01 2.423018E+00
3 0 03:00:30 9.407770E-03 2.306462E+04 2.306462E+04 1.131373E+16
(185,165,20) 5.332812E-05 9.564012E-05 1.524737E+00 3.294906E+00
4 0 03:00:40 1.466611E-02 2.306461E+04 2.306462E+04 1.131373E+16
(185,165,20) 1.145112E-04 2.117001E-04 1.830163E+00 3.864026E+00
5 0 03:00:50 1.994410E-02 2.306460E+04 2.306462E+04 1.131373E+16
(185,165,20) 1.993618E-04 3.810129E-04 1.729009E+00 4.163179E+00
6 0 03:01:00 2.522976E-02 2.306459E+04 2.306461E+04 1.131373E+16
(085,035,20) 5.850398E-04 2.208341E-04 1.295433E+00 4.310033E+00
7 0 03:01:10 3.083960E-02 2.306457E+04 2.306460E+04 1.131373E+16
(185,166,20) 1.166957E-03 7.500917E-04 1.074155E+00 5.084729E+00
8 0 03:01:20 3.714294E-02 2.306455E+04 2.306459E+04 1.131373E+16
(185,166,20) 1.558458E-03 9.072852E-04 9.010963E-01 6.331005E+00
9 0 03:01:30 4.434111E-02 2.306452E+04 2.306456E+04 1.131373E+16
(084,025,20) 1.487702E-03 6.875167E-05 9.478983E-01 7.402848E+00
10 0 03:01:40 5.240882E-02 2.306449E+04 2.306454E+04 1.131373E+16
(185,165,20) 6.468381E-04 1.710845E-03 1.215620E+00 8.248628E+00
11 0 03:01:50 6.117435E-02 2.306445E+04 2.306451E+04 1.131373E+16
(185,165,20) 6.712875E-04 1.987528E-03 1.281609E+00 8.876819E+00
12 0 03:02:00 7.045034E-02 2.306440E+04 2.306447E+04 1.131373E+16
(085,035,20) 2.105282E-03 3.678537E-04 1.280306E+00 9.346736E+00
13 0 03:02:10 8.012914E-02 2.306436E+04 2.306444E+04 1.131373E+16
(084,036,20) 2.289129E-03 2.694075E-03 1.308140E+00 9.747683E+00
14 0 03:02:20 9.020731E-02 2.306431E+04 2.306440E+04 1.131373E+16
(084,036,20) 2.524105E-03 2.958317E-03 1.164113E+00 1.017282E+01
15 0 03:02:30 1.007575E-01 2.306426E+04 2.306436E+04 1.131373E+16
(072,078,20) 1.260214E-03 7.331314E-04 9.014327E-01 1.069599E+01
16 0 03:02:40 1.118682E-01 2.306421E+04 2.306432E+04 1.131373E+16
(072,078,20) 1.418875E-03 8.514047E-04 9.218481E-01 1.135602E+01
17 0 03:02:50 1.235922E-01 2.306416E+04 2.306428E+04 1.131373E+16
(079,040,20) 3.397167E-03 3.012913E-03 8.764593E-01 1.215272E+01
18 0 03:03:00 1.359219E-01 2.306410E+04 2.306424E+04 1.131373E+16
(079,040,20) 3.586341E-03 3.186580E-03 8.709761E-01 1.305133E+01
19 0 03:03:10 1.487989E-01 2.306405E+04 2.306420E+04 1.131373E+16
(084,036,20) 3.741986E-03 4.064001E-03 1.009707E+00 1.399740E+01
20 0 03:03:20 1.621487E-01 2.306399E+04 2.306416E+04 1.131373E+16
(084,036,20) 3.984258E-03 4.296965E-03 1.204302E+00 1.493381E+01
21 0 03:03:30 1.759155E-01 2.306394E+04 2.306411E+04 1.131373E+16
(084,036,20) 4.220663E-03 4.543542E-03 1.195340E+00 1.581677E+01
22 0 03:03:40 1.900799E-01 2.306388E+04 2.306407E+04 1.131373E+16
(084,036,20) 4.451258E-03 4.801508E-03 9.934283E-01 1.662483E+01
23 0 03:03:50 2.046534E-01 2.306381E+04 2.306402E+04 1.131373E+16
(070,079,20) 1.063528E-03 1.024081E-03 7.611617E-01 1.736059E+01
24 0 03:04:00 2.196600E-01 2.306375E+04 2.306397E+04 1.131373E+16
(079,040,20) 4.985810E-03 4.319316E-03 8.383643E-01 1.804453E+01
25 0 03:04:10 2.351170E-01 2.306368E+04 2.306391E+04 1.131373E+16
(079,040,20) 5.251375E-03 4.513772E-03 8.581952E-01 1.870506E+01
26 0 03:04:20 2.510258E-01 2.306360E+04 2.306385E+04 1.131373E+16
(079,040,20) 5.510188E-03 4.705179E-03 7.382203E-01 1.936843E+01
27 0 03:04:30 2.673698E-01 2.306353E+04 2.306379E+04 1.131373E+16
(083,037,20) 5.755508E-03 5.568097E-03 8.896717E-01 2.005144E+01

Blowing-up: Saving latest model state into RESTART file

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

Elapsed CPU time (seconds):

Node # 0 CPU: 4.660
Node # 3 CPU: 4.690
Node # 6 CPU: 4.690
Node # 7 CPU: 4.700
Node # 1 CPU: 4.680
Node # 2 CPU: 4.680
Node # 4 CPU: 4.690
Node # 5 CPU: 4.690
Node # 8 CPU: 4.690
Node # 9 CPU: 4.670
Node # 10 CPU: 4.690
Node # 11 CPU: 4.690
Node # 12 CPU: 4.690
Node # 13 CPU: 4.690
Node # 14 CPU: 4.690
Node # 15 CPU: 4.690
Node # 16 CPU: 4.700
Node # 17 CPU: 4.680
Node # 18 CPU: 4.690
Node # 19 CPU: 4.700
Total: 93.750

Nonlinear model elapsed time profile, Grid: 01

Allocation and array initialization .............. 0.830 ( 0.8853 %)
Ocean state initialization ....................... 0.090 ( 0.0960 %)
Reading of input data ............................ 6.110 ( 6.5173 %)
Processing of input data ......................... 0.410 ( 0.4373 %)
Processing of output time averaged data .......... 0.930 ( 0.9920 %)
Computation of vertical boundary conditions ...... 0.150 ( 0.1600 %)
Computation of global information integrals ...... 1.040 ( 1.1093 %)
Writing of output data ........................... 5.000 ( 5.3333 %)
Model 2D kernel .................................. 13.710 (14.6240 %)
2D/3D coupling, vertical metrics ................. 2.020 ( 2.1547 %)
Omega vertical velocity .......................... 1.210 ( 1.2907 %)
Equation of state for seawater ................... 3.200 ( 3.4133 %)
GLS vertical mixing parameterization ............. 25.570 (27.2747 %)
3D equations right-side terms .................... 2.940 ( 3.1360 %)
3D equations predictor step ...................... 6.240 ( 6.6560 %)
Pressure gradient ................................ 2.240 ( 2.3893 %)
Harmonic mixing of tracers, geopotentials ........ 3.740 ( 3.9893 %)
Harmonic stress tensor, S-surfaces ............... 1.780 ( 1.8987 %)
Corrector time-step for 3D momentum .............. 3.030 ( 3.2320 %)
Corrector time-step for tracers .................. 2.390 ( 2.5493 %)
Total: 82.630 88.1387

Nonlinear model message Passage profile, Grid: 01

Message Passage: 2D halo exchanges ............... 3.360 ( 3.5840 %)
Message Passage: 3D halo exchanges ............... 2.900 ( 3.0933 %)
Message Passage: 4D halo exchanges ............... 0.920 ( 0.9813 %)
Message Passage: data broadcast .................. 11.990 (12.7893 %)
Message Passage: data reduction .................. 1.330 ( 1.4187 %)
Message Passage: data gathering .................. 1.860 ( 1.9840 %)
Message Passage: data scattering.................. 2.760 ( 2.9440 %)
Total: 25.120 26.7947

All percentages are with respect to total time = 93.750


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

Analytical header files used:

ROMS/Functionals/ana_btflux.h

ROMS/TOMS: DONE... Monday - June 26, 2017 - 2:28:23 AM

[kate]

This is the next step. Learn to view the ROMS output.

[shiyf]

This is the next step. Learn to view the ROMS output.

Thank you very much kate!
The problem has been solved which is caused by steep depth in some special regions(see figure below),

I was even more smooth the deep, and the model runs smoothly!