https://github.com/fschuch/Xcompact3d/blob/e99a2a77a12d77187a63bb87836d266deed53c2e/src/ibm.f90#L95 For consistency, we should force phi to zero inside the geometry, just like the lagrangian interpolator does.
Xcompact3d/src/ibm.f90
Line 95 in e99a2a7
For consistency, we should force phi to zero inside the geometry, just like the lagrangian interpolator does.