anchor.diff:  
           print*, '    max stage: ', maxval(domain%U(:,:,STG))
           print*, '    min stage: ', minval(domain%U(:,:,STG))
  +        print*, '    max elev: ', maxval(domain%U(:,:,ELV))
  +        print*, '    min elev: ', minval(domain%U(:,:,ELV))
   
           ! Ensure stage >= elevation
           domain%U(:,:,STG) = max(domain%U(:,:,STG), domain%U(:,:,ELV))
  
positive.code:     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)
      print*, 'skip_header_hazard_points_file: ', skip_header_hazard_points_file
  
negative.code:     character(charlen) :: timestepping_method, input_parameter_file, namelist_output_filename
  
      real(dp) :: approximate_writeout_frequency, final_time, manning_n, cliffs_minimum_allowed_depth
      real(dp):: timestep, global_ur(2), global_ll(2), global_lw(2), cfl, dx(2), linear_friction_coeff
      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
  