Merge pull request #109 from pmlmodelling/108_benthic_fauna_POC_diagn…

…ostic added excretion to diagnostic of flux from benthic fauna to POC
pmlmodelling · Jun 30, 2023 · 2e88016 · 2e88016
2 parents 004e1e2 + 6e824d9
commit 2e88016
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