c
      subroutine h2omsm(pot,r)
c
      implicit real*8 (a-h,o-z)
      dimension r(3)
      real*8 angtoau, cmm1toau
      data angtoau/0.529177249d0/,cmm1toau/4.5563352672d-6/
      external vsm2
c
      s3 = r(1)*angtoau
      s1 = r(2)*angtoau
      s2 = r(3)*angtoau
c
      pot = cmm1toau*vsm2(s3,s1,s2)
c
      return
      end