anchor.diff:      mask2 = .true.
   
       do i = 1, n1
  -        mask1(x1(2,i)) = .false.
  +        mask1(int(x1(2,i))) = .false.
       enddo
           
       do i = 1, n2
  -        mask2(x2(2,i)) = .false.
  +        mask2(int(x2(2,i))) = .false.
       enddo
   
       do m = 1, order
  
positive.code:     integer :: i, m, p1, p2
  
      double precision :: angular 
      double precision :: angular2
  
      double precision :: dcos1, dcos2, dsin1, dsin2
      double precision :: maxgausdist2
  
      integer :: pmax1
      integer :: pmax2
  
      double precision :: inv_width, diff, dist
      double precision :: c_width2, r_width2, r2
  
      double precision, dimension(order) :: s
  
negative.code: 
              if (r2 < maxgausdist2) then
  
                  d = exp(r2 * inv_width )  * ksi1(m_1) * ksi2(m_2) * &
                      & pd(int(x1(2,m_1)), int(x2(2,m_2)))
  
                  angular = a0 * a0
                  
  