{ "cells": [ { "cell_type": "markdown", "id": "c003e56e-b205-4138-a897-3cc36ac3eb17", "metadata": {}, "source": [ "# Expanding Intervals" ] }, { "cell_type": "code", "execution_count": 53, "id": "d07e22af-5cc4-4e83-8789-ece24e12b337", "metadata": {}, "outputs": [], "source": [ "import pystare\n", "import matplotlib as mpl\n", "import matplotlib.pyplot as plt\n", "import matplotlib.tri as tri\n", "import cartopy.crs as ccrs\n", "import numpy" ] }, { "cell_type": "code", "execution_count": 97, "id": "1d5a8a3d-940a-4cea-b384-8d4d15bac6de", "metadata": {}, "outputs": [], "source": [ "sids = numpy.array([2521902129312956427, 2521902404190863371, 2521899930289700875,\n", " 2521900067728654347, 2521900205167607819, 2521991876949573643,\n", " 2521992701583294475, 2521992839022247947, 2521992976461201419,\n", " 2521997099629805579, 2521889484929237003, 2521991739510620171,\n", " 2521902679068770315, 2521908726382723083, 2521894432731561995,\n", " 2521901579557142539, 2521989952804225035, 2521990365121085451,\n", " 2521991327193759755, 2521997237068759051, 2521997511946665995,\n", " 2521898418461212683, 2521898555900166155, 2521898693339119627,\n", " 2521889759807143947, 2521903091385630731, 2521901716996096011,\n", " 2521897456388538379, 2521902266751909899, 2521897868705398795,\n", " 2521898968217026571, 2521990227682131979, 2521889622368190475,\n", " 2521901167240282123, 2521887973100748811, 2521902816507723787,\n", " 2521892920903073803, 2521894020414701579, 2521890721879818251,\n", " 2521993526217015307, 2521901304679235595, 2521901854435049483,\n", " 2521902953946677259, 2521898006144352267, 2521890859318771723,\n", " 2521886461272260619, 2521888110539702283, 2521898830778073099,\n", " 2521903228824584203, 2521897044071677963, 2521903366263537675,\n", " 2521899105655980043, 2521903503702491147, 2521886048955400203,\n", " 2521991602071666699, 2521899380533886987, 2521890584440864779,\n", " 2521992014388527115, 2521886873589121035, 2521887011028074507,\n", " 2521889347490283531, 2521890447001911307, 2521886736150167563,\n", " 2521893745536794635, 2521893882975748107, 2521894157853655051,\n", " 2521897181510631435, 2521897731266445323, 2521899792850747403,\n", " 2521887698222841867, 2521900892362375179, 2521901442118189067,\n", " 2521906389920514059, 2521908588943769611, 2521901991874002955,\n", " 2521902541629816843, 2521909001260630027, 2521989540487364619,\n", " 2521897318949584907, 2521898143583305739, 2521898281022259211,\n", " 2521897593827491851, 2521993113900154891, 2521993251339108363,\n", " 2521993388778061835, 2521993663655968779, 2521886598711214091,\n", " 2521886186394353675])\n", "\n", "expanded_stare_test = numpy.array([0x22ff8a000000000b, 0x22ff8a200000000b,\n", " 0x22ff8a600000000b, 0x22ff8a800000000b,0x22ff8aa00000000b,\n", " 0x22ff8ac00000000b,0x22ff8ae00000000b,0x22ff8b800000000b,\n", " 0x22ff8bc00000000b,0x22ff8be00000000b,0x22ff8d000000000a,\n", " 0x22ff8e000000000a,0x22ff90400000000b,0x22ff91000000000a,\n", " 0x22ff91a00000000b,0x22ff940000000009,0x22ff96200000000b,\n", " 0x22ff96800000000a,0x22ff97800000000b,0x22ff97c00000000b,\n", " 0x22ff97e00000000b,0x22ff980000000009,0x22ff9c800000000b,\n", " 0x22ff9e800000000b,0x22ff9ea00000000b,0x22ff9ee00000000b,\n", " 0x22ffe8200000000b,0x22ffe8800000000b,0x22ffe8c00000000b,\n", " 0x22ffe8e00000000b,0x22ffe9c00000000b,0x22ffea000000000a,\n", " 0x22ffeb000000000a,0x22ffeb800000000a,0x22ffef000000000b,\n", " 0x22ffef200000000b,0x22ffef600000000b\n", " ])" ] }, { "cell_type": "code", "execution_count": 114, "id": "9a0d02f7-2129-43c0-9d2c-076a61807170", "metadata": {}, "outputs": [], "source": [ "compressed = pystare.to_compressed_range(sids)\n", "expanded = pystare.expand_intervals(compressed, -1, False) \n", "expanded_wrong = pystare.expand_intervals(compressed, -1, True) " ] }, { "cell_type": "code", "execution_count": 115, "id": "cd29402e-9567-46ec-9181-d3a182b82793", "metadata": {}, "outputs": [], "source": [ "def sids2triang(sids):\n", " lons, lats, intmat = pystare.triangulate_indices(sids)\n", " return tri.Triangulation(lons, lats, intmat)" ] }, { "cell_type": "code", "execution_count": 118, "id": "d8acc5fd-7202-46f6-adcd-12d228ebb195", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 1.0, 'Expanded SIDs (correct)')" ] }, "execution_count": 118, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVQAAAD3CAYAAABCbaxBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAB7z0lEQVR4nO2dd3hT1RvHP6d7UWZBZtkbWvbeIEMURWQKslVcP3GwEZAlKIq4GAICgoDiKntvWmbZe88uumeS8/vj3qRJm3RAEdLez/PkaXLvuTcn6Zv3rPd8XyGlRENDQ0Pj8XF42hXQ0NDQyC1oDlVDQ0Mjh9AcqoaGhkYOoTlUDQ0NjRxCc6gaGhoaOYTmUDU0NDRyiFzrUIUQY4UQi3K6bBbuJYUQFXPiXhoa9o4QYqAQYt/Trsd/hV04VPWfckoIES+EuC+E+FEIUSCja6SU06WUQ7Ny/+yUfRyEEDWEEFuEEBFCiEghxFEhRBf1XGshxG2zsruEEIlCiBghRLRadrQQwvVJ11Mj+wgh+gohjgghYoUQ94QQG4UQzZ92vZ51hBBDhBDnVTt/IITYIITIp55bKoSYqj4vq3ZWYtXHAyFEgBCiw9P9BJY88w5VCPER8AXwCZAfaAz4AluFEC42rnH672qYLf4FtgLPAUWB94HoDMq/K6XMBxQHPgJ6AxuEEOJJV1Qj6wghRgLfANOBYkAZ4Aeg21Os1rP8OwBACNEK5Tvro9p5NWB1JpcVkFJ6AX4ov6U/hRADn2hFs4OU8pl9AN5ALNAzzXEvIBQYrL6eBPwOrEBxUEPVYyvMrhkA3ADCgQnAdaC92fUr1OdlAQm8AdwEwoBxZvdpCBwEIoF7wHeAi9l5CVS08lmKqOcK2PisrYHbZq93AUPTlCkDxANdzepyRP3MD4A5T/t/ltceKI18LPBaBmVcURzuXfXxDeBq/n8HPgVCVJt6GegCXAQigLFm9zLa+mogBjgG+Jmdvw6MAk4CSYATSifkgGqzwUBrs/IDgavqva4B/dTjFYHdQJT6G1htdk1VFGcWAVzA7PcJFAb+UW0yCPgc2Gfje/kY+CuD720pMFV9bvxdOlm5xwPAQX09Crijfp4LQLv/1B6etkFmYqydAF3aL1E99wuwyszIUlRDdADcsXSS1VWjbw64AF+q5TNyqAvV+/iphllNPV9PNVAntew54H9m9bLlUAVwCQhQ61kszfnWZOJQ1eN7gC/U5weB/upzL6Dx0/6f5bVHRjZqVmYKcAhlVOKD4tw+N/u/64CJgDMwDKWzsBLIB9QAEoByaWy9h1r+YxRH6Kyevw6cAEqr9lsSpRPRRf1tdFBf+wCeKI6vinptcaCG+nwVME69xg1orh73BG4Bg9TfQB0Uh1tdPf8bsEYtVxPFudlyqC3UzzYZaIbayJidX0rmDrW8erwaUEWtWwmzayr8l/bwrA/5iwBhUkqdlXP31PNGDkop/5JSGqSUCWnK9gD+lVLuk1ImoxhvZiIGk6WUCVLKYJRW3Q9ASnlUSnlISqmTUl4H5gOtMvsgUvkPt0Ex+K+Ae0KIPUKISpldm4a7QCH1eQpQUQhRREoZK6U8lM17aTw+hbFto0b6AVOklCFSylAUB9Lf7HwKME1KmYLikIoAc6WUMVLKM8BZVPtTOSql/F0tPwfF4TU2O/+tlPKW+jt4Hdggpdyg/ja2ooxquqhlDUBNIYS7lPKe+n7GOvmiOKdEKaVxYakrcF1KuUT9DRwH/gBeE0I4Aq8CE6WUcVLK0ygdH6tIKfcC3YG6wHogXAgxR71PVrmr/i0E6FFGA9WFEM5SyutSyivZuNdj86w71DCgiI25oOLqeSO3MrhPCfPzUsp4lFY6I+6bPY9H6QEihKisTobfF0JEo8wBFbF2g7RIKW9LKd+VUlZAMdY4YFlWrjWjJMpQC2AIUBk4L4Q4LIToms17aTw+4di2USMlUKabjNxQj5nuIaXUq8+NnYEHZucTUO1PxdyWDShTBiWsnUexs9fURdBIIUQkykituJQyDugFvIXSwK8XQlRVr/sUZVQVJIQ4I4QYbHa/Rmnu1w9lXcAHpddq/v7mnzsdUsqNUsoXURxiN5QpiOwsEJdU/0ZIKS8D/0PpxYcIIX4TQpSwdeGT4Fl3qAdRhtvdzQ8KIbyAzsB2s8MZ9TjvAaXMrndH6Vk8Cj8C54FKUkpvYCyK4WULKeUt4HuUYVGWEEKURply2Kve45KUsg/KUPIL4HchhGd266LxWBht9OUMytxFcURGypDas3oUShufCCEcUGzb/H7mv4VbwHIpZQGzh6eUciaAlHKzlLIDSgflPMpUF1LK+1LKYVLKEsCbwA9qOOAtYHea+3lJKd9GmarQmddP/ayZovaetwM7yMZvAngFZe75gnqflVLK5ijft0T5XfxnPNMOVUoZhTI8mieE6CSEcBZClEWZo7kNLM/irX4HXhRCNFUjAybxCE5QJR/KvFOs2pq/nZWLhBAFhRCThRAVhRAOQogiwGCUubXMrvVQV0T/Rpno36Aef10I4aP2UiLV4oZsfyKNR0a10YnA90KIl9X/lbMQorMQYpZabBUwXgjho/7fJ6IsoD4q9YQQ3dVe8f9QHLotO1qBYvsdhRCOQgg3NUSvlBCimBCim9oIJ6GsMxgAhBCvCSGMnZCHKM7JgLIGUFkI0V/9nM5CiAZCiGpqL3sdMEn9HqqjLO5aRX3v3upvQwghGqJMn2XlN1FMCPEu8BkwRkppEEJUEUK0VUMLE1F69v/p7+GZdqgAUspZKL3AL1EcWSBKK9lOSpmUxXucAd5DmZ+6h2I4IShGlF0+BvqirCIuJPMwDyPJKJPk21A+x2n1/QdmcM13QogYlOHfNyhzVZ1UBwrKgsgZIUQsMBfobWX+WOMJI6X8ChgJjEfppd0C3gX+UotMRZm3PAmcQlmZn/oYb/k3ylD9IcpcbHd1PtVa3W6hDKXHmtXtE5TfvoNa77so00itSO0gNAACVdv6B/hASnlVShkDPI8SwncXZWrsC5S5S9TP7aUeXwosyeBzPERZhLuE8ptYAcyWUv6awTWRQog4lO+xC0p0xWL1nCswE2Uq8D7KyG1MBvfKcYS6GpanUKcMIlGG7deecnU0NLKMEGISShTJ60+7LhrpeeZ7qDmFEOJFdRjiidLbPYWy4q6hoaGRI+QZh4oy7DEGVldCGR7nve65hobGEyNPDvk1NDQ0ngR5qYeqoaGh8UR5psUTzOnUqZMMCwvLvGAe5ejRo5ullJ2edj00sodm1xljb3ZtNw41LCyMI0eOPO1qPLOo8Y0adoZm1xljb3atDfk1NDQ0cgjNoWpoaGjkEHYz5M8qey+FUt7Hi4RkHXFJeuKSdMQlG//qlL/q8ZCYJP4JvsvSQQ1oXaXo0666hoZNDAbJtnMPqFbcO50dG+07NklHvGr3Sw9cB+DajC5oeuT/HbnKoSbp9PT/OQiAMoU8yO/ujIeLI16uTni6OuHp6oinixMerk4U9XZl0T5lk9SYdadYN6IpxfO7P83qa2jYZMHeq8zceB6AysW8FHt2SbVpT1cnPFwd8XJxwskhdeA5+d+zfPZidc2p/kfkKoc6Z+tF6pQpwPGbkdyMiOfi1M64OFmf1VgVdBOAdSOacuR6BL3mH+LXoY0oXcjjv6yyhkamxCbp+DXwhsm2m5QvzORutgWZyo5eD8DJSc8zcHEQY/88zbSXa+LgoDnVJ02umUM9fD2CP4/dYdGA+hwZ3x6AyuM3kpiiT1c2Kj6FMetO0bZqUeqWKcjwlhUY0rwcvRcc4npY3H9ddQ2NDJkacJam5Yvw54hmTOhanV8O3uCTtcFWy87ZcgGAHR+1wtvNmWVDGnElNJaPfw9Gb9A28TxpcoVDjU3S8dGaYKa/UovCXq4U8XIleOLzAFSdsIm4JEsxdb8pWwD4+Y36pmNvNC3Lu20r0nvBIS6HxP53ldfQyIDt5x6w/0oY47tWA2BI83LM6F6LtUdv8/aKoxZl70Ym8O2OywxsWpbyPooetZerE78MakhIdBL/W32CFL2m7vgkyRUOddr6szQpX5j21YuZjuX3cOb05I4A1PhsM1EJirrZ11svAkoLnnZeqU/DMnzSsQp9Fx7i/P2MkpFqaDx5wmOTGLPuFF/28COfm7PpeJ+GZfimlz8bT9+n/8+BpuNNZ+4AYNJLNSzu4+7iyKI36hObmMK7K4+RrNOc6pPC7h3qjvMP2HsptQU3x8vViXNTlE0WfpO3cOZuFHO3X7JowdPyar1STOhanT4LDvH9zstPtO4aGraQUjLuz9O8XKckjcqnTy7xcp2S/PR6XfZeCuPl7/fz0RplCuDwuPZW7+fm7Mj8/sqIrPuP+7kflfjkKp+HsRtxlPr168u0O0rCYpN44du9zO1dh1ol8xMSk0RIdCIhMUk8iE4kNCaJkJgkbj+M5/D1h4/0vmvfakKDsoUyL/iUEUIclVLWz7ykxrOENbsG+P3obRbsucLf7zQnLllHSHQSITGKbYea2fnG0/et3DVrXJ/5wuNU/T/B3uzarlf560/dBkDvBdlL9ukgwK90AYrmc6VoPjflr7crRb2V50k6A91/OECZQh68tfwoP/WvZxdOVSN3cDM8no/VRadqEzdl69qKRb1Uu061Zx+jnXu7cvhaBD/uvkLAe82fRNXzPHbtUOf3r8eby48yqFlZJna1HWt34X4MHb/Zw6edqjBr0wUMEr7p5Y9vYev57IxhJ3s+bcPeS6G8tfwo8/rUoWlFu9pWrGGnlCnswUt+Jfgn+C5/vN2Uer4FbZbt+dNBgq5HsGxwQwYsDuJuZALbRlrPah4SnciXWy6w6I0GFnOyGjmHXc+hdqzxHEfGt+fglXBmbjqPtekLKSUdv9kDwIjWFbk6XUlH3mr2Li6HxKQr/+7KYwAcVUOvWlTy4ft+dXlv1XF2Xwx9Uh9FQ8OCb/vUUTsMRzh20/p01cEr4QRdj2Bub39aVvYh4L3mxCfrKTt6fbrfgpSST/84Sd+GZfAvXeA/+AR5E7t2qABFvFxZNawxBy6HM/nfs+kMyTgdYFzxd3AQXJuhONX2c/Zw9m7qav7pO1EEnLzHZy9Wp7CXq+l44/KFWTCgHiNXn2DbWfN06RoaT46ONZ5jdg8/hv1yhMCr4RbndHoDfRYeorCnC938ldT0NUvmZ/P/WgJQbswGi9/CqqBbhMUm8V67SkgpKTt6PVdCtfDAnMbuHSpAQU8XVgxtRPDtSMb+eRqDGsAceDWcwGsRfNPLHy/X1NkNIQTXZ76Ai5MDXb7dy4lbkUgp6TpvHwCDmpVL9x71fAuxeGADRq87ycZT97JVv9sP461uMNDQyIw2VYvybZ86vP3rMfZdStVNbf7FTgAOjW1nUb7Kc/nY+XFrQHGqBoPkRngcX265wNc9/XF2dDAF+OeETUYnplgdGeZV7HqVPy2xSToGLzlM6UIeTO9ekyrjN1HI04VjEzrYvu/UrYTFJpten52i9GRDopMIjU2yWF0NiU5i98UQU/n321YkNklPfLJRmCJVhCU+Sa+KWOiJVTcWXJrWGWfHJ9OG2dtqqIZCVuwalM7BiF+P8eVrfsQl63h35XF+HdqIZjbm9W8/jDc5Xb9S+XnJvySDm5UlMcVgElEx2mxsUqq9fvr7SQDebl3BJMBiy77D45Tfwa6PW1O2iPX1iMfF3uw6VzlUgPhkHcOXHWXfZaU19ytdgIdxyYTEJJKYkvMBzeNfqIaHi7nwiiLG4uHihJerE2+tOMqJW5E0KFuQAh4ufNe3Dq5OjjleD3szPA2FrNo1wLGbD3l9USDxyXr8Sxfgr3eaZVj+QXQijaZvB8DFyQGd3oCzo4Nin2aiKorQiiOB1yKIUJ1ky8o+tK3ig4e5CIurk0lsyN3FkTeXH6V9tWK806bi430JGWBvdm3Xq/zW8HBxYlSnquz7Thm+B9+KNJ3L5+aUGirlrYSWFPBwYfbmC6YyX73mR/e6JW1GDHT8eg8XHsTwXd86vLvyOH+fuMu/NkJQjt6I4MStSGb1qM3L/iV5b9Ux3lp+lB9fr4ebc847VY3cTd0yBfHJ58qN8Hj+HNE00/LFvN34c0RTXvnhAMk6A2cmd8TT1fpPPjFFT9UJm6jg44nOINlzMZShzcvRsrKP1fKL911DAG+2LP84Hyn3IaW0i0e9evVkVkhK0csuc/fIhXuuSN9RAdJ3VICMjE+2Wd5v8mbpOypApuj0stf8A9J3VID858Qdq2V3nHsgfUcFyC1n7ksppdx9IUT6jgqQrWbtSFdWpzdI31EBsvqEjaZjyTq9fHflMdlv4SEZn6TL0ufJKsAR+Qz8n7THk7FrKaU8eiNC1vt8q3wQnZDla6SUMjIu2fRbiEtKsVqm0rgN0ndUgNTpDVJKKbvM3WNh6+ZcehAt60zZIq+FxmarHo+Cvdl1rliUMue7HZco5u3GkOblOP956rZT41DGnD+P3yYyPoXVwxvj5OjAb8Ob0L5aMd5bdZw1R25ZlE3S6Rm09DDli3jSQdUMaFnZh9XDG3M9PJ56n2+1KN/my10AHJuYOn/r7OjAN738KertyhtLgkxzqxoamRGfrAgAfd6tBkXzuWXrWnNdi+oTNxOdmGJxflXQTZJ1BtaNaIqjKvG3/v0WNCxbiGHLjhBw8q6pbIrewIerg/no+cqULeLJhlP3aDJj+2N+utxDrnKox28+ZGXQLWa+WgshBG7OjlyYqjjVup9vJSQmdf9ybJKOD1cH06hcIYu90oveqM/L/iX49PeTLDt43XTcb7KiULU1TdB0o/KF+XNEU8Ljkqk8biMAm07f52ZEPEsHNUg3X+roIPiyhx/li3gy4OfAdMatoWGNGRvO41+6AJ1rFX+k6811LWpP2sJDtYORVsrSnDVvNaFt1aK8u/I4a9UOxrwdlyns5ULfhmUAOHErknuaLoCJXLMolZCsp8u3e/mkYxW6pDG6FL2BSqqzOzC6LSUKuJt2Q9lKETHuz1P8GniT0Z2rUtDDmVF/nMpw18qZu1G88O0+0+tqxb3Z+EELknUGQmMt91+HxCTxICqR1aqRvt+uEiM7VM7eF5IGe5u811DIyqLU7ouhjPnjJBv/15L87o+3wylJp6fKeGU7a9C4djScpvQuM0qV8t6q4/wbfJfX6pVi54UQNrzfAp98riTrDaY0LKZoATU9S2R8ChtP32fpoAaPlS3A3uw61zjUz/4+TVRCCt/0rmP1vN4gqTB2A6DIn60KusmWD1tSuVg+UxkpJZHxKarTS+R/v50whYYANCpXyCS4ktPD9aPj21tsJsgu9mZ4GgqZ2XVUfAqd5u5hdg8/mld69K3PKXoDYWoY4N3IBN7+9ZjpXHkfT7xcnUwhgpnpUOd3dyYuSYcQmFKxeLg4WqQZ2qJugHncnFb2Zte5YpW/94KDHLoawScdqzBn60WznmCiaiRJFuWN6U+e/3pPtt4n8FqE6bmHi6MpYsBHjRjYcuYBdyITTGWCxrWjiKer1dQT8ck6qk/cjF/pAjSvWJg+Cw+xYmijbM+PaeReElP0+E3ZQiFPF87fj2bPpVCLUU5IdCLRiY/fsF8NtcxSUcjTxSSq4mMmILT5zH2+eLU2+d2d8XR1sple6OiNCI7djGTjBy3yXC4ru3eoUfEpHLqqODrz8Ke0FPFywSefG+fupW417dOwNK0q++CTL1WVxzyc6eutF5m7/RI965dizZHb9GtUhmmv1LJ6/4dxySzZf50utZ5j3AvVaTZzBw2nbefytM44kN6oqk/cDMBfI5oihMDVyZHeCw6xcmhjnsuvOVUNZa4TICIumanrz1kt4+LooDg9tVG3UE/L54ZPPleKebtRyNPFtOAEqQJARraNbEXFotY1go0Mbp5+B2Fa4pJ0jFwTzNSXa+CT79FHXPaK3TvUz/45TduqRdlxPgSAy9M642RjN9LEv09z7l40B8e0pdWsXawKukWvBtbFIu5FJTB3+yUGNPFlSrea1CiRn8/+OUN8sp6ve/mnK19HXeX/oV89AILGtqPh9O1UHLcxXbLARXuvArDh/dQW/P12lXB2dKDXgoOsHNaYkgW0DKx5mZDoRLzdnahcLB8HroTzdusKjOpUNUfubZQGDBrXDh8vV8qN2UD7ObvZ+EELqhX3fqx7T99wjnq+BelU89EWz+wdu17lX3/yHidvR/F937oEqnuaK47baDXFw5XQWJYdvMGH7StTPL87F6d1ppCnCy9/vz+d8ARAkxlKOokpanbJN5qW5YtXa/Hn8TsMW2Y55/XFJiW9755P2piOFfV2MylWmScLDItNYur6c7xSpyTVS1ga79utKzCgSVl6zT/IzfD4R/pONOwfKVVlqEa+rBzWmE86VuHHXVf47O/Tj33vC/dj+P3obcZ0rkrRfG4IIdR5Tug8d6/FRpjssvNCCLsuhKZLwZKXsFuHGhKdyGf/nGZOL3/cXRwpZsOBgWKg7b7aDcAH7SuZjh+b0AHfwh70WnCIPWbSfJ+YteDm9GpQhrm9/dl69gF9VBWrWxHx/LjrCm+2LE+ZwpYpqAtbSRZoFMW21ssFJQnbm60q0HvBQa5qakB5EpMyVFtlS+c7bSoyMZNsp1lBylQpyzdbVTAdV5zqC+Rzc6Lb9/s5fD3C1i1s8jAumdF/nGR2j9p452GtVbtc5ZdSMnjpYWqVzM/I56tYlIuKTzFlNTVutRu89DA7zocQ/NnzVsNOus7by+k70SwcUJ8yhTzo+M0eRneuyltmRmfO5jP3eXP5UWqXys/J21GA9XQSsUk6QqITuR4ex+Clqb3aRuUK4eLkYFpVfRhvOxZ1/AvVGNoi8+199rYaqqGQdpX/RngcL3+/nzVvNqGSWQQKKIupY9adokut50xTS9mh1/yDBF6L4PTkjhbqa+Y0m7mDO5EJGQqvpEVKyburjlMsnxsTX6ye7XplhL3ZtV061FVBN/k18Abr3m5mdaUxNklHzc+URR+jkvkXr9aiV4MyNu/f48cDHLmRKuQ7o3stS6WpmCRC1RVWXQ7nN3cQmBYQjAsKq4KUGNUCHs6sGtY407ktezM8DQVzu9YbJL3mH6RTzedsNqJ/Hb/D/1afoEWlIiwf0ijL73Poaji9Fxzim17+vFynZIZljXoViwfWp23VYhmWBfj7xB3m7bhMwHvNc1yjwt7s2u4c6uWQWF776QBLBjUkv7uzZRiJKUwqkZsR8dyKSMj8xo9AQQ/ndL3KCV2rpzpENZ+PeS8g7arq4XHtba6Cjv3zFCsDb3JoTDuO3Ihg0j9nWTKwAbVK5bdZJ3szPA0Fc4f6464r7L4YwqI3GpCgSuWZy+YZA+f/OHaHIDWELyuJ9nR6AxXHbaSghzPH1SmozOj+w36O3Yzkh351022UMedeVAJdv93H0kENrdrn5H/PMKZzNZshVplhb3Ztdw41rWPKKiULuFOjhLcpnMQ8tKSotyuXH8TSd1FqjvPPu9Wgf5OyVu8lpaTcGGWTwLoRTen+wwHcnB04/3lnq+WH/nKEbeceEDzxeTxcHU27tg6NaZcuROpySAzt5+zhk45VTLJom8/cZ9yfp1gwoH667YFG7M3wNBSMdm3Mewbg5uxgktYzyuWZB827OTuy/NANACoV9Uq3HTotTWds525UYrpok8zo/3Mgey+FMaenH93rlkp3XkrJgMVBNChbiPfbVUp33rhD8ZfBDWllQ7UqM+zOrp+2OktWH0ZVnl0XQmSdKVvkdzsupRWmsWDb2fvSd1SA3Hjqnklp587DeKtlU3R66TsqQPpP3iyllHL0Hyel76gAuWD3FavlX5q310K559TtSNN7GAwGi7LHbkRI31EB8tdDN0zHjEpUvqMC5M3wONNxgyH1eFp2nHsg607ZIoOuhVutE3amyqM9LO1arzfI62GxMkWnT/e/TatMNnfbRfn6okNy/6VQ6TsqQDaati3dNUbWn7wrfUcFyD0XQ2yWyYi3lh9JZ79GfjlwTb703T6rdc4p7M2u7W6Vv1VlHzZ90II/j99hzpYLKN+5JUk6PUN+OUJ5H0861XyOK2pivqYzd3AjPC5d+ZazFGXzwLFKlMCM7rUY2rwc0zacY+62SxZl914KJfh2FD/0q4uHizKkr1kyP1s+TJ/Lx2CQvPLDAZwdBX0bpc7fOjoIU7LAFrN2mlbzBywOAuDUpPTDsjZVi/JNb3/eWn6UA5fD0p3XsG8cHAS+hT3TxVAn6wxUm7iJXReUOOtTt6P45cB1ZvWoTdOKRfjj7abcj06kupV00wnJekb8egy/UvlpUenReog/vl6PV+uWYuyfp1i875rp+JXQWL7Zdomve/rZjPvOi9jlN1HU243fhjdmy9kHzNyYPtupSRnqQ2UoZO7AlGynqeFI60/e425UIsuHNLQYDo3vWp332lbk620XmbFB2aWSrDPQ/+cgSuR3SzevVLlY+lw+xiHcqUkd030G82SBbb/aza+BN9h7KYyvXvOzmuJXpzdQqWg+3mpVgb6LAuk5/2DWvzANu8XFyYF1I5rSspIPiSl6PlxzgokvVqd4fmXjRz3fgjaznVZTnWxmyv6Z8VVPP2WDS8BZvt95GZ3ewMg1wfyvfSXK+2S8uyorJOsMRMYnc/thPHsvhdp1uKDdzaGa8zAumQGLg6jnW5DPXqyOEILVh28y6o9TrBvRNN18o5Spc58bP2iBb2EPZT99qfz8/a511f0fdl1m1qYL9G/sy78n7xIZn5LhbqxbEfG0UHu8AAOa+FKrZH5TfqrQmNTIgdCYJOKTHz1R2tXpXUw6AXY316QBZC8FyucBZ7kflch3feuk2yN/8UGMSZvi2owuLD1wncn/niXgvebULGl7MTM7zNx4np92XwGUudsf+tUlIcUyJ5UxB1Wc6XVaJapUdSrjeSkVkRV3Z0fuRydSu1R+/lF/j/Zm13btUAGiElIYuCSIqs95M6pTFfynbKVNFR+WDGpo816Vx20kWZ+6m+rslI4mFamQ6CQemImrhMYksffSkxtie7g4WjhVv9IF6FijGD5eSqSAUWOgoIeyF/tmeDwv/7CfVcMaU+W51DhFezM8DYWsOtQDV8L4cPUJNn3QkoKeLlbLXAuLMwmbA3StXZzv+tbNqaoCMHLNCdYduwNA+SKepoUz89xUSj41RyUflcUx81xWqTmtjCPD+buvsON8CKuGNbbbjoLdO1RIzXYapO7wKFXQndCYJJKsbEHNCWqXyq86ulSHZwyVKprPleZf7LCQQDv/eSeb8XkbTt1jxK/H+GVwQ0b9fpL70YmsHNaIphXSB1XrDZI+Cw7RvnpRhre03HRgb4anoZAVh7r93AOG/KKU8cnnSmga9TRbZCWk6lGYu+0SX2+7SJ+GpZnRvXaO3PP8/Wj6Lgzk73eaUbpQ6o5De7NruxdHAUWNfMGAevhPUQRKbj9U4k/zuTqZpPXMQ6WcHByYEnDWdP2ywQ1tJiOD1BjSub39+eC3E+R3d2bRGw2sll19+CYGqYRTlS/iif+UrVSdsImzUzqaFrGMGBcNapb0plVlHw6NbUf7ObvpuzCQpYMa0LpKUYvyi/ZeBQFDmmuJ0fISi/amLgYZnWkRL8Wui3lbJp38ed81rofH82O/nO2ZmvNB+0oU8nJhwl+niUvS820f6xrEWSVJp+fD1cGM7lzVwpnaJU87zCCrj6wkM0tI1smBiwPl8GWHZVKK7VAO89Ckzt/YTkYmpZQ/7rosfUcFyAv3o6WU0hSG1e27fenKRsYrydAGLg40HYtOSE2QFp1gmSzQeFyvtwy1eum7fdJ3VIDcdPqe6djZu1GyzpQtFmFW5mBn4SXaI+t2LaWUJ29Fynqfb7WZQFJKKW+Gx0nfUQFy+oazWbrn47Lm8E3pOypADl4S9Fj3mbnxnBz6y+F0IYdS2p9d2+Uqvy3cnB35qX89pIS3Vhy1EEgxMmOjsmK/91NFGWrDB6nJyP4NvmtRNiQmkZkbz9OzfimTsn+nms+xdFADTtyKpMOc3RbljdEFiwem9l7zuTlzdoqyyl9r0hYi45UMAEv3K72OgPeapxOg/vudZjQpX5g3lx/ln+C7agt+Ine04BqPRK1S+Vk+pCGfB5xl3bHbVssYF0PHdK72n9TptfqlmdenDtvPh9B7waNFnRy5HsHvR28zo3utXCFGnascKoCrkyPf96uLu4sjw5YdIcFswedWRDzzd1/lzZblLRyTMRnZe6tSk5EBpnw7s3r4WbxH6ypFWTWsMZdCYmmqZnycs/UiADs+apXOMDxcnEwZWP2nbOXSgxgm/XuWF2oXt7kCu2p4Y9pXK8b7q45TbcImShX04LV66XeraOQdqhX3ZuWwRszadIHf1KwTRt5RU5ocm9DB2qVPjBf9SrBwQH0OXY3gxXn7Mr/AjLgkHR+tDWbayzUp8hjpf54lcsWilDV0egOf/n6SO5EJLB7YAE9XJ9NcqK3JemMyss9frsmlBzEsO3iDg2PammL+0nL0RgSv/pjaMg9sWpaxXaopuXvMkvIZIwjuRCZYyARmh72ftsmwd2pvk/caCtm1a1BW819fFMibrcozoElZTt+Jouu8fUx6sToDm2Wuqv8k2HMxlAGLgyhXxNMUj50ZY9adQqc3MPs1P5tl7M2uc61DBWVVfNyfp7gUEounqxN7LoZybEIHCqlhJ1JaJuULiU7iozR6k+WKeBISnUjcY8SLZkSVYvko6u2Kj5dZDh9vN3y8XPFydWLgkiDeb1eJN5qWzfA+9mZ4GgqPYtegjLb6LjrEgMZlmaZuPHlSq/pZxahoVcTLlSOqNrEtdp4PYcLfp9n4QQurG1mM2Jtd52qHCsr2z/JqttPHwd3Z0SJvjzGPz6xNlnmszIPt02IM+h/esjwL9ihpUHZ93JqyRTytlh+z7hQpegNfZtCCG7E3w9NQeFS7BrgbmUDTmUpmiXNTOuHukrPSeY/C8ZsPeSUTsaCIuGQ6z93D3N51aFy+cIb3sze7zhVhUxkhBDQoW5DD1xWt0486VKZpxSJWk/IBJjHqN5r48svBG7zVqgKjO1vP5XPhfgyzNl1gbJeqtK9WjLZf7ab82A02napx0WBsl2qM6VyVcmM20PrLXWz9sGU6MeEd5x+w91IoGz9okRNfg0YuRKdXOkOfvVj9mXCmAHXKKFthu87bR9nR69OlkZZSMv6vU7zkVyJTZ2qP5LpFqbSsPnyLhBQ9F6cqreVXWy/i7uxI6UIe6ZzpsZsP2XE+hJndazG5W00+7VSFn3ZfYaKVXD5Spu7VH96yAuV9vEyRA+XHbkCfRoT67RVHATiuLhoYc/kAdPh6D2fuRpnKRsQlM2bdKb60sa9fQ0NvkHy09gTjX6jGoKc0b2oLW2JBAH+fuMvlkFg+SpNpI7eQqx3qzfB4Zm2+wJye/rg4OXB95gu4ODrQ5du9nEiTjExvkHT/4QAuTg70bqgoQ41oXZFJL1Zn2cEbfLTGcm6113wlp9SZyanCJ6ULeXBwTFsAKozdQIq6vfXU7Sg2nr7PlG41LLYNCiGUOjk58MK3+zhxK9LUgr9YO3e24Bo5w8K9V3EQgsHPmDM1Yk0s6G5kAp8HnGVOT/8cV/Z/Vsi1DtXYgo9oXcEUQwpwcVpnCno48/L3+02q54Cpt3nyM0vpvIHNyjGzey3+OHabt5YrvcyDV8IJuh7B3N7+eKbJzVM8vzuHxykT8pXGKckCX/xOCScZYEOw+uLUzhRWM7CO/fMUl0Ni+bhj7mzBNR6fc/eiWbDnKl++5mdzvv5ZoFwRT/aNSh21fbj6BIObl8sxsZZnkVw7h7oogxb8+MTnaTFrBz3nH2T5kIYkpRi4HBLLz2/Ut9py9m5YBjdnR/63+gR9Fx7iwJVwCnu60M0/NTePlFJJyqcKrEzpVoOJf5+h6gRFQq2+b0E6fr0n06R8q4JuMfmlGrm2Bdd4PJ7GJg+DQRKZkGKWYig1KiY76mmB1yKylQfLHsmVDvX8/Wjm77nK3+80s9mC7/yoNY1nbKf/z0GmY6fuRLH9fIhqKKm5qsznQw9cCQcgPC45W+lYzBMAOjkI06KYj7oPu7CnC/N2XAbgp91XaF6pCBVyQGtSI3fxzbZLlC6U8SaPxBS9YsOx6R2gKeFkTCJhsck5Xj8PF8fUSBg1KubErUiO34wE4JPfg/nqtdwrSp3rwqbuRSXQZMaOHH3vgh7O+ORz5eIDS+Hb7/rWsZDZMx/+J+sMVB6/kSJeLibDDZ74PPk9rC8ytflyF9fC4rgwtRN/H7/LV1svsHxII4vpioywt/ASDYXshE3V+3wr4XE55wSFwDIMUH34qPZsVFAr4uWCq9OjjZgSU/RUnbCJqs/l4693mjF8+VG8XB2Z27sOzllwqvZm17muh+rhnPqRWlb2wbeQh0VCPp98rhTzdqOQpwu/H73FqD9Okc/NiZhEHbN71Oa1+qWt3jc+WUf1iZvxL12ATztWoe+iQGZsOM/+0W2tlq83VVG+ChzbnhS9gaoTNuE3ZYvFxgIjW87c51pYHEsGNcDVyZGeDUrj7CTotyiQXwY1pHqJjFNIa+QN6voWZOvZBwB08y+RqqKmbgwp6q2MeLzdnJ6ZffHGKa8N77fAwUGwoH893l15jLdXHOP7fnUe2VE/szxtdZasPrKqyiOlkgTvw9+Oy54/HZCxiSlWyxiVoQapSjnvrjwmfUcFyGUHr1stnzYJX9C1cOk7KkDW/GxTurJrj9ySvqMC5GGzhHpJKXrTPR5EJZiOJyTrpO+oANnx693p7hMQfFfW+3yLDL71MNPPjJ2p8miP7Nu1lFJeDY2VTWdsl78cuJat654Gyw5el76jAtLZb1KKXr61/Ih8Y3GgTEjWWb9Yxd7sOldOZDg6CGa/5ke5Ip4MWBxEdGL6RSCjMtTPbyijiXl96tCnYWkm/HWaheouJiPG5GTr329uavkblC3EP+82IyZRZ5HLJyYxhY/XBtO8YhHqly1kuoeLkwOXpymxsA2nb+dupKLZakyutuH99AH8L9QuzvRXajFoyWGO3XyY7rxG3qNcEU9+G96YhXuvKvq4zygP45KZ8NdpOtYoRu1SBSzOuTg5MK9PHfK5OTPkl8PEJ+ueTiWfALnSoYLiVKe/UosaJbx5fVGgSTYPMGUyTasMNaN7bQY3U7KdfrtdKRMem8SUgLO85FeCGiUswz1qlypg2slkDGCuNUlx1MuHpE/B4uToYJGBdfqGcxgkVhfPknUG7kYmUNTbjY41n6P7Dwd49ccDj/u1aOQCShfyYPXwJqw4dIPvd15+2tWxSp3PlSmv+f0tpz+llMQkpnAzIp4+DUpz+NpDqk/c/MiiQc8auW5RKi1SSqatP8eBK+EsH9KQZL2BJjN2KFkcu9W0es2Xmy/w3c7LvNWqgikpWUbCE5dDYmg/Z4/p9Zsty+Pt7pxObSokJpHElMdLy2KrHvY2ea+h8Dh7+UOiE+mz8BAv1C7Bh+0rPTPzptmJfjHHmm3bm13nukWptAghGPdCNb7ccoE+Cw+ZVuoHNytH0LUIUyI+Y/yo8TVgcqaQPSOZv8f2UMzL1QmffK5cC4szHXvZvwQtKvmYBFeK5nOjgLuzqdd69EYEby4/pu3r17CgqLcbq99swuuLAknWGRjVqcpTd6pRCZbTa5+9WJ1i3spisFFRzTwaZuyfp0hKMfBVz8wFgOyBXO9QQXGqn3Ssyvc7Ux1ka7PskFmhfbViFmpT5pEDhb1cqDRuo0X5S9M62wwLmbXpPD/susKuj1vTbs5u/jpxl0HNyuFXukC6snFJOkauCWbqyzXxyZc7RHg1co4iXq6sGtaY/osDSQ4wMKFrtafqVCf/c4b+jX0Z2aEyr/8cyM2IeAY2LWu1TjsvhLD7Qigb/5eLOgpPe1Usq4/sroZaIzFFyTnV4osd8laE9dxMRowr8r8euiF9RwXIt1ccsVl2zLqT0ndUgLwXmSBDohNN1yampF/BNOb9mfLvGdOxWp9tShcVYGTsupPyw9XHM/1s2NlqqPbIObuWUola6fbdPjnuz5PpcpT9V2w4eVe2mrVDxiUpkTWRccnyJRt1iohNkg2nbZX7L4dmeE97s+tcuyhlDVcnR5YMasigZmXpveAQN8LjrJYb8auyZ//YhA70bVSGOT392HDqPgMWB6UrezkkhpWBN/n4+co8l18Z2hgVpaqM35Qur5VRwm9C1+qmYycndaRkAXd6/HSQA5fDTMd3Xghh14VQJr1U4/E+uEauJ7+7M8uHNOTC/RhGrzuZTu3sSRMSk8iEv88wp5c/Hi5OHLv5EHcXR1YMacj5e5Z1klIy/q/TdK1dwmq6dHsmTzlUI4OalePt1hXoveAQV0Itdz+dvhPFhlP3mfxSDVMAfve6pfihX132XAylh9lKu5TStBj1bttKpuMFPV0IVkVWqk7YRFySEhbywW/HAThqRc18/+i2VCmWj76LAtl5PoSHccmM+eMUs3vUxluT8NPIAvncnFk6qCG3IhL4eG0wOv3jLYBmFSklY/44Re8GpalbpiAA3X84wOh1J8nn5swvgxtyMyLeVKd/gu9y4UEMn+RCAaA86VAB+jXy5cMOlem78BAXH8QAighEVzXRWNqUI11qFWfxwPocufGQTqoylbHHenKSpUIVKD0Go7Rfjc82c+hqOH+fuMuErtUpbCMh2eYPW1LPtyCDlh6mzudb6VKrOE0r5q4WXOPJ4unqxJJBDQiPS+aD306YJCSfJKsP3+J+dCLvt0vtVJz/vBNzevqn1mlgQ0Jjkmj95S4+WXuSr3OphF+uD5vKjL9P3GHq+nMsHdSA0X+c4tSdqAzTSey/HEa/RYGm19NeqUnrKkVNIVIhMUmERicSGqtEDdyMiOdSSKzVe2XGX+80w9/KQpU17C28REPhSdl1Yoqed349hoOD4Lu+T26L55L915j871maVyxCiQJuxCXriU/SEZekJy5ZR3yyntgknXLMTIkqq/mv7M2u87xDBdhw6h4j1DS8oIQ2xSY9md0bRbxcqFbc21KQwrQX2427kQn0WxRIrZL5uReVyNJBDbKkH2lvhqeh8CTtOlln4P1Vx0nU6fnp9Xo53iPUGyQV1HxtZQt78FarCni4OuHp4oinqxOeLk54uDri5eqEh4sja4/c5o9jt/nj7aZZrou92XWeCJvKjI41nrN4HZuksypDZgyXMs+M6iDgyvQuNkNVjD3aub39+eC3E4TFJjO7hx/P5XdLV9ZgkLT5chfOjoJ/32vOptP3GLgkiEVvNMhyT1VDw4iLkwPf9a3DyDXBDP3lCAsG1MPDJed+8ov2XqVRuULM6lGb138OJDZJZ8p2kZYrobHM23EpW87UHsmzc6jmLNhzlcblCzFQnTf9qENlzk7pxK5P2rDmrSZ837cun71Yg7dbV8DLTTHIJYMasPGDFhhk+rw5RlL0BvotCqSYtyvd/Euatp02nrGdWxHx6cobswacmqTMvXaqWZwvXq3NkKWHOXojIl15DY3McHJ04Ote/jyX342BSw7n2Mjr3D1Fc/jL1/zwLezJb+pW2B92pd8Kq9MbGLkmmA87VKZ8Ltf4zfMO9dy9aBbuVQxj0ks1eKdNBb7aepEvNp1PVzYxRc+by49S9bl8tKlSlGrFvdk2shVg3ak2nr4dgP2jFIk/RwfBVdWptpi1k6tmEQbbzz3gUkgsiwZYZg1oV60Yc3r5M3zZUQ5dDc/ZD6+RJ3B0EMx6tTYVfLwY8HOgVbGg7GAta0DJAu6sfrMJfxy9zTfbLlr8Fn7YdQVvNyf6N/ZFb5CUHb2eiBzUdX2WyNMO1WgYYzpXpVRBxTA+6ViVj5+vzI+7rjDpnzMW5a0pQ1Us6sXuT1oDilM1xtr9dfwO4XHJrBrW2EKd3MEhNdtp2692c+F+DIkpeob8coQKPp60r14sXT1bVfZhXp86jPj1GPsuhaU7r6GRGQ4Ogumv1KR2qQLpxIKyyzfbLlGqYPqsAcW83fhteBM2nrrP7M0XkFJy6nYUvxy4zuwefgghTHHZ9rJ2k13y9KLUF5vOcyUklvn966WbA/153zU+DzhLz/qlmNXDjxWHbjD+r9M2V97vRCbQbKaSKSD4s+fxm7yFBmULsvatplbfW0pJuTEbLI5dntaZyIQUM20BJVrA+Hr9yXumsmlXSe1t8l5D4UkuShkxrsQbG/IZG8+z91IYK4Y0tBnCZ4sj1yN4+1dFV6KIjWsj4pJ5fVEgdcoU4NDVcN5vV4kXahVXIgCS1QiAJJ0SBaBGA8Ql6VkZdIMvXq1toepmb3adZxeljlyP4Pejt9n4QQurC0pDmpfD3dmRsX+e4m5kIvsuh9GhejGTM5XSMilfSEwiw1uWZ8Geqyat1YfxKfhN3pJOMMIWFdPoAWho5ATd65bielicyc7HdK6Kq9NF+iw8xIqhjSiaL/0CqTXiknR8tFbRlbDlTAEKebqwalhj/KYov4NP1p5k5JpgPFwc0638e7o44enqxIX7MVx4EEOK3j46eLbIkz3UU7ejePG7fdQs6Y1fqQI8yCApX05gSsqn5uoxhku5OjlazNWuGNKI5pVsB/J/teUCZ+9Gs+iN+ukaAXtryTUUcrqHuudiKB4ujiSmGJT4z2Sd2htUeoVKTKie2GSdacSzeGB92lZNP9WUlrF/niJZZ+DL17KmDHX85kP6LgykY41ifNXTH0cbCTPjk3V0mbuXUZ2q0rlWcYtz9mbXebKH+uJ3ym6o03eiOX0nOt35gh7Oplw9e9PMWa4e3thqUj6A0JgkGkzbRvtqxdh2Tsn9c3ZKR5uhKkZJwGszutB+zm5e/zmQXwY3pFVln3Rlj918yKqgW2z4oPlTl2jTeHYx7t5rWK4Q+d2dU2NC1bjQkgXcTa+NDnXw0iOsHt6YRuUL27xvRspQCcl6uv94gA3vW9pmnTIFOTy+PYOXHuaT34OZ3cPPqlOdvuEcdcsUTOdM7ZKnrc6S1UdOqfIYWX34pmw4bau8cD/aZpnbD+Ol76gAOemf03Lb2fvSd1SA7DJ3j83yRpUpKaWMSUwxvY5KSE5X9qvN56XvqAB5JSTGdOzFeXul76gAufn0PYuycUkpsvXsnXLDybs23xs7U+XRHk/GrhOSdXLAz4HyreVHZFKK3ma5//12XPqOCpBhMYmy9eyd0ndUgNxzMcRq2cyUoY7eiLCprialYr99Fx6U7648JpN1lnXadSFENp2xXUbGp/+NSGl/dp1nV/l71i/N2C7V6LcokLN30/dSAdMi02cv1qBdtWKsGNKIM3ejaTV7Z7qyxoiAA2oWVC9XJ85N6QRA7UlbLFZV70Ym8O2OywxsWtYiLu+fd5vTuHwhhi8/yr/Bd03HZ2w4T53SBXJHC67xRHFzdmTBgHqk6A2M+PUoSTp9ujLn7kXz5/E7jH+hGoW9XNn5cWtqlcxP/5+DTFlVjUiZuTJU3TIFuT7zBZvbWz1cnPj5jQZEJ6Tw3srjJOsUfYHI+GRG/X6SWT1qk989dwgA5VmHCtDNvySTX6rBgMWBnLwdaXFu5JoTABwel6oM1bxSEda82YQb4fHUU3PmAFwLi2Ppgeu837YiJQq4m467uzhyYariVP2nbCUsVskE0FR11NZk+X4b3oR2VYvy3qrjrD1yi90XQ9l+7gGfaRJ+GlnE1cmRH/rVw9nRgeHLjlpISEop6Tx3LwBDW5Q3Hf/3veY0KleIYcuOEHAytTHPKWUoo6PXGaTJ0Y//6zSdaz1Hs1wkAJSnHSooKlIzutdm0JLDHL2hZBa9cD+GdcfuMLZL1XQq+Q3LFeKvd5oRHpdMZXVVvo2q/j/y+fRG5+rkyCU122n9qdvo/7MirBI0rp1FuSSdnjuRCRy/+ZCeDUoD8MnvJ3ljcRCzX/PLNS24xn+DMbNoAQ9nBi9NzSzac/5BAJMSmjmr31Qa83dXHuf3o7e5F5XAlH/PPpYylJSSxBQ9EXHJhEQn8b/2lTh9J5oq4zex4dQ9RnWq+ugf8hkkT67yW2PXhRA+WhPM9/3q0nvBIcC2Io7BINl3OcxCcPrNluVJ0hkIiUkkJDrJpDaVkJJ+yJVdsqLMY2+roRoKT9qu9QbJqD9OcjM8nuEtyzN02RG+7VOHl/xKWC2fmKKn78JDHLsZCUD7akUZ3LwccUlKDKkxSiAuTfSA8bzxWLwxoiBZj4ODwNPFEQ8XJ7xcnYhMSOZBdBL+pQvw1zvNMqy/vdm15lDNSCvNlxPkc3XCx9uVIl6uBF1L3Y//ftuKNChXyKQ6ZZ6UD5Sh1txtF1n/foss9Q7szfA0FP4LuzYYJB+vDWbd8TuPdH3j8oXUuFFjtIAjHq5OeLmmOkkPi2gC9byLE+4ujrg4pQ6EpZS8seQw9coU5IP2lTJ4VwV7s+s8GTZlC/O5HCcHwWv1S1sk4zM+L+zpipODoPxYy51OV6Z3sRlrZxxqnZr0PLUmbeHbHZfZ+mFLKhXLl67s/ahEJv9zhiWDGuRqZR6N/wYHB8H07rVMDrV5xSJUKuaVLtlk0XyuFPBwZueFEAYvPUKpgu7cfphA0wpFLMSjH4cVgTeJik9mRJsKOXK/Z46nHWaQ1UdOh5dkxJHrEbLulC1yy5n7Nst0/Hq39B0VIBOSdabwKt9RATJFlz5U5eCVMOk7KkD+eey2lFJKg8FgKn/6TqRFWYPBIPv/HCi/3nohW3XGzsJLtMd/b9cGg0FOX39Wdvx6twyLSbRaJjFFJ31HBcg2X+6UUko5a9M56TsqQM7ceO6x3/9qaKysM2WLvGwWKpgZ9mbXeX5Ryhr1fAuyZFADxqw7xYZT99Kd33k+hPP3Y1ioKkOVLOBO4FhlkaniuI2msBBQpMt6LzhEQQ9nXq5TElDSWl+f+QIuTg688O0+TtyKNJVfEXiTyPhk3mlT8cl+SI08hxCC0Z2r8nz1YvRecIiQmMR0ZYzbprd+qKiofdKxKh91sC4WlB0UCb8TvN+2IhVysYSf5lBtULtUAZYNbshn/5zhL7O5pySdnkFLD1O+iCcdzJShinm7cURNvld5/EZTqIoxy2ng2PSJ+S5O7UxhTxde/n4/QdciuBYWx5wtF5jT0x9nR+1fo5HzCCEY+XwVXvIrQe/5h7gXlWA6t+bwLRJTDPzxdlOLqav32lVi/AvVWHrgOqN+P/lI7/vT7it4ujgxoEnZx/0IzzTarzYDqpfw5tehjZix8RxrjtwCoM4UJf50q6qDak4RL1dOTFRSSFedsIm1R25xLyqRFUMaWUzMm3NkfHsKebrQc/5B2ny5i/fbVaJi0dzbgms8G7zXrhK9G5am1/xD3H4YT1RCCp/+cZLWVXyo51swXfmhLcoz9eWarD5yi3dXHrNyR9ucvhPFkv3XmdWjtsXCa25EW+XPAldCY3l9USAOQnAnMoF5fepQqqC7RVI+4/OQmESuhsYRn/xo4VJXp3d5JKOzt9VQDYWnadegSPst2nuNO5FKT/XaDNvpfAB+P3qbj9cG07ZqURYPbJDp/RNT9Lw4bx/vtKlomvLKDvZm19oqfxao4OPFssEN6fC1kqLkvVXHs3W9f+kCVCrqZZaULzVBn08+V66GxvHG4iD+eLtprm/BNZ4tBjUrR8DJe9yJTGDbyFaZCu/0qFcKN2cH3l15nF7zD7L6zSYZlv9qywUqFfOim7/1uNfchuZQs8jao7dNzwt7unB0QgebZRfvu8aUgLOsG9GU7j8c4MStSH5+o75VMd/EFD0j15xgfNdqlC3i+UTqrqFhi5jEFB5EJzLtlZpZnmrqWrsErk6ODFt2hBfn7ePf95pbLXfoajj/BN9l4wct84xCmjaHmgUCr4bz1/E7HB3fnj9HNLXYdpqWiLhkpgSc5UW/EtQtU5CLU5Vtp/WmbiMkOv2q6pytF6ng48XL/tkfDmloPC5TA87RolIR+jXyzdZ1HaoXY9nghpy6E0Vbdeu1OTGJKXy0JpgZ3WtRyNMlh2r77KM51EyITdLx8e/BTH+lFoW9XKlTpiDr329Ost5A2dHrSTsHXVcVTZnXpw6g7Km+rO7lbzh9u2muCpQW/O8Td5j2Sq0804JrPDtsO/uAA1fDGPdC9Ue6vmVlH1YPb8zVsDgaTNtmce7zgLO0rFwkS8LVuQnNoWbC1ICzNC1fxCJ5Xo0S+dn6YUvAMtvp9A3nANj7aRuLezg5OpiynTabuYMb4XHEJKbw8VrFUeelFlzj2SA8Nomxf57iq9f88XJ99Jm/RuULs25EU0Jjkqg2QUliueXMfQ5djWD8Izpqe0ZzqBmw7ewD9l8JY8KL6Q2jUrF87Py4NaA41RvhcSzYc5U3W5U3pdY1xzzbaavZuxi05DDNKxahXbW81YJrPH2klIz78zSv1ClJw3KFHvt+dcsUJOC95iSk6Ck7ej3j/jrNVz390mW0yAvkvU+cRYwt+Lw+dfBydUJKSbLeYMrYGJ+sJzZJx8zutRi97hStZu8CIL+7M19uvpBOjceo0FO2sAfXw+M5cuMhSwc3fLofUiNP8ufxO1wLi2NuH39AcbDxyanZRxX1KJ3ZMfNMpakZS81/B/FJqWGCoTFJNCj7+I7aHtEcqg3qTVXmhD5cfcJkRA5C4OFqzNToaMrTU8HHkyuhcfiVLkBUQgqeLk4U8HBXVHjMVHnMr7scEvtYQy0NjUfhSmgsI9cEA9Dii53EJelISNHj6mRUi0pv36lKUo54uztRPL8bHi6O6ezby9WJ2w8TqPpcesGfvIL2i7bBkkENKOLpSkFPZ5N0ma3dTo9CMe+spe7V0MhJfAt58NVrftQpU8Akt+fu7GhTJS27WJvuyktoDtUGbaoUfdpV0NDIcZwcHXi1XqmnXY1ci7YopaGhoZFDaA5VQ0NDI4ewG3EUIcTTU5CwD8KklJ2ediU0sodm15liV3ZtNw5VQ0ND41lHG/JraGho5BCaQ9XQ0NDIITSHqqGhoZFD5AqHKoRoLYS4nXnJnL32SSOEqC6EOCKeYSkqIcQfQojOT7seuQnNnp8dhBCuQojzQgifrJTP1KEKIa4LIRKEELFmj+8ev6r2gRCimxDihBAiWggRJoTYIYQop56bJIRYYVZWCiHi1O8oXAixXQjR6zHe/nPgS/lsrxx+AUx92pXIKpo9a/acEUKIXUKIocbXUsokYDEwOivXZ3Wn1ItSym2ZF8tdCCEqAsuA7sAOwAt4HsgoYZSflPKyEKII0Bn4TghRVUo5OZvvXRxoA/R7pMqn3sdJSql7nHtkhJQySAjhLYSoL6W0lxAgzZ7t1J6t3DedfT8Bm18JnBBCjFUdrE0ea8gvhPhRCPGH2esv1FZMGIceQoixakt4XQjRz6zsC0KI42pLeUsIMcnsXFm1dXxDCHFTvX6c2Xl3IcRSIcRDIcRZwCJbmBCihDoUDRVCXBNCvJ/Va9PgD1yTUm6XCjFSyj+klDcz+26klGFSyuXA28AYIURh9f0HCiGuCiFi1LrZMrAOwDEpZaJ63SAhxL9mn+OSEGKt2etbQgh/9bkUQrwjhLgEXFKPDRNCXBZCRAgh/hFClDC7Vgoh3lLvGSmE+N44LBNCOAohvlL/B9eEEO+q5c0b413AC5l9J886mj3bJqftWb22tBBinfq5woU6UhBCOAghxgshbgghQoQQy4QQ+dN8l0OEEDeBHWod9gshvhZChAOThDJU/1L9vh8IIX4SQribvbd5T/2KEKKTEGIa0AKl0TCNXKSUt4GHQOPMvieklBk+gOtAexvnPICLwEC1ImFAKfVca0AHzAFcgVZAHFDF7HwtFKdeG3gAvKyeKwtIYCHgDvgBSUA19fxMYC9QCCgNnAZuq+ccgKPARMAFKA9cBTpmdq2Vz1ceSAS+RmldvdKcnwSsMHstgYppyjir30NnwBOINvsOigM1bLz3bOD7NHWJVD9fCeCG2Wcuj/IPdzCrx1b1M7oDbdX/TV31fzEP2JOm3gFAAaAMEAp0Us+9BZwFSgEFgW1qeSez60cC6zKzpWfhodnzM2PPjkCwWhdPwA1orp4bDFxW6+sFrAOWp/kul6nXuav/Lx3wHsqo21297z/q95IP+BeYod6jIRCF4uQdgJJAVfXcLmColfr/A7yfqX1l0QBjUX7Mxscws/ONgAiUH3gfs+Ot1Q/paXZsDTDBxvt8A3yd5ksrZXY+COitPr+K+oNXXw8n1QAbATfT3HsMsCSza23Uq7Fa71AUY1yKaohkwQDV4/dRhjqe6vf3KuCeyfe+EJiZ5tgtFKfYG1igfidVgUHAP2nq0dbs9c/ALLPXXkAKUNasfPM0/6fR6vMdwJtm59qT3qEOA3ZkxaE97Ydmz8+GPQNN1Do4WSm7HRhh9rqKaq9OZt9lebPzA82/I0CgNHYV0rzfNfX5fOP/xsp778K6Q/0VmJiZfWV1yP+ylLKA2WOh8YSUMlD9pwr1H2XOQyllnNnrGyi9K4QQjYQQO9XufhRKT6hImuvvmz2PR3EEqPe4lea+RnyBEurQNVIIEQmMBYpl4dp0SCkPSSl7Sil9UHotLYFxGV1jjhDCGfABItTvohfKZ70nhFgvhKhq49KHKC2rObtRftgt1ee7UHpKrdTX5ph/RmOP1viZYoFwlJbZSFa/a/PnRvKh/LDsBc2en749lwZuSOtznRb2qj53IvUzQ3o7NH/tgzLaOGr2nW1Sjxvf+4qNetoiSzb+2GFTQoh3UIZAd4FP05wuKIQwz41cRi0HykTvP0BpKWV+4CcUI84K91C+FPP7GrmF0hKZ/2DySSm7ZOHaDJFSHkYZftTM6jVAN5SeTZB6j81Syg4ow6PzKC23NU4CldMcMzrUFurz3dh2qNLs+V2UHyYA6v+kMHAnC/W/hzLcN1LaSplqKMM3u0ez50zJKXu+BZRJMxdvxMJeUT6TDmUaxVT9NNeYvw4DElCmH4zfWX4ppbEBuwVUsFHPtPc1kiUbf9xFqcooITOvA/2BT4W6MGLGZCGEixCiBdAVMC6k5ENp5RKFEA2Bvtl46zUoE+MFhRClUOZOjAQBMUKIUeqEvaMQoqYQokEWrk37+ZoLZTGnqPq6KvAScCizCgohCqkT9N8DX0gpw4UQxdTJcE+UObRYwGDjFluBukIIcyXq3ShzX+5SmSjfC3RCcY7HM6jOKmCQEMJfCOEKTAcCpZTXM/scKN/XB0KIkkKIAsAoK2VaAdbzatsRmj3b5gnYcxBKYzBTCOEphHATQjRTz60CPhRClBNCeKHY62obvdl0SCkNKI79a7PPWlII0VEt8jPK76GdUBbASpr1rB+gzN2af/aSKHOxmX5PWZ1zSkD5soyPP1G64EGoc21q2beBUygtfGvgNspwIgy4CfQ3K9sDpSsfg7Ig8h3q/A2p8yTm83S7UOc2ULrzy1C64GeBTzCbN0IZMqxCGWI9VL+I9lm5Ns1nr4kymf1A/dzXUeIunaXtOac4tWwEsBPoa3a+OIpTjFLffxdQPYPvfi3QK82xe6jzZ+rrI8DGNGWsLSa8hTLMiVC/71K2yqPMq01VnzuhTPCHA9eAD1Hms4zCOg1QVm+fyJxnTj80e3527Bml5/mXalthwLfqcQeURbhbKPOsK4CCGXyXA4F9ad7LDcURX0VZODuH2aIS8ApKrzkGZQHMuMjXBGVh8qFZfT4B5mTFvp6Y2pQQojXKP0eTB39EhBDVgV+AhvJJ/aOyiVB2Rf0kpfRVX/8B/Cyl3PB0a/Zk0ez58XkW7Tkz1BFdMNBSShmSWXktBcozjJQys7jCJ44au9cG2IKyKPAZSo8OACnlq0+pahp2xrNgz9lFKoH8thba0pEr9vJrPFEEMBllCHQcZeg08anWSEPjGUUTmNbQ0NDIIbQeqoaGhkYOYTdzqJ06dZJhYWFPuxrPLEePHt0s7Sj3joaCZtcZY292bTcONSwsjCNH7EXM6L9HKGpAGnaGZtcZY292rQ35NTQ0NHKIPO1QpZQcvRGBtjCnkdtITNFz8nbk065GnsNuhvxZpezo9ZQs4E6nms8Rn6wjNklPfJKO2CQd8cl64pJ0xCXriEvSE5uk7GQb2LQsk16q8ZRrrqFhmzuRCTSbuYOhzcuhM0jVni3tOy5ZR7xq10k6ZQfooTHteC6/WyZ318gpcpVDvRwSCyjGdykklo41iuHl6oSHixOero6m516uTrg6OVDn860AbDx9j+olvOlZ35ruh4bG0+erzRcAWLTvGmO7VKWAuwuerk54mOxa+evp6sSO8yF8+vtJWlfxoef8g6wc1ohSBT2e8ifIG+Qah5qiNzByzQmmvlyTlYE32XMxlP6NfelQvZjV8o2mKxkwLk3rzM2IeF5fFEiK3kC/Rr5Wy2toPC1O3Ipkz6VQNrzfgi7f7mX6hvOc/7wTbs6O6crGJ+v49PeT1PMtyNJBDVmy/xq95h9i5bBG+Bb2tHJ3jZwk18yhfr/zMgU9XOjXqAwbPmhBo3KFGLbsCP8G301X9t/guzyITmLl0EY4OzpQwceL34Y35oedV1i879pTqL2GhnUSkvWMXHOCSS/VoHoJb4InPg9A1QmbiE9OL75UfeJmAH5/qwkAg5qVY0SbCvRecIgrobH/XcXzKLnCoZ68HcmKQzeY1aM2Qs1Qu/rNJrSrWpT3Vh1nzZFU7dn4ZB3vrTpO3TIFaFoxNSLDt7Anq99szC8Hr/Pjruxqz2poPBm+2HSemiXy07W2kgIsv4czpycrKnTVJ24mOjHFVHbR3qsAbPyghel3ANCvkS8fPV+FvgsPceF+zH9Y+7yH3TvUxBQ9H65WWvBi3paT7z8PbEA3/xJ8+vtJlh28DqS24H+83TTdvUoV9GD18CasPXqLb7Zd1Fb/NZ4q+y6FsfnMfT7vZqn/7OXqxLkpSqx77UlbeBiXTFhsElPXn6N7nZJUK+6d7l496pVibJdq9FsUyOk7Uf9J/fMidrOXv379+tJaAPTkf88QHpvMt33q2Lx2zLpTrAq6SSFPFyLiktn4QQsq+HgRFpvEg+hEQmKSCIlJIlR9fvZeNCdvK0Z3eVpnnByf/XZHCHFUSln/addDI3vYsuuohBQ6f7OHma/WpmVlHytXQpJOT5XxmyyOXZvRhdgknWLT0UmExCQSqtp3SHQif51InQK7PvPZT1Rrb3Zt1w61y9y9nL0XTfXi3iTp9ITEJBGTmLMp6Ac1K8vErtUthlDPIvZmeBoK1uw6RW+g0jglAUKdMgUIiU4iNCaJZL0tMfxH4/e3mlC/bKEcvWdOY292bder/GfvRVv8BXB1cqCotytF87lRNJ+r8vB2Y7YadgLQvW5Jvuzhh4ODdSf5464rfLHpPGvfasLUgLNM+Ps0U16qabO8hkZOEhKTZHp+/Gak6XkBD2fVphXb9vF2JSFZz7KDqfnsNn7QwuqQ30jZ0esBWDa4IcOXH+W7vnVoWsGudnc+2zztlBRZfdSrV0+mxWAwyC82npMd5uySIdGJ6c4b+fzfM9J3VIC8FREnf9p1WfqOCpBj1p20WvZBVIL0HRUgP1l7QkopZXRCsnz1h/3y07XBUqc32HyPpw1wRD4D/yft8fh2LaWUoTGJsuPXu+WMDeekwWDb7nxHBUjfUQFSrzeYnl96EGO17Jh1J6XvqAB5LzJBSinlgcthsu6ULXL3hRCb93/a2JtdP/uTgxkghOCTjlV4oVYJei04yP2oxHRlboTHsWjfNUa0rkCpgh682aoCn3erwcrAm3zwW/q8dg2nbwdgVg8/APK5OfPL4IbciIjjk7XB6HJ42KWhYY0iXq6sGtaYfZdDmRJwFsW3WDJsmTJVcGJiBxwcBNdmKIlQ28/ZzTmzURvA5ZAYVgbe5KMOlU07p5pUKMz8/vX4cPUJtp97gMbjY9cOFRSn+kH7SvSsX5peCw5y+2G8xflWs3cB8Gmn1CwG/ZuUZVaP2vx94i5DfzlsOj7+r1OAsl3PHE9XJ5YMbEhobBIfrD5BiuZUNf4DCnq68OvQxhy/Gcn4v05jMKQ61eBbkWw9+4CpL9ekgIcLoPwWrs3ogrOjoPPcvZy4FQkoo9D2c/YA8F67ShbvUb9sIX4e2IBRf5xk0+l7/80Hy8XYvUM18larCgxsWpZe8w9xIzwOgDeXKy348Qkd0pXvWb808/rUYdu5EHovOMjlkFhWHLrJh+0rW9377O7iyMIB9UlM1jPi12MkpuhN5wzq3uqQ6ESuhcVx+k4UgVfD2XH+Af8G36Xs6PW8veLoE/rkGrmZ/O7OrBjaiEsPYvn0j5PoDRKDQdLt+/0IAa83ttzZJ4Tg0rQuFPRw5uXv9xN0LYIBi4MAODnpeavv4V+6AEsHNWTC32f4+8Qdi3MpegNR8SnKdu4HMZy4Fcn+y2FsOXOfP4/fpuzo9VZHhnkVu17lt8avgTf4bsdlRnWqyv9Wn+DzbjXo36SszfJbzz4wDZ0AmlUsrIabJBGVkGLzOgAXJwecHAQJKXrcnR1VnQBHk16Ah6sjD+NTCFZ7Cpv+14Kqz9leMHgc7G01VEMhq3Ydn6xj6C9HKOLlytl70VwOibW5/dRIy1k7uRmROmKrW6aAKUQwWZfxKKughzNxSXr0UuLp4oinqhNgfO7h4sQ2dZqgRglv/hzRDBennO+f2Ztd2/UqvzX6NfLF2dGB/60+AZChMwXoUL0Y/Rv7svyQslK6/3K46ZyLowM++VzVqAFldfVhfDIBJ5WhUbLOwPZP21CygLvNCADjquq8PnV4fVEQSwc1oGbJ/I/5KTXyGh4uTiwe2IBqEzchJSzoXy9DZwqw59M2JvsDOGYWMeDt5kRRbzd8vFLte+He1G3Xnq5OHBzTDlcnB6shg6fvRLHt3APGdK7KsZsPeWvFUX7oVzfTOuV2cp1DBUyt7ykbQ5y0fP5yTbrWLk6vBYco7OnCUStTBKDMRZUbo6Sfvzi1M5XHb6TFrJ0EjW1HUe/00wTvr1IWvY6Ob09hL1ecHR0YuCSIhQPqU6dMwUf5aBp5GCcHQQUfLx7GJfN8jeeydM31mS/Q86eDBF2PYF6fOrzoV8JquYNXwlm49xpze/tT0MOFAYuD6DJ3Lzs+bp2urJSSrvP2AfBmqwqk6A18uPoEw5YdYUH/+ri75F2nmmvmUI1cD4vjqy0X2DayJfncnLN8XaPyhflzRFPC45KprAZVp6X3gkMAnJ7cERcnBy5P6wwokQF3IhMsyp6+E8U/wXeZ2LU6hb1cAehU8zlm9ajN0F+OcPh6xKN8PI08zPw9Vynm7crhce2zdd2at5rQVtW1WGuma2FEpzfQZ6HSmejmX5KWlX1YPbwxV8PiqD91W7ryL/9wAMC0/dXZ0YG5vevg4+XKwCVBxCXl7OYaeyJXOVS9QTJyzQnea1uJikXzZfv6OmUKsv795iTrDZQdvd4iVOXQ1XACr0XwTS9/vFyVjr2TowNXpyuhKs1m7uB6mLIYZt6CD25ezuI92lYtxte9/Hlz+VEOXNGSs2lkjdN3oli87xqzM9iQkhGLBzbgRb8SfPL7SdP0lpFmX+wA4NDY1OiWRuULs25EU8Jik6g8PrWDsfdSKMG3IvmhX12Lnqijg2D2a36UK+LJgMVBFqIteYlc5VB/2n0FN2dHBjYt+8j3qFEiP9tGtgSg3JgNSCnR6Q30XnBIWTmtU9KivHn8X+svd3HpQYypBT87paPV92hZ2Yfv+9bl3ZXH2X0x9JHrqpE3SExRJPzGvVCNEgXcH/k+8/rUoU/D0kz46zQL9yjKVP+oUpa/qlKW5tQtU5CA95qTrFM6GMk6A/1/DqJkAXe61Cqe7v6ODoLpr9SiRglv+i8KJCo+7znVXLPKf+ZuFP1/DuLf95pT8jGMzsiN8DhTDKtPPldCY5K4OLVzupXMxBS9SYSix08HTcfzuTpRr2xBU8RAWGwSGfG4QhX2thqqoZCVVf7pG85xKyKeH/rVzRFNiSn/nmXx/mu82ao883dfpU6ZAvw5oplFGSklkfEphMYmceByGJP+PWs690Kt4oTFJhEam0RodBIxGQzxg8a1o2i+R0/BYm92nSscapJOz0vz9jO8ZXlerVfK5j0MBklkQgohMYkmR2d8Hmp8rqr0JJjFmeYURbyU1VQfVWNg7dHbpnM/vV6XTjXTt/pZxd4MT0MhM4caeDWc91YdZ9P/WlLI08VmuWSdwaZ6mrmdm+sE5BSuTg4mmy6az40CHs78dliZq61SLB8rhjbCJ5/rI93b3uw6V6zyGyXMPlobzEdrg3PknvlcnfByc+KeWdDyRx0qU6qQOz5ebqZQk/zuzo/Ua9h5PoS1R2+zoH89ShRwZ9DSwyTrJS/ZWIXVyHv8FnST0euU3Xt11fxnOYFxxGXO5JdqqEJCilP0yef6yCFQVdQ51yvTu/Ddjsv0XnCQlcMap9Mrzo3YvUMNuma5Wt62alGTypSPd6riVDFvN4p4uWYr+NgYw3d0fHvqTd3GV1svZhpMnRWSdHoGLT1MuSKepvCXFUMa0f/nQFJ0hgx72Rp5h2nrz5me1ylTgBL53c3iot0sHGABd+csL1aFxyZRb+o2uvmXoE7pAkz69yzBtyOZ09P/seu8+vBNknQG1o1oiqODsi3cxcmBXvMVp/o4c8D2gF071NgkHR+tPcHCAfWZuv4sN8LjGdSsLC0qWRfkzQ7G3FLr329OYS9Xgj97Hr/JW6g6YRNnJnfE0/XRvzr/yUpvY9vIVqZjVZ7Lx8phjRWnqjfQu2GZx/sAGnbN9nMPyO/hzF/vNqPdV7s5fjOS1VOb5MhupHpqKNTc3ooou5uzI6PXnSIuScf8/o8+uo5KSGHUH6doU8WHumZx1m+3rqA41QUHWTm0MaUL5d4MrHa9yj814CxNyhemQ/Vi7P6kDTVLetP/5yC2nn085Zzw2CSmBJzlRb8S1Cih7GrK7+7MGTWXT43PNme6LdUWqw/fJCFFzx9vKy24ORWLerFqWGPm7bhsStmikfcIj01izLpTfPWaHxV8vDg6Xok7rTx+o4WGxKMwNUBZXNo3qo3pWO+GZZjb25/NZx7w+qLAR7633+QtgBKilZYhzcsxvKWSLPCaGl6YG7Fbh7r93AP2XQ5jQtfqpmMB72Wc7TSrGFvweWnSqni6OnH+cyWY2W/yFiLikrN1X2ML3rqKD/V8re+UKlvEk9+GN2bh3qumpGsaeQcpJeP+PM3LdUrSqHxhAGWEZJbt9FED52+Gx1tIWZrTzb8kP71ej32Xw+j2/f5s3/vrrRcB2PFRK5trCv0b+/J+u4r0WXCIyyG5M1mgXTrUiLhkxqw7xZev+aXbDWUr22lWmbZeacH3ftrG6nk3Z0cuTFWcat3PtxISk3WlHWMLvsRKC25O6UJKssBfA2/y/c7LWb6/hv3z14k7XA2LZWSHyhbHzbOd1vhs8yMFzrecvROwlLI0p1PN51gyqAHBtyLpMGd3lu97LyqBudsvMbBpWcr7eGVYtleDMozqXIW+CwM5fz86w7L2iN05VKUFP0U3/xI0VlvwtJhnO12wJ+spoW+Gx7Nw7zXealUhw3keVydHLhm3nU7bzt00206tkZUW3JwSBdxZPbwx647d5uutWgbWvMDdyASmBpxjTk9/qwufabOd3ovK3O6MGOUjrUlZmtOmSlFWDm3EpZBYmszYnqV7N5mh7LSa9FKNLJV/pU4pJnStzuuLgnJdBla7i0M1rryXKeRBZHwy0VlIyteiUhGWD2mUaTnjvbMaZK83SCqMVcRS9nzShjKFrTvhe1EJNJmxgwFNfJmSJiVwZoTFJvH6okDaVC3Kpx2r2HTG9havp6FgtGvjyjtAqYLuWZLYg6wFzp+6HcWL3+3LVMrSnKM3Inj1x4N4ujhyRnXi1vhoTTB/HLvN4XHtsx1ruun0fcb/dYpFbzTAv3QBq2Xsza7trodq5GZEPNGJOlycHChV0J26ZQrQqcZzDGjiy8fPV2bWq7VNZfdeCjP1EG2R1RbcHEcHYdrL33L2Ti6HxFotZ2zBs+tMITUVxp6LoUxdf07rqeZSztxNHf7efphAss6At5sTlYp60axiYV6pU5I3W5Zn/AvV6GMWAZLZCMlgkLz4naIrkVVnClDPtxD/vtucuGR9Ol0LIxfux/DHsduM7VL1kQL3jWJBQ5Ye5khuEQt62kmtsvowJjMzGAxyyr9nZOdv9sjw2KR0Sb2MvPbjAek7KkDGJqbIr7ZckL6jAuS09Wetlj11O1L6jgqQvxy4ZvN+GWEwpCZIO3s3yuLcyNUnpO+oAPkgOuGR7m0kMi5ZvvTdPjn+z1NSbyVZIHaWzEx7WNq1lEqCyHZf7ZJfbT5vMzFfsk4vfUcFyHqfb5U6s8R8N8LirJbvMneP9B0VIOOTdFbPZ8a5e1Gm9zCvk7nNPy67L4TIulO2yAOXw9Kdsze7trseqhCC8S9Uo3UVH3ovOGh1UejQ1XCCrivKUJ6uTozsUJnRnauyYM9Vxv15yqKswZCqDDUgGy142jpdm9EFRwcll49Rod/Ygo/pXPWx9jODsiixYkhDzt2LZuyfp9AbtJ5qbqOotxu/DW/MlrMPmLnpPIo/saTZTGW0c3BM20xHSLsvhnLmbjQ/vV73kTVKqz7nzY6PlHjpcmM2mPJa9ZqvSFkaQwkfh5aVffiub13eXXmMPXYuFmR3c6hGpJR8u/0yf5+4w8phjU15oHR6AxXHbaSQpwvH0gzffzlwnc/+OcMrdUrydS9/ALrO28vpO9Gcm9IpR4Rx/SZvISohhTVvNqHnfEUsJStzsgnJekJjkgiNTd1zba4vEBqTeszI1eldTLtj7G2uSUPB2l7+h3HJ9F8cSH3fQnz2YnXTvPm/wXd5b9Vxfh3aiGYVi5jKS5kqfL7xgxZUK+5Nss5A5fEbKV3Inb2ftn3set6KiKfFLCVK4Nehjei3KJBv+9TJdKu0lJKohBSTRkZG9h1jth5i/M3Ym13brUM18uOuK6wKusnKYY0oVdCDJjO2cy8qkUvTOqeTIwMlsH7UH6dMqU8GLA56bGGStDSbuSOd4PSTwNxR25vhaSjYsuuohBTeWBxE9RLeTO1Wk0SdnuoTN+NfugB/vdMsXXkpJZXGbURnkPz1TjP6/xxITKKOK9O7pNtA8qgYF1efNPV9C/L7200B+7Nru3eooGwT/XnfNfo1LsOsTRdYMaQRzSultuBxSTq1hVRaw0V7rxJ8WwnXKFXQnX2jHr8FT4t5Lp+KRb0scvf45Evdi2187u3ulKVwqr9P3GHejssEvNfcIrTG3gxPQyEju45N0jF4yWHKFPbgd1WZ7NqMLiY7MRgkEfHJJvnIkJgkPv39pOn6tW81oUHZQjla36BrEaaRV0EPZ0oV9DBTmlLsWXmk2ndWtC8MBskbS4JoULYQ75ulurY3u84VDhVgyf5rTDbTbMwq1jROc4oePx7gyI2H/NivLp2tCPJml3tRCXT9dh9LBzWkVinLRH/2ZngaCpnZdXyyjuoTNz/SvR9XY9cWD+OSqaOqX+XUVNkvB67z5/E7/P5WE5zMRpb2Zte5xqFCaq/wnTYVaF7Rh6LeisqUVxohk3XHbjNyTfATacHTMmBxEHsuhvJ1Lz9eqfPoKlIGg2TA4iAalrNswY3Ym+FpKGTFrndeCGHQksMAfN+3LiULulM0n6tV9TTjb8C8J/skiE5MofYkZeff6ckd0/3GssOV0Fhe++kgv7/VJN1OK3uza7tb5c+I6zNfYG5vf9YcuU0+Nycq+Hil+0fHJKYwck0wTSsUfuLOFGDZ4IZ0qfUcH64OZlXQzUe+z/JDN4hJ0jGidYUcrJ2GPdCmSlEuTO1Eh+rFWHfsNlWfy0eJAu7pnKlxm/LWD1s+UWcK4O2WKhZU87PNj5zuRKc3MHJNMB+2r5TptlV7IFc5VFBEHj7vVoOBS4JM4Uvm1FJb1V+HZr5zKqf4oV89etQrxZh1p0yygNnhSmgs32y7yNc9/SyGQxp5B1cnR1Pe+2HLjpCQbKk6FRKdyOzNF+jdoDSVimU/QeWjYCEWNGUL4Zmk+bHG9zuv4O3mxOuNfXO6ek+FXPnr7FSzOF+8WpvBaXZg/LhL2de/5T9owdPy5Wt+ytbTgLPZEjzR6Q2MXH2CDztUzhUtuMajo6Rr9qeIlyuDllqma244Xdl3P9Nsh+B/gZuzIxenKroW9aZuIyQ662JBJ29HsvzQdWb38PvPf49PilzpUAHaVUtN13zwSjgh0Yl8sek8veqXpvJ/1IKnZUq3mrzVqgKzN1/gqy0XsnTN9zuv4O3uTP9c0oJrPB5Ojg58+ZofvoVS0zUbN6scGtMuk6ufDC5ODlw2igVN356lkMHEFD0frj7BxBdrmGLIcwO51qFC6g6Md1YeM7XgX/T4b1vwtIzuXJUP21dm3o7LJrFfW5y8Hcmyg7mrBdd4fBwdBDO616J6cW+azdjBr4E3+ahD5afqmJwcHbii7tpqNnMHN8IzFpH+YtN5qhX3znU51HK1QwVoUqEwz1cvBliqlD9NPmhfibFdqrJo3zXGrDtltUxqC149V7XgGjmDg4NgSrcaphTO71mJ/PivcXRQtmADtJq9y6aI9IHLYWw8dZ+pL2dfLOhZJ9c71Idxyey8EMLKoY3SqZQ/TYa3rMDn3WqwKugmH/x2PN15Ywvezb/kU6idhj3w94m7VCrqZdJIfRYw6loAtJ+zh7N3LUWkoxJS+OT3k3zRozYFPGynxbZXcrVDlVIy/q/TvFi7BE3N9j4/K/RvUpZZPWrz94m7DP0lNRZxfy5uwTVyhntRCXwecJave/nnSGB9TmJ0qk4Ogi7f7uWEWbTN5H/O0KaqD60qP34izWeRXO1Q/wm+y8UHMXzcscrTropNetYvzbd96rDt3AP6LDiktOBrg5n5aq1c2YJrPD4Gg+STtScZ2LQsNUvmz/yCp4AQgsvTu1DAw5mXv99P0LUINp66x/FbkYztUu1pV++JkWsd6r2oBKb8q7TgWdlL/DTQGyQh0YmUL+JJn4alOXg1HL/JW2hbrSitqxR92tXTeEZZfugGsUk63n6GN3kkJOuVpIADlE1OPecf5O1fj/FVTz88XOw6e32G5MpPJqXk099ztgWXUhJrEllRxCiMknoh0YmExiaZZMkeNcW0EWsqWRoakLrJ44+3m+bYJg+9QRIep8rrpZHVMxdeyWpalowo4O6ceSE7Jlc61J/3XSMiLpmX65Qk+FakagzmOoyJZvqMSTku1uzi6KAo8HgbVXgU5R3ltZvpXCEPF9OP4uzdaLp8u5fBzcqx9dx9SuR3Z1jL8jlaLw37JkVv4L2Vx012cehqeGqDbqaZa3SAkY+4HTQjvN2cKOrtZqEuVTSfG0W9LZ/nc01VT3v5+/2cuBXJ5Jdq0HdhIMuHNPzPdnP91+Q6h3rpQQxT158DMAniZkYhTxcL4zDJkXlbyux5PoYAREZIKeny7V4AJr5YnaEtytFvUSBJOj3vtn364TAazwaVxm0E4Oy9aGZtynxjiJODUGzXzAEaHZ7580KeLk9sVLT/chgnbkXyXd86dK1dAm93J/ouCmTZ4IZUK+79RN7zaZLrHGrFol6Mf6Ea3++8zMIB9an/HwigPC6v/ngAgLNTFLEJYwrpvosCSdZLPmxfSQvs1yB44vO8/etRpIRFb9R/Yg18TpGiN9BvUSDFvF3pWlsJ4H+lTilcHB3p/3MQSwc1eGYX1R6VXDdZJ4RgaIvyzO1dx7Tt9FnmwOUwjt1UWnDzyXpTfqEz9/li0wWr+YU08hb5PZxZPqQRZQp58MbiIGISc35In5M0Uncn7k8j4P5C7eJMf6UmA5cEcfzmw6dRtSdGrnOoRuwh8VeK3kDfNC24OcYU0vsuhzIl4KzmVDVM206rFffm9Z+DHlk270nz1/E7RMQl89vwxlYXz56v8Ryze/gx9JcjHM4tKaTJZQLT1jhyPYI3lx9l9mu1aVu12BOo2aNhnlhtbm9/IuKSzRYVjKutiYTFJltcZ0uF3d6EeDUUHtWupZRMXX+OQ1fDWT6kEYU8n52Y5ethcbT+chcAX73mZ1ooC7Ww7SRizdSy8rk5cWpS+gyq9mbXz/YkTA5Qv2whfh7YgKG/HGbqy7XoVPO5p10lAH5QpQQBPvjtRIZlC3m6EBGXnGEZjbyFMZ367M0X6LvwEMuHNMInn+vTrhaAyZkCfLQ22GY5N2cHElOUMKw3c0lES653qAD+pQuwdFBDBi45TIrewItPWeEmNknH7M2pq7QdaxRjfn/bjfC87ZcIvBbBssEN/4vqadgJQgg+6VgFFycHei84yMphjSnm/XSFdH7afcXi9T/vNqN2qQJWyyam6Ok6bx/vta2YazQrcu0calpqlszPiqEN+TzgLH+oGSSfWl0+U5KuXZvRhbm9/dl85gH9fw60Wvb0nSiWHrjOrB61ccihdMAauQchBP9rX5lX65Wi1/yD3P0P0pfbIiQmkZkbz/NavVJcn/kCXq5OvPTdfo7esD5HOnvzBaoUy5erJPzyjEMFqPqcNyuHNWb25gv89hj5nR6H+WoLvul/LRBC0M2/JD+9Xpe9l8Lo/sN+i7JGCb8JXatTooD706iuhp0wonVFXm/sS68FB7kVEf9U6tBwmrKqP/s1P0BJ3lc8vxuv/niQA5fDLMoeuBJGwMm7TH25Zq4KCcxTDhWUONXfhjdm3o7LLDt4/T9975CYRGZsPE+PeqWo+lxqUHOnmsVZMrABx25G0vHrPabjszdfoFIxL7r5554WXOPJMbRFeYa3KE/vBYe4FpaxwHNOM/Hv0wAcHGMZInVwTDsqFfWi76JAdl4IAZSMqZ+sPcnMV2tT8BlaTMsJ8pxDBShbxJPfhjdm4d6rLNp79T97X2ML/qXagpvTpmpRVg5txIUHMTSbucPUgk97uVauasE1niz9m5Tl/XYV6bPgEJceWBd4zmmuhsay7OAN/te+EsXzpx9JbR3ZCr/SBRi05DCbTt9n8j9naVXFhza5UAAoTyxKWaN0IQ9WD2+ibvE08E6bik/0/YzK/AdGtyUyPtlMS8BcYyAJZ0fBncgE+i4MZMnABrmuBdd48vRqUAYXJwf6LQrklye8xVNKSduvdgPwTpuK3ItKUMKjzGw6JCaRomoEwlsrjgKYUlDnNvKsQwXLLZ5JOsMT2+IZfCuSVeqcbdOZO7J83cnbUbSu4qP1UDWyzSt1SuHs6PDEt3gaY6khVWsgK1wLi8t1204hjztUSN3iWX/qNr7dfolrM7rkqAOTUtLt+9TFphGtK/Bpp6o2y49cc4J1x+6w8YMW/O+3EyTr9Xz8fBXNqWpkm661S+Dk4EDXeft4vXEZpr5cK0fvfyzNttH9o9tS0sbiqflGlp9er8vAJUEseqMB/qUL5GidnjZ5cg41LeYaj+XGbMCQg3J+xtQmwROf5722Fflh1xVmbjxvtez5+9GsO3aHMZ2rUq24N6uGN2bn+VCmbzinbTvVeCTuRSlhVCsO3eT7nZdz7L56g6T7DwdwcbLMdmorwqDHTwcBZajfqWZxvni1NkOWHrYZUmWv5HmHajBIPvk9mE86VmHvp0pW1PJjN6DTP56QLsDxmw/Zfj6EGd1rkd/DmY+er8InHavw0+4rfKauihqRUtLpG0XC781WihJ7IU8XVg5rROC1CCb/q+3l18gel0Ni+Xb7JXZ+3Jq3WlVg9uYLfLk5c9m/rPD818q86cnPnsfRQXBVdaotZu3kamisRdkDl8M4euMh3/apY1LIaletGHN6+TN82VEOXX22BYyyQ553qMsP3SA+Wc+bLctTupAHgWPbAVBx3EZSHsOpGgySV344gLOjoE/DMqbj77SpyMSu1fnl4A0+NtuW13N+agtuTgEPF1YMbcTJ25GM/fN0jvaeNXIvKXoDI9ecYOTzVShXxJPRnasyskNlvtt5mcn/nnmse287+4AroXEsHljflF7IwSyFdNuvdnPhfoypHn0XBVLEyzVdAH+ryj7M61OHEb8eY98lyzhVeyVPO9QrobHM3X6JOT39TYo4xbzdODK+PaBMsiem6B/p3h2/UeJJrQk+DG5ejumv1OL3o7cZ8etRDlwJ4/D1h3zTy9+qxqW3mzPLhjTiSmgsn/x+MsczDGjkPn7YeYUCHi683ii1MX+/XSXGv1CNJfuvM/qPk49038QUPUOXHaFSUa90YkPmKaQ7frOH03eiaDJDWYRNG59qpGnFIvz0ej0++O04O8+HPFKdniVyvdqULXR6A6/+dJAedUvSv0nZdOcj45Pxn7IVgHNTOmUrVe/2cw8Y8ssRFg2oT/vqthWu1h27zcg1qb1U/9IFTGEmKfqM/y9pVafsTZVHQyGn7Rrg5O1IBi89TMB7LXguf/q9/csP3WDCX6d50a8E8/rUyda9K43bQIpecnV6F4ut0HGmfGtKfrV3Vx43nXN0EFQq6kVITFKmIj/NKxZhxdBGptf2Ztd5dpX/h11XyO/uzOuNfa2ez+/uzN5P29Bi1k6qTdzE6uGNaVS+cKb3TUzRM+SXI5T38czQmQJ0r1uKr7dd5FaEsnBgnr88v7tzaioWNQVLXLKelYE3mfxSjax/UI08hXG78sQXa1h1pgB9GpQmOiGF2Zsv8G/wXZuSkGkpO3q96Xn5sRsyKGmJ3iA5r04BCIEp/Yq5fa8+cosH0Ul81zd7Dv5ZI0/2UDvP3cu5e9HZvi544vPk98g4a2Pl8RtJ1hm4Mr0LjlkUMzH2aGuU8Gb9+y2sltHpDbw2/yAv+ZVgULNy6c7bW0uuoZDTPVRzp5dVSuR348CYdhmWiYpPwW/KFtPrWa/WxrewhylhX9qpKmM9yhXx5FpYHMsGN6RlZR+r9z59J4o3Fgfx73vN02lW2Jtd58keaqNyhUwO9bV6pSyS8ZlnJnVzdrQwUL8pWzg6vj2FvazrTv4WdJNknYF1I5pm2ZmCsuK5YkgjXv85kFazd7L7kzbpyszfcxUPF0fesDI9oaGRlvI+njSrUMRqtt3Cnq5ExCXTYNo2Cng4czcqka7z9hLwnvXGHDA501OTnqfWpC18+sdJTk16nnxu6TsYC/co27k3ftCCasW9eem7fQxYHMTCAfXpkGbUZuxRj+9aLXcIAEkp7eJRr149mZPcjYyXrWfvlPO2X7RZZsHuK9J3VIA8dy9KJqXope+oAOk7KkDej0pIVzYyPln6jgqQAxcHPnKdAq+GS99RAbLe51ssjp+6HSnrTtki7zyMt3ktcEQ+A/8n7fF07VpKKdcduyXrT90qz9yJslnGaMtSSrn7Qoj0HRUgW8/eabXsV1suSN9RAfJqaKyUUsr4JJ3p+odxSRZlQ6ITpe+oAPnhb8ctjveaf0D6jgqQAcF3LY5//u8Z+faKI9JgMFh9b3uz6zy7yl88v7Lt9K8Td5mzJX0SvNCYJKZtOEf3OiWp+pw3Lk4OXJ7WGVCSj91+aBnA7DdZacEXD2zwyHVqWK4Qf73TjLDYZKqMV7bxJaboGbkmF7XgGk+cV+qUYtKLNRiwOJCTtyPTnZ/0jxI2dWC0svLesrIPq4c35lpYHPWnbrUoezcygW+3X2Jg07KUK+IJgLuLI+c/7wSA/5SthMUmmco3mLYNgDm9/C3u89vwJrSrWpR3Vh4z6REfvBLOvyfvMjUXCQDlWYcKZplFzz5g5qbzFk7VmmE4OTqYApibf7GT66pE2jfbLgKw46NWj20Y/qULsP795iTpDJQdvZ6vtlyggo8XL+cSRXON/wYls2gtBi05bLFF9FpYHEsPXOf9dpUsGuhG5Qvz54imhMUmU3l86p58o/bEpDQLoW7OjlxSOxj1p24jJDrR5Kj3j7YeIvXzwAZ0rV2cj9YGM3/3FT5eG8yM7rWeqXxYj0ueXJRKy8O4ZPovDqRB2UJM7FqdzwPOsXj/NZt7k6VM3Ze8bHBDBiwOYkATX6Z0q2nzPeKTdSYFHvOEZYrilBpuEpNEuJWwkmMTOmRqdPY2ea+h8CTtGmDnhRA+XhPMj6/Xo2G5QqY1AVsr+6fvRNF13j4Autctybpjdzg8rr3NfFUpeoOFKEqJ/G70a+xrSsQXEpNoUp2KT04f051ZhIG92bXmUFWiElJ4Y3EQLo4OBF2PoGf9UgxoUlYxCDNHaHoencjdqMQcr4cQ4OPlioeLI9fDlWmFrIS12JvhaSg8absG2HcpjPd/O050Qgo6g+Sfd5uRpDOoMnuJpky7Rrs2hjjlNO7OjhT1duWGatff9PLn5ToZj7zsza41h2pGbJLOlO8pu3i5OtGhejFTzKhPPsuV1XyuTtmaDvhkbTBOjoIZ3Wtnqby9GZ6Gwn9h1wBbzz5g2LJHe582VXx4Lr8bPmokjCl+1NuNIl4uuDo5IqUkSWcwbUW1RVhsEp3n7uXHfnWpX7ZQpu9tb3adJ8OmbOHl6sSWD1vy7spjVCvuzVev+Zm2pKblwv0YOn6zh9Gdq7Jwz1XC45Lp3aB0loL/M2PLmfsEXotgwwe2w1g0NLJDh+rFmNi1Ot/vvMzMV2unC18yp9f8gwRei2Dz/1rS8Zs97LwQyuVpnW3+FkDZdpqZM5VSMmbdKV6tWypLztQeydOLUtaoXCwf/7zbnIfxKbz/23GrAilSStNe/bdaVeDohA74Fvag14JD7L0U+ljvHxabxLi/TvNVTz+8rOzr19B4VAY3L8eSQQ0Ys+4kG07ds1om8Go4gdci+KaXP1Wey2chFmQuc/ko/H70Nrci4vmwQ6XHus+zjOZQreDm7MiC/vVI1hl4e8UxknSWk+l9Fyopn0+bKUPt/qQNNUp40//nILadffBI7yulZOy6U3SvW5IGubQF13i61C5VgF8GN+Szf87w94k7Fud0egO9FhyioIezaW7TXCyo8vhHFwu6/TCeGRvP83Uvf1ydsq6LYW9oDtUGbs6O/NCvHs6OguHLjpoMKehaBAevhvN1r/Q9yPXvt6Bh2UIMXXaE9Set9wAy4o9jd7gZEc/IDpVz5DNoaFijRon8rBjSiGnrz7HmyC3T8ZazdgIQOLa9RfkiXq6cmNgBgKoTNpFgZbU+IwwGycdrgxnWovwTzW/1LKA51AxwcXJgXp865Hd3ZvDSw8QkptBz/kHyuzvzSp1SVq9Z81YTWlfx4Z2Vx1h37HaW3+v2w3imbzjHnJ65uwXXeDao8lw+Vg5rzNdbL/Jr4A02nLrH3ahElg9piItTerdQwMOFk5OeB6DaxE3EJumy/F6L919Dp5cMb1k+x+r/rKI51ExwcnTg617+FM/vTq1Jym6ow+PaZ3jN0kENeaF2cUauCWZl4M1M38NgkHyy9iRDW5Sjeonc3YJrPDtULOrFqmGNmbPlIiN+PUbtUvlpUcm6gAkourxnpyjTXDU/20xUfEqm73HpQQw/7LrCVz39sqVvYa9oDjULODoI3myltK6TX6phtQVPy/d96/JavVKM/fMUi/ddy7DskgPXSdYbeLNlhRypr4ZGVilbxJOqxfMB8NeIZpmW93BxMm079ZuyhXCzbadpSdYZ+HDNCT7pWAXfwp45U+FnHM2hZoEUvYGP1gQz7ZWavNG0bJavm/2aH2808WVKwFmbCdIuPYjhux2XmJNHWnCNZ4sTtyK5cD+GoLHtLASjM8LN2ZGLU5Vtp/XUbafW+G7HJXy8XOndoHSO1fdZR3OoWWDejssU9nKhr1luqKwyuVtNU4K0r7ZYJkhL0Sst+Md5qAXXeHZISFaEdya/VJOi3tbFqG1hLhbUcPp27kYmWJw/fvMhK4Nu8cWrtXON8ElW0BxqJpy4FcnKwBvMegzDGN25Kh+2r8y8HZf5POCs6fi8HZcp4uX6SI5aQ+Nx+WLTeWqVzM8LtYs/0vVOjqkppJvO3MGNcEUsSHHUwUzpViPbjtre0RxqBhhb8EkvPb5hfNC+EmO7VOXnfdcYs+5UjjhqDY1HZd+lMDafuc+Ul2wL+mQF8xTSrWbv4nJILDM2nsOvVH661Ho0R23PaFtxMuCLTeepWSI/XWuXyLywFfQGSXyyjrgkPXHJOpqUL8ILtYqzKugmq4Ju8n3funmuBdd4+kQlpPDp78HM6lE705Q+tkjWGYhL0hGn2vcfbzfl1R8P0H7ObgCCP3s+J6tsN2gO1QYrA2+y9MB15vb2Z8Ope4rxJOmIS9YTl6QjXv1rNCjLc8qxRJ0eD2dHPF2d8HR1wsPF0SL3zqMOtTQ0HpUUvYEePx4gn5szyToDf5+4k2rLSXrik3XEqvYda2bLcWmOSYli16pNe7g60bRCYQ5cCQeUJJN5Ec2h2mDsn6cAZf+xp4sTHq6OeLk64eGiGFERL1c8XVVn6WLpMD1dHfF0ccLd2THLK6caGv8FtyLiuRQSS+1S+Vlx6IbJfo32XcDDhVIF3RU7V+3bw8VJtX31N+DqiIujgzZVZQXNodrg2owumsFo5DrK+3hptv0E0RalbKAZnEZuRbPtJ4fmUDU0NDRyCM2hamhoaOQQdpMCRQjx5PNE2DdhUspOT7sSGtlDs+tMsSu7thuHqqGhofGsow35NTQ0NHIIzaFqaGho5BCaQ9XQ0NDIITSHqqGhoZFDaA5VQ0NDI4f4P1Oji2B/I8n+AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "proj = ccrs.PlateCarree()\n", "transf = ccrs.PlateCarree()\n", "\n", "fig, ax = plt.subplots(nrows=2,ncols=2,subplot_kw={'projection':proj,'transform':transf})\n", "\n", "triang = sids2triang(sids)\n", "ax[0, 0].triplot(triang, '-',transform=transf, lw=1) \n", "ax[0, 0].set_title('Original SIDs')\n", "\n", "triang = sids2triang(compressed)\n", "ax[0, 1].triplot(triang, '-',transform=transf, lw=1) \n", "ax[0, 1].set_title('Compressed SIDs')\n", "\n", "triang = sids2triang(expanded_wrong)\n", "ax[1, 0].triplot(triang, '-',transform=transf, lw=1) \n", "ax[1, 0].set_title('Expanded SIDs (wrong)')\n", "\n", "triang = sids2triang(expanded)\n", "ax[1,1].triplot(triang, '-',transform=transf, lw=1) \n", "ax[1, 1].set_title('Expanded SIDs (correct)')" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.10" } }, "nbformat": 4, "nbformat_minor": 5 }