diff --git a/src/soilbiogeochem/SoilBiogeochemDecompCascadeMIMICSplusMod.F90 b/src/soilbiogeochem/SoilBiogeochemDecompCascadeMIMICSplusMod.F90
index df303399cc..36c03ff9c8 100644
--- a/src/soilbiogeochem/SoilBiogeochemDecompCascadeMIMICSplusMod.F90
+++ b/src/soilbiogeochem/SoilBiogeochemDecompCascadeMIMICSplusMod.F90
@@ -72,6 +72,7 @@ module SoilBiogeochemDecompCascadeMIMICSplusMod
   public :: readParams                          ! Read in parameters from params file
   public :: init_decompcascade_mimicsplus       ! Initialization
   public :: decomp_rates_mimicsplus             ! Figure out decomposition rates
+  public :: decomp_rates_after_FUN
   public :: calc_myc_mining_rates
   public :: calc_myc_mortality
   public :: calc_myc_roi
@@ -1764,10 +1765,69 @@ subroutine decomp_rates_mimicsplus(bounds, num_bgc_soilc, filter_bgc_soilc, &
 
 end subroutine decomp_rates_mimicsplus
 
-! add a new subroutine decomp_rates_after_FUN, will be called in FUN 
+
+
+subroutine decomp_rates_after_FUN (bounds, num_bgc_soilc, filter_bgc_soilc, &
+   num_soilp, filter_soilp, &
+   soilstate_inst, cnveg_carbonflux_inst, &
+   soilbiogeochem_carbonflux_inst, soilbiogeochem_nitrogenflux_inst, &
+   soilbiogeochem_carbonstate_inst, soilbiogeochem_nitrogenstate_inst,& 
+   soilbiogeochem_state_inst &
+   )
+! 
+! DESCRIPTION
+!! add a new subroutine decomp_rates_after_FUN, will be called in FUN 
 ! col variables, do we have corresponding patch variables
 ! patch variables have to be p2c
 
+! USES:
+   use SoilBiogeochemNitrogenFluxType  , only : soilbiogeochem_nitrogenflux_type
+   use SoilBiogeochemNitrogenStateType  , only : soilbiogeochem_nitrogenstate_type
+   use subgridAveMod                    , only : p2c
+
+! ARGUMENTS:
+   type(bounds_type)                    , intent(in)    :: bounds          
+   integer                              , intent(in)    :: num_soilp       ! number of soil patches in filter
+   integer                              , intent(in)    :: filter_soilp(:) ! filter for soil patches
+   integer                              , intent(in)    :: num_bgc_soilc       ! number of soil columns in filter
+   integer                              , intent(in)    :: filter_bgc_soilc(:) ! filter for soil columns
+   type(soilstate_type)                 , intent(in)    :: soilstate_inst
+   type(cnveg_carbonflux_type)          , intent(in)    :: cnveg_carbonflux_inst
+   type(soilbiogeochem_carbonflux_type) , intent(inout) :: soilbiogeochem_carbonflux_inst
+   type(soilbiogeochem_carbonstate_type), intent(in)    :: soilbiogeochem_carbonstate_inst
+   type(soilbiogeochem_carbonflux_type) , intent(inout) :: soilbiogeochem_nitrogenflux_inst
+   type(soilbiogeochem_carbonstate_type), intent(in)    :: soilbiogeochem_nitrogenstate_inst
+   type(soilbiogeochem_state_type)      , intent(in)    :: soilbiogeochem_state_inst
+
+
+   !LOCAL VARIABLES:
+   integer   :: fn                    ! number of values in pft filter
+   integer   :: fp                    ! lake filter pft index 
+   integer   :: fc                    ! lake filter column index
+   integer   :: p, c                  ! pft index
+   integer   :: g, l                  ! indices
+   integer   :: j, i, k               ! soil/snow level index
+
+   associate ( &
+      ! Carbon
+      decomp_cpools_vr => soilbiogeochem_carbonstate_inst%decomp_cpools_vr_col , &  ! Input: [real(r8) (:,:,:) ] (gC/m3)  vertically-resolved decomposing (litter, cwd, soil) C pools
+      pathfrac_decomp_cascade => soilbiogeochem_carbonflux_inst%pathfrac_decomp_cascade_col, &  ! Output: [real(r8) (:,:,:) ]  what fraction of C leaving a given pool passes through a given transition (frac)
+      rf_decomp_cascade       => soilbiogeochem_carbonflux_inst%rf_decomp_cascade_col    , & ! Input:  [real(r8)          (:,:,:) ]  respired fraction in decomposition step (frac)
+        
+      decomp_npools_vr      => soilbiogeochem_nitrogenstate_inst%decomp_npools_vr_col         , &
+
+      )
+
+      ! CALCULATE MORTALITY FLUXES
+
+
+      ! CALCULATE FLUXES CARBON FLUXES
+
+      !
+
+   end associate
+end subroutine decomp_rates_after_FUN
+
 
    
    subroutine calc_myc_roi(cpool_myc, npool_myc,cpool_somp,cpool_soma,cpool_somc, &
@@ -2241,6 +2301,7 @@ subroutine fun_fluxes_myc_update1 (cpool_myc, npool_myc, cpool_somp, cpool_soma,
 
   end subroutine fun_fluxes_myc_update1
   
+  
           
   !Moisture function, based on testbed code: https://github.com/wwieder/biogeochem_testbed/blob/957a5c634b9f2d0b4cdba0faa06b5a91216ace33/SOURCE_CODE/mimics_cycle.f90#L401-L419
   real(r8) function r_moist(h2osoi_liq,watsat, h2osoi_ice, dz) !As in testbed (and CLM) version of MIMICS