 Depth-dependent inclined porous-current equation

For porosity $P{z}^{\alpha -1}$ and permeability of a porous medium $K{z}^{\beta -1}$ above an inclined impermeable plane, stored liquid volume per area is $P{h}^{\alpha }/\alpha$. The depth-integrated Darcy flux is $\rho gK{h}^{\beta }(\sin \theta -\cos \theta h_x,-\cos \theta h_y)/(\beta \mu )$. Conservation of mass yields the displayed equation with $X=x\tan \theta$, $Y=y\tan \theta$ and $T=t\alpha \rho gK\sin \theta \tan \theta /(\beta \mu P)$.