Interannual simulation

[anu]

hello everyone,

I am getting error while interannual simulation as shown below. kindly help me out if anyone know how to solve it.

0 2003-01-01 00:00:00.00 5.757110E-04 1.520111E+04 1.520111E+04 7.950646E+15
(045,138,32) 1.858946E-02 1.080659E-02 0.000000E+00 6.403090E-01
DEF_AVG - creating average file, Grid 01: roms_avg.nc
GET_2DFLD - surface u-momentum stress, 2003-01-01 12:00:00.00
(Grid=01, Rec=0009498, Index=1, File: sustr_era_1990-2016.nc)
(Tmin= 0.0000 Tmax= 9860.5000) t = 4748.5000
(Min = -1.20152608E-04 Max = 1.49003681E-04) regrid = T
GET_2DFLD - surface v-momentum stress, 2003-01-01 12:00:00.00
(Grid=01, Rec=0009498, Index=1, File: svstr_era_1990-2016.nc)
(Tmin= 0.0000 Tmax= 9860.5000) t = 4748.5000
(Min = -6.81670388E-05 Max = 1.40685493E-04) regrid = T
GET_2DFLD - surface net heat flux, 2003-01-01 12:00:00.00
(Grid=01, Rec=0009498, Index=1, File: shflux_era_1990-2016.nc)
(Tmin= 0.0000 Tmax= 9860.5000) t = 4748.5000
(Min = -3.85925288E-05 Max = 3.96117364E-07) regrid = T
GET_2DFLD - surface net freswater flux, (E-P), 2003-01-01 12:00:00.00
(Grid=01, Rec=0009498, Index=1, File: swflux_era_1990-2016.nc)
(Tmin= 0.0000 Tmax= 9860.5000) t = 4748.5000
(Min = -1.45238460E-08 Max = 1.18113369E-09) regrid = T
1 2003-01-01 00:05:00.00 NaN NaN NaN NaN
(164,232,18) 2.716270E-03 1.262424E-03 6.567407E-02 4.635102E-01
Found Error: 01 Line: 298 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 298 Source: ROMS/Drivers/nl_ocean.h

Blowing-up: Saving latest model state into RESTART file

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

Elapsed CPU time (seconds):

Node # 3 CPU: 3.358
Node # 1 CPU: 3.454
Node # 0 CPU: 3.358
Node # 2 CPU: 3.408

[kate]

It blew up in the first step - this could be any number of things. The first thing I would try is a shorter barotropic timestep - i.e., a larger NDTFAST. Your surface inputs appear to be well behaved. What other inputs do you have?

[anu]

Hello Kate,
Thank you for your suggestion, we tried to reduce the NDTFAST and the NaN values were also reduced but not completely removed.

My Cppdef is as follows:
# define NO_LBC_ATT
#define UV_ADV
#define UV_COR
#define UV_QDRAG
#define DJ_GRADPS
#define SPLINES_VDIFF
#define SPLINES_VVISC
#define TS_U3HADVECTION
#define TS_C4VADVECTION
#define NONLIN_EOS
#define SALINITY
#define SOLVE3D
#define MASKING
#define CURVGRID
#define AVERAGES
#define MIX_S_UV
#define MIX_GEO_TS
#define ANA_BSFLUX
#define ANA_BTFLUX
#define TCLM_NUDGING
#define TCLIMATOLOGY
#undef TIDES
#undef TS_FIXED
#undef BODYFORCE /*use if applying stresses as a body force */
#define SOLAR_SOURCE
#define SHORTWAVE
#define BULF_FLUXES /*bulk formulation for surface heat fluxes*/
#undef QCORRECTION/*Activate linearized bulk formulation providing heat flux correction dQdSST for nudging towards model SST*/
#define DIURNAL_SRFLUX/*impose shortwave radiation local diurnal cycle*/
#undef SRELAXATION /*salinity relaxation as a freshwater flux */
#define EMINUSP

/*# define RADIATION_2D use if tangential phase speed in radiation conditions*/

# define UV_VIS2/* harmonic horizontal mixing */

# define TS_DIF2


# define GLS_MIXING
# if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# define RI_SPLINES
# endif

# undef DIAGNOSTICS_TS
# undef DIAGNOSTICS_UV
# undef DIAGNOSTICS_BIO
-----------------------------------------------------------------
My input files are:
1) surface forcing u-stress and v-stress, Net surface heat flux derived from ECMWF
2)Ocean boundary and climate files are derived from SODA

[kate]

Code: Select all

#define TCLM_NUDGING
#define TCLIMATOLOGY

This is from an older ROMS version.

Code: Select all

#define BULF_FLUXES /*bulk formulation for surface heat fluxes*/ 

You have misspelled BULK_FLUXES. If you spell it correctly, it will ask for winds instead of surface stress. ROMS will ignore random strings that it doesn't recognize so you have to check the output cppdefs to see if they match what you think you are doing.

But really, you need to look at the ROMS output 3-D fields to see what is going wrong. You can also try to reduce DT, but if that doesn't help, there's something else going on.

[anu]

Thank you Kate! The problem was fixed and the simmulation has begin. But again it is showing an error as shown below. Any help would be appreciated.

1373328 2003-01-21 12:00:00.00 1.912849E-03 1.520177E+04 1.520177E+04 7.950547E+15
(036,069,21) 8.275453E-03 1.443629E-02 2.116743E-01 1.571334E+00

Elapsed CPU time (seconds):

Node # 3 CPU: 13854.680
Node # 1 CPU: 13843.423
Node # 2 CPU: 13832.210
Node # 0 CPU: 13770.997
Total: 55301.311

Nonlinear model elapsed CPU time profile, Grid: 01

Allocation and array initialization .............. 0.441 ( 0.0008 %)
Ocean state initialization ....................... 0.152 ( 0.0003 %)
Reading of input data ............................ 5.892 ( 0.0107 %)
Processing of input data ......................... 125.662 ( 0.2272 %)
Processing of output time averaged data .......... 785.375 ( 1.4202 %)
Computation of vertical boundary conditions ...... 20.059 ( 0.0363 %)
Computation of global information integrals ...... 423.818 ( 0.7664 %)
Writing of output data ........................... 6.217 ( 0.0112 %)
Model 2D kernel .................................. 19957.787 (36.0892 %)
2D/3D coupling, vertical metrics ................. 786.063 ( 1.4214 %)
Omega vertical velocity .......................... 537.169 ( 0.9714 %)
Equation of state for seawater ................... 830.333 ( 1.5015 %)
GLS vertical mixing parameterization ............. 14723.375 (26.6239 %)
3D equations right-side terms .................... 2591.168 ( 4.6855 %)
3D equations predictor step ...................... 3981.610 ( 7.1998 %)
Pressure gradient ................................ 1150.535 ( 2.0805 %)
Harmonic mixing of tracers, geopotentials ........ 3971.411 ( 7.1814 %)
Harmonic stress tensor, S-surfaces ............... 1445.786 ( 2.6144 %)
Corrector time-step for 3D momentum .............. 1431.058 ( 2.5877 %)
Corrector time-step for tracers .................. 1938.871 ( 3.5060 %)
Reading model state vector ....................... 4.125 ( 0.0075 %)
Total: 54716.910 98.9432

Nonlinear model message Passage profile, Grid: 01

Message Passage: 2D halo exchanges ............... 2780.045 ( 5.0271 %)
Message Passage: 3D halo exchanges ............... 1925.272 ( 3.4814 %)
Message Passage: 4D halo exchanges ............... 583.991 ( 1.0560 %)
Message Passage: data broadcast .................. 8.847 ( 0.0160 %)
Message Passage: data reduction .................. 51.361 ( 0.0929 %)
Message Passage: data gathering .................. 1.008 ( 0.0018 %)
Message Passage: data scattering.................. 0.823 ( 0.0015 %)
Message Passage: synchronization barrier ......... 0.024 ( 0.0000 %)
Total: 5351.369 9.6767

Unique code regions profiled ..................... 54716.910 98.9432 %
Residual, non-profiled code ...................... 584.401 1.0568 %


All percentages are with respect to total time = 55301.311


>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

Dynamic and Automatic memory (MB) usage for Grid 01: 287x261x32 tiling: 2x2

tile Dynamic Automatic USAGE MPI-Buffers

0 247.52 40.13 287.65 21.40
1 245.84 40.13 285.98 21.40
2 245.68 40.13 285.81 21.40
3 244.02 40.13 284.15 21.40

TOTAL 983.06 160.53 1143.59 85.60

<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

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

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_srflux.h

Thanks
Anu

[kate]

I don't see an error in the output you show. That's just the normal ROMS report at the end of a run. How many timesteps did you ask it to do?

[anu]

Thanks for your reply kate.

If we are using retart file for the interannual simulation, can we need to change the NTIMES also.As we are running the model for one year using retart file we had given time steps for one year.

Thank you in advance
Anu

[kate]

It depends on your ROMS version. Older versions required you to change NTIMES to get a second year while newer ones don't.

[anu]

Thanks a lot for the reply kate, it is working.