anchor.diff:            end do
          end do
       end do
  +
       deallocate(x,w,p1,p2,p3,paux)
       return
     end subroutine gaussedges
  
positive.code:               sinal=-1.0D0
             endif
   
            !Determinando os valores para o node i 
             xx = mesh%v(i)%p(1)
             yy = mesh%v(i)%p(2)
             zz = mesh%v(i)%p(3)
  
negative.code:     deallocate(p1,p2,p3)
      return  
    end subroutine flux_olg
  
  
  
      subroutine flux_gas(nodes,mesh,z,time, uedges)  
      !----------------------------------------------------------------------------------------------
      !    Calculando o fluxo nas arestas
  