MIMICS into vertically resolved CTSM

src/biogeochem/CNFUNMod.F90

@@ -505,7 +505,6 @@ subroutine CNFUN(bounds,num_soilc, filter_soilc,num_soilp&
   !---------------------------------
   associate(ivt                    => patch%itype                                          , & ! Input:   [integer  (:) ]  p
          leafcn                 => pftcon%leafcn                                        , & ! Input:   leaf C:N (gC/gN)
-         lflitcn                => pftcon%lflitcn                                       , & ! Input:   leaf litter C:N (gC/gN)
          season_decid           => pftcon%season_decid                                  , & ! Input:   binary flag for seasonal
          ! -deciduous leaf habit (0 or 1)
          stress_decid           => pftcon%stress_decid                                  , & ! Input:   binary flag for stress

src/biogeochem/CNGRespMod.F90

@@ -60,7 +60,6 @@ subroutine CNGResp(num_soilp, filter_soilp, cnveg_carbonflux_inst, canopystate_i
          grperc                        =>    pftcon%grperc                                             , & ! Input:  growth respiration parameter
          grpnow                        =>    pftcon%grpnow                                             , & ! Input:  growth respiration parameter
          leafcn                        =>    pftcon%leafcn                                             , & ! Input:  leaf C:N (gC/gN)                         
-         frootcn                       =>    pftcon%frootcn                                            , & ! Input:  fine root C:N (gC/gN)                    
          livewdcn                      =>    pftcon%livewdcn                                           , & ! Input:  live wood (phloem and ray parenchyma) C:N (gC/gN)  
 
          laisun                        =>    canopystate_inst%laisun_patch                             , & ! Input:  [real(r8) (:)]  sunlit projected leaf area index      

src/biogeochem/CNGapMortalityMod.F90

@@ -134,7 +134,6 @@ subroutine CNGapMortality (bounds, num_soilc, filter_soilc, num_soilp, filter_so
          heatstress =>    dgvs_inst%heatstress_patch , & ! Input:  [real(r8) (:) ]    
 
          leafcn  	=>    pftcon%leafcn               , & ! Input:  [real(r8) (:)]  leaf C:N (gC/gN)                        
-         frootcn    =>    pftcon%frootcn              , & ! Input:  [real(r8) (:)]  fine root C:N (gC/gN)                  
          livewdcn   =>    pftcon%livewdcn             , & ! Input:  [real(r8) (:)]  live wood (phloem and ray parenchyma) C:N (gC/gN) 
          laisun     =>    canopystate_inst%laisun_patch  , & ! Input:  [real(r8) (:)   ]  sunlit projected leaf area index      
          laisha     =>    canopystate_inst%laisha_patch  , & ! Input:  [real(r8) (:)   ]  shaded projected leaf area index   

src/biogeochem/CNNDynamicsMod.F90

@@ -218,7 +218,7 @@ subroutine CNNFixation(num_soilc, filter_soilc, &
     !-----------------------------------------------------------------------
 
     associate(                                                                & 
-         cannsum_npp    => cnveg_carbonflux_inst%annsum_npp_col ,             & ! Input:  [real(r8) (:)]  nitrogen deposition rate (gN/m2/s)                
+         cannsum_npp    => cnveg_carbonflux_inst%annsum_npp_col ,             & ! Input: [real(r8) (:)]  (gC/m2/yr) annual sum of NPP, averaged from patch-level
          col_lag_npp    => cnveg_carbonflux_inst%lag_npp_col    ,             & ! Input: [real(r8) (:)]  (gC/m2/s) lagged net primary production           
 
          nfix_to_sminn  => soilbiogeochem_nitrogenflux_inst%nfix_to_sminn_col & ! Output: [real(r8) (:)]  symbiotic/asymbiotic N fixation to soil mineral N (gN/m2/s)

src/biogeochem/CNVegCarbonFluxType.F90

@@ -11,13 +11,13 @@ module CNVegCarbonFluxType
   use decompMod                          , only : bounds_type
   use SoilBiogeochemDecompCascadeConType , only : decomp_cascade_con
   use clm_varpar                         , only : ndecomp_cascade_transitions, ndecomp_pools
-  use clm_varpar                         , only : nlevdecomp_full, nlevdecomp, i_litr_min, i_litr_max
+  use clm_varpar                         , only : nlevdecomp_full, nlevdecomp, i_litr_min, i_litr_max, i_cwdl2
   use clm_varcon                         , only : spval, dzsoi_decomp
   use clm_varctl                         , only : use_cndv, use_c13, use_nitrif_denitrif, use_crop
   use clm_varctl                         , only : use_grainproduct
   use clm_varctl                         , only : iulog
   use landunit_varcon                    , only : istsoil, istcrop, istdlak 
-  use pftconMod                          , only : npcropmin
+  use pftconMod                          , only : npcropmin, pftcon
   use LandunitType                       , only : lun                
   use ColumnType                         , only : col                
   use PatchType                          , only : patch                
@@ -4004,7 +4004,7 @@ subroutine Summary_carbonflux(this, &
     use clm_varctl                         , only: nfix_timeconst, carbon_resp_opt
     use subgridAveMod                      , only: p2c, c2g
     use SoilBiogeochemDecompCascadeConType , only: decomp_cascade_con
-    use CNSharedParamsMod                  , only: use_fun
+    use CNSharedParamsMod                  , only: use_fun, CNParamsShareInst
     !
     ! !ARGUMENTS:
     class(cnveg_carbonflux_type)   :: this
@@ -4396,6 +4396,19 @@ subroutine Summary_carbonflux(this, &
             this%hrv_gresp_storage_to_litter_patch(p)      + &
             this%hrv_gresp_xfer_to_litter_patch(p)
 
+       if (use_mimics_decomp) then
+          ! Calculate ligninNratio for leaves and fine roots
+          ! TODO Declare ligninNratio vars wherever they need to be and
+          !      initialize if necessary for call p2c to work correctly
+          associate(ivt => patch%itype)  ! Input: [integer (:)] patch plant type
+            ligninNratio_leaf_patch(p) = pftcon%lf_flig(ivt(p)) * &
+                                         pftcon%lflitcn(ivt(p)) * &
+                                         this%leafc_to_litter_patch(p)
+            ligninNratio_froot_patch(p) = pftcon%fr_flig(ivt(p)) * &
+                                          pftcon%frootcn(ivt(p)) * &
+                                          this%frootc_to_litter_patch(p)
+          end associate
+       end if
     end do  ! end of patches loop
 
     !------------------------------------------------
@@ -4440,6 +4453,36 @@ subroutine Summary_carbonflux(this, &
          this%gpp_patch(bounds%begp:bounds%endp), &
          this%gpp_col(bounds%begc:bounds%endc))
 
+    if (use_mimics_decomp) then
+       call p2c(bounds, num_soilc, filter_soilc, &
+            ligninNratio_leaf_patch(bounds%begp:bounds%endp), &
+            ligninNratio_leaf_col(bounds%begc:bounds%endc))
+       call p2c(bounds, num_soilc, filter_soilc, &
+            ligninNratio_froot_patch(bounds%begp:bounds%endp), &
+            ligninNratio_froot_col(bounds%begc:bounds%endc))
+       call p2c(bounds, num_soilc, filter_soilc, &
+            this%leafc_to_litter_patch(bounds%begp:bounds%endp), &
+            leafc_to_litter_col(bounds%begc:bounds%endc))
+       call p2c(bounds, num_soilc, filter_soilc, &
+            this%frootc_to_litter_patch(bounds%begp:bounds%endp), &
+            frootc_to_litter_col(bounds%begc:bounds%endc))
+
+       ! Calculate ligninNratioAve
+       do fc = 1,num_soilc
+          c = filter_soilc(fc)
+          ! TODO Use cwd_flig in *CascadeBGC the same way for consistency?
+          ligninNratio_cwd = CNParamsShareInst%cwd_flig * &
+                             (this%cwdc_col(c) / this%cwdn_col(c)) * &
+                             this%decomp_cascade_ctransfer_col(c,i_cwdl2)
+          this%ligninNratioAvg_col(c) = &
+             (ligninNratio_leaf_col(c) + ligninNratio_froot_col(c) + &
+              ligninNratio_cwd) / &
+              max(1.0e-3_r8, leafc_to_litter_col(c) + &
+                             frootc_to_litter_col(c) + &
+                             this%decomp_cascade_ctransfer_col(c,i_cwdl2))
+       end do
+    end if
+
     ! this code is to calculate an exponentially-relaxed npp value for use in NDynamics code
 
     if ( trim(isotope) == 'bulk') then

src/main/clm_varpar.F90

@@ -9,7 +9,8 @@ module clm_varpar
   use shr_sys_mod  , only: shr_sys_abort
   use spmdMod      , only: masterproc
   use clm_varctl   , only: use_extralakelayers, use_vertsoilc
-  use clm_varctl   , only: use_century_decomp, use_c13, use_c14
+  use clm_varctl   , only: use_century_decomp, use_mimics_decomp
+  use clm_varctl   , only: use_c13, use_c14
   use clm_varctl   , only: iulog, use_crop, create_crop_landunit, irrigate
   use clm_varctl   , only: use_vichydro, rundef
   use clm_varctl   , only: soil_layerstruct_predefined
@@ -72,6 +73,8 @@ module clm_varpar
   integer, public :: i_litr_max = -9  ! max index of litter pools; overwritten in SoilBiogeochemDecompCascade*Mod
   integer, public :: i_met_lit = -9  ! index of metabolic litter pool; overwritten in SoilBiogeochemDecompCascade*Mod
   integer, public :: i_cwd      = -9  ! index of cwd pool; overwritten in SoilBiogeochemDecompCascade*Mod
+  ! TODO Not necessary but cleaner to use this copy of i_cwdl2 in *CascadeBGC?
+  integer, public :: i_cwdl2 = -9  ! index of cwd to l2 transition; overwritten in SoilBiogeochemDecompCascade*Mod
 
   integer, public :: ndecomp_pools_max
   integer, public :: ndecomp_pools
@@ -237,22 +240,28 @@ subroutine clm_varpar_init(actual_maxsoil_patches, actual_numcft)
     ! in InitAllocate (in SoilBiogeochemStateType) which is called earlier than
     ! init_decompcascade_bgc where they might have otherwise been derived on the
     ! fly. For reference, if they were determined in init_decompcascade_bgc:
-    ! ndecomp_pools would get the value of i_pas_som or i_cwd and
-    ! ndecomp_cascade_transitions would get the value of i_s3s1 or i_cwdl3
+    ! ndecomp_pools would get the value of i_oli_mic or i_cwd and
+    ! ndecomp_cascade_transitions would get the value of i_m2s3 or i_cwdl2
     ! depending on how use_fates is set.
     if ( use_fates ) then
        if (use_century_decomp) then
           ndecomp_pools = 6
           ndecomp_cascade_transitions = 8
-       else  ! TODO slevis: Currently for CN. MIMICS will get its own.
+       else if (use_mimics_decomp) then
+          ndecomp_pools = 7
+          ndecomp_cascade_transitions = 14
+       else  ! TODO slevis: Currently for CN (deprecated)
           ndecomp_pools = 7
           ndecomp_cascade_transitions = 7
        end if
     else
        if (use_century_decomp) then
           ndecomp_pools = 7
           ndecomp_cascade_transitions = 10
-       else  ! TODO slevis: Currently for CN. MIMICS will get its own.
+       else if (use_mimics_decomp) then
+          ndecomp_pools = 8
+          ndecomp_cascade_transitions = 15
+       else  ! TODO slevis: Currently for CN (deprecated)
           ndecomp_pools = 8
           ndecomp_cascade_transitions = 9
        end if

src/soilbiogeochem/SoilBiogeochemDecompCascadeConType.F90

@@ -27,6 +27,7 @@ module SoilBiogeochemDecompCascadeConType
      character(len=8)  , pointer  :: decomp_pool_name_history(:)       ! name of pool for history files
      character(len=20) , pointer  :: decomp_pool_name_long(:)          ! name of pool for netcdf long names
      character(len=8)  , pointer  :: decomp_pool_name_short(:)         ! name of pool for netcdf short names
+     logical           , pointer  :: is_microbe(:)                     ! TRUE => pool is a microbe pool
      logical           , pointer  :: is_litter(:)                      ! TRUE => pool is a litter pool
      logical           , pointer  :: is_soil(:)                        ! TRUE => pool is a soil pool
      logical           , pointer  :: is_cwd(:)                         ! TRUE => pool is a cwd pool
@@ -72,6 +73,7 @@ subroutine init_decomp_cascade_constants( use_century_decomp )
     allocate(decomp_cascade_con%decomp_pool_name_history(ibeg:ndecomp_pools))
     allocate(decomp_cascade_con%decomp_pool_name_long(ibeg:ndecomp_pools))
     allocate(decomp_cascade_con%decomp_pool_name_short(ibeg:ndecomp_pools))
+    allocate(decomp_cascade_con%is_microbe(ibeg:ndecomp_pools))
     allocate(decomp_cascade_con%is_litter(ibeg:ndecomp_pools))
     allocate(decomp_cascade_con%is_soil(ibeg:ndecomp_pools))
     allocate(decomp_cascade_con%is_cwd(ibeg:ndecomp_pools))
@@ -93,6 +95,7 @@ subroutine init_decomp_cascade_constants( use_century_decomp )
     decomp_cascade_con%decomp_pool_name_restart(ibeg:ndecomp_pools)       = ''
     decomp_cascade_con%decomp_pool_name_long(ibeg:ndecomp_pools)          = ''
     decomp_cascade_con%decomp_pool_name_short(ibeg:ndecomp_pools)         = ''
+    decomp_cascade_con%is_microbe(ibeg:ndecomp_pools)                     = .false.
     decomp_cascade_con%is_litter(ibeg:ndecomp_pools)                      = .false.
     decomp_cascade_con%is_soil(ibeg:ndecomp_pools)                        = .false.
     decomp_cascade_con%is_cwd(ibeg:ndecomp_pools)                         = .false.