rsum.F

Go to the documentation of this file.

      FUNCTION rsum(N,X,INCX)
      USE parkind1, ONLY : jpim, jprb
      IMPLICIT NONE
      INTEGER(KIND=JPIM) :: N
      INTEGER(KIND=JPIM) :: INCX
      REAL(KIND=JPRB) :: X(n*incx)
      REAL(KIND=JPRB) :: RSUM
      INTEGER(KIND=JPIM) :: I, J
      REAL(KIND=JPRB) :: ZSUM
!
      zsum=0.0_jprb
      i=1
      DO j=1,n
        zsum=zsum+x(i)
        i=i+incx
      ENDDO
      rsum=zsum
      RETURN
      ENDFUNCTION rsum