anchor.diff:      return
     end function f
   
  -  function velocity(p, time, uedges, mesh)
  +  function velocity(p, time)
       !-----------------------------------
       !  V - velocity in Cartesian coordinates
       !
  
positive.code:           p2=node(i)%edge(n)%xyz2(2,:) 
            do s=1,diml
               !Calculando o produto interno entre o vetor velocidade e o vetor normal unitario a edge 
               dot=dot_product(velocity(node(i)%G(s)%lpg(n,1:3), time, uedges, mesh),node(i)%G(s)%lvn(n,1:3))
               !Determinando os coeficientes da matriz de rotação para o node i 
               xx=mesh%v(i)%p(1)
               yy=mesh%v(i)%p(2)
  
negative.code:           allocate(node(i)%G(1))
            allocate(node(i)%G(1)%lwg(2*ngbr1))
            allocate(node(i)%G(1)%lpg(2*ngbr1,3))
            allocate(node(i)%G(1)%velocity_quadrature(2*ngbr1))
            allocate(node(i)%G(1)%lvn(2*ngbr1,3))
            allocate(node(i)%G(1)%upwind_donald(2*ngbr1))
  
  