anchor.diff:      print*, 'output_basedir: ', TRIM(output_basedir)
       print*, 'final_time: ', final_time
       print*, 'timestepping_method: ', TRIM(timestepping_method)
  +    print*, 'linear_solver_is_truely_linear: ', linear_solver_is_truely_linear
       print*, 'manning n: ', manning_n
       print*, 'cfl: ', cfl
       print*, 'hazard_points_file: ', TRIM(hazard_points_file)
  
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:     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 ####'
  
  