isrchfltpv.F

Go to the documentation of this file.

      FUNCTION isrchfltpv(N,ARRAY,INC,TARGET,NBITER)
      USE parkind1, ONLY : jpim, jprb
      IMPLICIT NONE
      INTEGER(KIND=JPIM) :: ISRCHFLTPV
      INTEGER(KIND=JPIM) :: N
      INTEGER(KIND=JPIM) :: INC
      REAL(KIND=JPRB) :: ARRAY(1+inc*(n-1))
      REAL(KIND=JPRB) :: TARGET
      INTEGER(KIND=JPIM) :: NBITER
      INTEGER   :: I
!
      IF (n.LE.0) THEN
        isrchfltpv=0
        RETURN
      ELSE
        IF (nbiter.EQ.3) THEN
        DO 3 i=1,n-nbiter
          IF ( array(1+inc*(i-1)).LT.TARGET ) THEN
             IF ( array(1+inc*(i)).LT.TARGET ) THEN
                IF ( array(1+inc*(i+1)).LT.TARGET ) THEN
                   IF ( array(1+inc*(i+2)).LT.TARGET ) THEN
                      isrchfltpv=i
                      RETURN
                   ENDIF
                ENDIF
             ENDIF
          ENDIF
 3      CONTINUE
        DO i=n-nbiter+1,n
          IF ( array(1+inc*(i-1)).LT.TARGET ) THEN
            isrchfltpv=i
            RETURN
          ENDIF
        ENDDO
        ELSEIF (nbiter.EQ.2) THEN
        DO 2 i=1,n-nbiter
          IF ( array(1+inc*(i-1)).LT.TARGET ) THEN
             IF ( array(1+inc*(i)).LT.TARGET ) THEN
                IF ( array(1+inc*(i+1)).LT.TARGET ) THEN
                      isrchfltpv=i
                      RETURN
                ENDIF
             ENDIF
          ENDIF
 2      CONTINUE
        DO i=n-nbiter+1,n
          IF ( array(1+inc*(i-1)).LT.TARGET ) THEN
            isrchfltpv=i
            RETURN
          ENDIF
        ENDDO
        ELSEIF (nbiter.EQ.1) THEN
        DO 1 i=1,n-nbiter
          IF ( array(1+inc*(i-1)).LT.TARGET ) THEN
             IF ( array(1+inc*(i)).LT.TARGET ) THEN
                      isrchfltpv=i
                      RETURN
             ENDIF
          ENDIF
 1      CONTINUE
        DO i=n-nbiter+1,n
          IF ( array(1+inc*(i-1)).LT.TARGET ) THEN
            isrchfltpv=i
            RETURN
          ENDIF
        ENDDO
        ELSEIF (nbiter.EQ.0) THEN
        DO 4 i=1,n-nbiter
          IF ( array(1+inc*(i-1)).LT.TARGET ) THEN
                      isrchfltpv=i
                      RETURN
          ENDIF
 4      CONTINUE
        ELSE
        CALL abor1('ISRCHFLTPV: ERROR WITH NBITER VALUE')
        ENDIF
        isrchfltpv=n+1
      ENDIF
      ENDFUNCTION isrchfltpv