Cleaned up boundary.c. Added new boundary condition identified as BcJoseph and uses input similar to BcNavier. Added new functions to boundary.c including is_interfacial, find_tcp, and boundary_tangential_gradient. NOTE: need to fix face_joseph so that the fluid velocity at the wall is accurately reflected using the discretization of Joseph's BC