c---------------------------------------------------
c parameters for qdiff4.f and subroutines
c---------------------------------------------------
	parameter (nclist = 1500)
	parameter (nrlist = 1500)
	parameter (nrmax = 1500)
c nrmax= maximum number of entries in radius file 
	parameter (ncmax = 1500)
c ncmax= maximum number of entries in charge file 
	parameter (ncrgmx = 20000)
c ncrgmx = maximum number of charged atoms 
	parameter (ngrid = 65)
c ngrid= maximum grid size 
	parameter (ngp = 274626)
c ngp = ngrid**3 = number of grid points 
	parameter (nhgp = 137313)
c nhgp = half the previous number 
	parameter (natmax = 25000)
c natmax = maximum number of atoms 
	parameter (ngcrg = 80000)
c ngcrg = maximum number of grid points that may be assigned charge (was 4226)
	parameter (nbgp = 8000)
c nbgp = number of points at one box boundary = ngrid**2
	parameter (nsp = 21130)
c nsp = maximum number of dielectric boundary points, appox.= 5*nbgp 
c---------------------------------------------------
c---------------------------------------------------
      dimension atnam(nrmax)		
c	atom names for radii table
      dimension rnam(nrmax)		
c	residue names for radii table
      dimension radt(nrmax)		
c	radius table 
      dimension irlink(nrlist)	
c	links for radius hash table
      dimension irnumb(nrlist)	
c	radii id numbers in hash table
      dimension catnam(ncmax)		
c	atom names   for charge table
      dimension cchn(ncmax)		
c	chain names for charge table
      dimension crnam(ncmax)		
c	residue names for charge table
      dimension crnum(ncmax)		
c	residue number for charge table
      dimension chrgvt(ncmax)		
c	charge table 
      dimension iclink(nclist)	
c	links for charge hash table
      dimension icnumb(nclist)	
c	charge entry id numbers in hash table
	dimension atmcrg(4,ncrgmx)	
c	atom charges for boundary condition
 	dimension phimap(ngrid,ngrid,ngrid),phimap1(nhgp),phimap2(nhgp)
	dimension phimap3(ngp)
 	integer iepsmp(ngrid,ngrid,ngrid,3),iepsmp2(ngrid,ngrid,ngrid,3)
	integer idebmap(ngrid,ngrid,ngrid)
	integer iatmap(ngrid,ngrid,ngrid),iatmap2(ngrid,ngrid,ngrid)
	dimension cgbp(5,10000),gval(ngcrg)
	dimension oldmid(3),oldmid1(3)
c---------------------------------------------------
c---------------------------------------------------
      character*1 chn,cchn
      character*3 crnam,rnam,res
      character*4 rnum,crnum
      character*6 atnam,catnam,atm
c---------------------------------------------------
c---------------------------------------------------
	common
     &	/link/  irlink,irnumb,iclink,icnumb,irtot,ictot
     &	/name/  atnam,rnam,catnam,cchn,crnam,crnum
     &	/value/ radt,chrgvt
     &	/maps/  phimap,phimap1,phimap2,phimap3,atmcrg
     &  /imaps/ iepsmp,iepsmp2,idebmap,iatmap,iatmap2
     &	/scale/ scale,oldmid,scale1,oldmid1,igrid,ibc,gten,cgbp,gval,rmmin,rmmax
     &	/log/   ihs,isen,isch
c---------------------------------------------------
c---------------------------------------------------