{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Python example" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## How to load the model" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import swami # SWAMI library\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "mcm = swami.MCM()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Single point" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Results are MCMOutput(dens=2.829354049223008e-10, temp=178.8555615350168, wmm=None, d_H=None, d_He=None, d_O=None, d_N2=None, d_O2=None, d_N=None, tinf=None, dens_unc=None, dens_std=1.749150709600944e-11, temp_std=9.60724715517875, xwind=-8.469296526223774, ywind=-9.41530132156296, xwind_std=28.02476670220009, ywind_std=21.61078819732512, alti=100.0, lati=3.0, longi=15.0, loct=12.0, doy=53.0, f107=70.0, f107m=69.0, kp1=1.0, kp2=1.0)\n", "\n", "Density is 2.829e-10 g/cm3\n", "Temperature is 178.86 K\n" ] } ], "source": [ "out = mcm.run(\n", " altitude=100,\n", " latitude=3,\n", " longitude=15,\n", " local_time=12,\n", " day_of_year=53,\n", " f107=70,\n", " f107m=69,\n", " kp1=1,\n", " kp2=1,\n", " get_winds=True,\n", " get_uncertainty=True\n", ")\n", "print(f\"Results are {out}\\n\")\n", "print(f\"Density is {out.dens:.3e} g/cm3\")\n", "print(f\"Temperature is {out.temp:.2f} K\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Altitude profile: temperature and density" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "altitudes = np.arange(0.0, 300, 10)\n", "temp = []\n", "dens = []\n", "\n", "for h in altitudes:\n", " out = mcm.run(\n", " altitude=h,\n", " latitude=3,\n", " longitude=15,\n", " local_time=12,\n", " day_of_year=53,\n", " f107=70,\n", " f107m=69,\n", " kp1=1,\n", " kp2=1,\n", " )\n", " dens.append(out.dens)\n", " temp.append(out.temp)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAELCAYAAAA2mZrgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAAA4S0lEQVR4nO3dd3xUVfr48c+TSqiBEAKEDglVWuhtQcWKvaGuig0LdndX2abfn+uu69pW116xInbEgoqE3pHeq/QSaoD05/fH3GRnIT2Z3DuT5/16zWtm7tzynORmnpxz7j1HVBVjjDEGIMztAIwxxniHJQVjjDEFLCkYY4wpYEnBGGNMAUsKxhhjClhSMMYYUyBgSUFEaojIfBFZKiIrReT/nOUNRORHEVnvPNf322asiGwQkbUicnagYjPGGFM4CdR9CiIiQC1VTReRSGAmcC9wKXBAVZ8QkYeB+qr6kIh0Aj4C+gBNgZ+AZFXNLeoYDRs21FatWgUk/tI4duwYtWrVcu34RfFiXF6MCUqOa9GiRftVNb4KQwLKd2579WdcFqFQBvB+OYo9r1U14A+gJrAY6AusBZo4y5sAa53XY4GxfttMBvoXt9+UlBR109SpU109flG8GJcXY1ItOS5goVbB38jJj/Kc2179GZdFKJRB1fvlKO68DmifgoiEi8gSYC/wo6rOAxJUdZeTkHYBjZzVE4Ftfptvd5YZY4ypIhGB3Ln6mn66i0gs8IWIdClmdSlsF6esJDIaGA2QkJBAampqJURaPunp6a4evyhejMuLMYG34qroue2lspRXKJQBgrwcRVUhKvsBPAL8Dms+CjgvxuXFmFSt+chrQqEMqt4vR3HndSCvPop3agiISAxwJrAGmAjc4Kx2A/CV83oiMFJEokWkNZAEzA9UfMYYY04VyOajJsA4EQnHd+nrBFWdJCJzgAkicjPwK3AFgKquFJEJwCogBxijxVx5ZIwxpvIFLCmo6jKgRyHL04AzitjmceDxQMVkjDGmeCF5R/O2A8f5fsUut8MwxhhXrNtzlCmr95Rr25BMCh8v2MaYD3/hWGaO26EYY0yV+2ThNu768JdybRuSSSGlVX1y85Sl2w65HYoxxlQ5VZDCLvIvhZBMCj1b1EcEFm496HYoxhhT5ZTCb/wqjZBMCvViIkluVMeSgjGmWvLVFMqXFkIyKYCvCemXrQfJzQvMgH/GGONVilpN4WS9WtbnaGYO6/YcdTsUY4ypUlqB9qMQTgoNAOtXMMZUT1ZTOEnzBjE0qhPNwi0H3A7FGGOqlKpan8LJRITerRqwYLMlBWNM9ZKnEGbNR6fq3ao+Ow9nsOPQCbdDMcaYKqMoYVZTOFXv1r5+BastGGOqkzy7ea1wHRrXpU50BPOtX8EYU434pqSxmsIpwsOEni3rW03BGFPNqPUpFKVP6was35vOwWNZbodijDFVIi/Pmo+K1LuV069gTUjGmGoiT62juUhdm9UjKiKM+daEZIypJnyXpFpSKFSNyHB6NI9lniUFY0w1oaqElfPbPeSTAkDfNnGs3HmYoxnZbodijDEBl2vNR8Xr27oBeWrjIBljqoc8hXBLCkXr2aI+EWFi/QrGmGohT9WuPipOTFQ4XZvVY96mNLdDMcaYgMvLs+ajEvVrE8ey7Yc5lpnjdijGGBNQuXlKeDnvXqtWSSEnT1lk/QrGmBBnSaEUUlr6+hXmWhOSMSbE5aoS4bWkICLNRWSqiKwWkZUicq+z/FER2SEiS5zHeX7bjBWRDSKyVkTOrsx4akVH0K15LHMsKRhjQlxFagoRlRyLvxzgQVVdLCJ1gEUi8qPz2bOq+pT/yiLSCRgJdAaaAj+JSLKq5lZWQP3aNOCVaZs4lplDrehAFt0YY9yTk+vB5iNV3aWqi53XR4HVQGIxm1wEjFfVTFXdDGwA+lRmTP3bNCQ3T20cJGNMSMtVDyYFfyLSCugBzHMW3SUiy0TkLRGp7yxLBLb5bbad4pNImaW0rE9kuFgTkjEmpHm1+QgAEakNfAbcp6pHRORl4DFAneengZsofEYILWR/o4HRAAkJCaSmppYpntZ1hR+WbKF/zJ4ybVeY9PT0Mh+/KngxLi/GBN6Kq6LntpfKUl6hUAZwvxwHDp2gTqSULwZVDdgDiAQmAw8U8XkrYIXzeiww1u+zyUD/4vafkpKiZfXMD2u19cOT9NDxrDJve7KpU6dWeB+B4MW4vBiTaslxAQs1gH8jRT3Kc2579WdcFqFQBlX3y3HOc9P15ncWFPl5ced1IK8+EuBNYLWqPuO3vInfapcAK5zXE4GRIhItIq2BJGB+Zcc1oG0ceYrd3WyMCVk5uXlERXiv+WggcB2wXESWOMv+CFwtIt3xNQ1tAW4DUNWVIjIBWIXvyqUxWolXHuXr0aI+NSLDmL0xjbM6N67s3RtjjOuyc/OIKOfY2QFLCqo6k8L7Cb4tZpvHgccDFRNAVEQYvVs1YPbG/YE8jDHGuCY7V4kML19SqDZ3NPsb2K4h6/aks+9optuhGGNMpcvOzSMy3MOXpHrNgLZxAFZbMMaEJF9SsJpCqXVuWo96MZHM2mBJwRgTejJz8oiOsKRQauFhwoC2ccxcvz//8ldjjAkJqkpGdi41IsPLtX21TAoAg5IasvNwBpv3H3M7FGOMqTTZuUqeQo1IqymUyaB2DQGYaU1IxpgQkpnju5Lfagpl1DKuFs0bxDBzvSUFY0zoyMjOAyDakkLZDWrXkDkb08jJzXM7FGOMqRQZ2b6agnU0l8OgdvEczcxh6fbDbodijDGVIj8pxFhNoewGtosjTGD6un1uh2KMMZUiPTMHgNrlnEisWieF2JpRdG0Wy4z1lhSMMaGhICnUsKRQLkOSGrJk2yEOH892OxRjjKmwY05SqBVlSaFchiTHk6c25IUxJjQczfAlhTpWUyifbs1jqRMdwXRrQjLGhID85qNa1qdQPpHhYQxoF8f0dTbkhTEm+BU0H0Xb1UflNiQ5nh2HTrBxnw15YYwJbodPZBMdEUZ0hCWFcvtNcjwAqWv3uhyJMcZUTNqxLOJqRZV7e0sKQLP6NWnXqDbT7H4FY0yQO3Asi7ja0eXe3pKCY2hyPPM2HeB4Vo7boRhjTLkdOJZFA6spVNzQ9o3Iys1jzsY0t0MxxphyS0u35qNK0bt1fWpGhZO61pqQjDHBK+1YptUUKkN0RDgD2saRum6vXZpqjAlKx7NyyMjOo0FtSwqVYmj7Rmw7YJemGmOC067DGQA0rluj3PuwpOBnWIdGAPy8Zo/LkRhjTNltO3AcgOYNapZ7H5YU/CTGxtChcR1+XmP3Kxhjgs+2gycAaF7fg0lBRJqLyFQRWS0iK0XkXmd5AxH5UUTWO8/1/bYZKyIbRGStiJwdqNiKc3qHRizYcpDDJ2zUVGNMcNl+4DhR4WE0quPN+xRygAdVtSPQDxgjIp2Ah4EpqpoETHHe43w2EugMnAO8JCLlu0+7Ak7v0IjcPLU5FowxQWfbweMk1o8hLEzKvY+AJQVV3aWqi53XR4HVQCJwETDOWW0ccLHz+iJgvKpmqupmYAPQJ1DxFaVHi/rE1ozk59XWhGSMCS7bDpygWf2YCu2jSvoURKQV0AOYBySo6i7wJQ6gkbNaIrDNb7PtzrIqFR4mDE2OJ3XdPnLz7NJUY0xwUFW2HTxeoU5mgPINuF0GIlIb+Ay4T1WPiBRZrSnsg1O+lUVkNDAaICEhgdTU1EqK9L+aaA4HjmXx1lc/k1S/6Bas9PT0gBy/orwYlxdjAm/FVdFz20tlKa9QKAO4U460E3kcOp5N2JHdpKZWYGQGVQ3YA4gEJgMP+C1bCzRxXjcB1jqvxwJj/dabDPQvbv8pKSkaCIdPZGnbsd/o379dVex6U6dODcjxK8qLcXkxJtWS4wIWagD/Rop6lOfc9urPuCxCoQyq7pTj+xW7tOVDk3TR1gMlrlvceR3Iq48EeBNYrarP+H00EbjBeX0D8JXf8pEiEi0irYEkYH6g4itO3RqR9GsTx0+r7H4FY0xwWLHjMOFhQqcmdSu0n0D2KQwErgNOF5ElzuM84AlguIisB4Y771HVlcAEYBXwPTBGVXMDGF+xhndKYOO+Y2zal+5WCMYYU2rLth8mqVFtakRW7KLNgPUpqOpMCu8nADijiG0eBx4PVExlcUbHRjwycSVTVu+lTXxtt8MxxpgiqSordhzm9A6NSl65BHZHcxGa1a9JxyZ1+dGakIwxHrfrcAZpx7Lo2qxehfdlSaEYwzs2YuHWA6SlZ7odijHGFGnxrwcB6JJYBUlBRMJF5EIRuUdEHsh/VPjIQeCszo3JU5hiN7IZYzxsyuq91K8ZyWlVkRSAr4FRQBxQx+8R8jo3rUtibAyTV+52OxRjjClUdm4eU1bv4YyOCUSEV7zxpzQdzc1UtWuFjxSERISzOzfm/blbSc/MoXZ0wO/1M8aYMpm36QBHMnI4q1NCpeyvNGnlOxE5q1KOFoTO6dKYrNw8UtdaE5Ixxnt+WLWbmMhwhiTHV8r+SpMU5gJfiMgJETkiIkdF5EilHD0IpLSsT1ytKL5fYU1IxhhvyctTfli5hyHJDSt8f0K+0iSFp4H+QE1VrauqdVS1YrfMBZHwMOGszglMXbOXjGzX7qUzxphTLN9xmN1HMjirU+NK22dpksJ6YIUzXka1dHbnxhzLymXm+v1uh2KMMQU+WbSNyHCplJvW8pWm53QXkCoi3wEFF+yfNJ5RSBvQtiH1YiL5dvkuzqykzhxjjKmI3YczmLBgO5enNKd+rahK229pksJm5xHlPKCQIa1DWVREGGd3TuC75bvJyM6ttLY7Y4wpr1embSRPlTuHtq3U/ZYmKYxT1S3+C0Skd6VGEQTOO60JExZuZ8b6/Qy32oIxxkV7j2bw0fxfubRnYoUn1TlZafoUPhORghnQRGQI8FalRhEEBrbzNSF9s2yn26EYY6q516dvIjs3jzuHtqv0fZcmKdwOfCkijZ2hr58Hzqv0SDwuMjyMczo35qfVdhWSMcY9+9MzeX/ur1zcPZFWDWtV+v5LTAqqugC4B/gBeBQYrqrbit0oRJ3ftQnpmTlMW7fP7VCMMdXU6zM2kZGTy5jTK7+WAMX0KYjI1/xvh3JN4DDwpoigqhcGJCIP6982jga1ovh66U7O7lx51wUbY0xpLNxygDdnbOaS7om0DdA8L8V1ND8VkCMGscjwMEZ0bcLHC7ZxNCPb7XCMMdXI/vRMxny4mMT6MTxyYeeAHafIpKCq0wJ21CB2cY9E3p2zle9X7KZyRhoxxpji5eYp947/hUPHs/nizj7Ui4kM2LFskp0y6tE8lpZxNflyyQ63QzHGVBPP/riOWRvSeOziLnRqGthRhiwplJGIcFH3RGZvTONgRp7b4RhjQtyU1Xv4z9QNXNWrOVf2ah7w4xWZFETkNRG5RESqxYQ6ZXFx96aowtxddmmqMSZwth04zv0fL6FTk7r830WB60fwV1xN4S2gG/CtiEwRkYdEpFuVROVxbeJr0615LLN35rgdijEmRB08lsXt7y8C4JXfplTZ8DpFJgVVnauqj6rqYOBK4FfgQRH5RUTeEpErqyRCj7qsZyLbjuaxYsdht0MxxoSY9XuOcvFLs1i/J51/j+xBi7jKHcqiOKXqU1DVNFX9SFWvV9UewItAUmBD87aLuiUSGQYfzf/V7VCMMSFk6pq9XPLSbI5l5vLR6H4Mq8RhsUujXB3NqrpIVR+v7GCCSb2akfRuHMFXS3ZyPMuakYwxFaOqvD59EzeNW0DLuJpMvGsgKS3rV3kcAbv6yGli2isiK/yWPSoiO0RkifM4z++zsSKyQUTWisjZgYqrMv2mWQTpmTlMWrbL7VCMMUEsMyeX33+6jMe/Xc05nRvzye39aRob40osgbwk9R3gnEKWP6uq3Z3HtwAi0gkYCXR2tnlJRDw/aUFy/TDaxtdivDUhGWPKaX96Jte8Po9PF23n3jOSePGantSMKs2sBoFRYlIQkZoi8hcRed15nyQiI0raTlWnAwdKGcdFwHhVzVTVzcAGoE8pt3WNiDCydwsW/3qIdXuOuh2OMSbIrNp5hIv+M4uVOw/zn2t6cP/wZMLCxNWYSlNTeBvfNJz9nffbgb9V4Jh3icgyp3kpv8EsEfAfeXW7s8zzLu2ZSGS4WIezMabUjmfl8OLUDVz+ymxy85RPbhvAiK5N3Q4LKN3Ma21V9SoRuRpAVU+ISHlT2cvAY/hGX30MeBq4CShsf4VO+Skio4HRAAkJCaSmppYzlIpLT09n+cI5pDQKY/y8LfSJ2UtMhLtZPj8uN38uhfFiTOCtuCp6bnupLOUVCmWAosuRk6ekbsth4sZsjmQp3eLDubGzkLbhF1I3VH2chVLVYh/AbCAGWOy8bwvML2k7Z91WwIqSPgPGAmP9PpsM9C9p/ykpKeqmqVOnqqrq4q0HtOVDk/SdWZtdjSdfflxe4sWYVEuOC1iopTjXK/tRnnPbqz/jsgiFMqieWo6c3Dz9dOE2HfjEFG350CS94pXZumBzmjvBafHndWlqCo8A3wPNReQDYCAwqjwJSESaqGr+pTqXAPlXJk0EPhSRZ4Cm+O6BmF+eY7ihR4v6dGseyzuzt3Bdv5autwkaY7xBVflh1R6e/mEt6/ak07lpXf52cRd+kxxP+RtcAqvEpKCqP4rIYqAfvmaee1V1f0nbichHwFCgoYhsx5dchopId3xNQ1uA25xjrBSRCcAqIAcYo6pBNbDQTQNbce/4JUxbv49h7av2ZhNjjPfM3rCfJyevZcm2Q7SJr8WL1/Tk3C6NPf9PY3Ezr/U8aVH+f/gtRKSFqi4ubseqenUhi98sZv3HgaC9Ie7cLk14vM5q3p61xZKCMdXY0m2H+NeCE6z8fh5N69Xgycu6cmnPRCLCg2NQ6uJqCk87zzWAXsBSfDWFrsA8YFBgQwsuURFhXNevJU//uI4Ne9Np1ygwU+UZY7xp/Z6jPPXDWiav3EOdSPjLiE5c27dFlQ1kV1mKm3ltGICIjAdGq+py530X4HdVE15wubpvC16YuoHXp2/in5d3dTscY0yA7TmSwbfLd/Ht8l0s2HKQ2tER3H9mMsls59xBrd0Or1xK09HcIT8hAKjqCqdfwJykYe1orunTgnfnbOGWwa1JSrCpKIwJNXuOZPDd8l18u3w3C7YeQBU6NK7D785K5pq+LWlQK4rU1OCdmbE0SWG1iLwBvI+vg/i3wOqARhXE7jkjic8WbeeJ79bw5qjebodjjKkEe49k8N2K3XyzfBcLtvw3Edx/ZjLnndYkpJqLS5MUbgTuAO513k/HdxOaKUSDWlHcMawtT36/ljkb0+jfNs7tkIwx5bD3aAbfr9jNN8t2Md9JBO0T6nDfGcmc37Ux7RqFZktAaS5JzQCedR6mFG4a2Jr352zlH9+t5ss7B3r+EjRjjM++o5l8v3I33yzbybzNvkSQnFCbe89I4vzTmlSLJuESk4KIbKaQISdUtU1AIgoBNSLDefCs9jz4yVK+XraTi7oHxTBOxlRL+9MzC2oE8zankafQrlFt7jk9ifO7NiG5GiQCf6VpPurl97oGcAXQIDDhhI6LeyTyxszN/GvyWs7p0pjoiOC6LM2YUKWqbD94gunr9/HNsl3M3eRLBG3ja3HX6UmMqIaJwF9pmo/STlr0nIjMBP4amJBCQ3iY8MfzOnDdm/N5avJa/nR+J7dDMqZaOnw8myXbD7F0m/PYfoj96VkAtImvxV3D2nFe1ya0T6jj2aEnqlJpmo/872wOw1dzqL5ptAwGJ8Vzff+WvD5jMz1a1Oe805q4HZIxIS0zJ5dVO484X/6HWbLtEJv3HwNABNrF12Zo+0Z0ax5Ln1YNSE6obYngJKVpPnra73UOsBm4MjDhhJ4/n9+J5TsO8/tPlpKcUDtkr1gwpqrl5Smb046xdNshlji1gFW7jpCd6+sCbVQnmu7NY7k8pRk9msfSpVk96taIdDlq7ytNUrhZVTf5LxCR4LxVzwVREWG8dG1PLnhhJre9t4iv7hpE7Wj3ptozJljtO5pZ8OW/1GkOOpKRA0CtqHBOa1aPmwa1pkfzWLo1j6VJPXfmOA52pfl2+hQ4eXC8T4GUyg8nNDWpF8MLV/fkt2/O4/efLOWla3taldWYYhzLzGHFjsMs3Z5fCzjMjkMnAF9/XfuEOozo1pTuzXwJoF2j2oTbpd+VorhRUjsAnYF6InKp30d18V2FZMqgf9s4Hj6nA49/u5rXZ2xi9JC2bodkjCfk5Oaxfm86S7Yd4rsVmTyxZDrr9hwlz7kQvln9GHq0iOXGga3o1jyWLk3rERNlV/MFSnE1hfbACCAWuMBv+VHg1gDGFLJuGdyaX7Yd5Inv1tC5aT0GtmvodkjGVLmjGdnM2ZjGwq0HWfLrIZbvOMyJbN/0KbUiIaV1Dc7qlED3FrF0bRZLw9rRLkdcvRQ3SupXwFci0l9V51RhTCFLRHjy8m6s35POTe8s4Okru3lmsm5jAiU3T1m+4zAz1u1jxvr9LP71IDl5SlR4GJ2a1uWq3s3p7vQDbFk+n2HD+rgdcrVWXPPRH1T1SeAaETllwhxVvSegkYWo2tERjB/dj9veW8RdH/7Cpn3HuPv0dtbHYELKzkMnmLF+H9PX72fWhv0cOp4NQJfEutw6pA1DkuLp2TL2lJs6t9rfgeuKaz7KHwl1YVUEUp3E1Y7mg1v7Mvaz5Tzz4zo27kvnn5d1DbrJOIzJdzwrh3mbDjB9va82sGFvOgAJdaM5s2MCg5MaMqhdQ+KsKcjzims++tp5eVxVP/H/TESuCGhU1UB0RDhPX9mNto1q86/Ja9l24DivXd/L2k9NUMjLU1btOsKM9fuZsX4fC7ccJCs3j+iIMPq2iWNk7+YMToq3m8OCUGkuSR0LfFKKZaaMRIQxw9rRumEt7v94CRe/OIs3b+hN+8Z2g5vxnr1HMgqSwMwN+wuGiujQuA6jBrZicFJDerdqYDXeIFdcn8K5wHlAoog87/dRXXx3NptKct5pTUiMjeGWdxdy2cuzeeGaHgxr38jtsEw1l5Gdy4ItB5ixfj/T1+1jze6jAMTVimJwUkMGJ8UzOKkhjeraFeqhpLiawk5gEXCh85zvKHB/IIOqjro1j+WrMQO5ZdxCbn5nAX8d0YlRA+3GcVN1VJV1e9ILOojnbUojMyePqPAwerWqz0PndGBIckM6Nq5rc4SEsOL6FJYCS0XkfVW1mkEVaBobwye39+fe8Ut49OtVbNp/jL+O6EREeJjboZkQlZaeycwN+wuahfYcyQR8w0hf07cFQ5Lj6du6ATWjbGiW6qK45qPlOJPrnNRRJICqatfAhlY91YqO4NXrUvjn92t4bfomNu8/xovX9rSBvEylyMrJY9HWg8xwrhJasfMwqhBbM5JB7RoyJCmeQUkNaRpr4wZVV8Wl/xFVFoX5H765GDrSpmEt/vzlCi57aTZv3tCbFnE13Q7NBBlVZdP+Y8xY52sSmrspjeNZuUSECT1b1ufB4ckMToqnS2I9GzvIAMU3H20tbLmIDASuAcYEKijjM7JPC1rE1eSO9xdz8UuzePW6FHq3sknvTPEOH89m1kZf5/CM9fsLBpJrFVeTy1OaMTgpnn5tGlDHap+mEKVqKBSR7vgSwZX45lP4vBTbvIWvtrFXVbs4yxoAHwOtgC3Alap60PlsLHAzkAvco6qTy1aU0DSgbUO+uHMAN49byLWvz+OJy07j0p7N3A7LeEh2bh5Ltx3i8/VZPLdyFsu2HyJPoU6NCAa2bcidw9oyuF281TRNqRTXp5AMjASuBtLwfZmLqg4r5b7fAf4DvOu37GFgiqo+ISIPO+8fEpFOzrE6A02Bn0QkWVVzy1iekNQmvjZf3DmA299fxAMTlrJp3zEeGJ5sV4BUY7+mHWf6+n1MX7ePORvTOJqZgwA9WsDdpycxJLkh3ZrF2kUKpsyKqymsAWYAF6jqBgARKfWlqKo6XURanbT4ImCo83ockAo85Cwfr6qZwGYR2QD0AWwgPkdszSjevakvf/lyBf+ZuoF1e47yj0tPs2EDqpm1u49y+/uLCqaYTIyNYUS3pgxJakje7jWcP3ygyxGaYFdcUrgM33/vU0Xke2A8viuPKiJBVXcBqOouEcm/QysRmOu33nZn2SlEZDQwGiAhIYHU1NQKhlR+6enpVX78c+KUsA5RTFi9h6FP7uH6TtH0bvy/v0Y34iqJF2MCb8VVmnP7p63ZbN6fxZXtI+nZKIKEmoJIGqSlkZ55zDNlKS8v/T4qIqjLoarFPoBawLXAJOA48DJwVknbOdu2Alb4vT900ucHnecXgd/6LX8TuKyk/aekpKibpk6d6tqx1+w6oiOen6EtH5qkYz5YpGnpmZ6IqyhejEm15LiAhVqKc72yH0Wd28/9uE5bPjRJs3Jyy1yWYBAKZVD1fjmKO69LbHBU1WOq+oGqjgCaAUvw9QWUxx4RaQLgPO91lm8Hmvut1wzfHdWmCO0b1+HzOwfw4PBkJq/czVnPTuP7FbvdDssE2MHjWdSpEUGk9RWYACnTmaWqB1T1VVU9vZzHmwjc4Ly+AfjKb/lIEYkWkdZAEjC/nMeoNiLDw7j7jCQm3jWIhLo1uP39Rdzz0S+kZ6nboZkAOXg8iwa1otwOw4SwgN27LiIf4etUbigi24FHgCeACSJyM/ArcAWAqq4UkQnAKnyD7Y1Ru/Ko1Do2qcuXYwby0tSNvPDzelIjICpxN2d1bux2aKaSHTiWRWxNSwomcAKWFFT1lNnaHGcUsf7jwOOBiifURYaHce+ZSQzvlMDtb89i9HuLuLh7Ux69sLN9iYSQQ8ezaVjbfp8mcKxhMsR0alqXv/avwb1nJDFp2S6GPzudSct25nfgmyB38LjVFExgWVIIQRFhwv3Dk/lyzEAS6kZz14e/cNM7C9h24LjboZkKUoUwm8nMBJAlhRDWJbEeX945kL+M6MS8zQcY/uw0Xp22kezcPLdDM+UUGS7k5NnvzwSOJYUQFxEexs2DWvPTA79hULt4/vHdGi54YSa//HrQ7dBMOUSGh1lSNwFlSaGaaBobwxs39OKV36Zw6Hg2l748m79+tYIjGdluh2bKIDI8jKwc6x8ygWNJoZo5p0tjfnxgCDf0b8V7c7cy/JlpfLd8l3VEB4nIcLGaggkoSwrVUJ0akTx6YWe+vHMgcbWiueODxdwybiHbD1pHtNdFhodZn4IJKEsK1Vi35rFMvGsgfzqvI7M3pnHuczP4YaUNleFldWpEsHrXUVbtPOJ2KCZEWVKo5iLCw7h1SBt+uH8IreNrMfq9RTz5/Rpy86w5yYv+cE4HosLDuOKV2Uxds7fkDYwpI0sKBoDmDWoy4bb+XN2nOS+lbuT6t+aRlp7pdljmJB2b1OWruwbSOr4WN49bwLjZW9wOyYQYSwqmQI3IcP5xaVeevKwrC7YcZMQLM1my7ZDbYZmTJNStwYTb+nN6hwQembiSRyeutJqdqTSWFMwpruzdnM/vGEB4mHDFK7N5f+5WuzrJY2pGRfDqdSncMqg178zewq3vLuSEXapqKoElBVOoLon1mHT3IAa0bcifv1zBg58s5USWDVzrJeFhwp9HdOLxS7owbd0+/j4vg52HTrgdlglylhRMkWJrRvH2qN7ce0YSX/yyg0tfns2vaXbZqtdc27clb4/qzf4TeVz04iyWbz/sdkgmiFlSMMUKcwbXe2tUb3YeOsGlL89ixQ770vGaIcnx/LlvDFHhYVz12hxS19qVSaZ8LCmYUhnWvhGf3TGA6IhwRr42l9kb97sdkjlJYp0wvhgzgFZxtbhl3EI+X7zd7ZBMELKkYEqtXaPafHpHf5rG1mDUWwv4fsUut0MyJ2lUpwYf39aPvm0a8MCEpbwybaNdJGDKxJKCKZMm9WKYcFt/uiTW5c4PFvPhvF/dDsmcpE6NSN4e1YcLuzXlie/W8P8mrSLPLlk1pWRJwZRZbM0oPrilH79JjuePXyznhSnr7b9Rj4mKCOO5q7pz86DWvD1rC3eP/4XMHLt6zJTMkoIpl5iocF67vheX9kjk6R/X8X9f23+jXhMWJvxlRCf+eF4Hvlm2i1FvLbCh0k2JLCmYcosMD+OpK7oV3EB138dLbFhnDxo9pC3PXtWNBVsOcOUrc9hvw5eYYlhSMBUSFib86fyOPHROByYu3cldHy4mK8cSg9dc0qMZb43qzZa0Y9zx/iL7HZkiWVIwFSYi3DG0LY9e0InJK/dw5weLrP3ag4Ykx/Pk5d1YsOUgj3690u1wjEdZUjCVZtTA1jx2UWd+Wr2X299bREa2JQavubBbU+4Y2pYP5/3K+3O3uh2O8SBLCqZSXde/FX+/5DSmrt3HbZYYPOl3Z7Xn9A6NeHTiSuZuSnM7HOMxriQFEdkiIstFZImILHSWNRCRH0VkvfNc343YTMVd07cFT17Wlenr9/lG77SB9DwlPEx4bmR3WsTV5M4PFts0rOZ/uFlTGKaq3VW1l/P+YWCKqiYBU5z3Jkhd2bs5T13ejZkb9nPzuAVk2rDOnlK3RiRvXN+L7Nw8Rr+7iONZOW6HZDzCS81HFwHjnNfjgIvdC8VUhstSmvHsld2ZuymNZxZlkJ5pXzxe0ia+Ni9c3YM1u4/w6ETreDY+4sadqCKyGTgIKPCqqr4mIodUNdZvnYOqekoTkoiMBkYDJCQkpIwfP76Koj5Veno6tWvXdu34RfFaXPN25fDq0gzaxIbzQEoNakaK2yEVKOlnNWzYsEV+tdmAqui5Xd7f+1MLMziUkcffBtUs87aVzWvnbnl5vRzFnteqWuUPoKnz3AhYCgwBDp20zsGS9pOSkqJumjp1qqvHL4oX43pq/I/a7o/f6AUvzNCDxzLdDqdAST8rYKG68DdSnnO7vL/3s5+dpje/s6Bc21Y2L5675eH1chR3XrvSfKSqO53nvcAXQB9gj4g0AXCebUD4EJKS4Js+cs3uo1z9+jzS7K5aT1BVth04TvMGMW6HYjyiypOCiNQSkTr5r4GzgBXAROAGZ7UbgK+qOjYTWKd3SODNG3qxeX86I1+by94jGW6HVO0dOJbFsaxcWjRwv+nIeIMbNYUEYKaILAXmA9+o6vfAE8BwEVkPDHfemxAzOCmed27sw45DJ7jqtbk2p7DLfj3guxzVkoLJV+VJQVU3qWo359FZVR93lqep6hmqmuQ8H6jq2EzV6Ncmjvdu7sP+o5lc+eocth2w6+Tdkp8UmltSMA4vXZJqqpGUlg344Na+HM3I4cpX59i8zy7Iy1Omrd0HQLP61qdgfCwpGNd0bRbLR7f2Q4DLX5nNV0t2uB1StZGemcNt7y/i8192MGpAK2pGRbgdkvEISwrGVZ2a1mXi3YPomhjLveOX8I9vV5Nrk/UE1LYDx7n85dlMWb2HRy7oxCMXdHI7JOMhlhSM6xrWjub9W/pyXb+WvDp9E6Pens/h4zZDWCAs2HKAi16cxY5DJ3jnxj7cOLA1It65mdC4z5KC8YSoiDAeu7gLT1x6GnM3pXHhizNZt+eo22GFlI8X/Mo1r88lNiaSL8cMZEhyvNshGQ+ypGA8ZWSfFowf3Y/jWblc8uIsJq/c7XZIQS8nN4/HJq3ioc+W069NHF/cOZC28d4dgsG4y5KC8ZyUlg34+q5BtEuow23vLeLZH9eRZ/0M5XL4RDY3jVvImzM3M2pAK94e1Zt6NSPdDst4mCUF40mN69Xg49H9uDylGf+esp6bxi2woTHKKD0zh2ten8vsDfv5+yWn8eiFnYkItz95Uzw7Q4xn1YgM51+Xd+Wxi7swe2Ma5/57BrM37nc7rKCQk5vHXR8uZs3uo7x+fS+u6dvC7ZBMkLCkYDxNRLiuX0u+GjOQOjUiuPaNeTzzw1pycvPcDs2zVJVHJq4kde0+HruoC8M6NHI7JBNELCmYoNCxSV2+vnsQl/dsxvM/b+Dq123cpKK8Nn0TH8z7lTuGtrUagikzSwomaNSMiuBfV3Tjuau6s2rnEc799wx+sKuT/sc3y3bxj+/WMKJrE35/Vnu3wzFByJKCCToX90jkm3sG07xBDKPfW8SjE1eSkZ3rdliuW38wl/snLKFXy/o8dUU3wsLspjRTdpYUTFBq1bAWn90xgJsHtead2Vu49KXZbNqX7nZYrtmy/xjPL84gMTaG167vRY3IcLdDMkHKkoIJWtER4fxlRCfevKEXuw6f4KIXZzF93T63w3LFu3O2cjQbzunSmPp2H4KpAEsKJuid0TGBSfcMJjE2hhvfWcB7c7e6HVKVGzOsLV3jw3k5dSO3vruQA8ey3A7JBClLCiYkJMbG8OkdAxiaHM9fvlzBoxNXVqvLVuNqR3N/z2j+OqIT09ft59x/T2f2Brunw5SdJQUTMmpHR/Da9b0K+hlueXchRzOqz2irIsJNg1rzxZgB1IqO4No35/HP79eQXY2So6k4SwompISHCX8Z0YnHL+nCjPX7uezl2dVuus/OTesx6e5BXNWrOS+nbuTyV+awNe2Y22GZIGFJwYSka/u2ZNyNfdh1OIOLX5zFoq0H3Q6pStWMiuCJy7ry0rU92bwvnfOfn8mXv9jMdqZklhRMyBqU1JAv7hxI7RoRXP36XCYu3el2SFXuvNOa8N19Q+jYpA73fbyEBz5ewt4jGW6HZTzMkoIJae0a1ebLOwfSvVks93+8hDkb09wOqcolxsbw0a39uO/MJL5csoNB/5zKQ58uY8Pe6ntfhymaJQUT8urXiuKNUb1oGVeTMR8uZkc1HDMpIjyM+85MZurvhnJV7+Z8uWQHZz4zjVvfXciirQfcDs94iCUFUy3UrRHJ69f3Ijsnj9HvLuREVvUcFqNlXC3fUOQPn849ZySxYMsBLnt5Dpe/PJsfV+2xyYyMJQVTfbSNr81zI7uzatcRxn6+DNXq+wUYVzuaB4YnM/vh03n0gk7sOpzBre8uZPiz05iwYBuZOdUzaRoPJgUROUdE1orIBhF52O14TGg5o2MCD5yZzJdLdvLmzM1uh+O6mlERjBrYmmm/H8q/R3YnOiKcP3y2jMH/nMor0zZypBrd52F8PJUURCQceBE4F+gEXC0indyNyoSau05vx7ldGvP3b1ezcr/9Rwy+PoeLuifyzT2DeO/mPiQn1OGJ79Yw4B8/849vV7PHrliqNjyVFIA+wAZV3aSqWcB44CKXYzIhRkR46opuJDWqw0tLM/g1rXrd3FYcEWFwUjzv39KXSXcPYliHRrw+YxOD/vkzv/9kqY2pVA1EuB3ASRKBbX7vtwN9/VcQkdHAaICEhARSU1OrLLiTpaenu3r8ongxLi/GdHP7PP4+V/l8ymy6N3L/T6Gi53YgfsaXNYHB9WKYvCWbqat2MLzBAaLCAzdPgxfPk/II6nKoqmcewBXAG37vrwNeKGr9lJQUddPUqVNdPX5RvBiXF2NSVZ3808/Ffg4sVBf+Fspzbgf6Z5yVkxvQ/at69zwpK6+Xo7jz2mvNR9uB5n7vmwHV7zZUU2UC+V9vqIkM99rXhQkEr/2WFwBJItJaRKKAkcBEl2Myxphqw/2GVD+qmiMidwGTgXDgLVVd6XJYxhhTbXgqKQCo6rfAt27HYYwx1ZHXmo+MMca4yJKCMcaYApYUjDHGFLCkYIwxpoBoEI8UKSL7gK0uhtAQ2O/i8Yvixbi8GBOUHFdLVY2vqmDylfPc9urPuCxCoQzg/XIUeV4HdVJwm4gsVNVebsdxMi/G5cWYwLtxlUcolCUUygDBXQ5rPjLGGFPAkoIxxpgClhQq5jW3AyiCF+PyYkzg3bjKIxTKEgplgCAuh/UpGGOMKWA1BWOMMQUsKRhjjClgScGYICAizUVkqoisFpGVInKvs7yBiPwoIuud5/p+24wVkQ0islZEznYv+v8lIuEi8ouITHLeB2MZYkXkUxFZ4/xO+gdjOQpjSaGSiUgbEXlTRD71W3axiLwuIl+JyFkux9dJRCaIyMsicrmbseQTkTAReVxEXhCRG9yOJ5+IdBSRV5w//jtcDicHeFBVOwL9gDEi0gl4GJiiqknAFOc9zmcjgc7AOcBLIhLuSuSnuhdY7fc+GMvwb+B7Ve0AdMNXnmAsx6mKmpKtOj6At4C9wIqTlp8DrAU2AA+Xcl+fFrKsPvCmm/EBDwKDndcTvfAzAy4B3gGeAc7w4O8yrCK/t0A8gK+A4U5ZmjjLmgBrnddjgbF+608G+nsg7mb4vjBPByY5y4KtDHWBzTgX6vgtD6pyFPWwmsL/egffl0YBJ6O/CJwLdAKudv7bPk1EJp30aFTC/v/s7MvN+N4DRorIv4C4CsRSmTG1B+ao6gNAZf1HXhlxISIXAjPxfZF5goi0AnoA84AEVd0F4Dznn4OJwDa/zbY7y9z2HPAHIM9vWbCVoQ2wD3jbaQZ7Q0RqEXzlKJTnJtlxk6pOd/7g/PUBNqjqJgARGQ9cpKr/AEaUZr8iIsATwHequtgD8Y1xviA/L28slRmTiGwHspy3uRWNqbLicvYzEZgoIt8AH1ZGbBUhIrWBz4D7VPWI79QqfNVClrl6/bmIjAD2quoiERlamk0KWeaFa+gjgJ7A3ao6T0T+jdNUVASvlqNQVlMoWZmyvIjEicgrQA8RGessvhs4E7hcRG53Ob5WIvIa8C7wr0qOpVwx4UtOZ4vIC8D0AMVU5rhEZKiIPC8ir+KB2QBFJBJfQvhAVfMT+h4RaeJ83gRfkxn4ytbcb/NmwM6qirUIA4ELRWQLMB44XUTeJ7jKAL64tqvqPOf9p/iSRLCVo1BWUyhZmbK8qqYBt5+07Hng+UqOK19Z49sCjA5QLPnKGtNx4ObAhVOgrHGlAqmBCqYsnNrmm8BqVX3G76OJwA34aqI34OtryF/+oYg8AzQFkoD5VRfxqVR1LL72dZyawu9U9bdOU2ZQlAFAVXeLyDYRaa+qa4EzgFXOI2jKURRLCiXzepb3YnxejAm8G1dpDASuA5aLyBJn2R/xfQFNEJGbgV+BKwBUdaWITMD3RZUDjFHVSmmaC4BgLMPdwAciEgVsAm7E1/ISbOU4hQ1zcRKnHXqSqnZx3kcA6/D9N7ADWABco6orLT7vxuTluIzxMutT8CMiHwFzgPYisl1EblbVHOAufJeRrQYmuJgQPBefF2PyclzGeJ3VFIwxxhSwmoIxxpgClhSMMcYUsKRgjDGmgCUFY0yJnJsylziP3SKyw+99lNvx+XNuOhwQ4GO8KiIDT1r2qIj8znldQ3wjpT4SyDgCwe5TMMaUyLkpszv4vvyAdFV9yq14RCTCuZqsMEOBdGB2GfYXXsZ7B/oCdxaxryh8d54vUtX/K8M+PcFqCsaYchGRFBGZJiKLRGSy3xAPqSLyrIhMF99cA71F5HPxzTPwN2edVuKbi2CciCwT3/DkNUux37+LyDTgXhG5QETmiW9Qup9EJMG5N+V24H6nFjNYRN4Rv2HiRSTdeR4qvjkqPsR3U2C4iPxLRBY4Md1WRLk7AuuKSCIR+IbwWK+qxY2H5FmWFAJERHKdk3KliCwVkQdEpFJ/3iJyu4hc77weJSJNy7GPq0XkT2XcZqyIXFvWY520j9tF5EURea8i+zGuEeAF4HJVTcE3VPnjfp9nqeoQ4BV8wz2MAboAo0Qkf3Te9sBrqtoVOALcKb7xnYrbb6yq/kZVn8Y3em0/Ve2B74v4D84wLq8Az6pqd1WdUUI5+gB/UtVO+IZaOayqvYHewK0i0rqQbc4Fvi9if38AclT1vhKO61nWfBQ4J1S1O4D4hmH+EKgHVFobo6q+4vd2FLCCsg/bcA5lH5fpLODKMm5zso+A3wDRFdyPcUc0vi/5H8U3Ums4sMvv84nO83JgZf6Q0iKyCd9QI4eAbao6y1nvfeAefF+2xe33Y7/XzYCPnZpEFL45Dspqvqrmb3cW0NWvVlEP3zhFJ+/3bHzDWhRmJtBfRJJVdV054nGd1RSqgKruxTcI3V3iU2g11anOpsp/p/n7QJy/DBF5QkRWOes/5Sx7VER+55zEvfCNxbJERM4XkS/yjy8iw0XklGGynX13BxaftLym+GZnWyYiHztV9F7OZ3WBKFXd51TXv3BqQktFZIBfs8AbIrLCKcOZIjLLaT7o4/xMDqvq1cBhEYmp9B+6CTTB92Xf3Xmcpqr+swpmOs95fq/z3+f/M3rynbNaiv0e83v9AvAfVT0NuA2oUUSsOTjfdc45798x7r8/wTccdv6xW6vqD/9TaF8TV6yqFvXP13TgPuC78tTcvcCSQhVxxvAPwzfxRnHV1B74TqpO+CbzGCgiDfDNTtbZqWr/7aR9fwosBK51aiffAh1FJN5Z5Ubg7ULC6gEs1VNva78TOOgc6zEgxe+zM/nvhDPPA9NUtRu+oYPzh4xoh2+6wq5AB+AaYBDwO3yDuCEiT4tvCO/6/HcuBRM8MoF4EekPvmG9RaRzGffRIn974Gp8/2WvLcN+6+Ebwwp8o5LmOwrU8Xu/hf+ewxcBkUXsbzJwh9OEhYgki2/yHH/DgKnFFUpVP8M3LP33IhJb3LpeZEmhauUP3XwWcL34Rruch28GtCTns/mqul1V84AlQCt87a0ZwBsicilwvLiDOF/y7wG/dU7K/sB3hax6ThHLB+Fro0VVVwDLitjmdOBlZ71cVT3sLN+sqsudMqzEN2+t4mtKaOWs/6CqjlbVi708YqQpUh5wOfBPEVmK71wt62Wgq4EbRGQZ0AB4WVWzyrDfR4FPRGQGsN9v+dfAJfkdzcDrwG9EZD6+q4aOnbInnzfwjWS6WERWAK9yahN7cf0JBZym3c/xTdBUVA3Gm9QDc4KG4gPfJXv+79sAafgSw2fA2YVsMxRn3lrn/X+AUc7raOA8fJPj/OwsexTfmPTgG/e/l9+2TYFF+Ka3fLKIGFOBuEKWfwUM83u/OH/f+P5Iw53X+4Dok7Zthd+8yPimxby8sM/sUX0fwXouOH8LkW7HEciH1RSqgNOM8wq+9k+ldNVU/+1rA/VU9Vt8TUvdC1ntf6rM6mvz3IlvXuh3CtlnPSBCfdefn2wmTkeyiHQCTnNedwbW6H//s5+CM6ey009St6gyGBMKVLWnqma7HUcg2dVHgRPjNA9F4uvoeg/InzHrDXz/KS12Or72ARcXs686wFdONVSA+wtZ5x3gFRE5AfRX1RPAB0C8qq4qZP3hwE9FHO8lYJxTrf8FX/PRYXztsf5V53uB18Q3qUguvgSxC2NKoL5LR7u4HYc5lQ2dHcJE5D/AL6r6ZiGfvQG8oapzC/ksHF8VOUNE2uKrESQD3wDXq3N5oTEm9FhSCFEisghfh9pwVc0saf2Ttq2D7wqLSHw1k4dUtbAOaWNMiLGkYIwxpoB1NBtjjClgScEYY0wBSwrGGGMKWFIwxhhTwJKCMcaYApYUjDHGFPj/1q+oeXN5mf4AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "f, ax = plt.subplots(1, 2, sharey=True)\n", "\n", "ax[0].plot(dens, altitudes)\n", "ax[0].set_xscale(\"log\")\n", "ax[0].set_xlabel(\"Density / g/cm³\")\n", "ax[0].set_ylabel(\"Altitude / km\")\n", "\n", "ax[1].plot(temp, altitudes)\n", "ax[1].set_xlabel(\"Temperature / K\")\n", "\n", "ax[0].grid(True)\n", "ax[1].grid(True)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Map at altitude" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Temperature and density at 160 km" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "tags": [] }, "outputs": [], "source": [ "lati = np.arange(-90, 90, 10)\n", "loct = np.arange(0, 24, 3)\n", "\n", "temp = np.zeros((len(lati), len(loct)))\n", "dens = np.zeros((len(lati), len(loct)))\n", "\n", "for i, lat in enumerate(lati):\n", " for j, lt in enumerate(loct):\n", " out = mcm.run(\n", " altitude=160,\n", " latitude=lat,\n", " longitude=15,\n", " local_time=lt,\n", " day_of_year=53,\n", " f107=70,\n", " f107m=69,\n", " kp1=1,\n", " kp2=1,\n", " )\n", " dens[i,j] = out.dens\n", " temp[i,j] = out.temp" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 0, 'Local time / h')" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA2wAAAFNCAYAAAB8A+zrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAABY+klEQVR4nO3debwkZX3v8c8Xhm0YEHBYZ0AGMqxeRSDEJRoUF8AF9UYvxAVRL+FexT0KMYlGQ2JcQ+JCJopLRBFRFL0jghpBoyDbyDagA4MwAwyOyOawDfO7f1Q11DTdffr0qb2+79frvM7pqurqp+uc009/+/fUU4oIzMzMzMzMrH42qLoBZmZmZmZmNpgDm5mZmZmZWU05sJmZmZmZmdWUA5uZmZmZmVlNObCZmZmZmZnVlAObmZmZmZlZTTmwWetJeqak66puh5mZmZnZdDmwWWkk3SjpPkn3SLpT0s8kHSep0L/DiPhJROzZ147nzmSfkv5C0ldmcP9fSdpjJm3o29+zJV2ZHtffSTpL0ry89m9mZsWQdG/ma13aT/Zuv6rq9k0ij352wsc9V9LzByz/gqR/yNzeV9Ktkt5ZbgvNJuPAZmV7cURsATwB+BDwHuBz1TZpIocDiye5o6TdgQ0i4lc5tuca4AURsRWwE/Br4DM57t/MzAoQEXN6X8BNJP1kb9lpVbevn6RZdXwMSZsDBwDnT7HdfsB/ASdFxMcmaqBZyRzYrBIRcVdEnA38L+BoSU8EkLSJpI9KuknSKkmnSNosXXewpBWS3inp9vTTsWN6+5R0uKRr0greSknvyt4v/fk/gV2A76SfXr5b0v+TdHy2fZKukPTSQW1PK4LPA84Zsv7dadtukfRGSSHpjzKbvJA07EnaTNLHJP1G0l2Sfpou2zW93zGSbpb0+7Qa+cdp2+6U9MnM8VwVEbdkHuNhIPuYZmbWIJI2kHSCpOvTkRNnSNomXTetPkLS6yT9t6R/S/uaayUdkln/OEmfS/uulZL+QdKGfff9hKQ7gPdL2l3Sj9J2rZZ0mqSt0u0H9bOP9MOZx3ykCifp/ZLOlPRlSXcDrxvVpiEOAf47Ih4YcUwPAn4A/HVEfHLYdmZ148BmlYqIXwArgGemi/4Z2APYjyRwzAP+LnOXHYDHpcvfAHxK0tbpus8Bf5lW8J4I/GjA472G9T/B/DDwReDVvW0kPTnd/7AK2kHADRGxun+FpEOBdwDPTdv/ZwPufzjw/9KfP0ryieDTgW2AdwPrMtv+CbCQJNj+C/DedN/7Aq+U9Mj+Je0i6U7gPuBdwIeHtN/MzOrvLcBLSfqRnYDfA5/q22bsPiLd9gZgLvA+4Ju9AEjSD64l6beeAjwfeOOA+24HnAQI+Ke0XXsDOwPvh6H97DiOAM4EtgJOG6NN/bJ96yAHkXzQ+vaI+OyYbTKrBQc2q4NbgG0kCfjfJC+md0TEPcA/Akdmtn0I+EBEPBQRi4F7gT0z6/aRtGVE/D4iLhvz8b8NLJS0ML39GuBrEfHgkO0fqZAN8Erg8xFxdUSsAf4+u1LSbOCPgfPTSt3rgbdGxMqIeDgiftb36eAHI+L+iDgX+APw1Yi4PSJWAj8h6cQAiIib0iGRc4G/Aa4d8/lbxSSdmlaNr8ppf+ekn7B/t2/5aZKuk3RV+pgb5fF4ZlaIvwTeGxEr0n7h/cCfa/3hgmP3EcDtwL+k/efXgOuAF0raHjgMeFtE/CEibgc+wfp97y0R8W8RsTYi7ouIZRFxXkQ8EBG/BT7O4A8op+PnEfGtiFgHbDlGm/odxuhTFZ4K3AV8b4btNCudA5vVwTzgDmBbYDZwafpm806ST8O2zWz7u4hYm7m9BpiT/vw/ST5h+42k8yU9bZwHTzvCM4BXpyHqKOA/R9xl1PlrOwE3Z27f3Lf+EOBnEXE/SbDaFLh+xGOtyvx834Dbc+gTEXeQfDL5bZVwroHl4gvAoTnu7yMkHzz0Ow3YC/gfwGaM/rTazKr1BOCsTH+4lGS4+/aZbabTR6yMiMjc/g1Jn/UEYCPg1sxj/TtJNa1nvb5M0naSTk+HKt4NfJmkT5uJ7GOM06Zse/4HcHdE9Pe5WZ8CLgbOy4zMMWsEBzarlKQ/JglsPwVWk3Qw+0bEVunX49ITsacUERdHxBEkL+jfIglhAzcdsOyLwKtIAtWaiPj5kPbuAOwIDKve3QrMz9zeuW99dsjGauB+YPch+5qJWSTHYcsC9m05i4gLSD60eER6jsg5ki6V9BNJe01jfz8E7hmwfHGkgF+w/t+qmdXLzcBhmf5wq4jYNK2eTWJeOpKlZxeSES43Aw8AczOPs2VE7JvZtr/f/Kd02ZMiYkuS0wo0Yvs/kHwgC0B6Ltq2fdtk7zNOm7KmGg4JSdh9Fclwze9Lcv9ojeHAZpWQtKWkFwGnA1+OiCvTYRD/AXxC0nbpdvMkvWCM/W0s6VWSHhcRDwF3k7w4D7IK2C27IA1o64CPMXV17Zy+TymzzgCOkbR3Ovzx7/rWPzJkI32+pwIfl7STpA0lPU3SJiMefyBJL5e0p5KT1LclGZ5yeVpts2ZaBBwfEQeQnJP46bx2nA6FfA1DJs4xs1o4BThJ0hMAJG0r6YgZ7G874C2SNpL0CpJzzxZHxK3AucDH0r55g/QDo1FDHLcgOSXhTiWXkPmrvvX9/eyvgE0lvTB9/fkbYGhfN0GbRp2qkN3vQ8ArSD4wXaxkZkmz2nNgs7J9R9I9JJ+evZckWByTWf8eYBlwYTrM4gc8eo7aVF4D3Jje7zgyE4n0+Sfgb9JhFu/KLP8SyVCxL494jJHT+UfE94B/JZkyeBnQq9Q9oGQmzHsj4qbMXd4FXEkyTOMOkklXJvm/nEfy5vuedH/rgJdNsB+rAUlzSCai+bqkJSRDgXZM1708PQet/+v703iITwMXRMRPcm+8meXlZOBs4Ny037yQZPKPSV1EMkHJapKJQ/48In6XrnstsDHJJWJ+TzL5x44j9vX3wP4k54T9P+CbfevX62cj4i7g/wKfBVaSVNxWMNpYbZL0OJLw+bMp9gdAen76y0lGuHxH6UzUZnWm4YUCs26R9Frg2Ij40yHrZwG3Abunnc84+9wbuIrkk8R3kAzveHdOTbYWkbQr8N2IeGI6VOe6iBj1hmmq/R0MvCsiXtS3/H0kExG8PK3ymlnLSXod8MZh/VuTSXolSfh8ZdVtMSuKK2xmPDJ74/8lGYY2zDbA304V1iS9LB2iuTVJxew76UQpNwKfz6nJ1mIRcTewPB22hBJPnul+Jb0ReAFwlMOambXEnSQzSJq1VqWBTdLbJV2dDuf5qqRNJW0j6TxJv06/eyYfK1R6jtxvScbcf2XYdulUyZ8ZY5d/me7vepLz6P5Pev8zImLpzFtsbSPpqyTDZ/dUcnH4N5CcHP8GSb8Eria5RtG4+/sJ8HXgkHR/vfNATyGZYe7nkpZI6j/Hctj+Rl52ID1/9Ir062fZcCnpUCWXElgm6YRxn4O5jzQbR0ScO2yiMLMypHMILMl83S3pbUouCL8ys/zwzH1OTPvF6zTOXA1VDYlMT1L9KbBPRNwn6QySc4P2Ae6IiA+lnfvWEfGeShppZmZIehbJBANfiognDlj/dGBpRPxe0mHA+yPiT5TMBPcr4Hkk56tcTFLdu6bE5jeS+0gzs+ZJ+72VJOebHkMyd8FH+7bZB/gqycXcdyKZr2GPiBg2WV7lQyJnAZul5wbNJple9giSKdZJv7+0mqaZmRkMvuxA3/qfRcTv05sX8ujlAg4ClkXEDemJ/qczjUqhuY80M2uYQ4DrI+I3I7Y5Ajg9vfD8cpJJ6g4atdPKAlt6HZGPklwP41bgrog4F9g+nc61N63rwIskmplZLb0B+F768zzWvxjuinSZTcF9pJlZIx1JUj3reXN6usCpmSHs0+4bZ+XbxvGljT4CWEBywujXJQ2bhn3Q/Y8FjgXQxhsfsNH27rPMrJkevHnF6ojov4jsRJ518Kbx+zumN5/IVVc+dDXJFNc9iyJi1AQ8A0l6Nklg681EpwGbeWriMeTZR24wa+MDNtvSfWTVHt6o6haYNc/9q6rtH2H8PlLSxsBLgBPTRZ8BPkjS732Q5Fq/r2eCvrGywAY8F1geEb8FkPRNkusOrZK0Y0TcKmlH4PZBd04P1CKATXbZOea98+0lNdvMLF/L3/bOUUMnpuX3d6zjm4vnTus+e+x86/0RceBMHlfSk0iusXRY5tpOK4CdM5vNJxnWZ1PLrY+c8/id44mHvq2cVtuU7p1X9dkoZs1x9UfeUWn/CNPqIw8DLouIVQC97wCS/gP4bnpz2n1jlYHtJuCp6XTq95GM+byE5GKKRwMfSr9/u7IWmtkjZq8c9IFQs6yZ5+JOESTtQnLh3NdExK8yqy4GFkpaQHIS9pHAX1TQxCZyH9lSc1Ymn/A7uJm1zlFkhkP2PlxLb76M5Lq8AGcDX5H0cZJJRxYCvxi148oCW0RcJOlM4DJgLXA5yaeBc4Az0mmtbwJeUVUbzbquDSEta9jzcZAbLb3swMHAXEkrgPcBGwFExCnA3wGPBz4tCWBtRBwYEWslvRn4PrAhcGpEXF3BU2gc95Ht5+Bm1h7ph2vPI7m0U8+HJe1HMtzxxt66iLg6nfn3GpLX9zeNmiESqq2wERHvI+n4sx4g+STRKjRr3h8e+Xntys0rbIlVoW1BbSqjnq/DHETEUVOsfyPwxiHrFpNMR2/T5D6yG3rBDRzezJoqItaQfHCZXfaaEdufBJw07v4rDWxWvWwwm2obB7f261pQG4ercmZWFlfdzGwQB7aOGCeYjbsPB7d2cUibjIOcmRXFwc3MshzYWiaPYDbOYzi0NZ+DWjEc5MwsLw5uZgYObI1VRjAb5/Ed3JrFIa06PvZmNimf52bWbQ5sNVd1MJuKg1szOCyYmbWDq25m3ePAVhN1D2ZT8TDJ+nFIMzNrLwc3s+5wYCtZ04PZKK621YODmplZdzi4mbWfA5vlzsGtGg5qZmbd5fPczNrLga1Eba6uDeJhksVzSDMzs36uupm1iwNbSboW1npcbSuGg5qZmU3Fwc2sHRzYStDVsJbl4DZzDmlmZjYJB7firJm/bui62St8vC0fDmwFc1izmXJQMzOzPMxZuc6hLUejwto46x3obFwObFYqn9dmZmZWHYe2fEwVxvLYhwOd9TiwFcjVNZspV9fMzMzqJY+wlsfjONB1hwOblc5VNjMzs+q4yja5ssLaOBzousOBrSCuro3m0DY1V9fMrCu2WH7fwOX3LNis5JZ0h0Pb9NUprI1jWHsd5JrHga0ADmtmZmbDDQto09nOYW7mHNrG17SwNsqa+esc2hrGgc0q4yqbmVk3jBvQ8ting5zlrU1hrcehrVkc2HLm6tr0OLQN5uGQZtZkRQS0PB7bYe6xXGUbro1BzZrJgc3MzMxmpMqANh2uyg3m0PZYXQhrrrI1hwNbjlxdm4yrbGZmzdKUgDYuV+Uc2rK6ENZ6HNqawYEtJw5rlhcPh7S6kXQq8CLg9oh44oD1ewGfB/YH3hsRH+1bvyFwCbAyIl5UQpMtZ20LaNPhqly3dCmsWT4k7Ql8LbNoN+DvgHnAi4EHgeuBYyLiTkm7AkuB69LtL4yI40Y9hiO11YIDr1mtfQE4dMT6O4C3AB8dsv6tJJ2TNcQWy+9b78seq43HZs7KboeVroa1rj7vvETEdRGxX0TsBxwArAHOAs4DnhgRTwJ+BZyYudv1vftMFdag4sAmaStJZ0q6VtJSSU+TtI2k8yT9Ov2+dZVtHIfDRj58HG1Sc1aue8yX5SciLiAJZcPW3x4RFwMP9a+TNB94IfDZ4lrYTmX2kQ5ok2vbMevq62fXQ0vXn3+ODiEJY7+JiHMjYm26/EJg/qQ7rbrCdjJwTkTsBTyZ5BPYE4AfRsRC4IfpbbNO8HDIwQYFsqnCmQNcbfwL8G7Av4TpK6yPdEDLX5uOY9deMx1WLEdHAl8dsPz1wPcytxdIulzS+ZKeOdVOKzuHTdKWwLOA1wFExIPAg5KOAA5ON/si8GPgPeW3cDyuCuXLE5B0UxlvDrKP0eYT6+9atymL791nmve6da6kSzILFkXEopm2RVLvvLdLJR080/11SRF9ZJsCRV31jrHPb2sOh7VHtX0Cksn6Rxi3j5S0MfAS1h/6iKT3AmuB03o7BHaJiN9JOgD4lqR9I+LuYS2octKR3YDfAp+X9GTgUpLzHLaPiFsBIuJWSdtV2MaRHNbMxlO3T2u7Et6mYXVEHFjAfp8BvETS4cCmwJaSvhwRry7gsdomtz5ygwfWOayVrA3BrQuzRjqs2ZjG7SMPAy6LiFW9BZKOJpm065CICICIeAB4IP35UknXA3uQTM41UJX/ibNIZhT7TEQ8BfgD0xjaIelYSZdIuuThex2c2sRBuFmmGq5Yt7DWryntbKKIODEi5kfEriTDRH7ksDa23PrIhx7ya2pVmh6U2/y66LA2mI/LjBxFZjikpENJRkC8JCLWZJZvm86ejKTdgIXADaN2XGVgWwGsiIiL0ttnknROqyTtCJB+v33QnSNiUUQcGBEHbjin/CF0DhXWRU0LY9PVtueTF0lfBX4O7ClphaQ3SDpO0nHp+h0krQDeAfxNus2WVba5BXLrIzfayMPMq9T00Gbd49A2fZJmA88DvplZ/ElgC+A8SUsknZIufxZwhaRfkry2HxcRQyf2ggqHREbEbZJulrRnRFxHMqvKNenX0cCH0u/frqqNVh2fy1Y/XQkzvefZ9qFA0xERR02x/jammP0qIn5Mcr6VjcF9ZLtssfy+Rg+PbBsHEstbWkF7fN+yPxqy7TeAb0xn/1VfOPt44LT0JL0bgGNIqn5nSHoDcBPwigrbN5Cra9Y1XQlrWV04f8Nqr5F9pA3m0FYPDmvjafsEJE1TaWCLiCXAoJP4Dim5KWNzWCuPq2z10MWw1uNqm1WpiX2kjebQVi2HtelxaKsP/xbMzKbQ5dBqZvlq2jltbXn9c1izJnNgmwZX18rnY16ttnTUefCxMLO8NC20NZ3D2uR87OrBgW1MDg5WhjXzouomPMIB5bHaODOmmVmbOXBYGziwWe05LFvdOLSZ2Uy5ylY8h7V8+DhWz4FtDA4M1jUOJFNztc3MZsqhrTgOGfny8ayWA5s1gkNzeRxCpsfHy8xmwqEtfw4X1jYObFNwUDCzqbjaZmYz4dCWH4e14vjYVseBbQSHNesah46Z8fEzs0k5tM2cA0XxfIyr4cBmjeEAXSyHjXz4OJrZpOoa2vy6ZlYtB7YhHA6sKnWa2t8m4yGSZmblcuWnPD7W5XNgs0ZxkC6Gw0UxfFzNbLrqWmWrMwcIazsHtgEcCqxLHCqK5WqbmU2XQ9v4HNaq4eNeLge2Pg5r9effkZmZWXnunVe/t4tr5q9zaKiYj3956vcfaGZmZtZxrrIN56BgXePAluHKjZmZmVl9OazVi38f5XBgSzmsNUvbf19lzRTpc6vMzKwpHA6sqxzYzMxK4HBsZk1Ul/PXHNbqy7+b4tXjv7BiZVZr9t5hFXvvsKq0xzMzM2uKWdevZNb1K6tuhtWMA0H9+XdUrFlVN6BLskFtUGhbetv2ZTan8WbN+wNrV25edTPMzCxnvdC2dvd5FbekWlssv497FmxWdTMq5SBg5sBWWnVtnKpa/zYOcFYkD9EzszoZVFnLLut6eOuiJoS1WTutqboJrL1ldtVNAJLf1+wVHrxXhE4HtjqFtXHu5wD3WK6ymRVP0qnAi4DbI+KJA9YLOBk4HFgDvC4iLkvXvR14IxDAlcAxEXF/WW23ZhhnGKTDW/mqPH+t7mGtDkGtZ5K21CXktYGkPYGvZRbtBvwd8KV0+a7AjcArI+L36X1OBN4APAy8JSK+P+oxHIMbpHf+m8+D64ayZoq08jS4qvkF4NAR6w8DFqZfxwKfAZA0D3gLcGAa9DYEjiy0pdYJvXPdfL5bOzmsFa+o51D3310RIuK6iNgvIvYDDiD54PIs4ATghxGxEPhhehtJ+5D0hfuS9K2flrThqMfobIWt7tW1SfbtCpyZFSEiLpC064hNjgC+FBEBXChpK0k7putmAZtJegiYDdxSbGutaWYautpeeevaeWx1fsPfhqBWho4PjTwEuD4ifiPpCODgdPkXgR8D7yHpM0+PiAeA5ZKWAQcBPx+208oDW5ooLwFWRsSLJG3DkPJh05RdBetqgPOwyOlrcKVnpC2W3/fIz116g5N199pN+dHqvaZ5rx/OlXRJZsGiiFg0jR3MA27O3F4BzIuISyR9FLgJuA84NyLOnWbjOq3NfSTMPKyN2l8bw5tVp41hbdZOazo1NHKy/hGm2UceCXw1/Xn7iLgVICJulbRdunwecGHmPivSZUPVIf6+FViauT2wfJinMqprdRiy2D+Esg5tMivCFsvvWy+s9ZbZ2FZHxIGZr+mENQANWBaStib5JHEBsBOwuaRXz7SxHVN6H9kWHjY5M1Wcv1bH6tqsnda0Mqz1eGjkWMbqIyVtDLwE+PoU+xvYZ466Q6WBTdJ84IXAZzOLjyApG5J+f2mej9mVsDaMA5y1yaCgZpVYAeycuT2fZOjjc4HlEfHbiHgI+Cbw9Ara10hV9JFlKjNIObzVXx3f4Lc5qGV15XmW4DDgsojovcFe1Ts9IP1+e7p8WJ85VNUVtn8B3g1k/0vXKx8C2w24X201LQS1JcCVefFzq1YvpI0T1OoY5lo6HPVs4LVKPBW4K339vgl4qqTZ6UySh7B+tchG+xda1kfWgYObTaXtVbVBini+dQzhBTuKR4dDQtI3Hp3+fDTw7czyIyVtImkByYRdvxi148rOYZPUmyL6UkkHT3D/Y0lmI2PDrbce6z5Fv6lvcuDp6T2Hrpz/Zs0xafiq2wn7VU6TPSlJXyU5cXqupBXA+4CNACLiFGAxyZT+y0hmxzomXXeRpDOBy4C1wOXAdIdbdlKefeSmGz8u38bloA6Byee7DVf261Sd3th3Lahlde2ctjxJmg08D/jLzOIPAWdIegPJB5ivAIiIqyWdAVxD0je+KSIeHrX/KicdeQbwEkmHA5sCW0r6Mmn5MD05L1s+XE86fnQRwCa77Oz5z3Pm4FYPa+YFs1cOGuo8uSZWeGZaKatbaGuaiDhqivUBvGnIuveRBDybntz6yC3n1O8aIWt3n1eL0NbTa4uDW/nqEta6HNSK1JUZIyNiDfD4vmW/IxlZMmj7k4CTxt1/ZUcwIk6MiPkRsSvJjCo/iohXM7x8OCOurk2mSc/LwyLbqY7DGs2KVnYfaQkPlyyXw1r9+FjUUx0j74eA50n6NUlp8UMVt2dKTQo1k2j787N68oQiZgM1ro9soi4HtyYO254JB5THyvuY1CWYN1nl12EDiIgfk1xMbmT5cFJFVl66Emb23mGVh0daaYoIah4WaU1VdB9ZlboNixzEQyWLU/WbeAe10fI+n60rQyOL4iNnY+tKOLVqtbmq1rVPrs3aoisVt7JeoxzWzKan9e8eXF3rFp/H1lweAmlmddeV4FakKsNaF6frnwkPjayPVgc2h7X8dfV5W7HKCmoOhGb10tShhg5tzeOgNhkft3podWArStdDS92fv6tszeIQZWZN1LbQVsZwyKoqLA4dM5Pn8XOVbTK1mHSkCH7T3m2z5v2BtSs3r7oZuSjiWmx14KBmZtYdVbxRd1CztnCFbZrqXl0qi49DM9XlotlVhrWqHtsTjpgN1tRhkWUpY3bbNr4+Oazly1W2arXvP5TiqmsOKeur+/FwlbWeXFkzM+uWMt+ge2KR4vi4VqeVga0IdQ8nVfFxsXHVaRbIurTDzKztyg5rVqy8jrGrbNPTunPYiqiqOJQ0V5vOZWsyByQzG6UJF9FuqyKHQ3YhrE3nPeLS27YvsCXlyfui2ja11gU2K9/eO6xqzYuQ5c9hzczaatb1K30OXsWqrKpN9wP9MgoATXo/tmb+Omav8GC/cbQqsLm6Vp06hzZX2apR96C2xfL7SjmZH9p5Qr+ZNV/Tq2tNqKqVbdy2zfQ9m6ts5fK7iBGK+Id8ztxrec7ca3Pfr1md1D2smVn9uFLVHg5r3ZDH78Dnso2nNYGtCTMCZoOaQ1u5mvD3McqaeVF1E8bmsGZmk3JoW19ZowCaxmFt5tr0XLqgFYFNG+Wfzsv4Q25bta3u//xND21N4LBmZjO1dvd5Dm4lKGo4ZJEVk6qm7N97h1W1f48ziTyek6ts5WhFYGuCUcGsbcHNzB7l89fMJtMLbg5vBq6qFaXtz68t/E5igLz/eMcNY20IbXX/x3eVrTiurplZURzc8tW06prDWr25ylY8B7Y+Vf9zutpmTeSwZmZl6FLVrWnnrzmsNVeXnmtTObAVbNLw1eTQVvd/fFfZ8uWwNpyHQ5oVpyvBrct8vlp5Zvqcq7weXhf43URGVUMhR92/ycHN2s9hrRsknSrpdklXDVkvSf8qaZmkKyTtny7fWdJ/SVoq6WpJby235dYFVVbdmhgYi/ggqYjqWlVhzSY3099Zk4dFStpK0pmSrk37vKdJ+pqkJenXjZKWpNvuKum+zLpTptp/qy6cPRN1C2v9+/rR6r1y25/5Ytpm0/QF4JPAl4asPwxYmH79CfCZ9Pta4J0RcZmkLYBLJZ0XEdcU32TromyAmnX9ygpbYjPhsFaNvXdYNeMLanfYycA5EfHnkjYGZkfE/+qtlPQx4K7M9tdHxH7j7twVtoZoWrWtCS98Hho5M66udUdEXADcMWKTI4AvReJCYCtJO0bErRFxWbqPe4ClQPNKEtZITR4yWdT5a02orpUd1ro6BHKYqodGNrHKJmlL4FnA5wAi4sGIuDOzXsArga9O+hgObNS7ujZo300KbtZODmtT69j5a/OAmzO3V9AXzCTtCjwFuKi8Zpl1a6KSsrUhrJnlYDfgt8DnJV0u6bOSssO4ngmsiohfZ5YtSLc9X9Izp3oAD4nMWVlhqgnDJJtQWm/S0Mg184LZK1V1M1oV1rZYfl/jZmKbyv0PbTTJ/91cSZdkbi+KiEXTuP+gP8x4ZKU0B/gG8LaIuHu6jTPLi4dMWo/D2nAzff82a6c1rL1l9sT3XzN/HbNX5P+h54T9I0zdR84C9geOj4iLJJ0MnAD8bbr+KNavrt0K7BIRv5N0APAtSfuO6h8rC2ySdiY5H2IHYB3Jkz9Z0jbA14BdgRuBV0bE74tqR5P/YXvhsO7BzdqjTWHN1rM6Ig6cwf1XADtnbs8HbgGQtBFJWDstIr45g8folDz7yHWbdKraO7ZeeOtKcMu76t/U6lqT3/eVqQkfupdoqj5yBbAiInojSM4kCWxImgW8HDigt3FEPAA8kP58qaTrgT2AbChcT5Wv4r2T0fcGngq8SdI+JE/whxGxEPhhersQTRoKWcfHHUcTXhi7ci7bnJXNGxfeVB0bDglwNvDadLbIpwJ3RcSt6bj9zwFLI+Lj1TaxcXLtI9tWSc6Th0xWz2GtnmZyvLo0zX9E3AbcLGnPdNEhQG9yrecC10bEit72kraVtGH6824kE3bdMOoxKquwRcStJCVBIuIeSb2T0Y8ADk43+yLwY+A9FTRxWqoOTa62zUyThkZWxdW17pL0VZLX5bmSVgDvAzYCiIhTgMXA4cAyYA1wTHrXZwCvAa7sTWcM/HVELC6t8Q1VRB95z4LN/H88hToMmSwiXNe5uuaw1l4zGRpZ1LDIAh0PnJbOEHkDj/aDR/LYyUaeBXxA0lrgYeC4iBg1sVc9zmHrOxl9+7SjIv2EdrsiHjPPf9yqw1pWE85tM7NmiYijplgfwJsGLP8pg89vs2nIs490aBtf14ZMjquJs/g5rE3OQyPHExFLgMcMm4yI1w1Y9g2SUwXGVnl0nfRkdEnHSrpE0iUP392NIW3jqlOAbJKuDI00s+bIo4986P71X9s8PHJ6PFyyOGVU1xzWqjWT33ETPxwoSqWBbcjJ6Ksk7Ziu3xG4fdB9I2JRRBwYEQduuOX0hrL5n7dcTTreDm02Ux08f80KklcfudGmHu6dhzKCW92HQzbtDXST3n/UmY9j9aqcJXLYyehnA0cDH0q/f7uC5o2trtUsD400M2uuovtID42cnKtt+Si6ulZWyBj1PtDvwywvVZ7DNvBkdJJO6AxJbwBuAl6R54N26VMCh7bJeAISM6uBwvtIh7b262rFvw5hbZz1g9T1fdtMzmWb6XXZrNpZIkedjH5ImW2ZVF2razZzDm3r85s6s3KV1Uc6tNm4mjIzZF3CWpn7rWvIs/zUYpbIsnSputbjKpuZmY3i0FYPXZkMxmEtf9Ntz6TvC11lq043a+U5qNs/q+XPE5CYWVd0JSx0SdcmG+lqWJtEG55D10xZYZP0HSD6Ft8FXAL8e0TcX0TD8tbF6lqPq2yT89BIm46uni/SZW3pI8GVNiteUdU1h7Xy+Lps1Rjn3cUNwL3Af6RfdwOrgD3S253jf9jpa3JgrlOlbc28/veFZlYx95E2Y3lXOLtUXXNYm0wVz6eM6+611TjnsD0lIp6Vuf0dSRdExLMkXV1UwyxfrrKZmRWiVX2kq2xWlCa/WW9bWJspV9nKN85HMNtK2qV3I/15bnrzwUJalbM8P33xP2031anKZma10vg+sp/PZ2s2V9dsXK6yNcc4/9XvBH4q6b8k/Rj4CfBXkjYHvlhk4yxfDpsz49BmZgO0so90aCtPF451EW/SPRSyWg7L5ZpySGRELJa0ENiL5Jow12ZOov6XAtuWC1fXzKwMnnCkm5reR47i4ZHdlld1zWGt3qo4ZcZT/E/flO8wJM0G/gp4c0QsAXaW9KKiG2bF6MKLT5FcZTOzrLb3kV2o/rRJFz44cljL36TP1VW28ozzn/15knH4T0tvrwD+obAW5cjVNSuCQ5uZZTS2jzQbpq7VNYc166pxAtvuEfFh4CGAiLiPZNiHNVRVL0Rt+iTGoc3MUq3vI11lK06ex7bt1TWHtWKVXWXz5CPTM85/94OSNiO9MKik3YEHCm2VWc35YtpmlupEH+nQ1h11rK45rFnXjXMdtvcB55CMyz8NeAbwuiIbZe3Ulmt2dC2sedIBs5E600d6EhJrM4c1q7NxZok8T9JlwFNJhnm8NSJWF96yGWrT8Durj66FNTMbral9pFWvzcMhm1hds8lnjJz0QtqeLXJ8QwObpP37Ft2aft9F0i4RcVlxzaoXf+piDmpmltXVPtJVtnar+4Wyi+L3eVZ3oypsH0u/bwocCPyS5NPDJwEXAX9abNPM6sFhzcwGcB9ptVC36lqeXF0zSwz9L4+IZ0fEs4HfAPtHxIERcQDwFGBZWQ00q5LDmllC0qGSrpO0TNIJA9ZvLeksSVdI+oWkJ2bWbSXpTEnXSloq6Wn992+aLveRnoBk5up4DOs22YgnGqmGZ4uczKB+TtL7Ja2UtCT9Ojyz/Ylpf3qdpBdMtf9xPpbZKyKu7N2IiKuA/SZ6Nsbhc66pugk2hrUrN3dYM0tJ2hD4FHAYsA9wlKR9+jb7a2BJRDwJeC1wcmbdycA5EbEX8GRgafGtLk0n+8g6Bg6z6XJYsxwN6+c+ERH7pV+LAdL+80hgX+BQ4NNpPzvUOIFtqaTPSjpY0p9J+g/a1dmarcdBzewxDgKWRcQNEfEgcDpwRN82+wA/BIiIa4FdJW0vaUvgWcDn0nUPRsSdpbW8eO4jbVrqONlIV6trNpiD7PRM0M8dAZweEQ9ExHKSURkHjXqMcf7TjwGuBt4KvA24Jl1mE+pila0JU/q7qmY21Dzg5sztFemyrF8CLweQdBDwBGA+sBvwW+Dzki5Pw02b/tE620e6ymZN5lCSvw4PixzVz705PVXgVElbp8vG6VPXM860/vcDn0i/rAUmmbK17RzUrC3ioQ0m+XueK+mSzO1FEbEoc1uDHqrv9oeAkyUtAa4ELgfWAhsB+wPHR8RFkk4GTgD+drqNrKOu95GeNbI6dZtsxNU1q7t4aINJLyMwVR85i8H93CeBD5L0lx8kmazq9YzXp65nnAtnN47/2W06HNbMWB0RB45YvwLYOXN7PnBLdoOIuJu0siRJwPL0azawIiIuSjc9k6QjM+ucOlYluziVv6trU5v0mmwtNU4f+Zh+LiIeCSTpcPnvZrYf2af2q9fHM2Yl8hBIs7FdDCyUtEDSxiQnS5+d3SCdIWvj9OYbgQsi4u6IuA24WdKe6bpDSIYNWkvUMYS0Xd2qa03isFasLg6LHNbPSdoxs9nLgKvSn88GjpS0iaQFwELgF6MeY+h/fDrd5FMmbv0MTTWFdFn8j91ODmpm44uItcCbge+TTKhxRkRcLek4Scelm+0NXC3pWpLZJN+a2cXxwGmSriCZQfEfS2t8QarsI+vSP1rzebIRs9wM6uc+LOnKdNmzgbcDRMTVwBkkH16eA7wpIh4etfNRQyKXA2+V9GSSk8m/B5wbEb+f2fOZWmYK6eeRlA0vlnR2RLTmU9nD51zD4nv7Z8W2MjismU1fOh3x4r5lp2R+/jnJp4SD7ruE5OLSbVJJH1nX/tHnsk3Nlcjq+UP46Zl0WOTeO6xqxGRzeRrSz71mxPYnASeNu/+hgS0iTieZupn0U8RDgW+mncUPSK41MLJ8NwOPTCGdPn5vCulGB7Yuzg5ZJw5q09eFN2BbLL/Pb6Rs2irsI2vbPzq0lcNT+VtbzdppzaSTgrTeWP/1EXF5RPxTRDwbeBHJFMZvLLBd057u0uqrDp+yOKyZWVFK7iPdP5pNwNU1a7Jpf0yTnkT+jYg4togGpaac7lLSsZIukXTJw3f/ocCmWJN5YpHEnJXdmwXMrAol9JFjTQed7SMfur+8PtLV6sHyOi5tra5ZfTno1kNdpxkaZwrpRRFxYEQcuOGWzXhDnj1nrarz19o+RWsvoDmomVlLjTUddLaP3GhTvxZavRU9HNKho3xdnC2ySHUNbFNOIT1KnkPw2h5wilaH4ZBtMnvloA/XzaxDZtQ/lsFVtvXVrbpmZs0z5YWz0wugvgrYLSI+IGkXYIcCJxwhItZK6k0hvSFwajoFZuNVOTNk28OnK2pWpTkr1/kNVQeV3Ue2uX+0cng4pFnzjPPu4tPA04Cj0tv3kEwpXKiIWBwRe0TE7unUl2ZmZnVTeh/ZhP7RVbaEq2tT83DI+vMxrN44rwB/EhFvAu4HSK8xs3GhrcqBh0Wur4rn4OGQ1hSeitxmoJF9pHVTXtU1s3H4PLb8jBPYHkqvKxMAkrYF/B9vteLhkGZWEfeRQ3S9ytbm5+831GblGiew/StwFrCdpJOAnwL/WGiraqjJVTZX12wSrjpNny+f0EnuI61QdZvKP08eDmk2niknHYmI0yRdChxCcv2Xl0bE0sJbloOlt21f+IuBWVkmnSHSIcKsOE3uI8twz4LN/OFPy7i61k3PmXtto4sXTTc0sEnaJnPzduCr2XURcUeRDbN8dKG65uGQZlY295E2St0mG6ljdc26Ye8dVk30vnDWTmtYe8vsAlrUTKMqbJeSjMkXsAvw+/TnrYCbgAVFN65ufrR6L5fXzVpqi+X3tfqcE8ud+8gxucpmVfD7NWuToR/dRMSCiNiN5FovL46IuRHxeOBFwDfLauBMdflcKpeuzcrnIajd0JY+0vLX5upansMhfcqK2fimPIcN+OOIOK53IyK+J+mDBbbJGszDIc2sY9xHjsFVNmuzw+dcU3UTpmXxvfuU+ngzGRZpiXE+vlkt6W8k7SrpCZLeC/yu6IbVVVOqVk1pp5lZw7mPtEe4utY9TQtrMHmbPcy0OuO8IhwFbEsybfG3gO3SZY3R5WGRZXJ1zSzhYZGd0vg+siw+R9R6PJ2/2fSMM63/HcBbS2iL5cTVtfbxlP7l8MQjNl3uI62nzdU1M6vWlIFN0n+RzIS1noh4TiEtKkie12TzbJGP5Spmu/hcE7PxtKWPLIvPZWsmD4e0mZr0PDZLjDPpyLsyP28K/E9gbTHNsZnqSnXNwyGt7uasXJfbJ+VWa+4jp6mNoc3VNTMr0jhDIi/tW/Tfks4vqD2FanOVrcqg5k9MrE08LNKmo019pNkgeVfXPJ1/sz1n7rWdKQ7UyZQf5UjaJvM1V9ILgB1KaFsh8gwXdfmD7VpY61p1bdLz16x6bTqHUNKhkq6TtEzSCQPWby3pLElXSPqFpCeOe98ma1sfWZY2fSji6lq91OnDdOsOSVtJOlPStZKWSnqapI+kt69I+8et0m13lXSfpCXp1ylT7X+cIZGXkozPF8kwj+XAG2bwnCxHXQtrNr42hQWrlqQNgU8BzwNWABdLOjsisnND/zWwJCJeJmmvdPtDxrxvk7mP7LA2Bc9BfO6a2dhOBs6JiD+XtDEwGzgPODEi1kr6Z+BE4D3p9tdHxH7j7nycj3P2jojdImJBRCyMiOcDF0/vOdRLW6psdanwlalr1TWrRp7n17QkOB8ELIuIGyLiQeB04Ii+bfYBfggQEdcCu0rafsz7Nlnr+siytD3sTEdXqmseDtltbf39S9oSeBbwOYCIeDAi7oyIcyOid07zhcD8SR9jnFeInw1Y9vNJH7Aumh7aqg5rHgpZjiqGQ7ZtMgDLxTzg5sztFemyrF8CLweQdBDwBJLOaZz7Nlkr+0gb7p4Fmz3ylYe6Tk7k6tp4Ft+7T9VNsOrtBvwW+LykyyV9VlL/m9bXA9/L3F6Qbnu+pGdO9QBDh0RK2oGkU91M0lNIhnsAbElS5rOMMich6WJYM2uyMmeM3ODBiYL+XEmXZG4viohFmduDdtg/lf2HgJMlLQGuBC4nGSI4zn0bx31kPuo8Y2TTKoB1r65ZvRw+55rOhc0NHoTZKybqi6fqI2cB+wPHR8RFkk4GTgD+FkDSe0n6w9PS7W8FdomI30k6APiWpH0j4u5hDRh1DtsLgNeRfEL68czye0jOVWi8PGeNLEvVYa0qXayuzURLhuFVqmOzRa6OiANHrF8B7Jy5PR+4JbtB2tEcAyBJJOdyLScJLyPv21Ct7yPLUofQVtX/uqtrZo0wTh+5IiIuSm+fSRLYkHQ08CLgkIgIgIh4AHgg/flSSdcDewCX9O+4Z2hgi4gvAl+U9D8j4hvjP6dmadJU/3UIa66umU2m4ddluxhYKGkBsBI4EviL7Abp7Fdr0vPU3ghcEBF3S5ryvk3UlT6ybdr6IUwTqmtN+4B8Ohbfuw+Hz2nLPEpT89T+64uI2yTdLGnPiLgOOAS4RtKhJJOM/FlEPPIpiKRtgTsi4mFJuwELgRtGPcaoIZGvjogvk5w4/o4Bjfv4gLs1UhNCWx3+MaoKa12trnk6f6uLdIarNwPfBzYETo2IqyUdl64/Bdgb+JKkh4FrSGdKHHbfKp5HnrrUR5ahiCpb3cNZXT/AcXVtMl0LbZPae4dVbf3w/3jgtHSGyBtIRpxcDGwCnJcMPOHCiDiOZIKSD0haCzwMHBcRd4za+aghkb13yXMGrGv8+QdN0uWwZlalvIdFNrnKFhGLgcV9y07J/Pxzkk8Jx7pvC7iPrIm6B7OiNaG6Vpaqr8Hm0NZdEbEE6B82+UdDtv0GMK2RGaOGRP57+uMPIuK/s+skPWM6D9JP0keAFwMPAtcDx0TEnem6E0k+mX0YeEtEfH8mjzWuulbZ6hDWqtTV6lpVqj6PxKwputZHlmGqKlubgllTP7ixqTm0WRHGecX4tzGXTcd5wBMj4knAr0guJIekfUjOb9gXOBT4dHrR1VLUbar/uoQ1V9fKN5PhkJ5wJF95h1j/flqnM31kGXqhLDt1fp5T6NdBnmEt7+paUcMh23z+2iBdm33RijfqHLanAU8Htu0bn78lyXkIE4uIczM3LwT+PP35COD0dPaU5ZKWkVx0tXPXtHFYc3XNzOrLfWRx2hTOiuShkPXmSttwLT6PrTCjPubZmGRs/ixgi8zX3TzaeeQheyG5iS6wGg/l92lVHapsDmsOa9ZurrK1QmP6SKuPOg+F9GQj+XOlzfIy6hy284HzJX0hIn4z3R1L+gGww4BV742Ib6fb9F9IbuwLrEo6FjgWYMOtt55u80aq6ny2ugQ18DDIKnl2yPrp2DXZbAxN6iM3np1vH2mT6eJQSOjecMh+ba20eWr/co2aJbJnTXoC9L7Apr2FEfGcUXeKiOeOWj/oQnKMcXHWzP4XAYsANtll51i7cnNmzfvDFE9lfGWHtjr90Vcd1lxdm5wrN83R5BkjbT217yPnPH5nz1ppVqG2hjYrzzjvFk4DrgUWAH8P3EhyXYGJZS4k95LsheSAs4EjJW2SXmR1IfCLcfdb5zf6wwLZj1bv5bCWUeffYdu1ZYbIWdevLGS/RRwfB+xWaEwfadXpanXNHrX43n08RNImNs4ryOMj4nPAQxFxfkS8HnjqDB/3kyRj/c+TtETSKQDpxVTPILno6jnAmyLi4ensOM83/HmHl/5gVregVnVYMw+H7CKHtsarfR/58EYzbI3NSJfDWteHQw7i0Jbw38b0jDMk8qH0+62SXkgy/GL+TB40IgZeSC5ddxJw0kz2n+fwyDyHRkK9QlpPXYKaq2sz4zf+xVXXenwumw3QiD7y3nkb+DWiAh72bIP0QpuHSdq4xnkl+QdJjwPeCbwL+CzwtiIbVTd1CTRFqMtzc1hzdS1PRQc3s4zG9JEOD83WtOqaTc3VNhvXlBW2iPhu+uNdwLMBJL2twDblos6TkNRBXYIaOKxZPhzSrApN7SOteHUeClmGNr1nKpKrbTaOSV9N3jH1JtXLOwjUKeTMRFueh1lPmWGtLZOzWKFq20e6ylaOuh9nV9fqp4vVNof68U36itKYsVsObeurW/tdXUvMdDjkTM9NcQgxy1Wt+8i6h4mmy/v4eihkd3gmSRtm0leVRl3TxaGtnrNAOqxZHqoYCumAa1OofR/p0NYMTRwKCa6czJRDm/Ub+oot6R5Jdw/4ugfYqcQ22gzVLaiBw1qWJxsphs9psyK1oY90aMtf3Y+pq2vN4dBmWUNfWSJii4jYcsDXFhExzuUAaqWrVbamtNMm1+WpuqsMZa6ydVvb+kibOQ+FtLzVfYjkc+ZeO+N9uBo7nnp/FJSzLoW2Og6B7HF1zfLgCprZzNW9ItQUdQ9rZfIb8PzVObRZOTr3St2F0FbHNvU4rK3PwyGbzVU2awOHtvZzda356l5ts2L5VToHdQpIdWqLFS+P4ZBNDB2urpnly6FtcnWvrjmstYtDWzd18hW6iCpP1UGpzkMge1xdSypq2S8rVhnBromB18zy4aC7Pg+HLEce1TYHv2bp7CtNm0Jb3YMadDesOaDlz9U1s2LcO28DB5BpKOJYubo2uTwmwGiatoQuB/2pdfqVuemhrQlVNehmWHNA65YuVNkkHSrpOknLJJ0wYP3jJH1H0i8lXS3pmL71G0q6XNJ3y2u1TcLBbbSijk/Tw5rfdFfD57bVg6StJJ0p6VpJSyU9TdI2ks6T9Ov0+9aZ7U9M+9PrJL1gqv13/hW5aaGtF9KaENSge2GtzKDW5en8rVySNgQ+BRwG7AMcJan/HcKbgGsi4snAwcDHJG2cWf9WYGkJzbWcOLQ9qhfSijomTZ4V0urBoa1yJwPnRMRewJNJ+rsTgB9GxELgh+lt0v7zSGBf4FDg02k/O5RfjWlGaGtSSOsiV9Ss5VW2g4BlEXFDRDwInA4c0bdNAFtIEjAHuANYCyBpPvBC4LPlNdny0PVqWxnPv4iw5upaNzU5tDX5b0jSlsCzgM8BRMSDEXEnST/5xXSzLwIvTX8+Ajg9Ih6IiOXAMpJ+dqjuvgr3qWtoa3JQ60J1raqglld1rYkhY+3u86puQhfNA27O3F6RLsv6JLA3cAtwJfDWiOj9of4L8G7AZYSG6lJoK7qaluWwZnnzEMlK7Ab8Fvh8OvT/s5I2B7aPiFsB0u/bpduP06euZ1b+bbaspbdtP9GLWVNDWk/bw5qraTbIFsvv454Fm1Xahg0fmijQz5V0Seb2oohYlLk96A8++m6/AFgCPAfYHThP0k9IPnW8PSIulXTwdBtm9dELMG0djl12KPUwSCvS4nv34fA511TdjFqZsH+EqfvIWcD+wPERcZGkk0mHPw4xTp+6Hge2jLUrN2fWvD/kvt/phLamBzVod1hzULOWWh0RB45YvwLYOXN7PkklLesY4EMREcAyScuBvYBnAC+RdDiwKbClpC9HxKvza76VqW3BrYrqYVFhzdU1y2paaNt7h1V1fR88Th+5IiIuSm+fSRLYVknaMSJulbQjcHtm+6n61PV0Z4zDmIoKG6P+AJs2kcgobQ1rdTpHrS1vktqqicNMx3AxsFDSgnQikSOBs/u2uQk4BEDS9sCewA0RcWJEzI+IXdP7/chhrR2aPEyyzGGPZenSFP42vkFDJPMcMtnFyyn0i4jbgJsl7ZkuOgS4hqSfPDpddjTw7fTns4EjJW0iaQGwEPjFqMdwhW2AsiptbQhoWW0Ma3UJaWZVioi1kt4MfB/YEDg1Iq6WdFy6/hTgg8AXJF1JMtzjPRGxurJGWymaVm2rQ0Bry1BIV9eapWnVtgY6Hjgt/VDzBpJRJxsAZ0h6A8mHmq8ASPvPM0hC3VrgTRHx8KidO7ANUWRoa6O2hTUHNbP1RcRiYHHfslMyP98CPH+KffwY+HEBzbOK3Ttvg9qGtjqEtB4PhSyGqzzj8WQkxYmIJcCgYZOHDNn+JOCkcfdfn1exGmpbCLGp1Wno4yB1fUNk62vpsEizkeo2xLBu7WlLWDMrUt0+DKiLVlTYNniwuH0XVWlrkzYE2zqHtKI4VJhZEaoeJlmnkNbTlmGQ4DfUZlWo9FVN0rskhaS5mWUnSlom6TpJLxh3X118w10HTQ5rvWqa/3asCA7ENlN59pFVKDM41XkSkSLDmqtrZt1Q2SubpJ2B55GchNdbtg/JLGL7AocCn5a0YTUtfFSTQ0mRmnpcmhrSPBzyUb54trVdk/rIUYoOUXUNaT1tC2uurlkZ/Hf2WFW+yn0CeDfrXyjuCOD0iHggIpYDy4CDxt1hkW/CmxpOitLE49HUoGbN5SqbzUDufWSV8gxVda6mZbVpGCT4TbRZlSp5tZP0EmBlRPyyb9U84ObM7RXpsrE5tBWvacfBQc3MmqTIPrJKMw1ZTQhpZfFQSLNuKWzSEUk/AHYYsOq9wF8zePrnQe+qY8AyJB0LHAuw0ZZbT9jK6fMkJM3RppDm4ZCPtXb3ecy6fuWU21Rti+X3cc+CzapuhtVMW/vIcWRD11SvbU0MaEVX1jwUcn2e0r+d9t5hVWsvhTWJwgJbRDx30HJJ/wNYAPxSEsB84DJJB5F8WrhzZvP5wC1D9r8IWASw2Q47r9dhzV4p1swb2IflosuhrWnVNWu3XiCbKriZ1U2VfWSdDJtRsolBDRzWzKwYpb8iRsSVEbFdROwaEbuSdED7R8RtwNnAkZI2kbQAWAj8YpLHKbq60sXg0qTn3Kbqmk1tUCWtDtU1s+kqq4+sm+x5aQ5rgzmsPZara9YVtXpVjIirgTOAa4BzgDdFxMPVtmq4JgWYmWrSc21bWPNwyPE4oFnbNa2P7BKHNeuyooKz/wYfVfmFs9NPELO3TwJOymPfRQ+NhG4Mj3RYs6ZwaLO2KbKPtHw4rFXD1TXrklpV2IpQxhv4JgWa6Wrzc2sCV9fawdP7m7WTw5qZlaHyCltbtKXS1uSA5uqamZmVoYxrrDmsDefqmnVNJwJbGUMjoZmhrckBzczMrCxlXgjbYc3MsjoR2MChraetAa2N1TUPhzQzq16ZQQ0c1szssToT2MpUp9DW1oBmZmZWlLJDWo/D2tQ8HNK6qFOBrawqG1QX2roY0NpYXTMzy8u6jR8NILNXtH6usRmpKqiBw5qZDdepwAbtC21dDGhd4OGQZlaEbCBxeEtUGdJ6HNbG4+qadZVfrQuWd6Bau3Lz9b66ztU1awpP7W91s2b+ulqElarU5fk7rFkbOEwXq3MVNii3ygYzq7Q5lJmZWZG6VHWrQ0DLclgbnwOBdVknAxvUO7Q5pI2nrdU1D4c0s6q0NbzVLaiBw5qZja+zga0Kw0KbA5qZmdVN08NbHUNaj8Pa9Li6Zl3X6cBWdpUNHM7y0tbqmrXXPQs2q7oJMyLpUOBkYEPgsxHxob71fwW8Kr05C9gb2DYi7pD0duCNQABXAsdExP2lNd5mrD/81DXA1Tmk9TisTY/DmjWBpBuBe4CHgbURcaCkrwF7pptsBdwZEftJ2hVYClyXrrswIo4btf9OBzaoJrSZDVPmcMh7FmzmiTBsLJI2BD4FPA9YAVws6eyIuKa3TUR8BPhIuv2LgbenYW0e8BZgn4i4T9IZwJHAF0p+GpajOgW4JoQ0qCaogcOaWYmeHRGrezci4n/1fpb0MeCuzLbXR8R+4+6484HNzKxoTa+uAQcByyLiBgBJpwNHANcM2f4o4KuZ27OAzSQ9BMwGbimwrVaBsgNcU0Jaj8OadcFz5l7Lj1bvVXUzakeSgFcCz5l0Hw5suMrWNB4OmR9X2YpXdljb4IF1k/xO50q6JHN7UUQsytyeB9ycub0C+JNBO5I0GzgUeDNARKyU9FHgJuA+4NyIOHe6DbRmKSLANS2k9TisTcbVNcvbhP0jTN1HQjLk/1xJAfx73/pnAqsi4teZZQskXQ7cDfxNRPxkVAMc2FIObVY1zw7ZPg2qrK2OiANHrB/0KcmwF8wXA/8dEXcASNqapBq3ALgT+LqkV0fEl2fQXmuYQWFrnBDX1JAG1QU1cFgzy9lUfSTAMyLiFknbAedJujYiLkjX9Y86uRXYJSJ+J+kA4FuS9o2Iu4ft3IHNGsXVNbNKrAB2ztyez/BhjUeyfsf0XGB5RPwWQNI3gacDDmwdN6wK1+SQ1uOwNjmHNWuiiLgl/X67pLNITiW4QNIs4OXAAZltHwAeSH++VNL1wB7AJY/Zcaqe0zxVxGHAuqhBVaBGadlxvRhYKGmBpI1JQtnZ/RtJehzwZ8C3M4tvAp4qaXY6jv8QktmxzNazZv66xoe1WTutcVibAYe15uvi71DS5pK26P0MPB+4Kl39XODaiFiR2X7bdDIvJO0GLARuGPUYrrCZmeWsZWGNiFgr6c3A90mm9T81Iq6WdFy6/pR005eRnKP2h8x9L5J0JnAZsBa4HOgf+2/WeFUGNXBYM6vQ9sBZyWeSzAK+EhHnpOv6R50APAv4gKS1JJcBOK53GsEwDmx9fC5bfbW5Alr1+WuefCQ/bQtrPRGxGFjct+yUvttfYMB0/RHxPuB9BTbPrFIOa2bdlc6g/OQh6143YNk3gG9M5zEc2AZwaDOzSbQ1rJnZYFUHNWhHWHN1zWw0n8NmjdDm6pqZmTVL1eeq9TisWR35d5o/B7YhHBCsLFUPh+xxdWhmfPzM2q0X0uoQ1MBhzaxLKgtsko6XdJ2kqyV9OLP8REnL0nUvqKp9ZmbjclizvLmPrIe6hbQehzWzbqnkHDZJzya5kOqTIuKB9CJzSNqHZDaVfYGdgB9I2iMiHq6inT6XrR5c7SyPJx+ZPoc1y1tT+si2qls46+ewZtY9VVXY/g/wofTCcUTE7enyI4DTI+KBiFgOLCO58FxlHBbMzKxkjekj26KulbR+bQhr1g0O5fmqKrDtATxT0kWSzpf0x+nyecDNme1WpMuso9oemOty/lqWK0bj87GyguTaR2qjdY0II2VrSkjraUtY8xt5s+krbEikpB8AOwxY9d70cbcGngr8MXBGeqXvQe/OB45JlHQscCzARltunUeTh/LQSOsaD40cziHN8lBmHzlr7uMeWT4snKy9Zfb4jW+opgSzfm0JauCwZjapwgJbRDx32DpJ/wf4ZkQE8AtJ64C5JJ8W7pzZdD5wy5D9LwIWAWy2w86FpymHtvK1vbpWdw5tj3JIs7yV2UduuvvUnVcbg1xTA1pPm4IaOKx10XPmXsuPVu9VdTNaoaohkd8CngMgaQ9gY2A1cDZwpKRNJC0AFgK/qKiNZoW7d169r6zR9aByz4LNOn8MrBLfoiZ9ZHbYYBOGDzalnVNxWDOzrEpmiQROBU6VdBXwIHB0+kni1ZLOAK4B1gJvqtPsV66ylcfVtfroWqXNAc1qoNZ9ZJ2qcU0PZv0c1KxtXGXLRyWBLSIeBF49ZN1JwEnltsjMRsmGmLaGNwc1q4um9pFlBLm2BbSetgU1cFgzy1NVFbbGcpWteK6u1dtUwaZJgc4hzax4g0LWuCGurQGtp41BDRzWbH2uss2cA9sEHNrMhmtCoHNQM6tW24PYOBzWzGxcDmxWK66utV9Vgc4hzczqoK1BDRzWbDhX2WbGgW1CrrLlz2HNIP9A56BmZnXQ5qAGDmtmRXJgmwGHNsvDvfM2YM7KdVU3ozHGCXQOaWZWF20PauCwZuNxlW1y9b4IlHWGq2uWF4c1M6sLhzUzy4MD2ww5aJiZmVnW3jusclgzG8B/M5PxkEirnEOvmZm1QRdCWo/feJuVxxW2HDhwmJmZdVdXKmo9Dms2E/77mT4Htpw4tE3Gx82sGSQdKuk6ScsknTBg/V9JWpJ+XSXpYUnbSNpZ0n9JWirpaklvraL9ZkXoWlADv9k2G0TSjZKuTPvAS9Jl75e0MtM3Hp7Z/sS0P71O0gum2r+HRFplHNYe5Zkirc4kbQh8CngesAK4WNLZEXFNb5uI+AjwkXT7FwNvj4g7JG0CvDMiLpO0BXCppPOy9zVrmq6FtB6HNctLS2eMfHZErO5b9omI+Gh2gaR9gCOBfYGdgB9I2iMiHh62Ywe2HHmaf5sJhzarsYOAZRFxA4Ck04EjgGGh6yjgqwARcStwa/rzPZKWAvNG3NesNroazPo5qJnl6gjg9Ih4AFguaRlJP/vzYXfwkMicuWo0tdkr5eM0xL3zNuDeef63tNqZB9ycub0iXfYYkmYDhwLfGLBuV+ApwEX5N9FsZnrDG7Nf5rBmxWnZ31YA50q6VNKxmeVvlnSFpFMlbZ0uG7tP7XGFrQCutA3mkDa+Xmhzxc2mSw88xKzrV073bnN7Y+5TiyJiUXa3A+4z7EXuxcB/R8Qd67VLmkMS4t4WEXdPt4E2uU03euiR8LH0tu0rbk09OIyNp2VvqK2GyhwaOWH/CFP3kQDPiIhbJG0HnCfpWuAzwAdJ+ssPAh8DXs/0+lTAga0wDm2PclCbnIOblWR1RBw4Yv0KYOfM7fnALUO2PZJ0OGSPpI1IwtppEfHNmTTUZmacoNKmUOdgNjmHNStLA85nm6qPJCJuSb/fLuks4KCIuKC3XtJ/AN9Nb06nTwUc2ArVCypdDW4Oavnx+W1WsYuBhZIWACtJQtlf9G8k6XHAnwGvziwT8DlgaUR8vJzm2kxMFXLqGugczvLhoGZVaEBoG0rS5sAG6XnamwPPBz4gacf0PG6AlwFXpT+fDXxF0sdJJh1ZCPxi1GM4sJWga9U2B7ViuNpmVYmItZLeDHwf2BA4NSKulnRcuv6UdNOXAedGxB8yd38G8BrgSklL0mV/HRGLy2m95a0OVTqHs2I4rFmVGhzatgfOSj6fZBbwlYg4R9J/StqPZLjjjcBfAqT95xkkk2+tBd40aobI3k6tBF0IbQ5q5XBwsyqkAWtx37JT+m5/AfhC37KfMni8vrVYXlU6B7PyOKyZTSadQfnJA5a/ZsR9TgJOGvcxHNhK1NbQ5qBWDQc3M2sqB7H6cFCzOmlwla1QDmwlGxZumhrkHNaq5+BmZmaTcFizOnJoeywHtpqYTvCpQ7hzUKsfT0xiZmbjclizOnNoW58DWwONG5aKCHYOavXmapuZmU3FYc2a4Dlzr+XsqhtREw5sLZZnsHNQaxYHNzMz6+egZtZMG1TxoJL2k3ShpCWSLpF0UGbdiZKWSbpO0guqaF/XzF6pKb+sme6dt8Ej4c3MmsF9pBXBYc2suaqqsH0Y+PuI+J6kw9PbB0vah+SCrPuSXEjuB5L2mOraBGY2ms9vqw8HaBuD+0jLjYOaWfNVFdgC2DL9+XHALenPRwCnR8QDwHJJy4CDgJ+X30SzdvEwyXI5mNkMuI+0XDismbVDVYHtbcD3JX2UZFjm09Pl84ALM9utSJeZWU4c3PLncGY5exvuI21CDmlm7VNYYJP0A2CHAaveCxwCvD0iviHplcDngOcCg06WGjgjhqRjgWMBNtpy61zabNYl/SHDAW40hzLLU5l95Gbbz8mlzVZfDmlm7VZYYIuI5w5bJ+lLwFvTm18HPpv+vALYObPpfB4dCtK//0XAIoDNdti5+guTmTWcA1zCwczKUGYfudVe27mPbBkHNLNuqWpI5C3AnwE/Bp4D/DpdfjbwFUkfJzmheiHwiyoaaNZ1bQ9wDmZWY+4j7TEc0sy6q6rA9r+BkyXNAu4nHbYREVdLOgO4BlgLvMmzX5nVQ1UBZyZB0aHMGsp9pDmgmdkjKglsEfFT4IAh604CTiq3RWZWVw5d1jXuI7vJAc3MhqmqwmZmZmbWaQ5pZjYOBzYzMzOzEjigmdkkHNjMzMzMCuKQZmYz5cBmZmZmlhMHNDPLmwObmZmZ2YQc0MysaA5sZmZmLbblrPvXCxU/Wr1Xha2pHwcuM6s7BzYzM7MOqXNA+dHqvWrdPjOzKvgCR2ZmNiVJh0q6TtIySScM2eZgSUskXS3p/MzyrSSdKelaSUslPa28lluTOKyZWRNJulHSlWkfeEm67CNpv3eFpLMkbZUu31XSfem2SySdMtX+XWEzM7ORJG0IfAp4HrACuFjS2RFxTWabrYBPA4dGxE2Stsvs4mTgnIj4c0kbA7PLa72ZmVkpnh0RqzO3zwNOjIi1kv4ZOBF4T7ru+ojYb9wdu8JmZmZTOQhYFhE3RMSDwOnAEX3b/AXwzYi4CSAibgeQtCXwLOBz6fIHI+LOshpuZmZWhYg4NyLWpjcvBOZPui8HNjMzm8o84ObM7RXpsqw9gK0l/VjSpZJemy7fDfgt8HlJl0v6rKTNi2+ymZlZaQI4N+3/jh2w/vXA9zK3F6R94vmSnjnVzhUReTW0MpLuAa6ruh195gKrp9yqPG7PaG7P1OrWpja15wkRsW0ejZB0TtqW6dgUuD9ze1FELMrs8xXACyLijent1wAHRcTxmW0+CRwIHAJsBvwceCGwJckni8+IiIsknQzcHRF/O+0nZxOpYR/Zpv/dotStTW7PaHVrD9SvTZO2p+r+EaboI9N97xQRt6SnA5wHHB8RF6Tr3kvSP748IkLSJsCciPidpAOAbwH7RsTdwxrQlnPYrouIA6tuRJakS+rUJrdnNLdnanVrk9szWEQcWsBuVwA7Z27PB24ZsM3qiPgD8AdJFwBPBn4CrIiIi9LtzgQGTlpihalVH1mX/5WeurUH6tcmt2e0urUH6temOrSnoP6xt+9b0u+3SzqL5FSCCyQdDbwIOCTSKllEPAA8kP58qaTrSUapXDJs/x4SaWZmU7kYWChpQTppyJHA2X3bfBt4pqRZkmYDfwIsjYjbgJsl7ZludwhwDWZmZi0gaXNJW/R+Bp4PXCXpUJJJRl4SEWsy22+bTuaFpN2AhcANox6jLRU2MzMrSDrD1ZuB7wMbAqdGxNWSjkvXnxIRS9PhJlcA64DPRsRV6S6OB05Lw94NwDHlPwszM7NCbA+cJQmSbPWViDhH0jJgE+C8dN2FEXEcyURcH5C0FngYOC4i7hj1AG0JbIum3qR0dWuT2zOa2zO1urXJ7SlRRCwGFvctO6Xv9keAjwy47xKS8ftWjbr9bbo9U6tbm9ye0erWHqhfm+rWntxExA0kpwD0L/+jIdt/A/jGdB6jFZOOmJmZmZmZtZHPYTMzMzMzM6upRgU2SYdKuk7SMkmPmWVMiX9N118haf8C27KzpP+StFTS1ZLeOmCbgyXdJWlJ+vV3RbUn85g3SroyfbzHzDZT8jHaM/Pcl0i6W9Lb+rYp9BhJOlXS7ZKuyizbRtJ5kn6dft96yH1H/r3l2J6PSLo2/X2cJWmrIfcd+bvNuU3vl7Qy83s5fMh9yzpGX8u05UZJS4bcN/djNOx/vcq/I7N+deof08erXR/p/nFgO9xHTr897h8f3af7x7JERCO+SE50v57kIqwbA78E9unb5nCSi9IJeCpwUYHt2RHYP/15C+BXA9pzMPDdko/TjcDcEetLO0YDfn+3kVxPo7RjRHJi5/7AVZllHwZOSH8+AfjnSf7ecmzP84FZ6c//PKg94/xuc27T+4F3jfE7LeUY9a3/GPB3ZR2jYf/rVf4d+ctf2a+69Y/p49Wuj3T/OPCx3UdOvz3uHx/dp/vHkr6aVGE7CFgWETdExIPA6cARfdscAXwpEhcCW0nasYjGRMStEXFZ+vM9wFJgXhGPlbPSjlGfQ4DrI+I3JTzWIyK5aGH/zDtHAF9Mf/4i8NIBdx3n7y2X9kTEuRGxNr15Ick1rkoz5BiNo7Rj1CNJwCuBr870cabRnmH/65X9HZn1qVX/CI3tIzvVP4L7yEnaMyb3j+4fc9WkwDYPuDlzewWPffEfZ5vcSdoVeApw0YDVT5P0S0nfk7Rv0W0BAjhX0qWSjh2wvpJjRHLdpmEvImUfo+0j4lZIXmyA7QZsU9Vxej3JJ7yDTPW7zdub0yEopw4ZzlDFMXomsCoifj1kfaHHqO9/vc5/R9Ytte0foVZ9pPvH8dT5ta0ufaT7xz7uH4vVpMCmAcv6p7gcZ5tcSZpDMjXn2yLi7r7Vl5EMcXgy8G/At4psS+oZEbE/cBjwJknP6ltfxTHaGHgJ8PUBq6s4RuOo4ji9F1gLnDZkk6l+t3n6DLA7sB9wK8kwi36lHyPgKEZ/eljYMZrif33o3QYs89S8lrda9o9Quz7S/WN+utxHun/s4/6xeE0KbCuAnTO35wO3TLBNbiRtRPIHelpEfLN/fUTcHRH3pj8vBjaSNLeo9qSPc0v6/XbgLJKSc1apxyh1GHBZRKzqX1HFMQJW9Ya5pN9vH7BN2X9LRwMvAl4VEQNfsMb43eYmIlZFxMMRsQ74jyGPVfYxmgW8HPjasG2KOkZD/tdr93dknVW7/hHq10e6fxxb7V7b6tRHun98zGO7fyxBkwLbxcBCSQvST6SOBM7u2+Zs4LVKPBW4q1eSzVs6VvhzwNKI+PiQbXZIt0PSQSTH+3dFtCd9jM0lbdH7meRE3av6NivtGGUM/dSn7GOUOhs4Ov35aODbA7YZ5+8tF5IOBd4DvCQi1gzZZpzfbZ5typ638bIhj1XaMUo9F7g2IlYMWlnUMRrxv16rvyPrtFr1j1C/PtL947TU6rWtbn2k+8f19uv+sSxRg5lPxv0imcHpVySzyrw3XXYccFz6s4BPpeuvBA4ssC1/SlK6vQJYkn4d3teeNwNXk8x8cyHw9IKPz27pY/0yfdxKj1H6eLNJOpjHZZaVdoxIOsJbgYdIPs15A/B44IfAr9Pv26Tb7gQsHvX3VlB7lpGM4+79HZ3S355hv9sC2/Sf6d/HFSQvoDtWeYzS5V/o/d1kti38GI34X6/s78hf/ur/GvR3VvFrf636yGGvDxUfo0r7x/Qx3EdOvz3uHx99DPePJX0pPWBmZmZmZmZWM00aEmlmZmZmZtYpDmxmZmZmZmY15cBmZmZmZmZWUw5sZmZmZmZmNeXAZmZmZmZmVlMObFZ7ku4tc9+StpL0fzO3d5J0ZkGPv5GkS8dpl5mZWZb7R7NucGAze6ytgEc6pIi4JSL+vKDH+lPgZwXt28zMLE9b4f7RrHQObNZIkvaTdKGkKySdJWnrdPkfSfqBpF9KukzS7pLmSPphevtKSUdMsfsPAbtLWiLpI5J2lXRVuv/XSfqWpO9IWi7pzZLeIenytD3bpNvtLukcSZdK+omkvYY81qHA94Y8x5PS53GhpO0nO1JmZtYl7h/N2seBzZrqS8B7IuJJwJXA+9LlpwGfiognA08HbgXuB14WEfsDzwY+Jkkj9n0CcH1E7BcRfzVg/ROBvwAOAk4C1kTEU4CfA69Nt1kEHB8RBwDvAj495LGeDfx4wPLNgQvT53EB8L9HtNfMzKzH/aNZy8yqugFm0yXpccBWEXF+uuiLwNclbQHMi4izACLi/nT7jYB/lPQsYB0wD9geuG3CJvxXRNwD3CPpLuA76fIrgSdJmkPSGX490+9tMuB57ATcERFrBjzGg8B3058vBZ43YVvNzKwj3D+atZMDm7XJsE8FXwVsCxwQEQ9JuhHYdAaP80Dm53WZ2+tI/qc2AO6MiP2m2M9hwPeHrHsoIiL9+WH8v2pmZpNz/2jWYB4SaY0TEXcBv5f0zHTRa4DzI+JuYIWklwJI2kTSbOBxwO1pZ/Rs4AlTPMQ9wBYzaN/dwHJJr0jbIUlPHrDp0PH5ZmZm0+X+0aydHNisCWZLWpH5egdwNPARSVcA+wEfSLd9DfCWdPnPgB1Ixu0fKOkSkk8Trx31YBHxO+C/JV0l6SMTtvlVwBsk/RK4GljvRG5JGwILI2JkW8zMzEZw/2jWAXq0qmxmZZH0p8CrI+K4qttiZmZWF+4fzR7Lgc3MzMzMzKymPCTSzMzMzMysphzYzMzMzMzMasqBzczMzMzMrKYc2MzMzMzMzGrKgc3MzMzMzKymHNjMzMzMzMxqyoHNzMzMzMyspv4/DUbNZbA/AqwAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "f, ax = plt.subplots(1, 2, figsize=(15, 5))\n", "\n", "lt, la = np.meshgrid(loct, lati)\n", "\n", "c = ax[0].contourf(lt, la, dens)\n", "f.colorbar(c, ax=ax[0])\n", "ax[0].set_title(\"Density / g/cm3\")\n", "ax[0].set_ylabel(\"Latitude / deg\")\n", "ax[0].set_xlabel(\"Local time / h\")\n", "\n", "c = ax[1].contourf(lt, la, temp)\n", "f.colorbar(c, ax=ax[1])\n", "ax[1].set_title(\"Temperature / K\")\n", "ax[1].set_ylabel(\"Latitude / deg\")\n", "ax[1].set_xlabel(\"Local time / h\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Winds at 80 km" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "tags": [] }, "outputs": [], "source": [ "lati = np.arange(-90, 90, 10)\n", "loct = np.arange(0, 24, 3)\n", "\n", "xwind = np.zeros((len(lati), len(loct)))\n", "ywind = np.zeros((len(lati), len(loct)))\n", "\n", "for i, lat in enumerate(lati):\n", " for j, lt in enumerate(loct):\n", " out = mcm.run(\n", " altitude=80,\n", " latitude=lat,\n", " longitude=15,\n", " local_time=lt,\n", " day_of_year=53,\n", " f107=70,\n", " f107m=69,\n", " kp1=1,\n", " kp2=1,\n", " get_winds=True\n", " )\n", " xwind[i,j] = out.xwind\n", " ywind[i,j] = out.ywind" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 0, 'Local time / h')" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA24AAAFNCAYAAAB49jzWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAABTnUlEQVR4nO3dedwkZX3u/8/FMoMiq4ADDDCIIFGjgkiMKwoYYoxkUY8mRkz4yc9Eo0aNQDgnmpyQoCYuJzEnmagJRhRxIRI3BBMwJoICooJgFAUcGTZlc2HGge/5o6qhp6f37qq676rr/Xo9r+fp5em+n36661tXfe+qUkRgZmZmZmZm6dqq6QGYmZmZmZnZeA5uZmZmZmZmiXNwMzMzMzMzS5yDm5mZmZmZWeIc3MzMzMzMzBLn4GZmZmZmZpY4BzczMzMzM7PEObhZK0l6iaTPL/gYT5H0jQV+PyQ9bJExmJlZd02qQ5L+SdKfTXPfJY7pjZLet+Bj/Kakz8z5u2vK+rrNImMwy5GDm1WmXDD/cMhXSPrjpsc3SUT8R0Q8vKrHl7RC0q2SHrSkx/sNSe9fxmOZmdn0JF0raaOk3Qauv7yseWvmedxZ6lDVNWuZIuKMiHhmVY8vaS9J65b4eH8k6c+X9Xhm83Jws8qUC+YH9X8BrwZuAv6h2dEl4anA5RHxwyU93rOATy7psczMbDbfAV7YuyDpZ4EHzPtg7igt5FnAp5f8eK6v1jgHN6uNpEOAtwEviIj15XV7STpH0g8kfUvSS/vu/0ZJZ0l6r6S7JF0p6bC+20+SdE1529cl/eqU4zhd0mvLn/cut4b+Xnn5YeVYJOmI/i125RbV10n6qqQ7JH1Q0nZ9t/+hpPWSbpD0O1MMZWQhKJ/rD8vn+pGkd0t6iKRPlX/v+ZJ26bv/VsDRwKclbSfpfZK+L+l2SV+S9JBpXhszM5vbPwMv7rt8HPDe/jtIWinpLyVdL+kmSX8n6QHlbUdIWifpREk3Av84pA4dIumysg58EOivQYP3/RlJF5R14EpJz+m77Z8kvVPSJ8rHuljSAX23v0PSdyXdKelSSU+Z5gWQdKGkXy9/fnJZX59VXj5K0uXlz5vtzlDe72WSvinptnJsKm/bunzNbpX0beCXphjKuPoakn6vfK67JP1vSQdI+kL5954laUXf/XcBDgK+IGk3SR8vX9MfSPqPsv6a1cJvNquFpJ2BDwN/FhEX9N30AWAdsBfwXODPJR3Zd/tzgDOBnYFzgL/pu+0a4CnATsCfAO+TtOcUw7kQOKL8+WnAt8vvUHTB/iMiYsTvPh84BtgfeDTwkvLvOwZ4HUV4OhA4aopxPAv4xJjbf718vIOAXwY+BfwRsBvFZ/eVffc9HPh2RNxKsbKwE7AP8GDgZcBPphiPmZnN7yJgxzIwbQ38D2BwX7A3USzTHws8DNgb6N91YBWwK7AfcEL/L5Zh4l8oAuKuwIco6sQWJG0L/CvwGWAP4PeBMyT1T6V8IUXt3AX4FnBq321fKse4K/B+4EP9GyrH6K+vT2XL+nrhmN99NvB44DEUtfYXyutfWt52CHAYxbrCSOXf/lTgvDF3OwZ4HPAE4PXAWuA3Kermo+jrnJbj+GxE3AO8lmKdZXfgIRQ1edT6gtnSObhZ5cqtZqcDVwBv7rt+H+DJwIkRcXdEXA68C/itvl//fER8slxg/jPFAh2AiPhQRNwQEfdGxAeBb1IEmEkuBJ5SbiV7ajmmJ5W3PY3xheX/lM/5A4qi+Njy+ucD/xgRV0TEj4A3jhuApIcC20bEuB3J/zoiboqI7wH/AVwcEV+OiA3A2RRFrOeXuH/r4k8pAtvDIuKeiLg0Iu4cNx4zM1uKXtftaOBq4Hu9G8pa+FLgDyLiBxFxF/DnwAv6fv9e4A0RsSEiBje4PQHYFnh7RPw0Ij5MEbCGeQLwIOC0iNgYEf8GfJzNA8lHI+KLEbEJOIP76xkR8b6I+H5EbIqIvwJWAtPsP3chmwe1v+i7PKm+nhYRt0fE9cC/s3l9fXtEfLesvX8xYQxPBb5Svr6jvCki7oyIKynWTT4TEd+OiDsoNpKOq697AvuV/4NxG3rNls7BzepwIsUWrOMGFnB7AT8YWLheR7EFsufGvp9/DGynct6/pBer2PH7dkm3l8+x2Y7hw0TENcAPKYrCUyiK2Q3llshJhWVwPL0Di+wFfHfg7xinvxCMclPfzz8Zcrn/oCb900L+GTgXOLOctvnmcgukmZlV65+B36CYjfHegdt2Bx4IXNpXtz5dXt9zS0TcPeKx9wK+N1BHR9WavYDvRsS9A/cdV1/vqymSXivpKhW7BdxOMYtjYn0FvgAcVE7PfyzFa7CPioO2HA58bszvLqu+TrM/2lT1tX83hPK2t1B0Jz8j6duSTprwPGZL5eBmlZJ0BHAK8NyIuH3g5huAXSXt0HfdvvRtoRzzuPtRHODkFcCDI2Jniq1mmnJoF1JMt1hRdrQupNhKugtw+ZSP0W89xRSLnn0n3H/SNMmpSVpFsQXwMoByK+CfRMQjgCdSTDF58ZiHMDOzJYiI6ygOUvIs4KMDN99KEQoeGRE7l187lQfuuu8hxjz8emDv3r5fpVG15gaKwLTVwH2nqa9Podjg+nxgl7K+3sEU9TUifgxcCrwKuCIiNgL/BbwGuKaczj+rxuorxdTNayPiFoCIuCsiXhsRD6XYheE1A7t3mFXKwc0qU+5vdibw6oj48uDtEfFdigX6X6g4oMajgeMppmxMsj1FgbulfK7fpui4TetCitDX2/p3AcU+AJ8vp2XO6izgJZIeIemBwBtG3bHcEf3w8jmX4VnAp3tbYSU9XdLPlvtY3EkxtWOev8nMzGZ3PPCMctr8fcru1z8Ab5O0B9x3gKxfGPIYw3wB2AS8UtI2kn6N0bsHXAz8CHi9pG3Ljai/TFGTJ9mhfJ5bgG1UnL5nxynHCPfX197slQsGLs/qLIq/eXV5oJCRXS5J+wMrI+LqOZ9r0GazYyQ9W8VBzERRX+/B9dVq5OBmVXopxc6779CW53L7u/I+LwTWUGwdPJtibv+4HYoBiIivA39FUchuAn4W+M8ZxnYhRXHqBbfPU0xhGTeNY9x4PgW8Hfg3imkU/zbm7kcCXxgzHWZWg9NCVlEcCOZO4CqKv3Whk6Wamdl0IuKaiLhkxM0nUtSIiyTdCZzPdPuOUXavfo1iGuZtFAc/Gezq9d/3OcAvUnT6/hZ48ZSB5lyK/bz+m2Ja4t1sPlVxksH6Onh5Vv9QjukrFDNLhv7NpWl2Q5jFYH09kOJ/9kOK9Y+/HTjgmlml5H0qzeol6W8pppD87RIeaxuK/QIOKHeqNjMz6yRJnwT+JiIWDm/lfnqXA3v5ACSWikY7bpL+QMW5Ra6Q9IFyutyuks4rz69xnvrOVWXWEpdTdBeXYVfgfzm0WZUkPbw8EFDv605Jr/byulqukWYzu4DiiJTLsBPwGoc2m6RcNn9R0lfKZfaflNe/RdLVKs7Je3Z5aqzFnqup96OkvSmmpz0iIn4i6SyKdvQjKI40eFp5tJ5dIuLERgZpZmabKfed/B7wc8DL8fK6Eq6RZmZ5KPd53D4iflgexfvzFAfo2RH4t4jYJOlNAIsur5vex20b4AHldK8HUuzndCzFOb8ov/9KM0MzM7MhjqQ4Otx1eHldNddIM7PEReGH5cVty6+IiM+U50kEuAhYvehzNRbcykOw/yVwPcWhXu+IiM8AD4mI9eV91gN7NDVGMzPbwguAD5Q/e3ldEddIM7N8SNpa0uXAzcB5EXHxwF1+h+KgPwvZZtEHmFc5L/9YYH/gduBDkl40w++fAJwAoG1XPG7F7vXXrq031v6UrXTPiqZHML9Y4anvtriN137v1ojYffI9J3vKEdvFbT+4d/Id+1z5tZ9eSXHkuJ61EbF28H6SVlAcqe7khQZpEy2zRm61zYrHbbfzbDUy9eVylcveFSs2Tb5TJnbcZrqDF9+5abuKR5K+jRsbWyWemTZOe8ra/G24YV2j9RGmq5HlqaQeW+7HdrakR0XEFQCSTqE4xcY0p7saq8l36VHAd3onNZT0UYqTBd8kac+IWF+eB+zmYb9cvmBrAbbbe5/Y72WvqWnYs9vhum6u3N+1X/sWLBv27V5aX7P6ls0uX7tuKctP63PdS066blmPddsP7uUjn9htpt85eN/1d0fEYVPc9ReByyLipvLyVMtrm8vSauT2u+8TP3PsH0x8wpyW2VUviweXe7k6etV0pzM778aDKx5J2nKrayuvT3zLyhL99x+/ptH6CDPVSCLidkkXAMcAV0g6Dng2cOQyDnTTZHC7HnhCebLin1DsN3EJxQkjjwNOK79/rLERLkkdxbDucJhTgV9E24PaPCsni6zQ5FYcbQsv5P5pkgDn0LLldUIqr5FdWY7P49p1u7cmvFm7dCm05ULS7sBPy9D2AIoNb2+SdAzFuRufFhE/XsZzNRbcIuJiSR+mOJniJuDLFFsHHwScJel4isL1vImPtSLuW8Hu6hvaBXg52hTUUlzpcOjLVxkgjgb+/76rT2PG5bVNZ5k1sp9rhZnZ0u0JnF4edXkr4KyI+LikbwErgfOKA09yUUS8bJEnanRCb0S8AXjDwNUbKLYszmXWFe+uBj0r5BzUUgxmVZr373XgW45ya+GDB677Pgssr228ZdXIe1Y4sM3DXTdLjddZ0xQRXwUOGXL9w5b9XPnsiVmRWVbc/YHJXy5BzSsLyzPNa+lwZ2ZddvSqqzu/n1vqvA5q4OA2E4e8vKQc0hzM0jLt/8MBz6x5dS7b3XVrPy/XLScObhVZtLA4+M0m5ZAGDmpt4YBn1j0Ob9Y0rxNaj4NboqoMInUsAFIPUlVzke82BzwzMzNbtlYEtzadLLMOXQ9VVXBQs3n03jdLO0mNmVXCXTdrirtt1q8VwQ0WX3H2lm+bVluL97Qnap3EO7ibWRs5vFndHNpsUGuC26J8fikbpU2FelnhbFnP4ZBnZmZN8fqb5cbBbQnc7WuXnINaHcFsmRzyzCwn7rpZXdxts2Ec3BKwjCLg8De/3IpwbuFsWRzyzMz7aJtZlzm4tYS7ftNLPah1NZgt0zSvocOdmc2jrV03n4Q7He622SgObgbUE2aaCocpFliHs+Y53JnZvNoa3qx5Dm02TiuC247b3D3zirBXyOpX97mtUiyqDmx5Gff/8jLEzNrEyzSz9LUiuM1jkRVoL9yqtUjASzGs9Ti0tcuo/6eXD2bVWXn9imT2c+vVoJTrTm95NLi88nKqsGb1LZ3aVcTy19ngtgiHvjSkXCz7ObB1iwOdWbVSCm+QV4CzLTm8WU4c3Go2aSXeC9f2cGCzfg50ZsuTWniDPAKcDefwZrlwcEuM96nJXxsD23N3vAyAD995aMMjaZ+jV13NhU0PwixDKYY3cIAzs+o4uGXEW+zT1ubANuryIAc7s/TEimh6CJVJNbyBA1xu3HWzHDi4tYC7dM3qQmBb5u853JnVrz/ctO1w472/xwHOFpVCeNuw78bWfUZteRzcWs6hrjoObNU8h4OdWbXaGuJS7r6BA1wuUghvZqO0IrjtsvWPvaV/Dg5183Fgq5Y/y2b16QWdtgS41MMb+OTdZja/VgS3ac27ctrFlUTvT7clB7Z0uGtntlxt6sLlEt7A3bdUuetmqepUcJuXA9/9hoWXLoS5toW2XAPbtEb9fW38TJotWxtCXA7hDRzgUubwZilycKvQtCvHua9MtjnMObC1y7i/P/fPoVkVcg5xuYQ3cIBLlcObpabR4CZpZ+BdwKOAAH4H+AbwQWANcC3w/Ii4rZkR1mPYymTuK5G5hzkHtu5xl85Sk1qNzHF/uJzCGzjApcjhzVLSdMftHcCnI+K5klYADwT+CPhsRJwm6STgJODEJgfZhDauROYQ5hzYbFAbP4uWjSRrZG5duNzCGzjAmdlwjQU3STsCTwVeAhARG4GNko4FjijvdjpwAR0MbqO0bSUylTDnwGazattncZLUuj9tl0uNzCXE5RjewAEuFe662TiS9gHeC6wC7gXWRsQ7+m5/HfAWYPeIuHWR52qy4/ZQ4BbgHyU9BrgUeBXwkIhYDxAR6yXt0eAYs9Gm6ZZ1hjkHNlu2Fge6JLs/LZZdjUx9KmWu4Q0c4FLg8GZjbAJeGxGXSdoBuFTSeRHx9TLUHQ1cv4wn2moZDzKnbYBDgf8bEYcAP6Io+lORdIKkSyRdctsP7q1qjJaItgUs657n7nhZtsG6r/vzbii6PxFxO3AsRdeH8vuvNDG+llpajbznrh9VNcbspBoqp3Xtut0dHhpUV3DOdQNDV0XE+oi4rPz5LuAqYO/y5rcBr6eYqbKwJoPbOmBdRFxcXv4wRZG6SdKeAOX3m4f9ckSsjYjDIuKwXXZt8s+wOlTVcUttHztrv16AyyzE9Xd/vizpXZK2Z6D7AyTT/WmBpdXIrXfYvpYBr7x+RfbBKBdtD2/eWOvwlitJa4BDgIslPQf4XkR8ZVmP39hUyYi4UdJ3JT08Ir4BHAl8vfw6Djit/P6xpsZoaXC4mt6H7zw0t0DQab3/1f9c4mPeds8D55iW+YndJF3Sd8XaiFjbd7nX/fn9iLhY0juYoftjs8ulRuYY1HKeMtmv7dMnj151dZL131Mm8zVffYQpaiQAkh4EfAR4NcX0yVOAZ84z1lGaPqrk7wNnlPtLfBv4bYou4FmSjqeYD/q8BseXrRbsT1Ob8248uFVb9xzebA63RsRhY24f1v05ibL7U+5rNbL7Y3NLtkbmGNjaqm0Brr8edz28bdh3oz9raZhUI5G0LUVoOyMiPirpZ4H9ga9IAlgNXCbp8Ii4cd6BNBrcIuJyYNgLcWTNQ7FE1bXAdngzGy2X7k/bpFYjvQKZtmvX7d6a8Gb3c3hLn4pk9m7gqoh4K0BEfI2+3QckXQsctuhRJb1zWAu522bg94EtXa/781XgscCfUwS2oyV9k+KoWac1NzyrShv3XWvb39PTxoOXpLpR1SHZ+jwJ+C3gGZIuL7+eVcUTNT1V0mykuqdHtK3rBu682fKk1v2xarU12HRFrtMnc6vBnjJpABHxeUAT7rNmGc/ljpslqak57SnOpV+UO29mNq02dtdG6cLf2ZbuW8qBzqcIsDo5uLWMV9JtGL8vzGycLgW2rsll+uSkcJZyeDOri4ObJafprlfTz29mVodeWHNg64Ycwluu3HWzuji4tYi7KsvTxvDm94eZgbtr/br2OqQa3qbtprnr5vDWdT44iSUlpcDkg5XkY1wobePfazarrgUUGy3XA5f0dP3cbtZtDm5mHdOG8DZL97B339z/ZrN5aePYg51ZR/m8b/nyUSa7y8GtJdowDS7FLWht7LpBXuFtWe/t/sfJ5W83s+qtvH5FZ6efpRDe5qmx7ro5vHWVg5slIcUFcI/DW73q2AjhEGdmVkghvLWJp0xalRzcWqAN3TZrRtPhLYX3rkOcmXVdU/u9LbJRNNWuW53cdeseH1XSGtf1BW/T6gxPH77z0M2+UpPquMysWl75LeTWKUp1NkydAXjDvhs7O9W3i9xxy1TOK5e5BbVUC8MyLbvrlPP708zMbFF1T5l0960bHNwykuvKcG5BrV8XQtugWUJcru9JMzNrh5SnTDYR3sAd5DZzcEtcbivGqS48Z9XFwDZMbu8/MzPLwzLrrMPb5tx9ay8HtwTltLKc6oJyEQ5t3eYDlJiZ2TI1Fd7A3be2cXBLRC5hrY1BrZ9Dm5mZWX5S7rpBc6cJcIBrFwe3BuUQ1lJeCC6bQ5uZmVm+HN5G8/TJdnBwq0kOIa0n5YVeVRzazMzM6tHlmts7VYC7bzYPB7clyymg9etiWOvpcgExMzNrk9S7bj1Nd996HOLy4uA2p1wDWr8cFmxVc2izfj4wiZlZ/hzepjd48m4HubQ5uE3QhoDWL4cFWV0c2szMzKxJKYS3fg5yaXNwK7UtoPVzWNuSQ5uZmVl75dJ1g/TCWz8HubQ0HtwkbQ1cAnwvIp4taVfgg8Aa4Frg+RFx27Ker80BrV8uC6smOLSZWS7qrpFmVauzBju8LZ+DXLMaD27Aq4CrgB3LyycBn42I0ySdVF4+cdYH7UpA65fLwqlJDm1mlplKaqSZpSeX8NbPBzqpV6PBTdJq4JeAU4HXlFcfCxxR/nw6cAETitJt9zywk0ENHNam5cBmZrlZVo0067Kcum6QZ3jrcTeuek133N4OvB7Yoe+6h0TEeoCIWC9pjyYGlrKcFkApcGizafiIkpagt7OEGrn1xkn3MGs3h7dmOMgt31ZNPbGkZwM3R8Slc/7+CZIukXTJj25rf1U678aD7/uy2fg1s2l0tWs/LUnXSvqapMslXVJet6uk8yR9s/y+S9PjbItl1shNP/kRO1wX932ZjdI7OXQbHb3q6vu+crBm9S2t/n+0iaT3SLpZ0hUD1/++pG9IulLSm5fxXE123J4EPEfSs4DtgB0lvQ+4SdKe5ZbEPYGbh/1yRKwF1gLs/cidW1uJHDqWY9TrmMsC3OoxGN7chdvC0yPi1r7L3t+qOkurkdvvvs9mNbIX3u7aTxUO33LTpZDQX/tTX8/q/V/a0IFrsX8C/gZ4b+8KSU+nmNr+6IjYsKwZhI0Ft4g4GTgZQNIRwOsi4kWS3gIcB5xWfv9YU2NsUuoLkrZwoJus9xp18TXpD3IOcUN5f6uK1FEjHeCsp87Qdt6NBydVT3IJcbkHuDZPk4yIz0laM3D17wKnRcSG8j5DN7LNqul93IY5DThL0vHA9cDzGh5PrVJeaHRJFwLdrO+1/vu36XWYlrtxBPAZSQH8fdnR8T7J9Vt6jXSA23JfnC7pUqdtkhxCXO4BrkMOAp4i6VTgboqNb19a9EGTCG4RcQHFlloi4vvAkU2OpwmpLiBsczkEujreS13uwvWkGuTu3LTdHO+BT+zW22+ttLYMZv2eFBE3lOHsPEnd/efXrK4a6QDXPU2FttS6bsOkHuJyOoBJKt22+eojTFkjB20D7AI8AXg8xQa3h0bEQrt3JRHcuizFhYHNro5Al+J7petduH6pBrkp3RoRh427Q0TcUH6/WdLZwOFMub+V5cUBrhvcaZteqiEuh+5bKqFtQRNr5BDrgI+WQe2Lku4FdgMW+uA5uDUkpQ++VWfaQNeG94ND3OYyD3KbkbQ9sFVE3FX+/EzgT4Fz8D7JreUA114phLYcum7DpFi/cwhwHfQvwDOACyQdBKwAbh37G1NwcFtACh9Wy1Pb3zueSrmlzIPcQ4CzJUFRN94fEZ+W9CU6vE9yV7Q9wHVt/7YUQlub9OpcCnU9tQDXkm7bRJI+QHGgrt0krQPeALwHeE95ioCNwHGLTpOEjga3FD5cZl3gLtxomwe5TzQ2jmlExLeBxwy5vpP7JHdV2wNcF6QW2nLtug3jANddEfHCETe9aNnP1YrgNv/OhmZWF4c4s3ZwgMtTaqGtrVLaH67JANeVblvdtmp6AGbWPefdeHDjBc3MFrPDdXFfiLO0pRza2lwLjl51dRIbKtesviXp94BNrxUdNzPLk7twZvnLuQPX5v3bvKKejlSmUdZ1CgF326rj4GZmSXCIM8tbzgGuTXIMbG3a122cFAKc93/Lm4ObWcdNu/Cuc2XAR6U0y5cDXDNyDGxdlUqAqyK8udtWLQc3s4bktrWrf7x1rSC4C2e2uK3vDna6ZgN3HLCy1ufd4bpweKtYm8JaV7pu/VIIcCnzPrRbcnCzkXor6m0qDHXLLZxNq4n3hkOc2WJ2umYDQK0BLuXuW877t7kut0tTAW7ZXbdldtsc2oZzcLMtDH6IHeDGa2s4m0YTXTjwVEqzRTjAFVZevyLr8NZGXey69Tt61dXuvuHQNo6DW0csI1zM+xi5Br4uB7J5eCqlmU2S2vRJhzdLTd3hra4jTc7irv3k8DaCg1vGUvugjTLPOKte8c/ltctV0yEuP59oegDWQU3s9wYOb2Zm83JwS4wDRcGvQ3t4qq1ZuhzeCrmFt2vX7e5laovlOmVyw74bfVTJinUyuDkUmNWvqf3hzCxNDm9mo9UZ3jxdMh+tCG4bN26T3BvOzEZziDNLR1NdN0gzvEHeR5w0s/ZqRXAzS9U0Uwa6voLgEGdmqXH3rTldP7JkvxynTHq6ZLUc3MymUOVCaNhjd3WFwSHOrBnuum3J4c1SkGN4s+o4uFmn5LIVyGHOBzUx6xKHN7PR6ghvqe7nZptzcLPs5BK+lq2rYc5dOLN6NNl1g7TDG6S3vPWRJc26x8HNGtHV8LVsXQtzDnFm7ZZqeAN336xZOU2Z9H5u1XFws7n4A5muroQ5hziz5Wu66wYOb5aX5+54GQAfvvPQyp8rp/Bm1WgsuEnaB3gvsAq4F1gbEe+QtCvwQWANcC3w/Ii4ralxLpPDjjWl7WGu6hCX2rx/a78u1sh+Dm+Wul5g679cR3irUor7udnmmuy4bQJeGxGXSdoBuFTSecBLgM9GxGmSTgJOAk5scJwzcTizaaSwUtLWMOeiYy3RWI1MoesG3Hfy3aaXlcM4vHXbYGjrv77q8JZL183TJavRWHCLiPXA+vLnuyRdBewNHAscUd7tdOACEgpufhPmqbcCkJJhY2p6BaWtYc4sN7nWyCqksKFrGIe3bhoV2vpvd3izqiSxj5ukNcAhwMXAQ8qCRUSsl7RH3eNxOEtbiiFsWQb/thRWVgY/D15RMatXEzUyla5bj8PbcD6yZH0mBbbB++Ya3jxdMm2NBzdJDwI+Arw6Iu6UplswSzoBOAFg6wfvPNNzOpilrc3BbFYpBzkHOLPqLaNGrly580zPmVJg6+fwZpYXT5dcvkaDm6RtKQrSGRHx0fLqmyTtWW5J3BO4edjvRsRaYC3Ayv1Xb7Z26zdJHhzQZpfSikv/58wrLWbLt6waueMOq6de2KYa2npS3e+tyQ1a7rrV48N3Hjp1162Og5TkMlXS4W25tmrqiVVsNnw3cFVEvLXvpnOA48qfjwM+NvGxNoqV16+478vSssN1MfTL5pPi6+fPXjdI2lrSlyV9vLy8q6TzJH2z/L5L02Nsi2XWyGncccDK5ENbvxSXg9DchmNPbatH7keNbErbN+5K+gNJV0q6QtIHJG1X1XM1FtyAJwG/BTxD0uXl17OA04CjJX0TOLq8bBlwQKtXiq+tN6C03quAq/oun0RxhMMDgc+Wl205aqmRuQW2QSnWGYe3dpsU3txtG27DvhtbGeAk7Q28EjgsIh4FbA28oKrna/Kokp8HRs11OLLOsVSl7mJS19SR1Ipkl6U0dXKQ94VrF0mrgV8CTgVeU17duSMc1qWOGplzYBuU2hTKpvZ787TJeswybdI219Kpk9sAD5D0U+CBwA1VPpHNKbUAk9p4rB4phzfwvnAt8nbg9cAOfdc1fhRgm12bAtuglAJcUxuvHN7qMSy8tWUqZdVHlmxTeIuI70n6S+B64CfAZyLiM1U9n4PbBA5DloOUVlbGcReuehs3bjNPwd1N0iV9l9eWB7cAQNKzgZsj4lJJRyw+SmtCmwPboJSWiT7qZHv1h7e6QluO0ySHaSK8zVkfYUyNLPftPhbYH7gd+JCkF0XE+xYd7zATg5ukfwUG08sdwCXA30fE3VUMrE4OZ9YWqXffetyFS86tEXHYmNufBDyn3MdqO2BHSe9jyiMctlkuNbJLoa1fKgGu7vDmrlt92tJla0JGnbdxNfIo4DsRcQuApI8CTwQqCW7THJzk28APgX8ov+4EbgIOKi8nb9RBM1Lcqdlms9M1G4Z+dVlu72kfzCR9EXFyRKyOiDUUO13/W0S8iIqOcJiZ5GvkPdulvzGnainU+7qXdT5YieWgBRtvrweeIOmB5dGAj2Tzg3gt1TRTJQ+JiKf2Xf5XSZ+LiKdKurKqgc2q6QWyVWPeEDb4e13b2pxL562fu3BZOg04S9LxFMXreQ2PpwlZ1EgrpNCBq7P75s5bu7RlmuSgjDpvW4iIiyV9GLgM2AR8mfIcmlWYJrjtLmnfiLgeQNK+wG7lbUmsXW2dxCgmGxVCuhYq+tXVHRv2PG1/3XMMbz3eFy5dEXEBxdEjiYjv05KjAC8g+RppW2o6wDm8mW0u8/D2BuANdTzXNMHttcDnJV1DcWji/YHfk7Q9xeGfbcCsYWTe8JJD8Eh52mIXunJNr5wsyl04y4BrZMaaXEY6vNks6u62VX1kyWFyDm91mRjcIuKTkg4EDqYoSlf37Wz99grHlrQUAsk8Y1h2OEnhdViWNnflcu6+9bgLZylyjWyHpgJcncs1hzfLgcPbeNMcVfKBFCdb3S8iXirpQEkPj4iPVz+85rUpmED7/p6qtakr14bwBu7CWVq6XiPbpskA5/BmVuh9FhzgtjTNVMl/BC4Ffr68vA74ENCqouRAY9PIPci1Jbz1uAtnCehEjeyaJgJcneENcIDLSFsPSjKJa/uWpgluB0TE/5D0QoCI+El5uMvsOJzZsuU4vTL3/d6GWWYXzlv4bEatqZG2pbqXl97vzczGmSa4bZT0AMoTjEo6AEg6ATmgWZNy6cq1rfvW4+BlNcuuRtpkTZ8uALq939t5Nx7M0auubnoYSehqt82Gmya4vQH4NLCPpDOAJwEvqXJQZm2y0zUbkg5v/doY5Mwq5hqZuVSXe3VPnexJJchNCixtDXYpBDWfvD1d0xxV8jxJlwFPoDhi1qsi4tbKRzYnd9ssRSmHt379QS7VlRmzlORWI7ssx2VanVMne1INcoPGBZycQl0KQa3HgS19I4ObpEMHrlpfft+3PNnoZdUNaz4ObZayXMJbj0Oc2Wg51sguadMyq4nw1q9/ZT7VEDco5W5dSkENHNZyM67j9lfl9+2Aw4CvUGxNfDRwMfDkaodm1j69jQs5BTjwlEqzIbKpkfesGP2ZHfxs56gLy6Omw1tPLt24SeoMdqkFtR4HtjyNDG4R8XQASWcCJ0TE18rLjwJeV8/wpudum+Ukt+7bIAc567rcauQoi3526wx+XV/OpBLe+rUlyA1aJNg5qFmVpjk4ycG9ggQQEVdIemx1Q5rd1nfnv8XQuieXo09Ow9MqrcOSr5FVqir45bAcaeLIvCmGt345TqucR6rhrMchrb2mCW5XSXoX8D6Kwx2/CLiq0lGZdVCO54Qbxt046xjXyAV4+TC71MNbT1u7calxSOuWaYLbbwO/C7yqvPw54P9WNiIzu08bwpy7cdZyrpEd1FuuNXU+zFzCWz8HucU5pNk0pwO4G3hb+WVmDcs5zLkbZ23jGtk9qRzQpc4TdVfBQW48hzQbZpqOm5klLtcw5yBnZrlrquvWk2P3bZiu7B83yAHNZuHgZtZSOR78xNMqzZYvVsR9K/a9Lo3NJ5Vu26C2hLeeRcJMyqHPIc0WNe4E3CcDn46IL9c4nv7nPwZ4B7A18K6IOK2JcXTRtCv4PgVDXnLryrkbZylrskYuUh9Hrdw70E2WamjraVt4m9ei4WhZwc8hzaowruP2HeBVkh5DcWLRTwGfiYjbqh6UpK2BdwJHA+uAL0k6JyK+XvVzd9msK/GT7u9gl76cwpy7cZaYRmpkVfXRgW4xTU+X7HF4W5wDVzrWrL6F65oeRGLGnYD7TOBMAEmHAMcAHy2LxvkUWxq/WNG4Dge+FRHfLp//TOBYwMFtyapcSR/32A516cohzLkbZ01rsEbWWh/HhYAuhbrUu239HN4sZylPdU3BVPu4lVNBvgz8haQdKbb0/X9AVcFtb+C7fZfXAT9X0XN1UtMr4g51eUl9f7lZg1xOK2GWvpprZDL10V26zaXSdQOHN8uHg9psZj44SUTcCXyk/KrKsCXfZmtakk4ATgBYuXLnCofSHqmtbFu++oNciu8rBzNrSg01cmJ9hM1r5NYP3rmioQzXxi5dKoFsFg5vZu2T6lEl1wH79F1eDdzQf4eIWAusBdhxh9VeSxsjxRVrawe/t8xqN7E+wuY1cuX+6dTINoa6lDm8mbVLqsHtS8CBkvYHvge8APiNZoeUF69Qm5m10sz1ccWKTZtNR0r14AueemlmNt7E4CZJwG8CD42IP5W0L7CqwgOTEBGbJL0COJficMfviYgrq3q+NnFgMzOrT901chn1cdw+JSmGOge6xaTYdfP/rnmpvSdsOtN03P4WuBd4BvCnwF0Uc/cfX+G4iIhPAp+s8jnawmHNLF857DsjaTvgc8BKirrx4Yh4g6RdgQ8Ca4BrgefXccqYxNReI6usj5MOFJBSsHOgm9404c2vW7ekGOhzVee5p6cJbj8XEYdK+jJARNwmyZ/uBDiwWZP8/ltMDoGtzwbgGRHxQ0nbAp+X9Cng14DPRsRpkk4CTgJObHKgDehUjcwh2DnQDdf1v9+21HtPOMDNr+5zT08T3H5aDirKAe5OsXUxGfdspy1WItt6SPkurCzfccDK1v7/zDILbABERAA/LC9uW34FxfnDjiivPx24gO4Ft+RrZJ1SnoY5bOXUYcbM3bcF1XpuzWmC2/8Bzgb2kHQq8Fzgf1YxmGUaFXByDQRdCGxmbZZjYOtXhpNLgYcB74yIiyU9JCLWA0TEekl7NDrIZmRZI5uQYrfO3Tmzgrtvc6v13JoTg1tEnCHpUuBIivPH/EpEXFXVgKqWU6BzWLNU+b05vboDmzZqnpXO3SRd0nd5bXk4+ftExD3AYyXtDJwt6VGLjbQd2lYjmzQs2DXVpXOgs66qu/t27brdazsJ95z1EcbXyKnOrbksI4NbudN5z83AB/pvi4gfVDWoJoxbEa071Hml2Cx/mXXYbo2Iw6a5Y0TcLukC4BjgJkl7lt22PSlqRSfkVCN33OZujl519WbXnXfjwQ2NZjYphTlwoEvVDtclc6rCsXKpC+6+bWFcjZzq3JrLMq7jdilFYhSwL3Bb+fPOwPXA/lUNKjV1dOkc1szaIZfCPItyv62flqHtAcBRwJuAc4DjgNPK7x9rbpS1y7pGDga5YVINd4Nhrul958D7zy1LLgFsXvP8fU3WFAe4qdR67umRwS0i9geQ9HfAOeXhh5H0ixRFu/OWEegc2MzaoY2Brc+ewOnlfm5bAWdFxMclfQE4S9LxFGHleU0Osk5dqJG5hLvUunI97s5Np+1hbVEphD0fvGS0us89Pc3BSR4fES/rG+CnJP3vqgbUBpMCncOa5czv3821PLABEBFfBQ4Zcv33Kfbt6rJO18hUw12qYQ7cnXNQq96k13ieuuXu22h1nnt6muB2q6T/CbyPYlrIi4DvVzqqlvIKr1l7dCGw2VRcIydIJdylOMWyp+3dOYe1tPT+H/MGOIe35kwT3F4IvIHicMcAnyuvMzPrHAc2G+AauQSjwl2VgS7lrlxPrt25tgS1RY5lkMPG+h2uC3ffMjPN6QB+ALyqhrGY3ccn4U5TDoWoCg5rNoprZLXqPhpmyl25nhTDXGpBLYX1h0ljSKWeLtp9Awe4Ok0MbpL+nSHnI4iIZ1QyIjNLSirFpWoOZzYP18h6OcgN11txrjPAVR3WUghfVUot2M3bfQNPn6zTNFMlX9f383bArwObqhmOWaHtC+wctC2wOZhZRVwjG9Qf5OreTy7FENdEgLNqNHFAu6bCW50n4c7dNFMlLx246j8lXVjReMysQW0Iaw5oVifXyHQ4xN0v9wDnjbf3638t6qjRDm9pm2aq5K59F7cCHgesqmxEZlYrhzWz+blGpqmpKZVdCnB37afk9m1ru7pCnMNbuqaZKnkpxfx9UUz/+A5wfJWDMrPq5R7YHNYsEa6RGairG9fFALds7rZNp+oQ5/CWpmmC289ExN39V0jKe43PrKNyDmsOapYo18jM1BHiUp5Gae1T1f5wPmBJeqYJbv8FHDpw3ReGXGdmCXJYM6uUa2TG6phSmVIXbsO+G5PuuuXabVtx9To2Hry66WFU0oVbJLzNy1230UYGN0mrgL2BB0g6hGIaCMCOwANrGJuZLSDXwOawZjlwjWynKrtxqXThUg9vuVhx9bqRl9sW4hY5UbenTC7XuI7bLwAvAVYDb+27/i7gjyoc08zuWTH/yp53rLU2cVjLj6eSZCubGmnzqSPEpdCFS0UO3bbBsDbpfikEOFhOiGsqvNnmRga3iDgdOF3Sr0fER2ocU60WXWF08LOmOazlxUGtHbpSI61Q1ZTKpgKcu27TmzasTfrdlEJcTuHNNjduquSLIuJ9wBpJrxm8PSLeOuTXOsedPmtKjoGta2HNhaq9XCPn99wdL5vp/h++M73dBZfdjUtlGmVTUuu2LRLWpnnMpkPcIuFtXg5vyzFuquT25fcHDbnNqWNBw1ZgHeZsktzCmoOatZhr5ASzBrRpHiflEJdbF85dt81VEdamea6mQty84c1HmmzWuKmSf1/+eH5E/Gf/bZKetMiTSnoL8MvARuAa4Lcj4vbytpMpzoFzD/DKiDh3kefKyeAHwUHOwGEtVS4+3eYaWVhWOJv3+VIKcjl24RYNb8s4CXfT3bY6A9u4528iwDm85Wea0wH8NVse1njYdbM4Dzg5IjZJehNwMnCipEcALwAeCewFnC/poIi4Z4Hnypa7ct2WU2DrQlhzobERWl8j6w5ns0g1yOXaheuKpsPaME114Rze8jJuH7efB54I7D4wf39HYOtFnjQiPtN38SLgueXPxwJnRsQG4DuSvgUcTnFOHMNhru0c1tLggmKTtLFGphzQppVakMshwDU5ZbLObluKYW2UukOcw1s+xnXcVlDM3d8G2KHv+ju5v4gsw+8AHyx/3puiSPWsK6+zMTzFMm8Oa81z8bA5ZFsj2xDQppVKkMtxGmUb5BTWRqkrxDm85WHcPm4XAhdK+qeIuG7WB5Z0PrBqyE2nRMTHyvucAmwCzuj92rChjHj8E4ATALbZaZdZh9dqbejKDS48mp4Dv0w5BTVwWDMbJqcaudOeD5h1eK3UdPfNrK0WCW82m2n2cftxuaP0I4HteldGxDPG/VJEHDXudknHAc8GjoyIXuFZB+zTd7fVwA0jHn8tsBZgu733ySuVNCD3rlxKQS634DWPNi6AHdSsIsnXyL0fufNmC/zBANP2DlxKgW3ZJ/J2t228Xocq985b06cPqJK7brOZJridQTFN49nAy4DjgFvG/sYEko4BTgSeFhE/7rvpHOD9kt5KseP1gcAXJz1erIih/3Qf5na0rge5LoSvWTmsmc0l+Ro5SVuDXEqBzYa744CVtWyIzTXA1RnYFl0v8pTJekwT3B4cEe+W9Kq+qSEXLvi8fwOsBM6TBHBRRLwsIq6UdBbwdYrpIS9f5GhZs74Juhz0cp9e6SA2H4c1s4VlWyNHyT3IObDVaxmnBKhLLgGuzR22URzepjNNcPtp+X29pF+imJax0DsqIh425rZTgVMXefx5zfKG6ULIy70rZ6O1LbB5YW8Nan2N7A9CqYa4HMKap0kOV1fXrV+qAa6pwLasjd+L7uvm8DbZNMHtzyTtBLyW4tw0OwKvrnJQOehiN89BLm9tC2uQd2DrPyrczEe2sJR0qkam1o3LIbBZulIJcF3ssI3i8DbexOAWER8vf7wDeDqApFdXOKZWGvUmzDnQ9QcBh7g0Oayloz+o5UbSPsB7KY6CeC+wNiLeIWlXiv271gDXAs+PiNuaGmcTul4jmwpyDmzLk8J6SBNdt35NBbg2BrZlHGGyTeFN0mOBv6M4eNUm4PciYu59k6fpuA3zGuDt8z6p3W/wjZnCAnQe7salw2EtHTmHtQGbgNdGxGWSdgAulXQe8BLgsxFxmqSTgJMoDqrRdZ2tkVUHOQc2q1JdAS6lwJbqMQJaFN7eDPxJRHxK0rPKy0fM+2DzBrf2rRkmoi1Hx3SQq1cbwxrkF9haFNQ2ExHrgfXlz3dJuorixM/Hcn8BOh24AAc3cI28z7KCXO6Bzfu3TdZ0161fVQEupcBWpWWd1y3H9d8hgmIKPcBOjDiFy7TmDW5eC69RG7pyDnLL57CWhraGtVEkrQEOAS4GHlKGOiJivaQ9mhxbQryAG2GWA53kHta6JKcjS85iWQGuK4HNhno1cK6kvwS2Ap64yIONDG6S7mJ48RHwgEWe1BbjINddbQ1rkFdgSzmsbb1xrs/TbpIu6bu8tjyB82YkPQj4CPDqiLizPFR9J7lGLm5UN86BrbtS6rr1mzfApR7Yqpwmuayu2zLNWR9hQo2UdD7FPuCDTgGOBP4gIj4i6fnAu4Gj5hkEjAluEbHDvA9atxUrNg1dkWrjVIJh2hbkHOI2l9qCb5lyCWspB7UluTUiDht3B0nbUoS2MyLio+XVN0nas+y27QncXPVAU5FTjcxFWwPbsqdJWnOmDXCpBzab2dgaGREjg5ik9wKvKi9+CHjXIgOZd6pkFiatbLU12OUe5NyNc1hLQQfC2tRUtNbeDVwVEW/tu+kc4DjgtPL7xxoYnlmntHXdpSfVrlu/UQEup8CW6kFJWugG4GkU+4A/A/jmIg/W6uA2ybgVszYtGHM/4ElXglybwxqkH9gc1MZ6EvBbwNckXV5e90cUge0sSccD1wPPa2Z4Ns6dm7bbrOtz9KqrGxyNWXvkFNSakOJ0yQa8FHiHpG2Au4ETFnmwTge3cdoe6nLuyrUpyLV9geaw1g4R8XlGHynxyDrHYotziLOeVGt/Dl03s2mU9fNxy3o8B7c5tDHUtSnIWbMc1szy4RBXrS6dBqCtR5ZsozqnSbrrtlwObks2uFKY8kJ2nJyDnNXPYc0sfw5xlpLBcOEO3GKa2KfNgW35HNwq1tYgBw5zXZZ6UAOHNbNF9EKcA9z8Uu+2LbuGV91tGxY8HOaGS+HAIw5t1XBwq1lbghw4zHWJg5pZNznAzadroa0pDnNphLR+DmzVcnBrWJuCHHiKZRvkENLAQc2sTg5w0+tiaEtp37Y2T7FMLaT1c2Crh4NbYtoe5MBhLhW5BLQeBzWz5p1348EOb2OkHNq6Wntz7cqlHNL6ObDVy8EtcW0LcuAwV7fcAlqPg5pZmtx9G66roS2lbtu0UgtzuYS0fg5szXBwy0wbgxw4zC1DrgGtx0HNLC8OcPframhrk7rCXI4hbZBDW3Mc3DLX1iAH3l9ulNwDWo+Dmlk7dD3AdTm05dhtm8WiYa4NIa2fA1vzHNxapktBDtof5toS0sBBzaztuhbglh3YIK/Q1lWjwlzbQlo/B7Z0OLi1XJuDHMwfbFIraG0KaD0Oambd1MYAV0VIG+TQlq+2hjYHtvQ4uHVM24PctNoYlJrmoGZm/XI9AmUdIW1QjqGt7dMku86hLU2tCG47bnP3ZsWhiYVurvpXtrsa4mw+DmpmNknK3bdU1hVyDG3WXg5saWs0uEl6HfAWYPeIuLW87mTgeOAe4JURce6sjzuqQKSykE6Vu3E2icOaWX2qqpFNaDrApVr/cw1t7ra1jwNbHhoLbpL2AY4Gru+77hHAC4BHAnsB50s6KCLuWcZzDisYqS7MUzBqJd2Brv0c0Mya1USNrEMdAS6Xup5raLN2cWDLS5Mdt7cBrwc+1nfdscCZEbEB+I6kbwGHA1+oahAOc7MbtlLvMJcvhzSzJCVRI6uyrACXa71eVs1sIrC529YeDm35aSS4SXoO8L2I+Iq02Ztmb+CivsvryutqNVhIci0MdXKYy4NDmln6Uq+RyzRLgGtLLc45tFk7OLDlq7LgJul8YNWQm04B/gh45rBfG3Ld0E07kk4ATgDYac8HzDlKq9Ka1bc4vCXCgc0sLXXWyJV77DDnKOszeATKtoS0QbmHNnfb8uFw1k6VBbeIOGrY9ZJ+Ftgf6G1JXA1cJulwiq2H+/TdfTVww4jHXwusBdj7kTtXuiQ5etXVrS0i1n4ObWbpqbNG7vDwVVmsbbe9zua8IdOBLT0OZt1U+1TJiPgasEfvsqRrgcMi4lZJ5wDvl/RWih2vDwS+WPcYzdrAgc0sP66RNonP09YNDmY2TFLncYuIKyWdBXwd2AS8PKejZZmlwqHNrH1cI/OVU7fNga0+Dmc2q8aDW0SsGbh8KnBqM6Mxy5sDm1m7uEZavyq7bQ5sy+dgZsvWeHAzs8U5sJmZpSnlA5I4rC3O4czq5OBmlfKRJavn0GZm42zc6FJvm3NgW4zDmjXFS3OzTDmwmdm0rl23u5cZDUit2+bAthgHNmvaVk0PwMxm5xUwq5Ok90i6WdIVfdftKuk8Sd8sv+/S5BhtMs9+6DaHtvnctZ/u+zKblaTnSbpS0r2SDhu47dGSvlDe/jVJ2016PAc3s4ysWX2LQ5s14Z+AYwauOwn4bEQcCHy2vGyJc3irT0rdNoe22Tis2RJdAfwa8Ln+KyVtA7wPeFlEPBI4AvjppAfzVEmzDDisWZMi4nOS1gxcfSxFoQE4HbgAOLG+UZm1X13nbDNPg7RqRMRVANIW769nAl+NiK+U9/v+NI/njptZ4hzaLFEPiYj1AOX3PSbc3xLhrlu3uNs2mjtr1qCDgJB0rqTLJL1+ml9yx80q5yNLzseBzeax9d3BTtdsmPXXdpN0Sd/ltRGxdonDssT4YCXVWkbN8xTJ6jiopW/DvhuX/phz1keYUCMlnQ+sGvJ7p0TEx0Y85jbAk4HHAz8GPivp0oj47LiBOLiZJcgrVFazWyPisMl328xNkvaMiPWS9gRurmJgVh2HN+sSh7X0VRHWlmRsjYyIo+Z4zHXAhRFxK4CkTwKHUuwzPpKnSpolxAcfsYycAxxX/nwcMGqroiXMsyGWz922dHgqZPo27Lvxvq+OORd4tKQHlgcqeRrw9Um/5I6bWSIc2CxVkj5AcSCS3SStA94AnAacJel44Hrgec2N0BbhzltaHNoW45CWvi6FNEm/Cvw1sDvwCUmXR8QvRMRtkt4KfAkI4JMR8YlJj+fgNsF5Nx7c9BCy5K24s/PKk6UqIl444qYjax2ImdkEO1wXDm8J6lJY6xcRZwNnj7jtfRSnBJhaZ4KbA9jsHL7qN+o1d6AzMzOznHQ1rFWpFcHtzk3bOZgtwAEtfeP+Rw51ZrYod/zbo8vTJK15DmvVakVws+k5pLWPu3RmZmb383TJejms1cfBrcUc0rrNgc7MzMyq4LDWDAe3FnBAs1l42qWZWXXWrL7FddlayWGteQ5umXExsCpN06XL+T3oYGpm1g2eLrk8DmzpcHBLVM4rx9Y+bXk/evqombWZD0xiy+KwliYHtwS0ZaXYLFeePmrmI0ta+7jrNhuHtfQ5uC2grYFr5fUran0+LygsZe7SmZlZW3kdLC+tCG4bN27T2hBVh7qDWtXP74WQ1cFdOjMbxQcosZR5PSlfrQhuNr2mQ1odRv2NXlBZXdylMzNLg6dLFrwO1A6NBTdJvw+8AtgEfCIiXl9efzJwPHAP8MqIOLepMbZBF4LatIa9Fl6QWZ28Bd6m5RppZsvg9Zx2aSS4SXo6cCzw6IjYIGmP8vpHAC8AHgnsBZwv6aCIuKeJcebIQW02DnNmlhrXSJuHjyhpPV6Paa+mOm6/C5wWERsAIuLm8vpjgTPL678j6VvA4cAXmhlm+hzUls9TLc2sYa6RZkvUlemSXk9pv60aet6DgKdIuljShZIeX16/N/DdvvutK68zikAx+GX18etvZjVprEZ6Ou/yeJ9Wq8OGfTfe92XtV1nHTdL5wKohN51SPu8uwBOAxwNnSXooMGxzyNDev6QTgBMAtn7wzksYcXqaDgbDpl10YYvVLNydq98snwv/HyxVKdZIB400bNh341z1/6795OmSY3Sl62btVllwi4ijRt0m6XeBj0ZEAF+UdC+wG8XWw3367roauGHE468F1gKs3H91a5ZUKYa1WW4fposLSge6zTX1vvb/wVKVUo10YLOu6K3DdHG9xNqhqX3c/gV4BnCBpIOAFcCtwDnA+yW9lWLH6wOBLzY0xlqkHtSaeI42L1Cb/n9bwYHOEvcv1FQjHdqqNe/53Nx1q5YDnOWqqeD2HuA9kq4ANgLHlVsWr5R0FvB1ikMgv3yao2Vpo7xCPKUcFugOetYUBzpLxFJr5DAObGYOcJafRoJbRGwEXjTitlOBU+sdUbvlENYWMerv84LYlsWBzupUZY10YMuHu271cYCzXDR2Am6rjhfYBR9cxarmQGc5cWjrDoe3+eR4AJOV169wzekQB7eW8AJ6Og5zVodxW8ldYK1uK1Zscmhr0Lz7ucH8XTebn7tvljIHt0w5qC2Pw5zVySthZlYHd90W4wBnKXJwy4QXvvVymDMzsxQs0nVzeFucA5ylxMEtYV7YpmXw/+GFuJmZTWOR6ZKWBgc4S4GDW0JSCWo7XbOh6SFs4Y4DVjY9hC24K2ddIukY4B3A1sC7IuK0hodkZlNw1225HOBsFpLeAvwyxaldrgF+OyJu77t9X4pTvLwxIv5y0uM5uDUslYVpimGt3zzjayLsOcxZG0naGngncDSwDviSpHMi4uvNjsysG3yQkvQ4wNmUzgNOjohNkt4EnAyc2Hf724BPTftgDm41SSWg9aQe1JZh8G9sqmvnMFe9SZ8vv94LOxz4VkR8G0DSmcCxFFsJzWwKTU6XdNetOjmeQsDqExGf6bt4EfDc3gVJvwJ8G/jRtI/n4LagXBaEXQhqk6QS5GD6942LwXI+Y369F7Y38N2+y+uAn2toLGadtGjXLZfwllKtnpa7bzal3wE+CCBpe4rO29HA66Z9AAe3EXJYuE3isDZeDsVhlvdhbgUjxc9YGwKe7t7IiqvXzfpru0m6pO/y2ohY2/+wQ34nvX+gmY2Vengbtt7Suy7FGj3I3be0zVkfYUKNlHQ+sGrI750SER8r73MKsAk4o7ztT4C3RcQPpenfM50MbikvtBbhoLaYYa9fDoWiJ6WQ19bPWE8bAt6AWyPisDG3rwP26bu8Grih2iGZ2aBl7OuWYnibZv0llwDn8NZKY2tkRBw17pclHQc8GzgyInofvp8DnivpzcDOwL2S7o6Ivxn3WK0LbqktjKrmsFat/tc39WIxi2k+J+MKT9c+Z/Na9HVOyJeAAyXtD3wPeAHwG80OyayblhXeIJ1l+R0HrJxqfSaXOuzwZj3lEZlPBJ4WET/uXR8RT+m7zxuBH04KbdCS4Lb1xnQWPnVxYKtfW0PcKP2fqRS30Fp9yqNhvQI4l+J0AO+JiCsbHpZZVpZ5YJJlHWUypQA3LrzlWHPrDG8rr1/Bhn031vJcNrO/AVYC55VTIi+KiJfN+2CtCG5d4sCWhi6HOFuuXF7biPgk8Mmmx2FmhWWeIiCVjXO9etqWGuvwZhHxsCnu88ZpH8/BLSMObWnKZd69mZk1o6rTACw7vEEaG5PaVE89bdKWaaumB2CT7XTNBoe2DPT+T/5fmZlZT9Xnblt2l+Wu/eSgsWR1hWGfpL39HNwS5hCQL//vzMysrhNuVzFFzuFtuRzebBkc3BLllf52cBfOzKyb6gptPVWFNwe4/Di8tZf3cUuMV/CnM+kEihsPXl3TSKbXlp2tzcwsTcvc561fSvu/5cz7u9miHNwS4cA22jxnue/9TooBDnxAkyr1F0WvZJhZE+rutvWrKrxBOkefzFld4c1HmWwnB7cEOLTdb56QNu3jpRji3IVbjlFFcNj1XukwS1cv8KxZfUvDI5lfk6Gtp7fC7u5bmhzebF4Obg3qcmBbdkCb5TlTDHDgLtys5i16DnNm6RkMO+PCT8qhLoXQ1q/q7ht4+Tkvhzebh4NbQ7oU2poIaeO4C5e3Kgqdw5xZM+YJOqmGutRCW0+V4Q0c4Mzq1Ehwk/RY4O+A7YBNwO9FxBfL204GjgfuAV4ZEec2McZl6kpISy2gTSOXLhx0O8Q1sTO3w5w1pe01ssqAM+mxqwp2qYa2nqrDG3j/t3m462azaqrj9mbgTyLiU5KeVV4+QtIjgBcAjwT2As6XdFBE3NPQOGfmkJan1AMcdC/EpXjkrVFj8sqKLVkra2QK4SbVbl0d6gpv4GXiLBzeRluz+haua3oQiWkquAWwY/nzTsAN5c/HAmdGxAbgO5K+BRwOfKH+IU7mkNY+qU+j7Gnr/nAphrVpuDtnS9aKGtmTQmCbxryhLpe/D+oJb+AANyuHty21fUPKvJoKbq8GzpX0lxQnAX9ief3ewEV991tXXte4roS0QV0KbYPchatHrmFtksG/yyswNoNXk1mNHCanQDNJm/6WKo84OcjTJ6fnc7zdz6FttMqCm6TzgVVDbjoFOBL4g4j4iKTnA+8GjgKGvWOHfuIlnQCcALBi+12WMmbobkAb1OXANii3LhzkFeJ6Rb2tBcsrLTZMnTVy5R47LGXMZpa/XmDPpfNmm6ssuEXEUaNuk/Re4FXlxQ8B7yp/Xgfs03fX1dw/RWTw8dcCawG2332fmdeMHNA256A2nRy6cJDnVMpRASfFQOcwZouqs0bu8PBVjbxhe1vN29StaoM6Om02u7prXcrTJq9dt7u7biM0NVXyBuBpwAXAM4BvltefA7xf0lspdrw+EPjiIk/kgLYlh7TF5daFyynADaor0DmMWUJqq5F1WLP6Foe3RDi0pSnFDZRNc3gbrqng9lLgHZK2Ae6mnM4REVdKOgv4OsUhkF8+7dGyHNBGc1CrVg5duFynUY4zTaBzGLNMLb1GNs3dt+Y5tKWpydCWctcNvLwYppHgFhGfBx434rZTgVNnebyt7/bKWT8HtWa4C5cGhzXL3bJrZEoc4OrnwGbjpB7ebHNNddxsiRzU0pNDiGtjF87M8uDpk9VzYEtfKlMkHd7y4eCWIQe1vOQ0ldIBzszq4u5bNVIKbJ4BMVoqoa3H4S0PDm6Jc0hrD3fhzMy25AC3uJTCmk2WWmjrcXhLn4NbYhzUusFdODOzzXn65Owc2GzZHN7S5uDWMAe1bsspwPVzmDOzKrj7Nh0Htnyl2m3r5/CWLge3mjmo2TA5TKPsN+70Gw51ZrYod9+Gc2DLWw6hrcfhLU0ObhVySLN55NCFG2dUqHOgM7NZOLzdz4HNmuDwlh4HtyVyULNlyj3ADXKXrp0kPQ94I/AzwOERcUnfbScDxwP3AK+MiHMbGaRlq+vhzYGtPXLqtvVzeFuMpP8NHAvcC9wMvCQibpB0NHAasALYCPxhRPzbpMdzcJuTQ5rVpW0BbhiHuqxdAfwa8Pf9V0p6BPAC4JHAXsD5kg6KiHvqH6LlrIvhzYGtXXINbT0Obwt5S0T8LwBJrwT+GHgZcCvwy2WIexRwLrD3pAdzcJuSg5o1rQsBbpi2Tb0cF1JzFBFXAUhbrJgcC5wZERuA70j6FnA48IV6R2ht0JXw5sDWPrmHth6Ht/lExJ19F7cHorz+y33XXwlsJ2llWTNHcnAbwiHNUpbbgUyq0rYA1EJ7Axf1XV7HFFsTzUZpc3hzYLMcOLzNR9KpwIuBO4CnD7nLrwNfnhTaABSR/1ntJd0FfKPpcQzYjaINmgqPZzyPZ7LUxtSm8ewXEUtZI5X06XIss9gOuLvv8tqIWNv3mOcDq4b83ikR8bHyPhcAr+vt4ybpncAXIuJ95eV3A5+MiI/MODZbUII1sk2f3aqkNiaPZ7zUxgPpjWne8TRdH2EJNbK838nAdhHxhr7rHgmcAzwzIq6ZNJC2dNy+ERGHNT2IfpIuSWlMHs94Hs9kqY3J4xkuIo6p4DGPmuPX1gH79F1eDdywnBHZjJKqkal8VnpSGw+kNyaPZ7zUxgPpjSmF8VRRH8vHnbZGvh/4BPAGAEmrgbOBF08T2gC2mmuEZmZm450DvEDSSkn7AwcCX2x4TGZmZrWRdGDfxecAV5fX70wR4k6OiP+c9vEc3MzMbG6SflXSOuDngU9IOhcgIq4EzgK+DnwaeLmPKGlmZh1zmqQrJH0VeCbwqvL6VwAPA/6XpMvLrz0mPVhbpkqunXyX2qU2Jo9nPI9nstTG5PEkICLOppjqMey2U4FT6x2RDZHae9PjmSy1MXk846U2HkhvTKmNpxYR8esjrv8z4M9mfbxWHJzEzMzMzMyszTxV0szMzMzMLHFZBTdJx0j6hqRvSTppyO2S9H/K278q6dAKx7KPpH+XdJWkKyW9ash9jpB0R9/c1T+uajx9z3mtpK+Vz3fJkNvrfI0e3ve3Xy7pTkmvHrhPpa+RpPdIulnSFX3X7SrpPEnfLL/vMuJ3x77fljiet0i6uvx/nF3usDrsd8f+b5c8pjdK+l7f/+VZI363rtfog31juVbS5SN+d+mv0ajPepPvI7NBKdXH8vmSq5Guj0PH4Ro5+3hcH+9/TNfHukVEFl/A1sA1wEOBFcBXgEcM3OdZwKcAAU8ALq5wPHsCh5Y/7wD895DxHAF8vObX6VpgtzG31/YaDfn/3UhxPo7aXiPgqcChwBV9170ZOKn8+STgTfO835Y4nmcC25Q/v2nYeKb53y55TG+kOCfXpP9pLa/RwO1/BfxxXa/RqM96k+8jf/mr/yu1+lg+X3I10vVx6HO7Rs4+HtfH+x/T9bHmr5w6bocD34qIb0fERuBM4NiB+xwLvDcKFwE7S9qzisFExPqIuKz8+S7gKmDvKp5ryWp7jQYcCVwTEdfV8Fz3iYjPAT8YuPpY4PTy59OBXxnyq9O835Yynoj4TERsKi9eRHG+q9qMeI2mUdtr1CNJwPOBDyz6PDOMZ9RnvbH3kdmApOojZFsjO1UfwTVynvFMyfXR9bESOQW3vYHv9l1ex5ZFYJr7LJ2kNcAhwMVDbv55SV+R9CkVZ0evWgCfkXSppBOG3N7IawS8gNELk7pfo4dExHooFjrAsMOvNvU6/Q7FFt9hJv1vl+0V5dSU94yY5tDEa/QU4KaI+OaI2yt9jQY+6ym/j6xbkq2PkFSNdH2cTsrLtlRqpOvjANfHeuQU3DTkusFDYk5zn6WS9CDgI8CrI+LOgZsvo5j68Bjgr4F/qXIspSdFxKHALwIvl/TUgdubeI1WUJx08ENDbm7iNZpGE6/TKcAm4IwRd5n0v12m/wscADwWWE8x/WJQ7a8R8ELGb02s7DWa8Fkf+WtDrvOhfG3ZkqyPkFyNdH1cni7XSNfHAa6P9ckpuK0D9um7vBq4YY77LI2kbSneqGdExEcHb4+IOyPih+XPnwS2lbRbVeMpn+eG8vvNFOdWOnzgLrW+RqVfBC6LiJsGb2jiNQJu6k1/Kb/fPOQ+db+XjgOeDfxmRAxdcE3xv12aiLgpIu6JiHuBfxjxXHW/RtsAvwZ8cNR9qnqNRnzWk3sfWWclVx8hvRrp+ji15JZtKdVI18ctntv1sUY5BbcvAQdK2r/cQvUC4JyB+5wDvFiFJwB39Fq1y1bOJX43cFVEvHXEfVaV90PS4RSv9/erGE/5HNtL2qH3M8UOvVcM3K2216jPyK1Adb9GpXOA48qfjwM+NuQ+07zflkLSMcCJwHMi4scj7jPN/3aZY+rfr+NXRzxXba9R6Sjg6ohYN+zGql6jMZ/1pN5H1mlJ1UdIr0a6Ps4kqWVbajXS9XGzx3V9rFskcISUab8ojvj03xRHoTmlvO5lwMvKnwW8s7z9a8BhFY7lyRQt3a8Cl5dfzxoYzyuAKymOlHMR8MSKX5+Hls/1lfJ5G32Nyud7IEWh2anvutpeI4qCuB74KcXWneOBBwOfBb5Zft+1vO9ewCfHvd8qGs+3KOZ5995Hfzc4nlH/2wrH9M/l++OrFAvSPZt8jcrr/6n3vum7b+Wv0ZjPemPvI3/5a/Br2Pus4WV/UjVy1PKh4deo0fpYPodr5OzjcX28/zlcH2v+UvnCmZmZmZmZWaJymippZmZmZmbWSQ5uZmZmZmZmiXNwMzMzMzMzS5yDm5mZmZmZWeIc3MzMzMzMzBLn4GbJk/TDOh9b0s6Sfq/v8l6SPlzR828r6dJpxmVmZtbP9dGsWxzczLa0M3BfYYqIGyLiuRU915OB/6rosc3MzJZpZ1wfzRrj4GZZkvRYSRdJ+qqksyXtUl7/MEnnS/qKpMskHSDpQZI+W17+mqRjJzz8acABki6X9BZJayRdUT7+SyT9i6R/lfQdSa+Q9BpJXy7Hs2t5vwMkfVrSpZL+Q9LBI57rGOBTI/7GU8u/4yJJD5nvlTIzsy5xfTRrLwc3y9V7gRMj4tHA14A3lNefAbwzIh4DPBFYD9wN/GpEHAo8HfgrSRrz2CcB10TEYyPiD4fc/ijgN4DDgVOBH0fEIcAXgBeX91kL/H5EPA54HfC3I57r6cAFQ67fHrio/Ds+B7x0zHjNzMx6XB/NWmqbpgdgNitJOwE7R8SF5VWnAx+StAOwd0ScDRARd5f33xb4c0lPBe4F9gYeAtw45xD+PSLuAu6SdAfwr+X1XwMeLelBFEXxQ331b+WQv2Mv4AcR8eMhz7ER+Hj586XA0XOO1czMOsL10azdHNysTUZtJfxNYHfgcRHxU0nXAtst8Dwb+n6+t+/yvRSfqa2A2yPisRMe5xeBc0fc9tOIiPLne/Bn1czM5uf6aNYCnipp2YmIO4DbJD2lvOq3gAsj4k5gnaRfAZC0UtIDgZ2Am8ui9HRgvwlPcRewwwLjuxP4jqTnleOQpMcMuevI+ftmZmazcn00azcHN8vBAyWt6/t6DXAc8BZJXwUeC/xped/fAl5ZXv9fwCqKef2HSbqEYuvi1eOeLCK+D/ynpCskvWXOMf8mcLykrwBXApvt8C1pa+DAiBg7FjMzszFcH806RPd3m82sLpKeDLwoIl7W9FjMzMxS4fpoNpqDm5mZmZmZWeI8VdLMzMzMzCxxDm5mZmZmZmaJc3AzMzMzMzNLnIObmZmZmZlZ4hzczMzMzMzMEufgZmZmZmZmljgHNzMzMzMzs8T9P3ePzOhtNRIEAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "f, ax = plt.subplots(1, 2, figsize=(15, 5))\n", "\n", "lt, la = np.meshgrid(loct, lati)\n", "\n", "c = ax[0].contourf(lt, la, xwind)\n", "f.colorbar(c, ax=ax[0])\n", "ax[0].set_title(\"Zonal wind / m/s\")\n", "ax[0].set_ylabel(\"Latitude / deg\")\n", "ax[0].set_xlabel(\"Local time / h\")\n", "\n", "c = ax[1].contourf(lt, la, ywind)\n", "f.colorbar(c, ax=ax[1])\n", "ax[1].set_title(\"Meridional wind / m/s\")\n", "ax[1].set_ylabel(\"Latitude / deg\")\n", "ax[1].set_xlabel(\"Local time / h\")" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "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.8.5" } }, "nbformat": 4, "nbformat_minor": 4 }