Skip to content

Instantly share code, notes, and snippets.

@d3v-null
Created July 16, 2024 07:24
Show Gist options
  • Save d3v-null/9262809b1f174ebac271411ba3f9cf51 to your computer and use it in GitHub Desktop.
Save d3v-null/9262809b1f174ebac271411ba3f9cf51 to your computer and use it in GitHub Desktop.
compare pyuvdata Van Vleck corrections with chebychev approximations, without, and raw
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"from pyuvdata import UVData\n",
"import os\n",
"from os import environ, path\n",
"\n",
"from time import time\n",
"\n",
"import cProfile\n",
"import re\n",
"\n",
"test_dir = environ['BIRLI_TEST_DIR']\n",
"obsid = 1090701368\n",
"vis_dir = path.join(test_dir, f'{obsid}/raw')\n",
"metafits = path.join(vis_dir, f'{obsid}.metafits')\n",
"raw_path = path.join(vis_dir, f'{obsid}_20140729203555_gpubox01_00.fits')\n",
"out_dir = '.'"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"some coarse channel files were not submitted\n"
]
}
],
"source": [
"uvd_raw = UVData()\n",
"uvd_raw.read_mwa_corr_fits(\n",
" [metafits, raw_path],\n",
" read_data=False,\n",
")\n",
"second_antpair = uvd_raw.get_antpairs()[1]\n",
"\n",
"nocorrect = {\n",
" 'use_aoflagger_flags': False,\n",
" 'remove_dig_gains': False,\n",
" 'remove_coarse_band': False,\n",
" 'correct_cable_len': False,\n",
" 'flag_small_auto_ants': False,\n",
" 'phase_to_pointing_center': False,\n",
" 'propagate_coarse_flags': False,\n",
" 'flag_init': False,\n",
" 'remove_flagged_ants': False,\n",
" 'fix_autos': False,\n",
"}\n",
"nocheck = {\n",
" 'run_check': False,\n",
" 'check_extra': False,\n",
" 'run_check_acceptability': False,\n",
" 'strict_uvw_antpos_check': False,\n",
" 'check_autos': False,\n",
"}\n",
"nofix = {\n",
" 'run_check': False,\n",
" 'check_extra': False,\n",
" 'run_check_acceptability': False,\n",
" 'strict_uvw_antpos_check': False,\n",
" 'check_autos': False,\n",
" 'fix_autos': False,\n",
"}\n",
"# uvd_raw.read_mwa_corr_fits(\n",
"# [metafits, raw_path],\n",
"# correct_van_vleck=False,\n",
"# **nocorrect\n",
"# )"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"print('profiling uvd_raw.read (no corrections)')\n",
"cProfile.run('uvd_raw.read_mwa_corr_fits( [metafits, raw_path], correct_van_vleck=False, **nocorrect )', 'read_mwa_corr_fits_raw.prof')\n",
"out_path_raw = path.join(out_dir, f'pyuvd_{obsid}_gpubox01_00.raw.uvh5')\n",
"uvd_raw.write_uvh5(out_path_raw, clobber=True, **nofix)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"approx 3.16s"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"print('profiling uvd_vvcheby.read (with cheby approx)')\n",
"cProfile.run('uvd_vvcheby.read_mwa_corr_fits([metafits, raw_path], correct_van_vleck=True, cheby_approx=True, **nocorrect)', 'read_mwa_corr_fits_vvcheby.prof')\n",
"out_path_vvcheby = path.join(out_dir, f'pyuvd_{obsid}_gpubox01_00.vvcheby.uvh5')\n",
"uvd_vvcheby.write_uvh5(out_path_vvcheby, clobber=True, **nofix)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"approx 11.2s"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"print('profiling uvd_vvnocheby.read (without cheby approx)')\n",
"cProfile.run('uvd_vvnocheby.read_mwa_corr_fits([metafits, raw_path], correct_van_vleck=True, cheby_approx=False, **nocorrect)', 'read_mwa_corr_fits_vvnocheby.prof')\n",
"out_path_vvnocheby = path.join(out_dir, f'pyuvd_{obsid}_gpubox01_00.vvnocheby.uvh5')\n",
"uvd_vvnocheby.write_uvh5(out_path_vvnocheby, clobber=True, **nofix)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"approx 2556.2s"
]
},
{
"cell_type": "code",
"execution_count": 27,
"metadata": {},
"outputs": [],
"source": [
"uvd_raw = UVData()\n",
"uvd_raw.read_uvh5(f'pyuvd_{obsid}_gpubox01_00.raw.uvh5')\n",
"uvd_vvcheby = UVData()\n",
"uvd_vvcheby.read_uvh5(f'pyuvd_{obsid}_gpubox01_00.vvcheby.uvh5')\n",
"uvd_vvnocheby = UVData()\n",
"uvd_vvnocheby.read_uvh5(f'pyuvd_{obsid}_gpubox01_00.vvnocheby.uvh5')"
]
},
{
"cell_type": "code",
"execution_count": 32,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0.5, 1.0, 'abs(vvcheby - vvnocheby) / abs(vvnocheby)')"
]
},
"execution_count": 32,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAiEAAAGzCAYAAAD5UcdSAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABltElEQVR4nO3deXxU1d0/8M/MJJPJNhPWhECAiMgimwYJQQSVaEC0RCki8hRKKVQlPGBcQYS6/bBYFNEIdSlolbL4VCzUxtLIYktACWBFgYIPmliYsJk9s97z+4MnI0MCzPfOJBPC5+1rXpjJ9849d+bO5Mw553u/BqWUAhEREVETM4a7AURERHR5YieEiIiIwoKdECIiIgoLdkKIiIgoLNgJISIiorBgJ4SIiIjCgp0QIiIiCgt2QoiIiCgs2AkhIiKisGAnJERWrlwJg8GAXbt2hewxS0pKYLFY8M9//jNkj0nB2bJlCwwGA95///2QPeavf/1rGAwGnDx5MiSPd+rUKcTGxuKjjz4KyePpwfdD4HhO0eWMnZBm7Omnn0Z6ejquv/76cDeFLiFt2rTBL3/5Szz55JPibdPS0vDAAw80QquCx/dD+ARzThFdCDshzdSJEyfw9ttv47777gt3U+gSdN9992H37t345JNPAt7m2LFj2LNnD0aPHt2ILdOH74fw03NOEV0MOyHN1LvvvouIiAjccccdF42tqalpghbRpaRXr17o06cPVq5cGfA2f/3rX2GxWHDzzTc3XsN0krwfqHHoOaeILoadkIv47rvv8MADD6BHjx6Ijo5GmzZtMG7cOHz77bcNxtfU1OBXv/oV2rRpA6vVikmTJuGHH37wi9m1axeysrLQtm1bREdHIzU1Fb/4xS/8YtavX4/09HTExcX53X/jjTeiT58+KCoqwrBhwxATE4O5c+cCAD788EOMHj0aycnJiIqKQrdu3fDMM8/A6/X6tl+6dClMJhPKysp89y1evBgGgwG5ubm++7xeL+Lj4/HYY4/pedrqycnJQVxcXIMdpgkTJiApKQlerxe33347rrjiigYfIyMjAwMHDvT9bDAYkJOTg/Xr16NPnz6IiorC1Vdfjfz8/Hrb7tmzB6NGjYLVakVcXBxGjBiBHTt21IsrKyvDgw8+iK5duyIqKgqdOnXCpEmT6s2ta5qG5557Dp06dYLFYsGIESNw+PDheo+3c+dOjBw5EjabDTExMRg+fPh51zScPHkSd999N6xWK9q0aYNZs2bB4XD4fj98+HD079+/wW179OiBrKwsv/tuueUWbNiwAYEWyv7LX/6Cm266CdHR0eeNaU7vB55Tzf+cIrooRRe0bt061b9/fzV//nz1+uuvq7lz56pWrVqpLl26qOrqal/cihUrFADVt29fdcMNN6ilS5eqGTNmKKPRqIYNG6Y0TVNKKVVaWqpatWqlrrrqKvXCCy+oN954Qz3xxBOqV69evsdyuVwqOjpa5ebm1mvP8OHDVVJSkmrXrp2aOXOm+t3vfqfWr1+vlFIqOztb3X333eqFF15Qy5YtU+PGjVMA1MMPP+zbfvfu3QqA2rBhg+++MWPGKKPRqAYOHOi77/PPP1cA1MaNG0PyPG7btk0BUGvXrvW7v7q6WsXGxqoZM2YopZR65513FAD12Wef+cV9++23CoB64YUXfPcBUP3791cdOnRQzzzzjFqyZIm64oorVExMjDp58qQvbt++fSo2NtYX9/zzz6vU1FQVFRWlduzY4YurrKxUffr0USaTSU2bNk0tW7ZMPfPMM+q6665Te/bsUUoptXnzZgVAXXPNNSotLU299NJL6te//rWKiYlRgwYN8mtzQUGBMpvNKiMjQy1evFi99NJLql+/fspsNqudO3f64hYsWOA7d+644w716quvqv/6r/9SANTPfvYzX9wbb7yhAKgvv/zSbz+fffaZAqDeeecdv/vffffdBuMb4nK5lNVqVa+++uoF45rT+4HnVPM+p4gCwU7IRdTU1NS7r7CwsN4btO5DNy0tTblcLt/9ixYtUgDUhx9+qJRS6oMPPlAA1Oeff37efR4+fFgBUK+88kq93w0fPlwBUMuXLw+orb/61a9UTEyMcjgcSimlvF6vslqt6tFHH1VKKaVpmmrTpo0aN26cMplMqrKyUiml1IsvvqiMRqP64YcfzttOCU3TVMeOHdXYsWP97l+7dq0CoLZt26aUUqq8vFxFRUWphx56yC9u0aJFymAwqO+++853HwBlNpvV4cOHffd98cUX9Z677OxsZTab1TfffOO77+jRoyo+Pl4NGzbMd9/8+fMVAPWnP/2pwfYr9eMfjF69eimn0+n7/csvv+z34axpmurevbvKysrybavUmdcoNTVV3XLLLb776v5g/OQnP/Hb5wMPPKAAqC+++EIppVRZWZmyWCzqscce84v77//+bxUbG6uqqqr87t++fbsCoNasWVPveM5VUFCgAKgjR45cMK45vR94Tp3RXM8pokBwOuYizh6adrvdOHXqFK688kokJCRg9+7d9eKnT5+OyMhI38/3338/IiIifKltCQkJAICNGzfC7XY3uM9Tp04BAFq1atXg76OiojBlypQLtrWyshInT57EDTfcgJqaGhw4cAAAYDQaMWTIEGzbtg0AsH//fpw6dQqPP/44lFIoLCwEAHz66afo06ePr73BMhgMGDduHD766CNUVVX57l+zZg06duyIoUOHAgCsVitGjRqFtWvX+g35rlmzBoMHD0bnzp39HjczMxPdunXz/dyvXz9YrVb87//+L4Az00p/+9vfkJ2d7Tck36FDB9x77734xz/+gYqKCgDA//zP/6B///648847G2z/2aZMmQKz2ez7+YYbbgAA33737t2LQ4cO4d5778WpU6dw8uRJnDx5EtXV1RgxYgS2bdsGTdP8HnPGjBl+P8+cORMAfOeOzWbDmDFj8Mc//tH33Hi9XqxZswbZ2dmIjY31277u/AkkTfOjjz5C79690bVr1wvGNaf3A8+p5n1OEQWCnZCLqK2txfz585GSkoKoqCi0bdsW7dq1Q1lZGcrLy+vFd+/e3e/nuLg4dOjQwTdnPnz4cIwdOxZPPfUU2rZtizFjxmDFihVwOp31HkudZ961Y8eOfh9Wdb766ivceeedsNlssFqtaNeuHf7rv/4LAPzaesMNN6CoqAi1tbX49NNP0aFDB1x77bXo378/Pv30UwDAP/7xD9+H4PmUl5fDbrf7bqdPn75g/Pjx41FbW4s///nPAICqqip89NFHGDdunN8H8vjx41FSUuLrEH3zzTcoKirC+PHj6z3muX9AgDMflHXrDk6cOIGamhr06NGjXlyvXr2gaRpKSkp8++nTp88Fj+F8+637cK7b76FDhwAAkydPRrt27fxub775JpxOZ73z59xzp1u3bjAajX7rLSZNmoTi4mLf6/T3v/8dpaWl+NnPflavjXXnz7l/7Bryl7/8JaCsmOb2fuA51XzPKaJARIS7Ac3dzJkzsWLFCsyePRsZGRmw2WwwGAy455576n3rCETdRYl27NiBDRs24OOPP8YvfvELLF68GDt27EBcXBzatGkDAPUW8NVpaOFgWVkZhg8fDqvViqeffhrdunWDxWLB7t278dhjj/m1dejQoXC73SgsLMSnn37q62zccMMN+PTTT3HgwAGcOHHiop2QWbNm4e233/b9PHz4cGzZsuW88YMHD0bXrl2xdu1a3HvvvdiwYQNqa2vr/SG44447EBMTg7Vr12LIkCFYu3YtjEYjxo0bV+8xTSZTg/s6XwcuVC6237rn+4UXXsCAAQMajD130fG5Gvqgz8rKQmJiIt59910MGzYM7777LpKSkpCZmVkvtu78adu27QX3c+TIERw4cADLli27YBzQ/N4PPKf8NZdziihQ7IRcxPvvv4/Jkydj8eLFvvscDodfdsnZDh06hJtuusn3c1VVFY4dO4bbbrvNL27w4MEYPHgwnnvuOaxatQoTJ07E6tWr8ctf/hKdO3dGdHQ0jhw5EnA7t2zZglOnTuFPf/oThg0b5ru/occYNGgQzGYzPv30U3z66ad45JFHAADDhg3DG2+8gYKCAt/PF/Loo4/6RlqA808fne3uu+/Gyy+/jIqKCqxZswZdu3bF4MGD/WJiY2Nx++23Y926dXjxxRexZs0a3HDDDUhOTr7o45+rXbt2iImJwcGDB+v97sCBAzAajUhJSQFw5lvivn37xPtoSN1wvtVqbfDDvCGHDh1Camqq7+fDhw9D0zS/KRKTyYR7770XK1euxG9+8xusX78e06ZNa/APWN1r36tXrwvu9y9/+QtsNptv+uJCmuP7gefU+YXrnCIKFKdjLsJkMtX7BvTKK6/4pb2e7fXXX/eb2162bBk8Hg9GjRoF4Mw3iXMfr+5bTd0QdGRkJAYOHCi65HXdB8bZj+1yufDaa6/Vi7VYLLjuuuvwxz/+EcXFxX4jIbW1tVi6dCm6deuGDh06XHCfvXv3RmZmpu+WlpZ20XaOHz8eTqcTb7/9NvLz83H33XefN+7o0aN488038cUXXzQ4bB4Ik8mEW2+9FR9++KHfEHRpaSlWrVqFoUOHwmq1AgDGjh2LL774Ah988EG9x5F+C05LS0O3bt3w29/+1m+9Qp0TJ07Uuy8vL8/v51deeQUAfOdOnZ/97Gf44Ycf8Ktf/QpVVVV+HcGzFRUVwWaz4eqrr75gWz/66CPceuutiIi4+HeS5vh+4Dl1RnM6p4gCxZGQi7j99tvxhz/8ATabDb1790ZhYSH+/ve/+4aIz+VyuTBixAjcfffdOHjwIF577TUMHToUP/nJTwAAb7/9Nl577TXceeed6NatGyorK/HGG2/AarX6fTscM2YMnnjiCVRUVPg+0C5kyJAhaNWqFSZPnoz//u//hsFgwB/+8IfzftDdcMMNeP7552Gz2dC3b18AQPv27dGjRw8cPHgQP//5z4XPVGCuvfZaXHnllXjiiSfgdDrP+4fgtttuQ3x8PB5++GGYTCaMHTtW9z6fffZZbNq0CUOHDsUDDzyAiIgI/O53v4PT6cSiRYt8cY888gjef/99jBs3Dr/4xS+QlpaG06dP489//jOWL19+3uspNMRoNOLNN9/EqFGjcPXVV2PKlCno2LEj/vOf/2Dz5s2wWq3YsGGD3zZHjhzBT37yE4wcORKFhYV49913ce+999bb7zXXXIM+ffpg3bp16NWrF6699toG27Bp0ybccccdF5y/r62txebNm7F8+fKAjqs5vh94TjWvc4pIpElzcS5BP/zwg5oyZYpq27atiouLU1lZWerAgQOqS5cuavLkyb64upTErVu3qunTp6tWrVqpuLg4NXHiRHXq1Clf3O7du9WECRNU586dVVRUlGrfvr26/fbb1a5du/z2W1paqiIiItQf/vAHv/uHDx+urr766gbb+s9//lMNHjxYRUdHq+TkZPXoo4+qjz/+WAFQmzdv9ov9y1/+ogCoUaNG+d3/y1/+UgFQb731lo5nKzBPPPGEAqCuvPLKC8ZNnDhRAVCZmZkN/h6A71oQZzv3tVHqzPOelZWl4uLiVExMjLrpppvU9u3b62176tQplZOTozp27KjMZrPq1KmTmjx5su8aEXXplOvWrfPb7siRIwqAWrFihd/9e/bsUXfddZdq06aNioqKUl26dFF33323Kigo8MXUpVN+/fXX6qc//amKj49XrVq1Ujk5Oaq2trbBY69Ldf1//+//Nfj7/fv3KwDq73//e4O/r7Nx40ZlMBhUaWnpBePqNLf3Qx2eU83nnCKSMCjFS981V1OnTsW///1v36p1ojovv/wyHnzwQXz77bcNZnPMnj0b27ZtQ1FR0QW/tT7wwAPYtWsXPvvss8Zsbkjw/dC4QnVOEUmwE9KMFRcX46qrrkJBQQErh5KPUgr9+/dHmzZtsHnz5nq/P3XqFLp06YK1a9fWWwB6rtdffx0pKSn11gg0R3w/NJ5QnlNEElwT0ox17tzZr84DXd6qq6vx5z//GZs3b8aXX36JDz/8sMG4Nm3aNLhwsSHTp08PZRMbFd8PodcY5xSRBEdCiC4R3377LVJTU5GQkIAHHngAzz33XLibRJc4nlMUbuyEEBERUVjwOiFEREQUFuyEEBERUVg02sLUvLw8vPDCC7Db7ejfvz9eeeUVDBo06KLbaZqGo0ePIj4+nmlgRER0QUopVFZWIjk5GUZj432vdjgccLlcQT+O2WyGxWIJQYtaiMa4+Mjq1auV2WxWv//979VXX32lpk2bphISEgK6IFJJSYkCwBtvvPHGG28B30pKShrjz5lSSqna2lqV1N4UknYmJSWd94Jxl6NGWZianp6O6667Dq+++iqAM6MbKSkpmDlzJh5//PELblteXo6EhAQMxW2IQGRgO9QxYmKKv3C1yXN5encV76PklvrVbi+k37BDovjHO+SL4gGga2SAz+n/8eo4PY563RcPOkucQbaPH7SGq41eiMXQcG2T89nvaifexw9arCj+uPvil+M/m1vJj/u0W3aeu3Q8t9Vesyj++6oE8T7cXlm7qpxR4n1IxUU5RfGVOtpUdUz2+iV8LXue2u2Rp90aD5WI4r2VOlJ7BZ87HrjxD3yEsrIy2Gw2+b4CUFFRAZvNhiNFXWCN1z/aUlGpITXtO5SXlwdUjuNyEPLpGJfLhaKiIsyZM8d3n9FoRGZmJgoLC+vFO51OX6EqAKisrPy/hkUiwtCInRCD7IMTEfLhM6NwyC0yVtamOB1vBmukbBs9nZAKr2wf8cJOiFuTH7dFuI8Yl/yPscMreztZ3LIOoUlHJ8Qs3Ac0+UeCyyM7byOU/I+xEnZCTKbG74RIPxL0tMkYLduJySx7niJMHlE8ABiFn52GQD/H/Qjer6puP40/fW+NNwbVCaH6Qv5snjx5El6vF4mJiX73JyYmwm6314tfuHAhbDab71ZXApuIiKg58Sot6Bv5C3uXbs6cOSgvL/fdSkpkQ31ERERNQYMK+qZHXl4eunbtCovFgvT09IvWelq3bh169uwJi8WCvn374qOPPvL7vVIK8+fPR4cOHRAdHY3MzEwcOuS/HOC5557DkCFDEBMTg4SEhHr7+OKLLzBhwgSkpKQgOjoavXr1wssvvyw+tpB3Qtq2bQuTyYTS0lK/+0tLS5GUlFQvPioqClar1e9GRETU3Ggh+E9qzZo1yM3NxYIFC7B79270798fWVlZOH78eIPx27dvx4QJEzB16lTs2bMH2dnZyM7Oxr59+3wxixYtwtKlS7F8+XLs3LkTsbGxyMrK8iuL4HK5MG7cONx///0N7qeoqAjt27fHu+++i6+++gpPPPEE5syZ41sLGqhGW5g6aNAgvPLKKwDOLEzt3LkzcnJyLrowtW4B0I0Y07hrQuLjRfGePqnifRSPjBHF977xsCj+oU4fi+IBYFCU7OUu9tSK9+EQrl0o02Rz5QlG2YJAADAKv4GUeOQL3ErcbUTxDqVnrlzmW0dbUXytcJEpIF+Y+oNT9r4AgJO1skW/HuG6IU01/nqCGqeO5/a47LgT9snW9CTurBTFA4DhwLeieE1PzRnJwlTlxhZ82KiLPev+Lh092CnohanJPb4XtVWa6DF+/HhUV1dj48aNvvsGDx6MAQMGYPny5VBKITk5GQ899BAefvhhAGcSQhITE7Fy5Urcc889fo+3cuVKzJ49G2VlZRdt64wZM7B//3588sknAR0b0EjTMbm5uXjjjTfw9ttvY//+/bj//vtRXV2NKVOmNMbuiIiIGp1XqaBvwJlOzdm3s5MzzlaX6JGZmem770KJHgBQWFjoFw8AWVlZvvgjR47Abrf7xdhsNqSnp5/3MQNVXl6O1q1bi7ZplIuVjR8/HidOnMD8+fNht9sxYMAA5Ofn11usSkREdKkIZl1H3fYA6iVgLFiwAL/+9a/rxV8o0ePAgQMN7sNut18wMaTu30CTRwK1fft2rFmzBn/5y19E2zXaFVNzcnKQk5PTWA9PRER0SSopKfGbjomKavyU8sa0b98+jBkzBgsWLMCtt94q2rbROiFEREQtiQYFbwhGQgJNwpAmegBAUlLSBePr/i0tLUWHDh38YgYMGBDwsdT5+uuvMWLECEyfPh3z5s0Tbx/2FF0iIqJLQVOn6JrNZqSlpaGgoODHNmgaCgoKkJGR0eA2GRkZfvEAsGnTJl98amoqkpKS/GIqKiqwc+fO8z7m+Xz11Ve46aabMHnyZDz33HOibetwJISIiKiZys3NxeTJkzFw4EAMGjQIS5Ys8Uv0mDRpEjp27IiFCxcCAGbNmoXhw4dj8eLFGD16NFavXo1du3bh9ddfB3DmyrKzZ8/Gs88+i+7duyM1NRVPPvkkkpOTkZ2d7dtvcXExTp8+jeLiYni9XuzduxcAcOWVVyIuLg779u3DzTffjKysLOTm5vrWk5hMJrRrF3jZC3ZCGpHR3bhpf3pSVas0WU/cZpQfQ2thbz/eWCOKd+sYDf1WmHJrMsjz+aXb1AhTW+ONjosHncMpvAx7lFFW9wcANMjOkRNKPgDrFabQSlNuvZqONH9jyK9uUI/BK2uXONNY4xU8Jc7OcNG7vdTFEj2Ki4v9qgcPGTIEq1atwrx58zB37lx0794d69evR58+fXwxjz76KKqrqzF9+nSUlZVh6NChyM/P96vuO3/+fLz99tu+n6+55hoAwObNm3HjjTfi/fffx4kTJ/Duu+/i3Xff9cV16dIF3377bcDH1yjXCQlGS7pOyPcjZDn+PW+RFbB7pvOfRfEAkGySvdx65j9Nwj9KlcJLGTdFJ0QP6XVCTnhk56CeTsjXNcmi+EhhoT8AqNVknani6lbifZyqlV1bxCu8TkhTdEL0XCekRljAznpAdo2epO3longAMPy7WBTfkq4TcmB/IuKDuE5IZaWGnr1KWcDuLFwTQkRERGHB6RgiIqIAeIPMjglm25aKnRAiIqIAeNWZWzDbkz92QoiIiAKg/d8tmO3JH9eEEBERUVhctiMh4qQgHRk4wuxIRBhl/WRptVoAMAqzHxw6UvhijLLqsHavRxQvTQkFgHamatk+dFRVtRsSRPFu4QlS6pVn+GjCdNgq6UkLoNJtuXjQWdxePeet7P1qMbtE8Q6P/LiV8BzRk4eohNlsRj2pY0IGHZ+FLYUGA7w6Pn/O3p78XbadECIiIglNnbkFsz3543QMERERhQVHQoiIiALgDXI6JphtWyp2QoiIiALATkjocTqGiIiIwoIjIURERAHQlEFX5tzZ25O/y7YT0iRpZsK0QukJGqmj0muNJi9QJvWdR5Zye1qTpXjGG2XplwDgFqaqunWkP0tZhBVrTcJCfwAQGyGrtHzaJSu6CAA1HnlhtsYmTbmVptsCgNPdBB+fzKZoVjgdE3qcjiEiIqKwuGxHQoiIiCS8MMIbxHf3xh+HvvSwE0JERBQAFeSaED3Tfi0dOyFEREQB4JqQ0OOaECIiIgoLjoQQEREFwKuM8Aoz7fy3D2FjWgh2QgKlpwSmkMMrrD7riRfvo0d0jSi+UpOnw5qEqcmtjQ5RfLkWJYoHAKMwnfm4V/7cOjTZ63faI0uHtUXIXjsAqPTI0p+dmjw12RIhSzX26PgQ9zpl25iEFakNwnMWALSIxl9mWOuVDd+Ls/Yv44q4emgwQAtiAkFjznU9nI4hIiKisOBICBERUQC4MDX02AkhIiIKQPBrQjgdcy5OxxAREVFYcCSEiIgoAGcWpgZRwI7TMfWwE0JERBQALcjLtjM7pr7LthOipHNzelLZhJfoTTDXiuIdSpYSCgDlwpRbl645TNlxS78dSNNtAX3PlZS08m7byMpGasmPkqPKRPEH3EnifUgvY+30Nv7HjqMJKtxGR8qqRVfVylPLlUlYiTuC37Tp0nLZdkKIiIgkuDA19NgJISIiCoAGIy9WFmLshBAREQXAqwzwBlEJN5htWyqm6BIREVFYcCSEiIgoAN4gs2O8nI6ph50QIiKiAGjKCC2IhakaF6bWc9l2QgzNsHqkntRTqRNe2RvIZpRXCpXOezqEqa160m2lFW5NOl6LeJOsGjCET+1JHVWTj7tk20QZZWmngDxF16jj26BXk523mjA+UkdF3AqHLOU2SpjSCwAOrZE/p/hHkcLssu2EEBERSXA6JvTYCSEiIgqAhuAyXBp/rPvSw+wYIiIiCguOhBAREQUg+IuV8Xv/udgJISIiCkDwl21nJ+RcfEaIiIgoLDgSEiCDt/GXFDm8sjTSzhE/iPeRaJIdR6lXlj4LADEGabqjLHXRoeSnrdUoS5/91t1WvI9Kr0UU/72rtSjeYnSL4gGgVWSNKN5okK/el1ZBbgrSdFi3MHUdAJRwgaJbmDash0Gacmvk91AJDYagzvfm+F4JN3ZCiIiIAsDpmNBjJ4SIiCgAwV8nhJ2Qc/EZISIiorDgSAgREVEANGUQlyk4d3vyx04IERFRALQgp2N4nZD6+IwQERFRWHAkJEBKR9VdaSFWaYqu3WuV7QBAp4jTovh4HVV0HcIhx+88rUTx3SNPieIB4D/eOFF8gkmW2qqHI1JeDVjqW4cs1fikU/Y8AUCVW1ZN1tMEGQIGYaqxWUcVXY8wfd3lkae7K1PzK3imLuPKu5oyQgvi/A1m25aKnRAiIqIAeGGAN4hrfQSzbUvFbhkREVEzlpeXh65du8JisSA9PR2fffbZBePXrVuHnj17wmKxoG/fvvjoo4/8fq+Uwvz589GhQwdER0cjMzMThw4d8ot57rnnMGTIEMTExCAhIaHB/RQXF2P06NGIiYlB+/bt8cgjj8DjkV0okJ0QIiKiANRNxwRzk1qzZg1yc3OxYMEC7N69G/3790dWVhaOHz/eYPz27dsxYcIETJ06FXv27EF2djays7Oxb98+X8yiRYuwdOlSLF++HDt37kRsbCyysrLgcPx4ZWmXy4Vx48bh/vvvb3A/Xq8Xo0ePhsvlwvbt2/H2229j5cqVmD9/vuj42AkhIiIKgBc/Tsnou8m9+OKLmDZtGqZMmYLevXtj+fLliImJwe9///sG419++WWMHDkSjzzyCHr16oVnnnkG1157LV599VUAZ0ZBlixZgnnz5mHMmDHo168f3nnnHRw9ehTr16/3Pc5TTz2FBx98EH379m1wP3/729/w9ddf491338WAAQMwatQoPPPMM8jLy4PL5Qr4+NgJISIiakIVFRV+N6fT2WCcy+VCUVERMjMzffcZjUZkZmaisLCwwW0KCwv94gEgKyvLF3/kyBHY7Xa/GJvNhvT09PM+5vn207dvXyQmJvrtp6KiAl999VXAj8NOCBERUQBCNR2TkpICm83muy1cuLDB/Z08eRJer9fvDz0AJCYmwm63N7iN3W6/YHzdv5LHlOzn7H0EQpwds23bNrzwwgsoKirCsWPH8MEHHyA7O9v3e6UUFixYgDfeeANlZWW4/vrrsWzZMnTv3l26q0bVFGlmzTEbK8YgSwv9X498NXcbY8O9+vOJN9aK4qt1VNF1K2E6pTAeABxK9tzWeGWprU4dx10rTPuOMsoWlQGA0SzLRa/2mMX7iDDJBrJrnLJ9SKvuAvrSeqVktZ91XEpAa/zq4C1JqArYlZSUwGr98RILUVGyz4KWRPxsVldXo3///sjLy2vw94EseCEiIrrUKBigBXFT/5eia7Va/W7n64S0bdsWJpMJpaWlfveXlpYiKSmpwW2SkpIuGF/3r+QxJfs5ex+BEHdCRo0ahWeffRZ33nlnvd8FuuCFiIiILsxsNiMtLQ0FBQW++zRNQ0FBATIyMhrcJiMjwy8eADZt2uSLT01NRVJSkl9MRUUFdu7ced7HPN9+vvzyS78snU2bNsFqtaJ3794BP05IL1Z2sQUv99xzT71tnE6n36KcioqKUDaJiIgoJEI1HSORm5uLyZMnY+DAgRg0aBCWLFmC6upqTJkyBQAwadIkdOzY0beuZNasWRg+fDgWL16M0aNHY/Xq1di1axdef/11AIDBYMDs2bPx7LPPonv37khNTcWTTz6J5ORkv6UVxcXFOH36NIqLi+H1erF3714AwJVXXom4uDjceuut6N27N372s59h0aJFsNvtmDdvHmbMmCGaXgppJ0TPgpeFCxfiqaeeCmUziIiIQi4cVXTHjx+PEydOYP78+bDb7RgwYADy8/N9f2eLi4thNP7YuRkyZAhWrVqFefPmYe7cuejevTvWr1+PPn36+GIeffRRVFdXY/r06SgrK8PQoUORn58Pi8Xii5k/fz7efvtt38/XXHMNAGDz5s248cYbYTKZsHHjRtx///3IyMhAbGwsJk+ejKefflp0fGG/bPucOXOQm5vr+7miogIpKSlhbBEREVHzkZOTg5ycnAZ/t2XLlnr3jRs3DuPGjTvv4xkMBjz99NMX7DCsXLkSK1euvGC7unTpUu9qrFIh7YScveClQ4cOvvtLS0sxYMCABreJioq6rFcGExHRpcELI7xBXNkimG1bqpB2Qs5e8FLX6ahb8HK+S79eMpqg7pDF5JbFG2TxAPC9V7ZNgo73jHSTBGFK7wlvrHAPQBtj41fFLXG3EcW3jqgSxR9zJ4jiAcCjyVKNnZr8I0FaRTdSR2VmKaOwiq5bR4Vb6dC6V5O/mQzCFHmDJrz0gJ7q4Dq2aSnCMR3T0ok/caqqqnD48GHfz0eOHMHevXvRunVrdO7cOaAFL0RERETiTsiuXbtw0003+X6uW88xefJkrFy5MqAFL0RERJcaDUZoQUypBLNtSyXuhNx4440XvNpoIAteiIiILjVeZYA3iCmVYLZtqdgtIyIiorAIe4ouERHRpYALU0OPnRAiIqIAqLMq4erdnvxdtp2QpkgzE2YJokZYXbS1SZ52ahE2ymaUpy6e9MpSMKVVcWOEKb0A4BDuI0FY2RcAIg2y4z7pjhfF67nkc7Qw7TtCR/psmStaFC9NnwUAk3Ab6TdOj1f+3FrMsudW17fgxv6YaoJq4i2JFwZ4g3hRgtm2pWK3jIiIiMLish0JISIiktBUcOs6pNeSuxywE0JERBQALcg1IcFs21LxGSEiIqKw4EgIERFRADQYoAWxuDSYbVsqdkKIiIgCwCumhh47IY1ImLEpTl087Y2R7QBArKFCFP+9R7wL/MfTShTfMULWpv94rKJ4ALAYZemU37raivchFWOSpRr/4JG/3tKU2+O1cfJ9GDRRvEtY2ReQvzciTLLjNkfIU5OliwwjjLLnCQAgLYorXfiop0lM66UQYieEiIgoAFyYGnrshBAREQVAQ5CXbeeakHrYLSMiIqKw4EgIERFRAFSQ2TGKIyH1sBNCREQUAFbRDT12QoiIiALAhamhd9l2QsRpZjqy0rRI2UYxES5RvJ6KjNJtYqR5xgDamKpF8XpSbqWkb35pRVwAiBVW9z3ulh13qwh51eRT7lhRfKzwHAQAjyZ7bmvcsqq7gDxFVwm/cTp1VNGV7qMpUnS1COFnAv8mUphdtp0QIiIiCU7HhB47IURERAHgZdtDj4NxREREFBYcCSEiIgoAp2NCj50QIiKiALATEnqcjiEiIqKwuGxHQgyGxu+RGl2yfdR4zKJ4r46c80STMJ1Sk6eqmoR5he1NVaJ4h5KfttVK9tyahJVhAaDaGyWKjzHK0mFLhSm9AODUIkXxsSZ5iq47QlYV12WWpTIDwKlaeQVhiUhh1V0AMBll53lVrez8ACD+mmiQlvbVkTV8OeNISOhdtp0QIiIiCXZCQo/TMURERBQWHAkhIiIKgEJw1/rQceHtFo+dECIiogBwOib02AkhIiIKADshocc1IURERBQWl+1IiLiKrh7CTq/F5BbFW40O2Q4AVApTbt06nqYKJUtFdCtZiqee1GSzsCpumbdxU0IBeUVjabVaACgTVqyt8chSegF5hVs93wad7sb9qNJ0PLfS/FaTSUc+rDxzWEbHYTfF5Q2aK46EhN5l2wkhIiKSYCck9DgdQ0RERGHBkRAiIqIAKGWACmI0I5htWyp2QoiIiAKgwRDUdUKC2bal4nQMERERhQVHQoiIiALAhamhx05IgAx6UnqFGXlmoywfT1oZFgBqlCyt1yJMvwQAq0FWJfWUavx02BijrE0WgyxdGgAcwoq1bk329osWpnADQGJUhSj+hCFOvA+PJkuxLhfGA0CEML3V45UN8kZGyHNha12y19tRK3+/QvpUSd+ul3G6rR5cExJ6nI4hIiKisOBICBERUQA4HRN67IQQEREFgNMxocdOCBERUQBUkCMh7ITUxzUhREREFBbshBAREQVAAVAqiJvO/ebl5aFr166wWCxIT0/HZ599dsH4devWoWfPnrBYLOjbty8++ugj/+NQCvPnz0eHDh0QHR2NzMxMHDp0yC/m9OnTmDhxIqxWKxISEjB16lRUVVX5xXz88ccYPHgw4uPj0a5dO4wdOxbffvut6Ng4HRMgpSOVTQmfXZcwddGh5BVPY4Qpt6eFaaR6tDHWiOLtXqt4H5WaRRTvElb2BYCkyHJR/ElPvCjeZJBXYXUK04ZNOlKyT7pklXr1DGdLK/WajLJ46eMD8qH1KIs8xbpGk1Wk1nEYYk1SgbyZ0mCAoYmvmLpmzRrk5uZi+fLlSE9Px5IlS5CVlYWDBw+iffv29eK3b9+OCRMmYOHChbj99tuxatUqZGdnY/fu3ejTpw8AYNGiRVi6dCnefvttpKam4sknn0RWVha+/vprWCxnPisnTpyIY8eOYdOmTXC73ZgyZQqmT5+OVatWAQCOHDmCMWPGIDc3F++99x7Ky8vx4IMP4q677sLu3bsDPj6OhBARETVTL774IqZNm4YpU6agd+/eWL58OWJiYvD73/++wfiXX34ZI0eOxCOPPIJevXrhmWeewbXXXotXX30VwJlO5JIlSzBv3jyMGTMG/fr1wzvvvIOjR49i/fr1AID9+/cjPz8fb775JtLT0zF06FC88sorWL16NY4ePQoAKCoqgtfrxbPPPotu3brh2muvxcMPP4y9e/fC7Q68w81OCBERUQDqsmOCuQFARUWF383pbPhiii6XC0VFRcjMzPTdZzQakZmZicLCwga3KSws9IsHgKysLF/8kSNHYLfb/WJsNhvS09N9MYWFhUhISMDAgQN9MZmZmTAajdi5cycAIC0tDUajEStWrIDX60V5eTn+8Ic/IDMzE5GRgY/AshNCREQUgLrrhARzA4CUlBTYbDbfbeHChQ3u7+TJk/B6vUhMTPS7PzExEXa7vcFt7Hb7BePr/r1YzLlTPREREWjdurUvJjU1FX/7298wd+5cREVFISEhAd9//z3Wrl170efxbOyEEBERNaGSkhKUl5f7bnPmzAl3k8TsdjumTZuGyZMn4/PPP8fWrVthNpvx05/+VLRuiAtTiYiIAlCX5RLM9gBgtVphtV58gX3btm1hMplQWlrqd39paSmSkpIa3CYpKemC8XX/lpaWokOHDn4xAwYM8MUcP37c7zE8Hg9Onz7t2z4vLw82mw2LFi3yxbz77rtISUnBzp07MXjw4IseH8CRECIiooCEak1IoMxmM9LS0lBQUOC7T9M0FBQUICMjo8FtMjIy/OIBYNOmTb741NRUJCUl+cVUVFRg586dvpiMjAyUlZWhqKjIF/PJJ59A0zSkp6cDAGpqamA0+nchTCaTr42B4khIoJrgQnceTVj5E/LKn+YmqJrpEvZt44UVa43S8sQAvErWpjYRVRcPOofdnSCKr/HqqKrayLw60mel1Z9r3I1/3F5NdhwGHe+LqEjZeev2ytO+pcQv32WcbnupyM3NxeTJkzFw4EAMGjQIS5YsQXV1NaZMmQIAmDRpEjp27OhbVzJr1iwMHz4cixcvxujRo7F69Wrs2rULr7/+OoAz5/rs2bPx7LPPonv37r4U3eTkZGRnZwMAevXqhZEjR2LatGlYvnw53G43cnJycM899yA5ORkAMHr0aLz00kt4+umnMWHCBFRWVmLu3Lno0qULrrnmmoCPj50QIiKiAISjdsz48eNx4sQJzJ8/H3a7HQMGDEB+fr5vYWlxcbHfiMSQIUOwatUqzJs3D3PnzkX37t2xfv163zVCAODRRx9FdXU1pk+fjrKyMgwdOhT5+fm+a4QAwHvvvYecnByMGDECRqMRY8eOxdKlS32/v/nmm7Fq1SosWrQIixYtQkxMDDIyMpCfn4/o6MCvHWRQzezKMxUVFbDZbLgRYxBhCDDNR8e3GGNcnCje2/cK8T5KbokVxfcacejiQWe5L3mLKB4ABkSVieLtOr69OYQX+pKOhJR4bKJ4PdyQH7d0JOQbR/0LDYVahUd2IbEKj+ziWABQIbxYWZlTFg8ADo/s+5LLI3v9Ikw6RteEoy16RkKqSmQX5kv4Sjbi1/7zClE8ABj+XSyK16rko4qSERqPcmMLPkR5eXlA6yz0qPu71GPV4zDFyN8jdbw1Thy89/lGbeulhiMhREREAQjVwlT6ERemEhERUVhwJISIiCgAZ0ZCglkTEsLGtBDshBAREQUgHAtTWzpRJ2ThwoX405/+hAMHDiA6OhpDhgzBb37zG/To0cMX43A48NBDD2H16tVwOp3IysrCa6+9Vu8SseEmTckzeOQL16RFTz06KrdKtTfJFsuWadXifUgXprqFs4JWo0MUDwCnNNlxn/LIFi4D8qrGMSaXKD7SIE/JLq5tLYo36qjUW+mWLdQzGeX7kG7jFqY/61mY6hEuNNWEKfgAYPAIP6ek37Tlh60rnZnofETviq1bt2LGjBnYsWOHr7zvrbfeiurqH/9QPfjgg9iwYQPWrVuHrVu34ujRo7jrrrtC3nAiIqKmpEJwI3+ikZD8/Hy/n1euXIn27dujqKgIw4YNQ3l5Od566y2sWrUKN998MwBgxYoV6NWrF3bs2BHwZVyJiIiaG07HhF5Q2THl5eUAgNatzwz5FhUVwe12+5UI7tmzJzp37nzessNOp7NeWWMiIiJq+XR3QjRNw+zZs3H99df7rsRmt9thNpuRkJDgF3uhssMLFy70K2mckpKit0lERESNh/MxIae7EzJjxgzs27cPq1evDqoBc+bM8StpXFJSEtTjERERNYpgi9dxOqYeXSm6OTk52LhxI7Zt24ZOnTr57k9KSoLL5UJZWZnfaMiFyg5HRUUhKkr/ZXCJiIiaAq+YGnqiTohSCjNnzsQHH3yALVu2IDU11e/3aWlpiIyMREFBAcaOHQsAOHjwIIqLi89bdjhcpCVzVIR80EiacWsUjtWZdaRsfu+R1XFIMMp77gdcstouXSNPi+IrNMvFg84RL0zrrTbKO8YnPbJaENLX+7RHlmYMAJHCCrfVHnmFW5u5VhRfLqw1AwBeYXprbJQs/dnpln8fixHu4+QP8eJ9qAjh55Q0fdYkf383s3JjdIkTvfNmzJiBVatW4cMPP0R8fLxvnYfNZkN0dDRsNhumTp2K3NxctG7dGlarFTNnzkRGRgYzY4iI6JLG7JjQE3VCli1bBgC48cYb/e5fsWIFfv7znwMAXnrpJV/Z37MvVkZERHRJC3ZdBzsh9YinYy7GYrEgLy8PeXl5uhtFRERELR9rxxAREQWAC1NDj50QIiKiQAR7rQ92QuoJ6oqpRERERHpxJCRQOipHGjyN0I6zxBidjbsDACe88n5qx4gyUXyZMOW2jUle2dfulaXPupX8rdE2QlZy4EBtsngfUtEmtyheT4quwyurHhwfKT9vHR7Z61HlaPzKvtVO2XNlEJe4bYIqupwfEGF2TOixE0JERBQo9ttCitMxREREFBYcCSEiIgoAp2NCj50QIiKiQDA7JuTYCSEiIgqI4f9uwWxPZ+OaECIiIgqLy3YkxCBNudV0pNcJi9y6NFnZ3RNeeVXOZJMsPTLSIE9d9Ap7+1bI2uTS0Xe2GGSpql4dc7fStF5bhKz6rJ4qutJKvRE6Xu+2Flll5m/K24r3IRVhkr35dBSLhlf4VBlN8ufWI6yiK65NIm+S/LOzJeF0TMhdtp0QIiIiEXZCQo7TMURERBQWHAkhIiIKhDLIp7zO3Z78sBNCREQUAFbRDT1OxxAREVFYcCSEiIgoEFyYGnKXbSdECcfFDDrG0aSFWDXhfKFJxxldqckGv8w6UjYbm0PJqrYCQLzRIYrXU0W3RpNVbq30yqoHx5tkxwAAp12ytF6Pkg+O1rhk1WQtEfLy0hUO2XMlvTy2R5PP1Uvfr7o09li1SX4M0s/OFoVrQkKO0zFEREQUFpftSAgREZGEQZ25BbM9+WMnhIiIKBBcExJy7IQQEREFgmtCQo5rQoiIiCgsOBJCREQUCE7HhBw7IQFSeipHCk84i0mWuljsbiPbAYA+5lOi+HJhZV89pFVxHZo8RddtkB1HmwhZZVgAqHTJ0kjdStYmzSs/B6NNLlH8CUeceB/S6s9eYZo4ANS6ZK95pLCKrtMt/yiMjJDtQ9dHiHAlo/CUArz8qyjCTkjIcTqGiIiIwoIjIURERIHgSEjIsRNCREQUCGbHhBynY4iIiCgsOBJCREQUAF4xNfQu206IQbpU3ahjGE24iSbcoGfUUdkOAEgXzxt1TGLajLKsgcPuGFF8mSaL1yPW6BRvE2mQHbeUUccnmEeYLmE118r3IcyOcXrlHztRkW7ZPtyybBprjLw4oDSjRtNRJM/o4mB1s8I1ISHHM5yIiIjCgp0QIiIiCovLdjqGiIhIwoAg14SErCUtBzshREREgWCKbshxOoaIiKgZy8vLQ9euXWGxWJCeno7PPvvsgvHr1q1Dz549YbFY0LdvX3z00Ud+v1dKYf78+ejQoQOio6ORmZmJQ4cO+cWcPn0aEydOhNVqRUJCAqZOnYqqqqp6j/Pb3/4WV111FaKiotCxY0c899xzomNjJ4SIiCgQKgQ3oTVr1iA3NxcLFizA7t270b9/f2RlZeH48eMNxm/fvh0TJkzA1KlTsWfPHmRnZyM7Oxv79u3zxSxatAhLly7F8uXLsXPnTsTGxiIrKwsOx49ZYhMnTsRXX32FTZs2YePGjdi2bRumT5/ut69Zs2bhzTffxG9/+1scOHAAf/7znzFo0CDR8RmUUs0qaaiiogI2mw03YgwiDAGm2emoDGWKjxfFe65OFe/jPzfFiuI7jSgWxU9N+VQUDwADhGm9NZp8xi7eKEunPOGNFsXrSdFNMNaI4v/tShTvo1KTHUeVV1bwTk/hvlNu2Tl43Cl7XwCAwyNrV5lT9jwBgFtY9E5a8M7lkZ/nJqMmiq91yF8/z3HZc9Vqn+x5avd5hSgeAAyHZZ9TWpW8GCQEf5Y8yo0t+BDl5eWwWq3yfQWg7u9Sl//3HIwW2fv2bJrDge/mPiFqa3p6Oq677jq8+uqrZx5D05CSkoKZM2fi8ccfrxc/fvx4VFdXY+PGjb77Bg8ejAEDBmD58uVQSiE5ORkPPfQQHn74YQBAeXk5EhMTsXLlStxzzz3Yv38/evfujc8//xwDBw4EAOTn5+O2227D999/j+TkZOzfvx/9+vXDvn370KNHD93PCUdCiIiImlBFRYXfzels+LpELpcLRUVFyMzM9N1nNBqRmZmJwsLCBrcpLCz0iweArKwsX/yRI0dgt9v9Ymw2G9LT030xhYWFSEhI8HVAACAzMxNGoxE7d+4EAGzYsAFXXHEFNm7ciNTUVHTt2hW//OUvcfr0adFzwU4IERFRAOqumBrMDQBSUlJgs9l8t4ULFza4v5MnT8Lr9SIx0X9kNjExEXa7vcFt7Hb7BePr/r1YTPv27f1+HxERgdatW/ti/vd//xffffcd1q1bh3feeQcrV65EUVERfvrTn17safR/XFE0ERHR5SpEV0wtKSnxm46JiooKqlnhoGkanE4n3nnnHVx11VUAgLfeegtpaWk4ePBgwFM0HAkhIiJqQlar1e92vk5I27ZtYTKZUFpa6nd/aWkpkpKSGtwmKSnpgvF1/14s5tyFrx6PB6dPn/bFdOjQAREREb4OCAD06tULAFBcHPi6IXZCiIiIAtHE2TFmsxlpaWkoKCjw3adpGgoKCpCRkdHgNhkZGX7xALBp0yZffGpqKpKSkvxiKioqsHPnTl9MRkYGysrKUFRU5Iv55JNPoGka0tPTAQDXX389PB4PvvnmG1/Mv//9bwBAly5dAj5GTscQEREFIBxVdHNzczF58mQMHDgQgwYNwpIlS1BdXY0pU6YAACZNmoSOHTv61pXMmjULw4cPx+LFizF69GisXr0au3btwuuvv36mDQYDZs+ejWeffRbdu3dHamoqnnzySSQnJyM7OxvAmRGNkSNHYtq0aVi+fDncbjdycnJwzz33IDk5GcCZharXXnstfvGLX2DJkiXQNA0zZszALbfc4jc6cjGXbSdEnJmsp4qucBfSKromHZOT0m2SIzzifZR4zOJtJKTptgBwSpOlqjaFEkdrUXxCpPy4Kz2ydEI9VZOlNB1Xjaxxys4pi1mWJq4nRTfCJKuabDbLB56l7z6Wim95xo8fjxMnTmD+/Pmw2+0YMGAA8vPzfQtLi4uLYTT+eG4NGTIEq1atwrx58zB37lx0794d69evR58+fXwxjz76KKqrqzF9+nSUlZVh6NChyM/Ph+Ws9OP33nsPOTk5GDFiBIxGI8aOHYulS5f6fm80GrFhwwbMnDkTw4YNQ2xsLEaNGoXFixeLju+yvU6IMS5OFO/t1028j/8Mk13PIjmzRBR/X8pWUTwgv05IvI7OV4nwuhFuYbl5PaSdkFMe2fkByK8Tsr86WRSvpxPyH0eCKN7plf8xrhF2Ok875Nd5qXTIFu5JOyGVtfJrP0RFyvbhdMuvE1Lzvew8bP2lrKPTdhevExKIur9LqU/9v6CvE3JkwdxGbeul5rIdCSEiIhIJUXYM/YidECIiogCEY01IS8fsGCIiIgoLjoQQEREFgtMxIcdOCBERUSCCnI5hJ6Q+dkICpSOJSJr0YTLIqnJajC7ZDgBECZNdTnl1pCYLWQyyRMTjXnnmisUgy2QwG2TplwDgFFa5jY1ouGjV+cSbHBcPOodbk52ELmE8IK+KG6Uj7bvaKHtuHcIqukZhRVxAnu3i8eiY/RZuoqS70NEkg45sRKLzYSeEiIgoEJyOCTl2QoiIiALBTkjIMTuGiIiIwoIjIURERAHgdUJCjyMhREREFBbshBAREVFYcDomUDqG0aQputJ0yjbGatkOAJzWZC95a6M8nbLcIzsOWfIsoOnoO2vi3EU5m0lWYK7cK0ttrdHk1YmlFWs9OlJ0E6JqRfGVblkxOgBQwuOQptwadFT2lSZx60l/hjBzWHwY8sxkeQXyloQLU0NO9Mm8bNky9OvXD1arFVarFRkZGfjrX//q+73D4cCMGTPQpk0bxMXFYezYsSgtLQ15o4mIiJpa3ZqQYG7kT9QJ6dSpE55//nkUFRVh165duPnmmzFmzBh89dVXAIAHH3wQGzZswLp167B161YcPXoUd911V6M0nIiIqMmpIG5Uj2hs/o477vD7+bnnnsOyZcuwY8cOdOrUCW+99RZWrVqFm2++GQCwYsUK9OrVCzt27MDgwYMbfEyn0wmn88crR1ZUVEiPgYiIiC5BuifKvV4vVq9ejerqamRkZKCoqAhutxuZmZm+mJ49e6Jz584oLCw87+MsXLgQNpvNd0tJSdHbJCIiosYTzCgIR0MaJO6EfPnll4iLi0NUVBTuu+8+fPDBB+jduzfsdjvMZjMSEhL84hMTE2G328/7eHPmzEF5ebnvVlJSIj4IIiKixsY1IaEnzo7p0aMH9u7di/Lycrz//vuYPHkytm7dqrsBUVFRiIqSr5YnIiKiS5u4E2I2m3HllVcCANLS0vD555/j5Zdfxvjx4+FyuVBWVuY3GlJaWoqkpKSQNThUpJUgDXrS0oSbRBplSX9lWoxsBwBsxh9E8eU60grjhdV9K4Wpp/FGWUooAHiFg34VmkW8D4cSVm4VniDSKssAEBshey2iTPKU7B+csvMwQsdxRJtlidxVDtkXmwiTvGqylK7is41dsNbEirgiTNENuaAvnqBpGpxOJ9LS0hAZGYmCggLf7w4ePIji4mJkZGQEuxsiIqKw4nRM6IlGQubMmYNRo0ahc+fOqKysxKpVq7BlyxZ8/PHHsNlsmDp1KnJzc9G6dWtYrVbMnDkTGRkZ582MISIiosuXqBNy/PhxTJo0CceOHYPNZkO/fv3w8ccf45ZbbgEAvPTSSzAajRg7diycTieysrLw2muvNUrDiYiImhSnY0JO1Al56623Lvh7i8WCvLw85OXlBdUoIiKiZoedkJBjATsiIiIKCxawIyIiCkCwi0u5MLU+dkIakbQAbZVLllb4H3cr2Q4AJEeUi+Kl6bMAEGuQpVNaDLInqkyTVZ/Vw2p0iLc5qmSvh1taZrkJSCs566GnmqxJ+OndFJ/1JmGlXqNR3ioV0bhHYvDoKKN7OeN0TMixE0JERBQIdkJCjmtCiIiIKCw4EkJERBQArgkJPXZCiIiIAsHpmJDjdAwRERGFBUdCiIiIAsDpmNC7bDshSk9VXCFN+Oy2jq4RxXeMlFXEBYAEYd6wpuRVNr3C0p/S+OSISlE8ABz1xIvipW0CgBhh9WCLUZbK7NBkVXoBwKOE1YPd8urBHk22D6OOT+LTtbJ2KeF5K40HAK/wuPV85Bg8wuMwCuN1lPY16ioH3EJwOibkOB1DREREYXHZjoQQERGJcCQk5NgJISIiCoDh/27BbE/+OB1DREREYcGRECIiokBwOibk2AkhIiIKAFN0Q4+dkEBpevLrpLuQbeBQ8pRNqUiDvMqmQ1h5t0yLEcWf8Db+zOoJj1W8jcUgS9Et98iqAdsiakXxABBtkqUBm6WlnwFIP0ZqhWmnABATJXtuHa7Gf29ERcqeK4dTR5s02XNl9Mo+pww68oab4vIGzRZHQkKOa0KIiIgoLDgSQkREFCiOZoQUOyFEREQB4JqQ0ON0DBEREYUFR0KIiIgCwYWpIcdOCBERUQA4HRN67IQESlidEgCk2a3SFF2LQZZ+CQDxBtkM3AllEu+jnUmWSnpKixXF6zluaTpzgklW0RgA7G6beBsJaUovANR6Gz9VVZrW6/bKUrL1iIzwiuKrHbK0cgBwCbOZIyNlbQIATyP/0dJTRZeXHqdQYieEiIgoEJyOCTkuTCUiIgpA3XRMMDc98vLy0LVrV1gsFqSnp+Ozzz67YPy6devQs2dPWCwW9O3bFx999JHf75VSmD9/Pjp06IDo6GhkZmbi0KFDfjGnT5/GxIkTYbVakZCQgKlTp6KqqqrB/R0+fBjx8fFISEgQHxs7IURERM3UmjVrkJubiwULFmD37t3o378/srKycPz48Qbjt2/fjgkTJmDq1KnYs2cPsrOzkZ2djX379vliFi1ahKVLl2L58uXYuXMnYmNjkZWVBYfD4YuZOHEivvrqK2zatAkbN27Etm3bMH369Hr7c7vdmDBhAm644QZdx8dOCBERUSBUCG5CL774IqZNm4YpU6agd+/eWL58OWJiYvD73/++wfiXX34ZI0eOxCOPPIJevXrhmWeewbXXXotXX331zCEohSVLlmDevHkYM2YM+vXrh3feeQdHjx7F+vXrAQD79+9Hfn4+3nzzTaSnp2Po0KF45ZVXsHr1ahw9etRvf/PmzUPPnj1x9913yw8O7IQQEREFJkSdkIqKCr+b0+lscHculwtFRUXIzMz03Wc0GpGZmYnCwsIGtyksLPSLB4CsrCxf/JEjR2C32/1ibDYb0tPTfTGFhYVISEjAwIEDfTGZmZkwGo3YuXOn775PPvkE69atQ15e3gWetAtjJ4SIiCgAoVoTkpKSApvN5rstXLiwwf2dPHkSXq8XiYmJfvcnJibCbrc3uI3dbr9gfN2/F4tp37693+8jIiLQunVrX8ypU6fw85//HCtXroTVKi/46Xtc3VteZvSksjX2Smgj5BVu7cIswQSjrHopAJQJq+jqSbmVqvTK01vF+9AsovhIg+zFKPPIU1vdmizF2mKSV9Etc8me2zhzw9/6LuREtSyN2+WRfbQZjfI3q9cr+w7ndsk/boVZ+xCeUrqq6FLwSkpK/P5wR0VFhbE1+kybNg333nsvhg0bFtTjcCSEiIgoECGajrFarX6383VC2rZtC5PJhNLSUr/7S0tLkZSU1OA2SUlJF4yv+/diMecufPV4PDh9+rQv5pNPPsFvf/tbREREICIiAlOnTkV5eTkiIiLOu16lIeyEEBERBcCgVNA3CbPZjLS0NBQUFPju0zQNBQUFyMjIaHCbjIwMv3gA2LRpky8+NTUVSUlJfjEVFRXYuXOnLyYjIwNlZWUoKiryxXzyySfQNA3p6ekAzqwb2bt3r+/29NNPIz4+Hnv37sWdd94Z8DFyOoaIiKiZys3NxeTJkzFw4EAMGjQIS5YsQXV1NaZMmQIAmDRpEjp27OhbVzJr1iwMHz4cixcvxujRo7F69Wrs2rULr7/+OgDAYDBg9uzZePbZZ9G9e3ekpqbiySefRHJyMrKzswEAvXr1wsiRIzFt2jQsX74cbrcbOTk5uOeee5CcnOyLOduuXbtgNBrRp08f0fGxE0JERBSIMFwxdfz48Thx4gTmz58Pu92OAQMGID8/37ewtLi4GEbjj5MaQ4YMwapVqzBv3jzMnTsX3bt3x/r16/06B48++iiqq6sxffp0lJWVYejQocjPz4fF8uO6tvfeew85OTkYMWIEjEYjxo4di6VLl+o/9vMwKNW8ViZVVFTAZrPhRoxBhCHAuhc6Fo0a4+JE8d4+V4j38f0I2WK6riO+FcXPTtkkigeAlIhyUbxZWgAH8oWpZVrjLxo94dG/ejtQJe7Wovgf3LLzQ8/C1Aq3bLGsHtKFqR5NPgvc2AtT9WiKhaneE7LXr/WXss/CdrtknwcAgG9KROHaea6yeUGCP0se5cYWfIjy8vKgsjQupO7v0jUTn4PJrP895XU5sOe9Jxq1rZcargkhIiKisOB0TCNSwi5ejVs2glCpawRB9s3HJT0IAF5hnU1phVuTjjHNeGFl36YYOWkbKfuGaNRReOK0SzZ6UuWWpwpaTLIU6x90jOgoaYXpSFmqscMt/yiMjpKlr2uafMRWWndXkx6GjoFwg57LFbQULGAXcuyEEBERBSCYInR125M/TscQERFRWHAkhIiIKBCcjgk5dkKIiIgCwOmY0GMnhIiIKBAcCQk5rgkhIiKisOBISCNSEbJur0E4VhdrlFcjldKE6bYA4Fayyq0JxhpRvDSlFwAcwnRmaYVbPRya7Diaok3WSId4G+k54tBxITGnMIU2wiS7yJ40BVjPNl6P7H0ByIfvjdJTREe6bTO7vmWT45RKaLETQkREFAildF1bxW978sPpGCIiIgoLjoQQEREFgNkxocdOCBERUSCYHRNynI4hIiKisOBICBERUQAM2plbMNuTP3ZCGpHRLUt/c3tlKXzHPfGieABoZ6oUxUsr4gLyVNITXtlxWAyyqq16tvEa5IOETmHKbblHljbcKkKWygwAEcJPPY+OqskVLoso3qvJ92Exy14/p1v2WkjT4wHA4RLuwyj/C6SM0nbJ3q9KR4ruZVxDl9MxjYDTMURERBQWHAkhIiIKALNjQo+dECIiokDwYmUhx04IERFRADgSEnpcE0JERERhwZEQIiKiQDA7JuSC6oQ8//zzmDNnDmbNmoUlS5YAABwOBx566CGsXr0aTqcTWVlZeO2115CYmBiK9oZPE+SlxZllVXHbRFSJ9xFj8IjiLToS20u9stTTWINLFB8pPAZAXnlXWgkYANpGyNKfjcJPJLvLKooHAKPw9atwxoj3oemoQCvlEaavS1NuI3Skz7rdsjYZxem2AITPbRO8FJc1TseEnu7pmM8//xy/+93v0K9fP7/7H3zwQWzYsAHr1q3D1q1bcfToUdx1111BN5SIiIhaFl2dkKqqKkycOBFvvPEGWrVq5bu/vLwcb731Fl588UXcfPPNSEtLw4oVK7B9+3bs2LEjZI0mIiJqcnXZMcHcyI+uTsiMGTMwevRoZGZm+t1fVFQEt9vtd3/Pnj3RuXNnFBYWNvhYTqcTFRUVfjciIqLmpm46Jpgb+ROvCVm9ejV2796Nzz//vN7v7HY7zGYzEhIS/O5PTEyE3W5v8PEWLlyIp556StoMIiIiusSJRkJKSkowa9YsvPfee7BYZPUizmfOnDkoLy/33UpKSkLyuERERCGlQnAjP6KRkKKiIhw/fhzXXnut7z6v14tt27bh1VdfxccffwyXy4WysjK/0ZDS0lIkJSU1+JhRUVGIiorS13oiIqImwuyY0BN1QkaMGIEvv/zS774pU6agZ8+eeOyxx5CSkoLIyEgUFBRg7NixAICDBw+iuLgYGRkZoWt1GBi0xj973Jos5S8Ssmq1AFCjZDNwlTpy/o4Lq+KWeWVpoZ0jT4viAcChmUXxMUZZujQA7Hcki+Klr7c1wiGKB4Bar+y4LSZ5heIaj2wfZpP8vHUaZOetUfhpL61gDQAmk7BCsUe+D+lfLS2C15+kS4vonR0fH48+ffr43RcbG4s2bdr47p86dSpyc3PRunVrWK1WzJw5ExkZGRg8eHDoWk1ERNTUNHXmFsz25CfkV0x96aWXYDQaMXbsWL+LlREREV3SeMXUkAu6E7Jlyxa/ny0WC/Ly8pCXlxfsQxMRETUbBgS5JiRkLWk5OIFIREREYcECdkRERIEI9qqnvGJqPeyEEBERBYApuqHHTkig9KxqFhbmFKcVQkfKn5CpCVZStRdWn/XqmFm1GGWVet3CVGYAaB8pKzlQ7Gwjind65bOnTVHhVpxariNFV5pC6xU+V0YdVXTNEbLjqK3ScT0k4etncsverwZ+M6cwYyeEiIgoEMyOCTl2QoiIiAJgUCqo0SOOPNXH7BgiIiIKC46EEBERBUKDeK1fve3JDzshREREAeB0TOhxOoaIiIjCgiMhgTI2fqqjR2v8PmGySZaqekKTnyJtTFWieIeKFO9DSlpF939d7RqpJT+yGGUVa6s0eYrnD65o8TZSsZGyc+pYpVW8jziLrKpxtVP2euupcOtVsrF1U6R8LF6TbiL8oq0M8s+1y/rS48yOCTl2QoiIiALBK6aGHDshREREAeAVU0OPa0KIiIgoLDgSQkREFAhOx4QcOyFEREQBMGhnbsFsT/44HUNERERhwZGQRiTt9drMDlF8mTdWtgMApd5yUXykjq67V9i3lVbqtRhkqa0AcApxonibqVa8jx88stejxitLI3UreRqp9Jz6wRkj3oeUV5Mnebo8sjRuaUVqg44Vg6oJKhQb3cJ9SMM5PSATpumYvLw8vPDCC7Db7ejfvz9eeeUVDBo06Lzx69atw5NPPolvv/0W3bt3x29+8xvcdtttZzVDYcGCBXjjjTdQVlaG66+/HsuWLUP37t19MadPn8bMmTOxYcMGGI1GjB07Fi+//DLi4s58lm7ZsgUvvfQSPvvsM1RUVKB79+545JFHMHHiRNGxcSSEiIgoECoEN6E1a9YgNzcXCxYswO7du9G/f39kZWXh+PHjDcZv374dEyZMwNSpU7Fnzx5kZ2cjOzsb+/bt88UsWrQIS5cuxfLly7Fz507ExsYiKysLDsePX1omTpyIr776Cps2bcLGjRuxbds2TJ8+3W8//fr1w//8z//gX//6F6ZMmYJJkyZh48aNouMzKNW8usIVFRWw2Wy4EWMQYQjw24+OC+4Y42TfjLWrU8X7+P4m2T66j/xGFD8uaZcoHgCuNh8VxesZCTmlyS6QpSlZX1jPSEixp7UoXs8ok3Qk5KRbdn44dVw4rtJjEcXrGQnRhF+/vy+3iffhbeQL+bl1XKzMZJK9Nxy1spEvAMD3sveS7ZDs4dt9LhsZBQB8UyIK16pkFy8EIBox8Cg3tuBDlJeXw2qVXwgvEL6/S9c9gYgI2XvqbB6PA1s+fw4lJSV+bY2KikJUVMMXI0xPT8d1112HV199FQCgaRpSUlIwc+ZMPP744/Xix48fj+rqar/OwODBgzFgwAAsX74cSikkJyfjoYcewsMPPwwAKC8vR2JiIlauXIl77rkH+/fvR+/evfH5559j4MCBAID8/Hzcdttt+P7775GcnNxgW0ePHo3ExET8/ve/D/g54UgIERFRAOpqxwRzA4CUlBTYbDbfbeHChQ3uz+VyoaioCJmZmb77jEYjMjMzUVhY2OA2hYWFfvEAkJWV5Ys/cuQI7Ha7X4zNZkN6erovprCwEAkJCb4OCABkZmbCaDRi586d531+ysvL0bq17Asf14QQEREFIkRrQhoaCWnIyZMn4fV6kZiY6Hd/YmIiDhw40OA2dru9wXi73e77fd19F4pp37693+8jIiLQunVrX8y51q5di88//xy/+93vGvz9+bATQkRE1ISsVmujTR2Fw+bNmzFlyhS88cYbuPrqq0XbcjqGiIgoEAqAFsRNOIjStm1bmEwmlJaW+t1fWlqKpKSkBrdJSkq6YHzdvxeLOXfhq8fjwenTp+vtd+vWrbjjjjvw0ksvYdKkSbIDBEdCmhWHV/Zy6Fmg2dYk2+Y7j3yhop5UUgk9VXfNBq8ovrGPAZCnkZa55RVxpYt+rWZ5avJpp2xBboxZft7WumXvDU24kDUyQnZ+AIDWBCm6njhpu2TnrTKxiq7E2es69G4vYTabkZaWhoKCAmRnZwM4szC1oKAAOTk5DW6TkZGBgoICzJ4923ffpk2bkJGRAQBITU1FUlISCgoKMGDAAABnFt7u3LkT999/v+8xysrKUFRUhLS0NADAJ598Ak3TkJ6e7nvcLVu24Pbbb8dvfvMbv8wZCXZCiIiIAqEQ5JoQ+Sa5ubmYPHkyBg4ciEGDBmHJkiWorq7GlClTAACTJk1Cx44dfYtbZ82aheHDh2Px4sUYPXo0Vq9ejV27duH1118HABgMBsyePRvPPvssunfvjtTUVDz55JNITk72dXR69eqFkSNHYtq0aVi+fDncbjdycnJwzz33+DJjNm/ejNtvvx2zZs3C2LFjfWtFzGazaHEqOyFERETN1Pjx43HixAnMnz8fdrsdAwYMQH5+vm9haXFxMYzGH0f+hgwZglWrVmHevHmYO3cuunfvjvXr16NPnz6+mEcffRTV1dWYPn06ysrKMHToUOTn58Ni+TH9+L333kNOTg5GjBjhu1jZ0qVLfb9/++23UVNTg4ULF/pl9wwfPhxbtmwJ+Ph4nZAANcV1QlJu+U4UP63Tp6J4AEi3yK4Tomc6Rs90iYT0CqsAUKbJjuM/7lbifVR5ZdcPOC28rsgxh3whm3Q6xqjjujDS6Zgyh3xaqbGnY/RMrUi30XOdEE+ZbJs2e2TTMW13V4jiAcDw72JRfEu6TsjN/R9DhKnhTJZAeLxOfPLFbxq1rZcajoQQEREFQkNwi2JYwK4eZscQERFRWHAkhIiIKABNnR1zOWAnJFCa/OTRhEsjnMIU3aSIMtkOIB8NTDLViPdh98rWX1QKa81YjPI0UofwxYgUpvQCgE34XB13xYvi25nlc+ulTtm8s3QNCQDUCivcRprkz221U7Y2wiNcExJhlI+TS+vNRETKj1vVNu5gtcHLP4oiYaqi25JxOoaIiIjCgiMhREREgeBISMixE0JERBQIdkJCjtMxREREFBYcCSEiIgoErxMScuyEEBERBYApuqHHTkiAlEk+cyXN8pRWVT3hlV/2t7XRIYqv1HEJdm8jz/JVa/LLJpcJ04YtRnml19Ne2WX6pa+3nsq+0srM0jYBQHykUxRf5pRftl1c5dYj3oWY0Sh7rpwO+cetFil9PWRf0/VU0TXqKJPRYnBNSMhxTQgRERGFBUdCiIiIAqEpQMdood/25IedECIiokBwOibkOB1DREREYcGRECIiooAEORICjoSci50QIiKiQHA6JuTYCQmQwSu/yox0/ZLT0/gvR7kwvbVayaqXAkCCsMrtCU2eaiwVY5SlkZp0LD6r0WTPlTQd1uOVp+jGRLhE8RUuefqsNA3YqeM4lJKlhTpcstTyCGkKMICoSFkat8cjn/3WvMJ0WOlpq+PiWYp/SCmE2AkhIiIKhKYQ1JQKs2PqYSeEiIgoEEo7cwtme/LD7BgiIiIKC46EEBERBYILU0OOnRAiIqJAcE1IyLETQkREFAiOhIQcOyEB0lMuQJhViASLLLW1wmuR7QBAj8jjovgyj6z6LABEGmSLryINspKnNTqq6HaM/EEU/x93K/E+ooSVdx2aLI00NkKWZgwAp12y108TVmEFAK+SLS2LNMoX5zndso8qc6TsnHJ75GnDtV5ZSrZHxz6ET60YS8tTuLETQkREFAiFIEdCQtaSFoOdECIiokBwOibkmKJLREREYcGRECIiokBoGnRd695vezobOyFERESB4HRMyHE6hoiIiMJCNBLy61//Gk899ZTffT169MCBAwcAAA6HAw899BBWr14Np9OJrKwsvPbaa0hMTAxdi8NEmm4LQNzFcwqrkSZFlMt2AHkKZkpEmXgfpzRZJVar0SGKNxvkFU8dSpYOe9obJ96HlFvJUjZPOuVtkr7eHk3+vUTPNo3NJUzptZhl6dUA4BV+KLicsnNQD2lKrzLIP9j0fBS2GBwJCTnxp8fVV1+NY8eO+W7/+Mc/fL978MEHsWHDBqxbtw5bt27F0aNHcdddd4W0wURERGGhqeBv5Ee8JiQiIgJJSUn17i8vL8dbb72FVatW4eabbwYArFixAr169cKOHTswePDg4FtLRERELYZ4JOTQoUNITk7GFVdcgYkTJ6K4uBgAUFRUBLfbjczMTF9sz5490blzZxQWFp738ZxOJyoqKvxuREREzY1SWtA38ifqhKSnp2PlypXIz8/HsmXLcOTIEdxwww2orKyE3W6H2WxGQkKC3zaJiYmw2+3nfcyFCxfCZrP5bikpKboOhIiIqFGpIKdiuCakHtF0zKhRo3z/369fP6Snp6NLly5Yu3YtoqNlixHrzJkzB7m5ub6fKyoq2BEhIqLmRwVZRZedkHqCWtaekJCAq666CocPH0ZSUhJcLhfKysr8YkpLSxtcQ1InKioKVqvV70ZEREQtX1AXK6uqqsI333yDn/3sZ0hLS0NkZCQKCgowduxYAMDBgwdRXFyMjIyMkDT2kiOc/jMKe9iVwlTYM9vI0mFtRnnlVim7xyaKb22qEu8jRljhVk/68wmPrAMdb5K9Fs5I+du13C2rtBwXKX+9jcIS08er5anG3kZOA5ZW6QUAj0fWJs2jIx1W+BmiI3udJDRN/qKcjWtC6hG98x5++GHccccd6NKlC44ePYoFCxbAZDJhwoQJsNlsmDp1KnJzc9G6dWtYrVbMnDkTGRkZzIwhIqJLH6djQk7UCfn+++8xYcIEnDp1Cu3atcPQoUOxY8cOtGvXDgDw0ksvwWg0YuzYsX4XKyMiIiI6l6gTsnr16gv+3mKxIC8vD3l5eUE1ioiIqLlRmgYVxHQMU3TrYwE7IiKiQHA6JuSaX9EHIiIiuixwJISIiCgQmgKE2WB+OBJSDzshAVIm+aCROL1OeHJHGjyyHQBoLUwLNekYeqzRokTx8aZaUbwbsuqzAFDmiRHFS9NtASBKmAbs9cSK4k3SnG8A8cKU22M18uN2a7LXI9IkzyOtdckq0JpMsufK5ZR/FJoihK+HWX7cHpPs/aciZGnABo1rFESUgvjaC/W2p7NxOoaIiIjCgiMhREREAVCaggpiOkZxJKQedkKIiIgCoTQENx3D6a9zsRNCREQUAI6EhB7XhBAREVFYNLuRkLqeogduwTVh5IWhjMolivd6ZFklAOB1yrIGPNWyTIaaSvlq+yqzbDhQV3aMW9Yur47XT6pGk7XJ4ZFnHmlG2TZOlyybxumRZwW53LLz3FMrL2DnEWbHeDzy89brkJ0jSsniNaf89YYwO0YTFrwDAE2WOAavS7YPj1dHgUrhZ6emZOf5mX0E/rnjgfv/Nmn8UQaPcgY1pVLXVvqRQTWz8aHvv/8eKSkp4W4GERFdQkpKStCpU6dGeWyHw4HU1FTY7fagHyspKQlHjhyBxSKrcN1SNbtOiKZpOHr0KOLj42Ew+H+bqaioQEpKCkpKSmC1yq9ncKnicfO4Lwc8bh63HkopVFZWIjk5GUZj460wcDgccLlko0ANMZvN7ICcpdlNxxiNxov2Zq1W62X1Zq3D47688LgvLzxu/Ww2W4hac34Wi4Wdh0bAhalEREQUFuyEEBERUVhcUp2QqKgoLFiwAFFRstoklzoeN4/7csDj5nHT5afZLUwlIiKiy8MlNRJCRERELQc7IURERBQW7IQQERFRWLATQkRERGHBTggRERGFxSXTCcnLy0PXrl1hsViQnp6Ozz77LNxNalS//vWvYTAY/G49e/YMd7NCbtu2bbjjjjuQnJwMg8GA9evX+/1eKYX58+ejQ4cOiI6ORmZmJg4dOhSexobQxY775z//eb3Xf+TIkeFpbAgtXLgQ1113HeLj49G+fXtkZ2fj4MGDfjEOhwMzZsxAmzZtEBcXh7Fjx6K0tDRMLQ6NQI77xhtvrPea33fffWFqcWgsW7YM/fr1810VNSMjA3/96199v2+JrzXJXBKdkDVr1iA3NxcLFizA7t270b9/f2RlZeH48ePhblqjuvrqq3Hs2DHf7R//+Ee4mxRy1dXV6N+/P/Ly8hr8/aJFi7B06VIsX74cO3fuRGxsLLKysuBwyKsaNycXO24AGDlypN/r/8c//rEJW9g4tm7dihkzZmDHjh3YtGkT3G43br31VlRXV/tiHnzwQWzYsAHr1q3D1q1bcfToUdx1111hbHXwAjluAJg2bZrfa75o0aIwtTg0OnXqhOeffx5FRUXYtWsXbr75ZowZMwZfffUVgJb5WpOQugQMGjRIzZgxw/ez1+tVycnJauHChWFsVeNasGCB6t+/f7ib0aQAqA8++MD3s6ZpKikpSb3wwgu++8rKylRUVJT64x//GIYWNo5zj1sppSZPnqzGjBkTlvY0pePHjysAauvWrUqpM69vZGSkWrdunS9m//79CoAqLCwMVzND7tzjVkqp4cOHq1mzZoWvUU2kVatW6s0337xsXmu6sGY/EuJyuVBUVITMzEzffUajEZmZmSgsLAxjyxrfoUOHkJycjCuuuAITJ05EcXFxuJvUpI4cOQK73e732ttsNqSnp7f41x4AtmzZgvbt26NHjx64//77cerUqXA3KeTKy8sBAK1btwYAFBUVwe12+73mPXv2ROfOnVvUa37ucdd577330LZtW/Tp0wdz5sxBTU1NOJrXKLxeL1avXo3q6mpkZGRcNq81XVizq6J7rpMnT8Lr9SIxMdHv/sTERBw4cCBMrWp86enpWLlyJXr06IFjx47hqaeewg033IB9+/YhPj4+3M1rEna7HQAafO3rftdSjRw5EnfddRdSU1PxzTffYO7cuRg1ahQKCwthMpnC3byQ0DQNs2fPxvXXX48+ffoAOPOam81mJCQk+MW2pNe8oeMGgHvvvRddunRBcnIy/vWvf+Gxxx7DwYMH8ac//SmMrQ3el19+iYyMDDgcDsTFxeGDDz5A7969sXfv3hb/WtPFNftOyOVq1KhRvv/v168f0tPT0aVLF6xduxZTp04NY8uoKdxzzz2+/+/bty/69euHbt26YcuWLRgxYkQYWxY6M2bMwL59+1rkWqcLOd9xT58+3ff/ffv2RYcOHTBixAh888036NatW1M3M2R69OiBvXv3ory8HO+//z4mT56MrVu3hrtZ1Ew0++mYtm3bwmQy1VsxXVpaiqSkpDC1quklJCTgqquuwuHDh8PdlCZT9/pe7q89AFxxxRVo27Zti3n9c3JysHHjRmzevBmdOnXy3Z+UlASXy4WysjK/+Jbymp/vuBuSnp4OAJf8a242m3HllVciLS0NCxcuRP/+/fHyyy+3+NeaAtPsOyFmsxlpaWkoKCjw3adpGgoKCpCRkRHGljWtqqoqfPPNN+jQoUO4m9JkUlNTkZSU5PfaV1RUYOfOnZfVaw8A33//PU6dOnXJv/5KKeTk5OCDDz7AJ598gtTUVL/fp6WlITIy0u81P3jwIIqLiy/p1/xix92QvXv3AsAl/5qfS9M0OJ3OFvtak1C4V8YGYvXq1SoqKkqtXLlSff3112r69OkqISFB2e32cDet0Tz00ENqy5Yt6siRI+qf//ynyszMVG3btlXHjx8Pd9NCqrKyUu3Zs0ft2bNHAVAvvvii2rNnj/ruu++UUko9//zzKiEhQX344YfqX//6lxozZoxKTU1VtbW1YW55cC503JWVlerhhx9WhYWF6siRI+rvf/+7uvbaa1X37t2Vw+EId9ODcv/99yubzaa2bNmijh075rvV1NT4Yu677z7VuXNn9cknn6hdu3apjIwMlZGREcZWB+9ix3348GH19NNPq127dqkjR46oDz/8UF1xxRVq2LBhYW55cB5//HG1detWdeTIEfWvf/1LPf7448pgMKi//e1vSqmW+VqTzCXRCVFKqVdeeUV17txZmc1mNWjQILVjx45wN6lRjR8/XnXo0EGZzWbVsWNHNX78eHX48OFwNyvkNm/erADUu02ePFkpdSZN98knn1SJiYkqKipKjRgxQh08eDC8jQ6BCx13TU2NuvXWW1W7du1UZGSk6tKli5o2bVqL6HQ3dMwA1IoVK3wxtbW16oEHHlCtWrVSMTEx6s4771THjh0LX6ND4GLHXVxcrIYNG6Zat26toqKi1JVXXqkeeeQRVV5eHt6GB+kXv/iF6tKlizKbzapdu3ZqxIgRvg6IUi3ztSYZg1JKNd24CxEREdEZzX5NCBEREbVM7IQQERFRWLATQkRERGHBTggRERGFBTshREREFBbshBAREVFYsBNCREREYcFOCBEREYUFOyFEREQUFuyEEBERUViwE0JERERh8f8BYsdHAxObD3kAAAAASUVORK5CYII=",
"text/plain": [
"<Figure size 640x480 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgcAAAGzCAYAAAC7ErTFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABf7ElEQVR4nO3deXgUVbo/8G93J+nOvq+QhCAQEEzAsAVFVgVUBGUEtwvuVwWvgo6KcxVlxosKig4i7qDjwqIDjqgoosCgAWVVRBAwkJCdJXu6O+k+vz+Y9I8iAertztIh38/z1POQzlt1TnUtnJxTbx2DUkqBiIiI6D+MrV0BIiIi8i5sHBAREZEGGwdERESkwcYBERERabBxQERERBpsHBAREZEGGwdERESkwcYBERERabBxQERERBrtunGwZMkSGAwGbN26tcm2mZubC4vFgu+//77JtumOTp064eqrr26y7R06dAgGgwHz5s1rsm3SSW3hWN1www2YOHFik23PHQaDAdOmTWvSbV555ZW46667mnSb3oDnFHmqXTcOmsPs2bMxYMAAXHLJJa1dFaIm8+ijj+KTTz7Brl27ROstWLAAoaGhqK2tbaaaue/777/H119/jUcffbS1q9IuuXtOUctg46AJlZSU4N1338U999zT2lUhalJ9+vRB37598cILL4jW+/zzz3HFFVfA19e3mWrmvrlz52LEiBHo0qVLa1elXXL3nKKWwcZBE3r//ffh4+ODsWPHtnZViJrcxIkT8c9//hOVlZW64qurq7FhwwZcddVVzVwzueLiYnz++efs1m5l0nOKWs552Tg4fPgw7rvvPqSmpsLf3x+RkZG4/vrrcejQoUbjq6ur8d///d+IjIxESEgIJk+ejBMnTmhitm7dilGjRiEqKgr+/v5ISUnB7bffrolZtWoVBgwYgKCgINdn06ZNQ1BQEKqrqxuUe+ONNyIuLg4OhwNXX301Onfu3Gj9MjMz0bdvX81n77//Pvr374+AgACEh4fjsssuw9dff91g3U2bNqF///6wWCzo3Lkz3nvvvQYxpaWlePDBB5GYmAiz2YwuXbrgueeeg9PpbLQ+8+fPR3JyMvz9/TFkyBDs3r3b9bvFixfDYDBgx44dDdb7v//7P5hMJuTl5TW6XYmtW7fCYDDg3XffbfC7r776CgaDAatXr8bHH38Mg8GADRs2NIh7/fXXYTAYXPW/9dZbERQUhLy8PIwfPx5BQUGIjo7Gww8/DIfDoVm3qqoKDz30kOs7S01Nxbx589DYJKfny7G6/PLLUVVVhbVr1zZa1unWrVsHm82GMWPGnDVu3rx5GDRoECIjI+Hv74+MjAx8/PHHZ4z/4IMPkJqaCovFgoyMDGzcuFHz+4qKCjz44IPo1KkTzGYzYmJicPnll2P79u2umM8//xx1dXUYOXKk6zOeU95/TlELUuehFStWqPT0dPXkk0+qN954Qz3++OMqPDxcJScnq6qqKlfc4sWLFQB10UUXqcGDB6u///3vaurUqcpoNKrLLrtMOZ1OpZRSRUVFKjw8XHXr1k3NnTtXvfnmm+ovf/mL6tGjh2tbdrtd+fv7qxkzZmjqsnHjRgVALV++XPN5VVWVCgwMVFOnTlVKKfXee+8pAOrHH3/UxB06dEgBUHPnznV99tRTTykAatCgQWru3Lnq5ZdfVjfddJN69NFHXTHJyckqNTVVxcbGqscff1y98sor6uKLL1YGg0Ht3r1bU4+0tDQVGRmpHn/8cfXaa6+pyZMnK4PBoB544AFXXHZ2tuu76tSpk3ruuefU008/rSIiIlR0dLQqLCxUSilVXl6u/P391UMPPdTguFx44YVq+PDhZz94Ap07d1ZXXnllg89vu+02FR4erux2u6qurlZBQUHqvvvuaxA3bNgw1bNnT9fPU6ZMURaLRfXs2VPdfvvtatGiRWrChAkKgHr11VddcU6nUw0fPlwZDAZ15513qldeeUWNHTtWAVAPPvigpozz6VjV1taeMb4x99xzj+rbt+854zp27Kjuu+8+9corr6gXX3xR9e/fXwFQq1ev1sQBUL169VJRUVFq9uzZ6rnnnlPJycnK399f/fLLL664m266Sfn5+akZM2aot956Sz333HNq7Nix6v3333fF3HnnnSoyMrJBXXhOefc5RS3nvGwcVFdXN/gsKytLAVDvvfee67P6xkFGRoay2+2uz59//nkFQH366adKKaVWrlypAKiffvrpjGUeOHBAAVALFizQfO50OlWHDh3UhAkTNJ8vX75cAVAbN25USilVVlamzGZzg4vk+eefVwaDQR0+fFgppdT+/fuV0WhU1157rXI4HA3KqpecnKzZvlJKFRcXNyjjr3/9qwoMDFS///67ZluPPfaYMplMKicnRyn1/28O/v7+6siRI664LVu2KABq+vTprs9uvPFGlZCQoKnf9u3bFQC1ePHiM3yDcjNnzlS+vr7q+PHjrs9sNpsKCwtTt99+u6Y+MTExqq6uzvVZQUGBMhqNavbs2a7PpkyZogBoPlNKqT59+qiMjAzXz6tWrVIA1N/+9jdN3J/+9CdlMBjUgQMHlFLn57Hq1q2bGjNmTIPPG5OUlKRmzZp1zrjTr1e73a569erV4D8SAAqA2rp1q+uzw4cPK4vFoq699lrXZ6Ghoa5G95lceumlmmNaj+fUSd56TlHLOS8bB6ey2+3q6NGjqqSkRIWFhWla4fWNg9dff12zTkVFhfLx8VH//d//rZRS6rvvvlMA1KxZszSNiFPVXySn/nVS78EHH1T+/v6qoqLC9dmECRNUhw4dNBf0+PHjVWJiouazjIwMlZmZ6fp57ty5CoDasWPHWfc7OTlZXXjhhQ0+T0tL09xI09LS1OjRo1VJSYlm+eabbzT7U39zuPHGGxtsc8CAASo1NdX185dffqkAqG+++cb12UMPPaT8/f1VeXn5WestsXPnTgVAvfXWW67PPvvsMwVAffXVV67P6m+8p9ZnwYIFCoDat2+f67P6G3lxcbGmnP/5n/9R4eHhrp/vvvtuZTKZGuxLfQO0voF4Ph6rAQMGqH79+p11f5RS6pdffmm0J+xcjh8/rkpKStS9996rwsLCNL8DoLkW6k2aNEkFBAS4/qNOTk5Wffv2VXl5eWcsp0ePHmrkyJENPuc55b3nFLWs8/KZg5qaGjz55JOuMbSoqChER0ejtLQUZWVlDeK7du2q+TkoKAjx8fGuZxSGDBmCCRMm4Omnn0ZUVBTGjRuHxYsXw2azNdiWamR8cNKkSaipqcG//vUvAEBlZSW++OILXH/99TAYDJq43NxcZGVlAQAOHjyIbdu2YdKkSa6YgwcPwmg04sILLzzn95CUlNTgs/DwcM3zFPv378eaNWsQHR2tWerHYouLizXrn/5dAUC3bt00z3NcfvnliI+PxwcffAAAcDqd+OijjzBu3DgEBwefsb52ux2FhYWa5fRx2VOlp6eje/fuWLZsmeuzZcuWISoqCsOHD3d9Nnr0aISGhjaI6927N7p166bZpsViQXR0tOaz07+zw4cPIyEhocG+9OjRw/V74Pw8VkopzTl7Jp9//jliY2MbPCvTmNWrV2PgwIGwWCyIiIhAdHQ0Fi1apOtard+n6upqlJSUAACef/557N69G4mJiejfvz+eeuop/PHHH43uy+l4TnnvOeWujRs3YuzYsUhISIDBYMCqVauarSwAeOqpp2AwGDRL9+7dm7XM5nBeNg7uv/9+PPPMM5g4cSKWL1+Or7/+GmvXrkVkZOQZH7I5G4PBgI8//hhZWVmYNm0a8vLycPvttyMjI8P1lG1kZCQANHiQEQAGDhyITp06Yfny5QCAzz77DDU1NZr/9AFg7NixCAgIcMUtX74cRqMR119/vbjOAGAymRr9/NSbotPpxOWXX461a9c2ukyYMMGtcm+66SZ88sknsFqt+O6775Cfn49bbrnlrOv98MMPiI+P1yy5ublnXWfSpEn47rvvcPToUdhsNvzrX//ChAkT4OPj44oxm80YP348Vq5cibq6OuTl5eH7779v8P3X1701tJVjdeLECURFRZ1zu1988QVGjx59zpv+v//9b1xzzTWwWCx49dVX8cUXX2Dt2rW46aabGv3PW4+JEyfijz/+wIIFC5CQkIC5c+eiZ8+e+PLLL10xkZGRjV6rAM8pbz2n3FVVVYX09HQsXLiw2co4Xc+ePVFQUOBaNm3a1GJlNxWfc4e0PR9//DGmTJmiyZ+1Wq0oLS1tNH7//v0YNmyY6+fKykoUFBTgyiuv1MQNHDgQAwcOxDPPPIMPP/wQN998M5YuXYo777wTSUlJ8Pf3R3Z2dqNlTJw4ES+//DLKy8uxbNkydOrUCQMHDtTEBAYG4uqrr8aKFSvw4osvYtmyZRg8eDASEhJcMRdccAGcTif27NmD3r17C7+Zhi644AJUVlZqnto+m/379zf47Pfff0enTp00n02ePBkvvPACPvvsM3z55ZeIjo7GqFGjzrrt9PT0Bk8tx8XFnXWdSZMm4emnn8Ynn3yC2NhYlJeX44Ybbmg07t1338W6devw22+/QSnV6I1cj+TkZHzzzTeoqKjQ/CW0d+9e1++B8+9Y1dXVITc3F9dcc81Zyy0tLcUPP/yg622Gn3zyCSwWC7766iuYzWbX54sXLxbtU0BAgOav8/j4eNx333247777UFxcjIsvvhjPPPOMK3Oie/fu+OSTTxotg+fUmbXWOeWJMWPGnDVjxmaz4S9/+Qs++ugjlJaWolevXnjuuecwdOhQt8v08fE5573L252XPQcmk6nBXx0LFiw4Yxf1G2+8oXmD26JFi1BXV+c6oU6cONFge/UXZv3Qgq+vL/r27XvGVzFPmjQJNpsN7777LtasWXPG/OpJkyYhPz8fb731Fnbt2tXgZjN+/HgYjUbMnj27QS+IO39pTZw4EVlZWfjqq68a/K60tBR1dXWaz1atWqVJRfrxxx+xZcuWBhdfWloa0tLS8NZbb+GTTz7BDTfcoPnLqzHh4eEYOXKkZrFYLGddp0ePHrjooouwbNkyLFu2DPHx8bjssssaxI0cORIRERGuuP79+yMlJeWs2z6TK6+8Eg6HA6+88orm8/nz58NgMLi+i/PtWO3ZswdWqxWDBg06az3rU+quuOKKc+6TyWSCwWDQXJuHDh06Y9dvVlaWJiUxNzcXn376Ka644gqYTCY4HI4GwxExMTFISEjQDANmZmbixIkTjQ438Jw6yZvOqeY0bdo0ZGVlYenSpfj5559x/fXXY/To0Y02hPTav38/EhIS0LlzZ9x8883Iyclpwhq3kFZ4zqHZTZ48WZlMJvXAAw+o119/Xd16662qY8eOKjIyUk2ZMsUVd3oq44IFC9S0adOU0WhUl156qevBwPnz56uuXbuqRx55RL3++utq3rx5KjU1VYWEhKg//vjDtb158+Yps9msysrKGq1Xly5dVHBwsAKgtm3b1mhMTU2NCg4OVsHBwcpkMqmioqIGMU888YQrlWnevHlqwYIFavLkyeqxxx5zxSQnJ6urrrqqwbpDhgxRQ4YMcf1cVVWlLr74YuXj46PuvPNOtWjRIjVv3jw1ZcoUFRgYqEpKSpRSjacyzZ49W0VERKjIyEiVn5/foKx58+a5njDfsmVLo/vbFP72t78po9GoAgIC1P3333/GuDvvvFMFBQUpg8GgXnjhhQa/r9/n082aNUudeqk4HA41bNgwZTAY1N13360WLlyoxo0b12ja2fl0rObNm6cCAgLO+VDp5MmT1dChQ88aU2/dunUKgBo8eLBatGiRevrpp1VMTIxKS0tTp9+ecIZURovFonbt2qWUUurEiRMqMDBQTZkyRb344ovqjTfeUBMnTlQANMe8sLBQ+fj4NHgYuR7PKe86p5oKALVy5UrXz4cPH1Ymk6nBw6sjRoxQM2fOdKuML774Qi1fvlzt2rVLrVmzRmVmZqqkpKQW28emcl42Dk6cOKFuu+02FRUVpYKCgtSoUaPU3r17VXJycqONgw0bNqi7775bhYeHq6CgIHXzzTerY8eOueK2b9+ubrzxRpWUlKTMZrOKiYlRV199tSalSqmT70Pw8fFR//jHPxqt11/+8hcFQHXp0uWs9b/55psVgEafpq73zjvvqD59+iiz2azCw8PVkCFD1Nq1a12/13tzUOpkdsbMmTNVly5dlJ+fn4qKinLdeOqzM+pvDnPnzlUvvPCCSkxMVGazWQ0ePNh1Yz5dQUGBMplMqlu3bmfdX0/t37/fdRPatGnTGePWrl2rACiDwaByc3Mb/F7vjVypk9/Z9OnTVUJCgvL19VVdu3ZVc+fO1WSa1DtfjtWAAQPULbfccsbfK3UynS4mJkY9//zzZ4071dtvv626du2qzGaz6t69u1q8eHGj3zkANXXqVPX++++74vv06aO+++47V4zNZlN//vOfVXp6ugoODlaBgYEqPT1d806Betdcc40aMWJEo3XiOeU951RTOr1xsHr1agVABQYGahYfHx81ceJEpZRSv/32m+tcONNy6jsmTnfixAkVEhKiyYBpC87LxkFruv3229Wll17a2tXwCiUlJcrHx6dBjjd5n3Mdqx07diiDwXDOFLr6lN5ff/21GWrZtDZu3KiMRmODHH9qGk11TjWl0xsHS5cuVSaTSe3du1ft379fsxQUFCilTjY4f/vtt7Mup6eqnq5v376anp224Lx8ILE1zZo1C926dcP333/f7mdmXLJkCRwOB/7rv/6rtatC53CuY/Xss8/iT3/6k66H4P7v//5PV6pdaxs8eDCuuOIKPP/883jzzTdbuzrnnaY8p5pLnz594HA4UFxcjMGDBzca4+fn51EqYmVlJQ4ePNjm7oMGpdzMFyI6g2+//RZ79uzBE088gWHDhuGf//xna1eJzoDHipqat51TlZWVOHDgAICTjYEXX3wRw4YNQ0REBJKSknDLLbfg+++/xwsvvIA+ffqgpKQE69atQ1pamluThj388MMYO3YskpOTkZ+fj1mzZmHnzp3Ys2dPg/ddeLXW7rqg88+QIUOUr6+vGjp0qOZVq+R9eKyoqXnbOVX/htvTl/rnz+x2u3ryySdVp06dlK+vr4qPj1fXXnut+vnnn90qb9KkSSo+Pl75+fmpDh06qEmTJrlef92WsOeAiIiINM7L9xwQERGR+9g4ICIiIo1my1ZYuHAh5s6di8LCQqSnp2PBggXo37//OddzOp3Iz89HcHBws07GQUREbZ9SChUVFUhISIDR2Hx/71qtVtjtdo+34+fnd843v3qF5niQYenSpcrPz0+988476tdff1V33XWXCgsLa/Rtf6fLzc095wsnuHDhwoULl1OXxl5C1VRqampUXIypSeoZFxenampqmq2uTaVZHkgcMGAA+vXr53pPuNPpRGJiIu6//3489thjZ123rKwMYWFh+PvGdPgH6ZvNbHtlsriOX+zvKYp3njCfO+g0ykf21Rrssp4SU428lZzyaYUo/vDoEHEZxtpzx2jizzwrc6N8K2XxAABhJ1TYAflfCHUW2fEoT5J13FV2cmNGUeHVHbX93DENC5GFB+U1nOr8XGqi/ETx/kWyMnJH+YviAcB8TLbjIYeEJzoAa4TsnKroJNt+YN65Y05nHSK7h9zfY724jL3V8bpj7VW1eHvMapSWliI0NFRclh7l5eUIDQ1F9rZkhAS73ztRXuFESsZhlJWVISREfm9tSU0+rGC327Ft2zbMnDnT9ZnRaMTIkSORlZXVIN5ms2kmRKmoOHni+QeZEKCzcWCGr7iexgBht461BRoHJtnNxqjkJ6mPSfY/t8mN7i9pz560cWByp2dP+B+Yj48bNwBf2Toms+zyM/q70TgQrmKS/R/8n0Jk4T4+8uFCH19ZxaRlGN04z01mWRk+vvLGgclPdk4ZhbvhzvE2BcjuIf5B8v9mzEb5Pb0lhqFDgo0eNQ7akibfy6NHj8LhcCA2NlbzeWxsLAoLCxvEz5kzB6Ghoa4lMTGxqatERETkMYdyery0Fa3eBJo5cybKyspcS25ubmtXiYiIqAEnlMdLW9HkwwpRUVEwmUwoKirSfF5UVIS4uLgG8WazGWazvMueiIioJTnhhCd/+3u2dstq8saBn58fMjIysG7dOowfPx7AyQcS161bh2nTpunezs/VibrHnWxO+fiU9AHD8F/knSzBR+pE8TUR+p6xqFeVIB9jyx4fLIoPypW3dP2PyS4Ah59sP2oD5ftdZ5Gt47+/RFxG/pgOshWEp1SPlxsOy51L9k0JoniDU37zCv+p6NxBpzg+IPbcQacxyi4lHLxeNvgemCs/pxzC8foT3WTXNwB0/KZMFB++V3YvlH5PAJAWLbs2Fs27VlxGaQ/99x2n1QpgpbgMOrtmec/BjBkzMGXKFPTt2xf9+/fHSy+9hKqqKtx2223NURwREVGzcygFhwcJfp6s29KapXEwadIklJSU4Mknn0RhYSF69+6NNWvWNHhIkYiIqK3w9LmBdv3MQb1p06aJhhGIiIjIOzRb44CIiOh84oSCgz0HREREVK89DSu0+nsOiIiIyLt4bc/B1zndYQrQl25YUyN/B6jyl73KtCJZ3o4yl8lSl6yRsnSqhH9XieIB4NhFAaL4wCL5K199y2XrnEiVHb/yC+Stb6Pwlcu/3yNLAQQAn2pZfGCebD9qukTJCgAQvk92LKTv8geAE/1kDxqHZNeIy6jqKJv7IO4H2fZ9bMJcSQDGWtnxO9FNnnK9/2ZZ6nHKv2QneuAR+fHeI3yw3DK2VFxGwL/Ddcc6bC33Ny6zFYiIiEjD+Z/Fk/XbCg4rEBERkQZ7DoiIiHRweJit4Mm6LY2NAyIiIh0c6uTiyfptBRsHREREOvCZAyIiImq3vLbn4IKIo/AN1JfidvC4PMUrNSlXFL/dnCQuo7Raln5lOSbrc8obGiiKB4DaYFkZgfKJAHG8u2zGS0uprD0duUs+g96JHrJ45UazOXqXLB2uPFF2+dmD5bP6GYR/qrgz02dFiiy+uJ/sugCAuM2y+MoOsgMYfES2fQCojJd9VyarvE85KFe2H0eGyWZZtMXKUzjD/h0kiq+6VJ5yXdtFf72cNfJ9cJcTBjggv0ZOXb+t8NrGARERkTdxqpOLJ+u3FRxWICIiIg32HBAREeng8HBYwZN1WxobB0RERDq0p8YBhxWIiIhIgz0HREREOjiVAU7lQbaCB+u2NK9tHPw9eTWCg/V1bHwVK08znL18oije5EYfS1Cu7NFUo/D1WfGf5YjiAcCeEiOKP3SVLDUKACzHZPHOStkFc+JC2fYBIHq77LsNKLCJy6hIkqVw2kNl2z8WIU9lTNgk2w+HRX6iW47Kjp8yupEu2VEWH7mnVrh9+a1QOjtoUI78u/UvkeWiKoOsDHuKPA3QcrXsAi89HCkuIyiuUneso1p+rbqLwwpERETUbnltzwEREZE3ccAIhwd/UzuasC7NjY0DIiIiHZSHzxwoPnNARER0fuEzB0RERNRuseeAiIhIB4cywuHOrGyu9ZuwMs3MaxsH17z8AExmfWl03SbuE2/fp1rWveNXIS4CtcJJE8t6ytKWjvZOlBUAIC5LdnaG7hcXAb9K2X44zLJj0eE7WZoaAPhWylK2jveQzxxYEyM8p8pl2w8slM8Gf6KrvplN6wUdlnd71obI4qN3yo9f7gjZrUqZfEXx1ij5XduvTPafRPzGUnEZ5d1kX25FV9l5bjTI97vmy1hRfPhIYW4zgPLfw3XHOq1W8fbd5YQBTg863J1oO60DDisQERGRhtf2HBAREXmT9vRAIhsHREREOnj+zAGHFYiIiKiNYs8BERGRDicfSPRg4iUOKxAREZ1fnB6+PrktZSt4bePAWKdgNOr7IvMqhdPbATBnytJrysoCxGWYjshmNAz9VTbrnrlMntpmC5Gd2P4nmr+M2gBZa/pYuvy0jdwlS20r7ywuAibh5HDSVFd3Zoo8mi47B6N3yN/+np8qO34VJfLjl7BJVi/pzI+2EvkN31gnuzZK+snvU9K37SaukcUHHLHLVgBw4CbZ7KOB38lnZdSZxQ4AcNg4Ot4cvLZxQERE5E3a0wOJbBwQERHp4ISx3bwEiY0DIiIiHRzKAIcHMyt6sm5L42ANERERabDngIiISAeHh9kKDg4rEBERnV+cyginBw8kOvlAoucevm85AoL1pfY9s3eMePsdQ8tE8ZVVspQwAIjaKTsRDE5ZapQyycev7NLdKBUXIZ6V0bdKtn1zmXy/zWXSVDhZWikA2CJk++1fJLvJVCXIUsgAAMJM1PzL5EX4Co9HyGF5+pyloFIUX91JNpth9FeHRfEA8Nv/JovifU/Iz9vI3bJ7yJHhsnMqME+eXulTJatT+Jh8cRknvkzQHWuQZ/iSDl7bOCAiIvImHFYgIiIiDSc8yziQv1Ku9TBbgYiIiDTYc0BERKSD5y9Bajt/j7NxQEREpIPnr09uO42DtlNTIiIiahFe23OwpbIzzNA3k17poTDx9ktDZVPiRWyUp5Edu0gWXxssi+/2D1l6FwBUXSpL8TrWS57SF3hE9kRudbzsAR+DfOJAlApnDuy0ukZcxolUf1G8PUy2/eJ+sngASFxbK4qvDZIf74JLZce7qJ/8Wuqwvk4U7/SVHe+iq1NE8QAQvF9Whl2eNYjiDFm8T7UsvrKL7PwAAD/hrJq3J20Sl7H5pi66Y+2Vdvy2UFyEW5wwwAlPHkjk65OJiIjOK/XDCp4sEnPmzEG/fv0QHByMmJgYjB8/Hvv27TvneitWrED37t1hsVhw0UUX4YsvvhDvKxsHREREOtS/58CTRWLDhg2YOnUqNm/ejLVr16K2thZXXHEFqqrO/Oa4H374ATfeeCPuuOMO7NixA+PHj8f48eOxe/duUdleO6xARETUnq1Zs0bz85IlSxATE4Nt27bhsssaf53pyy+/jNGjR+PPf/4zAOCvf/0r1q5di1deeQWvvfaa7rLZc0BERKSDUxk8XgCgvLxcs9hs+t4BXVZ28rX/ERERZ4zJysrCyJEjNZ+NGjUKWVlZon1l44CIiEgHp4dDCvXvOUhMTERoaKhrmTNnzrnLdjrx4IMP4pJLLkGvXr3OGFdYWIjY2FjNZ7GxsSgsLBTtK4cViIiIWlBubi5CQv5/5pjZfO4MnqlTp2L37t3YtEme/eEOr20cZL3WFyY/fVMIdrxJ1iICgJplcbIV3Jgvwx4rS78K+ENf6ma9w1fL0hLd0fE7eUpfSbospS/502PiMqSKB525G64xlUnyWThjv5XNPnfgTv0zzwFA7I+icABAzhjZJd5hvfzt7wH5svRHy1H5xeRTKsvRs3UJkG2/Rl6nOossLc0aJ7sfAEBAjuz4VSfLypCmJQKAb49yUfxLv48QlxHgpz/Fsq6q5aZl9HzK5pPrhoSEaBoH5zJt2jSsXr0aGzduRMeOHc8aGxcXh6KiIs1nRUVFiIuT/Z/HYQUiIiIdHDB4vEgopTBt2jSsXLkS3377LVJSzv0+jszMTKxbt07z2dq1a5GZmSkq22t7DoiIiNqzqVOn4sMPP8Snn36K4OBg13MDoaGh8Pc/2UM7efJkdOjQwfXcwgMPPIAhQ4bghRdewFVXXYWlS5di69ateOONN0Rls+eAiIhIh/phBU8WiUWLFqGsrAxDhw5FfHy8a1m2bJkrJicnBwUFBa6fBw0ahA8//BBvvPEG0tPT8fHHH2PVqlVnfYixMew5ICIi0sEBiIcGTl9fQqlzPwuzfv36Bp9df/31uP7664WlabHngIiIiDTYc0BERKRDU2UrtAXixsHGjRsxd+5cbNu2DQUFBVi5ciXGjx/v+r1SCrNmzcKbb76J0tJSXHLJJVi0aBG6du0qKsfoUDA69KUXHa+SpS0BgCFS1jVU179CXIb/Ltk0i9YYWRqZsVbevdXlvROi+CNXylIAAcAeIksLK+kvKyP0kDx1ySk80ysT5Bexsa8sVcgpnACxqJ8b3ZlO2bEI/EOWpgYARpvsPLeHyGd+3D8lShSfsEnWgRvwzc+ieAAwpCSK4jv888zvwz+T6l6ydNdDibLz1lQtP6fsdtnxC/5YON0sgBefWaA7tqrCiZHnDmsS7kyedPr6bYW4plVVVUhPT8fChY3Pkfn888/j73//O1577TVs2bIFgYGBGDVqFKxWq8eVJSIiai3qP1M2u7uoNjRls7jnYMyYMRgzZkyjv1NK4aWXXsL//u//Yty4cQCA9957D7GxsVi1ahVuuOEGz2pLREREza5J+ziys7NRWFiomfQhNDQUAwYMOOOkDzabrcEkFERERN6mfljBk6WtaNKa1r+gQTLpw5w5czQTUCQmysbxiIiIWkJTzcrYFrR6M2bmzJkoKytzLbm5ua1dJSIionatSVMZ6yd2KCoqQnx8vOvzoqIi9O7du9F1zGazrhmpiIiIWlP91MuerN9WNGnjICUlBXFxcVi3bp2rMVBeXo4tW7bg3nvvFW2rzt8A5aevC6ZTxHFpVRE0vuDcQac48E6quAxzhSw10SGblBGVifIuqupOspkcQw7JZ+mr6Ci7AAzCdLuaKOEXBSDmJ1kq6vGe8vSr4r6y/e60WpbBUxssv1zzhsjWOTJKnroqnbG0Tp55jND9skJsYbJj4XNJT1E8AFj2yu4hpZckictw+Mnifcpl94SO38rTKw/Ey2Zdtdwm+54AoLAuVHdsdZ30vYPu83RooC0NK4jvNpWVlThw4IDr5+zsbOzcuRMRERFISkrCgw8+iL/97W/o2rUrUlJS8MQTTyAhIUHzLgQiIiLyXuLGwdatWzFs2DDXzzNmzAAATJkyBUuWLMEjjzyCqqoq3H333SgtLcWll16KNWvWwGKxNF2tiYiIWpgTRjg9GBrwZN2WJm4cDB069KyTQRgMBsyePRuzZ8/2qGJERETexKEMcHgwNODJui2t7TRjiIiIqEVw4iUiIiId+EAiERERaSgPZ2VUbegNiV7bODie7oDRX1+KStnmTuLt+1hlLbhAozBfC0DhIFl80GHZiZP0aYmsAAA510SL4u1h8v2O+FWW/ljWRXjBSKczBGCNkKUmRm+Tz8JZmSRLE80fLEsJ8y+RH4vIn2Xr+NjkqauFA2XXUtf35a9I33+L7Ls1H5OdU04fYc4ggKp4WWqiQf7VojxFth9h+2THO//PdaJ4AAj3lc3saq2T/zfz2fE+umPtlXYAv4rLcIcDBjg8mDzJk3VbWttpxhAREVGL8NqeAyIiIm/iVJ49NyB831urYuOAiIhIB6eHzxx4sm5Lazs1JSIiohbBngMiIiIdnDDA6cFDhZ6s29LYOCAiItKhPb0h0WsbB4GHfWAy66veiEk/ire/Me8CUbz/rjBxGbYjspS76ljZ0yq5V8nSEgEganetKP5EV/kMiJbjsvSomjLZ6FZwrnwWtsoE2bEwZueLy4BTlm4XdER2vE12+dNMtf6ym1FtkPzmFfq7LL4mIVBcRsK/ZXmATh/ZdxVQaBPFA0BJH9n0kmEHZNceAFgjZNdGTYzs+Nl+l52zAJD0VqEo/reZkeIygp7Rn3pcVyeb3ZT08drGARERkTdpTw8ksnFARESkgxMevj65DT1z0HaaMURERNQi2HNARESkg/IwW0G1oZ4DNg6IiIh04KyMREREpMEHEr1AVac6GP31pcRlvdxPvP1u/71PFF+RVSouoyxFli5pFGbo1cTLp3mz5clS+qJ3ytOEckeaRfEBwqzB4C93y1YA4DPkQlF81aAu4jLqAmXpcyV9ZcfPfFw+G2XEHlkZlWHym5c9VBbvVyHfj6JLZN9txC7ZfhT3lc2QCQDRO2TfbXEfeVqwPVy235ajsr9Mw4RpqABQkRYjijdZ7OIyrnn9O92x1so6bOwvLoLOwWsbB0RERN6EwwpERESk0Z5en9x2BkCIiIioRbDngIiISAcOKxAREZFGe2occFiBiIiINLy258AUVAtjgL6UJ/PNJeLt/7g/RRQfMlGWngcA4ftlsxOWOmWHw+TGZGS1QbL4/Ess4jIc/sIZ9Hxlremc/0kXxQOAQZj16V8inwExME+2Hz7Vsra5qVaeulraRTjjZY68DIdZVobDT/7Xk3+e8Luyyo5f1C758S5LkdUp5LD8uzUelMUXCFM+rb3lN5G6atl9yijfbTgED+5JYj3VnnoOvLZxQERE5E3aU+OAwwpERESkwZ4DIiIiHRQ8e1eBfPCq9bBxQEREpEN7GlZg44CIiEiH9tQ44DMHREREpOG1PQeOWiOUXV/b5cgB2SxhAKDMsikQbWHiIlAdJZt9LuyALPWxLEV++Mq7yvKKur57Ql5G9zBRfGBulSg+e3yAKB4A/EplLfaQQ/KZ5I4M8xPF+1TL6pT8pnwKvdoJ3UTxJpt8VDR56RFR/N7/SRCXEXJA9l1VdhTOTnhQvt8OYXazQ5iyCwDWMNk6PebLjsWB58JE8QAQEiW7XoMsNnEZy3IydMc6qmwA1ovLcEd76jnw2sYBERGRN2lPjQMOKxAREZEGew6IiIh0UMoA5cFf/56s29LYOCAiItLBCYNH7znwZN2WxmEFIiIi0mDPARERkQ7t6YFE720cWE2AQV8qYFSn4+LN+5lkqYxFx2PFZZjLZWmDx3vIDof5uDz9KkA4u93hsRHiMpLWlInia8NkMz/6VMkvsA4bZOlXx3rJ0yV9K2X1sofJjl+RMC0RAMJ/l6WRFWTKZ+Es7ttBFB+XJZ+mr2Co7Hq1FMjSiKuj5Z2ooQdl+xH5Q4G4jN8ejhPFx2/wF8WHrQ4UxQNASV/ZeZvef4+4jH4hh3XHWivrsE1cgnva0zMHHFYgIiIiDe/tOSAiIvIiHFYgIiIijfY0rMDGARERkQ7Kw56DttQ44DMHREREpMGeAyIiIh0UACVPEtOs31Z4bePg0rR98AvSN8vdpg29xNu3HJV170TnytOv8sbI0q+6/EM2E2DeYFnaEgAYZRM/wimbaBAAcOCmYFF8tzeOiuI7vyefMfH3e2Xpdh3WC78oALYIX1G8jzD10SnbPADgeKps6sCEjdXiMhwBstvI8VT5SeVTJrutRuyVXXv5Q0XhJ8vYJet4Le0rS0sEgI5rZfstnRG1Jlreza0Ca0Xx/97VXVxG/IBy3bE2u6w+nnDCAAPfkEhERETtkdf2HBAREXkTZisQERGRhlMZYGgn7zngsAIRERFpsOeAiIhIB6U8zFZoQ+kKbBwQERHpwGcOvMCvx+JgqtGXhtXxO3na2dE0WV5YRUc3RmCcshPhaE9ZamJAkbwZahRm/fgfk6WEAUDhANmMeIXDo0XxtUHyC6zzoz+I4h1DLxaXYT4uu5xs4bL9cMiyEgEAvsLjXdFJPitjWWfZtSE9BwEgdqsslbgyXnYOmmXZtACAmihZfMSvVnEZB26SHY/OH8u+3MoOsu8JADp3KhbFl1vl51Ssr/5URquv/P5P5+a1jQMiIiJvwp4DIiIi0mC2AhEREWnUP5DoySK1ceNGjB07FgkJCTAYDFi1atVZ49evXw+DwdBgKSwsFJXLxgEREZGXqqqqQnp6OhYuXChab9++fSgoKHAtMTExovU5rEBERKTDyb/+PXnmQL7OmDFjMGbMGPF6MTExCAsLkxf4H+w5ICIi0qH+gURPFgAoLy/XLDabrcnr2rt3b8THx+Pyyy/H999/L15f1HMwZ84c/POf/8TevXvh7++PQYMG4bnnnkNqaqorxmq14qGHHsLSpUths9kwatQovPrqq4iNjRVV7MlunyMgWF+azV+SbxdtGwCMg06I4isr5ek4HVbKZp/zrZTNNpgzSt7xE75H1uqt6ChPdTLWycow2WTN6YpkefM7e06mKD5uizyFsyJFFu8wy9LzYrfItg8AlcIUXN8q+XerhKdIrUVehj1Qth+xP1aI4vMHy2YSBQCTcHLQ0m4B4jK6P7FXFF8+vJsoPmHdMVE8ABwO6iiKDxgkzxN9983RumMdNiuAb8VltKbExETNz7NmzcJTTz3VJNuOj4/Ha6+9hr59+8Jms+Gtt97C0KFDsWXLFlx8sf4UbdH/Lhs2bMDUqVPRr18/1NXV4fHHH8cVV1yBPXv2IDAwEAAwffp0fP7551ixYgVCQ0Mxbdo0XHfddW61XIiIiLyF+s/iyfoAkJubi5CQENfnZrMbLzE5g9TUVM0f7IMGDcLBgwcxf/58/OMf/9C9HVHjYM2aNZqflyxZgpiYGGzbtg2XXXYZysrK8Pbbb+PDDz/E8OHDAQCLFy9Gjx49sHnzZgwcOFBSHBERkddoqvcchISEaBoHza1///7YtGmTaB2PnjkoKysDAERERAAAtm3bhtraWowcOdIV0717dyQlJSErK6vRbdhstgbjL0RERNQ0du7cifj4eNE6bmcrOJ1OPPjgg7jkkkvQq1cvAEBhYSH8/PwaPCEZGxt7xhzLOXPm4Omnn3a3GkRERC2jqcYVBCorK3HgwAHXz9nZ2di5cyciIiKQlJSEmTNnIi8vD++99x4A4KWXXkJKSgp69uwJq9WKt956C99++y2+/vprUbluNw6mTp2K3bt3i7sqTjdz5kzMmDHD9XN5eXmDhzWIiIhanYfDCnBj3a1bt2LYsGGun+v/v5wyZQqWLFmCgoIC5OTkuH5vt9vx0EMPIS8vDwEBAUhLS8M333yj2YYebjUOpk2bhtWrV2Pjxo3o2PH/P7kaFxcHu92O0tJSTe9BUVER4uLiGt2W2Wxu0ocxiIiImkNrTNk8dOhQqLOsuGTJEs3PjzzyCB555BF5QacRNQ6UUrj//vuxcuVKrF+/Hikp2tytjIwM+Pr6Yt26dZgwYQKAk29pysnJQWamLJVs+rc3weivL32wxw2HRdsGgHKbLDXRsTlcXIblmGwWtpJ0abqkLBUOAPyPy9Zxp5FscMgeZQk5JMsJqwuQNyZ9K2RXpX+hfAa9+B9kx8/vhGy/bVGy1FgAsEXIDmD8i/J8ydq7ZQ8a+5XL75Che2WpiQdukaUmKoP8WvKpln23XRZmi8vIntZDFO8vnKn1+PhIUTwAdLniD1F8zV8a/8PwbA5fqX8/nFZP+vnpTESNg6lTp+LDDz/Ep59+iuDgYNdzBKGhofD390doaCjuuOMOzJgxAxEREQgJCcH999+PzMxMZioQEVGbxlkZz2DRokUATnZznGrx4sW49dZbAQDz58+H0WjEhAkTNC9BIiIiatOUwb3u1FPXbyPEwwrnYrFYsHDhQvEkEUREROQdOPESERGRDq3xQGJrYeOAiIhIj1Z4z0Fr4ayMREREpOG1PQevjVyCwGB9bZdbN8lnZezWsUgUf7yv/LXO1dlBovjQw3Wi+Mhf5TMHmn/4TRR/6M+9xWUkfl0lii9Nlc1WF5tVJooHgLxhoaL48H3iIpA3RNbWjtkqS3001sn/7AjfI1tHZabJy9gnm262OlaeklkwRHj8dsv22xopf1DMXCorwxEnTxtM/pds9lhbbKAo/ng/+X4XV8nua1UD5bNR9h+6R3dsbZUd8iRR9zBbgYiIiBpqQ0MDnuCwAhEREWmw54CIiEgHDisQERGRVjvKVmDjgIiISBfDfxZP1m8b+MwBERERaXhtz8Hzh0fDJ1Df7HsBv0pnMwQKtiWL4msuks2gBwD2vrJWotFhEsUH5soPX/lVslS1hH/L0yVLMmTpVHGbSkXx+28JEcUD8jRAh7/8uw05KDvex3rKth+9S94nWRskq5M9XJ5mWGeR/Y1R2lX+N0lgnmzffWyy+Kok+ayMoYdk8QdvlJ+3BuHlVxtTK4pP/Jf8L9kj18rSSsMuOyYu46LgPN2xVoNsnz3CYQUiIiLSaEeNAw4rEBERkQZ7DoiIiPTglM1ERER0qvY0KyOHFYiIiEiDPQdERER6tKMHEr22cXBRWD7MQb66Ym++bYt4+09/O14UbwqQzZgIAOE/6Kt/vap42XhUTaz8TIveKouv6CBLrwSAWtmkbajqJFsh+JC8w6syUfZdlaTJ02NDcmV5Z3X+su+2MFM+XukjmyATsR/Jp6PMvaeXKD78d3naYMElsvjgP2TnSNIaeTpcWYrs+g793Y3r9UfZrIz7HpLNgFicIdsHAAiLkM2KGvyGLPURAPKeCtMda7e1ZCpj+3nmgMMKREREpOG1PQdERETexKBOLp6s31awcUBERKQHnzkgIiIiDT5zQERERO0Vew6IiIj04LBC61ub0x2mAH2zMu6aK5tpEABC+8p23XZZjbiMYxmylC1ziSy1LW6zPL2yJF2WuhRQJD+bK5Nl6zh9ZF1t0duF+XkA4JTNFOnUd+ppWEpkM3earLJz0OCQp51F7baJ4vNul6UlAoBTmO1a1E/etRqxWxYfWCi7Nmoi5bdCy3HZ9V3SV1wEnL7honhlle23kmcqo+qXCFF80P2F4jJ+KknSHeuokp3jHmlHjQMOKxAREZGG1/YcEBEReZV21HPAxgEREZEezFYgIiKi9oo9B0RERDrwDYleINhihclf3zfp8A8Tb78ySXaUTL8LZxMCEFAh60LquK5CFF8b7CeKBwB7uGy/AwvERcB8TNYhdUz4gHxRX9nkMgAQmC+Lt0bKr+Ls8bIUB/8CYcedGz2SBQNldVLyhAgoo+y7shyT74jluGxSq/Jk2a0t/rtjongA2D9Z9tR+6O/iIhB2QJYBUxskO972MPl5bo+XTXTkdKMrfVLSNt2x1so66I/2UDt65oDDCkRERKTBxgERERFpeO2wAhERkTcxwMNnDpqsJs2PjQMiIiI9mMpIRERE7RV7DoiIiPRoR9kKXts4OP5zNIwWi67Y5L/uE28/b12qKN63Ut4dFL1TloZ0aKwsXTJ5daUoHgCit8ty1Yr7iYuAySqLt5TIvltDkPxYmGzC1FWrvIzOK0pF8XUh+s7vetnj5LNBdf5ENmHY8V7yNNGaaNl3ZXJjnhx7sKyTUzoZVN7ISNkKALq+UyyKz70mRlxG/mDZMbcmye45Ydvl6dDOFNnkTrUfxMrLeFT/8XaqFuwAb0eNAw4rEBERkYbX9hwQERF5E74hkYiIiLQ4rEBERETtFXsOiIiI9GhHPQdsHBAREenAZw68gNMXgM6sux93dBVv3xggO0oqtUpcRqGfLDUxeqds5rmjveUzRTossrSz0N/lZ3N5Z1l8xF7ZLG8FmfKpA6Wz20XtktUJABwBsrSwvMv8RfHm46JwAMAfE2RlBOTLUzh9qmXxlZ2c4jKCjsjOQ1u4LJfRFiYKBwDkXCtLTQzJlu+39H27BTGykWJbuGz7AOCwy77b4sGy1EcASLXon0a1ulZ23yR9vLZxQERE5FXa0euT2TggIiLSg88cEBER0ana0zMHTGUkIiIiDfYcEBER6cFhBSIiItLwcFiBjYMmsOlPCxGicya2iz77H/H2HcGytCLzTnnaoDVWVkZlgixFyKdKfqYFFcjSisqT5adI8hrZtIxHhspmJ4zeJU9dOjJclmYYtUueLnn0ItmTyPGbZemSxy6U1yl8j6xONdHiIlCdKDseQdnCKRMBHE2TnevxWbLv1j9fnqpccFmYKL7wUvn12uOZQ6L4Oossj/jocOEUqgBUteyeEHhIfg+Zfvw23bFOqxXAX8Rl0Nl5beOAiIjIq3BYgYiIiDTaUeOA2QpERESkwZ4DIiIiHfieAyIiImq32DggIiIiDa8dVvjN7otAu762iyWqRrx9w8/BoviwA/L0ucoaWcpWwhf6ZyIDgCNjE0TxAGBwyuoU/+1RcRnH+0SI4hO+l82YaI2Qn7Zhe2Xx1cLZ7QAg+JCsz7BgoGw/Er+Vp50VDJSliYYdlM8caLLJzqmgPPm1FJwrS8ksGCT7bqN3yu4HABCUL9sPh1mewpk38QJRvG+l7Bz0O2QWxQOAT5XsWFT3lJ+3T/X/l+7Ymso63DdLXIR7+EBi4xYtWoS0tDSEhIQgJCQEmZmZ+PLLL12/t1qtmDp1KiIjIxEUFIQJEyagqKioyStNRETU0uqfOfBkaStEjYOOHTvi2WefxbZt27B161YMHz4c48aNw6+//goAmD59Oj777DOsWLECGzZsQH5+Pq677rpmqTgREVGLUx4sbYio723s2LGan5955hksWrQImzdvRseOHfH222/jww8/xPDhwwEAixcvRo8ePbB582YMHDiw0W3abDbYbDbXz+Xl5dJ9ICIioibk9gOJDocDS5cuRVVVFTIzM7Ft2zbU1tZi5MiRrpju3bsjKSkJWVlZZ9zOnDlzEBoa6loSExPdrRIREVHz8aTXoI31HogbB7/88guCgoJgNptxzz33YOXKlbjwwgtRWFgIPz8/hIWFaeJjY2NRWFh4xu3NnDkTZWVlriU3N1e8E0RERM2NzxycRWpqKnbu3IktW7bg3nvvxZQpU7Bnzx63K2A2m10PONYvREREBGzcuBFjx45FQkICDAYDVq1adc511q9fj4svvhhmsxldunTBkiVLxOWKc8L8/PzQpUsXAEBGRgZ++uknvPzyy5g0aRLsdjtKS0s1vQdFRUWIi4sTV+y1oqHwrdQ3k57aI09DirjkzL0ZjSmIlk9XF/OjLC0s72pZamLYH7IZFgHA6StLQ8q9OkpchrTr7ER32WyDBqdsHwAgZofsu6rsKE+XTNxSKYoP/11WxsE/ydPOANk56FcuH2l0CieLNNXK/3yqs8iOuVM2CSfsgfL9Ppoh+24Dj4iLQGWyrIyEf8vij/YThQMAfISzMr6YuVxcxtaqFN2xNqtsBk6PtEIqY1VVFdLT03H77bfresA/OzsbV111Fe655x588MEHWLduHe68807Ex8dj1KhRusv1+D0HTqcTNpsNGRkZ8PX1xbp16zBhwgQAwL59+5CTk4PMzExPiyEiImpVTfX65NMfvDebzTCbG2/8jxkzBmPGjNFdxmuvvYaUlBS88MILAIAePXpg06ZNmD9/fvM1DmbOnIkxY8YgKSkJFRUV+PDDD7F+/Xp89dVXCA0NxR133IEZM2YgIiICISEhuP/++5GZmXnGTAUiIqL25vQH72fNmoWnnnqqSbadlZWlSQwAgFGjRuHBBx8UbUfUOCguLsbkyZNRUFCA0NBQpKWl4auvvsLll18OAJg/fz6MRiMmTJgAm82GUaNG4dVXXxVViIiIyCs10bBCbm6u5vm6M/UauKOwsBCxsbGaz2JjY1FeXo6amhr4+/vr2o6ocfD222+f9fcWiwULFy7EwoULJZslIiLyfk3UOGgLD99z4iUiIqLzRFxcXINpC4qKihASEqK71wDw4omXiIiIvElTPZDYnDIzM/HFF19oPlu7dq04McBrGwfDI36Df5DO6o08d8jpNm9NFcUbIm3nDjpNsPB9TpUd9bfqAKAkTX74fGXZdjDKJkwEANTJdgPhe4UzyVXJZ/Ur7C+bES/xG/nxLsyUpdQGFsj244J/yut0ZKjsYPjUyO9eRmFGbd4QeYel8hHWSzy5pDw9tsN6WXxpZ3ERCM6WfVeFA2X7EbFTFA4AqBpRIYqf8f0kcRnBO/SPwTtsVgCrxWW4pRVSGSsrK3HgwAHXz9nZ2di5cyciIiKQlJSEmTNnIi8vD++99x4A4J577sErr7yCRx55BLfffju+/fZbLF++HJ9//rmoXK9tHBAREXmVVmgcbN26FcOGDXP9PGPGDADAlClTsGTJEhQUFCAnJ8f1+5SUFHz++eeYPn06Xn75ZXTs2BFvvfWWKI0RYOOAiIjIaw0dOhRKnblV0djbD4cOHYodO3Z4VC4bB0RERDq0hWcOmgobB0RERHq0wrBCa2EqIxEREWmw54CIiEgHDit4gQ8fvRI+PhZdsaX3CfPzAIQclHWaGH4X5ucB+ONaWT5V8EHZ9q1uTJhY3UF2dvpUylO86gJkZYTvLhPFHx4bLooHAD9ZEfCtkOdw1sTIpgI01MnSK/0qvPNyPdFTdp7HbJGfU4EFsnxJh1l2fVfFy7/bgkHC/TCI8ysB4Qyk0pRPW5i889heECiKdycNvCZO/344rS34Py6HFYiIiKi98s4/RYiIiLxNO+o5YOOAiIhIBwPceZemdv22gsMKREREpMGeAyIiIj04rEBERESnYiqjF7jludW6Z2X8V0lv8fazL5elFZXtjRSX4VMlG2GK2yTLtyu9UDYLIACUdZGNJMVuqxWXYbLKvtvfbwsVxV+wrEoUDwBHRsjSr0pTZfEA4FcqTTuTbb9gkPxyDSiQ3Y0i9si/W4ND9l3Z5actSlNlaaLRO2XnoEO2eQDyGROls5UCgF+F7PiVyiabRWi2PL1SGWX7bbfpS0k/1bhRm3XH2iprsfApcRHuaUc9B3zmgIiIiDS8tueAiIjI67Shv/49wcYBERGRDu3pmQMOKxAREZEGew6IiIj0aEcPJLJxQEREpEN7Glbw2sbBJzcPgY/JrCu2aLA8zbB0sFUUbzC6cVSFM6o5QmT5VNK0RAAwCicbdCd9Lvw32XcVkCfbj/zL5GmG9lBZnapj5d+tb6Vwxssa2fbtofKXr5Z3kaWqxa0XVgpA1A+yWVH33RsjLiNyl2zfj6fKZrw0OkThAAD/IukMiPLjV2ORreMjzEQ12eQ7XhMvndlVfi2t/K237lhntRXAKnEZdHZe2zggIiLyKhxWICIiolO1p2EFZisQERGRBnsOiIiI9OCwAhEREWmwcUBERESnak/PHHht46BwSCRMfvpm87KHyLfvrJU9bmEwy49q1HZZfN5g2bRt9lD5jGpBObL9jt5VJy7Dr0w2k2OdWbbf5nL5fkfvkq1TlOErLsPp27yzMprkWYaI+kMWf3h8tLiMmjjZd9v1o2pxGTlXBIninb6y69UgTDsGgJpesninn/y8jdU/OSEAoCJZdn0f7yH/L8B8TBYfvVN+D9l462LdseUVToSLS6Bz8drGARERkVfhsAIRERGdyqAUDMr9/+E9WbelMZWRiIiINNhzQEREpAeHFYiIiOhU7SlbgcMKREREpOG1PQd33LUa/kH6qvd+zgDx9qv/kKVsBeTKZnkDgKO9Zc3EzqtkKV4+R2Wz4QHA8b5Ronhpeh4AlHfSl4JaryZGVkbAUXnzu7iPLDUxNFueduYQfldl3WTbN9rkxyLix2JRvKVUPsOp+TubKL44Q5aWCAA+wuxHqzAjM0qY6gq4cby7yo9fnSzLF0pYhOW4/FqSzi6ZO0pcBK7ZP1p3bG2VHcAb8kLcwWEFIiIiOhWHFYiIiKjdYs8BERGRHhxWICIiolO1p2EFNg6IiIj0aEc9B3zmgIiIiDS8tucg2xoNs4++9LOS7+PF2zd2lU1xV9VF3uQL2yVLn8seJ8tbMtqEeU4Aon6R7UdBpjyFM/w3WbxRNokjfGrkaWfSMkJXCKfUBFByW4YoXgmb5rVuzMLpCA8UxTt95Ol2tUGy87xOViUAgF+Z7LyN2i2bCbBwoHwWTv8iWZ0i9sjvIf7FshP3aH/Z9WqLcufvQ9l+jO2/Q1zChiNddMc6qmWptJ5qS0MDnvDaxgEREZFXUerk4sn6bQSHFYiIiEiDPQdEREQ6MFuBiIiItJitQERERO0Vew6IiIh0MDhPLp6s31Z4beNg5a4+MPrrm92v25Ac8fadEM4s9n1HcRlBeQ5RfHWsrCPHclyednY0TRbf6UuruIw6iyydqjbITxR/eIz8tDU4ZFdl5bg+4jJMdlm8b7ns+MVnyVO28ocEi+KrOsrvXqn/K8tdNV50kbiM6jjZdyVNr+y0qlQUDwA18bKczNKu8nTJY2myayN8l2z7pT3kx/uCT2Qn+vr8fuIy6gQTuzps8nuU2zisQERERO2V1/YcEBEReRNmKxAREZFWO3oJEhsHREREOrSnngM+c0BEREQa7DkgIiLSox1lK3jUOHj22Wcxc+ZMPPDAA3jppZcAAFarFQ899BCWLl0Km82GUaNG4dVXX0VsbKxo2yZLHYz++mZWO7Y0UVp1VHSSxYcdEBeBwgGyjhmfGlm6VkCRPA3J6SOrU+5w+cyPJmlmkTAjMzhb3uFlD5XFn+gqLgIdNlSJ4u0hbkxPKORTLYvv8bw8Lfjgw7LUxJjtshkTASAgR/bdOi2yW9v+m4UnCIDIn2Xx7lyvRrvs4jgxQJZmGLFZlioJAGWPVIrinTZhji+A7/u9qTu2osKJlLniItzCYQUdfvrpJ7z++utIS9Mmzk+fPh2fffYZVqxYgQ0bNiA/Px/XXXedxxUlIiKiluFW46CyshI333wz3nzzTYSHh7s+Lysrw9tvv40XX3wRw4cPR0ZGBhYvXowffvgBmzdvbrJKExERtbj6bAVPljbCrcbB1KlTcdVVV2HkyJGaz7dt24ba2lrN5927d0dSUhKysrIa3ZbNZkN5eblmISIi8jb1wwqeLG2F+JmDpUuXYvv27fjpp58a/K6wsBB+fn4ICwvTfB4bG4vCwsJGtzdnzhw8/fTT0moQERFRMxH1HOTm5uKBBx7ABx98AItF8PLrs5g5cybKyspcS25ubpNsl4iIqEmpJljaCFHPwbZt21BcXIyLL77Y9ZnD4cDGjRvxyiuv4KuvvoLdbkdpaamm96CoqAhxcXGNbtNsNsNsNrtXeyIiohbSnrIVRI2DESNG4JdfftF8dtttt6F79+549NFHkZiYCF9fX6xbtw4TJkwAAOzbtw85OTnIzMwUVWzkBb/DT+fMal8fvvjcQaepDZfNmFg8UFwEwnbLRm3iNh4Xxe+/NUwUDwDKKEunMtbJZ36M2S77bgsyZbM4Rv4q2z4AlCfJjoU9TFwEigbIUhODjsj2ozxZ3oiuiZHFV/STzz5qj5TtR0lveQZ1aLBsdsnyTrLHqTpskKdXSs+pwGJ5KqMtXLYfkT/IUhMN8ksJ4ZYaUXzd6ihxGRkF03XHOmusAJ4Ql0FnJzq7g4OD0atXL81ngYGBiIyMdH1+xx13YMaMGYiIiEBISAjuv/9+ZGZmYuBAN/53JSIi8hZOdXLxZP02osnfkDh//nwYjUZMmDBB8xIkIiKiNo1vSNRv/fr1mp8tFgsWLlyIhQsXerppIiIir2GAh88cNFlNmh8nXiIiIiINTrxERESkh6dvOWxDb0hk44CIiEgHpjJ6gRsjsxAYrG/U45taeSpj5y6Nv7HxTHKPhonLiPtBlh5VFyZ7sZSlRD4qZI2RpVN1ezVPXMbxzARRvE+1bCQu6ECZKB4ASi+IFMWH7Xcj7SxEth91/rLj53DjvWMG4W64k2YYv1FWiMEpz58rzpB9V11fOSTb/qhkUTwARP0iS+lz+MtSdgHAZJetc6yv7J4T+aP8eB8sjBbF/+8Dn4jLmBh0RHdseYUTHR4SF9GmLFy4EHPnzkVhYSHS09OxYMEC9O/fv9HYJUuW4LbbbtN8ZjabYbXKpsvlMwdERER6tMIbEpctW4YZM2Zg1qxZ2L59O9LT0zFq1CgUFxefcZ2QkBAUFBS4lsOHD4vLZeOAiIhIB4NSHi9SL774Iu666y7cdtttuPDCC/Haa68hICAA77zzzpnraTAgLi7OtcTGxorLZeOAiIioBZ0+E7HNZms0zm63Y9u2bZqZjo1GI0aOHHnGmY4BoLKyEsnJyUhMTMS4cePw66+/iuvIxgEREZEeziZYACQmJiI0NNS1zJkzp9Hijh49CofD0eAv/7PNdJyamop33nkHn376Kd5//304nU4MGjQIR47of44D8OIHEomIiLyJu0MDp64PnJzhOCQkxPV5U04+mJmZqZnLaNCgQejRowdef/11/PWvf9W9HTYOiIiIWlBISIimcXAmUVFRMJlMKCoq0nx+tpmOT+fr64s+ffrgwIEDojp6bePAF07onV/M4MbMgUXlslne/HYGics4MvLcMaeqlVUJkbvlKWFRP8vSzkqGdBCX4RSeVR2/k6WEFQyVpSUCgPmErLVfniwfcbMck5URsfOEKD5/eIQoHgDCfpcd79KubqTHhgtTUfPkaaJKmAX426OJovgeL8pTdvc8qu/mXM9ULf9ugw/K4gMO65vJtl7kLxWyAgD4/qlaFB9hqhSX0W/LbecO+g9HtQ3As+Iy3NLCcyv4+fkhIyMD69atw/jx4wEATqcT69atw7Rp03Rtw+Fw4JdffsGVV14pKttrGwdERERepRXekDhjxgxMmTIFffv2Rf/+/fHSSy+hqqrK9S6DyZMno0OHDq7nFmbPno2BAweiS5cuKC0txdy5c3H48GHceeedonLZOCAiItKhNd6QOGnSJJSUlODJJ59EYWEhevfujTVr1rgeUszJyYHR+P97pU6cOIG77roLhYWFCA8PR0ZGBn744QdceOGFonLZOCAiIvJi06ZNO+MwwukzI8+fPx/z58/3uEw2DoiIiPTgxEtERER0KoNTPl/J6eu3FXwJEhEREWl4bc/BFxVpMCt9aTmqlzwdBz+GisLtYfLuIL9SWYqXjyxDCOYyeSpjzhWyQ+5bIU8TDcqVfVf2MFn6lVU2KRwAIHm17Bw53kuYVwrAoTf39j8qusnOQctx+TkYdKTx17KeicEpfxnLsTThbJQB8ttOwkbZbIPm47L9rukiP6l8ymX5ld1elb2hDgD2zIwXxRutsr/3ahL8RfEAUG2THYtn/jpFXMbV08/8auDT2SprsV9cgps4rEBEREQaLfyeg9bEYQUiIiLSYM8BERGRDk01t0JbwMYBERGRHu3omQMOKxAREZEGew6IiIj0UAA8eVdB2+k48N7GQaxfGfz99FXPVmoRb9+WKkzxMsqPavJq4Yx43WRpRSe6CnPnAHR977go/tjF4eIyTshe4Y0TPWQpYZG/yI9FcT9ZamJZN3kZCf+WrWOslcUfTROFAwCUQXZtlF8gLyPpa6soPudy+fV61CK7VcXskH23/jllongAMNpls4MqX/nt1lQhuzb8ymVppbnXyNISASDaV7bOhfftFpdR49Cf3myXZ3S7jc8cEBERkZaCh88cNFlNmh2fOSAiIiIN9hwQERHp0Y6yFdg4ICIi0sMJQP5Gee36bQSHFYiIiEiDPQdEREQ6MFvBC7y27zKYAvTNEOcfUSPeftQ/AkTxRy+Sf1UHbpKl/IT9Kp3dThQOADhwiyw1sdO/5N9tRZKsYnE/1oriq2PkxyLumyJRvMMcJy7Dr1SWHlvZUZaK2mlluSgeAIoyZTM/hu6X37zyB8lScJO/FE4/CqCyoyz90WST5bcdvFGWlghA3L18LFN+TvlUe9KHfW7RG2QzogJAeWfZDJbbyuUzXsZemas7tq5Kdt15pB09c8BhBSIiItLw2p4DIiIir9KOeg7YOCAiItKjHTUOOKxAREREGuw5ICIi0qMdveeAjQMiIiIdmMroBWzZwTBa9KUvGdz4vkv6yJp/0Tvls5cFv5sviv99aqIoPuiIvAnrK5217XJ5vmRgnnC2wV6ydKqEjRWieAA4cFusKD7qZ3kTX5qaaA+SHYujfUNE8QBgE2boKZP8nFLCu0j2WFnqIwAEH5bF1+pMg67nzj0k8Igs3v+Y/B5S1UF2bYRky87bqnj5yLI9RTYLp+k3+Syc73X9SHdsRYUTPcQluInPHBAREVF75bU9B0RERF7FqdzrZjp1/TaCjQMiIiI9OKxARERE7RV7DoiIiHTxsOcAbafngI0DIiIiPdrRsILXNg4cIXVQ/vpSf4am7RVvf8vqi0TxhTfI0ncAoDpalpoY8ats+0f7yNPtgg7LRpKcsuw8AEBpD9kFkPiNLMUrd1SwKB4Akr+UzS7pUy4/3r9Pkc2AGLVTtn1bmDzN0Cm8wt053gEFsuOtTPIyqhJk+x58WFYnhxv7XRMjnEXVXz4DYlVn2YylJpusDGuU/D+r2DWyL+u2J1aJy/jBmqA7ttrqAFAsLoPOzmsbB0RERF7FqeDR0ACzFYiIiM4zynly8WT9NoLZCkRERKTBngMiIiI9+EAiERERafCZAyIiItJgz0Hr8z3mA6NFX/X8jPLZznz6nhDFW4/IZ8QzCKtlD5GlRsV/78bDLQaHKLyigzzvzKdKth+2UFkZ1ig39lt4UR7rHSYuInKXeBWR8N9laW0AUNFRdolXdZCnSwYUyb7b4oHyG2TEDtnjUSGH7aL46ljZLI4AkPT5cVG80yJPZawNkqXtVibKrg1llh+LE91lx+K9wwPFZRTujtEd67RaATTzxdcOeW3jgIiIyKsoeNhz0GQ1aXZsHBAREenRjoYVmMpIREREGuw5ICIi0sPpBODBi4ycbeclSGwcEBER6cFhBSIiImqvRD0HTz31FJ5++mnNZ6mpqdi79+SsiFarFQ899BCWLl0Km82GUaNG4dVXX0VsbKy4Ys6kGiBAXyurqk6ehlT1h2wGPUTKUqNOktXLt0rWqjTZ5a1Qa5g0bVBehvm4LB2uIkkW7yP8ngCgopNFFB+cKz/eJ7rJZqurDZLtd0lfeVv+ghWy2SVjPy8Ql3H08hRRfPBB+X7UBcjia4Nl57lfmWz7AFDRVXYPOd5DnhZslmVcI2a7LFU5Z5S8TnEDZedIyXr9MyzWS33vsO7YOqcN+qM9xJ6DM+vZsycKCgpcy6ZNm1y/mz59Oj777DOsWLECGzZsQH5+Pq677romrTAREVGrcCrPlzZC/MyBj48P4uLiGnxeVlaGt99+Gx9++CGGDx8OAFi8eDF69OiBzZs3Y+BA+YswiIiIqOWJew7279+PhIQEdO7cGTfffDNycnIAANu2bUNtbS1Gjhzpiu3evTuSkpKQlZV1xu3ZbDaUl5drFiIiIm+jlNPjpa0QNQ4GDBiAJUuWYM2aNVi0aBGys7MxePBgVFRUoLCwEH5+fggLC9OsExsbi8LCwjNuc86cOQgNDXUtiYmJbu0IERFRs1IeDim0oWcORMMKY8aMcf07LS0NAwYMQHJyMpYvXw5/f3+3KjBz5kzMmDHD9XN5eTkbCERE5H2Uh7MytqHGgUepjGFhYejWrRsOHDiAuLg42O12lJaWamKKiooafUahntlsRkhIiGYhIiKi1uPRS5AqKytx8OBB/Nd//RcyMjLg6+uLdevWYcKECQCAffv2IScnB5mZmeJtO4+ZgWp9qYDH44V5TgAsJbJ2kU+OLBUOAEq7y1qJfmXC2QzD5YfPFi6rU1CuuAgYa2VlHL9IFh+1Qz5zoEE41GfOFeaQAai59MyN4MYEFsj2O+HDGlE8AJhOVIvXkTqWJtsPHzeqlLRGttL+W2RpxD4V8r/oTDbZPSToiDspuLJzvWyALJXRfFA+U6T9bdl5bh8nP+C//U1/Gc4aK3CPuAj3OJ3ym8mp2tAzB6L/XR5++GGMHTsWycnJyM/Px6xZs2AymXDjjTciNDQUd9xxB2bMmIGIiAiEhITg/vvvR2ZmJjMViIio7WtHwwqixsGRI0dw44034tixY4iOjsall16KzZs3Izo6GgAwf/58GI1GTJgwQfMSJCIiImo7RI2DpUuXnvX3FosFCxcuxMKFCz2qFBERkbdRTieUB8MKbSmVkRMvERER6dGOhhU48RIRERFpsOeAiIhID6cCDO2j58BrGwcxFxyDKVBfOlKEWZ4qc0yYVhRYVCsuQ5lks/T5VDd/Sph0HUup/GQu6i+L7/qPClF8XZh8Fs4jw2TrHO0tn0nUJMw0FM/CeSBPVgCA6v6dRfFHhkeJyzDWydLt3EmPPXCjLJW44zey7Re6kVAVsUd2MVlj5OetNUJ2iw5ZKyvDr1KW+ggAZZ1lMzkGBNjEZdSZ9XdqO0zy7btNKQCepDK2ncYBhxWIiIhIw2t7DoiIiLyJciooD4YVVBvqOWDjgIiISA/lhGfDCkxlJCIiOq+0p54DPnNAREREGl7Xc1DfsnJU638CtdbHLi7HYbeK4uvq5NkKDpusC8lgl7UqDfLdhkPYHKwTTqIEAE7ZV4s6h+xp47o6eZ0cVtk6bnX+WWVP7TtqZaXUOeUHvK5WdjCcVtmT6CcLEe63G+ets0b4XQkvV6dV/ndSXZ3wHuLGteSwyW7RDuE9pE54DgKAwyY7RwyCe7mrDMGNqv7/ipb4q7xO2TwaGqiD/P+R1mJQXtbPceTIESQmJrZ2NYiIqA3Jzc1Fx44dm2XbVqsVKSkpKCws9HhbcXFxyM7OhsUin+m3JXld48DpdCI/Px/BwcEwGLR/kZSXlyMxMRG5ubkICQlppRq2PO4397s94H5zv92hlEJFRQUSEhJgNDbfSLnVaoXd7ka312n8/Py8vmEAeOGwgtFoPGfrLyQkpF1dRPW43+0L97t94X67LzQ0tIlqc2YWi6VN/KfeVPhAIhEREWmwcUBEREQabapxYDabMWvWLJjN8neUt2Xcb+53e8D95n6T9/C6BxKJiIiodbWpngMiIiJqfmwcEBERkQYbB0RERKTBxgERERFpsHFAREREGm2mcbBw4UJ06tQJFosFAwYMwI8//tjaVWpWTz31FAwGg2bp3r17a1eryW3cuBFjx45FQkICDAYDVq1apfm9UgpPPvkk4uPj4e/vj5EjR2L//v2tU9kmdK79vvXWWxsc/9GjR7dOZZvQnDlz0K9fPwQHByMmJgbjx4/Hvn37NDFWqxVTp05FZGQkgoKCMGHCBBQVFbVSjZuGnv0eOnRog2N+zz33tFKNm8aiRYuQlpbmegtiZmYmvvzyS9fvz8djfb5oE42DZcuWYcaMGZg1axa2b9+O9PR0jBo1CsXFxa1dtWbVs2dPFBQUuJZNmza1dpWaXFVVFdLT07Fw4cJGf//888/j73//O1577TVs2bIFgYGBGDVqFKxW4dSPXuZc+w0Ao0eP1hz/jz76qAVr2Dw2bNiAqVOnYvPmzVi7di1qa2txxRVXoKqqyhUzffp0fPbZZ1ixYgU2bNiA/Px8XHfdda1Ya8/p2W8AuOuuuzTH/Pnnn2+lGjeNjh074tlnn8W2bduwdetWDB8+HOPGjcOvv/4K4Pw81ucN1Qb0799fTZ061fWzw+FQCQkJas6cOa1Yq+Y1a9YslZ6e3trVaFEA1MqVK10/O51OFRcXp+bOnev6rLS0VJnNZvXRRx+1Qg2bx+n7rZRSU6ZMUePGjWuV+rSk4uJiBUBt2LBBKXXy+Pr6+qoVK1a4Yn777TcFQGVlZbVWNZvc6futlFJDhgxRDzzwQOtVqoWEh4ert956q90c67bK63sO7HY7tm3bhpEjR7o+MxqNGDlyJLKyslqxZs1v//79SEhIQOfOnXHzzTcjJyentavUorKzs1FYWKg59qGhoRgwYMB5f+wBYP369YiJiUFqairuvfdeHDt2rLWr1OTKysoAABEREQCAbdu2oba2VnPMu3fvjqSkpPPqmJ++3/U++OADREVFoVevXpg5cyaqq6tbo3rNwuFwYOnSpaiqqkJmZma7OdZtldfNyni6o0ePwuFwIDY2VvN5bGws9u7d20q1an4DBgzAkiVLkJqaioKCAjz99NMYPHgwdu/ejeDg4NauXouonzu9sWPfFPOqe7PRo0fjuuuuQ0pKCg4ePIjHH38cY8aMQVZWFkwmU2tXr0k4nU48+OCDuOSSS9CrVy8AJ4+5n58fwsLCNLHn0zFvbL8B4KabbkJycjISEhLw888/49FHH8W+ffvwz3/+sxVr67lffvkFmZmZsFqtCAoKwsqVK3HhhRdi586d5/2xbsu8vnHQXo0ZM8b177S0NAwYMADJyclYvnw57rjjjlasGbWEG264wfXviy66CGlpabjggguwfv16jBgxohVr1nSmTp2K3bt3n5fP0pzNmfb77rvvdv37oosuQnx8PEaMGIGDBw/iggsuaOlqNpnU1FTs3LkTZWVl+PjjjzFlyhRs2LChtatF5+D1wwpRUVEwmUwNnmAtKipCXFxcK9Wq5YWFhaFbt244cOBAa1elxdQf3/Z+7AGgc+fOiIqKOm+O/7Rp07B69Wp899136Nixo+vzuLg42O12lJaWauLPl2N+pv1uzIABAwCgzR9zPz8/dOnSBRkZGZgzZw7S09Px8ssvn/fHuq3z+saBn58fMjIysG7dOtdnTqcT69atQ2ZmZivWrGVVVlbi4MGDiI+Pb+2qtJiUlBTExcVpjn15eTm2bNnSro49ABw5cgTHjh1r88dfKYVp06Zh5cqV+Pbbb5GSkqL5fUZGBnx9fTXHfN++fcjJyWnTx/xc+92YnTt3AkCbP+anczqdsNls5+2xPm+09hOReixdulSZzWa1ZMkStWfPHnX33XersLAwVVhY2NpVazYPPfSQWr9+vcrOzlbff/+9GjlypIqKilLFxcWtXbUmVVFRoXbs2KF27NihAKgXX3xR7dixQx0+fFgppdSzzz6rwsLC1Keffqp+/vlnNW7cOJWSkqJqampaueaeOdt+V1RUqIcfflhlZWWp7Oxs9c0336iLL75Yde3aVVmt1tauukfuvfdeFRoaqtavX68KCgpcS3V1tSvmnnvuUUlJSerbb79VW7duVZmZmSozM7MVa+25c+33gQMH1OzZs9XWrVtVdna2+vTTT1Xnzp3VZZdd1so198xjjz2mNmzYoLKzs9XPP/+sHnvsMWUwGNTXX3+tlDo/j/X5ok00DpRSasGCBSopKUn5+fmp/v37q82bN7d2lZrVpEmTVHx8vPLz81MdOnRQkyZNUgcOHGjtajW57777TgFosEyZMkUpdTKd8YknnlCxsbHKbDarESNGqH379rVupZvA2fa7urpaXXHFFSo6Olr5+vqq5ORkddddd50XjeHG9hmAWrx4sSumpqZG3XfffSo8PFwFBASoa6+9VhUUFLRepZvAufY7JydHXXbZZSoiIkKZzWbVpUsX9ec//1mVlZW1bsU9dPvtt6vk5GTl5+enoqOj1YgRI1wNA6XOz2N9vjAopVTL9VMQERGRt/P6Zw6IiIioZbFxQERERBpsHBAREZEGGwdERESkwcYBERERabBxQERERBpsHBAREZEGGwdERESkwcYBERERabBxQERERBpsHBAREZHG/wNOgHJXF3BslgAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 640x480 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"from matplotlib import pyplot as plt\n",
"\n",
"data01_raw = uvd_raw.get_data(*second_antpair, 'xx')\n",
"data01_vvcheby = uvd_vvcheby.get_data(*second_antpair, 'xx')\n",
"data01_vvnocheby = uvd_vvnocheby.get_data(*second_antpair, 'xx')\n",
"\n",
"plt.imshow(np.abs(data01_raw - data01_vvnocheby) / np.abs(data01_vvnocheby), aspect='auto')\n",
"plt.colorbar()\n",
"plt.title(\"abs(raw - vvnocheby) / abs(vvnocheby)\")\n",
"\n",
"fig = plt.figure()\n",
"\n",
"plt.imshow(np.abs(data01_vvcheby - data01_vvnocheby) / np.abs(data01_vvnocheby), aspect='auto')\n",
"plt.colorbar()\n",
"plt.title(\"abs(vvcheby - vvnocheby) / abs(vvnocheby)\")"
]
}
],
"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.10.12"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment