An algorithm for imposing additional conditions on the solution of a boundary value elliptic problem at specified grid points in the interior of the computational domain. It effectively determines a modification to the right-hand side of the governing elliptic equation which will precisely satisfy the additional interior boundary conditions. Pragmatically it allows the inclusion of island and irregular coastal boundaries while retaining the use of fast and accurate elliptic solving routines at a modest additional computational expense.