fix balance

elisacw · Nov 1, 2024 · 86b46e5 · 86b46e5
1 parent b982c0c
commit 86b46e5
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 8 deletions.
diff --git 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
@@ -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