Made major modifications for ONR review. Navier BC has a new subgrid model. contact angle boundary conditions are now applied using the previous time steps contact angle. This formulation is not robust at all and only works for the withdrawing plate problem. Things run alright for a refine level of 5, but gives a segmentation fault for refine levels 6 and 7. current withdrawing plate results are closer to the experiment but it is not really a steady state result as the interface oscillates around an value.