anchor.diff:          final_time, timestepping_method, manning_n, &
           cfl, hazard_points_file, skip_header_hazard_points_file, record_max_U,&
           output_grid_timeseries, adaptive_computational_extents, negative_elevation_raster, &
  -        linear_solver_is_truely_linear
  +        linear_solver_is_truely_linear, grid_output_spatial_stride
   
       ! Predefine some variables that might not be in the input file
       manning_n = 0.0_dp
       linear_solver_is_truely_linear = .true.
  +    grid_output_spatial_stride = 1
   
       ! Read the input file -- the name of this file should be the first
       ! commandline argument
  
positive.code:     domain%record_max_U = record_max_U
      domain%output_basedir = output_basedir
      domain%linear_solver_is_truely_linear = linear_solver_is_truely_linear
      
      ! Allocate domain -- must have set timestepping method BEFORE this
      call domain%allocate_quantities(global_lw, global_nx, global_ll)
  
negative.code: 
      real(dp) :: approximate_writeout_frequency, final_time, manning_n
      real(dp):: timestep, global_ur(2), global_ll(2), global_lw(2), cfl, dx(2)
      integer(ip):: global_nx(2), skip_header_hazard_points_file, file_unit_temp, grid_output_spatial_stride
      character(len=charlen):: input_elevation_raster, input_stage_raster, &
          hazard_points_file, output_basedir
      logical:: record_max_U, output_grid_timeseries, adaptive_computational_extents, &
  