Created
May 21, 2024 14:28
-
-
Save bishtgautam/8782af4f41a5e6d2924659dd807a9f81 to your computer and use it in GitHub Desktop.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
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