数组越界，跪求大大帮忙， — 编程爱好者社区

主题：数组越界，跪求大大帮忙，

dingsheng1206 [专家分：0] 发布于 2010-09-12 20:28:00

废话不多说，上传程序吧，
    program input

c----------------------------------------------------------------------
c     Subroutine for loading or generating initial particle information

c     x-- coordinates of particles                                 [out]
c     vx-- velocities of particles                                 [out]
c     mass-- mass of particles                                     [out]
c     rho-- dnesities of particles                                 [out]
c     p-- pressure  of particles                                   [out]
c     u-- internal energy of particles                             [out]
c     itype-- types of particles                                   [out]
c     hsml-- smoothing lengths of particles                        [out]
c     ntotal-- total particle number                               [out]

      implicit none
      include 'param.inc'

      integer itype(maxn), ntotal
      double precision x(dim, maxn), vx(dim, maxn), mass(maxn),
     &                 p(maxn), u(maxn), hsml(maxn), rho(maxn)
      integer i, d,im

c     load initial particle information from external disk file

      call shear_cavity(x, vx, mass, rho, p, u,
     &                        itype, hsml, ntotal)


      end



      subroutine shear_cavity(x, vx, mass, rho, p, u,
     &                        itype, hsml, ntotal)

c----------------------------------------------------------------------
c     This subroutine is used to generate initial data for the
c     2 d shear driven cavity probem with Re = 1
c     x-- coordinates of particles                                 [out]
c     vx-- velocities of particles                                 [out]
c     mass-- mass of particles                                     [out]
c     rho-- dnesities of particles                                 [out]
c     p-- pressure  of particles                                   [out]
c     u-- internal energy of particles                             [out]
c     itype-- types of particles                                   [out]
c          =2   water
c     h-- smoothing lengths of particles                           [out]
c     ntotal-- total particle number                               [out]

      implicit none
      include 'param.inc'

      integer itype(maxn), ntotal
      double precision x(dim, maxn), vx(dim, maxn), mass(maxn),
     &     rho(maxn), p(maxn), u(maxn), hsml(maxn)
      integer i, j, d, m, n, mp, np, k
      double precision xl, yl, dx, dy

c     Giving mass and smoothing length as well as other data.

      m = 41
      n = 41
      mp = m-1
      np = n-1
      ntotal = mp * np
      xl = 1.e-3
      yl = 1.e-3
      dx = xl/mp
      dy = yl/np

      do i = 1, mp
      do j = 1, np
      k = j + (i-1)*np
      x(1, k) = (i-1)*dx + dx/2.
      x(2, k) = (j-1)*dy + dy/2.
        enddo
      enddo

      do i = 1, mp*np
      vx(1, i) = 0.
      vx(2, i) = 0.
        rho (i) = 1000.
        mass(i) = dx*dy*rho(i)
        p(i)= 0.
        u(i)=357.1
        itype(i) = 2
        hsml(i) = dx
      enddo

      end

居然报错  x(2, k) = (j-1)*dy + dy/2.   这行是数组越界了
我看了没问题啊，这里不是只要行循环时变量i和j在限定范围就可以了吗，谢谢帮忙，

本帖地址： http://bbs.pfan.cn/post/349248.html

回复列表（共6个回复）

沙发

jstzhurj [专家分：4680] 发布于 2010-09-12 21:58:00

include 'param.inc'
param.inc 没上传啊[em15]

板凳

dingsheng1206 [专家分：0] 发布于 2010-09-12 22:28:00

[quote]
include 'param.inc'
param.inc 没上传啊[em15][/quote]
它的我看了啊，谢谢你哦
c---------------------------------------------------------
c     Including file for parameters and constants used
c     in the entire SPH software packages.
c---------------------------------------------------------

c     dim : Dimension of the problem (1, 2 or 3)
      integer dim
      parameter ( dim = 2)

c     maxn    : Maximum number of particles
c     max_interation : Maximum number of interaction pairs
      integer maxn,max_interaction
      parameter ( maxn    = 12000    ,
     &            max_interaction = 100 * maxn )

c     Parameters for the computational geometry,
c     x_maxgeom : Upper limit of allowed x-regime
c     x_mingeom : Lower limit of allowed x-regime
c     y_maxgeom : Upper limit of allowed y-regime
c     y_mingeom : Lower limit of allowed y-regime
c     z_maxgeom : Upper limit of allowed z-regime
c     z_mingeom : Lower limit of allowed z-regime
      double precision x_maxgeom,x_mingeom,y_maxgeom,
     &                 y_mingeom,z_maxgeom,z_mingeom
      parameter ( x_maxgeom =  10.e0     ,
     &            x_mingeom = -10.e0     ,
     &            y_maxgeom =  10.e0     ,
     &            y_mingeom = -10.e0     ,
     &            z_maxgeom =  10.e0     ,
     &            z_mingeom = -10.e0     )

c     SPH algorithm for particle approximation (pa_sph)
c     pa_sph = 1 : (e.g. (p(i)+p(j))/(rho(i)*rho(j))
c              2 : (e.g. (p(i)/rho(i)**2+p(j)/rho(j)**2)
      integer pa_sph
      parameter(pa_sph = 2)

c     Nearest neighbor particle searching (nnps) method
c     nnps = 1 : Simplest and direct searching
c            2 : Sorting grid linked list
c            3 : Tree algorithm
      integer nnps
      parameter(nnps = 1 )

c     Smoothing length evolution (sle) algorithm
c     sle = 0 : Keep unchanged,
c           1 : h = fac * (m/rho)^(1/dim)
c           2 : dh/dt = (-1/dim)*(h/rho)*(drho/dt)
c           3 : Other approaches (e.g. h = h_0 * (rho_0/rho)**(1/dim) )

      integer sle
      parameter(sle = 0)

c     Smoothing kernel function
c     skf = 1, cubic spline kernel by W4 - Spline (Monaghan 1985)
c         = 2, Gauss kernel   (Gingold and Monaghan 1981)
c         = 3, Quintic kernel (Morris 1997)
      integer skf
      parameter(skf = 1)

3 楼

dingsheng1206 [专家分：0] 发布于 2010-09-12 22:28:00

[quote]
include 'param.inc'
param.inc 没上传啊[em15][/quote]
c     Switches for different senarios

c     summation_density = .TRUE. : Use density summation model in the code,
c                        .FALSE.: Use continuiity equation
c     average_velocity = .TRUE. : Monaghan treatment on average velocity,
c                       .FALSE.: No average treatment.
c     config_input = .TRUE. : Load initial configuration data,
c                   .FALSE.: Generate initial configuration.
c     virtual_part = .TRUE. : Use vritual particle,
c                   .FALSE.: No use of vritual particle.
c     vp_input = .TRUE. : Load virtual particle information,
c               .FALSE.: Generate virtual particle information.
c     visc = .true. : Consider viscosity,
c           .false.: No viscosity.
c     ex_force =.true. : Consider external force,
c               .false.: No external force.
c     visc_artificial = .true. : Consider artificial viscosity,
c                      .false.: No considering of artificial viscosity.
c     heat_artificial = .true. : Consider artificial heating,
c                      .false.: No considering of artificial heating.
c     self_gravity = .true. : Considering self_gravity,
c                    .false.: No considering of self_gravity
c     nor_density =  .true. : Density normalization by using CSPM,
c                    .false.: No normalization.
      logical summation_density, average_velocity, config_input,
     &        virtual_part, vp_input, visc, ex_force, heat_artificial,
     &        visc_artificial, self_gravity, nor_density
      parameter ( summation_density  = .true. )
      parameter ( average_velocity  = .false. )
      parameter ( config_input  = .false. )
      parameter ( virtual_part  = .false. )
      parameter ( vp_input  = .false.  )
      parameter ( visc  = .false.  )
      parameter ( ex_force  = .false.)
      parameter ( visc_artificial  = .true. )
      parameter ( heat_artificial  = .false. )
      parameter ( self_gravity  = .false. )
      parameter ( nor_density  = .false. )

c     Symmetry of the problem
c     nsym = 0 : no symmetry,
c          = 1 : axis symmetry,
c          = 2 : center symmetry.
      integer    nsym
      parameter ( nsym = 0)

c     Control parameters for output
c     int_stat = .true. : Print statistics about SPH particle interactions.
c                        including virtual particle information.
c     print_step: Print Timestep (On Screen)
c     save_step : Save Timestep    (To Disk File)
c     moni_particle: The particle number for information monitoring.
      logical int_stat
      parameter ( int_stat = .true. )
      integer print_step, save_step, moni_particle
      parameter ( print_step = 100 ,
     &            save_step = 500,
     &            moni_particle = 1600   )

      double precision pi
      parameter ( pi = 3.14159265358979323846 )

c     Simulation cases
c     shocktube = .true. : carry out shock tube simulation
c     shearcavity = .true. : carry out shear cavity simulation
      logical shocktube, shearcavity
      parameter ( shocktube  = .true. )
      parameter ( shearcavity  = .false. )

4 楼

jstzhurj [专家分：4680] 发布于 2010-09-12 23:09:00

cvf6.6 编译运行没有任何错误。

5 楼

dingsheng1206 [专家分：0] 发布于 2010-09-13 08:39:00

[quote]
cvf6.6 编译运行没有任何错误。[/quote]
谢谢，我的6.5，我再看看吧，

6 楼

dingsheng1206 [专家分：0] 发布于 2010-09-13 09:01:00

[quote]
cvf6.6 编译运行没有任何错误。[/quote]
你好，谢谢你，我刚才重新调了一下好了，谢谢你的帮助，

我来回复

您尚未登录，请登录后再回复。点此登录或注册

主题：数组越界，跪求大大帮忙，

回复列表（共6个回复）

我来回复

程序员工具箱 new

代码片段

本版新帖

主题：数组越界，跪求大大帮忙，

回复列表 （共6个回复）

我来回复

程序员工具箱 new

代码片段

本版新帖

回复列表（共6个回复）