added excretion to diagnostic of flux from benthic fauna to POC

pmlmodelling · Jun 26, 2023 · 6e824d9 · 6e824d9
1 parent 6e68f07
commit 6e824d9
Showing 1 changed file with 4 additions and 3 deletions.
diff --git 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