anchor.diff:          call setprob()
           call set_gauges(rest, nvar)
           call set_regions()
  -        call set_slices()
  +!        call set_slices()
   
           cflmax = 0.d0   ! otherwise use previously heckpointed val
   
  
positive.code: 
      use regions_module, only: set_regions
      use gauges_module, only: set_gauges, num_gauges
      use slices_module, only: set_slices, print_slices, num_slices
  
      implicit none
  
  
negative.code: 
          call set_gauges(rest, nvar)
          call set_regions()
  !        call set_slices()
      else
  
          open(outunit, file=outfile, status='unknown', form='formatted')
  