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

MIMICS-related refactor before introducing MIMICS #1340

Merged
merged 20 commits into from
Jun 3, 2021
Merged

MIMICS-related refactor before introducing MIMICS #1340

Show file tree
Hide file tree
Changes from 3 commits
Commits
Show all changes
20 commits
Select commit Hold shift + click to select a range
dfe3f71
Derive instead of hardwiring i_litr2, i_litr3, i_cwd
slevis-lmwg Apr 13, 2021
30c8b44
Reduce repetition of decomp_k lines in *CascadeBGC and *CascadeCN
slevis-lmwg Apr 14, 2021
0f60580
Moved initial_Cstock vars from namelist to new params file
slevis-lmwg Apr 18, 2021
22f3a3a
Rm sub setup_logic_century_soilbgcdecompcascade from CLMBuildNamelist.pm
slevis-lmwg Apr 18, 2021
b62b89d
Replacing explicit litter-pool lines of code w do-loops
slevis-lmwg Apr 28, 2021
2224191
Replacing explicit litter-pool lines of code w do-loops (contd)
slevis-lmwg Apr 29, 2021
55cf695
Replacing explicit litter-pool lines of code w do-loops (3rd commit)
slevis-lmwg Apr 29, 2021
8417cd8
Replacing explicit litter-pool lines of code w do-loops (4th commit)
slevis-lmwg Apr 30, 2021
69f4c5c
Rename i_soil* indices to descriptive names
slevis-lmwg Apr 30, 2021
1fed563
Merge tag 'ctsm5.1.dev040' into premimics_refactor
slevis-lmwg May 26, 2021
3122146
Revisions (part 1) in response to review
slevis-lmwg May 28, 2021
a63530d
Revisions (part 2) in response to review
slevis-lmwg May 28, 2021
1e49b3a
Merge tag 'ctsm5.1.dev042' into premimics_refactor
slevis-lmwg May 28, 2021
c5d0ee4
Updating to latest params files and including draft ChangeLog
slevis-lmwg May 28, 2021
8b9b418
Small update to ChangeLog
slevis-lmwg Jun 1, 2021
b998307
Another update to the ChangeLog
slevis-lmwg Jun 1, 2021
75a8e39
Add expected fails for the CN tests
ekluzek Jun 3, 2021
70bf79c
Update change files
ekluzek Jun 3, 2021
1fde421
Small note about memory increase in change file
ekluzek Jun 3, 2021
a1186e3
Note about new issue
ekluzek Jun 3, 2021
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
4 changes: 2 additions & 2 deletions bld/CLMBuildNamelist.pm
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3594,10 +3594,10 @@ sub setup_logic_century_soilbgcdecompcascade {

if ( (&value_is_true($nl->get_value('use_cn')) || &value_is_true($nl->get_value('use_fates'))) &&
&value_is_true($nl->get_value('use_century_decomp')) ) {
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'initial_Cstocks',
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl,
'use_cn' => $nl->get_value('use_cn'), 'use_fates' => $nl->get_value('use_fates'),
'use_century_decomp' => $nl->get_value('use_century_decomp') );
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'initial_Cstocks_depth',
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl,
slevis-lmwg marked this conversation as resolved.
Show resolved Hide resolved
'use_cn' => $nl->get_value('use_cn'), 'use_fates' => $nl->get_value('use_fates'),
'use_century_decomp' => $nl->get_value('use_century_decomp') );
}
Expand Down
14 changes: 0 additions & 14 deletions bld/namelist_files/namelist_defaults_ctsm.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1416,20 +1416,6 @@ lnd/clm2/surfdata_map/release-clm5.0.30/surfdata_ne0np4.CONUS.ne30x8_hist_78pfts
<leaf_mr_vcm phys="clm5_0" >0.015d00</leaf_mr_vcm>
<leaf_mr_vcm phys="clm4_5" >0.015d00</leaf_mr_vcm>

<!-- Initial Carbon stocks for BGC -->
<initial_Cstocks phys="clm4_5" use_cn=".true." use_fates=".false." use_century_decomp=".true." >20.0d00, 20.0d00, 20.0d00</initial_Cstocks>
ekluzek marked this conversation as resolved.
Show resolved Hide resolved
<initial_Cstocks phys="clm5_0" use_cn=".true." use_fates=".false." use_century_decomp=".true." >200.0d00, 200.0d00, 200.0d00</initial_Cstocks>
<initial_Cstocks phys="clm5_1" use_cn=".true." use_fates=".false." use_century_decomp=".true." >200.0d00, 200.0d00, 200.0d00</initial_Cstocks>
<initial_Cstocks phys="clm4_5" use_cn=".false." use_fates=".true." use_century_decomp=".true." >20.0d00, 20.0d00, 20.0d00</initial_Cstocks>
<initial_Cstocks phys="clm5_0" use_cn=".false." use_fates=".true." use_century_decomp=".true." >200.0d00, 200.0d00, 200.0d00</initial_Cstocks>
<initial_Cstocks phys="clm5_1" use_cn=".false." use_fates=".true." use_century_decomp=".true." >200.0d00, 200.0d00, 200.0d00</initial_Cstocks>

<initial_Cstocks_depth phys="clm5_1" use_cn=".true." use_fates=".false." use_century_decomp=".true." >1.50d00</initial_Cstocks_depth>
<initial_Cstocks_depth phys="clm5_0" use_cn=".true." use_fates=".false." use_century_decomp=".true." >1.50d00</initial_Cstocks_depth>
<initial_Cstocks_depth phys="clm4_5" use_cn=".true." use_fates=".false." use_century_decomp=".true." >0.3</initial_Cstocks_depth>
ekluzek marked this conversation as resolved.
Show resolved Hide resolved
<initial_Cstocks_depth phys="clm5_1" use_cn=".false." use_fates=".true." use_century_decomp=".true." >1.50d00</initial_Cstocks_depth>
<initial_Cstocks_depth phys="clm5_0" use_cn=".false." use_fates=".true." use_century_decomp=".true." >1.50d00</initial_Cstocks_depth>
<initial_Cstocks_depth phys="clm4_5" use_cn=".false." use_fates=".true." use_century_decomp=".true." >0.3</initial_Cstocks_depth>

<initial_vegC phys="clm5_1" use_cn=".true." mm_nuptake_opt=".true." >100.d00</initial_vegC>
<initial_vegC phys="clm5_0" use_cn=".true." mm_nuptake_opt=".true." >100.d00</initial_vegC>
Expand Down
10 changes: 0 additions & 10 deletions bld/namelist_files/namelist_definition_ctsm.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -306,16 +306,6 @@ Critical threshold of negative Nitrogen to die (abort when Nitrogen states are b
Critical threshold of negative Carbon to die (abort when Carbon states are below this value)
</entry>

<entry id="initial_Cstocks" type="real(3)" category="clm_physics"
group="century_soilbgcdecompcascade" >
Initial stocks of Carbon to use in soil organic matter pools for CENTURY decomposition
</entry>

<entry id="initial_Cstocks_depth" type="real(3)" category="clm_physics"
group="century_soilbgcdecompcascade" >
Soil depth to place initial stocks of Carbon in soil organic matter pools for CENTURY decomposition
</entry>

<entry id="freelivfix_slope_wET" type="real" category="clm_physics"
group="mineral_nitrogen_dynamics" >
Slope of free living Nitrogen fixation with annual ET
Expand Down
2 changes: 1 addition & 1 deletion src/biogeochem/CNC14DecayMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@ module CNC14DecayMod
! !USES:
use shr_kind_mod , only : r8 => shr_kind_r8
use clm_time_manager , only : get_step_size_real, get_days_per_year
use clm_varpar , only : ndecomp_cascade_transitions, nlevdecomp, ndecomp_pools
use clm_varpar , only : nlevdecomp, ndecomp_pools
ekluzek marked this conversation as resolved.
Show resolved Hide resolved
use clm_varcon , only : secspday
use clm_varctl , only : spinup_state
use decompMod , only : bounds_type
Expand Down
18 changes: 9 additions & 9 deletions src/biogeochem/CNCStateUpdate1Mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@ module CNCStateUpdate1Mod
use shr_log_mod , only : errMsg => shr_log_errMsg
use clm_varpar , only : ndecomp_cascade_transitions, nlevdecomp
use clm_time_manager , only : get_step_size_real
use clm_varpar , only : i_met_lit, i_cel_lit, i_lig_lit, i_cwd
use clm_varpar , only : i_met_lit, i_litr2, i_litr3, i_cwd
slevis-lmwg marked this conversation as resolved.
Show resolved Hide resolved
use pftconMod , only : npcropmin, nc3crop, pftcon
use abortutils , only : endrun
use decompMod , only : bounds_type
Expand Down Expand Up @@ -71,9 +71,9 @@ subroutine CStateUpdateDynPatch(bounds, num_soilc_with_inactive, filter_soilc_wi
c = filter_soilc_with_inactive(fc)
cs_soil%decomp_cpools_vr_col(c,j,i_met_lit) = cs_soil%decomp_cpools_vr_col(c,j,i_met_lit) + &
cf_veg%dwt_frootc_to_litr_met_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_cel_lit) = cs_soil%decomp_cpools_vr_col(c,j,i_cel_lit) + &
cs_soil%decomp_cpools_vr_col(c,j,i_litr2) = cs_soil%decomp_cpools_vr_col(c,j,i_litr2) + &
cf_veg%dwt_frootc_to_litr_cel_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_lig_lit) = cs_soil%decomp_cpools_vr_col(c,j,i_lig_lit) + &
cs_soil%decomp_cpools_vr_col(c,j,i_litr3) = cs_soil%decomp_cpools_vr_col(c,j,i_litr3) + &
cf_veg%dwt_frootc_to_litr_lig_c_col(c,j) * dt
ekluzek marked this conversation as resolved.
Show resolved Hide resolved
cs_soil%decomp_cpools_vr_col(c,j,i_cwd) = cs_soil%decomp_cpools_vr_col(c,j,i_cwd) + &
( cf_veg%dwt_livecrootc_to_cwdc_col(c,j) + cf_veg%dwt_deadcrootc_to_cwdc_col(c,j) ) * dt
ekluzek marked this conversation as resolved.
Show resolved Hide resolved
Expand All @@ -91,8 +91,8 @@ subroutine CStateUpdateDynPatch(bounds, num_soilc_with_inactive, filter_soilc_wi
! CStateUpdate1) if use_fates is true? Specifically, some portion or all of the fluxes
! from these updates in CStateUpdate1:
! cf_soil%decomp_cpools_sourcesink_col(c,j,i_met_lit) = cf_soil%FATES_c_to_litr_lab_c_col(c,j) * dt
! cf_soil%decomp_cpools_sourcesink_col(c,j,i_cel_lit) = cf_soil%FATES_c_to_litr_cel_c_col(c,j) * dt
! cf_soil%decomp_cpools_sourcesink_col(c,j,i_lig_lit) = cf_soil%FATES_c_to_litr_lig_c_col(c,j) * dt
! cf_soil%decomp_cpools_sourcesink_col(c,j,i_litr2) = cf_soil%FATES_c_to_litr_cel_c_col(c,j) * dt
! cf_soil%decomp_cpools_sourcesink_col(c,j,i_litr3) = cf_soil%FATES_c_to_litr_lig_c_col(c,j) * dt

end associate

Expand Down Expand Up @@ -197,9 +197,9 @@ subroutine CStateUpdate1( num_soilc, filter_soilc, num_soilp, filter_soilp, &
! phenology and dynamic land cover fluxes
cf_soil%decomp_cpools_sourcesink_col(c,j,i_met_lit) = &
cf_veg%phenology_c_to_litr_met_c_col(c,j) *dt
cf_soil%decomp_cpools_sourcesink_col(c,j,i_cel_lit) = &
cf_soil%decomp_cpools_sourcesink_col(c,j,i_litr2) = &
cf_veg%phenology_c_to_litr_cel_c_col(c,j) *dt
cf_soil%decomp_cpools_sourcesink_col(c,j,i_lig_lit) = &
cf_soil%decomp_cpools_sourcesink_col(c,j,i_litr3) = &
cf_veg%phenology_c_to_litr_lig_c_col(c,j) *dt

! NOTE(wjs, 2017-01-02) This used to be set to a non-zero value, but the
Expand All @@ -216,8 +216,8 @@ subroutine CStateUpdate1( num_soilc, filter_soilc, num_soilp, filter_soilp, &
! TODO(wjs, 2017-01-02) Should some portion or all of the following fluxes
! be moved to the updates in CStateUpdateDynPatch?
cf_soil%decomp_cpools_sourcesink_col(c,j,i_met_lit) = cf_soil%FATES_c_to_litr_lab_c_col(c,j) * dt
cf_soil%decomp_cpools_sourcesink_col(c,j,i_cel_lit) = cf_soil%FATES_c_to_litr_cel_c_col(c,j) * dt
cf_soil%decomp_cpools_sourcesink_col(c,j,i_lig_lit) = cf_soil%FATES_c_to_litr_lig_c_col(c,j) * dt
cf_soil%decomp_cpools_sourcesink_col(c,j,i_litr2) = cf_soil%FATES_c_to_litr_cel_c_col(c,j) * dt
cf_soil%decomp_cpools_sourcesink_col(c,j,i_litr3) = cf_soil%FATES_c_to_litr_lig_c_col(c,j) * dt
end do
end do
endif
Expand Down
18 changes: 9 additions & 9 deletions src/biogeochem/CNCStateUpdate2Mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -9,7 +9,7 @@ module CNCStateUpdate2Mod
use shr_log_mod , only : errMsg => shr_log_errMsg
use abortutils , only : endrun
use clm_time_manager , only : get_step_size_real
use clm_varpar , only : nlevdecomp, i_met_lit, i_cel_lit, i_lig_lit, i_cwd
use clm_varpar , only : nlevdecomp, i_met_lit, i_litr2, i_litr3, i_cwd
use CNvegCarbonStateType , only : cnveg_carbonstate_type
use CNVegCarbonFluxType , only : cnveg_carbonflux_type
use SoilBiogeochemCarbonStatetype , only : soilbiogeochem_carbonstate_type
Expand Down Expand Up @@ -66,10 +66,10 @@ subroutine CStateUpdate2(num_soilc, filter_soilc, num_soilp, filter_soilp, &
! column gap mortality fluxes
cs_soil%decomp_cpools_vr_col(c,j,i_met_lit) = &
cs_soil%decomp_cpools_vr_col(c,j,i_met_lit) + cf_veg%gap_mortality_c_to_litr_met_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_cel_lit) = &
cs_soil%decomp_cpools_vr_col(c,j,i_cel_lit) + cf_veg%gap_mortality_c_to_litr_cel_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_lig_lit) = &
cs_soil%decomp_cpools_vr_col(c,j,i_lig_lit) + cf_veg%gap_mortality_c_to_litr_lig_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_litr2) = &
cs_soil%decomp_cpools_vr_col(c,j,i_litr2) + cf_veg%gap_mortality_c_to_litr_cel_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_litr3) = &
cs_soil%decomp_cpools_vr_col(c,j,i_litr3) + cf_veg%gap_mortality_c_to_litr_lig_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_cwd) = &
cs_soil%decomp_cpools_vr_col(c,j,i_cwd) + cf_veg%gap_mortality_c_to_cwdc_col(c,j) * dt

Expand Down Expand Up @@ -172,10 +172,10 @@ subroutine CStateUpdate2h(num_soilc, filter_soilc, num_soilp, filter_soilp, &
! column harvest fluxes
cs_soil%decomp_cpools_vr_col(c,j,i_met_lit) = &
cs_soil%decomp_cpools_vr_col(c,j,i_met_lit) + cf_veg%harvest_c_to_litr_met_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_cel_lit) = &
cs_soil%decomp_cpools_vr_col(c,j,i_cel_lit) + cf_veg%harvest_c_to_litr_cel_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_lig_lit) = &
cs_soil%decomp_cpools_vr_col(c,j,i_lig_lit) + cf_veg%harvest_c_to_litr_lig_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_litr2) = &
cs_soil%decomp_cpools_vr_col(c,j,i_litr2) + cf_veg%harvest_c_to_litr_cel_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_litr3) = &
cs_soil%decomp_cpools_vr_col(c,j,i_litr3) + cf_veg%harvest_c_to_litr_lig_c_col(c,j) * dt
cs_soil%decomp_cpools_vr_col(c,j,i_cwd) = &
cs_soil%decomp_cpools_vr_col(c,j,i_cwd) + cf_veg%harvest_c_to_cwdc_col(c,j) * dt

Expand Down
6 changes: 3 additions & 3 deletions src/biogeochem/CNCStateUpdate3Mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -9,7 +9,7 @@ module CNCStateUpdate3Mod
use shr_log_mod , only : errMsg => shr_log_errMsg
use abortutils , only : endrun
use clm_time_manager , only : get_step_size_real
use clm_varpar , only : nlevdecomp, ndecomp_pools, i_cwd, i_met_lit, i_cel_lit, i_lig_lit
use clm_varpar , only : nlevdecomp, ndecomp_pools, i_cwd, i_met_lit, i_litr2, i_litr3
use CNVegCarbonStateType , only : cnveg_carbonstate_type
use CNVegCarbonFluxType , only : cnveg_carbonflux_type
use SoilBiogeochemCarbonStateType , only : soilbiogeochem_carbonstate_type
Expand Down Expand Up @@ -66,9 +66,9 @@ subroutine CStateUpdate3( num_soilc, filter_soilc, num_soilp, filter_soilp, &
! patch-level wood to column-level litter (uncombusted wood)
cs_soil%decomp_cpools_vr_col(c,j,i_met_lit) = cs_soil%decomp_cpools_vr_col(c,j,i_met_lit) + &
cf_veg%m_c_to_litr_met_fire_col(c,j)* dt
cs_soil%decomp_cpools_vr_col(c,j,i_cel_lit) = cs_soil%decomp_cpools_vr_col(c,j,i_cel_lit) + &
cs_soil%decomp_cpools_vr_col(c,j,i_litr2) = cs_soil%decomp_cpools_vr_col(c,j,i_litr2) + &
cf_veg%m_c_to_litr_cel_fire_col(c,j)* dt
cs_soil%decomp_cpools_vr_col(c,j,i_lig_lit) = cs_soil%decomp_cpools_vr_col(c,j,i_lig_lit) + &
cs_soil%decomp_cpools_vr_col(c,j,i_litr3) = cs_soil%decomp_cpools_vr_col(c,j,i_litr3) + &
cf_veg%m_c_to_litr_lig_fire_col(c,j)* dt
end do
end do
Expand Down
10 changes: 5 additions & 5 deletions src/biogeochem/CNNStateUpdate1Mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@ module CNNStateUpdate1Mod
use shr_kind_mod , only: r8 => shr_kind_r8
use clm_time_manager , only : get_step_size_real
use clm_varpar , only : nlevdecomp, ndecomp_pools, ndecomp_cascade_transitions
use clm_varpar , only : i_met_lit, i_cel_lit, i_lig_lit, i_cwd
use clm_varpar , only : i_met_lit, i_litr2, i_litr3, i_cwd
use clm_varctl , only : iulog, use_nitrif_denitrif
use clm_varcon , only : nitrif_n2o_loss_frac
use pftconMod , only : npcropmin, pftcon
Expand Down Expand Up @@ -67,9 +67,9 @@ subroutine NStateUpdateDynPatch(bounds, num_soilc_with_inactive, filter_soilc_wi
c = filter_soilc_with_inactive(fc)
ns_soil%decomp_npools_vr_col(c,j,i_met_lit) = ns_soil%decomp_npools_vr_col(c,j,i_met_lit) + &
nf_veg%dwt_frootn_to_litr_met_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_cel_lit) = ns_soil%decomp_npools_vr_col(c,j,i_cel_lit) + &
ns_soil%decomp_npools_vr_col(c,j,i_litr2) = ns_soil%decomp_npools_vr_col(c,j,i_litr2) + &
nf_veg%dwt_frootn_to_litr_cel_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_lig_lit) = ns_soil%decomp_npools_vr_col(c,j,i_lig_lit) + &
ns_soil%decomp_npools_vr_col(c,j,i_litr3) = ns_soil%decomp_npools_vr_col(c,j,i_litr3) + &
nf_veg%dwt_frootn_to_litr_lig_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_cwd) = ns_soil%decomp_npools_vr_col(c,j,i_cwd) + &
( nf_veg%dwt_livecrootn_to_cwdn_col(c,j) + nf_veg%dwt_deadcrootn_to_cwdn_col(c,j) ) * dt
Expand Down Expand Up @@ -132,10 +132,10 @@ subroutine NStateUpdate1(num_soilc, filter_soilc, num_soilp, filter_soilp, &
nf_soil%decomp_npools_sourcesink_col(c,j,i_met_lit) = &
nf_veg%phenology_n_to_litr_met_n_col(c,j) * dt

nf_soil%decomp_npools_sourcesink_col(c,j,i_cel_lit) = &
nf_soil%decomp_npools_sourcesink_col(c,j,i_litr2) = &
nf_veg%phenology_n_to_litr_cel_n_col(c,j) * dt

nf_soil%decomp_npools_sourcesink_col(c,j,i_lig_lit) = &
nf_soil%decomp_npools_sourcesink_col(c,j,i_litr3) = &
nf_veg%phenology_n_to_litr_lig_n_col(c,j) * dt

! NOTE(wjs, 2017-01-02) This used to be set to a non-zero value, but the
Expand Down
18 changes: 9 additions & 9 deletions src/biogeochem/CNNStateUpdate2Mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@ module CNNStateUpdate2Mod
use shr_kind_mod , only : r8 => shr_kind_r8
use clm_time_manager , only : get_step_size_real
use clm_varpar , only : nlevsoi, nlevdecomp
use clm_varpar , only : i_met_lit, i_cel_lit, i_lig_lit, i_cwd
use clm_varpar , only : i_met_lit, i_litr2, i_litr3, i_cwd
use clm_varctl , only : iulog
use CNVegNitrogenStateType , only : cnveg_nitrogenstate_type
use CNVegNitrogenFluxType , only : cnveg_nitrogenflux_type
Expand Down Expand Up @@ -66,10 +66,10 @@ subroutine NStateUpdate2(num_soilc, filter_soilc, num_soilp, filter_soilp, &

ns_soil%decomp_npools_vr_col(c,j,i_met_lit) = &
ns_soil%decomp_npools_vr_col(c,j,i_met_lit) + nf_veg%gap_mortality_n_to_litr_met_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_cel_lit) = &
ns_soil%decomp_npools_vr_col(c,j,i_cel_lit) + nf_veg%gap_mortality_n_to_litr_cel_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_lig_lit) = &
ns_soil%decomp_npools_vr_col(c,j,i_lig_lit) + nf_veg%gap_mortality_n_to_litr_lig_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_litr2) = &
ns_soil%decomp_npools_vr_col(c,j,i_litr2) + nf_veg%gap_mortality_n_to_litr_cel_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_litr3) = &
ns_soil%decomp_npools_vr_col(c,j,i_litr3) + nf_veg%gap_mortality_n_to_litr_lig_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_cwd) = &
ns_soil%decomp_npools_vr_col(c,j,i_cwd) + nf_veg%gap_mortality_n_to_cwdn_col(c,j) * dt
end do
Expand Down Expand Up @@ -171,10 +171,10 @@ subroutine NStateUpdate2h(num_soilc, filter_soilc, num_soilp, filter_soilp, &
c = filter_soilc(fc)
ns_soil%decomp_npools_vr_col(c,j,i_met_lit) = &
ns_soil%decomp_npools_vr_col(c,j,i_met_lit) + nf_veg%harvest_n_to_litr_met_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_cel_lit) = &
ns_soil%decomp_npools_vr_col(c,j,i_cel_lit) + nf_veg%harvest_n_to_litr_cel_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_lig_lit) = &
ns_soil%decomp_npools_vr_col(c,j,i_lig_lit) + nf_veg%harvest_n_to_litr_lig_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_litr2) = &
ns_soil%decomp_npools_vr_col(c,j,i_litr2) + nf_veg%harvest_n_to_litr_cel_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_litr3) = &
ns_soil%decomp_npools_vr_col(c,j,i_litr3) + nf_veg%harvest_n_to_litr_lig_n_col(c,j) * dt
ns_soil%decomp_npools_vr_col(c,j,i_cwd) = &
ns_soil%decomp_npools_vr_col(c,j,i_cwd) + nf_veg%harvest_n_to_cwdn_col(c,j) * dt
end do
Expand Down
6 changes: 3 additions & 3 deletions src/biogeochem/CNNStateUpdate3Mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -10,7 +10,7 @@ module CNNStateUpdate3Mod
use clm_varpar , only: nlevdecomp, ndecomp_pools
use clm_time_manager , only : get_step_size_real
use clm_varctl , only : iulog, use_nitrif_denitrif
use clm_varpar , only : i_cwd, i_met_lit, i_cel_lit, i_lig_lit
use clm_varpar , only : i_cwd, i_met_lit, i_litr2, i_litr3
use CNVegNitrogenStateType , only : cnveg_nitrogenstate_type
use CNVegNitrogenFluxType , only : cnveg_nitrogenflux_type
use SoilBiogeochemNitrogenStateType , only : soilbiogeochem_nitrogenstate_type
Expand Down Expand Up @@ -86,9 +86,9 @@ subroutine NStateUpdate3(num_soilc, filter_soilc, num_soilp, filter_soilp, &
! patch-level wood to column-level litter (uncombusted wood)
ns_soil%decomp_npools_vr_col(c,j,i_met_lit) = ns_soil%decomp_npools_vr_col(c,j,i_met_lit) + &
nf_veg%m_n_to_litr_met_fire_col(c,j)* dt
ns_soil%decomp_npools_vr_col(c,j,i_cel_lit) = ns_soil%decomp_npools_vr_col(c,j,i_cel_lit) + &
ns_soil%decomp_npools_vr_col(c,j,i_litr2) = ns_soil%decomp_npools_vr_col(c,j,i_litr2) + &
nf_veg%m_n_to_litr_cel_fire_col(c,j)* dt
ns_soil%decomp_npools_vr_col(c,j,i_lig_lit) = ns_soil%decomp_npools_vr_col(c,j,i_lig_lit) + &
ns_soil%decomp_npools_vr_col(c,j,i_litr3) = ns_soil%decomp_npools_vr_col(c,j,i_litr3) + &
nf_veg%m_n_to_litr_lig_fire_col(c,j)* dt
end do ! end of column loop
end do
Expand Down
Loading