Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Changes related soil organic carbon and soil organic flux #22

Merged
merged 2 commits into from
Nov 13, 2024
Merged

Changes related soil organic carbon and soil organic flux #22

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
7f72b39
Fixed overwrit ten changes to in carbon_module and command.f90
fgeter Nov 12, 2024
9d38f63
New total for tot_org
fgeter Nov 12, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
44 changes: 22 additions & 22 deletions src/carbon_module.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -14,9 +14,9 @@ module carbon_module
real :: part_DOC_para = 4000. ! |organic carbon partition coefficient 1000 to 1200 ! 500-2000 !replacing KOC=Liquid-solid partition coefficient for Microbial Biomass (10^3 m3 Mg-1)
real :: hlife_doc_para = 50. !days |DOC half life in groundwater, calculating DOC decay in groundwater ! 0-100
!Allocation of CO2 and Carbon transformation
real :: ABCO2_para_sur = 0.6 ! |Allocation from Microbial Biomass C pool to CO2; 0.6 (surface Litter), 0.85 - 0.68 x (CLAY+SILT) (all other layers) (Parton et al., 1993, 1994)
real :: ABCO2_para_sur = 0.6 ! |Allocation from Microbial Biomass C pool to CO2; 0.6 (surface Litter), 0.85 - 0.68 x (CLAY+SILT) (all other layers) (Parton et al., 1993, 1994)
real :: ABCO2_para_sub = 0. ! |Calculated -Allocation from Microbial Biomass C pool to CO2; 0.6 (surface Litter), 0.85 - 0.68 x (CLAY+SILT) (all other layers) (Parton et al., 1993, 1994)
real :: ABP_para_sur = 0.0 ! |Allocation from Biomass to passive Humus; 0 (surface Litter), 0.003 + 0.032 x SOL_CLAY (all other layers) (Parton et al., 1993, 1994)
real :: ABP_para_sur = 0.0 ! |Allocation from Biomass to passive Humus; 0 (surface Litter), 0.003 + 0.032 x SOL_CLAY (all other layers) (Parton et al., 1993, 1994)
real :: ABP_para_sub = 0.0 ! |Calculated - Allocation from Biomass to passive Humus; 0 (surface Litter), 0.003 + 0.032 x SOL_CLAY (all other layers) (Parton et al., 1993, 1994)
real :: ALMCO2_para_sur = 0.6 ! |Allocation from metabolic Litter to CO2; 0.6 (surface Litter), 0.55 (all other layers) (Parton et al., 1993, 1994)
real :: ALMCO2_para_sub = 0.55 ! |Allocation from metabolic Litter to CO2; 0.6 (surface Litter), 0.55 (all other layers) (Parton et al., 1993, 1994)
Expand All @@ -30,8 +30,8 @@ module carbon_module
!decomposition rates
real :: PRMT_51_para = 1.0 ! |COEF ADJUSTS MICROBIAL ACTIVITY FUNCTION IN TOP SOIL LAYER (0.1_1.),
real :: PRMT_45_para = 0.003 ! |COEF IN CENTURY EQ ALLOCATING SLOW TO PASSIVE HUMUS(0.001_0.05) ORIGINAL VALUE = 0.003, ASP=MAX(.001,PRMT_45-.00009*sol_clay(k,j)), ASP=MAX(.001,PRMT_45+.009*sol_clay(k,j)/100)
real :: BMR_para_sur = 0.0164 ! |Rate of transformation of microbial Biomass and associated products under optimal conditions (surface = 0.0164 d-1; all other layers = 0.02 d-1) (Parton et al., 1993, 1994)
real :: BMR_para_sub = 0.02 ! |Rate of transformation of microbial Biomass and associated products under optimal conditions (surface = 0.0164 d-1; all other layers = 0.02 d-1) (Parton et al., 1993, 1994)
real :: BMR_para_sur = 0.0164 ! |Rate of transformation of microbial Biomass and associated products under optimal conditions (surface = 0.0164 d-1; all other layers = 0.02 d-1) (Parton et al., 1993, 1994)
real :: BMR_para_sub = 0.02 ! |Rate of transformation of microbial Biomass and associated products under optimal conditions (surface = 0.0164 d-1; all other layers = 0.02 d-1) (Parton et al., 1993, 1994)
real :: HPR_para = 0.000012 ! |Rate of transformation of passive Humus under optimal conditions (subsurface layers = 0.000012 d-1) (Parton et al., 1993, 1994)
real :: HSR_para = 0.000548 ! |Rate of transformation of slow Humus under optimal conditions (all layers = 0.0005 d-1) (Parton et al., 1993, 1994; Vitousek et al., 1993)
real :: LMR_para_sur = 0.0405 ! |Rate of transformation of metabolic Litter under optimal conditions (surface = 0.0405 d-1; all other layers = 0.0507 d-1) (Parton et al., 1994)
Expand All @@ -54,7 +54,7 @@ module carbon_module
real :: microb_rate = 0. ! |rate of transformation of microbial biomass and associated products under optimal conditions
real :: meta_rate = 0. ! |rate of transformation of metabolic litter under optimal conditions
real :: str_rate = 0. ! |rate of potential transformation of structural litter under optimal conditions
real :: microb_top_rate = 0. ! |coef adjusts mocribial activity function in top soil layer
real :: microb_top_rate = 0. ! |coef adjusts mocribial activity function in top soil layer
real :: hs_hp = 0. ! |coef in century eq allocating slow to passive humus
real :: microb_koc = 0. !10^3 m^3 Mg-1 |liquid-solid partition coefficient for microbial biomass
real :: min_n_frac = 0. ! |fraction of mineral n sorbed to litter
Expand All @@ -65,24 +65,24 @@ module carbon_module

type organic_allocations
real :: abco2 = 0. ! |Fraction of decomposed microbial biomass allocated to CO2
real :: abl = 0. ! |Fraction of microbial biomass loss due to leaching
real :: abp = 0. ! |Fraction of decomposed microbial biomass allocated to passive humus
real :: almco2 = 0. ! |Fraction of decomposed metabolic litter allocated to CO2
real :: alslco2 = 0. ! |Fraction of decomposed lignin of structural litter allocated to CO2
real :: alslnco2 = 0. ! |Fraction of decomposed lignin of structural litter allocated to CO2
real :: apco2 = 0. ! |Fraction of decomposed passive humus allocated to CO2
real :: asco2 = 0. ! |Fraction of decomposed slow humus allocated to CO2
real :: asp = 0. ! |Fraction of decomposed slow humus allocated to passive
real :: abl = 0. ! |Fraction of microbial biomass loss due to leaching
real :: abp = 0. ! |Fraction of decomposed microbial biomass allocated to passive humus
real :: almco2 = 0. ! |Fraction of decomposed metabolic litter allocated to CO2
real :: alslco2 = 0. ! |Fraction of decomposed lignin of structural litter allocated to CO2
real :: alslnco2 = 0. ! |Fraction of decomposed lignin of structural litter allocated to CO2
real :: apco2 = 0. ! |Fraction of decomposed passive humus allocated to CO2
real :: asco2 = 0. ! |Fraction of decomposed slow humus allocated to CO2
real :: asp = 0. ! |Fraction of decomposed slow humus allocated to passive
end type organic_allocations
type (organic_allocations) :: org_allo
type (organic_allocations) :: org_alloz

type organic_controls
real :: cdg = 0. ! |soil temperature control on biological processes
real :: cs = 0. ! |combined factor controlling biological processes
real :: ox = 0. ! |oxygen control on biological processes
real :: sut = 0. ! |soil water control on biological processes
real :: x1 = 0. ! |tillage control on residue decomposition
real :: cdg = 0. ! |soil temperature control on biological processes
real :: cs = 0. ! |combined factor controlling biological processes
real :: ox = 0. ! |oxygen control on biological processes
real :: sut = 0. ! |soil water control on biological processes
real :: x1 = 0. ! |tillage control on residue decomposition
real :: xbmt = 0. ! |control on transformation of microbial biomass by soil texture and structure
real :: xlslf = 0. ! |control on potential transformation of structural litter by lignin fraction
end type organic_controls
Expand All @@ -101,15 +101,15 @@ module carbon_module
real :: cnr = 0. ! |c/n ratio of standing dead
real :: ncbm = 0. ! |n/c ratio of biomass
real :: nchp = 0. ! |n/c ratio of passive humus
real :: nchs = 0. ! |n/c ration of slow humus
real :: nchs = 0. ! |n/c ration of slow humus
end type organic_ratio
type (organic_ratio) :: org_ratio

type organic_transformations
real :: bmctp = 0. !kg ha-1 day-1 |potential transformation of C in microbial biomass
real :: bmntp = 0. !kg ha-1 day-1 |potential transformation of N in microbial biomass
real :: hsctp = 0. !kg ha-1 day-1 |potential transformation of C in slow humus
real :: hsntp = 0. !kg ha-1 day-1 |potential transformation of N in slow humus
real :: bmntp = 0. !kg ha-1 day-1 |potential transformation of N in microbial biomass
real :: hsctp = 0. !kg ha-1 day-1 |potential transformation of C in slow humus
real :: hsntp = 0. !kg ha-1 day-1 |potential transformation of N in slow humus
real :: hpctp = 0. !kg ha-1 day-1 |potential transformation of C in passive humus
real :: hpntp = 0. !kg ha-1 day-1 |potential transformation of N in passive humus
real :: lmctp = 0. !kg ha-1 day-1 |potential transformation of C in metabolic litter
Expand Down
13 changes: 10 additions & 3 deletions src/command.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -112,6 +112,11 @@ subroutine command
if (cs_db%num_tot > 0 .and. obcs_alloc(icmd).eq.1) then
obcs(icmd)%hin_sur(1) = obcs(icmd)%hin_sur(1) + frac_in * obcs(iob)%hd(3)
end if
! add to tile flow
ob(icmd)%hin_til = ob(icmd)%hin_til + frac_in * ob(iob)%hd(5)
if (cs_db%num_tot > 0 .and. obcs_alloc(icmd).eq.1) then
obcs(icmd)%hin_til(1) = obcs(icmd)%hin_til(1) + frac_in * obcs(iob)%hd(5)
end if
! add to lateral soil runon
ob(icmd)%hin_lat = ob(icmd)%hin_lat + frac_in * ob(iob)%hd(4)
if (cs_db%num_tot > 0 .and. obcs_alloc(icmd).eq.1) then
Expand Down Expand Up @@ -465,9 +470,11 @@ subroutine command
!if ((time%yrc == 2007 .AND. time%day == 213) .OR. (time%yrc == 2010 .AND. time%day == 319) &
! .OR.(time%yrc == 2011 .AND. time%day == 324)) then
if (ihru == 1) then
do nly = 1, soil(ihru)%nly
!soil1(ihru)%tot(nly)%c = soil1(ihru)%hact(nly)%c + soil1(ihru)%hsta(nly)%c + soil1(ihru)%microb(nly)%c
end do
if (bsn_cc%cswat /= 2) then
do nly = 1, soil(ihru)%nly
soil1(ihru)%tot(nly)%c = soil1(ihru)%hact(nly)%c + soil1(ihru)%hsta(nly)%c + soil1(ihru)%microb(nly)%c
end do
end if
write (9999,*) time%day, time%mo, time%day_mo, time%yrc, ob(ihru)%typ, ob(ihru)%name, &
(soil1(ihru)%tot(ly)%c/1000, ly = 1, soil(ihru)%nly)
end if
Expand Down
14 changes: 13 additions & 1 deletion src/soil_nutcarb_write.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -53,6 +53,15 @@ subroutine soil_nutcarb_write(out_freq)
soil1(j)%tot_org = soil_org_z
soil_prof_hact = soil_org_z
soil_prof_hsta = soil_org_z
soil_prof_hsta = soil_org_z
soil_prof_str = soil_org_z
soil_prof_lig = soil_org_z
soil_prof_meta = soil_org_z
soil_prof_man = soil_org_z
soil_prof_hs = soil_org_z
soil_prof_hp = soil_org_z
soil_prof_microb = soil_org_z
soil_prof_water = soil_org_z
do ly = 1, soil(j)%nly
soil_prof_hact = soil_prof_hact + soil1(j)%hact(ly)
soil_prof_hsta = soil_prof_hsta + soil1(j)%hsta(ly)
Expand All @@ -65,7 +74,10 @@ subroutine soil_nutcarb_write(out_freq)
soil_prof_microb = soil_prof_microb + soil1(j)%microb(ly)
soil_prof_water = soil_prof_water + soil1(j)%water(ly)
end do
soil1(j)%tot_org = soil_prof_hact + soil_prof_hsta + soil_prof_microb
! This a test of push
! soil1(j)%tot_org = soil_prof_hact + soil_prof_hsta + soil_prof_microb
soil1(j)%tot_org = soil_prof_hs + soil_prof_hp + soil_prof_microb + soil_prof_meta + &
soil_prof_str + soil_prof_lig

!write all organic carbon for the plant community
write (4560,*) time%day, time%mo, time%day_mo, time%yrc, j, ob(iob)%gis_id, ob(iob)%name, &
Expand Down