{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Modelling sun and sky irradiance"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "- This notebook explain how to use astk functions to simulate the natural lighting of a scenesky irradiance for diiferent conditions.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Conventions"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    " - Radiation fluxes are expressed in terms of irradiance of a horizontal unit placed at earth surface(total energetic flux passing through a square meter of ground). \n",
    " - The standard measurement of radiation at earth surface is the global horizontal irradiance (ghi, W.m-2 i.e. J.m-2.s-1), that captures radiations coming from the sun and from the sky in the shortwave domain (visible light)\n",
    " - The sun contribution is called direct normal irradiance (dni), whereas sky contribution is called diffuse horizontal irradiance (dhi)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Imports and setup location"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "from openalea.astk.data_access import montpellier_spring_2013\n",
    "from openalea.astk.sky_irradiance import sky_irradiance"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Specify a location on earth:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "Montpellier ={\n",
    "'longitude': 3.87,\n",
    "'latitude': 43.61,\n",
    "'altitude': 56,\n",
    "'timezone': 'Europe/Paris'}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Simulating clear sky conditions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "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>azimuth</th>\n",
       "      <th>zenith</th>\n",
       "      <th>elevation</th>\n",
       "      <th>ghi</th>\n",
       "      <th>dni</th>\n",
       "      <th>dhi</th>\n",
       "      <th>ppfd</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2000-06-21 07:00:00+02:00</th>\n",
       "      <td>65.222522</td>\n",
       "      <td>81.781574</td>\n",
       "      <td>8.218426</td>\n",
       "      <td>69.608228</td>\n",
       "      <td>133.668645</td>\n",
       "      <td>50.500666</td>\n",
       "      <td>167.643135</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 08:00:00+02:00</th>\n",
       "      <td>74.723897</td>\n",
       "      <td>71.647009</td>\n",
       "      <td>18.352991</td>\n",
       "      <td>241.741681</td>\n",
       "      <td>432.873390</td>\n",
       "      <td>105.442654</td>\n",
       "      <td>528.836249</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 09:00:00+02:00</th>\n",
       "      <td>84.230207</td>\n",
       "      <td>61.003100</td>\n",
       "      <td>28.996900</td>\n",
       "      <td>428.277428</td>\n",
       "      <td>708.490643</td>\n",
       "      <td>84.827871</td>\n",
       "      <td>897.324728</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 10:00:00+02:00</th>\n",
       "      <td>94.525066</td>\n",
       "      <td>50.169169</td>\n",
       "      <td>39.830831</td>\n",
       "      <td>603.268306</td>\n",
       "      <td>781.227518</td>\n",
       "      <td>102.874099</td>\n",
       "      <td>1233.792552</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 11:00:00+02:00</th>\n",
       "      <td>106.936113</td>\n",
       "      <td>39.517644</td>\n",
       "      <td>50.482356</td>\n",
       "      <td>751.612236</td>\n",
       "      <td>811.181337</td>\n",
       "      <td>125.843702</td>\n",
       "      <td>1515.360865</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 12:00:00+02:00</th>\n",
       "      <td>124.101947</td>\n",
       "      <td>29.712676</td>\n",
       "      <td>60.287324</td>\n",
       "      <td>862.071318</td>\n",
       "      <td>805.824816</td>\n",
       "      <td>162.194835</td>\n",
       "      <td>1723.771894</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 13:00:00+02:00</th>\n",
       "      <td>150.974668</td>\n",
       "      <td>22.305842</td>\n",
       "      <td>67.694158</td>\n",
       "      <td>926.637555</td>\n",
       "      <td>815.663749</td>\n",
       "      <td>172.009087</td>\n",
       "      <td>1845.290578</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 14:00:00+02:00</th>\n",
       "      <td>189.036867</td>\n",
       "      <td>20.362593</td>\n",
       "      <td>69.637407</td>\n",
       "      <td>940.709766</td>\n",
       "      <td>813.385068</td>\n",
       "      <td>178.153651</td>\n",
       "      <td>1871.754717</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 15:00:00+02:00</th>\n",
       "      <td>222.847350</td>\n",
       "      <td>25.217800</td>\n",
       "      <td>64.782200</td>\n",
       "      <td>903.290967</td>\n",
       "      <td>801.009537</td>\n",
       "      <td>178.621855</td>\n",
       "      <td>1801.370007</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 16:00:00+02:00</th>\n",
       "      <td>244.565120</td>\n",
       "      <td>34.007293</td>\n",
       "      <td>55.992707</td>\n",
       "      <td>817.036014</td>\n",
       "      <td>809.396716</td>\n",
       "      <td>146.073338</td>\n",
       "      <td>1638.894191</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 17:00:00+02:00</th>\n",
       "      <td>259.103202</td>\n",
       "      <td>44.316289</td>\n",
       "      <td>45.683711</td>\n",
       "      <td>688.127892</td>\n",
       "      <td>792.097153</td>\n",
       "      <td>121.387018</td>\n",
       "      <td>1395.154930</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 18:00:00+02:00</th>\n",
       "      <td>270.327385</td>\n",
       "      <td>55.107908</td>\n",
       "      <td>34.892092</td>\n",
       "      <td>526.051805</td>\n",
       "      <td>757.091487</td>\n",
       "      <td>92.970738</td>\n",
       "      <td>1086.065369</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 19:00:00+02:00</th>\n",
       "      <td>280.142111</td>\n",
       "      <td>65.894619</td>\n",
       "      <td>24.105381</td>\n",
       "      <td>343.602529</td>\n",
       "      <td>611.132818</td>\n",
       "      <td>94.005988</td>\n",
       "      <td>731.774188</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 20:00:00+02:00</th>\n",
       "      <td>289.566812</td>\n",
       "      <td>76.349471</td>\n",
       "      <td>13.650529</td>\n",
       "      <td>158.912224</td>\n",
       "      <td>318.904517</td>\n",
       "      <td>83.651016</td>\n",
       "      <td>359.130268</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 21:00:00+02:00</th>\n",
       "      <td>299.258754</td>\n",
       "      <td>86.099097</td>\n",
       "      <td>3.900903</td>\n",
       "      <td>15.621544</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>15.621544</td>\n",
       "      <td>41.965131</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                              azimuth     zenith  elevation         ghi  \\\n",
       "2000-06-21 07:00:00+02:00   65.222522  81.781574   8.218426   69.608228   \n",
       "2000-06-21 08:00:00+02:00   74.723897  71.647009  18.352991  241.741681   \n",
       "2000-06-21 09:00:00+02:00   84.230207  61.003100  28.996900  428.277428   \n",
       "2000-06-21 10:00:00+02:00   94.525066  50.169169  39.830831  603.268306   \n",
       "2000-06-21 11:00:00+02:00  106.936113  39.517644  50.482356  751.612236   \n",
       "2000-06-21 12:00:00+02:00  124.101947  29.712676  60.287324  862.071318   \n",
       "2000-06-21 13:00:00+02:00  150.974668  22.305842  67.694158  926.637555   \n",
       "2000-06-21 14:00:00+02:00  189.036867  20.362593  69.637407  940.709766   \n",
       "2000-06-21 15:00:00+02:00  222.847350  25.217800  64.782200  903.290967   \n",
       "2000-06-21 16:00:00+02:00  244.565120  34.007293  55.992707  817.036014   \n",
       "2000-06-21 17:00:00+02:00  259.103202  44.316289  45.683711  688.127892   \n",
       "2000-06-21 18:00:00+02:00  270.327385  55.107908  34.892092  526.051805   \n",
       "2000-06-21 19:00:00+02:00  280.142111  65.894619  24.105381  343.602529   \n",
       "2000-06-21 20:00:00+02:00  289.566812  76.349471  13.650529  158.912224   \n",
       "2000-06-21 21:00:00+02:00  299.258754  86.099097   3.900903   15.621544   \n",
       "\n",
       "                                  dni         dhi         ppfd  \n",
       "2000-06-21 07:00:00+02:00  133.668645   50.500666   167.643135  \n",
       "2000-06-21 08:00:00+02:00  432.873390  105.442654   528.836249  \n",
       "2000-06-21 09:00:00+02:00  708.490643   84.827871   897.324728  \n",
       "2000-06-21 10:00:00+02:00  781.227518  102.874099  1233.792552  \n",
       "2000-06-21 11:00:00+02:00  811.181337  125.843702  1515.360865  \n",
       "2000-06-21 12:00:00+02:00  805.824816  162.194835  1723.771894  \n",
       "2000-06-21 13:00:00+02:00  815.663749  172.009087  1845.290578  \n",
       "2000-06-21 14:00:00+02:00  813.385068  178.153651  1871.754717  \n",
       "2000-06-21 15:00:00+02:00  801.009537  178.621855  1801.370007  \n",
       "2000-06-21 16:00:00+02:00  809.396716  146.073338  1638.894191  \n",
       "2000-06-21 17:00:00+02:00  792.097153  121.387018  1395.154930  \n",
       "2000-06-21 18:00:00+02:00  757.091487   92.970738  1086.065369  \n",
       "2000-06-21 19:00:00+02:00  611.132818   94.005988   731.774188  \n",
       "2000-06-21 20:00:00+02:00  318.904517   83.651016   359.130268  \n",
       "2000-06-21 21:00:00+02:00    0.000000   15.621544    41.965131  "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    " sky_irr = sky_irradiance(daydate='2000-06-21', **Montpellier)\n",
    " sky_irr"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Simulating cloudy conditions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "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>azimuth</th>\n",
       "      <th>zenith</th>\n",
       "      <th>elevation</th>\n",
       "      <th>ghi</th>\n",
       "      <th>dni</th>\n",
       "      <th>dhi</th>\n",
       "      <th>ppfd</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2000-06-21 07:00:00+02:00</th>\n",
       "      <td>65.222522</td>\n",
       "      <td>81.781574</td>\n",
       "      <td>8.218426</td>\n",
       "      <td>13.921646</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>13.921646</td>\n",
       "      <td>37.808178</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 08:00:00+02:00</th>\n",
       "      <td>74.723897</td>\n",
       "      <td>71.647009</td>\n",
       "      <td>18.352991</td>\n",
       "      <td>48.348336</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>48.348336</td>\n",
       "      <td>120.629656</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 09:00:00+02:00</th>\n",
       "      <td>84.230207</td>\n",
       "      <td>61.003100</td>\n",
       "      <td>28.996900</td>\n",
       "      <td>85.655486</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>85.655486</td>\n",
       "      <td>205.795668</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 10:00:00+02:00</th>\n",
       "      <td>94.525066</td>\n",
       "      <td>50.169169</td>\n",
       "      <td>39.830831</td>\n",
       "      <td>120.653661</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>120.653661</td>\n",
       "      <td>283.847764</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 11:00:00+02:00</th>\n",
       "      <td>106.936113</td>\n",
       "      <td>39.517644</td>\n",
       "      <td>50.482356</td>\n",
       "      <td>150.322447</td>\n",
       "      <td>0.105614</td>\n",
       "      <td>150.240974</td>\n",
       "      <td>349.281690</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 12:00:00+02:00</th>\n",
       "      <td>124.101947</td>\n",
       "      <td>29.712676</td>\n",
       "      <td>60.287324</td>\n",
       "      <td>172.414264</td>\n",
       "      <td>0.618064</td>\n",
       "      <td>171.877462</td>\n",
       "      <td>397.754979</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 13:00:00+02:00</th>\n",
       "      <td>150.974668</td>\n",
       "      <td>22.305842</td>\n",
       "      <td>67.694158</td>\n",
       "      <td>185.327511</td>\n",
       "      <td>2.164173</td>\n",
       "      <td>183.325281</td>\n",
       "      <td>426.028282</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 14:00:00+02:00</th>\n",
       "      <td>189.036867</td>\n",
       "      <td>20.362593</td>\n",
       "      <td>69.637407</td>\n",
       "      <td>188.141953</td>\n",
       "      <td>2.292383</td>\n",
       "      <td>185.992823</td>\n",
       "      <td>432.186277</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 15:00:00+02:00</th>\n",
       "      <td>222.847350</td>\n",
       "      <td>25.217800</td>\n",
       "      <td>64.782200</td>\n",
       "      <td>180.658193</td>\n",
       "      <td>0.783678</td>\n",
       "      <td>179.949204</td>\n",
       "      <td>415.808807</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 16:00:00+02:00</th>\n",
       "      <td>244.565120</td>\n",
       "      <td>34.007293</td>\n",
       "      <td>55.992707</td>\n",
       "      <td>163.407203</td>\n",
       "      <td>0.422203</td>\n",
       "      <td>163.057211</td>\n",
       "      <td>378.010644</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 17:00:00+02:00</th>\n",
       "      <td>259.103202</td>\n",
       "      <td>44.316289</td>\n",
       "      <td>45.683711</td>\n",
       "      <td>137.625578</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>137.625578</td>\n",
       "      <td>321.337452</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 18:00:00+02:00</th>\n",
       "      <td>270.327385</td>\n",
       "      <td>55.107908</td>\n",
       "      <td>34.892092</td>\n",
       "      <td>105.210361</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>105.210361</td>\n",
       "      <td>249.555017</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 19:00:00+02:00</th>\n",
       "      <td>280.142111</td>\n",
       "      <td>65.894619</td>\n",
       "      <td>24.105381</td>\n",
       "      <td>68.720506</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>68.720506</td>\n",
       "      <td>167.479703</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 20:00:00+02:00</th>\n",
       "      <td>289.566812</td>\n",
       "      <td>76.349471</td>\n",
       "      <td>13.650529</td>\n",
       "      <td>31.782445</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>31.782445</td>\n",
       "      <td>81.596061</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-06-21 21:00:00+02:00</th>\n",
       "      <td>299.258754</td>\n",
       "      <td>86.099097</td>\n",
       "      <td>3.900903</td>\n",
       "      <td>3.124309</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>3.124309</td>\n",
       "      <td>9.353447</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                              azimuth     zenith  elevation         ghi  \\\n",
       "2000-06-21 07:00:00+02:00   65.222522  81.781574   8.218426   13.921646   \n",
       "2000-06-21 08:00:00+02:00   74.723897  71.647009  18.352991   48.348336   \n",
       "2000-06-21 09:00:00+02:00   84.230207  61.003100  28.996900   85.655486   \n",
       "2000-06-21 10:00:00+02:00   94.525066  50.169169  39.830831  120.653661   \n",
       "2000-06-21 11:00:00+02:00  106.936113  39.517644  50.482356  150.322447   \n",
       "2000-06-21 12:00:00+02:00  124.101947  29.712676  60.287324  172.414264   \n",
       "2000-06-21 13:00:00+02:00  150.974668  22.305842  67.694158  185.327511   \n",
       "2000-06-21 14:00:00+02:00  189.036867  20.362593  69.637407  188.141953   \n",
       "2000-06-21 15:00:00+02:00  222.847350  25.217800  64.782200  180.658193   \n",
       "2000-06-21 16:00:00+02:00  244.565120  34.007293  55.992707  163.407203   \n",
       "2000-06-21 17:00:00+02:00  259.103202  44.316289  45.683711  137.625578   \n",
       "2000-06-21 18:00:00+02:00  270.327385  55.107908  34.892092  105.210361   \n",
       "2000-06-21 19:00:00+02:00  280.142111  65.894619  24.105381   68.720506   \n",
       "2000-06-21 20:00:00+02:00  289.566812  76.349471  13.650529   31.782445   \n",
       "2000-06-21 21:00:00+02:00  299.258754  86.099097   3.900903    3.124309   \n",
       "\n",
       "                                dni         dhi        ppfd  \n",
       "2000-06-21 07:00:00+02:00  0.000000   13.921646   37.808178  \n",
       "2000-06-21 08:00:00+02:00  0.000000   48.348336  120.629656  \n",
       "2000-06-21 09:00:00+02:00  0.000000   85.655486  205.795668  \n",
       "2000-06-21 10:00:00+02:00  0.000000  120.653661  283.847764  \n",
       "2000-06-21 11:00:00+02:00  0.105614  150.240974  349.281690  \n",
       "2000-06-21 12:00:00+02:00  0.618064  171.877462  397.754979  \n",
       "2000-06-21 13:00:00+02:00  2.164173  183.325281  426.028282  \n",
       "2000-06-21 14:00:00+02:00  2.292383  185.992823  432.186277  \n",
       "2000-06-21 15:00:00+02:00  0.783678  179.949204  415.808807  \n",
       "2000-06-21 16:00:00+02:00  0.422203  163.057211  378.010644  \n",
       "2000-06-21 17:00:00+02:00  0.000000  137.625578  321.337452  \n",
       "2000-06-21 18:00:00+02:00  0.000000  105.210361  249.555017  \n",
       "2000-06-21 19:00:00+02:00  0.000000   68.720506  167.479703  \n",
       "2000-06-21 20:00:00+02:00  0.000000   31.782445   81.596061  \n",
       "2000-06-21 21:00:00+02:00  0.000000    3.124309    9.353447  "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    " sky_irr = sky_irradiance(daydate='2000-06-21', attenuation=0.2,**Montpellier)\n",
    " sky_irr"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Simulating actual irradiances"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In general, meteorological conditions are between overcast and clear sky conditions. The cursor between theses two extremes is a function of the ratio between actual irradiance (measured) and clear sky irradiance. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "get some meteorological data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "meteo_db = montpellier_spring_2013()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2013-05-26 00:00:00+02:00      0.000000\n",
       "2013-05-26 01:00:00+02:00      0.000000\n",
       "2013-05-26 02:00:00+02:00      0.000000\n",
       "2013-05-26 03:00:00+02:00      0.000000\n",
       "2013-05-26 04:00:00+02:00      0.000000\n",
       "2013-05-26 05:00:00+02:00      0.000000\n",
       "2013-05-26 06:00:00+02:00     31.615385\n",
       "2013-05-26 07:00:00+02:00    166.370370\n",
       "2013-05-26 08:00:00+02:00    335.425926\n",
       "2013-05-26 09:00:00+02:00    505.481481\n",
       "2013-05-26 10:00:00+02:00    660.689655\n",
       "2013-05-26 11:00:00+02:00    783.516667\n",
       "2013-05-26 12:00:00+02:00    865.300000\n",
       "2013-05-26 13:00:00+02:00    887.883333\n",
       "2013-05-26 14:00:00+02:00    843.233333\n",
       "2013-05-26 15:00:00+02:00    861.366667\n",
       "2013-05-26 16:00:00+02:00    700.766667\n",
       "2013-05-26 17:00:00+02:00    576.966667\n",
       "2013-05-26 18:00:00+02:00    415.766667\n",
       "2013-05-26 19:00:00+02:00    241.550000\n",
       "2013-05-26 20:00:00+02:00     68.566667\n",
       "2013-05-26 21:00:00+02:00      1.433333\n",
       "2013-05-26 22:00:00+02:00      0.000000\n",
       "2013-05-26 23:00:00+02:00      0.000000\n",
       "2013-05-27 00:00:00+02:00      0.000000\n",
       "Name: ghi, dtype: float64"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas\n",
    "when = pandas.date_range(start='2013-05-26', freq='h', periods=25,\n",
    "                             tz='Europe/Paris')\n",
    "observed = meteo_db.loc[when,'ghi']\n",
    "observed"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "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>azimuth</th>\n",
       "      <th>zenith</th>\n",
       "      <th>elevation</th>\n",
       "      <th>ghi</th>\n",
       "      <th>dni</th>\n",
       "      <th>dhi</th>\n",
       "      <th>ppfd</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2013-05-26 07:00:00+02:00</th>\n",
       "      <td>67.394678</td>\n",
       "      <td>82.902288</td>\n",
       "      <td>7.097712</td>\n",
       "      <td>166.370370</td>\n",
       "      <td>382.539320</td>\n",
       "      <td>119.103104</td>\n",
       "      <td>372.883646</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 08:00:00+02:00</th>\n",
       "      <td>76.998278</td>\n",
       "      <td>72.595513</td>\n",
       "      <td>17.404487</td>\n",
       "      <td>335.425926</td>\n",
       "      <td>708.145583</td>\n",
       "      <td>123.608595</td>\n",
       "      <td>712.692855</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 09:00:00+02:00</th>\n",
       "      <td>86.670875</td>\n",
       "      <td>61.831377</td>\n",
       "      <td>28.168623</td>\n",
       "      <td>505.481481</td>\n",
       "      <td>756.869728</td>\n",
       "      <td>148.187460</td>\n",
       "      <td>1043.076789</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 10:00:00+02:00</th>\n",
       "      <td>97.240796</td>\n",
       "      <td>50.953097</td>\n",
       "      <td>39.046903</td>\n",
       "      <td>660.689655</td>\n",
       "      <td>850.791955</td>\n",
       "      <td>124.727856</td>\n",
       "      <td>1339.890179</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 11:00:00+02:00</th>\n",
       "      <td>110.082221</td>\n",
       "      <td>40.367996</td>\n",
       "      <td>49.632004</td>\n",
       "      <td>783.516667</td>\n",
       "      <td>882.613527</td>\n",
       "      <td>111.053233</td>\n",
       "      <td>1573.703322</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 12:00:00+02:00</th>\n",
       "      <td>127.790259</td>\n",
       "      <td>30.809505</td>\n",
       "      <td>59.190495</td>\n",
       "      <td>865.300000</td>\n",
       "      <td>820.208961</td>\n",
       "      <td>160.843080</td>\n",
       "      <td>1729.864628</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 13:00:00+02:00</th>\n",
       "      <td>154.619856</td>\n",
       "      <td>23.898554</td>\n",
       "      <td>66.101446</td>\n",
       "      <td>887.883333</td>\n",
       "      <td>676.893420</td>\n",
       "      <td>269.023929</td>\n",
       "      <td>1775.511620</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 14:00:00+02:00</th>\n",
       "      <td>190.458560</td>\n",
       "      <td>22.437564</td>\n",
       "      <td>67.562436</td>\n",
       "      <td>843.233333</td>\n",
       "      <td>578.462848</td>\n",
       "      <td>308.562438</td>\n",
       "      <td>1695.371933</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 15:00:00+02:00</th>\n",
       "      <td>222.025552</td>\n",
       "      <td>27.340165</td>\n",
       "      <td>62.659835</td>\n",
       "      <td>861.366667</td>\n",
       "      <td>669.045851</td>\n",
       "      <td>267.056253</td>\n",
       "      <td>1725.259707</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 16:00:00+02:00</th>\n",
       "      <td>243.129279</td>\n",
       "      <td>36.027976</td>\n",
       "      <td>53.972024</td>\n",
       "      <td>700.766667</td>\n",
       "      <td>530.630224</td>\n",
       "      <td>271.630139</td>\n",
       "      <td>1425.687318</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 17:00:00+02:00</th>\n",
       "      <td>257.658296</td>\n",
       "      <td>46.287318</td>\n",
       "      <td>43.712682</td>\n",
       "      <td>576.966667</td>\n",
       "      <td>510.687303</td>\n",
       "      <td>224.060075</td>\n",
       "      <td>1188.519867</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 18:00:00+02:00</th>\n",
       "      <td>268.996379</td>\n",
       "      <td>57.092156</td>\n",
       "      <td>32.907844</td>\n",
       "      <td>415.766667</td>\n",
       "      <td>402.926728</td>\n",
       "      <td>196.860848</td>\n",
       "      <td>876.394984</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 19:00:00+02:00</th>\n",
       "      <td>278.938942</td>\n",
       "      <td>67.937211</td>\n",
       "      <td>22.062789</td>\n",
       "      <td>241.550000</td>\n",
       "      <td>248.774528</td>\n",
       "      <td>148.104704</td>\n",
       "      <td>530.208075</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 20:00:00+02:00</th>\n",
       "      <td>288.494733</td>\n",
       "      <td>78.476514</td>\n",
       "      <td>11.523486</td>\n",
       "      <td>68.566667</td>\n",
       "      <td>6.658009</td>\n",
       "      <td>67.236599</td>\n",
       "      <td>165.803474</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2013-05-26 21:00:00+02:00</th>\n",
       "      <td>298.338059</td>\n",
       "      <td>88.229752</td>\n",
       "      <td>1.770248</td>\n",
       "      <td>1.433333</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.433333</td>\n",
       "      <td>4.500906</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                              azimuth     zenith  elevation         ghi  \\\n",
       "2013-05-26 07:00:00+02:00   67.394678  82.902288   7.097712  166.370370   \n",
       "2013-05-26 08:00:00+02:00   76.998278  72.595513  17.404487  335.425926   \n",
       "2013-05-26 09:00:00+02:00   86.670875  61.831377  28.168623  505.481481   \n",
       "2013-05-26 10:00:00+02:00   97.240796  50.953097  39.046903  660.689655   \n",
       "2013-05-26 11:00:00+02:00  110.082221  40.367996  49.632004  783.516667   \n",
       "2013-05-26 12:00:00+02:00  127.790259  30.809505  59.190495  865.300000   \n",
       "2013-05-26 13:00:00+02:00  154.619856  23.898554  66.101446  887.883333   \n",
       "2013-05-26 14:00:00+02:00  190.458560  22.437564  67.562436  843.233333   \n",
       "2013-05-26 15:00:00+02:00  222.025552  27.340165  62.659835  861.366667   \n",
       "2013-05-26 16:00:00+02:00  243.129279  36.027976  53.972024  700.766667   \n",
       "2013-05-26 17:00:00+02:00  257.658296  46.287318  43.712682  576.966667   \n",
       "2013-05-26 18:00:00+02:00  268.996379  57.092156  32.907844  415.766667   \n",
       "2013-05-26 19:00:00+02:00  278.938942  67.937211  22.062789  241.550000   \n",
       "2013-05-26 20:00:00+02:00  288.494733  78.476514  11.523486   68.566667   \n",
       "2013-05-26 21:00:00+02:00  298.338059  88.229752   1.770248    1.433333   \n",
       "\n",
       "                                  dni         dhi         ppfd  \n",
       "2013-05-26 07:00:00+02:00  382.539320  119.103104   372.883646  \n",
       "2013-05-26 08:00:00+02:00  708.145583  123.608595   712.692855  \n",
       "2013-05-26 09:00:00+02:00  756.869728  148.187460  1043.076789  \n",
       "2013-05-26 10:00:00+02:00  850.791955  124.727856  1339.890179  \n",
       "2013-05-26 11:00:00+02:00  882.613527  111.053233  1573.703322  \n",
       "2013-05-26 12:00:00+02:00  820.208961  160.843080  1729.864628  \n",
       "2013-05-26 13:00:00+02:00  676.893420  269.023929  1775.511620  \n",
       "2013-05-26 14:00:00+02:00  578.462848  308.562438  1695.371933  \n",
       "2013-05-26 15:00:00+02:00  669.045851  267.056253  1725.259707  \n",
       "2013-05-26 16:00:00+02:00  530.630224  271.630139  1425.687318  \n",
       "2013-05-26 17:00:00+02:00  510.687303  224.060075  1188.519867  \n",
       "2013-05-26 18:00:00+02:00  402.926728  196.860848   876.394984  \n",
       "2013-05-26 19:00:00+02:00  248.774528  148.104704   530.208075  \n",
       "2013-05-26 20:00:00+02:00    6.658009   67.236599   165.803474  \n",
       "2013-05-26 21:00:00+02:00    0.000000    1.433333     4.500906  "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sky_irr = sky_irradiance(dates=observed.index, ghi=observed)\n",
    "sky_irr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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.10.14"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}