guniqi4.F

Go to the documentation of this file.

! Find unique integers in a sorted list
      SUBROUTINE guniqi4 (N,ORD,A)
      USE parkind1, ONLY : jpim, jprb
!      IMPLICIT INTEGER (A-Z)
      IMPLICIT NONE

      INTEGER(KIND=4) :: N
      INTEGER(KIND=4) :: A(n)
      INTEGER(KIND=4) :: ORD(n)
      INTEGER(KIND=4) :: I, J

      IF (n == 0) RETURN

      j = 1
      ord(1) = 1

      DO i = 2, n
        IF (a(i-1) /= a(i)) THEN
          j = j + 1
          ord(j) = i
        ENDIF
      ENDDO

      n = j

      ENDSUBROUTINE guniqi4