XH_source Ca*exp(-Ea/R_const/T)*log(max(tcs.p_wH2/p_eq(XH,T),1))*(XHmax-XH) "Absorption rate [mol H / mol M / s]"