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