{ "nbformat": 4, "nbformat_minor": 0, "metadata": { "colab": { "provenance": [] }, "kernelspec": { "name": "python3", "display_name": "Python 3" }, "language_info": { "name": "python" } }, "cells": [ { "cell_type": "code", "execution_count": null, "metadata": { "id": "VcEDUNH4G52w" }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "source": [ "a = 125_000.0\n", "dz = 1000.0\n", "n = int(a//dz)+1\n", "print(n)\n", "z = np.linspace(0.0,a,n)\n", "\n", "seg_por_ano = 365.0*24.0*60.0*60.0\n", "\n", "t_rift = 30_000_000.0 #duração da fase rifte em anos\n", "t = -t_rift\n", "\n", "beta = 2.0 #fator de estiramento litosférico beta\n", "kappa = 1.0E-6 #difusividade térmica\n", "\n", "G = np.log(beta)/(t_rift*seg_por_ano)\n", "\n", "dt = 1000.0 #passo de tempo em anos\n", "fac = dt*seg_por_ano*kappa/dz**2\n", "T1 = 1300.0\n", "T = z*T1/a\n", "\n", "tprint = t+5.0E6\n", "\n", "plt.plot(T,z/1000,\"g\",label=\"Temperatura inicial\")\n", "plt.ylim(a/1000,0)\n", "cor=\"r\"\n", "while t<=60.0E6:\n", " T[1:-1] += fac*(T[2:]-2*T[1:-1]+T[:-2]) + \\\n", " dt*seg_por_ano*G*z[1:-1]*(T[2:]-T[:-2])/(2*dz)\n", "\n", " t+=dt\n", "\n", " if(t>=tprint):\n", " plt.plot(T,z/1000,cor)\n", " tprint+=5.0E6\n", "\n", " if (t>0):\n", " G=0.0\n", " cor = \"k\"\n", "\n", "plt.plot(T,z/1000,\"b\",label=\"Temperatura final\")\n", "plt.legend()\n", "plt.ylabel(\"Profundidade (km)\")\n", "plt.xlabel(\"Temperatura ($^\\circ$C)\")" ], "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 323 }, "id": "U3f6fYLIHA0M", "outputId": "bd498a81-949b-4360-cd7a-d77e4efcb26e" }, "execution_count": null, "outputs": [ { "output_type": "stream", "name": "stdout", "text": [ "126\n" ] }, { "output_type": "execute_result", "data": { "text/plain": [ "Text(0.5, 0, 'Temperatura ($^\\\\circ$C)')" ] }, "metadata": {}, "execution_count": 2 }, { "output_type": "display_data", "data": { "text/plain": [ "
" ], "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEOCAYAAABmVAtTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOydd3gU1fvFz8zWbLLpAUIqCYGEVEjA0EIRSOgQeglFOtI7ijSlSBFRUQRBkPYVaRZERQR+FJUmXelFei9CQsqe3x8zGUJPIKHez/PMs2VmbmZmN3P23vO+75VIQiAQCAQCAJCf9QEIBAKB4PlBiIJAIBAINIQoCAQCgUBDiIJAIBAINIQoCAQCgUBDiIJAIBAINJ47UZAkKUGSpH2SJB2UJGnQsz4egUAgeJWQnqc8BUmSdAD2A6gK4ASAzQCakdz7TA9MIBAIXhGet55CKQAHSR4mmQrgfwDqPuNjEggEglcG/bM+gLvwAvBvltcnALyWdQNJkjoC6AgA9vb20cHBwTn+I1u3HgXgD+AWgJNwMt1E4dBQQJIe76gFAoHgBWLr1q0XSHrcb93zJgqPhOQ0ANMAICYmhlu2bHmcRmCUqyINEwBE4eqtP7F1W18UkDZi9TffILhBg9w9aIFAIHiOkCTp2IPWPW/DRycB+GR57a2+l7tIElL5K/7edQsmtFb/zHqc4TcIaTgIsiRhRIUKSD93Ltf/tEAgEDzPPG+isBlAkCRJhSRJMgJoCuC7vPpjwWGv4UbGTEwaPABGDAEQD2AviA8x/P92w5A/P16zt8fujz8G0tPz6jAEAoHgueG5EgWS6QC6AfgZwN8AFpLck5d/Uyfr0Gv0PPzfoXj0DQ+HDl+qh3AIQF9supmO8B49YDUYMKZUKSRv3ZqXhyMQCATPlOcqJDWnPLan8ACu3bqGUV+0AgYcwPib40DUBHAEwCAAC7Xt4kwmTOzeHTGDBwOurrn29wWC3CAtLQ0nTpxASkrKsz4UwTPGbDbD29sbBoPhjvclSdpKMuZ++whRuA//2zkfvw9rh+vLyuJLTAQQCeAPAH0BbNS2cwHQu0gR9Hj3XTglJgL6F863F7yEHDlyBFarFW5ubpBERN0rC0lcvHgR169fR6FChe5Y9zBReK6Gj54XmkY0R58v/8GlD67jf0VKoALaQPG/NwD4BkAAAOAygKH798O5SRPEmc1Y06QJuHv3MztugQAAUlJShCAIIEkS3NzcctxjFKLwAPyc/bC450YcmD4Czs3nYKd9EQRgKIAEKHbHB1D6CgrrMjJQaeFCuIaHY4inJ86OHQtcuvSMjl7wqiMEQQA83vdAiMJD0Mk6DIkbgkHvb0DLofnQp9S72IQgOGM2gB5QzOjeAIzaPlcAjDpzBgUGD0ZZNzcsL1cOGd9/L6KXBALBC4EQhWwQ6x2LdT13YNPgVujR7gy2evbAXETBiD+h9Bj2Amh4z34bAdTasAEudeqgn5MTjnXuDPz991M+eoHg6XLx4kVERUUhKioKBQoUgJeXl/Y6NTX1WR/eHaxZswYbN2589IaPSfv27bF378NLt02dOhVfffXVA9efOnUKDRvee3/JytGjRxEWFvZYx3gPJF/YJTo6mk+bBbsW0P1dR46sYmS6XseBcnVK2EGABNYTeI0AHriUArgwMJC3Pv6YvHz5qR+/4OVn7969z/oQNIYNG8bx48c/02NIS0t74LrHOb6HtfesOHLkCENDQ++77n7fBwBb+ID7qugp5JCmYU2xpdtO/NK8FIp1zkDr4K1IQXHUMXaFYkD/AWABlNpK97IJQONDh+DavTt6uLnhn5o1gZ9/BjIynto5CARPm61bt6JChQqIjo5GfHw8Tp8+DQCoWLEievfujZiYGISEhGDz5s1ITExEUFAQhgwZAkD5FRwcHIwWLVogJCQEDRs2xM2bNx/Zbq9evRATE4PJkyfj+++/x2uvvYbixYujSpUqOHv2LI4ePYqpU6di0qRJiIqKwrp169CmTRssWrRIO24HBwcASo+ifPnyqFOnDooVKwYAqFevHqKjoxEaGopp06bd97wrVqyIzAhJBwcHvP3224iMjERsbCzOnj0LABg+fDgmTJgAADh48CCqVKmCyMhIlChRAocOHbqjF3D06FGUL18eJUqUQIkSJfKklyNiKB8DP2c/rG69GmMCxiDCbTj6/u2MxT/NwYWM+UiQB2JHWk8A9QF8BGAUgKv3tHEDwMc2Gz7+8UdE/vgjejs7o1H79rB06AAUKfJ0T0jw0tLrp17YfmZ7rrYZVSAKHyZ8mO3tSaJ79+749ttv4eHhga+//hpvv/02Zs6cCQAwGo3YsmULJk+ejLp162Lr1q1wdXVFYGAgevfuDQDYt28fZsyYgbJly+KNN97Ap59+ip49ez603dTUVO2GfPnyZfzxxx+QJAlffPEFxo0bh4kTJ6Jz585wcHBAv379AAAzZsx44Hls27YNu3fv1sI7Z86cCVdXVyQnJ6NkyZJo0KAB3NzcHrj/jRs3EBsbi1GjRmHAgAGYPn26JnyZtGjRAoMGDUL9+vWRkpICm82Gc1nK7eTLlw8rV66E2WzGgQMH0KxZM+R2WL4QhcdEL+vxToV3UCWgClosaYE5vkfxy6ZgbP+/t7A5/yLUOd8dZ2x9AbwBYASAzwDc32zeAaDNlSvoPGECWk6YgK4RESjevTvQuDHg6Pj0TkogyANu3bqF3bt3o2rVqgCAjIwMeHp6auvr1KkDAAgPD0doaKi2LiAgAP/++y+cnZ3h4+ODsmXLAgBatmyJjz76CAkJCQ9tt0mTJtrzEydOoEmTJjh9+jRSU1PvidvPDqVKlbpjv48++ghLly4FAPz77784cODAQ0XBaDSiVq1aAIDo6GisXLnyjvXXr1/HyZMnUb9+fQBK4tndpKWloVu3bti+fTt0Oh3279+f4/N4FEIUnpDSPqWxvfN2dF/RHWEOX6F/RFGMXnwKp21tMa/ICnTY3wHJ+AhK6YwBAL59YFspAL4A8MXOnQjp0AHdu3ZF8/r14dSpE1CxIiCL0T5BzsjJL/q8giRCQ0Px+++/33e9yWQCAMiyrD3PfJ2uRu3dHVopSdIj27W3t9eed+/eHX369EGdOnWwZs0aDB8+/L776PV62Gw2AIDNZrvDGM/a3po1a/Drr7/i999/h8ViQcWKFR+ZD2AwGLTz0Ol02rnlhEmTJiF//vzYsWMHbDbbfYXjSRF3mVzA0eSI2fVmY0GDBZjmdQZ+Hf7DviZV0OLAN/ivQCsMCuwJGRkAlgFYDaDEI9v8G0DXtDTkW7gQrV5/HRu9vMBhw4AjR/L4bASC3MVkMuH8+fPazTstLQ179uSspNnx48e1/efPn49y5cqhaNGi2W736tWr8PLyAgDMnj1be99qteL69evaa39/f2xV65t99913SEtLe2B7Li4usFgs+Oeff/DHH3/k6Hzuh9Vqhbe3N5YtWwZA6WFleidZ/66npydkWcacOXOQkQdepBCFXKRpWFPs6LwDAf5RCA75FcNGVQFdnTHm0Ee4/nof1HEbCiAUwFYAs6HMKfRwUgHMAVD2zBkUHTkSHwYE4GL58sCcOcBdXxiB4HlElmUsWrQIAwcORGRkJKKionJskBYtWhRTpkxBSEgILl++jC5dusBoNGa73eHDh6NRo0aIjo6Gu7u79n7t2rWxdOlSzWju0KED1q5di8jISPz+++939A6ykpCQgPT0dISEhGDQoEGIjY3N0fk8iDlz5uCjjz5CREQEypQpgzNnztyxvmvXrpg9ezYiIyPxzz//PPD4ngRR+ygPSLelY8y6MRixdgQKWbzw65lq8Pv4K8BiwfHEtqgzzxc7bnUGYAMwEcA4AP9lu309lKyIjnZ2qNi8OaQ33gBKlxYzxwkAAH///TdCQkKe9WHkGkePHkWtWrWwW5SQeSzu930QtY+eMpkm9Lq265Bh0CHAcSY+nt4BjAiH78xJ2P7aUqxrOw75pW8BvAPgAIB2yO7HkQ7gawCVk5MROGMG3i9bFmeDgoD33wfUkDyBQCB4HIQo5CGZJnSL8BbocWQKyrVMw4UPRgE7dqDc/NE48/ZefFFlCMw4BMVi/gtAlWy1ndm/yyzs7XXoEBoMGoSfvb1hq1kTWLIEeM6yRwWCx8Hf31/0Ep4iQhTyGEeTI76q/xXmJ87H7ot7EZAyFosXvwfUqgW89x7anf0ON/53Er0CBkKGPYCVAH4A8Ojuv5TlMdPGTrDZUOinnzCqQQOcKlAA6N0b2Lkzr05PIBC8ZAhReEo0C2+GHZ13ILJAJBqu744WzU248c184NIlyM2aYlLNZFxeuAcJ1uEAygLYCeATAO4PbRe4UxwA4ITNhiEAfC5fRr3Jk/FjZCQyoqOBTz8FLl/Og7MTCAQvC0IUniL+zv5Y3Xo1RlYcia93f42wk2/hj59nAl27Ap98Ase+3bBiTnH8897XKKr7AkAnAAcB9ANguqe9zCEk212PhCIQNgA/ShJqAvDfuRMj33wTJwoUAJo3B1auBGw2CAQCQVaEKDxlsprQEiSUXVwdwxu6I33dWsBqBerVQ9Fdq/HPvqr4vtG7cMU6AOOhZC48vFIiszzKmYlu6o3/XHo6hgHwTU1FnUWL8EO1asjw9weGDweOHs2DMxUIBC8iQhSeEVlN6BFrRyDun4E48tsSYORIYOlSoGRJ1Erww8UjYRgb0Q9GXIcy69s6ACXvaS9z6Ej7QDN7AeoUoVmDVVfabKgNwP/cOYwYMQInChUCqlQBFiwAxLy+gidElM6+zbp16xAaGoqoqCicPHnykSWwH0bW4np5iRCFZ0hWE3rP+T2InFkSc+sWAnbsAMLCgHbtgLZtMXBRJ1xffRVJboMgoTCUWqtzAHhrbWm9hLtzFTJT6dWJu2UANjUL8kp6OoYD8AVQd+NGLG/eHBkFCgDdugHbtuXZeQtebtzc3LB9+3Zs374dnTt3Ru/evbXXRqPx0Q3kMg8rJ/E4opCT8hTz5s3D4MGDsX37dnh5ed1RgfV5RYjCc0Cz8GbY2XknIgtEImlpElrsfRdXf/oW+Pxz4K+/gPBwGDeuw1cnR+DExz8j1jAeylDSPijF9m5nNWYmI2aKg/aoputLOh0ydDrludqb0ANYmZ6OWgAKpaVh5NSpOBkdDRQvDnz8sZhWVPDEvIqls7/44gssXLgQ77zzDlq0aHFHCexZs2YhMTERCQkJCAoKwoABA7T9unTpgpiYGISGhmLYsGG5/VE8mgdNtPAiLM9ikp28JC0jjSPXjKRuhI7+H/pz/bH15KlTZMOGJECGh5N//EHeuMH1Hd6jl7RAndznJIHWBKR7JvWRZZkAaMh8BCipi8XOjgBo1emoU7d3NRqV/SSJ9ZycuAJghtFINm1KrlxJZmQ868skeARZJ1Xp2ZOsUCF3l549s38sw4YN47hx41i6dGmeO3eOJPm///2Pbdu2JUlWqFCBAwYMIEl++OGH9PT05KlTp5iSkkIvLy9euHCBR44cIQCuX7+eJNm2bVuOHz+eqampD223S5cu2nFcunSJNpuNJDl9+nT26dNHO76sk+y0bt2a33zzjfba3t6eJLl69WpaLBYePnxYW3fx4kWS5M2bNxkaGsoLFy7cc/5Z28s6Ec6XX37JQoUK8cqVK0xOTqavry+PHz9+R7vp6emsUKECd+zYoZ3T5s2bs3/xVcQkOy8wmSb0+jfWQ4KEuFlxGL7vc6R/vQBYtkz5xV66NDB4MMpO7IETJ+Lweem+sOA4gFkANgOIu6NNreJj5hs6nWJEA0hLTgYAyBYLMgAYZVkb87WXJKxMTkZ1AIEmE8YsW4azVasCgYHAu+8C//6b15dD8JKQtXR2VFQU3nvvPZw4cUJbf7/S2SaTSSudDeCe0tnr16/Hvn37Htru3aWz4+PjER4ejvHjx+e4IB9w/9LZmRPmZJbOzgmvv/46nJycYDabUaxYMRw7dgwAsHDhQpQoUQLFixfHnj17HjmdZ24jSmc/h8R6x2J75+3o9mM3jFg7AisPr8Tc+nNRaO9eYPBgZUhn6VLgs8/QceNEvPHnZvSs3w9TT/eADWsBLIZSpvuw1ubd4iDrdEjPyIAMIFWtEmlnb4/rN25AstlgliScT02FDkByRgbeSknBUJ0O9W/eRKehQ1F52DBICQmK71G7NvAMxooFj+bDZ185W5TOfgBZzzWzlPaRI0cwYcIEbN68GS4uLmjTpk2O231SRE/hOeUOE/rcHkROjcTco98BU6YA69cr4au1agFNm0Lv74spJ8fj7MzfUck8CkACgL0A3gfgCAlZjGjVZKYkKT0GnQ5pmb7D9euw2WxwsFjwn3qTdzWZcEUdv3XT6fDjtWuoAiDYxQUTN2zAxYYNAW9voH9/4J9/nt4FErwwiNLZ2efatWuwt7eHk5MTzp49ixUrVuRKuznhqYuCJEk+kiStliRpryRJeyRJ6qm+7ypJ0kpJkg6ojy5P+9ieR7JmQictTUKLJS1wtUSoYkBnhq+GhABffgn3No3x25V+2NrjEwRKC6AkvR0A0RGZH7WW4JblF1g6CZ1OhxTVgNbfvInk5GSY9Hpk2NvjFgAnkwkZaWm4kZICqyThZkYG+l27Bi+DAUlmMzZMmgSGhADlygGzZgE3bjzlKyV4XhGls7NPZGQkihcvjuDgYDRv3lwbMnuqPMhsyKsFgCeAEupzK4D9AIpBqR89SH1/EID3H9XWy2Y0P4y7TegNxzcoK/7+myxfXjGiK1Ui9+9X3j97lnMq96cVa1UzegeBSpoBLamPevXRaDBQkiTqdDpaDAYCoBtAWTWdPT09CYAmvZ5+JhMBUAcw0GqlndlMAAwvUIBT8uXjVYC0WslOncgtW0jV4BM8He5nLL7IZDVoBTnnuTeaSZ4muU19fh1Kqq4XgLpQZp6B+ljvaR/b88zdJnT5L8tj+JrhSC9SGFizBpg2TcktCA8HxowBXFzQctU4XNnhgj6+faGDFcBvAJYCCLoj+xkAkJYGktDrdEjJyIAsy0g1m2ED4Eziohrq55E/P06pY66+VivOXr+O5JQUeBuNuJKWhjfPnYOXnR06e3pix6xZQEwMUKKEMux15cpTvGICgeCxeJBaPI0FgD+A4wAcAVzJ8r6U9fVd+3QEsAXAFl9f31zQ0RePqylXmbQkiRgOlplRhocvqWFyp06RDRoovYaICPLPP5X3bTZe/t8yJlhGErhG4BaB8QScKElKGKtefTSpPQezyUSdTkdZluni4EAAdAFoVdd7FShAJycn5bmTE711OgKgoyQxzMODJjW0tXShQvzKz4/JAGlnR7ZqRa5bJ3oPecjL1lMQPBk57Sk8S0FwgDIvZaL6+spd6y8/qo1XafjofszfOZ+OYxzpOMaRc3fMvb1i2TLSy4uUJCWo/Pp15f3UVO4YOIlF5BkEMgicI9CRgEydelM3quJgp978LXZ21Ov1lCSJ7s7OBEAndWgJAPO5urJgwYJKjoODA4s5OChDTgDDnZzo7eGhDEU5ObF/ZCQP2dsrX7vgYHLiRPL8+Wdz8V5i9u7dq8XkC15tbDbbiyEKAAwAfgbQJ8t7+wB4qs89Aex7VDuvuiiQ5JHLR1huZjliONhicQteSb6irLh6lezaVREGX1/yhx9u73TxIufHD6Ij1qh+w3YCFZUeg16v9BhUcXBQb/72FgsNqu+QX73ROwAsqK53dnBgUFCQsq/RyCgPDzqr6wKMRpbw9aUsy5QkidXDw/l9cDDTATIzMW7VKpEYl0scPnyY58+fF8LwimOz2Xj+/Pk7Eu4yeZgoPPU5miUl4Hg2gEske2V5fzyAiyTHSpI0CIAryQEPagd4fudoftpknRPax8kH8xLnoYxPGWXlxo1Ahw7A3r1AkybA5MlA/vwAANvf/6B/rS8x+XAXZMAfwGLI0kDYeAgGgwFpaWkwyzJSbDZYAVwH4Gi1IuXWLaSmpsLL0xMnT5+GBUABKFkRFpMJAUFB2L9/P1JTUxHh7Y3/zp7F4bQ0uEgSQj09sf/WLZy7eBF+BQuis78/2u3eDY9r14DChZVjbdMGyJfvGVzJl4O0tDScOHHiqce3C54/zGYzvL29YVBrn2XysDman0UvoRyUX5A7AWxXlxoA3ACsgjJh8a9QREH0FHLAxuMbWejDQpRHyBy2ehjTMtKUFbdukSNHKr/KnZ3JL764Y0z/4uIfWcV+NIHrBFIIjKYsO93RczCrZTKcMnsGTk40qr6Br48PJYBmSWJRdejIoNMxKiKCzuqQUxFvb8a4uFBSo5ZKOzuzROHCypCV0ciWZcpwY2QkbQBpMJCNGpG//CJ6DwJBHoDnbfgotxYhCvdytwl95PKR2ysfFL5Kkunp3PzWFAbIc9UhpVMEWlGSdJo4yJJEkyoOzurwkpurK01qiGohf3/qZJkGSWJoFtM6KiyMvr6+BMCC+fOzgp8fHdV1oUYjqxYtSgfVzC4eEsIZVavyhouLcpyBgeTYseSZM0/9WgoELytCFF5BHmhCZ2SQn39OOjmRZjM5ejSZmnp7/dWr/KL6cNrjd1Uc/iTwGiVJ0gxng15PgyRRAuiiikM+Dw+a1XyFwoUL02AwUCdJjAA0AQgvWpTh4eFKlJLVyiqhoQxQBaWAJLG2vz+LFiqkRDq5uLBvjRo8WKqU8jXN7D38+qvoPQgET4gQhVeUI5ePsOyMsvea0CR58iSZmMh7wldV0vYfZOeAoZRxQhWHryjLPoqRrN7I7UwmymoCXKY4eBYooIlD0aJFaTaZKAGMAphPFYeifn4sXbq0JjSvlyzJ0qp5bQZY08OD1YoX10SoRlwcf6xfnxmZvYfChclx40i1OqZAIMgZQhReYR6YCZ3J0qW3w1d79CCvXbtj9ellv7KswyQCyarn8BYNBqsiCmrpbUeLRbuhO6ni4OPtTbPZTEmSWKxYMVrUbaJkmb6qOPjlz8/KlStr68qULMmaRYrQrK4vbzazeXQ08+fLRwAMDAjgBy1a8HJsrPLVNRrJZs3ItWtF3oNAkAOEKAgebEKTd4av+vjcGb5KkhkZXDt4GgvKS9Vew2EC9WkwKEazvb29MuRjVcTCqiaxZfoMZrOZsiwzPDyc1kzvQK9nUfXmX8DZmdUTEuih9hbCw8LYJDaWnqrJXUSS+EZICGOjopTcCYuFnRo14q7mzZVhMIAMCSEnTyYvX366F1YgeAERoiAg+QgTmiQ3bCCLFVO+Fk2a3Gvu/vcfJ1QdQxN2EiAlrCRQTMtfyPzF7+boqDwCdFDFIahwYU0cihcvrkUllTCZGKWKg6vFwjo1a7KwGpXk5+fHpBo1WNzVVVkP8A0vLzaNi9OGqCrFxXFJt25My/Qe7OzItm3JTZue2nUVCF40hCgI7mDeznmaCT1v57w7V966RY4YoQzNuLiQM2feMzSTfOAIG/mNo4SLBNIITKZO505Jkmg0GmkwGCjLMl3VnoMnQIskUZIkhgQH02g0UqfTMSYmhq7qDb+ExcLSqjg4GI2sX7MmS5YsqYiBqytbN27MmkFBlKDMHtfYamWPihXp66P4HL6+vhzbowcvtG5NZmZNR0eT06eT//339C6uQPACIERBcA9ZTeiWS1ryasrVOzfYu5csV075ilSuTB44cE8b+xetZqjdlwTSCZwn0IFmszKU5ODgQFmWaTKZ6Kz2IPzUMFW9OpRkNBqp1+sZGxurDR1FOzqyEpQqriadjonVq7NatWqKZ2E2s0WzZmxdrhwd1NDYijod3ypVihVjY7Vt2rdqxZ1vvUWGhSnH7+Sk+CWiJpBAQFKIguABpGWkccSaEdSN0LHQh4XuNaEzMsipU0lHx/uHr5Kkzcave86kk/R/6pDSVgJltLyDzJ6As7Mz7dUhn0A1gc1kMDAqKop6vZ5Go5Hly5dn/vz5CYAxrq6srkY26SWJ9atUYYMGDbReSGJiIns1bkxv1ewOBji0aFG2rV5dM8ArVarEpaNHM71ZM6Xnk5mf8c03956HQPAKIURB8FAeakKTSvjq/aqvZiHjxk12L/0hZfxLgJQxh4CndoPOFIf8+fLRZDBQD2hGs9XOjtHR0dTpdDSbzaxcuTILFCigiEO+fKwnyzRByZSuW74827RpQ6s6NBUfH88hPXqwuComHgAHeHjwnTp16OPtrZjdhQpx4ogRvDxsGOnvr5yHpyc5bBh54sRTucYCwfOEEAXBI7macpUtl7S8txx3VrKGr2atvpqFs7sOs2z+aVTKZVwjMIBGowP1ej31er1WbtvHx4eyJNEiSQzONJqtVpYqVYqSJNHe3p7x8fFaBdZSBQuyicFAe3XbGjExfLNrV+ZTw1VjY2P57rBhrBEWpoXHdjCbOTk+nuVfe02Lkur25pvc//nnZI0aynnodIrg/fabCGsVvDIIURBkm6wm9B2Z0JncXX11+fL7trNq2irmN6xQh5T+IZBARzUqydHRkRaLhZIksZC/vyIIsswi6g2/oJsbX1Nv5E5OTqxVq5YmDq/5+jLJzk6rwVQtIoJ9e/dmITUTOiQkhKNHjeIb8fE0yTIlgPUkiTPLlGGrWrVoNBopSRJr1arFlbNn09a3L+nqqvwrFCtGfvLJPbkaAsHLhhAFQY44fOkwy8woc/9M6Ew2bFByAwCl9PX9ahPZbBzR7CsasJ8AqcO3BApp4aj58uVTpv+0WOin1kbylmUGqjf8QC8vTRzc3NxYv359TRxKBwayvZMT3dVtKxUpwncGDWJERIQWjTRq1CgO7NiRrmoGdhmAM4sW5dDGjbUeRnh4OGd89hmTp00jY2KU87FayW7dlFpRAsFLiBAFQY7JakLfNxOaJFNSHhm+SpJXz19jQvBMKhnRyZQxgrLsoJnRXl5e2o2/oDoXdFGdjt7qDT8yKEgLT/X09GSTJk20OaPLhYSwa758LKBuG+fvz/fefptlypRRPAYPDw4fPpzjhg9nIdXXKAJwirs7pzVtygh1uMnDw4NDhw7lmeXLyaSk28b0668rw2bp6U/hqgsETwchCoLHJtOE1o3Qcfjq4fea0OSd4at3V1/NwuaVe+nn8J06pHSEQF3a2zvQYDDQYDBoN3o/Pz+6urgo+Qs6nVYzqXREBEWThvkAACAASURBVIsXL04A9Pf3Z4sWLbRopcrFi7OXnx+91G3LeHlxwtChjI+P14ahBg8ezOlTpzImIEAxvQGOMpm4rF491nr9dWaW8X7jjTe4a80actQoJcMbIP38yPffJy9cyNsLLhA8BYQoCJ6IrCZ02Rll782EJrMXvqry4aAVNMm71SGlFQSCtDwFFxcXuqiCEBISQjs7O+plmaV1OlrV/IUqsbEMDQ3VtmnTpo22f7XYWA4IDqaPKg6v5cvHj4YOZYMGDShJEu3s7NizZ09+/fXXjFeHphwA9pEkrqlcmV0aNNAipqpWrcqfli+nbdEismJF5d/FbCbbtSO3b8/biy4Q5CFCFAS5wkMzoTO5u/rqA8pNJCensX6ZrwlcIXCLeowhYNFu7t5qQT2DwcCIiAjKskyr0cjyOh2NAE2yzNqVKmlTgEZHR7NDhw50c3MjANaqVIlDoqLop4pDSVdXTh0yhK1ataJOp6PRaGSnTp24fPlyNq9blzpJogFgW4C/R0ZydMuWWs8lLCyMM2fOZMqWLWTHjkopDYCsUIFctIhMu0/vSSB4jhGiIMg1sprQLZe0vL8JTSrj8AULkrJM9up13/BVktz11ykGuv+gDikdp4QGNJvMWqRSYGCglvwWpo7/ezk4sHzmZD9GIxtVr05vNSehYsWK7Ny5s2ZmN6hRgyPLlqW/WoOppJMTZwwezE6dOmnlNlq1asVff/2V3Tp2pFmvpwSwPsB1BQtyVsuWDFd7JQUKFOCoUaN48eBBcvx4ZUgJUKKw3n+fvHgx7y68QJCLCFEQ5CrZMqFJ8soVskuX2zfOB4SvkuSn4zfSTq8U2tPjFwJFmM/Dg3q9nnZ2dlrIqZ+fnyYUoW5ujFVv9t729mxWp47WU6hduzY7d+6sldto2bgxx1SpQn9VTEpZrfyqf3/27NGDdnZ2lCSJTZo04Zo1azjkrbforFZ+rQzwZwcH/ty0KRMqVWJmldZu3brx0P79ivhlDi3Z2ZGdOpF79uTRlRcIcgchCoI8IVsmNEmuX//o8FUqQ0oNE74ncJnALRoxmoBF6wXkz59fy3SOjIzUwkrLeXoyXB0mCnN1ZVKDBrS3t6csy2zRogU7depEs9lMvV7PDm3acFzt2vTTKdOMvmax8OvevTlwwAAtGqp+/fpct24dJ4wfz4Lu7srwFMDFej2316nDNvXqaeU2GjVqxE2bNpE7dpBvvEGaTMp5VqtGrlghZokTPJcIURDkGdkyoclsh6+S5O7d5xhY8Gd1SOkoZdSjyWDUIo2CgoJob29PvV7P0qVLa8+r+/hoE/jEeXmxeWIijUYjzWYzu3Tpwnbt2tFgMNBkMrFnt26c0LgxfdU5G8qYzVzSrRvfeestLeu6du3a3LBhA6dNm8ZANY8iWJI4C+DRuDgObNJEG+aqUKECly9fTtvZs0rUkqen8u8VHEx+9hl540befQgCQQ4RoiDIc7JlQpPZqr6ayZSPttGs30uANGI5gUDmd3envb09dTodw8PDKUkSXVxcWLZsWeW5oyPr+/rSRY1USixShA1r16YkSXR2dubAgQPZsmVLyrJMBwcHvjV4MD9o1YreBoNyczeZ+GOXLnx36FCtXlONGjW4YcMGLliwgBGqv+Any5wC8Gx4OCe0aKH1ZsLCwjh79mzeun6dnDv3dkKcqys5eLCotSR4LhCiIHgq3D0n9D3luDO5O3x1zJgHhq8mJ6exbo2VBK4SSKYZwwmYtfIYLi4uLFq0KAGwcOHCLFGihGJQ+/iwvpcXjVDqILUrXpzV1FwELy8vvvvuu2zQoIGWNDdmzBhOaNeOBVRxqGowcFWHDhw9dKjmUyQkJHDjxo38/vvvGVuqlGI+63ScAPCStze/atGCYcWKadFTH3zwAa9fu0auW6dEZMkyqdeTLVqQW7fm4SchEDwcIQqCp8bd5bg3Ht/44I2zhq9GRj50trStW0/Tt+AaAqSMA9QjgWaDgb7qsE7hwoVvl8AoXZoBaoJamago1sisoCpJ7FW2LEtGRxMAixUrxg8++IBVq1bVxGLKJ5/w/U6d6KGKQy29nuvatOHYoUPprvoLmeLw22+/8XVVaNz0er4L8JKjI5c3acI4dX4HFxcXvvPOOzx37hx5+LASieXgoJxzxYrkd98J30Hw1BGiIHjqbDi+QTOhR6wZ8WATmrwzfPUB1VczGTduK436A0oeGRYR8GYBV1dtqCc6OpoWi4VGo5FVq1ali4sLJUlivYoVWVoNUy2q03HQ66+zsBrFFBcXx08//ZSx6o28SJEinDNnDt/r3JnOqjg01Om4OSmJ7w8Zcoc4/Pnnn9y4cSNr1qxJAHTU6zkE4AWjkb/Xq8d6quDY2dmxe/fuPHbsmBKVNWHC7WzpIkWUntPNm3nwSQgE9/JEogDAG0A/AN8C2Azg/wB8CqAmAPlR++flIkTh+SbbJjSZo/DVq1dTWLXKbwRuErhOC/oS0LNIoULU6/V0dHRkTEwMM2slxcfHU6fT0Wq1snX16iysTvZT0Wjk4Ph4LYqpcePG/Oyzz7Rs6ZiYGC5btozvdOpEB72eMsBWssxdTZty7KBB2rBSzZo1uWXLFm7btk0bknIwGDhAp+NZgHurVWObWrW08uGtW7fm3r17lSGz+fOVaUMB0t1dmePh3Llc/ywEgqw8tigA+BLALwB6ACgDoDCAMACJAD4GsBFA3MPayMtFiMKLQbZNaPLe8NWzZx+46dq1x5nP4w81t2EH9ShNi17PIHXoKDAwUPMbihcvznLlymnvt69Rg646HSWASQ4O7F2tGi0WCw0GA3v06MGPP/6YPur8z/Hx8Vy1ahX7vvEGzTodDQDflGXub9iQo/v21XopdevW5fbt27lr1y42bdqUkiTRYjCwr9HIMwCPlSnDHvXqaXkRiYmJ3LJlixKFtWYNWauWct5msyKQDzHhBYIn4UlEIewR640ACj9sm4fsqwPwF4Af1NeFAPwJ4CCArwEYH9WGEIUXhyOXj9yRCf1AE5pUwleHD1fCV11dyS+/fGD4akaGjYMG/UGd7gQB0h7TCbjS28Xljgl4XF1dKUkSq1evziJFimhhpK1ff50GSaI9wP6urmxbsSJlWaaTkxNHjx7N999/X7vpN2/enBs3bmTnli2pl2XaARwE8Gi9ehzZvbsWytqwYUPu2bOH//zzD5OSkijLMu0MBvayt+cpgOfCwjgkMVHbvlq1aly7dq1yQnv3KrWVjEZlzooGDe47051A8CQ8l54CgD4A5mcRhYUAmqrPpwLo8qg2hCi8WKRlpHH46uGUR8iPNqFJ5QZZtqzyNX399Yf+cj59+jpLllxNII0SztEerQiA4QEBNBgMtFqtjIuLo06no7OzMxs2bEhnZ2cl27llS9ZUo5a8AY719GR11Yz29/fnjBkzOHjwYNrZ2dFgMLBnz57ctGkTWyQmUgLoDHA0wBM1a/Kdjh1ptVopSRJbtmzJAwcOcP/+/WzTpo0y3ajBwF4uLjwF8GpAAMcmJmriVb58ef7000+02WzkqVNKCKuzMzVT+scfxexwglzhiUUBQC31V/0lANcAXAdwLTv7PqA9bwCrAFQG8AMACcAFAHp1fWkAPz+qHSEKLyYbjm+g/4f+2TOhMzKU5K9shK+S5KJF+2i17iBAmrCGOgTTqtezmFo4L2vYamhoKOvXr09Zlunq6spePXqwuJ+fkukM8CN/f0apZnSpUqW4ZMkStm/fnrIs09HRkaNHj+amTZtYWy3PXUCS+CnAU9WqcUCrVrSzs6NOp2O7du147NgxHjx4kG3btlXEwWhkr3z5eBrgjQIF+FHduvRW55WIiYnht99+q4jDtWvkxImkt7fy7xoeruQ/iCJ8gicgN0ThIIAIAFJ2ts9Ge4sARAOoqIqCO4CDWdb7ANj9gH07AtgCYIuvr28eXjZBXnIl+Ur2TWgyR+GraWkZTEpaS0m6ROAWnfAuATMLOzuzgBqeGhcXp03uU7NmTZYtW1bpWYSHc2C/fvRUh3aaA/wgMJBeavXWhg0b8qeffmLt2rW1fIQvv/ySa9euZTk1eqmwLPN/AE9WrswejRvTaDTSaDSyR48ePHPmDA8ePKj1HOxMJvb19eVZgCnOzpxWqxYD1ByMyMhIfvPNN8zIyCBv3SJnz1amDAVIf39yyhQRsSR4LHJDFFbnVqSR2uv4VH2eY1HIuoiewotPjkxoMtvVV0mlXEZAwDoCpE46SAdUoQQwxs+PBoOBjo6OjI+Pp9FopIODA5OSkrS8h/r16/PNzp1p0utpATgE4JDChWkxm2k0GtmvXz/+8MMP2oxwERER/Pnnn/nDDz9oVVWjdTr+CvBY+fJsX7u2NvXoW2+9xcuXL/PAgQNs1aoVZVmmxWzmwKAgXgCYZrFwdkICi6iGebFixbhgwQKmp6crPafvviNLl1b+fT08lLkrrjygWq1AcB9yQxRKAvgJwGDVC+gDoE929r1PW2MAnABwFMAZADcBzBPDR68ud5fjfqgJTd4Zvurnp4y1P4QJE7bRYDhMgLTK8wl40EOnY7haebVYsWJaTyE4OJitW7em2WymnZ0de/fuzcQ6dQiAvurwUFt/f0qSRHd3d06ZMoXz5s3TqrjGx8fzr7/+4ldffUXfzOglg4F/AdwXG8umlStrSW1jx47ljRs3uG/fPjZv3pySJNFqb893wsN5WZaZbjRyQZUq2tBXcHAw582bp4iDzUauXUtWr65cB0dHctCgBxYbFAiykhui8AuAJQBGABiWuWRn30e0WzGL0fzNXUZz10ftL0Th5SHHJjR5Z/hq8+YPDV+9fDmZFSqsJpBC4CJdpQ4EJEY6OjK/mxslSWJ8fDz9VE+hdu3aWkJaQEAAR40apdU9ipNlLgBYUa3YGhoayh9++IETJ07UzOv27dvzyJEjnDhxohb51Mps5jGAf0VHs4Y665unpyenTp3K1NRU7t69mw0bNlREw8mJo0uW5H8GAzN0Oi6sWJFhatRU0aJFb4sDSW7bRjZqpEQrmc1kjx7kv//m0icjeBnJDVF45FDO4yx3iUIAgE2qf/ENANOj9hei8PKRIxOavLP66iPCV0ly+fJDdHL6S0kH0G+gGcE0AyxbsCBlWaa7uzvr1q2rDSl16tSJISEhhJrBPGLECLq5ulKWJHaSZc7S6RigVkqtVasW//zzT/bq1YsGg4H29vYcOXIkT548yYEDB9JkMtGk17O/vT0vA/y/yEiWUYUmKCiICxcupM1m419//cVatWoRAPO5u3NSXByTzWZmSBIXlS3LcLXncI84/PMP2aaNUl/JYCDbtycPHcq9D0fw0pAbojAOQLXsbPs0FyEKLydXkq+wxeIW2TehyTvDVx9RfTUjw8a2bddpRrS7aQwBIwNkmaHqFJwlS5ZkXFycNrzUrVs3Wq1WGo1G9unTh126dKFOp6OLycQPJYlj9HpajUYaDAb27duX27ZtY2JiomZGz507l0eOHGHr1q0pSRJdLRZOcnRkCsDvihVjmDr8VLJkSa5atYok+fvvv2u1lbwLFuT0hASm2dszA+Ci115jWOHC2vEtXLhQMaRJ8uhRsmtXZW4HnY5s1UoRDIFAJTdE4ToAG4Bk5EJIam4tQhRebubumKuZ0PN3zn/0DjkMX929+xz9/NYTIPW6/XTVK+P95ezs6OzgQL1ezwYNGmiZzY0aNdKGd3x8fDhp0iRWUmdji3B25lKA7Y1GSpJED3d3Tps2jb/99huj1ZyHUqVKccOGDdy+fTurVKmiDE25u/MbFxemAZxVpAh91JyF+Ph4bt++nSS5atUqvqYONwUFBPB/iYnMcHRkBsCF0dEMUQUlPDycS5cuVUJZSSViq08f0mJRhpaaNhWzwglI5oIo3HfHXApPfZJFiMLLT45NaFKZs6B+fWYnfJUk33tvM/X64wRIV4f5lOBEd4Bxaiazn58fGzVqRIPBQBcXF/br14/h4eHazfvjjz/WopaaFSzIFQDLqYX0ShQvzrVr13L27NlaFdemTZvy6NGjXLFihTbvdNmAAP7h7s5kgBMCAuiiJsC1atWKx44do81m47fffqttXzwigj8lJdHm7Mx0gPOiolhEPYaYmBiuWLHitjicO6eY0A4Oijg0bkzu3JkLn47gRSU3egoj73otA5iXnX3zchGi8GrwWCY0SS5Zku3w1bNn/2NMzGoC6ZSk0yxgbavcYAEGqp5BtWrVWKZMGW2Yp3///nR0dKTRaOTAgQM5ePBgmkwm2tvZcWyhQpwD0Fud9rN5s2bct28f33nnHZrNZprNZg4dOpRXr17ltGnTtFnlmpUowWP58vESwP4+PjQZjTSZTBw4cCCvXLnC9PR0fvXVV/RXcxkqlS/PTV26kK6uTAP4ZUQE/dUhsLJly94un0GSFy6Qb79NWq3Kv36jRuSuXU/46QheRHJDFL4EMFh9boJSMXV4dvbNy0WIwqtFjk1oMsfhq3Pn7qWdnTLbm6P1V9rpfGkHsJpeT5NagbV169b08PDQoowaN25MQCmJMX36dNZRQ1iDfXz4fWAghwA0SRLtzWaOHjWK+/fvZ9OmTTW/Yd68ebx69SqHDBmiCcZbVarwmocHjwFMUm/ybm5unDx5Mm/dusWUlBROnjxZK+PdqH597u/bl3Rz4y2An4aG0lNdV61aNW7evPn2SV68SL7zzp3isHv3k304gheK3BAFCUqdosFqeGqv7OyX14sQhVePHGdCZ5I1fLVZs4eGryYnp7FGjdUEbhC4Qu/8QwlIDAZYSi27HR0dzcaNG1OSJHp5eXHEiBEMDg4mANarV48zZ87UJvppXLo0N/j4sB6U+aMLe3nxxx9/5Lp167SSG2XLluWWLVt47NgxtmjRggBYIH9+zmzcmBnu7twKsLJaqjsoKEjzDq5evcqhQ4dq81R3bd+eZ956i3Rz402A40NC6JalUN/ff/99+0QvXlR6DpnDSk2bCkP6FeGxRQFAiSzLawC2A5iS+d7D9n0aixCFV5e5O+bSOtqafROavF191WDIVvjqqlVH6eKylQBpsdtMF/sSlAEmAHQ3GKjX65mUlKT5CzVq1ODAgQNpsVhosVj43nvvcdiwYTSbzXRwcOD4Ro34g5sbi6jiULt8eR44cIAzZsygh4cHJUliu3btePbsWf7xxx/apD8loqL4f5060ebiwh8Ahjg4KPkScXFK6W2Sp0+f1iKiHBwcOPLtt/nf8OGkiwuvAhxWtCgdLBbKssx27drx36x5DBcuKJ6DxaIMtbVqJUJZX3KeRBRWP2T57WH7Po1FiMKrTVYTOmlJUvZMaDJH1VczMmxs1y4zfDWZ3l6fENDTR5IYr97cA/z82LFjR9rb29NisfDtt9/W8gzCwsL49ddfa69DixXjqs6d+b7FQnuAJlnm8O7defr0afbp04d6vZ5OTk6cPHkyU1NTOX/+fHp7exMAmyQm8lifPkxzdORnAD1MJkqSxNatW/PkyZMkyX379mmhsAULFuSMTz5h+vDhpJMTzwHsFRREo8FAs9nMAQMG8NKlS7dP9uxZsm9fJXJLryc7dRJJcC8peRJ99DwsQhQEj21C3x2+OnbsQ8NXd+w4Sy+vjQRIk2kPPfMnEACrAgyUJAJgYr16rFatGgGlmN348eO1cNb27dtz7ty5WsZ062bNuL1DBzZRjegAq5U/zJ7NvXv3auGqYWFhXLNmDW/cuKH1OMxmM4cNGMAbAwbwisXCgQCNskyLnR1HjhzJGzdukCTXr1+v9TTCw8P586JF5JAhpIMDjwJsFRBASZLo7OzMcePGMTk5+fbJnjyp5DkYDEquQ+/eYja4l4wn6Sm0fFjoKYBAAOUe1kZeLkIUBJmsP7Y+5yY0eW/4alZD9j4MGPA7ZfkMgTR6ec2mLFvobjCwHkA9QHdHR/bo0YOenp6UJImdO3dmjx49qNPp6OHhwS+++IKDBg3Swls/Hz+ev9Spw2C111E3KIhHtm/nkiVLNAFp1qwZT5w4wWPHjmmmtq+vL7+ZMYO2Pn142GhkQ1WYfAoW5IIFC2iz2Wiz2bhw4ULN20hISODudevI/v1JOzvu0OlYXRUtHx8fzpo163Z2NEkeOaJkSMuy4jsMG0ZezWZvTPBc8ySi0BPADgAzAbwJoDGAVgBGAlgLYDGAoIe1kZeLEAVBVh4rEzqTJUtIT0/lBti790PDVw8fvsygoP9Tkt70B+nn10wZ4zebWVy9ub8eG8ukpCRKkkRvb29OnjxZS0CrUqUKf/zxR1aoUIEAWLp0aW5Ztoxjo6JoAWgHcExCAi+fOsVhw4bRZDLRwcGBEyZMYGpqKtesWcOIiAgCYOXKlbl71SqySxeu0ekYpYpD2VKluHXrVpJkSkoKJ0yYQCcnJ8qyzC5duvDcrl3km2+SBgN/MxgYrYbERkVFceXKlXee8J49ygxwgDKP9IcfKv6M4IXliYaPoEybWRXAcACfA/gQQCcAvo/aN68XIQqC+5HVhM5WOe5MrlwhO3dmdsNXR4/eQp3uXwIZLFjwG5rNbrSaTGyi19MeoL3BwJ6dO7NYsWJa9M/YsWPp6OhIk8nEd999lzNmzKC7uzt1Oh379+/Pv5cuZX31Bh2s13N1v348tG+f5kkUK1aMa9asYVpaGqdMmUIXFxfqdDr27t2bV7ZvZ3pSEqdLEj0kiRLA9q1a8awaaXX+/Hl269aNOp2OTk5OnDBhAm/t20e2acMMSeJ8k4l+aqRS9erVufvuMNVNmxQPBlDmc/jqK2UYTvDCITwFwSvH4UuHWfqL0jnLhM5k3ToyOJjZCV89efIaw8PXECB1umMsXLgDATA2Xz5WUHsNr/n5sUe3bjSZTHR2dubEiRO1chlhYWFcsWIF27dvT0DJnl6+fDl/GDOGhUwmAmArJyeemz2b3337rZa01rJlS545c4bnz59nx44dKUkSCxQowDlz5tC2ezcv167N3uqQlpPZzMkTJjBNna1tz549TEhI0MJbv//+e9p27yYbNGAywPEWC53MZsqyzM6dO2uiovHLL2SJEsr1iYpSXgteKIQoCF5J0jLSOGz1sJyb0GSOw1c/+mg7DYYjBEhPz2V0cPCkndnMFvnz0xWgQZLYvW5dlitXTslErlSJU6dOpY+PDyVJYteuXblixQqtImvjxo156NAhvpWYSANAF4DTAwN5/ZdfOGTIEBoMBjo5OXHKlClMT0/n5s2btQl/4uLiuGvXLnLTJu4tXZpVVXEK8/Li6l9/1Y55+fLlLFq0KAGlXMfff/+t9AYqV+Z5gN2tVupkmVarlWPHjmVK1iGjjAxy/nylxwCQVauSaq0mwfOPEAXBK81jZUJncnf46sGDD9z0/PkbjI5eTSCDOt2/DA7uriS6BQayhp2dEglktXJwp050dHSknZ0dR40axW7dumlJcIsXL+a7775Lk8lEJycnTps2jbv++otxarnscgD3VKrEf5Yv1yqolixZklu3bmVGRganT59OV1dX6nQ69uvXj9evX6dt5UouLVyY/qo4NC1fnifUUNPU1FROmjSJTk5O1Ov17NOnD69cuUKuXElGR/NvgLWsViVCKiCAixcvvl1TiVTE84MPSBcXJQGubVvFvBc81whRELzyPJEJncPw1alTd9JoPKT2Gr6ns7MvjUYjk0qUoCdAGWCXyEjWVMNXY2JiOHfuXK3YXePGjblhwwZWrFiRAFi+fHnu3buXMz/7jK52djRAmR70Zps2nD9lCvPnz09ZltmjRw9eu3aN58+f14ajvLy8uGjRItoyMnhz3jwOc3WlCaCDLHN8165MVc/j3LlzbN++PSVJYr58+Thr1ixmpKWRCxaQAQH8BWCovb3Wy9l5d0G9S5eUHAejkbSzI4cOJf/77zE+KcHTIDfKXOQHMAPACvV1MQDtsrNvXi5CFAQ5Zc6OOTnPhM7kxAmyXj1mJ3z14sWbLFlyNYF06nQnGBramwBYPCyMiWqoaVGdjiNq16aHhwf1ej3ffvttDhs2jEajka6urpw1axa/+OILuri40Gg0cuTIkTxx4gSTGjVS/ACAq00mXu7Xj13atdN6G0uWLCGpzMcQFRWlZVsfOnSITEvjoVGjWEv1K0KtVq6dNUs77s2bN2v5DaVLl1YimG7dIj/6iGmurpwC0NVopCzL7Nq1Ky9cuHDniR8+rFRhBZRihLNmCTP6OSQ3RGGFGo66Q32tB7ArO/vm5SJEQfA4HLp06PFNaJJcvDjb4avTp+/Seg1eXt/TxcWXBoOBbWrXpq/RSAlgJzc3NlHnZcjMgM6sxpqQkMDNmzdr+QlhYWH8888/uXLlSgaoOQbtAF5yc+PvffsyQi25UbduXf77779MS0vjpEmT6ODgQLPZzDFjxii9gxs3+G1SEv3UENbWQUE8t2MHSTIjI4MzZ87USm907dpVyXy+coUcPJgXTCZ2lWXKkkRXFxdOnTr1zvwGUqk1VbKkcouJiSE3bMjpxyTIQ3JDFDarj39leW97dvbNy0WIguBxeSITmsxR+OrtXkMGdbrjjIhQew1RUWyq3vwDAY4JDWVBdSiof//+/OCDD2hvb08HBwdOmTKFy5Yto5eXF2VZZu/evXnu3DkOHDiQOp2O+Q0GLgR4KzCQ7ycl0c7OjlarlR9//DHT09P577//auUvQkNDuUG9Sd84fpxvlSqlmdnTqldnxuXLJMnLly+zW7dulGWZHh4enDVrluInHD9OtmrFHQDj9HoCyrwRGzfedQ0zMpSw1YIFlevUvLkom/GckBuisAaAG4Bt6utYAGuzs29eLkIUBE/K+mPr6TfJj7oROo5cM5LpGemP3ikrWcNXmzd/aDmIzz/fSYPhMAHSx+dburkpXkP7Nm1YyNmZANhJltlK9RaCg4O5ePFiVq1alQBYoUIFbtu2jV26dNGM399++43btm3TZner5+DAkwAPRUezmposFxsbq+UcfPvtVJ2bGgAAIABJREFUt1rpjS5duiimMsk9P//MOA8PAmAZvZ673n5b8022bdumDSnFxcVxT+bsbVu30lahAhcA9FLF4Y033uC5u6/B9etKNVaTibS3J0ePFslvz5jcEIUSADYAuKo+7gcQkZ1983IRoiDIDa4kX2Hzxc2J4WC5meV49PLRnDWQkqKUgMgMX50164Hhq+fP32Dx4msIZFCvP8KoqG5KhFJ0NJvXq0cALAzwA6uVvm5ulCSJffr04WeffUYnJyfa2dnxgw8+4KpVq1hYnaO5Y8eOvHjxIsf9P3vnHR1F9b7xd7Zls5tk03tCGhB6AIEQqhRRFAQbKkgTLPiVoqggCAoKgmJHmlQLKB2VIkoTFSlCAAtSpIceIJCEJLuf3x8zGUL5KcoGUec5555tcycz95zcZ+593vd5R47Ebrfj8vVlgr8/bhHer1OHkKAgrFYrgwYNIj8/n5ycHPr06YPJZCIqKorZs2cD4PF4mDxoECEWCxYR+gcHk/vxx+Dx4Ha7GT9+PEFBQVgsFvr166f6LHk8MHcuOUlJPCWCRVEICghg7Nix52tGF2PXrvOaTErKHyYHGig9eCX6SNMRKolIZRGxXmm/0mwGKRjwJoozoV3DXX9ehAbVDiIjQ/23atbsd8NX33prExbLHkTcJCbOJigoCh8fHx599FESoqJQRHhMhG7a03v58uWZP3++ntmckZHBpk2b6Nu3LyaTidjYWBYtWsT27dt1+4xmiYnssts5YrHQXqv1ULFiRb777jtAFZWLheg77rhDd1o9euQInbVw12QRvqxWDTTLjCNHjtCpUydEhMTERBYvXqze0Llz8Npr/OjnR+PipL2aNdm4ceOlN79kCZQvr45Tmzaqx5KBa4q/TAoicsfvtd/rey2aQQoGvI2SmdB/yo67GG43vPuuWtXM1xdGjIDCy+dFZGXlULHiSkTAZttO9erd9aif9vffj4hQzmzmLRHinU5MJhN9+/blvffeIzAwELvdzqhRo/jmm290K40uXbpw/PhxxowZg7+/P06Hg7fS03GL8LmfH7GBgSiKQp8+fTh79iwFBQW8/PLL6grD5WLChAl6HsJXX3xBikZKnUU43q6dnoOwbNkyypUrh4hq2KdnPR85gufhh5mmWW2YFIUnevcm52Ix/tw5NbTX4VDH6aWX1O8MXBNcDSlM1trnIpItqgHebBE5ISKf/V7fa9EMUjBQGrhYhP5u33d//iQlw1fT0n43fHXYsPWYTAcRKSQ19SP8/YNxOp306tWLuNhYTIrCkyYTD2o22xUrVGDhwoX6qqFevXps3bqV/v37YzKZiNEqu+3du1e3s2hQowa/ZmRwSoRHtZrTycnJLF++HIBff/1Vz4to0qSJGr4K5Obm0q9PH8yKQoQIn9hseLQchPz8fJ5//nlsNhtBQUFMnjz5fGLbpk2cyMjgYW3VEBcezoIFCy69+b17z5vtpabCsmV/fqwN/Gl4Q1P4QkSiSnyOEpElV9K3NJtBCgZKE1ctQsOF4atPPPH/JnTt2pVNYuLXiICv7xbS0tQazk2bNqVdu3aICFUDAnhXhBizGbPJxMCBA5k4cSKBgYH4+vry9ttvs2bNGn3V8OCDD5Kdnc3kyZN1PeLNhx/GXaECK0RI1kqLPvbYY+Tk5OB2uxk7diz+/v44HA7eeustXRfYuHEjNTQBvI0IByMj4YMPwO3mp59+ol69euqWVbNmOqHg8cDHH/NNeDiVNHK467bbyMrKunQAFi6EpCR1SurYEY4e/fNjbeCK4Q1S+Pmiz6aLv/s7mkEKBkobVy1CA2Rnq1XMRFSvoOJ9+Mugb99vUZRjiOSSlvYePj6+BAYG0rdvXyIiIrBaLAwMCaGDNslWr1CBL7/8kltuuUV/yt+2bZu+aoiLi2Pp0qXs37+fli1bqhFEDRqwc8gQzoaG0lsERYSEuDiWaU/pe/fu1c9Xv359fv31VwAKCwsZMWIEdpuNQLOZySJ4ateGNWtwu928++67OqG8/vrr53MXzpzhXL9+vGQ24yNCoK8vE8aOvdAuAyA3F559Vq36FhLyu4K9gauDN0jhHRFZIiKdtbZIRN6+kr6l2QxSMHCtUJwJ7RruYvqW6X/tJKtWnQ9fbd/+/w1fzcw8THj494iAy7WWSpVaqE/obdpw++23q7pDQgLv2u2EiWAzmxn5wguMGzcOPz8//P39mTRpEt999x2pmsBcvBqYNGkSAQEBOJ1Oxrz2Gp6nn+Zri4WyWhLb/x55hDNnzuDxeJg6daq+CnnjjTf0VcO2bdt0Y79bbDb2iUCnTnDwIHv37tW3terWraua7BVj+3a2NWyou8c2qVHj/KqiJLZuvVCwN+pFex3eij66Q0Re11rbK+33/5wrUERmicgvIvKziNQVkWARWSoi27XXoD86j0EKBq4lrlqEBjV8ddAgNXw1JASmTr3s07Db7aF9+5WI5KAo2aSnv4HFYiE6OppnnnmGwMBAHA4Hr9SpQxttkm1Qtiwrly3To49at27Nb7/9xhNPPIGiKKSkpPDNN9+wd+9ePfehefPm7P36a862aUNP7TzJ4eGs/vprAA4cOKBP8g0aNGCHFlHldrt58803cTgcBNhsTDKb8Tid8PLLePLy+OCDDwgODsbHx4eXX35Zt+3G48E9Zw5jg4LwF8FhsfD60KGXhq9eLNiPGgUXZ00b+MvwCil4s4nIVBHppr23aSQxUkT6ad/1E5ERf3QegxQMXGt4RYQG9Wm4bl31X7B58//3afjLL3fj57cZEYiKWkZy8g1qNFDnztyoWWPcWr8+byQmEiCq0d17/foxatQofHx8CAsLY+7cuaxYsYKEhARMJhP9+/cnPz+fsWPH4nQ6cblcah2GFStYUbYsCdqW0tMPPEBeXh4ej4cpU6bgcrnUFcaYMfrWz44dO2jYsCEiwm0RERwUgbJlYeFCsrKy9CzqOnXqXLhqOHOGfT16cKu2QqmfksKvv/xy6QDs2we33aaOU3q6GvZr4Krhje2jdBFZJyJnRKRARNwicvpK+l7mXC4R+U0uqv0sItuKxWxNyN72R+cySMHA3wWviNBuN4weff5peOTIy4av5uUV0rTpckQKMZv307DhID3j+amnnsJutxMaGsr4Bx+kkc2mbjXFxbFqzhyqV6+uh6ru37+fBx98EBGhWrVqbNmyhR07dugi8Z133snRw4c5PXo0D2lW35UDA9m4ZAnABSuMm266if1aeKrb7eb111/HbrcT5OfHjKgodWpp3RrPzp1Mnz5dXzWMGjXqglWBZ8sWppYrR6AIviYTb/Tvf+mqweNRazeEhKgurMOH/79hvgauDN4ghfUikiIiG0Utz9lFRIZfSd/LnCtNRNaKyBTtfO+JiFNETpY4Rin5+aL+D2nXsz4+Pr6Uh86Agf8fXhGhQX0avv129d+xenVYv/6yh02Z8qNmk+GmfPl5REaWwWq10qdPHz0J7cGOHXmpUSNsIkQoCvMeeIBnn34ak8lEYmIiX3/9NQsWLCAsLAwfHx9ee+01PVfBarUSGRnJwoUL4fRpFt59N5EiWEUY3qIFRZrW8O677+JwOAgMDOSjj84n+f3yyy96Hep7q1XjuMOhWo0//zxZu3bRunVrVehu2PBCLcHt5sCrr3Kr1YqI0KhMGXaVXFUU4/Dh8+GrtWurtS4M/CV4hRS0180lvtt4JX0vc64bRKRIROpon98UkaEXk4CIZP/RuYyVgoHrAV4RoT0emDULIiPV8NUnn7xs+Orhw2f0hDe7/ScyMh7UJ9oePXro2sH7o0ZRWSuO85jLxVfDhpGYmIjJZGLAgAHs37+fVq1a6WGv+/fvZ+PGjVSqVEkVnP/3P3Jzczm6di13RkermoLdzm/jxoHHw6+//qr7Id13331kayZ6hYWFvPjii6r+ERHBF40aqdNMUhKezz5j8uTJBAQE4Ofnx8SJEy+IQPIcOsTE9HT8tW2wif36XRqh5PHAjBmqnYjdDm+8YVhz/wV4gxRWaXv/07S9/z6i2Wj/2SYikSKyu8TnBqImxxnbRwb+sSgpQnec2/GvidCghq8+9JD6r5mYqFpCXAYDB36PohxFJJfatafgcKjawODBg4mPj8dsNvPcc8/RS4tWqiDC6owMutx5JyJqtbZt27Yxfvx4HA4HQUFBzJw5k7y8PPr0UV1cK1SowA8//KBGIj39NP4mE/4iTKtUCc+PP15AAHFxcXoiHKgmesWlRXu2bUtu2bLqPbVpw+4SBYTatGlziYHe7qlTaazVe7g9MZEju3ZdOgBZWdCqlXrOG2+EPXv+2nj/R+ENUigjIr4iEiAig0XkNRFJuZK+/8/5vhaR8tr750XkFa2VFJpH/tF5DFIwcD2hpAid9GbSXxehAVauhHLl1H/RDh0um8yVmXmY0NC1iEBIyLdUrar6FbVr1467tUI89erVY9qkSUQFBGAV4VWzmY/btCEoKAin08nEiRPZtm2bXt+5a9eu5OTksHTpUqKjo7Farbzyyiu43W52bdtG/aQkdXWgKGT36AGnTrF27VrKlSuHoig888wznNPsKnJzc+nZs6eqTVSqxOZevVTtxOHAPWIEr44Ygc1mIzIy8ryHkgb3qVO81qABPiJEmEwsfPHFS8fI44GJE8HPD1wuVXcwcEW4alLwdtN0hfUisllE5olIkKjW3F+JGpL6pYgE/9F5DFIwcD3CKyI0QF4eDBx4Pplr2rRLwlcLC920bbsCkXxMpsM0aTIck8lEUlISgwcPxt/fn4CAAMaPH08bzfKiuQjrw8K4Udsquuuuuzh8+DDPPvusvv20bt06jh07Rtu2bREtU/nAgQMUFRXxYv/+mBWFeBFWBwXB1KmcycnhoYceQkR1fN22bZt+jYsWLSIiIgKbzcYbgwbhKY4mqlKFTdOm6VtWvXr1Iv8iS+3N06ZRWRPPH69YkbxDhy4dp507z+c1tG8Pp/7iKu0/hL9MCiKyRZu4L9t+r++1aAYpGLhe4TURGmDLFjUcUwRuukm1oL4IM2duw8fnV0SgfPlPiYtLxmw28+STT+p7/507d+bNN9/E18eHULOZ+SK8nJSkb/+sWrWKFStWEBcXh9VqZeTIkRQVFTF+/Hh8fX0JDQ3l008/BWDNmjUkxcRgEuEFEYoyMmDTJubMmUNwsOrdVNIL6ciRI3q+Q8uWLTk8eTLExoKikNu9O//rrpoBpqWlXRi6CuRlZ9P7BjUUt6rVyo8TJlw6RoWFMGQImM2qXcaaNX99vP8DuBpSKKO1kVqrorURIvLy7/W9Fs0gBQPXO7wiQoOauPX22+pWicMBr756SVjm8eO5VKmyQhehGzbspkbzNGpEr169MJlMpKSkMHPmTKpVq6bu9/v48LWikOxyYTKZGDx4MIcPH9bzC5o3b05WVhY//fTT+T49e5Kfn8+pU6fo0L69KkJbLOxTFHj8cfZt3aon0N13332c0p7cPR4Pb7/9Nj4+PkRGRvLlggXQq5cqrEdF8Wn//oSGhuJwOJg0adIlIvPC114jzGzGV4T3GjfGc7kyqN98o1bCs1jUEF9DhL4svKEpXBJpJFoVtr+zGaRg4J+AkjWhO87tyOn803/9ZHv3nhdYa9aEH3645JBnn12DohxHJIcGDSbi6+sgODiYYcOGERsbi9VqZfjw4Tz++OPq03lICOtFeEAzyGvQoAF79uxh3Lhx2O12wsPDWbx4Mfn5+fTq1QsRoXr16voW0bRp03A6nQT7+DBfUSA8nKLJkxk6ZAgmk4nk5GTWlwizzczMJDU1FUVRGDhwIIXffac6yYpwoEULbtRKlLZv357Tpy8cq6xdu2gWH68Sjp8fp7744tIxys6Gu+5Sx+iWW363Gt5/Fd4ghU0iUq/E5wwxajQbMHDFKHQXMmjZIO+I0B4PfPIJRESo2yVPPw1nz15wyLp1BwkM/AERiI5eQaVKdRERunfvrucL3HLLLbz//vuEhITgsNuZkpjINBH8zGaCXS7mz5/P1q1bqay5oz7zzDMUFBTw6aefEhISgtPp5IMPPgBUP6TiRLnekZGcE4GGDVk1ZQqxsbHYbDbefvtt/en/zJkzdOnSRQ+nPbB7t1p7wm6nyN+fIa1bYzKZKFeuHJmZmRfcm9vt5qVu3TCLUFaETd2766VDLxijd99VS4DGxMDq1X99vP+F8AYp1BSRTBHZLSJ7NJKocSV9S7MZpGDgn4av93yti9BDVw796yI0wIkT8OCDFOcBsHTpBT+fO1dEs2bLESnEYtnDTTcN1LOZBw8erG/jfPzxx3qIaIf0dH5wuaih+SD1euQRTp48ycMPP4yIanK3e/du9u/fT4MGDRAtW/rs2bPk5+frq49aCQnscrnAYuHY//7HbZrr6h133KHnNIC6ynA4HISFhbFkyRLYvh0aNwYRVqSlERkWht1u57333rvk9lctWkS0w4FdhPcSEkBzc70AP/wAycnqdtKoUYbrqgavRR+JalHh+jN9SrMZpGDgn4jsvGzum3Uf8rzQYFKDqxOhAZYvV/2GRKBzZzh27IKfx47djNm8D5ECatX6mKAg9Sl/6NCh+jbOgAEDGDRokPp0npzM2jvuoJdGDDckJrJzxw5mzJiBv78/QUFBLFiwgMLCQgYMGICiKFSqVImftAzj2bNn43K5CHS5mNe0KYjgiYvj1S5dsFgsJCUl8UOJba+ffvqJypUroygKgwcPpqigAMaPh4AADvn60lSr8Na1a1dyc3MvuLfDhw/TrGpV9XeLhdyxYy+d+E+ehLZt1fG55x64nBbxH8PVCM0dtNcnLtd+r++1aAYpGPgnw2siNKi1CPr3V5+Iw8LUmP0Sk+OuXdlERX2HCAQHr+GGG9Qn9/vvv5/2mljcsGFDZs2aRXR0ND4+Pozp1485KSkEihBgNjPrtdfYvn27vk305JNPUlBQwJIlSwgLC8PhcOjbSTt37qRmzZrqcffeS0GlSiDCNw0aEBMZiY+PD+PHj9e3k86ePUvHjh0RUX2Vjh49quont9xCkQgD4+L06KSL7baLiooYqOVDVBfht1tvVYmgJDwetfynyQQVK15+VfEfwtWQwsPa6+DLtd/rey2aQQoG/unwqggNkJkJtWqp/9otW16Q6et2e2jXbgUieZhMB2nRYhiKopCamspLL72E0+kkNDSUGTNm6GU8291zD5sHDaK2Vgq0Z/XqnNy3jx49eujbSXv37uXAgQP6dtIjjzxCXl4e+fn5+nEN6tfnwIAB4HBwxOHgpvLl9TDZs5oe4vF4GD9+PDabjfj4eNatW6dO5lOmgMvF5zYbQVr29cKFCy+59U/nzcNltxMswtLISPj++0vH58sv1ZwPlwsWLbq6sf4H47pLXvNWM0jBwL8BXhWhQQ1fff11cDrV9sYbF9QimDHjF81Yr4i0tNmEh0dht9t54YUXLhCVX3zxRcxmMykpKaxdsoTe2m+1rFZ+e+cdZkyfjp+fHyEhISxZsoTCwkKefvppPYHtt99+A+DDDz/E4XAQERHBiunT9af/QZGRKIpCtWrVLnj6X7t2LfHx8fj4+DBx4kT1y/374ZZb2ClCNT8/FEXhxRdfvCRs9ddff6VSYiImEUaZTHgupyPs2gXVqoGiwCuv/Cd1hqtZKbz1e+33+l6LZpCCgX8TSorQV5UJXYzfflNDMkVUV9HNm/WfDh7M0WtCu1zrSU9vo1tkdOrUSQ9NnT17tr6dNH78eGYNG0aAyUSQCJ9Vr862JUt0PeD555/H7XYzf/58XC7XBU/0W7dupXz58pjNZl595RU8M2ZARAQLTSYCfXwICgpiUYkn96NHj9KsWTN95XHu3Dl18p4wgbNOJ+01R9U77riDnIs0gpycHO7UzP46ipDXqpUaploSZ87A3Xej6zAXZVL/23E1pNBJa+NFZLWIPK61VSIy9vf6XotmkIKBfxtKitBXnQkN6kT64YcQGqrqDQMGqPYZqNtJnTuvQuQsJtNhWrQYgclkonz58gwfPhyn00lYWBiffPKJXkehY8eObP7hB9JiYhARBpjNnO7Xj46aLnHzzTdz/PhxduzYQbVq1S4gi1OnTulJce3atePMvn3QvTs7RKhms6EoCsOGDdOf/kuuPOrXr8+hYouLXbvw1K/PKBFMIlSpWFFflRTD7XYz5IUXEBHSRTgUHw8bN3LRQTB4sDoNNmhwiUD/b4Y3QlLXiIilxGeriKy5kr6l2QxSMPBvhMfjuUCEnrFlxtWf9OhR6NhR/ZcvV0413NMwZ86v2Gw7EHGTljZX304aOnQoFStWRFEUBg0axKBBg1AUhcqVK5OZmcmD992HiNBUhMMJCYzt2ROr1UqZMmXYsGHDBeLxrbfeSnZ2Nh6PhxEjVPKpXLky27dvhxUrOJuczH1atNNdrVtf8PQ/ffp0fH19iY2NPZ8EV1QEI0eyxGwmUFEIdblYtWrVJbc9a9YsfH18iDObybTZVH3iYnz0kVq8JyXlPyNAe4MUtkkJgzpRDez+0Nq6tJtBCgb+zdh5Yifp76Ujzwud5na6ehEa4IsvVEtuEdWiW9tWyco6v50UFLSWOnXU7ZcOHTpwnzb5t2jRgk8++YSQkBD8/f2ZPXs2kyZNwm6zEWux8L0I3990E7HadtOUKVPweDyMHj0ai8VCSkoKW7du1S7jC4KDgwkMDFTzE3Jz8TzzDK8oCiYRqpYpc8HT/w8//EBcXBy+vr588skn5+9n40Z+TU6mvAhWk4nJ48dfcssbNmwgJioKP7OZhSLw2GOgObnqWL1aFaBDQ/8TvkneIIUuWtLaFFHrK/8mIp2upG9pNoMUDPzb4XURGtT99L59dc8hZs8G1O2kjh1X6dFJN9/8IoqiUKVKFYYMGYLNZiMuLo758+dTu3ZtRIS+ffuydu1aEhMSsJpMjDWbOeznRxMtuqhHjx6cO3eO1atXExkZidPpZNasWYAatlqlShVMJhOvvPKKum20YQOLk5JwiRBis7Fi7lz9sg8dOkSGZoHx/PPPnxeZc3M50bUrzbSVxjPdu19S0nP//v2kVauGSVEYIwL166uV3Eri11/VRDdfX9CM//6t8Er0kajFcW7XWuSV9ivNZpCCgf8KvJoJXYz163XPIdq2hQMHgOLopN8QKaR27Y8JDAzC39+fESNGkJCQgNVq5c033+TRRx9FRGjcuDE///yzHsbaNSqKHBGeiopCRK3pkJWVxYEDB3TH1oEDB+J2uzlz5gx33XWXvirJy8uDggK29e5NeREsIkx46CE9Qig/P18Xwu+77z71eA0Fc+bwiFac5+5atci9yPojJyeHW2+9FRGhv8WCJzb2Uu+ow4fhhhtU+5BJk65+jK9TeIsUYkT1PGpY3K60b2k1gxQM/JdwcSb0npNeqDZWUKB7DhEQAOPGgdvNnj0niYn5FhEID19N1aoN9Sf/4sm/Q4cOTJgwAV9fX6Kjo1m9ejXPPfecngW9NziYGYqCw2olOiqK77//nvz8fB58UC0h2rp1a06dOoXH42HIkCGICLVr1+bgwYPq/X7zDTcFBCAiPFG2LEVZWYCquQwfPlzPkyhZuc2zbx+vJiejiFA3LIyjFwnQhYWFet2HBxwOCux2tQxqSZw+Dc2bq9PjK69c/Rhfh/DG9tEIUX2PPheRT7W24Er6lmYzSMHAfw0ej4dpm6Z5V4QG1XPoxhvVKaFhQ9i2DbfbQ5s2KxApwGLZTfPmT+mhqk899RSKolC1alU+/fRTEhMTsVqtjB07lvnz5+Pv709YaCgrWrYkU4QEiwUfq1XXGd5++23MZjMVK1Zkx44dAMyZMweHw0FsbKxug1GYl8fjdVUzv9tsNnK0jGlQRWS73U5SUhK//PLL+XspKmJWu3b4iFDWamXnRYluHo+HoUOHqtFSLhdnROClly7MVzh3TrXEEFEjtv5luQzeEpp9ruTYa9kMUjDwX0WpiNDF5S0DA1V30ZdegoICxozZjMl0EJFcGjWajMPhICoqildffZWgoCACAwMvyILu3r07mZmZlC9fHovFwju9e3M0JYWm2p5/n4cfprCwkK+++org4GCCg4P1+s4bN24kPj4eh8PB3BJ6wjvPPotJhDQR9t9xhy6Qr1mzhrCwMIKDg1l9kRPq6rfeIlhRiBDhhyFDLrndCRMmYDKZqBsayonifIWSAnRREXTrpk6Tffr8q4jBG6SwSET8ruTYa9kMUjDwX8bFIvSafV6KmsnKOp/YVbUqrF3L1q1HdCvuhIQvSEhIxWKxMGjQID0fYfDgwfTv31/NDUhP5+eff9arrT3YuTNn+venp8mk+htVqsQJLZ+hQoUKWCwWxo0bp/35LGrXro2iKIwcOVIXlBcuWICfzUasCJnh4bor7M6dOylXrhw+Pj66iF2Mn1auJN7HB38RvmrZUs/RKMbs2bOx2WxUCQ8nSwSaNLnQN8njgZ491bF47LF/DTF4gxRmi8gOERlnZDQbMHB9oVREaIB58yA6Wo1S6tOHvGMnSU9fjgg4HFupV+9+PRHt3nvvRUS4/fbb9aI7kZGRrF69mgEDBiAiZGRkcGjVKiakpGAVoZzDwbYvvuDkyZP6KuOJJ56gqKiI3Nxc7rnnHkSEhx56iAKtXsKmTZuICQ/H32RiqYhauS03l2PHjpGRkYGiKLzzzjsX3Mb+3bupHBaGTYTZyckX+EEBfPnllzidTspGRLDXbIYqVVRbjWJ4PPDUU+p0+eij/wpi8AYpdLpcu5K+pdkMUjBgQIXX7biLcfIkPPII2hIBvviC/v3XIHISRTnGTTeN1LWF5557DrPZTGpqKgsWLCApKQmr1crEiRP55JNP9AS0DevWserJJwlVFAJF+LJ7dwrz8uipOZ22atWKnJwc3G63vvJo0aKFXtZz3759VKlUCYvJxDQR1fV040Zyc3P1AkIDBw68wBfpxIkT1C1fHpMIk/z8YNmyC27zm2++ISCljkcmAAAgAElEQVQggDIREex0OiEuDjQrcEAlgqef/tesGLwSfXQ9NoMUDBg4j1IToQFWrVIzoUWgUycWz8zEx+dXRIqoWXMmLlcgQUFBjBgxgtDQUAICAvjoo490e4yePXuydu3aCxLQdn37LZX8/bGIMDYuDjZuZPTo0ZhMJqpXr85+7Wn9vffew2w2U61aNf27kydP0qRJE0SEEf7+eKxWePVVCs+d06ObHn74YYpKGAGeOXOGm+rXR0R4S1FUo8ASk/v69esJCgoiLjKS7SEhajJbSadVjweefFIdg759/9HE4I2Vwm8isuvidiV9S7MZpGDAwKUoFREa1P34Z59VPZTCw8ka9yFxcd9oYatfk5paC5PJxDPPPENaWhqKojB06FB69+6t2mE0bcrPP/9MXS2aaMiQIZzMzuYWrT5DH0WhqH9/Fs6fj5+fHzExMXopzsWLF+Pn50dsbKyeFZ2fn69vW/VKSsItAs2b4zl4kH79+ulbW+dKiMf5+fm01XSO4SLQpcsFZngbN24kJCSEmIgIdsTFqS6zX311fgw8HnWlIAIvvOCdcf0b4A1SCCnRYkSkt4gMuZK+pdkMUjBg4PIoKCrguWXPYXrBRPKbyd4ToQE2bYKaNUEE922taN18ESKF2GzbadCgGyLCnXfeqWsCd999N+PGjcNms5GcnMzGjRt1T6T77ruP06dP01PLHWglQk5qKpvef5+YmBj8/f1VGwzUCTsyMpLAwEBWat5NbrebXr16ISLcX6uWmncQFgYLF/LKK68gotaiPlsika2wsJD7NeuO50WgXr0LspszMzMJCQkhLjqaXeXKqZFYCxacv3+3W41UEoG33vLeuF5DlMr2kYhs+Kt9vdUMUjBg4Pexavcq4l+P974IXVgIr76qWkIEBDDqrtEoylFETtKkyeuYTCaqVq1K//79MZlMpKWlMWfOHCIiIggICODzzz/n5ZdfRlEU0tPTOXToEO+88w4mk4kaVisHFYV9Dz1E1SpVMJvNTNKyi3fv3k1qaio+Pj7MmTMHULfNhg0bhojQsmFDzmpV3ujbl3GjR6MoCo0bN77AZK+oqIjOnTsjIgyyWPCUKQNbtui/b9y4keDgYBLLlGFftWrq6qhkZFNhIbRpo9ZkmOHFbbprBG+sFGqUaDeIyCMiknklfUuzGaRgwMAfo1QyoYuxY4caxinCt2mt8bVv1SJZZxMQEEhISAjDhg0jICCAsLAwZs2aRVpaGiaTiddff53Zs2fj6+tLfHw8W7Zs4fPPP8fpdBLndLJFhFPJyTSvVQsR4YUXXsDj8XDs2DHS09MxmUyMHTtWv5Rx48ahKAoN6tXjZNeu6vRWuzYfvP46ZrOZunXrcrJEuKnb7aZr166ICM85neDvD9qqBGDdunUEBARQLiWFQ8XWF9NLlE3Ny1MT/azWC7eY/gHwBiksL9GWilpfofyV9P1/ztdHRH4Uka0iMl1E7CKSKCLfa6GvH4uI7Y/OY5CCAQNXhlIVoT0eeO89cLk4bnORFPwpIhAW9jUpKTUwm80MHDiQlJQUrFYrY8aMoW3btroYvGbNGqKioggICGDJkiVs3LiR6OhoAhwOvgoPp0CEThUrXiAenzlzhpYtWyIiF1RgmzFjBhaLhZo1a3Js0iTVuiMwkNn9+mG1WqlVqxbZJQrulCSGIRER6sRfXO0NWL16Nb6+vqRVrcrJjAw1PLfkyiA7GypVUst7alrHPwF/mRREpJf2Wv/3jvszTdMkfhMRX+3zJyLSWXu9V/turIg8+kfnMkjBgIE/h1IToUE11GvbFrcIbV2DdJ2hTp0OqtfQAw/oEUN9+vThmWee0QXozZs3U7VqVcxmM+PGjWPv3r1UrlwZi8XCtCZN8IjQPzgYEaFNmzbk5uZSUFBAhw7quXv37q07o3722WfY7XYqVapE1nffQY0aIMKCtm2x2WzUrFmTEydO6Jftdrt1jWNk2bLqtPjcc3p00eLFi7FarTTIyCA3I0MljpJbSbt3Q2QklCkDxYWArnNcDSls0l5/+L3j/kzTSGGfiASLiEVEPhORFiJyTLRCPiJSV0SW/NG5DFIwYODPo9BdqIvQXs2ELsasWRARwatKUxQ5iqJk06iRamBXr149/cm8ZcuWjBkzBqvVSmpqKps2beKWW25BROjXrx8nTpzQSWRY16544uJ4S1FQRGhYvz7Z2dm43W49v6FLly4UFhYCsGzZMpxOJ+XKlWPfjh1q0pkIn1eqdFliKCoqol27dogIo9PT1anxwQdV7QB1BaIoCm1btaKobl11y+izz87f87p1qr6SkfGPKO15NaQwXUS2i8hZEdlcom0Rkc2/1/cPzttLRM6IyFER+VBEQkVkR4nf40Rk6x+dxyAFAwb+OkqK0C+ufNF7IjTA8ePQqROrJAG7kolIEdWrf4jN5kN8fDwDBgzAbDZTqVIlpk+fTnBwMKGhoaxcuZJHHnkEEeHee+/l9OnT3H+/mjn9WPfuFHXpwnQRrIpC1XLlyMrKwuPxMHjwYESEu+66Sw9BXb16Nf7+/iQlJbFnzx744ANwOPg8KAibtpVUUmMoKCiglVbb+f3WrdXp8fbbdWuMN998ExGhR7dueGrUUKOSVqw4f88zZ6p9One+7nMYrkpTELWOQqaIlLm4/VHf/+d8QSKyTETCRC3rOU9EOlwpKYjIQyKyXkTWx8fHl/LQGTDw70Z2Xjb3zroXeV5oOLmhd0VogMWLORxdllj5BBGIjlpKVFQiTqeTF154gcDAQMLCwpgxYwYpKSn4+Pjw0UcfMWLECESEhg0bcuzYMfr27auHuubNmcOSwEAcIiQHB7NTK6E5atQofQWSm5sLwPfff4/L5SIhIUGt5LZlC6SksMBkwmIykZGRcUFUUl5eHjfeeCNms5lPH35YnSIbNQItm/qpp1Sn2JGDB0OFCqo4vWHD+fsdNIh/QqjqVQvN6jnEJiKVtWa90n6XOc/dIjKxxOeOIjLG2D4yYODvQUkROvDlQD7e+rF3/8Dp07gf7UETGYgIOHw2UqlScxRF4YknniA5ORkfHx/GjRtHgwYN1MSy4cP58MMPsVqtVKxYkT179uiTfuPGjTm1axff3XgjwSJE22z8qEUNFUcgNW3alDNnzgBqpnJgYCAJCQns3r1bFYdvu41ZIpgVhRsbNdJJRL3c09SqVQu73c43gwer4ai1asHx47jdbj3/Yvb48aqOEB6uRmGBmsPQqpXa5yLX1usJV00KItJI1HKcK0VklSYU/6UiOyJSR9TII4eIKKKW93xcRGZeJDT3+KNzGaRgwID3UKoiNMCqVTwT1AWRHExygDo11W2ie+65h/qa/cSgQYP0LOXu3bvzxRdf4HK5iI6OJjMzkw8++ACLxUJaWhqHsrLYMnw4kYpCiAjrNHF46tSpmEwmGjZsyOnT6j0UE0NiYiJ79+5VJ+9Bg3hfBEWEW5s21U33AI4cOULZsmUJCgrix3feUbeKqlSBQ4fIzc0lPT0dX19f1s+apdphJCdDcbGf7Gz1c3T0pSU/rxN4gxQ2SIkQVBEpJ1eRvCYiL4jIL6KGpL4vIj4ikiQia0UNSZ0pV1C/wSAFAwa8i1LNhAbIzWVm64GY5TdEcqlTbhAiagW14ifwBx54QI9Muvnmm1mzZg0xMTEEBASwfPlyFi1ahMPhICUlhd9++40dq1aRYLfjL8LXjRrBsWPMmDEDs9lMvXr1dCO9devW4XK5SElJ4YBWepTZsxljs6kZ0S1aXFDbedeuXURERFCmTBkOfvwxOBzqltHBgxw6dIj4+HhiY2PJ+uwztXJd3bpQvOLYtEklkptuUgnoOoM3SOESUflqhGZvNYMUDBgoHZSqCA389PEyAkyrEYHUgNex2ewkJiby+OOPIyI0atSI1157DbPZTI0aNVi/fj0VKlTAZrMxc+ZMvv32W4KCgoiOjubHH39k3+7dlA8Lw1eEpSEh8NVXzJw5E4vFQnp6uk4M3333HX5+fqSmpnK4+Cl+0yaGu1yICP+7+eYL3FXXr1+P0+mkRo0anFm0SPVCKlcODh5k48aN+Pr6kpGRwbnp09Xs5rvvPk8CY8eqU+zLL3t17LwBb5DCZBF5T0Qaa22CiEy6kr6l2QxSMGCg9JCdl027me1KTYTOOZ5D2eC5aqKbMouwgEhcLhf9+vXDZrORmprKxIkTcTgcJCQksGbNGr1mwpgxY9iyZQtRUVGEhISwdu1aDh8+TJWUFHwUhc9FoF8/Zn/8MRaLhbp16+rEsGrVKnx9falatSrHjx8HwHPoEE9GRSEiDG3W7ILooc8++wyTyUTr1q0pWrFCJYbUVMjKYsaMGWpEUo8eMHKkOqUOHqx29HjgrrtUfWH9eq+O3dXCG6TgIyJPiMgcrfW5ku2d0m4GKRgwULrweDxM3TQVv2F+pSJCu90eWjVagIgbX1lHsrMcZrOZvn37EhQURHh4OFOmTNFLbi5btkyv5vbCCy+wfft2EhMT8ff3Z+XKlRw7dowaaWlYTSbmi0CdOswaPRqLxUK9evX0SKMlS5Zgs9moU6eO/p377FkeSEhARJhYvz6U0BjeeustRISnnnoKVq5Ut5IqV4YSkVHTpk49b5RXnNx24gTExkL58lDClO/vxtWGpJpF5Jc/Ou7vaAYpGDBwbVBShO48r7PXRejnBnyHyGlMsp+q5tqICF27diUhIQGHw8G4ceNISUnB19eXuXPn0qlTJ0SExx9/nL1795KamordbmfRokVkZ2dTu3ZtLGYzsx0OCAhgZp8+mM1mGjVqpDumzpkzB5PJRLNmzcjXEs4Kzp3jpqQkzCIsrF5dD0UF6NGjByLClClTVK8jHx+oUYPCY8do1KgRvr6+bNmwAdLT1dVEscHel1+qU23Pnl4ds6uBN1YK80Uk/kqOvZbNIAUDBq4dSluEnjlzG2bzPkTOUF3aICLc2rw5NWvWxGQyMXLkSG644QbMZjMTJkzgiSeeQETo0KEDBw4cIC0tDavVyty5czl58iR169bFYrEwS7Ou+KhpUxRFoXnz5uRpCWlTpkzRk96KC/KcPn2atLg4/ETYVL68WrMaNbmtSZMm2Gw21qxZA59/rm4NNWhAliZKV6hQgTPbt6u2F2XLnieVxx9Xp9tVq7w6Zn8V3iCFVSKSIyJficiC4nYlfUuzGaRgwMC1R2mK0JmZh3E6t6gCtDyNiFArKUmv4Pbkk0/q71988UWGDh2KiFrC8+DBg9SpUwez2cz06dM5deoUGRkZmM1mZt1+O4gwKS5OP744BLU4/+GRRx7RReYDBw4QGxpKjKKwPzYWtAS5Y8eOkZCQQHR0NFlZWaprqqLAbbfx5ZIlKIpCly5d1MnfbIY771S1hTNnIDFRJYoSORF/F7yVp3BJu5K+pdkMUjBg4O9BaYrQx4/n6hXdoq1TsYmFJIeDe7T6yx06dOA+rUhO7969efvttxERbrzxRg4cOEDDhg0xmUxMmzaN06dPk5GRgcViYc6AARAczGgfH0TUAj/Fq4PiSm1DhgzRryMzMxM/h4MaZjNnQkPhhx/07319fWnQoIFKLKNHq1Npt24MePZZRISPP/5YrTchAu+8o57wq6/Uz88+67Wx+qv4y6QgqqV1bxF5R0QeFi3j+HppBikYMPD3oTRF6MJCN/XrL0cE/H1WESQughWFBzWDvObNm/Poo48iIrRv357JkydjNpupU6cO+/bto6m2VTRx4kROnTpFeno6FouFBRMnQno6L4sgIjzcrRsejwePx6M7pU4sYZ1dHHl0h68vbj8/ffvngw8+QER44okn1AMHDgQRCgYPpk6dOrhcLvbu2QO33QY2G2zcqB7XqZO65fQ322xfDSl8LCIfaIQwT0Te/L3jr3UzSMGAgb8fO47voM6EOqUiQnfr9jUi57BafiXWWg4fER5NSsJsNnPDDTfoT/i33norM2bMwGazUa1aNXbv3k2LFi0QEcaPH8/JkyepVasWNpuNxZ9+Cr160V8jhv7/+x+gagY33XQTZrOZxYsX69fw2muvqdnWISGqE6r222OPPYaIqBXgPB51whdhx6hROJ1OmjVrhvvwYYiKUpPezp6Fo0chKAgaN/5bTfOuhhS2lHhvES9aaHujGaRgwMD1gYtF6O/3f++1c48a9QOKko2iHCYl+GYUEbrZ7fjabKSkpDBkyBC14lqDBnoltwoVKrBz507dinvcuHGcOHGCtLQ07HY7y5cvxzN9Og9bLIgIo3r0AODUqVOkpaXh5+fHpk2bAHVF1KVLF9XvKCFBffJfsID8/Hxq1aqFy+Vi165dcO6cOtnbbIzVjPNGjx4NS5eqU61GPowZo37+G8t4Xg0p/PB7n//uZpCCAQPXF0pLhP78851YLHsQOUu5Mr0REdqJEGS1EhEezssvv4zFYqF69erMmzcPPz8/kpOT2bZtm16hbdy4cRw9epSKFSvidDr59ttvKdq8mbv8/NQ8gw4dwONh//79xMTEEBsbq9th5OfnU6dOHZxOJ1srV1a3gGbPZteuXbhcLmrXrq3qC8ePQ0oKntBQmjdogNPpVAmjVy91ul26FIqKoFo1iIv723IXroYU3CJyWms5IlJU4v3p3+t7LZpBCgYMXH8oLRF669Yj+PltRsRNUqJa2+BmEWJNJgIcDkaMGIGvry+pqanMnz9ft8z++eef9RXDe++9x8GDB0lJSSEwMJBNmzaRf/gwTUJDsYiwuHlzyMtj06ZN+Pn5UbNmTT2vYf/+/URERFA2OZmTtWqpxDBnDrNmzVK3ofr3Vy/055/B5WJPpUr4+fnRvHlzPGfPqglsZcrA6dNqHQYRGDbMK2PzZ3HV0UfXazNIwYCB6xMej4cpG6d4XYQ+fjyX6OjvEIGoqE8QMVHXbqecCD5mMy8NHoy/vz9lypRh3rx5BAcHExsby9atW7n55ptRFIUpU6awe/duYmNjCQ8PZ9u2bZzKzqZaRAROETZUrgwHD/Lpp5+iKAp33nmnbpS3atUqLBYLbW67DU+dOioxzJ9Pt27dUBSF5cuXqxf62WegKLxTqxYiwtSpU+Gbb9Tw1cceU49p3Vqtx1DsrnoNYZCCAQMG/haUhgh97lwRNWosRwQCA5dhNvtRMTiYNBFMIgzq2pXg4GCioqKYPXs2oaGhREdHs3nzZpo3b47JZOKjjz7il19+ISwsjPj4ePbu3cvBgwcpExZGhKLwW0QErF+v5zAMLvYzAl5//XVEhFeGDoXatcFm48y8eZQtW5a4uLjz1dyGDMEtQt3EREJDQ1WfpZ49VWL49lv46ScwmaA4gukawiAFAwYM/G0oKCpg4FcDvS5Ct227QvVM8t2IwxFPXFgYGZoNdt+6dYmMjCQ0NJRPPvmEsLAwoqKi2LhxI40bN1YT2mbNYsOGDQQEBJCamsqxY8f46aefCPT3J9Vi4YTdjmfmTDp37oyIMEvzM/J4PNx5552YzWa+XrgQ0tLAbmfN2LGYzWY6duyoXqDbDbfcQqbFgtlspnv37urWUWys6ptUUKB6Jfn4wP79XhmTK4VBCgYMGPjbURoi9FNPfYdILhbLTlyuGoQEB9MkMhIR4dHQUGIjIwkMDGT69OmEh4cTFRXFhg0byMjIwGq18tlnn7Fy5Up8fHyoXbs2OTk5rFixAqvVSuOAAPJFyH/xRerWrYvT6WTz5s2AGqWUnJxMbGwsR3/+WdULAgIY2K0bIsK8efPUCzx2DOLi6ONyoSgKa9euhblz1an39ddh1y51C6o4MukawSAFAwYMXBcoDRH63XczUZTjKMphwsJa4nA4aJGWhojQ0WQiKTQUPz8/pk2bphPDunXrqFmzJj4+Pnz11VfMmzcPk8lEixYtKCgo0JPTOiUk4BHhYPv2REdHk5SUxIkTJwDYsGEDNpuNVq1a4dmzB+LiOBceTrUKFYiMjNSPY/VqTplMRNrtpKen49FWEPj7w6FD0K2bGuZ68OBVj8WVwiAFAwYMXDcoDRF64cKdWCx7ETlDZGRXLBYLLbXs5ztFKOdw4PD1ZdKkSfpW0vfff0/lypVxOp189913TJgwARG18pvH42Hw4MGICMMaNwYRvq1TB6vVyi233KLbY7zxxhuICG+//baqEQQH80NcHGazmc6dO5+/wBdeYLKWLPfBBx/AL7+oK4QHH1TrO5tM8OSTVz0OVwqDFAwYMHDdwdsidGbmYXx9f0KkkMhILdNZy1G4RVGoYDJht9kYN24cYWFhxMbGsmbNGj08NTMzUzfY69+/Px6PR/dYmvPoo2A2MyY+HhHh+eefB1SCa9myJT4+PmzdulUVkO12no2ORkRYunSpenGFhbjr1qWm2Ux8TAy5ubnQp48qOm/eDPfdB35+an3nawCDFAwYMHBdoqQInfJWylWL0AcOnCY4eD0iEBamRgnddNNNiAiNnU4qa6Gro994g+DgYOLj41m9ejUxMTFERETw66+/8vDDD+vZyLm5udSpUweHw8HGt97C4+tLRz8/FEVh0aJFABw6dIjw8HCqVq2q1mWYPZs8Ecr6+ZGUlKTnObB9O8s1M76Xhw9XE90CA+Hmm1WzPRG1ets1gEEKBgwYuK5RLEJbhlh4adVLVyVC5+ScIyFBrf8cFPQhIiZuvPFGFEUhPTKSqiLYFIW3nn2WwMBAkpKSWLZsGSEhISQkJLBnzx5atWqFyWRi/vz5ZGVlERsbS3x8PIcXLeJscDBVLBZCAgPZu3cvoBrniQhPP/20ehGvvsoybbtowIAB5y/unXe4VYRAh0PVHIpLeK5cCTfeqGY5FxZezVBeEQxSMGDAwHWP7Lxs7pl5j1dE6MJCNzfcsBwR8PNbhIiPXluhRnIy1cxmbCKMuvtu/P39KVeuHIsXL8bPz49KlSqxd+9eatWqha+vL2vXrmX9+vXY7Xbq16/PuU2b+CUiAj8RMqpU0esyPPTQQyiKwtdff62a3XXtygMiWM1mfvnlF/XC3G4ya9RAEeHZnj3V2grR0VC//vmopOJSnqUIgxQMGDDwj4C3RejWrVVisNu/QySAmjVrYrVaqZyaSjWnE6sII6pWxeFwULlyZebMmYPNZiMjI4Ndu3aRkJBAREQEu3fv5qOPPlJDXR99FPbsYboW+vpMu3YA5OTkkJiYSHJyMmfOnIFz5zhUqxYuEZqnp+sFfNi2jXtNJpwWC0eOHDlfj2HJEoiPhyZNrnYY/xAGKRgwYOAfhe3Ht+sidJd5Xcg5l/OXz/Xoo6sRKcBq3YpIBJUrV8Zut5NavjxVIyOxifBSUBB2Hx9q1qzJlClTUBSF2267jczMTFwuF5UrV+bUqVM8pbmfTpw4EQ4d4qHgYESEJUOHArBixQpEhF69eql/PCuLNwMCEBHmfvCBfk0//+9/KMWEkp8PMTHQoAG8+KI6LWuV3koLBikYMGDgHwdvitBDh65D5Awm02+YTGUpW7YsDoeD5ORkqiQlYRNhsMmExWSifv36eqhp165dWbp0KRaLhZtvvpm8vDyaNWuGzWZj7dq1nN27l0p2OxEiHJo6FVDrLCiKwurVqwEoXLmSSiIkORzka7WhycvjXj8//BSFYwcPwhtvqNPxvHlqeOozz1z1+P0eDFIwYMDAPxbeEqHfe2+LluR2CIvlBhISEnA6nSQkJFCpQgV8TCae1fyTWjRtyrNaac3nnnuO8ePHIyI8/vjjHDt2jDJlyhAXF8fRo0fZ8s032BWFloqCZ8ECcnJyKFOmDKmpqeRpJPCFViXuldtu069ny+jRanhrixaqhXZIiGqSd9ttqs5Q5L3a1xfDIAUDBgz8o+GtTOhPP92B2XwAkZPYbM2Ijo7G39+fuLg4UlNTsVssPKlFDd3VrJnuezRu3Dj69OmDiPDuu++yfv16fHx8aN68OUVFRbw9YgQiwltmMyxcyKJFi3RCAcDjoWVEBAEiHP3iC/16WoWFEawonNm9GwYNUqfk4lXDkiXeGLrL4m8hBRGZJCJHRGRrie+CRWSpiGzXXoO07xUReUtEdojIZhGpcSV/wyAFAwb+O/CWCP3tt/ux2XYgkovdfjfh4eEEBAQQExND2bJl8fXx4TGnExGhS7163HzzzZhMJubNm0fLli0xm80sW7ZMz4AePHiwmsTWvDl2ReFHmw2WLaN9+/ZYrVZ++uknAH789ltMIvR2uSBH1UhWv/8+IsI7N94Ihw+rdhePPAIul1res5Twd5FCQxGpcREpjBSRftr7fiIyQnvfUkQWaeSQLiLfX8nfMEjBgIH/HkpmQv9VEfqXX47hdG5FpBC7vTtBQUG4XC6ioqJITk7Gz+mkixZd1LNiRWpUr47D4WD58uVUrFiR4OBgduzYQefOnVEUhcWLF5OVlUVoSAjV7XYKHA4OL1xIUFAQDRs21COPut16K1YRftMiljweD+nh4SSLULRjh0oETid06KB6IxVrEF7G37Z9JCIJF5HCNhGJ0t5Hicg27f04Ebnvcsf9XjNIwYCB/yYKigoY8NUAlOcVUt5KYe3+tX/6HAcOnCYoaIMWsvoM/v7+uFwuIiMjiY+PJyAggLvKl0dEGBgZSUJsLBEREaxYsYKgoCAqV67MoUOHqFy5MqGhoezfv585c+ao20ZBQRAczDjNP2mqJkLv27cPH7OZTiJqIR5gxrvvIiIsaNoU1q5Vp+UePdTX+fO9Om7FuJ5I4WSJ90rxZxH5TETql/jtKxG54f8550Misl5E1sfHx5fKgBkwYOCfgZW7V+oi9LBVw/60CJ2dnUdUlFrJzcdnOHa7Ly6Xi4iICGJiYggODubmGjUQEYb6+hLo50eFChWYO3cuJpOJtm3b8uOPP+J0OmnQoAGFhYV07NgRs9nMupAQ3DExpNeoQXh4ONmar9GTvXtjEuHnsDDIzqagoIAYPz+aF4ei1qgBVapAcDA88EBpDNv1SQra52z+JCmUbMZKwYABAyVF6EaTG7H35N4/1T8vr8DrM0kAABe6SURBVJCUlFUaMYzFZrPjcrkIDw8nPDyciIgI6lWvjkmEFxUFq9lMkyZNeOWVVxARhgwZwvuaNjBw4ECys7OJiYmhYlIS+f7+rE9MRFEUevbsCcCRI0dw+vrSTkS1zQaGPP00IsK2Nm1gzBh1am7ZEoKC1GI8Xsb1RArG9pEBAwa8jotF6E+2fvKn+hcWuvUSn1brh1gsdgICAggPDyc4OJiYmBiqVa6M3WRigBad1LVzZ9q3b4+iKHz++ed06dIFRVFYtmyZHn307P33g9XKIzExmM1m1UkV6N+/P4oIP2q+R1lZWVhMJvoWu6ba7WrNBRFYscLr43U9kcIrFwnNI7X3t14kNK+9kvMbpGDAgIGS2H58O7Un1P5LIrTb7aFJE5UYLJYFmEy++Pv7ExYWhr+/P2XKlCEpKYlAu52HNWIY2q8faWlpuFwuMjMzKV++PNHR0Rw9epTOnTtjNptZ//zzHBUh0GajWbNmeDwejh49itPp5D6nE1JT4dw57mzZkhAR8nv0gHvvVVcJNhv07ev1cfq7oo+mi0iWiBSKyH4ReVBEQrStoe0i8qWIBGvHKiIyWkR2isiWK9k6wiAFAwYMXAZXK0LfeecKRMBs/gqTyR9/f3/Cw8NxOByULVuW8PBwYoODaaUoKCKMGzyY/2vvzuOjKs8Fjv+emclMJgnZyGIQSAJhcQNBwQCKWCxuLRbFpYVWQKyKVq73ai+gn2t7bxGs2tpeF7SIQouieN3qUkXEFZRNtoAJCXtYNYgJSyDJc/84J+MQAxLMMIl5vp/PfDjnzDvvPOclOU/Oec9539TUVO3WrZvOnz9f/X6/Dh48WL/88kvNysrS7t2768EJE/QvbiL55z//qaqqd955p3o8Hl3rDplde3YxOxBQnTXLOTx37+7M59zIonamEOmXJQVjzJG8v+F9bfendsfVCX399R8qVKnHs0AhWePj4zUzM1P9fr+eeuqpzmiqubna0+fTONBHfvMbFREdPny4Pvjgg6GH3F566SUF9N6JE/Xgz36mXUA7t22rBw8e1G3btmlsbKyObtdONT5eqzZv1pMzMvRSUL33XqejuUcP5zC9ZUujto0lBWNMi1S2r+y4O6H/4z8WKFSqx/OZQrrGxcVpZmamer1e7datm/p8Pu3Ts6ee7PNpW9DfXnyx8yDaww/roEGDNBgM6po1a/TKK6/UQCCgRZ99pq+6M7c97M7cdvPNN6s/JkZLfT7VUaN03Lhx6gXd3rat6siRqnFxzmF6xoxGbRdLCsaYFuv7dEI7A+ntU5HVClkaDAY1MzNTAe3Ro4fijpMUFNF80Is7dtSYmBh9/fXXtXXr1nrWWWfp+vXrNTExUQcOHKjVhYU6wOvVdJ9P9+zYoSUlJerxeHRcr16qIlrgnlk8BKr33OMcohMTG/3pZksKxpgW73g7oR966DOFr1WkRCFbA4HAtxLD5T/9qQJ6DWj7+Hht3769Tp8+XcGZee1R9wG1GTNm6KIHHnAecHOPX0OHDtWkxET9OilJ9dJL9czu3bV3TIzqRRc5TzXn5KhmZzdqW1hSMMYYPf5O6GnTVqlImcJmhU7q9/tDieGMM85QQC+77DIFdCxojIheetFFOmLECPV4PPrRRx9pfn6+pqena1lZmV7TpYvGg26fOlUXLFiggP7v4MGqoPeNGaOArhNxRkxNTHQO1Zsa9vzF0VhSMMaYMMfTCT17dqGK7FTYpnC6+nw+zcjIUEA7d+6sHo9H+/btqx4RHePeafS7ceM0Oztb8/LydP78+erxePSmm27SolWr1At6q9+vummT5ufna17Hjlqdnq7r+/ZVQCeD6tChzmEaVJ99ttH235KCMcbUUbavLDQn9LF2Qr/+eol6PFsVvlDoqR6PR9PS0hTQ9u3ba3x8vObl5WlyfLxeAuoF/eukSQrorbfeqmPHjlUR0YULF+oN11yjMaAb+vTRZ/7xDwX0jdGjVUHP7tpVeyUkqHburO4Tdaq1s7k1AksKxhhTj+PphJ47d4N6vZsUdivkq4hoqjstZ1pammZkZGhycrJ2bdtWc0Db+nz661/+0hn07tVX9aSTTtJevXrphg0b1O/z6SjQygcf1KysLL1k0CDVtDSd5A7EtwmcxJCUpHrOOY2235YUjDHmKBraCf3JJ6UaE7NeoVzhfAU0OTlZfT6fxsXFaW5urno8Hh3Qvbv6QC+Nj9e8nBzNzs4OzeL2t7/9TX9z663qFdHiQEDvcZ91KLn9dv3cvfz0sIhq796qXq9qINBo4yAdLSl4MMaYFi4vNY+PRn7EXefdxdPLnqbH4z1YVLroiOXPOacNCxcG8fu3A28AP+arr74iGAyiqmzcuJHOnTvz3vLlDMrP5429e/nJnj1s2rSJpUuX0q9fPyZMmMCYW24hxu9nYk0NNyxdisfj4YnqarrExdElKYlXUlJg82aorobKSigoiHxjHClbNIeXnSkYYxpbQzqhV63aqbGxnyscUHDuPoqLi9OEhIRQB7SIaO9TTlE/6LBWrRSc6T1FRMeOHatjx45Vr8ejJaA/69lT09PTtfLGG/UOj0djvF79GlRFnAs7Tz7ZKPuInSkYY8yx6Z/dn+U3LeeKU65gwrsTGDhjIJv3bK637GmnpbNqVQZxccXAi8AQ9u3bx6FDh0hKSqKoqIg2bdpQuHUrySkpLKyooENMDPdNmsSIESN45JFHuOKKK/DFxDApI4MbS0rYtWsXr3TtymU1NRyqrmYuQGYmeDywfHnE99+SgjHG1JESTGHWlbN46vKnWLx1Md2mdGN2wex6y3bsmMKaNW1JSCgEngeuobKykr1795KYmMjWrVuprq4mISWFtap0rapi3YYN+EWIj49n8uTJjB49mqfLyui8Zw/tW7Vi6ptv0m/AAFqJ8EZamnNTKlhSMMaYaBERRpw5gmU3LaNz685c/cLVjHplFBUHK75Vtn37JAoLc0hMXAXMBIZTVVVFeXk5CQkJVFdXs27dOk477TTeUGUQ8MS0afxq2DDefPNN8vPzAfjzGWcwsqKCOXPmsO2qqxioytv796M7dkBNDaxYEfH9tqRgjDFHcayd0G3atGLt2k4kJy8HpgMjUVXKy8sBSEpKoqCggOzsbD4NBskE5s6cSYcOHZg0aRLDhw9n6tq1XBYXh6oyY+dOLoqPZ+Pevayt/ZLdu2HXrojuryUFY4z5DjHeGP7woz/w3oj3qKyqpO+0vkz6cBLVNdWHlcvIiKek5BRSU5cC04AbANi/fz8VFRWkpKSwZcsWDgJpGRms3rOH3sEgq1evJicnh8rKSl7q3ZsBwPSpU7lwyBAA5gQC4PM5X/L55xHdV0sKxhhzjI6lEzo1NUhJyemkpS0EngDGAFBdXc3u3bsJBoMEAgFW7dzJGSkpvFhQQM+cHB599FEGDx7Mo0uXcnVyMsWlpezs359s4N3YWKejGSwpGGNMU1K3E7r7lO68sPqFw8okJ8dSUtKdzMwFOJNK3oaIAFBRUUFFRQXp6emsKS/H7/WiGzawc+dO0tLS2LNnD2X9+xME/v7yy1zQqhXzysupOXjQqbywMKL7Z0nBGGMaKLwTOi81j6tmX/WtTujExADFxWeRlTUf+AuqtyMieDweqqqq2LVrF8FgkITWrfkM6OXz8ewzz9C3b18eW7KEnwaDPP/225x37rnsrqlhZW3FEe5stqRgjDHHKS81j49HfcyEcyfU2wmdkOCnuLgXJ5/8MfAgqneiqni9Xnw+H+Xl5WzfuZOs9HTWVlVxcP9+WiUkUFpaSts+fSirqkLT0gB43+t1KrUzBWOMabpivDFMHDiRedfNq7cTOi4uhuLic2jX7kPgPlTHUV1dTU1NDTExMQQCAbbt2kW5x0Me8M6cOXTs2JG5u3aRCrzzzju093r5sDYp7NjxzXMLEWBJwRhjGsH5OecfsRM6NtZHUVEf2rf/ALgXuJuamhoOHTrEoUOHiI2NJTYY5HPAr0oisHzlSvp17syr27bRp00bPjx4EAVnDKSysojthyUFY4xpJEfrhI6N9bF2bT9yct4D/gf4LwBqamo4cOAA+/btIz4+nnivl89KSkhOTOSLxET2AQnADqCk9os21z/sRmOwpGCMMY0ovBO6U+tOXDX7Kq5/5XoqDlbg93spLDyP3Nx5wO+B34U+p6rs3buXL6qriQeSq6pYsHgx6YEAG7ZuBWB+beGNGyMWvyUFY4yJgPAnoZ9a9lSoE9rv91JUdD65ue8C9+CcNTg8Hg8+n49KETbs24fX4yGrbVs+rq6mlQgLaguuWROxuC0pGGNMhNT3JPTkjyYjHqWoaAC5uXOBu4GJgHMpqaqqihoRYoFkVYq2bOEAkOPz8UltxRF8gM2SgjHGRFj4k9Dj547nwr9fyLa9pRQVXUCHDu8AE3A6oB01NTUcAL5U5UBlJYmBAIcOHWIFUAGwbl3EYo1YUhCRaSKyU0RWhW27X0Q+F5EVIvKSiCSHvTdeRIpFpFBELopUXMYYEw21ndDTBk9jUekiuk/pzstFL1JY+CM6dHgbGA9MOuwzXveMQUVYB9QASwBKSyMWZyTPFJ4GLq6zbQ5wuqp2A4pwWgERORW4FjjN/cyjIuKNYGzGGHPCiQgje4w8rBP6xtdvYMnKfHJz3wLGAZND5atVOQCUHziAO8gFnwJ88UXEYoxYUlDVD4CyOtveVtUqd/UToK27fDkwS1UrVXU9UAz0jlRsxhgTTXU7oXs/eTazPkglJ+dN4D+B+w4r7wF8QDywCGDfvojFFs0+hVHAm+7yyUD4jbdb3G3fIiK/FpHFIrJ4V4THFTfGmEgJ74Q+UHWAfk/15Ybpy8jOfg34LeGJoQaoAiqBhQBVVRF7qjkqSUFE7sLZx5kN/ayqPqGqZ6vq2enp6Y0fnDHGnEC1ndBDug7hrnkTyL7nAbJOfom6iQGcg+YmYCeAO3lPYzvhSUFERgA/AYaphlJdKdAurFhbd5sxxvzgpQRTeG7oc0wbPI0l2xaz75ZRpGY8Q32JAdzO5m3bIhLLCU0KInIxzl4OVtXwi2KvAteKSEBEcoFOuGdJxhjTEoR3QndO60TZTcOITX4S55D5x8PKfgqwdm09tXx/kbwl9VlgAdBFRLaIyPXAw0ArYI6ILBORKQCqWgA8D6wG/gXcoqrVR6jaGGN+sELDcfefwIHbRuOJfwy4E7g/VOYdgIKCiHy/LyK1Aqr683o2P3mU8hOpfazPGGNasNrhuAd1HMSwxOGU3lUNlXcACvzWuXwUocl27IlmY4xpos7POZ+VY1Zw5fT3IeZhnDOGP3IA2L58eUS+05KCMcY0YSnBFGZf/TxTF8SB55vEcH+ELh9ZUjDGmCZORLj+rFGs2fZjnK7ZO/lT2LDbjfpdGsFp3SJNRHYBxzuweBoQuWfFI8tijw6LPTqaa+xNOe5sVa33Qa9mnRS+DxFZrKpnRzuO42GxR4fFHh3NNfbmGrddPjLGGBNiScEYY0xIS04KT0Q7gO/BYo8Oiz06mmvszTLuFtunYIwx5tta8pmCMcaYOiwpGGOMCWmRSUFELnbngi4WkXHRjqcuEWknIvNEZLWIFIjIWHd7qojMEZG17r8p7nYRkb+6+7NCRHpGOX6viHwmIq+567ki8qkb33Mi4ne3B9z1Yvf9nCjHnSwiL7jziK8RkT7NqM1vd39WVonIsyIS21Tb/Qjztze4nUXkOrf8WhG5LoqxN3ju+SZ9DFLVFvUCvEAJ0AHwA8uBU6MdV50Ys4Ce7nIrnPmsT8UZP3ecu30ccJ+7fCnOLHYC5AOfRjn+fweeAV5z158HrnWXpwA3u8tjgCnu8rXAc1GOezow2l32A8nNoc1xZilcDwTD2ntEU213oD/QE1gVtq1B7QykAuvcf1Pc5ZQoxT4I8LnL94XFfqp7fAkAue5xx9vUj0FRD+CE7zD0Ad4KWx8PjI92XN8R8yvAj4FCIMvdlgUUusuPAz8PKx8qF4VY2wJzgR8Br7m/zF+E/dKE2h94C+jjLvvcchKluJPcA6vU2d4c2rx2OttUtx1fAy5qyu0O5NQ5sDaonYGfA4+HbT+s3ImMvc57Q4CZ7vJhx5badm/qx6CWePnomOeDbgrcU/seOPNqZKpq7XRL24FMd7kp7dNDOLOC1LjrrYGvVLXKXQ+PLRS3+/4et3w05AK7gKfcS19TRSSeZtDmqloKPIAzU+M2nHZcQvNo91oNbecm0/51HMvc8001dqCF9ik0FyKSAPwf8G+q+nX4e+r8idGk7icWkZ8AO1V1SbRjOQ4+nMsCj6lqD2AvzmWMkKbY5gDu9ffLcRJbGyAeuDiqQX0PTbWdv4t8j7nnm5KWmBSaxXzQIhKDkxBmquqL7uYdIpLlvp+FO383TWef+gGDRWQDMAvnEtJfgGQRqZ3QKTy2UNzu+0nAlycy4DBbgC2q+qm7/gJOkmjqbQ5wIbBeVXep6iHgRZz/i+bQ7rUa2s5Nqf0bOvd8k4q9rpaYFBYBndw7M/w4HW2vRjmmw4iI4MxSt0ZV/xT21qtA7V0W1+H0NdRu/5V7p0Y+sCfsVPyEUdXxqtpWVXNw2vVdVR0GzAOGHiHu2v0Z6paPyl+Iqrod2CwiXdxNA3Gmh23Sbe7aBOSLSJz7s1Mbe5Nv9zANbee3gEEikuKeKQ1yt51w0vC555v2MSjanRrReOHc0VCEcwfAXdGOp574zsU5fV4BLHNfl+Jc950LrMWZpjXVLS/AI+7+rATObgL7MIBv7j7qgPPLUAzMBgLu9lh3vdh9v0OUYz4TWOy2+8s4d7U0izYHfg98DqwC/o5zx0uTbHfgWZy+j0M4Z2jXH08741y/L3ZfI6MYezFOH0Ht7+qUsPJ3ubEXApeEbW+yxyAb5sIYY0xIS7x8ZIwx5ggsKRhjjAmxpGCMMSbEkoIxxpgQSwrGGGNCLCkYY4wJsaRgzA+IiFwiIhNFxH63zXGxHxzT5IlIaxFZ5r62i0hp2Lo/2vHVR5y5GcZE+DuCIvK+iHjDNvfHeWK2T52yJ4nILBEpEZElIvKGiHQWEb+IfBA2HIZp4SwpmCZPVb9U1TNV9UyceQH+XLuuqgejFZc79MKRfoeSceYtaMw66xoFvKiq1WHbqoFhwJrwOoGXgPdUtaOqnoUzXHOm235zgWsaGqv5YbKkYJo9ERkuIgvdM4fHxZn5LcedDetpESkSkZkicqGIfOzO1NXb/WxtuZnizLb2gojEfUe9hSIyA2dIiXYi8rL713eBiPzaDWsy0NH97P3u58Jn67pDRH4XFsOx1FnXML4ZIwgAVb1bVa9S1bKwzRcAh1R1Sli55ar6obv6sluXMZYUTPMmIqfg/JXbzz2TqP1LGSAPeBDo6r5+gTOu1B3AhLBqugCPquopwNfAmO+ot5Nb/jRV3QiMcv/6Phu4TURa4wy7XeKezdx5DLtyLHWG77cfZ8yiDcdQ9+k48yscySqg1zHUY1oAu45omruBwFnAIucqCUGcYZc/wBlOeiWAiBQAc1VVRWQlzuxZtTar6sfu8j+A24ADR6l3o6p+Evb520RkiLvcDucAv72B+3EsdYYPb50GfNXA76iXqlaLyEERaaWq5Y1Rp2m+LCmY5k6A6ao6/rCNzox1lWGbasLWazj8Z7/uqJD6HfXuDVsfgDOfQR9V3Sci7+GMQlpXFYefmdct09A69x/he+pTwDdDaB9JACcRmhbOLh+Z5m4uMFREMgBEJFVEshtYR3sRqb1b5xfARw2oNwnY7R68u+JMLg9QDrQKK7cDyHDvpArgTMhyJEeqM0RVdwNeETmWxPAuEAjvmxCRbiJynrvcGvhCnQl6TAtnScE0a6q6GrgbeFtEVgBzcCZ2b4hC4BYRWYMzh8JjDaj3X4DP/exk4BM3ri+Bj0VklYjc7x5w/xtn/oI5OHMfHEm9ddbjbZw+kqNSZ3z8IcCF7i2pBcAkvrnEdQHw+nfVY1oGm0/BtGju5aDXVPX0KIfSYCLSE7hdVX/5Pet5ERinqkWNE5lpzuxMwZhmSlWXAvPqPLzWIO5dTC9bQjC17EzBGGNMiJ0pGGOMCbGkYIwxJsSSgjHGmBBLCsYYY0IsKRhjjAmxpGCMMSbEkoIxxpiQ/weQr336iSVJIQAAAABJRU5ErkJggg==\n" }, "metadata": { "needs_background": "light" } } ] }, { "cell_type": "code", "source": [ "def jarvis_heat(t_rift,beta):\n", " a = 125_000.0\n", " dz = 1000.0\n", " n = int(a//dz)+1\n", " print(n)\n", " z = np.linspace(0.0,a,n)\n", "\n", " seg_por_ano = 365.0*24.0*60.0*60.0\n", "\n", " t = -t_rift\n", "\n", " kappa = 1.0E-6 #difusividade térmica\n", "\n", " G = np.log(beta)/(t_rift*seg_por_ano)\n", "\n", " dt = 1000.0 #passo de tempo em anos\n", " fac = dt*seg_por_ano*kappa/dz**2\n", " T1 = 1300.0\n", " T = z*T1/a\n", "\n", " k = 4.0 #W/m/K (condutividade)\n", " qq = [-k*(T[0]-T[1])/dz]\n", " tt = [-t_rift]\n", "\n", " tprint = t+100_000.0\n", "\n", " cor=\"r\"\n", " while t<=60.0E6:\n", " T[1:-1] += fac*(T[2:]-2*T[1:-1]+T[:-2]) + \\\n", " dt*seg_por_ano*G*z[1:-1]*(T[2:]-T[:-2])/(2*dz)\n", "\n", " t+=dt\n", "\n", " if(t>=tprint):\n", " qq = np.append(qq,-k*(T[0]-T[1])/dz)\n", " tt = np.append(tt,t)\n", " tprint+=100_000.0\n", "\n", " if (t>0):\n", " G=0.0\n", " cor = \"k\"\n", "\n", " return(tt,qq)" ], "metadata": { "id": "1tIGdedqHfRj" }, "execution_count": null, "outputs": [] }, { "cell_type": "code", "source": [ "tt,qq20 = jarvis_heat(20.0E6,2.0)\n" ], "metadata": { "id": "A5AJH5-kSkfQ" }, "execution_count": null, "outputs": [] } ] }