ccommc.F

Go to the documentation of this file.

! Look for common elements in two sorted lists
      SUBROUTINE ccommc (KNC, KCOM1, KCOM2, KN1, CD1, KN2, CD2)
      USE parkind1, ONLY : jpim
      IMPLICIT NONE
      
      INTEGER(KIND=JPIM) :: KNC, KN1, KN2
      INTEGER(KIND=JPIM) :: KCOM1 (knc), KCOM2 (knc)

      CHARACTER(LEN=*) :: CD1 (kn1), CD2 (kn2)
      INTEGER(KIND=JPIM) :: I1, I2, IC

      i1 = 1
      i2 = 1
      ic = 0

      DO WHILE ((i1 <= kn1) .AND. (i2 <= kn2))
        IF (cd1(i1) == cd2(i2)) THEN
          ic = ic + 1
          kcom1(ic) = i1
          kcom2(ic) = i2
          i1 = i1 + 1
          i2 = i2 + 1
        ELSEIF (cd1(i1) > cd2(i2)) THEN
          i2 = i2 + 1
        ELSEIF (cd1(i1) < cd2(i2)) THEN
          i1 = i1 + 1
        ENDIF
      ENDDO
 
      knc = ic

      ENDSUBROUTINE ccommc