anchor.diff:      integer :: i, m, p1, p2
   
       double precision :: angular 
  -    double precision :: angular2
  -
  -    double precision :: dcos1, dcos2, dsin1, dsin2
  +    ! double precision :: angular2
  +
  +    ! double precision :: dcos1, dcos2, dsin1, dsin2
       double precision :: maxgausdist2
   
       integer :: pmax1
       integer :: pmax2
   
  -    double precision :: inv_width, diff, dist
  +    double precision :: inv_width !, diff, dist
       double precision :: c_width2, r_width2, r2
   
       double precision, dimension(order) :: s
  
positive.code:     double precision :: temp, sin1_temp, cos1_temp
  
  
      pmax1 = maxval(x1(2,:))
      pmax2 = maxval(x2(2,:))
  
      allocate(mask1(pmax1))
      allocate(mask2(pmax2))
  
negative.code:     mask2 = .true.
  
      do i = 1, n1
          mask1(int(x1(2,i))) = .false.
      enddo
          
      do i = 1, n2
          mask2(int(x2(2,i))) = .false.
      enddo
  
      do m = 1, order
  