mescla.isotopes.two_component_isotope¶
- mescla.isotopes.two_component_isotope(delta_sample, delta_1, carrier_1, delta_2, carrier_2)[source]¶
Fraction of end-member 1 from a solute isotope ratio, two components.
\[f_1 = \frac{C_2(\delta_2 - \delta_s)} {C_1(\delta_s - \delta_1) + C_2(\delta_2 - \delta_s)}\]The textbook formula. Note that it reduces to the familiar \((\delta_s - \delta_2)/(\delta_1 - \delta_2)\) only when \(C_1 = C_2\) – which is why applying the simple form to a solute isotope is wrong whenever the end-members differ in the carrier element.
Returned unclipped, like
two_component(): a value outside [0, 1] is information, not an error.