anchor.diff:  
       ! Number of quadrature points
       nquad = nint((order)/2.0D0)
  +    
  +    if (method=='G') then
  +      nquad = 1
  +    end if
   
       ! Maximum of neighbors in a Voronoi cell
       maxnnb = maxval(mesh%v(:)%nnb)
  
positive.code: !          phi%name=trim(adjustl(trim(transpname))//"_phi_t_"//trim(adjustl(trim(atime))))    
         end do
      end do
  !    print*, erro_Linf!, 'erro_Linf'
  !    print*, dsqrt (temp/contador)! , 'erro_L2'
  !    print*, contador 
  
negative.code:     deallocate(p1,p2,p3)
      return  
    end subroutine flux_olg
  
  
  
      subroutine flux_gas(nodes,mesh,z,time, uedges)  
      !----------------------------------------------------------------------------------------------
      !    Calculando o fluxo nas arestas
  