Index '7' of dimension 1 of array 'idsvar' above upper bound

Discussion of how to use ROMS on different regional and basin scale applications.

Moderators: arango, robertson

Post Reply
Message
Author
simaham
Posts: 8
Joined: Tue Oct 28, 2014 7:16 pm
Location: PTP

Index '7' of dimension 1 of array 'idsvar' above upper bound

#1 Unread post by simaham »

Dear ROMS users

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
and here is my ocean.in file

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

Post Reply