anchor.diff:      real*8, intent(out) :: cmax, cova
   
       ! internal variables
  -    real*8 ::    hsqd, maxnst, h, hr
  -    integer :: ir, is, ist, istart
  +    real*8 ::    hsqd, h, hr
  +    integer :: ir, is, ist, istart, maxnst
       
       !parameters
       real*8 :: DEG2RAD,EPSLON,PI,PMX
  
positive.code:     implicit none
  
  
  	!external references
  	real*8,external :: sqdist
  
      !input
      integer, intent(in) :: nxsup, nysup, nzsup
  
negative.code:     
    
      ! internal 
      integer, dimension (9) :: idrif
      
      
      ! constants
      real*8 :: EPSLON=0.000001,PMX = 999.0
      integer :: MAXROT
  