anchor.diff:         !call interpolation_g(mesh)
          !call flux_edges_g(nodes,mesh,0,0.0D0)  
   
  -      do i=1,nodes
  +       do i=1,nodes
             node(i)%phi_exa=node(i)%phi_new2
          enddo
   
  
positive.code:   
    
    
      subroutine flux_gas(nodes,mesh,z,time)  
      !----------------------------------------------------------------------------------------------
      !    Calculando o fluxo nas arestas
      !----------------------------------------------------------------------------------------------
  
      implicit none 
      integer,intent(in)    :: nodes
      type(grid_structure),intent(in)      :: mesh
      integer,intent(in),optional   :: z
      real(r8),intent(in),optional  :: time
      integer(i4)   :: i
      integer(i4)   :: j
      integer(i4)   :: k
  
negative.code:           node(i)%phi_exa=node(i)%phi_new2
         enddo
  
  
         !Avanço temporal 
         time=0.0D0
  
         do j=0,node(0)%ntime+1
            !Atualizando a solucao
            do i=1,nodes
               node(i)%phi_old=node(i)%phi_new2
            enddo
   
            !Plotando Inicial e Final
  !          if(j==0.or.j==node(0)%ntime)then   
  !             do i=1,nodes
  