anchor.diff:          end select
           md%domains(j)%boundary_function => boundary_function
       end do
  +
  +    ! Setup hazard points
  +    call md%set_point_gauges_from_csv("point_gauges.csv", skip_header=1_ip)
      
       ! Print the gravity-wave CFL limit, to guide timestepping
       do j = 1, size(md%domains)
  
positive.code:         boundary_information%boundary_data(nr - skip + extra,2) = &
              boundary_information%boundary_data(nr - skip + extra - 1, 2)
  
          call boundary_information%gauge4_ts_function%initialise(&
                  boundary_information%boundary_data(:,1), &
                  boundary_information%boundary_data(:,2))
          boundary_information%t0 = boundary_information%boundary_data(1,1)
  
negative.code:     type :: boundary_information_type
          character(charlen):: bc_file
          real(dp), allocatable:: boundary_data(:,:)
          type(linear_interpolator_type):: gauge_ts_function
          real(dp):: boundary_elev
          real(dp):: t0 = 0.0_dp
      end type
  