p_{wf}(t) = p_r(t) - \frac{q(t)}{J}

p_r(t) = p_i - \frac{ int_0^t q_i(\tau) d\tau - int_0^t q(\tau) d\tau  }{c_t \, V_{\phi}}