From 6e824d9bc2c72d6ff4fa0bf7efd34f3e78d19ea5 Mon Sep 17 00:00:00 2001 From: Rebecca Millington Date: Mon, 26 Jun 2023 14:15:27 +0100 Subject: [PATCH] added excretion to diagnostic of flux from benthic fauna to POC --- src/benthic_fauna.F90 | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/src/benthic_fauna.F90 b/src/benthic_fauna.F90 index 942430a..d30a2df 100644 --- a/src/benthic_fauna.F90 +++ b/src/benthic_fauna.F90 @@ -420,9 +420,10 @@ subroutine do_bottom(self,_ARGUMENTS_DO_BOTTOM_) _SET_HORIZONTAL_DIAGNOSTIC_(self%id_fYKIn,excess_n) _SET_HORIZONTAL_DIAGNOSTIC_(self%id_fYK1p,excess_p) - _SET_HORIZONTAL_DIAGNOSTIC_(self%id_fYQPc,mortflux*cP+excess_c) - _SET_HORIZONTAL_DIAGNOSTIC_(self%id_fYQPn,mortflux*cP*self%qnc) - _SET_HORIZONTAL_DIAGNOSTIC_(self%id_fYQPp,mortflux*cP*self%qpc) + ! fluxes from fauna to POC + _SET_HORIZONTAL_DIAGNOSTIC_(self%id_fYQPc,mortflux*cP+excess_c+fBTYc-nfBTYc) + _SET_HORIZONTAL_DIAGNOSTIC_(self%id_fYQPn,mortflux*cP*self%qnc+sum(grossfluxn)-sum(netfluxn)) + _SET_HORIZONTAL_DIAGNOSTIC_(self%id_fYQPp,mortflux*cP*self%qpc+sum(grossfluxp) - sum(netfluxp)) if (.not.legacy_ersem_compatibility) then ! Alkalinity contributions: +1 for NH4, -1 for PO4