SUBROUTINE MorsePotential(rau, Potential)
USE Numeric_Kinds_Module
IMPLICIT NONE
REAL(KIND=WP_Kind) r, rau, Potential
REAL(KIND=WP_Kind), PARAMETER:: De=333.758_WP_Kind, re=4.17005_WP_Kind, alpha=2.5_WP_Kind
!Convert the radial distance from bohr to Angstroms
r=rau*0.529176
Potential=De*(Exp(-(r-re)/alpha)-1.d0)**2
Potential=Potential/2.1947e+05

RETURN
ENDSUBROUTINE MorsePotential