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