diff --git a/docs/tutorials/PVSC 2021 Fig. 3.ipynb b/docs/tutorials/PVSC 2021 Fig. 3.ipynb deleted file mode 100644 index 0ec83f7b..00000000 --- a/docs/tutorials/PVSC 2021 Fig. 3.ipynb +++ /dev/null @@ -1,1841 +0,0 @@ -{ - "cells": [ - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "# PVSC Fig. 3" - ] - }, - { - "cell_type": "code", - "execution_count": 2, - "metadata": {}, - "outputs": [], - "source": [ - "import PV_ICE\n", - "import numpy as np\n", - "import pandas as pd\n", - "import os,sys\n", - "import matplotlib.pyplot as plt\n", - "plt.rcParams.update({'font.size': 22})\n", - "plt.rcParams['figure.figsize'] = (12, 8)" - ] - }, - { - "cell_type": "code", - "execution_count": 3, - "metadata": {}, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "Your simulation will be stored in C:\\Users\\sayala\\Documents\\GitHub\\CircularEconomy-MassFlowCalculator\\PV_ICE\\TEMP\n" - ] - } - ], - "source": [ - "import os\n", - "from pathlib import Path\n", - "\n", - "testfolder = str(Path().resolve().parent.parent / 'PV_ICE' / 'TEMP')\n", - "\n", - "print (\"Your simulation will be stored in %s\" % testfolder)\n" - ] - }, - { - "cell_type": "code", - "execution_count": 18, - "metadata": {}, - "outputs": [], - "source": [ - "years = list(range(2009,2051))\n" - ] - }, - { - "cell_type": "code", - "execution_count": 9, - "metadata": {}, - "outputs": [], - "source": [ - "USyearly = pd.read_csv(os.path.join(testfolder,'PVSC_USYearly_DataforPlot.csv'))\n" - ] - }, - { - "cell_type": "code", - "execution_count": 12, - "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "Index(['Unnamed: 0', 'Waste_glass_Reference.Mod',\n", - " 'Waste_silicon_Reference.Mod', 'Waste_silver_Reference.Mod',\n", - " 'Waste_copper_Reference.Mod', 'Waste_aluminum_Reference.Mod',\n", - " 'Waste_glass_95-by-35.Adv', 'Waste_silicon_95-by-35.Adv',\n", - " 'Waste_silver_95-by-35.Adv', 'Waste_copper_95-by-35.Adv',\n", - " 'Waste_aluminum_95-by-35.Adv', 'Waste_glass_95-by-35_Elec.Adv_DR',\n", - " 'Waste_silicon_95-by-35_Elec.Adv_DR',\n", - " 'Waste_silver_95-by-35_Elec.Adv_DR',\n", - " 'Waste_copper_95-by-35_Elec.Adv_DR',\n", - " 'Waste_aluminum_95-by-35_Elec.Adv_DR',\n", - " 'VirginStock_glass_Reference.Mod', 'VirginStock_silicon_Reference.Mod',\n", - " 'VirginStock_silver_Reference.Mod', 'VirginStock_copper_Reference.Mod',\n", - " 'VirginStock_aluminum_Reference.Mod', 'VirginStock_glass_95-by-35.Adv',\n", - " 'VirginStock_silicon_95-by-35.Adv', 'VirginStock_silver_95-by-35.Adv',\n", - " 'VirginStock_copper_95-by-35.Adv', 'VirginStock_aluminum_95-by-35.Adv',\n", - " 'VirginStock_glass_95-by-35_Elec.Adv_DR',\n", - " 'VirginStock_silicon_95-by-35_Elec.Adv_DR',\n", - " 'VirginStock_silver_95-by-35_Elec.Adv_DR',\n", - " 'VirginStock_copper_95-by-35_Elec.Adv_DR',\n", - " 'VirginStock_aluminum_95-by-35_Elec.Adv_DR'],\n", - " dtype='object')" - ] - }, - "execution_count": 12, - "metadata": {}, - "output_type": "execute_result" - } - ], - "source": [ - "USyearly.keys()" - ] - }, - { - "cell_type": "code", - "execution_count": 19, - "metadata": {}, - "outputs": [ - { - "data": { - "text/html": [ - "
\n", - "\n", - "\n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - "
Unnamed: 0Waste_glass_Reference.ModWaste_silicon_Reference.ModWaste_silver_Reference.ModWaste_copper_Reference.ModWaste_aluminum_Reference.ModWaste_glass_95-by-35.AdvWaste_silicon_95-by-35.AdvWaste_silver_95-by-35.AdvWaste_copper_95-by-35.Adv...VirginStock_glass_95-by-35.AdvVirginStock_silicon_95-by-35.AdvVirginStock_silver_95-by-35.AdvVirginStock_copper_95-by-35.AdvVirginStock_aluminum_95-by-35.AdvVirginStock_glass_95-by-35_Elec.Adv_DRVirginStock_silicon_95-by-35_Elec.Adv_DRVirginStock_silver_95-by-35_Elec.Adv_DRVirginStock_copper_95-by-35_Elec.Adv_DRVirginStock_aluminum_95-by-35_Elec.Adv_DR
000.0012230.0006284.285685e-072.017821e-080.000000e+000.0012230.0006284.285685e-072.017821e-08...8.025386e+031036.03900723.7642474.4853867.738673e+028.025386e+031.036039e+0323.7642474.4853867.738673e+02
110.0011240.0002862.009061e-071.964373e-087.984519e-120.0011240.0002862.009061e-071.964373e-08...7.713247e+03640.16298711.1574074.3717557.533675e+027.713247e+036.401630e+0211.1574074.3717557.533675e+02
220.0135940.0034132.024336e-062.375177e-074.722305e-100.0135940.0034132.024336e-062.375177e-07...9.326221e+047689.440032111.99710052.8272439.109097e+039.326221e+047.689440e+03111.99710052.8272439.109097e+03
330.0130840.0034231.558864e-062.286382e-075.549092e-090.0130840.0034231.558864e-062.286382e-07...8.976487e+047557.01503586.52493650.8411918.767506e+038.976487e+047.557015e+0386.52493650.8411918.767506e+03
440.0419590.0108473.640602e-067.628068e-073.735847e-080.0419590.0108473.640602e-067.628068e-07...3.015900e+0524561.064097203.098965169.7585822.924719e+043.015900e+052.456106e+04203.098965169.7585822.924719e+04
550.0410210.0103123.244131e-067.322020e-071.855378e-070.0410210.0103123.244131e-067.322020e-07...2.948416e+0523449.823276181.582129162.5722532.804289e+042.948416e+052.344982e+04181.582129162.5722532.804289e+04
660.0868050.0229376.211354e-061.783707e-067.133639e-070.0868050.0229376.211354e-061.783707e-06...6.551432e+0552589.010236347.870650396.2874536.340030e+046.551432e+055.258901e+04347.870650396.2874536.340030e+04
770.0801500.0202105.593786e-061.897091e-062.236392e-060.0801500.0202105.593786e-061.897091e-06...6.367041e+0549194.730224312.035735419.0250586.268391e+046.367041e+054.919473e+04312.035735419.0250586.268391e+04
880.1068700.0226797.980107e-063.157919e-066.008380e-060.1068700.0226797.980107e-063.157919e-06...8.960064e+0563216.628971442.471248694.6379348.945712e+048.960064e+056.321663e+04442.471248694.6379348.945712e+04
990.1090190.0177838.192585e-063.177004e-061.437589e-050.1090190.0177838.192585e-063.177004e-06...9.134817e+0555380.797254448.126733688.8030668.572974e+049.134817e+055.538080e+04448.126733688.8030668.572974e+04
10100.1958540.0291671.413817e-056.304032e-063.149377e-050.1949300.0290291.407329e-056.275057e-06...1.632957e+0694182.347348765.7833241356.7600121.523814e+051.645299e+069.489416e+04771.5710051367.0142121.535330e+05
11110.1774090.0248271.231834e-056.267344e-066.441594e-050.1765730.0247101.226383e-056.239381e-06...1.563819e+0684123.648855643.3912031309.3500481.476974e+051.575638e+068.475944e+04648.2538631319.2459301.488137e+05
12120.1874960.0259551.273085e-056.949641e-061.246754e-040.1863080.0257901.265886e-056.909687e-06...1.646076e+0686301.832023629.4433981385.7291341.554664e+051.667130e+068.740568e+04637.4943461403.4534191.574549e+05
13130.1851510.0253361.276931e-057.452795e-062.303997e-040.1839840.0251761.270379e-057.413310e-06...1.618034e+0682645.936987572.8630651369.5016261.528179e+051.638730e+068.370303e+04580.1903181387.0183511.547726e+05
14140.2467760.0311181.694638e-051.054243e-054.088584e-040.3885970.0491392.440697e-051.538162e-05...3.422864e+06168191.7901931135.5638652922.0811883.232780e+053.145681e+061.545716e+051043.6061092685.4515822.970990e+05
15150.2462020.0279461.839608e-051.212721e-056.991525e-040.3863490.0440502.524480e-051.695016e-05...3.382475e+06158970.3977521042.4334142912.2826383.194635e+053.108563e+061.460970e+05958.0173452676.4465172.935934e+05
16160.3793030.0410872.733941e-051.874639e-051.154571e-030.7100850.0773754.260399e-052.977600e-05...6.225730e+06284772.5752361811.8510805193.7380265.879995e+056.343689e+062.901681e+051846.1801105292.1434565.991403e+05
17170.3814530.0393903.128039e-052.216655e-051.844060e-030.7092490.0739054.556430e-053.274601e-05...6.169521e+06275074.2743711695.4513814981.7881965.826909e+056.286415e+062.802861e+051727.5749925076.1778265.937311e+05
18180.5014650.0498384.171634e-053.108794e-052.852477e-030.8847370.0887865.748299e-054.304799e-05...7.656610e+06333249.4221561986.3704665977.7439947.231415e+058.816110e+063.837159e+052287.1817996882.9996688.326524e+05
19190.5134610.0485614.902133e-053.881457e-054.279934e-030.8971740.0860146.372512e-055.026210e-05...7.665399e+06325531.2123021852.4596735721.4445867.178539e+058.826231e+063.748289e+052132.9918666587.8868738.265640e+05
20200.7205900.0659666.542370e-055.547309e-056.241336e-031.1548490.1070038.099552e-056.794914e-05...9.813826e+06410069.1872182219.2481287052.8552879.190513e+051.207558e+075.045760e+052730.7083468678.2953881.130861e+06
21210.7376740.0646257.666464e-057.055613e-058.866731e-031.1692730.1040909.108346e-058.248517e-05...9.753521e+06400950.6996452054.5721326739.0185709.134043e+051.200138e+074.933560e+052528.0802078292.1301231.123913e+06
22220.7950360.0687459.334149e-059.261527e-051.229902e-021.3415900.1187211.116072e-041.077438e-04...1.103668e+07453699.4366312324.8699147625.5981881.033571e+061.379303e+075.670082e+052905.4923079530.0458471.291699e+06
23230.8269510.0702151.121099e-041.201568e-041.667688e-021.3706810.1199291.303008e-041.352746e-04...1.097820e+07451295.5714492312.5519057585.1949861.028095e+061.371994e+075.640039e+052890.0979499479.5521861.284855e+06
24240.9052770.0755701.358183e-041.559843e-042.208345e-021.6447120.1431701.605935e-041.766693e-04...1.301877e+07535180.3370002742.3985228995.0964861.219192e+061.791196e+077.363309e+053773.14463112375.9547501.677433e+06
25250.9540530.0778761.618118e-041.980445e-042.845403e-021.6904310.1451761.865898e-042.190043e-04...1.295812e+07532687.8858142729.6265768953.2043661.213514e+061.782853e+077.329017e+053755.57227612318.3172381.669621e+06
26261.0506260.0842961.910774e-042.470661e-043.546196e-022.0547990.1760422.250153e-042.761605e-04...1.567521e+07644383.3684703301.98229410830.5372461.467967e+062.149434e+078.835975e+054527.77543914851.1518722.012920e+06
27271.1060090.0868622.181750e-042.976550e-044.252520e-022.1075620.1782882.524408e-043.281095e-04...1.561007e+07641706.2635663288.26413610785.5415411.461869e+062.140503e+078.799265e+054508.96469614789.4524342.004557e+06
28281.2486320.0973652.484255e-043.522708e-044.947745e-022.2689010.1903512.840788e-043.858802e-04...1.650213e+07678377.9275313476.17895711401.9041631.545410e+062.374752e+079.762230e+055002.41186916407.9644402.223930e+06
29291.2977420.0994442.671787e-044.017728e-045.561036e-022.3192730.1922783.042125e-044.400673e-04...1.644024e+07675834.1498483463.14400111359.1493671.539615e+062.365846e+079.725623e+054983.65385416346.4379492.215591e+06
30301.4013050.1076062.687268e-044.333743e-045.773818e-022.3594960.1941933.059051e-044.778077e-04...1.653371e+07679676.9588833482.83552011423.7377611.548370e+062.395775e+079.848659e+055046.70055216553.2317952.243619e+06
31311.4399550.1088722.779371e-044.798020e-046.241551e-022.4091550.1957143.191703e-045.374668e-04...1.647720e+07677354.5473623470.93490011384.7035991.543079e+062.387587e+079.815007e+055029.45630816496.6704892.235953e+06
32321.5232230.1151432.712952e-045.101453e-046.340878e-022.4544870.1973143.170890e-045.872917e-04...1.656200e+07680840.7326663488.79898911443.2980041.551021e+062.514029e+071.033479e+065295.80744717370.3050752.354365e+06
33331.5814390.1168812.854872e-045.818602e-047.018165e-022.5451140.2001503.413036e-046.915515e-04...1.650999e+07678704.0101663477.84988611407.3848291.546153e+062.506137e+071.030236e+065279.18730317315.7908262.346976e+06
34341.7964220.1314963.076460e-046.656371e-047.877043e-022.6985860.2063253.747291e-048.210428e-04...1.667850e+07685632.5044593513.35323211523.8361821.561937e+062.656348e+071.091986e+065595.61270418353.6695242.487650e+06
35351.8784860.1338843.189522e-047.535916e-048.860343e-022.8559730.2114714.078259e-049.800827e-04...1.663006e+07683642.9449133503.15822911490.3964631.557404e+062.648635e+071.088818e+065579.37542718300.4110832.480432e+06
36362.2011410.1510724.007223e-049.796619e-041.164729e-013.2681460.2309625.196237e-041.306075e-03...1.718313e+07706382.5925283619.68189711872.5953411.609208e+062.901691e+071.192849e+066112.45615420048.9215682.717424e+06
37372.4725020.1615414.852152e-041.233090e-031.483842e-013.6892970.2470056.457871e-041.696850e-03...1.713665e+07704476.5008653609.91460411840.5585171.604865e+062.893849e+071.189630e+066095.96239819994.8218742.710091e+06
38382.9351360.1878085.752492e-041.510840e-031.833232e-014.3880770.2840287.915824e-042.156751e-03...1.887736e+07776039.5296233976.62154513043.3612061.767893e+063.022863e+071.242671e+066367.76077520886.3234572.830925e+06
39393.2498340.1992386.546092e-041.793002e-032.202285e-014.9592500.3049889.393534e-042.668277e-03...1.882960e+07774082.6437413966.59397113010.4706531.763435e+063.015224e+071.239538e+066351.70362820833.6558422.823787e+06
40403.6336490.2159097.281388e-042.078377e-032.585010e-015.6660960.3340591.094822e-033.231128e-03...1.936071e+07795922.6070824078.50743013377.5479991.813188e+063.221011e+071.324142e+066785.23437922255.6413423.016522e+06
41413.9776720.2278307.933713e-042.362686e-032.982925e-016.3800550.3596731.252286e-033.833597e-03...1.936052e+07795922.6070824078.50743013377.5479991.813188e+063.220988e+071.324142e+066785.23437922255.6413423.016522e+06
\n", - "

42 rows × 31 columns

\n", - "
" - ], - "text/plain": [ - " Unnamed: 0 Waste_glass_Reference.Mod Waste_silicon_Reference.Mod \\\n", - "0 0 0.001223 0.000628 \n", - "1 1 0.001124 0.000286 \n", - "2 2 0.013594 0.003413 \n", - "3 3 0.013084 0.003423 \n", - "4 4 0.041959 0.010847 \n", - "5 5 0.041021 0.010312 \n", - "6 6 0.086805 0.022937 \n", - "7 7 0.080150 0.020210 \n", - "8 8 0.106870 0.022679 \n", - "9 9 0.109019 0.017783 \n", - "10 10 0.195854 0.029167 \n", - "11 11 0.177409 0.024827 \n", - "12 12 0.187496 0.025955 \n", - "13 13 0.185151 0.025336 \n", - "14 14 0.246776 0.031118 \n", - "15 15 0.246202 0.027946 \n", - "16 16 0.379303 0.041087 \n", - "17 17 0.381453 0.039390 \n", - "18 18 0.501465 0.049838 \n", - "19 19 0.513461 0.048561 \n", - "20 20 0.720590 0.065966 \n", - "21 21 0.737674 0.064625 \n", - "22 22 0.795036 0.068745 \n", - "23 23 0.826951 0.070215 \n", - "24 24 0.905277 0.075570 \n", - "25 25 0.954053 0.077876 \n", - "26 26 1.050626 0.084296 \n", - "27 27 1.106009 0.086862 \n", - "28 28 1.248632 0.097365 \n", - "29 29 1.297742 0.099444 \n", - "30 30 1.401305 0.107606 \n", - "31 31 1.439955 0.108872 \n", - "32 32 1.523223 0.115143 \n", - "33 33 1.581439 0.116881 \n", - "34 34 1.796422 0.131496 \n", - "35 35 1.878486 0.133884 \n", - "36 36 2.201141 0.151072 \n", - "37 37 2.472502 0.161541 \n", - "38 38 2.935136 0.187808 \n", - "39 39 3.249834 0.199238 \n", - "40 40 3.633649 0.215909 \n", - "41 41 3.977672 0.227830 \n", - "\n", - " Waste_silver_Reference.Mod Waste_copper_Reference.Mod \\\n", - "0 4.285685e-07 2.017821e-08 \n", - "1 2.009061e-07 1.964373e-08 \n", - "2 2.024336e-06 2.375177e-07 \n", - "3 1.558864e-06 2.286382e-07 \n", - "4 3.640602e-06 7.628068e-07 \n", - "5 3.244131e-06 7.322020e-07 \n", - "6 6.211354e-06 1.783707e-06 \n", - "7 5.593786e-06 1.897091e-06 \n", - "8 7.980107e-06 3.157919e-06 \n", - "9 8.192585e-06 3.177004e-06 \n", - "10 1.413817e-05 6.304032e-06 \n", - "11 1.231834e-05 6.267344e-06 \n", - "12 1.273085e-05 6.949641e-06 \n", - "13 1.276931e-05 7.452795e-06 \n", - "14 1.694638e-05 1.054243e-05 \n", - "15 1.839608e-05 1.212721e-05 \n", - "16 2.733941e-05 1.874639e-05 \n", - "17 3.128039e-05 2.216655e-05 \n", - "18 4.171634e-05 3.108794e-05 \n", - "19 4.902133e-05 3.881457e-05 \n", - "20 6.542370e-05 5.547309e-05 \n", - "21 7.666464e-05 7.055613e-05 \n", - "22 9.334149e-05 9.261527e-05 \n", - "23 1.121099e-04 1.201568e-04 \n", - "24 1.358183e-04 1.559843e-04 \n", - "25 1.618118e-04 1.980445e-04 \n", - "26 1.910774e-04 2.470661e-04 \n", - "27 2.181750e-04 2.976550e-04 \n", - "28 2.484255e-04 3.522708e-04 \n", - "29 2.671787e-04 4.017728e-04 \n", - "30 2.687268e-04 4.333743e-04 \n", - "31 2.779371e-04 4.798020e-04 \n", - "32 2.712952e-04 5.101453e-04 \n", - "33 2.854872e-04 5.818602e-04 \n", - "34 3.076460e-04 6.656371e-04 \n", - "35 3.189522e-04 7.535916e-04 \n", - "36 4.007223e-04 9.796619e-04 \n", - "37 4.852152e-04 1.233090e-03 \n", - "38 5.752492e-04 1.510840e-03 \n", - "39 6.546092e-04 1.793002e-03 \n", - "40 7.281388e-04 2.078377e-03 \n", - "41 7.933713e-04 2.362686e-03 \n", - "\n", - " Waste_aluminum_Reference.Mod Waste_glass_95-by-35.Adv \\\n", - "0 0.000000e+00 0.001223 \n", - "1 7.984519e-12 0.001124 \n", - "2 4.722305e-10 0.013594 \n", - "3 5.549092e-09 0.013084 \n", - "4 3.735847e-08 0.041959 \n", - "5 1.855378e-07 0.041021 \n", - "6 7.133639e-07 0.086805 \n", - "7 2.236392e-06 0.080150 \n", - "8 6.008380e-06 0.106870 \n", - "9 1.437589e-05 0.109019 \n", - "10 3.149377e-05 0.194930 \n", - "11 6.441594e-05 0.176573 \n", - "12 1.246754e-04 0.186308 \n", - "13 2.303997e-04 0.183984 \n", - "14 4.088584e-04 0.388597 \n", - "15 6.991525e-04 0.386349 \n", - "16 1.154571e-03 0.710085 \n", - "17 1.844060e-03 0.709249 \n", - "18 2.852477e-03 0.884737 \n", - "19 4.279934e-03 0.897174 \n", - "20 6.241336e-03 1.154849 \n", - "21 8.866731e-03 1.169273 \n", - "22 1.229902e-02 1.341590 \n", - "23 1.667688e-02 1.370681 \n", - "24 2.208345e-02 1.644712 \n", - "25 2.845403e-02 1.690431 \n", - "26 3.546196e-02 2.054799 \n", - "27 4.252520e-02 2.107562 \n", - "28 4.947745e-02 2.268901 \n", - "29 5.561036e-02 2.319273 \n", - "30 5.773818e-02 2.359496 \n", - "31 6.241551e-02 2.409155 \n", - "32 6.340878e-02 2.454487 \n", - "33 7.018165e-02 2.545114 \n", - "34 7.877043e-02 2.698586 \n", - "35 8.860343e-02 2.855973 \n", - "36 1.164729e-01 3.268146 \n", - "37 1.483842e-01 3.689297 \n", - "38 1.833232e-01 4.388077 \n", - "39 2.202285e-01 4.959250 \n", - "40 2.585010e-01 5.666096 \n", - "41 2.982925e-01 6.380055 \n", - "\n", - " Waste_silicon_95-by-35.Adv Waste_silver_95-by-35.Adv \\\n", - "0 0.000628 4.285685e-07 \n", - "1 0.000286 2.009061e-07 \n", - "2 0.003413 2.024336e-06 \n", - "3 0.003423 1.558864e-06 \n", - "4 0.010847 3.640602e-06 \n", - "5 0.010312 3.244131e-06 \n", - "6 0.022937 6.211354e-06 \n", - "7 0.020210 5.593786e-06 \n", - "8 0.022679 7.980107e-06 \n", - "9 0.017783 8.192585e-06 \n", - "10 0.029029 1.407329e-05 \n", - "11 0.024710 1.226383e-05 \n", - "12 0.025790 1.265886e-05 \n", - "13 0.025176 1.270379e-05 \n", - "14 0.049139 2.440697e-05 \n", - "15 0.044050 2.524480e-05 \n", - "16 0.077375 4.260399e-05 \n", - "17 0.073905 4.556430e-05 \n", - "18 0.088786 5.748299e-05 \n", - "19 0.086014 6.372512e-05 \n", - "20 0.107003 8.099552e-05 \n", - "21 0.104090 9.108346e-05 \n", - "22 0.118721 1.116072e-04 \n", - "23 0.119929 1.303008e-04 \n", - "24 0.143170 1.605935e-04 \n", - "25 0.145176 1.865898e-04 \n", - "26 0.176042 2.250153e-04 \n", - "27 0.178288 2.524408e-04 \n", - "28 0.190351 2.840788e-04 \n", - "29 0.192278 3.042125e-04 \n", - "30 0.194193 3.059051e-04 \n", - "31 0.195714 3.191703e-04 \n", - "32 0.197314 3.170890e-04 \n", - "33 0.200150 3.413036e-04 \n", - "34 0.206325 3.747291e-04 \n", - "35 0.211471 4.078259e-04 \n", - "36 0.230962 5.196237e-04 \n", - "37 0.247005 6.457871e-04 \n", - "38 0.284028 7.915824e-04 \n", - "39 0.304988 9.393534e-04 \n", - "40 0.334059 1.094822e-03 \n", - "41 0.359673 1.252286e-03 \n", - "\n", - " Waste_copper_95-by-35.Adv ... VirginStock_glass_95-by-35.Adv \\\n", - "0 2.017821e-08 ... 8.025386e+03 \n", - "1 1.964373e-08 ... 7.713247e+03 \n", - "2 2.375177e-07 ... 9.326221e+04 \n", - "3 2.286382e-07 ... 8.976487e+04 \n", - "4 7.628068e-07 ... 3.015900e+05 \n", - "5 7.322020e-07 ... 2.948416e+05 \n", - "6 1.783707e-06 ... 6.551432e+05 \n", - "7 1.897091e-06 ... 6.367041e+05 \n", - "8 3.157919e-06 ... 8.960064e+05 \n", - "9 3.177004e-06 ... 9.134817e+05 \n", - "10 6.275057e-06 ... 1.632957e+06 \n", - "11 6.239381e-06 ... 1.563819e+06 \n", - "12 6.909687e-06 ... 1.646076e+06 \n", - "13 7.413310e-06 ... 1.618034e+06 \n", - "14 1.538162e-05 ... 3.422864e+06 \n", - "15 1.695016e-05 ... 3.382475e+06 \n", - "16 2.977600e-05 ... 6.225730e+06 \n", - "17 3.274601e-05 ... 6.169521e+06 \n", - "18 4.304799e-05 ... 7.656610e+06 \n", - "19 5.026210e-05 ... 7.665399e+06 \n", - "20 6.794914e-05 ... 9.813826e+06 \n", - "21 8.248517e-05 ... 9.753521e+06 \n", - "22 1.077438e-04 ... 1.103668e+07 \n", - "23 1.352746e-04 ... 1.097820e+07 \n", - "24 1.766693e-04 ... 1.301877e+07 \n", - "25 2.190043e-04 ... 1.295812e+07 \n", - "26 2.761605e-04 ... 1.567521e+07 \n", - "27 3.281095e-04 ... 1.561007e+07 \n", - "28 3.858802e-04 ... 1.650213e+07 \n", - "29 4.400673e-04 ... 1.644024e+07 \n", - "30 4.778077e-04 ... 1.653371e+07 \n", - "31 5.374668e-04 ... 1.647720e+07 \n", - "32 5.872917e-04 ... 1.656200e+07 \n", - "33 6.915515e-04 ... 1.650999e+07 \n", - "34 8.210428e-04 ... 1.667850e+07 \n", - "35 9.800827e-04 ... 1.663006e+07 \n", - "36 1.306075e-03 ... 1.718313e+07 \n", - "37 1.696850e-03 ... 1.713665e+07 \n", - "38 2.156751e-03 ... 1.887736e+07 \n", - "39 2.668277e-03 ... 1.882960e+07 \n", - "40 3.231128e-03 ... 1.936071e+07 \n", - "41 3.833597e-03 ... 1.936052e+07 \n", - "\n", - " VirginStock_silicon_95-by-35.Adv VirginStock_silver_95-by-35.Adv \\\n", - "0 1036.039007 23.764247 \n", - "1 640.162987 11.157407 \n", - "2 7689.440032 111.997100 \n", - "3 7557.015035 86.524936 \n", - "4 24561.064097 203.098965 \n", - "5 23449.823276 181.582129 \n", - "6 52589.010236 347.870650 \n", - "7 49194.730224 312.035735 \n", - "8 63216.628971 442.471248 \n", - "9 55380.797254 448.126733 \n", - "10 94182.347348 765.783324 \n", - "11 84123.648855 643.391203 \n", - "12 86301.832023 629.443398 \n", - "13 82645.936987 572.863065 \n", - "14 168191.790193 1135.563865 \n", - "15 158970.397752 1042.433414 \n", - "16 284772.575236 1811.851080 \n", - "17 275074.274371 1695.451381 \n", - "18 333249.422156 1986.370466 \n", - "19 325531.212302 1852.459673 \n", - "20 410069.187218 2219.248128 \n", - "21 400950.699645 2054.572132 \n", - "22 453699.436631 2324.869914 \n", - "23 451295.571449 2312.551905 \n", - "24 535180.337000 2742.398522 \n", - "25 532687.885814 2729.626576 \n", - "26 644383.368470 3301.982294 \n", - "27 641706.263566 3288.264136 \n", - "28 678377.927531 3476.178957 \n", - "29 675834.149848 3463.144001 \n", - "30 679676.958883 3482.835520 \n", - "31 677354.547362 3470.934900 \n", - "32 680840.732666 3488.798989 \n", - "33 678704.010166 3477.849886 \n", - "34 685632.504459 3513.353232 \n", - "35 683642.944913 3503.158229 \n", - "36 706382.592528 3619.681897 \n", - "37 704476.500865 3609.914604 \n", - "38 776039.529623 3976.621545 \n", - "39 774082.643741 3966.593971 \n", - "40 795922.607082 4078.507430 \n", - "41 795922.607082 4078.507430 \n", - "\n", - " VirginStock_copper_95-by-35.Adv VirginStock_aluminum_95-by-35.Adv \\\n", - "0 4.485386 7.738673e+02 \n", - "1 4.371755 7.533675e+02 \n", - "2 52.827243 9.109097e+03 \n", - "3 50.841191 8.767506e+03 \n", - "4 169.758582 2.924719e+04 \n", - "5 162.572253 2.804289e+04 \n", - "6 396.287453 6.340030e+04 \n", - "7 419.025058 6.268391e+04 \n", - "8 694.637934 8.945712e+04 \n", - "9 688.803066 8.572974e+04 \n", - "10 1356.760012 1.523814e+05 \n", - "11 1309.350048 1.476974e+05 \n", - "12 1385.729134 1.554664e+05 \n", - "13 1369.501626 1.528179e+05 \n", - "14 2922.081188 3.232780e+05 \n", - "15 2912.282638 3.194635e+05 \n", - "16 5193.738026 5.879995e+05 \n", - "17 4981.788196 5.826909e+05 \n", - "18 5977.743994 7.231415e+05 \n", - "19 5721.444586 7.178539e+05 \n", - "20 7052.855287 9.190513e+05 \n", - "21 6739.018570 9.134043e+05 \n", - "22 7625.598188 1.033571e+06 \n", - "23 7585.194986 1.028095e+06 \n", - "24 8995.096486 1.219192e+06 \n", - "25 8953.204366 1.213514e+06 \n", - "26 10830.537246 1.467967e+06 \n", - "27 10785.541541 1.461869e+06 \n", - "28 11401.904163 1.545410e+06 \n", - "29 11359.149367 1.539615e+06 \n", - "30 11423.737761 1.548370e+06 \n", - "31 11384.703599 1.543079e+06 \n", - "32 11443.298004 1.551021e+06 \n", - "33 11407.384829 1.546153e+06 \n", - "34 11523.836182 1.561937e+06 \n", - "35 11490.396463 1.557404e+06 \n", - "36 11872.595341 1.609208e+06 \n", - "37 11840.558517 1.604865e+06 \n", - "38 13043.361206 1.767893e+06 \n", - "39 13010.470653 1.763435e+06 \n", - "40 13377.547999 1.813188e+06 \n", - "41 13377.547999 1.813188e+06 \n", - "\n", - " VirginStock_glass_95-by-35_Elec.Adv_DR \\\n", - "0 8.025386e+03 \n", - "1 7.713247e+03 \n", - "2 9.326221e+04 \n", - "3 8.976487e+04 \n", - "4 3.015900e+05 \n", - "5 2.948416e+05 \n", - "6 6.551432e+05 \n", - "7 6.367041e+05 \n", - "8 8.960064e+05 \n", - "9 9.134817e+05 \n", - "10 1.645299e+06 \n", - "11 1.575638e+06 \n", - "12 1.667130e+06 \n", - "13 1.638730e+06 \n", - "14 3.145681e+06 \n", - "15 3.108563e+06 \n", - "16 6.343689e+06 \n", - "17 6.286415e+06 \n", - "18 8.816110e+06 \n", - "19 8.826231e+06 \n", - "20 1.207558e+07 \n", - "21 1.200138e+07 \n", - "22 1.379303e+07 \n", - "23 1.371994e+07 \n", - "24 1.791196e+07 \n", - "25 1.782853e+07 \n", - "26 2.149434e+07 \n", - "27 2.140503e+07 \n", - "28 2.374752e+07 \n", - "29 2.365846e+07 \n", - "30 2.395775e+07 \n", - "31 2.387587e+07 \n", - "32 2.514029e+07 \n", - "33 2.506137e+07 \n", - "34 2.656348e+07 \n", - "35 2.648635e+07 \n", - "36 2.901691e+07 \n", - "37 2.893849e+07 \n", - "38 3.022863e+07 \n", - "39 3.015224e+07 \n", - "40 3.221011e+07 \n", - "41 3.220988e+07 \n", - "\n", - " VirginStock_silicon_95-by-35_Elec.Adv_DR \\\n", - "0 1.036039e+03 \n", - "1 6.401630e+02 \n", - "2 7.689440e+03 \n", - "3 7.557015e+03 \n", - "4 2.456106e+04 \n", - "5 2.344982e+04 \n", - "6 5.258901e+04 \n", - "7 4.919473e+04 \n", - "8 6.321663e+04 \n", - "9 5.538080e+04 \n", - "10 9.489416e+04 \n", - "11 8.475944e+04 \n", - "12 8.740568e+04 \n", - "13 8.370303e+04 \n", - "14 1.545716e+05 \n", - "15 1.460970e+05 \n", - "16 2.901681e+05 \n", - "17 2.802861e+05 \n", - "18 3.837159e+05 \n", - "19 3.748289e+05 \n", - "20 5.045760e+05 \n", - "21 4.933560e+05 \n", - "22 5.670082e+05 \n", - "23 5.640039e+05 \n", - "24 7.363309e+05 \n", - "25 7.329017e+05 \n", - "26 8.835975e+05 \n", - "27 8.799265e+05 \n", - "28 9.762230e+05 \n", - "29 9.725623e+05 \n", - "30 9.848659e+05 \n", - "31 9.815007e+05 \n", - "32 1.033479e+06 \n", - "33 1.030236e+06 \n", - "34 1.091986e+06 \n", - "35 1.088818e+06 \n", - "36 1.192849e+06 \n", - "37 1.189630e+06 \n", - "38 1.242671e+06 \n", - "39 1.239538e+06 \n", - "40 1.324142e+06 \n", - "41 1.324142e+06 \n", - "\n", - " VirginStock_silver_95-by-35_Elec.Adv_DR \\\n", - "0 23.764247 \n", - "1 11.157407 \n", - "2 111.997100 \n", - "3 86.524936 \n", - "4 203.098965 \n", - "5 181.582129 \n", - "6 347.870650 \n", - "7 312.035735 \n", - "8 442.471248 \n", - "9 448.126733 \n", - "10 771.571005 \n", - "11 648.253863 \n", - "12 637.494346 \n", - "13 580.190318 \n", - "14 1043.606109 \n", - "15 958.017345 \n", - "16 1846.180110 \n", - "17 1727.574992 \n", - "18 2287.181799 \n", - "19 2132.991866 \n", - "20 2730.708346 \n", - "21 2528.080207 \n", - "22 2905.492307 \n", - "23 2890.097949 \n", - "24 3773.144631 \n", - "25 3755.572276 \n", - "26 4527.775439 \n", - "27 4508.964696 \n", - "28 5002.411869 \n", - "29 4983.653854 \n", - "30 5046.700552 \n", - "31 5029.456308 \n", - "32 5295.807447 \n", - "33 5279.187303 \n", - "34 5595.612704 \n", - "35 5579.375427 \n", - "36 6112.456154 \n", - "37 6095.962398 \n", - "38 6367.760775 \n", - "39 6351.703628 \n", - "40 6785.234379 \n", - "41 6785.234379 \n", - "\n", - " VirginStock_copper_95-by-35_Elec.Adv_DR \\\n", - "0 4.485386 \n", - "1 4.371755 \n", - "2 52.827243 \n", - "3 50.841191 \n", - "4 169.758582 \n", - "5 162.572253 \n", - "6 396.287453 \n", - "7 419.025058 \n", - "8 694.637934 \n", - "9 688.803066 \n", - "10 1367.014212 \n", - "11 1319.245930 \n", - "12 1403.453419 \n", - "13 1387.018351 \n", - "14 2685.451582 \n", - "15 2676.446517 \n", - "16 5292.143456 \n", - "17 5076.177826 \n", - "18 6882.999668 \n", - "19 6587.886873 \n", - "20 8678.295388 \n", - "21 8292.130123 \n", - "22 9530.045847 \n", - "23 9479.552186 \n", - "24 12375.954750 \n", - "25 12318.317238 \n", - "26 14851.151872 \n", - "27 14789.452434 \n", - "28 16407.964440 \n", - "29 16346.437949 \n", - "30 16553.231795 \n", - "31 16496.670489 \n", - "32 17370.305075 \n", - "33 17315.790826 \n", - "34 18353.669524 \n", - "35 18300.411083 \n", - "36 20048.921568 \n", - "37 19994.821874 \n", - "38 20886.323457 \n", - "39 20833.655842 \n", - "40 22255.641342 \n", - "41 22255.641342 \n", - "\n", - " VirginStock_aluminum_95-by-35_Elec.Adv_DR \n", - "0 7.738673e+02 \n", - "1 7.533675e+02 \n", - "2 9.109097e+03 \n", - "3 8.767506e+03 \n", - "4 2.924719e+04 \n", - "5 2.804289e+04 \n", - "6 6.340030e+04 \n", - "7 6.268391e+04 \n", - "8 8.945712e+04 \n", - "9 8.572974e+04 \n", - "10 1.535330e+05 \n", - "11 1.488137e+05 \n", - "12 1.574549e+05 \n", - "13 1.547726e+05 \n", - "14 2.970990e+05 \n", - "15 2.935934e+05 \n", - "16 5.991403e+05 \n", - "17 5.937311e+05 \n", - "18 8.326524e+05 \n", - "19 8.265640e+05 \n", - "20 1.130861e+06 \n", - "21 1.123913e+06 \n", - "22 1.291699e+06 \n", - "23 1.284855e+06 \n", - "24 1.677433e+06 \n", - "25 1.669621e+06 \n", - "26 2.012920e+06 \n", - "27 2.004557e+06 \n", - "28 2.223930e+06 \n", - "29 2.215591e+06 \n", - "30 2.243619e+06 \n", - "31 2.235953e+06 \n", - "32 2.354365e+06 \n", - "33 2.346976e+06 \n", - "34 2.487650e+06 \n", - "35 2.480432e+06 \n", - "36 2.717424e+06 \n", - "37 2.710091e+06 \n", - "38 2.830925e+06 \n", - "39 2.823787e+06 \n", - "40 3.016522e+06 \n", - "41 3.016522e+06 \n", - "\n", - "[42 rows x 31 columns]" - ] - }, - "execution_count": 19, - "metadata": {}, - "output_type": "execute_result" - } - ], - "source": [ - "USyearly" - ] - }, - { - "cell_type": "code", - "execution_count": 21, - "metadata": {}, - "outputs": [ - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVAAAAE8CAYAAACW3un4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOzdeVxV1fr48c+ZOAdkdsLZ1JxSRAQFciCjq9nNsrRyuGp1tUkbtG5eLafq5q2sX19vWVbWtbRuo0ODGs6KE5UiDuSEqKjIPHOm9fuDOoogIpzDAX3erxcv2Puss9ZztvC4915rr6VRSimEEEJcNa27AxBCiPpKEqgQQlSTJFAhhKgmSaBCCFFNkkCFEKKaJIEKIUQ1SQIVFerfvz/bt28vs2/EiBEsX76cv//97yQnJ1e5rtTUVMaOHVvl8nl5eQwfPpzg4GC6dOnChAkTANizZw+xsbFVrudi48ePZ/Xq1ZWWadu2LQ8++KBje/bs2bz33nvVau9i06ZN45NPPqlxPaLukQQqKnTffffx1VdfObaLiorYunUrgwcP5sMPP6Rt27Zlylut1svW1bx5c5YsWVLlthcsWEC3bt1ISEjg4MGDPP3000DNEmhVbdq0iVOnTrm0DXHtkAQqKnTvvffy3XffObZ//PFHBg4ciMlkIjo6mkOHDpGcnExoaCgPPvggvXr1wmaz8fe//52uXbsyYsQIevbsSXJyMsnJyURERAClZ3UTJkygb9++tG/fnrVr15Zr+8yZMzRr1syxfdNNNwEwc+ZMPv74Y0JCQtiwYQNHjhyhX79+BAcHc99991FQUABAXFwcvXv3pkePHtx5551l6rbb7YwfP57XX3+9ws/9xBNP8NZbb5Xbn5SURExMDL169WLIkCFkZGQAsH37dvr27UtoaCgjR46kuLgYgIULF9KxY0eio6M5duwYADabjb/97W/cdNNNdO/enWXLllXtH0PUXUqIy+jfv7/asWOHUkqp+++/X3333XdKKaUGDBigDh48qI4fP650Op1KTExUSin1v//9T91///1KKaX27t2rtFqtOn78uDp+/Ljq06ePUkqpWbNmqdtuu01ZrVYVHx+voqKiyrW7a9cu1bBhQzVw4ED16quvqrS0NKWUUh9//LF6/vnnHeVuv/12tXz5cqWUUs8995x6+eWXVXFxsWrfvr06cOCAUkqpjIwMpZRS48aNUz/88IMaP368evXVVyv8vG3atFGZmZmqXbt2KjMzU82aNUstXLhQKaVUTEyMSklJUUop9d///ldNnTpVlZSUqH79+qmsrCyllFJz5sxRCxYsUKdOnVIdO3ZUubm5KicnR91www3q448/VvHx8ap///6O9rKzs6/uH0TUOfXiDDQ1NZXQ0FBMJtNlLxVXr15NdHQ00dHRNGvWjOXLl9dylNeeESNG8PXXX1NcXMyWLVsYPHhwuTJdunRxnCHGxcVx7733AhAcHEznzp0rrPeOO+5Ap9PRs2dPTpw4Ue718PBwDh8+zEMPPcTOnTuJiIhwnNldbM+ePdx1110AjBkzhq1bt5KUlESHDh3o0qULAIGBgY7yzz//PC1btmTatGmX/cyenp6MHz+ed99917EvLy+PuLg47rzzTkJCQnj99dc5efIkSUlJ7Nu3j+joaEJCQvjiiy9ITk5m9+7dxMTE4OPjg6+vL7fffjsA7dq1IyUlhaeeeor169fj5+d32ThE/aB3dwBVERgYyLp16xg2bNhlywwePNjxB96nTx9iYmJqK7xr1vDhw+nbty9RUVFER0djMpnKlWnQoIHjZ6UUGo2mzHZFjEYjAFqt9rL/IQYEBDB69GhGjx5Nt27dSExMLFemorYujeFiffr0YcuWLRQVFeHp6VlhGSi9jA8PD2fEiBGOOlu0aMGePXvKlEtISCA8PLzcbYiLb31cHFtAQAAJCQn8+OOP/Otf/2Lz5s3Mnj37snGIuq9enIGaTCYCAgIc20opHnvsMQYOHMgdd9xBVlaW47Vjx47RtGlTvL293RHqNSUoKIgWLVowa9Ys7rvvviuWj4qK4ttvvwVg3759JCUlVavduLg4cnNzATh//jwZGRk0b94cb29v8vLyHOV69OjBqlWrAPj888/p27cvnTt35vDhw462MzMzHeWHDx/O/fffzwMPPIDNZrts+4GBgQwdOpSlS5cC4OvrS0BAgCNRlpSUkJSUROfOnTl+/Dj79u0DSs9Ujx8/Tu/evYmNjSUvL4+8vDxH7396ejpKKe6//36mT59eLiGL+qdeJNBLff/997Ru3Zr169czadKkMkNNvv3220rPVMXVue+++0hOTmbQoEFXLDt8+HCMRiPdu3dn/vz5dO3aFV9f36tu8/Dhw0RFRREcHMwtt9zCSy+9RPPmzYmOjmbHjh2OTqT/+7//47XXXiM4OJhjx47x9NNPYzQa+e9//8vo0aPp0aMH48ePL1P3Y489Rs+ePZk4cWKlMUyZMoVz5845tj/77DP+/e9/06NHD3r16sWePXvw8PBg2bJlPPLII/To0YN+/fpx4sQJWrRowdNPP01YWBjDhw/n5ptvBuDUqVMMGDCAkJAQ/vGPf1R6K0HUDxp1ueusOig6OprY2FjefPNNli5dSkBAAFarlcjISEev6oABA/j2229p2LChm6O9PhUUFNCgQQMOHz7M0KFDOXjwoLtDEsJl6sU90Et16tSJsWPHMnXqVAAsFgsAZ8+excPDQ5KnG8XExFBYWAjAO++84+ZohHCtenEJb7FYiImJYe/evQwaNIigoCCSk5MZOHAgAwcO5KeffgJgxYoVjl5Z4R7bt29n79697N27l4EDB7o7HCFcyiWX8Dt37uSZZ55Bp9MRFhZWZmByamoqY8aMobi4mLlz50pvuRCi3nLJGWibNm1Yv349W7ZsIS0tzdFLCTBv3jxefvll1q5dy8svv+yK5oUQola4JIEGBQU5xgzq9Xp0Op3jtYSEBCIjI/H29sbHx6fMsBQhhKhPXNqJlJCQQHp6Ol27dnXss9lsjoHOfn5+ZGVl4ePjU+Z9ixYtYtGiRQDs37+f1q1buzJMN6s3gyCEuKZkZGSSnp5eozpclkAzMzOZNGkSX375ZZn9F5+N5ubm4u/vX+69EydOdIzTCw0NZdOmTa4K0+18fOQMXAh3CAsbWuM6XHIJb7VaGTNmDK+//jpBQUFlXgsODmb79u0UFBSQm5tbrYHWQghRF7gkgX711Vfs3r2b559/nujoaLZv387kyZMB+Mc//sGMGTOIiYlh+vTprmheCCFqRZ1/Ekku4YUQrhAWNpT4+Pga1VEvBtILIURdJAlUCCGqSRKoEEJUkyRQIYSoJkmgQghRTZJAhRCimiSBCiFENUkCFUKIapIEKoQQ1SQJVAghqkkSqBBCVJMkUCGEqCZJoEIIUU2SQIUQopokgQohRDVJAhVCiGqSBCqEENUkCVQIIapJEqgQQlSTJFAhhKgmSaBCCFFNLkmgqamphIaGYjKZsFqtZV4bP348ffr0ITo6mmXLlrmieSGEqBV6V1QaGBjIunXrGDZsWIWvL126lA4dOriiaSGEqDUuOQM1mUwEBARU+JpGo2Hs2LHceeednDhxwhXNCyFErXDJGWhl5s+fT2BgIFu3bmXq1Kl8/fXXtR2CEEI4Ra0n0MDAQAD69u3LtGnTKiyzaNEiFi1aBJTeT929e3etxVfbtFqbu0MQQlRTrSfQ3NxcfH19SUpKwt/fv8IyEydOZOLEiQCEhoYSHh5emyHWKh+fPHeHIISoJpckUIvFwu23387evXsZNGgQM2fOZOvWrcyYMYPRo0eTlZWFRqNh4cKFrmheCCFqhUsSqMFgIDY2tsy+AQMGALBq1SpXNCmEELWu1i/hhRDC1Y5kHuH9+Pf5PeP3Cl9v49/GKe1IAhVCXDN2nNrBwviFxB6NRVVSLs+cB3jXuD1JoEKIes2u7Px4+EcWxi/ktzN7arVtSaBCiDptzZE1vLP7HYw6Ix56D4w6Y+mX3ohBa2DbyW0kZ7vnoZxKE2j//v0r3K+UokWLFnzxxRcuCUoIIQBS81KZsnYKWUXZ7g6lQpUmUD8/v8v2mt95550uCUgIIaD00vzJn56ss8kTrvAsfGVDjmQ4khDClRbsXEDcye3uDqNSVboHOm/ePKZNm8bq1auZOXMmY8aM4cknn3R1bEKI69Qvqb8wf/v8ar1XA/Ty64H3QRN5x3MqLNO8WUt+IasGEZaqUgL9+eefmTZtGkuXLmXbtm1ERUVJAhVCuEReSR5P/PgEVvvVzRPhb/Slp6Y7GWvTiN+9t9KyBZ0VtTaMqbi4mCVLltCkSRMMBgOenp41blgIISryfOzzpOScdGzrtToCjQFYlIUSmxmzzVwmuXb0bU/rs81I+CSRDRnbajXWKiXQJUuWsGnTJubMmUNJSQmTJk1ydVxCiOvQ/xL/x/JDK8rs61sQwcbXSxOjUWfAz8OEp9EDg1GPwWTgSMpRfldHy7znxm7dCBz0FwqaNq2wnRa+vuyb+2GN461SAvX09CQ/P5+33noLpSob3y+EENVzLPMYL2x4ocy+m717s3H+hbPKEpuFkiILFBWUe3+DBl50GzSIjOhoDjdvXmlbVi8vp8RcpQR61113MWXKFHr06OGURoUQ159z+eeIT43HQ+fhGBTvofUo/a7z4MmfnqTAXOgo39G3A3tfT0Cv1xHYqBEWsxmL2YzZbMZstjjKdbipKw0HDyaxVy92enjU6meqUgJt27YtI0eOdHUsQohr1Oojq3l27bNkFlWt5zvA5I/1qxLy84vp/Y/n2NW7d7kyBqsVg93OEQ8Pjjg74CqqUgLNzMykV69edO3aFShd12jJkiUuDUwIUf8VmguZuWEmyxI/r/J7tBoNHQ+0Y2fir0SMfIAdFSRPAItej6XCV2pPlRLo4sWLXR2HEOIas+fMHib9NIljWcev6n397VFs/G4b3aMi2XX33S6KzjmqlEADAwNZuHAhR48epX379jz66KOujksIUU/ZlZ0FOxcwf/v8qx7LGe7bi82zt9GybRtOPP44dp3ORVE6R5US6OjRo3nggQe46667iI+PZ9SoUXz//feujk0IUQcdyTzCx799jFFvdHQCGXVGPHSlnUHLDy1n5+ldV11vW5/WHPm/wzTw9kb3/PPkmkwuiN65qpRAc3JyGDVqFACdOnXi/fffd2lQQoi66XzBecZ8O6bMQPeqau9zAy2tzVAahV1jx661Y9cobNiwa+zkr8onJyuP7i+9xN7GjV0QvfNVmkAXLFjA5MmTCQkJYcKECYSGhhIfH09wcHBtxSeEqCOKLEWMXT72qpOnh87AzaoPO/+9m6NFld8PjXj8MXZ06lSTMGtVpQn0u+++Y/Lkybz99tvEx8dz7NgxHnvsMcLCwmorPiFEHWCz23j0h0fZezbhqt7Xxa8jup80bNi+9Yplw+8Ywo6BA6sboltUmkCPHTvGzJkzy+xLTExk5cqVzJ0716WBCSHqjhnrZvDz0dgrF/yDp95EZHEYW1/ajtliw2Qy0nPcODK6dEFns6GzWtHabOhsNrQWC1q7nYRu3Vz4CVyj0gQaEBBATEzMVT++mZqayl//+lcOHDhAfn4+ev2FZhITE3n00UdRSrFw4UK5HSBEHfefXf9hScKnZfYFGP0w2y2YbWYsdmuZ14L9ulH8TQHr95aedXaPiiRz/Hi2BwbWWsy1pdIE6u/vf9llPSoTGBjIunXrGDZsWLnXXnzxRT7//HO0Wi2PP/44K1asqKAGIURdsPzgcl7d8mqZfQMt/Vg/fwsawMtgxM9gxOjpgc6kw2D0IDE5EbsdGjZuxA2PPEJ8SIh7gq8FlSbQOXPmVKtSk8mE6TJDEDIzM2nVqhVQ2rsvhKib4k7G8fSap8ssD3wLfVn/f1sAUECBpYQCSwkUXpjcQ6OB3vcM49A99xBfD4Yi1USlCXTevHn069ev3CW8UoqhQ4fyww8/XHWDdru9wp8vtmjRIhYtWgSU3g7YvXv3VbdTX2i1VzfQWAhnUCjSCtLILclFgwaNRuP4rtVo0KAlz5zHlGYXzj59NN7knc2n/+wH0Pv7lWZQpS58Ufrd5ulFsdGDgUVAkbs+YeVMOVYWOaGeShNoUVERt956a7n9SimCgoKq1aBWq63w54tNnDiRiRMnAhAaGkp4eHi12qoPfHzy3B2CuM7EHo3lxY0vciI7pcrv6e95M1te2Ub7Ll04/89/kuN5pTPLun1i0NnL6JR6Kk2gGzZscEojFwsMDOTUqVNotVr8/PycXr8QomIp2SnM3DiTtUd/vqr39W0QwZaXt9Guc2fSp00jR1akcKjSk0hXy2KxcPvtt7N3714GDRrEzJkz2bp1KzNmzGDOnDk88MADKKV45513XNG8EOIixdZi3tn1Du/sfodia8lVvTfCJ5ztL+2kXafOZPzzn2Q7aSLia4VLEqjBYCA2tuyYsQEDBgAQHBzM1q1XHlQrhKiaJXuX8OneTzHoDBi0htJn1P94Lt2oM7Ln3J6rulz/U7hvKL+88ittb+woyfMyqpRAP/nkE8aPH8/u3buZNWsWDz30EMOHD3d1bEKISqQXpjN17dSrGuD+p5YNmnNDckuS1x5Fb9SjM+rRexjQm/TojDq0eh2JuxJp3a4DmdOnS/K8jCol0M8++4zx48fzn//8h48++og777xTEqgQbhR7NJYpa6eQXphxVe/zN/kRkteNXW/Ec6owtdKy7Tp1JHP6dLIkeV5WlRJofn4+mzdvxtvbm2bNmtGgQQNXxyWEqECRpYg5G+eUezLoSow6D/powzj40UE2nt1GgwZeRD70INk33YTWbkdrtZY+Vmm1orVa0dhsHO3cWZLnFVQpgb711lv8+OOPvPjii5SUlHDPPfe4Oi4hrksZhRnoNDrHfUyd9sKEwglnE5j00ySOZJZdwtfTYKKzZ0ds2LFr7FiVBauyY8OKxW6hlaYFad+cY3NiHB4eBiJGPsDvgwezXU6EaqzK84G+9NJLHDlyhJdeeon777/f1XEJcV3JL8ln+vrpfH3gmzL7dRqtY8LifHN+uRneg/26UfhVLr/tu/wsSadIRafT0nvonZwcOpQd/v4u+QzXoyol0DfeeIPBgwfzr3/9iwkTJvDkk09e008HCVGbfjvzG0/8+ATJ2SfKvWZTdoosxRRZisvs9zP60CO9G1vf2oHBYOCm8DBsVis2iwWbxYL1j+8Wi4WW3bqRce+97KonkxTXJ1VKoHl5eaSkpKDT6YiMjJR7oEI4gVKKd3e/y2vbXis3o1Flwn17cWbJaTYf3U73yAgyH3yQ/ZXMdJTmjGBFhaqUQKdNm8aMGTN48cUXKS4uJiIiwtVxCXFNO5d/jid/epItKVUfE93IsyGdj7dn6/xdNGzciPAXZrD7Gp7pqD7QqCpO9llcXMz58+cdE4u0bt3apYH9KTQ0lE2bNtVKW+4gz8Jfe/JL8jmQfgC9Rl9mwTWDzoBRZyQ+NZ6pa6eSWZTleI9R50GkJZwj/0tCZ9Ch+3NsplGPzkOH3sPAyYOnyUzLpvc993Bg2DDyrvGZjlyps5cXuX+fSXx8fI3qqdIZ6Lx581i7di2HDh2iffv2GI3Gck8aCXE9yy/J5+djP7MyaSWbTmy6qkcmu/h1Qn1vZeOubQQ2DMRaYMZiNmOxWLBaL3Qate/ciUb/nMXOli1d8RFENVQpga5cuZK4uDiio6PZuHGj9MKL64JSitS8VPRaPUadEaO+9EurKZ1FLL8kn7VH17Lq91VXnTQBTHojESWly14ENmpCz5fm8luXLpcGgYfNhsFm47heX+fXSb/eVCmBenh4AODl5cXmzZs5cOCAS4MSwl1yi3PZdGIT65PXs/7Yes4Xppcro9fq8NB5YLFZrqrz52Jd/TpjX2lh0y/b6DNiBIl33cVvxgqmWNNoMOv1mPUumbZC1FCV/lXefvttSkpKmD9/PgsXLuSNN95wdVxCOE1qXioH0g6UnknqjWXOJk06E/mWfDYc38C64+vYfXr3FZOi1W7Daq94pmCT3khXr85453uDVoFWobRg1yjsWjtKY8eYaWTrnDha3tCOzm+8wY42bVzxsUUtqDSBpqSUzuASEBDAuXPnaNCgAc8++2ytBCZEdZ3NP0tcShxxJ+OIOxVX4fhKZwo0BdBZdyPqsJ1Dqw/xa9beSst7eBjo/eBDxMfEYJEzy3qt0n+9tm3bEh4ezk033eTYp5RCo9GwePFilwcnxMUsNgvjV4wnpzgHT70nngbPMt/NNjO7UndxNPNYjdtq7NmI9rob8Mg2gFZTejap04BOobSlZ5M6i47833LZt+4AcdZd+Ph40zEyEm2vXmS2aIHWbkdnsaC75Fnz9KAgtjdq5IQjItyt0gS6c+dOli9fzv79++nSpQv33nsvYWFhtRWbEGV88OsHbDi+0SV1exk86ezZEe/0BmTGp3Ngx++cV+Xvf16qTft2hD0wkrwePTjUujW/SCfPdaXSBBoeHu5Yj2jnzp1MmDCBsLAwPvjgg1oJTog/nc49zZvb36xyeU+9ifYN2hGY70/xsWIyD6Wj1WvRG/XoTQZ0Rh16Dy0aox5rgYVDm37n16K9GI0e3NijB30mTCC7UydOBgWhs9vRW6142Gzo/ugR11mtFHp5ccLfH9feIBB1WaUJND09nRUrVrBu3TpMJhMzZsxgyJAhtRWbEA4vbniRQksRfU0RaM1alB6Uzo5da8emtWPT2NAoLZ6ZRnL355AUd5jEogNotRra3ngjgW07o+x2bGYzNouFkgIzhWYzVkshHl5e9HhgJHmdO3O4TRsS5b6kqKJKf1OaNWtGz549ueWWW/Dw8CAhIYGEhNJZX+bOnVsrAQrx89GfWX1kDeE+vYibuxO7/fIPzwUEBtCmWzdC/jaOvPbtOd6qFcc8PKj5XVEhyqs0gcrTRsLdiixFvLDhBQI9A8j6Nge/Tz4hx2jEaLdjsFrx+OPLYLFg1es5HRhI1pWrFcIpKk2gfy4EJ4S7vL3jbU7mnKLPmTCK77+bI3/MkF6k01FkMLg5OnG907qq4meeeYZ+/frx1FNPldk/fvx4+vTpQ3R0NMuWLXNV8+IacCTzCO/98h7hPr2gpDV7O3Vyd0hClOGSBPrrr79SUFDAli1bMJvN5SZfXrp0KRs3bmTUqFGuaF5cI/4Z+098PHwo/MlMwr33ujscIcpxSQLdvn07MTExAMTExLBjxw7HaxqNhrFjx3LnnXdy4oQMABEV++bAN2w7GceNx9tjHfkQRX/MxyBEXeKSBJqdnY2vry8Afn5+ZGVduK0/f/584uLieP7555k6daormhf1XE5xDnM3zSXcJwyNsRtJMn2bqKNcMuDN39+f3NxcAHJzc/G/aBGrwD+WHujbty/Tpk2r8P2LFi1i0aJFAKSmpl7T6y9ptbYrF7oGZRVnYbNbAU251/LN+TzSaAYe+SaKIpoyQLrVhZOZcqwsckI9LkmgkZGRvP/++9x3333ExsYyfvx4x2u5ubn4+vqSlJRUJrFebOLEiUycOBEonZH+z6ehrkXX24z0Simmr5vOf/cvqbRcxPkIUqPGkRJ4ff4HI1yrs1cFUwdWg0sSaGhoKCaTiX79+tGjRw9at27NK6+8wowZMxg9ejRZWVloNBoWLlzoiuZFHWWz25i6dipf7v8KL4MnHtqK72t20ncEWx9SZMINUce57Jm1t99+u8z2jBkzAFi1apWrmhR1mMVmYfJPk1mZtIpQ3xCS3koiOzunwrKFUXr2PXNzLUcoxNWTh36Fy5VYS3jk+0dYe/RnojwjyPneQoO33sFuMqE0GhRc+NJo2KfVgqb8vVEh6hpJoMKliixFPLjiQbac2EI/WxTFewI4PHWMLFEhrgnyWyxcJq8kj7999zf2nN1DRHpvrMab2D3+L+4OSwinkQQqqi01L5W7vriLjMIMx/pCF9YcMpFVnEleSR7BSd0xhwzml+7d3R2yEE4lCVRUS3ZRNqO+GYXZYqahsSEWuwWzzUyeOY8Sawk2Zad5g2a0i+9Azh1jOdyihbtDFsLpJIGKq1ZkKWLc8nF4WD1o+EMQGo0B6x8TFVutVuxmC3oNeDdtwpnHHufcH0+lCXGtkQQqrorNbuPxHx6nsKAAn70d2D1tNDZZB0hcp1w2nZ24Nk2LnUZK+gn8DnZh132SPMX1Tc5ARZW9se0N4lPiaXw4hG3D7kVp5f9fcX2TBCqqZMneJaw6sIqglD5sGXqXDHQXAkmggtJOoe8OfUeRpQgPnQdGnREPfel3o87IydyTfBz/MUGpUWy5/U53hytEnSEJ9DpWYi3hs4TPWLBrAb46XwzosWLFpmxYlBWr3YpZmWlkaEjQ+f5svm2wu0MWok6RBHodstqtfJH4BW/vfJumNKXjse7YArpi9/BAY7OBxVL63WoFmw17YCCboyPdHbYQdY4k0GtMgbkArUaLUW9EqynbyWNXdr458A0Ldi0gsCSA9imh5NwYwda7Oso9TSGqQRJoHbT79G7+sfYfZJdk42nwxN/kj7/JnwDPAPyN/vgafckz55FZlElGYQYZRRlkFWVhtpkxao3YlR2LMqMADRo0Gg0GrQGtRktra2uap4WR2i2KI8HN3f1RhajXJIHWIVa7lbe2v8XOfTtpcjqMQKsendWMzsMGeitZehvZurMo7Wk0So/G5oHO1oxGttYE2D2wmTxRXl6gVJlLcY3NhhY7ymDgWGgYZ0IqXglACHF1JIHWEcezjjPr51nYj3tzttMYjgU3cXdIQogrkARaB3y+73N+2LGGQtWHHX1D5H6kEPWEJFAXsys7dmVHKYVCATh+NtvMvLxxGqkn9fzWcQxZXl5ujlYIcTUkgTqZzW4jNT+VlOwU0rPSyS8sRKGBPzpzSpVuazV6finpyYGQ1m6MWAhRXZJAq8Bmt3Es6xiFlkK0Wi06jQ6tRotWU/qzRqPhfN55MjIzyDWbsSsvivEh27cjeX6eldZ9wCuplj6FEMLZJIFehlKKUzmnOHYmmbSiIgqs3liV7o/bk3Y0gFZTuhSaBij08CbDtzM2mWBDiOuGyxLoM888Q3x8PKGhoWWWOE5MTOTRRywfQiYAACAASURBVB9FKcXChQsJDg52VQjVklGQwZGzKZwqKiDD5sk5nxbYPSUpCiHKc0kC/fXXXykoKGDLli089thj7N69m/DwcABefPFFPv/8c7RaLY8//jgrVqxwRQgOFpuFsznnOV9USLrZTKa5BKuyA3/eiSyl+WN53Uy7gQzPAPAOcGlcQoj6zyUJdPv27cTExAAQExPDjh07HAk0MzOTVq1aAZCTk1Ol+mx2G4VWCwVWK4U2G4U2K0U2G3alypTTaOCPjm6yLWbSSorJMFvJNphAoy/9MkpPtxDCOVySQLOzs2nfvj0Afn5+7N+/3/Ga3W6v8OeLLVq0iEWLFgFw4MABYm6NcUWYV+38+fM0btzYqXVW92FKV8RSn+MAieVyJJbycoFDhw7VuB6XJFB/f39yc3MByM3Nxd//wqOD2os6WbSX6XCZOHEiEydOBCAsLIz4+HhXhHnVJJa6GwdILJcjsVQsLCysxnW4pHckMjKSdevWARAbG0tERITjtcDAQE6dOkVqaip+fn6uaF4IIWqFSxJoaGgoJpOJfv36odVqad26Na+88goAc+bM4YEHHmDEiBHMmTPHFc0LIUSt0M2ePXu2Kyq+/fbbeeihhxgyZAje3t70798fgKZNm/LQQw/x8MMPExQUVKW6evXq5YoQq0ViKa+uxAESy+VILBWraSwapS7pyhZCCFElMkJcCCGqyW0JdOfOnURFRdGvXz+eeeYZAF5//XX69u3L6NGjsVgsWCwWIiMj8fb25siRI473XlrOnbH4+fkRHR1NdHQ0mZmZtRLL8ePH6devH/3792fUqFHYbDYAli5dSlRUFH/9618doyDcFUunTp0cx+XAgQO1Ekt6ejpRUVEMGDCAoUOHUlRUVGE5d8bijt+XP33zzTeOMdiVlavtONx1TCpq96r/hpSbnDlzRhUVFSmllBo1apTavHmzuv3225VSSs2bN099+eWXym63q7Nnz6px48apw4cPK6WUSktLK1fOXbEopdTNN99c4/avNpbMzEyVlZWllFJq+vTpauXKlcpsNqu+ffsqi8WivvjiC/Xaa6+5LRal3HNcrFarstlsSimlZs+erb788ku3/b5UFItS7jkuf3rggQdUZGSkUsr5f0fVjUMp9x2TS9utzt+Q285Ag4KCMJlMAOj1ehISEoiOjgYuPL2k0Who2rRpmfft2rWrXDl3xQJw8OBB+vXrx7Rp01BOuJ1clVgCAgIcY2v1ej06nY7ff/+d7t27o9fra/W4VBQLlD5x1r9/fx555BGKi4trJRadTucYW2yz2bjxxhvd9vtSUSzgnt8XgB9++IHbbrvNEZOzj0t14wD3HZNL263O35Db74EmJCSQnp6Ov78/vr6+QOmpdVZWVoXls7Ozq1SuNmIBOHz4MJs3byYrK4tVq1bVaiypqanExsbyl7/8xe3H5eJYALZu3crmzZtp06aN46my2ohl165dhIWFsX79em644Qa3HpdLYwH3/b7897//ZcyYMY7yrjouVxsHuO+YXNpudY6JWxNoZmYmkyZN4qOPPqr06aWLVbVcbcQCpQ8GaDQa7r77bhITE2stlpKSEsaNG8cHH3yAXq9363G5NBYoPS4Aw4YNq9Xj0rt3b+Lj4xk2bBiLFy9263G5NBZwz+/L+vXriYqKwsPDw/EeVxyX6sQB7vsburTd6hwTtyVQq9XKmDFjeP311wkKCiI8PJxNmzYB5Z9eulhVy9VGLAUFBY5Ok23btjme/6+NWCZOnMjjjz9O165dAejYsSOJiYnYbLZaPy6XxmI2mykpKQFq97iYzWZHeV9fXzw9Pd32+1JRLO76fUlMTGTlypUMHjyY/fv388ILLzj9uFQ3Dncdk4rardbfkBPv3V6VZcuWqUaNGqkBAwaoAQMGqLi4ODVv3jx18803q5EjR6qSkhKllFIjRoxQzZo1U1FRUWr58uVKKVVhOXfE8ttvv6mePXuqvn37qrFjxyqr1VorscTFxSlvb29HmW+//VYppdSSJUtUZGSkGjJkiMrOznZbLGfPnlU9e/ZU/fr1U0OHDlW5ubm1EsvOnTtV//79VXR0tLrnnntUQUGBUso9vy8VxeKu35eLXdxx4szjUt043HVMLtfu1f4NyUB6IYSoJrd3IgkhRH0lCVQIIapJEqgQQlSTJFAhhKgmSaBCCFFNkkCFEKKaJIEKIUQ1SQIVdc6fj+D9qWXLlmW2W7VqRXZ2dpXqSk5O5uuvv3Z6jEKAJFBRB0VERDhmwjl9+jSNGzcus+3t7V3lZ7clgQpXkgQq6pzIyEh27twJwI4dO/jb3/5GUlKSYzsiIoKBAwcSGhpKSEgIGzZsACAxMZFevXoREhJCz549yc7OZubMmaxZs4aQkBCWLl1KXl4eo0ePJjw8nD59+vDLL7+47XOKa0CNHzwVwgWCgoJUbm6uevbZZ9W2bdvUoEGDHNvvvfee4/n6U6dOqZ49eyqllJo0aZJavHixUkqpgoICZbFY1IYNG9T999/vqPe5555TK1asUEopdfToUdW7d+9a/mTiWlIvzkBTU1MdSyVbrdbLlluyZAm33nor0dHRnD59uhYjFM4WHh5OfHw8v/76K6GhoYSFhREfH8/OnTuJiIjg+eefJzg4mDvuuIODBw9is9mIjIzk9ddf54033uD8+fOOqfUutm7dOl588UVCQkK45557OHfunBs+nbhWlP8Nq4MCAwNZt24dw4YNu2yZ06dPs2nTJtatW1eLkQlXiYiIYOvWrVgsFkwmE3369GHr1q0kJSURHx9PcXExv/76K3q9noYNG2KxWBg1ahTh4eGsWrWKAQMG8NNPP5WrVynFTz/9RPPmzd3wqcS1pl6cgZpMJgICAhzbSikee+wxBg4cyB133EFWVhZr1qzBZrNx6623MnnyZMdcf6J+ioyMZPHixXTv3h2APn36sHjxYrp27Up+fj5NmzZFr9ezatUqx4Jgx48fp0OHDkyZMoX+/fuTlJSEt7c3eXl5jnpjYmJ45513HNsJCQm1+8HENaVeJNBLff/997Ru3Zr169czadIk3nvvPc6dO4fZbGbdunV4eXmxYsUKd4cpaqB3796cPHmS3r17A9CkSROg9Mx09OjRbNiwgfDwcFavXk3r1q0B+OKLL+jWrRshISEUFhYyePBggoODKSgocHQizZw5k9OnTxMcHEzXrl1ZunSp2z6jqP/q1Xyg0dHRxMbG8uabb7J06VICAgKwWq1ERkZyww03oNPpeOSRR1izZg3x8fHMmDHD3SELIa5h9eIe6KU6derE2LFjmTp1KgAWi4X9+/fzwQcfALBnzx7HIl5CCOEqtX4Jn5ycTNOmTYmOjnas4HglFouFmJgY9u7dy6BBgwgKCiI5OZmBAwcycOBAfvrpJ0JCQvD09CQ6Oprdu3czfPhwF38SIcT1rtYv4ZOTk3nhhRf47LPParNZIYRwOrd0Im3YsIF+/frx1ltvuaN5IYRwilpPoM2aNeP3339nw4YNxMbGyjASIUS9VeudSEajEaPRCMBf//pXEhMTCQ4OLlNm0aJFLFq0CID9+/c7hqlcm+rNIAghrikZGZmkp6fXqI5avweal5eHj48PAGPGjGHy5Mn06dPnsuVDQ0PZtGlTbYVX63x88q5cSAjhdGFhQ4mPj69RHbV+Cb9lyxZ69epFVFQUzZs3rzR5CiFEXVbrl/BDhgxhyJAhtd2sEEI4Xb0cSC+EqNy5cwVkZpa4O4w6w2DQ0qiRCX9/k1PrlQQqxDUoM7OEjh3bo9Pp3B2K2ymlKC4uITk5xekJtF5OJiKEuDJJnqU0Gg2ens5NnH+SBCqEcJrk5GTGjBnv7jBqjSRQIYSoJrkHKsQ1SrMxzul1quioMtslJSWMGDESs9mMv78/gwdfmCDoqaemsGdPAna7naVLPyEoKIh77rmPgoJCGjduxJdfLmPcuIdJTj6BVqtl3brVaLX165xOEqgQotqWL19JVFQk06Y9x2OPTSrz2quvvoyXlxexset4//0PefDBsTRq1Ijvv/8QpRQWi4VTp06xaVMsSik0Go2bPkX1SQIVQlTb8ePJBAd3AyAkpEeZ1157bT7r1m3AYrHQpUtnOnToQPfuNzF69Dh69erJlClPM27c3xgzZjxt2rTmpZdm17sz0PoVrRCiTrnhhrbs27cfgISEfY79GRkZbNy4mS1b1vPSS7NQSlFSUsIzzzzF0qX/ZfXqtZw7d46RI+/ns88+4fz58+zeXbPHKt1BzkCFENV2991DGTFiJIMG3YG3tzenT6cCEBAQgLd3AwYOHOQ4Qz1x4gQPP/woVquVdu1uwMPDg1tvHYzNZsPX14fu3bu586NUiyRQIa5Rl3b4uILRaOTbb79Er9fz2GOTGDgwmhkzpgHw/ffLy5XfsmV9me3Nm+v3MuSSQIUQNXLHHXeRn19Ahw7tiYyMcHc4tUoSqBCiRtas+cHdIbiNdCIJIUQ1SQIVQohqkgQqhBDV5LYE+uabb9K3b193NS+EcLKzZ8/yyivzAOjb9xYA5s17ndOnT7szLJdySydSSUkJe/fudUfTQggXCQoKcgxh+tO0ac+5KZrKrT+2/sqFqsAtZ6Affvgh48aNc0fTQlw3NBqj078uFRe3nT59+jJw4CAWL/6k3FR248f/nSNHjlBQUMDw4Q8wYEAMDz44AYDPPltGREQ/br45mr17S5c3j4zsz6RJTxESEs7q1WtcclzWH1vPq1tfdUpdtX4GarFY2LRpE0888QQzZ86ssMzFyxqnpqaye/fu2gyxVmm1NneHIK5Bfn4NXVJvXl5xme3vvvueWbNm0a9fP5KTk1mzZh15ecXYbHby8oqxWGzk55fwf//3Hv37D+DBBx/EbreTnV3AW28t4Oef15Kamsqzzz7HV199yfnz6UyZ8iyTJ1t49tnnuPnmAU6LvbjYzJp1v3Ik8zwjGjzOIj6ucZ21nkA//fRTRo0aVWmZiRMnMnHiRKB0WePw8PDaCM0tZFlj4QoHD2a6pF4fn7Izuz/zzBO8/PKrfP75MiZPfhyDQYePjwmdTouPjwmDQYe3t5GUlOM88cRjjvefPXuWdu3aEhjoQ2BgJwoK8vHxMdG0aRPatWsFQH5+Xrn2asLDqOfllKc4mnmMzo06Ad41rrPWL+GTkpJYuHAhgwcPZv/+/SxYsKC2QxBCOElAQADvvruAf//7FWbNmnvZcp06dWTHjp0A2O12GjduTHLyCSwWC8nJyfj5+QKUmdJOKeW0OJWCjKIMjmYec1qd4IYz0H//+9+On/v27cvkyZNrOwQhhJO8//4HfPvtCvLz87n//hHs25dYYbkJEx5m7NiH+PTTZbRv347FixcxadJj9Os3EK1WyzvvvO3SOFPzUimxOn+VUo1yZpp3gdDQUDZt2uTuMFxGLuGFKxw8mEmXLh3dHUadkFmYycncU6SdSGfc1tsBiGzemxPzbcTH12wKPRlIL4S4ZhWYCzmVV3YcaosGzTGt9XBK/TKZiBDimmS2WUjOTi5zL7WhZwANVhk5U3SeetmJJIQQrqYUpOSkYLVbHfs0aGi1pRm//3bcae1IAhVCXHPO5Z+lwFzg2NZoNBgsBvZsOuDUdiSBCiGuKXkl+ZwrSCuzr4HFC3OR2bHdMCjIKW1JAhVCOEVFk4nUNovNSkpOSpl9PnZv8jMunI2269SREicNn5QEKoRwioomE7ladru92u+t6L6nD97knc93bLe6oS3Z//wn+Z6eNQnTQXrhhbhGbdTEOb3OaFV2obq4uO0888xzNGjQgDFjRrJ+/UY+++wToHTWtb/+9W5+/vknAG69dRCrV3/PmjVree21N7FarcycOZ3BgwcREdGPnj174OnpyZtvvl6t2NIK0sg3X0iW3toG5J25sK3T6bC88AKZ3t40qVYL5UkCFUJU248/rubf/36F6OgBHD9+nPXrNzpeMxqNNG7cmJSUFOx2O61atUKn0/HGG/+P9evXYLfbuf32Oxk8eBDp6enMmDGNli1bViuOAnMB5/LPOrY9dSYKzxY6tvV6Pfj5cbbEr9qftSI1SqD9+/evcL9SihYtWvDFF1/UpHohRB33+OOP8PLLr/LRR58wefLj5V6/995hfPPNd9jtdu69927S09M5ePAQMTGlTwSlpaWhlKJJkybVTp5Wu40TOSn8OdrToNVjzbBit5fu0el06Jo3x3biZLXqr0yNEqifnx+rVq2q8LU777yzJlULIeqBPycTSU1N5eGHH6Fhw7LT6A0ZMpi77x6ORqNh0qTHMRgMdO9+E2vW/FB6SW2xoNFo0Gqr3x1zMuckFpsFAK1Ggz5fT1FJ6bR7Wq0WQ/PmFBtcc7Fdo1ovlzyv9JoQ4tpwpclEPD09CQgIwGAwYDSWTsg8ZcpT3HrrYDQaDV27duadd/6v2u2fzT9HbkmuY9vL6kV+XmmPu1arwdi8OUUehmrXfyVOmUxk3rx5TJs2jdWrVzNz5kzGjBnDk08+6Yz4ZDIRIarhephMJKMwg1O5F55z99Z4k3/2QqeRV7NmFHpd6G1PO3KccYWlTyF19vIi9+8z68ZkIj///DMAS5cuZdu2bXz66afOqFYIISqUU5zL6YuSp5fOk4JzF5Jng8aNyyRPV3FKAi0uLmbJkiU0adIEg8GAp5PGWAkhxKUKzIWk5JxwdBp56AyYz5v581q6gZ8fBb4+tRKLUxLokiVLsFqtzJkzh5KSEiZNmuSMaoUQooxiawnHs49j/yNb6jRatDlarJbStcVMXl4UNnTNelAVcUrXlKenJ/n5+bz11ltXnIY/MTGRiRMnotPp6NChA4sXLy4zjb8QQlTEbLNwLOsYNntpstRqtBgLjRQWFAFg8PDA0rQJqhbTiVPOQO+66y4aN25M//79GTBgAAMGXH4lvU6dOhEXF8eWLVsAanwTVwhx7bPZbRzPOn7RcCUtpkIjhbmlyVOn06FpFoStBsOhqsMpZ6Bt27Zl5MiRVSprMFwYUmA0GmnVqpUzQhBCXKPsSnE8O5li6x9jOy9JnhoNeDQLokhf+w9WOqXFzMxMevXqRdeuXYHSufeWLFly2fIrV65k+vTpdOzYsdzAW5B14YWoKT+/huXWcK+vLDYL3sofb4M/Wo0WjV2D3cuOr1fp61qDAbtWi5+18nqK7TAtSweAKcfKIifE5pRxoCdOnCi3r02bNld83+TJkxk4cCDDhg27bBkZByrE1Tt4MJMuXbs7v2JVdmVLu93OxImPc/jwEby8PJk69WleeGE2AC+/PJuYmFuJjr6NsLBQtm6N46GHxjFx4t8ZP/7veHl5kph4gOjo/sydO4vz58/z8MOPkJeXT5cunXj33QW8MHMW+37fT9qZNOa9+wrNvYIcZ55QOlypqj3udXYcaGBgIP/73//417/+xf/+9z8CAgIuW7ak5MI/gK+vrwx5EqIeW7FiFU2aNGbTplh++GEFs2e/zNq1P7B27Q/MnHlhnfhhw+5i27aNfPLJp5jNpRMb33xzFJs3r+PXX3/j9OnTzJv3Ov/85z/YsGEtPj4+bN++g9ySPNq0a8WiL98tnzybNKm14UqX45QEOnr0aFq2bMmUKVNo0aIFo0aNumzZ1atXOzqazp07x1/+8hdnhCCEcIPffz9MVFQEUPrcuUajwdfXF19fX3Q6naNcz54h6HQ62rRpTVpammMfQPfu3Th+PJmDBw8xbdoLREffxrp1G/g9+TBmm5mbenTFVGQsmzybNqXAp+aLwtWUUxJoTk4Oo0aNolOnTowePZrc3NzLlr3rrrvYtGkTmzZt4sMPP6zRJAJCCPfq1OlGduzYBZReztvtdnJzc8nNzcVmu3B/f+/eBGw2GydOpNCkSRPHPoDExP20bduGTp068uabr7Fx48/s2LWVkAE9APBUJgpzLkqeQUEUeDeorY9YqRp1Ii1YsIDJkycTEhLChAkTCA0NJT4+nuDgYGfFJ4Sow4YOvZNVq36kf/9b8fZuwMyZ0/nLX+5AKcXcuTMd5b766huefvpZHnxwLB4epWuyb9q0mXfffZ8BA/rRsmVLpk9/nokTHycnJxeL3cw/X/sHBo2eopzSzjCNBryCmlFQC49oVlWNOpEGDhzI+vXrgdLxnMeOHaNdu3aEhYU5LUDpRBLi6tWlyUSio28jNvan0kmN/zB+/N954YVpdOjQoVz5AnMhRzKPoNVoMOQYKCkyo9Fo8GzWjEJPU7XjcEUnUo3OQI8dO8bMmTPL7EtMTGTlypXMnTv3Mu8SQoiKKaU4lXsKAC+7F/lFBaXJs3kzCk3VT56uUqMEGhAQQExMzBUf3xRCXL82bvy53L5PPvmwwrLnC9MpthbjqTORf6Z0Xk+vxo0oqIPJE2qYQP39/S+7rIcQQlwNs9XCufyzaDUa7Fl2UODZoAEFPu4dqlSZGnWBz5kzx1lxCCGuc6fzTmNXCi97A0qKzKVLfjRu7O6wKlWjBDpv3jyUUo7hC39+2Ww27rjjDmfFKIS4xp3LTyO3JLf00v2PddyNTZpg1dXtYY41uoQvKiri1ltvLbdfKUVQUFBNqhZC1EMbN24iNnY9L79c9avTM3lnSStIQ6vRMHfKKzz37NTSSZHr0HCly6lRAt2wYYOz4hBCXIdS81I5X5AOlPa6P/fsVAwGA8WBgW6OrGpqf/4nIUSt0MwxOr1ONavsZCKpqamMHj0ei8VCcHA37rtvOADJycm88MJsPvvsEzZu3MTGjZuZPftFQkLCCQnpwS+//MqT0yax5OPPyEjL4P1P3wVN6fjQL2JXM/aOu+jRO4wdGzcx4qEHGTZuDONvG8InP/8IwMO3D+Wjn1by8O1DCYnozfZ167n3wfH8sm0bR/Yf5OmXZhMVM9Dpn/9SkkCFENXWqFEjfv75R/R6PWPGjOfw4SOVlj9z5izbt29h9ZY1TJ/6IsvWLGHjis2s+HQlf/vbaHQGPcXG0ieV/nLP3Tw6/XkeHXovw8aNuWydt909lEemPcdtN97Et/Fx2Gw2Xnnq2VpJoE65Q/vJJ58AsHv3boYMGcLXX3/tjGqFEHVcRkYGw4c/QHT0bWzdGkdq6hmAMsv0XDxOvEOHdqSVpGEMMNKuY1t88cHH4ENubi5GTxM23YVzuvZdu+BhNKK5ZL6MS8ed/1mubccbadikCU2aNSM3O9sVH7ccpyTQzz77DID//Oc/fPTRR8ybN88Z1Qoh6rhly77g7ruHsnHjz9x8cyQDBvQDwM/Pj7NnzwKwb18iUDqzvNlmJrs4Gx1a9HYD+ecL/qhJg/pjkpE/XbpWmlIKc0kJh/cfqLDc5ZK2KznlEj4/P5/Nmzfj7e1Ns2bNaNCgbsyUIoRwrYEDb2Hs2IdYvnxlmf3+/v60bt2amJjBdOjQnqCgIE7lnsJqt2HQGtAX6LGaS6eQ1+t1aH19MF9hSY6ho0cy/i9DuDmm/Mgfd3HKjPTbt2/nxx9/5IknniAgIID33nuPp556yhnxyWQiQlRDXZpMBC4MVTLqjNgz7FhKSheHM3l5Ym3StFbGe9bZGelzcnJ46aWXyM/PZ8qUKfTs2fOyZXfu3ElUVBT9+vXjmWeecUbzQog6LKMwg7SCNDx1nljPWRzJs4GfHyVBzer8YPnKOCXyN954A4B//etfjBkzhqlTp162bJs2bVi/fj1btmwhLS2Nffv2OSMEIUQdlFOcy+nc05h0JsxpJdhsduCPtYwaNazVNdxdwSn3QPPy8khJSUGn0xEZGVnpPdCLn1DS6/Vlpv0XQjiPUqpcR0xtKjQXkpJzAg+dEUuaBZvVjlarxRQUREEN5vWsDqUUdpzfseSUM9Bp06YxY8YMnnvuOYqLi4mIiLjiexISEkhPT3cshSyEcB6DQUtxccmVC7pIidXM8exk9FoDtnQbNqsNvUGPoWXLGk2KXB1KKazZeZywm51et1M6kQCKi4s5f/68Y/hA69atL1s2MzOTu+++my+//LLCZ+YvXhf+1KlTLFu2zBkh1kmyLrxwBZ3OgLe3r9vat9hL73NiA6VKl+NAp3fPJbuCAouZpKw0LPbSWwgmrY5FL86tcSeSUxLovHnzWLt2LYcOHaJ9+/YYjUZiY2MrLGu1Whk6dCizZs2iT58+V6xbeuGFqD9SslN4/MfHSc1JRfelltNJZ2jSLAiP2bM51bChu8NzqFO98CtXrmT9+vV07NiRLVu20LCSA/XVV1+xe/dunn/+eaKjo9m+fbszQhBCuIHZZmbLiS3M3jib/h/3J+KjSFKyUzCtMHI66QxNmwVhmDOnTiVPZ3JKJ9Kfq+x5eXmxefNmDhw4cNmyI0eOZOTIkc5oVgjhBkWWIr479B2xx2LZmrKVfHPp00RGnQfdfG+i+PsijiQcI6h5M3SzZ3O6nsysVB1OSaBvv/02JSUlzJ8/n4ULFzqGNQkhri0Hzh/g0e8f5UjmUXQaLR182tPY1pjCQ/n8vvEwibn7AWjWojna2bM5HRDg5ohdq0YJNCUlBShdXO7cuXM0aNCAZ5991imBCSHqlo9/+5iXNr9EW6+29D7ViyMbjpCUdpgkDpcp16xFczRz5nDa399NkdaeGiXQtm3bEh4ezk033eTY9+fYs8WLF9c4OCGE+2UXZTNl7RRij8USYe/Njld2YTWXHz3SpFkQbUJDOXXPPaT6+bkh0tpXowS6c+dOli9fzv79++nSpQv33nsvYWFhzopNCOFmu07v4okfn0Bjg86/dGTr+gudvo2bNqFNaCi6bt041aEDpxs2JM2NsbpDjRJoeHg44eHhQGkynTBhAmFhYXzwwQdOCU4I4R5FliLe/+V95sfNJ9wzjMT/JHIqKxWA3nffxdnBg0lp1Ijzbo7T3WqUQNPT01mxYgXr1q3DZDIxY8YMhgwZ4qzYhBC1wGwzcyDtAHvO7WHvub3sPbuXwxmH8TH6EHq6Jzs+3wWU3ttsNGkSu2680c0R1x01ElFXkAAADp9JREFUSqDNmjWjZ8+e3HLLLXh4eJCQkEBCQgIAc+fOdUqAQgjnKrYWs/PUTjYmb2T7qe0cSj+E2WZBr9XRukErGqlGRFoiOPLtEXYf+wWNBnrfey+J99zDmT+GLIpSNUqgl3vaSAhRtxzOOMzG5I1sPLGRHSd3UGwtpr1Pe5qYG9MrM5TspGySfz3BsYJkjpHseF+LVi0JePJJdt5wg/uCr8NqlEAHDBjgrDiEEDVktVs5nXua5OxkkrOTOZFzguTsZPal7eN0biqNPRvRTnsDwWe6cXTDMY6cPMoRjlZYl9HoQc9hw9g7dCin5azzsmRVTiHqKaUUO07tYPmh5Ww9uZWTOSex2m1oNRoaezamoT4Qb7s3N2S2wWuXJ0d2HeW8Si9XT7MWzWnaoQOe7dphbtmS882bcyowkB0y1eQVSQIVop7Zd24f3x36jpVJK8ktyaWLoQtNTzfB54w3WSeySDuWxrmiNM5VMKhIp9PSvmtXGvbsSUGXLhxv1YozJhNn3PA5rgWSQIWo4zKLMjl4/iA7Tu1gRdIKTuee5iZTV1okNSPr+yx2519+RiG9XkebG2+kcUgIxTfdxNG2bfndVLvzcV7LJIEKUUcopUhM28/B8wc4lHGIg+kHOXj+IGkF5/E3+dHO4/+3d/8xUd95Hsef8wNmEOTXimhtRe1aTVM3VfwFgtCr1dpYW7px01pTs9k9GpsmFy+XXKNNzmQvjam3f5jtZu+u215sYs+a1u3qtbvZRbOliqJWhfKrgij+4JcwA8MMzMx3Zt77BzqniK2OMF/Q98MQQ/KZ+b54+/m+Hb6/PjPJPJdB+4E2TrpPRV+XnpHOjHnzcGRnY83KwsjMxJuRgSszk6vJyZyzWm9zpFPdK22gSplMRPiy8Ut2Vu7E7XORZctiwsAE6BYmt2QhtRGuXunmFGeir5k8dQqzCgvx5uZSP2MGZ/R4pSm0gSplkohEOPDdAX5z/Dc4/IlwWLhS1saV2xyRfGTmDB4uLMS1YAHfTZtGp4nrHalB2kCVirNwJMwfG/7Ib0/8lpSBFIJ/DlB3vB4Ap9PBQzk5ZOTkkDhtGqGpU+nNyqI9K4tLTieXTM6ubmZKA21tbWXNmjXU1dXh9Xqx27WPq/tXRCI0uZqobq+mqqOKiosVpPvS8OzzUF/fAMATixdhf/55qubModk6ftdJf9CY0rkyMzM5ePAgJSUlZmxeqVHV3d/N3y78jTPtZ/i281tcXhepkVScfQ6MVoPOw53UX2rA4Uhk0drnca1cSc0wiyuqsc+UBup0OnHqpRTqPuIacPGnxj/xl6a/0NPdg+2Kjd7GHi5VXcbr8t00NnPSj8j75S/4rqCAEykpJiVWI2FM/u5847LGra2tnDhxwuREo0eXNR6/whLG4+/FE+gjZIQQfwqLPKuJhAUygSWQUJCA3eHA6kiExERC9gSCCQmIBYoMwG32T/FgcvaG+O8ReJ8x2UBLS0spLS0FBpc1vv7M0fuRLms8PrT2tdJwdfDazAuuC3S7u/C6vPSe9lD31zrCRgSr1cLsn/yEzCVLcM+dy+XJk/E6HMO8m/6naba5E4b7d7l7Y7KBKmWmUCTEqbZTVLZUcr7lPJ3dHfR39eNt8dJe107XJVd07MSJKTz51NNYc3M5N3cu3yUnm5hcxZspDdQwDFavXk1VVRWrVq3inXfeYcmSJWZEUQqAlp4Wjlw4Ql1jLc3nmmk700bT0XOEjUh0TFb2ZLKnz2LWoiJsDz1E3/TpNOTk8I1exP7AMqWBJiQk6LNEVVwZYYO2vjZava109HTgcrvodnXh6nHTcbWdlpqLNJY3YvhDWK0WZs2dy6KfriP04x/TM3kybZMmcTUx8YFfwkLdTH+FV/edAWOA2o5aGpsbaW5upvlCM5cvXab3Si/dl7rxuftvGp+QYOfRJ55g4c9eoX/OHJpzcmhyOmkyKb8aP7SBqnGlsa2Rrp4uvD4vvn4fAwMDeAe89Pf309vby7lzTTQ3NtNSfZFQYPBkzeSpU5g8YwbZ03OY9mQmkp5OKC2N/rQ0elNT6UpNpUFv5lAx0FmjxgWf38ev//PX/H7n76ONcaiUlGQenj2bSTOfZFLBC/RNm8aV7Gw6k5IeuOV2VXxoA1Vj3qmGU2zdupXutl4WbP4X/BMnYjidBBwO/E4n/QkJ+BwOvHY7DWaHVQ8UbaBqzApFQvxu1+947z/e44kXfkrDPz3LZf1VW40hOhvVmNR0uYktW7dwtcND1r9v52h2ttmRlLqFNlBlir6BPnq9vfgDfvx+P37/AP5AAH/AT3PzOf7rvffJKVnH2X9chujTidQYpQ1UxU0oHOLQ1wf55JNPKPvzQYxACKfTgcPpxJGURKLDQaLDQcaMGci//YpjaWlmR1bqe2kDVaPuSscV9n62l08/+ZTkjCySnn6apHUbMBwO/BYLfrMDKhUjbaDqntTU1+Dz+pDrf0QAEARPr4c//GEfdfVnyS56Cv+/buV8errJiZUaOdpA1V0LhULs/7/9/M+uXYQcSdiTkrBYLFisVrBYBr+sVmzJybiLn6F5w89pNju0UqNAG6i6Y54+D7s/3s2Xfy0j8cn5nH/jTVz69CH1ANMGqohEIrRcbCEcCmO1WbFardhsNqxWK1abFZ/Xx//u2UPV5cv0LVnKt5v/efBTplIPOG2gY5SIYBmlJnXx4kXOnDlDVXUVF9racBsGlswfwfWL1K8dx4xKSqJ54UI6nvqHUcmj1HilDXQM6ezspKysjGNnTtPm6yeC4LDbcdjtOO12kuwJOBMG/44gGJEIIRGMcJhwJIIhggACXG+9VixYrn0vFnAFDfzp6fQ9/DCXl+bRM2GCaT+vUuOdKQ108+bNnDx5kgULFrBz504zIowZtbW1fFVRQV13N612O62zZ3PxRV2tVKnxIO4N9NSpU/h8Pr7++ms2bdrEiRMnxuyaR5FIBJfLRUdHB13d3UREBo8PWizYbDbsdjt2ux2bzTZ4N41hEIhECF77CoTD0U+FFqsVrp2lvn62WiwWPu2tp3HuXML6VHOlxp24N9CjR4+yYsUKAFasWMGxY8fi0kB9Ph99fX14g0HCkcFlGiww2NBuGNcXDOL2++kJBOgLBvE7nXiTk/FNn37T+1lEsIpgi0SwihCYOBGJ4ZhlQ5rvhwcppcakuDfQnp4eHn30UQDS0tKora393vERETzBIMFwOPrJzohEMMJhIkPG3ti+BkIhvIaBzzDwhUL4bTYCiYnIDz3NJyEBbrg0xwqkXvsaDfP0MiCl4m6W08lILJYe9waanp6Ox+MBwOPxkD7MnSk3rgvfUF/PC6tWxTXj7Vy9epWsrCyzYwBjJ8tYyQGa5XY0y626gYaGEXh6rMTZN998I6WlpSIismnTJqmsrPze8bm5ufGIdUc0y63GSg4RzXI7mmV4I5El7s8JW7BgAU6nk8LCQqxWK4sXL453BKWUGhGmXMb0oF+6pJS6P9i2bdu2zewQPyQ3N9fsCFGa5VZjJQdoltvRLMO71ywWkaH37SmllLoTulaCUkrFyLQGWllZSX5+PoWFhWzevBmAHTt2UFBQwKuvvophGBiGQV5eHikpKTQ1NUVfO3ScmVnS0tIoLi6muLgYl8sVlyznz5+nsLCQ5cuXs379esLhwXXSd+/eTX5+PmvWrIleKmZWljlz5kTrUldXF5csXV1d5OfnU1RUxNq1axkYGBh2nJlZzJgv13322Wc88sgj0e9Hsi73ksOsmgy33bveh+75PH6M2traZGBgQERE1q9fL+Xl5bJ69WoREdm+fbvs3btXIpGItLe3y8aNG6WxsVFERDo7O28ZZ1YWEZFly5bd8/bvNovL5RK32y0iIlu2bJH9+/dLMBiUgoICMQxD9uzZI++++65pWUTMqUsoFJJwOCwiItu2bZO9e/eaNl+GyyJiTl2ue/nllyUvL09ERn4/ijWHiHk1GbrdWPYh0z6BTpkyBafTCYDdbqe6upri4mLg/2/xtFgsZA9Zzvb48eO3jDMrC0B9fT2FhYW89dZb0eUsRjtLRkZG9AaE6/finz17lnnz5mG32+Nal+GyALhcLpYvX87rr7+O33/vqx7dSZbrzzAFCIfDzJ4927T5MlwWMGe+AHzxxRc888wz0UwjXZdYc4B5NRm63Vj2IdOPgVZXV9PV1UV6ejqpqYM3TKalpeF2u4cd39PTc0fj4pEFoLGxkfLyctxuNwcOHIhrltbWVsrKyli5cqXpdbkxC8Dhw4cpLy8nJycneldZPLIcP36chQsXcujQIWbOnGlqXYZmAfPmy65du9iwYUN0/GjV5W5zgHk1GbrdWGpiagN1uVy8+eabfPDBB3d0iyfc2a2g8coCkJmZicVi4cUXX6SmpiZuWQKBABs3buT999/HbrebWpehWWCwLgAlJSVxrcvixYs5efIkJSUlfPjhh6bWZWgWMGe+HDp0iPz8fBITE6OvGY26xJIDzNuHhm43lpqY1kBDoRAbNmxgx44dTJkyhUWLFvHVV18BUFZWxtKlS4d93Z2Oi0cWn88XPWly5MiR6ENS4pGltLSUN954g8cffxyAxx57jJqaGsLhcNzrMjRLMBgkEAgA8a1LMBiMjk9NTSUpKcm0+TJcFrPmS01NDfv37+fZZ5+ltraWt99+e8TrEmsOs2oy3HZj2odG8NjtXfn4449l0qRJUlRUJEVFRVJRUSHbt2+XZcuWySuvvCKBQEBERNatWydTp06V/Px8+fzzz0VEhh1nRpbTp0/L/PnzpaCgQF577TUJhUJxyVJRUSEpKSnRMfv27RMRkY8++kjy8vLkueeek56eHtOytLe3y/z586WwsFDWrl0rHo8nLlkqKytl+fLlUlxcLC+99JL4fD4RMWe+DJfFrPlyoxtPnIxkXWLNYVZNbrfdu92H9EJ6pZSKkeknkZRSarzSBqqUUjHSBqqUUjHSBqqUUjHSBqqUUjHSBqqUUjHSBqqUUjHSBqqUUjH6O28zy/8S9glnAAAAAElFTkSuQmCC\n", - "text/plain": [ - "
" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], - "source": [ - "# USyearly.to_csv('PVSC_USYearly_DataforPlot.csv')\n", - "\n", - "\n", - "\n", - "plt.rcParams.update({'font.size': 8})\n", - "plt.rcParams['figure.figsize'] = (12, 8)\n", - " \n", - "keywords=['VirginStock_', 'Waste_']\n", - "SFScenarios = ['Reference.Mod', '95-by-35.Adv', '95-by-35_Elec.Adv_DR']\n", - "materials = ['glass', 'silicon', 'silver', 'copper', 'aluminum']\n", - "\n", - " \n", - "fig, axs = plt.subplots(2,1, figsize=(5, 5), facecolor='w', edgecolor='k')\n", - "fig.subplots_adjust(hspace = .6, wspace=.001)\n", - "fig.subplots_adjust(hspace = .3, wspace=.001)\n", - "\n", - "axs = axs.ravel()\n", - "i = 0\n", - "\n", - "### PLOT 1\n", - "\n", - "# Loop over Keywords\n", - "ii=0\n", - "keyw = keywords[ii]\n", - "# Loop over SF Scenarios\n", - "kk=0\n", - "\n", - "obj = SFScenarios[kk]\n", - "axs[i].yaxis.grid()\n", - "axs[i].axvspan(2000, 2018, facecolor='0.9', alpha=0.5)\n", - "axs[i].axvspan(2018, 2050.5, facecolor='yellow', alpha=0.1)\n", - "axs[i].plot([],[],color='c', label='glass', linewidth=5)\n", - "axs[i].plot([],[],color='k', label='silicon', linewidth=5)\n", - "axs[i].plot([],[],color='m', label='silver', linewidth=5)\n", - "axs[i].plot([],[],color='r', label='copper', linewidth=5)\n", - "axs[i].plot([],[],color='g', label='aluminum', linewidth=5)\n", - "\n", - "axs[i].stackplot(years, USyearly[keyw+materials[0]+'_'+SFScenarios[0]], \n", - " USyearly[keyw+materials[1]+'_'+SFScenarios[0]], \n", - " USyearly[keyw+materials[2]+'_'+SFScenarios[0]], \n", - " USyearly[keyw+materials[3]+'_'+SFScenarios[0]], \n", - " USyearly[keyw+materials[4]+'_'+SFScenarios[0]], \n", - " colors=['c','k','m','r', 'g'])\n", - "#axs[i].ylabel('Mass [Tons]')\n", - "axs[i].set_xlim([2010, 2050])\n", - "#axs[i].set_title(keyw+ ' ' + obj.name)\n", - "#axs[i].legend(materials)\n", - "\n", - "i += 1 \n", - "\n", - "### PLOT 2\n", - "\n", - "ii=1\n", - "keyw = keywords[ii]\n", - "# Loop over SF Scenarios\n", - "kk=0\n", - "\n", - "obj = SFScenarios[kk]\n", - "axs[i].yaxis.grid()\n", - "axs[i].axvspan(2000, 2018, facecolor='0.9', alpha=0.5)\n", - "axs[i].axvspan(2018, 2050.5, facecolor='yellow', alpha=0.1)\n", - "axs[i].plot([],[],color='c', label='glass', linewidth=5)\n", - "axs[i].plot([],[],color='k', label='silicon', linewidth=5)\n", - "axs[i].plot([],[],color='m', label='silver', linewidth=5)\n", - "axs[i].plot([],[],color='r', label='copper', linewidth=5)\n", - "axs[i].plot([],[],color='g', label='aluminum', linewidth=5)\n", - "\n", - "axs[i].stackplot(years, USyearly[keyw+materials[0]+'_'+SFScenarios[0]]*907185, \n", - " USyearly[keyw+materials[1]+'_'+SFScenarios[0]]*907185, \n", - " USyearly[keyw+materials[2]+'_'+SFScenarios[0]]*907185, \n", - " USyearly[keyw+materials[3]+'_'+SFScenarios[0]]*907185, \n", - " USyearly[keyw+materials[4]+'_'+SFScenarios[0]]*907185, \n", - " colors=['c','k','m','r', 'g'])\n", - "#axs[i].ylabel('Mass [Tons]')\n", - "axs[i].set_xlim([2010, 2050])\n", - "#axs[i].set_title(keyw+ ' ' + obj.name)\n", - "#axs[i].legend(materials)\n", - "\n", - "i += 1 \n", - "\n", - "\n", - "axs[0].set_ylim([0, 2e7])\n", - "axs[1].set_ylim([0, 5e6])\n", - "axs[0].set_title(\"Virgin Stock Needs\")\n", - "axs[1].set_title(\"Waste\")\n", - " \n", - "axs[0].set_ylabel('Mass [Tons]')\n", - "axs[1].set_ylabel('Mass [Tons]')\n", - "axs[1].legend(materials)\n", - "\n", - "fig.savefig(os.path.join(testfolder,'testplot.png'), dpi=600)\n" - ] - } - ], - "metadata": { - "kernelspec": { - "display_name": "Python 3", - "language": "python", - "name": "python3" - }, - "language_info": { - "codemirror_mode": { - "name": "ipython", - "version": 3 - }, - "file_extension": ".py", - "mimetype": "text/x-python", - "name": "python", - "nbconvert_exporter": "python", - "pygments_lexer": "ipython3", - "version": "3.7.3" - } - }, - "nbformat": 4, - "nbformat_minor": 2 -} diff --git a/docs/tutorials/PVSC 2021 Fig. 3.py b/docs/tutorials/PVSC 2021 Fig. 3.py deleted file mode 100644 index 3fbfdd82..00000000 --- a/docs/tutorials/PVSC 2021 Fig. 3.py +++ /dev/null @@ -1,147 +0,0 @@ -#!/usr/bin/env python -# coding: utf-8 - -# # PVSC Fig. 3 - -# In[2]: - - -import PV_ICE -import numpy as np -import pandas as pd -import os,sys -import matplotlib.pyplot as plt -plt.rcParams.update({'font.size': 22}) -plt.rcParams['figure.figsize'] = (12, 8) - - -# In[3]: - - -import os -from pathlib import Path - -testfolder = str(Path().resolve().parent.parent / 'PV_ICE' / 'TEMP') - -print ("Your simulation will be stored in %s" % testfolder) - - -# In[18]: - - -years = list(range(2009,2051)) - - -# In[9]: - - -USyearly = pd.read_csv(os.path.join(testfolder,'PVSC_USYearly_DataforPlot.csv')) - - -# In[12]: - - -USyearly.keys() - - -# In[19]: - - -USyearly - - -# In[21]: - - -# USyearly.to_csv('PVSC_USYearly_DataforPlot.csv') - - - -plt.rcParams.update({'font.size': 8}) -plt.rcParams['figure.figsize'] = (12, 8) - -keywords=['VirginStock_', 'Waste_'] -SFScenarios = ['Reference.Mod', '95-by-35.Adv', '95-by-35_Elec.Adv_DR'] -materials = ['glass', 'silicon', 'silver', 'copper', 'aluminum'] - - -fig, axs = plt.subplots(2,1, figsize=(5, 5), facecolor='w', edgecolor='k') -fig.subplots_adjust(hspace = .6, wspace=.001) -fig.subplots_adjust(hspace = .3, wspace=.001) - -axs = axs.ravel() -i = 0 - -### PLOT 1 - -# Loop over Keywords -ii=0 -keyw = keywords[ii] -# Loop over SF Scenarios -kk=0 - -obj = SFScenarios[kk] -axs[i].yaxis.grid() -axs[i].axvspan(2000, 2018, facecolor='0.9', alpha=0.5) -axs[i].axvspan(2018, 2050.5, facecolor='yellow', alpha=0.1) -axs[i].plot([],[],color='c', label='glass', linewidth=5) -axs[i].plot([],[],color='k', label='silicon', linewidth=5) -axs[i].plot([],[],color='m', label='silver', linewidth=5) -axs[i].plot([],[],color='r', label='copper', linewidth=5) -axs[i].plot([],[],color='g', label='aluminum', linewidth=5) - -axs[i].stackplot(years, USyearly[keyw+materials[0]+'_'+SFScenarios[0]], - USyearly[keyw+materials[1]+'_'+SFScenarios[0]], - USyearly[keyw+materials[2]+'_'+SFScenarios[0]], - USyearly[keyw+materials[3]+'_'+SFScenarios[0]], - USyearly[keyw+materials[4]+'_'+SFScenarios[0]], - colors=['c','k','m','r', 'g']) -#axs[i].ylabel('Mass [Tons]') -axs[i].set_xlim([2010, 2050]) -#axs[i].set_title(keyw+ ' ' + obj.name) -#axs[i].legend(materials) - -i += 1 - -### PLOT 2 - -ii=1 -keyw = keywords[ii] -# Loop over SF Scenarios -kk=0 - -obj = SFScenarios[kk] -axs[i].yaxis.grid() -axs[i].axvspan(2000, 2018, facecolor='0.9', alpha=0.5) -axs[i].axvspan(2018, 2050.5, facecolor='yellow', alpha=0.1) -axs[i].plot([],[],color='c', label='glass', linewidth=5) -axs[i].plot([],[],color='k', label='silicon', linewidth=5) -axs[i].plot([],[],color='m', label='silver', linewidth=5) -axs[i].plot([],[],color='r', label='copper', linewidth=5) -axs[i].plot([],[],color='g', label='aluminum', linewidth=5) - -axs[i].stackplot(years, USyearly[keyw+materials[0]+'_'+SFScenarios[0]]*907185, - USyearly[keyw+materials[1]+'_'+SFScenarios[0]]*907185, - USyearly[keyw+materials[2]+'_'+SFScenarios[0]]*907185, - USyearly[keyw+materials[3]+'_'+SFScenarios[0]]*907185, - USyearly[keyw+materials[4]+'_'+SFScenarios[0]]*907185, - colors=['c','k','m','r', 'g']) -#axs[i].ylabel('Mass [Tons]') -axs[i].set_xlim([2010, 2050]) -#axs[i].set_title(keyw+ ' ' + obj.name) -#axs[i].legend(materials) - -i += 1 - - -axs[0].set_ylim([0, 2e7]) -axs[1].set_ylim([0, 5e6]) -axs[0].set_title("Virgin Stock Needs") -axs[1].set_title("Waste") - -axs[0].set_ylabel('Mass [Tons]') -axs[1].set_ylabel('Mass [Tons]') -axs[1].legend(materials) - -fig.savefig(os.path.join(testfolder,'testplot.png'), dpi=600) -