!
SUBROUTINE findvasy(vasy, ri,ichanl)
USE Numeric_Kinds_Module
!==============================================================================
!this subroutine finds the value of the potential at ri
!==============================================================================
IMPLICIT NONE
INTEGER :: i, ichanl
REAL(dp), PARAMETER :: tol=1.d-6
REAL(dp) :: vasy, ri, v0, r1, r2, vl, vr

 v0=0.d0

 CALL diatomic_v(1.d30, vasy,ichanl)   
    DO i=0,10                         
       ri=2.d0**i 
       CALL diatomic_v(ri, v0,ichanl)
       IF(ABS(v0-vasy).le.tol)THEN     
          r1=2.d0**(i-1)
          CALL diatomic_v(r1, vl,ichanl)
          r2=2.d0**i
          CALL diatomic_v(r2, vr,ichanl)
          CALL findri(ri, v0, r1, vl, r2, vr, vasy,ichanl)
          EXIT
       ENDIF
    ENDDO

    RETURN
    END SUBROUTINE Findvasy