resTheta_r = (Theta_r*thk_r-bblTheta_r*bblEta_r)/resThk_r
bbl_tend = dVol * (bblTheta_d - resTheta_r) / deltaT
_RL Theta_r, bblTheta_d, bblTheta_r, resTheta_r