Skip to content

Commit

Permalink
Merge branch 'other_remin_bugfix' into development
Browse files Browse the repository at this point in the history
  • Loading branch information
mnlevy1981 committed Sep 10, 2021
2 parents e2061c3 + 08f59d4 commit 51fac36
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions src/marbl_ciso_interior_tendency_mod.F90
Original file line number Diff line number Diff line change
Expand Up @@ -1255,10 +1255,10 @@ subroutine compute_particulate_terms(k, domain, bot_flux_to_tend, tracer_local,
(0.06_r8 + 0.19_r8 * 0.99_r8**(O2_loc-NO3_loc))

flux_alt = POC_ciso%to_floor*1.0e-6_r8*spd*365.0_r8 ! convert to mmol/cm^2/year
other_remin = min(bot_flux_to_tend(1:k) * &
min(0.1_r8 + flux_alt,0.5_r8) * (POC_ciso%to_floor - POC_ciso%sed_loss(k)), &
bot_flux_to_tend(1:k) * (POC_ciso%to_floor - POC_ciso%sed_loss(k)) - &
sed_denitrif(1:k) * denitrif_C_N)
other_remin(1:k) = min(bot_flux_to_tend(1:k) * &
min(0.1_r8 + flux_alt,0.5_r8) * (POC_ciso%to_floor - POC_ciso%sed_loss(k)), &
bot_flux_to_tend(1:k) * (POC_ciso%to_floor - POC_ciso%sed_loss(k)) - &
sed_denitrif(1:k) * denitrif_C_N)

! if bottom water O2 is depleted, assume all remin is denitrif + other
where (O2_loc_col(1:k) < c1)
Expand Down

0 comments on commit 51fac36

Please sign in to comment.