my roms stop,but can not find error message

[shiyf]

hi,

I have run roms successfully first time, then i make nrrec = -1,and use ocean_rst file instead of inifile,
it stop,and do not output data,when i check the log message, i can not find any error!
I am very confused,Can you give me some advice?

log file:
Model Input Parameters: ROMS/TOMS version 3.6
Tuesday - February 25, 2014 - 9:03:29 PM
-----------------------------------------------------------------------------

SHIYF

Operating system : Linux
CPU/hardware : x86_64
Compiler system : ifort
Compiler command : /lustre/inspur/mpi/openmpi-intel/bin/mpif90
Compiler flags : -heap-arrays -fp-model precise -ip -O3 -free -free -free

Input Script : ocean_shiyf.in

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

Local Root : /lustre/home/shiyongfangWKJ/roms/trunk
Header Dir : /lustre/home/shiyongfangWKJ/roms/projects/roms_shiyf
Header file : shiyf.h
Analytical Dir: /lustre/home/shiyongfangWKJ/roms/projects/roms_shiyf

Resolution, Grid 01: 0459x0235x016, Parallel Nodes: 40, Tiling: 004x010


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

518400 ntimes Number of timesteps for 3-D equations.
60.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.
-1 nrrec Number of restart records to read from disk.
T LcycleRST Switch to recycle time-records in restart file.
43200 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).
43200 nHIS Number of timesteps between the writing fields
into history file.
518400 ndefHIS Number of timesteps between creation of new
history files.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
43200 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
518400 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.
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.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
2.0000E-02 Zos Surface roughness (m).
1 lmd_Jwt Jerlov water type.
1 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
7.0000E+00 theta_s S-coordinate surface control parameter.
2.0000E-01 theta_b S-coordinate bottom control parameter.
30.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)
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.
F 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.
F VolCons(4) NLM northern edge boundary volume conservation.
14.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).
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: roms_grid.nc
Input Nonlinear Initial File: ocean_rst.nc
Input Forcing File 01: roms_frc.nc
Input Boundary File: roms_bry.nc

Tile partition information for Grid 01: 0459x0235x0016 tiling: 004x010

tile Istr Iend Jstr Jend Npts

0 1 115 1 22 40480
1 116 230 1 22 40480
2 231 345 1 22 40480
3 346 459 1 22 40128
4 1 115 23 46 44160
5 116 230 23 46 44160
6 231 345 23 46 44160
7 346 459 23 46 43776
8 1 115 47 70 44160
9 116 230 47 70 44160
10 231 345 47 70 44160
11 346 459 47 70 43776
12 1 115 71 94 44160
13 116 230 71 94 44160
14 231 345 71 94 44160
15 346 459 71 94 43776
16 1 115 95 118 44160
17 116 230 95 118 44160
18 231 345 95 118 44160
19 346 459 95 118 43776
20 1 115 119 142 44160
21 116 230 119 142 44160
22 231 345 119 142 44160
23 346 459 119 142 43776
24 1 115 143 166 44160
25 116 230 143 166 44160
26 231 345 143 166 44160
27 346 459 143 166 43776
28 1 115 167 190 44160
29 116 230 167 190 44160
30 231 345 167 190 44160
31 346 459 167 190 43776
32 1 115 191 214 44160
33 116 230 191 214 44160
34 231 345 191 214 44160
35 346 459 191 214 43776
36 1 115 215 235 38640
37 116 230 215 235 38640
38 231 345 215 235 38640
39 346 459 215 235 38304

Tile minimum and maximum fractional grid coordinates:
(interior points only)

tile Xmin Xmax Ymin Ymax grid

0 0.50 115.50 0.50 22.50 RHO-points
1 115.50 230.50 0.50 22.50 RHO-points
2 230.50 345.50 0.50 22.50 RHO-points
3 345.50 459.50 0.50 22.50 RHO-points
4 0.50 115.50 22.50 46.50 RHO-points
5 115.50 230.50 22.50 46.50 RHO-points
6 230.50 345.50 22.50 46.50 RHO-points
7 345.50 459.50 22.50 46.50 RHO-points
8 0.50 115.50 46.50 70.50 RHO-points
9 115.50 230.50 46.50 70.50 RHO-points
10 230.50 345.50 46.50 70.50 RHO-points
11 345.50 459.50 46.50 70.50 RHO-points
12 0.50 115.50 70.50 94.50 RHO-points
13 115.50 230.50 70.50 94.50 RHO-points
14 230.50 345.50 70.50 94.50 RHO-points
15 345.50 459.50 70.50 94.50 RHO-points
16 0.50 115.50 94.50 118.50 RHO-points
17 115.50 230.50 94.50 118.50 RHO-points
18 230.50 345.50 94.50 118.50 RHO-points
19 345.50 459.50 94.50 118.50 RHO-points
20 0.50 115.50 118.50 142.50 RHO-points
21 115.50 230.50 118.50 142.50 RHO-points
22 230.50 345.50 118.50 142.50 RHO-points
23 345.50 459.50 118.50 142.50 RHO-points
24 0.50 115.50 142.50 166.50 RHO-points
25 115.50 230.50 142.50 166.50 RHO-points
26 230.50 345.50 142.50 166.50 RHO-points
27 345.50 459.50 142.50 166.50 RHO-points
28 0.50 115.50 166.50 190.50 RHO-points
29 115.50 230.50 166.50 190.50 RHO-points
30 230.50 345.50 166.50 190.50 RHO-points
31 345.50 459.50 166.50 190.50 RHO-points
32 0.50 115.50 190.50 214.50 RHO-points
33 115.50 230.50 190.50 214.50 RHO-points
34 230.50 345.50 190.50 214.50 RHO-points
35 345.50 459.50 190.50 214.50 RHO-points
36 0.50 115.50 214.50 235.50 RHO-points
37 115.50 230.50 214.50 235.50 RHO-points
38 230.50 345.50 214.50 235.50 RHO-points
39 345.50 459.50 214.50 235.50 RHO-points

0 1.00 115.50 0.50 22.50 U-points
1 115.50 230.50 0.50 22.50 U-points
2 230.50 345.50 0.50 22.50 U-points
3 345.50 459.00 0.50 22.50 U-points
4 1.00 115.50 22.50 46.50 U-points
5 115.50 230.50 22.50 46.50 U-points
6 230.50 345.50 22.50 46.50 U-points
7 345.50 459.00 22.50 46.50 U-points
8 1.00 115.50 46.50 70.50 U-points
9 115.50 230.50 46.50 70.50 U-points
10 230.50 345.50 46.50 70.50 U-points
11 345.50 459.00 46.50 70.50 U-points
12 1.00 115.50 70.50 94.50 U-points
13 115.50 230.50 70.50 94.50 U-points
14 230.50 345.50 70.50 94.50 U-points
15 345.50 459.00 70.50 94.50 U-points
16 1.00 115.50 94.50 118.50 U-points
17 115.50 230.50 94.50 118.50 U-points
18 230.50 345.50 94.50 118.50 U-points
19 345.50 459.00 94.50 118.50 U-points
20 1.00 115.50 118.50 142.50 U-points
21 115.50 230.50 118.50 142.50 U-points
22 230.50 345.50 118.50 142.50 U-points
23 345.50 459.00 118.50 142.50 U-points
24 1.00 115.50 142.50 166.50 U-points
25 115.50 230.50 142.50 166.50 U-points
26 230.50 345.50 142.50 166.50 U-points
27 345.50 459.00 142.50 166.50 U-points
28 1.00 115.50 166.50 190.50 U-points
29 115.50 230.50 166.50 190.50 U-points
30 230.50 345.50 166.50 190.50 U-points
31 345.50 459.00 166.50 190.50 U-points
32 1.00 115.50 190.50 214.50 U-points
33 115.50 230.50 190.50 214.50 U-points
34 230.50 345.50 190.50 214.50 U-points
35 345.50 459.00 190.50 214.50 U-points
36 1.00 115.50 214.50 235.50 U-points
37 115.50 230.50 214.50 235.50 U-points
38 230.50 345.50 214.50 235.50 U-points
39 345.50 459.00 214.50 235.50 U-points

0 0.50 115.50 1.00 22.50 V-points
1 115.50 230.50 1.00 22.50 V-points
2 230.50 345.50 1.00 22.50 V-points
3 345.50 459.50 1.00 22.50 V-points
4 0.50 115.50 22.50 46.50 V-points
5 115.50 230.50 22.50 46.50 V-points
6 230.50 345.50 22.50 46.50 V-points
7 345.50 459.50 22.50 46.50 V-points
8 0.50 115.50 46.50 70.50 V-points
9 115.50 230.50 46.50 70.50 V-points
10 230.50 345.50 46.50 70.50 V-points
11 345.50 459.50 46.50 70.50 V-points
12 0.50 115.50 70.50 94.50 V-points
13 115.50 230.50 70.50 94.50 V-points
14 230.50 345.50 70.50 94.50 V-points
15 345.50 459.50 70.50 94.50 V-points
16 0.50 115.50 94.50 118.50 V-points
17 115.50 230.50 94.50 118.50 V-points
18 230.50 345.50 94.50 118.50 V-points
19 345.50 459.50 94.50 118.50 V-points
20 0.50 115.50 118.50 142.50 V-points
21 115.50 230.50 118.50 142.50 V-points
22 230.50 345.50 118.50 142.50 V-points
23 345.50 459.50 118.50 142.50 V-points
24 0.50 115.50 142.50 166.50 V-points
25 115.50 230.50 142.50 166.50 V-points
26 230.50 345.50 142.50 166.50 V-points
27 345.50 459.50 142.50 166.50 V-points
28 0.50 115.50 166.50 190.50 V-points
29 115.50 230.50 166.50 190.50 V-points
30 230.50 345.50 166.50 190.50 V-points
31 345.50 459.50 166.50 190.50 V-points
32 0.50 115.50 190.50 214.50 V-points
33 115.50 230.50 190.50 214.50 V-points
34 230.50 345.50 190.50 214.50 V-points
35 345.50 459.50 190.50 214.50 V-points
36 0.50 115.50 214.50 235.00 V-points
37 115.50 230.50 214.50 235.00 V-points
38 230.50 345.50 214.50 235.00 V-points
39 345.50 459.50 214.50 235.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 254
HaloSizeJ(1) = 72
TileSide(1) = 121
TileSize(1) = 3630


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.
RST_SINGLE Single precision fields in restart NetCDF file.
SALINITY Using salinity.
SCORRECTION Surface salinity flux correction.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
SPLINES Conservative parabolic spline reconstruction.
TS_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
UV_ADV Advection of momentum.
UV_COR Coriolis term.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_C4VADVECTION Fourth-order centered vertical advection of momentum.
UV_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 # 4 (pid= 9874) is active.
Node # 8 (pid= 9878) is active.
Node # 0 (pid= 9870) is active.
Node # 9 (pid= 9879) is active.
Node # 7 (pid= 9877) is active.
Node # 1 (pid= 9871) is active.
Node # 39 (pid= 4196) is active.
Node # 31 (pid= 4176) is active.
Node # 32 (pid= 4189) is active.
Node # 33 (pid= 4190) is active.
Node # 36 (pid= 4193) is active.
Node # 14 (pid= 21829) is active.
Node # 16 (pid= 21831) is active.
Node # 10 (pid= 21825) is active.
Node # 11 (pid= 21826) is active.
Node # 24 (pid= 9146) is active.
Node # 12 (pid= 21827) is active.
Node # 25 (pid= 9147) is active.
Node # 13 (pid= 21828) is active.
Node # 28 (pid= 9150) is active.
Node # 15 (pid= 21830) is active.
Node # 23 (pid= 9145) is active.
Node # 17 (pid= 21832) is active.
Node # 20 (pid= 9142) is active.
Node # 6 (pid= 9876) is active.
Node # 22 (pid= 9144) is active.
Node # 30 (pid= 4175) is active.
Node # 38 (pid= 4195) is active.
Node # 5 (pid= 9875) is active.
Node # 37 (pid= 4194) is active.
Node # 21 (pid= 9143) is active.
Node # 29 (pid= 9151) is active.
Node # 2 (pid= 9872) is active.
Node # 26 (pid= 9148) is active.
Node # 34 (pid= 4191) is active.
Node # 18 (pid= 21833) is active.
Node # 3 (pid= 9873) is active.
Node # 27 (pid= 9149) is active.
Node # 35 (pid= 4192) is active.
Node # 19 (pid= 21834) is active.

INITIAL: Configuring and initializing forward nonlinear model ...


Vertical S-coordinate System:

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

16 0.0000000 0.0000000 0.000 0.000 0.000 0.000
15 -0.0625000 -0.0009137 -2.030 -1.875 -4.223 -6.416
14 -0.1250000 -0.0023093 -4.143 -3.750 -9.685 -15.227
13 -0.1875000 -0.0048224 -6.445 -5.625 -18.019 -29.592
12 -0.2500000 -0.0097619 -9.160 -7.500 -32.588 -56.016
11 -0.3125000 -0.0197716 -12.736 -9.375 -60.188 -107.640
10 -0.3750000 -0.0394986 -17.965 -11.250 -112.761 -207.558
9 -0.4375000 -0.0743313 -25.761 -13.125 -204.156 -382.551
8 -0.5000000 -0.1241359 -36.103 -15.000 -334.029 -631.955
7 -0.5625000 -0.1786344 -47.243 -16.875 -475.966 -904.688
6 -0.6250000 -0.2284617 -57.588 -18.750 -605.897 -1154.205
5 -0.6875000 -0.2764002 -67.613 -20.625 -730.974 -1394.334
4 -0.7500000 -0.3333248 -79.165 -22.500 -879.145 -1679.124
3 -0.8125000 -0.4130065 -94.586 -24.375 -1085.802 -2077.017
2 -0.8750000 -0.5326249 -116.796 -26.250 -1395.096 -2673.396
1 -0.9375000 -0.7162634 -149.890 -28.125 -1868.922 -3587.954
0 -1.0000000 -1.0000000 -200.000 -30.000 -2600.000 -5000.000

Time Splitting Weights: 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

Minimum X-grid spacing, DXmin = 2.33372303E+01 km
Maximum X-grid spacing, DXmax = 5.55664199E+01 km
Minimum Y-grid spacing, DYmin = 2.36170777E+01 km
Maximum Y-grid spacing, DYmax = 5.55654002E+01 km
Minimum Z-grid spacing, DZmin = 2.03032458E+00 m
Maximum Z-grid spacing, DZmax = 1.41204579E+03 m

Minimum barotropic Courant Number = 4.50937506E-03
Maximum barotropic Courant Number = 3.91802145E-02
Maximum Coriolis Courant Number = 7.91968582E-03

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

NLM: GET_STATE - Read state initial conditions, t = 360 00:00:00
(File: ocean_rst.nc, Rec=0002, Index=1)
- free-surface
(Min = -9.21388030E-01 Max = 9.34340298E-01)
- vertically integrated u-momentum component
(Min = -4.59825516E-01 Max = 8.06407392E-01)
- vertically integrated v-momentum component
(Min = -7.90307343E-01 Max = 7.75139153E-01)
- u-momentum component
(Min = -1.26890314E+00 Max = 1.62352765E+00)
- v-momentum component
(Min = -1.65238345E+00 Max = 1.39365602E+00)
- potential temperature
(Min = -1.82950630E+01 Max = 3.57503929E+01)
- salinity
(Min = 1.71523914E+01 Max = 4.05895386E+01)
- vertical viscosity coefficient
(Min = 0.00000000E+00 Max = 4.28411903E+01)
- temperature vertical diffusion coefficient
(Min = 0.00000000E+00 Max = 5.81563568E+01)
- salinity vertical diffusion coefficient
(Min = 0.00000000E+00 Max = 5.81563568E+01)
GET_2DFLD - surface u-momentum stress, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.68432195E-04 Max = 1.68611613E-04)
GET_2DFLD - surface v-momentum stress, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.30955021E-04 Max = 1.16607114E-04)
GET_2DFLD - solar shortwave radiation flux, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 7.14978554E-08 Max = 8.99074356E-05)
GET_2DFLD - surface net heat flux, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.22467813E-04 Max = 4.19740098E-05)
GET_2DFLD - surface net freswater flux, (E-P), t = 345 00:00:00
(Rec=0012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.01856720E-07 Max = 1.76678856E-07)
GET_2DFLD - sea surface salinity climatology, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.65167007E+01 Max = 4.11551018E+01)
GET_NGFLD - free-surface eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.73167685E-01 Max = 0.00000000E+00)
GET_NGFLD - free-surface southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.58965313E-01 Max = 7.52205844E-02)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.42503950E-02 Max = 3.94072602E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -5.20802142E-03 Max = 4.08989435E-03)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.78351520E-02 Max = 4.50582498E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.93825337E-02 Max = 3.05695320E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.85696218E-01 Max = 2.46939123E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -8.50273120E-02 Max = 5.58325617E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.58550369E-02 Max = 6.04781091E-03)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -8.92719773E-02 Max = 3.66097382E-02)
GET_NGFLD - temperature eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.76033770E+00 Max = 2.89773041E+01)
GET_NGFLD - salinity eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.33907192E+01 Max = 3.59697255E+01)
GET_NGFLD - temperature southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 9.74252701E-01 Max = 2.43833026E+01)
GET_NGFLD - salinity southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.42717653E+01 Max = 3.59174246E+01)

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


Initial basin volumes: TotVolume = 5.9777232065E+17 m3
MinVolume = 1.1279533624E+09 m3
MaxVolume = 4.3597847185E+12 m3
Max/Min = 3.8652171834E+03


NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00518401 - 00518400)

GET_2DFLD - surface u-momentum stress, t = 15 00:00:00
(Rec=0001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.40968276E-04 Max = 1.81788799E-04)
GET_2DFLD - surface v-momentum stress, t = 15 00:00:00
(Rec=0001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.30423325E-04 Max = 1.34807913E-04)
GET_2DFLD - solar shortwave radiation flux, t = 15 00:00:00
(Rec=0001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.81343152E-07 Max = 8.01640752E-05)
GET_2DFLD - surface net heat flux, t = 15 00:00:00
(Rec=0001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.28903545E-04 Max = 3.81661335E-05)
GET_2DFLD - surface net freswater flux, (E-P), t = 15 00:00:00
(Rec=0001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.56888076E-07 Max = 1.34596601E-07)
GET_2DFLD - sea surface salinity climatology, t = 15 00:00:00
(Rec=0001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.77527008E+01 Max = 4.05689011E+01)
GET_NGFLD - free-surface eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.50576890E-01 Max = 0.00000000E+00)
GET_NGFLD - free-surface southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.25712829E-01 Max = 9.00352100E-02)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -5.39048305E-03 Max = 3.10402664E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.49644176E-02 Max = 5.50461100E-03)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -6.44210545E-03 Max = 3.25296741E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.90746580E-02 Max = 4.79667278E-03)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.40916645E-01 Max = 2.03545932E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -6.73175687E-02 Max = 6.65828306E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -4.22836036E-02 Max = 2.35178151E-03)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -6.42212231E-02 Max = 2.90233429E-02)
GET_NGFLD - temperature eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.76033770E+00 Max = 2.83225178E+01)
GET_NGFLD - salinity eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.36186304E+01 Max = 3.60024570E+01)
GET_NGFLD - temperature southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 9.74252701E-01 Max = 2.58425064E+01)
GET_NGFLD - salinity southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.42761201E+01 Max = 3.58729397E+01)

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

518400 360 00:00:00 1.925411E-03 2.095307E+04 2.095307E+04 6.117299E+17
(149,064,16) 8.596160E-04 1.439332E-03 0.000000E+00 1.630741E+00

Elapsed CPU time (seconds):

Node # 1 CPU: 2.740
Node # 5 CPU: 2.623
Node # 8 CPU: 2.632
Node # 0 CPU: 2.244
Node # 4 CPU: 2.744
Node # 2 CPU: 2.753
Node # 6 CPU: 2.634
Node # 9 CPU: 2.341
Node # 3 CPU: 2.748
Node # 10 CPU: 2.208
Node # 24 CPU: 2.542
Node # 37 CPU: 2.558
Node # 13 CPU: 2.543
Node # 26 CPU: 2.563
Node # 34 CPU: 2.548
Node # 17 CPU: 2.526
Node # 21 CPU: 2.545
Node # 33 CPU: 2.538
Node # 30 CPU: 2.272
Node # 36 CPU: 2.517
Node # 39 CPU: 2.556

ROMS/TOMS - Output NetCDF summary for Grid 01:

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_hmixcoef.h
ROMS/Functionals/ana_nudgcoef.h
Node # 32 CPU: 2.531
Node # 16 CPU: 2.530
Node # 22 CPU: 2.552
Node # 35 CPU: 2.562
Node # 12 CPU: 2.543
Node # 20 CPU: 2.206
Node # 38 CPU: 2.562
Node # 14 CPU: 2.552

ROMS/TOMS: DONE... Tuesday - February 25, 2014 - 9:03:32 PM
Node # 28 CPU: 2.536
Node # 31 CPU: 2.529
Node # 11 CPU: 2.538
Node # 27 CPU: 2.557
Node # 19 CPU: 2.202
Node # 29 CPU: 2.217
Node # 15 CPU: 2.534
Node # 23 CPU: 2.546
Node # 18 CPU: 2.534
Node # 25 CPU: 2.558
Node # 7 CPU: 2.624

[kate]

You have told it to run some NTIMES steps from the beginning. It knows it already did that and doesn't need to run longer. I suggest you add to NTIMES for more timesteps.

[shiyf]

You have told it to run some NTIMES steps from the beginning. It knows it already did that and doesn't need to run longer. I suggest you add to NTIMES for more timesteps.

thanks a lot! I do not know if I understand correctly .when i run the first time making nrrec =0 for one year, NTIMES steps * dt equals one year, then i make nrrec = -1,and use ocean_rst file instead of inifile,
is the model starting from the second year? and my NTIMES steps*dt should equal The length of the two years? the output data continue to write in his avg file from 13 months to 24 months,and so on?

[kate]

Yes, that's right, try it out.

[shiyf]

Yes, that's right, try it out.

Thanks a lot,it is ok!