C
C
C     This code is distributed under the terms and conditions of the
C     CCP4 licence agreement as `Part 2' (Annex 2) software.
C     A copy of the CCP4 licence can be obtained by writing to the
C     CCP4 Secretary, Daresbury Laboratory, Warrington WA4 4AD, UK.
C
c+++ sphhrm.fh +++
c
c Spherical harmonic parameters
c  nprsph                     total number of parameters over all runs,
c                                = 0 none
c  nrnsph(maxrun)             number of spherical harmonic parameters 
c                              for this run, = 0 none
c                              The first constant is not counted
c  isphpt(maxrun)             pointer to parameter set, to allow runs to share
c  lmxsph(maxrun)             maximum l order 
c  lodsph(maxrun)             = 1 all = 2 even orders only (centrosymmetric)
c
c  sphclm(maxspt,maxrun)      coefficients (excluding constant term)
c
c  nplphi, npltheta           number of points for polar phi & theta
c  dplphi, dpltheta           interval on phi, theta (degrees)
c  nsphith(mxstheta, mxsphi, maxrun)  number of contributions to this bin
c  avphith(mxstheta, mxsphi, maxrun)  sum|average of this bin
c  sphmnmx(2,maxrun)          minimum & maximum values of correction
c                              applied in this run

c 
      common /sphhrm/
     $   avphith(mxstheta, mxsphi, maxrun),
     $   sphclm(maxspt,maxrun), sphmnmx(2,maxrun),
     $   dplphi, dpltheta,
     $   nprsph, nrnsph(maxrun), isphpt(maxrun), lmxsph(maxrun),
     $   lodsph(maxrun), mhpole(maxrun), nplphi, npltheta,
     $   nsphith(mxstheta, mxsphi, maxrun)
c
      real  avphith, sphmnmx
      real  sphclm, dplphi, dpltheta
      integer nprsph, nrnsph, isphpt, lmxsph, lodsph, mhpole,
     $   nplphi, npltheta, nsphith
c
      save /sphhrm/
c
c--- sphhrm.fh ---