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