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:     mask2 = .true.
  
      do i = 1, n1
          mask1(x1(2,i)) = .false.
      enddo
          
      do i = 1, n2
          mask2(x2(2,i)) = .false.
      enddo
  
      do m = 1, order
  
negative.code:                 do j_1 = 1, nj
  
                      l2dist = atomic_distl2(x1(i,i_1,:,:), x2(j,j_1,:,:), n1(i), n2(j), &
                          & ksi1(i,i_1,:), ksi2(j,j_1,:), sinp1(i,i_1,:,:,:), sinp2(j,j_1,:,:,:), &
                          & cosp1(i,i_1,:,:,:), cosp2(j,j_1,:,:,:), &
                          & t_width, r_width, c_width, d_width, cut_distance, order, pd)
  
                      l2dist = selfl21(i,i_1) + selfl22(j,j_1) - 2.0d0 * l2dist * pd(int(x1(i,i_1,2,1)),int(x2(j,j_1,2,1)))
  
                      if (abs(l2dist) < eps) l2dist = 0.0d0
  
  