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

User defined top snow layer #792

Merged
merged 23 commits into from
Sep 9, 2019
Merged

User defined top snow layer #792

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
23 commits
Select commit Hold shift + click to select a range
63b2c81
Changes from @olyson to fix #780
ekluzek Aug 16, 2019
20cbf74
Merge tag 'ctsm1.0.dev059' into snow_layers_user_defined
slevis-lmwg Aug 25, 2019
3e5856a
Allow users to specify snow layer structure in namelist
slevis-lmwg Aug 26, 2019
986cf49
Merge commit 'f2869fe202f71c4fd8c854fce47089630b4396c5' into clm50cos…
ekluzek Aug 28, 2019
7f67efc
Remove the commit to update the params file, which shouldn't be done …
ekluzek Aug 28, 2019
ee939a7
Merge tag 'ctsm1.0.dev061' into clm50coslatfirefix
ekluzek Sep 2, 2019
b2cf3a6
Changes for #787 for dwt_slash_flux
ekluzek Sep 2, 2019
f293334
Make sure urbantv is set for clm4_5 #175
ekluzek Sep 2, 2019
f49dac1
Revisions in response to @billsacks code review
slevis-lmwg Sep 2, 2019
47cd874
Merge tag 'ctsm1.0.dev061' into user_defined_top_snow_layer
slevis-lmwg Sep 2, 2019
ae415da
Merge tag 'ctsm1.0.dev062' into clm50coslatfirefix
ekluzek Sep 4, 2019
d31b2ee
Fix urbantv to point to different file for clm4_5 than clm5_0 as it s…
ekluzek Sep 4, 2019
ffe027c
Remove DWT_SLASH_* from output usermod at scales finer than gridcell,…
ekluzek Sep 6, 2019
ec4d91a
Update changelog
ekluzek Sep 6, 2019
5bcf4d6
Merge tag 'ctsm1.0.dev063' into user_defined_top_snow_layer
slevis-lmwg Sep 7, 2019
e00da7c
Revisions in response to code review
slevis-lmwg Sep 7, 2019
d7292f3
First draft of ChangeLog entry
slevis-lmwg Sep 7, 2019
e478244
Updated temporary test code checking bit-for-bit
slevis-lmwg Sep 7, 2019
db4bb80
Updated temporary test code checking bit-for-bit (part 2)
slevis-lmwg Sep 7, 2019
97cef78
Removing temporary test code
slevis-lmwg Sep 7, 2019
bd7f647
Updates needed for snowhydrology unit tests to pass
slevis-lmwg Sep 9, 2019
4c19afc
Updated ChangeLog/Sum
slevis-lmwg Sep 9, 2019
1d5ea12
Minor edit in ChangeLog
billsacks Sep 9, 2019
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
42 changes: 42 additions & 0 deletions bld/CLMBuildNamelist.pm
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3521,6 +3521,48 @@ sub setup_logic_snowpack {
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'reset_snow');
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'reset_snow_glc');
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'reset_snow_glc_ela');
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'snow_dzmin_1');
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'snow_dzmin_2');
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'snow_dzmax_l_1');
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'snow_dzmax_l_2');
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'snow_dzmax_u_1');
add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'snow_dzmax_u_2');

my $dzmin1 = $nl->get_value('snow_dzmin_1');
my $dzmin2 = $nl->get_value('snow_dzmin_2');
my $dzmax_l1 = $nl->get_value('snow_dzmax_l_1');
my $dzmax_l2 = $nl->get_value('snow_dzmax_l_2');
my $dzmax_u1 = $nl->get_value('snow_dzmax_u_1');
my $dzmax_u2 = $nl->get_value('snow_dzmax_u_2');

if ($dzmin1 != 0.01 || $dzmin2 != 0.015 || $dzmax_u1 != 0.02 || $dzmax_u2 != 0.05 || $dzmax_l1 != 0.03 || $dzmax_l2 != 0.07) {
$log->warning("Setting any of the following namelist variables to NON DEFAULT values remains untested as of Sep 6, 2019: snow_dzmin_1 & 2, snow_dzmax_u_1 & 2, snow_dzmax_l_1 & 2." );
$log->warning("Leave these variables unspecified in user_nl_clm in order to use the default values." );
}
if ($dzmin1 <= 0.0 || $dzmin2 <= 0.0 || $dzmax_u1 <= 0.0 || $dzmax_u2 <= 0.0 || $dzmax_l1 <= 0.0 || $dzmax_l2 <= 0.0) {
$log->fatal_error('One or more of the snow_dzmin_* and/or snow_dzmax_* were set incorrectly to be <= 0');
}
if ($dzmin2 <= $dzmin1) {
$log->fatal_error('snow_dzmin_2 was set incorrectly to be <= snow_dzmin_1');
}
if ($dzmax_l2 <= $dzmax_l1) {
$log->fatal_error('snow_dzmax_l_2 was set incorrectly to be <= snow_dzmax_l_1');
}
if ($dzmax_u2 <= $dzmax_u1) {
$log->fatal_error('snow_dzmax_u_2 was set incorrectly to be <= snow_dzmax_u_1');
}
if ($dzmin1 >= $dzmax_u1) {
$log->fatal_error('snow_dzmin_1 was set incorrectly to be >= snow_dzmax_u_1');
}
if ($dzmin2 >= $dzmax_u2) {
$log->fatal_error('snow_dzmin_2 was set incorrectly to be >= snow_dzmax_u_2');
}
if ($dzmax_u1 >= $dzmax_l1) {
$log->fatal_error('snow_dzmax_u_1 was set incorrectly to be >= snow_dzmax_l_1');
}
if ($dzmax_u2 >= $dzmax_l2) {
$log->fatal_error('snow_dzmax_u_2 was set incorrectly to be >= snow_dzmax_l_2');
}

if (remove_leading_and_trailing_quotes($nl->get_value('snow_overburden_compaction_method')) eq 'Vionnet2012') {
# overburden_compress_tfactor isn't used if we're using the Vionnet2012
Expand Down
48 changes: 27 additions & 21 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 @@ -298,6 +298,13 @@ attributes from the config_cache.xml file (with keys converted to upper-case).
<h2osno_max phys="clm5_0" structure="fast" >5000.0</h2osno_max>
<h2osno_max phys="clm4_5" >1000.0</h2osno_max>

<snow_dzmin_1>0.010d00</snow_dzmin_1>
<snow_dzmin_2>0.015d00</snow_dzmin_2>
<snow_dzmax_l_1>0.03d00</snow_dzmax_l_1>
<snow_dzmax_l_2>0.07d00</snow_dzmax_l_2>
<snow_dzmax_u_1>0.02d00</snow_dzmax_u_1>
<snow_dzmax_u_2>0.05d00</snow_dzmax_u_2>

<int_snow_max phys="clm5_0">2000.</int_snow_max>
<!-- For clm4_5, make this effectively unlimited -->
<int_snow_max phys="clm4_5">1.e30</int_snow_max>
Expand Down Expand Up @@ -1105,36 +1112,35 @@ lnd/clm2/surfdata_map/surfdata_ne120np4_78pfts_CMIP6_simyr1850_c170824.nc</fsurd
<popdensmapalgo use_cn=".true." hgrid="5x5_amazon" >nn</popdensmapalgo>

<!-- Urban time varying streams namelist defaults -->
<stream_year_first_urbantv phys="clm4_5" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm4_5" >2000</stream_year_last_urbantv>

<stream_year_first_urbantv phys="clm5_0" sim_year="2000" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" sim_year="2000" >2000</stream_year_last_urbantv>
<stream_year_first_urbantv sim_year="2000" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv sim_year="2000" >2000</stream_year_last_urbantv>

<stream_year_first_urbantv phys="clm5_0" sim_year="1850" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" sim_year="1850" >1850</stream_year_last_urbantv>
<stream_year_first_urbantv sim_year="1850" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv sim_year="1850" >1850</stream_year_last_urbantv>

<stream_year_first_urbantv phys="clm5_0" sim_year="1000" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" sim_year="1000" >2000</stream_year_last_urbantv>
<stream_year_first_urbantv sim_year="1000" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv sim_year="1000" >2000</stream_year_last_urbantv>

<stream_year_first_urbantv phys="clm5_0" sim_year="constant" sim_year_range="1000-1002" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" sim_year="constant" sim_year_range="1000-1002" >2000</stream_year_last_urbantv>
<stream_year_first_urbantv sim_year="constant" sim_year_range="1000-1002" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv sim_year="constant" sim_year_range="1000-1002" >2000</stream_year_last_urbantv>

<stream_year_first_urbantv phys="clm5_0" sim_year="constant" sim_year_range="1000-1004" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" sim_year="constant" sim_year_range="1000-1004" >2000</stream_year_last_urbantv>
<stream_year_first_urbantv sim_year="constant" sim_year_range="1000-1004" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv sim_year="constant" sim_year_range="1000-1004" >2000</stream_year_last_urbantv>

<stream_year_first_urbantv phys="clm5_0" sim_year="constant" sim_year_range="1850-2000" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" sim_year="constant" sim_year_range="1850-2000" >2106</stream_year_last_urbantv>
<stream_year_first_urbantv sim_year="constant" sim_year_range="1850-2000" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv sim_year="constant" sim_year_range="1850-2000" >2106</stream_year_last_urbantv>

<stream_year_first_urbantv phys="clm5_0" sim_year="constant" sim_year_range="1850-2100" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" sim_year="constant" sim_year_range="1850-2100" >2106</stream_year_last_urbantv>
<stream_year_first_urbantv sim_year="constant" sim_year_range="1850-2100" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv sim_year="constant" sim_year_range="1850-2100" >2106</stream_year_last_urbantv>

<stream_year_first_urbantv phys="clm5_0" sim_year="constant" sim_year_range="2000-2100" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" sim_year="constant" sim_year_range="2000-2100" >2106</stream_year_last_urbantv>
<stream_year_first_urbantv sim_year="constant" sim_year_range="2000-2100" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv sim_year="constant" sim_year_range="2000-2100" >2106</stream_year_last_urbantv>

<stream_fldfilename_urbantv phys="clm5_0" hgrid="0.9x1.25" >lnd/clm2/urbandata/CLM50_tbuildmax_Oleson_2016_0.9x1.25_simyr1849-2106_c160923.nc</stream_fldfilename_urbantv>
<stream_fldfilename_urbantv phys="clm5_0" hgrid="0.9x1.25"
>lnd/clm2/urbandata/CLM50_tbuildmax_Oleson_2016_0.9x1.25_simyr1849-2106_c160923.nc</stream_fldfilename_urbantv>

<stream_fldfilename_urbantv phys="clm4_5" hgrid="0.9x1.25" >lnd/clm2/urbandata/CLM45_tbuildmax_Oleson_2016_0.9x1.25_simyr1849-2106_c160923.nc</stream_fldfilename_urbantv>
<stream_fldfilename_urbantv phys="clm4_5" hgrid="0.9x1.25"
>lnd/clm2/urbandata/CLM45_tbuildmax_Oleson_2016_0.9x1.25_simyr1849-2106_c160923.nc</stream_fldfilename_urbantv>

<urbantvmapalgo >nn</urbantvmapalgo>

Expand Down
27 changes: 27 additions & 0 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 @@ -2440,6 +2440,33 @@ Changes in this value should possibly be accompanied by changes in:
glc_snow_persistence_max_days > 0.
</entry>

<entry id="snow_dzmin_1" type="real" category="clm_physics"
group="clm_snowhydrology_inparm" valid_values="" >
Min snow thickness of layer 1 (top snow layer); values other than the default 0.01 have not been tested as of Sep 6, 2019
</entry>
<entry id="snow_dzmin_2" type="real" category="clm_physics"
group="clm_snowhydrology_inparm" valid_values="" >
Min snow thickness of layer 2; values other than the default 0.015 have not been tested as of Sep 6, 2019; snow_dzmin of remaining layers is generated with the following recursive formula: dzmin(j) = dzmax_u(j-1) * 0.5_r8
</entry>

<entry id="snow_dzmax_l_1" type="real" category="clm_physics"
group="clm_snowhydrology_inparm" valid_values="" >
Max snow thickness of layer 1 (top snow layer) when no layers beneath; values other than the default 0.03 have not been tested as of Sep 6, 2019
</entry>
<entry id="snow_dzmax_l_2" type="real" category="clm_physics"
group="clm_snowhydrology_inparm" valid_values="" >
Max snow thickness of layer 2 when no layers beneath; values other than the default 0.07 have not been tested as of Sep 6, 2019; snow_dzmax_l of remaining layers is generated with the following recursive formula: dzmax_l(j) = dzmax_u(j) + dzmax_l(j-1)
</entry>

<entry id="snow_dzmax_u_1" type="real" category="clm_physics"
group="clm_snowhydrology_inparm" valid_values="" >
Max snow thickness of layer 1 (top snow layer) when layers beneath; values other than the default 0.02 have not been tested as of Sep 6, 2019
</entry>
<entry id="snow_dzmax_u_2" type="real" category="clm_physics"
group="clm_snowhydrology_inparm" valid_values="" >
Max snow thickness of layer 2 when layers beneath; values other than the default 0.05 have not been tested as of Sep 6, 2019; snow_dzmax_u of remaining layers is generated with the following recursive formula: dzmax_u(j) = 2._r8 * dzmax_u(j-1) + 0.01_r8
</entry>

<entry id="int_snow_max" type="real" category="clm_physics"
group="scf_swenson_lawrence_2012_inparm" valid_values="" >
Limit applied to integrated snowfall when determining changes in snow-covered fraction during melt
Expand Down
4 changes: 4 additions & 0 deletions bld/namelist_files/use_cases/1850-2100_rcp2.6_transient.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -40,6 +40,10 @@
<stream_year_last_popdens phys="clm5_0" cnfireson=".true." >2010</stream_year_last_popdens>
<model_year_align_popdens phys="clm5_0" cnfireson=".true." >1850</model_year_align_popdens>

<stream_year_first_urbantv phys="clm4_5" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm4_5" >2100</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm4_5" >1850</model_year_align_urbantv>

<stream_year_first_urbantv phys="clm5_0" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" >2100</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm5_0" >1850</model_year_align_urbantv>
Expand Down
4 changes: 4 additions & 0 deletions bld/namelist_files/use_cases/1850-2100_rcp4.5_transient.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -40,6 +40,10 @@
<stream_year_last_popdens phys="clm5_0" cnfireson=".true." >2010</stream_year_last_popdens>
<model_year_align_popdens phys="clm5_0" cnfireson=".true." >1850</model_year_align_popdens>

<stream_year_first_urbantv phys="clm4_5" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm4_5" >2100</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm4_5" >1850</model_year_align_urbantv>

<stream_year_first_urbantv phys="clm5_0" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" >2100</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm5_0" >1850</model_year_align_urbantv>
Expand Down
4 changes: 4 additions & 0 deletions bld/namelist_files/use_cases/1850-2100_rcp6_transient.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -41,6 +41,10 @@
<stream_year_last_popdens phys="clm5_0" cnfireson=".true." >2010</stream_year_last_popdens>
<model_year_align_popdens phys="clm5_0" cnfireson=".true." >1850</model_year_align_popdens>

<stream_year_first_urbantv phys="clm4_5" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm4_5" >2100</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm4_5" >1850</model_year_align_urbantv>

<stream_year_first_urbantv phys="clm5_0" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" >2100</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm5_0" >1850</model_year_align_urbantv>
Expand Down
4 changes: 4 additions & 0 deletions bld/namelist_files/use_cases/1850-2100_rcp8.5_transient.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -40,6 +40,10 @@
<stream_year_last_popdens phys="clm5_0" cnfireson=".true." >2010</stream_year_last_popdens>
<model_year_align_popdens phys="clm5_0" cnfireson=".true." >1850</model_year_align_popdens>

<stream_year_first_urbantv phys="clm4_5" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm4_5" >2100</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm4_5" >1850</model_year_align_urbantv>

<stream_year_first_urbantv phys="clm5_0" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" >2100</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm5_0" >1850</model_year_align_urbantv>
Expand Down
3 changes: 3 additions & 0 deletions bld/namelist_files/use_cases/1850_control.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -23,6 +23,9 @@
<stream_year_first_popdens phys="clm5_0" cnfireson=".true." >1850</stream_year_first_popdens>
<stream_year_last_popdens phys="clm5_0" cnfireson=".true." >1850</stream_year_last_popdens>

<stream_year_first_urbantv phys="clm4_5" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" >1850</stream_year_last_urbantv>

<stream_year_first_urbantv phys="clm5_0" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" >1850</stream_year_last_urbantv>

Expand Down
4 changes: 4 additions & 0 deletions bld/namelist_files/use_cases/2000-2100_rcp8.5_transient.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -39,6 +39,10 @@
<stream_year_last_popdens phys="clm5_0" cnfireson=".true." >2010</stream_year_last_popdens>
<model_year_align_popdens phys="clm5_0" cnfireson=".true." >2000</model_year_align_popdens>

<stream_year_first_urbantv phys="clm4_5" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm4_5" >2100</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm4_5" >2000</model_year_align_urbantv>

<stream_year_first_urbantv phys="clm5_0" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" >2100</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm5_0" >2000</model_year_align_urbantv>
Expand Down
3 changes: 3 additions & 0 deletions bld/namelist_files/use_cases/2000_control.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -24,6 +24,9 @@
<stream_year_first_popdens phys="clm5_0" cnfireson=".true." >2000</stream_year_first_popdens>
<stream_year_last_popdens phys="clm5_0" cnfireson=".true." >2000</stream_year_last_popdens>

<stream_year_first_urbantv phys="clm4_5" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm4_5" >2000</stream_year_last_urbantv>

<stream_year_first_urbantv phys="clm5_0" >2000</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" >2000</stream_year_last_urbantv>

Expand Down
4 changes: 4 additions & 0 deletions bld/namelist_files/use_cases/20thC_transient.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -34,6 +34,10 @@
<stream_year_last_popdens phys="clm5_0" cnfireson=".true." >2016</stream_year_last_popdens>
<model_year_align_popdens phys="clm5_0" cnfireson=".true." >1850</model_year_align_popdens>

<stream_year_first_urbantv phys="clm4_5" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm4_5" >2106</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm4_5" >1850</model_year_align_urbantv>

<stream_year_first_urbantv phys="clm5_0" >1850</stream_year_first_urbantv>
<stream_year_last_urbantv phys="clm5_0" >2106</stream_year_last_urbantv>
<model_year_align_urbantv phys="clm5_0" >1850</model_year_align_urbantv>
Expand Down
2 changes: 1 addition & 1 deletion cime_config/usermods_dirs/output_bgc/user_nl_clm
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@ hist_fincl2 += 'GPP', 'NPP', 'AGNPP', 'TOTVEGC', 'NPP_NUPTAKE', 'AR', 'HR', 'HTO

! h2 stream (monthly average, landunit-level)
! TOT_WOODPRODC:I, CROPPROD1C:I, and NEE are not available at the landunit level
hist_fincl3 += 'GPP', 'NPP', 'AR', 'HR', 'DWT_CONV_CFLUX_PATCH', 'WOOD_HARVESTC', 'DWT_WOOD_PRODUCTC_GAIN_PATCH', 'SLASH_HARVESTC', 'COL_FIRE_CLOSS', 'DWT_SLASH_CFLUX', 'FROOTC:I', 'HTOP'
hist_fincl3 += 'GPP', 'NPP', 'AR', 'HR', 'DWT_CONV_CFLUX_PATCH', 'WOOD_HARVESTC', 'DWT_WOOD_PRODUCTC_GAIN_PATCH', 'SLASH_HARVESTC', 'COL_FIRE_CLOSS', 'FROOTC:I', 'HTOP'

! h3 stream (yearly average, gridcell-level)
hist_fincl4 += 'SOILC_vr', 'SOILN_vr', 'CWDC_vr', 'LITR1C_vr', 'LITR2C_vr', 'LITR3C_vr', 'LITR1N_vr', 'LITR2N_vr', 'LITR3N_vr','CWDN_vr', 'TOTLITC:I', 'TOT_WOODPRODC:I', 'TOTSOMC:I','TOTVEGC:I'
Expand Down
Loading