{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# The RWTH-1 Borehole \n", "In this notebook, we'll take a look at temperature data from the RWTH-1 Borehole, drilled near the main building of the RWTH. Its initial purpose was to provide geothermal energy for space heating the newly built [Super-C](https://de.wikipedia.org/wiki/SuperC). Unfotunately, originally predicted **production** temperatures of 60 °C were not met. The water, pumped from a depth of about 2 km, has a temperature of just 31 °C when reaching the surface. \n", "But the borehole is not _cold_. Bottom Hole Temperatures (BHTs) of about 78.6 °C at a depth of about 2500 m suggest a normal average geothermal gradient at this location. \n", "\n", "In this notebook, we will visualize the average temperature gradient in the borehole. For this, we need at least information about temperatures at the surface and at depth: \n", "* T$_{surface}$ = 10.5 °C $\\pm$ 2 °C \n", "* T$_{2550 m}$ = 78.6 °C $\\pm$ 3 °C" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'depth [m]')" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVsAAAHoCAYAAADnkUw2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3XdUVWfWP/DvuY3ee0eki9hFbAjSRBC7MZZoNCoRzGTNZCaTyUrWvG8myaS8TgBRE01MLLFiV+wlllgjinSV3lV6597fH/44kycHxcK9XGB/1mK977PPcy+bO2R7OOc8++EUCoUChBBClErU3QkQQkhfQMWWEEJUgIotIYSoABVbQghRASq2hBCiAlRsCSFEBajYEkKIClCxJYQQFegRxfbgwYMICwtDcHAwtm7d2t3pEELIC5N0dwKdKS0txerVq5GYmAiZTIbXXnsNPj4+cHZ27u7UCCHkual9sb106RJGjRoFQ0NDAEBISAiSkpIQHR39zNfJ5XLU1dVBKpWC4zhVpEoI6QMUCgVaWlqgo6MDkej5Lw6ofbEtKyuDmZkZPzY3N8ft27c7fV1dXR0yMzOVmRohpA9zdXWFnp7ec89X+2Irl8uZM1OFQvFcZ6pSqRTAkw9EJpM91/fKysrCkSNH4O7uDhsbGz7+l7/8BRkZGfx4wIAB+Oyzz/g88vPzkZmZidDQULi5uT3X93pZKSkp8PLyUur36GnoM+kYfS5CXfGZNDc3IzMzk68xz0vti62lpSWuX7/Oj8vLy2Fubt7p69oLoUwmg4aGxnN9rwEDBuDixYtIS0uDpaUlxGIxAGDOnDmIiori5xUXFyMsLAzjx48HANjY2CArKwunTp3CgAED+Ncpy/P+PH0JfSYdo89FqKs+kxe9PKn2TyOMHj0aly9fxqNHj9DQ0IDjx4/zRa6rcRyH4OBgNDY2Iicnh4+PGDECo0ePZubGx8ejtbUVACASieDh4YHq6mpcvnxZKbkRQno2tS+2FhYWePfdd7Fw4UJMnToV4eHh8Pb2Vtr369evHxwdHZGdnY2WlhY+HhMTw/xLdv/+fRw6dIgfW1pawsjICBcvXkRzc7PS8iOE9ExqfxkBACIiIhAREaGy7zdp0iSkpqZCIvnvx+Pi4oKwsDAcPnyYj61fvx6hoaHQ1NQEx3Hw8vKCXC7H48ePYWFhobJ8CSHqr0cUW1UzNzeHvr4+UlNT8fuNLKKionDixAn+zLW8vBw///wzFi9eDAAwMDAAABQVFcHY2PiFL6CT7tfS0oKCggI0NjZ2dyqvTCKRIC0trbvTUCsv+ploamrC1ta2S/5bpmL7FJqamsjLy0N1dTUGDBgA4Mmlgjlz5mDz5s38vE2bNmHatGn8c8AKhQJ37txBVlYWZs+e3S25k5dXUFAAPT09ODo69vjns+vq6qCjo9PdaaiVF/lMFAoFHj58iIKCAvTr1++Vv7faX7PtThKJBDk5OaitreVjixcvhr6+Pj+uq6vDxo0b+THHcVAoFEhPT0dRUZFK8yWvrrGxESYmJj2+0JJXx3EcTExMuuyvHCq2zzBhwgSIxWLmGVt9fX3+skG7Xbt2obCwkB+7urpCLBbjyJEjKsuVdB0qtKRdV/4uULF9Bl1dXYwaNQrFxcWorKzk47Nnz4alpSU/bm1txdq1a/mxTCaDs7MzCgsLmUJNCOm7qNh2YuzYsdDU1GSKpoaGBlasWMHMS0pKYi689+vXD5qamjh58iRot3hCCBXbTmhoaCAyMpK/SdZu0qRJcHFxYWJxcXF8YRWLxRgyZAgGDRqE6upqleVLus6ZM2fg4eEBjuOU8uXh4YEzZ850949JVISeRngO7u7uMDExQUFBAeRyOYAnxTQmJgarVq3i5129ehW//vorfH19AQAmJiYAgLy8PAwYMOCFOgSR7rd8+XJkZWUp7f3T09OxfPnyZzZMunLlCtatWwepVIqCggIEBARAW1sbJ0+eBAB8++23SE1NRWxsLFpbW2Fra4v//d//hZGREU6cOIFt27ahsbERzc3N+PTTTzF06FAsWLAAAwcOxI0bN/Do0SN8+OGH8PPzU9rPSZ6g//qfk66uLm7evMk8YeDr64sRI0Yw8+Li4viCDABtbW24cOEC/x8HIS8qOTkZ//znP7Fnzx5s3boVxsbGSExMhJubG7Zv346vv/4aGzduxL59+zB27Fh89dVXkMvl2L17N9atW4cDBw5g6dKl+Pbbb/n3bGlpwY4dO/D3v/8d33zzTTf+dH0HFdvnpKGhgYaGBqSnp/PFlOM45swWADIzM3H06FF+LBaLoVAocP36deYRMqL+1q9fDw8PD6W9v4eHB9avX9/pPFdXV1hZWUFLSwtGRkb8X07W1tY4ffo0iouLsXDhQkRGRmLr1q3Izc2FSCTC119/jQsXLuCbb77B3r17UVdXx7/nuHHjADxZGfn7m79EeajYPieRSISgoCDU19cjLy+Pj3t4eCA4OJiZu3btWjQ1NTFzWlpacOzYMZXlS16dv78/v4pQGV+pqanw9/fvNI8/rl76fVc5uVyOoUOHYv/+/di/fz92796N2NhY1NXVYeHChSgoKMCIESOwYMEC5j3aO1/RY26qQ8X2Bbi6uvLtFNs7fgHA22+/zfRRKCkpwa5du/ixvr4+7OzscPfuXZSVlak0Z9K7eXt749atW3jw4AEAICEhAV988QXftW7FihXw8fHBiRMn0NbW1o2ZEiq2L4DjOISEhKCpqYn/5QYAW1tbzJo1i5n7/fffM08huLm5QSQS4cSJEyrLl/R+ZmZm+PTTT/GnP/0JERERuHv3Lv72t7/B3d0dbm5umDRpEiZPngwjIyNa0djNOEUvfQi0qamJ78re1Q2U23vW/n4HiMrKSkRGRjLXxRYuXMhc0y0vL4eRkRG8vb2hpaX1Ut/7xo0bGDZs2Etm3jt15WeSlpam1Ou0qkS9EYRe5jP54+/Ey9YWOrN9Cb6+voLdIgwNDfHGG28wse3bt6OkpIQfm5mZ8f0Weum/cYSQp6Bi+5I0NTVx+fJl1NfX87HXX3+d2ZyyubkZ69atY15XX1+Po0eP4tatWyrLlRDS/ajYviRdXV08fvyYeSBdU1MTy5cvZ+YdPnyYeTBeU1MTCoUCp0+fZm6yEUJ6Nyq2L8nAwAAjRoxAQUEBcyMsPDyc6X2pUCgQFxfHj0UiEdzd3VFbW4sLFy6oNGdCSPehYvsK/Pz8IJPJmCY1EokE0dHRzLxLly4xOwRbWFjA2NgYly9fRkNDg8ryJYR0Hyq2r0BLSwtjxoxBaWkpHj16xMfHjx+PwYMHM3NjY2OZlWeenp5obm7G+fPnVZozIaR7UCOaV+Tr64vm5mZ+Wxzgv8t433zzTT6WmpqKkydP8qvNDA0NMWrUKJiamqK1tZVZFEHUS3JyslKur0skEgwaNKjL35eoJzqzfUVSqRQBAQGCJZXe3t6CpZhr1qxhtkc3NTUFx3HIz89XSa7k5SjrRqa63SA9c+YMfvjhh2fOOXXqVKeNa3bu3IlDhw491/csLS3FW2+99cw5t2/fxpdffvlc7/eyAgICUFBQoNTvQcW2C4hEIrS2tuLatWvM87PR0dHMOvbCwkLs2bOHee3Dhw+xY8cOKrik26WkpHTaLGnixIl45513njnn5s2b/A7UnbGwsMB33333zDnZ2dl4+PDhc72fOqNi20U0NTVRWlrK/Ovo4OCAadOmMfM2bNjA/ELr6upCLpczncIIAZ70sv19A5n3338fiYmJKCgowNSpU/Hee+8hPDwcb7zxBt+56+DBgwgLC8PkyZPx8ccfo6WlBXV1dfjb3/6G6dOnIzIykj/rTExMxIIFCxAREYHo6Ghs374d27dvx549e1BaWoolS5Zg9uzZmDBhAn82m5iYiPfffx/Ak7PB//znP5g5cyYmT56MlJQUXLp0CadPn0ZsbCxOnjwJHx8f/ve9oKAAYWFhzM/Y3qO3/ef75JNPMHfuXAQEBGDPnj2orq5GbGwsTp8+jbVr16KtrQ2fffYZpk2bhilTpmDTpk38ZzVz5kxMnz4df/7znzF27FhUVFQAeLK6c+zYsWhpacH27dsxa9YshIeHY9q0abh//76S/tcTomLbRTw9PWFhYYHMzEym4cdbb73FLM2trKxktkKXyWRwcXFBcXEx7t69q9KcSc+Vnp6OxYsX49ChQ9DX18fBgwdRWlqKzz77DN9//z0OHz6MtrY2nDt3DmvXrsWAAQOQmJiIrVu3Yt26dfxfUqWlpdi7dy/i4+Px2muv4bXXXsOMGTNw6NAhhIeHY+fOnTh48CB+/PFH5iZwO0NDQ+zevRuvvfYa1q9fj9GjRyMgIACrVq1CYGAgJkyYgKSkJADAvn37MHXq1Gf+XCUlJdi2bRvWrl2LL774Avr6+li1ahUCAgIQFRWFnTt3AgD27t2L3bt349SpU/yTPjk5Ofjxxx/x9ddfIzQ0lP++x48fR1BQEJqamnD27Fls3rwZhw4dwoQJE7B169Yu+9+kM1RsuwjHcQgKCkJDQwNyc3P5uImJCebPn8/M3bp1K8rLy/mxg4MDtLW1qTMTeW4mJibw9PQE8KQnbVVVFX777TcMHTqU34z0k08+QWBgIC5duoTt27cjMjIS8+bNQ319Pb/QxtPTs8Obs0uWLIGVlRU2btyIf/3rX2hpaenwMcXO+uLOmDED+/fvBwAcOnQIkZGRz/y5xowZA47j4Orq2uH7Xb58GadPn0ZkZCRmzZqFkpIS/tHLfv36QU9PDwAwZcoUHD58mP++U6ZMga6uLj799FMcPnwYX3/9Nc6cOcOsAFU2KrZdqH///nBwcEBWVhZzI2z+/PkwNjbmx42NjUzXfLFYDHd3d1RVVeH27dsqzZmoL47jmHsAv/+d+n0DlPZ5EomE6U/7+PFjPHr0CHK5HF9++SXf83bnzp18kdTU1Ozwe3/++efYvHkzrK2tERUVBSMjow77eXTWF3fEiBEoKyvD8ePHYWtrCwsLi2f+zJ29X1tbG9577z3+Z9mxYwdmzpwp+Fm8vb35/55KS0sxZMgQFBcXY9GiRaipqcH48eMxbdo0lfYooWLbxUJDQ+Hj48OcLejo6AjuuO7fv59p02hlZYXRo0fzOzsQYmRkhPz8fDQ1NaGyshI3btx45vyBAwfi1q1b/F9NX331FU6dOoVRo0bh559/BgCUlZVhypQpKC4uFrxeLBbzT0hcvHgRS5YswaRJk/DgwQOUlpYy2z09i1gs5v9C4zgOU6dOxSeffILp06c/98/+tLxGjRqFnTt38teiX3/99af2GYmIiMDHH3+MyZMnAwDu3LkDW1tbLFq0CAMHDsTJkydV+pckFdsuZmlpiTFjxgg2d5w2bRrs7e35sVwuR3x8PD/mOA7GxsZoa2ujvqNqRlnPQHf2vi4uLvDz88PkyZPxzjvvdNpG0sLCAv/4xz+wZMkShIeHQ1NTE9OnT0d0dDQaGxv5m2nvvfce87vYbsSIETh48CA2b96M5cuX469//SvCw8OxZcsWeHl5PfejUaNHj8a6dev4a6aTJ09GQ0MDAgMDn+v1f+Tt7Y3k5GR89dVXeO211+Do6Ihp06ZhxowZmD59Onx8fDp83ZQpU5CWloYpU6YAeHKJQqFQICwsDNOmTUO/fv2U/rjX71E/WyU5ePAgqqqqmD6YJ0+e5O/kttuwYQOz2qygoACpqalYuXIl9PX1Be9L/WyFqJ9tx9Shn61cLsfPP/+MBw8e4MMPP+zWXADqZ9srKRQK3L9/n2kmPnHiRHh5eTHzYmNjmcsGhoaGaGlp4c8KCOnJoqOjsXv3brz99tvdnUq3o2KrJAEBARCLxUyTmo524719+zbOnj3Lj3V1dWFvb4/09PQOr6sR0pMkJCRg7969zA3ivoqKrZLo6urCx8cHRUVFqKqq4uNDhw7l7wS3i4+PZ5Zuurq6QiwW00KHbtJLr6yRl9CVvwtUbJVo7Nix0NDQQHp6OhOPiYlhbqDl5ubyzyICTx5/6d+/PwoKClBYWKiyfMmTx4cePnxIBZdAoVDg4cOHT3087kVRqykl0tTUREREhGD7cicnJ0RERDAF9ttvv0VYWBi/2szJyQmWlpaoqalRac59na2tLQoKCphFJz1Vc3MzsykpefHPRFNTE7a2tl3yvanYKtmAAQNgbm6OgoIC5jnF5cuXIykpCU1NTQCeNKTZunUrli5dCuDJs4V6enpoaGhAWVmZYINJohxSqZTZaaMnu3HjBrVw/IPu/EzoMoIK6OnpITk5mbnhZW5ujtdff52Z99NPPwnWn6empuKnn35iVg8RQnoeKrYqIJPJUFtbi/T0dObs9o033oCBgQE/rq+vF7SbMzY2Rl1dHc6dO6eyfAkhXY+KrQqIRCIEBgairq6O6Vurq6uLJUuWMHMTExORl5fHj83NzWFqaoqrV68yz+wSQnoWKrYq4ubmBmtra0ELxpkzZ8La2poft7W1ISEhgXmtp6cnWlpacPz4cZXlSwjpWlRsVYTjOISEhKCpqYlpQCOTyQSra06ePImUlBR+rK+vD1tbW6Snp1MLRkJ6KCq2KmRvb4+AgADY2Ngw8eDgYLi7uzOxPy7j9fDwwIQJE5htdgghPQcVWxUbN24c39y5nUgkQkxMDBO7efMmLl68yI81NDQgk8mgUCjw+PFjleRKCOk6VGy7gY6ODq5evcp0vvfx8YGvry8zLy4uTnDZ4Nq1a9i+fTutcCKkh6Fi2w20tLRQUVGBzMxMJh4TE8N0qL937x6/tUc7c3NzlJWVITk5WSW5EkK6BhXbbmBoaIjhw4cjPz+fWY7r6uqKSZMmMXPXr1+PxsZGfmxvbw8dHR2cPn2abpYR0oNQse0mfn5+kEqlgiY1UVFRkEql/Li0tBTbt2/nxyKRCB4eHqipqWGu6RJC1BsV226ira2N0aNHo7S0lFmia2VlhTlz5jBzN23axOw0amFhAWNjYyQnJ9PZLSE9BBXbbjR69Gj4+PjA0NCQiS9evJjfkhkAamtr8f333/NjjuMwZMgQ+Pr69oruVIT0BVRsu5FMJkNwcLCg5ZuBgQEWLVrExHbt2sVsBKmlpQWO4wTXfQkh6omKbTcTiUSQy+W4ceMG8zjXnDlzYGFhwY9bWlqwdu1a5rVtbW04e/YsDh48qLJ8CSEvh4qtGpBIJCguLmZ2ZdDU1MSKFSuYeUePHmVuqInFYlhYWCArK0ulWzITQl4cFVs14OXlBXNzc2RkZDA3vMLCwuDs7MzMjYuLY8YuLi6QSCS0Xxkhao6KrRrgOA5BQUFoaGhAbm4uHxeLxYJlvFeuXMGvv/7Kj2UyGVxcXFBUVITU1FSV5UwIeTFUbNWEs7Mz7O3tkZ2dzezKMHr0aAwfPpyZGxcXxzQhd3R0hJaWluC6LyFEfVCxVSOhoaEYMmQIJJL/bg3HcRxWrVrFzMvIyMCxY8f4sVgsxpgxYzBgwABm23RCiPqgYqtGrKys4Ofnx2xzDjxpHh4UFMTEEhIS0NzczI81NTWhUCiQmZnJxAkh6oGKrZrR0tJCSUkJMjIymPjKlSuZM97i4mLs3LmTmdPQ0IDjx4/jzJkzKsmVEPL8qNiqodbWVmRnZzN7jtna2mLGjBnMvB9++IFZ0KClpQUTExNcv36dFjoQomao2Kohf39/iEQiQQvGpUuXQkdHhx9XVVVh06ZNzBwPDw+0trYy13QJId2Piq0a0tPTg4+PDwoLC5kbXkZGRliwYAEzd/v27SgpKWFea29vj9TUVCZOCOleVGzV1Lhx46ChoSFowThv3jyYmpry46amJqxfv56Z4+rqCrFYjN9++00luRJCOkfFVk1pampi8uTJcHV1ZeJaWlpYtmwZEzt8+DCys7OZ1wYEBMDCwoLZeocQ0n2o2KqxgQMHwtvbW/Ao2JQpU+Do6MiP5XI54uPjmTntm0PeuXOHFjoQogao2Ko5AwMD3L17l7n+KpFIEB0dzcy7cOECbty4wcQqKipw9OhR3Lx5UyW5EkKejoqtmpNKpXj06BHS09OZJbp+fn4YNGgQMzc2NpY5izUxMYGenh7OnDmD1tZWleVMCBGiYqvmRCIRgoKCUFtbi/z8fD7e0TLeu3fv4tSpU8wcDw8P1NXV4dy5cyrLmRAiRMW2B3B3d4eVlRWysrKYFoyDBg2Cv78/Mzc+Pp5pZGNmZgZTU1NcvXoV9fX1KsuZEMKiYtsDcByHkJAQNDY24sGDB8yxlStXQiwW8+OCggIkJiYyr/Xw8EBbWxtSUlJUljMhhEXFtodwcHDA+PHjYWlpycQdHR0RGRnJxDZs2IDa2lp+bGBggMDAQEilUuaslxCiOlRsexB/f3/o6uoK4suWLYOmpiY/fvz4MbZs2cLMkUqlkMvluH37ttLzJIQIUbHtYerr63Hjxg1msYKpqSnmz5/PzNuyZQsqKiqYWE5ODg4dOoScnBxVpEoI+R0qtj1QaWkpsrKymNiCBQtgZGTEjxsbG/Htt98yc2xsbCCVSpGUlEQLHQhRMSq2PYy2tjaGDRuGvLw85rqsjo4Oli5dyszdv38/cxYrlUrh6uqK0tJS3LlzR1UpE0JAxbZH8vPzg1QqFTSpmT59Ouzs7PhxW1sb1qxZw8xxcHCAtrY2Tp06xTxGRghRLiq2PZCOjg58fX1RUlKCx48f83GpVIqVK1cyc8+cOYPk5GR+LBKJ4O7ujubmZmYnX0KIcnVLsV2wYAEmT56MyMhIREZGIjk5GQcPHkRYWBiCg4OxdetWfu6lS5cQERGB4OBgrF69ujvSVUtjxozB0KFDoaenx8QnTpyIAQMGMLE/LuO1srKCv78/ampqmCXAhBDlUXmxVSgUyMnJwf79+/kvS0tLrF69Gtu2bcO+ffuwY8cOZGdno7GxER988AESEhJw5MgRpKSk0LLT/08mk2Hy5MnQ0NBg4h0t401OTmY+N47jIBaL0draSgsdCFERlRfb+/fvAwDefPNNTJkyBVu2bMGlS5cwatQoGBoaQltbGyEhIUhKSsLt27fh4OAAOzs7SCQSREREICkpSdUpqy2RSASO43Dr1i3mzHXYsGEYO3YsMzc+Pl7QjObOnTs4cOAAKisrVZIvIX2ZyottdXU1fH19sWbNGmzatAnbt29HUVERzMzM+Dnm5uYoLS1FWVlZh3HyXxzHoaCgAEVFRUw8Ojqa6YObk5ODAwcOMHP69esHuVxO/4ARogKSzqd0rSFDhmDIkCH8eObMmfjss88QFRXFxxQKBTiOg1wuB8dxgviL6I1/Jv++b61CoYC+vj4yMjJgZWXFF1hnZ2dMnjwZBw8e5OeuX78ekyZNgpaWFgBAV1cXDg4OyMjIwNmzZwXXf3uSP/byJU/Q5yLUXZ+Jyovt9evX0dLSAl9fXwBPioWNjQ3Ky8v5OeXl5TA3N4elpWWH8Rfh5eUluK7Zk924cQPDhg1jYgYGBti2bRtyc3PRr18/Pr5ixQocP34cTU1NAICHDx9i69atzPO4rq6uKCgoQHZ2tuA53Z6io8+E0OfSka74TJqaml7qJE7llxFqamrwxRdfoKmpCbW1tdi7dy++/PJLXL58GY8ePUJDQwOOHz+O8ePHY9CgQXjw4AFyc3PR1taGQ4cOYfz48apOWe05OzvDzs4OWVlZzHVZCwsLvPbaa8zczZs3M4+LyWQyODs7o7a2lq7dEqJEKj+z9ff3R3JyMqZOnQq5XI7XX38dw4YNw7vvvouFCxeipaUFM2fOhLe3NwDg888/R0xMDJqamuDn54fQ0FBVp6z2OI5DaGgofvvtN8F+ZYsWLcK+ffv4LdHr6uqwYcMGvPfee/wcJycn9O/fH8XFxTAwMHjhSzWEkM5xil66SL79VL8vXEZo19DQgLS0NEHfgy1btuA///kPP5ZIJNi9ezdsbW2Zea2trTA0NISbm1vXJ65E9Odyx+hzEerKywgvWltoBVkvoqWlhfLyckGTmtmzZ8Pa2poft7a2CpbxAsC1a9ewf/9+6nlLiBJQse1lGhoakJWVxWyBI5PJmKc9AODEiRO4e/cuE3NxcUFDQwPOnDmjklwJ6Uuo2PYyAQEB4DgOmZmZTDwkJERweeCPy3hNTU1hZmaGa9euoaamRiX5EtJXULHtZfT19TFy5EgUFBSgurqaj4tEIsTExDBzb9y4gYsXLzIxDw8PtLa24vjx4yrJl5C+goptLzRu3DjIZDJBC8ZRo0bBx8eHicXFxTGtFvX19WFnZ4dHjx6hublZJfkS0hdQse2FtLS0MGnSJDg5OQmO/fHs9t69ezhy5AgTGzhwIIYOHYrCwkKl5klIX0LFtpcaPHgwhgwZInju1t3dHZMmTWJi69atQ2NjIz9uf01eXh7y8/OVnywhfQAV217M0NAQ6enpguY9UVFRkEql/Li0tBQ7duxg5sjlcvz66684cOAA7VdGSBegYtuLSSQSlJeXIy0tjWkSbm1tjVmzZjFzN23axK8yA56c3bq4uKCiogK3bt1SWc6E9FZUbHsxkUiEwMBA1NbWoqCggDn25ptvQldXlx/X1NTghx9+YObY2dlBV1cXp0+fFvTCJYS8GCq2vZyHhwcsLS2RmZnJPHVgaGiIRYsWMXN37NiB4uJifiwSieDh4YHa2lpcuHBBVSkT0itRse3lOI5DSEgIGhsbmW3NAeC1116DhYUFP25pacHatWuZOebm5jAzM0NZWRntV0bIK6Bi2wc4Ojpi9OjRgl7AmpqaWL58ORM7evQos/qM4ziMHDkSLi4utEsGIa+Aim0fERQUBBsbG0F88uTJ6N+/Pz9WKBSIi4tj5rTvmnHnzh08evRI6bkS0htRse1DDAwMcOvWLeaZWrFYLFjocPnyZVy5coWJNTU14cqVKzh69KhKciWkt6Fi24dIpVIUFRUJmtSMGTMGQ4cOZWKxsbHMNVoNDQ04OjoiOzsbeXl5KsmXkN6Eim0fYmRkhCFDhiA/Px+1tbV8nOM4rFq1ipmbkZEhaEbj7OwMqVRKZ7eEvAQqtn3MhAkTIBaLkZGRwcS9vLwQGBjIxBISEphmNDKZDC4uLigpKemVuxYTokxUbPsYXV1djBo1CsXFxYINHleuXAmxWMyPi4qrn14yAAAgAElEQVSKsHv3bmaOo6MjDA0NUVhYSMt4CXkBVGz7oLFjx8Lb2xs6OjpM3M7ODtOnT2diGzduZC45iEQijBkzBubm5rQbLyEvgIptH6ShoYHIyEhoamoKjr311lvQ1tbmx1VVVfjxxx+ZORzHobW1FRcuXEBTU5PS8yWkN6Bi20eJRCKIxWLcvn2buRxgbGyMhQsXMnO3bduGsrIyJlZVVYWbN2/i1KlTKsmXkJ6Oim0fJpfLkZeXh6KiIiY+b948mJiY8OOmpiasX7+emWNkZAQLCwvcvHmT6RZGCOkYFds+zNvbG6ampsjIyGCeqdXS0sKyZcuYuQcPHsS9e/eYmIeHB+RyOZKSklSSLyE9GRXbPozjOAQGBqK+vl6wUCEyMhIODg78WC6XIz4+npmjq6sLe3t7ZGRkMN3CCCFCVGz7OFdXV9jY2CArK4vpWSuRSBAdHc3M/eWXX3Dz5k3B6w0MDGhVGSGdoGLbx3Ech9DQUHh4eIDjOObYhAkT4O3tzcRiY2OZG2oaGhoYO3YspFIp6uvrVZIzIT0RFVsCW1tbBAcHQyKRMPGOlvGmpKR0+ARCc3MzTp06RT1vCXkKKrYEwJObYo8ePUJ2djYTHzx4MPz8/JhYQkKCYJuc0tJS3Lx5E9evX1d6roT0RFRsCa+mpgaZmZmCywHR0dHMluh5eXnYu3cvM8fGxgZ6eno4d+4cWlpaVJIvIT0JFVvCCwgIAABBC8Z+/fohMjKSiX333Xeoq6vjxxzHwdPTE/X19Th37pzykyWkh6FiS3gGBgYYMWIECgoKUF1dzRxbtmwZs7z30aNH2LJlCzPHzMwMZmZmuHr1KtNPgRBCxZb8gZ+fH2QymaAFo5mZGebNm8fEtmzZgoqKCibm4eEBIyMjFBYWKj1XQnoSKraEoaWlheDgYNjb2wuOLViwAIaGhvy4oaEBGzZsYObo6+tj5MiRqK+vp2u3hPwOFVsiMGzYMAwbNoy5KQY8WTG2dOlSJrZ3717BFukAUFdXR01qCPkdKrakQ0ZGRsjMzBR0+5oxYwZsbW35cVtbG9asWSN4fW5uLq5cuYIHDx4oPVdCegIqtqRDYrEYJSUlSEtLY1aMSaVSvP3228zcM2fO4Pbt20ysf//+kMlkOHr0KO3oQAio2JKnEIvFCAwMRE1NDQoKCphjgYGB8PT0ZGJ/XMYrlUrh6uqK8vJyJCcnqyRnQtQZFVvyVJ6enrCwsEBmZiba2tr4uEgkEizjvXXrFs6fP8/E7O3toaOjg9OnTzOvJ6QvomJLnorjOAQFBaGhoQG5ubnMseHDh2P06NFMLD4+nlnGKxKJ4OHhAWNjY5SXl6skZ0LUFRVb8kz9+/eHj48PTE1NBcdiYmKYTmEPHjzAoUOHmDmWlpYYOHAgSktLqUkN6dOo2JJOhYaGMk8gtHNxcUF4eDgTW79+PRobGwVzKyoq8MsvvygtR0LUHRVb8lyMjIxw584dQSFdvnw5NDQ0+HF5eTm2bdsmeP29e/fwyy+/4PHjx0rPlRB1RMWWPBeRSIS8vDxkZWUxcUtLS8yZM4eJ/fjjj6isrGRi7u7ukMvlOHr0qNJzJUQdUbElz8XY2BiDBw9GXl4e0+0LABYtWgR9fX1+XFdXh40bNzJzdHR04OjoiKysLMGjZIT0BVRsyXMLCAiAWCwWNKnR19fHm2++ycR27dolKKouLi6QSqV0dkv6JCq25Lnp6urCx8cHRUVFgssEs2bNgqWlJT9ubW1FQkICM0cmk8HNzQ0GBgbUgpH0OVRsyQsZO3YsPD09md62wJONH6OiopjY8ePHkZqaysQcHR3Rv39/5Ofn0zJe0qdQsSUvRFNTEzNmzICOjo7g2KRJk+Dq6srE4uLiBEVVoVAgOztbsC06Ib0ZFVvywkQiEWQyGe7evcsUUpFIhJiYGGbutWvXcPnyZcF73Lt3D6dOnUJzc7PS8yVEHVCxJS+lubkZDx48QHFxMRMfNWoURo4cycTi4uKY3gjt+5U1NDTg9OnTKsmXkO5GxZa8FG9vb5iYmCAjI4NZhstxnODsNisrS/AEgomJCSwsLHDjxg3BfmeE9EZUbMlLEYlECAwMRF1dHfLz85ljHh4eCA0NZWLr1q1DU1MTE3N3d0dbWxuOHTum9HwJ6W5UbMlLc3Nzg7W1taAFIwBERUVBKpXy45KSEuzcuZOZo6enBxcXF2hpadG1W9LrUbElL43jOISGhsLFxUVwzMbGBjNnzmRi33//PaqqqpiYq6srrKysaFUZ6fWo2JJXYmdnh0mTJkEikQiOLVmyhHlErKamBps2bRLMa2trw7Vr15Cdna3MVAnpVlRsySvT1tZGVVUV7t+/z8QNDQ2xaNEiJrZjxw6UlJQwMYVCgQcPHiApKYkWOpBei4ot6RKPHz9Geno6GhoamPjcuXNhbm7Oj5ubm7F27VpmjkQigZubGx4+fIgbN26oJF9CVI2KLekSEydOBABkZmYycU1NTSxbtoyJHTlyRDDPzs4Oenp6OHv2LLO1DiG9BRVb0iUMDAwwfPhw5Ofno6amhjkWHh4OJycnfqxQKBAfH8/M4TgOHh4eqKurw7lz51SSMyGqRMWWdBk/Pz9IpVKkp6czcYlEgujoaCZ26dIlXLt2jYmZmZnByckJIpGI9isjvQ4VW9JltLW1ERgYCGtra8GNrnHjxmHIkCFMLDY2VrD6zNPTEwYGBoKbaIT0dFRsSZcaOXIkfHx8IBaLmTjHcXjnnXeYWFpaGk6cOCF4j8bGRhw/fhwVFRVKzZUQVaJiS7qcsbEx7t27h/Lycibu5eXF30hrl5CQgJaWFiYml8uRn59POzqQXoWKLelyIpEIBQUFSE1NFVxOWLlyJXPWW1hYiD179jBztLS00K9fP9y/fx85OTmqSJkQpaNiS7qcWCxGYGAgampqUFhYyByzt7fH9OnTmdiGDRsE2+Q4OztDJpPRQgfSa1CxJUoxYMAAmJubIyMjQ9CkZunSpdDW1ubHlZWV+PHHH5k5UqkULi4uKC0txZ07d1SSMyHKRMWWKAXHcQgKCkJDQwNyc3OZYyYmJpg/fz4T27ZtG8rKypiYg4MDnJyc0NTURGe3pMejYkuUxtnZGcOGDYORkZHg2Pz582FiYsKPm5qa8O233zJzRCIRPD09IZFIBLv5EtLTULElShUeHg57e3tBXFtbG2+99RYTO3DggKCZDQBUVVVhz549aGxsVFqehCgbFVuidCYmJkhNTRXs1DB16lSmEMvlcsEyXgBobW1FYWEhTp06pfRcCVEWKrZEJR48eICsrCwmJpFIsHLlSiZ2/vx53Lp1i4kZGRnBysoKv/32G11OID0WFVuidCYmJhg8eDByc3NRV1fHHAsICMDAgQOZ2DfffCO4Iebu7g65XI6kpCSl50uIMlCxJSrh7+8PkUgkaK3IcRxWrVrFxO7cuYMzZ84wMR0dHTg4OCAjI0Pw7C4hPQEVW6ISenp68PHxQWFhoWAfsiFDhmD8+PFMLD4+XtDX1tXVFa6uroIWjoT0BFRsicqMGzcOrq6u0NDQEByLjo6GSPTfX8e8vDzs27ePmSOTyeDq6or6+nrB5QhC1B0VW6IympqamDNnDrMJZDsnJydMmTKFiX333Xeor68XzC0rK8POnTup5y3pUZRabGtraxEeHs5vU33p0iVEREQgODgYq1ev5uelpaVh+vTpCAkJwT/+8Q/+z8eioiLMmzcPoaGhiIqKorOZXkAkEkFTUxNpaWmCm2DLly9nznofPnyILVu2CN6jqakJeXl5gubjhKgzpRXb5ORkzJ07l+/a1NjYiA8++AAJCQk4cuQIUlJS+O1P3nvvPXz00Uc4duwYFAoFdu7cCQD45z//iddffx1JSUnw8vJCQkKCstIlKlRfX4979+4JGoSbmZnh9ddfZ2KbN2/Gw4cPmZiNjQ309fVx7tw5QXtGQtSV0ortzp078fHHH/M7q96+fRsODg6ws7ODRCJBREQEkpKSUFhYiMbGRgwePBgAMH36dCQlJaGlpQXXrl1DSEgIEyc93+DBg2FkZISMjAzBpYA33ngDBgYG/LihoQEbNmxg5rTv6NDQ0CB4aoEQdSVR1hv/61//YsZlZWUwMzPjx+bm5igtLRXEzczMUFpaisePH0NXVxcSiYSJv6iUlJSX/AnUV2/Y7rt///64fv068vPz4eDgwMd1dXWxdOlSfP3113wsMTERc+fOZVabmZqawtzcHNeuXUNAQECv+EyUgT4Xoe76TJRWbP9ILpeD4zh+rFAowHHcU+Pt//f3/jh+Hl5eXh3e/e6pbty4gWHDhnV3Gq9MoVCgsLAQWVlZsLW1ZRqKz5gxA9u3b+efp21ra8OaNWvw73//m3kPDw8PVFZWQiqV9orPpKv1lt+VrtQVn0lTU9NLncSp7GkES0tLZpuU8vJymJubC+IVFRUwNzeHsbExampq+F6o7fNJ78BxHEJCQuDo6Ci4USaTyRAVFcXETp06JfgF19PTg52dHTiOo2u3RO2prNgOGjQIDx48QG5uLtra2nDo0CGMHz8eNjY20NDQ4E/t9+/fj/Hjx0MqlWL48OE4cuQIAGDfvn2CB99Jz+bg4ICIiAhIpVLBseDgYHh4eDCx2NjYDvva5ubmYvfu3UrLk5CuoLJiq6Ghgc8//xwxMTEICwuDk5MTQkNDAQBfffUVPvvsM4SGhqK+vh4LFy4EAHz88cfYuXMnwsLCcP36dfzpT39SVbpERbS1tVFbWyvYa0wkEgmW8d68eRMXLlwQvEdjYyMyMzNx7949ZaZKyCtR+jXb06dP8/+/r68vDhw4IJjj7u7e4ZmJjY0NNm/erNT8SPcrKytDWloaLCwsoKWlxcdHjBiB0aNH49KlS3wsLi4Ovr6+/I1T4MnNtry8PCQlJeHtt99+qWv7hCgbrSAj3S4wMBAABC0YASAmJoYpnvfv38ehQ4eYORKJBK6urqioqMBvv/2m3GQJeUlUbEm3MzQ0xLBhw5CXlyfYZdfFxQVhYWFMbP369YJdG+zs7KCrq4szZ84INpgkRB1QsSVqwc/PD1KpFOnp6YJjUVFRkMlk/Li8vBw///wzM0ckEmHgwIHw9vamBuNELVGxJWpBR0cHAQEBMDc3FzxxYGlpiTlz5jCxTZs2CYqqiYkJjI2NkZ+fT01qiNqhYkvUxqhRo+Dr68sscGi3aNEi6Onp8eO6ujps3LhRME+hUODOnTs4fPiwUnMl5EVRsSVqxcTEBDk5OaioqGDiBgYGWLx4MRPbtWuXYNcGjuPQ2NiIW7du4dGjR0rPl5DnRcWWqBWRSIScnBykpqYKLifMmTMHFhYW/Li1tRVr164VvIe7uzsA8AtiCFEHVGyJWhGLxZg4cSKqq6sFZ60aGhqCZbxJSUmCm2ra2tpwdHTEvXv3kJubq/ScCXkeVGyJ2hk4cCDMzMyQmZkpuNE1adIkuLi4MLHY2FjBezg7O0MqleLYsWNKzZWQ50XFlqgdjuMQFBSE+vp6wZmpWCxGTEwME7t69Sp+/fVXJiaTyTB48GC4u7sLnt0lpDtQsSVqydnZGYMGDYK+vr7gmK+vL0aMGMHEYmNjBWfBFhYW0NHRQV5eHj0KRrodFVuiljiOw9SpU+Ho6NjhsT82qcnMzMTRo0cFc1taWvDLL78w/RUI6Q5UbIlaMzMzQ3p6Opqampi4h4cHgoODmdjatWsF8yQSCRobG3Hx4kXBMUJUiYotUWttbW24d+8esrOzBcfefvttpvtXSUkJdu3axcxp36+ssbERp06dUnq+hDwNFVui1kxNTeHt7Y2cnBzU19czx2xtbTFz5kwm9v3336O6upqJGRsbw9LSEjdv3kRVVZXScyakI1RsidoLCAiASCRCRkaG4NiSJUugo6PDj6urq7Fp0ybBPHd3d8jlcpw4cUKZqRLyVFRsidrT19fHyJEjUVhYKDhrNTIy4nf2aLd9+3aUlJQwMV1dXQwdOhS2trZobm5Wes6E/BEVW9IjjBs3Dv3792eu0bZ7/fXXYWpqyo+bm5uxfv16wTwrKytIpVLk5+crNVdCOvLUbXGOHz/+XG/wxzvChCiDlpYW5s6dizt37giag2tpaWHFihX45JNP+NihQ4cwb948ODs7M3Pr6upw/vx5hIWFwc3NTSW5EwI8o9j+7W9/g5eX1zNffPfuXSq2RGXEYjF0dHRw/fp1uLi4MNvlhIeHY+vWrXjw4AGAJ60W4+Li8M033zDvoampiebmZhw/fhwuLi4QieiPO6IaTy22Dg4OnW62OHXq1C5PiJBnqa6uRmZmJvT19WFpacnHJRIJoqOj8ec//5mPXbx4EdevX8fw4cP5mFgshpubG5KTk3H9+nWMHDlSpfmTvuup/6wnJiZ2+uLnmUNIVxoyZAiMjIyQnp4uWII7fvx4DB48mIl1tIzX1tYW+vr6OHfuHFpaWpSeMyHAM4pt+59XR44cQUhICHx9feHr68t30//9HEJURSQSITAwELW1tSgoKGCOdbSMNzU1FSdPnhTM8/DwQH19PS5evKj0nAkBnnEZod2XX36JDz/8EPb29qrIh5BOeXh4wNLSEpmZmbCxsWG20fH29oa/vz/OnDnDx9asWQN/f39IpVI+ZmZmhmHDhkFfXx9tbW0dbsVDSFfq9NTUxsYGEydOhIuLC/NFSHfhOA4hISGwtbXtsJvXypUrmeJZWFiIPXv2COZZWVmB4zgUFxcrNV9CgOcotlOnTsW///1vXL58GdeuXeO/COlOjo6OmDp1KrPF+R+P/d6GDRs67GtbXl6Obdu2oby8XGm5EgI8x2WEK1eu4Pz587hw4QITP3jwoNKSIuR56OjooL6+HmVlZXBwcGCOvfXWWzhy5AgaGhoAAJWVldi8ebNgWx1tbW00NTXh6NGjgpVohHSlTottamoqzp8/Dw0NDVXkQ8gLKS4uRnp6OiwsLKCpqcnHTU1NMX/+fHz33Xd8bOvWrZg1axaz2kxLSwtOTk7Izs7G/fv34eTkpNL8Sd/R6WUEU1NTtLa2qiIXQl5YYGAgFAoFMjMzBcfmz58PY2NjftzY2NjhMt7+/ftDJpMhKSlJsKMvIV2l02JrYWGByMhIfPTRR/jkk0/4L0LUgZGRERwcHJCfny+4Jqujo4O33nqLiR04cAA5OTlMTCqVwtXVFeXl5YKdegnpKp0WW3t7e0ydOhXm5uYwNDTkvwhRFy4uLhCLxR22YJw2bRrz2GJbWxvi4+MF8+zt7eHj44O2tjY6uyVK8dRrtgkJCXj77bcRHR391Be3zyGkO2loaGDChAkoLi6GQqFgeiZIJBK8/fbbeP/99/nY2bNncevWLWa1mUgkgpmZGVpaWlBRUQEzMzOV/gyk93tqsd2zZw9cXFye+a98YmIiFVuiFsaMGYOSkhIUFxcLnr2dOHEivLy8kJKSwsdiY2OxceNGpjADQG5uLk6ePIno6GhoaWmpJHfSNzy12FpbW+Onn3565outrKy6PCFCXpapqSmuXr0KTU1NmJiY8PH2ZbzLli3jY7dv38bZs2fh7+/PvIeuri7q6+tx8uRJREREqCx30vs9tdh21vGLEHWUnZ0NkUiEsWPHMmetQ4cOxbhx4/DLL7/wsfj4eIwbN45pSG5kZARra2skJydjzJgxzNMMhLwK6iRDeg2JRIKJEyeiqqoKRUVFguPR0dFM86Tc3Fzs379fMM/d3R0KhQJJSUlKzZf0LVRsSa/i7e0NU1NTZGRkCK7d9u/fH+Hh4Uzs22+/5VeZtdPW1oajoyOys7NRUVGh9JxJ30DFlvQqHMchKCgI9fX1yMvLExxfvnw5sxry4cOH2Lp1q2Ces7Mzxo0bh8rKSqXmS/qO5yq2paWlSE1Nxd27d/kvQtSVi4sLBgwYwGxx3s7CwgJz585lYj/99BMePXrExGQyGfT19VFdXU0Fl3SJTnsjrF69Gj/88AOznpzjOJw6dUqpiRHysjiOw8yZM5GRkdFhp69FixZh7969qKqqAgDU19djw4YN+Otf/yqYm5qail9++QXR0dHU85a8kk7PbA8ePIjTp08zX1RoSU9gaWmJrKwsNDc3M3FdXV0sWbKEie3Zs6fDLc719fVRWVmJK1euKDVX0vt1WmyNjIyYs1pCeorm5mZkZmYiKytLcGzmzJmwtrbmx21tbUhISBDMs7a2hoGBAS5cuCAo2oS8iKcW2/Zrs56envjkk0+QnJxM12xJj2JmZoaBAwciNzcX9fX1zDGZTCZY/XjixAlmlRnw5JKEp6cnGhoacPr0aaXnTHqvp16zjYmJYca//0Wja7akpwgICEBKSgoyMzMFO+8GBwdjy5YtTKev2NhYrF+/nlkQYWJiAgsLCyQnJyMgIKDD3SEI6cxTi217cS0pKYGlpSVzrKM/ywhRRwYGBhgxYgSuXLkCJycn6Ovr88dEIhFiYmKwcuVKPnbz5k1cvHgRY8eOZd7Hy8sLIpEIZWVlsLW1VVn+pPd46mWEyspKVFZWYtmyZaiqqkJlZSWqqqpQUVEhOOslRJ35+fnBwcGBWT3WzsfHB6NGjWJicXFxaGtrY2JaWlrQ0NBAaWkpampqlJov6Z2eWmz//Oc/Y9SoUcjMzOR/IX18fDBhwgQMGDBAlTkS8kq0tLSwYMECGBgYdHg8JiaGuWxw7949HD58WDBPoVDgypUr2L17t9JyJb3XU4vtxo0bkZ6ejmnTpiE9PZ3/SklJwddff63KHAl5ZWKxGHp6esjOzha0DXVzc0NoaCgTW79+PRobG5kYx3EwMjJCXl5eh9vwEPIsnT769dlnnyErKwubNm3CTz/9hNzcXFXkRUiXe/ToEdLT01FWViY4FhUVBalUyo9LS0uxfft2wTwnJydoaGjg+PHjtKMDeSGdFtvt27dj4cKFSE1Nxe3btzF37lwcOXJEFbkR0qWGDh0KQ0NDpKenCwqltbU15syZw8Q2bdokWKorkUjg5uaGhw8f4saNG0rPmfQenRbbH374Afv27cMXX3yBr776Crt37+7w4W9C1J1YLEZgYCBqampQUFAgOL548WLo6enx49raWvzwww+CeXZ2dtDT08PVq1fp7JY8t06Lra6uLiwsLPixtbU1PWdIeixPT09YWFggMzNT8MSBgYEBFi1axMR27twp6I3LcRyGDx+OkSNH4uHDh8pOmfQSnRbbMWPG4OOPP0ZmZibu3buH//u//4OjoyOtJCM9EsdxCA4OhqWlpaDYAsCcOXOYk4uWlhasXbtWME9HRwccxyE3N1dwI42QjnTa9evQoUMAwGwnAvz3cRlaSUZ6GicnJ1hYWCAjI0NwGUBTUxMrVqzAP//5Tz529OhRzJs3D+7u7szc1tZWnD9/Hnl5eZg2bZpKcic9V6fFltaDk95IR0cHTU1NKC0thZ2dHXMsLCwMW7duRXZ2Nh+Li4vDmjVrmHkSiQTGxsZISUnBuHHjqGETeaZOLyPU1dXhf/7nf/DGG2+gsrISH330Eerq6lSRGyFKlZ+fj5SUFMFlALFYjOjoaCZ25coV/Prrr4L3cHNzA/Dk7JeQZ+m02H7yySfQ09PDw4cPoaGhgdraWnz00UeqyI0QpQoMDIRcLu+w18eYMWMwbNgwJhYXFyfY10xLSwtOTk64f/8+cnJylJku6eE6LbZpaWl49913IZFIoKWlha+++gppaWmqyI0QpTI2NsaQIUOQl5cn2NGB4zisWrWKiWVkZODYsWOC9+nfvz9kMhkuX76s1HxJz9Zpsf1j8462trYOG3oQ0hP5+/tDLBZ3uPx2wIABCAoKYmJr164VNBGXSqUYM2YM3NzcqEkNeapOq+aIESPw5ZdforGxEb/88gtiYmLg4+OjitwIUTpdXV2MGzcO+vr6HS5QWLlyJSSS/95HLioqwq5duwTzdHR0oFAokJWV1eEjZYR0Wmz/8pe/QFtbG3p6eli9ejXc3Nw63BiPkJ5q/PjxGD9+fId/sdna2mLGjBlM7Pvvv+/wDLa2thZHjx7FxYsXlZYr6bk6LbZSqRQrV67Erl27kJiYiHfffRcaGhqqyI0QlTE3N0dJSYlgS3MAWLp0KbMtelVVFX788UfBPB0dHejp6eHSpUu00IEIPPU52wULFjA9Pv/op59+UkpChHQHhUKBtLQ0iMVijB07lvndNzIywoIFC7Bu3To+9vPPP2PWrFnMarP2/couXryIU6dOYfLkySr9GYh6e+qZ7fz58zFv3jyYm5tDW1sbCxYswKJFi2BkZAQHBwdV5kiI0kkkEkycOBFVVVUoLi4WHJ83bx6zaKGpqYkpvu2MjIxgZWWF3377TdAxjPRtTy22ISEhCAkJQV5eHhISEhAUFISAgACsXr2a2SCPkN7C29sbJiYmyMjI6PB52mXLljGxw4cPM6vM2rm7u0OhUOD69etKzZf0LJ1es338+DGampr4cV1dHaqqqpSaFCHdQSQSISgoCHV1dcjPzxccnzJlChwdHfmxXC5HfHy8YJ6Ojg78/f1hbGzM/LdD+rZOi214eDhmz56N2NhYfPPNN5g9ezZmz56titwIUTlXV1e4u7t3eBNYIpEIlvFeuHChwybiWlpa/HVgQoDnKLbvvPMO3nnnHVRXV6Ompgbvv/8+li5dqorcCFE5juMwZ84cODs7d3jcz88P3t7eTCw2NrbDZ3TLysqwf/9+KrgEwHMUW+DJGvIPP/wQH374Ifz8/JSdEyHdzsrKCvfv3xesFuM4Du+88w4Tu3v3boetRk1MTKClpYUTJ04IrgGTvofW3RLSgaamJqSlpXV4A2zQoEGYMGECE1uzZg1aW1uZmFgshpubGx4/foyrV68qM13SA1CxJaQDZmZm8PLyQk5ODhoaGgTHo6OjIRaL+XF+fj4SExMF82xsbGBgYIDz58+jpaVFqTkT9UbFlpCnCAgIAIAOm9Q4OjoiMjKSiX333XeCXhatC6kAACAASURBVM8cx8HDwwNNTU20jVQfR8WWkKcwNDTE8OHDkZ+f32EvhGXLlkFTU5MfP378GFu2bBHMMzU1xcSJEwGAmtT0YVRsCXkGPz8/2NjYdPi0gampKebPn8/ENm/ejIqKCsFcDQ0NtLW10YKgPoyKLSHPoK2tjcWLF8PQ0LDD4wsWLICRkRE/bmxsxLffftvh3Pv372PPnj0oKSlRSq5EvVGxJaQTYrEYhoaGePDggeAMV0dHR/Dc+f79+zvcIsfa2hoikQhJSUnKTJeoKaUX29raWoSHh6OgoAAA8Pe//x3BwcGIjIxEZGQkTpw4AQC4dOkSIiIiEBwcjNWrV/OvT0tLw/Tp0xESEoJ//OMfgsdrCFGFiooK3L17F2VlZYJj06dPh62tLT9ua2sT7MQLPNkm3cnJCbm5uR0+UkZ6N6UW2+TkZMydO5f5Vz4lJQVbtmzB/v37sX//fgQFBaGxsREffPABEhIScOTIEaSkpODcuXMAgPfeew8fffQRjh07BoVCgZ07dyozZUI6NHToUBgYGCA9PV1wdtve8/n3zpw5g+TkZMH79O/fHxoaGvzvM+k7lFpsd+7ciY8//hjm5uYAgIaGBhQVFeGDDz5AREQEYmNjIZfLcfv2bTg4OMDOzg4SiQQRERFISkpCYWEhGhsbMXjwYABPziDoTzDSHcRiMQIDA1FTU4PCwkLB8cDAQHh6ejKxjpbxSiQSuLq6orq6usP3Ib3XU5uHd4V//etfzLiiogKjRo3Cxx9/DD09PSxfvhy7d++GtrY2zMzM+Hnm5uYoLS1FWVkZEzczM0NpaekL5ZCSkvJqP4Qa6qjxSV+nis9EoVDAwMAAGRkZsLKyYhY1tO/Gu2LFCj6WnJyMc+fOCVab2dnZwcbGBqWlpS/8+/yi6HdFqLs+E6UW2z+ys7NjrmUtWLAA+/btQ0hICNMZX6FQgOM4yOXyDuMvwsvLq1dt43Pjxg0MGzasu9NQK6r8TIyMjHDu3Dm0tbUxxRYAhg8fjrFjx+LChQt8LD4+HmPHjmU2jRSJRPx+Z7q6unBzc1NKrvS7ItQVn0lTU9NLncSp9GmEjIwMHDt2jB8rFApIJBJYWlqivLycj5eXl8Pc3FwQr6io4C9JENId+vfvj9mzZz/1H/Do6Ghm48icnBwcPHiww7nJyclITExEfX29UnIl6kWlxVahUODTTz9FVVUVWlpasGPHDgQFBWHQoEF48OABcnNz0dbWhkOHDmH8+PGwsbGBhoYGf9q/f/9+jB8/XpUpEyKgq6uL1tbWDq+5Ojs7C/YeW79+fYf9FRwcHNDc3Mw/kUN6N5UWW3d3dyxbtgxz587F5MmT4eHhgfDwcGhoaODzzz9HTEwMwsLC4OTkhNDQUADAV199hc8++wyhoaGor6/HwoULVZkyIR26f/8+bt++3eFODCtWrGDOfCsqKrBt2zbBPENDQ9jY2OD27dsd7upLehdO0UufP2m/rkLXbHu/7vhMHj58iDVr1sDe3h4DBw4UHI+Li2O2O9fR0cG+ffuY1WYAUF9fj7Nnz6Jfv36Cpb+vin5XhLrymu2L1hZaQUbISzAxMcHgwYORl5cn6PQFAIsWLYKBgQE/rqurw4YNGwTztLW14ejoiPLyctTW1io1Z9K9qNgS8pL8/f0hEok6bMGop6eHxYsXM7E9e/bwKyl/z83NDePGjetwC3XSe1CxJeQl6enpYcyYMdDW1u5wNdjs2bNhZWXFj1tbWztcxisWi8FxHB49eoTc3Fyl5ky6DxVbQl6Bv78/JkyYIHjmFgBkMhmioqKY2IkTJ57aRPzXX39FYmIi9bztpajYEvKKLCwsUFZWhsePHwuOhYaGwtXVlYnFxcV1eCbs5OSE6upqXL58WWm5ku5DxZaQV6RQKHDnzh3cvXtXUERFIhFWrVrFxK5fv45Lly4J3sfKygqGhoa4ePEiGhsblZozUT0qtoS8IolEgoCAAFRWVnbYGHzUqFHw8fFhYnFxcYLLBRzHwdPTE42NjTh9+rRScyaqR8WWkC4wePBgGBsbIyMjA3K5XHA8JiaGGWdnZ+Po0aOCecbGxrC0tERRURFdu+1lqNgS0gVEIhGCgoJQW1uL/Px8wXF3d3dMmjSJia1du7bDywVDhgzB0KFDUVRUpLR8iepRsSWki7i5ucHZ2RlSqbTD41FRUcyx0v/X3p1GRXGueQD/Vy8sTbNDiyKoyI7YQTSJ1ygiICIao17HuMSLRiVGMZmTo0lMTjKZcWLMOOMd3MKo1yRqxn2Jxj3GRBY31Cg7sgvIToNs3XTXfMihx7pd5dWkGxp4fl8m9dSTuq+dzn8q1fW+b1UVDh48aNAnFovBsiyKi4t5d4YgvROFLSFGwjAMFixYAB8fH97zgwYNwpw5czi1r7/+GiqVyqBXq9UiOTlZcMUw0vtQ2BJiZO7u7igpKYFarTY4t2TJEsjlcv1xc3Mz9uzZY9AnFovh5eWFhw8fIjc316TjJd2DwpYQI2tra8P9+/dRUFBgcM7BwQFxcXGc2sGDB3mn6g4bNgxWVla4cOEC749upHehsCXEyBQKBYKCglBUVMS7ju3rr7/OWQRfo9Fgx44dBn1isRh+fn6or6/HrVu3TDpmYnoUtoSYQGRkJAAgPz/f4JyVlRXi4+M5tbNnz/IuaDN48GDY29ujpKSEduPt5ShsCTEBBwcHhIaGorS0lHfpxGnTpsHLy0t/zLIstmzZYtDHMAzGjRsHb29v1NbWmnTMxLQobAkxkbCwMAwYMIB3coJYLDaYxpuWloYbN24Y9IpEIuh0Oty7d4/WvO3FKGwJMREbGxssW7bMYHeGLuPGjcOoUaM4tcTERN4fw1pbW5GcnMzZMJX0LhS2hJiQWCyGo6Mj7zNXhmEM7m5zcnJw4cIFg+vIZDJ4eHggKyuLJjr0UhS2hJhYdXU17t+/j5qaGoNzI0aM0P+Y1mX79u287+h2LdXIt6YCMX8UtoSYWGhoKOzs7JCTk8P7RsHKlSs5i49XVFTgyJEjBn3W1tYYPnw4iouLUVhYaNIxE+OjsCXExMRiMSIiItDU1ITy8nKD8x4eHpg1axantnv3bt4fw4YPHw65XE5h2wtR2BLSDYKDg+Hq6oq8vDzeH8CWLVsGmUymP1apVJyt0LtIJBKEhYXB2dkZzc3NJh0zMS4KW0K6AcMwiIqKgoODAzQajcF5JycnvPHGG5zad999x/tjGMMw0Gq1SElJQWdnp8nGTIyLwpaQbuLj44N58+bB2tqa9/yCBQvg7OysP+7o6EBSUhJvb21tLa5fv47k5GSTjJUYH4UtId1ILpdDp9PxLjwjk8mwbNkyTu3UqVO8z2ddXFzg5OSEtLQ03vUXiPmhsCWkm+Xn5+Pu3bvo6OgwOPfaa6/B09NTf6zT6QSn8QYEBECtVuPixYsmHS8xDgpbQrpZZGQkdDod7yI1EonEYL+yq1ev4s6dOwa9jo6OGDRoEO7du4f6+nqTjZcYB4UtId3MxcUFSqUSJSUlaG1tNTg/ceJEjBw5klNLTEzkfUfX398fMpkMRUVFJhsvMQ4KW0J6QHh4OEQiEe8uDHzTeO/fv8+7vblMJsOECRPAMAzvYwliPihsCekBdnZ2GDt2LCwtLXnvWF944QWEhYVxatu2beN91YthGGg0Gvz8888mGy/54yhsCekhERERiIiI4EzVfdKqVasgEv3/v6KlpaU4fvw4b29ZWRnS0tKQmZlpkrGSP47ClpAeNGDAANTV1aGxsdHg3LBhwzBjxgxObefOnWhpaTHo9fT0hEwmw8WLF3nXzyU9j8KWkB529+5dZGZm8j5OWL58OSwtLfXH9fX12Ldvn0Ff135lKpUK169fN+l4ye9DYUtIDxKLxQgPD0dDQwMePXpkcN7V1RULFizg1Pbt28e7Rc6gQYNgb2+P5ORk3iUaSc+isCWkh4WEhMDR0RG5ubm8i9QsWrQIDg4O+uO2tjbs2rXLoI9hGAQGBkImk+Hhw4cmHTN5fhS2hPQwkUiEyMhIPH78mDck5XI5li5dyqkdP34cJSUlBr3Ozs4YO3YsVCqVycZLfh8KW0LMQEBAAIYNG8Z5++BJs2fPhru7u/5Yq9Vi27ZtgtdrbW3lXTuX9BwKW0LMAMMwWLRoEfz8/HjPS6VSrFy5klO7fPky7t+/z9v/4MED3L17l3fBG9IzKGwJMSPu7u54+PAh75q3kZGRCAwM5NSEpvH6+PhAJBLRfmVmhMKWEDPS1taGu3fv4sGDBwbnRCKRwTTeO3fu4JdffjHotbS0hLe3N8rKypCXl2ey8ZJnR2FLiBlRKBQIDAxEUVER2tvbDc6PHj0af/rTnzi1rVu38k7j9fLygqWlJc6fP89790u6F4UtIWama2tzoTvShIQEMAyjPy4qKsLp06cN+sRiMfz9/WFjY0NLMJoBCltCzIyjoyNCQkJQVlbGu8Ouj48PYmNjObWkpCTeO2EPDw8olUpUVFTQ3W0Po7AlxAyFh4fD2dmZ94cyAHjrrbdgYWGhP66pqcF3330neL3a2lrcvHnT6OMkz47ClhAzZGNjg/j4eM4GkE9yc3PD3LlzObVvvvmGd0EbAMjOzsalS5d4Fysn3YPClhAzJZFI4OzsjIcPH/I+Ali8eDHs7Oz0xy0tLdi9ezfvtfz9/aHRaHD+/HmTjZc8HYUtIWbs0aNHuHv3Lu/CM3Z2dliyZAmndvjwYd4pv/b29hg8eDAyMjJ4r0VMj8KWEDM2atQo2NraIicnh/fuds6cOXBzc9Mfd3Z2YseOHbzX6pqdRhMdegaFLSFmTCKRICIiAiqVChUVFQbnLS0tsWLFCk7t/PnzyM7ONui1traGj48PpFIp2traTDZmwo/ClhAzN3LkSLi4uAguwRgTEwNfX19O7WnTeH18fFBWVmay8RJ+FLaEmDmGYRAVFQW5XM67KLhIJEJCQgKndvPmTaSlpfFej2VZ5Ofn035l3YzClpBewMfHBwsXLoRMJuM9//LLL+PFF1/k1LZs2cK7HxnLssjOzsb58+dpv7JuRGFLSC/AMAxsbW3Bsiyqqqp4z//93W1+fj7vj2EikQj+/v5obm5GSkqKycZMuChsCelF8vLycOfOHXR0dBicCwgIQHR0NKf21Vdf8fa6ubnB0dERqampvNN8ifFR2BLSi0RERECr1fIuwQgAb7/9NiQSif740aNHOHTokEFf135lHR0duHTpksnGS/4fhS0hvYirqyuCg4NRUlLCO/XW3d0dc+bM4dT27NmDpqYmg15HR0cMGzYMAHiXaCTGRWFLSC8zadIkAMJLML755puwsbHRHzc1NWHPnj28vUFBQRg0aBDvO7zEuChsCell7O3t8dJLL0EsFvO+S+vg4IC4uDhO7eDBg3j06BHv9To7O3H9+nXa/tzEKGwJ6YWioqIQFRUluBvvvHnzoFAo9MdqtRpfffUVb69Wq0VeXh5N4zUxCltCeiGGYeDm5gaVSgWVSmVw3srKCsuXL+fUfvjhB+Tn5xv0WlhYwNvbGxUVFcjKyjLZmPs7CltCerFbt24hIyOD93HCtGnT4OXlpT9mWRZbtmzhvc6wYcNgbW2Nixcv8k4JJn8chS0hvZRYLMbEiRPR0NDAO9FBIpFg1apVnFpqairvjg1isRh+fn5obGzE9evXTTbm/ozClpBebNSoUXBwcEBubi7v3e348eMREhLCqSUmJvLevbq7u8PDwwPt7e20X5kJUNgS0ouJxWJERkaiubmZ920ChmGwevVqTi07OxsXL17k7VUqlXBwcEBNTY3JxtxfUdgS0ssFBgbCw8ND8G40ODhY/25ul+3btwtuJtnW1oYLFy6gubnZ6GPtzyhsCenlGIZBXFwc/P39BXtWrlwJsVisPy4vL8fRo0d5e9VqNQoKCmi/MiOjsCWkDxCJRPDw8EBlZSXvHeuQIUMwc+ZMTm3Xrl14/PixQa+trS08PT2RlZUlOBGCPD8KW0L6iNbWVqSnp6OgoID3/LJly2Btba0/bmxsxLfffsvb6+vrC5FIRBMdjIjClpA+QqFQwN/fH0VFRbzLJjo7O+ONN97g1Pbv38/7Y5iVlRWGDx+O0tJSwRXGyPOhsCWkD4mKioJOp+OdKQYACxcuhLOzs/64o6MD//M//8Pb6+XlBU9PT94ZauT5UdgS0oc4OTkhJCQEpaWlvM9jZTIZli1bxqmdPHkShYWFBr0SiQQjR46ESCTiXaKRPB8KW0L6mPDwcDg4OPDu0AAAr732Gjw9PfXHOp0OW7duFbxeU1MTTpw4QWve/kEUtoT0MXK5HCtWrOA8LniSRCLBypUrObVffvkFd+/e5e1vaWlBUVERfvnlF6OPtT+hsCWkD5JKpdBqtYKLgk+aNAnBwcGc2n//93/zToxQKBRwdnbG9evXeXeHIM+GwpaQPqqsrAy3b9/mfduAbxrv/fv38dNPP/H2BgYGQq1W48KFCyYbb19n0rDdunUrYmNjERsbiy+//BLAb6sOTZ8+HZMnT8bmzZv1vdnZ2Zg1axaio6Px0Ucf6Z8PVVRUYMGCBZgyZQpWrFiBlpYWUw6ZkD7Dw8MDcrkcOTk5vHesISEhGD9+PKe2detW3mez9vb2cHd3x/3791FXV2eyMfdlJgvb1NRUJCcn4/jx4zhx4gQyMzNx+vRprFu3Dtu3b8eZM2eQkZGBn3/+GQCwZs0afPLJJzh//jxYltXvCPrZZ59h/vz5OHfuHEaMGIHt27ebasiE9ClisRgRERFQqVSorKzk7UlISODs9lBaWooTJ07w9vr7+8PLywv19fUmGW9fZ7KwdXV1xQcffAALCwtIpVIMHz4cxcXFGDJkCDw8PCCRSDB9+nScO3cO5eXlaG9vxwsvvAAAmDVrFs6dOweNRoObN28iOjqaUyeEPJuRI0fC2dkZubm5vMsqenl5Yfr06Zzazp07eZ/NWltbw9/fH83NzbyTJsjTSf5xy+/j4+Oj/+vi4mKcPXsWCxcuhKurq76uUChQVVWF6upqTt3V1RVVVVVoaGiAXC6HRCLh1J9HRkbGH/yTmJ/09PSeHoLZoc+E3507d+Dl5YUHDx6go6ODM123S3x8PM6dO6d/Vayurg779+83eB+3S01NDQ4cOIARI0aAYRiTjt8Ueuq7YrKw7ZKfn4/4+HisXbsWYrEYxcXF+nMsy4JhGOh0Os4/tK561/990vP+wx0xYgQsLS3/0J/BnKSnpyM0NLSnh2FW6DPh1/W5sCyL5uZmFBQU8N7dKhQKzJ8/n7Pd+d69ezF79mw4OTkZ9Dc1NaGkpAShoaEGbzSYO2N8Vzo6On7XTZxJfyBLT09HXFwc3nvvPcycORNubm6cX0ZramqgUCgM6rW1tVAoFHByckJzczO0Wi2nnxDy7BiGgZ2dHRiGQXV1NW/PX/7yF9jb2+uPW1tbsXPnTt7eoUOHQiaT4dKlS/p/N8k/ZrKwraysxMqVK7Fp0ybExsYCAJRKJYqKilBSUgKtVovTp09jwoQJcHd3h6Wlpf72/uTJk5gwYQKkUilGjx6NM2fOAABOnDiBCRMmmGrIhPRpOTk5uH37NtRqtcE5uVyOpUuXcmrHjh1DaWmpQa9IJIK/vz+ampqQlpZmsvH2NSYL2927d6OjowNffPEFZsyYgRkzZuDYsWP44osvkJCQgKlTp8LLywtTpkwBAGzatAkbNmzAlClT0NraikWLFgEAPv30Uxw6dAhTp07FrVu38O6775pqyIT0aZMmTYJWqxVcxWv27Nlwd3fXH2u1Wmzbto23d+DAgXB0dERKSgr9WPaMTPbM9uOPP8bHH3/Me+777783qPn7++PIkSMGdXd3d+zdu9fo4yOkv1EoFBgxYgQyMzP1W5c/ycLCAitWrOD8e/vjjz8iIyMDI0aM4PR2TXRobm5GQ0MDBg4c2C1/ht6MZpAR0o9EREQAAPLy8njPT548GQEBAZxaYmIi76QIR0dHeHp64tGjR7RIzTOgsCWkH7G3t8eYMWOg0+l4A1QkEiEhIYFTu337NpKTkwWvWVJSgtOnTxt9rH0NhS0h/Ux0dDRiYmI4M8ee9OKLL2Ls2LGc2pYtWwTfPFCpVLh3757gojfkNxS2hPQzDMPAzc0NTU1NgouCJyQkcN5pLywsFLx79fX1hVgs1r81RPhR2BLST928eRMZGRm8jxN8fX0xdepUTi0pKYn3zQMLCwt4e3ujvLwcubm5Jhtvb0dhS0g/JBaLERYWhvr6esGJDitWrICFhYX+uLq6GgcOHODtHTZsGKysrHD+/HneWWqEwpaQfis0NBT29vaCSzC6ublh7ty5nNqePXvQ2Nho0CsWixEcHAw/Pz/ar0wAhS0h/ZRYLEZkZCSam5tRXl7O2xMXFwdbW1v9cUtLC3bv3s3bO2DAALi4uKCsrIw3vPs7CltC+rGgoCC4u7tDo9Hwnre3t8fixYs5tcOHDwuGM8uyuHfvHi5fvmz0sfZ2FLaE9GMMw2DJkiUICgoS7Jk7dy4GDBigP+7s7MSOHTsEr6dSqXD9+nXerdT7MwpbQvo5kUiEwYMHo7q6mvcO19LSEitWrODUzp07h5ycHN7rBQQEQKPR4Pz58yYZb29FYUsIQXt7O27cuIHCwkLe8zExMfD29ubUEhMTeXvt7Ozg4eGBzMxMwTcd+iMKW0IIXF1d4efnh8LCQv2ODU8Si8UG03hv3LiBa9eu8V7P19cXDMPg7NmzJhlvb0RhSwgBAERFRUGn0wkuUvOnP/0Jo0eP5tQSExN536u1traGUqnE0KFD0dbWZpLx9jYUtoQQAICzszNeeOEFlJaWoqWlxeA8wzBYvXo1p5aXlye4Cau7uzvkcjnvAuT9EYUtIUQvPDwccrlc8G40MDAQkydP5tR27NjB++gBANRqNX755RfcuXPH6GPtbShsCSF6tra2SEhIeOpef2+//bZ+x2vgty2wDh8+zNsrkUjQ1NSEy5cv9/s1bylsCSEcUqkULi4uePToEe/5wYMH489//jOn9re//Y13mm7XfmWPHz9+6pq4/QGFLSHEQGVlJW7duoXa2lre82+++SZsbGz0x01NTfj66695ewcMGAAnJyekpaX16x/LKGwJIQZGjRoFuVwuuEiNo6OjflPWLgcOHOC9G+7ar0ytVuPSpUsmG7O5o7AlhBiQSCSYNGkSGhsbBR8nzJ8/Hy4uLvpjtVqNpKQk3l4HBwcolUq4urr222e3FLaEEF5KpRJOTk7Izc0VfJc2Pj6eUzt9+rTgVukeHh6wsLDot9vnUNgSQniJRCJERUVBIpEIvto1ffp0DB06VH/Msiy2bNkieM3m5macOHGiX757S2FLCBHk5+eHuLg4zo9hT5JIJFi1ahWnlpKSglu3bvH2W1hY4PHjx4ITIfoyCltCiCCGYWBvbw+RSCT4ZkJYWBiUSiWnlpiYyPvDmoWFBXx8fFBZWYmMjAyTjNlcUdgSQv6hnJwcpKenQ61WG5xjGAbvvPMOp5aVlSX45sGQIUMgk8lw6dIlwe3R+yIKW0LIPzRx4kRoNBoUFBTwnh85ciTCw8M5tW3btvGujysWi+Hv7w+VSoWbN2+aZLzmiMKWEPIPDRgwAEFBQSgqKhKcmLBy5UqIxWL98cOHD3Hs2DHe3oEDB0KpVMLS0rLf7FdGYUsIeSaRkZEAgPz8fN7zQ4cOxWuvvcap7dy5k3d7HIZh4OHhAQD9ZoFxCltCyDNxcHBAaGgoNBqN4N3osmXLYG1trT9ubGzE3r17Ba9ZXV2N/fv394vtzylsCSHPbMqUKZg6dSpEIv7ocHFxwcKFCzm1/fv3C77JYGVl1W9eBaOwJYQ8M5FIhIEDB6KlpQXNzc28PQsXLoSTk5P+uL29XXAar1wuh6enJ3JyclBZWWmSMZsLCltCyHNhGAbXr19HRkYG7+MEGxsbLF26lFP7/vvvUVxczHs9X19fiEQinDlzxhTDNRsUtoSQ5yISiTB+/HjU1dWhpqaGt2fWrFn6H8AAQKvVYuvWrby9lpaW8Pb2xsOHDwV39+0LKGwJIc9t9OjRsLOzE1yCUSKRYOXKlZzalStXcPfuXd7reXl5ITQ0VHANhr6AwpYQ8tzEYjEiIiLQ1NSE8vJy3p6IiAgEBQVxakLTeMViMQYOHIj29nY0NjaaZMw9jcKWEPK7BAcHw83NTfBulG8a77179/Dzzz8LXrO4uBhff/0178yz3o7ClhDyuzAMg6VLl2LEiBGCPaNGjcL48eM5ta1btwouIG5lZQWVSvXUQO6tKGwJIb+bWCyGp6cn6uvrBQN01apVnPdyi4uL8f333/P2KhQKuLi44MaNG2hpaTHJmHsKhS0h5A9pb29Hamqq4JsEw4cPx7Rp0zi1pKQkwTUWAgMDodFocOHCBaOPtSdR2BJC/hAXFxf4+vqioKBA8PltfHw8LC0t9cd1dXXYv38/b6+dnR0GDx6MjIwMqFQqk4y5J1DYEkL+sKioKOh0OsFFagYMGIB58+Zxat9++y3q6+t5+/38/DBu3Dg0NDQYfaw9hcKWEPKHubi4QKlUoqSkBK2trbw9f/nLX2Bvb68/bm1txa5du3h7ra2tYW9vj/r6esHr9TYUtoQQowgPD4dMJhP8YcvW1hZLlizh1I4ePYqysjLBa2ZmZuLbb7/tE2veUtgSQozCzs4Oq1evxoABAwR75syZg0GDBumPtVottm/fLtgvk8lQVVWFX3/91ahj7QkUtoQQo7GwsICrq6vgkooWFhZYsWIFp3bx4kVkZmby9nt6esLGxgaXL1/u9fuVUdgSQoyqoqIC165dQ11dHe/56Oho+Pn5cWpC03hFIhH8/f3R3NyMlJQUk4y3u1DYEkKMKiQkBDY2NoKL1IhEIqxevZpTS09PFwxTNzc3ODo64vr164ITJ3oDCltCiFFJpVKEh4ejoaEBjx494u156aWX8PLLL3NqwLKyVAAAE0RJREFUW7Zs4X1UwDAMlEolxo0bJ/h4ojegsCWEGF1ISAgcHR2Rm5sLnU7H25OQkACGYfTHBQUF+OGHH3h75XI5LC0tUVFRgfb2dpOM2dQobAkhRicSiRAZGQmGYQTD0c/PD1OmTOHUkpKSBPtZlsW1a9dw/Phxo4+3O1DYEkJMIiAgAEuWLIFcLhfsWbFiBaRSqf64qqoKBw4c4O1lGAZyuRx5eXl4+PCh0cdrahS2hBCTYBgGDg4OkEgkgtNyBw0ahH/6p3/i1L7++mvBBcR9fHwglUpx9uxZo4/X1ChsCSEmlZWVhZs3bwouCL548WLO3e/jx4+xZ88e3l4LCwt4e3ujoqICWVlZJhmvqVDYEkJMKiwsDBqNBg8ePOA97+DggLi4OE7t0KFDqKio4O0fOnQorK2tcfXq1V41jZfClhBiUm5ubggMDERRUZHgj1+vv/46Z5qvRqPBjh07eHvFYjHGjBmDUaNG9aolGClsCSEmFxkZCQDIy8vjPW9lZYW33nqLUzt37hxyc3N5++3s7CASiVBcXNxr9iujsCWEmJyjoyNCQkLQ3t4u+J/+U6dOhbe3t/6YZVls2bJF8JoajQaXLl3Cjz/+aPTxmgKFLSGkW8TExGD69Omc/cieJBaLsWrVKk7t2rVruHbtGm+/VCqFjY0N0tPT0dTUZPTxGhuFLSGkW4jFYgwcOBBtbW14/Pgxb8+4ceMQGhrKqW3ZskVwFlpAQAA6Oztx/vx5o4/X2ChsCSHdhmEYXLt2DRkZGbyPExiGMVikJjc3VzBMbW1t4enpiezsbMF1GMwFhS0hpNuIRCL9gjJCi8oEBQUhKiqKU9uxYwfUajVvv6+vL0QiEa5evWr08RoThS0hpFuNHj0atra2gkswAsDbb78NsVisP66oqMDhw4d5e62srDBu3Dh4eXkJbo9uDihsCSHdSiKRICIiAiqVSnDigoeHB2bPns2p/e1vf0NzczNvv52dHQDgwYMHZjvRgcKWENLtRo4cCYVC8dSdc5cuXQqZTKY/VqlU+OabbwT7m5qacOrUKaSnpxt1rMZCYUsI6XYMw2DZsmUYOXKkYI+TkxMWLVrEqf3v//4vqqqqePttbW1hZWWFK1eumOWODhS2hJAeIZFI4OHhgcbGRsFwXLBgAZydnfXHHR0dSEpK4u1lGAb+/v5oaWnBzz//bJIx/xEUtoSQHtPR0YHk5GQUFhbynre2tkZ8fDyndvr0aRQUFPD2KxQKODs748aNG099RNETKGwJIT3GxcUFPj4+KCwsREdHB2/Pq6++iqFDh+qPdTodtm7dytvLMAwCAwOhVqtx48YNUwz5d6OwJYT0qKioKGi1WsElGCUSicE03qtXrwr+EGZvb4+JEyfC3t7erBapobAlhPQoV1dXBAcHo6SkRPA//cPCwgx+TEtMTBR8zUsul4NlWcEA7wkUtoSQHjdp0iRYWFgIvkfLMAzeeecdTi0zM/OpK35VVlbiyJEjKCkpMepYfy8KW0JIj7O3t8fq1avh5uYm2KNUKjFx4kRObdu2bYJvMjg7O0MsFuPs2bNmMdGBwpYQYhasra2hUChQV1cn2LNq1SrONN6ysjIcO3aMt1cqlcLX1xdVVVXIyMgw+nifF4UtIcRsPHr0CGlpaYK78Q4dOhSvvvoqp7Zz5060tLTw9g8ZMgQymQyXLl2CVqs1+nifB4UtIcRsKJVKyGQyZGdnC/6nf3x8PKysrPTHDQ0N2LdvH2+vSCRCQEAAmpubBbfY6S4UtoQQsyGVShEeHo6GhgbBabkuLi5YsGABp7Zv3z7BJRvd3NwwadKkHn8NzKRhu3XrVsTGxiI2NhZffvklAODDDz/E5MmTMWPGDMyYMQMXL14EAKSmpmL69OmYPHkyNm/erL9GdnY2Zs2ahejoaHz00UdmOeeZEGI8ISEhcHBwQG5uruDd7RtvvAFHR0f9cVtbG3bu3MnbyzAMrK2todVqe3QJRpOFbWpqKpKTk3H8+HGcOHECmZmZuHjxIjIyMrBv3z6cPHkSJ0+eRFRUFNrb27Fu3Tps374dZ86cQUZGhn5u85o1a/DJJ5/g/PnzYFkWhw4dMtWQCSFmQCwWIzIyElqtVvC9W7lcjqVLl3JqJ06cQHFxseB18/PzceXKFTQ2NhpzuM/MZGHr6uqKDz74ABYWFpBKpRg+fDgqKipQUVGBdevWYfr06UhMTIROp8O9e/cwZMgQeHh4QCKRYPr06Th37hzKy8vR3t6OF154AQAwa9YsnDt3zlRDJoSYicDAQLz55puwtbUV7Jk1axYGDx6sP9Zqtdi2bZtgv5ubG3Q6XY9liMRUF/bx8dH/dXFxMc6ePYv9+/fjxo0b+PTTT2Fra4v4+HgcOXIEMpkMrq6u+n6FQoGqqipUV1dz6q6uroLPcYSYwysfxmau63X2JPpM+PX2z0Wr1aK5uRkODg4G56RSKVauXIkPP/xQX/vpp5/w66+/QqlUGvTL5XIMGTIEubm5uHLlylOD3BRMFrZd8vPzER8fj7Vr18LLy4vz/3neeOMNnDhxAtHR0WAYRl9nWRYMw0Cn0/HWn8eIESNgaWn5x/8gZiI9Pd1g99H+jj4Tfn3hczl06BAePHiASZMmQSqVGpyPjIzE3r17kZWVpa8lJiZi165dvFnh6+uLhw8f4sGDBwaPIZ5VR0fH77qJM+kPZOnp6YiLi8N7772HmTNnGuySybIsJBIJ3NzcUFNTo6/X1NRAoVAY1Gtra6FQKEw5ZEKIGRk/fjw0Go3gkop8u/H++uuvguvZWlhYwNvbG9XV1aiurjb6eJ/GZGFbWVmJlStXYtOmTYiNjQXwW7h+/vnnUKlU0Gg0OHjwIKKioqBUKlFUVISSkhJotVqcPn0aEyZMgLu7OywtLfX/KXTy5ElMmDDBVEMmhJiZgQMHIiAgAIWFhWhvb+ftGT16NMaNG8epbd26VfDNpWHDhmHSpEmoq6vr1mm8Jgvb3bt3o6OjA1988YX+Na87d+5g+fLlmDdvHmJjYxEQEIBp06bB0tISX3zxBRISEjB16lR4eXlhypQpAIBNmzZhw4YNmDJlClpbWw22ySCE9G2RkZFgWRb5+fmCPQkJCZzHBsXFxTh16hRvr1gshlQqRVtbG8rKyow+XiEMaw4rNJhA13MVembb99Fnwq8vfS6nTp1CSUkJXn75ZcHfbT777DNOwLq4uOD48eOwtrbm7b9z5w4aGhrw7rvv8j4PFvJ7s4VmkBFCzF5MTAxeffVViETCkRUfH88Jv9raWnz33XeC/R4eHmhtbcVPP/1k1LEKMfnbCD2l64ZdrVb38EiMT2j7kP6MPhN+felzcXFxQXl5ObRaLe/dqpubG5YtW4ajR4/qaxcvXsTixYt5Q9rFxQWDBw9GRkYGxowZw9k2/Wm6MuV5Hwr02ccIzc3NyMvL6+lhEEL6KF9f3+d6V7fPhq1Op0NLSwukUulzv5tLCCFCWJaFRqOBjY3NUx9r/L0+G7aEEGJO6AcyQgjpBhS2hBDSDShsCSGkG1DYEkJIN6CwJYSQbkBhSwgh3YDClhBCugGFLSGEdAMKWzPGtzux0C7E/c3GjRvxwQcfAKAdmAHg8uXLmDVrFmJiYrB+/XoA9F05efKk/t+fjRs3Aujh7wpLzFJKSgo7d+5ctqOjg1Wr1eyiRYvYU6dOsWFhYWxpaSmr0WjYJUuWsFeuXOnpoXa71NRU9qWXXmLff/99lmVZNjY2lr1z5w7Lsiz74Ycfsvv37+/J4XW70tJS9pVXXmErKytZtVrNzps3j71y5Uq//q60trayY8aMYevq6liNRsP++c9/ZlNSUnr0u0J3tmaKb3fi4uJi3l2I+5PGxkZs3rwZb731FgDQDsz4bWWrqVOnws3NDVKpFJs3b4a1tXW//q5otVrodDq0tbWhs7MTnZ2dkEgkPfpd6bNLLPZ2fLsTL1y4kHcX4v7kk08+wT//8z+jsrISAIyyA3NvV1JSAqlUirfeeguVlZWYOHEifHx8+vV3RS6X45133kFMTAysra0xZswYSKXSHv2u0J2tmcvPz8eSJUuwdu1aeHh4/OHdhnuzw4cPY+DAgRg7dqy+ZowdmHs7rVaLtLQ0fP755zh48CDu3buHsrKyfv255OTk4OjRo/jpp59w9epViEQipKSk9OhnQne2Ziw9PR2rV6/GunXrEBsbixs3bvDuQtxfnDlzBjU1NZgxYwZUKhVaW1vBMEy/34HZxcUFY8eOhZOTE4Df9uw6d+4cxGKxvqe/fVeSk5MxduxYODs7A/jtkcHu3bt79LtCd7Zmim93YqFdiPuLPXv24PTp0zh58iRWr16NSZMmYcOGDf1+B+bw8HAkJyejqakJWq0WV69exZQpU/r1d8Xf3x+pqalobW0Fy7K4fPkyXnzxxR79rtCdrZl6cnfiLq+//rp+F+KOjg6EhYXpdyHuzzZt2oSPP/4Yjx8/RlBQUL/bgVmpVGLp0qWYP38+NBoNxo0bh3nz5sHLy6vffldeeeUVZGVlYdasWZBKpQgODsby5csRFRXVY98VWjycEEK6AT1GIISQbkBhSwgh3YDClhBCugGFLSGEdAMKW0II6Qb06hcxuvXr1+PmzZsAgIKCAri7u8PKygoAcPDgQf1fmyuWZbF48WIkJibCzs7OJP8bEyZMgLW1NZYvX46goCD8y7/8CxiGwfr16zF8+HAAwK1bt7Bjxw7U1dWhs7MT7u7uWLNmDby9vbF27VqkpKRg6tSp+Oijj0wyRmJk3bbkDemXwsPD2Xv37vX0MJ6LRqNhfX19WZVKZbL/jfHjx7NZWVksy7Lsp59+yhYXF7O5ubnshg0bWJZl2bS0NHbixIlsZmam/u85duwY++KLL7L19fUsy7Lsf/3Xf7Hr16832RiJcdGdLel2+fn5+Pd//3f9jKe4uDjMnDkTqamp2LJlC1xdXVFSUgKZTIalS5di7969KC4uRkxMDN5//32kpqYiMTERCoUCRUVFkMlk2LBhA7y8vKBWq/Hll18iPT0dWq0WQUFB+OijjyCXyzFhwgSEhoYiJycHa9asgU6nw65du6BWq1FfX4/Zs2cjISEBH374IQBgwYIF2LVrF+bMmYOkpCQEBAQA+O2uNCkpCTKZDIsXL4anpycqKyvx3XffoaioCP/5n/+J9vZ2iEQirF69GmFhYU/9PObOnYs1a9ZAJBLp16JNTEzEqlWrEBgYqO+bOXMmrK2todPpTPRPhphUT6c96dv+/s5WrVazMTExbHZ2NsuyLKtSqdjo6Gj23r17bEpKChsYGKg/FxcXx86bN49Vq9VsbW0tGxAQwNbW1rIpKSmsv78/m56ezrIsy+7du5edM2cOy7Is+9e//pX9j//4D1an07Esy7IbN25k/+3f/o1l2d/uJr/66iuWZVlWq9Wy8+fPZ0tLS1mWZdmKigrW39+fbWxsNLizffIu9Mnj4uJi1tfXl719+zbLsixbX1/PTp48mS0vL2dZlmUrKyvZ8ePHs5WVlQafy99f8+8FBwezhYWFT/1s6c62d6E7W9KtCgoKUFZWhvfff19fU6vVyM7OxuDBg+Hp6Ql/f38AgIeHB1xcXCCVSuHs7AyZTIbGxkYAQFBQEEaNGgUAmDNnDtavX4/m5mZcuXIFra2tuHr1KgBAo9FwFhsJDQ0FAIhEIiQlJeHKlSs4efIkHjx4AJZl0d7eDhsbm2f+80ilUiiVSgDA7du3UVNTgxUrVujPi0Qi5OXlwc3N7bk+J5FIRHewfQyFLelWOp0ODg4OOHnypL5WU1MDOzs7pKenw8LCgtMvkfB/RZ+sd4WSSCSCVqvFJ598gnHjxgEAHj9+DI1Go+/tCtLHjx9j5syZiI6ORmhoKGbPno2LFy+C5Zm9zjAMp/7k9aysrCASifTj8PX1xYEDB/Tnq6qq9KtxPQ+lUolff/1V/2NZl08//RRTpkzhLDNJegd69Yt0K29vb4hEIvzwww8AfttpYdq0acjJyXmu62RkZCA/Px/Ab284jBkzBjY2NnjllVewd+9eaDQaaLVarFu3Dn/9618N/v6ioiK0tbXhnXfeQXh4ONLS0tDZ2QmtVguxWAyGYfT7Uzk5OSEjIwPAb/t61dfX844pJCQEBQUF+lWlMjMzER0djbq6uuf6swHA22+/jcTERGRlZelrhw8fxo8//gg/P7/nvh7peXRnS7qVhYUFduzYgc8//xxfffUVOjs78d5770GpVCI1NfWZr6NQKLBp0yaUl5fD1dVVv6FfQkICNm7ciNdee03/A9natWsN/v7AwEC88soriImJgVQqhb+/P7y8vFBaWgp3d3dMnjwZ8+bNw/bt27FmzRp89tln2L9/P4KDg/U/lP09FxcXJCYmYsOGDVCr1WBZFps2bXruRwgA8NJLL+Gzzz7Dv/7rv6KtrQ1qtRpDhgzBN99887vulEnPo1W/SK+TmpqKjRs3ch5F9DZdbzQIBfez2Lx5M1pbW+k9216CHiMQ0kPeffddHD169Hf9vWvXrsWRI0eMPCJiSnRnSwgh3YDubAkhpBtQ2BJCSDegsCWEkG5AYUsIId2AwpYQQrrB/wFlwFlgQYkbzgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "\n", "import seaborn as sns\n", "sns.set_style('whitegrid')\n", "sns.set_context('notebook')\n", "\n", "temp = np.array([10.5, 78.6]) \n", "depth = np.array([0, 2550])\n", "err = np.array([2, 3])\n", "\n", "fig = plt.figure(figsize=[5,8])\n", "plt.fill_betweenx(depth,temp-err,temp+err, color=(.8,.8,.8), label='uncertainty interval')\n", "plt.plot(temp,depth,'-', linewidth=4, color='black', label='mean')\n", "plt.plot(temp-err,depth,'--', color='gray')\n", "plt.plot(temp+err,depth,'--', color='gray')\n", "plt.legend()\n", "plt.ylim(depth)\n", "plt.gca().invert_yaxis()\n", "plt.xlabel('Temperature [°C]')\n", "plt.ylabel('depth [m]')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "With this data, what is the average geothermal gradient? " ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "The mean geothermal gradient over the borehole is 26.71 K/km\n" ] } ], "source": [ "grad = (temp[1]-temp[0])/(depth[1]-depth[0])*1000\n", "print(\"The mean geothermal gradient over the borehole is {:.2f} K/km\".format(grad))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This is obviously just an approximation, and the real geothermal gradient over the depth might look different, i.e. not like a straight line. \n", "Reasons for this may be: \n", "* layers of different thermal conductivities \n", "* (vertical) advective heat transport \n", "\n", "Core measurements of the RWTH-1 well yielded an average thermal conductivity of the rocks of 2.99 W m${-1}$ K$^{-1}$ (using the geometric mean). If we neglect possible advective heat transport, we can estimate the specific heat flow in this area by applying _Fourier's Law of Heat Conduction_: \n", "\n", "$$q = -\\lambda \\nabla T$$\n", "\n", "As we assess just the vertical specific heat flow, $\\nabla T$ becomes $\\frac{\\partial T}{\\partial z}$." ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "The magnitude of the average specific heat flow is 79.85 mW/m²\n", "Considering heat flow direction, the average specific heat flow is -79.85 mW/m²\n" ] } ], "source": [ "tc = 2.99 # thermal conductivity in W/(mK) \n", "q = -tc * grad \n", "print(\"The magnitude of the average specific heat flow is {:.2f} mW/m²\"\n", " .format(np.abs(q)))\n", "\n", "print(\"Considering heat flow direction, the average specific heat flow is {:.2f} mW/m²\"\n", " .format(q))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "If we compare this value to the position on Aachen on a [Heat flow map](http://www.geni.org/globalenergy/library/renewable-energy-resources/europe/Geothermal/Geothermal%20heat%20-%20Potential_files/6-1-100.gif), we see that despite **a lot** of simplifications, the calculated specific heat flow is still in the right interval." ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.8" }, "latex_envs": { "LaTeX_envs_menu_present": true, "autocomplete": true, "bibliofile": "biblio.bib", "cite_by": "apalike", "current_citInitial": 1, "eqLabelWithNumbers": true, "eqNumInitial": 1, "hotkeys": { "equation": "Ctrl-E", "itemize": "Ctrl-I" }, "labels_anchors": false, "latex_user_defs": false, "report_style_numbering": false, "user_envs_cfg": false }, "toc": { "colors": { "hover_highlight": "#DAA520", "navigate_num": "#000000", "navigate_text": "#333333", "running_highlight": "#FF0000", "selected_highlight": "#FFD700", "sidebar_border": "#EEEEEE", "wrapper_background": "#FFFFFF" }, "moveMenuLeft": true, "nav_menu": { "height": "21px", "width": "255px" }, "navigate_menu": true, "number_sections": false, "sideBar": true, "threshold": 4, "toc_cell": false, "toc_section_display": "block", "toc_window_display": false, "widenNotebook": false } }, "nbformat": 4, "nbformat_minor": 2 }