anchor.diff:      md%domains(1)%timestepping_method = default_nonlinear_timestepping_method !'cliffs' !'rk2'
       md%domains(1)%cliffs_minimum_allowed_depth = 0.01_dp
       !md%domains(1)%theta = 4.0_dp
  +    !md%domains(1)%use_eddy_viscosity = .true.
   
       print*, 1, ' lw: ', md%domains(1)%lw, ' ll: ', md%domains(1)%lower_left, ' dx: ', md%domains(1)%dx, &
           ' nx: ', md%domains(1)%nx
  
positive.code:         real(dp) :: gauge_xy(3,8), leftmost_x(3)
          type(xyz_lines_type) :: polygons
          logical :: is_inside
          real(dp), parameter :: wall = 1.0_dp, initial_reservoir_stage = 0.4_dp, initial_downstream_stage = 0.02_dp
          
  
          ! These define the dam-wall and the building
  
negative.code: 
          ! Write gauges every time-step, but print and write grids less often
          call program_timer%timer_start('IO')
          call md%write_outputs_and_print_statistics(approximate_writeout_frequency=approximate_writeout_frequency,&
              write_grids_less_often=write_grids_and_print_every_nth_step, &
              print_less_often = write_grids_and_print_every_nth_step,&
              write_gauges_less_often= 1_ip)
  