project hrap grid.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Script to reproject the HRAP gridcell centers\n",
    "\n",
    "**geographic projection for HRAP grid:**\n",
    "\n",
    "'+proj=longlat +a=6371200 +b=6371200 +no_defs' \n",
    "\n",
    "**WGS 84 projection**\n",
    "\n",
    "'+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs'\n",
    "\n",
    "** Stereographic projection for the data from online**\n",
    "\n",
    "'+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-105 +k=1 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +to_meter=4762.5 +no_defs'\n",
    "\n",
    "I think this compresses the data too much, but if I remove the to_meter part of the argument then the dataset is spread out but the grid cells are not regular...\n",
    "\n",
    "\n",
    "Theodore Barnhart | tbarnhart@usgs.gov | August 28, 2017"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from pyproj import Proj,transform"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "header = ['cell','num','lat','lon']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "dat = pd.read_table('./data/hrap_lat_lon.txt',delim_whitespace=True,header=None, names=header)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "lats = dat.lat.as_matrix()\n",
    "lons = dat.lon.as_matrix()\n",
    "cells = dat.cell.as_matrix()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7fc20123bcc0>"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4cAAAJlCAYAAABpFoerAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xu8ZFdd5/3vr08IuQJt4iODTBJiY0CRUbsU5uU4plHi\nDEQugj4oykUwKs3zoBEvGQeIXAeRIer048jDBFDxwqAyDOgQMB11VC6nZLgIRtrAGblISDgiSUhC\n0r/5Y1cllZM651Sd2muvqu/+vF+v8+ruOrtq79P9y8pvrfVba0VmCgAAAADQb/tqPwAAAAAAoD46\nhwAAAAAAOocAAAAAADqHAAAAAADROQQAAAAAiM4hAAAAAEB0DgEAFUTExyPiqhmvfVpEZEScP8O1\nr4sImzOa5vl7chIR54/+zZ9W8RnOGT3DpbWeAQC6RucQABY0kcg+d4drMiLeuuW1q0avj7++FBGf\niojfjYiH7HLP/RFx8+h9P7DDdR/fco9bR6+9JiL++Rw/40kR8ayIuDIiPjt61n+MiPdGxMsj4kGz\nfhbuKiIujYjH1X6OPhp1AC+NiK+v/SwAsAxOqP0AANBzt0h65uj3J0s6KOnpkh4VEYPMvHqb9z1Z\n0omSPibpGZJ+c4d7fELSJaPfny7pfEk/NLrHQzPzup0eMCLOlfRWSQ+W9CeSXiXp05JOk/T1o896\nbkSclZmf3OmzJpwnqcQM3w9L+tECn1vSCyS9XtKbp3yv1N/TsvtTNf89fKnwfc5R8/f/cUn/a8v3\nNkbPcFvhZwCApUHnEADqui0zJzt2/39EfFjSL0l6tqT/Z5v3PUPSUUn/TdJlEfFVmfl321z7+S33\n+NWIuHb0+U+X9IrtHi4iTpb0NklfJem7M/MPplxzkqSf0C6dmIi4h6S1zLw5M2/Z6dq9yswvqXyH\nojOl/p66FBGnZ+YX5nlPZh6XdHOhR5r1GbL2MwBA1ygrBYDl88ejXx847ZsR8Y1qZuxeL+kNajpD\nT5/zHm8f/Xpgl+ueKelBkl4xrWMoSaPO3ssy81MTz3jpqIz1ayPiP0bEJ9Qk2g8ffX/qWrqIeGZE\n/E1E3BIRxyLiOZJi1h9q2prD8WsRce+I+NWIuHZUkvvnEfGwGT93/POcFxEvjYhPjJ7x/RHxqCnX\nPysiroiIT45KeT8dEb8ZEedMXHPOxLM+dbL8d+Kau/w9RcS7I+IzEXG3wd2I+M7R+3984rWIiB+L\niGFE3BQRX4iIoxFxaMafe7ze8ztGfwcbo5/7AxHxpCnXf3xULv0NEfH2iPi8pA9MfP/MiDgSEX8/\n+nv5+9Gfz9jyOVPXHM7780TEE0bf/8fR9VdHxC9HxImjzz46uvS1E3//V43eO3XNYUScEBE/ExEf\nHsXR9RHxBxHxdVuuu+P9EXFhNCXYN49i4RXT/g0BoDYaJgBozykRcWYLn/NVo18/t833nyHpRkm/\nl5k3RsTb1HQunj+acZnFuOO5Y0mppCeOfn3NjJ+71RskfVHSK9XMLH56uwtHnZpXSXq/pH8n6RRJ\nPyXp2j3ee6u3S/qspBdKOkPSxZL+MCLOmWNm6/VqOuO/qKas98clvTkivjozPz5x3XMlvUvSL6v5\nd3yImo72IyLi6zLz+tGz/KCk35D0Z5JePeP9j0j6N2pKfSc9RU0J5G9NvPYbkr5P0pskvVbSPdWU\nJL8jIr47M98y24+tl0s6VdKvqvl3fLqk346IkzLzdVuuPUvSlZL+q6TfU1N+rIi4t6S/UDMgcbmk\nv5L0DZJ+TM3fyzfP8O8w888TES9RE0cf1p2l0F8l6QmSnq+mdPWlo2terebfQJI+s8szvEHS90p6\nx+jv476SDkv6y4j41sx835brHyXpWZL+8+jnfqya+Ngc3R8Alkdm8sUXX3zxtcCXmjV8OcPXW7e8\n7ypJN0g6c/T1zyU9Ts36p5T0qCn3OklNZ+N1E689dnT9v51y/cclfWTiHg9Qk9j/o5pOzkN2+dmu\nV1OWuvX1tYnPHH+dPPH9S0fPdJWkE7Z5rqsm/nwfNR3eD0s6ZeL1+4/+jlLS+TP8W7xOo4rAra9J\n+v+2vP49o9d/ZIbPHf88b5UUE69/0+j1l225/tQpn/Hto2t/esvrOfnvucvf05epWaf6xi3XnT76\n+3vLxGuPH332RVuuPUHSupr1qjHtvhPXPm30GRuS7j3x+r1Hr31uy7/7OHafOeWzXjL63rO2vH54\n9PqLpvw39bS9/DySvnl07ZWSTtpyfUxcd7f7TFx3zuh7l0689sjRa7+7JQ4eqqZj/mdT3n+jpHO2\n3P9Dkj69W9zxxRdffHX9RVkpALTn1WqSx2lf2zlVzQzSZyX9b0l/oGZG6qmZ+YdTrv9uSfvVzCCN\nvU3N7NoPbXOPB03c4xo1sxfXSXpsZn5ol5/pXpL+acrrD574zPHX4SnXXZaZs2zocYGamcIjmXnT\n+MXM/ISamZo2vGrLn68c/Tq1fHcbv5SZd5R9ZuZ7JX1h62dk5o2SFBH7RuWsZ6qZEf28pJlKWafJ\nzM9J+u+SHhMR95n41hPV/P1NxsUPjJ7tzaNyzjNHz3Gf0Wecs/W5d/Crmfn5ief4vJqZsP1qOliT\nPqdmVm+rx6uJk60zpL+mJh4fv8szzPPzPHn06yWZeZd1gzmyy722M37Gl2yJgw+oGTj4VxHx5Vve\n8+acmFUeve+opPtGxGl7fA4AKIKyUgBoz0cz853TvhGx7bK5myV91+j3X6amNPCR2n5N+DPUJNif\niIjJ9YLvkPQ9EXFm3n330Y+r2cVTkm6V9KnMPLbDzzHpn9R0ELf6mO7s9P4LNWWW0/ztjPc5d/Tr\n30z53odn/IzdXDP5h8y8fvTvcsb0y3f/jJHPbf2MiHiEmtLFh6mZ7Z20f477TfPrakojv1d3drSe\noqZMcbLU9MFqZhR3KpP8Cs32b/SRKa+N/13O3fL632Xm7VOuf4Ck9a2DBZl5W0RcLekbd3mGeX6e\nB6qZtXv/Lp85rwdIOq7pfx8fUjOL/wA1/42OTYuZ60e/nqFmZhwAlgKdQwCo6/bJDmVEvElNgv/q\niPir0YzE+HsPkHRITVnadgn9D0i6bMtrN27XaZ3BhyT964h4QGZ+bPziaGbsnaPn2mlm8KYdvjdp\n3HueNqMz84Y0O9mmwzLv5+/6GRHxTZKukHRM0s+q6Uh/Uc3P9jtafDO4P1TT+XiKmjg5S9K3SfrP\nedfdTWN03ffv8Fm7zRyPzfPvMuu/+bzm+XlCZY4A2Ussbhcze/08ACiGziEALJHMPB7NDp0fVjMb\nd8HEt5+uJpn8YTVrBrd6sZqZxa2dw0W8SdK/VrOZys+1+LlbjY/heLDuLPfUxGur5PvVrMn8t5Md\n6og4VYvPGo5n2n5L0nOiOYPy+9TExeu3XPpRSV8t6V2Zuejs1NdI2rp5zfjfZdrM2DTXSDovIk6Y\nnD0c7dr51TN8zjw/z9VqNu15qKT37HDdvB3Iv5P0nWp+9g9s+d7XjH79mABgRbHmEACWTGZ+VM2O\nk4+MiH8lNWvX1GwO8sHMfE1mvmnrl6TflvSQ0cxVW16jptTzpyJiuzVhbcx+vEPN7NrhiDjljg+O\nuL92nilaRuOZoq1/L/9O0/+/e4OakuJ5jDuCT1Gz4+nVmfnuLdf8+uh+L5v2ARHxFXPc78dGu42O\n33tvST+qZpDiT2b8jDdL+nI1Aw2Tfnj0+tSjUibM8/OMd2x9aUTcc8q143+bcSdz1r//N49+vWTi\nMxQRD5H0GEn/MzM/O/WdALACmDkEgOX0UjUloj+vZpfLC9TsZvpfdnjP76nZVfMZkt7bxkNk5hcj\n4tFqSl1/f3QG3BWS/kHNWsQHSfq/1XSI/n6B+2xGxPPUzJb+RUT8upoNVn5UzYzRNyzyc3TsDyT9\nhJpjMl6tZp3nI9XMYk07OuRdkr4jIn5GzaZEmZm/s9MNMvN9EfHB0X3upabjufWaN0XEayU9O5qz\nMd86uv/9Jf1LNUdKbF0vuJ3rJL07Ii5X0+l9upojK545uYHQLn5BzQ6xR0bP8z41/67PUDPT9ws7\nvXmenycz3xMRL5f0M5KGEfG7amL2AWo27/lmNR3bD6vZ5OZZEXHT6LVrM3Pr7PX4Gd4REW+U9CRJ\n+yPirbrzKIubJf2/M/5dAMBSonMIAEsoM68eJ6ER8W1qEmhJ+v0d3vOhiPjb0Xt+IjO/2NKzXBMR\nB9XshvpEST+p5iiDG9Wsq3uNpP+SmVcveJ9XRsQNas4ffJmazuYvqtnh8/JFPrtLmfnnEfEESc+T\n9CI1M6LvVLMu8E+nvOVZas4u/Dk1G65IzdrE3bxezd/PcUm/uc2z/FBEHJV0kaRL1OyE+w9qzhi8\nZMYfSWo6Wd8q6dlqNn35qKQnZ+Zv7fiuuz7L5yPiW9QMeDxGTQfzM2p2PX1BznDW5Dw/T2b+bES8\nf/TMP61m1vHv1azZvGl0zRcj4klqSrIvU3Nu4p/o7qXNk548ut/T1JzfeePoPc/LzA/u+hcBAEts\nfM4PAADAXUTE09QcS3EoM6/q8L7frqZD/YOZObXjCwBoH2sOAQDAsrnf6Ndrqz4FAPRMlc5hRKxF\nxPtGtfqKiNdFxMci4n+Nvr5+m/c9NSI+Ovp6ardPDQAASoqIr4iIH1VTYvsFSX9Z+ZEAoFdqrTl8\njpoDZCcPVv6p0W57U0XEl0l6gaSBmq2nhxHxlszcLPqkAACgKw+W9Co1G8V81yzrEAEA7el85nC0\nLfmj1WxgMI/vlPSOzPzcqEP4DjVnGAEAgAIy83WZGV2tN8zMqzLz5Mw8mJmzHpEBAGhJjZnDy9Ts\nGnb6ltdfEhHPl/THkn42M2/Z8v2v1F23Sf/E6LW7iIiL1OxiplNPPfXggx70oLaee6qbbr1dN9xy\nm0675wk65cS1ovdCP137hVv0mX+6+Y4/f8W9TtL/dfrdju0CVg6xDVfENlwR26tpOBxel5lfPsu1\nnXYOI+JCNecHDSPi/IlvXaJmK+oTJb1azZbZL9z69ikfebetVjPz1aPP0GAwyPX19RaeHKhnuLGp\nJ7/mXfrSbcd1jxP26Q3PfLgOnr2/9mPBzHBjU++65no9/NwzOosvYhtd6Tq+iW24IrZXU0RszHxt\nl0dZRMTLJP2gpNsknaRmzeHvZ+YPTFxzvqTnZuaFW977fZLOz8wfGf351yRdlZm/vd396ByibTUS\n6Jr3RT+M/2d/623HdWLH/7MntlFarfgmttGFWgN7xPZqiYhhZg5mubbTmcPMvESjQ2onOoE/EBH/\nLDM/HREh6XGSPjTl7W+X9NKIGEfhBZrvAF9gITUT6INn76cBRjHvuuZ63XrbcR1P6Uu3Hde7rrme\n2IaNWvFNbKO0WnkJse1tWc45fENEfFDSByWdKenFkhQRg4h4jSRl5uckvUjSe0dfLxy9BnRiWoIB\nOHj4uWfoxBP2aS2ke5ywTw8/94zajwS0hviGK/ISlFDrKAuNdj67avT7R2xzzbqkZ078+XJJl3fw\neMDdjBOMcZ09CQZK6bpk5+DZ+/WGZz6cMiEUV6McjfhGF2rENnkJSuh0zWHXWHPoi7V/cFWzfBko\nidiGK9ZtY9kt7ZpDoA2s/YOzmuv/gJKIbbhi3TacLMuaQ2Bm1NjDGeuj0IXhxqaOHD2m4cZmZ/ck\ntuGK2IYTZg6xcqixR1dYHwVHNXc4JLZRGu02sBg6h1gYm2fAEeXLcEUJHFzRbsNZV/k2nUMshDN2\n4Ir1UehK1wNsVF/AFe02XHWZb9M5xEJoiNEFtgiHqxoDbFRfoCsMfADt6DLfpnNohAQajlgfBWe1\nBtiovkBpDHzAlXu+TefQBAk0XLE+Cs4YYIMrBj7gqA/5Np1DEyTQcEXyjK6wyyFcuc90AF3pQ75N\n59AEjTC6QPIMV+xyCFd9mOlAf7GutX10DgvheAe4IXmGMzbXgqs+zHSgn1jXWgadwwI43gGOSJ7R\nFUrggPYQ23DFutYy6BwWQBKNLlBKAUeUwMEZpflwxaCeD/vOIcEKR5RSwBUlcHBFaT5cMajnxbpz\neNOttxOssEQpBVwxuAZXVBXBFYN6Xqw7hzfccptOIFhhiAQaXaAEDs4ozQfaQWx7icys/QzFfM1D\nvyGPP/aldwRrlyUc6I8aCXTN+6IfapbAAaXVim/abZRGToJpImKYmYNZrrWeOTzlxDX9GiPQKIg1\nJHBFCRycUZoPR+QkaMO+2g9Q2sGz9+vwoQMELIqYlmAADsZlQmshyoRQ1HBjU0eOHtNwY7OzexLf\ncEROgjZYzxyif1hDAkes/YMrdjmEM3ISrCI6h7DB8Q5wRJkQnLHLIVyRk2BV0TmEDdaQwBFr/+CM\nmQ64IifBqqJzCBskGegCZUJwRfkyXNWIbdpurCrroywGg0Gur6/XfoxeYitlOGL7e7ji6BK4qhnb\ntN1YFhxlgapYIwVXlAnBFeXLcMW6VjgrMQBB5xCtI8mAK8qE0AVK4ID2ENtwVWoyhs4hWkdDjC6w\nPgqOONoBzmi3gfaUmoyhc9gDXTfGNMQojdJluKIEDq5ot+HMqeKDzqG5mqPQNMQohdJluKLyAq5o\nt+HKreKDzqE5GmN0geMd4IrKC7ii3Qba4VbxQefQHI0xSqsxYkYCjS5QeQFXtNtw5VTeWQudww6x\nEBuOON4Brqi8gCvabThyK++shc5hR1iIDVduI2bAGLENV8Q2HLmVd9ZC57AjjECjC8xOwxWxDVfE\nNtAOBj3aQeewIwQsSmN2Gq6IbbgituGMDb1WUy87h4zSwRGz03BFbMMVsQ1XbOi1unrXOWSUDq6Y\nnUYX2AkOrohtoD0MfKyu3nUOCVZ0hXIKuGEnOLgituGMgQ/Mo3edQ4IVXaCcAo7YCQ6uiG24YuAD\n8+pd55BgRReYoYYjBtfgitiGKwY+MK/edQ4lghXlkWigC5QuwxWxDUeUd2IVRGbWfoZiBoNBrq+v\n134MVFajMa55X/RDzc21gJKIbTiqGdfkI4iIYWYOZrm2lzOH6A92p4UrSpfhitiGI8o7sSr21X4A\noKRpjTHQtuHGpo4cPabhxmZn9xyXCq2FKBVCMcQ2XHUd28Q1VgUzh+gMtfZwxE5wcEVsw1WN2Cau\n0ZVF8206h+gESQZcUSoEV8Q2XNWKbeIapbWRb1NWik7ULO88ePZ+HT50gAYZRVAqBFfENlwR23DV\nRr7NzCE6QXknusIW+HBUoyyf2EYXiG2gPW3k2xxl0VM1GmO2UkZpbIEPR8Q1XBHbcLZMR6lxlAV2\nVHP9H40+SmILfDgiruGK2IarVT5KjTWHPcTxDugCW+DDEXENZxzvALRjlXNtZg57iPV/KI3daeGI\nuIYzjneAK45Smw+dw8pYiA1HbIEPR8Q1nHG8AxwxqDc/OocVrXI9MrCTVR4xA7ZDXMMZ8Q1HDOrN\nj85hRSzERheYnYYrji2BK9ptOKK8czXQOZzQddASsCiN2Wm4YtdluKLdhiPKO1cHncMRFmLDEbPT\ncEVswxWxDUeUd64OOocjLMSGI2an0QVKheCK2AbaQVyvjsjM2s9QzGAwyPX19ZmuHc8cjoO2yzIO\n9EeNRKPGPdEfNUvgiG2URGzDGflIv0TEMDMHs1zLzOEIJZ4ojTVScESpEFwR23BFPoKd7Kv9AMvk\n4Nn7dfjQAQIXRUxLNIBVNy4VWgtRKgQrxDZckY9gJ8wcopdYRwJXHO8ARxztAFfkI1g2rDlE77CO\nBK5qxjZQCnENV+Qj6AprDoEdsI4ErtgCH46Ia7giH8EyYs0heod1JOjCcGNTR44e03Bjs7N7Etvo\nQtexTVzDFbGNZURZKapjO2W4oVQIrmrFNnGN0mrFGLGNLlBWipXBdspwRKkQXNWKbeIaJdUc0CO2\nUdpwY1Nrp51x31mvp6wUVbGdMhxRKoQuULoMtINcBK7GAx9rp+7/ylnfw8whqmI7ZXSB4x3gpmbV\nBbGN0rpus8lF4Go88KGY/T10DnEHzpGCI0qX4YjSZbiq0WaTi6ALNc+01BybzNA5hCTq7eGLbfDh\niJkOuGJdKxzVrvZ42GX/+KlZ30PnEJJIoOGLJBpdoHQZjmrOdNBmw0ntao/bb7j+H2a9ns4hJNEY\noxuULsMRpctwVHumgzYbJbGudXt0DpcUo9BwQ+kyXFF5AUe1Zzr4bwilsK51Z3QOlxCj0HBEAo0u\nUAYHV8x0AO1gXevO6BwuIZJolEYCDUeUwcEVMx1wRT6yfOgcLiGCFiWRQMMVZXBwxUwHHJGPLCc6\nh7tgAw24IYGGKwbW4IrYhiPykeVE53AHbKABRyQZ6Aoba8ERg8ZAO8hHlhOdwx2w9g9dIIGGIzbW\ngiMGjeGKQQ+M0TncASMaKI0EGq4YXIMj4hqOGPTAJDqHO2BEA6WRaKAL7AYHR8Q10A5yEUyic7gL\nRjRQEokGSmM3ODgiruGMMy1RE51DYALr/+CG3eDgiLiGK860RG10DoER1v/BESPCcERcwxVnWqI2\nOofACDX3KI3d4OCKqgs4Yl0r+ojOIZYSDTLcsBscXFF1AUesa0Vf0TnE0qFBhiNmpuGK2IYj1rWi\nr+gcYunQIMMRM9NwRWzDEXGNvqJziKVDg4wusEYKjljXCle02UA3IjNrP0Mxg8Eg19fXaz/GSquR\naNS8L/qh5vo/oBTiGq6IbbjqKt+NiGFmDma5lplDbIsNNOCKNVJwRFzDFbENR8s66LGv9gNgeU1r\njIG2DTc2deToMQ03Nju757h0eS1E6TKK6Tq2iWt0gTYbaMey5tnMHK6QrkstWfuH0tiZFq5qxDZx\njdJos+GKI9TuROdwRZBowBE708JVrdgmrlESbTYcMehxV3QOVwSJBhwt66gZsChiG46Iazhi0OOu\n6ByuCBpkdIGtwuGI4x3giLgG2kGOfVccZbEHHO8AR8u6axawCOIajohrOKuR77rn2BxlURDHO8AV\nW4XDEXENR8Q1XNVc/8d/Q40qR1lExFpEvC8i3rrl9V+JiBu2ec89IuL1EfHBiPhIRFzSzdPe1bJu\nOwssiq3CURpb4MMVx5YA7SDPrq/WzOFzJH1E0r3GL0TEQNJ9dnjP90i6Z2Z+XUScIunDEfHbmfnx\nok+6BXXJ6AJrSeCG3eDgit3E4YrjHfqp885hRNxf0qMlvUTSxaPX1iS9QtL3S3r8Nm9NSadGxAmS\nTpZ0q6R/Kv7AW9AgozRKl+GI3eDgit3E4YgBvf6qMXN4maSflnT6xGvPlvSWzPx0RGz3vjdJeqyk\nT0s6RdJPZObntl4UERdJukiSzjrrrBYf+040yCiJtSRwxGgwXBHbcMSAXn912jmMiAslXZuZw4g4\nf/Ta/dSUjJ6/y9u/WdLtku4nab+kP4uId2bmNZMXZearJb1aanYrbfUHQO9QUgFXHFsCR5TkwxG5\nCLrU6VEWEfEyST8o6TZJJ6lZc3jL6Ovm0WVnSbomMw9see8RSe/KzN8Y/flySf8jM9+43f1KHWWB\nfqhZ3um+pTLqYht8OCKu4YhcBG2Y5yiLTncrzcxLMvP+mXmOpCdJujIz92fmfTPznNHrN23tGI78\nb0mPiMapkh4u6W86e3j0Ts0dsw6evV+HDx2gMUYR7AYHR8Q1HJGLoGtVjrKYVUQ8JiJeOPrjEUmn\nSfqQpPdKem1mfqDaw8EeW4WjCxzvAFcc7wBHxDXcdVpW2jXKSr3UKG+gpAIlUS4EV7Vim7hGScQ1\nVtU8ZaW1zjkE5lJzS2UaYpTCbnBwxfEOcERcow+WuqwUGGMtCRxRLoQuULoMtIO4Rh8wc4iVwJbK\n6ALHO8ANB1nDFceWAGXQOcSekETDDaXLcETpMhzVXK9NXKO02mtM6RxibiTRcFQziQZKoeoCjmiv\n4WoZzmulc4i50SjDEUk0ukDVBRx1Hde013C1DDk2nUPMjUYZXSCJhhuqLuCoRlzTXqMLNco7lyHH\npnO44liQDUck0XC0DCPCQNs43gGO+ryZF53DFcaCbLgiiUZpfR0RBtpGXMNRnzfzonO4wkig0QWS\naLjp84gwvFFNBEfkId2ic9giFmTDDUk0HPV5RBi+qCaCI/KQ7tE5bAkLsuGIJBqOGFiDI6qJ4Ig8\npHt0DlvCgmw4IolGaZTBwRFlcEA7iOvuRWbWfoZiBoNBrq+vd3Kv8czhOHhrHFoJbzWSjZr3hb9l\nOOwXaFvNuKa9Rmk1Yoy4XlxEDDNzMMu11jOHN916e2f3YiQaJbGWBI4og4MjyuDgimOm+mFf7Qco\n6ZrrbtBwY7Oz+x08e78OHzpAAKN105INoG3DjU0dOXqss3ZzXC60FqJcCMUQ10A7yEX6wXrmMFOM\nRKMIdqaFGzbVgiPiGq5Y14pSrDuHESJw0TqSDThiUy04Iq7hiOMdUJJ15/DcM08jcNE6kg04YkQY\njohrOGJdK0qy7hyecuJa7UeAIZINdKHrkiFGhFEax5YA7SAPQUkcZYGVx7bKcMMRD3BDTMMVx0xh\nFXCUBXqDbZXhiCMe4IaYhiOOmYIj66Ms4I9tleGIrfBRGsc7AIsjB4EjZg6x0qi7R2msk4IbdlyG\nK46ZAhZH5xCtIYmGG0qG4Igdl+GIQQ+46jq/pnOIVpBEwxHrpOCI2Q44YtADjmrk13QO0QqSaJRW\nY2aaJBpd4NgSOKLEE1hcjfyazqEhkmi4qbkrLUk0SmLHZTiixBOu+jDoQefQDEk0HNWcmSaJRklU\nXcARJZ5w1JdBDzqHZkii4YiZaXSBqgs4Iq6BdvRl0IPOoRkaZHSBNVJwQ9UFHBHXQHv6kmPTOSyM\nJBpuWCMFR1RdwBFxDVccn1YOncOCSKLhiDVScNSXEWH0C3ENRxyfVhadw4JIouGIZAOlMSIMR8Q1\n0A7y67LoHBZEEo0uULoMJ4wIwxFxDVdseOSnV51Dkmi4oXQZbhgRhiPiGo7Y8MhTbzqHJNFwRMKB\n0vpw4C/6h7gGFseGR5560zkkiUZplFbATV8O/EW/ENdwxaAH2tCbziEBjJIorYCjvhz4i34hruGI\nQQ+0pTedQwIYJVFaAUcMqsERcQ1HDHqgLb3pHEoEMMoh2UAX2FQLbjjeAY5YZoJVFplZ+xmKGQwG\nub6+Xvsvvln+AAAgAElEQVQx0LEajXLN+6Ifam6FD5RATMNRzbgmD8F2ImKYmYNZru3VzCH8cZYU\nXLGpFtwQ03DEMhOsun21HwBo07RGGXAwLhlaC1EyBAvENBwR11h1zBzCCjX36ALrpOCINa1wRFwD\n82HNIYqqkURTc4+SWCcFR8Q1HBHXcLSXPJc1h1gKNc/+o/FHKayTgiPiGo6Ia7jpIrdmzSGKYf0f\nShtubOrI0WMabmx2dk/Wk6A04hqOiGtgcV3k1swc9kjX5Zas/0NJNWemWU+CUohrOCKu4coxt6Zz\n2BM1GmYaZZTEduFwRFzDEXENR665NZ3DnqjVMNMooxRmpuGIuIYj4hqOXHNrOoc9QcOM0tguHI6I\nazgiroHFuebWHGVRQa2jFjjiAaWwXTgcEddwRFzDFcenbY+jLJZYzUaZEk+UwnbhcERcwxFxDUcc\nn9YejrLoGMc7oAtdbxnOduEojW3w4Yi4BtpBft0eZg475lqfjOXhunsW+ott8OGIuIarGqWW5Nft\n6X3nkEXZcOO6exb6i23w4Yi4hiMGPVZfrzuH1CfDEaNncENMwxFxDUcMeqy+XncOWZSN0mqUVjB6\nhtKouIAj4hpuKO/EXvS6c0gAoyR2poUjKi7giLiGG8o7sVe97hwSwCiJmWk4Iq7hiLiGG8o7sVe9\n7hxKBDDKYWYapVEyBEfENbA4Yhp7FZlZ+xmKGQwGub6+XvsxsCRqJBw17ol+qFm2TFyjFOIarshB\nUFNEDDNzMMu1vZ85RD+wngRuKBmCI+IajshBsEr21X4AoAvTEg5glY1LhtZClAyhmOHGpo4cPabh\nxmYn9yOu4YgcBKuEmUP0ArX3KI1t8OGmxmwHcY3SWNMK7IzOIaogkYYTSobgqFaJJ3GNUjjeAa7a\nzKvpHKJzJNJwwzb4cMRsB9ywphWO2s6r6RyicyTSKK3rmWmSaJRWoxSO2Q6URlsNLK7tvJrOYc9R\new83rJOCm5rHOzDbgVJoq+HIIa+mc9hj1N7DEeuk4IZqCziirYYbl7yazmGPUXsPR8xMozRK4YDF\nEddw45JX0znsMRpmlMY6KbihFA6u2EUcbhjI25vIzNrPUMxgMMj19fXajzGTGkl0zfvCX811UkAp\nR44e0yuvuFrHU1oL6eILztPhQwdqPxawENpruKkV08uaV0fEMDMHs1zLzOESYLMBOGKdFBy5jAwD\nk2iv4YY1rXtH53AJ0CjDEUk0ukApHNw47HYI1EZM7x1lpUtgPHM4DmDKOdA2ypbhiFI4uKkZ07TX\nKKlGfBHTd6KsdEGMRMMJZctwRdUF3LjsdghMqnnEAzE9PzqHWxDAcEMCjS5QCgdH7HYILI48ZLXQ\nOdyCAEZpJBtw43LwLzCJY0vgiIE87IbO4RYEMEoi2YAjSuHgiN0O4YaBPMyCzuEWBDBKItmAIwbV\n4Ii4hhsG8jALOodTEMAohWQDpdUoGWJQDaUR18DiyEEwC46yQK+xtTKccLQDHBHXcMQRU+gSR1kA\nM2BnWrhhQy04Iq7hhiOmsMz21X4AoJZpCQfQpuHGpo4cPabhxmYn9xuXDK2FKBlCMcQ1sBjyDywz\nZg7RW9TeoyR2poUj4hqOOGIKuBOdQywFNhuAG3amhSPiGm4Y8ICrydx6HnQOUR2193DEyDAcEddw\nw4AHHG3NrePEk0+d9b10DlEdmw2gC13PTjMyjNKouIAjSjyBxW3NrfedeMrps76XziHuhoYZbtiZ\nFm6ouIAjSjzhqMZA3tbc+vitN31h1vfSOcRd0DDDEbPTcENMwxElnnBTc3B6MrcevPiLN876XjqH\nuAsaZjhidhqlUXEBN8sw20FcY9XVHMjba25N5xB3QcOM0lgnBTdUXMDNssx2ENdYdauYV9M5XHJs\nogEnrJOCIyou4GYVZzuAWZBX747O4RJjEw24YZ0UHK3iyDCwE2IajsirZ0PncImRSKM01knBDWXL\ncMRsB7A48urZ0DlcYiTSKIl1UnBD2TIcMdsBR2x4tLzoHM6I0Wi4YZ0U3DAqDEfENdyw4dFyo3M4\nA0aj4YgRNLghpuGIuIYbNjxablU6hxGxJmld0icz88KJ139F0tMz87Rt3vdQSb8m6V6Sjkv6psy8\nufTzMmqHLrCmBG6IabihigiO2H8Ak2rNHD5H0kfUdPIkSRExkHSf7d4QESdI+k1JP5iZ74+IMyR9\nqfSDSgQxymNNCdwQ03BDFREcsf8AttrX9Q0j4v6SHi3pNROvrUl6haSf3uGtF0j6QGa+X5Iy8/rM\nvL3ks46Ng/jiC87r9H8G6I9ps9PAKiOm4YaYhqNacX3w7P06fOgAOfUSqjFzeJmaTuDpE689W9Jb\nMvPTEbHd+75aUkbE2yV9uaTfycxf2HpRRFwk6SJJOuuss1p7aEbtUBKz0yiNsiG4Iabhhh08sQwi\nM7u7WcSFkh6Vmc+KiPMlPVdNR+6Nks7PzNsi4oZpaw4j4rmSDkv6Jkk3SfpjSf8+M/94u/sNBoNc\nX18v8JPAWY3GueZ94a9WORwxjVKIabipWbZMXPuLiGFmDma5tuuZw2+R9JiIeJSkk9SsOfxrSbdI\nOjaaNTwlIo5l5oEt7/2EpD/JzOskKSL+UNI3qukkAq1gTQkccWwJ3BDTcMMOnlgWna45zMxLMvP+\nmXmOpCdJujIz92fmfTPznNHrN03pGErS2yU9NCJOGW1O822SPtzZw6MXWFMCR+OyobUQZUNo3XBj\nU0eOHtNwY7OzexLTcENMY1ks9TmHEfEYSYPMfH5mbkbEf5T0Xkkp6Q8z8211nxBuqL1HaWyFDycc\nZg1XHMWDvup0zWHXWHO4+mok0tTeo5SaZctACUeOHtMrr7hax1NaC+niC87T4UPTin+A1UFbDTfL\nvOYQmBnnpMFNzTUlQAlUW8ARbTXcxIknnzrrtXQOsbRonFEaW+HDDaVwcERbDSxm34mnnL77VQ06\nh5gZjTOc1JiZJpFGSVRbwBFtNbC447fe9IVZr6VziJnQOMMNW+HDDdUWcERbDUddT7jkrV+8cdZr\n6RxiJjTOcMPMNEqj2gJYHHENN8u+4RGdQ8yExhklcbwD3FBtAUe01cDilr3Kg87hCqJxhpOaI2jM\nTKMUqi3ghrYajmrk1Ms+4ULncMXQOMPNso+gAXux7P/zB+ZFWw03NTfxWuYJFzqHK4bGGaWxTgpu\nqLaAG2Y7gMXVzKmXecKFzuGKoXFGSayTghuqLeCG2Q64YnB6OdA5XBAHDsMJ66TghmoLuGG2A44Y\nnF4edA4XwIHDcMMoGkpjZBhuiGlgcQxOLw86hwtgRBolsU4KbhgZhhtiGq4Y9OgvOocLIJBRCuuk\n4IiRYbghpuGIQY9+o3O4AAIZpTArDUcMqMENMQ1HDHr0G53DBRHIKIGEA11gQy24IaaBxZGD9Ftk\nZu1nKGYwGOT6+nrtx4CBGuv/atwT/VGzdBkogZiGo1q5ADmIl4gYZuZglmuZOQR2wa60cETpMtwQ\n03DD/gOoYV/tBwCW3bSEA2jTcGNTR44e03Bjs7N7jsuG1kKUDaF1xDSwOPIP1LCnmcOI2CfpaySd\nIWk9M29s9amAHbC9MpzUnJlmrRRKIKbhqEaZJfkHapi7cxgRhyW9QE3HUJK+SdJfRcSbJV2Zmb/c\n4vMBd8H2ynBTsxSOsiGUQEzDDQMe6JO5OocR8cOSfknS5ZKukPTGiW//maQnSKJziGLYXhluGBmG\nG2IabhjwQJ/MO3N4saRXZubPRMTalu/9jaSfauexgOlIOlAaW+HDDTENNywvAcqZ6yiLiLhZ0qMy\n88pR5/BLkgaZ+VcRcb6k/5GZJ5V51PlxlEVZbK8MN2yFDzfENNzUimlyD5RWMsZKHmVxnaRztvne\neZI+OefnYUWxvTIcsRU+3BDTcMPyEjhapoG8eY+y+O+Snh8R5068lhFxpqSfkPTm1p4MS43tldGF\nrrfDZyt8lEZMww0xDSxumfLqeWcO/72kR0j6kKR3S0o1G9A8SNK1kl7Y6tNhaVF/j9LYmRZuiGm4\nIabhqO/HlszVOczM6yNiIOnHJX2npL8bfcZ/kvSqzPyn9h8Rs2DDAbihdAhuiGm4IabhhmNL9nDO\nYWZ+QdKLRl9YAjUDmcYZpSzTKBrQBmIabohpuOHYkj10DrF82HAApdUosVimUTT4IabhhpgGFseA\nxwydw4i4co7Py8z89gWeB3tAIKMkdqaFG2IabohpOGLAo45ZZg73qdl4Zuw8SfeV9HFJn5H0FWqO\nt/i0pKvbfbzVxPo/OGFmGm6IabghpuGGAY96du0cZub5499HxOMk/ZKkh2fmeyZef5ik3x19r9dY\n/wc3zEyjtK4H1IhplEZMA4thwKOeedccvkjS8yY7hpKUme+OiEslvVjSf2vp2VYSwYzSmJmGE7bC\nhxtiGo4Y8OiPeTuHD5T02W2+d62kA4s9zuojmFESM9Nww1b4cENMww0DHv0yb+fwY5J+RNIfTfne\nj6hZh7hUmGWBE2am4YYBNbghpuGGAY9+mbdz+POS3hARH5L0Jt25Ic0TJT1I0pPbfbzF3HTr7cyy\nwApJB0pjQA1uiGm4ocQTJUVm7n7V5BsivkNNJ3Eg6R6SviTpvZJekJl/3PoTLuCs874uT/ju/6Dj\nKa2FdPEF5+nwod5XvqJFNbZZrnFP9EPN3eGAEohpuKkV0+Qeqy0ihpk5mOXaeWcOlZnvlPTOiNgn\n6UxJ12Xm8Xk/pwun3fMEHWekA4Ww/g9uKFuGG2IabijxRGlzdw7HRh3Ca1t8ltadcuKafo3SDhRC\n0oGSaozSUjqEkohpuCGm4WiustKIeP4ul2RmvmixR2rPYDDI9fX12o+BDtQq73zya951RwNNuRLa\nUrMUjtIhlEBMww0xjVVSsqz00h2+N+5lLk3nEP1Qs7yTTQdQQs1ZaUqHUAIxDTfENFzN1TnMzH1b\nX4uI/ZK+S9JPSnpcS88FzIwGGm4oG4IbYhpuiGm42vOaw7HM3JT06xFxhqQjkh618FMBc6CBRmls\nhQ83xDTcENNAO+Y+ymLbD4p4hKS3ZOZprXxgC1hzWAfHO8AJW+HDDTENN8Q0HLWZ2xY9ymIHF0r6\nbIufhxXE8Q5ww660cENMww0xDTc1Bzzm6hxGxOVTXj5R0kMkfZ2kF7TxUFhdNNBwQ9kySuu68oGY\nhhtiGm5q5tPzzhw+QnfuSjp2s6QNSZdJen0bD4XVRQON0lhXAic1RoeJaZRGOw0spmY+3dqaw2XU\n9zWHtdbhsf4PpbCuBG6OHD2mV15xtY6ntBbSxRecp8OHDtR+LGDPaKfhaNX30yi25jAiniLpbZl5\n/ZTvfZmkCzPz1+f5TJRRs3Fm/R9KoWwZbqi2gBvaabjp234adzu3cBevlfRV23zvAaPvYwlMa5yB\nNg03NnXk6DENNzY7u+c4kV4LkUijdTVielwOd/EF5zHDgtbRTgOL61tOPe+aw9jhe6dKum2BZ7HG\nhgNwUnMUjXUlKIFqC7ihnYajGuWdfcupd+0cRsTXS/rGiZe+KyIesuWykyU9SdJHW3w2G2w4ADc1\ny4ZIpFECpXBwQzsNNwx4dGOWmcPH6s4jKlLSz21z3fWSntHGQ7mp1UDTOKOUvo2ioXtUW8ANMQ0s\nhgGPbszSObxM0uvUlJReI+m7Jb1vyzW3SPpMOm99ugAaaJTGtuFwQrUF3BDTcMSAh6ddO4eZ+XlJ\nn5ekiHiApE9n5q2lH6wkEmk46dsuWvBHtQXcENNww4CHr7k2pMnMjVIP0hUSabhhrRTcMDoMN8Q0\n3DDg4WuWDWlul/QvM/M9EXFczbrD7WRmzrsDaqdIpFESu2jBTY2YZnQYJRHTcEPugTbN0pF7oaRP\nTPx+pdcVEswohV204IbjHeCGmIYbcg+0bZY1hz8/8ftLiz5NBwhmlMIuWnBDpQXcENNwQ+6Btu2b\n5+KIuHy0Kc20750dEZe381hlHTx7vw4fOkBAo1XjWem1ELPSKGK4sakjR49puLHZyf2IaZRGTAOL\nIabRtpjn9InRmsOHZ+Z7pnzvoKT3ZOZai8+3kMFgkOvr67UfA5XUqMGvcU/0Q63SIWIapRDTcETu\ngWUUEcPMHMxy7V42j9muN3lfSV/cw+cBrWNXWrhhZzi4IabhhtwDDmbZrfTxkh4/8dLPR8R1Wy47\nWdK3Shq2+GzAnrGuBG7YTAulcaA1sBhyDziYZebwLDUdP6mZNfx6SbdsueYWSX8h6ZL2Hg3YO5IO\nlNZ1Is1mWiiJA63hiAEPYH7zrjn8mKTHZeb7yz1Se1hzuDyowYeTmtvhAyUcOXpMr7ziah1PaS2k\niy84T4cPHaj9WMCesaYVbhaJrWJrDjNz6k6lwE6owYcbSofghhkPuGFNK5x0mUvvZUMaRcR+SQ+U\ndNLW72Xmny76UPBCIo2SaozSkkijpBoxTYknSqKdBhbTZS49V+cwIk6SdLmk75UU21y2NEdZYDnQ\nQKOUmrPSJNIooWbJMjMeKIF2Go6c17POO3P4PEnnS3qqpN+QdFjSzZKeJumfSXpOi8+GQthIAy5q\nzkqTSKMEKi3ghnYabtw38Jq3c/gESS+U9DtqOofvzsy/kvTaiPivkv6NpD9q9xHRJtb/wQmz0ijN\neXQY6AIxDTfu61nn7RyeJemvM/P2iPiSpFMnvne5pNeK2cOlxqg0SusymWZWGiW5jw6jn6geAhbj\nPuAxb+fwekmnjX7/95L+haQ/G/35TEknt/RcvcGoNJzUSqZJNlCC++gw+ofqIThiwKNd83YO3yXp\nG9SUjv6epBdFxOmSbpP0k5L+Z7uP541RabhhZhpOGEyDG9pouGHAo33zdg5frqa0VJJeLOmAmjWI\na5L+UtKPtfdo/hiVhhuSaZTE6DDcUD0ELIYBj/bN1TnMzHVJ66Pff0HSEyLinpLuKembJb1V0kPb\nfkhXNNIoibPS4ITRYbiheghuOM/Sw7wzh3eTmbdIuiUi7i3paxd/pDpIpOGEs9LghtFhuKF6CE44\nz9LHwp1DByTScEMiDTeMDsMNMQ0nnGfpg86hSKThh6QDJVFpATfENNywnhV7RedQBDTKYyMNuKDS\nAm6IabhhPSsWsWvnMCLOnfGz7rvgs1RDQKMkNtKAEyot4IaYhhvWs2IRs8wcHpOUM1wXM163lAho\nlELigZIoHYIbYhpYDDGNRczSOXx68acAOkTiAReUDsENMQ03rGfFqtm1c5iZr+/iQYAukHjACaVD\ncENMwwnrWbGK9tV+AKBL0xKPLhw8e78OHzpAQ41WjWel10LMSsMCMQ0ntXIOYBHsVopeocQTJbEr\nLdwQ03DCshJgd5G5snvI7GowGOT6+nrtx8A2atTh17wvvNUsHwJKIKbhpFY8k3OgpFnjKyKGmTmY\n5TOZOUQV1OHDDbvSwg0xDSesZ4WbUrk0aw5RBXX4KG24sakjR49puLHZyf1YK4WSuo5niZhGWbTR\nwGJK5dLMHEISdfjwwq60cFKr0oKYRim00XDkkkvTOQSNNOxQPgQnNcs7iWmUQBsNN065NJ1D0EjD\nDjPTcEI8ww0xDTdOuTSdQ9BIo6gaO7UxM42SON4BTmijgcU55dIcZbGEajTUbLWMEtgKH26IaTgh\nnuGIo9LujqMsVljNjQeWLZCx+tgKH26IaTghnuGGo9IWx1EWS4YjHlASW4fDDTENJ8QzsBjy6MUx\nc7hknGqWsVycdtICJGIaXohnOHI53qFP6Bzugo0H4MJpJy1AIqbhhXiGGwY8VhOdwx2w/g9OGE2D\nG2IaTohnuGHAYzVV6RxGxJqkdUmfzMwLJ17/FUlPz8zTdnjvWZI+LOnSzPzFks/JQm2UxKw0nLAd\nPtzQRsMNJZ6YRa2Zw+dI+oike41fiIiBpPvM8N5XSfqjQs91FwQ1SmFWGk7YHQ5uaKPhhhJPzKrz\n3Uoj4v6SHi3pNROvrUl6haSf3uW9j5N0jaS/LvmMY+OgvviC8zj7B61iNy04IZ7hhpiGm1oxffDs\n/Tp86AA59AqpMXN4mZpO4OkTrz1b0lsy89MRMfVNEXGqpJ+R9EhJz93uwyPiIkkXSdJZZ5218MMy\niocSmJVGSZQOwQ0xDSyGmMasIjO7u1nEhZIelZnPiojz1XTyLpL0RknnZ+ZtEXHDtDWHEfGLkt6T\nmW+MiEsl3bDbmsPBYJDr6+ut/xzwUmOtVM37wlutcjjiGaUQ03BD3oGuRcQwMwezXNv1zOG3SHpM\nRDxK0klq1hz+taRbJB0bzRqeEhHHMvPAlvc+TNITI+IX1KxNPB4RN2fmf+ru8eGGtVJww+5wcENM\nwwl5B5Zdp2sOM/OSzLx/Zp4j6UmSrszM/Zl538w8Z/T6TVM6hsrMb5245jJJL6VjiEWxrgRuxqVD\nayFKh1DEcGNTR44e03Bjs5P7EdNwQt6BZbfU5xxGxGMkDTLz+bWfBZ6owUdpbIcPJ+x4CDesZwXu\nqtM1h11jzeFq6rqhpgYfpdQsHwJKOHL0mF55xdU6ntJaSBdfcJ4OH7pbsQ+wEljPCkfT4muZ1xwC\nO6o1Kk3jjBJqrZUCSmHWA05Yzwo3beTRdA6xVEimURLlQ3BSY/aBEk+URBsNLKaNPJrOIbZVI/Gg\noUYprJWCE3Y8hBvaaLhZ1TyaziGmqpV40FCjFMqH4IQqC7ihjYaTVc6j6RxiqpqJBw01SmBWGiVR\nDgcshpiGk1XOo+kcYioaaZTEWik4oRwObmij4YYBvNlxlMWKqNFQs9UySuB4B7jheAc4oY2GG44s\n4SgLOzXrlmsHM/ywVgpuVnmEGNiKNhpuWM86HzqHK4CGGk5IpFES5XBwQzkcsBhiej6Ula6A8czh\nOKgp8UCbKFmGC8rh4IZyOLipFVt9j2nKSgvrOsAYlUYplCzDCVUWcEM5HJxwPutqoHM4J5JpOCGZ\nRmldDqZROoTSKPEE9o6cYzXQOZwTgY2SSDzgpOvBNKosUBJHlsANOQemoXM4JwIbpZB4wE2NwTSq\nLFAKJZ5wQs6B7dA5nBOBjVJIPOCGwTQ4IZ7hhJwD26FzuAcENkog8UBJHPEAN2wOB+wdOQe2w1EW\nwDY44gEuOOIBbohpuCHnQEkcZQEsiF1p4YSNtOCGmIYTcg4sk321HwBYRtMSD6Atw41NHTl6TMON\nzU7uNy4fWgtRPoTWdR3PEjENL+QcWCbMHAJTUIuPUtghDk5qzngQ0yihRqklOQeWCZ1DLD020oAT\ndoiDk5rlncQ02sZgBxwNNza1dtoZ9531ejqHWGo1Nx0g8UAJjBDDCfEMJwx2wM04j147df9Xzvoe\nOodYamw6gNLYDh9OiGe46TKmGeyAm3EerZj9PXQOMTcaarhghzg4IZ7hpuuYZrADJdVcz6o5zi6k\nc4i50FDDCTPTcEI8w02NmGawAyXUXs/6sMv+8VOzvoejLDCXGtstHzx7vw4fOkBjjdaxHT6cEM9w\nQ0zDRc3jSg6evV+333D9P8x6PTOHmAtlniiJ9VJwQjzDCTuHA3u3SvlzzFGCunIGg0Gur6/Xfoxi\najTUNe8LbzV3pgXaRjzDCfEMNzVy2Zr5c0QMM3Mwy7XMHK4ojniAG9ZLwQnxDCfEM5ywedfOWHO4\nomrWLsPfcGNTR44e03Bjs7N7srYEpRDPcNN1TBPPcEIOvTNmDlfUKtUuY7XU3lGLkmW0iXiGmxox\nTTyjpK7LLcmhd0bnsCVsPAAXNcuHVqXkAquDeIabWjFNPKMEBjuWD53DFlC7DCeMqKG0LgfTiGe4\nIabhhMGO5UPnsAUs1EZJzErDSdeDacQzSqONhhNKPGHdObzp1ts7uQ+BjVKYlYabGoNpxDNKoY2G\nE0o8IZl3Dq+57gYNNzYJbKwsZqXhhsE0OKGNhhNKPCGZdw4zRWBjpZFIo6QaB/IymIZSasQzbTRK\nosQTNURm1n6GYk663wPzz//yPSQfaE2N5KPGPeGvVjkcUELNeKaNRgm1Ypp49hQRw8wczHKt9czh\nuWeeRmCjNawtgRPK4eCEI0vghhJP1LKv9gOUdMqJa7UfAUamNdTAqhqXD62FKB9C64Ybmzpy9JiG\nG5ud3I94hhtiGrVYzxzCG7X4cMJ2+HDBjodww/ps9AmdQ6wkkg84oWQZTiiHg5Oa61mJadRgXVYK\nX7VKPA+evV+HDx2gsUarKFmGE8rh4IT2GX3DzCFWEiWeKImSZTihZBlOaJ+BsqyPshgMBrm+vl77\nMXqBIx7ggu3D4YQjS+CE9hnYG46yQKdYLwUnrJeCE44sgRPaZ7hZxoEHOodYGMkHnFBChJIoiQP2\njniGk2Wt7KBziIXRWKMk1kvBBbssww3tM5x0Hc/LOrlC59AQjTVcULIMJ5TEwQntM5zUiOdlnVyh\nc2iGxhpOlnVUDauvxjqPZU0EgL2gfYaTGvG8rJMrdA7N0FijFJJpuKg5iLaMiQA8sJ4V2Lta8byM\nkyt0Ds3QWKMEkmk4qTmItoyJAFYf61nhhiVS9dA5LKjGTAvBjRJIpuGEQTS4YT0rnLBEqi46h4XU\n3J6W4EbbSKZREiPEcEIJPrAYlkjVReewEAIbpXWZgJBMoxRGiOGEEny4YbCjf3rTOWShNpzUWl9C\nwoG2MZAGJ5TgwwmDHf3Ui84hC7XhhoQaLhhIgxPiGU4Y7OinXnQOWagNNyQgKIX1f3BCPAN7R67R\nT5GZtZ+hmMFgkOvr63fMHI6Du8vNYeCvRj1+zfvCV82NtIC2Ec9wU+P/++QaHiJimJmDWa7txcwh\nI3kohV1p4YRyZTghnuGEzbvQlX21H6ArB8/er8OHDhDgaNW05ANoy3BjU0eOHtNwY7OT+41LiNZC\nlBChdcQzsHfkG+hKL2YO0R/sSgsXbKQFJ8Qz3JBvwBWdQ9gg+YATNtKCE+IZTsg34IzOIWyQfMAJ\no8RwQjzDCfkGnNE5hA2SD5TElvhwQjzDCSWeQHt6cZQF6mDLZbhgS3w4IZ7hpFY8k2+glBKxxVEW\nqFqatqEAAByUSURBVI4tl+GELfHhhHiGE0o84WQZBu96c5QFusWWyyiJLfHhhHiGi65jWSKe4WUZ\n8mdmDnuCeny4YJc4OCGe4aJmxRDxjBJqlA4vQ/5M57AHSD7ghBIiOCGe4aJmuTLxjLb1ebCDzmEP\nkHzAyTKMqgFtIZ7hgliGkz4PdtA57AEabJTElvhwUaOEiHhGKbTNwN71OXfmKIsKOOIBLpZhVy2g\nDcQynBDPcEPuvBiOslhiHPEAJ2yJDxfEMpwQz3BC7twtjrLo2DJsUQu0hS3EUQrHO8AJ8QzsHblz\nt5g57Fifa5hRFuul4IIdluGEeIYbjkfz1uvOIck0XNRcX9LXsguUww7LcEI8wwmDHf562zkkmYYT\n1pfACaPEcEI8wwmDHf562zkkmUZJlFzACVviwwnxDCfkG2hbb4+yGM8cjoObbZ7Rllqz0k5bLmN5\nsCU+nBDPcEK+gVlxlMUMGMlDKZRcwAlVFnBCPMMJ+QZK6G3nUCK4UQYlFyiJEiI4IZ6BvSOeUUJv\ny0rRHzXKHyi5QAmUEMEJ8Qw35BtYVpSVAiO1kg9mpVECJURwQjzDCfkGXOyr/QBASdOSD2BVjUuI\n1kKUEKF1w41NHTl6TMONzU7uRzzDCfkGXDBzCGvU46MktsSHCw62hpMapZbkG3BB5xCdIpmGC0qI\n4IQST7io2TaTb6CErnNnOofoDMk0nLAlPpww6wEXNdtm8g20rUbuTOcQnSGZhhOSaZRSoySOWQ+U\nwnElwN7VyJ3pHKIzNNgohWQaLmpVWEjMeqB9rGWFmz4MdtA57DHW/8EByTScUGEBJ6xlhZO+DHbQ\nOewp1v/BBck0SurDKDFQCvEMJ30Z7KBz2FMk1Cipy4Sa5AOl9GWUGP1BxRCcMHhXBp3DJUGAw0XX\nCTXJB0rpyygx+oGKIThh8K4cOodLgACHkxoJNckHSmAQDU6oGIITBu/KoXO4BAhwOCGhRimUxMEJ\nFUPA3hHP5URm1n6GYgaDQa6vr9d+jF2NZw7HAd7lbovwVuOIh5r3ha+au9ICbasVz7TNKKVGbBHP\ns4uIYWYOZrmWmcMpGJ2GA454gBNK4uCEiiE4YT2rFzqHWxDgcEEyjZIoiYMT4hnYO/INL3QOtyDA\nURJHPMABm2jBCfEMJzVKLck3vFTpHEbEmqR1SZ/MzAsnXv8VSU/PzNOmvOeRkv6DpBMl3SrppzLz\nyrafjQBHKRzxABeUxMEJ8QwXNavfyDd81Jo5fI6kj0i61/iFiBhIus8O77lO0ndl5qci4iGS3i7p\nK9t+MAIcpXDEA1wwiAYnxDNc1Kx+I9/w0XnnMCLuL+nRkl4i6eLRa2uSXiHp+yU9ftr7MvN9E3/8\na0knRcQ9M/OWtp+RAEcJJCAohU204KJGSRzxDBfkGWhD50dZRMSbJL1M0umSnpuZF0bEcyTty8xX\nRcQN08pKt3zGEyX9aGZ+x5TvXSTpIkk666yzDm5sbLT/Q2DlccQDXHDEA1wQy3DD8Q5YFkt7lEVE\nXCjp2swcRsT5o9fuJ+l7JJ0/42d8raSXS7pg2vcz89WSXi015xwu/tRwwxEPcMImWnBBLMMJu99j\nVe3r+H7fIukxEfFxSb8j6RFqSkQPSDo2ev2UiDg27c2jktQ/kPSUzPy7Tp4YdqYlIEBbhhubOnL0\nmIYbm53cb1xGtBaijAit6zKeiWU4IdfAqup05jAzL5F0iSSNZg6fO7lb6ej1GzLzwNb3RsR9JL1N\n0iWZ+ecdPC5MUZOPUtgSH07YYRlOOMsSmM1Sn3MYEY+RNMjM50t6tpoZxudFxPNGl1yQmddWe0C0\ngs004IIt8eGEHZbhgoE7OCmdN1frHGbmVZKumvL6aRO/f4ukt4x+/2JJL+7o8dARavLhhJFiOCGe\n4YKBO7joIm9e6plD+GMDApTErDRccMQDnFDiCexNF3kznUPcBQ02XDArDRfssAwnlHjCTZe5cxd5\nM51D3IEGG06YlYYLYhlOKPGEE8eNu+gc4g402HDCrDRKocICTohnYO8cN+6ic4g70GCjJNb/wQEV\nFnBCPAOLccyd6RwuMZJpuGD9H1xQYQEnxDOcsHFXO+gcLimSaThhzRRcOI4So7+IZ7hg46720Dlc\nUiTTKIk1JnBBhQWcEM/A3pA3t4fO4ZIimUYprDGBCyos4IR4hhMGoVcXncMZMZoHF6wxgQtGiuGE\neIYLBqFXG53DGTCaByeMrsEFsQwnxDNcMAi92ugczoDRPJTCzlpwQoUFnBDPcEGJJ+ZB53AGBDlK\nYGctOKHCAk6IZ7igxBPzonM4A4IcJTAjDSfEM5wQz3BBiSfmRedwRgQ52saMNEqijAguapTfE88o\nhbYZyy4ys/YzFDMYDHJ9fb32Y2CFdN1o10h64K9WSRzxjLbVLL8nntE22mbUEhHDzBzMci0zh8BI\nrbp8Gmq0jTIiuKhZ3kk8o220zVgF+2o/ALAspjXawCoalxGthSgjQquGG5s6cvSYhhubndyPWIYT\n4hmrgJlDYIS6fJTClvhwwK6HcEPbDNwdnUMsJc7/gwu2xIcLSuLghLYZTtrMm+kcYulw/h+csCU+\nXFBdASe0zXDRdt5M5xBLhwYbJbGNOFxQEgcntM3A3rSdN9M5xK5osOGCNVNwQUkcnNA2w8mq5810\nDrEjGmw4Yc0UXFBhASe0zXDhkDfTOcSOaLDhhFlplLLqI8VATcQzXDjkzXQOsSMabJTCjrRw4TBS\nDIzRNsMJA3fzi8ys/QzFDAaDXF9fr/0YrarRaNe4J7zV3JEWaNuRo8f0yiuu1vGU1kK6+ILzdPjQ\ngdqPBcyNthlOasXzMubNETHMzMEs1zJzuELYgAAuWC8FJw4jxYBE2wwvDiWeNdA5XCE02nBBMo2S\nOOIBLiiJA/aOeN4bykpXyHjmcBzklHugTV0nIctYdoHVR1kcXFASBzcsjaqHstKOMDoNF7U21CCG\n0TYqLOCCkjg4YWnU6qBzuEcEOZyQUKMUyuLgglgG9o48Y3XQOdwjghwlkYTAAUc8wAWxDCc1Si3J\nM1YHncM9IshRCkkIXFAWBxfEMlzUrHwjz1gNdA73iCBHKSQhcMEgGkqhugLYm5qVb+QZq4HO4QII\ncpRAEoJS2EQLDqiuAPaOHAO7oXMI7KBGXT5JCEpgEy24oLoCbrrMNcgxsBs6h8A2ap6XRhKCtrGJ\nFlww8wEnHCWFZUPnENgGyTRKYs0UXFCuDOwduQaWDZ1DrJQukxCSaZTCmim4oFwZTjjiAaBziBXS\ndRJCMo1SWDMFF8x6wAVHPMDN5GDHPOgcYmXUSEJIplECI8VwQSzDBUc8wMnWwY448eRTZ30vnUOs\nDJIQlMKaKbggluGCddnA3m0d7Nh34imnz/reyMySz1bVYDDI9fX12o9hqUZdfs37wlfNXWmBNhHL\ncFErlskxUEKNuBr/NzQe7PjbX3jC3xy/5aYHz/JeZg4xN454gBPWTMEFsQwXrMuGi2VZyzp48Rdv\nnPW9+0o+GDxNa7SBtgw3NnXk6DENNzY7ud+4lGgtRCkRWkUswwWxDOxNzZz54Nn7dfjQgbk7o8wc\nYm7U5aMUjniAC2IZLohluOFYtJ3ROTTRZaDTaKMUSonggliGC2IZTjgWbXd0Dg3UGtVbhQDHalnF\nETZgGmIZLohlOOFYtN3ROTTAJgQohW3x4YJYhgtiGdg7Bjt2R+ewAM7mgYOaO2yRfKBNxDJcEMtw\nUuOIBwY7dkfnsGUs3IYLZqThgliGC2IZLjgWbXlxlEXLam1Zu9ftaoHtsJU4SmFbfLggloG94Vi0\n5cXMYcso8UQprDOBA6or4IJYhhOWRGHMunN4062368jRY9QyY+WxzgQu2BYfLohluGCgA5OsO4fX\nXHeDXnnF1dQyY+WxzgQuGC2GC2IZLhjowCTrzmGmSKZRBOUXcEG5MlwQy3BBjoGaIjNrP0MxJ93v\ngfmVT71M9+h45hDeapV41tjyGd5q7hYHtIlYhgtyDJQQEcPMHMxyrfXM4blnnqbDF5xHoKNVlF/A\nBeXKcEEswwU5Bmqz7hyecuKaDh86UPsxYIbyC5RCKRFcEMtwQSyjb6zLSgeDQa6vr9d+DBRWoxSC\n8gu0jVIiuCCW4YJYhgvKStEbHPEAF5QSwQWxDBfEMvpoX+0HABYxreEGVtG4lGgtRCkRWjXc2NSR\no8c03Njs5H7EMlwQy3ARJ5586qzXMnOIlUZtPkphW3w44HBrOKFdBvZm34mnnD7rtXQO0Soabjig\nXBkuKIuDC9plOOk6Xz5+601fmPVaOodoDQ03XLAtPlxQXQEXtMtwUSNfzlu/eOOs19I5RGtouOGC\nhBqlUF0BFxzxAOzNsufLdA7RGhpulEJCDQdUV8AFa1nhhIGOu6JzaIyEGg5IqOFi2UeLgVmxlhUu\nGOi4OzqHpkio4YKEGqUwWgzsDbEMFwx03B2dQ1Mk1CiFhBoOGC2GEyqF4IIcoz46hx0h2OGAhBou\nGC2GCyqF4IIcYznQOewAwQ4XJNRwwQAaXFApBBfkGMuBzmEHCHa4IKFGKZTFwQWVQsDeEMvLITKz\n9jMUMxgMcn19vfZj3DFzOA72rko+4K/rJKTWPeGtVlkc0LZasUy7jBLIMXxExDAzB7Nc28uZQ0ao\n4YB1JnBBWRxcUCkEF+QY/dW7ziHBDhck1CiFsjg4qDEDQSzDBTlGf/Wuc0iwoxQSajhgAy04qDkQ\nTCyjBHIMdKV3nUOCHSWQUMMFZXFwUHMgmFhG28gx0KXedQ4JdpRAQg0XDKDBAXEMJ+QY6FLvOocS\nwY72kYigFDbQgosuY5k4hhNyDHSJoyxgie2X4YAjHuCCWIYTcgysGo6yQK+xIy1csIEWXBDLcEGO\nAXf7aj8A0LZpSQiwqOHGpo4cPabhxmZn9xyXEq2FKCVCq7qOZ2IZLsgx4I6ZQxTH9stYdWyLDyfs\nfAgn5BhAu+gcoiiSEDhgW3w4YedDuCDHgJNlWVdK5xBFkYTAASPFcEI8wwU5Blws06ZddA5RFEkI\nSmFbfLjguBI4qDHrQY4BF8u0aRdHWfQM2y/DwTKNsAGLIJbhoGYck2OgbbVy5Se/5l13DHS0/d8Q\nR1lgKrZfhotlGmEDFkEswwHrsuGCDeg4yqJX2H4ZLtgWH6X8n/buPtiuqj7j+Pcx4T3YRGgISZCk\nxFAgUx25FlEsAVsCKeIrTiqYKa2TMo0dKq3WFJFQX1rGcYaKb1haqQKD7yWIEkxrfCMZvClCEmIw\nNIghvBiMwg0mQPLrH3td3Zyc+3Lu3fvsc/Z5PjN7zt3rrLP32vNb59y19tp7bT/iweqinXXZ9djq\nosq28snHTmHpGXMqP9nhkcMe4mvzrSy+Z8rqwDMfWl20uy67HltduK3szmGl3KC2OvDlylYXnvnQ\n6qKKuux6bGVwW7n93DmsiBvUVhe+Z8rqwmeMrS5cl60O3FauhjuHFXGD2srS7rNsboRYWXzG2OrC\nddmsdW4rV8Odw4q4QW1l8D1TVhc+Y2x14bpsddLOEx1uK1fDncPEZ/WsDnzPlNWFzxhbXbguW114\noqPe4M4hPqtn9eGzbFYXrstWF67LVhee6Kg3uHOIz+pZeTwibXXhumx10O56DK7LVo4q6rJPdPQG\ndw5xZbdyeETa6sJ12eqgqnoMrstWrCp/k32io/5eUMVOJU2QdLekrzekXyNpYJjPLZO0RdJmSQuK\nKs9gZb/0rOPb+s/C6q3ZiLRZN3JdtjpwPba6qLIun3zsFJaeMcdt5RqrauTwEmAT8MLBBEl9wOSh\nPiDpRGARcBIwHVglaW5E7C2iQD6rZ0XziLSVxY8rsbrwzIdWB/5NtjpRRLR3h9JM4D+BDwGXRsS5\nkiYAq4C3AT+JiElNPrcMICL+Oa2vBJZHxJqh9tXX1xf9/f0lHIV1oyquz69in1ZvVV1O5LpsRaui\nLrseW9H8m2zdQNK6iOgbTd4qRg6vBt4DHJ5LeyewIiIekTTU52YAa3Pr21La80haAixJqwOSNo+7\nxMU7EthRdSF6iQ485LADpkyfC3oBxL5nd26/P5759a4WNuGYdZ9axmzCpCOmTThsygwERMQpV/9y\n+96BJx6tulwFqGW8am5cMatxXe5k/p4VrA312DHrLp0ar2NHm7GtnUNJ5wKPR8Q6SfNT2nTgfGD+\nSB9vkrbfsGdEfAb4zPhKWi5J/aPtvVtncMy6j2PWXRyv7uOYdR/HrPs4Zt2lDvFq98jhq4HzJC0E\nDia753AjsAfYkkYND5W0JSLmNHx2G3BMbn0msL38IpuZmZmZmdVfW2crjYhlETEzImaRTS7zPxEx\nJSKmRcSslP50k44hwApgkaSDJM0GXgLc1bbCm5mZmZmZ1VhHP+dQ0nlAX0S8PyI2SvoicB/wHLC0\nqJlKK9DRl71aU45Z93HMuovj1X0cs+7jmHUfx6y7dH282j5bqZmZmZmZmXWetl5WamZmZmZmZp3J\nnUMzMzMzMzNz57BIks6XtFHSPkl9ufQjJH1b0oCkjzd8ZrWkzZJ+lJapTbb7J5LWSVqfXs9sx/H0\ngrJilvItk7Ql5V1Q9rH0grHEK5dnhaQNQ7z3O5JulXRP2v5FZR1DrykrZun9+ek7uFHSd8oofy8q\nM2Ypzysk7ZX0lqLL3qtK/G28QNK9ablT0kvLOoZeUmK8JOljqe1xr6SXl3UMvWaM7cXbc+2KT0ua\n0GS7Hdf+6OgJabrQBuBNwLUN6buBy4F5aWl0QUT0D7PdHcDrImK7pHnASmBGAeW1kmIm6USyGXlP\nAqYDqyTN7eJJlDrFmOIl6U3AwDDbXQrcFxGvk/S7wGZJN0bEM8UUu6eVEjNJk4FPAmdHxENDnaSx\nMSnre0ZqHF1F9n/MilNWzLYCp0fETknnkE22cUohJe5tZcXrHLLZ/F9CFqdP4XgVZSwxe2tEPClJ\nwJfJnut+c0Oejmt/eOSwQBGxKSI2N0nfFRHfJ6tAY9nu3REx+EzHjcDBkg4aR1EtKStmwOuBmyNi\nT0RsBbYAfziOohpji5ekScClwAeH2zRwePoBnwT8gmxWZBunEmP2NuCrEfFQ2t7jBRW555UYM4C/\nAb4COF4FKitmEXFnROxMq2vJnjFt41Tid+z1wOcisxaYLOnoosrdy8YSs4h4Mv05ETiQrK2xXzY6\nrP3hzmFn+Gy6NOryVDmG82bg7ojY046C2ZBGitkM4Ge59W14tLcqHwA+Cjw9TJ6PAycA24H1wCUR\nsa8NZbPmRhOzucAUZZd5r5O0uD1FsyGMGDNJM4A3Ap9uV6FsWKP5nuX9JfDN8opjIxhNvNz26DCS\nVpKdDHuKbPSwUce1P3xZaYskrQKmNXnrsoi4ZQybvCAiHpZ0ONnZ1LcDnxti3yeRXY5z1hj207Mq\nilmzDqOfGzMKRcZL0suAORHxLkmzhsm6APgRcCZwHPAtSd/LnfWzYVQUs4nAycBrgUOANZLWRsT9\nreyvV1UUs6uBf4iIvSOfB7VGFcVsMP8ZZJ3D01rZTy+rKF5ue4xDCe1FImKBpIOBG8naGN9qyNJx\n7Q93DlsUEX9c8PYeTq9PSbqJ7NLD/TqHkmYCXwMWR8QDRZah7iqK2TbgmNz6TLKzQjaCguN1KnCy\npAfJfu+mSlodEfMb8l0E/EtkD37dImkr8PvAXQWWpbYqitk2YEdE7AJ2Sfou8FLAncNRqChmfcDN\nqWN4JLBQ0nMR8V8FlqW2KooZkv4AuA44JyKeKLAMtVbh76LbHmNUdHsxt93dklaQXfbb2DnsuPaH\nLyutkKSJko5Mfx8AnEt2w2tjvsnAbcCyiPhBe0tpeaONGbACWCTpIEmzyW4Od0ejzSLiUxExPSJm\nkZ3xvr9Z4wd4iGwECklHAccD/9euctpvtRCzW4DXpO/koWSTLmxqX0lt0GhjFhGzI2JWyvdl4K/d\nMazGaGMm6cXAV4G3e1S+Oi38Lq4AFqdZS18J/CoiHmljUS2RNGnwfk9JE4GFwI+bZO289kdEeClo\nIbuXYhuwB3gMWJl770Gym0wHUp4TgcOAdcC9ZBPN/CswIeU/D/in9Pf7gF1kw86Dy9Sqj7cOS1kx\nS+uXAQ8Am8nOuFZ+vN2+tBqvhs/OAjbk1i8GLk5/TwfuILvefwNwYdXHWpelrJil9XcD96WY/W3V\nx1qXpcyY5dKvB95S9bHWZSnxt/E6YGeu7dFf9bHWYSkxXgI+kdoe64G+qo+1LkurMQOOAn7Ib9uL\n1wATm8Ss49ofSgUzMzMzMzOzHubLSs3MzMzMzMydQzMzMzMzM3Pn0MzMzMzMzHDn0MzMzMzMzHDn\n0MzMzMzMzHDn0MzMakDSn0sKSXOqLkterlyzcmnLJZ1Zwr6uTw/JNjMzGxN3Ds3MzMpzG3AqkH8Q\n9RVA4Z1DMzOz8ZpYdQHMzMzqKiJ+Dvy86nKYmZmNhkcOzcysJ0i6UNI9knZL2iHp85KObsjzoKQb\nJC2StEnSLkn9kk5rsr1LUv7dku6S9Kq0fn0uz/MuK5UU6a3LUnpIWp7eWy1pdZP9PG+bKe21kv43\n7fsBSX81xDEfKukqSVslPZNeL5Pk//9mZrYfjxyamVntSVoCXAt8AVgGTAc+DJwi6eURMZDL/hrg\neOByYDfwAeDrkmZFxC/T9t4BXA38O/Al4DjgJmDyCEU5FVgDXJ/KA7CtxWM5AfgG0A8sAg4ClgOT\ngL25fBOBlcCJ6RjWA69Mx/Ui4O9a2a+ZmdWfO4dmZlZrkiaQdY5WR8SiXPqPge8BfwF8LPeRFwIv\ni4idKd+jwA+BhcBNadTtCuCbEfGO3PYeBb4yXFkiYq0kgIcjYu0YD+l9wFPAWRGxK+37TuABYHsu\n358BpwGnR8R3U9p/p/1fIemqiHh8jGUwM7Ma8mUlZmZWd8cDU4Eb84kR8X3gp8DpDfnXDHYMk/Xp\n9cXpdWZavtTwuVuA54oo8AhOBb4x2DEEiIifAT9oyHc22fHdKWni4ALcARxANopoZmb2Gx45NDOz\nuntRen2kyXuP5t4f9Iv8SkTsSaNtB6ekwfsUH2/It1fSjvEVdVSOBh5rkv4YMDu3PhU4Fnh2iO0c\nUXC5zMysy7lzaGZmdTfY2ZvW5L1pZPfutWKwkzk1n5guXz2yxW3l7Sa7pLVRY+f1EeCoJvka054A\ntgJvHWJ/D7ZSODMzqz9fVmpmZnW3mWxUbVE+UdKryEbWvtPi9ral5fyG9DcwupOuzwCHNEn/KTBX\n0oG5Mv4RcHhDvjXAQkmH5fIdA7y6Id/twDHAQET0N1naMcppZmZdxCOHZmZWJ2eniWHyfgW8H7hW\n0g3ADcAM4EPAT4DPtrKDiNgn6Urg3yRdR3bv4e8B70372jfCJu4D/lTS7cBOYHtEbAduBpYA/5Ee\nXTEbuDRtM++DZB3TOyR9BDgQuJL9LzW9EbiIbBKajwL3pLzHAecBb4iIp1s5djMzqzd3Ds3MrE6u\naZK2MSLmSXoaeDfZxDEDZI+DeE/DYyxGJSKukzQJeBdwIbABuAC4lf07c43eSTY76q1kj6G4Elge\nEd+WdDHw98CbgbvTtp83A2pEbJK0EPgI2aM5HgauIpuoZn4u37OSFpB1WpeQdTZ3kc1qehvZCKaZ\nmdlvKCJGzmVmZmbDkvQK4C5gcUR8vurymJmZtcqdQzMzsxZJmg0sJXtO4pPACcA/ko3GzfPlmmZm\n1o18WamZmVnrfg3MAxYDU8juHVwFvNcdQzMz61YeOTQzMzMzMzM/ysLMzMzMzMzcOTQzMzMzMzPc\nOTQzMzMzMzPcOTQzMzMzMzPcOTQzMzMzMzPg/wGu3rDD0PwIhgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fc203bc1278>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(15,10))\n",
    "plt.plot(lons,lats,'.',label='HRAP')\n",
    "plt.axis('equal')\n",
    "plt.xlim(-115,-114)\n",
    "plt.ylim(44,45)\n",
    "plt.xlabel('Longitude', fontsize=16)\n",
    "plt.ylabel('Latitude', fontsize=16)\n",
    "plt.title('HRAP Grid in native projection', fontsize = 18)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Reproject into WGS84"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "p1 = Proj('+proj=longlat +a=6371200 +b=6371200 +no_defs',preserve_units=True) # geographic projection for HRAP grid\n",
    "p2 = Proj('+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs',preserve_units=True) # WGS84 geographic projection"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "lonsWGS84,latsWGS84 = transform(p1,p2,lons,lats) # reproject the HRAP grid to WGS84"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7fc20028d630>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4cAAAJlCAYAAABpFoerAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XuYrFdV7/vfWB0WIRdkCR4jYIKwMHi/rFI4W7euhRr2\ngYggonhFNhiV5TloBD05W7mJoEK2cevaHhCDqGwV8YaAEiEreNsIXbBRLgaaQB8RMBIbIQlJSNY4\nf1RVUulUd1d1vfOdb/3e7+d58iSprup6EwYzY8w55pyRmQIAAAAA9NuB2g8AAAAAAKiP4hAAAAAA\nQHEIAAAAAKA4BAAAAACI4hAAAAAAIIpDAAAAAIAoDgEA2FVEfDAirprzvd8fERkRR+d4729EBPdJ\nAQA6g+IQAHokIo6Oi5en7/KejIjXbHvtqvHrkz8+HREfjojfi4gv3uM7D0XETePPfc8u7/vgtu+4\nZfzaSyPicxf4Zzw9Ip4aEVdGxL+On/XjEfHWiPj5iHjIvL9r1UTED43/3T15xs++c/yzf9rhs1eP\n/z0d2Pb6gyLisoh4Z0R8Yvy/y4cj4nUR8cMRcea2969FxPdGxF9HxEfH/9t/KCJORsRzI+Lu295/\nn4j4hYj4x4i4cfyZKyPiW+b4533l+J/pnfP9GwIA7Oa02g8AAFgZN0t6yviv7yHpiKQnSXpkRAwy\n8+odPvfdkg5K+oCkJ0v67V2+40OSLhn/9dmSjkr6z+Pv+NLM/NhuDxgRD5T0GklfIOlNkn5R0kck\nnSXpy8e/6+kRcW5m/vNuv2vK+ZJKrPD9gKQfavh3nhz/+ZikX9/2s6OSbpV0/4g4nJkbkx9ExDmS\nPl/Sn2bmqanXnyjpxePPvVLS/yvpRknnSPo6Sb8i6TGSHjH1Pf9D0rdL+htJl0raknSupK+W9P9I\n+m8axZIi4gxJfyvpcyX9mqS/l/SZkr5f0h9HxFMz81dn/YNGxIWSHifpU/P8iwEA7I3iEAAwr1sz\nc7qw+7WIeLekX5L0I5L+zx0+92SNipY/kXRZRDwoM9+/w3v/fdt3/GpEXDv+/U+S9MKdHi4i7iHp\ntZIeJOlbM/OPZrzndEk/pj2KvYi4m6S1zLwpM2/e7b37lZmflvTphn/n1RHxYY0Kwe2OSvodSd85\n/uuNqZ8dG//5qskLEfENki6X9E5Jj5xRTD9/XIx/x9RnjmhUGP5RZn7r9geIiM+W9O9TL32LpAdL\n+tHM/KWp971Eo4mCH5R0l+IwIs6S9N8lnZD06Bn/rACAfaCtFACwjDeO//zgWT+MiK/UaMXu5ZJe\noVEx9KQFv+P14z8f3uN9T5H0EEkvnFUYStK42HtBZn546hmfPW5N/KKI+K8R8SFJN0l62PjnM/cc\nRsRTxq2QN0fERkQ8TVLM+w81a8/h5LWI+IyI+NWIuHbclvk3EfHQOX/1VZLuFxG3/28SEZ+j0crg\nayUNdUcxOHF0/OeTU6/9wvjP377TKmtmXpOZL5h6afKdV+7w/n8ZF8UT9xz/+cPb3vrvkm4Y/zHL\nz2o0wf1TO/wcALAPrBwCQD+dERH3aeD3PGj853/b4edP1ijB/4PMvCEiXivpiRHxzOn2xT1MCo5d\nW0olfdv4zy+d8/du9wqNWhQv1Whl8SM7vTEiflSjltV3aNQqeYakZ0i6dp/fvd3rJf2rpOdKurek\niyW9LiIekJmf3OOzJyV9l0YF3/vGrx0d//lNkr5S0va9n0c1av98hyRFxAPG7/vLXdqFZ5msCD8+\nIl6RmVt7vP9KjVpWXxARN2jUVnpIo9Xde2lUBN5JRHy1RivJ35mZn4iYux4HAOyB4hAA+uk54z8W\nMlVQTvYcXjb++7vsIxy3cH6npFdl5mQF6OWSHqvRHrU/m/EVa1PfMdlz+CyNCojf2ePxvljSJzLz\nA9ueY02jgmPaDZm5fa/axyV9Y2beutuXRMSkaHmPpP+QmTeOX3+ZpH/c4xnn9bbMfOrUd75boz1/\n36XRHsDdTO87/LXxXx+VdHVmfnS8CvoTEfH5mfneqVXFP5kq2L9k/Of/tf2Xj/cJnrHt5ety5K0R\n8aeSvlnShyLibyX93fiPN07+XU1k5vsi4js0ak1+7dSP/kXSwzPzb7Z992njf6YrMvOVe/x7AAAs\niOIQAPrpJZJ+f4ef/cUOr5+p0WrWtI9IemJmvm7G+79Vo6Ls5VOvvVaj1bX/rNnF4UNmfMf7JX1X\nZu51IuU9JX10xutfIOkftr32DEkv2vbaZXsVhmMXaFQcnZgudjLzQxHxCkkXzfE79vKL2/5+0qY5\ns313Wma+f3wi6dGpl4/qjv2Efy3ptvFr79XsltJJu+cnZnzFcyX9+LbXPkt3rOw+TqO9gt87/t3f\nOH79kxHxnMy8dNtnP67RiuFLNSpG7zf+/X8SEd+Qme+Yeu8zNPp38NgZzwUAWBLFIQD00/sy8w2z\nfrBLm95NGq0ISaMTJb9P0jdp5/3rT9ao0PtQREzvF/wLjdoO7zPj9NEPanSKpyTdIunD06dq7uET\nuqOomfaB8XNK0pfprkXhxHvn/J4Hjv88a5Xw3XP+jr1cM/03mXnd+H+Xe8/5+ZOSvi8iztfo38vn\na7xSnJmfjIi3a7Sy+BLdURxeNfX5SVE469/niyX9+fivn6FRsTz9rJ/W6BTTXxkfEnRE0iM1OrDo\nRRHx4cz8HUmKiEdIep2kR2Xm5HcqIv5Qo3+/JyR97fi1w5KeKel5mXmnfz8AgGZQHAIA5nXbdEEZ\nEa/S6NqIl0TE2zLz76d+9nkaFR+hnYuu79EdbakTN+xUtM7hnZK+LiI+b7q1dNzS+obxc+22Mnjj\nLj+bNqmeZ5142sgGuMy8bY/v3stJjYr3o7qj0HvT1M/fpNEVIxq/5zqNVu8mJqu0Xz7j2d6n8V7G\n2OXeyvF7P6XRSuVfR8RJSVdoNGkwaRH+SY3+N//zbZ/7aET8lUZXmBzMzFs02gv6b5L+aNtkw2mS\nDo5fuyEzd9wrCgDYHaeVAgD2Zbw/bXJC5/bVuCeNX/8BSY+f8cfVGhUJTXrV+M9P2fVdy5scuvIF\nM34267UapvcdHpX0/m0njr5J0jkRcUyjVcU3Zebtxe64uH6bpK8drz424c3jP99v6rX7SToQs5er\nT9MoT5nkKudJuq+kd2lUnE7+uJ9Grabv0x17LAEA+8DKIQBg38YHivwPjVoYvzYz/zoiDmh0ifk/\nZObMk0Mj4oskPTsivioz39rQ47xU0lMlPSMi1ne4zqKJlb2/0OhU0+MR8bKpA2nur9GBMdVl5mZE\nfEDS12u0cvimbW/5K0mnJD17/PdXzfg1P6nRqamvjIj/Y/r6jyl3+vc5vj4jd2gFfsz4z9Ott+/W\nqDh9vEYH7kx+z+dJ+jqNYuim8ctP1+gE0+3+u0YtzxdrlxNmAQB7ozgEACzr+Rq1iD5H0jdotAft\ncyX9+i6f+QONCpMnS2qkOMzMT0XEozRqdf3D8amcV2h0SM09NTrs5js0Oozln5b4nq2I+GmNVkv/\nNiJ+U6MDan5Io9Wrr1jmn6NBJzU6+OccbbsSIjM/HhF/r1EBNnmvtr3nDRHxZI32GL43In5fozsS\nb5T02ePPXqBRQTYp4L5M0u9FxJs0Kjg/pNFBRg+V9O2SPqnRgTYTz5f0nyT9dkQc1ehAmvtL+mFJ\np2t0TcjtzzPrHzIiXiTp+sx81ayfAwDmR3EIAFhKZl4dEa+U9ISI+Hrd0S76h7t85p0R8d7xZ35s\nxrUS+32WayLiiEZF0bdpdOrlZ2h01+KGRquLv77g3X2zvufSiLheo9WqF2hUbL5Io8vbL1/mdzdo\nUhxKd105nLz25RqdMvquWb8gM39jvPfv/9Lo1NHHSzo4/sw7JB2X9FtTV5X8pUaH1HzT+Ls/W6PV\nxX+S9DJJL5xeVRxfffEfJP0XjU45vUijAvLvJP1cZl61n39wAMD+xNQWAwAAAABAT3EgDQAAAACg\nTnEYEWsR8faIeM34738jIj4QEf9r/Mddjs4ev++JEfG+8R9PbPepAQAAAMBXrT2HT5P0Ht35ct1n\n7LaZPCI+U9KzJA00ultqGBGvzsytok8KAAAAAD3Q+srh+KjvR2l0KMAiHiHpLzLz38YF4V9odMIZ\nAAAAAGBJNVYOL5P0E5LO3vb6z0bEMyW9UdL/nZk3b/v5/XTno8c/pDtfpCtJioiLNDrtTGeeeeaR\nhzzkIU0990w33nKbrr/5Vp1199N0xsG1ot+Ffrr2kzfrXz5x0+1//9n3PF3/29l3r/hEQDOIbbgi\ntuGK2F5Nw+HwY5n5WfO8t9XiMCIulHRtZg7H9xlNXKLRPVQHJb1Eo4t3n7v94zN+5V2OWs3Ml4x/\nhwaDQa6vrzfw5EA9w80tffdL36xP33pKdzvtgF7xlIfpyHmHaj8WzAw3t/Tma67Twx5479bii9hG\nW9qOb2Ibrojt1RQRm3O/t82rLCLiBZK+V9KtGl1ue09Jf5iZ3zP1nqOSnp6ZF2777HdKOpqZPzj+\n+xdLuiozf2en76M4RNNqJNA1vxf9MPmP/S23ntLBlv9jT2yjtFrxTWyjDbUm9ojt1RIRw8wczPPe\nVlcOM/MSjVYJp4vA74mIz8nMj0RESHqMpHfO+PjrJT0/IiZReMHkdwFtqJlAHznvEAMwinnzNdfp\nlltP6VRKn771lN58zXXENmzUim9iG6XVykuIbW9duefwFRHxD5L+QdJ9JD1PkiJiEBEvlaTM/DdJ\nPyPpreM/njt+DWjFrAQDcPCwB95bB087oLWQ7nbaAT3sgfeu/UhAY4hvuCIvQQm1rrJQZl4l6arx\nXz98h/esS3rK1N9fLunyFh4PuItJgjHpsyfBQCltt+wcOe+QXvGUh9EmhOJqtKMR32hDjdgmL0EJ\nre45bBt7Dn2x9w+uarYvAyUR23DFvm10XWf3HAJNYO8fnNXc/weURGzDFfu24aQrew6BudFjD2fs\nj0IbhptbOnFyQ8PNrda+k9iGK2IbTlg5xMqhxx5tYX8UHNU84ZDYRmmM28ByKA6xNA7PgCPal+GK\nFji4YtyGs7bybYpDLIU7duCK/VFoS9sTbHRfwBXjNly1mW9THGIpDMRoA0eEw1WNCTa6L9AWJj6A\nZrSZb1McGiGBhiP2R8FZrQk2ui9QGhMfcOWeb1McmiCBhiv2R8EZE2xwxcQHHPUh36Y4NEECDVck\nz2gLpxzClftKB9CWPuTbFIcmGITRBpJnuOKUQ7jqw0oH+ot9rc2jOCyE6x3ghuQZzjhcC676sNKB\nfmJfaxkUhwVwvQMckTyjLbTAAc0htuGKfa1lUBwWQBKNNtBKAUe0wMEZrflwxaSeD/vikGCFI1op\n4IoWOLiiNR+umNTzYl0c3njLbQQrLNFKAVdMrsEVXUVwxaSeF+vi8Pqbb9VpBCsMkUCjDbTAwRmt\n+UAziG0vkZm1n6GYL/zSr8hT3/L824O1zRYO9EeNBLrm96IfarbAAaXVim/GbZRGToJZImKYmYN5\n3mu9cnjGwTW9mBloFMQeEriiBQ7OaM2HI3ISNOFA7Qco7ch5h3T82GECFkXMSjAAB5M2obUQbUIo\nari5pRMnNzTc3GrtO4lvOCInQROsVw7RP+whgSP2/sEVpxzCGTkJVhHFIWxwvQMc0SYEZ5xyCFfk\nJFhVFIewwR4SOGLvH5yx0gFX5CRYVRSHsEGSgTbQJgRXtC/DVY3YZuzGqrK+ymIwGOT6+nrtx+gl\njlKGI46/hyuuLoGrmrHN2I2u4CoLVMUeKbiiTQiuaF+GK/a1wlmJCQiKQzSOJAOuaBNCG2iBA5pD\nbMNVqcUYikM0joEYbWB/FBxxtQOcMW4DzSm1GENx2ANtD8YMxCiN1mW4ogUOrhi34cyp44Pi0FzN\nWWgGYpRC6zJc0XkBV4zbcOXW8UFxaI7BGG3gege4ovMCrhi3gWa4dXxQHJpjMEZpNWbMSKDRBjov\n4IpxG66c2jtroThsERux4YjrHeCKzgu4YtyGI7f2zlooDlvCRmy4cpsxAyaIbbgituHIrb2zForD\nljADjTawOg1XxDZcEdtAM5j0aAbFYUsIWJTG6jRcEdtwRWzDGQd6raZeFofM0sERq9NwRWzDFbEN\nVxzotbp6VxwySwdXrE6jDZwEB1fENtAcJj5WV++KQ4IVbaGdAm44CQ6uiG04Y+IDi+hdcUiwog20\nU8ARJ8HBFbENV0x8YFG9Kw4JVrSBFWo4YnINrohtuGLiA4vqXXEoEawoj0QDbaB1Ga6IbTiivROr\nIDKz9jMUMxgMcn19vfZjoLIag3HN70U/1DxcCyiJ2IajmnFNPoKIGGbmYJ739nLlEP3B6bRwResy\nXBHbcER7J1bFgdoPAJQ0azAGmjbc3NKJkxsabm619p2TVqG1EK1CKIbYhqu2Y5u4xqpg5RCtodce\njjgJDq6IbbiqEdvENdqybL5NcYhWkGTAFa1CcEVsw1Wt2CauUVoT+TZtpWhFzfbOI+cd0vFjhxmQ\nUQStQnBFbMMVsQ1XTeTbrByiFbR3oi0cgQ9HNdryiW20gdgGmtNEvs1VFj1VYzDmKGWUxhH4cERc\nwxWxDWddukqNqyywq5r7/xj0URJH4MMRcQ1XxDZcrfJVauw57CGud0AbOAIfjohrOON6B6AZq5xr\ns3LYQ+z/Q2mcTgtHxDWccb0DXHGV2mIoDitjIzYccQQ+HBHXcMb1DnDEpN7iKA4rWuV+ZGA3qzxj\nBuyEuIYz4huOmNRbHMVhRWzERhtYnYYrri2BK8ZtOKK9czVQHE5pO2gJWJTG6jRcceoyXDFuwxHt\nnauD4nCMjdhwxOo0XBHbcEVswxHtnauD4nCMjdhwxOo02kCrEFwR20AziOvVEZlZ+xmKGQwGub6+\nPtd7JyuHk6Bts40D/VEj0ajxneiPmi1wxDZKIrbhjHykXyJimJmDed7LyuEYLZ4ojT1ScESrEFwR\n23BFPoLdHKj9AF1y5LxDOn7sMIGLImYlGsCqm7QKrYVoFYIVYhuuyEewG1YO0UvsI4ErrneAI652\ngCvyEXQNew7RO+wjgauasQ2UQlzDFfkI2sKeQ2AX7COBK47AhyPiGq7IR9BF7DlE77CPBG0Ybm7p\nxMkNDTe3WvtOYhttaDu2iWu4IrbRRbSVojqOU4YbWoXgqlZsE9corVaMEdtoA22lWBkcpwxHtArB\nVa3YJq5RUs0JPWIbpQ03t7R21r3Pmff9tJWiKo5ThiNahdAGWpeBZpCLwNVk4mPtzEP3m/czrByi\nKo5TRhu43gFuanZdENsore0xm1wEriYTH4r5P0NxiNtxjxQc0boMR7Quw1WNMZtcBG2oeaelFjhk\nhuIQkui3hy+OwYcjVjrgin2tcFS72+Ohl338w/N+huIQkkig4YskGm2gdRmOaq50MGbDSe1uj9uu\nv+6j876f4hCSGIzRDlqX4YjWZTiqvdLBmI2S2Ne6M4rDjmIWGm5oXYYrOi/gqPZKB/8fQinsa90d\nxWEHMQsNRyTQaANtcHDFSgfQDPa17o7isINIolEaCTQc0QYHV6x0wBX5SPdQHHYQQYuSSKDhijY4\nuGKlA47IR7qJ4nAPHKABNyTQcMXEGlwR23BEPtJNFIe74AANOCLJQFs4WAuOmDQGmkE+0k0Uh7tg\n7x/aQAINRxysBUdMGsMVkx6YoDjcBTMaKI0EGq6YXIMj4hqOmPTANIrDXTCjgdJINNAGToODI+Ia\naAa5CKZRHO6BGQ2URKKB0jgNDo6IazjjTkvURHEITGH/H9xwGhwcEddwxZ2WqI3iEBhj/x8cMSMM\nR8Q1XHGnJWqjOATG6LlHaZwGB1d0XcAR+1rRRxSH6CQGZLjhNDi4ousCjtjXir6iOETnMCDDESvT\ncEVswxH7WtFXFIfoHAZkOGJlGq6IbTgirtFXFIfoHAZktIE9UnDEvla4YswG2hGZWfsZihkMBrm+\nvl77MVZajUSj5veiH2ru/wNKIa7hitiGq7by3YgYZuZgnveycogdcYAGXLFHCo6Ia7gituGoq5Me\nB2o/ALpr1mAMNG24uaUTJzc03Nxq7TsnrctrIVqXUUzbsU1cow2M2UAzuppns3K4QtputWTvH0rj\nZFq4qhHbxDVKY8yGK65QuwPF4Yog0YAjTqaFq1qxTVyjJMZsOGLS484oDlcEiQYcdXXWDFgWsQ1H\nxDUcMelxZxSHK4IBGW3gqHA44noHOCKugWaQY98ZV1nsA9c7wFFXT80ClkFcwxFxDWc18l33HJur\nLAriege44qhwOCKu4Yi4hqua+//4/9BIlassImItIt4eEa/Z9vovR8T1O3zmbhHx8oj4h4h4T0Rc\n0s7T3llXj50FlsVR4SiNI/DhimtLgGaQZ9dXa+XwaZLeI+mekxciYiDpXrt85vGS7p6ZXxIRZ0h6\nd0T8TmZ+sOiTbkNfMtrAXhK44TQ4uOI0cbjieod+ar04jIj7S3qUpJ+VdPH4tTVJL5T0XZIeu8NH\nU9KZEXGapHtIukXSJ4o/8DYMyCiN1mU44jQ4uOI0cThiQq+/aqwcXibpJySdPfXaj0h6dWZ+JCJ2\n+tyrJH2LpI9IOkPSj2Xmv21/U0RcJOkiSTr33HMbfOw7MCCjJPaSwBGzwXBFbMMRE3r91WpxGBEX\nSro2M4cRcXT82n01ahk9usfHv1rSbZLuK+mQpL+KiDdk5jXTb8rMl0h6iTQ6rbTRfwD0Di0VcMW1\nJXBESz4ckYugTa1eZRERL5D0vZJulXS6RnsObx7/cdP4bedKuiYzD2/77AlJb87M3xr//eWS/jwz\nX7nT95W6ygL9ULO90/1IZdTFMfhwRFzDEbkImrDIVRatnlaamZdk5v0z8wGSniDpysw8lJnnZOYD\nxq/fuL0wHPv/JD08Rs6U9DBJ/9jaw6N3ap6YdeS8Qzp+7DCDMYrgNDg4Iq7hiFwEbatylcW8IuLR\nEfHc8d+ekHSWpHdKequkl2Xm31d7ONjjqHC0gesd4IrrHeCIuIa7VttK20ZbqZca7Q20VKAk2oXg\nqlZsE9coibjGqlqkrbTWPYfAQmoeqcxAjFI4DQ6uuN4Bjohr9EGn20qBCfaSwBHtQmgDrctAM4hr\n9AErh1gJHKmMNnC9A9xwkTVccW0JUAbFIfaFJBpuaF2GI1qX4ajmfm3iGqXV3mNKcYiFkUTDUc0k\nGiiFrgs4YryGqy7c10pxiIUxKMMRSTTaQNcFHLUd14zXcNWFHJviEAtjUEYbSKLhhq4LOKoR14zX\naEON9s4u5NgUhyuODdlwRBINR12YEQaaxvUOcNTnw7woDlcYG7LhiiQapfV1RhhoGnENR30+zIvi\ncIWRQKMNJNFw0+cZYXijmwiOyEPaRXHYIDZkww1JNBz1eUYYvugmgiPykPZRHDaEDdlwRBINR0ys\nwRHdRHBEHtI+isOGsCEbjkiiURptcHBEGxzQDOK6fZGZtZ+hmMFgkOvr661812TlcBK8NS6thLca\nyUbN74W/Llz2CzStZlwzXqO0GjFGXC8vIoaZOZjnvdYrhzfecltr38VMNEpiLwkc0QYHR7TBwRXX\nTPXDgdoPUNI1H7tew82t1r7vyHmHdPzYYQIYjZuVbABNG25u6cTJjdbGzUm70FqIdiEUQ1wDzSAX\n6QfrlcNMMRONIjiZFm44VAuOiGu4Yl8rSrEuDiNE4KJxJBtwxKFacERcwxHXO6Ak6+Lwgfc5i8BF\n40g24IgZYTgiruGIfa0oybo4POPgWu1HgCGSDbSh7ZYhZoRRGteWAM0gD0FJXGWBlcexynDDFQ9w\nQ0zDFddMYRVwlQV6g2OV4YgrHuCGmIYjrpmCI+urLOCPY5XhiKPwURrXOwDLIweBI1YOsdLou0dp\n7JOCG05chiuumQKWR3GIxpBEww0tQ3DEictwxKQHXLWdX1McohEk0XDEPik4YrUDjpj0gKMa+TXF\nIRpBEo3SaqxMk0SjDVxbAke0eALLq5FfUxwaIomGm5qn0pJEoyROXIYjWjzhqg+THhSHZkii4ajm\nyjRJNEqi6wKOaPGEo75MelAcmiGJhiNWptEGui7giLgGmtGXSQ+KQzMMyGgDe6Tghq4LOCKugeb0\nJcemOCyMJBpu2CMFR3RdwBFxDVdcn1YOxWFBJNFwxB4pOOrLjDD6hbiGI65PK4visCCSaDgi2UBp\nzAjDEXENNIP8uiyKw4JIotEGWpfhhBlhOCKu4YoDj/z0qjgkiYYbWpfhhhlhOCKu4YgDjzz1pjgk\niYYjEg6U1ocLf9E/xDWwPA488tSb4pAkGqXRWgE3fbnwF/1CXMMVkx5oQm+KQwIYJdFaAUd9ufAX\n/UJcwxGTHmhKb4pDAhgl0VoBR0yqwRFxDUdMeqApvSkOJQIY5ZBsoA0cqgU3XO8AR2wzwSqLzKz9\nDMUMBoNcX1+v/RhoWY1Bueb3oh9qHoUPlEBMw1HNuCYPwU4iYpiZg3ne26uVQ/jjLim44lAtuCGm\n4YhtJlh1B2o/ANCkWYMy4GDSMrQWomUIFohpOCKusepYOYQVeu7RBvZJwRF7WuGIuAYWw55DFFUj\niabnHiWxTwqOiGs4Iq7haD95LnsO0Qk17/5j8Ecp7JOCI+IajohruGkjt2bPIYph/x9KG25u6cTJ\nDQ03t1r7TvaToDTiGo6Ia2B5beTWrBz2SNvtluz/Q0k1V6bZT4JSiGs4Iq7hyjG3pjjsiRoDM4My\nSuK4cDgiruGIuIYj19ya4rAnag3MDMoohZVpOCKu4Yi4hiPX3JrisCcYmFEax4XDEXENR8Q1sDzX\n3JqrLCqoddUCVzygFI4LhyPiGo6Ia7ji+rSdcZVFh9UclGnxRCkcFw5HxDUcEddwxPVpzeEqi5Zx\nvQPa0PaR4RwXjtI4Bh+OiGugGeTXzWHlsGWu/cnoDtfTs9BfHIMPR8Q1XNVotSS/bk7vi0M2ZcON\n6+lZ6C+OwYcj4hqOmPRYfb0uDulPhiNmz+CGmIYj4hqOmPRYfb0uDtmUjdJqtFYwe4bS6LiAI+Ia\nbmjvxH70ujgkgFESJ9PCER0XcERcww3tndivXheHBDBKYmUajohrOCKu4Yb2TuxXr4tDiQBGOaxM\nozRahuD+74XYAAAgAElEQVSIuAaWR0xjvyIzaz9DMYPBINfX12s/BjqiRsJR4zvRDzXblolrlEJc\nwxU5CGqKiGFmDuZ5b+9XDtEP7CeBG1qG4Ii4hiNyEKySA7UfAGjDrIQDWGWTlqG1EC1DKGa4uaUT\nJzc03Nxq5fuIazgiB8EqYeUQvUDvPUrjGHy4qbHaQVyjNPa0ArujOEQVJNJwQssQHNVq8SSuUQrX\nO8BVk3k1xSFaRyINNxyDD0esdsANe1rhqOm8muIQrSORRmltr0yTRKO0Gq1wrHagNMZqYHlN59UU\nhz1H7z3csE8Kbmpe78BqB0phrIYjh7ya4rDH6L2HI/ZJwQ3dFnDEWA03Lnk1xWGP0XsPR6xMozRa\n4YDlEddw45JXUxz2GAMzSmOfFNzQCgdXnCION0zk7U9kZu1nKGYwGOT6+nrtx5hLjSS65vfCX819\nUkApJ05u6NIrrtaplNZCuviC83X82OHajwUshfEabmrFdFfz6ogYZuZgnveyctgBHDYAR+yTgiOX\nmWFgGuM13LCndf8oDjuAQRmOSKLRBlrh4MbhtEOgNmJ6/2gr7YDJyuEkgGnnQNNoW4YjWuHgpmZM\nM16jpBrxRUzfgbbSJTETDSe0LcMVXRdw43LaITCt5hUPxPTiKA63IYDhhgQabaAVDo447RBYHnnI\naqE43IYARmkkG3DjcvEvMI1rS+CIiTzsheJwGwIYJZFswBGtcHDEaYdww0Qe5kFxuA0BjJJINuCI\nSTU4Iq7hhok8zIPicAYCGKWQbKC0Gi1DTKqhNOIaWB45CObBVRboNY5WhhOudoAj4hqOuGIKbeIq\nC2AOnEwLNxyoBUfENdxwxRS67EDtBwBqmZVwAE0abm7pxMkNDTe3Wvm+ScvQWoiWIRRDXAPLIf9A\nl7FyiN6i9x4lcTItHBHXcMQVU8AdKA7RCRw2ADecTAtHxDXcMOEBV9O59SIoDlEdvfdwxMwwHBHX\ncMOEBxxtz63j4D3OnPezFIeojsMG0Ia2V6eZGUZpdFzAES2ewPK259YHDp5x9ryfpTjEXTAwww0n\n08INHRdwRIsnHNWYyNueW5+65cZPzvtZikPcCQMzHLE6DTfENBzR4gk3NSenp3PrwfM+dcO8n6U4\nxJ0wMMMRq9MojY4LuOnCagdxjVVXcyJvv7k1xSHuhIEZpbFPCm7ouICbrqx2ENdYdauYV1McdhyH\naMAJ+6TgiI4LuFnF1Q5gHuTVe6M47DAO0YAb9knB0SrODAO7IabhiLx6PhSHHUYijdLYJwU3tC3D\nEasdwPLIq+dDcdhhJNIoiX1ScEPbMhyx2gFHHHjUXRSHc2I2Gm7YJwU3zArDEXENNxx41G0Uh3Ng\nNhqOmEGDG2IajohruOHAo26rUhxGxJqkdUn/nJkXTr3+y5KelJln7fC5L5X0Ykn3lHRK0ldl5k2l\nn5dZO7SBPSVwQ0zDDV1EcMT5A5hWa+XwaZLeo1GRJ0mKiIGke+30gYg4TdJvS/rezHxHRNxb0qdL\nP6hEEKM89pTADTENN3QRwRHnD2C7A21/YUTcX9KjJL106rU1SS+U9BO7fPQCSX+fme+QpMy8LjNv\nK/msE5MgvviC81v9jwH6Y9bqNLDKiGm4IabhqFZcHznvkI4fO0xO3UE1Vg4v06gIPHvqtR+R9OrM\n/EhE7PS5z5eUEfF6SZ8l6Xcz8xe2vykiLpJ0kSSde+65jT00s3YoidVplEbbENwQ03DDCZ7ogsjM\n9r4s4kJJj8zMp0bEUUlP16iQe6Wko5l5a0RcP2vPYUQ8XdJxSV8l6UZJb5T0U5n5xp2+bzAY5Pr6\neoF/EjirMTjX/F74q9UOR0yjFGIabmq2LRPX/iJimJmDed7b9srh10h6dEQ8UtLpGu05fJekmyVt\njFcNz4iIjcw8vO2zH5L0psz8mCRFxOskfaVGRSLQCPaUwBHXlsANMQ03nOCJrmh1z2FmXpKZ98/M\nB0h6gqQrM/NQZp6TmQ8Yv37jjMJQkl4v6Usj4ozx4TRfL+ndrT08eoE9JXA0aRtaC9E2hMYNN7d0\n4uSGhptbrX0nMQ03xDS6otP3HEbEoyUNMvOZmbkVEf9V0lslpaTXZeZr6z4h3NB7j9I4Ch9OuMwa\nrriKB33V6p7DtrHncPXVSKTpvUcpNduWgRJOnNzQpVdcrVMprYV08QXn6/ixWc0/wOpgrIabLu85\nBObGPWlwU3NPCVAC3RZwxFgNN3HwHmfO+16KQ3QWgzNK4yh8uKEVDo4Yq4HlHDh4xtl7v2uE4hBz\nY3CGkxor0yTSKIluCzhirAaWd+qWGz8573spDjEXBme44Sh8uKHbAo4Yq+Go7QWXvOVTN8z7XopD\nzIXBGW5YmUZpdFsAyyOu4abrBx5RHGIuDM4oiesd4IZuCzhirAaW1/UuD4rDFcTgDCc1Z9BYmUYp\ndFvADWM1HNXIqbu+4EJxuGIYnOGm6zNowH50/T/+wKIYq+Gm5iFeXV5woThcMQzOKI19UnBDtwXc\nsNoBLK9mTt3lBReKwxXD4IyS2CcFN3RbwA2rHXDF5HQ3UBwuiQuH4YR9UnBDtwXcsNoBR0xOdwfF\n4RK4cBhumEVDacwMww0xDSyPyenuoDhcAjPSKIl9UnDDzDDcENNwxaRHf1EcLoFARinsk4IjZobh\nhpiGIyY9+o3icAkEMkphVRqOmFCDG2Iajpj06DeKwyURyCiBhANt4EAtuCGmgeWRg/RbZGbtZyhm\nMBjk+vp67ceAgRr7/2p8J/qjZusyUAIxDUe1cgFyEC8RMczMwTzvZeUQ2AOn0sIRrctwQ0zDDecP\noIYDtR8A6LpZCQfQpOHmlk6c3NBwc6u175y0Da2FaBtC44hpYHnkH6hhXyuHEXFA0hdKurek9cy8\nodGnAnbB8cpwUnNlmr1SKIGYhqMabZbkH6hh4eIwIo5LepZGhaEkfZWkt0XEH0u6MjP/W4PPB9wJ\nxyvDTc1WONqGUAIxDTdMeKBPFioOI+IHJP2SpMslXSHplVM//itJj5NEcYhiOF4ZbpgZhhtiGm6Y\n8ECfLLpyeLGkSzPzJyNibdvP/lHSM5p5LGA2kg6UxlH4cENMww3bS4ByFrrKIiJukvTIzLxyXBx+\nWtIgM98WEUcl/Xlmnl7mURfHVRZlcbwy3HAUPtwQ03BTK6bJPVBayRgreZXFxyQ9YIefnS/pnxf8\nfVhRHK8MRxyFDzfENNywvQSOujSRt+hVFn8q6ZkR8cCp1zIi7iPpxyT9cWNPhk7jeGW0oe3j8DkK\nH6UR03BDTAPL61JevejK4U9Jerikd0r6O0mp0QE0D5F0raTnNvp06Cz671EaJ9PCDTENN8Q0HPX9\n2pKFisPMvC4iBpJ+VNIjJL1//Dt+RdIvZuYnmn9EzIMDB+CG1iG4IabhhpiGG64t2cc9h5n5SUk/\nM/4DHVAzkBmcUUqXZtGAJhDTcENMww3XluyjOET3cOAASqvRYtGlWTT4IabhhpgGlseExxzFYURc\nucDvy8z8hiWeB/tAIKMkTqaFG2IabohpOGLCo455Vg4PaHTwzMT5ks6R9EFJ/yLpszW63uIjkq5u\n9vFWE/v/4ISVabghpuGGmIYbJjzq2bM4zMyjk7+OiMdI+iVJD8vMt0y9/lBJvzf+Wa+x/w9uWJlG\naW1PqBHTKI2YBpbDhEc9i+45/BlJPz1dGEpSZv5dRDxb0vMk/UlDz7aSCGaUxso0nHAUPtwQ03DE\nhEd/LFocPljSv+7ws2slHV7ucVYfwYySWJmGG47ChxtiGm6Y8OiXRYvDD0j6QUl/NuNnP6jRPsRO\nYZUFTliZhhsm1OCGmIYbJjz6ZdHi8DmSXhER75T0Kt1xIM23SXqIpO9u9vGWc+Mtt7HKAiskHSiN\nCTW4IabhhhZPlBSZufe7pj8Q8Y0aFYkDSXeT9GlJb5X0rMx8Y+NPuIRzz/+SPO1bf06nUloL6eIL\nztfxY73vfEWDahyzXOM70Q81T4cDSiCm4aZWTJN7rLaIGGbmYJ73LrpyqMx8g6Q3RMQBSfeR9LHM\nPLXo72nDWXc/TaeY6UAh7P+DG9qW4YaYhhtaPFHawsXhxLggvLbBZ2ncGQfX9GJaO1AISQdKqjFL\nS+sQSiKm4YaYhqOF2koj4pl7vCUz82eWe6TmDAaDXF9fr/0YaEGt9s7vfumbbx+gaVdCU2q2wtE6\nhBKIabghprFKSraVPnuXn02qzM4Uh+iHmu2dHDqAEmquStM6hBKIabghpuFqoeIwMw9sfy0iDkn6\nZkk/LukxDT0XMDcGaLihbQhuiGm4Iabhat97Dicyc0vSb0bEvSWdkPTIpZ8KWAADNErjKHy4Iabh\nhpgGmrHwVRY7/qKIh0t6dWae1cgvbAB7Duvgegc44Sh8uCGm4YaYhqMmc9uiV1ns4kJJ/9rg78MK\n4noHuOFUWrghpuGGmIabmhMeCxWHEXH5jJcPSvpiSV8i6VlNPBRWFwM03NC2jNLa7nwgpuGGmIab\nmvn0oiuHD9cdp5JO3CRpU9Jlkl7exENhdTFAozT2lcBJjdlhYhqlMU4Dy6mZTze257CL+r7nsNY+\nPPb/oRT2lcDNiZMbuvSKq3UqpbWQLr7gfB0/drj2YwH7xjgNR6t+nkaxPYcR8X2SXpuZ18342WdK\nujAzf3OR34kyag7O7P9DKbQtww3dFnDDOA03fTtP4y73Fu7hZZIetMPPPm/8c3TArMEZaNJwc0sn\nTm5ouLnV2ndOEum1EIk0GlcjpiftcBdfcD4rLGgc4zSwvL7l1IvuOYxdfnampFuXeBZrHDgAJzVn\n0dhXghLotoAbxmk4qtHe2beces/iMCK+XNJXTr30zRHxxdvedg9JT5D0vgafzQYHDsBNzbYhEmmU\nQCsc3DBOww0THu2YZ+XwW3THFRUp6b/s8L7rJD25iYdyU2uAZnBGKX2bRUP76LaAG2IaWA4THu2Y\npzi8TNJvaNRSeo2kb5X09m3vuVnSv6Tz0adLYIBGaRwbDid0W8ANMQ1HTHh42rM4zMx/l/TvkhQR\nnyfpI5l5S+kHK4lEGk76dooW/NFtATfENNww4eFroQNpMnOz1IO0hUQabtgrBTfMDsMNMQ03THj4\nmudAmtsk/e+Z+ZaIOKXRvsOdZGYuegJqq0ikURKnaMFNjZhmdhglEdNwQ+6BJs1TyD1X0oem/nql\n9xUSzCiFU7Tghusd4IaYhhtyDzRtnj2Hz5n662cXfZoWEMwohVO04IZOC7ghpuGG3ANNO7DImyPi\n8vGhNLN+dl5EXN7MY5V15LxDOn7sMAGNRk1WpddCrEqjiOHmlk6c3NBwc6uV7yOmURoxDSyHmEbT\nYpHbJ8Z7Dh+WmW+Z8bMjkt6SmWsNPt9SBoNBrq+v134MVFKjB7/Gd6IfarUOEdMohZiGI3IPdFFE\nDDNzMM9793N4zE7V5DmSPrWP3wc0jlNp4YaT4eCGmIYbcg84mOe00sdKeuzUS8+JiI9te9s9JP1H\nScMGnw3YN/aVwA2HaaE0LrQGlkPuAQfzrByeq1HhJ41WDb9c0s3b3nOzpL+VdElzjwbsH0kHSms7\nkeYwLZTEhdZwxIQHsLhF9xx+QNJjMvMd5R6pOew57A568OGk5nH4QAknTm7o0iuu1qmU1kK6+ILz\ndfzY4dqPBewbe1rhZpnYKrbnMDNnnlQK7IYefLihdQhuWPGAG/a0wkmbufR+DqRRRByS9GBJp2//\nWWb+5bIPBS8k0iipxiwtiTRKqhHTtHiiJMZpYDlt5tILFYcRcbqkyyV9u6TY4W2ducoC3cAAjVJq\nrkqTSKOEmi3LrHigBMZpOHLez7royuFPSzoq6YmSfkvScUk3Sfp+SZ8j6WkNPhsK4SANuKi5Kk0i\njRLotIAbxmm4cT/Aa9Hi8HGSnivpdzUqDv8uM98m6WUR8fuS/pOkP2v2EdEk9v/BCavSKM15dhho\nAzENN+77WRctDs+V9K7MvC0iPi3pzKmfXS7pZWL1sNOYlUZpbSbTrEqjJPfZYfQT3UPActwnPBYt\nDq+TdNb4r/9J0pdJ+qvx399H0j0aeq7eYFYaTmol0yQbKMF9dhj9Q/cQHDHh0axFi8M3S/oKjVpH\n/0DSz0TE2ZJulfTjkv662cfzxqw03LAyDSdMpsENYzTcMOHRvEWLw5/XqLVUkp4n6bBGexDXJP1P\nST/c3KP5Y1YabkimURKzw3BD9xCwHCY8mrdQcZiZ65LWx3/9SUmPi4i7S7q7pK+W9BpJX9r0Q7pi\nkEZJ3JUGJ8wOww3dQ3DDfZYeFl05vIvMvFnSzRHxGZK+aPlHqoNEGk64Kw1umB2GG7qH4IT7LH0s\nXRw6IJGGGxJpuGF2GG6IaTjhPksfFIcikYYfkg6URKcF3BDTcMN+VuwXxaEIaJTHQRpwQacF3BDT\ncMN+Vixjz+IwIh445+86Z8lnqYaARkkcpAEndFrADTENN+xnxTLmWTnckJRzvC/mfF8nEdAohcQD\nJdE6BDfENLAcYhrLmKc4fFLxpwBaROIBF7QOwQ0xDTfsZ8Wq2bM4zMyXt/EgQBtIPOCE1iG4Iabh\nhP2sWEUHaj8A0KZZiUcbjpx3SMePHWagRqMmq9JrIValYYGYhpNaOQewDE4rRa/Q4omSOJUWbohp\nOGFbCbC3yFzZM2T2NBgMcn19vfZjYAc1+vBrfi+81WwfAkogpuGkVjyTc6CkeeMrIoaZOZjnd7Jy\niCrow4cbTqWFG2IaTtjPCjelcmn2HKIK+vBR2nBzSydObmi4udXK97FXCiW1Hc8SMY2yGKOB5ZTK\npVk5hCT68OGFU2nhpFanBTGNUhij4cgll6Y4BIM07NA+BCc12zuJaZTAGA03Trk0xSEYpGGHlWk4\nIZ7hhpiGG6dcmuIQDNIoqsZJbaxMoySud4ATxmhgeU65NFdZdFCNgZqjllECR+HDDTENJ8QzHHFV\n2l1xlcUKq3nwQNcCGauPo/DhhpiGE+IZbrgqbXlcZdExXPGAkjg6HG6IaTghnoHlkEcvj5XDjnHq\nWUa3OJ2kBUjENLwQz3Dkcr1Dn1Ac7oGDB+DC6SQtQCKm4YV4hhsmPFYTxeEu2P8HJ8ymwQ0xDSfE\nM9ww4bGaqhSHEbEmaV3SP2fmhVOv/7KkJ2XmWbt89lxJ75b07Mx8UcnnZKM2SmJVGk44Dh9uGKPh\nhhZPzKPWyuHTJL1H0j0nL0TEQNK95vjsL0r6s0LPdScENUphVRpOOB0Obhij4YYWT8yr9dNKI+L+\nkh4l6aVTr61JeqGkn9jjs4+RdI2kd5V8xolJUF98wfnc/YNGcZoWnBDPcENMw02tmD5y3iEdP3aY\nHHqF1Fg5vEyjIvDsqdd+RNKrM/MjETHzQxFxpqSflPRNkp6+0y+PiIskXSRJ55577tIPyyweSmBV\nGiXROgQ3xDSwHGIa84rMbO/LIi6U9MjMfGpEHNWoyLtI0islHc3MWyPi+ll7DiPiRZLekpmvjIhn\nS7p+rz2Hg8Eg19fXG//ngJcae6Vqfi+81WqHI55RCjENN+QdaFtEDDNzMM972145/BpJj46IR0o6\nXaM9h++SdLOkjfGq4RkRsZGZh7d99qGSvi0ifkGjvYmnIuKmzPyV9h4fbtgrBTecDgc3xDSckHeg\n61rdc5iZl2Tm/TPzAZKeIOnKzDyUmedk5gPGr984ozBUZv7HqfdcJun5FIZYFvtK4GbSOrQWonUI\nRQw3t3Ti5IaGm1utfB8xDSfkHei6Tt9zGBGPljTIzGfWfhZ4ogcfpXEcPpxw4iHcsJ8VuLNW9xy2\njT2Hq6ntgZoefJRSs30IKOHEyQ1desXVOpXSWkgXX3C+jh+7S7MPsBLYzwpHs+Kry3sOgV3VmpVm\ncEYJtfZKAaWw6gEn7GeFmybyaIpDdArJNEqifQhOaqw+0OKJkhijgeU0kUdTHGJHNRIPBmqUwl4p\nOOHEQ7hhjIabVc2jKQ4xU63Eg4EapdA+BCd0WcANYzScrHIeTXGImWomHgzUKIFVaZREOxywHGIa\nTlY5j6Y4xEwM0iiJvVJwQjsc3DBGww0TePPjKosVUWOg5qhllMD1DnDD9Q5wwhgNN1xZwlUWdmr2\nLdcOZvhhrxTcrPIMMbAdYzTcsJ91MRSHK4CBGk5IpFES7XBwQzscsBxiejG0la6AycrhJKhp8UCT\naFmGC9rh4IZ2OLipFVt9j2naSgtrO8CYlUYptCzDCV0WcEM7HJxwP+tqoDhcEMk0nJBMo7Q2J9No\nHUJptHgC+0fOsRooDhdEYKMkEg84aXsyjS4LlMSVJXBDzoFZKA4XRGCjFBIPuKkxmUaXBUqhxRNO\nyDmwE4rDBRHYKIXEA26YTIMT4hlOyDmwE4rDfSCwUQKJB0riige44XA4YP/IObATrrIAdsAVD3DB\nFQ9wQ0zDDTkHSuIqC2BJnEoLJxykBTfENJyQc6BLDtR+AKCLZiUeQFOGm1s6cXJDw82tVr5v0j60\nFqJ9CI1rO54lYhpeyDnQJawcAjPQi49SOCEOTmqueBDTKKFGqyU5B7qE4hCdx0EacMIJcXBSs72T\nmEbTmOyAo+HmltbOuvc5876f4hCdVvPQARIPlMAMMZwQz3DCZAfcTPLotTMP3W/ez1AcotM4dACl\ncRw+nBDPcNNmTDPZATeTPFox/2coDrEwBmq44IQ4OCGe4abtmGayAyXV3M+qBe4upDjEQhio4YSV\naTghnuGmRkwz2YESau9nfehlH//wvJ/hKgsspMZxy0fOO6Tjxw4zWKNxHIcPJ8Qz3BDTcFHzupIj\n5x3Sbddf99F538/KIRZCmydKYr8UnBDPcMLJ4cD+rVL+HAu0oK6cwWCQ6+vrtR+jmBoDdc3vhbea\nJ9MCTSOe4YR4hpsauWzN/Dkihpk5mOe9rByuKK54gBv2S8EJ8QwnxDOccHjX7thzuKJq9i7D33Bz\nSydObmi4udXad7K3BKUQz3DTdkwTz3BCDr07Vg5X1Cr1LmO11D5Ri5ZlNIl4hpsaMU08o6S22y3J\noXdHcdgQDh6Ai5rtQ6vScoHVQTzDTa2YJp5RApMd3UNx2AB6l+GEGTWU1uZkGvEMN8Q0nDDZ0T0U\nhw1gozZKYlUaTtqeTCOeURpjNJzQ4gnr4vDGW25r5XsIbJTCqjTc1JhMI55RCmM0nNDiCcm8OLzm\nY9druLlFYGNlsSoNN0ymwQljNJzQ4gnJvDjMFIGNlUYijZJqXMjLZBpKqRHPjNEoiRZP1BCZWfsZ\nijn9vg/Ov/mfbyH5QGNqJB81vhP+arXDASXUjGfGaJRQK6aJZ08RMczMwTzvtV45fOB9ziKw0Rj2\nlsAJ7XBwwpUlcEOLJ2o5UPsBSjrj4FrtR4CRWQM1sKom7UNrIdqH0Ljh5pZOnNzQcHOrle8jnuGG\nmEYt1iuH8EYvPpxwHD5ccOIh3LA/G31CcYiVRPIBJ7QswwntcHBScz8rMY0arNtK4atWi+eR8w7p\n+LHDDNZoFC3LcEI7HJwwPqNvWDnESqLFEyXRsgwntCzDCeMzUJb1VRaDwSDX19drP0YvcMUDXHB8\nOJxwZQmcMD4D+8NVFmgV+6XghP1ScMKVJXDC+Aw3XZx4oDjE0kg+4IQWIpRESxywf8QznHS1s4Pi\nEEtjsEZJ7JeCC05ZhhvGZzhpO567urhCcWiIwRouaFmGE1ri4ITxGU5qxHNXF1coDs0wWMNJV2fV\nsPpq7PPoaiIA7AfjM5zUiOeuLq5QHJphsEYpJNNwUXMSrYuJADywnxXYv1rx3MXFFYpDMwzWKIFk\nGk5qTqJ1MRHA6mM/K9ywRaoeisOCaqy0ENwogWQaTphEgxv2s8IJW6TqojgspObxtAQ3mkYyjZKY\nIYYTWvCB5bBFqi6Kw0IIbJTWZgJCMo1SmCGGE1rw4YbJjv7pTXHIRm04qbW/hIQDTWMiDU5owYcT\nJjv6qRfFIRu14YaEGi6YSIMT4hlOmOzop14Uh2zUhhsSEJTC/j84IZ6B/SPX6KfIzNrPUMxgMMj1\n9fXbVw4nwd3m4TDwV6Mfv+b3wlfNg7SAphHPcFPjv/vkGh4iYpiZg3ne24uVQ2byUAqn0sIJ7cpw\nQjzDCYd3oS0Haj9AW46cd0jHjx0mwNGoWckH0JTh5pZOnNzQcHOrle+btBCthWghQuOIZ2D/yDfQ\nll6sHKI/OJUWLjhIC06IZ7gh34ArikPYIPmAEw7SghPiGU7IN+CM4hA2SD7ghFliOCGe4YR8A84o\nDmGD5AMlcSQ+nBDPcEKLJ9CcXlxlgTo4chkuOBIfTohnOKkVz+QbKKVEbHGVBarjyGU44Uh8OCGe\n4YQWTzjpwuRdb66yQLs4chklcSQ+nBDPcNF2LEvEM7x0IX9m5bAn6MeHC06JgxPiGS5qdgwRzyih\nRutwF/JnisMeIPmAE1qI4IR4houa7crEM5rW58kOisMeIPmAky7MqgFNIZ7hgliGkz5PdlAc9gAD\nNkriSHy4qNFCRDyjFMZmYP/6nDtzlUUFXPEAF104VQtoArEMJ8Qz3JA7L4erLDqMKx7ghCPx4YJY\nhhPiGU7IndvFVRYt68IRtUBTOEIcpXC9A5wQz8D+kTu3i5XDlvW5hxllsV8KLjhhGU6IZ7jhejRv\nvS4OSabhoub+kr62XaAcTliGE+IZTpjs8Nfb4pBkGk7YXwInzBLDCfEMJ0x2+OttcUgyjZJouYAT\njsSHE+IZTsg30LTeXmUxWTmcBDfHPKMptValnY5cRndwJD6cEM9wQr6BeXGVxRyYyUMptFzACV0W\ncEI8wwn5BkrobXEoEdwog5YLlEQLEZwQz8D+Ec8oobdtpeiPGu0PtFygBFqI4IR4hhvyDXQVbaXA\nWK3kg1VplEALEZwQz3BCvgEXB2o/AFDSrOQDWFWTFqK1EC1EaNxwc0snTm5ouLnVyvcRz3BCvgEX\nrBzCGv34KIkj8eGCi63hpEarJfkGXFAcolUk03BBCxGc0OIJFzXHZvINlNB27kxxiNaQTMMJR+LD\nCe/DF+0AABpPSURBVKsecFFzbCbfQNNq5M4Uh2gNyTSckEyjlBotcax6oBSuKwH2r0buTHGI1jBg\noxSSabio1WEhseqB5rGXFW76MNlBcdhj7P+DA5JpOKHDAk7YywonfZnsoDjsKfb/wQXJNErqwywx\nUArxDCd9meygOOwpEmqU1GZCTfKBUvoyS4z+oGMITpi8K4PisCMIcLhoO6Em+UApfZklRj/QMQQn\nTN6VQ3HYAQQ4nNRIqEk+UAKTaHBCxxCcMHlXDsVhBxDgcEJCjVJoiYMTOoaA/SOey4nMrP0MxQwG\ng1xfX6/9GHuarBxOArzN0xbhrcYVDzW/F75qnkoLNK1WPDM2o5QasUU8zy8ihpk5mOe9rBzOwOw0\nHHDFA5zQEgcndAzBCftZvVAcbkOAwwXJNEqiJQ5OiGdg/8g3vFAcbkOAoySueIADDtGCE+IZTmq0\nWpJveKlSHEbEmqR1Sf+cmRdOvf7Lkp6UmWfN+Mw3Sfo5SQcl3SLpGZl5ZdPPRoCjFK54gAta4uCE\neIaLmt1v5Bs+aq0cPk3SeyTdc/JCRAwk3WuXz3xM0jdn5ocj4oslvV7S/Zp+MAIcpXDFA1wwiQYn\nxDNc1Ox+I9/w0XpxGBH3l/QoST8r6eLxa2uSXijpuyQ9dtbnMvPtU3/7LkmnR8TdM/Pmpp+RAEcJ\nJCAohUO04KJGSxzxDBfkGWhC61dZRMSrJL1A0tmSnp6ZF0bE0yQdyMxfjIjrZ7WVbvsd3ybphzLz\nG2f87CJJF0nSueeee2Rzc7P5fwisPK54gAuueIALYhluuN4BXdHZqywi4kJJ12bmMCKOjl+7r6TH\nSzo65+/4Ikk/L+mCWT/PzJdIeok0uudw+aeGG654gBMO0YILYhlOOP0eq+pAy9/3NZIeHREflPS7\nkh6uUYvoYUkb49fPiIiNWR8et6T+kaTvy8z3t/LEsDMrAQGaMtzc0omTGxpubrXyfZM2orUQbURo\nXJvxTCzDCbkGVlWrK4eZeYmkSyRpvHL49OnTSsevX5+Zh7d/NiLuJem1ki7JzL9p4XFhip58lMKR\n+HDCCctwwl2WwHw6fc9hRDxa0iAznynpRzRaYfzpiPjp8VsuyMxrqz0gGsFhGnDBkfhwwgnLcMHE\nHZyUzpurFYeZeZWkq2a8ftbUX79a0qvHf/08Sc9r6fHQEnry4YSZYjghnuGCiTu4aCNv7vTKIfxx\nAAFKYlUaLrjiAU5o8QT2p428meIQd8KADResSsMFJyzDCS2ecNNm7txG3kxxiNsxYMMJq9JwQSzD\nCS2ecOJ4cBfFIW7HgA0nrEqjFDos4IR4BvbP8eAuikPcjgEbJbH/Dw7osIAT4hlYjmPuTHHYYSTT\ncMH+P7igwwJOiGc44eCuZlAcdhTJNJywZwouHGeJ0V/EM1xwcFdzKA47imQaJbHHBC7osIAT4hnY\nH/Lm5lAcdhTJNEphjwlc0GEBJ8QznDAJvbooDufEbB5csMcELpgphhPiGS6YhF5tFIdzYDYPTphd\ngwtiGU6IZ7hgEnq1URzOgdk8lMLJWnBChwWcEM9wQYsnFkFxOAeCHCVwshac0GEBJ8QzXNDiiUVR\nHM6BIEcJrEjDCfEMJ8QzXNDiiUVRHM6JIEfTWJFGSbQRwUWN9nviGaUwNqPrIjNrP0Mxg8Eg19fX\naz8GVkjbg3aNpAf+arXEEc9oWs32e+IZTWNsRi0RMczMwTzvZeUQGKvVl89AjabRRgQXNds7iWc0\njbEZq+BA7QcAumLWoA2sokkb0VqINiI0ari5pRMnNzTc3Grl+4hlOCGesQpYOQTG6MtHKRyJDwec\negg3jM3AXVEcopO4/w8uOBIfLmiJgxPGZjhpMm+mOETncP8fnHAkPlzQXQEnjM1w0XTeTHGIzmHA\nRkkcIw4XtMTBCWMzsD9N580Uh9gTAzZcsGcKLmiJgxPGZjhZ9byZ4hC7YsCGE/ZMwQUdFnDC2AwX\nDnkzxSF2xYANJ6xKo5RVnykGaiKe4cIhb6Y4xK4YsFEKJ9LChcNMMTDB2AwnTNwtLjKz9jMUMxgM\ncn19vfZjNKrGoF3jO+Gt5om0QNNOnNzQpVdcrVMprYV08QXn6/ixw7UfC1gYYzOc1IrnLubNETHM\nzME872XlcIVwAAFcsF8KThxmigGJsRleHFo8a6A4XCEM2nBBMo2SuOIBLmiJA/aPeN4f2kpXyGTl\ncBLktHugSW0nIV1su8Dqoy0OLmiJgxu2RtVDW2lLmJ2Gi1oHahDDaBodFnBBSxycsDVqdVAc7hNB\nDick1CiFtji4IJaB/SPPWB0Uh/tEkKMkkhA44IoHuCCW4aRGqyV5xuqgONwnghylkITABW1xcEEs\nw0XNzjfyjNVAcbhPBDlKIQmBCybRUArdFcD+1Ox8I89YDRSHSyDIUQJJCErhEC04oLsC2D9yDOyF\n4hDYRY2+fJIQlMAhWnBBdwXctJlrkGNgLxSHwA5q3pdGEoKmcYgWXLDyASdcJYWuoTgEdkAyjZLY\nMwUXtCsD+0euga6hOMRKaTMJIZlGKeyZggvaleGEKx4AikOskLaTEJJplMKeKbhg1QMuuOIBbqYn\nOxZBcYiVUSMJIZlGCcwUwwWxDBdc8QAn2yc74uA9zpz3sxSHWBkkISiFPVNwQSzDBfuygf3bPtlx\n4OAZZ8/72cjMks9W1WAwyPX19dqPYalGX37N74WvmqfSAk0iluGiViyTY6CEGnE1+f/QZLLjvb/w\nuH88dfONXzDPZ1k5xMK44gFO2DMFF8QyXLAvGy66spd18LxP3TDvZw+UfDB4mjVoA00Zbm7pxMkN\nDTe3Wvm+SSvRWohWIjSKWIYLYhnYn5o585HzDun4scMLF6OsHGJh9OWjFK54gAtiGS6IZbjhWrTd\nURyaaDPQGbRRCq1EcEEswwWxDCdci7Y3ikMDtWb1ViHAsVpWcYYNmIVYhgtiGU64Fm1vFIcGOIQA\npXAsPlwQy3BBLAP7x2TH3igOC+BuHjioecIWyQeaRCzDBbEMJzWueGCyY28Uhw1j4zZcsCINF8Qy\nXBDLcMG1aN3FVRYNq3Vk7X6PqwV2wlHiKIVj8eGCWAb2h2vRuouVw4bR4olS2GcCB3RXwAWxDCds\nicKEdXF44y236cTJDXqZsfLYZwIXHIsPF8QyXDDRgWnWxeE1H7tel15xNb3MWHnsM4ELZovhgliG\nCyY6MM26OMwUyTSKoP0CLmhXhgtiGS7IMVBTZGbtZyjm9Ps+OO/3xMt0t5ZXDuGtVotnjSOf4a3m\naXFAk4hluCDHQAkRMczMwTzvtV45fOB9ztLxC84n0NEo2i/ggnZluCCW4YIcA7VZF4dnHFzT8WOH\naz8GzNB+gVJoJYILYhkuiGX0jXVb6WAwyPX19dqPgcJqtELQfoGm0UoEF8QyXBDLcEFbKXqDKx7g\nglYiuCCW4YJYRh8dqP0AwDJmDdzAKpq0Eq2FaCVCo4abWzpxckPDza1Wvo9YhgtiGS7i4D3OnPe9\nrBxipdGbj1I4Fh8OuNwaThiXgf05cPCMs+d9L8UhGsXADQe0K8MFbXFwwbgMJ23ny6duufGT876X\n4hCNYeCGC47Fhwu6K+CCcRkuauTLecunbpj3vRSHaAwDN1yQUKMUuivggisegP3per5McYjGMHCj\nFBJqOKC7Ai7YywonTHTcGcWhMRJqOCChhouuzxYD82IvK1ww0XFXFIemSKjhgoQapTBbDOwPsQwX\nTHTcFcWhKRJqlEJCDQfMFsMJnUJwQY5RH8VhSwh2OCChhgtmi+GCTiG4IMfoBorDFhDscEFCDRdM\noMEFnUJwQY7RDRSHLSDY4YKEGqXQFgcXdAoB+0Msd0NkZu1nKGYwGOT6+nrtx7h95XAS7G21fMBf\n20lIre+Et1ptcUDTasUy4zJKIMfwERHDzBzM895erhwyQw0H7DOBC9ri4IJOIbggx+iv3hWHBDtc\nkFCjFNri4KDGCgSxDBfkGP3Vu+KQYEcpJNRwwAFacFBzIphYRgnkGGhL74pDgh0lkFDDBW1xcFBz\nIphYRtPIMdCm3hWHBDtKIKGGCybQ4IA4hhNyDLSpd8WhRLCjeSQiKIUDtOCizVgmjuGEHANt4ioL\nWOL4ZTjgige4IJbhhBwDq4arLNBrnEgLFxygBRfEMlyQY8DdgdoPADRtVhICLGu4uaUTJzc03Nxq\n7TsnrURrIVqJ0Ki245lYhgtyDLhj5RDFcfwyVh3H4sMJJx/CCTkG0CyKQxRFEgIHHIsPJ5x8CBfk\nGHDSlX2lFIcoiiQEDpgphhPiGS7IMeCiS4d2URyiKJIQlMKx+HDBdSVwUGPVgxwDLrp0aBdXWfQM\nxy/DQZdm2IBlEMtwUDOOyTHQtFq58ne/9M23T3Q0/f8hrrLATBy/DBddmmEDlkEswwH7suGCA+i4\nyqJXOH4ZLjgWH6VwxQNctBnLxDFc1MyVj5x3SMePHa4+2cHKYY/Qm49S2DMFB5x8CBdtxzJxDBfk\nyhSHVZFQwwHtynDByYdwUSOWiWOUQK7cPorDSkio4YI9U3DBjDFcEMtwQK5cB8VhJSTUKKXtWTaS\nEJTCjDFcEMvA4siV66A4rISEGiWwZwoumDGGC2IZTtqc6CBXroPicIxZPThgzxRcMGMMF8QyXHDQ\nUT9QHIpZPfhglg0uiGW4IJbhgoOO+oHiUMzqoRxWpOGCWIaDtuNYIpZRRo1YZqKjHygORbCjDFak\n4YJYhoNacSwRy2hWzTGZiQ5/B2p8aUSsRcTbI+I1217/5Yi4fpfPXRIRGxFxdUQ8oqnnmQT7xRec\n3+p/LOBt1oo0sIqIZTggjuGiZiwfOe+Qjh87TK5srNbK4dMkvUfSPScvRMRA0r12+kBEfKGkJ0j6\nIkn3lfSGiPj8zLytiQdiVg9NY0Ua/397dx4sR1WGcfj3mrAHJYIhJEES2QRSSmkUQZSICxARXICK\nAilRKlKipeJGZAvulGUVCi4oKipQKIsSZAmihDUpvBEhCRgMBjGExWAUEkyA5POPPlebydxl7p2e\nnul5n6quuX3mzOmv65uZe86cXori25VYVfjKh1YF/k62KlFEtHaD0gTgp8BXgFMi4nBJI4CbgA8A\nf4mIUXVeNwsgIr6W1ucCsyNifl/bmjJlSvT09BSwF9aJyjg+v4xtWrWVdTiR38vWbGW8l/0+tmbz\nd7J1AkkLI2LKYOqWMXN4LvA5YNtc2ceAORHxqKS+XjceWJBbX5HKXkDSTGBmWl0jaemwI26+HYBV\nZQfRTbT5VttsNnrcHqAXQWx8bvXKB+LZ/6xtoAnnrPNUMmcjRm0/dsQ2o8cjICL2O/dfKzesefKx\nsuNqgkrmq+KGlbMKv5fbmT9nTdaC97Fz1lnaNV+7DLZiSweHkg4HnoiIhZKmprJxwNHA1IFeXqds\nk2nPiPgB8IPhRVosST2DHb1be3DOOo9z1lmcr87jnHUe56zzOGedpQr5avXM4RuBIyRNA7YkO+dw\nCbAeWJZmDbeWtCwidqt57Qpg59z6BGBl8SGbmZmZmZlVX0uvVhoRsyJiQkRMJLu4zO8jYnREjI2I\nian8mToDQ4A5wHRJW0iaBOwO3NWy4M3MzMzMzCqsre9zKOkIYEpEnBkRSyT9ErgPeB44uVlXKi1B\nWx/2anU5Z53HOesszlfncc46j3PWeZyzztLx+Wr51UrNzMzMzMys/bT0sFIzMzMzMzNrTx4cmpmZ\nmZmZmQeHzSTpaElLJG2UNCVXvr2kmyWtkXR+zWvmSVoq6U9pGVOn3bdLWihpUXo8uBX70w2Kylmq\nN0vSslT3kKL3pRsMJV+5OnMkLe7juZdIukbSPan9E4rah25TVM7S81PTZ3CJpFuKiL8bFZmzVOd1\nkjZIOqrZsXerAr8bj5V0b1rulPTqovahmxSYL0n6dup73CvpNUXtQ7cZYn/xhly/4vuSRtRpt+36\nH219QZoOtBh4L3BBTfk64AxgclpqHRsRPf20uwp4V0SslDQZmAuMb0K8VlDOJO1NdkXefYBxwE2S\n9ujgiyi1iyHlS9J7gTX9tHsycF9EvEvSy4Clki6JiGebE3ZXKyRnkrYDvgscGhEP9/UjjQ1JUZ8z\nUufoHLL/Y9Y8ReVsOXBQRKyWdBjZxTb2a0rE3a2ofB1GdjX/3cny9D2cr2YZSs6OiYinJAm4guy+\n7pfV1Gm7/odnDpsoIu6PiKV1ytdGxO1kb6ChtHt3RPTe03EJsKWkLYYRqiVF5Qw4ErgsItZHxHJg\nGfD6YYRqDC1fkkYBpwBf7q9pYNv0BT4K+CfZVZFtmArM2QeAqyLi4dTeE00KuesVmDOAjwNXAs5X\nExWVs4i4MyJWp9UFZPeYtmEq8DN2JPCzyCwAtpO0U7Pi7mZDyVlEPJX+HAlsTtbX2KQabdb/8OCw\nPfwkHRp1Rnpz9Od9wN0Rsb4VgVmfBsrZeODvufUVeLa3LF8Cvgk800+d84G9gJXAIuATEbGxBbFZ\nfYPJ2R7AaGWHeS+UNKM1oVkfBsyZpPHAe4Dvtyoo69dgPmd5HwauLy4cG8Bg8uW+R5uRNJfsx7Cn\nyWYPa7Vd/8OHlTZI0k3A2DpPnRYRVw+hyWMj4hFJ25L9mno88LM+tr0P2eE47xjCdrpWSTmrN2D0\nfWMGoZn5krQvsFtEfErSxH6qHgL8CTgY2BX4raTbcr/6WT9KytlI4LXAW4GtgPmSFkTEA41sr1uV\nlLNzgc9HxIaBfwe1WiXlrLf+W8gGhwc2sp1uVlK+3PcYhgL6i0TEIZK2BC4h62P8tqZK2/U/PDhs\nUES8rcntPZIen5Z0Kdmhh5sMDiVNAH4FzIiIB5sZQ9WVlLMVwM659QlkvwrZAJqcr/2B10p6iOz7\nboykeRExtabeCcDXI7vx6zJJy4FXAnc1MZbKKilnK4BVEbEWWCvpVuDVgAeHg1BSzqYAl6WB4Q7A\nNEnPR8SvmxhLZZWUMyS9CrgQOCwinmxiDJVW4vei+x5D1Oz+Yq7ddZLmkB32Wzs4bLv+hw8rLZGk\nkZJ2SH9vBhxOdsJrbb3tgGuBWRFxR2ujtLzB5gyYA0yXtIWkSWQnh3ug0WIR8b2IGBcRE8l+8X6g\nXucHeJhsBgpJOwJ7An9tVZz2fw3k7GrgTekzuTXZRRfub12k1muwOYuISRExMdW7AvioB4blGGzO\nJL0cuAo43rPy5Wnge3EOMCNdtfQNwL8j4tEWhmqJpFG953tKGglMA/5cp2r79T8iwkuTFrJzKVYA\n64HHgbm55x4iO8l0TaqzN7ANsBC4l+xCM98CRqT6RwBfTH+fDqwlm3buXcaUvb9VWIrKWVo/DXgQ\nWEr2i2vp+9vpS6P5qnntRGBxbv0k4KT09zjgRrLj/RcDx5W9r1VZispZWv8scF/K2SfL3teqLEXm\nLFd+EXBU2ftalaXA78YLgdW5vkdP2ftahaXAfAn4Tup7LAKmlL2vVVkazRmwI/AH/t9fPA8YWSdn\nbdf/UArMzMzMzMzMupgPKzUzMzMzMzMPDs3MzMzMzMyDQzMzMzMzM8ODQzMzMzMzM8ODQzMzMzMz\nM8ODQzMzqwBJH5QUknYrO5a8XFwTc2WzJR1cwLYuSjfJNjMzGxIPDs3MzIpzLbA/kL8R9VlA0weH\nZmZmwzWy7ADMzMyqKiL+Afyj7DjMzMwGwzOHZmbWFSQdJ+keSeskrZL0c0k71dR5SNLFkqZLul/S\nWkk9kg6s094nUv11ku6SdEBavyhX5wWHlUqK9NRpqTwkzU7PzZM0r852XtBmKnurpD+mbT8o6SN9\n7PPWks6RtFzSs+nxNEn+/29mZpvwzKGZmVWepJnABcAvgFnAOOCrwH6SXhMRa3LV3wTsCZwBrAO+\nBPxG0sSI+Fdq70TgXOBHwOXArsClwHYDhLI/MB+4KMUDsKLBfdkLuA7oAaYDWwCzgVHAhly9kcBc\nYO+0D4uAN6T9einw6Ua2a2Zm1efBoZmZVZqkEWSDo3kRMT1X/mfgNuBDwLdzL3kxsG9ErE71HgP+\nAEwDLk2zbmcB10fEibn2HgOu7C+WiFggCeCRiFgwxF06HXgaeEdErE3bvhN4EFiZq/d+4EDgoIi4\nNZX9Lm3/LEnnRMQTQ4zBzMwqyIeVmJlZ1e0JjAEuyRdGxO3A34CDaurP7x0YJovS48vT44S0XF7z\nuquB55sR8AD2B67rHRgCRMTfgTtq6h1Ktn93ShrZuwA3ApuRzSKamZn9j2cOzcys6l6aHh+t89xj\nued7/TO/EhHr02zblqmo9zzFJ2rqbZC0anihDspOwON1yh8HJuXWxwC7AM/10c72TY7LzMw6nAeH\nZmZWdb2DvbF1nhtLdu5eI3oHmWPyhenw1R0abCtvHdkhrbVqB6+PAjvWqVdb9iSwHDimj+091Ehw\nZmZWfT6s1MzMqm4p2aza9HyhpAPIZtZuabC9FWk5uqb83QzuR9dnga3qlP8N2EPS5rkY3wxsW1Nv\nPjBN0ja5ejsDb6ypdwOwM7AmInrqLK2Y5TQzsw7imUMzM6uSQ9OFYfL+DZwJXCDpYuBiYDzwFeAv\nwE8a2UBEbJR0NvBDSReSnXv4CuDUtK2NAzRxH/BOSTcAq4GVEbESuAyYCfw43bpiEnBKajPvy2QD\n0xslfQPYHDibTQ81vQQ4gewiNN8E7kl1dwWOAN4dEc80su9mZlZtHhyamVmVnFenbElETJb0DPBZ\nsgvHrCG7HcTnam5jMSgRcaGkUcCngOOAxcCxwDVsOpir9TGyq6NeQ3YbirOB2RFxs6STgM8A7wPu\nTm2/4AqoEXG/pGnAN8huzfEIcA7ZhWqm5uo9J+kQskHrTLLB5lqyq5peSzaDaWZm9j+KiIFrmZmZ\nWb8kvQ64C5gRET8vOx4zM7NGeXBoZmbWIEmTgJPJ7pP4FLAX8AWy2bjJPlzTzMw6kQ8rNTMza9x/\ngMnADGA02bmDNwGnemBoZmadyjOHZmZmZmZm5ltZmJmZmZmZmQeHZmZmZmZmhgeHZmZmZmZmhgeH\nZmZmZmZmhgeHZmZmZmZmBvwXNhMxBFWzHBIAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fc22bc726d8>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(15,10))\n",
    "plt.plot(lonsWGS84,latsWGS84,'.',label='HRAP')\n",
    "plt.axis('equal')\n",
    "plt.xlim(-115,-114)\n",
    "plt.ylim(44,45)\n",
    "plt.xlabel('Longitude', fontsize=16)\n",
    "plt.ylabel('Latitude', fontsize=16)\n",
    "plt.title('HRAP Grid in WGS84', fontsize = 18)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Move to a projected grid"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# +to_meter=4762.5\n",
    "p3 = Proj('+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-105 +k=1 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +no_defs',\n",
    "          preserve_units=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "x,y = transform(p2,p3,lonsWGS84,latsWGS84) # reproject the HRAP grid to WGS84"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7fc2001aaa20>"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA6UAAAJlCAYAAAAivlRGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XuYZFV97//3d6YZJSTCOHhBUWAiRCXR6DTYOYkG1IN4\nSTA5muAhEQmIIeqJOZqoMUe8xdvRoP4OmpCBiGaiISiIxgRRQX85sYVuNSoYwjgyMoEIDCMaFYam\nv+ePvUrLpvo20zOrevX79Tz1TPeqVXuvb9Wumvr0XnvvyEwkSZIkSaphVe0BSJIkSZJWLkOpJEmS\nJKkaQ6kkSZIkqRpDqSRJkiSpGkOpJEmSJKkaQ6kkSZIkqRpDqSRpj4qI6yPiigX2fV5EZEQcs4C+\n740Ir2u2AuypbWjYDcM2HhGvKc/noTXHIalthlJJWoSIOKZ8QXvZHH0yIj42o+2K0t673RURN0bE\n30bEz86zzrURcUd53G/N0e/6GevYWdo2RsRDFlHjvSPi9yLi0xFxSxnrtyPiqoh4S0Q8fKHLWq4i\n4tER8YGI2Fye+1sj4ssR8RcR8Zi+fgeUL+3HVByutFsi4pkR8Zra45C0co3UHoAkrSB3AqeVn/cF\nNgCnAE+LiNHMvHaWx50ErAG+AZwK/PUc69gGvLL8/FPAMcDvlHU8KjNvnWuAEbEe+BjwCOAzwFnA\nTcBPAj9flvWyiHhoZv77XMvq8zPAntjb83zgd5d6oRHxDOBi4BbgfcBm4ADg4cCvA9cBXyzdDwDO\nLD9fsdRj0S55P/BBYGftgSyBPbKND/BM4GTgNQPuewPwZrrPL0naIwylkrT3TGVmf6D8y4i4Bngn\n8CLgxbM87lTgcuAjwDsi4qcz8+uz9L19xjreExE3l+WfAvzv2QYXEfsCfw/8NPDrmXnRgD73Bv6A\neUJmROwDrM7MOzJzj3yZzcy7gLv2wKLfBPwAOCozt/XfUeq67x5Y50AR8VOZ+d29tb5dMWxjzMy7\ngbtrj6NfRASwX2b+52Ietwe38cWMYQqYqjkGSe1z+q4k1fWp8u/hg+6MiMfS7aE8H9hE9wX1lEWu\n49Ly78Pm6Xca3d7A/z0okAKUkPmmzLyxb4y9Y86OjIg/i4htwB3AWLl/4PGAEXFaRPxrRNxZpsn+\nPhALLWrQ8Xa9tojYPyLeExE3l+m3/zciHrfARR8OXDszkJb678rMb5V1HUO39xrgzL5p09fPGNNv\nRsQ/RcR3I+L7EfH5iHjWgHqyjP9Jpf9/Ah/tu3//Mn16c3nObilTjNcPWNaBEXF2RNxQpnHfUH5f\nN6DvoRHxoYj4TkTcHhEfiYjDBr1uc40xIh4UEW+PiC9FxI7yvF8TES+PiNUzltM77vPJZfvZWmr6\nckScOPhlgYh4eET8fXkub4+ICyPigbMs+5gZ7Wsi4o/K+L5fHj8RES+abX19j72iPB/ry/Nze3m+\nLpr5/MePpvg/LyJeWP7wdAfwsr4+T4iIy8pyfhARX4iIUwesd+AxpRFxUNm+v1le3xsj4pyIuP+A\nvveJiD+NiK+V12R7ee1O7NVGt5e09/r2bs8rbQOPKS3bzfsj4lvltft6RLwxIn5iRr/e43+m3L+t\n9P+XiHjafM+9pJXBPaWStGt+IiIOXILl/HT597ZZ7j8V+B7wocz8XkT8PXByRLw6M6cXuI5e4J1z\n6i7QC0obF7jcmTbR7WF8O92e1Jtm6xgRL6GbGvwvwB8DPwH8IXDzLq57pkvppt++DlgH/E/g4xFx\n6AL26n0dODIi/ktm/vMc/b5Gt9f4LOAi4MOl/Yd7wyLiDcCrgH8E/hcwDfwa8HcR8aLMPHvGMkeB\n/wb8Jd0fInrL2R/4Z+ChwHnA1cBBwO8Bn49u+vfWGX0fVvp+AXgMcAbwxIg4uvcclJD6/wMPAP68\n1PR4uj3z+81S98AxAo+im958Ed1zuA/wVLqpn+uBFwxY1lvKet5Dt82cAnwgIu6dme+d0ffBdFOk\nL6LbVh5dlnkf4LhZxkqpcw3dNnEM8Am6KfB3AD9Xxvx/5np8sR/d83Il3RT5w+me/7GIeExm/seM\n/i+h2/b+EvgP4IYyll8pNfwH3Xvlu8CJwMaIWJ+Zr5qnlocCn6Ob0n8u3XP9MLrX99iyLdxe+h4A\n/BNwJHAh3fO8mm57eAbdNOc/pdtJ8Xjgt/tWNeu2HxGHlOdh/7LMf6N7bl8J/GJEPKnsYe13Pt0f\n1d5Wxv4S4OKIOCIzr5+rZkkrQGZ68+bNm7cF3ui+eOUCbh+b8bgr6MLKgeX2ELrjuK4v/Z82YF33\npgur7+1rO6H0f+qA/tfThYreOg6j+5L/bbovgz87T23b6ab/zmxf3bfM3m3fvvtfU8Z0BTAyy7iu\n6Pv9ALqgfQ3wE33tB5fnKIFjFvBavLf7b+yebcC7Z7Q/u7S/YAHLfRZdeEzgy3Rh7XeAQwf0PbT0\ne82A+x5b7nvjgPsuBr4D/FRfW2/befKA/u+kC/yPntF+SFlO/zbyp2U5vzej7wtL++v72t5a2k6a\n0bfXfsWM9rnGuC8QA9rfTzed9qC+tueV5WwF9u9r37+03TZjG7u+9P+NGcs+u7Q/fMCyj+lr+6M5\nXotVC9gmriiPf8eM9l8r7X/e13ZMabsNuP+A99JWuvfkg/ra1wD/tzxPh8+zjX+E7o83B89oH6Wb\nZvuavrZ3l7GcPlfdg9bTd99ryjIO7WvbxIDPLbrDAxI4dcDjP9a/fQBHlfY3zff8e/Pmrf2b03cl\nadecA/zXWW6z2Y9u790twDfp9pasAU7OzI8P6P/rwFp+fG/U39N9If2dWdbx8L51bKHbU3YrcEJm\nfnWemu5DF3BmekTfMnu3Fw7o9468596RQY6j2zN6dmZ+v9eY3XTZTQt4/EKcNeP3T5d/B06T7peZ\nFwJPoNuz9BC6vXHnAt8oUzfvt8AxnET3pfv86KbT/vAGXEJ3IqpfmPGYf8nMT/Y3RESUZX0W+PcZ\ny/keMM6P7yn8NbrX6JwZy/4Lum3h1/rafoVuj/YHZvR92xx13WOMAJn5g8zMMuY1EXHfMsZL6fbE\njQ5Y1nuy7NUry7id7o8Aa+nCXb8bM/OCGW2913W+qeknATvo9pzPHPdCZxxAt9e3/7EXAdfS/YFp\npvdl5sw9/xsoe7uzbwp8Zu6kC3Sr6P7wNFDZC/4Muu3njhnbwvV0J+U6rvRdRbcH9mt0e2t/zCLr\n7h/DKuBXgS8O+Nx6Ez+aDTDTO3vbR1n/VXR7ied9T0pqn9N3JWnXXDfoizlAlyEGuoMuBEB3spzn\n0oXY2f5AeCpduNgWEf1fui8Dnh0RB+Y9z6Z7Pd0ZO6E7++iNmbl5jjr6fYcumM70DX4Uth/N7IHl\n3xa4nt4xeP864L5rFriM+Wzp/yUzt5fX5R7HVA6Smf8E/FMJhIcDx9JN1fxVuqmfT1nAYh5Bd4zs\noDp7HjDj90HP4f3oxn0c3fYwSH/AOAyYmPkHgsyciohr6fbg9ve9cmZAycybI+Lbs6xr4OscESPA\nK+i264dxz+OD1w542NcGtPW2gZnHym6Z2ZFu7z7M/7oeDnwpM++Yp99cvp33nKILXQ3PjIj9MvN7\nfe2DnqfDyr9XD7iv90ejexwj3Odn6D4vTi23QXrP04F0z/k/9ofBJXA/urNx36OGzLwtIm5icA2D\nXr/bWOB7UlLbDKWStPfc3R9kI+JCuilt50TEFzLzy333HUYXhILZw95vAe+Y0fa92cLyAnwVeEJE\nHJaZvRP4UL5of7KMa649od+f475+vbAy6Ivygk90NJfszsA617oXupyke/7/LSLOp/siflxEHJwD\nToQ0YF1Jd1zlbOOZ+cV+0HPYG/Mn6Y7BrG221/nP6M4g/bd0U4hvpps2/li6cQ/648titoG5zqi7\nkNd1d4PZbI+fbd1zvZa7qvf4v+bHZ1D0+8GMvkt9OaZdrWFJ3pOS2mQolaRKMnM6ujPOXkO397F/\nCuYpdF/Wnk93/NlMb6DbUzIzlO6O3rTV0+hOzrOn9C5n8wh+NP2SvrahlJl3RMSX6PYCPZjumrBz\nfeG/Djge+GZmDtojuFC30G0D91ngHxy2AD8TESP9e0vLnswj+PE9VtcDD4uIVf17S8tZXA9Y5Dh/\nG/hsZv7Y2XNn7OWf6ZF0U1H79baBQXvWdtW/AY+IiHvlrl+iaG1EPHDA3tKHAzfP2Es6m962f+SA\n+x5Z/p2r7s1029yaBWwLt9BNWf75BYxrMcH1Zrppt/eoISLW0p2A60uLWJ4keUypJNWUmdcBfwP8\n14j4JfjhMVvPA76SmRsz88KZN7pjAH82Io5awuFspJtq+ocRMeiYMFiavRqX0e3NeWH/5SMi4mDg\nvy/B8ndLRBwfA+Zgl2NJf5HuZDLXlebemXYHXbv0/eXfN8aMS6KU5d3j8h2DlLC4CTg6BlxKZsCy\nLqabYnnajG7PL+39l/v5KF2IeM6Mvi9j8e5mxvYREfvRnaF4NmeU4yR7/fcHfpcuhH9mF8Ywm010\nU1n/ZOYdg17rObxixmN/jW5K7cULfPwX6I4nPyX6LmUT3fVv/5AuHH5ktgdn5nbg48CvR8TYzPuj\nc7/Sd5ruc+KRMfhyM/11/2dpm/cavGW5HwUeExHHz7j7FXTfLQdeUkqSZuOeUkmq7410U3FfCzyJ\nbo/pQ+hOrjObD9Gd1fJU4KqlGERm/iAink43pfjD0V2/8BN0l664D90eod+kCx837MZ6dkTE/6Lb\nO/zPEfE+uhMf/S5d2HvM7tSxBC4Ebo6Ij9HtxZ6i2zv623THgL4uM2+DHx6ruhk4MSK+DnyLbgr1\nRzPzqog4k+51/VJE/B1wI10I3AA8je5EVwvxKrpAfEFEXEB3cqOddGfffRowSfeHDOjOnPts4Ozo\nrnP7Rbrn9FS6k/K8tW+5b6H7Q8BfRcTRdH+U+KWyrltZ3B60C4EXRMTf0k01fgDdCbm2z/GYW+ku\naXMeXaA9he5EQKf1nwRrCbyT7njuPyl/yPkE3THeR9KFyicvYBm30oXBB9Gdjbd3SZhv0b0X55WZ\nd0d3XdSLgKsi4hy6vY6/SXdd3zeWP1TN5Qy6y7x8trx3vkgXBNfTnSTpfX3j+RPgiXSXmzmuPC7o\ntocRfnQJmHHgRcC7y2Wn7gI+3z+Nf4Y/pjvO/OKIeDfdHtwnlDo+y+xTiyVpIEOpJFWWmdeWoHFi\nRPwyPzqByYfneMxXI+LfymP+IDN/MFvfRY5lS0RsoAsTzwJeSneZju/RffHcCJybmdfu5nreHhH/\nSXf90DfRhdy3AbfTnTG4plPojgN9Et2X9p+kOyHLF4CXZOaHZvQ/ie5sv2+kC9db6fYkkZmvi4hJ\n4H/QXZdxP7rpj18Ffn+hA8rM2yPiF+lej9+gCx9TdFOI/4m+a8v29X0t3YmZTqELTn8OnJl912nN\nzFvLHvq3073mSXctzmPp/tixmO3qf9IFrN74bqA7A/BVlGOSB3g53fUxX0QXYq+juzzN3yxivfPK\nzJ0llL2ULoS/kS6UXgf81QIX8z26gHcW3Vl4g+76sy/NzFmvyTtgLB+NiCfRBcY/pPvDxNeA52fm\nvNcIzswbynv05XTP82+VWm6g2+4u6Ou7IyJ+gS5E/jrdWXG/S/fHlv+vb7EfoAuqJ9L9QWMV3XYz\nMJRm5taIeBzd2Yx/i26q9za69/IbFngWbkn6oVjaE7JJkqTlLiLW0e0Z/IvM/N09sPzn0YXBYzPz\niqVe/lIrswYOzcxD9/J63w88JzPdiSCpaR5TKknSChYR+w5ofnn597K9ORbdw4Po9qxLUtOGKpRG\nxIsj4tqIuDoi3trX/sqI2Fzue0pf+/UR8ZWI+FJETPS13zciLouI68q/a0t7RMS7yrK+XI616T3m\n5NL/uog4ua99Q1nH5vJYT10uSWrJP0TE+eX/4JdExEfpppX+Mws/gY+WUET8l4h4Hd1xmpfWHo8k\n7WlDMx0kIo6lOzbiUZl5Z+9MghHxSLpjHI6k+4vhJyPiiL5r0B074OLxrwA+lZlvjohXlN9fTneM\n0OHl9jjgPcDjytnmzgRG6Y6nmYyISzJzR+lzOt1JAD5Od3r/f9gjT4IkSXvfR4HnAs8E9qU7NvDt\nwGvnuN6r9qzTgWcAf0d3HKwkNW1ojiktJ/k4Z+Z1tyLilQCZ+aby+6XAazLzcxFxPTA6M5RGxLXA\nMZl5U0QcBFyRmT8TEX9Rfv5Af7/eLTNfUNr/gu7MelcAl2fmw0v7c/r7SZIkSZJ2z9DsKaW7oPfj\nI+JP6c4i97LMvIruAuXjff22lTbo9mp+IiKS7mQM55T2B/TOhFeCae/6bQ/mxy9j0FvWXO3bZln3\nPUTE6XR/3WS//fbb8PCHP3whdUuSJElScyYnJ2/NzPvN12+vhtKI+CTwwAF3vaqMZS3ddbqOorsW\n23oGX6i9t3v3FzPzxhI6L4uIf83Mz841hFmWtdj2gUooPgdgdHQ0JyYmZusqSZIkSU2LiK0L6bdX\nQ2lmznpx6og4A/hwdvOJr4yIaeBAur2TD+nrejDdxcfJzN6/N0fERcDRdBdt/lZEHNQ3fbd35rrZ\nlrWNbgpvf/sVpf3gQeuWJEmSJO2+YTr77sV0F6UmIo6gu5j0rcAldBeHv1dEHEZ3kqIrI2K/iPip\n0n8/4Di6i5FTHtM7g+7JwEf62p9bzsI7BtxepvleChwXEWvLmXqPAy4t9303IsbKWXef27csSZIk\nSdJuGqZjSs8DzouIrwI7gZPLXtOry0mQrgGmgBdm5t0R8QDgonKFlhHgbzLzH8uy3kw3/fdU4JvA\ns0v7x4GnAZuB7wOnAGTmbRHxeuCq0u91mXlb+fkM4L10ZyT8BzzzriRJkiQtmaE5+25rPKZUkiRJ\n0koWEZOZOTpfv2GavitJkiRJWmEMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJ\nkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmS\npGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKk\nagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRq\nDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoM\npZIkSZKkagylQ2py6w7Ovnwzk1t31B7KkmipnpZqAetZDlqsSdLK5Wea9ia3t+VhpPYAdE+TW3dw\n0sZxdk5Ns2ZkFZtOG2PDIWtrD2uXtVRPS7WA9SwHrdU0uXUH41u2M7Z+3bKuo6elelqqBdqrpxWt\nfaZBm9taKzW5vS0fhtIhNL5lOzunpplOuGtqmvEt25f1RtdSPS3VAtazHLRUU2tfDlqqp6VaoL16\noJ0voi19pkG721orNbm9LR9O3x1CY+vXsWZkFasD9hlZxdj6dbWHtFtaqqelWsB6loOWahr05WA5\na6melmqB9urpfRF9+yeu5aSN48t6GmJLn2nQ3rYGbdXk9rZ8uKd0CG04ZC2bThtr4i+i0FY9LdUC\n1rMctFRT78vBXVPTTXw5aKmelmqB9uppaW9PS59p0N62Bm3V5Pa2fERm1h5Dk0ZHR3NiYqL2MCRp\nqLQyBbGnpXpaqgXaqqe3p7T3RbSlKXstaGlb62mxplYst9cmIiYzc3TefobSPcNQKkmSlspy+yIq\nSbDwUOr0XUmSpCG34ZC1hlFJzfJER5IkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJ\nkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmS\npGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKk\nagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRq\nDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoM\npZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagyl\nkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJInJrTs4+/LNTG7dUXsoS6K1eqS9xfeO\nahipPQANNrl1B+NbtjO2fh0bDllbezi7zXqGV0u1gPUMO+sZTpNbd3DSxnF2Tk2zZmQVm04bs54h\n0sp21q+1mlqpp7X3DrTz2kBbtcxkKB1CrX0gWM/waqkWsJ5hZz3Da3zLdnZOTTOdcNfUNONbti/b\nWqCtelraznpaq6mlelp670Bbr01LtQzi9N0hNOgDYTmznuHVUi1gPcPOeobX2Pp1rBlZxeqAfUZW\nMbZ+Xe0h7ZaW6mlpO+tpraaW6mnpvQNtvTYt1TKIe0qHUO8D4a6p6SY+EKxneLVUC1jPsLOe4bXh\nkLVsOm2smWlhLdXT0nbW01pNLdXT0nsH2nptWqplkMjM2mNo0ujoaE5MTOzy41ubM249w6ulWsB6\nhp31SIvX4nbWWk2t1dOSll6b5VhLRExm5ui8/Qyle8buhlJJkiRJWs4WGko9plSSJEmSVI2hVJIk\nSZJUjaFUkiRJklSNoVSSJEmSVI2hVJIkSZJUjaFUkiRJklTNUIXSiHhxRFwbEVdHxFv72l8ZEZvL\nfU+Z8ZjVEfHFiPhYX9thEfH5iLguIv42ItaU9nuV3zeX+w+dbx0RcXxp2xwRr9iT9UuSJEnSSjM0\noTQijgVOAB6VmUcCbyvtjwROBI4EjgfeHRGr+x76+8DXZizuLcBZmXk4sAM4tbSfCuzIzIcBZ5V+\ns66jrOds4KnAI4HnlL6SJEmSpCUwNKEUOAN4c2beCZCZN5f2E4APZuadmfkNYDNwNEBEHAw8HdjY\nW0hEBPBE4MLSdD7wzL5lnV9+vhB4Uuk/2zqOBjZn5pbM3Al8sPSVJEmSJC2BYQqlRwCPL9NqPxMR\nR5X2BwM39PXbVtoA3gH8ETDdd/864NuZOTWg/w+XVe6/vfSfbR1zrfseIuL0iJiIiIlbbrll/ool\nSZIkaYUb2Zsri4hPAg8ccNeryljWAmPAUcAFEbEeiAH9MyKeAdycmZMRcUz/agb1n+e+2doHhfYc\n0NbdkXkOcA7A6OjorP0kSZIkSZ29Gkoz88mz3RcRZwAfzswEroyIaeBAur2TD+nrejBwI/CrwK9G\nxNOAewP3iYi/Bn4bOCAiRsre0F5/+pa1LSJGgP2B2+ZYB3O0S5IkSZJ20zBN372Y7lhQIuIIYA1w\nK3AJcGI5c+5hwOHAlZn5ysw8ODMPpTtJ0acz87dKqL0ceFZZ7snAR8rPl5TfKfd/uvQfuA7gKuDw\ncjbfNWU9l+y5p0CSJEmSVpa9uqd0HucB50XEV4GdwMklMF4dERcA1wBTwAsz8+55lvVy4IMR8Qbg\ni8C5pf1c4P0RsZluD+mJAJk56zoi4kXApcBq4LzMvHrJKpYkSZKkFS663KelNjo6mhMTE7WHIUmS\nJElVRMRkZo7O12+Ypu9KkiRJklYYQ6kkSZIkqRpDqSRJkiSpGkOpJEmSJKkaQ6kkSZIkqRpDqSRJ\nkiSpGkOpJEmSJKkaQ6kkSZIkqRpDqSRJkiSpGkOpJEmSJKkaQ6kkSZIkqRpDqSRJkiSpGkOpJEmS\nJKkaQ6kkSZIkqRpDqSRJkiSpGkOpJEmSJKkaQ6kkSZIkqRpDqSRJkiSpGkOpJEmSJKkaQ6kkSZIk\nqRpDqSRJkiSpGkOpJEmSJKkaQ6kkSZIkqRpDqSRJkiSpGkOpJEmSJKkaQ6kkSZIkqRpDqSRJkiSp\nGkOpJEmSJKkaQ6kkSZIkqRpDqSRJkiSpGkPpkJrcuoOzL9/M5NYdtYeyJKxneLVUS09LNbVUC7RV\nT0u1SHtbS++flmppVUuvUUu19BupPQDd0+TWHZy0cZydU9OsGVnFptPG2HDI2trD2mXWM7xaqqWn\npZpaqgXaqqelWnomt+5gfMt2xtavW/a19LRUU2u1tPL+aamWnpa2NWjrNWqplpncUzqExrdsZ+fU\nNNMJd01NM75le+0h7RbrGV4t1dLTUk0t1QJt1dNSLfCjLzpv/8S1nLRxvIm/wLdUU0u1QFvvn5Zq\ngfa2NWjrNWqplpkMpUNobP061oysYnXAPiOrGFu/rvaQdov1DK+WaulpqaaWaoG26mmpFmjzi05L\nNbVUC7T1/mmpFmhvW4O2XqOWapkpMrP2GJo0OjqaExMTu/z4FqdOWM9waqmWnpZqaqkWaKue1mo5\naeM4d01Ns08jU8JaqqmlWnpae/+0VEtr2xq09xotp1oiYjIzR+ftZyjdM3Y3lEqStDctty86C9FS\nTS3VouHmtqalZCitzFAqSZIkaSVbaCj1mFJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWG\nUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZS\nSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJ\nkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmS\nJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIk\nSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJ\nUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlLtFZNbd3D25ZuZ3Lqj\n9lCWREt2GlCvAAAgAElEQVT1tFRLT4s1SXtDS++dlmqB9urRcGtpe2uplpaN1B6AZje5dQfjW7Yz\ntn4dGw5ZW3s4u2xy6w5O2jjOzqlp1oysYtNpY9YzJFqqpae1mlr5HOhpqZ6WaoG23jst1QJt1tPa\ne6e1elrZ3lqqpae17a3HUDqkWnoTjW/Zzs6paaYT7pqaZnzL9mVbC7RVT0u19LRUU0ufA9BWPS3V\n0tPSe6elWqCtelp777RWD7S1vbVUC7S5vfU4fXdIDXoTLVdj69exZmQVqwP2GVnF2Pp1tYe0W1qq\np6VaelqqqaXPAWirnpZq6WnpvdNSLdBWPa29d1qrB9ra3lqqBdrc3nrcUzqkem+iu6aml/2baMMh\na9l02lgzUw1aqqelWnpaqqmlzwFoq56Waulp6b3TUi3QVj2tvXdaqwfa2t5aqgXa3N56IjNrj6FJ\no6OjOTExsVvLaHXOuKSFa+1zoKV6WqpF2ptae++0Vo+G23Lb3iJiMjNH5+1nKN0zliKUSpIkSdJy\ntdBQ6jGlkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJ\nkqRqDKWSJEmSpGoMpZIkSZKkagylkiRJkqRqDKWSJEmSpGqGKpRGxIsj4tqIuDoi3trX/sqI2Fzu\ne0ppu3dEXBkR/1L6v7av/2ER8fmIuC4i/jYi1pT2e5XfN5f7D51rHaX9+NK2OSJesTeeB0mSJEla\nKYYmlEbEscAJwKMy80jgbaX9kcCJwJHA8cC7I2I1cCfwxMx8NPDzwPERMVYW9xbgrMw8HNgBnFra\nTwV2ZObDgLNKv1nXUdZzNvBU4JHAc0pfSZIkSdISGJpQCpwBvDkz7wTIzJtL+wnABzPzzsz8BrAZ\nODo7/1n67FNuGREBPBG4sNx3PvDMvmWdX36+EHhS6T9wHeW2OTO3ZOZO4IOlryRJkiRpCQxTKD0C\neHyZVvuZiDiqtD8YuKGv37bSRtmb+SXgZuCyzPw8sA74dmZOzezfv6xy/+2l/2zrmHXdg0TE6REx\nERETt9xyy6KKlyRJkqSVaGRvriwiPgk8cMBdrypjWQuMAUcBF0TEeiAG9E+AzLwb+PmIOAC4KCJ+\nFvjWbP3nWNZs7YNCew5oo4znHOAcgNHR0Vn7SZIkSZI6ezWUZuaTZ7svIs4APpyZCVwZEdPAgXR7\nJx/S1/Vg4MYZy/12RFxBdzzo24EDImKk7A3t799b1raIGAH2B26bZx1zrluSJEmStOuGafruxXTH\nghIRRwBrgFuBS4ATy5lzDwMOpwut9yt7SImIfYEnA/9aQu3lwLPKck8GPlJ+vqT8Trn/06X/wHUA\nVwGHl7P5rqE7GdIle+wZkCRJkqQVZq/uKZ3HecB5EfFVYCdwcgmMV0fEBcA1wBTwwsy8OyIOAs4v\nZ8hdBVyQmR8ry3o58MGIeAPwReDc0n4u8P6I2Ey3h/REgMwcuA6AiHgRcCmwGjgvM6/es0+DJEmS\nJK0c0eU+LbXR0dGcmJioPQxJkiRJqiIiJjNzdL5+wzR9V5IkSZK0whhKJUmSJEnVGEolSZIkSdUY\nSiVJkiRJ1RhKJUmSJEnVGEolSZIkSdUYSiVJkiRJ1RhKJUmSJEnVGEolSZIkSdUYSiVJkiRJ1RhK\nJUmSJEnVGEolSZIkSdUYSiVJkiRJ1RhKJUmSJEnVGEolSZIkSdUYSiVJkiRJ1RhKJUmSJEnVGEol\nSZIkSdWMLKZzRJw3x93TwO3AJPDhzLxjdwYmSZIkSWrfokIpcCywP3AAMAXcChxYlvPt0ucPgK9H\nxLGZuW2pBipJkiRJas9ip+/+d7q9of8NuHdmPgi4N/Bs4DvAM4DHlbY3LeE4JUmSJEkNWuye0rOA\nt2TmRb2GzJwGPhQR9wfekZlHR8SbgFcv4TglSZIkSQ1a7J7SRwNfn+W+rwM/W36+Bli7q4OSJEmS\nJK0Miw2l/wE8a5b7ng18q/x8H2DHrg5KkiRJkrQyLHb67juAsyLiQcCFwM3A/ekC6dOAl5R+jwe+\nuFSDlIbN5NYdjG/Zztj6dWw4xEkBkiRJw8jvbMvDokJpZr4zIr5Hd7zo0/vu2gY8PzPPLb+fDfxg\naYa4crX0JmqtlpM2jrNzapo1I6vYdNrYsq6ppdemp6WaWqoF2qqnpVp6WqqppVpa1Nrr01I9rdXS\n0nc2aOv16bfYPaVk5saIOBc4GDgIuAnYlpnZ1+f6JRvhCtXSm6ilWgDGt2xn59Q00wl3TU0zvmX7\nsq2ntdcG2qqppVqgrXpaqqWnpZpaqqWnpS+irb0+LdXTUi3Q1nc2aO/16bfYY0oByM4NmXll+Tfn\nf5QWY9CbaLlqqRaAsfXrWDOyitUB+4ysYmz9utpD2mWtvTbQVk0t1QJt1dNSLT0t1dRSLfCjL6Jv\n/8S1nLRxnMmty/u0Ha29Pi3V01It0NZ3Nmjv9em36D2lEXEfuuNHH0p3PdJ+mZmvX4qBrXS9N9Fd\nU9PL/k3UUi0AGw5Zy6bTxpr4i3Vrrw20VVNLtUBb9bRUS09LNbVUC7S3t6e116elelqqBdr6zgbt\nvT79YjE7OSPiF4GPAgfM0iUzc/VSDGy5Gx0dzYmJid1aRmtTdVqppTUtvjYt1dRSLdBWPS3V0tNS\nTa3VctLG8R9+EW1hyl5Lrw+0VU9LtbRoub0+ETGZmaPz9ltkKL0KWA08H/hKZu7c9SG2bSlCqSRJ\nEiy/L6KSBAsPpYudvvsI4Dcyc3LXhiVJkqTF2nDIWsOopGYt9kRH3wTutScGIkmSJElaeRYbSl8L\nvKKc7EiSJEmSpN2y2Om7zwAeAHwjIj4H3Dbj/szMk5dkZJIkSZKk5i02lP4SkMB3gCMH3O/1SiVJ\nkiRJC7aoUJqZh+2pgUiSJEmSVp7FHlMqSZIkSdKSmXdPaUQ8FLgpM+8qP88pM7+5JCOTJEmSJDVv\nIdN3vwH8AnAlcD3zHze6ejfHJEmSJElaIRYSSn8H+Hrfz57MSJIkSZK0JOYNpZl5ft/P792jo5Ek\nSZIkrSie6EiSJEmSVM1ir1NKRPwy8BzgocC9Z9ydmfmkpRiYJEmSJKl9iwqlEfEC4D3AduA64M6Z\nXZZoXJIkSZKkFWCxe0pfCvwN8DuZuXMPjEeSJEmStIIs9pjSBwN/ZSCVJEmSJC2FxYbSSWD9nhiI\nJEmSJGnlWWwo/R/ASyLiCXtiMJIkSZKklWXeY0oj4gYg+5r2By6PiO8DO2Z0z8w8ZAnHJ0mSJElq\n2EJOdPQpfjyUSpIkSZK0JOYNpZn5vL0wDkmSJEnSCrSoY0oj4tUR8aBZ7jsoIl69NMOSJEmSJK0E\niz3R0ZnAwbPc96ByvyRJkiRJC7LYUBpz3LcWuHM3xiJJkiRJWmEWcvbdY4An9jW9ICKeMaPbvsDT\ngauXbmiSJEmSpNYt5Oy7vwz8Sfk5gVMG9NkJXEN3HVNJkiRJkhZk3um7mfnazFyVmavopu+O9X7v\nu907Mx+bmZ/b80OWJEmSJLViwceURsQa4F14zVJJkiRJ0hJZcCjNzJ3A8+mOH5UkSZIkabct9uy7\nXwJ+bk8MRJIkSZK08iw2lL4UeFlEPCMi5ro8jCRJkiRJ81rI2Xf7/R2wP/ARYCoibubHjzHNzDxk\nqQYnSZIkSWrbYkPpp/BER5IkSZKkJbKoUJqZz9tD45AkSZIkrUCLPaZUkiRJkqQls+hQGhE/FxEX\nRsQtETEVETdHxAUR4Vl5JUmSJEmLsqjpuxFxFPAZ4AfAJcB/AA8EfgV4ekQ8ITMnl3yUkiRJkqQm\nLfZER28Cvgo8KTO/22uMiJ8CPlnuP27phidJkiRJatlip++OAW/qD6QA5fe3AL+wVAOTJEmSJLVv\nsaF0vsvBeLkYSZIkSdKCLTaUfh744zJd94ciYj/g5cD4Ug1MkiRJktS+xR5T+sfAFcDWiPgYcBPd\niY6eDuwLHLOUg5MkSZIktW1RoTQzr4yIMeDVwFOA+wK3AZ8GXp+ZX1n6IUqSJEmSWrXYPaVk5peB\nZ+2BsUiSJEmSVph5Q2lEvHoxC8zM1+36cCRJkiRJK8lC9pS+ZgF9+s+6ayiVJGmZmdy6g/Et2xlb\nv44Nh6ytPZzd1lo90t7ie0c1LCSU7jPP/T8PvIHuGNPrdntEAtr7QLCe4dVSLT0t1WQtw62Vmia3\n7uCkjePsnJpmzcgqNp02Zj1DppVtraelelqrxffOcGutnp55Q2lm3j2oPSIOp9sr+mzg34HTgb9a\n0tGtUK19IFjP8Gqplp6WarKW4dZSTeNbtrNzaprphLumphnfsn3Z1gLt1dPStgZt1dNSLeB7Z9i1\nVk+/xV6nlIh4SERsBK4Gngi8DDg8MzfOFmC1OIM+EJYz6xleLdXS01JN1jLcWqppbP061oysYnXA\nPiOrGFu/rvaQdktr9bS0rUFb9bRUC/jeGXat1dNvwWffjYj7A6+i2yN6B91e0rMy83t7aGwrVu8D\n4a6p6SY+EKxneLVUS09LNVnLcGuppg2HrGXTaWPNTAlrrZ6WtjVoq56WagHfO8OutXr6RWbO3SFi\nf+DlwIuBAP4P8JbM3LHnh7d8jY6O5sTExC4/vrX54tYzvFqqpaelmqxluLVYk4ZTa9taS/W0VEuL\nWnt9lls9ETGZmaPz9ltAKL0N2B/4BN0JjW6aq39mblnEOJu1u6FUkiRJkpazhYbShUzfPaD8+xTg\nuAX0X72APpIkSZIkLSiUnrLHRyFJkiRJWpEWckmY8/fGQCRJkiRJK8+iLwkjSZIkSdJSMZRKkiRJ\nkqoxlEqSJEmSqjGUSpIkSZKqMZRKkiRJkqoxlEqSJEmSqjGUSpIkSZKqMZRKkiRJkqoxlEqSJEmS\nqjGUSpIkSZKqMZRKkiRJkqoxlEqSJEmSqjGUSpIkSZKqGapQGhEvjohrI+LqiHhrX/srI2Jzue8p\npe0hEXF5RHyt9P/9vv73jYjLIuK68u/a0h4R8a6yrC9HxGP7HnNy6X9dRJzc174hIr5SHvOuiIi9\n82xIkiRJUvuGJpRGxLHACcCjMvNI4G2l/ZHAicCRwPHAuyNiNTAFvDQzHwGMAS8sfQFeAXwqMw8H\nPlV+B3gqcHi5nQ68p6zjvsCZwOOAo4Eze0G29Dm973HH75EnQJIkSZJWoKEJpcAZwJsz806AzLy5\ntJ8AfDAz78zMbwCbgaMz86bM/ELp+13ga8CD+x5zfvn5fOCZfe3vy844cEBEHAQ8BbgsM2/LzB3A\nZcDx5b77ZObnMjOB9/UtS5IkSZK0m4YplB4BPD4iPh8Rn4mIo0r7g4Eb+vpt40fhE4CIOBR4DPD5\n0vSAzLwJoPx7/3mWNVf7trnWPWMcp0fERERM3HLLLXMWK0mSJEmCkb25soj4JPDAAXe9qoxlLd1U\n3KOACyJiPTDoGM7sW+ZPAh8CXpKZ35lvCLMsa7HtA2XmOcA5AKOjo7P2kyRJkiR19moozcwnz3Zf\nRJwBfLhMk70yIqaBA+n2Tj6kr+vBwI3lMfvQBdJNmfnhvj7fioiDMvOmMgW3NxV4tmVtA46Z0X5F\naT940LolSZIkSbtvmKbvXgw8ESAijgDWALcClwAnRsS9IuIwupMNXVnOgnsu8LXM/LMZy7oE6J1B\n92TgI33tzy1n4R0Dbi/Tey8FjouIteUER8cBl5b7vhsRY2V9z+1bliRJkiRpN+3VPaXzOA84LyK+\nCuwETi57Ta+OiAuAa+jOuPvCzLw7In4J+G3gKxHxpbKMP87MjwNvppv+eyrwTeDZ5f6PA0+jO1nS\n94FTADLztoh4PXBV6fe6zLyt/HwG8F5gX+Afyk2SJEmStASiy31aaqOjozkxMVF7GJIkSZJURURM\nZubofP2GafquJEmSJGmFMZRKkiRJkqoxlEqSJEmSqjGUSpIkSZKqMZRKkiRJkqoxlEqSJEmSqjGU\nSpIkSZKqMZRKkiRJkqoxlEqSJEmSqjGUSpIkSZKqMZRKkiRJkqoxlEqSJEmSqjGUSpIkSZKqMZRK\nkiRJkqoxlEqSJEmSqjGUSpIkSZKqMZRKkiRJkqoxlEqSJEmSqjGUSpIkSZKqMZRKkiRJkqoxlA6p\nya07OPvyzUxu3VF7KEvCeoZXS7X0tFSTtQy3FmuS9rTW3jet1aPh1ur2NlJ7ALqnya07OGnjODun\nplkzsopNp42x4ZC1tYe1y6xneLVUS09LNVnLcGutpsmtOxjfsp2x9euWdR091jOcWnzftFQPtLOt\n9bRUT4vbW497SofQ+Jbt7JyaZjrhrqlpxrdsrz2k3WI9w6ulWnpaqslahltLNfW+6Lz9E9dy0sbx\nZf8XeOsZXi29b6C9elra1qC9elrb3voZSofQ2Pp1rBlZxeqAfUZWMbZ+Xe0h7RbrGV4t1dLTUk3W\nMtxaqqm1LzrWM7xaet9Ae/W0tK1Be/W0tr31c/ruENpwyFo2nTbWzFQD6xleLdXS01JN1jLcWqqp\n90XnrqnpJr7oWM/waul9A+3V09K2Bu3V09r21i8ys/YYmjQ6OpoTExO1hyFJ0oK0dNwVWI+0q1rb\n1lqrZ7mJiMnMHJ23n6F0zzCUSpIkSVrJFhpKPaZUkiRJklSNoVSSJEmSVI2hVJIkSZJUjaFUkiRJ\nklSNoVSSJEmSVI2hVJIkSZJUjaFUkiRJklSNoVSSJEmSVI2hVJIkSZJUjaFUkiRJklSNoVSSJEmS\nVI2hVJIkSZJUjaFUkiRJklSNoVSSJEmSVI2hVJIkSZJUjaFUkiRJklSNoVSSJEmSVI2hVJIkSZJU\njaFUkiRJklSNoVSSJEmSVI2hVJIkSZJUjaFUkiRJklSNoVSSJEmSVI2hVJIkSZJUjaFUkiRJklSN\noVSSJEmSVI2hVJIkSZJUjaFUkiRJklSNoVSSJEmSVI2hVJIkSZJUjaFUkiRJklSNoVSSJEmSVI2h\nVJIkSZJUjaFUkiRJklSNoVSSJEmSVI2hVJIkSZJUjaFUkiRJklSNoVSSJEmSVI2hVJIkSZJUjaFU\nkiRJklSNoVSSJEmSVI2hVJIkSZJUjaFUkiRJklSNoVSSJEmSVI2hVJIkSZJUjaFUe8Xk1h2cfflm\nJrfuqD2UJdFSPS3VAu3VI2nXtPZZ0FI9LdXSotZen9bqadVI7QFosMmtOxjfsp2x9evYcMja2sPZ\nLZNbd3DSxnF2Tk2zZmQVm04bW9Y1tVRPS7VAe/VAW58FYD3DrLVaWvosaKmelmrp18r7p7XXp7V6\noJ1tbSZD6RBq7Q00vmU7O6emmU64a2qa8S3brWdItFQLtFdPa58F1jO8WqoF2vssaKmelmrpaen9\n09rr01o9LW1rMzl9dwgNegMtZ2Pr17FmZBWrA/YZWcXY+nW1h7RbWqqnpVqgvXpa+yywnuHVUi3Q\n3mdBS/W0VEtPS++f1l6f1uppaVubyT2lQ6j3BrprarqJN9CGQ9ay6bSxZqYatFRPS7VAe/W09llg\nPcOrpVqgvc+CluppqZaelt4/rb0+rdXT0rY2U2Rm7TE0aXR0NCcmJnb58a3OF5e0OK19FljP8Gqp\nFmlv8/2jvWW5bWsRMZmZo/P2M5TuGbsbSiVJkiRpOVtoKPWYUkmSJElSNYZSSZIkSVI1hlJJkiRJ\nUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElS\nNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1\nhlJJkiRJUjWGUkmSJElSNYZSSZIkSVI1hlJJkiRJUjVDFUoj4sURcW1EXB0Rb+1rf2VEbC73PaWv\n/byIuDkivjpjOfeNiMsi4rry79rSHhHxrrKsL0fEY/sec3Lpf11EnNzXviEivlIe866IiD37LEiS\nJEnSyjE0oTQijgVOAB6VmUcCbyvtjwROBI4EjgfeHRGry8PeW9pmegXwqcw8HPhU+R3gqcDh5XY6\n8J6yjvsCZwKPA44GzuwF2dLn9L7HDVqfJEmSJGkXDE0oBc4A3pyZdwJk5s2l/QTgg5l5Z2Z+A9hM\nFxzJzM8Ctw1Y1gnA+eXn84Fn9rW/LzvjwAERcRDwFOCyzLwtM3cAlwHHl/vuk5mfy8wE3te3LEmS\nJEnSbhqmUHoE8PiI+HxEfCYijirtDwZu6Ou3rbTN5QGZeRNA+ff+8yxrrvZtC113RJweERMRMXHL\nLbfMM0RJkiRJ0sjeXFlEfBJ44IC7XlXGshYYA44CLoiI9cCgYzhzV4cwy7IW2z5QZp4DnAMwOjq6\nq2OUJEmSpBVjr4bSzHzybPdFxBnAh8s02SsjYho4kG7v5EP6uh4M3DjPqr4VEQdl5k1lCm5vKvBs\ny9oGHDOj/YrSfvAi1y1JkiRJWqBhmr57MfBEgIg4AlgD3ApcApwYEfeKiMPoTjZ05TzLugTonUH3\nZOAjfe3PLWfhHQNuL9N7LwWOi4i15QRHxwGXlvu+GxFj5ay7z+1bliRJkiRpN+3VPaXzOA84r1ze\nZSdwctlrenVEXABcA0wBL8zMuwEi4gN0ezgPjIhtwJmZeS7wZrrpv6cC3wSeXdbxceBpdCdL+j5w\nCkBm3hYRrweuKv1el5m9EyidQXeW332Bfyg3SZIkSdISiC73aamNjo7mxMRE7WFIkqT/1979x0pW\nlgcc/z67lxXrL5bFn2Bgt0VbNBRhSy+xbSxYQbRFG0wwJKCVklI1tthUkMZKa1M1VaopxhrAYIsF\npaJEYhEE28a44F4FBBG4rKxsQFfWxWpaWLb36R/nnTpc5i57792975n3fj/JyZx5z5k555lnzjnz\nzHnnjCSpioiYysz1TzZfn7rvSpIkSZKWGYtSSZIkSVI1FqWSJEmSpGosSiVJkiRJ1ViUSpIkSZKq\nsSiVJEmSJFVjUSpJkiRJqsaiVJIkSZJUjUWpJEmSJKkai1JJkiRJUjUWpZIkSZKkaixKJUmSJEnV\nWJRKkiRJkqqxKJUkSZIkVWNR2lNTm7dz4Y3TTG3eXntVJFXkvkBSi9y3SRo2UXsF9ERTm7dz6kUb\n2LFzhlUTK7jsjEmOOnh17dWStMRa3BdMbd7Ohk3bmFy3xlh6qKWYWoqlNa3t21p8r7UYUwtazotF\naQ9t2LSNHTtnmEl4bOcMGzZtG/s3XmsbUUvxtBQLtBVPa/uClj6IthTLQEsxtRQLtLVfg7b2ba29\n16C9mFrZflrLy2wWpT00uW4NqyZW8NjOGfaZWMHkujW1V2lRWtuIWoqnpVigvXha2xe09EG0pVgG\nWoqppVha269BW/u2lt5rAy3F1NL201JeRrEo7aGjDl7NZWdMNvGtDrS3EbUUT0uxQHvxtLYvaOmD\naEuxDLQUU0uxtLZfg7b2bS291wZaiqml7aelvIxiUdpTRx28emw3mtla24haiqelWKC9eKCtfUFL\nH0RbimWgpZhaiqXF/Rq0s29r6b020FJMLW0/LeVllMjM2uvQpPXr1+fGjRtrr0ZvtNKff6CleFqK\nBdqLR5Lcr0kL5/ZTV0RMZeb6J53PonTvsCiVJEmStJztblHq/5RKkiRJkqqxKJUkSZIkVWNRKkmS\nJEmqxqJUkiRJklSNRakkSZIkqRqLUkmSJElSNRalkiRJkqRqLEolSZIkSdVYlEqSJEmSqrEolSRJ\nkiRVY1EqSZIkSarGolSSJEmSVI1FqSRJkiSpGotSSZIkSVI1FqWSJEmSpGosSiVJkiRJ1ViUSpIk\nSZKqsSiVJEmSJFVjUSpJkiRJqsaiVJIkSZJUjUWpJEmSJKkai1JJkiRJUjUWpZIkSZKkaixKJUmS\nJEnVWJRKkiRJkqqxKJUkSZIkVWNRKkmSJEmqxqJUkiRJklSNRakkSZIkqRqLUkmSJElSNRalkiRJ\nkqRqLEolSZIkSdVYlEqSJEmSqrEolSRJkiRVY1EqSZIkSarGolSSJEmSVI1FqSRJkiSpGotSSZIk\nSVI1FqWSJEmSpGosSiVJkiRJ1ViUSpIkSZKqsSiVJEmSJFVjUSpJkiRJqsaiVJIkSZJUjUWpJElq\nztTm7Vx44zRTm7fXXpU9orV4JGnYRO0V0GhTm7ezYdM2Jtet4aiDV9denUUznv5qKRYwnr5rKZ6W\nYhloJaapzds59aIN7Ng5w6qJFVx2xqTx9Ewr77WB1uJpSWu5aS2eAYvSHmrt4GM8/dVSLGA8fddS\nPC3FMtBSTBs2bWPHzhlmEh7bOcOGTdvGNhZoL56W3mvQZjytFD0t5qaleIbZfbeHRh18xpnx9FdL\nsYDx9F1L8bQUy0BLMU2uW8OqiRWsDNhnYgWT69bUXqVFaS2elt5r0FY8g6LnQ1++i1Mv2jD23cVb\nyg20F88wz5T20ODg89jOmSYOPsbTXy3FAsbTdy3F01IsAy3FdNTBq7nsjMlmzva0Fk9L7zVoK57W\nzsq3lBtoL55hkZm116FJ69evz40bNy748S11nQDj6bOWYgHj6buW4mkploEWY1I/tfZeayWewZnS\nQdHTQvfQVnIzMG7xRMRUZq5/0vksSveOxRalkiRJ0lIbt6JH/ba7RanddyVJkiQBXXdxi1EtNS90\nJEmSJEmqxqJUkiRJklSNRakkSZIkqRqLUkmSJElSNRalkiRJkqRqLEolSZIkSdVYlEqSJEmSqrEo\nlSRJkiRVY1EqSZIkSarGolSSJEmSVI1FqSRJkiSpGotSSZIkSVI1FqWSJEmSpGosSiVJkiRJ1ViU\nSpIkSZKqsSiVJEmSJFVjUSpJkiRJqsaiVJIkSZJUjUWpJEmSJKkai1JJkiRJUjUWpZIkSZKkaixK\nJUmSJEnVWJRKkiRJkqqxKJUkSZIkVWNRKkmSJEmqxqJUkiRJklRNr4rSiHh7RNwVEXdExAeH2s+N\niOky7fih9hNK23REnDPUvjYiboqIeyLiiohYVdqfUu5Pl+mHLHQZkiRJkqTF601RGhG/DZwEHJ6Z\nLzq6Z+QAAAubSURBVAH+rrQfBpwCvAQ4AfhYRKyMiJXAhcCrgcOAN5Z5AT4AXJCZhwLbgbeU9rcA\n2zPzl4ALynwLXYYkSZIkaZF6U5QCZwHvz8xHATJza2k/Cbg8Mx/NzO8B08DRZZjOzE2ZuQO4HDgp\nIgI4FriyPP5S4HVDz3VpGb8SOK7MP69l7KX4JUmSJGnZmai9AkNeBPxmRPwN8AjwZ5n5DeBAYMPQ\nfFtKG8D9s9p/HVgDPJyZO0fMf+DgMZm5MyJ+Uuaf7zJGiogzgTPL3Z9FxF27CriCA4CHaq+E5sWc\njRfzNX7M2fgxZ+PHnI0X8zV++pyzg3dnpiUtSiPieuB5IyadV9ZlNTAJ/BrwmYhYB8SI+ZPRZ3lz\nF/Ozi2nzXcZImfkJ4BNzTa8tIjZm5vra66HdZ87Gi/kaP+Zs/Jiz8WPOxov5Gj8t5GxJi9LMfOVc\n0yLiLOBzmZnAzRExQ1f1bwFeODTrQcADZXxU+0PAfhExUc6WDs8/eK4tETEBPAv48QKWIUmSJEna\nA/r0m9LP0/0WlIh4EbCKrsC8GjilXDl3LXAocDPwDeDQcqXdVXQXKrq6FLU3AieX5z0d+EIZv7rc\np0y/ocw/r2XstVdAkiRJkpaZPv2m9BLgkoi4HdgBnF4Kxjsi4jPAd4CdwFsz838BIuJtwLXASuCS\nzLyjPNe7gMsj4n3At4CLS/vFwD9FxDTdGdJTADJzIcsYR73tWqw5mbPxYr7GjzkbP+Zs/Jiz8WK+\nxs/Y5yy6uk+SJEmSpKXXp+67kiRJkqRlxqJUkiRJklSNRekYi4g3RMQdETETEetnTTs3IqYj4q6I\nOH6o/YTSNh0R5wy1r42ImyLinoi4olzYiXLxpyvK/DdFxCFLFd9yNleetDQi4pKI2Fp+4z5o2z8i\nrivbyHURsbq0R0R8tOTqtog4cugxp5f574mI04faj4qIb5fHfDQiRv0tlXZTRLwwIm6MiDvLPvEd\npd2c9VRE7BsRN0fErSVn55f2eR+L5nu808JFxMqI+FZEfLHcN189FhH3lf3WLRGxsbS5X+yxiNgv\nIq6MiO+WY9oxyyZnmekwpgPwK8CLga8C64faDwNuBZ4CrAXupbtQ08oyvo7u6sa3AoeVx3wGOKWM\nfxw4q4z/MfDxMn4KcEXtuFsfdpUnhyXLwW8BRwK3D7V9EDinjJ8DfKCMnwh8ie7/jieBm0r7/sCm\ncru6jK8u024GjimP+RLw6toxj/MAPB84sow/A7i77AfNWU+H8jo+vYzvA9xUcjGvY9FCjncOi8rb\n2cCngS+W++arxwNwH3DArDb3iz0egEuBM8r4KmC/5ZIzz5SOscy8MzPvGjHpJODyzHw0M78HTANH\nl2E6Mzdl5g7gcuCk8i3JscCV5fGXAq8beq5Ly/iVwHG9+lalTSPzVHmdlpXM/A+6K3QPG94WZm8j\nn8rOBrr/SX4+cDxwXWb+ODO3A9cBJ5Rpz8zMr2d3hPjU0HNpATLzwcz8Zhn/KXAncCDmrLfKa/+z\ncnefMiTzPxbN63i3l8NqWkQcBLwGuKjcX8hnB/NVn/vFnoqIZ9J9KX4xQGbuyMyHWSY5syht04HA\n/UP3t5S2udrXAA9n5s5Z7Y97rjL9J2V+7T1z5Ul1PTczH4SuCAKeU9rnu70dWMZnt2sPKN0EX0Z3\n5s2c9VjpCnoLsJXuQ9O9zP9YNN9cauH+HvhzYKbcX8hnB/O1tBL4ckRMRcSZpc39Yn+tA34EfLJ0\nk78oIp7GMslZn/6nVCNExPXA80ZMOi8zvzDXw0a0JaO/hMhdzL+r59Le42s+XubK13zbtUgR8XTg\nX4E/ycz/2kWnDnPWA9n9H/gREbEfcBXdT1KeMFu5nW9u5jreaQEi4rXA1syciohXDJpHzGq++uXl\nmflARDwHuC4ivruLed0v1jdB99Oht2fmTRHxEbruunNpKmeeKe25zHxlZr50xDBXQQrdNx8vHLp/\nEPDALtofojvlPzGr/XHPVaY/iyd2a9SeNVeeVNcPS9cXyu3W0j7f7W1LGZ/drkWIiH3oCtLLMvNz\npdmcjYHSPe2rdL+Jmu+xaL651MK8HPi9iLiPrmvtsXRnTs1Xj2XmA+V2K90XP0fjfrHPtgBbMvOm\ncv9KuiJ1WeTMorRNVwOnlKvfrQUOpfth8zeAQ6O7Wt4quosPXF36ld8InFwefzrwhaHnGly162Tg\nhjK/9p6Reaq8Tnr8tjB7GzmtXAVvEvhJ6V5zLfCqiFhdrpT3KuDaMu2nETFZfmN12tBzaQHK63gx\ncGdmfnhokjnrqYh4djlDSkQ8FXgl3W+B53ssmtfxbu9H1qbMPDczD8rMQ+heyxsy81TMV29FxNMi\n4hmDcbr92e24X+ytzPwBcH9EvLg0HQd8h+WSs11dBcmh3wPwerpvPR4Ffkj3hhtMO4/u9zl3MXRl\nLbordd1dpp031L6O7sAwDXwWeEpp37fcny7T19WOezkMc+XJYcle/38BHgQeK9vYW+h+D/UV4J5y\nu3+ZN4ALS66+zeOvhP0HZduZBt481L6e7sPBvcA/AFE75nEegN+g64J0G3BLGU40Z/0dgMOBb5Wc\n3Q68p7TP+1g03+Odw6Jz9wp+fvVd89XToeTm1jLcMXhN3S/2ewCOADaWfePn6a6euyxyFmUFJUmS\nJElacnbflSRJkiRVY1EqSZIkSarGolSSJEmSVI1FqSRJkiSpGotSSZIkSVI1FqWSJC2BiHhTROQc\nw8N7eFlHRMR7I2L/EdMyIt67J5cnSdJiTNReAUmSlpk30P3/7bCde3gZRwB/Cfwz8ONZ044ZsXxJ\nkqqxKJUkaWndkpnTtRaemRtqLVuSpFHsvitJUk9ExLMj4h8j4u6I+O+IuD8iPh0RB86a70URcVVE\nbI2IRyLi+xHx2YiYiIg3AZ8ss94z1EX4kPLYx3XfLd18MyIOjYhrIuJnEbE5It4TEStmLffIiPjP\nssz7I+LdEXF+ROTefF0kSW3zTKkkSUtrZUTMPv7OZOYMsD/wCHAu8CPgBcA7ga9FxC9n5iNl/i8C\nDwNnAQ8BBwIn0n3ZfA3wPuAveHxX4QefZL2uoitmLwB+FzgfuL+0EREHAF8BHgBOA3YAfwocMq/o\nJUmaxaJUkqSl9d0RbdcAr83Mu4B3DBojYiXwNeD7wKuBq0pxeChwUmZePfQcny63P4qIe8v4fLoK\nfygzB2dYr4+IY4E38vOzrmcDTwOOz8wtZf2uBe7bzeeXJGkki1JJkpbW63nihYb+/+q7EXEW8EfA\nL9IVgQMvLrfbgE3A+yPiucBXM/OePbBe18y6fzvwsqH7k8DXBwUpQGb+T0RcA7x5DyxfkrRM+ZtS\nSZKW1u2ZuXHWMA0QEW8HPgZcD/w+cDRdMQiwL0BmJvA7wEbgb4G7I2JTKWYXY/ZVeh8dLLN4PrB1\nxON+uMjlSpKWOc+USpLUH6cAX8nMdw4aImLt7JkycxNwWkQE8KvA24CPRcR9mfmlvbRuDwLPGdH+\n3L20PEnSMuGZUkmS+uMXgMdmtc3ZNTY7t9D93hPgpeX20XL71D24bhuAYyLioEFDRDwVeM0eXIYk\naRnyTKkkSUvriHKxotk2Av8GvCsi3g3cDBwLnDw8U0QcDnwEuAKYBlYCbwJ2AjeU2b5Tbt8aEZfS\nFbq3ZeaORaz3h+mu9nttRJxPV/ieXW79SxhJ0oJZlEqStLQ+O0f7s4G/Avaj+6uVfYF/B46nu7DR\nwA/orsZ7NnAQ3V/IfJvu6r1TAJl5a/kv0jOBP6TrGbWWRVwpNzMfiojjgI8Cn6K74NLHgQPo/iJG\nkqQFie56CZIkSfNT/rLmm8BDmXlc7fWRJI0nz5RKkqTdEhF/TddleDOwBjgDOBw4seZ6SZLGm0Wp\nJEnaXQm8B3hBGb8NeN1evOKvJGkZsPuuJEmSJKka/xJGkiRJklSNRakkSZIkqRqLUkmSJElSNRal\nkiRJkqRqLEolSZIkSdX8H7ZO93F/gDabAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fc2001db358>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(15,10))\n",
    "plt.plot(x,y,'.',label='HRAP')\n",
    "plt.axis('equal')\n",
    "plt.xlim(100,50000)\n",
    "plt.ylim(-6000000,-6050000)\n",
    "plt.xlabel('Easting', fontsize=16)\n",
    "plt.ylabel('Northing', fontsize=16)\n",
    "plt.title('HRAP Grid in Stereographic projection', fontsize = 18)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "4643.3413071078248"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x[-1] - x[-2]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "4808.8484406499192"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x[1] - x[0]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "tmp = pd.DataFrame({'cell':cells,'x':x,'y':y})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "tmp.to_csv('./data/hrap_x_y_stere.csv',index=False,header=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Notes:\n",
    "\n",
    "Plotting the X,Y data using Qgis it looks like removing the +to_meter field from the projection puts the data on grid that covers the US nicely and allows the dataset to be rasterized with a cell size of 4762.5, the raster was indexed using the cell number so that it will allow the manual re-mapping of the Stage IV data to the projected HRAP grid without too much of an issue."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}