Because of Matlab derailing for values of k3 and k4 used in the paper (0.1) but working for lower values, 0.01, we want to know what values are physically relevant. Sulpis et al. (2022) propose different values of k and n depending on calcite saturation state.
Do these values (and intervals they consider) apply to us? For this we need to know what Ω_calcite values we see in the system so we need to track its value two_factors_low_lim - 1) ** self.n1 (Eq. 45).
This is preferable to calculating them each time from raw output of Ca and CO3 concentrations because we want to explore the effect of changing k3, k4 and n1 quickly.
Because of Matlab derailing for values of
k3andk4used in the paper (0.1) but working for lower values, 0.01, we want to know what values are physically relevant. Sulpis et al. (2022) propose different values ofkandndepending on calcite saturation state. Do these values (and intervals they consider) apply to us? For this we need to know what Ω_calcite values we see in the system so we need to track its valuetwo_factors_low_lim - 1) ** self.n1(Eq. 45). This is preferable to calculating them each time from raw output of Ca and CO3 concentrations because we want to explore the effect of changingk3,k4andn1quickly.