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:     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:     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 ####'
  
  