求助， [复制链接]

各位楼主，请帮我一会!
我在学习ABAQUS软件，我按照ABAQUS软件在道路工程中的应用 动态模拟 但做到 动态荷载分布子程序 （vdload30.for）就做不了。求助帮我看一下这个程序那地方不对。多谢

C User subroutine VDLOAD
      subroutine vdload (
C Read only (unmodifiable) variables -
     *      nblock, ndim, stepTime, totalTime,  
     *      amplitude, curCoords, velocity, dircos,  
     *      jltyp, sname,
C Write only (modifiable) variable -
     *      value )
C
      include 'vaba_param.inc'
      parameter (zini=4.5,vel=30,dlen=0.166666667,pressure=0.7d6)
C
      dimension curCoords(nblock,ndim),  
     *      velocity(nblock,ndim),

*      dircos(nblock,ndim,ndim),  
     *      value(nblock)
      character*80 sname
C      -------------------------------------------------------------
C      distan为steptime时间内荷载移动的距离（车速vel为30km/h）；
C      zmax和zmin分别为荷载的上下边界（两者相距dlen为0.166666667m）；
C      -------------------------------------------------------------
      distan=vel*stepTime
      zc=zini-distan
      zmax=zc
      zmin=zmax-dlen
      do 100 k=1, nblock
       if(curCoords(k,3).lt.zmax.and.curCoords(k,3).ge.zmin) then
             value(k)=pressure
       else
             value(k)=0
      end if
100  continue
      return
      end