subroutine srhel(u,v,p,nvp,cx,cy,ip1,ip2,isrh)
c     This routine calculates the storm relative helicity
c     integrated from pressure level p(ip1) to p(ip2)
c
c     Input:
c        u(p) - zonal wind component (m/s)
c        v(p) - meridional wind component (m/s)
c        p    - Pressure levels (low pressure to high pressure)
c
c        nvp  - No of pressure levels
c        cx   - zonal component of storm motion (m/s)
c        cy   - meridional component of storm motion (m/s)
c        ip1  - Index of first pressure level for integration
c        ip2  - Index of second pressure level for integration
c        isrh  - Storm relative helicity (m**2/s**2) * 10
c               (missing value = 9999)
c   
      dimension u(nvp),v(nvp),p(nvp)
c
      isrh = 9999
c   
c     Check indices
      if (ip1 .ge. ip2) return
      if (ip1 .lt. 1 .or. ip1 .gt. nvp) return
      if (ip2 .lt. 1 .or. ip2 .gt. nvp) return
c
      srh = 0.0
      do i=ip1,ip2-1
         srh = srh + (0.5*(u(i+1)+u(i)) - cx)*(v(i+1)-v(i))
     +             - (0.5*(v(i+1)+v(i)) - cy)*(u(i+1)-u(i))
      enddo
c
      isrh = nint(1.944*1.944*srh)
c
      return
      end