Skip to content

Instantly share code, notes, and snippets.

Show Gist options
  • Save bishtgautam/035ad3e622a80359bb43f66b720de699 to your computer and use it in GitHub Desktop.
Save bishtgautam/035ad3e622a80359bb43f66b720de699 to your computer and use it in GitHub Desktop.
diff --git a/components/elm/src/biogeochem/CarbonStateUpdate1Mod.F90 b/components/elm/src/biogeochem/CarbonStateUpdate1Mod.F90
index 057fd0e56b..8ee384c6ba 100644
--- a/components/elm/src/biogeochem/CarbonStateUpdate1Mod.F90
+++ b/components/elm/src/biogeochem/CarbonStateUpdate1Mod.F90
@@ -4,6 +4,7 @@ module CarbonStateUpdate1Mod
! Module for carbon state variable update, non-mortality fluxes.
!
! !USES:
+ use elm_varctl , only : iulog
use shr_kind_mod , only : r8 => shr_kind_r8
use decompMod , only : bounds_type
use elm_varpar , only : ndecomp_cascade_transitions, nlevdecomp
@@ -277,8 +278,16 @@ contains
if (woody(ivt(p)) == 1._r8) then
veg_cs%livestemc(p) = veg_cs%livestemc(p) + veg_cf%livestemc_xfer_to_livestemc(p)*dt
veg_cs%livestemc_xfer(p) = veg_cs%livestemc_xfer(p) - veg_cf%livestemc_xfer_to_livestemc(p)*dt
+ if (p == 997671) then
+ write(iulog,*)''
+ write(iulog,*)'0.',veg_cs%deadstemc(p)
+ endif
veg_cs%deadstemc(p) = veg_cs%deadstemc(p) + veg_cf%deadstemc_xfer_to_deadstemc(p)*dt
veg_cs%deadstemc_xfer(p) = veg_cs%deadstemc_xfer(p) - veg_cf%deadstemc_xfer_to_deadstemc(p)*dt
+ if (p == 997671) then
+ write(iulog,*) '1.', &
+ veg_cs%deadstemc(p),veg_cs%deadstemc_xfer(p),veg_cf%deadstemc_xfer_to_deadstemc(p)*dt
+ endif
veg_cs%livecrootc(p) = veg_cs%livecrootc(p) + veg_cf%livecrootc_xfer_to_livecrootc(p)*dt
veg_cs%livecrootc_xfer(p) = veg_cs%livecrootc_xfer(p) - veg_cf%livecrootc_xfer_to_livecrootc(p)*dt
veg_cs%deadcrootc(p) = veg_cs%deadcrootc(p) + veg_cf%deadcrootc_xfer_to_deadcrootc(p)*dt
@@ -300,6 +309,10 @@ contains
if (woody(ivt(p)) == 1._r8) then
veg_cs%livestemc(p) = veg_cs%livestemc(p) - veg_cf%livestemc_to_deadstemc(p)*dt
veg_cs%deadstemc(p) = veg_cs%deadstemc(p) + veg_cf%livestemc_to_deadstemc(p)*dt
+ if (p == 997671) then
+ write(iulog,*) '2.', &
+ veg_cs%deadstemc(p),veg_cf%livestemc_to_deadstemc(p)*dt
+ endif
veg_cs%livecrootc(p) = veg_cs%livecrootc(p) - veg_cf%livecrootc_to_deadcrootc(p)*dt
veg_cs%deadcrootc(p) = veg_cs%deadcrootc(p) + veg_cf%livecrootc_to_deadcrootc(p)*dt
end if
@@ -362,6 +375,10 @@ contains
veg_cs%livestemc_storage(p) = veg_cs%livestemc_storage(p) + veg_cf%cpool_to_livestemc_storage(p)*dt
veg_cs%cpool(p) = veg_cs%cpool(p) - veg_cf%cpool_to_deadstemc(p)*dt
veg_cs%deadstemc(p) = veg_cs%deadstemc(p) + veg_cf%cpool_to_deadstemc(p)*dt
+ if (p == 997671) then
+ write(iulog,*) '3.', &
+ veg_cs%deadstemc(p),veg_cf%cpool_to_deadstemc(p)*dt
+ endif
veg_cs%cpool(p) = veg_cs%cpool(p) - veg_cf%cpool_to_deadstemc_storage(p)*dt
veg_cs%deadstemc_storage(p) = veg_cs%deadstemc_storage(p) + veg_cf%cpool_to_deadstemc_storage(p)*dt
veg_cs%cpool(p) = veg_cs%cpool(p) - veg_cf%cpool_to_livecrootc(p)*dt
diff --git a/components/elm/src/biogeochem/CarbonStateUpdate2Mod.F90 b/components/elm/src/biogeochem/CarbonStateUpdate2Mod.F90
index 48174fb8dd..60592bd175 100644
--- a/components/elm/src/biogeochem/CarbonStateUpdate2Mod.F90
+++ b/components/elm/src/biogeochem/CarbonStateUpdate2Mod.F90
@@ -5,6 +5,7 @@ module CarbonStateUpdate2Mod
! Module for carbon state variable update, mortality fluxes.
!
! !USES:
+ use elm_varctl , only : iulog
use shr_kind_mod , only : r8 => shr_kind_r8
use shr_log_mod , only : errMsg => shr_log_errMsg
use abortutils , only : endrun
@@ -90,6 +91,10 @@ contains
veg_cs%frootc(p) = veg_cs%frootc(p) - veg_cf%m_frootc_to_litter(p) * dt
veg_cs%livestemc(p) = veg_cs%livestemc(p) - veg_cf%m_livestemc_to_litter(p) * dt
veg_cs%deadstemc(p) = veg_cs%deadstemc(p) - veg_cf%m_deadstemc_to_litter(p) * dt
+ if (p == 997671) then
+ write(iulog,*) '4.', &
+ veg_cs%deadstemc(p),veg_cf%m_deadstemc_to_litter(p)*dt
+ endif
veg_cs%livecrootc(p) = veg_cs%livecrootc(p) - veg_cf%m_livecrootc_to_litter(p) * dt
veg_cs%deadcrootc(p) = veg_cs%deadcrootc(p) - veg_cf%m_deadcrootc_to_litter(p) * dt
! storage pools
@@ -181,7 +186,15 @@ contains
veg_cs%frootc(p) = veg_cs%frootc(p) - veg_cf%hrv_frootc_to_litter(p) * dt
veg_cs%livestemc(p) = veg_cs%livestemc(p) - veg_cf%hrv_livestemc_to_litter(p) * dt
veg_cs%deadstemc(p) = veg_cs%deadstemc(p) - veg_cf%hrv_deadstemc_to_prod10c(p) * dt
+ if (p == 997671) then
+ write(iulog,*) '5.', &
+ veg_cs%deadstemc(p),veg_cf%hrv_deadstemc_to_prod10c(p)*dt
+ endif
veg_cs%deadstemc(p) = veg_cs%deadstemc(p) - veg_cf%hrv_deadstemc_to_prod100c(p) * dt
+ if (p == 997671) then
+ write(iulog,*) '6.', &
+ veg_cs%deadstemc(p),veg_cf%hrv_deadstemc_to_prod100c(p)*dt
+ endif
veg_cs%livecrootc(p) = veg_cs%livecrootc(p) - veg_cf%hrv_livecrootc_to_litter(p) * dt
veg_cs%deadcrootc(p) = veg_cs%deadcrootc(p) - veg_cf%hrv_deadcrootc_to_litter(p) * dt
diff --git a/components/elm/src/biogeochem/CarbonStateUpdate3Mod.F90 b/components/elm/src/biogeochem/CarbonStateUpdate3Mod.F90
index f8d8c9e196..bbe86ec3d4 100644
--- a/components/elm/src/biogeochem/CarbonStateUpdate3Mod.F90
+++ b/components/elm/src/biogeochem/CarbonStateUpdate3Mod.F90
@@ -5,6 +5,7 @@ module CarbonStateUpdate3Mod
! Module for carbon state variable update, mortality fluxes.
!
! !USES:
+ use elm_varctl , only : iulog
use shr_kind_mod , only : r8 => shr_kind_r8
use shr_log_mod , only : errMsg => shr_log_errMsg
use abortutils , only : endrun
@@ -115,7 +116,16 @@ contains
veg_cs%livestemc(p) = veg_cs%livestemc(p) - veg_cf%m_livestemc_to_fire(p) * dt
veg_cs%livestemc(p) = veg_cs%livestemc(p) - veg_cf%m_livestemc_to_litter_fire(p) * dt
veg_cs%deadstemc(p) = veg_cs%deadstemc(p) - veg_cf%m_deadstemc_to_fire(p) * dt
+ if (p == 997671) then
+ write(iulog,*) '7.', &
+ veg_cs%deadstemc(p),veg_cf%m_deadstemc_to_fire(p)*dt
+ endif
+
veg_cs%deadstemc(p) = veg_cs%deadstemc(p) - veg_cf%m_deadstemc_to_litter_fire(p) * dt
+ if (p == 997671) then
+ write(iulog,*) '8.', &
+ veg_cs%deadstemc(p),veg_cf%m_deadstemc_to_litter_fire(p)*dt
+ endif
veg_cs%livecrootc(p) = veg_cs%livecrootc(p) - veg_cf%m_livecrootc_to_fire(p) * dt
veg_cs%livecrootc(p) = veg_cs%livecrootc(p) - veg_cf%m_livecrootc_to_litter_fire(p)* dt
veg_cs%deadcrootc(p) = veg_cs%deadcrootc(p) - veg_cf%m_deadcrootc_to_fire(p) * dt
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment