\phi c_t \partial_t (P + \tau_w \partial_t P ) = \nabla (  \alpha \nabla (P + \tau_p \partial_t P ))