anchor.diff:      
       ! High order advection scheme variable
       character (len=6):: advmtd
  -
  -    ! Reconstruction interpolation method for high order adv scheme
  -    character (len=8) :: reconadvmtd = "lsqhxe"
   
   contains
   
  
positive.code:     type(scalar_field), intent(in):: h  !General
  
      !Velocities (defined on edges - only normal component)
      type(scalar_field), intent(in):: u  !General
  
      !Temperature (defined on voronoi centers)
      type(scalar_field), intent(in):: htheta  !General
  
negative.code:     !Use routines from the interpolation pack
      use interpack, only: &
      plot_scalarfield, &
      plot_cart_vectorfield
  
      !Use routines from the SWM
      use swm_data      !Everything
  