The implementation of Josephs Boundary condition seems to work for both squares and cylinders in some, but not all cases. If the solid boundary happens to lie exactly on the grid boundary, then max residual warnings occur. Max residuals are significanly higher than the tolerance (10^2 vs a tolerance of 10^-6). I think this is because it has difficulty computing the triple contact point intersection and limiting it to be within the cell. Also floating point exceptions can sometimes occur if multiple boxes are used. This does not occur in all case so I'm not exactly sure what is causing this