diff --git a/notebooks/examp_plotting_split_axes_vertically.ipynb b/notebooks/examp_plotting_split_axes_vertically.ipynb new file mode 100644 index 0000000000000000000000000000000000000000..3200ca9b5b8c25d1f279704b9ff7e1b690109a0e --- /dev/null +++ b/notebooks/examp_plotting_split_axes_vertically.ipynb @@ -0,0 +1,211 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "id": "7f0d3233", + "metadata": {}, + "outputs": [], + "source": [ + "# Jupyter Notebook with widget matplotlib plots\n", + "#%matplotlib notebook\n", + "# Jupyter Lab with widget matplotlib plots\n", + "#%matplotlib widget \n", + "# with Jupyter and Jupyter Lab but without widget matplotlib plots\n", + "%matplotlib inline\n", + "%load_ext autoreload\n", + "%autoreload 2" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "id": "3f7a1734", + "metadata": {}, + "outputs": [], + "source": [ + "import numpy as np\n", + "import matplotlib.pyplot as plt\n", + "from netCDF4 import Dataset\n", + "import sys\n", + "import cartopy\n", + "import pyicon as pyic\n", + "import matplotlib" + ] + }, + { + "cell_type": "markdown", + "id": "c6f1dd69", + "metadata": {}, + "source": [ + "Generate some example data" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "id": "9d11da1f", + "metadata": {}, + "outputs": [], + "source": [ + "depth = np.linspace(0,6000,100)\n", + "lat = np.linspace(-90,90,180)\n", + "\n", + "L, D = np.meshgrid(lat, depth)\n", + "data = np.cos(L*np.pi/180.)*np.cos(D/6000)" + ] + }, + { + "cell_type": "markdown", + "id": "f4a595a4", + "metadata": {}, + "source": [ + "Example where pyicon is only used for splitting the axes:" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "id": "a34c9ca0", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdoAAAEWCAYAAADBzlZgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAmFklEQVR4nO3df5xcdX3v8debRJJgSPkRAyGJTYQIBa4oRKSiVQyWEJGoD2lji6YWm2KBorZKqPdWe1t6sXqtehG5KVCxUmlEKKlG+RGliJUf4TchIAvRsCQQAkXAkkDCp3+c75LJZnZ2dmfOnHNm3s/HYx47c35+dnb2vM/5njPfo4jAzMzM8rFL0QWYmZl1MwetmZlZjhy0ZmZmOXLQmpmZ5chBa2ZmliMHrZmZWY4ctGYlIenVkp6TNKbBNCHpgE7WZWatcdCaFUjSzyUdCxAR6yJiYkRsS+Oul/SRYis0s1Y5aM3MzHLkoDUriKR/Al4N/FtqMv5UahoeK+kc4K3AeWnceXXmHyfpC5LWSXpc0gWSJnT69zCzxhy0ZgWJiA8C64B3R8REYFnNuE8DPwZOT83Jp9dZxOeA1wKvBw4ApgF/mXfdZjYyDlqzCpIk4I+Aj0fEUxHxLPC3wMJiKzOzwcYWXYCZjcqrgN2A27LMBUDAkFcsm1kxHLRmxWp0+6xG4zYBzwOHRMSj7S3JzNrJTcdmxXoceM1Ix0XES8A/AH8vaQqApGmSjsulSjMbNQetWbH+D/A/JT0NvH/QuC8D75f0n5K+Umfes4A+4CZJzwDXAQfmWayZjZx843czM7P8+IjWzMwsR10ftJLmSXpAUp+kJUXXY2ZmvaWrm45T5+w/A94J9AO3Ah+IiPsKLczMzHpGtx/RHgn0RcTDEfECcBmwoOCazMysh3T792inAY/UvO4H3lQ7gaTFwGKAXcbsesSE3ae0tYBtOfY8+9L4/FsjJo7bkvs66pnyimcKWa9Vy8YXJxWy3ue2jMt9Hbts1vATjcKY59u/zF893b8pIl7VyjKOO+aV8eRT21qq47a7t1wdEfNaWkgOuj1o631Sd0iniFgKLAWYuOeMeP0xZ7a1gCcPya+jns0Hbc5t2QPecsBDua+jntP3XVnIeq1azntsbiHrvbFv/9zXMf7+8bksd+/VrYVZPT+58pO/aHUZTz61jVuufnVLyxgz9cHJrdaRh25vOu4HZtS8ng6sL6gWMzPrQd0etLcCsyXNkrQrWYfrywuuyczMekhXNx1HxFZJpwNXk3W2fnFErC64LDMz6yFdHbQAEbECWFF0HWZm1pu6venYzMysUA5aMzOzHDlozczMcuSgNTOzniDpYkkbJd3bYJq3S7pT0mpJ/96O9TpozcysV3wdGLLnKEl7AOcDJ0bEIcBJ7Vipg9bMzHpCRNwAPNVgkt8DroiIdWn6je1Yr4PWzMy6xWRJq2oei0c4/2uBPSVdL+k2SR9qR1Fd/z1aMzPrGZsiYk4L848FjgDmAhOAn0q6KSJ+1kpRDlozM7NMP1lY/wr4laQbgMPI7ms+am46NjMzy1wFvFXSWEm7kd1WdU2rC/URrZmZ9QRJ3wLeTnYutx/4DPAKgIi4ICLWSPoBcDfwEnBhRAz5VaBmOWjNzKwnRMQHmpjm88Dn27leNx2bmZnlyEFrZmaWIwetmZlZjhy0ZmZmOXLQmpmZ5chBa2ZmliMHrZmZWY4ctGZmZjly0JqZmeXIQWtmZpYjB62ZmVmOHLRmZmY5ctCamZnlyEFrZmaWIwetmZn1BEkXS9ooqe49ZiX9vqS70+M/JB3WjvU6aM3MrFd8HZjXYPxa4G0R8Trgr4Gl7Vipb/xuZmY9ISJukDSzwfj/qHl5EzC9Het10JqZWbeYLGlVzeulETHao9JTgO+3oSYHrZmZdY1NETGn1YVIOoYsaN/SekkOWjMzs5dJeh1wIXB8RDzZjmVW7mIoSWMk3SHpu+n1XpKulfRg+rlnzbRnS+qT9ICk44qr2szMyk7Sq4ErgA9GxM/atdzKBS1wJrCm5vUSYGVEzAZWptdIOhhYCBxCdpXZ+ZLGdLhWMzMrCUnfAn4KHCipX9Ipkk6VdGqa5C+Bvcny4s5B53tHrVJNx5KmA+8CzgE+kQYvAN6enl8CXA+clYZfFhFbgLWS+oAjyd5kMzPrMRHxgWHGfwT4SLvXW7Uj2i8BnwJeqhm2T0RsAEg/p6Th04BHaqbrT8N2IGmxpFWSVm3d8lwuRZuZWe+qTNBKOgHYGBG3NTtLnWGx04CIpRExJyLmjB03saUazczMBqtS0/HRwImS5gPjgUmSvgk8LmlqRGyQNBXYmKbvB2bUzD8dWN/Ris3MrOdV5og2Is6OiOkRMZPsIqcfRsTJwHJgUZpsEXBVer4cWChpnKRZwGzglg6XbWZmPa5KR7RDORdYJukUYB1wEkBErJa0DLgP2AqcFhHbiivTzMx6USWDNiKuJ7u6mPSF4rlDTHcO2RXKZtZlznus7r+9WelUMmjNzKy7/CqCW7e8WHQZuajMOVozM7MqctCamZnlyEFrZmaWIwetmZlZjhy0ZmZmOXLQmpmZ5chBa2ZmPUHSvHR/8j5JS+qM/zVJ/ybpLkmrJX24Het10JqZWddL9yP/KnA8cDDwgXTf8lqnAfdFxGFkt1/9v5J2bXXdDlozM+sFRwJ9EfFwRLwAXEZ23/JaAewuScBE4CmyLnxb4qA1M7Ne0Mw9ys8DfoPsTm/3AGdGxEu0yF0wmplZ4Ta+OKkN/VdfNFnSqpoBSyNiaXrezD3KjwPuBN4B7A9cK+nHEfFMK1U5aM3MrFtsiog5Q4xr5h7lHwbOjYgA+iStBQ6ixVusuunYzMx6wa3AbEmz0gVOC8nuW15rHelucJL2AQ4EHm51xT6iNTOzrhcRWyWdDlwNjAEuTvctPzWNvwD4a+Drku4ha2o+KyI2tbpuB62ZmfWEiFgBrBg07IKa5+uB3273et10bGZmliMHrZVS61cfmpmVg4M2Z3uv3lZ0CWZmViAHrZmZWY6GvBhK0t1NzP9ERLiNz8zMbAiNrjoeA8xvMF7s/B0kMzMzq9EoaP84In7RaGZJf9LmeszMzLrKkEEbETcON3Mz05iZmQ3nuS3juLFv/6LLyMWwF0NJOkHSHZKekvSMpGcltdTBspmZWa9opmeoLwHvA+5JHS2bmZlZk5r5es8jwL0OWTMzs5Fr5oj2U8AKSf8ObBkYGBFfzK0qMzOzLtFM0J4DPAeMB3bNtxwzM7Pu0kzQ7hURbb+bgZmZWS9o5hztdZIctGZmHTL+/vFFl9CVJM2T9ICkPklLGkz3RknbJL2/HettJmhPA34g6Xl/vcfMel23ftez20kaA3wVOB44GPiApIOHmO5zZDeIb4thgzYido+IXSJiQkRMSq8ntauAkZD0cUmrJd0r6VuSxkvaS9K1kh5MP/esmf7stOfygKTjiqjZzMxK4UigLyIejogXgMuABXWmOwP4DrCxXSseMmgl7TvczM1M0y6SpgF/CsyJiEPJ+mJeCCwBVkbEbGBlek3aU1kIHALMA85PeypmZtadJktaVfNYXDNuGtnXVQf0p2EvSznzXuCCdhbV6GKoFcDhw8zfzDTtNBaYIOlFYDdgPXA28PY0/hLgeuAssj2VyyJiC7BWUh/ZHs1PO1ivmZk1YZfNase56U0RMWeIcaozbHD/EF8CzoqIbVK9yUenUdAeNsy5WAEdO1cbEY9K+gKwDngeuCYirpG0T0RsSNNskDQlzTINuKlmETvtvQCkPZ6BvZ7nfnLlJx9oU8mTgU0AXNmmJbZue01N+nk+ddQasqZv5r/uRkb8XnWAa9rBRUON8PvUvHbV9ettWEbe+oEZNa+nkx2s1ZoDXJZCdjIwX9LWiPjXVlbc6KYCpWpmTedeFwCzgKeBb0s6udEsdYbt1LtVRCwFlrajxh1WLq1qsGdVCNfUvDLW5Zqa45qaV9a6cnIrMFvSLOBRslOLv1c7QUTMGngu6evAd1sNWWjuquOyOBZYGxFPRMSLwBXAm4HHJU0FSD8HTmA3s/diZmY9ICK2AqeTXU28BlgWEaslnSrp1DzX3UyHFWWxDjhK0m5kTcdzgVXAr4BFwLnp51Vp+uXAP0v6IrAfMBu4pdNFm5lZOUTECrJri2qH1b3wKSL+oF3rrUzQRsTNki4Hbge2AneQNflOBJZJOoUsjE9K06+WtAy4L01/WkRs62DJbW+ObgPX1Lwy1uWamuOamlfWurqKmrkpT/pazD7UBHNErMuxLjMz6yET9psRM0/5REvLuP9vPnFbGc85D3tEK+kM4DPA48BLaXAAr8uxLjMzs67QTNPxmcCBEfFk3sWYmZl1m2Zv/P7LvAtphrtgNDOzqhnyiFbSQGP5w8D1kr5HgTd+r+mC8eCIeD5d6LSQrHPolRFxbrobwxLgrEFdMO5Hdhei13b4gigzM+txjZqOd08/16XHrmy/8fvwV1Dlw10wmpl1oTHPw96ru/M4qFHPUH8FIOmkiPh27ThJJ+VdWJ16cu+CcZcxux4xYfcpgydp2bYJbV8kAC+Nz39/Z+K4LcNPlJMpr/DdGG1oG18s5CZiPLdlXO7r2GVz+/rZHWzM8+1f5q+e7t8UEa9q/5K7QzMXQ50NfLuJYbnqRBeME/ecEa8/5szWix3kyUPy6c1y80Gbc1lurbcc8FDu6xjK6fuuLGzdVm7nPTa3sHV34n60ed74PY+jxp9c+clftH2hXaTROdrjgfnANElfqRk1iawDiE57uQvGVN8OXTCmo1l3wWhmZqXS6Krj9WRdHG4Gbqt5LAeKuIL35S4Yld1aYS5Zf5XLybpehJ27YFwoaVzqRNpdMFZMkUctZmbt0ugc7V3AXZL+mawZ9iCyptcH0t3pO6qCXTB2hRv79i+0+djMrOqaOUf7TuD/Aw+RBe4sSX8cEd/PtbI6IuIzZL1U1dpCdnRbb/pzgHPyrsvMzGwozQTtF4FjIqIPQNL+wPeAjget7Wj8/eM7ckGUmZmNXjM9Q20cCNnkYbZfcGRm1lE+d2+jJWle6imwL3VwNHi8JH0ljb9b0uHtWG8zR7SrJa0AlpGdoz0JuFXS+wAi4op2FGJmZpaXdBe6r5KdDu0ny7HlEXFfzWTHk104Oxt4E/C19LMlzRzRjie7c8/byHpgegLYC3g3cEKrBZg14qMXM2uTI4G+iHg4XdB7GVnfDLUWAN+IzE3AHulroy0Z9og2Ij7c6krMzLpBJzqrsJZMlrSq5vXS1CkRZD0DPlIzrp+dj1brTTMN2NBKUc3cj/a1ZIfP+0TEoZJeB5wYEX/TyorNzMwGjPmvF9j99kdbXcymBjd+b6a3wKZ6FBypZpqO/4Gsy8UXASLibrK74piZdZRPJVgLmuktMJceBZsJ2t0iYnCPSkV0wWgFKbq5zBtXM2uDW4HZkmZJ2pXsgHH5oGmWAx9KVx8fBfxy4KY1rWjmquNN6buzASDp/bTYXm1mZtZJEbFV0unA1cAY4OLUg+CpafwFwAqyPv77gP8C2nKNUjNBexpZV4cHSXoUWAs0umuOmZmNUp537ul1EbGCLExrh11Q8zzIMq+tmrnq+GHgWEmvBHaJiGfbXYSZ2XB8CsGqqtFt8j4xxHAAIuKLOdVktpPzHpvr+9NaoYq+VsGqq9ER7e7p54HAG9l+0vjdwA15FmVmZtYtGt0m768AJF0DHD7QZCzps8C3O1KdmRluNm6XvVf7TqFFaObrPa8Gau8/+wIwM5dqbMQ6deGEm83MzEanmaD9J+AWSZ+V9BngZuCSfMsy25mPasysipq56vgcSd8H3poGfTgi7si3LDMz6yXxwotsXddfdBm5aOZ7tETE7cDtOddiZraTMrRk+NSJtaKZpmOz0ijDRtfMbCRKF7SSZkj6kaQ1klZLOjMN30vStZIeTD/3rJnnbEl9kh6QdFzN8CMk3ZPGfUUDXwK2UfFevZnZyJUuaMluWPBnEfEbwFHAaZIOBpYAKyNiNrAyvSaNWwgcAswDzpc0Ji3ra8BiYHZ6zOvkLzLAl9SbjU6vtWC4+8XuVLqgjYgN6Zww6bu7a8huvLuA7Vc7XwK8Jz1fAFwWEVsiYi1ZZ9BHSpoKTIqIn6b+K79RM49VWK9tfM2s2koXtLUkzQTeQPaVon0GbleUfk5Jk00DHqmZrT8Nm5aeDx5uZhXgHSrrFqUNWkkTge8AH4uIZxpNWmdYNBg+eD2LJa2StGrrludGV2zBerG5yRth6xRfm2CtKmXQSnoFWcheGhFXpMGPp+Zg0s+NaXg/MKNm9unA+jR8ep3hO4iIpRExJyLmjB03sb2/SBfyRsfMulGjC27rTDtG0h2SvtvMsksXtOnK4IuANYPuELQcWJSeLwKuqhm+UNI4SbPILnq6JTUvPyvpqLTMD9XMY13AR7Xdy39bK0DdC26HcCbZ9UNNKV3QAkcDHwTeIenO9JgPnAu8U9KDwDvTayJiNbAMuA/4AXBaRAxc5vtR4EKyC6QeAr7f0d/EzMyqYqgLbncgaTrwLrJsaUpTPUN1UkTcSP3zqwB1d3Mj4hzgnDrDVwGHtq86M8tbmY5mfaqkciZLWlXzemlELG1y3h0uuJU0ZYjpvgR8iu23kh1W6YLWbCR8Q3jrFr14UWMONkXEnKFGSroO2LfOqE83s3BJJwAbI+I2SW9vtigHrY3YjX3785YDHiq6DOtCZTqate4TEccONU7S45KmpqPZ2gtuax0NnJhOZ44HJkn6ZkSc3Gi9ZTxHazYi3jibWRsMdcHtyyLi7IiYHhEzyXok/OFwIQsO2q7hZierurLtMHXb+Vl3BTusuhfcStpP0opWFuygta5Qto20mVVLRDwZEXMjYnb6+VQavj4i5teZ/vqIOKGZZTtobVS6bW/fiuUdJetmDlrrGt5Ym1kZOWitqzhsq6eMf7NOt9j4Govu5qDtkG68EMHNx9aqMoasWbs5aLuI94oz3nibWZk4aK0rOWzLr6x/I7fUWLs5aK0l3ijZaJQ1ZM3y4KC1ruWNuZmVgYPWuprDtnzK/DcpooXG11Z0PwetmXVMmUPWLC8O2i5TxN5x2c/TeuNuZkVy0FpPcNgWr+x/g7LvMLaiG7/HXyUOWmuLKmykyr6h72Z+762XOWitp3iD33l+z4fmC6F6g4PWzHJTlZCtQouM5UvSXpKulfRg+rnnENN9XNJqSfdK+pakYfeWHLQd1KnzJEXtJVdlY1WVjX/V+X22ilkCrIyI2cDK9HoHkqYBfwrMiYhDgTHAwuEW7KC1nuQQyJffX6ugBcAl6fklwHuGmG4sMEHSWGA3YP1wC3bQWs9yGOSjau9rVVpirCmTJa2qeSwewbz7RMQGgPRzyuAJIuJR4AvAOmAD8MuIuGa4BTtora2qttGqWiiUnd/P5vlCqFxsiog5NY+ltSMlXZfOrQ5+LGhm4em87QJgFrAf8EpJJw8339jR/CZWfuPvH8/mgzYXXUYlDITD6fuuLLiSaqtiyFZtx9BaExHHDjVO0uOSpkbEBklTgY11JjsWWBsRT6R5rgDeDHyz0Xp9RGttV9WNVxWDoiz83lkXWA4sSs8XAVfVmWYdcJSk3SQJmAusGW7BDlqzGg6Mkavqe1bVHcKRcq9QTTsXeKekB4F3ptdI2k/SCoCIuBm4HLgduIcsQ5fWX9x2bjq2XNzYtz9vOeChossYlfMem+tm5CZUNWDLwOdnyyciniQ7Qh08fD0wv+b1Z4DPjGTZPqI1q8Mh0pjfH7PmOWg7rJPNOEXvNVe9ae68x+Y6UOrohvek6p9NqxYHrdkwuiFY2sE7Hmaj46C1XHXLkUOvh0w3/e5FfyaLbmmyzvPFUF3O36dtr177zm03BaxZURy0lrsqX4E8lG4P3G4N2KKPZq03OWjNWtBtgdutAQsOWSuOg9Y6ohuPamtVPXC7OWDLpNPnZ91ZRTk4aHuAz9N2Tm1glT10eylcfTRrRXLQFmDv1dt48pAxRZfRcd1+VDtYGUO3l8LVrCwctGYdMDjgOhW8DtbyHM36az29y0HbI8rSfNxrR7VDGSoARxvADlSz8ipl0EoaA6wCHo2IEyTtBfwLMBP4OfA7EfGfadqzgVOAbcCfRsTVafgRwNeBCcAK4MyIiM7+JlaPw3ZoDsz2KsvRrPW2svYMdSY73uNvCbAyImYDK9NrJB0MLAQOAeYB56eQBvgasBiYnR7zOlO6mZVBmULWzca9rXRBK2k68C7gwprBC4BL0vNLgPfUDL8sIrZExFqgDzhS0lRgUkT8NB3FfqNmnp5Vpn/2Mm0EzbqRv9ozMpJOkrRa0kuS5gwxzQxJP5K0Jk17ZjPLLl3QAl8CPgW8VDNsn4jYAJB+TknDpwGP1EzXn4ZNS88HD9+JpMWSVklatXXLc235Baw5DlvLiz9bNgr3Au8DbmgwzVbgzyLiN4CjgNNSy2pDpQpaSScAGyPitmZnqTMsGgzfeWDE0oiYExFzxo6b2ORqW+e9TbN8lC1ky9SSZEOLiDUR8cAw02yIiNvT82fJTnHWPYirVaqgBY4GTpT0c+Ay4B2Svgk8npqDST83pun7gRk1808H1qfh0+sM73ll+6cv20bRzCpt8kALZXoszmtFkmYCbwBuHm7aUgVtRJwdEdMjYibZRU4/jIiTgeXAojTZIuCq9Hw5sFDSOEmzyC56uiU1Lz8r6ShJAj5UM4+VjMPW2sWfpZ63aaCFMj2W1o6UdJ2ke+s8FoxkJZImAt8BPhYRzww3fSm/3lPHucAySacA64CTACJitaRlwH1kbeenRcRAm+xH2f71nu+nh5WUv/JjrSpjyJatBanXRcSxrS5D0ivIQvbSiLiimXlKG7QRcT1wfXr+JFD3C4YRcQ5wTp3hq4BD86uwusrSeYVZu5QxZIvka0DykVpILwLWRMQXm52vVE3H1tu8sbTRKOvnxkez1SLpvZL6gd8EvidpoPOj/SStSJMdDXyQ7PqhO9Nj/nDLLu0RbS8o8uYCZT2qdROymRUhIq4ErqwzfD0wPz2/kfrfamnIR7RWOmU9QrHy8WfFqsBBa6XkDagNp8yfETcbWy0HbQ8r+8agzBtSK5Y/G1YlDlorNW9QbbCyfyaK3oH1Fcfl46DtcUVvFJpR9g2rdY4/C1ZFDtqCee+zOd7AWhU+A1XYcbXOc9BaZTYOVdjQWj78t7cqc9BapXiD23v8N7eqc9AaUJ2jWvCGt5dU6W9dhv8hn4oqJwdtCfifY+SqtAG20fHf2LqFg9ZeVoY98pHwhrg73di3f+X+tlX737HOctBapVVxo2xD89/SupGD1nZQ1T1zb6Crr6p/w7L8z/gUVHk5aK1rVHVD3euq3CpRlpC1cnPQlkSZ9karvPGo8ka7F/lvZWUh6SRJqyW9JGlOg+n2kHS5pPslrZH0m8Mt20FrXckb8HLrhh2iKu+QWl33Au8Dbhhmui8DP4iIg4DDgDXDLdhBa3V1w0akGzbm3ch/k/YrU4tYVUXEmoh4oNE0kiYBvwVclOZ5ISKeHm7ZDlobUjeELXjDXhbdtOPTLf8bXWiypFU1j8VtXv5rgCeAf5R0h6QLJb1yuJnGtrkIa8Heq7fx5CFjii6jKw1s4N9ywEMFV9J7uiVcBzhkS21TRDQ6v3odsG+dUZ+OiKuaWP5Y4HDgjIi4WdKXgSXA/xpuJrMhjb9/PJsP2lx0GW3jwO2sbgtZq7aIOLbFRfQD/RFxc3p9OVnQNuSgtWF1W9iCAzdv3RqwZTya9fnZzomIxyQ9IunAdD53LnDfcPM5aEvGzced5cBtn24N1wFlDFlrH0nvBf4f8Crge5LujIjjJO0HXBgR89OkZwCXStoVeBj48HDLdtBaU7rxqLZWbUg4dEem2wPWekNEXAlcWWf4emB+zes7gSHPA9fjoLWmdXvYDvBRbnN6KWB9NGutcNCaDcFHuTvrpXAdUOaQ9fnZanDQllCZz9P2ylHtYL0cur0Yrmbt5KC1EevVsB3Q7aHrYN2uzEezVh0OWhuVXg/bAfVCqWrh62Ctr+wh62bj6nDQllSZm48HOGzrK3P4OlSbU/aQtWpx0Jp1wHAB144gdoi2h0PW2s1Bay3xUW17OCTLoSoh62bjavHde0qsKv9MVdk4mTXiz7HlxUFrbeGNlFWZP7+WJwettY03VlZFVfvcVqWly7Zz0JZc1f6pqrbRst7mz6t1goPW2s4bL6uCKn5Oq7bjbRkHbQVU8Z+rihsx6x3+fFonlTJoJe0h6XJJ90taI+k3Je0l6VpJD6afe9ZMf7akPkkPSDquZvgRku5J474iScX8Rr3JGzMrI38urR5Jn0+Zc7ekKyXtMcR081LW9Ela0syySxm0wJeBH0TEQcBhwBpgCbAyImYDK9NrJB0MLAQOAeYB50sa6FLpa8BiYHZ6zOvkL9FOVTyqhWyj5g2blUWVP4tV3QZUyLXAoRHxOuBnwNmDJ0jZ8lXgeOBg4AMpgxoqXdBKmgT8FnARQES8EBFPAwuAS9JklwDvSc8XAJdFxJaIWAv0AUdKmgpMioifRkQA36iZxzqsyhs46w7+DFojEXFNRGxNL28CpteZ7EigLyIejogXgMvIMqihMvYM9RrgCeAfJR0G3AacCewTERsAImKDpClp+mlkb8qA/jTsxfR88PAdSFpMdtQL8NxPrvzkA236PSYDm9q0rMyVLS+h/TW1row1QTnrck3NcU3Na1ddv97qAp7lP6++Li6f3OJixktaVfN6aUQsHcVy/hD4lzrDpwGP1LzuB9403MLKGLRjgcOBMyLiZklfJjUTD6HeeddoMHzHAdkfYTR/iIYkrYqIOe1ebitcU/PKWJdrao5ral6Z6oqI3E/tSboO2LfOqE9HxFVpmk8DW4FL6y2izrCdcmWwMgZtP9AfETen15eTBe3jkqamo9mpwMaa6WfUzD8dWJ+GT68z3MzMelBEHNtovKRFwAnA3HTKcbCh8qah0p2jjYjHgEckHZgGzQXuA5YDi9KwRcBV6flyYKGkcZJmkV30dEtqZn5W0lHpauMP1cxjZmb2MknzgLOAEyPiv4aY7FZgtqRZknYluxB3+XDLLuMRLcAZwKXpF3kY+DDZTsEySacA64CTACJitaRlZGG8FTgtIgYuz/so8HVgAvD99OiUtjdHt4Fral4Z63JNzXFNzStrXUU4DxgHXJu+CXpTRJwqaT/gwoiYHxFbJZ0OXA2MAS6OiNXDLVj1j47NzMysHUrXdGxmZtZNHLRmZmY5ctC2maTXS7pJ0p2SVkk6smZc3a4iO1TXGWm9qyX9XRlqSuv/c0khaXLNsEJqatQFW8F/uxF3+ZZDDTMk/Sh1ibpa0plp+JBdo3awtjGS7pD03RLVNKJuZDtU08fT3+5eSd+SNL7omnpGRPjRxgdwDXB8ej4fuD49Pxi4i+xk+yzgIWBMh2o6BrgOGJdeTym6prT+GWQXFfwCmFx0TcBvA2PT888BnytBTWPS+l4D7JrqOLhTf6OaOqYCh6fnu5N1UXcw8HfAkjR8ycB71uHaPgH8M/Dd9LoMNV0CfCQ93xXYo8i6yDpaWAtMSK+XAX9QhveqFx4+om2/ACal57/G9u9Y1e0qskM1fRQ4NyK2AETEwHeQi6wJ4O+BT7HjF74LqymG7oKtyPdpVF2+tVtEbIiI29PzZ8n6H5/G0F2jdoSk6cC7gAtrBhdd00i7ke2UscAESWOB3ci2TUXX1BMctO33MeDzkh4BvsD2jqnrdd21U5eQOXkt8FZJN0v6d0lvLLomSScCj0bEXYNGFfk+1fpDtn8drMiayvJ+vEzSTOANwM0M6hoVmNJg1jx8iWxn7aWaYUXXVNuN7B2SLpT0yiLriohHybZH64ANwC8j4poia+olZf0ebak16saLrIONj0fEdyT9Dtle7bGMsuuuNtU0FtgTOAp4I9n3kV9TcE1/QdZUu9NsRdUUQ3fBlmtNwyhy3TuRNBH4DvCxiHhGBd55UtIJwMaIuE3S2wsrZGcj7UY2d+nc6wKyUx9PA9+WdHKRNfUSB+0oRINuvCR9g+wmCADfZnuT1qi67mpTTR8FrojsRMwtkl4i60y8kJok/Q+yf/i70oZ6OnB7unCssPcp1VavC7ZcaxpGkevegaRXkIXspRFxRRo8VNeonXA0cKKk+cB4YJKkbxZcE4y8G9lOOBZYGxFPAEi6AnhzwTX1DDcdt9964G3p+TuAB9Pzul1Fdqimf021IOm1ZBdnbCqqpoi4JyKmRMTMiJhJtmE6PLLuNwt7nzR0F2xF/u1G1eVbuynbI7oIWBMRX6wZNVTXqLmLiLMjYnr6DC0EfhgRJxdZU6prpN3IdsI64ChJu6W/5Vyy8+yFvle9wke07fdHwJfTBQebSbfgi8ZdRebtYuBiSfcCLwCL0tFakTXVVfD7VLcLtiJrilF2+ZaDo4EPAvdIujMN+wvgXOp0jVqwMtTUdDeynZCasC8Hbif7DN9B1v3ixKJq6iXugtHMzCxHbjo2MzPLkYPWzMwsRw5aMzOzHDlozczMcuSgNTMzy5GD1ixnkp4bZvwekv6k5vV+6asYA3eDmj+KdX5W0p+PvFozazcHrVnx9gBeDtqIWB8R708vX092FygzqygHrVmHSJooaaWk2yXdI2ngLjznAvsru4fx5yXNTPcM3RX438DvpnG/O/hINU03Mz3/tLL71l4HHFgzzf6SfiDpNkk/lnRQ535rM3PPUGadsxl4b+qMfzJwk6TlZP3gHhoRr4eX745DRLwg6S+BORFxehr32XoLlnQEWTeEbyD7v74duC2NXgqcGhEPSnoTcD6pS04zy5+D1qxzBPytpN8iu63bNGCfNi37rcCVA300pwAfuNvOm8nu1jIw7bg2rdPMmuCgNeuc3wdeBRwRES9K+jnZXWdGYis7nvKpnb9ef6q7AE8PHC2bWef5HK1Z5/wa2f1TX5R0DPDrafizwO5DzDN43M/J7nWKpMPJbjcIcAPwXkkTJO0OvBsgIp4B1ko6Kc0jSYe171cys+E4aM0651JgjqRVZEe39wNExJPAT9KFTZ8fNM+PgIMHLoYiux/sXukOOh8FfpaWcTvwL8CdaZof1yzj94FTJN0FrCa7AbiZdYjv3mNmZpYjH9GamZnlyEFrZmaWIwetmZlZjhy0ZmZmOXLQmpmZ5chBa2ZmliMHrZmZWY7+G/hlZqG/VNy4AAAAAElFTkSuQmCC\n", + "text/plain": [ + "<Figure size 432x288 with 3 Axes>" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "split_ax = True # switch between splitting axes (using pyicon) and not splitting axes (without pyicon)\n", + "\n", + "# --- make figure and origional axes\n", + "fig = plt.figure()\n", + "ax = fig.add_subplot()\n", + "\n", + "# --- add colorbar (similar to make_axes_locatable method, which cannot be used here)\n", + "pos_axes = ax.get_position()\n", + "pad = 0.05 # space between axes and colorbar\n", + "size = 0.05 # width of colorbar\n", + "pos_cb = [pos_axes.x0+pos_axes.width*(1+pad), pos_axes.y0, pos_axes.width*size, pos_axes.height]\n", + "cax = fig.add_subplot(label='cb')\n", + "cax.set_position(pos_cb)\n", + "\n", + "# --- do contourf plot\n", + "if split_ax:\n", + " ax = pyic.split_axes_vertically(ax, frac=0.4)\n", + " for axsp in ax.axs: # iterates over upper and lower axes\n", + " Cf = axsp.contourf(lat, depth, data, np.arange(-2,2,0.2), vmin=-0.2, vmax=1)\n", + "else:\n", + " Cf = ax.contourf(lat, depth, data, np.arange(-2,2,0.2), vmin=-0.2, vmax=1)\n", + " \n", + "# --- fill colorbar\n", + "fig.colorbar(Cf, cax=cax)\n", + "\n", + "# --- makes axes nice\n", + "ax.set_facecolor('0.7')\n", + "if split_ax:\n", + " ax.set_ylim([6000, 800], [800, 0])\n", + " ax.set_yticks([6000, 4000, 2000, 800], [800, 400, 0])\n", + "else:\n", + " ax.set_ylim(6000, 0)\n", + " ax.set_xlim(-90, 90)\n", + "ax.set_title('title')\n", + "ax.set_xlabel('latitude')\n", + "ax.set_ylabel('depth [m]')" + ] + }, + { + "cell_type": "markdown", + "id": "9ab5c2c1", + "metadata": {}, + "source": [ + "Example with using pyicon also for generating axes and plotting:" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "id": "8d70590a", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiIAAAEOCAYAAAC0Mnz6AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABjxElEQVR4nO29fbA+aVnf+b2e57z8fjNAABEcZkhgtwYjUHGz4EjWjYuiYUTLMbvLZnwJqLhTsKBk3wLE3VhbKVKTzW5KUqtSU4hAxTg7aygZ3VFEsqjZBWF4McAgMgUKIxOR+IYzv9/vnOd5rv2j++6++vS37/vu7ueln3OuT9VvTs/91vfp03331dfbLaoKx3Ecx3GcXTDb9QQcx3Ecx7m4uCDiOI7jOM7OcEHEcRzHcZyd4YKI4ziO4zg7wwURx3Ecx3F2hgsijuM4juPsDBdEHMdxHMfZGVsXRETksoj8uojMO+qPROQ3RORg23NzHMdxHGe77EIj8oMA3qGqS1apqicA3gPg72x1Vo7jOI7jbJ1dCCLfC+CdIvIYEXmPiHxYRD4mIreZNr9QtnMcx3Ec5xwj20zxLiJHAD6nql9Vml6uU9U/F5EnAXg/gJtVVUuzzb9T1a/c2uQcx3Ecx9k62/bDeBKAPy2PBcA/FpFvBLACcCOAp6AQQJYiciIij1XVL295jo7jOI7jbIltCyJXAFwqj78XwFcCeK6qnorI75k6ADgGcHW703Mcx3EcZ5ts1UdEVf8EwFxELgH4SwC+WAoh3wTgr4R2IvIVAP5IVU+3OT/HcRzHcbbLLkJkPwHg0wBOUQglfwPARwH8jmnzTQDu2/7UHMdxHMdhiMhbAHwHCiXCc0i9AHgjgBcDeBTA96vqh1PjblUjUjqhPhvA/QD+Kgp/kZeq6g+p6teo6u+VTb8HwF3bnJvjOI7jOFHeCuDWSP23Abi5/HcHgJ/KGXTb4bu3AHgAwC8CWAK4G4AN2w2RNb+gqp/a8twcx3Ecx+lAVX8DwB9HmtwG4O1a8H4AjxeRG1LjblsQuRHA51X1LWVCs4fKsgpVPVHVt295Xo7jOI7jjONGAJ83/996xzO27SMipKyVyERE7kCh1sFsdvzcS9d9VVE+M93n0ioLx42yst1M6rJZ6EvKGu1m7TKpysxU7PmkXT8TMs6ZuqKvlnUrU7Zql5WXLNQ16m1eGF2ZY+1Zpu2yFWm3Iudr9A1Vpowdp/pU57NF2uoaxmuUVe3IOUi7xiEpo13aQ58Zu923fefz+uw+qfEmAFsAQmGsrnHIyvr0YeeTZp2tt2Xhf9Jltro9uNQLQHtedqDZ2Q5mfSNzOLPwNNt39ZFZq0/jl5nN2mXkd6rGGVCm5ptYdVaWtetDXbOsbrcKz3p9tqpsZR7Ietmpy5YrbdQBzTWoGoeVLe060S6r5mXbhbVq1S5DpN2VRz/3pVSOrcfKXE8GPP0n0L8AYK0Rd6lqHzeJrHf8WbYtiDwE4Gnm/28C8IWzjcpf/C4AuP6xT9dnPe9HAQDz6w6rNgePOeosOzBll68vji9dbpcdN8qKvpeO60tyuTy+/nJdFuovm3bXX6rHuXRUPByXD+qtdK47LI4vHdQP0eXy+JJpdzQvst4fzq5UZeG4USbtsoOyDItrVRlOTfTzaVGui5NWWfUTAEL9VVN20i7T00i7kzrYSa8tip9X64z++qipv0rqT5dl32Wr3cqUrU4LqWRRngMAluU4y8Wq1W5p2lVlp+12ALBatOsXZfeF2ZwglFkZrWpXn66qX9rFadXuy45pmcbbsfFiZZtgRvStsbIZeR82yli7RP08fGCYsoOD5k9bb8vCo2nL5odFw5l5lmdlWaizZQAwL9eKRlnZf36pfv4PSLvZcVEvl+pJSCg7rPvKJdKuXAfFnEPMuoWjct06OqrLLh2XY7fLWDscmLLD4+ZPABLqTRkOyywNB3XZQi9Xx6er4viUla3aZSfL+ve7Wj6cV8zzf7U8fvS0fnCvlO2untTtHrlarEtXzDpx1Rw/cmXRWX/lkXpdvXalHOeRep27SsoW5Tq4+Iu67zJSFn5+5AOv+n0kOIHiZ3BzqlmL78XvzlT1eb071mS948+ybdPMBwHcLCLPKH1Bbgdw75bn4DiO4zjnmtms/781cC+Al0rB8wH8mao+nOq0VY2Iqi5E5NUA3gVgDuAtqvqJbc7BmShzo9Hz7DGOk2bOtOCOUzBIsEhoS0Xk5wC8AMCTROQhAD8G4BAAVPVNKNJuvBjAgyjCd38g57RbzyOiqvfBc4Q4juM4zsZYk4ajgap+d6JeAbyq77i7SGjmOI7jOM6GENQ+TvuACyKO4ziOc87YhEZkU7gg4jiO4zjnDBdEHMdxHMfZGS6IOE6M2ea8/UMCOyzi7QLzw1kjV8imsYvDtvJ5XDT2YQGWTUW8bPDZcvaLfXgOAi6IOBeabQohjuM428IFEcdxHMdxdoYLIo7jOI7j7AwXRBzHcRzH2Q1S77O0D7gg4jiO4zjnDNeIOI7jOI6zM1wQcRqsVNON9pWxd3tQH/pGd44znDFq+H16Yw3gXK+/HQj268/qgojjOI7jnDNcEHEcx3EcZ2e4IOI4juM4zs5wQcRxHMdxnJ3hgojjOI7jODvjYL7rGeTjgohzodnkpnfb3uAu93ybnFfuV9i2v9b26evQcdbBPt3zLog4Fxrf9M5xnPOICyKO4ziO4+yM2cxTvDuO4ziOswsEEN9rxnEcx3GcXSGuEXEcx3EcZ1e4IOI4juM4zs6YHe6Pt6oLIhcF2Z+bMsYYKV8abuTriZYJQ24jPNfZHOuOMJARA+7Tl2ySc7Lu7CPuI+I4juM4zk4odt91QcRxHMdxnB2xT5o1F0Qcx3Ec55zhphnHcRzHcXaEjPJT2jYuiDiO4zjOecITmjl7xx5JzueJdUXcxMbZ9sZ7jNTtta7bz2/jHeAXfbK4j4jjOI7jODvD84g4juM4jrMzPHzXcRzHcZydsU+mmdG6GxGZi8hHROSXyv9/ooi8W0Q+Xf58gmn7ehF5UEQ+JSIvGntux3Ecx3HayFx6/9sV6zAivQbAJ83/vw7Ae1T1ZgDvKf8fIvIsALcDeDaAWwH8pIjM13B+x3Ecx3EMMpPe/3bFKEFERG4C8O0A3myKbwPwtvL4bQC+y5TfrarXVPWzAB4EcEvqHLrSMVPMZrVUrJb557pybbHB2Yzg4Lg+lln9ry+zWf0vVsb6EOTS5mROmc2qf/VUpPrH+xQP3sFx/3nNpP63DlKXdZOs67y7nv/aIm9G/G1TCzq7J9m9uzVynuU+z38MthYtrvUfZwtcPs73WOj73hjD/LrDXu0vkkbkxwH8fTR3EHuKqj4MAOXPJ5flNwL4vGn3UFkWZap2rj4361aZ6MOtV5e7ngJlcW2a83Kcc4/9aJoQU/3IXD562qu9FXJz/2WNK3Jr6V7xoIi8jtT/JRH5RRH5bRH5hIj8QGrMwYKIiHwHgC+q6odyu5AyKkqKyB0icr+I3L84/fLQKTqO4zjOxUM2Y5op3Sl+AsC3AXgWgO8u3S4srwLwgKp+LYAXAPjfReQoNu6Yz/pvAPCdIvJiAJcAPE5E/gWAPxSRG1T1YRG5AcAXy/YPAXia6X8TgC+wgVX1LgB3AcD1j336dvRejuM4jnNOmG8mj8gtAB5U1c8AgIjcjcLt4gHTRgE8VkQEwGMA/DGAqJpp8ExV9fWqepOqPh2FE+q/VtXvA3AvgJeVzV4G4J3l8b0AbheRYxF5BoCbAXxg6Pkdx3Ecx+EM9BE5CNaI8t8dZ4bNcbH4PwB8DQpFw8cAvEZVo3mdN+HocCeAe0Tk5QA+B+AlAKCqnxCRe1BITgsAr1JVN9A7juM4zpoZ6F+5UNXnxYYlZWetFi8C8FEA3wzgPwTwbhH5TVX9865B1yKIqOp7Aby3PP73AF7Y0e4NAN6wjnM6juM4jtPBZqJgclwsfgDAnaqqAB4Ukc8C+KuIWEAmGvpxPtlSJPL5ZF5aEU/Xu3NbI2Qt4ZSeu0kdaxfKlgkdYGqTuqre1LH7KjaHdW2ON2Qzu1gf+wEXbdfDoNx3DkPabTTscb4RO/+5x9fajUWcfhDAzaV7xR+gcMv4njNtPodCGfGbIvIUAF8N4DOxQV0QcRzHcZzzxgaEWFVdiMirAbwLwBzAW0q3i1eU9W8C8I8AvFVEPobClPNaVf1SbFwXRBzHcRznvLGhHFyqeh+A+86UvckcfwHA3+ozpgsijuM4jnPO2GWm1L64IOI4juM45wkRYDN5RDaCCyKO4ziOc44QTHd7FIYLImtgpRNy0e6zwV0VyrHmOZBwA7k077/fzJoepNgDOT+YYbkowkSsKlNWxbHdoGy1ZVd8FrWRG9Gy7vPuirN7ru1mDnaTuuK4j9p7LS+Ede6wuCmyw4ymc4NNau1eN3sUceWCiLMVprrpXRBCHMdxzhOuEXEcx3EcZ3e4s6rjOI7jODvDBRHHcRzHcXaFm2Ycx3Ecx9kdh/NdzyAbF0Qcx3Ec5zwhntDs3LGJsM3lmofUMiRua7de9s5hA2aUucEdUz3uSh3JNo1jZSJ2fuu9CXI3rhuywd2QOayj3RCa17j7fLsMUx517/YNy+zzTGzxoqgN413To7DudXXbIftrxU0zjuM4juPsDNeIOI7jOI6zK9xZ1XEcx3Gc3eGZVR3HcRzH2RlumnEcx3EcZ1e4acbJZsymS4oBqrfUhlObCi+w0rk9DmPPTfjGaeY4sXYEIX3tpndDCNMfEn0ypm9jnPLXYg7+qQiZtc0hcoukbp8p7Oe27mCRUaGTuX3pc9TxnI2h75owYFO7QWtZybneuG4Mh26acZy9wDe9cxzn3OF5RBzHcRzH2Sm7TJTTExdEHMdxHOdcIS6IOI7jOI6zQ1wQcRzHcRxnZ3jUjOM4juM4O8M1ItNmkxsZjRladdb42Vk2IDxuY0zgZh+ygVjwKJ+t6nap+JnwqzZCYpfd7WxIbCxMtusS5rbNDcvNnUNuKG8yLLdn1KetW0e7Tc6h2ae4h1JRClPapLHBBJ5hSwjlTa6DOnzem9zLbjIb5U3s7xrjQgoijuM4jnOuOdif1/v+zNRxHMdxnDxcI+I4juM4zs5wQcRxHMdxnJ0xBf+jTFwQcRzHcZzzhmtEnMkRi7QRIjkzaXpIiMK6GCHdV9EMCzLs4Qyr0yJMxEYwCIuEiHjD50ad2HZLEnHT1TZr7LypJtnG+pX75xwWxbLedvYeqO6LHvfjqD0/NvlV23enQjYXtnY06vfnZXiuEOyVIDJqpiLy34rIJ0Tk4yLycyJySUSeKCLvFpFPlz+fYNq/XkQeFJFPiciLxk/fccYRhBDHcZzzhMis979dMfjMInIjgB8B8DxVfQ6AOYDbAbwOwHtU9WYA7yn/HyLyrLL+2QBuBfCTIjIfN33HcRzHcZqUe830/bcjxp75AMBlETkAcB2ALwC4DcDbyvq3Afiu8vg2AHer6jVV/SyABwHcMvL8juM4juOc5WDe/9+OGCyIqOofAPjfAHwOwMMA/kxVfxXAU1T14bLNwwCeXHa5EcDnzRAPlWWO4ziO46yTPdKIDHZWLX0/bgPwDAB/CuD/EpHvi3UhZdSlTkTuAHBH+b/X7n/vHR8fOs8LxpMAfGnXk9gj/Hrl49eqH369+uHXK5+vzmq1R86qY6JmvgXAZ1X1jwBARN4B4D8B8IcicoOqPiwiNwD4Ytn+IQBPM/1vQmHKaaGqdwG4qxz3flV93oh5Xhj8WvXDr1c+fq364derH3698hGR+7Ma7pEgMmamnwPwfBG5TkQEwAsBfBLAvQBeVrZ5GYB3lsf3ArhdRI5F5BkAbgbwgRHndxzHcRyHsSHTjIjcWka+Pigir+to8wIR+WgZVfvrqTEHa0RU9bdE5OcBfBhFhoaPoNBiPAbAPSLychTCykvK9p8QkXsAPFC2f5WqJjIpOI7jOI7Tmw1oRMpI158A8K0orBwfFJF7VfUB0+bxAH4SwK2q+jkReTIdzDAqoZmq/hiAHztTfA2FdoS1fwOAN/Q8zV0DpnZR8WvVD79e+fi16odfr3749con71ptJhneLQAeVNXPAICI3I3CV/QB0+Z7ALxDVT8HAKr6xdYoZ5h8ZtXSX8TJwK9VP/x65ePXqh9+vfrh1yuf7Gt1sJHXO4t+/fozbZ4J4FBE3gvgsQDeqKpvjw06eUHEcRzHcZweDE/xfnDGGfauM4JPTvTrAYDnorCMXAbwPhF5v6r+budJh8zUcRzHcZwJMyxl+yIRvZQT/foQgC+p6iMAHhGR3wDwtQA6BRHfa8ZxHMdxzhsy6/8vzQcB3CwizxCRIxTbttx7ps07AfxNETkQketQmG4+GRvU95pxHMdxnHPFZvaaUdUFgFcDeBcK4eKeMiL2FSLyirLNJwH8CoB/iyJFx5tVNZqUdKxpJuw1c4p6r5nXA3hBWf82AO8F8FqYvWYAfFZEwl4z7xs5B8dxHMdxLLKRqBmo6n0A7jtT9qYz//9PAfzT3DHH5BH5AxEJe81cAfCrqvqrItLYa8bEEN8I4P1miM69ZmyK99ns+LmXrvuqotyGI82lVRaOG2Vlu5n5o8xIO1Y2D2W2r5lC6DOz40joY6dalqE9jh0v9BGsTJm2y8pjkXYZAECD75DxIdKyfmX9ilZn2ptjNeOFPrbdKtaXlAHQMI6dw6rdLhzbKVT1tm91Om01s+3qy2HL2Djt8yXLuqdP6xtVkbLOItYn0je7OtE3SWTNSy6HpIGQOjoOqydlbE0W8uzRelMZL7N924M3zjFrl0m9ALTaNc43Q6usOrYftWRNo33oODNSZtvNuvvaSVTztxNrn1dNH9VZu6w8Vm33aSxLdJko/mdJ1omVKVuVnVamc2Ocsu3Sri2kDysL59Flu52SviDtwvGVRz/3JVX9SqQY5iOyEya514xN8X79Y5+uz3rejwIA5tcdVm0OHnPUWXZgyi5fXxxfulyXHV8+KOuO6rJLoV19Sa4vyy4f1xakS8dF/fVmvMvHdZ9LR0XbSwf1TXD9YSgz45T1l027cJrD2ZWqLBynyg7KY1lcq8qwOKmPT4tyPb1q6q91tsOJKbsayk6rIr12tVVW9SFlq0cXdd+r9nhZjrcwfQoJZHUt3k5Pi3aLq3VevNWiOF6asuWiHO+0lmyW5TiNMnMcyleLdv3CTDXIYwuTmi/Ur1bxsnC8XGqrzB6nymhf8mSxPjl1OcS0uqyuKdB397F1sTIAmIePDlJvIxlZWXg0bd9QPz+sC2fl8zozZaG+UWbWhFA+N8/6/NK8HK9eEw7KMjHjSDmOXKrbzcLYR/F2cqkom11nftGjI3N82Fkmx5fa7S4dt8c5NGUHZdlBXSaHlzrbqWm3WF2ujk/L49PMsmvm2btSPq9XzXN7tXw4HzldmrKy3UlddqVcEx65Uq9fV816c6U80SNXTf2Vov6aKbvySLHmXbuyMO1Oy7q63eLR4njxF/Vau4yUhZ8f+cCrfh85XARBBBvca8YZQR913JgbtWdoWOOrjGi2bFns41wy26XnE+bf/w1sf3X2Ag/1to6VsfEYbJyu+qpdeZmsQNJ3Dl3CSe6fPiWAxNr3FU56zaGn0JSLjOncGCfvGWaaX6Yh7sW4C5DZbjMmA+cM8/0Jih0z02qvGRSmmRcCuB/AIyj2mLkT7b1m/qWI/DMAT0XmXjPKPu16sBrZv4tHrpw2tCLROWxmCg0qdebsEmarq90N7ULDEuyzxWQNK7dcmldajSGEhVZP23Wzeb2wrRbN9gCKDQU6mB8fYFV+Kdk+spLGT6A2wbF7KiWc9G23ybGHtMtljJAytTkwc231wp/3f+E37slwDlKW6rN2hkhzMaEjd01YnmA1u5RuOJK+6+/1lw8bWpG1zGENL4H5dYeVViSLi6AR2dZeM4Ok+i2QK4Rsm6gQskPGCCGbZHU6zXk5znlnG0LIENYthKyLXkIIsFeap33Ya2aSDJFwG45RzKOx93jG5jzmnsuVnFNCYaU7X88DkDLDMCGVmVzCl61uQKitzB1Ensk34dh55d0XfU00AP8y7GOSyelL2+XeNgM+4BqO5CM+0te9PxhzhG/Wt08Y++ha2wdZyjmHtYuxpq9uu5aNG0cbP3v13YbqeptcBI3IthhrmomxWsbHjt3MV64tGk6qg+fQ8O6OeIubslV5zJKwrLpMM32l48TKLBLmlTeOXDerHU3n1shfvvWsJBWO2QeA7Xva/cZMLdxB5T2f16YZZsKhEVgr7lwaiPmIWIb4i/Tt25gDuSS5wkku6xI6Uv4guULHGB8R5nfBqMyGqxVmh92pkbKFiXni4oTnwz4zYWzbtzJnripn1hSyBpPLmQG7h1tdpVqRVeY6GI5XOtLT+gyXjw8qx9XGvCLvg9S7ZB30Ns1sYPfdTTF5QWSqppl1CCGbYLKmGfJgTwE3zewfU11fY0LILskVQrbNVE0zTAiZAn1NM1GhcmJM8w7dM5ZTVemlbkRW39tZdUuC4gi1da4wa1XowYyjiS+d6svdhs6eqQPyNSdDyI2GyY2uGTSHzGiYIfXZc+j5MZ9rpbCEe4mZW1J9+tYVJ9qGs2qmGsuSu3bk1O2Qya7d62Ci15wx+Zn2Mc3oSkebcq5eaUvDj1xtl3VJzSvVtfh/nMXGz1dl2i5bzS4htm+AyKz61xuSClgOj9r1ly612/X5KptL+W9W/4u1M8hMIDNp5HEIHBzPIXMp/kXadRH6zMw/BsuWHMrYztwHB3W9iFT/xmRljs2h66U85l/fOXTBfjdWxq7ROuZg/7bh783WlCUxC65Ol1Wf6j7rcN4K7RqQ+7muM89C2U7JHJSsVQ1tZN8baIiUuDhtrzE2f1G5Ps30pN13g1w2OVvCOm1zRFXttqRBCu8rmwsrEHJiWVi7KJvZa2YjTF4jsm3TjE1oFrj+Urts26YZm7ysKpN22bYfbj0l5ztpl23bNLMk51su2gs3a7dJFuR0YzUifdn2+XLJ9eeQLUcDsDXIJjmryra8JgiZg5C1atumGbHJywKHxBdE2i/bTXJl0TbDXrnGyra8VhGTi01oFmvXjQDzaZoKGZMXRKbAcogH9jZyhxh3VaWuq4awoq8jd0hO/dl2LImZPW6MR96UJPFZoPHVedouq3KBUGdUo7WZ1X+0ylHRRLTEtG0NE8iSlK3aZWfrhhIbc5dziNXlOp4OmUOuQJNtjrH3Ek0cFh+I5R5hY9cT63hWWu1m7XbMgbXPBe77h0x9RUfG67V+rYEha/KQtX8y7JFpxgURx3EcxzlPCFwQmSKpGHHN/CQckz9kjO+IjbMPWVSTu4nZrxGqCWFfYz0/HTcYwpCbR8TWhUyVK6L9YM6oQ+Zjv2JDKC8Lk7V3VK5DqSXkF5nNEmHmmZqOXM3JGIaE5bK+rF0wyfTRpqScVM+2G5IxlbVPObOyjKpbySNCJzPkIlKb2fBzGML6Nia3yJicTX3W+Nh7Y+d5SVwQcVLUu0SSHRqHCB2BIZEyuQx5y5THwlTG9pjUU+HDtAvp3uOJzYDYfjLz43m1sV0jnfuqnQSNRdDUOzC3c4vkmkUsY5KJpYSO3HbrSmjW10QyRBhKCS+5c4kJDnSXbmKuSZGbxCw1HhOIYzuPd0pum/rA2MAaVG1hYdZG9oG363f/pNjgR+K6cUHEudCsiAOr4zjO3uMakf0lV522c7XbGXprTizMhLMuZ9UxkC++xlUPYb0mjJFrRyJfiz0c5KkDK9WSFD9t2vdczQNvZ7/I26Yg1jdXKzMmuytjXY6iNkz3bFmfW7Pvuek9kkmXWSd2TzYgYerCNIabZIzXcU90TWnd18W+rv2dXJS9ZvaVdafj3faNGTzMGw9y6p7LTUBE/UYybb9DdOxh8bUZwWIp3ukCTqJm2AuFjDc/nmNVCjIsgqZhCiJ/51ja98ZUiY/IkDU85jfCxs4t2wRjTDMxf5A+lgYWtMXax4QObgKZtcpS/bnw0fN5s30a0TWZf8ghzyubD1sn+iY5MzTM0VuIoAlsYu3eRrr3LFwj4jj7wSqyX43jOM6+orP9eb3vz0x7kBsBk2JQhEx56nVlV10F7UfmV4IaKZhmT82VklO5RXK9+Fkug7n5+8Ry9MzJF58VHGKOfmTcZrRL2xmV5R6ZWTMMMfEELcms4dBItBUkt0hVN4E8IkPm0zeVTFfZtvOI1O3iJhOWzp3lBIm16yKZPyQQ7vtc00zD6TuiDuoilntkzNphyqo1KrFE2jVvtQYtScOpdcAztw7tybreTVnnmpjpK8beCiKb3JX3LGPvndq7m5W1Pb7tQ5f7+GkqW/+YfSAibwor7AwJieWJmMKYkcRmADXDhKs5a5hUusebH85ouu4YTGCxviLUBEJCes+274P1G1ENwtCQcdYzn1jf/OjPuD9IyjRD55MwyVTnJpEo6yKYcbJDdlNzmBHhZMi8w+884vkfEyGTXLMIdm2MRR6ytXYI285CvM732jZNXGPZW0HEcdZBXyHEcRxn+ohrRPaFMcnJpkBaC5L5ZWWTncW+cIbo4ulnbG7ehXhukfD1qqftPqlIBZqLYWXrZ+XPdtp3+9VCv6RX3RqKRG6yQQTtiGpCKzNRZ9W+jqkpZ9Vcf8xU8rJYOnfmDN0nUiY3f0jsfqY0npkBKqTcdn23hUisRUO0I5tiyBo/xUiaKV3TFBdaEBnDtgWS3CRnth199Id4u+eGOgx4q1ULLNlfwwoYIIt99RewEQNhQWDCCRE6DuZzLMjGV1W7hF9JMMnwrJr9I2lySYX51u2G+zkNmUOqnm1YNyahGR0nYY6ZZQqrLHlZjK4kZ1HzSyxk17Yj50kJNklyn+EqljrTbyQldCS+1Fnysm0wpY/MdTAmM+22cUHEudDEhBDHcZx9xX1EJsQ2ktQ0vbH77SuzJE5VQ1K8B+m38+ONfLkEJzVNtKtYk2mGpWlv9icOrKGd/YKcMTVD+xxBk9Fn/5naDNPuQ51VV20TjvWYoxqMSCRNitw07VYDEXNqHWuuiWvlu7UgXX3HOKum6mO74OZGuwzZVyYZKTPr1pLQXXUj8ysn1vx59jjGANNM3yi9Pl/ssRTvuWk7xqzTnWNONAmaAljp/rze92emW2IjCW6od/fw8YLtb2VsgFT2TdlxgzIgtfndOt4UzHYN1KG8dvFlDqQsLJeUUVNPVWeEodJk0UhoRgQM1gfMR4Scr/mC6t4cb6wpNz8stz3HMRE3jJTQkVs2yEckchvPiBkj16eDJS9L7vNCfsGkKYUJKjHByD4zdF8ZVjZCsuO2tu72jXbxG8yuZWN8G9a91tJzTNAfhLHaIx+R5ExF5Gki8v+IyCdF5BMi8pqy/Iki8m4R+XT58wmmz+tF5EER+ZSIvMiUP1dEPlbW/XNhq5bjbBFPaOY4znlEddb7367I0YgsAPz3qvphEXksgA+JyLsBfD+A96jqnSLyOgCvA/BaEXkWgNsBPBvAUwH8mog8U1WXAH4KwB0A3g/gPgC3Avjldf9SjCFpd5cTkHzD10HqK0FTic9yI2RYu6quQ6uxDlhOhMY5Qhr2trNqMw172YfsP5Ob9r3RJ5HkbBmcVUkq+NxImkZykRGXdUiEDNOSjGGIA2tv5VqHFiQ3QiZWlpu8jNFnr5mgzYhFypyZmOlLTJfrIubx2zWfQGYkTXKtqtqtR0uyLvq+D3ad6v1c+Yio6sMAHi6PvywinwRwI4DbALygbPY2AO8F8Nqy/G5VvQbgsyLyIIBbROT3ADxOVd8HACLydgDfhbGCyBYyzlhV3CbVcmHo3IcuuddMtpd7QuUa9RsxC27wOekTQpht7y59MdhGeAk1dxWWy/aauVSbZpiA0ZhPJKS32S72govvEZNiHfLf1MJ31xSA1TFOtxkmHYLbPYlBWVQTPiJsg7vsRGsxM0zH81b5dKSEjtBuzP4yCdNM7te4XRtXzIFrDWxrvTcn2ciw5zaPiIg8HcBfB/BbAJ5SCilQ1YdF5MllsxtRaDwCD5Vlp+Xx2XJ2njtQaE5wdPzEPlN0nF64acZxnPPIFLRIuWQLIiLyGAD/CsDfU9U/j7h3sAqNlLcLVe8CcBcAXP/Yp69dJF2HlLvskJpZiuFlZor3GE1v8eKYfSR1aknW8eUy9pOVUH3pJZKXcWfVthkG5Gs3lva9od1gX8g2lToxwwQtitWg0MRnYZdeM/0VM9ckzDS5jqnroOtDbUzgRV9fyD5+z3VZXFsRM8mkNFyxHB5Wg5JM5x7bV6ZRFnOizdfKZJOb1GXEepL6Ul/1zCOSSvEeXX9N4bpM8dNwZpXzl0dERA5RCCE/q6rvKIv/UERuKLUhNwD4Yln+EICnme43AfhCWX4TKd8Z29yvZii54busXWMDvOogsdovSbu+wgkru3QJODlp1zPCgjyLRM8AJqS3bYZJZVsNx/PjA6wWy0bf4tRtsw8L6Y2pzlNJpmap+y+cJyKQAF2+H+26vpveDfH3yG03KCgrMwgkHdGS5/sRC9Vlycu6/EKiWVRzI2UaSc7W8Dx21TMq00ymZMmED1sWTNADwnenyhTfJfvkI5L865aRLT8N4JOq+s9M1b0AXlYevwzAO0357SJyLCLPAHAzgA+UZpwvi8jzyzFfavo4550ghEyMIIQ4juOcJ1Z60Pvfrsg58zcA+LsAPiYiHy3L/gGAOwHcIyIvB/A5AC8BAFX9hIjcA+ABFBE3ryojZgDglQDeCuAyCifV9UbMLIfb+3elTkuaa5j5JeJ13mkXpI5mmc5ngTGJHlKhDnPztyOajtDOmlLqT+S2GYalfacaj1XdjplpUrlFgtrCOrquzrS32C+nmANrgzVF1URPscGEZqM0In38nmO5QIgGg0bIjMgZMiPnAJCdzp3mpKmej4TmZNgFyytj5KZ2J2VsjbJrWexLPmWG2Taj3hsj3ldJdPpaJEtO1My/Qbdh4IUdfd4A4A2k/H4Az+kzwW3DPKbHCimx7aqzxzAPJ3tMw8N9urqMw9mVooztO5PIiKisXe7iFOPoiGtF2AIaXgrLdhnbfwYrO9d2ttVYkrP54Sy6Ay/LvNo0wxABKQg+7MXI9oNp/F9/c802/Ua6GCOIsLLsIJCU/wYpq+qIMMGEjlTmVAY1zaSyqGbuK1OHtyecZcZAbV7xtaMuzPMLSTlTrgaYFtjH3JjlexPvg20x5PrtCs+ses4IQsjkmKhpJiaEOI7j7CvnSiOyj+gGE8kwqZh5Xo/Zw6BhrimlWqoFIREy1lO68YGVm7xMyKdoLLkZ+YxtaFhSHobl4DRnCDtfI0lYecx26bXRLiR/SCrygGku+kbSMBEn6cDa+L/IPbQFc00f1qIRybQqdEXF5DqmjoqQiSQ5S/1tqbnGwhxYM81NdXt2wcxcpf285ub1abQbkdqdRXMkI2nKFXCUdkPJOr1ljccm302tc21oYRCRWwG8EcVr6c2qemdHu69Dkcrj76jqz8fGPJeCiGWbf/gh5JprUhvhhX0FVnoZh1KaZpj4MiRCJiw6KXNNNGrmGDg57a63MJv6adv+UJlDWOgvizKw4bSlsHAwrxOaMX8Re//MyNj1y6DdtyFuEYElBTGYcZjEs2XhpLcgkvsOJWaYVAbTZMZUEm6bu5ldLFS3U9CoosDyhCb+TMRNOL2fx1R9Y01gvh/d5ho2f7YWde2Fwja4Y/RNgbBL1vEeml93iOWjHamg2Tk3oBERkTmAnwDwrSgiYT8oIveq6gOk3T8B8K6ccSfwLeWskyCETI6T/Adom3hCM8dx9oE+QghQaJH6/svgFgAPqupnVPUEwN0osqmf5YdRpPz4Iqlrce41IoHVFlLBd567pwOVbRfLI5LyPm/QN3lZbh6RIaS8EsMnPvnCpCaXhrmmHCeR5IyVpRKana3rrE/06UsvZ9bAlrUkUY1IphY/1wyT/JuNMa+kyqIJzfi8osnLUgnXounhU8/RmohpP1i7RFmvaL9Qr0HjG7//1xEYsC52+c4BBmtEDkTkfvP/d5UJRgM3Avi8+f+HAHy9HUBEbgTwtwF8M4CvyzrpkJluk00kiqltgat2WUKFFrL0PXL1FJePDhp9m2ObTpnOy7FTpzaACjfdqV7C4exqowxALcikFouwkK0y2w0xzVgi4b1sMW+aXMpj++KPJTmzv1PIajqfQ8leMyzzalCvskgaS622r2+CmN9IH7LNNYzYyUe+v3JdC3LLYmaYTr+eTH+Qes+gdt+YOcb2paG6zBzTVV+1I/dzMqSX+XZELuhY00w9mf7tSmjUTMeLMnejz77WDrsmx3z9jg7nuHKy6B5nydb79rtk3fQxzSgESx308bNQ1edF6tmgZ3/hHwfwWlVdRjKwN5i8IJJyANsVQQiZGkEImRwTNc2om2YcxzHEhJBd0ts0sxl5qCtzuuV5AO4uhZAnAXixiCxU9Re6Bp3m23QibNIZKjcZT+yrgEXIdO01E1IsN8S62K6arF3CIz/mpU935LUQx0FN7ivT1lrQJGekjOV2SKYBzzS55DqwjhGBBplrYoyUx9aTRyTTBNLlrMoSlWU6po6JkIlqLcxxdvIyC0uCVmlOUmql8hwsUsa2TYUrRU24Em1XpXY3tyjThKT2RanXwfaNypJCrospmHiGoZua7wcB3FxmTf8DALcD+J7GmVWfEY5F5K0AfikmhAB7IIhMMYc/UEjNuVqRXKEj1i4VNRPsrqeNqBmmNo0vHJXgIGRxIi/iJEFYOL4EPSV7zVQLH0ntabOthnNTkwsz17TNMMJML/MZ9DSy10xH5tWz7ezymOs30r209qNLwb1N+gsicaFulBlmgD9I3BRU/wKxCBldKeRw3llPI1+SJh5S1jOLqi4XkMMj9II9/6YsmsgsM5vqgZzgVC+320aiZnKFjjHCyeWjg0lqRfpGzWzi1amqCxF5NYpomDmAt5TZ1F9R1r9pyLiTF0TcNNOPqUbNVELIxAhCiLM/THVNCELI1OgthGwJJoRMgSkKIcAQ08xmPkZU9T4A950powKIqn5/zpjTfJvuGXT/g4Y42r1wNm+WdruqngzRx/s8fF3QLxn7FRU0ALke8sl9Kpbd7Vptz44TN6UwR1GaHp45uibU6blbrVdOjqYs14G1Oq05HqIdCXNo7GNDW3YvTI1bIHMSuQEaudqPRn1mFAs1wwxwTI39vbMjZHrcS7EddtNjD3BMZWW9zTCkXSKJWW5q99ROsXHtR94Lt+GYukc5SIayT7/b5AWRKZhmrlxb4Lj82gmL9CNXF7h83O8LaDkgfDdIIDbGu8oySPaVOdW6zNpdZxJJbsYiZCwsfC8miFy6VKd0r14ER8Y0k1Ipt5N/VeYXspkdS3KmJKSXmmHmB7XDKkl4xiJprLiQ6zcyOwx9ZyZ3yapqX0XkmD7rdqOdBxNCItxgXRGgsZDlzlBX0re3GcaaUg6JIJLtD1L0ZREyhRkmmBDDc2TLrNDRLqPJyyImHKGmGfLsLRbA0SEsDdNMriCSG77PImnsucu1aiYLLFZNzcyBXMPpqq0VqdY3lgQtM3x3mfkSvnR8gKvXFo2xjw/nuHKtWCd2ta8MM8P0Nc1MPJdng8kLIlNQwx4TlWtfIWQTsH1lJrHXDNlXZgqmGWaG2XbUDEugNvXsv7uEPf9DNqFbN5XAkSjbOmeEEGAappmzQggAKoRsmyCEWIIQskuYwDEV08wmmLwgsgmUxZCXLwNdUzz4kORlrF3McZXuNWO+qWdM3ZlM506+kpi5JpL/I+W539xDJtOcQ5KcVeYXpjlZkXa5kTR23rS+vwNr1Z7s3Nul+4jV2nkxrSEz12yb1I63VVlsh9xMLUhRPiNlEW1KwjG1d4RMIgKLJuhjWpJU8jJWRuqrZ67LrBPLR5LSfkadVeOmF5banTrjG2Jr45AkZ0MI74PGO2LZfpdMQYsPAIuJzCOHCymIOI7jOM55Zo/kkH0QRLRSXTdU2MtVq0wriZVlvtvcX2VJzpeils7bXz9McmfaD5ZttTOPSNCi2C++6sD2IRqMzOx4vbOtApU/SLZmxWoCqs3E4rlFKn+RzJDervp5eY2XxryS68AaejCfkpRnyLr9RlLalHWNTesj1yt7szqiBbHl+blF4n4j8+BfkgjFTWVWBdP4sBDcedufJbpDbu5z1sdRnEHWBLoRZonatYP4n1VrUWa2aEta+9HPcXW5wfcCy9qqRHOSfK8NNN26aeYcM0SwieUH4fHuQspq2EMbHFPtWthQd/bdVZeZSRvmmu6FiPbpXAAjUTVGmKidBM0vmJtbpDxmkTQ4rlO8p6J0qq4JB9ZZ2An51JaVGD8CvuFeXOyYE5NLJeRYE9WO/E5Se+r0FTrY2Ez4sH1m5hozc86MCBgxx9TcCJlOB9wROUNo8jJqzkS/dhYmYKT6UBNO+29GP5BKupKY8fWtfT/3X1fjbOOjdVsoXBBxnL3BU7w7jnMe2Sd5ygWRCEPUd7l9YtIq+5hNqi4rc03tFJbqQ3OKVJXEMS1VFlUPm688q7qNqplJtlWrJYnlFskM6WUhu4DRdNj6s+dFrb5fNhzaV63xql/D/k/Z12ovVjY0eD4vfyerzo2bbpozSGtJ1uXUGjPJpLQfbIyU9qMa56BdnzTnkLHnLMyXOpnmmWuSobp9c4YkQ96J+ST2bNnjIc86I2RlTuQOofmPkvlGijs6peiLrqsDTOhD+kwF14jsIWPUcv39QgDuG8L6tMtiKd67fESCGjSaTwSIR8isyOKUEDroORJmmJhAIybteyy3SFOoCD4iZtEvzTRyfFCF9dKEZ4dWsGungq/Ga/zO4aWWyNcRptwos/Z11r/dqzJbmVWavjLKObL7lb3ku8w7KfNLdbqYcNIR+VLXM/NJt2+HPV86KVm30MEEjEbG1EiETGfCskha+FE5Q/r4g5ydf1ef2AdGZkKzRoRMZNfdriRmfaNmGEOSnOXCfD+myISn1sIFEedC4yneHcc5j7hGZE9YDXDoy9WcNKXm8sB8AMRi25Nl4QPfDDgjWpLsnCJUDds2pSQ3wstN+55thkmNU14TkluE7fabG0kDGKfXujbbgTX8JjPzJb2KCDxtBXQOMYMM15LEevJT5Gk+ukjlBanbtWeUu/tuOqomOLDO22UDHFOjETJdaeuZeZKZfXpHw4zQnBSTKH+ScdgGd2zDzEReIq2ypHZrRor6vCyqqTKqOQmJkzPX7EFakgkmJfQ8IhNkjDqtV4RMpG1+6JnEy0La5FRIL6u3kTThBZ3YkZfSN+27hfiLaHIBJb4f83a0i5p021W7WCTNfA6ckHGCzwCsKSiUGUpn11kjq+bK/LccOmwRkNDApISSyqej4TcSjuPCSTXGBqNr0lEz/YQO1q4ruoaNE657oyyMY/5mlRBj/47MHyTblyRPeKFmpgGhujR52dn2Z4/7pnO3EOEkRObFQnaBtO/aGKEkpw7oYU4fYbrZtdkmN839FLgwgojjUE48asZxnHOGuo/IpBiVpn1NUnNdX389LDNj4FkZd1adle3qLwuWU6SRbEjJF34g5cAa8n8MSvtOjjU1TigjTq3GgZXlDOkdSYNaE8IcWBuah+AcW5dEE5rlmmsa/c0cYvcaS1QWIm+ArqiZ9jgsMoefL+/rOZbErKuMaT9i7Wx50gwTcy6ds3uXaT+s5oSMl4iQYTlD1uKYmmvCscexdO32OPH3ponKIvlDUjuFD1kbmRYgqqVek7Z7k30Hn9M1IvsH23+GEf64XWqv2PbS/CHS1nGqrJoze+Abwkk7lJcKKkMiZNiOvbEEaQlBRMx4fUN6G+r0kIU3M5JG5nVCM2mEGDb72nqbe6QyldRniyY0Y+Ya++JnJhcWDWPHYf4gLCw3N1TXCi9j6LvXjCU7AobtNUPMMCyhGQ+xNWWRnXSbJpy2UJHMxspMQaNCdTOFk0QCMvqs50bSlKyIGSYVsptKYtZ3bcxday3LyHiNcSa4rwzDBRHH2RM8oZnjOOeRCctILc6lIJLrgDcs/8d66nPVjrHkZsmEZh05RepxyrbJhEbtLzBlqaA18uXUqRGJaT9SuUW0XXZYhSPV0yeRNI3cEGH6NuqE7dgb6qxmpbyGYswsteuoccALTq3mfCQBPJAZNzO38wravKye3ISzLobsNcO0H6x9cq+ZiBlmRjQYjfMFp1aW/yPhmBq0GjTfSEefKtX/2swweU6tSZNL7r4yZJ3oygsCnNV+5GpH6juarYP5ppt1BBDEyX2XbGvrBYVrRHZGtgCyhnYpMzq7MfPD0ch4ao/L/7HuDZFImqK+KJ8RQaWxEV7MvMLMNcy+PMROTYUOc5EPylvVXvhQzyJpDu3CVtaZKVbmmOOD+piE94pdNEnis3BkF9cglAwx1zChpBHlEnxbyE3ChJOUf0mfRGY5fRkpPw82XmxPGiZ8AAPMMNaXpIpYIXOwAkbmBnYNk8vRvFWWHSET7vtBoboR4cSSmzGZCCdsg7tGQrNIhAyLjgGsyYXVdbdPliV9+aLVa3lvWDYtlCwWLog4zl7gphnHcc4jrhHZY/okKmOwPWZYWe5Ou9whq+iTSmjGnFkbZWRvCIk5pmZrSVbtMmZSseVG+8FzixAzTFSbQpw0c1PBwzimWn0FySNSO7qaskxzTRXZQnOCAKtgUmrcPyQvCHNMJQ6sbOdeyshEZmfn1SjLjaBhe81ETC/NPnHzCjXDRPKIMDNMbp6QxrnHJCBL5QyJ9e1yUB2Tzp2sHcxkzNads3VAOnfIoMSPJblrcjXGyHfAVOizR9qu2XtBJNfGnb65SPrT7L55asCG6q88Tdq7u71Ihz7cR8Q83KY+7DHD6hv7z6wj26omFrHcLKpK+uRG0rA9acheMzg2UTNkXxlQ4cT8eqG9jaSpxoiba1Kve2a6kVn3RngsURkP6R1ujukiZqYZkmGV9WGmF2bOGWSGqQSMGSkjgg3xB+m110zPCJlRobpdCctm3c91bhZVtq8M2+CORdJYGkJJxCSTbcomis7cdbpffbdGNdtvZEMCw74ITMA5EEQcZwxumnEc5zzippkJoSM0HUP6ZmtHMiV71j6WU6ShFpW2Y2ryM5yRa66JqXpTCc0aDnql6ryRX6OsPzC37IppSZgJh2gHQj6IulVTKKlyUhjNQ/mzkSwN3eNQc43NQRB28wWnNjW0TTdUS5Iy4UT2n0lFuwwhriUZp/2gfWPaikwzTNKsE3NMZeYY2zZXq2Hv8XAfSvv5GBRdk4yQS/TpCdtpN2auAeL5QZrtuvs2yrI11rnvA3LeEX03hWtE9piUxzP19xhUVixONBMgFTqIOp34gKyM8MF80tlGeMn9Z6opJxasFSsLqvFjYHHars9dQMMcWSSNhQgnwUxjr1f12x3OagGEbnpHzCtHZpyT7rG5uaamemHa36mZJvLsb1eRCt/lwkmoq4/XvUCmsq3mCh1V+1zTC1ALAnT/lkzB4YiUMbNOKslZboSMvYdHRchEBAhmPrX1qeRlVZ2558i+MrEN7iwrtu4kQthjwkluNtXcspRvxRQ3uGO4IOI4ZwlCyMRw04zjOOcRF0S2QMzBJ+X8E3NM5e3z1XhjzDC11N/+MrRCOHPcCpoQ+xs1zTTEdBPMOWz/GWZeaTihSqudxhIfzc2tplarEbx24w6sIQV8MpKGaUnCGFiYKUQcWAG++y6Lhim/oOXEtss111S/iJnjqlXf1ApoOf/2DUTziBBzxqqRCn496dwZsxEOrDQJGOtL9oYZZIY5ImUsuobsvptM8R7uyYN5uyw7QibPnNkoC89cl7NqbvIy4vQa21cmFUkTWDUcXa0WuPiZSmI2JLV7VbYmcw0j5sCaejety3FVAZwu9+cja28FkXUzSpU3IrQs98HKNdecri7jQK4BaAoYLEKG7T9T1ef6ftDwXW2XzWfAatE9DrWVm4U7PNzZkTQJgSX4jRzPgVJI0COjgj8pzTBWxR42wmtE0pTTyjXXNPZ+qf64VVkjWVq1IvcXTmL+IFRgWRNDMqv2Fjq6/EZipptcMwwTJg5Z3/bGeumkfamyiC+JfRb6hup2CiKR5zozi2pqXxkaSUOEEgvfVybWrrsuVdZnk7x1CSpbQXW6cyO4IHLOCELI5Fgt0m12gZtmHMc5h7ggMkHYl+Em4rez1YCxRGV0h0miYjddK7UoUZ8W5W0HsVXZNjvte26SMyVmFpZbxNanHPRikTTJPWmIA+thOUbdCmKEkjqqxp6vHZXBIl8qc42552odVlu7YaGmG2M+ydWSgOQRqdoTzckmiGpHBmg/qnZMC2LGZOOw3XKzI27YeA0TTjgH0VrYOWZqSfIjZMgzk7mXzJA+uenc2b4yKXNN7rYWKU1H3/V3XbB3ybb2lenCw3d3xDoEizHht5ZUNj/mplKbYRJzJA9leLitacZG0Mwigkp2tlVmcrHJkGbF7dSIxAjtDo6AZakVYdEwjbBcZkppt6vNHUTwifiKAEBwF5HjGbBYlONYk0u4NuZXCeez0TCH7QRjPHkZeXlUv6dNvpYw3YSDhHASzkOjZpjwMXbRzBRoqHDSV+hIRM0k/TfWnOSMRsUwfxASlts0M8XaxZ+F7AgYcxye1+zkZQOyqLJ9ZVi7lfl4iW9Sx4+rsnKYvtlUu867SXNM3/fV7PgAq2v5muV90ojEjXUGEZmLyEdE5JfK/3+iiLxbRD5d/nyCaft6EXlQRD4lIi8y5c8VkY+Vdf9caIyoM4bJmmaWEzXNLCY6L8dxHEMfIQQoBJG+/3ZFH43IawB8EsDjyv9/HYD3qOqdIvK68v9fKyLPAnA7gGcDeCqAXxORZ6rqEsBPAbgDwPsB3AfgVgC/vJbfpAOmHqNSbm47KjUn+kT+wHa8mNNrs10w4ZB2RAsC1JqOcWnfE3vIhDmmcosocUJlX5M0kiZT+2Gf2ZgDa0PDYqJqws63jXTVoZ35VXruNWOpNCyJqBkQx9TGHz84phJti/0zVjBv+g0kNGtEtETORzU0A6JmqNPrGvaaoWYdWxaLgLHH7B5PaEmyI2Sq55GYh5gptNGHaVHaWpIx6dzZvjJdzqjL6haPm1yWEXONJbYm01TwHeONeW/swkzD3idTJUsQEZGbAHw7gDcA+O/K4tsAvKA8fhuA9wJ4bVl+t6peA/BZEXkQwC0i8nsAHqeq7yvHfDuA70JCENEJ+BKym+3KtQWOD9kqX/ZJCCT1MfP9sA9gt9DBWKzmOJgtO+stdRI0I5xEk5yRBa0RNbNqtdOwMM7mwLLMJZIKDY5FzbDIl1V7rkm/kerOnwOLZau+9iGxSckyzTWhzCZDq+ZvzCzzcv72/mKmG+IPgpQ/SGwxTIXs5n4ZZQovUUHDkil0JPd0YfWNBHcjzDCH5F6KRbbY45Q/SG4EWUzAoP4e5j4cEqqbmbwsFQ2TC4uGYdGDtG8lYPQ3qYT648M5rpUfEVNKXtbfNLPByayZXI3IjwP4+wAea8qeoqoPA4CqPiwiTy7Lb0Sh8Qg8VJadlsdny6OMyC68UWJCyC7JFUK2znKaCc0qIcRxHAeohJCp0dc0s1jsjySSFERE5DsAfFFVPyQiL8gYk30maaScnfMOFCYcHB0/MXoyqvJqqJ5LyT1TQk4nQxsuaafziKDVLu641dacsCRmtlwbuvpFu09ukrPc3CKxMqA20yQSmlUq7EY7oiWpnF8RJ/iGdKWML+uDZgToYa4JZpGGmSX0iad9t1RjWg1GrmNq+IpvJC8jjyC7d8dE0qS0JCmtR9WOaD/YGKkkaEHzkExKlmmGYeZAFhXDUrcnc4FExknlDJnFzDWkrGscljMkM3lZVZdwTK00sZlaXlvPzDmpPCJR0/iA9Zy9I7LHIWbRTZltzlvUzDcA+E4ReTGASwAeJyL/AsAfisgNpTbkBgBfLNs/BOBppv9NAL5Qlt9Eyluo6l0A7gKA667/K/aOMq3Wq5Hoe0NeubbAEdGKxFSDjXb0weprrmn3vbYEjubELlsKIMxvpLFwRJKcqfERkViyJOuZH44Pj6GLk+4+TJhgIbiJaJj65W6+HmKJz46O6vvKOq6G85iyfHNNmdCskQgrrJr131Yr05IVKtqmm4YwQcwKleBjb8ew4KWSl20wfDcQjZTpasfmNUToyE1ylmuGYfvBJKNcmHkokm01JeRkR8gEYXQBOTgi9d19mh8iecnLYv4g7GPnYAaclPcpFzDGrKGkPTWN14SomqPDOU6IVmSrCc2IXWUqUTMiciuAN6J4Cb9ZVe88U/+9KNw0AOAvALxSVX87NmbS8KGqr1fVm1T16SicUP+1qn4fgHsBvKxs9jIA7yyP7wVwu4gci8gzANwM4AOlGefLIvL8MlrmpaZP9y89UdMME0KmwFHMQXCHVELI1NgnQ6oDoEOwmQKzaT57lRAyMU4mmoKcCSFToI8QoigEq77/UojIHMBPAPg2AM8C8N1lgIrlswD+M1X9awD+EUqlQowxeUTuBHCPiLwcwOcAvAQAVPUTInIPgAdQ6P5fVUbMAMArAbwVwGUUTqobjZjpoq+kqAlNjB0vGvli20UjadrHTOpfmiRntF0ij0j4wrGmhmiSMxvdEdtDZj4nZbYdcSRVItixxGjJ5GURB1Z2tzMtiC0fYq4JKv1GQrP24saiZuwhdXANf2j75R7O03B6be++SyNoAmNfCpkCcNTUknRWJdoB25Zoi2iismhCs0wzDNOCNMwxbTNMtmNqqiw8S6kIGbbXjH022TPMcoZUz388eVksnTvbV4ZFyjTr88wwlrCuNrQfsYibpBmm/3MxlfwdG5rHLQAeVNXPAICI3I0iQOWB0EBV/z/T/v1oWkIovQQRVX0viugYqOq/B/DCjnZvQBFhc7b8fgDP6XXOketjzOeDh161hY6TkyWOjpq+Jlc7TDPVODbcNmqjNMezIGC0H6IVETquLlY4mjfNNFcXSk0zdNO7KnzX+DLEkpyBLGIspHd+hLP7yogp0xkRVBrq4fACtgt7+OVTviRoUV25OaBVUrXyIvYwzUTNNY2EZmGvGTuHYMNvJy+TQ9iYxbq+SoJmxqkyNlmho1Snz9v3mVAfKlK2ic3vEuaVQFQ4GSJ0sPBemtDMmkpi/hkJM0zCNEMTleVmTK3mYIWFiFlUV8D8sGxmymZEKMkM1c1OXhYxw8zlFIvVMYD6Fj+YScI0A1MWWRsT74hKOCECydHRDCfl3lBhbT+Yz6JaEbapXSq6pm/ysvl1h1g+2nTw722aGeYjciAi95v/v6t0lQjcCODz5v8fAvD1kfFejgyFw+Qzq07BNHN01F6kp2CaOSKL+iRMM2xfmQnsNaMsqdoUTDMTChHcBzaZlj4bJiBNwTQzP2yXzXa/zAchxDIF08zJSXsOUzDNnBVCgP5RMwM1IgtVfV6kvk/QyTehEET+09RJd3+HrhOWN33McMt+jqeNKbB2xFmqqUJs5wzJJZb2HTAqUmKusbkAWJKzpZZfWI19VyLq5sa+MpFcBUBtkrGqr8pZVdtlydTtJBomzNXc7tXI9tnONM1Qcw2JkBH7hww75DKTSq8U720tSbjp2F4z1Fl1yy/yZORLgLzIh6R4Tzu1Chl7hBmG5PqQOdOSJM7HkpzFHFMb2pRI1AxL4W7bJnKGhOc/mbysZzr3IV/s1GydWFfpOKswRqqdGXvdHwyb/AjSjZlmuoJRGojIXwPwZgDfVlpPopwvQcRxHMdxHCw2s7P4BwHcXAai/AGKAJbvsQ1E5C8DeAeAv6uqv5sz6IURRLJzfYzICdJ1vtyx6TiRELXUjrzMcTUZvkuyrUpkw7xkbpHwZWjmIuQLrKE+Dl8KqZTyMQfWVNp3lm+EaUdy/UbM142UGg4lKt5GvpF60qY+nh9EZuQeKkOCG/kIiJVAYl9gm9z0LmGyoP4eZNxkGHA0s6oNy+3rn8H8RmxZZDxbn7uZXa5WY9bWnNDnkTijNsdO+H4QzSlzYB2Szp2V9Q3LZQzSkmQ6sKbCgHPqtsEQR9vkmKoLEXk1gHehMD28pQxQeUVZ/yYA/xDAVwD4yTLNQ8rcs7+CSCwJzJCkZOtql90nKmDEzTWx/WfsH9T20Ygwkdp9l+1TIyjNNTP7Ao6YWVgkDWDMOaZsnnlbxppZYYFEzVAHViaUpEwzlaNrewo2uia0azi1hj9QI3dI27zSiHYJQg5JXsaFFNuubbKs0763uw4h6b/RN7U783nqNM0QoYM5oebu35Jphqnzl8SjZqhAw/KIWAEj/M0aZWTXXBo1M2+XEXNOI5JuVdy0qf1ieu++m/CZ5uugrY+16y+wDPk4XMd7I/VuWmdys00JQqp6H4q94mzZm8zxDwH4oT5j7q0g4jjrgDqwOo7j7Dm71sj04UILIixUl7fr77BEVYMk7fCy0a78ab5s2Q6TuVqShkaEqU3LT21rDmC7767K+kZukZA91GpJaKZGUmY34So1IdroE77SrVNrxIGVZWBlX58NLUmpbWg4sFqTS/kzZa452x6oTCWNtPbBodQ60ca0JGaOSrQajbDcWBZVpqFIpX1fFxHtSK6ZhfZhDqzoof1g2UpjWouE5qTShHSaZlh9z4ypiSyq1bNnw7Bjz6M5ppoOljMksftuLIvqSu0aE1+rmPaDroMhkj1zrc3WaCzbfZN9phB9dwYXRCZErqorLNxDBA32ILB49yHkqySlVWY9Eqo+LMlZwzTTvZg0F522wBJSwAvLE8J21wUAPTgzU0BpbpGwqLK8JalImvInjbipd9+lQknKXFNJj2ZelaBiF6du20fDL8S8ZOsU8FboCEKaacde2rE9mLadmTRhXgkko2tYtIsl6neRipCJ9E2ZYVi7lDkn1DP/DeYPQtKwN+dAhJ1YdA2sGZYIGMzk0lgTyDoRSV5m4aaX1PrWGiYbmoOEfgj2F1Ry84Nsaj+ZGFPaOTjFuRdEHCeK777rOM45xDUiWyRbIjVfrLrM88xTotLTEblKqAqRqhXjX6x9077bY+qEiraWZGluDSm1EI0yhFwAtqzULDAzzIpoMoBaw8EibXKxskQs2+qCtOuImgnaEWquMRvXMSdVulNwpTkxfYPphrWD0ZQQLYn941YaApZtNeGsWrHJqBlDfk6RiPajywQS037kplJPOJlSM0xu1AyrZ86lzDHVahlZ6vZYtlVmroF1JLWawKLtEqTMtmM5Q2ja9+bPYpxQxrXGsSW9saZF1tCcfVO6sO+K5trfT/uRG7XSN+tq1phQLPfoI2vvBRFLX/XXkGiXVF+WdjgudOR5gadDdZvt28fFT7qVdyKdO0/nHBYis0izSJqU+aQy3ZgXa5hDrt8Ie782BAxiZkGiPmSob4TahNTtTDhp+4NwQSRTOOnoU10SssgpiZqR1H1NTD2jGBA1kxQ2zpYxX4vO+tyyts9G1PeDCSyJvWaa5ppgmiFCQq4/CN27hkXXtJ9vADxCJmaaGZD2nSUvi0XFWHJDelP+IHUZGW/EO6DrfLls1FyjrhFxnP2BOZ46juPsOe4jsmPGSJpNCbmU5pczWr8OmNNr/kZ4IGVtB1VbzpKcyYBN75h3vVbmmjqBxrwyvZDcIUDlrMrgyb8Iy3Y0TOPWrkwzRJNhYapUYoZpOLXOiKajr2xjr0fjj0+ib2i74CRM6om2QYkJZxNEc4qk9mqh9ZFol9Q4uU6oTMuQ1G5E2tm2M9KHmSQb7cpnlDmmWm0KdVZtm2vss8mf4RHaj0jysnQZMzPHNc1sp9110BUlWb0Pxph9tigcuEZkQmzC/gZwc4wtT+8r0z2v1P4zdV8SFUOSnAHAjAgqdDGJhPTare6DDbnRDu12KynazRrRLuaXoVEBxG+EJT7LjZph+9QEoeRoVq9uycyq5U8jGIRspUlfkpg/CNtRuKs+HFu7f0pQOUMjhHiTxDKqxgSNrnZMqGBtB/iNcDMMEXxyM6um9osJf79U1AwTWGL+IKRsZQXnhv9GW8AIfiC0XSpUlyQvq4WKdlnXGhhNaDYgO2q2aXyDL+5NvYdiuCDiOPvCBOP/HcdxxrIL4WcoF0YQUZKkhqngxjgsAVxdyJLxsHbVsf3YzU5oBlIWV4cGjYmYXCBSfhHNTCjKqnRCtVoN5l0ftCOKdjud2XYEmvbdOLDOQhQLwcoS1dcnaWdNODGNBxDPGdLwKwlaIOLUSpKX9Ymaqa5DV320jFzlcBMk9n5ZOylNR25Zthkmof2I7luU2Gsm11xDTS6JZGNESxLu+9Suuewc4ZmzeULYDrosQqbpcFruvtvQkhR9bNTMkpl/IzlBupxVU6abqoyttUTTsSRakjFmndR7IxzvIneIxTUi54SoSs/eZAOuYsxcw5KSJbMM0odXTJ9mO3tMMyFa4aTy/TCLOU1yFsw1y2g7YWaYlHll2e6ruUJHrP7wqBaC6AvddI7tNUMEjMYmc32FE1ufMtcwgYW9tHel/UkJPrkRMqz9uoSOqC+J9c8QUhYRDIBaOMgUTvjGdcQfpGHWaZtrWATMiggT3LxC1gQWSZcdqgvTjq1f8fWNtVtX5EvVNyJUjB17F7izquPsC+qmGcdxzh8LzyNy/mkmvTFfBdQJqlkHoDK/ZO8/Yz+aI9qPVB6RVG6R8NUzo19Etcf9rMqp0d6TRtRuO1sg1qxjjuucCAkHVpIrJFySxsZ1YRz7tRjq2Yd5Q0uQMt3EEpWxFO/tE1pH3qiWBOgwrxDtB+vLyNWI5LbLNfGM0ojkJzSj2g/aJ6LVGGSuYWXEDEMdU5nDrM0FQrQpEcfU5u7a7URkNmqG7SFT7w1jomuY9oOacGI5Q1gZWmW2vJm8rPiZ2lfmbJ3tmwoqyE1Ati/skwbnXAoiDSedcn8NZcICCdVttKsy5I1Tz8VMPNntTChk3yRn9jglvEgldLQzq9oImWW5UDVCf0EWwyqkj99qlVab1hpWIXa2HTVjQxur38SuKaGemWtmc2AVvhzafh7NORCpMDuzagjFbf/dmXACdIQGh6YxIYXNGYgLBLntuhgjlPT0IaEhtl19cze4y8ysSrOaVu14BlPeh4TlxvaaSWRbDf4gDUEjmF6IoAFwASM81zwsl2RWzQ7VRavM0iWUVGXrWEMz6cqmyt8H7ZDeykeEvoe2Jxy4IOI4+8Jqf9SXjuM4Oaju3lm2D9MXRBTV15rdI6aSLJfmCzI4CW5hWukU720zDEvQY1WNB8xcE3FMzU1y1uzTzi2ysB/D87ZjGt2ltzxeklvImmGk2pOi7ehqx6EOrI2v3citapUkcxJdE+pT5hr7JRrkk4bTIYmaGZXinWhJrKkhtCXaGWW5QHKjZxhMm9KH3MRomdEyjfvhbLvOPCI9tR+5OUNSzqPUlMjyiBgH8IiWBCx5WSJRWcwxtWlSsfWF9oPtIZPaVTc8twtros7MGcJzHuWac9ov1yVZf5cpc82aHV2HwN9hpTalcRGHmYxcI7LH8EQ45mELwlDHlYtvL01MKYkHgWZbJe+vQZlVFZ19hEXDWAGjCtWNZ1YNv/JSDzCPpRwlAQwU7Vjs68kWwzFzTddmgsFh1Zpuwq9lz7eMRc0kMqLmCieWSGZVYX1yBRZ6rpH28UzTTFTASJYRYYH1GSJ0sHbUt4MkE4sJFaacmmFSm9kxs08ww5hFiAkVVRlpB8SzqFqhg0XXxLM755pr2mWNepJF1ZK7z0t8TSbtSTbVPuebCvsyT8AFEWdLRIWQXeJRM47jnEP2yfn2QgsiLMlZLl37EbArGt2RV5S0M2XCTD1Mu9FdZssbZpjy625m+oT6YKIB6q8jIQ6nqRTvMQfW4phoVoJzbOPLsHXqWkMR0YwAxlzDnFWZUytgtBWLdj3JDwKaH2SElgTg+UHOjgfUn5DMxJEbSTPv0DL0ZUiETFWXUIsN2lemp/aDJg6z82KmkrizarYZhmlgSA6Sau8nkqgsd48YgJtugrOqdUxlETLBJNM0r4S6PJNxroOqPc41r7CddhmrEe8A22eK/hhTnFMXF0YQsX+UWlKML5pDPLBj6rtUFlWabXUWs5eavsE2ulQclouvXRBmlTmHRd+YOUYiaZiPyCpxCwVzzkLnmMlp2Tnvtmus9atICO6SCAsWZq4JmpD5Qd1fGydsjxP62BdT6LNu4aQxfzIOS3LGSPmIjPUNOUvKV6T3/jNEMOisX7PQEfP9YHu/2IyozK8kZYaJJDmzmYlXRFjgAkRbOGGmGxYNw5Ocpcwr3QIGE06YX0hznLjQ0ddck7uOj01oFt4vuxYEVqeuEXF2xOEGd1MdQyWETI3lRE1GjuM4I/C9ZnaN+ULMjaRhqdvZx7USibtpiyNfD0z7waT+iDSfn+TMlsVNN5UZxnwF5kbSVGnfSZkSR7cGRiiJRdU0EqiFj9NVZvIy+scz5yjHUWaOKU7eHoelgmdakhBxQyOBMrUklj4akz51OfV9GWWaydR4dNZnaj8CqcgXuvMtcx4t15hGX+LgmjLDkDmw/CA0DTt5jqjGI+GYGitLRcgErcfY3CIseVldR9bLjuibWLtAMzFlaWZJab0HmHOqMbfpt7F0jchOqFRh7cSeFHsTkcSddbuUD0h25Iutk3ZZJNtqw5cklNn3U9nn2lJxMAu+H+36lBmGRdKEBebALK5hUzxmrmkmQwqL5hHmsmjVL+UgdG5hHVyrPmzzvFTyMprQLOyLg0qYUOZQq8QMY802TKNSRdzYix0mnSmcWPoIKrG+Q+r7si4fEdY+JajEhI5E+GtcMGAJxoxQyzams31iZhjiN2LNMDzpWFs4CZExyX1jyDHvU8+hFjDqqeZHyHQLC51ZoCMm8ZQvCSsbFXGTaV4Z4l/S13QzO5z1MLeoa0Sc3XEwUdNMEEImh0fNOI6zB/Ty+dDd+6j04UIKItUfiGo32mYdVr8yydXsOGwHR67pKH8STQfdT2FGvjLM13Ht1GqdUdE65o5kdRmLpGF9l7lqpwQzFGYaG1WzInvVnG1f/I/RagRYrhD21V/lDqmLpKFOL7UkykwpZpzwZcva2TIh84o5ulr6aExadefIWZVpOlIOp6wd25OGmmbmrXbC2sWSnAE80oZFw8xYvo7u/CAsXXtDS0I0JykH17ospdUofqYjZBApa6+HXeeLro3MbJ2KwiEv6Xrt5s9TbZaPCwRTEQBcI7JGIpFXRb252LElMGn3y/SiDpyeLHB4dNCqn5M+MVWjlWfi4WhtVeOsEfpb/PYny1WlFWnab1dlH7vINfsW7YpCs7WNMeHUZWGcXL+RpR5jVplmrK25LXTUJhzr7d9qVgsl1lwTi66xt3sQEmYHqEJ0V20zjBhhImq6YRE3TDhp9I1E4Vj6CCqxvjG2lNAs2i5X0GjUJ6JcWF8qTBBzTfD9mCXMNZVvCmnXqI9Hw8SECebbkRRO2MZ0uT4iCZNLrrkmrDvJLNBkfcsN313mmnUSJpzA4dEBTk8WnfWp/oFcQSC3XT/TDMY/01tk8oII+7iZAkEImRpTNc3MpmqaWU10Xo7j7IQghEyNvuG4U9HM5DDNt+kaCX8MIdr+lLNqqJ+Tsi5JmOUoiasJpV1mThjVkhhzzWIZtCS2b9tMk1SRlofNe74774p1YB3CvIygSWo/pLvOJlWrtPPM98N+pVZ7zXTsuBsz3dgdcmkkTaaWJFxPFoXTmHcPjUmrb+bilZsKvvM8mRqRWLtcjUfXeLR/T+0HM7k0nEzbJpX0XjMhv4nRfkhb+xE0IksS5dI0n8S0JG0zDCuz5cwx1ZpcglbjlJph4ibcqEl42R4PiGs/LLkajlS0TOy8qbFz6oDdCAWeR2SNdJlmYhEy9o8u5DdkqjBlvh0lXaaZg3KXumwzjCmbBZ+ORkBErhmGjFfWXz1RHB221au1Tde8tKvF0MyBmGFC35TfSCzx2UKPK6EjlQQtJnRYQSvU2+iaEPJr85ZUAkSXYyoTJoKmpCGFtgWMSuAZIpxUgyQEkT6CylnY+VLkqnRzzTGW2FzHCBqd/UcIHZUAwfxLbPbTiH+JKWd+UCwahoW/p4WT7uRkMyyr+nSisraAsSTPeixCxq4xPEKm+HkwF5yUL0u2cV2+uQatvmlfknZZeC8cHMyxWCxb9Wf7AlzAiJlaUgJJrL6vacZ9RNbIVE0zQQiZGkEImRrzqSY086gZ55yzjDh/75KTiX6xByFkavTWcHgekQli/yiZL+vgRd3IeUUiblieEVqWincPmvpEMjSaW4Smgjd9iKZjUWlbUmnf88YLGpbUwjc3kS+VSpmotpj2g9XbuhnaD5+U9pWGn0q4BTSh+m9EvjAn1BFakqoyYXqJpZu3xBxi+xC+5rv65pphcvumNB3RssTfbx3aj2S+kUheEtQaQLs3zKqKaGk7l6bMMMzkwvONdEfFFGPmOaZGU66T9SRXc0LXuTPHZ8toPpJcE00kUiZVFivvZMeCgGtE1kgqaqbRtueFpxlRmRnm2hIHZ4SXk2uLqmxu/UFifiWz9sNmhQmQJGeVuYY8tCerJQ4Oyhd5WXbl2hLHh+WC1vAXKVWW9qVdmWuG3LDdfiOAVE6zq2pxnVdaESaozIx5JSaUcB8RO4cg5Jjka6Ev5pVQIsYnohIc2AZ3ybDc4cKJJVtQqcoyBZZGnx0tjCnBJSVgxNqx0FrWZ4zQQduRjfAAsyFdXMBYkQgZFg0TM8M0w3yJcBI2oZRV3aexcZ02fhbHbX8Q3q5dFttDphkhU/w8mM9w7XTZaFfUd695uXvIdAk5Z8sODmdYlJqGEJZry9g7hY03ZKfb2PtKDufQ06Zmprdpxp1V18cUTDNnhZCusm0ThBBLEEJ2CYvcmYJpRtBWuVKn1i2jE5iDMx4doynaIFM1zVw73b0JZEFe7Kxs25wVQoCephl1jcikiDmm2g9vrskgETCZEnlzD4P2i7lyQm1oScqfxOTCHFOZ5mRh1IGzmf1iCj9ZJA1TzZrfb4QDK7vFrFBS5TBomFcKzUSuucZmba3zKTATjfl6RsirQpxagVp7wEw3a9KSVM1ZXwMVVHI1J436Hb0oU18TmeaaqMaja5x1az/YfjCRCBiAR7RUphS73UFPM0xzp912zhC6I++aHVPTa0e7XVijUrvcphxYWfKyWDp3tq9MrlmnKO+OtLH0dWDdFK4RcRzHcRxnZzCtylTJEkRE5PEA3gzgOSi+XX8QwKcA/J8Ang7g9wD8V6r6J2X71wN4OYosDD+iqu8qy58L4K0ALgO4D8BrVPt4geTBsq3SsgEOpSxz6mxAn6ZvSLPeZludEcm9ckwlZfaLYLEwczxsb3BXazXMHIj2I+bAmtIWBuuR1UbQ9PD2oznMK+E3IsQ/RUqNwlzbfa1pJubUCtQ+JA1txTq0JJbKLyTurMpzomRqThjbyriYGd6b1HREy0iG01zNyQjtRzMMttsHBKg1HA1HUrJJHcuOGvMHWWVrTtq+HfY4nR+k+/lPb2ZXnsusRVUobpezaiR3Esuiaumr6VByjlSfMVqQbWpGNnUuEbkVwBtRvNrerKp3nqmXsv7FAB4F8P2q+uHYmLkakTcC+BVV/S9F5AjAdQD+AYD3qOqdIvI6AK8D8FoReRaA2wE8G8BTAfyaiDxTVZcAfgrAHQDej0IQuRXAL2fOISlMsKRlucScTFm7rh15q/0KrMmFlTFhojK5CClrP8gztiMvMesAzeRBcWJOqG2smwoz/1QLn7GpzHOFkoS5pspLYhbfKrmZfT+VL2jm1LrSI8ykXR/GWbtwYkk4sFJBpaojAgvry8hNcjaWmM9E0oE1shNvV9/eKd6HCx3M9LIiQkWjHkToyDThMDNMyoF1Wd6TTYdSmGNiNsl0TM12al22zxHL6wH0yBnC2jEhhuSIiu8102WayTfj9GFTppxNmGZEZA7gJwB8K4CHAHxQRO5V1QdMs28DcHP57+tRvPe/PjZu8m0jIo8D8I0AfhoAVPVEVf8UwG0A3lY2exuA7yqPbwNwt6peU9XPAngQwC0icgOAx6nq+0otyNtNH8fZCbNtvZQdx3G2yWrV/1+aWwA8qKqfUdUTAHejeOdbbgPwdi14P4DHl+//TnI0Iv8BgD8C8DMi8rUAPgTgNQCeoqoPA4CqPiwiTy7b34hC4xF4qCw7LY/Plm8W4nCanfmO5ASxzqX1KUwoKPl6DU5SMadVO3ZubhEaDmzKqINrQoPR14GVh/7W1+Og/BJttCPaEaYZUWZ6SZhUQpbVFewXaRmy28hgGn6Yvx0x3axdS2KJheza+iFhudFdd43eb92W0dwwt9zcIay+q13UNLM57UfM9GLLGynZM9O51+1YmC9LD9/WRtpHz2Y9PSUajHU7pmY7nibCd1nOEJZF9ewYAA+trbUpxMSZ0Jwks6OyNXEHG9BtyFn1RgCfN///ENraDtbmRgAPdw2aI4gcAPiPAfywqv6WiLwRhRmmC7YSaaS8PYDIHShMOABw7UP/7ys/njFPB3gSgC/tehJ7hF+vfPxa9cOvVz/8euXz1akGq9XJlz78/v/mugFjr0TkfvP/d6nqXeb/c97j2e/6QI4g8hCAh1T1t8r//3kUgsgfisgNpTbkBgBfNO2fZvrfBOALZflNpLw94+IXvwsAROR+VX1exjwvPH6t+uHXKx+/Vv3w69UPv175nBEUKKr6lRs6fdf7vW+bBkkfEVX9dwA+LyJBCnshgAcA3AvgZWXZywC8szy+F8DtInIsIs9A4bDygdKM82UReX7pVftS08dxHMdxnGnzQQA3i8gzysCV21G88y33AnipFDwfwJ8FN44ucqNmfhjAz5Yn/gyAH0AhxNwjIi8H8DkALwEAVf2EiNyDQlhZAHhVGTEDAK9EHb77y+gRMeM4juM4zu5Q1YWIvBrAu1A4m72lfOe/oqx/E4qI2BejCFR5FIW8EEU2kMZjrYjIHWdsVE4Hfq364dcrH79W/fDr1Q+/Xvmcx2s1eUHEcRzHcZzzyzR3aXIcx3Ec50IwSUFERL5WRN4nIh8TkV8sk6qFuteLyIMi8ikRedEu5zklROSHy2vyCRH5X025Xy+DiPwjEfm3IvJREflVEXmqqfNrRRCRW8tr8mCZRdkpEZFLIvIBEfnt8tn7X8ryJ4rIu0Xk0+XPJ+x6rlNBRB4vIj8vIr8jIp8Ukb/h14sjIq8RkY+X99bfK8vO3bWapGlGRD4I4H9Q1V8XkR8E8AxV/Z/L9PE/hyK721MB/BqAZxpn2AuJiHwTgB8F8O2qek1EnqyqX/Tr1UZEHqeqf14e/wiAZ6nqK/xaccqUzr8Lk9IZwHefSel8YSkjAK9X1b8QkUMA/wZFwsf/HMAfmy0wnqCqr93lXKeCiLwNwG+q6pvPbBni18sgIs9Bkbn0FgAnAH4FRcDHf41zdq0mqRFBkbDlN8rjdwP4L8pjmj5+B/ObGq8EcKeqXgMAVQ05Xfx6nSEIISXXo06049eKk5PS+cJSprH+i/J/D8t/iu4tMC40A7YMuch8DYD3q+qjqroA8OsA/jbO4bWaqiDycQDfWR6/BHVylK7UsRedZwL4myLyWyLy6yLydWW5Xy+CiLxBRD4P4HsB/MOy2K8Vx69LAhGZi8hHUSR1fHeZ/LGxBQaAJ0eGuEjYLUM+IiJvFpHr4deL8XEA3ygiXyEi16EIiX0azuG12pkgIiK/Vtq+zv67DcAPAniViHwIwGNRqKWAAaljzwuJ63UA4AkAng/gf0SR30VwQa9X4lpBVX9UVZ8G4GcBvDp0I0Od+2uVgV+XBKq6VNX/CEUGyVtKlbrDCVuG/JSq/nUAjyC+ZciFRVU/CeCfoLAK/AqA3wbKzbTOGbkJzdaOqn5LosnfAgAReSaAby/LeqeOPS/ErpeIvBLAO8pdjT8gxe5uT8IFvV4Z91bgXwL4vwH8GC7otcrAr0smqvqnIvJeALeiewuMi07fLUMuNKr60yjNWCLyj1Fcv3N3rSZpmpFyJ18RmQH4nwC8qayi6eN3M8tJ8QsAvhmoBLcjFBtI+fU6g4jcbP73OwH8Tnns14qTk9L5wiIiXykijy+PLwP4FhT3VNcWGBeaAVuGXGjMu/Avo3CA/jmcw2u1M41Igu8WkVeVx+8A8DNAMn38ReYtAN4iIh9HYcZ6Wakd8evV5s5yEVwB+H0AITWxXytCV0rnHU9rStwA4G1ldNEMwD2q+ksi8j6QLTAcAD22DHHwr0TkKwCcoliT/kRE7sQ5u1aTDN91HMdxHOdiMEnTjOM4juM4FwMXRBzHcRzH2RkuiDiO4ziOszNcEHEcx3EcZ2e4IOI4juM4zs5wQcRxHMdxnJ3hgojjOI7jODvDBRHHcRzHcXbG/w+L7aUHxuMaIgAAAABJRU5ErkJggg==\n", + "text/plain": [ + "<Figure size 566.928x283.464 with 3 Axes>" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "split_ax = True # switch between splitting axes (using pyicon) and not splitting axes (without pyicon)\n", + "\n", + "# --- generate axes\n", + "hca, hcb = pyic.arrange_axes(1, 1, asp=0.5, fig_size_fac=2, plot_cb=True, \n", + " sharex=True, sharey=True,\n", + " )\n", + "ii=-1\n", + "\n", + "# --- do plot\n", + "ii+=1; ax=hca[ii]; cax=hcb[ii]\n", + "if split_ax:\n", + " ax = pyic.split_axes_vertically(ax)\n", + " hca[ii] = ax # exchange old with new axes object\n", + " for axsp in ax.axs:\n", + " pyic.shade(lat, depth, data, ax=axsp, cax=cax, clim=[0,1])\n", + "else:\n", + " pyic.shade(lat, depth, data, ax=ax, cax=cax, clim=[0,1])\n", + "\n", + "for ax in hca:\n", + " ax.set_facecolor('0.7')\n", + " if split_ax:\n", + " ax.set_ylim([6000, 800], [800, 0])\n", + " ax.set_yticks([6000, 4000, 2000, 800], [800, 400, 0])\n", + " else:\n", + " ax.set_ylim(6000,0)\n", + " ax.set_yticks([6000,4000,2000])\n", + " ax.set_xticks(np.arange(-90,100,30))" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "6759f5c2", + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.8.2" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +}