{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Read GMT *.grd file\n", "\n", "The GRD file is prepared as follows:\n", "\n", "```bash\n", "cd ~/data/test/FernandinaSenDT128/mintpy/geo\n", "save_gmt.py geo_velocity.h5\n", "```" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Go to directory /Users/yunjunz/data/test/FernandinaSenDT128/mintpy/geo\n" ] } ], "source": [ "%matplotlib inline\n", "import os\n", "import numpy as np\n", "from matplotlib import pyplot as plt\n", "\n", "work_dir = os.path.expanduser('~/data/test/FernandinaSenDT128/mintpy/geo')\n", "os.chdir(work_dir)\n", "print('Go to directory', work_dir)\n", "\n", "fname = os.path.join(work_dir, 'geo_velocity.grd')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Use `gdal`" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUoAAAD8CAYAAAARze3ZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO29f5Ac53nf+Xm5s9glsEsMgSV2CS6oJQmaNExFsIWTWBJ9R0e0xci6UM7JPiWlO6ZOLjlVUcW6xBVRl6qLczlX6e5iy3/kzhfZcVkV2ZFVTmSpWIpokSdGoWXSBmVIAkFQXIorYQnsggtwxV2CA2AW7/3x9oN+5sXbv+bHzszu+62amp7ut99+u6f729/nx/u+xlpLREREREQ2rut3AyIiIiIGHZEoIyIiIgoQiTIiIiKiAJEoIyIiIgoQiTIiIiKiAJEoIyIiIgrQM6I0xjxojHnBGDNvjHmkV8eJiIiI6DVML/IojTEjwPeAnwUWgb8C/q619kTXDxYRERHRY/RKUb4DmLfWft9aewn4PPBQj44VERER0VPUelTvLcAp9XsReGdW4ampKTs3N9ejpvQGz54qLuPj7Qe6346IiIju4Nlnn12x1t4U2tYrojSBdS02vjHmo8BHAW699VaOHj3ao6Z0H+Z/bm+/ZzPW20+33ZSIiIguwRjzg6xtvSLKRUDrp1ngtC5grf0M8BmAI0eObOsO5yHijeQZETE46BVR/hVwpzHmNuAV4EPA3+vRsbYkslRrJNCIiM1HT4jSWts0xnwMeAwYAX7fWvtcL4612WjX7O718SOBRkT0Dr1SlFhrvwJ8pVf1R7QiEmhERO/QM6KMGAxE/2dEROeIRLkNEdVnREQ1RKKMuIqoPiMiwohEGZGLqD4jIiJRVkK/I96DhEigEdsJkSgjuopovkdsRUSijOg5ovqMGHZEoozoG6L6jBgWRKIsieif3BxE9RkxiIhEGTEUiAQa0U9EoowYahQp/UikEd1AJMoSiGb38CIq0YhuIBJlxLZEJNCIKohEGRGhECPxESFEooyIKEBUnxEdEaUxZgFYAzaAprX2iDFmD/DHwBywAPyStfa1zprZP0T/ZEQWIoFuH3RDUf6MtXZF/X4EeMJa+yljzCPJ70904TgREUOBaL5vPfTC9H4IuD9Z/izwJJEoI7Y5ovocbnRKlBb4M2OMBf5NMrPitLX2DIC19owxZl+njYyI2KqI6nM40ClRvttaezohw68ZY06W3dGf13sQEf2TEf1AVJ+Dh46I0lp7Ovk+a4z5IvAOYNkYc3OiJm8GzmbsG+f1joiogEig/UPbRGmM2QVcZ61dS5Z/DvjfgC8DDwOfSr6/1I2GRkREhBEJtPfoRFFOA180xkg9f2St/aox5q+ALxhjPgL8EPjFzpu5+Yhmd8Sww7+HI3G2j7aJ0lr7feBtgfXngPd00qhBgP10JMuIrYW8+zmSaD5iz5wMRJKM2E6I5ns+IlFGRERkIg5j5xCJMoCoJiMiymG7KNFIlBEREV3HVkukj0QZERGxKbhKnu8HGsAUsA52CEK/kSg9RLM7IqKHuD/5rgErwDiYJ3DEWUu+x8G+tz/Ny0IkyoiIbqMGNPvdiAHFk8n3g8l3A3etxpPfCSOZL6my4Ah2Fuzbe9y+DESijIjoNpq0PlmRNK/FOilB1pNlIU249pqtAIcvYx4bddd2QpVv9l6BRqKMiOgWsp4mf70QqU8GoXJbFU8BR3AkmZjb4rPk0UD5OWAxIUm5do10s3ksWUiuYbf9npEoFaJ/MqI08p4ceZClTBYhFj19WftvFRzFkaOQ3irhc/1gsg1SotSEGYD582T7OE6NziYb1tsz3yNRRkRUQS2w3PTW6d9bleS6ha/i/I8SzHkqUGbc+z1Barr75caTbdqcJ1meuAzro5gXcOTZTOuxD+U3MxJlRERZhEhSlkO+tW6S5CKpKtpqeLJg+xIpCY7j1KUQqw4ETZFGzkWhapJdGnXfDWDGwoqBWQvHDebb+U24ruy5RERsK2jzuBb4rdFr1SjBi9WiglsUj+MU4gpOAa7irof4NSFVkqGAUBP3oiHZdx1Hkg33vfvBpcJrG4kyQfRPRgBh1eiTZWifXttmUziy2K4QdeirRPFrisJcoVVNriTr15OP9h8n5PqjkzPM/Dffzz18JMqICA0/tacMATYD+0V0F4/i/JnjwDFS8pOI+SrXukCEJB8lVaLjpD7OBlfTjJaeuz338JEoIyIEvjIsIj4hyIjNw1eBKbhu7g2YAcZhz72vpGpS/o8lXGDo8eT3seT7Sa++9aTsUv5hC4nSGPP7xpizxpjjat0eY8zXjDEvJt83qm2fNMbMG2NeMMYMWEekiIgCFJnZefvp714gKlaH43ClsePq/3T+5C1u/SKOKB8nJUaNp3Bk+jnSoA84l0bBtS2jKP+AtMOR4BHgCWvtncATyW+MMYeADwE/kezz/xhjRkoco6+I/sltDB1NbYccQ6TaD0Lz8wv9ANRWQ7Pm0n0WSP2Ud1McQRcs4Py+8r93SpTW2m8A573VDwGfTZY/C3xArf+8tfaitfZlYB43M2NExGBByHFYkGXiFz3kW5Es7weWjEv3uUetX8won4X15FsCQTlo10c5ba09A5B870vW3wKcUuUWk3XXwBjzUWPMUWPM0VdffbXNZkREFEAUg/7WKiJLhXWKbhNTVjS+DLaaL/VJ0oTzhWR5nGp5pg1g7jLMXna+zgJ0O5hjAuuCc3Zbaz9jrT1irT1y0003dbkZ5RHN7iHFeODjp/N4I9KUCs50E90mSx1ZLxto2koEqTGDi1gLOTYoDMi0YA5YGHWqtAH2nfnF2yXKZWPMzQDJ99lk/SJwQJWbBU63eYyI7Q6tAifU9wRhBaiVo78+L2Ecsket6QT9JN6tZGqH0ISb3vbD9H6o4EqZ+N9fddenQl5qu0T5ZeDhZPlh4Etq/YeMMWPGmNuAO4G/bPMYEdsVWg36SjGrfBnTOWubKLWtTi5bCZ+D1XN1J83quICOJr4HuDYEnWD98Ztao94lXmhl0oP+PfAXwF3GmEVjzEeATwE/a4x5EfjZ5DfW2ueALwAncBlP/9Bau1HcjIhtD20ut0NYIUVYlvwkUXmrmalbnPgvH7vBmdBLpMO1gSPJmlrWeD/ObK9xtc94kdkNYKwNuhA3FUeOHLFHjx7d9ONG/+SAoIxizEIzUCZEeLrrWl493cKgkO6gtKNXuBf2fPAVzj91S+t0Eln/swzrtpKusj/vvo0xz1prj4R2iz1zIvqHIpO5jCKUMr6S9MuE1gt6QSb9VnNbUSGH8DScX9zn1OHdl9P+3CHIwMBybcS/WQKRKCP6g05UZJVjFCnNMm3p5PgRvcdvJ8OnPTWaH9BZIu0TrodrK4Ft+1dGs3uTocmozAASIZO6TP15gZmyidnbQYltNTyZfN9H6ucODXAiw7FNUOm/3rZEGdFjFEWY88hS38BF88qUPW4RQjmKw0qYw9rubsAnydBEbzXSXjklEYlykCDD1+sHVf+h4k+RoaWkjMAnmJoqKybGOOEHSSu4vActpNh0m3Xb8lAmQJN3/G4jPglbB3n+an97ScTbo1eQnK5V9Ts0irKYAAJ/WeeGyR887pUJQdaPe9/+PiET1zdj/XrzVF6ZOyorKh1CUaCmHZQ1ubphhvfDlO+GIh5WF0TWpG4dnk8kyk6Rld0vV1bSEXTZLFVU5s8segiyiKUokOGjHeVWYnCBTGQdrxtEVbQ97xjDSjrdOJ5/7vL/DjprlE0Zq4AY9e4EQnx+mktIhen1oeUyqiyr3lC5snWWQVa6jY+QO8BH0bastJayqUJVy5cNGmX9f2Ux6OTiI/QfVAyA9AVPkj/Umtd2yaEswrD9ff3HVHGRXBSRV9noa1Y9ZfepcrMX+TSzymUdo4q6LSLpZmA5y3UQIsat+gR0qoTz1H2ea2ZQ8CRuODbIvmcqtH1QT7OnqJQaVGL0456gUyLrZnl/3zKE26nqCJGcvxx6ULP8mSGXRDttDEVSq+xLm8dtF+2eZ5lzHCb2yHvZl0A0vbNQ51qS3Mwbo6y52+7DXqbu0LH8B6+bD33Zm1krgirXqUitF12HrMDXVlWnw96758nkuwvnEIkyC536pNpFXlAjy/dY5BvNQ41qo31XzWvM2i+vXCiHssilUOSTLasE87bnKek8n3OVY/QCvTheP9RxO8hpn32ofDXbjigr98hpV33l1ZenQOSBzlNNITLxy+fdyH4KRYNyXbmqnHsWkeQFbPDKlD2e7q+b9SKBa18GWfWXOW7V4FIV9JLYOq3fv68GXUV3SRVvO6IsRK8nmvdJzic2X5n45X2V5aucsjdGSLX520LriwJQ/jq/bUXt0ftmvQz8Yxe1K+8YWei1FeFbAqH/v1fH9JerwHdhDDqepvW5qPKMKAzDqW4uQuRTJnjRjeNWqb9d08cvLz11mrT24AntF/Lt+eqUjHJFSiTPtzhDdpezrIe26Pih4+YRcbuo8j9txn2m0e4xqt6r/cbTavk+92X/u2pVtDuv968bY14xxhxLPu9T2wZ2Xu+2BsLoxc3rq8IQ+WT56qoco0zAQsrkDTellW8ZBVrGn1h0LA2ZjlTXIWiQTb5Z6/LalrV+s32PmyFhNtP3PuRod15vgE9baw8nn6/A8M7rvemoShTdvJnzSExmswu1R3yZWsE1A+U7hW9KhwgqK/qsv3XX0DLTSVRtY6iebiusMuQ8aBh0lflU8qmIduf1zkKc17tTVH0ouvUQ1WhVlXkBI73d//imu69s2zm/UMDKb49P4qI0/aklekFmm4Fe+0u7Wdegk2Ub6CSY8zFjzHcS0/zGZN3WmNe7yDTuNqqajHl+tnaOBfn9tP1pX5tqfagtIYe5JtK8AFWorUWmvLSl5n2LqvTN8zx/aDsouv7dCoAMC1nCliPLdonyd4A7gMPAGeA3k/UDO6/3wA7Um/cGDkVD9T5Vbu4if6eQWCNQRs9W55NfJwjVFVrnE7W/HFonuaHav1rkJuhUzef5OOW7WwSyWRH5LQj76er7tEWU1tpla+2GtfYK8Luk5vUiwz6vd17UFzbvzeurMr9dRQ+dVmNFCspXe/62rH31MdpMu7gGIX9klkLN2keWQz7VrH38OqtiM03wPP9tu8hyq7SDLUiwbRGlMeZm9fMXAImID/+83iXn0Aj657LKtYOQWgmRmHzyHvo8xarrylJsvsosGwUPmdih39q3mRdhzjPFtWr0ybEsmXRKPHnXwj9G1Xrz7rduKvuIIAr/smRe7/uBKWPMIvDPgfuNMYdxZvUC8Cvg5vU2xsi83k2GdV7vPPVYVVl2wyelydAntjzlV1YdSv6kv12fq79tlZSUyqrdvGvhq0h/yDZNzD6ha7+ongpAbw+dnz5ulq+1U3O5l+rKvxcHhfAGpR1dxLaY17st/2TVUYOKfIZVfYr+vnm/y+zjtyNPkYR+631kugdN4Jps/aknoNy5y1QYUt5XrSGC9OuVyL1uW5aPlcDvEPJeIFXqKVN3p+V62YaqdQ2gCZ7nn8yb13sAT2WAUIbcyipM/61f9spnqaAqD1gn5UNmNLSayg2vbNmRzv3gi14X8rH6vrlQPXhlitRu3rUJ7Z91vHZQhaC6GQjqJYalnRURiTILq1QbpLcbZpquR9BOnUUmtyayIkLMM9fXVV3iZ/TzFiGf1HxCDKlGf1l+5/lmtaotIsc8X2uZ69/uf17FZO626uumuV6mvk4sqgHAEDe9HDpKC6r65/o+vbwAS+gYnTysoWPodqxybUK5XvbV4DrXdmvUZdZpJS/tH8zyJWYtFynG0H6hevxrFVLxeX7XdtFJXf1SX908bpG6z1u/iWgnLUgwAM0fcFQhyzLmYF7QpezvrLpCClIrLP9cskgVWn19PsZJ+19P4IhzBZjz6pP9x7199fbQNcsjztBLxfdn5r2o2iHLXpmT/TK9+2keD6myjMOs5SE0vayPUFAh7ybMUktVUYb0fLXaxClBMZn1p5nxwfuWsqtJPVJuBlhS9Wn/3kmvLbWSHwLl9Tq89ZJkXiM89mQROYf+l1AAKFSmHVS9D/rhG+02htSHOYTcvslYJY2Ay3JRQKbKVe3EV1REyHDtQBeaFDUB6fQZHcAYD6wTJQmOHGdwXQ1mVPn1ZJuQaJEprcktj1R98tYk6dcpfkr/4cxzZ5T1X+q6OkUv/aDtHq/XGDJlOURNrY6udVuUnEExMzW5hEzadlBWkeryZQMUq7i2CuELecoYjxNcqyDlnFZI+0z7bRQFt5LULSb5SeBuXDDMD6ro/uFSxzqtgyXXaCVOn/xCJnSD8FBx/vFq3nLWf+dfS70+j2iq3AvdIquiF3evjtsu+sA6nfgnIZre5aHzAvMe5G4pjE4hZKfnw9Gms5jhopQXk99CqvIt5yvmusYKjgxFSS6SkohOixWTXAhRq0AhOE2iE7Sa0NqkhlaFKesnvPV55nzI9NbXLbRely+j5EMIuTOKUMb012W7VVdECyJRtousG65T0gw9qFUDQEKIkKpGUV3rOIJbwpGckFgNR3RSzwqpP3OFlFgFMzjVKOuElOdxQ6WsJ2WO4hTmYnK8xaS+cVLS1HP2SNu1f3hCfY/nfIQ4J7xvTary0ggRL+STaDcQ8r3mqdAi32noPshaPygYpLaUxJY2vXuKMn4e38zzkRWQ8cky5GMLLevf/hBoIVfBEo6slpIydRxBzSXL66qOGVKVKfVLnceAT+EGRD2WlHsQ+BPgl0lJFxxRzibrHsUNzb+EI9ippE4x48UfvKKOO6POSQg0dG19k79Ga9v9a5jltigys6sgz8QP/ad5JFrk2+6mlaPbkeW+yGuDvPyq9NIaMAxhk4cMRTdWnqlXdCOitvt+OB19FhNbbtQpXA99fRzxW0pgRshlCkdiMziCE2Uo7Z1Ltn0sqedw8j0L/H3VrruT9jxIapY/qPa5R7VlFUeOuu3SRiFRPYK5r8w0QaLW1bzlda+Mvo5lyLKIrIpQhuTyyuUdt9uqzVfcTa69dnrZV846KNjtUfEL0Kl/ErYwUfZ8/MmQAgk9UP7DlLUdWh9mv0yRqe0fS4hyilQZSj11HME1cCQoJqpOJAdHiuNJHSuk5rBPUBM40jtI6qsU/6X4HFeSfU4mv4WY55J14JTlvcD7gcdJifrpZNv9qi3arPahVbN+YUDrdZDgnL6GIfjbJACVpfi7gc0kwRDyXgKhe1m/aPIU8CaTZLewZYmy5+jWldM3l05nkfpDUyv4N7H0otGpQEIkkqLTVL+16prDEZsOhhwnNcHnccR6ONk2r9oo39LVU/sWV4AjOOWqyfJgsr2myospfi+pn3Mm2Xcdp0al7X+AI+UajjizRl+XoFWWj0+X1QgRRMgK8E3bqmRZ1lLwUeYYnSrdLEWbd893oqp7yELdUJOwRYlyU0Yzz7qJ9Q1eZFbLNr1dj8oD4bQXTa6ibHQe5KpXr5isfk7lBGluaJ1W5SXEJj5DSfk5qH7fgyOzJRzZ6ai1KFnxbYoy1V0ptSKU9CBJY5pQx2rg/J8P4si3nhxbiFai7jULTRM2/TRp+rmVoZeRjzxi0/7QMuhEEZYl5BABlT1uL9Rs3n7+SFMDiC1JlJuCLGd5yLTOMlU0sojXDzzIvpokNVnq7Us4wplTdfllm6Rm9QROzUk+5FSybgF4gFS5LtJqZt9DGhSSB1QHXbTilNxK3U59PpKXKcpTzv/XSN0GEkU/Rur7vAdYNTB1GZo1qJlW87ipPqHcTn09sv5bDf9/LbOP3lejinLMU6955m+n5noVxVxGJWYJhQFFmXm9Dxhjvm6Med4Y85wx5leT9XuMMV8zxryYfN+o9hnYub27Bv3gCUK+Rb3s+yP94IH+nqD1QRRf2yqtZqUoNQlyaFUoBKVThISk6mqdRJeljeIDlG6K0pf7KVIyfRRHWHUcaeHVJ3VJSo4Q9kyyrFN9pLxsP5y0fcY64puyqVKdcW3d/aGllOCfBP41zpd5dBRqzbQt0gb9rbf5191H2Ye4LImUUa9Zx83yc1dtQzvohonst99/fgYYZfIom8A/sdb+OM6L9A+T+bsfAZ6w1t4JPJH83vpze3diloRIVW6WUB6hNhUhNaH9IITUKyavvP11fqEEakLtEQUoAZc6sH4ZmstpPQtJmf/39TToAinRrQIfpjUaLXX7+Y5yDlM4k53ke67B+D3nYarBxOyr7Jk7zXh9jdH6GtRt6u98EC6s73R33BzubhPTH+Dp0cTvasPEmKVi/PV5L77QvlBMKGXNy3YIRJ+rvi9CFkk76ITU5GUr9/cC14qFAUbhJbPWnsHNtIi1ds0Y8zxuCtqHcO50gM/i3uufQM3tDbxsjJnHTT72F91ufAh9n21RKzpoNZ/F36jNI21yNr39tZKE1CSuq+21wH5aJUqbxGwWMhQ/oUTEhbyaJOlBo1CbhlULGLftXoAbUnM8UXdXe+ccU+1ZVfVpn6S0fQIYt+yZO83qSp2dExcAqO9aZcfuM1xiB01GmNy9xiRrnJo4wI8WZmACrpt6g43mCNQbcO8IrIym0fNVnB9zBWiY1peNTk3R184nOd+8zTM7/Tr0eqnTR1V1lkckofbq9aFt/SAm+e/F762T/3tkencrkAMVm2iMmQN+EngGmE5IFGvtGWPMvqTYLTgjSLBIYG5vY8xHgY8C3HrrrVXb3R9kBVb0VfRTVrTZ66sbuWGFPBZIzWW9XUxXvxdMU633k6kbXtmGt79s177FeZxvUJRrE1g0rk2SSlQjjWxLz56DpCk+0va5ZJ8ZUsV61QS+zHh9jand59jBJfZPn+YCO9mfTNi5zD4OsMwak2wwwkV2cLGxg933LXGxsYPGyo1cN3GB8YkLNFZudOY5o+mxFkmDSDMWxk3r+JplCDHL15xXviwB5ZFo1rHL1FeFrHsN/wUE6YtKXDsw0AEcjdJEaYyZAP4D8HFr7evGhKbwdkUD666ZmMda+xngM+DmzCnbjr4i1LvAJz3/ivpR5pDjX9Rg3VvfJH3AZT8xu32FVFf7iEqVeiSRXLokisKS+mZJE8khJZmlZWDarZd8xikcMc4lZfU+kjs5pbYp9TA68zojtSaTu9eZZpl9LDPGJS6wkzqrTLPMYY6xwQgnOMQGI5xmPwCXdo+xwQiTY3Bx9zobGyNsNEeozaxwsTHG5doG19U2uPL0LrivAU+PuwBP00DNJudvWl8c2kWhswZCChNvGW85T3WGyufBv4c6Jbpu19cutJUj96pkRAw4YZYiSmPMKI4k/9Ba+x+T1cvGmJsTNXkzcDZZv0if5vbeNLM7RHZZ27MeoCbX1iOm8QypepV1496+8nueNDqtj63JVQhygtSf2PQ+ElzR3RaZdl8rOMV4GGdeH8FFm4+pttZxieJLan9p/4xlfOo1pnafYz+nuZ43Ocg8TUa4xBg7uMg0Z9nHMiNssJcV3sU3+To/wwFOcYJDzLHAWfZxkTEALo647x1jlxipbXCx1qSxvpPx+8/TeHSPa8eRy9AYTU1wHfXXnzImNh1sK7M9D2VIOOuYefdpL5GnZkOjPFV1R2wyCptmnHT8t8Dz1trfUpu+DDyM6+X7MPAltf6PjDG/BexnGOf2LoJvhoVuRv3tBwhCD6YsC/Gsk5q8cO00DvJ7jvTNLGSwRKrqNAnqXipCGvq3Nt+XAud9HOeVXsW5CerJ94I6h7vV8YDrZt9g7/Q5plnmECe4ngsc4gRrTHKJMc6xl5/mG4mZXeMApzjECS6xgw/yJ7zEHRziBM/wTvZzmlMc4BI7WGOS10bqbFBjpLZBrbYBi+M0ZmHiA6+y/vmbkuvWgMaYI8txy9W0oSLfYVVSKUOiWf7EzUI/TfGi44bcWh2gm/5JKMfh7wb+B+C7xphjybr/BUeQXzDGfAT4IfCLsIXm9g5Bp3Zk+YV88y0UMNDrQ/tCqvwk/Uf7BnW0UPyJkuDdII0ki1nv9/aRfccD9YkClbbcTatqXcKRuSSCH0yOJ3Uv4UzecRifPc9du7/HTi5QZ5W9nOMg87xGnUnWuMQl6qxSY4MbWWWOBe5gnttPJiy9Dm+983lO776Jwxzji/wCcyzwEnewnKjdi4xxaWQH55f2wzhcV9ug2RxpyeMcra9xeeWGdIW8JEKo4mvMIr0832O7JJl1HxXtE2qDX2evUKXuLpJkL1BIlNbapwj7HQHek7HPbwC/0UG7BhO+cgxdvbx1+sbJ6poob1YdCdfbmjgVuURrFFkPfKEfdvFprqpt2uT0P1KXBF5Ovg5MOl/fQnLsJ1+E+p2OIJ/CjQA0hSNVuDqoRmN1D9+u3cvb3vn0VT/kd3gr+znDJGvAGod4nlMcYI6Xue/FbzknzTO4V++Pw+hpeMuPvcrUbefYO7bCN3k3dV7jWY4wQpM32cmJi4dg3sAsXFnaRWN8l2vH/CgchMurk05ZNkegNpq+JPRUFvq/0oQWshD0/6X3KYNOTXC93C7h5v2mzXq3OLbMeJSb4p/USdqhBymkRvygTyj6revzAz5ixso2HZTRx9XtaXh1Qqtp7ROEfz5CyuuXcVlfBlYvOGJcWAbuhD8FPp/s/1RS/igpmTbd9dpz5BXmeJk5XmYnFzjIS1xiBwc4xRTnGKHJBiMu4v0j4CvAf8Yt/2ecH/QJ2HX0CrefWuKn+QYHeYnDHGOMS2wwwt6xc/CkU7Cjs68zMfvq1eT18fqaI/nGGKyPXnvO+j+S6wmtLzKfJPX6vP+81wjdU92qt8r6EEJk2+12ZqDbZjcMtPt0ACHd+rQqw1sO+S5D5rn/cIrK8c0rIcSG+q19jk21LCPraLM91Ea/y15Dfa+r34zi+i5eAHYmdU+nx5RA0wSOT+9J9jtK2gcbOMN+DnCKHVziv/DTvJ1n2aDG9VzgDPsZYYMNavBG0sazyWcPjij347adg7e88Sobd/91EvhpcoJDPM8hvv+xBo31nczccpqlx26Hgy6AtNGscV39DQCu1HYkwR1S1Rwyw4sUVRl/ZGhdL5+2bqtA/37W63t1rAHGEDRxgBAyvbUK87fpmypksun1WWbcAml6TxNncgs5ySATEsjxlaLOnRTS9CPWFXAAACAASURBVEletulgkE6V4c3kx3loXA+NJkxMw4pNktebMDWaugskCJV0q9w5coEJ1qgnDtabOc03eRdHOMoak+zkAtMsc+fZRdgAbgX2Ad/Fmd8AP44jztPA63D7d5c48HNLnN09zSXGWKXOW25Z4Nwbezm7PO3asGi4NDHG3PQC3//Bj3Hd+CWuG7/EFekHrglSzllHxfPIICso46/L8l36Qbx2sRkmci8Yol/BrA6wZUzvTUOW+exvlyBKFjHqsipK3GI+r+DMXRk89zhpis8K6TQOsu+K9y0O8lBbZdnPX2vpWXQBJ+sE1wPLCcEY4HlgtDVg1OTqiOYzd32fETa4xBgvcBer1KmxwRwvc469TLPMHcw7Er0IjABjwG5gGmwTnj0Lrz+NU5Yv4sz8H8Lon8HbOcpOLjDJepKUPsaVxg4mHnwVJmBueoGzb+xLCHKEK+s7nRkeCozk+SM1sqwEQdYLMVTGL59lyod8yf0imm4dd8gk2pA1N4y+dVvMUpj+Nq3QNCHKNr/sEimRreOIZ4V0jpoFXLRZDy4hI+mIObyutgtx+w+5319cTFEhWq5XDT2PI8Z3c9UU5xCwnJrjT5OOJDQLS6/sZ6K+xsiuDTYYYQcXmeLc1XxJgLcde/GqWc3LuGSyH7kjml3usH98EWZfhL8FqVk+AW85+ir1I6tMssZtLPDmxE727X2ZF595G7PvfJHvf+knmH3oRa7f9SZrP5qg0RxJrsVoK1nq6H8ZAtTf/nZBlrWRV7e/rlPl1W2Tv922ZF2bHqAX/kmIirI6NOlp9ddirnrbQ+ZulqqBdCQgXb/0fBHTVvsiZVgyMc3FP7nKtcde937759ASjT8PvJ4sX49Tl99Klq3bVktIclHtP4NL9j7uGOjsG/uumt6n2c8Yl1hmmovsoHEnsDc5xI9wWbgN4Ifw4iknpP/OiDvi11506zmNG4ZlHO7n60yzDMDFxg7OvLGf2Xe+yBwLvO2hp1n88zu5dHEHAKMTb7rj6HQoff55xFZENiGfs17fjhrslEz6LYOK3BdDhEiUVSG+Pv8B8JWiIC/Sp8vKkGYrOMUoKlJ6uohyFP/fOOnI5JLf6JOrHoRD2qJNbT260Lgq14S01+mbpGQ5CtyQNNDAzA1OQb4fR9T34BLSDwILozAO60tTrK/UOXHuEGfYzyV2cJEdnGMvC9zGsV1v47m7b3e9ffbj1GLip7wB+O/vhOYGvHUEHtgD7IEXn8GR6hNc7eK4g0tM7T7Hu3Z9k+Vz05zYOMQqdajB3NgCG82aG0gD3KhCfgJ/FtHpdVnbfOSZ592K/GaZ4z4Z16hGzt1GllruttrtMYakmdnYVLNbR6X1Hx26ir5aDKFGmkQudU/Rehzdo0YGrZVtizhzewbnyxRzW/yekAZ0mrSSoY76ar/r1bbK3fwKMIkjzGmcttub5k0ewaXi3JO0oUnqGlgBJgzMbHB56QbmOcj+vad5kp/hAR5ngTl2coE7eIlX95/lplvXXbeFk8A5mN6bfO/Didm9wG64816cL/M9uJ45bDDJGm/lu2wwwrv2fpNz7OX4E/8Vb3nPSb79g7ezZ/Ysqyt1Z3bLi8Tv861P2//fxH2RhdALUr57QQplzNes7WXcAWXr6hS+0BhQDD1RbiqyhulCLftqwidMXU7ISk/R4EejpT/1MdJAjhChRJqfxA2BdjIpO48jrQmvXm1+C3xCvhoNfh24jCPJJo4kE0UpIwHdTeuQb1Okw69JN8p1GB2/RP0nzrL2owkusJNplnmBu3gvj7GDS+zkAgvMMXXkOKYBvBXnsxzHmdmncZHw/cBtuG4OY9BwYpK7eIFz7OUEh1hh79X8Shpw7o29XDd+ifNP3cLoPa9zRTIGlmgdSESgXxiaMGWdvgd8hMgrjySLXqR52Kx9OtlvCyESZbvwo9/+fDQhc9xXLrJvaC6ZJq05ilJO6pAhzhaS5UdxKY+QDnt2kHQebT3cm6Qa+cSu27aayDlG0wK1vY5cDpIGk2SkIYCvko5EdJCrBHx55QZenb8BxmFjtsbE3jUusoO/5jD7OcMaE+znDN/Zc5EDnOLG8QbmH+DMcGm35FLeCa9PjzLS3GB5bB8H+CGTrHGCQ1xkjJ28yTeX38Wh6RP83M9/ie/wVtaXplw7Tt7g2rSQtO8o6eRqci2yVKOfwZBlUvrIU36dEJD/5Fatq4z7oB3l2Uf0KpADQ3UZBhD6ZvO7B4ZMrpBJJwpFKz8xlxeT30dwD/c9pAnhx5PyEu2ewZHjSRyZ3Zdsly6MQo4S5BFyhtaJxcTvuQrO1r2cbEh6sc6Sjkz+960br3IhKS/9wnUXSLk2SXfLnRMXWOA27uIFzjLNTt5kJxc4xa2M0OQcU9T3rDK9Z5n958+zsmeC+o8cW57bvTvxb06xd2SF/W8sMdKE07v3U2eVKVa4wE6uLOziwPQpTrOfpWdvZ/zgeRrze9y1eYDURSDf8j/IoMh+wA5v2bcMyqAT0vHJNuu+6iZC599L5Cn1AcBQE2Xf0oL8B0UTJLSmmWilqXvECLQvEZxKlADOHOn0B4dxD/o4bs7rw1w7DFsTZ4IvkQ6Lpk3jUDdGXwXPqLavALXRVNXNJvsnxLxn7jTna/uAUdfORVJlLMP+10gH8Z1wUzjsGLvEMvuY5iwnOMQITS4xxgK3UWeVVeqcZR/f2FPnLr7Hud0XeIG72M9pLibdH99y+lUau2H0FNy4ezUZ4HeM//SD/5bb3/kcF9jJOfZy+9uf49S5AzAOs7/6IovP3ZmOvakJUVwaoT74Ap2kHnK95EFeWFnmub+/H1wqu22z4LsrOkEZf2ufMdRE2RfoN21WcAe1XkdXdd6iHj9S19WkVZ0dTMo9TWswR+daiom+gCNHKTuOM8lruGi0qFHxbc6Sqko/Od03y4XwZBShe2ByZI0L9etpNPe0XgdxEUhgSs6tASO1JnVe4yzT1NjgAjt5iYMc4gTf5F0c4BR7WeHNZBCNSdZYY5J1JrnEDu7gJXbyJq9Pj3LD05dhH6wwxRqTHOMwN73lDN9/7Ceov/c1dnApvfYTsPjnd7aOgKRfaDo1KkR++mVY9aH2r6NGltIMvYyz6u43yeQdfxDa1wVEoqwKMRFEfejvJtdO/aAfRFk/ocpKDxxJ9JYAzHEckS2R5k4ukk4PK8T4AZzCXEzqn0+O8SCOBN+P88UJwc2TkmiN1tGHpP+29BuHNBcT0vm7kzb/4A/vdoS+QNq3u6n2l2HYJHUJuNQYY6F5G9N7l7nATtaY5DT7WaXODpzSvI0F6qxykTFOcYANRjjMMUbYYIQmL3EHPzbyAuy+zOu3jzLFCmtM0mSEu3iByfeusZM3OcWt7OQCY+MXeevbn3Km+GO3u3MW18WMui5+ipRvfvo+yqLAjg9dPxQTSBnzt1t+zl6RWdl65Tlq0/zupX8SIlFWQ5bq0v5FMW/1HC2SHC6QPtGiJLWyO4kjQnDkt0Q6Go9+mCXZ/Mmk/AqOHOWGOwY8tQwPTqd1zZAGd+QGFoKWhHZRg5LDKcpTCH6edBoIaYOoMTlPIVlISTIZN/PK0i6urMJi4wZOH36DHeMXWR2vs3/sNDu4xCFOsEqdJiPU2GCEDeqsssw+Jlnju/wNDnGCDWpcPgDzIwf5Lm9lgxHey2OsUmcnF1hjkjEucu6NvYzUNvjWC/clE47hyFuyCJ7EqW15cWnzWpCn5uS7jA9Sk6/eP3ScskqsKLhU9gmvqvx64R8dYDbqZF7vXzfGvGKMOZZ83qf22Zrzemt1pc0zUReQTsOpVaMsSzk9Fe26+p7DKUBtEorpug6wnJKdRJYfIPUbCklLsOWXp10gSPqEN3C/F0kH+NU+ulm1/0zSnrnk2DM4/+d9pMEeOad7k7bUSX2cWrGKStZBo3saXGmOXE0ElykeltnH9bgeNBfYySp1zrGXlzjIKQ5wln0sMMcGIzyz+6d4iTtYYI4THOIUB3iTnVxkjIu43jj1Xatuatt5HFEukhLin9KqJvH+j5A7Rb9g8h7sLOLKgx8I1ERchFCZKoGmKuXbOUaZNmiLbMBQ5jSbuHm9v2WMmQSeNcZ8Ldn2aWvtv9KFvXm99wOPG2N+bEuMch6KZvvLs8lvHSjRV1lIRwd7JnBEJ+RzmPRB1YR7cDo1f+/Dqc+7Sbs2ruDM7F8jDax8lVZyq6tjLqn6F0iDOWKCC6ncS5rwPkvrKDu6b7iQ8RzpqEaSAzqDG9AiGVXoutoGV1Z2cXl8DOqwOl6HEbjA9VeHZTvAKZqMcAcvscYky0wzxkV+mv/CAnOcZZoV9vINfppJ1vged7GXc5xjLxcZ4/Tyfibra1xev961e4V0WtsJ4IPJOcxzbRBHrpNGXsAlazmrfBZpNihXV1b9VRRelhugSj3dNtkHVFUWNitnXu8sbMq83n2fv1uQ9Tb21YB8iwkuD4T4BP1IeUOVFXNaSFbM4Q/jHnpJ9anhXlH34QbVXcSZ4w2cH/MIqQ9T3AKLuIDP/cl6UZ+LSb0P4gjuIOmQb5JutK6OKwEomSxEpxuJK0GIdR2uPL4rIXfD5YM3cH7uBs7P7GNiapU7dr3UEog5wSEusYN9LPMzfJ1jHOYF7nKmNVMc4BQL3Mb1XGCZfW7otYt1doxf5EfzM6lbQPrAz5POPqktBPkfJlRbs/Iq/f89y/zOMlF9943A99FVMZ9D9bZDZJsZfPGP1YaPstf+SajI39683u8GPmaM+R9xOuafWGtfYyvP6w1hh3PWTa8VaNNb1o79CVofODH91tX2+0h7k8wDf5CUezIpI77CJdy/cRBHjFodHk/2EXP8KVLz+BhObR0nDcrMk+ZzTpAO93YwKXdQnav4K6WMBKIk7UbORVSmKFxR4IvAVI2JXWtXe+qcY+/VqPdeznGWaY5yBICzuGnkT3GANSavztC4xiRvvnE966uTyajmSXsWSZPL7yadh1yGrpNeT5A9+ZiG/z9nEVoe0YXcOHnHyUNVQu0U3SBTeQ4kTWsiv3g/UfrSBub1/h3gX+JGT/iXwG8C/xNbeV5vaCW9rO01wkpCR76zHPv+NolQJ8GQq/5FcA+7+Ns+SJoK9ECyXco2VR1zOKL4PRxh6mj646TdEr9KSqhztAafhFyESKWHjpC/5FFKQGhV1TGOI/IJnIthhdS/umBYWr2d2ttd3+1THLiqLE+znx1c5BQHuIOXOM3NnOIAY1xijUnWNiZdX25w406um1SFSwR+njRopYNj8n9JyhW09gPX/4+G73rJ2q4JMcusDt0PPpGGtoXq8n8XtVOjakCnE8Ic974HGKWIMjSvt7V2WW3/XdxjCpswr/dAmd15N37o5vXVRyhgAK35h2IGN3HkOIczicWkFVP5bnUMIYk68DmcIr03qeMozox+EPjXpGb0QdW+JeDRZNzJWVpnghSTWrpRyraZpD0nSYlzPWlXk1T13kOqcEWBjif7TMHiKwegOQLNGuNTr9FY38memXOcn7+bl+pv8AzvZOfEBQDWVyedv3N9J6wmvYREocsxx9VHUrAkU0BMconY65He/YBKiGzynqCQWZlFfnn75hFqWRSp2nZIrxOS1JbVEKCQKLPm9TbG3Jz4LwF+gTQzbevP691uvpf/0PmOe//hEN+amK4SjKmTzrt9N2kajx+NlqHbPkCaAnQsqfcR3KttKtnWvAwTo62RevmWaPUKaTdKCeBAGriRc5K2ioKVwI+Y4xJZ918YS0m9E+NXCb8xtQcacH7lFmjClfVdbvi2xWSmxUW4ImpQ53JK76YZUoKWYI7OSxVlOUvqu5XrXfZ/1m6VssSZ5bssQrv7Fe2zmX7JLmIz/JNQTlFmzev9d40xh3Fm9QLwK7DF5/XWKHooivxV2k8ZKit+QVF5J3GEdQ8pAQiZLpBGymfU/klXw6se46dxhHAQR5hTSb014L7RtHvkeFLf0zth/WWo3QbHX4eDNziy0QpQVOYqaRRdiEknuE+ofYSEdbqQTrcRE11IW4jLD7z46lnIUncIOE6qlmdJze8lVd89SXlRm1K/nzeb9X8Xbff/4xAphUzlUEBmM4m1VyirqgcIhUSZM6/3V3L22Zrzemv4pOdfySyfUJZvUj/8WpkJoc5xLZGI//GI2k9ShiSIIgR2P+nAv3W13+M4ApnBKcw5YH4Rnpp15vqx25J0mhtaz2VWtfs4qcl6MDmeKN05Wgnv6oAbCaQt/rWSa6oTwEPKrYl7UYjfVOqpk46sJIp7jjTSLYnwuq+8PgaESbwKyqhMQcgd4y93At2WTk3mTvYbROIugXb+/r5iYPyTkO9rzCsfIk9t7s7RGoTQydBahYp5qB/mg6TmekOtW8SRxOdxvkmp7x8k38eS7/mX3QHuI01qXyAd7OIorVNVCIlNkY5MNEvq+5P1ElFeUvtq9ab9d1pRCQHrnMYmqek/o35LYEgU9EHSqD3qW66zdmeQtEECOpocG97vItLRLzvIvkf8ffNItRvEEvK3FqHbRK3r6zDKvVlmN8SpIDpHltNdb9c3iL8sDyukKlIUjigweeC0YhLCmVJ1nCRNHJ9S+66oOmZwZCejDzWT35LIThPuvzPtYQSw+nLaI0jqqKnPLNcOECxmbE1tE7O4SRo48c/Lv0aS2yjXSpNbjTS3U8z+x5Pr0EjOa5GUCBdJXy7iI9U+SLluojSFpDXySKbmbdcvRV8l+v5pvc0nlH6oL78N7dZRtDwkGDpFOdAoUpWhB0hUolaEen9RmvoBEtNSHjZRW3M44tBvbyFKIeAHSBVYDZePWccFfI4DM3e2JsOvAlO3OTKU3EMhmXXSnE0hHfGtos5L8jvF7BUi9QM6QsxaTQoJTpEOtCu5pGI2J9Fy5kkHA1kijczXSSPxktsp0OpQpwQ1aU34F/hKMs/lklembDCwiGDaeYKLiKoKkbVb1wDnTIYwVIpyoMxuSAmnigKQVBRRUwdpzZUUE1Kg1ZY8zEIiclxJ7xGyFQLSvXkOk/oDpdxTpIP8Pk6a5qMVaz3ZV8i5TivRa3UpSlhMVSGDBXVuct4S4NHXS1+3RVL3gbwUpLz4O+VcxNwXk38KR6gfUOfwAK3jd/rHFH+uVsAhgqvyEvTdLDXvW5cJuWRCZf3lkPnuWy/+p1topy5xuQwZhoooBxK6C6IgdMOL6hsnTUWR3/KQirmcdSNpv5cmGq0shbggjfBOkE4NsYrrzbMCHF90Zf8URyxLpOk7zWSfRVoH4xXTVNRjnTQFZ0pdEzkX+S0EO0GrKS3nJOtQ+yyo82smbZtP1sm+UqcQ5UHSiLY29UnOUa6PmNlCjNqf6hOaDqxlqcQiZRciyaxyeXXnEXGn5m1ZIu2UcNeLixRhM/2TEImy+8h6ECSAIDeYTo+BdBAGbQpq+DemJlrZNsG1xLBIqnwlt1FSie6eTfMshaylTVLHLIqQL6dJ7eDSanT7dc6n9tdNkQZdat63jmhr8tTuAlGA96h9RfWJH1R6B4nPtIkjVSmvX2SyjyZC+X/0uYv61Qnr3XZWFdVXlZSqtq9Iafrkm9eWLAUbOsbqNXsPNKKPspvI81FK17mDye85UlNPD8SQpUy1z0vKC7HMqN8CMd/FFL6PVKHpXMfDpIGgpUWozzqiOYIj2flkGeDgaHqeDdLpJvyAiCZKUbsSUBJ1qglJR+6ljhlVVs5BlsW9oF8s4sutkY4Cf1jVL2rXv5667fr4EtDR1122deOp0X5k/xid1JmHKseqok6rtr1Jq594CDA0inLg/JMh6GCEfIv/UXq1CHmIWQitKivPXwWtKkibsuPefuKbq5OOMzlHOtLQLE4dggviLAH8nlsvvW+WLrhlP6ItpCbHl2MKsWk/piYbWd/wPpIorpVHw6tHFJ3fHm1+y76Sdyq/xWcbgvaT+n7Gmrde2tUuqWW9SEN+Rv93FWWZd/+E6shSfmUUYhkl2mVsttkNUVF2HzofElof3qVknQRVGtfsnW12a+d9SOn45cfVN+q3kI74+xaS+ngd+PU0xWYOOLgzJTjJ79T+ulmvTp33KcSljwutwRnta9MEJOt11NknRz9HUZ9jnWvTf7ICJSGC9IlCX/uy0eoQsgIwWe0K7ePvm/cE+/dI6PjdJLoekmO/EYmyFxCylGX94AvR+AjdvH4AQH/rB8dfrwmqodYJuejue3cn6yduSP2RD5Km3fgKUroritleU+v9NCZNklKf3wNG150FP2Ah+2Q98Dq6X4Qa6f+i26tN8ywTuZdPT+hYWWowry1F5Byqt12SDBF8qNwQEurQmN5DB0mOlpui4f0u4/cJffJUgr5JQyajtEFM1Vm1v+QbalNeotc+UdZoJVHtYxSyFFNauwK0mR56GVRRN77y9pWlPu8s+ETrXz99nInAfr2WGaEXZFUCCq3z6+mWzzVUdwi+v7gC+mF2w5AoyqHwT+bBN9/wlkMkUcbHlPdghAhZ6hTXgPYhav+hTjGS/eT3LKk69NOYRIEtkAZSIDXLNbHq9uh9dZt9ovBNYtkWug5+gCgE3cNHXxNN7kK2muT9tvVSIWk/bNE5k7G+zLo8c79dhOoZsiCOYCiIcqih+2dn3SRlFFOZ8vo7z0z0e8Do3Eg9QIQmGykrJrZOhJfy4lIQc17I0Yfu3qjblnXOWaooaxtcS+IhiE9V6tOBHa2U/S6WIVXsv8Cq+g5DrhR/e6dkVsZMz7pn8v6DIgxhgrmPSJSbAVEkfv/ekNIsujl9ZKks39T2ywsxSj6nENeMKitpObqdU6qsr7D83Ea9Tfv9ZgPn5iu/LPLUhJFFqtrPWAQZfESOKf9RyOyFa32ZcsxQ2/32ZxFVGb9niCxDxwot+8jblqVWi65lXru3AAb+NIbe7BaIsvSR5UgvUlm6rO4Z5Pslfd+Vb/bqAIuORosprAlOeuZoIslTKSFCCJX3yUiXD5FiGcVWFvJikLb5eZV5Sr2MKZylDjWKzifLH9sO4eW1r0hxFrWzneNXQL/8k1DitI0x48A3gLGk/J9Ya/+5MWYP8MekQzH8UjK5GMaYTwIfATaAf2StfawnrR825CmHIgUQUiaSj9nkWiIW8tSpSKKYtCLUy/oY2k+HKqvJI0t1ZSHkm8z7Haq36BjtBFrylF4ZhV+WEPLUWkgttnOMPFR5wVRVlWVelEOMMqdyEfib1tr1ZO6cp4wx/wn4O8AT1tpPGWMewU0u8IktPa93p1ggjTQXPYRCbroXDlyb7iPRdU0qsl0rRDH9ffhmsazzy4gv09/uHze0fxax6m2d5CdqVH048whBrrWfw1m2zrLkkWVSFx0rpOqLVGZZbCGS6wYKL4e11pK66UeTj8XN331/sv6zuKEWPsEmzes9lJgn7cbo5w5q1ScK0SdHSPsj633GuZZ8JeDim8myLD2GIExg/u+sBzmPNPP8aVmmZJXhx8r4+vJ8mX475HfIx1mWjNpB6FqWQZ469ZezjlsVvTj/IUCpS2WMGQGexT3i/7e19hljzLRMLmatPWOM2ZcULzWvd6njbhX/pOB+XK8XPYCD9v81vXW+31EP2iDkOKN++71hfIiyhNYeO3KMoruhiulVVLaMjzIE3ac7azuqTBmiDrUvTwX3kiyL1lU9bpE7oZ168nyYsn0L+Seh5CVLzObDxpg68EVjzD05xUvN622M+SjwUYBbb721TDO2BmQEHxmMVpvYTVpTKfQNJze47po4TmvOo0+yvqqTffx0l6zlohtflG+ZtKciMzRru65DUpOyjiXXJ3SMLNWW5f4o+9Ko6suTfcqSVTukk7dPleMKfB95mX3y3EpDiEo9c6y1Mprhg8CyMeZmcFPXAmeTYouUmNfbWvsZa+0Ra+2Rm266qY2mDykWz7kr1CAdUWfd+2jilA+kilAi1ZL0HRoCTMr65Oqb4z7J+nVkEYvU76c8FcH3k2YRpN8uuJaQs1wB4956qP6gl0UVNaz38f/bMscog7J1lqlHvkP3QRV0su+AoJAojTE3JUoSY8z1uLGiT+Lm7344KfYw8KVk+cvAh4wxY8aY29iK83p3go/vhfnL6YC4ojC10pSxGPXAGvpbTO4sE1ug03tkJCG9rspDnlU+1B2txrXHCLVJly+CvFD8kYwEWf7SUNvLnHcVctDR/HYVW1n0guyz4N93ndYzxChzWW8GPpv4Ka8DvmCtfdQY8xfAF4wxHwF+CPwidG9e7y3nn2zB63ByryM7XyXVaDXD/Te7lAkFPXwzG8LqSpfN+5bj5qEo+OL7WYtIMuuOlFQnKeO7FcqYqPqc2vXVhZDX8yTPjeG3pyyKzrXMf+efexk3S2g/f11WQK0D9Ns/CSVuFWvtd4CfDKw/B7wnY5+tP693J/j4Xvjtc7C0NyUOPfJOVnK69JkuIknfL+eX85dD5fN8bz7RZD0oZYkr67dGnXDbs/YvUosa/rl0SqT+da/y4vHLt4ui4xQp8Kz9i9rmv2gh7eE1xIijB/ULH9/rvnUakE88csPV1SfLv5WlAvzAjl82RKRlTfFQHVk+QVkOjQ9Z9Xhl1lepL4/ou11/1kvO/930PmVfKp36Pbvh4yx6iQ0hBrL5W9vs9iBEBq3msh+k8U1o3+fmP1xiwvvdEMv6B0PrQ6a5vy0PRYPolkEeaXTTnPaP1w55ZF3DdtvRTfTiWhUdb4gRFWW/IdFo+ZYpHDRRarKsBT6hgW8bge2dPhi+Kgydi1Ykvoqt2oaqqTeb9eCXOY6+BmVUe9W2V/E/dqPOrLrLHqvNycQGwT8JkSj7Dz+9RqfvZBFj6OObjKGUoW4hi8DyHqK8dmQ9oINi7/jtq6LGfL9uVlpQt4ivarvK7KvvM39dWcwUFxlkRKLsJz7OtZODQVgt+v5LSMkw5FfLIkn/pu8GGeURYp7qzdpnkFHFh1vF35d13SDb3ZDXlixiq0q63fI3xmBOd7Gt/JMAv0wa6RZ1qU1vuPYmH1fb9IMonKbiCgAADS1JREFUQZs615JhluLLS/sIIS9olEUMWaZ6mQe3XaLpNqqYy1WvaV59Za5RkTuiirlcxnqpCiH6xWq7DYrZDQNIlNsSHyYlrbyE7JAy87cVDXThI0vJ+vDzIPMUUEixdvKA5aWqhJY3C0Uugypt6iRSLdvKKt2sdZ0SYh5kHvghRSTKQcGDtM4GWOZtnkUi3TKn80zqTh/svGOF9m+HCDshz7xzr1JHlWtQxvfZTlvKXMdeKvKVpP4HeniMHmOYvENbHzLUiEx6Ba2k6Sfu+jd9nh9Qj2tZo1it5W3TZfx2VCGGsvCPU0Utd5KOo/2MPtl0y2Wg0a0XnD6+T/jtuAXK3AshyL06pBOKaQyUotx2/sksyKjkviLJG3i3CCHSyLv5y5q7zcD6POQp1bzzaCfvr1NzXJNKlWsn+1RZ323Idc47h7Log1tjkPyTEBXl4KLIV1TFQe+jyo1fpFCqIkudFXVz68TX2Smqvgyy0AlZtVtPkUujl+jW+Q4ABkpRRngITQ4G3SfFdszITvx2ofqz5tMZVrTjjqhSvhcBl25ji5AkDBBRRrM7AzLUWlUiyVJgIUVYxoxsJ8gQioSXMemHHdptos+/zDlWyb3sdlS9G/BHi9oiGBiijChAtybf0ujGg5ZHgO2Y61vhAdMvIE18VZVit1RjVf9q0f5F5fQcT21g0PyTMPjiPQJa54kpG9ToxFwvizIqqWoQZgv5tboaxRb4Jn07Abmy7crKU81q01b53wIoM8L5uDHmL40x3zbGPGeM+RfJ+l83xrxijDmWfN6n9vmkMWbeGPOCMea9vTyBbYP14iKZ0Gku3byZh9lE7BRlXhCdKLkst0jIldKvbIBO2jBkKHN6WfN6A3zaWvuvdOF25vV+9lR7jd920PmVWWgnlWVQMCzthHI+226+mMr6k7uJshaD/t6iKFSU1iE0r3cWrs7rba19GTeb9Ts6bmlEucTdQb1hB7Vd7aAZ+IS29+qYIfRKVealF20TkoSSwRxjzIgx5hhupsWvWWufSTZ9zBjzHWPM7xtjbkzW3QJojbhIm/N6Ryjo+WJ6DT9w1K0HYZgeKD+ZfhiUVZXc2nZ6RoU6C3T5fhzEQA6UJEpr7Ya19jBu6tl3JPN6/w5wB3AYOAP8ZlK89LzexpijxpijvPlqW43fVmj3QWwnkNLLIbF6obZ6iayeSEVlBwWhjIN2yM3vJCCfIR8+rSzantfbWrucEOgV4HdJzetFKs7rzfXbaF7vdiFT2kLvgiOhCHWvHv5BJBXornk7SCgiyTLnJmk/WR0htjDantfbGHOzKvYLwPFkOc7r3WuUTQ/qVP34xNltotjsh63oenSjF1O3UfUaZbWv7LkVdZ2tmu61RVDmlLPm9f53xpjDOLN6AfgV6N683hE5kOh3XiJzOz1pihRlL0ltMx7ATtWyvt69VNp+/bWcbe3U3c5+m/BCG1T/JICxNi+AvUmNmD5i+XtH+92M4YOOgveCZHr1cPgP/maj26k7vUBW/mRZsiybLB5aX8Uv20X0myiNMc9aa4+EtsUujMOMVdJE9Co3dBmV2Gv1CJtHkn70utcPfzfqDpm/ZQN6RUnqZfbVvsiISJRDjyaOMEOBnqz0liK/42Y8HJt1jJAC64cLIe/4WXmY7bS5KomGjtkHguy3mizCQBDl2w8M/oUaCqyS3ugrpCQq8BVV6APhkW+6jW7V6aus0LlstirKuq55Jm1ZtdgptsOgJD3AQMWvssgyDsFWAX6fcJ3rJt/j6lsejHVaA0SCXhBNlcE9spBFMpupHn2UzSTI+t1PkowEmYuBIsoshAg0kmcJ+GNZ4v3216/SSqJSxie2bgVhOqlns/I9y6Lfx89CiIj9YeD6jGGwJoeCKEOI6rNH0CQ6QfaoRTVSX9YUKenJLJJNtSzlQwM75EVksx7iIsLeSkO1VUHonPMGSdmO16hNDC1RZiESaBeRN7SbkFQtUE7/9s19n9xCaSlawYYeZn3sEFl2iwDy6ilqYy9R1v/ay2uzzbDliDIL0XzvM/RUFtrklztQT7VbRvH4ZKnr2gxspiorm8KV56ONBNkRtg1RhhDVZx+QZQZqFarnMm+nrm5gUIglS1FnmdniDgm5Obb1094Z4qULIBJonyGKcyL57tVd2gvTNERMVYJWWcGXJu4FUuQO0cf1lwcQwxDIgYG/jIOFaL5vMtZxd6gQRFmlWYSsCH63VWS7qT86IKb37WQ6kIiOEImyQ0T12WM0Cac5TXjl9ARseXUVLVdFEdHm5XjqgJDvu90m4zwOCyJR9ghFJkUk0g4RUldFZNluhDpEdn7vn3GKyS0vTWdQfKKbiGExuyESZd8QlWiPoE30duGrTa32soIjUQFuaUSiHDBEP2gXoElLfH1FuZv+8oD1XonoL0oTZTJw71HgFWvt+40xe4A/BuZwA/f+krX2taTsJ4GPABvAP7LWPtbldm8rRPXZAXwl6I9FmRWRjuQYoVBFUf4q8DxwQ/L7EeAJa+2njDGPJL8/0c683hHtIRJoGwgN+hGx6Rgm/ySUvE2MMbPAzwO/AfzjZPVDwP3J8mdxk459AjWvN/CyMWYeN/HYX3St1RG5iOZ7RER3UfZ9+tvAPwUm1bppa+0ZAGvtGWPMvmT9LcDTqtwicV7vviNPfdpPRyKNiMhDIVEaY94PnLXWPmuMub9EnaXn9QY+CnDrrbeWqDaiFxAC9Yk0EmdERIoyivLdwN82xrwPF0O8wRjzOWDZGHNzoiZvBs4m5RcpOa838BmAI0eO9H+Gs4gWRPM9olcYNv8klCBKa+0ngU8CJIry16y1HzbG/F/Aw8Cnku8vJbt8GfgjY8xv4YI5cV7vLYIYPIrYrugk5vcp4AvGmI8APwR+EeK83tsRUX1GbHUMxLzeR44csUePxnm9twMigW5vDLLZnTevd8wii9hU5D0okUQjBhWRKCMGBjHyHjGoiEQZMbCIwaOIQUEkyoihQxzCbjgxyP7JIkSijNhyiEo0otuIRBmxbRAJNKJdRKKM2PaIeaC9xzCb3RCJMiIiiKg+IzQiUUZEVEAk0O2JSJQREV1ANN+3NiJRRkT0CLEXksOw+ychEmVERF8QFehwIRJlRMSAIJLn4CISZUTEAGPYg0dbweyGSJQREUOJ6P/cXESijIjYYogmfPcRiTIiYhsgkmdniEQZEbFNMez+z81EJMqIiIgWdItAt0ogBwZkzhxjzKvAG8BKv9vSAaaI7e83hv0cYvv7i7dYa28KbRgIogQwxhzNmthnGBDb338M+znE9g8urut3AyIiIiIGHZEoIyIiIgowSET5mX43oEPE9vcfw34Osf0DioHxUUZEREQMKgZJUUZEREQMJPpOlMaYB40xLxhj5o0xj/S7PSEYY37fGHPWGHNcrdtjjPmaMebF5PtGte2Tyfm8YIx5b39ancIYc8AY83VjzPPGmOeMMb+arB+mcxg3xvylMebbyTn8i2T90JwDgDFmxBjz18aYR5PfQ9N+Y8yCMea7xphjxpijybqhaX9HsNb27QOMAC8BtwM7gG8Dh/rZpox2/tfATwHH1br/E3gkWX4E+D+S5UPJeYwBtyXnN9Ln9t8M/FSyPAl8L2nnMJ2DASaS5VHgGeDeYTqHpF3/GPgj4NEhvI8WgClv3dC0v5NPvxXlO4B5a+33rbWXgM8DD/W5TdfAWvsN4Ly3+iHgs8nyZ4EPqPWft9ZetNa+DMzjzrNvsNaesdZ+K1leA54HbmG4zsFaa9eTn6PJxzJE52CMmQV+Hvg9tXpo2p+BYW9/KfSbKG8BTqnfi8m6YcC0tfYMOCIC9iXrB/qcjDFzwE/iFNlQnUNith4DzgJfs9YO2zn8NvBPgStq3TC13wJ/Zox51hjz0WTdMLW/bfS7r7cJrBv2MPzAnpMxZgL4D8DHrbWvGxNqqisaWNf3c7DWbgCHjTF14IvGmHtyig/UORhj3g+ctdY+a4y5v8wugXX9/g/eba09bYzZB3zNGHMyp+wgtr9t9FtRLgIH1O9Z4HSf2lIVy8aYmwGS77PJ+oE8J2PMKI4k/9Ba+x+T1UN1DgJr7SrwJPAgw3MO7wb+tjFmAedi+pvGmM8xPO3HWns6+T4LfBFnSg9N+ztBv4nyr4A7jTG3GWN2AB8CvtznNpXFl4GHk+WHgS+p9R8yxowZY24D7gT+sg/tuwrjpOO/BZ631v6W2jRM53BToiQxxlwPPACcZEjOwVr7SWvtrLV2Dnef/3/W2g8zJO03xuwyxkzKMvBzwHGGpP0do9/RJOB9uCjsS8A/63d7Mtr474EzwGXcm/IjwF7gCeDF5HuPKv/PkvN5AfhbA9D++3Bmz3eAY8nnfUN2Dn8D+OvkHI4D/2uyfmjOQbXrftKo91C0H5eZ8u3k85w8q8PS/k4/sWdORERERAH6bXpHREREDDwiUUZEREQUIBJlRERERAEiUUZEREQUIBJlRERERAEiUUZEREQUIBJlRERERAEiUUZEREQU4P8H9kwL1UWLWmgAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# read\n", "import gdal\n", "ds = gdal.Open(fname, gdal.GA_ReadOnly)\n", "data = np.array(ds.ReadAsArray(), dtype=np.float32)\n", "\n", "# plot\n", "plt.figure()\n", "plt.imshow(data, cmap='jet')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Use `scipy.io.netcdf_file`" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUoAAAD8CAYAAAARze3ZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO29f5Ac53nf+Xm5s9glsEsMgSV2CS6oJQmaNExFsIWTWBJ9R0e0xci6UM7JPiWlO6ZOLjlVUcW6xBVRl6qLczlX6e5iy3/kzhfZcVkV2ZFVTmSpWIpokSdGoWXSBmVIAkFQXIorYQnsggtwxV2CA2AW7/3x9oN+5sXbv+bHzszu+62amp7ut99+u6f729/nx/u+xlpLREREREQ2rut3AyIiIiIGHZEoIyIiIgoQiTIiIiKiAJEoIyIiIgoQiTIiIiKiAJEoIyIiIgrQM6I0xjxojHnBGDNvjHmkV8eJiIiI6DVML/IojTEjwPeAnwUWgb8C/q619kTXDxYRERHRY/RKUb4DmLfWft9aewn4PPBQj44VERER0VPUelTvLcAp9XsReGdW4ampKTs3N9ejpvQGz54qLuPj7Qe6346IiIju4Nlnn12x1t4U2tYrojSBdS02vjHmo8BHAW699VaOHj3ao6Z0H+Z/bm+/ZzPW20+33ZSIiIguwRjzg6xtvSLKRUDrp1ngtC5grf0M8BmAI0eObOsO5yHijeQZETE46BVR/hVwpzHmNuAV4EPA3+vRsbYkslRrJNCIiM1HT4jSWts0xnwMeAwYAX7fWvtcL4612WjX7O718SOBRkT0Dr1SlFhrvwJ8pVf1R7QiEmhERO/QM6KMGAxE/2dEROeIRLkNEdVnREQ1RKKMuIqoPiMiwohEGZGLqD4jIiJRVkK/I96DhEigEdsJkSgjuopovkdsRUSijOg5ovqMGHZEoozoG6L6jBgWRKIsieif3BxE9RkxiIhEGTEUiAQa0U9EoowYahQp/UikEd1AJMoSiGb38CIq0YhuIBJlxLZEJNCIKohEGRGhECPxESFEooyIKEBUnxEdEaUxZgFYAzaAprX2iDFmD/DHwBywAPyStfa1zprZP0T/ZEQWIoFuH3RDUf6MtXZF/X4EeMJa+yljzCPJ70904TgREUOBaL5vPfTC9H4IuD9Z/izwJJEoI7Y5ovocbnRKlBb4M2OMBf5NMrPitLX2DIC19owxZl+njYyI2KqI6nM40ClRvttaezohw68ZY06W3dGf13sQEf2TEf1AVJ+Dh46I0lp7Ovk+a4z5IvAOYNkYc3OiJm8GzmbsG+f1joiogEig/UPbRGmM2QVcZ61dS5Z/DvjfgC8DDwOfSr6/1I2GRkREhBEJtPfoRFFOA180xkg9f2St/aox5q+ALxhjPgL8EPjFzpu5+Yhmd8Sww7+HI3G2j7aJ0lr7feBtgfXngPd00qhBgP10JMuIrYW8+zmSaD5iz5wMRJKM2E6I5ns+IlFGRERkIg5j5xCJMoCoJiMiymG7KNFIlBEREV3HVkukj0QZERGxKbhKnu8HGsAUsA52CEK/kSg9RLM7IqKHuD/5rgErwDiYJ3DEWUu+x8G+tz/Ny0IkyoiIbqMGNPvdiAHFk8n3g8l3A3etxpPfCSOZL6my4Ah2Fuzbe9y+DESijIjoNpq0PlmRNK/FOilB1pNlIU249pqtAIcvYx4bddd2QpVv9l6BRqKMiOgWsp4mf70QqU8GoXJbFU8BR3AkmZjb4rPk0UD5OWAxIUm5do10s3ksWUiuYbf9npEoFaJ/MqI08p4ceZClTBYhFj19WftvFRzFkaOQ3irhc/1gsg1SotSEGYD582T7OE6NziYb1tsz3yNRRkRUQS2w3PTW6d9bleS6ha/i/I8SzHkqUGbc+z1Barr75caTbdqcJ1meuAzro5gXcOTZTOuxD+U3MxJlRERZhEhSlkO+tW6S5CKpKtpqeLJg+xIpCY7j1KUQqw4ETZFGzkWhapJdGnXfDWDGwoqBWQvHDebb+U24ruy5RERsK2jzuBb4rdFr1SjBi9WiglsUj+MU4gpOAa7irof4NSFVkqGAUBP3oiHZdx1Hkg33vfvBpcJrG4kyQfRPRgBh1eiTZWifXttmUziy2K4QdeirRPFrisJcoVVNriTr15OP9h8n5PqjkzPM/Dffzz18JMqICA0/tacMATYD+0V0F4/i/JnjwDFS8pOI+SrXukCEJB8lVaLjpD7OBlfTjJaeuz338JEoIyIEvjIsIj4hyIjNw1eBKbhu7g2YAcZhz72vpGpS/o8lXGDo8eT3seT7Sa++9aTsUv5hC4nSGPP7xpizxpjjat0eY8zXjDEvJt83qm2fNMbMG2NeMMYMWEekiIgCFJnZefvp714gKlaH43ClsePq/3T+5C1u/SKOKB8nJUaNp3Bk+jnSoA84l0bBtS2jKP+AtMOR4BHgCWvtncATyW+MMYeADwE/kezz/xhjRkoco6+I/sltDB1NbYccQ6TaD0Lz8wv9ANRWQ7Pm0n0WSP2Ud1McQRcs4Py+8r93SpTW2m8A573VDwGfTZY/C3xArf+8tfaitfZlYB43M2NExGBByHFYkGXiFz3kW5Es7weWjEv3uUetX8won4X15FsCQTlo10c5ba09A5B870vW3wKcUuUWk3XXwBjzUWPMUWPM0VdffbXNZkREFEAUg/7WKiJLhXWKbhNTVjS+DLaaL/VJ0oTzhWR5nGp5pg1g7jLMXna+zgJ0O5hjAuuCc3Zbaz9jrT1irT1y0003dbkZ5RHN7iHFeODjp/N4I9KUCs50E90mSx1ZLxto2koEqTGDi1gLOTYoDMi0YA5YGHWqtAH2nfnF2yXKZWPMzQDJ99lk/SJwQJWbBU63eYyI7Q6tAifU9wRhBaiVo78+L2Ecsket6QT9JN6tZGqH0ISb3vbD9H6o4EqZ+N9fddenQl5qu0T5ZeDhZPlh4Etq/YeMMWPGmNuAO4G/bPMYEdsVWg36SjGrfBnTOWubKLWtTi5bCZ+D1XN1J83quICOJr4HuDYEnWD98Ztao94lXmhl0oP+PfAXwF3GmEVjzEeATwE/a4x5EfjZ5DfW2ueALwAncBlP/9Bau1HcjIhtD20ut0NYIUVYlvwkUXmrmalbnPgvH7vBmdBLpMO1gSPJmlrWeD/ObK9xtc94kdkNYKwNuhA3FUeOHLFHjx7d9ONG/+SAoIxizEIzUCZEeLrrWl493cKgkO6gtKNXuBf2fPAVzj91S+t0Eln/swzrtpKusj/vvo0xz1prj4R2iz1zIvqHIpO5jCKUMr6S9MuE1gt6QSb9VnNbUSGH8DScX9zn1OHdl9P+3CHIwMBybcS/WQKRKCP6g05UZJVjFCnNMm3p5PgRvcdvJ8OnPTWaH9BZIu0TrodrK4Ft+1dGs3uTocmozAASIZO6TP15gZmyidnbQYltNTyZfN9H6ucODXAiw7FNUOm/3rZEGdFjFEWY88hS38BF88qUPW4RQjmKw0qYw9rubsAnydBEbzXSXjklEYlykCDD1+sHVf+h4k+RoaWkjMAnmJoqKybGOOEHSSu4vActpNh0m3Xb8lAmQJN3/G4jPglbB3n+an97ScTbo1eQnK5V9Ts0irKYAAJ/WeeGyR887pUJQdaPe9/+PiET1zdj/XrzVF6ZOyorKh1CUaCmHZQ1ubphhvfDlO+GIh5WF0TWpG4dnk8kyk6Rld0vV1bSEXTZLFVU5s8segiyiKUokOGjHeVWYnCBTGQdrxtEVbQ97xjDSjrdOJ5/7vL/DjprlE0Zq4AY9e4EQnx+mktIhen1oeUyqiyr3lC5snWWQVa6jY+QO8BH0bastJayqUJVy5cNGmX9f2Ux6OTiI/QfVAyA9AVPkj/Umtd2yaEswrD9ff3HVHGRXBSRV9noa1Y9ZfepcrMX+TSzymUdo4q6LSLpZmA5y3UQIsat+gR0qoTz1H2ea2ZQ8CRuODbIvmcqtH1QT7OnqJQaVGL0456gUyLrZnl/3zKE26nqCJGcvxx6ULP8mSGXRDttDEVSq+xLm8dtF+2eZ5lzHCb2yHvZl0A0vbNQ51qS3Mwbo6y52+7DXqbu0LH8B6+bD33Zm1krgirXqUitF12HrMDXVlWnw96758nkuwvnEIkyC536pNpFXlAjy/dY5BvNQ41qo31XzWvM2i+vXCiHssilUOSTLasE87bnKek8n3OVY/QCvTheP9RxO8hpn32ofDXbjigr98hpV33l1ZenQOSBzlNNITLxy+fdyH4KRYNyXbmqnHsWkeQFbPDKlD2e7q+b9SKBa18GWfWXOW7V4FIV9JLYOq3fv68GXUV3SRVvO6IsRK8nmvdJzic2X5n45X2V5aucsjdGSLX520LriwJQ/jq/bUXt0ftmvQz8Yxe1K+8YWei1FeFbAqH/v1fH9JerwHdhDDqepvW5qPKMKAzDqW4uQuRTJnjRjeNWqb9d08cvLz11mrT24AntF/Lt+eqUjHJFSiTPtzhDdpezrIe26Pih4+YRcbuo8j9txn2m0e4xqt6r/cbTavk+92X/u2pVtDuv968bY14xxhxLPu9T2wZ2Xu+2BsLoxc3rq8IQ+WT56qoco0zAQsrkDTellW8ZBVrGn1h0LA2ZjlTXIWiQTb5Z6/LalrV+s32PmyFhNtP3PuRod15vgE9baw8nn6/A8M7rvemoShTdvJnzSExmswu1R3yZWsE1A+U7hW9KhwgqK/qsv3XX0DLTSVRtY6iebiusMuQ8aBh0lflU8qmIduf1zkKc17tTVH0ouvUQ1WhVlXkBI73d//imu69s2zm/UMDKb49P4qI0/aklekFmm4Fe+0u7Wdegk2Ub6CSY8zFjzHcS0/zGZN3WmNe7yDTuNqqajHl+tnaOBfn9tP1pX5tqfagtIYe5JtK8AFWorUWmvLSl5n2LqvTN8zx/aDsouv7dCoAMC1nCliPLdonyd4A7gMPAGeA3k/UDO6/3wA7Um/cGDkVD9T5Vbu4if6eQWCNQRs9W55NfJwjVFVrnE7W/HFonuaHav1rkJuhUzef5OOW7WwSyWRH5LQj76er7tEWU1tpla+2GtfYK8Luk5vUiwz6vd17UFzbvzeurMr9dRQ+dVmNFCspXe/62rH31MdpMu7gGIX9klkLN2keWQz7VrH38OqtiM03wPP9tu8hyq7SDLUiwbRGlMeZm9fMXAImID/+83iXn0Aj657LKtYOQWgmRmHzyHvo8xarrylJsvsosGwUPmdih39q3mRdhzjPFtWr0ybEsmXRKPHnXwj9G1Xrz7rduKvuIIAr/smRe7/uBKWPMIvDPgfuNMYdxZvUC8Cvg5vU2xsi83k2GdV7vPPVYVVl2wyelydAntjzlV1YdSv6kv12fq79tlZSUyqrdvGvhq0h/yDZNzD6ha7+ongpAbw+dnz5ulq+1U3O5l+rKvxcHhfAGpR1dxLaY17st/2TVUYOKfIZVfYr+vnm/y+zjtyNPkYR+631kugdN4Jps/aknoNy5y1QYUt5XrSGC9OuVyL1uW5aPlcDvEPJeIFXqKVN3p+V62YaqdQ2gCZ7nn8yb13sAT2WAUIbcyipM/61f9spnqaAqD1gn5UNmNLSayg2vbNmRzv3gi14X8rH6vrlQPXhlitRu3rUJ7Z91vHZQhaC6GQjqJYalnRURiTILq1QbpLcbZpquR9BOnUUmtyayIkLMM9fXVV3iZ/TzFiGf1HxCDKlGf1l+5/lmtaotIsc8X2uZ69/uf17FZO626uumuV6mvk4sqgHAEDe9HDpKC6r65/o+vbwAS+gYnTysoWPodqxybUK5XvbV4DrXdmvUZdZpJS/tH8zyJWYtFynG0H6hevxrFVLxeX7XdtFJXf1SX908bpG6z1u/iWgnLUgwAM0fcFQhyzLmYF7QpezvrLpCClIrLP9cskgVWn19PsZJ+19P4IhzBZjz6pP9x7199fbQNcsjztBLxfdn5r2o2iHLXpmT/TK9+2keD6myjMOs5SE0vayPUFAh7ybMUktVUYb0fLXaxClBMZn1p5nxwfuWsqtJPVJuBlhS9Wn/3kmvLbWSHwLl9Tq89ZJkXiM89mQROYf+l1AAKFSmHVS9D/rhG+02htSHOYTcvslYJY2Ay3JRQKbKVe3EV1REyHDtQBeaFDUB6fQZHcAYD6wTJQmOHGdwXQ1mVPn1ZJuQaJEprcktj1R98tYk6dcpfkr/4cxzZ5T1X+q6OkUv/aDtHq/XGDJlOURNrY6udVuUnEExMzW5hEzadlBWkeryZQMUq7i2CuELecoYjxNcqyDlnFZI+0z7bRQFt5LULSb5SeBuXDDMD6ro/uFSxzqtgyXXaCVOn/xCJnSD8FBx/vFq3nLWf+dfS70+j2iq3AvdIquiF3evjtsu+sA6nfgnIZre5aHzAvMe5G4pjE4hZKfnw9Gms5jhopQXk99CqvIt5yvmusYKjgxFSS6SkohOixWTXAhRq0AhOE2iE7Sa0NqkhlaFKesnvPV55nzI9NbXLbRely+j5EMIuTOKUMb012W7VVdECyJRtousG65T0gw9qFUDQEKIkKpGUV3rOIJbwpGckFgNR3RSzwqpP3OFlFgFMzjVKOuElOdxQ6WsJ2WO4hTmYnK8xaS+cVLS1HP2SNu1f3hCfY/nfIQ4J7xvTary0ggRL+STaDcQ8r3mqdAi32noPshaPygYpLaUxJY2vXuKMn4e38zzkRWQ8cky5GMLLevf/hBoIVfBEo6slpIydRxBzSXL66qOGVKVKfVLnceAT+EGRD2WlHsQ+BPgl0lJFxxRzibrHsUNzb+EI9ippE4x48UfvKKOO6POSQg0dG19k79Ga9v9a5jltigys6sgz8QP/ad5JFrk2+6mlaPbkeW+yGuDvPyq9NIaMAxhk4cMRTdWnqlXdCOitvt+OB19FhNbbtQpXA99fRzxW0pgRshlCkdiMziCE2Uo7Z1Ltn0sqedw8j0L/H3VrruT9jxIapY/qPa5R7VlFUeOuu3SRiFRPYK5r8w0QaLW1bzlda+Mvo5lyLKIrIpQhuTyyuUdt9uqzVfcTa69dnrZV846KNjtUfEL0Kl/ErYwUfZ8/MmQAgk9UP7DlLUdWh9mv0yRqe0fS4hyilQZSj11HME1cCQoJqpOJAdHiuNJHSuk5rBPUBM40jtI6qsU/6X4HFeSfU4mv4WY55J14JTlvcD7gcdJifrpZNv9qi3arPahVbN+YUDrdZDgnL6GIfjbJACVpfi7gc0kwRDyXgKhe1m/aPIU8CaTZLewZYmy5+jWldM3l05nkfpDUyv4N7H0otGpQEIkkqLTVL+16prDEZsOhhwnNcHnccR6ONk2r9oo39LVU/sWV4AjOOWqyfJgsr2myospfi+pn3Mm2Xcdp0al7X+AI+UajjizRl+XoFWWj0+X1QgRRMgK8E3bqmRZ1lLwUeYYnSrdLEWbd893oqp7yELdUJOwRYlyU0Yzz7qJ9Q1eZFbLNr1dj8oD4bQXTa6ibHQe5KpXr5isfk7lBGluaJ1W5SXEJj5DSfk5qH7fgyOzJRzZ6ai1KFnxbYoy1V0ptSKU9CBJY5pQx2rg/J8P4si3nhxbiFai7jULTRM2/TRp+rmVoZeRjzxi0/7QMuhEEZYl5BABlT1uL9Rs3n7+SFMDiC1JlJuCLGd5yLTOMlU0sojXDzzIvpokNVnq7Us4wplTdfllm6Rm9QROzUk+5FSybgF4gFS5LtJqZt9DGhSSB1QHXbTilNxK3U59PpKXKcpTzv/XSN0GEkU/Rur7vAdYNTB1GZo1qJlW87ipPqHcTn09sv5bDf9/LbOP3lejinLMU6955m+n5noVxVxGJWYJhQFFmXm9Dxhjvm6Med4Y85wx5leT9XuMMV8zxryYfN+o9hnYub27Bv3gCUK+Rb3s+yP94IH+nqD1QRRf2yqtZqUoNQlyaFUoBKVThISk6mqdRJeljeIDlG6K0pf7KVIyfRRHWHUcaeHVJ3VJSo4Q9kyyrFN9pLxsP5y0fcY64puyqVKdcW3d/aGllOCfBP41zpd5dBRqzbQt0gb9rbf5191H2Ye4LImUUa9Zx83yc1dtQzvohonst99/fgYYZfIom8A/sdb+OM6L9A+T+bsfAZ6w1t4JPJH83vpze3diloRIVW6WUB6hNhUhNaH9IITUKyavvP11fqEEakLtEQUoAZc6sH4ZmstpPQtJmf/39TToAinRrQIfpjUaLXX7+Y5yDlM4k53ke67B+D3nYarBxOyr7Jk7zXh9jdH6GtRt6u98EC6s73R33BzubhPTH+Dp0cTvasPEmKVi/PV5L77QvlBMKGXNy3YIRJ+rvi9CFkk76ITU5GUr9/cC14qFAUbhJbPWnsHNtIi1ds0Y8zxuCtqHcO50gM/i3uufQM3tDbxsjJnHTT72F91ufAh9n21RKzpoNZ/F36jNI21yNr39tZKE1CSuq+21wH5aJUqbxGwWMhQ/oUTEhbyaJOlBo1CbhlULGLftXoAbUnM8UXdXe+ccU+1ZVfVpn6S0fQIYt+yZO83qSp2dExcAqO9aZcfuM1xiB01GmNy9xiRrnJo4wI8WZmACrpt6g43mCNQbcO8IrIym0fNVnB9zBWiY1peNTk3R184nOd+8zTM7/Tr0eqnTR1V1lkckofbq9aFt/SAm+e/F762T/3tkencrkAMVm2iMmQN+EngGmE5IFGvtGWPMvqTYLTgjSLBIYG5vY8xHgY8C3HrrrVXb3R9kBVb0VfRTVrTZ66sbuWGFPBZIzWW9XUxXvxdMU633k6kbXtmGt79s177FeZxvUJRrE1g0rk2SSlQjjWxLz56DpCk+0va5ZJ8ZUsV61QS+zHh9jand59jBJfZPn+YCO9mfTNi5zD4OsMwak2wwwkV2cLGxg933LXGxsYPGyo1cN3GB8YkLNFZudOY5o+mxFkmDSDMWxk3r+JplCDHL15xXviwB5ZFo1rHL1FeFrHsN/wUE6YtKXDsw0AEcjdJEaYyZAP4D8HFr7evGhKbwdkUD666ZmMda+xngM+DmzCnbjr4i1LvAJz3/ivpR5pDjX9Rg3VvfJH3AZT8xu32FVFf7iEqVeiSRXLokisKS+mZJE8khJZmlZWDarZd8xikcMc4lZfU+kjs5pbYp9TA68zojtSaTu9eZZpl9LDPGJS6wkzqrTLPMYY6xwQgnOMQGI5xmPwCXdo+xwQiTY3Bx9zobGyNsNEeozaxwsTHG5doG19U2uPL0LrivAU+PuwBP00DNJudvWl8c2kWhswZCChNvGW85T3WGyufBv4c6Jbpu19cutJUj96pkRAw4YZYiSmPMKI4k/9Ba+x+T1cvGmJsTNXkzcDZZv0if5vbeNLM7RHZZ27MeoCbX1iOm8QypepV1496+8nueNDqtj63JVQhygtSf2PQ+ElzR3RaZdl8rOMV4GGdeH8FFm4+pttZxieJLan9p/4xlfOo1pnafYz+nuZ43Ocg8TUa4xBg7uMg0Z9nHMiNssJcV3sU3+To/wwFOcYJDzLHAWfZxkTEALo647x1jlxipbXCx1qSxvpPx+8/TeHSPa8eRy9AYTU1wHfXXnzImNh1sK7M9D2VIOOuYefdpL5GnZkOjPFV1R2wyCptmnHT8t8Dz1trfUpu+DDyM6+X7MPAltf6PjDG/BexnGOf2LoJvhoVuRv3tBwhCD6YsC/Gsk5q8cO00DvJ7jvTNLGSwRKrqNAnqXipCGvq3Nt+XAud9HOeVXsW5CerJ94I6h7vV8YDrZt9g7/Q5plnmECe4ngsc4gRrTHKJMc6xl5/mG4mZXeMApzjECS6xgw/yJ7zEHRziBM/wTvZzmlMc4BI7WGOS10bqbFBjpLZBrbYBi+M0ZmHiA6+y/vmbkuvWgMaYI8txy9W0oSLfYVVSKUOiWf7EzUI/TfGi44bcWh2gm/5JKMfh7wb+B+C7xphjybr/BUeQXzDGfAT4IfCLsIXm9g5Bp3Zk+YV88y0UMNDrQ/tCqvwk/Uf7BnW0UPyJkuDdII0ki1nv9/aRfccD9YkClbbcTatqXcKRuSSCH0yOJ3Uv4UzecRifPc9du7/HTi5QZ5W9nOMg87xGnUnWuMQl6qxSY4MbWWWOBe5gnttPJiy9Dm+983lO776Jwxzji/wCcyzwEnewnKjdi4xxaWQH55f2wzhcV9ug2RxpyeMcra9xeeWGdIW8JEKo4mvMIr0832O7JJl1HxXtE2qDX2evUKXuLpJkL1BIlNbapwj7HQHek7HPbwC/0UG7BhO+cgxdvbx1+sbJ6poob1YdCdfbmjgVuURrFFkPfKEfdvFprqpt2uT0P1KXBF5Ovg5MOl/fQnLsJ1+E+p2OIJ/CjQA0hSNVuDqoRmN1D9+u3cvb3vn0VT/kd3gr+znDJGvAGod4nlMcYI6Xue/FbzknzTO4V++Pw+hpeMuPvcrUbefYO7bCN3k3dV7jWY4wQpM32cmJi4dg3sAsXFnaRWN8l2vH/CgchMurk05ZNkegNpq+JPRUFvq/0oQWshD0/6X3KYNOTXC93C7h5v2mzXq3OLbMeJSb4p/USdqhBymkRvygTyj6revzAz5ixso2HZTRx9XtaXh1Qqtp7ROEfz5CyuuXcVlfBlYvOGJcWAbuhD8FPp/s/1RS/igpmTbd9dpz5BXmeJk5XmYnFzjIS1xiBwc4xRTnGKHJBiMu4v0j4CvAf8Yt/2ecH/QJ2HX0CrefWuKn+QYHeYnDHGOMS2wwwt6xc/CkU7Cjs68zMfvq1eT18fqaI/nGGKyPXnvO+j+S6wmtLzKfJPX6vP+81wjdU92qt8r6EEJk2+12ZqDbZjcMtPt0ACHd+rQqw1sO+S5D5rn/cIrK8c0rIcSG+q19jk21LCPraLM91Ea/y15Dfa+r34zi+i5eAHYmdU+nx5RA0wSOT+9J9jtK2gcbOMN+DnCKHVziv/DTvJ1n2aDG9VzgDPsZYYMNavBG0sazyWcPjij347adg7e88Sobd/91EvhpcoJDPM8hvv+xBo31nczccpqlx26Hgy6AtNGscV39DQCu1HYkwR1S1Rwyw4sUVRl/ZGhdL5+2bqtA/37W63t1rAHGEDRxgBAyvbUK87fpmypksun1WWbcAml6TxNncgs5ySATEsjxlaLOnRTS9CPWFXAAACAASURBVEletulgkE6V4c3kx3loXA+NJkxMw4pNktebMDWaugskCJV0q9w5coEJ1qgnDtabOc03eRdHOMoak+zkAtMsc+fZRdgAbgX2Ad/Fmd8AP44jztPA63D7d5c48HNLnN09zSXGWKXOW25Z4Nwbezm7PO3asGi4NDHG3PQC3//Bj3Hd+CWuG7/EFekHrglSzllHxfPIICso46/L8l36Qbx2sRkmci8Yol/BrA6wZUzvTUOW+exvlyBKFjHqsipK3GI+r+DMXRk89zhpis8K6TQOsu+K9y0O8lBbZdnPX2vpWXQBJ+sE1wPLCcEY4HlgtDVg1OTqiOYzd32fETa4xBgvcBer1KmxwRwvc469TLPMHcw7Er0IjABjwG5gGmwTnj0Lrz+NU5Yv4sz8H8Lon8HbOcpOLjDJepKUPsaVxg4mHnwVJmBueoGzb+xLCHKEK+s7nRkeCozk+SM1sqwEQdYLMVTGL59lyod8yf0imm4dd8gk2pA1N4y+dVvMUpj+Nq3QNCHKNr/sEimRreOIZ4V0jpoFXLRZDy4hI+mIObyutgtx+w+5319cTFEhWq5XDT2PI8Z3c9UU5xCwnJrjT5OOJDQLS6/sZ6K+xsiuDTYYYQcXmeLc1XxJgLcde/GqWc3LuGSyH7kjml3usH98EWZfhL8FqVk+AW85+ir1I6tMssZtLPDmxE727X2ZF595G7PvfJHvf+knmH3oRa7f9SZrP5qg0RxJrsVoK1nq6H8ZAtTf/nZBlrWRV7e/rlPl1W2Tv922ZF2bHqAX/kmIirI6NOlp9ddirnrbQ+ZulqqBdCQgXb/0fBHTVvsiZVgyMc3FP7nKtcde937759ASjT8PvJ4sX49Tl99Klq3bVktIclHtP4NL9j7uGOjsG/uumt6n2c8Yl1hmmovsoHEnsDc5xI9wWbgN4Ifw4iknpP/OiDvi11506zmNG4ZlHO7n60yzDMDFxg7OvLGf2Xe+yBwLvO2hp1n88zu5dHEHAKMTb7rj6HQoff55xFZENiGfs17fjhrslEz6LYOK3BdDhEiUVSG+Pv8B8JWiIC/Sp8vKkGYrOMUoKlJ6uohyFP/fOOnI5JLf6JOrHoRD2qJNbT260Lgq14S01+mbpGQ5CtyQNNDAzA1OQb4fR9T34BLSDwILozAO60tTrK/UOXHuEGfYzyV2cJEdnGMvC9zGsV1v47m7b3e9ffbj1GLip7wB+O/vhOYGvHUEHtgD7IEXn8GR6hNc7eK4g0tM7T7Hu3Z9k+Vz05zYOMQqdajB3NgCG82aG0gD3KhCfgJ/FtHpdVnbfOSZ592K/GaZ4z4Z16hGzt1GllruttrtMYakmdnYVLNbR6X1Hx26ir5aDKFGmkQudU/Rehzdo0YGrZVtizhzewbnyxRzW/yekAZ0mrSSoY76ar/r1bbK3fwKMIkjzGmcttub5k0ewaXi3JO0oUnqGlgBJgzMbHB56QbmOcj+vad5kp/hAR5ngTl2coE7eIlX95/lplvXXbeFk8A5mN6bfO/Didm9wG64816cL/M9uJ45bDDJGm/lu2wwwrv2fpNz7OX4E/8Vb3nPSb79g7ezZ/Ysqyt1Z3bLi8Tv861P2//fxH2RhdALUr57QQplzNes7WXcAWXr6hS+0BhQDD1RbiqyhulCLftqwidMXU7ISk/R4EejpT/1MdJAjhChRJqfxA2BdjIpO48jrQmvXm1+C3xCvhoNfh24jCPJJo4kE0UpIwHdTeuQb1Okw69JN8p1GB2/RP0nzrL2owkusJNplnmBu3gvj7GDS+zkAgvMMXXkOKYBvBXnsxzHmdmncZHw/cBtuG4OY9BwYpK7eIFz7OUEh1hh79X8Shpw7o29XDd+ifNP3cLoPa9zRTIGlmgdSESgXxiaMGWdvgd8hMgrjySLXqR52Kx9OtlvCyESZbvwo9/+fDQhc9xXLrJvaC6ZJq05ilJO6pAhzhaS5UdxKY+QDnt2kHQebT3cm6Qa+cSu27aayDlG0wK1vY5cDpIGk2SkIYCvko5EdJCrBHx55QZenb8BxmFjtsbE3jUusoO/5jD7OcMaE+znDN/Zc5EDnOLG8QbmH+DMcGm35FLeCa9PjzLS3GB5bB8H+CGTrHGCQ1xkjJ28yTeX38Wh6RP83M9/ie/wVtaXplw7Tt7g2rSQtO8o6eRqci2yVKOfwZBlUvrIU36dEJD/5Fatq4z7oB3l2Uf0KpADQ3UZBhD6ZvO7B4ZMrpBJJwpFKz8xlxeT30dwD/c9pAnhx5PyEu2ewZHjSRyZ3Zdsly6MQo4S5BFyhtaJxcTvuQrO1r2cbEh6sc6Sjkz+960br3IhKS/9wnUXSLk2SXfLnRMXWOA27uIFzjLNTt5kJxc4xa2M0OQcU9T3rDK9Z5n958+zsmeC+o8cW57bvTvxb06xd2SF/W8sMdKE07v3U2eVKVa4wE6uLOziwPQpTrOfpWdvZ/zgeRrze9y1eYDURSDf8j/IoMh+wA5v2bcMyqAT0vHJNuu+6iZC599L5Cn1AcBQE2Xf0oL8B0UTJLSmmWilqXvECLQvEZxKlADOHOn0B4dxD/o4bs7rw1w7DFsTZ4IvkQ6Lpk3jUDdGXwXPqLavALXRVNXNJvsnxLxn7jTna/uAUdfORVJlLMP+10gH8Z1wUzjsGLvEMvuY5iwnOMQITS4xxgK3UWeVVeqcZR/f2FPnLr7Hud0XeIG72M9pLibdH99y+lUau2H0FNy4ezUZ4HeM//SD/5bb3/kcF9jJOfZy+9uf49S5AzAOs7/6IovP3ZmOvakJUVwaoT74Ap2kHnK95EFeWFnmub+/H1wqu22z4LsrOkEZf2ufMdRE2RfoN21WcAe1XkdXdd6iHj9S19WkVZ0dTMo9TWswR+daiom+gCNHKTuOM8lruGi0qFHxbc6Sqko/Od03y4XwZBShe2ByZI0L9etpNPe0XgdxEUhgSs6tASO1JnVe4yzT1NjgAjt5iYMc4gTf5F0c4BR7WeHNZBCNSdZYY5J1JrnEDu7gJXbyJq9Pj3LD05dhH6wwxRqTHOMwN73lDN9/7Ceov/c1dnApvfYTsPjnd7aOgKRfaDo1KkR++mVY9aH2r6NGltIMvYyz6u43yeQdfxDa1wVEoqwKMRFEfejvJtdO/aAfRFk/ocpKDxxJ9JYAzHEckS2R5k4ukk4PK8T4AZzCXEzqn0+O8SCOBN+P88UJwc2TkmiN1tGHpP+29BuHNBcT0vm7kzb/4A/vdoS+QNq3u6n2l2HYJHUJuNQYY6F5G9N7l7nATtaY5DT7WaXODpzSvI0F6qxykTFOcYANRjjMMUbYYIQmL3EHPzbyAuy+zOu3jzLFCmtM0mSEu3iByfeusZM3OcWt7OQCY+MXeevbn3Km+GO3u3MW18WMui5+ipRvfvo+yqLAjg9dPxQTSBnzt1t+zl6RWdl65Tlq0/zupX8SIlFWQ5bq0v5FMW/1HC2SHC6QPtGiJLWyO4kjQnDkt0Q6Go9+mCXZ/Mmk/AqOHOWGOwY8tQwPTqd1zZAGd+QGFoKWhHZRg5LDKcpTCH6edBoIaYOoMTlPIVlISTIZN/PK0i6urMJi4wZOH36DHeMXWR2vs3/sNDu4xCFOsEqdJiPU2GCEDeqsssw+Jlnju/wNDnGCDWpcPgDzIwf5Lm9lgxHey2OsUmcnF1hjkjEucu6NvYzUNvjWC/clE47hyFuyCJ7EqW15cWnzWpCn5uS7jA9Sk6/eP3ScskqsKLhU9gmvqvx64R8dYDbqZF7vXzfGvGKMOZZ83qf22Zrzemt1pc0zUReQTsOpVaMsSzk9Fe26+p7DKUBtEorpug6wnJKdRJYfIPUbCklLsOWXp10gSPqEN3C/F0kH+NU+ulm1/0zSnrnk2DM4/+d9pMEeOad7k7bUSX2cWrGKStZBo3saXGmOXE0ElykeltnH9bgeNBfYySp1zrGXlzjIKQ5wln0sMMcGIzyz+6d4iTtYYI4THOIUB3iTnVxkjIu43jj1Xatuatt5HFEukhLin9KqJvH+j5A7Rb9g8h7sLOLKgx8I1ERchFCZKoGmKuXbOUaZNmiLbMBQ5jSbuHm9v2WMmQSeNcZ8Ldn2aWvtv9KFvXm99wOPG2N+bEuMch6KZvvLs8lvHSjRV1lIRwd7JnBEJ+RzmPRB1YR7cDo1f+/Dqc+7Sbs2ruDM7F8jDax8lVZyq6tjLqn6F0iDOWKCC6ncS5rwPkvrKDu6b7iQ8RzpqEaSAzqDG9AiGVXoutoGV1Z2cXl8DOqwOl6HEbjA9VeHZTvAKZqMcAcvscYky0wzxkV+mv/CAnOcZZoV9vINfppJ1vged7GXc5xjLxcZ4/Tyfibra1xev961e4V0WtsJ4IPJOcxzbRBHrpNGXsAlazmrfBZpNihXV1b9VRRelhugSj3dNtkHVFUWNitnXu8sbMq83n2fv1uQ9Tb21YB8iwkuD4T4BP1IeUOVFXNaSFbM4Q/jHnpJ9anhXlH34QbVXcSZ4w2cH/MIqQ9T3AKLuIDP/cl6UZ+LSb0P4gjuIOmQb5JutK6OKwEomSxEpxuJK0GIdR2uPL4rIXfD5YM3cH7uBs7P7GNiapU7dr3UEog5wSEusYN9LPMzfJ1jHOYF7nKmNVMc4BQL3Mb1XGCZfW7otYt1doxf5EfzM6lbQPrAz5POPqktBPkfJlRbs/Iq/f89y/zOMlF9943A99FVMZ9D9bZDZJsZfPGP1YaPstf+SajI39683u8GPmaM+R9xOuafWGtfYyvP6w1hh3PWTa8VaNNb1o79CVofODH91tX2+0h7k8wDf5CUezIpI77CJdy/cRBHjFodHk/2EXP8KVLz+BhObR0nDcrMk+ZzTpAO93YwKXdQnav4K6WMBKIk7UbORVSmKFxR4IvAVI2JXWtXe+qcY+/VqPdeznGWaY5yBICzuGnkT3GANSavztC4xiRvvnE966uTyajmSXsWSZPL7yadh1yGrpNeT5A9+ZiG/z9nEVoe0YXcOHnHyUNVQu0U3SBTeQ4kTWsiv3g/UfrSBub1/h3gX+JGT/iXwG8C/xNbeV5vaCW9rO01wkpCR76zHPv+NolQJ8GQq/5FcA+7+Ns+SJoK9ECyXco2VR1zOKL4PRxh6mj646TdEr9KSqhztAafhFyESKWHjpC/5FFKQGhV1TGOI/IJnIthhdS/umBYWr2d2ttd3+1THLiqLE+znx1c5BQHuIOXOM3NnOIAY1xijUnWNiZdX25w406um1SFSwR+njRopYNj8n9JyhW09gPX/4+G73rJ2q4JMcusDt0PPpGGtoXq8n8XtVOjakCnE8Ic974HGKWIMjSvt7V2WW3/XdxjCpswr/dAmd15N37o5vXVRyhgAK35h2IGN3HkOIczicWkFVP5bnUMIYk68DmcIr03qeMozox+EPjXpGb0QdW+JeDRZNzJWVpnghSTWrpRyraZpD0nSYlzPWlXk1T13kOqcEWBjif7TMHiKwegOQLNGuNTr9FY38memXOcn7+bl+pv8AzvZOfEBQDWVyedv3N9J6wmvYREocsxx9VHUrAkU0BMconY65He/YBKiGzynqCQWZlFfnn75hFqWRSp2nZIrxOS1JbVEKCQKLPm9TbG3Jz4LwF+gTQzbevP691uvpf/0PmOe//hEN+amK4SjKmTzrt9N2kajx+NlqHbPkCaAnQsqfcR3KttKtnWvAwTo62RevmWaPUKaTdKCeBAGriRc5K2ioKVwI+Y4xJZ918YS0m9E+NXCb8xtQcacH7lFmjClfVdbvi2xWSmxUW4ImpQ53JK76YZUoKWYI7OSxVlOUvqu5XrXfZ/1m6VssSZ5bssQrv7Fe2zmX7JLmIz/JNQTlFmzev9d40xh3Fm9QLwK7DF5/XWKHooivxV2k8ZKit+QVF5J3GEdQ8pAQiZLpBGymfU/klXw6se46dxhHAQR5hTSb014L7RtHvkeFLf0zth/WWo3QbHX4eDNziy0QpQVOYqaRRdiEknuE+ofYSEdbqQTrcRE11IW4jLD7z46lnIUncIOE6qlmdJze8lVd89SXlRm1K/nzeb9X8Xbff/4xAphUzlUEBmM4m1VyirqgcIhUSZM6/3V3L22Zrzemv4pOdfySyfUJZvUj/8WpkJoc5xLZGI//GI2k9ShiSIIgR2P+nAv3W13+M4ApnBKcw5YH4Rnpp15vqx25J0mhtaz2VWtfs4qcl6MDmeKN05Wgnv6oAbCaQt/rWSa6oTwEPKrYl7UYjfVOqpk46sJIp7jjTSLYnwuq+8PgaESbwKyqhMQcgd4y93At2WTk3mTvYbROIugXb+/r5iYPyTkO9rzCsfIk9t7s7RGoTQydBahYp5qB/mg6TmekOtW8SRxOdxvkmp7x8k38eS7/mX3QHuI01qXyAd7OIorVNVCIlNkY5MNEvq+5P1ElFeUvtq9ab9d1pRCQHrnMYmqek/o35LYEgU9EHSqD3qW66zdmeQtEECOpocG97vItLRLzvIvkf8ffNItRvEEvK3FqHbRK3r6zDKvVlmN8SpIDpHltNdb9c3iL8sDyukKlIUjigweeC0YhLCmVJ1nCRNHJ9S+66oOmZwZCejDzWT35LIThPuvzPtYQSw+nLaI0jqqKnPLNcOECxmbE1tE7O4SRo48c/Lv0aS2yjXSpNbjTS3U8z+x5Pr0EjOa5GUCBdJXy7iI9U+SLluojSFpDXySKbmbdcvRV8l+v5pvc0nlH6oL78N7dZRtDwkGDpFOdAoUpWhB0hUolaEen9RmvoBEtNSHjZRW3M44tBvbyFKIeAHSBVYDZePWccFfI4DM3e2JsOvAlO3OTKU3EMhmXXSnE0hHfGtos5L8jvF7BUi9QM6QsxaTQoJTpEOtCu5pGI2J9Fy5kkHA1kijczXSSPxktsp0OpQpwQ1aU34F/hKMs/lklembDCwiGDaeYKLiKoKkbVb1wDnTIYwVIpyoMxuSAmnigKQVBRRUwdpzZUUE1Kg1ZY8zEIiclxJ7xGyFQLSvXkOk/oDpdxTpIP8Pk6a5qMVaz3ZV8i5TivRa3UpSlhMVSGDBXVuct4S4NHXS1+3RVL3gbwUpLz4O+VcxNwXk38KR6gfUOfwAK3jd/rHFH+uVsAhgqvyEvTdLDXvW5cJuWRCZf3lkPnuWy/+p1topy5xuQwZhoooBxK6C6IgdMOL6hsnTUWR3/KQirmcdSNpv5cmGq0shbggjfBOkE4NsYrrzbMCHF90Zf8URyxLpOk7zWSfRVoH4xXTVNRjnTQFZ0pdEzkX+S0EO0GrKS3nJOtQ+yyo82smbZtP1sm+UqcQ5UHSiLY29UnOUa6PmNlCjNqf6hOaDqxlqcQiZRciyaxyeXXnEXGn5m1ZIu2UcNeLixRhM/2TEImy+8h6ECSAIDeYTo+BdBAGbQpq+DemJlrZNsG1xLBIqnwlt1FSie6eTfMshaylTVLHLIqQL6dJ7eDSanT7dc6n9tdNkQZdat63jmhr8tTuAlGA96h9RfWJH1R6B4nPtIkjVSmvX2SyjyZC+X/0uYv61Qnr3XZWFdVXlZSqtq9Iafrkm9eWLAUbOsbqNXsPNKKPspvI81FK17mDye85UlNPD8SQpUy1z0vKC7HMqN8CMd/FFL6PVKHpXMfDpIGgpUWozzqiOYIj2flkGeDgaHqeDdLpJvyAiCZKUbsSUBJ1qglJR+6ljhlVVs5BlsW9oF8s4sutkY4Cf1jVL2rXv5667fr4EtDR1122deOp0X5k/xid1JmHKseqok6rtr1Jq594CDA0inLg/JMh6GCEfIv/UXq1CHmIWQitKivPXwWtKkibsuPefuKbq5OOMzlHOtLQLE4dggviLAH8nlsvvW+WLrhlP6ItpCbHl2MKsWk/piYbWd/wPpIorpVHw6tHFJ3fHm1+y76Sdyq/xWcbgvaT+n7Gmrde2tUuqWW9SEN+Rv93FWWZd/+E6shSfmUUYhkl2mVsttkNUVF2HzofElof3qVknQRVGtfsnW12a+d9SOn45cfVN+q3kI74+xaS+ngd+PU0xWYOOLgzJTjJ79T+ulmvTp33KcSljwutwRnta9MEJOt11NknRz9HUZ9jnWvTf7ICJSGC9IlCX/uy0eoQsgIwWe0K7ePvm/cE+/dI6PjdJLoekmO/EYmyFxCylGX94AvR+AjdvH4AQH/rB8dfrwmqodYJuejue3cn6yduSP2RD5Km3fgKUroritleU+v9NCZNklKf3wNG150FP2Ah+2Q98Dq6X4Qa6f+i26tN8ywTuZdPT+hYWWowry1F5Byqt12SDBF8qNwQEurQmN5DB0mOlpui4f0u4/cJffJUgr5JQyajtEFM1Vm1v+QbalNeotc+UdZoJVHtYxSyFFNauwK0mR56GVRRN77y9pWlPu8s+ETrXz99nInAfr2WGaEXZFUCCq3z6+mWzzVUdwi+v7gC+mF2w5AoyqHwT+bBN9/wlkMkUcbHlPdghAhZ6hTXgPYhav+hTjGS/eT3LKk69NOYRIEtkAZSIDXLNbHq9uh9dZt9ovBNYtkWug5+gCgE3cNHXxNN7kK2muT9tvVSIWk/bNE5k7G+zLo8c79dhOoZsiCOYCiIcqih+2dn3SRlFFOZ8vo7z0z0e8Do3Eg9QIQmGykrJrZOhJfy4lIQc17I0Yfu3qjblnXOWaooaxtcS+IhiE9V6tOBHa2U/S6WIVXsv8Cq+g5DrhR/e6dkVsZMz7pn8v6DIgxhgrmPSJSbAVEkfv/ekNIsujl9ZKks39T2ywsxSj6nENeMKitpObqdU6qsr7D83Ea9Tfv9ZgPn5iu/LPLUhJFFqtrPWAQZfESOKf9RyOyFa32ZcsxQ2/32ZxFVGb9niCxDxwot+8jblqVWi65lXru3AAb+NIbe7BaIsvSR5UgvUlm6rO4Z5Pslfd+Vb/bqAIuORosprAlOeuZoIslTKSFCCJX3yUiXD5FiGcVWFvJikLb5eZV5Sr2MKZylDjWKzifLH9sO4eW1r0hxFrWzneNXQL/8k1DitI0x48A3gLGk/J9Ya/+5MWYP8MekQzH8UjK5GMaYTwIfATaAf2StfawnrR825CmHIgUQUiaSj9nkWiIW8tSpSKKYtCLUy/oY2k+HKqvJI0t1ZSHkm8z7Haq36BjtBFrylF4ZhV+WEPLUWkgttnOMPFR5wVRVlWVelEOMMqdyEfib1tr1ZO6cp4wx/wn4O8AT1tpPGWMewU0u8IktPa93p1ggjTQXPYRCbroXDlyb7iPRdU0qsl0rRDH9ffhmsazzy4gv09/uHze0fxax6m2d5CdqVH048whBrrWfw1m2zrLkkWVSFx0rpOqLVGZZbCGS6wYKL4e11pK66UeTj8XN331/sv6zuKEWPsEmzes9lJgn7cbo5w5q1ScK0SdHSPsj633GuZZ8JeDim8myLD2GIExg/u+sBzmPNPP8aVmmZJXhx8r4+vJ8mX475HfIx1mWjNpB6FqWQZ469ZezjlsVvTj/IUCpS2WMGQGexT3i/7e19hljzLRMLmatPWOM2ZcULzWvd6njbhX/pOB+XK8XPYCD9v81vXW+31EP2iDkOKN++71hfIiyhNYeO3KMoruhiulVVLaMjzIE3ac7azuqTBmiDrUvTwX3kiyL1lU9bpE7oZ168nyYsn0L+Seh5CVLzObDxpg68EVjzD05xUvN622M+SjwUYBbb721TDO2BmQEHxmMVpvYTVpTKfQNJze47po4TmvOo0+yvqqTffx0l6zlohtflG+ZtKciMzRru65DUpOyjiXXJ3SMLNWW5f4o+9Ko6suTfcqSVTukk7dPleMKfB95mX3y3EpDiEo9c6y1Mprhg8CyMeZmcFPXAmeTYouUmNfbWvsZa+0Ra+2Rm266qY2mDykWz7kr1CAdUWfd+2jilA+kilAi1ZL0HRoCTMr65Oqb4z7J+nVkEYvU76c8FcH3k2YRpN8uuJaQs1wB4956qP6gl0UVNaz38f/bMscog7J1lqlHvkP3QRV0su+AoJAojTE3JUoSY8z1uLGiT+Lm7344KfYw8KVk+cvAh4wxY8aY29iK83p3go/vhfnL6YC4ojC10pSxGPXAGvpbTO4sE1ug03tkJCG9rspDnlU+1B2txrXHCLVJly+CvFD8kYwEWf7SUNvLnHcVctDR/HYVW1n0guyz4N93ndYzxChzWW8GPpv4Ka8DvmCtfdQY8xfAF4wxHwF+CPwidG9e7y3nn2zB63ByryM7XyXVaDXD/Te7lAkFPXwzG8LqSpfN+5bj5qEo+OL7WYtIMuuOlFQnKeO7FcqYqPqc2vXVhZDX8yTPjeG3pyyKzrXMf+efexk3S2g/f11WQK0D9Ns/CSVuFWvtd4CfDKw/B7wnY5+tP693J/j4Xvjtc7C0NyUOPfJOVnK69JkuIknfL+eX85dD5fN8bz7RZD0oZYkr67dGnXDbs/YvUosa/rl0SqT+da/y4vHLt4ui4xQp8Kz9i9rmv2gh7eE1xIijB/ULH9/rvnUakE88csPV1SfLv5WlAvzAjl82RKRlTfFQHVk+QVkOjQ9Z9Xhl1lepL4/ou11/1kvO/930PmVfKp36Pbvh4yx6iQ0hBrL5W9vs9iBEBq3msh+k8U1o3+fmP1xiwvvdEMv6B0PrQ6a5vy0PRYPolkEeaXTTnPaP1w55ZF3DdtvRTfTiWhUdb4gRFWW/IdFo+ZYpHDRRarKsBT6hgW8bge2dPhi+Kgydi1Ykvoqt2oaqqTeb9eCXOY6+BmVUe9W2V/E/dqPOrLrLHqvNycQGwT8JkSj7Dz+9RqfvZBFj6OObjKGUoW4hi8DyHqK8dmQ9oINi7/jtq6LGfL9uVlpQt4ivarvK7KvvM39dWcwUFxlkRKLsJz7OtZODQVgt+v5LSMkw5FfLIkn/pu8GGeURYp7qzdpnkFHFh1vF35d13SDb3ZDXlixiq0q63fI3xmBOd7Gt/JMAv0wa6RZ1qU1vuPYmH1fb9IMonKbiCgAADS1JREFUQZs615JhluLLS/sIIS9olEUMWaZ6mQe3XaLpNqqYy1WvaV59Za5RkTuiirlcxnqpCiH6xWq7DYrZDQNIlNsSHyYlrbyE7JAy87cVDXThI0vJ+vDzIPMUUEixdvKA5aWqhJY3C0Uugypt6iRSLdvKKt2sdZ0SYh5kHvghRSTKQcGDtM4GWOZtnkUi3TKn80zqTh/svGOF9m+HCDshz7xzr1JHlWtQxvfZTlvKXMdeKvKVpP4HeniMHmOYvENbHzLUiEx6Ba2k6Sfu+jd9nh9Qj2tZo1it5W3TZfx2VCGGsvCPU0Utd5KOo/2MPtl0y2Wg0a0XnD6+T/jtuAXK3AshyL06pBOKaQyUotx2/sksyKjkviLJG3i3CCHSyLv5y5q7zcD6POQp1bzzaCfvr1NzXJNKlWsn+1RZ323Idc47h7Log1tjkPyTEBXl4KLIV1TFQe+jyo1fpFCqIkudFXVz68TX2Smqvgyy0AlZtVtPkUujl+jW+Q4ABkpRRngITQ4G3SfFdszITvx2ofqz5tMZVrTjjqhSvhcBl25ji5AkDBBRRrM7AzLUWlUiyVJgIUVYxoxsJ8gQioSXMemHHdptos+/zDlWyb3sdlS9G/BHi9oiGBiijChAtybf0ujGg5ZHgO2Y61vhAdMvIE18VZVit1RjVf9q0f5F5fQcT21g0PyTMPjiPQJa54kpG9ToxFwvizIqqWoQZgv5tboaxRb4Jn07Abmy7crKU81q01b53wIoM8L5uDHmL40x3zbGPGeM+RfJ+l83xrxijDmWfN6n9vmkMWbeGPOCMea9vTyBbYP14iKZ0Gku3byZh9lE7BRlXhCdKLkst0jIldKvbIBO2jBkKHN6WfN6A3zaWvuvdOF25vV+9lR7jd920PmVWWgnlWVQMCzthHI+226+mMr6k7uJshaD/t6iKFSU1iE0r3cWrs7rba19GTeb9Ts6bmlEucTdQb1hB7Vd7aAZ+IS29+qYIfRKVealF20TkoSSwRxjzIgx5hhupsWvWWufSTZ9zBjzHWPM7xtjbkzW3QJojbhIm/N6Ryjo+WJ6DT9w1K0HYZgeKD+ZfhiUVZXc2nZ6RoU6C3T5fhzEQA6UJEpr7Ya19jBu6tl3JPN6/w5wB3AYOAP8ZlK89LzexpijxpijvPlqW43fVmj3QWwnkNLLIbF6obZ6iayeSEVlBwWhjIN2yM3vJCCfIR8+rSzantfbWrucEOgV4HdJzetFKs7rzfXbaF7vdiFT2kLvgiOhCHWvHv5BJBXornk7SCgiyTLnJmk/WR0htjDantfbGHOzKvYLwPFkOc7r3WuUTQ/qVP34xNltotjsh63oenSjF1O3UfUaZbWv7LkVdZ2tmu61RVDmlLPm9f53xpjDOLN6AfgV6N683hE5kOh3XiJzOz1pihRlL0ltMx7ATtWyvt69VNp+/bWcbe3U3c5+m/BCG1T/JICxNi+AvUmNmD5i+XtH+92M4YOOgveCZHr1cPgP/maj26k7vUBW/mRZsiybLB5aX8Uv20X0myiNMc9aa4+EtsUujMOMVdJE9Co3dBmV2Gv1CJtHkn70utcPfzfqDpm/ZQN6RUnqZfbVvsiISJRDjyaOMEOBnqz0liK/42Y8HJt1jJAC64cLIe/4WXmY7bS5KomGjtkHguy3mizCQBDl2w8M/oUaCqyS3ugrpCQq8BVV6APhkW+6jW7V6aus0LlstirKuq55Jm1ZtdgptsOgJD3AQMWvssgyDsFWAX6fcJ3rJt/j6lsejHVaA0SCXhBNlcE9spBFMpupHn2UzSTI+t1PkowEmYuBIsoshAg0kmcJ+GNZ4v3216/SSqJSxie2bgVhOqlns/I9y6Lfx89CiIj9YeD6jGGwJoeCKEOI6rNH0CQ6QfaoRTVSX9YUKenJLJJNtSzlQwM75EVksx7iIsLeSkO1VUHonPMGSdmO16hNDC1RZiESaBeRN7SbkFQtUE7/9s19n9xCaSlawYYeZn3sEFl2iwDy6ilqYy9R1v/ay2uzzbDliDIL0XzvM/RUFtrklztQT7VbRvH4ZKnr2gxspiorm8KV56ONBNkRtg1RhhDVZx+QZQZqFarnMm+nrm5gUIglS1FnmdniDgm5Obb1094Z4qULIBJonyGKcyL57tVd2gvTNERMVYJWWcGXJu4FUuQO0cf1lwcQwxDIgYG/jIOFaL5vMtZxd6gQRFmlWYSsCH63VWS7qT86IKb37WQ6kIiOEImyQ0T12WM0Cac5TXjl9ARseXUVLVdFEdHm5XjqgJDvu90m4zwOCyJR9ghFJkUk0g4RUldFZNluhDpEdn7vn3GKyS0vTWdQfKKbiGExuyESZd8QlWiPoE30duGrTa32soIjUQFuaUSiHDBEP2gXoElLfH1FuZv+8oD1XonoL0oTZTJw71HgFWvt+40xe4A/BuZwA/f+krX2taTsJ4GPABvAP7LWPtbldm8rRPXZAXwl6I9FmRWRjuQYoVBFUf4q8DxwQ/L7EeAJa+2njDGPJL8/0c683hHtIRJoGwgN+hGx6Rgm/ySUvE2MMbPAzwO/AfzjZPVDwP3J8mdxk459AjWvN/CyMWYeN/HYX3St1RG5iOZ7RER3UfZ9+tvAPwUm1bppa+0ZAGvtGWPMvmT9LcDTqtwicV7vviNPfdpPRyKNiMhDIVEaY94PnLXWPmuMub9EnaXn9QY+CnDrrbeWqDaiFxAC9Yk0EmdERIoyivLdwN82xrwPF0O8wRjzOWDZGHNzoiZvBs4m5RcpOa838BmAI0eO9H+Gs4gWRPM9olcYNv8klCBKa+0ngU8CJIry16y1HzbG/F/Aw8Cnku8vJbt8GfgjY8xv4YI5cV7vLYIYPIrYrugk5vcp4AvGmI8APwR+EeK83tsRUX1GbHUMxLzeR44csUePxnm9twMigW5vDLLZnTevd8wii9hU5D0okUQjBhWRKCMGBjHyHjGoiEQZMbCIwaOIQUEkyoihQxzCbjgxyP7JIkSijNhyiEo0otuIRBmxbRAJNKJdRKKM2PaIeaC9xzCb3RCJMiIiiKg+IzQiUUZEVEAk0O2JSJQREV1ANN+3NiJRRkT0CLEXksOw+ychEmVERF8QFehwIRJlRMSAIJLn4CISZUTEAGPYg0dbweyGSJQREUOJ6P/cXESijIjYYogmfPcRiTIiYhsgkmdniEQZEbFNMez+z81EJMqIiIgWdItAt0ogBwZkzhxjzKvAG8BKv9vSAaaI7e83hv0cYvv7i7dYa28KbRgIogQwxhzNmthnGBDb338M+znE9g8urut3AyIiIiIGHZEoIyIiIgowSET5mX43oEPE9vcfw34Osf0DioHxUUZEREQMKgZJUUZEREQMJPpOlMaYB40xLxhj5o0xj/S7PSEYY37fGHPWGHNcrdtjjPmaMebF5PtGte2Tyfm8YIx5b39ancIYc8AY83VjzPPGmOeMMb+arB+mcxg3xvylMebbyTn8i2T90JwDgDFmxBjz18aYR5PfQ9N+Y8yCMea7xphjxpijybqhaX9HsNb27QOMAC8BtwM7gG8Dh/rZpox2/tfATwHH1br/E3gkWX4E+D+S5UPJeYwBtyXnN9Ln9t8M/FSyPAl8L2nnMJ2DASaS5VHgGeDeYTqHpF3/GPgj4NEhvI8WgClv3dC0v5NPvxXlO4B5a+33rbWXgM8DD/W5TdfAWvsN4Ly3+iHgs8nyZ4EPqPWft9ZetNa+DMzjzrNvsNaesdZ+K1leA54HbmG4zsFaa9eTn6PJxzJE52CMmQV+Hvg9tXpo2p+BYW9/KfSbKG8BTqnfi8m6YcC0tfYMOCIC9iXrB/qcjDFzwE/iFNlQnUNith4DzgJfs9YO2zn8NvBPgStq3TC13wJ/Zox51hjz0WTdMLW/bfS7r7cJrBv2MPzAnpMxZgL4D8DHrbWvGxNqqisaWNf3c7DWbgCHjTF14IvGmHtyig/UORhj3g+ctdY+a4y5v8wugXX9/g/eba09bYzZB3zNGHMyp+wgtr9t9FtRLgIH1O9Z4HSf2lIVy8aYmwGS77PJ+oE8J2PMKI4k/9Ba+x+T1UN1DgJr7SrwJPAgw3MO7wb+tjFmAedi+pvGmM8xPO3HWns6+T4LfBFnSg9N+ztBv4nyr4A7jTG3GWN2AB8CvtznNpXFl4GHk+WHgS+p9R8yxowZY24D7gT+sg/tuwrjpOO/BZ631v6W2jRM53BToiQxxlwPPACcZEjOwVr7SWvtrLV2Dnef/3/W2g8zJO03xuwyxkzKMvBzwHGGpP0do9/RJOB9uCjsS8A/63d7Mtr474EzwGXcm/IjwF7gCeDF5HuPKv/PkvN5AfhbA9D++3Bmz3eAY8nnfUN2Dn8D+OvkHI4D/2uyfmjOQbXrftKo91C0H5eZ8u3k85w8q8PS/k4/sWdORERERAH6bXpHREREDDwiUUZEREQUIBJlRERERAEiUUZEREQUIBJlRERERAEiUUZEREQUIBJlRERERAEiUUZEREQU4P8H9kwL1UWLWmgAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# read\n", "from scipy.io import netcdf_file as netcdf\n", "with netcdf(fname, 'r') as f:\n", " width, length = f.variables['dimension'][:]\n", " data = np.array(f.variables['z'][:].reshape(length, width), dtype=np.float32)\n", "\n", "# plot\n", "plt.figure()\n", "plt.imshow(data, cmap='jet')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.7" } }, "nbformat": 4, "nbformat_minor": 2 }