Skip to content

Commit 370c1a4

Browse files
arnoldjjmsarnoldjjms
authored
Merge pull request #22 from fgeter/main
Changes related soil organic carbon and soil organic flux
2 parents f9b1200 + 9d38f63 commit 370c1a4

File tree

3 files changed

+45
-26
lines changed

3 files changed

+45
-26
lines changed

‎src/carbon_module.f90‎

Lines changed: 22 additions & 22 deletions
Original file line numberDiff line numberDiff line change
@@ -14,9 +14,9 @@ module carbon_module
1414
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)
1515
real :: hlife_doc_para = 50. !days |DOC half life in groundwater, calculating DOC decay in groundwater ! 0-100
1616
!Allocation of CO2 and Carbon transformation
17-
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)
17+
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)
1818
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)
19-
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)
19+
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)
2020
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)
2121
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)
2222
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)
@@ -30,8 +30,8 @@ module carbon_module
3030
!decomposition rates
3131
real :: PRMT_51_para = 1.0 ! |COEF ADJUSTS MICROBIAL ACTIVITY FUNCTION IN TOP SOIL LAYER (0.1_1.),
3232
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)
33-
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)
34-
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)
33+
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)
34+
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)
3535
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)
3636
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)
3737
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)
@@ -54,7 +54,7 @@ module carbon_module
5454
real :: microb_rate = 0. ! |rate of transformation of microbial biomass and associated products under optimal conditions
5555
real :: meta_rate = 0. ! |rate of transformation of metabolic litter under optimal conditions
5656
real :: str_rate = 0. ! |rate of potential transformation of structural litter under optimal conditions
57-
real :: microb_top_rate = 0. ! |coef adjusts mocribial activity function in top soil layer
57+
real :: microb_top_rate = 0. ! |coef adjusts mocribial activity function in top soil layer
5858
real :: hs_hp = 0. ! |coef in century eq allocating slow to passive humus
5959
real :: microb_koc = 0. !10^3 m^3 Mg-1 |liquid-solid partition coefficient for microbial biomass
6060
real :: min_n_frac = 0. ! |fraction of mineral n sorbed to litter
@@ -65,24 +65,24 @@ module carbon_module
6565

6666
type organic_allocations
6767
real :: abco2 = 0. ! |Fraction of decomposed microbial biomass allocated to CO2
68-
real :: abl = 0. ! |Fraction of microbial biomass loss due to leaching
69-
real :: abp = 0. ! |Fraction of decomposed microbial biomass allocated to passive humus
70-
real :: almco2 = 0. ! |Fraction of decomposed metabolic litter allocated to CO2
71-
real :: alslco2 = 0. ! |Fraction of decomposed lignin of structural litter allocated to CO2
72-
real :: alslnco2 = 0. ! |Fraction of decomposed lignin of structural litter allocated to CO2
73-
real :: apco2 = 0. ! |Fraction of decomposed passive humus allocated to CO2
74-
real :: asco2 = 0. ! |Fraction of decomposed slow humus allocated to CO2
75-
real :: asp = 0. ! |Fraction of decomposed slow humus allocated to passive
68+
real :: abl = 0. ! |Fraction of microbial biomass loss due to leaching
69+
real :: abp = 0. ! |Fraction of decomposed microbial biomass allocated to passive humus
70+
real :: almco2 = 0. ! |Fraction of decomposed metabolic litter allocated to CO2
71+
real :: alslco2 = 0. ! |Fraction of decomposed lignin of structural litter allocated to CO2
72+
real :: alslnco2 = 0. ! |Fraction of decomposed lignin of structural litter allocated to CO2
73+
real :: apco2 = 0. ! |Fraction of decomposed passive humus allocated to CO2
74+
real :: asco2 = 0. ! |Fraction of decomposed slow humus allocated to CO2
75+
real :: asp = 0. ! |Fraction of decomposed slow humus allocated to passive
7676
end type organic_allocations
7777
type (organic_allocations) :: org_allo
7878
type (organic_allocations) :: org_alloz
7979

8080
type organic_controls
81-
real :: cdg = 0. ! |soil temperature control on biological processes
82-
real :: cs = 0. ! |combined factor controlling biological processes
83-
real :: ox = 0. ! |oxygen control on biological processes
84-
real :: sut = 0. ! |soil water control on biological processes
85-
real :: x1 = 0. ! |tillage control on residue decomposition
81+
real :: cdg = 0. ! |soil temperature control on biological processes
82+
real :: cs = 0. ! |combined factor controlling biological processes
83+
real :: ox = 0. ! |oxygen control on biological processes
84+
real :: sut = 0. ! |soil water control on biological processes
85+
real :: x1 = 0. ! |tillage control on residue decomposition
8686
real :: xbmt = 0. ! |control on transformation of microbial biomass by soil texture and structure
8787
real :: xlslf = 0. ! |control on potential transformation of structural litter by lignin fraction
8888
end type organic_controls
@@ -101,15 +101,15 @@ module carbon_module
101101
real :: cnr = 0. ! |c/n ratio of standing dead
102102
real :: ncbm = 0. ! |n/c ratio of biomass
103103
real :: nchp = 0. ! |n/c ratio of passive humus
104-
real :: nchs = 0. ! |n/c ration of slow humus
104+
real :: nchs = 0. ! |n/c ration of slow humus
105105
end type organic_ratio
106106
type (organic_ratio) :: org_ratio
107107

108108
type organic_transformations
109109
real :: bmctp = 0. !kg ha-1 day-1 |potential transformation of C in microbial biomass
110-
real :: bmntp = 0. !kg ha-1 day-1 |potential transformation of N in microbial biomass
111-
real :: hsctp = 0. !kg ha-1 day-1 |potential transformation of C in slow humus
112-
real :: hsntp = 0. !kg ha-1 day-1 |potential transformation of N in slow humus
110+
real :: bmntp = 0. !kg ha-1 day-1 |potential transformation of N in microbial biomass
111+
real :: hsctp = 0. !kg ha-1 day-1 |potential transformation of C in slow humus
112+
real :: hsntp = 0. !kg ha-1 day-1 |potential transformation of N in slow humus
113113
real :: hpctp = 0. !kg ha-1 day-1 |potential transformation of C in passive humus
114114
real :: hpntp = 0. !kg ha-1 day-1 |potential transformation of N in passive humus
115115
real :: lmctp = 0. !kg ha-1 day-1 |potential transformation of C in metabolic litter

‎src/command.f90‎

Lines changed: 10 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -112,6 +112,11 @@ subroutine command
112112
if (cs_db%num_tot > 0 .and. obcs_alloc(icmd).eq.1) then
113113
obcs(icmd)%hin_sur(1) = obcs(icmd)%hin_sur(1) + frac_in * obcs(iob)%hd(3)
114114
end if
115+
! add to tile flow
116+
ob(icmd)%hin_til = ob(icmd)%hin_til + frac_in * ob(iob)%hd(5)
117+
if (cs_db%num_tot > 0 .and. obcs_alloc(icmd).eq.1) then
118+
obcs(icmd)%hin_til(1) = obcs(icmd)%hin_til(1) + frac_in * obcs(iob)%hd(5)
119+
end if
115120
! add to lateral soil runon
116121
ob(icmd)%hin_lat = ob(icmd)%hin_lat + frac_in * ob(iob)%hd(4)
117122
if (cs_db%num_tot > 0 .and. obcs_alloc(icmd).eq.1) then
@@ -465,9 +470,11 @@ subroutine command
465470
!if ((time%yrc == 2007 .AND. time%day == 213) .OR. (time%yrc == 2010 .AND. time%day == 319) &
466471
! .OR.(time%yrc == 2011 .AND. time%day == 324)) then
467472
if (ihru == 1) then
468-
do nly = 1, soil(ihru)%nly
469-
!soil1(ihru)%tot(nly)%c = soil1(ihru)%hact(nly)%c + soil1(ihru)%hsta(nly)%c + soil1(ihru)%microb(nly)%c
470-
end do
473+
if (bsn_cc%cswat /= 2) then
474+
do nly = 1, soil(ihru)%nly
475+
soil1(ihru)%tot(nly)%c = soil1(ihru)%hact(nly)%c + soil1(ihru)%hsta(nly)%c + soil1(ihru)%microb(nly)%c
476+
end do
477+
end if
471478
write (9999,*) time%day, time%mo, time%day_mo, time%yrc, ob(ihru)%typ, ob(ihru)%name, &
472479
(soil1(ihru)%tot(ly)%c/1000, ly = 1, soil(ihru)%nly)
473480
end if

‎src/soil_nutcarb_write.f90‎

Lines changed: 13 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -53,6 +53,15 @@ subroutine soil_nutcarb_write(out_freq)
5353
soil1(j)%tot_org = soil_org_z
5454
soil_prof_hact = soil_org_z
5555
soil_prof_hsta = soil_org_z
56+
soil_prof_hsta = soil_org_z
57+
soil_prof_str = soil_org_z
58+
soil_prof_lig = soil_org_z
59+
soil_prof_meta = soil_org_z
60+
soil_prof_man = soil_org_z
61+
soil_prof_hs = soil_org_z
62+
soil_prof_hp = soil_org_z
63+
soil_prof_microb = soil_org_z
64+
soil_prof_water = soil_org_z
5665
do ly = 1, soil(j)%nly
5766
soil_prof_hact = soil_prof_hact + soil1(j)%hact(ly)
5867
soil_prof_hsta = soil_prof_hsta + soil1(j)%hsta(ly)
@@ -65,7 +74,10 @@ subroutine soil_nutcarb_write(out_freq)
6574
soil_prof_microb = soil_prof_microb + soil1(j)%microb(ly)
6675
soil_prof_water = soil_prof_water + soil1(j)%water(ly)
6776
end do
68-
soil1(j)%tot_org = soil_prof_hact + soil_prof_hsta + soil_prof_microb
77+
! This a test of push
78+
! soil1(j)%tot_org = soil_prof_hact + soil_prof_hsta + soil_prof_microb
79+
soil1(j)%tot_org = soil_prof_hs + soil_prof_hp + soil_prof_microb + soil_prof_meta + &
80+
soil_prof_str + soil_prof_lig
6981

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

0 commit comments

Comments
 (0)