* move the contact angle functions into twophaseCoeff.cpp