Computer program and input data set for performing comparative temporal analyses of heating infrastructure scenarios

The program developed for performing the comparative study of district heating and conventional heating infrastructures for the different feedstocks - heating oil, natural gas, forest residue and willow crop, is presented here. The program is developed in R programing language. The primary file is named as 'RF_temporalmodel_soildepth100.Rmd'. The 'rmarkdown' approach to the code is followed for distributing the code into different compartments. The 'myfunctions.R' contains function definitions for functions called in the primary 'rmarkdown' program code.
To successfully run the 'RF_temporalmodel_soildepth100.Rmd' code file, the variables saving the directory file paths need to be updated. Comments are provided at the specific code lines needing update within the code file.

Two cases of land use change for willow crop cultivation needs separate parameterization. Once the code chunk named ‘load_lib_func’ is run, analyst need to set the land use flag change variable ‘set_W_luc_case’ type of either ‘croplands_to_willow_soildepth_100cm_Mg_C_ha_yr’ or ‘pasture_hay_grasslands_to_willow_100cm_Mg_C_ha_yr’. Both lines of code are written, the needed code line be commented out to mark it not executable. In addition to setting the ‘set_W_luc_case’ variable, analyst need to enable the corresponding code line under the code chunk named ‘willow LUC cases’, and update the two variables ‘Wloc_DH_em_RF_2’ and ‘Wloc_DH_em_CRF_2’. Again, analyst can comment out the corresponding lines to reflect the land use change under study. Once these two variables are updated the analyst can run all the code chunks from ‘load_lib_func’ up to ‘set_W_luc_case’. Next, analyst can enable the alternative land use change case and re-run the code lines. To understand the algorithm of the executed code, the analyst is referred to the supplementary information A (SI-A) document.
At this stage, all the data matrices needing to plot the analyses graph for radiative forcing and cumulative forcing are populated. The code chunks beyond ‘willow LUC cases’ facilitate plotting the graphs, some of which are presented in the journal manuscript. The analyst can follow the code, as the commands are self-explanatory.

Please reach out for any question or clarification regarding the code and data.

本文介绍了一款用于对比研究不同原料（取暖用油、天然气、森林残留物及柳树种植园）的区域供暖与传统供暖基础设施的程序。该程序基于R语言开发，主文件名为`RF_temporalmodel_soildepth100.Rmd`。代码采用rmarkdown框架进行模块化分发，`myfunctions.R`文件包含主rmarkdown程序中调用的所有函数定义。

若要成功运行`RF_temporalmodel_soildepth100.Rmd`代码文件，需更新存储目录文件路径的变量。代码文件中已在需要修改的特定代码行添加注释说明。

针对柳树种植的两种土地利用变化场景，需进行单独的参数化设置。运行名为`load_lib_func`的代码块后，分析人员需设置土地利用变化标志变量`set_W_luc_case`，其可选类型为`croplands_to_willow_soildepth_100cm_Mg_C_ha_yr`或`pasture_hay_grasslands_to_willow_100cm_Mg_C_ha_yr`。两段相关代码均已预置，只需将无需执行的代码行注释即可。除设置`set_W_luc_case`变量外，分析人员还需启用`willow LUC cases`代码块下对应的代码行，并更新`Wloc_DH_em_RF_2`与`Wloc_DH_em_CRF_2`两个变量。同样，分析人员可通过注释对应代码行来标识当前研究的土地利用变化场景。完成这两个变量的更新后，分析人员可运行从`load_lib_func`到`set_W_luc_case`的所有代码块。随后，可启用另一种土地利用变化场景并重新运行对应代码行。若需了解所执行代码的算法逻辑，可参阅补充材料A（SI-A）文档。

至此，用于绘制辐射强迫与累积强迫分析图表所需的所有数据矩阵均已填充。`willow LUC cases`之后的代码块用于绘制图表，其中部分图表已发表于期刊论文中。代码中的命令均具备自解释性，分析人员可直接参考代码进行操作。

如有任何关于代码与数据的疑问或需要澄清之处，请随时联系。