anchor.diff:      domain%linear_solver_is_truely_linear = linear_solver_is_truely_linear
       domain%nc_grid_output%spatial_stride = grid_output_spatial_stride
   
  +    !domain%friction_type = 'chezy'
  +    domain%linear_friction_coeff = linear_friction_coeff
  +
       !! Optionally suppress limiting with rk2
       !!domain%theta = 4.0_dp
       
  
positive.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)
      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
  
negative.code:     namelist /MODELCONFIG/ &
          input_elevation_raster, input_stage_raster, global_ll, &
          global_ur, global_nx, approximate_writeout_frequency, output_basedir, &
          final_time, timestepping_method, manning_n, linear_friction_coeff, &
          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, cliffs_minimum_allowed_depth
  
      ! Predefine some variables that might not be in the input file
      manning_n = 0.0_dp
      linear_friction_coeff = 0.0_dp
      linear_solver_is_truely_linear = .true.
      grid_output_spatial_stride = 1
      cliffs_minimum_allowed_depth = 1.0e-03_dp
  