lambda a, _: p0_14C_with_decay_disc(a) * 31 * 12 / ALPHA_14C,
(0, 45000000),
np.array([0]))
soil_14C_by_integration_disc = soil_14C_by_integration_disc.y.reshape(
-1)[-1] * ALPHA_14C
soil_14C_by_integration_cont, soil_14C_by_integration_disc
# +
# Delta 14C in Soil in t0 = 1920 AD
print("Soil Delta 14C at t0 = 1920")
print("%2.2f ‰" % (
(soil_14C_by_integration_disc / soil_by_integration_disc / ALPHA_14C - 1) *
1000))
# -
# ## Or use CARDAMOMlib
# +
start_values_14C_disc = CARDAMOMlib.compute_start_values_14C(
dmr, nr_time_steps)
start_values_Delta_14C_disc = (
start_values_14C_disc / start_values_disc.values / ALPHA_14C - 1) * 1000
with np.printoptions(precision=2, suppress=True):
print("Initial Delta14C (t0 = 1920)\n")
print(pwc_mr.model.state_variables)
print(start_values_Delta_14C_disc, "‰")
# -
np.array([0]))
soil_14C_by_integration = soil_14C_by_integration.y.reshape(-1)[-1] * ALPHA_14C
soil_14C_by_integration
# +
# Delta 14C in Soil in t0 = 1920 AD
print("Soil Delta 14C at t0 = 1920")
print("%2.2f ‰" %
((soil_14C_by_integration / soil_by_integration / ALPHA_14C - 1) * 1000))
# -
# ## Or use CARDAMOMlib
# +
start_values_14C = CARDAMOMlib.compute_start_values_14C(pwc_mr, nr_time_steps)
start_values_Delta_14C = (start_values_14C / start_values.values / ALPHA_14C -
1) * 1000
system_Delta_14C = (
start_values_14C.sum() / start_values.values.sum() / ALPHA_14C - 1) * 1000
with np.printoptions(precision=2, suppress=True):
print("Initial Delta14C (t0 = 1920)\n")
print(pwc_mr.model.state_variables)
print(start_values_Delta_14C, "‰")
print("\nInitial Delta14C total system (t0 = 1920)\n")
print("%2.2f ‰" % system_Delta_14C)
# -
ax.axvline(mean_age_vector_disc[nr],
ymax=pool_age_density_disc(mean_age_vector_disc[nr]) /
ax.get_ylim()[-1],
ls="--",
c="red")
ax.legend()
plt.suptitle("(Fake) equilibrium age distributions")
plt.tight_layout()
# -
# ## Initial Delta 14C from (fake) initial age distributions
# +
start_values_14C_m = CARDAMOMlib.compute_start_values_14C(
pwc_mr_m, nr_time_steps_m)
start_values_Delta_14C_m = (
start_values_14C_m / start_values_m.values / ALPHA_14C - 1) * 1000
system_Delta_14C = (start_values_14C_m.sum() / start_values_m.values.sum() /
ALPHA_14C - 1) * 1000
with np.printoptions(precision=2, suppress=True):
print("Monthly initial Delta14C (t0 = 1920)\n")
print(pwc_mr_m.model.state_variables)
print(start_values_Delta_14C_m, "‰")
print("System (1920): %2.2f ‰" % system_Delta_14C)
print("\n")
start_values_14C_y0 = CARDAMOMlib.compute_start_values_14C(
pwc_mr_y0, nr_time_steps_y)
start_values_Delta_14C_y0 = (
start_values_14C_y0 / start_values_y0.values / ALPHA_14C - 1) * 1000