*dk h2odp
      subroutine h2odp(dpl,rij)
C
C     H2O - dipolmoment
C     CALCULATED BY STAEMMLER + PALMA;  FIT BY SCHINKE
C
C     ABSTAENDE IN ANGSTROM,  DIPOLE MOMENT IN DEBYE
C
      implicit real*8 (a-h,o-z)
C
      dimension r(3),re(3),alpha(3),c(100),phi(100)
      dimension rij(3)
      real*8 angtoau, autodby
      data angtoau/0.529177249d0/, autodby/2.541750d0/
C
C     EINLESEN DER DATEN
C
      data norder/ 50/
      data re/      1.0000d0,      1.0000d0,      1.0000d0/
      data alpha/      0.2000d0,      0.2000d0,      0.2000d0/
      data gamma/      0.0000d0/
      data c(  1)/  0.99074524d+00/
      data c(  2)/ -0.77176762d+00/
      data c(  3)/  0.69985137d+00/
      data c(  4)/ -0.63681282d+01/
      data c(  5)/  0.64115665e-01/
      data c(  6)/ -0.12161895d+00/
      data c(  7)/  0.76708290d+01/
      data c(  8)/ -0.62791600d+01/
      data c(  9)/ -0.74103365e-01/
      data c( 10)/  0.52642150d+01/
      data c( 11)/  0.10310915d+00/
      data c( 12)/ -0.61293702d+01/
      data c( 13)/  0.55889839d+01/
      data c( 14)/  0.12266109d+03/
      data c( 15)/ -0.27575632d+00/
      data c( 16)/ -0.90865798d+02/
      data c( 17)/  0.33107891d+02/
      data c( 18)/ -0.90123009d+00/
      data c( 19)/ -0.26998948d+01/
      data c( 20)/ -0.19625944d+00/
      data c( 21)/ -0.23340421d+01/
      data c( 22)/  0.26326061d+02/
      data c( 23)/ -0.26909940d+03/
      data c( 24)/ -0.48099992d+01/
      data c( 25)/  0.17902869d+03/
      data c( 26)/ -0.34551909d+01/
      data c( 27)/ -0.59593047d+02/
      data c( 28)/  0.10379498d+03/
      data c( 29)/  0.11567520d+02/
      data c( 30)/ -0.35853770d+02/
      data c( 31)/  0.50357697d+01/
      data c( 32)/  0.49106000d+01/
      data c( 33)/  0.14507240d+01/
      data c( 34)/ -0.45981721d+00/
      data c( 35)/  0.18736500d+03/
      data c( 36)/  0.20127856d+03/
      data c( 37)/ -0.52398738d+03/
      data c( 38)/  0.52610133d+03/
      data c( 39)/  0.43101052d+02/
      data c( 40)/ -0.18386779d+02/
      data c( 41)/ -0.17411415d+03/
      data c( 42)/ -0.23899011d+02/
      data c( 43)/  0.75039407d+02/
      data c( 44)/  0.52807753d+02/
      data c( 45)/ -0.31026072d+01/
      data c( 46)/ -0.21710015d+01/
      data c( 47)/ -0.20880270d+01/
      data c( 48)/ -0.51271737d+01/
      data c( 49)/  0.22135686d+00/
      data c( 50)/  0.14258218d+00/
C
c  The following conventions exist for the internuclear distances:
c     r(1) = r(OH),  r(2) = r(HH),  r(3) = r(OH)
c  Convert from aph to staemmler-schinke coordinate convention
c  Convert distances from au to angstroms
      r(1) = rij(3)*angtoau
      r(2) = rij(1)*angtoau
      r(3) = rij(2)*angtoau
C
      call funcdp(r,phi,norder,re,alpha)
C
      dpl=0.0
      do 1 i=1,norder
 1    dpl=dpl+c(i)*phi(i)
C
      dpl=dpl/autodby
C
      return
      end
