|
89 | 89 | "import pylab\n", |
90 | 90 | "import matplotlib.dates as mdates\n", |
91 | 91 | "from mpl_toolkits.axes_grid1.inset_locator import inset_axes\n", |
| 92 | + "import gc # Import the garbage collection module\n", |
92 | 93 | "\n", |
93 | 94 | "import pyet\n", |
94 | 95 | "import spotpy\n", |
|
120 | 121 | "source": [ |
121 | 122 | "#edit default values for plots\n", |
122 | 123 | "params = {\n", |
123 | | - " \"font.family\": \"Arial\",\n", |
124 | 124 | " \"legend.fontsize\": \"15\",\n", |
125 | 125 | " \"axes.labelsize\": \"16\",\n", |
126 | 126 | " \"xtick.labelsize\": \"15\",\n", |
|
773 | 773 | "fig.supxlabel(\"PET$_{benchmark}$ [mm/day]\", x=0.475, fontsize=16)\n", |
774 | 774 | "fig.supylabel(\"PET$_{pyet}$ [mm/day]\", fontsize=16)\n", |
775 | 775 | "fig.subplots_adjust(wspace=0.05, hspace=0.05, left=0.05)\n", |
776 | | - "plt.tight_layout()\n", |
| 776 | + "plt.tight_layout();\n", |
777 | 777 | "#fig.savefig(\"figure1.png\", dpi=600, bbox_inches=\"tight\")" |
778 | 778 | ] |
779 | 779 | }, |
|
872 | 872 | "\n", |
873 | 873 | "plt.tight_layout()\n", |
874 | 874 | "\n", |
875 | | - "axs[2].legend(loc=(1.05,0.1), ncol=2, bbox_transform=axs[0].transAxes, fontsize=8)\n", |
| 875 | + "axs[2].legend(loc=(1.05,0.1), ncol=2, bbox_transform=axs[0].transAxes, fontsize=8);\n", |
876 | 876 | "#fig.savefig(\"figure2.png\", dpi=600, bbox_inches=\"tight\")" |
877 | 877 | ] |
878 | 878 | }, |
|
1007 | 1007 | " cbar.set_label(\"PET [mm/day]\", labelpad=10)\n", |
1008 | 1008 | " plt.subplots_adjust(hspace=0.02, wspace=0.1)\n", |
1009 | 1009 | "\n", |
1010 | | - "plt.subplots_adjust(hspace=0.02, wspace=0.01)\n", |
| 1010 | + "plt.subplots_adjust(hspace=0.02, wspace=0.01);\n", |
1011 | 1011 | "\n", |
1012 | 1012 | "#fig.savefig(\"figure3.png\", dpi=300, bbox_inches=\"tight\")" |
1013 | 1013 | ] |
|
1287 | 1287 | "axs[1,5].set_yticklabels([\"\",800,1000,\"\"])\n", |
1288 | 1288 | "axs[0,5].legend(loc=(-3.65,1.05), ncol=7, bbox_transform=axs[1,5].transAxes)\n", |
1289 | 1289 | "clb = fig.colorbar(density, orientation=\"horizontal\", ax=axs[0, 0], cax = axs[0, 0].inset_axes([0.04, 1.2, 1, 0.05]))\n", |
1290 | | - "clb.ax.set_title(\"Number of points per pixel\")\n", |
| 1290 | + "clb.ax.set_title(\"Number of points per pixel\");\n", |
1291 | 1291 | "#fig.savefig(\"figure4.png\", dpi=600, bbox_inches=\"tight\")" |
1292 | 1292 | ] |
1293 | 1293 | }, |
|
1336 | 1336 | "source": [ |
1337 | 1337 | "# Load CO2 data for each RCP scenario\n", |
1338 | 1338 | "rcp_co2 = pd.DataFrame()\n", |
1339 | | - "rcp_co2[\"rcp_26\"] = pd.read_csv(\"data/example_10/co2_conc/RCP3PD_MIDYR_CONC.dat\", skiprows=38, \n", |
| 1339 | + "rcp_co2[\"rcp_26\"] = pd.read_csv(\"data/example_10/co2_conc/RCP3PD_MIDYR_CONC.DAT\", skiprows=38, \n", |
1340 | 1340 | " delim_whitespace=True, index_col=\"YEARS\").loc[\"2020\":\"2100\", \"CO2\"]\n", |
1341 | | - "rcp_co2[\"rcp_45\"] = pd.read_csv(\"data/example_10/co2_conc/RCP45_MIDYR_CONC.dat\", skiprows=38, \n", |
| 1341 | + "rcp_co2[\"rcp_45\"] = pd.read_csv(\"data/example_10/co2_conc/RCP45_MIDYR_CONC.DAT\", skiprows=38, \n", |
1342 | 1342 | " delim_whitespace=True, index_col=\"YEARS\").loc[\"2020\":\"2100\", \"CO2\"]\n", |
1343 | | - "rcp_co2[\"rcp_60\"] = pd.read_csv(\"data/example_10/co2_conc/RCP6_MIDYR_CONC.dat\", skiprows=38, \n", |
| 1343 | + "rcp_co2[\"rcp_60\"] = pd.read_csv(\"data/example_10/co2_conc/RCP6_MIDYR_CONC.DAT\", skiprows=38, \n", |
1344 | 1344 | " delim_whitespace=True, index_col=\"YEARS\").loc[\"2020\":\"2100\", \"CO2\"]\n", |
1345 | | - "rcp_co2[\"rcp_85\"] = pd.read_csv(\"data/example_10/co2_conc/RCP85_MIDYR_CONC.dat\", skiprows=38, \n", |
| 1345 | + "rcp_co2[\"rcp_85\"] = pd.read_csv(\"data/example_10/co2_conc/RCP85_MIDYR_CONC.DAT\", skiprows=38, \n", |
1346 | 1346 | " delim_whitespace=True, index_col=\"YEARS\").loc[\"2020\":\"2100\", \"CO2\"]" |
1347 | 1347 | ] |
1348 | 1348 | }, |
|
1416 | 1416 | "dpet_rcp_etco2 = pd.DataFrame(index=rcp_temp.index, columns=[\"rcp_26\", \"rcp_45\", \"rcp_60\", \"rcp_85\"])\n", |
1417 | 1417 | "dpet_rcp_etco2_5th = pd.DataFrame(index=rcp_temp.index, columns=[\"rcp_26\", \"rcp_45\", \"rcp_60\", \"rcp_85\"])\n", |
1418 | 1418 | "dpet_rcp_etco2_95th = pd.DataFrame(index=rcp_temp.index, columns=[\"rcp_26\", \"rcp_45\", \"rcp_60\", \"rcp_85\"])\n", |
1419 | | - "for year in rcp_temp.index:\n", |
| 1419 | + "for year in rcp_temp.index: \n", |
1420 | 1420 | " for rcp in [\"rcp_26\", \"rcp_45\", \"rcp_60\", \"rcp_85\"]:\n", |
1421 | 1421 | " df_rcp_et = pd.DataFrame()\n", |
1422 | 1422 | " df_rcp_etco2 = pd.DataFrame()\n", |
1423 | | - " for i in np.arange(0, len(methods)):\n", |
| 1423 | + " for i in np.arange(0, len(methods[:2])): # only for two methods to reduce processing for RTD to handle!!!!!!!!!!!!\n", |
1424 | 1424 | " input1 = input_rcp(rcp_temp.loc[year, rcp])\n", |
1425 | 1425 | " df_rcp_et[methods[i]] = simulate(sollutions2[i], methods[i], input1[i])\n", |
1426 | 1426 | " df_rcp_etco2[methods[i]] = simulate(sollutions2[i], methods[i], \n", |
|
1430 | 1430 | " dpet_rcp_et_95th.loc[year, rcp] = df_rcp_et.resample(\"y\").mean().mean().quantile(0.95)\n", |
1431 | 1431 | " dpet_rcp_etco2.loc[year, rcp] = df_rcp_etco2.resample(\"y\").mean().mean().mean()\n", |
1432 | 1432 | " dpet_rcp_etco2_5th.loc[year, rcp] = df_rcp_etco2.resample(\"y\").mean().mean().quantile(0.05)\n", |
1433 | | - " dpet_rcp_etco2_95th.loc[year, rcp] = df_rcp_etco2.resample(\"y\").mean().mean().quantile(0.95)" |
| 1433 | + " dpet_rcp_etco2_95th.loc[year, rcp] = df_rcp_etco2.resample(\"y\").mean().mean().quantile(0.95)\n", |
| 1434 | + " # Call gc.collect() after processing each RCP scenario for a year\n", |
| 1435 | + " gc.collect()" |
1434 | 1436 | ] |
1435 | 1437 | }, |
1436 | 1438 | { |
|
1487 | 1489 | "\n", |
1488 | 1490 | "for i, letter in enumerate([\"a\", \"b\", \"c\", \"d\"]):\n", |
1489 | 1491 | " axs[i].text(0.85, 0.9, \"({})\".format(letter), transform=axs[i].transAxes, fontsize=12)\n", |
1490 | | - " axs[i].tick_params(axis='both', which='major', labelsize=13)\n", |
| 1492 | + " axs[i].tick_params(axis='both', which='major', labelsize=13);\n", |
1491 | 1493 | "\n", |
1492 | 1494 | "#fig.savefig(\"figure5.png\", dpi=600, bbox_inches=\"tight\")" |
1493 | 1495 | ] |
|
1529 | 1531 | "pet_df.plot() # daily PET [mm/day]\n", |
1530 | 1532 | "pet_df.boxplot() # boxplot PET[mm/day]\n", |
1531 | 1533 | "pet_df.cumsum().plot() # cummulative PET [mm]\n", |
1532 | | - "plt.scatter(pyet_makkink, pet_knmi) # plot Makkink pyet vs KNMI" |
| 1534 | + "plt.scatter(pyet_makkink, pet_knmi) # plot Makkink pyet vs KNMI;" |
1533 | 1535 | ] |
1534 | 1536 | }, |
1535 | 1537 | { |
|
0 commit comments