anchor.diff:      real(r8),allocatable  :: p2(:)      
       real(r8),allocatable  :: p3(:)      
       real(r8),allocatable  :: vn(:)      
  +    real(r8) :: dir_i(1:3), dir_j(1:3), rot_dir_i(1:3), rot_dir_j(1:3)
   
       allocate (p1(3),p2(3),p3(3),vn(3))
       FEPS = 1.0D-8
  
positive.code: 
      ! Number of quadrature points
      nquad = nint((order)/2.0D0)
  
      ! Maximum of neighbors in a Voronoi cell
      maxnnb = maxval(mesh%v(:)%nnb)
  
negative.code:     aaxx = 0.0D0
      erro_Linf = 0.0D0
      contador = 0  
  
      do i = 1,nodes
           node(i)%S(z)%flux = 0.0D0
           jend = node(i)%ngbr(1)%numberngbr
  