Takes a grid patch with mesh widths hx and hy, of dimensions nrow by ncol, and uses the boundary conditions to set the solution values on any piece of the patch which extends outside the physical domain.
User-specified boundary conditions must be included in this routine.