{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## 課題\n", "\n", "### 標準大気\n", "\n", "縦軸を気圧(hPa)にして[標準大気](https://ja.wikipedia.org/wiki/%E5%9B%BD%E9%9A%9B%E6%A8%99%E6%BA%96%E5%A4%A7%E6%B0%97)を描画したい。`import matplotlib.pyplot as plt`となっているものとする。\n", "\n", "1. まず対応する気圧をリストにして`plot()`を用いて描いてみよ。どのような問題があるか。\n", "2. `plt.plot()`の代わりに`pltsemilog()`を用いるとどう変化するか。依然残る不自然な点を指摘せよ。\n", "3. `plt.gca().invert_yaxis()`を追加するとどのようになるか。\n", "4. `plt.title(), plt.xlabel(), plt.ylabel()`を用いてタイトルと軸ラベルを付加せよ。" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAETCAYAAADKy1riAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4XGeZ///3rV4sy1Yduch2XOJYchISJYSQRpoL31BC\nNkvf0ALXfnd/LCzssvwWlrYBloVlKQsEAmHZULIJJIS4Jk6cTnBIHEnuiXs0sizZcpFltfv7xzmS\nx00e29KMZvR5XddcPm3O3Ee259bzPOfcj7k7IiIi8chIdgAiIpI6lDRERCRuShoiIhI3JQ0REYmb\nkoaIiMRNSUNEROKmpCGjgpm5mZUNwXm+YGbfG4qYkiHV45fkU9KQEcPM8s3ss8mOY6iZWbmZ/X/J\njkNkKChpyEhSCdye7CCGQQ1wc7KDEBkKShqScGaWYWZ3mtkzZvasmX0o3PUzIGJmT5nZx8JjP2Nm\nL4fH/dDMLNx+t5n9q5k9ZGarzGylmZWG+8zM7jCzejNbYmYfP+bzLw3P91z4vpnh9mvM7HEz+1sz\nW2FmmWZ2gZk9bWZPmtm9wJRBrusDZrY6PP4+M8szszHAd4HXhdf11vBznjGz/zGzF81skZldG8ay\n1sxujznnFDP7fXjOP/X/XMJ9bw23PWVmPzezseF2N7O/MrMnzKzRzL52TKilZvar8GfwZzObEXPO\n94Y/l6fM7A9mNjHcflv4M7/DzH4TbrvEzB4Lj33azK44nX8HkqLcXS+9EvoCLgIawuUC4JJweSqw\nJea4fOBjQF64fg9wXbh8N7ABKAnX7wM+HS5/EHgWyA3XPwI4UBauLwAuDJdvAH4RLl8DbAU+FK5n\nAZuAt4brJUA98L0TXJOFn9Mfz7/GnOca4PGYY68BuoGLwvWngeeAHGAycAjIJfil7iXgr8LjxgKr\ngZvD9ZeBunD5SiAjXHbgWzE/3+eB28L1LwC7ganh+jeA74fLl4U/tzHh+nuBh8Pl24BXgHnhejGw\nDpgRrs8Kf3YFyf73pdfwvrIQSbw1QKuZ/Ra4C1h0kuM6gTxgiZllEnyhLovZ/zN3bwuXG8P9ECSF\nn7n74XD9LuDOmPftBT5nZhUEX9QdMfvygJ+Gy+cSfAk+CODubWb2O4LkcRR3dzPrBO41s1ygHPjV\nID+Dje7+53B5E9Do7l3AdjPrI+iqKwSq3P3n4WfsM7MfArcC/T+7n5vZnQSJry/m/P8evqfDzH4K\nLCRItAD3u/uWcHkN8PZw+W3h5/4hbNBlAKUx52xz96Xh8huACuAn4bEAh4HpBIlV0pSShiScu3cC\nV5vZG4F/BN4DvPsEh74L+BBwg7tHzeyrBL/R9+uOPS1Hultjl4ldDr/QlwHvdPeHzexc4Ecxx7a4\ne38Vz2PPwwnW+897OfBvwJXuvsnMPgpUnejYUOwXvJ9gPfZajvs4AHf/TzO7j6AV8KKZvdHdd5zg\n/H3HvP9kP7dM4Hfu/vcnibk5ZjkT2Ozu15zkWElTGtOQhDOzUjO73t2fJvjCWxju6gbGhMcYQXfM\nq2HCKCX4jTg7jo/4PfBhMysI1/+ZI1+cOeHrj+H6hwc553pgr5ndGsY0hSDBnchYIAq8En7uu2PO\ne+x1xWsDsMvM3hm+dwzBjQL3h+tvA3a7+78SdFvVxbz3o+ExuQQ/45O15mI9ALyvf4zDzCaa2X+d\n5NingYlm9vb+zzGzH5tZ8Wlcn6QgJQ1JhhzgdjN7FlgKfCrc3gSsNrMXgPcD/w24ma0Cfgl8n6Dv\nfFDu/j/AQ8Afw89oA1rCffuBjwNPmNlTBN0zY83suP8L7t5LkKj+xsyeDz//lyf52GXAC8CLwMME\nXUf9sa4GssNYbjxV/Md8/luA94axPgb81N3vDeOdCTxpZk8TdOUtjXn7QTNbCqwCHnP3u+P4vKeB\nvwP+18yeIehe+++THLsXeDPwCTP7I7ASeM7d2+O9PklNdqQlLiLpwMwcKHf33cmORdKPWhoiIhI3\ntTRERCRuammIiEjclDRERCRuafGcRllZmU+dOjXZYYiIpJQXXnhht7uXn8570iJpTJ06lVWrViU7\nDBGRlGJmW0/3PeqeEhGRuClpiIhI3JQ0REQkbimdNMzsJjO7s71dlQtERBIhpZOGuz/k7rcXF6tG\nmohIIqR00hARkcRS0hARkbgpaciQ+dnTm3ngxZ3s6+w+9cEikpJS+uE+M7sJuGnGjBnJDmXUc3d+\n+cdtbNx1gJzMDK6cWcb82gg3zKlkXEFOssMTkSGSFlVu6+rqXE+EJ19fn/Pi9r0srm9icUOUnXsP\nkZVhXD6jjAW1EW6cU0npmNxkhykiITN7wd3rTn1kzHuUNGQ4uDsv72hncUOUxQ1NbG3tIMPgsnNK\nWTC3ink1lVQU5SU7TJFRTUlDRiR3Z03TPpY0RHm4volXWw5iBpdMKWHB3AjzayNUFecnO0yRUUdJ\nQ0Y8d2fjrgMsqm9iSUOUddH9ALyuehwLa6uYXxthcklBkqMUGR1GXdKIGQj/yMaNG5MdjpyBV1oO\nsKQhyqL6Jhpf2wfA3InFLJgbYUFtFdPKCpMcoUj6GnVJo59aGulhW2sHixuaWNQQZfX2vQDMjhSx\ncG4VC+dGmFFRlOQIRdKLkoakjZ17D7GkIcri+iZWbd0DwMyKMSyojbBgbhWzI0WYWZKjFEltoy5p\nqHtqdGje18nSxqAL6/nNbfQ5TCsrDBJIbRW1E8cqgYicgVGXNPqppTF67D5wmGWNzSxuaOKZV1rp\n7XMmjc9n4dxgEP3CSePIyFACEYmHkoaMKnsOdrF8bTOL65t4atNuunudquI85tdGWDi3iourxyuB\niAxi1CUNdU9Jv/ZD3Ty6tplF9VGe2NhCV08f5UW5zK+JsGBuhEunlpCVqVJrIrFGXdLop5aGxDpw\nuIcV63axuL6Jx9bvorO7j5LCHObVVLKgtoo3TC8lWwlERElD5FgdXT2sXN/CooYoK9Y2c7Crl+L8\nbG6YU8nCuRHeOKOM3KzMZIcpkhSjLmmoe0pOR2d3L09u3M3i+iaWr21mf2cPRblZXD+nkvm1Ea6e\nVU5ethKIjB5plTTMbDzwJWCKu79lsGPV0pDT1dXTx9OvBAlk2Zpm9nZ0U5CTybWzK1hQW8WbZpdT\nkJPSMweInFK6JY1sIBf4nrvfdpJj1NKQs9bd28dzr7ayuCHK0oYorQe7yMvO4JpZFSyYG+Ha2RUU\n5WUnO0yRITeikoaZ3QLcClzm7tUx228FPgVkAo+7+9+f4jx3nyxp9FNLQ4ZKb5/z/OY2ljQEc4Ls\n2n+YnMwMrppVxoLaKq4/r5LiAiUQSQ8jLWlcDTQCDe4eCbdNAZYBlwL7gF8D9wK7gItj3v6gu28O\n33PCpGFmtwO3A1RXV1+8devWYbkOGb36+pw/b9vDovooSxqaeK29k+xM4/LpZSycG+GGORFKCjUr\noaSuEZU0Bj7ALBqTND5KMEbx2XD9WuAD7v6+Qd6vloYknbuzekc7i+ubWNTQxPa2Q2RmGJedU8KC\n2irm1UQoL9KshJJaziRpJHqkrxSIxqw3ARUJjkHktJkZF04ex4WTx/GZBbNpfG0fixuaWFwf5Z8f\naOBzDzZwydQSFtZGmF9bRaRYsxJKekp0S+NDwDR3/+dw/U0ELY33n81nqKUhyeLubGgOJpVa3NDE\nhuYDAFw8ZTwLaoNZCSeN16RSMjKlQvdUFfAIweD4fjP7BfCAu99/Np+hpCEjxaZdB1jS0MSi+ihr\nmoJJpS6YVMyCuVUsqI0wpVSTSsnIMeKTRrj+HoK7p7qAJ939U2f7GUoaMhJtbT3I4nBOkNU72gGY\nUzWWhXODLqwZFWOSHKGMdiMyaSSCkoaMdDv2dASTSjVEeSGcVGpW5RgW1FaxcG4VsyrHaE4QSTgl\nDZEUEG3vDLqwGqL8aUsb7nBOWeHAvOg1EzSplCSGkoZIitm1v3NgUqnnXm2jt8+ZXJLPwtoqFsyt\n4oJJxUogMmyUNERSWNvBLpavibKoPsrTm3bT0+dMKM5jfm0VC+dGuEiTSskQU9IQSRPtHd08sjZo\ngTyxYTddvX1UFOUyP5wX/dJpJWQqgchZUtIQSUP7O7vDSaWiPLZ+F4d7+igbk8ONNREW1Ea47BxN\nKiVnRklDJM0dPNzD4+tbWNzQxIp1u+jo6mVcQTY3zglmJXzjjDJyspRAJD5KGiKjSGd3Lys3tLCk\nIcoja5rZf7iHorwsbjivkgVzq7hyZpkmlZJBjbqkofk0RAKHe3p5etNuFtVHWb6mmfZD3RTmZHLt\neZUsrI1w9bmaVEqON+qSRj+1NESO6O7t49lXWlnc0MTSxmbawkml3nRuBQvmVnHt7ArG5CqBiJJG\nssMQGXF6evt4fksbi+ujLGmM0rL/MDlZGVw1s5yFcyNcd14lxfmaVGq0UtIQkZPqHZhUqoklDVGa\nwkmlrpgRzEp4w5xKxmtSqVFFSUNE4tLX57y0Yy9LGqIsqm9ix55gUqnLp5cyvzbCvJoIZWM0qVS6\nU9IQkdPm7jTsDCaVWlTfxJbWDjIMLp0WzEo4vzZC5VhNKpWOlDRE5Ky4O+ui+1lc38Tihigbdx3A\nDC6uHs+CuUECmTguP9lhyhBR0hCRIbWxeT+Lwy6sddH9AFwweRwLw3Im1aWalTCVjbqkoec0RBJn\n8+6DA/Oi1+8MJpWqmTCWheGshOeUa1KpVDPqkkY/tTREEmt7WzCp1KKGJl7ctheAcyuLWDA3wsK5\nVcys0KRSqUBJQ0QSrqn9UDArYX2UP20NJpWaXl7IgtoqFsyNMKdKk0qNVEoaIpJUu/Z1srQxmNb2\nuVdb6XOYUloQJJDaCOdrUqkRRUlDREaM1gOHWbammcUNUZ4JJ5WaOC6fBbURFsyN8LrJmlQq2UZd\n0tBAuEhq2NvRxfI1zSxpiPLkxmBSqcqxuQMtkLqpmlQqGUZd0uinloZI6tjX2c2KtbtYVN/Eyg0t\nA5NKzasJBtFfP62ELE0qlRBKGiKSUg4e7uGx9cGshCvW7eJQdy/jC7K5cU7QhXX5dE0qNZyUNEQk\nZR3qCiaVWtzQxKNrd3HgcA9j87K4fk4lC2uruEKTSg05JQ0RSQud3bGTSkXZ19nDmNwsrjuvggW1\nEa6eVUF+jhLI2Rp1SUMD4SLpr6unj2de2c2ShihLG6Ps6egmPzuTa2dXML82wrWzKyjUpFJnZNQl\njX5qaYiMDj29ffxxcxuL6oNZCXcfOExuVgZXzypnQTip1Ng8TSoVLyUNERk1evucVVvaWNwQZUlD\nlOi+TnIyM7hiZhkLaiPcMKeScQWaVGowoy5pqHtKRCCYVOrF7XsHSrrv3HuIrAzjDdNLWTi3ihvn\nVFKqSaWOM+qSRj+1NESkn7tTv7OdRfVRFjc0sTWcVOr100q5+aKJ/EXd5GSHOGIoaYiIxHB31jbt\nZ3FDEw+8tJPtbYdY/fkbKS7QuAecWdJI6VsOYrqnkh2KiIxAZsacCWOZM2Es5UW5fP7BRnr6+pId\nVkpL6Uct3f0hd7+9uLg42aGIiIwKKZ00REQksVI6aZjZTWZ2Z3t7e7JDEREZFUZk0jCzXDP7qZn9\n2MyuPNlx6p4SEUmsEZk0gPcCdwO3Ax862UFqaYiIJFbCkoaZ3WJm95rZtmO232pmz5vZC2b2zXDz\n+cAqD+4H7j3ZOdXSEBFJrES2NFqAvwYGnus3synAl4EbgDpgkpm9g+BW4EEfIDGz281slZmtamlp\nGb6oRURkQMKShruvdPfdx2yeD9zv7u1hq+JHwNuAl4GLzSyDkzxL4u53unudu9eVl5cPa+wiktr2\ndnSxLro/2WGkhWQ/3FcKRGPWm4AK4BfAfxKMZ9yVhLhEJMVF2ztZtiYop/7cq2309jnnlBdSpCq4\nZyXZSaMZmBazHgGa3b0D+EhyQhKRVPVqywGWNjaztDHKS9v3AjC9vJCPXnUO82oinD+pGDNLcpSp\nLdlJYxHwiJl93d33Ax8EHkhyTCKSItydxtf2DUzQtHHXAQDOn1TMp+edy7yaCDMqxiQ5yvSS1KTh\n7k1mdgfwhJl1AU+6+/3JjElERrbePudPW9pY2hhlWWMzO/ceIjPDuHRqCe95fTU31kSYMC4/2WGm\nrYQnDXePHLN+D3BPouMQkdTR2d3LM6/sZmlDM8vXNtN2sIucrAyumlnGx6+fyfXnVVJSqAmXEuG0\nk4aZZbt793AEIyLSb39nN4+tb2FpY5TH1+3iYFcvRblZXHteBfNqIlw9q1xzgyfBKX/iZvYXwCeA\ncUAPsAV4y/CGJSKj0e4Dh3lkTTCQ/fSmVrp6+ygbk8tbLpzIvJpKLp9eRk7WSC1kMTrEk6Y/BVxD\nMGh9M/CFYYxHREaZ7W0dA+MTq7a20ecwuSSf979hCvNqI1xUPZ7MDN3xNFLEkzSyCFoYBuwDLhzW\niEQkrbk7G5oPsLQxuOOp8bV9AMyOFPG3185kXk2E86qKdGvsCBVP0vgv4FbgEeA5jn4YT0TklPr6\nnJd27A0SRUOULa0dmMFF1eP57MLZzKuJMKW0MNlhShxOmTTc/S4AC9L+Q8Da4Q5KRFJfd28ff3y1\njSWNTSxf00zzvsNkZRhvmF7Kh688hxvnVFIxNi/ZYcppOmnSMLM5wFeBduBT7r4LWJ2owEQk9Rzq\n6mXlhhaWNUZ5dN0u2g91k5+dyTXnljOvJsKbZldQnK8yHqlssJbGd4GvATOArwMfSEhEp8HMbgJu\nmjFjRrJDERm12ju6eXRdcMfTyg0tdHb3Ma4gm+vPq2ReTSVXzSonLzsz2WHKEBksaWS4+3JguZm9\nK1EBnQ53fwh4qK6uTnWqRBKoeV8nyxqjLG1s5rlXW+npcyJj8/jLusnMq4lw6bQSsjJ1a2w6Gixp\nxM5n0TPcgYjIyLZ598GBO55e3BYUAzynvJCP9BcDnFhMhm6NTXuDJY0rzOw1glttS2KW3d0nJCQ6\nEUma/mKA/YliQ3NQDHDuxGI+deMs5tdGmFFRlOQoJdFOmjTcXYVcREaZ3j5n1Za2gfLiO/ceIsPg\n0mkl/MtNc7ixJsJEFQMc1eIq3GJm44Cx/evuvm2Qw0UkhRzu6eWZTa0sbYyyfE0zrWExwCtnlPHx\n62Zy3XkVlI7JTXaYMkLEU3vq+wRzeEcJu6eAq4Y5LhEZRgcO9/DYul1BMcD1LRw43MOY3CzeNLuC\neTWVXHNuBWNUDFBOIJ5/FRe7+6xhj0REhlXrgcM8sraZpY3NPLVpN109fZQW5nDTBVXcWBPh8uml\n5Gbp1lgZXDxJY62ZFYUz640oek5DZHA79nSwLByf+NOWoBjgpPH5vO+yKcyriXDxFBUDlNNj7n7i\nHWa/IuiKGg/MBp7t3+fu705IdHGqq6vzVatWJTsMkaRzdzbuOsDShihL10Rp2BkUAzy3soh5NZXM\nq40wp2qsigEKAGb2grvXnc57Bmtp/PAs4xGRBOjrc1bv2Dtwx9Pm3QcBuKh6HP+0ICgGOLVMxQBl\naAx2y+3KRAYiIvHrLwa4tDHKsjXRo4oBfvCKadw4p5JKFQOUYaDbI0RSxKGuXp7YGEx/+ujaoBhg\nXnYG18yqYF5tJdeeW0lxgYoByvBK6aShgXBJd+2HulmxrpmlDc2s3NDCoe5eivOzuS6cJ/uqmeXk\n5+iOJ0mckw6EDxxglg38NTAFuBvA3V8e9shOgwbCJZ3s2tfJ0jXNLGuM8uwrQTHAyrG53Dgnwvza\noBhgtooByhAY6oHwfj8DngcuBb4I3A9cf/rhicjJbIktBrh9L+4wrayQD195DvNqKrlg0jgVA5QR\nIZ6kMdHdv2Nmb3P3djNL6S4tkZHA3VnTtC+446khyvrm4DGo2olj+eT1s5hXG2FmxRjdGisjTjwJ\noM/MrgAws5lA9/CGFD+NaUgq6e1zXti6Z6BFsWNPUAywbmoJn/s/c7hxTiWTSwqSHabIoOIZ05gF\n3AXUAi8DH3H3DQmILW4a05CR6nBPL8+80sqysBjg7gNd5GRmcMXMMubVVHL9eZUqBihJM1xjGhXu\nfuUZxiQy6hw43MPj63extLGZx9bt4sDhHgpzMsNigBGuObecojzdGiupKZ6k8SkzW+XuncMejUiK\najvYxSNrgieyn4wpBvjmuVXMr41w+QwVA5T0EFfBQuAJM3sC6AJw988Oa1Rx0piGJNPOvYdY1hhl\nScORYoATx+Xz3tdPYV5NJXVTS1QMUNJOPEljXfgacdz9IeChurq6jyQ7FhkdNu3az5KGKEsbm6nf\n2Q7ArMox/N83zWBeTYSaCSoGKOntlEnD3X+eiEBERiJ3Z/WO9oE7nl5tCYoBvq56HJ8JiwFOUzFA\nGUXimblvM0GJdMI/D7v7nGGNKk7qnpLh0NPbx/Ob28JE0Ux0XydZGcZl55TygcuncsOcCJFiFQOU\n0Sme7qnYb+S3AyNmFj91T8lQ6ezu5YkNLSxtbObRdc3s7QiKAV41s5x/qD2X62arGKAIxNc91Ruz\nep+ZPQLcMXwhiSRG+6FuHlu3iyUN0YFigGPzsrj+vEpurIlw9SwVAxQ5VjzdU7fHrEaA/OEL5/So\ne0rOxOPrd3HXU5sHigFWFOXyjosnMq8mwmXnlKoYoMgg4umeqopZbgduHqZYjmJm44EvAVPc/S0n\nOkbdU3K6HnxpJ5/4zUtMGJfPh66YxrzaCBeqGKBI3OJJGne4ezeAmdUSJI5EOAD8E/C9kx2gloac\njv6EccnUEn72gUsoyFHtTZHTFU87/C4zm2RmHwU+Cfz2TD/MzG4xs3vNbNsx2281s+fN7AUz+yaA\nu3e7+4HBzufuD7n77cXFxWcakowSD7wYJIxLpylhiJyNeJLGZHffAZzv7h8ExpzF57UQTOiU07/B\nzKYAXwZuAOqASWb2jlOdyMxuN7NVZraqpaXlLEKSdPe7F3fwyXtf4vXTSvnZbZcqYYichXiSxiEz\n+zWw1MwmAHvP9MPcfaW77z5m83zgfndv96Dk7o+At8Vxrjvdvc7d68rLy880JElz97+wg0/eu5rL\nzinlp7ddoruhRM5SPL9y3QJMc/dGM8sjaCkMpVIgGrPeBFQM8WfIKHTfCzv49H2reeP0Mn78/jol\nDJEhEE9LYxaw1cxKga8A04c4hmaOThKRcBsA7n7bEH+ejAL3rtrOp+9bzRUzyvjJXylhiAyVeJLG\nt8LjPgP8HvjGEMewCHi7mRWF6x8EHhziz5BR5N4/becf73+ZK2YELYy8bCUMkaEST/dUPpAJuLs/\nYWYHhzIAd28yszsIyq93AU+6+/1D+Rkyevz6+W185rf1XDWrnDvfd7EShsgQiydp3Efwm//7zGwG\n8OzZfqi7R45Zvwe452zPK6Pbr57fxj/9tp6rZ5XzIyUMkWERT+2pbwLfNLMsd+8BRsQETCKxfvnH\nbXz2d/W86dxyfvBeJQyR4XLKMQ0ze52ZPQ3Um9mVZvaJBMQlErf/eW4rn/1dPdfOruCHamGIDKt4\nBsL/g+C226i7P0mCak+JxOMXz23lnx9o4LrZFfzgvRdpHm6RYRZP0sh09yaOTMSkym4yIvz3s1v4\n3AMNXH9eBf+lhCGSEPEMhL9sZj8BSs3sy8CfhzkmkVP6+TNb+JffN3LDnEq+/+6LyMlSOXORRIgn\nafwNwbMTh4BtwF3DGpHIKfzs6c188aE13Dinku8pYYgkVDxJ40vu/rlhj+QMqDT66PPTpzbzpT+s\nYV5NJd99lxKGSKLF8z+u28yuH/ZIzoBKo48uP3nyVb70hzXMr4mohSGSJPG0NG4FPmFmfUAXwZPh\nE4Y3LJGj/eTJV/nKw2tZUBvhO+96naZkFUmSeB7uq01EICInc+cTr3DHonW8eW4V337nhUoYIkl0\nyqRhZjkE5dBrgEbg+/3Tv4oMtx+tfIWvLl7Hm8+v4j//8kKylDBEkiqe/4H/A0wimOa1EvjFsEYk\nEvrB40HCuOmCCUoYIiNEPGMak9z91nB5sZk9PozxiADw/cc28Y2l63nLBRP41q0XKGGIjBDx/E/c\nZGaTAMysBFg/vCHJaNefMN56oRKGyEgTT0ujmqBY4VpgBtBsZs8S3EV1+bBGdwp6TiP9fPfRjXxz\n+Qbe/rqJ/PtfXEBmhqrWiIwk5u6DH2A25WT73H3rkEd0Burq6nzVqlXJDkPO0nce3ci3lm/g5tdN\n5BtKGCLDzsxecPe603lPPLfcjojEIOnt249s4NuPbOTmiybyjVuUMERGqni6p0SG1X8s38B/PrqR\nWy6exNffcb4ShsgIpqQhSePu/McjG/nOoxv5i4sn8TUlDJERT0lDksLd+dbyDXx3xSZurZvE124+\nnwwlDJERL6XvZTSzm8zszvb29mSHIqfB3fn3Zev57opNvPOSyUoYIikkpZOGqtymHnfnG0vX8/3H\nXuFdl07mjrfPVcIQSSHqnpKEcXe+vmQ9P1z5Cu9+fTVfeWutEoZIilHSkIRwd762ZB0/Wvkq73l9\nNV9WwhBJSSndPaUxjdTg7nx1cZAw3ntZNV95mxKGSKpK6aShMY3U8J1HN3HnE6/y/jdM4ctvrcVM\nCUMkVaV00pDUsKi+iUunlfDFt9QoYYikOCUNSYiSghwlDJE0kNJJQ2MaIiKJldJJQ2MaIiKJldJJ\nQ0REEiulk4a6p0REEiulk4a6p0REEiulk4ZaGiIiiZXSSUMtDRGRxBqRtafMrAj4DnAYeNbdf57k\nkEREhJHb0pgL/Ie7fwy44WQHqXtKRCSxEpY0zOwWM7vXzLYds/1WM3vezF4ws28CuPsz7v6ymV0I\nHDrZOdU9JSKSWIlsabQAfw3k9G8wsynAlwlaE3XAJDN7R7jvcuDjwN+d6GRmdruZrTKzVS0tLcMd\nu4iIkMCk4e4r3X33MZvnA/e7e7u7O/Aj4G1mVgG8B/iQux88yfnudPc6d68rLy8f3uBFRARI/phG\nKRCNWW8CKoAPAFXAnWb2g2QEJiIix0v23VPNwLSY9QjQ7O5fT1I8IiIyiGS3NBYBbw9vsQX4IPBg\nEuMREZFBJLWl4e5NZnYH8ISZdQFPuvv9yYxJREROLuFJw90jx6zfA9yT6DhEROT0Jbt7SkREUoiS\nhoiIxE2ZihYoAAAPQklEQVRJQ0RE4qakISIicVPSEBGRuClpiIhI3FI6aag0uohIYqV00lBpdBGR\nxEp27SkZBTIyjBXrdjH/209QXVJAdUkBU0oLqC4tZEpJARPH55OdmdK/v4iMGkoaMuw+93/O45E1\nu9jWdpDNuw+yckMLh3v6BvZnGEwYlx8kkpLCI0kl/LMoLzuJ0YtILCUNGXaXTy/j8ullA+t9fU7L\ngcNsbe1ga+tBtrV1sK2tg62tHSxtjNJ2sOuo948vyB5olUwpLWBySUG4XEhFUS4ZGZboSxIZtZQ0\nJOEyMozKsXlUjs3j0mklx+3f39kdJJLWDraGyWR7Wwcvbt/DH15+jT4/cmxuVsZAEqkuPfJndUkh\nk0vyyc3KTOCViaQ/JQ0ZcYrysqmZUEzNhONvcOju7WPnnkNBy6Stg22tB9naGrRUnn21lY6u3oFj\nzaBqbF6QVEqDlknsmMq4gpzjzi8ig1PSkJSSnZnB1LJCppYVHrfP3dl9oIttbUcSSX9rZcW6FnYf\n2HHU8WPzso4kkqNaKQVUFeeTqW4vkeOkdNIws5uAm2bMmJHsUGQEMDPKi3IpL8rl4inHd3t1dPUM\njJ1sC5PK1rYOGl9rZ2ljlJ6Yfq+czAwmjc8fSCLV4RjKlNICJo8vID9H3V4yOpm7n/qoEa6urs5X\nrVqV7DAkhfX09tHU3jmQVLa2HWR7TILZf7jnqOMrinKPv9srbK2UFOZgplaKjHxm9oK7153Oe1K6\npSEyVLIygwH1ySUFvPGYhqu7s6ejO0woBwe6vLa1dfD0pt3cv6/zqOPH5GbF3OFVcGRMpaSQCePy\nyNIzKZLClDRETsHMKCnMoaQwhwsnjztuf2d3L9tjbhvuTy4bd+1nxfpddMU8k5KZYUwceCblyPMo\n1SVB11dhrv5Lysimf6EiZykvO5OZlUXMrCw6bl9fnxPd1zlw2/DWmEH6h+ub2NvRfdTxZWNyYm4h\nLjzqVuLyolx1e0nSKWmIDKOMDGPCuHwmjMvnDdNLj9vffqg77O46eORur9YO/rRlDw+ufo3YIcf8\n7EyqS2K6u2IG6SeNLyAnS91eMvxSOmno7ilJdcX52cydVMzcScc/k3K4p5edew6Fz6PEdn8d5KlN\nLXR2H12Kpar4SLdXdTiG0j+mUpyvUiwyNHT3lEgKcnda9h8eeGJ+W1iOpT/BtB5TimVcQfZAl1d1\nST5TSgqDxFJaQGVRnkqxjFK6e0pklDAzKsbmUTE2j0umnrgUy/a2QwMPOm5tC8ZUVm/fy6L6Jnpj\nn0nJymDy+PzjnpifUhp0e+Vl65kUOUJJQyQNFeVlM2dCNnMmjD1uX3dvH6/tPXTkqfnwbq+trR08\nd0wpFoDI2LyBwfgjtxAHg/TjCrI1OD/KpHTS0JiGyOnLzswIn24/cSmW1oNdR+726h+kb+1g5YYW\ndu0/fNTxRXlZMbcOFx41OD9hnEqxpCONaYhI3Dq6etjedui4kvbb2jrYsaeD7t4j3yfZmcak8QVH\nPejYP0hfXVJAQU5K/86aFjSmISLDqiAni3MjRZwbOf6ZlN4+p6n90FFPzPffTvzitj3s7zy6FEt5\nUW4wON9/t1dMWZayMSrFMlIpaYjIkMjMCFoWk8YXcPkx+9yd9kPdA4Pyx5a0/+2LO486vjAnc+B5\nlOqYBx2nlAbdXpoeOHlSOmloTEMkNZgZ4wpyGFeQwwUnKcWyY8+R7q7+hPJKy0EeW99yXCmWCePy\nBm4bro6dhKu0kDEqxTKsNKYhIiNaX5/TvL/zuDlS+p9N2XNMKZaSwpwjtw2XxNztVVpAhUqxHEVj\nGiKSdjIyjKrifKqK87nsnONLsezr7B4ovxIMzgddX6u27OGh1UdPD5yXncHk8Se426u0gEnjNT1w\nPFI6aah7SkTG5mVTO7GY2onHl2Lp6ulj596Yu71ajzw1//SmVg51Hz098ITi/IFbhqtj6ntNKSmk\nuEClWEDdUyIySrk7LQcOD7RS+p+a708wuw8cXYqlOD/7yMONMc+mVJcWUDU2NUuxjLruKbU0RORM\nmRkVRXlUFOVRd4JSLAcO9xyZvTGmpH3DznaWNpxgeuCS/JhbiAuPeoI+nUqxjMiWhplNAO4AMoHf\nuftvBzteLQ0RSaSe3j5e2xtODxw+MR9bluXAMdMDV47NPfpur4E/CxmfxFIs6dTSmAl8BtgF3A0M\nmjRERBIpKzMjSAClBVxB2VH73J22g10x3V1HWitPbmyhed8xpVj6pwcuLYi5hTgYpK8qHnnTAycs\naZjZLcCtwGXuXh2z/VbgUwStisfd/e/dfWW4Lxf41SDnVPeUiIwoZkbpmFxKx+RyUfX44/Yf6upl\n+56YlknrQba2dbA+up9H1jYfVYolK8OYOD4/5hbiwqMeekzG9MAJ654ys6uBRqDB3SPhtinAMuBS\nYB/wa+Bed7/fzOYRdFE96u7/MNi51T0lIumgd2B64IPHl2NpPci+Y0qxlI3J5YY5lXz15rln9Hkj\nunsqpvUQu3k+cL+7t4f7fgR8wMwed/elZrYcuPdE5zOz24HbAaqrq090iIhISsnMMCaOy2fiuHwu\nn378/r0dXUcVidzW2sGk8fkJjTHZYxqlQDRmvQmoAGrN7J1AD/D7E73R3e8E7oSgpTHMcYqIJF1/\nKZbzJx1fiiVRkp00moFpMesRoDlslaxMTkgiInIyyR6WXwS83cz66yx/EHgwifGIiMggktrScPcm\nM7sDeMLMuoAn3f3+ZMYkIiInl/Ck0X/nVMz6PcA9iY5DREROX7K7p0REJIUoaYiISNyUNEREJG5K\nGiIiErcRWeX2dJlZC7A12XEAZcDuZAeRQKPtemH0XbOuN72d6+5Fpz7siGQ/3Dck3L082TEAmNmq\n063jkspG2/XC6LtmXW96M7PTLtqn7ikREYmbkoaIiMRNSWNo3ZnsABJstF0vjL5r1vWmt9O+3rQY\nCBcRkcRQS0NEROKmpCEiInFT0jhLZlZgZveY2dNmtsrMvhiz71Yze97MXjCzbyYzzqFkZleb2Qoz\ne9TMfmNm48PtaXm9/Syw3My+ELMtLa/ZzP7VzJ4xsz+Z2editqfr9abldcUKr/FZM3vSzO4Nv7su\nMLOVZvacmT3U/395UO6u11m8gL8FvhkuZwHPA5cAU4D1QDFgwG+AdyQ73iG43mLgcaAwXJ8G5Kbr\n9R5z7X8H/AL4QrieltcMvBm4K1zOBJ4Dzk/j603L6zrmGkuAVUB+uP4N4OPAWuDCcNtfA9891bnU\n0jh7W4D+SXozgb3htoH5zz34G/kR8LZkBDjEFgD1wF1m9hTwVnc/TPpeLwBmNofgGn8aszktr9nd\nHyb4AumXAXSSptdL+l7XAHdvA65w90PhpiyCv9M97v5SuO0nBL8wDCotnghPBDO7Fvj8CXa9E7jK\nzP4RmA78wN1bzOxk85+nhEGudwVwOXA90AE8bGYvc/L53lPGINf8PuAHwG0Ev5X2S+lrHuzftLtH\nzWwiwS2Zd7r7BjO7hRS+3kGk9N9jvNy908zygK8T9A40EHPd7t5lZqfMCUoacXL3FQRfmEcJ+3tf\ncPdfh+vfM7P9nGT+80TEOhQGud7bgYfdfU+4/iBwMSl+vTDoNX8VuMfdN5tZbNJI6Ws+2fUCmNk1\nwKeAT7r7+nBzSl/vINL1uo5iZpOAHwPfcffFZjadmORoZrlA16nOo+6ps1dL0C3VL5egxZGu858v\nA24IB9EMuAF4ifS9XoBrgYVm9gDwFeCdZvYPpOk1m9ls4JPAzTEJA9L0eknf6xoQtjDuBm5398UA\n7v4KMMbMasPD3gcsPtW51NI4e/8/cKeZfYxgbOMl4L/d/VA6zn/u7lvM7NsEySMbWO7uywHS8XoB\n3P31/cvhb+DXuPu/hevpeM0fJvjFZ1nwewEA33L336fj9bp7Uzpe1zGuB84DfhHzd7qCoMv1x2bW\nB7QCf3WqE+mJcBERiZu6p0REJG5KGiIiEjclDRERiZuShoiIxE1JQ0RE4qakISnLzK5MdgwnYmZ5\nZjas80yb2VQze+4M3zvFzO4zs8fN7I9m9m/xPAksAkoaktp+kewATuIy4G+SHcSJhA95/Z6gyOY1\nBLFmAm9JZlySOpQ0JCWFJegj4W/LF5rZuWG59pVm9suwJAJmttHMfmBmq83s9rCU+5/N7JPh/mvM\n7H4z+11YBvwHYQn0DDP7r7Bc9rP9T82a2d1m9nELysIXmdmHw/P9ycz+Mgzvi8B8M3s8fE80Ju6v\nmNlt4XKjmX3UzO4K1z8flqh+xszmn+Car7Og/P7DwMditheG1/WEmT1iZueE228Jj3/czP5gZu8g\nKEj3nLs/C+CBv3f33w7l34+ksWSX7NVLrzN9AVtilp8CrgqXPw78XbjcDVQD5QR1daYTlHrZGe6/\nBtgMjCcoi30fcDNQBNwaHnMVYcloglIMP4j53PeF5ysAno05590xx0Rjlr8C3BYurwPeHy5fDzxC\n8IvcGODPQHbM+wx4FZgcrr+T4Mu//5yfD5cvBB4Il7eG15VDUJk4F/g08Ilk/93plbov9WNKuqgF\nvhSWSMghmNcEoNXdtwGYWZsH9XYws9hSCH/0IwUYVwA1BHOG3GhmfxOeb13M8YvDYzOAqcByoI/g\nC/p05AEPhMsXElTQXRGzbyJBmX2AMuCAu28P12PHMy4EKsOqtRAkB4D9QCHQSzBXRA5BIrn6NOMU\nGaCkIaksO2a5Hviou683s3EEX+bxOt/MCty9g6BVcS/wfuCgu19lZtcTtCj69VcCPR94K/AGgpbJ\nE+F2J/iC7udmlk+QWOYD3zvBuV4mSATvd3c3s8uB12KO2w0Umtk0d99MMK8JMe/d7O4/DhPZG8Pt\nvwB+F37Gv7j7fjN7CPi8mV3p7k8CmNlnCVptv4zvxyWjmcY0JJWtsWDqytkERfa+b2YrgdPtn99F\nULTtWYKWyW+Bh4ErzGwJQXIYd4L3rSUoob0C+AKwJRxLWQtcamYPmFkm8G/AM8D/AqtPFIC7LyOY\nPe4ZCya3ejtB11r/fgc+APxveI3FMW+/A7jGzJ4AniTojgOoIkhUh4E3mVm1B5PwLAQ+Ho7/PEmQ\n8H4T909LRjUVLJRRLaxa+zF3f2eyYxlKZjaGoOVzPXCQIHE1ufvXkhqYpDx1T4mkp4MEZfofBHrC\ndSUMOWtqaYiISNw0piEiInFT0hARkbgpaYiISNyUNEREJG5KGiIiErf/B4eVModGsP/qAAAAAElF\nTkSuQmCC\n", "text/plain": [ "<matplotlib.figure.Figure at 0x175156e00f0>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "T = [15.0, -56.5, -56.5, -44.5, -2.5, -2.5, -58.5, -86.2]\n", "p = []\n", "plt.plot(T, p)\n", "#plt.semilogy(T, p)\n", "#plt.gca().invert_yaxis()\n", "#plt.title('')\n", "#plt.xlabel('')\n", "#plt.ylabel('')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 天気予報\n", "\n", "天気予報のノートブックの`weather.py`モジュールを使う。\n", "\n", "1. このモジュールには,重大な不具合を発見せよ。\n", "2. バクを修正するにはどのようにすればよいか。ヒントバグの該当箇所をf-string `f'city:06d'`で置き換える。\n", " * [f-string](https://docs.python.jp/3/reference/lexical_analysis.html#f-strings)\n", " * [formatspec](https://docs.python.jp/3/library/string.html#formatspec)\n", "3. 出身地に最も近い地点の天気予報を表示せよ。" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import weather" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "道北 稚内 の天気\n", " 北海道付近は、20日夜から21日にかけて、オホーツク海から北海道の\n", "南東海上に中心を移す高気圧に覆われますが、21日夜には日本海から気圧\n", "の谷が近づく見込みです。\n", "\n", " 宗谷地方の20日15時の天気は、晴れ又は曇りとなっています。\n", "\n", " 20日夜は、晴れのち曇りでしょう。\n", "\n", " 21日は、晴れの見込みです。\n", "\n", " 海の波の高さは、20日夜は1メートルでしょう。21日は1メートルの\n", "ち1.5メートルとやや高くなる見込みです。\n", "\n", "今日 (2017-08-20)晴のち曇 -- --\n", "明日 (2017-08-21)晴れ 19 25\n", "明後日(2017-08-22)曇時々晴 -- --\n" ] } ], "source": [ "weather.forecast(weather.loc['稚内'])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 海面水温\n", "\n", "1. 生まれた年月の海面水温偏差を図示せよ。ただし気候値は授業で作成した1961~1990年とする。\n", "2. 1961~1990年のNINO.3の海面水温5か月移動平均値において,エルニーニョに加えてラニーニャ期間も同定し,`plt.axvspan()`を用いてエルニーニョとラニーニャの期間を分かりやすく図示せよ。" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAASoAAADuCAYAAACONEo5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXeYG9XZ9n9nNOq72up1r4CNjW1wwTZuFION6byAaQkk\nL50A70cgEELABBKSAAlJSAIJhBICCS3UEFoMNmC6MSYYDAZsA+5lu7Srcr4/RiPNjGZGs+tdew26\nr2suSTPnnDkz0tx6nuc8RUgpKaGEEkroyVB29gRKKKGEEoqhRFQllFBCj0eJqEoooYQejxJRlVBC\nCT0eJaIqoYQSejxKRFVCCSX0eJSIqoQSSujxKBFVCSWU0ONRIqoSSiihx0Pd2RMooYQSuhbjZ5bL\npm1pT21X/jf+rJTy0G6e0najRFQllPA1Q9O2NL9+fHdPbY/a7f3abp5Ol6Ck+pVQQgmuEEJcJoRY\nLIR4VwhxpxAiYDl+kRDiTSHEUiHEpYb9BwohXsseu9faryMoEVUJJZTgCCFELVABTJNSjgMiwNGG\n49OAk4HpwCTgGCHERCFEGXAXcIKUchKwDriws/MoEVUJJXyzUSuEeNuwnW08KKXcLKW8Ukops+QT\nA/5raHIEcJeUsl1K2Q7ciUZk04DFUsovs+1uA47p7CRLNqoSSviaIY2gKRP22nyzlHJisUZCiPuA\n2cANwEeGQzXAa4bP64DJ2f3rLfvrvE7KipJEVUIJJRSFlPJUYDAwBTjdcGgDZgLqk93ntL9TKBFV\nCSWU4AghxD5CiNMBpJStwMdApaHJ48BpQgi/EMKHRmJPAK8Ck4UQfbPtzsi27RRKRFVCCSW4YQUw\nNWu/ehkYAtwuhHhJCNFHSvk2GjG9CbwOPCmlfFtKmQDOA54SQiwGBgG3dHYSJRvVLgYhRBUwEuiH\nZtgsF0LEYgfsf20m0UamLUEm0YZsayOTSCCT7YhAECUYRAkFEcEgSiiEEgzS+NLCq6SUTUAj0AR8\nBXwopazfeVdYQk+ClDIOnGNz6ABDm5uAm2z6vgBM6Ip5lIiqB0IIIYCBwEghxJ5lU6b8JrlxI+0b\nNiCCQfx1daiVlRrhhEJZ8gmiVtfkySgY0t4HAsj2dnQSk4k2MlkSi82adZ1MJHKfU/X1JDduRI3F\nqAwE3m1vb1/a1NS0FPgQzYD6pSwl2e/xSEuFpnRoZ0+jS1Eiqh4AIUQ/YFbZvvv+tX3dOkQggBIK\n4e9dh7+uN4G6OqJ7j8Vf1xtfRQyNx7oHUkrSDQ0kN2wYJzduHFe+YSPJjRtJbthAJpHA7/e3Al+m\nUqn7gXullJ9122RKKCGLElHtBGT9UfaPzZz5VPzjj1HCYUJ77EFw0CBiU6fir+uFEva8vNxlkKkU\nyU2bED6VQL9+hHbbDaHmfyKprVtpeW9ZpPWjj4bLjRuvySQS1wT69SM8fDiNCxfOBV6WUrbs8ImX\n8LVHiah2AIQQKnBcZO+x/0iuXw9CgJSgKPQ68UQCAwcglJ2/rtH8zjtsfuBBAJSyMjLNzeDzEZs2\nDbW6ioaXXiJd35A/BiTXrSO5bh3AvwGEEETGjaP13XenAG9LKb1Fx5ZQggtKRNWNEEJMEEI8A9QC\ntL63DCUaBSkJDBxAbMYM/FWVjv23Pfss8RUfI9Np+l14gUm6MSKTTJLatIn2jRtp/2otCAjU9cZX\nUYG/tgalrAzF7y863/LJk0knEtQ/+xzh4XvQsuRdSKdpXLTIfL4sSTmh9d13QVsBIjx8DxKfrDxU\nSvls0QmU0CXIIGjNBHf2NLoUJaLqYgghqhVFObW8vPyi6urq/lu3bg0DXHXVVVx88cUsW7aMAw44\ngPYvvmTtTTdRNnEiSjhMfMUKZDqNEg7ji0QI7z2W+MpPkek0/tpasEhcya1byTQ10/jaa7QsXYpa\nWcHhU6cxavJkWltbWbBgAcs++ABfNEqmtRXh92vqpJT4e/dGra4i0KcPgX79iK9YQbq1lcRnn+vS\nEYlPVmpG+bIy0lu2mM599NFH8/Of/5w//elP/Pa3v3W9H4lPVgI8EwqFNieTyZ9lMpl7pJTbuu6O\nl/BNgCgt4mw/hBAKMKuiouLCtra22XPnzk2fd955kYMOOgifz1fQ/p133uHhhx9myJAh3HDDDXz2\n2Wf85S9/YfTo0TQ0NLB27VpuuOEGVjY20ufcc/BFIqb+DYsWsfWxx4nFYpxzzjlcccUVVFVVOc4v\nlUrR0NDAli1bEEKwYsUKVq1axZIlS1iyZAkzZ85k6NChTJ48mVAoRHV1NQsWLGDlypVUVlYSDod5\n8803WbRoEZs3b6a1tRUAn89HOp3X7ITfj0wmi94vVVU/T6VSFwDPSCkz3u5yCV4xcHRMfv+hKZ7a\nfn/U8+94CaHZ2SgR1XZACDEwGAyeo6rqOf379w9eeOGF5aeccgrV1dWdGk9KyaJFixg+fDjRaJSK\nigoAfJWVlE2YQNWhc2j/6ivOqapm/vz5ANx6662ce+65XXZNXrBgwQLuuOMOotEojY2NPPTPf3L/\nvffS2NjI8uXLuf3223Nktnr1ai688EKeeOIJfNXVpLdu1QYRAhEKQSLxqpTyh1LKV3boRXyNUSKq\nEgAQQgwoLy+/Np1On3zaaaeJs88+Ozhu3LhOj/fxxx/zxRdf8Le//Y27774bgIsuuohJkybxrW99\nK9fuwAMP5MUXX2TOnDm8/fbbbNmyhZEjR7J8+fLtvaQdgra2NkKhvH+PiESQWUJTFOWLTCYzS0r5\nyc6a39cFJaL6hkMI0be8vHx+Op0+7bzzzvP98Ic/DNTWbl+CxHg8TiQSIThoMNGJ44kMHEb7p6vZ\n9uoCkvXbCFb34brLL+ass86ipaWF0047jQULFpjGyGQy3epb1dX485//zDnn2Dg7KwpCynVSyiOk\nlEt2/My+Hui/V6U854GZntrOH/PkLkFUO39NfBeAEKKurKzslnA4/Ol3vvOd73722Wfhm266yZGk\n2tvbue6665g2bRqJRKLg+IYNG3j44YepmXoQdfvPxN+nN30vvZDyQ6biH9yP6EH70fenl1M9739o\n27KOyy67jIceeoh4PJ6zIQ0YMIBTTz2VBQsW7FIkBXD22WcTHDKkYL/apw4pZV/gHUVRPhVCjN7h\nkyuhR6K06ucCIURNJBK5MhwOn/vtb39bueqqq4L9+vUztWloaOCuu+7iRzf+hvja1QwZMoRVq1YB\nMGHCBAKBAJlMhldeeYWjfnAFLSuWk2lrIzRkKIEhg/AlJLXfOgkAJZ7/31DSPmLT9qNq/DRC69L8\n7sUvOPvsPbj66qt59NFHCQTcs7o2Nzfzox/9iBNPPJH99tuPjz/+mJqaGnr16tW1N6mTaFyxgpkz\nZ7Lcp1B1wjE0/OtZ2tetJ7VWS2EkpRyGEO+XTZ9Ay6tLqksrhd9slIjKBkKISCQS+XE4HP6/k046\nSZk/f35o0KBBAKxfv56amhrWrVvH2JNPpenddwjvtgfxtasBciQ1adIklq7+AlVVUcrLUaJRysaP\no/fsc1BrawocPEVKk4qkWqiKf7XgYeqXvQFAS0sLd9xxB/PmzcNJolu8eDHTpk0DoKysjCuvvJKF\nCxciwiGmjp/AlClTCAaD1NbWMn78eMaNG0csFtv+G9cBBAIBnn/+eQ4//HDe+evf6Pf9I/BFQmSS\nKVo/3ciWe58h8dFqWl5dArBV+FVIpWPZIOoSvmEo2agsEELsX1ZWdv+sWbOqrr322vC//vUvFixY\nwKpVq1i1bRuprE+Rr7KS8IgRKMEgTW++icyqeKHBQ2nftJGykaOJTNyb+IqPKZs4Ef+gPgBIv8y+\naqvyIqkgksKRqERKEH9mMRtf+hdDxh7FtvqVbFu5hAEDBjBs2DAuvvhi3n33Xfbff38OOuggfv3r\nX3PJJZcggn5kWxK1uozQ7gMI7zeB0KgRtC3/kvZPv0CKNOnmJtpXf0Vy3QYi40ZTftBU1l772x2q\nSra1tdHrwAOJf/oxVYdNJDZtNGq1RppNS1ex4Y+Pkt5iSuZwrpTyTztsgrsg+uxVLU+7f5antjfu\n8/AuYaMqEVUWQohYeXn5bwOBwLybb745cuE999D02muER48kOm4ffL2rSG3YSP1TzyOEQrqhUXPQ\nDIVIbd1K+T4TqDz4EAK9tKSGGZ2QVJ2Y8vdZJynQiEqJKyhJYepjbK/EFfxNArUVMn7Y+NYLtH65\nkvY1X4KUZJrjXH755dTU1HD99dfTpEiqj5hEzeTBhAfV0toeJpVSSMX9yJRAJM3SXLqllebXXqd5\n0RtkWhMokTAi4OfXP7ySCy64oHtuuAVPP/00J/70WuIfvE/N8ftTeei+pKWfdGuCzXc/Q/PCd3Jt\no9HonS0tLf8npXR3kf+GokRUX1P4fL7Dw+HwPUceeWRZr169gn/4yx1Ex+9D+ZwD8NfmfaLav/iK\nbQ89Saq+gdt+/gumTZvG3rMOJrlpI8N++ksyfmmruhlJR/ucl6Z0klJSkMkq4hm/JBM2kplZ4spE\n0ojUFhoXvM22B7UVQBEOERjcn7aPPgUgMKg3o287A4CWeNBEVHYQqkSoGWRLPZl4O8n6BJv/+BAn\nzzmM+fPns9tuu3Xm1nYYH374IeOPPpLEJ5/iH9CL6IQ9ic6YBAK2PfAs8eWfc+Ckye1vvPFGfUtL\ny7ellM/tkIntQuhqohJCzAMuBlJouc+/k832qR+/EDgVkMC7wAVSyowQ4hq0gg66SHy/lPLPHbwc\n7RzfZKISQvSKxWJ/ikQic+bOnRu558EHiIwfS8VhB6NWVuTayUyGpmcX0bBgIX+6+WZOOeUUfvaz\nn/GLm28mPHw4dSecDBF/AbkUnM+i3omUyJGU9hnSIUiHMzly08eRfkkq0UjzwjdoWvAqmcZ8koLY\noTOITh2H2qs3za+8hb9/H0IjB+OPaAOnkz5kSjGRlNAlvew+nah8/rynebIhScO/F9P0wmvUXXAm\n667/TSfvdMcgpWT9+vV88cUXPPTQQ/z6T7cS3GMQ0f3GQlph69+fIrr33oQ/+rg1kUg82djYeL6U\ncusOmdwugK4kKiFENfAcMENKGRdC3Ah8IaX8Xfb4XsDNwFwpZVoI8RBa+p8nhBD3AD+TUn68vdf0\njSQqIYTw+XwnBYPB26ZPnx557rnn1PCE0VQcO4dgTR8TychMhi33P0xy/QaWv7iAzz77jCMvPB8y\nUHPS8QQjtQghTOSir94ZCQlAMYyb8dvf90w4Y5K4AJpeeYOt9z1sahcaNYTAsCGUTR2PWtc7f20G\ntc6oYhbcAw9EBZCK+2l+dQn1D/4btXcda197s9Oe951FY2Mjjz32GOf84lraPvsCMhlkMsWTTz7J\nv/71r7a//vWv8UQicWY6nX5kh06sh6IbJKpQNrUwQoibgeVSytsdjj8K/FFK+bwQ4kW0hIujgDXA\nxVLKzZ25pm8cUQkhYrFY7IHa2toZt99+e3TWrPwX2vdnPyCi9s3bijIZtjz4CMmNG/jwxQXsOXsW\nmXiC6PixVEyZQUBqOaMyap54dDKy2qiM0lMq7KwiSn8mRyLEtTjBjX+8i/j7y4kdegDRieNRB9eZ\nDN46ORmlL9trTwrHY/k5WMhNJllz1o9RKsoIDBlIasNm7vvt75k3b57rON0FKSX33nsvp5+uFUKp\nOOxgHvvBlZx++umt9fX1LzQ2Np76Tbdd1Y2qkcf/ba6ntrdOuG81YCSPP9upZ0KIEPBLIAh8z5q+\nRwhRCfwRWCGl/El23w1o0tX72QIRh0opT+7MNX2jHD6FEHuUlZUtO+GEEw5cvnx5NBjMp8KoOHI2\n5am+KCmNbJSkYNtjT9Hy/jL2HTiA3UaNIrjbUPpe8X9UHXBIjqQAlBSocYEad1b39HYZ1d4FQScp\n0KQco5pWd/53GXzrjVQdfTjBun6oTX6UVl92xbD4VyiSImfnyr132JRWn2lsiUp43F5kGppJbdpK\n9alHcur5Z9Pru0ewdeuO17aEEJx22mksW7YMpSxCw9MvcMjJ83j++ecjRx111OyysrKlQoihO3xi\nuy42SyknGjY7khoAPIoWRH6uDUmNBh4AfquTFICU8jIp5fvZjw+hVVLuFL4xROXz+Q6NRCJLbrrp\npoF33HFH8JZbbmH69OkA9L3ge/SeeCiQN2jXv7yQhlcWkmlv44O6Gvpd+QOqD5uLv9mfIySR0jZf\nHFNf0MhOb+eLK7n3Gb80kZcOjSiU3KaRhdCM7YbN2DcvQWVsxhK5/iIlchKdz/DZuBlhWnFM+ej1\nvdOoOvkoUms30Pzy2/S+5GTiH3xGryEDaWxs7ND30FUYM2YMzRs3Ezt0Gqn1mxgxYgRXXHFF6Gc/\n+9nQSCSyVAjhTfcpwRVZSepu4Gwp5b9tjvcCfgPMk1K+YdgvhBDXCSF0Y+9coNNhUV97h08hhAiH\nwz+MxWJXPfnkk+F9992XyN5jiC/TqlIPueRKArEapJpXydrXrWfLU4/jq6yg7n9PIzRgiCZlZV0N\nRdb4rSQ1d4F02ExSAL4ESNVilzK0MUlaWdcEK2HofY0rgrmxMKiVSXvnUeO59XGyLcxjGc4tVZlT\nEfOvCrGZM4hOm4jwgy+q0ufSU1h/8z/otc9Iak6dzcor/kDEko6muxEOh2n49yu88cYbTJkyhWnT\nprFt2zZl9OjRsWOPPfbJUCj047a2tptLBSm2CwejVT2612BuWAAcBJwEHA8MBR43HL9fSvlnIcR/\ngReFEM1AA3BWZyfxtbZRCSEi5eXl9w0cOHD2M888E4lEIvQbO5b2tWsBGHDNlYT9NWRUSMXy0mym\nPk77F18SGrEHasJnWpXTyUmHVAtJJE8IZhjdD4wkpI8tVXM7p7EKSMsoARlWEb3AaSyjemr1AdMN\n7xtv+Rutb2QzNwiYNnUat956K2PGjPE+gS5CJpOhqakplxrn888/Z9q0aW3btm1bnEgkjvwm5XLv\noI1ql/Cj+tqqfkKIQeXl5e/OnTt3zttvvx1pbW2lbvDgHEn1vfJiwv4akuUZE0kB+KIRygbvSaDZ\nhy+RV+90kpJqfrOSipLKS1xWKDlJTJja6gSYk9RS+c1pHH3zJTT7mK5qGseE/HvrZh0r9zlrnzOp\nhkmz4V6zoSnUnjWPquMOwN+7CiS8+uqrjB0/jsCgvjz99NNFv6OuhKIoOZICGDp0KC+99FIQOABF\nNGeX0b8RyEhBSyroadtV8LUkKiHEjEgksuyqq67a7R//+Ef4P//5DyPHjCHT0gJC0OfccyiLDKS9\nOm2SFkRSoDb68DcpOYKCPEHpJGWElVSE5VV/74tnCclCZjnCs0lp7kYwVjKykp6SzJOrYpN00zqe\nFTph6XYtJa7kCEsnK6GGqDh2NgN/8/8YcNW3CQzoBak0yS/Wc/jhh/OXv/zF+QRdjKG3/IrBf7nB\ntG/48OEsXbpUkJEA/xVCHLnDJlRCl+Jrp/r5fL65kUjk4UceeSQyc+ZM+s08hPr330YmElTMOZje\nUw9DqtBebZai1EZf7mG3U++sMBKNVX0CbNVFfRyjmmc3Zm6MLMHYkZjdcf2zlZh0krWDnbrpBN3z\n3hivaPS/UtUM9S++y4bfP5br89FHHzFixIjig28nht7yq7x7R1Jh1bmX5o499thjHHvssQAIIb6b\nyWTu7vYJ7UTUjqyVh99ztKe2f5185y6h+n2tjOlZknrkueeeCw8dOpSqocNJNm1DJhJExo6l776H\nkQxljccGg7G/SckRCtgbx6HQnpRRNS9yO2RSmiomVUhaitYWE2ONUlDGX0icViIytjces/Zxktqk\nam+wz49TGINo50yaSimUzZhAuiVJzeKVlJeXU1+/Y6rDf37hJY7HjjnmGH75y19y+eWXI6W8Swgx\nQ0p5xg6ZWAldgq+N6mckqfb2dgYO3xO1poZ0i+b7N+Cwb5E21PQUKW3pPrBNwRcHNRu5lCzXVv8y\nfpst+yDrK4TpcIZMJF2w6Q9zRi0cqxh00rFT1/TjxnZqq7YpSVAToKQNW/E6C+ax7Qz3FpKSfoPn\nfFaasiJ60FRWiST/+7//y+TJkzs2iW7CZZddRn19PUcffTQ+n+90n8/3tXVfyCBIpP2etl0FXwvV\nTyepZ599NnzUT66nfvEiwkN2o/m/Sxly8Y8JRatzqldGzT+Qug0HIBXR4uyM6o2uRoBmv/LFFZNn\nufRLCKdND6tMKRD3uYbR+AqTfmaP2ZOLnTSVe02bPxv7ZHyFap+dVOWm/jmpe0Duuo1hN+mk5k2f\n+GQtwbuf5PPPP99pGUiH3PNLANRwknTSx+enXgHA888/z+zZswF+IKW8aadMrhtRPbKXPOTO4zy1\nfXDqn0qq346ATlKPPvpoePbZ55JqaqB8+lTqn9HqXbZvXE9guBabJlLgzzpn6qpQKmJW44xhLAJA\nTWdTo2gPoJGkpD+DYomPSwOSfIksq79UfuUPw373a3Qip2J9je2sdizjPjCTlNF9wghTeE8W+rWr\naoZUKi+gB3fvy1fxZpYuXcr2FL7YLmRDkFLkYxqH3vIrACorK6mvr79RCFEmpbxm50ywBK/YpVU/\nn883NxwOP3LSSSeFDzvtW8iYH//gfjS+/DKREaMACPbRUgcrSQg05AkqGdPUvGS5zLkoZCJpRHkS\nJaRJSSZJyZ9V9YyZDexCYVIKdpkTdJIyrsrpqpsd2RjbWNU643G7fnZwa6uvGhrn2hGo2fukZknb\n508jhCC89x785z//6dBYXYlV514KYe2PZvUZlwGaLSsTSbM266YCzPf5fIfstEmW4Am7LFH5fL65\noVDony0tLeE77riD6P7jaH1tCc0vv06mpRUyaQZdeCn+WGXuoU5FoK0a2urStPVOkurdjqxrg+p2\nCKfxl7fnHjSjlCRUTcWTsZStETmd9JFO+kjF/RD35QOQjZuFpMBZqjK9N9ibrPYrV2nKIH25kZTp\nXBYbVTGbmkwpOVVPh6pmcoTlr45y7TMPuY7R3Vh1+uUIVeYkKYDVZ1xGOBzmsce01UkRDT73dbZZ\nfR2wSxKVEGL/YDD4z9bW1lDs0Fn0u+lH1P/jmdzxuu9+h/7fOYdgn34mG1SyPOuBHk7npCYn2OVw\ngrwUJZJa4HAm4cttxH320pSTp7rFXpQjoiw5qQl3w3pXQXc2Nfp56bBmG9UDpvObYlL5Av4UAX8K\nVc1QdfA+tC75kDfeeIOdirjPNmvE0UcfTeXxc0g3xSEUeEEIsd9OmF2XIyMFrSm/p21XwS5HVEKI\nIaFQ6Kk2CPU68ztUHX0oTc+9DEDvS05mwI3XEB0zRlPnwjIXi5csl1kbVCbn96NLTTL7oBklI514\njJApkTOUi1Q+00Buiyuawd3geW71pyoGq2rnphba7Xca004i06ETptGo7uwzZd5Au29GsgKNsNTq\nGDXHTud3v/td8QvvRqw691ICW322x7Y9pP3BZVoTqNXli4UQ/Xfk3Erwhl3KmC6EiIbD4UWJTKas\n5qijiI4dTeMLL9P03MtERw2gcuqeJDZGyISzT2KjSjokClQYoxRgzCEuyWYsMAXo5n/gSlaly8Eh\n6Nfd49v8ajrmsILnVaLSfa7c+rk5fhpjEc2xfoUrfdJATDKloCvKAcP4AX+KyJhhPPizv3F7a+sO\nD1o24uMfXVy0Tc3hE2l8fcWXQoiwngiuhJ6BXYaosmkj1sbj8VjdKacRnbg3bZ+tYesjjwMw6tpj\n2ZpSIGz2OLeSlJZjKZOr/mIkH51ocitglv1g9kg3fjbCyWPcety0L124z3wd7se9wur4aXTb0D7n\nScqayM+oKtupzemkj3ZVJeDP35TAgFoiew1h7ty5vPTSSz2yWGo8HiccDhP/chuB3pX4+9fFhRBK\nKetCz8Euo/pF9t47I6WM1Rx1LOV770Nq7SbW/ua3AIz/zQmkI+Ukm/JFOaVNniXNGTOTK6qg54ny\nJWxICnJBv3qsnp53Sod1Bc/O2F3MfcC4kqf3Mb7qLhS6P5Rxs4O+v6C9jU9VYV/7jAmASV22pis2\nwqoCCkWhzwXH8OYXK/jDH/7gfPKdiFBICx1ofHEpFSfPpW3VeipPPMLZgFnCDscuQVRCiPmt772H\nWlVN9eQZJNeu54ubfgHA7j87GYbvwdatZUCWoBpV/Fv8JjUtHUvl7Cr6yp1ikI6MXuTFYBenp8ON\npKzGcmsfKJR2ciRj2axtjciRlbGPDdEZpSn79C7uz6q22GDe0kkfrfEA7cn8TRKKwqDLjuPSn19N\n3TGTbMvc72wsXLgQgDXf+xU1555O078XsquuBJaM6TsBQogRwDUAfc84hw0P/Z01v7sRgAE//ha+\nPfdk69ayrP+Sgn+LH3+ToTR6UpAJmx0VjYnmMi5BxVZ3glz/IvmenKQstxg9LzCRmA1ZGSWmjM9h\nfxGJzAo7XzE36GSVSikmsgr2r2HEr08n/tkGBs2b3qExdwRmzpxJsP9AANre+4za009FRKIvCCGG\n7eSp7VQIIY4XQjwohFjjcPwlyzYpu/9AIcRrQog3hRD3CiECdv29okcTlRCiQolGPgKom3cKX/7u\n1zQufQsUhdhhMwiOGUFb3J8LWwls9RWEp6QNVV0Ak7uBLUFZfJ6cYMxPZbRFGSUm62obYOst7gTv\nKpsziVlJSocxn5ZdsVQ3fzGr75QTUtlFC520ZKScgZcdT+M7n/GnP/W8YseJL9fgKytny98fItRv\nEFWzDibQp8+nQoiynT23nYhNwPmAE9EEpZQHGLY3s/frLuAEKeUktFqAF27PJHqsMV0IIQJ7DKlv\n/2QVoSFDafngfSJ9B6MM7kVy2xYqj59LKq7kVu2CG302huJCHyDdiK4fN6qHOkkZYwC1dhjaFBrL\njTCSk13mAzvYHfci8WRc+MJKTlYXhI5CphRXvzNrW52PdXuWZrtSCVREGTz/FC744aWoqsoZZ/Ss\nJAbJxgYqx09m2wOPUnf8SbR99RX+uj5N31TjupRyIWC7CCKEUIFKIcSDQB9gIZr2Mw1YLKX8Mtv0\nNuAe4FcFg3hEj5WoFEU5o/2TVQAkVn1OYtUqIqNH0fre+1R/6xjtKY1r1VLURo2k0mHzg2jKKR73\n5UjKrriLN/DJAAAgAElEQVQC2Kt7xVQlO+O5J38pF7KzGsCtx4rBK0l5yUFlhNUeVawtmKUwXbIS\ntXUM/Ml3Oe+qy7niiis6NoluhhCCta+8SHz158Q//pi6I46jff06FEXxFuXbA5CRgkTK72kDaoUQ\nbxu2sztwqjLgJeBs4ACgL3AmUAOsN7RbB9RtzzX1SKISQvTOZDK365/DA4bgL6+k/s3F9D7vLALh\nWnyNaq4yC2he59bAWmN5KGsVFz20BfKOmUZ1TycLnfz0LJzGzUhQaqsmTTmRhA4vTpwZH6RC3knJ\nTr3z0teYGnl74FQm3s7QrquCgf61DP75mdx0x63Unngg/X78HXqKwBKNRqkaO5mWpUtR/H7qjpuH\niJY9lK0a/HVD0XJZTpBS1kspz8u+ZoB/opXE2oCZmPpk93UaPZKoovvsvV4py5sF4l+uom7/wxl0\nwaWEq/rnvb8NKpxxtc5IWD6DG4Ix9g4KPcdz/bMuAcb8VWDw3HZSBW1sQ8ZVOmMf62pe7rwh7dz6\nZifROUl5dmRVTOUzkpW5FFfeIdYJxvqDhaE1dptCKu7PEZavPEL/H55MujnO5nueJtCrijvvvJN0\nuohT2Q7A/HmHkWqoxxeHyIChlI3dm/JxE7fs7Hn1JAgh+gghfiTyeuGhaCWxXgUmCyH6ZvefATy+\nPefqcTYqn893jCgrI9OcL3a720VX449VklbRcnZgNgLbwerEaX1vhDU3ekftOFbHTp2ETG18eYnL\nTq1LkScm4/k768xjNZzbttmOb1/akVpRdwa9T96GFRral9AZfUm2Q+t7Kznv+vl87+afs+HVd4jF\nYp2f4HZi1KhRtNdrvCRSUDvrMFbfciM+n+/QdDr9TJHuX2sIIV5CK5W1AU39W5ItibUUrdJyUghx\nHvCUEKINWAlcu13n7CniNoAQotIXq9iWbmwAoHzyFPocNQ9f3CwZGL2ojdCN4VanTacc5WCWpLyU\nvjKSmhV6OhY725ITjNKPdRXOmL/KS1iOcTywEF4H0g0DuKUc1mGtJ1iMqEx9DbnWjUglYPOfH+Fb\noydz2223eR6vq9HW1kY4EmXPS2+gvUYj2MQHK9jw+IOkGrbFpJRNO21yRVA2vI/c+w+ne2q7ePYN\npcR5HUVswuRtbRvXkW5soO//XUSk3xBEwsbOY/PA2ZGUDq8SklVKs9puivlPZfw5gc8zUeltrSSl\n7ZeAdl1uRGt3LitJGV/N47tDJBWbSsz2FgO7/U7kpUlXCliISvh81Jx+FH+58hYOO+wwjjjiCK67\n7jpmzJjBjBkz8Pt3jJNiMBjEF4qQbNgKNTVk/JLIiBFEhu2B4gs0Aj0vFigLKQWJVI96tLcbPcZG\nJYSY1bpyBW1frCGy91hCgwcXSDMFpaos+Z70B91rfnKrNNUZOHl8F2tjPGbnz2T0ENedUq2bPm+3\nuXdEvZN2iQBtiKZDklO2RL3t+VLCtCqo+2cpkRC155/Isd86hQ8++IBrrrmGg484jOjgvqxZY+t3\n2C0oH7gnjSvey/nmpcKSXnOOonnFMoQQPc9r9WuMHkFUQoioWlX9Qq8jjqX6qCOp+/a3s8bdwrb6\ng2dMo2IkKaeHvOCcxVQ+F2nKaIsyFiOVqmaAT4fNq4ZOcJK6hE2corlfvniElbDcgqXtxtDvl9nZ\n066IgzR5qduRlb7CapeTSycs4wb2YTgypRAcOoTK4+cwfvYszj//fGIzxhAdtzvz5893v7DtwMCz\nL6Ru7rH0OvwYVq1axYsP3M6WtxaSam/J2zzLI/Se8z8EaupeFkKEigxZQhehR8iHoUFDm2U6RWTs\nGELZB8jXoj10dvXzdHuL6ZglNQlgcug0FnVwIymvS/VGo7eth7vhvR152KmzUJiSRr8Gu/15WNLN\nFJ17YVwf2NulvIbQuPmnGccu2O/QR6iS8v0n0b7qK+549glSG7Yw8KxZ/P2BJ5h6++2cddZZnubl\nhpaWFv74xz8SjUYZMGAAG598hHRjI5HhIxk6dCgAweo6Nv37CepOPFn7HlSIjhlL47IlVIwcF6cH\nq4BfJ+x0ohJC9BZ+P/2/e27uAdLr4emwMzLr+yFfcMEO1j5GuBnmDTPM2Yh0ojGpbo6rj/a5qozn\nNs3TQrbGh96U1cDmOtOWPpqPlyjoaz6/PUltD6QqHcnKKGXlcs47qIQAMrveWX3qkWy+9X7aP/2S\n1Tc/Rc3svTn/8ou5bv1rrLnqzu2a79VXX80f//UUvnCY9nXrSW7QXH16H38y4X3HEP/4E5oXv0Hb\nko1Ex4wlOmqv3G+h1yFHsvrO3yCEqJZSbt2uiZRQFDudqKom779ekia421AyyJzKpz/I6axwbXyI\nswpJwX4j7B74/LHsmA4kZe2bUzct41jJxYh09qE1EqSbZ3hBsjrLQ29djbNDngBEjrw6CqtK5yjx\nGEjGREK2di5pamunGlrb5l0e/NSedzrx/77Ppt/8jS3PvYe/dyVrf/0QA3wBvvxR51YG161bx2//\ncge9zv8uob4Dka0pWt9agk8NoqQF5buPITJmLwKD+7Hlrw+y6eEHCF94KUq55jIRqKqhbMRo1Gj5\nFnqYVCURJFPe4jF3FexUG5UQol/De29SdcDBuX2+hD2RZHKlrLQHOpOtCpMvXVVccigwSNsY3Z0k\nM+N8cmNthxTiRlJWuB2zbe9wP/L2p8LNK6zGcTfScR0na4cTKa1eYu5z0rwprVrok7//IHzVFVx8\n8cVUEyQybgTrfnEv5QfuS3t7e4fPP/K7pxHZZzTBQQO0+URUyqZNpnysVtpLSWn3KzJjPP4BfVGi\nUbY+/0wut3zGD7XTD2HbksUIIWo7dRNK8IydSlSVE6d/VTFuEmrWsU9JmlW+3ApeNgOC/m8vLcSl\nI29jMT+AOrFZNyu8EIKTcd4LrCt2+jUWGrLtV//AmRj0+2PKFOGSCK8jsBrA9Xl4JSk7tc8N1vsg\nkgKak6S3NnDH8iXcc889tLzzEVXz5pDaXE/dvDkdD79JZlBCEa1Ih/EP0GBPFEmB4ofKY2eRbmmm\n6b9LaV21kmSjVqbeX1lNbNQ+VO934KaOnbyEjmKnEZUQYmDjB0uomn5QgaSiv+ZWpAyFBYwPpPXB\nhMIHwfrQW7fOwkh0dqt0BdlFbcjNKkm5rZqZxrZpZyWTjhCJeWzF8bPX+TnN2fhqknSN35EDkQVq\netPvJ5eT2rSFQw89lODug9j20HOE9x5O65KPiE4Yw5Yt3iNcnv/FL2h5cwm0S9M1ad9Jds7Z7zA6\neRThccORQrLuwXv5/FfXUr/0DdJhqNr/YOqXvo4QonfH7kgJHcFOI6rYvvutqdh3Cr7K8tw+48Ov\nVy3ORNIdTt5mRWekCaN9yYs0psfK6ZsxplCHsa+TJOWmihnHd1KVOkskpvO4kJ71/EbYzdscP9jx\neRnH9NfV0uv8b1N32ZnIeALZ1k7jUy/T+5IzkG3tDL/Ie+D/vvvuiz9WScuiN/DFDYUqVH3xJSvh\nZ+dfc9qRKOEQ6ZYmlGCIjY8+oM2popLYPhOpnDJzfeFZdg4yUtCeVD1tuwp2ClEJIYY0f/AeVdMP\nyNl60uF8euBUOO/L44WkrNKVk6RVjLCcHkCvsAY96ygImu6kSrY9EuH2kJdJhXQ4p1Xt7iyc5pkn\nZYXQ8D2o+8FZ+Pv1pn9FFVvueYzYgTOof+wZqmbP9hzUXDllJq0rPszmILNkffVr9lKlVavvqERC\n1J73bQCiQ0dQM+dIpNSkr6qZs2h87y2EEP06feEluGKnEFX5+EmfV06ahqyNmg3l+qperv6eNJdV\n90AgHfGadoOX5X0jWdi1sxJUAdEY/JacVNmC87pIX3YSV2dgNx8325mT4V5vawe7udnuM6q3KS1d\nj0gqkPERGT+G9X5B++dfkq5vof8l3yfx6adUzZxJJlP8d/D0D/+P5NbNpvFFSuTUPyWJltY6rpFV\ncNhAep1yCuHhI9jy7JM0vP4KAGosRvnEyVRMnvZV0ZOW0CnscKISQkxvWvo25TNnFPzgRUpkpSkz\nSQk1U1ix2MUHx/NcDGNaS7CD2dPd6jpgByf10O7B9hL0a4TrymYX+EB1Si0rItFZSVR/tRKS8b7n\nychZhdVXA8tnTie0+26k6xvYct8DNL/6Bn1OO4Pk+vWUjR5NfX296/z22msvUo31pFtbcvMwvmb8\nWrZXtdGHaFQRSYWyiROomDyFYJ8BbHrqUWSWECsOPoCm999FCDHQ9aQldAo7nKgqJk59OTJsBKJW\nK0apl7DSoRdiMEtShQbe7ZWcCozdLvanzsJOPTOvxHXuGrxIYF5XDE32p9xrx34WXvyi3JxBjXCK\nDLAL9VECAaqOmUvf718IPh/bFjxP67vv0u/M81CrK+m9z1iampyTHPj9fiLDhtPy0XLT/oxfZhdz\n8iqgTqBSStoat9Ln+FPw1/RCppJkVPCHYkTHjaPq4Dk7LhjxG4QdSlRCCLV5xX+pOeoo2+P6Cl8x\n2MWg2ZUaLzh/EanBqMYV823KjWk1KLvYjwpq5dkQgheSMK/E6Y6RHbddORGXPq5bQHFn4HVe1iR+\n+qJKJpJGxlLIWArC+T+4wB6D6Puj/wdA2+YNWgaG/zmOQL9+1E2cQDLpnBv695dfQstHHxTsN0vE\n5BIvbn34Udb88qes/v0NlO01xrSaWzZ1Eo1L3kQIsXNjaKW5sIbbtqtgR8/0UH9lNeqAOkRKFHgs\nS38hyWhBqjYP1HauBJrOYfynttiQ7FQux1WvTszJLlDXiMaXX2XTX+8vIBW79m4OncWM3MbrcVpJ\n1OP5nGxJ27vaaDcn0+espJ3f8n9sgX596Hf596k88jCtrRDUnHgcwqdQd8yRjuc46KCDaP38UzI2\ntQZ1qcqIxtdeBUAJhti2aAFf3Hwj6XicjF8SHNgfXzgMcFDnr7rnQQgxL1v26h0hxK8sx060lMv6\nSgjx/7LHrhFCLDUc60g+dhN2KFGVjRzzZPmkSbnPVudEr6qQF5Kyz43k7uRpth110BtcdScGq8rn\nxWgOIFQ/sr2dttVrWH3eD2hbsQrZ3G5oV3ysYufQ5w/2Bu2OwI2s7MZ2UvX0/V7mo1+/f0hfRMSf\nk4aEolB9/LE0LnqF5cuX2/YdOHAg0VF7seWZfzkGfjetXM7GZ56g/tVF+Gt7EZs+jZqjjsFXHiO1\ndSsb/n4PMpNBpARlUyYRnTDu+aKT3kUghBgMXAccAkwEBgghcoUupJQP6KWygMOBLwG93sFQYJ6h\nlJbnfOxW7DCiEkLUtX6+kvCksZoHcNgQoW+xOen2KSdpqrvhGH7itDTv4qhY6ICasX3v1Ld8v0nU\nnfkdfLEYvlg5G/54O2t+cCWb/nJfscuwdY4t3Nw94XXoktT2OsqCsxuHXTvwZug3XY9hjmq/GqL7\njefee+917Lvm+WdpXraU5JYtBepz4/KlbP7nw3z/4APZ+szT9D7rTGr+538o229f+v/gEtSaGuIr\nPmb9rbchUymiE8fRuvxDhBBVXu9HD8ehwCNSyoZsubA/Acc4tL0M+KOUsiX7eRBwsRBiYbYIaadD\njXYYUdXOPmpDdK/RiFgwt68gLKOjMW2G4gLW/Y59OrDEbtvfqpp6cDNwW+HzIg2JlMBfXsWga65h\n0E+uoWzqfrS+vZTG199AptKexnGy4RndD4zztZUKO2D7sjqfGp1grbnElFTWZ8kmh5aS1GIBlVat\nPFomkd/ydRoNDptZsjJeU3LtBn7z+MOO862srCQ2fTqbH34IafDBkqok1dZKaI/dufrqqwn060uq\ntTF33FdWRu+zzsQXi5H49FM2P/JPlMoI4VEjqD7xmJ2WUUHKwoSEThvFy2V5Kn2VJeajAOM/6Fto\nxLU/8AJwS2evaYcQlRBCNC59g8j+E00/JB1O0lTBOIaHzLlEk9EI62F1zMF248WvyUpSbmMYw4Ds\nyKJYaJAOJRik9qTj8NVUse2vj7Dpd3d5inMrcO/o6J+CQzxiMVhtWnY56PUyZcJAXHrbHJFly53p\nQcr5zSn2MT+/2NwDSK7dwPr1zs7jm556CoSg/j//MV1faMQw4p99CsD03XYn3dZikrr8/XpRNmVf\nABIff8zq836AryxK8+K3PN2fHoBi5bK8lr46B7hfSpn7hqWUl0kp389+fAitlFansKMkqkkynSYw\ncjDg3cHRi9pneuAtJFUwnjGoWV9Fsqhi9l7txVQ08yqkdX5uK5HFYF3Ny4V4BLQK28l1GyCT98g3\nnbuLVuy82O3sFhisfmrm9hQUytDJCsyEphOYGhd5worni8nqHutOCxKhPXej7KApDNxnDIsWLbKd\nv6qqLH/6GZpefIX2dRty1yv8fjItrQDMnj2b+LLlBfehcu5sBt10PbWnnwpA08LFpBubEELs43jD\ndh08DRwrhNBj3f4X+9JXZwI5/VpouE4IUZHdNRetlFansEOCfYL9B74eHDUMxQ9QKBHZFhBQM0De\nRmVHAJ2BkVRy+zpQlEp/+PLG88LxzO+dihs4q7xWgrZbHa089jBSm7ZQNmsKwucztM3Yui+A/XUa\n4/ncUIyk3D6DIVOpaiYhuwo6dmmUcwSmfdLaG86n28/soVB59CwCA/pw0NFHUHHUbDbd9SCKYia1\noUOHUnXcEay/6Q9Ex49FCQVpfmMJFbNnAbBmzRqaF79JzbdPAAz3RCrUP/4USnk+bjXd3Ex02sR3\nc5PdRSGlXCeEuB5YJIRoB16WUj6il8ySUq4XQkwE6qWU6w39pBDiv8CL2VJaDUCn07J2O1EJIcIA\n0UMnO7bpqrAX8E5gZgIpVp5cZOdotIUUj0XUz2EsCaUXMPAyVzsJSX84ImNGaW0ciFBrb31G7Py2\nigcWu8FpVc4acVB43ExITvnjrWXJzN+CZdykKPB/y2c9FUQm7IXau4bNtz1AsK4XV15wIfPnzydf\nPxM23/MAjbc0csIJJzBixAjO+91tjBw5EoDaWntbsPRLGhe9mpNsAUhnaH3rPYQQqlEd2hUhpbwP\ns+2J7Cqf/v5ttBVBa78HgAe6Yg7dXtdPCHGB8PtvGXT7ddpnNePqG2W1UcmU6FKfKX0OHYFxvsZi\nm/q83Maz1q2DPFnZ2eG0/R5UXosB2XrMTkKyzW7gkisqJy06lscqdCGwk2rsQmbAnPfeDna1EyFP\nanbVd+xKnhlzmglVkklCatUGtt7/KEowxEfPPsugQYOcJ5LFmjVrGDpqJANuvMZEbgDJ9RtpWriY\niiPnsPGPdxCbewANjz5L8qsNU6WUrxUdvAsRHDpA9r3mQk9tV3/nh6W6fgDRyRNvycQ1Hb+jBKH1\n6R6SsiMQO6STvqxzYZZYVK2fwHw9XsfrKJzUQicp1Kv/mN2xjkq20q/lcnJLk6wft88BLwpkPLcC\nr9b9drNVACUlCgp2SFUACpIMSkohMKAffS46j4Zn/sPQkSM549RTOe+88xg3bpzj9Q4aNAgRDJJa\nvxF/X3P6KX+fOqq+pUVc9LniewC0ffIZkUljF7OLq389Ad1uTE+uW0/5nBmmfd6M5JkOE5vZa7mw\nv5WkVDVTsFlhJCCnsYuRlLFmnRus7hZ2JO1kmHfzVC963mKLGl2omhuhx9SZzpVLP51/NdZCNKqI\nWs4o85ZbRbRsRqM7ZF0ofD6q5syh/6WX8vDatUycOYPY/tP55JNPuPrqqwnvvjv9zjCv1qcbGml+\n/R1a3nyXxhcWks4a2o3Qv6PwXnuQWP5Jl9yrbzq6VaISQlSLYJDgsOJiteMYWVXRSzvIk4ZODFay\nMpKUHfT91jiozkpMRoJyI6uOOrYKVXryF3M+biVxZ4O+l9XDYmqlG6wGdt1+ZVfIFWzsVob9SlL7\nbCwQoh3TCl5Y08YAqJUVVMw5kOiUCTS++goj9tqLwOABlE+fysaHH+T+WQdwyimnaJ3SaRqfe5Ga\nmhoaM2maX3ub3t8/D19UC7I33sfgiKG0/+E+hBDlPbkE/K6A7lb9DgruPhThVxFq5x50L7YqI0np\nRKOqGdugSyeC8gJjX+PYOgFZySyd9BUl2e3xvHciK5NfWoGUZbPy5yC5uc2tMPZQOJJjZ8JynGom\nKilIO/xq05Y6kE7zcKrwo1bEqDp6LpWHHwI+H0pKIdirL6eddTYjR45k3LhxWrI8KRFCIKUkNnM/\n6h97mppTjzefS83gU1WCw/qT+PDzmcC/OnwTSsihW4mqbOrkh9QB3rzm3WxRXh5mI0npMJJVRwhK\n7+NFXbOTmKzk1Fky6oqFBKu7ghGdGdtOIoHsw1/gK2XO72QHoyTllFc+1zYpihbWUMiTlV3164zf\nPt2MyQVEzXcM9OtLr6OOY9JBs1i59F0GDx6cM6QLISjffyqb/vxXx/mERu9BcEifp9iRdiq5fX+A\nPRHdaqNqXvwGSqzctM/kP1QkjYjXWD8vqmFH0RGSKixJLkwb5K/VqbQ52BNHd/3girpGFDlvMd8p\na954u3xfdiEzRtgFCVtTOlvHzeWR8vgX7BjbaCDX6Lh9qDjwAIbuvjuHH344q1atIpXSJv/JFdeR\nrm+0jQ7w+dNE9xlK6/srvU2mBEd0G1EJIY4CCO89svCYneHX1l3BuxSUTvo859cplo9HVTP4/Omc\nKudmnzKqpvrmSkZFijDYGcs7EtPY1YZvpz8SJzcEuzQw1kypRlLpLJxSPxvH9VKUA4ye9873LjZj\nOgN/81NeScUZNmpP/H4/Tz75JNFoFBEMkG7IxwAKg/khukcf0tuaS/nUtxPdKVGdCuALa6WOjdKG\nXaknrU1+n2m1zea4ndSik1UqpdAW99MW9+dW3IzkZJSW3MjKOo9ifYzxZ8WqwxSz29jGBNpcs90c\n3D47jeVlbDCvJBYrNGFMQGgKdLaQVa4StUXCslMZXdVI3WfKhZSKIec7ZiDcXNoYn4/KeYfS95qL\nAbjhhhsACO42hLZPPtP6Z/v4/GnavtxE25oNhEcPBTiYEjqNbrNRiXBonowXJiNTWn0F2RtB97rW\nVvjsJKn88cJ0uXo/cCYRI93IlEKajq/k2RnQtfFEjqCgc8bjztqiOpoltFgaZ6/xgfnvrdAobbVd\nORnZ9bAacJauiqWC6Qi2Nz2NDn/vSgb87se8Nv8WXn31VdRYhUmikikF/Gk2/2MBza8tp+LwqZTN\n3OcewNmY1ZWQwvP3uKug267GSFL6CpJOMopJmjLm7LYngnxbazoPQ5UUa3ycRTow2pCs6EhKVq8+\nUV0Fu6IW3ZEmuLPj2flqFUuJbKo4ZKOeuW1u53Gt+OOQILHo9Vn6638ovlgZVfPmcuC84/jRYUeS\n2rLN1C+VgObXtGR9jQveoX1Njyn7t0uiWyQqIUQu6ZQpQDepOEobbsvbxuO6ROApy6eLO4NVmjKS\nldW1wYmY8qqsxd5WJIzEnDLFY1bTDpBIdzloup+zUMLSYf1u9awHYL4XXrKDWqsqG8N8FENtPrvx\n83PV74+zmmztp11DxhSREJmyNw3/Wsj1C5+h7ZPP831Tgvj7H+Y/x9to/3ITQghFSrnjv5yvAbpL\nojoOoM5SuVb/Ibqt9pjtUNYfUj6TQjHP9WI+V1YYJSVnFU9xlcwg/29dsKn5Byzfrkj6GJvEcD0B\nbil6HPPMW/KuAybblbWvdRz7ij6WFD02SQ31cZ36WGGU8O1SSRt99hS/JDhiGG0rV5NubaV9Rb4A\nTaY1RWTK3vmBtQSHnfd8/oaju56APQGCY4bldugPm53NwVhUAChCBA4xbh5CbrzE+dlkQHScj1s6\nlsJ5W43KxmSB3RsY7gX5XF3ei6Ea+9kfKx6gbD2n0zgFBvwccbjn+rKdg4f77ZhaOqWYfj+RcXvR\nvnIVIhSg/t//Id3QxNb7HwfFT3pzA6Cljxk2bFgcKFwCL8ETuoWognvufhWQc4zLrfSlCg2nxsKU\nVliN1I6ez2rG9L6jcYK6B7k9GRW6UVj9o5xgfej1RH3WbJ9eUCwLqNvD7rWkWGdVxo70s0qWTmPZ\nXWex7KLWJINucLrvXuxXRik7tOduVJ9+HOmmZpJfrOWrS35K84LF1D/8NAT8HHHEETQ1NTFnzhyf\nEGJP14G7CrLQDabYKnRPR7cQVdtHmoNbLg1K1jZlrUBshdfVMo0onKWu3HidCJcxjmscvyBguJiz\nqkVi8lJ30GmcYgHDrqt4rvmyum6lsTNkpY1V6Hdm28dEWO7Sl/7qLBU5+KUVkSbtFnuEKimbMYmy\nKeNBCJRYmTZWJkNwt8Hcd999LFu2jNGjRwfKy8vH206oh8OtXFb2+EXZ40uFEJca9h8ohHgte+xe\nIUSgs3PY6cYPt38+u7xKRp8qq9+VnUpnlbDsMifY9bOSlN0cOoJimR3s+zgXWPWS8riwjfu57cbv\nSGrmYjYrNxuUcYyOZHKwIx231DVuRAgUXHtBu+yfpBXlc2aSrm8g09isjaMIGp54nlgsRt++fRk5\nciR+v3+s44X1UBQrlyWEmAacDExHy4l+jBBiohCiDLgLOEFKOQmtKIS3JFk26HKiEtaMYmj/mt5L\nIymO/3j58YxhJ2aisSMdncTcMie4EUmBp7iHQgyeKh53QOKzm1sx0rL2tTu3PcHZP7R21+pEbF7u\nkRHGNsUI2om07VRh63432JG80/yNkrdQJWpVBYNvvTF/PJlk1apVuc977rknra2tw9j1UKxc1hHA\nXVLKdillO3AncDQwDVgspfwy2+42nMtsFUV3uCf019/k7DgdHKBYRkkNCoSds2eCfUqXgN9uyTF/\nG8ypgjOmH6MV1mISXbE6VyxbqDF9jZPq65UAzddnT2QavF2Xo92nWLI/VeZ+I8bzdiSGM+er5yLt\nuams3gzs2hgSZ/eYfj+5nGT9ZhLLV3D33XdzzTXXANCnTx+AgBCiVkq52dNF7RjUCiHeNnz+s6US\nTbFyWTXAa5bjkz306xC6g6jGAkQm7tM5FcmDL40XD2M9m4JOTH6HNDPJlPbg675T1vzmxcp3eSUr\nR+14lqkAACAASURBVI/7IqRilRCN+ba85upygzcV1PkeuI1nJXm3kl1O98ZcBKPYXAtDrdznaj8v\nt3bFyMpfV4vavxoySX75z7/niEoIwbBhw+IffPDBnsArnibYWUjv9l6y5bJcjm9Aq3isw1ouy6mc\nltcyW57QHTaqvgDVJ5qlPLfYK5OnehdmCwj4U/jVNOWhNipD8dwWUlMF7YCCLJ9eQmy6y7XAyeZW\neP4d5z9opx7n1UR79dKqquqvxnzzuvpu3fTjxlcrjO2t57HCbVXU+Zo79/36qipIb2s07Rs2bJgC\nDOzUgDsPxcplPQ6cJoTwCyF8wOnAE8CrwGQhRN9suzOwL7PlCd0hUUUjk/bBVxYtOGCM7dJhrf2m\n5wuyeizbSVNavqbCCejSlF9NUxmKE1GTRNR8WsjWlJ+tiQiJlP3lG7N8GtUtL3DKnLkjCWVnwEli\ndIPVvmhnOzRKtUY42RvTubm4q+qdQTH10jS/sgiZZnOa4srKSh9Qbt+jZ8JDuay3hRBPAG8CKeAf\n2ao0CCHOA54SQrQBK4FrOzuPLicqIUS5WlvleDwXid6BH400JDszrugorT6NuCySh67yhdRUjqRC\nPnP+2q1oal97UrsFulSlf9bH0ckqleqYWgGFJOVUNqujsHrLW+GmojmlZnaakx0BWdM9W8/rXMvQ\nmA8+47q4YdxvPY8xSaLZ5qh9d/oVFUQ2WFR161yl5TvuaOZTI4H5qitINzaTTCbx+7VcyVVVVX52\nMaICvJTLugm4yabfC8CErphDlxNV2ZyZPxWhoGmfkWjAW6oOvZ9pHNUc2yVSgkxKkE76CnJHGW1S\nVpLSYSQl43snOKZWcVUd7L3hnQjCS6EIr/O0lyoU1znZZSl1kpbspE1n14cMufKhDiRlXejQvxP9\nPNY+umpvhmFhxHYmzivGujRWDJ7iTBUFEfCzfv16Bg7UtL3Kykq/ECJWtPN2QsiuzTrRE9DlRCUT\nbSjVldp7A6k4QSelYjdWD5vIGTbJkMH9RxOyqHw6WlP+nNqXSilF0xQbnfvcDMJWeC3N1ZHSXXbn\ncPL5sktlbOdPpsOo6nbGUG8lIGvMpNu5rYse+iKHDmMfY1urvVGDy8+6mKoJyJR55dcNbvdIiYRy\nmUABYrEY4XC4xnXAEmzR9UTV2oaSTZZnJBWZtTMZCamjyc2M0etC1cjL+kMy/pATKT+tKXMpk9bs\nvmQqnxHULq+6Y3I8F9HfzZbiZXzrg+Mlg4N+XrsVNqdVNSeVy6jqdoSsjAQUCbfn9gf8mlRkXU21\nO7cVeUlJNf2ZGKUonaRCapJEym+zSKK9WqVlO0lMI0atnVsKbLtVYd0mZuyT3LiFdEMTgwcPzu0r\nLy8nEAhUu154CbbocqLKJNqIf/gxkdGjUcL2dbqVVD7eL5dLiEK1z2y0LLR/eDHeGolKk6Q0aUr/\n8RoryFjJw4kc9B9lMRXATVKyShxuD25n7FmdXa2yIys3GNUy/RqsUlExUirmQlKMpIyvXuAsiQGo\npLNqqpOKqn9vOVsdFJBV/L2PiOw7FkXJf8exWAyfz1fpeaIl5NDlRJVcu4HU5i3EZswgWFWYJrpY\nQn8lKUiH7UIx8n4rTsZc0zwM4nsiS1aJlJozoFtTEneUDDpDUlbJzamklxE7MkmfFXb32cnzX5dk\njWSjv3ey/9k53+oEkkip+NV0wWKHua1GTrp635ryeyKskJpybZeyZIA1Suka8lJeKmX+A5UpH0KV\npLfU890J00zjlpeXI4SoKDrBEgrQ5USV2rwFgEzCnIZYN6hbi00aYZS0oPNSQXtSNah/+QGNJOWE\n7qhoY4WdXcyLrWxnwc5oXtzbn4LjRsIy9rFKSR2B0QYZMRBXa+7PyV6qt4N+/ki4nXY1T5DlobZc\nG6M0Dub7oEmgEp8/TmLFpxx45v+Zxo/FYkgpu92YjiwuEOxq6Na6fm6pUKxuClZjujWbonHlyYuU\nYfcvbiSpYlWLPadfMaywufkSWYuU2pGlV9XT65z0eWn7DPYuD2O4kabVNcAoRVmN4MXgpO7pYzmR\noP4HZPWRA02yMpJVSE2Z/rCM742qo5HUzAb7pElq0+1ZeptUSksK4POnSQPpdRvINLdy7LHHmuZV\nXl5OOp0uc7zgEhzR9X5UkTCyNY5QVTLb4vjKovng0KQguw7YZeezSymsP0hWsupIbnQvZFVshU2H\nVcW0U5+6S8Wzm5dOWvosOlPkwkhkyZTPkXD8arrD5OXkiOt0Pp2YOnoOJwnO5NqS9cWzHneSqgBS\nDdsI9a822acAIpEImUwm3KGJlgB0y6pfHICNt96FbI0z6Pe/AH/eWdPN0dOtzpuT06FxlclNrbPu\nt7O/eKlwbBcfZqzs4uQgaVf23eoBb52PV7ilpHHap0Gbq1NJejdo9zMvVRjJSFe99QfeK1kVIyjj\nefTXzpCU9XwdMcTrMF6f0dQAKi1rNhMc1KugTzweR1GUwtJMJRRFt6h+/qH9Sa7fgq8iRmrLNi1I\nk2ysVU6qysMurEYkFVDTuQff6vWsqhmi4bacYTSR8tNEXgw3otjqWodCZFxyUxmlF6vkZLXzGOdR\nEAbiqpZ20LdJL4ZhLUCRKz1mr666SZ/6/PR6idb5F7NZQf7hdpLG7IjPro2urhkJy+qSUowArXYs\nq6RlHE9vq68cJlJqwSJCywerqZg2quA8jY2NKIrS7DqZLoIo2ajcERy+GxVHH8C2+54iuXYDzYvf\novHZBdSecxKRfccVlap0iKQg4xc52w9g+2CZbAyq6sm2YgdzlL6zr5R5jh0jDSeJJeBPmVQJp/jC\nYiRllfKs761pcranWo2xjmKKvLNnsQUB48qnkaycoN8XO5tjwG+1PXk3nNvBidCMq8ZOMJGbIlCC\nhXNpampCUZTGggMlFEWXL3EpgSCZlnaCI7QcYY3PLuD6669n6/1P0rrkPcA+TUvTR8to/eLTgv35\nunzmVx0Jl9Udp6wIXuDqcW5ImWtN+OYFnVnd6whJ2aEz+bHtahjm88tbqysrOQlLJxWrmmaU0vT3\n7UnVdXNDe1JzN2lKBEmkVMdNn0NHbWU6jOPYwSoRNr35CWtueKSgXVNTE1LKElF1Al2f4TMQYMud\nDxJ/f0Vu30UXXcQbC15i29+fouGFBba5rOuff4F1v/8DrV9+po1jqZdnJAHrQ6uTlNHbXCcpL2qI\nW4UZr+Wq3AKQi8HugTT2t8vdXqwEuylne8GKavFSXXY+ZqakgpY8525kZSQda1my7dmM904nIish\ndYag7FYI7cZxIq+qo6fZ7m9sbCSdTm+zPViCK7qeqDJavF9681Zth+ojGo0yceJEVr+/nObF77D5\n7n+QTidMklW/i7R0yutu+T2ZdnfjpkwJ2uJ+GprCbGmK5rbWeMdyxxudPm0N50UIyipJefGUd5Km\njA+f4/ks5NmZmn+mQqAWCdBaJsx6zA1GskrF/bTF/bTGA7TGA7lr06Ux67mczus2H3CWymD7JCiw\nl6KK2dL0z+nWJPPnzy9o29TURDKZ3NrpSX2D0fVEFQoR3mcv1AF9AIhOzWd56N+/P/WffIZMtLP1\n0ScATQ3M+CWEfEQnakU6Gl5ZmB/PJo+1SCpkEr7cA9EW93fI9QAwPThOamVnbThu0pR1nl4ISp9b\nnpzM5Y68klW+Qkvx63IrH+YY65hUClRB42as7mO3AQXt9a0YqenQpTk325YdihncvZAUQLw+SdOr\nyzj//PMLjjU2NhKPx7udqLTsCd62XQVdbkxvfv31FiA6+O5f0L5mLRt/dRfVpx7D1vseAyAajbJ2\n4cvUDBpI2/ovCfYZAEBy02Za3l4CQPuaL8jE4xAOOp0mlw42F7KQ3a+ThNcfqGk8S4krr0nSOgo7\nYrKuEtqlW9HmaV11FLaqdGfgREwdKkKh2+7Ip3XJj+WSZtqwaGLtYzxmnIsbWenGenBeNYRCgrKu\nQrr5iIE5bjCRUmn58EsCg/pRV1eYHry+vj6ZyWQaHAfbhZAt4nI9cBAQBG7M5q0ytvk5sJ9h177A\nNCnl0mziPSMuk1K+6XS+7nBP2Eg2x7KvuoKK4w6j/qGnWLZsGWPHatWCKioqKJu8Ly1vvQv7qbS8\n9CaNixZpE+pTR+23T2b1D66k9sxvEZ0+2uSd7gTjA1UsHMVOmrKDEkqbVhy9Si5OD5CTpGVnA9LO\n6Zysz5rx1Eq0HUFxQ33+3up5pSQZx3MaS5o5HQOzZOf2PRjdPrzmUHf6DbgRj9X4D+6e83bIZHwg\n7O/ntm3bUkBThwbsuTgF2AOYgpYM8HUhxAIp5Tq9gZTyCv29EGIC8CMp5dLsrqCU0khiruiOwLan\nfNWVyJSg5eV32HrnA5QfPJkZl19garT4N7fQ8tZStj34eI6kACoPOwQlqElSm+/4GzLTsYfPyVDr\nBtv6gYYH0+6BsFvpK/bAe7HLOI3lWJknZV/xdnsq4rgZ6T31dyAvq5uE1dbmWtDVZvHAqBLawel7\n3177lT6G3T6ZUqHNXqeqr6//OhHVEWgVa2R2JfNh4DCX9jcAlwEIIVSgUgjxoBBikRDiumy+dUd0\nB1E9rfauBUBEtLCm8Jg9SHxkdj0YNWoUH775FomVnxEZPwaAyuPmEpk6FumX1Jx0PABb73mM5Obm\nDjk62v14vZCWbmOxPWbILurmjrC9Qc3mSs0dI4vOGtjN5xdFPntYAXUgILNdzb6suNf5W+dVjPzt\n3B2cVgrd2ri1bU+qZBoToNh/Z6v+f3tfHiZHVbX/3qrqnu6ZyUzWyZA9hC0EEELYVyUIyKYBTRQR\nUVQUf3zqp/iBivp9LghuiIBEwiJIEpZgCBCWJARCIMtkY8lOQpbJNjOZme7pvarO74/q6r5Vfau6\numfJZKj3eeqZ6epb995e6u1zzz3nPR9/rANoLPriehaDGWMN3PFt/knG2GcYY4vtB0ooh8UYuwjA\nbiIySaAawGIA3wZwIYyCMDe5TbI7ln4b1L0HwBRCcMxwyANqEB5Xh8yeJot+NAA0NDRArumHgVO/\ngEHf+TKYnP8CVJ8+CW2vvI7Yu2vQ8eYKHPG7HyEwZCDA8tebSw9RPpubQ1sYTOnBz8MUKiiTZLe2\nTKIRnRe1d4JT8GZX+aPMMbwnX7un6ZgQkUyxGC7+efP1eV3KOvnQnHIqRSJ6buft14rO82qx4Unj\ncXDmc9i1a1dOghgAiAhbt24NAdhQ9EV1FlSSo9y1XBYRLQKwyH6eMfYECsth7XDo5icAfsH12Qbg\nu1xfcwBcA6O4qRDdYVHt0uMJUDqO4PB6DL/nZ6B0BkyRceT0uwAABw8eRL9zz0JtbS0olYaWiVtI\nCgCYomDEr36G+t/+BCwYwP57HsbOb/0caos1DEXkD3FK+uV/bU3/lMhnYu5a2XekjHkVr0VXbsUZ\nkU/KuQS5TWKYs1Dsh/uYTPh/sTZ5K8nZcvIyvh3lBKZ2xop1Cix1Cj61t+f9WlpGhlQhIThyKLZs\n2WLpr6mpCZqmaQCayp5s78JcGCWwwBirBDAFwHx7I8bYWAAjiGgld66eMXYHV1X9UgCr3QbrcqIi\nIl2pH4y2l/J+J7U1Cj2exIF7HsLBgwfx/e9/H2q0DVdO/SIqzzwFB594rjDDP0BgkoRg9QCMfOA3\nqPvBNwEAe356F9pefh1aJL/Uz5OVWMJF5KR2Iql8G2e54XLJyO0aYYHTonFc1CUWlpclJt/GbV48\n0TCVeUuXsrXzGnZhJVBxqIOTnI6XYFKnawAUkFU6o+St5lAF2tutm3sbNmxAVVXV9mxZ9L6A5wDs\nyVZZfhPAXdnSWiczxmZx7a4B8Irt2v0wln+rGWNLYMipTIcLuiUpmQWCiMxdjAFfuNQYZFANas4/\nCfGNjRgz7SqoB5rQ78JTkRpVBz2hQW06iPia91F5yom2kID8r3ZF3TAM/83t2HvX39A+7zWkd+1G\n3XduBGDfPbIqAvBw+uW1RnFz5AWrL8qtCKYJryWnnFBoseSVJwpqBubSeAqtq+6AV4ICYCUerk6j\nCF4qXzte67B8te8C8ygWvCqSpRbBXAryqh1M0SHXVKOpyWo4bdy4EZqmvV+008MEWcL9b8H5tQCm\ncY9FZbQIwB3ZwxO6Y+kHSTH8SLkquNWVGHLLF1H/85uQ/HAzEus2I1AbQv9LT0N8+RpUnX4Kmu79\nN7RIK8wqurzz2rwJ5Pr+UIYOgjywPxJr1yMTaRFucRfKtYgd1E6Vc024bbM7QVTh2EslE3Neot0x\nx+sc5l6upeXkjBefK395V+qcisExjcghUNRpWW+iGJHZ07TMa3IuAx0FWlQffPBBOhKJrCr6YroC\nfsCnNyQ3bf07gFw8gvnBKwNrMeovt4JaDyIwxNC4r73qAsRXfggWCmLPHfdh+D235arYAMhKlHA5\nZkyCdrANAKC2RqAMshY7zX9hnYMN+YRiEZxUBkp13HqBaBlqXz4VszjcLC1LX9ngUCfpF1Ef1nPl\nE5KX1+E2F2+OddH8nPM43drwmumW89x3McjlwVsJjwqIau3atQki2iiYoA8P6BaLCrYa8xYZ3IwM\ntTUKKRFBMKBiwJVnouLo0WCSBL0jjvaXFufa5oL8OOug7rbv4Ij/MyzOA3c9ALWlycH0Fyftlnpz\nFvYr/jUWVfO1PxZZW9bxWUkk1ZkwBMe4LAcLqTutJjtE1qAo5EFkAdrPiRK43cI43Kxnz/FaGpD3\nExvYtGmTgp7Y8euj6C6iWgJFAakq2l9egvTOPZabe+ftD2PXH58HADBJwsCvXIYxM27HkG9fjepz\n87mBwshmSULgiDrU/8ow2DK79xe2cbwJvQVRFvp8nOOrnMjKKQUGyJMYb03lHMhZx7KIpDpDTF5R\nyrLRi6Pc3r7YNeVuDoh3HZ2IzabO6rJsNCHaoLHvAObSrtIaKiry6V+xWAytra1BAB+X9eJ8dM/S\nj4hSgbohyOxpQ+tT8wHMx4DrrkT1uZMghSvAKgLoWLsNqVT+xtO1AKrOPyN385pwWqZVjBmB0Y/d\nJR7fJa3D63lz2WRZLtliqHJtOUdqZ7TP+ZuYBNak03xdI7q55V9XxmCZc3Wy+kTn+GtEbUudn5OV\n51WYkR/PkawE15ouAN6SshRrTasWotq8eTOqq6v3tLa2lucf8NFtFhWC9UcgtW0PpBojOr313/PQ\nOvslAEDo+HEAgNjmfGCraWUYDvh80B8gJh63DH6hSd9FycXuOWmFTvzegK4KYzBRqiWVm4dCuaOz\n/YnCIArCHGzn7eM4EZ09TafUz5FsRLV69WowxtaV1EknwAiQNG/H4YJuu5MSH3zw/eT6LaB0BgsW\nLMDA66dAa4uCVAalrg6KoiCywIgBM0nK9OGUGzDpVT/Kek1pKSpA4S6T6ItdDkl5cTa77baZu5j2\no6uRs4AEpFNKH3bi8vpZ2EnKhJRhBfr7BdcKCcvBqrZ9rm6wCh2qCIfzxWZefPHFWGtr6wuuHfhw\nRbcRFRG9nly/FcEjR6G1tRW/Ou1C6MkUmEJQBg+EcvQYxFZvBhHlpILNbHeDrLgIcJcvEg/rTSyO\nJi8I7HSwNOy7ZaK5OBGWW/KsdYu8MHUkd9OWaGW4EZIbaZVLZsLlXtZy64wF50ZWos0GEyZB6SX6\n19wSpO3hKcU2UMzVgJ5JIRQydq51XcfChQtlAK97mpgPIbpzbbKF0mmk1m/Bvn37cOKJJwKq8asT\nHDUMWqQDIODgrAXCi5lSnlUgJyRHRYGuQjHLTURmhUJx3qKvxTLCVjIo9v64EXVXQURMTueKEVo5\nqUCAd5LyMl7+OefPR7SDS7qOTONeTJgwAQCwbt06MMYOEtHOsibnA0D3WlSkDDZUFP7rp7dhyJAh\n0DriAICKEaOhtrTi2ks/h7b/vIWOBkNf3Us0MB8QKkq70QOE1O7daP73LCCRz91z7K9MQiu2VV5s\nm9xsY58/UH6ktv094R/z/3f1stDNeirVynLyJxVDOSRVsATk/VkedgaBwsIheiQGFgqhf38jTvD1\n118nXdcLcuB8lIZu9famd+66FQCuveJKjBs3DnosgeTWRqACCB13NK666irU/eh6NP17oZCknEhJ\nlGtnRrKTQtASHehY3oCOdxtc58dLjgjPl+k0NvvwsqwQjQs4La2cfVBOBFUwhmg3jtu0cOq/lOVc\nuSToShyC51znUKL/rGAnsgziI1UCQQHUfMT63Llzox0dHS+W3Fln0Acj07t7W2o2APzud79DKBTC\nU/98GM0PPgmtPYrA0CG45ZnHseO306G2RKAejBjBoIlAbu3PJwCb63+eoMSCdoTgcWMBAO2vLALS\n9qhvcXpKgQ65ykC6jsTmLdBTqZJetN2/Va6ygfV15QnJqzSL4/xsfRQjllJIx5LS5OE6+4+C6RAX\nOca9Wls8QRUjK7f+7DmgTqk6fJEQqSoM0nREo1Ekk0k0NDSEALzhOgkfRdGtREVEBwDgqKOOAgBM\nnToVVWecjJZHn4EybAgyew8gFAqhctLxiLxp5GuKvghuu4FW4jKsKlYlo/7XP4R2sBXtLxk+TMc4\nKodf69i6ddj1m9+g6ZEnkNq5yzlmx+Muk+MWegmkJSp0USoKib/4AcBiWZWCzi4vi+3ilQKvVhYp\nhQnggPg7ZK4ETCe7pMpQBg3Etm3bsHTpUlRVVW0loj6hk34o0S0Bn25oev4V9Dt6HDoWLAVk40v4\nxh//gTM+cyGqzzsDUpWxrWt3VJq+APOL4RRYad5YweH1GPTNqQgMEooOGm1tMTcm9HQara++Bq2t\nHbUXXojQ8DFGG1uQYLFUE7ddKcdf8U5aS8I+VbHCgPcwkMIdyq7erCDFqKCtB8hCTl59T+b19uBR\nlnFOQxLtFJJC0Cvd0pzMwF+RfJBhBSoDB2L79u1455131Hg8/h9PL8CHK3qCqCwJysFgEAfe+8Bw\nNmYTNydNmoTQ8UehY+kq1H7ubCFJif7nl4gm8lHihKqzjPJbplwLgKLJuOntu9F4318AAIOuvQY1\nZ59tPOFSx7SUm9bNQtAD1pvNOoZkeR125Bz0DjclT1alJFDn32PKFXQArESQG6OEjID8NVzkPEdW\npYCP6+Ln5jRPHpIKAMaYetgbcbOMBB223dRsKhQFCFIwgGg0irlz58ZTqdSrJb2YLoBZLqsvoSdC\np+8AgFWr8goXtbW1aGxsxPJ3382de+r2XyOxVpxc7rYbWEyyo1ToySTk6moEhtQB8XTB824O3qJ9\nZ3clnW5EEYl1ZUS5CbslUEqJeZE/y42ovQSfUsC7T6ngWq9Syrb3UcowSCqgK8bnQrYlrqc+bYnO\n5mtU29sRDoexY8eOAIBlnjv04YhuJyoiilZWVi6eP3++5ZsybNgwnH766bnHkydPRnr3PiS3GOEm\n9iDQrkKxX/jQ+KNQ/7VvYNAVV6H2nPM8RToXg6gPEWnpNkuAj5MSWisO/iThPHPOYOvrL+X9Ffms\nSkWpZGWPYOf9TK5ifE7vmZp30ttJyq0fXhtNlPCclwLSkNl/AMuXL9eDweDTRFT4a9cHwBirYIzd\nmq0iM9OhjcIYa7YVhghmn7uVMbaCMbaWMfbjYuP1SDJaPB7/5YwZM2JuKqxVVVV4/ol/48Cf/4XY\nSkMNIxhQEQyoBZKy/C4LIHbAF0StWyLd80sNUf5Z8KjRqDrueDDF28rYycLy+mvPE5bhyC0kKDsJ\nFT4WSyTbE7ztVYx5eV179LyaCAitVFdCFMzTyTFvh52sLIRlD3J1WuJ6yDQAjPdcC+u595svKuGc\nqM7vAhYueZlCSG3bBnnQADz55JOJaDT6gLCjvgEVwEYAvwfg9Gs9EsBrRHQhd6QZY+cA+DKAcwGc\nDuDzjDHHAhNADxEVgCUtLS3RlStXuja66qqrsHLxEhz45wuILP0gd970qZg3VUHqAheQl97RiKYH\nnoKaFdcrF6X6ScztdCcHPQ8nK8FOUoA4aLMzZbpYRgISMvSkXJDaU2wZzRNhYbCt7mmexSAKDrWT\nT2eWw3brttjmiBssGyfZ1xpf+T5uuOQyxOPxZgDuX/jeAddyWU4gIo2IXgOQcGk2BkAdY2w+Y2wJ\nY8yUKL4CwKNElM5anI8AuNptvB7Z9SMiqqioePAf//jHHaeffnrIre3EiRNR/9MbsefOB9B/4g+g\n1FQCANREIHeTuN+UDPEV7yG+4j0Mvf07CI06CmZJrXLhtPwTkZmUYS4u79Ig0movVi3a+tj6qy8W\nivNIcI5FWPObF2a7cn2FhZaWu3xNMRToj6n5XT7+R8Frn0U3BjIq4ivfx97Kwcl4PP73Q1bIgQAp\n43lo13JZjLHPALhT8NQ0ItonOM8jDqN+3+9hFHNYxBh7D0ZNwHe5dnsBnOHWUY/pkKTT6Udmz57N\n4vF40baNt/8NNZeeg51/mZcTJiOV5awAHmYxCPNLVDGuHoO/NxUAsP/3D6GjocGS1mKimASICUkV\nHwAct9HN8/agRdGveckpJh4q4DhqiLs4tQuDa62H41gBHQjnJXrsfdnnBYhveNMis1hjYa3kYFPz\ncJJ2trxmlxg2sW58/ntiWS5mSTrx3kYodXVYuHAhU1X1cc8T78UgokW2pZt5FCMpENFyIvpN1vpq\nB7AQwKkwqtDYawIWKmBy6DGiIqLGioqKFc8995yn9s1zFiC1swlNzy0FkF1OZJctxUqOV5/9KQy5\n9SsAgJYZs9E279WCiPRi/iOTVJhqHHLCOMzHkku4An+9ndyMPjsfg2S1srj3IyFbb9IigZx8X8Wk\nkgvnQBBZfXa45zfqtmVjIck5pfZ42VEU+absn71zjURb/iZHUvzcTEQXLcdXzr8QoVBoMRH1lfp9\nZYMxdo653GOMVcCoirwGhlT51xhjgWwp9xsAuMrg9KiyW2tr67333XdftHhLIBAIYMYf/4qDT72G\nPb97MnfeLuFinLPKcQBA1eknYNgffgB5UH+0v7wARIVBit6ilI2/elbI38yRMslKtCwUkVMp+FrN\nggAAIABJREFU45a7o+YYfe/BWuIJqtRCFaXoh7nlKZarQ9ZdEKU7OYUxpBv3I717H1asWBFta2u7\nr8cn20tgq+u3AcAUxthKGEvA6UT0ARE1wCCmFTDCN+Zlzzn325PLaMZYMBwON7///vv9xo0b5+ma\nNWvWYOLEiRjxt59DDtRanstvGeeTavl8QFIlEBEQBVhAyX/ZbA5UPs/MCSbpMBv5kGJsc7u1zZGd\n4uwfEb4uT9ZKIUk7XVtOoGexcfPjiVRY3ecmgt0P6ap8IdxUKJRoKV52rNCiAgotX3vUOv8eN8+Y\nC2gS9IZ1bfF4fAgRFbG5uw/Vg0bSCZf90FPb5f/+71VuPqregh61qIgoLcvyEzNmzPD8IZ5yyin4\n4Q9/iObps6HDmhwsqkDCaz8BADQJLFC4Z5D/VfeW/6VnCYkE2w8i/5Wd0OxkVgpKdU6XS1L2svel\noFBvy5rE67QTaLfsSJWyvkhb9aASYsZycyo7P5GEy3PLDxw3vhaNIb7iPUw7ZkKGiB4+lCTVV9Hj\not4dHR333X///Zlo1NMKEABw99134+oJE9H093+BMs7fAXNZaI8YdgNPVvZdPD4o03zOTlamz8p+\nFENX5crZwwV6YhlVjDi9lIi39+VmfXkNc3DVsy+luo5tM0TK5gs6bUJEXnkH4U+Nx8yZM7VEIvGg\np0F8lIQeJyoi2sgYe+3ee+/1vA5RFAVPPvkkpEAQzY/Pgi6rFh8H/8spcoJ6nptSSEwi8GRVTNtH\nZIGJYNevAmzhBR5FBY2/pZMUX96rs+lIpcZReVnm5fsuvuvZFZAyzPG7wP8Y6vEkoguX4/qTTldl\nWX6eiLZ1++Q+gehRH1VuUMaOra6uXrt79+5QbW1t8QuySCaTqD3uaFR/+kxUn3taYb+2mm1ucJLA\nFc7XdgPxO3pMtZKUzlfP5UhKzznlvaR/WH+93QjIKTCTX+p5Wco59ePUzim2y37eLVnavhNYkEMo\n2JkUzcUTwRXxV/HJy8Jln8ivGNbQ/vIipLftAdZvS8bj8eOJaHvRyXQzfB9VF4GINkmSNO/Pf/5z\nSWv5UCiEt5+Zg/bnX7UsAc1fb6+7ZU6a6k4R0XZSEQd62voqYkkV06my9GVui5dg5ZjkJLKURId4\nXG/lv/j4KSc4Ve1x+nGxkxRfabrUXUkvcIqn4j97+y6gnkgiMn8ppp54WkaW5dm9gaQA5AI+vRyH\nCw5Z4blIJHLHn/70J7W1tbWk60477TQow4YivsZIseF3+ZglBcU5iLJYeSZrQnDxJaAecLekvDrS\nXa0sAVkVIy57tWanKjlOY5nXmoe57FLCmZJJgylU6CC3ScbwbfP/GyTllqTuJSHb+OvB8jLfZ1vy\ns/0HRa/UEHltKULHjsWsWbO0aDT6i6Kd+ygbh4yoiGirLMtz7r777pJ3SGb+9m60z1sIPZUu+DW3\nS+HmnOUCZQK7AJyrXrgtSdYqtpYnK9Fyj4dTCISnmC6OrIQR+rnDeffNqU8ePKnYn3cjDadKO3y/\nuf8FgZ98G5MUzfHMBPVggLekrTuGnZZo5khKdN6ElGHQOuKILngH1x4/MaMoypNEtKtTg/twxSHx\nUeUGZ2xMZWXlhh07doQGZyvWeEX12adCqq7EwOuuKPSXJLKWhE3z2vjLHK0k+83itjyxaHyXGIog\nWjralQLs8+8MLOW7nCK47dHaBaEBedLgYSaJC4mJkz8xya/oXG3+ODmgoTJsVUqJJ4JCv5vTHIy/\n7snITiQFFBKVHtbR/NTTIJ1Aqz9MJhKJcUS0p+iL6yFUDxxJJ138A09t3336x76PqhiI6GNZlmfd\nddddJesRbnx6LmJLV0HrKJT7ETllnWJqvFRzMa0yy3VKPqTBXN55XeY5FS7I/V+keo1oR9OtZJcT\ncu2KyLa4kRRgc5oL5mC06ZpwjHRGgeLis/IWX1WYQuNk0YqCPpM7tyP5/iZMO/aklCzLj/Qmkuqr\nOKREBQDRaPTOBx54QNu/3zUnsQAjRoxAcORwpD74CHpStiQs88TExxfx4mcmTDKSQprjlzzXhy0h\nlY+/cjp42NNq7NVW7HXlih1SXIYUl4VkxFsSpVa8sb43YpIyYVhThXl8nY0T45e26YySO8zlH7/8\nLE2h1JqnZ3EVFEttUghaWEXbsy+j38XnYubMmXpHR8f/lvCyega+M73rQUS7FEV5/Pbbb0+Weu3v\nb74F8XXrc7/gokhmHvmM/MLoaP6vcJ7ZmoHlwkJQglxAV9IS7Ax6JbX8/B1IWCj9YoYUOJMUv+Sz\n9+NYL1EQWe4WbW7qZakC0rLPpRg6U2Y+Zz1Xakiu+wB6LI5LpH5JSZL+QUSl/cL6KAuHnKgAIBqN\n/s/s2bPjb775ZknXTZs2DYl1H0BrjgpuCvvyj3O8cvIh5k3pVpLLhGOJcsVZLreYAF8x0iq2TOTF\n+kqSi7GVnRf97wS7X6pAmcIiXcMLAIqd38XkZOxkZc6BF1IUxVSVGvBbOPd8XiYFCJRO4OCTL+B3\nP/wJXnzxxUgsFhPpNPnoBvR4uSwRiKhNluUbr7vuupmx0eMqlf4DcODF58GY+01TX1+Pr0/7Mv6z\ncg1qJl9gCVUAuJ3A7Jc458sIaMaXWyEg+8XmnbNOQYuAe7hCwetSxDljTjCLDQifK3Mp5TbfPLlL\nQNZLaJ4xCcbutLaTAr/RIPL/AdYlZK4fpzm5/FA4lUorN4peTLC2nUwu7KT92UWoOHYc7rnnnngs\nFrueiDrKGthHyegVFhUAaJr2QjQaXdD2zltoefVFTJkyBYmEm8qpgWuuuQbxDzcYZaYcSIH/8ps+\nDdN64pcbopggvk1XlDsvBtES0UtbrxA74h2WizYtdV5D3f5eO/l+OhsyYIKfAz8/a1xZ6daUnWSd\nkp/V5iZ0vN2Aa489MZ1IJOZlZXh7JRgR5Iy343BBryEqAIhEIt8CgP4XfRYLNn2E+jPOx80334wP\nP/zQ8ZqxY8dCbTkIJWE4l+2/rm6BiXJAswQv2tubhxLO5JztIq1wO2HZixSIdgbddgd5gnJrZ8rG\nlAsngrUmd0sFigaFqS9UQN4idDZank9gdpNd9gp7uozj3FNAy4xnUTXpFDz77LPxaDT6vZIH89Ep\n9CqiIqIDjLE/tS58DUMu+zzSzfvw0EMP4YQTTnC8pqqqCpRKQYppCLZKYBGl4Mvutjsk2j2ynzMt\nMCWcyVlY9sPNynJKdi4lrMGpvBbvDxNZRQW+O8EOoDNhZWPSsqqh9tADL+RkD1I1rbLcktrsP6Lk\n1FtLgSgo1Ys15ab0apWdYYi8vhjQGfp/tDOWTCa/SUQHS5pkHwVj7LeMsXcYYysZYwWR+VkFz4cZ\nY29nC0fcxD232Hacbr+eR68iKgAgop+Eho1A64olGDE197rQ/5zzhO1HjRoFpbIfUjt3IdgOjqyM\nL7CTI9YOM+LZjIJOrv4QUtPugu3w/JLRfhSSFpBfUuhmaSYXhQY7YdkTmXPnS7SinEIUnMjVkqAr\n3EG0VsgplsIiSt0x/VtSQoISkSEnJIslZ1znlocoXuoXIymLSoUq9q3xSG39GJGFS3DjxZdk2tvb\n39I0bY7rAJ8QMMYuB1BPRGcDOBPA5Yyxk2zNvg1gJxGdC+B8AHcwxgZkn6uwabCvcBuvNxIVJXft\nGBVdvRJpimPktQZZtb/zNhYtWiS85vb/dzPaP2iAHjCSgwNRKRdfpCdlIVmJIpv59IzGv8/Dhpv/\nifY5i6FG4jnCqgynURE2PM/25aFJWuYy0UoAOvRKzTjCeo64tLAu1LwSqS0UlzAuLBbhpfyUca6I\nZcRFyzsRk1d/FMsYn4+UMD6TXI09fo4cEdlTgRx9keXI+lRqwiISpErQIh1ofugp/PW3v8Pjjz+e\nikQi3yhpgMMD5ZbLegkAvwSWANhDjB6EUYEGMPZpVAAaY0wB0J8x9nS2gOn/ZbXTHdHriAowYqsG\nTr4M++fMQsWnjsPR3/0F5H41uOSaayEKYbj++usR3fI+NFnPJQfzv9B6UkYmGkQqEbAQFl+A0x6f\nc+ZjNwAAdj62BLvveRZEhIBi+LGCAdXi8zKXjDxxmVYWX7zAXl1Fr9QK8gjdAkYBa0Q8IIicdil8\n4MXpLyIsO0k5XutRboVfduVIO0CuZFms73JCEfRKDVLIYf+RZdD84Ez0O3EiHnjggVg6nb7VS+WV\nXgECpLTu6UC2XBZ3TOe7Yox9RrBMW8wYqyeiFGNsOAz98+lEtNkyDSKdiDKMsfEAXgTw30QUgVE6\nazEMi+tCAEcAuAku6JVEBQBNLzwrSxUVaF40HxgxAMd895cYcslVuOhzl2PZsmWWtmPHjkWgZiBi\nH76Xk1uRVCtZsYwEigagJgKWYpvmbpKdxPSqGpz/wvfQ/+SRiK7bgW23PYLouo9zZMX7u0RLSZ6s\n3GKI+KTpYoVJ7X4vS4Z/QcyS1erpqrxBHvZKOJ6uyTqw+TirYgTVVek3gNXqFJGU+ZpaZ70KkIzr\nTjw+tWfPngZN0x7rskkcRnArl8UYuxDAQwB+REQPi65njE0F8GsA1xHRvGyfbUT03exfHcAcGBWT\nHdFriYqI9MTWLXUda1cj+t4aUJBhwJhPYejnp+G8Sy7Fnj359CrGGJ5/9CHsX/Iy9Ew6JwVskpWU\nIyyWddiKdwbtS0SlugKn3DMF/U8egfimPdhy+1OIbTbGJd2j9pUHRUqRjyi3TMwe/A1tf85cvhT1\nFXkoL1UMbqqjxXISWUYySIojacC7r8sJpVhTudceLswRNIkr8to7iDesxx+/fys99thjLZFI5AuH\nrJhoLwVj7DgAPwIwhYg2ObS5EMClAKby+ZCMsXrG2B0sHyh5KYDVbuP1ioBPJxBRE2PslKZ5c9YE\nb6xD5YDhqDlyAjKT9mPMiSej4Y0FSKfTaGpqgqqqyBxsRry1EdWBsdABkGJu81urF+sBVnBDmBaS\nlpGRVhQks5ZTSMnguN98CatumI5aFsLWn89Cv1OPQqRhK/qddzIGfPWqnFYTD3tIhEgbCuAsK1t9\n5WJqBp7evxItEcdyW3xgKIy5Wqw11QgNsWz1c5aemx6UXS2hlODNUgnKOq74cXLjdrTPfROzHn4E\nN9xwQzwej3+WiEoTTftk4CYA4wC8xgVm/xnATgD/Q0TTAHwfwFgAb3BtboNR6r4awGrGWAeAtQAs\nS047ejVRAQARrT1i2tew59+PYMyNP4RSVY1BZ38GcmUVJp59LrRYvkjEgIs+i+DY0dB05Kwqp1gj\nUvNkZe70AUZmvqpKyKh5354UBMZ++9PYt2A9tj+3GvPnz8fPJ7yK1mcWYfCNl0ND/gbjI7D55aHZ\nr9nWLg7HFPeb1C49Y15X+Lq8a6vz/VkkWTwQFk+sJkmZ0fPlxnZ5mnsZvig7TMuJj7IHAC2WQNOD\nz6D22s/h5ptvTiSTyeuIyDmI7xMMIvoxgB87PD0t2+Zaly7uyB6e0OuJCgD2zvoXG3jBRdT4n8cx\n4vqbwWQZNaedgX4TJyGyYhm0RBytDUsQHFoPJknQJeuaVhwrU3jO3PULBmBxrIcUFdENe5BpjmL4\n8OG48cYb8b9j/xfQrfFXlmVjjvz4/oG0ouTaahnZGojqYoGJUOzGFmlAFeqvZ3fTbIKDBnEVKqHm\nkqL5lCM1T1AiXS3LnAWWojHP4nIw5ZKUSH/eGEvKjS9JaTQ//DTCpxyLUas2xTbF43/RNG1uWQMe\nYjAdpqO8z+CwICoAaH1rkVJ5zHHqgdfmou7KKcZJRUbNeedATgCDz7sYesDIITN2zKz5cV5SXExH\nuR1JVUHVaceg7f3dCI8cgvVvrcDe1hYM+8lUBAMqAoqWtcCMt9NOUGa/GVXmnlMc5UlUVcr5zERL\nIidVzq6ESVYiyAlu08EWhe8E4bIva8V1VbiBV9jzN+VQGs3/+A+0VAbX1B+dnLt87pJEIvHLbhnc\nR1k4bIiKiDTG2IDA4CGtbauXof/EMwFkl3b9ADkJHHhpDuLbtyA04ThUjh+P8NFHW9QF+Cq3QN6q\nMCwh460wSYUnrdTeVmz92UwMvfZMJPa249Sbp0LqF0ZgcE1B+4CiIaQYZJRUrW+vExHaZUt4AhPF\ne9mVM/nXknu/XAIleQvGCSK/Em85meWkSiUn0TjFlpydBT+GPf2H9cug/dkFSO9uwgiEac6cOc3x\nePxL2d0oH70Ehw1RAca2JmNsfMurL28IDq1H5fAxAPIWVGj4SLSteBtaQxwdDatwxC3fQ2hgvWFZ\nZcwt8cJYIFF14EDOma4isa8JLCBj/7P5sAgWUKAMrkVGZQVkFVIyub9JNVBAWHxb08qyk5UJ3rLK\nLY9cVAnM18O3B8Q3K/+c9XFxXXdNtJx2kfN1Q2cIykm5gUeBI90cL6whunA5Im9/gMDw0di6ag2D\nsUPlvTqujx5Brw1PcAIRbUQ8fsW+Jx9FvLXR8lz4vEkY9rVvgYjw9S99Cc0zZ0GTucDMSJ6MCkq/\nA0glAogngkhnFIszHQDOP/tc6LqOAZNPxoCrz8Hw334HUtCILjXb2q/hkVHlguf5x6LiBXwgKWBV\ncigFdmniXGQ4pw6aUw1NGCEEcvavaxkvpxxBB7G/roBTWo8XFdNcPF02ROHgE8+j5dG5kIfUI75q\nDRhj04jonS6b7KECEaSM7uk4XHDYERUAaJr20sCrr8a+6dORbGq0BD5WnHwchnx5Kh6bNROpnTux\nb/p0ix9FzuYBAoWWhUhhIakqSOxowbBhw8AYQ/qjXRh81hhUjK7PRbQDVtJJqgFLH/xzJmGJSC2g\naAgomoWw7FHvfIqO8L1xcMIXqH/a1EGBvEVkFqwoVwPLDuFmRgl9d0ZOuTDK3ngcW7EC0YXLAQDJ\n99cDwIO6rs8ueQAfPYLDkqgA4MC/nmADpxhkldm51wiSzObPhU8aj7pvGWlZgZH1lpQTADlLgpcy\nscO0qiKrt2PXrGX42c9+hk2bNiG2vQXxd95HbT9DK8uefpNRZUSTFWjqqEZTRzVaolUWzW/7YVpw\nbuB9VjyZkspy6UHmYZdkcaugY/7vZvVINhKz9OOBPIr5sIqRUGe11y1zyZKU2rgFLdOfAwtXIKsU\n+FMi8qVbejEOKx+VHQcefYLV4XraN306hn73WwiMOQJmrb7QuLEY87c/GkuTjHHDaE6KmwkZFIYl\nPMDEzgdfxexHnsCECRNw5513ov78IzH+xlOhhxLIqDLiiSAXllCo5c2HIDiVmhI5zE3ysucm2lHg\n38kAolu71OVXKSXoTdgVTUv3VXWelJx8VuaST921CXt+/QikflXQO2IA4X+I6O5OD+yjW3FYExWQ\nJ6v9D/4Tdf/vJgRHDCuQKeHBazeZu4AU0AGVQU0EoIQzUFUJiqIjlWRI723DlVdeCQB4tmIDEhEV\nem1/JNVAziEuIhAAWe0lJoxcN2E6y+25hiaclnIirXIRihGUqNimPWCTjzLPtXMZszNFMLoavOM8\ntWETmh6YjdApxyK5drNJUn84tDP04QWHPVEBebI6cN/DObICkLOuAFhyzEzwO2GAcfOriYAlGDJ8\n4pGoHjsMNROGI1gVAKm6xQfFp944EZZBVrKFsEwriycilQtAtF9vzjc3dziThb34gwkv1lGOyMt0\ngPcWkrJYVmENHYveROvcJQiOGI7kui1gwDd0okcP8TS7BwSwtJMq/eGJPkFUgJWszGWgHcXIKp+m\nIuWKDxxx23XQ9uwDbdiAxI5mHHnrZ3PhBnaHOE9EfI3BXOVm6DDdghp3jfGXS2Xh0lNYRsqRUjnk\nIbKOnOSHRXBq35MoJ+wh91mHNaTWv4+Wf81H1TmnIfZuAxjhel3Xn+yu+froevQZogKsy8Ah37wB\nFceMzj1n+k+8WFZAnnSYokAZNQKVx9ZhcDaYM6k6hyLoSblwRy6sFewwFlRyMc8LVA1MkhJZRyZK\nzbOzB1narU8niBzzXglEVHW4WJtizzn1QQpBC2UQnTsfHUvWIXTCMYi90wAGfM8nqcMPh+2unxMO\nPPoEo47YFQf++Sg63lppeY5f1phb8yZy+uB8+6yzGzCc27FEBaLJioL4J0XJB1LaC6Hm+jdvqIS1\nSIJ9V05YEDTgLmEM5Amq3JCCUkiq3PgoSx6gwDfWmf7s0AIptMx4Con3twGyjOSGj8CAqbquP1jS\nID56BfqURWVC07SXGGPj2xe+sSHd2IiBX7gKTM4uv7idKd7CyjngI0rOyW6mmlhX+wpUVUJlOF32\n/KS4LEz4NSG08BwCK3mYqS1OZKVbnOFWR7xXgrLD024g91pFtQ7LWdKKlqSm5akihgN/ewwEILO7\nyfjsNe1Y3aZA2VfBiCCly6il1ovR5ywqE0S0UT3QNEBtbsH++6dD64jln7PpjzOV5cT1gHycVd76\nkXLKoGZIQTQStsRPGcGY5g1pK1oquBGlrKCfE0T19wpeo0AJVFQIwjwKl7zFLaPORJc7abUXI7fO\nBJpm9A7s/9tDoEQK6U0fQxk6BHosXmOXyfVxeKHPEhVg5AYm1m9UKsaMxt4//hWp3daUG9HyI/d/\nLhAxT1h6UraUe0pGK3JhBbndO873Yy1QYFVyKCBLQVkrz6+TL+KgGIUSzEMkZ2yfR2dQLFXGribq\nJP3Cz4MvA+alqIWJZGsj9v3hXmjxJNK79qL602chs2O37OfudT0YY2HG2GOMsSWMsfcYY/8laBNk\njM3IltRazRibzD13K2NsBWNsLWPMSdcqhz659ONBRBoANuTrX6X990/HoC9+AVUTT84/X2Q7nmWT\nmZ0lSaSC+CdTIsXNIuKd1/xYoscllZHn+i13t67cEAUnLapyqsOU+ny6pREH/jIDelsUkBgY4bro\ngrefKmlgH6VgEoCFRPR1xlgYwHrG2FNE1MS1+QmANiI6O1sEYjFj7ITstV8GcG623SLG2GIianAa\nrE9bVDyaHnuS6bHYxIMvvITWeS8XaJ67/XLnfDQC8rHXqisFXomk3Dw3T/2KHP9d6Nj2qsxQbPkp\n8r2ZfUWWvoO9P/+rQVKyBOh0kq7rPkl1I4hoCRE9kX1YD2APgDZbsytgFH8AETUCeBcGOV0B4FEi\nShNRGsAjAK52G6/PW1Q8iGgNY6wutWPngX1/vR+DvvIlBOuHerrWtKyM/61yK4Vte5b/3fLknMiQ\nv6bTEsIuScdeotlLJUV+AyT58Ta0Pvk8AEAZOgjq/pZBn/hKxkRg3p3pgxljvCUznS+ZxRj7DIA7\nBddNA5AG8AqMclffIqKMrc0gAHyJsb0A6rLn37WdP8Ntkp8oogJgFoyQB35xirbv3gdQc+F5qJ38\n6dyuYEF7l5uXJ6TuKEfVFfCaXlMuSXVmfPscAGe/oWgzQG2PYP8fjGiDfpecj+irb0l+tZiS0UxE\nk5yeJKJFAMSVfw2czhgbAWABY+wKItrKPbcfBjFFso/rs+fM87Cdd8QnZunHg4j0lqefY3osNjr5\n8cfY88e/IrlnV6cisEXloYpf0/XLOc9jdxFJddYZ7xRbZQkdEezSJnd9jAN/fQiBkcMA4JjIK28y\nn6R6DoyxLzPGTCtoD4AWGJVleMxFtrAoY2wojNLvS7Pnv8YYC2QrJN8Ao4ipIz5xFhUPItrJGJMG\nfX2afuD+Gag+axJqL/8spGCgR0hE5FDvTnSlgF1Xgs8acANTGTQtidb5ryDesBZ6tGMqET3jE9Qh\nwTIA9zPGBsIweF4gorWMsZORL5f1NwAzGGPLYaSn3kJEKQANjLEXAKyAUeZ9lpsjHQCY/xkbYIzV\nVZ76qf3pXY0Y9NUvInT0kT0zbhcTVU+TUVckIRebs+lD69i5CS2znkHF2DGINaweTEQtnR68D6I2\nfASdPfZGT21f2fD7VW5Lv96CT7RFxYOIDgBgsix/vvmRp54Pnzge/a+8BHI/uzXbu9BbraSugpRh\nSCfa0frKfCQ2bQa1Ry7vaF718qGeV68GEZC2+7UPb/hEZYOmaf9hjA1gsty659f3oPrcM1Az+QLI\n1VVdPlZnranDkaRKmbPW0YGDi95AdMVyVJ92GrTWtloiihS/0kdfg09UAhBRGwDGGBulx+I79vzq\nblR/+izUfvoCSJXhwvZCxYNi5agOP5IRwYsiQqnQYjG0L34T0XffRdUpJ0OPJ0a0L36zsfiVPvoq\nfKJyARHthEFYY7W29m2Nv/wD+l10NmouOB9SOOR+ra2Ap1tMU19AVxCUmo4jsvgtRJcsRdWJJ0GP\nx0dH3l66swum5+Mwh09UHkBE22EQ1tHqvpbNjXfehZrJ56PfmWdBj8WhDBuUb+tSTLO7rSh75LaX\nsINSNaG6AvYxtWgHIsveRfSNt1E54Xjosfi4yLvvbuvWSfRl6L6P6hMNItoCg7DGp3fvXd/4f3dB\nj8UBAENu/rqh165pUAb0BwL2/L/iEeKdmlu2eEVJ1zjMyas2ernyLCY6tryPtudegbpnPyonngS9\nI3ZMdNmKLSV36qPPwyeqMkBEG2AQ1vEV48Z+mPpoO5ofnwUAkMJh6IkEqiadjPCJx6Ni7ChXR3wp\nssAF13okJq/EU+zazoICBE1PomPxUrTPewOUSiM4eiQAfCm2at0zXTaQjz4Hn6g6ASJaD4OwKiiR\nvGrAgAH/lU6lT736C1NYbW1txaaNm7DyiWcQjccBTUN4wnhUjB2DwJDBCB17NOTKynxfLkJ6wrFL\naJ9XgHBXZyivwKf1WhGxqYkkmu+fgdRHH4OFKlBx5BhQMnUXgF+ntu9Iljyoj08c/IDPLgZjbGRF\nRcU3AoHA9+rq6ipvueWWqsmTJ7NIJILGxkY0NDTg/pdeQvKjj6AMHIjA4MFAUIFUFUKwbigCY4Yj\nNHIMmNQ92U1eCiV4rc2nJ5JgigIWyP/emYSodcRw8D8vIrbCIgfdBOAsIvqo3Pn7KI7a4FA6u/7L\nntq+suvewyLg0yeqbgJjTAJwQW1t7feTyeTlJ5xwQurqq6+uvvjii6VJkyaBMYY1a9YMU5TmAAAE\nl0lEQVRgx44dSKVSaGlpwc9mz0Ry+8eQKipQd9MNkPv1cx/DRiiddYIndm2HFA6DqQBlMpBrayFX\nVoIqJDCWTQ5OZ9CxbAUS6zcgsX4jAECuqYEWiSB8/HiEjjwSrS++ZO/6QRjpE/6XrQdQG6ijswdd\n66ntK/sf9InKhwHGWDWA8yorKy8PBAJXpNPp+gsuuCB95ZVX9rv44otx1FFH5YlA1yFnlRxG/+UP\njqoOlv6LEJSWSCC+eQMYY1Db2qF3JJBpOgCt1ZAPYhUVoEzaUEBV8/IgTFFA2cf9r7gMTJLQ+sJL\nYMEgxo0aha1bt6KqqgqxWF7mWZKkNIBtuq5vAvA1P0Cz5+ETlY8uAWOsHsDk2traqzKZzMVVVVWB\nyy67TPrc5z4XnjhxIogIJ13yWaitbag8aQJCRx+FwKh6KEMG5ZaETnFZ2oF27J/+CCpGjoBcWwO9\nI4boO8sAWUb4uGMRGDwQP7rwIhxzzDEYNWoUACAWi6GqqgrRaBREhEQigeHDh6Ourg5r1qzBW2+9\nBVmW0dTUhEgkAiLCsmXL9FQqlVZVNa4oymvt7e3zACwion0Fk/LRo/CJykeXgxmm1HjG2OQBAwZ8\nPpPJnJBIJAYMHz48PnLkSDkYDFYu2bSBZRr3AgAG3/RVBOqHIlA32OIbMkHRFHbe9nPLuWHDhmHz\n5s2oqvKWBtTR0YGNGzdi48aNWL9+vf7cc89ltm7dGtR1nUmSpCqK0phOp+cBeArAMn9J17vgE5WP\nHkFWg/oYAMcpijKhpqZmUjqdPjkejw+VZZlkWWbpdFoKhUKora3F4MGD0b9/fyxZsgQAMGXKFMyZ\nM8fS57Jly5BOpxGJRBCNRi1/29raMq2trZmtW7dmNm3aFIhEIoHq6urdkiStj0Qiq1RVbQWwAcCb\nROTv0vVy+ETl45AiKzI2GsB4ACMBDGOMDVcUZbgkSYNSqdQkxlg6FAod1HW9UlXVsKZpAUmS1Jqa\nmkZZliOMsSgRtWua1ppOpw8mEomDRNQOIAqgEQYh7SCi3ilZ6qMoapUhdFY/VwnyHF5tm3FYEJUf\nR3UYIVtRZ1v28OHjkCFr9T8IYByAWgAziOheW5vfAziLO3UagHOyAnuLbV3eRkQrnMbzicqHDx/l\noGi5LCK63fyfMXYqgDuIaG32VAURnQWP8InKhw8fJYOIlgBYkn3oVC6Lx90Avg0AjDEFQH/G2NPZ\na98E8KvsikEIn6h8+Ohr0HVQMuW1dXeWyzL7uAjAbi4joRrAYgC3w6hQMx1GEYiHnCbpE5UPH59s\ndGe5LBM/AfALrs82AN81HzPG5gC4Bi5E9Yksl+XDh4/OwWO5LDDGxgIYQUQruXP1jLE7mJmOAVwK\nYLXbeD5R+fDhoxwsA/BLxtiy7P/zzXJZjLFZXLtrYCwPeeyHQWqrGWNLYJTSmg4X+HFUPnz0MTDG\nXgEw2GPzZiK6tDvn0xXwicqHDx+9Hv7Sz4cPH70ePlH58OGj18MnKh8+fPR6+ETlw4ePXg+fqHz4\n8NHr4ROVDx8+ej18ovLhw0evh09UPnz46PXwicqHDx+9Hv8fSSaSvvGHHXcAAAAASUVORK5CYII=\n", "text/plain": [ "<matplotlib.figure.Figure at 0x175182c7358>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import netCDF4 as nc4\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import cartopy.crs as ccrs\n", "from cartopy.util import add_cyclic_point\n", "\n", "from datetime import date\n", "dates = np.array([date(y, m+1, 15) for y in range(1961, 1991) for m in range(12)])\n", "birthday = date(1969, 12, 15)\n", "i = np.where(dates==birthday)\n", "\n", "f = nc4.Dataset('HadISST1_SST_1961-1990.nc')\n", "# 90年生まれより若い人はhadisst1.readtxt()を使う。\n", "lon =f.variables['lon'][:]\n", "lat = f.variables['lat'][:]\n", "sst = f.variables['sst'][i]\n", "\n", "fc = nc4.Dataset('HadISST1_SST_1961-1990_clim.nc')\n", "\n", "sst_clim = fc.variables['sst'][birthday.month-1,]\n", "\n", "dsst, lonc = add_cyclic_point(sst[0,]-sst_clim, lon)\n", "\n", "ax = plt.axes(projection=ccrs.Orthographic(central_longitude=210))\n", "\n", "plt.contourf(lonc, lat, dsst, 60, transform=ccrs.PlateCarree(), vmin = -4, vmax = 4)\n", "plt.colorbar()\n", "\n", "ax.coastlines()\n", "\n", "plt.show()\n", "\n", "f.close()\n", "fc.close()" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "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.1" } }, "nbformat": 4, "nbformat_minor": 2 }