blowing up: Diag Speed Trouble

[jitraporn13]

Dear All,

I'm trying to run ROMs by forcing with u,v,p components
It was blowing up since the beginning.
The error:

Diag Speed Trouble
Diag speed ijk

I don't know why the max speed is so big

Please help me to figure out the problem.

Thank you so much

Here is part of my log file
____________________________________________

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

3600 ntimes Number of timesteps for 3-D equations.
150.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.
T LcycleRST Switch to recycle time-records in restart file.
30 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).
30 nHIS Number of timesteps between the writing fields
into history file.
0 nQCK Number of timesteps between the writing fields
into quicksave file.
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.
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.0000E+01 blk_ZQ Height (m) of surface air humidity measurement.
1.0000E+01 blk_ZT Height (m) of surface air temperature measurement.
1.0000E+01 blk_ZW Height (m) of surface winds measurement.
1.0000E-01 Dcrit Minimum depth for wetting and drying (m).
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
0.0000E+00 theta_s S-coordinate surface control parameter.
1.0000E-01 theta_b S-coordinate bottom control parameter.
100.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
5844.000 dstart Time-stamp assigned to model initialization (days).
20000101.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
9.0000E+01 Tnudg_SSS Nudging/relaxation time scale (days)
for sea surface salinity.
0.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
0.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
0.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
0.0000E+00 obcfac Factor between passive and active
open boundary conditions.
F VolCons(1) NLM western edge boundary volume conservation.
F VolCons(2) NLM southern edge boundary volume conservation.
F VolCons(3) NLM eastern edge boundary volume conservation.
F VolCons(4) NLM northern edge boundary volume conservation.
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.
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

Output/Input Files:

Output Restart File: SE_rst.nc
Output History File: SE_his.nc
Input Grid File: ../../../easygrid_v2/SE_Asia_grd.nc
Input Nonlinear Initial File: ../../../easygrid_v2/SE_ASIA2_ini.nc
Input Forcing File 01: ../../../Documents/nesting/matlab/forcing/utest.nc
Input Forcing File 02: ../../../Documents/nesting/matlab/forcing/vtest.nc
Input Forcing File 03: ../../../Documents/nesting/matlab/forcing/ptest.nc
Input Forcing File 04: ../../../Documents/nesting/matlab/forcing/SE_Tair_era.nc
Input Forcing File 05: ../../../Documents/nesting/matlab/forcing/SE_Qair_era.nc
Input Forcing File 06: ../../../Documents/nesting/matlab/forcing/SE_lwrad_era.nc
Input Forcing File 07: ../../../Documents/nesting/matlab/forcing/SE_rain_era.nc

Tile partition information for Grid 01: 0258x0221x0030 tiling: 001x001

tile Istr Iend Jstr Jend Npts

Number of tracers: 2
0 1 258 1 221 1710540

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 259.50 0.50 221.50 RHO-points

0 0.00 259.00 0.50 221.50 U-points

0 0.50 259.50 0.00 221.00 V-points

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

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

zeta 1 Closed Closed Closed Closed

ubar 1 Closed Closed Closed Closed

vbar 1 Closed Closed Closed Closed

u 1 Closed Closed Closed Closed

v 1 Closed Closed Closed Closed

temp 1 Closed Closed Closed Closed

salt 1 Closed Closed Closed Closed

tke 1 Closed Closed Closed Closed

Activated C-preprocessing Options:

SE2 SE Asia
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_SRFLUX Analytical kinematic shortwave radiation flux.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ASSUMED_SHAPE Using assumed-shape arrays.
BULK_FLUXES Surface bulk fluxes parameterization.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
GLS_MIXING Generic Length-Scale turbulence closure.
LIMIT_BSTRESS Limit bottom stress to maintain bottom velocity direction.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
LIMIT_BSTRESS Limit bottom stress to maintain bottom velocity direction.
MASKING Land/Sea masking.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
OUT_DOUBLE Double precision output fields in NetCDF files.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
RI_SPLINES Parabolic Spline Reconstruction for Richardson Number.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
SPLINES_VDIFF Parabolic Spline Reconstruction for Vertical Diffusion.
SPLINES_VVISC Parabolic Spline Reconstruction for Vertical Viscosity.
TS_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection 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.
VAR_RHO_2D Variable density barotropic mode.
WET_DRY Wetting and drying activated.

Process Information:

Thread # 0 (pid= 61980) 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

30 0.0000000 0.0000000 0.000 0.000 0.000 0.000
29 -0.0333333 -0.0011675 -0.128 -1.725 -6.319 -9.585
28 -0.0666667 -0.0046693 -0.257 -3.567 -18.842 -31.790
27 -0.1000000 -0.0105031 -0.386 -5.525 -37.563 -66.601
26 -0.1333333 -0.0186649 -0.516 -7.600 -62.470 -113.997
25 -0.1666667 -0.0291493 -0.646 -9.791 -93.550 -173.950
24 -0.2000000 -0.0419494 -0.776 -12.097 -130.784 -246.421
23 -0.2333333 -0.0570566 -0.906 -14.519 -174.150 -331.364
22 -0.2666667 -0.0744609 -1.037 -17.056 -223.620 -428.725
21 -0.3000000 -0.0941507 -1.168 -19.708 -279.164 -538.442
20 -0.3333333 -0.1161130 -1.300 -22.472 -340.748 -660.443
19 -0.3666667 -0.1403332 -1.432 -25.350 -408.332 -794.651
18 -0.4000000 -0.1667954 -1.564 -28.340 -481.874 -940.979
17 -0.4333333 -0.1954820 -1.697 -31.441 -561.329 -1099.333
16 -0.4666667 -0.2263742 -1.830 -34.652 -646.644 -1269.609
15 -0.5000000 -0.2594516 -1.963 -37.973 -737.768 -1451.699
14 -0.5333333 -0.2946926 -2.097 -41.401 -834.641 -1645.485
13 -0.5666667 -0.3320740 -2.231 -44.937 -937.202 -1850.842
12 -0.6000000 -0.3715715 -2.365 -48.579 -1045.387 -2067.639
11 -0.6333333 -0.4131593 -2.499 -52.325 -1159.128 -2295.736
10 -0.6666667 -0.4568104 -2.634 -56.174 -1278.352 -2534.987
9 -0.7000000 -0.5024966 -2.770 -60.125 -1402.984 -2785.239
8 -0.7333333 -0.5501882 -2.905 -64.176 -1532.946 -3046.333
7 -0.7666667 -0.5998546 -3.041 -68.326 -1668.156 -3318.103
6 -0.8000000 -0.6514640 -3.177 -72.573 -1808.530 -3600.376
5 -0.8333333 -0.7049834 -3.314 -76.916 -1953.980 -3892.975
4 -0.8666667 -0.7603786 -3.450 -81.352 -2104.415 -4195.714
3 -0.9000000 -0.8176145 -3.587 -85.881 -2259.741 -4508.404
2 -0.9333333 -0.8766550 -3.725 -90.499 -2419.864 -4830.850
1 -0.9666667 -0.9374628 -3.862 -95.206 -2584.684 -5162.850
0 -1.0000000 -1.0000000 -4.000 -100.000 -2754.099 -5504.199

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

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

Minimum X-grid spacing, DXmin = 1.50000000E+01 km
Maximum X-grid spacing, DXmax = 1.50000333E+01 km
Minimum Y-grid spacing, DYmin = 1.50000000E+01 km
Maximum Y-grid spacing, DYmax = 1.55908600E+01 km
Minimum Z-grid spacing, DZmin = 1.28384748E-01 m
Maximum Z-grid spacing, DZmax = 3.41348659E+02 m

Minimum barotropic Courant Number = 1.94485509E-02
Maximum barotropic Courant Number = 6.57079840E-01
Maximum Coriolis Courant Number = 9.20089964E-03


NLM: GET_STATE - Read state initial conditions, t = 5844 00:00:00
(Grid 01, File: SE_ASIA2_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 = 9.00000000E+00 Max = 9.00000000E+00)
- salinity
(Min = 3.00000000E+01 Max = 3.00000000E+01)
GET_2DFLD - surface u-wind component, t = 5844 00:00:00
(Rec=0000001, Index=1, File: utest.nc)
(Tmin= 5844.0000 Tmax= 5875.0000)
(Min = -1.22636754E+01 Max = 4.81609346E+00)
GET_2DFLD - surface v-wind component, t = 5844 00:00:00
(Rec=0000001, Index=1, File: vtest.nc)
(Tmin= 5844.0000 Tmax= 5875.0000)
(Min = -1.21289526E+01 Max = 3.68782660E+00)
GET_2DFLD - surface air pressure, t = 5844 00:00:00
(Rec=0000001, Index=1, File: ptest.nc)
(Tmin= 5844.0000 Tmax= 5875.0000)
(Min = 1.01185653E+01 Max = 1.03030460E+01)
GET_2DFLD - surface air temperature, t = 5844 00:00:00
(Rec=0000001, Index=1, File: SE_Tair_era.nc)
(Tmin= 5844.0000 Tmax= 5875.0000)
(Min = 2.56414107E+03 Max = 2.08692103E+04)
GET_2DFLD - surface air relative humidity, t = 5844 00:00:00
(Rec=0000001, Index=1, File: SE_Qair_era.nc)
(Tmin= 5844.0000 Tmax= 5875.0000)
(Min = 5.00000000E+01 Max = 5.00000000E+01)
GET_2DFLD - net longwave radiation flux, t = 5843 22:30:00
(Rec=0000001, Index=1, File: SE_lwrad_era.nc)
(Tmin= 5843.9375 Tmax= 5874.9375)
(Min = 6.78818394E-07 Max = 7.42758226E-07)
GET_2DFLD - rain fall rate, t = 5843 22:30:00
(Rec=0000001, Index=1, File: SE_rain_era.nc)
(Tmin= 5843.9375 Tmax= 5874.9375)
(Min = -2.73776874E+00 Max = -1.20218290E+00)

Basin information for Grid 01:

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

Initial basin volumes: TotVolume = 1.3989106538E+16 m3
MinVolume = 3.4621784284E+07 m3
MaxVolume = 7.6841808058E+10 m3
Max/Min = 2.2194641220E+03

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

GET_2DFLD - surface u-wind component, t = 5844 12:00:00
(Rec=0000002, Index=2, File: utest.nc)
(Tmin= 5844.0000 Tmax= 5875.0000)
(Min = -1.19747591E+01 Max = 3.37829314E+00)
GET_2DFLD - surface v-wind component, t = 5844 12:00:00
(Rec=0000002, Index=2, File: vtest.nc)
(Tmin= 5844.0000 Tmax= 5875.0000)
(Min = -1.27883484E+01 Max = 2.67607881E+00)
GET_2DFLD - surface air pressure, t = 5844 12:00:00
(Rec=0000002, Index=2, File: ptest.nc)
(Tmin= 5844.0000 Tmax= 5875.0000)
(Min = 1.01103045E+01 Max = 1.02726652E+01)
GET_2DFLD - surface air temperature, t = 5844 12:00:00
(Rec=0000002, Index=2, File: SE_Tair_era.nc)
(Tmin= 5844.0000 Tmax= 5875.0000)
(Min = 3.66556213E+03 Max = 2.04005658E+04)
GET_2DFLD - surface air relative humidity, t = 5844 12:00:00
(Rec=0000002, Index=2, File: SE_Qair_era.nc)
(Tmin= 5844.0000 Tmax= 5875.0000)
(Min = 5.00000000E+01 Max = 5.00000000E+01)
GET_2DFLD - net longwave radiation flux, t = 5844 10:30:00
(Rec=0000002, Index=2, File: SE_lwrad_era.nc)
(Tmin= 5843.9375 Tmax= 5874.9375)
(Min = -3.83092890E-08 Max = 6.42174955E-07)
GET_2DFLD - rain fall rate, t = 5844 10:30:00
(Rec=0000002, Index=2, File: SE_rain_era.nc)
(Tmin= 5843.9375 Tmax= 5874.9375)
(Min = -3.03388892E+03 Max = -2.52202799E+03)

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

0 5844 00:00:00 0.000000E+00 1.508919E+04 1.508919E+04 1.455726E+16
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: SE_his.nc
WRT_HIS - wrote history fields (Index=1,1) in record = 0000001
1 5844 00:02:30 2.761005E+01 1.572569E+04 1.575330E+04 1.398958E+16
(104,113,30) 1.410010E+01 1.268761E+01 0.000000E+00 1.898574E+03
DIAG speed trouble 1898.5744590090233
DIAG speed ijk 1509.5531256078064 103 113 30 1

Blowing-up: Saving latest model state into RESTART file

WRT_RST - wrote re-start fields (Index=2,2) in record = 0000001

Elapsed CPU time (seconds):

Thread # 0 CPU: 121.579
Total: 121.579

Nonlinear ocean model elapsed time profile, Grid: 01

Allocation and array initialization .............. 13.727 (11.2906 %)
Ocean state initialization ....................... 2.178 ( 1.7912 %)
Reading of input data ............................ 8.801 ( 7.2387 %)
Processing of input data ......................... 0.102 ( 0.0839 %)
Computation of vertical boundary conditions ...... 0.075 ( 0.0617 %)
Computation of global information integrals ...... 4.167 ( 3.4277 %)
Writing of output data ........................... 3.613 ( 2.9721 %)
Model 2D kernel .................................. 5.005 ( 4.1169 %)
2D/3D coupling, vertical metrics ................. 3.077 ( 2.5305 %)
Omega vertical velocity .......................... 1.427 ( 1.1735 %)
Equation of state for seawater ................... 7.637 ( 6.2814 %)
Atmosphere-Ocean bulk flux parameterization ...... 0.697 ( 0.5732 %)
GLS vertical mixing parameterization ............. 26.115 (21.4800 %)
3D equations right-side terms .................... 5.254 ( 4.3217 %)
3D equations predictor step ...................... 9.365 ( 7.7025 %)
Pressure gradient ................................ 5.271 ( 4.3356 %)
Corrector time-step for 3D momentum .............. 7.315 ( 6.0166 %)
Corrector time-step for tracers .................. 8.237 ( 6.7748 %)
Total: 112.062 92.1725

All percentages are with respect to total time = 121.579


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/Functionals/ana_srflux.h
ROMS/Functionals/ana_stflux.h

[kate]

DIAG speed ijk 1509.5531256078064 103 113 30 1

It is at least telling you to look at point (i,j,k) = (103,113,30)

GET_2DFLD - surface air pressure, t = 5844 00:00:00
(Rec=0000001, Index=1, File: ptest.nc)
(Tmin= 5844.0000 Tmax= 5875.0000)
(Min = 1.01185653E+01 Max = 1.03030460E+01)

What units is your pressure in? ROMS wants something around 1000. Maybe you're using my varinfo.dat for converting from Pascals to millibars? Better check your Qair and Tair too.