{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "5ec17927",
   "metadata": {},
   "source": [
    "# Section 1.1 - Orbit Insertion DV"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "7263c18b",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from astropy.time import Time\n",
    "from AMAT.arrival import Arrival"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "bb17efc2",
   "metadata": {},
   "outputs": [],
   "source": [
    "arrival1 = Arrival()\n",
    "arrival1.set_vinf_vec_from_lambert_arc(lastFlybyPlanet='EARTH',\n",
    "                                      arrivalPlanet='MARS',\n",
    "                                      lastFlybyDate=Time(\"2024-10-01 00:00:00\", scale='tdb'),\n",
    "                                      arrivalDate=Time(\"2025-09-03 00:00:00\", scale='tdb'),\n",
    "                                      ephem_file='../../../spice-data/de432s.bsp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "4f884957",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Arrival v_inf vector, ICRF: [-0.99752618  1.63758919  1.73292484] km/s\n",
      "Arrival VINF MAG: 2.58 km/s.\n",
      "Arrival Declination: 6.03 deg.\n"
     ]
    }
   ],
   "source": [
    "print(\"Arrival v_inf vector, ICRF: \"+str(arrival1.v_inf_vec)+\" km/s\")\n",
    "print(\"Arrival VINF MAG: \"+str(round(arrival1.v_inf_mag, 2))+\" km/s.\")\n",
    "print(\"Arrival Declination: \"+str(round(arrival1.declination, 2))+\" deg.\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "ee0e39bc",
   "metadata": {},
   "outputs": [],
   "source": [
    "arrival2 = Arrival()\n",
    "arrival2.set_vinf_vec_from_lambert_arc(lastFlybyPlanet='EARTH',\n",
    "                                      arrivalPlanet='VENUS',\n",
    "                                      lastFlybyDate=Time(\"2010-05-10 00:00:00\", scale='tdb'),\n",
    "                                      arrivalDate=Time(\"2010-12-06 00:00:00\", scale='tdb'),\n",
    "                                      ephem_file='../../../spice-data/de432s.bsp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "2bcd54e1",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Arrival v_inf vector, ICRF: [-3.26907095  0.67649492 -1.06977471] km/s\n",
      "Arrival VINF MAG: 3.51 km/s.\n",
      "Arrival Declination: -21.93 deg.\n"
     ]
    }
   ],
   "source": [
    "print(\"Arrival v_inf vector, ICRF: \"+str(arrival2.v_inf_vec)+\" km/s\")\n",
    "print(\"Arrival VINF MAG: \"+str(round(arrival2.v_inf_mag, 2))+\" km/s.\")\n",
    "print(\"Arrival Declination: \"+str(round(arrival2.declination, 2))+\" deg.\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "d0bff732",
   "metadata": {},
   "outputs": [],
   "source": [
    "from AMAT.approach import Approach"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "7988c941",
   "metadata": {},
   "outputs": [],
   "source": [
    "approach1 = Approach(\"MARS\", v_inf_vec_icrf_kms=arrival1.v_inf_vec,\n",
    "                            rp=(3389.5+200)*1e3, psi=1.135*np.pi)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "3b2aeb68",
   "metadata": {},
   "outputs": [],
   "source": [
    "approach2 = Approach(\"VENUS\", v_inf_vec_icrf_kms=arrival2.v_inf_vec,\n",
    "                            rp=(6051.8+200)*1e3, psi=np.pi)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "b1ad1fee",
   "metadata": {},
   "outputs": [],
   "source": [
    "from AMAT.orbiter import PropulsiveOrbiter"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "405fd7b3",
   "metadata": {},
   "outputs": [],
   "source": [
    "apo_altitudes = np.linspace(200, 60E3, 1001)\n",
    "mar_OI_DV = np.array([])\n",
    "ven_OI_DV = np.array([])\n",
    "\n",
    "for alt in apo_altitudes:\n",
    "    orbiter1 = PropulsiveOrbiter(approach=approach1, apoapsis_alt_km=alt)\n",
    "    orbiter2 = PropulsiveOrbiter(approach=approach2, apoapsis_alt_km=alt)\n",
    "    \n",
    "    mar_OI_DV = np.append(mar_OI_DV, orbiter1.DV_OI_mag)\n",
    "    ven_OI_DV = np.append(ven_OI_DV, orbiter2.DV_OI_mag)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "980a292a",
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "a51d07b7",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "The PostScript backend does not support transparency; partially transparent artists will be rendered opaque.\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAksAAAE4CAYAAABL47fwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAACaFklEQVR4nOydd1wUx/vHP4CgqGADFE9FJcYebDEaa9RoNLHEFjXGGkss0Z81lkSjGBFEFFAEFbHSRIliJSJiQxBEIIgi0nuRJhzlbn5/3JeNK3W5O27vmPfrtS9vd2Znn/s8y97j7MwzaoQQAgqFQqFQKBRKpagr2gAKhUKhUCgUPkODJQqFQqFQKJRqoMEShUKhUCgUSjXQYIlCoVAoFAqlGmiwRKFQKBQKhVINNFiiUCgUCoVCqQYaLFEoFAqFQqFUQyNFG6AMNGvWDEKhEBoaGjAwMFC0ORQKhUKhUGRAeno6RCIRmjRpgvfv31dZT40mpawZDQ0NiMViRZtBoVAoFApFDqirq0MkElVZTnuWakF5sKSurg5DQ0OZt19WVoZGjagragvViztUM25QvVQH6kvu8E0zedqTkpICsVgMDQ2N6isSSo0IBAICgAgEArm0f+jQIbm0q6pQvbhDNeMG1Ut1oL7kDt80k6c9tf19p6/hakGHDh2QlJQEgUCAxMREmbdfUFCA5s2by7xdVYXqxR2qGTeoXqoD9SV3+KaZPO2p7e87nQ3HAw4fPqxoE5QKqhd3qGbcoHqpDtSX3OGbZnywh/Ys1QJ59yxRKBQKhUKpf2jPkgIQiUQQCoWcNysrqzqdp2pbdTMRPsTU1FTOnlQ9qGbcoHqpDtSX3OGbZnywh/Ys1YKaIk9CCFJTU5GTk1On9kUiUc0j8RsILVu2RLt27aCmplZlnczMTOjp6dWjVcoP1YwbVC/VgfqSO3zTTJ721LZniT9zA5WY8kDJwMAATZs2rfaHvjKysrLQpk0bOVmnHBBCUFhYiPT0dACoNkXDpUuXsHLlyvoyTSWgmnGD6qU6UF9yh2+a8cEeGixJiUgkYgKlugY8rVq1QpMmTWRsmfKhra0NQJJR1cDAoMretiFDhtSnWSoB1YwbVC/VgfqSO3zTjA/20DFLUlJaWgoAaNq0aZ3bKCkpkZU5Sk+5juW6VgYdZM8dqhk3qF6qA/Uld/imGR/socGSjOD66o1SObXRkQ6z4w7VjBtUL9WB+pI7fNJs586duH//frXrttUHvAyWcnNzMXfuXLRo0QLGxsY4e/ZshTrff/891NTUmM3f358pKy4uxsqVK7FmzRrMmTOHVVaOm5sb5s6di+XLlyt8pL2WlpZCr69sdOzYUdEmKB1UM25QvVQH6kvu8EWzpKQk/PXXXzh48CBGjhypUFt4GSzt3r0b8+fPh7e3Nzp16oTFixfj5cuXTPmrV6+gqakJKysrWFlZwcHBgfVOc/HixWjdujVsbW1hZWWFqVOnIj4+nim/desWtm/fDicnJzg4OCAgIABHjhyp1+/4IQUFBQq7tjLy5MkTRZugdFDNuEH1Uh2oL7nDF83c3d2ZXq7Jkycr1hi5LbhSR/Ly8sjr16+Z/ZiYGAKAuLu7M8fWrFlD0tLSKj0/ICCAACBRUVHMsQkTJpAlS5Yw+7179yY7duxg9p2dnYmuri4pKCiotM3q1o4pKioiERERpKioqPZf8iNKS0vrfK6qURs9MzMz69Ei1YBqxg2ql+pAfckdvmg2ZMgQAoAAIBEREXK5Rm3XhuNdz5KOjg66devG7Ldp0wa6uroYPXo0AMkKwadOncLEiRNx6NAhFBcXs86/dOkSdHR08MknnzDHBg4cCHd3d4jFYrx8+RL//vsvBgwYwCrPy8vDzZs35fvlqqB8unx9cffuXXTv3h3a2trIysqqtM748eOhq6uLixcv1jpZZH1hZ2enaBOUDqoZN6heqgP1JXf4oFlcXBwzhKZt27bo2bOnQu2RKli6ceMGbty4gdDQUACAo6Mj+vXrh9mzZ8ssALh48SLOnz/PJKQKCwvDxIkTkZ6ejo0bN2LIkCHIyMhg6gcEBFSYwm9gYID8/HxER0cjICAAAFh1DAwMAADPnz+Xic1cad++fb1eb+zYsZg5cyYAwMHBoUL5q1ev8OjRI3z22WeYN28e7xJm7ty5U9EmKB1UM25QvVQH6kvu8EEzNzc35vOvv/6qQEskSBUsTZ8+HdHR0ejcuTPu37+Pn3/+Gfr6+hg2bBi2bdsmlWExMTHYtGkTVq9eDVdXV2Z6/fjx4+Hh4YH4+Hg4ODggIiICS5cuZc5LT09H69atWW3p6OgAkGQBLQ/iPqzzYXl1pKSkoEOHDqzt9OnTUn1PAEhOTpa6Da5oampi3rx5sLOzQ1lZGavs+PHjmDlzJho14mcaLkUPyFdGqGbcoHqpDtSX3OGDZq6ursznmn6b6wOpgqUFCxZg7dq10NXVxbZt29ClSxdcu3YN69atk7q3pF27dpg7dy6mTZuGCxcuYO/evaxyNTU1LFu2DDY2Nrh27RqSkpKY4+XJDcspf42kqanJTE3/sM6H5dUhFouRlJTE2nx9fSESiZCdnY33798jNzcXZf37gwgEELVvD9KhA0SGhsD//iUdOkDUvj3EAgGztfv882rrEoEA4g/OqapuiYkJiouLkZWVhfz8fOTn5yM7OxtCoRCpqakQi8VMYJafn48lS5YgLS0Nzs7OyM3Nxbt375CRkYHU1FQIBAImQL1y5QqWL1+OjRs3Yvjw4fj3338RHx+P3377Df3798exY8fQsmVLhIWF4cSJEzh48CCmT5+OwYMHIzs7G3l5eSgoKEBmZiaKi4uRkpICQgiSk5NBCEFKSgqKi4uRmZmJ9+/fQygU4saNG4iKioKlpSWKioqYP15TU1PMnTsXDg4OCAoKwo0bN+Dp6YmwsDDY2toiJyeHVTc3Nxc2NjYIDw/HlStXcOvWLQQGBuLEiRNITU1l1RUKhbCwsEBUVBRcXFzg6+uLhw8f4ty5c4iLi8P+/fshEolY5yQkJMDJyQmPHz/G3bt34ebmhlevXsHS0hKFhYWsuhkZGbC3t0dwcDCuX7+Oq1evIjQ0FEePHkV2djarbl5eHqytrREREYHLly/j9u3bCAwMxKlTp5CcnMyqW1xcDHNzc0RHR8PZ2Rm+vr548OABzp8/j9jYWJiZmWH16tWscxITE3H69Gn4+/vD29sb7u7uiIyMhJWVFQoKClh1MzMzcfz4cYSEhMDLywvXrl1DSEgI7OzskJWVxaqbn5+Pw4cPIyIiAh4eHrhz5w6ePn0KR0fHCnaXlpbiwIEDePv2LS5cuAA/Pz/4+fnhwoULePv2LQ4cOIDS0lLWOcnJyXB0dMTTp09x584deHh4ICIiAocPH0Z+fj6rblZWFuzs7BASEoJr167By8sLISEhOH78ODIzM1l1CwoKYGVlhcjISLi7u6Nv377w9/fH6dOnkZiYyKpbVlYGMzMzxMbG4vz583jw4AF8fX3h7OyM6OhomJubo7i4uILdp06dQmBgIG7fvo3Lly8jIiIC1tbWyMvLY9XNzs7G0aNHERoaiqtXr+L69esIDg6Gvb09MjIyWHULCwthaWmJV69ewc3NDXfv3sXjx4/h5OSEhIQEVl2RSIT9+/cjLi4O586dw8OHD+Hr6wsXFxdERUXBwsICQqGQdU5qaipOnDiBwMBA3Lp1C1euXEF4eDhsbGyQm5vLqpuTkwNbW1uEhYXB09MTN27cQFBQEBwcHJCWlsaqW1RUBEtLS0RFRcHV1RU+Pj549OgRzp49i/j4eOzbtw9isRimpqYQi8XYt28f4uPjcfbsWTx69Ag+Pj5wdXWt8hmRlpbGPCN69OhBnxG1eEaUlZUx5xQXFyv0GfHmzRsEBQUBAAQCARYtWiS3Z0StUxJIMzBq69athBDJAGk1NTXi6urKlH399dfSNM1i+PDhpF+/fpWWlZWVEQMDA/Lo0SNCCCHjxo0jffr0YdU5dOgQAUBiY2PJ+fPnCQDy7Nkzpjw7O5sAILt37670GuUDwNTV1YlAIGBtx44dqzggWSAgBKj/rYYBah+ya9cuEhMTQ+bPn0+GDx/OHLe1tSX37t0jO3bsIKNGjSKEEDJw4EBy//59QgghkyZNIgcPHiQlJSXk7NmzREdHh9y9e5ecOnWKZGdnk5EjRzJtbd++vdb2lFObAd5WVlac223oUM24QfVSHagvuaNozUxNTZmB3ebm5nK1p7YDvKV6z9K1a1cMHjwYISEhmD59OmbPno2oqCiYm5vj7t270jTNYsaMGXB2dq60TENDA0ZGRujQoQMAwMTEhBlDVU5qairatWuHTp06wcTEBAB7UHVqaioAYPDgwdXaYWhoWCGTqFAoRExMDLtiu3Y1f6kPIITIJqklx+sCwNq1a/HFF1/g+fPn6N+/P/z8/LB69Wr8888/TB1ra2sMHDgQL168QGZmJgoKCqCpqYkOHTqgVatWGDNmDFM3ISEBK1euxJ49e7Bu3Trpv1MljB8/Xi7tqjJUM25QvVQH6kvuKFqzD8crzZ49W+EJKQEpX8MtX74cV65cwdOnT3Hp0iUAkmUqfv75Zzx+/FgmBgKSgGTEiBGVluXn56Njx47o1KkTAGD+/PlIT09nBTChoaGYPn061NTU0KdPH5iYmLASVYaGhkJPT4+ZcSc1z54BiYm13nLCwznVr3J79oyzqYMHD8bgwYNhbW2Nu3fvsgKfcgwNDfH7778jKysLPXv2ZPJelCcE/RBPT08EBASgS5cu8PDwqJt+NfBhzi1K7aCacYPqpTpQX3JH0ZpdvnwZu3fvxty5c2FkZKRwewCOwdK+ffvw+vVr1jGBQID+/fsz+7169cIXX3yBL774ok4GpaamwtLSkpnhlpKSgqtXrzIDxletWoWjR49CJBIhPz8fW7ZswYEDB5jz+/Xrh6lTp+LKlSsAJD0dL168YA043717N/7++2+IxWIAwJkzZ/Dnn39WGOtUX6irKzaDw9q1a+Hi4gJbW1vMnz+fVUYIwVdffYVffvml0kDqY3R1dfHs2TOYmZlhzZo1iI6Olrm95QPyKbWHasYNqpfqQH3JHUVrZmxsjF27duHixYu8sAcAOL2Gmz17NjPI69NPP8WsWbPQvXt3mRqUkZGBo0ePYu/evRgzZgw6deoENzc3Zqq/jo4OfvvtN1haWmLUqFEwNTWFQCBgtXHu3DmsWbMGGzduRFpaGq5fv868pgOAadOmIS0tDYsXL0bjxo0xYcIErFq1SqbfgwuKmJovEomYWXCzZ8/Gpk2b0KVLFzRr1oxVnp2djbi4OGRkZEBLSwsRERFo27Yt03P3cQ6mgwcPwtbWFmvXrsWZM2fkssZQixYtZN6mqkM14wbVS3WgvuQO3zTjgz2cgqVu3bph+/btAIA3b97A3d0dERERTODUo0cPqQ3q27cv3r59W2X5gQMHWD1JlaGjo4MzZ85UW2fFihVYsWJFnWyUNUKhEM2bN6+36/n5+eH69esoLCzEr7/+CiMjI6xcuRI//fQTAMDX1xc3b95ETEwM/vnnHyxatAjffPMNFi1ahO+++w4XLlzArFmzcOHCBWbmwZIlSwAAx44dQ35+Pnr16oX58+ezkoPKin///bfOPZcNFaoZN6heqgP1JXf4phkf7FEjMvivf3R0NNzd3fHvv/+iW7dumDVrlsKzbcqSDh06ICkpCQKBoMoB3l26dEGTJk3q1H5JSQldTPd/1EbP5OTkek/kqexQzbhB9VIdqC+5oyjNNm3ahGbNmmHBggUwNjauF3uq+33/EJkMljE2NsZvv/2Gc+fOYf78+bh69SoWLFggi6YbBHxIuKVMODo6KtoEpYNqxg2ql+pAfckdRWhWnqdrz549GDp0KCtZMh98KJOepcrw8PDAjBkz5NF0vSPvniXKf1A9KRQKpeFx4sQJLF++HIBkeZMjR47Uy3XrpWcpKSkJa9euxYQJEzBmzBhmGzVqFGsJEkr1KGK5E2WGD6n4lQ2qGTeoXqoD9SV3FKHZ2bNnmc8fv5nigw+l6lnq378/srKyMHz4cDRu3Jg5TgjBvXv3EBcXJxMjFY28e5bEYrHC0wfwhdroWVpaWuPSNBQ2VDNuUL1UB+pL7tS3Zm/evEG3bt0ASNIPhYeHs3L4ydOeeulZevv2LR4+fIiLFy/i9OnTzObk5ARbW1tpmm5QpKWlKdoEpeLQoUOKNkHpoJpxg+qlOlBfcqe+NTt58iTzeeHChRWSHfPBh1IFS0uXLmWSR35Mr169pGm6QdG6dWtFm6BUzJo1S9EmKB1UM25QvVQH6kvu1KdmJSUlOH36NADJYvaLFi1SqD1VIVWwZG5ujsuXLzMrAZdvvr6+2LVrl6xsVHkKCgoUbYJS8eTJE0WboHRQzbhB9VIdqC+5U5+aXb16lVmr9fvvv4eBgYFC7akKqRbS/fPPP7F//36YmZlVyNSspqaG8+fPS2VcQ4HmWOJGx44dFW2C0kE14wbVS3WgvuROfWrm4ODAfC6fDfcxfPChVMHSkSNHsG3bNvz000+sddVEIhGsrKykNo5CoVAoFIpqkpubi+DgYACSfI1fffWVgi2qGqlew/Xr1w//93//hx49esDIyIjZunbtih07dsjKRpWnpKSkXq936dIl6Ovro3HjxvDw8GCVWVlZQUtLixXt842EhARFm6B0UM24QfVSHagvuVNfmrVo0QIJCQk4d+4cTE1Nq5wVzgcfShUsHTt2DDdv3qy0zM/PT5qmGxT1uS4cAMycORPbt2+Huro6Jk6cyCqbMWMGVq5cWWV3KB8YOnSook1QOqhm3KB6qQ7Ul9ypT820tbUxf/58zJkzhxf2VIVUwdKWLVuwb98+VkLK8qSUfFmkVhnIzs6u92suWrQIampqcHd3Zx13dXXFypUr690eLnxsM6VmqGbcoHqpDtSX3OGbZnywR6pgqV27dtDR0YFAIGC9hjMyMqr33hJlpm3btvV+zVatWuGHH36Avb096/jLly/RrVs3mJmZYceOHRg8eDCuX7+O3Nxc/Pnnnxg0aBDu3LkDY2Nj9O3bF9nZ2YiLi8PUqVMxevRoAMCDBw/Qtm1bODk5AQDCw8Oxd+9ebN26Fdra2hAKhVLZvmHDBqnOb4hQzbhB9VIdqC+5Ux+aVZV2qDL44EOpgqU1a9bA2dkZ586dYyWlPHv2LM6cOSMrG5WSQ4cOoUOHDjVuU6ZMQWpqKuvcKVOm1OpcaRN1rVy5Ek+ePEFoaCgA4PHjxxg6dCgOHDiAYcOGYd++ffjtt9/w448/QkNDA5999hnevn0LoVCIV69eQV1dHe7u7jAyMsK0adOYdkeMGIGePXsy+3/++SdWrFiBAwcOYPv27VLZDAAHDhyQuo2GBtWMG1Qv1YH6kjvy1uzFixdo37495syZg6CgIIXbUxukmg03YMCAKsvGjBkjTdNKT15eHpKSkmqs17FjR7Rv3551LCMjo1bn5uXl1dk+APjiiy/Qv39/HD9+HMeOHYObmxv27t2LgQMHQltbG0FBQSgoKMCQIUOQlpaGli1bQldXF1OmTAEA9O3bl8k+/nHG1Q9p3rw5Zs2aBVtbW/zyyy9Sp0rYuXOnVOc3RKhm3KB6qQ7Ul9yRt2Y2NjYoKyuDq6srRo0ahYEDByrUntpAFySTE7q6uhAIBDVu+vr6FRbS1dfXr9W5urq6Utu5cuVKXLhwAWlpaRCLxdDR0UFCQgLmzZuH9evXY+fOnbh16xaMjY0rBESNGjWCWCyu8RqHDh2CgYEBTExMsH37dohEIqls5sOiisoG1YwbVC/VgfqSO/LULDMzExcuXAAgmQ33008/KdSe2iJVzxKlajZs2FDr96wfpw64evWqPEyqlB9//BGbN2/GrFmzsG/fPgCAoaEhPDw8sGbNGgBASkoKysrKqm1HXV29yiAoNzcX7u7u8PHxwQ8//IAhQ4ZgyZIldbZZmnMbKlQzblC9VAfqS+7IU7OTJ08y41aXLFlSq/HNfPAh7VniAdK+TpOGZs2aYf78+UhNTcWIESMAAHPnzsXWrVthbW0NHx8f/PnnnxAIBBCLxRUytZfvt23bFq9evUJ2djbu3buHhIQEZGRkoKysDIcPH0ZZWRnGjBmDBQsWVGiDK7du3ZLq/IYI1YwbVC/VgfqSO/LSrKysDMeOHQMgGbqxevVqhdrDBbn1LL148QKtW7fmRZpyvtOkSROFXv+XX35h+en3339HRkYGdu3aBWNjY5w5cwaFhYVwd3dHamoqrl+/js6dOyMoKAgxMTGIjY3FuHHjYGJigt69e+PAgQPo0aMHYmNjkZaWhpCQEIwbNw7fffcd1NTUatXtWh29e/eW9is3OKhm3KB6qQ7Ul9yRl2aenp5Mgslvv/0WxsbGCrWHC3IJlmbNmoXHjx9j48aNEIvF2LRpkzwuozJIO4ZHWvr06YM+ffow+02aNIGDg0OFLN52dnaws7Nj9sPCwljl3t7ezOcFCxYwn319fWVqb25urkzbawhQzbhB9VIdqC+5Iw/NCCGwtLRk9teuXatQe7gil9dwCQkJaNy4MTZs2IC5c+fCy8tLHpdRGWozSJryH/n5+Yo2QemgmnGD6qU6UF9yRx6aPXr0CP7+/gAkM6m//vprhdrDFbn0LPn5+THjUgQCAVq2bCmPy6gMin4Np2x8mMOJUjuoZtygeqkO1JfckYdm5ubmzOfNmzdXm26mPuzhCqeepWvXrtWqnpaWFho3bszsN2vWjJtVDQxFDvBWRu7cuaNoE5QOqhk3qF6qA/Uld+Sh2Q8//IC+ffuiY8eO1a4DV1/2cEWNcJia1K5dOxw/fpyVrbkh0KFDByQlJUEgECAxMZFVJhQKERMTgy5dutS5h0gkEkFDQ0MWpio9tdEzPz8fOjo69WyZckM14wbVS3WgvuSOvDQjhCAxMZHzxC95+rC63/cP4dSz1Lx5c9y9exdjxozBn3/+ifj4eKkNpYDJgk2pHUeOHFG0CUoH1YwbVC/VgfqSO/LSTE1NrU4z5PngQ049Sz4+PswyJn5+fnByckJ2djZ++OEHTJ8+nfXqTZWoTc+SkZERmjZtqiALVYfCwkLExcVJ1VNHoVAoFEptkEvP0ofrvY0cORKOjo44f/48CgsLMX36dKxZs6ZWi+KpElpaWlBXV0dycjJyc3NRVFQEoVDIaYuPj+d8jqptRUVFyM3NRXJyMtTV1atdP44Pqe+VDaoZN6heqgP1JXdkpVlWVhZmzZqFwMBAXtgjDZx6lmrCxcUFy5YtQ9euXbFkyRKsW7euTu3k5uZi5cqVuHHjBvT09LBr1y5W3p64uDhs2bIF7du3R0ZGBg4ePIh27dox5YQQ7Ny5E+np6SgsLMScOXMwefJk1jV8fHxw/PhxGBgYQFNTExYWFmjUqPLJgTVFniUlJUhJSUFhYWGdvi8ds/QfTZs2haGhYbXBUlZWFtq0aVOPVik/VDNuUL1UB+pL7shKs507dzLLaNnY2DBLaCnKnsqobc8SiJQUFRWR06dPky+++IKoq6uTxo0bkzlz5pD79+/Xuc3169cTLy8v8vTpUzJ69Giirq5OIiIiCCGEvH//nhgbG5Pbt28TQgjx8PAggwYNImVlZcz527ZtI3PnzmXqd+zYkQQGBjLlISEhpG3btiQjI4MQQsjatWvJ+vXrq7RHIBAQAEQgEFRZRywWk5KSElJUVMR5c3R0rNN5qraVlJQQsVhc4/1x7NixGutQ2FDNuEH1Uh2oL7kjC82ysrKIjo4OAUA0NTVJXFycQu2pitr8vhNCCKdgKSAggPn877//kl9//ZW0atWKqKmpka5duxIzMzOSnp5eN4v/R15eHnn9+jWzHxMTQwAQd3d3Qggh5ubmrC9VVlZGmjVrRs6ePUsIISQxMZFoaWkRb29vps6KFSvImDFjmP1JkyaRH3/8kdl/8uQJ0dDQIDExMZXaVFsx68rz58/l0q6qQvXiDtWMG1Qv1YH6kjuy0Oz3338nAAgAsnz5coXbUxW1/X3nNGZp8+bNsLS0xPDhw9G3b18cO3YMo0aNws2bNxEdHY2tW7dCX1+/rr1hAAAdHR1069aN2W/Tpg10dXUxevRoAMClS5cwYMAAplxDQwMmJiZwdXUFAFy9ehUlJSWsOgMHDsS9e/eQnp6OvLw83Llzh1Xer18/EEJw6dIlqWyvK+Vr5VBqB9WLO1QzblC9VAfqS+5Iq1lOTg4zg61Ro0bYtm2bQu2RBZwyePv5+eHBgwdo3749/vjjDyxbtgzt27eXl20AgIsXL+L8+fPQ09ODSCRCUFBQhYVYDQwMEBAQAAAICAiAuro6WrVqxSonhODFixfQ0tJCWVkZ6/1nkyZNoKuri+fPn1drS0pKCjp06MA6tmHDBmzYsEGq78glkymF6lUXqGbcoHqpDtSX3JFWsyNHjjDJlhctWoTOnTsr1B5ZwKlnqV27dnB3d0dcXBx27dol10ApJiYGmzZtwurVq+Hq6oqSkhJkZ2dDJBKhdevWrLo6OjrIzMwEAKSnp6Nly5YsccuTWWVmZiI9PR0Aqm2jKsRiMZKSklibl5cXkpOTcerUKQQGBuL27du4fPkyIiIiYG1tjby8PGYkv6mpKbKzs3H06FGEhobi6tWruH79OvLz82Fvb4+MjAxW3cLCQlhaWuLVq1dwc3PD3bt38fjxYzg5OSEhIYFVVyQSYf/+/YiLi8O5c+fw8OFD+Pr6wsXFBVFRUbCwsIBQKGSdk5qaihMnTiAwMBC3bt3ClStXEB4eDhsbG+Tm5rLq5uTkwNbWFmFhYfD09MSNGzcQFBQEBwcHpKWlseoWFRXB0tISUVFRcHV1hY+PDx49eoSzZ88iPj4e+/btg1gshqmpKcRiMfbt24f4+HicPXsWjx49go+PD1xdXREVFQVLS0sUFRWx2tfW1oaDgwOCgoJw48YNeHp6IiwsDLa2tsjJyWHVzc3NhY2NDcLDw3HlyhXcunULgYGBOHHiBFJTU1l1hUIhLCwsEBUVBRcXF/j6+uLhw4c4d+4c4uLisH//fohEItY5CQkJcHJywuPHj3H37l24ubnh1atXsLS0RGFhIatuRkYG7O3tERwcjOvXr+Pq1asIDQ3F0aNHkZ2dzaqbl5cHa2trRERE4PLly7h9+zYCAwNx6tQpJCcns+oWFxfD3Nwc0dHRcHZ2hq+vLx48eIDz588jNjYWZmZmaNeuHeucxMREnD59Gv7+/vD29oa7uzsiIyNhZWWFgoICVt3MzEwcP34cISEh8PLywrVr1xASEgI7OztkZWWx6ubn5+Pw4cOIiIiAh4cH7ty5g6dPn8LR0bGC3aWlpThw4ADevn2LCxcuwM/PD35+frhw4QLevn2LAwcOoLS0lHVOcnIyHB0d8fTpU9y5cwceHh6IiIjA4cOHkZ+fz6qblZUFOzs7hISE4Nq1a/Dy8kJISAiOHz+OzMxMVt2CggJYWVkhMjIS7u7uSExMhL+/P06fPo3ExERW3bKyMpiZmSE2Nhbnz5/HgwcP4OvrC2dnZ0RHR8Pc3BzFxcUV7JbmGREcHEyfERyeEWlpacwzIiYmhj4javGMKCsrY865efNmnZ8Rrq6usLCwAACoq6tj5cqVUj8jAMjtGfH+/XvUCi7v9k6ePCnFm0FuFBYWkmfPnpEZM2YQAGTnzp0kIyODACA7duxg1Z03bx5p1qwZIYSQb7/9tsK7xzt37jDjntzd3QkA1pgmQghp3749+fbbbyu1pfydprq6OhEIBKzN0tJS6u9qZ2cndRsNCaoXd6hm3KB6qQ7Ul9yRRrNdu3YxY5UWL16scHtqorZjlji9hlu6dGmt66alpaFt27Zcmmehra2NgQMH4tKlSxgxYgS8vLywZ88eaGlpITc3l1W3oKCAGStlaGiI+/fvVygHAH19fSY9QHVtVIWhoWH1UwvryMyZM2XepipD9eIO1YwbVC/VgfqSO3XVLD09HZaWlgAkY5V27NihUHtkCafXcDUhFotx7do1TJs2DZ06dZJZuzNmzICWlhbU1NTQt29f5lVaOampqRg8eDAAwMTEBAUFBaycR6mpqWjUqBEGDBiAnj17QlNTk9VGYWEh8vLymDbqm+PHjyvkusoK1Ys7VDNuUL1UB+pL7tRVs+zsbPTp0wcAsGzZMhgbGyvUHlkik6SUb968gaOjI86cOYPU1FS0atUKGhoaMlvzzMzMDJmZmTh48CAOHz4MKysrxMXFAQDKysrQunVrnDp1CrNmzUJGRgYEAgFu376Nr776CgDwyy+/ID4+HtevXwcATJs2DS1atMCZM2cAAE+fPsXIkSMRHR1dYQA3wCFpFYVCoVAoDRhCCK5evYrBgwfD0NBQ0ebUiFyWO/kQoVCI8+fPY/To0ejevTvMzc3Rt29fODs7IykpiTU1nwupqamwtLRERkYGAMkMtKtXrzJTD5cuXQpCCIKDgwEAly9fRq9evTBjxgwAkldtq1atwpUrVwBIViu+efMm9u7dy1xj586duHv3LtP7dObMGaxbt67SQKk+4EMqd2WC6sUdqhk3qF6qA/Uld6TRTE1NDVOnTpVpoMQHH3LuWQoODsbJkyfh7OyM3NxcGBsbY9GiRVi4cCEr2Jg0aRJu3LjB2aCwsDBMnToV2dnZGDNmDDp16oRNmzax2n758iX++OMPdO7cGVlZWbCwsGClAigtLcWGDRugpqaGd+/eYdGiRRg7dizrOp6ennB1dYWenh5atGiBPXv2QF298thR3j1LBQUFaN68uczbVVWoXtyhmnGD6qU6UF9yh2+aydMeuSx3MnPmTKKmpkaaNWtGFi5cSHx9fausO3HiRC5N8xp5Z/A+dOiQXNpVVahe3KGacYPqpTpQX3KHq2a//fYbuXv3rpyska8P5ZLB293dHbdu3cLEiRNRUFAAoVAoXUhHAQBMnDhR0SYoFVQv7lDNuEH1Uh2oL7nDRTNfX1+YmZlh7NixWLJkicLtkRecxyyNHz8e7u7uOHLkCB4/foyRI0fi4MGDyM7Olod9DYKwsDBFm6BUUL24QzXjBtVLdaC+5E5tNROJRFi/fj2zP2LECIXaI0/qPMBbIBDgzz//xL1799C5c2fMnz8fCxcuxJMnT2RpX4OgZcuWijZBqaB6cYdqxg2ql+pAfcmd2mrm6OiIFy9eAAAGDBiAhQsXKtQeecIpKWVlaGhoYObMmZg5cyZevXoFOzs7bNu2DfHx8bKwr0FQvhwLpXZQvbhDNeMG1Ut1oL7kTm00y83NZSWdPHz4cJWTpOrDHnkj02/WvXt3HD58GDdv3sTOnTtl2bRK8/LlS0WboFRQvbhDNeMG1Ut1oL7kTm00K1/LDgBmzZolt1dwtbVH3sgkKaWqI+/UAYmJiQrL8aSMUL24QzXjBtVLdaC+5E5Nmr158wa9evVCaWkpGjdujMjISHTu3Flh9kiD3JNSUmSHk5OTok1QKqhe3KGacYPqpTpQX3KnJs02bNiA0tJSAMCmTZvkGijVxp76gPYs1QK63AmFQqFQKMDff/+NadOmAZAsLv/69WteJbDkisJ7lrKysuTVtMrBh1TuygTViztUM25QvVQH6kvuVKdZnz598M033wAArKys6iVQ4oMPpe5ZIoQgLS0NJSUlzDGRSARra2tYWVlJbSAfkHfPUllZGRo1knpiYoOB6sUdqhk3qF6qA/Uld2rSjBACX19fjB49Gmpqagq3RxrqpWfJyckJrVu3hkAgQJcuXZjN2NgY1tbW0jTdoDh48KCiTVAqqF7coZpxg+qlOlBfcqcmzdTU1PDVV1/VS6BUG3vqA6lCtfXr12P69OmYOXMmmjVrxhwnhMDR0VFq4xoKc+bMUbQJSgXViztUM25QvVQH6kvufKyZWCyGUChE06ZNeWGPIpCqZ6lr164wMzPDpEmTMGrUKGYbPXo09uzZIysbVZ6HDx8q2gSlgurFHaoZN6heqgP1JXc+1szR0RG9evWCl5cXL+xRBFIFSydPnoSHh0elZQEBAdI03aAwMjJStAlKBdWLO1QzblC9VAfqS+58qFlSUhI2btyIuLg4TJ48Gc+fP1eoPYpCqtdwO3fuxOvXr+Hq6spKcy4WixESEoLZs2dLbWBDQCQSKdoEpYLqxR2qGTeoXqoD9SV3yjUjhGDlypXIy8sDACxatAj9+/dXmD2KRKpgqV27digoKICRkRErWCotLcWbN2+kNq6hkJKSomgTlAqqF3eoZtygeqkO1JfcKdfs4sWLzKu3du3a4dChQwq1R5FIFSytWbMGBgYGlaYh9/HxkabpBsXgwYMVbYJSQfXiDtWMG1Qv1YH6kjuDBw9GWloafv31V+bY8ePH0apVK4XZo2ikGrM0YMCAKtdrGTNmjDRNNyiqGvdFqRyqF3eoZtygeqkO1JfcuXTpElavXo3s7GwAktloU6dOVZg9fPCh1EkpxWIxLly4gMDAQGhqamLs2LGYNGmSrOzjBfJOSllcXIzGjRvLvF1VherFHaoZN6heqgP1JXdOnz6NJUuWAAD09PQQEREBfX19hdkjTx/WS1LKlJQUmJiYYOHChTh+/DhcXFwwY8YMDBs2DLm5udI03aCwsLBQtAlKBdWLO1QzblC9VAfqS27ExcXhl19+YfaPHTum0EAJ4IcPpepZmjVrFgBgx44dMDExYbJ5Xr16FXfu3IGtra1srFQwdCFdCoVCoTQEbt++jVmzZiE/Px8LFizAmTNnFG2SXKmXniWhUAh3d3f069ePlfZ8ypQp9ZYGXRXgwyKBygTViztUM25QvVQH6ktuTJgwAcuWLcPcuXN5s2wZH3wo1Wy4nj17Vnq8uLgYwcHB0jTdoCh/N0ypHVQv7lDNuEH1Uh2oL7mzceNGtG/fXtFmMPDBh1L1LDVq1Ah2dnZIS0uDUChEXFwczpw5g6FDh6Jbt26yslHluXnzpqJNUCqoXtyhmnGD6qU6UF/WTFlZGT4ckcM3zfhgj1TB0u7du+Hn5wdDQ0M0a9YMXbt2xeLFi6Gvr48jR47IykaV57PPPlO0CUoF1Ys7VDNuUL1UB+rLmlm2bBkWLlyI9+/fA+CfZnywR6rXcFpaWnB2dsb27dvh6+sLQggGDx6MIUOG4N27d7KyUeUpz2VBqR1UL+5QzbhB9VIdqC+rx8nJCU5OTgCAyMhI+Pv7804zPtjDKVhKTk5Gy5Yt0bRpUwBAfHw8AKBFixashFUxMTGwtraGlZWVDE1VXcqjeUrtoHpxh2rGDaqX6kB9WTXh4eFYtWoVs79u3Tqoq6vzTjM+2MMpWOrbty8mTpyI8+fPAwBGjx6NuLi4CvUIIVBTU6PBUi3p0aOHok1QKqhe3KGacYPqpTpQX1ZOQUEBZs2ahaKiIgCSV3E//vgjAP5pxgd7OI1ZunXrFszMzJj9VatW4a+//oK3tzd8fHyYzdvbG7Nnz66zUampqZgxYwZatGiBTz75BHZ2dhXqfP/991BTU2M2f39/pqy4uBgrV67EmjVrMGfOHFZZOW5ubpg7dy6WL1+u8GmJ//zzj0Kvr2xQvbhDNeMG1Ut1oL6sCCEEv/zyCyIjIwEAJiYmrHHGfNOMD/ZIlZQyJycHpaWllWb3TExMrHLduJqYNGkSBg8ejG7duuH06dO4e/cuzp49i59++gkA8OrVK/z+++/48ssvAQDNmjXDsmXLmPPnzZuHzp0746+//kJKSgr69euHwMBAdOrUCYAk6FuzZg3+/fdfNG7cGFOmTMHYsWOxbt26Su2Rd1LKvLw86OrqyrxdVYXqxR2qGTeoXqoD9WVFjh07htWrVwMAmjdvjqCgIHz66adMOd80k6c9tf59J1Jw69atSo/fuXOHpKen16nN8PBw4u7uzuyXlpaS7t27kxEjRjDH1qxZQ9LS0io9PyAggAAgUVFRzLEJEyaQJUuWMPu9e/cmO3bsYPadnZ2Jrq4uKSgoqLRNgUBAABCBQFCn71QTe/fulUu7qgrViztUM25QvVQH6ks29+/fJ40aNSIACADi4uJSoQ7fNJOnPbX9fZcqWPrll18qPf7u3TsydOjQOrWZnJxMxGIx69j69etJ7969mXJtbW0yYMAAYmlpSYRCIavuli1biI6ODuvY9u3biY6ODhGJRCQiIoIAIB4eHkz569evCQBWkPYh8g6WKBQKhUKRN/Hx8URfX58JlDZu3KhokxRObX/fOedZCg8Px5w5czBmzBhcvXoVY8aMqbANHDgQCQkJnLvDAMDQ0LDCUillZWUYMmQIACAsLAwTJ05Eeno6Nm7ciCFDhiAjI4OpGxAQgDZt2rDONzAwQH5+PqKjoxEQEAAArDoGBgYAgOfPn1drW0pKCjp06MDaDh06VKfv+SGKHjOlbFC9uEM14wbVS3WgvvwPXV1d5rf066+/Zo1B/hC+acYHezgHS3369IGDgwO0tbXRvHlzGBkZsbbOnTtj0qRJuHbtmsyMvHfvHjZs2AAAGD9+PDw8PBAfHw8HBwdERERg6dKlTN309HS0bt2adb6Ojg4AIDMzE+np6QDAqvNheXWIxWIkJSWxNi8vLyQnJ+PUqVMIDAzE7du3cfnyZURERMDa2hp5eXmMo01NTZGdnY2jR48iNDQUV69exfXr1zFs2DDY29sjIyODVbewsBCWlpZ49eoV3NzccPfuXTx+/BhOTk5ISEhg1RWJRNi/fz/i4uJw7tw5PHz4EL6+vnBxcUFUVBQsLCwgFApZ56SmpuLEiRMIDAzErVu3cOXKFYSHh8PGxga5ubmsujk5ObC1tUVYWBg8PT1x48YNBAUFwcHBAWlpaay6RUVFsLS0RFRUFFxdXeHj44NHjx7h7NmziI+Px759+yAWi2FqagqxWIx9+/YhPj4eZ8+exaNHj+Dj4wNXV1dERUXB0tISRUVFrPZnzpwJBwcHBAUF4caNG/D09ERYWBhsbW2Rk5PDqpubmwsbGxuEh4fjypUruHXrFgIDA3HixAmkpqay6gqFQlhYWCAqKgouLi7w9fXFw4cPce7cOcTFxWH//v0QiUSscxISEuDk5ITHjx/j7t27cHNzw6tXr2BpaYnCwkJW3YyMDNjb2yM4OBjXr1/H1atXERoaiqNHjyI7O5tVNy8vD9bW1oiIiMDly5dx+/ZtBAYG4tSpU0hOTmbVLS4uhrm5OaKjo+Hs7AxfX188ePAA58+fR2xsLMzMzFgTGUxNTZGYmIjTp0/D398f3t7ecHd3R2RkJKysrFBQUMCqm5mZiePHjyMkJAReXl64du0aQkJCYGdnh6ysLFbd/Px8HD58GBEREfDw8MCdO3fw9OlTODo6VrC7tLQUBw4cwNu3b3HhwgX4+fnBz88PFy5cwNu3b3HgwAGUlpayzklOToajoyOePn2KO3fuwMPDAxERETh8+DDy8/NZdbOysmBnZ4eQkBBcu3YNXl5eCAkJwfHjx5GZmcmqW1BQACsrK0RGRsLd3R09e/aEv78/Tp8+jcTERFbdsrIymJmZITY2FufPn8eDBw/g6+sLZ2dnREdHw9zcHMXFxRXsluYZERwcTJ8RHJ4RaWlpzDOiW7du9Bnxv2dEZmYmfvjhB6xcuRLr1q2Di4sL84woKytjziksLOTVM2LSpElye0bUOi1BXbuuhEIhcXBwqOvpteby5ctk3bp1VZbb29sTACQxMZEQIhmPNGzYMFYdBwcHAoAEBgYSCwuLCmOaSkpKCACyevXqSq9R3k2nrq5OBAIBa7O0tJT6O9ra2krdRkOC6sUdqhk3qF6qA/Uld/immTztkdtruHIaN26Md+/ewcPDo65N1Mi7d+9w/vx5mJubV1ln6dKlMDAwYPI9GRoaIjc3l1WnoKAAAKCvrw9DQ0MAYNX5sLw6DA0NkZiYyNrKe7ykYcSIEVK30ZCgenGHasYNqpfq0NB9eePGDcTGxnI6h2+a8cEeqdaGO3DgAEJCQiotI3XPSAAAEIlE2LJlC6ytraGlpVVlPQ0NDRgZGTFpCkxMTJhXbeWkpqaiXbt26NSpE0xMTACAVSc1NRUAMHjwYKlsritcb+SGDtWLO1QzblC9VIeG7MvHjx9j+vTpGDJkCJ49e1br8/imGR/skSpYsrOzg7GxcaVl0i6ku23bNqxatQoCgYA59nEQBAD5+fno2LEjk0Np/vz5SE9PR0xMDFMnNDQU06dPh5qaGvr06QMTExNWosrQ0FDo6elh9OjRUtlcVzQ0NBRyXWWF6sUdqhk3qF6qQ0P15Zs3bzBlyhQUFxcjLS0Njo6OtT6Xb5rxwR6pFtJ1cnJCVFQUbGxs0KJFC+Z4cXExgoKCsH79+jq1u2PHDojFYqSlpeHWrVsoKyuDn58fRo8eDS8vL/Tu3RsrV65EYWEhtmzZggMHDjDn9uvXD1OnTsWVK1ewYcMGJCQk4MWLFzhx4gRTZ/fu3di9ezd27doFdXV1nDlzBn/++Se0tbXrrIU0lL8apNQOqhd3qGbcoHqpDg3Rl1lZWZg0aRKysrIAAGPHjuXUgcE3zfhgj1TBkqGhITIzM9GzZ0+oq//XSSUWi5lFdrly8OBB/PXXXwAAS0tL5nirVq2wb98+3L9/H7/99hssLS0xatQomJqasnqfAODcuXNYs2YNNm7ciLS0NFy/fp2VTXzatGlIS0vD4sWL0bhxY0yYMIG1mGB9ExgYiAEDBijs+soG1Ys7VDNuUL1Uh4bmy6KiIkybNg1RUVEAgN69e+PSpUvQ1NSsdRt804wP9ki13ElwcDCaNWuG7t27Vyi7desWvvnmG6mM4wvyXu4kIyOjxsHllP+genGHasYNqpfq0JB8WVpaiu+//x7Xr18HALRr1w7+/v4wMjLi1A7fNJOnPbX9fZdqzFK3bt1gaWmJRYsWAZB0/Zmbm8PHx0dlAqX6wN7eXtEmKBVUL+5QzbhB9VIdGoovxWIxlixZwgRKOjo68PLy4hwoAfzTjA/2SNWzVO6YkSNHwt3dnTk+Y8YMTJ8+HT/++KNMjFQ08u5ZolAoFAqlrhBCsH79elhbWwOQpPa5deuWwiYtKRP10rP0/PlzvHr1CoMGDWIdHz9+PHbt2iVN0w0KPqRyVyaoXtyhmnGD6qU6NBRftmvXDgCgrq4OV1dXqQIlvmnGB3ukGuD95ZdfomXLlhWOP3jwoNJp/pTKkUViy4YE1Ys7VDNuUL1Uh4bgSzU1NWzbtg2tW7dG48aNMXXqVKna45tmfLBHqp4lAwMDpKSkMAvf5uTkYPPmzbh48SK+//57mRjYELCzs1O0CUoF1Ys7VDNuUL1Uh4bkyxUrVjBjiKWBb5rxwR6pepa2bt2KLVu24Nq1a3ByckJsbCyEQiG+//572NjYyMpGlee7775TtAlKBdWLO1QzblC9VAdV9aWdnR0EAgGmTJki87b5phkf7JGqZ6lJkyawtrbG/fv3YWlpiYsXL+LVq1fw8PCArq6urGxUeV68eKFoE5QKqhd3qGbcoHqpDqroSzs7O6xatQozZ87E33//LfP2+aYZH+yRKli6ceMGbty4gZycHEycOBHZ2dmYNWsWZs+eTccscaBNmzaKNkGpoHpxh2rGDaqX6qBqvnRwcGCSKJeWliIoKEjm1+CbZnywR6pgafr06YiOjkbnzp1x//59LFu2DPr6+hg2bBi2bdsmKxtVHkUts6KsUL24QzXjBtVLdVAlXx45cgQrVqxg9n/77Tf8+eefMr8O3zTjgz1SBUsLFizA2rVroauri23btqFz5864du0a1q1bh/bt28vKRpXn9evXijZBqaB6cYdqxg2ql+qgCr4khGDv3r2s9VY3b96Mv/76i5lgJUv4phkf7JFqgHfr1q0BAC4uLvD394eLiwuaNGkCAHj69Kn01jUQxo4dq2gTlAqqF3eoZtygeqkOyu5LQgg2b97MWiv1jz/+wO7du+USKAH804wP9kjVs9S1a1cMHjwYCxYswPTp0zF79mxERUVh2bJluHv3rqxsVHnOnDmjaBOUCqoXd6hm3KB6qQ7K7EuRSIQVK1awAiVLS0v8+eefcguUAP5pxgd7pFruBACSkpKQnp6O/v37AwAiIiKQn58PAPjiiy+kt5AH0OVOKBQKhVLfhIeH4/PPP4dQKISamhrs7e2xbNkyRZulUtTLcicAIBAImEAJAHr16oUvvvgCn3zyibRNNxj4kMpdmaB6cYdqxg2ql+qgzL7s06cPLl++jKZNm+LixYv1FijxTTM+2CN1zxIhBGlpaSgpKWGOiUQiWFtbw8rKSmoD+YC8e5ZEIhE0NDRk3q6qQvXiDtWMG1Qv1UEVfJmeng4DA4N6ux7fNJOnPfXSs+Tk5ITWrVtDIBCgS5cuzGZsbMysfkypGXNzc0WboFRQvbhDNeMG1Ut1UCZfhoaGYs+ePfi4D6M+AyWAf5rxwR6pZsOtX78e06dPx8yZM9GsWTPmOCEEjo6OUhvXUJg3b56iTVAqqF7coZpxg+qlOiiLL//55x9Mnz4d+fn5aN68uUIXj+WbZnywR+rZcGZmZpg0aRJGjRrFbKNHj8aePXtkZaPK4+fnp2gTlAqqF3eoZtygeqkOyuBLe3t7TJw4kZkc5e7ujtLSUoXZwzfN+GCPVMHSyZMn4eHhUWlZQECANE03KLp06aJoE5QKqhd3qGbcoHqpDnz2ZWlpKVavXo2VK1eirKwMADBlyhTcvXsXmpqaCrOLb5rxwR6pXsPt3LkTr1+/hqurK9TV/4u7xGIxQkJCMHv2bKkNbAiU/5FQagfViztUM25QvVQHvvoyKysLs2bNwr1795hjGzZsgLm5ucIHV/NNMz7YI1Ww1K5dOxQUFMDIyIgVLJWWluLNmzdSG9dQSE1NVbQJSgXViztUM25QvVQHPvry33//xZQpU/D27VsAgJaWFuzt7bFo0SLFGvY/+KYZH+yRKlhas2YNDAwM0KFDhwplly9flqbpBsXAgQMVbYJSQfXiDtWMG1Qv1YFvvnzw4AEmTZqEgoICAEDbtm1x+fJlfPnllwq27D/4phkf7OE0Zik5ORmFhYXMvp6eHsRiMeLj41lbTEwMHjx4IHNjVRVPT09Fm6BUUL24QzXjBtVLdeCbL/v06QN9fX0AQP/+/REYGMirQAngn2Z8sIdTUso2bdpg4sSJOH/+PADJbLi4uLgK9QghUFNTg0gkkp2lCkTeSSmFQiGzADGlZqhe3KGacYPqpTrw0ZfPnz+HjY0NbGxsWGl3+ALfNJOnPXJJSnnr1i2YmZkx+6tWrcJff/0Fb29v+Pj4MJu3tzcd3M2BgwcPKtoEpYLqxR2qGTeoXqqDon3p4+OD5ORk1rH+/fvD0dGRl4ESoHjNWJSU4NCBA4q2QrrlTnJyclBaWsp0KX5IYmJipWOZlBG6kC6FQqFQuCASibBv3z78+eef+PLLL+Hj46PQdABKy8mTwB9/AD//DKxaBbRrJ9Pm62W5k5YtW1YaKJUbQKkdfFgkUJmgenGHasYNqpfqoAhfxsfHY/To0di1axfEYjEePnwIJyenerejrvDm/icEOHYMSEkB9u4FEhIUZorUC+nKg9TUVKxevRr//PMP9PX1sXHjRvzyyy9MeVxcHLZs2YL27dsjIyMDBw8eRLsPok1CCHbu3In09HQUFhZizpw5mDx5MusaPj4+OH78OAwMDKCpqQkLCws0alT55EB59yylpqay7KdUD9WLO1QzblC9VIf69qWbmxuWL1+O3NxcAIC6ujp2796NHTt2sFLs8Bne3P9PnwJDhkg+DxoEBAbK/BL10rMkL5YsWYK+ffvi2LFj6Ny5M1atWoVz584BAAoLCzF27FgsXboUVlZWmD59OiZPnswaTL5jxw7ExMTgxIkTOHHiBFavXo1nz54x5S9evMC8efNw7Ngx2NraQiQSYfPmzfX+Pcu5du2awq6tjFC9uEM14wbVS3WoL18WFBRg6dKl+OGHH5hAycjICH5+fvj999+VJlACeHT/29n99/mDDhOFQHhGeHg4cXd3Z/ZLS0tJ9+7dyYgRIwghhJibmxOBQMCUl5WVkWbNmpGzZ88SQghJTEwkWlpaxNvbm6mzYsUKMmbMGGZ/0qRJ5Mcff2T2nzx5QjQ0NEhMTEylNgkEAgKAdV1ZEhAQIJd2VRWqF3eoZtygeqkO9eHLR48ekW7duhEAzDZnzhzy7t07uV9bHvDi/s/MJKRxY0IAUqqjQ8j793K5TG1/36UKdW/fvl3pcW9vb2RkZNSpzdatW2PGjBnMfqNGjTBx4kRkZ2cDAC5duoQBAwYw5RoaGjAxMYGrqysA4OrVqygpKWHVGThwIO7du4f09HTk5eXhzp07rPJ+/fqBEIJLly7VyWZpycrKUsh1lRWqF3eoZtygeqkO8vZlbGwsRo4ciaioKABA8+bNcebMGVy8eBEtW7aU67XlBS/ufycnoLgYAJA4dizQtKlCzZEqWPr7778rPf75559j6tSpdWrT0NAQampqrGNlZWUYMmQIRCIRgoKC0KZNG1a5gYEBnj9/DkCygK+6ujpatWrFKieE4MWLF3j+/DnKyspYbTRp0gS6urpMG/VNUVGRQq6rrFC9uEM14wbVS3WQty87d+6M1atXAwC++OILPH/+HAsWLKjwO6ZMKPz+F4tZr+Cixo5VoDESOC93Eh4eDlNTU6Snp+P169eIjIysUCcuLg4lJSUyMRAA7t27Bzc3N2RnZ0MkEqF169asch0dHWRmZgIA0tPT0bJlS9aNqqOjAwDIzMxkBnFX10ZVpKSkVJjlt2HDBmzYsKFuX+x/dOvWTarzGxpUL+5QzbhB9VIdZO3L4uJiaGpqssYg/fXXX+jevTtWrFih8EVwZYHC7/9r14DoaMnnsWNhOHq0Qs0B6tCz1KdPHzg4OEBbWxvNmzeHkZERa+vcuTMmTZokswFiV65cwbhx49CrVy8mANLW1mbVEYlETP4KNTW1SssBQFNTs1ZtVIVYLEZSUhJr8/LyQnJyMk6dOoXAwEDcvn0bly9fRkREBKytrZGXl8dMwzQ1NUV2djaOHj2K0NBQXL16FdevX8fZs2dhb2+PjIwMVt3CwkJYWlri1atXcHNzw927d/H48WM4OTkhISGBVVckEmH//v2Ii4vDuXPn8PDhQ/j6+sLFxQVRUVGwsLCAUChknZOamooTJ04gMDAQt27dwpUrVxAeHg4bGxvk5uay6ubk5MDW1hZhYWHw9PTEjRs3EBQUBAcHB6SlpbHqFhUVwdLSElFRUXB1dYWPjw8ePXqEs2fPIj4+Hvv27YNYLIapqSnEYjH27duH+Ph4nD17Fo8ePYKPjw9cXV0RFRUFS0tLFBUVsdr/+++/4eDggKCgINy4cQOenp4ICwuDra0tcnJyWHVzc3NhY2OD8PBwXLlyBbdu3UJgYCBOnDiB1NRUVl2hUAgLCwtERUXBxcUFvr6+ePjwIc6dO4e4uDjs378fIpGIdU5CQgKcnJzw+PFj3L17F25ubnj16hUsLS1RWFjIqpuRkQF7e3sEBwfj+vXruHr1KkJDQ3H06FFkZ2ez6ubl5cHa2hoRERG4fPkybt++jcDAQJw6dQrJycmsusXFxTA3N0d0dDScnZ3h6+uLBw8e4Pz584iNjYWZmRn++ecf1jmJiYk4ffo0/P394e3tDXd3d0RGRsLKygoFBQWsupmZmTh+/DhCQkLg5eWFa9euISQkBHZ2dsjKymLVzc/Px+HDhxEREQEPDw/cuXMHT58+haOjYwW7S0tLceDAAbx9+xYXLlyAn58f/Pz8cOHCBbx9+xYHDhxAaWkp65zk5GQ4Ojri6dOnuHPnDjw8PBAREYHDhw8jPz+fVTcrKwt2dnYICQnBtWvX4OXlhZCQEBw/fhyZmZmsugUFBbCyskJkZCTc3d1hZ2cHf39/nD59GomJiay6ZWVlMDMzQ2xsLM6fP48HDx7A19cXzs7OiI6Ohrm5OYqLiyvYLc0zIjg4mD4jODwj0tLSmGeEra2tzJ4RPj4+MDIywu7du1nPiMuXL+Pbb7+Fubm50j4jysrKmHN+++03xT4jPkiK6d6xI9zd3eX2jHj//j1qRV0HRQmFQuLg4FDX02tFdnY2mT59OikuLiaEECIWi4mWlhZZs2YNq96UKVNI586dCSGE/Pzzz6R58+as8suXLxMAxNfXlzx8+JAAIJcuXWLV0dXVJYsWLarUjvIBYOrq6kQgELA2S0tLqb9nTk6O1G00JKhe3KGacYPqpTrIwpfZ2dlk+fLlzODtpk2bkqioKBlYx08Uev8/eUKIJMMSIb16ESISydUeuQ/wbty4MZYtW1Zl+bt37+raNABJT8+WLVtgbW0NLS0tAJJeo759+yI9PZ1VNzU1FYMHDwYAmJiYoKCggLXgb2pqKho1aoQBAwagZ8+e0NTUZLVRWFiIvLw8po2qMDQ0RGJiImuT9hUcANjY2EjdRkOC6sUdqhk3qF6qgzS+JITAxcUFPXv2hIODA3O8V69eKCsrk4V5vESh9/+HS61s3Aioq/Pj75FLBJaUlETefzB9Ly4urtLt7du3ZP369XWK8srZvHkzCQ4OZh1LS0sjVlZWpFOnTsyx0tJSoqOjQ9zc3AghhKSnpxNNTU3i4+PD1Fm5ciWZNGkSsz916lSyYMECZt/f359oaWmRhISESm2Rd+oACoVCofCL6OhoMmHCBFY6gObNm5MjR46QsrIyRZunmkRFEaKmJulVateOEKFQ7peUS89S3759sXz5cmZ/9OjR6NKlS4XN2NgY1tbWdQ7gduzYAbFYjLS0NNy6dQteXl7YsmULnj17hqVLl4IQguDgYADA5cuX0atXLybdgL6+PlatWoUrV64AAPLz83Hz5k3s3buXaX/nzp24e/cu0/t05swZrFu3TmFLtPAmtbySQPXiDtWMG1Qv1YGrLwsLC7Fr1y707t2blR7n+++/x8uXL/Hrr7+qxCDu6lDY/b9vn+QFHAD8+ivQuLFi7fkATsudBAYGwtDQkAkqDh48CJFIhM8//5x185SVleHkyZNwdnbmbNDBgwcrzabdqlUrpKWlQVNTEy9fvsQff/yBzp07IysrCxYWFqxUAKWlpdiwYQPU1NTw7t07LFq0CGM/mnro6ekJV1dX6OnpoUWLFtizZ0+VGVblvdxJTk6O0ubjUARUL+5QzbhB9VIduPpy//792L59O7PfoUMH2Nra1jkdjjKikPv/zRugRw9AJAJatQJiYoAWLeRuT61/36Xpvnr37h1JT0+vtKyqV1rKiLxfw9nY2MilXVWF6sUdqhk3qF6qA1dfFhQUEIFAQBo1akQ2bNhA8vLy5GQZf1HI/b9gwX8Du/furTd7avv7zjnP0oe0bNkS4eHh2Lt3L96+fYtmzZph/PjxWLhwocJeaSkjo0aNUrQJSgXViztUM25QvVSH6nyZlZWFp0+fYtKkScyxZs2a4fz582jXrh169OhRHybyjnq//1+/Bs6fl3xu3VryCk6R9lSCVBm83d3d0b9/f5w7dw5isRgAcPToUQwfPhwFBQUyMbAhEF2efItSK6he3KGacYPqpTpU5suioiKYmZnB2NgY06dPR0xMDKt89OjRDTZQAhRw/+/dK8naDUhmwOnqKtaeSpCqZ2nbtm0YPnw4PD090eJ/7xYB4NatW9i1axcsLS2lNrAhUJ4agVI7qF7coZpxg+qlOnzoS5FIhLNnz+KPP/5gjU/ZsWMHLl68qAjzeEm93v9hYUC59m3aAGvXKtaeKpCqZ6mgoADbt29nBUoA8M033yA/P18qwxoSbdu2VbQJSgXViztUM25QvVSHtm3bghCCmzdvol+/fliyZAkTKKmrq+Pnn3+GhYWFgq3kF/V2/xMi6Ukq71XasgX43/JkCrGnGqQKlvbt24fw8PAKx0UiUaVrxlEqJygoSNEmKBVUL+5QzbhB9VIdXFxc8NVXX2HSpEms36vvvvsOoaGhOHHiBAQCgQIt5B/1dv/fugV4e0s+GxlVGKtU7/ZUA6fXcMuXL6+QtfTp06cICwtjHUtMTORFJKgsNKQpqbKA6sUdqhk3qF6qgZeXFw5+mBEawOeffw4LCwteDBrmK/Vy/5eVAZs2/bdvZgY0aaI4e2qAU89SRkYGfHx88ObNG8TExCAmJgYGBgbM5/KttLQUu3fvlpPJqseJEycUbYJSQfXiDtWMG1Qv1WDChAlo1aoVAKBbt25wcXHB06dPaaBUA/Vy/586BURESD5/8QXwww+KtacGOCWl/Oeff9CpUyd8+umn8rSJd8g7KSWFQqFQpOP58+fw8/PDunXrWMc9PT2Rk5OD+fPno1EjqeY0UWRFZqYkAWVWlmT/0SPgyy8VYkptf9859SyNGzeOFSiZm5vDw8Oj7lZSAPAjlbsyQfXiDtWMG1Qv5eHBgweYOHEiBgwYgP/7v/9DRHlvxf8IDw/HokWLaKDEAbnf/1u3/hcozZlTY6DEh79HTj1LH9OmTRusWrWKte5aOYQQqKmpSWUcX5B3z1JRURG0tbVl3q6qQvXiDtWMG1QvflM+u+2vv/7Co0ePWGXLly+Hvb09s099yR25avbgATBypOSzri4QGQkYGirMHrn0LH2MnZ0djI2NKy07cuSINE03KI4dO6ZoE5QKqhd3qGbcoHrxE5FIBFdXV/Tv3x/ffvstK1Dq3Lkzjh07VuG3h/qSO3LTrKQEWLnyv/2//qoxUJKrPRyQql/SyckJUVFRsLGxYeVaKi4uRlBQENavXy+tfQ2CKVOmKNoEpYLqxR2qGTeoXvzD09MTmzZtqpDNuVevXti2bRt++OEHaGpqVjiP+pI7ctPM3Py/Qd2ff84OnBRhDwek6lkyNDREq1at0KdPHxgZGTHbJ598An19fVnZqPIEBwcr2gSlgurFHaoZN6he/OTDQGnw4MHw9PREWFgY5s+fX2mgBFBf1gW5aBYSAuzZI/msrg4cPw5oaCjOHo5I1bO0evVqNG3atNI1dG7duiVN0w0KGlhyg+rFHaoZN6heisXf3x+NGzdG//79mWOTJ0+GsbExunTpgm3btuGrr76q1bhY6kvuyFyz4mJg4UKgtFSy/9tvwIABirOnDkgVLA0YMABisRjnzp1DYGAgNDU1MXbsWEyaNAnffPONrGxUeRo3bqxoE5QKqhd3qGbcoHrVP8XFxbh8+TKsra3h7++Pb775Bjdv3mTKNTQ0EBwcDN2PFlmtCepL7shcsz17gNBQyefPPgP++EOx9tQBqV7DpaSkwMTEBAsXLsTx48fh4uKCGTNmYNiwYcjNzZWVjSoPH1ZUViaoXtyhmnGD6lV/xMbGYtu2bejYsSPmzZsHf39/AJK3Ex+nAeAaKAHUl3VBppo9fCjJzg0AmprA2bMAx+CHDz6UKlj69ddf0aNHDwQHB6O4uBhJSUkoKirC1q1bsWPHDlnZqPKMHj1a0SYoFVQv7lDNuEH1ki8ikQjXr1/Hd999h65du8LMzAwZGRlMed++feHo6IiuXbtKfS3qS+7ITLOsLGDu3P8Wyv3jD8DERHH2SIFUwZJQKIS7uzv69evHenc8ZcoUlcmxVB+cO3dO0SYoFVQv7lDNuEH1kh95eXn45JNP8N133+H69esoT/WnqamJH374Ab6+vnjx4gUWL16MJlWsFcYF6kvuyEQzQiTjlMpzF40cKRmrpCh7pESqpJRbtmyBubl5hePFxcUYM2ZMhWRhyoq8k1KKxWKoq0sVtzYoqF7coZpxg+olOypLUDx69Gjcv38fANCpUyesWLECS5culcsC7NSX3JGJZocOARs3Sj7r6UlmwwkEirOnCuolKWWjRo1gZ2eHtLQ0CIVCxMXF4cyZMxg6dCi6desmTdMNir/++kvRJigVVC/uUM24QfWSnn///RcbN27EF198AXH5a5j/sXr1akycOBHXrl3D27dvsX37drkESgD1ZV2QWrN794AtW/7bP3euzoGSTOyRAVL1LJWUlGDhwoVwdXVl/udACMHXX38NNzc3VqJKZYb2LPELqhd3qGbcoHrVjZycHLi6usLR0REBAQHMcW9vb4wbN04hNlFfckcqzWJiJAkny9d+27ZNkqlbUfbUQL30LGlpacHZ2RkvXrzA4cOHYWVlhcePH+P27dsqEyjVB/v371e0CUoF1Ys7VDNuUL1qj1AohIeHB6ZPn462bdti5cqVrEBJS0sLL1++VJh91JfcqbNmBQXA1Kn/BUoTJwKVrB1bb/bIEKl6lhoK8u5Zio+PR6dOnWTerqpC9eIO1YwbVK/asWXLFjg4OFSaKqZ///5YsmQJ5s2bh9atWyvAOgnUl9ypk2ZiMTBrFnD5smS/e3fA3x9o2VIx9tSSeulZosgGX19fRZugVFC9uEM14wbVq3a8f/+eFSi1bdsW69atQ3BwMIKDg7FmzRqFBkoA9WVd4KwZIZLB3OWBUosWwN9/yyRQqpM9coAGSzzA2NhY0SYoFVQv7lDNuEH1kkAIQXBwMLZv344ePXogNTWVVf7jjz9CR0cHCxcuxO3bt5GYmIjDhw+zlilRNNSX3OGsmaUlcPiw5LOGBuDsLOlZUpQ9ckCq5U4osqG4uFjRJigVVC/uUM240ZD1IoTg2bNnuHTpEi5duoS3b98yZVeuXMEvv/zC7A8dOhSpqalo2rSpIkytFQ3Zl3WFk2bnzgGbN/+3b28vGaukKHvkhFTBUn5+PjZu3Iji4mKcOXMGWVlZOHXqFAYNGoQxY8bIykaV58PMtZSaoXpxh2rGjYaml1gsRkBAABMgxcXFVaijrq5eYdkJNTU1XgdKQMPzpSyotWZeXsCSJf/t790LLF2qOHvkiFTB0rp163D9+nWMHDkSANCmTRts2bIFM2bMQEpKCn788UeZGKnqDOCw+jKF6lUXqGbcaGh6zZs3D66urhWOa2ho4KuvvsKsWbMwbdo0GBgYKMA66WhovpQFtdLs5k1gxgygrEyyv2oVIKdlzvjgQ6nGLD1//hyvXr3CoEGDWMfHjx+PXbt2SWUYIJmSeuTIEQwZMqTS8u+//x5qamrMVr4AIyDptlu5ciXWrFmDOXPmsMrKcXNzw9y5c7F8+XKYmppKbW9duXr1qsKurYxQvbhDNeOGquqVkZGBs2fPoqz8B+5/DB8+nPncqFEjTJgwASdOnEBqaiq8vb2xfPlypQyUANX1pTypUbM7d4DvvwdKSiT7c+YA1taAnJY544UPiRSsWrWKEEKImZkZ6/iPP/5IdHR0pGmalJaWEkdHR9KrVy9iZGRUoTwyMpLMmjWLWFlZESsrK+Lg4MAqnzt3Ltm2bRshhJDk5GRiYGBA4uLimPKbN28SY2NjIhQKCSGETJ48mRw+fLhSWwQCAQFABAKBVN+pKgoLC+XSrqpC9eIO1YwbqqKXWCwmkZGRxNzcnAwfPpyoq6sTAMTPz49VLzY2lnz33XfE0dGRZGVlKcha+aAqvqxPqtXszh1CmjQhRDIHjpBZswgpLVWcPVJS2993qYKl3bt3k+TkZHLgwAFCCCHv3r0jmzZtImpqamTBggXSNM2wdevWSoOlNWvWkLS0tErPCQgIIABIVFQUc2zChAlkyZIlzH7v3r3Jjh07mH1nZ2eiq6tLCgoKKrQn72Bp7969cmlXVaF6cYdqxg1l1is/P59cvXqV/PLLL6RLly4EQIVt06ZNijaz3lBmXyqKKjVzdydEU/O/QGn6dEJKShRnjwyol2CpqKiIrF27lnTu3Jn07NmTaGtrEzU1NTJ9+nSSm5srTdMMu3btqhAsJScnE21tbTJgwABiaWnJ9A6Vs2XLlgo9W9u3byc6OjpEJBKRiIgIAoB4eHgw5a9fvyYAiLu7ewUb5B0sUSgUiixYuHAh0dTUrDRAAkB69OhBtmzZQoKCghRtKkXZcHAgRF39v0Dp++8JKS5WtFVSU9vfd6nGLDVp0gTW1ta4f/8+LC0tcfHiRURGRsLDwwO6urrSNF0tYWFhmDhxItLT07Fx40YMGTKENVo+ICAAbdq0YZ1jYGCA/Px8REdHM6n4P6xT/j7++fPnVV43JSUFHTp0YG2HDh2S+vsocryUMkL14g7VjBt81ysnJwf//PNPheNNmzZFaWkps6+pqYkxY8bAwsICr169wsuXL3HgwAFeDJitL/juSz7C0owQ4MABYPlySZZuAFi8GHBzA7S06t8eBSFVsPT+/XsAQKdOnTBx4kRMmzYNn376qUwMq47x48fDw8MD8fHxcHBwQEREBJZ+MF0xPT29QtZYHR0dAEBmZibS09MBgFXnw/KqEIvFSEpKYm1eXl5ITk7GqVOnEBgYiNu3b+Py5cuIiIiAtbU18vLyGEebmpoiOzsbR48eRWhoKK5evYrr16/jiy++gL29PTIyMlh1CwsLYWlpiVevXsHNzQ13797F48eP4eTkhISEBFZdkUiE/fv3Iy4uDufOncPDhw/h6+sLFxcXREVFwcLCAkKhkHVOamoqTpw4gcDAQNy6dQtXrlxBeHg4bGxskJuby6qbk5MDW1tbhIWFwdPTEzdu3EBQUBAcHByQlpbGqltUVARLS0tERUXB1dUVPj4+ePToEc6ePYv4+Hjs27cPYrEYpqamEIvF2LdvH+Lj43H27Fk8evQIPj4+cHV1RVRUFCwtLVFUVMRqf9q0aXBwcEBQUBBu3LgBT09PhIWFwdbWFjk5Oay6ubm5sLGxQXh4OK5cuYJbt24hMDCQGbz6YV2hUAgLCwtERUXBxcUFvr6+ePjwIc6dO4e4uDjs378fIpGIdU5CQgKcnJzw+PFj3L17F25ubnj16hUsLS1RWFjIqpuRkQF7e3sEBwfj+vXruHr1KkJDQ3H06FFkZ2ez6ubl5cHa2hoRERG4fPkybt++jcDAQJw6dQrJycmsusXFxTA3N0d0dDScnZ3h6+uLBw8e4Pz584iNjYWZmRkWL17MOicxMRGnT5+Gv78/vL294e7ujsjISFhZWaGgoIBVNzMzE8ePH0dISAi8vLxw7do1hISEwM7ODllZWay6+fn5OHz4MCIiIuDh4YE7d+7g6dOncHR0rGB3aWkpDhw4gLdv3+LChQvw8/ODn58fLly4gLdv3+LAgQMoLS1lnZOcnAxHR0c8ffoUd+7cgYeHByIiInD48GHk5+ez6mZlZcHOzg4hISG4du0avLy8EBISguPHjyMzM5NVt6CgAFZWVoiMjIS7uzs+/fRT+Pv74/Tp00hMTGTVLSsrg5mZGWJjY3H+/Hk8ePAAvr6+cHZ2RnR0NMzNzVFcXFzBbmmeEY8ePcK6devw66+/okOHDmjTpg2+/vprxMbGsp4RhoaGMDQ0xFdffYWTJ09i69atuHv3LoRCIYyNjRvMMyItLY15RnTp0oU+I2rxjCgrK2POKSgoQGJiIpxOnEDajBnAb78xv4PPRo9GwZEjMDUzq7dnxPjx4+X2jCiPY2pEmu6rn376qdIxPrKkstdwH2Nvb08AkMTEREKIZDzSsGHDWHUcHBwIABIYGEgsLCwqjGkqKSkhAMjq1asrtF/eTaeurk4EAgFrs7S0lPo72tvbS91GQ4LqxR2qGTcUrVdpaSl5/Pgx2bt3Lxk9ejTR0tKq9LXauXPnWOeJRCIiFosVZDU/UbQvlRF7e3tCsrMJGTv2v9duACF//UWIAu4vefqwtq/hpMqzJBQK8dtvv6FJkyaYN2+ewlLcL126FL///jvi4uIgEAhgaGhYIS1/QUEBAEBfXx+GhoYAwFrT6MPyqjA0NJTLQroDBw6UeZuqDNWLO1QzbihSr9LSUhgaGiKrfOX2Sujbty/Gjh2LPn36sI6rq9MVrD6G3vvc+dLAABg6FHj1SnJASws4eRL46SeF2MMHH0oVLB07dgx6enrIz8/HxYsXYWdnh0GDBmHevHlo3ry5rGysEQ0NDRgZGaFDhw4AABMTE4SGhrLqpKamol27dujUqRNMTEwAgHkdV14OAIMHD64nq/8jLS2t3q+pzFC9uEM144a89SopKUFQUBD8/PwgFouxbds2pkxTUxO9e/eGn58fc6xLly4YO3Ysxo4di6+++gpt27aVq32qBL33OXL1KrrPnw+Uv57S0wM8PYFhwxRmEh98KFWwpKenB0Ay3mfFihVITk7GmjVrsHHjRsybNw/29vYyMbIm8vPz0bFjR3Tq1AkAMH/+fFhaWiImJgZdunQBAISGhmL69OlQU1NDnz59YGJiAn9/f0z83xo2oaGh0NPTw+jRo+vF5g8pKU/sRakVVC/uUM24IWu9CgsL4e/vz4y58Pf3R1FREQDJc/S3336D2gcJ/WbMmAGBQICxY8dizJgxzHOMwh1679eSsjJJBm5zc2iWH+vVS7KkiYLvP174UJp3fU+ePCGEEHLjxg0ybdo0oqmpSQwNDcnOnTtJfHy8NE0zbN++vcK7xF9++YXY2tqSsrIykpeXR1auXMkaf0QIIVOnTmXGE8XHxxNDQ0OSkJDAlF+5coWYmJgQkUhECCFk4sSJ5OjRo5XaIO/UAaGhoXJpV1WhenGHasYNWegVExNDNm/eTIYMGUIaNWpU5XR+ACQiIkIGVlMqg977tSA5mZCRI9njk2bOJERGKYCkRZ4+rJcxS99//z0AyeusMWPGwNnZGdOmTYOGhoY0zTJ4eHjg6tWrSElJweHDhzFv3jwYGBhAR0cHv/32GywtLTFq1CiYmppCIBCwzj137hzTy5WWlobr168zr+kAYNq0aUhLS8PixYvRuHFjTJgwAatWrZKJ3Vy5f/8++vbtq5BrKyNUL+5QzbjBRa/S0lKEhIRAX18fnTt3Zo4XFhbCwsKi0nOMjIwwYsQIjBw5EiNHjqyXWcQNFXrv18DffwPLlgHl6XcaNYLf1KkY6eYmt+VLuMIHH6oRQkhdT27WrBlWrVqFFStW4JNPPpGlXbyiQ4cOSEpKgkAgkMsA75ycHLRs2VLm7aoqVC/uUM24UZ1eqampePLkCbM9e/YMQqEQv//+O/bs2cPUE4vFaNOmDXJyctCjRw+MHDkSI0aMwIgRI2BkZFRP34RC7/0qyM8H/u//gFOn/jvWoQPg5oacnj15pZk8fVjb33eppk7Y2NjAwsKi0kDp3bt30jTdoLC1tVW0CUoF1Ys7VDNufKhXcHAwLCwsMHv2bHTp0gWGhoaYPn06LCws8PDhQwiFQgDAkydPWG2oq6vj+vXrSEtLw8uXL2Fvb4/58+fTQKmeofd+JTx6BJiYsAOlKVOA4GBg6FDeacYHezj1LCUnJ6Nly5Zo2rQpACA+Pr7SeiKRCNbW1rCyspKNlQpG3j1LFApF8bx//x5hYWEYMmQI6/i6detgbW1d5Xldu3bF0KFDMXbsWCxevFjeZlIodScvD9i+HTh2TDIyCQCaNweOHJFk5ebJa7f6RC49S3379sXy5cuZ/dGjR6NLly4VNmNj42ofLhQ2fEjlrkxQvbhDNWNTVFSEgIAAHDt2DEuWLMFnn30GXV1dDB06lJVFGAA+//xz5nPTpk0xcuRIbN26FZ6enkhNTUV0dDTOnz9PAyWeQu/9/3HlCtCzJ3D06H+B0rBhwIsXwJIlrECJb5rxwR5OPUuBgYEwNDRkBkofPHgQIpEIn3/+OWtQd1lZGU6ePAlnZ2fZW6wA5N2zlJubixYtWsi8XVWF6sWdhq5ZUVERbG1tERISgpCQEERGRkJcvs7VR/z9998YNWoUo1dSUhLu3LmDzz//HD169ECjRlLNi6HUMw393kd8PPDrr5KB3OU0bQrs3QusWwdUMiGLb5rJ055a/75LM+Xu3bt3JD09vdKyD6fpKzvyTh1gbW0tl3ZVFaoXdxqCZiKRiERFRRF3d3fi7e3NKisrKyPa2tpVTt3X0NAg/fr1Iz///DMJCAhoEHo1FBqsLwsKCPnjD0K0tdkpASZNIiQmptpT+aaZPO2pl9QB1Y1O/3CaPqV6vvrqK0WboFRQvbijaprl5OTg33//RXh4OMLCwvDixQu8ePEC+fn5AIDJkydj3LhxTH0NDQ189tlnePr0KZMh28TEBIMGDcKgQYNgYmICbW1tpv6HnynKjard+zUiFgMXL0oWv01K+u9427aAtTUwa1aNY5P4phkf7JEqWLpx4wYASWD02WefwdHREdbW1vj0009ha2sLAwMDmRip6kRFRVVY44lSNVQv7qiKZg4ODti7d2+Nr8NDQkIqHLOyskLTpk3Rs2dPaGlpVXu+quhFaWC+9PEBtm0DAgL+O9aokeQ13O+/A7Wcfs83zfhgj1TBUvn02eHDh+P+/fv4+eefMXbsWAwbNgzbtm3DqQ+nJVKqhP4vlhtUL+7wXbPS0lJERUUhPDyctd27d4+VcLZRo0ZVBkpGRkbo168fsxFCWEuIDB06tNb28F0vSu1pEL588gTYuVMSLH3I5MnAwYMAx6SnfNOMD/ZIFSwtWLAAa9euBQBs27YNXbp0wbVr19CkSRP8/vvvMjGwIdCmTRtFm6BUUL24wyfNSktL4ezsjMjISLx8+RKRkZF48+YNysrKKtQNDw9nBUt9+/ZFy5Yt0bdvX/Tp04fZ+vbti1atWsnMRj7pRZEOlfZlSIikx8jLi328Tx/A0hIYP75OzfJNMz7YI1Ww1Lp1awCAi4sL/P394eLigiZNmgAAnj59Kr11DYSQkBDW9GRK9VC9uFOfmhFCkJqaygRCAoEAU6dOZco1NDSwcuVKZiHZqtDW1q6w2vigQYOQnZ3N6jGSB/QeUx1U0pePHgFmZhWDJGNj4M8/gTlzKp3lVlv4phkf7JEqWOratSsGDx6MkJAQTJ8+HbNnz0ZUVBTMzc1x9+5dWdmo8kyePFnRJigVVC/uyEOzgoICREVFsbbIyEhERkYiNzeXqTdlyhRWsKSuro7u3bsz44oaN26M7t27o0ePHqweoy5dulRYZ1LeQVI59B5THVTGl4QAt24B+/cDDx6wyzp0AP74A1i0CNDUlPpSfNOMF/ZIO+0uMTGRBAcHM/v//vsv8ff3J48fP5a2ad4g79QBe/fulUu7qgrVizt11ez9+/fkxYsX5NKlS0QkErHKNm7cWOVU/A+3bt26VWj32rVrxMvLi0RHR5OysrI62SZP6D2mOii9L4uKCDl9mhATE3YKAICQDh0IOXJEUkeG8E0zedpT2993Tkkpq1re5GNEIhGOHDmCw4cPcwrc+Apd7oSiyhQWFiImJqZCL1FUVBSSPph6nJCQwEoJYm9vj5UrV1Zor3PnzujRowd69uyJHj16oHfv3hg2bFi9fBcKRWVITATs7AAHByAzk13WowewdSswbx5Qw8xOSvXIJSmlkZERUVdXr3FTU1Mj6urqUsR6/IL2LPELqhc3ysrKyKZNm4ivry+JjIxklcXHx9eqdwgA8fHxYZ0bFBREFi9eTP766y/i7u5OQkJCyPv37+vzq8kNeo+pDkrlS7GYkHv3CJk9mxANjYo9SZ9/ToiHByEf9fLKGr5ppnQ9S+bm5lBXV8eAAQMqjCX4ELrcCTeEQiEzMJ5SM1SviuTm5uLt27eIiYnB27dvmS0mJgaxsbEoKSkBAGzduhVmZmbMeSKRCE2aNKl0Jpqenh66devGbHPnzoWxsXG9fSdFQu8x1UEpfJmcDJw5A5w6BURHs8saNZIkkvz1V+CLL+plsVu+aSZPe2r7+85pgPeyZcsgEomgp6dXY90ePXpwabpBY2Njg82bNyvaDKWhoelFCEF2djbi4uIQFxeH/Px8LFiwgFVn3rx5TJLY6nj79i1rX0NDA9OnT4eWlhY++eQTJjD65JNPZDoVX9loaPeYKsNbX5aWAjduACdPSv79eK1CAwNgxQpg5Uqgfft6NY1vmvHBHk7BUm0fnkVFRWjbtm2dDGqITJs2TdEmKBWqqldKSgp8fHwQHx/PBEZxcXGIj4/H+/fvmXrNmzfHTz/9xJoZ1rVr10rbbNasGbp27QoDAwP07du30sSMrq6usv8ySo6q3mMNEV75UiyWTPu/eBFwdweysirWGTsWWLoUmD4daNy4/m0EzzQDP+yRy/LZT58+RWJiIubPny+P5lWOoKAgdOvWTdFmKA3KpBchBFlZWUhKSkJiYiKzxcXFYePGjTAxMWHqhoSE1OpvpqCgAO/evWPynAHAqFGjUFhYiK5du6Jr167o0qULunbtCn19faipqcHFxQVz5syRy3dURZTpHqNUj8J9SQjw4oUkQHJxARISKtbp0AFYvFiydelS/zZ+hMI1+wg+2CNVsKSurl5l3pPPP/+cBku1pF27doo2Qangi15isRhFRUVo1qwZ6/jSpUsRHR3NBEbFxcWVnj9+/HhWsGRkZFShTpMmTWBkZIROnTrByMiI2T5e22zmzJmYOXNmlbbyRTNlgeqlOijEl2Ix8OwZcOWKZHv1qmIdbW1gyhRg4UJJpm0pkkjKGr7d/3ywR6pgady4cZg7dy7U1dWZYyKRCJ6entU+uClsGjWSSwefylIfegmFQqSkpCAlJQVJSUkVeoYSExORlJSE6dOnV3iNde/ePcTExNR4jbi4ONZ+586dYWlpyQqMynuGpIXeY9ygeqkO9ebL0lLAz08SHHl6Ah+k3WDQ0AAmTJBM+Z86FWjevH5s4wjf7n8+2COVBZaWlujbt2+F461bt67VjwVFQkxMDIYPH65oM5QGafQqLCxESkoKkpOTmWAoJSUF+/fvZwUl27Ztq1WesMpmT3To0AExMTFo2bIlBAIBOnTowNoEAgGMjIzQuXNn1nlNmzbFhg0b6vS9aoLeY9ygeqkOcvVlRgZw5w5w86ZkkPa7dxXrqKkBw4dLAqSZM4FaTJBSNHy7//lgj1TBUmWBEgB06tQJS5cuxf/93/9J03yDYeTIkYo2Qan4WC9CCPLz89GoUSM0bdqUOR4dHY1du3axgqMPl+H4kM2bN7MWazQ0NKzWBj09PXTo0KHSWZ+urq7Q0dFBcx79r5HeY9ygeqkOMvWlSAQEBkqCo5s3Ja/aKsu+o6UFjBsHfP+95FWbgYHsbKgH+Hb/88EeqYKlJUuWVDhWVFQEPz8/tK/nqY7KzMWLF7Ft2zZFm8FbHj9+jMTERKSlpSEtLQ23b9+GoaEhs5+WlgahUIhTp06x7snS0lJcuHChVtdISUlhBUv9+vXDrFmzYGhoiPbt27N6htq3bw9tbe0q26op0FIE9B7jBtVLdZDKl4QAb94A9+4BPj6AtzeQnV15XR0d4NtvJQHSN98Aurp1N1rB8O3+54M9nJJSfoyOjg4GDBjAGrOkpaUFY2NjbNmypcJrBmVF3kkpRSJRtUk+VQWxWIysrCxWkJOWlob09HTmc48ePWBpack6r3v37nj9+nWN7ZuammLHjh3Mfm5uLlq2bAlAMoW+PPAxNDSs8PmLL77gVU+QrGko95isoHqpDpx9GRsrCY7Kt+qe+Z99JgmMJk4EvvxSZZYe4dv9L0975JKU8mMOHz6MpUuXStMEBcD+/fuxc+dORZvBmZKSEmRmZiIjI4P5t/zzzp07ofnB6tempqbYvXs3RCJRtW1mfrwGEiQzISoLltTU1NCmTRu0bdsWbdu2hUAgYJXr6uoiMjIS7du3h46OTh2/pWqgrPeYoqB6qQ7V+pIQ4PVr4PFj4MEDSXAUG1t1Y7q6wNdfS4Kjb74BPnrmqAp8u//5YI9UPUsAEB4ejpMnTyI6OhqNGzfG8OHDsWTJEugqcRfkxzSEhXQJIXj//j0r4OnduzdrOvvr16+xcOFCpryq8T8AkJyczHoddejQIWzcuLFGO7p06VIhy/S5c+eQlpYGAwMDJjBq27Yt9PX1eTFLgkKhKAmFhZJxRo8fS5JDPnlSeWLIcrS1JYOzv/pKsg0cCHzwn0CK8lMvPUtnzpzBzz//DEIIOnToAAMDAxw6dAgHDx6Et7c3evbsKU3zDQZTU1OZRs2lpaV49+4dsrKyoK6uju7du7PKt23bhjdv3iArKwvZ2dnIyspCZmYmhEIhq56dnR1rVXkNDQ34+/vXyoaMjAxWsNSlSxcMGDCAFex8vBkYGLDGDZXz008/sfZlrVdDgGrGDaqXCiAWA69fw3PnTkwTCCSB0fPnQCXrIDJoaUlep5UHR4MHKyyLtiLh2/3PB3uk6lnq1KkTjI2N4eDgwMqu+eDBAxw+fBgeHh4yMVLRyLtnKSEhAR07dqxwXCwWIzc3F9nZ2RUWML1+/Tpu3rzJCnjKP+fl5TH1xo8fj9u3b7PO7du3L8LDw2u0a8+ePfj999+Z/by8PLRo0QItW7aEvr4+9PX1oaenx3z+cH/YsGFy612sSi9K1VDNuEH1UjLEYskCtM+e/bcFBwMFBdWf17q1JDgq3wYPlvQmNXD4dv/L05566VkSi8WwsrKqkIZ8xIgRuHnzpjRNQygUwt7eHs7OzhV6M+Li4rBlyxa0b98eGRkZOHjwICvDJyEEO3fuRHp6OgoLCzFnzhxMnjyZ1YaPjw+OHz8OAwMDaGpqwsLCQiGvdAIDA/Hbb7+hXbt2yM7Oxrt375Cdnc18Fv9vccXCwkLWDCx/f38cPXq0xvazKuli/rD3RlNTE23atGEFPeWfR48ezTpPR0cHJSUlrLFIiuDu3btYtGiRQm1QNqhm3KB68ZjiYuDlSyA0VLI9fw4EBQHVDAtg6NkTGDbsv+Do008leZAoLPh2//PBHqkHeD98+BD9+vWrUFZUVMTa9/LywnfffVerdsvKyuDs7AwHBwfWAqKAJGgYO3Ysjh07hvHjx+Py5cuYPHky/P39mdHyO3bsQGxsLC5evIjCwkL06NEDhoaGGDRoEADgxYsXmDdvHsLDw6Gnp4dff/0VmzdvhpWVVR1UkI63b9/Cx8enxnrZ2dmsAcwfv65SV1dHq1at0KZNG7Ru3Rpt2rRBmzZt8Omnn1Zo6/Tp01BXV0ebNm3QrFmzWmeIVlNTU3igBKDS70SpHqoZN6hePIAQSRbs8qCofIuMlOQ7qolOnYBBgxCnpwejGTOAQYMkPUmUGuHb/c8He6QKlry9vRESEoJHjx4xvR4ikQj//vsvgP/yMJWWlsLHxwdJlaV/r8yoRo2wePFivHr1Ci4uLqyyo0ePQigUYvz48QCAqVOnYsGCBbh48SJ++uknJCUlwdLSEtevXwcgyYo8adIkbN26FXfv3gUAbN++HePGjYPe/zKpzps3D8OHD8e6devqPd1B60r+eHV1dZlgpzzw+TigmTVrFr788kumvEWLFqwUDtXRhQcLNUrDx4E4pWaoZtygetUjhEim5798KdkiI4GICCAsrPKM2JXRvr0kGCrfBg5kEkG+uXsXRmPHyvELqB58u//5YI9UwVJeXh7S09Ohra3N+jEvn6ZdvuRJWVkZCmp6d1wJTZo0qXDs0qVLGDBgALOvoaEBExMTuLq64qeffsLVq1dRUlLCqjNw4EA4ODggPT0dTZo0wZ07d3DgwAGmvF+/fiCE4NKlS9i0aRNnO6Vh8ODB2LdvH2bNmoXWrVujRYsWtXodKBAIKkyVbyhU9mqRUj1UM25QveRAaakkwWN5QPRhcPTRG4Qq0dSUvEr77DP2Vk0iWOpL7vBNMz7YI1WwtHr1auzcuRO9e/eusa6Dg4M0lwIg6bUKCgqqMDvKwMAAAQEBAICAgADmldSH5YQQvHjxAlpaWigrK2O9xmrSpAl0dXXx/Pnzaq+fkpKCDh06sI5t2LBBqvW8WrRogRkzZlQY90WpGhMTE0WboHRQzbhB9aojJSWSPEVv3lTcYmKqn4n2Me3bVwyKunfnnPiR+pI7fNOMD/bU7r1NFQwfPrzKQOnjmXA///yzNJcCIBm3IxKJKry60tHRYZIZpqeno2XLlpX2dGVmZiI9PR1AxddfH7ZRFWKxmFmBvnzz8vJCcnIyTp06hcDAQNy+fRuXL19GREQErK2tkZeXB1NTUwCS6Y/Z2dk4evQoQkNDcfXqVVy/fh3Hjx+Hvb09MjIyWHULCwthaWmJV69ewc3NDXfv3sXjx4/h5OSEhIQEVl2RSIT9+/cjLi4O586dw8OHD+Hr6wsXFxdERUXBwsICQqGQdU5qaipOnDiBwMBA3Lp1C1euXEF4eDhsbGyQm5vLqpuTkwNbW1uEhYXB09MTN27cQFBQEBwcHJCWlsaqW1RUBEtLS0RFRcHV1RU+Pj549OgRzp49i/j4eOzbtw9isRimpqYQi8XYt28f4uPjcfbsWTx69Ag+Pj5wdXVFVFQULC0tUVRUxGq/fDxbUFAQbty4AU9PT4SFhcHW1hY5OTmsurm5ubCxsUF4eDiuXLmCW7duITAwECdOnEBqaiqrrlAohIWFBaKiouDi4gJfX188fPgQ586dQ1xcHPbv3w+RSMQ6JyEhAU5OTnj8+DHu3r0LNzc3vHr1CpaWligsLGTVzcjIgL29PYKDg3H9+nVcvXoVoaGhOHr0KLKzs1l18/LyYG1tjYiICFy+fBm3b99GYGAgTp06heTkZFbd4uJimJubIzo6Gs7OzvD19cWDBw9w/vx5xMbGwszMDH///TfrnMTERJw+fRr+/v7w9vaGu7s7IiMjYWVlhYKCAlbdzMxMHD9+HCEhIfDy8sK1a9cQEhICOzs7ZGVlserm5+fj8OHDiIiIgIeHB+7cuYOnT5/C0dGxgt2lpaU4cOAA3r59iwsXLsDPzw9+fn64cOEC3r59iwMHDqC0tJR1TnJyMhwdHfH06VPcuXMHHh4eiIiIwOHDh5Gfn8+qm5WVBTs7O4SEhODatWvw8vJCSEgIjh8/jszMTFbdgoICWFlZITIyEu7u7jh06BD8/f1x+vRpJCYmsuqWlZXBzMwMsbGxOH/+PB48eABfX184OzsjOjoa5ubmKC4urmC3NM+I4OBg3jwjIgICcGHbNrx3dob3xInAqlV4+8knEHXuDHGTJpKA5ttvgXXrABsbyRpqUVGVBkpETQ0FBgbI/OILJM2di8dLlyL18mVY7tgBcUICTIcNg3j/fuyLjUV8ixY46+JSq2dEWloa84ywsLCgz4haPCPKysqYc9avX8+rZ8SZM2fk9oz4eFx0VUiVOiApKQlmZmZ4/fo1SktLmeMikQgvXrxATk5OXZsGAOzevRtOTk6I/V9G1czMTOjr62PHjh3MFwaAH3/8EX///TcKCgrw3XffISQkhDUF0NvbG+PHj4e7uzsAyXgfb29vjBs3jqkjEAjQv39/eHl5VbCjfGqhurp6hXW/pO1ZAiSD1j9cAJZSPVQv7lDNuNGg9crNlfQOxcWx/y3/XJdXIk2bAp98AvToIXmNVr59+ilQyXALWdKgfVlH+KaZPO2pl9QB3333HbKysjB8+HA0/iBxFyGECXBkSZs2baClpVUhc3RBQQH09fUBSBYxvX//foVyAKyMz9W1URWGhoZyybN06NAhhSfcUiaoXtyhmnFDZfUqKQGSkyWzzBITJVtCAjsYqut/cnV1gW7dJEHRJ58Axsb/fW7XTmFT9FXWl3KEb5rxwR6pgqW3b98iLCwMnTp1qlB27do1aZquFDU1NfTt25d5lVZOamoqBg8eDEDybvPkyZOsSDQ1NRWNGjXCgAEDUFpaCk1NTVYbhYWFyMvLY9qobxR9EygbVC/uUM24oZR6vX/PDoISE9n7SUlAWlrd21dXBzp2BIyMgM6dJduHwVGbNrzMWaSUvlQwfNOMD/ZINWZp6dKlyMjIqLSsV69e0jRdJfPnz2clqSwrK8PLly8xc+ZMAMAPP/wATU1NPH36lKkTGhqK8ePHQ0dHB61bt8akSZNYbYSFhUFLS6tC4sr64sNXipSaoXpxh2rGDd7oRYjktVd4OPDPP8D588DBg8CmTcD8+cDYsUDv3kCrVkDz5pLxQmPHAgsXAjt2AMeOAVevSrJZ1xQoNWoEdO0qWeZj8WJg927AyQnw9ZUMzhYKJb1P9+8DZ84Af/4psWHIEEBPj5eBEsAjXyoRfNOMD/ZINWaprKwMu3btwoQJE1jHxWIxTp48ifPnz0tl3I4dO3DmzBnWq6/8/Hz07t0bnp6eGDBgANzc3HDo0CE8fvyYyTO0fv16iMViWFtbIz8/H3379sXly5eZdALPnj3DtGnT8Pr1azRt2hSrVq1C8+bNYW5uXqkd8l7uJCMjo8ZXgJT/oHpxh2rGDbnqJRZL8gdlZPy3paZKtpSU/z6npkoCnA/Gg9YZDQ3J9PoOHQCBQPLvh5uRkWT22f8S+6oS9N7nDt80k6c9tf19lypY+v3337Fv3z6oqanh42bU1NQgqk2W1Srw8PDA7t27ERERAUtLS8ybNw8G/0sy9vLlS/zxxx/o3LkzsrKyYGFhwUoFUFpaig0bNkBNTQ3v3r3DokWLMPajpGSenp5wdXWFnp4eWrRogT179lSZ1FHewZK9vT1WrFgh83ZVFaoXd6hm3OCkV1mZpPfnw+Cnui0rq3YZqGuLtrYkEOrYkR0Iffi5bVuVDIRqA733ucM3zeRpT70ES7q6uli7di1++ukn1rplIpEIVlZWsLGxqWvTvELewVJwcDAriSaleqhe3KGa1QJCgLw8IDsbLx8/Rk8DAyA7WxLcZGdX3DIzJcHPu3eSc2WJurokA3W7dpJAqF27ilv58ebNefsKjA/Qe587fNNMnvbUy2y4fv364f/+7/+YZUM+ZMeOHdI03aBISUlRtAlKBdWLOw1GM0KAwkLJjK7cXMm/5du7d9UHP9nZTI9PT3nYpq0N6OtXvn0cAOnpNdieIFnTYO59GcI3zfhgj1TB0rFjx3Dz5s0KGbUBwM/PD7Nnz5am+QaDNK8rGyJUL+4ojWYikaR358Ng58Ogp6bPubncskRLQ/Pm/wU7BgZVB0LlW7Nm9WMXhYXS3Ps8gm+a8cEeqYKlLVu24O3btzh9+jTruEgkQmhoKA2Wakl9L96r7FC9uCNXzUQioKAAyM+XbHl57H+r+lxVeX2jqytZjb58a9MGmYRAr1s31jFWnVatOC+7QVEM9HnBHb5pxgd7pAqW2rVrh4yMDAgEAtbiryKRCG/fvpXauIbCgwcP8NlnnynaDKWB6sUdRrOyMkk+nvfvJQHOx58rO/b+/X+BTGXBTm0XQZUXGhpAixZAy5b//VvZ5w+DnfIAqGVLyeKsH+F69ChWr15dr1+DIh/o84I7fNOMD/ZwGuBdntCxnODgYOjq6uKTTz6pUPfOnTsYP368bKxUMPIe4J2dnV1hrTpK1ai0XiIRUFQk2QoL//v3488cgx1xfj7UCwuB4mJFf0M2amqAjs5/W3XBTlXHmjWT+QBnlb7HGhjUl9zhm2bytEcuA7zd3d0xb948Zr+60enJyclcmm7QHDt2jBcZSpWFetdLJJIk5BMKKwYxlQUztSmr6lhJiVy+glTZZys0pi55dVUe4JR/ruxYVZ/L/23aVNIez6B/k6oD9SV3+KYZH+zh1LNkZGSE2bNns9IEVEZJSQlcXV0RExMjtYF8QN49S5QaEIslPSLlAcuHW1FR5cdr2ricV18DhmWNurqk16VZM8lg5I8/V3asuvLyAEdbm05Vp1AoKoFcepYSEhJgaWlZq7pq9GFaa0xNTRUeNQOQTLsuLZUEJiUlkn+r+1zbenU9vzxYkVNvi8LR1pb0rDRt+t/n2h6rTYDTpAkT1PDmHlMSqF6qA/Uld/imGR/s4dSz9Msvv+Dw4cNo3LhxtfUKCwuxbds2HDlyRGoD+YDcepbi44GcHBS8e4fmWlqSoKB8Ky2tfr+udWoKXBoqamqSQKRJk8q38rLGjf8LXOoS5JT/+0EgUx/k5eVBV1e33q6n7FC9VAfqS+7wTTN52iOXnqVffvmlxkAJAJo2bYqlS5dyabphsnw5cPs2mivaDj6gqSkJRBo3lkzJbty4ysDlVXw8upuYVB3Y1BTwVLY1aqTSr5acnJzw66+/KtoMpYHqpTpQX3KHb5rxwR5OwRKXqXuKnuanFCgiT0t5QPJhUFLbz/I6R0uLU6AiiogAevWSo0iqx7hx4xRtglJB9VIdqC+5wzfN+GCPVHmWKFIyfjzQrh2iExJg3LOnpHdFS4u9fXysLnXK9zU0VKL3JDIyEr1osMQJqhk3qF6qA/Uld/imGR/socGSIlmzBgDw5vZtGE+YoGBjlIdmdNkIzlDNuEH1Uh2oL7nDN834YA//Epw0QPiU/EsZoHpxh2rGDaqX6kB9yR2+acYHe2iwxANCQ0MVbYJSQfXiDtWMG1Qv1YH6kjt804wP9tBgiQdMnDhR0SYoFVQv7lDNuEH1Uh2oL7nDN834YA8NlhTMoUOHMH/+fBw6dEjRpigFVC/uUM24QfVSHagvucM3zfhiD6eklA0VeS53QpdS4QbViztUM25QvVQH6kvu8E0zedtT2/ZpzxKFQqFQKBRKNdBgiUKhUCgUCqUaaLBEoVAoFAqFUg10zFIt0NLSQmlpKdTV1WFoaCjTtlNSUiAWi+XStipC9eIO1YwbVC/VgfqSO3zTTN72lLevqamJkmoWk6fBUi3Q0NCAWCxWtBkUCoVCoVDkgLq6OkQiUZXldLmTWtCkSRMIhUJoaGjAwMBA0eZQKBQKhUKRAenp6RCJRGjSpEm19WjPEoVCoVAoFEo10AHeFAqFQqFQKNVAgyUKhUKhUCiUaqBjlmTAnTt3cOvWLXTr1g2jRo1Cr1698O7dO6xduxZ6enpITU3Fnj178Omnn7LOs7GxQWBgINTV1TF8+HD8/PPPrPIXL17A1NQU7du3R2FhIaysrNC8eXOmvLi4GOvWrUOjRo2QmZmJ9evXY8iQIfXynetCcXExOnfujNTUVACSlaQTEhLQtGlTxMXFYcuWLWjfvj0yMjJw8OBBtGvXjjmXEIKdO3ciPT0dhYWFmDNnDiZPnsxq38fHB8ePH4eBgQE0NTVhYWGBRo3+u8Vr4xM+kp2djR49esDNzQ2jR48GULvv0tDur9zcXKxcuRI3btyAnp4edu3ahQULFjDl9B5TLoRCIezt7eHs7Ax/f/8K5d9//z08PT2Z/SdPnjD3Z23uXTc3N1y5cgU6Ojro1KkTdu7cySqv6X5RBhT9zK2O+vh7lSmEUmfEYjHZsGEDWbZsGSktLWUdHzZsGLG3tyeEEBIUFEQ6depE8vLymDr29vbkyy+/JIQQIhKJSP/+/YmnpydTnpSURNq2bUsiIiIIIYRYWlqSadOmsa4/d+5csm3bNkIIIcnJycTAwIDExcXJ58vKAHt7e/Lbb78RKysrYmVlRW7cuEEIIeT9+/fE2NiY3L59mxBCiIeHBxk0aBApKytjzt22bRuZO3cuU79jx44kMDCQKQ8JCSFt27YlGRkZhBBC1q5dS9avX8+U18YnfGXp0qUEALl37x4hhN5fVbF+/Xri5eVFnj59SkaPHk3U1dWZ70fvMeWitLSUODo6kl69ehEjI6MK5ZGRkWTWrFnMs8TBwYFVXtO9e/PmTWJsbEyEQiEhhJDJkyeTw4cPM+W1uV+UAUU+c2tC3n+vsoYGS1Kwfft2MmbMmAp/QG5ubkRTU5MUFhYyx7p370727NlDCCGksLCQ6OnpkRMnTjDl+/fvJ8bGxkQsFhNCCFm1ahUZNmwYU56SkkIAED8/P0IIIQEBAQQAiYqKYupMmDCBLFmyRPZfVAaIRCIyZ86cSsvMzc2JQCBg9svKykizZs3I2bNnCSGEJCYmEi0tLeLt7c3UWbFiBRkzZgyzP2nSJPLjjz8y+0+ePCEaGhokJiaGEFKzT/iKn58fWbx4MStYovdXRfLy8sjr16+Z/ZiYGAKAuLu7E0LoPaasbN26tdJgac2aNSQtLa3Sc2pz7/bu3Zvs2LGD2Xd2dia6urqkoKCAEFLz/aIMKPqZWx318fcqa2iwVEf8/f2JpqZmpf/Tnj17Nunbty/r2Lx580jv3r0JIYTcuHGDACBBQUFM+Z07dwgAEhgYSMRiMWnbti1Zu3Ytq4327duT1atXE0II2bJlC9HR0WGVb9++nejo6BCRSCST7yhLLl26RLS1tcm4ceOIh4cHq2zw4MFk8uTJrGNffvkl+fbbbwkhhBw7dowAIFlZWUy5g4MDUVNTI2lpaSQ3N5c0atSIWFpaMuVFRUVEXV2dWFhYEEJq9gkfKS4uJtOmTSORkZGsYIneXzWTl5dHdHV1mf/10ntMOdm1a1eFYCk5OZloa2uTAQMGEEtLS6Z3qJya7t2IiAgCgPUcev36NevHuqb7RRlQ9DOXC7L+e5UHdIB3Hdm7dy8GDx4MJycnjBs3DsOGDYOPjw8AICAgAG3atGHVNzAwwMuXL1FcXIyAgAAAYNUpz9/0/PlzxMfHIy0trdI2nj9/Xu018vPzER0dLdsvKwMSExMxevRoBAQEYMaMGfjpp58gFoshEokQFBRU43dVV1dHq1atWOWEELx48QLPnz9HWVkZq40mTZpAV1e3Rr3KfcJHDhw4gPXr16Nx48as4/T+qpmLFy/i/Pnz0NPTo/eYihEWFoaJEyciPT0dGzduxJAhQ5CRkcGU13Tv1vT3UZv7RRlQ9DOXC7L+e5UHNFiqA+/fv4e3tzfat2+PdevW4ebNmxAIBJg0aRJiYmKQnp6O1q1bs87R0dGBWCxGdnY20tPTAYBVR0dHBwCQmZlZaXl5nczMTACo8hrlbfCNdevW4caNG0hOTsaqVatw/vx5HD58GNnZ2RCJRDV+15YtW0JNTY1VDkivV7lP+EZUVBRSU1MxatSoCmX0/qqamJgYbNq0CatXr4arqytKSkroPaZijB8/Hh4eHoiPj4eDgwMiIiKwdOlSpryme7emv4/a3C/KgKKfubVBXn+v8oDOhqsDb9++RUlJCX744Qe0aNECALBv3z64u7vDxsYGampq0NbWZp1TnkZdU1OTcfCHdWoqL6+jqakJADVeg680a9YMR48exbt373D8+HFm9oMsvqsq6fXHH3/g6NGjlZbR+6tq2rVrh7lz5yI2NhYXLlxAly5dsG7dOgD0HlM11NTUsGzZMhBCsGLFCiQlJUEgENTL34cyoahnbm2Q99+rLKE9S3WgoKAAAFjTrLt164aOHTvi9evXMDQ0RG5uboVzNDQ00Lp1a2YxwA/rlLepr69faXl5HX19fQCo8hrlbfCd7du3IzY2Fm3atIGWlpZU31Uavcp9wifOnDmDb7/9tkq76P1VNdra2hg4cCAuXbqE4cOHw8vLi95jKs7SpUthYGCAuLg4AHXz5YfltblflJH6fubWBnn9vcoDGizVgQ4dOgBAha71tm3bonXr1jAxMWG6KctJTU3FwIEDoa6uDhMTEwBg1SnPgzF48GAYGhrCwMCg0jYGDx4MAFVeo127dujUqZMMvqV86dq1Kzp06AA1NTX07du3xu9aUFCAwsJCVnmjRo0wYMAA9OzZE5qamqw2CgsLkZeXV6Ne5T7hE05OTvjpp5+gpqYGNTU1dOnSBQDw1VdfYfTo0fT+qiUzZsyAlpYWvcdUHA0NDRgZGTHP5Zru3Zr+Pmpzvygj9f3M5Yos/17lAf0LrgMdO3bEZ599hidPnrCOZ2dnY+zYsZg/fz5evHjBGtQZGhqKmTNnAgAmTJgAfX19VqK10NBQ9OrVCz179oS6ujrmzp3LKk9LS0N6ejpmzJgBAJg/fz7S09MRExPDamP69Oms97h85dmzZ1iyZAkAyXf58LuWlZXh5cuXjF4//PADNDU18fTpU6ZOaGgoxo8fDx0dHbRu3RqTJk1itREWFgYtLS0mSVlNPuETJ06cwPPnz5nt+vXrzPGTJ0/S+6uWCIVCjBgxAgC9x1SZ/Px8dOzYkQnia7p3+/TpAxMTkwp/H3p6ekzS15ruF2Wkvp+5XJHl36tckMscuwbAzZs3iZ6eHjPVMTAwkPTu3ZuUlJSQsrIyMnDgQGa65tOnT0mXLl2YHB6EEGJlZcVMjSwtLSX9+vUj165dY8rj4+OJnp4eSU5OJoQQcuDAATJr1iyWDVOnTmWmbsbHxxNDQ0OSkJAgvy9dR969e0emTZvG5MR48+YNWbRoESkuLiaESKaNduzYkZnq7urqSr744gvWFPV169YxU93z8vKIkZERa2p8YGAgEQgE5P3794QQQn755ReyefNmprw2PuEr5TlIylMH0PurIikpKeTgwYMkPT2dECKZXj506FCSmZlJCKH3mLKyfft2Vr4dQiS629rakrKyMpKXl0dWrlzJyqlESM337pUrV4iJiQnj/4kTJ5KjR48y5bW5X/gMH5651VEff6+yRo0QQuQThqk+bm5ucHV1Rb9+/RAfH49du3YxXcHJycnYsGEDOnXqhJSUFOzZs4d5nVLOn3/+ieTkZBQXF+Obb77BnDlzWOUPHz6EtbU1OnTogJKSElhYWLAGteXn52PNmjXQ09NDWloaNm7ciP79+8v/i3NEKBRiypQp8PPzQ69evTBp0iT8/vvvrCnxL1++xB9//IHOnTsjKysLFhYWrKmjpaWl2LBhA9TU1PDu3TssWrQIY8eOZV3H09MTrq6u0NPTQ4sWLbBnzx7W64/a+ISPxMbGokuXLrh37x7zP196f7EJCwvD1KlTkZ2djTFjxqBTp07YtGkT8/cI0HtM2fDw8MDu3bsREREBS0tLzJs3DwYGBti6dSuOHTsGfX19jBo1CqamphAIBKxza3Pv2tvb4/Hjx2jcuDF69+7NDCwup6b7hc/w5ZlbFfX19ypLaLBEoVAoFAqFUg10zBKFQqFQKBRKNdBgiUKhUCgUCqUaaLBEoVAoFAqFUg00WKJQKBQKhUKpBhosUSgUCoVCoVQDDZYoFAqFQqFQqoEGSxQKhUKhUCjVQIMlCoVCoVAolGqgwRKFQqFQKBRKNdBgiUJRAby8vDBjxgwsWrQIP/74IyZNmgR3d/caz8vLy4O1tTV69uwJX1/fKusFBASgdevWSExMlKHVikUoFEIgEMDDw0Om7SYlJWHnzp0wMDCo9+vHx8dj27Zt0NfXl0l7Xl5emDx5Mn7++WeZtEehKCuNFG0AhUKRjs2bN+Off/6Bt7c39PT0AAAZGRmYMGECHjx4AGtr6yrPFYvFMDAwQGRkZLXXMDQ0xDfffIOWLVuyzn358iV69+4tk+9R3zRu3BjffPMNunbtKtN21dTUkJ+fj4yMDNbx/Px8ZGVloXPnztVePywsDH379q3TtTU1NSEUCpGZmVmn8z9m8ODB+PXXX5VmTTQKRV7QYIlCUWJOnTqFgwcPIiQkhAmUAEBfXx+Ojo7o378/+vbti2XLllV6fsuWLfH555/XeJ2OHTvi4sWLrGNubm4QCoVKGyypqanh1KlTMm+3ffv2+OyzzyocP3r0KIYMGcIES5VdPzY2FkePHsXx48frdG1DQ8M6B1qVYWBgwFrclEJpqNDXcBSKkiIWi/H7779jwIABMDExqVDer18/fPbZZ9i1axdEIlGV7WhoaHC+dlhYGFauXMn5vIbCx5revXsXf/zxR7XnvHv3DjNmzIBQKJTq2rVZ9V2R7VEoygj9K6BQlJRnz54hJSUFgwYNqrLOgAEDkJKSgqCgIAQFBWHJkiWYOHEiXF1d0bp1a5iZmTF1CwoKMHfuXOjo6ODrr79GQkICAMmPuIWFBbp06YLY2Fjk5eXh1KlTyM3NxZkzZ7By5Urk5ORUuLZQKMTKlSuxdetWzJgxA+PGjUNiYiIIIbhx4wa+/fZb7NmzBzY2NmjTpg26deuGmzdvstpwcHDAunXrsHjxYgwfPhxPnjyRSfulpaW4ePEihg4dijNnzjBtmpqaYseOHZg3bx40NDQQHh5eqa5VXbsyEhMTcebMGZSWluLQoUNYvXp1pdc/fvw4UlJS8OTJE6xcuRIBAQFYsWIF1NTUmPFkLi4uaNGiBRYtWsS0X1hYiDVr1mDNmjWYN28eXFxcKthgb2+PDRs2YMSIERg7dmyNr12rIiUlBV27dsV3330HFxcXJCYmYseOHWjbti3evXuHH3/8Ebq6upg0aRKKiopw6tQpCAQCGBoawtvbu07XpFB4AaFQKEqJq6srAUC2b99eZZ0tW7YQAMTV1ZVERESQQYMGka5du5JTp06R7du3Ezc3NxITE0MAkKlTp5IrV64QR0dH0qJFCzJo0CAiFotJWloa2b9/PwFAYmJimLYBkNOnT1d57d27d5OhQ4cSQggRiUREIBCQzZs3k9LSUvLs2TPSrFkzMmrUKHL69Gny5MkT0rt3b9KsWTOSkpJCCCHEzs6OzJw5k2nPxsaGNG3alERFRUndfmFhIfH29mZ9h8ePH5Pvv/+eud769etJWFgYp+9WzunTp8mHj9dyje/du0cIIZVenxBCRo0aRRYuXMjsR0dHs84jhJARI0aw6nz//ffkwIEDhBBCxGIxmTBhAuvaTk5OxNPTk7F1/PjxpFu3bkQsFlf63T7mQ5vc3d3Jxo0bmbKEhASyatUqAoDY2tqS2NhY4uXlRQD8fzv3FtLkG8cB/KvNs2nmsZI0EmVm4UgpsIOaYAamF1EIpi7MCy2ptEi6UJeCUIJheVOQmRYdYHaRVmikdaHZSjPmIbOwSG1qUpuWh/3+F9LLpnOZ9k+t3+dqe9/n/T3P8+5iP54TxcTEkFwuJ5VKRbt27SJfX98Z1cfYQsRrlhhb5LRa7bT3zMzMAExMC4nFYojFYjQ3N+PAgQNCmXfv3gEAjhw5guDgYAAT62mkUikUCgUCAgJmtK5pMn9/f7i5uQnf3dzc0NfXB5FIhI0bN8LR0RGbNm0SRknOnTuHsLAwXLt2DampqcjMzERBQYHw/MGDByGTyZCbm4vLly/PKf6xY8cQEhKi197u7m5UVVWhuroaO3bsQEpKCiwtLX+pbzNlZWU1pX5DDE2B6V6rqamBXC5HSUkJgInfLTo6Gvfv3xfKyGQySKVStLS0AJhYfzY+Pg6VSvXTHXu6bt68iYaGBpw9e1a45u7uDolEAgBISUkBAHh4eMDFxQXe3t6Ijo4GAOzevRuHDh2acV2MLTScLDG2SK1ZswYA0NPTM22ZgYEBAMDatWsBTPzR2tvb/zR2REQEAKCzsxMBAQGzWtcUFRUFtVqNgoICDA4OQqPR6CV2JiYmsLKyEr5v374dS5YswZs3b6BUKvHp0ye9tlpYWMDPzw8NDQ1zjg9MXVe0c+dO+Pn5ISwsDOHh4cjLy4OXl9es+jYTs3mnk1VUVMDR0RG2trbCNd0Eb2hoCJ2dnUhMTNRL7n6VUqlEWVmZwak0kWjq34juewcmfruxsbFZ18/YfOM1S4wtUhKJBK6urnj69Om0ZZqamuDh4WFwd5YxS5cuBQDY2NjMun0vX74UEo+srCy4uroaLS8SieDg4AALCwsQEQCgt7dXr4ybm5swWjaX+IZYW1vj8ePHuHDhAp4/f47AwEDcu3fvt/Tt/6JWq/Hly5dpE7Xv378DABQKhd51jUaD4eHhGdfj6+uLyMhIxMbG/rZjCRhbTDhZYmyREolEyMjIgFKpRF1d3ZT77e3tqKurQ05Ozi/vaHr9+jWsra0RFBQ06/YlJydj69atEIvF05bR3aU3MjKCgYEBBAUFwcfHB/b29lMOyuzv70doaOic4xtSXV0NjUaD5ORkKJVKiMViFBUVzbpvukxMTGZVztzcHAD0EhutViskRz4+PhgdHUVNTc2UWFqtFg4ODlixYgVOnz6N0dFR4d6lS5dm3CbdZ0xNTZGQkCAks4z9KzhZYmwRS01NRXx8PBITE/UOQezt7UVMTAzS0tIQGxsrXNdqtcJow2Rv374FABAR8vLykJubKxxC+SPp0E0+bGxs0NbWhsbGRnR1dU2J9+3bNzx48ADt7e0oKSlBZ2cnent7UVVVJZRpbm4WPt+4cQN+fn6IioqCpaUlTpw4gdu3bwvTZj09PXj16hXS0tLmHN9Qn1QqlXCAp5OTE0JDQ+Ht7W3wXf2s7smxf4zQtbW1obq6Gmq1etp32tHRgf7+ftTW1sLV1RW2trYoKyuDUqlEfn4+Pnz4gPb2dnR0dGD//v1wdHREcnIyWltbMTQ0JLShpqYGGo0GJ0+eRH19PbZt24aioiIcPXoUg4ODsLS0xMOHD+Hp6TnlDC1d4+PjGB8fx/Lly1FaWorKykqcOXNG7/7kfugmdACE5MrYERaMLWjzu76cMTZXWq2Wrl69ShERESSVSikuLo4iIyPpzp07euUqKipo1apVZGZmRvn5+TQ2NkZEEzukiouLKSAggKKioighIYGuXLkiPNfV1UX79u0jAJSamkp9fX1ENLEjzM7OjjIyMgy26+7du+Ts7ExeXl5UXl5Ox48fp5UrV9KjR4+IiMjDw4O2bNlC6enplJ6eTnv27KH379/r9Ss7O5s2bNhASUlJFB8fT01NTb8l/vDwMMlkMgJAYWFh1NLSQtevXycAFB0dTadOnaKkpCRSq9W/3LfW1lYKDQ0lACSTyejz589ERJSYmEjLli2j8+fPG6z/R1wnJyeKiIgQ6r548SLZ2dnR+vXr6cWLFxQeHk5xcXH07NkzIiKqr68niURC5ubmFB4eTtnZ2RQYGEilpaU0MjJCWq2WZDIZubi4kJOTE6WlpdHo6CgREcnlcnJ2dqbNmzcb7GdxcTHZ2trS6tWrqbKykj5+/Eju7u5kampKmZmZ9OTJEwoJCSEAlJOTQyqVigoLC8nU1JQkEgnV1taSQqGg4OBgAkDZ2dmk0WgM1sXYQmZCxOOpjLE/z9PTEwkJCcjKylqU8f8mhw8fRmFh4Xw3g7EFi6fhGGPsH3br1i3s3bt3vpvB2ILGRwcwxubF2NiY3qLjxRb/b9DY2Ih169bB19d3vpvC2ILGI0uMsT/q69evyM/PR3d3N8rLyw3u5FrI8f8m/v7+nCgxNgO8ZokxxhhjzAgeWWKMMcYYM4KTJcYYY4wxIzhZYowxxhgzgpMlxhhjjDEjOFlijDHGGDOCkyXGGGOMMSM4WWKMMcYYM4KTJcYYY4wxI/4DLwmJ4eXg+aUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 650x325 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig = plt.figure()\n",
    "fig.set_size_inches([6.5, 3.25])\n",
    "\n",
    "plt.rc('font',family='Times New Roman')\n",
    "params = {'mathtext.default': 'regular' }\n",
    "plt.rcParams.update(params)\n",
    "\n",
    "plt.plot(apo_altitudes, mar_OI_DV, linestyle='solid', color='r', linewidth=2.0, label=\"Mars\")\n",
    "plt.plot(apo_altitudes, ven_OI_DV, linestyle='dashed', color='k',  linewidth=2.0, label=\"Venus\")\n",
    "\n",
    "plt.ylabel('Impulsive orbit insertion ' + r'$\\Delta V$' + ', m/s',fontsize=12)\n",
    "plt.xlabel('Orbit apoapis altitude, km',fontsize=12)\n",
    "plt.xticks(np.array([300, 5E3, 15E3, 30E3, 40E3, 50E3, 60E3]),fontsize=12)\n",
    "plt.yticks(fontsize=12)\n",
    "plt.grid('on',linestyle='-', linewidth=0.2)\n",
    "\n",
    "ax=plt.gca()\n",
    "ax.xaxis.set_tick_params(direction='in', which='both')\n",
    "ax.yaxis.set_tick_params(direction='in', which='both')\n",
    "ax.xaxis.set_tick_params(width=2, length=4)\n",
    "ax.yaxis.set_tick_params(width=2, length=4)\n",
    "ax.xaxis.set_tick_params(width=1, length=4, which='minor')\n",
    "ax.yaxis.set_tick_params(width=1, length=4, which='minor')\n",
    "ax.xaxis.grid(which='major', color='k', linestyle='dotted', linewidth=0.5)\n",
    "ax.xaxis.grid(which='minor', color='k', linestyle='dotted', linewidth=0.0)\n",
    "ax.yaxis.grid(which='major', color='k', linestyle='dotted', linewidth=0.5)\n",
    "ax.yaxis.grid(which='minor', color='k', linestyle='dotted', linewidth=0.0)\n",
    "ax.invert_xaxis()\n",
    "\n",
    "\n",
    "for axis in ['top', 'bottom', 'left', 'right']:\n",
    "    ax.spines[axis].set_linewidth(2)\n",
    "\n",
    "plt.legend(loc='upper left', fontsize=10, framealpha=0.8)\n",
    "\n",
    "plt.savefig('../../../data/mdpi-aerospace/smallsat-mission-concepts/orbit-insertion-DV.png', dpi= 300,bbox_inches='tight')\n",
    "plt.savefig('../../../data/mdpi-aerospace/smallsat-mission-concepts/orbit-insertion-DV.pdf', dpi=300,bbox_inches='tight')\n",
    "plt.savefig('../../../data/mdpi-aerospace/smallsat-mission-concepts/orbit-insertion-DV.eps', dpi=300,bbox_inches='tight')\n",
    "\n",
    "plt.show()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.16"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}