Skip to content

Instantly share code, notes, and snippets.

@flamingbear
Last active August 29, 2015 14:14
Show Gist options
  • Save flamingbear/0c2a617d1e64a753a62f to your computer and use it in GitHub Desktop.
Save flamingbear/0c2a617d1e64a753a62f to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"metadata": {
"name": "",
"signature": "sha256:524be89bf6206cf8844e047497122bc2a6552ae3899fa6310c87f088965d8930"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"import netCDF4 as nc\n",
"import numpy as np\n",
"from netCDF4 import Dataset"
],
"language": "python",
"outputs": [],
"prompt_number": 67
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"file = '/sharehome/savoie/projects/nsidc0081/ancillary/PS_N25km_LOCImask_01.nc'"
],
"language": "python",
"outputs": [],
"prompt_number": 68
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"rt = Dataset(file, 'r')"
],
"language": "python",
"outputs": [],
"prompt_number": 69
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"print rt.variables['loci_flag']"
],
"language": "python",
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"<type 'netCDF4.Variable'>\n",
"int8 loci_flag(rows, cols)\n",
" long_name: Land, lake, ocean, coast, possible ice classification mask for sea ice data.\n",
" valid_range: [0 4]\n",
" flag_values: [0 1 2 3 4]\n",
" flag_meanings: ocean possible_ice coast land lake\n",
"unlimited dimensions: \n",
"current shape = (448, 304)\n",
"filling off\n",
"\n"
]
}
],
"prompt_number": 70
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" imshow(rt.variables['loci_flag'][:])\n",
" loci = rt.variables['loci_flag'][:]"
],
"language": "python",
"outputs": [
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAMIAAAEACAYAAADhtqk6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXdcU2fbgK8DggwVBwLiqgNUBBFZDkhCcNLauqq14+34\n3g6t9a2t1Q4rOLqHbW1rq91Dq63WWhUVCQloVQQRZYpatzhwy5Bxvj8CgchICAkGzPX7RXOe86yE\n3OdZ9xBEUcSChbsdqzvdAQsWzAGLIFiwgEUQLFgALIJgwQJgEQQLFgCLIFiwAJhIEARBGC0IQpYg\nCDmCIMw1RRsWLBgTwdjnCIIgWAPZwHDgNLAXmCqKYqZRG7JgwYiYYkQIAg6LonhMFMVi4DfgARO0\nY8GC0TCFIHQGTla5PlWeZsGC2WIKQbDobFhocrQwQZ2nga5VrruiHhU0CIJgERYLJkUURaE++U0h\nCEmAhyAI9wBngCnA1OrZIk3QdFWUgMzs2vCaXEbcuGW4XLoIwMDkN0n9vq6Buf5tVODcT8Rrchmf\nDfsT30PpWvd+9LgPldMAvh9so9WGJLKUoP+VkrDQmj2fWNdad2RpIVHL3tO7L1GbIOre+n+G+lDR\nhjCj/mWNPjUSRbEEmAFsBTKA1ZYdo0oy1lixxkFi8nae+yoXxSe/o+q4uJoQADyes5FPspYSm6ig\nm7SM2EQFsYkKVB0XMyQ1qU4haI6YYkRAFMVoINoUdTcHXhjnw+8LvOgzrozUGcZfpi3NO82MVd9A\nTt352ty4jjwxgScdQJ6oTvvDayQfrA+qlnfk4iLmjtzJrRJr3o8ZWq/RoClgEkEwD+4x6zYk80tZ\nLNiapA3XvAv1yi/zUP+/pt8opsgHa9JbtoHgWaU42BYT7fQelAvL6A7Kevepog1T0pA2jH6gplej\ngiCafo1g3vSdUIaLj4ibn8iaccZ7Hrn4iCjeX0P/I1n1LnvT3gFZykskfW6NJLKUt0dsY9j+RKP1\nrbEQZpjHYtmCHmSts8LFtxRRBDc/kdyUev3dasTVV0Tx7mq8jmQbVN6xIB/VwI/YnRiAPDEB9je4\nS00Gi9JdIzFlQwmxiQokkaWAekSgDJy6i7j5129U7jJEpNeoMs11p0Ei//vmNIp3fjNYCCpwKCxQ\nC8FdhmVEaCRkQjryxATkHRNQJIbS/0gWrnkXGLhU1/ZpJV1DRH74OI5uuadpeauIp+99kgM/WaGI\nWkXfYzlw1MQfohljWSOYGJ9Hy3j+oQM8e+wvrfQZts9w8KwLh9Zb6TUtkkSW8p3fL/Q6dUyTdsrV\nnZv2DvQ5dtjY3W7SWNYIZsY3m/cxPm8H7Y9drnZPFCE+Ur+9eklUKSrnxbedz0OXc2eM0U0LWNYI\nJuWz1QNJ9exPnlO7avfiVa7VC3iU7/916QI9e6rfd+9O1Bil6TppAbCMCCblRILAxB1yOj8zk/bB\n7YiyX0DY3h3qmxMnQhtXWL8ecnPVaY88Uq2OyBHzKstU4ZSrOwUt7fA4YVkYGAOLIBiTQYNAJoN9\n+3DqLnKlRxgAlwsAJcgZzucjn+X5fcvV+bOzK4UAYOVKkEjUIwIgD99KVMpbNTa1/JQXwtErLOhj\nEQRjYBGEhiKTQUAAtGoFwEC/JPa3kXG1luwztn2NOFJAUqok7VI37ZuHDqkFAejR8zDfnX+61maf\n7ppJQW87OG6Ez2DBIggNIihILQhVuHjBRWexF7Z9hYtrLqR+pX3j6aehc2d69sohzkFOt7Onaq4A\n6Jp72pAeW6gFy2K5oXzwASQlQXIyfv57OZXRVr9yogC+vpXXnp7QWW3I91TPb+oUAgvGxyIIDSEx\nEW7ehL59aZUdR8rYTRAXB1drmxhV4uKaC+PHqy/69IGHHwagV+9D/Ofsr6bstYUasEyNGoj/tDKS\nP/yQG1UTnZxqzT9h0ipk1+JZfvY5dUJUlNb93r0P0TXJMu1pbCwjQgM5kVA/ZbmsTG+mHPujxnu9\nPbL55uSzxuiWhXpiEYQGciGtiiAMGgQP1O25JiPdh2kDljLfe361e709si2nxXcIy9TImNx/v17Z\n1v3xEId8+mmleXhks+L4c6bolQU9sIwIxsLfv17Z+/TTtiP+1X0qXc6dNWaPLNQDiyAYg2nTYOxY\nvbP380rj+evLcHWz/PDNBcvUqKEEBoJrDQp0ddDPK42LNu05l9vJRJ2yUF8sI0JDqOFkuS76ex/g\nt0cmsOzATK10zz6ZdL5gGR3uJBZBaAgyGTg66p09PW0AFy90xOXSBR6M2cAA3xQADmX343RHy+hw\nJ7EIQkNwcKh3kS9Pz+B8e2cTdMZCQ7AIgqEEVXeCpQ8Z6T6ca9/RyJ2x0FAsgqAPc+bA669DmzaV\naWFhBlcX1UX7MM3DMxP3C7m15LbQGFgEQRctWsD774OtLVy7pk4LDoaCAoOrXPfHQxzw7M+qllNx\ndTvLM/d8TefzlsXyncQiCHXx+utgYwNvvAGArXO5i8YxY6B9+wZVHen2Jl5Hs3FzO0NYobKBHbXQ\nUCyCUBe2tjB3LtjYMGv22wgAr75q9GainKrrHVloXCwHarVRPgpUkJwUTNGM8DvUGQumxjIi1IRE\nop4SlTNkaALxSuMKwfp1U1g1ZgKyDnEcOtSPD/xe1Nw73LUHUu8Y4gJDjNqmhdqxjAg1cdtp8ZBh\nKtraXiZaqZ92qb7Mivlco2YxJ3sJbk+cocu5M7x25R0unHRlXqfFvB3wBtKknUZt10J1LIJQE7e5\nwVy06m329RtAbMtR3CpqabRmggfvYMP6BzXX//lhtdb9o0c8OPp/3SkThBp9G1kwHpapUU2oVNWS\nQlL2MGKIcYMAVRWC2wkfoW7rqW9/RR6dQGyw6cNN3c1YBKEmpNIakzemjcfJ6UqjdOHfo7017+eP\neJPwPfGN0u7dikUQakKo3Q553v/Na5QuHD2i9oMaOWIeC1IWN0qbdzMWQagJq8qvZciweGxKSjTX\ns3/6gvtlaxulGwtGvF6ry0dzJS4wRPNqSlgWyzVRVgZWVgwLVRJzdAw2JcVat+dav0eM/RgKCuqv\nfaov4SOimZ/yjsnqNxaReaFIR4kkxAiUloKwrgTs1O7ub4WHMOpA01jkWwThdqpMi3YmyJg/93U+\n+Fb75Hdo6l5GDd7E+rjaF7sNoWevHL7LfcYkdRuDMisrVP5DUUYLCFEJqMpCABHZvQKS5L20KFWP\noLcy9Ykaah5YBKEqMpneFmdzhffZ5nAv+fnGHRV69T5EnJ2crmfN08lXqZU1C88MRpi+AyJDkOwJ\nRb63PObabQE4bYtvNX4HDcSyRqiKUlntDGHhL2/XmHXwgSTsHW4avQv/12OFWTv4jRs4BOHdXRAV\nwoKOOyqFoImjUxAEQfhOEIRzgiAcrJLWXhCEGEEQDgmCsE0QhLZV7r0mCEKOIAhZgiCMNFXHjY4g\nqN0vlk+NhoUqiY2QYF9UWGsRiU+cUbswOmIDr6V8aNQ6jUmMXwg7ZLuQ7RhClHPTmPvriz4jwvfA\n6NvSXgViRFH0BGLLrxEEwQuYAniVl/lSEISmMeqIonpEKGdngoxi0ab2/IDUJ9ZozTcFB18iAmGq\nYGTNUOVD549UFMUE4PZoePcDP5a//xEYV/7+AWCVKIrFoigeAw4Dhtk03gmqCALA6OhYClra1Zo9\n/qDcaE3/956vzd7B18iUBKTJ/9zpbpgEQxfLrqIonit/fw6ocOzjDuyuku8U0NnANhoXQYDI+oW8\nVRlBEDz7ZPJ096+Zve/TBtdlwXAavGskiqKojptce5aGtmFqJFGlxC/Q/iqGhSpZ2Ho+9om1rxHq\ni1f/Azg7X6R/izT6WGfj0+Ig/Y4eotO+c7oLWzAphgrCOUEQ3ERRzBUEoRNwvjz9NNC1Sr4u5Wk1\noKzy/p7yl/FwcBHJP6+fy/b4KGtAVC+Wy+MV7EyQsUA2H0mLMbSocrJcwQSf1eTF1eyNQiJV0Lbd\nJeJVcq5cVpt0Lhv1XyYeXU/H9DxDPo6FOlAeAmVOw+oQRFH3A1sQhHuAv0VR9Cm/fh/IE0XxPUEQ\nXgXaiqL4avlieSXqdUFnYDvQW7ytEfUIUr9pSKNxW+CO2XMXVTtQ2zUggPBdqmonyxKpgi5WJ/ny\n2EyOdL6HK20qA4bIE5vHNmNTQJgBoijWK3CFzhFBEIRVgBRwFgThJDAfeBdYIwjC/wHHgMkAoihm\nCIKwBsgASoDptwuB2SOK1ZTuxvhs4K+sidgWq1Ut3hPn1qxeIZQx2nYLTtevMSjrQGP01oKR0GtE\nMHqj5joiSKW1+iuaNfttPv5BbcfcruQSV660qzHfE//3Fd//Nc1kXbSgG0NGhKaxx9+EWLjHojLd\nFLEIQlVqGA2GDFMbxCxco1a12OUbSGFhzWcLT/33S7NWj7BQOxZBqKAWY5yKneH5k18HYKdvEIWF\n9jXmlZ2yLIibKk1GEJy9TLiWkclqPUz7Z4e22ebib9VTn8CgXdXyhp20mFM2VZqEIHTsL3Ixo35h\nXPXC2lq9XapD9Vo+fAvv/hrFP75BFBapp0Vz2r/H+LA1mjy/jJ5iiYjZhGkSgnAh3QRCAGoHv3Uw\n57UFSKQKtqTdj21xMTv9gigqXx+E5aif/hKpgpWjH+SRpDV1VWXBzLm7DXOKinRmUaVXerhTJanf\nrx1zPx32XmbdwSkm65qFxqVJjAgmJV57Xl+xOB4WquTt7xZp3Xutxbu0bFlInGdoo3XPQuNg9iNC\nv4llZK41kbxaW6v9nFYhJkKGIIpIk3ZiXVamdW/Y/j3Y2RlPCc+C+WD2gmAyIYAa1whhiQlY1XHa\nHuoXx/t/WNy4Nzfu3qmRvT289ppWUqhUga5l+d9pE+o037TQNDH7EcEkjBgBw4ZVS5bKYhHSm5aO\noAXjcPeNCDJZjUJg4e7m7hKENm1qPTwTBBFZM3FNYqH+NJ+pkZ8fHD4ML78MGzdCq1aQnw8dOqht\nDAoL6zxBlkhiCbcYz9y1NA9BcHKCBx6ovL7vvnpXIZXGQroR+2ShSWGWgnDfihI2Pl1H19q3h5kz\n4fJlOH68mne6+mJtXcKCL99tUB0WmjZmJwid/EUO/V3L0kUmU091fHzU1+3aqV8NwMqqlO0jZJDU\noGosNHHMThAADm2oQRA6dtTbQa++2NoWsU0ebgnWZ8H8do06BYp43q+t2oCLCzz/vNHbunWrJfML\nowj33mL0uk3J4a497nQXmh1mNyLs+6oG2XRxMVl78crhAEhlMZq0uPRRWIlltRW5oyQMGszDmatR\ndAvD48TRO92dZoPZCYIWrq4wrXE8QlQIREyE1ChCUBE6qaawsOm9+uJ8JQ/XvAu1lj/StQcrnacQ\n0mIHqhIpccUyQlvsYFfZYE6d7MajnX9iD00rPJM5Y5aC4DG2jJy/rWp1rWIq5kW+wfAvGmZuGRcY\nwndFT/FL9JMA5Pj1pPfJfzX3p/l+xtqDkxkWomRmwecaQcm+pzd7C9x49Jz6+nDXHszfoB06Kp7h\nDeqbhdoxP0Fwd6fPcyXk/H0e+va9072pkz0+/uTb2ROYnkJWDw9eFj4gPlr7xzqhbC0zBi0lqTSA\n7FJP4mPV99evm8J6ppAc6sugzAP0OXYYR9fKwCO9T/7L0BBVNZvpCqStqseCtmA4Zufgq9NsKf5T\nL7PxA687IggS2Xa91wjOnCfvYkf8A/Zw+nRXcs+617u9Pn0zcHWrbuv87xEPTp7sXmu5TeNGELFj\ne73buxswicvHRuf0acCh0YVAKt3OfMeFDI9WUeYiYFWqu8xGnzHMLXufmzcduXGjtUHtZmd5kZ3l\npZXWucsJTp/qVme5BblRDHDNoMu5M5zp6MYHnV/ifpsNNa5JLOjG/LZPnQ4hLarfPL2N05V6tyOR\nKhgui2ZLRDgSqYJIhwXIExOIHS2hRakeUgAIwJjgDSQnBXPjumGCUBMjR21mXuTrdeZJ3D2MsFsK\nZg98h7igED7Z/gry6ASk3jFk9OpjtL7cLZjd1IiHH2berz+weEHNQfxux6ntZb4f+RgT1mysO5/T\nFfp5pbF7VwgODvncdHCsZ6+1meCzmq177iU/v2H11EXXbsc5eaL26VFNuLjm8tugB+/qkaFZ+D6d\nNW094+I2653/6pV2LMyJqjPP7LmL8PVLYvcu9XZjfr4DC559RSuPMkC3jcI/vkHEBYYwzud3/oyb\nbFIhADh5ojvduh+rV5nz59w0I8NzvktN07FmiNmNCJ2WP8299/3JN8tn6KynbdvLDPBN4Um77/j0\n/Iu0aXONGzdasy85UK9+LHp9NvOWf0SJtTW254uJHS0FRKzEMsrKYyBGFVT2c8+uUIqKWupVd33x\n8MgmJ0c9penZK4cuXU4S6RjFPWdOIL8cx/Hjhp0mPxP+Bf7Wyfx662EAPi15iYFZB3WUato0i8Vy\np+sKFiS8RXyfcA5l99O618H5AnkXK6PUfDzxeZ78axUA41v9jdPJa5zp6EZndAfls7MrYN7yjwBo\nUVrK62/OQ75QvTaxsiqjrKzxBkvPPpm8OTiS7Id7s3jB2/zc41GG7q+M3t29/78GCcK8yNdZvOBt\nJLJYelkdQURo9kJgKGYnCOTkcK1Pa/p0zqgmCM7OF1gTNJGcUg9WFj7CwOw0zT2nG9cAcL+Qy9PP\nLWXFVy/U2UxhoT1Rs+cwNHkvbxW8TvzCyv3/xhQCgEPZ/Xgsu9JT3mslbzHSL4Z80YGePXKIX2FY\n0MJ45XBiIyR0uHIJ30MWY4u6MLupEcDADffxmPNPvLztM630Ds4XUHiGMUDHH/WGgyMyl1iSk4KN\n0t/bqXjSmjOfvfAML6xacae7cUdoFotlADZsYOqJNfTpm6FJkkgVxHnIdAoBwMZ+o6ka6NPOroDY\nCAnrx9xL9APhdZTUzboxY83enUt/79S7VggMxTwFAeh08RzPdltG337qH36/Fhn45GToKKVGFASS\n9g7WXL817WXkiQk8sHczo3cqeOkVw6PavHdlLm+8+7HB5Zs7JdbmN9vWB7MUhP3fWPGi3X+Zte9z\nlGVhTJj0G7JOCr3LC1WmexKpgsEHjGd+tmeX+Wt8pqf5Mt33U+ICQzjo4aW7gJ6UWLdA6h2Dyn9o\ntXsq/6EogkKJGdA0/cKapSAArP7CncyenrjmXWDZgZk8tOVPvcuKVaLfpB30ZWjqXgCKbG2JCwxh\n9z+S2oo2G5bFzmTy3nUstTeOQdPmgBGM8tmIDSVIk/+pdn94rIpboi1jUuKM0l5jY3bjWCd/EY/7\nyojfEEBGzxT6HT2Ey6Xa9fZrYmrSH7hGSLhY5syqoqlQvmM496EFfPrxq1p5JVL1SBOvMmxnxpy5\neLEj0i5xcAAuObVjdvd3eazlz7S9fhW/emyj3rKxwV4oINJ6AZL91SMFRU2by+yyhYz+OpaEQUMo\nsbYG4IbYipEH4hAFATszX1eZnSBMff4UD9+TTMDXnVFdkDGRvw2qpyLA92TUI0mhbUt2/xOKRKrA\nXshnjsP7mnxlghVh0m3Eqxq2kDY3nDuex/1CLssHPcnnuS9wUOHH9zxDt+7H+LPfA4iCgH9GqlaZ\nuMAQ7IsKtaaTrz4SxVdfzCLfUR1bWuU/tNqosDMhjLjAOMYoYigqqgy2OHjIDgKDdvHZr3NM+Ekb\njtltn3YKEPn1CwXPbX6JQ8J4ls54msxMb75IfbHB7RbZ2iKIoiZweAUl1tbYnCtpcP1NiTZOVxBF\ngXsDNvBMy+UA3LR3YOy6rdjbFzBi8Gb+OjiJWza2vPpIpEYQbtnY4njhJsVtbbTqk/ffSpxqJBKp\ngjkO72EvFADwcf5LXCzsSNzJ4Y2222bI9qnZCQKAJKoUt4c8WSO7Cc89B8CqRyfUa51QH8oEK2T9\nYkiIbz7TI+eO51kdOIkFN6PIu9SB9IO+9a7jvrA/6T0oi08+eg17+wLWycbSwrqEMVu2awShwiR1\n+FYVQ4YmoMwZTotS7YeKIihUM0I3BiYRBEEQugI/AS6ACCwXRfEzQRDaA6uB7sAxYLIoilfKy7wG\nPAWUAjNFUdx2W511CoIGNzeQSqFfP/p5pZFx3qc+n61eiIKAzCum2UyPOnY8z+qAiWSX9mHatm8a\nVNeQYfHsTRxCcPBOBAFsrG4RmzEaWf8Y4pXq72trhJzw5His9VRhNyWmOlArBmaJotgfGAw8LwhC\nP+BVIEYURU8gtvwaQRC8gCmAFzAa+FIQBMN2p3Jz8fZOpWPHc2RmePOs7+cGVaMPgigilcWarP7G\n5sIFF2ac/IJVRQ/T3ydVd4E62LVTQkmxDTt3yLCyLuFNh0XMe+4N4lXhbI+QEhshYWRinEYIdvgN\n5paNjY5azQudP1BRFHNFUdxf/v4GkAl0Bu4HfizP9iMwrvz9A8AqURSLRVE8BhwGggzuYWYGil4y\nlo58jqySfjzj+4XBVd1thLgmoEoPx8dnf4PratFCva6KVw5HHp3AjvgwFGMkhCfGa017EgYNIUwZ\nz6uPRDW4zcakXk9qQRDuAfyAPYCrKIrnym+dA1zL37sDp6oUO4VacAwibaUVGb36MmPf13yVP411\n+ydzwLO/odXViWzvDi3VjKaOgPE+S0mJ+gkfLttCbIQEVXo4YTXM+0eoFJQU27Bo9VtGa7sx0FsQ\nBEFoBawF/ieK4vWq90T1QqOub72Ge8oqr2N1tr1gpTqMq/OVPMb7/q6XvpEhKANDqOfU0qwREbjY\nrgOFhfZGqW/oMBXSoduRJyYgCgLKwMpT9rjAEKTeMZSWNv6OvPIQRG2qfBmCXr0WBMEGtRD8LIri\n+vLkc4IguImimCsIQifgfHn6aaBrleJdytNuQ6ZfD11dwcuLuG4hhO3dwYrL0/UrZwCyvTtYLIjN\nQhj6eaWx+EQkf/QYx+EcT6PU2cKmmAXvvcf20NGIohUJW8KQSNXrqqpubF6eswjH7/KN0qY+yDzV\nrwoWRNe/Dp2CIAiCAHwLZIii+EmVWxuAx4H3yv9fXyV9pSAIH6OeEnkAiRiKiwsXLrgiT0/g1/9M\n4uHNaw2uqi5EQWBhfvOJlpmZ4Y3MW0HGtgFGq7PCG2DF/7e/b8roMyIMAx4FDgiCkFKe9hrwLrBG\nEIT/o3z7FEAUxQxBENYAGUAJMF1syGHF+fNc+FwF/U2zLqjgzWmvo1rYPP6oFWSkGU8Imjs6BUEU\nxR3Uvpao8ZcjiuLbgHEsV0QRryllhMrjTDIalFpZIffaqmWhZqHhxAWGEFUQSZT9gibhUcNstU8B\n9fpg0iSk8vN8lVq36aWhlFlZN5vhvS7atc/TK5+9fT6xERKCgqtrmOqLjU0J8ugE4pXDicpfQGFL\nO92F7jDmLQiAdycFX6b+z2T1x/mZv32BMbh8qYNe+QoKHFiQP59328/lzzH1j0X3yqsLtTR8X7L/\nyOw1T8FMdY2q4v1wGQeHLjJJP7YGyhkd3XxOk42JnV0hscMkXLR15oFo/f1M1cTlHm1pe/2qkXqm\nm+Zjs1wFifSc7kwGohjY/A10DKWw0I53yl7j/r3RxEZIkEgVDB1WtwfuYSEqWras/vSP62P+Vmtm\nLQiuA0Smd2tYvILaKGhpx/vv6Dcq3a1sjBvPfT5/0rL4Fqr0cFSH1e5htkbINQZNFUikClQ5NSss\nvnv5NQrMfJ1g1lOjP7L2MXG7YYY5uihoaYfD6QKT1G3O/Ef+LamX/EjdP0jvMja2txgyVP1AqtgF\nUgYMo8yq8jkqT0wg3n8I4bHxlJRob0ZapkYNwHWgiDLezSR1F7ewYXiPphVA0Bi0a3eJDfsm0K37\nv8x88X1at7mmV7niW7bEK4drFO6k3jHItyawMP9NBFFEEEVU/kN5s2hhNSEAGO26yTIi1NiojhHB\nzU9Esfg3+v17yCTtv/zkIj7+YJ5J6jZ3JFIFqfsH4TswmWvXnNifEmCcemXb69yGlsi2E334fhwK\nTT8KN4sRoZO/yJTnT5tMCJqanryxOZDqh0S2nR69DhtNCKB2VYvN40ZgZ1dIvHI4r4wz3J+UqTGr\nEcE9UEQRuZI+xw6bpN0iW1vC7tnOrn/MfxejMfAbtJeUffp5DjcUiWw7ibtDKCy001xPb7mMKcnr\nTNZmk7ZZlkSW4h4osuqY6fTYrzm2xum4fvPi5s7YB9bSs+dhPl0yt9HatLfPp6BA7QnjTF83Ol00\nzdZ4k5waSSJLmbqlhKgIJWeSBLLv6W2yttrcvM6GMWNMVn9TYv++AFL2GW9qpIs/Rj/A28+8BMCq\n0ZNMJgSGckdHBElUKV2HipzcKRB0z2lEUeDD/O9M2vYtGxtkPWLv+ulRG6crXLvattHaGzXmbwoK\n7JFIFRpnahv/HU/rmzeM3laTChTSTSLy6DPHWf1Ld45us+Kb11R4nDh6p7pz19GYQgCwNXosoL2o\nvubV2iSCYAh3bGr0H1UxT6/7GVnXo/z4aVyjCYFtcTHbT4xqlLbuGhYtgjLdcamr4h+whzY3r+vO\n2EjcMUH4SWKDIiiUeXkr6+X86WLb9g023k/28mX2XNMo8t2VSCSwYwfE1q3AGBC0C4lUwaSwVWRn\neTFr5LuN1EHd3LE1wpH13el56nijtw2Qb2eP46nGs6k1Vwb6JRnvLGHxYigpgdBQsLMDQYDiYrWD\nNiAwaBfKM+GaA7UpPr+wJu4Rfh71EFuLRxFuE8sTyb8apStNao1wp4QAwKGwgK1j5Sy4HMk/O6R3\nrB93GqMJQXy8WggAEm4b3f38CBqRhvJUuJbv09niR2xuPZbHtv4GgNUTxTyBcQTBEO749umdoLiF\nTY1CECpVIKnysrYu5dU3tA36gwbvrKZ52VQZ6GeEACoqFShq/z6CB+9EdZsQAOzt78eN62001wv2\nLdb4Ub0T3JWCsNVPTgtrbR+dgiAS6RilVjcuf8UOl/DO14twcNCeRq0/OVGjo99UeUT+A6rLcvwG\n7TWsAlEEpRLi6g4Msi1/So0WaovWVq7RZr74PmuHj0UencCRrobFk24ozVYQEr1rVzMWRJGSUmvN\ntZVVGW+++Trhty3aa4oMAzDaZSOCKPK+1RxGR2wwTocbESenK8wq/gRE6u/ZT6VST4WUSvWrDl48\nU0ybG9WTNhSqAAAgAElEQVR3hhK9/bhxo7Xm+rNP5jD74y8Y4LuPdtcu168/RsLsAoUYgxsOjmxM\nc67V4erI/XEIfiL3op4alZVZsWjRW8hGq70t1OZ1IXjwTlq2LAJAHp2Ag8NN8vMdNXGdW7e+xvUq\nw73ZIsC1Vq15ccR77PtG/S1JpApu5juQXCUIYzViY6uvAWpBElnKe39/WC090XsQ8r0q+vTNoFWr\nG8Sr1EY+9/Q4zCfxc2h3rfHsFqrSLEeEVvk3mWu/gz0+/jXetykpJuyA9h9UFK2QRydo9O23BYbV\nWHZrTgSq9HBGjfmb/t4HAAg7Fs/ksF956ukvjftBTIUIbW5c55Ptc1n/YAQS2XYQypjT4f26y+kh\nBJJI9ZRzc7ePsS2+Ve3++1av4B+QiPJC+fSz3I/qp6o7JwTQDEeEDX1CyDoo8ET7TKKfPUbwjJrz\nzX/0dbWPviqEShVEOUZprsN9tyAZHEv+W2pFsT27h7Fl9HDGJW1i756hXCr3DPHRxdlkZ3txI65y\nuHdqe5mrV9oZ86MZnTY3rvNAXDQ3A1rxyJY1xNfspkpvov9R4SAUcX2kHY4pNW9PvyJ+SHB6sua6\n4gypMS3YaqIZjggCN7Nv8PmhfpCVx7sFunWKrKzKiImQkqCSI9u7A3liAvLEBOa3WMji9yr9lG0e\nNxxREPj4sWlcutSBgMDdODreJDkpmBvX1UIw0C+Jx5/8mrZtrwDgM6DSJbtEqsC543n6eaVpdqYA\nPDwzjfkF6OR+/3UIiFxzbI3UOwaPgsOsGj2pQXWGzitlxP4EJMm7uLeOyJrBB5NrvXcnaXaCEHYy\nhYAHHbGKSoC+HfAJgH98tXXu8+3s+fC9N7G2LtX8GFuUlRIbIcGqygHjW0VvaJkeLjwXRVvUP/Dp\nL3yEc8fz3LzpCMAgf/Xuy/6UABYkVRqgHDwwEIlUQWyEhM+LZiAAq9tMRpUezrKC6Xw0ciZKa7lm\nF8q980mdn7GqEBnCz4qnCDugYmyPddy43obZwvu4Cud46RUdKvAhtW9vhkaWEjPQNIqMcYEhJt9a\nbXZTo9b5Nwg7kQKxoSRPSaaDVT6DU2t+CoWFbmNLxn3Ihm1nQWEkcWkjte5LhiqIUUZornfvCkFO\nApey29HyVhGPevygubcvWS1sbdpcJaXfALqfO8rxYz3wGbCfdacfpMOVS6yMmIik63Z84tQjgNfR\nbLzIBsD9Qi5ywnll4NusEh7h9Klu1forkSqIdIzSTCfeHfEy0bciDAqNe+1qWy0FODnDQZcXaaH2\nw9qw1H8owfhexJe2CCUvWP1522/rj6+JQgI0uxEB1Itl+bFkxOcHodxdPTaAQ2EBs+cuQkDEuqyM\n+a0XolKF1xjgYciw6u5kRrtsZEzvDWzdcy/+AXtwdT2ruXftmhNPbP0VSZUwVGt7PgCAz6EM4pXD\nWT1qfK19/2D/6zzS9+dq6Y89/g2q9HAtvaxHz6ziPttKLx8vhH9surONkpI6F8txvkMZvU9p1CZ/\njZhE3nOVbWbeNN2aq1kKAoBjQT4LOu7gVUftrdCNAaMY47OBf3ZKEMufYOF745FIFNWeZ/OWf0RK\nciBDhmr/ABJ3D8PbJ5X8fEdmO3/AEr+ZWj/AHj0Pa57S66wnMPHoXwD4HM7kt8BJfHlqBlLvGHKd\nXar1Oy4whJXZj9Gj52EkUgXduh9DIlXw+EVt4Tjl6s4jHX9i462xvD9yFm5uZ3ju+te0a3+x/l/W\njTpUoVXlTr106KTtWGxNqbV1nXnqyyOb/6DDshAkiaFIEkPp52i6M4ZmNzXShcJbwpYPx2JnV0BB\nK/VukJVYhjJ9RLW8I703sdb+AUakKEmIGMy4hI2ac4Ivl74MqMMzPZS0liuTWxOvktOmjXr34/ix\nXpp6Trt0osOVS1xs1wF5YgLe7R8irXdf3C6e12rv60lP8NxX3zNi5Ga+Ofss3dJPMW9gJHJbBXFF\nYYRTOTp1OXeGFXbP4nn8CGdc3DjdvzNeqYcMD9L1wQcwfTrs3QuFhbB7d+W9nTshqO4weGIplFpZ\n6xVVM9Hbj6C0FJ35AF4o3QGJMMN3Cb1bHcYX00yN7pj2qdiAAJl5bdvT4colg8r+HTia+6OjadGi\nhFi5BMm+XTXmUwUMo0CwY/TeWMb4bCBftGfPrlCKilpq8gQE7kaZG06i9yDu37EZD49spob9yJyP\nlmry+Pim0K7tJSIdF5BV2pfpySuqtZV1jwfPtvqSM6e74u5+ml/y/kPX3BqCDOnBuNDV/PXnZIPK\natDj1Lgmir5+p8azg6pUHKj5B+7W20Jtuu+nrD04hW29w/VaIzRJm2VDOO3SSet6Xz/9A2KM2h/L\nkKEJlJS0YPGt2n0bSZN2MnpvLKN9/mZL3FjilcO1hAAgae9gInr9hTw6Qb19KoiaAH6+A5Np2+4S\n7dpexr3zKeSJCRohSO/Vl88GTdPshoSdVtLK8QZxtmGo0sMNFgKjIZPBiy/Wq8jrxbdQ+AzVme/e\n9Ghu3mxFvHI4G/pG6MwPoDwr5/x5VyJdTRfRqEkKQp5TO1L7eANwpbUTgzIP6F22qoVajDICqXcM\nqoCa/4AjvTexNU7tGn3W7Ldp3foaL7z4PrEREuzt1QdG8crh/DVpDG6dzpCSHMiUrevp2u04Tk5X\nGTBgP5GOUcjytBfc/Y9kkdB2mOYkO/esO5s3jeO0i3u9vwtz4W0bW+QH6xdT4eOLL3PNsXWdeVaN\nmUBmhrfm+suWoazsZPyt1CYpCHEX3fnzudMsEUNZEtedbyY8ZlA98yJf53WHd5hfuIAw721a94Z7\nR2ttncrT4pH4K3h/zXwEUcTevtJjmyCK3LjRmp9GT+V8O2euXXXiCbvvNbs8k46s1+S90M6ZuMAQ\n4lXaDnM9PLLpfP6MQZ+jKrKbRnKa3LYtDBlSryIjzs6hxLruZedG7wgcHdXToX3JQbw44r0681dl\nwbmFyDqfw62r8bdpm6QgPFqWDEcuc/WFBDh2hf+uq77dqIuRsk0MS0wk3n8ow0LUu0ZVid+p/UOV\nZexgY9p4Pn/ov8ijE7h0qQMSqQJHx5vcv3YLXw37L48l/Uayly9Xr7blqa0/sy78Pi45tdWsZ571\n/RyZrQJ5dALt2+cRGyHRvOJswuhy7iwNReloRFf3Mln98ouwfUDdT+vgg8mMCfpbs8u2f18AV1vp\nVlQcP3E1vofS8Tp6qF6mvfrSJAXB8/gRJImh+GwORfp+93qfOiZ7+bI17T5G740lasV7vLN4IfHx\n6u1OZcAwygTtr6Xq9mlRy5YEBf+jURZ79yn1XFrZLZRrjq15Zu1PdO2mtr6buPpvZK0U5Dm145d7\nH2R57PNkpPsAavVnr6PZGnWOzucbLgRgxBEBoGVLGKp73q9BgLHhUgrs6o7r/PvBh/lQmE2r1tdI\n2RfIX33ubWBHG06TFARQK2tNPKr+EV25LLB84n/0Lhu6r3JrsCJ0lEo1nLjAUORb45n37BuUllV+\nNYIgMrJLNFLvGOa98xFbLtzLkANq667+h7NwdT3LN8tnMGuE2hj9z+4PaFQmDqb6MaHrGlYcf0ar\nD5kZ3uR2cOGss6vB30FNvLjvC1Y+OtF4FdZDEEqLoKQQ4vrWocpdTmBaCvcGqG05Ht+6Eql3DN/7\nP1rtoXahnXbIq8weniZRt2iS26c1cbVVGwRRrLeLkOIWNtjm1r3lB1QemAkiIFJQ4IC9fX4157c/\njZrK/dmbcbpxjZS+PoSnK7l8uX21+lY+NpHcc+6ICLy0b2m1+w3lt9HjWZSyQDMCGcytW/B2/QOk\n3vj2AxwL6naQ8NWDTzBt2ffV0iWy7aw7NYVfIibx4mdfa6VnZ/XnXG4nJLLt3GeziVdSPqlWvkkZ\n7xsbpxvG8WkqCCKCIDJkWDy7/wmltNSav8eM5r69WwFoW3yZq1WcY/kH7OF9l1cAGJewkf9sXcVA\nvyQ+7jeLdteuUJv6zZqbk/lz30NG6XNN/FE0qeFCAGpBMIA4z8Hcl1q3uod/Riqt21zj+jXtNUK8\ncjgyn1jOrqw8HeznlUa7dpeo0IOJVw7H8d6bvEJ1QTCEJjs1MhY2JcVsi6g0whFFga1j5Kiyw7G2\nLmHTmJEaIfjgiRkaIQgM2sXq0RNQnZPT+uYNFuTP1/zo96cEII9OICxDxeVL1UeDCmpSsTAWX6b+\nD6/+Bxte0T+GhZl9e8NQbtnY1pknMD2F+/zX13hPRCDvYkf1hUpFptKF7Kz+IIBnn0xiIySs+Hea\nQX2riSYjCOm9+jL/YqhJ5ocjEpW8MX8eoVIFUul2ZEk72Ok3mNLSFghVpo6Lvq5UU57T/j0mJ/1J\neq++hCXH84zdCtaH3EdshIQ2TmpV7Ss1TIkqWL9uCmGCAql3DB8NmskZF+NGB3K5dBHlreqxzhqL\nXR9YI82ZQ5FtyzrzrTz4OI/If6iWnn7QV/0mLU3tIECpJEvpzLncTnj2yTTqBgPoEARBEOwEQdgj\nCMJ+QRDSBEGIKk9vLwhCjCAIhwRB2CYIQtsqZV4TBCFHEIQsQRBG1lp5Pel/JAtZBCiXXOLn+xqo\nQlADi5a9TXx6OFJpLPI+2wjbHk9pqTViuerxvT7rteyRK06Q78vYzIqQJ3kkaQ3yxAScrl/jiSeX\n6/UDzMrsT7xyOLO3fcrTPZYZ/TMB5OU5G15461aDRwSA3R9aI8t6RWfA8V8OPMmj8u9p21atVFdh\nb9GhwwVIL1epyMiAZctApWLjhgl86Gfc2Nt1CoIoioVAmCiKA4GBwGhBEIKBV4EYURQ9gdjyawRB\n8AKmAF7AaOBLQRCMNurIExPwfqQDlzYZ70lQgSCKxAWGoIwPJyFeTlmpNa/Ne5OIvTEkDBrM9p2V\n7uQ3jR+p0VyVeisYm1GpyL/E5kXCzqhYUjpL84fVh82bxtVqY20oHS9fRFkkp79PqlHrrQ+7P7Ym\n2kP3KP7zgadQeEs129Kq9HAUnjKcn5dDx46VGePiQKXilZhPkHrH8J7fS0aZJej8kYqiWLH0twVs\nUC9X7gd+LE//ERhX/v4BYJUoisWiKB4DDkOtziQM4sHjCcg8DFO404UoCLzpsEhzyLVoxTuAWkhs\nbCoXjUUFLZmUpFatni1+SKv8mwA8POBHTpZ2I+rQQq62bkPbdqbpZ334o8e4ymlGfSgshF01KyTW\nl7ju+v0E/LIOah2WDcjJwLljuYZuv364TO5bXmEcrFpFvGwHr8Z8hDw6gbf9XmlQH3UKgiAIVoIg\n7AfOAdtEUUwEXEVRrIj0cA6o2Ax3B05VKX4KwxWDa8VUVkryxASGJ8ZrDrmsS0tJGDQYuTIB/4BE\nTb4JW/4mz0ltJBJ8cJ8m/aXSJeqn2VU53xc9ybF/e1Vroy6i7I2rVJbWux+T/l1/R0cEgMRPrRoW\nVdPFBTIzGffZeQY8Xu51O1tt2UdUFIBmhDYUfUaEsvKpURcgWBAE79vuqzfW66iipsSoTZUvpWni\nBhoFEYF1w8cSEfSXJi0o+B/mR7xRLW9AutpQ3+n6NWQ976wXvPTefQk7omRi598MGxGMyJ5PrInu\nbZg985q2D+I8XQrTp7P8fyM48ONtP1kPD0j8neXLzzIrvjVRmwzro97zd1EUrwJxwCjgnCAIbgCC\nIHQCKixMTgNdqxTrUp5Wjah7K18yT0O63jhI9u2ibckVFn1ZuWM0t/27nEnrxu8B47hp71BjOb+s\nA3d0atT/cBarAyfxUMvVd6wPFQS/WIrslGF+Vn1yMonrJWP61FXQrx+43ba7JpOB+xBO9PiJCS/2\nJ8pAbQ1du0bOFTtCgiDYAyOATGAD8Hh5tseBis3gDcBDgiDYCoLQA/AAEmnCVEyNfAakEFzuAHjM\ngRj+PDiZaUnLkbnF8pv/RB4a8LOWSrFf1sE77hvVJyeD/SUDeejhH3VnriA/X22e+a7xYhfs+cQa\nZRf9PG+n9vFm/sVQ8tpWbj17H87ki9QXiY2Q8MLv+dCp3B7FwwM6d+atJz9j8Yi5Wqoz9UXXiNAJ\nUAiCkIr6B71NFMXNwLvACEEQDgHy8mtEUcwA1gAZqH0iTBfvhA6HEUkYNITiYhscW91AdUKOKl3b\ns3PS3sEsK3qOTXsf4Forbd36v3Y8SM9eOY3dZQ0dL+cxsMV+Dh4cqH+hggKdjn0NoT4jgtUvaZo1\nWFU0C+mz5buGOervNqxMyRspOrz06aBOFQtRFA8C1bzpiqJ4CWp2iyaK4ttA/ZVTzARlwDBKRGuG\nJ6u1OF//ZgmDI5K0djN2DQhAGRhC3nvqbb145XAGBSTS5kLDQiEdzunDKVd3upxruF1CBdOSv2Fy\n23W88MhHrPr1ibozb94MiXd2APfNTiNvZSietahaq3a4AeXqNJ6e9Op9yCjfV5M5WW4s4oJCUfgN\nAyB+0BCiA4ZX03+P2LeN199bAqjtlgH2JQVVGxEA1rpNoH37PL3a7u2RbVQhuNCuA1LvGE67uLF0\n90v4DNBhMG/CwXv+EN3qEE8MWIHUOwZBFLWmRlr07Qt9+qinR4cOcUTlyCnXhlv2WQThNmR7dxB2\nYBdxgSEogkKJ2BJDfh369Unl3qNrWw8Imn9087DLr8wZ+BYn3Yyz49zxch5ubmcJy1Yxoeuauvxz\nqamvIY4RSenrw1//u078/64jj05AZh+L1DtG8zrXXj36RvaNhGvX6OSufmD0sldYRoT6UmRryy0b\nG+ICQ0gYNJhSa2uGe6vtDJQB6lFAEEV2BAQjj05gcHIyLVqUsMNvsNYf5cptzn0lUgU3brQmckj1\nLVXf7DSmz/hIr/7954fV/Jr5GD93erjhHxb1iJCb64639wEQ9fhT55surtwXvW34w6t2jRtREMB3\nIKSmwqVLpB0cSLxyuOY1zWcpL/h+jOqiDKRSzl7vA1FR9H6krVGcHTQbNezaKLCzZ7ePP1EFkdwq\nskMQRHZFh9KiRQlDQ5QaewK5EE6oNJaEreGIWwS2RISzd6AfJdEtGB+7kfz8mrdJV42exENJa+lk\ncxoRddyBqsExljz6nMYHki56e2Tj7n6a7/59mkfdVtKtgX/gjpfzWOo8kwE5Gazwf4JnUqvr/mvx\nxRcNak8XD8qHsDpW4N5jO3AsuKl1b1DmAZzaXi73LAtrpo5j2rYV5OWpR4KszP6Ey2OY8fk30Bc8\n7odOnRTISlRG6VuzFoTfA8bxweVX2Btd3cqqpKSFllGNKAqa620RYYxIVGJXWkiLFiW1CgGole+m\nDPgF92unyD3rztV2bbQE4fO/X9K4j9eFu/tpxthuJsgzkW6JxnHpMiAng7y27Vl6dqZR6msoU8IH\nEzA9iJfG78XDOpeAjEpv4aTup2KSsk4Yj8JTxrJBT/PVFy/St186XkeycXU7y7ncTjx9z9e8ss84\ntgjQTAXhhoMjm71GMmXLunqVk8hiCZUoGPGlElCHjmrRoljLI/btPLRlrXbCbbt+ku5Kjh7x0Kv9\n963m8Ic4kQX58+nt9kSDR4QK1vZ8gIPb/GrPcONGpWvHRiDpSyse/jKY1u4iP8d00jw4rn5Uvoi5\ndIkF2YvwPH6EZTmzUHkNJ+9iRx7K+p0Hxv3Oiq9fIKxIadQ+Ncs1wksR79RbCADileHIknZouSF/\n88Xq8/7a+Gn01GrzVdk9+h+qRdnPR0AkXjmcE5266F2uLuICQ3SPBh9+qHb12MhcPyMwrn8Q8uBw\n5MHhXPm3XBB++YVI5+e52FY9kko7xeHiepa+/dJY8fULuLkZb2etgmZjs1yVVQGTeHjL7waVlci2\na6ZIW8aEs6RwlsbJVwVBwf+QuEc93fIP2IOj400td+23838DvuI7xbN696FnrxziHOR0O3tKd+Y6\n+HLQ0zy/bXndmf7+G5LVbvOdvUS8Hiyj4BLsXWpch76G0P+F9ix99Hf6H8niQjtnlrV6ll7WRxhz\neSt9/639oPKucfmoi7Bjhrs0qbpuUAyUIAnWfqKvGzMW5alwxoWpBU0QRDb9e3+dvna+PfAcHp5Z\nOtvu7ZGNRKrgqZ7fNEgIlg56Dql3jG4hqIL/c2XMTziLquNiVAM+ZK3ScIMcY5G+9BLy4HDCrqxG\n4STjzVPvMGvf53UKgaGYtSDkOruQ0bNPvcvNl9Xu01QXMRFS5s9/FSurUt5/J5Ktu7X9c8qz4rEv\nKuTPg5OJjZDwXsc5GnuEqqT09eFym0ojf2kX3WoLvT2yUaWHN0hdINWzPzO3LavmXaM2Wt3IZrLK\nA2Xgx7ywSu2b1b6oEPlx45wwVwQXbAgZ0+OZ+cZDyKwUJrPzNmtB+LIomNVzL5Heq6/eZVb6P8iK\nr14wuM1FN+ezePHblJWppwa3/6BGuW1C6h3D2oD7NXYLt7Ov3wDGn/qLmSGV4VVXHJjO0/K6o25G\nFSwEIKdbT4P7r6+j3A9ensFPo6fS+uWprD74KK0KqguzPoTOr/uHHr/ASFOsDh3IyuzPb90eNE59\nt2G2u0ar3EMQfj4EIug727vm2JqP8mY3qN3bfZJWIJFtp6jQnj271IvoeIZzwcMZ58vV1Sf29fPl\neEIPnJz8uOTUlvZX1bvj+gb39jhx1MDe60e37sfwz0jl4zQpkj5xcJuzi4KWdoxO0e9hkrCwkdYS\nGRkseekrHjrxh0mqN9sRIfOKE8Kf2bi80Afvw7rn1wBtbl7nA5fZDPJv+LA+LFRJ5Ly5mut45XD2\n7B6mlWf+yJp3lAZlptK27WUOpA7iUpt2HPDszzO+X7A89vkG96suUvr66BUb4cSBdsinvULKCmt+\nO1jdgXKBnT17PjHeD9yure48OpHJCLuhxO3iOd15DcAsBeFMRzeEvw4h/seb54t36C5QBXliAptu\njm1wHwRBZPE7bzNkWDyxERJsbas7ukraO5jVAROqpVdGSQCli4SJF9ayInZ6ne1F3Lu+QaFXz3Z0\nZerF33RnjIqCJUtg3z56jSljarG2N4h8O3vGJBlXYAuv6M5TJ15e9JVewDXvglH6UxNmKQgCIqKf\nG1MeKjaovNvF8zw7/dMG9WFHfBilpdYMC1EhT0wgVirBtmWRVp69iUO0bGUrDHP8M1P5eJJ6avH0\n1h84nKN7wW+oF4vTLp14aeB7yMriyM7yqj3jpUtwWdurRnyUNRs/cGLWxNeZXDSLuMAQZKkvk/ip\n8UaD4FmlLC6s9ELo0LGOzLURFsbqtpNxyzuvO6+BmKUgdLpwDr/JbfA6mm1wHQu3L9adSQdDQ1SE\nHVAvhkNS9qAIlSKRKrCzqxIbofzZv2zykxonwGCYqaaqRf1cur/n9xJhJQqWbJ/Doex+tWfMyYHP\nPoNPqz8cbpwR+KSTDb+/3AZ5cLjRzw+2+X6K3a0iAl8oJTZRweg5BthsKJWoWhvR3X0NmJ0g/OoW\nwvxukxivMNAK24jMb72QiL0xmuth+/egSg/nrWmVSnRfFk7jN/+JTFvzPVui7+dnf7U/UxFB/1V+\nOXE2Mr3znnJ155tjz5BzqJYdtQsX4GJ5hM2DRnD9aCBxPYIYkrqXuRMTkScmIHvKAEFIT2fmx8+T\n6tnf+B0sx+x2jTp1E/C/cGfdj1Qg2697fRKvHM7NgNbc67CV+RHzWJL4MmX+ViiK5dXUtfUhp1tP\nvXaNfnJ/hMMxtUy51qxRe4a7g8w6V8wXvWxQ9Azm078qDRaD0vbh6NKJm+fr6X4lI4PIc+NYb6Ko\nmmY3IsgTE4xycuhYkK+xHjMEiVRBi7Ka98iHpO7VxFADSE4K5r6ef+Jx4ihHDnsQ5y7hJ8X/1bvN\nwzl9arfM0sXZs2rFuaioOy4EAEtcbbh1A8Zc3KPl8jH4YDJj5lwnNlFBbKKCqUv09AaYkcFfj7Ug\npa8RPHzXgNkJgrFwLMgnMKj+ntokUoXG7WBtMYOHHEhi1GDtqVu8cjizz31EmzZX+fH7Z+vl7rEC\nfRbVFYSVKasnmsDovqGM8QthzvCXWNe/Mo7173ZLNIeRK20+0whFbKJCcxLdf2oZkshSnPs1ji5c\ns1S6q+Biuw50zKlfJPrSjlZY6fmdjA7cwNZo7a3agX5J7E/Rz3VJTYwavZEtSbq3f0f5b2Tb1nIn\nPmfPwtdf113gDhP0v1I6d8knL9+BzV2WkNHTk77HcmqMs6wICiXnfHuePfYXz1o/x/JpakeKfcaX\noZyyXOdZgkXp7jYcCgrqPSrMn/a6Xvk+efS5akIA6tgITk6Gb5wfOeKp02b5x7FTKoXgzBmzF4IK\n/nylNUmfW6MICuH+l8Yi2zWLR0pmVgsmKE9M4Nljas+CVe2s+04ou7sO1IyFQ2E+yjPhBAUbX5My\nPlUdXKSmdUjViDr15XBO3YJwvFMXFsYvqkxQKg1uqzFJ/NSagBnqac/9PoNxDxTZt9yKlS+24y+3\n2t1BPt9DxbNp6vOkBa41BxUxBs1aEAAcCgtQngonaPDOaveGhhhmlbVrQABbdqttFCq8WBiTKKvI\nWu91P3tK22mYpxn7y6xC4AulqAZ+hOKPjbj4qIWggo9XD+JKa6cay/nkZCJrmQnAQ19NuDu1T42F\nAATXIAj/7JBWS1u4TLdvsiEHkrQCjlcgkSro0uWEQX2sSsy2CHb4BeuX2b3hPn1MzQv/FqMa8CEO\nhQXV1Ej8ni4j6pVk2l6/WmPZtN79mDlOvVOUtc7K6FFIK7grBOGVSQtZ+skcnfnmRb6h90L5zf9o\n2zwM8t/L38fG1ZK7/kRR86iQMGgw27dF1HjPXFnaw0bLG7bk8eu0clN/zykrrBjvFcgx967Vyh30\n8CJs1oNcSFcvFPqML6PTRdOoWTR7Qchzasfnn9YdRMLKqgypdDuLvtDfU+XgA0k4OKjPEiRSBcqL\nctrcvI7EwzhbmLExYzR+lI65dzNKnXcSscqq93e7JSh++5tW7pUPnTgHX/b39dGKfvP5UQkXMyvL\nqaswzS6n2Z0sG5N8O3vGOG+EI9rpC954hZCUPZprq7IyZEnVp051MfhAErHBEvLtHbSMc6SeClbG\nPX/ysnIAABIISURBVF5HyTrIyQEHB+isXixXGAXJ7onjPwN/YOH+RUSV3TZSNIHF8qbUnUQkbNdK\nC0xLobX7WG6U2+H/edaXlxYOouswkU9fgKhoGfFR2npPzz6Widtp04wIZiMIVZ8EnS6co++xhp8u\nz37gLdofzdMytgewKSmp08ZYXwY3QG26GsuXq7dCAUaNgiFDNLeOH+vJT+KTlA204t+jvbXLHTLj\nKCvlbGkThH1gIcEHk3EorFxbSR+5ym9J6h22v59U/xSTPvmZnW5BWkIwcU0JEwdlMtUAzyT6YhZT\no5NunVFuBlXwDlTBO1g1/7omlGxknvr/tN76m2tWEHYknnil9o7RGNkGXvvWeI6hjMKGDZVCAOpo\nlj//DD/9pEk6frwHb22Pqi4IzzzTSJ00nKU9bJAHh/NK6CxN2pdTnuK3WdW3mff36c+332v7gVo7\nuQWLXu3PuQ6mi0ttFoLQNfc0CzruoCzxCQCEfblYRSWgCt6BEKn+/4+XLnKwdx2qxjUQ10tCQYGD\nZo1g27KI121M67H+mbU/0bXbcd0Z166F1avhq68q/f1X5cgROKqHyaa7e5MQhttZNLJmv00/7/Wt\nNiUCtaNuU+pAmIUgVOBx7AhiZC2hQjMusvbF8xz0qMP45DYWbntL67qlbZHW2sBUSHvr4dTr/HnI\nzITc3JoFoYKoKPjhB7VCXWysOu12r3Tu7mp36U2Ezx/6L7kpNWtA3C4EXpPLeO6rXJSTvjKpYY7Z\nrBEAHjufAB1BkajealNtEuHgBYR15QY6WXmkesrxydFPuzKuVyhUWWpI/BWQZuxeG0BGBpyrh6rA\nsWPqF0BC+domJwf++1/1+1WrKqNMmjF7P7dC6jyP5P61P38vH9EWkOmPpvH8iT+h3EdCrrMLbibY\nQjUrQahAnpjAWWdXVKlOCOurLAb7dmCAnkIAcG/aNiRSBQMGJuPmkMsbKz42QW+1Sfby5USZju1O\nLy91yNTzDfiDnjqlHi08PZvEghnq7z1v8voSJl+PJ71XX5adlHDwjAvnYgWeeTQLvy7qUdTjxFGj\nxEcwa+3T+a3HYvX431pp7b8MZWZZ/XZ8ClraacU9MzWPDvieXxVP6M64bFn9Roa7EK8pZVxIEzSH\narfTO6KMuCe/1RKGZqd92q91DVqcBghuYwoBwKziT3DSxx7BIgQ6yVhtVasQABzebIXs2/9rcJQh\nsxSEsx1diQsMYeqZBCSJofRdX3k83z6i0x3smZG5PWawBYM4ssWKn0saZtxvVoLwi2sI8y+Gsuya\nP6rgHUTmhZK5H3xaXqDs9/GI4zx5dJNhXq4bE//MVBTeMt1Warm5jdOhu4A3HvJkcQfDQ26ZlSAA\nCElnsXp6s/p9ZAIXnklgzQ8CM4oT6POxaexVTcGgrAMovKV1Bx2fMsUyKhiR777rxfFO1ZX39MEs\nFssn3TqjLOvBkcmVXiOqnieI/TtCxkWsfstg4qeuem+f3klS+vrgl3WQ1D7eyLOUXMqrJXxUVpZa\nm2zVqsbtYDNFtSsG6ZB/6r1Y1mv7VBAEayAJOCWK4lhBENoDq4HuwDFgsiiKV8rzvgY8BZQCM0VR\n3FZX3SfcOvPtipYIW9VCIPq74TatN9MKK4XiYJt+rI1SG2esnSkiLu3PgEOmcethDD4uC+Xafw/g\n9Mk96tPQup41TeggzNyRRJbS/exJg8rqOzX6H5BB5Z/0VSBGFEVPILb8GkEQvIApgBcwGvhSEIRa\n2zjRqQvfLrdF2KpWJRADO/HwglZaQgDgcziTst/HA1A2tT9xh+rvL6gxkfa/hmTNAHqeOsbA7DTG\nDfyDDh10OBF42DghZe9WFv92CFXHxXQ/YyJBEAShCxABfENl6Oz7gR/L3/8IVFikPACsEkWxWBTF\nY8BhIKi2ur851Bth27+aa/dpvehz7HC1fOfbOyP8rvaIbRWVgPhIfz6hdjvXO41/ZqqWduu3B57j\nk/vqdgKMpycE1fpVWaiFXmPKiE1U8MbFhk0t9ZkaLQFeAaq6GnAVRbFiE/wcUGE/5w5UtWY/BdS4\nwasKGIbVYKXmWgxyJ6zjWaiir5br7EJmD0/OnwVhTeUocaWd2mvDkqWhzBL+v71zj67pyuP453eJ\nekZEUrS0MSWkEu80HnnchIrX8uiU0a5OtZ0aM2Z1lCqqnSZRnUV1Jp3p6qhB1Wi1oymWVjNE5IF6\nEyLeSqsUISiCkuz545zIw725j5wbkpzPWnflOGff394/9/7uPmc/ft+KL6e+ZzhpjIpmTSEyroBP\ngz+i5daKzyyXGwgiMhg4q5TaJSJWW2WUUkrKV8Cwee3hUyduPxC3bA1RTX8i8PvSO2j+dSMM+d1m\nZG/pdOCFq0ZiGbSUiytOwPDyPKhiREXBkiV3uxX3PIFDCpnzRrrW6xo0J+moR+gFDBGRgUBdwFtE\nFgNnRKS5Uuq0iLQAihbNnARKjl+11M/dwcfzfkDQNrq3vWA7GcO4+7ZwfnZT/nupN5ZRWioPu6tT\nPURaaDjx1+J4te67NJCrRG9zTa8BYE9gB6Zvmu64YGAgWK1VYtfZ3aLdsELWPTWPB7YWz8GkH4L0\nCu7jcnr4VESigEn6qNE7wHml1CwRmQr4KKWm6g/LS9CeCx4E1gJtVJlKXMl0tyewA8v6aSNETfJ+\nz/gl/yYpIIK947MZOrcVXQ4Ym+k5vXtvlAjx1+KIDk8hYcasUtcjrWt5ue4/8JGLdoNiZ1BHLjX0\nJv5aHPH1Enhy6wr7w6e2iI+vgAfVk/a/LuSDKek8+t0hj2S6czUQXlFKDdGHT5cCD3Hn8Ok0tOHT\nW8B4pdRqG7acCoT9vwpk941m7M9S+DWHlwqKv3hZ7YLpfNC4NdVpoeEk/M9KpkpAlTOG8FiPjezJ\n6sr16/WItK4lvl4CwO2gSAsNZ/i3X91O8hXSMYsXYuYyY/FbnD/v51xj5s/XVpea3Gb4klssy3vb\ncUE8HAhG4igQlrQIhy7NOBh/kITHPLuBoMBiIa1LOP16RRPxZgFxg2wn/Xo3fxLJ6UMAmDTlLbZu\n7k1mRszt65FWbXN6SRXOkE67SL/Sh8kBf2VB6h9ca5jZK5Ri9+ovnF6C704g3HP7EU7e34KDH5+B\nGYdhWCBx5yNQCkb2yCX4yAHO+fjy/lFNMMLHF6xBF92+PSqoVYuYM69Ra00hqZtTtV91OzqEMaxn\nsHU5q9KHkZkRw7CoJLZv60F+fn3gThlagBXyBJMCZrIw1Y2tlK1awQn3xsSrIxnX2tOs6TmaeWiX\n2j3ZI3zaPILDT5YZFn3Uj8KRQYQ81YS9S/KQBP02KaAx3hM78qJ3Ft5XXFdjyejemyg7qVwyu/bk\nL/lvEt9I65KHbljF5cvFo8iR1rV2hb0XxD6LtLxF4rbJnPyxFXl5LjwjHD2qbd43KUX7JwpJGznX\n4Q61anNrBFpq8IxVettyziFJxRKzTT6IICQULpyD7IFawKgXO9Eo2JtJtY2dV6h3JZ/r1+u59J6F\nsc/w3I5PAchr3IQwyyaXtA9YvRo2ua7tUBNIXJbDy6fK11quVoFQkrO+/qXSuTS9mEcnfa3Rskci\n2L30ErJwDy1Oj2Zs0iJ7Ztwis2tPpqu/kJoywKnykVHrWPnDcBpfLlaS3B3Y4bYSzm+2JXEu10Fa\nkiNH4JNP3G5zdScyvoDPgubzQK7tZezVNhAccamRNzuCOhmStKssmV170jctg5s3vZwqvyj2aZ7d\nYX+6/4+d/smHqQ5U7c0ewSGzk/Yz6fRSm9eq3VZNZ2l8+WePBAHAhq49uHnTy6bgeFmWjxhUbhAA\nnG3ih7+/g+nQ2FhtYs3ELunn2zou5ALVIhA8SfjOzXh53SQiKtVh2cTcCUQFp7Co29OkhYaXSmMJ\nkN0miHPn/MnNdSK1udUK4ZU7i16VWDW2NluDuxpmzwyEMqSFhpPRvdft4157tpEwcTL5VxsyeWpC\nKcHx7qGb8b+/+Nc9M70vmYlNea5nW2KmziAmeT3zu40mKjiFtNBw5nQbg7+/55JU1TTi9zn33OYM\n99w8gqf5xcuLjZ1ti3AkXIsjI7kvFksh4ZHryEzuy9TX32T6gpmsbx/Ne4lTEVHcuF6XyKh19O/x\nFdNmJdqu6MIF2L2bxSGD8b2ezojpMcSs8Gfc1Q/Z3yGYfTkdHTe2b19tLuF7J1JI1lCc1aV2RI0K\nhJUhfZj5dU82jbWXaGoDxPWhV+8MUBaSB/QlceMEVneJ4ZttQ9kYEcZZL3/ey5/AZN+ZDFqQQn5c\nPWYk2MinmpUFQOaYU5BXG0JCOJNzi5zo9uz7xokgKOL557VFeOZCvDtIHlebPWHPEJlQwKwrH9Pq\ntPvL2GtMICT37sPQLo7vuSMi15F+oB8WVQjAfaE3iN62gV+8vBCl8OEisT2/ZvDMNdqEWkKZCbWS\nWsd6MBRxPimHJL+RrjfeatV6hWPHHBataZzcInzWvzbBn0cyDfc351S7QEgLDadh/lVCc3bdPnfB\n24eBTgQBQGRUKpb9hWR268WM1ChSwuoQGW/l5lXY+WVvbnS1QnJdSE4m851sCPPS9hEA7NwJP/9s\n23B2NjnZFmhveya6XI4fN4PAAa+PCiTi2x6sPhoI/Mdh+bJUuXmEcz6++F3Mu+P84rYD8X5Q8fcv\nurFjroUBr15mXPg28vLr4Vv/GsubR/L+w87NBUTGF7Dh7VoU3ixz4ZVXoFGj4n/n52uZ92bPhm7d\nYIeLwiFBQcXpXA4e1DQSrNbiodPvvtOS/ipVnATYxC6tHy/kWIoFSKj+E2qjHnqdQR0P89uviydT\nLjX0JnrDeHbNt9CwhaLr2ELyc2H7B8XPAmETCtiS6FoSWqdo0EDba3z5MmzfXnF70dHQsiXUqQMp\nKfBDxVU6ax7VPBDmPDSMnSdaME/NYX7rITzip/UME+dGk/VRNRoJbtMGbtwwV5+6jeuBUKWeEdIK\nOjBk4CFYBWm1QxgTVqWa7zxH7szkYeJZqszPaM4j7fliuPbFT3p8CEtiHQXBcY+3yayj+tRRJQLh\nrK8fIxeOoGl7xbxF7RjRrosT7zru6WaZdVSjOqpEIJzx9WffUgvNOyubQnMmJhWlSgRCs7xcHh1Z\nSM7nVaK5JlWQuzZqVOmVmtQoqsTwqYnJvYZ5r2FighkIJibAXQgEEekvIgdE5LCITDHQ7nER2SMi\nu0Rkq37OV0RSROSQiKwRER8XbX4kImdEJLvEObs2ReQ13a8DItLPTfvxIvKj7scuERlQ4ppL9vX3\ntBKRNBHJEZG9IvJnD/hhrw7DfBGRuiKyRUSy9DriDfVDKVVpL6AWmmZCAOAFZAFBBtk+BviWOfcO\nMFk/ngLMdNFmBNAFyHZkE00cJUv3K0D30+KG/Thgoo2yLtvX39cc6KwfNwQOAkEG+2GvDqN9qa//\nrY0mPxBmlB+V3SM8BhxRSh1XSt0EPkcTFzGKsiMF9gRNnEIptR4oK41piEhKOfbhTj/csq/XcVop\nlaUfXwH2oyVoNtIPe3UY7Uu+flgH7QuujPKjsgPhQaDkSjK7QiJuoIC1IrJdRMbo5+wJmlSE8kRS\nSmburYhvL4nIbhFZUKKrr7B9EQlA64G24CE/StRRJBhjmC8iYhGRLL29a5RSW43yo7IDwZNjtb2V\nUl2AAcCfRKSUtpTS+ktD63fCpjv1zQFaA52Bn4C/GWFfRBoCX6JlKC+VG9MoP/Q6kvQ6rmCwL0qp\nQqVUZzTdjTARCS5z3W0/KjsQygqJtKJ01LqNUuon/W8usBytGzwjIs0BygiaVAR7Np0WSSkPpdRZ\npYOmW1fUnbttX0S80IJgsVJqhSf8KFHHJ0V1eMIX3e4lIA2INcqPyg6E7UBbEQkQkTpoCpwrK2pU\nROqLSCP9uAHQD8jWbY/Wi40GVti24BL2bK4ERolIHRFpDbTFbm5t++gfZhHD0fxw276ICLAA2KeU\nes8Tftirw0hfRMSv6NZKROoBj6M9ixjjhyujKEa80G5dDqI9vLxmkM3WaCMEWcDeIruAL5pqzyFg\nDZqyjyt2PwNOAb+gPds8X55NYJru1wEg1g37L6BtuN0D7NY/1Gbu2tffEw4U6v83u/RXf4P9sFXH\nACN9AUKAnbqtbOANR5+xK3WYSyxMTDBnlk1MADMQTEwAMxBMTAAzEExMADMQTEwAMxBMTAAzEExM\nADMQTEwA+D9OCUNImq7kawAAAABJRU5ErkJggg==\n",
"text": [
"<matplotlib.figure.Figure at 0x110232310>"
]
}
],
"prompt_number": 86
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"np.unique(rt.variables['loci_flag'][:])"
],
"language": "python",
"outputs": [
{
"output_type": "pyout",
"prompt_number": 85,
"text": [
"array([0, 1, 2, 3, 4], dtype=int8)"
]
}
],
"prompt_number": 85
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"imshow((loci == 1) | (loci == 4))"
],
"language": "python",
"outputs": [
{
"output_type": "pyout",
"prompt_number": 91,
"text": [
"<matplotlib.image.AxesImage at 0x112d69110>"
]
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAMIAAAEACAYAAADhtqk6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXecXFX5/99PeiM9pBdINmWTkEJ6m0EgCSoEERXQ7w9R\nKQYUEJBEyrqIFKWELxAEbHxFQAQMoAIJkZlN732zySaShPTeCWnn98fce7l35k6/d3Z297xfr33t\nzC3nnLs7nznnPOc5zyNKKTSa6k6Nim6ARpMPaCFoNGghaDSAFoJGA2ghaDSAFoJGA/gkBBEZLyJl\nIlIuIvf6UYdG4yXi9TqCiNQE1gGXANuARcC1Sqm1nlak0XiIHz3CEGCDUmqTUuoU8AYwwYd6NBrP\n8EMI7YHPbO+3Gsc0mrzFDyFonw1NpaOWD2VuAzra3nck0itYiIgWi8ZXlFKSzvV+CGExUCAiXYDt\nwHeAa2MvK/KhajshIKjrqJZ1FKd9p+dCUEqdFpHbgI+AmsAftMVIk+/40SOglPoA+MCPsjUaP6jC\nK8tddB26jpTRQtB16Dqo0kLQaFJHC0GjQQtBowG0EDQaQAtBowG0EPKKAsoB6MBWzue/AHRmc0U2\nqdrgy4KaxklfVtGc/YQJWMeu4h3OZQ/TmMBO2gBwHa/F3Ptnvp+rZlZrtBA8ZCBLCRBmKQMBCBJy\nnA8SYioT2UMrAMroYYkA4DWuI0AJ7Q0fxZlczGY6u9Y1jPk04wAfcJkPT1L90ELIkgBhBrGYRhwF\nYBkDYgRgZyJTmcpEyuhFbU46zpVTQIASAP7L+cxmVNxy5jMs+8ZrLPQcIQuGsJAgIUsEAK3Yk/S+\niUwlQJgV9HMc/xG/pz1b2UhX/sL/eN5eTXx0j5AlT3A3F/EJAAJ0Mya8yVH0Y4UlhgLKac82AMK+\nuytrotE9QhYsZAjHaEgP1hEiyPt8nRBBDtEk6b27aM2VTAOgO+utifIGuvEZHXxttyYW3SN4wJPc\n5XjfhENxr32Lq9lDK0YxG4CiqE0kG+jmfQM1SdE9Qo7pxVp2c67ruXIKWMDQHLdIA1oInjKQpVzB\newmv6c0aerOGd7gq5lw5BX41TZMEPTTykMt5P6XrruYterDOcaycAhYx2I9maVJA9wgecSFL0rq+\njF6O969xnZfN0aSJ7hE84BZ+R2t2pXx9KYWU0otdtKE1O31smSZVdI+QJYNYnJYIICIEQIsgj9A9\nQhYMZhEBwilfv5o+vM03Y46vo4eXzdJkgO4RsiBIiIYcS/n6Pqx2uGCYq8rRE2dN7tFCyIIGHE/7\nnnR6EE3u0ELIkCEszOi+3qzxuCUaL9BCSIF7+C2TeZRzOGIdS+RqnYxoMayne8ZlabxBCyEJtTjN\nb7mHOpzkCOcAMJQFfE79jMu8mrcAmMaV7KQNIe1tWuFoISRgMo9Si9P8gkcAqMsXAIznQ5qzP6uy\n+7IKgJ20YQdts2uoJmu0EBJQh5Pcy+PU5hRPcDcAk3jc83oGstTzMjXpodcR4mD2AiaDWJzVvECT\n3+gewYXRzKI2p6z3cxnhuQiu4h0AyuhJd9bTju2O89rMmlu0EFyI/tDPYSQzuNTzem7hRa7hDXqw\njht52To+lul05DMuZqbndWrc0UMjFxTO9FvHacBcRjCGEmvC7AXzGGZt14TY3Wpd2chMLvasPk18\ndI/ggj0QVyrHM8UugmjMHqiI4hiBaLxHC8GFeOPzeQxPaWO+F5ghHwFe4fqc1Fmd0UJwQRKkip7C\nHTlpQ1c2ApGQj5tykm2meqOF4EINzlqv5zAy5vxHjMtJO/7EDXFDPmq8RU+WXThLDWpwllmM5j98\nJeb8fIYxhhLq87lvbZjBpWyhk2/le0WAMCWMYRSzqcFZanOK08bHai4jOEmdCm5hamghRGEfFo1m\nFnMZwQnqxVwXJsB4PvSlDRvpylxG+FK2lwQJEyBk/c0qs8+UFoKNIKGUF7IWMJQxlGS0JyERG+jG\nX/mup2V6zUV8whhKCBPgIYoSzKgqD3qOYCNEMGYNwa03MDlOA8/bkO9xT2tyhlHMJkTQ+HtVDZIK\nQUT+KCK7RGSV7VhzEZkhIutFZLqINLWdmywi5SJSJiJj/Wq41wiKIoqtbn4WoymmKOE9Xsch+jdf\nZSvtPS3TS2pxmsk8yq94wPM1lYomlR7hT8D4qGOTgBlKqe7ATOM9IlIIfAcoNO6ZKiKVotdRiOOf\nO5pZ1ORMwnsWMsSz+itLgK+Hub+im+ALST+kSqlZwIGow1cArxivXwGuNF5PAF5XSp1SSm0CNoCH\nnxafiZ7s3c/DCa8fzCLf6s5HTlfhKWWmT9ZaKWUG89kFtDZetwPm267bCnnc19sQFA/yUFr3ZLpv\n2c46ehAmoDfnVDBZS1wppUQk0ZypUs6n4q0hZMMaerOXluzmXPbRkl1xomJrck+mQtglIm2UUjtF\npC2w2zi+Dehou66DccyFkO11F+On4lAIxRRZDm6jmcUZasadFA5hIS3Z63ouTIADNCNAmGbGqPI5\nbmMfLfxpfLVnk/GTOZkK4T3geuBx4/c02/HXROQpIkOiAog3fghmWHXuCBJiAUNdTahuG3XCBDhM\nYyurZnSONI1fdMH5RZr+pqZUzKevA3OBHiLymYjcADwGXCoi64GvGO9RSpUCbwKlwAfARKVUpRoa\nRa8jAIxgruP9UBa4ulcoRGe8qaQk7RGUUtfGOXVJnOsfgagNv5WEAOEYz9N7jc36K+jHMRoCkRzH\nbgQJsZGuHKaxvw3VeE6lsPFXJrbqRICVEi0EG27jftMN2+wNAOpxwvX+3/MjX9ql8R8tBIN4m3HM\nOYM96nU8IRzRQ6JKS9VdKkyDAOG44VrMNLAmdzAFiLhXRC+oHTZCQmoqH9VaCDU5k9SNAuBjLrGG\nRvWMKBYhguyjBZfxAQBPReVa1lQuqr0QEvEokxnBXMd2zbrGsMgMAhwmwGIGcZRG/jVU4zvVWgip\nbCO0O8MNZx4Aj3Mv4K33qaZiqfaT5VmMcbw3J8ezGB0jlHkM5wvqxp0sayov1bpHqMkZRlPiOPYQ\nDya8J9GONU3lpVr3CMnmCG4sYBgHaZr8Qk2lotr2CPU4YblPmKSy/XAew/xqkqYCqZZCuISPGcmc\nmOOVYZeYxh+q3dAoQNhVBJrqTbUSwjkcibuC7OZ+rak+VJmh0QCWsYFu/Iyn+CdfpxFHOU4DWrAP\nhXCCegmz3lS18CSa9KgSQmjCIa7gPev91/ln2mVoIVRvKqUQmrOfn/AsB2jGZjpnPaw5Q02PWqap\nrFQqIQQI04J9Vo7iZhywNsdnyllqVNmgVZrUqTRCaMUezzNbnqQOjzLZ0zI1lZNKYTVqxR4mMtXz\ncutwkgBhRjPL87I1lYtK0SOca4VN8h6zl7Fnycn3ifMP+CN/5AcV3YwqRV4LoTW7uIXf5aQuUxDJ\nnO5ySSe2sIVOdGYzXdjEFjrRju105DPG8VHOUlhVB/JaCF7PCZKRLAx8qvRnORN417XMQkoJEGY2\no1hFX8e5dmxnO+2s9zfwJ0/ao0lO3gmhHdtRCDtoS0/KKro5aeG29/kbTCPMGNqxnRbss85fxTtc\nxTsOodhFADCbUTF7pk020tXTtld3pCIC0UWCBrt/+wYJcZjGFFBeIUIIGZnBUuFWnqcle1nMINqz\njbbsSLu+MnqykzYxx7uykY58Fvc+r3qvqkkxSqm0FpfyrkfYRnvO4UjORWAK4EEeSlkIz3MrQUI0\n5BjncDSjentSFvOsW+lAB7YmvC96jtCDdayjR0Zt0OShEMopSPub9RBNaMKhtO4JE+AMNZnLSEYb\nifEg/clyiKAVQdsrPmIcW+mQsNxhzKcF+1jEYMop4BresNqT71avfCTvhHAtr/M68cKtxnKQpjzD\n7Uk/jIdoQimFDGcex2ng2HuQyT6EISwkQNjzrJoAP+QPAGyhE53YEve6AsopoNxxrBdrtRAyIO8W\n1JZwYVrXN+Wgw+HOjce5l2UMsKJQZPvhHcoCLuMDX0RgpxNb2EzntO5pzS6KKCZAmN6s8allVY+8\n6xEChDmSYsS4gzRlBf1YxgAADtOYRhzlQpY4rovekglQRLFjwvkgD7kOi+x5l4czj7pGgC+vKafA\n+nbfSFe20sHqqb7Pn+nM5rTKM61TgmI77Sz/LN1buJN3QggRZAdtWUcPerDOcW4vLR1Zap7hduv1\ne1xhvY4WghvR0Sh+xQPW8OosNRwrzX6zjh68wTVARKCv8j3H+U10SVsIgJUBKESQAzSLG99Vk4dC\nKKcAiJgV3YTwPLfSgv30ZWXcMn7HLUlXpM3YRDU4y2hmOez/uRQBRCw+9jlOgDD/5XxqcZpPOc/R\nK6VDkJA2s6ZI3gkB4HLed805bPYG+2iecIK7i9YsZhCDWJywnkytPfZca34QNPLaZ4MWQHrkpRDe\n53Ig0iuYNvYwgZStO4046tisc4J6PM691OULvqBuVh/ixyO51fOa1fSp6CZUOvLOamQnTJC19AJg\nD63SuteeDNyMVfoFdQH4Lfdk3KYAYV97A03FkNdC2Elr3uTbvMXVrKF3RmV4bSUZZphg85k+rKaQ\nUl/KznS+ku/ktRBMMhUBYJkN7URnyayKfIu/cyvP04fVnpRXi9OMZA5n43xk7uPXGYXQzBfyzunO\npD/LWU7/rOppwHH6sooFDAUicY1+xlOOa8weo6p+09knzRewkpVc4Es9tTjN6agpZx1OphR633uq\ngNOdyU7a0JMyyuiZcRnHaWCJACI9QZgAp6jtSP5hUtXEsNeYVzVnPwHCXMBKvsE/2Exn/sz3Uy6n\nIce4lef5DT+Pe81oZvEJFzGJxxyLjvMYzkKG5H3g5LwVQht2MoF3KaaIQkoppTDrMhPt6KpqIgBo\nyR7XiX1nNnMHUxAUYQIsZaB1rohiPqc+YQKOLxE79/Frfs19jmNnqEkRxYQJMIdRnDI+WsOZx1AW\n5P1uurwVwufU50Vusf6RftvFSxjDmKhcCZWZPbRiKhMJEqI5+2PmSqa37uW8z+W87zhXn88Zz4cI\nitX04S6e5HPquw5/TMxdd9Em7nkM9+6hfCTpZFlEOorIJyKyRkRWi8hPjePNRWSGiKwXkeki0tR2\nz2QRKReRMhEZm0nDTN9603z6Lf6eSTEp8wkXVUk/nFX0dTUYpMI4PuIunmQuI6nNaQKE+Qr/YRaj\nAQgapuQiivk19zGTi71sek5JxWp0CrhTKdUbGAbcKiK9gEnADKVUd2Cm8R4RKQS+AxQC44GpIpKR\ndWonrVlNX/bQikJKffemrEph4Vuxh6t5m76sitkbnS4jmEMtTjGK2ZyhJiWMoQZnCRDiIUMK8axJ\nlYWkrVdK7VRKLTdeHwXWAu2BK4BXjMteAa40Xk8AXldKnVJKbQI2QOZZ90rpxVQmMpWJtGKPZ+bA\n6kAZPQkRzFoIgDUkChKiiGLGUEIxRa5ufPfxiCNBe2UgrTmCiHQBBgALgNZKqV3GqV1Aa+N1O2C+\n7batRISTFXtoRYggt/Gcry4ECqkyXppePkUtTgPwH75iDY3cmMRj1OKUlZe6spCyEESkEfA2cLtS\n6ojIl2ZapZSKrA3ExeVcyPa6i/Hjjt1q5Pc4vqqIALC8rbzKAjqHkQlFECBcQYtqm4yfzElJCCJS\nm4gI/qKUmmYc3iUibZRSO0WkLVjh6LYBHW23dzCORRFMqYGt2e0QghfdfCKqSo9QSiGr6U0L9sVs\n58yUU9Tmfh5mFqMRFAFKCBlfTHZv2UTrDf7QBecXafqm8KQryxL56n8F2KeUutN2/DfGscdFZBLQ\nVCk1yZgsv0ZkXtAe+BjopmwVpbKybNKXVYxiNuey23cTaoAwAcJVQggAa+hD7wqYU/2Gn/M59XNe\n75f4s7I8EvgesFJElhnHJgOPAW+KyA+J9EvfBlBKlYrIm0ApcBqYqLLw49hNa8IEfHMiMxFUziPr\n+U1FiKCyklQISqnZxLcuXRLnnkeAR7Jo15dlGb/N9QQ/cItQp8meAOFKszaT18bf1uzim7wNZOeB\nmozqIIL9NE/pus+pTzFFcd0rUuE0tSmi2JOddrkir4UAkc05XvgZxSPX+5MriubsT+m6+nxOgDAf\nMp7HMtiN9xiTuNPm4VtZXCzy1g07F9TkDPfzcEU3Iy8xt7fW4SSTeTSrsnK/fzr9yXLe9wh+Upk3\nkvhNPU4wjPmcpA7FFBlepbGu63ZmM8raDhtdVr6Tt96nuSDbb7qqzjg+QlDMY7jlh/Uxl1CDs4yh\nxOG6bgZXsO8VNwkS4kPG56rZGVGthVAdeZcJdGYz/Vme0vVjmc5FhJjDCCDygT9LDWMaHIy53q2X\nzXcRQDUeGl3MzIpuQs45QDMChNlCZ6ZwB4dpnNJ9tTlpWYDMuKoP8lDMZqYAYcsnyc44PvKk/X5S\nLSfLDTjOPfy2wuqvSMIErP3gjTnMAJYlvykFQgQTmkpzG65eT5Y1SejHCkIE2UhXz0QA8ddiiini\nBPUIEqIF+zyrz2uqXY9wKTOqRTiXVFjKQAay1Nc6QgQZxnzLchQiyGIG+eymrXuEpGgRRHiXCTkZ\nqgQJWSL4nPoECXE3T/heb7pUOyE8yuSKbkJeMJClvvcGdn7Dzy337Ce4O2f1pkq1E8JJ6jDXMAVW\nZxpzOKchbPqyigBhiiliEIvzLnyOXkeopqSbfDFbLuMDwDmpzifP1GrXIwDM4NKKbkKV4lc8kHYU\ni8UM8qk1mVEthQBfhorXZM8YSpjNKP7DVxJet5AhhAnwb75KT8pox/YctTA51XZo5JZgsLqxjAGe\nrCUECPNr7uM0tVAI9TiBoDhFbWsusJAhfMBl1j0K4UZe5mnu5Hz+y6ecxyGaZN2WTKl26wh2LmYm\no5hd0c2o9JQwhk+4yPXcnTzNWnq5+hv9jKc4hyMA/Jnvp51KNz56HSFl3EQQtrbvR37OUDMm2O18\nhuXVJC8bzLS82RAmEFcEEPl7xXO6M0UAeCiCzKiWQ6PanOIMNR3HFBLjTWm+P04DR3Jx0/MySCjv\nzICp8h5XsIwBKCSj9QSFUELACucSj3g71G7mRev1FO4AYnNf55JqKQRwugufpQa/4oGU7x3HR3zE\nOEIEOUZDvsq//WiibxyiidUb2JMupkKYADU4yxlqZtUzNuKo9foOpgCwgn4Zl5ctVVYIg1nkmqIW\nIoGqwgSsoVENzvIgD/EQDyYscz7DrB1YRRRbPYWZ1/kI5zi6+/wl8uFvzzYrOXuYAA047rqxxmQm\nFzObUVnXfjdPUEZPjtLIinTxKedV6PCoygohnghMouP8C8rKxRAiyCxGu9rG7cOlFuxjCAs5RkP+\nxddZR/eY1FT5ScRAso32FFNkDe/i7TAz8UIEg1nEYgZZf8d8iUBe5YTQiKP0YB1LuJBbeZ7nudX1\nOrd9tNG5nEcxm9mMsuYHw5hPiCBfUJehLLAiQwQI05N1fI1/WvcepClNOejhk/nHEi7kLp7MSeiV\nZF9QFUWVsxoJ0JYdBAjTkr0MYWHSe85Sg4d4MGbiO5tR3MevrffFFFkj6ubsZxGDOUZDBrGYRsaQ\naBkD+BM3WDnD7Mn7wgTYQytKKbQsUwDr6Z7FE6ePXewBwhylUV46wuWSKieEIzSyLDp7aRk3aPC9\nPO6Y8CkkxmIxitmOrYdjmc4JY47wLD9hD62sPABLuBCAASxjC52sey5gJWECFFNkfQD/zrcsy9OL\n3MzrXGtFitiWQgR9u4gyYQLvcgdTCBDmHI4QIMwxGvJEkkTsI5mTcZ35TpUbGgHWN9xEpiYMRhsm\nwCxGW1EZShjjOD+bUY69zcOZx3DmUUwR+2lONzZY58xJp7kPeDOd6cImVnKB4xs4eky8kzaOc9vo\nwBjCdGCra3vt95dTQAHlGZlwm3DIMRRKZVhUEcGRe7OGq3kL8Dc+UpXrEUyO0ZBFDGYoC1zPP869\nlukw0X5bt1g+4/jIWnZbzCB2WTlSIu7NdzDF8YG1R5lL9oErpxthlwnkK1wfI6JttHcMq/7Gd3xb\n7DtNLU8my+liigAiIUD9okr2CCZuFok6nGQQi6nP51Yy7LPUiPsBupAlzGWEY2fbMObzLD8hSIgQ\nQQTl8LH/L+dbr//BNxzlmZkuFRK3zgBh/sv5fEZHurCJTXRhk0siFTOE/cvcxOW8Rxk9qcPJpH+X\naI7SyGHXtxMmYDxLeusNXvE2V7Pax7i3JlVaCG7U5hSXMsMKaWji9qEcxWye5C5OU4sZXOrwjfkJ\nz1rXHTXmJQHC1tCoS4IMLolMhkUUM52x1orseXzKp5xn/bZjb3O2ZsjIUPIFFjKYE9RzBAGey8iU\njA5+YIqgF2vZT3NH7+slVXZoFI9T1AZwjb8TzTyGc5pajGAuAcIx37am1QjgLp5iKQN5mjsZwDJL\nCN9gmjUcasneuHUFCHMdrxEm4HBLMD/80SKIxkzHmwmNOMrdPEEDjjn2GJucpDazk4R7TIe7eSKt\neU0hpQQJ+SYCqIZCMLdq1uJ0Us/TM9RkJHO4lBkECVGXLxznB7OIgJFrOGI+/XL4sJz+HKAZB2hq\nWYL20tI6f66VaQtu4iWO0ZDXuC4vFpiChLidZ3wpeyJTaWgILtVV+N6s4Vx2cwErfWkTVEMhwJc7\n1C5mZsJvplHM5hI+BiJDhyOcwxTuoJgiyxoVJEQxReygLQMNy9EWOnGIJpbFaD8tYspuwT4rWXdb\ndvA1/uX1Y+Y9qfYK9mxJPSmjDTs9b0ulFEIB5fyUZ+nHCibwbsblFFPEbEYTIMwYShznRjPLYTo9\nRW3mMdzaPBJtlj2HozzNz4CIaXI5/a1v930uSTqCUR+Ccgoyfg47Xg0fmnKQYY4swd4wlYnWcPJC\nltDeLc9kHFZyAWX0dJicvaJSTpbLKeA6XuNKprlaU1JhJhdTkzMohNmMinHLHs0sx/uT1LHG7qZP\nUpgAg1hMEcU8xV0coREQEcLtPBNj9+7NGnqxlt6sYS29eIFbMmp7Irw0MQYJMZ9hnpVnEiZAQ44R\nIMxAlqa0iPg23/S8HXYqpRAg8m3ejAMcoFlG99tt4kUU8wkXxSyomUSHf1nAUGuzSYigNUcwhbCF\nTnRiC0UUs5J+/IMrAadN3C9zpJcTyrp8wQjmeh7+xvQ3MmMrlTCmQrdpQiUdGpmYIsgmNImZOsq0\n7DzIQ9TkDDVtKaUUwqXMsCbG0TuudtGam3jJ6ub/xA2Wy8QFrLAW3+z08ilL6C5ae7oCO5x5npUV\njWn+Nd09KtJJsVrvWTYxhzqJsPskQSTXmBnC0M7T3OkIt/5T/pdmHIgpr5gi2rALUL6MeSHSA/Vm\nTVZlnKSOr9EB3f729sjZ9vMhgtZkOUSQcgrYTjuXUv3Js1ytUAgKYQ4jGcFcanKGR/iFtf5wB1Mc\nPdBiBvEvvgZENqPfydMsYwDvcYVxhfv/oy+r4joEekEhpVmLALBW3/3iMI1pzGHHMTNydlt2WMdK\nKeQAzSx/pyAhjtEwjhDSRwuBSICqB/gVEHEsM3eqDWUBjzHJsYnHFIEZo+c4DQAcbggDWMYAlnGA\nZq69QS4opZA19M5aDH5nxQwT4HLedz1nLkCGCdCLtQgKhbCOHrzBNZ62o1IJwa8E1vb9CPY9vNFp\nkOy7z0IELRPqXTzJS9zsGPM24VBCEVzFOxSylp20Zj3d2UFbLx8JgLe4mr20zOsAA0sZiEK4gvcc\nx/uyCoA19CFEkN2cyxhKaM0uFjLE83YknCyLSD0RWSAiy0VktYj80jjeXERmiMh6EZkuIk1t90wW\nkXIRKRORsV42NkzASkDuNWYUizABLuITHuBhhxCGM891JfTHvMCL3GxZjCAyWU5FsD1ZS5AQN/ES\ng1jszYNEkU1yjumMzUnA5MhQcoK1mcncb7GPFqwxcmyXUsjvuIUwAa7gPcewyQsS9ghKqRMicpFS\n6riI1AJmi8gHwDeBGUqp34jIvcAkYJKIFALfAQqB9sDHItJdKeVZVu8QQV/H1/bFtYe531pfsE+K\niymiPp8DEVPgUZsIBrCMZQwgRJB+rEjZEvI1/uVLPNC3+SYKsb5h85Vl9GcZ/R3HQgS5lefZS0v2\n0Mo6BhG3lBBBNtAtpXWIZCQdGimlzIA+dYDaRHZ+XwHWV94rQIiIGCYAryulTgGbRGQDMAS8W6Lc\nRwv2ubgseEUJY1y/zU9R23K6q8cJa1hk//AOZCmNOcwE3uVdJnCAZhW+b7kF+zISwQnq+T4/SAVT\nAL1Yi0IooychguygHevoblmV/sAP2UqHjOtJuo4gIjVEZDmwC5iulFoItFZKmUuYu8BaxWkHjq1V\nW8EDueYQtzg/v+ARxwc+XtzUpQwkRJB3mUB/lnMen6ZVd6IIEpmyjxa+WqdyQSv2sJZejsXTdcaG\nJHPNJNvdc0mFoJQ6q5TqD3QAhopIn6jzChK2Is65kO1nUwpNrTh+yz0Or9AFDHXsOnNjeVQ3XxHc\nylQChPN+WJSIN/k2QUL8mBdcd/cVUM4iWlKDmXz5eUqflFeWlVKHgE+AccAuEWkDICJtwfIp3gZ0\ntN3WwTjmQtD20yWNJueeU9R25P0KEaSAckcYSDcydf/wiueZyJoc7O7ym6lMJEyAnpTF+FIFCNOW\nenyfzXz5eUqfZFajlqZFSETqA5cCa4H3gOuNy64Hphmv3wOuEZE6InIeUAAVtLXJQ37BI6zkAisA\nsLmD6wb+xFf5N404am3et5MPwYLbsJPXuC7l6z+nAWECFZo/ws3cW0qh5V9mumEXUE57tvER4/gj\nP8iqzmST5bbAKyJSk4ho/qaU+reIzAfeFJEfEhnXfBtAKVUqIm8CpcBpYKKqCB8Oj6nNKY7SyNWU\nOJhFRmyjJVZIF5MV9KMvq+jKxlw1NYadtIlxMU/EcepX+OYgMwROPEyXFNN1PZtJskky8+kqYKDL\n8f3AJXHueQR4JOuW5RHxnNjMlc8gIZZ4YPqMjHe9jQS3l5a8w1Ws5AK+y18TXvsBX2VhHkSie5af\npHRdd9azgW6e1FmpvU9zgdve5tt5xjLbmR/cC+MsiL3K99jvsjHHDa8259gxhxkb6OaIuudGRXfd\n/QxP3USHE1UgAAARZ0lEQVT0pIwerKMtO1hPd2s9J1u0EFyoYXPBtod8dMM0eSaeD6TmCLmIwY6g\nYV6wkzbcxnMx7iNuRO+ayzWb6cxBmlJEMVfxjiNti8k0ruQwjS2XFK/mYZXK18gPRjOLGpx1/EHP\nUoMiinmY+y2HO/s/I3qRLEyARhylPdtcVzmf5baUXL2LKGYb7TlBPU/GvRCZLK+mT0p29mOGA2FF\n0Y8VhAkwhpIYk+8+Wlj/i/4sZykDuZbXHXnZsqHaCCFAmLp8gUIYwVwrcptpmzYd+sygWQ/zAGeo\nQS1Ocw+/jWsqfYK7OUZDbuIlBFyFYI+BlAjTvz5ImFf5bqaP6sAUeHP2J23HVCZ6UmemmL5GEJmX\n3cZzlq9UL9byd75lfaEIKuVYsalQpYXQgOOMocQ17GMtTjsWaARlvbfnDT5NraTrBReyhO20p42L\nI9j3eDXp4pvJdtpRToEvVpt89kA1sWfMKaCc57iNVuxlIs+zll5AZKjXhp2ECbLdQ4/dKiuEhhxz\nLIKlQoggJYxxjKVPUythMLBkdazigpTNp5GFug2eu5s3Z3/CmEBHaRR3v3ZFsJ/mluFgDy0ppZCW\n7OUWfsc0ruRmXvRUBFBFJ8ut2ZW2CMA9QG86C0tmOBc7K5JYauwMZhEqTjuyIVl5T3JXXiXw+Cvf\ndbiPl1LILlqzll7czIs6nEuqmHFzMmEMJdYH52Hud43ts4Ch1nBrMYM4RsOEw5l/cBXf4J2kdZtJ\nCTfSNf2Gu9CSvdzK8wmv+Sdf96Qur9lHC1qyl720dLiJvMjNWgipYt8fkC72b8+anGEOI61odwCP\nM4kT1OUAzRjPh65paaNZSV96s5rurE94nTlZ3pjlItG57KaXsemnMmOGyAwQZi29rHmCH1RJIWQT\n5OohHkQhPMCvmMyjMR9yM2POAoY6IkYnYzV9UhKCF0OUH/NCyteupzuDWOzLpiCvCBMw8guFeDPi\nzeM5eT1H6McKfszv0rqnEUe5Jc177AQI8wC/iol3ZDKO6QQIJ13RvJ5XHPsRVtE3abQ2L0SQaqDc\nYop4mju5hjfyWgQmLdhHL9b6liwkb3uENuykJ+tId+E/29XReGbGEEHqcYJhRsArM/hvPLqwiUM0\ncYRzTze5t1+Y+YwP05ilsa5keclaejGbUdUvP0JECGtd0ygl4p98LcYLNBNmMdrxQQ8Sipk4J8p3\ncJCm9GOF9b4Pq6NCPvrDN/hH0muacpABLAPw5G+VCwKEfc2PkLc9Qk/KWEE/SjOYIP2Tr7vuD0gP\n4X4eZg4j+ZhLmMyjMYlCBrOIEsYktFI14wCX8QEFlCeszQwSlg3f4W9Jr/EzIZ9fFFLq6z51yOMe\nYQdtmWYEz82EbN0FRlNCTc5YwYIfZTJfGBNlk0TplJ7hdiAS8jGZCICs8iN0Zz3X8AY9KYt7zX6a\nV/iOuUzxc5Jskrc9QrYrq2b0g2yYzShHlLvHmESQEENZEJNeCYhxuks30l07tqcVwrA92xhDSUrW\nqHR2qeUbYQKcy252c65vdeRdj9CO7Un/sbkiRNAhBPOYfbV5EIutjJSX8YEjjme6kRXSSZoBkXFz\nvL/VHlpZdvjKHsViDb1j8lV4Td71CF4FdfWC6OQhbgQJ0Yij/Iuv8Xt+xOW8Twlj6MImX2MadWBr\n3CHX3/kWpUaEuKqCn4tpkIc9gpdkY5dPNjSzp46y51juxgaacIgrrXgGqZPKXCIeO2lDmADFFFU5\nEQC+P1Pe9QhespAhaQfNChNIyQ06TIDxfGi9DxKiMUc4TGNu4E8cpGnaPUI6QnDbuFPRm+4rM1W6\nR7Cnc02VVD9MCxgasztqIEusSXqmwyK3vRNu2Nc0dtCWF7k5o/o0Eaq0EIC0Q4inM8G9jA9ijg1g\nWVb5wLqluHdhHB8BkTnVS9yUcX2aCFVeCB9wWVrOcaliriG4zUOyyenWLYXh0fd41XqdD0HEqgJV\nXggAHzLeNU2qPbNmupjOeH4E7k2WwM++XyGbCbbmS6qFEABXIYxidsyxVB3johOOQ+Tb2YvoE2OZ\nnvK17diedX2aaiKEphzkDqYkvS4dP5zoaGxLuNBTq80I5sY9l45QNKlRpc2nJrfzTMLzZ6mR0eb1\n4zSgAccdJtdlDKCDI0VEZlzKDGpzCtDzgFxQ5YVgt/WbeOWB+VvuiTm2lIFxs0Qmo5wCGnDccrUw\n5yFd2ESYAJvowkjmOO7RIvGGKi2EluxlP80dm+3zlZe50XIvGct0x4S5C5usgFbn81/HfeuNzDGa\n7MgbIQQJWW4KO2jLOnpYkecy/dY7TgOChHiJmywhzOBSz9rsFe9zucPHajpj2UhXFML/8BcAOrOZ\n/8f/xdx7Iy/zMjfmrK1VlbwRQogg6+nOjbxMW3bQlh3W0CBIiD20SnuPQQOOU5/PrTnCF9TNSbrU\nLXSiE1sSXvMOV3Ga2hygqesiXqohXdqxnZt4mZe0GLIir6xGZvfvRiv2JI3RE020i8VJ6mTctnRY\nxoCk1+ymNWvpyU7aJEw2XkwRr3A9YQL8h68AsfOCtmynB+uya3Q1J296BIgE0N1Ge8sCEyRkxeiB\nxHuE3YiOWZoP6VLBjNyW+iaTTXRhk5FnbhajgUi+gx/yBwDe4Fory6QmM/JKCNG0Zpdj+2G6TnTH\njXxgy+nviLTsN8lcLAopzXrH1VY6UEwR3VmvJ8wekFdDo2jK6Ol435K99GF1WmWECOZUBBAJk5KM\nH/OCJzF6tAi8Ia+FsMslxmW2iaVzwTIGpCQ+P8OTaNIjr4Wwk9YUU8QL/Ng6Vtn339pp41PUNk36\n5JUQ2rGdAGEGspQiigkQphdr2c25/J1vxQyV8plnuD1pr7BT9wh5g1REGmQRUbi4OZhCiI7MsIbe\nvMXVuWqep9gXCqMpoxdhxvgS5rx6U4xSKq34mnljNWrPNn7E7633dlv5bs4lQJg+rOZ5bq2I5mVM\niCB9WU1zW+ILk56sJeL4rXiDa3LfOI1FSkIQkZrAYmCrUupyEWkO/A3oDGwCvq2UOmhcOxn4AXAG\n+KlSKqnP8Pf4K12NtKrbaUcJY1hHD8c1ZhK523iO57gtpYerKPqznABhK9pdokDGiaLTaXJHqnOE\n24FSvvyPTgJmKKW6AzON94hIIfAdoBAYD0wVkYR1fI9XHSJ4mRtjRACRWD0QmSxfEJV6NN9YTn+b\nCCBMMGnszmt53e9maRKQVAgi0gH4KvB7vsycfQXwivH6FbCClE4AXldKnVJKbQI2QPzd86OY7Ui0\nl2hjSyGlQGTMvZK+KecByAdW0jdpL9ad9b5s+9SkRipDo6eBe8CxStRaKWXa/naBZf5oB47Y6Vsh\nfiLci5lpvd5GeyuTYjSt2ENv1ljv7bl27SlJKzv5FOWvupFQCCLydWC3UmqZiATdrlFKqYgVKC5x\nz5kT4q10ZIOLt2U/VjCCuZzLbsfx17iO63iNfqyoUkIIEK7UwXorM8l6hBHAFSLyVaAe0FhE/gLs\nEpE2SqmdItIWrE/qNqCj7f4OxjEXQo60fW6soB8r6Edv1lhJNipiR1aAMPMYnpX3aipzgALKPc+x\nXD3YZPxkTsrrCCISAO42rEa/AfYppR4XkUlAU6XUJGOy/BqReUF74GOgm4qqJN46QjzMoZC5xbIF\n+/kur/K//DTlMtLF/ECadZuECLKAoZygXspl/ITnXM2n8aiMyTzyC//XEcwP9GPAmyLyQwzzKYBS\nqlRE3iRiYToNTIwWQbq0ZjdhAo58B/to7psIHuQhy5/JLT1rPU5wO89QjxOECLp+e9/BFMsDtTn7\neZbbuI3nHEm0E9GBrZ6EhdGkTl6tLJu0ZQc7aMuVTMsqa46XjGAulzIDgMe5l2HMd6wY2/dQmKyg\nH9O4kr6s4qoUEo7b0b1CNlTilWU7kcjShymjp/Vhs3/zmscO0jRnk+W5jEAhjDXSy37EOAazyNr8\n49Z7TONKLmBlSgn+ounIZ3zmmG5p/CQvhRAiyE28RBt2Wsd6s4Y19CZMAEFZYggSsjbf+MFI5jCH\nkcCXoRiHMZ9hzCdEMG52e3NBrSsb2U9zmrM/5To30lWLIMfk5dDIxPyQtWKPtaAGkY3vq+hLMw7w\nU/4XiMQT2kZ7z/MG38vjrvnSEjGFOxwRsa/jtbRilE5nbN5sK62cpD80ymshpEJTDhIgTH+W+zau\nts8PkpEs0chEptKKPQnL2Eg3XuW76TRR4yB9IeTVfoRMOEhT3mWCr5PLeMMfN5L1SKnER93I+SnX\np/GGSi+EXFCbUyktphVTlHS/sqCSpr41J+Sa3KGFkAKnqJ3Sam+AMAHCCaNYtGRv0qERRHoht7D1\nGn/QQkiBR/gFDTnGY0wiTMBaVV7EYEdIliAhmnCIQzQhSIgiimnGAevbvSV7PUmErvGevDSfVhQB\nwgQJWWHig4T4NfdxmlocoyF38jQg1OWENSmOdsEwMdc4OvGZ1Zu0ZC+lFLKG3g5v2nhczEw+oyOb\n6ezlY2pc0D2CjTABFMIcRiIoHuZ+hjGf2pxiLiN4jElM4Q6HZSjeJH05/fmU86z1gL6soiZnAFIS\ngcn3+TNBPV/wHd0jRBE2Rvom0XnWTlDX6gncFtTsVqHoVe819M5oEhwgxGY68SnnpX2vJjW0EKJI\nNCkexnyChKnLCT5kPKvo61jlXsrAuFYjMx5TOqZYk0100SLwmUq/oJZLfsZTnMMR6/1xGqAQnuBu\nLmQJS7gwrfJ6UUZrw41kPd3ZTjtH+JdPOY9ZjEYhVhBgTSpUEae7fOUpfhZzrCHHCBLiKOekXd5a\nerLWJWjZRrpSh5OECLKFThm1VZMeukfIQ7qxkS+oox3vMkb3CFUCt/3bGn+pwubTTboOXUfKaCHo\nOnQdVGkhaDSpo4Wg0VChViONxj8qxQ41jSbf0EMjjQYtBI0GqAAhiMh4ESkTkXIRudfDcjeJyEoR\nWSYiC41jzUVkhoisF5HpIpJWnlkR+aOI7BKRVbZjccsUkcnGc5WJyNgMy/+liGw1nmOZiFyWafnG\nPR1F5BMRWSMiq0Xkpz48R7w6PHsWEaknIgtEZLlRxy89fQ6lVM5+iET73QB0AWoDy4FeHpX9KdA8\n6thvgJ8br+8FHkuzzNHAAGBVsjKJJEdZbjxXF+M5a2RQfhHwM5dr0y7fuK8N0N943QhYB/Ty+Dni\n1eH1szQwftcikn5gqFfPkeseYQiwQSm1SSl1CniDSHIRr4i2FMRLaJISSqlZwIEUy0wrSUqC8iH2\nOTIq36hjp1JqufH6KLCWSIBmL58jXh1eP8tx42UdIh9w5dVz5FoI7YHPbO8TJhJJEwV8LCKLReRG\n41i8hCbZkChJylbbddk8209EZIWI/MHW1Wddvoh0IdIDLcCn57DVYSaM8exZRKSGiCw32jtdKbXQ\nq+fItRD8tNWOVEoNAC4DbhWR0Y6KI/2lp/WnUGYm9b0AnAf0B3YAT3pRvog0At4GbldKHbGf8+o5\njDreMuo4isfPopQ6q5TqTyTvxlAR6RN1PuPnyLUQohOJdMSp2oxRSu0wfu8B/kGkG9wlIm0AohKa\nZEO8MtNIkhIfpdRuZUAkb53ZnWdcvojUJiKCvyilpvnxHLY6XjXr8ONZjHIPAZ8A47x6jlwLYTFQ\nICJdRKQOkQyc72VbqIg0EJFzjNcNgbHAKqPs643LrgemuZeQFvHKfA+4RkTqiMh5QAGwMN3CjX+m\nyTfASiGaUfkiIsAfgFKl1BQ/niNeHV4+i4i0NIdWIlIfuJTIXMSb50jHiuLFD5Ghyzoik5fJHpV5\nHhELwXJgtVku0JxI1p71wHQimX3SKfd1YDtwksjc5oZEZQK/MJ6rDBiXQfk/AP4PWAmsMP6prTMt\n37hnFHDW+NssM37Ge/wcbnVc5uWzAH2BpUZZq4D7k/2P06lDu1hoNOiVZY0G0ELQaAAtBI0G0ELQ\naAAtBI0G0ELQaAAtBI0G0ELQaAD4/9U5yXX65Y7QAAAAAElFTkSuQmCC\n",
"text": [
"<matplotlib.figure.Figure at 0x1129a9fd0>"
]
}
],
"prompt_number": 91
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" clim_file = '/sharehome/savoie/projects/nsidc0081/ancillary/icelakeclim.nsidc0622.N01_v1.byte'\n",
" clim_mask = np.fromfile(clim_file, dtype=np.uint8).reshape(rt.variables['loci_flag'][:].shape)"
],
"language": "python",
"outputs": [],
"prompt_number": 92
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"imshow(clim_mask)"
],
"language": "python",
"outputs": [
{
"output_type": "pyout",
"prompt_number": 96,
"text": [
"<matplotlib.image.AxesImage at 0x111137ed0>"
]
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAMIAAAEACAYAAADhtqk6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXecXFX5/99PeiM9pBdINmWTkEJ6m0EgCSoEERXQ7w9R\nKQYUEJBEyrqIFKWELxAEbHxFQAQMoAIJkZlN732zySaShPTeCWnn98fce7l35k6/d3Z297xfr33t\nzC3nnLs7nznnPOc5zyNKKTSa6k6Nim6ARpMPaCFoNGghaDSAFoJGA2ghaDSAFoJGA/gkBBEZLyJl\nIlIuIvf6UYdG4yXi9TqCiNQE1gGXANuARcC1Sqm1nlak0XiIHz3CEGCDUmqTUuoU8AYwwYd6NBrP\n8EMI7YHPbO+3Gsc0mrzFDyFonw1NpaOWD2VuAzra3nck0itYiIgWi8ZXlFKSzvV+CGExUCAiXYDt\nwHeAa2MvK/KhajshIKjrqJZ1FKd9p+dCUEqdFpHbgI+AmsAftMVIk+/40SOglPoA+MCPsjUaP6jC\nK8tddB26jpTRQtB16Dqo0kLQaFJHC0GjQQtBowG0EDQaQAtBowG0EPKKAsoB6MBWzue/AHRmc0U2\nqdrgy4KaxklfVtGc/YQJWMeu4h3OZQ/TmMBO2gBwHa/F3Ptnvp+rZlZrtBA8ZCBLCRBmKQMBCBJy\nnA8SYioT2UMrAMroYYkA4DWuI0AJ7Q0fxZlczGY6u9Y1jPk04wAfcJkPT1L90ELIkgBhBrGYRhwF\nYBkDYgRgZyJTmcpEyuhFbU46zpVTQIASAP7L+cxmVNxy5jMs+8ZrLPQcIQuGsJAgIUsEAK3Yk/S+\niUwlQJgV9HMc/xG/pz1b2UhX/sL/eN5eTXx0j5AlT3A3F/EJAAJ0Mya8yVH0Y4UlhgLKac82AMK+\nuytrotE9QhYsZAjHaEgP1hEiyPt8nRBBDtEk6b27aM2VTAOgO+utifIGuvEZHXxttyYW3SN4wJPc\n5XjfhENxr32Lq9lDK0YxG4CiqE0kG+jmfQM1SdE9Qo7pxVp2c67ruXIKWMDQHLdIA1oInjKQpVzB\newmv6c0aerOGd7gq5lw5BX41TZMEPTTykMt5P6XrruYterDOcaycAhYx2I9maVJA9wgecSFL0rq+\njF6O969xnZfN0aSJ7hE84BZ+R2t2pXx9KYWU0otdtKE1O31smSZVdI+QJYNYnJYIICIEQIsgj9A9\nQhYMZhEBwilfv5o+vM03Y46vo4eXzdJkgO4RsiBIiIYcS/n6Pqx2uGCYq8rRE2dN7tFCyIIGHE/7\nnnR6EE3u0ELIkCEszOi+3qzxuCUaL9BCSIF7+C2TeZRzOGIdS+RqnYxoMayne8ZlabxBCyEJtTjN\nb7mHOpzkCOcAMJQFfE79jMu8mrcAmMaV7KQNIe1tWuFoISRgMo9Si9P8gkcAqMsXAIznQ5qzP6uy\n+7IKgJ20YQdts2uoJmu0EBJQh5Pcy+PU5hRPcDcAk3jc83oGstTzMjXpodcR4mD2AiaDWJzVvECT\n3+gewYXRzKI2p6z3cxnhuQiu4h0AyuhJd9bTju2O89rMmlu0EFyI/tDPYSQzuNTzem7hRa7hDXqw\njht52To+lul05DMuZqbndWrc0UMjFxTO9FvHacBcRjCGEmvC7AXzGGZt14TY3Wpd2chMLvasPk18\ndI/ggj0QVyrHM8UugmjMHqiI4hiBaLxHC8GFeOPzeQxPaWO+F5ghHwFe4fqc1Fmd0UJwQRKkip7C\nHTlpQ1c2ApGQj5tykm2meqOF4EINzlqv5zAy5vxHjMtJO/7EDXFDPmq8RU+WXThLDWpwllmM5j98\nJeb8fIYxhhLq87lvbZjBpWyhk2/le0WAMCWMYRSzqcFZanOK08bHai4jOEmdCm5hamghRGEfFo1m\nFnMZwQnqxVwXJsB4PvSlDRvpylxG+FK2lwQJEyBk/c0qs8+UFoKNIKGUF7IWMJQxlGS0JyERG+jG\nX/mup2V6zUV8whhKCBPgIYoSzKgqD3qOYCNEMGYNwa03MDlOA8/bkO9xT2tyhlHMJkTQ+HtVDZIK\nQUT+KCK7RGSV7VhzEZkhIutFZLqINLWdmywi5SJSJiJj/Wq41wiKIoqtbn4WoymmKOE9Xsch+jdf\nZSvtPS3TS2pxmsk8yq94wPM1lYomlR7hT8D4qGOTgBlKqe7ATOM9IlIIfAcoNO6ZKiKVotdRiOOf\nO5pZ1ORMwnsWMsSz+itLgK+Hub+im+ALST+kSqlZwIGow1cArxivXwGuNF5PAF5XSp1SSm0CNoCH\nnxafiZ7s3c/DCa8fzCLf6s5HTlfhKWWmT9ZaKWUG89kFtDZetwPm267bCnnc19sQFA/yUFr3ZLpv\n2c46ehAmoDfnVDBZS1wppUQk0ZypUs6n4q0hZMMaerOXluzmXPbRkl1xomJrck+mQtglIm2UUjtF\npC2w2zi+Dehou66DccyFkO11F+On4lAIxRRZDm6jmcUZasadFA5hIS3Z63ouTIADNCNAmGbGqPI5\nbmMfLfxpfLVnk/GTOZkK4T3geuBx4/c02/HXROQpIkOiAog3fghmWHXuCBJiAUNdTahuG3XCBDhM\nYyurZnSONI1fdMH5RZr+pqZUzKevA3OBHiLymYjcADwGXCoi64GvGO9RSpUCbwKlwAfARKVUpRoa\nRa8jAIxgruP9UBa4ulcoRGe8qaQk7RGUUtfGOXVJnOsfgagNv5WEAOEYz9N7jc36K+jHMRoCkRzH\nbgQJsZGuHKaxvw3VeE6lsPFXJrbqRICVEi0EG27jftMN2+wNAOpxwvX+3/MjX9ql8R8tBIN4m3HM\nOYM96nU8IRzRQ6JKS9VdKkyDAOG44VrMNLAmdzAFiLhXRC+oHTZCQmoqH9VaCDU5k9SNAuBjLrGG\nRvWMKBYhguyjBZfxAQBPReVa1lQuqr0QEvEokxnBXMd2zbrGsMgMAhwmwGIGcZRG/jVU4zvVWgip\nbCO0O8MNZx4Aj3Mv4K33qaZiqfaT5VmMcbw3J8ezGB0jlHkM5wvqxp0sayov1bpHqMkZRlPiOPYQ\nDya8J9GONU3lpVr3CMnmCG4sYBgHaZr8Qk2lotr2CPU4YblPmKSy/XAew/xqkqYCqZZCuISPGcmc\nmOOVYZeYxh+q3dAoQNhVBJrqTbUSwjkcibuC7OZ+rak+VJmh0QCWsYFu/Iyn+CdfpxFHOU4DWrAP\nhXCCegmz3lS18CSa9KgSQmjCIa7gPev91/ln2mVoIVRvKqUQmrOfn/AsB2jGZjpnPaw5Q02PWqap\nrFQqIQQI04J9Vo7iZhywNsdnyllqVNmgVZrUqTRCaMUezzNbnqQOjzLZ0zI1lZNKYTVqxR4mMtXz\ncutwkgBhRjPL87I1lYtK0SOca4VN8h6zl7Fnycn3ifMP+CN/5AcV3YwqRV4LoTW7uIXf5aQuUxDJ\nnO5ySSe2sIVOdGYzXdjEFjrRju105DPG8VHOUlhVB/JaCF7PCZKRLAx8qvRnORN417XMQkoJEGY2\no1hFX8e5dmxnO+2s9zfwJ0/ao0lO3gmhHdtRCDtoS0/KKro5aeG29/kbTCPMGNqxnRbss85fxTtc\nxTsOodhFADCbUTF7pk020tXTtld3pCIC0UWCBrt/+wYJcZjGFFBeIUIIGZnBUuFWnqcle1nMINqz\njbbsSLu+MnqykzYxx7uykY58Fvc+r3qvqkkxSqm0FpfyrkfYRnvO4UjORWAK4EEeSlkIz3MrQUI0\n5BjncDSjentSFvOsW+lAB7YmvC96jtCDdayjR0Zt0OShEMopSPub9RBNaMKhtO4JE+AMNZnLSEYb\nifEg/clyiKAVQdsrPmIcW+mQsNxhzKcF+1jEYMop4BresNqT71avfCTvhHAtr/M68cKtxnKQpjzD\n7Uk/jIdoQimFDGcex2ng2HuQyT6EISwkQNjzrJoAP+QPAGyhE53YEve6AsopoNxxrBdrtRAyIO8W\n1JZwYVrXN+Wgw+HOjce5l2UMsKJQZPvhHcoCLuMDX0RgpxNb2EzntO5pzS6KKCZAmN6s8allVY+8\n6xEChDmSYsS4gzRlBf1YxgAADtOYRhzlQpY4rovekglQRLFjwvkgD7kOi+x5l4czj7pGgC+vKafA\n+nbfSFe20sHqqb7Pn+nM5rTKM61TgmI77Sz/LN1buJN3QggRZAdtWUcPerDOcW4vLR1Zap7hduv1\ne1xhvY4WghvR0Sh+xQPW8OosNRwrzX6zjh68wTVARKCv8j3H+U10SVsIgJUBKESQAzSLG99Vk4dC\nKKcAiJgV3YTwPLfSgv30ZWXcMn7HLUlXpM3YRDU4y2hmOez/uRQBRCw+9jlOgDD/5XxqcZpPOc/R\nK6VDkJA2s6ZI3gkB4HLed805bPYG+2iecIK7i9YsZhCDWJywnkytPfZca34QNPLaZ4MWQHrkpRDe\n53Ig0iuYNvYwgZStO4046tisc4J6PM691OULvqBuVh/ixyO51fOa1fSp6CZUOvLOamQnTJC19AJg\nD63SuteeDNyMVfoFdQH4Lfdk3KYAYV97A03FkNdC2Elr3uTbvMXVrKF3RmV4bSUZZphg85k+rKaQ\nUl/KznS+ku/ktRBMMhUBYJkN7URnyayKfIu/cyvP04fVnpRXi9OMZA5n43xk7uPXGYXQzBfyzunO\npD/LWU7/rOppwHH6sooFDAUicY1+xlOOa8weo6p+09knzRewkpVc4Es9tTjN6agpZx1OphR633uq\ngNOdyU7a0JMyyuiZcRnHaWCJACI9QZgAp6jtSP5hUtXEsNeYVzVnPwHCXMBKvsE/2Exn/sz3Uy6n\nIce4lef5DT+Pe81oZvEJFzGJxxyLjvMYzkKG5H3g5LwVQht2MoF3KaaIQkoppTDrMhPt6KpqIgBo\nyR7XiX1nNnMHUxAUYQIsZaB1rohiPqc+YQKOLxE79/Frfs19jmNnqEkRxYQJMIdRnDI+WsOZx1AW\n5P1uurwVwufU50Vusf6RftvFSxjDmKhcCZWZPbRiKhMJEqI5+2PmSqa37uW8z+W87zhXn88Zz4cI\nitX04S6e5HPquw5/TMxdd9Em7nkM9+6hfCTpZFlEOorIJyKyRkRWi8hPjePNRWSGiKwXkeki0tR2\nz2QRKReRMhEZm0nDTN9603z6Lf6eSTEp8wkXVUk/nFX0dTUYpMI4PuIunmQuI6nNaQKE+Qr/YRaj\nAQgapuQiivk19zGTi71sek5JxWp0CrhTKdUbGAbcKiK9gEnADKVUd2Cm8R4RKQS+AxQC44GpIpKR\ndWonrVlNX/bQikJKffemrEph4Vuxh6t5m76sitkbnS4jmEMtTjGK2ZyhJiWMoQZnCRDiIUMK8axJ\nlYWkrVdK7VRKLTdeHwXWAu2BK4BXjMteAa40Xk8AXldKnVJKbQI2QOZZ90rpxVQmMpWJtGKPZ+bA\n6kAZPQkRzFoIgDUkChKiiGLGUEIxRa5ufPfxiCNBe2UgrTmCiHQBBgALgNZKqV3GqV1Aa+N1O2C+\n7batRISTFXtoRYggt/Gcry4ECqkyXppePkUtTgPwH75iDY3cmMRj1OKUlZe6spCyEESkEfA2cLtS\n6ojIl2ZapZSKrA3ExeVcyPa6i/Hjjt1q5Pc4vqqIALC8rbzKAjqHkQlFECBcQYtqm4yfzElJCCJS\nm4gI/qKUmmYc3iUibZRSO0WkLVjh6LYBHW23dzCORRFMqYGt2e0QghfdfCKqSo9QSiGr6U0L9sVs\n58yUU9Tmfh5mFqMRFAFKCBlfTHZv2UTrDf7QBecXafqm8KQryxL56n8F2KeUutN2/DfGscdFZBLQ\nVCk1yZgsv0ZkXtAe+BjopmwVpbKybNKXVYxiNuey23cTaoAwAcJVQggAa+hD7wqYU/2Gn/M59XNe\n75f4s7I8EvgesFJElhnHJgOPAW+KyA+J9EvfBlBKlYrIm0ApcBqYqLLw49hNa8IEfHMiMxFUziPr\n+U1FiKCyklQISqnZxLcuXRLnnkeAR7Jo15dlGb/N9QQ/cItQp8meAOFKszaT18bf1uzim7wNZOeB\nmozqIIL9NE/pus+pTzFFcd0rUuE0tSmi2JOddrkir4UAkc05XvgZxSPX+5MriubsT+m6+nxOgDAf\nMp7HMtiN9xiTuNPm4VtZXCzy1g07F9TkDPfzcEU3Iy8xt7fW4SSTeTSrsnK/fzr9yXLe9wh+Upk3\nkvhNPU4wjPmcpA7FFBlepbGu63ZmM8raDhtdVr6Tt96nuSDbb7qqzjg+QlDMY7jlh/Uxl1CDs4yh\nxOG6bgZXsO8VNwkS4kPG56rZGVGthVAdeZcJdGYz/Vme0vVjmc5FhJjDCCDygT9LDWMaHIy53q2X\nzXcRQDUeGl3MzIpuQs45QDMChNlCZ6ZwB4dpnNJ9tTlpWYDMuKoP8lDMZqYAYcsnyc44PvKk/X5S\nLSfLDTjOPfy2wuqvSMIErP3gjTnMAJYlvykFQgQTmkpzG65eT5Y1SejHCkIE2UhXz0QA8ddiiini\nBPUIEqIF+zyrz2uqXY9wKTOqRTiXVFjKQAay1Nc6QgQZxnzLchQiyGIG+eymrXuEpGgRRHiXCTkZ\nqgQJWSL4nPoECXE3T/heb7pUOyE8yuSKbkJeMJClvvcGdn7Dzy337Ce4O2f1pkq1E8JJ6jDXMAVW\nZxpzOKchbPqyigBhiiliEIvzLnyOXkeopqSbfDFbLuMDwDmpzifP1GrXIwDM4NKKbkKV4lc8kHYU\ni8UM8qk1mVEthQBfhorXZM8YSpjNKP7DVxJet5AhhAnwb75KT8pox/YctTA51XZo5JZgsLqxjAGe\nrCUECPNr7uM0tVAI9TiBoDhFbWsusJAhfMBl1j0K4UZe5mnu5Hz+y6ecxyGaZN2WTKl26wh2LmYm\no5hd0c2o9JQwhk+4yPXcnTzNWnq5+hv9jKc4hyMA/Jnvp51KNz56HSFl3EQQtrbvR37OUDMm2O18\nhuXVJC8bzLS82RAmEFcEEPl7xXO6M0UAeCiCzKiWQ6PanOIMNR3HFBLjTWm+P04DR3Jx0/MySCjv\nzICp8h5XsIwBKCSj9QSFUELACucSj3g71G7mRev1FO4AYnNf55JqKQRwugufpQa/4oGU7x3HR3zE\nOEIEOUZDvsq//WiibxyiidUb2JMupkKYADU4yxlqZtUzNuKo9foOpgCwgn4Zl5ctVVYIg1nkmqIW\nIoGqwgSsoVENzvIgD/EQDyYscz7DrB1YRRRbPYWZ1/kI5zi6+/wl8uFvzzYrOXuYAA047rqxxmQm\nFzObUVnXfjdPUEZPjtLIinTxKedV6PCoygohnghMouP8C8rKxRAiyCxGu9rG7cOlFuxjCAs5RkP+\nxddZR/eY1FT5ScRAso32FFNkDe/i7TAz8UIEg1nEYgZZf8d8iUBe5YTQiKP0YB1LuJBbeZ7nudX1\nOrd9tNG5nEcxm9mMsuYHw5hPiCBfUJehLLAiQwQI05N1fI1/WvcepClNOejhk/nHEi7kLp7MSeiV\nZF9QFUWVsxoJ0JYdBAjTkr0MYWHSe85Sg4d4MGbiO5tR3MevrffFFFkj6ubsZxGDOUZDBrGYRsaQ\naBkD+BM3WDnD7Mn7wgTYQytKKbQsUwDr6Z7FE6ePXewBwhylUV46wuWSKieEIzSyLDp7aRk3aPC9\nPO6Y8CkkxmIxitmOrYdjmc4JY47wLD9hD62sPABLuBCAASxjC52sey5gJWECFFNkfQD/zrcsy9OL\n3MzrXGtFitiWQgR9u4gyYQLvcgdTCBDmHI4QIMwxGvJEkkTsI5mTcZ35TpUbGgHWN9xEpiYMRhsm\nwCxGW1EZShjjOD+bUY69zcOZx3DmUUwR+2lONzZY58xJp7kPeDOd6cImVnKB4xs4eky8kzaOc9vo\nwBjCdGCra3vt95dTQAHlGZlwm3DIMRRKZVhUEcGRe7OGq3kL8Dc+UpXrEUyO0ZBFDGYoC1zPP869\nlukw0X5bt1g+4/jIWnZbzCB2WTlSIu7NdzDF8YG1R5lL9oErpxthlwnkK1wfI6JttHcMq/7Gd3xb\n7DtNLU8my+liigAiIUD9okr2CCZuFok6nGQQi6nP51Yy7LPUiPsBupAlzGWEY2fbMObzLD8hSIgQ\nQQTl8LH/L+dbr//BNxzlmZkuFRK3zgBh/sv5fEZHurCJTXRhk0siFTOE/cvcxOW8Rxk9qcPJpH+X\naI7SyGHXtxMmYDxLeusNXvE2V7Pax7i3JlVaCG7U5hSXMsMKaWji9qEcxWye5C5OU4sZXOrwjfkJ\nz1rXHTXmJQHC1tCoS4IMLolMhkUUM52x1orseXzKp5xn/bZjb3O2ZsjIUPIFFjKYE9RzBAGey8iU\njA5+YIqgF2vZT3NH7+slVXZoFI9T1AZwjb8TzTyGc5pajGAuAcIx37am1QjgLp5iKQN5mjsZwDJL\nCN9gmjUcasneuHUFCHMdrxEm4HBLMD/80SKIxkzHmwmNOMrdPEEDjjn2GJucpDazk4R7TIe7eSKt\neU0hpQQJ+SYCqIZCMLdq1uJ0Us/TM9RkJHO4lBkECVGXLxznB7OIgJFrOGI+/XL4sJz+HKAZB2hq\nWYL20tI6f66VaQtu4iWO0ZDXuC4vFpiChLidZ3wpeyJTaWgILtVV+N6s4Vx2cwErfWkTVEMhwJc7\n1C5mZsJvplHM5hI+BiJDhyOcwxTuoJgiyxoVJEQxReygLQMNy9EWOnGIJpbFaD8tYspuwT4rWXdb\ndvA1/uX1Y+Y9qfYK9mxJPSmjDTs9b0ulFEIB5fyUZ+nHCibwbsblFFPEbEYTIMwYShznRjPLYTo9\nRW3mMdzaPBJtlj2HozzNz4CIaXI5/a1v930uSTqCUR+Ccgoyfg47Xg0fmnKQYY4swd4wlYnWcPJC\nltDeLc9kHFZyAWX0dJicvaJSTpbLKeA6XuNKprlaU1JhJhdTkzMohNmMinHLHs0sx/uT1LHG7qZP\nUpgAg1hMEcU8xV0coREQEcLtPBNj9+7NGnqxlt6sYS29eIFbMmp7Irw0MQYJMZ9hnpVnEiZAQ44R\nIMxAlqa0iPg23/S8HXYqpRAg8m3ejAMcoFlG99tt4kUU8wkXxSyomUSHf1nAUGuzSYigNUcwhbCF\nTnRiC0UUs5J+/IMrAadN3C9zpJcTyrp8wQjmeh7+xvQ3MmMrlTCmQrdpQiUdGpmYIsgmNImZOsq0\n7DzIQ9TkDDVtKaUUwqXMsCbG0TuudtGam3jJ6ub/xA2Wy8QFrLAW3+z08ilL6C5ae7oCO5x5npUV\njWn+Nd09KtJJsVrvWTYxhzqJsPskQSTXmBnC0M7T3OkIt/5T/pdmHIgpr5gi2rALUL6MeSHSA/Vm\nTVZlnKSOr9EB3f729sjZ9vMhgtZkOUSQcgrYTjuXUv3Js1ytUAgKYQ4jGcFcanKGR/iFtf5wB1Mc\nPdBiBvEvvgZENqPfydMsYwDvcYVxhfv/oy+r4joEekEhpVmLALBW3/3iMI1pzGHHMTNydlt2WMdK\nKeQAzSx/pyAhjtEwjhDSRwuBSICqB/gVEHEsM3eqDWUBjzHJsYnHFIEZo+c4DQAcbggDWMYAlnGA\nZq69QS4opZA19M5aDH5nxQwT4HLedz1nLkCGCdCLtQgKhbCOHrzBNZ62o1IJwa8E1vb9CPY9vNFp\nkOy7z0IELRPqXTzJS9zsGPM24VBCEVzFOxSylp20Zj3d2UFbLx8JgLe4mr20zOsAA0sZiEK4gvcc\nx/uyCoA19CFEkN2cyxhKaM0uFjLE83YknCyLSD0RWSAiy0VktYj80jjeXERmiMh6EZkuIk1t90wW\nkXIRKRORsV42NkzASkDuNWYUizABLuITHuBhhxCGM891JfTHvMCL3GxZjCAyWU5FsD1ZS5AQN/ES\ng1jszYNEkU1yjumMzUnA5MhQcoK1mcncb7GPFqwxcmyXUsjvuIUwAa7gPcewyQsS9ghKqRMicpFS\n6riI1AJmi8gHwDeBGUqp34jIvcAkYJKIFALfAQqB9sDHItJdKeVZVu8QQV/H1/bFtYe531pfsE+K\niymiPp8DEVPgUZsIBrCMZQwgRJB+rEjZEvI1/uVLPNC3+SYKsb5h85Vl9GcZ/R3HQgS5lefZS0v2\n0Mo6BhG3lBBBNtAtpXWIZCQdGimlzIA+dYDaRHZ+XwHWV94rQIiIGCYAryulTgGbRGQDMAS8W6Lc\nRwv2ubgseEUJY1y/zU9R23K6q8cJa1hk//AOZCmNOcwE3uVdJnCAZhW+b7kF+zISwQnq+T4/SAVT\nAL1Yi0IooychguygHevoblmV/sAP2UqHjOtJuo4gIjVEZDmwC5iulFoItFZKmUuYu8BaxWkHjq1V\nW8EDueYQtzg/v+ARxwc+XtzUpQwkRJB3mUB/lnMen6ZVd6IIEpmyjxa+WqdyQSv2sJZejsXTdcaG\nJHPNJNvdc0mFoJQ6q5TqD3QAhopIn6jzChK2Is65kO1nUwpNrTh+yz0Or9AFDHXsOnNjeVQ3XxHc\nylQChPN+WJSIN/k2QUL8mBdcd/cVUM4iWlKDmXz5eUqflFeWlVKHgE+AccAuEWkDICJtwfIp3gZ0\ntN3WwTjmQtD20yWNJueeU9R25P0KEaSAckcYSDcydf/wiueZyJoc7O7ym6lMJEyAnpTF+FIFCNOW\nenyfzXz5eUqfZFajlqZFSETqA5cCa4H3gOuNy64Hphmv3wOuEZE6InIeUAAVtLXJQ37BI6zkAisA\nsLmD6wb+xFf5N404am3et5MPwYLbsJPXuC7l6z+nAWECFZo/ws3cW0qh5V9mumEXUE57tvER4/gj\nP8iqzmST5bbAKyJSk4ho/qaU+reIzAfeFJEfEhnXfBtAKVUqIm8CpcBpYKKqCB8Oj6nNKY7SyNWU\nOJhFRmyjJVZIF5MV9KMvq+jKxlw1NYadtIlxMU/EcepX+OYgMwROPEyXFNN1PZtJskky8+kqYKDL\n8f3AJXHueQR4JOuW5RHxnNjMlc8gIZZ4YPqMjHe9jQS3l5a8w1Ws5AK+y18TXvsBX2VhHkSie5af\npHRdd9azgW6e1FmpvU9zgdve5tt5xjLbmR/cC+MsiL3K99jvsjHHDa8259gxhxkb6OaIuudGRXfd\n/QxP3USHE1UgAAARZ0lEQVT0pIwerKMtO1hPd2s9J1u0EFyoYXPBtod8dMM0eSaeD6TmCLmIwY6g\nYV6wkzbcxnMx7iNuRO+ayzWb6cxBmlJEMVfxjiNti8k0ruQwjS2XFK/mYZXK18gPRjOLGpx1/EHP\nUoMiinmY+y2HO/s/I3qRLEyARhylPdtcVzmf5baUXL2LKGYb7TlBPU/GvRCZLK+mT0p29mOGA2FF\n0Y8VhAkwhpIYk+8+Wlj/i/4sZykDuZbXHXnZsqHaCCFAmLp8gUIYwVwrcptpmzYd+sygWQ/zAGeo\nQS1Ocw+/jWsqfYK7OUZDbuIlBFyFYI+BlAjTvz5ImFf5bqaP6sAUeHP2J23HVCZ6UmemmL5GEJmX\n3cZzlq9UL9byd75lfaEIKuVYsalQpYXQgOOMocQ17GMtTjsWaARlvbfnDT5NraTrBReyhO20p42L\nI9j3eDXp4pvJdtpRToEvVpt89kA1sWfMKaCc57iNVuxlIs+zll5AZKjXhp2ECbLdQ4/dKiuEhhxz\nLIKlQoggJYxxjKVPUythMLBkdazigpTNp5GFug2eu5s3Z3/CmEBHaRR3v3ZFsJ/mluFgDy0ppZCW\n7OUWfsc0ruRmXvRUBFBFJ8ut2ZW2CMA9QG86C0tmOBc7K5JYauwMZhEqTjuyIVl5T3JXXiXw+Cvf\ndbiPl1LILlqzll7czIs6nEuqmHFzMmEMJdYH52Hud43ts4Ch1nBrMYM4RsOEw5l/cBXf4J2kdZtJ\nCTfSNf2Gu9CSvdzK8wmv+Sdf96Qur9lHC1qyl720dLiJvMjNWgipYt8fkC72b8+anGEOI61odwCP\nM4kT1OUAzRjPh65paaNZSV96s5rurE94nTlZ3pjlItG57KaXsemnMmOGyAwQZi29rHmCH1RJIWQT\n5OohHkQhPMCvmMyjMR9yM2POAoY6IkYnYzV9UhKCF0OUH/NCyteupzuDWOzLpiCvCBMw8guFeDPi\nzeM5eT1H6McKfszv0rqnEUe5Jc177AQI8wC/iol3ZDKO6QQIJ13RvJ5XHPsRVtE3abQ2L0SQaqDc\nYop4mju5hjfyWgQmLdhHL9b6liwkb3uENuykJ+tId+E/29XReGbGEEHqcYJhRsArM/hvPLqwiUM0\ncYRzTze5t1+Y+YwP05ilsa5keclaejGbUdUvP0JECGtd0ygl4p98LcYLNBNmMdrxQQ8Sipk4J8p3\ncJCm9GOF9b4Pq6NCPvrDN/hH0muacpABLAPw5G+VCwKEfc2PkLc9Qk/KWEE/SjOYIP2Tr7vuD0gP\n4X4eZg4j+ZhLmMyjMYlCBrOIEsYktFI14wCX8QEFlCeszQwSlg3f4W9Jr/EzIZ9fFFLq6z51yOMe\nYQdtmWYEz82EbN0FRlNCTc5YwYIfZTJfGBNlk0TplJ7hdiAS8jGZCICs8iN0Zz3X8AY9KYt7zX6a\nV/iOuUzxc5Jskrc9QrYrq2b0g2yYzShHlLvHmESQEENZEJNeCYhxuks30l07tqcVwrA92xhDSUrW\nqHR2qeUbYQKcy252c65vdeRdj9CO7Un/sbkiRNAhBPOYfbV5EIutjJSX8YEjjme6kRXSSZoBkXFz\nvL/VHlpZdvjKHsViDb1j8lV4Td71CF4FdfWC6OQhbgQJ0Yij/Iuv8Xt+xOW8Twlj6MImX2MadWBr\n3CHX3/kWpUaEuKqCn4tpkIc9gpdkY5dPNjSzp46y51juxgaacIgrrXgGqZPKXCIeO2lDmADFFFU5\nEQC+P1Pe9QhespAhaQfNChNIyQ06TIDxfGi9DxKiMUc4TGNu4E8cpGnaPUI6QnDbuFPRm+4rM1W6\nR7Cnc02VVD9MCxgasztqIEusSXqmwyK3vRNu2Nc0dtCWF7k5o/o0Eaq0EIC0Q4inM8G9jA9ijg1g\nWVb5wLqluHdhHB8BkTnVS9yUcX2aCFVeCB9wWVrOcaliriG4zUOyyenWLYXh0fd41XqdD0HEqgJV\nXggAHzLeNU2qPbNmupjOeH4E7k2WwM++XyGbCbbmS6qFEABXIYxidsyxVB3johOOQ+Tb2YvoE2OZ\nnvK17diedX2aaiKEphzkDqYkvS4dP5zoaGxLuNBTq80I5sY9l45QNKlRpc2nJrfzTMLzZ6mR0eb1\n4zSgAccdJtdlDKCDI0VEZlzKDGpzCtDzgFxQ5YVgt/WbeOWB+VvuiTm2lIFxs0Qmo5wCGnDccrUw\n5yFd2ESYAJvowkjmOO7RIvGGKi2EluxlP80dm+3zlZe50XIvGct0x4S5C5usgFbn81/HfeuNzDGa\n7MgbIQQJWW4KO2jLOnpYkecy/dY7TgOChHiJmywhzOBSz9rsFe9zucPHajpj2UhXFML/8BcAOrOZ\n/8f/xdx7Iy/zMjfmrK1VlbwRQogg6+nOjbxMW3bQlh3W0CBIiD20SnuPQQOOU5/PrTnCF9TNSbrU\nLXSiE1sSXvMOV3Ga2hygqesiXqohXdqxnZt4mZe0GLIir6xGZvfvRiv2JI3RE020i8VJ6mTctnRY\nxoCk1+ymNWvpyU7aJEw2XkwRr3A9YQL8h68AsfOCtmynB+uya3Q1J296BIgE0N1Ge8sCEyRkxeiB\nxHuE3YiOWZoP6VLBjNyW+iaTTXRhk5FnbhajgUi+gx/yBwDe4Fory6QmM/JKCNG0Zpdj+2G6TnTH\njXxgy+nviLTsN8lcLAopzXrH1VY6UEwR3VmvJ8wekFdDo2jK6Ol435K99GF1WmWECOZUBBAJk5KM\nH/OCJzF6tAi8Ia+FsMslxmW2iaVzwTIGpCQ+P8OTaNIjr4Wwk9YUU8QL/Ng6Vtn339pp41PUNk36\n5JUQ2rGdAGEGspQiigkQphdr2c25/J1vxQyV8plnuD1pr7BT9wh5g1REGmQRUbi4OZhCiI7MsIbe\nvMXVuWqep9gXCqMpoxdhxvgS5rx6U4xSKq34mnljNWrPNn7E7633dlv5bs4lQJg+rOZ5bq2I5mVM\niCB9WU1zW+ILk56sJeL4rXiDa3LfOI1FSkIQkZrAYmCrUupyEWkO/A3oDGwCvq2UOmhcOxn4AXAG\n+KlSKqnP8Pf4K12NtKrbaUcJY1hHD8c1ZhK523iO57gtpYerKPqznABhK9pdokDGiaLTaXJHqnOE\n24FSvvyPTgJmKKW6AzON94hIIfAdoBAYD0wVkYR1fI9XHSJ4mRtjRACRWD0QmSxfEJV6NN9YTn+b\nCCBMMGnszmt53e9maRKQVAgi0gH4KvB7vsycfQXwivH6FbCClE4AXldKnVJKbQI2QPzd86OY7Ui0\nl2hjSyGlQGTMvZK+KecByAdW0jdpL9ad9b5s+9SkRipDo6eBe8CxStRaKWXa/naBZf5oB47Y6Vsh\nfiLci5lpvd5GeyuTYjSt2ENv1ljv7bl27SlJKzv5FOWvupFQCCLydWC3UmqZiATdrlFKqYgVKC5x\nz5kT4q10ZIOLt2U/VjCCuZzLbsfx17iO63iNfqyoUkIIEK7UwXorM8l6hBHAFSLyVaAe0FhE/gLs\nEpE2SqmdItIWrE/qNqCj7f4OxjEXQo60fW6soB8r6Edv1lhJNipiR1aAMPMYnpX3aipzgALKPc+x\nXD3YZPxkTsrrCCISAO42rEa/AfYppR4XkUlAU6XUJGOy/BqReUF74GOgm4qqJN46QjzMoZC5xbIF\n+/kur/K//DTlMtLF/ECadZuECLKAoZygXspl/ITnXM2n8aiMyTzyC//XEcwP9GPAmyLyQwzzKYBS\nqlRE3iRiYToNTIwWQbq0ZjdhAo58B/to7psIHuQhy5/JLT1rPU5wO89QjxOECLp+e9/BFMsDtTn7\neZbbuI3nHEm0E9GBrZ6EhdGkTl6tLJu0ZQc7aMuVTMsqa46XjGAulzIDgMe5l2HMd6wY2/dQmKyg\nH9O4kr6s4qoUEo7b0b1CNlTilWU7kcjShymjp/Vhs3/zmscO0jRnk+W5jEAhjDXSy37EOAazyNr8\n49Z7TONKLmBlSgn+ounIZ3zmmG5p/CQvhRAiyE28RBt2Wsd6s4Y19CZMAEFZYggSsjbf+MFI5jCH\nkcCXoRiHMZ9hzCdEMG52e3NBrSsb2U9zmrM/5To30lWLIMfk5dDIxPyQtWKPtaAGkY3vq+hLMw7w\nU/4XiMQT2kZ7z/MG38vjrvnSEjGFOxwRsa/jtbRilE5nbN5sK62cpD80ymshpEJTDhIgTH+W+zau\nts8PkpEs0chEptKKPQnL2Eg3XuW76TRR4yB9IeTVfoRMOEhT3mWCr5PLeMMfN5L1SKnER93I+SnX\np/GGSi+EXFCbUyktphVTlHS/sqCSpr41J+Sa3KGFkAKnqJ3Sam+AMAHCCaNYtGRv0qERRHoht7D1\nGn/QQkiBR/gFDTnGY0wiTMBaVV7EYEdIliAhmnCIQzQhSIgiimnGAevbvSV7PUmErvGevDSfVhQB\nwgQJWWHig4T4NfdxmlocoyF38jQg1OWENSmOdsEwMdc4OvGZ1Zu0ZC+lFLKG3g5v2nhczEw+oyOb\n6ezlY2pc0D2CjTABFMIcRiIoHuZ+hjGf2pxiLiN4jElM4Q6HZSjeJH05/fmU86z1gL6soiZnAFIS\ngcn3+TNBPV/wHd0jRBE2Rvom0XnWTlDX6gncFtTsVqHoVe819M5oEhwgxGY68SnnpX2vJjW0EKJI\nNCkexnyChKnLCT5kPKvo61jlXsrAuFYjMx5TOqZYk0100SLwmUq/oJZLfsZTnMMR6/1xGqAQnuBu\nLmQJS7gwrfJ6UUZrw41kPd3ZTjtH+JdPOY9ZjEYhVhBgTSpUEae7fOUpfhZzrCHHCBLiKOekXd5a\nerLWJWjZRrpSh5OECLKFThm1VZMeukfIQ7qxkS+oox3vMkb3CFUCt/3bGn+pwubTTboOXUfKaCHo\nOnQdVGkhaDSpo4Wg0VChViONxj8qxQ41jSbf0EMjjQYtBI0GqAAhiMh4ESkTkXIRudfDcjeJyEoR\nWSYiC41jzUVkhoisF5HpIpJWnlkR+aOI7BKRVbZjccsUkcnGc5WJyNgMy/+liGw1nmOZiFyWafnG\nPR1F5BMRWSMiq0Xkpz48R7w6PHsWEaknIgtEZLlRxy89fQ6lVM5+iET73QB0AWoDy4FeHpX9KdA8\n6thvgJ8br+8FHkuzzNHAAGBVsjKJJEdZbjxXF+M5a2RQfhHwM5dr0y7fuK8N0N943QhYB/Ty+Dni\n1eH1szQwftcikn5gqFfPkeseYQiwQSm1SSl1CniDSHIRr4i2FMRLaJISSqlZwIEUy0wrSUqC8iH2\nOTIq36hjp1JqufH6KLCWSIBmL58jXh1eP8tx42UdIh9w5dVz5FoI7YHPbO8TJhJJEwV8LCKLReRG\n41i8hCbZkChJylbbddk8209EZIWI/MHW1Wddvoh0IdIDLcCn57DVYSaM8exZRKSGiCw32jtdKbXQ\nq+fItRD8tNWOVEoNAC4DbhWR0Y6KI/2lp/WnUGYm9b0AnAf0B3YAT3pRvog0At4GbldKHbGf8+o5\njDreMuo4isfPopQ6q5TqTyTvxlAR6RN1PuPnyLUQohOJdMSp2oxRSu0wfu8B/kGkG9wlIm0AohKa\nZEO8MtNIkhIfpdRuZUAkb53ZnWdcvojUJiKCvyilpvnxHLY6XjXr8ONZjHIPAZ8A47x6jlwLYTFQ\nICJdRKQOkQyc72VbqIg0EJFzjNcNgbHAKqPs643LrgemuZeQFvHKfA+4RkTqiMh5QAGwMN3CjX+m\nyTfASiGaUfkiIsAfgFKl1BQ/niNeHV4+i4i0NIdWIlIfuJTIXMSb50jHiuLFD5Ghyzoik5fJHpV5\nHhELwXJgtVku0JxI1p71wHQimX3SKfd1YDtwksjc5oZEZQK/MJ6rDBiXQfk/AP4PWAmsMP6prTMt\n37hnFHDW+NssM37Ge/wcbnVc5uWzAH2BpUZZq4D7k/2P06lDu1hoNOiVZY0G0ELQaAAtBI0G0ELQ\naAAtBI0G0ELQaAAtBI0G0ELQaAD4/9U5yXX65Y7QAAAAAElFTkSuQmCC\n",
"text": [
"<matplotlib.figure.Figure at 0x1130e4790>"
]
}
],
"prompt_number": 96
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"imshow(((loci == 1) | (loci == 4)) - clim_mask)"
],
"language": "python",
"outputs": [
{
"output_type": "pyout",
"prompt_number": 99,
"text": [
"<matplotlib.image.AxesImage at 0x1137e01d0>"
]
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAMIAAAEACAYAAADhtqk6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAC51JREFUeJzt3V+IXOd9h/Hna8nCSd1WCMNKlgQrqAxWUpBv5NI0xIVG\nlkMruTeWAwWRmN44bUxDWkulVNadE2ibK/emTlFdolQ0RKjQUCnGMrmxVRNtLGctS4IueJVoVUoI\nMaFUwr9ezKt4tN3VzozOrLbW84FBZ86e877noHn2zI7+nFQV0p3urtt9ANJKYAgShiABhiABhiAB\nhiABYwohya4k55JcSPLsOOaQupSu/xwhySrgHeB3gEvAvwOfraq3O51I6tA4rgg7gItVNVNVV4Fv\nAnvGMI/UmXGEsBF4t+/5bFsnrVjjCMG/s6H/d1aPYcxLwOa+55vpXRV+IYmxaKyqKsNsP44Q3gC2\nJpkEfgTsBT77fzc7OIap+50CHnGOO3KOQ0Pv2XkIVXUtyR8B/wasAl70EyOtdOO4IlBV3wG+M46x\npXH4EP/J8qRzOMfADME5nIMPdQjS4AxBwhAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAk\nwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAkwBAk\nwBAkwBAkwBAkwBAkwBAkwBAkwBAkYIAQknw9yVySs33r1iU5meR8khNJ1vZ97UCSC0nOJdk5rgOX\nujTIFeHvgV3z1u0HTlbVA8DL7TlJtgF7gW1tnxeSeNXRirfki7Sqvgf8ZN7q3cDhtnwYeLwt7wGO\nVNXVqpoBLgI7ujlUaXxG/W49UVVzbXkOmGjL9wOzfdvNAhtHnENaNrf8tqWqCqibbXKrc0jjtnrE\n/eaSrK+qy0k2AFfa+kvA5r7tNrV1CzjVtzzZHtIoZtpjdKOGcBzYB3yl/Xqsb/03kvw1vbdEW4HT\nCw/xyIhTS/NNcuM30leHHmHJEJIcAT4F3JfkXeAvgeeBo0meopfiEwBVNZ3kKDANXAOebm+dpBUt\nt+N1mqTg4LLPqzvFIaoqw+zhZ/wShiABhiABhiABhiABhiABhiABhiABhiABhiABhiABhiABhiAB\nhiABhiABhiABhiABhiABhiABhiABhiABhiABhiABhiABhiABhiABhiABhiABhiABhiABhiABhiAB\nhiABhiABhiABhiABhiABhiABhiABhiABhiABA4SQZHOSV5L8MMlbSb7Y1q9LcjLJ+SQnkqzt2+dA\nkgtJziXZOc4TkLowyBXhKvAnVfUx4DeALyR5ENgPnKyqB4CX23OSbAP2AtuAXcALSbzyaEVb8gVa\nVZeraqotvwe8DWwEdgOH22aHgcfb8h7gSFVdraoZ4CKwo+Pjljo11HfqJJPAQ8DrwERVzbUvzQET\nbfl+YLZvt1l64Ugr1upBN0xyL/At4Jmq+lmSX3ytqipJ3WT3Bb52qm95sj2kUcy0x+gGCiHJ3fQi\neKmqjrXVc0nWV9XlJBuAK239JWBz3+6b2rp5HhnxkKX5JrnxG+mrQ48wyKdGAV4Epqvqa31fOg7s\na8v7gGN9659MsibJFmArcHroI5OW0SBXhE8AfwC8meRMW3cAeB44muQpetelJwCqajrJUWAauAY8\nXVU3e9sk3Xa5Ha/R3s8TB5d9Xt0pDlFVWXq7D/j5voQhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAh\nSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAh\nSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSIAhSMASISS5J8nrSaaSvJXkubZ+XZKT\nSc4nOZFkbd8+B5JcSHIuyc4xH7/UiZuGUFX/Dfx2VW0HtgO7kjwM7AdOVtUDwMvtOUm2AXuBbcAu\n4IUkXnW04i35Iq2qn7fFNcDdQAG7gcNt/WHg8ba8BzhSVVeraga4COzo8oClcVgyhCR3JZkC5oAT\nVXUamKiqubbJHDDRlu8HZvt2nwU2dni80lisXmqDqnof2J7kV4FvJ/n4vK9XkrrZEAuvPtW3PNke\n0ihm2mN0S4ZwXVX9NMkrwKPAXJL1VXU5yQbgStvsErC5b7dNbd0CHhnleKUFTHLjN9JXhx5hqU+N\n7rv+iVCSjwCfBt4GjgP72mb7gGNt+TjwZJI1SbYAW4HTQx+VtMyWuiJsAA4nWUUvmn+qqn9N8hpw\nNMlT9K5JTwBU1XSSo8A0cA14uqpu9rZJWhFyO16nvZ8pDi77vLpTHKKqMswefsYvYQgSYAgSYAgS\nYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgS\nYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSYAgSMGAISVYl\nOZPkX9rzdUlOJjmf5ESStX3bHkhyIcm5JDvHdeBSlwa9IjwDTAPVnu8HTlbVA8DL7TlJtgF7gW3A\nLuCFJF51tOIt+SJNsgn4DPB3QNrq3cDhtnwYeLwt7wGOVNXVqpoBLgI7ujxgaRwG+W79N8CfAu/3\nrZuoqrm2PAdMtOX7gdm+7WaBjbd6kNK43TSEJL8LXKmqM3xwNbhBVRUfvGVacJPRD09aHquX+Ppv\nAruTfAa4B/iVJC8Bc0nWV9XlJBuAK237S8Dmvv03tXULONW3PNke0ihm2mN06X1DH2DD5FPAl6vq\n95J8FfivqvpKkv3A2qra335Y/ga9nws2At8Ffq3mTZKk4OAtHbi0uENU1YLvYBaz1BVhvusv6OeB\no0meopfiEwBVNZ3kKL1PmK4BT8+PQFqJBr4idDqpVwSN1fBXBD/jlzAECTAECTAECTAECTAECTAE\nCTAECTAECTAECTAECTAECTAECTAECTAECTAECfhQhzDjHM4xMENwDufgQx2CNDhDkLit/3hfGp9h\n//H+bQlBWml8ayRhCBJwG0JIsqvdRORCkmc7HHcmyZvthian27pFb2gy4JhfTzKX5Gzfus5ukrLI\n+M8lmW3ncSbJY6OO3/bZnOSVJD9M8laSL47hPBabo7NzSXJPkteTTLU5nuv0PKpq2R7AKnr3TJgE\n7gamgAc7Gvs/gHXz1n0V+LO2/Czw/JBjfhJ4CDi71Jj0bo4y1c5rsp3nXSOMfxD40gLbDj1+2289\nsL0t3wu8AzzY8XksNkfX5/LR9utq4DXg4a7OY7mvCDuAi1U1U1VXgW/Su7lIV+Z/UrDYDU0GUlXf\nA34y4JhD3yRlkfFh4f+Cf6SbsFTV5aqaasvvAW/T+w+auzyPxebo+lx+3hbX0HuBV1fnsdwhbATe\n7Xve5Y1ECvhukjeS/GFbt9gNTW7Fctwk5Y+T/CDJi32X+lseP8kkvSvQ64zpPPrmeK2t6uxcktyV\nZKod74mqOt3VeSx3COP8rPYTVfUQ8BjwhSSfvGHiWvKGJkMbYMxR5vtbYAuwHfgx8FddjJ/kXuBb\nwDNV9bMbBunoPNoc/9zmeI+Oz6Wq3q+q7fTuu/Fwko/P+/rI57HcIcy/kchmbqx2ZFX14/brfwLf\npncZnEuyHmDeDU1uxWJjDnGTlMVV1ZVq6N237vrlfOTxk9xNL4KXqurYOM6jb45/vD7HOM6ljftT\n4BXg0a7OY7lDeAPYmmQyyRp6d+A8fquDJvlokl9uy78E7ATOtrH3tc32AccWHmEoi415HHgyyZok\nW4CtwOlhB2+/mdf9Pr3zGHn8JAFeBKar6mvjOI/F5ujyXJLcd/2tVZKPAJ+m97NIN+cxzKcoXTzo\nvXV5h94PLwc6GnMLvU8IpoC3ro8LrKN3157zwAl6d/YZZtwjwI+A/6H3s83nbjYm8OftvM4Bj44w\n/ueBfwDeBH7QflMnRh2/7fNb9G4EOQWcaY9dHZ/HQnM81uW5AL8OfL+NdRb4i6V+j4eZw79iIeGf\nLEuAIUiAIUiAIUiAIUiAIUiAIUiAIUgA/C+g1lWlEtiunQAAAABJRU5ErkJggg==\n",
"text": [
"<matplotlib.figure.Figure at 0x1132da690>"
]
}
],
"prompt_number": 99
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"outputs": []
}
]
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment