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:     integer(i4) :: nlines
      integer(i4) :: ncolumns
      real(r8):: time
  
      !Numero de vizinhos e vizinhos de cada no
      integer(i4),allocatable   :: nbsv(:,:)   
  
negative.code: 
    !Logical for plots or not
    logical:: plots
  
  
    !Logical for moist shallow water model
    logical:: moistswm
  
    !Logical for moist shallow water model initial condition
    integer:: moistswm_ic
  
    ! Reconstruction interpolation method for high order adv scheme in the moist swm
    character (len=8) :: reconadvmtd = "lsqhxe"
    !character (len=8) :: reconadvmtd = "lsqtrc"
  
    !Plotsteps - will plot at every plotsteps timesteps
    integer (i4):: plotsteps
  