{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# import libraries\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import micom2reg as mr \n",
    "from netCDF4 import Dataset"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "# read gx1 sample data (3d temperature)\n",
    "nc = Dataset('datasample_gx1.nc')\n",
    "T_gx1 = nc.variables['templvl'][:]\n",
    "depth = nc.variables['depth'][:]\n",
    "nc.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "# interpolate data to 1x1 deg lonlat grid\n",
    "T, lon, lat = mr.micom2reg(T_gx1,'gridsample_gx1.nc',1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.QuadMesh at 0x3245c8a58>"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOy9ebxlR1ku/LxVa+29z9Td6QxkIKGTG8IMwg2DYSYmGCeuXCIB5IuIRsQLGOCLDAIRwxQgDRdFjXARFQVR7wXUAGEQ+C4aCWEQgQQIIQFRkSlJ9zl7r7Xq/f54q2pV1Vp7OOfs031OZz2/X/fZe+011Jqeeut5hyJmRocOHTp0ODKhDncDOnTo0KHD1qEj+Q4dOnQ4gtGRfIcOHTocwehIvkOHDh2OYHQk36FDhw5HMLLD3YAQPerzAEuHuxkdOnTosKNwG77/n8x8bNtv24rkB1jCg+nsw92MDh06dNhR+BD/5TfG/dbJNR06dOhwBKMj+Q4dOnQ4gtGRfIcOHTocwehIvkOHDh2OYHQk36FDhw5HMLZVdE2HDh06HAm42rzbfz5HnX8YW9KRfIcOHbYRrjbvjkgxJEtg6wkzPH567M3sM8WhJP65kDwRXQzglwAwgH8G8DQAJwB4J4C9AK4D8FRmHs3jeOMbYtUnNvXnNqS/s9n4sdr23aFDh6kYR6LzINfN7GNe5D7rMdJObd4dwKY1eSI6CcCzAZzJzPcGoAFcAOA1APYz810BfB/A0zd7rA4dOnQ4krEVFv685JoMwAIRFQAWAXwbwGMAPNn+/nYAlwL4val7IgXSOvhKMzWADfv13ed0e7+ck2WpVT6LNb7DLPbDPWTs0GGzFvK05/VQWOBbganWe6hQbACbJnlm/hYRvQ7AzQBWAXwQwKcB/ICZS7vaNwGc1LY9EV0E4CIAGGARpDXU8hLADBCB8gwwDGYDmHgWq/d/7w8ntu0cdT66ea8EHaF3OFKxU8kdmPG9HEfukeQ8fvNNkzwRHQXgcQBOBfADAO8GcF7Lqq3NYOYrAVwJALuzY1nv3gUsL8oJKBJiNxVoOAKMAa+ugQvpOx47eAq4qgAAHyze2RHZJrCdogE6HJlIn6v1kPOkZ/Icdf6OJvqpGOP/CxUPTDDy5xEn/2MAvs7M32HmAsBfAzgLwB4icp3InQH869Q9EQFaAaMCGI2AsgQyDSwMgD27gJVl0O5doIWBdAABHrvw1DmcSgdgZ1tGHTaPQ3n/z1HnTyXwO7TRYQmeFMk/rUXtyDNQnoF6uSgfEzAPTf5mAA8hokWIXHM2gGsBfBTAEyARNhcCeM/UPRkDHo5AWQb0e34xKyXdUd4HFnqglSVQUQLDEbC6Bh4OxbonteO0cmBrPOqbRafh37GQ3u9xz+S0Z7UzEOYP71ckJZ8jmYbBo8lBi/PQ5K8hor+EhEmWAD4DkV/+FsA7iegyu+ytmz3WkYqOPDscbmxGSkn3czgMhFn2f0ftgIh5+7gmd+fH8lnH/JxY8UoD/VyseE0i5RCBiSKRiSoGqgpUGvB3vgs+sAozGrVb9DZyx/WMbBhg4/8eqZjnw911SB22AodyNLudyH6mqBrE1ny0zPouP3jwTz7NzGe27WZ7ZbwyJKrGoTIgIllEgQYfED7nAIychsqPg1odQd16G8yBg6LpA4DWoF4P1Mut7q/lb1UBZSVyT1VJFE9Vyec5kP7hkGG20wPcYedgO0qGW4V5Z7Ruth0NjCH3iOhdUMoM2F4krwjIsojQxXJXYtFnCqwDix4ACIACiBnck21JKailJaAoAGPJ2q9P4sx1RE8EKivQaAQMR0L4o5GQfXgRN0D6hyIF+1BjPUPxLmJn+yO8RxspIdDWOUx6Lsfp/NO2m7U9OwWTCN5Z6KS1cCKR/axiYxcAFIkP8+D4Y20vkicC8kwIWGtwpoFcgzMhec4UQBCiT2GE+ClToFwDFYPcqIBZ/oXLHJilMwCE2IvSJ0iRDtcLkqwmEL57UDf6QE4ixsNteYzDLO3ayHW5I1mXhwPT7ttGnrf1bLPe/d+RngdSEmkYkntI+rDGbyvxJ9heJN+hQ4dDhtBJOi8JY5Klvt0IelIhtCMJ28vx2r8Tn3XCk+sFSgOZBg9ycD+HyTU4c3JNvRrZU6CKAWOtdWbAAGRMbcGH5+os+7IChkOgMrKsLMFlKd8BwGn09u9WOmmP5Adtq7DdiONIwbwSlcL9zSLrzfOY68Ghzg1ohZVqSGvR3bUCOUtdqdqyz6xtTrVV//5b3riDHK+G60QnUwElQEMCa5FuQARWJFp8sJlz2qrSgEohdocPfOa3xh7yvLteAqwNY2Inkpj9sjwk5A50BL9R7FQdd6e2uw2zyijjLPpxFvV69PvNXK/t9O55gvcRgEaIvqokUdSR/hSJJsT2InlATsSl67oomp6WMEpFMNb5ytYBTQyArRXPhA/+40vXdbirvnI5Hrt8obX8g+gaYEtIvXNGHlpMqk++WYfhetuxFevPyzreimdxUinhacdbb3s2Iwmdm1+ADxbvXPd2s2LdbVItJO6+K1esjOtlU9SY7UXyWoGXF6JhCGuC6eeoFjJUfYmZZxJS1yMDNRIi/tDHXrThw37g9rc3lm3kYdmqF7nD5jCPuuXbVVfeyoktthJt13OWNsxVRw/DEueEc9T5UT4OzxjmGMFwlAvkrXkjOUHIMi/fuCCVSdhWJM+ZQnHMcr2AAJMrmIzAmcg0VDHUyAjBDyuoUbUlbZlm8W2XWNsOhxbdSKwdG30P5v3+TL0/LTHoW2nFtx7fKQRBXRrf4WhtizLGkXxsa9Z4aAX0ejCLffBgMo13E3l36NBhW+FQd57nHfuMueznHP3Euexn3thelrwilMsZWNl5PVwUjRtRsXwkAyljoAkf+NTLtrxdbdZGZ8HfMXE45x+9I2DWazhrbobbH9msd7WyDPR6YgmXFbA23Hg7g2gYWlyU42Q5KM+g9uwGAPBwKMlKi4uSA/SDH4LXhjCra35foWREvR4oy2T+jKK0kX2mlmoqI6ydZTCLfZilHqr+DpJrQECViw7vyd2FRzKgSha5pjBgrfDhj75w7K52SiJRh52BNvLpSP3wYL3aPQBQvw914vEoj1nxkSvqwAjqNkkVPftRr8SH/366X8+Rexv03U/H6PgVFLsyVD0JDlElwEoMU1UygGOQHRSJ+aMf+I32YzzsFTAZQa+V0N87AP7e96X8ei+Xul4LA5jlBSH4BQ2TTRZkthXJs3O0ZpDMVrJWu4HEv1udqlzU+NjfXtK6j3N7TwaXxaFrdIcOM6Dz4TSx0U5yvc5X0loI/rhdKBezOjJvMYM6auD178ec/Wp85MMvmKkNlOVQgz5o1wrMsUdh9cQljHZrGC3h3e4YVU8MVKoYrMgSveAhP38FAOAf//S5eOAvXgEmUTPyfQuyPQOqWkL/+0ej9901UFEJPw5yT+7iq5zc1m1F8uJolRMFM8gAesRyYey1GUfuAPDju3+xdXn3YnXYKLarzrqTMc8RUBvhh8tIa6iFBYxO3I1qoMG58qQoCZUKqgK0ta4f+ZOXt3LMY5cvBCmJmtEnHI/q+L0YHtVHsaRR9S3REoEkktuTdK1GECoFVD1CuRCz8n/95f2Rc7RYkJEGGYZhQnl8hoPHLUMXDCqjTX0i6CRsK5JnK9foEUMVQu7KBs987H3/70z7UMtL4IMyBPvA2ju2qqkd7ojYxKQ0naFx6OF0eNq9C6an48KHqs6aNwTwsgaV0xlTH3M0hnc9HmvH5GCNKPMeLNa6scsDV6IPcQl9jOHnKrcb2E6BjNTocomeZOyc15n97ir2zpAUNReSJ6I9AN4C4N72nH4RwPUA3gVgH4CbAPwcM39/Hsfr0GESNuOPOTe/wH+mLJ9bmzocgolD5jTqetR5l+PvrxqvGOw0zMuSfyOA9zPzE4ioB2ARwIsAfJiZX01ELwDwAgDtngYHO8zRI3GwArNb8AAApUC9HJTv2thZdNjRGEcibXp46EC7unoXzu09uTUxRg368sGWtTBFeUgnmHns8oVijS4MAABXfetNh+zY88Q8CH7maCbnGCUlkS3Liz5iz/2r13XLSKaaGArvPPrc10CviTbyoY+/WFZdGID37kGxouOMe9h9aEKVy1+q4lFBHTHo1qfYokfwOYkodHIPM4Oqus1Gkx1NTLbmN12gjIh2AfgcgNM42BkRXQ/gUcz8bSI6AcDfM/PdJu1r6ZiT+d7nXYxsTbJZAYBKnqjDhzjvTr8qwxfbDHPgYF2PBsAHVv9kfSfXYcdgQyTiklFskkk4005UwzsArw3BZVGXvlgnptassdZotvcoYGkB5j+/F5E87xUDRsp4sIQBQspzbGfMy4qfeT4De2+VDWdcu+9dYFzknkKj/hUAL4uoQv45kterJdTN/wYsLmB0l2OwelwfrG1ACADWIreYHDCZjdwpOOoApE3uO3mi9k7asOCigfdBxufU1l6J2rn2bc/b0gJlpwH4DoC3EdH9AHwawHMA3ImZvw0AluiPa9uYiC4CcBEA9BaPgqoYrKSX0sPYYjrnoZeBAtJu1Knp9WqNaoMvYYftga2oMpg65BrZjzrWbaWekX3eXE3vhQGozGTieDbrLl43tRRvEp6njtkL5Dm4l09MXXzs/V82sRDfHQ72fpjRCPy9H0AP7wyT1XTnal41rHpNqDSh6gth978vcfR8wrEwgwzVQAOUELF1sJpMIgOJIbNdc0zwTNZBq5qWvW9CuN+0I7DHIgOoAuKIHdchBJgHyWcAHgDgWXZS7zdCpJmZwMxXQib+xvLRJ4ufQQNsNlBTYmXJzxxFNv7VhVMe0tTlDtsGkeVnk1fCeiI+VTwoCkW2HoifpAGIRojIM9CC3YFhP5vYPKQcUgTq9yVZZ2EgE+XYDobDuY7dBrZk9nbHPGr/TKtBNE6W46qCvn2Iqq8htELCxAwoV/k2IFZn8ZeLCqY3qDsEVzfLQKL/2CZw9ghln1D1IARuAA34kEhP1o7cVULaXP+lqm5HtC5DooCGDD2CzReSbcwUFp8HyX8TwDeZ+Rr7/S8hJP/vRHRCINf8x0YP8Kgffw1UxZic19XhSECXYLRxnHPWZQCAqz/5m4e5JVuDcJKTDrNj0yTPzP9GRLcQ0d2Y+XoAZwP4ov13IYBX27/vmbozA6gR++n9TE5Qo9pKMbkCcoVySZodxrSec9ZloOWBDGkNoNcKiWtdWMD7b33bZk+zwyHAPAm+jQxq7d1I5IxWoNxO7m4tdrKTNfiJGsJyrol04y3/Xg4aFaBkQniwmSrl6N27pXZKppGddCL49gOgxYWxVjwrJc+4PXbowHPW6LkPeTmAFjnzCMe48iPnqPOh/vU/ked3QrGrL7lPmmyCpQ1PDCx5KBKLX4n1DiOWs14zYBINnBVQDRSKRYVykVAOgrBJgsTOO6kllGVCgcIEVryp/wFinRtt20Aiz2QHGZl1DP/Dnz3P7+ah579+4nWZV3TNswC8w0bW3AjgadI8/AURPR3AzQBme4NJPNPuZGE91awIpAgf/sgLcc5DL8NoTw8A8JhzXi2rATADeYnVqEK1dwkfvKEj952CrSb4c3t2xjFLuJRnQJZZacaFSlBM4Kk+r4L1gEDi0UAvF70ekM6gqoCyApcleHW1MTG8XloEmyb50+4V8S1l9dye9aT19V9P7lQn4rh20kg+/tgjXuEjQ7YLDke9n6vNu3FufoEoASffCeVKDyav769LvgSEdNXIgMo6IUmkFJlpzoG1QrmgUCwSigWRmOUHIerSBmU1ROegTEtE8FGUjt1fqMlrSaTSNpn/4f/ttfjE/5HIQ1VMlgjnQvLM/FkAbZ7ds9e7LzKSBOULlCkheD0Useqcsy4DCOj9UM62WpCrawbavwzsXtoORw42mIh0bu/Jsl2wvVjxVm93RB2Se0j64XGdJu+WA7FmnwFR58AyvST1JN7e3Ho74Apa9fuy3HY2CCZpZiJvrXsdHvDLfXSGs+7JHRcAKOQinPuQl29ri349pL8ZqYarCtV/fhdqOEL/6KPASwOYQQbT035aURde6coQUBh5GBJ8Rp7gy0XA9OwPVD9q4XbO+eoMV0fuHuSjJGsdPki0IiBwuHIjIIWmvBbbao7XlT135vs/4jn2ZOueVI0YelXCmdTIyLytfhrAuHbDdrNcOsyGqRbdOkg+JANP8m43NgsSvVykGUX1rPeO6LVCGEoZSTbp+xJa+6HFrZXEXjNLcakw2stJQpkWSSbdhyP6TMWEn1jvbDspvw6kA1ClARjQq3UNp+1M9CkmRSDNZd+koHo90KAPOuE4lHuXpBaMppqACbGEE8D0FEbLCsM9CqNlIXkfBunsABv14soctMkyPjImWM87Z1HvT9nH93NvuHjseRHRzpnjlSpG1Vee3JWNHJilQlyHnYtxL7F/4TcatZIQvLOaKSTxEC5WXseWvRslkrFlX8cZR1bfd+TNRDI3sf3NSy8qWD/4G5F6aMT43+Gtdm/ZByn6AFD1Ncgwsttr/9ZOwlbOwhXum6oK+rvfh84UWC+C+xouEclfU11fW7KROMbdzsr+M2J5mxzxukAcPWPibUKSl41QSzSwlvsMBD8N24vkO3SYA1wUxjn6iXOf3m2n4khL1b+j4oFPk8qVrhP5pz967tRtthXJ3+2ux6NY1tBDI1ln1oofV3e5w5GPWSZ9boN3tALwU6s5R6uLpiGqE6BCWcb/pupoFucEZQUYE6Wte1nGfoaz4sO49sBC9+uFaHGi1jHWsRwTF7gKrPhgl6wIo72D2vIH8IifeS0+/t51lAk5gqF6PUBrcFGCbj2ATClg9wCmp2pHrL0fRgdVJg17KVmPGPlBQmWAkq0133NOU/t8GKoteKvXk6q/A4jloeBecpDodL9f3w+qgI1UU9pWJP+lm/4d97k7IVsT58ZHZyxn0OHIQlpXpgHn2UolnHDOTMBPpuzkF3K6OzBZh9canGkJs1MqImvR5VV8aOcuCCWYgMxFukHsP2qRX9z6bj8mcrq67eCdvW3kHqLM6zZPc85tR2y0g5+G8874DdDyknT4CwP5ZwzUaglWubCiIgmF5brDNTkA1KURTBBRA9SRMSZnQFuR3Ui9GVIEJVG1cvsI4szleB9em7f6vCHpEJzjt1wkUAVfl/7+z9w/9d5uK5IH5MKt7dH4p7dPH4Z0uAMimvxY7BrSup3wA7IPI2hcRms04737XVtS11r+BiRfhzLaY7jQOvc1tc7tuiZXiHR0oGGZ14QdxFfr+nihRZ6mw3NyWDKI9wnyhuVZT5SY6k++63nYCVivw7XRKZCtAX/0UTKrEgD0eqCVZfHPKAXONTi3nTqH5rT4CJW9F6ykbIFkxApXsZaYdud4ZcXgjMG59dsEljyXgCoJJrDm04xX9zh5R65tisT2w5c0AGwJhRk6721F8vfYd6cdRe6Pvf/LxBEHAKX8ff+/vOIwtugIgiXtNkvOz26vtcS793I7D6aJwxydlW7DFr2FrSSyxb3kXppxpK6UONz83zAZieqhNTDeAesbSzY8r8WBCnjLPK5xAm/1t5E/gGBE0HT2UVi2iYM/QWfwo08Wsg+Tau4ouOqG1+C8u73Aj9bYdvhtnbTMSgc7Kx1JfZpcZnwyGcBZ7XQ1GeromNyG6xoCV+RHckYxqCQfhUv2kXUE74uehWGY2pK8tepZyfGVgeREHILaNR06HBoEEydDa1CWSXXGXs9OcmwnOk71dpVE0pAjeZtbEWaTapFX2G7jLDgJV6yjLigZZkfEj2AZrOwSzFscEjeAeHlkyc9gwYe85EheB4RvOFnfyjfbJ3J6rnCW/7n5BQh9MRj05a/FVddLEuVjH/hbtaSm7IinZMBa20bDyir1hWZlM1L7sIXMAM5YiD5nIDeg3GnyTsN3uReyK0UAB1Z+W7imn2EqA5SL6CmtNV8CRtVRN3TleOO4I/k548fv9xKf/v7+f+6s+vWiq13TYavxqPMu9x1db/rq2xr3euH+qet0JL9BPPInLke2W3KXpQa1ARUVqOhKHG8UY2f2cRYZKVCe+Xh37zjr9aRjLUsfFeGlGCDOJHVww/RMRTo5Kzc5MknmtLPerQXPtl74WEvYWc4qsMjdtoHTtSG9pLq7G/aH64VWe2Tx1dE9UZo8W2krsBJra14WPfBpV+BTb9s5EumsoF5PJgzJM1Ceg3ctoTh6EeVAR/dutHeA7EAJkxE4dx5viNMzU/GEHkrkmaonVnw5AKoBg3sMkzGQMZAzoA1I22grkAzy3P0pZRhoyE4sYiD3Y0rJYOPCajKx/E0PcYTOBNyhSf7+v7bfl+n83BvXn2xQLmW1c68w0GsV9FoJGnZEv16cd/Jz5IMjdMAX94okml5PImWc43RhANPPrW9EOl0nuzSIPal2wVpZaSZYh1xEjJNp4KMponDFENT8EmnpqibtNuklLWSVdg6N6Jnku8hHbmf25WdXc6W5bqMe+hGIDxz4YwDAeSf8GrC0CFYKamSgMiUTeyiSckQlo1zM6nvDLMnVJA5Wk4lkN1qRf1UfKJdEpqkGDNMzQuy5ATRDZQakLMHbcukMkprGtpOgksA6dPCy7aj91xjp/Vb2vmaucNnka7GjSX4zJH2v39gfDdXu/lvJsMde6C9f2r5fzgjFkoTGUQlkBJnwpFAgvQPj1bYBMjtBBvIM5jvflXIAkDBI6vXEene1XpQCBn1UuxckegVICLSFjet3yjtT/UxBurbQQ7J1RJ9a30BigbfBtsEkxabqNrrPhDCaoq164dTOJYzSUJb00xop4Tb273W/t/FMyu2Ghzzl9f4+XvPHz8WDLrwC+tGnQQ8Zg+8M4YofIqO6ZrzzsRiGsqGJbBimp1D1CKNlWbdcBMpFIfdiRcidewzqGZA2UJmB1kb4gBimUhL1yACDQfYeEwOsGVQFRE+Aq28fkn0jaiq8n3UYfjj5XSt2NMlvFsOj6hdwvWBNGO5SqAaAHgHL3zIoFzQ+cvXM86V0CMB7d4lU0tcoV3rA6cfhI1e/QOY4DerNwEo0PMhhFnsoVnq1pR1Y5I30crSTpid7Z2kjiXBIHKJRjZKUkEMkckpbydlw1p+238Z3Hm0XEFEhqzDuOlyn0Y47CNaO6YsTmghGA6ZHfo5UAUGVdo5XZlR9QjkgFEs2mmZg5ZlFhlmsgL6Bzish98xAEYOI61QKIxeZw5tIYs2DpfAiVWLhsyN6hv0eIDRMfEdOdVCXq1g5ATuC5E9/raTyhg/n154X64hnXLYfN/zmbFbJPV66f9JMahHu/rL9+PJvHTnWznbE2Y9+1YYy+Tp0ACT2P50d6Uf+x/4d71SdF7YlyZ/xiv0SPqYgPV7LOqe94Qrc+LvPxd0vrSWbU98ksb9ff9b02N9qAY1wpUgTC7rTe75YpJwvvqIm+0/87+fjPs/dj+FeoFpg3HYXaUTXKUzHuT/62zC5QrWQoRoo5ADWTlwB4Jxb9Y357hPvB1UCu25ag751CB5kKJdzmEzB5ISqV8eupxMiyw6T76mF7FPXA8eo+ykcNjutPrUOAmvcr5uiTaIJ49zD/bRY+q1IpJqoTnlY3TBcxx0rzNbc4XjgL17hDQSTk72PctajXSTWMkKfiFi+YcSiHoofQxUS7lj1COWC6O/FUuBsdRr8oILuV8jyCpk20MqAgiGTMQolwQ7H0KK71H8ZALSdSrBVn4nxtWcLt51++X7/u5lisc6N5IlIA7gWwLeY+aeI6FQA7wSwF8B1AJ7KzKON7Dt68IMX+dQ3vR5ZX0mGWfDQnvqm108kepPV+2iQgHs5wtrPYzDcC5TLDCiGAZDfKjs98+lX4Nq3HnkRC/PAOWdd1uCt1WN7GO2SJ95NijBaIfzI/9jvH9DbTx4AGACIyTiUOhrRJ4lkkU6q3JBM0k7CJamExJw03hOHO2ZiOETtmbQ8JPUZdPhUBQgzJn2Z2yA3zB/DOv+mzQu6E3DWz70eWBaHqNFk49RtNUgtxOxJPfGfSEVJuUZ6BFAh2aisRZqp+nZffQZrBjsHqzboLRbI8wqZFpJX9kEzIFRGRYTfIG6g7nBV+HD5/8bKdDf9Ws1pX70kNiTpkvF8N89b/RwAXwKwy35/DYD9zPxOIvp9AE8H8HuTdvCFf/t36aFyewFU8uBT8yKUKzazrBwvMu578+v855ue+Xxc/9KLcbffto5Wu9n1L5GLdo+X7m/OxA7gi69sWuc3/ObFOO0NVkqyD9HyLUdwyMImcfZjXiWzM/Y0RntylAOF0Ypon1Uf3vKU2hyyTbFE9XMA+MmO01ntgZjEI6RE10Km4bJJev64aJe27Rtoe9/dMoWxncLYbeD9dbEl79Llg5rm/jyzmgR3Ku7zvP2oFoA99ns5IJgeUOUS2lgsIzL8ZOIhSIQL2b9AbWkbgFwkjLahkBnbiBgGKYbODLK8giLGQr9ApsSCd2AmO0GInfDI3Q9ia1CIPwCwx3fHdsTedqI07of1YS4kT0R3BvCTAF4B4LlERAAeA8CVAnw7gEsxheQdOKtvSqtDi9G4ACZnX7MjxL4/eB3a4Eg9xZdevj6pRRXixIEl+dHuCW/nEYzzTn4OoBXM0btApQHdtgoAuOrG1+HsR70SpqehAJQrOYpljYPHaO/M4iBD01mYvdtkGGp6wcw7LOSf1vzwaLv0Y+SatmiWSRLJ2A4k7GDGWNkTI2MmWe3pui3rhdfAW/IKnuRDa95b8bklsx2OH5wu0W0ml2ekGjBMLnHrrFli1a2zk5SNfFHiHIV1lAI19xIJoYdOVEUMpQy0YmRKHKy9rGxY71WYERs9AIDcDCllwDrolZOb2Rbu2njGN4B5WfJvAHAJgBX7/WgAP2Bma4/hmwBOatuQiC4CcBEA6KOO8hXcxkYtAI0XxEFe3JYroscsnxNMJh2MyQlrxwA3vKjT5NswPCq3EUmE4S6xujggJFBN8p97w8W428v3CyFlDBiJfjBBGvi0FyDysUyRZlLrvg1tRD+e/JsNW1e0TPpsjyF3xxe1RU82Th5AZS32gOTd9dzJlvxwLzA8voAaVFBJZ6UAT+BkSV0RQ+swAqYmcCL2hC0ljGQ9+S7bKrD/DMh2hgnGDvnBdUcAACAASURBVAHJdxiyLCR5b8kDMiqIWsv189vycJC1EDY6X050TTYDIvopAP/BzJ8OF7es2vo6MvOVzHwmM5+pl5Y22xyPcRb8Xd56+dyOMQ5eCrqD4LxTfn3u+wwltg5bBxessFPwX6644nA3YcdhHpb8QwH8DBH9BMQztgti2e8hosxa83cG8K9T90QA5yx6vP3uMWlYHmqzabflss8CTX+eqAaBV12LE+cOi7KCOjhCtWcRvCTZp2c/+lV+6sb7PXs/RrudNWm3sSn4puXeDY+tbFwxgMoWcyopmDqNoucidEACaJdqnPYdSoFjpJYGZtTM423iRoR6/8RjjTsO139czHQcQy26MxvyVQ7dSMZoa8VnqN8xiEF00688f8YTOnz4L6+/AiCgPG0Ve1ZWsdgrvHVtWGQTZpFQQmhib6FHFnlgoavEUne/pzDWWnfHVGAYiPUuf4HKtFjz/gu8TOTXGUtL8kBvdnazTZM8M78QwAsBgIgeBeD5zPwUIno3gCdAImwuBPCemfanGWOHzSGZW4eJ/052WB+8VPvedjkAahD8Xf7X5fH+AHzj6RuboMR76VmOb8wdT5O/6uY3+M8/fq8Xy6xEu3vIVuPyDiazkQsDe48Dh1fIel+95GLs++PX+GG3OLUIprIhcYbqe2+ooUun5C8/yI+1ZMPxb+FfhzY5aJyPaJ3gdF8zyD31xnL+7Dq1cIYhO7wnw1HyDFA7IL0j0r1rkHfiG7+4/Sfp6Z92K45ZOYClfISeqiLppGQlf21MoSNj1UbkaBK7Sq55SOZenkl+S+GIW/R9AzePQDjvgLvmZOUadp7aiWQPnPqOV0nZYkPSiZf2XZjy/G1lINVvAHgnEV0G4DMA3jp1i4DcJ1kW+/7wtfVb4h5wxaDM1C+de7BndDDte9vl/sW66Rdmn27wxmfXoUt3fZXINF954R1Tkz/nrMtAS30UKzmKFd0grtvOqKAPKvG7+N/qh/vrz5Kw0/u892XIBjmUrkmeDUmqeGUtV9c5MNVBCkyuAGjT+0nwDrdWMm8bNYaZhQHaKib49dZD+KnOHpF+y45ce1g6RGelO7KPSN+9/CZoXGg8Kev7Ui3XY5via897Ln7y48/GYlZgoIua4EEojY6IPiR49ze10OvfatHbcNBBUDuZR2SPptUu1a3FiTvJWmd7P+t1yN/DNrSGY84A4ran+DDhzDPP5GuvvXbiOqdfvh/lnrJOQQ+bnxs75E+0nfQFtBez8bJugOQ71DjnrMswOqqHA8dnqPqE3u1yPW87hTDcw6iWXfIB4lEZgJUTbkOeVciUweooh1IGma0BxEyoDKGsNIwdCvtb6eWKmpDDl8bdU2ohzWi9aHndvokvVrDL6UPvcRpj/HNbO/3+nWHjSN0ROQfFrQwBlR0hhZa8+93KVJzLDEaUGS+TqazCjU9+0fjzPYy45/+5FEctreKUle9jJRsiU5Ulc4WClZA8yC8L4Uh8mgwDwG/viDzsNNLRgpNoKqNQGYXC/mUGyqr2bKcU66cUtM8xB/fHhFJPYqiwteDZELhQ0Sj2G79yyaeZ+cy2a7czUyIKBQyqujB/ZYeoAETDCq8q1xeRIP8xxwSfEE6HjYEzQjVQGO0im0xCzSgON1J7RjxSu897X+Y/715cRUaSRVixfYEqjZ6uIsvJ+Jclace4SAXEL5hJLCn/N91PQvjcSu7tbYl2GDUobJv7wA2Sb7Py2BKBC5V0lj0zZERLNiLDW/T1JCFRiGgQXhi/M9sXK9kQC9bxZaBQEdelgVnZc2uvAquSuryR1d6CVNsPO4J6H/VnBQaUkXtDVftzEnx2z59bxnb0EG6YWvukGF9/ygsbbaVfGS+17UyS73DE4EHvf5FNKukf7qZ0AHD3v345AODLj3/pYW5JjXu959LD3YQdjR1H8tWCATJJdCCbccYgO6mE9HTeg506OIyVaVIdMuhxv/7zzV6yw2z40MdfjIc9XkIfi13A6okGZsEAfSP3yuroN/0/L8CjPvx8S+5SqmDXwhp6ukRPV+gpscTc8HhUaZ9dOM7qSp1mKVz8svtbGaqjItqkHtTWvnyul6UWf7pdvbD9N2+tJxLNOEEnlaTYKJmEIrDqwXZ6CuviEP+DVDsUJ23TR1F/B0hLVud2wwOvehGAHo5buR0nLN2KpWwITWyTjwwq1lBkpH4NmYb1PRkmknbapJ74Nxr7/BExtKvz7tsXXnf2+xF1AX4OAHnOAGYGcS1HMhMItU8KAG580vrltB1H8mlNmtP+/JWS2ZY4z9zFOP0vLgMgLwaIQXZSgHD8yulL0GHDYCVZwGoI8ApACyV6i1KQpirlBTr7o8+FE4J3D1aRKYOMjM8odGFppVGAAVTGrS9XSOyhzhq+6OHLaZhQWU2VSI7fRuThPtLfQ5JP2zRJe3WgQJLxOTLBefhoIiCSk1wEBgNgI+sYGybJzFafl+Wk5fLKpCtOt+VY420YOuSThrYLHvKBFwIg3Pu4f8NRvQNY0AVyqmCCWFsNOQ9FtWEwCV6igXPMGpgJmWGTiF22Z6cWWX6JHbGcPJc6kA0N18+M8m0KeMmux0aBCPjaBRvzl+w4kk/hyPyu776s1Zn61Z/7zej76X9xmUQ12Yt9wxNeIh3FZnOHOwAQHR6wxmIlV7WXVdDawOSEpb7oqUv5CIOsREamoZUC9mVUky30tt8aYW8gGKpJXrF0GFlg9YUvcVPvjy39Np9AinRUIBpt7QCOJqtKCD7cLoy3jnwI9nxIkR+dsv1MzEKCDKvZB3p9EsYnU9TVy/KsQp4JWT7sQ5fg//uxrU8cnIYz9n4Hx/VvQ98WM0qvuSKDsGKdI9GqxaI3TJbUlRgFhKgcgbLx0K4jSB20LprHgYilvA2xf8ZCsgeSDt4/FzKoTYldE3vHrW+zIYAYX3vii6deq3HY8STv8JXzfxNn/OVvT3WepqQPbGwI1KEdVY9QLEuJVrNg0FsosTwYYjEfIVMGA11ioEv0VBkRXPryGmKodeRzhxaaAUXx0ynZu+O1OdA8gQcdQGzNU9IpTA6p9BZjkDo/TVpyxyUK5SQWUjAKhiXu2hglFRLdaNSOAhRMIOPI/qIEMe98dbKnSJ9aG/SzMiq8dThx1MJBHNO7HYtqFIRL1lZ3HQJZNTvcRHapgmfCEX20H46te0DZdWNiT0sbwO0zIHsAtUUfbK/g7ivE2EiIPoyxN+7ezUFhOGJI/nDhLm95rc+q3QnJJB06zIqf/Piz8beP+J9z2ddd3nq5EJZiqH6FweLIV3MEgFGlsTrMvfZ8wxNeMpfjdjjCSP5wPxin/fkr7/CjgtEuwuqdGNXuEoM9azhm5QD29FexmBXo6Vqe8cNqZ+UkBksdDtcOk1jcLitOEXycc2hZpY4zZ91H+wr+6uQ4LiY6tfgdxoY+Br/NasEzkwz3jYKbH9Ttk8iAXAglGYmrBsSit9KNszCJXfgwN9tDLisTUNpA2ckvtDJYyArv/P7Z//tr+N8P/d2p7R6He73nUiz2RqgLA8+OXfkQfVUiVy4k0sosbKwWJajsvY6RjEZ8Ys1k3b7ej1Xaubbaldf/ObLyQ2u/tuCDPE04C5795xRkR1bMMlpjQ6gqteniZMARRvKHA3qpgNIGppKH6L7veyk+/9MvP8ytOnw4cBJgjhlh1+5VHL10AEcNVrGYjdBTJTRxw0GqxxCfpjrWuUpeCsOqfhmd44o4GG7Hw3fDSn6julOIUt8p+Js4bd1fR/yVqcnCYVwCk0Nbynx8PvVxlJWVQk2XONZxoQBj6s9spRupUw6RazjelgK5aFyVRk2MjAx6qsJyPsLITh/1tE89zV7zunN+y5l/NPGcUxx1/G2+TUoxBlmJlf4aALkXi0xY7g9xcFRP2veuH/19vODzT5BrCPb31UBDhSROsTQz5irDsIIGoy0Cx0k29T5lG0WxbBNmyzpZEGg+KzDKP1fjyF0hkGycLGgUylLB2CCFeSSndSRv8diP/bq/YVc94o0zb9cfFMizCpVRXst80PtfhF5Wom9DAte7z52Kfb//OvTPOIhjF9awZ2EVe3prWM7XvNMstHi882sWJO+kpioaAbQ5yqIX2Tp2Q8tdWVu9zpJsIXe3fjAaiCJ6EvIOa4yH3yc6iBGTg7FWvIK8/ZVRXrcNdVxH7jDSgbmytD7BiYFQ5wUSC95+V8oW71JSJ30hL7Crt4aVbOgd4sNNTCN1/K5bfdRUWGogIyPJbUxYqzIYJiyRQV/Ls/LTn3gW3vfwNwEA1kyOgSp8vZkcFaqwml2budtK/LEW7+CdseF9sr4Ld0zX7jT5yX9PnOlOp28QPWodPk6WCpzrRoHUxqNpUnQk34L1aJFf+tmX4X5/8xIMesVE0nriPzwD7/rR359XE7ctTtzzQ+wdrGJRj7CgC0/wHlOGy0CcmTiOIN00sCLrxCQ+DlXYwXjHrPJOWbfvcBgeHtPtPip4BW6uNwFjwzuTKCDDEiIJZcYSvQFBKSNWowulZEqKn4VRPXEddfddK7HQ+7pET0meQl+JkZLZEdUPioXoPJ/zmSfhjff/84nnCgC7esPou7tmmaqQkUFp78XIZjRnyvhr+JzPPAkLur6m2kdEqZjY2wg9+N2FXBpWrVE30q5m/Zq03W3PRDiiDB3+PrQyIPqoeXb9yigYQygryeyuCjE+5kXwQEfyHu7h2wg+91O/DUDCzpgJPV1iOR9hoMtWkjpi0Tc4celWrGRivacvRBzR0LS+3PoaDEX15A2pfu/2J38DbbYlJDK0vjKqOwXFymr3xtYjieOnDVHjeVAk5ASEVQzbyaERLdRWGyXlGxf54o5nNWAOYrErKxmkRN8WsueQErqU2JXfXAlerYxPRFPWb5JRhUU1Qq5K7M4O4rvFCqpk1DMNzaqPkhPhOhD3+2qVozQapVIYmSwZlSnkVF93kI9MtzJafQ+chR+NHGBQcIaqpbZNKBc2ZUGK9tVG8A3LfoZkLMOSr1FW7p+GqQimVF72nSc6kh+DeUYWODzxH54ROfzCh+KvznrzXI91KLHvT15ta5/s4OmGOqwbz7ruKXjTA95xuJvRYQo6kre46hFvxE9+/NnRslAXnBW5rXF9sMyxqyeOpdR6aMOjPiwFu/7+7J01I5JkFMtjpBcLHNe/DTkFdb4h1fkcXIaiDiw1RQa59WcoWMsyjY4AJHkFsbXshuIVE1Sgj7ZZaxIrXVnrj+AU0kwK2zessFSKcY5HJ2GMrz8etBFkrd8gjt9em9I0yzGH1rxLtoGVZAwArawVaYvFEzXbTBTH7zt5JpzaLp0wI1PGR0DlVN+DXJXeil5QIxgmDDnHsJqNOmJN22U2V8iownI2QkYVNAx2ZQqrVe63y5Xxv6XXyFnzVaCje0kHpmHNG6jW9y+1vOtyCU3pLz3OOBnRjb7aYFjkmdKIBV+UGmWpUZUKxlaVZENzD8XuSL4FG5VYJmUIPuWaX2445jIS/fFg1Ru73XbDvj96DVb2HkSeVdjVXwOwB3v2HsCugXRo/YAUCmgYjpNWAKDviaYmnZwq6BYJJ0WOClVAmBXcUBowFBA+JROWMNkoBuVJPxv7MiqE5Wldh5SWrE337zsXu56Tj0qjYTyBqIYMkCIK1QM80bvwPGohGifDAM2QzdSHkP6W2W0XsxFyq8fngZSSW7Id5D/E98pl6WAsnvOZJwHARH1+0nWTjqRCjgoDVbRun6NCwRoKDG3JHUD0ObyWjugdwVdMkRYfJkO1SSxpZx3JdhSvlzqU007AAKgrlNYhrq48sbHEDiaZEGQL0JF8gL99xP/ET3/iWVuiox/Tux2rlsyddeeswj35KpYzcVC94PNPwKvv+5cbPs7P/t9f8xELGRns6a/iHQ/+ww3t624vl0lQrn/pxTjtz14JnRsAebTOfU76Nnbla9iVrWEpG2JRuTKw5InPkakj/5DM3QsULRtX5gCQl9r/Lg44bS18hcoTfr2v2AJzWJhw3gYKCsb/DfeVjjCc1ViPKizB25GFlGdwIZ0agEEGoLTrKnuN4sgjmyKvOYr+SaOAHEG7z+5aRssSTTy8vqHjWJHxtWFyZXynq5NOYXd2EAVrDEyBg1UPxZTqoePaYFjIt4BGjsoTeL2dHLsCoTCZz2QeJwiG1yUl+NJMlhEd0Ycjrvo6JbNMBYZa2kmHUmyKMJPa2Bh4UzmCl3VmigTdADZN8kR0MoA/BnA85Am/kpnfSER7AbwLwD4ANwH4OWb+/maPt9UIrR4A+O+ffOam9fJLPnc+gMzXwS5YeysgJ7GcdmWrnjBcbPBmyD7ERjqOM1653z+up7/2CuAk+Xz0nW7FMYsHsLd/ELvzVSxnI/SpsEkrJTQYayb3Q+mcCv8Ct5F3KNs4tNX91lR5a0wDdv+Vj233Fv44Y8g5GmcYLVSsvJWYWo6NtnI4OqtD+yTeXaFiFrL3y5Uv2VAaLUTuyMVHDAWEEBB/+NehneSrBoFL+5vhnspLM+ydod66TkZDPSrRoxIVKwyoQF8VWJwyqXFKjrJMYcgKfVMi1+2130PkaYQWALjchwQVFErW63YQj4MKDQqgQeptCDtkDj5Xxso1lYonBwEk5lUnnvc5YR6WfAngecx8HRGtAPg0EV0N4BcAfJiZX01ELwDwAsiUgNsa73nY7+C/f/KZM3vKZ0UYbaJM/cDnVFnrqUIfJYYmi2OA14nFbIQlPYIi4wl4o/ur+oAr2XDM3ttx7OLtWM5HWMyGWNYjLGVDDJS87DKkl3UL0p7UXaTMNKQavUNI8FL8qYJhFRHQuLC4jcJZcyJZBEP7lkianOLORSxC5aWkjGBJR4jdhfEZJvRVCTezUTgbkdtP3KYWWSewkN13RcZq6jV5p20Pt9V+H8Y/h65Tduefwln5OVcYULvEEp6HM2hc+9y1c07bS7/wuKRtsdzVPG9Xdjq20J0F3wY3Kgj3ker0ci2STNckoxrJMdLRVYh6cm8rz3BM7pK3QICbJnAOnVIbNj1AYOZvM/N19vNtAL4EsfseB+DtdrW3A/hvmz3W4cL5n/zVw92EqXjVF38CL//Cz4z9/Zeu/YVD15gOHTpsG8xVkyeifQDuD+AaAHdi5m8D0hEQ0XFjtrkIwEUAcMopp8yzORvGX531ZjzpHy8CAK/nnffx52woa/Ul//yz6Nuu1A/9Fby2q6i2mgCx+A8a0e4v/cLjcOm93zP1GL/z5cfATb5xt6V/x0AVKFijYI01k+M/R8vrbne5bLxl8cAH3YC1KsOe3ir6qsSCHqGvSn8OklDE0FQhpxKLiekwiyUPoDWiRpOxenttPetERsgbW03HJOvfJBJSerx0WW1ZKqmRQuIncGWNJXOylm8yxM69rCV+uy3Mthm3Hzs03agpVyZyarchvNbuGQxlNbevcITl5Cvflhnu61+d9WY88R+eMXbd1IoH4uikcc+OIoNL7/U3AOQda1/Htd3uoyHnBfsOatsABj6T2u2CmnHzcZvr2kZhgptzsKZzDZByGbAyciBm8HTlakOYG8kT0TKAvwLw68x8axraNQ7MfCWAKwGZyHte7dksem6iYFK4tVj/1HQv/vzjo+gEB0VV5IBMXyRFBouoX/IXf/7xeMV9/3rscf7g+kcC0BhQicVshKHJI4IvbVQCADz1ml/Cnzz4LbOdwLFDLCzWeuuuXByrzjEHWKejASpSWEMOTQbH5bciR4uGOgHT6ro4QnKEDzQ7hHHhjGOP2TKIHdeJjIPvKDxBmMZ3TXVEkGFCFXjXwhBQAN4ZH+rJbYk40kZudQo7P48n47ADGgOv2U+Q1qSzsqRnOwRH9Gsmx6VfeBxur8TQeN393hVt+64f/X0h+uB7G8KoKSDu1Nqkshd//vF2CUWSpEJcwCycSUoF+9AUSjj1vXMyjSuf4cnedgBhPfqJUk2LlONKSSAog29gQEaBZ+TM9WIuJE9EOYTg38HMjpH+nYhOsFb8CQD+Yx7HOlR4+4Peigv/6enRskd/5Hn46GNeP3G753/uiegHOmVbaKAKLPdQM5WHhwAFHKz6U4nqLTc8HE5xy6lEwZkn+MKIc7eyFsWscc0Ou3at4ujFA9jbX8ViJtUAXWidC68LrTqnRWvU0S7T0BoLnyIJOVBw160m/vbtatJoQ+wIpMn7akE4uvBOYft5anSQWxaYldHyJILGIY1KCr87/XwyUW9cna1cKYGA6B0Kjp+tZ1738wCANz/gT/2yaSU9Wq+FnRBEBx2V93dMOReVPlukojh4ly8RhtVWPncCnoElt6LW6Wui13Y/zXrzYQftspW1AqSqnAJQgRXZYnfKH45aZ4HfPOYRXUMA3grgS8x8RfDTewFcCODV9u903WGbQdk49owMbi/WF8seRieElhKAyKEVyjYgiQcGgEU9nPmlrKAwNLmEm1kLPrQSFRgLusCVZ759yp5qLPeHtpZIhQVdYFEX6FPhzyscgYgjMbSczGwEPgMiWQFNwmu8zEB03m2QULwwamb97apD6Vw7VO0UDtdrsaRDqcgEZNKGNtkk3Hck+03o+JxzOO18U0NiWuccWfQWBWt/LddrTKRkHbdHeaKPciMmjkraw1ydcaBceQnUpO+rlraQvY4s/ZroFXFUr6auhUTIlEFpVDO/gQhk5HkzICgj2xWlhlEGtIlOeBLmYck/FMBTAfwzEX3WLnsRhNz/goieDuBmAOfP4VjbGhd/9oINEcZ6IRb8+vGwD12CUZl56/AfH/uqeTarQ4cIz/nMk6JoItehrcfY6LB5EG/REGEjOPPMM/naa6893M2IcP4nfxUGYs0DzWHncz7zJChi9JWEP/ZViUU1amTvhdZYTnXpXQdnfVRQNo6+rphXsMYL7/l3fl1H8iMrzwxNjrVAi0/jhH9QLOJ7w0V8f7iwLpJ/0j9ehJ6qsLd3AMt6iEU99HJNiLCtK2rNO/GA2qpsy1L05z5FJgmt9bT+fNuIwVmjacEpIFRlY0fitOF/23ppmGPlrfJgVBO0L0ycStcbVwLX6ewig7XE6c9YVM+1ddx5pm2bhLRQ1w+rRf/8uWfP/dZG8pmqIiknxMWfvaCWA925J+c97V45tD1rhqlxrmFpjLCdYbvDZKnS+ldKo/11DfMYSqMwsnKpW9426cyolO2LUmNYiL19/eNfOtO5pSCiTzPzmW2/dRmvU/Dus34P53/yV3F0//Zo+bOue4qtSFi/FMta4sZTIhw3rHaISBGybgENcIY1bt6iipWXaArWqED+RdBkUBjZXsoKEHJVYU9vFbt6a7i96OPW0WyO5D9/yJV45nU/j5wqLGrJZnUp7r4tQW6o+x7KBgpS7yW8HmknEV2LCTJPqjWHdUqi9fxQu40M6oSlWdsUIorPTwmRYLMs60SdkIRzVDZ+Xvnv0iKRetJzC30e4+Qpk5xjeOx0eXqtwmvtrmVK4G0IOykTdN5+P3aiDSkzYSLydMbAU6/5JRhQIxtbW+mjjtwizBLpPc5QyKlKrkedGd2WKe0ioAD4JLbKOb+hUAA+Y7mnyyRmvp5gJMpgbklqAwCtpNpmmSvcaqPj7v7XL8eXN0j049CR/Ax491m/h4s/ewGAulaHe/C8M9Jq1bVjcrxemlpf4QvsyqpqNlDMyLk9UqV27hifJASI1Z9bbb8w2r+IPVX6cNDwRX7Yhy7BQiajjqsftb9xnDc/4E/x4s8/HgOS8xxQEVm1jhzcRA6GCbcZKRpwp/yHwTVofwnHkfq4GaNSkmvTqqVtNP7Fx2yO4TbEvofWI/sIGhctE7YJgK+xE7c39mkAQWc15lpI2WHjOxag1qzbtPfGMtcRtTgP4zNqj0RyvpjK/nOZ3KFl3NyOUAZJTE+55pcBICJ7r/v782w9fQDw1n6I1CntjtYoiQE0CB9swyetZl8Y2FBHIXuX0JaP8aG4gIfS3Q92JRaC62uf2XodQllp/ODApGIbG0dH8jPCWSKOrDJVRZEmUp6ggAZHckwbuTsL10FHso2QdU4Vcq6Lcf3ZVx+MJ59+TdQmRQZ5YlVrMCoi5FyhIOsQyyTrcphk8Z7z9xcjjDJvm9jk0i88DrmSY/WojKzz2oaX9g+okJIGQZij6+xmIe1xpDyuIwhHDUBMlO5FHQd3nWdBSFazFVGz2yVknrvY67ZjJ5yROuaBFnmHlG/bOGki3TaUIgpHf/bW+MgcrqXDdKKNWMoQch8akQ1drZj0XYmvQU3Kwyrzz2NbeG8jFj3Yf7RPW7FyXCRbZAik4a6ylv1NooeaMfDueWZxzAZ5CSEqEHqqjEJeXadWBdfEZb+761Swwp7eKsyebRxCeUfAG+//54EVDz8cLRgAAyt6zZP+OGIH4AmyGefdfKlzci9ZU8MOIwb8UBuVjfCQeOyc63ovBWtkVKHPJZa0FEP7j9Vl7O6tYTErJsaqu3Oqz8VpxICzXCsm9KgCFDBA4c9JBx3CJIL0Fv86Jm5RAVkaNC33Wptvkr1uud4palkHY/dTH7/5W45myn+bFd8WIRS20R3byTvGyjuVvc/uGQk7XVei17XZj7Tss1IFfp+wRENu9++Wh5Ezqa7vRgeOxEqjvY4NwMech2HCiiQZLCcpNzw0WaOAmLEWdXpdXfy727+rfPmMTz/V7p+D8F72NZUMGfRVgT4Z7+9y9Y8A1DOL2ZGIH5EkHRxgiX5MmGoWvENO1ks7JheK7Hx2bgTkoupmlQ3Xg62J2emwY3HRtRcCsNm2LdmIHTp02Fq8/StnzXV/XXTNnLD/S+eirwr0qGzIMrEcU0sYfllQJMptBzgnWt0PP+60z+LN1z862bap8QOxNuusF58oZa2Grxw4Dj1VYtkWNHOlkPuqxNG9222bGANVYFmvYUAFcirFYgcaUoxMs9Ys6yqO19KffxvSEc84rT1Fm6WVLh8XEeO0W29pWkt40v7btPO4vn27czO9R20TlETbjYkoiZKGgpFKXTddwZeQzAAAIABJREFUBZJKEDViPzvLXJ6DLHoe0uO5Z8hZtnL8UCuXZaEF31aSIMz8diUXXFtKo1Fw7d9y64f17J384a7luAzwC//p6ShZSdSLbdvdV/7NBw605XX4ka4tZ7xmcqxWeasV7zBuekp3Xq798lucsNZXBRbVKOIJA8KIZTLzFS3zMgyowBNP/1Tr8dvQRddsMd745R+DInhJw9WiCZ2OKam3adENArT7dC/4u7/2XwHsamQ+DlTR6uA0AcmHL7d7+ffmB+HqiIfbODgHck4VelQip9JGP9TtdQQeyjnhMXrB8HNcYlN8DeIXJv09hRojxUT1bFrklrBWvCP6caURopr0JMQQdUJu/74CZdymcXLMOKQyDRAP/92xXLkEF0HifABhBmrl2yzkr2BgSHnHvmIDxcY66jNP7MZ3GHVd9vCapGGG/lwD8gPizlo0fIayyUR968zPgk7E6dXOoeqcom2hlCne/qC3jv3tVV/8CcD50IAomsmwQq6rqLMbmqyVHCVaLZi1KkDYxvDdcf96VKJPBXq2xpPjgR6VVkKTa3rQrL+MyiR0JD8nuBvpCNdbrmMIHZhMei4r020jnn9gjz4IAL6ImcsudQ+N7Mtup0xktVX2ZXUkcP/lb+DLqycEE0QYH52Tpskr63RyIZ7peTn07HnnqPDTp30e77nxR6L12rJW5Vq06/apRZ8Sce711NrZNVZnDwyzKNQTLmwvrSLV3H5aWeM6vHC6E2220g+pc8925gH5Ow3dkb4fYUAhJxcFxajsKDItQyD3lAFVk3daxndi+6hl2RiURqOERhbMQJVTHZklxyYMdOXnJHCYpVjfJIS6d5h/4DrBgrV9DtlHqLl3I/RDtGV9O6R1gFxEWk6iwzuCH1AdpJGTQc/6W9Y4w5pa29R5pujkmk3iD65/JCrIJApLamhvWk2aDm3k1jb0d8vbhv7hy/ndctlv5x6i0MEZWdbQ3sIcWcJfM2LrHpvdii+tnWRffrHmtHVUAfKw9pU7twoDGkUjh/D8XFsc1jjHAWuVhG3yck+LJe0JJ8A4ScYfO4l2mFXqkW3j9d0VH2czzkLebrqASDJKtmuWE9hAeQHUowe/HLEMET5L4UhOkpeySCpw8o2bjcnJNG0WfdqW9dbFcdUys5Bww9DclmihgSrwYlt5ciN4xb/8lN+361xCidVdi6HJscaZfxfCMiQuQibsKEKk+Q1O3nQkn0NIf6BK9FAhJ4OcGLl9BodMGMFZ9PKOPvguX596bp1cs8VwBC+JULV0oZOePrV8nRUGxBKAnySD2RM00B7Kp6x15oaDsiwe2hobMz7iDD0rAS1lQ9+eS+55FfZ/6VxoW8kyHCorYvsSmOhBTeHONSXYARUwoKgTCC2XkKCF6NtJuo1c28h5vWUlxsVgp5Tlrmab5j8ObqKKUNKoks47tfzHdfztbVfeig+XNc+F5QxcWCkb39EWnHm5oF6foBVDs/GdhORcMFTLeYdJT+l5jYOydWMqzsQ44XrSkjCzWkaXmR+5bgYvvtffiGwD+M4MRka8uSqQw6BPdiTOJdZISDYMiwYFE/1Ywyp+L53U5EbYdoIVVOhRZeWZCj1UGFCFPjF6RMit1NZjxogr3DbjSGoWbCuSv230BXziptPlYtiX/4Gn3HR4G3UHwAe+fk8Adz7czejQoUMLrrv5FDzglJs3vP22Ivk2XP31e0iKvpUqHnfaZ6dscejw9q+cBU0aA1WMteJTnTl08Im146z21FknFkMPFSoiP8x0w2xn5QxUgUUaRg7OBrw8wjYBJrPV92oN/eJ7fBBvueHhjVEAbDucTJMjTjh5zL7r8bGbzrD7j6MpljBCZWtkazAGVGJAFXKSCRM0ATqwDOsr4CI5JqPNqgwRbl85C5ybv80CscbH/5ZOL1fY6AyRRjJ7zJZM10BuSUshjE0CY+UTmVzymZMRQusxTaRq1gxiACVcbL1kLpfWJV3HthfG1jtypYYTTKrwCLRb9WkcugajUBWWSaJgFtXIS4QFayyajUwN00RYA8rhD65/pDiwlbxTWq+JHm/k2PVo1j7XVpZt87eFkmQtncZSpPP/5MTiVI6uJdAjwhIk58RFHv3zLXe20lktn41z9KfYViRPcDPTyMV0L8d2hiPAaQktbZqzI0AFhZxMq9NOw0aKWAfakh5Fv6ek6xC/0K6+iHQaitZwS3k0qmRyj0U1sudUz/rkomecROOO9Zh91wOQTrhHQLMzM+hR/UC7hzq3xC4k784/zWx0JBeez3SdfRKpt6+T6uTyV5E7B8ZaI+NTPo+g4ZLQ2u6bM0raag/JunH4I9BWskGkitBZ7uYduM0M/Hq+Y2Y75yxKgOtOIgyrjI8tOrsrLnabGfgs1vp6xdE2k+AlPtShqG1JRaHM5j73tUSeuDmDB1TUviZkUGq93fLs+JW7fczHpksWeX1MMajKSIZZUkP/3oWGXPiO54FUG08KJH9zsDdyDIBw6igFYEAEhQoFGRSW0DeKbcWiGowlKurePYj4aHP4HE782VcfDAXJ8hyXsg/ISxzKT0Bt6SrYEDFrzTs9OZy1KCpoRBXSlPi2TFXR4MOJrptOrJOz7+Lh+74aLVtRq1jj3L9gaVjjyIZhnnPql+JztKQdtsc98D0Yr7M7gs+JPLHr5OGtydnFEDd/89cD1FgW7YuBosUCr0A4wE5vjevia7BP81fE3pIKKxeGFlTq+HQInZ3yvSbWtuqTbTHtABplIty6EglSj5w0C9EUrOEq2o8rauYcrAdNHwdNz9c5CkMIx2VAh/HybYXH2tZ3VRz9vn3kifEO2D4VWNSjelJ4+5vUdzfAHDXqNlx410/iT77yEHlGiaGpgGKDHleNEdKAilpjb7HQnUGT+oniZz01bJrQ7r1x8fRRdJgbJU33D20rkg/h6rC4dOy0Dsh2gPOapxEtbqgGNEMm097dTVLgEyMs6QOIiH9yBcEkGoUo3oZrZ+6k1PyfPu3zeN+N9/U1apyVAshL3aMKj9x3w8Rr0jZVHVAPSevoFXHiOZI2AbkbjLfEZV1BAca4sylAkQUUDnXddwBS/yXYzn3W1pHYLA/cJPhJhA7UVn2bJVwFpCnf4ySk1Ap3UGSCjgG2amndaYRliusImTghzhG7I+y0vECYzOQMmZTYXTJTWHZ3EtxzrogBPcKCqnz5gTCePOyYelShp6stl2rDEGSXyOjKULjRrHOiOsUhtNrT97Cd3AWpcbJeGdEfiwiNCWQTbDnJE9GPA3gjZPT9FmZ+9VYfcyvxnht/BCMWXftIwye/cRqA9U/63aFDh+2LLWUqItIAfhfAOQC+CeBTRPReZv7itG0lSUESVgYoMC2c7FDgfTfeF17fRj2cTEOkolT5ZDgn24ZFm1APxxg+4cZZ937+UdTFxmKNL8wojH/321Bt0YfV9z5x0+kAEMk2u9QaXGJNeM0ffpdY2nE459QvecdreL4hdKI/Gm95xJb4NAveH6NlmTuic5COQudUar3PgDTMMbXgpY3NfIYwljxcXqEZY56WCHb7rB1rzSqQfn+BfNFWGTLcp29rMjJJJZciHTH4UUEzEzgsSiaTaDTLGaTbjis+p0iytheto7XOnM68Q/lQ4MmnX+NHsj4hKvAMSYKidSzD+tBIQTMniXaEEQOax4cFx36e+p1P5cEQqSwUjlInYavN0QcB+Coz3wgARPROAI8DMJXkm3U5Dq9cc/XX74GwLK+U3q38MM6nKkdk3syoDMndoZZu7L4hD4gQtCwTGWdy9uc4wg/liGlx3nI+60soAlz8tXwuoLCiRvGxmMfq52ls+ySCbx4X/hhhEpIJpBrf6VCT7Ke9IDGhNwl+/HY18Y0sIYeSTN3+Znw8gCCaBTLHqX+h42SnaFkLufu2Jufs9PexBI92XT5MfNLEMCQ+jGnE7v1Q1jDqqRKLusCCGtkJaeoJdypWGHLun6dfOuMTjXZsFSJ51Upebdzja9MzsIYs6sDiZ68+/7YEtjBPJoyaiqu+msjR6xzs4f4nYatJ/iQAtwTfvwngweEKRHQRgIsA4PiTdMPyclEDvW3gd+1R6ROKAPj05Jrkq4jkgfYkn2lZlYDtBLgmyWkTO4SI64i3XzixruVB+dTN+/DAU27CWXe5Edd849QJrYpx9dfvYV8Kdwxl71UlpOEzGOtzDSNX6vZOL4ca6ptt1r6k19uVmL0PAqivQe6tr4C8xzi5G1m29jn0LSUDWKIfN0NVva1tKDfrC/ljR5NkuGOMt/ClvbYksLXqXORN6qJ0xAnU0TJShrp+z1z9mGkTi8usR5XfTin2x43OJ6j741P9iZFRheVshD4V2J2tQvkEOWnjAdP3PoP1ZtHOAz9x6hds3oggfE+AliKC1lnfltRmIDNJiZETH6cxgksNWar9Q+6eS6SdcE0VsMi0qKetJvm2JyU6XWa+EsCVAHCP+/a5rSJeiI/ddMZU5988UROZsk7VOuxwkUboUWVnTTKezFPoMQ7CNrIPXy5H9G5oNjGKp9W5Ey+b1DGE66SFuWZFBcKKDcNsRKxYOGu7sGc/CDNzXRuS/abx9O6zJmAUhZ45Z5TbX0v7A4IPh8lSg78m0NARXm8rRG9sUEBK9CAb8kr2Jbd19pWz6icED4STmITOblf3XFEl0wdSTSJyLrG0E55fm0QQyj/hCEDmCs4ip/G4UU6c5wFkkXRU7zuN+89UFVnvfVseQ5HE6IflBA4nHntqLTQ4STMtaOcwzpKujYL4eQ7rCgEIMtpNXaMqXL8lwCLqNGYISNnqq/lNACcH3+8M4F+nbdRGWG1lYA8lnBwTYokKn9wzCW2kBdSWaBzRHocOunVzS/TT6qLMqjuPS7RxGcbX3XyK3fckCzUu2/DIfTf47erjMCZJRKfd+du48ZsnwNjPDjd984TWEgVpawZEPt4+hIvSGYfKyljOCi4gYZ06jVQISFOhLvQFSJEwH3UCKUPhz9u9zKRscbh6HtwwAiZELYPY6x5Y/v7Z91a+kwLdhDHsJ4xJ04aqYFttJwZxI4Do+EQAxErMkkuaPi/hKCBHPY9rm1kXTuLhinuFz42z3N2cxdsJD9/31cjn5GrqRJVgw/c/HF1OSFB0v4XJjgAiSdCv2xKLvx5sNWt+CsBdiehUIuoBuADAeze7U9e7dpgfPnfzyRN/v+Ybp0b/2iDROR06dNhO2FJLnplLIvofAD4ACaH8X8z8L5O2UcQN/cpHhATLD4Vs85Gb7iZZmzBYso7ESQkPbdmVk6SH9RTTqsBQzD4hQqGWPup14rkla2tJrp+zktzwOHQMA0L09zvlFruvyf3/p27eh0HLUDEsyiXx+bWEFTun5O8XbzkJVUNFBvZZq/5fv3WitbrZR+f48wd8YpVGnWAFoCHZVODIclUcjKhswomx242LfNBs589Fncchlm+9LK0CqViir+rJvdWYibwpHq4z/FA8nLDbWZLhhBOhg9CPAtAc1TntfNwosK9KGBtZJde4OTJx5+Cvj4s6CnRjh1xJpExaybG+tsqWLNG+6iXgJhdpn8D+cKCNZ5x1X8FE+S+zyJ2NXBKuLfsKkuGeFqprC+iYFVsufjHz3wFoFowYg7QuuBuKylDQeIfGoYKCwUBJzZXQgarQTHZoRIkEL8l6QwKb6xByCqNJnNxQSzk5uJZzfEJUHYpZkZs8uD1hCahrZJggHV9R+xAxDyQDhxHqqA1YaSB0OPs1qZmJes03Tm2UVT3xJFH3rr/lRDlHknN0JRIUI5i0I3R4O4kl/gvYRCzXGYKRW2euC+HMwbWzOEg2kVolVGcpOw3fdgQ60ezrZTL/riPk1rlnKQ6xVHZfFeo6MCHRu+2jsMTAGKqgPOGEmbJtUDa71LCq5SB3HVucxCEBG1bRvKRhjfVwPoK2srxDIzMwOf0+PFZhMizqIZ5/jw+0tvlw45H7bsDfff3e0iHaYAOFuF7QpMCXcbJqWInWoS1az2GWcgfbqp78Pe7b5z9634neWTdNezpoeltqzX/spjNsin+JRVVF5F6TSIyIhDle5hBa/0B7av+4/afHAOJsudQKBWLn4gg6sjTSME/ZjrDGmS1LEJdkCEPt2hx7YfEkV1J4QFWrQ9rFsxfW8ecwrn72p27e52twO4dteD9iB1dN7inhy7WqLXz3PSylMO7aFi11atz5uvK14/TVtmkCx2XJ1jHy8RMQlfNt/FZny7p9hG2M95s4XxHH2IeI9edmDoRbnhK6q63urPnQoVyw9sX21mw8vDuH3Ga/AjL6bCsqtp3gJsYZ2FmfUkybGWxcDtA4p24bLz7oLt/YWfXkw0I/7rtDBSGgdHabrUBaVTElFIM0Mqa2ANt62NByVVQnBbXFj1dg2wnU5DWO3MPfJp4PMXpcYUShtR3snWUdmchghHue/C188ZaTYsdpUs8mTcgIOwHDhBFJ0ozmOBkshAb7UcYkPPCUm/C5m0+2z4BYy4uq7hBDeaxihiZxLoYlFNoIv74+bj9i3cv51Ra+SC/toyaFuhJlRS4uvcLIlg0A4Eeivr47UDsq02gJHyNfty+06uUcwsikICqoRd5M9ytyn6pDTW3Ia1s4pFyb5rvoljtyl99jggdiEgvb6eZxddMPDoNKk2FHdPFnL8D+H3lno03bBY877bP4s68+GGuUY5da9cvrUVZtwLQR+rhQ0aJ1KYKonTjxchy2FckbENY4l4fXDYGIoeH+MsDKk9BW4+H7voprvnGqzQ4EDDmNNrT4woiXWI8dN1GEhvGk11aRz8fUMwAK9Ai0dwhAWBumjirRVMs1YXJVj5NM2Ba4ttzz5G9Fy0WLLz2hudFBYa3ASLt1oV52UomayMKsXFsXxZ6m8wlMQugTGZBCXeLWdcaJNJNcv7i8cU36xmr/4frKdhRO0ilCOcfp94mcE2ZA5lRhjTNL7vCdhAuNdGQfTvEo10jVRA9LyFa+DENp2yaSSSNwqmBbl2Qln4P1razktmvGvStbLyfpHO22VZjwxbXcpGEwoiyy6l1bXXvd/MNS9TPzyVlDM5/JQg4Fnnz6Nb6SpUyUE18/3yHbWzcpgz/Mu5g0RSigMYs+v61IvsP2wle+eaL/XHdG851kOMUt3zoBAHDySeJ4/eItJ/nfhGDW4a3u0OEQ4Xe+/BisbK/oT49tRfLMZFPANbS1JhQbGCIMUMow05ZVPVRwevYIKqpp7qzY8Hu8XR3hEGbBaTIwIBSQnjmUpCpXk8JKPk7SCa33IjhOmkFba/w1nNUJqmPC3QikgonOIYwcatPQAeBYNfTtqJzDEVLioaJal473GUcftGFRVb6GDyBRNWvMwdnV+4Ntf26tbu/4bvFtjCtp3AZlI3Rak6iASD4JI3OKoJXOiq9j6wlA6SWcCsZPEOOseqC26iOLnpQ4eIPIG4XEyg5lGDR1+hCpNS/ttxZ9IN1IS2I5Dqit+RDig9D+/vpfWUORQQFtR602Nd/JOKyiiTdcxJFig6HJfWmDoU2MeuZ1P483P+BP///23j3atqSsD/19VXOtvc/pbpqn8ujm8hA0xOsgiEiQkIA0BIy3YyKPkdyr8RoZCkSDEhA7uZfEcIeidA8yjDgw6asYDY9ormT44OEjRgURuICNvBpBwPaKqN109zl7rTmrvvtH1Vf1Vc2aa6199j7n7HN6/sbYY68151xzVtWs+tVX36sm63Yx8aMfeSb2Sfz9w6rlgMsVSE4L3VClVdcUx6rMmPW1nm1Dwi9xskg+qmuAGDkIFwxosQ7i2qWJ5C2f/Go855HvO6/lkkHXq++avCd/V7ia5UAYuY9hpbcEJ9K3iTjKQCJ5Zu3KadKSvp3DWnv91EFAoLErZiv1QK5HNkD2bJrXAIgRwGWn3aSTL7/nAkugmeexHUACekbuk7F+NbHra+oAqqlNSbapxwACiFMglKhwjNLZg3xyv8wBVaFfB0LNBKvJfi12p4rIR66W9TWUPW+SRw/n5X/L8yb0g1y2nCIht/lU8BbQcC6IKhwb1YdZpZZdLYMgl3+T9kON7pYSMXvGHc+uUOcDr/rw38OeiZNlbCvJ16+REp1xFvZq9Sag0yZwuj4lQ2Svgu92F3RPFMmLTt5GV7HQ8SiFi0vE4Jot7uR93O3Pr+oAQJx0+tg5s9TlYVLiKWBax1b7u4oHROgQNpG6+Nv2UMRZGN2y37uLqxkfpX9EAhTp30aCMUDU6Ss3T0X6yYDbKHdpUBYbQ+mp04rME7fJ8LtQd/HQWW7omFJ+sUPY9JxQpzyRlb/rORgatUQ/hUJC54njGBP7pGSfyh6IvuUmm/N+I9tGKrLPSbCCflv09baW+pTffCbr8TWFgX+C6PX9Etk3iD631zgSU87XCdIKKL1/D5UegFTchmpfAw4b1pg4Lthg5QJN/ePf+w787Nf+5PgZFwmvvOV6LGI9ZJOdsAeBLRLY1d5OmxwMRu6Ssf0WKcNtXv3vmp3zxJH8Gb8XgiFkazoaYFiCd2RJGELFD45p38ddIFKWuMX1bIslmZS33ldTQwxMGhZ5s5FNG4GE+2UPBg9SBJ8J3zKjT9I+0McJQG+wUfvz60AuDxSbbdQDu9ipaoL8tNEouFH6Yl/XdF2LFAnJs+UqAv6Sg2GzmbUzrkqCRMzwRFgk742sdpFAKJcMpFvUMnV9Gq6V5fnSCF+3nW6nHOznEfb3NCkGREv2wfgZK608cPT+qdpAV3vFaFJvHWv5zE8lJ2t5fxTbCVbJtsL5oEZI5+JE3COotHqy6Mk2vXHE8+a0WaG3Fl8c9vGF9cnc58CxSbtZAdFtlvWGL9IXzGT7FpgQFOpYA0mHvEsK9gsXVTRjxowZMy44TpYkzxT2tCSDniw89YBBNFaGmEHHBgdRir8QvvJPe9jH8Lt//IjgsmeGGCi0CLptzhssOwpZIoMuvUymVm8iUSSo4nG0XFDblAizNqXf9wg6fUCFPVeG3ODKVrpo6oyY2yLnmoniJvTquZwhYZW4ONYZOEWlIpGsAvGqCeqkhkG70pmLZO+Um6lFdpEUd0ggBzqJJN7vGABok0SejeESDKVXFrUaS8qt1VXhN8r4nIrgo7skKv92o1YN4jrpJ6X5bWl5JVVxKG97U+7ayKqP1ffXqwSdm76FsFduqdo0cTxlKTVsGG+Ii0079kzYGlDwDb/13filp/y7jXW9ELjhQ/8Ai8hNYkuQlbrswQuU/v6b9gLQSGq1Wl+v7XfxueH7Zmn+ZJE8hEANlgj7PBr2WFBoKKNSGhjyuMqePe9GV0HKnBc9fIzvcEALLMjhdncaAJJeXXKxlHm/y4hGAADbMDFEi7wsyVopheuQ6aSXU8bamvTDM8Z1mfJ42WYQbZ0beeEwov+6BBRRIlfHITXww1S2yXy/TKjy2Ac/5DZ89jMPDUt8cErnLHDR3UViF0JKAkrpeIGxqqUVhbzNgK7tDxLsVW8SIfdp7QmrU+5OZXMsiF98zGN2QgmemiL68Pt2EJNEwI5SEjfq3SL21s5V47YZ65oLLx/O/V57i+yxLQi/Z4sFhXTKSf3BJvXP21enRvW72JByipFc8uBvSusMjNtzFz09gEKllTY736J2PFEkz6CQbpQsVljAmazbtSYHhmxLnnW+IMZP0R9eYVZ42sM+hps//uQioEeQB03OS1LvctVzeHE9bPS2yYO9CIQg3zwu1xebiSfdcw6GaaH07Ck9ffajVKX13xotDxzJLXMVAfvUYRENe1c86I/x57c9ZNJTBQiJxgIhU3Hd4x76GQAhCAscjM4AUkpgTZkGcXNvlmjhWC4ee3+0kpC1zoXzpUtofS7r4McEn3YQQiC61m5ShYQPJKnep0CjkuhDXfNqcdumJS0ICWnUNqM6+ViLiPR2hzqXvH6OnNPChKRP7sglG5MBwxlJEdHBkMcBd8mN8rApdi8EpB17xBV29PyTnPgtMq/7SZ3SQkPzgSn4gLM9Y4te/kSRvGODO91+ufyRjkwMYJ2u26ce3/iID533Mkn63BzNGl8Um+TSGQYD5WW31KcaeELw9UD3aSlri9/r6EAh8ho2DRxfhJZvw9TkEHadiUvBOGkt4Te6VWrc2xD24krmigf9cTr+gAf/ydRPAACL6NbXs0cr4OlrHvpp/P5nHgbvRXId752bJ0EerWA0cbdy7gimopTTb3lzTu96hyAAhVpxrIYTqb2dNgAoDabiBROI348Sm6U6qTw2myaAOrdOi9i3bSTSsy2l+JG0ql1XRZJn9OywoPDeFxzyzIcYkq6Ihl2Qw30XZ3DK9jjrFviHv/tCAMDPP+nHJ+t1IeDYBDI3SBHlve9CHXxpeK2JXO/CBZSTYY1WX5OJUnL1b8KJInnPhLvcXsppsW/CkBAVhui913xhQ8uSV098Aa0EWiGc249SBdQSvE4gla5pejpkMhMPGqDhqjk1KUxATxQ2SQNe7VFbhtNoSX7KO0Z7lpym0KU0we+Cez/4s/iL267ZeI2PwXJ9nAzz5ult1UkNreuUDJqtQaU9qbTEbCBBX2VAXkuC189rlSeXRbaPC89y0Evy4GGeA6s4EWftPWOxeVOdKT18OJdJuUXwmtzHKp1AVknQqe+Z+n/lfglGxwZr6tCRw0AWvXHYw5DI3kQ709XdWVxpV/Ag3HZwbxycAL/54BcfJyXfwcf3IRugyEbnQG6jVpsBmvjbqyWNYv9cCrzYbdicBJi9aybx3z/96J02wfixjz7tApRmxowZNb7ht777YhfhksCRJHki+hEA34igR/kkgG9j5tvjuVcA+HYEDcd3M/PWxNCODe4c9rGMVvfBZh2hg8FVhgpj0/nAf//0o5MEBYjBL0baTjSXYxNVTJUaoFjqtpfP9cw98q3343l4Kghik5qmXvLphFELGkJYucnqgJ4tQA6SXREA6pS+sseqRMBaMP7crZuG1V1wvwd/buP5u1VcgrTjAS+Ud1E7tL8MYqPkV9/bl9ipAAAgAElEQVSyWWi1iaTYyOcoraqKoDVgpPbQ5QiZKMdJ6oCgqnOpfcfScssHPScq8zF3vJ2U0sP9jme8lGmJlacXm0Jy95XqYUqqNy6ovQx5LI3DnhkwGIsubhEYvGuG6DO/xoIG9EuLu90yGWGf+uvfBwD4jae95ljquAt+9CPPBHBV+u5B4T1GST6oc6mQ4LV6RtQ5wLi9Wu3Ygqxe9UbpUziquuYdAF4Rd4D6YQCvAPByInoMwlZ/fx3AgwG8k4gezcwb1xWeDe4elnA2uiOyQW/KFAJ6I+3XfvTp+J6veOcRq5ARtrULS8FF3Jxb8pc7EO5uEK5D1s2HHCNj/VhbD18a6/Lxlp6u8oaYUMHp3N5T9xSIC1YIWLLwpofMbS4Gdi3YYd8MISANHvvIOzEBiPrzAPGouYrOnwbwaQ/7GH75U19ZTGZr7pJOOwWVKVWGbvfRZh08nhAk1YQQqejTxbup3v1oanLRE8sBL4pzRhlKw6QSyp0FmozWQD+MkTXbkjYH40wZ5zdBVA41MSV1xIjAavLPKsnOBKLft2Ef2C4SfDI027Bj0mm7xlXdAdauwxf78x/xPoV905fquB1ULYCy36Fsr/DdVG04vZm6xjaD9JFGJDO/XX19N4Bvjp+vB/BGZl4B+BQR3QrgCQDetel+HsCBC7vFdMbBoMNgxSUpdNLTZr1zOO9hEDIu7mGfeshmGbXLXkvXrSP1dkVpYd/sO7vNBUuOadcqXd6Wt4QOlRZfZOcNHIVUDfumj2RvAA9YkzfpWCCn9rVE2KfoG08Xxk7y7Iffgv/2R19VtLmPfaNHmJzrCbTOIVTEN1SEqV1ZNRyCrtWR3jbRpxxD+f7lhiB1+gt5ivxe4j8AFJNLiq6uBrkWIlqb2zumYtUoddzmr70Lptw0wz2nJU9NWPq6QQlNxjMGY7D2Fkvj0BmHwfaB+KIhyFrGabPGA/e+iD0z4HR3Gmt/4VM/SgK1Xce9uF4DXjlZlNBSft1e+ngL55XkK/zvAN4UPz8EgfQFn4vHRiCiFwB4AQDsf+lVOHBdcLfyFoYYAxusKLhRnbULXGlXuLI7SOT7yluuxyu/8hePpQL7NGCNHGqtcdBISfDajz49huu3pPft0lI98JrBKA2SDudLIqvdzKb8scWYZqIHzcI49GyxZ3r0FHJwCAHuUzR8s8eCQ8K4vZgz/bRZwiiCMQ88v/vtanzjIz6E//ZHX1Ucc3HiXLNt7s6jr/FRum5NmPUuRho2ef/owBSVNKoR9LZRjQKPntW9uLxXK9hPYjEATHrWbMK2ICcpv/aqqX87vuduZRhLriVxDd6gMx5rsuiMx+AtXBfOL8hh31icNiucNmv0Nky6GC6sNP/P3v+Pcb/lmKgD4QfnC69UfGXOfgT35kj44XqoCUCNp5SYr5Ts5dhhsJXkieidAB7YOHUDM/9ivOYGAAOAn5WfNa5vTjfM/HoArweAe335l/LgQwczkTwHE1+8D0Qv1us9Mxxagt6ER11zG97/mYcqXSmlzS3CFnWBPNqeNZReok4OVWz3pga8lqzq7dfq3xXfKx/bllsWoIi+Ip/0u3i/FQXy3qOw9FyYoJsPKWB90tMv2aHnAT159DE47Yxf414P/sw5tvbRIe6zv/hHj03kLaSoI4p1rAKA1H8AFO9knHK5RfJKHeRLSW7SHqA8evT3cC1Be0TVCaiaE0RUMYW4ie3+8eKD3lIlTSUWC3YZIaxz980Qz7Ckg1ffa6LiqMd23mAd1Tea2DrjsPAOV9qDoGqU7TjPw6p+G5xqtyL9byR6nT02wweOMC7wGww8cdDNK6I3kfvqNhNGPRei30ryzPz0TeeJ6FsB/D0AX895w9jPAbhWXXYNgNvq386YcRS85ZNfjeX5s8HPmFHgZR98DoBLY6cqjaN61/xdAC8H8LeZ+Yw69VYAP0dENyIYXh8F4D3b7segpKcLVuOgrum8T58HNlj7DqdMj85s9g89V9ReE5tmTUkTG3K0t/yjx540ddizhI9PhTrX0vsm/1vBpmAdLen2HKIMe7bYYwtvDBYUVGUwwRC55AEL7mA9w5PDHRjw5deejDn7wC9iPbrkzdTDBsk67q/aWlFtCvG3E7ltal19q4130Zu23pONBt8FQl752ossrdrSuaCR1zjfkeDT+7+O9+fV12rVg4Z858pAa8Bgl8edob245afDyi5Ujphx1O5x45W3XI8zbglJdN5FSVzXq8aCXOxfNgcrUsifZGQMF0FhJhmx5RXq754jL1bSPO8ozR9VJ/9jCPvBvYOC+9e7mfk7mfnDRPRmAH+IoMZ50TbPGoGEcPfxPxFjiFFxa29xMHQ46BZY2TX2bPAC+abfeRH+69f9+yNWZRxtKCobjXd86q/huod/JF6fB1mhBojkIR1dLylbBC+ddioPCBAGubhT1m5ZpcU+E/hWcMh/4pkwkI0GTMLCuOL34rXiyODui7A8nsJzHvk+/MwnngjEdm9t1qAJHajdACubSMOYVuZUR6EiSPvmVvdtYTKhW5V0yiInvNLBTy03T3HtnVIVieHdEIdU2NyN9nJN7VSolUwO7OOQi6lF5OHZ21MNbPISqQk+qW0QJuWD6GbZmSU643DGL9NYOy7X0BZu/viTY1+6d64HYr6gysNKo8j/Xk0CBhR3aMsBZtYwOhAGH9xsTXKt9On7EN/h4A0GmGivyFajbWR/VO+aL9tw7lUAXnUu9x20TlryNlDWXQ5sgpHG+WN50W/71GNwlTnAFfFWki8eiFGGHCItZdeqmz7yDPRscdrI9aLv7dLvQwSjDNCSeGqC733stFpvqtwxpUOsfFdY4R3nzlFL6Nq/W6MlfXRkMZi4M5GhFHWo4WBCtOd5HFjnghUvEtH08R31aju6ekIV1NL9rp4LdYKoWseu7y3XT2UVlPvIOUk2pd1bQx0q8pat+yjft/bwKa6HT5HKweQuO1PFxH/aOAhE8jcponkqNmUXQaLQLWN6pSMEL//1teJi2ZHHnd0+PFNKBDbw8UvyP/bRp6WxLR5oveoGQfgLEa9AO9J8Yd1m25rasMXDoKfgANGRKwQ3A4Zhg8HbIOV7FEQf7r8ZJyqtgQar2V9eOhGjh8XadViZDgt7vOqanGnQpsAbWbpLOP3t7oriN1py1FumWUWSktis2LYuSo6b/JPrgAoJlR7705Z+yjhExzdgDFEVJs8bfPAwEompZ4sDHzJuXhH3eD0p+KeP/h/4Dx//WykQZeUXaTOZTd5I+tiUcVHyygCZpOR9tyS5qSC32ke+eEaRCjonnBJyaZYr5XUJSbwsZ2+wTQFxi5gCwif/KsYCsT/KpBhX0PWTZcUn+8HqeST0ez0OstRvkvtgSdpF+4gE78vJ1gBwPgh0a2exNhZ3R2+ang1Wvgvkdx4gwoJLgtN4JeM577JWT9LhfNtjDihjZxwbdGSxx0NaoUiQmSVOnob5f5DqRZWzDedXoTVjxowZMy4qTpwkL3rI1rKNvUj0cYZn01Q/HBaSK97B4MB3af9Y0Y+KDu2Alzjjl3Fzg7yU1nr3VrTptsyNue7qt2SS1J72VK2k+PS79EyT9J+DkopayaF0eluRUkSK37MDjAvulb236E2HMzHVxB10eud2vVD4p4/+H/jxjz0VBxykeHG1bWVPrAPNNhmrtVRW5vRmDNG4m/YsBYrnhWe1fdHrPqJdXk2UyvuoBy/dNMtAN0k/rTfe0MbhepWQ/pMk14oumBwTecR7Sprc2pXS1NJ8xIJ8lCpL1dUQSy0qz1QmHre19D+xeUufJXGr9BYHboGzfpFsduuorz5uuGj/SipbKvfvlWSEYdtMA6NWXNrGYtV+x6Ey+hljm5CobsWO1LGF5wG9sVj7Dn1UVa/RJYneM21tgxNF8gSGNR7eEzjqAVsv3nCYCJw34NhZvu4dL8fvXPfD5/RcyVPi2eAv3BX4S3dlykeziGkUtJeGLHuzAUhFgTbymdSDXH8XQ9cCWr1TbsasVTICnbtC2kZ+oyMLh7jcZUX4YtAWghd9p4nRhtJpOu+wog57dkDn3UZVwMXGyi8K9ZLYL1oG6pYu14DRmZB/SNp2AEYkb8HwalIN58akXud3KZ7V0s2r++sUsq2grFYucSnzIuZ6yTuD+UTWmuhrPb8Qd/Deib71iuTEECv6eSE2zxS83DxCDvykgqAwTjkQk2dbELsIGjIhauOhR1TfADCRxGxU2Zx1iyQIDj7sI/uWJ71u1EZHgXai2AQXyd4TJQuYCHzb0n4vEGwqyWBPhAU79JRz4CzYhbQpMb11Rw5D3Dt34JDBMwSOXkKphonKGZ6IAUVOzASKgQKGGc4bwGR95FN//fsOnajo5o8/GdcuLO70V+BOfwq3u9M445ew4LRBbyt6NBhCS73dUSBSUvIGiURduPYRAwi53QupqdJv1gET4sFQ6zyBkLKgo7ApixhzOhO6bEcenXFY+S4QYCRCAHjh+/9X/Pjj/tOR631c0IZVmbx6WZ1ECWmIq6ChkgAl8G4ZJ8CRV4QMWDUBaKKujYradqKP19fXz5D7J5JvBPvoHaRkEtDBWx4ETyZJ9gsgSJCciX6bq6XWQxu4pHuWRGpG9UlPBt6FIB+pq086+xwcZFTbAUAHj4FCQNDUQjfo6QlEIeXGYAwO3AId+dFq9rjw6j98FhBX5zrhWB186GBTugLDuW42TX7hHWzyxAHKiFj5rdhZerZpAyBYwHgTmtQiZwTwttmfNE4UyQN5wPlI8MyZ4IE84/vYuMwEEGNph0233Yjb/Wnc1t8Hd7l9HPgFJOL2wC+SEawebGLkk+UVMOH/PCGJaWhyl3u3pHcASpLyDQPsbp1eG7WdNxgorKCsCZGEXfTe6IxH5wPRd17OHY+K7LgxNdn6aFAeoofC2tuirWRylTB6vQNXRh5I9STQUs/VRvFdSd4kkhcPmywR1huj5JKVkn/vbHLfW5CDM4QFgspxgQF+gzukTJQ2SvUm5nVPaglF+OH6IMnqVYCJgkgrAlzq5imsMGUyNrH/wYd3lb1sAE+ymrDAIIKHP1R/3xUhbfge8n4ClMaidtnM7tHZr13DREO2TIYpQrlCK0dVeR+OO4RFV24TXCvjSZlD4d15dKE836j9P2vLO5QEtXbnVpUFDfjL4Uqc8cu4OXhM+O9JXeMKaQkIHTxLbNPLcq1vrfe5rDGVnlQkll55EtREIlJqa6IhYhgGPLHuG4XE711MkepNUJlF4ltGYu+iPrAj2+ywFxs5I2LZBpIlUXSXyV4RpXaZ0D0ThqifnsImKX9UFqUuK68fS67hmIFWwXXGwSuS15C+JGRTrwBCyg9GTw6LuLPSwgzJy0s/X6sURXrVbQe002rI91R+IXagSMZl0oQxnZxLhDqiHL/PHMY/M8ImMQOADmGSRmlzOm7odA96jMvxnk2sE2eiZUJPNsU4hJ3VZILKthONFrHrd9DyDjPkYUDo4OBoN0+6E03y5xs3f/zJWJwst+/LGtf/9ouT5AYc3/ZtN3zoH+D0hU9GOOMywuW8+c+JI3mRgLQUz9xemnkQiPPKBQCe/M6X4bef/uqtz7n5409On8/4vUJ/unJd8sP1bLBnhqSTFoikLTO99teV5faUDleuGfvdtqSosV+8XFurarTBVeoiz3IisQLJzgEgqcMAwDkbjbJBql8YB28JnfdY2mCQXZoBhrZLD4In/OoPYGEdlnbAvaq0H8exV+fLPvgcLCira1zQUocgMWXbSO1VSfmhEZDUBxotNYv2RhIJvIX6neh7bEJSC0UjuEjn6XnF/bVdIVzXkVcbb4QMrgty2OMsZerfhHuJhJoTt9VtMLX5dI3RBhm10wHaY7m4Jqpskk7eA+gAcsBg83vcNax/V+isn1JO6T892xR42Ee1no2qp545rfYHCl5Xi8Ljpr3Z9iZbnm53sQdoLhikPDvYJk4WyfN4cNQV0J42RW7ttNzbjJs+8gzsmR77FDxbeu7KgAWEIIS7h72Q2x6EtbHB+ySStzT0GbfEwdAlPa94qiyNiwM0GyrriMSWvlbu3xrMNcnnc1ScK1Qwot9ULqlifK11n83XkSI9ADhFQDu4hP7Nt39/MIwr/FXczQcArlqssN8FO8p1v/kSvOPv3LT1nlOQvPdpf01vi2CZte+SqgYYE03KBimBT8geS0A5WXYbVFW1AXY8Me8QuKLevaTcnYK+r/aUkv4aiD4bzxfGJwJqTTa6fFOk3op8FQHHqfD8qWhs6ZeibtFCHafP2egqKhvngu3IsQEz7zRZ7IpX3nI9rrbB7daBCrWt9CUx3ovaTwQdUV+uqStUZsB4Ex+tuk3tPBL0SltROFaqb/tkX+pS9P8mnCiSl1nSeRPcIxmTRDTywlH4+t/4XvzaU2/E8971nTjdrXClXeOUXePq7mwKVwbGUWeeCb2P5BDPrZ3F2tnkZihYe4svrvfQe1usNMR4ubRC9N2ItA9ch7XL+WoAVL/hBtHTSO9e/6/zgKRjFaHLoKoHiZaMRLLto+HLm6Cj90O+RtoZCC6s+jlWUj4wYegXuHu9TO+QCDi7XuDUsk/5h3ZdgdUQ8hpg07uryV0Tim5T/V8w1faF3h6bJfKpibaFTRPmurVyaJRdR4R3xmPf5v7aGZf+y6bPQvw1pshHk3TS60Pr87PtSBNhHdMh70HeSe+Cl5MIH85Ht2j1rqRuzhNMPO+ApqfYueClH3werrTAHe5UGo93ub1kaNUr9VadtKEcyI4jOQ4le2YJ9DuvV1PNncAqQU/KIu7Rl5SfvGfC3etleIG+JCmBSOw+qmkY2bi4MFlCed67vnPjs8KsnUPhe7ZY8aLIDyMva+0tDvpOvbgQViwTki6ndNjeW6yia6KuX+8tzvaLNIlFF+skhYlqo6t+55lwMJS/m1q9lIOk9E4Soq+9laQN9TPFfU0MTC2iecqv/Yv4ueycgzNp0DAD3leSoQlGpQOzwDJK9Ifdr/ObfudFuOaUwcoFSWftO6xchwO3CIEyvm2Qzu1QupYyU5qkrclqkhZZS6yBuPWOnARUe2hCmmrvluqBVH9roUVyJnrYiEuseKPIKkRI38h4mVAnttQs2i1UP9/XZK+IXK9mhrSKjP+VMJfeQyL6aHiNbcPegBzQOxvbK5Th95/1fzXbZhNe8oHnA0AM/APuGE5h8BZ3uSUGb3HGLdPYl6yTtTFd9ynNC8VntN/f+HdlAsNW+8rneqJ0lxzJe8LZ1SIt3equTYgeIiZa8w1AlAlsaUMHXsacNqe7Fe7VHWDPDNgzgUgcgn87gJSTJW1GktIcexhSLl6xJOvUwUS6MKgl5CQBRf02EAhZyN55g3XfjVYnRMAKcSXQBelUfiMD4ux6kbx+5J6hLcqOpAleJp3srZAHDqIzgwRGEXGaPJjjnqNR4POGirYYvbs0wWXi1JOJLpe0Q6obE+576szWZWcLZ+O7XLkOa9/FyMiwUhKpslYb6YlKyuzYFB5VV586SHpUjoOqJXDUbSCo610H9GlMeZFp76jDoI+posUt1o6IPk8AWr3QerctAm+fKz9rTyZpv3q8aCEgqRI9VeQubRfuL2MQAD7wDT94qHZp4a5hiYEt7uz38YWD06nP6JV1qm8UckQ9K8cEtQAkbQ6M1bXpNyrTaH2/1opdt6WLLqd1/64x566ZceLwxLe94mIXYcaMywYnSpJnJqxWMduDlm6UlEmG4ZnBnDNEWrVEktl3aQYsopFUIse8o5yhkAkrXhQ52sXA0vLEEKne+eyPLVKJGDGdUjHVMCbrVp0neG+yRA0kqbqHxbrvYIwvVgDMhFXfJclG1DxEnKTx+qlyfzeYsAmDM6WPdmxPYxkkq4J4L2sI3oRc/r0PUqEBw5pSutRtoCUxaWO9WtEqoqDiMCAfwrv/gq8oJJ9d4JlwZ7+flq5agh+cbUrotZqL1fvQ7/sv7roC97niTFlHn1UQ6b3u4JFTr6x0W+jfaJvFJmxrI6pWZZY8rOH4XssI5+xQwCMjby2da120Pl6nzSjGhlLHyDHpI9mwWq4uxdgavhMk/oeZsDrkqqaF2/tThVH4wC1w+2ofX1ztJ6k4qOwYnS2N1JI6urRzjVdtsjeBtHvd53T7Th0Lx5HatdV+wiObcLJI3hPcShdJ6zMCIQVS8mBLYPYgCjsbtXSiQQ0TVA46elCMRCvX5aCYhhfAwOMBnSNx20vsWvecBq7PgzN1Dj8e/EDQkFC03hNlIneDCSQvpNUg7Jog2AN+bcNvFNkh3hOG4RwAMvH3DDKAtx6DmkC0OscaX3TY1oCVc61XqQ2FgEzchMH4dO7RPz9ein/8H/4rfPkv/BsQgEXn8MB7LXBmWBT2EbEFDM6kSaceiPWgSwQfJ1qtWtPvyzGhH8pAuRp1v6j17oVr8MTvpuwsLehrqXE8EDsn4rcmk76oc8S+ICSv1Q61qkDrh1tEro9rFaEmp1oVk8g8dZxS+EmJCX12BLj1uf+y2R6b8Nhf+lcwxHjIVVemYwMb3Lnew12rPax7m57bm5jWNwpb0l9stMdN2baknp4pDC/1e/1uNgV6SjvU9x61W2z3bSx/LCRPRC8F8CMAHsDMX6CwTdRrATwbwBkA/4SZ37/1RgxgpaRNec8EgBhMAFmGNwbGeviuHDTWBJdGkeYl/LnW9dbuhm2f9bFfMjCta9XnNOHpqqXPyjVMOnj4fb64lssphi+zj2QEqm6qvpK0XZiQ0BukqqQyjQnfGwbZOKGQSQRf/kdaWQHhe3OgHgJEYVKRSWYKX/bmfwureuxd670kXYtU45KUY2J75YFQTC5pJRTLL+1KgOsC2b3rm34IX/EL/wZXnFoFXagz6AebJwTd3qOBme+v0WqfVh+Qdh71jUbbhQ8lidTXpAkcJfHoz2KDaJVZS5KahOoVbE1ItWSuBRQZB0XHVWMgFyDe04/z3O+Cp/zav8Bf3X0aNqrz//zMldjveljjcTAEz6/VuoNzBjX5praMh63xWQdevXNpK459r/i9qmI4XhWylvSlnesJL7VX/M8AuL5ZiSOTPBFdC+A6AJ9Rh5+FsK/rowB8LYDXxf+bwQSsDSpDf2gXGQeGActwnYG3DL9wGAYD23ms+y6WidHZLG3WagS5pl4+jesWznc2qypKCTYbj7RXQKGKqZefdX2BYjYu1Df6Uh86T+vFcmtghAqEf55QTwjl53gdKPjERxLM/UkRvHxPz5i4b+s8qQurpvCO0ioiPUMRipAU+yAJDs7g7tUSgzNJEteSoax4ikkxPlqEhtFgiV9rSb0fbJK+6ufkuleDsp4ENGrBoDVpSxsU39GETgfQvqA9SVNxPEwAelwUZURD0ImXpbHVIKZM6lU7q2vTealy3Q5MsQ8eLQDq7GoBNxgslwP8fjSiOovVukPf26TO1OM2NWFsk8FyrpNqg3qy4oGKuk3JLukRUUsBm722ijbT7VW31ZY2OQ7D600AXoby1VwP4A0c8G4A9yaiBx3Ds2bMuGB45BsP7543Y8ZJw5EkeSL6XwD8CTN/kMr1x0MAfFZ9/1w89qeNe7wAwAsAwN7nPjADgRyCsOdRzlLEYCJIHim2DL+0cJYxdB4ruVRmV08gV0pIbILKB4azXrrWcbMybhKDOg+7dLCdj/rosUGtlt69j6oCX87ytT66tdQvJZk8m4908VrSmJLm5ZifuL7+TTAIBOGpcavmc5SwJueTtCznlcTCBrntEVUmQNgbwIxvTiYkePJs4E2QpFznMQw2tnOW2EMbxRv6sfSTKpEfGvoZqTJLX15k3by84ykVUFoxMOXnorEsh1p5pQONho4/Jt2OExiZ3BuSXRYy63ei30OW6tMtRysWqYPui9MSZyFtFpK7PH6DBF8gdEiOu3c88kdvxCdf+r2TVz/sZ34IxnqAgFOnr8DZu5fwgwH3Bv6UARHQWRcM9C6kS/aDie+vep+qaGQ5GYJRl1+1EbnQD6ReVLVBbpdcPe4Af8rBR5Vp0Yaeyt9Iv62f3cBWkieidwJ4YOPUDQB+AMAzWj9rHGsWhZlfD+D1ALB/zbVs1kiNQb5ujHLZw4bgu0D2IKsaBrlhPFBslUkAW8AvGJKqQiaV4mXIMyzDLxnDgtEvfJwgYnWqSaVUXVAiVo7/U10oziuVvq+J1vIM+dh4adu4h+4gFca6+lg1XbdWp9wEUh+EPCOxyySN6MmTrhWVQ7o+P4hjThVaGXzqhS/Fw37qh+EcwRHKdo7vHMgCQtlmqlz1uTipferF34cv+5Eb4fYZvD9gGLKqz3uKRIDR8lyIQYQK0jMelRNgap5JckfqJ8UvNFHUbT01abeulc+67UtOK5ljSjjQzxgRWNm+VJdd3WckOHB1XIpgAB4IdMis4mfu3AOf7YCeYPvwPg8MYBcuCAyDge8N4IIXGnkEFafuQ1IOwzBbUvsC+feJ5OVePn6uzyHwUu8tuOPIaXKexu3Y+P0UtpI8Mz+9WQmi/xnAwwGIFH8NgPcT0RMQJPdr1eXXALht27PAgFnnDq1fftER4mfDgFkDWRQp70WpcVE0UsgUmgnFOGA0MOIg8JbAnUwmFqzIiSlfJ4ZhAOWm9wCMfuHyCOm0QmZKGk5Sbj2A9eBRx7fN5AV2HWjxcz04tz0rk4S6MSEMGCJkvTAVbVaIvLoNVFkA4BE33QhcYXIb15JNg2DGdaRiYpfndGcIX/6DNwH74bs76OANl6uLov2ilMdBckurRk3+umuOJuRpsmittEZEqc7rSZgt2v2nNSFUfXcbfdXvP/2uSfT6M+XPdRnqe9fn1LhiUwltU+W8YwG/F/cG8ARzYGDiUp+tgVtYuD5K7o6AIfwPQqEICiWZhoPSdyohq54cR1I8ivqTQxJASfEPeYLbI7g9TfJ5kij7Apr9vMY5q2uY+Q8AfIl8J6JPA3h89K55K4AXE9EbEQyudzDzSFVTgxiwK5SdoXho+VH8FO8AACAASURBVL8tpalz9bVT8I1BGEncUuxcqTPHjiLSqXwmiiqk/CfPzy8p1006rPwuGTptJH55ZmNgjghA2kuk2HogihRdW2D0AGw9o9XO25AMp8gfqPqYSIWUMTcWIV1DmewFUr+VCdHOSoKZJJkWhJQ5V/3Wl78kEDwQyMBTEDgI4I7hLYJxLBIAdwzY8ABykSAmnjlq46nybZDKp96Tfob0Md+FstZzSLFqK8YY7ZKWfPLZutwF+cj11fetQsnEM9J42PL7R954I7AE6MAkKXhxZyBnvwT8QOCzeeVPkdxl5Z+JPfxmtOITdXLNPbFPFxOfaoN0TDjB5T8Tn2sGwJwCTB/fSd2mDb47srrmHPHLCO6TtyK4UH7bTr/yQHcGhdSbJF5qNCpnct44KajP+r5Fg9eDLpL76Fo1K+tr5XpvEXRrFoW0qTuLfM7SvPqtBbijSPYoOkwhJVN5Dl7VQ+vqWK4lcFdOTqldpibIVjtuQWvpOJLuFfGPhNn6uJ7s5JijJAnVG28dVmr8xA+8ZFT2NLes44A1FEg+SZQMXgIxri6raFrxC4p8R2Vr9MlN2DiYVf80A8JKtV7tQPUV1fel0tuevw0b1S4bSCrfAKm/Tl4DdX4DzCq3MzlCd3e+nzUEGkwuW5S6C4KPzy/fXRiTZthSn3oyahGy8E4kezPE97YO/+1ZJA5I109MlNve27GRPDM/TH1mAC86rnvPmHG+8ZgbbgKW26+bcbLxFf/nTcDVF7sUJwsnKuKVoiSf1BvIUq6oNMLBCelAY5t0qq4hz22pr5aWq+eGWZ6h9ckikfsuS85JkvccIl8HwAycln1+AbglwS+C7l/+uEO5jFZL1qJNat2dUg8RIxuq1Epjm+RWS5et9tsoQahzYz29um/1PL30TWqwpMZCMMSkFQvBtIxw9cquITGzCbaYx/zLm0L72GDfSW0qUn1sXxpimeTduvzak6qsWum1JPfmalTKp4t8SKla15eqldwo7kRUBUpyrZ2PRiupiWfVdZk83xqz6l1v7C+orgVwy6tfginYM3kVg7iysatQX9NHO56hQkpO5d606ontYqbUNbU0P1F2OSYcEiR5TisEu462QFnNm/Zv8yp9sikAnECSX9zNSn3C4PgytN46XV8PnMMiLec4vcCxPlv3Pk7PTYOEGSTH1cTkLcEvs15NSJ08hxfqABo41dEvKZCHDWTvIvEnFQuQDGq+1v3ryU8ZdPTyXUi/nhy0SipUIt9Ln9dtNhqQupNtWmprvepE50y2C30vNamJwVyuHelH0XiPjYmbq56fJsLYns1xQ1FFE70vWE8wLTLfMNHo463+23o+t9p5AonY6/uLiqC2E1VlS6/RoFDvbZqkpupNaJzXk49W0U30Df0Mu9pcbxMnZCFoGgJx0hBsbHaN3Len+qGuZ0WqQsajd7qNh+rnMLLg5wIPmT6fDGpf2irobsPJInkHLO/yqgJceMGwzG6VdA2MpZVthqQk7fjQuPlEdZ8kZhG0gzPpTqpIXu5lDcGv8yRhBgYNnFYN5Dl2QgZbghnKero9A7cfiF/0+34RJ40OMXwaYIfyxXtVtvjZ9DKouSBQpvzM0eCSslSrhbIRY+ebaLvciOVv9P3T/9pOMH7UzqhJc6S/tHmCMH04z1ImTTwNaY0cQB3APYBCAKiua9ShuSpqSYQbK6fu1fjclCArSbXlCNAsC+VrinLXZN14Rk1Cozrr8dp4uVz1x6lVSY2P/OBLwupMGTWTBN/HsaknCenDaVlW1R1hjCZJXnTyrfaSW1bvcUrSFo7LuvkgAMoYlvHHJgp/agVe2/c24WSRvGd0Zzy0xAxIZSkaJGnkdSKNVd6M0m/ToZEUxZmYFfkVMGIt59G51OF8Jnm5J8AwicgoztIe5Dkdy2Xg5FFCzoM7YU4D6gBjI/Ezw9vgAVBIIa2lnPqcOo42/nK8oIHcuVBIz9K5RpNC4/+oPI0JQp4lapBiUhkVqlnUnVC0ByEEfPns0ZAkdy77yNSqMaluWsv9Vnn1XLBJaqx/16xM4171O6jacEToDXIvrlNlIe2yWJH6RkJvfJ6aAFtOFdSoq9zrfT/ZVtX8jRfdBPKAPRWbxFV193EssKoTpL2qcZlOclGHmuRzHUvOaq546/Pqt6l/RUHQOM5OBXEMekuApXG/9Js7zQkjecCeDS2YGo6Uh0NnYiduVHQELgfBSBpRDasniMItT0XYFlIvxfJy+XufZ/065JE8gwYheUo6waTXD0WG6R3ADHPQwe9buKUJE1tHGPYNfJdVWCNylbqWzRDUNYxiEsrSHBfn0s+E2A0Fb5JEwnJc/em2aXkkJfvAeDIVCSp5JMkfyuuOBPXudZnMUE4s2tuhHqSpOAwYvXxGNXAbfbHQcW+SuKe+y/WN98zVeQCFJ5C+56Rkrco++pnS709NEE0VVWMSCfcZV05PGlpSH3mW7dgXugNVR1GhWvGGIthexkAed+qpqqx6jIdjvqNiRZxItqqWFlBzJfMjasE0KQxkte98UCMrAheJPoxB2rk9jiN3zYwZM2ZcMHztt9w4OvbV33HuG8Ff7jhRkjyYYVeumOGCJB8kXxq8OpbP52vjbdLGqdX9Rz7V1UwuqhNR30SpPauIqtlTvHL0jFzp6XPdAPI+qAnYg4lAVq/n4vKMGdQ72PUAezehMyHwh62Fu3IRJPuFiVIvtQ3TVN53k+SePrtsL9BSbCG9dxSt/hIBHJ9NnJ/fCrgSw5J6ptxbGp8tJUned1Tcr1BDNLCLDlSkuuQZIauL4p0hSXjNe5KKyK7VIhNS/LbybdSn6nMNdUxtP5HP5EvJd6QOqaVO9byRPUOO++p8oW5Q5zZJ7S1pt7ESaKnMDE+/f4FdRY+1WuMSHSHS+3UI401J1KNnR+nc9kGiFpWIX9roqMHpulF5a5UPKU7S7TJaTWX+EUkekTO0VgHGhM+20dEbOFEkT55hzqzjF0qVKpKfaXdFTeZyvWpkiSitURpzMqkHdQLnxpVyRKJNz5Ci5MTa5QMaxiG5NxyDnEs6eH3vgugHDwwuTAxDUGHZO/bA+wv4Uwv4vQ5+aRLpCuGGMnKuH6olp1qiFoafVO/xxBXKGf78gsIk0xG85WQMFy8AbdTTA4hqktckEQsk+n/fcapPIn7lUjZu2+q9phP5uExmdd6Q4j0W7VKqDkgzofyjBoFuQEv3vpnka2EHeeIz5QQs9QQhG+gkDcBUObaVp0HA9e8KlWejLpN1Vu3bUqGm3+ubMUpvN/0cz9mLrHquj6rG5V1B7RjUIA2Cj8e1rc6uXCDc3gGOgU70kWrs6/EuJByvKwy6VTuM6iD8ELLhZR6KzzCK6Iko6A13wIkieXgGHfTJa4FahI48q5GQY/xjmz+n3zQmgDxeI/nFFUJq3MGnGRRA+Ryo+021cepoNfn7TOBirDWK6I1JEwFkUvAcnjs44O4zoJWFPbuE2evAewv4ZQdemChdm1FnSuQdv4+gJOtiwlP1T21oCORMyAS4MKAFwTPCBisI2fn0AK5JvZxUkKQligNMv1uxQ7hFdCldSvzAOAS/aSxtSLmF8dQxCikc1TXVRFCe5/J7o20nvUBqMp3aDaS+d/xSrKx0VDTJuRzxWp9v1jUeaxrHebp++R48+v3URNGqd3MS8ONj+Xg48KTnvga/++bvK+pCDBiX40/k3oaA7iB8MWuOq7lSAqfqOUmaXjtQ70AuahgGTI9vIEvYgy0F1G17Osq9hOS1UOg94CJ3aFvfJUny7IH1ujoWO0PVSEESNoA1obLWhMal6k8IWkvM+t4+Epp8di6QsMtePiMQgaLxQ08krXJnY2wk7DiRQDoNEUjXw1C+Dgjfh/gbKafzoN6GDri04K4DuiA9JKOMKgd5JSHUTa6vlT6j2wQI3ijGAJ0JhmMxWvngNWQ9gs94i1CTyygKUq8lKDP4NLDgOKqoTEglvWfg9g3c0oTcIxOGRd05tPsZTRD9FL0mqb9F6FIX9cwWYU9K6Oo3tVFvhIZBD4Tk8aVzHyUPtKTG42zIrhwV6rpuksR1XZpqpmpCnCTxFoFrNJ5dTgCZ3Fp1+f2bvxd/8x+9BjQg+Z6n++j7NVYLRTFEFkyqE5c5IY0Hn8ZuXTYi6XiKjFNBSiGxXEHy+L+MQ+ED5xFzl4/vvQEni+RnzJgxY0c86bmvSZ+1VD+jxMkieQbQ75AsOs2I0WHc2jhDm0K6JiLAB3UIIVijSmlELYuiFA8XZ8y0RKqttVHcjRIzUSU51zpDr54hM/BqHfTseia3tlyCaYl+GIK6poupLZnDd16BXAcyA9DZoAc0JqutBErPV+sQCUgrIE7PzuoTIEiGQJT6owRiHMDsw4oGnFUGyhgkKpnSyMsjV1NiBPfS3pdSExFMZ2CXHdyBhd8LUr1vxUsoaP9+tqx02JTe+6bApGR8re0LwMjtbqRPVuemQIxs6Ja2mLTjlHXUdihOGUxjW0Q7hlsE1ZpRUn2t1knlYBSSbXJXbJUb42tHNpa6HqotdrY/TCCtrEjVOeJJz30NrAvqRjaEehUmNq+mgRUoxy2Fth2Nbb0S13DquyHAufzb6DSSqhk1CulIvRrQqiC9+i8kfN6pvQQni+ShCl835JT+SS+dvM3EayPh26iDZ8p6/AqkdOBJTcEcCH5wWVUi5TDh3uiq5ktOxYqogfhSsqoF6zW4H4LKRtdPvIhUIVlNMmQW4043DCFayvuQWY8IZE2phpHfpyUfjweANSCxa8gxbe8QgoxtFQKClNcRhFSpfB4wJvVqOR1sAR4kxqbBJz1kmIQMTBcmc7/s4E9Fo3NnQoBcRfaiuvBd8I3OcQW8leCy0U2pmrShzpWDq0ny2yDqEZlMxJMilSOqOobcSHlHNJNJrlITSjyJ2bPREwrRa0kZrZUBt2Ucl/IV7anGTEtfPqr3IdVWrd+FkyJcVIc53CDxY3zfAjO0Hsrl82oHhEicQeCIKs66PPW4EcIF8tgXzx2iwAMmCJpSH/KceUFNVCMnjlo4FA7ZVJ4JnCySZ4zJXeD9mOiVsaJIO0AEDEHPTc4AtmGore+jpanRDB7JPjws6v/t+Pc15D7JiBImDF73gHNgXVfnQKayyGvit0qKr8vvXejp0oG8T947I+hOI1CdLuxoFIyeqSxi8BEDdfpdrmcusxB+ozmS50LZ3oXBe/CgQYxcLtRNIoINwXYd7DJ4GPHCwu13QUdts3FdBr3pCNyLJNvW5Y48c1hL8Vy6l4pBrvWua7IqdLHqmiRNctkeLgsYSegYfHYysJS9NrTwIv1LupIN/YAXBtQZcIewZV0jgE7KNDKI66rIdVPqX/V73Z4tEh+tnrZEaqabGCpXIYndVdFimxmHciLW5a4n8zTRcyGI6MDFYJ/jvLr3HPpk0iZMcZKPfp+EevWfG4LGzarJXpP7pFQ/ZeTIOFkkD0zO5qNr6mUOkMk+zBZhBjU25B9XhFm7QgJqeSSS7jCoF+uDRK1mbXIesJVPf9VpmX2SytPvmcF9P11PUdlocpfVQ6uNPDJJswe8CfU3atKrB5NIBGkFZPLI8wjk7RgwkdAcAPLjDqkl/GoCnVxJMI+NwMzRRU2tmGQlpTu246DOO1jBHCyBRQc6vQde2KiqoCThhlWFKVxLiyYQotPkL68gGokz0au8QyJdax7ZIE1N9jNtkBdi1ySv2oJiiC13AA3hniO1SDIGGxhrQt9zDO8M0BGQPG2UIDRBzMWxbWisyqZQCB4T1zUnGZfTfoT3pqThqiySIwpAmVhNVkp1nWNZUsQ3hxUUuehVc9BnRwm9Ck6qmKrAur/WQmOKVK1W+o16hHIpI2u9sriQ6hoi+mcAXozgXPRLzPyyePwVAL4dgSK+m5nfdk4P0ISupfnWUqqW5iWIQRrX2rRRdeGtoyXtpI+PUrdzZQN7D7g4HuqVRT2rMo+90rTXipQTcSAnvXxF7jVhFhPClo6SH9wmfQOELGfI7eQBks6rVU7pN+r5Skesj6VS6XeSysJJ2koqpH5AsRytJReRWpKk5IC+g2EO9ggb3q/YI4KnkVFSfilx5gmBctoKKU/Sk08Hp2xFbI+iHYq6IasGZdKTvid1F4lRbE7OJzVk06OLKHg+OZ8KYOBBTGBXqdaoJviGumlHHqltCqNgnw2rgKJNJh+QbxDUcyLIxdOy0nIeZu2yoKJWOeRcfq8oJ19S/YtEkBA72KBshPUKzZggSLa4pKVu9lGAQ9XPJ+paSOyta3eQ4oEjkjwRPRXA9QC+iplXRPQl8fhjADwfwF8H8GAA7ySiRzOzm77bjBkzZsw4bhxVkv8uAD/EzCsAYObPx+PXA3hjPP4pIroVwBMAvGunu9aSbqGS2XGeqGdcIKx1W+of7YPO3JbeW2XcUhaemmmJkv69qZ4R6U3OyfN0vRr+8KGtfENPOKHTk9UOEcLyZErkqstv1Eete6b2fw3xGpIypz8tpfux1FLUJ/6mj8m3rQF1HWBMMDp3BtSb7G0kki8wluZ0sBwqdYr3pWeS6M03to2SOlurHCm/lEPqukkl4HyxGk2eGzb3IY6qAxoAGLmnCeodWd1ofXbrVe+qARipeZS6Y8u9atVWS603/pH8J2gxMdk0+hCRSusBtIpxNrUu26t2RlV9GWPa6cIzxOlBO0IU4865cR9vSd1E+VrNKa0+PmUv3IQtbXhUkn80gL9FRK8CcADgpcz8+wAeAuDd6rrPxWON8tELALwAAPbNleBhBxfKbWgZdDS5tLBtKVST7BbdIoBM5KMTMvjzIA3qGaVWKlQhokBulHmK6DU0aWoC0f89kEZmq2MVHWlictvU2Wp9ZCqzGoj62fr5hZ4+D8Aw+HpgoLC0FlWRNUDXgTob3Ge1S6nUl3k6iZ8m+VpltIuQEetIqc407i+qLmBfqWpifWXy15NxIn8T6myCIEBSNxvrVrvSpolAf54gqPQ9/rSlC/aNc617KVUOVe+Vpdz6PqpPNydJoqQySwb73gGrHlitgHWPwn6WytBQk7ZUoUDRLqQdLOq6pXtvmx3VeRMbI/6m8J4jU5a9HgM1n+wyQWIHkieidwJ4YOPUDfH39wHwRABfA+DNRPQIHEJOYObXA3g9AFxt7887+clvwwbi1dL1SJpOJ2jcoEA5+FrYUUdW3GOTXh0oBslIetfXaFdSfU8txQhB6cmhrquWZjS2dahdOpyuc13+FibIPZxjMFzuaEJ4ScodkmRfTJraACb31qu2mnzrFVA9sW5qi9aqK5Fk65kT/adBmMWzdJ2sDa604hKrdc+ip9YJrqZ0wDWJj8raONdqg7p91fmiFafONxwlUpi/rH6iDp3Xffiv+42+XyRsWi5LD7lWn2/10YZjxVbUE0ByUx1L9AyMbXYtW1hVvklhMmIryTPz06fOEdF3AfiFuHH3e4jIA7g/guR+rbr0GgC3bXsWgKAqkftPFH5SDdIqY0UsyaCpjrVJtnHMTpF7fBHazXGnwk2QhfdjF02KKxEt1egBNAzhd2nQm9SBeBgmB9qh0GgnqiXzTRPWLu9tYnBp7yR9bTHQXCyPeBYRAc4gpMBUREgmB5bJM2qpWpN6a2m9re2KibNBpC0iqlFLb5U6oXxefOfWhrrFeo7UabLaaQkM9QQ0JbFu8vpoYROBTgkt+rxVKxVxr9WOEZHwWVKFNNrGXHEaWCxCG3VdMGrLe5axMlX2CSTJuwUtjat3XUxAtY89AN6wKqDkeqk86FByZgtHVdf8PwCeBuA3iejRCPvdfwHAWwH8HBHdiGB4fRSA92y/XSkpnROZA6U6JA7ukV5NX7uD6mVcVEU+LdVCC80yTpBFrR6S70XyMmVHiL736TfSoZzPnUB3yFZUGFCqjVBNkpLDRhez2PEa2Emvv2l5q9tjk11kSgBgD9JupPXEl5b8LhFh/rEieJk0WxNM3adaHkhwk25+uqw1xvesVG16ma9dTpXqgRZd2a/1O9w2OdWripEKZ7yiKs5NYWpcTK2YhQcWXanGaJFmizAr33Re9/jVP3tdesSzHvjCeJ0aX/UEUb+/hsp3KnBRyllMPOwzie+yCqjAXp7lwiTtMDmONY5K8jcDuJmIbgGwBvCtUar/MBG9GcAfIrhWvmj2rJkxY8ZJwLMe8J3Tq/LLEMSHkV7PM6429+Mn7j279MPGhNqmJaUIUqpiM9Zj15LMlPFlCjIzOw9erZK0XN5zQhJuqImKMrauEcgSVSTzYYhLVRekgygl6OV5sfRrSQ4tI2g8TrUU0/KoUWXfphcMt29IhzrMe1S+ht6xZT8ZPahaBW1qZw2V5ZOHnHZipCedenb93Na5KacA0gF7UadeG/20FJ8yE7r2cn2bznZq/BSpOMYSfVPdUNSz6mdkzlFqrSRm3ecaz9TXk7WJxN92108DAK6zzwNF7zVz1VUTqy+FXeNPWmVnDz5Y5bEa679JFbMLpsYkGcLbVz/3PmZ+fOt3JyrilQGkqE3TIEnBJrIUaONXrVvUcK7w6KGuCzq7Wo0jJDu4TLI6gIdM8ULlGGKHI2ubS+vpxuCSdJwDr1bFoG6ROFdLOOnYgrcd/CwA4DrznHDe2nA+/t6v1/FzSTBFZ9X1i8/gWu1Rl+8QaHVmrs/taguI5DlKFQEUbTy57N+Cow5cAMGWoNswpuNg67MtASiJV6njapIblRGNNq1Rj58aKnfKVuGhrlzrWTti8n3rsrSEE+fxtrM/07ynP3s22+YaalMis53M67YSFZbz4L4P6tMt7+WwGE980gZHNLxeaLBnUMNrqcDIvdGEHV/MInxvGataulUZLP2QfNTZR79onXws5q4R6U5Ilvb30iW/+lf/oVnUZyyeH4x/tjEYjCRNm5ZMk2GJI9G3JBwyACzevv654tx19nnhXBwgutO/w7+lWd5R2QG8vX/j5DXXmefEiUVJUmmCpvTb6+zz8rPdmybvBTKNgZGJovCTLib1sqMXK5opUqlJa5fVziYcRmKtJNwwUGPlnHgKTRDvOeh3typLiUd8XE8MRyKsKRtQq35TYf+ey1VLw8bEnvH2dYPc2YNdHE/9kMYNWZttGJH0OZLnyKmgLktKSuYzL8T/x0nuU0gr9y3v/0Spa+5l7sdP7J5ZLlMr1U3LOBKWZxa0XIyuK1Qb9eAQaElAGR7JmCzhSeeKZXvb3W84cn1lCRmKsNt70JK5SOX3JMhE0Az6wpZl/RR2IWft+z76+UR/ugCYVGucCw7TZpO3mFpRTTg+jK4rV/B6hSWryVadt40JWbm2ypQIXyZWa8o0IxqGgH7IBlWlLg3lvbB8Ku39juFNk+qae471YcaMGTPugThZ6hqOkvYwXoZNgZjC7DlK5O9GhsniUb5UL4TvShVjKPoch6Xc+ZCaW2qLpsQRClSoP+6p2KRmus4+r1zCtgzIUMvcBHt41zZt97AoJblzkYgnnnsYyfBQUvwRpPaNKhbByK6yWS3ZQktaHtkDRL23c9pM/QBtx+KxKlBL93K9tYUt7kJL7jV2ef7JInm0Cu2KJXhq9HhMq02esfxHY/VHw9oflncnU92xi658RhutSTP1iRYpK90vqcTyLaEg/3hCXwwLtOxJu5J3w3AdDh9B/dK698R9NwbhbLB9JNRqVWDk8DDyAqo9igSVw8GUAfNYx0rlX5+Fvuq644jIPwYcRs17onTyj3/84/m9733vxS7GjHsInrF4fmmHAabtNhF1sNgk6nG1wSNlcqBObT15Dtg2WdQr28nnNqX4hh1sR0lXT8BF2zb08FKuc1nNTq6Qd8Ex2CqaqAO3ikcq21K8blO9iejScKGcMeNCQgbNMxbPT8b7ENVoknsasE3KnYi81epD5nIS0ct/aHXh6OaHqM1m7CTxTUnsDdVLuCRKvRNq0bf3b9xKrkFSduk+yYiqiZ8J1BnQogN9yf231uPYcT4M6crrrTycJ73aW+5cMZP8jBmIenq941Ct3hG0gr8aaTEK75A6ViNfVBD9tvJtuWC3302sDsJEU9kw6rrWwXoxPoVi4NGv/H8/Pnq+qFR2kqTFzXFkM4nvwNqdgpF2xkRg1XnBDioz4NxWKdswe9fMuEdD++/PmHE5YpbkZ9zjcJ15zihlABlq65tbBkWU6phaZTNS1eySbbAhQZ6z50YjhiBLjbZUoyTdeQzGYvEqiZHnNv5uuRgnTyMK3mfL5dYi7WIklRiIFKCE44lHmcSFimfYEKh3IbzlZpKfcY9CLbmbRZcHYSuToeTWUdlREySfi/YaqQystZ56J+I+CvlEXe9UVLHgXD1TnvWl35VTGkvqj2Miy8vSs6wi+G3v5XxgJvkZ9ziQtaBOUmC00gaUhtEi3H8qqRsqAj+nPC+bcRJI8FdUut5LCbrtjuRpcxiovnIxyF0wk/yMexSaKTOAUkoXFQWTykQpZF05Tu+iZpnKibOJ9CNBXExyuFwxNVkeN/mflHc3G15nzJgx4zLGkSR5InosgJ8AsI+wOcgLmfk9REQAXgvg2QDOAPgnzPz+oxZ2xoyj4lx8j5956n9LUbI7ZRjcks3ynp5k7qSilvDPRbI/CSq1GkdV17wawL9m5l8homfH738HwLMQtvx7FICvBfC6+H/GjEsOU3nJZ1zeOImEfS44qrqGAdwrfr4aebPu6wG8gQPeDeDeRPSgIz5rxowZM2YcEkeV5P85gLcR0Y8iTBhPiscfAuCz6rrPxWN/Wt+AiF4A4AUA8NCHPvSIxZkxY8aMGRpbSZ6I3gnggY1TNwD4egAvYeafJ6LnAviPAJ4OoBWj3VRkMvPrAbweCAnKdiz3jBkzZszYAVtJnpmfPnWOiN4A4Hvi17cAkD3wPgfgWnXpNciqnBkzZsyYcYFwVJ38bQD+dvz8NACfiJ/fCuBbKOCJAO5g5pGqZsaMGTNmnF8cVSf/HQBeS0QdgANE3TqAX0Zwn7wVwYXy2474nBkzZsyYcQ443kTzCwAABVBJREFUEskz828D+OrGcQbwoqPce8aMGTNmHB1zxOuMGTNmXMaYSX7GjBkzLmPMJD9jxowZlzFmkp8xY8aMyxgzyc+YMWPGZYyZ5GfMmDHjMgZxaxf5iwQi+nMAf3yxy3EOuD+AL1zsQpwj5rJfeFyq5Qbmsl8M7FLu/4mZH9A6caJI/lIFEb2XmR9/sctxLpjLfuFxqZYbmMt+MXDUcs/qmhkzZsy4jDGT/IwZM2ZcxphJ/njw+otdgCNgLvuFx6VabmAu+8XAkco96+RnzJgx4zLGLMnPmDFjxmWMmeRnzJgx4zLGTPLnACL6NBH9ARF9gIjeG4/dl4jeQUSfiP/vc7HLCQBEdDMRfZ6IblHHmmWNm7z8OyK6lYg+RESPO2HlfiUR/Uls9w8Q0bPVuVfEcn+MiJ55cUqdynItEf0GEX2EiD5MRN8Tj5/odt9Q7hPf7kS0T0TvIaIPxrL/63j84UT0e7HN30REy3h8L36/NZ5/2Akr908R0adUmz82Hj98X2Hm+e+QfwA+DeD+1bFXA/j++Pn7AfzwxS5nLMtTADwOwC3byoqw0cuvIOzR+0QAv3fCyv1KAC9tXPsYAB8EsAfg4QA+CcBexLI/CMDj4uerAHw8lvFEt/uGcp/4do9td2X8vADwe7Et3wzg+fH4TwD4rvj5hQB+In5+PoA3nbBy/xSAb25cf+i+Mkvyx4frAfx0/PzTAP7+RSxLAjP/FoC/rA5PlfV6AG/ggHcDuDcRPejClLTERLmncD2ANzLzipk/hbAj2RPOW+G2gJn/lJnfHz/fCeAjAB6CE97uG8o9hRPT7rHt7opfF/GPEbYl/S/xeN3m8i7+C4CvJyK6QMVN2FDuKRy6r8wkf25gAG8novcRkWx5+KUc97GN/7/kopVuO6bK+hAAn1XXfQ6bB/nFwIvjMvVmpRI7seWOaoC/gSChXTLtXpUbuATanYgsEX0AwOcBvANhZXE7Mw+N8qWyx/N3ALjfhS1xQF1uZpY2f1Vs85uIaC8eO3SbzyR/bvg6Zn4cgGcBeBERPeViF+iY0JJkTpKP7esAPBLAYwH8KYDXxOMnstxEdCWAnwfwz5n5i5subRy7aOVvlPuSaHdmdsz8WADXIKwo/lrrsvj/xJS9LjcRfSWAVwD4CgBfA+C+AF4eLz90uWeSPwcw823x/+cB/FeEDvVnsmyK/z9/8Uq4FVNl/RyAa9V11wC47QKXbRLM/GdxQHgAP4msGjhx5SaiBQJR/iwz/0I8fOLbvVXuS6ndAYCZbwfwmwg663sTkexlrcuXyh7PX43d1YPnBarcfzeqzpiZVwD+bxyhzWeSPySI6Aoiuko+A3gGgFsAvBXAt8bLvhXAL16cEu6EqbK+FcC3RAv+EwHcIeqFk4BK9/hNCO0OhHI/P3pMPBzAowC850KXTxB1u/8RwEeY+UZ16kS3+1S5L4V2J6IHENG94+dTAJ6OYFP4DQDfHC+r21zexTcD+HWOls0LiYlyf1QJA4RgR9Btfri+cjEsypfyH4BHIHgUfBDAhwHcEI/fD8CvAfhE/H/fi13WWK7/jLDE7hGkgG+fKivCUvDfI+gy/wDA409YuX8mlutDsbM/SF1/Qyz3xwA86yK3+ZMRltAfAvCB+Pfsk97uG8p94tsdwFcB+H9jGW8B8H/E449AmHhuBfAWAHvx+H78fms8/4gTVu5fj21+C4D/hOyBc+i+Mqc1mDFjxozLGLO6ZsaMGTMuY8wkP2PGjBmXMWaSnzFjxozLGDPJz5gxY8ZljJnkZ8yYMeMyxkzyM2bMmHEZYyb5GTNmzLiM8f8DsWp/r33CfkoAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# plot map of surface temperature\n",
    "plt.pcolormesh(lon,lat,T[0,0,:,:])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.QuadMesh at 0x32558a278>"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAD8CAYAAACYebj1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO2de7RkVX3nP7+qurdfPLoBjUrDiGvaZDAagy2amSSDiNDtZNlmxsygM5FEs5gwmsnDRECSmKgkoolmXONjMZFEMzoNMRqIARpQTGatKNA+UBCBHjDS4oNH8+zue29V/eaPs/euc3adqlv31r23btPfz1q1Tp199tnnd3ad2r+zf/u398/cHSGEEAKgMWkBhBBCrB6kFIQQQiSkFIQQQiSkFIQQQiSkFIQQQiSkFIQQQiRWnVIws21mdqeZ7TGzCyYtjxBCHE7YapqnYGZN4C7gFcBe4Bbgte7+zYkKJoQQhwmrradwKrDH3e9x91lgJ7BjwjIJIcRhQ2vSAmQcD9xX2t8LvCTPZGbnAucCbNiw4UU/9mM/tqxC3fXle5a8TDOLX6rbSqbsS54n7pbLCN89pgW17424H46n9LC1/jRKx8r7WOhdNhwLaU3rFttGsZ0K20ZMD+e0rFPs08XwivgN6nut3XCj3SBIsV+ti0H93VRFeCo/yhT3y3kKefrTmylvVWYLuTxKUN0AFqucVqjAWG435OqEbTuc1A75Ot6gE7575d5727Y3U95ia+kHm69O0u84CgPK9LS1vnv3WDcxPZXRS6977EcSx3vlpOumC1fT6Vb3rVu69WwbHo3BlVciiR7yNmfCs/Hkgb68W045af4Cx2Sh7dTj7HvQ3Z+Wp682pVD3iPT9PO5+KXApwNatW3337t3LLRevaP6nkgDdwRkjNrwTZs3iz2zTU5V94hZSA26p5W1U8zSy9EYT1k4XIk4XP203btcV287a4tz2uka2NdrrQkOzriiuvZ5wTthfV/wU3TXdcA3HpkOjv6YNwJo1cwCsC9tnHfEYAOtbxX6rUSiFI1sHOaI1W6SVFAXAEa2DABzdPFBJ399dA8AP5o5KjWEjHFvTKK6/tlFcZyqUGZmyDmvCsQ2NmSKvzVXOWWuzlXOjIpiyNlOhtdhgsfxwLHs88yfjoDfohMd6Y7j34xprQxnFPTzRLeT5Qac4/ki4z4e6G3i8sy6UUzwnB7tTlbp4NBx/aPYIAB5rr+Vgp8jTzRVJ2G8E2XtKcnTlMKjMdrfR+x4UVLtbbOe6xX3OdYrtbNi2Ow26IY8NkGFQurvRCee2O0EhtsN2rijf58L/Yib812YL+ZoHGjSLR4zw06dts3gEaMV2PVcedTKGH73RKTIddW/xe07fUXq/3QvXfP9DgwtZBJU2CXiJLUzx3OCf+ue69NWmFPYCJ5T2NwP3T0iWCtd3Lu/7EVLDX6ckBimOcI6HBoDwEDpz4bBh00XDHht/j3+MeE47exMJr1u2dg20i8bR1hUNTyO8Ppnnr1OtkB6K6DTSm87UE0V5nTW5kggNQauQqzsN7fXFOd4MjdWaorE6cFRxP48+Xpw8PV3IddyRTxZylf5lG6aKSjhq+mDltrpeVaxHTxX3va45m47FBvex9tpK3tiIx15B7Lkc0TyYGv14bH2juP6zpvYV8mRKA2/RCb/nXPj9Hu8W13ukswGAOw8+E4CH5jZUrrfG5pgKyiAqpSh7bOCPbu0v6ipo4dseexYAex46jmazKOdpG54o5A91FRvc6WZR9rHTT4ZrtEuNcvidvPqulfZjcuz4jaAcUp6Y1bJ0oBVV4yjG6WaXuXYz9SJyJTAo3cxTrzR7pHu9qryHEuql610sKqPsb9FJ/4cg3swI9xDqoNssvjx+YvEfOHq2+B1b33lghEJG58zp1xWXjS+Hqa1p1p9QpvyyOlefZbUphVuALWZ2EvBd4GzgdZMVqcf1ncsBeEXjFyrp1hzyY4QfIeWJtoTu4D+gzxUNaFICsY89yCkgHj84A+FcC+ear6/kaUQFEp4jCw1HY6ZLd03407TCn2h/sZ2bCefOhf0jepee6gTFsSH+m0IZjxSKrbup+Fe1G8V1HnqiaDQ3bdhPKwjx5FyRNzZmx6zZH2QtjsdGdN9scS8zrRZHtaoKJCcqi/gX6Lgx022xvzOVegAbWoVsT1A08A/OHVlcv1U0sEc3CzmObB5gOrQSndDA3D+3CYC/+8ELAPjW3ccX1zsQlH6zuMbRJz7K5qMeBeDR2eI6e793TJEnvMnaVDdUXVCwTxZ/y+YTzdT1eLyxqXJ/qSPUrTb4nQ1dpp9VyP/sYx8GeookNty5klgMuXIoK4W+HsmAnknTnNlOk0bD6XSGy5QrBy/dg2U6rhHrMdvSjPtGtxVkmYvm1lBGeGC6oWWMPYdRzEmRoOvpToXC2m26jzzGWet+EYBdB/5q9MJqsEYm7DBl0BhSr4eCUnD3tpm9GdhFcaeXufvtExarj+u7fw30um8eGviy5k7fW+GNfCo8Kc15Xp/KymLQDzpb/JqptxFx7ymOkAcr3q6tHUwiwbzUXV9sbTY0dkdM9XoN0+FNs3jhYerJ8AcK3fRGu6ccwosxrdi7iC/sQUl09xWFtI8pGuDY8D26fx2bNoTGP/zjDrSLOnrQi0KPW/tk7e0/0V6TGp4jp4Yrh6hQHmtPJ7NK5OG5QskcGRRMVDTfmz262FJsuzSYC2+WT3aK+7n1geIt8KG7jgNgzZOZHTtsDzx8DN9aXzToHhqiZmyIpmIjFcwbIb01U2rwOtVtsJKlRqyRHeeRBn5/ody+vbbYzhxXHJz+keJZODb01tYGk17sbaxvzSVFvRCTUsyfK4NcYcQGPb7hd7qNnqLKGvthZqO6awN4KLcbFGVqM6NyaPbGwjy++MS/bD7GFtKjcmhE5VAnU2yjM5EPHlf8x6Zum+uN5c2105v+dbOfHFzoEJIlYUEnjf4isKqUAoC7Xw1cPWk5FkXoxlmzWR0bgJ6JJ5qPYt5sH+hTBmlMIaavC2XHHkWZbEA5mpPS+EOwqzY64Q90RNHINfcD66qNZiu8vcUeBAdDgzHX6zk0Zvt7DwCNYN+dOzpcJ/QcZo+K9wAPPl40/kevLxrjNc1C1plO8Vh+f3/RqB23NvYceo3N/k61dxFNS5HYQMVGvKwQej2Q0N0PpqcHZorr/WB/cTMPBflmDkzRORje3h8vtlOPFueuDTopiVbTQLQORrNbdZC/OxX2W9WTUy+u3VMCCY9lhjzZe0GZ+JY79UQwI323uK9908U2VA2d9aGR3tBhw7FFXZ+w6RGg1MsYYfwhPzaoxxDd4MsNf8/kE+z+C1BKyTEg/seCErBG7IHF/0TsMfQURPxN4phCGruOHfqoFOLxYWJl5rjZo0JPcOPRvZe0/EVuMUzHrsjCFPeorDaX1EOK6zuXJ5PSoYxPjWCLFEIcFqy6nsJTAe967y2omb2550QPo2AXdDq9N4D8nCzdplrV9Gajv4cSewhpJK36ptLYX7wJdtdO00xjF6GMNWE7E3oIzWh/DWaAGZh+IpiaHgs9hjAYHT2XmgdDjyH0JNqh9zEz22TqiOLa+54oMsdB1alWIWMcoH1ytnilXTs1F8Rq0wpvgbPBgyf2BKYb1furt58XMu0P5qp7Hj62uM7eoqew9sHieKsY2+XI2fIYTHXrmQ06Un6Djzbm3CSR0ltWSaemt5Ekjz2H3JWyNPCbvnar58SeQzf0MprxhfPJOA7UYPahwmR251HFDzb9tKIH9qyNxbjIprXFfqzniqkoDZdVewj5bxC9kWY7TQ7MFELkZqFoAorjA7kJqo6+Hkk4t5vGFoKoTe+Zc1q9cYZyHsvMSCl9hBf95A4e5Dj4nONY88BDRVqr96CcdcQ5AOx64mPzF1rCgjMH3fgQ5k4kgUX2JKQUlpDe2AIlj4D6zlgyCWW2PrNGoRjK5Wampj6i8uiUTFD5GEYrUyCRMNZgs+2eh1KUJeqg9AcKjVdUDk3DwoMZG5zpx8PYwlTVgykqi7kjCnlmNzaYOzq4Kk6HCwVXV5rVgdehxP9fbDSil+5UcV+tYC83g04YEzn4ZGF6av6w2K59oChk0+Ph3OCFFe/JOjVyWKyLYrevQS/5u/cUR3xBoHJuNONERerRi8Z6eZMvfiwiVlVm1vAm/Y9cVBxRSWTKwoOHTXeqpEAOhvkPjxXK4dtHFl5k3zm6qJRNm4pxiePWF9uWdZObaiS6nkZz4Gy72EbPqQMzUz031uBG2o2D7+G3b8ZB+DTgFeRrdGk1s2OBPG9uTqJppXGGqkKOjb5nv1EaeO6yoEFngANPa7E2eANyYPgY2Eg0srYj/m+73Wq+BY4LpeIXKZYocaibkXx6av5MQojDAvUUlgPvQjN49wxyV01uZdafnka5RjMjVc8Nx+biK2Ro8OMrUOxBxHxmWBiwTmaH6L4a3pDjxDfiW2RId7PSrOde7wF63hrRk2Y6eiftC2+GDxqzRwbTUvBg6gZzVc+lL5QdXQnrzCrRzJB54Vg0s8RO1xwE6xNrg2zRB70VegaNMKU4mQiGvmhFU0XYS7O9+89Nb6Pt6htsvJ84aNwMg/ZxXkh3upcn0psoVb1O7m9fR7yveL+xrNiD6azp1ZtnLpnNMAGsE+ad7HswTLBbt7HIsKZLI5j9LPxejfSGXvVGiqahbqeBB3OidxuVPB5k7bablRuMvYB2w5nJegLR/BhNUZ3Q64hl+Fyv4pN5KPYYYu8seoZlL+OppzDHwOcin/4Raa8tpTSW8T18WNl5L2IIUgpLSOwtnDl19oQlEUKsJuJs5m2bfmXCksyPzEdCCCES6iksE3FWsrUGVHFaBS6OaHVTepq7EAfGutXJOQMHnOuI8xSiXSGas3LPpYb1Bppng+y5LSKaQabjOk3gaaJQ9G4KtxHNG+1qGc1gb5k6YGlQOvrrd6arg34xvZtPNDL6FjRrZGajaPpK5qRuKa0d08J+7smTp1My41jVvOCZCS+ZHbrltGBOiaaSzHIYfRKmutGrK5iR1vZ7Ls1nuqgjeUwN8JyJZrNGp3ftZMIL86SiOTCuGdTdH2e/xxnc0A0OA3FSXjtfUDGXy2u+53lz85iVkq2SJT3KFiZXxnuZytY3sk7JEyzeZ+YokMqu8UKyES0x+bOxZOT/y1EmpqWJtfMPPqunIIZPhRdCHFZIKQghhEhIKQghhEhoTGG58G4xkyqNGQRb3lScfTPAVdWsN2O5EQyesYwwPhAXwovL/9aaf3LbYZzJfHCmUqZP9R6BNIYQXxXCOEdjJizAN9WsHPaSrNHuH9eZyV1Ue5O6osus0wq2Xw926lYcQ2hWz+2bNVxjQ03jAKUxhMr1vJQnHsuN0ZEa+3bvvqpleJo9O8AfEbDoghnXOIr3lTJk143Z5kq29NJ4SrXweI3euXm5fYu81YYtyWUOMoTHJXdVTeM+8TFulGbypmdigOzl+81cPtOkssz9OK63Vfl9B8wy740rhW27mm5eGpvIxhbCQrk98jq00uOXm/YXN1ds1SGlsAxcN7czfT8UXNCEECvPtfv+fNIi1CLzkRBCiIR6CitEWrxu0KzDsmtq7J/m6xfF0J3xnLgcb7tmCe2cfEG8YEayULa3Gnia5plNj41LErezBbiaTbwVYyyEQ3EWa7yHbP3/hFlvEbK4plI41IyX6VRNM71lpGv66ZnrYpSxYioqmZIqeQdRdoOM7qxJ/ihsZjbK3FrxkjtpDGEZTU7h3EEusd61fjfS3GUyM4NQ466bm2gsc5vty1+50Sxvya0TqpbL3G2UTLa4erlPZfkAS663VZNhItZNySSUzETt/mOV9DoX5Di7OnN3jufEWeVeMo9BMJ+N8HdbdYzgihqRUlhmYhdx+4+cN2FJBhMbdiHE8rJaTUZl1BoIIYRIqKewUnRHtFk0Gr3ZzXksVs+mUsbVTeO2vHR2LGOQaSmakWbmsBnwdWt6oV7nW7QrTb7uYnOlxfEAYq8jebxUTSTJZOBeMklEU0/oyseZv5nHT88cMdh7ps+jqFzt7tU8o1K+XPSqGnS9mN7tT0/L7McFB8k8beLs8JLnVDJvxNvJZwlnVeEV4cqJC7zvfBVm5tkvmWTymfDx2YjxG2LI1s5Ur6C0OF80I2WPYB4bwuZKZqJoNoqmrWQ+ykyIubmMkqkuPrbJFJV5UNX9JQZ4Ia0qFmA2iqinIAqFIIQQSCkIIYQoIaUghBAioTGFlSYPnTc0MEY2w3g+mo1e+Y3oY5etwhrHGOL+7Cw2Owvr1+GN7HEYIFu0GXuHNA6RbMzRbTW6PSZX1RhmseRymLlV9o0UVD1ie/k6PtrKkFCxqfbZ1AetzJlRngGbu+mSzZLOyy6PZTQz234nxZbMBEw/u6Wd3phCdaXVQS6qQ++rz30329YdS+kD3Hq9xrU2nRvuIa5WGmay2/qSzT6rv3wWdL/bbmmcIR9L6FR/k3wsofxb5T9rlCMNW8VwyDE44ZAVUj0bH1lxFjF+UId6CivENQ98hGse+Mikxahn/bpJSyCEWCVIKQghhEjIfLTSjOKaGs02+UJ3lqV3sz4vlGwteb8/myVdnrX85H5oNpOrpE9nATliep3JZsSAH8nkFLvyWG/2c+jEezTFpJPCtlHdLWSpvcy8pqDi5AH7Q87NXWoHmVdq3V6tmrW3oGA8JTMJpRnQHkxIZfNR7u6ZXcNstDqoFFISNf89+8xH1fSKSWbQrPK8jDRL22gH99S40F5aaHBAgJ5yjOo8lvZgk1eWXnKHTqRHPVxwrvpQRLdX71n0+kx3A01ty8UylT9WT8HMfsHMbjezrpltzY5daGZ7zOxOMzurlL4tpO0xswtK6SeZ2U1mdreZXW5m0+PIJhbAoBVbhRCHHeOaj24D/j3wj+VEMzsZOBt4HrAN+JCZNc2sCXwQ2A6cDLw25AW4BHi/u28B9gFvHFM2IYQQC2Qs85G73wGlLlePHcBOd58B7jWzPcCp4dged78nnLcT2GFmdwCnA68LeT4G/AHw4XHkW9UM80LqTX2tPzfOYI7Emcxz7X43jmRqyq5bvkaKzVwE4rVG4XLhaY37RvWcIa8SFu8nmkSGzT5OpoIBZqR4apS9vIAauRmlKlutqatPgOru0DjHo/bULTMF1ZQR7y8tIpjNhk4eRk1LMRjyheLSDNvcvGI+MOZC7wL1ctV6H2XnDE3PzUU13j4AjZCvdbD4tNcanWAXKMdCrr2HUpmNLHZGn0mrT45Q73O9733mqmzxwjibPsWgaAz32iqu2++h9eI3vI9bLvstlgLff6D4Ej0J44KV0WFkvhUJ5mG5BpqPB+4r7e8NaYPSjwUecfd2ll6LmZ1rZrvNbPcDDzywpIIvN6vaC0mIFaa9dqEDIGK5mbenYGY3AM+oOXSRu1856LSaNKdeCfmQ/LW4+6XApQBbt25dzSuPCCHEIcW8SsHdz1hEuXuBE0r7m4H7w/e69AeBjWbWCr2Fcn4hhBArxHK5pF4FfNLM3gc8C9gC3EzRI9hiZicB36UYjH6du7uZ3Qi8BtgJnAMM6oUsKz+7473p+z9e+TvLf8E6t7Jkwx/RuheD8DSnoR0MrXNz9eUPclmF6tgEYOHxSO6CSebSeMF8tvto703jBqGIbm8nrY7ayGzS+Szdiptn9bq5HTd3Zx1pjGFA2UWxXnts2DhEr7woS/1sZO9zPY6Fe1qtM91Yso8PKKvskTrqLY8wljD4eOk3y/Nmtv483TqeZiMnG/6geM41MvfNWM5Wxq0bSwBoznlaQTVf3TYFfGpnD1AoLK5aC4NdTpd9RnO+8nF0945jDSGgl1mjz63d5+aPEDSuS+rPm9le4KeAvzezXYWMfjtwBfBN4FrgTe7eCb2ANwO7gDuAK0JegPOB3wqD0scCHx1HttXONT/4MNf84Kk7ji6EODQZ1/voM8BnBhy7GLi4Jv1q4Oqa9HvoeSgJIYSYAJrRnGFz9V3CZSN3TW02e+ag+c6pM4m0Qv87lpG61GEbYjMPLTd3dYuBSOJ0wri4XafX3R4663mxRNMENSabgbM5g8xxEbaqH2v9ZUYyAQ2aqZ2VsZDbD+c2sljUlanPKa5wbi6K51SvW10Qr16YvlnZQ2WsZhpmXho0S3gQZbfSWHA3CJcvgth3bslclbuXlmc9AzRnw/5cN2xrZjRHwv+/mQXoYUPPjBTrPF9ob756PfWX3gfAzX+5NK6pfSTX8sJOttiWTGsfCSGESEgpCCGESMh8lJE8GMw4bft7+MI1b12ZC8f1h0Yxv0QZh8wWTuXEPLFruS6sQDYTzEjluM7xnJgWF6hvBW+G4JHhU6XrdqIsWRl9M1LrpvhWZ+kOzVs+7DVmhfnqbQHmlZyRYhREcg+jEc4ZGKvBy55YmfdW7rlVMwk+LZ433wznGubzQkr56swwC7BbNOZCHOo0+7dq9hsoqvebjdIM507V2yh6GvXNgK6UV/V8i6a8xmyROa0hsN4qHkjl8vqqOXlBeb+H2SpGPYV5OG37e5a1fHkhrRwjjR2IFSVvYMXyc+1jf8G1j/3FwONSCkIIIRJSCkIIIRIaUxjEcgfIyGktYExhHKJtc00IPhtWRk0zoQHWhmMz4VgMLhzGGpJde6rZswEPsAn3uarW3d6QlUUXTBo7sGy7xNeZT4yym+Jw79k++32cPW1OZXwhlVc5J+StGTOJ4w0pZ/4TjFMPi/l71A0r5SupponE2dhCTfCfRorJHMuojiVYitUcn98aoQfNSs4mR/fGFhppRdduFhhoEEUwoENniTb1FERVIQghDmukFIQQQiRkPhrAyEFVlhr3ki0g+rplsZkX496WL4QXy5iehmmgU+otRFMWoZ+cL8AViwQ8+upZNH/FGdrRJbaUuSSHN8rfh5iWYHjPOzMX9S2AV7ruomYdL5RutX6HPkf5bOhANb5zNqN4BBF65VVdXAfd/zgzuodfPz9Qvb51SfeXQk93MjflIWbcPnNR3M9cU20BZqP8eAr8FJIbbe8d6wyvk4m1IWMipbBKuObb7wdg+7N/c+Uv3pH5SIhDhehOuu2oX17UefMh85EQQoiElIIQQoiElIIQQoiElIIQQoiElIIQQoiElIIQQoiEXFIzesvd9tJeduYl3Hjd+RMQJujsPJpaZFgEtoUQl+3OideN2zjzOUWLc2ymXciVlrcIj9SA141epDBL8xPSPIX5XlHqlhPIq2TQHA73tGzEiriP1/nGD2Cpp030lzdgCe08qPvQ5b7r76d2ToINKbeSD+ibAjN6vZV/azdozfSey3JR5SWsx6brNLKlN3rzPwY/ewkzfubn/wSA//uZ3x5fnmVAPYVVRpyvcMgwX+jQ1cRyryslJsY/XfGWSYvwlOEQ+kcLIYRYbqQUhBBCJKQUhBBCJKQUhBBCJMZSCmb2XjP7lpl93cw+Y2YbS8cuNLM9ZnanmZ1VSt8W0vaY2QWl9JPM7CYzu9vMLjez6XFkE0IIsXDG7SlcD/y4u78AuAu4EMDMTgbOBp4HbAM+ZGZNM2sCHwS2AycDrw15AS4B3u/uW4B9wBvHlE0IIcQCGUspuPt17h49jb8EbA7fdwA73X3G3e8F9gCnhs8ed7/H3WeBncAOMzPgdOBT4fyPAa8eRzYhhBALZyknr70BuDx8P55CSUT2hjSA+7L0lwDHAo+UFEw5fx9mdi5wLsCJJ544tuAVsvktjfY8AViXg6kQuSYFjllC//pBwXa6pRvPAsWkIC3TUwwkxgPKYzOE2Wo+VZ0gV8RMri/KB80nGDDHbhTMKQVOKcnwFGNYYJe++80m2JnZ4Lwp04DdyjNaLXfUsmoZOt9sQHzl/C87LFDPfLF2YgzsGGyn5jZTGf3BtovkUj37GM9wn2zd4W2TNRb3zj+vUjCzG4Bn1By6yN2vDHkuopib+Il4Wk1+p75n4kPy1+LulwKXAmzduvUQjW+0iuiqCsWhzz/83e8AcNr290xYktXFqMF1IvMqBXc/Y9hxMzsH+Dng5e5JJe8FTihl2wzcH77XpT8IbDSzVugtlPMLIYRYIcb1PtoGnA+8yt33lw5dBZxtZmvM7CRgC3AzcAuwJXgaTVMMRl8VlMmNwGvC+ecAV44jmxBCiIUz7pjC/wTWANcHe+SX3P1X3f12M7sC+CaFWelN7t4BMLM3A7soLMSXufvtoazzgZ1m9i7gq8BHx5RNCCHEAhlLKbj7vxxy7GLg4pr0q4Gra9LvofBOEkIIMSE0o1kIIURCSkEIIURCSkEIIURCSkEIIURCSkEIIURCSkEIIURCSkEIIURCSkEIIURCSkEIIURCSkEIIURiKeMpPKWwsOBrc/8cANte+PsAXPu1d6ysIKPEUgiyelw/fZCqLy+/HtbOt263uEbd8tldL47lcQia1TgL5PEPzDD3Ii5CPDfEpbCwIL2H+/JWA8OL/W5xESOcu4hXFh9pkf4g5mpfMXyJ4z0UMSWGlFuKO2AMiIMwINaF5zELStkse4YHxVeov17YeDXPOL/dWOd2ga7jFkSrkTnFc8jqxBxoO501S/su7jMz9emLLE9KYQQaB9vzZ1pCrrnrEgC2/6sLF37yKAqh79iAx2cUhVSjEGBIoJwS3goBeFIwH0Y+VyyCBbQSYwUgWsS5Sxnw6AvXvLWyf9q2S5aucEqyLkDmsiL6x6t+Z0nk2PXExwA4c+rsJSkvIvOREEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIxFhKwczeaWZfN7Ovmdl1ZvaskG5m9gEz2xOOn1I65xwzuzt8zimlv8jMvhHO+YCZlsoUQoiVZtyewnvd/QXu/kLgs8Dvh/TtwJbwORf4MICZHQO8HXgJcCrwdjPbFM75cMgbz9s2pmxCCCEWyFjxFNz9sdLuBnortu8APu5F5I0vmdlGM3smcBpwvbs/DGBm1wPbzOwLwFHu/sWQ/nHg1cA148i3GKxT3EIrBNdpPnoAOsOCEUyQPHbBMBU/4FgMzGOzHWiO+I4Qq6NRWiR+UMeuEwP0hJPiYvTNXv7YKewF7wnldguEIWIAABIDSURBVKtxFurIYy/YokOL9AIrjcLEYj6M8BP1yTZg/f++GAal8/qP5XLEYEk15aXnsnrugq43iOWq9vjTjxKAyKvbQY9c+XnyUWKTrBLGDrJjZhcDrwceBV4Wko8H7itl2xvShqXvrUkfdM1zKXoVnHjiiePdwCrmmjv+GFhksJ2FXuvOd7P95Lct+3WEWGm+cO35A4+97KzBAXhu3HX+kgXouXHXYBnG5bq5nUA12E5MWwzzvneY2Q1mdlvNZweAu1/k7icAnwDeHE+rKWpQIMBh6bW4+6XuvtXdtz7taU+b7xaEEEKMyLw9BXc/Y8SyPgn8PcWYwV7ghNKxzcD9If20LP0LIX1zTX4hhBAryLjeR1tKu68CvhW+XwW8PnghvRR41N2/B+wCzjSzTWGA+UxgVzj2uJm9NHgdvR64chzZhBBCLJxxxxTebWY/SjH0+M/Ar4b0q4FXAnuA/cAvA7j7w2b2TuCWkO8dcdAZOA/4S2AdxQDzig8yCyHE4c643kf/YUC6A28acOwy4LKa9N3Aj48jjxBCiPHQjGYhhBAJKQUhhBAJKQUhhBAJKQUhhBAJKQUhhBAJKQUhhBAJKQUhhBAJKQUhhBAJKQUhhBAJKQUhhBCJseMpPOXZ92ixPerIycpBLyDOSqnycrAWc8enmsVOuxR0yL02wE4MMJLKyALYWAfohnJicJ0QVMdSAJcYzaRafjkgi+F4s4F16wOa2DIFSBovmE/4kgf1ye+zfB/x/rxaJ7WBa+LzEbYeV6YfFOelrqwBQXWSbDUBe/L7youoC9zjuWgDZEx14UVgH2/W51suUmCdtveC59T8jt3pFRZsGZBSyPj85y6o7G9/xn8D4Jq7libYxmrmmm/+UWV/2/Muqhzb9hO/V3vetV9/18BjK8HnP3cBL3/ZH1f2Iy8/7Y/qThGHMP90+Vv4qdf96ZKUNV/wmxigZ1gwnkp51y1fMJ35GCewThmZj4QQQiSkFIQQQiSkFIQQQiSkFIQQQiSkFIQQQiSkFIQQQiSkFIQQQiSkFIQQQiSkFIQQQiSkFIQQQiSkFIQQQiSkFIQQQiSWRCmY2W+bmZvZcWHfzOwDZrbHzL5uZqeU8p5jZneHzzml9BeZ2TfCOR8wq1l6UwghxLIytlIwsxOAVwDfKSVvB7aEz7nAh0PeY4C3Ay8BTgXebmabwjkfDnnjedvGlU0IIcTCWIqewvuBt0JlgfkdwMe94EvARjN7JnAWcL27P+zu+4DrgW3h2FHu/kV3d+DjwKuXQDYhhBALYKx4Cmb2KuC77n5rZu05HrivtL83pA1L31uTPui651L0KjjxxBPHuINDjHFUeIw1sxSvATEoTqc/0Ex3esAjlV+39Lx4bimMQVo6MW9IjjGGOqWsDTjjZy/GAG/231xzppe57zoDsDz4zQqR5Iv3OyRgSwpYY4MD6FgnBrupv590vYb3p/UVNkKdLLTaSpeqDRYEeCsGYFqeYEmin3mVgpndADyj5tBFwNuAM+tOq0nzRaTX4u6XApcCbN26dVn/wdd8/0PLWfxoMtxRBJDZ9vyL5sm5tFx7+8XV/VvfWcjxgt/tzxuOnXnqO5ZfsAF87sYLJ3bt5WDY/RzuwYO++Mm3rOj15gvG81RiXqXg7mfUpZvZ84GTgNhL2Ax8xcxOpXjTP6GUfTNwf0g/LUv/QkjfXJNfCCHECrJoY4K7f8Pdn+7uz3b3Z1M07Ke4+/eBq4DXBy+klwKPuvv3gF3AmWa2KQwwnwnsCsceN7OXBq+j1wNXjnlvQgghFshyxWi+GnglsAfYD/wygLs/bGbvBG4J+d7h7g+H7+cBfwmsA64JHyGEECvIkimF0FuI3x1404B8lwGX1aTvBn58qeQRQgixcDSjWQghREJKQQghREJKQQghREJKQQghREJKQQghREJKQQghREJKQQghREJKQQghREJKQQghREJKQQghRGK51j4SE8A6Du0QPyBt28U2X4++Ed4HNqxbXqHia8eQWAYxfkFcy9+6+Wro1fgK3mzQODiXjnbWT80vxoG5efOMGm8Blj/mQmfD9OJO7Hpf/IS+UAgxXkV+biXGxeBj5ZNHqbNUV1GOeP2aKixft7NOzdMkMJ9QQJGlYuvWrb579+5JizFRYpyFcZRCjNmwoOs+r7huHncBSnEVRlAKkaQURngmba4XPGfX7j+YN/+2F/7+yNcfheVWCtd+bXFxKeriLAxSCv0ZV5dS+NwX3jZv2WLxmNmX3X1rni7zkRBCiISUghBCiISUghBCiISUghBCiISUghBCiISUghBCiISUghBCiISUghBCiISUghBCiISUghBCiISUghBCiISUghBCiMRYSsHM/sDMvmtmXwufV5aOXWhme8zsTjM7q5S+LaTtMbMLSuknmdlNZna3mV1uZotcJlIIIcRiWYqewvvd/YXhczWAmZ0MnA08D9gGfMjMmmbWBD4IbAdOBl4b8gJcEsraAuwD3rgEsgkhhFgAy2U+2gHsdPcZd78X2AOcGj573P0ed58FdgI7zMyA04FPhfM/Brx6mWQTQggxgKWIYvFmM3s9sBt4i7vvA44HvlTKszekAdyXpb8EOBZ4xN3bNfnFchDjKzz+5OLLmGriZpz1k28HYNdX/7Avi3XSIvopzVsNrF2N7+CtZR7eGiVGw0JiJByKcUhiHIPuoOO9A9YXJyEEQgo/kzeLL+WAPvPFVmjMDbpwTzaao8e0EMvDvErBzG4AnlFz6CLgw8A7KZ6YdwJ/CryBmqBOIU/dP9+H5B8k07nAuQAnnnjiEOkPD679Rn+Qm5ztm//7CkjS47qbi6A2Z724X1Fc98Xf60u/7ou/tyJyPdUpB6apC7iz0nz+cxfw8tMXHsBJTI55lYK7nzFKQWb2v4DPht29wAmlw5uB+8P3uvQHgY1m1gq9hXL+OpkuBS6FIvLaKPIJIYSYn3G9j55Z2v154Lbw/SrgbDNbY2YnAVuAm4FbgC3B02iaYjD6Ki9igt4IvCacfw5w5TiyCSGEWDjjjim8x8xeSGHq+TbwXwHc/XYzuwL4JtAG3uTuHQAzezOwC2gCl7n77aGs84GdZvYu4KvAR8eUTQghxAIZSym4+y8OOXYx0GfsDm6rV9ek30PhnSSEEGJCaEazEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIhJSCEEKIxFIE2RFCHAK09u1f0vLmjjtiQfmbjx+sJgwIVNTeuG6xIoklwPxQjCBVYuvWrb579+5Ji3HYEqOuQX3ktdXItp/oD+hz7a3vXFgZL/jdgceu/fq7hh6PeVaabc+/aEnLGyW4U+X6eb0PaHsmUTeHI2b2ZXffmqfLfCSEECIhpSCEECIhpSCEECIhpSCEECIhpSCEECIhpSCEECIhpSCEECIhpSCEECIhpSCEECIhpSCEECIhpSCEECIhpSCEECIxtlIws18zszvN7HYze08p/UIz2xOOnVVK3xbS9pjZBaX0k8zsJjO728wuN7PpcWUTQgixMMZSCmb2MmAH8AJ3fx7wJyH9ZOBs4HnANuBDZtY0sybwQWA7cDLw2pAX4BLg/e6+BdgHvHEc2YQQQiyccXsK5wHvdvcZAHf/YUjfAex09xl3vxfYA5waPnvc/R53nwV2AjvMzIDTgU+F8z8GvHpM2YQQQiyQcZXCc4GfCWaffzCzF4f044H7Svn2hrRB6ccCj7h7O0uvxczONbPdZrb7gQceGPMWhBBCROaNvGZmNwDPqDl0UTh/E/BS4MXAFWb2HMBq8jv1SsiH5K/F3S8FLoUiyM4w+YXIWWhAndoy5gkEsxoDxSw0KM6SX38J6l0sP/MqBXc/Y9AxMzsP+LQX4dtuNrMucBzFm/4JpaybgfvD97r0B4GNZtYKvYVyfiGEECvEuOajv6UYC8DMngtMUzTwVwFnm9kaMzsJ2ALcDNwCbAmeRtMUg9FXBaVyI/CaUO45wJVjyiaEEGKBzNtTmIfLgMvM7DZgFjgnNPC3m9kVwDeBNvAmd+8AmNmbgV1AE7jM3W8PZZ0P7DSzdwFfBT46pmxCCCEWyFhKIXgQ/ZcBxy4G+oyY7n41cHVN+j0U3klCCCEmhGY0CyGESEgpCCGESEgpCCGESEgpCCGESEgpCCGESEgpCCGESEgpCCGESEgpCCGESEgpCCGESEgpCCGESEgpCCGESFixft2hi5k9APzzmMUcR7G662pjNcq1GmWC1SnXapQJVqdckml0lkquf+HuT8sTD3mlsBSY2W533zppOXJWo1yrUSZYnXKtRplgdcolmUZnueWS+UgIIURCSkEIIURCSqHg0kkLMIDVKNdqlAlWp1yrUSZYnXJJptFZVrk0piCEECKhnoIQQoiElIIQQojEYa0UzOyFZvYlM/uame02s1NDupnZB8xsj5l93cxOWWG5fs3M7jSz283sPaX0C4NMd5rZWSspU0mG3zYzN7Pjwv7E6srM3mtm3wrX/YyZbSwdm2hdmdm2cO09ZnbBSl8/yHCCmd1oZneEZ+nXQ/oxZna9md0dtpsmIFvTzL5qZp8N+yeZ2U1BpsvNbHoCMm00s0+FZ+oOM/upSdeVmf1m+O1uM7P/Y2Zrl72u3P2w/QDXAdvD91cCXyh9vwYw4KXATSso08uAG4A1Yf/pYXsycCuwBjgJ+H9Ac4Xr6wRgF8VkweNWQV2dCbTC90uAS1ZDXQHNcM3nANNBlpNX8rcKcjwTOCV8PxK4K9TNe4ALQvoFsd5WWLbfAj4JfDbsXwGcHb5/BDhvAjJ9DPiV8H0a2DjJugKOB+4F1pXq6JeWu64O654C4MBR4fvRwP3h+w7g417wJWCjmT1zhWQ6D3i3u88AuPsPSzLtdPcZd78X2AOcukIyRd4PvJWi3iITqyt3v87d22H3S8DmkkyTrKtTgT3ufo+7zwI7g0wrirt/z92/Er4/DtxB0dDsoGgACdtXr6RcZrYZ+HfAn4d9A04HPjVBmY4Cfhb4KIC7z7r7I0y4roAWsM7MWsB64Hssc10d7krhN4D3mtl9wJ8AF4b044H7Svn2hrSV4LnAz4Tu4T+Y2YtXgUyY2auA77r7rdmhicpV4g0UPRaYvEyTvn4fZvZs4CeBm4AfcffvQaE4gKevsDh/RvFy0Q37xwKPlBT8JOrrOcADwF8Es9afm9kGJlhX7v5dinbpOxTK4FHgyyxzXbWWsrDViJndADyj5tBFwMuB33T3vzGz/0jxlnAGhSkkZ8l8d+eRqQVsojDFvBi4wsyes9wyjSDX2yjMNX2nLadcw2Ry9ytDnouANvCJlZBpBCZ9/QpmdgTwN8BvuPtjxYv5xGT5OeCH7v5lMzstJtdkXen6agGnAL/m7jeZ2f+gMBdNjDB+sYPCBPoI8NfA9pqsS1pXT3ml4O5nDDpmZh8Hfj3s/jWhO0uhfU8oZd1Mz7S03DKdB3zaC4PhzWbWpVgAa1llGiaXmT2f4sG8NTQom4GvhIH5idVVkO0c4OeAl4c6Y7llGoFJXz9hZlMUCuET7v7pkPwDM3umu38vmPp+OLiEJeffAK8ys1cCaynMt39GYXZshTfgSdTXXmCvu98U9j9FoRQmWVdnAPe6+wMAZvZp4F+zzHV1uJuP7gf+bfh+OnB3+H4V8PrgWfNS4NHYhVwB/jbIgpk9l2LA68Eg09lmtsbMTgK2ADevhEDu/g13f7q7P9vdn03xBzrF3b/PBOvKzLYB5wOvcvf9pUMTq6vALcCW4CUyDZwdZFpRgq3+o8Ad7v6+0qGrgHPC93OAK1dKJne/0N03h+fobODz7v6fgRuB10xCpiDX94H7zOxHQ9LLgW8ywbqiMBu91MzWh98yyrS8dbVSI+mr8QP8NIWN7lYKW+uLQroBH6TwIPkGsHUFZZoG/jdwG/AV4PTSsYuCTHcSvKYmVG/fpud9NMm62kNhu/9a+HxktdQVhVfWXUGGiyb0O/00hWnh66U6eiWFDf9zFC9BnwOOmZB8p9HzPnoOheLeQ9FrXzMBeV4I7A719bcUZtyJ1hXwh8C3QnvwVxQedctaV1rmQgghROJwNx8JIYQoIaUghBAiIaUghBAiIaUghBAiIaUghBAiIaUghBAiIaUghBAi8f8BRT1DG+/omV8AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# plot meridional temperature section through Atlantic\n",
    "plt.pcolormesh(lat[:,0],-depth,T[0,:,:,330])"
   ]
  }
 ],
 "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.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
