From 86b46e5ef673040e9ae13d62932f0163c85aeb39 Mon Sep 17 00:00:00 2001 From: mvdebolskiy Date: Fri, 1 Nov 2024 18:22:48 +0100 Subject: [PATCH] fix balance --- src/biogeochem/CNFUNMIMICSplusMod.F90 | 7 +++---- src/soilbiogeochem/SoilBiogeochemDecompMod.F90 | 6 ++---- 2 files changed, 5 insertions(+), 8 deletions(-) diff --git a/src/biogeochem/CNFUNMIMICSplusMod.F90 b/src/biogeochem/CNFUNMIMICSplusMod.F90 index 8154636096..47a87be44f 100644 --- a/src/biogeochem/CNFUNMIMICSplusMod.F90 +++ b/src/biogeochem/CNFUNMIMICSplusMod.F90 @@ -1100,12 +1100,11 @@ subroutine CNFUNMIMICSplus (bounds, num_soilc, filter_soilc, num_soilp ,filter_s Nuptake(p) = Necm(p) + Nam(p) + Nnonmyc_no3(p) + Nnonmyc_nh4(p) + Nfix(p) + Nretrans(p) !---------------------------Extra Respiration Fluxes--------------------! - soilc_change(p) = npp_Necm(p) + npp_Nam(p) + npp_Nnonmyc_no3(p) + npp_Nnonmyc_nh4(p) + npp_Nfix(p) + npp_Nretrans(p) - !soilc_change(p) = soilc_change(p) + soilc_change(p) = npp_Nnonmyc_no3(p) + npp_Nnonmyc_nh4(p) + npp_Nfix(p) + npp_Nretrans(p) ! mycorrhiza fluxes are used separately + soilc_change(p) = soilc_change(p) npp_burnedoff(p) = burned_off_carbon/dt npp_Nuptake(p) = soilc_change(p) - npp_Nuptake(p) = npp_Nuptake(p) * (1.0_r8 - grperc(ivt(p))) - soilc_change(p) = soilc_change(p) * grperc(ivt(p)) + npp_Nuptake(p) = npp_Necm(p) + npp_Nam(p) + npp_Nnonmyc_no3(p) + npp_Nnonmyc_nh4(p) + npp_Nfix(p) + npp_Nretrans(p) ! how much carbon goes to growth of tissues? !npp_growth(p) = (Nuptake(p)- free_retransn_to_npool(p))*plantCN(p)+(excess_carbon_acc/dt) !does not include gresp, since this is calculated from growth diff --git a/src/soilbiogeochem/SoilBiogeochemDecompMod.F90 b/src/soilbiogeochem/SoilBiogeochemDecompMod.F90 index dd1dc15766..a0f3200ff2 100644 --- a/src/soilbiogeochem/SoilBiogeochemDecompMod.F90 +++ b/src/soilbiogeochem/SoilBiogeochemDecompMod.F90 @@ -326,11 +326,9 @@ subroutine SoilBiogeochemDecomp (bounds, num_bgc_soilc, filter_bgc_soilc, end do ! transitions ! fluxes that are not part of the cascade. decomp_cpools_vr(c,j,i_ecm_myc) = decomp_cpools_vr(c,j,i_ecm_myc) + c_ecm_growth_vr(c,j) * dt - decomp_npools_vr(c,j,i_ecm_myc) = decomp_npools_vr(c,j,i_ecm_myc) + (n_ecm_growth_vr(c,j) + & - sminno3_to_ecm_vr(c,j) + sminnh4_to_ecm_vr(c,j)) * dt + decomp_npools_vr(c,j,i_ecm_myc) = decomp_npools_vr(c,j,i_ecm_myc) + n_ecm_growth_vr(c,j) * dt decomp_cpools_vr(c,j,i_am_myc) = decomp_cpools_vr(c,j,i_am_myc) + c_am_growth_vr(c,j) * dt - decomp_npools_vr(c,j,i_am_myc) = decomp_npools_vr(c,j,i_am_myc) + (n_ecm_growth_vr(c,j) + & - sminno3_to_am_vr(c,j) + sminnh4_to_am_vr(c,j)) * dt + decomp_npools_vr(c,j,i_am_myc) = decomp_npools_vr(c,j,i_am_myc) + n_ecm_growth_vr(c,j) * dt decomp_cpools_vr(c,j,i_avl_som) = decomp_cpools_vr(c,j,i_avl_som) + c_ecm_enz_vr(c,j) * dt end do ! layer enddo !column