c
c
c     ###################################################
c     ##  COPYRIGHT (C)  1992  by  Jay William Ponder  ##
c     ##              All Rights Reserved              ##
c     ###################################################
c
c     ###############################################################
c     ##                                                           ##
c     ##  module polar  --  polarization & induced dipole moments  ##        
c     ##                                                           ##
c     ###############################################################
c
c
c     npolar    total number of polarizable sites in the system
c     ipolar    number of the atom for each polarizable site
c     jpolar    index into polarization parameter matrix for each atom
c     polarity  dipole polarizability for each atom site (Ang**3)
c     thole     Thole polarization damping value for each atom
c     tholed    Thole direct polarization damping value for each atom
c     pdamp     value of polarizability scale factor for each atom
c     thlval    Thole damping parameter value for each atom type pair
c     thdval    alternate Thole direct damping value for atom type pair
c     udir      direct induced dipole components for each atom site
c     udirp     direct induced dipoles in field used for energy terms
c     udirs     direct GK or PB induced dipoles for each atom site
c     udirps    direct induced dipoles in field used for GK or PB energy
c     uind      mutual induced dipole components for each atom site
c     uinp      mutual induced dipoles in field used for energy terms
c     uinds     mutual GK or PB induced dipoles for each atom site
c     uinps     mutual induced dipoles in field used for GK or PB energy
c     uexact    exact SCF induced dipoles to full numerical precision
c     douind    flag to allow induced dipoles at each atom site
c
c
      module polar
      implicit none
      integer npolar
      integer, allocatable :: ipolar(:)
      integer, allocatable :: jpolar(:)
      real*8, allocatable :: polarity(:)
      real*8, allocatable :: thole(:)
      real*8, allocatable :: tholed(:)
      real*8, allocatable :: pdamp(:)
      real*8, allocatable :: thlval(:,:)
      real*8, allocatable :: thdval(:,:)
      real*8, allocatable :: udir(:,:)
      real*8, allocatable :: udirp(:,:)
      real*8, allocatable :: udirs(:,:)
      real*8, allocatable :: udirps(:,:)
      real*8, allocatable :: uind(:,:)
      real*8, allocatable :: uinp(:,:)
      real*8, allocatable :: uinds(:,:)
      real*8, allocatable :: uinps(:,:)
      real*8, allocatable :: uexact(:,:)
      logical, allocatable :: douind(:)
      save
      end