anchor.diff:      print*, 'output_grid_timeseries: ', output_grid_timeseries
       print*, 'adaptive_computational_extents: ', adaptive_computational_extents
       print*, 'negative_elevation_raster: ', negative_elevation_raster
  +    print*, 'grid_output_spatial_stride: ', grid_output_spatial_stride
       print*, ''
       print*, '#### END INPUT ####'
   
  
positive.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
      character(len=charlen):: input_elevation_raster, input_stage_raster, &
          hazard_points_file, output_basedir
      logical:: record_max_U, output_grid_timeseries, adaptive_computational_extents, &
  
negative.code:         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, 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
  