SUBROUTINE PotDiatomic(r, hstep, v0, v1, v2, jvalue, mu)
USE Numeric_Kinds_Module
USE Numbers_Module
!USE OneDim_Parms_Module
USE FileUnits_Common_Module
USE FileUnits_OneDim_Module
!USE CriticalDist_Module
USE QState_Module
USE DiatomicPot_Module
IMPLICIT NONE
CHARACTER(LEN=4), PARAMETER :: PESName='Vli2'
INTEGER jvalue
REAL(KIND=WP_Kind) r, mu, hstep, h, v0, v1, v2, cent0, cent1, cent2, pot, potph, potmh, rvals(3)

h=1.0d-3
rvals(1)=1.d+6
rvals(2)=1.d+6
rvals(3)=r
CALL Surface(pot,rvals)

rvals(3)=r+h
CALL surface(potph,rvals)

rvals(3)=r-h
CALL surface(potmh,rvals)

CALL CenTerm(r, jvalue, mu, cent0, cent1, cent2)
v0=pot+cent0
v1=(potph-potmh)/h+cent1
v2=(potph-Two*pot+potmh)/(Two*h**2)+cent2

RETURN
END SUBROUTINE PotDiatomic