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

CALL DiatomicPES_SplineFit(r,v0, v1, v2)
CALL CenTerm(r, jvalue, mu, cent0, cent1, cent2)
v0=v0+cent0
v1=v1+cent1
v2=v2+cent2

RETURN
END SUBROUTINE PotDiatomic