Bugfix: update operator

PON%remin is PON%sflux_in * dzr_loc + ..., not PON%sflux_in + dzr_loc + ...
mnlevy1981 · Jan 9, 2025 · 01f796f · 01f796f
1 parent 2d313de
commit 01f796f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/marbl_interior_tendency_mod.F90 b/src/marbl_interior_tendency_mod.F90
@@ -3294,7 +3294,7 @@ subroutine compute_particulate_terms(k, domain, bot_flux_to_tend, p_remin_scalef
 
         if (PON%sflux_out(k) < c0) then
            PON%sflux_out(k) = c0
-           PON%remin(k) = PON%sflux_in(k) + dzr_loc + &
+           PON%remin(k) = PON%sflux_in(k) * dzr_loc + &
               (c1 - PON%gamma) * PON%prod(k)
         end if