From fe449c60c28ebe631fae876eae85a488091d5921 Mon Sep 17 00:00:00 2001 From: Keith Lindsay Date: Tue, 19 May 2026 15:42:20 -0600 Subject: [PATCH] Correct interior tendency alt_co2 pH bounds --- src/marbl_interior_tendency_mod.F90 | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/src/marbl_interior_tendency_mod.F90 b/src/marbl_interior_tendency_mod.F90 index 698c80af..98058227 100644 --- a/src/marbl_interior_tendency_mod.F90 +++ b/src/marbl_interior_tendency_mod.F90 @@ -1122,16 +1122,17 @@ subroutine compute_carbonate_chemistry(domain, temperature, press_bar, & call marbl_status_log%log_warning_trace('marbl_co2calc_interior() with dic', subname) end if + ph_prev_col = ph + do k=1,dkm - ph_prev_col(k) = pH(k) if (ph_prev_alt_co2_col(k) /= c0) then if (domain%delta_z(k) > 1._r8 * unit_system%cm2len) then - ph_lower_bound(k) = ph_prev_col(k) - del_ph - ph_upper_bound(k) = ph_prev_col(k) + del_ph + ph_lower_bound(k) = ph_prev_alt_co2_col(k) - del_ph + ph_upper_bound(k) = ph_prev_alt_co2_col(k) + del_ph else ! Double initial bracket width for very thin layers - ph_lower_bound(k) = ph_prev_col(k) - 2._r8 * del_ph - ph_upper_bound(k) = ph_prev_col(k) + 2._r8 * del_ph + ph_lower_bound(k) = ph_prev_alt_co2_col(k) - 2._r8 * del_ph + ph_upper_bound(k) = ph_prev_alt_co2_col(k) + 2._r8 * del_ph end if else ph_lower_bound(k) = phlo_3d_init