From 4abf0d5ac4832bc83f00b433dcf5a7edc36e432e Mon Sep 17 00:00:00 2001 From: jlinmod <134327817+jlinmod@users.noreply.github.com> Date: Thu, 18 Jul 2024 17:48:44 -0400 Subject: [PATCH 1/2] Update marbl_ciso_surface_flux_mod.F90 comment and eps_aq_g equation fix --- src/marbl_ciso_surface_flux_mod.F90 | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/src/marbl_ciso_surface_flux_mod.F90 b/src/marbl_ciso_surface_flux_mod.F90 index 862de6c3..39bf3bfc 100644 --- a/src/marbl_ciso_surface_flux_mod.F90 +++ b/src/marbl_ciso_surface_flux_mod.F90 @@ -77,16 +77,16 @@ subroutine marbl_ciso_surface_flux_compute(& real (r8), dimension(num_elements) :: & eps_aq_g_surf, & ! equilibrium fractionation (CO2_gaseous <-> CO2_aq) - alpha_aq_g_surf, & ! alpha_xxx_g_surf => eps = ( alpa -1 ) * 1000 + alpha_aq_g_surf, & ! alpha_xxx_g_surf => eps = ( alpha -1 ) * 1000 eps_dic_g_surf, & ! equilibrium fractionation between total DIC and gaseous CO2 - alpha_dic_g_surf, & ! alpha_xxx_g_surf => eps = ( alpa -1 ) * 1000 + alpha_dic_g_surf, & ! alpha_xxx_g_surf => eps = ( alpha -1 ) * 1000 frac_co3, & ! carbonate fraction fCO3 = [CO3--]/DIC alpha_aq_g_surf_14c, & ! for 14C, with fractionation being twice as large for 14C than for 13C alpha_dic_g_surf_14c ! for 14C, with fractionation being twice as large for 14C than for 13C ! local parameters for 13C, Zhang et al, 1995, Geochim. et Cosmochim. Acta, 59 (1), 107-114 real(r8) :: & - alpha_k, & ! eps = ( alpa -1 ) * 1000 + alpha_k, & ! eps = ( alpha -1 ) * 1000 alpha_k_14c ! for 14C, with fractionation being twice as large for 14C than for 13C ! kinetic fraction during gas transfer (per mil) (air-sea CO2 exchange) @@ -150,9 +150,9 @@ subroutine marbl_ciso_surface_flux_compute(& !----------------------------------------------------------------------- ! individal discrimination factor of each species with respect to - ! gaseous CO2, temperature dependent, based on Zhang et al. 95 + ! gaseous CO2, temperature dependent, based on Zhang et al. 1995 !----------------------------------------------------------------------- - eps_aq_g_surf(:) = 0.0049_r8 * sst(:) - 1.31_r8 + eps_aq_g_surf(:) = -0.0049_r8 * sst(:) - 1.31_r8 !----------------------------------------------------------------------- ! compute the equilibrium discrimination factor between dic and From 3968707effed812d667c7686e100bdf1a0de3bf6 Mon Sep 17 00:00:00 2001 From: jlinmod <134327817+jlinmod@users.noreply.github.com> Date: Thu, 19 Sep 2024 08:54:57 -0400 Subject: [PATCH 2/2] Update marbl_ciso_surface_flux_mod.F90 9-19-24 --- src/marbl_ciso_surface_flux_mod.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/marbl_ciso_surface_flux_mod.F90 b/src/marbl_ciso_surface_flux_mod.F90 index 39bf3bfc..32583d89 100644 --- a/src/marbl_ciso_surface_flux_mod.F90 +++ b/src/marbl_ciso_surface_flux_mod.F90 @@ -152,7 +152,7 @@ subroutine marbl_ciso_surface_flux_compute(& ! individal discrimination factor of each species with respect to ! gaseous CO2, temperature dependent, based on Zhang et al. 1995 !----------------------------------------------------------------------- - eps_aq_g_surf(:) = -0.0049_r8 * sst(:) - 1.31_r8 + eps_aq_g_surf(:) = 0.0049_r8 * sst(:) - 1.31_r8 !----------------------------------------------------------------------- ! compute the equilibrium discrimination factor between dic and