I meet with blowing up at the beginning of the model, can someone help me to find the problem?
Thanks a lot.
Physical Parameters, Grid: 01
=============================
1296000 ntimes Number of timesteps for 3-D equations.
1200.000 dt Timestep size (s) for 3-D equations.
120 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.
1 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).
72 nHIS Number of timesteps between the writing fields
into history file.
72 ndefHIS Number of timesteps between creation of new
history files.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
72 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
72 ndefAVG Number of timesteps between creation of new
time-averaged file.
2.0000E+01 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
2.0000E+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.
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.
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).
1 lmd_Jwt Jerlov water type.
1 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
3.0000E+00 theta_s S-coordinate surface control parameter.
3.0000E-01 theta_b S-coordinate bottom control parameter.
50.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)
9.0000E+01 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
9.0000E+01 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
9.0000E+01 Znudg Nudging/relaxation time scale (days)
for free-surface.
9.0000E+01 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
9.0000E+01 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
1.0000E+01 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).
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(idu3dE) Write out 3D U-wastward component at RHO-points.
T Hout(idv3dN) Write out 3D V-northward component at RHO-points.
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(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(idu3dE) Write out averaged 3D U-eastward at RHO-points.
T Aout(idv3dN) Write out averaged 3D V-northward at RHO-points.
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
T Aout(idHsbl) Write out averaged depth of surface boundary layer.
Output/Input Files:
Output Restart File: /public/userspace/ljd/IO_d10/IO_rst_d10.nc
Prefix for History Files: /public/userspace/ljd/IO_d10/IO_his_d10
Prefix for Averages Files: /public/userspace/ljd/IO_d10/IO_avg_d10
Input Grid File: IO_grd_d10.nc
Input Nonlinear Initial File: IO_ini_d10.nc
Input Forcing File 01: IO_frc_d10_01.nc
IO_frc_d10_02.nc
IO_frc_d10_03.nc
IO_frc_d10_04.nc
IO_frc_d10_05.nc
IO_frc_d10_06.nc
IO_frc_d10_07.nc
IO_frc_d10_08.nc
IO_frc_d10_09.nc
IO_frc_d10_10.nc
IO_frc_d10_11.nc
IO_frc_d10_12.nc
Input Climatology File: IO_clm_d10_TSUVZ.nc
Input Boundary File: IO_bry_d10.nc
Tile partition information for Grid 01: 0799x0629x0040 tiling: 016x012
tile Istr Iend Jstr Jend Npts
0 1 50 1 50 100000
1 51 100 1 50 100000
2 101 150 1 50 100000
3 151 200 1 50 100000
4 201 250 1 50 100000
5 251 300 1 50 100000
6 301 350 1 50 100000
..........
Maximum halo size in XI and ETA directions:
HaloSizeI(1) = 124
HaloSizeJ(1) = 130
TileSide(1) = 59
TileSize(1) = 3304
Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Closed Chapman Imp Chapman Imp Closed
ubar 1 Closed Flather Flather Closed
vbar 1 Closed Flather Flather Closed
u 1 Closed Rad + Nud Rad + Nud Closed
v 1 Closed Rad + Nud Rad + Nud Closed
temp 1 Closed Rad + Nud Rad + Nud Closed
salt 1 Closed Rad + Nud Rad + Nud Closed
Activated C-preprocessing Options:
IO_D10 Control run for the Indian Ocean, 1/10*1/10 deg
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.
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_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.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
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_QDRAG Quadratic bottom stress.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
Vertical S-coordinate System, Grid 01:
.........
level S-coord Cs-curve Z at hmin at hc half way at hmax
40 0.0000000 0.0000000 0.000 0.000 0.000 0.000
39 -0.0250000 -0.0076503 -1.250 -1.250 -26.133 -51.016
38 -0.0500000 -0.0156736 -2.500 -2.500 -53.479 -104.459
37 -0.0750000 -0.0241410 -3.750 -3.750 -82.270 -160.790
36 -0.1000000 -0.0331262 -5.000 -5.000 -112.745 -220.490
35 -0.1250000 -0.0427053 -6.250 -6.250 -145.151 -284.053
34 -0.1500000 -0.0529559 -7.500 -7.500 -179.742 -351.984
33 -0.1750000 -0.0639559 -8.750 -8.750 -216.770 -424.791
32 -0.2000000 -0.0757823 -10.000 -10.000 -256.486 -502.972
31 -0.2250000 -0.0885086 -11.250 -11.250 -299.129 -587.009
30 -0.2500000 -0.1022034 -12.500 -12.500 -344.922 -677.345
29 -0.2750000 -0.1169269 -13.750 -13.750 -394.061 -774.373
28 -0.3000000 -0.1327288 -15.000 -15.000 -446.708 -878.416
27 -0.3250000 -0.1496453 -16.250 -16.250 -502.980 -989.710
26 -0.3500000 -0.1676968 -17.500 -17.500 -562.943 -1108.387
25 -0.3750000 -0.1868861 -18.750 -18.750 -626.608 -1234.466
24 -0.4000000 -0.2071978 -20.000 -20.000 -693.923 -1367.845
23 -0.4250000 -0.2285984 -21.250 -21.250 -764.779 -1508.309
22 -0.4500000 -0.2510387 -22.500 -22.500 -839.018 -1655.535
21 -0.4750000 -0.2744569 -23.750 -23.750 -916.437 -1809.124
20 -0.5000000 -0.2987836 -25.000 -25.000 -996.811 -1968.621
19 -0.5250000 -0.3239476 -26.250 -26.250 -1079.908 -2133.566
18 -0.5500000 -0.3498824 -27.500 -27.500 -1165.512 -2303.525
17 -0.5750000 -0.3765328 -28.750 -28.750 -1253.444 -2478.138
16 -0.6000000 -0.4038607 -30.000 -30.000 -1343.580 -2657.159
15 -0.6250000 -0.4318501 -31.250 -31.250 -1435.867 -2840.484
14 -0.6500000 -0.4605110 -32.500 -32.500 -1530.338 -3028.176
13 -0.6750000 -0.4898810 -33.750 -33.750 -1627.116 -3220.482
12 -0.7000000 -0.5200268 -35.000 -35.000 -1726.417 -3417.833
11 -0.7250000 -0.5510432 -36.250 -36.250 -1828.549 -3620.848
10 -0.7500000 -0.5830519 -37.500 -37.500 -1933.909 -3830.319
9 -0.7750000 -0.6162000 -38.750 -38.750 -2042.975 -4047.201
8 -0.8000000 -0.6506573 -40.000 -40.000 -2156.300 -4272.599
7 -0.8250000 -0.6866143 -41.250 -41.250 -2274.502 -4507.754
6 -0.8500000 -0.7242808 -42.500 -42.500 -2398.264 -4754.029
5 -0.8750000 -0.7638835 -43.750 -43.750 -2528.324 -5012.899
4 -0.9000000 -0.8056655 -45.000 -45.000 -2665.473 -5285.946
3 -0.9250000 -0.8498858 -46.250 -46.250 -2810.552 -5574.854
2 -0.9500000 -0.8968188 -47.500 -47.500 -2964.454 -5881.408
1 -0.9750000 -0.9467546 -48.750 -48.750 -3128.123 -6207.496
0 -1.0000000 -1.0000000 -50.000 -50.000 -3302.557 -6555.114
Time Splitting Weights for Grid 01: ndtfast = 120 nfast = 162
==================================
..........
ndtfast, nfast = 120 162 nfast/ndtfast = 1.35000
Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):
1.000000000000 1.016453155986 0.508226577993 1.000000000000 1.000000000000
Power filter parameters, Fgamma, gamma = 0.28400 0.26033
Metrics information for Grid 01:
===============================
Minimum X-grid spacing, DXmin = 9.61983482E+00 km
Maximum X-grid spacing, DXmax = 1.11132834E+01 km
Minimum Y-grid spacing, DYmin = 9.62858029E+00 km
Maximum Y-grid spacing, DYmax = 1.11132783E+01 km
Minimum Z-grid spacing, DZmin = 1.25000000E+00 m
Maximum Z-grid spacing, DZmax = 3.47617667E+02 m
Minimum barotropic Courant Number = 2.81833239E-02
Maximum barotropic Courant Number = 3.42277905E-01
Maximum Coriolis Courant Number = 8.73901467E-02
NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: IO_ini_d10.nc, Rec=0001, Index=1)
- free-surface
(Min = -1.65856037E-01 Max = 1.09070785E+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 = -7.12925634E-01 Max = 7.55215659E-01)
- v-momentum component
(Min = -7.72852753E-01 Max = 7.65704231E-01)
- potential temperature
(Min = 5.25969674E-02 Max = 2.98913774E+01)
- salinity
(Min = 2.12743881E+00 Max = 4.07141128E+01)
GET_2DFLD - surface u-momentum stress, t = 30 00:00:00
(Rec=0000030, Index=1, File: IO_frc_d10_01.nc)
(Tmin= 1.0000 Tmax= 30.0000)
(Min = -3.34676881E-04 Max = 4.07706802E-04)
GET_2DFLD - surface v-momentum stress, t = 30 00:00:00
(Rec=0000030, Index=1, File: IO_frc_d10_01.nc)
(Tmin= 1.0000 Tmax= 30.0000)
(Min = -2.24989736E-04 Max = 2.88494194E-04)
GET_2DFLD - solar shortwave radiation flux, t = 30 00:00:00
(Rec=0000030, Index=1, File: IO_frc_d10_01.nc)
(Tmin= 1.0000 Tmax= 30.0000)
(Min = 2.57790513E-05 Max = 7.66878176E-05)
GET_2DFLD - surface net heat flux, t = 30 00:00:00
(Rec=0000030, Index=1, File: IO_frc_d10_01.nc)
(Tmin= 1.0000 Tmax= 30.0000)
(Min = -5.53345135E-05 Max = 3.85001502E-05)
GET_2DFLD - surface net freswater flux, (E-P), t = 30 00:00:00
(Rec=0000030, Index=1, File: IO_frc_d10_01.nc)
(Tmin= 1.0000 Tmax= 30.0000)
(Min = -1.88254621E-07 Max = 9.84993660E-08)
GET_NGFLD - free-surface eastern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 9.16177432E-02 Max = 8.75259395E-01)
GET_NGFLD - free-surface southern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 7.29231921E-02 Max = 7.42595725E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = -2.83704406E-01 Max = 7.63229013E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = -2.64802093E-01 Max = 6.51729572E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = -7.15777892E-01 Max = 1.10803898E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = -3.84782651E-01 Max = 9.38375501E-02)
GET_NGFLD - temperature eastern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 6.58606760E-01 Max = 2.96212242E+01)
GET_NGFLD - salinity eastern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 2.48901541E+01 Max = 3.57507704E+01)
GET_NGFLD - temperature southern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 4.04895103E-01 Max = 2.52432386E+01)
GET_NGFLD - salinity southern boundary condition, t = 356 00:00:00
(Rec=0000072, Index=2, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 3.43807778E+01 Max = 3.58496595E+01)
GET_3DFLD - potential temperature climatology, t = 271 00:00:00
(Rec=0000004, Index=1, File: IO_clm_d10_TSUVZ.nc)
(Tmin= 1.0000 Tmax= 271.0000)
(Min = 5.98539002E-02 Max = 3.09193401E+01)
GET_3DFLD - salinity climatology, t = 271 00:00:00
(Rec=0000004, Index=1, File: IO_clm_d10_TSUVZ.nc)
(Tmin= 1.0000 Tmax= 271.0000)
(Min = 2.12883019E+00 Max = 4.07080536E+01)
Basin information for Grid 01:
Maximum grid stiffness ratios: rx0 = 2.087711E-01 (Beckmann and Haidvogel)
rx1 = 9.239405E+00 (Haney)
Initial basin volumes: TotVolume = 1.3422006839E+17 m3
MinVolume = 1.1592918981E+08 m3
MaxVolume = 4.1624495376E+10 m3
Max/Min = 3.5905103317E+02
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 01296000)
GET_2DFLD - surface u-momentum stress, t = 1 00:00:00
(Rec=0000001, Index=2, File: IO_frc_d10_01.nc)
(Tmin= 1.0000 Tmax= 30.0000)
(Min = -3.95970791E-04 Max = 3.27253677E-04)
GET_2DFLD - surface v-momentum stress, t = 1 00:00:00
(Rec=0000001, Index=2, File: IO_frc_d10_01.nc)
(Tmin= 1.0000 Tmax= 30.0000)
(Min = -3.27691439E-04 Max = 2.14956658E-04)
GET_2DFLD - solar shortwave radiation flux, t = 1 00:00:00
(Rec=0000001, Index=2, File: IO_frc_d10_01.nc)
(Tmin= 1.0000 Tmax= 30.0000)
(Min = 2.25483297E-05 Max = 7.95320454E-05)
GET_2DFLD - surface net heat flux, t = 1 00:00:00
(Rec=0000001, Index=2, File: IO_frc_d10_01.nc)
(Tmin= 1.0000 Tmax= 30.0000)
(Min = -5.94075649E-05 Max = 4.05880867E-05)
GET_2DFLD - surface net freswater flux, (E-P), t = 1 00:00:00
(Rec=0000001, Index=2, File: IO_frc_d10_01.nc)
(Tmin= 1.0000 Tmax= 30.0000)
(Min = -1.70621406E-07 Max = 9.32263919E-08)
GET_NGFLD - free-surface eastern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 1.09218201E-01 Max = 8.56462430E-01)
GET_NGFLD - free-surface southern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 8.11142240E-02 Max = 7.40974584E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = -2.72014740E-01 Max = 8.39977933E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = -2.63939739E-01 Max = 5.91890203E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = -6.70640665E-01 Max = 1.22306183E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = -3.81899509E-01 Max = 8.58668860E-02)
GET_NGFLD - temperature eastern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 6.58755465E-01 Max = 2.94668588E+01)
GET_NGFLD - salinity eastern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 2.48780853E+01 Max = 3.57520208E+01)
GET_NGFLD - temperature southern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 4.08595854E-01 Max = 2.55525390E+01)
GET_NGFLD - salinity southern boundary condition, t = 1 00:00:00
(Rec=0000001, Index=1, File: IO_bry_d10.nc)
(Tmin= 1.0000 Tmax= 356.0000)
(Min = 3.43756636E+01 Max = 3.58334744E+01)
GET_3DFLD - potential temperature climatology, t = 1 00:00:00
(Rec=0000001, Index=2, File: IO_clm_d10_TSUVZ.nc)
(Tmin= 1.0000 Tmax= 271.0000)
(Min = 5.79100773E-02 Max = 2.97284260E+01)
GET_3DFLD - salinity climatology, t = 1 00:00:00
(Rec=0000001, Index=2, File: IO_clm_d10_TSUVZ.nc)
(Tmin= 1.0000 Tmax= 271.0000)
(Min = 2.12750411E+00 Max = 4.07142143E+01)
STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed
0 0 00:00:00 7.816705E-04 2.102642E+04 2.102642E+04 1.351798E+17
(114,286,36) 7.585973E-02 7.725372E-02 0.000000E+00 1.001434E+00
DEF_HIS - creating history file, Grid 01: /public/userspace/ljd/IO_d10/IO_his_d10_0001.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
1 0 00:20:00 7.853795E-04 2.102642E+04 2.102642E+04 1.351798E+17
(427,355,22) 9.046983E-03 1.832398E-02 3.615468E-01 1.001033E+00
WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001
2 0 00:40:00 7.884234E-04 2.102643E+04 2.102643E+04 1.351798E+17
(093,251,40) 5.191574E-03 3.447813E-03 1.963723E+00 1.002399E+00
WRT_RST - wrote re-start fields (Index=1,1) into time record = 0000002
3 0 01:00:00 7.909288E-04 2.102643E+04 2.102643E+04 1.351798E+17
(017,011,40) 2.232034E-02 7.210791E-03 1.492642E+00 1.002914E+00
WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001
4 0 01:20:00 7.933181E-04 2.102643E+04 2.102643E+04 1.351798E+17
(017,011,40) 1.075373E-02 1.031695E-03 1.974933E+00 1.005651E+00
WRT_RST - wrote re-start fields (Index=1,1) into time record = 0000002
5 0 01:40:00 7.955499E-04 2.102643E+04 2.102643E+04 1.351798E+17
(093,238,40) 4.780234E-03 2.868451E-04 1.005678E+00 1.005552E+00
WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001
6 0 02:00:00 7.977047E-04 2.102643E+04 2.102643E+04 1.351798E+17
(090,250,40) 1.649237E-02 1.200453E-02 1.148647E+00 1.003872E+00
WRT_RST - wrote re-start fields (Index=1,1) into time record = 0000002
7 0 02:20:00 7.999087E-04 2.102643E+04 2.102643E+04 1.351798E+17
(090,250,40) 4.511523E-03 1.395434E-02 1.543594E+00 1.001765E+00
WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001
8 0 02:40:00 8.023086E-04 2.102644E+04 2.102644E+04 1.351798E+17
(090,250,40) 1.885090E-02 9.810953E-04 1.732885E+00 1.006152E+00
WRT_RST - wrote re-start fields (Index=1,1) into time record = 0000002
9 0 03:00:00 8.047665E-04 2.102644E+04 2.102644E+04 1.351798E+17
(109,294,40) 1.139387E-02 2.705044E-03 9.080376E-01 1.003519E+00
WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001
10 0 03:20:00 8.074709E-04 2.102644E+04 2.102644E+04 1.351798E+17
(088,512,06) 6.383354E-02 5.553126E-02 9.530703E-01 1.004710E+00
WRT_RST - wrote re-start fields (Index=1,1) into time record = 0000002
11 0 03:40:00 8.110126E-04 2.102644E+04 2.102644E+04 1.351798E+17
(088,512,07) 6.624754E-02 5.991228E-02 1.224564E+00 1.019313E+00
WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001
12 0 04:00:00 8.154429E-04 2.102644E+04 2.102644E+04 1.351798E+17
(088,512,07) 7.194483E-02 5.467485E-02 1.655688E+00 1.168393E+00
WRT_RST - wrote re-start fields (Index=1,1) into time record = 0000002
13 0 04:20:00 8.207404E-04 2.102644E+04 2.102644E+04 1.351798E+17
(090,512,11) 5.065749E-03 1.447481E-02 2.383933E+00 1.304853E+00
WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001
14 0 04:40:00 8.288798E-04 2.102645E+04 2.102645E+04 1.351798E+17
(090,512,12) 1.696763E-02 4.843554E-02 3.511587E+00 1.659252E+00
Blowing-up: Saving latest model state into RESTART file
WRT_RST - wrote re-start fields (Index=1,1) into time record = 0000003
........
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 = 00000014
Analytical header files used:
ROMS/Functionals/ana_btflux.h
ROMS/TOMS: DONE... Thursday - April 3, 2014 - 12:23:49 PM
real 452.72
user 5805.97
sys 1272.95
+ exit 0
Blowing up at the beginning
Re: Blowing up at the beginning
Have you tried a smaller timestep? Does it run longer?
Have you looked at the restart file to see how and where it's blowing up?
Have you looked at the restart file to see how and where it's blowing up?