Application of harmonic coordinates to 2D interface problems on regular grids