anchor.diff:          call read_csv_into_array(polygon_coords, polygon_filename)
           inside_point_counter = 0
           do j = 1, domain%nx(2)
  -
  -            call points_in_poly(polygon_coords(1,:), polygon_coords(2,:), x(:,j), y(:,j), is_inside_poly)
  +            y = domain%y(j)
  +            call points_in_poly(polygon_coords(1,:), polygon_coords(2,:), x, y, is_inside_poly)
               inside_point_counter = inside_point_counter + count(is_inside_poly)
   
               do i = 1, domain%nx(1) 
  
positive.code: 
              do j = 1, domain%nx(2)
                  ! Find points in the polygon in column j
                  call points_in_poly(polygon_coords(1,:), polygon_coords(2,:), x(:,j), y(:,j), is_inside_poly)
                  inside_point_counter = inside_point_counter + count(is_inside_poly)
  
                  do i = 1, domain%nx(1) 
  
negative.code: 
      use global_mod, only: dp, ip, force_double, charlen, wall_elevation, pi
      use domain_mod, only: domain_type, STG, UH, VH, ELV
      use read_raster_mod, only: multi_raster_type
      use which_mod, only: which
      use ragged_array_mod, only: ragged_array_2d_ip_type
      use file_io_mod, only: read_character_file, read_csv_into_array
  