Skip to content

Instantly share code, notes, and snippets.

@angus-g
Created February 8, 2017 20:16
Show Gist options
  • Save angus-g/474b25a669c1df9b13036861cca77aa2 to your computer and use it in GitHub Desktop.
Save angus-g/474b25a669c1df9b13036861cca77aa2 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"%matplotlib inline\n",
"\n",
"import numpy as np\n",
"from netCDF4 import Dataset\n",
"import seaborn as sns\n",
"import matplotlib\n",
"import matplotlib.pyplot as plt\n",
"import matplotlib.cm as cm"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"sns.set_style('whitegrid', {'grid.color': '.9',\n",
" 'legend.frameon': True})"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.colorbar.Colorbar at 0x7ff18f0d2b38>"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhoAAAGACAYAAADrrMueAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3X10lNW9N/zvvCaTVxKSAAkIIbypWEFERUpBg4fncDze\nreA761n4cMQ+2MfevpYWldN6d2kofbk9lceWY3ms9rTa5U3L3eMRQY9HURSNWsHGRrEghrckJCGZ\nmWRenz+Qa//2JtdMEjLDMPP9rMVae+bas699XXMl7Fzz+/3GEY/H4yAiIiJKAeeZngARERFlLy40\niIiIKGW40CAiIqKU4UKDiIiIUoYLDSIiIkoZLjSIiIgoZdxnegKJNDY2nukpEBERDcisWbPO9BQy\nUkYvNABgzj89fqanoPEUllrtrt/9k7ZtzcR/VNvCsbTNiYiIzqzb3n3+TE8hY/GjEyIiIkoZLjSI\niIgoZbjQICIiopTJ+BiNTBOPRtUDp376XA5HmmdDRESU2XhHg4iIiFKGCw0iIiJKGS40iIiIKGW4\n0CAiIqKUYTDoIMUiIdX2lWrbqvPV6Sx0qTXckb6I7XguET8ajQ98HkN9nZ1ofGCDJNqXnJMMjHUZ\nMbLB4ZgwERGdFXhHg4iIiFKGCw0iIiJKGS40iIiIKGW40CAiIqKU4UKDiIiIUoZZJ4MUj6kS5AFP\nibZt5gVVVjuvxGu197x9UOsnMzxkdkZPRP9q+VBs8NkZZkKHXTZJosSPRBkosQT95LHIrBuvU087\naQupLJyusBpxoJkviTChhYgos/COBhEREaUMFxpERESUMlxoEBERUcpwoUFEREQpk9Zg0D179mDV\nqlUYP348AGDKlCl48MEH0zmF0yaDQQ/79dLiNXMnWu0RdTVWu+rCQ1q/WEwFQLo86i0IHO3Q+oW6\ne9VrjCjHWEjNo/d4n2p39Gr9Ir1qjpGgaof7olq/YFTNSQahmgGpIT1eVVPuVevWqppiq10yVg+a\n/etbLVb7sJifGQw7UPJV5nwHGmAqA1kTvd6u32D2K8eQ/VIdyGqWgpfkvoej30ADme3OxUBfY85D\n7neg5fIHY6DzTeV7ac7drtz/QOdgHod8nQzkNvcr+w01kPtM/SxQeqV1oREIBLBo0SKsWbMmnbsl\nIiKiMyStCw2/35/O3REREdEQBYNBrF69Gu3t7QgEArjjjjswY8YM3H///eju7sbo0aOxfv16eL3e\nhOOkNUYjEAigsbERt956K5YtW4adO3emc/dEREQ0QK+88gqmT5+OZ555Bv/yL/+CdevWYd26dViy\nZAmee+451NTUYMuWLUnHSetCY9q0abj99tuxadMm/PCHP8T3vvc9hEKh5C8kIiKitPqHf/gH3Hbb\nbQCAw4cPY9SoUdi1axeuvPJKAEB9fT127NiRdJy0fnRSV1eHuro6AMD48eNRUVGBI0eOYNy4cemc\nBhERUcZ59X//B4qrq5J3HEZerxcXXHBBwj7XXXcd2tra8Mtf/hK33HIL8vPzAQDl5eVoa2tLuo+0\nLjQ2b96Mrq4uLF++HO3t7Whvb8eoUaPSOYVh9eGRHu3x34lMk7yqSqtdOXq0/kKnSzULimFL9HO4\nXPb9hFhQj4OJdbVb7fDx41a7t71L6xdsPy76qTEC7fp4fSLDJW6Ehlecp97L6isutdqe2vO1fnnr\n/qfVHvmOysgJiLEHQ07DzHaQ2TSJMiYGGv1ul60QjJrZOanLOhlotsNAj3EwBpqRcDzSfxZTvtPM\nGFGPw6dkP/Rfqj/RGH6b7CmTnK9dJhFw6i3fgWYn9WrZL4PPThloZgkA5Dvty/3bzVHuN9F5l2MX\nufWzkSg7bSh6RTae/AoDANgXCJ/2+GeD4uoqbLx4SVr3edu7zyft8/vf/x4fffQR7r77brjE/0Xx\neByOBD8/J6V1obFw4ULcd999eOmllxCJRLB27dqkQSRERES5Yqip16mwe/dujBw5EtXV1Tj//PMR\ni8Xg8/kQDAbh8/nQ1taGqqrkd2DSutAoLi7GE088kc5dEhERnTUS3WFLt/fffx8tLS347ne/i7a2\nNvj9ftTX1+Pll1/G1VdfjW3btmH+/PlJx0m60Ojo6MDOnTvR0nKiyFJ1dTXmzJmD8vLy0z8KIiIi\nsmTSHY0bb7wR3/3ud3HzzTcjFAph7dq1OP/883HPPfdg06ZNqK2txeLFi5OOY7vQ8Pv9aGhowKuv\nvoqLL74Yo7+MM2hqakJDQwMWLFiA73znOygsLBy+oyIiIsphmXRHw+v14sc//vEpzz/99NODGsd2\nobFy5UosWbIEa9eu1YI/ACAajWLz5s1YuXIlfvOb3wxqh9nkrb8d0x7/owj6dBaNUBucRiCnKGPu\n8KgYlXhULwuekNujhvfmW21XWaXWzVGnoonznOrtLnLowVaOkAr6jHa0ivZRrV+su1O1Q3q5c+9E\nEfR57jyr2evVA17P/R8VVrv8Nxut9oH/+ljrJ0umx6L25cnjMvAuZH8OZfCqw/izwSGCz2IJxrAT\n8uvBanLuJpco1S5Ly0eHWIJdGkrp8xOPBz++Hgyqn8+usDqHMlDQZwT5yTGCxntsF9hqBiV6RMBi\nR0iWtLc/KG+CxH55LGZwpV2w5UCDkBO9zi7Q9sSc+m8D+vmQ5zfRe5wokFO+rtSjfncV+vT/LqLi\n+wjCMfv3bqDkVxCcU6P/zti3+6jZPWtl0h2N4WK70Fi3bh1qamr63eZyubB06VJcfvnlKZsYERFR\nrsmkOxrDxXahcXKRsXPnTjz77LPo7u5GXKx0f/WrX6G6ujr1MyQiIsoROXVH46SHHnoIt99++1ld\n74KIiIjOjKQLjQkTJmDp0qXpmAsREVFOy6mPTk66/vrr8cADD2DGjBlwu1X3r3/96ymdGBERUa5J\n6xeQpUnShcYvfvELFBQUoK9PlYd2OBxcaAB4c/dh7bHriv6zThx5Pv2FEfVFco6CUvsdxETmglN/\nq2JeX7/b4nE9+jsuskvicgwz6yRPpSk784qstpnFEjcyTTSlqkJcROwrZISgd42YarVL7v6J1Z59\n/Xtav9BnH1ltWUo9FtEzPPo6uq328X36e3L88/7r8OeV5PX7PABEgvbljl1eUX5XHFefUT49liDs\nXmadyMj9SK+eqZIo08ZOTIxnZtbIuZvZOWY5ebs5yAwfd756j819+Q6p8vwyAye/LF/rJ8cwj1/O\nySnm7i30aP3k40Kx30Bb0HY8t89+7pInX/+5014nMkPCRtaRzEKKR+V7ov/cyfdBzs9lpMWYr5MK\nqwqstrdIZbGZ77F8HA2L687IkIqKjKGiKvV7oaBC/z0WFu+XOYbd9WSS11eFmNPoGcZH9TmVdZKD\ndzQKCwvx1FNPpWMuREREOS0ng0G/9rWvYdeuXad8dOJ0ZuMNHiIiojMnJ+9o/PznP0cweOIWpMPh\nsL6trampKeWTIyIiyiU5dUcjHA7D4/GgsbHR9u7FyT5ERER0+nLqjsaKFSvwz//8z5g4cWK/2/fu\n3Yvvf//7+PWvf52yyWW6fe9/pD123vp3VjuWp0roxn16wKcjrILUYl4VbBU3SpU7YgnKaYu+Mtgy\nnqAEtaw0HDHKDmshfyq2DD6j3HNh5z41v64j+nz7VBlzZ58Kyjsa0y+zplbVr0CUOJ5Y9hWtX+mc\nGVbbLUtVG6Wlg6J0sRl4Wl2g9j2iVwWUOfu6tX5y7jLg1ZFXAL2jOh+Rlk/VHJr3aN16vmiFnbwy\nFWwbOh6w2tHeUH/dTxFPECQaDduXPtfLrOvBi3FRQtqR4GNR2c9bLAKIPfp7fHzfIavd26HOZ9GY\nEq2f06te19fh17bJoEQZNFpUowcoe0vUe3Ssab/VNoNBneJPRRlc6kzwJ6S32Ahe9Xn77SffR0A/\nZhkYbO4r3Nv/+2UGvCYKBi05R5X0zx+pzm/ErwduR8T1FTquznWoW+8ng0ZLa1WAd1FNhdYvLMYP\ndevHn+ga1fuJYFhxbY083/h/58kPBjReNsipOxpr1qzB3XffjerqasydO9eqAnrw4EHs2LEDBw8e\nxLp169I2USIiomyXU3c0pk6dis2bN+OVV17Bjh078NprryEej2PMmDG49tprUV9fz4BQIiKiYZRT\ndzSAE8Gf9fX1qK+vT9d8iIiIclbOLTSIiIgofbLxo5OUfvbR3NyMhQsX4plnngEAtLe3Y8WKFbj+\n+utx5513IhQaWPAbERERnZ1SdkcjEAjg4Ycfxpw5c6zn1q1bhyVLlmDx4sVoaGjAli1bzuovbOvY\np2cadBTfZrWPiWjy4+16ZHm+R5W/9nep0tUxI2PEkyAGJiwitHsj9hHeMrkkLCLBw7GBlQiWWSEA\nUF1cY7UnVY/UtrnbPrPazkCH1X7TSMAIhPvPpjnq1xeeRSIzQM7D5bRf8UeN4zrYrfoWeFQmkMep\nZwJ5fKqfp1BkZ0Afr1OUWo5VT7DaXSMXaP3ke3J+ZZG2bVKR2uZpVZkr8ZBRMrsvKLb19tsGgHjQ\nJmOmQM/wcIhUdFnS/ZTXefVMC62fKP/unTLTdjznK9vUtpA6Z4Vj9GvGO7Lcaoc7O7Vt8nX5Y1RJ\narlfAHDkq/PrfuFZq921t0Wfk1cdv7dYZaq4vPa/Bn2VZdpjT1lZv/36juoXeU+LepyoBLnM/nC4\n1DWeN0K/ZuTrnMYYJdMmW233mAlWO+bXM6vigeNWO9Kh5hc82qH1i4nMpZIpKvvDO0nPCovKrwUw\nryf5NQEiQ87MaJKZJpKnps545t/67ZeNsvGjk6R3NP74xz/immuuwRVXXIEFCxZg/vz5WLBgQdKB\nvV4vNm7ciKoqlR61a9cuXHnllQCA+vp67NixY+gzJyIiyjIuhyOt/9JhQJVBH330UYwePTpZV31g\nt1srWQ4Afr8f+fkn/loqLy9HW1v/X3hFRESUi7LxjkbShcbEiRMxa9asYdmZrCJ6spQ5ERERnZCN\nwaC2C42dO3cCAKZNm4Yf//jHuPTSS+ESnyHK2IuBKiwsRDAYhM/nQ1tbm/axChERUa7LqTsaGzZs\n0B5/8IEqAetwOIa00Jg3bx5efvllXH311di2bRvmz58/6DEySdwoEf7eIVV2++M21Q6G9H5eUdZb\nbjMDGQcT9DiQfrIdMgJIzZLkJ7mNOVSWqEDW6aOKtW1zxqpgsZZuFVD2yZEvtH5/PqCC/uS5GDPC\np/XziqC3PPfQEqTKi1TJ6CJRxtr8qyFqU7rdPM/meRuIL7r0IM+3ferOXlm+Cq4N2gTJAkBZsXqN\nUWUdfRH1Oo84Zx1Bvcx4e0AF244dq5/ry8aqwFF5zD6Pft4jYuev7Vfv47jK87R+l35VBK+KoFbM\n0YO/ww71nhR2fq5tcxxXJeNROMJq7naO0/q1BdRxzr9uldUu2PkHfbx8UTK9WI1nBr863KLM+Kha\nbVusUAWzxkU5+hIR1AsA+Z/YlMw2fmdEjSDKk1wjx+hzEkGUjjz9vXNMu9xqBwpUefb8WJ/WT34t\ngK9bnVvvF81aP/l+OS9YYLVb8/U/DEeIKuEuMR4AOKLqWos7xX8zDuPnOC5+nlziGje+tqEyT/2R\n29pn/3OSDXLqjsbTTz8NAHjrrbdw2WWXadu2b9+edOA9e/agoaEBLS0tcLvd2Lp1K9avX497770X\nmzZtQm1tLRYvXnya0yciIsoezlxaaHzxxRc4cOAAGhoasHr1auvLuvr6+vDII49g4cKFCQeePn26\ntViR+nuOiIiIAEcWfnZiu9BobW3FCy+8gJaWFjz++OPW806nEzfccENaJkdERJRLEn2b8NnKdqEx\nc+ZMzJw5E/Pnz09694KIiIhOn1nULRskTW+tra3Ft7/9bXz66adwOByYOnUqvvWtb6G2tjbZS4mI\niGgQcuqjk5Puu+8+3HTTTbjzzjsRj8fx3nvv4d5778Xzzz+fjvmdVV7/TEWQ7/xEFSNzGJkbMpsk\nIrIYYgPMJDHFbTImAGi1SmS/eIJ9yXk4jbl78tQl81HLcW3bvs7yfre989ERfXwxDxn41NauZ2e4\nvWpl7xKrfPM8JToutyhd7haZK4kyeuQ2M+tEPk7UT/J59TLuxTL7RWQThCJ6NL0c3+vWx7DrJwVD\nkX6fB4CPvXp56nf3qzLUMgPJzPY52q0yGQ51qvdrbHmB1q/tPJVRVlGgsjje+0Av1S3P2/lVennv\nC8eMVfv1q2N55KW/av1k5tJHteoavOPiv9P6QWQ/yOyRPk+h1k2+k83teubGnk/VeZMl/eeMm6b1\nm3a+nhlijW1kXTgPfKzaBSqLK16tj6dlqxjBgp/FVQZNi8h8M7/CoDRfjV9eql5TVaoXY3SG1TE3\nx9R52r1PLxFfKn4XTCqv0LbJxJ2B/lYTpxO93Xp211fPUVkomz85NsARz0459dHJSaWlpbjuuuus\nx5MmTcKWLVtSOikiIiLKDkkXGrNnz8arr76KSy65BPF4HDt37rTa8XgczgRf/EVEREQDZ37xXDZI\nutB4/PHHEY2eWiBlw4YNcDgcaGpqSsnEiIiIck1OfnSye/du3rUgIiJKg5wMBu3s7MSGDRvQ3t6O\nn/70p9i2bRtmzJiBysrKZC/NOW/8VQW6HRblmV1GupIMDpXBi4kCNE2xBAGgp0vO45RAVlni+kiP\ntm3/UfU4cFwFlLUf1gMPw72qPLXDqYIcPXlerZ9LBPmZ87BjVtWTr5PjJQq8lQGwpwSe2rzOfN5u\nvwDg9fX/Y5doDHlc0ageKGdeXyf19Yb7fR4ACkv0stt2AcBuowR5SJR/7ulU7+PxLj1oskuUO2/v\nUe3OVr/WTwYX/3uR/v6PGqHP8aSPP9a/9VnO/W8t6lorzqvT+l1crQIK3UF1jKEefe7vtnRZ7af+\n6zNt22EREOkVc39lhh5Q+T+/fr7VluXtu/v0927MhaqOt1tcut1h/VqQAbrmJRjqVYGysnx1jxEM\nLB8f9YuvQSjVv0qgwKvOU9sxFfDbZ5Tfb+lT779Z7r4or/9r3AyaliXz5e+0gFGOf/bfT1EPPnmr\n37GzRTamtyY9ogceeAC1tbU4fPiw9dzq1atTOikiIqJc5HQ50vovLceUrIPf78ctt9wCt/vECvWq\nq65CKBRK8ioiIiIaLIfTkdZ/6ZD0o5NwOIxwOGzVY2hvb0cwGEzyKiIiIhosZxZ+dJJ0oXHzzTdj\n6dKlaG1txTe/+U3s3r0ba9asScfciIiIckpOBoNeffXVmD17Nnbv3g2Hw4Ef/OAHqKqqSsfciIiI\nckpOLTT+8Ic/9Pv8m2++iVgshmuvvTZlkzpbHWhW0fA9Rw5YbXd+ke1r4rFTa5QMhN3rhjqeHZkV\nAgBOUVvY7dOPyy8yTfpEFLq/9fMB7ctp7Evu25zHQMnXybkPlTy/cuxYxD5uydyv09P/PMz3zm7u\nCfuJz1zDvQHbOfV06mW3oyKjQI7hMbIHev0q06CvS2VZBXuq9fG71Mer0YjKJuj1258nt1Gq/ain\n//e881Brv88DQLhPlTH/5YvN2ranRLZPWanKaBkzQi8XHgyp8xvXEy0QEh8bh/vU/HbvOar1u+YL\nVYK/RGTPzBYl0gFgxWxVZr2iQJSmN/6viYrHkaieueEVnct8HqttZmAd9ff/82lmeNQUq/nKMeTY\nANAmMouOGRlO8rHLYf8fp0dmVol2d5+eMTPmqxepB49ld9ZJTn108sYbbwAAOjo60NzcjAsvvBCx\nWAwffPABLrroIi40iIiIhllO3dH40Y9+BAC455578OKLL6Kg4MSXJvX09OCBBx5Iz+yIiIhyiPlF\nltkg6T2aAwcOWIsMACgqKsLBgwdTOikiIiLKDkmDQevq6nDTTTdhxowZAIAPP/wQ48ePT/nEiIiI\nck02VgZNutB45JFH8MYbb6C5uRnxeBwrV67EvHnzBjR4c3MzVq1aheXLl2PZsmV4+OGH8f7776Ow\n8EQg2ooVK7BgwYLTOoBMcujD1/t93u319fs8MPzBoKbYaQaHJgrQdHn1EtGBQlW6ONqnguaCHUds\nx3e4ZCDj0AI+B2qgwaXy2xPjsdiw9gPsAzsTBXnKdqJ+UjTUa/RT8zADeWPhkOgnysIb/SIhEQzp\nP97v8wAQ7FDXgpxvuFcvWy+3eYygaRlELYNtA+0tsBMR4we79cBLT766XtsL1HvwN6NEvAyMPd7W\noW3r3LdHzE8F1EaCNVq/Dpsp7v9rhfZYfm3B7LqRVvucigKt3+Ryta8JZfrPnd23EfQZpeo/aVPl\n34+KwO0RBXqQZ3SMGjBfnBszaFQGbJqlxbXxxLae3ojtNq/YV2dADy7NO/8S8WiD7b6yQbq/VM3+\nN9fwSbrQAIC5c+di7ty5gxo4EAjg4Ycfxpw5c7TnfvjDH+Lcc88d3CyJiIhyQE4Fg54ur9eLjRs3\nYuPGjdZzfr8/wSuIiIhyW6Z9dPKTn/wEb7/9NsLhMG677Ta8++67g/5kImULDbfbbX0/ykl+vx+P\nPfYYuru7MWrUKDz44IMYMWJEqqZARER0Vkn3RyeJvPPOO2hqasKzzz6Lzs5OXHPNNZg7d+6gP5lI\n69LpxhtvxD333INnnnkGU6dOxWOPPZbO3RMREWW0TPpStZkzZ+JnP/sZAKCkpAThcBjd3d2DPqaU\n3dHoz1VXXWW16+vrsXbt2nTunoiIKKNlUmVQ+cnE73//e8yfPx+tra2D/mQirQuNVatW4Xvf+x7G\njh2LXbt2YfLkyencfcqF/V39Ph8J9vT7PDD8JcPTycx2cHWrCH2ZJWAe45k65qGWMR9uiTJIBvIa\nk11Z9ERjO7v1bAq7vq48PWNKZqfI99gswd5rk2lkZkHFo+qxyyjNLrOaYhGVhRA2fp7kMctMm94u\nvVS5zLQJFpTCjjyWkPEzHZXXdUDNvdemrLwpFNDH6z6q5vTZh6rt9uZp/XxF6nHZKL18fJ4oDV5V\novr1RfR8gs8PqH33+tX5zC/Us056LlIZNLWVal8tHXpm0Vt729VxHNO3ObQy9vKahK1J49R7cqhT\nHy86d5L9C7NMJgaDbt++Hc899xw2bdqEt99+G7W1tZg0aRJ++ctf4rHHHsNDDz2U8PUpW2js2bMH\nDQ0NaGlpgdvtxtatW7Fs2TLcddddyMvLQ2FhIR555JFU7Z6IiOisk2nBoK+//jo2bNiAJ598EiUl\nJUP6ZCJlC43p06fj6aefPuX5RYsWpWqXREREZzWz7s6Z1N3djUcffRRPPfUUyspOfGHhUD6ZSOtH\nJ0RERGQvk2I0XnjhBXR1deGuu+6ynrv22msH/ckEFxpEREQZIpM+Ornhhhtwww03nPL8N77xjUGN\nw4VGGpzNAZ+JmMcV6bUPes0EmfI+DGUeA33NQPtFjZLhdswgTztmgGai0ur2czLK3dsEUScaT87X\nfL1DBMD25bWp541AWzlGpM8IShTnTb7ODHJ1iTLzsp/TbZT7FuM7eo7BjgwhPfKJ/dcCeEWQq9MI\nUJVB6fIcuoyvSPhP0d4lgkt7OvWS9p2tqgCjv1X/ok2zJP1JMpjYdOicaVbbk6f/1xTIn2m1vSLQ\nNJSg9PnZKpMWGsMl+46IiIiIMgbvaBAREWWITAoGHS5caBAREWUI+Y3W2YILDSIiogyRjTEaXGgQ\nERFlCCc/OiEiOtVAs1MSGY5S9cOddZNoPLnN/JqBqMgE0bJOjEwQp00/MxNmoMfV1yWzafT/sOKx\nmGj3X7Ye0M+NzEgxs8pkefaBvv+ylDygl3Rv/XiX1fYWl2n92oJXqm3MOjnrcKFBRESUIbjQICIi\nopRh1gkRERGlDO9oEBERUcpwoUFElAWGuxz9QANZEwWhOt1e220DZQZ2DoUsJy9rOsSj+jHJAFDz\neO0CWxOdp4jN8wDQE1KPZTBoNsqkL1UbLlxoEBERZQjGaBAREVHKZONHJ9l3RERERJQxeEeDiIgo\nQ2TjHQ0uNIiIiDIEYzTOgJ3/eseZngIREZ1B0cN7rfaSN587gzNJPSe/vTW9Zs2adaanQERElDb8\n6ISIiIhShgsNIiIiShnGaBAREVHK8I4GERERpQwXGkRERJQy/OiEiIiIUmY4vhgv03ChQURElCm4\n0CAiIqKUycKPTrLviIiIiChj8I4GERFRhnCwBDkRERGlDGM00quxsfFMT4GIiGhAhuX7ubjQSL85\n//T4mZ6Cxlc22mq3PXKJtu2/z/nvVjsaT9uUiIjoDLvt3eeHZRzW0SAiIqLU4R0NIiIiShkuNIiI\niChV+NEJIR6LWm2H23MGZ0JERFmHdzSIiIgoZbjQICIiolRhwS4iIiJKHcZoEBERUcpk4Ucn2bd0\nIiIioozBOxqDFOkLWm3niEpt2zifykKRlUGP9EVsx3M5VDud1USjcfudDcc85HG5HA5tWyjGsqlE\nRP1xZOEdDS40iIiIMgVjNIiIiChVeEeDiIiIUifdC41Y6nfBhQYREVGmSPdHJ1xoZJ5YJGS1I6XV\n2rbLzq2w2r6yfKv953cP2Y7ndapAya6w/o4PJWjSDPK0C+w0+9lda0MNGs0XxyWPEQD8UbU3ecyp\nDlAlIsp0aS/YZZ+rAAD4yU9+grfffhvhcBi33XYbLrnkEtx///3o7u7G6NGjsX79eni93oRjcKFB\nRESUKTIoRuOdd95BU1MTnn32WXR2duKaa67BnDlzsGTJEixevBgNDQ3YsmULli5dmnCc7AtvJSIi\nOls5Xen9l8DMmTPxs5/9DABQUlKCcDiMt956C1deeSUAoL6+Hjt27Eh+SKd/VoiIiGg4OJzOtP5L\nxO12o7CwEADw+9//HvPnz0cwGER+/onQgPLycrS1tSU9Jn50QkRElCky6KOTk7Zv347nnnsOmzZt\nwuuvv249H4/H4TAKMvYnrQuNPXv2YNWqVRg/fjwAYMqUKXjwwQfTOQUiIqLM5cisDxpef/11bNiw\nAU8++SRKSkpQWFiIYDAIn8+HtrY2VFVVJR0jrQuNQCCARYsWYc2aNenc7bCSWSfdeeXatnPm1Vrt\n0roaq105vUXrJ6OK3fkqWtd/uF3r19uhyp1HQ3peSDQctdp9XX2qfTyk9Yv0qpDiUCBstXsi+njB\nqMz+UM8dC5ASAAAgAElEQVSbmS+Jsj+K3GplW12Sp54fU6T1O7i3w2p/HlDzk3NIJDyI7JREmSyS\nzIyRY5ivN8upn2TOXY7hMl4ix8jWcuwDPe/6awbWL9H5HMp+UyGdWVLy2k30lQYDPTd2167dtT+Y\nsYeyX3Nb1sughUZ3dzceffRRPPXUUygrKwMAzJs3Dy+//DKuvvpqbNu2DfPnz086TloXGn6/P527\nIyIioiF64YUX0NXVhbvuust67tFHH8Xq1auxadMm1NbWYvHixUnHSfsdjcbGRtx6660Ih8O44447\nMGfOnHROgYiIKGPF035Hw/520Q033IAbbrjhlOeffvrpQe0hrQuNadOm4fbbb8eiRYuwf/9+LF++\nHFu3bk1a7IOIiCgnpH2hEU3e5TSldaFRV1eHuro6AMD48eNRUVGBI0eOYNy4cemcBhERUWYaQBbH\n2SatC43Nmzejq6sLy5cvR3t7O9rb2zFq1Kh0TmFYHTACL8fWjrHavjGj+m0DgMOj7uA4C4rVhkRp\nTQNMeYoHjmuPo10qwDTYqoIwg0c7tX697V1WO+xXwaXBjl69n3gcMwIgR05WwbFjr5hptfOnTNf6\nfbbp36x2yfa/qbF7k9TCtSEDKs3gSn2bet5r/NFgFwxqjmcXHJeofLxZgl2OYQbl6v36/4Uz1FLt\nZoDdUCQK2JP8NoG9ngS/RM0gX7t95Ru5/3JbrzjvQw1QlHMcaACkebR2Qb6nBlerxy5H/0GdJvN6\nKnSp8+ET7UQBygOdu/wqATm2KdFxJSKPuTemZlLh1f9rauruQ87g18SfnoULF+K+++7DSy+9hEgk\ngrVr1/JjEyIioi+lP0Yj9dK60CguLsYTTzyRzl0SERGdPXJxodHR0YGdO3eipeVELYjq6mrMmTMH\n5eXlSV5JREREg5JLCw2/34+Ghga8+uqruPjiizF69GgAQFNTExoaGrBgwQJ85zvfseqgExER0WnK\npYXGypUrsWTJEqxduxYulx6IGI1GsXnzZqxcuRK/+c1vUj5JIiKiXJBTMRrr1q1DTU1Nv9tcLheW\nLl2Kyy+/PGUTOxu8d0jP8JgyqsJqu0pHDmwQd4Jg2JjKb3YYiz0Zmezw5qunK0Zr3bw+ldWS71L7\nKjN25Qipqq2R9sNWOyraABDrVpkr0V496yavdqoab7oqSxsr0Pc2afz5Vru0boPV3rdtj74vkSZi\nZrhIcZlpENJzwuM2aRgOI6zfISLq42Jfdq83mdk5IVEi3uvR3zuXVz2WZeHNSH2PTfR5LGHWyeDL\nTCfuZ98xUUaGXUl7X4J0ikTl7uXLzOwHmYVxTLz/Q8l8MPdlZnjYSZTtpD9v/7qB7tc8heXiesrP\nV7/SzSyukJaRYz9X+bjUo851ods+8y0cG1iGi3mu5Ta533ETRmj9mhoP2e476+TSQuPkImPnzp14\n9tln0d3djbi4KH71q1+huro69TMkIiLKFblYR+Ohhx7C7bffflbXuyAiIjor5NIdjZMmTJiApUuX\npmMuRERElGWSLjSuv/56PPDAA5gxYwbcbtX961//ekonRkRElGtyKhj0pF/84hcoKChAX58qAetw\nOLjQALD9L0e0x7dUq0BM9+hz1AazfLgI8kShCpR0xI1IMfE47tTfqrgnT4yvtpkhaPKi1S5g42J2\neH1W25VXpNplVfqAET0AVCOOJSoCGWMO/fi7SmutdtG311vtS25s0vqF//quGq/jqBovrAe59XV2\nW+2Opv3ath4RsCsDQPNK8mAnkqAUusMmSK/vuP15cRn1zmXgaaKS7jJoVAalmv2cYjwzGFbfr5q7\nGeRqF/Rq7isqSq27PPa/EAPtwX7HzivRg5/luYgZc4+KgFqXCKg13zunOL9FbWq/4QTvo5yTGRjs\nEQGVTu/ASv+bc7fbt9lPBmzKoGGXsV8ZoGzOqWSsCvj2lanAcP+RgO2c5HjRkPkeqznK8QpG+mAn\n5A9rj+2uV6cRyCv3FROvGT3D+Kg+l4JBc7EEeWFhIZ566ql0zIWIiCi3ZeEdjaRH9LWvfQ27du1C\nKBRCLBaz/hEREdEwczjT+y8Nkt7R+PnPf45g8MTtSIfDgXg8DofDgaampiSvJCIiokHJwjsatguN\ncDgMj8eDxsZGOG0+MzrZh4iIiE5fNgaD2h7RihUr8Nlnn9kuMvbu3YsVK1akbGJEREQ5J5c+Olmz\nZg3uvvtuVFdXY+7cuVYV0IMHD2LHjh04ePAg1q1bl5ZJZqr339cjod3zLrTa0bKxVjvu0aO1HRGV\nwRMVGR5mxWSnyCGJG9XiZN+wLC1slBOWj+SmeILyzA6xL1+lvt/8zs9Vv079+B1hFfHv7FMlzQ9G\ni7V+H7epaPgCkblQXTxJ61d8qSppLpM9eiN6jFAwIs6TcViTS9UdN19rs9oQ06P/tYwf8f7AZX/H\nLrR3t9Xu+uADbVtfR4/V9pQUaNvc+SrzorddZcXEzWwSr7vfbWY/KRq2z7SQ4lGjVLuIu3LI8vYu\nM/tBvc7ptT83waOiVL3ItMgvs/8SxrBfL+Mus3+0rBNjDFe+ykKR+zWzh+wyhkxunzou+V6Z5PsQ\n9vdp20L+/rOQzAwP2c9bqPbl9tl/qm1m+5RNUb9rCqpU5lf353pWXKhb/UzKc21mjESC6ryVjC21\n2kXn2Bdt7Ovs0R7r16t9JpQks8kqvqL/LsCTHyBn5FJl0KlTp2Lz5s145ZVXsGPHDrz22muIx+MY\nM2YMrr32WtTX19ve7SAiIqIhyMKPThIGgzocDtTX16O+vj5d8yEiIspZ2RijkTTrhIiIiNKECw0i\nIiJKFTMeLxukdKHR3NyMVatWYfny5Vi2bBna29tx//33o7u7G6NHj8b69evh9doHW2W6lvff0B53\nV65V27pVYFN3lx5sVSBSgrva/bDjSrCwDYtyvVERAWnEgiIsSw2LbbEEwaBOcaHnu/VJjCkaY7Un\nja/WtrmPftLveK/t79AedwbV+XCJAL2W43pAXZEotWzOw45ZSfuL4yrorTR/nNX2GOPJc+3MV3My\nz1OXCDCMTlSl1HvG/oPWT57DyUbp5rqICqKVAbUOtxFcKUrXx/wiaDSkB03Kx/E+0Y7o150MgDW3\nxftUIK/DI34mjaDZqF+Ve3ePHG3b79h7H6rxxMmVwYqAHngqgxUBPWDRN1IFJeafM0Efw6vKZB//\n8/tW2wxQlPuSgYdOj/5rsHDMSKvtLdUDmaW4OOZwt17uOyCCUrVy3716kGiwvctq55WpfXkK8mEn\nr6xIe1wyc7bVdo2dYrV9n7yv9Yt1tVvtSI86N71iDoB+3ksn1aj9Tpmh9ZPXUEx8RYC5TZPo6xhE\n21N7vtbN6/yt1Q6Zv+Qo4yVdaPzxj3/Ek08+ie7ubsTjcatg16uvvprwdYFAAA8//DDmzJljPbdu\n3TosWbIEixcvRkNDA7Zs2cJvhiUiIvpSgr8Bz1oDqgz66KOPYvTo0cm6arxeLzZu3IiNGzdaz+3a\ntQvf//73AQD19fV45plnuNAgIiL6UqK7zWerpAuNiRMnYtasWYMf2O3WvlYeAPx+P/LzT9wSLC8v\nR1tb26DHJSIiylbZt8xIsNDYuXMnAGDatGn48Y9/jEsvvRQuUbxHfiQyULJc+cmPYIiIiOiEbAxB\nsV1obNiwQXv8gah86HA4hrTQKCwsRDAYhM/nQ1tbG6qqqgY9BhERUbZKVLX5bGW70Hj66acBAG+9\n9RYuu+wybdv27duHtLN58+bh5ZdfxtVXX41t27Zh/vz5QxonU/R2tWqP3zusouYbW1Qkt1kWXGZa\nBEV55kT9TGbfwfYLGWW8Izb98ozsjPIilZEwtVKPfp9eNdlq7+tUket7WvRSyPtFpo1bHGNpgZ6B\n5BNlpwd6Lsx+VSWqPHWpKC3tNO6myc9FZRaPef7keZP7SvTe7T2mZ1OU+1QJ7TLfuVY70KtnbkRE\ntkJJvshIMRK1AmH1uuMeFe1/pFvP4nF51Jxm1pZq2y4ao97LPnGM5l1HmZ3TGlCZG5UF+q+Sstj/\nZ7XdlSpzIfyV/0MfL6yyNQo7Dmjb4p3qunGMUOWvj5bq5anlFCsnfkW9vuOw1g+ykrEsue4zyqIX\njlDd8vSsE1lm3xFW5zevS8+68LW2oD9RkfkBAJHjKpvIXaCyk5ylI7V+MiPDWab/gdZ5njqnXX3q\nuMZfeo7Wz92jfl/l+4+puR/cp++qp1NtO+9Sq3286jytn/zd4O4yvo4gIjKjRF2IuFO/TmTpf1mo\nKlKqZ7RdWKp+jt/p0LOusk1O3dH44osvcODAATQ0NGD16tXWKquvrw+PPPIIFi5cmHDgPXv2oKGh\nAS0tLXC73di6dSvWr1+Pe++9F5s2bUJtbS0WL148vEdDRER0FsvCdYb9QqO1tRUvvPACWlpa8Pjj\nj1vPO51O3HDDDUkHnj59unVXROrvOSIiIsqxOxozZ87EzJkzMX/+/KR3L4iIiOj05VSMxkm1tbX4\n9re/jU8//RQOhwNTp07Ft771LdTW1iZ7KREREQ1CLHmXs07ShcZ9992Hm266CXfeeSfi8Tjee+89\n3HvvvXj++efTMb+zyn/tVYFeO/6qAq8cTjOgTj0OiWBQcyWbKP13uFe98QHer/PkqUvmwxF6meTd\nY0qs9v52FeS3+1M9AC4W6f9HyevTL0e3x9VvP/PYE81djpHntQ8uTRRQatcv0fNyDJ9XP44CbR4q\nAC4as/8V43UPPjA2GIrY9pPvDwC83KzOvQwMNoOBA+J6PeZXwZBF+Xr59L8/b4l6jQhW3fmiXqZe\nnoubL5qgbZtSpx63hlS/5z/Sg4tlQHV93XirfeG4Cq1fPE8EL8ufLeN6CjpUtO27h/RA3t6ICtis\nKFDHPK1uitavuHyvGl4EObp7jmn9nPuarLZ3kgpkjRaP0vo5RNBsrKhS2/ZZp3ofZIn81jz9uqsq\nHGu1y0rUeSouH6/18wZV+fQjxROt9pt/69T6eURk8KRyPXjVl9f/NWr+9Mhe8usD2o4GtX6zrlBz\nfOd//bXfsbNFFt7QSL7QKC0txXXXXWc9njRpErZs2ZLSSREREeWibIzRSPpNVbNnz8arr76KQCAA\nv9+P7du345JLLkE8HkcswV9gREREREnvaDz++OOIRqOnPL9hwwY4HA40NTX18yoiIiIarJwMBt29\nezeczoF9RTcRERENXTZ+TpB0odHZ2YkNGzagvb0dP/3pT7Ft2zbMmDEDlZWVyV6ak974RH1R3MG/\nnQiqqqkrx5HPVaXQ6gkjcHDficCqiupiHP1yW+W4ErQeUJUCq84ptbbJ9sm+J8ccdU5pv+1E28x+\nVWNLcPSL4/2221q6rX5jJpZZrxszIh9/+as63mljSvAf756oiHhebRnebzpRLdHtceHo5yqQrKK6\nBEe/+PK4xpZa7bGTR+LwPtVv7OSROPTlORxTW4bD+zsHNXcAqK4tsx6Pm1iGli/HP2diGT7/TAW9\n1UwYgQNfPp4wqdzads7EMuz7VAXwjZtYhs+/fHzOpHLrNdUTRqBFjHfOpHKr39TzKrFXBMRecF4V\nPt57Ytv5kyvQ9OW2KRPL0PSJ6nfu5JFWv69MrcSe5hPn+sJplda5BYCZ51ZZj78ytRJ//rI9pa4c\nH/1FBSWff16l9fjyi2uwe686rsvPq8IbXwZZXnpuFXZ9fKLfvOmjsLNZjXHhhHI0fjmP2ppia34X\nnz8Kf9mvjv/vzxuFX725DwBw4+xx+H+3nwgCnTFxJN744KDV76pLxuG13ScqeN58UQ3u3bzH2rbl\n/5yOe7buBwCsvmIi7v/Tibunl9aNxKZ//9jqd+NVk/G7bSfGr6+rwN3Pnvi6hJdvnYxvvaKuz3/5\n+yI0vHXiGL8zZzQ27j4x99uml+FXH6nr7qbpVbj/hWYAwPUza7Dmt+rrFx687itY/eQ7AIB//X8u\nxzd/8RYA4NXvLcB/36Lm9OQ8Dz73TQAAjOv9HPvyTlTprO05hkNVM61+Ffua4Lj8y2+wPtqM3poZ\nAADP8UPo9KmA0LLw39BVeiLTrzjag9a4rGbah4Ivg567eiNW++Tjk6oKvfjkyyDgS2qKsefoifac\nUuDzuKoUW4sOHCutO/EgFscXx0PWth4ReFvmc+LQl9VnJ5UX4K0v1Dm8YkIZ3jt04vfGRWOK8cHh\nHgDAhaOL8P4h9fvkojHF2H30RMDteZWF+Eur35rrX1p7VD8AExZ9+eWeDAY96zjiSe7TrFq1CnPn\nzsWf/vQn/Pa3v8W2bdvwu9/9Dk8++WTKJ9fY2Ig5//R48o4Z4msr/i+r3SJ+idtlTwD2GRjJpOur\nhM1S3U6RheA1otp9RapMcF+vKoV9TPxiAYCYTbSTeZ5cYl9m5s5AyTGdCcaQc5L9zLnKDBc5p4SZ\nL179jqDM3JHnNxrVrwU5vuyX6L2X8wj3nfqR50nFZT7tsfw1ILOdzLlHQmqO8j0uKNLropeVqoyk\nji5VMtp/XC+LnifKwheJcvEAMHmUKv8tM1z+8rFe+l+Sx3Xu+DJt21ixTZbSdxnX+PtiUfxnY19y\n/vLczJg+Wuu39u9UOX6ZTZHv0vdV6Vb/iUc9BWo/Yf1akJdub0R//9uDajEhFxY9CbKO5GKkxjjv\npeLnep/IaJH/8Z+Yh5qjxziuUptr3Lx2ZeaK1BYIaY+vfe8Jq7369n/r9zVn2m3vPj+kbzqXGhsb\nUT5+SvKOw+jY/ubTnncyST8T8fv9uOWWW6yvfL/qqqsQCoWSvIqIiIgGK57mf+mQ9KOTcDiMcDhs\n/ZXT3t6OYDCY5FVEREQ0WNmY3pp0oXHzzTdj6dKlaG1txTe/+U3s3r0ba9asScfciIiIcko2xmgk\nXWhcffXVmD17Nnbv3g2Hw4Ef/OAHqKqqSvYyIiIiGqRYFn5/q+1C4w9/+EO/z7/55puIxWK49tpr\nUzaps5UMAO34/G9W251f1F93AEA8Fu23DQAOp8t220DGGw5yDgDgdKsgOrdPPy5PntoW7lNxPL0d\nh7V+MZs5Op32QbNyv6ZExyxfl2iMgZL7kucmFrGPWzL368rz9dvPHEOOP5S5R0P6R5xyvEC3HpQZ\nE4GoMgjVDNDtC6oA0EhQBQcWjBih9ess9Pb7ml5/r9ZPXjOdefqvo6MtKnsoJiIq/cf1MWTwarhX\nlDtv1cusuz0qJK2gWAVAevP1/crA2FCfHlDZK4IUo8fVOXv9tb9p/a7asc9qF4lS/fdf9xWt3zfO\nVWXSZfn4qPFnbSBsv80rAjFlkKcZDNoWUO9DOKbe/yN+/bqrKysQ/dQxlvn0MvMyYLPbOE/mYzte\nEQzqEdddmzGnwpmXiUeZGQw6XDLtjkZzczNWrVqF5cuXY9myZXj44Yfx/vvvo7DwRObTihUrsGDB\ngoRj2C403njjDQBAR0cHmpubceGFFyIWi+GDDz7ARRddxIUGERHRMMukGI1AIICHH34Yc+bM0Z77\n4Q9/iHPPPXfA49guNH70ox8BAO655x68+OKLKCg4scrt6enBAw88MNR5ExERkY1MuqPh9XqxceNG\nbNy40XrO7/cneEX/ksZoHDhwwFpkAEBRUREOHjyY4BVEREQ0FJkUo+F2u63SFif5/X489thj6O7u\nxqhRo/Dggw9ihPHR6SnjJNtRXV0dbrrpJsyYcaJi3Ycffojx48cneRURERFlmxtvvBG1tbWYNGkS\nfvnLX+Kxxx7DQw89lPA1SRcajzzyCN544w00NzcjHo9j5cqVmDdv3oAmNBxBJGeTg39+y2pHjEA8\nt7f/AMBUB3naBV4OhgzSlAGFLm++1k8Gh0b79OPv7WrDQDhc/QeEJgoUTXSMdnM3H9sFeSbiML4D\nKC4C5+R4TrceRCcDO8331e59NoNBB3o9REMqcFIel6ewRB+vny9OtPqK9zUsAkBlsGn3IcBXpqpj\ndot9hXv1ipKxsAr0k++3t6BU65dXqr7mQAaentjfp1bbJX62vIX6GPKY5fXZU6gHMssqtFFR8bL7\n8Odav94uVSnUW1yubZPXmrwme46o5+/+H5+gsm6S9fhP01UGX6cI1gQAlwiOnFip5rts1litX74I\ncpUBoAUeFz5pV7e5G0XF05BRkbjUCPS0nhcBunkuJwJhdVxhEaAbNQILtMBW0e4yjlFu83n1n7uj\nMmD5vAv6nV82yqSPTvpz1VVXWe36+nqsXbs26WuSLjQAYO7cuZg7d+6gJjNcQSTZwG6RcTZI9B/8\nQJ3uIiMVBrqYGCi5yBj6GMO76JSLjKHy+OwzpiS5yEhELjKGSi4yEhnoeywXGUM10J8TuchIxDXE\nkvuSXGQkYrfIMMlFRqodNbKickm6vl5iqFatWoXvfe97GDt2LHbt2oXJkycnfc2AFhpDMVxBJERE\nRLkimkFf37pnzx40NDSgpaUFbrcbW7duxbJly3DXXXchLy8PhYWFeOSRR5KOk7KFxnAFkRAREeWK\ndN/RSHQ/b/r06Xj66adPeX7RokXDto9hd+ONN+Kee+7BM888g6lTp+Kxxx5L5+6JiIgyWjQeT+u/\ndEjrQuOqq67CpEknPqOsr69Hc3NzOndPRESU0WLxeFr/pUPKPjrpz1CCSM4mQaPU9klhf1eaZ5Ie\nicqTy4DFROW5zxRz7kPJOhmuffc3h4G+Jtnr7ISMa9JuDLNcugzm1N9jPZtAbovavMbUl6cHDcuf\nJxnYGjEymhy9KvYr1NNhO768Pr3FZbb95DFGevW4Mnkt94kg50QZSDLrqHWvHsj6XwdarLZZMl56\nX3yNwWetegZOsSihLjNXgr1GCfKD3VY7IgI7vUbp9y5RMn7ORJVZY5YFf2tvu9U+JMYG9Nv/5vha\nP5F1UjtWZUId6dQDmf2F6v8KUXEd0cyOmxySTIrRGC4pW2gMVxAJERFRrsj0rJOhSNlCY7iCSIiI\niHJFuuIm0imtH50QERGRvUz6UrXhwoUGERFRhjCrrGYDLjRoyMzAvkTBbJkmUVDiUMu9D8e+h/M1\niQz0vRpoIG9vV+uQglwTzUmWO5djmOMNJXg3bJQ0l+QxJzp+syy8HRlQa87dlaBqsOwrg3ff+0/7\nsvXmVx9I8lhkyXmnRz+Ot0Jq22eH9SDPYI8ao+OIOoc9bUe1fn09x/qfg1EZVs69ddKFVtvh0Cuj\nHuxRga1eo2pqMMsiQrMxRiOt6a1ElJ1SnamTa1K92B0KucjIFNm2yMhWvKNBRESUIbJx7cSFBhER\nUYbIxo9OuNAgIiLKEAwGJSIiopThHQ0ion6kInhxoBkvQ9l3pNc+62SghlJa38y6cHpUiXNngoBa\nGWzbc2Sfbb+BnotEJfePfPKR1e46Umm1zePt61aZJYmymGQ5dnMMWZ7+6F/ettqybDsAHOz+qtX2\nuVQOQzCaeUGzp4sxGkRERJQyvKNBREREKRNjjAYRERGlCj86ISIiopThRydERDQoiYJGZRBlouqq\ncluifsMRlBvpU3MKdhzRxxfBl1pJ8wRl4WXb7CcfJzoXbQG1r0IRDFroc+JAMGxzJGentH97qyN5\nl9PFhQYRESUVz8AMj2xbZABnIEYjDd8ewIUGERFRhsjGGA1+qRoRERGlDO9oEBERZQgGgxIREVHK\npD0YNA0yfqGx81/vONNTICKiMyly2GpOfv13Z3AiqccvVUuzWbNmnekpEBERpQ0XGkRERJQyXGgQ\nERFRynChQURERCnDhQYRERGlDBcaRERElDJcaBAREVHKcKFBREREKcOFBhEREaVMNi40+KVqRERE\nlDK8o0FERJQhIll4R4MLDSIiogyRjR+dZPRCo7Gx8UxPgYiIaECG4/u5uNA4Ay5f+cSZnoKtI7dX\nao9X3/5vZ2gmRER0Jt327vPDMg6/Jp6IiIhShnc0iIiIKGW40CAiIqKU4UKDNE6Pfvp8LofVjopr\nJZSFF86ZJk61dq6JiM5m0VgszXt0JO9ymrjQoLOOK/U/F0REZ0T672hwoUFERJQz+NEJERERpUw2\nVgbld50QERFliGgsntZ/yTQ3N2PhwoV45plnAADt7e1YsWIFrr/+etx5550IhUJJx+BCg4iIiE4R\nCATw8MMPY86cOdZz69atw5IlS/Dcc8+hpqYGW7ZsSToOPzo5Dd4xY7XHU4ryrLas7rYvEB7QeIkq\nwrkc9gE7iV4nMzKGUnHO5xraWnQomTaJ5ifnYQaD9kQGFqUtz4U5BjNX6HR4nfY/n8w6o8HIpBgN\nr9eLjRs3YuPGjdZzu3btwve//30AQH19PZ555hksXbo04ThcaBAREWWITFpouN1uuN36MsHv9yM/\nPx8AUF5ejra2tuTjpGR2RERENGiZtNDoj8fjsdrxeByOBHfbT+JCg4iIKENk+kKjsLAQwWAQPp8P\nbW1tqKqqSvoaBoMSERFliHgsntZ/gzVv3jy8/PLLAIBt27Zh/vz5SV/DOxqnwTvxfO3xtNpSq+3O\nV6fW1dSu9Sv32q/vDvdG+30+UTVMGQxpBp7JMEm5LVHg5TifR3tc6lHzlUGTLcGI1q9XlM71JZiw\nx+ZWWzjBnCq8+qUqh2/uUelV5nHJID15/KcG1/a/b7OfPobelwGlZ6dEP1uJ3lN5bSUKBpXXpDme\n3b55LeWuWAbd0dizZw8aGhrQ0tICt9uNrVu3Yv369bj33nuxadMm1NbWYvHixUnH4UIjg9gtMs4k\nucjIFJlYgpz/MRDRcIgPITswVaZPn46nn376lOf7ey4RLjSIiIgyxFA+zsh0XGgQERFliEz66GS4\ncKFBRESUIeLp/pb4NOBCg4iIKENkUozGcEnrQmPPnj1YtWoVxo8fDwCYMmUKHnzwwXROYVhFqiZr\nj6tn11jtEZNUu/qSVq1fweiRVtvlUW/BsY/3a/1iIRUc6vS6tG3ufK/V7trfYbU7PuvU+gWOBVW/\nsBrPjJKXJb49eWpfpeNLtX4OEYk59khA27b/ULfVDkbtI+1rfOqY5TycRoaHP6LmO7KyUPUzokEP\nfEX6rfQAACAASURBVKwq00Xj+jY5vsthnzESsvkrIlHpdzl2MBqz3ZYo00COn6hUdaJsl6Ewj2so\n5en115vjD2y/sl9wGCJq9ffbvl+icvR2/U487v99SDRGomuhyK1+7uQ1lOhcZEq2k5zHQOcw1Ayf\nXMKPTk5TIBDAokWLsGbNmnTuloiI6KyQjcGgac1d9Pv96dwdERERnWFpv6PR2NiIW2+9FeFwGHfc\ncYf29bNERES5LBvvaKR1oTFt2jTcfvvtWLRoEfbv34/ly5dj69at8Hq9yV9MRESU5WIMBj09dXV1\nqKurAwCMHz8eFRUVOHLkCMaNG5fOaQyb1nih9rioptJq549UQZTeYr2fp6RAtUepYy+dMcN+Z259\nMeYsLLbaDrEt2nFU63f8z+9b7bC/12rHwnr58PyRJba77uvosX2dVNOigl4793eJ/Ya1fp5CVeK8\ndKzab/7IYtjxlujnsLROBdvu+r9/Y7WPR/SgzAoRRCu3Fbr0Tw1lIJ4MWDSD1zpEQK0si348old1\nLXHrwbuSHFMGA3aF9bnLwEM5d3NOcr4DDeo0y8DbZdSZAapy33Lbqf1UxxJxjImCQeW5BewDNhMF\nDcrS9/lO+0+G5ftlBkbLa6PXOC4z6FfNz/64EpUnv7A032p/Kkrpt4XMc6HmYY4nz728nkx25y1R\n2X4pUWB0InoArf0Y8tyOytP/a5JfM5DteEfjNG3evBldXV1Yvnw52tvb0d7ejlGjRqVzCjRIcpGR\nKeQig4gom3ChcZoWLlyI++67Dy+99BIikQjWrl3Lj02IiIi+xPTW01RcXIwnnnginbskIiI6a+Rk\nwa6Ojg7s3LkTLS0tAIDq6mrMmTMH5eXlKZ8cERFRLsmpEuR+vx8NDQ149dVXcfHFF2P06NEAgKam\nJjQ0NGDBggX4zne+g8LCQrshiIiIaBBy6qOTlStXYsmSJVi7di1cLj2CPhqNYvPmzVi5ciV+85vf\n2IyQ/Zra9BLcs6rKrLZ3pCozHo/qEeRSPCIyMpz6eXa4PbAVU8temYEi2wAwcmydeuBS48XyjAWi\nQ0SriyW1I9xnO4X4sUPa4+5dr6v9nq9Kn0eNTJVCUYLdO2qM2pc3X+vnLFX9nBVjrfah0ilav0t/\n+h9We//nx7VtlQXqmKMycyNBdL4scX48oGfMnBNX49mVjwb0ku55xUbGkBjfna9+BPuO65H1IbHv\nYyH7LAn5OFEZ80SZC3YZCeZ48rF8jdlPju8Vp7rSyCaQZef3+Y3jtzmWROW5R+ern6HqkjxtW0y8\nrlH87JqZEBeU+6x2njHG0aOq6GBLUF3XiTI35HTl/ADgopsvtNr5z+622geD+nU30PL05V77bKdE\nY0h2WSfm84kygQZaJt/uGppQO0Lr19yo/67JZjkVDLpu3TrU1PQf3e9yubB06VJcfvnlKZsYERFR\nrsmphcbJRcbOnTvx7LPPoru7WwtS+dWvfoXq6urUz5CIiChH5GTBroceegi33347610QERGlWE7d\n0ThpwoQJWLp0aTrmQkRERFkm6ULj+uuvxwMPPIAZM2bA7Vbdv/71r6d0YmeD1z9r1x5fXqFSfl2j\nzlEbYkYwqAj6dFcacTCibyyoAs8cRkBurHKi1Q4XqCBUMzfK1X2k/21O/a2Pi8cO0S/u8mqvi+Wr\nYNN4hZoDABSdc4Ha1vyW2Jc+d3eVKrsux9MCUgFERqgA0OYucQ77Ytj+WZv18LYHrrPaeY9t0caI\niiDKqCjxHTeCNx0ieFMGa+KI/o3DMqCwoEIFDY4s1AM+3T51PvNK9G1e0dcpgvf6juuBt3KOZQdV\nhVanUT7dIeYrj9f8y0huk0GoZt9YgmDLiAiADIvgTfM1cfG4dLwqM19QUaD1087Tbr18fjQkgndF\nRGlPR6/eT+y6Ypza14R6PWjY5VH7cv32PastjwnQA0AnXX2etq1OvCcf/dsHVrvJeO9kMLC8nC7+\nx8lav3PufdBqB47ea7VLPu3Q+gXaVHC1OV9PkQpQdnnU9eQ0IjLDvep18v2R59YcQ15bLiPQVPaL\nGT9P5jVqR75OXoNVF1TpHRkMelZLutD4xS9+gYKCAvT1qR8kh8PBhUYqmAuSTJCBSd1ykUFElE1y\nKr31pMLCQjz11FPpmAsREVFOy8bKoEnvb33ta1/Drl27EAqFEIvFrH9EREQ0vOKxeFr/pUPSOxo/\n//nPEQye+HzQ4XAgHo/D4XCgqakp5ZMjIiLKJTn10Uk4HIbH40FjYyOczv5vfJzsQ0RERKcvnomx\neqfJdqGxYsUK/PM//zMmTpzY7/a9e/fi+9//Pn7961+nbHKZ7qXGFu3xg3//FasdH6ui1eNGNoUj\nJiL3RcZI3MgEcUT1MsRSzKPKdfeI6Py+iP6xVvGICVY7P66yBOIuPROio09klogFdaLg8e6gvq9Q\nVJUNHj/zH612Xs8RrV/EqzIPovkqS8Bh/IB94RcR/q09sNN0yT9Z7YvWlGnbYl0qMygeUcff296l\n9fMUqgySuPho0OxnRtdbrzGeL6qptNp542q1be5q9TgeVnPqfvMVrV9Q7LtojDoOh/GmeIv7/76h\nsD+oPY72yiwRfb7ymD2F6tpyGH9k9HWq9+HI7lb1emO8yvMqxHgqi6OkdozWr0CU7S8Y+bExfxWA\n7spXf9D0HOrW+kXD6ropqlLnYtS8S7R+8lgur1PZXp3NB7R+oW6VaVRUU6FtyxuhsqS+IrJYen+1\nS+tXPkrNQ2ZuTPjGQq1frFCV2a+7YZHV9vzpVa1fb6d6L3uNrJv8MvV+uTzqGN0+/Y9AmdUkM1fM\nrBOZdSOvNXldAIAzwR+ZWhaXOO+JPnaP9qr5jTzf+H/nyQ+QK3JqobFmzRrcfffdqK6uxty5c60q\noAcPHsSOHTtw8OBBrFu3Lm0TJSIiynY5tdCYOnUqNm/ejFdeeQU7duzAa6+9hng8jjFjxuDaa69F\nfX297UcqRERENHiJvoTzbJUwGNThcKC+vh719fXpmg8REVHOyqk7GkRERJReXGgMUnNzM1atWoXl\ny5dj2bJlaG9vx/3334/u7m6MHj0a69evh9frTT5Qhtr7zp+1x+G7V1ntlm4VyBk10pWcUMFhflla\nG2HImDqPUUJY6upVAWtdfSK41AjKKxWlpkvzVPBWIBzQ+rUFRKCgUTCmNL//oK/uPr0Uclgc5xG/\nek25b4TWLxBQx7yvQ1X5vKquXOs3qlDNvViUPy4eWYjmdnX8u1o61RymfkMbQwbHTipXAZ8e41O/\noCjJ7BXnvaLjE71joBP9KtCPMVqiSigHvPq2j9tVMN9fRJBr3d9dpPWb26PKZIcPqHmYZesd+Sq4\nVpaqD3/erM9JBMY6C0tgJy5K30cCekCpu6jIapfW/U29xrzuJqk5uorV8XvG6WXBHSPUecob974+\n33ZVdlrO9/jH+nsSOq6u5bJzx6v9fvU6rZ/7qDiH4isC8qddqM/JbR/kGPOrQNRR1WpfM43zVPGV\nSeo1IfVz4p14vj6g/IqAy6+1mmPa9eq3wVZVkryvQw+MlgG2ThGg6vTqv97D4jzJgFdT/shSqy0D\naN0l+jXjyFM/T2ZVY4dXBI6KryAwg4tlEHI8pH4uvBOmaf3Kvc9a7WOh7PuPWMrGhUbKgiwCgQAe\nfvhhzJkzx3pu3bp1WLJkCZ577jnU1NRgy5YtCUbIPTYJDWeU3SLjTJKLDCIiymxJFxp//OMfcc01\n1+CKK67AggULMH/+fCxYsCDpwF6vFxs3bkRVlfqLZdeuXbjyyisBAPX19dixY8fQZ05ERJRl4rFo\nWv+lw4Aqgz766KMYPXr04AZ2u7VvewUAv9+P/PwTt9TKy8vR1sYvxyIiIjopGz86SbrQmDhxImbN\nmjUsO5NVRE+WMiciIqITYrm00Ni5cycAYNq0afjxj3+MSy+9FC4RaCZjLwaqsLAQwWAQPp8PbW1t\n2scqREREuS6n7mhs2LBBe/zBB6oErMPhGNJCY968eXj55Zdx9dVXY9u2bZg/f/6gx8gkx7/Qo/rf\nP6yCFP98WEWnO407NzKroytoX2bc5VSvMzNXgiLyWm4LGSXIi0TWiU9kbgSNyO2eXhUZHxDb5BzM\nx+a+pLEiw8McQ76uvUdlu3zeqUfuL52uPq6rLVPj7Tmil6CW3j2olwzvCqjz23hQHf+s6lKtX0CU\nsfaI8skVBRO0fqVVagz5voaN9ycUUo9dET07p03MSe73L0aZ9SMF0632qHMv6vc1ANAt3q/jYuwD\nLj07Zeo0lTGyaJKe4SPfoWBEzT0Q1t/jMUXqehpZ+b+stswQAQBXqSqtLUuuHxurlwWPip+F8kvP\n0bY5o+racIjy8eWX6Mcf+8sbqh1Q10Ysr1jrd7BCZZcU16j3sSDQqvVz9qiPdON5enl3l/g6gbhL\n3aGdMHaS1i9y+HOrLc+Nlo0BPeuir0Cds4KrbtD65TU3Wu1Yj5755B4zQY3vVll8Dp8+d2/HUavt\nE21Hvt7Pc47IDBL/6cV9etZJ3C2OJRax3SbPE9x6lqH8OgZHVLV7Ciq1fvPHq5/XzZ8cQzbLqYXG\n008/DQB46623cNlll2nbtm/fnnTgPXv2oKGhAS0tLXC73di6dSvWr1+Pe++9F5s2bUJtbS0WL158\nmtMnIiLKHjlVGfSLL77AgQMH0NDQgNWrVyP+5V8efX19eOSRR7Bw4UK7lwIApk+fbi1WpP6eIyIi\nohy7o9Ha2ooXXngBLS0tePzxx63nnU4nbrjhBruXERER0RDl1EJj5syZmDlzJubPn5/07gURERGd\nvpxaaJxUW1uLb3/72/j000/hcDgwdepUfOtb30JtbW2yl2Y984L4z72qxPNrf1UBZjIIE9CDN82g\nzLgRVGgnKsqIRhMEZbrcKnitwKeCsvqM/UZEgGHMLJkugjnltogRKCj7fVykgr7MYNA+Ubpczl3O\nFQCOiUDR20Sg4FviPANAgTi/Y0p9sLO/XZVg/t+NLdq2kAiGnXyOKpm9eLpeP0YG8spgULNse5cY\n73ivHvD7/n5VTrp5rwpsG1GpB+WVFqj3S14zBcb1JHX3Rmy37WlRgbKvG+dw/MgCszsA4ONDeuDt\nRDFHr1sFdlaM04P8FtSqYNOaYnUcL/5Fr51TJq7JfOP933tMnTeXONc1JUVav6umzrbanl413xB0\nL32mzrUsTd9jlNK/froKhsxzG4HMolS9fE9qqvTjdxfqZedPikeM4G+vGsMDNadjxRO0biO+ooIh\nvUc/1bZFq1QgqgwolXMFgMJxIsg7oK7BuFd/77vd6j2W34LQHdJ/3oPyZ9dhf55CYdWOhfU5uRwq\naNQjfk98euC41m/69TOt9uYfvoxsJgOEs0XShcZ9992Hm266CXfeeSfi8Tjee+893HvvvXj++efT\nMb+cMtBFRjo5nax1QkSULjl5R6O0tBTXXae+nGjSpEn8jhIiIqIUyMaFRtLvOpk9ezZeffVVBAIB\n+P1+bN++HZdccgni8ThiWXiLh4iIiIZP0jsajz/+OKL95PVu2LABDocDTU1NKZkYERFRrsmpEuQn\n7d69G05nyr5NnoiIiL6USQW79uzZg1WrVmH8+PEAgClTpuDBBx8c9DhJFxqdnZ3YsGED2tvb8dOf\n/hTbtm3DjBkzUFlZmeylOeeNT1REfctnKqrb7dUXanHxiZPM4oglyB4xySyHRK9zikh+lyitbQae\nyiwWyWEEgw50jnK/3jz9Mgt09/W7rc/IzvAfV/2+Wqei6fd+qmdMuD0qC6OzRi+TXCxKsMvS6m0H\n9WyK9v2qZPT+3SqD4NMDE7V+/+2ycejPhwf0stByX2Zm0WdNKiNJnt9ev378XcV6JsNJid4Dcwyp\npEJlFxw+6te27RbZL6LK9imZRX87pLIB8n39Z8UAwO4W1S8oMpr+/LFe7lsyM3cCXer99+Sr9zgv\n36P3+8b5Vvu/Vaj3+1kjw6Xhqff63W/ZqKJ+nweAj1r07IdPD6jMHXmeLpw4UusnM4OWzbraap9/\n7F2tXyxflUk/FFDn2nyLQ26VxVIxYba2LSw+AT/Yo97/FvHzAwDFXnVuxpfq89XGE9klraKkfbPI\n2gL0jJRS4z2RDomfd7N8vsw0yhPtA8bXEVzxtcvVg6zPOsmchUYgEMCiRYuwZs2a0xon6a2KBx54\nALW1tTh8+LD13OrVq09rp0RERHSqeCya1n+J+P3+hNsHKulCw+/345ZbboHbfWI1fNVVVyEUMjPU\niYiI6HRl0kIjEAigsbERt956K5YtW2Z9q/tgJf3oJBwOIxwOw/FlQZb29nYEg8EkryIiIqLByqSP\nTqZNm4bbb78dixYtwv79+7F8+XJs3boVXm//H+vaSbrQuPnmm7F06VK0trbim9/8Jnbv3n3an9cQ\nERGRzuv1Yue/3pH2fdqpq6tDXV0dAGD8+PGoqKjAkSNHMG5c/7FqdhzxeDxpOcojR45g9+7dcDgc\nuOCCC1BVVTWonQxVY2MjLl/5RFr2NRwmzb/Gah87sN9qu/Ptg81ikYF9DGWucuXrEqVDOZ0qKM3h\ntC9dHRXjOY1+8nWJ+tm9xjz+QLsq/+0rG93v8wBQXK1KK3vyVanijn16SnX3ob1qTm79h6Zy2mVW\n25WnypOb57OvSwUpyvPpMsYrqqqx2pMvGGW1Q0b0Xmer+myz9YAevNq5b48a36uOS84PAApGqn3J\n+Ub77O8ohnt7rLb5/uSLc53ny9O2mWXnrfH69OvT5VZ/mxSNUPM1X19YosaXgcfHO/S5yzLuvQH7\nnwWnCGQ2S9XX1Kly51OqVTDwzvf06+noXvUzKa8TT2Gp1m/sZBUoaQbD7v9zs9XOK1UB8ebPcX6h\nCryVwcof3jVB63eoUH2Vw79/ol8nUoUoR3/RmGJtm7z0ekVQ9wdG+XgZbCtLv08q10uQh8QYR8TX\nAPTZBIwDpwZ57jmogmjNnw1Jfj2DLLl/uLNX6/fIDDXGHef8o+14Z9Jt7z6PWbNmnelpDKvNmzej\nq6sLy5cvR3t7O6677jq8+OKLw3dH4w9/+EO/z7/55puIxWK49tprBzdjOislWpwQEVH2WrhwIe67\n7z689NJLiEQiWLt27aAXGUCChcYbb7wBAOjo6EBzczMuvPBCxGIxfPDBB7jooou40CAiIspixcXF\neOKJ0/9UwXah8aMf/QgAcM899+DFF19EQcGJ22s9PT3/f3v3H9PUucYB/AstRVYuImNW1BvGCFMj\n3mBGTBwCTmZYokuWTZNtIRnL4pa0srGhC4tjMtAAxiwKy2ZmNrZgYqgsms0tWWTLMtEpCzFT2Lwk\nzOv4eR2dF0pbKD96/zD2vOdgawVPaTnfz1/l9PSc97znVJ+cPs9z8O677856x0RERDT/3bW8tbu7\n2xtkAEBsbCz6+vpUHRQRERHND3etOklNTcULL7yAjIwMAMDly5e97UiJiIiI/LlroFFVVYVz586h\ns7MTHo8Hr776KrKzswPaeGdnJ8xmMwoLC1FQUIDKykpcunQJRqMRAPDKK69g48aNszqAUNJ7+aL3\ntdshZV1HxQRWdaJMvAy0nnpy3He2foQusGROX/31dVHyxB9xX4FuO2qBUfb3mF1qzy4ev9spb/es\nE6pQpoT9TrrllQs6g7xaQ/T3tV+9r8UqDGV1iihSL2W/Kxt6jwrVKUN9UrVCdGyCbD2xgkRZJTIh\njN/tEFtay+dz3Md70yuQpFFOCPtSnh/Xzf96X+sV50TcprgvZTWF+J5DUa0hGhYqMsRtKyuLxPMw\nNiRvGS6OUTznyvMdFZ3hfd0jVG6MOuQts8V9i5Umyuup/9qd2/YD8vMvHteEYhvjDuk7Lx5jhy5H\ntl5zh3ROLgvtzQ2Kyhqd0Kq+Y0BeTSISKzx+VbTFF9vExwsVHpkp8mtX3EaiURr7PxfK531c2F6U\n4lEFaYul4x8XqlX6huTVJD1CFdL/hHbnfzvk7dPHEtNB4euugQYAZGVlISsr65427HQ6UVlZifXr\n18uW7d+/H6tWrbq3URIREVFYUu2xrAaDAUePHpX13LhffdOJiIgoPAR0R2NGG9brvc9Huc3hcKC2\nthZ2ux0mkwllZWWIj4/3sQUiIiIKd6rd0biT559/HiUlJTh27BhWrFiB2traYO6eiIiIgky1Oxp3\nsnnzZu/rvLw87N27N5i7V52YbCdSJtSF0kNz7kaZKDcTU4pkVV8JoMqkSTEZ0jMlJZQF2rZduW+P\n2ErdTwKtXkjkFBNDleNw2aQyb8df3bL1xARY8TPTxufnWMYCvE58bUOZXCrOr5iE6o/yWhW3OTZy\nU7m69J79b2F8UpKfv2t/wk/SrPwaku/3ry4p8VD8jDK5VDzm8VHpZ9zpycrS2JVt4SeEz8la8yu+\nJ+Jci9vYc/o32Xp2h3Rctn6xfbw8uTJhiXSMnX3ypOkx14S0PSG50jEsT7wUE1sNQgvyaz3y7U0J\nyaARwjgWGOXfhQihffy/khfJ3ktOlLc1v21S0aq+S0hsTRTa1g8Oy5NB/xzmE8PDWVDvaJjNZvT0\n9AAAWltbkZaWFszdExERUZCpdkejvb0dNTU16O3thV6vx3fffYeCggK8+eabiI6OhtFoRFVVlVq7\nJyIiohCgWqCRnp6OhoaGacvz8/PV2iURERGFmKD+dEJERETawkCDiIiIVBPUqhOtCqcqEzX4q6xQ\nVqSIJt1S1vyUMIf+qlj8UVZhBPKeWDEBAJPCviZcI8rVpfXGfFfrBHo93Et1zb3uZ6bXpPg5f+3z\nR8VqnwD35a/CxR+x4ke8Tny11Qfk58ffeso27uK1569Vu3jtiq+vtP7pcxz+Wvo77VIVhiFa/s/2\nyP+kyo2hP6WqFuX8ia3QYxYt8b6+Ke8K7/O8Gv4hb1UuVsb857cbPseuj/J9XI7BAe/rRcuW+1yv\n44bv7xqFPt7RICIiItUw0CAiIiLVMNAgIiIi1TDQICIiItUwGZTmlJgAqEwG9NUm/H4kMvoz7ifJ\nU0yUE7enTLzTWgKw8njvx/HP9nwFmgzr7xEB/s5roG3cI4TreLjn37L3IqOkBE0xKVU5JqesLb5B\n9p5baJk+qmi7LhuHcCxiK3UlsY17pPAZ5X5FekUbd2Xr9tuU58Rtl9rJu53SfCr3dfXGCmnbQnf2\nSXlHcwpRvKNBREREqmGgQURERKphoEFERESqYaBBREREqmGgQURERKph1QnNKX+VAb5acKtR0eGv\n+sWXQFtkk7pmez3MtFW7vxbsvrYxNnLT53p+OYel7Slapk+6fbe79zWOQCtmpsQ/FG31xWP2V6ml\n91GBAiha+vs5jus2p/d1lM+1KFTxjgYRERGphoEGERERqYaBBhEREamGgQYRERGphoEGERERqSbC\n4/GEbLf4tra2uR4CERFRQB577LG5HkJICulAg4iIiMIbfzohIiIi1TDQICIiItUw0CAiIiLVMNAg\nIiIi1TDQICIiItXwoWohqLOzE2azGYWFhSgoKIDNZsPbb78Nu92OJUuW4ODBgzAYDHM9zKD44IMP\ncPHiRYyPj2PHjh1Yt26dJufC5XKhtLQUNpsNTqcTFosFGRkZmpyL20ZHR7FlyxZYLBbk5uZqdi7a\n29thNpuRnJwMAHj00UdhNps1Ox9ff/016uvr4fF48MYbb2DNmjWanYtQwTsaIcbpdKKyshLr16/3\nLjtw4ACee+45WK1WLFu2DF999dUcjjB4fvnlF/z+++9obGzEZ599hqqqKs3OxQ8//ID09HQcO3YM\ndXV1OHDggGbn4raPP/4Y8fHxALT7HQFu/ZuRn5+PhoYGNDQ0oKysTLPz4XA4UF9fj+PHj+PIkSNo\nbm7W7FyEEgYaIcZgMODo0aNYvHixd1lrays2bdoEAMjLy0NLS8tcDS+o1q5di0OHDgEA4uLiMD4+\njgsXLmhyLrZs2YIdO3YAAAYGBmAymTR7XQBAV1cXurq6sHHjRgDa/Y4At/5zVdLqfLS0tCA3NxfR\n0dEwmUzYt2+fZucilDDQCDF6vR4LFiyQLXM4HN5lCQkJGBwcnIuhBZ1er4fRaAQAnDhxArm5uXC5\nXJqci9u2b9+OXbt2oaysTLPXBXDrDkZpaan3by3PhdPpRFtbG15++WUUFBTg559/1ux89Pf3w+Vy\nYefOnXjxxRc1PRehhDkaYSAqKsr72uPxICIiYg5HE3zNzc2wWq2or6/H2bNnvcu1OBcnTpxAR0cH\n3nrrLeh0Ou9yLc3FqVOnkJmZieXLl3uXafk7snLlSrz22mvIz8/H9evXUVhYCLHhs5bmw+12o6en\nB4cPH0Z3dzcKCws1+z0JJQw0woDRaITL5UJMTAwGBwdlP6vMd2fPnsVHH32ETz/9FHFxcZqdiytX\nruDBBx/E0qVLsXr1akxNTSEmJkaTc/Hjjz+ip6cHZ86cwcDAAAwGA6KjozU5FwCQmpqK1NRUAEBy\ncjISExNx48YNTc7HQw89hIyMDOh0Ojz88MOIjY1FZGSkJucilPCnkzCQnZ2N77//HgBw5swZ5Obm\nzvGIgsNut6O6uhqffPIJFi1aBEC7c3Hp0iV88cUXAIDBwUE4HA488cQTmpyLQ4cOoampCVarFdu3\nb4fZbNbsXADAyZMn8fnnnwMAbDYbbDYbtm3bpsn5ePzxx3HhwgV4PB7YbDZNf09CCR+qFmLa29tR\nU1OD3t5e6PV6mEwmHDx4ELt27YLT6URKSgqqq6uh18//m1GNjY2oq6tDSkqKd1l1dTVKS0s1Nxdu\ntxvvvPMO+vv74Xa7YbFYsHr1apSUlGhuLkR1dXVYtmwZNmzYoNm5sNvt2L17N4aHhzExMQGLxYJV\nq1Zpdj4aGxtx+vRpOBwOWCwWrFmzRrNzESoYaBAREZFq+NMJERERqYaBBhEREamGgQYRERGphoEG\nERERqYaBBhEREamGgQbRPHb+/Hns3LkT3d3dyMnJmdE2iouL8dNPP93nkRGRVjDQIJqnHA4HwLzr\nggAAAlVJREFUysvLUVFRMau2y++//z4qKysxMjJyH0dHRFrBQIMoDNXX16OsrAwA8Mcff+Cpp56a\nFgg0NTUhJycHCQkJsuUDAwPYunUrrl69it27d6Ompgavv/46Nm3ahG+++QYlJSV4+umnUV5eDgBY\nuHAh8vLy0NTUFJRjI6L5hYEGURh66aWXcO3aNbS1taGiogIVFRWIjY2VrdPS0oKsrCzZspGRERQV\nFaG8vBwrV66ETqeDzWZDbW0tioqKUFFRgb1798JqteLLL7+E0+kEcKu1s/hAOyKiQDHQIApDkZGR\n2L9/P4qLi5GWloZ169ZNW6evrw9JSUnevycnJ1FUVIStW7ciMzPTu3zt2rUAAJPJhEceeQRxcXGI\niYlBXFwc7HY7ACApKQm9vb0qHxURzUcMNIjC1NDQEB544AH09fXd8X3l0wWGhoaQnp4Oq9XqvVMB\nQPbcB+UzIPiEAiKaLQYaRGFobGwM7733Ho4cOQKDwYBTp05NWycpKQn9/f3evxMSElBSUoInn3wS\n+/btu6f99ff3Y+nSpbMeNxFpDwMNojB0+PBhbN68GSkpKdizZw/q6uowMDAgW2fDhg04d+7ctM8W\nFRWhq6sL3377bcD7O3/+PLKzs2c9biLSHj69lWieGhkZwTPPPAOr1Tqt8uReDA0NYdu2bTh58uS0\nhFMiorthoEE0j7W0tOD48eP48MMPZ9xLo7i4GM8+++yMG34RkbYx0CAiIiLVMEeDiIiIVMNAg4iI\niFTDQIOIiIhUw0CDiIiIVMNAg4iIiFTDQIOIiIhU83+xZaPBYFNIwQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7ff18f2190b8>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"path = 'lock_exchange_kh0.01/prog_{}.nc'\n",
"cmap = cm.RdBu_r\n",
"\n",
"fig = plt.figure()\n",
"\n",
"ax = plt.subplot(311)\n",
"d = Dataset(path.format('plm'), 'r')\n",
"pcol = plt.pcolormesh(d.variables['xh'][:],\n",
" d.variables['zl'][:],\n",
" d.variables['temp'][-1,:,0,:],\n",
" linewidth=0, rasterized=True, cmap=cmap)\n",
"plt.axis('tight')\n",
"ax.invert_yaxis()\n",
"ax.xaxis.set_ticklabels([])\n",
"plt.ylabel('depth (m)')\n",
"d.close()\n",
"\n",
"ax = plt.subplot(312)\n",
"d = Dataset(path.format('ppmh3'), 'r')\n",
"pcol = plt.pcolormesh(d.variables['xh'][:],\n",
" d.variables['zl'][:],\n",
" d.variables['temp'][-1,:,0,:],\n",
" linewidth=0, rasterized=True, cmap=cmap)\n",
"plt.axis('tight')\n",
"ax.invert_yaxis()\n",
"ax.xaxis.set_ticklabels([])\n",
"plt.ylabel('depth (m)')\n",
"d.close()\n",
"\n",
"ax = plt.subplot(313)\n",
"d = Dataset(path.format('ppm'), 'r')\n",
"pcol = plt.pcolormesh(d.variables['xh'][:],\n",
" d.variables['zl'][:],\n",
" d.variables['temp'][-1,:,0,:],\n",
" linewidth=0, rasterized=True, cmap=cmap)\n",
"plt.axis('tight')\n",
"ax.invert_yaxis()\n",
"plt.ylabel('depth (m)')\n",
"plt.xlabel('x (km)')\n",
"d.close()\n",
"\n",
"plt.tight_layout()\n",
"\n",
"fig.subplots_adjust(right=0.82)\n",
"cbar_ax = fig.add_axes([0.85, 0.2, 0.05, 0.7])\n",
"fig.colorbar(pcol, cax=cbar_ax)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.colorbar.Colorbar at 0x7ff18ee83358>"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhoAAAGACAYAAADrrMueAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt4VdWdP/73ueTkHiBAAkEbQxDwK05BBEWKQQKT+SHj\nWEDxkt93YinQX/BSBSwtN4U6EMRLURgrVcZipwX1h6UtUwQcvhoMYqNMSUsnI1YKCSA5hpCcc5Jz\n/f4R2XutTc4ll71zcvb79Tx5nrXPWmftdTa3D3t/PutYQqFQCEREREQ6sPb2AoiIiChxMdAgIiIi\n3TDQICIiIt0w0CAiIiLdMNAgIiIi3TDQICIiIt3Ye3sBkVRXV/f2EoiIiGIyfvz43l5CXIrrQAMA\nJn13iy7z9vvGddLx+WduVdpr/u5/S30X2gK6rIGIiBLDgj+83dtLiFt8dEJERES6YaBBREREumGg\nQURERLqJ+xwNvQR9Xuk4lJSqtJMsFqOXQ0RElJB4R4OIiIh0w0CDiIiIdMNAg4iIiHTDQIOIiIh0\nY9pk0OazJ5E6YIhyHOifp7Qz7HL8ZRNyQwMh3ZdGRESUMEx7R0MMMoiIiEgfpg00iIiISH8MNIiI\niEg3DDSIiIhINww0iIiISDemrTrxupqk4y9DGUr7uoGp8tigWmpyvs2vtD0sQSEiIoqIdzSIiIhI\nNww0iIiISDcMNIiIiEg3DDSIiIhIN4Ymg9bU1KC8vBz5+fkAgJEjR2LVqlVGLkER8HpgsdqU44/r\nm5X2xAdvkcb+nfi+Vq/SPvfpKWlc3bHzSvvPl9RxDV6/NE5MLo24RuaaEhFRH2dooOF2u1FSUoIV\nK1YYedoOiUEGERER6cPQRycul8vI0xEREVEvMzTQcLvdqK6uxoMPPojS0lJUVVUZeXoiIiIymKGP\nTkaPHo1FixahpKQEp06dQllZGfbt2weHw2HkMoiIiCgKj8eD5cuXw+l0wu12Y/HixRg7diyeeOIJ\nNDc3Y8iQIdi0aVPUf8MNDTQKCwtRWFgIAMjPz8egQYNw/vx5XH311UYug4iIKO4c+s1/IDMvx9Bz\nOhwO3HDDDR32vffeexgzZgwWLFiAuro6fOc738HYsWMxZ84czJw5ExUVFdizZw/mzp0b8RyGBhq7\nd+9GU1MTysrK4HQ64XQ6kZuba+QSFKFgQDr+9R/PKu0759wj9QVdakVKyK9Wk+Tc3iyN+zthzn/w\n+5S2/6sL0rjmv6nVKU0n65T2xVPytuhNpy4p7bpGj9R32qNWsrT4g0o7EOpaRUu497HyhYjIGJl5\nOdh20xxDz7ngD2+H7bvjjjuU9rlz55Cbm4ujR4/iqaeeAgAUFxfjjTfeiK9AY/r06Vi2bBneffdd\n+P1+rFmzho9NiIiIvmaz9PYKrnT33XejoaEBr7zyCh544AGkpKQAALKzs9HQ0BD1/YYGGpmZmXj5\n5ZeNPCUREVGfYbPEX6Tx5ptv4k9/+hMef/xx2Gzq1hChUAiWGNYbNdBobGxEVVUV6urab/Hn5eVh\n0qRJyM7O7sayiYiISCue7mgcP34cAwcORF5eHq6//noEg0GkpqbC4/EgNTUVDQ0NyMmJnlMStrzV\n5XJh9erV+Kd/+iccOHAAjY2NaGxsxMGDB3HXXXdh9erV3BeDiIioB9ksFkN/Ivn000/x+uuvAwAa\nGhrgcrlw++234+DBgwCA/fv3o6ioKOpnCntHY+HChZgzZw7WrFkj3SoBgEAggN27d2PhwoX4xS9+\nEfUkfcHHf1CTMlv+abrUlxJSE0AhJHxavW5pnCUobDUuJFcmBbzSuDSvmtg5xN+qvsXbKo0L+eT3\nSX2t6rkDTU6l7Tt3WhrXUqcmonq+bFTarU45kbXlSzVobKlvUdpfnZeDyXOt6mf0aDJFY00o9Qnj\nxPdot2YX3xeQ3tPhaTol1qTZyHN0fx1ERKJ4uqNx77334oc//CHuv/9+eL1erFmzBtdffz2WLFmC\n7du3o6CgADNnzow6T9hAY+PGjRg2bFiHfTabDXPnzsWtt97a9U9AREREknjK0XA4HHj22WeveH3H\njh2dmidsoHE5yKiqqsLOnTvR3NyMkPC/wNdeew15eXmdOhkRERGFF093NHpK1GTQ1atXY9GiRb22\n3wURERH1XVEDjWuuuSbqZhxERETUffH06KSnRA007rnnHqxcuRJjx46F3a4Ov+uuu3RdGBERkdkY\n+k2nBokaaPz0pz9FWloa2tralNcsFkvCBRpn//yp0n7nL2OlviEZyUrbZlWjzX7J6dK4lCT1t0iS\nMM7hkCPU5FR1nEN4IJek+R0mnksb5Vqh5stYhKqW5IBPGpfsUytZLH7h11CskAGAULDjtoZVqJgJ\nNp6X+vwX1Modv/Oc0m5rlCtc2i62dNj2fHlRGtd8Vu1rFiphWhvl6pzWQPgt2MNVrkSqcBH7tJUl\nYp92Dnk+tS/81Yxc/RJu7T0xrid0pXKHlTpEkZnyjkZ6erpSR0tERET6MWUy6G233YajR49e8ejE\nak3EGzxERES9x5R3NF566SV4PO23yy0Wi7K3+YkTJ3RfHBERkZmY6o6Gz+dDUlISqqurw969uDyG\niIiIus9UdzTmz5+PJ598EsOHD++w/+TJk3jqqafw85//XLfFGcl1Qd26e+vv/iL1OVLVy2Szq0FX\ncrJ8+VId6lbtaUK7X5pDGpcpvK9fmhqoie8HAIdwrhSbHOwl24VzCVmkGQ55TWlJqeocdjV51ar5\nzZwkhNEW4VTa3/RJQmKro99oeU0jxMRWYZwmRE8V5hSnt0RILhSTVy1tcnKpzaVuwY6mL6U+cXv2\nYLOabBpqlbdWD7jUOcXk1VbnJWmc+5w6n9upzuFrkbeLDwqJoiFNBmRQSF4NeoVEVl9AGhfwBoS2\nOs7fKify+j3qsS+ojot1S/dItMmb4baCjzQu1uTaWNcUKbk23Bqu7Ivtfb2VQGtk0mxPbMcfeX5d\np084prqjsWLFCjz++OPIy8vD5MmTlV1A6+vrUVlZifr6emzcuNGwhRIRESU6U93RGDVqFHbv3o33\n3nsPlZWVeP/99xEKhTB06FDMnj0bxcXFTAglIiLqQaa6owG0J38WFxejuLjYqPUQERGZlukCDSIi\nIjJOIj460fXZR21tLaZPn4433ngDAOB0OjF//nzcc889eOSRR+D1eqPMQERERH2Zbnc03G431q1b\nh0mTJimvbdy4EXPmzMHMmTNRUVGBPXv2xM0XtoWCaob/yY+rpT6rXa0asVptHb6uPbYI42zJqdI4\nm7DxmT1JnE+O+2xCpYnFKke5Nrt6LM4htrXH4nusmioWqzC/XViHQ7OmzBR17RkpcmlzhlBNI1bd\naOcQj6Vt1rWfUThOEapsUjTzZTgGqudNypH60oao70saps6XbJOvk3g5rGGqYrTSbGK1j7ymfsnC\nr4m3ReqztKnHFp+6pbtF2N5d2xcUqmKCLfJW7YFGtdLG52xQ2t5mubIm6FWrU0JBuXZDrIQJiVUx\nPrnCRZwjIPQFWtukcQFxXKu6Lb4vQsWMtupGrNYR16clVu5InyOorfYR5hMqerRrCoap9gGAYJgK\nja5uER9u6/ueqBiKJNY1RarA6co6emI7/ljf0+KXz9ZXql8S8dFJ1Dsav/71r3HnnXfi9ttvx9Sp\nU1FUVISpU6dGndjhcGDbtm3IyVH/4j969CimTZsGACguLkZlZWXXV05ERJRgbBaLoT9GiGln0A0b\nNmDIkCGdm9hul7YsBwCXy4WUlBQAQHZ2NhoaGjp6KxERkSkl4h2NqIHG8OHDMX78+B45mbiL6OWt\nzImIiKhdIiaDhg00qqqqAACjR4/Gs88+i5tvvhk24bm2mHsRq/T0dHg8HqSmpqKhoUF6rEJERGR2\nprqjsXXrVun42LFjSttisXQp0JgyZQoOHjyIWbNmYf/+/SgqKur0HEZoOfdFTOPEhM+OjpXXNYmH\n1nDjYpyvvc/a4Thtgqp47sjzddynXauUGKs5l5j0KiXDah6haRNbldc1kbyYoComyto0iax2IRFT\nmwwrjo00h7jtunhe7VrFNYrb0WekaLZ+F5JhtUmzqcIak6XEWzlp2CFsGW+zDhbamqTZfsJ27wPV\n+fprzjskI1lpD81MlvrEZNZkm9iWzyVuJy+2rX45GVRMZLV6mjpsA0DIoya5hgJyMqhETF4NapJG\n/WqyacjbqrZ9clWb1Ce2PXLSrN+trj3QqtlaXlijVfizFdSsXUy2jZRcKybN+l3qmnxuOTHYJ/SJ\nyaqRWDS/dmJyrZh4G9DMJ21pLyTKiq+3z9HxZ4zkyu34O16T9jOKSbnBMNvbA/I/0p+1+KS+k66+\nUeVoqjsaO3bsAAAcOXIEt9xyi9R34MCBqBPX1NSgoqICdXV1sNvt2LdvHzZt2oSlS5di+/btKCgo\nwMyZM7u5fCIiosSh/R6qRBA20Dhz5gxOnz6NiooKLF++HKGvo8i2tjasX78e06dPjzjxmDFjlGBF\n1NFrREREdOWdqEQQNtC4cOEC9u7di7q6OmzZskV53Wq1Yt68eYYsjoiIyEysZgo0xo0bh3HjxqGo\nqCjq3QsiIiLqPost8b6sNGp5a0FBAR599FF89tlnsFgsGDVqFB566CEUFBQYsT4iIiLTMNWjk8uW\nLVuG++67D4888ghCoRA++eQTLF26FG+//bYR6+sVQX/fyE6+LFI1iZ5zdKWKBbiyCicW4Sp1tOe6\nsnKn4+qcSHNEej1cFY8txu3oI6030nWJde3htrfXHieldLU6x9Jhn0XznzBrmMoi7db3GWlpSjt/\nYLrUlz9I7csRKmaykuW/tjKzhK3vk8Rt4eVxYpVMqvAZxfe0v0+ougnIlQsxC3VcJWMJypUbEP6u\nsQjnsvjc0jCLT63quWKOcLS/KOKaAh1X6gBAUKjCCbWq6wi6LsnTie/TVgKFqyDSbH0fEj6/363O\n5/fIVUxyn/oebVWQVfg1H/Tun6W+k0fqOl5TnDHVo5PL+vXrh7vvvls5HjFiBPbs2aProoiIiCgx\nRA00JkyYgEOHDmHixIkIhUKoqqpS2qFQCFZr4j1PIiIi6g2WBPw3NWqgsWXLFgQ6uA22detWWCwW\nnDhxQpeFERERmY0pH50cP36cdy2IiIgMYMpk0IsXL2Lr1q1wOp14/vnnsX//fowdOxaDBw+O9lYy\nSCgY25bEes9hBrEmzfZEgm5Pr8GWFD5BNdz7Yr2N29XPKybNHknSbGkv9NkcqR2+B5CTaMW+WLe+\n1xKTZtOz5K3aM7PVdVw9WE1eHZyZIo8TtqQXt6fP0CSyigmrA1LV5Nfc9P7yfGlCwq9m58gg1G24\nhd25odmdG+LHt4tJvZrL4hASdsNtOa99X6y7WWqHWYW/d+zCNvYWv5ygagmTNIuAJnFfSIDtP/JN\nqetXH72mvk1zbeJJIpa3Rv1EK1euREFBAc6dO6e8tnz5cl0XRUREZEZWm8XQH0M+U7QBLpcLDzzw\nAOxf/+9gxowZ8Hr7VvknERFRX2CxWgz9MULURyc+nw8+n0/51kqn0wmPxxPlXURERNRZ2r1mEkHU\nQOP+++/H3LlzceHCBXzve9/D8ePHsWLFCiPWRkREZCqmTAadNWsWJkyYgOPHj8NisWDt2rXIyckx\nYm1ERESmYqpA45133unw9Q8//BDBYBCzZ8/WbVFE8SrW6px4rOLRJPLHhXjcPj/SOLFy5xOrWO2S\nFPZ9PfEZrcJ5U/sPkfqyBg9U2mJVTEq6vKa0VPVYrIrJzpAra4b2Uyto8oT2wDS52idJuMUfDIUv\n4xArUrT/hiYLW8GnSNvCp0njHLYMdQ7hctrs8oQBYR0FJd+R+vJSdijt054ubi1vAFM9Ojl8+DAA\noLGxEbW1tfjmN7+JYDCIY8eO4cYbb2SgQURE1MNMdUfjmWeeAQAsWbIEv//975H29ZcftbS0YOXK\nlcasjoiIyETCfSFhXxb1Hs3p06eVIAMAMjIyUF9fr+uiiIiIKDFETQYtLCzEfffdh7FjxwIA/vjH\nPyI/P1/3hREREZlNIu4MGjXQWL9+PQ4fPoza2lqEQiEsXLgQU6ZMiWny2tpalJeXo6ysDKWlpVi3\nbh0+/fRTpKe3b907f/58TJ06tVsfgIj6rr62fX48JNReOlMrHX8ZJtlUm4Ta3aRU7fuTMwco7dQB\ncoJq2sBcpZ0iJJFqE1SThQTVfv3VxNOrsuVkUDFBVdzSPdUhr8nrDyrt6YWDpL4b+qlJr/GdDGrs\no5Ng9CHdFjXQAIDJkydj8uTJnZrY7XZj3bp1mDRpkvTa008/jeuuu65zqyQiIjKBREwG1e0ejcPh\nwLZt26Q9N1wul16nIyIi6vMsNquhP0aI6Y5Glya225XvR7nM5XJh8+bNaG5uRm5uLlatWoX+/fuH\nmYGIiMhcjH50YgTdAo2O3HvvvSgoKMCIESPwyiuvYPPmzVi9erWRSyAiIopbRn3RWayee+45fPTR\nR/D5fFiwYAH+8Ic/dDrX0tBAY8aMGUq7uLgYa9asMfL0REREcS2edgb9+OOPceLECezcuRMXL17E\nnXfeicmTJ3c619LQQKO8vBw/+tGPcNVVV+Ho0aO49tprjTw9EVHCCVd1o3c1jr+1RWm7LpwOO64n\ntmC3OdSt1dNzrlba2d+Q/w1Jy1IrSy7c2ib1zRmnVsLgP1oQr+IpGXTcuHF44YUXAABZWVnw+Xxo\nbm7u9Dy6BRo1NTWoqKhAXV0d7HY79u3bh9LSUjz22GNITk5Geno61q9fr9fpiYiI+px42kdDzLV8\n8803UVRUhAsXLnQ611K3QGPMmDHYsWPHFa+XlJTodUoiIqI+zWKNn0DjsgMHDmDXrl3Yvn07Pvro\no07nWhr66ISIiIjCi6ccDQD44IMPsHXrVrz66qvIysrqUq5lfH0iIiIiE4unfTSam5uxYcMGvPLK\nKxgwoH0n2PLycpw5cwYAYs615B0NIiLSTU8kpYqJp01/O6G0tduxJ2dmK+3/TJ0h9T108wj14D9O\ndntNeomnHI29e/eiqakJjz32mPLa7NmzO51ryUCDiIiIrjBv3jzMmzfvite//e1vd2oeBhpERERx\nIh6TQbuLgQYREVGcsNi6v+9IvGGgQUREFCfiKUejpzDQICIiihNWPjohIiKKD9qKFq+rSWk3npe3\nGc8uHi0c7dNzWd3COxpERESkGwYaREREpBtWnRAREZFueEeDiIiIdMNAg4iIKE6JyaGe5mapL+nq\nkUrbZlFfD4R0X1anxNuXqvUEBhpERERxgjkaREREpJtEfHSSeJ+IiIiI4gbvaBAREcWJRLyjwUCD\niIgoTjBHoxdU/Wxxby+BiIj6uP8S2t/5+O1eW0c0Vn57q7HGjx/f20sgIiIyDB+dEBERkW4YaBAR\nEZFumKNBREREuuEdDSIiItINAw0iIiLSDR+dEBERkW4sVpa3EhERkV4YaBAREZFuEvDRSeJ9IiIi\nIoobvKNBREQUJyzcgpyIiIh0wxwNY1VXV/f2EoiIiGLSI9/PxUDDeJO+u0WXeYdNmCkdf3rTJ0p7\nxff+XeoLhHRZAhERJYgFf+iZb4TlPhpERESkH97RICIiIt0w0CAiIiK98NFJAqn7eC+GjpuuHNvT\nU3pxNUREROAdjUQiBhlERERxgYEGERER6YUbdhEREZF+mKNBREREuknARyeJFzoRERFR3DDtHQ1v\n81fScdKwQqXtsFqkPg+3BiUiIgNYEvCOhmkDDSIiorjDHA0iIiLSC+9oEBERkX6MDjSC+p+CgQYR\nEVG8MPrRCQMN/Tg/+wRpA/OUY+t1dyjt3GT5sjT6Akq7xa/+qjBHlIiIepLhG3b59T+FaQMNMcgg\nIiKKC8zRICIiIt3EWaDx3HPP4aOPPoLP58OCBQswceJEPPHEE2hubsaQIUOwadMmOByOiHMw0CAi\nIooT8fQ18R9//DFOnDiBnTt34uLFi7jzzjsxadIkzJkzBzNnzkRFRQX27NmDuXPnRpwnfj4RERGR\n2Vltxv5EMG7cOLzwwgsAgKysLPh8Phw5cgTTpk0DABQXF6OysjLqRzL0jkZNTQ3Ky8uRn58PABg5\nciRWrVpl5BKIiIjilyV+/v9vt9tht7eHCW+++SaKiorw3nvvISUlBQCQnZ2NhoaG6PPoukoNt9uN\nkpISrFixwsjTdsjruiQdn00arLTv/u54qS+5f6bSbqm7oLRPHz4tjftzXYvS/sLtVc8VjK08hVUs\nREQmF0eBxmUHDhzArl27sH37dnzwwQfK66FQCBaLJcI72xn6iVwul5GnIyIiom744IMPsHXrVvzs\nZz9DVlYW0tPT4fF4AAANDQ3IycmJOoehgYbb7UZ1dTUefPBBlJaWoqqqysjTExERxbWQxWroTyTN\nzc3YsGEDXnnlFQwYMAAAMGXKFBw8eBAAsH//fhQVFUX9TIY+Ohk9ejQWLVqEkpISnDp1CmVlZdi3\nb1/U0hgiIiJTMPzRSSBsz969e9HU1ITHHntMeW3Dhg1Yvnw5tm/fjoKCAsycOTPqGQwNNAoLC1FY\n2P517Pn5+Rg0aBDOnz+Pq6++2shlEBERxacYch6MMm/ePMybN++K13fs2NGpeQwNNHbv3o2mpiaU\nlZXB6XTC6XQiNzfXyCUo/K0t0rfk/edfG5X2fffeLY0Vxw32qHkm+XddlMZ9y68mgIbaWpV2y+lz\n0rimk3VKu/GvX6ntz+X5zjS4lfYXbp88h0/cCl3NIrXF+JtUfE/7cUxvIyIiPcXRPho9xdBAY/r0\n6Vi2bBneffdd+P1+rFmzptcemyTiV/ESEVHfFi1voi8yNNDIzMzEyy+/bOQpiYiI+g4zBhqNjY2o\nqqpCXV377f68vDxMmjQJ2dnZui+OiIjIVMwUaLhcLlRUVODQoUO46aabMGTIEADAiRMnUFFRgalT\np+IHP/gB0tPTDVssERFRQjNToLFw4ULMmTMHa9asgc0m5zMEAgHs3r0bCxcuxC9+8QvdF0lERGQG\npsrR2LhxI4YNG9Zhn81mw9y5c3HrrbfqtjC9hYJy7fD29z9X2vf8f3dIfRahQsPib1PatrYWzTi1\nEgTC/Nn+VmncoDa1ciXk9ajtNo80DkF1vpBQ0QIAgUZ1K3Rv3SmlfemLs9I411mn0vY41fO6G+Rz\nub5U+85dUKtd6jx+aVxrjNupS2vVVLiIW7KL1S5drYTRvk9PrM4hIl2ZKdC4HGRUVVVh586daG5u\nRkj4C/21115DXl6e/iskIiIyizjaR6OnRE0GXb16NRYtWtRr+10QERGZhpnuaFx2zTXXYO7cuUas\nhYiIiBJM1EDjnnvuwcqVKzF27Fjle+kB4K677tJ1YURERGZjqmTQy376058iLS0NbW1qEqTFYkm4\nQOOzP3ymtHefKJT67Db1Fz4tSa3A6ZecKo1LsavjUpI6bgNAcqr6DC7JqrbtVvnZnE04tGn6HEJn\nspCgmuOTkzwtPiERVUhQtQTlJE+IiaxC2+Jrk4dd/FJp+8/+VerznlcTUT0X1C3dvZfc8rhmISn1\ny0tKu7leTq69dKZZaTc1q+uIlJAZa+KpNxg+8VScI9ZxQYQXKVk11s/SE+Ni1dPJtUygJeoEM25B\nnp6ejtdff92ItRAREZlbAt7RiPqJbrvtNhw9ehRerxfBYFD5ISIioh5msRr7Y4CodzReeukleDzt\nt+MtFgtCoRAsFgtOnDih++KIiIhMJQHvaIQNNHw+H5KSklBdXQ1rmGdGl8cQERFR9yViMmjYTzR/\n/nx8/vnnYYOMkydPYv78+botjIiIyHTM9OhkxYoVePzxx5GXl4fJkycru4DW19ejsrIS9fX12Lhx\noyGLNMJXn/+X0v7Jb74RdlxSsl1oy98B43Cox6lCu3+aQxqXIcyRmaK2HXb5Fz01zHwAkCnMkeFQ\n22JVDACk2NUvvRMrXJJschWL1SJWwqjrcKTI41KGXau0k/O/JfWJlTAO4VwpEappxPNqN8SzCFUy\nVq9aqWJtuSCPu3hOafsbzkl9Ibda1RJ0qVUsAY9cCdN2Ue1ra1SrX9znnNI415fqOK/Lp55HU1oR\nDATD9onHAZ/6GQNeeVv8gDfY4Ti/Zlt4r9AnVslEqpiJJNYqmVi3j+9KFQ8QvpKnq1Ux4T5/Z6qC\n9NzuPh633JfP2yunNR8z7Qw6atQo7N69G++99x4qKyvx/vvvIxQKYejQoZg9ezaKi4vD3u0gIiKi\nLkjARycRk0EtFguKi4tRXFxs1HqIiIhMKxFzNKJWnRAREZFBGGgQERGRXkJmytHoCbW1tSgvL0dZ\nWRlKS0vhdDrxxBNPoLm5GUOGDMGmTZvgcDiiT2QAf6uaAHjy42qpLxQQkhKT1PXa7PLaLVY1EdMq\n9Nk0W5VLfcL3x1htciRrFZIorZpEUbuQ9GkXtjh3JMu/pOL7pPm0yaBSn/qe1BR5PjF5tZ8myTVV\nWFOakLyq3T5dPBbb2mRY8VicO8PRTxqX4chWzzv0BqlPTHoVk1xTNOcSk1LFXwar5g+9mOQ6QLju\n/VM0SbgBdSt4a5u8tbpF2Cbe4hXa/lZpXNAlJLK61STUYPNFeVyTmrDa1vCV0va55O3oAz41iTQU\nCL/pXlAYF/TKiafiHEGfL6Zxfo86zt8qjxMTWwM+eU3a5Fh17RESb4Vk0+AVSbgdJ+iKibbaNfna\nNAm6XdgWXn49/LhIibFCXnAHCbXdy9Ls6prEX62eTtC9clz3t+PXXjcyTtRA49e//jVeffVVNDc3\nIxQKKRt2HTp0KOL73G431q1bh0mTJimvbdy4EXPmzMHMmTNRUVGBPXv28JthiYiIvtZLRUW6imln\n0A0bNmDIkCGdmtjhcGDbtm3Ytm2b8trRo0fx1FNPAQCKi4vxxhtvMNAgIiL6WjABI42ogcbw4cMx\nfvz4zk9st0tfKw8ALpcLKSkpAIDs7Gw0NDR0el4iIqJElXhhRoRAo6qqCgAwevRoPPvss7j55pth\ns6nPocVHIrEStyu//AiGiIiI2iViKknYQGPr1q3S8bFjx5S2xWLpUqCRnp4Oj8eD1NRUNDQ0ICcn\np9NzEBFS4GBHAAAgAElEQVQRJaqQmR6d7NixAwBw5MgR3HLLLVLfgQMHunSyKVOm4ODBg5g1axb2\n79+PoqKiLs2jt+b6k2H7xMqSWPu0r1tsHY+zdmJu8dgiVFNEHhdhflvH47TvESttbA5NNY1QkRNp\nDqk6R6jisGiqU8Q7XnaHWGUjzydW3dg01STinDabtcPXr1iHuC26pqRdnF9cR0aE6pxUh9znENbh\nEB4vJtuz5HH2/kJb+IxJmiqeHKFyJ08dl50qVwUNzUxW2rkZcl+6cA2TxaojTXWSQ6riEa5t0CeN\nkypr2oTt4z1N0ji41QqaUJtcdYNgmKoTv3yukDhO6Au1yVU34vtC3lZhnHxecXt6v0vu87e2dbim\nK9YoVcKo7aCm2kes1hGrfbQVQz5hHT6XvAZttU5n1ydudd8+X8fb4mu3vvcJW/BrK3e0FT/h1xES\n2uHX5BOuYazb7Lf45TlqLsX2a9fbTHVH48yZMzh9+jQqKiqwfPlyJcpqa2vD+vXrMX369IgT19TU\noKKiAnV1dbDb7di3bx82bdqEpUuXYvv27SgoKMDMmTN79tMQERH1YQkYZ4QPNC5cuIC9e/eirq4O\nW7ZsUV63Wq2YN29e1InHjBmj3BURdfQaERERmeyOxrhx4zBu3DgUFRVFvXtBRERE3WeqHI3LCgoK\n8Oijj+Kzzz6DxWLBqFGj8NBDD6GgoMCI9REREZlG5zNu4l/UQGPZsmW477778MgjjyAUCuGTTz7B\n0qVL8fbbbxuxvl4RCpOEFq0vHkRK+Ozp+WI9V1fXFGvSbMwJrxGSZsPNd8W5wyS8arejF7eZ7+q5\npM8V5lpox4nnvSJptgsJtTZbhORaYZy2VD3W74US15HTP1vqyx+YrrSHZauJx9ok1wzhc2UKW/Cn\nJcmLSBPOJSbkpmvGZTiEPrvmcwW8HXyKKIS/MyxBv6ZPPbb4vUJbTlyUtqr3uqQ+bXJsLEI+4XNo\nk2vDJM2KW+Jrj4OtbqlP3u6948RYbV/AKySXtsrXWUqGdQvtS/K1EF06Iyce1/z2f8KOjScJeEMj\neqDRr18/3H333crxiBEjsGfPHl0XRUREZEaJmKMR9f8dEyZMwKFDh+B2u+FyuXDgwAFMnDgRoVAI\nwWAi3uQhIiKinhL1jsaWLVsQCFz5uGDr1q2wWCw4ceKELgsjIiIyG1Mmgx4/fhxWa4wPXImIiKjL\nEvE5QdRA4+LFi9i6dSucTieef/557N+/H2PHjsXgwYONWB91gZiwarHalONw7WjjgkKSmtXuUI5t\njlQEhCQ18Vhs21My4G9tUcaJx+HaAJCU3g8+V9MV7eTMbLQ1f6WMS+k3GK1NFwAAqQOGwNN4TukT\nj9MG5sHtrAcApA++Gq4Lp5Vx4nFG7jVoOf/FFW0AyLpqJC6dqQUA9PvGdUp7wDVj0PhFjTJu4Igb\n4fzsEwDAoJETlHZHfQ21H1/RBoDBo2/Bhb8cAQDkXD8ZX/7p8BVtAMi94TacP/4+AGDouOk491//\nCQC46qbpqPvkP5Vx+bdMx98+bj8efut0fHFE7RsxpRgnDx8CAIwsmoba//MeAOC6adPxl0PquOum\nTcOJ99r7xvz9dPzpwNftGcWo2X9QGXdDSTGO7zsYsQ0A42b9PY7tbT/++/vvwP5f/V7p++7iOdj+\ns70AgJVPzMH6534DAHhmxRz8oOI3yrgtq+7Cw+vbx/3bk3fiO0+3z/GrJ2fi/h+ruxi/81QJvv1U\n+/Hv1v497lj9LgDg0Pp/wNQfquf9w7MzcdOS9vlO/OQOjFnyH0rfnzYWY+KP26/90ZWTUbyl/dfx\n4OIbMfff/6yMe+v+/4Xy334OANg6Mx9PHKxrX/vtuVhb5VTGrb65H577r/aExiXXJ+Fn7W/Bgm8A\n/39DpjJuTpYHHwavBgBMxl9wIvMGpW904yf4n4E3AgCudcrt/86+URk36iv1eOT5Izg5tH3H58LT\nH+CLb9ymjMv//CDqR7RvaTD0z7+F84Z/BAAMOPILtH6rVBnn2LcVljseaj94eyPs316i9PneegaO\nOe3HbbueQfI9ywAArb+qQMq8Zco4z79XIPXepQCAlp+vR2Zpe9/Fnz2NAf+szvfl1qeRu7B93JkX\nfoyrH2pvf/4vT6HwsceUcSeffx6FS9r7Pn1sFa5/+F6lD79dh+yvE4e/8sZvUn8C3tCAJRTlPk15\neTkmT56M3/72t/jlL3+J/fv341e/+hVeffVV3RdXXV2NSd/dEn0gUS/Qu+qmKyJu/S5UrkSq3Inl\n9fa+7t/pjLjdfbJaaSJW9VgjVPiIVUHacV35ddBuaZ+elaK0Mweo6xswIEUalyOMyxAqYTI1W9WL\nW9f3S1W/dHJgmmb7+Ax1+/hsYRwA2ML8MkRKKhTfY9NUDIlby4tbzqdqKnCSxaojTTWeJRTm/+Xa\n18NU5Fh88tbvYhWOxa/2WTxyJQyEX2Pv53+Sutb+478obT0CjQV/eLtL33Quqq6uRnb+yB5aUWy+\nOlXb7XVHE/VvCpfLhQceeED5yvcZM2bA6+1CiRcRERFFFDL4xwhRH534fD74fD6lTt7pdMLj8UR5\nFxEREXVWIpa3Rg007r//fsydOxcXLlzA9773PRw/fhwrVqwwYm1ERESmkog5GlEDjVmzZmHChAk4\nfvw4LBYL1q5di5ycHCPWRkREZCrBBPz+1rCBxjvvvNPh6x9++CGCwSBmz56t26KI+oJYt6OP923r\n41VXkjf1TrwNl7xqS4ot8fSKBN0wCbVXbNWtqRATJfcbpLQzc/OVdkb/VGlcSpqaRJokJKgmp8r/\nDGQK44YKc+QPTJPGDRUTXh3ymqyaBNPLbNbwiadicmlaUoY0Lt3RT2mnJKvjktPl+cS7AYOnjJX6\nxmQ9q7Tfb5C3TI8nprqjcfhwewlXY2Mjamtr8c1vfhPBYBDHjh3DjTfeyECDiIioh5kqR+OZZ54B\nACxZsgS///3vkZbWHs22tLRg5cqVxqyOiIjIROLtjkZtbS3Ky8tRVlaG0tJSrFu3Dp9++inS09u/\n9HD+/PmYOnVqxDmi5micPn1aCTIAICMjA/X19d1bOREREV0hnnI03G431q1bh0mTJkmvPf3007ju\nuutinidqoFFYWIj77rsPY8e2P+/64x//iPz8/CjvIiIior7M4XBg27Zt2LZtm/Kay+Xq9DxRA431\n69fj8OHDqK2tRSgUwsKFCzFlypSYJu+JWy5EZE5dSaLtrcRbf2v0MXoRt9y/KGyDD3Rtx9dIwiWl\niru4Au3b/V+WOmCI0k5Ol5M8k4VdTlMz5YTa9Cx1N9Q8IRE1f2C6NG5YtnrugJDgMDRT3q11xHh1\nHe/v+xzxKp4endjtdmWzzstcLhc2b96M5uZm5ObmYtWqVejfv3/keWI52eTJkzF58uROLbCnbrkQ\nEVHn9XQFjpFBnBhkRCIGGYkiGE+RRgfuvfdeFBQUYMSIEXjllVewefNmrF69OuJ7dPta1su3XMQ9\nN7pyy4WIiMgsAkFjfzprxowZGDFiBACguLgYtbW1Ud+jW6Bht9uRkiLfurp8y6W0tBRLlizBxYsX\n9To9ERFRnxMMhQz96azy8nKcOXMGAHD06FFce+21Ud8T06OTntKVWy5ERERmETD40Umkuw01NTWo\nqKhAXV0d7HY79u3bh9LSUjz22GNITk5Geno61q9fH/UchgYaM2bMUNrFxcVYs2aNkacnIiKKa/GU\nozFmzBjs2LHjitdLSko6NY+hgUZ5eTl+9KMf4aqrror5lgsREXVepOTNnkjsFOcI+r1SX5OrSW3/\n7US3zyWyp6iVK/2/IRcWpA3MVdrDCrOlvl8VCWPjuOqkK3kT8U63QKOnbrkQERGZRTzd0egpugUa\nPXXLhYiIyCyMztEwgqGPToiIiCg8U32pGhERERkrkICRBgMNIiLqM/ytLUq7ofZjqU/cDdXTOEHq\ny15ws9K2WX6ntANx9u96IuZo6LZhFxERERHvaBAREcWJeLvD0hMYaBAREcWJRHx0wkCDiIgoTjAZ\nlIiIiHTDOxpERERxStwW3e2sk/psI+9R2g6rRWl74iwpIs6W0yMYaBAREcUJ3tEgIiIi3QSZo0FE\nRER64aMTIiIi0g0fnRAREfUB/laXdBzof5XSju9kUIPXY4k+pLsYaBAREcUJw3M0bNGHdBcDDSIi\nojgRZzdYegS/VI2IiIh0wzsaREREcYLJoERERKQbw5NBDRD3gUbVzxb39hKIiKiPq/nvz5T2PVVv\n9uJKIuOXqhls/Pjxvb0EIiIiwzDQICIiIt0w0CAiIiLdMNAgIiIi3TDQICIiIt0w0CAiIiLdMNAg\nIiIi3TDQICIiIt0w0CAiIiLdJGKgwS9VIyIiIt3wjgYREVGc8CfgHQ0GGkRERHEiER+dxHWgUV1d\n3dtLICIiiklPfD8XA41eMOm7W3SZd+K9/690/E/rFint0x6fLuckIqLE9M9H3+qRefg18URERKQb\n3tEgIiIi3TDQICIiIt0w0EggR3+1A9+8617l2GG19OJqiIiIgEAwaPAZ9f+3z7SBhhhkEBERxQPj\n72gw0CAiIjINPjohIiIi3XBnUCIiItJNIt7R4JeqERERUYdqa2sxffp0vPHGGwAAp9OJ+fPn4557\n7sEjjzwCr9cbdQ7TBhpfnW2WfgYn25Qfh9Ui/dgsUH6IiIj0EgiGDP2JxO12Y926dZg0aZLy2saN\nGzFnzhzs2rULw4YNw549e6J+JtMGGkRERPEmngINh8OBbdu2IScnR3nt6NGjmDZtGgCguLgYlZWV\nUT8TczSIiIjiRDzlaNjtdtjtcpjgcrmQkpICAMjOzkZDQ0P0eXRZHREREXVaPAUaHUlKSlLaoVAI\nFkv0nAIGGkRERHEiFOeBRnp6OjweD1JTU9HQ0CA9VgnHtIHG6Y9+h6yrRirHI8YPUdrej89KY5t8\n6pawX7jVDNtAfP9+ICKiPiYY54HGlClTcPDgQcyaNQv79+9HUVFR1PeYNtAQgwwiIqJ4EArFT6BR\nU1ODiooK1NXVwW63Y9++fdi0aROWLl2K7du3o6CgADNnzow6j2kDDSIiongTT49OxowZgx07dlzx\nekevRcJAg4iIKE7E+6OTrmCgQUREFCdCRn9LvAEYaBAREcWJeMrR6CmGBho1NTUoLy9Hfn4+AGDk\nyJFYtWqVkUtQuBvqpePrH5mmtNMGVkl9SenJSrvuyBml/d9/uySNa/GroagroLYDmt84XuHWmDfC\nbTJWtRARmQsfnXST2+1GSUkJVqxYYeRpiYiI+oR4SgbtKYZ+14nL5TLydERERNTLDL+jUV1djQcf\nfBA+nw+LFy+WvhWOiIjIzBLxjoahgcbo0aOxaNEilJSU4NSpUygrK8O+ffvgcDiMXAYREVFcCjIZ\ntHsKCwtRWFgIAMjPz8egQYNw/vx5XH311UYuAwDgb22RjpP+n5VKu0AzNhRUEztzxqk7it5ok588\nBbw+pd34l1NK2/VlszSu7VKb0v7qfxqlvr995VH7vAGlrU0a9YU6TiiNlEAqJqUy0ZSIKP7wjkY3\n7d69G01NTSgrK4PT6YTT6URubq6RS4hrYpBBRETmw0Cjm6ZPn45ly5bh3Xffhd/vx5o1a/jYhIiI\n6Gssb+2mzMxMvPzyy0aekoiIqM8w5YZdjY2NqKqqQl1dHQAgLy8PkyZNQnZ2tu6LIyIiMhNTbUHu\ncrlQUVGBQ4cO4aabbsKQIUMAACdOnEBFRQWmTp2KH/zgB0hPTzdssURERInMVI9OFi5ciDlz5mDN\nmjWw2WxSXyAQwO7du7Fw4UL84he/0H2RRvjdGbXCY1ZxadhxqW1CBUkgIPVZ/K1KO+vmL5V2qNUt\njQt61I3LQt5WhOX3Ks2Wv/5N6nKd+0od5lKTSCNVuLQ2que6qNk+/VyrXz2XsJW6ttqlNcbt0yNv\nrc7qFyKijpgqGXTjxo0YNmxYh302mw1z587FrbfeqtvCiIiIzMZUgcblIKOqqgo7d+5Ec3OzlKTy\n2muvIS8vT/8VEhERmYQpN+xavXo1Fi1axP0uiIiIdGaqOxqXXXPNNZg7d64RayEiIqIEEzXQuOee\ne7By5UqMHTsWdrs6/K677tJ1YUb7l7ePK23ft8dIfTaL2h6QOkRppyXJSbJZGer16TfwBqXtECfQ\nHGv7rAF1G3OLX03k7D9JTvIcICSeWnxtwnvk5NKgS0j69PsQi2Cr/C27/vovlLa38aLSbruoSTxt\nVLd197eqa/K55DV5nOr8l86oc2gTVJt84bdgb/GLCaXhk0vl7dnVtraCLNwckZJaIwn08O1PJs0S\nmYMp72j89Kc/RVpaGtrahH/MLJaECzTigRhkxAsxyIgXYpBBRJRITFXeell6ejpef/11I9ZCRERk\naom4M6g12oDbbrsNR48ehdfrRTAYVH6IiIioZ4WCIUN/jBD1jsZLL70Ej6d9QyiLxYJQKASLxYIT\nJ07ovjgiIiIzMdWjE5/Ph6SkJFRXV8Nq7fjGx+UxRERE1H2hYCD6oD4mbKAxf/58PPnkkxg+fHiH\n/SdPnsRTTz2Fn//857otzki1/+eQ0l7b7A07LiVdDaxsNjkAc6Sql7N/vxSlnZkiX+aMFHWOjGS1\nL80hV7H0S1PHDUx3SH1pSVnqHA5hjiR5TWn91TlT7Gqf1SJXu9itFmGc2k7+X/J84kdO0VTMZInz\nh6meAYCQJeoTu/Y53I3qeVsuKO1gk1MaF2hWK2GCl+Q+serG16xuBe8Ttm0HgLaLasVM89/OK+2L\np5qkceI27iGhFMQnbOEOAEGv+pdFUFMyEgizxbu2wkWsXJHHScPCzqGtfIm1ciVSxUy4h6aR3iOe\ntzNr6unKnfDnMeQ0vS7W62nT/N0QzyJ9pr70OUSmCjRWrFiBxx9/HHl5eZg8ebKyC2h9fT0qKytR\nX1+PjRs3GrZQIiKiRGeqQGPUqFHYvXs33nvvPVRWVuL9999HKBTC0KFDMXv2bBQXF4d9pEJERESd\nFwqYKNAA2pM/i4uLUVxcbNR6iIiITMtUdzSIiIjIWAw0Oqm2thbl5eUoKytDaWkpnE4nnnjiCTQ3\nN2PIkCHYtGkTHA5H9IkM4Gk8p7RPHT0k9QWFX3ibXV2vxSonb1rtSUJbHWdNkj+jOIc4TjunLTlV\naSelpMnnErIy7UISqV2zLbpNSOwUk1ctVjlRSuwTk1q1czqEc/VPkyuOUoWk1FRhnMMW2yM2m2ZN\nYnJsv7RBQs8gZAkJtWlZ6ri0gfK5ku1Cn5AoKybGAkCS8BiwVUjW9Gn2jBmUpv56DU5TP++AZHk+\nq0dNZLW2ylu1W7weoa1uxx70yFu/B93q+0JCUmuw5aI0zn2mXm2f+0ppa7eIjyQUUD9nwKcmtopJ\nrYCc9Cr2BTTjAr5gh31+TdKs36Mee13yzrjimrQJtfI4YWv5UGzJsOESbbXjtH1ysm3H83VGuDm0\nHzfWBN2urKMnkoT1TuoNN39fTfiMJBEDDd2SLNxuN9atW4dJkyYpr23cuBFz5szBrl27MGzYMOzZ\ns0ev0/dJ2sAlHmgDl3ggBhlERBTfogYav/71r3HnnXfi9ttvx9SpU1FUVISpU6dGndjhcGDbtm3I\nyclRXjt69CimTZsGACguLkZlZWXXV05ERJRgQsGAoT9GiGln0A0bNmDIkCHRhsoT2+3St70CgMvl\nQkpK+/4S2dnZaGho6NScREREiSwRH51EDTSGDx+O8ePH98jJxF1EL29lTkRERO2CZgo0qqqqAACj\nR4/Gs88+i5tvvhk2m/q8Xsy9iFV6ejo8Hg9SU1PR0NAgPVYhIiIyO1Pd0di6dat0fOzYMaVtsVi6\nFGhMmTIFBw8exKxZs7B//34UFRV1eg4jtDV/FbbPF7ZH1hOJneIckeaL2GfruM8a43yRKmuu6Evq\nuCLninFhzn3luTquztFW8Uh9msoVsbpGvIOmHReOtjonSdgyXtyOPlmToGp3qOlPycnyHzOxIidZ\nuJ4Oe7Y0zmEfrLS1FTkSIV53DFXPe1V2qjRs+MB0pf0NYYt8AMgUK4aS1HM5NOdNtouVO2pfkmZz\ncnHbeYtf3bZdW4EjbjMfaPxS6gu1qlvGWxzyeqVxXmFb+DZPh68DQMjvE8apfQGPWxoXaFW/gkC7\nVb3P1drhOC1LmEorsZIGAPzCHIHWtg5fBwC/xyf0yZU7AWFPerFPWwkkr0/9tQtpt8gXq4mE6qFY\nK5AAwNemHse6LX6kSiBx2/3I48JXuJwWrmE8V6uYKtDYsWMHAODIkSO45ZZbpL4DBw5EnbimpgYV\nFRWoq6uD3W7Hvn37sGnTJixduhTbt29HQUEBZs6c2c3lExERJQ5T7Qx65swZnD59GhUVFVi+fDlC\nX0eRbW1tWL9+PaZPnx5x4jFjxijBiqij14iIiMhkdzQuXLiAvXv3oq6uDlu2bFFet1qtmDdvniGL\nIyIiMhNTBRrjxo3DuHHjUFRUFPXuBREREXWfqQKNywoKCvDoo4/is88+g8ViwahRo/DQQw+hoKDA\niPX1WT3xm6Uv/4bTc5fTrs4d6X3hrnWk92iTUkV2h5qIKSbQdlWsycBi8q92e3txTTaHnCgqfpZI\nybXi1vdin1WzpbtVSLbTJtTKa1f77Enydeo3aKjSHjmsn9LOHyhvx99P2Ao/U9gWPsMhX7O0JFvU\nNgCkCFvVpyfJnytTmDNLSIbVpiD6hCxHsc+uTS4Wr2FI/T0oJtO2H6tJqBZfbH2WoJw0GlZITlBF\nQE2aFLfFF5NzATnZNihskQ8AwWZ1m3x/S4twKvlcYnKsX0iGFZNuAcAvHIsJutpxHqe8jb/otXdq\nw/bFE+01SgRRA41ly5bhvvvuwyOPPIJQKIRPPvkES5cuxdtvv23E+oiIiEyjL/8HM5yogUa/fv1w\n9913K8cjRozgd5QQERHpIBEDjajfdTJhwgQcOnQIbrcbLpcLBw4cwMSJExEKhRBMwFs8RERE1HOi\n3tHYsmULAh3U9W7duhUWiwUnTpzQZWFERERmY6otyC87fvw4rFbdvk2eiIiIvhZPG3bV1NSgvLwc\n+fn5AICRI0di1apVnZ4naqBx8eJFbN26FU6nE88//zz279+PsWPHYvDgwdHeSiam53PGeHmGGfB6\nwvb5XE2GraMrVTh6VgX11PziHNVhqmLax1k7fM8V42zht8UPd17tdvninEnpaiVM/5z+0rj+g9XK\nmKzMZKWdpqmEcQjVOhnCNvb9UzUVOEJlzYA0+XMNSFXXoa2gEfmECo+gsI23VbMdt1SRkym0s7Vr\nF7aqt8lzpAufyy70af/LKn4tgFjgk66tYvIK29EH1O3ZLW1ylYnVo/6501bdXHfwfyvtE81y5U48\niZe/3wDA7XajpKQEK1as6NY8UW9VrFy5EgUFBTh37pzy2vLly7t1UiIiIrpSKBgw9CcSlyt8uXBn\nRA00XC4XHnjgAdjt7Tc/ZsyYAa83/JcJERERUdfEU6DhdrtRXV2NBx98EKWlpcq3undW1EcnPp8P\nPp9PucXldDrh8YS/ZUxERERdE0+PTkaPHo1FixahpKQEp06dQllZGfbt2weHI/xmhR2JGmjcf//9\nmDt3Li5cuIDvfe97OH78eLef1xAREZHM4XCg6meLDT9nOIWFhSgsLAQA5OfnY9CgQTh//jyuvvrq\nTp0jaqAxa9YsTJgwAcePH4fFYsHatWuRk5PTqZMQkX668j+gePpfUyz8rdHHGCFcEum5CONiTYwN\nt5W8lk2T5CpucS8mq6b0kxP20wbmKu3UDDVBNTlV/mfAJmwzbxeSVx0p8rhU4TgnK0Xqu2qAusV9\nqiPCZxG2YBcTY9Md8rnErd9ThWTVzOQsaVyGY4A6tyZB9fa5o5V2zWvH1HGW8FvkG+2GG27o7SVI\ndu/ejaamJpSVlcHpdMLpdCI3Nzf6GzXCBhrvvPNOh69/+OGHCAaDmD17dqdPRkRERH3D9OnTsWzZ\nMrz77rvw+/1Ys2ZNpx+bABECjcOHDwMAGhsbUVtbi29+85sIBoM4duwYbrzxRgYaRERECSwzMxMv\nv/xyt+cJG2g888wzAIAlS5bg97//PdLS2uvCW1pasHLlym6fmIiIiBJf1PLW06dPK0EGAGRkZKC+\nvl7XRREREVFiiJoMWlhYiPvuuw9jx44FAPzxj39UtiMlIiIiiiRqoLF+/XocPnwYtbW1CIVCWLhw\nIaZMmRLT5LW1tSgvL0dZWRlKS0uxbt06fPrpp0hPTwcAzJ8/H1OnTu3WByAiMkqs1Tp6VvX4YhzX\ncv4L6Thc9UtXt4uXtmpP0lTCCH32lHSlnTZwmDQudcAgpZ0ibK2eminP50hW/6lKSlbnzhWqWwCg\nnzBHIBiS+l5eUKb2vfp9pW2Ln6KThBU10ACAyZMnY/LkyZ2a2O12Y926dZg0aZL02tNPP43rrruu\nc6skIiKiPkm3r2V1OBzYtm2btOdGT+2bTkRERH1DTHc0ujSx3a58P8plLpcLmzdvRnNzM3Jzc7Fq\n1Sr0798/zAxERETU1+l2R6Mj9957L5YsWYI33ngDo0aNwubNm408PRERERlMtzsaHZkxY4bSLi4u\nxpo1a4w8PRGRaYVLUO2JxNWAV/6iTTFR1OtqUtpuZ33YcZFetyerSZ/iNuuOzAHSuEHDR4VdY+ud\nJWH7SF+G3tEoLy/HmTNnAABHjx7Ftddea+TpiYiIyGC63dGoqalBRUUF6urqYLfbsW/fPpSWluKx\nxx5DcnIy0tPTsX79er1OT0RERHFAt0BjzJgx2LFjxxWvl5Tw9hUREZFZGProhIiIiMyFgQYRERHp\nxtCqEyIiSnzd3apdW3Xi9Xs7HNfW/JV83NSgtK32JKmv8m83x7Qm6nm8o0FERES6YaBBREREumGg\nQURERLphoEFERES6YTIoERHFFW2SaLityrXjxORQ7Xv+83/URFGH1dLdJVIn8I4GERER6YaBBhER\nEV4kUbEAAAedSURBVOmGgQYRERHphoEGERER6YaBBhEREemGVSdERBTXYt3SPJL/On1RaY+1qVUn\nnkCo23NTZLyjQURERLphoEFERES6YaBBREREumGgQURERLphoEFERES6sYRCobhNua2uru7tJRAR\nEcVk/Pjxvb2EuBTXgQYRERH1bXx0QkRERLphoEFERES6YaBBREREumGgQURERLphoEFERES64Zeq\nxaHa2lqUl5ejrKwMpaWlcDqdeOKJJ9Dc3IwhQ4Zg06ZNcDgcvb1MQzz33HP46KOP4PP5sGDBAkyc\nONGU18Lj8WD58uVwOp1wu91YvHgxxo4da8prcVlrayvuuOMOLF68GEVFRaa9FjU1NSgvL0d+fj4A\nYOTIkSgvLzft9fjNb36D7du3IxQK4dFHH8UNN9xg2msRL3hHI8643W6sW7cOkyZNUl7buHEj5syZ\ng127dmHYsGHYs2dPL67QOB9//DFOnDiBnTt34rXXXsP69etNey3ee+89jBkzBm+88QZefPFFbNy4\n0bTX4rJ//dd/Rf/+/QGY988I0P53RklJCXbs2IEdO3Zg1apVpr0eLpcL27dvxy9/+Uu8/PLLOHDg\ngGmvRTxhoBFnHA4Htm3bhpycHOW1o0ePYtq0aQCA4uJiVFZW9tbyDDVu3Di88MILAICsrCz4fD4c\nOXLElNfijjvuwIIFCwAA586dQ25urml/XwDAyZMncfLkSUydOhWAef+MAO3/uGqZ9XpUVlaiqKgI\nycnJyM3NxY9//GPTXot4wkAjztjtdqSkpEivuVwu5bXs7Gw0NDT0xtIMZ7fbkZ6eDgB48803UVRU\nBI/HY8prcdndd9+NpUuXYtWqVab9fQG038FYvny5cmzma+F2u1FdXY0HH3wQpaWlqKqqMu31OHv2\nLDweDx566CHcf//9pr4W8YQ5Gn1AUlKS0g6FQrBYLL24GuMdOHAAu3btwvbt2/HBBx8or5vxWrz5\n5pv405/+hMcffxw2m0153UzX4p133sFNN92Eq666SnnNzH9GRo8ejUWLFqGkpASnTp1CWVkZxA2f\nzXQ9vF4vzpw5g5/85Cc4ffo0ysrKTPvnJJ4w0OgD0tPT4fF4kJqaioaGBumxSqL74IMPsHXrVrz6\n6qvIysoy7bU4fvw4Bg4ciLy8PFx//fUIBoNITU015bU4dOgQzpw5g/379+PcuXNwOBxITk425bUA\ngMLCQhQWFgIA8vPzMWjQIHz55ZemvB6DBw/G2LFjYbPZcM011yAjIwNWq9WU1yKe8NFJHzBlyhQc\nPHgQALB//34UFRX18oqM0dzcjA0bNuCVV17BgAEDAJj3Wnz66ad4/fXXAQANDQ1wuVy4/fbbTXkt\nXnjhBbz11lvYtWsX7r77bpSXl5v2WgDA7t278W//9m8AAKfTCafTiblz55ryetx66604cuQIQqEQ\nnE6nqf+cxBN+qVqcqampQUVFBerq6mC325Gbm4tNmzZh6dKlcLvdKCgowIYNG2C3J/7NqJ07d+LF\nF19EQUGB8tqGDRuwfPly010Lr9eLH/7whzh79iy8Xi8WL16M66+/HkuWLDHdtRC9+OKLGDZsGL71\nrW+Z9lo0Nzdj2bJluHTpEvx+PxYvXozrrrvOtNdj586d+O1vfwuXy4XFixfjhhtuMO21iBcMNIiI\niEg3fHRCREREumGgQURERLphoEFERES6YaBBREREumGgQURERLphoEGUwD788EM89NBDOH36NG67\n7bYuzfH9738f77//fg+vjIjMgoEGUYJyuVx48sknsXbt2m5tu/zUU09h3bp1aGlp6cHVEZFZMNAg\n6oO2b9+OVatWAQA+//xz/MM//MMVgcBbb72F2267DdnZ2dLr586dw6xZs/CXv/wFy5YtQ0VFBR55\n5BFMmzYNv/vd77BkyRL84z/+I5588kkAQL9+/VBcXIy33nrLkM9GRImFgQZRH/TP//zP+Otf/4rq\n6mqsXbsWa9euRUZGhjSmsrISkydPll5raWnBww8/jCeffBKjR4+GzWaD0+nE5s2b8fDDD2Pt2rVY\ns2YNdu3ahbfffhtutxtA+9bO4hfaERHFioEGUR9ktVrx9NNP4/vf/z6uvfZaTJw48Yox9fX1GDp0\nqHIcCATw8MMPY9asWbjpppuU18eNGwcAyM3NxfDhw5GVlYXU1FRkZWWhubkZADB06FDU1dXp/KmI\nKBEx0CDqo5qampCWlob6+voO+7XfLtDU1IQxY8Zg165dyp0KANL3Pmi/A4LfUEBE3cVAg6gPamtr\nw+rVq/Hyyy/D4XDgnXfeuWLM0KFDcfbsWeU4OzsbS5YswfTp0/HjH/+4U+c7e/Ys8vLyur1uIjIf\nBhpEfdBPfvITzJgxAwUFBVixYgVefPFFnDt3ThrzrW99C4cPH77ivQ8//DBOnjyJvXv3xny+Dz/8\nEFOmTOn2uonIfPjtrUQJqqWlBXfddRd27dp1ReVJZzQ1NWHu3LnYvXv3FQmnRETRMNAgSmCVlZX4\n5S9/iZdeeqnLe2l8//vfx+zZs7u84RcRmRsDDSIiItINczSIiIhINww0iIiISDcMNIiIiEg3DDSI\niIhINww0iIiISDcMNIiIiEg3/xchhVzsxIDVygAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7ff18f00cd30>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"path = 'lock_exchange_kh100/prog_{}.nc'\n",
"cmap = cm.RdBu_r\n",
"\n",
"fig = plt.figure()\n",
"\n",
"ax = plt.subplot(311)\n",
"d = Dataset(path.format('plm'), 'r')\n",
"pcol = plt.pcolormesh(d.variables['xh'][:],\n",
" d.variables['zl'][:],\n",
" d.variables['temp'][-1,:,0,:],\n",
" linewidth=0, rasterized=True, cmap=cmap)\n",
"plt.axis('tight')\n",
"ax.invert_yaxis()\n",
"ax.xaxis.set_ticklabels([])\n",
"plt.ylabel('depth (m)')\n",
"d.close()\n",
"\n",
"ax = plt.subplot(312)\n",
"d = Dataset(path.format('ppmh3'), 'r')\n",
"pcol = plt.pcolormesh(d.variables['xh'][:],\n",
" d.variables['zl'][:],\n",
" d.variables['temp'][-1,:,0,:],\n",
" linewidth=0, rasterized=True, cmap=cmap)\n",
"plt.axis('tight')\n",
"ax.invert_yaxis()\n",
"ax.xaxis.set_ticklabels([])\n",
"plt.ylabel('depth (m)')\n",
"d.close()\n",
"\n",
"ax = plt.subplot(313)\n",
"d = Dataset(path.format('ppm'), 'r')\n",
"pcol = plt.pcolormesh(d.variables['xh'][:],\n",
" d.variables['zl'][:],\n",
" d.variables['temp'][-1,:,0,:],\n",
" linewidth=0, rasterized=True, cmap=cmap)\n",
"plt.axis('tight')\n",
"ax.invert_yaxis()\n",
"plt.ylabel('depth (m)')\n",
"plt.xlabel('x (km)')\n",
"d.close()\n",
"\n",
"plt.tight_layout()\n",
"\n",
"fig.subplots_adjust(right=0.82)\n",
"cbar_ax = fig.add_axes([0.85, 0.2, 0.05, 0.7])\n",
"fig.colorbar(pcol, cax=cbar_ax)"
]
}
],
"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.0"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment