roms.in file output error...

Message

rakeshthp · #1 Unread post by **rakeshthp** » Mon Dec 31, 2007 4:04 am

Hello everybody..

After installing intel fortran compiler roms is working properly.. properly in the sense, the given example is only working.. To work with indian coast, i just changed the lon_min, lon_mat, lat_min and lat_max, and got the values of LLm and MMm values and entered them in param.h file.. and after running it i'm not getting appropriate output.. some error is coming.. what may be the problem..?? i'm getting following output after running
./roms roms.in

Southern Benguela
480 ntimes Total number of timesteps for 3D equations.
5400.00 dt Timestep [sec] for 3D equations
60 ndtfast Number of 2D timesteps within each 3D step.
1 ninfo Number of timesteps between runtime diagnostics.

6.000E+00 theta_s S-coordinate surface control parameter.
0.000E+00 theta_b S-coordinate bottom control parameter.
1.000E+01 Tcline S-coordinate surface/bottom layer width used in
vertical coordinate stretching, meters.
Grid File: ROMS_FILES/roms_grd.nc
Forcing Data File: ROMS_FILES/roms_frc.nc

WARNING: Unrecognized keyword: bulk_forcing --> DISREGARDED.

Climatology File: ROMS_FILES/roms_clm.nc

WARNING: Unrecognized keyword: boundary --> DISREGARDED.

Initial State File: ROMS_FILES/roms_ini.nc Record: 1
Restart File: ROMS_FILES/roms_rst.nc nrst = 480 rec/file: -1
History File: ROMS_FILES/roms_his.nc Create new: T nwrt = 480 rec/file = 0
1 ntsavg Starting timestep for the accumulation of output
time-averaged data.
48 navg Number of timesteps between writing of time-averaged
data into averages file.
Averages File: ROMS_FILES/roms_avg.nc rec/file = 0

Fields to be saved in history file: (T/F)
T write zeta free-surface.
F write UBAR 2D U-momentum component.
F write VBAR 2D V-momentum component.
F write U 3D U-momentum component.
F write V 3D V-momentum component.
F write T(1) Tracer of index 1.
F write T(2) Tracer of index 2.

F write RHO Density anomaly.
F write Omega Omega vertical velocity.
F write W True vertical velocity.
F write Akv Vertical viscosity.
F write Akt Vertical diffusivity for temperature.
F write Aks Vertical diffusivity for salinity.
F write Hbl Depth of KPP-model boundary layer.
F write Bostr Bottom Stress.

Fields to be saved in averages file: (T/F)
T write zeta free-surface.
T write UBAR 2D U-momentum component.
T write VBAR 2D V-momentum component.
T write U 3D U-momentum component.
T write V 3D V-momentum component.
T write T(1) Tracer of index 1.
T write T(2) Tracer of index 2.

F write RHO Density anomaly
T write Omega Omega vertical velocity.
T write W True vertical velocity.
F write Akv Vertical viscosity
T write Akt Vertical diffusivity for temperature.
F write Aks Vertical diffusivity for salinity.
T write Hbl Depth of KPP-model boundary layer
T write Bostr Bottom Stress.
1025.0000 rho0 Boussinesq approximation mean density, kg/m3.
0.000E+00 visc2 Horizontal Laplacian mixing coefficient [m2/s]
for momentum.
0.000E+00 tnu2(1) Horizontal Laplacian mixing coefficient (m2/s)
for tracer 1.
0.000E+00 tnu2(2) Horizontal Laplacian mixing coefficient (m2/s)
for tracer 2.
0.000E+00 rdrg Linear bottom drag coefficient (m/si).
0.000E+00 rdrg2 Quadratic bottom drag coefficient.
1.000E-02 Zob Bottom roughness for logarithmic law (m).
1.000E-04 Cdb_min Minimum bottom drag coefficient.
1.000E-01 Cdb_max Maximum bottom drag coefficient.

1.00 gamma2 Slipperiness parameter: free-slip +1, or no-slip -1.
1.00E+05 x_sponge Thickness of sponge and/or nudging layer (m)
800.00 v_sponge Viscosity in sponge layer (m2/s)
1.157E-05 tauT_in Nudging coefficients [sec^-1]
3.215E-08 tauT_out Nudging coefficients [sec^-1]
1.157E-06 tauM_in Nudging coefficients [sec^-1]
3.215E-08 tauM_out Nudging coefficients [sec^-1]

WARNING: Unrecognized keyword: diagnostics --> DISREGARDED.

WARNING: Unrecognized keyword: diag_avg --> DISREGARDED.

WARNING: Unrecognized keyword: diagnostics_bio --> DISREGARDED.

WARNING: Unrecognized keyword: diagbio_avg --> DISREGARDED.

WARNING: Unrecognized keyword: sediments --> DISREGARDED.

WARNING: Unrecognized keyword: sediment_history_fields --> DISREGARDED.

WARNING: Unrecognized keyword: bbl_history_fields --> DISREGARDED.

WARNING: Unrecognized keyword: floats --> DISREGARDED.

WARNING: Unrecognized keyword: float_fields --> DISREGARDED.

WARNING: Unrecognized keyword: stations --> DISREGARDED.

WARNING: Unrecognized keyword: station_fields --> DISREGARDED.

WARNING: Unrecognized keyword: psource --> DISREGARDED.

Activated C-preprocessing Options:

REGIONAL
BENGUELA
OBC_EAST
OBC_WEST
OBC_NORTH
OBC_SOUTH
SOLVE3D
UV_COR
UV_ADV
CURVGRID
SPHERICAL
MASKING
AVERAGES
AVERAGES_K
SALINITY
NONLIN_EOS
SPLIT_EOS
QCORRECTION
SFLX_CORR
DIURNAL_SRFLUX
SPONGE
CLIMATOLOGY
ZCLIMATOLOGY
M2CLIMATOLOGY
M3CLIMATOLOGY
TCLIMATOLOGY
ZNUDGING
M2NUDGING
M3NUDGING
TNUDGING
ANA_BSFLUX
ANA_BTFLUX
UV_VIS2
MIX_GP_UV
TS_DIF2
MIX_GP_TS
CLIMAT_TS_MIXH
LMD_MIXING
LMD_SKPP
LMD_BKPP
LMD_RIMIX
LMD_CONVEC
OBC_M2FLATHER
OBC_M3ORLANSKI
OBC_TORLANSKI
M2FILTER_COSINE

Linux 2.6.22.5-31-default i686
NUMBER OF THREADS: 1 BLOCKING: 1 x 1.
-11

Spherical grid detected.

-10

hmin hmax grdmin grdmax Cu_min Cu_max
75.000000 4521.358610 .104455839E+05 .109436150E+05 0.31910429 2.44980675
volume=2.937838787677007500000E+15 open_cross=9.547022086081026077271E+09

-9

Vertical S-coordinate System:

level S-coord Cs-curve at_hmin over_slope at_hmax

32 0.0000000 0.0000000 0.000 0.000 0.000
31 -0.0312500 -0.0009350 -0.373 -2.452 -4.531
30 -0.0625000 -0.0019030 -0.749 -4.979 -9.210
29 -0.0937500 -0.0029380 -1.128 -7.660 -14.192
28 -0.1250000 -0.0040767 -1.515 -10.578 -19.641
27 -0.1562500 -0.0053591 -1.911 -13.825 -25.739
26 -0.1875000 -0.0068304 -2.319 -17.504 -32.689
25 -0.2187500 -0.0085426 -2.743 -21.734 -40.726
24 -0.2500000 -0.0105560 -3.186 -26.654 -50.122
23 -0.2812500 -0.0129416 -3.654 -32.425 -61.197
22 -0.3125000 -0.0157835 -4.151 -39.240 -74.330
21 -0.3437500 -0.0191819 -4.684 -47.329 -89.974
20 -0.3750000 -0.0232566 -5.262 -56.965 -108.669
19 -0.4062500 -0.0281514 -5.892 -68.478 -131.064
18 -0.4375000 -0.0340388 -6.588 -82.262 -157.936
17 -0.4687500 -0.0411263 -7.361 -98.792 -190.223
16 -0.5000000 -0.0496640 -8.228 -118.640 -229.052
15 -0.5312500 -0.0599527 -9.209 -142.495 -275.781
14 -0.5625000 -0.0723554 -10.328 -171.187 -332.046
13 -0.5937500 -0.0873092 -11.613 -205.717 -399.821
12 -0.6250000 -0.1053416 -13.097 -247.290 -481.484
11 -0.6562500 -0.1270882 -14.823 -297.363 -579.903
10 -0.6875000 -0.1533158 -16.841 -357.689 -698.538
9 -0.7187500 -0.1849493 -19.209 -430.385 -841.560
8 -0.7500000 -0.2231040 -22.002 -518.002 -1014.002
7 -0.7812500 -0.2691252 -25.306 -623.619 -1221.933
6 -0.8125000 -0.3246355 -29.226 -750.949 -1472.672
5 -0.8437500 -0.3915923 -33.891 -904.471 -1775.051
4 -0.8750000 -0.4723564 -39.453 -1089.586 -2139.719
3 -0.9062500 -0.5697755 -46.098 -1312.811 -2579.524
2 -0.9375000 -0.6872846 -54.048 -1582.005 -3109.962
1 -0.9687500 -0.8290268 -63.574 -1906.649 -3749.725
0 -1.0000000 -1.0000000 -75.000 -2298.179 -4521.359

Time splitting: ndtfast = 60 nfast = 89

-8
-7

Maximum grid stiffness ratios: rx0 =0.2360003877 rx1 = 2.5917040544

GET_INITIAL -- Processing data for time = 0.000 record = 1
-6
-5
-4
-3
-2
GET_TCLIMA -- Read climatology of tracer 1 for time = 345.0
GET_TCLIMA -- Read climatology of tracer 1 for time = 15.00
GET_TCLIMA -- Read climatology of tracer 2 for time = 345.0
GET_TCLIMA -- Read climatology of tracer 2 for time = 15.00
GET_UCLIMA -- Read momentum climatology for time = 345.0
GET_UCLIMA -- Read momentum climatology for time = 15.00
GET_SSH - Read SSH climatology for time = 345.0
GET_SSH - Read SSH climatology for time = 15.00
GET_SMFLUX -- Read surface momentum stresses for time = 345.0
GET_SMFLUX -- Read surface momentum stresses for time = 15.00
GET_STFLUX -- Read surface flux of tracer 1 for time = 345.0
GET_STFLUX -- Read surface flux of tracer 1 for time = 15.00
GET_SST -- Read SST and dQdSST fields for time = 345.0
GET_SST -- Read SST and dQdSST fields for time = 15.00
GET_STFLUX -- Read surface flux of tracer 2 for time = 345.0
GET_STFLUX -- Read surface flux of tracer 2 for time = 15.00
GET_SSS -- Read SSS fields for time = 345.0
GET_SSS -- Read SSS fields for time = 15.00
GET_SRFLUX -- Read solar shortwave radiation for time = 345.0
GET_SRFLUX -- Read solar shortwave radiation for time = 15.00
-1
DEF_HIS/AVG - Created new netCDF file 'ROMS_FILES/roms_his.nc'.
WRT_GRID -- wrote grid data into file 'ROMS_FILES/roms_his.nc'.
WRT_HIS -- wrote history fields into time record = 1 / 1
0

MAIN: started time-steping.

STEP time[DAYS] KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME trd
0 0.00000 0.000000000E+00 2.0677142E+01 2.0677142E+01 2.9378388E+15 0

******************************************
ABNORMAL JOB END
USER TERMINATED
******************************************
VMAX = 0.295E+03 IMAX,JMAX = 52 13
IINT,IEXT = 1 8

MAIN - number of records written into history file(s): 1
number of records written into restart file(s): 0
number of records written into averages file(s): 0

MAIN: Abnormal termination: BLOWUP

So i wanted to know what all files need to be edited to get proper output..? Whether we need to edit roms.in file..?? if so then what part of that file must be edited..??

andres · #2 Unread post by **andres** » Wed Jan 02, 2008 12:52 pm

Hi,

A couple of questions:

Do you have a close boundary?
(Could you post a figure of your domain? )

What is the spatial resolution you are using?

The first thing I always try is to decrease the time step by half until it works. If it never works, then check other details. Please try this first.

Andres

ldebreu · #3 Unread post by **ldebreu** » Wed Jan 02, 2008 6:59 pm

Hello,

Your time step is too large.
As a rule of thumb, the value of Cu_max should be less than 0.8

> hmin hmax grdmin grdmax Cu_min Cu_max
> 75.000000 4521.358610 .104455839E+05 .109436150E+05 0.31910429 2.44980675

It is currently equal to 2.449...
So you have to divide the baroclinic time step (dt in roms.in) by at least
a factor of 3

Laurent

kate · #4 Unread post by **kate** » Wed Jan 02, 2008 8:04 pm

Also, it looks like you are attempting to use several cpp options available in the Rutgers ROMS which the Agrif ROMS is not recognizing. It is kind enough to let you know that it doesn't know about them, though: FLOATS, SEDIMENT, etc. Are any of these critical to your needs?

rakeshthp · #5 Unread post by **rakeshthp** » Thu Jan 03, 2008 4:26 am

hi..

thx for your reply.. but i really dont know about it.. i'm novice.. i dont have any knowledge of oceanography.. came in this field just three months back. am MCA student.. so i'm jjust trying to implement n run the model.. thats it.. what all you said about baroclonic, n floats n sediments all are going above my head... first of al i just need to run this model without any error.. next i want to understand what are all these fields present in that roms.in file..

anyways, i reduced the tiimestep n its working fine.. but i didnt understand the concept of Cu min which must be < 0.8 or something.

Thanks a lot..