anchor.diff:          real(dp) :: expected_grid(10,10)
           real(dp) :: lx(3), ly(3), lz(3)
           integer(ip) :: i
  -
  -        lower_left = [0.0, 0.0]
  -        upper_right = [10.0, 10.0]
  +        real(dp), parameter :: eps = 1.0e-03_dp
  +
  +        lower_left = [0.0_dp, 0.0_dp]
  +        upper_right = [10.0_dp, 10.0_dp]
   
           x = [1.6_dp, 2.2_dp, 3.6_dp, 8.7_dp, 9.9_dp]
           y = [0.5_dp, 1.3_dp, 2.6_dp, 3.7_dp, 4.8_dp]
  
positive.code:     !! provides routines to burn these into grids. For example, this can be used to ensure that 
      !! an elevation grid includes a continuous representation of a breakwall or riverwall.
      !!
      use global_mod, only: ip, dp, charlen
      use logging_mod
      use stop_mod
      use file_io_mod, only: read_csv_into_array
  
negative.code:             print*, 'FAIL'
          end if
  
          !print*, 'grid(1:6, 4): ', grid(1:6, 4)
          !print*, 'expected_grid(1:6, 4): ', expected_grid(1:6, 4) 
          !print*, 'grid(6, 5:10): ', grid(6, 5:10)
          !print*, 'expected_grid(6, 5:10): ', expected_grid(6, 5:10) 
  
      end subroutine
  
  end module
  