changed some water grid into land grid,the model blow up

[fanzhongya]

hello everyone

I had a problem.Please help me.THANKS!

I run a real case in Changjiang Estuary.The model can run well in pictureI and pictureII.However,only changed some water grid into land grid as described in pictureIII,the model blow up abnormally.

the record file as followed:

Process Information:

Node # 0 (pid= 7529) is active.
Node # 2 (pid= 7524) is active.

Model Input Parameters: ROMS/TOMS version 3.2
Thursday - June 24, 2010 - 11:01:58 AM
-----------------------------------------------------------------------------
Node # 5 (pid= 7526) is active.
Node # 3 (pid= 7523) is active.
Node # 14 (pid= 4698) is active.
Node # 30 (pid= 28175) is active.
Node # 12 (pid= 4695) is active.
Node # 26 (pid= 28173) is active.
Node # 25 (pid= 28174) is active.
Node # 7 (pid= 7528) is active.
Node # 16 (pid= 6493) is active.
Node # 8 (pid= 4693) is active.
Node # 10 (pid= 4692) is active.
Node # 22 (pid= 6499) is active.
Node # 9 (pid= 4691) is active.
Node # 1 (pid= 7522) is active.
Node # 19 (pid= 6495) is active.
Node # 20 (pid= 6497) is active.
Node # 24 (pid= 28170) is active.
Node # 23 (pid= 6500) is active.
Node # 11 (pid= 4694) is active.
Node # 21 (pid= 6498) is active.
Node # 18 (pid= 6496) is active.
Node # 15 (pid= 4696) is active.
Node # 31 (pid= 28177) is active.
Node # 17 (pid= 6494) is active.
Node # 4 (pid= 7525) is active.
Node # 13 (pid= 4697) is active.
Node # 28 (pid= 28172) is active.
Node # 27 (pid= 28171) is active.
Node # 29 (pid= 28176) is active.
Node # 6 (pid= 7527) is active.

Yangtz River Estuary

Operating system : Linux
CPU/hardware : x86_64
Compiler system : ifort
Compiler command : /usr/mpi/intel/mvapich2-1.2p1/bin/mpif90
Compiler flags : -ip -O3 -xW -free

Input Script : ./cjk_in/cjk0402.in

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

Local Root : /export/home/pxding/fan/cjk
Header Dir : ./ROMS/Include
Header file : cjk.h
Analytical Dir: /export/home/pxding/fan/cjk/ROMS/Functionals

Resolution, Grid 01: 0270x0232x010, Parallel Nodes: 32, Tiling: 004x008


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

57600 ntimes Number of timesteps for 3-D equations.
60.000 dt Timestep size (s) for 3-D equations.
5 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.
129599 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).
60 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.
1 ntsDIA Starting timestep for the accumulation of output
time-averaged diagnostics data.
14400 nDIA Number of timesteps between the writing of
time-averaged data into diagnostics file.
1.0000E+01 tnu2(01) Horizontal, harmonic mixing coefficient (m2/s)
for tracer 01: temp
1.0000E+01 tnu2(02) Horizontal, harmonic mixing coefficient (m2/s)
for tracer 02: salt
0.0000E+00 visc2 Horizontal, harmonic mixing coefficient (m2/s)
for momentum.
2.5000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
2.5000E-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.
0.000 gls_p GLS stability exponent.
1.000 gls_m GLS turbulent kinetic energy exponent.
1.000 gls_n GLS turbulent length scale exponent.
5.0000E-06 gls_Kmin GLS minimum value of turbulent kinetic energy.
5.0000E-06 gls_Pmin GLS minimum value of dissipation.
5.5440E-01 gls_cmu0 GLS stability coefficient.
9.0000E-01 gls_c1 GLS shear production coefficient.
5.2000E-01 gls_c2 GLS dissipation coefficient.
2.5000E+00 gls_c3m GLS stable buoyancy production coefficient.
1.0000E+00 gls_c3p GLS unstable buoyancy production coefficient.
1.9600E+00 gls_sigk GLS constant Schmidt number for TKE.
1.9600E+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.
0.0000E+00 rdrg Linear bottom drag coefficient (m/s).
0.0000E+00 rdrg2 Quadratic bottom drag coefficient.
3.0000E-04 Zob Bottom roughness (m).
2.0000E-02 Zos Surface roughness (m).
5.0000E-01 Dcrit Minimum depth for wetting and drying (m).
1 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
6.0000E-06 theta_s S-coordinate surface control parameter.
1.0000E+00 theta_b S-coordinate bottom control parameter.
0.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1024.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
0.000 dstart Time-stamp assigned to model initialization (days).
0.000 tide_start Reference time origin for tidal forcing (days).
20040120.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.
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).
T Hout(idFsur) Write out free-surface.
T Hout(idUvel) Write out 3D U-momentum component.
T Hout(idVvel) Write out 3D V-momentum component.
T Hout(idTvar) Write out tracer 01: temp
T Hout(idTvar) Write out tracer 02: salt

Output/Input Files:

Output Restart File: cjk_in/Output/cjk_rst.nc
Output History File: cjk_in/Output/yz_0402_4_his.nc
Output Averages File: cjk_in/Output/cjk_avg.nc
Output Diagnostics File: cjk_in/Output/cjk_dia.nc
Input Grid File: cjk_in/Input/yz_3.2_4.nc
Input Nonlinear Initial File: cjk_in/Input/cjk_0402_4_ini.nc
Input Forcing File 01: cjk_in/Input/cjk_0402_tides.nc
Input Forcing File 02: cjk_in/Input/yz0402_rvr.nc
Input Forcing File 03: cjk_in/Input/cjk_0402_smf.nc

Tile partition information for Grid 01: 0270x0232x0010 tiling: 004x008

tile Istr Iend Jstr Jend Npts

0 1 67 1 29 19430
1 68 135 1 29 19720
2 136 203 1 29 19720
3 204 270 1 29 19430
4 1 67 30 58 19430
5 68 135 30 58 19720
6 136 203 30 58 19720
7 204 270 30 58 19430
8 1 67 59 87 19430
9 68 135 59 87 19720
10 136 203 59 87 19720
11 204 270 59 87 19430
12 1 67 88 116 19430
13 68 135 88 116 19720
14 136 203 88 116 19720
15 204 270 88 116 19430
16 1 67 117 145 19430
17 68 135 117 145 19720
18 136 203 117 145 19720
19 204 270 117 145 19430
20 1 67 146 174 19430
21 68 135 146 174 19720
22 136 203 146 174 19720
23 204 270 146 174 19430
24 1 67 175 203 19430
25 68 135 175 203 19720
26 136 203 175 203 19720
27 204 270 175 203 19430
28 1 67 204 232 19430
29 68 135 204 232 19720
30 136 203 204 232 19720
31 204 270 204 232 19430

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 67.50 0.50 29.50 RHO-points
1 67.50 135.50 0.50 29.50 RHO-points
2 135.50 203.50 0.50 29.50 RHO-points
3 203.50 270.50 0.50 29.50 RHO-points
4 0.50 67.50 29.50 58.50 RHO-points
5 67.50 135.50 29.50 58.50 RHO-points
6 135.50 203.50 29.50 58.50 RHO-points
7 203.50 270.50 29.50 58.50 RHO-points
8 0.50 67.50 58.50 87.50 RHO-points
9 67.50 135.50 58.50 87.50 RHO-points
10 135.50 203.50 58.50 87.50 RHO-points
11 203.50 270.50 58.50 87.50 RHO-points
12 0.50 67.50 87.50 116.50 RHO-points
13 67.50 135.50 87.50 116.50 RHO-points
14 135.50 203.50 87.50 116.50 RHO-points
15 203.50 270.50 87.50 116.50 RHO-points
16 0.50 67.50 116.50 145.50 RHO-points
17 67.50 135.50 116.50 145.50 RHO-points
18 135.50 203.50 116.50 145.50 RHO-points
19 203.50 270.50 116.50 145.50 RHO-points
20 0.50 67.50 145.50 174.50 RHO-points
21 67.50 135.50 145.50 174.50 RHO-points
22 135.50 203.50 145.50 174.50 RHO-points
23 203.50 270.50 145.50 174.50 RHO-points
24 0.50 67.50 174.50 203.50 RHO-points
25 67.50 135.50 174.50 203.50 RHO-points
26 135.50 203.50 174.50 203.50 RHO-points
27 203.50 270.50 174.50 203.50 RHO-points
28 0.50 67.50 203.50 232.50 RHO-points
29 67.50 135.50 203.50 232.50 RHO-points
30 135.50 203.50 203.50 232.50 RHO-points
31 203.50 270.50 203.50 232.50 RHO-points

0 1.00 67.50 0.50 29.50 U-points
1 67.50 135.50 0.50 29.50 U-points
2 135.50 203.50 0.50 29.50 U-points
3 203.50 270.00 0.50 29.50 U-points
4 1.00 67.50 29.50 58.50 U-points
5 67.50 135.50 29.50 58.50 U-points
6 135.50 203.50 29.50 58.50 U-points
7 203.50 270.00 29.50 58.50 U-points
8 1.00 67.50 58.50 87.50 U-points
9 67.50 135.50 58.50 87.50 U-points
10 135.50 203.50 58.50 87.50 U-points
11 203.50 270.00 58.50 87.50 U-points
12 1.00 67.50 87.50 116.50 U-points
13 67.50 135.50 87.50 116.50 U-points
14 135.50 203.50 87.50 116.50 U-points
15 203.50 270.00 87.50 116.50 U-points
16 1.00 67.50 116.50 145.50 U-points
17 67.50 135.50 116.50 145.50 U-points
18 135.50 203.50 116.50 145.50 U-points
19 203.50 270.00 116.50 145.50 U-points
20 1.00 67.50 145.50 174.50 U-points
21 67.50 135.50 145.50 174.50 U-points
22 135.50 203.50 145.50 174.50 U-points
23 203.50 270.00 145.50 174.50 U-points
24 1.00 67.50 174.50 203.50 U-points
25 67.50 135.50 174.50 203.50 U-points
26 135.50 203.50 174.50 203.50 U-points
27 203.50 270.00 174.50 203.50 U-points
28 1.00 67.50 203.50 232.50 U-points
29 67.50 135.50 203.50 232.50 U-points
30 135.50 203.50 203.50 232.50 U-points
31 203.50 270.00 203.50 232.50 U-points

0 0.50 67.50 1.00 29.50 V-points
1 67.50 135.50 1.00 29.50 V-points
2 135.50 203.50 1.00 29.50 V-points
3 203.50 270.50 1.00 29.50 V-points
4 0.50 67.50 29.50 58.50 V-points
5 67.50 135.50 29.50 58.50 V-points
6 135.50 203.50 29.50 58.50 V-points
7 203.50 270.50 29.50 58.50 V-points
8 0.50 67.50 58.50 87.50 V-points
9 67.50 135.50 58.50 87.50 V-points
10 135.50 203.50 58.50 87.50 V-points
11 203.50 270.50 58.50 87.50 V-points
12 0.50 67.50 87.50 116.50 V-points
13 67.50 135.50 87.50 116.50 V-points
14 135.50 203.50 87.50 116.50 V-points
15 203.50 270.50 87.50 116.50 V-points
16 0.50 67.50 116.50 145.50 V-points
17 67.50 135.50 116.50 145.50 V-points
18 135.50 203.50 116.50 145.50 V-points
19 203.50 270.50 116.50 145.50 V-points
20 0.50 67.50 145.50 174.50 V-points
21 67.50 135.50 145.50 174.50 V-points
22 135.50 203.50 145.50 174.50 V-points
23 203.50 270.50 145.50 174.50 V-points
24 0.50 67.50 174.50 203.50 V-points
25 67.50 135.50 174.50 203.50 V-points
26 135.50 203.50 174.50 203.50 V-points
27 203.50 270.50 174.50 203.50 V-points
28 0.50 67.50 203.50 232.00 V-points
29 67.50 135.50 203.50 232.00 V-points
30 135.50 203.50 203.50 232.00 V-points
31 203.50 270.50 203.50 232.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 246
HaloSizeJ(1) = 129
TileSide(1) = 76
TileSize(1) = 2812


Activated C-preprocessing Options:

cjk_calculate Yangtz River Estuary
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_FSOBC Analytical free-surface boundary conditions.
ANA_M2OBC Analytical 2D momentum boundary conditions.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ANA_STFLUX Analytical kinematic surface temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged fields.
CURVGRID Orthogonal curvilinear grid.
DIAGNOSTICS_UV Computing and writing momentum diagnostic terms.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EAST_FSCLAMPED Eastern edge, free-surface, Clamped condition.
EAST_M2REDUCED Eastern edge, 2D momentum, reduced-physics condition.
FSOBC_REDUCED Using free-surface data in reduced physics condtions
GLS_MIXING Generic Length-Scale turbulence closure.
INLINE_2DIO Processing 3D IO level by level to reduce memory needs.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MIX_S_TS Mixing of tracers along constant S-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.
NORTH_FSCLAMPED Northern edge, free-surface, Clamped condition.
NORTH_M2REDUCED Northern edge, 2D momentum, reduced-physics condition.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
OUT_DOUBLE Double precision output fields in NetCDF files.
PERFECT_RESTART Processing perfect restart variables.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
RAMP_TIDES Ramping tidal forcing for one day.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
SOUTH_FSCLAMPED Southern edge, free-surface, Clamped condition.
SOUTH_M2REDUCED Southern edge, 2D momentum, reduced-physics condition.
SPLINES Conservative parabolic spline reconstruction.
SSH_TIDES Add tidal elevation to SSH climatology.
THREE_GHOST Using three Ghost Points in halo regions.
TS_MPDATA Recursive flux corrected MPDATA 3D advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
TS_PSOURCE Tracers point sources and sinks.
UV_ADV Advection of momentum.
UV_COR Coriolis term.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_C4VADVECTION Fourth-order centered vertical advection of momentum.
UV_LOGDRAG Logarithmic bottom stress.
UV_PSOURCE Mass point sources and sinks.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
WESTERN_WALL Wall boundary at Western edge.
WET_DRY Wetting and drying activated.

INITIAL: Configuring and initializing forward nonlinear model ...


Vertical S-coordinate System:

level S-coord Cs-curve at_hmin over_slope at_hmax

10 0.0000000 0.0000000 0.000 0.000 0.000
9 -0.1000000 -0.1000000 -0.050 -4.037 -8.023
8 -0.2000000 -0.2000000 -0.100 -8.074 -16.047
7 -0.3000000 -0.3000000 -0.151 -12.110 -24.070
6 -0.4000000 -0.4000000 -0.201 -16.147 -32.094
5 -0.5000000 -0.5000000 -0.251 -20.184 -40.117
4 -0.6000000 -0.6000000 -0.301 -24.221 -48.140
3 -0.7000000 -0.7000000 -0.351 -28.258 -56.164
2 -0.8000000 -0.8000000 -0.402 -32.294 -64.187
1 -0.9000000 -0.9000000 -0.452 -36.331 -72.210
0 -1.0000000 -1.0000000 -0.502 -40.368 -80.234

Time Splitting Weights: ndtfast = 5 nfast = 7

Primary Secondary Accumulated to Current Step

1 0.0119306869089341 0.2000000000000000 0.0119306869089341 0.2000000000000000
2 0.0475049776008205 0.1976138626182131 0.0594356645097547 0.3976138626182131
3 0.1047629417620785 0.1881128670980491 0.1641986062718331 0.5857267297162622
4 0.1760826281731955 0.1671602787456334 0.3402812344450286 0.7528870084618955
5 0.2416469636634478 0.1319437531109943 0.5819281981084765 0.8848307615728898
6 0.2602974116474964 0.0836143603783047 0.8422256097559728 0.9684451219511945
7 0.1577743902440273 0.0315548780488054 1.0000000000000000 0.9999999999999999

ndtfast, nfast = 5 7 nfast/ndtfast = 1.40000

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

1.000000000000 1.084198606272 0.542099303136 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.00000

Minimum X-grid spacing, DXmin = 1.20649429E-02 km
Maximum X-grid spacing, DXmax = 1.26877789E+01 km
Minimum Y-grid spacing, DYmin = 3.82419589E-02 km
Maximum Y-grid spacing, DYmax = 1.04691316E+01 km
Minimum Z-grid spacing, DZmin = 5.02027180E-02 m
Maximum Z-grid spacing, DZmax = 8.02338420E+00 m

Minimum barotropic Courant Number = 1.53122431E-02
Maximum barotropic Courant Number = 7.93292092E+00
Maximum Coriolis Courant Number = 4.83202820E-03


NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(File: cjk_0402_4_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 = 1.50000000E+01 Max = 1.50000000E+01)
- salinity
(Min = 0.00000000E+00 Max = 3.40000000E+01)
GET_NGFLD - tidal period
(Min = 4.30820424E+04 Max = 9.67260960E+04)
GET_2DFLD - tidal elevation amplitude
(Min = 0.00000000E+00 Max = 2.67858389E+00)
GET_2DFLD - tidal elevation phase angle
(Min = 2.53042579E-06 Max = 6.28313850E+00)
GET_NGFLD - river runoff XI-positions at RHO-points
(Min = 1.00000000E+00 Max = 1.00000000E+00)
GET_NGFLD - river runoff ETA-positions at RHO-points
(Min = 1.15000000E+02 Max = 1.22000000E+02)
GET_NGFLD - river runoff direction
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - river runoff mass transport vertical profile
(Min = 2.00000000E-02 Max = 1.50000000E-01)
GET_NGFLD - river runoff flag
(Min = 2.00000000E+00 Max = 2.00000000E+00)
GET_NGFLD - river runoff mass transport, t = 0 00:00:00
(Rec=0001, Index=2, File: yz0402_rvr.nc)
(Tmin= 0.0000 Tmax= 59.5000)
(Min = 3.50000000E+03 Max = 3.50000000E+03)
GET_NGFLD - river runoff salinity, t = 0 00:00:00
(Rec=0001, Index=2, File: yz0402_rvr.nc)
(Tmin= 0.0000 Tmax= 59.5000)
(Min = 1.00000000E+00 Max = 1.00000000E+00)
GET_2DFLD - surface u-momentum stress, t = 0 00:00:00
(Rec=0001, Index=1, File: cjk_0402_smf.nc)
(Tmin= 0.0000 Tmax= 157.0000)
(Min = -8.10309795E-05 Max = 3.76189468E-05)
GET_2DFLD - surface v-momentum stress, t = 0 00:00:00
(Rec=0001, Index=1, File: cjk_0402_smf.nc)
(Tmin= 0.0000 Tmax= 157.0000)
(Min = 2.75237187E-05 Max = 2.15437673E-04)

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


Initial basin volumes: TotVolume = 6.38128429852600E+12 m3
MinVolume = 3.11297470446051E+03 m3
MaxVolume = 9.75363793296542E+08 m3
Max/Min = 3.13322106954151E+05

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00057600)
GET_NGFLD - river runoff mass transport, t = 0 12:00:00
(Rec=0002, Index=1, File: yz0402_rvr.nc)
(Tmin= 0.0000 Tmax= 59.5000)
(Min = 1.25000000E+03 Max = 1.25000000E+03)
GET_NGFLD - river runoff salinity, t = 0 12:00:00
(Rec=0002, Index=1, File: yz0402_rvr.nc)
(Tmin= 0.0000 Tmax= 59.5000)
(Min = 1.00000000E+00 Max = 1.00000000E+00)
GET_2DFLD - surface u-momentum stress, t = 1 00:00:00
(Rec=0002, Index=2, File: cjk_0402_smf.nc)
(Tmin= 0.0000 Tmax= 157.0000)
(Min = -1.03412595E-04 Max = 3.15749988E-05)
GET_2DFLD - surface v-momentum stress, t = 1 00:00:00
(Rec=0002, Index=2, File: cjk_0402_smf.nc)
(Tmin= 0.0000 Tmax= 157.0000)
(Min = 7.95569603E-06 Max = 1.50063234E-04)

STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME

0 0 00:00:00 0.000000E+00 2.359258E+02 2.359258E+02 6.402577E+12
DEF_HIS - creating history file: cjk_in/Output/yz_0402_4_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file: cjk_in/Output/cjk_avg.nc
DEF_DIAGS - creating diagnostics file: cjk_in/Output/cjk_dia.nc
1 0 00:01:00 4.761448+158 -Infinity -Infinity -3.213876E+60

Elapsed CPU time (seconds):

Node # 0 CPU: 1.438
Node # 1 CPU: 1.679
Node # 2 CPU: 1.682
Node # 6 CPU: 1.680
Node # 3 CPU: 1.679
Node # 13 CPU: 1.656
Node # 5 CPU: 1.681
Node # 11 CPU: 1.656
Node # 8 CPU: 1.656
Node # 7 CPU: 1.681
Node # 16 CPU: 1.679
Node # 14 CPU: 1.676
Node # 10 CPU: 1.668
Node # 9 CPU: 1.668
Node # 21 CPU: 1.678
Node # 22 CPU: 1.678
Node # 30 CPU: 1.680
Node # 29 CPU: 1.680
Node # 19 CPU: 1.675
Node # 17 CPU: 1.678
Node # 15 CPU: 1.672
Node # 31 CPU: 1.680
Node # 28 CPU: 1.681
Node # 12 CPU: 1.672
Node # 27 CPU: 1.678
Node # 24 CPU: 1.682
Node # 23 CPU: 1.677
Node # 25 CPU: 1.680
Node # 26 CPU: 1.680
Node # 18 CPU: 1.678
Node # 20 CPU: 1.675
Node # 4 CPU: 1.622
Total: 53.320

Nonlinear model elapsed time profile:

Initialization ................................... 23.142 (43.4018 %)
Reading of input data ............................ 7.677 (14.3985 %)
Processing of input data ......................... 0.165 ( 0.3094 %)
Computation of vertical boundary conditions ...... 0.017 ( 0.0319 %)
Computation of global information integrals ...... 0.151 ( 0.2832 %)
Writing of output data ........................... 13.461 (25.2452 %)
Model 2D kernel .................................. 2.402 ( 4.5044 %)
Tidal forcing .................................... 0.043 ( 0.0806 %)
2D/3D coupling, vertical metrics ................. 0.327 ( 0.6132 %)
Omega vertical velocity .......................... 0.257 ( 0.4820 %)
Equation of state for seawater ................... 0.512 ( 0.9601 %)
GLS vertical mixing parameterization ............. 1.297 ( 2.4322 %)
3D equations right-side terms .................... 0.339 ( 0.6357 %)
3D equations predictor step ...................... 0.601 ( 1.1271 %)
Pressure gradient ................................ 0.161 ( 0.3019 %)
Harmonic mixing of tracers, S-surfaces ........... 0.073 ( 0.1369 %)
Harmonic stress tensor, S-surfaces ............... 0.115 ( 0.2157 %)
Corrector time-step for 3D momentum .............. 0.825 ( 1.5471 %)
Corrector time-step for tracers .................. 1.312 ( 2.4604 %)
Total: 52.876 99.1673

Nonlinear model message Passage profile:

Message Passage: 2D halo exchanges ............... 1.528 ( 2.8655 %)
Message Passage: 3D halo exchanges ............... 1.473 ( 2.7623 %)
Message Passage: 4D halo exchanges ............... 0.631 ( 1.1833 %)
Message Passage: data broadcast .................. 26.207 (49.1497 %)
Message Passage: data reduction .................. 0.146 ( 0.2738 %)
Message Passage: data gathering .................. 3.609 ( 6.7680 %)
Message Passage: data scattering.................. 8.186 (15.3529 %)
Message Passage: boundary data gathering ......... 0.032 ( 0.0600 %)
Total: 41.811 78.4155

All percentages are with respect to total time = 53.320

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

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_fsobc.h
ROMS/Functionals/ana_m2obc.h
ROMS/Functionals/ana_stflux.h

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


MAIN: Abnormal termination: BLOWUP.

[kate]

That's quite a grid you have there!

As for why it's blowing up in the very first timestep, it's hard to say from here. I attack this sort of thing in the debugger. Having it go bad so very fast makes it easier to see what's going on - I can get it to blow up after running for a year. You know you have some infinities - where did they come from? Where did you divide by zero? You don't have zero depth anywhere, do you?