anchor.diff:  
     !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
  
positive.code:        !call interpolation_g(mesh)
         !call flux_edges_g(nodes,mesh,0,0.0D0)  
  
        do i=1,nodes
            node(i)%phi_exa=node(i)%phi_new2
         enddo
  
  
negative.code:              call flux_gas(nodes,mesh,k-1,time)
               call rungekutta(nodes,mesh,k) 
            enddo
            
         
            time=time+node(0)%dt
            print*, time
         enddo
  