Recipes/Cross-slope_section.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "tags": []
   },
   "source": [
    "# Cross-slope section\n",
    "\n",
    "We compute a cross-slope section of fields through gridded data in ACCESS-OM2-01 using the  `metpy.interpolate.cross_section` function; see https://unidata.github.io/MetPy/latest/examples/cross_section.html#sphx-glr-examples-cross-section-py."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Load modules"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import cosima_cookbook as cc\n",
    "from dask.distributed import Client\n",
    "import numpy as np\n",
    "import xarray as xr\n",
    "import xgcm\n",
    "\n",
    "from metpy.interpolate import cross_section\n",
    "\n",
    "import cartopy.crs as ccrs\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib.path as mpath\n",
    "import cmocean as cm\n",
    "from matplotlib.colors import BoundaryNorm\n",
    "from matplotlib.ticker import (MultipleLocator, AutoMinorLocator)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Start a cluster with multiple cores"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "    <div style=\"width: 24px; height: 24px; background-color: #e1e1e1; border: 3px solid #9D9D9D; border-radius: 5px; position: absolute;\"> </div>\n",
       "    <div style=\"margin-left: 48px;\">\n",
       "        <h3 style=\"margin-bottom: 0px;\">Client</h3>\n",
       "        <p style=\"color: #9D9D9D; margin-bottom: 0px;\">Client-3e5e5d98-6737-11ef-a3f8-0000018ffe80</p>\n",
       "        <table style=\"width: 100%; text-align: left;\">\n",
       "\n",
       "        <tr>\n",
       "        \n",
       "            <td style=\"text-align: left;\"><strong>Connection method:</strong> Cluster object</td>\n",
       "            <td style=\"text-align: left;\"><strong>Cluster type:</strong> distributed.LocalCluster</td>\n",
       "        \n",
       "        </tr>\n",
       "\n",
       "        \n",
       "            <tr>\n",
       "                <td style=\"text-align: left;\">\n",
       "                    <strong>Dashboard: </strong> <a href=\"/proxy/8787/status\" target=\"_blank\">/proxy/8787/status</a>\n",
       "                </td>\n",
       "                <td style=\"text-align: left;\"></td>\n",
       "            </tr>\n",
       "        \n",
       "\n",
       "        </table>\n",
       "\n",
       "        \n",
       "            <button style=\"margin-bottom: 12px;\" data-commandlinker-command=\"dask:populate-and-launch-layout\" data-commandlinker-args='{\"url\": \"/proxy/8787/status\" }'>\n",
       "                Launch dashboard in JupyterLab\n",
       "            </button>\n",
       "        \n",
       "\n",
       "        \n",
       "            <details>\n",
       "            <summary style=\"margin-bottom: 20px;\"><h3 style=\"display: inline;\">Cluster Info</h3></summary>\n",
       "            <div class=\"jp-RenderedHTMLCommon jp-RenderedHTML jp-mod-trusted jp-OutputArea-output\">\n",
       "    <div style=\"width: 24px; height: 24px; background-color: #e1e1e1; border: 3px solid #9D9D9D; border-radius: 5px; position: absolute;\">\n",
       "    </div>\n",
       "    <div style=\"margin-left: 48px;\">\n",
       "        <h3 style=\"margin-bottom: 0px; margin-top: 0px;\">LocalCluster</h3>\n",
       "        <p style=\"color: #9D9D9D; margin-bottom: 0px;\">30678711</p>\n",
       "        <table style=\"width: 100%; text-align: left;\">\n",
       "            <tr>\n",
       "                <td style=\"text-align: left;\">\n",
       "                    <strong>Dashboard:</strong> <a href=\"/proxy/8787/status\" target=\"_blank\">/proxy/8787/status</a>\n",
       "                </td>\n",
       "                <td style=\"text-align: left;\">\n",
       "                    <strong>Workers:</strong> 8\n",
       "                </td>\n",
       "            </tr>\n",
       "            <tr>\n",
       "                <td style=\"text-align: left;\">\n",
       "                    <strong>Total threads:</strong> 48\n",
       "                </td>\n",
       "                <td style=\"text-align: left;\">\n",
       "                    <strong>Total memory:</strong> 188.56 GiB\n",
       "                </td>\n",
       "            </tr>\n",
       "            \n",
       "            <tr>\n",
       "    <td style=\"text-align: left;\"><strong>Status:</strong> running</td>\n",
       "    <td style=\"text-align: left;\"><strong>Using processes:</strong> True</td>\n",
       "</tr>\n",
       "\n",
       "            \n",
       "        </table>\n",
       "\n",
       "        <details>\n",
       "            <summary style=\"margin-bottom: 20px;\">\n",
       "                <h3 style=\"display: inline;\">Scheduler Info</h3>\n",
       "            </summary>\n",
       "\n",
       "            <div style=\"\">\n",
       "    <div>\n",
       "        <div style=\"width: 24px; height: 24px; background-color: #FFF7E5; border: 3px solid #FF6132; border-radius: 5px; position: absolute;\"> </div>\n",
       "        <div style=\"margin-left: 48px;\">\n",
       "            <h3 style=\"margin-bottom: 0px;\">Scheduler</h3>\n",
       "            <p style=\"color: #9D9D9D; margin-bottom: 0px;\">Scheduler-78bfd87c-be80-495c-af15-80018501222b</p>\n",
       "            <table style=\"width: 100%; text-align: left;\">\n",
       "                <tr>\n",
       "                    <td style=\"text-align: left;\">\n",
       "                        <strong>Comm:</strong> tcp://127.0.0.1:35469\n",
       "                    </td>\n",
       "                    <td style=\"text-align: left;\">\n",
       "                        <strong>Workers:</strong> 8\n",
       "                    </td>\n",
       "                </tr>\n",
       "                <tr>\n",
       "                    <td style=\"text-align: left;\">\n",
       "                        <strong>Dashboard:</strong> <a href=\"/proxy/8787/status\" target=\"_blank\">/proxy/8787/status</a>\n",
       "                    </td>\n",
       "                    <td style=\"text-align: left;\">\n",
       "                        <strong>Total threads:</strong> 48\n",
       "                    </td>\n",
       "                </tr>\n",
       "                <tr>\n",
       "                    <td style=\"text-align: left;\">\n",
       "                        <strong>Started:</strong> Just now\n",
       "                    </td>\n",
       "                    <td style=\"text-align: left;\">\n",
       "                        <strong>Total memory:</strong> 188.56 GiB\n",
       "                    </td>\n",
       "                </tr>\n",
       "            </table>\n",
       "        </div>\n",
       "    </div>\n",
       "\n",
       "    <details style=\"margin-left: 48px;\">\n",
       "        <summary style=\"margin-bottom: 20px;\">\n",
       "            <h3 style=\"display: inline;\">Workers</h3>\n",
       "        </summary>\n",
       "\n",
       "        \n",
       "        <div style=\"margin-bottom: 20px;\">\n",
       "            <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n",
       "            <div style=\"margin-left: 48px;\">\n",
       "            <details>\n",
       "                <summary>\n",
       "                    <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 0</h4>\n",
       "                </summary>\n",
       "                <table style=\"width: 100%; text-align: left;\">\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Comm: </strong> tcp://127.0.0.1:36959\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Total threads: </strong> 6\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Dashboard: </strong> <a href=\"/proxy/42865/status\" target=\"_blank\">/proxy/42865/status</a>\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Memory: </strong> 23.57 GiB\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Nanny: </strong> tcp://127.0.0.1:38621\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\"></td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td colspan=\"2\" style=\"text-align: left;\">\n",
       "                            <strong>Local directory: </strong> /jobfs/123969626.gadi-pbs/dask-scratch-space/worker-91nversz\n",
       "                        </td>\n",
       "                    </tr>\n",
       "\n",
       "                    \n",
       "\n",
       "                    \n",
       "\n",
       "                </table>\n",
       "            </details>\n",
       "            </div>\n",
       "        </div>\n",
       "        \n",
       "        <div style=\"margin-bottom: 20px;\">\n",
       "            <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n",
       "            <div style=\"margin-left: 48px;\">\n",
       "            <details>\n",
       "                <summary>\n",
       "                    <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 1</h4>\n",
       "                </summary>\n",
       "                <table style=\"width: 100%; text-align: left;\">\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Comm: </strong> tcp://127.0.0.1:41035\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Total threads: </strong> 6\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Dashboard: </strong> <a href=\"/proxy/34597/status\" target=\"_blank\">/proxy/34597/status</a>\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Memory: </strong> 23.57 GiB\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Nanny: </strong> tcp://127.0.0.1:35061\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\"></td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td colspan=\"2\" style=\"text-align: left;\">\n",
       "                            <strong>Local directory: </strong> /jobfs/123969626.gadi-pbs/dask-scratch-space/worker-ucnwqn4f\n",
       "                        </td>\n",
       "                    </tr>\n",
       "\n",
       "                    \n",
       "\n",
       "                    \n",
       "\n",
       "                </table>\n",
       "            </details>\n",
       "            </div>\n",
       "        </div>\n",
       "        \n",
       "        <div style=\"margin-bottom: 20px;\">\n",
       "            <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n",
       "            <div style=\"margin-left: 48px;\">\n",
       "            <details>\n",
       "                <summary>\n",
       "                    <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 2</h4>\n",
       "                </summary>\n",
       "                <table style=\"width: 100%; text-align: left;\">\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Comm: </strong> tcp://127.0.0.1:37911\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Total threads: </strong> 6\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Dashboard: </strong> <a href=\"/proxy/34817/status\" target=\"_blank\">/proxy/34817/status</a>\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Memory: </strong> 23.57 GiB\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Nanny: </strong> tcp://127.0.0.1:38633\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\"></td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td colspan=\"2\" style=\"text-align: left;\">\n",
       "                            <strong>Local directory: </strong> /jobfs/123969626.gadi-pbs/dask-scratch-space/worker-593hvh83\n",
       "                        </td>\n",
       "                    </tr>\n",
       "\n",
       "                    \n",
       "\n",
       "                    \n",
       "\n",
       "                </table>\n",
       "            </details>\n",
       "            </div>\n",
       "        </div>\n",
       "        \n",
       "        <div style=\"margin-bottom: 20px;\">\n",
       "            <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n",
       "            <div style=\"margin-left: 48px;\">\n",
       "            <details>\n",
       "                <summary>\n",
       "                    <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 3</h4>\n",
       "                </summary>\n",
       "                <table style=\"width: 100%; text-align: left;\">\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Comm: </strong> tcp://127.0.0.1:39315\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Total threads: </strong> 6\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Dashboard: </strong> <a href=\"/proxy/41551/status\" target=\"_blank\">/proxy/41551/status</a>\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Memory: </strong> 23.57 GiB\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Nanny: </strong> tcp://127.0.0.1:35795\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\"></td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td colspan=\"2\" style=\"text-align: left;\">\n",
       "                            <strong>Local directory: </strong> /jobfs/123969626.gadi-pbs/dask-scratch-space/worker-kvaoqbbl\n",
       "                        </td>\n",
       "                    </tr>\n",
       "\n",
       "                    \n",
       "\n",
       "                    \n",
       "\n",
       "                </table>\n",
       "            </details>\n",
       "            </div>\n",
       "        </div>\n",
       "        \n",
       "        <div style=\"margin-bottom: 20px;\">\n",
       "            <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n",
       "            <div style=\"margin-left: 48px;\">\n",
       "            <details>\n",
       "                <summary>\n",
       "                    <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 4</h4>\n",
       "                </summary>\n",
       "                <table style=\"width: 100%; text-align: left;\">\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Comm: </strong> tcp://127.0.0.1:42405\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Total threads: </strong> 6\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Dashboard: </strong> <a href=\"/proxy/34047/status\" target=\"_blank\">/proxy/34047/status</a>\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Memory: </strong> 23.57 GiB\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Nanny: </strong> tcp://127.0.0.1:35757\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\"></td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td colspan=\"2\" style=\"text-align: left;\">\n",
       "                            <strong>Local directory: </strong> /jobfs/123969626.gadi-pbs/dask-scratch-space/worker-zzvxvcp2\n",
       "                        </td>\n",
       "                    </tr>\n",
       "\n",
       "                    \n",
       "\n",
       "                    \n",
       "\n",
       "                </table>\n",
       "            </details>\n",
       "            </div>\n",
       "        </div>\n",
       "        \n",
       "        <div style=\"margin-bottom: 20px;\">\n",
       "            <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n",
       "            <div style=\"margin-left: 48px;\">\n",
       "            <details>\n",
       "                <summary>\n",
       "                    <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 5</h4>\n",
       "                </summary>\n",
       "                <table style=\"width: 100%; text-align: left;\">\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Comm: </strong> tcp://127.0.0.1:36415\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Total threads: </strong> 6\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Dashboard: </strong> <a href=\"/proxy/41767/status\" target=\"_blank\">/proxy/41767/status</a>\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Memory: </strong> 23.57 GiB\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Nanny: </strong> tcp://127.0.0.1:33299\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\"></td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td colspan=\"2\" style=\"text-align: left;\">\n",
       "                            <strong>Local directory: </strong> /jobfs/123969626.gadi-pbs/dask-scratch-space/worker-t98rw48y\n",
       "                        </td>\n",
       "                    </tr>\n",
       "\n",
       "                    \n",
       "\n",
       "                    \n",
       "\n",
       "                </table>\n",
       "            </details>\n",
       "            </div>\n",
       "        </div>\n",
       "        \n",
       "        <div style=\"margin-bottom: 20px;\">\n",
       "            <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n",
       "            <div style=\"margin-left: 48px;\">\n",
       "            <details>\n",
       "                <summary>\n",
       "                    <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 6</h4>\n",
       "                </summary>\n",
       "                <table style=\"width: 100%; text-align: left;\">\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Comm: </strong> tcp://127.0.0.1:33347\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Total threads: </strong> 6\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Dashboard: </strong> <a href=\"/proxy/45505/status\" target=\"_blank\">/proxy/45505/status</a>\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Memory: </strong> 23.57 GiB\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Nanny: </strong> tcp://127.0.0.1:43813\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\"></td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td colspan=\"2\" style=\"text-align: left;\">\n",
       "                            <strong>Local directory: </strong> /jobfs/123969626.gadi-pbs/dask-scratch-space/worker-ssawjy8_\n",
       "                        </td>\n",
       "                    </tr>\n",
       "\n",
       "                    \n",
       "\n",
       "                    \n",
       "\n",
       "                </table>\n",
       "            </details>\n",
       "            </div>\n",
       "        </div>\n",
       "        \n",
       "        <div style=\"margin-bottom: 20px;\">\n",
       "            <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n",
       "            <div style=\"margin-left: 48px;\">\n",
       "            <details>\n",
       "                <summary>\n",
       "                    <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 7</h4>\n",
       "                </summary>\n",
       "                <table style=\"width: 100%; text-align: left;\">\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Comm: </strong> tcp://127.0.0.1:46213\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Total threads: </strong> 6\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Dashboard: </strong> <a href=\"/proxy/46165/status\" target=\"_blank\">/proxy/46165/status</a>\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Memory: </strong> 23.57 GiB\n",
       "                        </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td style=\"text-align: left;\">\n",
       "                            <strong>Nanny: </strong> tcp://127.0.0.1:36317\n",
       "                        </td>\n",
       "                        <td style=\"text-align: left;\"></td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <td colspan=\"2\" style=\"text-align: left;\">\n",
       "                            <strong>Local directory: </strong> /jobfs/123969626.gadi-pbs/dask-scratch-space/worker-2bknnw3x\n",
       "                        </td>\n",
       "                    </tr>\n",
       "\n",
       "                    \n",
       "\n",
       "                    \n",
       "\n",
       "                </table>\n",
       "            </details>\n",
       "            </div>\n",
       "        </div>\n",
       "        \n",
       "\n",
       "    </details>\n",
       "</div>\n",
       "\n",
       "        </details>\n",
       "    </div>\n",
       "</div>\n",
       "            </details>\n",
       "        \n",
       "\n",
       "    </div>\n",
       "</div>"
      ],
      "text/plain": [
       "<Client: 'tcp://127.0.0.1:35469' processes=8 threads=48, memory=188.56 GiB>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "client = Client(threads_per_worker = 1)\n",
    "client"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Nominate a database from which to load the data and define an experiment"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load database\n",
    "session = cc.database.create_session()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "depth = 3000 # metres"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Pick shelf coordinates"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "shelf_coord = (-62,   -60 + 100)\n",
    "deep_coord  = (-56.5, -59 + 100)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Load velocity and bathymetry data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Select data in the southern part of the Southern Ocean\n",
    "lat_slice = slice(-80, -59)\n",
    "\n",
    "# We use a time-mean from some output from the RYF run\n",
    "expt = '01deg_jra55v13_ryf9091'\n",
    "start_time = '1950-01-31 00:00:00'\n",
    "end_time = '1959-12-31 00:00:00'\n",
    "\n",
    "# Import bathymetry\n",
    "hu = cc.querying.getvar(expt, 'hu', session, n=1)\n",
    "hu = hu.sel(yu_ocean=lat_slice)\n",
    "hu = hu.load()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "## importing density - ignore warnings\n",
    "pot_rho_2 = cc.querying.getvar(expt, 'pot_rho_2', session,\n",
    "                               start_time=start_time, end_time=end_time, frequency='1 monthly',\n",
    "                               attrs={'cell_methods': 'time: mean'},\n",
    "                               chunks={})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "#importing potential temperature - ignore warnings\n",
    "pot_temp =cc.querying.getvar(expt=expt, variable='pot_temp', \n",
    "                             session=session, frequency='1 monthly',\n",
    "                             attrs={'cell_methods': 'time: mean'},\n",
    "                             start_time=start_time, end_time=end_time,\n",
    "                             chunks={})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<xgcm.Grid>\n",
       "Y Axis (not periodic, boundary=None):\n",
       "  * center   yt_ocean --> outer\n",
       "  * outer    yu_ocean --> center\n",
       "X Axis (periodic, boundary=None):\n",
       "  * center   xt_ocean --> right\n",
       "  * right    xu_ocean --> center"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Load model grid information with depth\n",
    "path_to_folder = '/g/data/ik11/outputs/access-om2-01/01deg_jra55v13_ryf9091/output000/ocean/'\n",
    "grid = xr.open_mfdataset(path_to_folder+'ocean_grid.nc', combine='by_coords')\n",
    "\n",
    "# Give information on the grid: location of u (momentum) and t (tracer) points on B-grid \n",
    "ds = xr.merge([pot_temp.sel(yt_ocean=lat_slice).sel(st_ocean=slice(0, depth)), grid.sel(yu_ocean=lat_slice).sel(yt_ocean=lat_slice)])\n",
    "ds.coords['xt_ocean'].attrs.update(axis='X')\n",
    "ds.coords['xu_ocean'].attrs.update(axis='X', c_grid_axis_shift=0.5)\n",
    "ds.coords['yt_ocean'].attrs.update(axis='Y')\n",
    "ds.coords['yu_ocean'].attrs.update(axis='Y', c_grid_axis_shift=0.5)\n",
    "\n",
    "grid_depth = xgcm.Grid(ds, periodic=['X'])\n",
    "grid_depth"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Take time mean for density and temperature"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "pot_rho_2 = pot_rho_2.mean(dim = 'time')\n",
    "pot_rho_2 = pot_rho_2.sel(yt_ocean = lat_slice).sel(st_ocean=slice(0, depth))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "pot_temp = pot_temp.mean(dim = 'time')\n",
    "pot_temp = pot_temp.sel(yt_ocean = lat_slice).sel(st_ocean=slice(0, depth))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we make the cross sections using the cross_section function of metpy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "#choosing number of steps in cross section\n",
    "step_no = 400\n",
    "\n",
    "# Create dataset\n",
    "ds_pot_rho_2 = xr.Dataset({\"pot_rho_2\": pot_rho_2, \"lat\": pot_rho_2.yt_ocean, \"lon\": pot_rho_2.xt_ocean})\n",
    "\n",
    "# Interpolate to xu_ocean and yu_ocean\n",
    "# Rename coordinate names\n",
    "ds_pot_rho_2 = ds_pot_rho_2.rename({'xt_ocean': 'x', 'yt_ocean': 'y'})\n",
    "\n",
    "# Convert latitude from ACCESS-OM2 default range of [-280, 80] to [-180, 180] which is what metpy expects.\n",
    "ds_pot_rho_2['x'] = ds_pot_rho_2['x'] + 100\n",
    "ds_pot_rho_2['x'].attrs = pot_rho_2['xt_ocean'].attrs\n",
    "\n",
    "# MetPy parsing\n",
    "pot_rho_2_parsed = ds_pot_rho_2.metpy.parse_cf('pot_rho_2', coordinates={'y': 'y', 'x': 'x'})\n",
    "\n",
    "pot_rho_2_section = cross_section(pot_rho_2_parsed,\n",
    "                                  start=(shelf_coord[0], shelf_coord[1]),\n",
    "                                  end=(deep_coord[0], deep_coord[1]),\n",
    "                                  steps = step_no,\n",
    "                                  interp_type='linear')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create dataset\n",
    "ds_pot_temp= xr.Dataset({\"pot_temp\": pot_temp, \"lat\": pot_temp.yt_ocean, \"lon\": pot_temp.xt_ocean})\n",
    "\n",
    "# Interpolate to xu_ocean and yu_ocean\n",
    "# Rename coordinate names\n",
    "ds_pot_temp = ds_pot_temp.rename({'xt_ocean': 'x', 'yt_ocean': 'y'})\n",
    "\n",
    "# Convert longitude from ACCESS-OM2 range of [-280, 80] to [-180, 180] which is what metpy expects.\n",
    "ds_pot_temp['x'] = ds_pot_temp['x'] + 100\n",
    "ds_pot_temp['x'].attrs = pot_temp['xt_ocean'].attrs\n",
    "\n",
    "# MetPy parsing\n",
    "pot_temp_parsed = ds_pot_temp.metpy.parse_cf('pot_temp', coordinates={'y': 'y', 'x': 'x'})\n",
    "pot_temp_section = cross_section(pot_temp_parsed,\n",
    "                                 start = (shelf_coord[0], shelf_coord[1]),\n",
    "                                 end = (deep_coord[0], deep_coord[1]),\n",
    "                                 steps = step_no,\n",
    "                                 interp_type = 'linear')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "# load variables\n",
    "pot_temp_section.load()\n",
    "pot_rho_2_section.load();"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Finally calculate the distance along the transect (for plotting purposes)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Define number of points you want to interpolate\n",
    "step_no = 400\n",
    "\n",
    "# Radius of the Earth\n",
    "Rearth = 6371 # km\n",
    "\n",
    "# Difference between points in lat/lon space\n",
    "dlon = deep_coord[1] - shelf_coord[1]\n",
    "dlat = deep_coord[0] - shelf_coord[0]\n",
    "\n",
    "# Calculate distance in km between the two end points\n",
    "distance_endpoints = Rearth * np.deg2rad(np.sqrt(dlat**2 + (dlon * np.cos(np.deg2rad(np.mean([shelf_coord[0], deep_coord[0]]))))**2))\n",
    "\n",
    "# Create array with length of step_no\n",
    "distance_in_km = np.linspace(0, distance_endpoints, step_no)\n",
    "\n",
    "# Repeat by the number of depth levels\n",
    "distance_in_km = np.tile(distance_in_km, (len(pot_temp_section.st_ocean), 1))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Plotting cross-slope section of potential temperature"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAk4AAALrCAYAAAAIiAioAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOydd5gc1ZW33wqde/JoNBqFUZYQkhASoEA0GRsMxtjGERwBs+tlbX9rnHHAeDH2YnvBNmtjjAnGGExOAiTAILKQhEBCEkoz0mhy7FhV5/ujuns6T4/SKNT7PP1U1617b93qrvCrc889VxERwcHBwcHBwcHBYUjUkW6Ag4ODg4ODg8PBgiOcHBwcHBwcHBxKxBFODg4ODg4ODg4l4ggnBwcHBwcHB4cScYSTg4ODg4ODg0OJOMLJwcHBwcHBwaFEHOHk4ODg4ODg4FAijnBycHBwcHBwcCgRRzg5ODg4ODg4OJSII5wOYlavXs3nP/95Jk2ahNfrJRgMMn/+fK6//no6OztHunl7lVNOOYVTTjllpJtxQHHXXXdx44035t2mKArXXHPNfm3PcHnmmWc45phjCAQCKIrCAw88UDT/Cy+8gMfjYevWrSXVLyL8+c9/5rjjjiMQCFBeXs78+fN58MEHU3l27tzJ9773PRYvXkxtbS3l5eUsWLCAW265BdM0Sz6W3/72t8ycOROPx8OkSZP40Y9+RDweL7l8Idrb2/mP//gPJk6ciMfjYfTo0Zxzzjk513d/fz9XXXUVDQ0NeL1e5s2bx9/+9reS99Pa2sqll15KbW0tfr+fxYsX88wzz2TkicfjTJkypeA5l82OHTu45ppreOutt0pux3AYzvlwzTXXoChKzsfr9ebNX+rvno98+1EUhZ///OfDPsZ0Jk6cmLfeyy+/PCfvq6++yllnnUVZWRnBYJAPfOADvPjiiyXva1+cD4cU4nBQcsstt4iu63LkkUfKTTfdJMuWLZOnnnpKfvazn8mkSZPkggsuGOkm7lVOPvlkOfnkk0e6GQcUH/rQh6SxsTHvthUrVsj27dv3b4OGgWVZUl1dLYsWLZKnn35aVqxYIZ2dnUXzz58/X6688sqS93HZZZeJx+ORq6++Wp5++ml54okn5Be/+IXcddddqTwPP/ywjB8/Xr773e/Ko48+Kk899ZT853/+p6iqKp///OdL2s9Pf/pTURRFvv3tb8uyZcvk+uuvF7fbLV/+8pdLbms+mpubZfLkyTJ9+nT54x//KM8995zcd9998m//9m+yc+fOjLxnnHGGVFZWyu9//3t59tln5Utf+pIAcueddw65n0gkIrNnz5Zx48bJHXfcIU899ZScf/75ouu6LF++PCPvbbfdJlVVVdLe3j5kva+99poA8uc//3lYx10Kwz0ffvjDHwogTzzxhKxYsSL1eeWVV3LyDud3zwcgF110UcZ+VqxYIc3NzcM+znQaGxvl+OOPz6n3/fffz8j36quvisfjkRNPPFH++c9/yv333y+LFi0Sj8cjL7300pD72Vfnw6GEI5wOQl566SXRNE3OPvtsiUQiOduj0ag8+OCDResIhUL7qnn7BEc45VJMOB3oNDU1CSD//d//XVL+xx57TABZt25dSfn/+c9/CiD33HNP0XydnZ0Si8Vy0q+88koBZNu2bUXLt7e3i9frla985SsZ6ddee60oiiJr164tqb35OP/882Xs2LFFBaWIyKOPPipAhiAUscVUQ0ODGIZRtPxNN90kQMZDNR6Py6xZs+S4447LyBuNRqW6ulquvfbaIdu/L4XTcM+HpHBqa2sbMm+pv3shgGEJ/FJpbGyUD33oQ0PmO+uss2T06NEyMDCQSuvt7ZXa2lpZsmTJkOX31flwKOEIp4OQc889V3RdH/KmniR5wd13330yb9488Xg88q1vfUtERNasWSMf/vCHpbKyUjwejxx11FFy2223ZZQ3TVN+8pOfyPTp08Xr9UpFRYXMmTNHbrzxxlSe1tZW+fKXvyzjxo0Tt9udukiXLl06ZPs2bdokn/jEJ2TMmDHidrulrq5OTj31VFm5cmUqTz7h1NHRIVdccYU0NDSIy+WSSZMmyXe+850cMZm8kf3+97+XadOmidvtliOOOELuvvvunLbs3LlTvvKVr8jYsWPF5XLJxIkT5ZprrpF4PD7kcTzzzDNy8sknS3V1tXi9Xhk/frxceOGFGTewaDQqP/nJT2TGjBmp3+nSSy+V1tbWnPruvPNOWbRokQQCAQkEAnLUUUfJH//4x9TvAeR80o/5hz/8YUZ9pfzXy5YtSz2Ev/Od78iYMWOkrKxMTjvttJIfUi+88IKceuqpEgwGxefzyeLFi+WRRx5JbU8+xNI/QwnA8847T4499tiS9i8icsopp8jEiRNLzp/NX/7yl5yHRz7uuOMOAWTFihUZ6Tt27BBgtx8omzdvFkVR5Jprrhky75e+9CUJBoM55+hdd90lgLz44otFy59++ukyY8aMnPSf/exnAkhTU1NG+hVXXCGNjY1iWVbBOpPnUfYn/Zx88MEHZdGiReLz+SQYDMrpp59ekkVEZPjnQ6nCaTi/eyFGWjgFg0H5xCc+kZN+4YUXCiA7duwoWn5fnA+HGo5wOsgwDEP8fr8sXLiw5DKNjY0yZswYmTx5stx6662ybNkyefXVV2XdunVSVlYmU6ZMkdtvv10effRR+eQnP5ljCbjuuutE0zT54Q9/KM8884w88cQTcuONN2bcXM466ywZNWqU3HLLLbJ8+XJ54IEH5Ac/+IH87W9/G7J9M2bMkKlTp8pf//rXlFn8G9/4hixbtiyVJ1s4hcNhmTt3rgQCAbnhhhvkqaeeku9///ui67p88IMfzKgfkPHjx8usWbPk7rvvloceekjOPvtsAeTee+9N5du5c6eMHz9eGhsb5Q9/+IM8/fTT8pOf/EQ8Ho9ceumlRY9h8+bN4vV65YwzzpAHHnhAli9fLnfeead89rOfla6uLhGxBejZZ58tgUBAfvSjH8nSpUvlj3/8o4wdO1ZmzZqVYQX8/ve/L4BceOGFcu+998pTTz0lv/rVr+T73/++iIisXbtWjj/+eKmvr88w26cfc/pDqtT/OvnAmzhxonz605+WRx99VO6++26ZMGGCTJs2bUjrxfLly8XlcsmCBQvknnvukQceeEDOPPNMURQldS5s375d7r//fgHk3//932XFihXy5ptvFqwzGo2Kz+eT//qv/yq67yTxeFw8Ho985CMfkV/+8pcyYcIEUVVVJk2aJL/4xS9KusFfcsklouv6kF0QV199tQDS39+fs622tlY++clPltTmbG6//XYB5JZbbpGLL75YAoGAeDweOfnkk3PExaJFi/KKiLffflsA+cMf/lB0X/X19fKxj30sJ/2RRx4RQJ588smM9HvuuUcAWb16dcE6e3p65M9//rMA8r3vfS91fia7j++8804B5Mwzz5QHHnhA7rnnHlmwYIG43W554YUXirZ3uOeDyKBwqq+vF1VVpa6uTj772c/K1q1bM/IN53cvBCBVVVXi9XrF7XbL/Pnz5dZbby25rYVobGyUsrIyCQaDouu6HHHEEXLDDTfkXJNut1s+97nP5ZRPXu/Z/2c2++J8ONRwhNNBRktLiwBy8cUXl1ymsbFRNE2T9evXZ6RffPHF4vF4cixX55xzjvj9funu7hYR28I1b968ovsIBoNy1VVXldymJO3t7QJkWK/ykS2cfv/73wsgf//73zPy/fd//7cA8tRTT6XSAPH5fNLS0pJKMwxDZs6cKVOnTk2lXXbZZRIMBnNupjfccIMARbtd/vGPfwggb731VsE8d999twBy3333ZaQnuzRuvvlmERF5//33RdM0+fSnP12wLpHiXXXZwqnU/zopnLLF59///ve8lpVsFi1aJHV1ddLX15dKMwwj5TORFC2bN28WQH7xi18UrU9E5JVXXhGgJBEuYgtgQMrLy2XcuHHyl7/8RZ555hm5/PLLBZDvfOc7Rcs/+eSToqqq/Od//ueQ+/ryl78sHo8n77bp06fLmWeeWVKbs7nuuutSx3D++efLE088Iffdd5/MnTtXvF6vrFq1KpV32rRpctZZZ+XUkbR6/exnPyu6L5fLJZdddllO+ksvvZS3C3DDhg0CyO9+97ui9RbqqjNNUxoaGmTOnDlimmYqva+vT+rq6obsThru+SBiC6Jrr71WHnvsMXn22Wfl5z//uVRXV8vo0aMzLCjD+d0L8alPfUruvPNOef755+Uf//iHnHPOOSkBuSd89atflVtvvVWee+45eeCBB+TTn/60APKZz3wmI9+8efNk+vTpGb9tPB6XyZMn5/0/s9lX58OhhCOcDjJ2VzgdffTROel1dXU5D0iRwTeIxx9/XEREfvzjH4uiKHLFFVfIE088IT09PTllTj31VKmsrJSf/OQnsmLFihy/EcuyJB6PZ3yS6VOmTJGxY8fKL3/5S3nzzTczLvgk2cLp4x//uAQCgRzrwa5duwRIdUWK2CLi3HPPzakz+RaafAseO3asnHfeeTntXLt2bYawycfGjRvF7XbLcccdJ7fddpts2rQpJ8+nP/1pqayslFgslrOP+vp6+fjHPy4iIn/4wx9K6iYajnAq9b9OCqff//73GfnWrVs35MOqv79fFEWRr371qznbkoL23XffFZHhCaekv9Kzzz6bkW6aZsZvmHzzbm5uTnUNZQu9Cy64QLxeb4awS+eNN96QiooKWbJkSV7/wWy+/OUvi9frzbtt+vTpeQVNOoZhZBxD8ty/9tprBZBZs2ZlWBR27Nghfr8/Q1RPmzZNzj777Jy6k8LpuuuuK9oGl8sll19+eU568kGZ3aXd09NTkhAoJJzeeecdAeT666/PKXPFFVeIqqoZ3dvZDPd8KMQrr7wiqqrK1772tVTacH734ZB0r8jXJZ9OofOhEP/2b/8mQIbF9k9/+pMAcsUVV0hTU5Ns27ZNvvjFL4qmaSUJzn11PhxKOOEIDjKSw0M3b948rHJjxozJSevo6Mib3tDQkNoO8O1vf5sbbriBl19+mXPOOYeamhpOO+00Xn/99VSZe+65h0suuYQ//vGPLF68mOrqaj73uc/R0tICwF/+8hdcLlfGB+yhu8888wxnnXUW119/PfPnz2fUqFF87Wtfo6+vr+DxdHR0UF9fj6IoGel1dXXoup5qe5L6+vqcOpJpyby7du3i4YcfzmnnkUceCdhDlAsxZcoUnn76aerq6rjyyiuZMmUKU6ZM4de//nUqz65du+ju7sbtdufso6WlJVV/W1sbAOPGjSu4v+FS6n+dpKamJmPd4/EAEA6HC+6jq6sLERnWfkohuc/soeM//vGPM37DKVOmAFBVVYWiKJSXl7No0aKMMueccw6RSIR33nknZz8rV67kjDPOYNq0aTz22GOpYy5GTU0NkUiEUCiUs62zs5Pq6uqi5adMmZJxDD/+8Y9T9QKcfvrpaJqWyj9mzBiOOuoo3nzzzYw25Ptdk0Pnh2rDcMsn/4di50IxkvsqdJ5YlkVXV1fB8sM9Hwpx3HHHMX36dF5++eVU2nB+9+Hwmc98BsMwMu6Z+Sh0PhSrF8g4hi984Qv8/Oc/569//Svjxo1jwoQJvPPOO3zzm98EYOzYsUXr3N/nw8GIPtINcBgemqZx2mmn8fjjj9PU1FTywzVbYIB9gezcuTMnfceOHYAt0gB0XefrX/86X//61+nu7ubpp5/mO9/5DmeddRbbt2/H7/dTW1vLjTfeyI033si2bdt46KGHuPrqq2ltbeWJJ57gvPPO47XXXsvbtsbGRv70pz8B8N577/H3v/+da665hlgsxu9///u8ZWpqanjllVcQkYxja21txTCMVNuTJAVcvrTkzbK2tpa5c+dy7bXX5t1n8uFfiBNPPJETTzwR0zR5/fXX+e1vf8tVV13F6NGjufjii6mtraWmpoYnnngib/mysjIARo0aBUBTUxPjx48vus9SKfW/3hOqqqpQVXWv7ydZJjuGzle+8hXOPffc1HpS6Ph8PqZNm5b3PxcRAFQ1851x5cqVnH766TQ2NvLUU09RUVFRUtvmzJkDwJo1a1i4cGEqPSmEZ8+eXbT8ww8/TDQaTa0nz7G5c+cWLCMiGe2fM2cOd999N4ZhoOuDt/Q1a9YADNmGOXPmpPKmU6h88n/Y3XMmeb0VOk9UVaWqqqpg+eGeD8XI/i2H87sPh0LnXTaFzofh1vutb32Lq666ig0bNlBWVkZjYyOXXXYZgUCABQsWFK1zf58PByUjZ+xy2F3SwxFEo9Gc7bFYTB566KHUeqHRGJ/85CfF6/XmxBf50Ic+lOH3ko8bb7xxSL+fCy64QEaNGlXKIeUwb968DIfX7K66ZHfW/fffn1HuF7/4hQAZo/ko4uM0ZcqUVNqXvvQlaWho2O1hyNl0d3cLIP/v//0/ERkcgfXyyy8XLbd582bRNE0++9nPFs134YUXSl1dXd5tZHXVlfpfJ7vq0p3mk22ihKHlixcvlvr6+gxHd9M0Zc6cObvt47R9+3YB5Ne//vWQeZN8+9vfzjui7MMf/rAEg8GM9q1cuVKqq6tl7ty5w45H09HRIV6vN6dr47rrrtujcASmacq4ceNk5syZGV1Gzc3N4vP55Itf/GIqLTk0P7sL5uyzzy4pHMHNN9+cc17G43E58sgj8w5CefHFFwUYMuTJ6tWr83Zxm6YpY8eOlXnz5mV0tff390tdXZ0cf/zxRevdnfMhHytWrBBVVTN8M4fzuw+HD37wg+JyuUoKhzAcrrjiiiF9K0VEtm7dKhUVFSX5oe6r8+FQwhFOBynJAJizZ8+Wm266SZYvXy5Lly6V66+/XqZOnZoRALOQcEqOtJo+fbrccccd8thjj6UcDtP9D84991y5+uqr5R//+Ic899xzcvvtt8vEiROlsbFRYrGYdHd3y9FHHy2/+MUv5OGHH5bly5fLL37xC/F6vfKpT32q6HGsWrVKTjzxRPnNb34jjz/+uDzzzDPy3e9+V1RVzXDiLTSqrqysTH71q1/J0qVL5Yc//KG4XK5hjapLf9js2LFDGhsbZebMmXLzzTfLM888I48++qjcdNNN8qEPfahoQMnf/e538rGPfUxuu+02efbZZ+Wxxx6Tiy66KGMUimEYcs4550h1dbX86Ec/kscff1yefvppue222+SSSy7JEIHJUXUXXXSR3HffffL000/Lb37zG/nBD36QypP00br55pvllVdekddeey3jmPONqhvqv95T4ZQcVbdw4UK599575cEHH5SzzjorY1Rden2lCCcRkcmTJw9rhFpHR4dMmDBBGhoa5E9/+pM8+eST8uUvf1kAueGGG1L51q1bJzU1NVJdXS0PP/xwTnDBdJ+U5cuXi6Zp8qMf/ShjX8kAmN/5zndS577H48kJgJk85ksuuaSkY7j33ntFURT50Ic+JI888ojcc889Mnv2bKmoqJCNGzdm5D3jjDOkqqpKbrnlFnn22WdTx3rHHXdk5PvCF74gmqbJli1bUmmRSESOPPJIGT9+vNx5552ydOlS+chHPpI34KGIyC9/+UvRNC01WrQQAwMD4vP55Pjjj5dly5bJa6+9lhLuyVF1H/zgB+XBBx+Uv//973LssceWNKpOZPjnw9y5c+X666+Xhx9+WJYuXSrXXnutVFZWSkNDQ87w/FJ/93znw/XXXy+XXnqp/PWvf5Vly5bJPffcI2eeeaYAOSEOhnM+3HnnnfLRj35Ubr31VnnmmWfkvvvuk4svvliAnBG/a9askWuuuUYeeeQRWbp0qdxwww1SW1srxxxzTI5v3/48Hw4lHOF0EPPWW2/JJZdcIhMmTBC32y2BQECOPvpo+cEPfpBxwy8W/2PNmjVy3nnnSUVFhbjdbjnqqKNyHo6//OUvZcmSJVJbWytut1smTJggX/ziF1MXWyQSkcsvv1zmzp0r5eXl4vP5ZMaMGfLDH/6wqJOniO3Mfemll8rMmTMlEAhIMBiUuXPnyv/8z/9kvPEViuN0+eWXy5gxY0TXdWlsbJRvf/vbBeM43XzzzTJlyhRxuVwyc+bMvFGV29ra5Gtf+5pMmjRJXC6XVFdXy4IFC+S73/1u3iHnSVasWCEf+chHpLGxUTwej9TU1MjJJ5+cYfkTsd/cbrjhBjnqqKPE6/VKMBiUmTNnymWXXSYbNmzIyHv77bfLsccem8p39NFHZ/w3nZ2dctFFF0llZaUoiiLpBuRs4SRS2n+9p8JJZDCOUyAQEJ/PJ4sWLZKHH344b32lCqfvf//7UlVVVZLDdpJt27bJxRdfLFVVVeJ2u2Xu3Lk5w8KTQ+YLfdKPN/nbZP+uIiK//vWvZfr06anr44c//GHOAIk1a9YIIFdffXXJx/DAAw+kzoGKigr58Ic/nNeK1dfXJ1/72tekvr4+daz54pRdcsklAsjmzZsz0ltaWuRzn/tcKgbZokWLCsZgO/HEE+W8884rqf133323zJw5U1wuV85v98ADD8jChQvF6/VKIBCQ0047bciYU0mGez5cfPHFMnXqVAkEAuJyuaSxsVEuv/zygjGNSvnd850PDz30kJxwwgkyatQo0XVdysrK5MQTT8z7XwznfFixYoWcdtppUl9fLy6XS/x+vxx77LFy88035ziQr1+/Xk466SSprq4Wt9stU6dOle9973t571/7+3w4VHCEk8MhT1I4ORy8NDc3i9vtHtYQ9AONm266SQKBQEaX8cHGxo0bRVGUjHAfI4FzPhwYHCjnw/7GGVXn4OBwwNPQ0MBVV13Ftddei2VZI92c3WLZsmV87WtfY/To0SPdlN3mpz/9KaeddhpnnHHGiLbDOR8ODA6U82F/44yqc3BwOCj43ve+h9/vp7m5ea+NNtyf3HvvvSPdhD3CMAymTJnCt7/97ZFuCuCcDyPNgXY+7E8UkcR4RgcHBwcHBwcHh6IcVl11N998M5MmTcLr9bJgwQJeeOGFkW6Sg4ODg4ODw0HEYSOc7rnnHq666iq++93vsnLlSk488UTOOecctm3bNtJNc3BwcHBwcDhIOGy66hYuXMj8+fP53e9+l0o74ogjuOCCC7juuutGsGUODg4ODg4OBwuHhXN4LBbjjTfe4Oqrr85IP/PMM3nppZfylolGoxmh7y3LorOzk5qamrzTlzg4ODg4ODiMPCJCX18fDQ0Nuz1NTjEOC+HU3t6OaZo5wz5Hjx6ddz4rgOuuu44f/ehH+6N5Dg4ODg4ODnuZ7du379XJ0pMcFsIpSbalSLImiE3n29/+Nl//+tdT6z09PUyYMIGzf/EQLl8AXbXLb3n+QVrXvUmkp8OeGX72cRxx5idw+8vo2PQ265bey+lf+ykA65/5Jy3vreL4z16Fr6IaRUxUTUcTgz9/7aNMPGoRc8+8kJYNa9m5fjUnfPqrVIwawz0/vIJpx57Esed9CoCX77+NltUvMn/2EbiCldSe8BFcwdxJMc1h9MKa1vB7bHenzN7EiIZpe+t52t5eQTzUl0pX3V7Kx0+jYuIsqiYegTtYPqx6NXVoi6JW4LxRi5QtVG92upbnBSk3T+Z6vv0m26hmtTU553t6mUJ5k1mSedO3a8ltiTQtK09ye/o1lspTYL/JY08/HD3xxji4TUnUm13HYJnkNhX7i2Bf6088eD+PP/wgv/79H3G53Tzz5OP87a+38d0fXcuUKVP5zS+vp2n7Nn7x65vo6+3m7dWraG5qorOjnffWvcspp57OhRd9DID+3h7uuuOvPPDP+ykrC7Jw8RL+3399C0VRUsesKomJWBMTxKb/fsmvSvI6lcx4REpyPT09+5rOKVPCNSkF4h4ppb+ZS77zP7t86gAz0yV9PfE9u77kYViJL+lHlbztiAiWafLCv/7Frf93C+9v2gjA5ClT+Pev/z9mzz0KRVEyfrL0sgDJX0JS23N/v+yUPbntFbu1ZG/KvhYL5YPCzsrFekcKtaWUDpVSzpRC+64t95dQes/o7e1l/PjxqYnT9zaHhXCqra1F07Qc61Jra2vB4GMejyfv7Nqqx4/q8aOoCgow/ezPMP3sz6CpCvHwAJrbg8vlAmD0kcey+ZWnefTaKykfM4FIdwezzryIstF2zJHkQ8Sjq0T7+zjxkq8TrKpl7Kxj+L+vnENveyt1jdPobmmicd5idF8AgLgRp2bGfPSAGzUepuuFfxCccSy+8TPRy6pTJ2z2n1tM6Ogl5ClUJomxn4WU5g0w7uQLGXvi+fRueYeu91bSu/VdjFAf/VvfpX/ruzQ/B+7yagL1jXjKq/FU1OKpqMFbOQrdF0TzeNE1beidJfdZovgZKj3ftnx5h8pTUplCAiptNSmMcoRMlnBKJyfPEAIqfx473VVAHOVrk0vLFE5JcZR+eKl6Uo8ZWzi1tbZSXVNDeWUVqqqyauWbNIwbz5x581EUmDx9Ojuam9B1nXHjJzBu/IRUvb/6+c/42113cMkXvkgoFOJPf7yF5qYm7nvoYTa+9x6maRIMBtE0LfHAFjRV4Z21a/nB97/HC88/z7nnnsctf/wjuq6jKgotLS289tqrVFZWUl9Xx+jRo6mssIW+kk/gpMRUcQGVTt56hokUElX50ocSTHnKJIVTIcGUfmvp7e1l9aq3ePWVV3jqySdo2bkTgMrKSj73hS/xyc98Dl3XU6LIzNCdmb9bSjgVEWqpvNk/eQm3u6FESL7bQ6Ei+YRUofqLCZtCgqaoqBvqOIpvzqC8PDCM3HvGvnKrOSyEk9vtZsGCBSxdupSPfOQjqfSlS5dy/vnnD6su0xLUPALBsARcPiwhtV1TFY699Fv0bt9Iz44tVIyqZ8ysBUQjYTateIr4QB9HnfsZ+vv7mXD0ErpaW/BW1GCZJoqqonl8xAyLUG8X1eOnYJgWiqLQ07qThulzGPWBD9H32mPEO5rpW/sifWtfRNHd6OU16OU1qIEqVH85qjeI5guiePwoujvvyTQc0TPSlqZcVComz6Fi8hxELEK7ttOz+W1633+bUFsTsd5OYr2d+YsqCprHj+71o/sC6N4Aui+IyxdAc3tRXR50t4fNa99kzfNP0rlzO2dfdjVTF5yA6vKguty8fP9faHrnDUSEM75yNWNnHgXA+hXPsOnV5fZdR4Sjz/kY42YelRIN4VCId154klg4hKKqzFpyKuW1haMIx2NRwj2dBCtrcCVEvWnEAQVNz72Uk/vJtjwmb8Bmvtf4nDugkrE9n4AarCf59EnuKKuOtDxa1uMh9aZvJdVQ7jmWbHcyb1IUWUpunYPVpZ6IqKpC664WqmtGpfwemrdvY+68+Rn7KSsvp62tFZ9/IhvWr+PxRx6kZccOuro6+ejHP4GI8OD999He1s7X/+tq6kaPob5+jL0b00QRUBPHaQkEyyv49W9/yx9vuYX169ejaRqWZaFqGi+88DxXf+tbTJ8+nVgsRm1tLd/4xjdYtHAhVsIiniF8UqIjS0Al0/OIpIKiZ3coZFUqkqcUwWSaFqFQiN6+Pnp6euju7s74tO7axZYtW9iyeTOtrbsy6igvL+fCj32CT33uUioqKrCwf51iwmYowZTvFjeUULLS5FbK0pkjtuwERbEf6un7URUyoqAnz9HkL5xso2maaWUUtLQXP0Wxt8ct+1mhJ+4L6b98tnhUUtdV4WMb6gyy9o0+OWA5LIQTwNe//nU++9nPcswxx7B48WJuueUWtm3bxuWXXz6sej7o2oilldMvLgbERci0lxHFXhpoKWGRfHD5x07DP3Yabl0lFI2jay5Gzz2BaG8n4biJpnmYecYneOq33ydYMxpPsJyaSUdQPmk2LS07qZs2h/fXv8PoaXMY6Gpj15YNTDvtIkJaAHXhR3E1r8fcsQ6rfTtixIh37iTeuTP/AagauLyIywu62/5oLkR1gabb3xNLVB1LTUvLzqMmP1rBV5JSur3sfKXf3Iess6KBwLwGAvPORDGihHdtJda5k3hfJ/HeTuK9HcR7OxAjBiKYkQHMyADR7raCVfb39zOjcSyrezrY9vyDhN+1BxU0NTXR39/PmDFjMAyT+390ObPnzcdfVkF/fz8SjmCJYFoWD/33N5i76GTKqmpQNY01r7xALBLB7fURDYeId7VSO7YRzaWjaC5UTUfVdBRVQ/e4WfPCUl785185+0vf4MjjT0fVXTz2+/9m+7pV+MsqsCyL6QtP5rTPfhVvYNBEnfy9Ci0B4sm8VnGLk5Z2h8zOY2VZkVICJ+3csFL7tDK2WSlLU1JwpO+HjDLJbSlrVaIBZtqDSyTTCqUk5u1t3bWLuUcvIByN4Xa7ad3Vwqj6MZgiKAId7R14fX50txcRoaKyislTprFp40Y8Hi8fOP1MLGDDxg309vbwxz/8jjdff51FS5Zw1Tf/i1G1tQjpOlRSDqqmJdTXj8FKtMYUob2jg4s/+Sl+8tOfZnRBCqCo9gMxQ/gkhVEiLceaJCVcb0UUQDwep6l5B6ZlYhgGccMgHo9jxA0Mw/7E43Fi8cQyFkusx4jGjMR6jFgsRjQWtwfZxGLEolE7LRolGo0RjUaIRCKEQiH6+/vta2UYrgVjGsay4NhjWbh4CSecciperxfTgriVpxsuo6uuuEBKCak8Nqf0eq797n/x7JOP07KjiQeWvcyU6TMA6Ghv4+p/+wrN27fi8Xj54X//ivnHLQLgxut+zLKnnkDTVBRF4cv//p+cc/6FqToffeCf3Pyr61MuJP931z8YM3Zczq310k9cSEd7O6qqEggE+P61/82Rc+2XtbadLXzj3y5jR9N23B4PP7ruBo5bcnz+rtWsAyshy5AcLgOnDhvh9IlPfIKOjg5+/OMfs3PnTmbPns1jjz1GY2PjsOpxKYJPi1FJLO/2iKXRJy56xE2/5aVP3IQsD2H0lFVHNwXcAbx1QcIxA4DyqfM45Vs307NtPWY8Rt0Rx9ITioGngvGLzualO37NhONOp6vpfVzlNbjHzqSjL4KiKBjlk6B8ElgmargHBjrRQp2okT7U6ABqzP4olgGWCdEBlOhARrsLne6ldmKJoiGqlrM0E99R1MF0RYVkuqqBqiKKmlpX1OR3dbCsomK5fMTLxwyr+wuAika0qoloiePxJttsxpFoGDUexoqFsSIDSGJpRfrBjCHxGGLGcceiYMZR31mHq6waLVCJGDGam3cwbdpUKioqANi6dSs97a24sPAAHu/gL7h2ZzPtm97GrKqitbWVzl27mDNnjr2xKkBk21qatq3N/F0TN9Hu7m7a2trwuN00v/wkofWvoigKu9atpb66mvr6ertAywaW3/B1UBQUTUdVNRRNR9G01HdbkGm2KEssNV1H1Vyouo6q6+i6/V3TXWi6C93tRnO5cSWWmsuNx+O101wudJcbl8dOd7s8aC4XbrfHtpym/SV6og8tXqB7z5XPl0pNvhFnibqEnrASU25mlsl8eVFRsUwIh8NMmj4TVXdhWILPHyQUCmFa9oPj7dVvMWfefIIVFcQtobJ2FOec/xHOOf8jfPXzn+XpJ5/k81+5nN6eHt5es4Zf/Pq3fPsHP+bqb17FnX+9ncuuuBK3e9CqK4BhmOi6QnNzMxMnTsK0JGFZUNm5s4VVq1bx8COPMHfuXMaMGYM70dWfJMMvSs28IiVxhRb1bRpGV13Lzl2c8+ELSs6/t9F0nfKycioqK6moqKA8sayuqWX8hEYmTp7MhImTKCsrT72gWkDMlJKsR9nCaLBMZjvy/ZzpYur0cy/gkq/+B5d+5BwsJGXV/eVPfsD8hYv5/d338/Zbb/LNr1zC46+sRlEULr3yKv79Wz8AoL11Fx86YQGLTz6dYFk5a1ev5Lc3/Jzb//k4ldXVRMJhFEUhbklOF9pv/3wn/oDd3fXsE4/xX1/7Kg88+y8ArvvxD1h4/Elc/h/f4J23V3PF5y5m2etvo+VznEyQqr/IKTSUHEqdo4dHdKPDRzgBfPWrX+WrX/3qHtVx785RVAU9VLhMynWDct2kTDMI6iY+1cKrmngxGUUE6E2VMwT6LZ1+00WvqdNv6fQaeiJNIyoqoCFVswBoao8BMUzDwqxfgG9qmE1vvIq/fgqjz/4E723oSD1UrURfiZW4+sWqACoQS0ADyyvgBcWMo0sMl0TRzBgaBpoVRxUDDdNeioGKObie+O7TQcMWXooYtghLuyErYqKkmZD3Bd1qDRu8x+fdVsgpWykisjK3+RKfxE3Ak5kv+VvHXMvoHncG5uQFAMSee57o4qvoCFahWHFcbTfRXduIfsTxaIrQs2U13e+vJNbXQd2s41Fmn0rI42XHpr9QOXUh2+I6VjxKsGo0ZTV19u9rmSAmWAaqmJixCO+vfocjFhxPz6vLsdxBxFuBWIbdJZHy8Uk7HhHEiGOmbEkjg6LpqLobRXeh6m40txtVd6O73KguN7rbY4svrxfN5cHtsbtGXR6vvc3tweuzt3m9PnSXG6/Ph+724HG7UVR1UGyl/Z/ZAuzxu27lqX/cwbaN62lqbubee+7mP6+9kXM/+2WevP9vVI6fis/vZ9Vbb3HRl/+D9qhCuKON8orKlI9V70CItp4+OsIGons4+4KPMWHmUcR0P8efdjZPP/4IG7bvZMzYcYhYtjM4CqZpoesmO3e1MmfBcXRH7etEiZtMmDaTLdubuPXWP9PX18ulX76Ccz9suw+kO5jn/K776MU+hE4wWIama+i6jq7bolrXXbhcOprmQtd13G43LrcLV0I4u1y2uHa73bhcbtyexNLtQU+kudz2x+322B+vB6/Phz9Yhj8QJBAIors9GQ7dSbGS7h7QZ0FfX3yw6yqxjCfvg1np6aIoJbZSeTKPP59zeD6qps8jBMRNix29UbTuCABPPfIgf3r6dbZ2RyibOIuy6lE88cxyjjpuCeCmpzeCZVns2LELVdfZ2RchqHj4w02/4cIvfJWwy0+oN4yqqvZ5NzD4gj4ooN309EQQy2J7Sysur4+OUBwR4alHH+bBFT+jIxxn9JQjqKmr55nnnufYxScUPJZCTuhFtNaQZfNxKBmjDivhtDdo7YOOmJZ488v0F3IpFhUukxqvSY3boMZtUOUyKHeZ6ApUagaVmpG33pilEDJUBgyVkKkSMhQipkLYUgh7FernLmJg1hIGokIoBOHYAMn3ANOw6xTLzFyaZt70QXTEUgB3zrbs9Y9N7aXGlyuMDGvwY1oKhgimZTtkmgKWZd+MLEuwBCyxEt8FSX7S1i2xEEsQy/ZUEEtQFEFF6DGF7vaNQ/1Fecl+Uy8FNa1MUjgZ0Qi9ra1EZCOSaGvXjhY0bwiAcDhOvDuE0RpBUTVM6jCrZmHJVvpCIZRwJarhIxQxiTQ1oVc0gGnQ0rSS8rkX4Koan2Pm61l7D54jP0pz3TFE1dfYVXUiPXVHABAJdrF513a2tA6gB2upnP9RPJUNKGLZfhRi2d1fYtlpYqIq9lJRBEVMEAtVbLGmJbcl1tXEd0VMFMtIiOY4WLZ4Vqw4ihj2ummAFc/oPhLTwDQNiLJPUDQXqsuFqrvQdHdiOZimu9y2FU3VmH/yGRy15BTCoX4G+vt47MH78Hp9xEzhhu99k3g8xmkf/TRbWtpo7nqZh2+7mUgoRM3oetxeH7va2qmfNZ/1OzuYvugUHr/7Nrb2RPEZOq+9+RaW7qFH3Ji90ZTYSYoul2bS3tmJGqhkV38MRVHQFIVFZ57P8Wfbfpcb163h61/8DHWNU5g644i8Dxo1691/dx5GxbSBVjaKvz31r7xdVell8wmMpBO2aVlYpkkkHscyDWJxEzPRzWeaBrF4HNM0iMTj9PTHMbtbMYxmDMPuEjQtk3g8jmkY9se0uw0t08SyDIy4gWUaTDnmRGrGNuYKJStTFKW3NdvfLymkdsdvU1MVDEto7otgdocY6OnCFKEDD909YQCCo8awbtMWqmccDcCye2/j2Xv/SndbC1/8wQ30KV76+qK8v2EDFfXjeeST5xMJDTDvhFO56IpvZIzETB+E8Zvv/gdvv2a7C/zg93fRHorT22W/TId1H9GQ/cJUVT+WTVu2Memo3B6SwoJp6JOqFA8MTVH4v2u/g2UaoCh4vD4u/X/XDF3wIMARTsMk3B9FjWmpt1slbRkB+oCmRKeQqtidQgpCucuiwmNR6TYpd5lUeix76bYIuAS3KrjdJpXu0qw2piWEDQjHIRSTxNIiFIdwzEysm4TjQjguhGJCv2EQjont7VFATImZX0CFYhZeDVwq6NrghaOr9ieRu0BrFQaNvcN3VDUswTDtN8qjzD7ilhA3B9PTlzFD7HXLXjeT263BpWnZ4i5u2vXETbtM3CzgFGqag9a9eJSBjhbCMXvIuZhx+nZsRHEHAYh27cDQyok3b7QdezUNtHIYNYfomnvpXvccWvVkzHgM0fxo1bYAMvufp2v1w7gmnoSiaoP762ki1rQayz+GyKbXMCP99De/SyRqoqgaypgF+BrtMrGdq9n1zK8JzP8cim53k0GuaEyup6crmg7oqTRNd2fm1TTQQE2kp29TXYN5q6QdnxLFxO5eNUXBRMcQBcRCU2zBZVsv7WVSfClmDFUMMGO2KDPt7lHFiqOYcTDjtjAzB61oYsYxzTgmDMu25gLWL3sYsI2LUyc0ALDjrRfZ8daL9n8yMIARDrO9rZl4PM6oUaN48o+/4slEHdvff5/PnzIXTdNxe9xMm3EEv7vmG4RDYaKxGLWj69E0DdM0cXs87Ni2heVPPMyq11agKCooCqqqoqgqiqKiqgrtra3c+vubqKkbbZ8/kFJHye+JsQYkbY2DWkBS6aS/mEjmOmIlvifLJF5W0stZdrplWUjykyxnJdNNLMvEMpPfbcGU+4K2b+jyjqYyVpYSQ9kiaHA5KOSzBVL2esqdopSQJKpCNG6yels3W80O4gO9WAKvb+5KiY/OgRhbOkL4t/cAULvkQj6+5EK6tm/knl99C3XCXLzBCvrDEd5evZqLrrkFRYS/XXMFsVv/yHEf/lRe38QPfv1nfBBY9czD3HTt9/jKDX9moCeMBWzuCqXy9kVitIdjbOkOp9pdUDBlJRcNrVKiYt+xdRNW4pni9vnZ1hMpqdyBjiOchuCmm27ipptuSo1kOGG8gu6BSNwiHLeIGArhuEXUVAjHhb6IhZllMlBVhTBgjwVRsbKGIGgYlLuFoFsI6LaQCrjA5xJ8uv3x6xBwg98leHQFTVUIuiHoBgLZoqTw32qJEIlDKG4lBJctqkIxez0UtwjHhIGomRBi9raH35eULBLLQFMSIkoFHROXJmhYuLS0dBV0xRr8roGuSMa6ptjfXZrt5OtKfR+8wHVVQVfB69r3tl7TkoSgEuKGEDOttO/Co++qzJ8SxFvmImZYrN42AX/0Parr5xGORNk0sIvqGcdjqBGi0Rgx0Qbf0qP9GOF+rIEe8FZhKSrxAfuGaqEj0V4Y6MkQNFZfK6J5CK9/3E6I9hPZ+AxaZAC1otHOm4xw7x+DFe4h2t2C4smNX5JPMKW2JUbmqFl5Bpe5IiwpopJlFVXjiOB26t2DMbXSiYtGDDcxcdlL3Biql7gaIK54iCv2uoELS/dgoieEQuKBkVDoioItuhJdxqpqW8NUzEwrmJhoWLY1TEzbBxAzbd1AFXs7Sf+/hFUN0yAQNAhY9nfSlwkmT57MhAkTiMfjKIqCW1OI9vUw0NdHd3c3vsRPtWHDBjo6OohGo/xr6eP4fD7mzp2LZVk0NTXh8/nweDz09/ejqQq9O7fSv2t73t/woEVR7XNH1ezzRdVQEgNLFE23fRk1DRJ+jpLh66ghqIlBKGrKR7Ip6qNpW5ftkgCYRtJvKVMMWWniKJlXsgRTqV10kBZHTFUwDGHbzj58Zk9qX2+/swVvWRUAbc3bCB4dYM2WzlQZAEWpwfJV8dyzyxl11ElQVod7+hLe2NaLpip4py9h9ZtvoMw5O0c4pYs6beoJbPz5t3nlnW24A2VYAive3pzaf9PWrYxVKnl9S9dgmaFCnBQRRUNZo7K3T/rwZYPO56rKyu3dfHb+3g9Iub9xhNMQXHnllVx55ZX09vZSUVHBJ46pJBAoHoeiL2ywoztGU1eM5q4YO3ri7OyO0RmSDJtM8oFjAu0xaE/fVqBrSVE1NEXwu7AFlUvwu1S8ie9ezRZaXl3w6SSWyTT7ove7we/WYBjhNEQgZMBATGEgnvjEJPEdexm1GIhDd9giYiSuFwGsPN2IRboG7e92N5KuCjpmQoSZCVEmuFTQSIoyQVcELbHUVdCSy4RISwq01Lpq34CSQk9NuzlpqoI30z+XpUuXsnbtWvp7enhz2X243W6uvPJKjhpzLg8++CAdq7ZgWRZnnnYKCxdOA+CRRx4hFo+hajqxuIE5eTzHnHAUluqmo2MRq155loC1EUvRaY/tpLxxLmq5n4H+fsKhfsRXi1YzHa1meqodsbf/hmvqOSjeCvvtP9oH7gCq5sLq2W6PhnT5iv+XRSwCqXfzRB61gLUKwDLiibRBUdWpqWD5cKkWLkVwqSYe1URVwKWYuAgTUMLkkHauAGDYzuBx9ITQchHHTRyXLbLSlobqwcCFoQZTYivbIpx0jlUKpGek6Zk+Y6ngn3bU25TYUhTLFmNioapC3DJRxMInFn4sTDFRFWHqQoupCEaoj1h/F7FQH4yfBvEI8dBjhAZ6iHeFcPvKOfpT38RTWZXRnzbovJvlKKykW3LtdVVRkqapwe0KtnDJLqOoiW12mkXSwqXaeVQlkUexHfPtMfQoioqJkrCaqaCqoGhY2PnNhMgRlIQQsi1qqS60tMBKg36ZyUPM7G6zjORIyjzixxToDA/WIZl1WVkiSREL1Ywlup/jtj9nwo9Ttwb9O1Wx7F9DLFQS3xG7qxtJdHFj+7EZEaqaX6UstBlFgbETp9H/xI2MXXAS3a3NWN0tTDZbUTY8TV9XG2VVtYgIAz0dDGxbS33/0QTXPsbk0bXsWHE/M3z2SN/ulx+kbtxEvCsftC2Mlm2BjseihOIRvP4gCGzf8h5erxflrceJAxMmTmHN/32PIxcsoXPXTsLtzZS1rCXUkhx0kvgt8vUOSCJdSWXL3JiDkvdroTyW5mIz9YUyHlQcNpP87ilJ4fS9G+/A7/fjc+v43Spet4bfreF1q/hcGnq2vTONqGHR0hOnqSvOto4YW7ribOmIEYrl/gXFhFOpeZMMRn1OCChtUGBlfkgJMF+W4BoulkAobvs2DsRhIGrRH4OesEVvFHrCJr0RoS8qxMx8wgmsocSWmVumsDCz8qYn0zTFtnLpmoKmJC1nim0F08CtK7ZFTbESgs0iGPDj0sCI9CNGFE0R6uvr8Ll1PLrCjubt9Pf3Yya6+qZOnYrfPxg1d9OmTWzatAnTNBk3blxqhN22bdvYtGkTx594MuGYxUDUIByziBrw9EN3MeXoU1F9VfSF47y/4n4MI24/tDQP+tjjUIP5Y0EVO08yuuTSKCacClmjcpaaikux8KgGXk1wqwZe1cKdEFUuxcCtWLgVA1fioym7d1uyUGwRpbgSSzeG4sJU3Qlx5cZU3BiaG1NxIZoHU3FhqImQHIpSXDiRaXFIMmgNS24jo+ygpWGwTI64yxMmInt9uG/8w6GYn89QXVzpaWYhq05SDA1HOMUNVMtAsWKoloEqdletZhkJERRHNeP20rLT7QEvCWFk2UtN4qiDrwV7zMsvv8y2bduIRCJ4PHbQ4wsvvJBIJMLzzz9Pf38/qqqyePHiVJDlZ555hr6+PttvSVWZM2dOxqju1157jebmZhRFob6+nuOOOw5FUdi+fTvbt29nyZIlDAwMsGzZMkzTRFEUvF4vxx13HJWVlQBEIhFeeOEF+vr60DSNxYsXU1dXt9eOe0+wNDdNcz/DM187cZ/vK/m87unpobx8eDNHlIIjnEok+UdUfOBbCf+R3AeKAgS8OrWVPsZV+2isK2NstZdxNT7qK9ypodjZ7OyxBdT2rjg7+4RdvSbtYdvhWk2USd6Ys9+c07el4sxl3fCz37Ih/80/3zrYsWg8qoVPtUcO+jQTr2rhUxJLbXBEoU+z8OQJXliMmKUwkHCMT35Chkp/XCVkKvRHFUKGQl8s8e6XMs2nOSFL5s168I01lkjPdJjPSCvBCpadZpWQV1UEl2Li0QSPhr3UwauaeDTwurBFrA5eHTwJ8erRS3O+zCYctwjFFQZi0B8jZ9kXtcXs7lCKc33ONaFpedPT0woJM01VbEFFHLdi4FbiuFULlxLHjWEvFRMXcVtsYaApe/ZgFMDEhaHomIoLS3VhYscyMxUdK/ERNfE9LbyGKBqWoqHomh0iQEuE3cgIv6Em0u3vqqpmdEVmi63B32RwvWDU55JHj+4eYloZAwwUyx4wIGZiIIGVNqggMTLU7jYdXLetc8l0e8QoVlIYGaltakLwqJaxV8VOEgsVU9Ht/1axxw1bip4YP6xhoSJoWIqCPSxFsdMUNeFVpqTsNvYiMW1Maj3N6kciJpeSK5yFhJ9kskyyzqz/WMmqb9CIM2g3UtLKC4Png2QWIGVJzNwDtkExNz2boeI02edaolUZjwEBRSXkG83L3zu9aB17g30tnJyuumFiRAZQdKOgw213BHr7Q7zfNOgDorrcqAqMqQ7SUO2lsa6MSbVeJtf5GFWmM6bCxZgKF4vT6jMtoXXAYke/sKPPoiUEOwcsQom/THcNijDdZe/b7U4sE0IqufSk1gfbnNyWesvNiqeTrz89m1QEbEvI6HwRC10MNMuw3w7FQLfiuCSObsZwWXFcVgyXxNGwUo7xVUM4xgsQQSckLsKiE5LB772mRlh0+uOKfaPLervN5/+QbcpP70JIT09/c87uFkg75Izt2XWA7bwcB/qyhkTbKyCRlBeZLbgUMyEcTDyKgUc1UkuvYuBVDbxqHJ9qW2h8LhWfC2qKTAUVt1T6TDd9ppueuJs+w0Ov6abX8BCxtJRlbrD9hf+TfeEEnNEVCEQSn+xt+cqoYuIihk4ctxJHF1tU6RJHJ7GOHZIjua4TT1gj7MdeMh0Jsw+e2TlYie6upPiStPXBrjY1EcDQfqimvqOkPWPTH6r2w0tBEg8vIfVYTziKKwlHcPvEFTsOlAyKI3ub3fWIZZG3a2c/Yikqoui2mFVtUWsl15MiV9FT20xVx8SFpeoY6FhpeUjrNix0HedtQ9a1XUysJik4U00eAVJqSJViQriU8ACFyhetdy+Ibyg9LuCBjiOc9hOWQEtPjJaeGKua7dd+zeMj6FGZUh+ksVpnbIVOQ4VOXVDF71YZU6YxpgwWjBmsJ2YKHRHojEJHBNojQltMaDuQBisoqt1ForoxsqxPueZ8A7fE0UxbWLnFQLdiuDBwWYZtXZA4bnucFj4MfIqR+zKU8EmyvBARjZClExKNAVMjZGn0xTV6TBc98WTMrL1zI9g3KMRFJy46Yrnz5si8iQsuxcKnxlMff/ZSswWWS7WoViNUuyKDkUATxC2VPsNFr+G2440ZbjqjLrriHgzZi9N27CMsRSOKjyg+wnmsrFDAGiuCiomu2F06Lkw0iaOT8IFRTLu7BzPlF6OKbQ1RE+JCxRz0iUlaZdJ8ZQaXmagJsXKwIYrtv5QMapsMViuqlkq3Eo7dVtIqp+opy5yoScudllpaqo6FvW4kLHtJ8SNpo3GzX3gk34tInu2HNWLbzuxzFnuZWE/a1VRLGLSxZX7UhP+TkrS5yaAYH7S/5fefSp7zFgrb9bH77ZD3JY5wGiZGLIJiWgV9QgC0Am/iyTJ68g1ZdxMy4b1O4b1O2x7h9tkKoK7KxZigwoQqnXq/QkNAocYDbk1hTADGpBy77dPSEgiL0CcuooqLAVxYLi+GouFNWJpcaf5X3kSXX7L70JVlaSoWwTn75aOUoanZ8VOS97J8vhLZwenipj0M2oxHkWgYIxZBYhHMaBiJRyEeReIRiMdQFcGvmPjVxH/gyqkeS9UwdR+G7sV0+Yhp9vew4rG7ZYrEdkmmGQXypMrK0HXklrVK3m92enaaBfRakhaCFSzTxGOGcRsD+MwQXiOEzwzjs0J4rIgtqtxRqt2ZQZcECCs++tUAfUqAfjVAvxIghNd2HC5gZcv3wCr0cMtmdx52pQc7LdxVnZ5mqQoWYGZ3neTpUsv2YcrrvyRiBxVMOJFrCTHlUmVw3bJQFdvKo9keyPb1JsnvaWlJstdTx5fWxZNmsVJVhaSTd9L5O+nkrSQEEYoCqmZPWJyI7o+ipkbEWaVM7ZIg5zxPj6tU6PrJCiWwJ9fE8K7Bwufd3j4nIdcnTRELt2LavloJq72WCEzswkBLdHfa4t1Ek8TghGR6wrk9Jehl0Ml98CCHfRh7BVPVkfoZI7PzvYwjnA5QemPQ2ylsjQAI/oALFWFMmUalSxjjVynXLSp0i2qX7VcUUEwCmKQ6N0yIoxK1PIRVD1H8RFV3Th/6wYCiKKguD7g8qb5zJa1rLW4lAmjGo2BEMaMRMGJIPIISj6LEI6jxMKoRs28ysX5csf6c/ZiKTlzzENO9RDUfET1ARA/Yc/Yd5IiiEtEDDKh+ush8EChi4bXCeIwQPiuM1wzjt0IEzH48EsMvYfxmmLq0sZ8mKgOKnz7FT7dSTqdaaYsph/wkRqWh2VMAi5qYQU9XU34pwuDDVpIP1bR5HJPbsrtO9uaQ8pLyOq6xpSOC24risSL2x7SXrqT7gsRxWfZSk/2jagQl5bMlia7glA9XKj1hT0r6cKW2p/lyke6jleb7lXHa2CvmbsTwO1A5+J8G+xmxTFAGg7wpVq4DrKiFInXbJAMGujz2MmllAvAH7TR/uT3nR12l/SAaVWYvaxLbNb+LAWAACHl1NDEpVw10M4pX4mhGFNWM4cLCZYUJWmEwulE0Hc1fjjdYie4L4tKTzrj2/lO+TokTP/3GmTzts6eCSF5H+S6LoZwJi41NyJ3BPLNMegdHMk/2NAvZ6YZhEgn125+BASID9vdoqJ94NGK/3RkGXiNzLj+3L4C3rBJfVSW+skpcwQo0ffDySVnJsnpd8kZY3gPrW6F608sOVSY7snJmmcw8cUuwomGifZ2YfV3E+7ow+7oxB7rRLJNy6adc+hlLq/0m6/ZBZR1KeR1WWR14gxn7yX7Tz7bcZVvsCh1jvvVStg1nnsOhxEi+0W56uoWpyLLYNreeeSVllElde9ltyzyWfD4phcTUcKbN2BOKxUrKuSaK/bcFNhWqP9+5lNo2DOtUMUQEiYWRgR6sUC9WuBcJ92OF+pBIvz2FwnDQdJTEVEWK7kpEyE+sa3rWx55XUtF0FNWefzI9ZpaiaqnJwpPWw/T7cpGB4HudOftvV/sURzgNQXYAzKRwGs4oo2R3XjIas+a24+x4EoIpKZZgUDCNqbbzjKuyPX3rEwKqMlGmKi3QUEUiXoA/4SSeXLoUi3g0gsTCRPt7iQz0IaaB0ddJf18niqpSVl5BRXUtVZUV9jQQWc7ialp/tZIM/pd8K0qKx9RY4qxl9nfIuYHkzO6eB8n2sEy+gaen50tLL5u+rCgHygdHwySWhmnR399PX18ffX19dHV309nZycDAALGw/eltbU7VXV5eTnV1DXV1dYyur8fjHbS2FLvVZt+HcwL25SmcnZQtIjO2peotXiaf8Cw0QaqIHbRuUFhZ9PX00NPVTndbG207m+jYtRMrFoa2rUjbVhTA4w9QVTeG2rGNjG6cgua2z++kwCw0r1gxITgoUrN/t7QuoCGedfkesqUKiHwPmhwhk9XVnT5dRm6eZL35uxPzDcbNFU75y+Ztf4EX//0loIqxO4EoR2L/8ViU1uZttO9ooreznd7OduKxwvMKKYqCNxDEHyzHFyzDHyzD7fXh9fpwe324vV7cHh9ujxfd7bYjyhc4vOxpdzL3U6DMHv5U2fUe7gZHJxxBiSSHN7qO+hyK5i4YswZAS1iUdK/tiKR5bBHk8trTcnir7CBggXJvYjk4o+zeFE7uxN06uVQRIgN9hHq7CfX2YBqDY9Pdbg81o0Yxqq7OntizqHBKPvUOIOGULZCyyw7Gasgpk5MnVdbOG41G6erqoquzk47OTro6OwmFQjntrK6poaFhLA0NDQTLywta2w5O4WQv09/es7dF4zE6d+2kdcd2OnY209XakhG2QVFVahomMGbSdKrHTUJ3uR3hhCOc0jmQhVN4oJ+dWzaxY8tGOnY2I1n3LkVRCJRXEqysIlhRRaC8gkB5JYGyCnzBMjv8xDCO4VAVTjNH7/3wANk44Qgc9hqqquIvq8BfVoE+DiLhAQa6u+jubCcWi7KzuYmWHc1U19Qypr4en7/IuPbDCI/HQ319PfX19akbRiQSobOzg/b2dlp27qS7u5vOjg46Ozp4e81qAsEgDQ1jGTt+HNXVNUN2WR4K6LqLurETqG0YD0AsZtDV1kLbju20bN1If1cH7U1baG/agqpp1IydSN3EadSOm5Thx+PgcKBgGgZNm9azfeM62ndkToMTrKiibnwjVaPqKa+qIVhZhaY5j9TDAcfiVCJJBes++vMoWu5Ep+kWJz1hYUpanHSfbWlKWpz81aMAKE9YlYKVg108tVV2WmONXXZcIk9DwjpV4bEvzIq0cN4VHtv65NEzLUypZXKkXNqrctIXQzWimJZFR0cHLW0dhMKDcQ1GV5UzoaYs5QeVtB6lLE+JZcpqZGVaojK25bNGZeVNkh1LqBBKvodtdhdqtrUqT2CVHItWoeArGXkyhdBAJEZzRy/b23to6erLeHOtDPqYMXYUk+qrB3/LfPsttv8h2ij5rG/Zbc22sKWTbcXLWqbKpP++Wdsksa2Q9aqjvZ1333mH9evX0d3VlarG7w9wxOw5zJpzFGWJt8P81rDkeqbFLN8dLHvfudsLbChQ31CU+qavZGzLGlGV7adUdH+F/JWG0bYi9Q93v8Nhd6wfxYrk/MZF2lho39nJpmmyetUqXl6xgoGBwUEkDQ0NTJs+nalTp1FVXV1kP6UdZLFsSikn4hAW+1Is+vua7PuN1z+Mub52Eydy+AHCoSycALDsaUH6+3rY2dZFZ699s9BUhbE15YypKkt1NzjCiaJ3vLhhsqPTFlHb2rpTXUouTWXymBpmjKujIuA97IQT2N1rIsKuXbvYsH49695dy0D/4INp4uSpzJl3NA0TGlMPQEc45dufI5wK1bcnwklEWL/uXV544QV6ursBKCsv56ij5nHErFlUVFSU1F5HOA3iCKfDmOQf4VnwpUwfp6xZ5WEwYrgrYP9hScHkCtgXXVl1ZWJpi6KK6sFJWRtr7ZOqMRH+eULCx6kuYNdZnfBxCrgH9xdIRBHPFkyexFJPPGIUc9CnSTEiiaU9JQmJbYppr/cOhHm/I8RAzBY2LlVhUoVOnV9LCSdJ+kglpyyJZ66n50mNMEwtE3PHmXlGHmaLqT2JUF2CE/+w8g6zLTFReb/fYmOfRZ8xmD7GpzC3UqPaM3jnzBcTbLBthRxTisxdmFVmMD13e2rf2fUlz3PdlawkvUIgj49Yoq68vmOJsA7JbaYFG7Zu582169nctDOVrba6koXz5nLkjCnorqx9FxJ1eQWnkndbXvFYwMdtb1L0YVjCQ27IB2EpD8r99UAulKdY2SF8IvPue3Ciu7zpeduap46mXe08+dKb7Gy3raFBv5cTjz6SedMnZkwGXYhhiZRSXg6H82jeA4FUqN0FX+x2e0eJ2Somztu79ebB8XFyGBEqfC7mjS2nrTfM1u4oUVN4ryvOrpDJtHIFn37o++zsDdyawswKjRnlKrsiwnu9FjvCws6wsDNsMCFgC6gy1+H5e2qayszJjcycMpH2rm7eeHs9q9dtoL2zm0effZ7nXn6dY+bN5ujZR+Dz+Yau0MFhmPSHwjz7ylus3rAFALdLZ/HcmSycMwO3S98jUeJwaOIIJ4eCKIpCXcBFrV9nR1+cbT1ReqIWb7bBhDKVsZ7ipmaHQRRFod6nUO9T6Y8La7pNtgwI2waE7QMGU8tUZlerh7Ugra2q5KyTFnPywgWsfOc9Xl31Nv0DIZa/9CovvvYm8448guOOnltyd4mDQzEsy+L1dzbw3GtriCas5UdNn8Spx80l4HMCuToUxumqK5Gk6c97zGV2YLKs2EwZ4Qjc9kXnTnTNufwVifUyIK2LLtEdN3bUYJ9v0rdpQnJbwrcp2VVXnvBxSn/A5nbVJUyipu2/pCT9mJLdcmlpya65ZJddqjsvsV1ig87i4ajBprCLbtM+1qASZ7rehyexH4kn6s/TVZfqzsv2g0qmp3fZ5fGVsssO3U1Wqn9U0Tqyo1iWuG24++/GzTtqDbuUhC+cWBwp7UyS3hyfjrz+XOnb83QlZKcl68gbgyyr+y6n6y7ly6fmKZO1LbmeCMtBnjLZ3YmKnplX0TRMy+Lt7W28/F4zrT12QFJFUZg1oY5FMxsZO6o6eWC5x1OoGy/HkWiw7JC+biX5vu1G98aeWDSGU3ZPut8Kndf5Hh/Zw/RL6Zor4PuY93oqdO8o5b6RyLOto5cn391Ba599f6sv83LWzHrGVhYZSTyMazuvC8JQ7I1Js/fBxNt5GY4LRBaeUz+3FxuSH6erboTJDoB5OONVhVn+GK1RYXPMR7+4WBWvZLLSQ61yIM0yfHBQSYwl1k7a8LFWraFL8bJKqaNZgsy3WglgDF3JIYymqhzVOJq5E8fw/q4uXlq3nS2tXazduou1W3cxoa6KRTMnMn18/V6bvd3h0KY/GufZ9TtY02z7MXl1jVOmj2ZeQ8UBEcfK4eDAsTiVSFLBBhb/G4ruSTmA53uLT46icyctTWX2m3EyUnjS0lRZay8n1wVTZZNO4WMrbKvU2ERAzFEJi1PSuuRLm5YhaX3yJNKSI+VSVqV4JGPd/p4QOkamU3jSwiTRSMY6gCSsUhKPE0Fjg1QxoNoWsWqjl8ZwMy7MDMtWtsN48u3PiicczBMWnPQ3y1RalnXHMou8hWaVLbSesa1APVYJdeTsxxp+mYxtQFPlWDaNmoSlamiWyZS29xnbvaPoiKJi5LNCFUrPtmgl8ySXanI9j+UrO28yj6ppGekAqksvWibHgR0yrFC7Qiavtcd4tyuWGhlX7laZP8rHUaP9g9fFENawvE71Wn6n+hwrWT4KOeLvJ0qxxqYocN4XrSNnwEa+87vAoI4cC1AJVqQC1qOMPAUsTZLHWhW3hDe74ZVOIZY4b+aUwfEVFv60v2pvWKyLMZTFeiQpdL8YVh0lxGMr//yP93g/Q+FYnBwOSLyYzIo3sUOrolmrplMvpy/gY3KkiYo04eRQGgowvruZmoEO1tXPoNtfyXujp9EWrGXWznfxpI2IPJwZ7dc4d2KQkxss3myPsqo9Qm/MYnnzAP/aMcCMai/z6nyMK/ceFkFHHYpjibC2F17sFPoTBtx6D5xaA2O8sJ/m1HU4xHCE015m3Kgyjp45nvVN3TT1jnRr9i0qMM7sotLoY5N7DBHVw3rfRMaoPsaGd2ZM1+JQGv54hKO3r6K5soFNoybTFajitYkLOHLHOqrC3SPdvAOGMrfKyeOCLGkI8G5nhDd2hWkNm6ztiLC2I0KNV2dunZ/ZdcGM0B0OhwciwoZ+4YUOoTPxHlemw4lVMDPoDGrZmwgQRSOETgQ3ETRiqMTQiKOmvqsIJ7FzyPoOBhzhtJc4ZkY93/z4ccwYX5NKW9fcy3Nv72LVToP1OwdGsHX7lqBEmR3dylatljZ3NTu99fTqZUwZ2IIHx1IyXBRgXPcOqkLdvN0wiwFPgJXj5zKpfQsTO7ftdtfdoYhLVZhb62NOjZeWsMVbbSHe7YjSETFYtq2X57b3MqXKx1Gjg0yu8hWcp83h0EBE2Dxg8a92i5aEl4FXhUXVCvMqQHM8U3YLCxjARS9uBhQXIVyEFJ0BXITQsUoYFKEdQmEdHOE0TFTdjaK7URMBAVXdzdFT6/jtv52K161hmBYbd/YxdUwZM8eWM3Os3b86EDVZuzPEY293silsn0CVfruO5MS9YMdPAgi6Myfs9Sbu+Mko4O60KODJUXRD+jaZaV1o2b5NEVvYJX2SUr5O6T5O0bC9TAa1TI6YM2KowMR4F+WeKraUTWJAD/B22UzGGRuo6d9pRzok17cpn9+SFTNy0uy8mb5Pef2IhvA1Sl+3hsxb2Kcquw05flFFygyuZ01Sm+d4pr+/nW0zj6Zj7CQ2j5pEh+Jj0uqX0eN73h2qFHGoVrNH5OWZ2TbbPym7bLYfU768SpYflKrmKVNodGBW/QFN5XjgWI/CRtPNu4aHVsvFhs4wGzrD+BSLGR6DI9xxanQrpz1qVgDSfO3Pdwx58xyACm3IEaElXE/ZefOVscz8Pk3F/P725FoX00KAnZaLN/RKWhN+l7pYzDJ6OTLahbvXoiunht0c/VaAogFsSyq/++dMvgmEd3f/cVWjzxOk3x2wl54AA25/8YCYInjMGF4jituM4TYNdMvAZRm4TAOXFcdlGfREetn7Hkf7H0c47SFuXeWazy3C69Z4fs0OfnzXa4QsNxV+F2ccO4UFkyqYP7GcMp/OcRPLWDAhyNVPNtHSf2iOmKqOdhEId7G5agZ9nkq21cykx1fD+NZ3cVmO9Wm4aJbJpHdeJ9jZxrZZ8+mtrefdxWcwZeWL+Pu6R7p5ByRuRZilR5mlR+kSnXcND+8ZbsKi8lbEzVsRN6M0k5meONO9Jv4DT+M4lIgAzYqHlWo5u3RbMGkizDD7mBPvwoeFUFwwHs4IEHL76fZV0B2opNdbRtiVP9CsZpkEYiF88TC+eASvEcFvxvAaEbxG7LByzXCE0x7y6dOOYFxtGa3dYb53+8uEYyYur5ueUJyH3mzloTdb8fk9TB7l44cXTKKuzM3C8UEefLd7pJu+z3BbMaZ3rGFXYCzN5RPp8Y9iYHwF49vepbynZaSbd1BSu2MLgd5ONh19PFF/GesWnkrjO29Qs2PrSDftgKZaszheC7PYF2WbobMu7mZLXKfN1GgLafwrJIx1WUz3GEz1CT5HRB0UiMAWvKyijLbE3KGqCNOtfuYZPfgxsRzBlIMAA+4AXf4Kuv2VdPsqiOvunHzeeISyaD/B6ABlsX6C0RBeI5IbX+4AtKzuDxzhNExU3YWiu9A8PmrKvHzh7NkA3Pz4exiqD5d3MByBx2eHEnD73TQNmNy3posrloxm4YQgS7eFqPDbJ2yFf7CrLui2/5JkF11ymeya05PLtG4WNTE0ZDB4Zeb8c6SCXKbNVZcdfqBAF12yew5yu+hS2xLrVnKZ6GobFdlMoGcXW0fPJuIpY3P9UVR4amjYsRYlmgiamcdUb8byhyoo2j2Qk2f43QVWgW6D9DLJ7rXBvMl1yVt/endcoTxW1oyzmWXS8u4aYNyWB2k56WQGxo9ny5yFdOll1L36KsT3z/CgfF12hbZld/flzZM4j7NvwGrersGsSXFToRLy5FVz2xIEjgFmqzrby6vZGqim0xOgKa7RFNdY1ieMjvYxIdzDuEgvnsR1ldttWfhhMdSDpJTh2nuDUobV706g133VJZ19HRWqy1RUtpfXsLGugX6vHbpFtUwa23cxpaUZXzxGv2nRT3Gyr7lsisUFK/Qf5zsPB8vkP3eHU0+xrvXBsrn1mppOd81oumvG0FtdR9yTaVFSTYNgbydlPe2U9XYS6OtCN+IZbbSAUJH97q/z+kDBEU5DUCwA5hfPPhK/R2fNlk6eemvHkHW9tq2fLy+qY3y5i7rA4fHT+2L9TN/+CjsrJ9FWM5meqnEMBGtp2PIGwb62kW7eQYcWi9Hw9FI65h1N59FH03PELKI1NdQ/swzXwKE7AGFv4rUMpve1Mb2vjQG3l23+Srb5quh2+2jxltPiLec1EUZH+xkf6WF8vA/f/orI7JCXiO7i/Zp6NlePJp7wL9VNg0ntLUxq2YHHcNwA0om5PXTXNtA1qoHeqjok3ZfPNAh2d1De00Z5TzuBvi5Ux2l+WBweT+894Morr+TKK69MBdRK4nFpfPDYiQD87tG1JdU1ELN4d1eYOWP8zKxxs7WYhD+EUBHGtL5HeV8r28fMIeYNsm3a8VS3bmLU9rdRnYfSsFCA2rdW4m1vp+Wkk4jUjWbbRy5g9HPPE9y+faSbd1ARNGPM6mtlVl8rvbqH7b5Ktvsr6Xb5aPGW0eIt4zURRsVDjI32MS7aSwWGM7JxPyBAu7+MzTX17KioTjkn+2IRJre10NjRgsuycqxUhytxl4euhkY668bSX1GTEXPBO9BHZftOKjtaCPZ0oIpV1HrsUBxHOO0mJ85uIOhzsaNjgDc2taEVcKjLZn2rLZymVHnYGjq8AkUGwt1M3fg8LWNm0Vkzkc66KfSVj2bMljfx97SOdPMOOoJN25nw0EPsPOUUoqNGsfOsM6lcvZra199A2ccRkA9Fyo0oR/btYvZAG32am22+Cpq8FXS6/bS5A7S5A7xVVk+ZGWVsrJ+GWD918RD6YeQUuz+IqRrbK2vZXFlHn3dw7rjqgV6mtO9kTFeHI1wTCNBTW0/72Mn01I7JsCwFejupattBVdsOfKG+kWvkIYgjnIbJtge+Q3l5Obc/sgKAp1btQvcE0dL6jd1+e9JWdyK0gDfhw+Qrc7MlMZpuWrWHVd12/gpfuo+T7dNU5rGXHj0z/EC+cARKPMuXyTIyl2bmOqT5KSWGtaemU0ku8w3FT1qGkj5USb+h1BQIuT5CSX+n1BDluEH91rfwt22nZdIC4t4g22aeRPnOjdRuXoVmxnOnY0n6IA3hx5S+76HCBeTPW9x/KX/eTL+lpO9Edl3paYV8Pgr5eeTbb6pse5jav/6DnpOPp3/BXLrnzmWgto7qh55E7+ktWF++uvY2pbzRFvbnGNpnYjg+JYPppfuW1CY+Ea+H9rpRtI+uo6u6mj7Nwzqfh3W+GhTLoqK7m+qOTqo7OqkY6D9guz3yXQPpSBG/n8LnUDE/qfxl8tVlahpd9aNoH1dPx9h6JDG0XzUMRm1tpn7jFoLd9vnctg/P2/z+csX974qF6sjJm8f3Lrds8fpVTSXu89M9ZSrdU6dhBAan7PJ2tFG5ZTNlTdtwhe0uDRMI7WOfreHUeyjgCKfdQESIGSaRmMHSVcOLhLqpM4olQo1Pw6dD+NCMSjAkwZ5dTFrzFK3j5tBTP4XeMVMZqG5g9IbX8LU63U3DQTEtKp99Ac/2JjrPPpV4Qz2tl15M5dLl+N55b6Sbd9DjjUQZt62J8U3NGJpGZ00NHaNq6aypJuLz0V1dTXd1Ne9PAz0ep6ajk9qOdmrbO3DHHd+bQliqSlf9KFonjKWroQ4rLQ5SoKuH0e9vY9TWJvT4YXqTzEKAUEMD3TOOoH/suNT8iFo0QsX7m6jctBFvf8/INvIwwRFOu4GiKHznL6/icalYmn/oAmmEDaEtZDI6oFPrVdnef/h2qWimwej336CsfRu7Ji8g7i9nx+yTKduxkZoNb6Ca0aErcUjh27CZ0bv+Tue5ZxAbO4au884iMm0y5Y8vQw1Hhq7AYUh006SutZW61lYECPv9dNZU01lTTVd1NYbLxa760eyqHw1AsK+P6s5Oqju7qOruRt+LARcPRkRR6K6rpm18Ax3j6jHcg0Phvf0D1GzfSe32HQS6epzuuASiKPRNnETHnLnEKqtS6f5dLVRufI+ybVtRk1Z+x29pv+AIpz0gGrdw7Uaw2NaQweiATpXn8BZOSfy9bUxY+SQdjXPoHjuDvoaphGoaqH33FQJtjvVpOOi9fdTeeR99i4+hb8mxhGdOIzqugYpHn8a7yYn5tDdRAH8ohD8UYtz2JgToq6qkvbaW9toa+srL6S8ro7+sjG2NjSiWRXlfH1WdXVT2dFPR03tYWKRiHjdddbV01dXSWV+H4RkUS+5whFHbdjBqWzOB7t5U97+DLZh6p0yl86ijiJfb8bbVWIyK9zdSteE9PL2OdWmkcITTbuLy2v3KyZhNbt+gj5Mn4bPkD9o3CF+Zvawss+M69Rv2W0F9UKcpqlLhGfwbyhLfk/Gckj5NepZvk2rmxmRK+TKZmdOppHyT0nycBtOs/EsjqyzpPkFmxrZicZayp1QpNAWKYlrUbnyTQNt2WmcsJO4vZ9e8D+DftZWqd15BD/fn+Dal+1IVqjc7NtPuxGQqFosp26fJjJkZ6xk+TkPEqsn2m0pvd26ZrHZkbfc+9RLaqg30ffyDmHU1dH3ifNyvr8H/6HLUiG3JM/P44hRyH8mX91BF26MZYLsJsIUAYPg89I+tY2BsHf1j64hXBOmpqKAnbXSuq7cfX2sXvrZOe9nehRbL7JraV7/9cFyFhtOGuK4TGTuKyITRhMePJl5bmbFdDUfwv99M4L1teJvbUEToBrpLb85uk22Qyf6v8xlsknkKlVWVZCyytNh6BeI2JfMU86VSXBoDR86g57gFmBUJwRQOU75yNWWr3kYzDSwgnG8/BXyZdseXaqjyQ5XJqaOEOFQHC45wGgE6IvbDrsJ16JxIewtfTxvjXn2Urolz6J4wi9DoRsK1YynfuIryjauc0WLDQN/RSuX/3sHAmScQWTKf2DFziE+bSOCfS3Gvf3+km3fIo4ejVG7cTuVG22oaC/oZGDuKgTGjCNfXEK0qJ14eJF4epHfq+FQ5d3cf3o5uPN19eLr6cHV04+nqOyDPfculEx1VRXRUFbHaSiJ11cSryzOGwiOCu60b37adeLe14N3RjnIYCfFSEUUhNHsmvUuOxSwvA0AdCFH+5luUrVmLmvT1OoSdrg8WHOE0BMUCYO4uHYlJfstdQ2Q8TFEtk5r33yLQtJGOWYuIVNfTM2MBA2OnUL36RTytTSPdxIMGxTAIPrYcz9oN9H30TKzaavovvRD36vV4Hn4Gtc8Jmrm/cPeHcK/fStV6u8vUdLsIj6oiXFdlL0dVEy8PEKssI1ZZllnYtHD39OHp7LHFVG8/rr4QeiiCPhBG24cO1AKYXjdGwI9R5iNeWZ4SS/GqskyRlEDv7sPb1IqvaRfe7bvQIodX6JXhIEBk6iR6P3A8RrXtw6T1D1D22kqCa95Bk8PbL+5AxBFOQ1AoAOae0BOzhZOuKricl4eCuPu7qX/1CfpHT6Rz1nEYwUpal3wI/7b1VL79MmrkMIkguhdwbW2m4je3Ez79eCLHLyA2dwax6RPxPvkC7pffciwAI4AWixNsbiXYPBjDzPC6CddWEa2pIFoRJFpdTqS6AsvjJlZdQay6gnwReZR4PCGiIujhCGo0jhaLo8YN1LiBEjdQDQPVMFEMExQFU1URTUM0FUvX7O+6hul1Y/o8GH4fRtD+iF74UaH1hfC0deJp78bT2omrpRMt7AzsKIV4bTU9p51IdOIEANRQmPJX36Rs1dsoyZd1x+H7gMMRTiOAKRA2LHy6im83nMsPJxQgsPN9fG1NdE1fQH/jTEITZhCpb6R87asENr+D4gQgLAklbuB//Dncb73LwEfOwBw/hsj5pxNbMBvfg0+jbx9eaA2HvY8eiVHWtIuypl2pNEMEI+gjWlVBtKaCWEWQWHnQFjUBH5bbhbhcxCtcxCvKitS+Z6jhCHpfGFdvP562LtxtXXjautCzRJITyHtoTJ+XvhMWMjBvth1WwDApe/0tyl55A910wi8c6DjCaTdJBrxMOoV70oJYJp3Bk8uaCi8AdeXJpYewCT4dGoJudPdg2aSjeDLwpVe3TVKpSX6TLx/pczNlT+qbdAJPvLEkA19KWpmcwJepIJZZZuH06VD209Qo+YLqqUaM6jUvEtj+Hp1zjydeUUv30ScxMOkIKle9iLu1OaNsMafwUvef11l7iECXhRzO0+spGEQzjyN4thN4bl2JOtIsRtmOvMnqUulNLbj+9w5YeBTmOSdhjatn4MrPoL62Bh57DiWt+y77IVjMSbiUB2apD9XdcYje0wf23nTC3pO25G3HQBx29QKDo0xVwA2IS0MCPgj6kYAX8XvB40I8bnDpWC4dXDq4ddA1+yPY9wfDSiwTH9NECccgHEXpD0N/CKUvBH0hFNOyRRy2Y/LeZDgO+YUMMIXq2B2H7/T0ofKk1o2h8yi6TnjJ0YQ/sAjxJiaBX/sewSdfwNXbiwlYWY7kSZKO4PmCww4VlDNfMMrhBOe0sPLWW2h/+TiUOhwd4TRChAyo9oBHlUPqhNrXeLrbGL38fvonzqLniGOJV46i7eQL8DZtpGLVS2h9zhDdUlBE0Fa8hbrmPcwPnox17BysY+fAnOkoT7+E8sIbg10FDgc0StxE6e6H7v682x0L0IFBfMoEBi44A6s24cfUvIvgE8/h3pLw2XS65A4aHOE0QoQN+27mVvf+29uhjiJC2ea1+LZvoHfWsQxMmkVk3FSi9Y0E175GcP1KGIaF6XBG6Q+h//1xrJffwjz/dGTCGOTcDyCL5qE+sgzWbBjpJjo4HNRYAR/hc08ldvQsAJTefgJP/QvPyrWH1BD9wwlHOI0Q0cRz3aU4wml30WIRqt56geD7a+maewKxurH0HbWE8MQZlL/yLO62HSPdxIMGddtOlP/9K8b8I5FzToLaKqxLL4SN2+ChZ6F519CVODg4pBAgNm8WkfM+gAT8YAnel1fiX/oiatQZZXgw4winESKWsJ/rKji+zXuGq7eT2uceINw4g565SzAqaug882P4Nr5N8I0XUGPOdCOloAior7+NrF6PnLoIOflYmDoBvn4pvPkOPPECtHeNdDMdHA54rOpKoh85A3P6RAC0na0E7n8Kt/MCckjgCKfdZNttny247cRfLAMGI4UnncLHVdmO5GMqvPg8Agg+XaFKH3QOD7jtYXa+LKdwT8JPTzHsESypqODkRgpXEg7fSspJPOk0nu7onR0xPNOfRfZSsD01MXGnaSWdC9WM+lPraV1r2Y6MyTWzSPebAvi3rse7Yws9cxcTmnwk4amziYybQtnKf+Hd+DYKpTmJD7apdB+fQjPB720KOYUXiwJupvJkb5ecvABE45iPPQ8vvYXyoZNRFxwJ82chR82AV1ZjPvUi9PSX7Dhe7KcZ6TLpDHVmlFLH3ozGvbdOqQMl6nspDuDDcfweKgp4vkgvg07bUjQ9va4cJ/BUnuztIJqKnHwccsYScLkgbqA//RKe51/Dsiwi2XXK4H7ETEYgt+tPOYMn7kep8zPtvpTMY5H/njlI7tk9VBkrTxk1769aGsOJMn6g4winIdgXATABEi5OaI65aa+ixqNUvbEc3/vv0nPMKRiVtfQuPoPwlCMpe20ZelvLSDfx4KG7F7nzYcxlr9gC6ogpaEuORj1mNtaLb2I+/TIMOB3NDg4AMnEs1sfOhtG1ACgbtuC+/ynUjm6UPZrGx+FAwxFOQ7AvAmBCmnByrqd9gqd9J6OevIf+aXPpn7OQeF0DnR/8FL53VxJ481+oxqE/uepeY0crxi33okwej/ahk1Anj0f7wEK8i+dhLH8NY/lrEHECHjocprhdcM5JWCcsAFWBvgG0h5ehrnwHx/f70MQRTiNEUjipjsVpn6GIReDdN/FuWU/f0ScQmTST8Kz5RMdPofzFJ3E1bxvpJh5UyPvbMX57J8rMybaAGleP6+wT0E9cgLHsVcznXwfH6dXhcGLqBPj4OZAIMaC8shrlkWedKWYOcRzhNEIk4yI6byT7Hi08QMW/Hse7aS29i07HKqug+6yP4Xv7DQJvvOCELhgmsu59jPXvI3Nm4DrnRNT6Wlznnoz+geOILX+V+PNvOBYoh0MbrwflvFNQlxxtr3f1ot77BMp7m+11p2vukMYRTiOEmbI4OewvPDu3Uf3AX+g/9mQiM44iPOcYYuMmEVz+qOP7NFwErFXria5+D23+LPSzlqDW1eD50Mm4P7CQ+POvE3vudQg5PlAOhxbKnOnoHz0DJTm9zYtvwqPPocQcK9PhgiJygAy5OMBJ+jj19PRQXl6+W3X89U07Qmyt34WGMN0bQYCotyr1huJPzPqbHFUXSKxrhj0eQ0ku44ND7HNG2hmZo+wkMRw/Oc1Keh4rOVQ/sT44FUvCB8gYLCNZeVJ1JNOT05GkzdRuxezvVsK5PrlNsqZEMWODZbK3perImj4lfdRb9tQquVOUZNYZHTuJ3sVnYPkCYFn41ryG99UXUEyz4NQodj2Z28x41hQv2aPe0tqRnH7FjJlZeUqfcsXKGmpVbMqVnKlWUttz0wuNkMtNJwdTBBQFff4RuM88HrXedo6VSJTYv1YSXf4q0jtQtC1568xJKz1vIRtiqfsvVO/eKjtU+VLrGG6d+5JS/TWLja4baqRcKXmyR8Glo5K5rbQpVxLTXVWW4bnoTPS50wGQ1k6Mvz+B+v62nP3ljsQrvg6gJqdlSWzMHlWXJH10WqFpWIpNwZIzanmIqVfsPIXbUIhCeRY/+9yQZfeUvfG8LoZjcRohkvc3x6A7MniaN1P9wG30LTyV6OQjCB+1kOiEqQSeexytuWmkm3fwIYLxxjsYb76DNncG7jOPRxs3Gs/pi3Cfcgzx194m+swrmLs6RrqlDg7DQ1Fwnzgf77kno3g9iGliPPMKsvRFe34/5yZ+2OEIpxEi841YcK6+/Y8ajVDx/GNE3l9P/5LTsapq6Dv/03jeXIH3lRdGunkHJwLmqvWEV61HmTUVz+mL0KeMx714Hu7F84itXk/4qRWYW52o7g4HPtqYUQQ+9UH0iWMBMDc3Eb/nCaSl3RkRfRjjCKcRQ8HCNh0rIohzEY4Ynq0bce3cTt+xpxCbOYfogiUYDRPwP/5P1P6+kW7eQYuxdiPG2o1ok8biOW0R+uxpuOfOwD13BvH3thJZ+hLmu5tHupkODrnoGr6zjsd7xmIUTUPCEaIPL8d4aSWa85J72OMIpyHYVwEwwY7cqmKBWIC21+t3KB01FiW4/DFi2zbRf/LZmGPG0f/JL+Jb+gj6Jmei2z3B3NxM6I/3odZV4z59Ee5jZ+Oa3ohreiO+pl2EnnmZ6BvvDkaxd3AYQVwzJ1L28bPR6qoBiK9+j/A/nkLtSbxEObrpsMdxDi+Rve0cDjDFG8ONScwVwNLt6Vkc5/D94xyemTfhCB2364sHygmdeQHm6DEAuN98Bc+/luV1HHecw4c3BYopglpZhvfUhXiWzEPxuO30rl7Cz71O5MW3kHDEcQ4vsY7h1rkvOdidw9XKMsovPB3vglkAWN19DPzjKazV7+WtN3fKFcc5vJQ8h4JzuCOcSmRv/BFPrm8FoNpnC6dqJYxmxsDjR3H7AfDp9snmSwgmt5WYdy5mj0pKzUNnDMbJSYmoxHxzKcGUFDTRxJDwtGjZqW2J+iSeOZ9dScIpKVYSebJFUb40M54psgbT8winLLGTnW7lFUFm8TIZQiY7z6DYEVVjYOHJROccA4Da1oL/iQdR29sTeYYnoDK2WZK3jmwhld7G7DnqsoVUxn6y5q8rRRSVIq7y1VWs/uJlBMXvxXPCfLwnHYNaEbSPIxoj8soawstfw2ztLKm+/S2cRmIevVLqOdDYk3nnSslbqmBSM/LkKaNrBE9dSODs41E8bsSyiD3/BpHHnodIrGTBlG8k3lCCSVXSBU1+wZSdnk5SMBUSStkiKW9eNVn/0OKokBjKLpuPZNnjnnh2yLx7ijOq7hBGlMTJ5nRRHHAolkngxWdwbd/CwAc+iDWqnv6Lv4D3uaW41qwc6eYdEkgoQuSpl4g8+wr6/Fn4T12IPrYO30kL8J20gOjajYSXv4757iacAPsO+wLPUTMo+8hp6KPsyN/x97cT+vuTsKNthFvmcCDjCKcRJCWcxBFOByrubZtQ7voj4TPOxZgwmcjpH8RonIzniYdRok507L2CYRJ9ZQ3RV9bgmj4R3weOwX3kNDxHTsVz5FSM1k5Cz79B+OVVSNj5zR32HNf4eiouOh3PtEYAzJ4+Bh5chvH62yDFuxMdHBzhNIIMWpz2vuO5w95DHejH/8DfiB19HJElH8CYNhOzZhS+B/6O0ua8me5N4u9tIf7eFtTaKnwnLcC7aC56XTXlF51B8LyTCb+8mv7lr2M48aAcdgO1Ikj5h0/Bv3AuiqogsTgDz75K+KmXIBZ3Qgw4lIQjnEYQURIj6SwTEUFx3nIOWBTAs/JVtO3bCJ37UaS6htCnP4/n4X86o+72AVZ7FwP3P83AI8/hPvZI/Ccdg2tsHYGTjyFw8jFE3tnEwHNvEHl7Izhumg5DoHjdlJ++mLLTFqJ67QEJoVfX0P/gMqzuPsfC5DAsHOE0goiiDoa+tEzQnL/jQEfbtZPAXbcSOvejWOMmEP3oxVj/Wo763PMojiPO3icWJ/yvlYT/tRL39Eb8pxyLZ850vLOm4J01BaOjm/4X3qT/pbew+kIj3VqHAwzFpVN20gLKzz4eLWgPwIm930T3P5YS37LDsTA57BbOk3okURRE0VDEdITTQYQaDuH7+x3EPnAm8aOPIX7iB1BHN+B66J8w4Dy89xWx97YSe28rVFcSOGk+gSXz0GsqqbzgVCrOPZnQW+vo/9dKouu3jHRTHUYaXaNsyVFUnHMCeqU9qire0k7vQ8uJvLV+hBvncLDjhCMYgvQAmO+9994eDW9cu7MXAE8i5IBbU7DC/RCPoHt9uH1BPIn4TR7FdhhXk2EIoollItRAKmYT5MZtygolkAw1IHlCC6TiNyVjIhlZdaT5X2WnZYcwyI7nBHniNg0RlgByYzsNFZYAMkMTZObJH54g/Xt2jCcrEWcpO3yAnZYZDiA6Yw6R084G3YXS3YX7/r+j7mrJH8epQIynZCiDopMK54nxlFmnlVsmK8ZTsaH4uxOOYOgwAeSwNyf5NXUN7/wj8J+4APeksal0o7WTgRffIvTKaqwCkwvnG45RaiiBYuED9iQ203BCDexOjKd9QSldXKXEZhpO3qJhB3SNwKK5VJx9AnpNBQBmZw/9j71A+JXVaFK8ruIxmYrnLSUcQaFYTXZa/vADhUIOpOctFJupWEym7NAC2ev5QgwUDEeQZz+DZTLrWfDw0oJ59xZOHKcDhL3xR+QTThILI5EBVN2Nt6zCEU4cXMJJTMEcNdr2e6qsAsPAtfQJ1DdeR8ERTtnsVeGUlkcfNxrv8UfjP3Y2qs8OJiumReSdTYRWrGJgzQZ7QtYEjnDaOxwowklx6ZQtmUf5mYvRqxOCqbuPvidfIvLSytR/P1TsJ0c45eYrlieV9zATTk7f0Eij2n+BmMYQGR0OVLS2Xfhu+z+iHzofc9oM4uecizqhEdejj0A4PNLNOywwmnbR87cn6P3nM/jmzyJw/Dzck8fhmzMN35xpVA6ECb3xDgOvrCH2ftNIN9dhL6GW+Sk/6RjKTjkm5cNkdPfS99QKQi++BXHD8WNy2Os4wmmk0eyRdSJWwqowdARWhwMPJRrBc/89xI5ZhHHKaVhHziE2pgHtb3ehOiEL9hsSjRNasYrQilXo9TX4F87Ff9xstKpyyk5aQNlJC4i3dtLyy9sxe/tHurkOu4leXUHVmYsJLpmH6rZnYoi3ddH39Mv0v/QWGKYzUs5hn+EIpxFGUVRQNcQyscw4zl9y8KIArldXoDZvJ3a+HbLA+NJX0B/4J+q774x08w47jJYOeh9cRu9Dy9GnNxJYOAf/vJlgWY5oOkhxN9RRfeYiyo6bjZJ46YxuaaZn6ctEVq5zQlM47Becp/QBgKq7MGMmphEHfCPdHIc9RGtuwvPHW4hfeBHWpMkYH78Y9V8voD79NIozvc7+R4Tous1E122m6+7HU07DDgcJqkJgzjQqP3AsgSMmp5LD6zbT88SLRBKjKB1bvcP+whFOBwCarmPGwDIcP6dDBSUcwnXXHcRPPQ1r8fFYJ5yINIxFu/fv4Fg7RgyJxYnvbB/pZjiUgF5ZRtXxR1Nx4tG4Eg7fYlkMrFxH99IVGFt3jnALHQ5XHOF0AKDqdh+9ZcSdCOKHEIpY6E89idnUhHn+BcjkyRiXX4F2999Qtm8f6eY5OBx4KArB2VOpPPFognOmp0ZkGX0D9L60iu7lryNdPYAzn5zDyOEIp/1Imce+CbgTwzzciSGcLlUlqiiICJoZw+Vyocbt0VhKzA6omBOGIBkagGGEIUjvJsqeHy97XdVy07PT1My59pLDZjOqSZxiFpnWNI2EWFQT4QK0tOH0iTAEmjtRNiv8gKhWThnFzApHkMhjJW68g2XTwiskyitZ9aeG/2aFJ7C3JfetJNpm+1Qkh+kmQwCYiTKKpqBtXId12x+JXfgJpLYW88tfQnvqCdQVL6MwOOy31PAHGelaMkxB2rDmrDLJdS153Mn1NH+QQqEEiocYUBJpQ4UlSE9QCuQhCyVtW6H68++/WH3F2urKehAn87oSyWlnW579lRbKYDhli9c3sqKhlJFqhYRNTllFwTtlHIH5swjOPwK9siy1KbJhG/0vvElo5btgmOiAlnWfKRTXqVB6vm25YRByj2N3whEUCj+gurWMdSgcdiA75EBmmWS9aoEyhcMRZOfNDimQHUYgX9uy6yhGvvoOVhzhNATpATD3FYqioLtcxGMx4nEDl8u1z/blMDKoHe24b72F+LnnY806EvOcD6GMb0R78AEIRUa6eQ4OI0bdJR+mbNHc1LrZH6L/lTX0/2slZovTrXowI5qGFSxHysrB48G1ZeNIN2mv4AinIbjyyiu58sorUwG19hW6bgunWDyO33EQPyRRYjFc99+LsW0r5hlnIrNnYzQ0oP3tbyjNO0a6eQ77iH4xeNRoo1mizFPLOFWvSW3rkBibzBAKCuNUL3WKPQGtKYKJ4EI55Lvuw+s24587ndCq9YTefJfwus2pgJWHjo3i0EV0HbOuDqtmFGZVDVZlNVZZORIsR/yBVD4lHKLi1t+MYEv3Ho5wOkBwuVyEgXg8PmReh4MXBdBeeRll+3aMiz4O1dWYX/kK6qOPwopXRrp5DvsAFYXT9RpeMrsJYab8GDskxutmLyEx0VBoNiMs0ioZpbjZJVFesrrpFQMBJqo+jlcrcSU6XC0RNkqINonhVVTGKl7qFU/JbTqQfCn7Xl9L32tvg2k5fksHOKLpmKPqsEY3YI6uxxw9Bqu6dtBtIx/xGGpfH2p/L3KI/L+OcDpA0BPdc6ZlETcMSr8FOhyMqDua0X93s+00fuSRWOefD6PrUR5+BCxndOWhhA8Vv6LhRUVN801aafZhInxYrwPgIaOVdVY/o7Rq6hQPF+n1qbwPGq28ZfVxrFoJwOvSQ6vEqMLFLoniRs0QTklh9LrVw1tWL2XouBSFiFjMUcuYo5bRKlGaJYoHBQ8aZYpGNS50ZT/beYx95wbhsGcIIDW1GJOmYk6agjl2POi5skEJDaC2t6J1daB2daD29qD296L09aJEIyg4Pk4O+wBVVXF7PMSiUSKRCB6PNnQhh4MaJRJBu+dvWMefiHXG6bBoITJ6NPz1TpR+J2TBoUYIkxrFnbL0tEiUWWowtb1acRFHiIqFR1ExRVIWmAFMYgn39E2WbWlapFZSm+jaA9sKlXJGTiwXKOXM0oIIsF3CvCRdxBP1bJQQ66wB6hU3BkJQdI5Wy6lyOsgOa0TXMSdMwpw2DXPSVKSiMmO7MtCPtqsFtXUn2i77o0YGRni4wv7lgBZO11xzDT/60Y8y0kaPHk1LSwtgv1X96Ec/4pZbbqGrq4uFCxdy0003ceSRR6byR6NRvvnNb3L33XcTDoc57bTTuPnmmxk3btx+PZZS8Hi8xKJRwpEo5W7fAWNKd9h3KID6/PPQ0oL1iY/DpInwtSuR2+9EaXLmVDuUiIrgTbPmxLAIKoMvSBoKJrZvE0CTRHjF7GaLRFigljNftScr3Swhgug0S4T3ZIAxeGhUfCnRlI6iKPgT3Xs9YjBO8TFFsed0M0U4Si1jgeoEBD3cEY+H+MxZmNNnYjZOgvQBSoaBtn0r2uaNuLa9j9rVmSOSShlVdyhxQAsngCOPPJKnn346ta5pgzea66+/nl/96lfcdtttTJ8+nZ/+9KecccYZrF+/nrIye0jrVVddxcMPP8zf/vY3ampq+MY3vsG5557LG2+8kVHXgYDb40FRFCzLImaYeFwH/N/jsJdQ33sP+d+bkc99FupGwRVfQf75ILzy+kg3zWEvEcXCzaAlKSIW3jTrTggTHyp64rHUqPpoVH30isHjZhvdYuBTdGJYtEgUnxJAQ2Gl9KKgMFEpPKjEEIt1Vj/HqpVUKPZDMZKoR7MUanBRr3hw7e9uOocRQzQNa8o0zLlzsaZOz+iCU3q60TZvRH9/E9r2LSgJ39tC4QgONw74J7Ou69TX1+ekiwg33ngj3/3ud7nwwgsB+Mtf/sLo0aO56667uOyyy+jp6eFPf/oTf/3rXzn99NMBuOOOOxg/fjxPP/00Z5111n49lqDLvim5EiefnozBYURTeXwujVDMIByJ4BUVxbRP2FT8pkRYBCXdDyYZVymxFDNzPYlkx2pKJ1/cpvR0BiPHSCIttZ7cbiXyJuNGpVeTONWyYzApZmacJRiMZyKJOEpqdhynZHyltDKSnSc7RlNiv6SJUcvKKpP8beNGRnp6fCo1p0wyrpK9bsYTZbLiOtl5MuMrWZZk5FEHupFbb8E4/0LkiCPg4xehTpqI+sgjSDQRm8uSvHXl3U/qd5KcPOllAbSsPMn3TSsn7lG+2E9DxUgqHJMpN2++fNnxlbLzDieuUv5YTaXkLV6meBsBlLhQr7kIJM6nGtVFD3Gma7YFaKcR5QS9kgpt8BwVEXyKmwpLp1OJMU33o1sKDaqHE/VKXIrKC0YXm2SAyZoXv6JlOH6nfJ2MPqpVF9N1H8HE/qcpfnZaUbqJs1XCzFQCzNPKhuWgXYrMKlZfoedwoZhMpeQprWxm3tw4TsqQeQvFaoLC8ZqyYzWll8mO05QdoymfVWcwPlT+PMnYTMm6BTDHN2LMPBJjxhHgHRTbakcbrvfewb1lI2pn2+AZrYPqyQyRo6j5/7jsGFGHKge8cNqwYQMNDQ14PB4WLlzIz372MyZPnszmzZtpaWnhzDPPTOX1eDycfPLJvPTSS1x22WW88cYbxOPxjDwNDQ3Mnj2bl156qahwikajRKODgqa3t3ffHGAWPrfLFk4GVLrlsOo3dgAlGkX7291YJ56EdeqpyLHHYI5tQLnzLpSOzpFunsNu8Fy8k5VmHx0S55/xVoLofNY9htP0apYbXYQMC0HwKCqTVftBtt2KMEpx4VU0TBG2WhFmavbQ7vGqlzaJpaxDvoQVK3mvyNfF/7LZw0l6FRXK4C1/thZktmb7WHVbcf4v1sQE1cuoNL8ph0MDq6yc+Jx5GEfOzfBZUvp6bbG0/m3U9tZDzol7X3FAC6eFCxdy++23M336dHbt2sVPf/pTlixZwtq1a1N+TqNHj84oM3r0aLZu3QpAS0sLbrebqqqqnDzJ8oW47rrrcvyr9gduXUVTFEwRIib493sLHEYaRQTt+edQmpswL/oYNDQg//5vcN/9sGrNSDfPYZicqFdxrF5BVCz6xSQkJq5E3KajtCDNVpQYwgWuUfgUDQt40+xlhxVNWXQW6hVMU+27wRwtyKPxdh6Pt1Oh6GyyQizQyvGiEhWLOBbBhEBSFIXNZggNhSmqv6DfZKXqwkCwKGyBczi4EEXBmDIdY958zMlTB0MGRCPo69/FvX4tWvO2opEEHPJzQAunc845J/V9zpw5LF68mClTpvCXv/yFRYsWAblvV6XEJyklz7e//W2+/vWvp9Z7e3sZP378cA9h2CiKgs+t0x+NEzIE/wH9DznsS9RNm+B/b8L8xMdh4kTk05+CyS/Dw4+C6cT7OlhQEw7afkWjXDIv6DlaGXO0spwy57vqUt/jYqGhpLqEvIrG8Xol71thImIxX7NHzgE0m2HWWQOco9em7nHLjS7maEECaY7olgjLjU4CikZA0WmzYoxVvRmj9BwOTqyKCsy5R2McdTSUlafStW1b0FevRN+4HsUw0vyVnH6N4XJQPZYDgQBz5sxhw4YNXHDBBYBtVRozZkwqT2tra8oKVV9fTywWo6urK8Pq1NraypIlS4ruy+Px4PGMTDQln8cWThETTK20eaEcDk2U3l60P92Keeqp8IFTYPEimDAebr0duntGunkO+4F8Dtv1qod6NTdu02TNz2Qt0079SXd9yuE8lR/wKCrdYtAiMbyofMw12glAeRBjja7HOPEkrBlHQPJ/DA2gv70a95qVqF1OV//e4qAy0kWjUd59913GjBnDpEmTqK+vZ+nSpantsViM5557LiWKFixYgMvlysizc+dO3n777SGF00ji0jQSfuSEzIPqL3LYByiWhfrkUyh/uhX6+2HsWLjq32BS40g3zeEAoZgF3atoOUEtNUXheL2Ks1y1nO+q4yxXLT7lwBpl7FAa1uh6Yh+7mNiXL8eaOQsUBXXz+7ge+Af+m2/Es2ypI5r2Mge0xemb3/wm5513HhMmTKC1tZWf/vSn9Pb2cskll6AoCldddRU/+9nPmDZtGtOmTeNnP/sZfr+fT33qUwBUVFTwxS9+kW984xvU1NRQXV3NN7/5TebMmZMaZXegEtAVumPCgKUSFMsxpjqgbNiI/OYmuPSzKA0NyOVfhvsfgpdeHummOTg47GekvBzjgx9EZiXiFoqgvfM2+ksvoLa3AYVHvznsGQe0cGpqauKTn/wk7e3tjBo1ikWLFvHyyy/T2Gi/af/Xf/0X4XCYr371q6kAmE899VQqhhPA//zP/6DrOh//+MdTATBvu+22Ay6GUzY+HXpiYIhCTBS8I90ghwMCpbsbuen3yMcuQpk3F+XjFyLjxyF//ycYzlQtDg6HDfE4MnkKiIX6zlr0559D6+4Y6VYdFigiRQKZOKTo7e2loqKCnp4eysvLhy6Qh86+EFA8jlMqblMsRFfUImSAX7WoUQbsDMlYQ2ZarKREGTESy3hiW3I9Gd/JSHMoTsV+yoxhlIzBlIwBlBHXKbuenLhRibrS4jhJTtwjq2i6Xa1ZUtmicZyK1J/aT4E4TlZWHKeM/ZQYx6lofKUCcZysPGVytqXHcfrASfDBs1FUFdm2HetPf4XOrt2K41QoT2lxnHK35dueL0+pdRXPu3/2U0qZvVE290w9cHHiOI1sHCeZMxva2tA72/PmyY7jlLEtTxuGypvKswdxnGbd/vCQefaUvfG8LobjQHMA49ftkzNsKQVvxA6HMcueR/7wJ2RgAGXCeNT/ugqOnDnSrXJwcNhPqO+sRW1rHelmHHYc0F11hxpePfFGIckI1VlRwWEw6rZl4gZcCHFU+i0XFUoUZOj3USXRDZm0OKRHgE0iWb06g1HA7eHISsoipeaUScbQk0SkcEXNilSeFhhEybJCJSN6S1ZE72Q65EYTtxL1p+eBwYjiGfVlWaOKWZxyo4AnyrpdmevxQUtdTt7EMlmX5s5viYLC1qhClqh89UiWZcjauhn57U3Ipz6JMmE82uVfRF32HDz2JJKynGVZkbIsUulp2fvN/t20PNawVB1FrFODJKOMUzRvoajd+fPmy5U/mnlumXxRx4fab+llSik73LoOJEoZ8VvKSL1iFp+h8+YvW8xKNVSZ9P1nW5QGLTVq3nTItSzlWp4yrUnp9RSKBj4c61F2nsy2ZVmlsvIWqiP9uHL3P7QN5lAKrHnoHMk+4qabbmLWrFkce+yx+33figJlmv3w68eNdRDdUB32H0p3N/zuD/DCi/b6B06Gf78CamtGtmEODg4OhyCOcBqCK6+8knfeeYfXXnttRPbvU0x0LCxUBnCC0znkRzFNlIcfgb/8FQmFUCaMR/nmVXDsgpFumoODg8MhhSOcDnDSrU59eJwJERyKoqx9B375a2Tj+yheD+qnP4H6hc+C3zd0YQcHBweHIXGE00GAXzHRsDBR6cLviCeH4nT3wO//D+uRxxHTRDl6Lur/+w+YMG6kW+bg4OBw0OMIp4MARYEqwoAwoHjpc6I6OQyFCDy9DPmf/0XaO1Bqa1D/80qUU04Y6ZY5ODg4HNQ4wukgwacYVBIBoEfxM6CMzDx6DgcZTc1Y19+IvLUGRddRP3o+2hVfgIB/6LIODg4ODjk4wukgokyJERRbPHUo5URxjXCLHA4KwhGsP92Ode8/kbiBOmcW+ne/gTJj2ki3zMHBweGgwxFOBxmVhPBJDBSFNrWCOAf21DEOBw7y/EtYv/wN0rILpbIC139egfbR80B3ziEHh32BIRb/u2Mr/2xtYVVf70g3x2Ev4QingwwFqKYft8SxFJVWtRLT+RsdSqV5J8bPf435/EsAaGd8AP3bX0cZP3aEG+bgcGixPRrh3ze9S9SyiIvws80beaGrk0hySilntrODFidy+BDcdNNN3HTTTZimOXTmIUhFDE/OR2clQnGnRQNXkt+TSyUhipIRW0VQgVrpp1XKMRSNVlcto63ulHwajOidiCCemLtOsQYtC6ltWXPWDc5hl7i406JmJ8tkz3OXXYa0iNTJNMWVjMadv6ykzYm3J/PblVqHXSZ/npyI4q7By6TQ/HaF9mfGB0O0F4oqnjH/XNZxDRVdPF+E8pzo4ln7AYFHHsbauAE+9lHUsWNwf/c/kaeeQZ58NjUfYsZ+0//T7G0F5sBzZeXLzJP/oZEdfRzyRf8uLfp4sbzFyyS35a+rlLK5dRxaEcPzMVQU8WLbC0UKH84cdaVsz55fLkmhqOB5txWYby69jOZS6e0JMykY4DvTZ6BoCn6vzrKeTvw+F0uqalP5BvdTvN7s+eLsthSfzy4V3Twt4ndupPDi6+nlk7NS5GwvISq4WiDq+MHIoXMk+4iRDoBZCA1hFH2oYhFXXLSr5U6YAofhsfZd+MX/IKvfRtE01HPOtMMWjHfCFjg47ClbwiGiaS8cFzY0UON281pXN92xWJGSDgc6jnA6iNGxGEUfiggRxUOXVumIJ4fhMTAAf7kD6y93Iv39KGPHoH7z31E+ci64ncEHDg6lICLctH0Ld+xs5ol2e9Ld02treamri1U9Pal85zfU82pnFwPGnvdgOIwcjnA6yHFjUmP1gAj9WoA+NTjSTXI4GFm5Crnul1ivr0RRVdRTT0b9zjdRZh8x0i1zcDig2R4Jc/l7a4lYFpN8Pm5p3s5fmrZTrrv48oQJ/GzjhlTeSYEA43xeVqeJKYeDD8fH6RDAT4xK6adbKaNbK0e3YvisyEg3y+FgY2AA+ctdmK+9gfrxC1FqqtG/+kWsN1Zh3P0P6B8Y6RY6OBxwrPv/7d15fBT1/T/w18zeuTb3BSHghSBI5VBQK4LIYQGPtmr9FaW1aG3AUrXWoxatVdR6tJXUr22t1NseotaDigcgRRQRCooiKDe5j92ce8x8fn/skezubHY21ybZ1/PxyGMzM5/PZ96fmU3yzuzM59PagmkZmfjB8BEAgEavF78/tB8lNhu+N2w43q+vx4rPv8APSkfADRXH2ttxQhr/wR3MeMVpiEgXbUhTmgFJQp0pBx6JOTF10+49UO99EOrb6yEUBfKkCTCt+AWkiRMSHRnRgHPM7cKO5qbgcp7JjONTUvBObS0A4P4xY2E1yPjTgQP4zedf4tLhw3BiOhOnwYx/XYcICUCW4oBHMsElW1BrykWBcgwy73qi7nB7IF55HcrW7TBeeRmkYcUwXXMV1E92wvvCvwBnU+w2iJLAnKxcfOBowG8PfIUiixU7mpxYWjoKd+37ElsaGjA1Kws3n3gimrxeWIwGpHDctEGPV5yGEAlArrceBuGFVzah3pzLtIl65vBR37hPr7/lu/o08VSYVvwC8pmnJzoyon6139WGdo1hOQotFtxWejyOT0lFs+LFspKRODUjAxPSMzqGPwCQaTIxaRoieMVpiDFARY67DtXmfLQZU9GkupHuqUl0WDSYKQqUf6+F+slOGK68DHJpCYxXXg7p9EnwPvMPiGPViY6QqE/9ueYoNrU0YsWw43CcITVkrCghBEZYbShNTQ2u29PSjBq3G8UW35yiUpRxqmhw4hWnGMrLyzF27FhMmTIl0aHoZhFuZHkbAAAOUybajKkxahDFJo4eg/f+38P7r1ch3G4YTj4R5hU3w3jxBRy6gIYkVQg8XnME9YoHfxx5Mk6wpkQMsNk5KWpWvHji6GHc+eWXmJ+fj2E2W3+HTP2AiVMMA3UAzFjSlBakepyAJKE+pQge2ZzokGgoUFWo69bD8+sHoOz6HJLJCNOC2bD85lYYppyW6OiIepVbqPjK1YZrc4ch3WDErtZmfNLsRIPXo1leCCDFYMBj48ZjTl5+P0dL/YUf1SVCF5dthRQYIt8/JUp4Adn3GXtwuhYAQglMy+I7nYFpW7LVVnhVi+9m8bThKHBVwYDQqU6kwJQn/mk2JP8vhMAUKJKxI+ESYVOrSGrotCPhbfrqRNlflPWh7YVN6RIWR+dpQOKZniVAjTFdSmB96HQmXe8nMMVKYNmgxJ5KRo3Slu97fdOzdJ7CRO/0LJ3rhE+BEm1/QW3NEE88CfXUcZAuXgA5Nxvm666CmDMdyj9fhThwKKLfoe2HT9MS/W688HOnNYVLtDa0pnABOk+XEvmzGGvqk67q6m1Du73BLdq0KR3b46+rZ8qVWNOnhG6TNct0nkbFK1SY/PciyaqKAqcVVSYV/2qowO6WJtiNJpibZZyTm4PZefmQDBL+U1WNU3MyMTwrDT/MS+to3z8lisFsCNmv1jQq0aZPCZ82RTZETtMSMW1KF1OkBPcta0+10tWUK9HairY8VDFxGsIkALlqI6qkHHglI2rMech3V/NJO+o9Oz+F+HwPxIxzIJ03A9JxI2G8+XqoW7dDefVNoK4+0RES6fZ4xWG0qAqyTCb8sGg4BIBWRcH7jfUwyBL+OHY8WlQFb9RU438OJ87OzkGLV4FLVZDFj6uTRnKkh0nMAIE8Tx1kocAjm1FrzoWqY8JRIt08HuCtd6DefT/UDz6CUFXIU06D8Vc3Q75kAZDGe+xoYPuqvRXL93+BBq8HF+fk49mqY1hbVwOzLGNBXj42NNTD4fXCrapIMxqRazZDEQIpBgPyLRYsLCpCqpHXIZIFE6ckYIIX+d46SEKFS7aizpTDa07U+xxOiOf+Ae99v4P6+ZeQTEYYZk2H6Te/hOHCC4DUlERHSBTBI1T819mAhdl5uKXkOIyypuDi3AJkGk1QhcAZ9izMy82DW1XxTl0tjrW347XqKhRaLRBCQAyRj1pJP6bIScIsPMh116LWnIt2gw315lxku2t57Yl635FjUB79E9SxoyEvmAu5tASGebMgzzgb6obNUN5eDzRyAE1KLK8QMAnAJMn4dk4hUg0GCCHw6LFDeKOhBvVeD9bW1+Ca4aW4smg4/utswDanA/+pq8E52Tn47rBhADjUQDJi4pRErMKFHE8dak25aDWmwSC8yHTVJTosGqLE7j1Qdu+BMm4MDPPnQB4xHIY5MyGfcya8a9+Bsm6D72M+on7255qjqPN6YDcZcW3+cKT6b6L+ur0NpVYbXj91EoySjF8f2Ic36qqxZNgInJudg7OzsuGB7yM6Sl78qC7J2NR2ZHt8N+w2mTLRYuScSdS3xM7P4L33YXjK/wL14GFINitMF3/LN4TBmad3+ZQpUW8RQsAtVKysOoAqrxuXZxfgqNuFhysOYl9bKwDgeFsKLszJh9H/dPNZ9kxUuV3BNoySxKSJmDjFMhgHwIwlVW1FuscBAKi35MElWxIcESUDsWs3vPf9Dt4nnoGoq4eUnQXTD74H84qfQz5tfKLDoyEu8JGaSZJwVU4RRlpsWDHsOKgQ2NzUgBqPO6T8UVc73qyrxTfSMhIRLg1gTJxiGKwDYMZi9zTA6m0FJBm11kJ4JX5qS/1ACKhbP4Hrlyvh+ccrEC2tkIcVwbLsalh++TPIp4xOdIQ0xDQpHWPetSgKvna3weUfW8wiy5ifmYe97a045GoDANR63Hiy4gju3L8P38zMwvy8goTETQMXE6ckJQHIcdfApLigykbUpg6HIvHtQP3E64Xy1nq4bvsNvK+vg2h3QT6uFJYbr4P5F8sgn3xCoiOkIeCJhmNYcmwPvnL7kqIsowmnWtOwquZwsMy4lDTkmcx4x+G7hcEqy3ALFb857kRcmMukiSLxL2USkyGQ214Bg+qB12BGXcowjvFE/au1Dd6X30D7L+6G9z/rITweGEYfD8vNS2G+eSnkk09MdIQ0SP27qRbHPG6cbEnB3x3VaPXPOPDDnGI4FC/+1VAVLHu+PQceofrGaTIYcW3xCBSYeQsDaePnM0nOKBTktVegyjYMbqMN9eklyGk6FLsiUW9qaobnxZfh+c97MH1rFgznTIPh5BNgOPkEKPv2w/vaOig7PvNNBkakw8zULJxuy0CB0YwfH9uDd5obsDArH2ZJxk0Fpbjj2FewGAwYabHib7UVOC01A+YoU4sQdSYJjt6li9PphN1uh8PhQEZG390s6G6s9n3jn28O/s/nA/PPBedtEx3zdwXnrRNh89iFzwfXqU6wvv+1XRhQY8yFkCSkKC3IdtX4rj3pmEOuYw47NXS5c/0o89pFlEPkPHaIMr+dntg62uz+/Haqjvnmwue105ojL9a8dkKJHWO0ee1C56rTN6+dGjJXXdh8cFHmt+tMKKFlwvff0e/IefSi7TfYVqYd0vkzIE07HZLJN5WFOFYJ5e31EFu3A4oSMQdetDnstGLUmrswXDxz4kXU7cav1cE0Z12seek6C59LrjMpygR24fPMac6fFlYmfN45wHdMTUYZHzgb8VjFIawYeQJG2mwwGY14r74OX7W3YkeTA9Ozc3B58bCIeeI6fx/YZjCFzukW2K9s7njariOW8Lnqos87J5uNEfH76oTW7TzHXPj8cpF1o88/J8eo09W6aHPXyTqeOMy67r6YZXqqr/9e84oTAQCswo1cbx1qjDloNaRCNinI9NTzgztKjEYHxD9ehlj7DqQZ34R09lRIxYUwXnk5xMJ5UN/bBGXjf4G29kRHSgOYQZIghMC0jExsdjbgb1XHcPeoE1HjduGszCycZ8pFq6Ig1cQ/haQfr0tSkE24gmM8NZsy4DRlJjYgoqYmiFffgPqre6GseQ2i0QEp0w7Dxd+CaeUKGC69GMjLTXSUNAjcOHwUqt0u3PLVHlz52f/wkbMRADguE8WNaTaFSFVboXpkNJqy4DRlQlIVZHgaEx0WJbv2dqjr1kN9931IU06D4bxzIA0rhmHmNyGfe5ZvjKi3N0D9Yl+iI6UBJjB+k1eoaFEUWGUDHh8zHiOstgRHRoMVE6cYysvLUV5eDkVRYhceItKVZggBOMxZcFhyIEEEB8wkSihFgdjyMbxbPgZGnwDDzHMgjx8LacI4mCeMg3rkGJT3NkH5cBvgcsduj5LGc9UVuCAnD1cUFGvex0OkF28O12mo3xyudQO2w5QZ/Lgus70Gaa4GzXZ5czhvDg+p2xs3h2uIqBs4TgX5MMw4G/LUKZCsvkfIRWsblA+2QnnnfYiaWt4c3ssGy83hnderQsBg1L6xO3ATN28O583hevCKE0WV4WmEEAJN5iw0WvMgKx6keJsTHRZRqKpqKC+8BM9Lr8Nw1hkwnHsW5II8GM87B4YZZ0Pdvgve/6yHum9/oiOlBOoqgSOKBxMnikoCYHfXQ0BGs9mO+pRCGFqOwKK2Jjo0okht7VDe3gDlnY2Qx54Ew8xzYDh1LAyTJsAwaQLUrw/Cu24DlI//B+i44kREpIWJE3VJApDproVXNqLdmIq6lGHIazoIk8r7R2iAEgLqZ3ugfrYH3uJCGGZNh2HqJMjHlcJ87ZUQlzbC8+5/4V2/GXA2JTpaIhpkeIccxSQByGmvgsnbBlU2oDa9BAonBaZBQByrhPepF9F+013wrHnDN5xBVibM3/4WbI/cBfNPFkMeexLAj3GISCf+9SNdZAjkth5FTeoIeA1m1KSXIM+xHwaRPE8b0iDW1Azvv9+C9413YJjyDRjO+yYMx4+E8YyJMJ4xEWp1Lbwbt0DZ9BFEQ2OioyWiAYyJE+lmECpyW46gOq0EXoMFtRkjkOc4CAlMnmiQUBQoW7bB+9+tkEYMh/HcaTBOmww5Pxfm78yHuOQCqJ9+Ae+mj6B+vBNwexIdcVIxjDkRyuFjQBvvo6SBi4kTxcUovMhrOoya9BHwGFNQmzECOY1fQ8bgeZSaCADEoSPwPPUPeF54GYbJ34Bxun9i4VPHwnDqWJjb2uH9ZBe8H/8Pys7POS5UH5Gy7DCefTpM06dCLsyH+6U34X3lzUSHRRQVE6cY+nsATGHyj2YbGIfIEBiXKHSspsD4Tv5KIXVE+HhO/mXReeyn4PhNgfGcREhZKXycp07rLKqCfKUeVVIu3KZUNGSfgFxXVce8dj2cGFgK2xbZlsY4Tl2M8aS5vXMsMWLsydhPeurqGftJdXv9IYaX7Z2xnzoOSddjQHUWPh5UtLGgtOt0PX5TV7EFy6jh40d1NclvjPGiPt0JfLoTSm4OMGUS5NMnQsrNgemsKTCdNQXC7YH44kuoOz6F+unuqDeVdxVDZEy9/8+GnvGpenPwx2jjMHUWPiYTAMBkgjxhHORpUyCNOSk4LpBob4chxQxDhiXqfsLjDx8zqWN97DGZgsv+uoG2DebIsZJijc1k6DTfXfj4SbJ/ourwsZg690X2148YMym4/8hxl6KN09QxTlTkuEqx6oSX62pdtPGcQssMvSltmDjFUFZWhrKysuCAWuRjFh7kuWtQY85Du8GGWnM+ct3VnBSYBrfaOqiv/Qfqa/+BNHIEpNNOhfyNcZDyciGdegrkU08BAKj7D0Ls3gP18z0Q+w8BSTSzQLfJMqSTT4Q85TTI3xgPydYx5Ym692uoH3wE8clOwO3WTraIBggmTtRtFuFGrqcWNaY8tBtTUIc85LhrmDzRkCAOHII4cAjqmteA4sJg4iSPHAF5VCkwqhSGb82GcLkg9u2HumcvxOf7IA4dibhamrTSUiGfOgbS+LGQTzk5JFkSdfVQP9wG9cOPgZq6BAZJFB8mTtQjVtWFXFc1ai35aDOmog5ATnslkycaWo5VQj1yDOob6wB7BuRTxkAacyLkk0+ClJ4G6ZSTIZ9yMowARFs71K/2Q923H+LgEagHDwNNSTLivkGGVFoCwyknQx53su+qXaePc4SzCeon/4P68Q6Ir/YHJ+AlGkyYOFGP2dS2kOSp1lqI3PYqSLxhnIYihxPq5g+BzR9CkSRIRYW+j6BGnwDpxOMhpdhgGDcGhnFjglVEQyPUg0cgDh2BeugoxLEKiKq6jnsLB6v0NMjHlUIeVQrpuFLIx5VCsoTen6QeOgLx6edQd34GcfBwaJ+ZONEgxMSJeoUveapCnSUf7f7kiVeeaMgTwpcEHauA+u5GqCogDS+CfOLxkEeOgFQ6HFJhPqSsTBiyMoFvjOuo6nZDVNVCVFZDrayGqKqBqKqBWl0DNLUkrk9ajAZIuTmQ8nMhDyuCPKoE0sgRkHOzI4qK5haoe/ZB/fRzKJ9+AbnJmYCAifoOEyfqNTa1PXjlqd2YglpbEXKaj3CoAkoeQkAcPgbl8LGO0c0sZkjDh0EuHQ5pxDDII/zJlNkMqaQYKClG+HNHor0dotEJ4XBCNDgg6hsh6ht8r40OCIfTl1x5vegVBhlITYGUng4pNwtSXi7kogJIhXm+G+OzMzWfoBKqClFR5btZfv8h38eTFVWhV5V4ozcNMUycqFdZ1XbktR1Dja0ILoMN9SnFyGk9yitPlLxcboiv9kP5an/HOkkCsrMhFeRBKsyH7H+VCvIg52RBslohFVqBwvwumxYuN9DSCtHaCtHSBtHSCrS1+29OFwAkQJZ8r5LkuwfJZATMZkhpKb5kKTUVUoqty/0A/mSuus53dezAIagHDvtuhHe5enR4iAYbJk7U6yyqC3ltFaixFaPdlIoGWyGyWo4xeSIKEAKiuhaiuhbY9Xno2PsmE6ScLEgZ6ZAyMyBl2SFlZ/mu+mRnQbJnABnpkIwGSBaz74pWdmbPQ1JVoLUNoq4eorYeakU1hP8jRLWmFnBG3uCuZxwnoqGGiVMM/T0A5lBhUV3Iaa9CrbUQreYMmL1tSHM1JDosooHP4/ElLJXV0ctIEmC1Qkq1ASkpkNJSIKWmQFitkGxW33ZJ8n1kJgSgCgACUFTfoLYut+/qVEsrRHNL8PtBf7M6UT9g4hQDB8DsPpvSiuy2SrSa0pHqakx0OERDhxBAWxtEWxuA+uBdhP09cjhRMmLiRH0qxdMEm6eJQxMQEdGQwH89qM/xLggiIhoqmDgRERER6cTEiYiIiEgn3uM0wKgGEwBAkv2nRqiBDf5l/71Chk43gare0LJKYNkUur7TEzPCv04Kbgt9jdgOdExcGv7kTbS2wusDgKr4i6ghy8FXACLwBGP4RKnBupHbhapolg2U6WhTiSwTbEM7JtE5jvBt4bGGxdi5buDG3eCrqmquVzs9wRmtjOLxaq7vst2I9ZFPikaPKTTmkDrBsiKkbEebHe8XVQl974TfzCyUyHvhAnXCy6pqeFtadbVvltYqq7WuyzZU/ffthfdbDz03eve1eG4kl2MMTSDJkdvlsPbDhzfQGu4gUCewTZbD68gh8XRuI2Kbv67BbAjZ3nmwz8D+ZLMxZJtsNvn3H7q9c5lAe4awusH9dOp/sF2DQbNMeF2tPocPUirJhkAnIuogvKwhrIxmndB1WoOiRqXV3iDFK05EREREOjFxIiIiItKJiVMM5eXlGDt2LKZMmZLoUIiIiCjBmDjFUFZWht27d2Pr1q2JDoUAVNZy9HEiIkocJk40aHz06R78YMXDePGtjYkOhYiIkhQTJxoUhBA4fdxoXD53Oq6842Fs3/N1okMiIqIkxMSJBjxFUSBJEpzNrXj2jXdx1fzzcPywwkSHRURESYiJEw1oqqrCYDDA7fHgwp/dheK8HNz/0x8gIy0l0aEREVES4gCYNGDVNjiQm2UHACy8fgUKc7Lw4A0/QlZGGgDgs68Pwe3xorq+EbPPmJDIUImIKEkwcaIByeV24yf3PIpzp0zA5h2fIdVmxUM3LkFuph1CVfDc2vX49ePP4exvnIL/7d2Ptz/cgQeuvyrRYRMR0RDHxIkGJIvZjHuu/wHm/PhWuD1eHHjzaRj9Hyw/88Z7WHb//+GZu2/At86eAiEEzvrhzdiyaw+mjh+d2MCJiGhI4z1OMXAAzMQ5ccQwbPrbI8hKT8Ptjz4JANi0/VNcd285/vZrX9IEAJIkISMtBS3trkSGS0RESYCJUwwcADOxivNy8O4TD+D4kmIAwO+fewV3XnsFFk4/I1jm0Rf/jebWNsyYNA4AsPdwRUJiJSKioY+JEw14eVmZuOY7FwAAUm1WnP2NU4LbnnlzPTbt+BzLv7cAsiyjzeXCS+9uxg/uXpWocImIaAjjPU4DjFsRAACD5FuWJQMAQDL6X4VvO4TaUUk1+Nf5tkmy17+ohpYNqdPFtk5tiU7rpahl1YiyEXWC+xOadSUd+xGKF15Jxuo3N2LamVNx3+NP4f2Pd+LC887C9G+eBYM9G2mqgtuWXo0f3fFbfPeXv8O/Hv5V6P5UpdOulZDYRKdtIWXVjjiEooRtU/x11a7Xa21TwtoP1ukcoz82Jew1ynpfM4quuqrSuV/hxzr2foL7iyij0ffwdsPaUcP3r3YVW+h7KLyur354mdBlrX6EtxutbldtRNRVtet2JVociSQFfiF1QZa1y0iG6P+fy2HtBvYTXqdzOcm/HzmsTHhdSZYjygW2ySajf9kQuj5Qx9zxpzHQTqCMwb9NCisbWPbt06C5LRhbWIy+7w2Byv7XQNmw9Z3513VuJ7SNyDpRy0Zb7hxDxP7juPaiFf8gldArThs3bsSCBQtQXFwMSZLw8ssvh2wXQuDOO+9EcXExbDYbzj33XHz22WchZVwuF5YtW4bc3FykpqZi4cKFOHLkSEiZhoYGLFq0CHa7HXa7HYsWLUJjY2Mf9456myRJeO7hu3CkqgaX/vQOPP/a27jxB5fhu3NmoCA3G16vgo8/2wMA+POdNyA3y47n33wvwVETEdFQktDEqaWlBRMmTMCqVdofqzzwwAN4+OGHsWrVKmzduhWFhYU4//zz0dTUFCyzfPlyrFmzBi+88AI2bdqE5uZmzJ8/H4rS8V/7FVdcgR07dmDt2rVYu3YtduzYgUWLFvV5/6hvvP6n3+Kv996KLS/+H2ZOnYgsezoAwKsoWL1mLZb95g8AgP/t+RoNzuZEhkpERENMQj+qmzdvHubNm6e5TQiB3/3ud7j99ttxySWXAAD+9re/oaCgAM899xyuvfZaOBwOPPHEE3j66acxa9YsAMAzzzyDkpISvP3225gzZw4+//xzrF27Flu2bMEZZ/huKP7zn/+MadOmYc+ePRg9mo+vD0ZpqSmAqmDP/kNQVYExx5fCajHj/huuwSkLf4gGhxMAsOQS3/tr/db/YcLo45DFEceJiKgHBuzN4fv370dlZSVmz54dXGexWDB9+nRs3rwZALBt2zZ4PJ6QMsXFxRg3blywzAcffAC73R5MmgBg6tSpsNvtwTI0eL2xYQt+eu8fIPz3TqWm2DDtG2Nx7WULsPmZ38Pkv5fhcFUNLvnZr3G4siaR4RIR0SA3YBOnyspKAEBBQUHI+oKCguC2yspKmM1mZGVldVkmPz8/ov38/PxgGS0ulwtOpzPkiwaeny2+FJNOGY0LrrkZR6tqcPdjT6OxqRmnHF8KWZaDCdWi+bNwyayzceUdD2LvoaMJjpqIiAarAZs4BUhS6FMXQoiIdeHCy2iVj9XOypUrgzeT2+12lJSUxBk59ZeVN1yDaaeNwy0PPY433v8Q1122EDaLJXiOvV7f/W7LvnchLjnvLEy/+mZ8dYRjPRERUfwG7HAEhYWFAHxXjIqKioLrq6urg1ehCgsL4Xa70dDQEHLVqbq6GmeeeWawTFVVVUT7NTU1EVezOrv11ltxww03BJedTieTpwHsVz+5Ci2tbTAZJJhNpo5H/YWA0WjAJ5/vwzsffoI17/wXI4sL8Om+AziuOPJKJBERUVcG7BWnUaNGobCwEOvWrQuuc7vd2LBhQzApmjRpEkwmU0iZiooKfPrpp8Ey06ZNg8PhwEcffRQs8+GHH8LhcATLaLFYLMjIyAj5ooEtNcUGs8kUsu69j3bgwdX/wPylv8SBY1X4wYWzsflvD+PCc6clKEoiIhrMEnrFqbm5Gfv27Qsu79+/Hzt27EB2djZGjBiB5cuX495778WJJ56IE088Effeey9SUlJwxRVXAADsdjuuvvpq3HjjjcjJyUF2djZuuukmjB8/PviU3ZgxYzB37lwsWbIEjz/+OADgmmuuwfz58/lE3RC28eOdePLl/+B/e77CaSefgD/9ajnmT58aHGRSVVXEHtKPiIgoVEITp48//hgzZswILgc+GrvqqquwevVq3HzzzWhra8NPfvITNDQ04IwzzsBbb72F9PT0YJ1HHnkERqMRl156Kdra2nDeeedh9erVMHQa6fTZZ5/F9ddfH3z6buHChVHHjqKhwWg0oLGpGc+uvAUjivKRarOGbJdlGUJVoKrqwL3sSkREA44kRPgcGKTF6XTCbrfD4XD06cd2jpY2AJ2nXPFPI+Bf1ppyRVK9/nUibDn2lCvRp1HR2E+MKVci1neuE8eUK1Hbiza9iVYZRYGiKDBIki85CkwNoDHlSktLKxqamtHoaMIpx4+IaItTrnDKFS2cciUUp1zhlCtd8rdrnvZt/XW6qa//Xg/Ym8OJekqWZUCIjqTJr/MTlbu/PoRFv3wIJ48cjp179+PaS+ag7LsXJCJcIiIaBJg40ZAlSVLwCldtgwP1ziacVNLxhOam7Z/hpkeewPlnfAP3LbsKXx85hoU33IPxJ5TinNNOSVTYREQ0gPH2DkoKAgIP/u2f+OpIBSRJQk2DAzc98gTmnT0Z9y27CgAwqrgA3zztFJiMQ2cWbyIi6l1MnGIoLy/H2LFjMWXKlESHQj2Ql5WJ3/5sCY4f7rviVHZvOY4vKcKKa64Illn34Q78e+NHyExLBQA4mlsTEisREQ1cTJxiKCsrw+7du7F169ZEh0I9ZE/3JURCCGTZ03DLD74b3LZl1xf49V9exO0//C7GjCrB/mNVuOQX9+MXq55OVLhERDQA8R6nAUbxP4Ujgg+R+J+U89/MHHhwxSB1fJwk+T9aCjxxJ4RvEMjwp+20nnqL+uSdRlkR/iReRBuRTwMJrafmNPYjNLZFPJHX0WjUGKVY+1MVeDxeVDQ04R8bPsaE007D+o+24/r7HscVC2bhusXfQ21dPW557HlMGDsan3y5Hz9d9RwevW1paCzhT8Cp4U+FhT1916mM7ifzutoW7cm88PqdYunqCblYT8+FP6nXVdnwp/a0n2ALfWovfL3WOjVG2eCx0Igleh2NpwSjPC2n9ZRgV+2Ebu/+k3LR4ulL4U+shdPzlF34k3ERT3Np7Cf8qTOtutGehIv2xJqk9VRdlDKBp+3kTk+RRT6JF1bHaAp0plPH/GVM5pBtwX5oPPUW68k4zSfbotSRIp6Uizz2kWW6uEUhytNzUZ+2G+J4xYmSjslkxDMP/grvbfkE31l2O8rufAjXXH4hll/puwK1YtVqVNc34Lc3XIN3nvgtqusbcMejTyY4aiIiGgh4xYmSUmZGOt576g+oqquHoqoYOawoeKVm8cVz8dOVq/B//3gNy664CC8+eAfe2vwx3B4PzEn6HxYREfnwihMlLbPZhJKiAowcVoSmllY0tfhuBj/j1DH45yN34qW338e23XsBALPPnBwxDx4RESUfJk5EAB59+p/4fz+/O7g8vDAPBdlZUDTumyEiouTFxIkIwG0/vhL52Zm48hf34GhVLZ577R28uekjpKemJDo0IiIaQJg4xcBxnJLHX37zC2RlpOPGB/6IO8tX44lf34Qxx42IeJqLiIiSF28Oj6GsrAxlZWXBSQNpaPv97dejtaUFbo8X9lQbhMZcd0RElLyYOBGFSbFZkWLzjQsUGD+LiIgI4Ed1RERERLoxcSIiIiLSiYkTERERkU5MnIiIiIh0YuJEREREpBMTJyIiIiKdmDjFwAEwiYiIKICJUwxlZWXYvXs3tm7dmuhQiIiIKMGYOBERERHpxMSJiIiI+oQKCR7ZBJdsTXQovYZTrgwwXlUAQHCqDzkw44fwr/cvKpII1gkvK/nzYdlg8b2io2yQUENfI9Z3o47Wsv/7YGv+CXOlLuuI0DpCTx1Vs05AsK7WhL2BvobVkaP1t/M6VfEv+pZVSHAY7bC7633H3b8d8E3hohlDsI3I7R3rQrcF2wprQysmhLUhtI5BeJmI9sPa6hxbcDkQW+hrV2WC68OWVUWJqBOtbFf7UWOU1aobbb966sRTJmAgTyIda55GyRB9e/g2SaOtiDJhy7J/uXPdQBkpbJtsMGi2JZuMMetKsiGwQ/9rp/2ZzCHbgrEElo2m0LqdtxnC2gu2obGf8BiCbYUfo7DtIXFrbIP2sY8oK3VxrrvaBkBIMhRIcMMIDwzwwAAFMrwwQPHXlYWKYWjssp3BgokTUS8RAGrNuXDJVngkI3Jd1bykS0RDkgoJLhjRLpnRDhM8UvR0QhIqDFC1/oUflJg4EfUSCYDd40CN2QyXwYZaSwFy245pX/EjIhpEBBBMlFwwwQUjEDYJulEoMMMLIxQYoQZfDUPsdyATJ6JeZBFu5LlrUGPOg8tgRY2tGLltFTBg4H4cQ0SkRQDwwIgWyYZWmKFIoR/vGYUCCzywwgMLvEMuQYqGiRNRL7MIN/Jclai1FMDtT57y2isgQ4ldmYgowTwwoFWyoEWywtvpIzhJqLDBDavwJUtG/z+EIsY9UEMNE6cYysvLUV5eDiX8ZlmiLlhUN/LaK1FjKYDHYEG1bRhymw/DKLyJDo2IKIIKCS2SBS2SDW7JFFwvCQEb3EgRLtjghtRFG8kiudLEbuAAmNRdZuFBfttRGFQPvLIJ1Wkj4JHNiQ6LiCjIDSPq5QwclXPRIGf4kiYhYBVuZKtODFNrkSuakMKkKYhXnIj6kEl4UdB2FDXWIt+Vp9QS5DUfhllpT3RoRJSkBIA2yYImKQUu2RJcbxJepIo2pApX2H2ZTJk6Y+JE1McMQkF+2zHUWAvhNtpQk16CvCYmT0TUvwSAZjkFTkM6lMC9S0LAJlxIF22wwMMUSQcmTkT9QIaK3JYjqE0dBrcxBbVpw5HXdBAm1ZPo0IhoiBMAWuQUODolTLJQkSZakaa2+m7yTrIbvHuCiRNRP5EhkNtyFDWpJfAYrahNK0F+00EY+LQdEfUBAaBdtsJhzIRH9t3wLQsFdqUJqWjnTc7dxONG1I9kCOQ2H4ZBcUMxmFGbVgKV/+kRUS9zSWbUmPNQa86DRzZBFioyvQ4Ue6qQrrbwj38P8NgR9TODUJDXfBiy6vVdeUofAZV3FhBRL3BLJtSYclFtKfBNrCsE0pUmFHsqkaE2cyaDXsDEiSgBjKoHuc2HIakK3KZU1KWXQDB5IqJu8koG1JmyUWUuQLvBBgiBVG8zilwVyFKcTJh6EROnGMrLyzF27FhMmTIl0aHQEGNWXMhtPgJJqHCZ09GQXsJfbUQUFwHAabSj0joMrYZUQJJgU1pR6K5EtrcBRt5D2euYOMXAATCpL1mUNuQ0HQKEijZrJhrThjF5IiJd2gw2VKaUwGHOgpBkmFUXClyVyPXUwcRZCvoMn6ojSjCrpwXZzUdRnzYcrbYcyIoH9pbKRIdFRAOURzKi0ZKLdmMqAEBWvcj0NCBFtPMD/37AxIloAEhxO6E2H0Vj+nA0pxVCFgrSW2sSHRYRDSAqZDSZs9FkygQkCRACaR4H7F6H7x4m2ZDoEJMCEyeiASK1vR6qZIAzrQjO9GGQhYKUFiZPRMlOAGg1ZcBhzYUq+/5sW72tyHTVwiQ8TJj6WVyJ0549e/D888/j/fffx4EDB9Da2oq8vDycdtppmDNnDr797W/DYrHEboiiCs4OJHx3uij+G17Cb0aTpM4XZH2FZP8qgxQog5D1IfUDLfrHEJL9haWwspKIvONGCDWsjH85bH3nfnQsq50ijlVX1awDNbKsFF4/WFeE1tXYT7QYRBexRexPDd1fRBzwjdSr2V5gWVVghwrhdaLJmIHG9BLI1lSkqq2AqgTLAIBQAssaxy1QRtUuI8La0iyjaJQJWxbh+9ZqN7xcF+2F7DekjKpdtqv9BNsL71f4sdA4flHKxlrfmarVjzh1FVtfkGT9t7zKhuh/qCWD3PWyxn7iqSMFkoTgq+yvExZTeLnO7YRtC18vmTpNxh2+LXy/RnNouc6xaMSgGYdvJcIKAQA8MKBeSoNL8g1gaRQKMtEKq8EDKTUjrE7YL3D/frTumxTRxo/TM65cko89p6v327dvx/nnn48JEyZg48aNmDJlCpYvX467774b3//+9yGEwO23347i4mLcf//9cLlcfR030ZAkAbB7HUj1NgOShHpTNlplW6LDIqJ+pgJolFJQIWXCJZkgCQG7aEUhHLBxTrmE0nXF6aKLLsLPf/5zvPjii8jOzo5a7oMPPsAjjzyChx56CLfddluvBUmUTCQAWd4GCEhoNaaizpQDCBUpSmuiQyOiftAKMxqkdCiS74qUTbiQiTbfnHKUcLoSp71798JsNscsN23aNEybNg1ut7vHgRElMwlAtrcekIBWQyrqzHmAuwYpalOiQyOiPuKFjAY5HW2S75YXg1CQJVqQAnf0j9ao3+lKnPQkTT0pP5CVl5ejvLwcSi/cs0AUDwlAtqceQKfkSVWRorQkNjAi6lUCQItkRYOU5kuQhEAG2pAhWjnY4gDUrafqPvroI6xfvx7V1dVQw25gfPjhh3slsIGirKwMZWVlcDqdsNvtiQ6HkkwweRICrcY01FkLgPYq2BRnokMjol6gQkKdbEebbAUAmIUHOaoTJolD4Q5UcSdO9957L375y19i9OjRKCgoCHm6Swq/o5+IekwCkO2uBYBg8pQlBFI9/NiOaDBzwYRagx2KZASEQKZoRrpo89/4zWtNA1XcidPvf/97/PWvf8XixYv7IBwi0hJIniShosWUgQZbIYQkI629IdGhEVGcBACnnAaHIQOQJBiFF7mqE2ZwmpTBIO7ESZZlnHXWWX0RCxF1QQKQ5aqBJFQ0mzPRaCuAKiRkuOoTHRoR6eSFAfXGbLj8H82lqG3IVp3BsfRo4Iv7WuDPfvYzlJeX90UsRBSDBCDTXYf09joAgDMlH422PE4MTDQItMlWVFoK4ZKtkISKbMWBHNU/XQoNGnFfcbrpppvwrW99C8cffzzGjh0Lk8kUsv2ll17qteCIKJIEwO6qgyxUOGx5aLbmQBYqMvzJFBENPK2yzTcmmyTBrLqR462HibcxDUpxJ07Lli3De++9hxkzZiAnJ4c3hBMlSLq7AZJQ0ZhSAKctD5JQke7iPU9EA03npClFaUGO0ui/AZxzzA1GcSdOTz31FP71r3/hW9/6Vl/EQ0RxSHM1QJVkOG15cKQUQBYqUtt5zxPRQNFm6JQ0eVuQ7a2Pa15AGnjiPnvZ2dk4/vjj+yKWAam8vBxjx47FlClTEh0Kkab09jqk+T+ma0gpRIuZ440RDQTtshW15vzglaZsbz3nmBsC4k6c7rzzTqxYsQKtrckxb1ZZWRl2796NrVu3JjoUIk0SAHtbDVLbGwBJQkPaMLQyeSJKKJdsQa3FlzTZlFZke5g0DRVxf1T3hz/8AV999RUKCgowcuTIiJvDP/nkk14Ljoj0kQBktlUBkoQWSybq04YhS6hIcTUmOjSipOOWzai1FUFIMixKG3I8dUyahpC4E6eLLrqoD8Igop6SAGS2VgJCRYs1Gw3pJRCQkNrGp+2I+otHNqHGWgRVMsCstCPXVc17moaYuBOnFStW9EUcRNQLJACZLRUAgBZrNhozSgCAyRNRP/BKJtSmDocqG2FSXMhzVXOMpiGoT9JgIfhGIUqUQPKU2uqb364xowRNKXmJDYpoiPNKRtSkDYcim2BU3MhrOwYZaqLDoj6gK3EaM2YMnnvuObjd7i7L7d27F9dddx3uv//+XgmOiLpHAmBvOYa01moAgDN9GBxpRfzfl6gPqJBRlzosmDTltx+DgUnTkKXro7ry8nL84he/QFlZGWbPno3JkyejuLgYVqsVDQ0N2L17NzZt2oTdu3dj6dKl+MlPftLXcRNRDBKAjJZKyIoXzvRiNKcWQJFNyGw8AIkpFFGvEADq04rhMVggq17kthyBQebP11AmiTg+V9u8eTNefPFFbNy4EQcOHEBbWxtyc3Nx2mmnYc6cOfj+97+PzMzMPgw3cZxOJ+x2OxwOBzIyMvpsP0cbWrrcLnfxaEZgEHc5uCxFrRO+KryM1ojwEWUitkfWCV8lhb/dhOpfr/HfWfi6KHU168da1lqnRolFs67OfnS137Btwbpqp/Xd2U9EvxQ0yymoN2QCkgSL2o4cdx0MqiekTEhzihIZi0Y53+60ywqNshH1w+soeuqExRoeY3frRKsb7/bO+9HqT18K9KuXb0aWDDpGuJajlAlbr3mjdHjdrupEKRuMMVBWIx5JDi3TsRy2vnN/w8pIndoXABoMdjRLNkhCIA9OWKAAUlgfg7+cfetF+HZfw/EtR2snStmIWOKpE2t/cTDnDu9xG7H09d/ruG4OP/PMM3HmmWf2ehADWXl5OcrLy6H09y9Aol6UprbCoHpQa8qFS7ai2pyPXFcVTMKb6NCIBq0mKQXNkg0AkI1mX9JEQx6fkYyBA2DSUGFVXch3V8MgvPDKJlRbi9AuWxMdFtGg1CpZ0GjwXc3IVFuQAk+MGjRUMHEiSiJm4UGBqwpm1QVVMqDGUoBmYzrveCKKQ7tkRq2cCQBIE21IR1tiA6J+xcSJKMkYoCLPXYMUb7NvihZzDhoseRAc25goJrdkRI3RN2mvTW1HlmjhT06SYeJElIRkCGS7a2F31wNCoMWUgSrbMHiluMfEJUoaXsmAGmOubyoV4UaO2sikKQkxcSJKUhKADK8Tea4qyEKBx2BBZcpwtBpSEx0a0YDjhQHVpnwokgEm1YM8pYF/QJNUt/69VFUV+/btQ3V1NdSwx3rPOeecXgmMiPqHVW1HQeth1FkL4DbYUGcrRKpsRWZ7Df+bJgKgwIAacx4U2Qij8CLfWwu5q7FhaEiLO3HasmULrrjiChw8eDBiahVJkvjYPtEgZBQK8tuOwWHORpM5Cy2WTLiNVmQ3H4NJ7XrGAKKhzCsZUGMphFc2wSC8yPfU+kcF1zG2FQ1JcSdOP/7xjzF58mS8/vrrKCoq0hwokYgGHwlAprseFqUd9ZY8eAxWVGeMhL21CqluB68+UdLxSkZUWwqhyCYYVC/yPdUwSnwGNdnFnTjt3bsX//znP3HCCSf0RTwDDgfApGRjU1pR0HwQ9bZCuEypaEwtgtuYgqzWCiZPlDS8khHVtmG+j+dUD/LcNTBqjQpOSSfud8AZZ5yBffv29UUsAxIHwKRkZBAKcluPwt5aDQiBVosdtWklUPlHg5KAIhlQYyvyJ01u5LdX+pImIui84rRz587g98uWLcONN96IyspKjB8/HiaTKaTsqaee2rsRElFCSADSXfUwKS7UpRXDZUpFdcZxyG06CKPKUZJpaFIho8ZaBK9shkH1IM9VBQMU8J4mCtCVOH3jG9+AJEkhN4P/8Ic/DH4f2Mabw4mGHqu3BXlNh1CXNhxeowXV9uOQ03QIFi9HS6ahRUBCXWoxPAYLZFVBXlsFjFIXE0JTUtKVOO3fv7+v4yCiAcysuJDvPIDatOHwGG2oyRiJrOZjSHU7Eh0aUa8QAOpTCuEypkASKvLaj8EkPIDEK00UStcNC6WlpcGvgwcPYtiwYSHrSktLMWzYMBw8eDCunW/cuBELFixAcXExJEnCyy+/HLJ98eLFkCQp5Gvq1KkhZVwuF5YtW4bc3FykpqZi4cKFOHLkSEiZhoYGLFq0CHa7HXa7HYsWLUJjY2NcsRIlO4NQkOfYD5vLCUgyGtKHw2HL5zx3NOgJAI0phWgzpQNCRW57JcwchoOiiPtOzxkzZqC+vj5ivcPhwIwZM+Jqq6WlBRMmTMCqVauilpk7dy4qKiqCX2+88UbI9uXLl2PNmjV44YUXsGnTJjQ3N2P+/PkhHxleccUV2LFjB9auXYu1a9dix44dWLRoUVyxEpF/qpbmw0hvrQEANKXkoT6jFCrHUKZBzGnLQ4slExACOa2VsCr8GJqii3s4gsC9TOHq6uqQmhrfVA3z5s3DvHnzuixjsVhQWFiouc3hcOCJJ57A008/jVmzZgEAnnnmGZSUlODtt9/GnDlz8Pnnn2Pt2rXYsmULzjjjDADAn//8Z0ybNg179uzB6NGj44qZKNlJAOxt1TAqLjSkFaPdYkdNlhk5jfthVPhfOg0uTlsumqw5AICstirYvM2AgR/PUXS6E6dLLrkEgO9G8MWLF8NisQS3KYqCnTt34swzz+z1ANevX4/8/HxkZmZi+vTpuOeee5Cfnw8A2LZtGzweD2bPnh0sX1xcjHHjxmHz5s2YM2cOPvjgA9jt9mDSBABTp06F3W7H5s2boyZOLpcLLpcruOx0Onu9b0SDWarbAaPTjbq0EniNNlRnn4RsxwFY3c2JDo1Il2ZrNpwpBQAAe2s1Uj38PU+x6U6c7HY7AN8Vp/T0dNhstuA2s9mMqVOnYsmSJb0a3Lx58/Dd734XpaWl2L9/P+644w7MnDkT27Ztg8ViQWVlJcxmM7KyskLqFRQUoLKyEgBQWVkZTLQ6y8/PD5bRsnLlStx111292h+iocbibUN+417UZYyEx5SCuszjkdF8DGn+j/KIBqoWazYaU4sAAOlttUh31QMyP3Km2HQnTk8++SQAYOTIkbjpppvi/liuOy677LLg9+PGjcPkyZNRWlqK119/PXgFTEv4x4laHy1G+8gx4NZbb8UNN9wQXHY6nSgpKYm3C0RDnkH1Iq/xKzSmFaPVlgNn+jC4TanIbNwPWfBRbhp4WqxZaEwfDgBIa6tFRnttgiOiwSTue5xWrFgBAKiursaePXsgSRJOOukkzas6va2oqAilpaXYu3cvAKCwsBButxsNDQ0hV52qq6uDHxsWFhaiqqoqoq2amhoUFBRE3ZfFYgn5OJKIopMgkOk8DJOnFY70YWi3ZqImdwyyG76Cydue6PCIglqs2WhMCyRNdbC3VkGSeU8T6Rd34uR0OlFWVoYXXngh+OSawWDAZZddhvLy8uBHen2hrq4Ohw8fRlGR7/LqpEmTYDKZsG7dOlx66aUAgIqKCnz66ad44IEHAADTpk2Dw+HARx99hNNPPx0A8OGHH8LhcPTJPVk91XmQUS2KxubglTP/NuFflPwrtOrEuiCtdTFOjnKBrmO1iFou/OpeRxnZvz0yIjmsofCYpE7HKvhd2BUOKfyKh9YVkLBjLrpqI0Z74XWhRu4vWkya596/LbxOcD+a/YkSf3gsXe0vvK6O/UBVYAdg9dah1pgFxWhFbe4Y2L0OpCnNvjZVJaJOSJNdDaAbHn94WyGhhe9H7Xq75v5ilNE4t72hy2PQC6Se3Pis56OsKElIl8lJWLsRZbXq+utE9CdQ1v8qdWq72ZCGRoPvb1Sa2oIsswLJnBecf04Efg91/n0U+MXjb0cNL6NRJ6KdiF9ekcdRhK+Ldqy7mvYo2jYdUyVF7L+b7SSDuI/Cj370I3z44Yd47bXX0NjYCIfDgddeew0ff/xx3Pc4NTc3Y8eOHdixYwcA30CbO3bswKFDh9Dc3IybbroJH3zwAQ4cOID169djwYIFyM3NxcUXXwzAd9/V1VdfjRtvvBHvvPMOtm/fju9///sYP3588Cm7MWPGYO7cuViyZAm2bNmCLVu2YMmSJZg/fz6fqCPqAxbhRqGnBlalDUKS0WjKQo05Dx4p7v/TiHqNU05DvT9pSldbkKU2cdJq6hZJxLrEESY1NRX/+c9/cPbZZ4esf//99zF37ly0tLTobmv9+vWaYz9dddVVeOyxx3DRRRdh+/btaGxsRFFREWbMmIG777475F6j9vZ2/PznP8dzzz2HtrY2nHfeefjjH/8YUqa+vh7XX389Xn31VQDAwoULsWrVKmRmZuqO1el0wm63w+FwICMjQ3e9eB2pj/+JpGhXc7q4hUvHFafIyrGvOEUvF/2Kk3YbvjL6rzgF9cIVpy7biNVeD644dXkFKGqdvrniFFFX5xWnkEVVRbMhDQ6jHUKSIQkVGR4H0r2OjvPNK04ReMWpd684QZbRaMhAkyEdgC9pygwkTWFXhvRccYq8msQrTvGwpGf2uI1Y+vrvddyJ04gRI/D6669j/PjxIet37tyJCy64IGLU7qGCiRMTJ836TJw6hCdB/v15JQPqjdlwGawAAJPqRpa7FhbVzcRJAxOn3kucVEhoMGWj1ZACAMhUnMgQrZ13FPLKxCkGJk4AuvFR3S9/+UvccMMNqKioCK6rrKzEz3/+c9xxxx29GhwRDX5GoSDPU4Nsdx1kocAjm1FtKUKDKZsjjlOfUSCjxpznS5qEQI63PjRpIuqmuG86eOyxx7Bv3z6UlpZixIgRAIBDhw7BYrGgpqYGjz/+eLDsJ5980nuREtGgJQFIVVthbWtBozkbrcY0NJsy0GpIhd1dh1Qv7zeh3uOWzKi15EORjJCFilxvHazCHfVqGFE84k6cLrrooj4IY+AqLy9HeXl5yNx3RNQ9BqjIcdci1duMBnM2vLIZDdZ8NKmZsLvqYVWcTKCoR1oNKag350JIMoyqB3neepjgTXRYNITEfY9TsuI9TrzHSbM+73HqEOUep2hlBHyPhzvNWVAl35UAk7cddlctLN7WyPcE73HqFUP1HicBwGnOgtOcDQCwKO3I9dTC0PmXTXg7vMcp+v672U4sSXmPEwA0NjbiL3/5C2699VbU19cD8H0sd/To0V4NjoiGLglAuseBopaDyHDXQxIqPEYralOHoyZ1OFwGW8w2iABAhYQ6a0EwaUrzOJDnqYEMXheg3hf3R3U7d+7ErFmzYLfbceDAASxZsgTZ2dlYs2YNDh48iKeeeqov4iSiIUqGgN3dgDS3A06THc3mTLiNKahJS4HF24KM9npYlLZEh0kDlCIZUGsthNtgBYRAtrsOqUoz72eiPhP3FacbbrgBixcvxt69e2G1WoPr582bh40bN/ZqcESUPAxQkdlei6KmA0h1NQJCwGVMRU1aCWpSh6PdmMLrBxTCI5tRZRsGt8EKWSjIbzvmS5qI+lDcV5y2bt0a8uRcwLBhw1BZWdkrQRFR8jIIL7Laq5HuqkeTJRstZjtcxhS40kfA6mlGZmsVjKon0WFSgrllC2rShkNIBhhVD3LbKmASHl5poj4X9xUnq9UKp9MZsX7Pnj3Iy8vrlaCIiIz+BKqwaX/wClS7KQ2VGaPgsOZC8Pm7pOWRLahN9SVNZqUd+a1HfEkTUT+IO3G68MIL8etf/xoej+9NKkkSDh06hFtuuQXf/va3ez1AIkpugQSqwLkfFk8LIMlosuWiMvMEtJnSEh0e9TOPbEZNeglU2QCTtw15bcdgQN884UikJe7E6cEHH0RNTQ3y8/PR1taG6dOn44QTTkB6ejruueeevoiRiAgm1Y3c5sPIbj4Kg+qBYjCjLqMU9WnDgsMZ0NDmNlhQkz4CqmyEyduOvJajfHKO+l3c9zhlZGRg06ZNePfdd/HJJ59AVVVMnDgRs2bN6ov4Eo4DYBINHBKAFE8TrI5mOK05aLbmoNWSiXZTKrKaK2DzNCU6ROoDAkCLORONKfmAJPuSptYjkKECYNJM/YsDYOrEATA5AKZmfQ6A2SHOATA168Q5ya/LaENDajG8Rt8TvjaXA5ktFZAVd5d1OQBmNyVgAEyvbEJjahHazekAAKu7CdmtFR2/57qY5NfXlhy5rWNHIa8cADMGDoAJIM4rTqqqYvXq1XjppZdw4MABSJKEUaNG4Tvf+Q4WLVqk+ceWiKivWLxtKHB8DWdKHpqsuWiz2OEypcLefBQ2l4O3jw9iAkCzNQfOlHzfH3UhYG+rQZqr3ndee+GPOFF36H7nCSGwcOFC/OhHP8LRo0cxfvx4nHLKKTh48CAWL16Miy++uC/jJCLSJEHA3lqNfMfXMHrbocpGNGSUoi5jJLwGc6LDo25wGVNQZT8ejtRCCEmG2dOKAud+pAeSJqIE0n3FafXq1di4cSPeeecdzJgxI2Tbu+++i4suughPPfUUrrzyyl4PkogoFrPS7rv6ZMtFky0XLksGqsxpSGurQ3pzJWTB+xQHOkU2wpFahDZrFgBAVr2wt1YhxdPEhIkGDN1XnJ5//nncdtttEUkTAMycORO33HILnn322V4NjogoHhJ8H+fkN+yFxe0EJBnNKXmozB2DppR8qPzzOyCpkgxnSgGqsk/2JU1CILW9HoWN+5DqauRZowFFd+K0c+dOzJ07N+r2efPm4X//+1+vBEVE1BMmxYUcxwHkOPbD6G2DkI1wphejKncsmm0cPHOgEJDQlJKHquyT0ZRa4P9YrgX5jq+R1VLBq4Q0IOn+qK6+vh4FBQVRtxcUFKChoaFXgiIi6ikJviewLO4mtJntcKYVQjFY4MgYjubUfKQ1VSClrQ4SxwHqdwJAqy0bTamFUPz3oRm97choqYTV7YTMaVNoANOdOCmKAqMxenGDwQCv19srQRER9RYJQEp7A2ztjWixZaMptQCKwQxHZima0ouQ1lKNlNYayFpDLlCvEgBarf5zYLQAAAyKG+mtVUhpb+B1QBoUdCdOQggsXrwYFotFc7vL5eq1oAaS/h4Ac3i29hQSh/zjO2l9tjosK7XLNo82tESsi90b33/h0cZuAqKPE6UVY9SyWmM+xYih6zGmwsaDCXtkWdbYYazxoYTGH9SIcZWChWOPlSRijc2ksU6El1GjjLekVTfGOFVdxRKx3y7aiWijq3GOosSk9TbRf6y72J9/TCY7gHSlBs0iFU2GdCgGM5wZw9GUXow0pRnp3qaO6Tv0jPUUZT+DXneu+EQdv0mGCgnNcgqaDGlQJN+fHVkoyFBbkC5aIdlMgC2/c6XwRiLajRh3KHzMpMC4S1p1orxqjuMUbF/SLitr1IkSW9SYtfYTsb4HYzHFMXSD6MawQsk2GqTuATB/8IMf6GrwySef7FFAA1V/DYAZTVeJU7RkK0ArcdKr9xIn7cJdJ07aZXo2OGf8iZOuRCNa2e4Matlle/EnTrEG+NQVSz8lTlp6M3EKKQ6gWbKhyZAOr2wKxpWqtCBNaYZZaY/eXhz7GZR6KXFySya0GNPRItuCf9RloSBDaUaa2gJZ7yCPPUicNOswceqijb5NnFJs1rjbj9eAGQBzqCZERJScJABpSgtSlRa0y1Y4jRlwyxa0GNPQYkyDWWlHqtICm7cVBh3XaMlHAGiVU9BsTINb7viEwqR6kKY2I1VtjX+SVKIBJO656oiIhhIJgE1th83dDpdkRpMxHW2yDW6DFW6DFQ2mbJhVF2xKG2xKK4zCw3txwggALtmCVjkFbQZbx6TLQiBFbUOa2gKLcPO40ZDAxImIyM8i3LB46qBARotsQ5shFW6DJZhEOZAFg+qFVWmFRXXDpLpgStJEyisZ0C7b0G7wfXX+mMggvEj1+j7yNPLyEg0xTJyIiMIYoCLD60SG1+lLEAy+Kyku2QpFNqJFzkDwzkEhYBIemFQ3jIobJuGGUfXAqHogD5GhDgQARTLCLZvhkq1oN1jhlUOns5GFApvShhS1FRbV1SmZZOZEQwsTJyKiLhiFgjRvE9K8TVAhBRMHj2yBWzZDSDI8khke2RzxG1UWCgyqB0bhhVH1wiA8MKpe33qhQBbKgEuuVEjwSma4ZV+fPLI52M8QQsCsumBV2mBT2mCSlKS88kbJh4kTEZFOMgRsahtsahuAjisxHtkEj2SCRzLCK5vhlU1QJYPvy2CAp4s2JaH6EiihdiRTQoEMFZIQkIUKCSpkoYZ9r0DSMQa6ACAgQ5HkYExKILZO3yuB7+UofxaEgEl1w6y6YVXbYFXaIaPTk4wSB62k5MDEiYiomyTAdzVJ8cKGtpDhCFTI8MpGeCUjFNnkfzVCkYz+JMUISBKE5EtqFKDLBEuTECFJVSDZUiXJnyz5kqKoj7hHIQsFJtUdTJRMqjtp7+UiCsfEKYb+HgCTiIYGGSrMqhtmuDVHnPVdCZI0rgD5rgwJyFAlGUKSoMLgS4YkGSrkjmRIkiDgv1oUI56OK1uBjwnVkI8MDf4vo+rRHCiWiHyYOMVQVlaGsrKy4IBaRES9QQIgQUAWXhhF/NNV+ZKuQCIVSLJkCMi+JAkdiZLvalQ891LxYzeiaJg4ERENQhKE7yoRFAyw+8uJhjQ+J0pERESkExMnIiIiIp2YOBERERHpxMSJiIiISCcmTkREREQ6MXEiIiIi0omJUwzl5eUYO3YspkyZkuhQiIiIKMGYOMVQVlaG3bt3Y+vWrYkOhYiIiBKMiRMRERGRTkyciIiIiHRi4kRERESkExMnIiIiIp2YOBERERHpZEx0AKTPiOy0btdVROTU6RqrQkhSoG7s9rvKvqVAQ2HTt8tSZNnwfUdvX4S1rRFT2KaOxcgORZQNa1eWNHroXxcegdxFx8LDlcJPglAj64SvCyyH1RXx1NVYFtHKBKix24+IQetNFqX9iFi7ajdGW93dj1YfowSkr9xQ0sXPWlSy72ck/GgJrZ+n4H7krpc1tkW0F1gOxKzRRrCOLIfVkSPq6G1fs19h60TYcdR6K6kx3l9dbVWjbdT6UYy6n8j1UdvthpTeaypheMWJiIiISCcmTjFwAEwiIiIKYOIUAwfAJCIiogAmTkREREQ6MXEiIiIi0omJExEREZFOTJyIiIiIdGLiRERERKQTEyciIiIinZg4EREREenExCkGDoBJREREAUycYuAAmERERBTAxImIiIhIJyZORERERDoxcSIiIiLSiYkTERERkU5MnIiIiIh0YuJEREREpBMTJyIiIiKdmDjFwAEwiYiIKICJUwwcAJOIiIgCmDgRERER6cTEiYiIiEgnJk5EREREOjFxIiIiItLJmOgAqO8pqohYF7lGu4Ckp/2wZVnqVEto70nS0XC0rF4KVo5sW47Sblf7i9yPCNuPRp1o++niyIbXCW9fljR67F8XLCkFyobW1QpVCj/2Qg3bHrrcVVmt8xhcE6tdHfuJaL2rOgFq/P3p8n0f5b0ara0hTeu9qLtutB+O2G2KaGXk6D8bsZZD2gx8H4gxvExwOfrPfvjbRPWvECHrAoUDdUJXaPxKjvre1CrbEYu+96yed6zet79m+z2oOxjxihMRERGRTkyciIiIiHRi4hQDB8AkIiKiACZOMXAATCIiIgpg4kRERESkExMnIiIiIp2YOBERERHplNDEaeXKlZgyZQrS09ORn5+Piy66CHv27AkpI4TAnXfeieLiYthsNpx77rn47LPPQsq4XC4sW7YMubm5SE1NxcKFC3HkyJGQMg0NDVi0aBHsdjvsdjsWLVqExsbGvu4iERERDSEJTZw2bNiAsrIybNmyBevWrYPX68Xs2bPR0tISLPPAAw/g4YcfxqpVq7B161YUFhbi/PPPR1NTU7DM8uXLsWbNGrzwwgvYtGkTmpubMX/+fChKx9CMV1xxBXbs2IG1a9di7dq12LFjBxYtWtSv/SUiIqLBTRJ6hx7tBzU1NcjPz8eGDRtwzjnnQAiB4uJiLF++HL/4xS8A+K4uFRQU4P7778e1114Lh8OBvLw8PP3007jssssAAMeOHUNJSQneeOMNzJkzB59//jnGjh2LLVu24IwzzgAAbNmyBdOmTcMXX3yB0aNHx4zN6XTCbrfD4XAgIyOj7w5CH9hf2xSxTu9J1zNyeLjwUa012+2VkcO19h3//np1P9F3o2Pk8Oh1wzcNhJHDdbcbz3701AnoxsjhXeLI4R04cnj0GOMZOTxYJ3RFMo8cPiwrtfs70qmv/14PqHucHA4HACA7OxsAsH//flRWVmL27NnBMhaLBdOnT8fmzZsBANu2bYPH4wkpU1xcjHHjxgXLfPDBB7Db7cGkCQCmTp0Ku90eLBPO5XLB6XSGfBEREVFyGzCJkxACN9xwA84++2yMGzcOAFBZWQkAKCgoCClbUFAQ3FZZWQmz2YysrKwuy+Tn50fsMz8/P1gm3MqVK4P3Q9ntdpSUlPSsg0RERDToDZjEaenSpdi5cyeef/75iG3hH2kIIbr8GEWrjFb5rtq59dZb4XA4gl+HDx/W0w0iIiIawgZE4rRs2TK8+uqreO+99zB8+PDg+sLCQgCIuCpUXV0dvApVWFgIt9uNhoaGLstUVVVF7LempibialaAxWJBRkZGyBcRERElt4QmTkIILF26FC+99BLeffddjBo1KmT7qFGjUFhYiHXr1gXXud1ubNiwAWeeeSYAYNKkSTCZTCFlKioq8OmnnwbLTJs2DQ6HAx999FGwzIcffgiHwxEsQ0RERBSLMZE7Lysrw3PPPYdXXnkF6enpwStLdrsdNpsNkiRh+fLluPfee3HiiSfixBNPxL333ouUlBRcccUVwbJXX301brzxRuTk5CA7Oxs33XQTxo8fj1mzZgEAxowZg7lz52LJkiV4/PHHAQDXXHMN5s+fr+uJOiIiIiIgwYnTY489BgA499xzQ9Y/+eSTWLx4MQDg5ptvRltbG37yk5+goaEBZ5xxBt566y2kp6cHyz/yyCMwGo249NJL0dbWhvPOOw+rV6+GwWAIlnn22Wdx/fXXB5++W7hwIVatWtW3HSQiIqIhZUCN4zSQcRwn/TiOk746HMepi20cxymxOI5T9Bg5jlNk+xzHiYiIiIi0MHEiIiIi0omJUwzl5eUYO3YspkyZkuhQiIiIKMGYOMVQVlaG3bt3Y+vWrYkOhYiIiBKMiRMRERGRTkyciIiIiHRi4kRERESkExMnIiIiIp0SOnI49Y9Ruekxy+yriRwks7uEnuE1w4poDTenRKna1UCRwfZ6ZYDN6P3ozkCbkfsNbV9rwM3oA22KmOUi25NDy2oMNihFfBNoP3rHwjfFGngztKyegSljlNEYgDBW3V7ZbzLqxsCYUQe1jKddzRFeQ+tE7EdjMEub1Ro7lhjqm1p97fqXNQez9L8nw99BgR+NruqE643BK1Udo1vGM/5lPANeDkW84kRERESkExMnIiIiIp2YOMXAATCJiIgogIlTDBwAk4iIiAKYOBERERHpxMSJiIiISCcmTkREREQ6MXEiIiIi0omJExEREZFOTJyIiIiIdGLiRERERKQTE6cYOAAmERERBTBxioEDYBIREVEAEyciIiIinZg4EREREenExImIiIhIJyZORERERDoxcSIiIiLSiYkTERERkU5MnIiIiIh0YuIUAwfAJCIiogAmTjFwAEwiIiIKYOJEREREpBMTJyIiIiKdmDgRERER6cTEiYiIiEgnJk5EREREOjFxIiIiItLJmOgAaGAQovfaUrrRliTF075vB3JXdaLE0LmKEqOyrCOoeOKO9V9KV21F62uXdaIcBKmLStGPqa+trrobXrdjP5E9D5aVtI+KFHVBa7/RC0TbJOl5wws1dhnqEOVchhNx/NDoOU1qWKHAkqoG2ug4jzbde44uOz1Fc/2h+uaOGESM2ERouc7UsJ9bfceg6+262oj2SzPOdpIBrzjFwAEwiYiIKICJUwwcAJOIiIgCmDgRERER6cTEiYiIiEgnJk5EREREOjFxIiIiItKJiRMRERGRTkyciIiIiHRi4kRERESkExOnGDgAJhEREQUwcYqBA2ASERFRABMnIiIiIp2YOBERERHpxMSJiIiISCcmTkREREQ6MXEiIiIi0omJExEREZFOTJyIiIiIdGLiFAMHwCQiIqIAJk4xcABMIiIiCmDiRERERKQTEyciIiIinZg4EREREenExImIiIhIJyZORERERDoxcSIiIiLSiYkTERERkU5MnGLgAJhEREQUIAkhRKKDGAycTifsdjscDgcyMjISHU6v21PtTHQIcZMh9Uo7UjeakXXWiadpWUcgemKN9d+QpKMRvf3ztRd/HHpiiDeWrorpaSOemKK1qUb5bRrP8Rxqoh0TPbrz5yl8f50Xc9JTuh9MDLsqHB37DItBDVuhqJHrFZ19DW9Lq92etKMnjp6c01kn5nW/sk59/feaV5yIiIiIdGLiRERERKQTEyciIiIinZg4EREREenExImIiIhIJyZORERERDoxcSIiIiLSiYlTDBwAk4iIiAKYOMVQVlaG3bt3Y+vWrYkOhYiIiBKMiRMRERGRTkyciIiIiHRi4kRERESkU0ITp5UrV2LKlClIT09Hfn4+LrroIuzZsyekzOLFiyFJUsjX1KlTQ8q4XC4sW7YMubm5SE1NxcKFC3HkyJGQMg0NDVi0aBHsdjvsdjsWLVqExsbGvu4iERERDSEJTZw2bNiAsrIybNmyBevWrYPX68Xs2bPR0tISUm7u3LmoqKgIfr3xxhsh25cvX441a9bghRdewKZNm9Dc3Iz58+dDUZRgmSuuuAI7duzA2rVrsXbtWuzYsQOLFi3ql34SERHR0GBM5M7Xrl0bsvzkk08iPz8f27ZtwznnnBNcb7FYUFhYqNmGw+HAE088gaeffhqzZs0CADzzzDMoKSnB22+/jTlz5uDzzz/H2rVrsWXLFpxxxhkAgD//+c+YNm0a9uzZg9GjR/dRD4mIiGgoGVD3ODkcDgBAdnZ2yPr169cjPz8fJ510EpYsWYLq6urgtm3btsHj8WD27NnBdcXFxRg3bhw2b94MAPjggw9gt9uDSRMATJ06FXa7PVgmnMvlgtPpDPkiIiKi5DZgEichBG644QacffbZGDduXHD9vHnz8Oyzz+Ldd9/FQw89hK1bt2LmzJlwuVwAgMrKSpjNZmRlZYW0V1BQgMrKymCZ/Pz8iH3m5+cHy4RbuXJl8H4ou92OkpKS3uoqERERDVIJ/aius6VLl2Lnzp3YtGlTyPrLLrss+P24ceMwefJklJaW4vXXX8cll1wStT0hBCRJCi53/j5amc5uvfVW3HDDDcFlp9PJ5ImIiCjJDYgrTsuWLcOrr76K9957D8OHD++ybFFREUpLS7F3714AQGFhIdxuNxoaGkLKVVdXo6CgIFimqqoqoq2amppgmXAWiwUZGRkhX0RERJTcEpo4CSGwdOlSvPTSS3j33XcxatSomHXq6upw+PBhFBUVAQAmTZoEk8mEdevWBctUVFTg008/xZlnngkAmDZtGhwOBz766KNgmQ8//BAOhyNYhoiIiCiWhH5UV1ZWhueeew6vvPIK0tPTg/cb2e122Gw2NDc3484778S3v/1tFBUV4cCBA7jtttuQm5uLiy++OFj26quvxo033oicnBxkZ2fjpptuwvjx44NP2Y0ZMwZz587FkiVL8PjjjwMArrnmGsyfP59P1BEREZFuCU2cHnvsMQDAueeeG7L+ySefxOLFi2EwGLBr1y489dRTaGxsRFFREWbMmIEXX3wR6enpwfKPPPIIjEYjLr30UrS1teG8887D6tWrYTAYgmWeffZZXH/99cGn7xYuXIhVq1b1fSeJiIhoyJCEECLRQQwGTqcTdrsdDodjSN7vtKd68A23IEP7xv54RXk+oOt966wTT9OyjkD0xBrr8/doD0SExhJ7Px3txR+HnhjijaWrYnraiCemaG2qUX6bxnM8h5pox0SP7vx5Ct9f58Wc9JTuBxPDrgpHxz7DYlDDVihq5HpFZ1/D29Jqtyft6ImjJ+d01ol53a+sU1//vR4QN4cTERERDQZMnIiIiIh0YuIUQ3l5OcaOHYspU6YkOhQiIiJKMCZOMZSVlWH37t3YunVrokMhIiKiBGPiRERERKQTEyciIiIinQbMXHWUWKPztR/Z7GqYgmhPrXZ+qrsvB7tQ0LPGg3F2oxnFX6c7QyJEe+q9q8eA5Rixdm5SibH/njwa3xvDIehtJ2r7uoYWiNGG5lrtgxvvMAX9aagOc9CTx90T5YijPfi9EiX+rocS6M7QC/rqRIsn3rbUwXhi+gCvOBERERHpxMSJiIiISCcmTkREREQ6MXEiIiIi0omJUwwcAJOIiIgCmDjFwAEwiYiIKICJExEREZFOTJyIiIiIdGLiRERERKQTEyciIiIinZg4EREREenExImIiIhIJyZORERERDoxcYqBA2ASERFRABOnGDgAJhEREQUwcSIiIiLSiYkTERERkU5MnIiIiIh0YuJEREREpBMTJyIiIiKdmDgRERER6cTEiYiIiEgnJk4xcABMIiIiCmDiFAMHwCQiIqIAJk5EREREOjFxIiIiItKJiRMRERGRTkyciIiIiHRi4kRERESkExMnIiIiIp2MiQ6ABrbR+Rk9qv9FlbPL7aroUfM90419y1LostKNRqTYRSIIob+WFKOoImLHHN7PjkC62G+g/Zitd9F+N8TqL9C9/xClYMP98ybtjWOi51j0lfBjrCZov/n21H7ac6jPa5oj1hn8J9UQ5eQadJwwuYsyqo6fZUDfz3ywbBy/lOMpO5TwilMMHACTiIiIApg4xcABMImIiCiAiRMRERGRTkyciIiIiHRi4kRERESkExMnIiIiIp2YOBERERHpxMSJiIiISCcmTkREREQ6MXGKgQNgEhERUQATpxg4ACYREREFMHEiIiIi0omJExEREZFOTJyIiIiIdGLiRERERKQTEyciIiIinZg4EREREenExImIiIhIJyZOMXAATCIiIgpg4hQDB8AkIiKiACZORERERDoxcSIiIiLSiYkTERERkU5MnIiIiIh0YuJEREREpBMTJyIiIiKdmDgRERER6cTEKQYOgElEREQBTJxi4ACYREREFMDEiYiIiEgnY6IDGCyEEAAAp9OZ4EgGl+amro+XKvopkF4iSz1vQ+pGGzLiqNQLMXann/FUkbtzEHqw3+78hyj1Yox6JOq91VvCj7GaoP1aJaWf9hyqvaUpYp3Bf1LlKCfXoOOEdfWzogp9v0AVneUAQI3jl7LSjV/g/fE3NLAPEUe/48HESae6ujoAQElJSYIjISIiGpxu78d91dXVwW6393q7TJx0ys7OBgAcOnQo7hMxZcqUbt0jNRjqOZ1OlJSU4PDhw8jIyOjz/Q2WeslwXLpbNxmODY9L79XjcdGWiOPSk7r9Wc/hcGDEiBHBv9u9jYmTTrLsuyBst9vjfpMaDIa46wymegCQkZHB46JhKB+XntYdyseGx6V36wE8LtH053HpSd1EHJvA3+3expvD+0FZWdmQrtddg6V/PC59U7c/9zdY6nXXYOkfj0vv1uuuRPzMD5Zjo4ck+uruqSHG6XTCbrfD4XB0O/sdinhctPG4RMdjo43HRRuPizYel+j6+tjwipNOFosFK1asgMViSXQoAwqPizYel+h4bLTxuGjjcdHG4xJdXx8bXnEiIiIi0olXnIiIiIh0YuJEREREpBMTJyIiIiKdmDgRERER6cTESYc//vGPGDVqFKxWKyZNmoT3338/0SH1qzvvvBOSJIV8FRYWBrcLIXDnnXeiuLgYNpsN5557Lj777LMERtx3Nm7ciAULFqC4uBiSJOHll18O2a7nWLhcLixbtgy5ublITU3FwoULceTIkX7sRe+LdVwWL14c8R6aOnVqSJmheFxWrlyJKVOmID09Hfn5+bjooouwZ8+ekDLJ+J7Rc1yS8T3z2GOP4dRTTw0Oajlt2jS8+eabwe3J+F4BYh+X/n6vMHGK4cUXX8Ty5ctx++23Y/v27fjmN7+JefPm4dChQ4kOrV+dcsopqKioCH7t2rUruO2BBx7Aww8/jFWrVmHr1q0oLCzE+eefj6amyEkvB7uWlhZMmDABq1at0tyu51gsX74ca9aswQsvvIBNmzahubkZ8+fPh6IkZnLS3hDruADA3LlzQ95Db7zxRsj2oXhcNmzYgLKyMmzZsgXr1q2D1+vF7Nmz0dLSEiyTjO8ZPccFSL73zPDhw3Hffffh448/xscff4yZM2fiwgsvDCZHyfheAWIfF6Cf3yuCunT66aeLH//4xyHrTj75ZHHLLbckKKL+t2LFCjFhwgTNbaqqisLCQnHfffcF17W3twu73S7+7//+r58iTAwAYs2aNcFlPceisbFRmEwm8cILLwTLHD16VMiyLNauXdtvsfel8OMihBBXXXWVuPDCC6PWSYbjIoQQ1dXVAoDYsGGDEILvmYDw4yIE3zMBWVlZ4i9/+QvfK2ECx0WI/n+v8IpTF9xuN7Zt24bZs2eHrJ89ezY2b96coKgSY+/evSguLsaoUaNw+eWX4+uvvwYA7N+/H5WVlSHHyGKxYPr06Ul3jPQci23btsHj8YSUKS4uxrhx44b88Vq/fj3y8/Nx0kknYcmSJaiurg5uS5bj4nA4AHRMGs73jE/4cQlI5veMoih44YUX0NLSgmnTpvG94hd+XAL6873CSX67UFtbC0VRUFBQELK+oKAAlZWVCYqq/51xxhl46qmncNJJJ6Gqqgq/+c1vcOaZZ+Kzzz4LHgetY3Tw4MFEhJsweo5FZWUlzGYzsrKyIsoM5ffUvHnz8N3vfhelpaXYv38/7rjjDsycORPbtm2DxWJJiuMihMANN9yAs88+G+PGjQPA9wygfVyA5H3P7Nq1C9OmTUN7ezvS0tKwZs0ajB07NvgHPlnfK9GOC9D/7xUmTjpIkhSyLISIWDeUzZs3L/j9+PHjMW3aNBx//PH429/+FrwBL9mPUWfdORZD/Xhddtllwe/HjRuHyZMno7S0FK+//jouueSSqPWG0nFZunQpdu7ciU2bNkVsS+b3TLTjkqzvmdGjR2PHjh1obGzEv/71L1x11VXYsGFDcHuyvleiHZexY8f2+3uFH9V1ITc3FwaDISIjra6ujsj6k0lqairGjx+PvXv3Bp+u4zGCrmNRWFgIt9uNhoaGqGWSQVFREUpLS7F3714AQ/+4LFu2DK+++iree+89DB8+PLg+2d8z0Y6LlmR5z5jNZpxwwgmYPHkyVq5ciQkTJuD3v/990r9Xoh0XLX39XmHi1AWz2YxJkyZh3bp1IevXrVuHM888M0FRJZ7L5cLnn3+OoqIijBo1CoWFhSHHyO12Y8OGDUl3jPQci0mTJsFkMoWUqaiowKeffppUx6uurg6HDx9GUVERgKF7XIQQWLp0KV566SW8++67GDVqVMj2ZH3PxDouWpLlPRNOCAGXy5W075VoAsdFS5+/V+K+nTzJvPDCC8JkMoknnnhC7N69WyxfvlykpqaKAwcOJDq0fnPjjTeK9evXi6+//lps2bJFzJ8/X6SnpwePwX333Sfsdrt46aWXxK5du8T3vvc9UVRUJJxOZ4Ij731NTU1i+/btYvv27QKAePjhh8X27dvFwYMHhRD6jsWPf/xjMXz4cPH222+LTz75RMycOVNMmDBBeL3eRHWrx7o6Lk1NTeLGG28UmzdvFvv37xfvvfeemDZtmhg2bNiQPy7XXXedsNvtYv369aKioiL41draGiyTjO+ZWMclWd8zt956q9i4caPYv3+/2Llzp7jtttuELMvirbfeEkIk53tFiK6PSyLeK0ycdCgvLxelpaXCbDaLiRMnhjwymwwuu+wyUVRUJEwmkyguLhaXXHKJ+Oyzz4LbVVUVK1asEIWFhcJisYhzzjlH7Nq1K4ER95333ntPAIj4uuqqq4QQ+o5FW1ubWLp0qcjOzhY2m03Mnz9fHDp0KAG96T1dHZfW1lYxe/ZskZeXJ0wmkxgxYoS46qqrIvo8FI+L1jEBIJ588slgmWR8z8Q6Lsn6nvnhD38Y/FuTl5cnzjvvvGDSJERyvleE6Pq4JOK9IgkhRPzXqYiIiIiSD+9xIiIiItKJiRMRERGRTkyciIiIiHRi4kRERESkExMnIiIiIp2YOBERERHpxMSJiIiISCcmTkREREQ6MXEion4jSRJefvnlRIfRJ9xuN0444QT897//BQAcOHAAkiRhx44dvbqfVatWYeHChb3aJhHpx8SJiHpk8eLFkCQJkiTBZDKhoKAA559/Pv76179CVdWQshUVFZg3b56udgdbkvWnP/0JpaWlOOuss/p0P0uWLMHWrVuxadOmPt0PEWlj4kREPTZ37lxUVFTgwIEDePPNNzFjxgz89Kc/xfz58+H1eoPlCgsLYbFYEhhp33n00Ufxox/9qM/3Y7FYcMUVV+DRRx/t830RUSQmTkTUYxaLBYWFhRg2bBgmTpyI2267Da+88grefPNNrF69Oliu81Ukt9uNpUuXoqioCFarFSNHjsTKlSsBACNHjgQAXHzxxZAkKbj81Vdf4cILL0RBQQHS0tIwZcoUvP322yGxjBw5Evfeey9++MMfIj09HSNGjMCf/vSnkDJHjhzB5ZdfjuzsbKSmpmLy5Mn48MMPg9v//e9/Y9KkSbBarTjuuONw1113hSSA4T755BPs27cP3/rWt6KWUVUVS5YswUknnYSDBw8Gj8fjjz+O+fPnIyUlBWPGjMEHH3yAffv24dxzz0VqaiqmTZuGr776KqSthQsX4uWXX0ZbW1vU/RFR32DiRER9YubMmZgwYQJeeuklze1/+MMf8Oqrr+Lvf/879uzZg2eeeSaYIG3duhUA8OSTT6KioiK43NzcjAsuuABvv/02tm/fjjlz5mDBggU4dOhQSNsPPfQQJk+ejO3bt+MnP/kJrrvuOnzxxRfBNqZPn45jx47h1Vdfxf/+9z/cfPPNwY8V//Of/+D73/8+rr/+euzevRuPP/44Vq9ejXvuuSdqXzdu3IiTTjoJGRkZmtvdbjcuvfRSfPzxx9i0aRNKS0uD2+6++25ceeWV2LFjB04++WRcccUVuPbaa3Hrrbfi448/BgAsXbo0pL3JkyfD4/Hgo48+ihoTEfURQUTUA1dddZW48MILNbdddtllYsyYMcFlAGLNmjVCCCGWLVsmZs6cKVRV1azbuWxXxo4dKx599NHgcmlpqfj+978fXFZVVeTn54vHHntMCCHE448/LtLT00VdXZ1me9/85jfFvffeG7Lu6aefFkVFRVFj+OlPfypmzpwZsm7//v0CgHj//ffFrFmzxFlnnSUaGxsj+vjLX/4yuPzBBx8IAOKJJ54Irnv++eeF1WqN2GdWVpZYvXp11JiIqG/wihMR9RkhBCRJ0ty2ePFi7NixA6NHj8b111+Pt956K2Z7LS0tuPnmmzF27FhkZmYiLS0NX3zxRcQVp1NPPTX4vSRJKCwsRHV1NQBgx44dOO2005Cdna25j23btuHXv/410tLSgl9LlixBRUUFWltbNeu0tbXBarVqbvve976H5uZmvPXWW7Db7RHbO8daUFAAABg/fnzIuvb2djidzpB6NpstajxE1HeYOBFRn/n8888xatQozW0TJ07E/v37cffdd6OtrQ2XXnopvvOd73TZ3s9//nP861//wj333IP3338fO3bswPjx4+F2u0PKmUymkGVJkoIfxdlsti73oaoq7rrrLuzYsSP4tWvXLuzduzdqcpSbm4uGhgbNbRdccAF27tyJLVu2aG7vHGsgydRaF/6EYn19PfLy8rrsCxH1PmOiAyCioendd9/Frl278LOf/SxqmYyMDFx22WW47LLL8J3vfAdz585FfX09srOzYTKZoChKSPn3338fixcvxsUXXwzAd7/SgQMH4orr1FNPxV/+8pfgfsJNnDgRe/bswQknnKC7zdNOOw2PPfaY5hW26667DuPGjcPChQvx+uuvY/r06XHFq+Wrr75Ce3s7TjvttB63RUTxYeJERD3mcrlQWVkJRVFQVVWFtWvXYuXKlZg/fz6uvPJKzTqPPPIIioqK8I1vfAOyLOMf//gHCgsLkZmZCcD3dNw777yDs846CxaLBVlZWTjhhBPw0ksvYcGCBZAkCXfccUfElZhYvve97+Hee+/FRRddhJUrV6KoqAjbt29HcXExpk2bhl/96leYP38+SkpK8N3vfheyLGPnzp3YtWsXfvOb32i2OWPGDLS0tOCzzz7DuHHjIrYvW7YMiqJg/vz5ePPNN3H22WfHFXO4999/H8cddxyOP/74HrVDRPHjR3VE1GNr165FUVERRo4ciblz5+K9997DH/7wB7zyyiswGAyaddLS0nD//fdj8uTJmDJlCg4cOIA33ngDsuz7tfTQQw9h3bp1KCkpCV5ZeeSRR5CVlYUzzzwTCxYswJw5czBx4sS4YjWbzXjrrbeQn5+PCy64AOPHj8d9990XjHPOnDl47bXXsG7dOkyZMgVTp07Fww8/HPIkXLicnBxccsklePbZZ6OWWb58Oe666y5ccMEF2Lx5c1wxh3v++eexZMmSHrVBRN0jCSFEooMgIhrsdu3ahVmzZmHfvn1IT0/vs/18+umnOO+88/Dll19q3mxORH2LV5yIiHrB+PHj8cADD8R9z1W8jh07hqeeeopJE1GC8IoTERERkU684kRERESkExMnIiIiIp2YOBERERHpxMSJiIiISCcmTkREREQ6MXEiIiIi0omJExEREZFOTJyIiIiIdGLiRERERKQTEyciIiIinZg4EREREenExImIiIhIJyZORERERDoxcSIiIiLSiYkTERERkU5MnIiIiIh0YuJEREREpBMTJyIiIiKdmDgRERER6cTEiYiIiEgnJk5EREREOjFxIiIiItKJiRMRERGRTkyciIiIiHRi4kRERESkExMnIiIiIp2YOBERERHpxMSJiIiISCcmTkREREQ6MXEiIiIi0omJExEREZFOTJyIiIiIdGLiRERERKQTEyciIiIinZg4EREREenExImIiIhIJyZORERERDoxcSIiIiLSiYkTERERkU5MnIiIiIh0YuJEREREpBMTJyIiIiKdmDgRERER6cTEiYiIiEgnJk5EREREOjFxIiIiItKJiRMRERGRTkyciIiIiHRi4kRERESkExMnIiIiIp2YOBERERHpxMSJiIiISCcmTkREREQ6MXEiIiIi0omJExEREZFOTJyIiIiIdGLiRERERKQTEyciIiIinZg4EREREenExImIiIhIJyZORERERDoxcSIiIiLSiYkTERERkU5MnIiIiIh0YuJEREREpBMTJyIiIiKdmDgRERER6WRMdACDSXt7O9xud6LDICIioi6YzWZYrdY+aZuJk07t7e2wpWcD3rZEh0JERERdKCwsxP79+/skeWLipJPb7Qa8bTCOvRSSyXciJNkQfJUk//cG36efcqdtoa+y9rK/Pgza9WTZoFE37DUYQ7R9BdqS/Mv+V6ljufP3ACCHLQfXB5cR0oYcWA5rW5YkIFobUmhb0eILL28Ia8cgS8F1gVejv82Isv6dhZcPvkqRbfvai9xHtLIdr9BcL4fXlyT43wWa23yvodsD6w2B9eHHJmy7JHXEH/icPnwfgXjlsDYjtwfOh38ZHctycH+h2xBsM2Sxo63AMjriDSkvdWwP3xYgS1Hqhu8rsF4I3zdC9S+rwe8RsS10OdqrFLFeo51YddVodaO1rfhXB+orvi8AQlFC2/S/Cv/2QFvBZTW0fGT9TvsKthG+/7CyUdsIrSdU4S+vdqxTwl/9dUSU7f421Ih6akS9wDo12r7U0FcEyodvD9YToe2potMhCiujiJD1auBcKuHxi7B6gX50tBNZRoTE17Ec2EfY+kA9f5tK8DX0+5AyQETZ0NdAm7HKdWwTUdvUXo5oG4AbKp6tPAq3283EaUAwmCAZzADCEic5PGnRmziFvXaZOMVqW18MXSUmEQlS3IlTlPKSRtvdTZyiJihaiVO0svoSJ632otbp7mvnxCksMYpIanqYOMmdEyf/uuiJU5QYepA4BZdjJDPdSZw62opSN8q+4kucuk6YAklB7MRJjdpG/ImTRjITWNaZtEQkUlETJyWsvqq9Tqus3sSpU+ISLXFSleh1Oq+Pmjh1qhc9cep6HzETp05JnCr7kwLZv86/HL5eDSSNkn+7pIYsR6wPJIgQEPCvQ8c6Xx0pdDnwKvzrAwlLMBGJnThFLCMskUK09dGXlWDc8dcNWQ5kXn2IN4cTERER6cTEiYiIiEgnJk5EREREOjFxIiIiItKJiRMRERGRTkyciIiIiHRi4kRERESkExMnIiIiIp2YOBERERHpxMSJiIiISCcmTkREREQ6MXEiIiIi0omJExEREZFOTJyIiIiIdGLiRERERKQTEyciIiIinZg4EREREenExImIiIhIJyZORERERDoxcSIiIiLSiYkTERERkU7GRAcw6CgeCNng+150elUD38v+F/9yxKusvSz5loW/HUkOfYVsgOQvG74t+OpvQzIE1muXhyz5l/2vkuSPWQp+H9gmwpaD64PLCGsD2m1LUsR+5bAygbbkKPGFlw/EoHZ6VcPWBeIJLCv+V4N/Z4bAcvhr2D4NnerrLdvxCs31cnh9SYL/DGlu872Gbg+sNwTWS2H7CNsuSR3xB/5rCt9HIF45rM3I7YHz4V9Gx7Ic3F/oNgTbDFnsaCuwjI54Q8pLHdvDtwXIUpS64fsKrBfC941Q/ctq8HtEbAtdjvYqRazXaCdWXTVa3WhtK/7VgfqK7wuAUJTQNv2vwr890FZwWQ0tH1m/076CbYTvP6xs1DZC6wlV+MurHeuU8Fd/HRFlu78NNaKeGlEvsE6Nti819BWB8uHbg/VEaHuq6HSIwsooImS9GjiXSnj8IqxeoB8d7USWESHxdSz79xFYr4auD8SgBF9Dvw8pA0SUDX0NtBmrXMc2EbVN7eWItgG44T9XfYSJk05msxmFhYWo3P33RIdCREREXSgsLITZbO6TtiURSFkppvb2drjd7oj1TqcTJSUlOHz4MDIyMhIQWeIka9+Ttd8A+56MfU/WfgPJ2/fB3m+z2Qyr1donbfOKUxysVmuXJyIjI2NQvsF6Q7L2PVn7DbDvydj3ZO03kLx9T9Z+d4U3hxMRERHpxMSJiIiISCcmTr3AYrFgxYoVsFgsiQ6l3yVr35O13wD7nox9T9Z+A8nb92Tttx68OZyIiIhIJ15xIiIiItKJiRMRERGRTkyciIiIiHRi4kRERESkU1InTitXrsSUKVOQnp6O/Px8XHTRRdizZ09IGUmSNL9++9vfAgDq6+uxbNkyjB49GikpKRgxYgSuv/56OByOmPv/4x//iFGjRsFqtWLSpEl4//33+6SfWhLZ9zvvvDOizcLCwj7ra2e90W8AuPbaa3H88cfDZrMhLy8PF154Ib744ouY+x/s5xzoXt8Tec6B3ut7gBAC8+bNgyRJePnll2PuP1HnPZH9Hirn/Nxzz43Yfvnll8fc/2A/593pd6LPeX9J6sRpw4YNKCsrw5YtW7Bu3Tp4vV7Mnj0bLS0twTIVFRUhX3/9618hSRK+/e1vAwCOHTuGY8eO4cEHH8SuXbuwevVqrF27FldffXWX+37xxRexfPly3H777di+fTu++c1vYt68eTh06FCf9jkgkX0HgFNOOSWk7V27dvVZXzvrjX4DwKRJk/Dkk0/i888/x3/+8x8IITB79mwogYlMNQyFcw50r+9A4s450Ht9D/jd734XnEQ4lkSe90T2Gxg653zJkiUh5R5//PEu9z1Uznm8/QYSe877jaCg6upqAUBs2LAhapkLL7xQzJw5s8t2/v73vwuz2Sw8Hk/UMqeffrr48Y9/HLLu5JNPFrfcckt8QfeS/uz7ihUrxIQJE7obaq/qrX7/73//EwDEvn37opYZqudcT98H0jkXomd937Fjhxg+fLioqKgQAMSaNWu63NdAOu/92e+hcs6nT58ufvrTn8a1r6FwzrvT74F2zvtKUl9xChf4iCk7O1tze1VVFV5//fWYV1QcDgcyMjJgNGpPBeh2u7Ft2zbMnj07ZP3s2bOxefPmbkTec/3V94C9e/eiuLgYo0aNwuWXX46vv/66e4H3UG/0u6WlBU8++SRGjRqFkpISzTJD9Zzr6XvAQDnnQPf73traiu9973tYtWqVro8gBtp5769+BwyFcw4Azz77LHJzc3HKKafgpptuQlNTU9T9DJVzDsTX74CBdM77TKIzt4FCVVWxYMECcfbZZ0ctc//994usrCzR1tYWtUxtba0YMWKEuP3226OWOXr0qAAg/vvf/4asv+eee8RJJ50Uf/A91J99F0KIN954Q/zzn/8UO3fuFOvWrRPTp08XBQUFora2ttt96I6e9ru8vFykpqYKAOLkk0/u8orLUDvn8fRdiIFzzoXoWd+vueYacfXVVweXEePKy0A67/3ZbyGGzjn/05/+JNatWyd27dolnn/+eTFy5Egxa9asqO0MlXMeb7+FGFjnvC8xcfL7yU9+IkpLS8Xhw4ejlhk9erRYunRp1O0Oh0OcccYZYu7cucLtdkctF/jB2rx5c8j63/zmN2L06NHxB99D/dl3Lc3NzaKgoEA89NBDcdXrqZ72u7GxUXz55Zdiw4YNYsGCBWLixIlRE8uhds7j6buWRJ1zIbrf91deeUWccMIJoqmpKbhOb+I0EM57f/Zby2A851o+/vhjAUBs27ZNc/tQOOdaYvVbSyLPeV9i4iSEWLp0qRg+fLj4+uuvo5bZuHGjACB27Nihud3pdIpp06aJ8847L+YfEJfLJQwGg3jppZdC1l9//fXinHPOib8DPdDffY9m1qxZEfcE9KXe6HdnLpdLpKSkiOeeey7q9qF0zjuL1fdo+vucC9Gzvv/0pz8VkiQJg8EQ/AIgZFkW06dP12xroJz3/u53NIPtnGtRVVWYTCbxwgsvaG4fCudcS6x+R5OIc97XqZE1fwAABrdJREFUkjpxUlVVlJWVieLiYvHll192Wfaqq64SkyZN0tzmcDjE1KlTxfTp00VLS4uufZ9++uniuuuuC1k3ZsyYfrt5MJF9D9fe3i6GDRsm7rrrrm7Vj0dv9Tucy+USNptNPPnkk1HLDJVzHk5P38P15zkXonf6XlFRIXbt2hXyBUD8/ve/7/KPUyLPeyL7HW4wnnMtgf53dbP1YD/nWvT0O1x/n/P+ktSJ03XXXSfsdrtYv369qKioCH61traGlHM4HCIlJUU89thjEW04nU5xxhlniPHjx4t9+/aFtOP1eoPlZs6cKR599NHg8gsvvCBMJpN44oknxO7du8Xy5ctFamqqOHDgQN91uJNE9v3GG28U69evF19//bXYsmWLmD9/vkhPT++XvvdGv7/66itx7733io8//lgcPHhQbN68WVx44YUiOztbVFVVBcsNxXPe3b4n8pwL0Tt916L1kdVAOu+J7PdQOOf79u0Td911l9i6davYv3+/eP3118XJJ58sTjvttAH7+z2R/U70Oe8vSZ04AdD8Cv/P+fHHHxc2m000NjZGtPHee+9FbWf//v3BcqWlpWLFihUhdcvLy0Vpaakwm81i4sSJcWXyPZXIvl922WWiqKhImEwmUVxcLC655BLx2Wef9VFPQ/VGv48ePSrmzZsn8vPzhclkEsOHDxdXXHGF+OKLL0LKDcVz3t2+J/KcC9E7fY/WbngCMZDOeyL7PRTO+aFDh8Q555wjsrOzhdlsFscff7y4/vrrRV1dXUi5oXbOu9vvRJ/z/iIJIYSep++IiIiIkh3HcSIiIiLSiYkTERERkU5MnIiIiIh0YuJEREREpBMTJyIiIiKdmDgRERER6cTEiYiIiEgnJk5EREREOjFxIiIiItKJiRMRDUq///3vMWrUKKSkpOCiiy6Cw+GIWraurg75+fk4cOBAxDa3240bb7wRN910E9xud7di+c53voOHH364W3WJaHBh4kREg85tt92GVatW4W9/+xs2bdqE7du346677opafuXKlViwYAFGjhwZse3FF1/ExIkTcdZZZ+Hpp5+O2L548WJcdNFFIev++c9/wmq14oEHHgAA/OpXv8I999wDp9PZo34R0cDHxImIBpWtW7fi/vvvx4svvohzzjkHEydOxLXXXovXXntNs3xbWxueeOIJ/OhHP9LcrqoqDAYDFEWBqqox9/+Xv/wF/+///T+sWrUKN998MwDg1FNPxciRI/Hss892v2NENCgwcSKiQeXBBx/EzJkzMXHixOC6vLw81NbWapZ/8803YTQaMW3aNM3tl19+ObZs2YL//ve/uPLKK7vc9wMPPIClS5fiueeei0jEFi5ciOeffz7O3hDRYGNMdABERHq5XC78+9//xoMPPhiyvq2tDXa7XbPOxo0bMXny5KhtWiwW/O53v4u571tuuQXl5eV47bXXMGvWrIjtp59+OlauXAmXywWLxRKzPSIanJg4EdGg8cknn6CtrQ033nhj8GMyAPB4PJgxY4ZmnQMHDqC4uLhH+33zzTfxyiuv4J133sHMmTM1ywwbNgwulwuVlZUoLS3t0f6IaOBi4kREg8aXX34Jq9WKXbt2haxfuHAhzjrrLM06bW1tsFqtPdrvqaeeitraWvzqV7/ClClTkJ6eHlHGZrMBAFpbW3u0LyIa2HiPExENGk6nE/n5+TjhhBOCX2azGV988QW+/e1va9bJzc1FQ0NDj/Y7bNgwbNiwARUVFZg7dy6ampoiytTX1wPw3W9FREMXEyciGjRyc3PhdDohhAiuu+eee3DBBRdg7NixmnVOO+007N69u8f7HjFiBDZs2IDq6mrMnj07YuiBTz/9FMOHD0dubm6P90VEAxcTJyIaNGbOnIn29nbcd999OHDgAO699168+uqreOyxx6LWmTNnDj777LMeX3UCgOHDh2P9+vWoq6vD7NmzQwbdfP/99zF79uwe74OIBjYmTkQ0aBQUFGD16tV47LHHMHbsWGzevBmbNm1CSUlJ1Drjx4/H5MmT8fe//71XYgh8bNfY2Ijzzz8fjY2NaG9vx5o1a7BkyZJe2QcRDVyS6HzNm4hoCHrjjTdw00034dNPP4Us9/7/i+Xl5XjllVfw1ltv9XrbRDSw8Kk6IhryLrjgAuzduxdHjx7t8upUd5lMJjz66KO93i4RDTy84kRERESkE+9xIiIiItKJiRMRERGRTkyciIiIiHRi4kRERESkExMnIiIiIp2YOBERERHpxMSJiIiISCcmTkREREQ6MXEiIiIi0omJExEREZFO/x9bcSvs4NCXfAAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 600x800 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, axs = plt.subplots(figsize=(6, 8), sharex = True)\n",
    "\n",
    "pot_temp_section_top = pot_temp_section.sel(st_ocean=slice(0, 3000))\n",
    "\n",
    "axs.set_title('Cross-slope section of (%.1f,%.1f) to (%.1f,%.1f)' % (shelf_coord[0], shelf_coord[1]-100, deep_coord[0], deep_coord[1]-100))\n",
    "\n",
    "cmesh = axs.pcolormesh(distance_in_km,\n",
    "                       pot_temp_section_top.st_ocean.values,\n",
    "                       pot_temp_section_top.values,\n",
    "                       cmap = 'RdBu_r')\n",
    "\n",
    "# Colorbar\n",
    "cbar = plt.colorbar(cmesh, orientation = 'horizontal')\n",
    "cbar.set_label(r'$\\theta$ ($^\\circ$K)')\n",
    "\n",
    "cs = axs.contour(distance_in_km[0, :],\n",
    "                 pot_rho_2_section.st_ocean,\n",
    "                 pot_rho_2_section,\n",
    "                 cmap = 'gray')\n",
    "\n",
    "axs.clabel(cs, cs.levels,\n",
    "           fontsize = 8,\n",
    "           colors = 'k',\n",
    "           inline = True,\n",
    "           use_clabeltext = True)\n",
    "\n",
    "# Axes\n",
    "plt.gca().invert_yaxis()\n",
    "plt.xlim([0, 350])\n",
    "plt.xlabel('Distance (km)')\n",
    "plt.ylabel('Depth (m)')\n",
    "axs.xaxis.set_minor_locator(MultipleLocator(10))\n",
    "axs.yaxis.set_minor_locator(MultipleLocator(50))\n",
    "plt.tight_layout();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "client.close()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [conda env:analysis3-24.01]",
   "language": "python",
   "name": "conda-env-analysis3-24.01-py"
  },
  "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.10.13"
  },
  "other_supplementary_files": [
   "images/topographic_gradient_sketch.png"
  ]
 },
 "nbformat": 4,
 "nbformat_minor": 4
}