#include <misc.h>
#include <preproc.h>
module CNC13StateUpdate1Mod 1,1
#if (defined CN) && (defined C13)
!-----------------------------------------------------------------------
!BOP
!
! !MODULE: C13StateUpdate1Mod
!
! !DESCRIPTION:
! Module for carbon state variable update, non-mortality fluxes.
!
! !USES:
use shr_kind_mod, only: r8 => shr_kind_r8
implicit none
save
private
!
! !PUBLIC MEMBER FUNCTIONS:
public:: C13StateUpdate1,C13StateUpdate0
!
! !REVISION HISTORY:
! 4/21/2005: Created by Peter Thornton and Neil Suits - copied from
! CNCStateUpdate1 for C13 state variables.
!
!EOP
!-----------------------------------------------------------------------
contains
!-----------------------------------------------------------------------
!BOP
!
! !IROUTINE: C13StateUpdate0
!
! !INTERFACE:
subroutine C13StateUpdate0(num_soilp, filter_soilp) 1,3
!
! !DESCRIPTION:
! On the radiation time step, update cpool carbon state
!
! !USES:
use clmtype
use clm_time_manager, only: get_step_size
!
! !ARGUMENTS:
implicit none
integer, intent(in) :: num_soilp ! number of soil pfts in filter
integer, intent(in) :: filter_soilp(:) ! filter for soil pfts
!
! !CALLED FROM:
! subroutine CNEcosystemDyn
!
! !REVISION HISTORY:
! 7/1/05: Created by Peter Thornton
!
! !LOCAL VARIABLES:
! local pointers to implicit in arrays
real(r8), pointer :: psnshade_to_cpool(:)
real(r8), pointer :: psnsun_to_cpool(:)
!
! local pointers to implicit in/out arrays
real(r8), pointer :: cpool(:) ! (gC/m2) temporary photosynthate C pool
! !OTHER LOCAL VARIABLES:
integer :: p ! indices
integer :: fp ! lake filter indices
real(r8):: dt ! radiation time step (seconds)
!
!EOP
!-----------------------------------------------------------------------
! assign local pointers at the pft level
cpool => clm3%g%l%c%p%pc13s%cpool
psnshade_to_cpool => clm3%g%l%c%p%pc13f%psnshade_to_cpool
psnsun_to_cpool => clm3%g%l%c%p%pc13f%psnsun_to_cpool
! set time steps
dt = real( get_step_size(), r8 )
! pft loop
do fp = 1,num_soilp
p = filter_soilp(fp)
! gross photosynthesis fluxes
cpool(p) = cpool(p) + psnsun_to_cpool(p)*dt
cpool(p) = cpool(p) + psnshade_to_cpool(p)*dt
end do
end subroutine C13StateUpdate0
!-----------------------------------------------------------------------
!-----------------------------------------------------------------------
!BOP
!
! !IROUTINE: C13StateUpdate1
!
! !INTERFACE:
subroutine C13StateUpdate1(num_soilc, filter_soilc, num_soilp, filter_soilp) 1,3
!
! !DESCRIPTION:
! On the radiation time step, update all the prognostic carbon state
! variables (except for gap-phase mortality and fire fluxes)
!
! !USES:
use clmtype
use clm_time_manager, only: get_step_size
!
! !ARGUMENTS:
implicit none
integer, intent(in) :: num_soilc ! number of soil columns filter
integer, intent(in) :: filter_soilc(:) ! filter for soil columns
integer, intent(in) :: num_soilp ! number of soil pfts in filter
integer, intent(in) :: filter_soilp(:) ! filter for soil pfts
!
! !CALLED FROM:
! subroutine CNEcosystemDyn
!
! !REVISION HISTORY:
! 8/1/03: Created by Peter Thornton
! 12/5/03, Peter Thornton: Added livewood turnover fluxes
!
! !LOCAL VARIABLES:
! local pointers to implicit in arrays
!
real(r8), pointer :: woody(:) ! binary flag for woody lifeform (1=woody, 0=not woody)
real(r8), pointer :: cwdc_to_litr2c(:)
real(r8), pointer :: cwdc_to_litr3c(:)
real(r8), pointer :: frootc_to_litr1c(:)
real(r8), pointer :: frootc_to_litr2c(:)
real(r8), pointer :: frootc_to_litr3c(:)
real(r8), pointer :: leafc_to_litr1c(:)
real(r8), pointer :: leafc_to_litr2c(:)
real(r8), pointer :: leafc_to_litr3c(:)
real(r8), pointer :: litr1_hr(:)
real(r8), pointer :: litr1c_to_soil1c(:)
real(r8), pointer :: litr2_hr(:)
real(r8), pointer :: litr2c_to_soil2c(:)
real(r8), pointer :: litr3_hr(:)
real(r8), pointer :: litr3c_to_soil3c(:)
real(r8), pointer :: soil1_hr(:)
real(r8), pointer :: soil1c_to_soil2c(:)
real(r8), pointer :: soil2_hr(:)
real(r8), pointer :: soil2c_to_soil3c(:)
real(r8), pointer :: soil3_hr(:)
real(r8), pointer :: soil3c_to_soil4c(:)
real(r8), pointer :: soil4_hr(:)
real(r8), pointer :: col_ctrunc(:) ! (gC/m2) column-level sink for C truncation
integer , pointer :: ivt(:) ! pft vegetation type
real(r8), pointer :: deadcrootc_xfer_to_deadcrootc(:)
real(r8), pointer :: deadstemc_xfer_to_deadstemc(:)
real(r8), pointer :: frootc_xfer_to_frootc(:)
real(r8), pointer :: leafc_xfer_to_leafc(:)
real(r8), pointer :: livecrootc_xfer_to_livecrootc(:)
real(r8), pointer :: livestemc_xfer_to_livestemc(:)
real(r8), pointer :: cpool_to_xsmrpool(:)
real(r8), pointer :: cpool_to_deadcrootc(:)
real(r8), pointer :: cpool_to_deadcrootc_storage(:)
real(r8), pointer :: cpool_to_deadstemc(:)
real(r8), pointer :: cpool_to_deadstemc_storage(:)
real(r8), pointer :: cpool_to_frootc(:)
real(r8), pointer :: cpool_to_frootc_storage(:)
real(r8), pointer :: cpool_to_gresp_storage(:)
real(r8), pointer :: cpool_to_leafc(:)
real(r8), pointer :: cpool_to_leafc_storage(:)
real(r8), pointer :: cpool_to_livecrootc(:)
real(r8), pointer :: cpool_to_livecrootc_storage(:)
real(r8), pointer :: cpool_to_livestemc(:)
real(r8), pointer :: cpool_to_livestemc_storage(:)
real(r8), pointer :: deadcrootc_storage_to_xfer(:)
real(r8), pointer :: deadstemc_storage_to_xfer(:)
real(r8), pointer :: frootc_storage_to_xfer(:)
real(r8), pointer :: frootc_to_litter(:)
real(r8), pointer :: gresp_storage_to_xfer(:)
real(r8), pointer :: leafc_storage_to_xfer(:)
real(r8), pointer :: leafc_to_litter(:)
real(r8), pointer :: livecrootc_storage_to_xfer(:)
real(r8), pointer :: livecrootc_to_deadcrootc(:)
real(r8), pointer :: livestemc_storage_to_xfer(:)
real(r8), pointer :: livestemc_to_deadstemc(:)
real(r8), pointer :: livestem_mr(:)
real(r8), pointer :: froot_mr(:)
real(r8), pointer :: leaf_mr(:)
real(r8), pointer :: livecroot_mr(:)
real(r8), pointer :: livestem_curmr(:)
real(r8), pointer :: froot_curmr(:)
real(r8), pointer :: leaf_curmr(:)
real(r8), pointer :: livecroot_curmr(:)
real(r8), pointer :: livestem_xsmr(:)
real(r8), pointer :: froot_xsmr(:)
real(r8), pointer :: leaf_xsmr(:)
real(r8), pointer :: livecroot_xsmr(:)
real(r8), pointer :: cpool_deadcroot_gr(:)
real(r8), pointer :: cpool_deadcroot_storage_gr(:)
real(r8), pointer :: cpool_deadstem_gr(:)
real(r8), pointer :: cpool_deadstem_storage_gr(:)
real(r8), pointer :: cpool_froot_gr(:)
real(r8), pointer :: cpool_froot_storage_gr(:)
real(r8), pointer :: cpool_leaf_gr(:)
real(r8), pointer :: cpool_leaf_storage_gr(:)
real(r8), pointer :: cpool_livecroot_gr(:)
real(r8), pointer :: cpool_livecroot_storage_gr(:)
real(r8), pointer :: cpool_livestem_gr(:)
real(r8), pointer :: cpool_livestem_storage_gr(:)
real(r8), pointer :: transfer_deadcroot_gr(:)
real(r8), pointer :: transfer_deadstem_gr(:)
real(r8), pointer :: transfer_froot_gr(:)
real(r8), pointer :: transfer_leaf_gr(:)
real(r8), pointer :: transfer_livecroot_gr(:)
real(r8), pointer :: transfer_livestem_gr(:)
!
! local pointers to implicit in/out arrays
real(r8), pointer :: cwdc(:) ! (gC/m2) coarse woody debris C
real(r8), pointer :: litr1c(:) ! (gC/m2) litter labile C
real(r8), pointer :: litr2c(:) ! (gC/m2) litter cellulose C
real(r8), pointer :: litr3c(:) ! (gC/m2) litter lignin C
real(r8), pointer :: soil1c(:) ! (gC/m2) soil organic matter C (fast pool)
real(r8), pointer :: soil2c(:) ! (gC/m2) soil organic matter C (medium pool)
real(r8), pointer :: soil3c(:) ! (gC/m2) soil organic matter C (slow pool)
real(r8), pointer :: soil4c(:) ! (gC/m2) soil organic matter C (slowest pool)
real(r8), pointer :: cpool(:) ! (gC/m2) temporary photosynthate C pool
real(r8), pointer :: xsmrpool(:) ! (gC/m2) execss maint resp C pool
real(r8), pointer :: deadcrootc(:) ! (gC/m2) dead coarse root C
real(r8), pointer :: deadcrootc_storage(:) ! (gC/m2) dead coarse root C storage
real(r8), pointer :: deadcrootc_xfer(:) ! (gC/m2) dead coarse root C transfer
real(r8), pointer :: deadstemc(:) ! (gC/m2) dead stem C
real(r8), pointer :: deadstemc_storage(:) ! (gC/m2) dead stem C storage
real(r8), pointer :: deadstemc_xfer(:) ! (gC/m2) dead stem C transfer
real(r8), pointer :: frootc(:) ! (gC/m2) fine root C
real(r8), pointer :: frootc_storage(:) ! (gC/m2) fine root C storage
real(r8), pointer :: frootc_xfer(:) ! (gC/m2) fine root C transfer
real(r8), pointer :: gresp_storage(:) ! (gC/m2) growth respiration storage
real(r8), pointer :: gresp_xfer(:) ! (gC/m2) growth respiration transfer
real(r8), pointer :: leafc(:) ! (gC/m2) leaf C
real(r8), pointer :: leafc_storage(:) ! (gC/m2) leaf C storage
real(r8), pointer :: leafc_xfer(:) ! (gC/m2) leaf C transfer
real(r8), pointer :: livecrootc(:) ! (gC/m2) live coarse root C
real(r8), pointer :: livecrootc_storage(:) ! (gC/m2) live coarse root C storage
real(r8), pointer :: livecrootc_xfer(:) ! (gC/m2) live coarse root C transfer
real(r8), pointer :: livestemc(:) ! (gC/m2) live stem C
real(r8), pointer :: livestemc_storage(:) ! (gC/m2) live stem C storage
real(r8), pointer :: livestemc_xfer(:) ! (gC/m2) live stem C transfer
real(r8), pointer :: pft_ctrunc(:) ! (gC/m2) pft-level sink for C truncation
! local pointers for dynamic landcover fluxes and states
real(r8), pointer :: dwt_seedc_to_leaf(:)
real(r8), pointer :: dwt_seedc_to_deadstem(:)
real(r8), pointer :: dwt_frootc_to_litr1c(:)
real(r8), pointer :: dwt_frootc_to_litr2c(:)
real(r8), pointer :: dwt_frootc_to_litr3c(:)
real(r8), pointer :: dwt_livecrootc_to_cwdc(:)
real(r8), pointer :: dwt_deadcrootc_to_cwdc(:)
real(r8), pointer :: seedc(:)
!
! !OTHER LOCAL VARIABLES:
integer :: c,p ! indices
integer :: fp,fc ! lake filter indices
real(r8):: dt ! radiation time step (seconds)
!
!EOP
!-----------------------------------------------------------------------
! assign local pointers
woody => pftcon%woody
! assign local pointers at the column level
cwdc_to_litr2c => clm3%g%l%c%cc13f%cwdc_to_litr2c
cwdc_to_litr3c => clm3%g%l%c%cc13f%cwdc_to_litr3c
frootc_to_litr1c => clm3%g%l%c%cc13f%frootc_to_litr1c
frootc_to_litr2c => clm3%g%l%c%cc13f%frootc_to_litr2c
frootc_to_litr3c => clm3%g%l%c%cc13f%frootc_to_litr3c
leafc_to_litr1c => clm3%g%l%c%cc13f%leafc_to_litr1c
leafc_to_litr2c => clm3%g%l%c%cc13f%leafc_to_litr2c
leafc_to_litr3c => clm3%g%l%c%cc13f%leafc_to_litr3c
litr1_hr => clm3%g%l%c%cc13f%litr1_hr
litr1c_to_soil1c => clm3%g%l%c%cc13f%litr1c_to_soil1c
litr2_hr => clm3%g%l%c%cc13f%litr2_hr
litr2c_to_soil2c => clm3%g%l%c%cc13f%litr2c_to_soil2c
litr3_hr => clm3%g%l%c%cc13f%litr3_hr
litr3c_to_soil3c => clm3%g%l%c%cc13f%litr3c_to_soil3c
soil1_hr => clm3%g%l%c%cc13f%soil1_hr
soil1c_to_soil2c => clm3%g%l%c%cc13f%soil1c_to_soil2c
soil2_hr => clm3%g%l%c%cc13f%soil2_hr
soil2c_to_soil3c => clm3%g%l%c%cc13f%soil2c_to_soil3c
soil3_hr => clm3%g%l%c%cc13f%soil3_hr
soil3c_to_soil4c => clm3%g%l%c%cc13f%soil3c_to_soil4c
soil4_hr => clm3%g%l%c%cc13f%soil4_hr
col_ctrunc => clm3%g%l%c%cc13s%col_ctrunc
cwdc => clm3%g%l%c%cc13s%cwdc
litr1c => clm3%g%l%c%cc13s%litr1c
litr2c => clm3%g%l%c%cc13s%litr2c
litr3c => clm3%g%l%c%cc13s%litr3c
soil1c => clm3%g%l%c%cc13s%soil1c
soil2c => clm3%g%l%c%cc13s%soil2c
soil3c => clm3%g%l%c%cc13s%soil3c
soil4c => clm3%g%l%c%cc13s%soil4c
! new pointers for dynamic landcover
dwt_seedc_to_leaf => clm3%g%l%c%cc13f%dwt_seedc_to_leaf
dwt_seedc_to_deadstem => clm3%g%l%c%cc13f%dwt_seedc_to_deadstem
dwt_frootc_to_litr1c => clm3%g%l%c%cc13f%dwt_frootc_to_litr1c
dwt_frootc_to_litr2c => clm3%g%l%c%cc13f%dwt_frootc_to_litr2c
dwt_frootc_to_litr3c => clm3%g%l%c%cc13f%dwt_frootc_to_litr3c
dwt_livecrootc_to_cwdc => clm3%g%l%c%cc13f%dwt_livecrootc_to_cwdc
dwt_deadcrootc_to_cwdc => clm3%g%l%c%cc13f%dwt_deadcrootc_to_cwdc
seedc => clm3%g%l%c%cc13s%seedc
! assign local pointers at the pft level
ivt => clm3%g%l%c%p%itype
cpool_deadcroot_gr => clm3%g%l%c%p%pc13f%cpool_deadcroot_gr
cpool_deadcroot_storage_gr => clm3%g%l%c%p%pc13f%cpool_deadcroot_storage_gr
cpool_deadstem_gr => clm3%g%l%c%p%pc13f%cpool_deadstem_gr
cpool_deadstem_storage_gr => clm3%g%l%c%p%pc13f%cpool_deadstem_storage_gr
cpool_froot_gr => clm3%g%l%c%p%pc13f%cpool_froot_gr
cpool_froot_storage_gr => clm3%g%l%c%p%pc13f%cpool_froot_storage_gr
cpool_leaf_gr => clm3%g%l%c%p%pc13f%cpool_leaf_gr
cpool_leaf_storage_gr => clm3%g%l%c%p%pc13f%cpool_leaf_storage_gr
cpool_livecroot_gr => clm3%g%l%c%p%pc13f%cpool_livecroot_gr
cpool_livecroot_storage_gr => clm3%g%l%c%p%pc13f%cpool_livecroot_storage_gr
cpool_livestem_gr => clm3%g%l%c%p%pc13f%cpool_livestem_gr
cpool_livestem_storage_gr => clm3%g%l%c%p%pc13f%cpool_livestem_storage_gr
cpool_to_xsmrpool => clm3%g%l%c%p%pc13f%cpool_to_xsmrpool
cpool_to_deadcrootc => clm3%g%l%c%p%pc13f%cpool_to_deadcrootc
cpool_to_deadcrootc_storage => clm3%g%l%c%p%pc13f%cpool_to_deadcrootc_storage
cpool_to_deadstemc => clm3%g%l%c%p%pc13f%cpool_to_deadstemc
cpool_to_deadstemc_storage => clm3%g%l%c%p%pc13f%cpool_to_deadstemc_storage
cpool_to_frootc => clm3%g%l%c%p%pc13f%cpool_to_frootc
cpool_to_frootc_storage => clm3%g%l%c%p%pc13f%cpool_to_frootc_storage
cpool_to_gresp_storage => clm3%g%l%c%p%pc13f%cpool_to_gresp_storage
cpool_to_leafc => clm3%g%l%c%p%pc13f%cpool_to_leafc
cpool_to_leafc_storage => clm3%g%l%c%p%pc13f%cpool_to_leafc_storage
cpool_to_livecrootc => clm3%g%l%c%p%pc13f%cpool_to_livecrootc
cpool_to_livecrootc_storage => clm3%g%l%c%p%pc13f%cpool_to_livecrootc_storage
cpool_to_livestemc => clm3%g%l%c%p%pc13f%cpool_to_livestemc
cpool_to_livestemc_storage => clm3%g%l%c%p%pc13f%cpool_to_livestemc_storage
deadcrootc_storage_to_xfer => clm3%g%l%c%p%pc13f%deadcrootc_storage_to_xfer
deadcrootc_xfer_to_deadcrootc => clm3%g%l%c%p%pc13f%deadcrootc_xfer_to_deadcrootc
deadstemc_storage_to_xfer => clm3%g%l%c%p%pc13f%deadstemc_storage_to_xfer
deadstemc_xfer_to_deadstemc => clm3%g%l%c%p%pc13f%deadstemc_xfer_to_deadstemc
froot_mr => clm3%g%l%c%p%pc13f%froot_mr
froot_curmr => clm3%g%l%c%p%pc13f%froot_curmr
froot_xsmr => clm3%g%l%c%p%pc13f%froot_xsmr
frootc_storage_to_xfer => clm3%g%l%c%p%pc13f%frootc_storage_to_xfer
frootc_to_litter => clm3%g%l%c%p%pc13f%frootc_to_litter
frootc_xfer_to_frootc => clm3%g%l%c%p%pc13f%frootc_xfer_to_frootc
gresp_storage_to_xfer => clm3%g%l%c%p%pc13f%gresp_storage_to_xfer
leaf_mr => clm3%g%l%c%p%pc13f%leaf_mr
leaf_curmr => clm3%g%l%c%p%pc13f%leaf_curmr
leaf_xsmr => clm3%g%l%c%p%pc13f%leaf_xsmr
leafc_storage_to_xfer => clm3%g%l%c%p%pc13f%leafc_storage_to_xfer
leafc_to_litter => clm3%g%l%c%p%pc13f%leafc_to_litter
leafc_xfer_to_leafc => clm3%g%l%c%p%pc13f%leafc_xfer_to_leafc
livecroot_mr => clm3%g%l%c%p%pc13f%livecroot_mr
livecroot_curmr => clm3%g%l%c%p%pc13f%livecroot_curmr
livecroot_xsmr => clm3%g%l%c%p%pc13f%livecroot_xsmr
livecrootc_storage_to_xfer => clm3%g%l%c%p%pc13f%livecrootc_storage_to_xfer
livecrootc_to_deadcrootc => clm3%g%l%c%p%pc13f%livecrootc_to_deadcrootc
livecrootc_xfer_to_livecrootc => clm3%g%l%c%p%pc13f%livecrootc_xfer_to_livecrootc
livestem_mr => clm3%g%l%c%p%pc13f%livestem_mr
livestem_curmr => clm3%g%l%c%p%pc13f%livestem_curmr
livestem_xsmr => clm3%g%l%c%p%pc13f%livestem_xsmr
livestemc_storage_to_xfer => clm3%g%l%c%p%pc13f%livestemc_storage_to_xfer
livestemc_to_deadstemc => clm3%g%l%c%p%pc13f%livestemc_to_deadstemc
livestemc_xfer_to_livestemc => clm3%g%l%c%p%pc13f%livestemc_xfer_to_livestemc
transfer_deadcroot_gr => clm3%g%l%c%p%pc13f%transfer_deadcroot_gr
transfer_deadstem_gr => clm3%g%l%c%p%pc13f%transfer_deadstem_gr
transfer_froot_gr => clm3%g%l%c%p%pc13f%transfer_froot_gr
transfer_leaf_gr => clm3%g%l%c%p%pc13f%transfer_leaf_gr
transfer_livecroot_gr => clm3%g%l%c%p%pc13f%transfer_livecroot_gr
transfer_livestem_gr => clm3%g%l%c%p%pc13f%transfer_livestem_gr
cpool => clm3%g%l%c%p%pc13s%cpool
xsmrpool => clm3%g%l%c%p%pc13s%xsmrpool
deadcrootc => clm3%g%l%c%p%pc13s%deadcrootc
deadcrootc_storage => clm3%g%l%c%p%pc13s%deadcrootc_storage
deadcrootc_xfer => clm3%g%l%c%p%pc13s%deadcrootc_xfer
deadstemc => clm3%g%l%c%p%pc13s%deadstemc
deadstemc_storage => clm3%g%l%c%p%pc13s%deadstemc_storage
deadstemc_xfer => clm3%g%l%c%p%pc13s%deadstemc_xfer
frootc => clm3%g%l%c%p%pc13s%frootc
frootc_storage => clm3%g%l%c%p%pc13s%frootc_storage
frootc_xfer => clm3%g%l%c%p%pc13s%frootc_xfer
gresp_storage => clm3%g%l%c%p%pc13s%gresp_storage
gresp_xfer => clm3%g%l%c%p%pc13s%gresp_xfer
leafc => clm3%g%l%c%p%pc13s%leafc
leafc_storage => clm3%g%l%c%p%pc13s%leafc_storage
leafc_xfer => clm3%g%l%c%p%pc13s%leafc_xfer
livecrootc => clm3%g%l%c%p%pc13s%livecrootc
livecrootc_storage => clm3%g%l%c%p%pc13s%livecrootc_storage
livecrootc_xfer => clm3%g%l%c%p%pc13s%livecrootc_xfer
livestemc => clm3%g%l%c%p%pc13s%livestemc
livestemc_storage => clm3%g%l%c%p%pc13s%livestemc_storage
livestemc_xfer => clm3%g%l%c%p%pc13s%livestemc_xfer
pft_ctrunc => clm3%g%l%c%p%pc13s%pft_ctrunc
! set time steps
dt = real( get_step_size(), r8 )
! column loop
do fc = 1,num_soilc
c = filter_soilc(fc)
! column level fluxes
! plant to litter fluxes
! leaf litter
litr1c(c) = litr1c(c) + leafc_to_litr1c(c)*dt
litr2c(c) = litr2c(c) + leafc_to_litr2c(c)*dt
litr3c(c) = litr3c(c) + leafc_to_litr3c(c)*dt
! fine root litter
litr1c(c) = litr1c(c) + frootc_to_litr1c(c)*dt
litr2c(c) = litr2c(c) + frootc_to_litr2c(c)*dt
litr3c(c) = litr3c(c) + frootc_to_litr3c(c)*dt
! seeding fluxes, from dynamic landcover
seedc(c) = seedc(c) - dwt_seedc_to_leaf(c) * dt
seedc(c) = seedc(c) - dwt_seedc_to_deadstem(c) * dt
! fluxes into litter and CWD, from dynamic landcover
litr1c(c) = litr1c(c) + dwt_frootc_to_litr1c(c)*dt
litr2c(c) = litr2c(c) + dwt_frootc_to_litr2c(c)*dt
litr3c(c) = litr3c(c) + dwt_frootc_to_litr3c(c)*dt
cwdc(c) = cwdc(c) + dwt_livecrootc_to_cwdc(c)*dt
cwdc(c) = cwdc(c) + dwt_deadcrootc_to_cwdc(c)*dt
! litter and SOM HR fluxes
litr1c(c) = litr1c(c) - litr1_hr(c)*dt
litr2c(c) = litr2c(c) - litr2_hr(c)*dt
litr3c(c) = litr3c(c) - litr3_hr(c)*dt
soil1c(c) = soil1c(c) - soil1_hr(c)*dt
soil2c(c) = soil2c(c) - soil2_hr(c)*dt
soil3c(c) = soil3c(c) - soil3_hr(c)*dt
soil4c(c) = soil4c(c) - soil4_hr(c)*dt
! CWD to litter fluxes
cwdc(c) = cwdc(c) - cwdc_to_litr2c(c)*dt
litr2c(c) = litr2c(c) + cwdc_to_litr2c(c)*dt
cwdc(c) = cwdc(c) - cwdc_to_litr3c(c)*dt
litr3c(c) = litr3c(c) + cwdc_to_litr3c(c)*dt
! litter to SOM fluxes
litr1c(c) = litr1c(c) - litr1c_to_soil1c(c)*dt
soil1c(c) = soil1c(c) + litr1c_to_soil1c(c)*dt
litr2c(c) = litr2c(c) - litr2c_to_soil2c(c)*dt
soil2c(c) = soil2c(c) + litr2c_to_soil2c(c)*dt
litr3c(c) = litr3c(c) - litr3c_to_soil3c(c)*dt
soil3c(c) = soil3c(c) + litr3c_to_soil3c(c)*dt
! SOM to SOM fluxes
soil1c(c) = soil1c(c) - soil1c_to_soil2c(c)*dt
soil2c(c) = soil2c(c) + soil1c_to_soil2c(c)*dt
soil2c(c) = soil2c(c) - soil2c_to_soil3c(c)*dt
soil3c(c) = soil3c(c) + soil2c_to_soil3c(c)*dt
soil3c(c) = soil3c(c) - soil3c_to_soil4c(c)*dt
soil4c(c) = soil4c(c) + soil3c_to_soil4c(c)*dt
end do ! end of columns loop
! pft loop
do fp = 1,num_soilp
p = filter_soilp(fp)
! phenology: transfer growth fluxes
leafc(p) = leafc(p) + leafc_xfer_to_leafc(p)*dt
leafc_xfer(p) = leafc_xfer(p) - leafc_xfer_to_leafc(p)*dt
frootc(p) = frootc(p) + frootc_xfer_to_frootc(p)*dt
frootc_xfer(p) = frootc_xfer(p) - frootc_xfer_to_frootc(p)*dt
if (woody(ivt(p)) == 1._r8) then
livestemc(p) = livestemc(p) + livestemc_xfer_to_livestemc(p)*dt
livestemc_xfer(p) = livestemc_xfer(p) - livestemc_xfer_to_livestemc(p)*dt
deadstemc(p) = deadstemc(p) + deadstemc_xfer_to_deadstemc(p)*dt
deadstemc_xfer(p) = deadstemc_xfer(p) - deadstemc_xfer_to_deadstemc(p)*dt
livecrootc(p) = livecrootc(p) + livecrootc_xfer_to_livecrootc(p)*dt
livecrootc_xfer(p) = livecrootc_xfer(p) - livecrootc_xfer_to_livecrootc(p)*dt
deadcrootc(p) = deadcrootc(p) + deadcrootc_xfer_to_deadcrootc(p)*dt
deadcrootc_xfer(p) = deadcrootc_xfer(p) - deadcrootc_xfer_to_deadcrootc(p)*dt
end if
! phenology: litterfall fluxes
leafc(p) = leafc(p) - leafc_to_litter(p)*dt
frootc(p) = frootc(p) - frootc_to_litter(p)*dt
! livewood turnover fluxes
if (woody(ivt(p)) == 1._r8) then
livestemc(p) = livestemc(p) - livestemc_to_deadstemc(p)*dt
deadstemc(p) = deadstemc(p) + livestemc_to_deadstemc(p)*dt
livecrootc(p) = livecrootc(p) - livecrootc_to_deadcrootc(p)*dt
deadcrootc(p) = deadcrootc(p) + livecrootc_to_deadcrootc(p)*dt
end if
! maintenance respiration fluxes
cpool(p) = cpool(p) - cpool_to_xsmrpool(p)*dt
cpool(p) = cpool(p) - leaf_curmr(p)*dt
cpool(p) = cpool(p) - froot_curmr(p)*dt
if (woody(ivt(p)) == 1._r8) then
cpool(p) = cpool(p) - livestem_curmr(p)*dt
cpool(p) = cpool(p) - livecroot_curmr(p)*dt
end if
! maintenance respiration fluxes
xsmrpool(p) = xsmrpool(p) + cpool_to_xsmrpool(p)*dt
xsmrpool(p) = xsmrpool(p) - leaf_xsmr(p)*dt
xsmrpool(p) = xsmrpool(p) - froot_xsmr(p)*dt
if (woody(ivt(p)) == 1._r8) then
xsmrpool(p) = xsmrpool(p) - livestem_xsmr(p)*dt
xsmrpool(p) = xsmrpool(p) - livecroot_xsmr(p)*dt
end if
! allocation fluxes
cpool(p) = cpool(p) - cpool_to_leafc(p)*dt
leafc(p) = leafc(p) + cpool_to_leafc(p)*dt
cpool(p) = cpool(p) - cpool_to_leafc_storage(p)*dt
leafc_storage(p) = leafc_storage(p) + cpool_to_leafc_storage(p)*dt
cpool(p) = cpool(p) - cpool_to_frootc(p)*dt
frootc(p) = frootc(p) + cpool_to_frootc(p)*dt
cpool(p) = cpool(p) - cpool_to_frootc_storage(p)*dt
frootc_storage(p) = frootc_storage(p) + cpool_to_frootc_storage(p)*dt
if (woody(ivt(p)) == 1._r8) then
cpool(p) = cpool(p) - cpool_to_livestemc(p)*dt
livestemc(p) = livestemc(p) + cpool_to_livestemc(p)*dt
cpool(p) = cpool(p) - cpool_to_livestemc_storage(p)*dt
livestemc_storage(p) = livestemc_storage(p) + cpool_to_livestemc_storage(p)*dt
cpool(p) = cpool(p) - cpool_to_deadstemc(p)*dt
deadstemc(p) = deadstemc(p) + cpool_to_deadstemc(p)*dt
cpool(p) = cpool(p) - cpool_to_deadstemc_storage(p)*dt
deadstemc_storage(p) = deadstemc_storage(p) + cpool_to_deadstemc_storage(p)*dt
cpool(p) = cpool(p) - cpool_to_livecrootc(p)*dt
livecrootc(p) = livecrootc(p) + cpool_to_livecrootc(p)*dt
cpool(p) = cpool(p) - cpool_to_livecrootc_storage(p)*dt
livecrootc_storage(p) = livecrootc_storage(p) + cpool_to_livecrootc_storage(p)*dt
cpool(p) = cpool(p) - cpool_to_deadcrootc(p)*dt
deadcrootc(p) = deadcrootc(p) + cpool_to_deadcrootc(p)*dt
cpool(p) = cpool(p) - cpool_to_deadcrootc_storage(p)*dt
deadcrootc_storage(p) = deadcrootc_storage(p) + cpool_to_deadcrootc_storage(p)*dt
end if
! growth respiration fluxes for current growth
cpool(p) = cpool(p) - cpool_leaf_gr(p)*dt
cpool(p) = cpool(p) - cpool_froot_gr(p)*dt
if (woody(ivt(p)) == 1._r8) then
cpool(p) = cpool(p) - cpool_livestem_gr(p)*dt
cpool(p) = cpool(p) - cpool_deadstem_gr(p)*dt
cpool(p) = cpool(p) - cpool_livecroot_gr(p)*dt
cpool(p) = cpool(p) - cpool_deadcroot_gr(p)*dt
end if
! growth respiration for transfer growth
gresp_xfer(p) = gresp_xfer(p) - transfer_leaf_gr(p)*dt
gresp_xfer(p) = gresp_xfer(p) - transfer_froot_gr(p)*dt
if (woody(ivt(p)) == 1._r8) then
gresp_xfer(p) = gresp_xfer(p) - transfer_livestem_gr(p)*dt
gresp_xfer(p) = gresp_xfer(p) - transfer_deadstem_gr(p)*dt
gresp_xfer(p) = gresp_xfer(p) - transfer_livecroot_gr(p)*dt
gresp_xfer(p) = gresp_xfer(p) - transfer_deadcroot_gr(p)*dt
end if
! growth respiration at time of storage
cpool(p) = cpool(p) - cpool_leaf_storage_gr(p)*dt
cpool(p) = cpool(p) - cpool_froot_storage_gr(p)*dt
if (woody(ivt(p)) == 1._r8) then
cpool(p) = cpool(p) - cpool_livestem_storage_gr(p)*dt
cpool(p) = cpool(p) - cpool_deadstem_storage_gr(p)*dt
cpool(p) = cpool(p) - cpool_livecroot_storage_gr(p)*dt
cpool(p) = cpool(p) - cpool_deadcroot_storage_gr(p)*dt
end if
! growth respiration stored for release during transfer growth
cpool(p) = cpool(p) - cpool_to_gresp_storage(p)*dt
gresp_storage(p) = gresp_storage(p) + cpool_to_gresp_storage(p)*dt
! move storage pools into transfer pools
leafc_storage(p) = leafc_storage(p) - leafc_storage_to_xfer(p)*dt
leafc_xfer(p) = leafc_xfer(p) + leafc_storage_to_xfer(p)*dt
frootc_storage(p) = frootc_storage(p) - frootc_storage_to_xfer(p)*dt
frootc_xfer(p) = frootc_xfer(p) + frootc_storage_to_xfer(p)*dt
if (woody(ivt(p)) == 1._r8) then
livestemc_storage(p) = livestemc_storage(p) - livestemc_storage_to_xfer(p)*dt
livestemc_xfer(p) = livestemc_xfer(p) + livestemc_storage_to_xfer(p)*dt
deadstemc_storage(p) = deadstemc_storage(p) - deadstemc_storage_to_xfer(p)*dt
deadstemc_xfer(p) = deadstemc_xfer(p) + deadstemc_storage_to_xfer(p)*dt
livecrootc_storage(p) = livecrootc_storage(p) - livecrootc_storage_to_xfer(p)*dt
livecrootc_xfer(p) = livecrootc_xfer(p) + livecrootc_storage_to_xfer(p)*dt
deadcrootc_storage(p) = deadcrootc_storage(p) - deadcrootc_storage_to_xfer(p)*dt
deadcrootc_xfer(p) = deadcrootc_xfer(p) + deadcrootc_storage_to_xfer(p)*dt
gresp_storage(p) = gresp_storage(p) - gresp_storage_to_xfer(p)*dt
gresp_xfer(p) = gresp_xfer(p) + gresp_storage_to_xfer(p)*dt
end if
end do ! end of pft loop
end subroutine C13StateUpdate1
!-----------------------------------------------------------------------
#endif
end module CNC13StateUpdate1Mod