Skip to content

Instantly share code, notes, and snippets.

@wtbarnes
Created July 16, 2025 19:38
Show Gist options
  • Save wtbarnes/7e6889ad2656c83d5f918db09ed919cf to your computer and use it in GitHub Desktop.
Save wtbarnes/7e6889ad2656c83d5f918db09ed919cf to your computer and use it in GitHub Desktop.
Comparing free-free radiative loss between fiasco and CHIANTI IDL across different versions of CHIANTI.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"id": "382bc1ee-0333-42ec-9d65-e421a8573c45",
"metadata": {},
"source": [
"# Free-free Radiative Loss Testing\n",
"\n",
"Trying to understand differences in free-free radiative losses in IDL and fiasco across different versions of CHIANTI."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "2b578b1c-e4e0-409b-9968-c48b072818ab",
"metadata": {},
"outputs": [],
"source": [
"import astropy.units as u\n",
"import fiasco\n",
"import fiasco.io\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"\n",
"from fiasco.util import build_hdf5_dbase, get_chianti_catalog\n",
"from fiasco.tests.idl.helpers import setup_idl_environment"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "e29a6772-9ed8-4c77-9682-476c62502237",
"metadata": {},
"outputs": [],
"source": [
"ascii_dbase_v8 = '/Users/wtbarnes/Documents/projects/chianti/database/CHIANTI_v8.0.7/dbase/'\n",
"ascii_dbase_v9 = '/Users/wtbarnes/Documents/projects/chianti/database/CHIANTI_v9.0.1/dbase/'\n",
"ascii_dbase_v11 = '/Users/wtbarnes/Documents/projects/chianti/database/CHIANTI_v11.0.2/dbase/'"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "056fb86b-5a9f-4056-b094-d632c6ccd3db",
"metadata": {},
"outputs": [],
"source": [
"ff_rad_loss_script = \"\"\"\n",
"abund_file = FILEPATH('{{abundance}}.abund', ROOT_DIR=!xuvtop, SUBDIR='abundance')\n",
"ioneq_file = FILEPATH('{{ionization_fraction}}.ioneq', ROOT_DIR=!xuvtop, SUBDIR='ioneq')\n",
"ff_rad_loss, temperature, free_free_radiative_loss, abund_file=abund_file, ioneq_file=ioneq_file\n",
"\"\"\""
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "67c5346f-6cc1-4563-abdb-4704a46e593d",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"HEAD is now at 6489b99 Adding files for CHIANTI v8.0.9\n",
"v8.0.9\n"
]
}
],
"source": [
"!cd /Users/wtbarnes/Documents/projects/chianti/chianti-idl-github/ && git checkout v8.0.9 && git describe --exact-match --tags"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "d91e7b03-719c-46e2-a0df-7599eda48fa9",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"IDL 8.9.0 (darwin x86_64 m64).\n",
"(c) 2023, L3Harris Geospatial Solutions, Inc.\n",
"\n",
"Licensed for use by: ASA - Goddard Space Flight Center - Code 600, Science &\n",
" Exploration Directorate, Greenbelt - Greason\n",
"License: 881395\n",
"% Compiled module: HISSW_PROCEDURE.\n",
"% Compiled module: FF_RAD_LOSS.\n",
"% Compiled module: SETUP_ELEMENTS.\n",
"% Compiled module: READ_ABUND.\n",
"% Compiled module: READ_IONEQ.\n",
"% Compiled module: TRIM.\n",
"% Compiled module: READ_IP.\n",
"% Compiled module: READ_GFFINT.\n",
"\n",
"\n"
]
}
],
"source": [
"idl_res_v8 = setup_idl_environment(\n",
" ascii_dbase_v8,\n",
" '/Users/wtbarnes/Documents/projects/chianti/chianti-idl-github/',\n",
" '/Applications/harris/idl/',\n",
").run(\n",
" ff_rad_loss_script,\n",
" args={'abundance': 'sun_coronal_1992_feldman_ext',\n",
" 'ionization_fraction': 'chianti'},\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "2203a4cc-4f5f-490b-af6b-7d907ad68e2b",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Previous HEAD position was 6489b99 Adding files for CHIANTI v8.0.9\n",
"HEAD is now at 1f6b0eb Adding files for CHIANTI v9.4.3\n",
"v9.4.3\n"
]
}
],
"source": [
"!cd /Users/wtbarnes/Documents/projects/chianti/chianti-idl-github/ && git checkout v9.4.3 && git describe --exact-match --tags"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "b41e8fd9-c4bc-4f2c-9947-98863aa24b41",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"IDL 8.9.0 (darwin x86_64 m64).\n",
"(c) 2023, L3Harris Geospatial Solutions, Inc.\n",
"\n",
"Licensed for use by: ASA - Goddard Space Flight Center - Code 600, Science &\n",
" Exploration Directorate, Greenbelt - Greason\n",
"License: 881395\n",
"% Compiled module: HISSW_PROCEDURE.\n",
"% Compiled module: FF_RAD_LOSS.\n",
"% Compiled module: READ_ABUND.\n",
"% Compiled module: READ_IONEQ.\n",
"% Compiled module: TRIM.\n",
"% Compiled module: FF_GAUNT_SUTHERLAND.\n",
"% Compiled module: READ_GFFINT.\n",
"\n",
"\n"
]
}
],
"source": [
"idl_res_v9 = setup_idl_environment(\n",
" ascii_dbase_v9,\n",
" '/Users/wtbarnes/Documents/projects/chianti/chianti-idl-github/',\n",
" '/Applications/harris/idl/',\n",
").run(\n",
" ff_rad_loss_script,\n",
" args={'abundance': 'sun_coronal_1992_feldman_ext',\n",
" 'ionization_fraction': 'chianti'},\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "296d10e5-1071-423e-9019-95034419701c",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Previous HEAD position was 1f6b0eb Adding files for CHIANTI v9.4.3\n",
"HEAD is now at 5efffae updated VERSION file\n",
"v11.0.1\n"
]
}
],
"source": [
"!cd /Users/wtbarnes/Documents/projects/chianti/chianti-idl-github/ && git checkout v11.0.1 && git describe --exact-match --tags"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "c91ede7b-0892-4caf-9976-59f6f032cc96",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"IDL 8.9.0 (darwin x86_64 m64).\n",
"(c) 2023, L3Harris Geospatial Solutions, Inc.\n",
"\n",
"Licensed for use by: ASA - Goddard Space Flight Center - Code 600, Science &\n",
" Exploration Directorate, Greenbelt - Greason\n",
"License: 881395\n",
"% Compiled module: HISSW_PROCEDURE.\n",
"% Compiled module: FF_RAD_LOSS.\n",
"% Compiled module: READ_ABUND.\n",
"% Compiled module: READ_IONEQ.\n",
"% Compiled module: TRIM.\n",
"% Compiled module: FF_GAUNT_SUTHERLAND.\n",
"% Compiled module: READ_GFFINT.\n",
"% Program caused arithmetic error: Floating underflow\n",
"\n",
"\n"
]
}
],
"source": [
"idl_res_v11 = setup_idl_environment(\n",
" ascii_dbase_v11,\n",
" '/Users/wtbarnes/Documents/projects/chianti/chianti-idl-github/',\n",
" '/Applications/harris/idl/',\n",
").run(\n",
" ff_rad_loss_script,\n",
" args={'abundance': 'archive/sun_coronal_1992_feldman_ext',\n",
" 'ionization_fraction': 'chianti'},\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "f8008de0-cefc-4246-9895-1469f0d0bd26",
"metadata": {},
"outputs": [],
"source": [
"hdf5_dbase_v8 = '/Users/wtbarnes/Documents/projects/chianti/database/CHIANTI_v8.0.7/dbase.h5'\n",
"hdf5_dbase_v9 = '/Users/wtbarnes/Documents/projects/chianti/database/CHIANTI_v9.0.1/dbase.h5'\n",
"hdf5_dbase_v11 = '/Users/wtbarnes/Documents/projects/chianti/database/CHIANTI_v11.0.2/dbase.h5'"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "55b46bbd-ad17-4291-8fb8-aca6e60825f9",
"metadata": {},
"outputs": [],
"source": [
"col_v8 = fiasco.IonCollection(*[fiasco.Ion(ion, idl_res_v8['temperature']*u.K,\n",
" abundace='sun_coronal_1992_feldman_ext',\n",
" ionization_fraction='chianti',\n",
" hdf5_dbase_root=hdf5_dbase_v8)\n",
" for ion in fiasco.list_ions(hdf5_dbase_root=hdf5_dbase_v8)])\n",
"col_v9 = fiasco.IonCollection(*[fiasco.Ion(ion, idl_res_v9['temperature']*u.K,\n",
" abundace='sun_coronal_1992_feldman_ext',\n",
" ionization_fraction='chianti',\n",
" hdf5_dbase_root=hdf5_dbase_v9)\n",
" for ion in fiasco.list_ions(hdf5_dbase_root=hdf5_dbase_v9)])\n",
"col_v11 = fiasco.IonCollection(*[fiasco.Ion(ion, idl_res_v11['temperature']*u.K,\n",
" abundace='archive/sun_coronal_1992_feldman_ext',\n",
" ionization_fraction='chianti',\n",
" hdf5_dbase_root=hdf5_dbase_v11)\n",
" for ion in fiasco.list_ions(hdf5_dbase_root=hdf5_dbase_v11)])"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "2d7b177a-1290-467f-87d0-b8517d51d946",
"metadata": {},
"outputs": [],
"source": [
"ff_rad_loss_v8 = col_v8.free_free_radiative_loss()\n",
"ff_rad_loss_v9 = col_v9.free_free_radiative_loss()\n",
"ff_rad_loss_v11 = col_v11.free_free_radiative_loss()"
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "1a73775e-d83b-407f-aaaf-639baf6616d5",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x127a1b020>"
]
},
"execution_count": 16,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAApAAAAGhCAYAAADWe4TFAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAAPYQAAD2EBqD+naQAA5xhJREFUeJzsnQWcFPUbxp/bve6ju7u7SzollG5EBAv5A4pYhCKKCoqBgaLSCAiCIB2S0t3dcd2x+/88v2WWu+Oavd29u/eL4+7OzM3+dmd25p03ntfBaDQaIQiCIAiCIAhpRJfWFQVBEARBEASBiAEpCIIgCIIgpAsxIAVBEARBEIR0IQakIAiCIAiCkC7EgBQEQRAEQRDShRiQgiAIgiAIQroQA1IQBEEQBEFIF47pW13QiI2NxeHDh5E/f37odGKHC4IgCEJWwGAw4O7du6hZsyYcHcUMyijyzWUQGo/16tWz9TAEQRAEQcgA+/fvR926dW09jCyLGJAZhJ5H7QAsWLCgrYcjCIIgCEIauH37tnIAaddxIWOIAZlBtLA1jcciRYrYejiCIAiCIKQDST97OuTbEwRBEARBENKFGJCCIAiCIAhCuhADUhAEQRAEQUgX2T4H8vr16xg4cCDu3bunyvXfe+899OzZM9VlliIuLg4xMTEW3aaQtXBycoJer7f1MARBEATBYmR7A5KG4axZs1CjRg3cuXMHtWvXRseOHeHh4ZHisqfFaDSqbQYGBlrkcwhZG19fXxQoUAAODg62HoogCIIgPDXZ3oBklbQms8MLeJ48eeDv76+MxJSWPS2a8ZgvXz64u7uL4ZBD4Y1EeHi48nITkXwSBEEQsgM2NyB37NiBGTNm4ODBg0qbaeXKlejWrVuCdb755hu1Do2y6tWrY/bs2RkS8eZ7MKRctGjRdC1LL9yOZjzmzp37qbcnZG3c3NzUI41IHhMSzhYEQRCyOjY3IMPCwpRROGzYMPTo0eOJ5UuWLMH//vc/zJkzB/Xr11ch53bt2uHs2bPqYkwYgmZrwcRs2LABhQoVUs/pWRw0aBB+/PHHJ9ZLaZlGVFSUmjRCQkKSXVfLeaTnURDiHws8NsSAFARBELI6NjcgO3TooKbk+OKLL/Diiy9i6NCh6jUNybVr1+Lnn3/GhAkT1LwjR46k+B40/OjV5PqNGjVK87L4fPzxx5g8eXK6PpuErQUNORYEQRCE7IRdy/hER0er0HLr1q0TKMfz9Z49e9KcgzZkyBC0bNlSVVyndVli3n77bQQFBZmnU6dOZfBTCYIgCIIgZG3s2oB88OCByidM3K+Sr5kPmRZ27dqlwuB//vmnCnVzOn78eKrLEuPi4gJvb2/z5OXlZYFPKAiCIAiCkPWweQg7s2nSpAkMBkO6lyUHC3o40TsqCIIgCIKQE7FrDyRldVhwcPfu3QTz+ZqyO7bglVdeUeHrbdu2IbvBcH78Cni+Zu4eJ4ph0/Pbpk0blX+a2PAuUaKEKnDKTP777z+0atVKaSr6+fmpYqqjR49m6nsKgmAbmGIUGhqKq1ev4ty5cyrqREkszheyP6HhYfhp/1rcCPK39VCErGhAOjs7K3HvzZs3m+fRcOHrhg0b2nRsOYX27dsreaUrV65g3bp1eOaZZzB69Gh07tw5ycr3zIIXEo6lWLFi2LdvH/7991+VRkAjUjr9CELWhQYhVTUCAgLM8/744w+4urqq3zhvTsuXL680VKnRy+sCnQuUcmNx5eeff47169er85SQ9Y3GDxdPwNCvG6DjonpYv/9DzN0519bDEuw1hE3D4MKFC+bXly9fVlXVuXLlUsYCJXwGDx6MOnXqqBMGvVyU/tGqsq1NRkLYPEFGxMTBFrg56Z+qApi5n5q3t3DhwqhVqxYaNGigPIHz5s3D8OHDn3qM9C7wAnH69GlUqFDBPH/mzJn4+uuvcfHiRZw5c0bJLU2ZMsWs1fnBBx+gWrVqykNRpkyZpx6HIAjWh92/aAD++uuvSk6N0FDUzrE0GHke4rWC51LeuD58+FBNjEpo8DzHcwQjFITeSpFSs39i4+Kw7MS/2LhzMi653cVDRx2gSgx0KIkouPgfsvUQBXs1IA8cOKC8Who0GAmNRhoovXv3xv379/H++++rEAYLXXiySVxYY80QNqcbN26kWXScxmOl9/+BLTg1pR3cnS27m1m1Tu3OFStWWMSALFeunLpBWLBgAaZOnWqez9f9+vVTz2lgUpR97ty5mDhxoiqu4vOKFSsqD4UgCPYPHQT8Xb/55pvKMCSVK1dWKUHxU5WaNWumoh78zdOYpHHI6BOdB1TBoKHIG8+TJ0+aJxqZmvFIKA9H8f4ffvgBTZs2tcnnFZImNiYGn6+ehoPhQTgfvRexjiGoCldlPHrGGVAh3BmVIvR4oWAF5Gr0oq2HK9irAdmiRYtUc1peffVVNQn2Az2Fx44ds9j2+vfvr7yNmgHJiwMlnObPn69eM5TFiwxzNLV1ypYti3/++Uf1NBcEwX65du0a3nvvPfz+++/qfM8bv+eee04t4/xp06aZDUpCozFxS1lKuPE8wKlIkSIq+vD888+bl8dPZeFznj9ocMbPl9+5c6cKddO4FCUN6xuNX675FKdursMV1wDcc9KhZqA3Yn1D4GEwoGBsJKrfc8GL+UoiV7NeQPn2gIvsI3tGrrxWCiPTE2ir984MeBGwpDh2nz59MG7cOOzdu1eFyOmlYLhcC2lHRETghRdeQOPGjbFo0SLlgfzss8/QqVMnFcbS2gUKgmA/sKUrmzB8+eWX5k5e1PFlDqOGj4+PRd6LhX7xnzNKxFa5vNGMnxbDdrn0VnIc3bt3x7PPPou8efNaZAzCk0bjN3/PwrHrq3HV1R93nR6Hp90MBvgZYjDr7kM0zF8H7s17ARU6AW5+th62kEbEgLRCDiQNLUuHkW0N8xVLlixpse3RS8DQ+MKFC5UBycdRo0aZl/M1Q1oUkKcnQpvHauxVq1YpA1QQBPuAN5jfffed8i4y3KxFmz799FPUrVvXKmNgOJvGYXyYenPixAmcP39edTTjxPMJb0xpTDLCYcnzWk6EqQb/nD+C347/iTv3duKBxwOz0ehqMKBiuBPKh+sxzDcPCjbpCVTsCniKAZ8VyV5WjZ3mQGY3tmzZogTXx4wZY9HtMozN3Ki+ffvi0qVLCYxCJsTzRB/f66m9Tq+WpyAImQfzDikBRtUGLceRhiPDxlZr6RnuD9w7BQReB4JvAsG31PRB4Qd4/428OHXDESsPP8DK40E4dDNahbY5MQe/WpnCeLZ6XnTp8wLq9HjZdMOqzjFGQCd97JPzNP606XscvLgMRkM49vlEmhZ4AMWjAL9YHcqEO2GIuzeK13oeqNwN8C5k62ELT4kYkEKKMOzE4iWGjJnkzgImhqQo46NVTGrcvHnzib7kxYsXV17CtNCjRw/ldeTEwqpChR6fYKg/OX78eGW8v/baa8ponD59usp/jF+EJQiC7WBOMgsgea6gDA8NR/6erZanfGET8NcbQND1ZFehCVvZGahcH3i3viuuBTnjzxIf4s9Vq1TI+9iFm2r6cPlrKFDgI6xevRp1CzsCP3cA8lcGClQFCtUASjQFcpViiAk5lZ82fId955bghssD3HB2UJ5GhqZ9o5xROzYY7cLC0dC9BHyr9QQqdwdyiXc3OyEGpJAiNBipv8YLAA1BhoC++uordZHQQskazEnkFB8mzQ8YMEBVStMrMWnSpGTfi0ntXbp0wdKlS5VYeXyYC/nXX39h8uTJSgOU712zZk3z+ARBsB1M6Xn77bfxxRdfqNdVqlRRucp8zDSO/wEcWQDUGmQyTgjz5zTj0beYycDzLvxoKgR45AWc3EyTo6uaihnj8HqBqnh99GglDfT3LzPw19//YP1/51U7XZVDeWk1EBuB5Rv3ICByN7qUc0R+Tx3gUwwo1Rwo1QIo0ypH5O/tuHwS89e9jZtOl3DtkdFIs9zJaESFcD3KRrhghKMDClcfAVTpAeR5nIMqZC8cjCLrn+EcSOoTXr9+XVUExicyMlLJVTCXhnfhOR2GnynHwZAW86ByInJMCNkZarJSo5VQMYOeR4sXtoXeA9xyAfpHfo9Nk4B/ZwLV+wLd55jmxUYBNw6YPIVujyV9MgLP8ex0pXI2DXGA/yU0bNUZe4+dw3cDK2NkmTuAIQYGoxE6eiF1TkDZtkDV54ByHQDn7KNBOX/rr9h29wIOhe1HjOMtVAl2xwnvcDjSaIxgeNoZAx0MKFe9G1DlOSBfJbv2zGopaEldv1OC1/4ZM2aoqFz16tUxe/ZspU+dHMuWLVN5wMzf543IJ598onRPk2LkyJH4/vvvVaHXG2+8gayAeCDTieRApp+tW7eqApmcajwKQnaH6SXMIWSXqq5du1puw8w9PP8PcHAecH4j0G8pULa1aVnlHoBHPqBsm8frO7oAJRpb5K0pK2Qu+NHpYcxdBl16D4bBdRW6fvInkNsbuLoH382agTkrd6JX+Sj0e/AXSp9dCzh7AjUHAB0+QVZl4fb52H5iHm453cEVFwfUCfRCjG+IMhp9EI6uD1zQ3xiNilW7mDyNhWrZtdH4tCxZskTlyM6ZMwf169dXTU3YCY1dlPLly/fE+rt371b5/FrKF4s+WaR16NChJzzzVAagAkn8tK2sgHggM+EORrxNQmLkmBCyG2zwQDkerTDGotJecTHAiRUmD+P904/nNxsPtHwX9kTbtm2xceNG8+sGxd3Qv5IRvYeMQt5eppA+eJmlF1PzntopS/9dgi1Hf8Itp9u47PJ4X+qNRjQMckXbmFtoAXf4Vepu8jQWqcdqRmTV6/epU6dUhzUNyjtxSgoajbyhoF4xMRgMahvMyZ8wYcIT67MJCnVI16xZY55HhRE2Q6ERGr92gNtm/jBl6eh9FA+kIAiCkC2hF4VV1a+//jreeecdNc8ixmNsNHDoV2D3V0DgNdM8Zy+gzhCg1mC7zKejZ4qFNvQwbdq0CXuvRmDvVWDMxtl4bsUtVUTUrIgBDqtfA5qMAeq+AHvi4M2L+PHQnzhxaxOCXG8AnpzrAJ3RiPKROpQOd0bvmDDUqNTY5PUt0STbVKNXqlTpiVSMpPL0mc5AYXrm+WrodDqlJUppuaTgfK2zngY9ln/++af5NY3QgQMHKg8+1QqyGmJACoIgCOmChhLlev744w+MHTv26b3q9NCdXQdseEflGipY8NJgFFDnhafOZ8xMWFzIokJOzI1bvHixaoTANr00LjlVKuKLl6tGYlChA/CyAwPyz71/4p+DcxAT5499PhGmma5AqUgHOBuNKB3uit7RwahZoZXJaGShkP6xUHt2ISkPZFKwmIpKJIlbKOfPnx9nzpxJ8m94LCS1PudrMCeSBaq8EcuKiAFpBSFxQRCE7AT1WpmOQY/KUxuPd08B/7wNXNpmeu2Z3xSqZg4hq6WzEGyIoIUgKWlGMXW2Yz11IxCv3gDe27MYo88UUJ4pr8hbwMMLQPmOVskdXPvfGqz57xvc0d/ABe4yT6h8xtwRriiFICW50yTSiMIVOpqMRlaVM6c0G0PlD29vb5u898GDB1WHJnrzraaPamHEgEwnUkQjCEJOzXl0d3c396ju2bPn020wLhbY+B6wbw5gNAB6Z6DhK0DTsdmiBzJz3VhVy4r03377TeXOnTt37nGVLavIz/4NlG4FdPrMJDlkYU7evY5vV72N28ajOE+j8VF7cQejEeUiHVAq3A1DIu6jUtlmQNMepirybFQ9bimY66vX65W+aXzu3r2boNd6fDg/pfVZdEYvfrFixczL6eWkR58FOqzctnfEgBQEQRBSJCgoSHkbGeJj+79cuXI9/UaZR8cuMTQeK3QG2n6YLYWm2eubhRYvv/wyli9frvqD+3h5Afkqwnh+IxasXIfnLvwLt5bjgcavP7XXb8ORjfj7wn/YHXIAEQ7nUSnaA+cf2ePlIoBSEa7oHhGARiXrAo16AOXbZwuDPTNhRX7t2rWxefNmVUmt5S9u3rxZyVYlBfWKuTx+QQyLrTifMPeROZTx4W+M84cOHYqsgBiQgiAIQooKArxoHj58WMmVsLd1hg1IVlfHRACu3qawbedZQM1Bj6V5sjH0YPXq1evxjFbvY929Qhj4wYuYsj0aJwxT4XxsCdDly3RLEW05thXLd32BOw5Xcd7VgHpBnojwDVPLnHWhaOfviq7hgWhavDpQvwdQoZNd55XaI0w7YJ5rnTp1lPbjrFmzVJW1ZuyxMxvzKSnbQyhp1bx5c3z++eequpq5scyL/eGHH9RyaiNzio+Tk5PyUJYvXx5ZATEgBUEQhCShNA89Z9u2bVP5YmwGUKZMmYxtzP8ysPwFU1eYXr+ZDEj3XDnCeEwOo3dBJQPXtVVNOHufBh6eB+Z1Ahq9ZpIrSsEbuePEDizd+TnuOFxWRqPBXcujc4CDLhrjHgagbXgkChZtCNTrAVR8FvBIaLAIaYeyPEzjeP/991UhTI0aNVQnNK1Q5tq1awm6szVq1EhV5r/77ruYOHGiEhJnBXamdmeyMqIDmUGyow4kWw0yvKLJDPD1r7/+qp6zUoxeh2rVqilxVC6L/2Nhq8LM1q9iOICq/sePH1d5WLwb/Oijj6zXZ/cpyKrHhJCz+fHHHzFixAj1W+fFkj3pM8y1vSbjyNkDGPmvqdWgoDp1MRzqqY9VVejH1/+Khcdj8H7PWnDrMxfI/1hq5rL/Xcz5bxV2XV+PCP05ROseF1+UjjSiVLgrOocFoGWhsqZCmEpdAa+ElcBCxjvRCAmx/yuvnZHTqrDbt2+PX375RSX3MgGYFxG65infQe0zaxlvbCnGFlDUnGNCOsVX2fqJ40rcf1sQhKfnv//+M+d38UbtqYxHUqwB8PzPQOHagI9ctDVYmKRh6DIbL7y3Ef8dP48VZ/Zj7vHGKNLlRXxx6zzCY29in1ck3T6AE/MZ9YjRxSmjsWNYINrmLwE0eGQ0+jyWphGEzCLrScjbGFZgUzuKIZ2cAJPmmZPB3I5atWopV/yqVatUKGvevHkWeQ9WJlLGILGeFqsVS5curZ5TS43eT4YPGEJjbgmrG2nMh4SEWGQcgiA81r17/vnn1Y0yWxO+9dZb6d9IdBiw6hWTTI+GMm7EeEwOenrfmToDefLmwrmHBjSbG4gvpn2JDg9347BnOAqFu6NiVDRG+wfio3u3sNqYC7PqDEDblzcDwzcBDV8W41GwGuKBtAbMEogJt817O7lbXGOMfa3ZSH7FihUYPnz4U2+vXLlyKjGZ4rtTp041z+frfv36qedRUVFPhH7d3NxUaJh6WtJnWxAsA736/fv3VzldvFljGku6depC7wHznwPuHAOu7wdG7bH7Nn625lawPz5b+R5uhe9CgQ/yIXaJAYE7AzF7fzT+uuqPxoNzY5TnHTTOXQioPdwUos5tusEWBFsgv2hrQONxmo2apE+8Zco5sjAVKlTAsWPHLLY9XrCok6YZkPRK0jCkCK8mb8Cqt0WLFqlKRiYxT5kyRS27ffu2xcYhCDmdyZMnY8OGDeoGjTeJlKFJd7HM792BgMuAex5TVbEYj0ly4vpJLDm8EdsCDiLQcBzlwtxwzitOXZobDSgMx0p5sHPBZVy5G4u7X15Aq08nodGI8XBwzH5dYYSsh4SwhQzB2itLquf36dNHCafu3bvX7H1kyJyGKmnbti1mzJih8h4ZVqfXkjmRJH4xjyAIGYfixh9++KG5gKZq1arp28CdE8DP7UzGI4tkXtgAFG+UOYPNopy6fhr/+7Ef+n5XE4M39cbtywsRiCOALg5xjkFoHuiCKTdjsDbCEat6DsTp3ZtU/mlEVAyGjX4HA1uURciKsYCBhqYg2A65LbRWGJmeQFu9dyZw+vRpVVFsKZhnydA4ZQ8aNGigHkeNGvWEDteYMWOUx5H9Z2lwsrl9qVKW7+AgCDmN0NBQpa7Am0M+MiqQLq7uBhb2AaKCgHyVgQHLAe+CmTXcLMXF25fw7dpJuB19HOfcYhDl7AA4c4kDop0iMTLA1EqwtHsBOFTuAFTpARSsodKPWEO9fn0z1TeZKhQLdl3FwbNfYnWBlijbqJOtP5qQgxED0hrQU5cJYWRbsWXLFiWlQ2POkvCCxR67lAm6dOmS8komhl7PQoVM6QAMZ1OKgZ5KQRCeDrZVY9iardWYLpIuLmwGFvcDYiOBYg2BvotzvFB1YEQofjiwFn9f/BuxhgMIctKp6mkajYWiDSgT7oomYRHo6WSAY9UBJqORFepJRHYYZeHNctOmTdHn+a44c9cf9ToNUMWFjM4Igi0QA1JIERavMN8wvowPlfY7d+6slPfjQ2mdI0eOJJhXvHhx5S1MCz169FBeR07PPPOM2VDUYAibskI8mTI3a/r06Vi6dKnq8CAIwtNBTz47ZVy9ejV9eY+Xtj82Hsu1B57/Jcf2U7527zq+XvMBAsJOYY97NBwcY2gvomKkM9ydo1E6zBVNwyPRSx8NxyqdTEZjkXq0ENO0/SZNmuC/wyfw3HPPYc+ePUpKrW31wqYOP4VqZPrnE4T4iAGZTnKaDiQNxoIFCyq9RxqCrL7+6quvlIh34txD6jEm1mT8/fffMWDAACU0zrDYpEmTkn0vdrro0qWLMgp//vnnJ5ZTOoh6dDRqOQ7KCXXo0MGCn1YQcnY+M5UO0tVGjWHrRX0eGY8dTB1mHFVsNsdw2/8Ovlz1Hm6GH8I590iE63WAN1A61BNRLhFoGxaOhqExqKt3hWOVtqbqaeaFshd4BuD5eOvWreo8/MYLfYF57YAIf6DX7zm6q49gfaQTTQbJjp1oMrPTAnt+0gDMqXI7ckwI9gp7+bLNGtNH0tUYgPI8rLaODgXKtAb6LEyx9V52IiwqCp+tmoYL99fivHsEwuLdTOeNNaBcmCs6hAajiz4GuopdTEZjiaaWr0aPDAKWDkbshS34cGcMxnw0Bz6NBlv2PbIh0onGMogHUsh0eLfMApmcajwKgr2yfft21RCA0QSqGrC/b5rZ+bnJeCzZHOg9P9sbj/eDH2DB/r+x/s5e3Iw+gFIROlzyjFJiJnmU0eiCuuGxGBIXBMfKzYHK3YFSLQB9JkruuPoA/ZZiXNea+HLbKazr9QL2LI+Crv6IzHtPQXiEGJBCptOpUyc1CYJgXzRr1kylmdy6dSt9xiNhW8Jt04EWEwAnN2RH/EP8MevPybgavBvn3cJROdQNN32iAD0Q4hyLhkHuqB0RgyExgXCp2NGU01i6pXWNaUdnDJn6O/7Y1xxvN4mDbt14k2ey2TiLN5EQhPiIASkIgpBDYe4jc5TTTGzUY+OIyhJtH3eOyi4EhgZi5qrJuBK4CxfdwhDEnEZPLtEhwDkSfYJDleROzThH6Ms3NIWny7SyqRFdo1YtnLt6B+77ZgE7PgW2fojYMH84dpgmRqSQaYgBKQiCkMM4efKk6m/v6+ubvt7W8zoDZduavI7ZyDCJionBwqPb8MfZNYgN34xb1Gl8ZDT6xhlQPswZ1cONeCHqDtwrtDOFp/k92FG1ubuHB9DyHcA9F64sfhMdBn+Kr/93Ga0mLMpW+0qwH8SAFARByGEFXd27d0dISAj++usv1Yc+TZxZC9w6BARcAWoPyfIi4aHhYfhi1WRcf7AL+9ziYHQKU/OrRLsgTB+pjMZqEUYMi7gHr/KtTeHpsu0AF2VZ2i8NRuHjT5fizINt6PTeEqyIikbHD5aLESlYHDEgBUEQchDsN3/+/HklB8Pq6zRTrZfJC5mvUpY1Hmk0zlw9FZcebMMltxD4O5okdyoGe+OmLg6twyNQN9SApg8i4VOulcnTWL494OKFrMRXi9bjbvtGWLXtELpNXYnF0c+ix4er0qw3KQhpQQxIQRCEHMKxY8fw6aefqudff/112gTDqfSmea/qDEVWIyYuDl+u+xonryzGZbdgPKTRqOxBHTzjDKgQ7oxmYUHoH3APzsxlbNINqNDRVOGcRXFxccGyDXsxsHMzLNmwF72mr8FvTkPQb/Jvth6akI0QA1IQBCEHwG5SI0aMQGxsrAphs/NTmloU7ppl6i7jkQdZhfDICPy+by3W3tqLy+F7UCQqDje8whIYjZXCgeHhN+BXpjnQpAdQvmO2ar/o5OSEBX//C7cerTBv9XYMmDofcWXaYODAgbYempBNEANSEAQhB8AOWvv27YO3tzdmz56d+h88vAj8MdQkCbPrS7uvuKbR+MWfH+LSgy245BaM0uEuuOwdo65yAcY41A12RYVI4MWw6/Ar0wxozPB0R1V0kl1hm9e5K7fA9eWXMef771U3MIrF9+3b19ZDE7IBYkCmk5zWylAQhKwP+1tPnDhRPf/kk09UBXaKRIWY+lvTeCxSF2j5LuzVaPxs5RRcfrjtifC0u2sUOgdFoENEGOpHxcKFgueVGZ7unK2NxsRQJP6bb79FnMGAH3/8EQMHDoDTmVV4fvJiWw9NyOJIRm06eeWVV3Dq1Cls27YN2Q3enXbr1i3Ba+rEcWI4JH/+/GjTpo3qU20wGBL8LXtdz5o1K1PH9/rrr6N27doqvycp0WNWl3LMVatWVXfZ8T+LIORU2K325ZdfRlhYGJo0aaLC2CnC3/aKl4D7ZwDPAqYey3bUZSYqNhZzD2xG5wVj0efXelgWvQYHvEKV8ehhMKB2iCMG3nPC7zfv42O/GmjWZgZcxp0HBq4Aag3KUcZjfCNyzpw5GNKrC+LiDOj74RKs+uFjWw9LyOKIB1JIkfbt2+OXX35R+VN3797F+vXrMXr0aPzxxx9YvXp1+nrnWoBhw4apMByLARLDMbq5uSlDc/ny5VYdlyDYK0uWLMHff/8NZ2dn/PDDD8qYSJHtnwBn1wJ6Z6DPAruouA4KC8XMP6fguv+/OOKsQ7RrkJpfNcYFnk4RKM+cxggHvBB2A7nLNAUa0dPYKUcai8nB/f7TwpWICayPBRsOouerH2Bl4WrSJUzIMGJACilCb1+BAgXUc4a9atWqhQYNGqBVq1aqh+7w4cOf+j3OnTuH8uXL4/Tp06hQoYJ5/syZM1Wl6MWLF9Xrr776Sj3ev38/SQPSw8MD3333nXq+a9cuBAYGPvXYBCEr8/DhQ3VDRd555x1UrFgx5T84twHYPt30vMuXQJE0akRmAoFhIfhi5WRcC/wXl+NJ7lQL9MJlJwOeCQ9H/dAgNH4QidzMaWzUVYzGNOREzlu7F7EDBqgbi48++ggd2reHTq+39dCELIgYkFYKIUXERtjkvd0c3VQI2pK0bNkS1atXx4oVKyxiQJYrV06JGS9YsEBp1Gnwdb9+/Z56+4KQUxk3bpy64apUqRImTJiQ8sqB14AVL5qe1x0O1LD+by80Kgpfrv8O568vwiW3UAQkIbnTKDIEPwfdhUvpZ4DGjyR33PysPtasCqNG7H/Om/b/vfIidAt7AvVHAuXa2npoQhZDDEgrQOOx/sL6Nnnvff32wd3J8u226ClMyguYUfr376+8jZoBSa/kwYMHMX/+fIu9hyDkJA4cOKCiBLyBZPEEQ9gp9rheOhiIDAQK1QLaTbPaOO8HBWDe7hXY+OAobsXsR/5oA+558YZbBy+2EdQkdyJuwa/sM0CTrkD5DtlKcsfaMKd98uTJwNaPgYubgau7Efr8YnhWaGHroQlZCDEghQx7VS3p2ezTp4/yluzdu1eFyOl9ZLg8fkhbEIS0w4IzGo4XLlxAo0aNUl75n3dMbQpdfYFev2Z60czdwAf4ctUU3Azdi4tuYSgb4YjbngY46IFI5zg0CHJDhUijMhp9yrUEmtJobJ+lxb3tkmbjgJsHMfP3NZj1eWvs3PwPitVqZetRCVkEMSCtFEamJ9BW750ZMF+xZMmSFtse8ywZGl+4cKEyIPk4atQoi21fEHIavMFLU4rJsWXAfz+anvf4EfAtlinjuR1Ao/ED3Ardj0vu4QjS6wDVVlqHWy6x6BEYgXYRYagTCziXbQNU6gaUa5fl2ghmKfROCO/8HX4Yvw7XAqOw7IO+GPv7QcC3qK1HJmQBxIC00ok8M8LItmLLli04fvw4xowZY9HtMoz95ptvKpHbS5cuKa+kIAjp48SJEyhWrJgSDE8VSvZQJJw0HWfxPLiHYeGYe2AD1l/5B16Ru3DJ1WjOafRleDrMGVUijHgh6iG8KrQFKnUFyrQBXJRlKVgBd9+82LB9L1a92wWvVPQHFjwPDFsveaVCqogBKaRIVFQU7ty5k0DG5+OPP0bnzp0xaNCgBOvevHkTR44cSTCvePHi8PNL24mIrdXodeT0zDPPoFChQgmWMxQXGhqqxhMREWF+LxYIaPld1OikyLu/vz9CQkLM6ySlGykI2Q3+Prp0odZfHNasWYNq1aql/AeU9Bm6Ftg7xxTOtADXHtzF7NXv4174IZx10iPMlS0EgQIxbvBzDEO5cGdUizBiaHQgvCq2Ayo+C5RpDThnn5vsrEbRclXx6o//Aj+1Ufqf0fP7QTfwDzi6eth6aIIdIwakkCI0GAsWLKgq92gIsvqacjqDBw9+Qk/us88+U1N8WO03YMAAJTROke9JkyYl+15eXl7q4rd06VIlVp4YhuO2b99ufl2zZk31ePnyZbV90rFjR9V1I/E6zNkUhOzO9evX1e+Sx3uaU0yYV9jirad638t3b+HbNZNwJ/IQLrhFIJThaS+gdqAXrjjGoXV4OBqHPkT9ID3cK3UGKj0LlHoGcHJ9qvcVLIhPEaD/MgR/2xY9Pt6IUktq4ft/TsFBJH6EZHAwypU1Q9y4cQNFixZVJ+wiRYo80RGFRg1P4K6ucoIMDw9H7ty5sW7dOrRokTOr/OSYEKwFO87wJoqe+RTzHllxTbmeDBbDPQgNx8z1s3Ht7rLHRuMjcscaUDbcBU3Do9EfYdBTn5E5jSWbAY4pVIMLNmfDvE/RYdhbMBiBKYNa4L1ftyInXb+FtCMeSCHT2bp1qyqQyanGoyBYEwrqp2g8Bl4H/hoNxISZRLerPJeu8PTPOxdje9Bp3DcchG+sAUGe0SqnkUZjuXAX1AyPwwvGcDhXZiFMV6BEU0Avl5qsQtshb+Lrc2fw8se/4P3ftqFw0cEY9uGvth6WYIdk+1817zAGDhyIe/fuqTDse++9h549e6pl7FTSunVrxMbGqokt+l588ZGQrmAx2CpL2mUJQubBvGFKXr322muptypkqLLlO8CFzSavYCrc8r+vqqdvh/+H827hKBOpxwMPo5LccTbGommgJ6pGRuMFhyg4V+5gMhqLNwJ0EvrMqoya9jNuXL+CafO3YsS031C0dAW0Gfq2rYcl2BnZPoR9+/ZtVfzBIgoWX1AbjSLVvEtnojmLRNzd3VXYp0qVKkp8l+HW1JAQtpAe5JgQMotly5ahV69eynA8dOiQylNOE6zATsbYfBAciC/+nIybIbuUTqOS3HlE/hgDWoRGokNECGo654aOBiOnovWT3Z6Q9TAaDBj4TAUs2HEe3l6e2LV7j7pGZgckhG0Zsr0HkgUgnDStwTx58qgKXRqQ7AtK45HQkKQtnc3taUEQshG8EL700kvq+dtvv52y8Xh9P5Cv0mOJnETGXkRMDH49uB0rz6+FV+hWnHWPM+s0+j0KT9eIMOAFnQFuVXoDlbsBheuI0ZhNcdDpMPfvg7jerh127Nqjokhs9KBdTwXB5r/8HTt2qMpbSrZQL/HPP/98Yp1vvvlGVdnSc1O/fn3s378/Q+/F1nj0OvLOQ4NhbJ50eRcyfvx4ZWAKgiDYOwaDQSkbBAQEqF7yH3zwQfIr+18G5j8H/NAcCLphnh0VFY1Pl3+MYV83Qbsf2+Kbs2Nwy7AJjrGu8I4zoF6wI4bf1WNtoAN+KtcBr/b7CW5jTwDtpwFF64nxmM1x8fDCytVrUK5cOVy7dg3Ptm+FsID7th6WYCfY3APJ0DENuGHDhikdwMQsWbIE//vf/zBnzhxlPM6aNQvt2rXD2bNnkS9fPrUOw9PMYUzMhg0bzFqC9DpSt5CtveLj6+uLo0ePqjA33//5559H/vz5M+3zCoIgWAKeCzdv3qyiKOwZz/7GSRIbDSx/AYgKVh5Io2d+/LD+J/x3YT6uu9zHLWeT5E69QE+cNBjQKiwczcIeokFULvhU7QhU7g4Urp3ham0ha5MrVy6sXbsWDerWwoFjp9G/bS0s330J+uSONyHHYHMDskOHDmpKji+++EIVtgwdOlS9piHJg5k6gRMmTFDzEotXJ4bh6W7duqn1k+sJS6ORhuzOnTuVEZnUNjhpUKRaEATBFhw7dkyFrLVzZPny5ZNfectU1e84xskD74VG4dTcmrjs4mDuCONiMKJihB5VoyLwtX8c3Cr1Bqr0kPC0YKZMmTJY9dMMtOozEqsO3MCECW9ixuczbT0swcbY9dmBHUUYdmaltAYTxfl6z549adoGcxoZ5qGMDKux40Ovo2YIBgUFqXB6cididl/x8fExTynKZAiCIGRi1KZ3797q/MiOUCNGjEh23Ut7FgC7v1LP3/XyxN8+D5TxqDcaUTncAd0eOGPJ/Tj8Xrw13uj1DdzGnAI6TJfwtPAEjZ97CfO+mq6ef/bFLMybN8/WQxJyugcyJR48eKByFhOHlPn6zJkzadrGrl27VBicLb20/Ep2R6lataoS2+XJVyueoQQG5ycF7/YZSo/ftk+MSEEQrM3rr7+uzn8sZmAkhrnj8Tly5Qx+/mcSomJO4NOHN9W8Rcp4dEHtYBfkiwH6xoSiZoXWJg1I0WkU0kifUW/h9J1wdQ2tW6eOKZ+WslBCjiTbnzWaNGmiks2Tol69eqmGvzVcXFzUxIIeTrz7F9IGPcAsVkqqQEoQhLSzcOFCs9G4YMEC5M2bV82/eOcGvl3zHu5EH8VZt2hEuwDfB9yHj9GAa3pnRDsA/9wOQKGyHUzh6dKtpCOMkCFYrPXGKy/Bb8e7wLqdwIitYkTmUOw6RsGKaErtMNQcH76mJI8teOWVV3Dq1Cls27YN2dHQ44WJk7Ozs8p7mTJlSpIFSklx5coV9bdpNcrtuTiBqQxubm6qYn/MmDFKx1EQbMmFCxcwcuRI9fzdd99FjfoN8Mn2FRg8uyX6rWuHDU4HcMwjBlE6B7zmH4yGkZGIhQOKFW+Mwe2+QaH/nQOe+xEo30GMRyHDMI3Mz8cLuHcGCLuHPR91QvCD27YelmAD7NoDSSOGwt+sNGQRDKE3ka9fffVVWw8vW9K+fXv88ssvqmDo77//VgYzqzu1hP3sDj08LLail4cFVxSd1wxrFisIgi1gik3//v1VznaRcoVwpMBVNF/cFEZ9FCo4eCJcp0NeajWGuaBTaBg6Rwerv3NkPmN9k9EpCBbDxQvotxiLXq6LwUuOoc2muli9/wr0jnZtUgjZzQMZGhqqPFaa14rdOvicmlOEeYeU3vn1119x+vRp1bKLSeRaVbZgWRimp3e3ePHi6rtmwdLq1avVd+7t7Y0//vgjwfoMS1OUnRc2dlkhNWvWVAZX4t7Xn332mcrbYqcfGqYxMTHmZdSyo8ySn5+fkiVhZf758+fNy5mwTcmlf/75BxUrVoSnp6cydtlpKK1MnDhRSUElhtX39LSS3bt3o3HjxujXr5/SHm3bti369u2bYe1RQXha7gcF4N35/4NPgyC4F3GB+0ueyBdzRhmPuWPjUDn2Pl69o8d6x3KY0/AVdPFwgwOMQMUuQD2TyLggWBzfYijb92OwSZFH1F1Eb/3U1iMSrIzNbxfYOvCZZ54xv9YKVQYPHqyMBlYb3r9/H++//75qRUjNx/Xr19tMq/FpciDDo5MPBescHODqpLfouu7OT797GcZ9+PChMhL79OmjvJPxZY60115eXsrIYl7ppk2bULlyZeVB1ti6dasyHvnIUBz3K/el1nucXj4ajDRWaai+9dZb6Nixo0oX0PTtwsPDlRHKBG6GUQYMGIBx48apXLC0QA8Oq+kvXryI0qVLq3knT55UkijLly9Xr+l1pKae9lkuXbqkPLGJK/gFITN5GBKML/78ADeCduOSeygCeZWu5YSSNcrAz2iER7ABP925j9oFG8KxWQ+gQifAPRew6hUg4DLgXRjo8pVoNwqZSp0uw/Df/NuofOoTOOz6BChaA6jQ0dbDEnKKAUkvVWrtAxmutpeQNT1nnLRemumh0vv/JLvsmfJ58cvQeubXtaduQkRMXJLr1i+ZC0teamh+3eSTrfAPe9KgvTK9EzIK9wlTBejxY3U6GT58uDKw6PWjMXjv3j1lXNFgJFpCPz2MiXNU6Vn8+uuvVU5rhQoVVFssbp8GpGY4smJe0+mkUcjvlx7Onj17qnn0WFIHVDP+eExonsO0QKOW3kaGqd977z3z+9AryXxPQs8jq/9ZfMXvgPmfzDuj91IQMpPgyEj8cnAL1l5cB//IfTA4hyM8JhaGewbky+eE8uHOqBrugBfyF4NXs+ceG40aMRHAgwtsQgf0+CHhMkHIJKr0fgdYGwD89yOMy4fjVueFKFw9YfRJyJ7YPIQt2Bdr1qxR4WG2jWQYmZ7CSZMmqWX0yNEIYzoBoaeOoe5mzZqlul3+HY1HDc0AJUxNcHR0TBBephHKQhYu02BoWzMeE28jrdALSQOS0EBctGiRmqfB4qhp06bh22+/xaFDh7BixQolXD916tR0vY8gpIW7Qf6Y+NtoDP66Pkb+Uhc/nnsTt43bEOUSgXLBTgj8+gaufnARrx/3wdzGY/HGq5vgNfhPoNbAJw1EJzdgyFpg0CqgRBNbfSQhJ9L+Y4QVaIA+i+6jfvO2uH3xlK1HJOQED2RW42lC2KemtEsxLB2fg++1TvO6/771OAXgaWE6wXfffafCz2wDScMuPvRC8vOz0ITha+aiJtahS4rEbdb4N8nJK6VnG6l5rxPDfEaGx2kcRkRE4Pr168pI1qBnkuFqfk5CXVDmf1Iv9J133lGhc0F4Gm7538dXqyfjVth+XHQLQzDD06orDFAx2AOVHe6hVUQMyvnWR9/cMThw/xpqv/o7EO/mKVmo51iqeaZ/BkFIgN4JcV3n4PikirgZFIUe7Zpg69FrcPXwtPXIhExEroZWlPFhTmJyU/ycRkutmxGY68hwbrFixZ4wHgnzDinA/tVXX6nvgbmqGlrOI8Xf0wOLYhgq3rdvn3ke8y7Z79zSYu1FihRB8+bNVeiaU5s2bcw91bU8y8RGouY5Ta+xKggaD0LD8fG2lRg0uy16/tkCax2247BnhDIefeIMqBfsiKH3XPBTnnKY1OJzNH39DPK/uBwb9xxVqR3xPe9PsOdbYP1EIEakpgTb4V2wJFYtXwZfVwfsvRiAV54zpQFlJ+g8YXElI3T169dPtbhy2bJlKmWL69MZwZSv+DC6x+W87jLNi0Wr8a+D9o54IIV0wYO8R48eGD9+vKpQpkGmQUOMRTcscuJ8/mjY9jE1ypYti65du6p8yO+//14V5NDDWbhwYTXf0jBkTTFcepFnzkzYz7VLly5KroeV5DxBsOCHXknOjx+CF4TUuHj7Jr77exIuR4XhousZxOljUMbBUxmNvnEGlA9zRpUoR7xQrBq8GnYHyrRSYWhKaFENgfCYS647liLkDrBpEhAXBRSqCVQz5QsLgi0o27gLln75HtqPnIKf/zmK2u+OwMsf/YjsADvasciXefi8NsyaNQvt2rVTjo74TggNKnow4sXCTbYcZeoU5QgZ/apSpYpap1y5cqo2oFSpUioixusRr6u87mg1BXaNUcgQ169f562VekxMRESE8dSpU+oxKzF48GBj165dU11v8+bN6rMvXbr0iWU//vijsWjRokadTmds3rx5stsdPXq0eTnx9/c3Dhw40Ojj42N0c3MztmvXznju3Dnz8l9++UUti8/KlSvVODS2bt2qXl++fDnF8QcEBBhdXFyM7u7uxpCQkATLYmJijJMmTTKWLl3a6Orqqj7Lyy+/rP7maciqx4SQPk7fuGIc8+NQY/+vaxnr/1zJWGVeFePwL+uox+Y/VTR+MKuEcdaMusbg5S8bjec3Go0xUQn+PjY21tiqVSvjqFGj0n6snFlnNK582Wg0GDLnQwlCOvn0xdbqXOzoqDdu377dmJWu38lRr1494yuvvGJ+HRcXZyxUqJDx448/TnL9Xr16GTt16pRgXv369Y0vvfRSsu8RFBSkxrVp0yZjVsCB/7O1EZtVcyApB8McuvheOMKuJdSzpC4ivXDZDcrosDvLrVu3Ekj12BrmZLIAJr70j72Q3Y+JnMyZO7cwZ/XbeBBzHOfdo5Sot0buWAOqh7ijr0Mo6pZ9FvrKXYGi9QFd0t7s6dOnK9F+FowdPnxYeSgEIathjItFv949sXj5n8o7R7m+9KqWZCaaigqvFYx0JW5ZnBhe7/mbpA6y1tSEMIWLbXpXrVqFxDANjB7LN954wzyPkS8qixw9ejTJ92Bq2Icffqg8kOzEZ+9ICNuKMj5ZHeYHUsKHF7mXXnrJroxHwvwSGpD2ZjwK2Y8zt2/jl6ObsO/G3wh0OAkXJwPCXVhMZuoIUzbMBTVjXTGsUms4V+0KFKqVqibj3r17VYtCMnv27JSNx1uHAY98gM/ji58g2AsOekfM/W0BTl9opIylHs92wo6du+Dm+ahazE5InGNPA09THYkPpd2Y259Yfzp//vw4c+ZMktumbnVS63N+YuUTaizz+kplkY0bN2YJ45GIASmkmU8//RQfffSRku2xx9aGTFgWhMxi3/njWLBlGh7EnUGoYzQuO+mAR/cq5UPc4GEwoFacB4ZV6wh91W5A3vJpFvIOCgpS+VK8SFEVIMVOW1GhwNJBQGQQ0H85ULSuhT6hIFgOeuxWrlyJOjWr48CR4xjVowl++edImlQ7rEVSHkhr88wzz6juezRS2XWvV69eqpAmqbxKe0MMSCHN8M4sqbszQciu7D13HAtpNBpO45xbLKJctYufDrUC3VEDd9DSvRSq1XoODhU7A34l0v0ezCKiR//KlSuqwpOFZCleZDe+DwReA3yKAfkqZPzDCUImw5SdpZ+PRdsXJ+HXjcdQZ/ZXePX10bAXWLDJzmepQY8gC9ru3r2bYP7du3efaJqhwflpWV9TPuHUoEEDVVQ6d+5cu3TSJEYMSEEQhHicu38fP+xfg3sXf8VJrweIdtOMOQcUijagdLgbGuhzoX/9ntBXfhbwerq2qmzZygpPXqAobJ+icsGlbcCBuabnXb8GXOwrJCgIiWn1wgf49PxNjPvkR4wZOw7Va9ZC06ZNkZVgulbt2rVV9zQtB9JgMKjXyXXJa9iwoVoePweS4WnOTwlul0oMWQExIK0oJC4Ign1y6NJZ/LpxCu5GBeGM1zXE6YworvdAtM4BhaMNKBXuhkbO+dC3bj/oK3cC3Pws8r5MltcuQOx2RA9EskQGA6seXazqDhfBcCHL8L+Pv8eBqyHYunWrybseF6PEx7MSLIhh0UydOnVUV7ZZs2apJhNausmgQYNUOJyyPWT06NFKc/jzzz9XrXsXL16siol++OEHtZx/y5SwZ599VuU+MoRN2+LmzZvm9r32jhiQ6SQnF9EIQnbi7M1rmPP3+7gfcwxn3aIQ6apDg2hnZTwWi4lBi5h7KBBRGX0bDIK+YgfAxbJdNSiez4sOk+dbtGiBN998M+U/2PgeEHQd8C0OtJ5s0bEIQmZCo/Gnn35C4IN7KHziG2DJbKDPIiALdfZibvL9+/fx/vvvq0KYGjVqKM1jrVDm2rVrCZpQNGrUSGk/sjBu4sSJKjTNCmxNA5IRBxbgsDUwjUe2761bty527typWv9mBUTGJ4NoBmROlPER0o8cE/bBnaBAfL58Au5G/Iez7pEJJHcKxBhQOdQDQ/OVRLUaA+BQtrWpv3QmQU8FLyzMwTp27JjqK58sFzYD83uYnrPftfS6FrIi904DP7QAYiNxq/LLKNTT5K2zp+u3kHbEAykIQrYmNCoac//bhG1nF+K68RjcHAwI9IwvueOKui5FMLTpCOjLtQYcM78Sk1WXlAwh1H5L0XiMCgFWv256Xn+kGI9C1iVfRaDzTPwxbTiGTJuOmWfC8eJ7X9p6VEIGEQNSEIRsR2RUNL5cMxPnbv2JAOdgXHTWw+BoCrZUDHaDK+JQw7EQXmw+As7l2gOO1tM0ZdBn5MiRiImJUQn5DGOnyOYpQPANU4V3q/etNUxByBxq9MM5t28RFrMLK+f/gOGjXoNDnjK2HpWQAcSAFDKdIUOGKLV+5n8IQmYaZj9v/h37Ts/FNdcHuOmsA7w1yR1PFHG8j2e8q6JFm2FwLNPGKp7G5PLB2M1p3LhxqUv2XNsH7H/US7jLl4Czh9XGKQiZxYQfN6GQrjYGFL0Oh6UDgeGb5NjOgmSdDFbBKoYeL2acKFtAXaopU6aoZP+0QB07/i3Dc1kVeoX4mUuXLq1yFatXr64SpQX75cD1q3j5+0Ho/kM1zLo5A3u8/ZXx6GIwokaoHs8H5cOMxqPw0YhjaD1wBRwrdLaZ8ajBhHq2P0tRLDg2Clj9Gk1joMYAoFQLaw5REDINnbMrhny5CY6UwLp3Sh3nRoPB1sMS0ol4INNJdpfxad++veopTR0qtgZkxTlbA2YFUVNLwIq5+fPnq44AFSpUwD///IPu3btj9+7dqFmzpq2HJzziwIXjWLp3EU6G78Z154coYHDHbXdAZzSiQqQOpaL8MLRGT5RrMMTi1dNPU0jFYhlKgKSJ038BD84CHnmBtlMze3iCYF28CwK9fkX03M5445P58Fh5DTMWbLL1qIR0IB7IdEKDiu2Ptm3blv4/jg5LfoqJTMe6EWlbNwOwlROV8pnUP2rUKLRu3RqrV69WmlWsFmUz+fgwLE0l/ZCQEFVhTGho0RNJaZL4fPbZZ0rvinIF/B7p7dMICAhQuWB+fn6qBVaHDh1w/vz5BGLLvr6+yqCrWLEiPD09lbHL3txphRWv9evXf2I+vYz0OhKGFrlex44dUapUKfUd8Dm1vATbcuXebYz7cTD6f1sTI3f2xYPAdbjm8hBGByCXLhDtAn3wXd5eWDJ4Oz4evR3lmr9qN8ZjfJ1H7VhLlarPA71+N4Wu3XNl9vAEwfoUb4QdeQbiuwMx+GzhZiz5WuSpshLigbQm0wolv6xsW6B/vF7OM8oAMeFJr1u8CTB07ePXs6oC4Q+fXG9SEJ4WNzc3PHz4UBmJbPhO7+Tzzz9vXq69Zkuo/fv3K+/Kpk2blI4Vw+AaFJCl8chHiidTU4s6Wi+++KI5fE6DkcYqDdW33npLGW401ukBJdTLoxFKI496WwMGDFB5ZAsWLEjTZ+nfv7+STrl48aIKUZOTJ08qr9Dy5cvVa3peE8vs8Dv4999/n/q7FNJPQEgIPl8xETdC9uCsewRCmdeoDisHBDnFYGCYK56t8Bwq1HsB8Mxn1/mZ9+7dU49Vq1ZN+x9WejYzhyUINqf1a7MwYds2TF95DC+MnYzKdZqiSoOWth6WkAbEgBSShBc6tmGix++115iHBQwfPlyJo9LrR2OQF0SGuWkwkrx586pHehgT9/ukZ/Hrr79W4qkMDVOZn9unAakZjrt27VLbJzQKqdNFD6emyk+P5Zw5c8zGHzt4pNmbAyijlt5Giru+99575vehV5L5nqRdu3b44osv0KxZM/U+HOOKFSsQFxdngW9VSAvRsXFYdHADtp74BYGxZ3DR1Qiojn065Ik1oHyYG+p6V8aQ5yZBn8fk9bZ36JFnWsTLL7+ceirE5R1AvkqARx5rDU8QbIeDAz5c8C8O1iqKjWeC0KNrZ+w/eQW+eez3hlAwIQakNZl4K/llDvqEr8dfSGHdRJkHbxyHpVizZo0KD9NYY0/Ofv36YdKkSWoZvYs0wqicP2HCBJUryFA3ja3U4N/ReNSgAXr8uGncp0+fhqOjY4LwMo3Q8uXLq2UaDG1rxqO2DRqx6YFeyJ9//lkZkDSS2XuYLao0vvzyS2XU0sjlRZ/vx1ZV/Bsh8zAYjPhm7UycuLIcl50icdstWp2dqka4wdspDBXCXVDZuRRe6TwJLoVMnRyyIqkaj6H3gSUDlIcVw9abdPMEIZujd/PCwtWbULteA5y/F4FB/Xriz/VbE3R2EeyPdO+dQ4cOmS/8hJWE1DJj3lh2LSyxGJQpSG5yck3Huom6YyS3XgZ45plnVBU1vYIRERHKWGT4WoNeSOYjauFrGlcpypA8QgtDa/BvaKCmh6S2kd5GSn379sXZs2fVcczCGHYiYDhdg15Uej2Z83n16lXVaooGNfMhBcvCfTd/828YNbsZuvxUFT/4/4Ld3sEoEOUCN4MBrcKN6ORbEWta/Ii5rx7C/0b8keWMx0uXLikPOo+zNBHhD/gUA3yLAXnKZfbwBMFuyFO2DlbM+xYuLs74a+MO1SdayGYeyJdeekl5n5jHw5Mj8+JYpbps2TKVo8YG40LWhcaiFs5NCuYdsmcvu2cwP5HN5TW0nMf0hntZFEOpoH379plD2My7pKFXqVIlWBK2rWKDe4auaSC3adMmSSkV5kEWLlxYeWKZH9mrVy+LjiMns/HINqza+SFuOt/FBd43Ka1GBzgajagYrkcFp7z4puUX8CpWX4W3srKBzPMlUzyYW8s0jVTJWx4YsQ0IuwfoEkUlBCGbU7vrCMyZ46wcE+zUVLtWLXTs1MnWwxIsZUCeO3dOFT8QGo0MXzKnjPlrNCazuwGZ3WV8UoO5jD169MD48ePRtm3bBH1EaYix4IS6iZxPI8zHxydNmnhdu3ZVoWMKK7MghzcpNOA439IwjM2TE/fhzJkzEyyjEXvz5k11jPOR4Xt6Smk0Cxnn3O0bWLzzaxwN2IqHxjg89I5S8x2MRpSPdECJ6LwY0GAEqtftnaWNxvjw/Ejjkb+DdFXx6x0B7xQK7gQhG8OCSp6Hme8+sHc3XL50Gd75pF91tghh865aCz3y5MhKWcKChwcPHiC781QyPtmEF154QRlfw4YNSzCfeYz0TNIILFSoULqMP4bDa9eujc6dO6Nhw4bqOGOBTuKwdUpwnzCsTUHzlGDVOD2c9Jgz/SKxVh+1IOn5pGedRiwrsCkhJKSPczcu4s3vnkP/b2vg4+WdsSxiLc65huOhWxSqhTqhVZAvvir+MpaNOIoZr29B9Xp9so3xGBoaas6tpYYqb5JS5L+fgO2fPinnJQg5kFmff4aOlX0xf+pL8M5b2NbDEZLBwZjOJLKWLVsqY5H6gDQkaEwx5Ll9+3YVzkzt4p1duHHjhvoemNsU3wunGSGXL19WuoiJJWGyA5TRGTNmDG7dupVAqsfW0AidNm1aAukfeyG7HxMad/zv4evlY3E9+iTOuEch/FESvKvBgOIRTqjlWgBdawxB5Ro9s3WIlkbj9OnT1f6mVBQ988kSfBv4ui4QHQJ0/x6o3seaQxUE+ySWhXTOVr9+C5kYwmaImiFAFhq888475nw5Ckxr+WtC9oQeO0r48MLI3C57Mh4JPZY0IO3NeMzuhEdFY9m2n7Dr7AKc9AxEsFmrUYd8MQaUCfdEvfxNMKTfh9C7pGBIZROY5qOFrFnVn6LxSDa8azIeC9cGqkqurSAoMsl4zMlcvHhROVr4yHMT087WrVuHYsWKKaWUTDcgq1WrlqAKW2PGjBkJZFqE7Menn36qKuOY92qPrQ2ZcyZYh6jIKHy9/B1cDTyO46438cDRAVV07gjW6+AXa0C5cDdU9a2NV/p/AUe3jCkCZEUY0Hn99ddV8RXTe5iSkSJXdwMn2N3JAej0OSCyJYIgZAKMErPDW+PGjbFjh6nKnQbk0aNHMXfu3Ce6zGWKAUmXL/PMNLcvu4+wiIY5YyNGjEj3AISsAwtKNE1IIedhiDPgpz+n4dDNtbjsFoxbzjrUifVUxqNnnAF+ehf0jSuH13p+Di/vnCkCTFkziu/TO887/BQlrgxxwLq3TM9rDwYKSa91QRAyBxamfvjhhyo3m4Wq8dMS2eQjI6TbgKSwNA3FgQMH4s6dO0oGha5PyqLw9fvvv5+hgQiCYJ8etWXrv8P2c4twzS0AV1wcAFVYr4OLwQg3GPFWnu54/pkxcHX3Q06GslBvvPGGek6VgpTksBRHFgB3jgEu3sAz71pnkIIg5EiOHz+unH2JoRcyowXQ6TYgT5w4oTqSkKVLl6JKlSpKwmfDhg0YOXKkGJCCkA04dWY3Vu2ZjSMRp3HaPRZGXwezVmP5CD2KG4tiROcPUbqYSdJLMKXxUHyeyfmppnhEBgObH7XhbP4m4GlqAyoIgpAZUEmENQws7IvP4cOHldqIVQxI5va4uLiYZXyeffZZ9Zyt3zg4QRCyJkePb8P8bdMRavDHXp9wxOocAA+gdKROGY7F4gqif8sJqF2xpa2Hanfw3PfJJ5+o55999lmC7k1JsvMzIOw+kKs0UO8l6wxSEIQcS58+ffDWW2+pWgGtExydf+PGjcOgQYOsY0AyXE2Bz06dOmHjxo2YOnWqmk9JF/YvFgQh63Dj2gnM+esdXMVV5WmM8jZ5Gf0i3eChj0E9r9ro3uJlVCle29ZDtWu0TlzUMGXrwhTxvwTs/c70vN00qTYVBCHToUIJdawZIWG3ONat8JFpidQ+tooBybtsCiwzXEPdx+rVq6v5bNOlhbYFQbBfwoPu4dtlb+FU5DGccY9EiKdW+euAgjEGlI7wxbN1X0eHBo97hAspw8prpvOw81KqveE3vAfERQOlWwLl2llriIIg5GCcnZ3x448/qjRD5kOy2UHNmjVTb3JgSQOyRYsWKuEyODhYtbXTYGGNu7s7sjs5vZWhkDWJDvHHxm3fYNeNv7HPOQh5Ir1xyovHsA65Yw0oG+GJGoXaYuSzH0DvmO7TgvComjFVLm0HzqwBHPRAu4+zTecdQRCyBvRAcrIEGbpSUO8xNjZWtXgj5cuXR4kSJZAToAuYk6ZkL6Stt2lgYKASnxeshzEmAstWf4TtNzbhsnsICoV7YZ9vGH/BiHaJRO0QF5T3rY/RvT6Few7SarQkLCrMmzcv8ufPn7Y/2P+D6bHuC0C+Cpk6NkEQBI3nnntORYmZB5lY3/m///7LkI5yulVrw8LCVA/kggULKkFpTux7zLaGzAESsrahx/AbJ7q7KUMyZcoUdbOQFtjGkn975MgRZFXYdo4/NN4Q8bOw81JiKMLapUsXddxzHbsyjONisXPTdxjzZTM893M9TA1dhR2+YbjurEOISyiqxebF6GKjsGHIXsx79QDeHvCNGI8ZhPlD7MrF3wlVKNLE878AHT4Fmk/I7OEJgiAkuG6xuUFiKC7OZVbxQFKEkormf/31l1I0J/REMgdo7Nix+O67R8nhQpakffv2qtVRVFSUag1IbytbA9pj55nMgDdBpUqVUoUQ7Ped3E0Uc395I9WjRw/YHKMRDy/swp87Z2Fz5AVVDBOrZHdUs3uUj9ShqLE4hnechkrFq9p6tNmG+/fvq77mjMjUrp3GIiMWzNSXqmtBEKwLcx6Taj/M6ztTEjNCuj2Qy5cvV21vaLV6e3uriVYtkzMz0gonJxEeE57sFBUXleZ1I2Mj07RuRqBEU4ECBVC8eHGMGjUKrVu3VgVSNJq4rxPvY3rfKFkSEhJi1pdiYi49c8yXjQ/lTei5ZrU+DVNKQmkEBAQoKQHm1TKXlsfX+fPnzcvnzZundKzY5aNixYrw9PRUxm56pKMmTpyI+vXrPzGfxiA9raRu3bqqQIySB5pcVWI4Nir6s5jMloTdPoWv5vTGm99WRZedL2GW/iyidEbEOjigRJQRzUPy4YsaM7Bs5DF8MeovMR4tDH8ne/fuVd24UlWguHMciHt8vAuCIFiTqlWrYsmSJU/MX7x4sarItooHkh6apPJ9qGYuIeyUqb/wSeNFo2nhpvi29bfm1y2WtkBEbESS69bJXwe/tP/F/Lr98vYIiAp4Yr3jg5/sWZ5e3Nzc8PDhQ2Uk0qiid/L55583L9deszUSL6TMsaA+KOWe4t/tbN26VRmPfLxw4QJ69+6NGjVq4MUXXzSHz2kw0liloco8Dd6YnDp1St0hER5fNEJ///136HQ6DBgwQGlYsQtSWmC48eOPP1aN5EuXLm0OWR87dkzdGGUJwv2x5M9J2H5/L857hOGuqwMKh3siRB8Gv1g9CujK4tlCLTG4zShbjzRHwBulcuXKpbxSRAAwrzPgkQcYuBLwLWat4QmCICjee+89FTHj9U8r+Nu8eTMWLVqUofzHDBmQ1Dn74IMP8Ntvv6nwjdbCa/LkyWqZkH1a2PHgosfvtddeU/OGDx+ORo0aKa8fjcF79+6pMDcNRsJiAkJvDL0z8aFnkf02Ge6j6Dx1RLl9GpCa4UhRU26f0ChkkRI9nJquHj2W1CDVjL9XX33V7DlMCzRq6W1kOyf+mLT3oVcy1bZztiQuBsd3zcNvh3/HJZcAnHNjWwEu0MHdYECpWG8Mq/4BelRrC71Ob+vRZnvobeexTC86b3ZS5f5ZQOcI6JwA74x1fBAEQXgamLfP6yn1IBlJpHOoWrVq6vrdvHlz6xiQX375Jdq1a4ciRYqYNSCPHj2qwn1pTiTPoezrty/ZZYkv/Nt6bUt2XZ1DwsyD9c+th6VYs2aNCg/TWKNSPUVGJ02apJbRu0gj7Ndff1WN2efPn69C3SykSg3+HY1HDRqg1KIip0+fhqOjY4LwMo1QVvdzmQZD25rxqG2DRmx6oBfy559/VgYkjWTefTGv1x6JunkE/2z9AhsfnMYenwhEPcpr1BmNqBChQ1HHinit66conq+4rYeao2BhFbXU+FvhTU+qFGsAvH4YCLkNiIEvCIKNoOOGk6VItwFJsVx6jOi5OXPmjJrXt29fdWGmRSskj7uTu83XTY1nnnlGFUIx/MwqYxp28aEXkjqYNCAZvh46dGjqwsmPEnXjo7VSSg9JbYNGYHrgscrw+KFDh5Tn/Pr16yqcbjeE3scfqybj8N3D+M/HH7cdHaH3cQA/eZFoA0pE50ePJm+iTfX2th5pjsTf31+lURDNM58mXL1NkyAIgg2hhjUdL4mvv8WKFbOODiQ9QVrumsalS5cwcuRI8UJmcZjrmFI4l3mHb775Jr766iuVn8huRBpaziPlTdIDi2IoFbRv3z5zCJt5l2fPns1wcm9y0HNOdz1vgGhAtmnTRuXv2pTYKJza9QvmHV6Ai66BKkRd0jMOt3WOcDM4oKK+EpoUaYVhzYaKyLeNYYEVKxYZ+unVq1fqoev7Z4CKz4pguCAINoWOPyqH7N69O8F8OmHojEnvdZs4WjIviDltQvaGuYxMxB0/fjzatm2rDDINGmL0Qq9fv17NZ44sW7ulBlspde3aVd2UfP/996oghx7OwoULq/mWht5y5vHyTmzmzJkJlnEeDWPt+c2bN5WuJcP6mmFNOQQWAmlcvnxZrZMrV66038UZjQi/uANfrf8Upx1u4pR7DCL9dGbpHbc4J/QvOBivtB0JL5fs3+EpK3Dnzh1140SmTp2qCrlSZPMUU9eZxm8AbSZbZ5CCIAhJwEJVRhSZesP0r7REDlND3BlCuqFoPAtReDcTHx6cvMCysIU5Yk2bNsW2bcnncsaH4fDRo0ejc+fOynBjXiULdBKHrVOC78UQPA26lDojsWqcBTjMyezWrVuCZbdu3VIyRBoMV3Ki11L7LAcOHFDvo6HlUNIbS7mhFAm+hcMbvseGa3/jmIMTjnlp0i46FIo2oFRUfnRpOA4daz8p+CrYFlbwUwmAucBMSE+RGwcftSzUATX6W2uIgiAISUInx8GDB1URq6VwMKY3iSwZWEhTq1atDLlBsyJaK0Pm0MX3wpHIyEhlxFAXUatUz05QRoci2zS2khImtRU0QllhFl/6xy4wxCHg4S2cvnQBK49+hvUu19TskqEeuO8WggrhHqhSqBPeePYdCVHbKdeuXVOect7cbNy4Uemjpshv3YBLW4Hq/YDu0lxBELLK9Tu7UrduXRVxa9KkicW2KVcrIc3Q+0IJn+nTp+Oll16yK+OR0GNJA9IujEejEXGRwbgXfBeRDjEIizMgQmeA0YF3bUCJuIJoXKEPBtZsj0K5Ctl6tEIqfPTRR8p4pDh+q1atUl758k6T8UjZnhYJ+84KgiDYgk8++UTVL/AaSVHxxNfJNEmSZdSA1LqLJIe9iojzDmPgwIGq6oghVsq3aLqC8cfOQg7O1yoshSdh03VeSBletsfWhhkVQ7UosZHwD7yNkLgIROoMiNWbfjPGOMDRCLjq8uCPzj+gXB7RA8xK3kd6twnTM1LMHWJAZ8tU0/PagwG/5FMpBEEQrIUWNUl8A2yVIprEuWJZBRqN1G1j1xMmwbNnLTucsNpYg0ZRgwYNbDrOrAD1IDVNSCEehljEhPkjKOIhghGLKJ0DDMpwdIAORrgZ9HDUe0Hvqce7/b/PlmkN2Rl63KmLyrxX5vWmyPkNwPV9gKMb0Gy8tYYoCIKQIuwCZ2nSbECyajUrwmojToTdUfLkyaO03DQDkqXt1LNkUvyJEydsPFohy2A0wvAoRB2FWITrDDDoaDTqlPSOEUa4OLgjv18RODk6qbzYkIdhth61kIFcqblz56btHEhdtc2PvI/1XgS8EnZjEgRBsBUZ7TaTEqnoUGQ+O3bsUMYbRavpRmWrncRQuJpVtfTcsFsJey5nBFYg0U3L5FkN9lJmdWVmYKH6JMGeiI1C0MOruHbvDC6EXsdDfQxC9UY4G/TQG3Xw0vuhoF9plM5XGUXyllTGI5FjIevmDTH3kSffVE/Ap/4E7h4HnL2AJmOsNURBEIQ0sXPnTqXlTL1lStRpRbH//vsvsqQBGRYWploi0khMiiVLliiZFN79s3sI12Urxfgt7BieZoecxBOrhDXodRw0aBB++OEH87xVq1ahXLlyakqNqKgoJSCsTdS9TA4tOdVe80KFdGI0IDbsIW7dPYOL/hdwwxiKEL0BMQ6mELW7wQFerrlQPk8FFPMrBDcnlyc2oR0LdlHgI6QZtuDkzS1lqVL1Pm7/xPS84SuAey6rjE8QBCEtLF++XNlO1GqmLUWbhgQFBanCGpvK+FgCeiBXrlyZIN+SHkeWn3/99dfqNdvv0IPINmIUm04L/KLYcYRC1Syo0WAhCPs5Uw+Q4tDMcxo7dmySFwvm/k2e/KQYcHIyAKxWDgwMVOLa7NxjCdFOwcrERiEi9AFCYsMQ4aBHtO5x6ycXA+Dk4ILcXgXg7JR8NTp/XjQeecPj6+trTqcQsg48f1BxIMXf8IkVwB9DAVcfYPQxwM3XmkMUBCEd5EQZn5o1ayr5PTrS2KyD0oulSpXC4cOH0aFDB1Ujkq1kfBg6Ytg5fsUvuz+wmmjPnj1p2gYv4FRgb9myZQLjkTB0rYWvKQDNHMjkPA0cgyYYTej+TanNHvMtSXxPqZBFchtjwhEcGQIDjIjUPb6/YhW1Hjq4OnnC0dUTcYjDzSBTGCA1aDxqx4SQtXBxedKjnKz3scErYjwKgmB3sDUwFVQSw25xdHZlhKcyIFkYkJkVpQ8ePFA5i/nz508wn69Z+JIWdu3apcLg7F2r5Vcy5k8dpPReROJfSBjGTgl6K+htogeSnk3BzokMxsHt32Hd3QM45x6BQL0OXgYDYuIc4avPh8Ylu6NPtVZwcUy/9iXD1vRyC1kHCu7yPNO7d+/U911MOFCsIRD+EKj/krWGKAiCkGbowGAL3sRd2pj/SE+kVQxIhpApezNnzhzcvXsX586dU29OfUUOjG3u7AmqrnPMqUEvZVpgriYnekfTAi8+YjzYMUE3sfiPt7Ah4iSOeMYgxsUBiAN8og0oHu6H3s0/Quuqlq9eE+wXRg3eeecdREREKCMyVeFwF0+gyyyg3TTAWfqWC4JgfzCFj+2Cf/75Z+XgYo0II7ksJKb9ZhUD8sMPP8Svv/6qRKU5IA0WrVBv0ZIGJCV3aHzRUI0PX9sqHPjKK6+oScuhELIo/pexZc2bWBJwFbt9IwAvznRA6UgjSjlUxpvPf4UCuRJ6voWcAaMqEydOxO7du1XqS5oR41EQBDuFNSN0pvGGmHn5DGczqkoDkjUlVjEgf/vtN1XJzEGMHDnSPJ/V0WkNK6cVJq5T+Hvz5s3mwhp+AXz96quvWvS9hBxC2ANsWzUOC4J3Ya+bK3K5ukBnNKJyuCOq5umMN/tPkn7UORy29Hr33XdTX5GRjX8mAjX6AQWrWWNogiAI6YapgEzno/Nr/PjxKpTNwmHWcXh6eiKjpPtKyeKRMmXKPDGfhl1Gcv34IfhhNC5fvowjR44gV65cKFasmCpcGTx4MOrUqYN69eopLyelf4YOHQpbkN4QtmAnRIdhyeLRWBNyAI4w4oCnK/RGoIBbJbxa4Xn0rPesrUcoZDVOrwb2fQccXQj874x4IAVBsEsYyW3bti1Onz6tCjpTKgDOVB1IvjHFKBPzxx9/qDLx9HLgwAH1d9rf0mDkc60amkns7E/N19R7pHG5fv36JwprrAUt+FOnTmHbtm02eX8h/ZzY9QuGf98MHxv24ohnHI54xKFmWFn82Hollgz6TYxHQUFtV8pZ/PXXX2kTfs9fGajyvKnyWoxHQcj2pLepybJly1ChQgW1Pgt3//77b/MyOtzeeustNZ+d8ag3S4md+PrVloRphpcuXbLsRo3p5M8//zT6+PgYp0+fbnR3dzfOmDHDOHz4cKOzs7Nxw4YNxpzC9evXeYVRj4KdEh5gfP/L1sbmP1UyVplXRU19vqlu/G7tl7YemWCHfPrpp+o3XbZsWWNsbGza/9BgyMxhCYJgB9fvxYsXKzvn559/Np48edL44osvGn19fY13795Ncv1du3YZ9Xq9Oq+cOnXK+O677xqdnJyMx48fV8sDAwONrVu3Ni5ZssR45swZ4549e4z16tUz1q5d25gZrFu3zlijRg3jX3/9Zbx165YxKCgowZQRMiQkTg/klClTlBAlQ9C1atVSHkK6SHMKOVGINCuxfcNn+OHC7zjmYarALxhjQGPXNvhgwCxbD02wQ1hxXbJkSVWgxypFW6XICIJgn9fv9DY16d27t0q3W7NmjXlegwYNVCSVKjZJ8d9//6lUvatXr6oUPktCDW2N+E0RaALyNfMk00uGqgWaNm2KjRs3IiciOZD2jSE8AD8t6o8FhpuIc42Fo9EBdUN98E6vRSiez7I/SCH7QKORxiNP2uwVmyIXNgMnlgNN/gfkeTIfXBCErJO2El/TObHe89M0NdmzZ0+C5iOErQQ1PeqkYFtBGnPMU7Q0W7dutfg2021ADh8+XJ1gW7RogZyIyPjYL1s2fY3fz32DAx6mw7paUC40KtkNrwwZa+uhCXYMLw6ffGLqJMOcpBT7lTNgw64z1/cBbn5Au4+sN1BBECxK4mKSDz74QLUttkRTkzt37iS5fnItA9mYheefvn37KiUIS9O8eXPbG5D3799H+/btkTdvXvTp0wf9+/dXLllBsBWxESF4+8dnsdPrHipG+8HZLRgt3NtiWr9P4JKSMSAIAObPn69CWdSWHTZsWMorX/nXZDzqXYBGGdNOEwTBPmBBbOHChdPetjSTiImJQa9evVQ4+bvvvsu092H64ffff6+KaVjgw8/OznxM32HTlUyvwl61ahVu376tlMsZr6dOY+XKlTFt2jRcuXIl3QMQhKdh7tJx6PtrI6z3eYAwnQ5RzqH4rP5P+Lz3F2I8CqlCr8LHH3+sno8dOzb11qw7PzM91hoIeElvc0HIynh5eSlvnzYlZ0BmpKlJgQIF0rS+Zjwy75GpgZnhfSTLly9XIXQ3NzccOnQIUVFR5rA57beMkG4Dkvj5+WHEiBFKyoYfmm0AacUmpQ+Z3WD+I93eOTWEby/s2v8Hhn1dG1+Fr8cZN8DZYETT4Dz4fvBuPFOpga2HJ2QRli5dqnRoqTsbvzFCktw4AFzaBugcgcajrTVEQRBsTPymJhqGR01NGjZsmOTfcH789QkNxPjra8bj+fPnsWnTJuTOnTvTPgO7CLJ458cff0yQptO4cWNlUGaEp2q5wQ9PHcd9+/Yp76OttBmtieRA2pbAgOuYPH8IdnveQbgX738cUD1Mj84VXkeflqmEHwUhHrwAfPSRKYdxzJgxqXdk2PHI+1itD+ArBVmCkJNIranJoEGDVEhYi2iw7zTzDj///HN06tQJixcvVvYSO/lp9tPzzz+vjDdWajMaouVH8oaWRqslOXv2rGpfmBgfHx8EBgZaz4BkNc/ChQuVS5Qn4R49eqgvIF19YwUhHVy9fhQ/bXwbW3AV3novhOt0KBllRAPXFpj4sklWQRDSw+rVq3Hy5EkVMkq1NertY8C5dYCDDmgyxlpDFATBTqAsD2tAKFlIQ69GjRoJmppcu3YtgVROo0aNlJ3EtqgTJ05E2bJlVQU2Bb21rn48B5HEdSS0sSwd5WTonNEWCqHH599//0WpUqWsY0DSwvb391eFNLSku3TpYrPEUyH7s2TdDGy4uBh6gwP2+DBnQwc35zB0iKyESYN+gbubdAAR0g+T1RnSIdRxS1U2498vTI+Vu4t0jyDkUHijmdzN5rYkutP17NlTTUlBQy4DMtwZ5sUXX1ReUUqWUSqIHW8oNTRu3DhV02IVA5Il7vxCMkOnSBBIYNA9zFj0Es7jIk67GQEvwM1gQKUQDzQq3B4jOr4LNxfLuveFnMU///yjdN3c3d3xxhtvpLzygwvAyUfabdR+FARByGJQ7JwR41atWiE8PFyFs+n8owHJm2irGJC0YnMyIiSeeazbMQ+rj36P0+7BeOhmCgU4GI2oHKFHZZeGeGvEN3Bylspq4ekpV66ckuyhHBkrLFNkF7sXGYFy7YECpvCTIAiCvXPs2DEVMmdonV7Hd955B+PHj1ehbHYRZEFwqrnfKZCmVobMcZw3b57KFeLzlFixYgVyAtLK0DJERYVhycZPsPHWWsRF6XHcO0LN94s1oGK4F1pWGIbebUbYephCNkVr45UsQTeBL6sDhhhg2AagWH1rDk8QhEwgp1y/9Xq9kl3Mly+fynOk9KIlK73T5IFklY52kqURmeIJVxDSwJ9b5mD9yXmIjNPhoF8I4MJItTPKRzigVFwJjOn9HQrmeSzwKgiZQarnsj3fmIzH4k3EeBQEIUvh6+uLy5cvKwOSSjkMYVuSNBmQv/zyi/k5PZGCkNHcxs8XjcI543mcdjPA6E0JHpO3sbahOLrWHY0WNdvZephCNoaSYzNnzlSVkVo1ZIqUaQncPAA0lcprQRCyFs8995ySEipYsKC6WaYEEb2SScHuNJmeA0mpHoapExfRsCF5t27dsGXLlnQPQsjebN67GH/89xVOuweZcxup31gxwgFFUArf9P4NPp6Zo74vCPGZMmUK/v77b1U8w2rEVCnTGijdyhpDEwRBsChUymHaIXMeX3/9dVXDwu47liLdBiRL1ZMqIGEjcPZZFARiNBjw978/4M8zv+IWQnFN2Yc6U25jmBealu6HAZ1et/UwhRwGRX49PDyULluakZQdQRCyaBFN27ZtlewiVSco42MTA5IDid+AXFNMJ1RQp6Bm/Kbk2RWpwk6Z8LBgfPj7IJw3XMEZjziV21gx2APlDGEoHVcSY3p/i4J5sm/SsmDfVKtWTbUvTJUDPwOh94F6LwLuuawxNEEQBItSs2ZNcxHN9u3bLW63pNmApFI6Y+ickuo4wwbds2fPRnZHWhkmzf371zF9yQs44XITt9x0aBDoARdDKOpE50LH6i+hS9N+Unwl2G+1dXxiIoFt04HQu4B3IaDWwMweniAIQvYsoiEcBE/CLAXfv3+/0k/TYM9GDjC55Ewh+xIU/ACTf+uLw2638ED1ptbBK84AV7hgQaufUb5YJVsPURDQr18/pXfGNmSp3vjpnYD2HwNHFgHVeltriIIgCNmziKZ48eLq0dIWrJA1iYuNxfs/9cYB/RncemQ45o41oEpYHozsNAtVSte09RAFQXHixAksXrxYPU+16wzR6YEqz5kmQRCELMoP9lZEEz8Pks3DE8fUn332WUuMS7BjVm79FvPPf8/eHLjlrINPnAE1Q/PhrV6/oEiBhI3aBcHWaD2vn3/+eVSuXNnWwxEEQbAaLKAhNi2iie/m7N69O44fP65colojGy2/iAU1QvZk/toZ2HZ5JfZ5mYS/C4e7onGQI15u9yWqla9n6+EJwhOcOXPGXDRD7ccUYXTl925A2bZAnWGAs7t1BikIgpDJxNfztpkBSQu2ZMmS2Lx5s3pkPuTDhw8xduxYfPbZZxYfoGB7tv/3F37bNxkHPSNR2+gMvdGIxtF58Uq7mahUqoathycIKXofeZPbtWtXVK9ePeWVz/wFXN4O3DoC1BpkrSEKgiBkCpndhjrdBuSePXuUWHiePHlUg25OTZo0UfpqjLEfPnw43YMQ7JOrt8/j02VDccArAOEqz9EBUfo4fF5xKlrV727r4QlCipw+fRoLFy5Uz1k8kyKMpOx4dANcfwTgKsL2giBkbXzitaHmc0uTbgOSIWothk4j8tatWyhfvrwqsjl79iyyOzlBBzIyMhzvze2Jg65XcN/HVCBTMsqIBk4NMPHln2w9PEFIc9cZzftYq1atlFe+sBm4cwxwcgfqj7LWEAVBEKwStraLEDb7xx49elSFr+vXr49PP/1Uyfiw2ocSP9md7K4DuXLrHKw6/gsO+oQrwzEv+1RHlsDUF5bB1VVywoSswcmTJ7FkyRL1fNKkSan/wc5H3kfmPnrkzuTRCYIgWJ8HDx4oPUh6JUuUKIHcuXNb14BkInpYWJj5Dr9z585o2rSpGoh2whayHpFR4fhgQW+s012Gi7sj8sQaUCnMD2/2/AXFC5a19fAEIV1MnjxZeR+Z98MmCClyZRdwbQ+gdwYavmqtIQqCIFjthnrUqFHYtWtXgvnUiPz2229RoUIF6xiQ7dq1Mz8vU6aMqnL09/eHn5+fdBrJoqze+jPWHZuDnd7hah/WNnjhtZYzUblMHVsPTRDSDRUili1blnbv444ZpseaAwDvgpk8OkEQBOvBttM0FNn85YsvvlDGIm+uKcX4448/olmzZkorl81grKYDGZ9cuaRXbFZl4pznsM35DEJ8dWgY4IkaxZrilcGPLqiCkAXRjMaePXuiatWqKa98fT9waSugcwQap0FkXBAEIQsxc+ZMVaNC76Orq2sCfUh6JVkEzXVYCJ0pBmRml4ILtuGN2a2x2fuuynUsFWlEm0oD0avtK7YeliBkGN5V8xxET/oHH3yQ+h+w5zWp0Q/wM3XbEgRByC5s3LgREyZMSGA8ari5uWH8+PGqliXTDMjMLgUXrM+3Sydih9cdJc1TP8QV0/qtQr5chWw9LEF4KipWrIhVq1aprgupdp25cQC4uBlw0ANN/metIQqCIFgNNn9JSYWC/bEz0gc7zQZkZpeCC9bl5MVDWBW0CjHOOlQJ02HOiF1wdHa29bAE4anhjS7bqaapparmfazeF8hVMtPHJgiCYG1CQkJU9Dg5KMsYGhqaoW1bJAdSyDrExcbik/Uv4JanTlVaj24yU4xHIVsQERGhQjJp4uZB4MJGk/ex2djMHpogCIJNjcikQtgkODjY3JI6UwzImjVrprnC+tChQxkaiGAdZv3xOsINLnA3xKCjUzM0qNba1kMShKdm37596NSpk8r1GTduXOp/kKcc0HoSEPYAyJX99WsFQciZGI1GlCtXLsXlGVXQSZMB2a1bN/PzyMhIpRtUqVIlNGzYUM3bu3ev0hl6+eWXMzQIwTr8e3gNFkfsQKS3AzrFlMX4od/ZekiCYBHmzp2Lhw8fqvNQmnDxApqMyexhCYIg2JStW7dm2rYdjOn0XQ4fPhwFCxbE1KlTE8xnxeP169fx888/I6e0Mrx48aL6zEWKFIG9ExcTg24/18QVVwdUiXLCb0P2wMnZxdbDEgSLEBsbi0WLFqFx48apd8TiKU80awUhx6J1kssq1297Jd0GJKuwDxw4gLJlE3YnOX/+vKrmCQoKQk4gqx2AsxaOwbbQDcgb4Ynx3b5HueKpdOcQhOwIcx/XjAGaTwAqdLT1aARBsAFZ7fptr+jS+wdMUk/cDockFqkU7IuTD3fhoosOUY7RYjwK2QZeAJhWk2Z2fgHcPgqcWpWZwxIEQcj2pLsK+4033lDq5SyWqVevnjmBnaHr9957LzPGKDwlkZHhuOjO/uU6FNWXtvVwBMEiMHgyYMAAlUqycOFC1ZIrVZ6dDeQuDdQZZo0hCoIgZFvSbUCyypE5Rl9++SXmz59vFu+lPmSvXr0yY4zCU/LNH2/hvqMO7gYDRnWXNoVC9mDLli3YsWMHXFxcUs971HDPBbSZktlDEwRByPZkSAeShqIYi1mH0wH7AG+gQrgziuQXwWQhe3gf33//ffX8pZdeSj2PKeyhyXiU4hlBEASLIELiOSB8feFR+LqYU8LCJ0HIqmzatAm7d+9WedeMiqSIIQ6Y1xFwzw10/Ua6zgiCkCPo0aNHmtddsWJF5huQcXFxmDlzJpYuXYpr164pOZv4+Pv7p3sQQubx1dL/4aGjDh4GA17uJuFrIXt4HydNmqSejxw5UsmKpcixJcD9M4CrL+DmZ51BCoIg2Biq5sQ/b65cuVLNo2IOOXjwIAIDA9NlaD6VATl58mT89NNPGDt2LN5991288847uHLlCv78809zSEmwH64GX0BFnSPcDXoUzFfc1sMRhKdm48aNZu/jm2++mfLKsVHA1o9Nzykc7uZrlTEKgiDYGtamaLz11lsq9XDOnDnQ6/VmhyAbwKTUK9uiMj4LFizAjz/+qAxIR0dH9O3bVxmUNB7ZkUawH6Kjo3DC7TZOe8TimZLdbT0cQbC+9/HAz0DQNcCrIFBvhHUGKQiCYGdQKYdtXjXjkfD5//73vww3gEm3AXnnzh1UrVpVPff09DQLh3fu3Blr167N0CCEzGHVjh/h76iDV5wBzz/zuq2HIwgW8T7u2bMnbd7HcH9g23TT8+ZvAs7uVhmjIAiCPXbrOnPmzBPzOc9gMFgnhM1qx9u3b6NYsWIoXbo0NmzYgFq1auG///5TchqC/fDvieUo7eiJ/I56eLh72Xo4gmAx7yO1aFP1Pm7/BIgMBPJVAmoOss4gBUEQ7JChQ4fihRdeULq58TW8p0+frpZZxYDs3r07Nm/ejPr16+O1115TQr5z585VBTVjxozJ0CAEyxMY/ACHve4hwFGHujG1bD0cQXhqeLOaZu/j/bPA/h9Nz9tNA/QiOCEIQs7ls88+Q4ECBfD5558rJyDhTfj48eNVSqJVemEnhnmPTGhnb+wuXbogp2DvvTQ/+nkoFusPwDvOgNU9tiK3bz5bD0kQMgxPU40aNVLnG96ofvHFFyn/wfzngQsbgfIdgb6LrDVMQRCyAPZ+/c5sgoOD1WNGi2c0nvq2vEGDBmoiBw4cMJeH2ws8QAYOHIh79+6poh+2W+zZs6d5eYkSJdSXqNPp4Ofnh61btyI7cDXslBIPLxvuIsajkOWh0sPp06fT5n08v9FkPOqcgLYfWmuIgiAIdp8HuW3bNhXG7tevn5p369YtZQOxpiXTDcjQ0FBVuePm5maed+TIEWWY/f3336os3J6g0Thr1izUqFFDFQDVrl0bHTt2hIeHh3kdelAz8uXZM0GOEerRz5jL1kMRhKemZMmSyojcv3+/CsMkS1wM8M9E0/P6L5n6XguCIORwrl69ivbt26t0w6ioKLRp0wZeXl745JNP1GvK+2RaFTY9eQ0bNlQilJxY+h0eHo5BgwapfEgaZDTE7A3G+Gk8El548uTJkyPEzu85mwz5wr4VbD0UQbAIvr6+aNu2beqyPQ/OmbrONBtvraEJgiDYNaNHj1YR4oCAgAQOQK2uJSOk2YBkomVkZCS+/PJLNGnSRD02b95cuT7pDl28eLEyJNPLjh07VO5koUKF4ODgoATJE/PNN9+oUDPDV3wPeiEyAlXX6SFl7oMG35Ofo27dukrjMjtw4OQOPHA07do29fvbejiC8FQw7zHNqdpBN/irBlq+K6LhgiAIj9i5c6dq/uLs7Iz40La6efMmMoJjegw99kpkviPVzOnN69+/P9544w08DWFhYahevTqGDRuWZDudJUuWKG8n3as0HhmObteuHc6ePYt8+Uy5ffQwMrafVNUmDVNCryO9pRRBj8+///6LwoULq6qk1q1bK43LatWqPbEtung5aYSEhMBe2XxgoXrMG2tA9XINbT0cQcgwPO/wBq9p06YqPzm+CG6StJ0KVO8D5ClvrSEKgiDYPdR6TCrFkAVFDGVnqgF59+5dlYdEaLi5u7ujQ4cOeFq4jZS2w2rLF1980axTREOSguVUTp8wYYI5BzMlaPh169ZNrc9KzvjQeNRC3cyNPHToUJIG5Mcff6zaOGYFohGO4qEeyB+bysVWEOwc3igy3FKpUqXUjUeN/JUze1iCIAhZCqb/0AH3ww8/mKOvrGn54IMPlO2T6Z1oWKkc/3liV6iliY6OVmFnegbjvy9fUw8uLTD0NWTIELRs2VJVYyf2fmqeRH6RW7ZsQeXKSV983n77bdV1R5tOnToFe+Ve9HVc9QxDLs/HjdQFISvCm8fLly+nfPNmiAPWvQU8OG/NoQmCIGQZqP+4a9cudTPOdERWYWvhaxbSZKoHkoZYuXLllNWqGVw1a9ZMYFQSSxaoPHjwQLlc8+fPn2A+XyfVkicp+IUxDE6vopZf+fvvv6tQNb2qTCAlfB9erJgLmRTsshO/046mo2SP3DEGqMei3uVsPRRBeGoS//6f4NCvwL45wPE/gDEnASdXaw1NEAQhS0C9y6NHjyp7iI+04diZhqmI8YtqMsWA/OWXX5AVYcFPcn0eS5Uqpb7I9MCCHk70jtojsdHR8It0Qp0Id5Qobqo+F4SsxsmTJ9VNGpUfUqVkc6BsW6B0KzEeBUHINHjtnzFjhpIErF69OmbPnm1uC5gUy5YtUxKHlCBjsxV6+uKHi1lXwrQ8RlrpfDt8+LBZNSazZA1pMHKyyPbSuuLgwYNhbSi5w7wnegrjw9cpasFlIq+88oqaNCV7e2PDviXY4xMJvTEC0yo9Y+vhCEKGYMrIX3/9pU64qQqHU+ux31KGSaw1PEEQchhpKeiND2UN+/btq+onOnfujIULF6paDNZZVKlSxZxGRycXC5MZAc1MaEs1a9YMy5cvR65cuRLYUyw2zoiGd7pyIK0Ncywp/B1fo4jeRL5Ok2ciB/LfmfXqsWCMEQXzFbf1cAQh3fAES+OR6TE84SZL7GNVBDC1JlE6jSAIgqWIX9DLPMI5c+aoYmIW9CYFpQ4p3E0JxIoVK2Lq1KmoVasWvv76a/M6rMt4//33E9R5ZBZMQ2RBMbUgGeFJvCwj2PyMyzg8q6i1SmomzPM51dIJLX5K7/z666+qldmoUaOU1a5VZdvChc2Dp0WLFrBH7kdcV495YpxsPRRByBA80RLevTPvOkmYlvJLB2D1a0B49m8MIAiC5WERLVNltCm+VN/TFvTu2bPnCcOQHsu0FgBbGtav0PtI3W064FatWpVgWZY0INk/m8U4nDSDkc9plZPevXvjs88+U6+ZG0Djcv369akn1mcSDF+zApv9JO2RUH2oevSMe9yqURCyCsxJZrEbT2gUvU2WIwuAmweBEytN7QsFQRDSCZ1BWnc9Tgw3p7eg986dO0n+DeenZ/3Mhl5GhrHpGaVNRdvqww8/zLD3MUO9sC0NPXmpfYBXX31VTULqPHRicY8DcjkVtPVQBCHdTJkyRT3y5FahQjJtOCODgM2PZH1avAV42eZmUhCErA2dQZoWNImvtJKdGTFihCrq6dmzp2rWkFEy7IGkS5fJo0l1gBFsQ3BIAG49ilxXKdHM1sMRhHRx7NgxVZVI7yMrF5Nl+6dA2H0gd1mg3kvWHKIgCNkIdmBhO2ZtSs6AzEhBb4ECBeyqALh48eIJmjE888wzqk3s9eumtDerGJDh4eFKO4jJoxTd1nIVX3vtNUyfPh3ZHXvOgVy14ydE6xzgYjCie4vhth6OIKQLhlMI74r5G0uS+2dNmo+k/XTAMXObGQiCIGSkoLdhw4YJ1icbN260WQEw60ty586dYF6ZMmWUdNClS5esY0BSXoN5SswBdHV9rLnGZFGWuWd37DkHMizaH15RzqgQ6glXV3dbD0cQ0gyrAv/44w/1PNncR6a6rJ8AGGKBch2AsplfuSgIgpCWgt5BgwYp+0hj9OjRql6DHWDY+GTSpEmq5iN+Oh61H1nXoXW2Y1SXr62ZJ0k7jt5Jq+RAMsGdhmKDBg0SVO7QG3nx4sUMDUKwDFcDTyPEJRqeEONRyHqV18yFfu6551SXqCQ5uw64uAXQOwPtPrL2EAVByMEwL/v+/fuqoJcGXo0aNRIU9DIaG78zX6NGjZT2I2+IJ06cqHIOaT9pGpBk9erVCRRl+vTpox7Zn5oG59NCvcdz586pELyfn1+K1dYZ6SKYbgOSX2BSopm0xDNaCi5YhjvRtwFXoLCb6D8KWcv7uHTpUvU82dxHaj7+8+juvuGrJvFwQRAEK5JSQe+2JKKSTMfhlBxDhgxRU2Yxc+ZMledJKHxuadJtQFKEcu3atSrnkWhG408//ZQjxL3tuZWhY3QMGkS5Im8hMSCFrAPv6Ol97NGjh2oPliTMewy4AngWAJqOtfYQBUEQshyD43UQzIxuguk2IKdNm4YOHTqomD0rsKkpxOds27N9+3Zkd+y1leHV2+ex3ysaBgcH9Clay9bDEYQ0ERMTo25CGfrRJHyeIOwhsONz0/NW7wEunlYdoyAIQlYkODg4zeuyCj3TDUj2bWSSJyuumau0YcMG1Z6H6urJ5i4Jmc5fO35SxqNPnAHNa3ax9XAEIU04OTmp4pkrV66gRIkSSa+0fToQFQTkrwpU72vtIQqCIGRJfH19U00tZPSH62SkF3aGhMRLly6tqpEE++HKg2OAJ1AgWgdHZ5E2EbIWyRqPhjiTdA9p9yGge6xjJgiCICTP1q1bkZlkyIBktfUvv/yitIOYmMmimnXr1qFYsWKqGluwPsEGUwWVb+xjaSVBsFd41ztjxgxVdcjzRrLQYBy0Cri+DyjWwJpDFARByNI0b948U7efbh1I5jkyVL1v3z7VmDs01NR7mdqQLD3P7tirkHiQY4R69IKvrYciCKmyadMmvPXWW+pcEhISkvLKDMGI8SgIgvDUsBkMdSnZ+Sv+ZBUDcsKECapjBBXVqc6u0bJlS9UWJ7tjr0Lid51N+QtFfJPpHywIdgR1ydhKixpomsxEAgwGYNdXQESALYYnCIKQrbh//z46d+6szreMFNesWTPBZBUD8vjx4+jevfsT8xnGfvDgQYYGITwdF64dh7/elCjbul4/Ww9HEFKFJ6wtW7aoMHaSHF8KbHwP+L45EBdr7eEJgiBkK9544w0EBgaq6LGbm5sSQWdXHQqcU9DcKjmQrOq5ffs2SpYsmWA++ykWLlw4Q4MQno6b9y/BYATyRzqiernsr8UpZK8q7CTxKgDkrQhU6wXoM5SqLQiCIDyCN+yrVq1SWt6UTWP7wjZt2ij5no8//hidOnVCpnsgmfTO3CW28mHpNxuK79q1C+PGjVO9IAXrc/XuSTjoHOCuF0+NYN/8/vvvGDt2bOpts0q1AEbtMnWdEQRBEJ4KdgvUugiyrSFD2oR56IcOHcrQNnUZERKvUKGCEtFmAQ0LSpo1a6b6PrLno2B9bgdeUo9+RqnAFuw7eZs51F988YUKnaQKK7AdRZJKEAThaSlfvjzOnjVJorHj1/fff4+bN29izpw5KFiwYIa26Zhe6Q16Hr/66ivVfoz5kDQimc/EOHpOwB5bGd66fx51HHXwjpEOHYL9wq5Vt27dUqGTUaNGJb3SjhmA3hmo9xLgJDdEgiAIlmD06NEq/ZBQMad9+/ZYsGCBKoaeN29ehrbpYKRVmEYYrnZ1dcXJkydzjMGYHForw+vXr6NIkSI2Hcugb2risGcsmgf54evXd9h0LIKQFCywYwMCttZiGHvAgAFPrhR4HZhdG4iLAvovB8q2tsVQBUHI5tjT9dvWcj7U4aUqRqaHsJl4ScPx4cOHGXozIXMIe5T76OnoZ+uhCEKSUPqLxmONGjXQr18ySgFbp5mMx+JNgDKtrD1EQRCEHIO7u7tqQ51R45Gku7yRPbDHjx+P7777DlWqVMnwGwuWI8DJoO4F8nkXt/VQBOEJ2LHq22+/Vc8//fRTdSP6BHdOAEcXmZ63mWISDxcEQRAsAoPNf/zxh2pveO/ePRVRjs+KFSsy34BkpTVdn0zCZOycekLxSbW6UrAogcEP8PCRBmSVUo1tPRxBeAIW18XExKBt27ZKNiJJNk3iKQ6o1A0oUtvaQxQEQcj2OpDff/+9auCQP39+paLztKTbgGTva8F+2H5wFQwODnAyGtGsVhdbD0cQEnDgwAEsWrRInaw++eSTpFe6tB24sBHQOQKt3rf2EAVBELI9v//+u/IyduzY0WLbTJMB+b///Q9Tp06Fh4eHEhCnZI+jo4j72gNnrv2nHvPEGuHq6m7r4QhCgpAJ011I//79Vf7jEzCMsuGR/FedYUDu0lYepSAIQvbHx8cHpUqVsug201REM3v2bCXXQ+j+lDC1/fAwzFSW7xObbklPQchU1q1bp3rGu7i4qCKaJDm+DLhzDHDxBpq/Ze0hCoIg5AgmTZqEyZMnIyIiwmLbTJMbsUSJEkr7kTlM9Crs2bNHKZknBUXFszP2pgPp4u4Mo8EIP4OvrYciCGbi4uLw5ptvquevv/660n58gphIYMtU0/MmbwAeGa8GFARBEJKnV69eKp2I3Who0yVuI5uRbjRpMiBnzJiBkSNHqn6JzGXq3r17kutxGS8c2ZlXXnlFTZqOlK0JiLkPB2cH5HbLbeuhCIKZ5cuXK71Y3mi+/fbbSa+0/3sg6DrgVQion4ywuCAIgvDUDB48GAcPHlQavFYtounWrZuaGMZm4222w9F6Kgq2JcAQoh7zuBWy9VAEwcxzzz2n2hXyhjLJaEW4P7Djc9Pzlu8CzpK/KwiCkFmsXbsW//zzD5o0aWKxbaarEsbT01NpCLGQRopo7APnmDjUjtXB0TmhnJIg2BK9Xq8kv5Jl73dAVBCQvwpQvY81hyYIgpDjKFq0qHIAWpJ0V140b97cbDx26tTJ3FtRsD6x0dE46RGDg54G5PEWD6Rge4KCgpRObKo0HQu0ngy0+wjQ6a0xNEEQhBzL559/rvLSr1y5YrFtPlXp7o4dOyxa0SOkj8NndiLiUVePJrW62no4gqBEw9nu9K+//kp5RSdXU+FMqRbWGpogCEKOZcCAASqCXLp0aXh5eSFXrlwJpowgcegszH+nN6pHv1gDihcsa+vhCDkcKhMwx+bWrVtKMzZJQu8BbrkAvZx6BEEQrEVmNIF5qrM4pTkSl4IL1uNWwAXAjQak9A0WbA9bmx4/flx5H1u2bPnkCkYjsHQQEBEI9PgeKFjdFsMUBEHIUcTExGD79u147733VA2LXYSwT5w4YRdSNtnRk3Pnzp1U1wuOeageveLEmyPYBxQNf/7555NeGHgVuH8GCLgCeOS19tAEQRByJE5OTkpazdJkyIDcuXOniqezpeHNmzfNfRb//fdfS48vx7F48WL4+vripZdeSnXdCISpR7c4ZyuMTBBS7joTGxub8kp+JYDXDwO95wNS9CUIgmA1KMX4559/WnSb6XZd0YodOHCg6m1L5fKoqChz9eW0adPw999/W3SAOQ32qmRhEo10g8EA3aMimaSIcjCoRzeIhp5gO1hM17FjR1StWhUHDhxQoexkcfMDyra25vAEQRByPGXLlsWUKVOwa9cu1K5d+4k8dXYMy3QDkj1t58yZozTe6C3TaNy4cfL9brMRmd3KsFatWkpvMyAgQOWTVa+efJ5YqHO0amNYvnD9TBmLIKQGW5u+9ZaphzUFapM0HiODgCu7gPId2K7K+oMUBEHI4cydO1dFN9mNhlN82JXGKgYku9Ak1e/ax8cHgYGByO5kditDamzSGGc1K5NeUzIgH+rj4KDToYwUIwg2YvXq1di7dy/c3d1VgnaS/DsL+PcLoOYAoOs31h6iIAhCjufy5cu2z4EsUKAALly48MR85j8y/Co8PRRrJzQgkyMo1B/+jqbdV75YLauNTRA02KZw4sSJ6vkbb7yBggULPrlS0E1g77em5+U7WXmEgiAIQlKRI05WNyBffPFFjB49Gvv27VNuT2q+LViwAOPGjcOoUaOeekDCYwOSuWXJ7eR/di1CpQgj6gU7olgB0YAUrM/8+fNx6tQp1et6/PjxSa+0bRoQGwkUb2wKYQuCIAg24bffflO56m5ubmqqVq2aKoDOKOkOYU+YMEEVd7Rq1Uq1LGM4m9IdNCBfe+21DA9EeEydOnXUzn3w4IG6QFeuXPmJdc7dPIBTbg4ooo+GTi+t4ATrwhzgSZMmmc8JzK15gntngCMLTc/ZtlDyHwVBEGzCF198odKMXn31VZUmp0WOR44cqWyNMWPGZL4BSa/jO++8ozwODGWHhoaiUqVKqvBDsAwsRKBE0ubNm1UYOykDMiD8NuAFeMeK8ShYn59++kn1VGVKC09ISbJ5CmA0ABW7AEXrWnuIgiAIwiNmz56N7777ThVAazz77LPKvqAzICMGZIaFxK9du4br168rdyiNR0vE04W050GGGYPVo3ucdAISrAsjD1OnTlXPeUfLAponuLYXOLsWcNABLd+3/iAFQRAEM7dv31aOqcRwHpdlhHQbkA8fPlTh63LlyintN+2NX3jhBYwdOzZDgxBSNiCTMs4jHSLVo5vR1epjE3I2X3/9teqUVKJECQwfPvzJFXi8bvzA9LzmQCBvOauPURAEQXhMmTJlsHTpUiRmyZIlSiMyI6Q7hE03J9vi0ANZsWJF8/zevXvjf//7Hz7//PMMDURISL169VRu6d27d3Hu3DmUL18+wfIwR1PXD3edj41GKORE2DBg+vTp6jnDHknqPp5dB1zfCzi6AS0mWH+QgiAIQgImT56s7DQW52o5kBQVZ6pcUoZlphiQGzZsUBqFRYoUSTCfFuzVq1czNAjhSVxdXdGgQQPlgeSU2IAMcoxTDuTcngn3gz0RHWvA9YBwXH0YhssPwnEjIBwxcabuORpOeh0qF/LB87Uff47r/uHI7+0K50cyRYJ9JWJT5J43j2xn+gSGOGDzZNPzBqOkZaEgCIId8Nxzzyn1nJkzZ5pbGvI8vn//ftSsWdM6BmRYWFiSOU/+/v7KYyZYNoy9e/duJZUUn+ioSMTBVNFarnBt2BMGgxEbT9/FTzsv4eDVABjSkBrbtlJ+swHJcH27WTsQGROHreNaoHhuU7ulwPBoeLg4KoNTsB0UuqdCAHMg9UlV/x9dBNw/A7j6Ao1H22KIgiAIQhKwhSHl1yxFug3Ipk2bKi0hLYmeVdmU9fn000/xzDPPWGxggildgG3iEhvsV+6cwT0nByDOgLYN+8CeeG/VCSzYd8382t1Zr4zAknncUTSXO9ycEhodUbEGlM33uII/IDzG/Lygj5v5+Sfrz2LNsVtoXi4vWlXMh3aVC8DdOd2Hr/CUsGhmxIgRyJcv35MLY6OArR+bnjcdC7glIe0jCIIgZAvSfQWmocgimgMHDigtuDfffBMnT55UHkjG0wXLkaS2HjUgrx1Rj3mNgJeHfeVAdq1RGH8dvYUBDYqjb71iKOLnpm4y0kouD2ecnNwOD8OiE4Swj90IREhkLNYcu62mAt5n8XbHCni2eqF0bV94evLnz5/0gugwoFQL4PJ2oN6L1h6WIAiCkAidTpfqNZLLY2NNdRWZakBWqVJFFXWwEtPLy0vpQPbo0UP1h06ylZlgEejl5YFArj84qx5zxdneAzdv12XExBnxYjNTG8u6Jfywd2Krp/IO8mDO45kwHeKvV5vgyI1AbDl9DysP38TNwAiMXnwE8/dexQddKqNKYfsypLMb/L3XqFEDTZo0SX4l91xAt2+A6HDA6bH3WBAEQbANK1euTHbZnj178NVXXyn7IiM4GNMh4BgTE4P27dtjzpw5GS77tjbUqhw4cCDu3bun8rcYguvZs2eCBuPDhg1T1c7M6dq7dy88PEx5dylx48YNFC1aVG0/cUGRJfn777/VmKtXr46ff/5ZzRvzVSvccr6LPJHe+OaN3bAVey4+RN8f96oGI6tfaYKqRaxjxDE/cu6/l/H1lguIiIlT709v57udKkpYOxOgYDhlu/j7P3LkiDoWBUEQsirWun7bK2fPnlUdxP766y/0798fU6ZMQfHixdO9nXRVJFC+59ixY8hK0GicNWuWagnICvI33nhDFQJpDBkyRH15XM5qZ3srBOL4Dx06hC1btpjnhTiEqDaGMXpWYtuOBqVyYXDD4nizXQVUKexttfd1ddLjlWfKYMu45iqEzVughfuuoce3u1XFt2B5RYChQ4cq3dckjceIQGDlSFPrQkEQBMEuYUHuiy++qBrAMGRNh8Cvv/6aIeORpLukldIdc+fORVaBYXWG3gjbruXJk0flaxLmbtIoZmEQyZUrlzLY7AnqNS1atEhVY2tE6aPVoyusHya88iAM/mHR5lDzpGcrY1SL0jbJQ2SRzVd9a2Lhi/VVyPvMnRB0mf0vtp69Z/WxZGf4u/n+++/V3WqS7J5tqr5e/oJJRFwQBEGwK/1eFuRSTJx2D7UfeT5nSuLTkG4DklYr+ynWqVMHL730khIPjz+lF4padunSBYUKmYohNH2i+HzzzTeq6wU9IfXr11e6RRnh4MGDiIuLU65rcv78edWGke9fq1YtTJs2DfYGw+l9+vRR3w/h+INjTZXK7nrr5v1tOnUXXb7+F68vOoy4R/o89lDA0qh0Hqx5rQlqFvNFcGQshs37TxXdCJZFy8F9gqrPAxU6Ay3e5gFh7WEJgiAIKRQ+lypVCmvWrDE7ozSn2dOSbgPyxIkTythiAQ2LaQ4fPmye6A5NLwwnMyxGIzEp2GaHhukHH3ygQrlct127diqnUYMeRlrSiaf4+on0OrKJ+A8//JDAGN65cye+/fZblUy6ceNGNSVFVFQUgoODzVNISAisDdNVabRvm3UJsaGxyOVpHZFmGoufbziL4b8dUJXQzDsMjUx/xVZmUsDHFYtHNMCABsXQvWZhVJWimqfm+PHj6NatW+ppK/kqAn0WABU7W2togiAIVie9zqxly5ahQoUKan2GjVnTkPia/v7776tIKfV1W7durRxbloS5jpGRkcr7yHA1i56TmjJCmuO1ly5dQsmSJbF161ZYkg4dOqgppc4XjNkzB4uwgGft2rWqoIRfDEnNcKXxxwsh14/fTLxw4cLKk6p5JJnjxW21adPmiW18/PHHqhWQLWH7SFZUBftHIurjy7g3vLEykjXvpKXhwX38ZhBm/HMWO88/UPOGNCqBiR0r2mWXGBdHPT7sVlUZvJpn9F5IJDafvoc+dYvahbc0K8HirVWrVql2hUm2umK4Wr5TIYfA8yEjHF4ujtDpTMf90euBKtpRLr8X6pfKreaFRsVi7s7L6qfBxgc8VzrrHdRzN2e9yuGmHi6f85FauSz+o4SZPZ5XhYTOLNogNB5nzZqlnFksSElKF5eevr59+yrboXPnzli4cKGyQ+gI00LH9A6yCpqGHe0rnnO5TdZk0Oi0BHScZda1L81V2KxQvn37tvmLYk9FfvBkNeEyMhgHB2Ug8Usm1JmkiPYff/xhnkcGDx6MwMBAdXFLDX68fv36qVaA7N0bH3og69atqwpUfHx80LVrV+Xh485OygjlpHHz5k1UqlRJHTysSNd2EMfMalXmUsYvyNEKd3iXoYUCuR7X53cb/2B5GHAPkTFhcHZxggOMMBhjERMTi9g4B2zavh4jR45BbOBjD2CBggVQu25tVK1UVX0Ohr31znoVni9WtBhaNGuh1gsICcDhY4eVQVC6TGn13fI74OfimDg2GIHTd4Kx9tA1bDx9G7dDDHCAB1yddPioa2U0Ke0Co4MRub1yw0nvpLYbHBqM0OhQ9Rm0HFKG2iMiItT3El8InXdCfE+OIa9XXvO6/iH+iDXGIrd3bjjrTf2VQ8NDERoVCidnJ5WrSig3EB4erp7Hr5bnZ+B2uV5uz9zQOejUdof9vBVbz97G4KZVMLWrqV1TeEQ4giKD1LpaL2ceJ9wuHzlejpvPtf3J9XJ55IJep1fz+V1Gx0XDx9MHbo8kayKjIhEUHmTe98ZH/yh1ldS+5+Tk6ARfD1846kzfW2BwIKLiouDl4QV3Z9P3xjEEhAVAp9fBxfXxMRURHqH2F793rStM/GOKY3PSPd5HHK+7m7t5u/y+AkIDzPueY1XbjYjAwQMH0b5Ve7XswJED6nhR29Xp1djUeNeMQVh4JCIbvQrXvCXg4eKRYH8m3vfcLvcfvxt+R9xHRoMRUZGm48/Tw1PN0/ZRWFSY+t593E3eZP4tt8vvP/G+53ZcnV3N+5NjCIsIg8FogLent/mY4vEXEhmi9j2/d23f+wf7q3W5XY6BE/en+i07OZqPP+04IV6eXnB1dH0coYgIVp+L29XOB9yfPK61Yyq5cwS3q50jPNw9zPuI6wWGBap9z+3y+9e2G2OISfV8wn3K8fK74773cvVKsO8ddA7w9Xx8/AWFBiE6NjrBMcV1+V2YzxGJ9qeriyu83bzN33tgaCDYKIvHn7bd4LBgRMVEwdnFWW2X41LniEjTOYKf2YH/HBzU8cBl8fc9XweFBal95OP1+LgOCw9DeHS4Wjf+OYJjS+kcwe1qTTA43lhDLHy9fNVxEhtnwKkb93H4+m1cfBCFywFRuBUQiduBEQgNC8ea1xqjdKHc6u9nbTqH77eeQ+/ahTCxcxV4u3vjQWg06n64CcZY7s84ODi5wwGPjp84HquRcNDp4eBomqfGHB2J7wbWRvOKhdX3vOLwTczecArNy+bG250qq+3ymBw1/xB0ceFwcjLAy90Tnq6mc7Mj4hAXFwY3Fyd4uLnB2UmvjFZDTBQqFvBC4Tw+6nsPjIjG9fvBcNEZUSyPF7w9vNU+uh8ShdCwYBiMMfB0d4eXmyf0OgcYmTIVHqj+Nsl9n+g40a4l/M1p+ygkLASRMZFqXe344/7k8UdSOkfwGqMz6tTvNqlriXacpLbvE6+bkSpsGo20Fyhppr1X0aJF8dprr5mdWfGhjcTfNEPHGmxPzIgpjVD+5un8GTt2LMaNG2f6/QUFKZtq3rx5KnXN7jGmEQcHB+Pdu3fNrz09PY0XL140WhIOZ+XKlebXN2/eVPN2796dYL3x48cb69Wrl6Zt7ty5U429evXq5unYsWPm5X///bexSpUqxsqVKxvHjBmT5rFev35djY3TvXv3zPM//PBDNW/48OEJ1nd3d1fzL1++bJ43c+ZMNa9fv34J1vX0dlHzy3xUxlhlXhU1FRpSSM3zqullLPdZOaNfCz+ja1FX8xiSmwoWL2TebuHiuc3zS71TSm232GvFTPMcYXQr6aa2n7dzXqNLYdMYyr5e0vjy/IPGM7eDjRs2bFDz+L7/HPnZvN2qdSqo+UVfLmoeb8l3Sqp5zvmdzfM4eVbzVPMLv1DY/PeHDx9W8xx9HY2LTi8yz2/ZuaXpMwwoaP77stPLqnk6N12C7fo29lXz8/fKb3wY8TDhPtLDOGnHTPN2+wztp+bn7ZrX/PcVv6lo/m4q/1TZPD9P+zxqHh8vBFxQfx8dHW1ed+KGiebtjp84Xs3L1TJXgrHx/Tm//Mzy5nkcJ+dx3Ptv7398nHiZjpOhvw81z5s9e7aa513HO8F2+X1xfunJpc3z+L1yHr/n9ZfXm7dRqITp+On4eUfzvKVLl6p57uXdE2w3/nHF8XFe8TeKq9c8RtQ+enDBaJzkZ6xW3DSGehMe/x63bt2q5rkUckmwXY/KHmp+kRFFzPNKvV9KzXPK7WT89si35m206tjKdPy9WNY878SJE2qe3kufYLs+DXzU/AJ9C5jnlZtRTs1zcHYwfrjnQ/M2BgwZoObn65HPPI+/X+3zxt9u7jam3wt/D9q8St9XMq/72t+vmbfx9ttvq3n8G/8If/N8bd0KX1Uwb4PvzXl+zfwSvB/HyvnPz3v+iXMEP+PFgMfnW09fzxTPEfG3y++W85vPaG7++/nz56t53Cfxj7/CpU3HT4m3Spj/XjtHuJdJeJzwWFC/lzcrm/8+wTni8j/m+VXrV7Wrc4S2j27cuGE+RwxeMdXYc85uY7l3/jZ61XkmQ+cI7qPgiGjjm0sPmtdt+d1I48jfDxiH/bLfWKnLgAyfIyKiY43F31pj1Lmazs0l3+2sXnPK1WZkhs4R2j6q/P56o1Nu0/Fe7H+NzdvN03VCiueIEuMeHydPnCOMRmOPb3cZPUuYfuPFX65mbDx9s7HJJ5uNVUd8kaZzxJIzS4z79+83/X3x4sb4dO3aVc3/4YcfnjhH5MmTJ8G6vMZyPn9P8a8Np06dMgYFBZmnyMhIY1JERUUZ9Xp9AvuEDBo0yPjss88ak6Jo0aLm99N4//33jdWqVVPPaT9xDDy249OsWTPj66+/bswK2FfJcSZA4eOURDJTC6EnlQPBiXf61sDBaFSJqvREajjncUbhIYXV87NjzyLmYQyKtcyFgKAmMESHI+b+VUTfOQ+dmw7Ouf0eb8xo+h507jo4uCRyaccCEZdNd24hhx/ndwYd9MdbY/OgWAEvXDseb/3gm/H+1vR36SYuBnjkxTR/3kc9vi1BfLd9Pq/HnruMFthYcmypvpcFQg4ZHa8h6tHvRQfk65pEy0KydRrdKYAzPQomT8JTY7TNd53TyanfdVB4DN5e8Ti/d/fFh4h+YFLocHwUos4IXq5OKpXm00evG5bOjQ/b1FbP3zxfEKeSETNIDZ2DA6Z2q4LhX+vAq0/x3O5oWaak0sXdc88XT5Nclg456HRxJyhSNZogcQbgRoDpWhEZEgl7gFHE+LDWInGkkjx48EB5MBNHXPPnz48zZ5KWL7tz506S63O+tlybl9w69k66Qtj8UHnzmsKOLKJhcj3j9vYcws4sNBd4ZoSwQwLuqb7C7u6e0DPM4aBDTEw0oiLC8dXKj7Diemk4xPjBCAcYomOg1+lQpkheFCxQTOXVOOuM8Im6Cjd9HNzylsWgZpXVdkPvXsKug4cQERWDukU9kN/DATGR4Th/7SaOnL2KqPBw+PsHYPeRs1i35xRiYh8b3qNHj8b06dMRvXMWdAFX4dHkJTgUMskjhZ36BzG7ZsMFUXA0RKpOJHFREYgMC4IuOgTu7Nv9iKhYozqROLu4wfX9OyqHToUbV40F/C/Ds9VY6Eo2M4cxo2OiTWEH54QhCl70kgtPubmY2idy3dCwUBUuY3iU4b/w6FjUnrQWYQyJMCzj6oKGpfOgebk8qF/UU0kDuXu4q3AR30PbR9wuQ3VaaJshYS2UqoXyOYbwqHDoHfVqzNo2GELko6ubqwoZquMkKlp9Nh4nHm6msCkJCglS4TSG9FydTMcExxASEQJHvSPcXB+HkRi+Y0hPC/ckPqY8VRjKFF7ieKN4TLm5m0PN/L5UeFSnU/M1Wj/TGrv+3YWhw4bim++/MW+Xn4/j9wm9DqcfTGkRIQPWIyJP6SfCmAw3Ev5+44fGuUwbL+fHxMao/Uxyeecyh4Q5TwtNatvVwo3xQ9jcBrcbGxerQs387h8tUKFQfu8MTWr7iKEslRbxKIypVjUaERgSiDhjnHm8/F65j7R0Cxdnl1RD2Co07uikwo3a+YCh8TiDabupnSO04yTxvg8OD1b7k+FG7ThJKYTN7cY/n6hUB+4LV3cV5tf2fUh4iAo1c7vaccIQdkxcjCkdRe+YIDSZXBiTnyF++oLa94+2q4WwGcZkaJzrmsONcQbTPnJwgJu7m/rOtf3J7cbf90p9giFWmFIS4oewI2Ii0hTGDAoJx4W7QahaLLfaLq98Dadtwm1/f3i56dG4fFE0LVcAjUrnRgEPHaJiI01h1Efnnvj7PvH+5MQxeLl5mVIzjEb1XXLfM+XD2dGUQsHvMTQyNNnrA88R2v6Mf5xwu+b0hRBTyN3T3dN8nHBdplAkvpZoaTlaqDn+OY2vvdwfb5fnnvCYKLi6uMPFyVXlkkdFx5hSB3j+cnVFnNGI2DgjIiMi4OXqiDw+nmo7zPu8/jAUDnHRKJ7Hy5zmcvZOCAJDTOkLjs4ucHPm7+vR8RcWrI4TVzd3tS84zqjICBT0cUHxfL7qc3MbDkaHTAthM9eQtRAa3CdJaUGz1oDrMa+xYcOG5vlvvvmm0o/et2/fE3/D92JuI/MgNViwy1oKNi7htijTx23H7+LXq1cv9VmZc2nvpNkDyZ1L0W3ty+UOGTly5BNdW1asWGGxwXEH1K5dW2kWaQYkDxC+fvXVV2EPxD+RaGPW8rDik1R3Gx7I2oEfHy+/J70+Ts7uCDa6YcHV5giL9YCLPhgTOjZAjaK+qFTIWxWQJKTaE9vwzF8K7TqaWg7GPwCq1Xxy7YCAAGW4L1iwQP1AvvzyS1Vpz6qyxAnDHpXaAZzij1fpVD7yMkYGAeH+QPhDOIbeAYJvm7yWj743dXG8+S/w8DwQ+5p5GypfL55ho9ABLvG8iebP4fbkoaxycbwSCpy7Ozti+9vtsGjfdSz+7xpuB0Vi6xl/NZGy+TzVd1qNU2EfVCjoBW+vhMnM3N/Mw0pMciefpNZ1cmOukkfa1nVyQi6nXE/Mp2H0xBj0Lo+++ITwopsYnqDz+ORJMI/5wDQeeQxPnjQZbo4mg0E9arbDmke/vco94FWmIUxZTY/h/qQxmBg3ryR0S50AX7cnP0dS+577M83bJUl8D7wIJk5O5/70847nqdd48mes8HD2SPO+T2q87k6Jjmltvs+T81VOr4+pOCM+vt5PfmfqR5fEpt283ZLc90l9ZuYtPrky4Onq+cRsl0TtRlPa9zQonlwZcHN2S9P+5HaTGi+Nbf5LbR+dvh2MXnMOwMlRh/0TS6h9ztPP9Oerw8fdSf3WHfUJC1g8kvgyk9pH6reRaDa3n+TvM5njxNXnyYNV++0lJqnt8veaxzlPmva9uulJYtM89/gk4U3N45308Zpwm0AejyePkfIFvABOSWJyRqWFpK6fSRWZ8ByR1nU1R5i3d+pNMKgfzWOQhl987t69q3Ryk4LzU1pfe+S8+AYkX2va1dnGgKTXL7GguCVggcGFCxcStBZkJTRFvYsVK6aqnvjerJauV6+eqnzi3ZpWlW1trB3Cjg+r9PI4n0IciqOp3woMbfz4zsbS+Pn5qep3TvT0sh0kNTtp0PMmgcnEaYKeH488piklun8P3DwIFK71eN5/P5kqfWsPBR55QyxBPi9XjG5dFq+3KoNTt4Ox9cw9bD17H4evBeD8vVA1LTt4Q63rpHdAyTwe+GVoPRT2NZ112YebN1QFvOmlyT5Vk/xM7777rnrOYjJNnSAB1/cD59YBDnrgmXesP0hBSAMhkTG4dD8M1YuajK0y+Tzh4qSHi6NO/X6L5zYZGc9USCZFQxCe0pnVsGFDtZzd7zQoE6h5MBm9pRHJdTSDkRKB9GaOGjUKWQJbJ2FqCfeJp8GDBycoIihWrJjR2dlZFc/s3bvXaGu0JFw+WpP+s2saK8+ta3z1q2ZWfd/Tp08by5cvb0qcLlFCFZJkKiH3jMaPChuNH3gbjd81NhojAjP3/YxG48PQKOM/J24bZ6w/Yxw4d5+x+uR/zMnkTGDXeHPZUTVv1sZz5nl3gyKMn284a1x24Lpx78UHxluB4ca4OIMxK8GCMpUE7+ZmvHXr1pMrGAxG4y+dTPvkz1dsMURBSJFzd4KN7648bqz03jpjw2mbjLHxfoNXH4Rlud+kYD/X78WLFxtdXFyM8+bNU8U3I0aMMPr6+hrv3Lmjlg8cONA4YcIE8/q7du0yOjo6Gj/77DN1/fzggw+MTk5OxuPHj5vXmT59utrGqlWrVHEvC4NKlixpjIiIMGYFbF5E06JFi1QTeGnh20vI2tpERMfhn5N3UL9ULpWfF+kYAwd9DNz1SbvNMwuKoVI09YUXXlB3VEmF3i2Kmx/QZhKw5SPgznFg68dAh+mZ7uFtW7mAmgiPSyZ9X/cPV9ptGjEGg/JMFvZ7HAdiG8WvNicUgKW8RvFcHiiRxx0l8nigdF5P5Qnh5O2ayd9fOuFnpQYZeeWVVxKEVMxc2gZc2QlQEqf5W9YfpCAkAWV3Np66i1/3XMHeS6ZUFK25wN3gSBR6FDkoljv1UKwgJAdlee7fv6+Ev1kPUqNGDaxfv95cBEOd5vjduqg5Te1HRnUmTpyoaiXYaS9++0DmUDKiOmLECFXXwaJfbtNSGpB2U0QjJCQjOlIZYdeFB+j/0z4U9HHF7gkt0eGnKrjprEN/Y0NMGPK4q46tOHDggFLYTyqvxyJc3AL83t0UMn1pB1Dg6Xp3WgqDwagSyikOTE7cDMKCfVdxzT8c1/0jVJhMa/eYFPm9XZT48A8D66jCJ1vDE1v37t2VdiibBmjFcmZ4mvixJXDrEFB/VKYb84KQGv5h0Vi0/xoW7L2KW0Gmql7WqLWplB+DGpZQxTDSPECw5fU7u2NzD2RWw9o5kLEGI6oV8UHZfF7qZBjwaI8VzVcBtoa9xelBZn4q8yQpYm5xSrcEKj4LnF4N/D0OGMr8O9tfFNiJgjXWGlUK++DjHo9LkWLiDMp7eeVBGK48DMPlB2G4eD8UF+6F4m5wlJpoX8Y3HtljnLlbr7Uqi1rFkijqyCS0dlpatf0TxiM5tcpkPDp7Ak3HWm1sgpAYdn/5fe9VrD56C9GPlCIYPWC3qf4NiptzlQVByFzEgEwnDO9x0u5gMpvm5fKqiRf5UxcPI/yRi7x2pZawNVTNp8ueCcbx5RUsTrtpwIVNwLU9wLGlQPXesHfomWTxDafEBEfGKEMyKCImwfx9lx8qw/LVlmXM82iAhkfHoWJB0w1EZsDtfvfdd6qtFrsiPEFcLLBlqul5w1cBz7RXTwqCpVh3/Da+3XZRtVfVYM/7wY1KoHO1ggnSTARByHzEgMwi8CJ/+Mxm9dzdYECFErYv82/ZsiX27NmjQgBaTiSr6Kk9VrmySXvSIvgWBZqNAzZPATa+B5TvALimLr1grzD/MbGHkTcIvwypp1IWahR9vOzHnZewYN81lTf5bPVC6FK9UJJG6dNCPbJkdVXPrgUeXgDccwMNX7H4e+dEuL+phUdtOiHl70m7cbrlHwz/4FCU8HFC8/J58WyNQqhU8JHmZlwMIikZJgiPZJ8oVSUpDJmLGJB23inBydFBaReSO4HX4RdrhLvBfn4UiQ1FJgwvWrQIPXv2VP09a9WqZe6V+lTQ83V4AeB/Edg2HWg/DdkJnuio58kpPgajEc6OOuWx/GLjOTWVz++FFuXzqotoneK51PKMQgMm1f3DFII+C4GYiCxtuNsLTH+5ffu2WZRaeBIK/odGxsLT1dF8/qufx4iqHQvBw9lRpZAgyh9XrjwumhGE+DAqxmLApHSZBcsgBqQd50D+vOsyvtl6AaNalMbYtuXh4uqEgFgHFIiyT+1BGiPsckCvwdKlS9VEUVdNw5PN6EuXLq0U/SnMmq67Q0cXoOOnwPzngH1zgJr9gfwW9HLaKcyrfLtjRWw4eVflfNFDefZuiJq+33EJHs56NCiVG1WL+KBKIR9ULuyt9CnT8t3SA0Y9z9atWyvDP9kcVm6rQifLf7gcCLXj6KWn0V6oUCF1cRMvCdvcGdT3wHZ95EFoFBxDo+Dh4ogiflI9LaQdXn94fWbFNH9rrH6OXx0tWA6pwrbjKq4+P+xRshQfdquCAQ2KY8pv/bHMeAz1Iz3x00t7YK+wxeVHH32EdevWISTkcV/t+PDCyQsov0Ou27RpU7OQKtulJSsTtLg/cGYNULk70HMechoBYdHYeeEBtp29hx3n7uNB6JM3Mt/2r4WOVQuacyhZwMOOEJqciQZD1hTFzZ07tzrRsitDAiKDTbKsrplQHJVDYQcvftfFixfP3LzhLCK/ExIZq3KBQ6JiUczPDT7uJm9RdGwcgiNj4evmlK3E+gXrQQ//1atXlWB3YlkcqcK2DOKBtFOiYuNw+Fqgek4PEwmKfqDalXnr7TuMWK1aNdXHkx5JNpqnsj4nVm3zB3vv3j11h3jlyhU10SujQa8lK4LHjx+PMWPGPLnxxqNNBuTFrWymCzzq45pT8PNwVrmQnCglxE46+y774+StIJy6Fay66FQs+Pj4+PvEbXy6/ix61CyML3rXMFeIv7PyOIr4VsS7s3+DqyECEXCCZ7x8M8WOGcDh34EOnwLVetni42Zb0usR4X0+ZaGi4wyqR716jDOq+UoR+ZEbgLuPBVzUKeWjhxbutQM41shYGo0xCI6IVWHq+IRFx0HrFOjsqEcez5z12xYsi3gdMx8xIO2Uo9eDEBVrQB5PF5TOayqaCA33R0VnI/TRWeOHwTAdcyQ5DRs27IkcMN4FUnyVBqcG7xi5jA3mk6RQLcDZyyQnE3zLVGCTQ6FhQPkgThqRMXFwjuex8XRxVDmTpfM97lN7KzACSw+YWjUCpr7FP360WQmfUwKlaC53FPV1xusXNyBvRABOBzjA+X4oivq5P1W+pZC+kC6NQs37Rm/c1Ydh6d5OhQLecH5kQLL6n8Ynw8LWqFhWBmOMAWHRsQiL4hSH2Hg3i4TjYFEZ+1G7yrElCFkKMSDtlH2XHqrH+iVzmb1C/s6ROOPmgOKGrH1nzvA1Q3icWP0bH3ZDocGp9Qsl/v7+ytBU/ULZE/u1g4BnPrvQg7Q3EhsGFFTmFB894jCyYQE8jHZUwufUq7wdFKEu9hfvU6/SZKgsxDtoqTuMjevcgHXbsXFMM5TNbwpzU4dv/YnbyhPau24xNY9Gwl9Hb8HX3VkJpbMTSF5PFwlBppPbgREqBzCftyvye5tCb26P9iu/S+dHEwvsmDOofgUO6j+lLRrzyDvJEDE9kfHTHxgu5n7RjhOuy/3m4qhXNwf6DHgrNe+o0kZ99Ju8HxL5SOs0YYYUl9OA9XZ1hJerk9yQCEIWRgxIOy2i2Xv5kQFZyuQhIiF63r3r4Odqap2UHaFx2adPH/PrU6dOoX379iocfvToUVV8A6/s+/mtweol8/HJO++o3NMZL7+s5lGQmUYkjckbAaZuOny8E5wXpUKicC84Cvm8HucRnbsTgl0XHqJ2PDmiO8GRmLDieIL3oj1CLzrbPhbP5Y5iuT3UI9s7Fs/tgdweObeIhHl+zAHkVMjXVYVtiZOjToWlNZFsQkOLRVJPE46m4UYD0+NRVTOh8cibCA2GvV0cGQLXPTIITUYfp4cPH6BJ3ZpYt2Un6lYtb/4bpk3Q810qr6fyeBMaojQe+Xfuzno1n+9P4XzNyEwOprUwb+3w4cOmm8YcAq8p5cqVwx9//KEKDwXB3pHbv3RCEXEaNdu2bcu09+CF4+DVAPW8fklT/iMJfHTeL5KvLHIKWrEBQ9oMgyeo+WIOZKKQmJAyERERynBk39X43yUNFBp0jcvkQW+3/zCuZXHM6lMTi0c0xJaxLXBicjsVZtToU68oZvWugXZVHvdkd9Q5oFWFfKhZzFeFwvmaBsu9kCiVz/vnkVuqX/jYZUfx3Hd7UOfDTag6aYMq8tE4fC1ACUZnJFybFfra77/sr7yAl++Hqf7pbHnJ0DKNSA0WjlQs4K1SCeLztLmMNOSpIUpDToPGO2VyNM8jPZKhUbEICI/Gw9Ao3A9h16RIdXPx0UfT0LxNB3jkflSgdeWK+vvyBbxRvagfihXKj7Zt2yrDj2FptuqsXMhbGZb0pvJ9ExuPQ4YMUYVcWYl58+bB19c3wWt+D5yYtuPn56cUJ6ZMmaKaLaT18/LmmdJnb72VdJ/5yZMnY8CAAep5iRIlMGvWLPMyvtbGwCJEvu7Vqxe2bNmSYBvaPjty5AgykwULFqB69epmKR2eux8+NDlFhOyDeCDtkOM3A1U40c/dCWUf5a5dvHYCYY+SgmuVfwY5BcoALV68WJ2Q//rrL8yePRuvv/468NcbwMmVJn3CEgnD4ELyfP/998oYZwXi8OHDn1yBxUl/DANylQJe3gc4Jq2hVpmSQYUSVmfTAJ07pK75NYt8HoRF4W5QFK4HhOPqQ05h5sfbwZHKWGG4W2PJf9ex+L/rGNO6HEa3Nt0oXfcPx5t/HEMRv/+3dx7wTdTvH3+6aaEUCpVdlAIyZS8VkL1B+SkiqIj8EEVU4A8KogxlCSiKiOiPpQLKHiKCTNkbZe9d9oaWQsf9X58nXLikSZukyWU979frmia5XL73TXL3uWeG0+N5s/NvAq702OgIh1yuegG37sELt7n00vqjV2nbqesUEx5Ag+s9RuGpqRQYHMTiLTJbsNFyB/R0+UeFh/AC4PJG3DUWxCri2gxWRCwJCYm0aNZ0mjFvMeXJYfqdWLb8LypfrhydPx/Pv81mzZpx8pxWZPk6OXPmpMOHD/NFGS7ONm3aRCNGjKCpU6fSxo0bueKELXTs2JG7Qe3fvz9djV1UTejXr5/V10Kwdu3a1ZigOH36dC7R9fnnn9OAAQNIL7C/r7/+Oo0dO5ZatWpF8fHx9Pbbb/PY5s+fr9s4BNcjAtIDQeke1fqoWh12HDB0ocmGLjRFK5E/ATfWmDFj+OSE7GyU/Kl0/zZR0k2iU+tFQNpIQkICn9TUWNOwsEfCjYFiQKcfUKKJVfFoK/juwu2NBXUqzYHbEwk9aqFoAKtblaK5ufOOCnqJb34YE6wFVtPiMTmoYmwudqVXLpqbHs8T4XaX+G/bztDSfRdp1+kbLJC1xOSI4NqdBaKyUXTOHCZi0Twr2RYQC6luAwIQ2dmw8mljYS1tVzvnKtgOluxmXwswd9VSypYtjFo1qpvuuXyPxVDBggV4we8Ucc1q5QVUVdi3b1+63zOEBax1P/30Ez+mfmZr1qxh6xk4ceIEV2LAdlDLb+LEiSax0fPmzeOKDceOHWMr13vvvWfSihPbeeutt/j5OXPmsGUQ9U7xmLPB+PPnN1jjMZbSpUvzPkIEfvjhhyzmbAFjxPzhohnCTwXVKyAqEc5jDZThUscQGxtLderU4bFgjl588UV68slHYQeO8uOPP9LgwYM5AVKb5dymTRsuBzZlyhTuToa55wt9Ig5H6NatG33xxRdZfn/BsxAB6YGgLIt5/OOpSwf5Nncqzuv+V1m/R48etHLlSlq8eDHHSO78YxrlqPG2IStbsAnE7qKEEg7ocKWlY+9soot7icJyEtXp6/LxQOTAvanl3XrFedGCLPKv2lXg+MwTV+5yzB0688BShjJGWGZuPcPrRmcPpYpFcnFm+nv1i3Msn6uAYEMy0b742zTshXJG0fbvuZtcoxNEhgVTtSeiqXaJvFS7RAwVigxi61DO8NB0lsYyA5fbPYbvOlSmFk8ZXMrL91+id2fu4sS7Wd0eCa1nv1hD1xNMY7ZPjbSvMPz69eupSpUqma4H9ymAFQxuS7hdt2/fzgXrAdzbqBMLS9Rjjz1GBw8e5NqvsNSB6OhoYwUGWM0gSCEe8f8rr7zCYhAt6lASDC5aiJmXX36ZLX7du3dnEaP9bn/55ZcsxD7++GOOLXznnXeobt26ThFTmYH9g0URosqmjk8PQdMFzLcWHPeee+45tnTawwcffMD7D+slhGxWQYcxCHUI/QYNGhiTHJctW0ZLly7l+xD5mG/chzUaxxzMffPmzbP8/oJnIQLSw0AM0s6H7bm08Y93km5Q7mCFIlO8OwM7K1f4OBAjrubIkSM0cMIc+uqrr9w9LK8BBd1HjRrF/8Mika5QO9oUrnpo8Xi2F1H2R989d4P4ubaVC6dzDyPJBy7iXWducszw3vhbLJRWH7pM++JvUe9GJY3rf7pwH8caoqsTStsAxPWdv5nESR4Qf9gmXLcpyGBOU+j2vWSOAURGNG6RJASB+lmbcvx6uM+/XX2M3/P1WkWpQhGDy/b5ioVY9EI44r20bnYUEvdGUF4rMzcsXLcQKzly5GARlC9fPmrSpAmLQ1VA4n8IuGLFihkF5/37942WMy2IB2zRwiB0IURhzYOALFWqFP/2IWBgSQdIPkFs+ujRo00EJEQLhCVAbCHcqhA/eghIgLHit4f4PwhKW8A8Y761QADCymcvEOR4X1y0OANYSCEKZ86caRSQEIdIbqxXzxBaBQsqYiAh7PF9R8crWGNxASv4FiIgPSwLGyc+LqgbHkKl8j/qDBIWEUo30gKoeKppNxF/AtYFzD2C0NFvG9YJKRZrG+PGjeOTGE60aiC+CVt/ILp9jihnIaKa75CnA1GGmEssTcsVMBbfRzH1PedusQDUsvLgJbpwK4neePpRSaPF/5ynEX8esut9kRz0WZtHFzWv1SzK1aQgLFVqFMvDi70c+KyJ3a/R1vxsUjYfb8M8UWXDR/Wcknxl3s1D5emnn+bfIUIkIAzRRADiESDuDZZICD6sA+EBEWcL2vqwcMUCWLMgymC5NBdUEC5ILNFa+7TbUN3M2IYlIHrgalVBJy21Q5ajqIlq9oRVQFRr+6TDQvv333/T5MmTHR6DM8M6YFXF5zphwgQOg8G8wSukHosh5GH5xIUqLiBQ1xehR4iDdHQfBM9EBKQDWdhY1FZIzgZuOtR8q/Z4tEnWJXehCSaKCvTvtnK4+oUb5+LFi7Rtyc9Uk3YR5YkzdKgRrFqGILbBoEGD2AVoQsI1ovUPrbn1BhCFeOdFCmoZVorNzYs5g1qV4eSdYnkfucxheSwSHc7Z0UhagyhF3UTcBgcGcmJLTGTYoyVHGCfxaOmlsXJmFUtxifagxjA6e7sAFqYbNwyVIcyBYCxTpgxf4JknzsDyBJGxYMECzjJOTk7meDxb0FrJVQGk7Vpl7zbU7VjbRuvWrTlZT6VQoUKUVSB0cbzC3NgKXMIxMTEmQhbz68j5BheN6AmNsBVngc8UovSPP/5gyzLc7dqLAsRZQ8xDNKoiHsmQEONDhw41XgwI3o8ISA+jVYWC1LRcfrqZmGzy+O3U2wYBGfIoLtIfwUkIIhInrYWLFlLNImuI8pUXAZkBOLhDROIkBLdSOtaNIrp/yzCPFR7V4PQlVCulltdqPc6LkDmVKlWymggCYRMXF2fxOVysdOrUiV3Xao1XNU4S4DFYDO0FSSrI9tWC+7Cw2xpraCkJJV0/+CwASycsrvCY2OMpQdIR5jur7mvwzTff8Hs7s1QSLNFt27ZlyyNCChAOULnyo1h0WE/NL1LVz8SkDJvg9YiA9EAQ+A+Lh5bkB4lUKg0lNaTuIQ6mEJCLNhygka8Q0aW9RAlXibLndffQPA5YH9RYUcSRpTu5XjtOtH2S4f/Gn/tdb3HBNuCK7N+/P1shEQdnDygXBcEHzEUfsnWXL1/OJXBgpYuKss3DgmxrWL8Qc4mLImT+jh8/nt2q7gDCCF4RtYwPxjN8+HDen5EjR5qsi9qQ5nUYse+qhREWPTUDG/GDsEAiHjQzEGuJMcDKe/LkSRb8kyZNYotg8eKmiWmYb3MQY5ouNjoDN3bLli05M9w8JAYWSri4v//+e6MLu2fPnhwXa2s5I8E7EAHpJVwKe0BnQwOpQqqpC80fgQUSV7iHjhylI0pFKhlwgujkOqJybd09NI8DJ7G7d+9yBi2sBulYOZgoLYWoeEOiOP+pLyrYR/ny5dnKhLI82jhBW0AWNeIk4ZrVuogBhAaaMqDzCr6n2jI+GaGOBXF2EFtwi6IOosXqAjqAOEWMAS5yuKxhlYPlFbGA5pnT2F+thRF06dKFxR6EJwSm6uZH7COSkrQWPmtgLrDAqotYz5o1a9KqVauMyS1atN2+tKWCIFjh7sbngKxva9SvX58TdCBEO3ToYPIcPgOIWQh6CH2ENWB9KePjewQoYlN2CDUGEj+6woVNM0RdwdNTytKdoEDqHdWeOj+vX1FYTwUdL3Cg/aZ9Sap5azFRlTeIWn3j7mF5FEj0KleuHB09epT++usvatSokekKZ7YSTWlMFBBI9PZGonxl3DVUvwFZqbAO4SRtLSnFU0HMG+La4GK1xyWLUwxEJLKhe/fu7dIxejuwpqLSBMrgANRShKjTy7IK4YgLTdTgtNfS7E2/Nb3P376KWCC9gDMXjrF4BOWfTF/I1x9ZsmQJX2nT4T+Jfl1ssEAKJmB+UHMPMVTpxCOuG//6xPB/pVdFPAqZgpI6uBhBZxFbEzoQQoGi2HCtdu7c2eVj9PYLPlh6UTxdBReA2uLprga1GyFevV08CvogAtIL2Lp3Bd+GpilUMa66u4fjEbB4BEWfIQoIIrp+gujmWaJczs+M92Zw5W0xcQZ9xEu1ILp11pB5LQg2gFg2e0ANQmRwo4OJiJLMj2nolKPFFV1zMgJ1NAXBVkRAelgdSEucQncQ7kKj+GUXmoy4dV+h60Fx9ETKEaLzu0VAPmTFihUc+5SuZI9KUDDRsz2JanbPcstCQbCGREgJgu8iVZjtBDUgUSgVgdB6cf2uobVXzlT5uLT8/PPPXC+t1+8P+yRfMu25668gEB8xogi8R5ePDBHxKAiCIDiAKBIvICktiaJTFApPFYOxFhSoRcmKc7dTDZaOiyIgAWLUkCGJshko4mxC4nWiKU2Jjq501/AEQRAEH0AUiRcQHB5E14MDqEKafxcRN0fti10i+CLRTy0N9SAFLgGCpBmLYRYbxhKd2Uy0cpChbI/UfRQEQRAcQASkF3An9RZ/UrlD7e+v68ug5hrKg9C9hwXEb54huneTKNy0nZo/YrUgc+3eRChGX6KRiEdBEATBYcSF7QXcVhL5Nnd4PncPxTMJz01JEYUoDW7sS/vJX0G7OPQczjBxITw3UZNhRMWsFwkWBEEQhMwQAekFBKQlU6l7CqUkSRtDS6CLQ97BJ2hnXE+iPJZ78vo6ly5d4q4XKAKMgs/pSLhmqP0oCIIgCE5ABKQXcDosjQ6FB1CuSLFAWgJFihPu3actt2KIIvOTPzJgwABuH4aWcM2bN09f8/GXNkTTWhjqZQqCA1y7do3rOp46dcql74PtIzzFvF+0r4OYZbRx3LFjh7uHIgg2IQLSwzl/6TTdUrvQxD3t7uF4JGp/3a1bt5I/snPnTpoyZQr/P27cuPRt5v6ZSYRaoihzlE3iQwXHGDZsGLVp08bYq1oVeuqSJ08eLh+1e/dum7eJvsnPP/88eRPTpk3j/s7a++ocBAUFccF0HJPQmxvtVm3dXxQS79OnD3300UcWnx8yZAi9+uqr/D8+g6+//tr4HO6rYwgPD+f77dq1o9WrV7tFnKNWcunSpXks6AuOkmuC7yEC0sPZfsBQbiVYUahKqTruHo5nC8iN64h2TCFKTSF/AfGOcF3jtmPHjunbniGpaNUQw/91+hJFSCa/YD+JiYk0efJkDhcxZ+XKlXThwgVavnw53b17l5o1a0Y3b94kfyJnzpw8B+ixvGnTJu4gA9FUsWJFOn/eUMfXFvAb3rBhA+3fnz6WGy1JW7dubfW1EKwYw+HDh/m9IXIbNmzIwl9Pvv/+e+rfvz8NHjyY9wPCF/WTf//9d13HIbgeEZAezqmLhgNJLulCYxXUOwTHTp2lq7M/ILp2lPyFWbNm0caNGykiIoJGjhyZfoW1I4gSrhDlKUFUvZs7hijYyoME+xftxRL+x2PJ9zLfrgM9klFTtGbNmumeg+Uxf/78HD4xZswYjseFNwCCBr2czYGo+vTTT1lg/PTTTyyMVOuZtkHDiRMnuJsSvtso2YUC+VrmzZtHZcuW5XHB4vbll1+aPI/Hhg8fTm+++SZFRkZSbGwst1R0BRg75qBAgQJseYPQhpCEoP7www9t3g6sl8888wz3D9dy9uxZFmNNmza1+lrsI8aA/axTpw7vK+Z54MCBLCqdAbZZsGBBSkszjceHZRrzDH755Rfq1q0bt1AtVqwYtW/fngX1F1984ZQxCJ6DCEgHTPNlypSh557TJ4v16u14vo1MDdDl/bwRHHThJgHb0soTpSaTv8Sk9erVi//HFX/hwoVNV0Bh9W0PT5jNR0vXGU9neEH7l0Maqw7+x2PTXzTd7tfl07/OTtavX09VqlTJdD24LNV4PgiKgwcP0vbt243Pw729Z88e6ty5M7tr4WaFKILlDMvTTz9tEteLdeBuLVmyJL3yyiuUkpJiDNvAayFO9u7dy2IUYgnuZC0QlRC2eN/u3bvTO++84zQxlRmIF4VFcfHixZSammrXBTHmWwu2gXMOLJ32oHonINKdwUsvvcTHnTVr1hgfu379Oi1btoz3FaD7VbZs2dJ9L7Zt28aNHwTfQQSkh7cyTEi+wbcRqVKzLyNUy8iWiIZEBZ4if+C9997jBCJYPHCiNQEZ10v7GGo+lmljKBouCA5y+vRptjxlBNzWn3/+OeXIkYNFEC5omjRpwuWlVPB/3bp12TKF9SAsYEGE5QwL4gBV8J1u0aIFi0e4QTGGY8eO8XNfffUVNWjQgEUjnkdsYY8ePWj06NEmY0JCGYRj8eLFObYwb968JuLH1ZQqVYqT2yC6bAXzjH21x31tDXSkcmbiEy7WEaIwc+ZM42Nz587leYW1GOAznzRpEot8iFckBeE+xOPVq1edMg7BM5BC4h5OcloK5UlRKCw1xN1D8fg4SLjD/CWRBu67X3/9lYP2sd/mV/y0Z7ah40xIBFGT4e4apmAPH9seK2ckSNOqslQrwzYCzOwCPbPeoenevXvpv2MPgdUQiVsJCQksDBFWkS+foWJE165d2RIJwYd1IDzGjh1rc6tSFbiGweXLl1mUwbIJt6kWuH6RWAJrH34X5ttQ3czYhiVmzJjBrleVP//8k2rXrk1ZQa3Jive2FYhqxJyq3L59m/7++2+OQXV0DPa8f2bA0ojPdcKECSz+MW+wBKvJexD1uLDFRT3eG9+FTp060ahRo9In+AlejQhIDycgPICuBQdQjQj/LE9jbyIN3CRpd65QYGQM+So4Ab799tv8P6wq1apVM10h6RbRX588SpyJMnNtC55JaPasvT4o2LA4e7tEbGG6ccPgDTEHghFhPYiF1GYng1atWrHIQIF7WBdhhUKrTVsICXl00awKIPPYO3u2oW7H2jZg4VOPI6BQoUKUVSB04XbG3NgKXMIxMTEmQhbzW6RIEbvfH5bPK1eu0BNPPEHOAp8phCHqzeLYA3e79qIAAhhVIX744QeOh4X4R+wkYjS1+yV4PyIgPZy7aXf5Nqe0McyQ8uXL84ELbrQj/Z+gUiOOE/lg3UwcuOGSgysI+4wA+XSsHUmUcJkoT3GiWj3cMUzBx6hUqRJNnz7d4nMQNnFxlgv4BwcHs/UJrmsISFiq1DhJgMfsiQ9UQdgGkse04D7c2ar10V4gcLA480IPFleU7bHH8rZv3z6eb6372tzaaivffPMNv7czSyXBEo2GBbA8IqQA8eeVK1e2KN7VuGwkBbVs2VIskD6GCEgP507Afb6NFgtkhuBghSB/lMDYei6VSl3a65MCEmVCEMOFEzNc17DumHBhD9HWHwz/NxsliTOCU0BcGxK1YIVEHJw9/Pe//2XBB8xFHzKlUf4HiS2w0lnt4W7G//3f/7H1CzGXyPZFhvb48ePZrequCzu4bXGLi1iMBxng2B/z6gioDWlehxH7rloYYdHDfgEkDcECmS7G2QKItcQYYOU9efIkC37EHo4YMYJjQLVYSiRCRru5xTYjNzYEITLD1dqUKkeOHGFPEKy5+L4gfAGiGMcrwbcQAenhBKem0pNJafQgWLLXMgMxNxCQW86lUidkIBdvSL4GTjI4aK9bt87ESsGkPCBa+A6RkmpInCnewF3DFHwMWLthZZo9e7ZJnKAtlChRguMk4ZrVuogBYumQkIhMaZS8wcWRWqg8I9SxwAIPsQU3KcoGIZnGHSBOEWOAixwua1jlYHlFFrR55jT21/y3i7I/EHsQnhCYqpsfsY9INrJk4TMHc4EFVl3EeuJ4uGrVKmNyixZYgs1BqSAIVri78TlkVGmkfv36nKADIdqhQweT52BRRvY7noMgxfujpJEtn6vgXQQoapSvYLclCCdz/OjSlU9xIrWmlKW7QYH0YZ7X6LWWttcT80eQDYgyE5XyB9Ku8Z2J/jOJ/IoHiUQrPiXav5Co+xaiHBJv5GkkJSWxdQgnaWtJKZ4KYt769u3L1iR7XJE4xUBEIvSid+/eLh2jtwNrKmpefvzxx3z//fffZ1Gnl2UVwhHuadTgtNfS7E2/Nb3O376OWCA9mNPnj7F4BOWLP+Pu4Xg8uCr+a8pwqnZ4hKF1n4+AEzDq16FsCQSyVUIjiFp8SVRvgHScEZwOSuocPXqU4uPjbU7oQAIH4t/gWkXtR8E6qJ0JS69a2xWgEHu67lIuBAXjIV69XTwK+iAWSAfR4wpmwaofaeC5bylEUWhbh13SicYWbsUTjS1DFBBkKGkS4l1WHkvAVQfLBOIecQJP5wpC4XTsrwSoezzebIF0BLh0kcGNZA5zV6cguBKxQLoesUB6MKcuHODb3CnSxtBmchYkCo8muned6MpBooJmcYJeCFxKSBooWrSo5TiiNcOIzu0gavMdUe6i7hiiIFhE7BOC4LuIgPRgriXEE2UjyiFtDG3m8JEjNPnvQIpITKLBrff6hICE5RE9hi2SeJ1o2ySiB3eILu4RASkIgiDogvi8PJiEBzf5NiJN2hjaCmKtRv95gibvTjb0gvZSkLGKWCTERWUIYh3fXkfUYBBR6VZ6DU8QBEHwc3xeQCLGAeUIUMkfba3mzJljfA5lBipWrGhcUOB24cKF5CmkkNrGUNzXtoJyIN1ebECf1wujtHO7yBu5f/8+vfDCC1y/DTX0MiW6GFFtyW4VBEEQ9MPnXdhw/6E/KgQirFMoNt28eXPKnj071+pSC7qiBhniyxo1akSeQkA24jaGNSMMfWCFzMHnOnHiD0TfVia69C9RcpJXJdIgZgy9g1HnEfXjUDYlHWjFtqQnUZnWPlnrUhAEQfB8fN4CieKuEI8AxVWREQj3oDmLFy/mMikQIJ7CnbREvo0KkzaGdgGLXPYYotQHRBdMOz54OigEjPZnuPCZN28el/VIx9oRRLt+IvrtVaK7l90xTEEQBMEC0Bfo1AMDAHrDo0g8DFSZZYy/++673JEIheP/85//cB9xLagJCgMYuo+pmob8XUDC0oLm7AULFuSSD5ZcyN999x1bB5GKj04GaJPkCDt37uQq+ZZqmKmlUjyJOwFJfBudXdoY2kOaotDO5BI07R8IyD3kLUyZMoWGDh3K///www/UsKEF6+I/vxKtG2X4v8UYohyP6TxKQRAEwRoQj+gWtmLFClqyZAlrnLfeeosyArU/f//9dw6xQ/eh8+fPc/UNc+Cd8iidoriZpUuXKgMGDFDmz5+Peg/KggULTJ7/7bfflNDQUGXKlCnK/v37la5duyq5cuVSLl26ZFynQoUKStmyZdMt8fHxxnWuXbumlClTRtm4cWO6Mdy6dUuJiYlR7t27Z/O4z549y+PFrat4eUJZ5T8TyyjfzuzjsvfwRU6dOsWfTVBQkHLz5k3FG1ixYoUSHBzM4/7kk08sr7R/oaIMyaMog3IqyopBeg9RcAI4xhw4cMCuY42ncPXqVT5Onjx50qXvg+3jd7B7927Fn7h//75StGhRZfv27e4eis//1tTzN57H+V9dkpKSsvSeBw4c4O1qP8M///xTCQgIMNEjWnCOCgkJUebMmWN87ODBg7ydzZs3p1t/0KBBrHk8AbdbIJs1a8ZWFyQNWAKN2NEvFV0MkAgzceJEioiIYGuNCuIY0V7LfIFVU01KeP7556lfv37ck9WcRYsWUePGjTMs7IttoN+puqBxvas5mU2hw9kCKTqXxEDaA+ollixZkq3NaM3l6eBqFS4LtCxDsWX09E3H5glEszsRpSUTlW1LVH+gO4Yq+DHDhg2jNm3aGGuRnjp1ir1G6gL3G46ju3fvtnmb6F2NY7M3MW3aNHZNau+rcxAUFMRdXOApw+8Yfa1t3V/0sO7Tpw999NFHFp8fMmQIvfrqq/w/PgPE9qvgvjoGJIPifrt27Wj16tUm21A/MzX23xVcuHCBj2M4BqPlZc+ePa0e89Rxa/dFT6ApoqKijAsSF7PC5s2b+buBZE4VeJIwD1u3brXqGU1OTjbxOJUqVYpiY2N5e56M2wVkRqCECSZXO7H4IHDf1olFUgJ+tGhz99prr1lcxxb3Nb5Y2i8avniu5NjpvZT4sLNIxZJ1XPpevoiaDPXXX3+RJ4M4F7SIw0XJs88+yxdGOKCaJMws+5hoeX98m4mqdjH0+JauM4KOJCYm0uTJkzmey5yVK1eyaFi+fDnHesEocPOmoQSZv4B4N8wBOpxs2rSJXZY///wzx6rBHWmP+3PDhg0ssCwZOlq3bm31tRCsGAOqi+C9IWRwroTw1xMYW2JiYuiTTz7hvt7Wvk/FihWjkSNHcm6Cuzhw4ACLfHXp3x/HWce5ePEiPfaYaVgR4tmjo6P5OWuvwcWD9qIE5MuXz+prPAWPPgtdvXqVrUiYSEcnduPGjTRr1iyOrVTL9ezd+6hPMr40iKls0qRJhtvBF0v7RcMXz5XsOryOb8PSFHoy1vKPULAOLCFgxfxfiPbOJU/k3r17bNE5ffo0FS9enBYsWMAB0kaQQT63M9GW7wz3Gw429LoOlLqgvkhicqLdS0paivH1+B+PJaUkZbpdR3ok47tZs2bNdM/B8ggRAKsLCt7jogjWFgga9HI2B8fgTz/9lAYPHkw//fQTCyPVerZ27VrjeidOnKB69eqxxwlCxNxogCSzsmXL8rhgyfryyy9Nnsdjw4cP57ixyMhItuj8+OOP5AowdswBkjZLly7NQhtCEoL6ww8/tHk7sF4+88wz3D/cvBwdRGXTpk2tvhb7iDFgP+vUqcP7inlGYh5EpTPANuHZS8OFrQYcxzDP6ryjdeXrr7/OxhZLVKtWjUaPHk3t27c3PebpDOYM4l9drI0F3kuttd3ScujQIfI3fL6MD6w65l92LfiCm2c7WQJfLCxI6MGSaYHnLHL20kG+zZWqUFBIiEvfyxdB7c+gwAA6dukundy8iJ4o/yJ5EvhO4gCLEy2uTv/44w+uEGDSYea3DkRnNhMFhhA9P4HoqXbuHLLgYmrMrGH3a8bUHUNNHjdc/K46s4r6/N2HquarSlObTjWu03ReU7px/4bJ6/Z2enQRbQvr16/nDNDMgPsU4PgIQQG36/bt21kwALi39+zZQ/Pnz2dLzcGDB9n6PnWqYbz4LagWuwEDBrAgLVGiBP//yiuv0LFjx9iiA88UXLQQofAeQax1796dxSw8TioQlZ9//jkX5Z87dy698847VLduXS7h5mqwf7AowqsAQwjc27ZQvXp1nm/zKiE4pkHk2MMHH3zA+w+Rbo+QtcZLL71E7733HocGoWqJmnW8bNkyvsjwVdBKVvu9skSxYsVYwF++bFoZA6FJmCNrllY8jt8LrPZaKyR0iTuts15vgcQJFT86c4HnzolFqj2sj9orZVdwLeEC30amevRH5LHgQFurqsFyu+JWMfI0YF2B6y8kJIQtj4gXMnL9BNHkRgbxGJaT6NV5Ih4FtwIruRpTbg2cACFWUIYEIqhw4cLs2VHFIcD/EHA42WI9CE5cmON4jgWuPBXEAyK8A78NCFGMAQJSjY2HgIGFDc/j5N6jRw+2amlBzV8IS1j4EVuIc4qecdGIZUO8/LVr12x+DeYZ+2qP+9oaEOQQsoh9dAawkCJEAaXGVCDMMa+wFvsqcMnjs8xoCQ0NpVq1avHvABc4KohDhcEAcbGWwIUZzgOrVq0yPgaL8ZkzZ3h7noxHWyDxgWByMbFq4DE+CNzHwcKXuZd8iyicKFwEpMM0avECbdj2D63YdoAyLqKgPzihbdmyhZO94G4y4cBiomvHiKKKEHWYTZTPtfG2gmewtYPlIPuMCA16JLgaxDbgbQQGmB4zlv1nmVPCLawlGSIxEbHpCQkJLAwRMqSGHSEBEpZICD6sA+ExduxYm94TncNU4BoGsO7gZA3LJdymWuD6RTKG1tqn3YbqZja3EKnMmDGDunXrZrz/559/Uu3atSkrIAZffW9bgahGjKAKLLQo7YIYVEfHYM/7ZwasqvhcJ0yYwOIf8wZXND5ff6d06dIcZoD5QcIvkmOgVTA/6gVYfHw8X/wgThUXWvCCIuShd+/eLPhh/ICVF+JRGzKCiyeERCB8D79HNREK+RjaCy+/EpCYEPWqEpw8eZInBhOJWA5MaqdOnTi+BpONAwQOVMjK9mUeKA/bGKa5Lz7EF+IgBw0axBcc9riQXAm+u2qxerjRLLrSnvmAKDWZqPJrRJGe7cIQnEdESESWXh8cGMyLs7cLYGG6ccPUDa4CwYiTGNzH5okAqPELkQErO05yOKG++KJt4SSwyqioAiijcKTMtqFux9o2YOHTWokKFSpEWQVCF4IAc2MrcHfC4qUVsphfS/WLMwOWzytXrtATTzxBzgKfKUQpwm4QmgB3u60XBf7AjBkzWDRCJEJUI9t83LhxxufxG4CFUXuRgPlT10USEiz3EOha0NYWFxIqlSpVMmomtTKC3wnIHTt2mJi+IRgBRCPKIyC+BT8ABAJDeSMAG/EW5ok1eqFXDGRAtgBuY1gre8ZuI8E6uOiIisrJJ76dc76i6u0ttAXUEfz4EUP066+/GuOHjBxZTvREXUPbRZws67p3rIKgBSer6dOnW3wOwiYuLs7ic4hXxLEcrmsISFhi1DhJgMdwceeIpQcJklpwH+5sRy8UkVCBxVnA0gmLK7xn9ljn4JVQxYHqvja3ttoKklnw3s4slQRLNIpcQyjB+IOL4MqVKztt+95OdHS0iYvfHIg91TKtnVNVW1jD1WFzXikgERhsPpnmQM17issaMZBYUK7BkStCW7mjJPCttDF0HJy8GlR5kuav3k5/zfrR7QISB3NcDP3vf/8zFZA7fyL6/X1DX+v2vxIFu8cdIQjWgEUElShwMYY4OHuA5QSCD5iLPpxMUf4HFhlY6axl7VpKaoD1CzGXMDIgQ3v8+PHprDZ6gXMYDBy4RQwcxoMMcOwPStVoQRUP8zqM2Hf1fAKLHvZLTcCABRLxoJmBWEuMARYuWKUg+CdNmsQl6BAyo8VSVjYy2s0tthm5sVu2bMmZ4WptSi3q/sHDiGMe7uNiQS1/BwOMWskE/8Oti3UQF2s+VsGDcXclc2/F1Z1o2v7wlFJuWjnlx4VWupIINvH96CH8OdUpGqIoyVnrMpBV0OVgyJAhSmJioukTJ9crytD8ivJnP0VJS3PX8AQd8OZONNWrV1cmTpzoUMeY2rVrc3cwcy5fvqw0atRIyZEjB29rzZo1Frd748YN4/Mqc+fO5e5i6OIRGxurjB492mTb6OoyduxYk8fQwQOdPLLC1KlTlaioKJP7hiKtxB1H8Bzm6rPPPuPuJlo6depkXFe7dOnShZ/ftGkTd1pTjxErV65UChcunG4M5vuG++q20LkN89GuXTtl9erVJq9T59bSgnOZ+rx2ni2RmpqqFChQgNc9fvx4uuctbR9jzGwcdevWVfTsROPKTnL+QAD+uFvEehNaF/bx48e5PheyDZ1Nu4nlKC1AoXq5WtG7L5tewQq2c+L4cYorXpxCAomu7V1FkWXq6/r+sNggLizTIPZrx4miixnc14LPkpSUxNYhxKRl1PnKE0HMW9++fdnFao9LFqcYlOJBNrQaoiRYBtZU1LxE2SHw/vvvsxVSL8sqMtThnkaVCHstzd70W1M9iK46f/sLkjblgWV8UuGCCDO0MYzJLV/urFAsLo5aVylEvWuF0v2T+raFglsGbja4n9Jdp237H9ElTbeJPHEiHgWPBiV10GEF32tbgfsSrmW4Vn098TGrwChRvnx56tWrl/ExFGJH7Uq9QC1HiFdvF4+Cn8RACuk5eGYPJT28wq9SStoYZpVF3/QlWjGQ6Hb69mCuAi3F0D4TVmpkoOKgbMzEXDeaaPVQooi8RN03E+UwbX0lCJ6Kpb7GGYEahMjgRgcTESUZgxhBtP/TAsGuJ+Z1NAUhI0RAeiB7HrYxDE9Lo7jYR3XMBAcp8rCW1tmt8Ke53NKHQvdIkjly5AiXokIhWRaPeG8Ix/VjDCtW70qU/VG5DkHwNSRCShB8F3FheyDnrhzh2yj7q1sIlihYkRJSQ2npvxcofs+jOlquADE1KEuF+m+IrUFMEdfowon0r08eicdGnxE910/c1oIgCIJXIgLSTpBAg1IEKD/kKm5IG0PnEhxGL/0eTC1m3qP+vbq77G0QG4uuHKp4hOURnTkIhYuX9iHaPN6wYrPRhmLhgiAIguCliELxwCSaxNTbfBue6v7OKb7CZ4MHU96IAKoecY4o8brTt4+6b2h7huw+tFrbtGkTZ55SchLRvDeJtk9CeXiiVuOIanhaY0VBEARBsA8RkB5IipJGeZMVCkuVgtLOomqbt+jMiBrUo2oA0Q7HespaA4V+EfOIFmRohbZhwwZDUWAI1V9eINq/gCgwhKjt/4iqdHLqewuCIAiCOxAB6YlkU+hqSADlz5H1XqzCQwICKPy5Xo9K6KTcd0qCAHqcoocumts3bdqU+25zwsyN00RTmhCd2UQUlpPo1XlET72U9f0QBEEQBA9ABKQHxkDeUQxN1nNnkwxdp1L2BVIiC9KC7eeoTaNnuUBvVsMZPvjgA94O2nktXryYsmfPTnR+N9HkRkRXjxDlLET05jKiYnWdthuCIAiC4G5EQHpgDOTdAIN1LDp7AZe9h18SFEJ3nv2U/vtXMC1eu4OmTJmSpc3B8ojabV9//TX9/PPPj/rIbptEdPcSUb5yRF1WEOUr65zxC4IbuXbtGtd1PHXqlLuHItjJsmXLqGLFipSGhD5BcBIiID2Q4LQ0KpmURin3s2YhE9KTs0YHGjRkKP8/cOBAunPnjs2vTU1Npd27dxvvw2WNll+wQpq0Kmw+iqjG20SdlxJFSRiC4BsMGzaM2rRpYyhLRcRCEt/7f/75x+S+ukRGRlLZsmX5ovvo0aMm25o2bRq3+HQVM2bM4JaAERERVKBAAXrzzTdZAGcV7DsuGLX31f0NDw/n++3ateMKDFrM58rZIHQGFSAw5/nz56ePPvrIxMOCYxUucDEvlkC7v5UrV7JhBOO8efMmP67ex4L2lVFRUVSpUiX68MMPuVmClsGDB7NIdQXJycn02WefUVxcHLclxGcLUSy4FxGQHkbKgwd0IozoSLZAypu7iLuH45O8/fbbVLx4cS743aVLl3QHQmst2SpXrswHaWRaqxQqVIjoVjzRqs8N5XpAaHaiZl8QZYty5W4Igm4kJibS5MmT+feSGRAi+E39+++/NHz4cC5rhRM+RI4ebNy4kV5//XUe6/79+2nOnDm0bds26tq1q0veD8IG+3v48GH2REAYN2zYkAW3HmCemzdvziIRF7izZs3icJp+/fqZrPfGG29wzLY5e/bsoRs3blDdutbDbLBv58+fp+3bt7M4xWeMNot79+4lPUCHnh9++IG+/fZb9gDiGP7CCy+YXNAL+iMC0sPYd2wH3Q80WLOqlq3v7uH4JKEhIfRV56f5f5xcICY//fRTun3bUD5JBVfwsDoCtGPD1T2ufvft2/dopdQUomnNDQXC1TqPguAACQkJdi9aKxP+x2NI6Mpsu470SA4LC6OaNR92dcoAJJHBCoYaqLBYQmygOgEEnfp7chS0RCxYsGA6VyzeB1ZGtaQWLIHvv/8+W9aeffZZ6tatG4tIV6Ba/dB1qk6dOjxGHE/g4YDwygroZgXr36FDh0weHzt2LFvjAATjU089xe+HYxmE4KhRozheX+thadWqFe3YsYPbq2pZtGiR0UJpDYQuYB9LlixJ7du3Z5EeExPjlD7daPOK74c5uOiAOAe//PILrwehjO8V3hf/f/nll1l+f8FxREB6GHuPb+DbiLQ0KlqgpLuH45sEBFCrEoH09xsRVPOpEmxdGTp0KB+Y3nvvPWrbti0nSsH9hat7w0sC2O2Ggy8OtkaCgolq9SDKV56odCv37ZPg9eTIkcPuBX3WVfA/HmvWrJnJdiGmzF9nL+vXr6cqVao4tF9wfSLM4/Tp07Rz507KCi+99BK7otHhSQXls+DO7NixI9+vVasWd4SC6EWlBHga5s6dy4JDL7C/eG+Is6wAwVa1atV0rmfc79ChA/9///59vrDVAnd6UlKSyXxD4ObLl48/Sy2wVkKA2wO2DysghOTly5cpK+Bzg7jXCltYjmEZzWwfUTJNcB8iID2M+GvH+DaXhD+6ljp9qU77nrTpn8M0f/58evLJJ/nENH78eD4Rw+2GuButBQFX99G5cxP9O4vopKFfOVO1C9F/VxJFP+GefREEFwPxB8ufo6C4PshqAk7u3LlZIM+cOdP4GMQhPARoIQqeeeYZFlgvv/wyJ7nBcobYPVjk9CI6OtppCUcQWL/++quJVRLCUBXMTZo04cYFWAcW3vj4eKPlzjw8B58hPksVrAuhZn7RoednijhZWBu1nyk+P1glccxV9/Grr77iWFpYn1esWMHHbVvCjwTXIQLSw8r43Eq8xLfZ0+SjcSn5yhA1GcaWRcTS7Nu+nqZ+8RFbDr755hu2aJw8eZJPQkzyPaJ/fzPUdlzwFtGiHkT37xqeCwwkCjG9OhYEe7l7967dC767Kvgfj6GwvRac4M1fZy9wi5tbgOwB1jhgkmzmIBBO8+bNY6uUKjbgVoWlEyBGDr9juHQhtPBbxhzAYpaRiFGts46IKWv77Iz9xb5h/Fu2bDHuL+KxVQHXuHFjGj16NO8fwgxgtVStreqcaK128LhorY9w8TuS0OTsz1QVkNguxLAqkAGOyejshX3GRUGPHj2oc+fO6fZP0Jdgnd/P60FGIRYkUnC3ESeT8LCNYYS0MdSP1BQKXvBfeuPeOqJWvYnyFyEKDSIKvEh07qKhk8y/vxIl3TKsHxRKVPk17rEtCM6Ca4hmgeDgYF6cvV0ACx8SLRwFFn2AmMSsglg+iIw//viDqlWrxi5ZxASqjBgxgq2Qffv25fuID8QcoNUoQlWQlW0O3N3wOKgiK6vAm4HEO2fsLyyo9evXZ4GFGFTcmsce9u7dm3r16sUWOVhpITj79+/PYTla4O5H7KJWQKIcWVY+UzUrPyu88sornJyza9cuvlhBCILx4p2Ix7xw4UJ2y2NuYUlFkpD5/gn6IgLSA9sYxnAbQxEnuqGkEeWKNdwiGcYaUbFElV8nqtSRKKfj7jxB8DZQumX69OkOvRYuR2T/QkxhO1kFllDEKcMSd+zYMQ4/gUVOBRY2cyEdFBRkYjUzp2jRouRMYDGDdez55593yvZgjUPpHAgtlA6DVdIcWALVMANY8GDg0M4LxBfiDNXPAJZoxJJ+//33do8HIg/JQkga0gpSRylcuDAn/+AzxbYbNWrEIQCWPntUvoDYhxUaJZME9yEC0sOIzpWX7qdcoyK5JZ5ON4JDiVqNIypcnejg70T37xA9uGNwUcN1XaQaUZU3iIrVIwoUy7DgfyAGDRYtWCFh4coIWIguXrzIQg4VC1A3EUkSsBiqQg4gXs+8LiJcsKVLl7ZJULVs2ZKTLdAFytxCiZI9EEYYN6xyPXv2pOrVq2cpjtMayHTG/kLUIOwFQnvSpElsCVVj+FQsZWXDfZ5RBjSAYIbVEQtiPc33Ay5sJPdBtCI2cOTIkTR79myT+YYLHPOLJCMA1z7c3bZYEJEoAwGKfUVYALK8r169yu+lBeLP/DNFlrqaMZ7ZZzpo0CB68OCBiUUZbN26leM1UWcSt6g5iQsTiGrBjSiCQ5w9exaXsnzrCtLS0lyyXUEQ3Me9e/eUAwcO8K23Ub16dWXixInG+ydPnuRj4O7du03uq0tERIRSunRppXv37srRo0dNtjV16lSTddUlLi7OprGkpqYqBQoU4NccP3483fPjxo1TypQpo4SHh/N6HTt2VM6dO5flOShatKgyduxYk/vq2ENDQ5XY2FilXbt2yurVq01eZz432sXWcwi2i/WnTJmS7rl69eopUVFRSrZs2ZQaNWooS5cuTbfOW2+9pXTr1s14/9VXX1UGDBhgss6aNWv4PW7cuGFyH0tAQIASGRmpVKhQQenbt69y4cIFk9cOGjTI4v41aNDApv3De4aFhfH35s6dOybPrV27lr9LeD5PnjzKa6+9psTHxzv8W3P1+dtfCMAfdwpYb0WNgUSsBszvgiAImQErDqxUcOdmJSnFHcCCiLhCWBUlecG7gLUQrn7UgcR3DzVDUdIHCVewzPrbb03O385BXNiCIAhCprRo0YLLqMCF6IoEQsF1IKlmwoQJxqQeJNMg6QZJSILgKHIZ6WFlfARBEDwVxBK6UjyeOXMmw8LpeN7X0JYQMl+s9a62FxQj12Y1I0EF7QGdUYInM5Aln9FnKngvYoH0sDI+giAI/gqSQ8yTMMyf9zW0JYTMgZvZ24F4zegzFbwXEZCCIAiCR4DyO+aZy76Os0sIeRqoq+lvn6m/IC5sQRAEnZHcRUFwLfIbcz0iIAVBEHRCrfenbScnCILzUX9jmdXYFBxHXNiCIAg6gcLO6DuMwswgIiJCl0QGQfAnyyPEI35j+K1pi6kLzkUEpCAIgo6gtzFQRaQgCM4H4lH9rQmuQQSkIAiCjsDiWKBAAS6lYi37VhAEx4HbWiyPrkcEpCAIghvACU5OcoIgeCuSRCMIgiAIgiDYhQhIQRAEQRAEwS5EQNqJtDIUBEEQBMHfCVCk2qZDoCcrOghs27aNA+IFQRAEQfB8Lly4QNWrV6fTp09TbGysu4fjtUgSjYNcunSJb/ElFARBEATB+87jIiAdRyyQDpKSkkK7d+/mZvf169entWvXWlwPrm7z58wf096/c+cOu8gPHDhAkZGRLt4L62N01ettWTejdaw9Z8s8u3uuZZ5lnm19TubZ/nVlnp3z+szWd9Y8W3pcr3lOS0tj8VipUiXuvy44hsycg+BLV61aNf4/NDSUChcubHE9S8+ZP6a9f/v2bb4tVKgQ5cyZ04V7kPEYXfV6W9a1dz6tPZ7ZY3rPtcyzzLOtz8k827+uzLNzXp/Z+s6aZ0uP6znPYnnMOpJE4wTeffddu54zfyyj1+tBVt/fntfbsq6982ntcVsf0wuZZ32QedYHmWffm2db1nfWPHviuVCwD3Fhexi46oqKiqJbt27pZoH0V2Su9UHmWR9knvVB5lkfZJ49H7FAehhhYWE0aNAgvhVci8y1Psg864PMsz7IPOuDzLPnIxZIQRAEQRAEwS7EAikIgiAIgiDYhQhIQRAEQRAEwS5EQAqCIAiCIAh2IQJSEARBEARBsAsRkIIgCIIgCIJdiID0QhITE6lo0aLUp08fdw/FZ3n88cfpqaeeoooVK1K9evXcPRyf5eTJkzy/aFlWvnx5SkhIcPeQfI7Dhw/z91hdwsPDaeHChe4els8yduxYKlu2LH+n33//fZJCJ65hzJgxPM/lypWj6dOnu3s4fom0MvRChg0bRjVr1nT3MHyeTZs2UY4cOdw9DJ/mjTfeoKFDh1Lt2rXp+vXrUvPNBTz55JP0zz//8P93797li6NGjRq5e1g+yZUrV2j8+PG0f/9+CgkJoTp16tCWLVuoVq1a7h6aT7F3716aOXMm7dy5kwU6LkJbtmxJuXLlcvfQ/AqxQHoZR48epUOHDlGzZs3cPRRByBLqSRbiEURHR3OPecF1LF68mBo0aEDZs2d391B8lpSUFEpKSqLk5GReHnvsMXcPyec4ePAgi/Js2bKxRb1ChQq0bNkydw/L7xABqSPr1q2jVq1aUcGCBSkgIMCiG+m7775jCwF+GDVq1KBt27aZPA+39YgRI3QctX/OM15Xt25dqlatGs2YMUPH0fvPPONiCBZebKNy5co0fPhwnffAf77PKrNnz6aXX35Zh1H751zHxMTwMTo2Npa30bBhQ4qLi9N5L3x/nuG2Xrt2Ld28eZNu3LjB/8fHx+u8F4IISB1BfBeulPDDsMSsWbOod+/e3L5p165dvG6TJk3o8uXL/PyiRYuoZMmSvAium2ewYcMGdo/AYgNhs2fPHh33wD/mGZaa9evX04QJE2jz5s20YsUKXgTnf5/V3sIIy2jevLlOI/e/uYaYWbJkCZ06dYoFDeYbYklw7jyr8aX169entm3bckhXUFCQznshIH5AcAOY+gULFpg8Vr16deXdd9813k9NTVUKFiyojBgxgu/369dPKVy4sFK0aFElT548Ss6cOZUhQ4boPnZfn2dz+vTpo0ydOtXlY/W3ed60aZPSuHFj4/OjRo3iRXDN9/nnn39WOnbsqNtY/XGuZ8+erXTv3t34PL7PX3zxhY6j9s9jdJcuXZQlS5a4fKyCKWKB9BAePHjAFi+4PFQCAwP5PqwzAK7rs2fP8tUtMtC6du1KAwcOdOOofXOecXV8584dY9LB6tWrOdtPcO48IzwAFgVYbdLS0thSU7p0aTeO2jfnWUXc166f6yJFirDVETGQqamp7FpFEpPg/O+0ao1ElQG4t2GhFPRFItY9hKtXr/IBJ1++fCaP4z6SZgT95vnSpUv0wgsv8P9YF0IdYkdw7jwjYQbhAchUhSGicePGnEkpOP+4cevWLT7Jzps3zw2j9J+5hisVIQKVKlVi0YOEpdatW7tpxL79nW7Tpg1/r5EQNnXqVEnAcwMy415c/kRwDcWKFaN///3X3cPwC1BNQCoKuJ6oqCi+MBL0KbOGRXAt5hZ2QX/Ehe0h5M2bl4OAzQ/yuJ8/f363jcvXkHnWB5lnfZB51g+Za32QefYeREB6CKGhoVSlShVatWqV8THEheG+FKF1HjLP+iDzrA8yz/ohc60PMs/eg7iwdQQJGceOHTNp44YOESigjLphKFvQqVMnqlq1KlWvXp2+/vprTujo3LmzW8ftbcg864PMsz7IPOuHzLU+yDz7CGZZ2YILWbNmDZcsMF86depkXOfbb79VYmNjldDQUC5lsGXLFreO2RuRedYHmWd9kHnWD5lrfZB59g0C8MfdIlYQBEEQBEHwHiQGUhAEQRAEQbALEZCCIAiCIAiCXYiAFARBEARBEOxCBKQgCIIgCIJgFyIgBUEQBEEQBLsQASkIgiAIgiDYhQhIQRAEQRAEwS5EQAqCIAiCIAh2IQJSEARBEARBsAsRkIIgCIIgCIJdiIAUBEEQBEEQ7EIEpCAIgiAIgkD28P/Fb7p752jpmgAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 640x480 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"fig = plt.figure()\n",
"ax = fig.add_subplot()\n",
"l1, = ax.plot(idl_res_v8['temperature'], idl_res_v8['free_free_radiative_loss'], label='IDL, v8')\n",
"l2, = ax.plot(idl_res_v9['temperature'], idl_res_v9['free_free_radiative_loss'], label='IDL, v9')\n",
"l3, = ax.plot(idl_res_v11['temperature'], idl_res_v11['free_free_radiative_loss'], label='IDL, v11')\n",
"ax.plot(col_v8.temperature, ff_rad_loss_v8, color=l1.get_color(),ls='--', label='Python, v8')\n",
"ax.plot(col_v9.temperature, ff_rad_loss_v9, color=l2.get_color(),ls='--', label='Python, v9')\n",
"ax.plot(col_v11.temperature, ff_rad_loss_v11, color=l3.get_color(),ls='--', label='Python, v11')\n",
"ax2 = ax.twinx()\n",
"diff_v8 = (ff_rad_loss_v8.value - idl_res_v8['free_free_radiative_loss'])/idl_res_v8['free_free_radiative_loss']\n",
"diff_v9 = (ff_rad_loss_v9.value - idl_res_v9['free_free_radiative_loss'])/idl_res_v9['free_free_radiative_loss']\n",
"diff_v11 = (ff_rad_loss_v11.value - idl_res_v11['free_free_radiative_loss'])/idl_res_v11['free_free_radiative_loss']\n",
"diff_v8_v9 = ((idl_res_v8['free_free_radiative_loss']-idl_res_v9['free_free_radiative_loss'])\n",
" /idl_res_v9['free_free_radiative_loss'])\n",
"ax2.plot(col_v8.temperature, diff_v8, label='(Python - IDL)/IDL, v8', color=l1.get_color(), ls='-.')\n",
"ax2.plot(col_v9.temperature, diff_v9, label='(Python - IDL)/IDL, v9', color=l2.get_color(), ls='-.')\n",
"ax2.plot(col_v11.temperature, diff_v11, label='(Python - IDL)/IDL, v11', color=l3.get_color(), ls='-.')\n",
"ax2.plot(col_v9.temperature, diff_v8_v9, label='(IDL_v8 - IDL_v9)/IDL_v9', color='k', ls='-.')\n",
"ax2.axhline(y=0, color='k', ls=':')\n",
"ax.set_xscale('log')\n",
"ax.set_yscale('log')\n",
"ax.set_ylabel('Free-free Radiative Loss')\n",
"ax2.set_ylabel('Normalized Difference')\n",
"ax.legend(loc=2)\n",
"ax2.legend(loc=4)\n",
"#ax2.set_ylim(-0.001,0.001)"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "e0237996-078b-43b2-810a-5048f08a29af",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python [conda env:fiasco-dev]",
"language": "python",
"name": "conda-env-fiasco-dev-py"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.3"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment