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