when I run Roms I got memory error as following"
Program received signal SIGSEGV: Segmentation fault - invalid memory reference.
Backtrace for this error:
#0 0x7F3599E25E08
#1 0x7F3599E24F90
#2 0x7F359976C4AF
#3 0x4ECA9F in __analytical_mod_MOD_ana_btflux
#4 0x472D03 in set_data_tile_
#5 0x475447 in set_data_
#6 0x438690 in main3d_
Segmentation fault (core dumped)
So I enabled use_debug option then after compiling and running it I got the following error
]At line 1260 of file mod_ncparam.f90
Fortran runtime error: Index '7' of dimension 1 of array 'idsvar' above upper bound of 6
any idea what is the "idsvar"?
could you please help me to solvbe it
here is my include file
Code: Select all
** Application flag: CURRENT_WIND
** Input script: ocean_current_wind.in
*/
#define UV_ADV /* use to turn ON or OFF advection terms */
#define UV_COR /* use to turn ON or OFF Coriolis term */
#define UV_QDRAG /* use to turn ON or OFF quadratic bottom friction */
#define UV_VIS4 /* use to turn ON or OFF harmonic horizontal mixing */
#define MIX_S_UV /* momentum mixing on s-surfaces */
#define TS_U3HADVECTION /* use if 3rd-order upstream horiz. advection */
#define TS_C4VADVECTION /* use if 4th-order centered vertical advection */
#define SOLVE3D /* use if solving 3D primitive equations */
#define ANA_SMFLUX /* use if analytical surface momentum stress */
#define ANA_STFLUX /* use if analytical surface temperature flux */
#define ANA_SSFLUX /* use if analytical surface salinity flux */
#define ANA_BSFLUX /* use if analytical bottom salinity flux */
#define ANA_BTFLUX /* use if analytical bottom temperature flux */
#define MASKING /* use if analytical masking is enabled */
#define SOUTH_FSCHAPMAN /* use if free-surface Chapman condition*/
#define SOUTH_M2FLATHER /* use if 2D momentum Flather condition*/
#define SOUTH_TRADIATION /* use if tracers radiation condition*/
#define WEST_FSCHAPMAN /* use if free-surface Chapman condition*/
#define WEST_M2FLATHER /* use if 2D momentum Flather condition*/
#define WEST_TRADIATION
#define ANA_FSOBC /* use if analytical free-surface boundary conditions*/
#define ANA_M2OBC /* use if analytical 2D momentum boundary conditions*/
#undef BULK_FLUXES
Code: Select all
! Application title.
!
TITLE = Current Wind
!
! C-preprocessing Flag.
!
MyAppCPP = CURRENT_WIND
!
! Input variable information file name. This file needs to be processed
! first so all information arrays can be initialized properly.
!
VARNAME = /home/user1/SIMA/Ocean/ROMS/External/varinfo.dat
!
! Grid dimension parameters. See notes below in the Glossary for how to set
! these parameters correctly.
!
Ngrids == 1
Lm == 197 ! Number of I-direction INTERIOR RHO-points
Mm == 224 ! Number of J-direction INTERIOR RHO-points
N == 10 ! Number of vertical levels
! Domain decomposition parameters for serial, distributed-memory or
! shared-memory configurations used to determine tile horizontal range
! indices (Istr,Iend) and (Jstr,Jend), [1:Ngrids].
!
NtileI == 4 ! I-direction partition
NtileJ == 1 ! J-direction partition
! Time-Stepping parameters.
!
NTIMES == 864864 ! 10 days
DT == 100
NDTFAST == 20
! Number of eigenvalues (NEV) and eigenvectors (NCV) to compute for the
! Lanczos/Arnoldi problem in the Generalized Stability Theory (GST)
! analysis. NCV must be greater than NEV (see documentation below).
!
NEV = 2 ! Number of eigenvalues
NCV = 10 ! Number of eigenvectors
!
! Input/Output parameters.
!
NRREC == 0
LcycleRST == T
NRST == 360 ! Every 1 hour
NSTA == 360 ! Every 1 hour
NFLT == 360 ! Every 1 hour
NINFO == 360 ! Every 1 hour
!
! Output history, average, diagnostic files parameters.
!
LDEFOUT == T
NHIS == 60 ! Every 10 minutes
NDEFHIS == 0
NTSAVG == 1
NAVG == 60 ! Every 10 minutes
NDEFAVG == 0
NTSDIA == 1
NDIA == 60 ! Every 10 minutes
NDEFDIA == 0
!
! Harmonic/biharmonic horizontal diffusion of tracer: [1:NAT+NPT,Ngrids].
!
TNU2 == 20.0d0 20.0d0 ! m2/s
TNU4 == 2*0.0d0 ! m4/s
!
! Harmononic/biharmonic, horizontal viscosity coefficient: [Ngrids].
!
VISC2 == 100.0d0 ! m2/s
VISC4 == 0.0d0 ! m4/s
!
! Vertical mixing coefficients for active tracers: [1:NAT+NPT,Ngrids]
!
AKT_BAK == 1.0d-6 1.0d-6 ! m2/s
! Set vertical, terrain-following coordinates transformation equation and
! stretching function (see below for details), [1:Ngrids].
Vtransform == 2 ! transformation equation
Vstretching == 1 ! stretching function
! Vertical S-coordinates parameters, [1:Ngrids].
!
THETA_S == 7.0d0 ! 0 < THETA_S < 20
THETA_B == 0.1d0 ! 0 < THETA_B < 1
TCLINE == 5.0d0 ! m
! Time-stamp assigned for model initialization, reference time
! origin for tidal forcing, and model reference time for output
! NetCDF units attribute.
!
DSTART = 40639.0d0 ! days
TIDE_START = 0.0d0 ! days
TIME_REF = 19000101.00 ! yyyymmdd.dd
! Logical switches (TRUE/FALSE) to specify the state surface forcing
! variable whose stochastic optimals is required.
!
SOstate(isUstr) == F ! surface u-stress
SOstate(isVstr) == F ! surface v-stress
Hout(inert) == F ! inert passive tracers
Hout(idBott) == F F F F F F F F F F F F F F F F
! Input NetCDF file names, [1:Ngrids].
!
GRDNAME == /home/user1/SIMA/Ocean/ROMS/ROMS_Grid2.nc
ININAME == /home/user1/SIMA/Ocean/ROMS/ROMS_initial.nc
! ITLNAME == ocean_itl.nc
! IRPNAME == ocean_irp.nc
! IADNAME == ocean_iad.nc
! CLMNAME == ocean_clm.nc
! BRYNAME == ocean_bry.nc
! FWDNAME == ocean_fwd.nc
! ADSNAME == ocean_ads.nc
NFFILES == 1 ! number of forcing files
!
FRCNAME == ECMWF_atm_2014.nc ! forcing file 1, grid 1
! Output NetCDF file names, [1:Ngrids].
!
GSTNAME == ocean_gst.nc
RSTNAME == ocean_rst.nc
HISNAME == ocean_his.nc
TLMNAME == ocean_tlm.nc
TLFNAME == ocean_tlf.nc
ADJNAME == ocean_adj.nc
AVGNAME == ocean_avg.nc
DIANAME == ocean_dia.nc
STANAME == ocean_sta.nc
FLTNAME == ocean_flt.nc