Including grazing
(Wroblewski et al., 1975; Denman and Platt, 1975)
     Pt = kPxx + mP - R[1-exp(-lP)]
     - RlP
     Lc=pÖ[k/(m-Rl)]
But this is only for small times…