anchor.diff:          domain%forcing_subroutine => apply_rainfall_forcing
           domain%forcing_context_cptr = c_loc(rainfall_region_indices)
           call domain%store_forcing()
  +
  +        deallocate(x,y, is_inside_poly)
   
       end subroutine
   
  
positive.code:         domain%U(domain%nx(1), :, ELV) = domain%U(domain%nx(1)-1, :, ELV)
          domain%U(:, domain%nx(2), ELV) = domain%U(:, domain%nx(2)-1, ELV) 
  
  
          !
          ! Set friction
          !
  
negative.code:         class(domain_type), target, intent(inout):: domain
  
          integer(ip):: i, j, k
          real(dp), allocatable:: x(:), y(:)
  
          ! For reading background elevation
          character(len=charlen):: input_elevation_file
          type(multi_raster_type) :: elevation_data
  
          ! Add this elevation to all points in houses/ csv polygons
          real(dp), parameter:: house_height = 3.0_dp
  