docs/tutorials/14 - Advanced topics - Cement Pavers albedo example.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 14 - Advanced topics - Cement Pavers albedo example\n",
    "\n",
    "This journal creates a paver underneath the single-axis trackers, and evaluates the improvement for one day -- June 17th with and without the pavers for a location in Davis, CA.\n",
    "\n",
    "![Paver](../images_wiki/AdvancedJournals/Pavers.PNG)\n",
    "\n",
    "Measurements:\n",
    "![Paver](../images_wiki/AdvancedJournals/Pavers_Geometry.PNG)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Your simulation will be stored in C:\\Users\\sayala\\Documents\\GitHub\\bifacial_radiance\\bifacial_radiance\\TEMP\\Tutorial_14\n"
     ]
    }
   ],
   "source": [
    "import os\n",
    "from pathlib import Path\n",
    "import pandas as pd\n",
    "\n",
    "testfolder = str(Path().resolve().parent.parent / 'bifacial_radiance' / 'TEMP' / 'Tutorial_14')\n",
    "if not os.path.exists(testfolder):\n",
    "    os.makedirs(testfolder)\n",
    "    \n",
    "print (\"Your simulation will be stored in %s\" % testfolder)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from bifacial_radiance import *   \n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "path = C:\\Users\\sayala\\Documents\\GitHub\\bifacial_radiance\\bifacial_radiance\\TEMP\\Tutorial_14\n",
      "Loading albedo, 1 value(s), 0.200 avg\n",
      "1 nonzero albedo values.\n",
      "Getting weather file: USA_CA_Sacramento.Metro.AP.724839_TMY3.epw\n",
      " ... OK!\n",
      "8760 line in WeatherFile. Assuming this is a standard hourly WeatherFile for the year for purposes of saving Gencumulativesky temporary weather files in EPW folder.\n",
      "Coercing year to 2021\n",
      "Saving file EPWs\\metdata_temp.csv, # points: 8760\n",
      "Calculating Sun position for Metdata that is right-labeled  with a delta of -30 mins. i.e. 12 is 11:30 sunpos\n",
      "\n",
      "Module Name: test-module\n",
      "Module test-module updated in module.json\n",
      "Pre-existing .rad file objects\\test-module.rad will be overwritten\n",
      "\n",
      "Module was shifted by 0.078 in X to avoid sensors on air\n",
      "This is a Cell-Level detailed module with Packaging Factor of 0.85 %\n",
      "Module test-module updated in module.json\n",
      "Pre-existing .rad file objects\\test-module.rad will be overwritten\n",
      "\n"
     ]
    }
   ],
   "source": [
    "simulationname = 'tutorial_14'\n",
    "\n",
    "#Location:\n",
    "lat = 38.5449 # Davis, CA\n",
    "lon = -121.7405 # Davis, CA\n",
    "# MakeModule Parameters\n",
    "moduletype='test-module'\n",
    "numpanels = 1  # AgriPV site has 3 modules along the y direction (N-S since we are facing it to the south) .\n",
    "x = 0.95  \n",
    "y = 1.838\n",
    "xgap = 0.02# Leaving 2 centimeters between modules on x direction\n",
    "ygap = 0.0 # 1 - up \n",
    "zgap = 0.06 # gap between modules and torquetube.\n",
    "\n",
    "# Other default values:\n",
    "\n",
    "# TorqueTube Parameters\n",
    "axisofrotationTorqueTube=True\n",
    "torqueTube = False\n",
    "cellLevelModule = True\n",
    "\n",
    "numcellsx = 6\n",
    "numcellsy = 10\n",
    "xcell = 0.156\n",
    "ycell = 0.158\n",
    "xcellgap = 0.015\n",
    "ycellgap = 0.015\n",
    "\n",
    "sensorsy = numcellsy   # one sensor per cell\n",
    "\n",
    "cellLevelModuleParams = {'numcellsx': numcellsx, 'numcellsy':numcellsy, \n",
    "                         'xcell': xcell, 'ycell': ycell, 'xcellgap': xcellgap, 'ycellgap': ycellgap}\n",
    "\n",
    "# SceneDict Parameters\n",
    "gcr = 0.33 # m\n",
    "albedo = 0.2  #'grass'     # ground albedo\n",
    "hub_height = 1.237 # m  \n",
    "nMods = 20 # six modules per row.\n",
    "nRows = 3  # 3 row\n",
    "\n",
    "azimuth_ang = 90 # Facing east \n",
    "\n",
    "\n",
    "demo = RadianceObj(simulationname,path = testfolder)  # Create a RadianceObj 'object'\n",
    "demo.setGround(albedo) #\n",
    "epwfile = demo.getEPW(lat, lon) \n",
    "metdata = demo.readWeatherFile(epwfile, coerce_year=2021) # read in the EPW weather data from above\n",
    "mymodule=demo.makeModule(name=moduletype,x=x,y=y,numpanels = numpanels, xgap=xgap, ygap=ygap)\n",
    "mymodule.addCellModule(numcellsx=numcellsx, numcellsy=numcellsy,\n",
    "                       xcell=xcell, ycell=ycell, xcellgap=xcellgap, ycellgap=ycellgap)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Material exists, overwriting...\n",
      "\n",
      "Added material sw_chantillywhite to file ground.rad\n"
     ]
    }
   ],
   "source": [
    "description = 'Sherman Williams \"Chantilly White\" acrylic paint'\n",
    "materialpav = 'sw_chantillywhite'\n",
    "Rrefl = 0.5\n",
    "Grefl = 0.5 \n",
    "Brefl = 0.5\n",
    "demo.addMaterial(material=materialpav, Rrefl=Rrefl, Grefl=Grefl, Brefl=Brefl, comment=description)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Simulation without Pavers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Created tutorial_14.oct\n"
     ]
    }
   ],
   "source": [
    "timeindex = metdata.datetime.index(pd.to_datetime('2021-06-17 12:0:0 -8'))  # Davis, CA is TZ -8\n",
    "demo.gendaylit(timeindex)  \n",
    "    \n",
    "tilt = demo.getSingleTimestampTrackerAngle(metdata, timeindex=timeindex, gcr=gcr, \n",
    "                                   azimuth=180, axis_tilt=0, \n",
    "                                   limit_angle=60, backtrack=True)\n",
    "# create a scene with all the variables\n",
    "sceneDict = {'tilt':tilt,'gcr': gcr,'hub_height':hub_height,'azimuth':azimuth_ang, 'module_type':moduletype, 'nMods': nMods, 'nRows': nRows}  \n",
    "scene = demo.makeScene(module=mymodule, sceneDict=sceneDict) #makeScene creates a .rad file with 20 modules per row, 7 rows.\n",
    "octfile = demo.makeOct(demo.getfilelist())  # makeOct combines all of the ground, sky and object fil|es into a .oct file."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Linescan in process: tutorial_14_noPavers_Front\n",
      "Linescan in process: tutorial_14_noPavers_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers.csv\n",
      "Simulation without Pavers Finished\n"
     ]
    }
   ],
   "source": [
    "analysis = AnalysisObj(octfile, demo.name)  # return an analysis object including the scan dimensions for back irradiance\n",
    "frontscan, backscan = analysis.moduleAnalysis(scene, sensorsy=sensorsy)\n",
    "analysis.analysis(octfile, simulationname+\"_noPavers\", frontscan, backscan)  # compare the back vs front irradiance  \n",
    "print(\"Simulation without Pavers Finished\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Looping on the day"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_1_Front\n",
      "Linescan in process: tutorial_14_noPavers_1_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_1.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_2_Front\n",
      "Linescan in process: tutorial_14_noPavers_2_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_2.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_3_Front\n",
      "Linescan in process: tutorial_14_noPavers_3_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_3.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_4_Front\n",
      "Linescan in process: tutorial_14_noPavers_4_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_4.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_5_Front\n",
      "Linescan in process: tutorial_14_noPavers_5_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_5.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_6_Front\n",
      "Linescan in process: tutorial_14_noPavers_6_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_6.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_7_Front\n",
      "Linescan in process: tutorial_14_noPavers_7_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_7.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_8_Front\n",
      "Linescan in process: tutorial_14_noPavers_8_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_8.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_9_Front\n",
      "Linescan in process: tutorial_14_noPavers_9_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_9.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_10_Front\n",
      "Linescan in process: tutorial_14_noPavers_10_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_10.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_11_Front\n",
      "Linescan in process: tutorial_14_noPavers_11_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_11.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_noPavers_12_Front\n",
      "Linescan in process: tutorial_14_noPavers_12_Back\n",
      "Saved: results\\irr_tutorial_14_noPavers_12.csv\n"
     ]
    }
   ],
   "source": [
    "j=0\n",
    "starttimeindex = metdata.datetime.index(pd.to_datetime('2021-06-17 7:0:0 -8'))\n",
    "endtimeindex = metdata.datetime.index(pd.to_datetime('2021-06-17 19:0:0 -8'))\n",
    "for timess in range (starttimeindex, endtimeindex):\n",
    "    j+=1\n",
    "    demo.gendaylit(timess)\n",
    "    tilt = demo.getSingleTimestampTrackerAngle(metdata, timeindex=timess, gcr=gcr, \n",
    "                                       azimuth=180, axis_tilt=0, \n",
    "                                       limit_angle=60, backtrack=True)\n",
    "    # create a scene with all the variables\n",
    "    sceneDict = {'tilt':tilt,'gcr': gcr,'hub_height':hub_height,'azimuth':azimuth_ang, 'module_type':moduletype, 'nMods': nMods, 'nRows': nRows}  \n",
    "    scene = demo.makeScene(module=mymodule, sceneDict=sceneDict) #makeScene creates a .rad file with 20 modules per row, 7 rows.\n",
    "    octfile = demo.makeOct(demo.getfilelist())  # makeOct combines all of the ground, sky and object fil|es into a .oct file\n",
    "    frontscan, backscan = analysis.moduleAnalysis(scene, sensorsy=sensorsy)\n",
    "    analysis.analysis(octfile, simulationname+\"_noPavers_\"+str(j), frontscan, backscan)  # compare the back vs front irradiance  \n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Simulation With Pavers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "demo.gendaylit(timeindex)\n",
    "tilt = demo.getSingleTimestampTrackerAngle(metdata, timeindex=timeindex, gcr=gcr, \n",
    "                                   azimuth=180, axis_tilt=0, \n",
    "                                   limit_angle=60, backtrack=True)\n",
    "# create a scene with all the variables\n",
    "sceneDict = {'tilt':tilt,'gcr': gcr,'hub_height':hub_height,'azimuth':azimuth_ang, 'module_type':moduletype, 'nMods': nMods, 'nRows': nRows}  \n",
    "scene = demo.makeScene(module=mymodule, sceneDict=sceneDict) #makeScene creates a .rad file with 20 modules per row, 7 rows."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Custom Object Name objects\\PAVER0.rad\n",
      "\n",
      "Custom Object Name objects\\PAVER1.rad\n",
      "\n",
      "Custom Object Name objects\\PAVER2.rad\n"
     ]
    }
   ],
   "source": [
    "torquetubelength = demo.module.scenex*(nMods) \n",
    "pitch = demo.module.sceney/gcr\n",
    "startpitch = -pitch * (nRows-1)/2\n",
    "p_w = 0.947 # m\n",
    "p_h = 0.092 # m\n",
    "p_w2 = 0.187 # m\n",
    "p_h2 = 0.184 # m\n",
    "offset_w1y = -(p_w/2)+(p_w2/2)\n",
    "offset_w2y = (p_w/2)-(p_w2/2)\n",
    "\n",
    "customObjects = []\n",
    "for i in range (0, nRows):    \n",
    "    name='PAVER'+str(i)\n",
    "    text='! genbox {} paver{} {} {} {} | xform -t {} {} 0 | xform -t {} 0 0'.format(materialpav, i, \n",
    "                                    p_w, torquetubelength, p_h, \n",
    "                                    -p_w/2, (-torquetubelength+demo.module.sceney)/2.0,\n",
    "                                    startpitch+pitch*i)\n",
    "    text += '\\r\\n! genbox {} paverS1{} {} {} {} | xform -t {} {} 0 | xform -t {} 0 0'.format(materialpav, i,\n",
    "                                    p_w2, torquetubelength, p_h2, \n",
    "                                    -p_w2/2+offset_w1y, (-torquetubelength+demo.module.sceney)/2.0,\n",
    "                                    startpitch+pitch*i)\n",
    "    text += '\\r\\n! genbox {} paverS2{} {} {} {} | xform -t {} {} 0 | xform -t {} 0 0'.format(materialpav, i,\n",
    "                                    p_w2, torquetubelength, p_h2, \n",
    "                                    -p_w2/2+offset_w2y, (-torquetubelength+demo.module.sceney)/2.0,\n",
    "                                    startpitch+pitch*i)\n",
    "\n",
    "    customObject = demo.makeCustomObject(name,text)\n",
    "    customObjects.append(customObject)\n",
    "    demo.appendtoScene(radfile=scene.radfiles, customObject=customObject, text=\"!xform -rz 0\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Created tutorial_14.oct\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "'tutorial_14.oct'"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "demo.makeOct()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You can view the geometry generated in the terminal with:\n",
    "\n",
    "### rvu -vf views\\front.vp -e .01 -pe 0.01 -vp -5 -14 1 -vd 0 0.9946 -0.1040 tutorial_14.oct"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "## Comment the ! line below to run rvu from the Jupyter notebook instead of your terminal.\n",
    "## Simulation will stop until you close the rvu window\n",
    "\n",
    "#!rvu -vf views\\front.vp -e .01 -pe 0.01 -vp -5 -14 1 -vd 0 0.9946 -0.1040 tutorial_14.oct\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Linescan in process: tutorial_14_WITHPavers_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers.csv\n",
      "Simulation WITH Pavers Finished\n"
     ]
    }
   ],
   "source": [
    "analysis = AnalysisObj(octfile, demo.name)  # return an analysis object including the scan dimensions for back irradiance\n",
    "frontscan, backscan = analysis.moduleAnalysis(scene, sensorsy=sensorsy)\n",
    "analysis.analysis(octfile, simulationname+\"_WITHPavers\", frontscan, backscan)  # compare the back vs front irradiance  \n",
    "print(\"Simulation WITH Pavers Finished\")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## LOOP WITH PAVERS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_1_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_1_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_1.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_2_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_2_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_2.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_3_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_3_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_3.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_4_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_4_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_4.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_5_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_5_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_5.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_6_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_6_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_6.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_7_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_7_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_7.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_8_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_8_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_8.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_9_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_9_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_9.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_10_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_10_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_10.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_11_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_11_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_11.csv\n",
      "Created tutorial_14.oct\n",
      "Linescan in process: tutorial_14_WITHPavers_12_Front\n",
      "Linescan in process: tutorial_14_WITHPavers_12_Back\n",
      "Saved: results\\irr_tutorial_14_WITHPavers_12.csv\n"
     ]
    }
   ],
   "source": [
    "j=0\n",
    "for timess in range (starttimeindex, endtimeindex):\n",
    "    j+=1\n",
    "    demo.gendaylit(timess)\n",
    "    tilt = demo.getSingleTimestampTrackerAngle(metdata, timeindex=timess, gcr=gcr, \n",
    "                                       azimuth=180, axis_tilt=0, \n",
    "                                       limit_angle=60, backtrack=True)\n",
    "    # create a scene with all the variables\n",
    "    sceneDict = {'tilt':tilt,'gcr': gcr,'hub_height':hub_height,'azimuth':azimuth_ang, 'module_type':moduletype, 'nMods': nMods, 'nRows': nRows}  \n",
    "    scene = demo.makeScene(mymodule, sceneDict=sceneDict) #makeScene creates a .rad file with 20 modules per row, 7 rows.\n",
    "    # Appending Pavers here\n",
    "    demo.appendtoScene(radfile=scene.radfiles, customObject=customObjects[0], text=\"!xform -rz 0\")\n",
    "    demo.appendtoScene(radfile=scene.radfiles, customObject=customObjects[1], text=\"!xform -rz 0\")\n",
    "    demo.appendtoScene(radfile=scene.radfiles, customObject=customObjects[2], text=\"!xform -rz 0\")\n",
    "    octfile = demo.makeOct(demo.getfilelist())  # makeOct combines all of the ground, sky and object fil|es into a .oct file\n",
    "    frontscan, backscan = analysis.moduleAnalysis(scene, sensorsy=sensorsy)\n",
    "    analysis.analysis(octfile, simulationname+\"_WITHPavers_\"+str(j), frontscan, backscan)  # compare the back vs front irradiance  \n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# RESULTS ANALYSIS NOON"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_0 = load.read1Result(os.path.join(testfolder, 'results', 'irr_tutorial_14_noPavers.csv'))\n",
    "df_w = load.read1Result(os.path.join(testfolder, 'results', 'irr_tutorial_14_WITHPavers.csv'))                        "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>z</th>\n",
       "      <th>rearZ</th>\n",
       "      <th>mattype</th>\n",
       "      <th>rearMat</th>\n",
       "      <th>Wm2Front</th>\n",
       "      <th>Wm2Back</th>\n",
       "      <th>Back/FrontRatio</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0.772445</td>\n",
       "      <td>4.729862e-17</td>\n",
       "      <td>1.137844</td>\n",
       "      <td>1.116107</td>\n",
       "      <td>a9.1.a2.0.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.0.0.cellPVmodule.2310</td>\n",
       "      <td>944.8593</td>\n",
       "      <td>116.66480</td>\n",
       "      <td>0.123473</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0.601509</td>\n",
       "      <td>3.683183e-17</td>\n",
       "      <td>1.164490</td>\n",
       "      <td>1.142752</td>\n",
       "      <td>a9.1.a2.1.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.1.0.cellPVmodule.2310</td>\n",
       "      <td>944.9186</td>\n",
       "      <td>108.86900</td>\n",
       "      <td>0.115215</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>0.430574</td>\n",
       "      <td>2.636503e-17</td>\n",
       "      <td>1.191136</td>\n",
       "      <td>1.169398</td>\n",
       "      <td>a9.1.a2.2.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.2.0.cellPVmodule.2310</td>\n",
       "      <td>944.9760</td>\n",
       "      <td>101.76430</td>\n",
       "      <td>0.107690</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>0.259638</td>\n",
       "      <td>1.589824e-17</td>\n",
       "      <td>1.217781</td>\n",
       "      <td>1.196044</td>\n",
       "      <td>a9.1.a2.3.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.3.0.cellPVmodule.2310</td>\n",
       "      <td>945.0333</td>\n",
       "      <td>97.60841</td>\n",
       "      <td>0.103286</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0.088702</td>\n",
       "      <td>5.431448e-18</td>\n",
       "      <td>1.244427</td>\n",
       "      <td>1.222689</td>\n",
       "      <td>a9.1.a2.4.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.4.0.cellPVmodule.2310</td>\n",
       "      <td>945.0907</td>\n",
       "      <td>96.65984</td>\n",
       "      <td>0.102276</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>-0.082233</td>\n",
       "      <td>-5.035345e-18</td>\n",
       "      <td>1.271072</td>\n",
       "      <td>1.249335</td>\n",
       "      <td>a9.1.a2.5.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.5.0.cellPVmodule.2310</td>\n",
       "      <td>945.1480</td>\n",
       "      <td>95.05819</td>\n",
       "      <td>0.100575</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>-0.253169</td>\n",
       "      <td>-1.550214e-17</td>\n",
       "      <td>1.297718</td>\n",
       "      <td>1.275980</td>\n",
       "      <td>a9.1.a2.6.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.6.0.cellPVmodule.2310</td>\n",
       "      <td>945.2054</td>\n",
       "      <td>96.83778</td>\n",
       "      <td>0.102451</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>-0.424105</td>\n",
       "      <td>-2.596893e-17</td>\n",
       "      <td>1.324363</td>\n",
       "      <td>1.302626</td>\n",
       "      <td>a9.1.a2.7.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.7.0.cellPVmodule.2310</td>\n",
       "      <td>945.2628</td>\n",
       "      <td>102.45900</td>\n",
       "      <td>0.108392</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>-0.595040</td>\n",
       "      <td>-3.643572e-17</td>\n",
       "      <td>1.351009</td>\n",
       "      <td>1.329271</td>\n",
       "      <td>a9.1.a2.8.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.8.0.cellPVmodule.2310</td>\n",
       "      <td>946.0780</td>\n",
       "      <td>105.59260</td>\n",
       "      <td>0.111611</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>-0.765976</td>\n",
       "      <td>-4.690252e-17</td>\n",
       "      <td>1.377654</td>\n",
       "      <td>1.355917</td>\n",
       "      <td>a9.1.a2.9.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.9.0.cellPVmodule.2310</td>\n",
       "      <td>946.0906</td>\n",
       "      <td>111.11640</td>\n",
       "      <td>0.117448</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          x             y         z     rearZ                        mattype  \\\n",
       "0  0.772445  4.729862e-17  1.137844  1.116107  a9.1.a2.0.0.cellPVmodule.6457   \n",
       "1  0.601509  3.683183e-17  1.164490  1.142752  a9.1.a2.1.0.cellPVmodule.6457   \n",
       "2  0.430574  2.636503e-17  1.191136  1.169398  a9.1.a2.2.0.cellPVmodule.6457   \n",
       "3  0.259638  1.589824e-17  1.217781  1.196044  a9.1.a2.3.0.cellPVmodule.6457   \n",
       "4  0.088702  5.431448e-18  1.244427  1.222689  a9.1.a2.4.0.cellPVmodule.6457   \n",
       "5 -0.082233 -5.035345e-18  1.271072  1.249335  a9.1.a2.5.0.cellPVmodule.6457   \n",
       "6 -0.253169 -1.550214e-17  1.297718  1.275980  a9.1.a2.6.0.cellPVmodule.6457   \n",
       "7 -0.424105 -2.596893e-17  1.324363  1.302626  a9.1.a2.7.0.cellPVmodule.6457   \n",
       "8 -0.595040 -3.643572e-17  1.351009  1.329271  a9.1.a2.8.0.cellPVmodule.6457   \n",
       "9 -0.765976 -4.690252e-17  1.377654  1.355917  a9.1.a2.9.0.cellPVmodule.6457   \n",
       "\n",
       "                         rearMat  Wm2Front    Wm2Back  Back/FrontRatio  \n",
       "0  a9.1.a2.0.0.cellPVmodule.2310  944.8593  116.66480         0.123473  \n",
       "1  a9.1.a2.1.0.cellPVmodule.2310  944.9186  108.86900         0.115215  \n",
       "2  a9.1.a2.2.0.cellPVmodule.2310  944.9760  101.76430         0.107690  \n",
       "3  a9.1.a2.3.0.cellPVmodule.2310  945.0333   97.60841         0.103286  \n",
       "4  a9.1.a2.4.0.cellPVmodule.2310  945.0907   96.65984         0.102276  \n",
       "5  a9.1.a2.5.0.cellPVmodule.2310  945.1480   95.05819         0.100575  \n",
       "6  a9.1.a2.6.0.cellPVmodule.2310  945.2054   96.83778         0.102451  \n",
       "7  a9.1.a2.7.0.cellPVmodule.2310  945.2628  102.45900         0.108392  \n",
       "8  a9.1.a2.8.0.cellPVmodule.2310  946.0780  105.59260         0.111611  \n",
       "9  a9.1.a2.9.0.cellPVmodule.2310  946.0906  111.11640         0.117448  "
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>z</th>\n",
       "      <th>rearZ</th>\n",
       "      <th>mattype</th>\n",
       "      <th>rearMat</th>\n",
       "      <th>Wm2Front</th>\n",
       "      <th>Wm2Back</th>\n",
       "      <th>Back/FrontRatio</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0.772445</td>\n",
       "      <td>4.729862e-17</td>\n",
       "      <td>1.137844</td>\n",
       "      <td>1.116107</td>\n",
       "      <td>a9.1.a2.0.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.0.0.cellPVmodule.2310</td>\n",
       "      <td>944.7910</td>\n",
       "      <td>127.7906</td>\n",
       "      <td>0.135258</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0.601509</td>\n",
       "      <td>3.683183e-17</td>\n",
       "      <td>1.164490</td>\n",
       "      <td>1.142752</td>\n",
       "      <td>a9.1.a2.1.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.1.0.cellPVmodule.2310</td>\n",
       "      <td>944.8566</td>\n",
       "      <td>121.8282</td>\n",
       "      <td>0.128938</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>0.430574</td>\n",
       "      <td>2.636503e-17</td>\n",
       "      <td>1.191136</td>\n",
       "      <td>1.169398</td>\n",
       "      <td>a9.1.a2.2.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.2.0.cellPVmodule.2310</td>\n",
       "      <td>944.9208</td>\n",
       "      <td>115.7133</td>\n",
       "      <td>0.122458</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>0.259638</td>\n",
       "      <td>1.589824e-17</td>\n",
       "      <td>1.217781</td>\n",
       "      <td>1.196044</td>\n",
       "      <td>a9.1.a2.3.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.3.0.cellPVmodule.2310</td>\n",
       "      <td>944.9850</td>\n",
       "      <td>111.3911</td>\n",
       "      <td>0.117876</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0.088702</td>\n",
       "      <td>5.431448e-18</td>\n",
       "      <td>1.244427</td>\n",
       "      <td>1.222689</td>\n",
       "      <td>a9.1.a2.4.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.4.0.cellPVmodule.2310</td>\n",
       "      <td>945.0492</td>\n",
       "      <td>105.4847</td>\n",
       "      <td>0.111618</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>-0.082233</td>\n",
       "      <td>-5.035345e-18</td>\n",
       "      <td>1.271072</td>\n",
       "      <td>1.249335</td>\n",
       "      <td>a9.1.a2.5.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.5.0.cellPVmodule.2310</td>\n",
       "      <td>945.4636</td>\n",
       "      <td>109.0423</td>\n",
       "      <td>0.115332</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>-0.253169</td>\n",
       "      <td>-1.550214e-17</td>\n",
       "      <td>1.297718</td>\n",
       "      <td>1.275980</td>\n",
       "      <td>a9.1.a2.6.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.6.0.cellPVmodule.2310</td>\n",
       "      <td>945.4807</td>\n",
       "      <td>109.1814</td>\n",
       "      <td>0.115477</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>-0.424105</td>\n",
       "      <td>-2.596893e-17</td>\n",
       "      <td>1.324363</td>\n",
       "      <td>1.302626</td>\n",
       "      <td>a9.1.a2.7.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.7.0.cellPVmodule.2310</td>\n",
       "      <td>945.4962</td>\n",
       "      <td>110.7382</td>\n",
       "      <td>0.117122</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>-0.595040</td>\n",
       "      <td>-3.643572e-17</td>\n",
       "      <td>1.351009</td>\n",
       "      <td>1.329271</td>\n",
       "      <td>a9.1.a2.8.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.8.0.cellPVmodule.2310</td>\n",
       "      <td>945.5118</td>\n",
       "      <td>114.4013</td>\n",
       "      <td>0.120994</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>-0.765976</td>\n",
       "      <td>-4.690252e-17</td>\n",
       "      <td>1.377654</td>\n",
       "      <td>1.355917</td>\n",
       "      <td>a9.1.a2.9.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.9.0.cellPVmodule.2310</td>\n",
       "      <td>945.5274</td>\n",
       "      <td>118.3969</td>\n",
       "      <td>0.125218</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          x             y         z     rearZ                        mattype  \\\n",
       "0  0.772445  4.729862e-17  1.137844  1.116107  a9.1.a2.0.0.cellPVmodule.6457   \n",
       "1  0.601509  3.683183e-17  1.164490  1.142752  a9.1.a2.1.0.cellPVmodule.6457   \n",
       "2  0.430574  2.636503e-17  1.191136  1.169398  a9.1.a2.2.0.cellPVmodule.6457   \n",
       "3  0.259638  1.589824e-17  1.217781  1.196044  a9.1.a2.3.0.cellPVmodule.6457   \n",
       "4  0.088702  5.431448e-18  1.244427  1.222689  a9.1.a2.4.0.cellPVmodule.6457   \n",
       "5 -0.082233 -5.035345e-18  1.271072  1.249335  a9.1.a2.5.0.cellPVmodule.6457   \n",
       "6 -0.253169 -1.550214e-17  1.297718  1.275980  a9.1.a2.6.0.cellPVmodule.6457   \n",
       "7 -0.424105 -2.596893e-17  1.324363  1.302626  a9.1.a2.7.0.cellPVmodule.6457   \n",
       "8 -0.595040 -3.643572e-17  1.351009  1.329271  a9.1.a2.8.0.cellPVmodule.6457   \n",
       "9 -0.765976 -4.690252e-17  1.377654  1.355917  a9.1.a2.9.0.cellPVmodule.6457   \n",
       "\n",
       "                         rearMat  Wm2Front   Wm2Back  Back/FrontRatio  \n",
       "0  a9.1.a2.0.0.cellPVmodule.2310  944.7910  127.7906         0.135258  \n",
       "1  a9.1.a2.1.0.cellPVmodule.2310  944.8566  121.8282         0.128938  \n",
       "2  a9.1.a2.2.0.cellPVmodule.2310  944.9208  115.7133         0.122458  \n",
       "3  a9.1.a2.3.0.cellPVmodule.2310  944.9850  111.3911         0.117876  \n",
       "4  a9.1.a2.4.0.cellPVmodule.2310  945.0492  105.4847         0.111618  \n",
       "5  a9.1.a2.5.0.cellPVmodule.2310  945.4636  109.0423         0.115332  \n",
       "6  a9.1.a2.6.0.cellPVmodule.2310  945.4807  109.1814         0.115477  \n",
       "7  a9.1.a2.7.0.cellPVmodule.2310  945.4962  110.7382         0.117122  \n",
       "8  a9.1.a2.8.0.cellPVmodule.2310  945.5118  114.4013         0.120994  \n",
       "9  a9.1.a2.9.0.cellPVmodule.2310  945.5274  118.3969         0.125218  "
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_w"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Improvement in Rear Irradiance"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "10.8"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "round((df_w['Wm2Back'].mean()-df_0['Wm2Back'].mean())*100/df_0['Wm2Back'].mean(),1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# RESULT ANALYSIS DAY"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_0 = load.read1Result(os.path.join(testfolder, 'results', 'irr_tutorial_14_noPavers_1.csv'))\n",
    "df_w = load.read1Result(os.path.join(testfolder, 'results', 'irr_tutorial_14_WITHPavers_1.csv'))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>z</th>\n",
       "      <th>rearZ</th>\n",
       "      <th>mattype</th>\n",
       "      <th>rearMat</th>\n",
       "      <th>Wm2Front</th>\n",
       "      <th>Wm2Back</th>\n",
       "      <th>Back/FrontRatio</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0.407436</td>\n",
       "      <td>2.494829e-17</td>\n",
       "      <td>0.573299</td>\n",
       "      <td>0.562299</td>\n",
       "      <td>a9.1.a2.0.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.0.0.cellPVmodule.2310</td>\n",
       "      <td>579.5359</td>\n",
       "      <td>27.76016</td>\n",
       "      <td>0.047901</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0.320937</td>\n",
       "      <td>1.965169e-17</td>\n",
       "      <td>0.723122</td>\n",
       "      <td>0.712122</td>\n",
       "      <td>a9.1.a2.1.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.1.0.cellPVmodule.2310</td>\n",
       "      <td>583.5565</td>\n",
       "      <td>26.74639</td>\n",
       "      <td>0.045833</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>0.234436</td>\n",
       "      <td>1.435510e-17</td>\n",
       "      <td>0.872944</td>\n",
       "      <td>0.861944</td>\n",
       "      <td>a9.1.a2.2.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.2.0.cellPVmodule.2310</td>\n",
       "      <td>591.8787</td>\n",
       "      <td>26.89933</td>\n",
       "      <td>0.045447</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>0.147936</td>\n",
       "      <td>9.058500e-18</td>\n",
       "      <td>1.022766</td>\n",
       "      <td>1.011766</td>\n",
       "      <td>a9.1.a2.3.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.3.0.cellPVmodule.2310</td>\n",
       "      <td>596.0308</td>\n",
       "      <td>27.00327</td>\n",
       "      <td>0.045305</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0.061437</td>\n",
       "      <td>3.761903e-18</td>\n",
       "      <td>1.172589</td>\n",
       "      <td>1.161589</td>\n",
       "      <td>a9.1.a2.4.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.4.0.cellPVmodule.2310</td>\n",
       "      <td>601.7504</td>\n",
       "      <td>27.09849</td>\n",
       "      <td>0.045033</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>-0.025063</td>\n",
       "      <td>-1.534695e-18</td>\n",
       "      <td>1.322411</td>\n",
       "      <td>1.311411</td>\n",
       "      <td>a9.1.a2.5.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.5.0.cellPVmodule.2310</td>\n",
       "      <td>605.6486</td>\n",
       "      <td>26.53406</td>\n",
       "      <td>0.043811</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>-0.111563</td>\n",
       "      <td>-6.831292e-18</td>\n",
       "      <td>1.472234</td>\n",
       "      <td>1.461234</td>\n",
       "      <td>a9.1.a2.6.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.6.0.cellPVmodule.2310</td>\n",
       "      <td>612.0669</td>\n",
       "      <td>27.14638</td>\n",
       "      <td>0.044352</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>-0.198064</td>\n",
       "      <td>-1.212789e-17</td>\n",
       "      <td>1.622056</td>\n",
       "      <td>1.611056</td>\n",
       "      <td>a9.1.a2.7.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.7.0.cellPVmodule.2310</td>\n",
       "      <td>615.5109</td>\n",
       "      <td>27.55200</td>\n",
       "      <td>0.044763</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>-0.284564</td>\n",
       "      <td>-1.742449e-17</td>\n",
       "      <td>1.771878</td>\n",
       "      <td>1.760878</td>\n",
       "      <td>a9.1.a2.8.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.8.0.cellPVmodule.2310</td>\n",
       "      <td>619.5745</td>\n",
       "      <td>28.53212</td>\n",
       "      <td>0.046051</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>-0.371063</td>\n",
       "      <td>-2.272108e-17</td>\n",
       "      <td>1.921701</td>\n",
       "      <td>1.910701</td>\n",
       "      <td>a9.1.a2.9.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.9.0.cellPVmodule.2310</td>\n",
       "      <td>622.4969</td>\n",
       "      <td>29.04205</td>\n",
       "      <td>0.046654</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          x             y         z     rearZ                        mattype  \\\n",
       "0  0.407436  2.494829e-17  0.573299  0.562299  a9.1.a2.0.0.cellPVmodule.6457   \n",
       "1  0.320937  1.965169e-17  0.723122  0.712122  a9.1.a2.1.0.cellPVmodule.6457   \n",
       "2  0.234436  1.435510e-17  0.872944  0.861944  a9.1.a2.2.0.cellPVmodule.6457   \n",
       "3  0.147936  9.058500e-18  1.022766  1.011766  a9.1.a2.3.0.cellPVmodule.6457   \n",
       "4  0.061437  3.761903e-18  1.172589  1.161589  a9.1.a2.4.0.cellPVmodule.6457   \n",
       "5 -0.025063 -1.534695e-18  1.322411  1.311411  a9.1.a2.5.0.cellPVmodule.6457   \n",
       "6 -0.111563 -6.831292e-18  1.472234  1.461234  a9.1.a2.6.0.cellPVmodule.6457   \n",
       "7 -0.198064 -1.212789e-17  1.622056  1.611056  a9.1.a2.7.0.cellPVmodule.6457   \n",
       "8 -0.284564 -1.742449e-17  1.771878  1.760878  a9.1.a2.8.0.cellPVmodule.6457   \n",
       "9 -0.371063 -2.272108e-17  1.921701  1.910701  a9.1.a2.9.0.cellPVmodule.6457   \n",
       "\n",
       "                         rearMat  Wm2Front   Wm2Back  Back/FrontRatio  \n",
       "0  a9.1.a2.0.0.cellPVmodule.2310  579.5359  27.76016         0.047901  \n",
       "1  a9.1.a2.1.0.cellPVmodule.2310  583.5565  26.74639         0.045833  \n",
       "2  a9.1.a2.2.0.cellPVmodule.2310  591.8787  26.89933         0.045447  \n",
       "3  a9.1.a2.3.0.cellPVmodule.2310  596.0308  27.00327         0.045305  \n",
       "4  a9.1.a2.4.0.cellPVmodule.2310  601.7504  27.09849         0.045033  \n",
       "5  a9.1.a2.5.0.cellPVmodule.2310  605.6486  26.53406         0.043811  \n",
       "6  a9.1.a2.6.0.cellPVmodule.2310  612.0669  27.14638         0.044352  \n",
       "7  a9.1.a2.7.0.cellPVmodule.2310  615.5109  27.55200         0.044763  \n",
       "8  a9.1.a2.8.0.cellPVmodule.2310  619.5745  28.53212         0.046051  \n",
       "9  a9.1.a2.9.0.cellPVmodule.2310  622.4969  29.04205         0.046654  "
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_w"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>z</th>\n",
       "      <th>rearZ</th>\n",
       "      <th>mattype</th>\n",
       "      <th>rearMat</th>\n",
       "      <th>Wm2Front</th>\n",
       "      <th>Wm2Back</th>\n",
       "      <th>Back/FrontRatio</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0.407436</td>\n",
       "      <td>2.494829e-17</td>\n",
       "      <td>0.573299</td>\n",
       "      <td>0.562299</td>\n",
       "      <td>a9.1.a2.0.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.0.0.cellPVmodule.2310</td>\n",
       "      <td>579.1415</td>\n",
       "      <td>20.88081</td>\n",
       "      <td>0.036055</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0.320937</td>\n",
       "      <td>1.965169e-17</td>\n",
       "      <td>0.723122</td>\n",
       "      <td>0.712122</td>\n",
       "      <td>a9.1.a2.1.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.1.0.cellPVmodule.2310</td>\n",
       "      <td>582.8452</td>\n",
       "      <td>21.34479</td>\n",
       "      <td>0.036622</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>0.234436</td>\n",
       "      <td>1.435510e-17</td>\n",
       "      <td>0.872944</td>\n",
       "      <td>0.861944</td>\n",
       "      <td>a9.1.a2.2.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.2.0.cellPVmodule.2310</td>\n",
       "      <td>591.8007</td>\n",
       "      <td>20.74030</td>\n",
       "      <td>0.035046</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>0.147936</td>\n",
       "      <td>9.058500e-18</td>\n",
       "      <td>1.022766</td>\n",
       "      <td>1.011766</td>\n",
       "      <td>a9.1.a2.3.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.3.0.cellPVmodule.2310</td>\n",
       "      <td>595.6917</td>\n",
       "      <td>21.40979</td>\n",
       "      <td>0.035941</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0.061437</td>\n",
       "      <td>3.761903e-18</td>\n",
       "      <td>1.172589</td>\n",
       "      <td>1.161589</td>\n",
       "      <td>a9.1.a2.4.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.4.0.cellPVmodule.2310</td>\n",
       "      <td>602.5026</td>\n",
       "      <td>22.87528</td>\n",
       "      <td>0.037967</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>-0.025063</td>\n",
       "      <td>-1.534695e-18</td>\n",
       "      <td>1.322411</td>\n",
       "      <td>1.311411</td>\n",
       "      <td>a9.1.a2.5.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.5.0.cellPVmodule.2310</td>\n",
       "      <td>606.3221</td>\n",
       "      <td>23.75233</td>\n",
       "      <td>0.039174</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>-0.111563</td>\n",
       "      <td>-6.831292e-18</td>\n",
       "      <td>1.472234</td>\n",
       "      <td>1.461234</td>\n",
       "      <td>a9.1.a2.6.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.6.0.cellPVmodule.2310</td>\n",
       "      <td>612.8148</td>\n",
       "      <td>24.00990</td>\n",
       "      <td>0.039180</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>-0.198064</td>\n",
       "      <td>-1.212789e-17</td>\n",
       "      <td>1.622056</td>\n",
       "      <td>1.611056</td>\n",
       "      <td>a9.1.a2.7.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.7.0.cellPVmodule.2310</td>\n",
       "      <td>616.6312</td>\n",
       "      <td>25.03161</td>\n",
       "      <td>0.040594</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>-0.284564</td>\n",
       "      <td>-1.742449e-17</td>\n",
       "      <td>1.771878</td>\n",
       "      <td>1.760878</td>\n",
       "      <td>a9.1.a2.8.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.8.0.cellPVmodule.2310</td>\n",
       "      <td>619.6736</td>\n",
       "      <td>26.37031</td>\n",
       "      <td>0.042555</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>-0.371063</td>\n",
       "      <td>-2.272108e-17</td>\n",
       "      <td>1.921701</td>\n",
       "      <td>1.910701</td>\n",
       "      <td>a9.1.a2.9.0.cellPVmodule.6457</td>\n",
       "      <td>a9.1.a2.9.0.cellPVmodule.2310</td>\n",
       "      <td>622.9435</td>\n",
       "      <td>27.44040</td>\n",
       "      <td>0.044050</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          x             y         z     rearZ                        mattype  \\\n",
       "0  0.407436  2.494829e-17  0.573299  0.562299  a9.1.a2.0.0.cellPVmodule.6457   \n",
       "1  0.320937  1.965169e-17  0.723122  0.712122  a9.1.a2.1.0.cellPVmodule.6457   \n",
       "2  0.234436  1.435510e-17  0.872944  0.861944  a9.1.a2.2.0.cellPVmodule.6457   \n",
       "3  0.147936  9.058500e-18  1.022766  1.011766  a9.1.a2.3.0.cellPVmodule.6457   \n",
       "4  0.061437  3.761903e-18  1.172589  1.161589  a9.1.a2.4.0.cellPVmodule.6457   \n",
       "5 -0.025063 -1.534695e-18  1.322411  1.311411  a9.1.a2.5.0.cellPVmodule.6457   \n",
       "6 -0.111563 -6.831292e-18  1.472234  1.461234  a9.1.a2.6.0.cellPVmodule.6457   \n",
       "7 -0.198064 -1.212789e-17  1.622056  1.611056  a9.1.a2.7.0.cellPVmodule.6457   \n",
       "8 -0.284564 -1.742449e-17  1.771878  1.760878  a9.1.a2.8.0.cellPVmodule.6457   \n",
       "9 -0.371063 -2.272108e-17  1.921701  1.910701  a9.1.a2.9.0.cellPVmodule.6457   \n",
       "\n",
       "                         rearMat  Wm2Front   Wm2Back  Back/FrontRatio  \n",
       "0  a9.1.a2.0.0.cellPVmodule.2310  579.1415  20.88081         0.036055  \n",
       "1  a9.1.a2.1.0.cellPVmodule.2310  582.8452  21.34479         0.036622  \n",
       "2  a9.1.a2.2.0.cellPVmodule.2310  591.8007  20.74030         0.035046  \n",
       "3  a9.1.a2.3.0.cellPVmodule.2310  595.6917  21.40979         0.035941  \n",
       "4  a9.1.a2.4.0.cellPVmodule.2310  602.5026  22.87528         0.037967  \n",
       "5  a9.1.a2.5.0.cellPVmodule.2310  606.3221  23.75233         0.039174  \n",
       "6  a9.1.a2.6.0.cellPVmodule.2310  612.8148  24.00990         0.039180  \n",
       "7  a9.1.a2.7.0.cellPVmodule.2310  616.6312  25.03161         0.040594  \n",
       "8  a9.1.a2.8.0.cellPVmodule.2310  619.6736  26.37031         0.042555  \n",
       "9  a9.1.a2.9.0.cellPVmodule.2310  622.9435  27.44040         0.044050  "
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "17.3"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "round((df_w['Wm2Back'].mean()-df_0['Wm2Back'].mean())*100/df_0['Wm2Back'].mean(),1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "17.3\n",
      "39.4\n",
      "33.5\n",
      "22.8\n",
      "8.7\n",
      "9.7\n",
      "10.0\n",
      "9.7\n",
      "17.2\n",
      "32.3\n",
      "37.8\n",
      "39.9\n"
     ]
    }
   ],
   "source": [
    "average_back_d0=[]\n",
    "average_back_dw=[]\n",
    "average_front = []\n",
    "hourly_rearirradiance_comparison = []\n",
    "\n",
    "timessimulated = endtimeindex-starttimeindex\n",
    "\n",
    "for i in range (1, timessimulated+1):\n",
    "    df_0 = load.read1Result(os.path.join(testfolder, 'results', 'irr_tutorial_14_noPavers_'+str(i)+'.csv'))\n",
    "    df_w = load.read1Result(os.path.join(testfolder, 'results', 'irr_tutorial_14_WITHPavers_'+str(i)+'.csv'))\n",
    "    print(round((df_w['Wm2Back'].mean()-df_0['Wm2Back'].mean())*100/df_0['Wm2Back'].mean(),1))\n",
    "    hourly_rearirradiance_comparison.append(round((df_w['Wm2Back'].mean()-df_0['Wm2Back'].mean())*100/df_0['Wm2Back'].mean(),1))\n",
    "    average_back_d0.append(df_0['Wm2Back'].mean())\n",
    "    average_back_dw.append(df_w['Wm2Back'].mean())\n",
    "    average_front.append(df_0['Wm2Front'].mean())\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Increase in rear irradiance:  20.1\n"
     ]
    }
   ],
   "source": [
    "print(\"Increase in rear irradiance: \", round((sum(average_back_dw)-sum(average_back_d0))*100/sum(average_back_d0),1))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "BG no Pavers:  8.3\n",
      "BG with Pavers:  10.0\n"
     ]
    }
   ],
   "source": [
    "print(\"BG no Pavers: \", round(sum(average_back_d0)*100/sum(average_front),1))\n",
    "print(\"BG with Pavers: \", round(sum(average_back_dw)*100/sum(average_front),1))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "\n",
    "#metdata.datetime[starttime].hour # 7\n",
    "#metdata.datetime[endtimeindex].hour # 17\n",
    "xax= [7, 8, 9, 10, 11, 12,13,14,15,16,17,18]  # Lazy way to get the x axis...\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 0, 'Hour')"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEGCAYAAABPdROvAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAA4k0lEQVR4nO3deXycdbX48c/JniZNk7RZuiXd0iXdaaBFFttCK6JQEFCgcBG5wEVA5aoo15+IuFxEES+KIiiCsnhRURBRKW0KlKWlCzRJm7SlS9q0TdKmmSRNs5/fHzPlhtrSSWb5zkzO+/Wa18w8med5zjTpnHm+y/mKqmKMMcb0V5zrAIwxxkQ3SyTGGGMCYonEGGNMQCyRGGOMCYglEmOMMQFJcB1AuA0bNkzHjBnjOgxjjIkq69atO6CqOcf72YBLJGPGjGHt2rWuwzDGmKgiIrtO9DNr2jLGGBMQSyTGGGMCYonEGGNMQCyRGGOMCUhUJRIRiReRDSLygu95togsE5Gtvvss1zEaY8xAE1WJBPgisLnX868Dy1W1CFjue26MMSaMoiaRiMgo4BPAr3ptXgI87nv8OHBRmMMyxpgBL2oSCfAT4Hagp9e2PFXdB+C7z3UQlzHGRLT2rm6eWl3NP8r3h+T4UTEhUUQ+CdSp6joRmd+P/W8AbgAoKCgIbnDGGBOhWju6eHrNbh55dTv7m9q4cOYIzpuWH/TzREUiAc4ALhSR84EUIENEngBqRWS4qu4TkeFA3fF2VtWHgYcBSkpKbCUvY0xMa2rr5Hdv7uLXq3bQcLiDuWOz+eFlMzhzwrCQnC8qEomq3gHcAeC7IvmKql4lIj8ErgHu8d0/5yrGYPO0drKnsZWpI4a4DsUYEyUaDnfw6KodPP7mTprbupg/KYdbFkygZEx2SM8bFYnkQ9wDPCMi1wHVwGWO4wmanyzfwm/f3MUzN85jTmFo/wiMMdGttqmNh1/dzlOrq2nr6ua8qfncvGAC00aG54to1CUSVV0JrPQ9Pgic4zKeUHl3dyPdPcqtT23gxS+eReagJNchGWMizO6GVn7xynv8ce0eulVZMnMEN80fT1He4LDGEXWJZCDo6u5h074mTh83lLW7GvjqHzfy8NVzEBHXoRljIsC2uhZ+vnIbz72zl3gRLi0ZxX+cPZ6CoYOcxGOJJAJtP3CYts4eLisZxTlTcvnu3zbz2Bs7ufaMsa5DM8Y4VF7j4ecrt/H38v0kJ8RxzeljuOHsceQPSXEaV9gTiYikAW2q2h3uc0eLsj0eAKaNHEJRbjpvvneQ/36xkpLCbKaPss53Ywaadbsa+NmKbZRW1TM4OYHPzx/P584Yy9D0ZNehAWFIJCISB1wOLAVOBdqBZBGpB14EHlbVraGOI5qU1XhITYxnfE46IsKPLpvJ+Q+8xi1Pr+eFW89kcEqi6xCNMSGmqry+7SA/K93KW9sbyBqUyFcWT+Tq08cwJDWyPgPCcUVSCryMd/huuar2gLfgIrAAuEdE/qyqT4QhlqhQsddD8YgM4uO8fSJZaUk8cMVsLn/4Lf7rz+U8cPks6y8xJkapKss31/HT0m28u7uRvIxk/t8npnDl3AIGJUVmb0Q4ojpXVTuP3aiqDcCfgD+JSGSlV4e6e5SKvU1cNmfUB7afOiab284t4kcvbeGM8UO5/DSboW9MLOnuUV4s28eDpduo3N/MqKxUvnfxNC6dM4rkhHjX4X2okCeS4yUREflPVf2x7/EkVa0KdRzRYseBw7R2dB93/PdN8yfw1vYG7vprBbMLspiUH94hfsaY4Ovs7uHPG2p4aOV7bD9wmPE5adx32UwunDWCxPjoKIcY1uskEckE7gcmiUgbsBG4Drg2nHFEsvKa/+toP1Z8nPDjz8zk/P95jVueWs/zt5xJalJkf1Mxxhxfe1c3//v2bn75ynZqGo9QPDyDny89hfOm5hMXF11N12FNd6raqKrXAt8GVgNFwLPhjCHSldd4SE6Ioyg3/bg/zx2cwk8+M5tt9S3c9XxFmKMzxgRDT4/yhac3cOdzFeQPSeE3nz2Vv33hTM6fPjzqkgi4m0fyiqq2AescnT9ildV4mDI8g4QPuaQ9s2gYn58/ngdL3+MjE4ayZNbIMEZojAnUz1du458VtXzj/Cn8+1ljo37wjKsGuPtE5FER+bWIfNVRDBGnx9fRPm1kxklfe9u5EykpzOK/ni1jx4HDYYjOGBMMpVV13LdsCxfNGhETSQQcJRJVvVlVPwd8Cch0EUMk2tXQSkt7F9P9KLSWEB/HA1fMJjEhjlueWk97l83vNCbS7TxwmC8+vYEp+Rn896dmxEQSAUeJRES+LCLn4F1bxIb++pT5Otr9LR0/IjOVH106k4q9Tfz3i5WhDM0YE6DWji7+44l1xMUJv7x6TkwNlHHVtPVPYDjwFWCzoxgiTkWNh6T4OCb2oXLnucV5fO6MsTz2xs6QLaNpjAmMqnL7HzeypbaZn14xm9HZboorhoqrRPIZYAmQD7zhKIaIU1bjYVL+YJIS+vZr+frHJzNj1BBu/+O77DnUGqLojDH99chr23lh4z6++rHJnFWU4zqcoHOVSDJV9TLgeuALjmKIKKpKeY2nXwvRJCXE8dMrZtOjcOvTG+js7glBhMaY/nh92wHu+Xsl50/P5z8+Os51OCHhKpF0iMgpgAJpjmKIKLsbjtDU5l9H+/EUDk3jnkums6G6kfte2hLk6Iwx/bHnUCu3PLWe8Tnp3HvpzJjpXD+Wq0TyDbwrGz4M/K+jGCJK2fsz2k8+9PdEPjljBFecVsBDr7zHyqq6YIVmjOmHts5u/uOJdXT1KA//WwnpyZFZcDEYXCWSp4AZgOCtADzgle/1kBAnAdfP+tYFxUzKG8yXn3mX2qa2IEVnjOkLVeW//lxGeU0TP/nMLMYOi+2GF1eJpEJVr1bVzwI/cBRDRCmv8TAxb3DAVT5TEuN5cOlsWju6+dLv36G7R4MUoTHGX799cxfPrq/hS+cWcc6UPNfhhFzYE4mIPAKcKSJfFZHJQEO4Y4g0qkpZjaff/SPHmpA7mLuXTOXN7Qf52YptQTmmMcY/a3Y08J0XNnHulFy+sLDIdThhEfZEoqrXA+fjHfZ7HvBQuGOINDWNR2hs7WRaEJfRvXTOKC6ePZL/Wb6Ft7YfDNpxjTEntt/TxuefXEdB9iB+/JlZUVmAsT9cNW3dpKqvq+pP8JaVH9DeLx0/ov8d7ccSEb5z0TTGDE3ji7/fwMGW9qAd2xjzr9q7vJ3rRzq6+eXVc8gYQEtihzWRiEimiPwGuExEPi8iZwBf82O/FBFZIyLvikiFiHzbt/0uEakRkXd8t/ND/R5Cobymifg4Ycrw4CUSgPTkBH565WwOtXby5T+8S4/1lxgTMnc9v4l3djfyo8tmUtSH6hSxwNV6JHfhXY9kIv6tR9IOLFTVmcAs4DwRmef72f2qOst3ezEEYYdcWY2Hotx0UhKDX3tn6oghfPMTU1hZVc+vVm0P+vGNMfD0mmqeXlPNTfPH8/Hpw12HE3YuOtvvBs4GxgGrVPWvJ9tHvVp8TxN9t5j4eh3IjHZ/XTWvkI9Py+fef1SxofpQyM5jzEC0ofoQ33qugrOKhvGVxZNch+OEi872O4EHgGbgEt8orpMSkXgReQeoA5ap6mrfj24RkY2+9U2yTrDvDSKyVkTW1tfXB+FdBM/+pjYOHu4Iav/IsUSEey6ZQf6QFG59egOeI50hO5cxA0l9czs3PbGevCHJ/PSK2cQPkM71Y7nqbL8R2KCq9/hGcZ2Uqnar6ixgFHCaiEwDfgGMx9vctQ+47wT7PqyqJapakpMTWQXTyvZ4O9qnB3HE1vEMSU3kp1fMZr+nja//aSOqMXFBZ4wznd093PzkehqPdPDLq0rIHJTkOiRnXCWSR4GbROSHIjKrLzuqaiOwEjhPVWt9CaYHeAQ4LdiBhlr53ibihKB3tB/P7IIsbj9vEn8v388Tb+0K+fmMiWXf+9tm1uxs4AeXzKA4hC0K0cBVIvkC3vXiE/A2c30oEckRkUzf41TgXKBSRHr3al0MlAc/1NAqr/EwPiedQUnhqcPz72eOY/6kHL7zt81U7PWE5ZzGxJpn1+/hsTd2ct2ZY1kya6TrcJxzlUjew7s64nOqerYfrx8OlIrIRuBtvH0kLwD3ikiZb/sC4LaQRRwi5UGc0e6PuDjhvstmkjUokVuf2sDh9q6wnduYWFBe4+GOZ8uYNy6bOz4+2XU4EcFZrS1gBXCdiLx9sher6kZVna2qM1R1mqre7dt+tapO922/UFX3hTrwYKpraqOuuZ2pYUwkAEPTk/mfy2ez8+Bh/t9fyq2/xBg/NRzu4MbfrSM7LYmfXXkKCfGuPkIji6t/hfF4m7UeBq51FINzR0vHh/OK5Kh544byxXMm8ucNNfxx3Z6wn9+YaNPV3cOtT6+nvqWdh66aw7D0ZNchRQwnBfJV9XEX54005TVNiOCso+6WhRN4a/tB7nyugtkFmUzIHVizcY3pix/+s4rXtx3k3ktnMHN0putwIoqz6zIRGebq3JGirMbD2GFpzha8iY8TfnL5LAYlxXPzkxto6+x2Eocxke6FjXv55avbuWpeAZ8uGe06nIjjsoHvUYfnjggVe8Pb0X48eRkp3PfpmVTVNnP/Mlui15hjVe5v4qt/2Micwizu/ORU1+FEJJeJZGBOAfU50NLOPk8b00a4TSQA8yflcuHMETy1pprWDhvFZcxRntZObvzdOtJTEvj50lNISrDO9eNx+a8yoIcKvV863vEVyVFXn15Ic1sXf313r+tQjIkIPT3Kl/53A3sbj/CLpaeQl5HiOqSIZVckjhxNJFNHRsaM2JLCLCbmpfPk6mrXoRgTEX7y8hZKq+q584KplIzJdh1ORHOZSO5weG7nymo8jBk6KGIWvxERrppXyMY9HjbuaXQdjjFO/bNiPw+s2MZlc0Zx1dwC1+FEPCeJRETyVDXqypkEU3lNU8Q0ax110eyRpCbG8+RbdlViBq5tdS18+Zl3mTFqCN+5aBoiA7rxxC+urkiucXTeiHDocAc1jUciLpFkpCRy0ewRPPdujZWaNwNSc1snN/5uLckJcTx01ZyQLDYXi1wlkiUicouIDMhVYMr3upvRfjJL5xbS1tnDn9fbbHczsKgqX37mXXYebOVnV57CiMxU1yFFDVeJ5FPANuBiEfmVoxicOVoaZWoElp6eNnIIM0dn8sTqaqvBZQaU1TsaeGlTLbd/bBKnjx/qOpyo4iSR+NYR+YdvYat/dxGDS+U1HkZnp0bsQjhL5xawra6FNTsaXIdiTNiUVtaRGC9caZ3rfeaqs/0uEfm7iPxKRG52EYNL5TVNEdmsddQFM0aQkZLAEzYU2AwgpVV1nDomm8ERMpIymrhq2soE3gK+BwyofhJPayfVDa1MjYAZ7SeSmhTPJXNG8Y/yfRxoaXcdjjEht+dQK1tqW1gwKdd1KFHJVSJpAOKBOt/jAaMigjvae1s6t5DObuWZtbtdh2JMyK2sqgdgweQcx5FEJ1d9JHcDD+FdZndArfdaFmGlUU5kQm4688Zl89Tqanp6rNPdxLaVVXWMykplfE6661Cikqs+kgeB76vqdcAmFzG4Ur63iZGZqWSnRWZHe29XzStkz6EjvLK13nUoxoRMW2c3r287yMLJuTb5sJ9cNW11ANt9jxc4isGJ8hoP0yKkvtbJLC7OZ1h6ss10NzFtzY4GjnR2W/9IAFwlklZgiIgkAgNmrF1TWyc7DhyOiNLx/khKiOMzp45iRWUtNY1HXIdjTEiUVtWRnBDHvHE2d6S/XCWSXwDvAQ8CTzmKIew27W0CYNqo6EgkAJefWoAC/7vGrkpMbFpZVc/p44eSmmTlUPrLVSK5UlV/rqo3qOoLjmIIu/fXIImSKxKA0dmDWDApl9+/vZvO7h7X4RgTVDsOHGbHgcPWrBUgq7UVRuU1HvIzUsgZnOw6lD5ZOreAuuZ2Xt5U6zoUY4JqZVUdgCWSAEVFrS0RSRGRNSLyrohUiMi3fduzRWSZiGz13WeFOvBAlEVRR3tv8yflMjIz1Ra9MjGntKqecTlpFAwd5DqUqBYttbbagYWqOhOYBZwnIvOArwPLVbUIWO57HpFa2rvYfuBwxM8fOZ74OOGK00azatsBdhw47DocY4KitaOLt7YftKuRIAh5IhGRJb3raYnIahHZ7rtd6s8x1KvF9zTRd1NgCfC4b/vjwEXBizy4Nu9rQjXyZ7SfyKdPHU1CnPDU6l2uQzEmKN587yAdXT2WSIIgHFcktwPP93qeDJwKzAdu8vcgIhIvIu/gLauyTFVXA3mqug/Ad3/cvwgRuUFE1orI2vp6N5PryvZEx4z2E8kdnMLHpubzh3V7aOvsdh2OMQErrapjUFI8p46N6BbxqBCORJKkqr0LNq1S1YOqWg2k+XsQVe1W1VnAKOA0EZnWh30fVtUSVS3JyXFTS6d8r4ecwcnkZaQ4OX8wLJ1XQGNrJy+W7XMdijEBUVVKK+s5Y8IwkhNs2G+gwpFIPpDuVfWWXk/7/Kmuqo3ASuA8oFZEhgP47uv6HWWIldd4mBaBC1n1xenjhjIuJ40n3rLmLRPdttW1UNN4xJq1giQciWS1iFx/7EYRuRFY488BRCRHRDJ9j1OBc4FKvE1mR9d/vwZ4LhgBB9uRjm621bVEbf/IUSLC0rmFrK9ufH9ypTHRqNQ37Hf+JKv2GwzhSCS3AdeKSKmI3Oe7rQQ+C3zJz2MMB0pFZCPwNt4+kheAe4BFIrIVWOR7HnE27WuiR6O3f6S3S04ZSXJCHE9ap7uJYqWV9UzOH2zrsgdJQqhPoKp1wEdEZCEw1bf5b6q6og/H2AjMPs72g8A5QQk0hMqjpHS8PzIHJXHBzBH8ZUMNd5w/hfTkkP8JGRNUzW2dvL2zgX8/a5zrUGLGST8FRCTbj+P0+PouTsiXOPxOHrGkvMbD0LQkhg+J3o723pbOLeCP6/bwlw01XDWv0HU4xvTJ69sO0NWjLLBmraDx5+vkXt/twwr1xzOAqvj2VVmNh6kjh8TMWgezRmcydUQGT66uZuncgph5X2ZgKK2sZ3BKAqcU2rDfYPGnj2Szqo5T1bEnugEHT7SziJwuA/iTpq2zm611LUyPwtIoJ3K0033zvibWVze6DscYv6kqpVV1nF2UQ2K8qwpRsceff8nTA3zNNcA6Efm9iHxWRPL9Cy02VO5vprtHo6rirz+WzBpBenKCdbqbqLJpXxN1ze02WivITppIVLUtkNeo6n+o6inAXXjnlDwmIm+KyPdF5GwRienZQNGyRntfpSUncPHskbywcR+HDne4DscYv6ys8la2+KglkqDq87WdiMwTkRUi8rqIXOTvfqpaqar3q+p5wEJgFXAZsLqvMUSTihoPmYMSGZUVe8MMl84roKOrhz+t3+M6FGP8sqKyjukjh5A7ODYGvkSKkyaS4zRF/SdwId6Z5d/pz0lV9Yiqvqiqt6pqSX+OES3KajxMGxE7He29Tc7PoKQwiydXV9PTo67DMeZDHTrcwYbqQzZaKwT8uSJ5SES+KSJHU3gjcCXwGcCmN3+I9q5uttQ2x1yzVm9L5xWw48Bh3tx+wvEWxkSEV7fW06Mwf7KVRQk2f/pILgLeAV4QkavxzkbvAQYRwWXbI8GW/S10dmtULmblr49PG07WoESrv2Ui3sqqerLTkpg5KtN1KDHHrz4SVf0r8DEgE3gWqFLVB1S1XzXZReQ/ez2O2eV2j3a0R3uNrQ+TkhjPZSWjeWlTLbVNJx2XYYwT3T3KK1vq+ejEHOLjYq+Z2TV/+kguFJFVeGellwOX410i92kRGd+Xk4lIpoj8BrhURD4vImcSwasaBqp8r4fBKQkUZMf2Mp5XnlZAd4/yzNu7T/5iYxzYuKeRhsMdNuw3RPy5Ivku3quRS4AfqGqjqv4ncCfwvb6czLfvtcC38Y7WKsJ7hROTymO4o723McPSOKtoGE+vqabbOt1NBCqtqidO4OwiSySh4E8i8eC9CrmcXut9qOpWVb28n+d9RVXXqepvfM1mMaezu4fKfc1MHxW7zVq9LZ1byF5PG6WVEbskjBnAVlbVMbsgi6y0JNehxCR/EsnFeDvWu/CO1gqG+0TkURH5tYh8NUjHjChbapvp6O5hapQvZuWvc6fkkpeRzBM2091EmPrmdjbu8diw3xA6adFGVT0A/DSYJ1XVmwFEZDAx2kdSUeMdGR3LHe29JcTHcfmpBTywYiu7G1oZHeP9QiZ6vLLFOyZovq2GGDL+dLavD8Zrjnn9l0XkHCAFSOzLvtGirMZDenICY4b6vSx91Lv8tNHEifDUmmrXoRjzvtKqOnIHJw+Y1gEX/CkjP8W3MuGJCNDXr93/BGYBi4HNfdw3KpTVeCgekUHcABpqOHxIKudMzuWZt3dz27kTSUqw6qrGra7uHl7dUs/Hp+XH/KAXl/xJJJP9eE13H8/7Gd9xW4FH+7hvxOvq7mHzvqYBuejT0nmFvLSpln9U7OfCmSNch2MGuPXVjTS3dbHAmrVCyp8+klD0nmaq6mUikgTcD9wcgnM4s62+hfaunpie0X4iZ00YRkH2IJ58a5clEuNcaVUdCXHCGUXDXIcS01y1PXSIyCmAAjHXiVA+wDrae4uLE66cW8DqHQ1srW12HY4Z4Eor6ygZk0VGSkx2xUaMfiUSEckL8LzfAM4BHgb+N8BjRZzyGg+DkuIZOyzddShOXDZnFEnxcTy52jrdjTv7PEeo3N9szVph0N8rkmsCPO9TwAy8HfULAjxWxCmr8VA8PGPA1vQZmp7Mx6fn86f1e2jt6HIdjhmgji5itcCq/YZcfxPJEhG5JYCCixWqerWqfhb4wcleLCKjRaRURDaLSIWIfNG3/S4RqRGRd3y38/sZT9B09yib9jbFdOl4f1w1r5Dmti5eeHef61DMAFVaWcfIzFSKcgdmy0A49TeRfArYBlwkIr/qy44i8ghwpoh8VUQmAw1+7NYFfFlVpwDzgJtFpNj3s/tVdZbv9mJfYgmF7fUtHOnsHvCJpKQwi4l56TbT3TjR3tXN69sOMH9Sjg37DYP+JpKbgC8CE4ENfdlRVa8HzgfewLvK4kN+7LNPVdf7HjfjnXsyso8xh0X53tgvHe8PEeGqeYVs3ONh455G1+GYAWbtzkMc7ui2/pEw6W8iyQTewlsZuD/NWzep6uuq+hO8w3/9JiJjgNn831rvt4jIRl/trqx+xBJUZXuaSEmMY3xOzA1G67OLZo8kNTGeJ9+yTncTXqWVdSTFx/GRCUNdhzIg9DeRNADxeKsB+9M0BXxgPZLLfOuRnAF8rQ/7pwN/Ar6kqk3AL4DxeGfJ7wPuO8F+N4jIWhFZW1/fr7W4/Fa+18OU4RkkxNus7oyURC6aPYLn392L50in63DMAFJaVcfccdkMSvJnzrUJVL8+7VT1brxNUg/gLTPv735H1yO5C+8VxUT8XI9ERBLxJpEnVfVZ3/FqVbVbVXuAR4DTTnDeh1W1RFVLcnJCVwG0x9fRPtCbtXpbOreQI53d/Hn9HtehmAGi+mAr79UftmatMOrvPJIHge+r6nXApj7uezdwNjAOWOXPeiTi7S37NbBZVX/ca/vwXi+7GO8Kjs7sOHiYlvYupo2wRHLUtJFDmDk6kydXV6Nqi16Z0Fu5xbsmjg37DZ/+tr90ANt9j/s0D0RV78R7JdMMXOIbxXUyZwBXAwuPGep7r4iU+YpKLgBu60sswVbuW6N9oI/YOtbSuQVsrWthzQ6/W0GN6bfSyjrGDB3E2GHWTxku/W1AbAWG+JqbCvqx/43AL1X1H/68WFVX4Z28eCznw317K6/xkJQQR1GejVvv7YIZI/juC5t4cnU1c8dZ56cJnbbObt547yBXnNafjyXTX/29IvkF8B7wIN5Z6n31KHCTiPxQRGb1M4aIU17TxJT8wSRaR/sHpCbFc8mcUfy9fB8HWtpdh2Ni2JvbD9Le1WPNWmHW30+8K1X156p6g6q+0I/9v4D3aigBbzNX1FNVyvd6mGrNWse1dG4hnd3KH9Zap7sJnZWVdaQmxjN3bLbrUAYUVyVS3sO7OuJzqnp2P48RUXYdbKW5rctGbJ3AhNx05o3L5qk1u+jpsU53E3yqSmlVPWdMGEpKYrzrcAaUQEukXNzXEik+FcAK4DoRebufMUQUm9F+clfNK2R3wxFe3RrauTxmYHqv/jDVDa22NrsD/epsV9Va4B++W3+MBw7hLSN/qJ/HiChlNR4S48U62j/E4uJ8hqUn8+TqavvPboJuZZV32O/8SaGbK2aO76RXJCKyRERu7vV8tYhs990u7c9JVfVxVX1eVV9TVadzP4KloqaJSfmDSU6wS+oTSUqI4zOnjmL55lr2Nh5xHY6JMaVVdUzMS2dU1iDXoQw4/jRt3Q483+t5MnAqMB9v8cZ+EZGYWftSVSmr8dhERD9cfmoBCvz+7d2uQzExpKW9izU7Gmw2uyP+JJIkVe39v36Vqh5U1WoCWyb30QD2jSh7Dh3Bc6TTJiL6YXT2IBZMyuX3a6rp7O5xHY6JEa9vO0Bnt1qTqSP+JJIPVNRV1Vt6PQ2kMTJmFgk4OqPdOtr9s3RuAXXN7bxSZZ3uJjhWVtWRnpxAyRjnBcAHJH8SyWoRuf7YjSJyI7AmgHPHzBjQshoPCXHCpPzBrkOJCmdPzCEjJYGXNu13HYqJAapKaWU9ZxUNs8nAjvgzaus24C8iciWw3rdtDt6+kosCOHfsXJHsbaIob7CNXfdTYnwcCybnsnxzHd09OmDXtjfBUbm/mf1NbdY/4tBJ07eq1qnqR4DvADt9t7tV9XTfMOD+uiOAfSOGqlJe42HaiAzXoUSVRcV5HDzcwfrqmBj9bRwq9Q37/agN+3XG73kkqroC7yTCgIlIXqwM+93naaPhcAfTR1n/SF98dGIOSfFxLNtUy6ljrJyF6b+VlfVMHZFBXkaK61AGLFcNitc4Om/QlVnp+H4ZnJLI6eOH8lLFflunxPSbp7WTddWHrFnLMVeJJNBaXRGjvMZDnMCUfGva6qtFxXnsPNjKtroW16GYKPXatnq6e5QFk61ZyyVXiSTQWl0Ro7zGQ1HuYFKTrKO9rxYV5wHw0qZAutrMQFZaWU/moERmjbZhvy753UciIsnAJcCY3vv51m/vkyDU6ooI3hntTZw9MWYm6YdVXkYKM0dn8tKmWm5eMMF1OCbK9PQor2yp4+yiHBv551hfrkieA5YAXcDhXrc+E5G7ROTvIvKr3nW8ok1dczsHWtptImIAFhfn8e7uRmqb2lyHYqJM+V4PB1o6rFkrAvSl+u8oVT0vSOfNBN4CfovjddYDUbbHOtoDtag4jx/+s4qXN9eydG6h63BMFCmtrEcEzi6yROJaX65I3hCR6UE6bwMQD9T5Hkel8r0eRKB4uHW091dRbjpjhg5imfWTmD4qrapj5qhMhqYnuw5lwOtLIjkTWC8iVSKyUUTKRGRjf07q61d5CO8yu57+HCMSlNd4GJ+TTlpyv5Z1MYCIsKg4jze2HaSlvct1OCZKHGxp5909jTbsN0L05RPwPLxlTQIe9C8iDwJpqvpZEflYoMdzpazGw+njhroOI+otKs7nkdd28EpVPZ+YMdx1OCYKvLq1HlWsfyRC+LOw1SrfwwqgDCj33Sp89/3RAWz3PV7Qz2M4VdfcRm1Tu/WPBMGcwiyy05JYZkUcjZ9KK+sZlp5kawBFCH9qbZ3pux+sqhm9boNVtb+dA63AEBFJBApO9mIRGS0ipSKyWUQqROSLvu3ZIrJMRLb67sM2mLyipgmwjvZgiI8Tzpmcy4rKOlujxJxUd4/yypZ6Pjoxlzgb9hsRXE1I/AXwHvAg8JQfr+8CvqyqU4B5wM0iUgx8HViuqkXAct/zsDi6BslUK9YYFIuK82hq865yZ8yHeWf3ITxHOq1ZK4K46iW+UlXv9ffFqroP2Od73Cwim4GReOe1zPe97HFgJfC1oEZ6AmU1HsYNS2NwSmI4ThfzzirKISXRW8TxjAk2wdOcWGllPfFxwlk27DdiRF2tLREZA8wGVgN5viRzNNkcdwiHiNwgImtFZG19fXBW5Suv8TDVmrWCJjUpnrOKcqyIozmp0qo65hRmMSTVvsRFin4nEhE5wzf6qj/6VWtLRNKBPwFfUtUmf/dT1YdVtURVS3JyAv8Wc7Clnb2eNqaPtGatYFpUnMdeTxsVe/3+1ZoBprbJ+/dhw34jS5+atkRkFnAl8GmgFpgM9LnESX9qbfk65v8EPKmqz/o214rIcFXdJyLD8U5wDLly3wedjRgJrnMm5xInsGxTrQ1iMMf1SpW3RcH6RyKLP8N/J4rInSJSCfwKOAjMV9W5+DErXUSW9K6nJSKrRWS773apP0GKiAC/Bjar6o97/eh5/m9tk2vw1gMLufc72u3DLqiGpiczpzDLqgGbEyqtqmP4kBQm5Q12HYrpxZ+mrUrgE8ClvuahH6jqTt/P/GnMvh3vB/5RycCpeDvJb/IzzjOAq4GFIvKO73Y+cA+wSES2Aot8z0OuvMZDQfYga6MNgcXF+Wze18TuhlbXoZgI09ndw2tbDzB/Ui7e75YmUviTSC7Bu077MhH5nYhc4Gtm8leSqu7u9XyVqh5U1WogzZ8DqOoqVRVVnaGqs3y3F33HOUdVi3z3YRk7Wr7XYxV/Q+ToGiUvb7arEvNBa3ceoqW9iwW2NnvE8WdC4p9V9TPABLx9GjcCe0TkN4A/vc0fmCSoqrf0ehp1fxGNrR3sbjhibfghMmZYGkW56VbE0fyLlVV1JMaLDQ+PQH6P2lLVw6r6pKp+EpiCtwx8mR+7rhaR64/dKCI3Amv8jjRClL8/o91GbIXK4ql5rN7RQGNrh+tQTARZUVnH3LFDrUhqBOrX8F9VbVDVX6qqP3WybgOu9ZU4uc93Wwl8FvhSf87vUvle3xokNmIrZBYV59Pdo5RWhWUQnokCuxta2VrXwnxr1opIIZ+QqKp1qvoR4Dt4+1p2Aner6um+YcBRpazGw8jMVLLSklyHErNmjBxC7uBka94y71u55eiwX5s/EonCdo2oqiuAFeE6X6hU1FhHe6jFxQnnFufxlw01tHV2k5IY7zok49jKyjoKsgcxbphf43NMmLkqkRKVmto62XmwlemjLJGE2uLiPFo7unnzvYOuQzGOtXV28/p7B1gwKceG/UYoJ4lERPJdnDdQVvE3fE4fP5T05ASbnGhYvaOBts4e5luzVsRydUXyoqPzBsTWIAmf5IR4Pjoxh5c319LTY0UcB7LSyjqSE+JsNdII5iqRROX1aVmNh+FDUhiWnuw6lAFh8dQ86pvbeWdPo+tQjEMrq+r4yPih1lcWwVwlkkccnTcg5Xs9djUSRvMn5ZIQJzZ6awDbceAwOw+22mitCOckkajqz12cNxAt7V3sOHDYRmyF0ZDUROaOy7ZEMoCVVnrnEs2faIkkktmoLT9V1HhQtRnt4bZoSh7b6lrYXt/iOhTjwIrKOsbnpFEwdJDrUMyHsETip1MKs/jbF87ktLHW4RdOi6Z6B/jZVcnAs3ZnA6u2HeDCmSNdh2JOwhKJnxLj45g6YgjpVucnrEZmpjJ1RIYlkgGmu0f51vMVDB+SwvVnj3UdjjkJV/NInhWR34rIYyJyr4sYTPRYVJzHuupD1De3uw7FhMkza3dTsbeJO86fwqAk+/IW6VxdkVSo6r+p6meBHziKwUSJxcX5qMKKSrsqGQg8Rzr54T+rOG1MNhfMGO46HOOHsCcSEXkEOFNEvioik/FjuV4zsE0ZPpiRmanWvDVA/M/LWznU2sGdFxRbSZQoEfZEoqrXA+cDbwDnAQ+FOwYTXUSERcV5vLb1AK0dXa7DMSG0ra6Z3765k8tPLbA5W1HEVdPWTar6uqr+BLjfUQwmiiwuzqO9q4dXtxxwHYoJEVXl23/dRGpSPF9ZPNF1OKYPwppIRCTTt0TvZSLyeRE5A/haOGMw0enUsdkMSU205q0Y9vLmOl7beoDbzp3IUCtDFFXCOhxCVRvxrpb4MeAAMAN4NpwxmOiUGB/Hwsm5rKispau7h4R4G7keS9q7uvnu3zZRlJvO1acXug7H9JGLzva7gbOBccAqVf1ruGMw0WlRcR6HWjtZu+uQ61BMkP161Q52HWzlzguKSbQvCVHHRWf7ncADQDNwiW8UlzEndfbEHJIS4qx5K8bUNrXxsxXbWFScx1lFtiZ7NHKV+m8ENqjqPb5RXCclIo+KSJ2IlPfadpeI1IjIO77b+SGL2DiXnpzAGeOHsmxTLaq2Rkms+MHfK+nqUb75iWLXoZh+cpVIHgVuEpEfisgsP/d5DO9w4WPdr6qzfLeoXDDL+G9RcT7VDa1U1Ta7DsUEwfrqQzy7oYbrzxprhRmjmKtE8gW8Hf0JeJu5TkpVX8UmLw54507xlhNfVmHNW9Gup0e56/kK8jKS+fz8Ca7DMQFwlUjeA1KA51T17ACPdYuIbPQ1fWUd7wUicoOIrBWRtfX19QGezriUm5HC7IJMlm22RBLt/rhuDxv3eLjj41NIs2KoUc1ZrS1gBXCdiLwdwHF+AYwHZgH7gPuO9yJVfVhVS1S1JCfHOvOi3aLiPDbu8bDPc8R1KKafmto6ufeflcwpzGLJrBGuwzEBcpVIxuNt1noYuLa/B1HVWlXtVtUevMv3nhak+EwEW1ycB8DLNnorav10+VYOHu7grgumWj2tGOBqqd3HVfV5VX1NVctPvsfxiUjv0qAXA/0+loke43PSGTssjZcskUSl9+pb+M3rO/n0nNFMH2X1tGKBs5k/IjKsj69/GngTmCQie0TkOuBeESkTkY3AAuC2EIRqIoyIsLg4j7e2H6SprdN1OKaPvvPCJlIT4/nqeZNch2KCxOUU0kf78mJVvUJVh6tqoqqOUtVfq+rVqjpdVWeo6oWqui9UwZrIsqg4j85uZWWVDZ6IJisqa1lZVc8Xzy1imNXTihkuE4k1jJp+m12QxdC0JJvlHkU6unr4zgubGZeTxr+dPsZ1OCaIXCYSm5ps+i0+Tjh3Sh4rK+vo6OpxHY7xw29e38GOA4e585PFJCVYPa1YYlckJmotKs6jub2Lt7YfdB2KOYm65jZ+umIb50zOZf6kXNfhmCBzmUjucHhuEwPOLBpGamK8NW9FgXv/UUV7Vzf/75NWTysWOUkkIpIXyLBfYwBSEuM5q2iYFXGMcO/sbuSP6/bwuTPHMnZYmutwTAi4uiK5xtF5TYxZPDWf/U1tlNV4XIdijuNoPa2cwcncurDIdTgmRFwlkiUicouI2EByE5CFk3OJE6x5K0I9u6GGd3Y38vXzJpNu9bRilqtE8ilgG3CxiPzKUQwmBmSnJVEyJtsSSQRqae/iB/+oZNboTC6ePdJ1OCaEXJVIqVXVf/gWtvp3FzGY2LG4OI/K/c1UH2x1HYrp5acrtlLf3M5dF04lLs4GacYyV53td4nI30XkVyJys4sYTOxYXJwPwEub9juOxBy148BhHl21g0vnjGLW6EzX4ZgQc9W0lQm8BXwPsH4SE5CCoYOYlDfYmrciyHdf2ERyQjy3Wz2tAcFVImkA4oE6bNVDEwSLivN4e2cDhw53uA5lwFtZVcfyyjpuXTiB3MEprsMxYeCqj+Ru4CG8y+zauE0TsMVT8+hRWF5Z5zqUAa2jq4e7X9jE2GFpXHvGWNfhmDBx1UfyIPB9Vb0O2OQiBhNbpo8cQn5GCsusn8Sp3765k+31h/nmJ6dYPa0BxNVvugPY7nu8wFEMJoaICOcW5/LqlgO0dXa7DmdAqm9u539e3sr8STksnJznOhwTRq4SSSswREQSgQJHMZgYs7g4nyOd3azaesB1KAPSj/5ZxZHObr5p9bQGHFeJ5BfAe8CDwFOOYjAxZt64oQxOTrDRWw6U7fHwzLrdXHvGGMbnpLsOx4SZq5oFV6rqvY7ObWJUUkIcH52Uw/LKWrp7lHibBBcWqspdf61gaFoSt55j9bQGIqu1ZWLKouI8DrR08M7uQ65DGTCee2cv63Yd4vaPTSYjJdF1OMYBq7VlYsqCybkkxgsvVVjzVjgcbu/iv/++mRmjhnDpnFGuwzGOWK0tE1MyUhKZN26o9ZOEyYOl26htaudbF1g9rYEs5IlERJb0rqclIqtFZLvvdmmoz28GnkXFeWw/cJhtdS2uQ4lpuw4e5lev7eBTs0cypzDLdTjGoXBckdwOPN/reTJwKjAfuCkM5zcDzLlTvHMYrIhjaH33b5tJiBe+9vHJrkMxjoUjkSSp6u5ez1ep6kFVrQb8XndTRB4VkToRKe+1LVtElonIVt+9fS0yjMhMZfrIIda8FUKvba1n2aZablk4gbwMq6c10IUjkXzgw11Vb+n1NKcPx3kMOO+YbV8HlqtqEbDc99wYFhXn8c7uRuqa2lyHElLdPUpzWydd3T1hO2dndw/f/usmCocO4rozrZ6WCc88ktUicr2qPtJ7o4jcCKzx9yCq+qqIjDlm8xK8TWQAjwMrga/1O1ITMxYV5/HjZVt4eXMdV86N7OIJqkpbZw+NRzrwHOmksbUTz5FOPK2d/7rNd2ts7aSxtYPm9i5UvccZnJLAkNREMgcleu9Tk8jwPc9MTXz/Zxm+nw3xbR+UFI+I/x3lv3tzF9vqWnjk30pITogP0b+KiSbhSCS3AX8RkSuB9b5tc/D2lVwU4LHzVHUfgKruE5Hc471IRG4AbgAoKIjsDxUTHJPzBzM6O5Vlm/aHPZHUNbVR03iERl8y6J0IGo900PSB597XdHzIFUV8nPgSgzcJZKclMW5YGkNSExkyKIn05HhaO7ppbO30HvuIN8ns9zS9n3g6u/WEx0+M9x5/yPvJJqnX4w/epyTGc//LWziraBjnTjnufzczAIUjkdSr6kdEZCEw1bftb6q64ugLRERU9cR/6QFS1YeBhwFKSkpCdh4TOUSERVPyeWL1Llrau0hPDs2fenePUrm/ifW7DrF21yHW7TrEnkNHjvva9OSED3xAT8hN/+AVQq+rh4xeH97pyQl9umI4lqrS2tF9zJVNxweTWa+kV9fcxpbaZjxHOmlu6/qX4yXECd+6oDigmExsCUciKRWRPwHPHZM8koAzgWuAUrx9IH1VKyLDfVcjw/EulGUM4G3eevT1Hby6pZ7zpw8PyjFb2rvYUH2ItTsPsb76EBuqG2lp937Y5mUkU1KYzbVnjPVeMbzfxORNDInxbub/ighpyQmkJScwIjO1T/t2dffQ3Nb1frJpbO0gOy2JCbmDQxStiUbhSCTnAZ8DnhaRsUAjkIJ3hcSXgPtV9Z1+Hvt5vInoHt/9c4EGa2LHqWOyyByUyLJNtf1KJKrKnkNHWOe70li76xBV+5voUYgTmJSfwcWzR1IyJotTCrIYlZUac9/SE+LjyEpLIistyXUoJoKFPJGoahvwc+DnvrLxw4AjqtrYl+OIyNN4O9aHicge4Ft4E8gzInIdUA1cFsTQTZRLiI9j4eRclm+uo7O756RXBJ3dPVTsbfIljgbW7TpEbVM74G2Wml2QyeKFRZSMyWLW6EwGW10pY4AwV/9V1U5gXz/3veIEPzqn/xGZWLe4OI9n19fw9s4GPjJ+2Ad+1tjawXpfM9W6XYd4d08jbZ3eTu9RWanMGzeUksIs5hRmMyl/sFUTNuYEXJWRNyYszirKISkhjpcqasnPSGHtrkPvd4wfLaGSECdMHZHBlacVMqcwi5IxWTbJzpg+sERiYlpacgJnTRjGY2/s5LE3dgIwJDWROYVZXOyrETVzVCapSTYfwpj+skRiYt4tCyeQm5HCjFFDKCnMYnxOulWqNSaILJGYmDe7IIvZBVaGzZhQcbWwlTHGmBhhicQYY0xALJEYY4wJiCUSY4wxAbFEYowxJiCWSIwxxgTEEokxxpiAWCIxxhgTEAnhelIRSUTqgV2u4/DTMOCA6yBCKJbfn7236BXL7y+Q91aoqjnH+8GASyTRRETWqmqJ6zhCJZbfn7236BXL7y9U782atowxxgTEEokxxpiAWCKJbA+7DiDEYvn92XuLXrH8/kLy3qyPxBhjTEDsisQYY0xALJEYY4wJiCWSCCUit4lIhYiUi8jTIhK1i4iLyKMiUici5b22ZYvIMhHZ6ruP2pWnTvD+figilSKyUUT+LCKZDkPst+O9t14/+4qIqIgMcxFboE703kTkVhGp8v3/u9dVfIE6wd/lLBF5S0TeEZG1InJaMM5liSQCichI4AtAiapOA+KBy91GFZDHgPOO2fZ1YLmqFgHLfc+j1WP86/tbBkxT1RnAFuCOcAcVJI/xr+8NERkNLAKqwx1QED3GMe9NRBYAS4AZqjoV+JGDuILlMf71d3cv8G1VnQXc6XseMEskkSsBSBWRBGAQsNdxPP2mqq8CDcdsXgI87nv8OHBROGMKpuO9P1V9SVW7fE/fAkaFPbAgOMHvDuB+4HYgakfrnOC93QTco6rtvtfUhT2wIDnB+1Mgw/d4CEH6XLFEEoFUtQbvN6FqYB/gUdWX3EYVdHmqug/Ad5/rOJ5Q+hzwd9dBBIuIXAjUqOq7rmMJgYnAWSKyWkReEZFTXQcUZF8Cfigiu/F+xgTlStkSSQTy9RcsAcYCI4A0EbnKbVSmP0TkG0AX8KTrWIJBRAYB38DbLBKLEoAsYB7wVeAZERG3IQXVTcBtqjoauA34dTAOaokkMp0L7FDVelXtBJ4FPuI4pmCrFZHhAL77qG1COBERuQb4JLBUY2fC1ni8X3DeFZGdeJvs1otIvtOogmcP8Kx6rQF68BY6jBXX4P08AfgDYJ3tMawamCcig3zfhs4BNjuOKdiex/tHje/+OYexBJ2InAd8DbhQVVtdxxMsqlqmqrmqOkZVx+D94D1FVfc7Di1Y/gIsBBCRiUASsVUJeC/wUd/jhcDWYBzUEkkEUtXVwB+B9UAZ3t9T1JZtEJGngTeBSSKyR0SuA+4BFonIVryjf+5xGWMgTvD+fgYMBpb5hlo+5DTIfjrBe4sJJ3hvjwLjfENmfw9cE61Xkyd4f9cD94nIu8D3gRuCcq4o/TcyxhgTIeyKxBhjTEAskRhjjAmIJRJjjDEBsURijDEmIJZIjDHGBMQSiTEhJiItxzz/rIj8zFU8xgSbJRJjopSIxLuOwRiwRGKMUyJSKCLLfeuWLBeRAt/2x0Tk0l6va/HdzxeRUhF5Cu9kVWOcS3AdgDEDQKqIvNPreTbeEjHgnQH/W1V9XEQ+BzzAyUvqn4Z3rZMdwQ7UmP6wRGJM6B3xLSQEePtIgBLf09OBT/ke/w7/FhpaY0nERBJr2jImshytWdSF7/+nr3BnUq/XHA53UMZ8GEskxrj1Bv+3jPJSYJXv8U5gju/xEiAxvGEZ4z9LJMa49QXgWhHZCFwNfNG3/RHgoyKyBpiLXYWYCGbVf40xxgTErkiMMcYExBKJMcaYgFgiMcYYExBLJMYYYwJiicQYY0xALJEYY4wJiCUSY4wxAfn/B5O0Bz4Yoe0AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(xax,hourly_rearirradiance_comparison)\n",
    "plt.ylabel('$\\Delta$ in G$_{rear}$ [%] \\n(G$_{rear-with}$ - G$_{rear-without}$ / G$_{rear-without}$)')\n",
    "plt.xlabel('Hour')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}