This restriction of V to places where Wstag is positive is correct and may help time stepping.  It does help diagnostic output.
forgot this: it is correct to have P=overburden when there is no water