C
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 ( ywater_level = 91.75d0, 
     *     water_density = 1000.0d0,
     *     gravity = 9.81d0)
C
      dimension curCoords(nblock,ndim), 
     *     velocity(nblock,ndim),
     *     dircos(nblock,ndim,ndim), 
     *     value(nblock)
      character*80 sname
*
C
C     USER SUBROUTINE  TO APPLY THE WATER 
C     HYDROSTATIC PRESSURE ON KOYNA CONCRETE DAM. 
C     
      do k = 1, nblock
        depth = ywater_level - curCoords(k,2)
        value(k) = amplitude * 
     *       depth * gravity * water_density 
     
      end do
*     
      return
      end