Skip to content

Instantly share code, notes, and snippets.

@wtbarnes
Created November 19, 2018 23:57
Show Gist options
  • Save wtbarnes/2d271e91c74dbbfd5d11c8f383decbef to your computer and use it in GitHub Desktop.
Save wtbarnes/2d271e91c74dbbfd5d11c8f383decbef to your computer and use it in GitHub Desktop.
Create a matplotlib funcanimation using a SunPy Map
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Animate a Map\n",
"Animate a simulated map over time for a single channel and single heating frequency."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import os\n",
"import sys\n",
"import glob\n",
"\n",
"import numpy as np\n",
"import matplotlib\n",
"import matplotlib.colors\n",
"import matplotlib.pyplot as plt\n",
"from sunpy.map import Map\n",
"import astropy.units as u\n",
"from astropy.coordinates import SkyCoord\n",
"from astropy.visualization import ImageNormalize,AsinhStretch,SqrtStretch\n",
"from IPython.display import HTML\n",
"\n",
"sys.path.append('../paper')\n",
"from formatting import get_figsize\n",
"\n",
"matplotlib.rc_file('../paper/matplotlibrc')\n",
"\n",
"import warnings\n",
"warnings.filterwarnings('ignore',category=UserWarning,)\n",
"\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 92,
"metadata": {},
"outputs": [],
"source": [
"fits_files = sorted(glob.glob('/storage-home/w/wtb2/data/timelag_synthesis_v2/intermediate_frequency/nei/SDO_AIA/193/*.fits'))"
]
},
{
"cell_type": "code",
"execution_count": 93,
"metadata": {},
"outputs": [],
"source": [
"clipped_fits_files = fits_files[500:2501:3]"
]
},
{
"cell_type": "code",
"execution_count": 94,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(224.446,232.74,'$193\\\\,\\\\,\\\\mathrm{\\\\AA}$')"
]
},
"execution_count": 94,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAnkAAAJiCAYAAAC7GNT5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvd2OHNeZrvmt/M+snyyWKLUNDcYuSoMGGtNoUy2ezMkATakvwDTlC9im5H3ebNr7AvamrL4Ac+QLGFGiz1sS+wKmKWkDA2xgAxLpHmzNuG2JrCLrJ//XHKzvXX8RmVWVlWQVk+8DFCIjMtaKFZFRVZHv+r73M9ZaIYQQQgghy0XltAdACCGEEEIWDx/yCCGEEEKWED7kEUIIIYQsIXzII4QQQghZQviQRwghhBCyhNROewDLhDHmfz/tMRBCCCHkheIba+3/KHvD0EJlcRhjeDEJIYQQ8kyx1pqy7ZyuJYQQQghZQviQRwghhBCyhPAhjxBCCCFkCeFDHiGEEELIEsKHPEIIIYSQJYQPeYQQQgghSwgf8gghhBBClhA+5BFCCCGELCF8yCOEEEIIWUL4kEcIIYQQsoTwIY8QQgghZAnhQx4hhBBCyBLChzxCCCGEkCWED3mEEEIIIUsIH/IIIYQQQpYQPuQRQgghhCwhfMgjhBBCCFlCaqc9AEIIIYSQF5VutytvvvmmiIjcvXt3oX1TySOEEEIIOSXeeecduXv3rty9e1euXLmy0L75kEcIIYQQsoQYa+1pj2FpMMbwYhJCCCHkyCxiutZaa8q2MyaPEEIIIS8UV65cke3t7dKHqosXL8p7770n33zzjbz00kvy0UcfyVdffZXsc/nyZXn77bfl+++/l9dee02++OIL+fDDD+cay87Ojrzxxhvy+eefz9V+FnzII6V0u11599135ZNPPpEHDx7I1taW/OxnP5Mvv/wy+aW4ePGivPnmm3L//n25cOGC3L9/v/BLs6h9CCGEkJNy+fJl+fDDD+Xq1aul7924cUP+/u//3m+7d++eXL16VR48eOD3ERH51a9+leyzsbEhH3zwwVxjunTpkrz22mvyy1/+cq72U7HW8mdBPyJiF/2ztbVlv/7664X3e5Tjxjx8+NBeuXKlsM+nn36abLt9+7bd2tpa+D784Q9/+MMf/pzkZ2try/72t7+1165ds19//bW9fPlyYZ+y7VeuXLE3b97067dv3y7sc/Pmzbn/V29tbfn/ufOe27TnEiZenHHee+89uX///qkc+6233pKNjQ25cOGCbG5uyp07dwpju3XrVrLt1q1b8v777y98H0IIIeQkPHjwQH75y19OnVbtdrvy2muvFf7nfvnll3Ljxo1k29tvv11ov729Pde4NjY25MGDB/LFF18sPLv21NWvZfqRp/DN4969e8k3iGf1s7W1VfotJ/75+uuvC2pbt9tNvo0sah/+8Ic//OEPfxb1U6bYXbx40VprbbfbTbbj/1G+Pe/v+vXrc43l4sWLVkTs9evX7e3bt+fqY+pzyWk/GC3Tz9O4Ea21hz5sPY2fwx7yZt301lq7tbW1sH2e9bnzhz/84Q9/lvun7CEP/4/y/zuYSp32/+jatWtzP5zhAS8+/jz9THsuYeLFGeTy5cty9epVuXDhgoiIXL16Va5evSq3bt0qZPg8TTY2Nrx0vLm5KQ8fPvRTtpubmyLisoLKwNgXsQ+CXaeBLCmM+e2335YbN25M7ZMQQgjJ2dnZkc8++0zeeOON5P8O/lfl/4+uXLnip23feeedhRz/iy++kGvXrs2dqZvDh7wzCJyvr1+/Lpubm0fOtrl9+7ZsbGwc+Tjvvffe1Aeohw8fyubmZnKj3b59W0RE7ty5c6TjLGqfWVy+fLmQBl8WK0EIIYQcBmLEIWh0u13/3sOHD5N979y5I3fu3JFutyv37t2Ta9eunViIuXXrlly9epUPeS8Cb7/99rF8cxbxTQLs7OwUbrJbt24lN/9Z4MKFC7KxsZE85H322WenOCJCCCHPKw8ePJCrV6/6Wazt7W2fiDEtCXJnZ0du3bold+/e9bNcR+HKlSty6dIl+fnPf+63YTaq2+0uZjbqtOPYlulHFhwzcFrxeNN+EJfQ7XanBqjG417UPrPG1O127cOHD+2nn35qr127xhg+/vCHP/zhz6E/0yxUyn7wf2rWPvj/eJz/2dP2vXfvnr127dqxzmfacwktVM4oFy9eFJH5S5yclOvXrxe2QaqGWbGIFL61QNq+f//+wvaZxc7Ojmxtbclnn33m7WYwrUwIIYSclDfffFM+/vhjERHZ2tqShw8fytbWVum+Rw1B6na7helf8NFHH8l7770332AzOF17Rnnrrbfkiy++8OtHkW4XFZO3tbUlv/nNb3y1C4AHsfv378vOzo588803heNtbm7Ko0ePfLtF7TONixcvyldffSUffPCBfPDBB9LtduXjjz9eaOAqIYSQF4NPP/1UPv744+T/x40bN5JY73v37hUe0JCc8eWXXx7pOG+99dbU0KdPPvlEfvOb38jW1tah/wMPgw95Z5RLly7JvXv3/Pq77757aLmURcXkPXjwQN59993CzYUHTzxsfv755/Lmm28mgaZ5/b1F7TONt956S7a3t/1Yd3Z25MaNG0mMAyGEEBKzublZKopsb28n/3uuX78ut27d8v9jHjx4UBr3fePGDXn//feP/FB26dKlqQ95Dx48kG+++Ubee++9pHTaXJx2HNsy/cgC4wVu377t5+QRu7bI/g/7uXLlShLf1u127b179wqePnk5sk8//bTQbhH7TPu5fv26/e1vf1uIc8hLsPGHP/zhD39e7J9ut2tv3rxpb9++ba219uuvv7a//e1vk/8XW1tb9vr16/bmzZv25s2bU2Pjrl27Zm/evOn/Bx01hu7y5cv23r171lo7tc3NmzctuHnz5pH+F057LjH6cEIWgDFmYRfz4sWL8utf/1o++uijgkXIs+LKlSty4cIFeemll2RjY6P0W8rFixflrbfekvv378uFCxfkyy+/LIx1UfuUce3aNbl3755cuHBBtre3/Tezs5QBTAghhDxNrLWmbDsf8hbIIh/yCCGEEEKOwrSHPGbXEkIIIYQsIXzII4QQQghZQviQRwghhBCyhPAhjxBCCCFkCeFDHiGEEELIEsKHPEIIIYSQJYQPeYQQQgghSwgf8gghhBBClhA+5BFCCCGELCF8yCOEEEIIWUL4kEcIIYQQsoTwIY8QQgghZAnhQx4hhBBCyBLChzxCCCGEkCWED3mEEEIIIUsIH/IIIYQQQpYQPuQRQgghhCwhfMgjhBBCCFlC+JBHCCGEELKE8CGPEEIIIWQJ4UMeIYQQQsgSwoc8QgghhJAlhA95hBBCCCFLCB/yCCGEEEKWED7kEUIIIYQsIXzII4QQQghZQviQRwghhBCyhPAhjxBCCCFkCeFDHiGEEELIEsKHPEIIIYSQJYQPeYQQQgghSwgf8gghhBBClhA+5BFCCCGELCF8yCOEEEIIWUL4kEcIIYQQsoTwIY8QQgghZAnhQx4hhBBCyBLChzxCCCGEkCWED3mEEEIIIUsIH/IIIYQQQpYQPuQRQgghhCwhfMgjhBBCCFlC+JBHCCGEELKE8CGPEEIIIWQJ4UMeIYQQQsgSwoc8QgghhJAlhA95hBBCCCFLCB/yCCGEEEKWED7kEUIIIYQsIXzII4QQQghZQviQRwghhBCyhPAhjxBCCCFkCeFDHiGEEELIEsKHPEIIIYSQJYQPeYQQQgghSwgf8gghhBBClhA+5BFCCCGELCG10x4AIYQQIiLSNSdr38za961b7tijH3OePgg5q1DJI4QQQghZQoy1/HqyKIwxvJiEkBcCKGBQvuqZAjac8dcw37d2BAVvdIK/rrPGAjCmfCw47lH6oOpHTgtrbelvEZU8QgghhJAlhEreAqGSRwhZZuL4tWkKHpilfPWz9w7ra1Z/s9oAqHNHUQOnqYplbWuHqH9xm92JW1LlI08DKnmEEEIIIS8QVPIWCJU8Qgg5faZlzM5S/XI1blacIBS6XMmrZm3GttgmX1LhI4uASh4hhBBCyAsElbwFQiWPEEKeL6ZlCecqXZmylyt3s5S8/L18n/4kbHsydss/TYr7E1IGlTxCCCGEkBcIPuQRQgghhCwhnK5dIJyuJeT54hX9mhtbejAAnpTxSiSJrFXdclrCRdnUbGXKviCe2h3qNC2mcHu6ZJIGmQanawkhhBBCXiCo5C0QKnmEnG0QZL9RTbcfRAHuTzPYPbf2EKEqswxA5WvrsqXLWLWrZ9tmKXtQ9aDoTWy6HQofEjREmKTxokMljxBCCCHkBYJK3gKhkkfI2eRVVe66NbeEejJQ9eOPw7Dv01DWqOC9mJTF8TWh9uk6FL5KyT0ymXKP5EqfiEhPVb1dXX47FvICQSWPEEIIIeQFgkreAqGSR8jZ4dUo7u58wy3r2dfaJ6rg/ffBYo/9SnYcZO9SvSPgR6oqr+uyGd2v+X2aq3xlqt8IGbmq4EHZ29Plv43mHys5+1DJI4QQQgh5gaid9gAIIWSRQMGDeici0tG/dFBAoHqMFqSsIeZuNfvaPKSCR6bglTVd/ij6b5yre1D2an7pbrhY0Wvovi1tu6IBfWuq5K2oar0TKXqM21t+qOQRQgghhCwhfMgjhBBCCFlCOF1LCFkK8mna1Xp4D9NbE000gwVFbw4D2bKp2VoW8gxzZRrUkqOSJEboa9zTm3ovd/x/bHcDNyIn5VrBisW9QGhCo+radKLjdPX1dzqVy/t1+aCSRwghhBCyhNBCZYHQQoWQZ0+udqyrktfK5TURGaiE91gtU/7Qc8tZiRF5KbRWyVfjUaYMMqCdPA3yez1Wq3G/I/GimvmsjDURYxSpdQPYregNDEshKnvPH7RQIYQQQgh5gWBMHiHkucTHxqm6AVWjU3dv1KKvsKOs0Pu0clExMDQ+p38lW9X0/aSklL7epfJBniLfZiXLfhTF123U3U0dlGy3rOt9C2WvGv1eNHWfDpRBbbs2cH2di0zCF20YTp4NVPIIIYQQQpYQKnmEkOeSpip5LW8C6zY0S2LxkFWLMKVKFme3Gil7aL6CfjMFD5m5/Ui1QzYtTY/JsyTOyMXr13R5ruFuxpU6YvXcej3KyIWq10A8n96/bW3TqYcbekXbfz8oHpucXajkEUIIIYQsIVTyCCHPJW39igqlreFjj9yyzDigYpxCUa9o/FLJX0CofHmReMTxDWeUREOc4DyKXu6/F2fxNvU1RJhcTURM4DA6bp8l1V5IvtHMWNHlX6qit6Hxdp0oI7etrxGvB0Wvob8Xjej3AzF+bVX0Vvpu/VtV9nifnU2o5BFCCCGELCFU8gghzyWInUMWLWKNKpnaVUY9U8ZiMnuxqRm5cegfVMW2riPGD8oaVLVmSRvE/qGaQTWLG4yP7RW8zIcPx6GvGclBVuyres+cj34x1rP7BUo3FL1aEr+HmD63RIxfp5Z67P1x6JtQ3TsDUMkjhBBCCFlC+JBHCCGEELKEcLqWEPJck0+vgnh6tTDVmn29xTRV2t4mbTHLhanePDGjbF8sMesVTw83q+mykc0dD6JpNUzP7qltxSNdcnqWHBVvpHwQtv2lbhvpjWstbih3czeimAS8zm1X2mqzAgNlLEXCFC7v19ODSh4hhBBCyBJCJY8Q8lxzlBJl+T5Q7qDo5cqea+P2QUm0inGdtIu7+v6x77AQ0O6WsVgXjm2SPqDg7UYB7P+uwfNUQsgiQVLGj/x96+690cRtWG0WlTwkY6BcmjdSVkVvpREZKGvG0ZpmBp3TexqWPzAR70e/n0zWWCxU8gghhBBClhAqeYSQ5xKoAT4OboLSZSZZj/cBME6GGlGLAvsQngcz5dxCxWaqnYjIIIttyqlmdi8xUO4Qd7ejasc3w+K+hDwNUKLMK2t6T25Ejt9rqup1Gu4mbmrps3otXTajUmhe3RtMdOm27w3Te34/sgRa1THsslTgQqCSRwghhBCyhFDJI4Q8l8AA2JcZ06UxaVZsjC9ZpgpeC2pEtajkQRFElu144tYHqm6MoixCqHEYS56JCwWvzOAYagYLv5PTBjGff+q55WuRwnZO7/vuKI3XW9Gae0HJCzc54vdC3J4qeqpSHyBGL1IM9/X+7+pYNnQM21C6qewdCyp5hBBCCCFLCJU8QshzCVSHH2ZKHtSyOGYOqhli8Zrw+Mrii0REKpWsPJo27qlyByVvEKkceZkxlCiD9x2UvFhdhPrX07ZnScHr6rkjuTK6PAXy0m1UWk6X7ozPSuR4n08cF/pqFq8H9W2o61D0WiVKXojfc7+Abc3A7ejN0xuGQflt2n9Lfy+wXI1+177Nfu9IkefqIc8YsyEi7+rqJRH5L9baL6P3/1FE7ovIBRH5/KjvlRznyPsSQgghhJxFnquHPBF531r7noiIMeaCiHxhjNmy1m4bYz6W6KHPGPOZiLytr6e+l3OcfQkhp8+u/4YPNaAoZeSxeFAWVlqq5DVC5Eo1qzwx1EC74cjJBlAu4ooU6L+tf1FXG6gOkMb19SO1zmfTPmMF71VkFJcoPtiGEKnRDNVnmL1HBe/pc5hKJxLUVwBhDZ/tRvQePt/8swSxfx2yXYcaOzrIlL01/T1ZiwbQbqQxeSETFwq6KnpxdRcdzIEq53u6bOr2RknWORW96Tw3MXn6UPcN1q2198WpbVD23sgUt/vGmLeO8F7OcfYlhBBCCDmTPDcPeeK+gLxfsv0lfQjbzrZvi8jbs97LOzrOvoQQQgghZ5nnZrrWWvulMeZvs81viMgNSRVo8L24uL1Z7+UcZ19CyBngiU7VnNNlmV0JttV9ILjb0G5Vk6VISLyYaJaEVngSa8e6TPsUEenonBhsJRrZXOgBprjGIRvkkW572lNNr+g417Jp2rKp2HyaNp/Ga0fn3II9DKaqZ5SqApjJa2dtsaweYToSjEumlPH6YFLcfxpnOWEkn57NE2GmTbPOIr418RlWj3D9MT3by5YjvY9REi0OY1gdwUBZl0jOaLj3cyNlEZFGXZc1d4BaVZf+97p40rgOLPtX5HlS8iRLpHhXXFLE5yKyOaPZrPdOsi8hhBBCyJnluVHyYjTL9qq1FtOoD2fsPuu9k+xLCDkD4Nv7X6gi1qxqgHYkS0DVq+rXWiReeCWvXff7wgx5oB4plUoqtSGZYrURtiHAHDYS6APB432VmLb7oc3TLlsGBS9XzUCZopcra11tC6GzHf3HaKFgfWYPM/KWNkVTanwOMJquTFGO4u3T9smP57bZZNswU3ZwS0Rik+8HVjZ9jB+KLexkSiQR9NPLVMwycvWtTOkEhyVPYNma3kVBqc3bxuC6NKvpMiY3786Nv/H+wShuk34eY5XFYSwOC6NYyYPtSkUvOJT1amWSLEWKql5ff7/Ooip7WpzJhzzNaL2AdWvta9ku74vI1Wh9W4pTrS+JS8yY9V7OcfYlhBBCCDmzTH3IM8b8sxSTEObFiIi11v78KDtHCl3ZuP5RnJXKtq6/Ya393BiTT7VuiMjHs94rOe6R9yWEnC1Q5LyZxeaJFEucQR2o6U6NRpAuLNSHUSrLQJ0LJdHCATrNNLZI3VZkoH08Vtnm/+4d54wOB3YoeXxcDNxh6tmyLHYRChsuR7uemtlCsRQJ5xpsYtx2f/30Gkyiix/2wTpKxqXbY6C+muzcsG8U5uiNq8eZmojzQl/xcYbe3NetI6YMY63561VUhnGLoA1UrlzditsAvDfOljHT1Dgob/USdTHvdx4w1lqJjIp7o0ypndYPrtPBMG0DZW/WfVXDvaj7WBtO2tpJMpbzsCV6yir588QsJc8c9aHsKBhjPl1AHz8TkS9F5KFO2V4QkTd12+f6wIe4vQsaryez3jPGvCGSxPvN6ocQQggh5Llg1kPeoqs8nKg/9ckrU9Sg+l0TkV/rfpfEZd3KEd77uTi17r0j7EsIOaMgxu2vfMZgLDE4FQAKzySTH1KVKFUvoPoYA2VBVbtqUX3APiM9EMo1bQ+OdSoFoNh19S82YuRmqXJQYRpeeUzVOGQY10rir3xhed0Hy9goulIpv04hKxmqWqTkITZLr89I5abhyCbrsdIWYrLcOj6rMiVvkil5UAorJlWFYqA49jTIbDBKDaxrlWJbqEwYgzfILil5589jSmzhyI813jY7trByDCUvj0OMj5OrfbOya2vZsVqZbFx2ftNUPtwT/WG6LhIU8prG1uJ3KvweRuPN7vGOtnlFrz+zbWc85FlrfzWroTFmS0QuWmt/n22z1to/HLe/w1Dz46m3oE7f4oHsk2O8d+Oo+xJCCCGEPC/MlXhhjPmJiLwmIh8aY+5ba/+riIi19oEx5qIxZgPbCCHkWfKtqmavJluhtLg1rxypVBLH33l1RuWTsVeX3BIKXqzsoA2UKJRm2tfl7jG88P5Ss3bXQsKvz2qFetKEvxi876C8ReoK4uigwqF0GzKKEYdoImlkmsKZq3bxvvkSbcIy6t9gH2Qh49pOdOn2i5WdytT+3D7VRGFKVbiKpCogrlOqBqFfzfSsQg1M2zZKMkDRzyhT8nw8YiRJVjNFspIFGQ4jWc0rnOP0vUJs6Qz1LN+3bDvGV+y3mP1cy5TUPM4x/5xibHbsfH0U/X7sTZBFW+wn3zdcFyi2bnu9ZAwvKvP65L1prb1jrd3MH+astV9JlBlLCCGEEEKePc+VGTIhhBBCCDka8/rkvX7I+5dE5PeH7EMIIQvHG6FGyQ7/s07fdGCTgQB5rdXU74c5IEw77qur67665I5mTLlieg3mx3vqqAtj2DJrE5gVn9O/wptNt1xXO5ZOZCsRLEyQCJFPxWIKNnxvr+tcLmxictuYqs6H2WhKMZ++xrXAPvF0br7vaJQmWuTTuCLF6Tpct3E23RknXqC/MHWYTjGW7TvMEjgwVYqEjHiq3WRTk5huxHQh2sbTkJVK2rZusvUqBlWc4q366f70pkiu7ThNRJllMePb2/Sc8ynYPEnE9ZtOb+fEFir59QD59SojHz8+J2PT90XCdUBoRZ6YEjsbIUElN2iep9zbsjKvkveRMeafjTF/E280xvydMeZfReSzkw+NEEIIIYTMy1xKnrX2K2PMByJyRzNqwbaIvGOt/ZeFjI4QQuYkLm0Eo+RBphj1femyoE4MVSp4suekhMf7KMWklhpegYnajFMFD+bHUBZWVDn66+gv7qomVmy0XD/rWn9sVZedVpCbkCzRbLoOkDRRh92ESlPVyOeiUkm/w8O2ZDJGkoNNtosEVW6oEf/jcbpPrDZ59Q3WI8PUFmWSKTDxMXP7k9zEuIxcbSpTtyaZUoT1agUKX9H+Bu/lbYNlh3vRsKENjo224fzSMcf3VcUrYegX65nDsYjUa+VqZUh6kAKzTIlFwmc3ShI80vdmJXbgXKaZOoNYGS5+vunnHJJBbNTeLQf572ym1sWvVZD35eVmlYx70Zi7rJkaBL+uD3kXROS+tfbBwkZGCCGEEELmZu6HPGPMj0XkZyKyba39nW77OxG5Z619vJDREULIMUGsW5mNAmLkdvUrf3sf8XbFeCgoeDsHKJ3k3m/VUisVkaDk7WpMHo4Ds1qodqtRnB2Uu67KfBvqmbKy4pbtdvjz3Gy5bfWG2wblzvh4MRjGRooIlBuVa/JlsIgpU/LSeLsQZ+d3DbFTui9sY6YZEsfbckVnWhxZ3j7tt9gmxG2lbbwZsipvjWp4PzdXnqZmxTGZUPByu5DcViSOUwtm0Tiuxhpm1irxNh8niHW/TI9X+t4hZtXx61y9LDMezq1T8tjL3NhaJNwbIR4QS20zQ7nFe7khdGzgjNejbEkCc8XkGWMuisjnInJeXLUIERHRadpL+gBICCGEEEJOiXmVvHesta+LiBhjLsdvWGvvGmN+KiJ/OOHYCCHkULqqMLyamQjH3/j3VVnTMDt5eAAlxkkLayXpeL3c0FhLMEGdq0fxWBBHNBHXh1etqJy40XbLbid8r4Zyt9F1A19fb4mISKvj1qHaiQTlrlLVZVZWazJJ4+3c+ZenA3v1ZpyaPcfv5bFzYRmrcuk1y01+oSjFCtjYqz5uPVd0phn5JuPP1LoyJS8v4VVR1azq1ajQZlq5MeDLmkXXaVobfxw992aJgXLFG0KnbeOsVV86r5Yab1f99cfFDmOqQkWspiXooPKWmVPnY8g/51n74l6B6otlryRTPcQ9pvF1gxKT52mfYV6CTaS8rJuIyC7LmXnmza69F70u+3XcnLNfQgghhBCyAOZV8uIHu7LvNIf56BFCyEI4r3/FzqnPHOLe4ngfKAjfq3fe/zhwy57KAxuDsHNL++shfk8VvG1dIu6nzPsO287rWFZUXYSCt9kNtco2zznlbk0VvHbHNarWiuXGvGI3dhKInWSmbraoDlVU0anjz7yPs4Li404wjsnLve18TF6WZZvuiyGk6g/W43hHqDz7GhMJX0FkO08r3yUS1LOg5JnCvlD38jguLJF1GStxZepeDPqcR/2L1U5NjC74yYVdIlUOanFWri4ofBoXWnYTot9JquBB4avF6mLmoQi8d2DFFLaF/jWr2meoj3V7OA9kqvtx67I/Sj9nKO0iQQ3P4+ww7FrJZ0emM6+SZ4wxv9DXNtr4E2PMP4vIpyceGSGEEEIImZt5ffLuGGOuGWMeiojVJaZo36VPHiHkaYMs2o7GK9WyLMVYnGhiH932bz7L1i03opi8df2rCAFqR/f9Ngtx68b9Q8FToU5t7KSFShVN9b6LMmbb7TRjFspXv+fkxuEgyBuItTMa6NaAX14TfbhlLSrlkCsvdoKMWfXA0+UkUedSPzy8F9pEMX+Z316I/UoVpFHUpt9353Sg6cdQ9vqD1Gsv9svD51jJFTBkX0ZK4TDzQcQlPBgiFqwoEdYO8X/LM3XdvmmbXGXM+0j7SQ+E7WUegVq0xJ8PFL3mBEph6MuPd5Squvg8alV3gLgqCh4BQhWUPIu3mL0bqGj/6cDLqmegqVcmNbsZannsffdIt/1pSlxd/Hu3mt0T+DXeKQ7hheUkPnkfGmNui8iboj554uxTdhY1OEIIIYQQMh/zWqisG2PWRcRaa++fOXZ0AAAgAElEQVRaaz8UkS0R+ZVm1hJCCCGEkFNkXiXvN+Ji8T4Tkd8bYz4S55f3KxHZNMb81Fr7+wWNkRBCPJiuOad/vXIbhWGJcwimi1a1zas6FYQp2J0o8Fvi1zPYiIrco2wZ7FtaNUwdp1YnsV0JSodVDtz0LKZIDw6GugwDwZQopng3znVERKSmJ9ReaYuISKPV9G1gt5InQvhyZiPXP5I53GskWqRTuqPhKFnG780q5RWfV9x+oNO2gwGWqUFz3CZPHPDH91PJsZmzTsvCzkOngTEdHMqRhX5q1XR62feVTf3G54nPNZodT/pHmzjxIjdKDm2KSRrB+kWSseVjzKfk3TjT8fppc4wjcgxGskToD4kY6COaDs6ScnD/9nTOdXfP3be7URbFAcr7jdPrkZdnixMoykzMY+Kp2J1ylyASMe9D3mfW2jsiIlrW7Gcicg6VLowxVxY0PkIIIYQQMgfzPuQ9il6/JSJfZaXMHgkhhDwFEGzdUhUFSRXBaNUk6zFQ/TZVccsTMY7Cj/Sv5kuNsG0lU/BggoshHPSdxFOFH4sERaSpJ+BLPw1RfixWXFK1qaYR+c22U+7aqysiItJot32bai398+4VHm+cDLUuVueg2A113e1TrQ50GWQoqHK5gXKuHJpK/DnU9HzU3kWTTlpQEJHgUZKFUKkcruTl2/LybLjGsdoIS5FcgYTqmpd2ExGpFQyH0zZ5QolIUOegHOJSzlLj0CZX5eqaKFEvsUPJEy7y9Wokm+HzrNdhr1It7OPHn9nq5Are9hO3/uQgyGswFM9/FfNyZvGlX9Xf5x/pNq0qODURg8xmXguVmJ+JyP+ZbTu3gH4JIYQQQsiczKvknTPG/J2I/K2IXBKRqyJ+6vYfReTWYoZHCCEpKFl0PtseSmNp/M8MGwUoeq05vuYi/m4l+Bp7BU/FFD+Gngp3ubWHSFBGEGPY6bgO22qzEis8UJlaGpPXbDkZEWqdL3dWDUFilVzJUwXPTgvwmrZNROykloxZRMTA1sUijk77z5S9uA32gYsIlKR8SJVIMcz7mxUDWKm4869D6WymY5hkMWHxGHK1FMeBclUWJ5ibCWOflip5iDWM30OboJ5BFSzejFAmJ5k3C44Xt8n7MZkNSpliGKxS9DyyzyM+7qDvbmYozVAtEf849HGIoX9v4p0t8buKZSwcwuEF2/A72kab6FLk5ctonVJkLiVP4/GMiGyLyJa19rEx5qI4Ve++OFsVQgghhBBySsztkyfuYa4LXzxr7VfGmG19/WARgyOEkBx8W/elqrIMO2yPwqG8KpCHGo3m+Obfz1QIkbjgultCXKppPBoUvlaUOtjSWLzNTRdHd+6lNRERaahKl5Y1S+PdanX3pxuqTVDRwqDwOjcphquwN/SNYt3wD2HSdLF+E43RGw5cTF6t3w/79t22PPM2xNW5pbGxApbG6+XbwzkXFU9vyFxSYs2PKVOzEFuGPtA2Pn4xGzVVwNBnrGrl48+VNSwbjWqhDfaF2oi4RKzHY8DnGoyn0/NNYwtde8RrVn18nQ/+kwJZsF+hdFmUeY3xIzMaINO41UCsYbifcMkQz4gydvuD9ERi1X2WAp8DE/I+FbypzOuT9xMReUNEfqevRcQ/3G3E2wghhBBCyLNnXiXvTWvt70TkTv6GKno/FZH/eqKREULIDPZUZMA3/7yg+UFJNl5bv9YiKbE3R8bedypkrCaChpad8mWb3NaOqhttDTTqrga15uXzTsF75YcuT6173i0brVbhmD6WLVNeoML52LwoDq9aqyfbCmXHNH7NRMpLpZru6z31hqraDYKSN+z3dNsgWQ5V7RsOhro9ZBT7jNxxmkWbl0+TQZzxCzUr9WfD+pQwwuQ88qzbPMZNpFjKK4+ziz0OsS3fB2octpcpeV6NzVTHOA7Rj8HiPXz+WV+Rkod+oPJiGeLspmfxilczy1VNN6Y0/hDjhyLts5RLjgPPwT31fny0qzGL+nHvR5J6/vuM31GULGOW7fFYRHYtIYQQQgg5Y8yr5L1+yPuXRIQVLwghTw0odbliN+ubflfVAFSrgKL3alS54NsjuujHlTY2Wq6jdZUK1zquw46a+TV0Z2TOioisrTslr9lxyh0UPHjdxZmyQX1D5mRazcLHkUVtqnUX21drtJL1ChQ+Ve2qtWD455W8XNFDbNho4PcdD51iN+ofuOUAyp4qeb2DZOneQxwffPhSX76hxvkNeuE4/QPXX08rg+RZnnGcWu4Jl2en5ssYnwk6SRW3sjZB8UpVv2o1VRdjJSz/rGymyk2iGEP0W4hV9Iuit18lyxLOM2fzsceDCOeYxnaWXSfE+uFehqKXxxzG44aC2vC/NO6zHE/cZxh7WlZV+O3p7yEUvTyTlhyNeZW8j4wx/2yM+Zt4ozHm74wx/yqu3BkhhBBCCDkl5lLyNO7uAxG5o954YFtE3rHW/stCRkcIIYQQQuZibgsVa+3nIvK6PuRdEJH7tE4hhDwr5gnAhv3KD3QOY0X/AsZllc5rv7sa+/9NyBtwbdUE+eVOaPTDc66jlzbc1Ofamls2m4f/icXUZG9/P9kelyirN92Uaw1LnXo12VRcvI5p2kZ71fWhy2rD9Vutu/cxRevaw1S5nqyDyTgqy6ZTt5i2HQ/ddO3YT9vqNG4/TNeGfbXtCMkZ/WQ57PV8m4G+xnKYWbfAssW91ulf3TYa5su0FFsZhWnUzIw5JiQh4LqlU6Zl+GnU6qzEizRZIrTJjxMnzaTb/NR+Zu4c2+x4mxtMFes9H6ZPo2lUnZZtturJuIO9S1aDTcK1C8kybt+VobvnYKCclDXTDIt9vdVasFvRGfydY5QgJHM+5Gm1iw1r7e/1wY4Pd4QQQgghZ4h5lbxfisj/tciBEELI0wYlxKDgdZtpOTKRYHJcVwXhR5klC9quNEOjNd3Y7ToT4e5GR0REGk2nekBZQtJAvG3vsVPwYCsCpWolMqKFKlNraJJEnohhoAYF5Q1KHhS85uqmrq+LSKTk1ZrRcarp0iCoHgH5UZkuVfUmUPCypVf6omQNbJuMhsl7dowEjFTpc6/TbVD7xj55Yxjtq8og3ssSPKBcpUkO02xEdIlkhBmGuyZLcgiJEkVFLyTJQMkrsbKppJ+nydb99hm2KLaQRAG7mvgzVIUNptd678H+xgxHhX1zdREWN8GsuqgU5tY1NbUaQmJSNZLy+qrktdTNHG8Nte0rkYJPO5XDmTfx4jNr7T9Ne1N98gghhBBCyCkxr5L3kT7I3bfWlpkevy20UCHkheQV/ep4Fr9lI55uTU2KVxowcA37DNWNFYreOoQPFTVg1jqIDFyhVHgTWVhq1KrJuo1iwYbDNH4MKl+tDquQYDzcUEUqWFyofQgUPV+UPih53iql1kiWUPBqrTVdb0dtmtqP/muwUGBGyXFFRCYjN75Jrv7pcoLjjWM1SM8Dit447bcBS43oOHgP+05TAd22obaZJG2CuoW+SvqfpDdsrpIlJeMytS+3uPGxcma6jhLUulrSNt7mlVlTbmg8ia8t4hszxROl6XBt4nP3CucwDTwNti5h35EqbOMs3nE4LF+PXw8GeC+9xnV4GEV6kzFq1zOByl6M2yNHZ96HvAfiHHs2jDG5iG3ERWv+x5MMjBBCCCGEzM9JHvKuTsumNcbcnn9IhJDnjVciwQKlw17R9bOg6GF8a6rktfQvHxIao8RGr5a0dNwdle6g4MGk9fv9cGJtTflrNlw8Wl1NX6HcIf4qztAMJcncst5wgwtlqep+36q+RqkyLI1XfFDWLLTJ4+okU7Ow7lU7CUpepdrQXdw5Gih5UXatmaFSxe9PZu2HmD8ok7osa5HvM8nWk31mBdBlbYKaqOeaxdPlcYnxvvE2EZGKz06GOhdlJ+dprj62TdtEn53Pbq5m/6K9woaYw6jMXG9Pl/var1sfHLhlrnLGr71C7DOKUUou7IvM5H7P9TPoaxxfruiNQhsod9iG8ma5wfQoKhmHmLyDAbJs3RK/d/3ZHy3JmDcm71rZA54xZl1f/pf5h0QIIYQQQk7K3GbI+TZjTFdE3jTGbIibri3sQwhZThK17gwodyKpungOyl0W14M4n3r0htaY9++FWCAnITxRMetREFGksoPMQ42DUkVko+fUjpZKh9XYCy33SfNxfKrkNUK5sVrDKWzwy6u3XPauj7Pz8Xex+qdZsyjP5dUfpzYarxbNUP9yTCJ56nloe6haiB9D/1F2rX/Px+9BXRqXLkUiVc7HvWlbHC9S7WzJtrhNYT8pxgWG00v95tL2RU+4eF98LvFxp11T700YlZdDqbn4s3FjRAyd3mfD4Cfo1VztD+eDcnNQ7YaDONsZ8Xoa7zjG9Z/ueTfycXYjXaYKXqzkDVTJg6IHXzwsEdcK9U5EZK/v9t3uuW0P9ffskYYf7lDJOxbzKnkeY8y6MebHInJORL4Rkfsi8vOT9ksIIYQQQuZnXjPkroh8LCJvicgj8WWTZUNEvhCRGwsZHSGEHMKrKpasVqfvA6EOoT+jTG1s1CIlT183x8icnCRtxqpu9CPh50BVhkd7bqeqKlRQMDpazB2KnohIs5kOuKYSYqOFLNKQOVmINVNVqKqVL2rNjq4HzzsobLmqBBXLZ8cOo4oU2T6FeLRYNcszbr1XXz1pGytYiHezkzTmrxBvF1fWQAapTePGooGEl3nc4RRiJc9kylcYa5axHGe25mMoxNcdQ8lTBa4a+RVWfDWSVMnDZ4Z9x5H657OoszbwFxz0DnQc4TzyWDw/pkqaLZy3c22nxNdFWec9VeX2DtQTUtdDvJ3bby9S8qCU7+rH8e2Mj7KbKfNU+YrMq+T9SkRuWGsrIvKutXZTfyr63sOFjZAQQgghhBybeR/y7kVxeefiN6y1d082JEIIIYQQclLmtVCJRdEvjDE/tdbG5scXRKTMJJkQQhbCj/Sv14bOTqlriSBmfFiSADLJpmsxA1WPpmubDZ121Gmnhi6bNbdsj9PjiYi0MMWrS1hCPNl3O/e0RFO7F+aeOm03lTf208KpgXKtEabdkIThrTl80gOSNNq6HiVe6DZM4frpwDxhIsJPm8K0uDCdGlWH94bA2YXO1sumRvNECICp2cpkxlR1fryyY+Vjy9ZtMu1cPh/oEyJKxpqfs983nyI/gn1MsL+Jpmu9KXU1Gf+4kpojJ1PIfiocCRduehZJO7U67qHIMLuK6XnY6aSl1pJEoSluxLiUuI8PojiGx3vu9fc694pkim3N/UAyxVFsljA1uxpdUvza+hli7YfTtoF5lTxjjNkyxnwkItsi8p+MMf8hSsJg4gUhhBBCyCkyr4XKHWPMFXFJF9+LyDsi8qmIfChO5Xt7YSMkhJAIJFq8rMLHSj0tTTZQRaESCTRlqp5IKIzeiGS5ei1NNqhA+YJJsv7VjIWNlaaWSWtWtD9VVaAqeuuIoFDBXmKSqYpQT2CKLBKVR/NqTbHYvEhWGssH/9fSdahNXnVqRu3TxAuv6I1TqxP3Ok28yJW2gsIXj9MrUeXlusr29buYWrJ9VumwfCxHMUuu5GqjN5OOEzzSxBSvUGbXNlbygqqYqYDYNzalhsrqzaJHSRvJzid+Pc0SRrxSHI5TuH+yz6EyHBXew7VDUhEsVPY1+2hnN/T55yfu9R/3XJv/R+1QjqO05QpeJLp7BW+XCt5U5p2uFWvtHRG5o6uPReR1Y0zXWruzkJERQgghhJC5mfqQZ4z5B2vtPx2ns1kPePP0Rwg5G+Db9LP+pozjno/+Up1X4WmznapnYKhKXi+yZVCxwcfkrajzRLtRVPJq6rcCkSMvfYa4u0Y0ptWWe3NF4+zQH2wl9nuIzQtjOlA7iUrFKSBVPW6tpqXRogOg5Fk1i6GCauOXkfFwXL4qpobrVRKTFxrnMW1QjoqxcrnK5+1PSmL2DiuFNpOszJiPaSspHTbVriRXBeP3pih4uYIlEhkx523nOT9v+xKu7STrH9d01N/T5a6IiIwHwf4Gnz32gXUKbFIQ+1dvtX0bxOSNh2r5o0bJZffOUFW9gwO3796eW+7suuVjLfP33W5o++/77tz+W2QcflxUJBf8Oke/1meiZOJZZ9YduegpV07hEkIIIYQ8I2ZN1xpNrFgERoJhMiHkOQGlwdq6XNVv0bMMSuchV+y6KjKph7B06uHPx1rLve52nAoB9QzZf0NNnR1EgXiIiRvpuDXETdZX3AHarSjjEEawyNwbp/Il1MBaVAqtpYogMnNbzTSGLmQmBrUGY4FChHJQfc1OPNgP8kez5S5IvYkyZoOk/4mP0YvUoDFKVrl9m6r61NtrunSlxmvN1XDukbmuyNFi2PJ9JTMkTkqI2SnGxtl6Eht2iDpWtq8pqH4maxOZ/EK5m4zL25YwLd6waKA8vQ8/7rK4xKwdYiPHg30RERn13XIclTWDkhffAyJBwWuoghfH5I0G7vVAXH9Q8oYDd7z93X2/7+Nt9/q7790xv9t2+zzaSzNnH0aq3R+Dp/WRyQ2OqdadjKkPedbav3+WAyGEEEIIIYtj7sQLQshygW/QzeibdDvLaBstOCbvrzS+7pwu1/XgiJVrqzLWjuLuoNytrjh1q61yH7JSoQaNIwWu33fqRk9j4xAr10D/7ciTTrNroazV61DnxoV+QSXL0m00qklfiLeLxZq+eudBEUQfGNugHxSZ3oFTWBotJ5P4GL26WyJLchQVn+/vPdH+tDh8z8VqNdpOucsVPRGResu9l/u0lWa/5mXL/HKStBUb2hrEEkKdsWk5M1+WrHg0z1QvvGhbHpuXx9sl5cZmnaPMVvSmja2g6ElR3fOKoR9TrbBvoc0hPoPxsaHgofRdiJkMMjxeDw6cyru37cLqt7/bFhGR7/702O/7x393St7/+727x/6ssXfbqtw91tt1O1L5jxrDG6t3+PtDBW8xnCAKlhBCCCGEnFX4kEcIIYQQsoRwupaQF5yyaVpwkE2Z9E8wXYsyZD9ohW0vr7jvmd2OW66208QFJER0OmE6tdNx008tXTaa7r0asimQjDAOg+9pEsPuEzctheQGTKe2OyHhoKHWJWNt325rEPq+m6aC+WvZtK2f2lXH5HpdLVUaaXJIPAbMEuZTvJi2FRE52HPjr2l/MEr25c4qmKIrTuN5K42R2mQcuGnbav2RG2NrJYwfU7k6bYvSaL5sWpSYYbxRb26oq9sFU5fRmHQ6FcbSFp8V3vfJG9GNl02XGptZmiSztdn0r7jlTDUDSR8lpsfx8ZOp12w6uJCI4ZMq5ss3DAk7esxsehaWKpMo8WKs5zwe9tPlwO0z2He2KwdPgtPZ3rabln38vbsXHv7Zvffdn92+f/4+9P+nbTcf+yedpn0EY2Odpj2OIfGsvzkn+RtDilDJI4QQQghZQqjkEfKC4795L/gbNOxX/kLFn1c67mv7y2vhu+VL6+5P0Pqa2wlJFEiEaLV1+0qQ/1od9xp2IlC1oGJBBRkNg3/D/hOnXhlV0pp99x4MhztR/+gPSlpv36kZVVXY+j1VxCKLFqhKUO4wbqiCSH6IC77X66nVRXVKAoZISBwxjw90X1X0NAGj0W7r+cVlzVBiSw2UdQwjVXaGPRdI39sNyk5V26CoPZS9hiZp1JqdsG89Tc4oKG4lViQF0+B0IbYkoWCapclkkqp28euCupgnXpQobF4ZzM9D2yCRoWwsxtSyNnpPlpxPQfXLbGTi8fv8FCSsjPq6VOPjyAy5v+tUuf6eS5YYqGLb23Xr+zvuc37yKHzeOw/de9sP3b47O+7e2H7s7vGHT8J13D5wY1DvY1HnlCMpeFDuNvRy5Ilc8YzBPCXPWM5sOlTyCCGEEEKWECp5hJCF8qp+W/+himOvrLiv23/RdX9uXtoIcV3nzrmd1tbdsrPqFCmodc22U4ugVImI1JtuG+LRavXUwBdmsP39/WgbFB6Y7krSf7MT+odKBiWwjpg/Vfh6B05NGQ5i4+FJso83L1YlD8fD++61U2PGI8SuJafhLVxERIZDGCTDBNkpLxVVBmFwGxvdwjoDSl6sRImIjCcogxXirqD6+Ni5rBRWXBILsXw4TrXWSNoYb+HRDGPSfYyWVCsoaiWlyoJSmO6bx/WJREogVLFp5sgRk0kaG+n7h+qH84na5jYufqyFUmgl6hxMizO1cZbqlyt3wwNnj3Pw+KHf92DnexER2X/slLqDJ0903S2fbLs4uyc74ffiyRN3L+/vu3t570CthtTeZxB5Jg1QLlCHDauUXEWL7VCg3K3osp5dfvR1FGum3CSZHA0qeYQQQgghS8jcSp4x5sci8jMR2bbW/k63XRaRf7XWPp7RlBCyZMTfsv9ChRsoeD/YcH9mXt50b7z0Uojr6m6u6LIrIiKddRf71Vpx22tNp+jFChXixWoN9x5iz6CUIBYpZjwaaRunJNV12VpdTfqMGQ2cytFSo9j26oFuT0s/uW3uda7YYYmYqkY/mBUjxm+o8XbI5vXxe5G6gfcQB7i320/6tb6sWdHotr3mzI7rLXfdEbcHBW4SFaM3moE76Dl1bzR06o/suHivapS9W6mlGb6I50O/UFjxGbrXmq2r6l6IG6yn61EWr1f/MqVrWvyd7py20WWefSsSSoblbQUKXkkJM6jFaOvLsWUxf3E5OK/G+ZJk/fQ8on29gofSdIijVNUV93h/b9e36e1pDN6e6/9AS5Lt7x7o0rXt9SIFeoIycnrKvqRfUTZDiGiuuuF3f1UvWzd6qkC1wCY+BsTiTdI+a9ElfuUQ6SnOvmUs3uHMpeQZYy6KyOcicl5ENrDdWntXRC7pAyAhhBBCCDkl5lXy3rHWvi7i1TuPtfauMeanIvKHE46NEPKc8GoUFvcysmjX0xg8xN9BvRMROffyORERWT9/XkREVs+9JCIi9bbGe3mlJ441S5U8qCdQRrAeq1qS+Y41V5A1ulLoH8pW8BsbJEv0Ox4FBWjYd/vGZcVcv+lYUD7KnZt772DPKSz9HhTCojwRYsr0eBqj91izIREfOOiHMQ16bkxrm26fzvp6cs4+kzaLaRSJFEJd4vyGvRC/N/HxjVCDcI3TOEEofSJRPKWqe7V6qv55FbAR1D983j7WL1PJ4ri3vKSX/1wzZQoKmYiIhSrn1bi0j3A/hTZQ1rwaN03JixTDkd4/yGoe9g60D8ToFcuzhXttpH0Mk/X4HvfxprpEvCZ8JBE3Gl+K2LdRRGSoElulJP4N2zqqyp3X7ciUhWrXiZ4qGsgYL/Tnzg8K30qJIpcrhkNdp3p3POaNybsXvS675Jtz9ksIIYQQQhbAvEpe/GBXlvPy+pz9EkKeI17TQhQ/XAl/Bs6vua/nL3Xdm92uU2LWui4mbG1jze+7tum+D66f/wsREWl3Vcnz8WPFzNBqHbF4aeF1SBRQgRrtEPvnfd9UQWquuBjAaqQYgRDrNUnW82zIOLMSis6w79SZyQjqELzp3Pu93SeF4yEuyiuIo5IYsyltQqyea4MqGiIhE7e358a0ds7Fb61suAib1qr7HOpRzByUx4Zm0Xp1TrdD0XPnNEiWiFHE+HOFz/WTVuaAmjlL/YPahxjAsqoevn9/v6TL3PsuVuW8GlZQ8tLjQD0TCRnJOHd8dkFxLYnJ0+NAjfPXbzhK+ojJY+NCZZP0eomIVJFtrpeu0XL3/KDljoNM7+punCWsSrDeP5XseJGto6w29PdLt61PUn2npnJdI7psuYKHWLwypXAaUPD+VLw85AjMq+QZY8wv9LWNNv7EGPPPIvLpiUdGCCGEEELmZi4lz1p7xxhzzRjzUESsLjFF+6619l8WNkJCCCGEEHJs5rZQsdZ+aIy5LSJvisgFEbkvIvestTuzWxJCnndgm/AD9cY9vxomBTZRqmzdzRutrbudVtfd9OnqRtfvu7LhEi9aa24K0U+jomSWN9SNAvExTYupULWksN7UNp22FQnTvg0t04VyXbDuKCuNlZvwmsyWI7b0QJtR39lYoHTYRJMzhr29QpuxTtthmhPLSjWz9JA4ySDdHqZt3XIQJX6gFBqmbfeeuGnbdS3xtrbprjkSMkTCFHeYToVpcXF6cKzTg0NM1yI5o49pyHT61m1Lpya99cskTfSIy78Fo2dMweaGw8WwcExr+unO3BIknmofp9PjlSnTtXGSA87Dnw/Mti360HsmLjOXzVHm5xxTzaa1YcWDa+CNwBux1QxKqVW0f7XiqYcpdpE0OQdl/jCEiZ9id+vNehgz7FU6Ou7h2K2Ps2nbEvcVGWdTrZi2RWXAfvQ+Spz1mWixEE5U8UIf6O7qj8cY81Nr7e9P0jchhBBCCJmfQx/y5vC82xCR90WED3mELCk/VpHspY5TDborQfVYXXHqT6fjVIb2ilPlOuurugzKUXPVKXcwyfXKh8oBUO3i0lhQlcYIflf1DGoZrE3ioHVvxQJlqqSElO8fwfso05WZ8Vaz7fGYRgOnKo56u7ru1DMkjsAuQ0Skv+eSMKDSeAUJKlSk/Pga9lMUPWwfjcMbsMNAMsbenjv27hN33TbUJHd9M5hHd9bcZ9TUMnJQiqAoxcpUrQlDYxg/o5RbquTFtjK5AoblWE2FR5oAEKtNlYoqnJXypIZJpCQFhXC6Shb3kb925wgrmNQSJgb9egPrTKqC2ohkB5FIjUPSRDW919P7NVXwgtVMvXR92jjjfoNhdhgrrvc4G39dfVHaNvSJy4zrjZJnqO6HLmJhD8ogtvW0jVZRk239mLcjMZXK3WKZ+ZBnjNkSkW9E5JH+xGyJSD41C2PkfF9CCCGEEPIMOUzJeygin1hr34k3asWLTa1wIdl7l7UdIWTJ+EsVr15Rw+OuKnmddlDG2m33Z6XVdipDe8WpQm0tIdZaDUpeHQqeKmGIYUOsXJniBnVssO++Y/aeuO+UUPKgssS2KxVf/ktjp2CAW2Jamx8zt3GBgldtdIptsv5gu0ytvV0AACAASURBVBIMnIMiCeNf4+PH0vNMrEdmlJuKiYUrr7ho4NO+VoPfUxlld9cpbN2dYNC81nUVKVfW3Lm19LNrdtw6zIxFgpqE8eeqH9bHURvYkCAeEWof4vhg6hzblUyymC+vznnxKXo/U5ugUEHF8gpxtRgr56+tV01TxbAStzGp2mqzOEEocTAiFhGp+ri61NQZxKshDjFdekVVdy4rY+etWlRB7e1rWTONxdzfC4bWKHE2GqVKXq1avN+CigyFE/F7ri2UvdgBSMNCZV99UHZVuftOxd1vD3cLIidkpoWKxtxdK3lrq+wBT9vcFafyEUIIIYSQU+LQmLw5s2W352hDCDmjoGj4K5pNu7niNqxrLB7i8ERE2qrgtTpOwWmvOhWovabmu+1Q1gzqWGxOK1IsPj8eBLUpV/D6e24d8W5QBaUexWrhOKoUjs3077emCuUmK2aPIvdYT8ppofxamo2aB8/FKmGu0sxSdmrqQAtlKqhPkrQ1phjQBFUL8XpQ9noaq7e7F+LfVrbddV5ddTGFq2su+HJlzX3wnbXw2TXb7r16Ky1RlhsRx1mq2GcClU9Vp0YLGcYoIRcreeWZuF65ipQ+H6eHWLksRg/W/dUZhsrAx+RV06xe916aweo/50zJq0UxcyE7ODt2SdygV3d16WMAs3JmseKJ18jS7h+4rFqUzdvfdUvEZoqEeM1RVkPMZ/dWinUPisoqsmzdem8Yfi8Qg/dYD/nvuqSx8bNjXjPk1w55/405+yWEEEIIIQtgXguVz7WyxT+KyANr7WMRV/FCRH4tIrcWND5CyBnghxpWda7tvheuaZXyFa1G3ulESp5m0yKuq6MKXnPFLWszSolJVlB+PFJlJ1LyvIK36yYMBlroHcpaHh8nErJqkXmblygzkbIXx/KVjbFqUQA+nLOJXpe1mYyL3nfI+A2Zkmn2ZVkGaCXL/MwVvFj9K6p8bh1KTG+AZVCDdjVe79GOu+6rHaf+rK66rOHVtV2/76qqe+0VLRXXcesop4WYvNhbL/egC3588ENUpS9WqOBB57NDs4zQsgxaH1+XZtmW+eYFVSwtxxZKrRUVuPxaFkuiIYM2apN593mV0abKZBk+xhBqXRbDKBIykvu9gS7d+sGB22cwcP0PI6XNq6JZ1jbOr1Ka4e23JGPE9jhRFyIxFbzTY96KF18ZYz4QkTsishX9MdoWkWuseEEIIYQQcrqcpOLF5yLyujHmgrhEi/vW2gcLGxkh5NRBNu35tvsit9rSbNqWUyh8/F07eMZ1VNlpq+daa81l0zY0Fi/2WkMsXp7BivWJKnnDg6AgDfadv9zgAL542gfUrqxv148qeJmXGNSU2GMs9r+LqagaV6iIIUUlD9U4vBKp8YJpFm+alVpvIKatGEUTqiNo2ymFHMrwqgzErbxofOStB3Vvsu/Gv7PrxtvOlD0RkTWoe6tu/Cur+rlD2Ws7RQ/Knkg4Ryh2eTyi/xxi9U/VMK+4+Xi1VHkTCcpZrprlxG2Q4VupjJL30IdXVqN7xKt9AjWuXGGNjwN1EooklEPE0I3iOMRMrQz7un2g0kGdi18Ph6liB7/EED4YxpSrvLPup6DUZTGeI1TA0PVIkNSEbip4p8i8MXkea+19a+3d+AHPGPOLk/ZLCCGEEELmZy4lT2PvprEpIjdE5HdzjYgQQgghhJyYYz/kqRHyF1JeBeOCuAoZtFAh5Dnl1cjh4bzmSHR1uraj07XNhpZtqrtlIyrfhGk6mB83O26aFokG3phWwpRqraHTbEh6gGXEwE0P9nWKViRM08Ls1U5JuDDRcTBNOsqmTTGWpFxXNl2H/jB1PCkpm5abN+O80Gak54HED9cegf1ausonKpQYQI9gdYHpQfSRJWREYxprYLyfpq2kc3GhLFjYhjwCTMHtaeX43Z5ru7MXdu48duey0u7r0iXArPiydmqG3QnT302d1m/pPVJvuvVabhRcVqLLjzdNwIg/Bz/9npUKy6dtk3JnWXIGyKdcbWI8rNOo3oC4fD4y7hNTrvk07UAdg+NSbvic8yWmYrE+jqbaMX2aT+1jerXip2TjhJ7SYfvrN4mnm7Wfod4b/WG67GEZ2bHs0ez41JlHybsvIv+HtfaXZW+y4gUhhBBCyOlz7Ic8a+2OMebGjPfvGmN+KiJfnWhkhJBT4XyUd9BtuW/0K01V8OpuHea8flkP6lNDTXIbvsi9W/fK0ajv9y1YWyjjodunv+/KbA0iJW/YT1Uxr9ZAyStRZrw9BVQgb6Xi1mObj0ISSGacjCSKGLSpNtw5B+sUlPHqF9pab8WSlbAqsVDJy3TlhrRlgfN5UH21Yku3x6BbqEK4hAeq0mApIrKnH+OOJmm0GyNdunOE6ttuhWuLkndYNpu6bGnyia7X68XEi8wVpTxpBtYlsD2pFJMmRFK7khESH0ap7FRQ/yLJc5SpcmOvsKaf0zhS+MZelXNtkChRZm0yzJS6YHWSXoMyckXSo1JeteRzz+8FqH/xfYZ7YpwlXPR1mZcuExE5YMLFqTNX4sWcVTAIIYQQQsgzYm4LlTKMMevi4vIuicjvF9k3IeTpgtJlq5EbSEuVu3ott4gQXdfYvCgmr6FlrupNDejzxsZaUixSsxCnh32sKl+IYRvsqZJ3sO/bjFTJy1U4r57pfpWo7NjEK2vjZAxe9Yj2hbJm64203yyOz5a18ca2qflxMPINNhkFq468KH1kpRIMdFM1JRSNl0NBH1ByKiWKTojjSt/UUDCv8IgEVQmxWFB2Hqt8U991fbTq4ZxbjVQRRmxnq+mUt6ZfhvsJanFlSkxhrORVvcKsyqq3P0nbxjF0iHPDtlwJ89dkHCt5aaxcrrCOM9UreS/bB6rdKNlXSjmedQ7iNNPziAnKY7qPv4+jNrmaiPPQkEJvfPwwuvVpnXL6zJt4cW/GLl+KyLW5R0QIIYQQQk7MvIkXd6y17yx6MISQ02NNQ9patSAT1Kvl8Vu5KW8tiqGqN13mpFe3Jql6FhdVz4u1wzx4oObHyKQd9kJZM2TVgjx2ykCFKMnQhPHtaNBP2iaZsojXg1GzPw/E8xWNjQeZiXN+waD6xbFgeTxiKPGFZZTxW0vj9HxWpyoluUo0bdthYNgNKLe6oYHMyihsDSFsI3+cdDvi93b74fi1iqqvemq4v7BsN1KFTySofkUlLx2ziEhV+6n5smImWeaKqEgxUzUv8ZWrW3H78ZTr7z+X6NJj+JOsv7LqbPnvW9Unek+X8CoFlU/HIumBhyWS3jQlOD4cxo1zRUW1XTXQ/l5/Lb9lRu2Z4tgxeRqPNzXxghBCCCGEnD7z1q59kG8zxqxbax+ffEiEkNNgVf8aROFQUp2ingBkM9YaIZAPfm+INUM5J6hxScxTC9moiEvq6b4uBg+ZtLF6N+z3k/59PF/meVfmnwYFD/2VKSOIifOedsiuRamyWdm1WhKtUqsX9hHJYvKgDGZZolBFY3UUSl41S40cDLOYsJIYqikWbqXxfLn/HoZQ0+3x/TDUbN2qKjdQ+aCEqUAmg3FRJqpkapNXDFUIbdbCoBHuWfVegDhOOuYYKF81r+yVtxUpUeWyDoPnXvE440m5GpeXkou35ePFWMr2xTmXqYk5dkobXGMIbJPk807789e4RDmEP96BKnd7utzWX4d/C7c2OUPMlV07pWzZa8aYK8aYXxhj/u6E4yKEEEIIISdg3uza1/MN1tqvRL3x1CePEEIIIYScEvM+5E3VjI0xPxaRt4UWKoQ8V7R0qqwWzRvlAd8Ifse0YV1NkBvN4KBcU+sRmAhj2hMlxcqmUb1ZsO4bpmv7yVJEZJhNtebTtLGxMUBSA6aO4wSIfEyYyq3W02lZO5k+XYv5taqee1062gjJFZWkD5FwPfL5O4y/0QxTvnjdVOPcfj+1dUFpqbLpvDwJILfCKCP/3L39SiVulM9flvdbKZnn9GPxU7v5FG/YVz2WpVZJOw59FLfhFkbuCqYhcV62pM20aVqcc7Vk3suWjCE+fnzqZaXCShtJsLnx/ReOl05zx8fK7U/y8mZlhGl4tEWCSrhfUeJu58Dt87Dnlo84TXumOdJDntqmvCUiPxeRi7rtH8t2FVe79uqiBkgIIYQQQo7PkR7yoqnYD4wx74vIxFr766c6MkLIM6GbBYnHioW3SMkC2EM5MzXurQfVKS8hZTMT4bKkBKhjsEzp76t1iip4g1jJ67t9kajgC9VnSl6szgWT4nQs/viRrQsUtuqgn+zjjZRHqdFxfOyxmjjn5sj5OESCqlgIftdrWY/U0WbbvR6ob0W97sYAZdWb70ZqTZ5UkicYjEsUvWCci8+7XD1zr23Sxt8rmSpUK1HAgnrm9h1lySFlRs35PkdR8vx743Sss46F9XFBETuCAzHGgRfR2MbZBFh+3coSYCbaxmZqbJ5UEffjjzcpJtbk+1WzC5WXtetFZex2em4Q2/prAdNjGh6fbeZJvPjP4rzyCCGEEELIGeXYMXnqk/fhrH2MMT+21v5h3kERQp4dq/pVryyOCN/0g4KXKnm+bFRZ3FVWSszH9UUxc0GNU5sVWKZksXiDXlDVUFDeK3fjtBg9lLFaibqYq2aTkvpRUNig7sEeJVckY/KyZd5IWWPxEBMYq5wT3yaND8S+9UYck+eUvGbTnXuj4dRGlJvD5Y/tUipmWgxbqtrE6p81qQqXf67xvrltSMFex9fgkgIwrK5kylsZuSo3Lf6trJ983zK7ElDzcYjl4xiVxLSVxRuKiBfYyuLv0CZYnui+sdI3Kb/u40zROwq5Qh+PGf2qO4r0YZOiKt2TQTjQIw0h/U5F751jjIGcHnNZqByB959Sv4QQQggh5AgcquQZY74WkY8Rg2eMuSwin85qIjOybwkhZ5PSeCWYymYKXsUbriJmKEhIXs1CeTC8pwoCsm9FQnwe9oUcNM5i52xJua6qscl7E4s21WRsyTjz2DyvakX9o2QYxl8pl3iS+Lox1D/XpubjD1VlVJPkahSP6EuUjbL0RKhcUck3lDqD+bTPalb34EYdWbZRqbUR1LjS7j2xsuOvwyRX9Ipt86zNaepSuUqn2Zv4VzFFpSvbdhwV67DMWZGg4HmV7wj9Gq8ulg9mWtatiIjV+zaIdUWzZcTv5fGToxlKHq4zltNKEuK+EBEZauzdvv769fS9J7r+XRR+yti755Oj3M83ROR2tP5QXO3a6pSfiojceSqjJYQQQgghR+JQJc9aeydb/8oYc+2QZqxtS8hzRpk6EPzFUp884JWwSMnz6lsWrOU916K4tGqWaQv1ChmyYRlUrbwMmMnGlHu7xWMIGbKpeharZnmsXH68vM+4P5RNqzX6yfn5ZaRi5j5/hbJs0fGg4CEGEopqs+nG3Wq69YN++BwOBurrNyVzdUat+4KiN2vfeahkKhbUrTyDNh3T0fsfTVEXKyUK3GhiSvfN28RZwng5TSWdpbjl/Zb51+Xxh6G/dN/Y07IYe5eOCSrdfqTOIeZuT38ddlUI/jYNEyXPMXPF5GnyRYIxZj16/8FJBkUIIYQQQk7GXBUvjDG/sNb+Ltv8mjHmgoicE5H71tp/OfHoCCFPnWGmFpSFGeXq2LSlSKxEqTpm1A9uhqyBfaGoITMWClktyjSFshaqcKTfVcsyaRHjNxoOk6UfdxyXlql9lZouK0GFy0Esnhm4Y4/UL6/WaCXnFyuXOFebxSFWS86rWk1j8qBsIjavpYpesx4kmCcHbgmvM1z+PIu0NiO11ajCVpaZGfbRMWbKkUmF3CNRNpSp3nclbcK+aUxhPu5Zx8mzj8NxQyMNhSyoJLnyFmfkTssOPsqY8n3wmcXqIl7jeiPerqcq3b7eB9tRwZY/qqrHTNnlhbVrCSGEEEKWkHktVKY+90e1awkhhBBCyCnB2rWEvOD0YYVRMp3k7TF8iaTyZTpdm5oT5wkGKGEmUiwRlk/bVrNlWb+FMmqY/sytSSSURBsNBslx4unavFyaNzrOpomtjbIEMMUHI2U9R5xftQ4j6Gi6Vl+bzJ8kT8Bw40yta9AGCRiwUmk2QhtYaGDKEIH3wWpDx2SK39kx9YelT7yJ9jW+n3Qd+8LwOE78GGcJBLOsRkA+zZlbhcwyUkZix1HahCnWdBnej6dey/tFm4H+MvWjBIYh7FC0m2rWthqNqZ4lT4TPozhNm49/ME7HsKu/Yv+uv3a0QnmxYO1aQgghhJAlZJ6YvP8sIu+UvWGMWbfWPj7ZkAghz5KdTCmJlZehqj8DlSHquhzqcqzywTgqzYXyXyAoX2qTMQzqHQyHTbWW7GumGBDH5BYtuRVJHPHvj60K3mioyRR60olqVnHvwb4l9KvKnj/PyADav86UvKGzUqno+cXJFEguwbEnmfVMbOvik0yQgJEperBWaTZCm1bDjbdWcf1C6YGiBFUuUW4zlUnzOUKiQUkJMahKUA6rlVSWK0u8mGY1EhIXYouT9L38uGXJILmFSq64lSWboA2O1xunx42p6O9FPVPUhtk1jhxtfOmwGq4trpsum+GjK4xzmmIYXycoeDAy3lEhm3YoLzYLq11rjOmKyJvGmA23m/39AsZHCCGEEELmYN7sWo/6423q6jfiLFR+JSJ8yCPkOcKrEONIHciUvKoGGSEWrN93sgFi3USmGw17FaukhFgVShcUvUoka8jsEmXizYtTJS82UIYql+PjBKMxQyUrGCfnZc1iA+hMVfSKYb+nx68XziOPXfRjQdm2EgkMil1uAI0x1yNpqdVw25oqHY21vzw2LFabnqjqgzjNVe1uteo2tKKPBWpf3ceLuX0a1XL1KX8dkyt4g0h9mqaohbi1ok0Jzg3nWox/C22wLVfh9vS4ByUxbBpiKS2cOxS2bIzxOvZtZModlmXqImIXczUTsX696DpRuSNlzOuT1xWRj8UlYzwSRLeKbIjIF8KKF4QQQgghp8q8St6vROSGlji7Epc+M8ZcFlfflhDyHPFYlYCNYZAfVlSKQGxeQyWSkUoLw4FrNOgFJQ8qli8HBpUuK2EmIl4dwz5ST+PRSkuUKbm5MuLsxqOxdhXe9+XR6ukYxl79i9VFjWHT2MHctLhSYracK3nIqh32D9LzyRRKty0N7IJyGPcPBa9ad+dRqxX7cdsjJU8lorbG5kH52tfPGWpXL1Kq8szLHahCJerQqzqEXM1qVNJ4tThurZZljebxdJNMgYu3jTOVTEoUtlydzJU1jLFWoiiOsuuxq0vErHZL2rR0WTXly5j8euQZs2mpNT2PTLGDurg9TsdGyDTm9cm7pxm3Im561mOtvXuyIRFCCCGEkJMyr5IXf3/4whjz0yzR4oKI/Nf5h0UIedZAHXg5CqkbZfFQuaA2VulkOAgZs8O+yyidjDWDtdoUkZBNmsSl5SW8VO2r1V0JMShvsQJXreFYKJeWSjrofxKpc1DLanVkueqJ9AeFfaGkoU2lOkxPvlquopWNYTyEoqexefViabSg8mWxedF1qlT0Oviyb8hGTo9XjSSktgbQrajU1qm7fh9qwF2e7Xlcvs1Uvlzpaup6O5ISinFpOu4jjGGWrxyoHEHtm0au8vl4RH2/FZ1HnhmLZT6mgvoYb/NqXRorKSKyq7+D3+mSih2Zl3mVPGOM2TLGfCQi2yLyn4wx/8EYs64VL36+qAESQgghhJDjM5eSZ629Y4y5Ii7p4ntxvnmfirNWscKyZoQ8d0At2I+UvP6oqIqJFBUkxOaJiAwOXBwaMkuDKtdI2rh+UC4Bqpkqbk0X7dRotUUkKGIiIet1NNCsUZUba1lVjkqUrehFuBqyd+FNp20iJdKH1/ms3XH6Rt5pybY8htCWZMpCwas13HWB0tk3vfLjScgYxtJnAsPvLzpus+nOdbXjxr+isXnILN0tibODGjePclRokylWZfxI/wOtqDhaps7l3n0grwohEhTnPDZuMCnfHrfHttaUKhzxcaapitMydUVCXB0UVCyRvctKFORpMLeFiiZbIOHisYi8bozpqo8eIYQQQgg5RRbhk9cVp+R1dX1bRO5ba//lpH0TQgghhJD5OPFDXlwBwxizJSKfiMjfLKJvQsizZyearn1J7VRgijwap9O3SLwY9MN0Z2/fTTdi2rbR7oiISL3ulvGUIsqZAW8a3HDTtfWWlgcbxdOpqZUJrE5yK5LU4kSSY2PatlZ3c2iTpCxbOm+Wl1rLjzNrX28EXcU0a2ThYtw0bUNLoNWbLkGl3ujpOYcx5WXNkHiB6WF8DvHQYKey0nH7djvuOnWeuPP7Vl1vTjOo/99wr+nylZJLC+NhJETkU7t5aTGR6abLZYkQZabNhzHMplbRb25Lsx2bFTN5gpwC8yZelGKtfWCt/Vtx07eEEEIIIeSUeFpq2+dPqV9CyFMmLot0XtWertpudNRlttnQpbq2DvpB/jvYcwpeb29PRERaq2siIlJvOSUvNkWGOXBQvlRha7qECzspUdpsmmRQgdGxWqtAKBwOiuqftxrxSl690D9MlWWKKoe2ZWbIBfPmbL0S2a9UapqIookRzZ67brCgwTJuF/pRc+RqmkACc2eRkJTRaLi26yqBdZt6rvtnT1o6UvKBfjy5GbNIsbxYbmwMmpFqB4uXViZ59DJD5cgj3Jd9ozpHzjoLVfIi7j+lfgkhhBBCyBE4VMkzxvyDtfafjtkvv9+QpQdKwioUhcgjNzdH9UXPVYV4pMuzbpvwnSp56yprtFH0Hkqe2nL0IyWvt+8UqINdKHm7IiJSVzuUWMmrZgpete7i0qreJLmkDFhWCq06cMcbZ2XIKtWebzMtzi5X50RCWbMcxABWasU/m7nNCo5XOG50PrCUwViaHafkIZZxEreFUqhjwHjrjdQUud8PbapV9zqPzdtccZ/VD9RDZSdUpHuuyM2Yj0X8H+qM/w4SchKOouRdmqPfjTnaEEIIIYSQBXGUmLyrxpi3jtnvhoj8xznGQ8iZIi7VtKEizIYKUau67Gj6X7sedkb9+Hrm4AqRaKhpePuDICk81kCfP6oA9W8jOXWglnR1TCiN1dKs1IZKlo1GkFP6PaeoHey5Ru1dp+Q12k7Ji0t7wSgZip3JFL1ao0T9w3v7T0REZKCxbGPNUs3NhUXS+LaYciVP4wAzc+WQ1ZvGEYoENQ5tba7oeSUuVvKaybI5cNerv++u12gQJDaUWoOa2Gi6c2x3tK0aH+9HTtYjjZdETF6rpVm2Gpt3vuXG+t+fUyWPEHI4R3nIuy8i7x+jTyMi1+cbDiGEEEIIWQRHecj7xFr74XE6Nca8Med4CDkTvKYi0LmonnxXU/LWW27Z0SrlTVXwapFqh8zGaiVdz8tdDUchIGiggXtbGv/21z2n/uzo+veqpv23clHqqQJPtc6BG1OrnsbktZpRWbOB+7OC2Lz9Jy42r95y3ne1RtPvi9fQ3KCSVetuX2Sg1lqroU2zo0un8kHRG/b2RURkpDF6sWJYqbn3ELeXZ79WI9XPx8KNUinVBrO9Qv9QAuHZNxm7tj4zt1L0yavreUDds9pmcOCu1yhSH6FE+pg8HW9nbcUtO8616vHj0CYv7xbKnGmWrd7HP4qybKEen6S8GSHk7HBoTJ619lfH7dRa+8v5hkMIIYQQQhbB07JQIYQQQgghpwhLjxEioZzSD3Um8XzbzVed64TvQev6uqNeKc0GDGmL9ZCCKa7oPum+WI9nbxHoj4D5wcBNhfbUFmO/59Z/0gtTvNv7btv3OuX2wM1KLtyaBdN2/59OGbdraQJGqxGuU6Mx1KWbOqxjulaTBVC+S0SkoQbJ45ZrY63rF1OYtaabjqxE5c/qOnWLZa25LSIig303ZYlp21rjwLfBtPCo704A06m5ObJI+OwwtZuTl0aLx1dRQ2YkgeA8YLsSn0cFCSR6jth3pAkYGKse1LVBsgdsV3Qat6UJGDWdlo7BPYe2raY7zpq6/27UQ9LMgd43yCFq6ud+1q1+CCHlUMkjhBBCCFlCqOSRFxokWPzAxfDL+RX3vWdTHY7XVsKvSFsVPJjLBoWkmFSBbViiTa2OPjQQvxrZcOi+wX4jLVWFgvVQ+EREej2nSO3uOuXof33ilKRHanT7/V6InP+zJk2cJHFjqqVKZKHSqI/03NyYqnqu9bpLkICiJyLSUIPkWlOTMpBM0dRMD1X2KrUoWaOpCl5rTZduvdHeERGRwYE7zqgfVK1R36l6UMmgtI1HmigxCqpdXY8Na5ZJZl8CYlUOSiGSMSa69MkPum9soVKp1vV8nJJX1SST8dCNEYrkLIYtt2+j6fqKVeXxOM2awL1Y13sQ6ut69NkNM8UOXfyJNiuEPJdQySOEEEIIWUKo5JEXkr92wpH8YMWpGy+tuu8759bUMHbNKSMwkBUpxtFhPVfpRETq9Vqyrabr1UzBiwvWm0oaO5XbfIC4VNZo6FSm/oHalew69Wl/T9f3g0K1t+de/2+7rs22qn2P9l1/j2DZEoWiPVGRBzFZsNaAwrOrZs6P94MaBCNoDDtXMxutOCYPtipOxYLFSIhx0+2Njm9TVQWv3u66Pjrn3LXQ8mmjvi57u77NqK+2JAM1TvaKXl+3h/g3vIaVyeDAKWreFsUUlVsYGsO0GZ/RZKQSmElj6dym1C7GqKLXGg91zCGmEPhYQu1/qHF79QaUvPC9faQKZG7qjCHgc+pEJt4D/WAzEVBe0/VvysMUCSFnlLmVPGPMj40x/2CM+UW07bIxZn0xQyOEEEIIIfMyl5JnjLkoIh+LyCfxdmvtXX3Q+8Za+4cFjI+QhfFXQUCS/2nNqRcvrzs549y6U0K660456nTcer0eK22pcoc4KCzjWDOvrKhcAsUOyzI1KBzIJPtW87ZRbJhXdLQEFkxzURILSp97z22DSfHerlOBdp/oUuP69qLSWLuabtkb6HLoJB2UZatkWcQiIsORlmzTbOB63ck/rX3tX7NtRcI1Q6ZqHLPmG1WU0gAAIABJREFU+tX4x1pUCk1LnVV9bJ77XllXRc+OkeEalLDxAEbJ+8k64t+g6IkE1Q/Zuv1dZO86ZXAycedlopg8KHn11koy/rEqeVbb1BqtQht/Po00yxaqndtW0fGnMYVDVfuabdcXSpiJhNhNKHlB0dM+kSUcfdVvpJff42P1qOQR8lwx73TtO9ba10Wcehe/oQ96PxWRP5xwbIQQQgghZE7mfci7F70uK3yzOWe/hCwcxJH9cCXITVDwXta6Zd2uU0JW15zS0mxp5mMcZ6fqXLPl2jTabt9GE5mVoWRVrtj5klhKXrg+3ub7QEmsGuL70ng1PYAbm8a2TcZjmQb6R7H7Qc+pQv19p2r19nX9IKhag/4wW7q2yOrFcpinZUrI9EQ5LVyL0TCMsb/vlKiDhsuI9Yqej0dE7GLkL6eqnqniemjpM1X2Ku1zhbFMxqqojYfJ+mTkzrVM9RvuOwWv9+Q7t3z8Z3cNDnZ1jEECg0JXa6WZsvh8J3rcWJGsIzu4vSEiIo2O+7MJRbIaZxSr2tffe+TGhgxiVfbaq+647U7ov99Py7JNsmxtrFciGbaWBfDoLlI17kU3UmxZ8oyQs8+8MXnxr3fJfJO8Pme/hBBCCCFkAcyr5BljzC+stb+T6IHPGPMTEXlffwg5E/wvmpj58lr4TnN+wykem5tOIVnvumVnVdU5zQBNM0E17qqJWCpVkKDgRYpIrppZXWIdits4Ut5yFQ6qVk234/200kKqGIbYv1ph3zzeDTFfqOwwGmg8X1TpwWeJ+mxRVe405q9/oJmovWCkhjhA7xFXQVxdmmksEtRR8SqfKobqUVetQ+ELKinOA9UxEO/m+/z/2Xu3JTuuMzt3zHWuE1AAKfbRbQlShK8cu5ukXqBF9QO0SekFWmT7AUxY+wGsjW6/ACU/wJZI9b62CfYDWGyprxzhCBNq2S1LpkigANRhnee+mP+Yp5WratWqAlBYNb6Iirkyc2bmzFyJQtaY/z/+NjNyU3xcx9Q916p+5XleV8qunY0tE3cYfPcGR18CAIZPfhuu2ZS9SZa9W3vesRoH13OsXAay7OCdV0K7+3vWxyqDZIok4w+7Tz8HAIwPw5gYS7h9Mxxr72bKQqZiyu8h+i1aPOXU2nmmJreqP9mn8+ZsWyHEy8FaL3ne+586577nnHsIwFvLKdp3vfd/f2kjFEIIIYQQ52Ztnzzv/Y+ccz8B8CaAOwAeAPjUe//4sgYnhBBCCCHWY10LlT/13v+jvdB9Yj9CXEl+34yOv7Kfpspu3w7Tsvu3w/Ta3k0zot0J07b9rdAyoQHIDXuZHFCGtOaWF7M47WglqmhtMimna+d5skXytgjnsanXaYflwVjkPk27Mtmj3tbpNiRC2L/2OJXbLQ2IWVqMU5hAmhp1VUkvjnXK8mCTNF3La+e1LdjF5IbA3BandLvF2EhuVjw6fGTDtClvm7JkEoWfh3vc230t7sMp0FYnXGOrnrbtZ/aeNss7HQULlS4TIuL0algeHz3Mrnlo4w7HZ1IFp1k5TZxPmXf6Zuq8bdO1268WQ+pkY+LUc9fsYk7sPLR72b75zwCAvVtpHz5rTJahQfZ0Gr6XsVndzLJHpZ6uJcyrUbKFEC8X6yZe/OhSRyGEEEIIIS6Vdadrv+6c+wGAzwD8xHv/5BLHJMSlwNJlr+wF9eSV20mVu/VKUFhu3g7Kx/ZeUEr62yFwvWdKHgvPA0kli4HsUcVi4kJSs0YnIXFgPCzNiRn8TiWstlYBkro1n5XJG60220zJs3N3KkWPCRKtLIkiWrB0yr5RrYuqWp7YwfJiVP+sfFel7OWKZJ0IUduhuIbSXjVLjwFgNjWV1GxEaE6c7FFMLZ0mKxgmVnTNOqUzCLYl7U56Jmravd1i3NHOxfbJEztmo8NyXyZKmPKWjI5zK5h+MZbT6O28VoyBauXIEjC29sIx9m6lY/G5PHwc1D4qeeNJqeTl0DGItirmp4xZw3MqhLj6rPuS94b3/pdALGX2NQAPlHAhhBBCCHE1WDe79pfZ508AwDl30ypd3AZwX2XNxIvm1a2gRrxyMyhX+7eSvcTeflBatm+Edms3tP3toLgwPi23IKHqNh2xtFRQSqjaTYYpbowK3mRMdYl2FsvH66O6R1XMymfRHqU1XRjTfMYSWKW9SrRSyRS3SWWqHGMMLf6NfYuyadWAqby1Woxxs2OcUp6NsYs8bmFsbJYidczf3NQ6lgXjcj4mKnbT0TE3hCbG6qXvYzq0+LqdYHDcNQPizuCm3YP0bFBtixYtpv4xVpEmzFTpwj7tYgyMzaOC16TkLdi5rABVRMcYQ7uXLKe2fTMpeVSPaWTN+0Ylj6Jyt52+u16nuIwYo9cg+gkhXgLWjclr4lsA/hrA3wB47xKPK4QQQgghzsm62bV/6b3/OzM//msA3wHwMwAfeO//4jIHKMS6vLIT/oa5cSOoHjt7W3Hb1m5QbgZVDF53UMbi5bFgzCSlajatYvGo2gGZCXIsXwZrl8fiJeNhKlVl6TNnskqnk7J4k3Jn2amsNh/Vv/R3HNUyGhDT1JltitFb3IeKUceyXGf98j41mRUvKHnzRbUvKoOuU447Kkp2T7LM5aTYUZEaF32o/jHrFgCmwxC/1z2mkhcUPJZCa/eSkpcbFofT8Xw0eS7jKvNtMWau2odlzVyWuYwpFeGgRHrLps2zasnM+k6OrLTaUTBkZnYt71tvK11Hf8cycnuPi+FSweOj0s7+1E+xeKHz0CS8p8ur5QkhrjDrxuR95Jz7zD5/AOBr8scTQgghhLg6rPuSdwDgPSVaiKvIa6ZM3DQlb28vqE1bO0nJ6/YYU0YlrIxli2pWptZEb7tZ5XU3K8tHhQOU/m/Os2+l1mVZqdzGclST6Kk3L47ZzqQXfk4ZrCiXM1pUAq2sWH8Q7gHLtXVM4cuzd12lCNKXL6mAA9s3+doxxq/OnI3HylS/Tm9Q7BPj9aq4u1zJ42fG3tWqaLtDhSxT8sYhbpLxeyxJxjJkjKFrGncsfWZqXGNgZSwrZ5nLpgyyRFq7u8WDZ9fBWD9f7NuyvrkP48RKrI2srNnJwW+s/a1tP14Yeyp5VwwR3U74sNULx59mNcuYcXs0Du1Du4WfL9ouCiFeAtaNyfuBXvCEEEIIIa4u62bX/u1p251zX1V2rXhR/LGJMjd3gzKytWNVCPpJbcrVKiCLhzNVzrVMJcr82moFz89LeSNXz6h8eRZ4NzVuMmYVAvN6GyeFKm4zc7LRiIrecs+4s8i7cr9uN4yt3w///AcDy9DsMfs2/VqISmGlAvb6VAFDS2UPaKoMUle16GR9S0WQKl/tw5erdT4FOAKo4vWQVfQYJyUveRqGdV3LvJ1b28qUvLZlDMebR1/EmOkb2vzZiJ+pttoxmpRCUvv7eV8+TzwPAIyPg5I3fBqqbIwOw3JeCSSMMVc8y4zuTic889uD0E5MtTsepvNOTNV7PArtf09DEEK8hFxmdm3OvWd0XCGEEEIIsQJnvuQ55/6HVbfg8recc7NTfuYA3n6moxZCCCGEEKeyynTtXQAPsuWHAH7qvf/Osh2ccz+56MCEWJdXzAR5ZydM/Q22OH2YpmijsbHZn7gqOj1NCWaJETYVNp+X1ia+svQAsuLwZop8fBim1Y7YHoUpupOTNL12MrRC8jZ9NpqYea1Nq83j+dK1zsuZywgvJ7fHoOltx9p+11oLwN/qh/vT66Wder2wrttlG7b1+l1ryySO/HPbpnbbNjXO+19O1/aa225zYkz4XE61x+9jTmub8J22sqlMTmuy5fQtTZY7/WQ94u1z0xQrkBsfZ88GDZhtinViJdfw9Mti33KKtzmbgWMcHSXDguFTm649OrTxh3HzXvK+5dPaYzPp5rRt177L7e3w3dEU+WS0OF17XM6ACyFeUs58yfPe/7Ra/oVz7ntn7Hb3QqMSQgghhBAXYt3EiwVPPOfcDe/9E9v+y8W9hHg+3NwyZWprMZGAMEC9Ni2uyZMcYqB8laQxGQX15uTwOPalcnf4NKgpjx8HRe/gSej75Cjse5gFvR+OSgPa6RLbipZb/My2R5XOLrmddebndhWkQTsXWml05nmSgx0/GjGbGleZL+e2Lq12qdjx/nM5V/Lq8mtx32jVEuxE2r2U2EHjZX43sbyZJVqwpVoHALPpcbmNfc1Spbu1G/v2tvbCuE3RY+k2Knh1C6RSZ84SOWbRsiW0w8Ngwjw6ehL3iSoibWTsGr1dz/Aw9T15Gj6PrXQen0Eack/GpVE3kCl59n3yu6Ma22otJvDwmRvJMkWIjWCtxAvn3F81rP66c+7fOOf+yjn35xcclxBCCCGEuADrmiF/o17hvf8FgF8AoezZRQYlxDrcNGFiu0+rkLK8Vq5yzKalkS6VkVqFymPBYjmzMZW7YHT75GFQWR5+8TT2ffgwqCgPH4e+jw7Dvo9PLObJ4u1GDeWiTGiBhVBh22LnBl2a2OYF5a2cVZdxdaHdNhWTaiaQLFM6HapmpaITlZ5e2oexd7SfYUt1rmOKW26GvMwomfeyKJvmagWP1iPb1rJ8WoqPq8uO0RZlMgzfB2PZRkfp+5iMGJNHlc8sbMxEuHN8mI3/cXHuqCZ2rYSbjSUfx4Lxsy0zBm9yEsb29IvPY5/RcTh330rrbd/cL87XydTLwe6ejaFbjH8yGhXHoqqcf2Z8KBU9qrO1ggukcmbTBr9nIcTLx7oWKkt/BTjnvgrg22seVwghhBBCXAIrKXnOuT8D8BaA7wL4M1v3flNXAJ8BeOeyBijEqvyJCR+DXlnqax6NiCexL9fVxrpUoWZVJi2QYpwOD4LS8+VvgzHtb/53WP7N5yex728Pwv6Ph2V8XR07t58lcO6YCrdnMYU7W0EN2tkuVbmtraQgMbuV6hvbwVa4jv52KuXW3zJlzdSgWHLtVLNiZsoyQ7l533ZWqoxKV62ARVW0yczZ7j8VsKSalcfI+1AlY/wbqhg9xrzl66ZjqrGm6B0f2enT3611CTcqasvapmtk/Ob4OKiJh4++AAA8/t0XcZ+jJ0E9HFi5ParJO7duA0j3DwD627vFOUemPDJGb3gU7sHR06O4D8vjUbU+ObZM7+NwXw6Pw/qnWVwok71PFJMnxEaw0kteNhX7t865ewDm3vvvP9ORCSGEEEKItVknJu8/AFjqkSfEi+KGhYW126XXHWOScpJSVyt5ZUmuPI7v6aOQIfn5r4Ma87/+Z1j+n/8nKCS/fZJlNtrHgcXMvbpjyt1OUKFuWLu3mxSw3d1w7p3doIRt7QTVhqpcb8BYtxT/RlWuVWWuxri4QVKDGCOX1DhTxMzvrY6PA5IqxzJdrvIRzH3f4j5V7FrH1Lg8vpGkjOWqjBzPY8t56TJXravLjdVlzvJr5TXy+x8eWRzfcVJhZ/adR0XPYhWTh99iHGL9PZCJKW1Hj4Pae5xlYI+GLGcWlnsWC8hj5Ooiv7OZqZOMyWMs3tODoBgePExKHsvjsaTesZnffXkQ7hPV5v99mM7zhYX0fS4lT4iN4Nwxed77x977HznnvlYnWNi6r17W4IQQQgghxHqsa6HypwBeB/Cf7DOA6I93K18nhBBCCCGeP+taqLzpvf9PAH5ab7CKGH8J4B8vNDIhzgmnRpncEC0iJpzWy6ZTR2G6azotg9M5NcdpWpYlA4AvfvsIAPDP/ytMq/3z52HbkTnH3t5OfzNxWvaVm2GabX8/TGHu3QhTmDt7YQpzaydNpzIAf2CWGt1BaUHCKbt8OjVNQ5bJEywL1s6sR5ZN09ZG0Pm0Kq1LWkyssCnEmMjAJIdsitTXNdYajhv7cuo1TkOOimPEBI/ceDhO5YZxc7p2ZuXAJkPaiaSyZjx+muKdWh+zIDlJfSfjqZ2H46atTuiT7GPS9XCaPNrD2Phn0XZnsrDPzl6r2IfP6dFBCANgok9+zXwuR7bt8ZeWBPS7MF37fz5P08FHJ3Zus0Wh8faXR2H5n639LOWnCCE2jHUtVIQQQgghxBXm0syQK74J4O/WPLYQlwIDzpNql5UQe2pF4M2NeHubxrEs5h7+aRwdJoXnd58H24onh+F4u2Z18gevBtXsVuaHcut2UONu3Aomtjs3dgAk49seVbp+suFI5sGWuNCtymlVSQ9AnrAwL/pQeWtlxsMxWaJKUGD5L24vDKCtL0t7Uf3jPjT5nY7TfeJ4a6uRZKGSjYk2KFHJGxfreaxcvawVwVRmrryeWWGhUppf875QJe0NksqbmzUD2T2tyoDNs/JvmE4bt3FsPG+u5NXfJ5W84bEpkpnlTzLvRrHPyZGVyzsI7W++TPv87jCcm2XyhnaJTK74dYMRtxBis1hXyfuxc+4/O+f+r3ylc+7PnXM/A/DxxYcmhBBCCCHWZS0lz+Lu/hbAT51zX8s2HQD4jvf+7y9ldEKsgQl4mExMyTgJ6sbxcSr59MWXQS15aoawW/2ghBwehj49qymWq39UBl97NShut26FGLpbr94AAOzduhH77uyHElUsR9Ub0NR3eYmvtK4skUVq1QtIChVjzmIsm6lZrexYVOFSebZR0XJ9rprVqiLsGIx3o0XIZJzuLRUq2rjQVJjXlx8/G1wYA8vL2X2hIpkbNEfj5Ooetk75dTZvVeXT7LhUUvtbKTaScXoLhthLYg3DJsYqhr4pDnR61q7xvg9Pwj08fGLxdo+TOkrFucvydVu0WQnbh7ad8aEA8DsroTdaYofy2il/4o/suI9V3kyIl5p1p2vhvb8P4Bv2kncHwAPLrn1mOOf2Abxri98E8APv/c9tG9f/BMBtAO957+9m+74P4IGN9T73W3KelfsKIYQQQlxF1n7JMz+8twEcWKYtnHN/DuBT7/2TSxndIve89+/Zue4A+Afn3Ne89wcA9gHcA/ABwgtarJ/rnPsQ5Qvhx1hSX/c8fcXVJMY4mQp3YkrekydJbfo/j8K6h0ehz1Y3qE9jU/9u3QhKz95eMrz9ymuhtNTt14JKt/+VVwAAu1aGarB3M/btbYUYvJidWseRNZgJz6v4unmVucpYMypv+Wduq2PPqpOGvlVmaVSuTIVqt5PEE9W4yuQ3HoMKXm7ca8odr6NtfZOKlmellhnDyVS4LCmWlzWL+8S4w1IZ5H3kPQHSPaxjGNlnPExZqeOT8Hm6TNEz8qzkqZkTU8Ebjyxm0eJBeW/zY3iL2xtbNu/RYRlf98VBGv/JOPTdHYTvhs8njb9PTMkbZqGBs+oRsKp52G4Wiot9juyxfLzoKy2EeIlY1yfvzwDcB/AqwssVAMCmab/5LAyR7aXus+xcDxBe5qjgHXjvHYBb3vuv23byeqXGPXDOvbXkVOfpK4QQQghxJVlXyfuO9/4bAOCc+1a+wXv/ifnk/dMFx1ZDpe5vqvWvVOc/yJftBa1YZ8vfRnhRXauvuDowtojJjxSVqOSxPTxOsgQ9w44noXPf/iX0LObp1q2gwP3hv0iP16t/9PsAgP3fC+3Ora+EfXaCgperTTG+zjJJUzxd6ek2y5S8uK7KDqVaxzg4KkxAUtJY5irG6DGrM1eOqrgxZm/WZd/yDNDoH2jr2m167FlcXMzmTftQ/WMGMZU9KnAxvg+Znx9jFas+Hbun0acPKdOX97vVrvz/YtmzLHYxlj4zZW1KVXS4MCZ+d9Pu0I7THJs3yzJqncUm8l7WmbO851T08r5Umo+OQvvYsrcfH6e+x6bkUdEbW8Zsx5S8A5Peno6TfHdoj9aE/x6sNTEQ7TJZOIzB9vmVFDwhNoJ1s2s/zT43hebeXvO4SzF17Y1q9evIMnmdc+865952zt1zzr1uq/exyJcI8XY15+krhBBCCHFlWVfJy1/sGv4ePNNHb72TZtOolmhx3xJAYJ85RfuRc+4z59wbON8L56W/nIpnj4XToVP9yTKL6ollGU7SY8tqBre2woc/uBX+KfzxH4ZYun/5jT8I6+98Ne5z6w//JQBg53bY1tsJfxNQSWL2KpDH0VENKis6TKvqDPnnqNjFeLuJrQ/L42FencGyais/Np8M1VATN8VMVlPjTKXr9pOq1RsEtazLqhvdMoYuKnkNGbkLPnmdUtELx2OfQdXHPAJNVeMyALS6g2LfdqdvfdNxAcBnSh6/h+k4ZK7OrJ2cPLV7kPvwhc8zy4iOHn7TMkYyj43kvZz2TCHslXGJvE/IVNPZrHxgl1UKAQAmeQ9NgT6ZlH0PTcH7XRrSog/ekixbIcTmsq6S55xzf2WffbbyT51z/xnAf7nwyE4/+T6Ad7z3MSGiisEDzM4FwMNzHPo8fYUQQgghrizr+uT91Dn3PefcQwDeWqpg717UJ88yWuMUqff+61WXewDeyfrfAfAP3vtbWZ8HAL5ubT0N+4qtrzk4R18hhBBCiCvLRXzyfuSc+wmAN2E+eQj2KY8vOqhcoasxD7t7TLCw2LsDAHerrvsAPvPe33fO1dOw+wA+bDjvyn3F1YGB5JwRY7Up2qFwujafDtsxP4lbu2Fq7o/+IEzT/vFXXwOQpmlf/eq/ivvc+L3g+93fDckYTACI03ejw9g3Bvrb9OxkGMp/jY+f2vJx0QLJuoPTsVMmVdg0IW06ZvmUX1VOq7YpyUtx1VOsvGFcpn1JNy+1FpMnmAjRbOYcj4k0xRr7MlGisj4Jx7Up1x7NokvLmVTWLN/H+sRpW+6TpnQBYD5LFiScrm2PaY9yVB4/H5NNHTMRZjI6sX1sinfE6fJsbpTXbveh3WFpt3IKObdQcU2BLgDa9p11s8yIljNbINv9xLIo+Kwf2WNwoFJlQoiMdadrAQDe+8fe+0+89z+y9sIveKfhnHsbwM8BPHTO7dsL3ps2Vbuf9dsHcMd7/0NbdT9LxIBtu299X6+2Le0rhBBCCPGycKaSZ7F3n3rv/9GWv4VggrwKBwB+zH0vgk3JNilqVP1+aCrfAUIWbq4Gfg/A9+0Y30Sp+n0X4QXxvRX6ipcAFoefmtoxNSWvnalat/fC3zevvRoUpN//o/A3Am1Sbv7+vwAA7H3lT+I+g5thW6cfVD9EQ+Og4PlZUtWmI1PlLLB/dBT8waOSVyVXhHWVOfGkKlVGxS1Th/g5Jk1YgkSnSpTI+yR7l1LJi4bEpmSF/XvN+7JEmSlheSk2qnt18kQsa5YlUXSo4LGtkiiS6pj+HuVxW12qf1t2/Eo1yxIvWh2aKlfqX6dUEgGgPQzf56QVvismXsTkCpaDa7CymVeJFq1W+Xc0n81w3NLih0IzL7mTmRZ3KpsginwmVkd7FJUhE0LkrDJd+9cA/l8AfFHbR3j5+WCFffcB/N/Ouf/qvf+P6w0xYGrdkgmO6I9Xe+jl2/iy9lG17e6qfYUQQgghXhbOfMnz3r9ZrXqAUPbrp6uexDn3b847MCFWhSrGhOavtoIxeaTXTX8j7N8IatDtV7YBADdfCUpeLFG2G5ZzhWfBbHdKRSfEd02GKSaPn6ncpVg8s/CYlubFOYynS3F1QcZJalqTasYyYKW5cJPJr1sWDOZaRb9in1a5bZlKFz53y33bpS1KruS14rp+0ce16jJw6btMxy/bpP61imPm+9NWZdbbKfrkFirsMz4O0Se0vWHpM8ZOjk5O4j4pfrIsK0dlL9r5TNP3PY19i0uNcZS58tzvhM+MyRubOl2XLhNCiJxzJ154738B4Ber9nfO/RjAz857HiGEEEIIsT5rZ9eeg0dQSTDxDDmJ6kZohxOWfAotFbydrfS47+wE1WdrmyW4GLNlKpOpXcySBYDJyRPbFLZFg13GcFn8HZCyZml6XCt2MdM0U8C6/UFx/KjcVX2LOK8YV1erfWV2av55QcmrlnNVq469i6XEulWMW6bOuSpOb0FpaxhTiu3rLBlTN9unVAjZRlWOWbZZHB9MyYsZtyw3Z6rdtPU06xq+K3530yp+ch5NkdN3ys/LzKjruEoA6DADtxPG1rWSeiyt1+9mJdBm/J7D8tDCDeenGCgLIcSFXvKcczcAvIVkoXLfe/8k7+O9/+uLnEMIIYQQQpyftV/ynHPfQzAlfoBQKeI2gA+dc3cvmmQhxDpMLXNxNqdqEtYP+kEx2d5Oj/tgQFXJCsnPqfRYOTJTccZHB3Efxmal+C6Wtzop9gnHC8oOFbWulchimxSsXJVrFeui4lWpdTlxLPP5wrbiGNlYojJYH79S1ULf0hevVSl3rRhDl6tz5a8VVx0fjdcRa63ZPlQd7V4UPnllVm1U8trN8Xzh+GVMHuMpZ+MQT0lVFgBmmXobrjEcr7dVZVU3qWiWEV2rpcxXzr/vVpseg/VYFw/baYfnqTeyUnQutENL6O44KXpCiEXW8smzRIp97/1t7/2b3vu/sLYdNru/vNxhCiGEEEKI87Cuknfbe/+3TRu893/rnPt3FxiTEOeCSbM9Mw/b6oW/XW7shsf75o0Qq7Wzk+LGen179JmZO7aMSmZOWjZsrlC1xymbEsjir0y1yyWZ7iCoPt2t3WIfV6tyuYwTFSL6ptVxXfNiuXEM6URhzN2UYUr/O2YMt6oYuSZPuqjg1X1bpdJW4EtV0fNvyXif0vhrlS/F8TF7l5m56TpaC156LTu8KXCm0pUVL0xttaokrE4ys+8099Sjd9/Wza8AAHpbewBS1ZLhYVB3O08exX3ah+apZ0oes2rnVbWS7jB567EqBrNpqez1emH91lbyXTw+Dvs/PbTsYLvFvZEpfBeytRdCbCrr/mr48oztj87YLoQQQgghniHrvuTV9V1rbq15XCGEEEIIcQmsO137S+fcDwD8B+999B6wbNt7aC4/JsQzYWB/quz0wrTX7Rvhsf7Kq8Ho+MZNS3ropcc92ZSEnWc2rTY2g9vxcZjOa2eF6zldWtuTtKOdSD41WiYFpIMwUcIMlbNSaJwy5LZ6uraRJfYnTIjIxxRLh3Vp1dIqrmthyjdcZHFc78sEldOItijtMpEkT4yV0PxiAAAgAElEQVRo1fYnnIqtp4OzaWE/nxYt7xunZGNSRZZMMZ8Mi771+Xs76e9Wt9fmiYrj0By5O/jc9k3X0en1re/YdqX5cTjfZBiO0e4cpfNUZeU63Y7tE76HyTjd414v7M9H4cQSMGiS3G0p8UIIschaL3ne+0+cc7cA/Mo59yVCvdh9BIXvrvf+7y9xjEIIIYQQ4pysbaHivf8IwEfOuW8h+eR96r1/fFmDE2IV+iby3NwOH165HZSqV1+7AQDY3guKXisrE8UyU74qKE9FLxalz1S0ZIfCklimltFOpDuIfaNpcKXkzaOFR1B8XKaItea9om+dCNGUGFEnckQFrL1othxNhG2feG0cky+VxBxfqX6UlOrzAbmCZ0kTdl+SHUuTcXJpgswxzGZBnfOZAjevbVBMaWNLa5t5dh1R4bSxdLfCs9HbDpElncGN1NfGR6WQSRrt3u8W7gvpDrZtTKYqWjsdhfGPuocL+8T7b9fc7XftvFaWbzRe6HtyEsbU74Xj9zrhGs0lCH+U8oTw6wZhVghxvVjrJc859+cIFip/573/BMAnlzssIYQQQghxEdZV8v4awH+9zIEIsS4DUy92LDiPMXi7+8G+ZGs3KG+5YfB0UpamYjxUpxdUHJYY62UWKP2d/dCnH1SbqJZF5Sr9c2rHGLPSLDjG21Vt6EQ7kVJxWzAPzixKFgx5fbMpMpDUsFk8t8W2NShfaUg2huj4QnPiMvavlVucLLFbSeXGkqIaFbtJqdglg+njYhlIKmiMr6vuAcfEuEgA6AzC98jYu/5usEfpbr8S9ummvtGSZWoxfXHZlMPtYw4+O/5OObao5B3bWKjsZvGIZrLc7Z/YNYZ9abeS0+mGc7fNJqjbMbsVi8kbWHuzk+6FlDwhxLrZtR+fVtVCZshCCCGEEC+WdZW8H9uL3APv/T82bP82gL9bf1hCrI7Vc0fX1AwaHfcGQT3pmjqXF5Qn3gxp+9tBndu+cRMAsHXDFJ+dm+k8Zoob1bnK0LjIGqWa1A8KUqvKaI3nP0V5i32ovDEWrSE+rd5WZ+oWx6lKiMUM2sr4OP+8YIJcq5jZtae4PRoCm1pnqlw+JsbRTUdH1h6Xy1YqLt8nZjUzJtJMnrum1qV4u/24T3f7tm0zNdbadnd7Yfy5iTKwWFqtO9gr1gNAxxS8+nug2TKVPo4VSGrxyDK5kxE31cv0PUdB1eJKO6boJSUvbM8q9+Hr9jV8tigMCiGuCWtbqABwAPadWyia6BDqCPzbiwxMCCGEEEKsz0Ve8t7x3v+yaaNz7ifrD0mI8zExxWI6C39vzKb0NwuxZpNx6V0GJDWLCt7uraD07Nx+DQAw2KPysxf3qePr6ozWdpZd2zG1h8pRu7dT9I3jKJQ2xshVqpy3Ml30hZum0liMAZuxlBcVpRU871yL8XT01GMW7GJ8XU3trUeVDkiedFNbF8c4WfSvi8qdKXbzKqvZtcKYO4UCFr6z7nZQWXtsLd6OmbJUUYF0/9tdlnTrFdcxn2Tj5720dVT24vccv8vMJ4/fXSwzNy2uNZWSy+M2ed+Zpc3s2jCmYSdl5LZaZYY1lb22KXtdlvTLYvL47+I1exQ+P1s0FkJsGOvG5H1v2Que8YM1jyuEEEIIIS6Bdc2QfwHEChdvwnzyaILM7UIIIYQQ4sWwthmyc+57AD6wxQcA7jjnHgH4K+/9/3cZgxNiFR5bYPnT4zAf9fRJCHbfPggV9wYWwO4y647BTpg+29oL06p7r/5+2Gc/TNf2tsOU30JZMiBZdkSHE9qJZCXEultFm8p1lf/kfBbkP+N0bDUtS3uRWbTnSPvUNiLJrqS7MP5Ujq3T2Ldpn5g8EUuH2bRqNSXLaVcAmJw8sXU2bk7TTjkNnQygOd0by8zZFCbNhXs2Xc6pWSCbnrXECiZVpCSXreI688+cTp2Nj2wsZcJEsa6ymlm4x9n0fA335T3lNXd6mRXMoEyO4VR1u2sm0p28/Buna4uhYDYvQ6I7meE3S52xSpuma4W4fqw1Xeuc+zcA3gHwhve+5b3/hve+hZBV+2/NLFkIIYQQQrwg1lXyvum9/4t6pff+5wD+wjn3AwCqXyueC7SIuPM0KCK/+8JsK8wwdncvKDPbu0l52bsV1KDtG6Gs1WAvtLWCV9h9UO2xBAYGznM5Kj5YLL1FFpS83NiY52LwfmWvEpW2TJFs+VJNcpWdS5E4UclA9fFjskCe2FFZnEyGIRlgau3EFLzJ8CjuQwWPylQ8fUxQSWXNaC1C02kmulCto4kxbUsAoG2KXacf1rV7VrauXZaFK74Pmi1TJY3JLQ2JKrWtTWV7s2BWjUVrnKjCRtuYMjHDdrJx2/GihQ2V1sXydYRJRuNpaCe2nAt7louBTrmrEOIasW7ixc/O2P7pmscVQgghhBCXwNoxeWdwkC84577qvf+nZ3QuIQAAX54EGePRk6DS7O4ElabbC4/57s2d2LfbDzFyVJWoxDDuDQ1WJCl+LPSdtlkKLcSAlXFdZVxdVJssNo9lz1ymPtXxe7Q0SSbGi0FVC+uqEl/59lpdivF109LahDF0QKbYxbayPJmV8Wo5LbME6fTDdUXrk6xUXDQwZnzdVml/kmIZk/UMFbtWpbbO5iyNVtrK5NfoK5No3p9aicvPGVVRlqqjOtqwT60Uxns9prqZGS033LNibA2WPzMqeOaPMpqE5aF9DdNMyps1H14IcY1Yu6zZstJltv6zavW9Nc8jhBBCCCHWYF0l75/QXO0CsJzDbBsrYAjxTNnrheCj7UFQYLa2wuO9Y7F4g+0Uv8bi8LNJUFaGTx+G9SdBsapjoIDMoLeKpZqyrFamgM2iEji1fcyImGqQqXadQcoabXXK+LoYx1WZJJexf0vKaTEjd5zFyg1DtvH4OAjtk5On1if0ncY2mRXz/lCBmpu61KpMios4O1tHhTPF1TGWLimqVDjr8m/tLFM5XHNSVHnNUZ2jMXMVB5krefHeVepZimHMTKpjzCXV1ubs41z9i/GMlYLHWEaqo81qrMVIzhgTWbYAMDVj76kZfVPJG0Ylz2LzssNT1Jvqt68Q15Z1lbwHAL7uvW83/LTsJy4D+OkljlkIIYQQQpzBukreWRUvalQBQzwz/pWJSH9yO6gxf/B7QR169bUQ37W3H1Si3iCpQ/NZUIZOngZVa2yxZqlsVKnWASnbMfnLhbZrGaK52kS/t7hsChXLnPV2vgIAaPdvpD5VqS0fFSuLlcNxsQzkMV9HxfLEVDuqdwAwOnoUth0/sT62z7SMU8vhfehYPB0VOy6nNqlzVPBimTSqcw1l0xY9+3gPLFuU92CSq4sW31bH11FFq1TUJuqSdHkW8sJ3Z2NKil75PeVjWIhv5FiaYv/s3jKucTph+beJLafxs1QflbyZnXpsXxkVPCl5Qoicy6h48Ras4gWA+977J8v6CyGEEEKI58NFK17cQ3i5ewjgNoAPnXN3vff/8ZLGJ8SZvBaEI9y+GZSWm/thxfaeZXP2S3UIAMYnrMZA77vSQ46qXV51gJ/TtnC+3vaoWA8kvzcSM0It7s5Vnm5AHmvGqhKlSjcdBVVucvwoXYd9pmLHGDBWoMgrUTDmjnF2VJeozrW7/aLNP1OpiypdrdZl6lydEZt85VrFecOCVdSw+DMqdr6q+sGxh/FbnxhnZ8pnzEpl5uyihFXH4FG1a4rJq2PxoqJnff0sz34t4wJjprVnDGPYZ56PiWqlqX3TsV3rhPGVWRyixULG+D1rW6d44I1teAeLAq0Q4ppwkYoX+9772977N733f2FtO2xuzrwVQgghhBDPh3UTL2577/+2aYOt//r6QxJCCCGEEBdl3enaL8/Y/uiM7UJcGgwwn9uH2TTMT42HZjzLKbLMFqXNgu+t8u+cOA0Wp2RzE952sY3L/e0wRZqb17LkWbcqy8XpzVg6KwvEjwkWVTJFnLIclcbEQJqmjRYd1ZRlnjgSS7axvJiNsZ56jWPE8qSJegr2NDiVWVufhM+2rprmTGXhmIyQEi+iAXM8QWlLMq/Kw+Wk0mr9Yvx54gWtU1qVhQpqWx2/OF2bTmSJO1USx7ywgmGCBe8LjZQtiSMzNq4vpW3ztF17PPvW5okXQ/v8WIkXQlxb1lbyzth+a83jCiGEEEKIS2BdJe+XzrkfAPgP3vvo0WDZtvcAfHgZgxNiFX5rMfkPnwQlZPcLsxGZBEWk26VNRpYY0Q2PfrtjSg4VPRNr2u3FfVqt5r+Jun0zGR4uqk11UgAVKZr/5upQbQGyjHamDrWsHBhiQgEVNjtuprQlWxhLJFgw+y3LhYXd2wvHsQsqr6vJeLgq8UUFLy/tNY0l1prLjs0aysul66lV2CoBI5O/mPjQ6pZJIbzWQr2sypctPU+RQFLKZU3HBQBkCSTxGWEJvIZrTIdnwgWvJ7T9TvhOB/ab/CTzSzlZ9F0WQlwz1rVQ+cQ5dwvAr5xzXyLUqt1HUPjueu///hLHKIQQQgghzsnaFire+48AfOSc+xaST96n3vvHlzU4IVbhM3Or+MbToIQMekEtO7Gq7b1ukD16vRRfNzDpYzAIikvH1D7G6rU782IZQFT5/Ly0sRiPFpW3bo8xebRxCS1jtLhvbiJMRa3X27G+/ca2yVA3mfuW/6RPMwReUOfSTuljrY7Nyhi6aEzcZFbMeLppWXptVvQtzYNTTGF53ubx1zFy/F7MtqRBkaQK2qpa19B3Id6wOn4+tqjq8fugOlrZpOSK33xaWaVwWzSCTsefzUoLFcbk9bvWThctVUaKxRPi2rP2Sx7x3n8C4JN8nXPuq977f7rosYUQQgghxHpc+CVvCfcAfPcZHVuIgtcqwWViqsbJMCgh04a6Tv1+ePQZc9ftheWk6DGTNpNGKsPkqMZZfF9eNq3TNxPhXmkwzCzb/l4oa9bffS3u0+4zAzdku9I4udU1FbC1+j/XeTQTzuIEp6XRcGxjRuti/FvM+I3xdNaybJepdrMs1mwhzq7K+KUKGI7HzFIbQ9W3vtdAHm9YKnkxU7bTLVogZRBTOWWsnFsSd5cTVbrqvvF+5evqsaTtZUwmAEzH4XPMrjXlbtpQ1mw6KTNueelU9DrWtp3kOyFE4szsWufc/7AkCy5/yzk3O+VnDuDtZzpqIYQQQghxKqtIA3cR4u3IQwA/9d5/Z9kOzrmfXHRgQqzK14LQhRtb4W+WQS+0fYvBOy0mb3snKGzbu0HpYQk0euDlYhE99ZhlW5c+G+yk+Lrtm8FFaLAbsl97O6HtWjYsFbz+zT+O+3QH++e46tOh6ucz9Y/q0rwqGTaPMXRU5xri66jY2bZUiouefnnZsTKLtvaQy1W5mFlqbcqmreLesi8ixR/y/jO+zuIrq9JrANCxmMja76+VlzOL46OvX1DUXD3+qDZOF9ZF5bEuy2bl5iYnyeNwMrTSc1bObDIK7XgY2pPjlLE8Gtn9mVHJM+XO/kxn28n+bO/zlkncE+LacuZLnvf+p9XyL6xu7WncvdCohBBCCCHEhVjXQuWsDNqb6xxXiFXJ4/D2B0Gy2O6HlVuDoM5sb4XHu2tKHuPwAKDXD6pPfysoOoOdoPr0t03xMXUu98ZrdRi3Rw+0UvXjcjjOru3Dagnr+o6vxyxWz0jK0Wxk1TGskkZso5JnqpOpc0Aee1e2debsfJJUp9O87YCzlLzSJ89F/7+kuDlmQHf7jW2TkpcUvMoL0L6XQm1k7CDCdbh51adS6/JtdaWOyckTAMDo6MDa9KszKnmjcE+p5I1OQjs8TjF/Q8sUn85KWS4peqHttZPiudcOfT+XX54Q15Zn9T/P95/RcYUQQgghxAqcquQ55/4MwKfnPKaDokCEEEIIIV4oZ03XPgDwiff+L85zUCVeiGfN72V131naicHn3Q6nZ9vWcto2TfmxnFkq9UWzXJuSNTNjTs2GdbRBMTuOHqcAzbKjmNo1k10mOzBRIQbgh+m7VjdNKcbz9G8U+64Cp2fnEyZCHNr5nqY+I64LbT1NW0/BhuNV07XR8qSyPslsUeqkCV+ZCBeGwJX5MY2YU9kxlp/L7FCq6VmWKmMCRp1ckX+OlinLDI+zMXgaEHN1w/jTLjbdbPcpTtMePgQADJ8+CstH6fuYjGihEu4dy+INj60dpns6ndIMGTbu0MaEC3u0B9lv9F1+TocRQlwzTv1fxGLv3lvjuEq8EEIIIYR4gaySXfvLpvXOua8i+OEdeO//k637FoCfLdtHiMtiu+HJZVD6dFZGmlPB62WJF51OaZ1B49mZqSpU9nxn8USLFh6m5LVT36gydWi2G6QXGg9TYWsdJ7UpqkGmwtEMecFYNy87xjJjU5YHC0oeLU/yJIqUcEGz4jKJgmpcaYZs6l5VfizapFjCRJM6V2/zDQpYSqygPU1pZNyuVDoAcFHday5RRtUu34ffTWuJgleWKGuWvurx599DSrQIiRXDp18AAE6efGnLQckbDzOrmWn5zE3GpuiNbHmSl5ez64iWKaZAW9u1hIteJyVebHfCTkxSUgKGENePtRIvLFbvPoBXAURzLytx9k17ARRCCCGEEC+Idcuafcd7/w0gqncR7/0nzrm/BPBPFxybEAvcNKEiN32dWqmn0SS0w1GQLBjH5FjMfZBZnJh1CsuZkbHZWMQ4snmmpszLIvG1yW9vazd+jrFkVfmsujQa4+LCJlP7TNFrtW0bVSfGijWY8NLIOFp4sJ3kJsWn26FQefN52TFanNQxeFUZsiaomvF+UWNyWexiLEVGpa0yNE6WJ5mpM9W4GLfXrfqWcXc5UY1jibKG7zIZGzdLX3Gf7D4x1pIK3vHB59aG5fGxqajT/N7OinZmajKf2xwqdlTwZla+jFX3uD1X8gb2WVYqQlxf1rVQyTNumzJpb695XCGEEEIIcQmsq+TlL3auYfs31jyuECuRix1Wux1DU/KOh2HFaMTMzbA+V+22doPpMcuYUVmbm1rHjMc8DivGUE1LA1+qUZ1+ypR1UZmyjE9m27bKf3K5KhczVKNy1xwD1qTkxTJaMa5uWmy3hfK4Ua20Y9i+eXZtypC1fZmNzFhD37HVi2XHippwS9bXZsetKkOWMXmFGXIVT1crerFv1m9BwavWFzF5Tfeu6MvYyaSSjs3s+Pjx70L7KLQnT0OMHp+Z/D5RGZ5RyZuVSjHVOQBoW8zdfO6K47BPx9S6bmaGzCp+/efrwy2EuEKs+8/fOef+yj77bOWfOuf+M4D/cuGRCSGEEEKItVm3rNlPnXPfc849BOCt5RTtu977v7+0EQqR8dj+pDhOYlb0ButMw8aTcWgPrdPNIePHkjLW6YadBlbGjOoSMx2pvORQYUlKXqly5cXuqdzVJbZijF5DlmdS6sqyWqnD8qCqFNsW1Cy/RA2sdqraVtliUclbjHIrzx8+V8c7jcoXj7GLKat2UcmLu54RM+ewuD0+A/W9blDyaqWTfajgTU6S593wafDDYwze8eOQTctSZSl7ON2TWhmmshfj7zJVjutcjMELbbfD9WFFnlneloInxLVn3elaeO9/ZKbHbwK4g2Cc/OkKdW2FEEIIIcQzZu2XPCCaJX9iP0IIIYQQ4opw5kuec+6vaHa8KuvsI8R5eJj51dL0tWfTW7RSOToJ02tPn4Zpz5tHw7jPntlVcP6LZcxY1qy2SwldmXRgZdL6YeqVCRd54kVnsFu03cFeOH5vJxyLth/tVK6LSRnLSm41TbymGVdatFT7zLISZZZY0Z6GKepU1qw0VJ7nZc2qxI5LIbun9ZTrMjuU4rq4T8wjmRdt3N4wXVxPvfomA+h4zVVCihlAj22adnR4EPc5fhyma0+eHFiflJQBAG3HxInMQDlaqJT3oN3m1G5TTlugY30s6gATy0QaNxgoj2SdIsS1ZRUl7x3n3McAHq14zNsA3gGglzwhhBBCiBfEKi9530aIt1sVh2bRQYhL49eZu8Wr9nnbFIuJlTc7HIYVB0+CMrX7KJX42tkLakxvYOXHrIxZp0sbDpYuSwH/7MMyZt1BUO7623u2nMyQO1WiRatLc1+eb9GwNyp53LasnFlTebCoDPZ4MOubK0fhPlCZooFyNEWuSqSF3StVq1bL4tiyMflKBa2Vtny/yqYkGR136o5nj6k+NvJSZaWBcTKNXrRLiTY0M5YdC/djOgrPz+gohB2fPEl/9x4/Dgre6Pi4OF+7U1vmpLEuJvBYwogrbVKa6FjCRadTqn6Hx+k6aDM0lJInxLVllZe8+wDuAjjI1jkA7wO419D/FoD/5+JDE0IIIYQQ67LKS96H3vtf1Cudc//De//Lhv6/dM59ePGhCbEa5n2M8YylnkwBMSuVx0ehw+5Bisnb3glKXqeXYuIAYGsvqHK9LZolD+K2dpe2KH1ry5g8xpEBWDACrtWspIhl1iPVdUU1y1S5uL1BuXJtmghvlfsWQ2CJMpYzM1WxMyjX52XNTNU7SzUrFcMyzs1XJcQKVY7jPKuEWJMBdMM227Aw1lr9i+XaoqqZxSFa6bGpWaVMR6EdnxwCAIaHT0L79Enc5+Tw0I4Xjk9FOCp5jMnLrHmm0a6HcYK+sc0uKap7VPDYzma85rQPlbyJ5lWEuLac6aTkvf/ReQ+6zj5CCCGEEOLyuJCFihBXgUcmjuxEkcSUEB9Uj+4oLD89SirK44OgzvR6SY3JoWktVTtgMSaPhr15hmwcwWyx9BUAtKko0eQ325eKVGsh47QsieYaMnKpxrGtVcDi+NZnNjm24/H4Fg+XxeQxLjCqcysoeTxnjH/j+lhSLFcv5+W6KuavNojOPy+LzavX55/rGDyWcGPcHZCUu8nwCAAwPrH2OLTDw6ACD49TjCdNtKncMQObsZ2+KpcHANNx+DwZl4oey5uVSh7j9cJyMkd2xfZpFuJIZXskJU+Ia8uZSp5z7s+XbVpjHyGEEEII8RxYyUIFQFOZstP+PnxvyT5CXDqfm/jDLNuu/enScuERndj6URacdHISFJ3DpyEOzdGTrMpcdVkZqlrJW4zzSioNFaIYi0XlyDJZmWXbzuL4YuYtt1l8ne+yNNqOjWMr22dgfbdtWxYXWGPqnq8ycOss3tM4U9HLjreKgLSQgXsOVW5Zn8aMWXoBVrF4s0lYz8xZICl3o+MQZzc6sli8I1P0hoxdTPegZb51VPLa9qzw+ibjcB6WOQufbZ0pedMJ1cUyRi+HHnrtqmbZ1Er6jTKfvJFd/mMpeUJcW1Z5yXvPOfedhvX7zrnvL9ln/wJjEkIIIYQQF2SVl7wHaLZKWYYD8O/WG44Q60Nn/+0lysV0ljawQsBwGJSddptZpKFPjIua5f5pVIwqHzj6mhWed6xaYTFTVJm4npmbWXwdVT1m6VLJa/eT/x5QxQCaktdU3WEZVNpShY1OcVw/T2qgM7Wsbv2sip0rlDZm09YZxQ3ZtfNSdaszc5tUudqHb5miyixZYDH2Lip446DKTUy9A5KCN7SMWXrfjYdWEcS+/9wDL2bTWhtj5EzBYwWM4VGK0RyeWAWNURjneByudTpdjMljDF5dDWNaVbo4Gad9RtktE0JcT1Z5yfvovNmyzrk7a45HCCGEEEJcAqtYqPz78x50nX2EEEIIIcTlIQsVsTGwfBNnuVqxPFRYznMqaB47sWD34dCC3y3onesZFA8A0wnLXJVB+5wW5NQiAPSrKUSWN6uncZGb8FpShmtb8H43LHeY0NFUFqwyBI5WKpXtSrFPnCJtnlbNWZjadZNie5oiHTesq0yRq2UgJaRwG6dVF5aL6drK7NjTPNjOG6dm05g4LRvbaJMS2tFxPl17XLScpuUUfrvTLlogTdfG89kzwmOcHIZ2eJwSL0ZDmzKe0l5lViznxsacpu12y+nasU3Pcrp2mCUXje323LRHTQkYQlw/Vg/kEUIIIYQQLw1S8sTG8GsTe16tFItOq1T0gBTUTkVvPC5VFC4zKB4AJgyMt2B6qjXTcVBnBuOk0sxMMerv3AQAdLdC8gQVvVa7wazYVL7aDDmqaVUSRzhPUKCivUqnLFVWJINUyRmxzBkVttmialarfvN6eVYqb/nxYtJEXC7b/FzR2oRWJ1WiRDGmhUQOXxwjJlmMUxk7KnbjoalzJ2y5PvWlckezYibc0E4n2qR0Fn99zqxsGY9LBe/kKBx/eJKpi7RMMYVwUi277IGty5jx0vkcU8kbZ4LnVMqdENceKXlCCCGEEBuIlDyxcRxaaNyNXlneLIdxe1TyWq1S2aNKNM6kEX6molcre7nR7XhYxnz1d4Oi1zNFr9M3exRT9oBkneIz9a0YMxWw8VG2jqoVS5SF2LAUS9eg5LmyhBjVOVTxfWFVFfc2K5W2uZVAa1bnqBCWimEeu7hgmbLE6mSWxS7W0NKG+1JZnQxzg+NSuRtFBc9Uu0muRJZKaozBs7i7bo/3OD1XVPD4DCxT8PjMAEmxizF5lXUKVTsA6PXaxbo6bm9WxaMCQMeGt2uH6du2usyZYvWE2Fyk5AkhhBBCbCBS8sTG8ZmJMjet3e5S7UiSBRW7qUkgrerPnTmVkXlSdcbjMl5vWmXg7g2Tkrd1Uip5E8vmHFDR294L7VYyOu70Q2myqBBR2avNf/M0YSpg3mLXJma2W8etNZCya6vSYg37LDUtrlS6vE/sGzNoG8qNVXF0zFRmFiwVvVz9SydiVq0pYabgTUdhX6p3QJ4pa+bHsZTY4nFb9jB0usxyLmPwXCxVl54NqrnD43D840MzPz5h/OZiqbL4DFYKHjNnqd4BQK/fKcY2peI58+W+mWg9sN0Zm9exdo+Pl7X7mXD5q4bbLIR4eZGSJ4QQQgixgUjJExvLlxbGtWUSxqCbVJS+fe52TMmrPPWo+k2zFEWuGzE2zzIaR6Mgf9D3DABu7JexXozZYlsrVQDQ27ZM1Qa/OttgbdrObFo4qkGneNItlAyrYuSalLwlJcSWlR/Lx+pev40AACAASURBVOmrkmfxGEVMHu9LreCVXnf58dNY5tW+ZQmxPGN2MqKiNi325Reee94x5o7ed612GSPJbNvxKH3fo2PLprUYvJNjfv+lglcqeTZ+U+P47FHB6w+S916v37Vzl6ofn8FxQypt247Xtz/lt23ZrPai6tfO1TspeUJsFFLyhBBCCCE2ECl5YmNhfNENE1wGnaR2dNtWvaD6M4eVBaJSMk2q2cTUkiGzFC1Gj8penok7q7If80xMIFewsmzLqnIG49Oo9rW7QbVr97biPp2exfG1y4oLMXO2wb8uxrlVnnSxb6bkJTWuVOXmVXzdvLEiRan+1cv5NVOxq+PquEz1LJyrHFOd2cp4u1mWMTurMmapzjHurtNL948KHn3xeB6qsszIZfwdkGLwRsyinfCaeUvK+LvQp1TlooJn8XdU74D0XM6m5bN2PLSKLRNmkmfXaI8cY/Oo4DHUj31H2VcnhNgspOQJIYQQQmwgeskTQgghhNhANF0rNp4nNiO6mxVv74+bp2t71aznJAtoH45LKxY35nqWpUpTcZyC4yxtXQqLVhj51OjCNK0lIdBapdMLxsntXjJQ7tjULadwWzRDbjBUXrA2icbDw+L8ZVmzebGNVifROoUJDFkySG31ks5bHiucm9O0NJQ2+5PhsFifT9fWx4/JCLNZ0eb5I9HQOH4PnK5dTK6IFilMsLCxcHr25IhJFilpZjQM31k9TRuvc7b8GWm3S8uUbs+mkLNkEN5nmimfnITzcbqWiRfTbL6WiRdsB+aZQo/lMUv6LcnzEUK8/EjJE0IIIYTYQKTkiY3nxJSKPMDchBB0WlTcmCgR1lNom86SMkK15NikD8byd00qGWVK4dwH1YcqTVRnes22HOHclTWIqVyd/ontO7DlPPHC1lWKHpM0ckUvWZgwwaM50aPJrHgVa5OF8yzYroSWyRbhOKaAMWnCWip40fIkk8aiChptb0LLexrbzBmY66KSasunGRuPTky5sxJlx4elPQqtc4DFcnj1+kmDGTKTKZhoQcsUJoPkxxrbuY6OwtiOjs22Z0IFDwtQsevYfRjYb3sOgfscKvFCiI1FSp4QQgghxAYiJU9sPJ+bYvFqplhsmcJCMaYdDWKtODxLQmX7TGwfqoDHUbkL7dORy/pO7LhlvBVtMRgT5vJ6apXNSoo5Y+mwRePhWlnr0G6F8XvdFL+XYs3K2LtkcdIUM0dlzWLkllibNBk4xxi8eXPMHJCZRNem0bPlFjStNsuOlVYnrlL28vtZW9iQ2hYFSMbGR08rBe+EcXeLptH18Tl+WqbQVJvKLgD0++EZ2NoK15FKlzkbW7pPLI92fBzGcDIq1WSKlp1cvYyxeKHlsziysVDZ/rWUPCE2Fil5QgghhBAbiJQ8cW0YZWIT45HYUjRpT0tT2aLCV4xlMuXOFL3HJnwNWlkm7pQF40OnTieoQt1uqeDlClDS3FBsm0VlarQwqDpjlW2nWgaAlsXn1Yqd9+XF+sIMueFGnEEdk3favjG+zmLlYnxdi6bU87IfsgxZZsa2yr9Va+Uw/0zljrF+VPCYMQsk5e7Y4t8YezebldeRx/xRzWXMHfvWGbSDQfqVu7PTC+u2e3YL2ravmWxnZfIYi3d8YiX1qjJm7croGAC6NEG2c/NroGf34xQaKYTYUKTkCSGEEEJsIFLyxLXhaRZ7dNMUPCp29L6jojdvEJ+4rhXVt7Di0I6Vxzalc4UP7dbY2sPyoA2q3NbeXjiPKVStduVVl6l/VOVadZxdLEeWBuXa5T/32r+uLllW7E+FrWMZoJWXXOGtVyl484XjN/jwLSlVxr5lTF4Vz+ipnln2LjN2x0kJG1um7OjEMmNNJWOs2/AkyVrDoZVJqzzvUqgfM3IXFU+u4z6dTplBS/UOSApe17JpeQ/GozAWqncAcMxs2nFZIs5VsXid7M92qnuMxWPpM6rMXywmRgshNgwpeUIIIYQQG4he8oQQQgghNhBN14prw+dZ4sUf2HTavGon0dR20XJj7plMEZYHFtjeL2dvi3N1LZ6/95jTjmG6kNPD08zFdpltiKvMf3MWplyrJAefmfy22jTZLadGF8qdFdOpNPFl+bJybJxS9sUUMqdr58W4fdP4q3Xz+tobzIrrsXBqdzKmsbJNwVoZsvA53Pc4PTss7VDy76G+h0yw4FA5JZtP16YydmZ+3Q33hYkWtEnpb6XpWhpke4YKsGTZcWmXAgBD8zuZzsrp4Gid0uZybhsTWj5rJ5MqYWj1XBohxEuKlDwhhBBCiA1ESp64liQFj5YppkxVil7TPsTEGuyZipIrhYTJGLfN/aT3lIqbmRhn56GSl6tWOVED8t24zlNJq8tpUT1rKAcWkwSmLGtGta5UBcPnZiXvNJuUqBRWZdp8k7VJta42R6blyXyWK55MsCgVPCYsUBEbZ2XHaIMytsQFKnf1fQNyFbFcX9ujlPe2NL2mcjeIRscsWZZ5nPBZqxS8kxOONbvmebO6SAWv3VpUSbnPkSl4h+PQ/nq80FUIsaFIyRNCCCGE2ECk5IlrySwqeau1+WeaIVM82TY15Y8yAY4K3k3GRdm+tK94csLOKe7K++NijHVcGpWjbr8f+9AY2M8Xi9rX+Oo4Myp507KkWFHWbFqXVCvVvtr65LRts9miUhjj26KSR9PiadHOsvpyU/tMBYzlv6iE1kbE+bYUO8cti4pnUilRteW+3UyVo0UKY/Co3DHurk1vk+zr4fhHFh9I6xaqjbl6x2eta8dpt8ryexxjbpJMyxQ6sfzWQhQViyfE9UFKnhBCCCHEBiIlT4iMWsHLki6jglf3jdm22Z9Mr9m2LhWYVtmX+1JtAYBHT4OS026fWNsqWpIrYZ1eiNRrd4L606IpcUOJMpoI18dZZkQMAFMzFK4NhmPs3HS2sA/XMY4u9mX2cIOJcFQX4z7Mei2PFdaVqtyiOrdIUkVLlS6ZF+dj4j7lvlTumDlL9Q4AevaZyh1LlEUVltnUk6RIjkZU8CymcMz4x8Xroaly/QubfUcWvzccp+s4GoXPD4eh/ZXMj4W4dkjJE0IIIYTYQKTkiWtJu86c9M3L4yz7dTJv7ks62TG3WuU6no/lp/q2oZsNhDFYT4+C5NL9Mih6VHGa/PL6VQZrq4qd6+SqX7drxwkqU+3HhyoGMP+8kAVbZbayBRbj6GbTSjFsuH91WbC4bxVnl/eNXnGVj10dQ9d4/Ia4PdK274T3nQoeM2dr1S70LZU7wvFP4/1anvFbK3jdrEYZ1/Fa2XdsD2XMpB2l7/vLk7Duf44WLlEIcU2QkieEEEIIsYFIyRPXkjq7ltQZtKMUQoXhrNy3VgOnDQoVlbxYmcD+rOrZhq1eOghVGgpQh1aUvvvQFD1TlOoYvbCPZfxa3F27wWsvVqloN1eV4L5+nmL34vEsi5fxda41K/ZtGkutkjX5z9Wq3OKxVluXr286/6yqFJEyZC3WLVPNqNjFGDzG2dlyVO0avOmYPTuNWcJlJvAo8+6LWbSV7x7vUzt7wNrVw8ZKHWwPLVv74XH63n+jbFohrj1S8oQQQgghNhC95AkhhBBCbCCarhXXknq6trZHsZm0mGwBAI9tpm1STX91G2btyLQ6D2f4OOPay3bOA+2BNH13ZNO2D78MZsksT1Z1DsfIjJLrvrP2pNzWLpMFYt9OVjatSrhg8sZCqbIGM2ROXdYJF0VCRLzvpS0KabJFOc0qpbyetI7TsvVUaD0lC6Tp2FZlYeNa5fR2fh11kslkUho1c5p2MlluBVOPO5+W5zVznxOLHXhyFNovDsP6z0/SNf86CzUQQlxPpOQJIYQQQmwgUvLEtYQK22QxP6FYn9WIh4kl5wpkpyly3/6cokI4jwkAuZJXqky1Ye/RUVDiWq2juE9tIry9F9r+VmWPgsyCJSY7BFWuVvRcJoG1TNXrNJgrF8d3eZJA2D9aqcxq1W9RAaOXcru6t0xGadcbsJjIwXE3q39Uyex4ptalZJZ21hfFOGvDZj8v73l+rUywmDDhIiZIlEkW+fFr0+XaEgZICh5V3UdPwrPw+ZNw3N8chWP99/HitQshri9S8oQQQgghNhApeeJaUit5dYwel4d5TN4aVhSfc38Lh9u2slM71m7300G3LJyuZ/FjVHSo/kxMzTk8THJNbRcyr1SzJkNgxtW1Lc6uZfYorba1rTwWrLRXYRm1ej0tVoC85Nm0Gltp1Bz60GrEYtiqMmYLhs3V57CiXN9k67KgkrXKe5wTlcdZaepc32OONYw/fKYtCpW32sy5yXy5HhvhsQDg+CR8rhW8//U0HO+/yfBYCNGAlDwhhBBCiA1ESp64llDJY8wdvWZrBe/gkjIUo6Jn2Y/9dlgxyLJrB2aMPOiXpscMF0sZlkkNOj6muW7IvK2zXWeZ2tTfCioQ1bio6FHJi4peZoZcZfLG7FTrWyt6QFII67H4OKZkCJxUv1Lla1K86jGQOmYOS+IHi32rGmh5fN2sMjJOKiMzjReVPMbe1crdYgzj4vWwD79XHuvoJN2nx4fhXF9a+9sjKXhCiLORkieEEEIIsYFIyRPXEmbKsuwYW4pkVPAuuyQUFb2uJch2WklBsoTP6JfHGC16vLHclnOLXmsnUfUJtawoIM2KuLHQpzcwRa8b/vlTlWs3KXntdtHWcW9NMXNxf8a9+TImr8jejTF9qyt5JGW7lqplHvOX+pb7zOel2liocuPS4y4uV3F2+RhTJm5YTrejvE+5sFerfiOTlY/pgZeVKHtknx8NQ9/fSMETQqyAlDwhhBBCiA1ESp64llCh260EIyp5z7qoO6sRdA7TiVqm0LXctOi7sx3+mVLR63QW/zajGjQclopeGWtmHnrjoOR1e+aB16WCZ8pbHl9XrYstPelaZSWJgirurY6ZK9ZV+zi3qHzV+9TVOGJ836yMnQPyqhvNXnezIr6uWcGr1bpVqH0Mx5nx4thi705G4XxHFgh6aGrd4yy1+4klVH9h7edpkxBCLEVKnhBCCCHEBqKXPCGEEEKIDUTTteJawwSMXftzZ/KMp2lrfpXPzJqxbcs1+7Zw2pZlz4BFSw7acDARI08O4LTjYBSma/t9Ttcy+aG0bgHyadpWsS0mYjRM10ZT3yXmxIVBc1XqrDZxTsuL+8eEC1qa0OqkKqcWOpf7zqvpWiZZACkJY5mBcbMNChqPSzuUOqkCAE5Gts7aw1HY94m1Tyfp+I/tOfn1JVn6CCGuB1LyhBBCCCE2ECl54lrDBIv+c1bwmqCq16kUvVo5olkykFS9us84mvMuKnksl9XrWQKGKXlU9PLEjpiM0aqVvNLmpSg7Vql7tbLXlHhRW5w0lWVL+4S2NnyO5sVVabQmosVMNDZetKWpFby67Fi+nYod7y3v/9CSKo5jUkU6D5W7Y5OPD025O7Tn4Gmm2inRQgixDlLyhBBCCCE2ECl5QgAYUdEzseamtc/aSqWJz0zRqRU9ciNTkLb6tDQpFT2qTLlCNY5qU2hpyVJbsxRKnh23XSl3CzF6rTxO0MZSqXxNNiu1UpdsVha6LvRJ10gLlTJGLz92rTguxORl9+msWDweNt+Hil2t3B2N2JpaN07HXKbcSbUTQlwWUvKEEEIIITYQKXlCICl2N6vlF8l/N+Pbtil6c2+KVUM5LcbptSvFquxrGZ+mQHXGoU9U8NpU7RYzZXk8buP6TmdRyavVvVrBa8pOPQ91/F5dHoz3JD9PPYZa6SyVPMbklX1nUf0LLZVRIGXKUrk7GZfK3ZGpdo+zjNkvTMG7Cs+aEGIzkZInhBBCCLGBSMkTIuMqqir/zYrRzxiL5pOCRFXrJmMKe+XfbbmARVFvbkrUxDz1WpbdSXu8Trshvs6Vy5ZAG5XDdsM+MX6vapsroDVvayyXZtSxeXUMXdm33CcpeLwX6Z7yvkxnpXLH5bFtH2amivzMTFmzKcQx1Tr53AkhXgBS8oQQQgghNhApeUK8JDBGD4e5YmXxY7Zqbyt8qBU9ICl1jDWjMuWnZQxbUb3CPtbKXVPfGm66yD6nqX8L3npUKqtlYFGNm1Zq5jT3E+S2qNyhWB7a8jjbx5JqwYIWFpqnjFkhxAtFSp4QQgghxAailzwhhBBCiA1E07VCvGTEaVsgTt1OWeLL5m13bSqx101/x3G6Nhoc09S3SihoKiXG6dNO2xfLTa7Fdf5DVQ2s2N50LiCzbGlxebFPmq4tjxutVbIp0noK1vyTMZ039J3XbejDKVm242wfViubVpczar48IYR4LkjJE0IIIYTYQKTkCfESQ1VvSHPfealU7Q5S30GvTGbIbU+ApISNMzmKyiAVsNaUCluZrNFqkNrqVXViRL5uEdrFLParEyvqPlzOjaDHUbljWx6/yX2FXsfMr+DyaIlq18Su/Rm9y2Nm+9Qq31W07xFCvNxIyRNCCCGE2ECk5AmxAfzKbD1ODkPLGL1JZvOxOwt/0233WcYsrKei12s4Lu1DYtzerFTAUrxdOk8dR9eug/JOgV1rNW6amQhThauVu7qdzhf34brJCpYmsyUqX8eVLZBUvU51qVzfpPpJuRNCPGuk5AkhhBBCbCBS8oTYIGi6+7kpev86k5BesxRSqnJbZpjc74btNDzudxfNkMfWsnzXPMbqhfVNMW2tqORZPB/LnWVqV4rpK5d9VPJCO81OwHX1OU/L6m255vY8zM6hvNUK3qGNWeqdEOJ5IiVPCCGEEGIDcct8qsT5cc7pZoory7803f4Pt0K7Pwhy1o7F6G316IWXZC4qa7WXXmrD9jz2b7pEaSNNKlqn1bytKb5u1ePm/eqYvBjzt8K/2EkVx8fs2mE2psmS4zCDVgqeEOJZ4r1vnJ+QkieEEEIIsYHoJU8IIYQQYgPRdO0loula8TLxr8wz5VUzTL5ZTd8CwMCSMGpj45gQUZULC+tQrKuNiHMz5EUrlmbyqdfzTAPX+0yWTNcuO2beh9O0T+16Pl/BhkUIIZ4Hmq4VQgghhLhGSMm7RKTkiZeZf22K3n4/rdvrhz8Oqej1zO232y4VvvzXCFW+OkmjVvaAxTJjpDY2njUoeWclTeRV22rFblyVJlulRNkJEy+UTCGEuGJIyRNCCCGEuEZIybtEpOSJTeBm9vfgn5iqd8Pi93Z7pbLHNrddYWwcFT2WJuNyHr/Hz8n0uG7LMmrAoqXJKiXKyKxS7pqUvGV2KERKnhDiqiElTwghhBDiGiEl7xKRkic2HWbk7lfK3lZWCq3bDi3VvdNMin2MlfPFcq365WbIjOkbM7P3HAbHdZ9a2QNKk+PTaFL8pPIJIV4EUvKEEEIIIa4RUvIuESl54rrxR6ba3e6mdbv2eWCZuAMrp9ZuMTM39WV2bss1Z+vWCh+wqPJR2aszcnN4nLpMGlXAXP1rWpfTlIlbr+NyrfaNzvEbQmqgEGJVpOQJIYQQQlwjpORdIlLyhEhQ5btpSt62tVud1KdncXs960u1r21/fjKer664ASzG79Vee019ayWPcX1NFTWWVcVoqpKxTPWrY/7yeL/DKvbPLAnRbbhW9pW6J4RoQkqeEEIIIcQ1Qi95QgghhBAbSOfsLkIIcX5+PStbjEKTmy3vt8P84w37TdS35b5N39bTuUCayiWciq2ndFv5Ck77LoSnhA2nzPRm07dlokdTqbVldjFNhs1MyqiTMc4yYxZCiFWRkieEEEIIsYFIyRNCPFfy5IHHU/swLftQ7aPSt5MpeV3707Qd7VfK9Z0qaaP+nHOagteqjl+rdaeZLtfHiGPLFD0mWFDJ+9y23VwyViGEOC9S8oQQQgghNhApeUKIKwfVvmVK32m8Zn+6brWS1Dao1L/2KWpZ64w+57FQOe08y5BNihDispCSJ4QQQgixgcgM+RKRGbIQIue0+DopdkKIy0JmyEIIIYQQ1wgpeZeIlDwhhBBCPG+k5AkhhBBCXCP0kieEEEIIsYHoJU8IIYQQYgN5qXzynHP7AN4FcADgDQAfeu/vZ9vfB/AAwB0A9733P19lW8N5Vu4rhBBCCHEl8d6/ND8A7mWf7wDwAPZt+UMAr2fbP84+L93WcI6V+zbs6/WjH/3oRz/60Y9+nufPsveSl2269l3n3FsA4L1/YOvuWPt6pbg9YN8zttWcp68QQgghxJXkZXvJe4PTs845vtzxJeyg6nsA4NunbasPfp6+QgghhBBXmZfqJS9T7wDgPQB3vfcHAPYbun+JoPKdtq3mPH2FEEIIIa4sL1XiBRAVvLcB3PHe37XVt0/Z5bRtF+krhBBCCHFleele8kzN+xvn3B3n3GcIWbYPT9nltG0X6SuEEEIIcWW5ktO1zrmPnXOf8SdbH6dT7WXvAMD3ra2nWl9BsEE5bVvNefoKIYQQQlxZrqSS571flhTxMYC6Ptu+9/6+c66eat2H+egt29Zw3pX7CiGEEEJcZa6kkreEBwDuVuvuIL2A3XfOvZ5vy4ySl25zzr1ebTvtOEIIIYQQLwXOTHxfCkzNex2p4sXH3vuPbNs+wtTtzwB8E8CP6Xd3xrZ7CGrge2f1XWF8L8/NFEIIIcRG4L2vZzkBvGQveVcdveQJIYQQ4nmz7CXvZZquFUIIIYQQK6KXPCGEEEKIDUQveUIIIYQQG4he8oQQQgghNhC95AkhhBBCbCB6yRNCCCGE2ED0kieEEEIIsYHoJU8IIYQQYgPRS54QQgghxAailzwhhBBCiA1EL3lCCCGEEBuIXvKEEEIIITYQveQJIYQQQmwgeskTQgghhNhA9JInhBBCCLGB6CVPCCGEEGID0UueEEIIIcQGopc8IYQQQogNRC95QgghhBAbiF7yhBBCCCE2EL3kCSGEEEJsIHrJE0IIIYTYQPSSJ4QQQgixgeglTwghhBBiA9FLnhBCCCHEBqKXPCGEEEKIDUQveUIIIYQQG4he8oQQQgghNhC95AkhhBBCbCB6yRNCCCGE2ED0kieEEEIIsYE47/2LHoMQQgghhLhkpOQJIYQQQmwgeskTQgghhNhA9JInhBBCCLGB6CVPCCGEEGID0UueEEIIIcQGopc8IYQQQogNRC95QgghhBAbSOdFD0AI59w+gHcBHAB4A8CH3vv7tu1d6/YTALcBvOe9v5vt+z6ABwDuALjvvf/5KedZua8QF+G0Z9q2L30W9UyLq0j2TAPANwH8gM+bfk9fYbz3+tHPC/0BcC/7fAeAB7Bvy+/bsgfwGYA7Wd8PAbyeLX98yjlW7qsf/Vz054xneumzqGdaP1f1B8AH2ec7AB7p9/TV/9F0rbgKvOucewsAvPcPbN0daw+89w7ALe/917PtQPhlkP+V94DHaeA8fYW4KKc906c9i3qmxZXDOXcH4eUNQHymHyApe/o9fUXRS564Crzh0/Qs/yPMf0nAe3+QL9s//GKdLX+7Pvh5+gpxSTQ+06c9i3qmxRVmH8C9hvWv5Av6PX310EueeOFUf/W9B+Bu/svCOfeuc+5t59w959zrtnq/4VBfIqklOefpK8SFOeWZPu1Z1DMtriSmrr1RrX4dwMdc0O/pq4kSL8SVwNSOtxFiOe5mm+5n/2F+5Jz7zDn3BkJw76qcp68Ql8KSZ/q0Z1HPtLiy+DI56F2E381MJtLv6SuKXvLElcB+QfyNc+6Oc+4zhOmug0oRAYJ8/x0AD89x+PP0FeJSaHqmcfqzqGdaXHksy/Yd732cRtXv6auLpmvFc8M597H9hfeZ/afH9VGmt18WBwC+b/85PqoO8wDA161PLe+/giqWzzhPXyFW5rzPNE5/FvVMixfOsmc64x6Ad7L++j19hZGSJ54b+V9+xIJtPwbgqk38x363Yf1n3vv7zrla3t9HSMGvz7tyXyHOw3mf6dOeRT3T4irQ9EwT87C7x5hpi707gH5PX1mk5IkXzQMs/oK4g/Cf3gNkf9mZOnLHe/9DW3U/C/CFbWNG4+vVtqV9hbhklj7T9vm0Z1HPtLiSOOfeBvBzAA+dc/v27L2p39NXG+eD4aAQLwxTPvgX4RsIBpgf2ba6csA9xn/Ytu8D+BmCA/uPfXJgv4egnLx3Vl8hLpsVnullz62eaXHlqH3yMr5tCpx+T19R9JInhBBCCLGBaLpWCCGEEGID0UueEEIIIcQGopc8IYQQQogNRC95QgghhBAbiF7yhBBCCCE2EL3kCSGEEEJsIHrJE0IIIYTYQPSSJ4QQQgixgeglTwghhBBiA9FLnhBCCCHEBqKXPCGEEEKIDaTzogcghBBCrItz7mP7+LH3/m9e6GAugHPuHoDXAcB7/+0XPByxIeglT4gNwDm3D+CO9/7nL3osQjRhz+ht7/2DSz70gff+nUs+5nPHe38XAJxzH77osYjNQdO1QjxHnHP3nHOfOee8c+4D59xb2bZ3nXMf27aPnXPvrnjMOwA+BPAP1fq3nXOfXe4VXD7OuX9Y9Vov4VzP/Z5U3/mHzrn3n+f5T+N53Q97wfsRgM/seRVCPAek5AnxHPHe37X/VO9579+rtv3QOfcpwsvae6sqHtbv2845X216AOCjyxj3M+bHAD69zAM6515fomo+93tSfedXTXFqvB+n3L+18N4fAHjHOffoso4phDgbveQJcbU4sPbhRQ9k/0lf+enbZxRH9RYarv1luSfPi1PuR+P9uwQu/FwLIVZH07VCiE3kuy96AC85un9CbAB6yRNiQ3HOvZ7H/F1FbIxvX9Y4nXP7eZbikvNd6XvyPKnvxwr3b9859/pZ64QQVwNN1wrxkuCceztb/CaAHy+Lm7JA9+8DeBuAq7a9BeAOwtTZbQAPvfcf2bbXEQLk7wD4GsK0Hc/3s1P6fQfAtwF84L2/f9Z5Mu7YOD8FcH/Va7XzfxcAkwY+te3fyfrcs4+fWbzjwj2xJIjvA9gH8E52fXcQYiMfAPgez32e72BdVvx+Xvfeu2w8xfeTHet9pBAAIE2XfhfA97D4jJx2/+4AuIfydDrzUAAABjFJREFU/u3bunedc280fEdvIdxDZG3TNV/qfW043g8sLnDVZ3fZ88Vj857u2753s+Ov++wKcfl47/WjH/08xx8A7wLwAN5v+PnAtu037PNute4zBNsULvuGc/lq+W2EBIB83fv5sRH+8/MA3q76fdzQ7xH72XHeXvU81bV9sOq1Irw4fNxwXafei1PuyVsILzJ1v3r8Z34HZ3znj1bot8r3s2/fz1tnXNf7+bH4jPG7O89zc9Y2ew5er+5p/R29ZeO+U61f+75a3w/PeBZebxjLac/u0ufLttX3/QPYv9dLenY/XHat+tHPeX9e+AD0o5/r9nPaf/hIL1j72br9pv4ICsr72bJv6OOr5UeoXiCb1i851usNY/P5f+7nPU92Pz5Y9VrtP83XG4779mnjP+Pamo5Zv1id+R2s852v+/00jLe4BlvO78myl9kzn5vz3r+m+5ldR/5Cc6H7an3rl7wPGtY1XuOSMS59vux7fL/a9rZdx2U9u3rJ08+l/Wi6Voirz5tAnMbL+RLlVNyp2P4PvU0rVTxAeAlYai/ivf+5c47juZ+vv8TznHqtNmXYZPr8rYZ15+EDAO/ZD8//k1XHdYHzwo67b+c4z307y2LnAcJ0L9nHJYz1LPgd4ezxAc/mvt4FcNumRd9EmKI+cM7t1/e24dk99fmy7+mXzrlvI3hT3vdpKv2t067lGT67QixFL3lCXH3uIPznf79aXy+vcpxl/3E+tO1n8WCFfhc5z6nXuixp4hL+k/whgF/CXvIQ/jPOx3BZ38EyGAu38n1b8jKYcw/hheeHtvxdAD9Yd4Dn4A6w0vjY97Lv622EF7AfeO9/CADOuR/Z+rPG1Phs8vny3h84576GoOi9B+AD59xHPvgfvqhnV4ilKLtWiKtPrchc5Dj7S7bdxmrKyyoKzUXOc9a1nrsklmWQLhsPgPhC8qll+u43nOeyvoNl8JwX/X5yHgK4a5VU3kV46Tm3EfQq96/iPON8Fvf1Y4Tp/yaT57Ou49Sxm0n0gff+b7z3bwC4BWDf7u+lP7tCXBS95Alx9fkUiFl5Bee0A/kUYRqr6T+613GGemLnP2hQKi7zPKdeqw/VPQ6arvuUcll3VlSVOGX7VsM1XtZ3sEA27gt9Pw285b3/uff+h/azrmK0cP8axhjvffYdNdmq1Ptd6n21cd1BOdWen/fN0170Vni+3sqfM7svdwF8Hc/m2RXiQuglT4irBf8DioqA/UfyPYTpt4j9Z7FyObDsP6T6OO8js4DI1udWELD9vnfZ52nY96xrfadhe/0f64PsP86VlCJTft5s6n9Z38ESPkR4eT7v99P0snK76rPuC+hp9+/nyF7q7By1ikV7lqaxxHFf9n3N7lEcix3r53behbi8Bs56vu5W/W8j2Ndc1rMrxKXhvK/LXQohnhXmPUY7ho9QenO9i/CfwFsIqs2HjCmy7W8heHJ9hjAV93PvPf8z5nE/Qoi7+v/bu6OjhIEgjOPfdZCxBOxAxxLswBkqUEqgBuwgtkAJ0Qq0hpTA2MH6cHt4JGHIQAb0+P+eIECSu9zDktvsfXe2Lf1OQtrPnX+nktR2p7ZCXAf31s+zUqz31dhuHbFUY+3NP+vu4+BxsnbfW7aW7762Dnz+6Zta69ciWyhO3X2b2ftAPy2tsz6wX5+6u33see35Ta14TWfeVynISCU5KkmPtltfcLDf9lzr1re9KI6blbe3UswzTDaKwU59qD+G+q/Tt3Nlfa9Yd67yY6c8uNSOFGC1kj78PJb5WDimX7Pfri1bE9j3tVBcEznt46azbczY7Y0vH6tf+s07raTtH4RRbRkxdnfaA5yCIA9ATwjBzCwc/uYkx+oFeTiNB3BL7RbpTU+9rnRkft5fVFpQVFp7cFlM1wK4tLOU9rgyC8U7wdt+NbPW78g9i7VpgatAkAfg7EIIdQhh7W8fFKfRMJ1Gcep/yFz0N3AVqJMHYCvLtZMHYb28tYnUik8qPikmrVMrbEKec7fxhzby6zdTllsJoGzk5AEA/q0QQuMvGzN7vejJnMAf+kllisjJwyQI8gAAAApETh4AAECBCPIAAAAKRJAHAABQIII8AACAAhHkAQAAFIggDwAAoEA/5iS1xahZ3rcAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 720x720 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"fig = plt.figure(figsize=(10,10))\n",
"m = Map(clipped_fits_files[0])\n",
"blc = SkyCoord(Tx=-375*u.arcsec,Ty=-310*u.arcsec,frame=m.coordinate_frame)\n",
"trc = SkyCoord(Tx=-225*u.arcsec,Ty=-160*u.arcsec,frame=m.coordinate_frame)\n",
"m = m.submap(blc,trc)\n",
"ax = fig.gca(projection=m)\n",
"norm = ImageNormalize(vmin=50,vmax=1200,stretch=AsinhStretch())\n",
"im = m.plot(axes=ax,title=False,annotate=False,norm=norm)\n",
"ax.grid(alpha=0)\n",
"lon,lat = ax.coords\n",
"lon.set_axislabel('Helioprojective Longitude [arcsec]', fontsize=20)\n",
"lat.set_axislabel('Helioprojective Latitude [arcsec]', fontsize=20)\n",
"lon.set_ticklabel(fontsize=15)\n",
"lat.set_ticklabel(fontsize=15)\n",
"xpix,ypix = m.world_to_pixel(\n",
" SkyCoord(Tx=-370*u.arcsec,Ty=-170*u.arcsec,frame=m.coordinate_frame)\n",
")\n",
"time_text = ax.text(xpix.value,ypix.value,f'$t={m.meta[\"t_obs\"]:.0f}$ s',color='w',fontsize=20)\n",
"xpix,ypix = m.world_to_pixel(\n",
" SkyCoord(Tx=-240*u.arcsec,Ty=-170*u.arcsec,frame=m.coordinate_frame)\n",
")\n",
"ax.text(xpix.value,ypix.value,f'${m.meta[\"wavelnth\"]:.0f}\\,\\,\\mathrm{{\\AA}}$',color='w',fontsize=20)"
]
},
{
"cell_type": "code",
"execution_count": 95,
"metadata": {},
"outputs": [],
"source": [
"def update_fig(i):\n",
" m = Map(clipped_fits_files[i]).submap(blc,trc)\n",
" im.set_data(m.data)\n",
" time_text.set_text(f'$t={m.meta[\"t_obs\"]:.0f}$ s')\n",
" return im,time_text"
]
},
{
"cell_type": "code",
"execution_count": 96,
"metadata": {},
"outputs": [],
"source": [
"interval = 30/len(clipped_fits_files)*1e3"
]
},
{
"cell_type": "code",
"execution_count": 97,
"metadata": {},
"outputs": [],
"source": [
"ani = matplotlib.animation.FuncAnimation(\n",
" fig,\n",
" update_fig,\n",
" frames=len(clipped_fits_files),\n",
" interval=interval,\n",
" blit=True,\n",
" repeat=True,\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 98,
"metadata": {},
"outputs": [],
"source": [
"ani.save(f'noaa11158_aia{m.meta[\"wavelnth\"]:.0f}.mp4',dpi=85)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python [conda env:synthesized-timelags]",
"language": "python",
"name": "conda-env-synthesized-timelags-py"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.6.6"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment