{
"cells": [
{
"cell_type": "markdown",
"id": "0667fa3c",
"metadata": {},
"source": [
"## Importing custom source data\n",
"\n",
"In this example we illustrate how to incorporate a source with a custom spatial dependence in a Tidy3D simulation. We will illustrate this both using data taken from a monitor from another Tidy3D simulation and using a simple array defined from scratch."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "590771ce",
"metadata": {
"execution": {
"iopub.execute_input": "2023-02-03T05:03:51.676941Z",
"iopub.status.busy": "2023-02-03T05:03:51.676653Z",
"iopub.status.idle": "2023-02-03T05:03:52.355750Z",
"shell.execute_reply": "2023-02-03T05:03:52.355406Z"
}
},
"outputs": [
{
"data": {
"text/html": [
"
[23:03:52] WARNING This version of Tidy3D was pip installed from the 'tidy3d-beta' repository on __init__.py:103\n",
" PyPI. Future releases will be uploaded to the 'tidy3d' repository. From now on, \n",
" please use 'pip install tidy3d' instead. \n",
"
\n"
],
"text/plain": [
"\u001b[2;36m[23:03:52]\u001b[0m\u001b[2;36m \u001b[0m\u001b[31mWARNING \u001b[0m This version of Tidy3D was pip installed from the \u001b[32m'tidy3d-beta'\u001b[0m repository on \u001b]8;id=969208;file:///Users/twhughes/Documents/Flexcompute/tidy3d-docs/tidy3d/tidy3d/__init__.py\u001b\\\u001b[2m__init__.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=676435;file:///Users/twhughes/Documents/Flexcompute/tidy3d-docs/tidy3d/tidy3d/__init__.py#103\u001b\\\u001b[2m103\u001b[0m\u001b]8;;\u001b\\\n",
"\u001b[2;36m \u001b[0m PyPI. Future releases will be uploaded to the \u001b[32m'tidy3d'\u001b[0m repository. From now on, \u001b[2m \u001b[0m\n",
"\u001b[2;36m \u001b[0m please use \u001b[32m'pip install tidy3d'\u001b[0m instead. \u001b[2m \u001b[0m\n"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
"
INFO Using client version: 1.9.0rc1 __init__.py:121\n",
"
\n"
],
"text/plain": [
"\u001b[2;36m \u001b[0m\u001b[2;36m \u001b[0m\u001b[34mINFO \u001b[0m Using client version: \u001b[1;36m1.9\u001b[0m.0rc1 \u001b]8;id=959706;file:///Users/twhughes/Documents/Flexcompute/tidy3d-docs/tidy3d/tidy3d/__init__.py\u001b\\\u001b[2m__init__.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=116020;file:///Users/twhughes/Documents/Flexcompute/tidy3d-docs/tidy3d/tidy3d/__init__.py#121\u001b\\\u001b[2m121\u001b[0m\u001b]8;;\u001b\\\n"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# standard python imports\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"\n",
"# tidy3D import\n",
"import tidy3d as td\n",
"from tidy3d import web\n"
]
},
{
"cell_type": "markdown",
"id": "2f985512",
"metadata": {},
"source": [
"### Starting simulation with an in-built source\n",
"\n",
"We will first run a simulation with a [GaussianBeam](https://docs.flexcompute.com/projects/tidy3d/en/v1.9.0rc2/_autosummary/tidy3d.GaussianBeam.html) source propagating in empty space. The definition of the source below is for a converging Gaussian beam (negative `waist_distance` parameter)."
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "cdd2a272",
"metadata": {
"execution": {
"iopub.execute_input": "2023-02-03T05:03:52.357660Z",
"iopub.status.busy": "2023-02-03T05:03:52.357500Z",
"iopub.status.idle": "2023-02-03T05:03:52.360107Z",
"shell.execute_reply": "2023-02-03T05:03:52.359821Z"
}
},
"outputs": [],
"source": [
"# Free-space wavelength (in um) and frequency (in Hz)\n",
"lambda0 = 1.55\n",
"freq0 = td.C_0 / lambda0\n",
"fwidth = freq0 / 10\n",
"\n",
"# Simulation size and run time\n",
"sim_size = [10, 10, 10]\n",
"run_time = 20 / fwidth\n",
"\n",
"# Grid specification\n",
"grid_spec = td.GridSpec.auto(wavelength=lambda0)\n",
"\n",
"# In-built GaussianBeam source\n",
"pulse = td.GaussianPulse(freq0=freq0, fwidth=fwidth)\n",
"waist_radius = 2\n",
"src_pos = -4\n",
"gaussian_source = td.GaussianBeam(\n",
" source_time=pulse,\n",
" center=(src_pos, 0, 0),\n",
" size=(0, 10, 10),\n",
" waist_radius=waist_radius,\n",
" waist_distance=-8,\n",
" direction=\"+\",\n",
")\n"
]
},
{
"cell_type": "markdown",
"id": "635a5b42",
"metadata": {},
"source": [
"We will use monitors to record the fields at the beam waist along the propagation direction, and inject those in another simulation using a [CustomFieldSource](https://docs.flexcompute.com/projects/tidy3d/en/v1.9.0rc2/_autosummary/tidy3d.CustomFieldSource.html). The best way to achieve this is to use a monitor with a slightly nonzero size in the propagation direction. This is because of the details associated with the Yee grid used in the FDTD algorithm. A monitor with a zero size will interpolate the fields to the exact monitor location along the zero-size dimension, which loses a bit of information related to the numerical grid. On the other hand, a monitor with a small but nonzero size along the propagation direction will store the fields exactly as they are in the simulation. We will illustrate the difference in the results below."
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "d21c681b",
"metadata": {
"execution": {
"iopub.execute_input": "2023-02-03T05:03:52.361572Z",
"iopub.status.busy": "2023-02-03T05:03:52.361474Z",
"iopub.status.idle": "2023-02-03T05:03:52.364346Z",
"shell.execute_reply": "2023-02-03T05:03:52.364076Z"
}
},
"outputs": [],
"source": [
"# Monitor for propagation\n",
"mnt_xy = td.FieldMonitor(\n",
" center=(0, 0, 0), size=(10, 10, 0), freqs=[freq0], name=\"field_xy\"\n",
")\n",
"\n",
"# Monitors for forward and backward flux\n",
"mnt_flux_pos = src_pos - 0.5\n",
"mnt_flux_f = td.FluxMonitor(\n",
" center=(-mnt_flux_pos, 0, 0), size=(0, td.inf, td.inf), freqs=[freq0], name=\"flux_f\"\n",
")\n",
"mnt_flux_b = td.FluxMonitor(\n",
" center=(mnt_flux_pos, 0, 0), size=(0, td.inf, td.inf), freqs=[freq0], name=\"flux_b\"\n",
")\n",
"\n",
"# Monitor to be used as custom source, 0D along x\n",
"mnt_yz_1 = td.FieldMonitor(\n",
" center=(-src_pos, 0, 0), size=(0, 8, 8), freqs=[freq0], name=\"yz_zero_size_x\"\n",
")\n",
"\n",
"# Monitor to be used as custom source, small nonzero along x\n",
"mnt_yz_2 = td.FieldMonitor(\n",
" center=(-src_pos, 0, 0), size=(1e-5, 8, 8), freqs=[freq0], name=\"yz_nonzero_size_x\"\n",
")\n"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "9773710c",
"metadata": {
"execution": {
"iopub.execute_input": "2023-02-03T05:03:52.365821Z",
"iopub.status.busy": "2023-02-03T05:03:52.365706Z",
"iopub.status.idle": "2023-02-03T05:03:52.479347Z",
"shell.execute_reply": "2023-02-03T05:03:52.479083Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 1,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcYAAAHWCAYAAADttCmyAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABsoUlEQVR4nO29eZhcZZn3/z2n1k51VbpT3YEOiUDiEhCYQBLwZRlkyAAOLoyKjgIDyATBoGCYUfLixaZDVFAZWcV5DV6IENGBcUBZZBFFUQTGH2sQAc2eVHV3qrrTtZ7n90dxautazvKcpVPfz3XVBV05VXW67/PcnzrnPM99K0IIAUIIIYQAAFSvd4AQQgjxExQjIYQQUgfFSAghhNRBMRJCCCF1UIyEEEJIHRQjIYQQUgfFSAghhNRBMRJCCCF1UIyEEEJIHRQjITOUxx9/HIqi4PHHH/d6VwjZo6AYCfE5N910E2677Tavd8MSP/zhD3Hdddd5vRsAAE3T8PWvfx37778/otEoDjnkENx5552GXz8+Po5zzz0Xw8PDiMViOO644/Dss8+23PanP/0pDjvsMESjUbztbW/D5ZdfjlKpJOtXIU4jCCG+5t3vfrc49thjpz1fLpfF1NSUKJfL7u+UQU4++WSx7777er0bQgghLrnkEgFArFy5Utx6663i5JNPFgDEnXfe2fW15XJZHHnkkSIWi4krrrhC3HDDDeLAAw8U8XhcvPrqqw3b/uxnPxOKoojjjjtO3HrrreKzn/2sUFVVnHfeeU79akQyFCOZEUxMTHi9C57RTowzAb+IcdOmTSIUColVq1ZVn9M0TRxzzDFi/vz5olQqdXz9+vXrBQBx9913V5/bsWOHGBgYEJ/4xCcatj3wwAPF3/zN34hisVh97tJLLxWKooiXX35Z0m9EnIRiJK6zadMm8alPfUqMjIyIcDgs9ttvP3HeeeeJfD4vhBBi3bp1AoB4/PHHxfnnny+Gh4fFwMBA9fU33nijOPDAA0U4HBYjIyPiM5/5jBgbG2v4jFdffVV8+MMfFnvttZeIRCJin332ER//+MfF+Ph4dZuHHnpIHHXUUWL27NkiFouJd77znWLNmjVd99/I63K5nLjsssvEokWLRDgcFvPnzxf/9m//JnK53LT3u/3228Xy5ctFX1+fGBgYEMccc4x48MEHhRBC7LvvvgJAw0OX5GOPPSYAiMcee6zh/X70ox+Jww47TESjUZFMJsVpp50mNm3a1LDNmWeeKWKxmNi0aZP40Ic+JGKxmBgaGhIXX3xxV0kIIcS9994r/uEf/qEaw4ULF4qrrrqq4bXHHnvstH3vJMkzzzxz2vb64/LLL++6T5248cYbBQDx4osvNjz/wx/+UAAQv/rVrzq+/tRTTxV77bXXtLPzc889V8yaNasa1xdffFEAEDfeeGPDdps3bxYAxJe//GVbvwdxh6CLV20JwZYtW3D44YdX79csXrwYmzdvxo9//GPs3r0b4XC4uu1nPvMZDA8P47LLLsPk5CQA4IorrsCVV16JFStW4Pzzz8eGDRtw88034+mnn8aTTz6JUCiEQqGAE088Efl8Hp/97Gex9957Y/PmzbjvvvswPj6O2bNn48UXX8T73/9+HHLIIbjqqqsQiUTw2muv4cknn+y4/0Zep2kaPvjBD+LXv/41zj33XBxwwAF4/vnn8a1vfQuvvvoq7r333uq2V155Ja644goceeSRuOqqqxAOh/G73/0Ojz76KE444QRcd911+OxnP4v+/n5ceumlAIC99tqr7f7ddtttOPvss7F8+XKsXbsW27dvx3/8x3/gySefxHPPPYeBgYHqtuVyGSeeeCKOOOIIXHvttfjFL36Bb3zjG1i0aBHOP//8jn+H2267Df39/Vi9ejX6+/vx6KOP4rLLLkMmk8E111wDALj00kuxa9cubNq0Cd/61rcAAP39/W3f89Of/jRWrFjR8NwDDzyAO+64A3Pnzq0+l0qlOu6bTjweRyQSAQA899xziMViOOCAAxq2Ofzww6v/fvTRR7d9r+eeew6HHXYYVLVxWsbhhx+OW2+9Fa+++ioOPvhgPPfccwCAZcuWNWw3b948zJ8/v/rvxOd4bWbSW/zzP/+zUFVVPP3009P+TdM0IUTtjPHoo49uOAPZsWOHCIfD4oQTTmj45n7DDTcIAOJ73/ueEEKI5557btplr2a+9a1vCQBi586dpvbfyOtuv/12oarqtLOQW265RQAQTz75pBBCiD/96U9CVVXxj//4j9PORPS/hRDtL6U2nzEWCgUxd+5ccdBBB4mpqanqdvfdd58AIC677LLqc/rZ2VVXXdXwnoceeqhYunRp5z+CEGL37t3Tnvv0pz/dcPYkhL1LqX/605/E7Nmzxd///d83HAdoc1bZ/Fi3bl3DfixcuHDaZ0xOTgoA4pJLLum4L7FYTHzqU5+a9vz9998vAIgHHnhACCHENddcIwCIv/71r9O2Xb58uXjPe95j9NcnHsJZqcQ1NE3Dvffeiw984APTvlEDgKIoDT+vXLkSgUCg+vMvfvELFAoFXHTRRQ3f3FeuXIlEIoH7778fADB79mwAwIMPPojdu3e33Bf9zOm///u/oWma4d/ByOvuvvtuHHDAAVi8eDFSqVT18Xd/93cAgMceewwAcO+990LTNFx22WXTzkSa/xZG+MMf/oAdO3bgM5/5DKLRaPX5k08+GYsXL67+feo577zzGn4+5phj8Prrr3f9rL6+vur/Z7NZpFIpHHPMMdi9ezdeeeUV0/vezOTkJP7xH/8Rg4ODuPPOOxuOg4cfftjQ48QTT6y+Zmpqqnr2WI/+d5qamuq4P0Zfr/+33bbdPof4A15KJa6xc+dOZDIZHHTQQYa233///Rt+/stf/gIAeNe73tXwfDgcxsKFC6v/vv/++2P16tX45je/iTvuuAPHHHMMPvjBD+L000+vSvPjH/84/vM//xP/8i//gksuuQTHH388PvzhD+OjH/3oNEnVY+R1f/rTn/Dyyy9jeHi45Xvs2LEDAPDnP/8ZqqriwAMPNPT36Ea7vw8ALF68GL/+9a8bnotGo9P2cXBwEGNjY10/68UXX8SXvvQlPProo8hkMg3/tmvXLrO7Po2VK1fiz3/+M37zm98gmUw2/Fvz5VYj9PX1IZ/PT3s+l8tV/13G6/X/ttu22+cQf0AxEt9iJ4l84xvfwFlnnYX//u//xkMPPYTPfe5zWLt2LZ566inMnz8ffX19eOKJJ/DYY4/h/vvvxwMPPID169fj7/7u7/DQQw81nKE071O312mahoMPPhjf/OY3W77HggULLP9eMmn3O3ZjfHwcxx57LBKJBK666iosWrQI0WgUzz77LL74xS+aOgNvxX/8x3/gzjvvxA9+8AMsWbJk2r9v27bN0PvMnj27egyNjIzgsccegxCi4Wx869atACr3ADsxMjJS3bae5tePjIxUn2+O89atW6v3NIm/4aVU4hrDw8NIJBJ44YUXLL1+3333BQBs2LCh4flCoYA33nij+u86Bx98ML70pS/hiSeewK9+9Sts3rwZt9xyS/XfVVXF8ccfj29+85t46aWX8O///u949NFHq5c629HtdYsWLcLo6CiOP/54rFixYtpDP6NbtGgRNE3DSy+91PHzjF5Wbff30Z9r/vtY5fHHH0c6ncZtt92GCy+8EO9///uxYsUKDA4OTtvW7CXhX/3qV/jXf/1XXHTRRTjttNNabjMyMmLosX79+uprlixZgt27d+Pll19ueK/f/e531X/vxJIlS/Dss89Ok/7vfvc7zJo1C+985zsb3ucPf/hDw3ZbtmzBpk2bun4O8QcUI3ENVVVxyimn4H/+53+mJQ4AEEJ0fP2KFSsQDofx7W9/u2Hb//f//h927dqFk08+GQCQyWSmVRk5+OCDoapq9RLX6OjotPfXk1ary2A6Rl73sY99DJs3b8Z3v/vdadtOTU1VZ9iecsopUFUVV1111bSEW//7xWIxjI+Pt90nnWXLlmHu3Lm45ZZbGn6Hn//853j55Zerfx+76Gea9ftYKBRw0003Tds2FosZvrS6detWfOxjH8PRRx9dndnaCiv3GD/0oQ8hFAo17KMQArfccgv22WcfHHnkkQ378corr6BYLFaf++hHP4rt27fjv/7rv6rPpVIp3H333fjABz5Qvaf47ne/G4sXL8att96Kcrlc3fbmm2+Goij46Ec/auhvQbyFl1KJq1x99dV46KGHcOyxx1aXMmzduhV33303fv3rXzcsJ2hmeHgYa9aswZVXXomTTjoJH/zgB7FhwwbcdNNNWL58OU4//XQAwKOPPooLLrgAp556Kt75zneiVCrh9ttvRyAQwEc+8hEAwFVXXYUnnngCJ598Mvbdd1/s2LEDN910E+bPn99x2r6R151xxhn40Y9+hPPOOw+PPfYYjjrqKJTLZbzyyiv40Y9+hAcffBDLli3D29/+dlx66aX48pe/jGOOOQYf/vCHEYlE8PTTT2PevHlYu3YtAGDp0qW4+eab8ZWvfAVvf/vbMXfu3OpEnnpCoRC+9rWv4eyzz8axxx6LT3ziE9XlGvvttx8+//nPWw1bA0ceeSQGBwdx5pln4nOf+xwURcHtt9/e8ovN0qVLsX79eqxevRrLly9Hf38/PvCBD7R838997nPYuXMnvvCFL+Cuu+5q+LdDDjkEhxxyCABr9xjnz5+Piy66CNdccw2KxSKWL1+Oe++9F7/61a9wxx13NFxWXrNmDb7//e/jjTfewH777QegIsb3vOc9OPvss/HSSy9haGgIN910E8rlMq688sqGz7rmmmvwwQ9+ECeccAL+6Z/+CS+88AJuuOEG/Mu//Mu05SLEp3g4I5b0KH/5y1/EP//zP4vh4WERiUTEwoULxapVq6Yt8G+1pEOIyvKMxYsXi1AoJPbaay9x/vnnNyzwf/3118WnPvUpsWjRIhGNRsWcOXPEcccdJ37xi19Ut3nkkUfEhz70ITFv3jwRDofFvHnzxCc+8Ylp5b2aMfq6QqEgvva1r4l3v/vdIhKJiMHBQbF06VJx5ZVXil27djVs+73vfU8ceuih1e2OPfZY8fDDD1f/fdu2beLkk08W8Xjc0AL/9evXV99vzpw5HRf4N3P55ZcLI2nhySefFO95z3tEX1+fmDdvnvjCF74gHnzwwWn7MzExIT75yU+KgYGBrgv8WxUE0B92F/gLUSnrdvXVV4t9991XhMNh8e53v1v84Ac/mLadvpTljTfeaHh+dHRUnHPOOSKZTIpZs2aJY489tu0xes8994glS5aISCQi5s+fL770pS+JQqFg+3cg7qAI0eX6FSGEENJD8B4jIYQQUgfFSAghhNRBMRJCCCF1UIyEEEJIHRQjIYQQUgfFSAghhNTRUwv8NU3Dli1bEI/HLXUvIIQQMjMRQiCbzWLevHkdGwUAPSbGLVu2+KaAMyGEEPfZuHEj5s+f33GbnhJjPB4HAHzhC19o6NxQKBSq5azC4bAnZ5Plcrla3zMYDBrqfKDv89DQEEKhUPcPEUUg19T9PDoEKKHq++3cuRNApZ+c3qLJbXbt2lWt9Tk8POxJPAqFQrU+6cDAAMLhsOv7wHjUYDxqGI5Hl/HejWKxiFSq8vrmz7CSr2QjhEChUKjuX7djYmpqCl//+terHuhET4lRD25fXx/6+/sBVAo/K4oCVVWrLWkikUjXU22ZlEolFAqFBrmFw2EEg53Do2kaNE1DIpEwlii0AhCabHyuLwGoYWiahnQ6jVgshkgkglwuB1VVDR1EMslmswgEApgzZw7y+TyKxSKSyaSr8SgUCpicnKwmvmKxiNmzZ7uajBmPGoxHDVPx6DDejaD/3VVVbXh/q/lKJpqmIZ/PIxAIQFEUaJqGYDDYskF0M0a+2PX05Jt8Po9SqYRgMIi+vj6Ew2EIIZDP5233lDOKfpApioJoNIpoNApFUVAoFKZ1iHAKfdAXi0UMDQ0hmUwikUggk8kgm826sg9AZdBnMhkkEgkkk0kMDQ2hWCwinU67Fo9CoYBUKoVQKIShoaHq2Xgqlap+O3UaxqMG41HDD/HwS77K5/MQQiAcDqOvrw/BYBClUqljZxwz9KwY66Wof8sIBoOuyrH+INPPUlVVRSQSce1gax70+rfweDzu6uCvH/T6t/BwOOzq4K9Pwvq3cFVVkUwmXUvGjEcNxqOGH+Lhl3xVL0X9LDUSiUiVY0+KsVgsTpOijltybHWQ6bh1sGlCID06Pm3Q67g1+FsNeh23Bn+rJKzjVjJul4R1GI8KjEcNt+Lhi3zVRoo6MuXYk2Lsdj3aaTl2Osh0nD7YNCGQzpRRLJVaDnodpwd/p0Gv4/Tg75SEdZxOxt2SsA7jUYHxqOF0PMrlsvf5qosUdWTJsSfFqAexE07J0YgUm/dT9sFWlWJZYGjOYNdJDE4NfiODXsepwW8kCes4lYyNJmEdxqMC41HDqXhomuZ9vjIoRR0ZcuxJMRpa2gD5cjQjRR3ZB1uDFOMBhMPG/hayB7+ZQa8je/CbScI6spOx2SSsw3hUYDxqyI6HfsvJ03xlUoo6duXYk2I0gyw5WpGijqyDTdO0RimGzIVf1uC3Muh1ZA1+K0lYR1YytpqEdRiPCoxHjYZ4jI5Ds9iHvlAoYGxszPN8ZUWKOnbkSDEawK4c7UhRx+7BpmlaZaKNRSnq2B38dga9jt1kbCcJ69hNxnaTsA7jUYHxqFGNR6mEdKZsWo56PPSF+1YKOsjIV3akqFMvx2KxaPh1FKNBrMpRhhR1mg+2crls6HXVQV8q2ZKijtXBL2PQ61hNxjKSsI7VZCwrCeswHhUYjxrhcBhDcwZRLAtTcqyPx8DAgNR8ZVSOsqSoo8vRzBc2itEEZuUoU4o6zQdbt31oGPRzBm1LUcfs4Jc56HXMJmOZSVjHbDKWnYR1GI8KjEeNcDiEoXjAsBydiocZOcqWoo7py8FSPrWHMCpHJ6SoU3+wdbpEMH3QG5toYxSjg9+JQa9jNBk7Meh1jCZjp5KwDuNRgfGoEQ6pjXL0KB5G5OiUFHWMTroEKEZLdJOjk1LUqT/YxsbGpg1+pwe9TrfB7+Sg1+mWjJ0c9DrdkjHjUYPxqOFKPOrlODruWTw6ydFpKZqFYrRIOzm6IUUdRVEQCAQQDAYbBr9bg16n3eB3Y9DrtEvGbgx6nXbJmPFgPDyPhy7HUsnTeLSSo9+kCPRYdw3Z6AEsFArI5/MIhUKuSVFHVVUMDAwgm80ilUohmUwik8m4Nuh19IGdyWSqz7k16HX0ZJxKpZBOp5FIJJBOp10Z9Dp6Mk6n04wH41HFF/EIqRiaM4jU6C5P4xGJRJDP56tfVIrFoq+kCPCM0Tb6maP+rQcwf6PXLvrgrz9zdHPQ69R/M3Z70Ovoybh+yrnbbZIYjxqMRw1/xCPki3jolcf0q21+kiJAMRJCCCENUIw20e8ptvoW5Bb6PZNSXUFwN/vV6dTfM/GiXx1Qu2ein6mUmu6puAHjUYPxqOGPeBR9EY/mq2tu9nM0wowS4+bNm3H66acjmUyir68PBx98MP7whz94tj/NE23c7ucIVA6y8fFa66hIJOJqvzqd5okEbverA6ZPJIhEIq43c22e2MF4MB6AT+JR1JAaHfM8Hs0Tbdzs52iUGSPGsbExHHXUUQiFQvj5z3+Ol156Cd/4xjcwODjoyf60m33qphyFECiXyw3fhAH3+tXptJtd5+bgbze7zs1mru1mOzIejIfn8ShqSGXLCDXdU3Q7Hq1mn7rd7NgIM0aMX/va17BgwQKsW7cOhx9+OPbff3+ccMIJWLRokev70m1JhhtyrD/IBgent45ya/B3m3LuxuDvNuXcjcHfbQkA41GD8ajhSjx0KQYUJOdML/PmVjw6LcnwmxxnjBh/+tOfYtmyZTj11FMxd+5cHHroofjud7/b8TX5fL46A0x/2MXoOkUn5Vh/kAWDwbYVHaYPfuNFdI1gdB2Wk4Pf6DosJwe/0XVxTidjxqMC41GjQYqJgGfxMLJO0Wk57pFFxF9//XXcfPPNeMc73oEHH3wQ559/Pj73uc/h+9//ftvXrF27FrNnz64+FixYYGsfzC7ed0KOzQdZt31oGPyjYygU5RzwZhcnOzH4zS5OdmLwm10s7lQyZjwqMB41CoVioxS7dMlwKh5mFu87JUez+XfGiFHTNBx22GG4+uqrceihh+Lcc8/FypUrccstt7R9zZo1a7Br167qY+PGjZY/32pFG5lybD7IAoGAoddVB38wiFS2bFuOVit2yBz8Vit2yBz8ViuoyE7GjEcFxqNGoVCoTLQxKEWd+niMj08vH2cGqxVtZMsxn8+jVCrtmUXER0ZGcOCBBzY8d8ABB+Cvf/1r29dEIpHq1Gj9YQW7Zd5kyNFu2SRVVZGcM4BQQLElR7tlrGQMfrtlrGQkY7tlxWQlY8ajAuNRoxqPYNCUFHXql3KUy2UIC82OZeQrGXLUpdjpllPLz7f0aR5w1FFHYcOGDQ3Pvfrqq9h3330d/VxZtU/tyFFWLUFVVZFMBCzLUVZtRzuDX1ZtRzvJWFatTbvJmPGowHjUaIjHnAHTUtQJh8MYHBz0PF/ZkWO9FPU15oY/29TWHvL5z38eTz31FK6++mq89tpr+OEPf4hbb70Vq1atMv1eRm/Cyi4IbkWOsgvsqorSKEeDE3JkFzy2MvhlFzy2koxlF6C2mowZjwqMRw3Z8QiFQggGg97mK4tytCNFYAaJcfny5bjnnntw55134qCDDsKXv/xlXHfddTjttNNMv1d95YV2ONUlw4wcnao63yDH0ektq5pxqguAmcHvVBcAM8nYqa4MZpMx41GB8ajhVDxUVfU+X5mUo10pAjNIjADw/ve/H88//zxyuRxefvllrFy50tL7qKqKUqnUVo5Ot44yIkenW7FU5djUsqoZp1vjGBn8TrfGMZKMnW5VZDQZMx4VGI8aTscjEAh4n68MylGGFIEZJkZZ6JcIWsnRrX6KneToVn8yVaks+G03+N3qF9dp8LvVL65TMnarf1+3ZMx4VGA8argVD1/kqy5ylCVFoEfFCKBa27Rejm42GQZaH2xuN+1sN/jdbKIKtB78bjZRBVonY7eb2jIeNRiPGn6Ih1/yVSs5ypQi0OONivU/YKlUgqZpEEK42mQYaGx2nMvlqs+72Z+suZlrJBJBLpdzvV9cfTPX+ubPbvaLq2+um0qlAGBaLVqnYTxqMB41/BAPv+Sr+mbHxWIRmqZJkyLQw2eMOroEdTG63WQYaPwm5lUna33wCyGQy+UQjUZdb6IKVAZ/NBpFLpeDEML1JqpAZZAnk0kUi0UUi0Ukk0nXm9oyHjUYjxp+iIdf8lUkEoEQApqmNbT9k/L+0t5phlIqlRoWsJqppycLTdMaPlf/BuQ2k5OT1f/Xv425jf5NuNU+uYWmaQ11dTOZDOPRYp/cgvGo4Zd4+CFf1e+DEEJqbdWeFmP9PcW+vr62E3KcpP4afSQSqX4LcrvZcf09k5GREdf71QGN90xGRkY8aeZafw9reHgYw8PDrvar02E8KjAeNfwSDz/kq/p7in19fdJrq/asGFtNtGk1IcdJWt249qLZcfNEArf71QGtJxK43cy11cQON/vV6TAeFRiPGn6Jhx/yVfNEGycKj/ekGDVNazv71C05dprN5ebBls1Otpxd5+bg7zS7zq3B32m2o5vJuN1sR8aD8fA6Hn7IV+1mn8qWY0+KsVgsdpx96rQcjUxxduNgy+4uI5OdaDu7zo3Bb2TKudOD38gSADeScbclAIxHDcajhtPx0HOQ1/mq25IMmXLsSTEC6Dr71Ck5mln34+TBlt1dRmZKQyLe33F2nZOD38w6LKcGv5l1cU4mY6Pr4hiPGoxHDSfj4Yd8ZXSdoiw59qQYjTT4BeTL0cpiWCcOtqoU+1TE47Gu2zsx+K0sTpY9+K0sFnciGZtdLM541GA8ajgRD73tlJf5yuzifRly7EkxKiZasciSo50KETIPtgYpzjLW6BiQO/jtVOyQNfjtVFCRmYytVlBhPGowHjUa42F9KYemaRgfH/c8X1mtaGNXjj0pRrPYlaOMskkyDrZsdtKSFHVkDH4ZZazsJmMZZcVkJGO7ZcUYjxqMR41qPLITyO4um369Hg9dSIGA+VwhI1/Zbh3VJEcz+0AxGsSqHGXWEmw+2Mx01s5ms5WJNhalqGNn8Mus7Wg1GcustWknGcuqtcl41GA8asTjcSTi/chMaabkWB+PwcFBW/tgR47SumTUydFM8RaK0QRm5ehEgV0rB1t10Mf7bUlRx8rgd6Lgsdlk7EQBaivJWHYBasajBuNRIx6PIdGnGpZjczxCoZDtfbCSr2QXBLdSLo5iNIlROTpZdb7+YCuXyx0PtsZB332ijVHMDH4nuwAYTcZOdmUwk4yd6srAeNRgPGrEZwUMydHJeJiRo2wp6ugNlw1vL+2Te4hucnSjFUv9wTY+Pt7yYHO6NY6Rwe9Ga5xuydiNVkVGkjHjUYHxqOFKPOrl2GJCjhvxMCJHp6SoY2bSJcVokXZydLM/WSAQqO5D8+B3q19cp8HvZr+4dsnYzf59nZIx41GB8fAoHrocsxOexaOTHJ2Woll6uh+jXer7OQJAKBRysWmnACCgKgoGB2ZjVyaDdDqF5JxBTE7uRiabRSIeR7x/FiDeuoQitMqj4W202r9bRFWA5JwBpEfHkErtxFByDgAglR6tDPo5A1AVYftzulH5XTVkMrsAoSEWm4X06Fhl0CfnIBwKOL4P4VAAQ8lBpNKj3ePhEIxHjZ6OR9N4j/cpQHiW8XgIDRCVPCOL+n6O+XwekUgExWLRV1IEKEbb1MuxVCpBURR3pFjOAZoABBAqpzE0SyCVncLWzZX2PIk+FfHAOLB7vO5lZaAwPv29FAkTcgAk+wTSxTJ2bt8MAAgHFST7ilBzU7bf3yjxAIBwGZnxnciMA4oCDMUDCJd2APK60nQkDGBoltY9Hg7CeNTo2Xi0GO/xsADCMBaPkgBEESgDUPsAGL8U2Yl6OU5NTVWf84sUgR4Vo97cUhaBQKBhAane+Ng5BKAJaEJACKBYEggFFQRVIF+sfLsLhxQUSk3f9IQAyk0Hd0kAirxvhH1hBbmC9tb/qyiV39pfFwmHFGiTlc+MhBRAafG3cBrFQDxcgPF4i16MR5vxHg6phuJRLFUaEWsAoGmQJUYA0y4bBwIBxztzmFneRjFKeC/9TFFRFAghUCgUnD1jFEJ3I4SmIZ+bQF4oEJqCvoCAJoBMNo9wQEPD/WYhAK3pa2FxAuhyU1oIgaJoXAMUUkLTbmZrmoJCWUH0rRPQ3buBUkBAVd1LQEIAhbKKiFq5hFUuK5jI5hAMuNtItVRWu8fDYRiPGj0ZjxbjXRQnUCgHDMWjpCnQRKSiw7LWNU+YQT+RUFW1kl+KRQSDQVMTZMxCMXZheHgYiUTC9vsIITA2NoZSqYTBwUGEQiFMTk5iYmIC/f39iMXkLY9o/GANyG1DsSSQz03gH961E7P7DFyXKheA3JbG56LzgEDnG+5FrYgdpZ0Nz80NDiOk2l/nRAhxCIvjXWfXVBA/e21fREIqQv17A4qcyUHNObJYLGJsbAzBYBCDg4OOyTGTyRjetifFGAqFbM++0mdzaZqGuXPnVt9Pv7+YyWQQDAadmfEmykCpcvDkhYLZfSXMmf83QDdRlXYDY03fDAeXAMFZHV82VZ7ClqlXGp5L9C1GX6DP7J4TQtzC4niHVgRG/wAACAYUhIIKwuGQlLkI2WwWu3fvxsDAQDU3hsPh6qzdTCbj2AxdMwULelKMduk2xVkPuP4Nxanp4JNTZRS1t75dqaHuYlSD0w9uNdj9daIEoTYdKkY+jxDiHVbHu0N0WiKjL61JpVJIp9OOL1/pBtcxmsTouh+nm4dmd5cxkdMQcvF+ESGEWMHIulGnm0+zVqpDmF0M65Qcs9kJZKY09EdV1ycxEEKIGcwUU3BKjoVCAePj44a3pxgNYrVChGw5VrpkZJHoUxHrs3/NnxBCnMJKhSHZctQrDJlZKUAxGsBu2SRZcqx1yYhL6ZJBCCGdKJWt36qxU3ZPlhzry+4NDAwYfh3F2AVZtQTtyrHxIOu3tA+EEGKGYllgcsp8qToZtWjtyrG5Fi2LiHfB6E1Y2QV2rcrRrYLHhBBSTyigYCKnIZudMPwamfnKqhztFmjvSTGOj493bR7qVNV5s3KkFAkhXhEMKOiPqshks57lK7NylNG1pCfFGAwGOzYPdboVi1E5UoqEEK+J9QWQiMc9zVdG5SirlVdPinFgYKBt81C3+pN1kyOlSAjxC/F4v+f5qpscZfa37EkxKorSsnmom007gfZypBQJIX7DD/mqnRxlN33u2ZJwemftdDqNVCqFZDKJTCbjmhR1msvH6f9PKRJC/IYf8lVz+bhEIoF0Oi1NikAPixGoyTGVSiGVSgGodN5wS4o6zQcbpUgI8St+yFe6HHfu3Cn1TFGnJy+lEkIIIe3oaTHq9xRLpVL18mmn2apOUX+N3snC44QQYhc/5Cv9nqJ+5lgqlaTWVu1ZMTZPtIlEIi0n5DhN841rp7tyEEKIVfyQr5on2kQiEemFx3tSjEKIlrNP9XuObsmx3WwuypEQ4jf8kK/azT6VXXi8J8U4Pj7edvapW3LsNsWZciSE+IVsdsLzfNVtSYZMOfakGOvvKbbCaTkaXfdDORJCvGZyqlxpdedhvjK6TlGWHHtSjAMDA12XZDglR7OLYSlHQohXlMoCEzmt0urOo3xldvG+DDn2pBhDoZCh7WTL0WqFCMqREOIFxbJAf1Q11epOZr6yWtHGrhxnrBi/+tWvQlEUXHTRRY5+jiw52i2b1HiwGW8BQwghVgkFFMT6zDdFlyFHu2XemuUohPGmyzNSjE8//TS+853v4JBDDnHl8+zKUVYtwerBls0iu9t881BCCDFDMGC8uW8zduQoq/ZpvRzHx8cNv27GiXFiYgKnnXYavvvd72JwcNC1z7UqR9kFduPxeKUFzJRmqbM2IYS4hRU5yi4IXl8EwCgzToyrVq3CySefjBUrVnTdNp/PI5PJNDzsYFaOTlWdj8f7kehTMZHTUCrPuBASQnoIM3KULUWdcDiMgYEBw9vPqCLid911F5599lk8/fTThrZfu3YtrrzySqn70NyVo92yD6dbscRnBVAqCxQK1i91EEJmALmdwOafNj63zweB6LA3+2OB5sLjrXKiU1LUMTrpEphBZ4wbN27EhRdeiDvuuAPRaNTQa9asWYNdu3ZVHxs3bpSyL93OHN3qTxbrCyCkGr+hTAghXtHpzNFpKZplxpwxPvPMM9ixYwcOO+yw6nPlchlPPPEEbrjhBuTzeQQCjbOnIpEIIpGII/vT7szR7SbDwYCcormEEOI0rc4c/SZFYAaJ8fjjj8fzzz/f8NzZZ5+NxYsX44tf/OI0KbpBsxwjkQhyuRz7KRJCSBvq5VgoFJDP530lRWAGiTEej+Oggw5qeC4WiyGZTE573k10OW7duhW5XA7RaJRSJISQDuhnirlcDgB8JUVgBt1j9DOTk5PV/8/n8673cySEkJmEfqaoU59D/cCMOWNsxeOPP+71LjTcU4zFYl1nqxJCSC/TfE9xcnKy42xVL+AZow2aJ9q43c8RADSNyzUIITODVhNt/FgLmmK0SLvZp27KsVjUUChTjIQQ/9Np9qnf5EgxWqDbkgw35FgoahibLEOhFwkhPsfIkgyn5bjHFxH3EqPrFJ2UY6FQQCpbRjCgIMx1jIQQH2NmnaJTctQ0bc8uIu4lZhfvOyHHQqGAVHoUoYCCwf4AzxgJIb7FyuJ92XLUNA3pdHrPLiLuFVYr2siUY8NBlghAoRUJIQ5isr9vA3Yq2siSoy7FYrFoqog4xWgAu2XeZMix4SCbMwiVUiSEOEyhpKFYNG9HGWXe7MqxXopDQ0N7ZhFxmZi5CSur9qkdOfqxliAhZM9HURSMTZY9y1dW5dgsRbNrynsyw46Pj0MzcI1AdkFwK3KkFAkhXhEOKggGFKTSo57lK7NytCtFoEfFWCqVkE6nO8rRqS4ZZuRIKRJCvERRgMH+gOf5yqgcZUgR6FExDgwMoFgstpWj062jjMiRUiSE+AFFUZCcM+h5vuomR1lSBHpUjKFQCENDQy3l6FY/xU5ypBQJIX7CL/mqnRxlShHoUTECQDgcniZHt5sMtzrYKEVCiB/xS75qlqNsKQIzvLuGXXQ5plIpbN26FQBcbzJc3+x4586d1f2iFAkhfsMv+aq+2XEmk4GiKFI7GvV85g2Hw4hEItWfY7GY6/ugqioSiUT150QiQSkSQnyJX/JVfa6ORCJS2/z1fPbNZrPI5XKIRqNQFKXrbFUnKBQKSKfTCIVCCIVCSKfTbHZMCPElfshX+uVTRVEQjUaRy+Wk1lbtaTHW31NMJpNtJ+Q4Sf01+qGhoWqFBrf6ORJCiFH8kK+a7ykmk0nphcd7VoytJtq0mpDjJK1uXHvR7JgQQrrhh3zVbqKN7MLjPSnGycnJtrNP3ZJjp9lclCMhxE/4IV91m30qU449K8ZOs0+dlqORKc6UIyHEDxSLWqXVnYf5yuiSDFly7EkxxmKxrksynJKjmXU/lCMhxEs0DRibLHuar8yuU5Qhx54VoxFky9HKYljKkRDiFYWShmCgUhLOi3xldfG+XTn2pBjNIEuOdipEUI6EEC9QFAWD/QFP8pXdijZ25EgxGsCuHGWUTWo42NKjKFhoHkoIIWYIBxUoFpqi25WjrDJv9XKcnJw0/DqK0SBW5SizlmDDwZYtW+qsTQghRrHgxCpW5Si79qkuR4rRIczK0YkCu6qqVlrABCqdtTXNxpFLCCEOYlaOThQEBypyNFPuk2I0iVE5Oll1XlVVJBMBBAMKCmWKkRDiX4zK0Skp6lCMDtNNjm60YlHfuilu51IHIYS4QTc5Oi1Fs1CMFmknRzf7kymKgnCA9xkJIf6nnRz9JkWAYrRFsxzz+bzrTTt5xkgImSk0yzGfz/tOigDFaBtdjvqZYjAYZJNhQghpgy7HYDBYPXP0kxQBipEQQghpgGK0iX6mqJ85lkolT5odE0LITEC/p1gqlapnin6r6EUx2qB5ok0kEnG92bEQjn8EIYRIoXmiTSQS8WW5S4rRIu1mn7rZ7FgIgUKZISSE+J92s0/9WAuaWdUC3ZZkuCFHTQiMTZR5xkgI8T3dlmT4TY4Uo0mMrlN0Uo6apiGdKaNUFggHaEZCiH8xuk7RaTmyVqpDmF2874QcNU1DenQMxbLAYCwAVaUYCSH+xOzifafkmM1mKUYnsFrRRqYcGw6yeAChEMNHCHEOO7dqrFa0kS3HbDaLTCbDWqmysVvmTYYcGw6y5ByEKUVCiMMUSgLCgh3tlnmTJUddiolEgmKUiazap3bk6MdagoSQPR/x1iQ/L/KVXTnWSzEej5v7bFNb7yEYvdYsuyC4FTlSioQQrwgHVZTKAunRMU/ylVU52pEi0MNizGazHbdxqkuGGTlSioQQL1FVYDAW8DRfmZWjXSkCPSrGWCyGTCbTVo5Ot44yIkdKkRDiB0IhFUPJOZ7mK6NylCFFoIfFmEgkWsrRrX6KneRIKRJC/IQf8lU3OcqSIjCDxLh27VosX74c8Xgcc+fOxSmnnIINGzZYfr94PD5Njm42GQZaH2yUIiHEj/ghX7WTo0wpAjNIjL/85S+xatUqPPXUU3j44YdRLBZxwgknmFq02Uy9HNPptOtNhoHGgy2VSiGVSlGKhBBf4od81SzHdDotVYoAEJTyLi7wwAMPNPx82223Ye7cuXjmmWfwt3/7t5bfNx6Po1AoIJfLAYAnTYbD4TCSySRSqRQAUIqEEN/ih3yly3Hr1q3I5XKIRqPSpAjMoDPGZnbt2gUAmDNnTttt8vk8MplMw6OZQqGAfD5f/dnOGahVNE1r2LdMJsN+joQQX+KXfFWfq/P5vNTaqjNSjJqm4aKLLsJRRx2Fgw46qO12a9euxezZs6uPBQsWNPx7/T3FkZGRthNynKT+Gv3w8DCGh4dd7edICCFG8Uu+qr+nODIyIr226owU46pVq/DCCy/grrvu6rjdmjVrsGvXrupj48aN1X9rNdGm1YQcJ2l149rNfo6EEGIUv+Sr5ok2ThQen3FivOCCC3Dffffhsccew/z58ztuG4lEkEgkGh4AqjeOW020cUuOnWZz+VWOoqRBaOzmQUiv4Zd81W72qWw5zhgxCiFwwQUX4J577sGjjz6K/fff3/J7jY+Pd5x96rQcjUxx9pscCztyeOP/PofxR7d6uh+EEHcRolISzut81W1Jhkw5zhgxrlq1Cj/4wQ/wwx/+EPF4HNu2bcO2bdswNTVl+r2CwWDX2adOydHMuh8/yXH3C2Mo7ypg4tlRS9X2CSEzDyGAsYmy5/nK6DpFWXKcMWK8+eabsWvXLrz3ve/FyMhI9bF+/XrT7zUwMGBoSYZsOVpZDOsXOU48NwqhCRS27EZxe86TfSCEuEuhJFAqi0qrO4/yldnF+zLkOGPWMco8S1EUxfC2eiD06clW18rYqRChH2z6YtbknAFXv9FMbs8gv3E3AokQtKkyJl8YA6xfySaEzBCEEBiMBezlKxtrw61WtNHlqBduGRoaMvW5M+aM0UvsnjnKKJvU8E1sdAyai5czd/x/m6EVylCjASiKwsuphPQI4aCKkIWm6DLOHG23jmo6cywWi8Zfa/rTehSrcpRZS7DhYMuUXZPTtj/8FVAqZ9pqXwCFLbsxuXV6sQRCyJ6FnSJgduQorUtGnRzHx8eNv87yJ/YgZuXoRIHdcDhcaQFTrnTWdtqNE9t3IfPXUQSilavuSkSFVtCw4//b7OwHE0JmPFbkKL0g+FtyDAaN3zmkGE1iVI5OVp0Ph8MYigdQKgsUys6GcNsfN6KcL0GJVj5HURQoioJtT/+Fl1MJIV0xI0fZUtRRVRUDAwPGt5f2yT1ENzm60YolHFIxGAs4fsa4+enXgbdkqKP2BZDdPI7s1nFnP5wQskdgRI5OSVHHzKRLitEi7eToZn+yUEhFOOCcGSe278L4X9II9YUanlciKsqFErb/cWObVxJCSCOd5Oi0FM1CMdqgWY5eNBlWVefEqF9G1e8v6uiXUzf9/nVeTiWEGKaVHP0mRWAGrWP0K/XrHDOZDBRF2WP6KVYuo751CaLJf8FZIWQ2jSG7dRyJeYPe7CAhZMZRv85x69ZKiUk/SRHgGaMUYrFY9f8jkcgeIcXaZdTWv0sgEuTlVEKIJcLhMCKRSPXn+hzqByhGm+iXTxVFQTQaRS6Xc7Wfo1Pol1GD0VDLf+flVEKIVbLZLHK5HKLRKBRF8bwWdDMUow2a7ykmk0lPmh07QcNl1DYEZ4Wrl1MJIcQI9fcUk8mkL2pBN0MxWqTdRBu3mx2XHFjH2O0yqk6Ql1MJISZoNdHGL40S6qEYLdBt9qlbcpycKqOoGV+bY5Rul1F1eDmVEGKUTrNP/SZHitEkRpdkOC3H7O4yJnIaQg4s1zByGVWHl1MJId0wsiTDaTmyiLhDmF2n6JQcs9kJZKY09EdVBANyDx6jl1F1eDmVENIJM+sUnZJjoVBgEXEnsLp4X7Ycs9ksMtksEn0qYn0B2+/XjNHLqDq8nEoIaYeVxfuy5VgoFJBKpVhEXDZ2K9rIkmP1IIvHEZ8lX4qAucuoOrycSsieSals/cuunYo2suSoSzEUCrGIuExklXmzK8fGg6zf0j50w+xlVB1eTiVkz6RYFpicKpt+nYwyb3blWC/FZDLJIuLdMHoTVnbtU6tydKuWoNnLqDq8nErInkkooGAipyGbnTD8Gpn5yqocm6Womuy43JNiHB8fR6FQ6LiNUwXBzcrRzQK7Vi6j6vByKiF7HsGAgv6oikw261m+MitHu1IEelSMwWAQqVSqrRyd7pJhVI5uSlG/jKoGAijliw2PVmeBWrHcsI0CoJQr8nIqIXsYsb4AEvG4p/nKqBxlSBHoUTEODAwgFAq1lKNbraO6ydHtVizlQhmBUACKqqBcKFcfpVwRxWy+YVtR0pDP5Bq2KxfLCEZDKBfN348ghPibeLzf83zVTY6ypAj0aNspRVGQTCaRTqeRSqWqAnS7n2J9y6r6n73oTzZ7wRwc/+WPQCs3iu33Nz2Ksb/sRP13KKEBgVAA/+eiExCb27h/fYP+qpJPCJGDH/JVfcuqdDpdFaBMKQI9KkYAUFW1QY7JZBKZTMbVJsPA9INN/38v+pP1zZkuNTXQ/gCLDs5CbDjh5C4RQnyEH/JVsxwTiQTS6bQ0KQI9LEagJsdUKoVUKgUAGB4edr2fYvPB5remnYQQouOHfKXLcefOnVLPFHV68h4jIYQQ0o6eFqN+T7FUKlUvn3aareoU9dfo95R+joSQPRM/5Cv9nqJ+5lgqlaTWVu1ZMTZPtIlEIkgmk21nqzpF841rt/s5EkKIUfyQr5on2kQiEemFx3tSjEKIlrNP9XuObsmx3WwuypEQ4jf8kK/azT6VXXi8J8U4Pj7edvapW3LsNsWZciSE+IVsdsLzfNVtSYZMOfakGOvvKbbCaTkaXfdDORJCvGZyqlxpdedhvjK6TlGWHHtSjAMDA12XZDglR7OLYSlHQohXlMoCEzmt0urOo3xldvG+DDn2pBhDIWPdI2TL0WqFCMqREOIFxbJAf1Q11epOZr6yWtHGrhx7UoxmkCVHu2WTGg824y1gCCHEKqGAglif+aboMuRot8xbsxzNtMSjGA1gV46yaglWD7ZsFtndLNZNCHGWYMB8CzodO3KUVfu0Xo7j4+OGX0cxGsSqHGUX2I3H45UWMFOapc7ahBDiFlbkKLsgeH0RAKNQjCYwK0enqs7H4/1I9KmYyGkolRlCQoh/MSNH2VLUCYfDGBgYMLw9s6pJjMrR6VYs8VkB9EdVFDXrlzoIIcQNjMjRKSnqGJ10CVCMlugmR7f6k8X6Agipxm8oE0KIV3SSo9NSNAvFaJF2cnS7yXAwIKdoLiGEOE0rOfpNikCP92O0S3Oz40gkglwux36Kkvjj2Eb8ZXIU793rnUiE+rzeHUKIBOr7ORYKBeTzeV9JEaAYbaPLcevWrcjlcohGo5SiJP7lqduxafcYBsKz8Hd7vQsnzns3JUnIHkA8HkehUEAulwMAX0kRoBilMDk5Wf3/fD6PQqHQteQc6Y4mBEJqAPlyEf+18Tncs/F/MTvcR0kSMsPRzxR1JicnfXVCQTHapP6eYiwWq15W7VSknBgnoKiYE4kBAEpaGZlijpIkZAbTfE9xcnISmUwGAHwjR4rRBq0m2tTfc3RDjloPLdcIqgFKkpAZTKuJNvX3HAF/yNE/F3UNcuONN2K//fZDNBrFEUccgd///vee7Ee72aduNjsuFjUUyr0jxnp0Se4zawBzo/3Vy63n/e4OHPHAV7Hq9z/ETzf9EZnilNe7SghB59mnfmuUMKPOGNevX4/Vq1fjlltuwRFHHIHrrrsOJ554IjZs2IC5c+e6th/dlmQ0z1Z14syxUNQwNllGxKWvNgLAlt1Z5MpFBEtApBTAyY/dgLHZzi0XmSzloSrdxc8zSUL8jZElGU6fOZopIm5ajGeeeSbOOecc/O3f/q3Zl9rmm9/8JlauXImzzz4bAHDLLbfg/vvvx/e+9z1ccsklruyD0XWKTsqxUCgglS0jGFAQVt1axygwVSpCEwIBoUAIgVy5iKmyc/VaA4qC2SZFRkkS4i/MrFN0So6appkqIm5ajLt27cKKFSuw77774uyzz8aZZ56JffbZx+zbmKZQKOCZZ57BmjVrqs+pqooVK1bgt7/9reOfD5hfvO+EHAuFAlLpUYQCChKzVORduFIoNAFRFAiUgQAUhIQKVVExOxSFCPu38o5RSf793gegL8iJUoTIxsrifdly1DQN6XTa2SLi9957LzZv3ozzzz8f69evx3777Yf3ve99+PGPf4xisWj27QyTSqVQLpex1157NTy/1157Ydu2bS1fk8/nkclkGh5WsVrRRuY9x4aDLBGAYuAyo1365sSgBlRoU2X0FVVEiyqCZQX5sEBhBrmk/p5kMhLDzlwWt7/xFP712R/j5V1bvd49QnyJyf6+DdipaCPrnqMuxWKxaKqIuKV7jMPDw1i9ejVWr16NZ599FuvWrcMZZ5yB/v5+nH766fjMZz6Dd7zjHVbeWipr167FlVdeaft97JZ5k3Hm2HCQzRmAmptC5c6fsxz2qb/F2058F17L/Rnf+N2D0IRALBiCNiuC3bMc/3hplIWGbDGPyVJl7dRgeBb+z/BCvG/eQXjX7L093jtC/EmhpCFUBMx+B5ZR5s3umWO9FIeGhjA1Zfzymq3JN1u3bsXDDz+Mhx9+GIFAAP/wD/+A559/HgceeCC+/vWv4/Of/7ydt29gaGgIgUAA27dvb3h++/bt2Hvv1oltzZo1WL16dfXnTCaDBQsWmLoJK6v2qR05TjvIFPcuXwajIczeL4nI7m3Y9qcSpkpFRAMlhNQy1IKzZ6yKoiAejFg+M26WYX8wgr8fOQDvm3cQjt97MeZGEzJ3l5A9DkVRMDZZRiheQDhi7J68zNqnVuXYLMVwOOysGIvFIn76059i3bp1eOihh3DIIYfgoosuwic/+UkkEpVEc8899+BTn/qUVDGGw2EsXboUjzzyCE455RQAlV/+kUcewQUXXNDyNZFIBJFIZNrz4+PjSCQSXQMmuyC4FTm2PMiE+w2KFSh45+wkXt2VBgCoUBy/lJsvl6AJDQNh46emlCEh8ggHFZQ1Ban0KIaGhq3lK5uYlWMrKZrFtBhHRkagaRo+8YlP4Pe//z2WLFkybZvjjjvO1PVco6xevRpnnnkmli1bhsMPPxzXXXcdJicnq7NUjVIqlZBOpzsGzqkuGWbk6Keq86qiYO2yFSiLyk2HxdHF6As4O6vzPQ9+FWOF3V23owwJcQZFAQb7A8gUQ57mK6NylCFFwIIYv/Wtb+HUU09FNBptu83AwADeeOMNSzvUiY9//OPYuXMnLrvsMmzbtg1LlizBAw88MG1CTjcGBgZQKBTaytHp1lFG5OgnKeqoigJVCQAAwoEgwgFnl8EqaH9GShkS4g6KoiA5ZxDp0XFP81U3OcqSImBBjGeccYblD5PBBRdc0PbSqVFCoRASiQRSqdQ0ObrVT7GTHP0oRT9AGRLiDX7JV+3kKFOKwAyrfCOTcDiMoaGhBjnqxWzd6qfY6mADQCnWoQkN44UpypAQj/FLvmqWo968QZYUgR4WI9Aox61bK2vZ3G4yXH+w7dy5s7pflGLl/kZOKyGsBilDQnyAX/JVvRwzmQwURZFaXaynxQhUghqJRKoNM2OxmOv7oKpq9dIuAEMzZnuB65f9E96cTFOGhPgIv+SrWCxWPWuMRCJS61H3fPbNZrPI5XKIRqNQFAXpdBqanXIPFtAnAoVCIYRCIaTTaUe7cswU/s/wInxiv8MpRUJ8hB/ylX5PUVEURKNR5HI5qV05elqM9RNtkskkhoaGUCwWXZVj/Y3roaEhDA0NudKyihBCzOKHfNU80SaZTEpvWdWzYmw1+1S/5+iWHFvN5nKznyMhhBjFD/mq3exT2f0ce1KMnWafuiXHTlOcKUdCiJ/wQ77qtiRDphx7VoydZp86LUcj634oR0KIHygWtUqrOw/zldF1irLk2JNijMViXZdkOCVHM4thKUdCiJdoGipFxD3MV2YX78uQY8+K0Qiy5WilQgTlSAjxikJJQzBQKQnnRb6yWtHGrhx7UoxmkCVHO2WTKEdCiBcoioLB/oAn+cpumTc7cqQYDWBXjjJqCTYcbOlRFIrurrUkhPQe4aC19nJ25Sir9mm9HCcnJw2/jmI0iFU5yiyw23CwZcsoUo6EEAex03LVqhxlFwTX5UgxOoRZOTpRdV5VVSTnDCIUqHTW1jRnmwUTQohVzMpRthR14vG4qXKfFKNJjMrRyVYsqqoimQggGFBQKFOMhBD/YlSOTklRh2J0mG5ydKM/mfrWTXE7lzoIIcQNusnRaSmahWK0SDs5utm0U1EUhAO8z0gI8T/t5Og3KQIUoy2a5ZjP511v2skzRkLITKFZjvl83ndSBChG2+hy1M8Ug8EgmwwTQkgbdDkGg8HqmaOfpAhQjIQQQkgDFKNN9DNF/cyxVCp50uyYEEJmAvo9xVKpVD1T9FtFL4rRBs0TbSKRiOvNjoVw/CMIIUQKzRNtIpGIL8tdUowWaTf71M1mx0IIFMoMISHE/7SbferHWtDMqhbotiTDDTlqQmBsoswzRkKI7+m2JMNvcqQYTWJ0naKTctQ0DelMGaWyQDhAMxJC/IvRdYpOy5G1Uh3C7OJ9J+SoaRrSo2MolgUGYwGoKsVICPEnZhfvOyXHbDZLMTqB1Yo2MuXYcJDFAwiFGD5CiHPYuVVjtaKNbDlms1lkMhnWSpWN3TJvMuTYcJAl5yBMKRJCHKZQEhAW7Gi3zJssOepSTCQSFKNMZNU+tSNHP9YSJITs+Yi3Jvl5ka/syrFeivF43Nxnm9p6D8HotWbZBcGtyJFSJIR4RTioolQWSI+OeZKvrMrRjhSBHhZjNpvtuI1TXTLMyJFSJIR4iaoCg7GAp/nKrBztShHoUTHGYjFkMpm2cnS6dZQROVKKhBA/EAqpGErO8TRfGZWjDCkCPSzGRCLRUo5u9VPsJEdKkRDiJ/yQr7rJUZYUgR4VIwDE4/FpcnSzyTDQ+mCjFAkhfsQP+aqdHGVKEQCCtt9hBqP/ATOZDAqFAvL5vKtNhoHawZZKpZBKpQCgoeo8IYT4BT/kK12O6XQaqVQKkUgEuVxOmhSBHj5j1InH44hGo8jlchBCeNJkOBwOI5lMolgsolgsIplMUoqEEF/ih3yly1EIgVwuh2g0Kk2KAMVYPVPUMVM2SBaapiGTyVR/zmQy7OdICPElfslX9bk6n89Lra3a02Ksv6c4MjLSdkKOk9Rfox8eHsbw8LCr/RwJIcQofslX9fcUR0ZGpNdW7Vkxtppo02pCjpO0unHtZj9HQggxil/yVfNEGycKj/ekGIvFYtvZp27JsdNsLsqREOIn/JKv2s0+lS3HnhTj+Ph4x9mnTsvRyBRnypEQ4geEqJSE8zpfdVuSIVOOPSnGYDDYdfapU3I0s+6HciSEeIkQwNhE2fN8ZXSdoiw59qQYBwYGDC3JkC1HK4thKUdCiFcUSgKlsqi0uvMoX5ldvC9Djj0pRkVRDG8rS452KkRQjoQQLxBCYDAW8CxfWa1oY1eOPSlGs9iVo4yySQ0H2+gYNDuttQkhxADhoIqQhaboMuRou3VUkxyLxaLx15r+NA948803cc4552D//fdHX18fFi1ahMsvv1zqgs5uWJWjzFqCDQdbpmypszYhhBjFThEwO3KU1iWjTo7j4+OGXzcjaqW+8sor0DQN3/nOd/D2t78dL7zwAlauXInJyUlce+21ru1HfW3V+p/b4USB3XA4jKHkHKR2bMbYRBl9M+KrDSGkF6mvrZpOpw2V3JRdEFyXo5mqZjNCjCeddBJOOumk6s8LFy7Ehg0bcPPNN7sqRsC4HJ2sOh8OhzEUD2DHrhIKgmYkhPgXM3KULUUdVVUxMDBgfHtpn+wyu3btwpw5czz57G6XVd1oxRIOqRiMBcCrqYQQv2PksqpTUtQxM+lyRpwxNvPaa6/h+uuv73q2mM/nGwqE1xe+tUu7M0c3+5OFQirCAZqREOJ/Op05Oi1Fs3h6xnjJJZdAUZSOj1deeaXhNZs3b8ZJJ52EU089FStXruz4/mvXrsXs2bOrjwULFkjd/+YzRy+aDKsqxUgImRm0OnP0mxQBj88YL774Ypx11lkdt1m4cGH1/7ds2YLjjjsORx55JG699dau779mzRqsXr26+nMmk3FEjvp7ZzIZKIrCJsOEENKG+jPHrVu3AoCvpAh4LEa9bYkRNm/ejOOOOw5Lly7FunXrDFWuiUQiiEQidnezK7FYrHpJNRKJUIqEENKBcDiMSCSCXC4HoJJD/cSMmHyzefNmvPe978Xb3vY2XHvttdi5cye2bduGbdu2eb1r1cuniqIgGo0il8u52s+REEJmGtlsFrlcDtFoFIqi+K6i14yYfPPwww/jtddew2uvvYb58+c3/JuXi9xb3VPUr5cD3dc5EkJIr9F8T1HvjWt0naMbeL8HBjjrrLMghGj58Ip2E23cbnZcKs+IEBJCSMuJNn6sBc2saoFus0/dkuPkVBlFzfjaHEII8YpOs0/9JkeK0SRGl2Q4Lcfs7jImchpCXK5BCPE5RpZkOC3HPa6IuF8wu07RKTlmsxPITGnoj6oIBry/7EAIIe0ws07RKTkWCgVTRcQpRoNYXbwvW47ZbBaZbBaJPhWxvoDt9yOEEKewsnhfthz1yT3BoPG5phSjAexWtJElx+pBFo8jPotSJIQ4S6ls/VaNnYo2suSoSzEUCvVGEXG3kFXmza4cGw+yfkv7QAghZiiWBSanyqZfJ6PMm1051ksxmUyaKiLek2I0ehNWdu1Tq3L0Yy1BQsieTyigYCKnIZudMPwamfnKqhybpWh2bWRPinF8fByFQqHjNk4VBDcrR0qREOIVwYCC/qiKTDbrWb4yK0e7UgR6VIzBYBCpVKqtHJ3ukmFUjpQiIcRrYn0BJOJxT/OVUTnKkCLQo2IcGBhAKBRqKUe3Wkd1kyOlSAjxC/F4v+f5qpscZUkR6FExKoqCZDI5TY5u91NsJ0dKkRDiN/yQr9rJUaYUgRlSRNwJVFVFMplEOp1GKpVCMplEJpNxtckw0NjPUYdSJIT4ET/kq/p+jul0GolEAul0WpoUgR4WI1CTYyqVQiqVAlDpEel2P8Xmg41SJIT4FT/kK12OO3fulHqmqNOTl1IJIYSQdvS0GPV7iqVSqXr5tNNsVaeov0bvZssqQggxix/ylX5PUT9zLJVKUmur9qwYmyfaRCKRlhNynKb5xrXb/RwJIcQofshXzRNtIpGI9MLjPSlGIUTL2af6PUe35NhuNhflSAjxG37IV+1mn8ouPN6TYhwfH287+9QtOXab4kw5EkL8QjY74Xm+6rYkQ6Yce1KM9fcUW+G0HI2u+6EcCSFeMzlVrrS68zBfGV2nKEuOPSnGgYGBrksynJKj2cWwlCMhxCtKZYGJnFZpdedRvjK7eF+GHHtSjKFQyNB2suVotUIE5UgI8YJiWaA/qppqdSczX1mtaGNXjj0pRjPIkqPdskmNB5vxFjCEEGKVUEBBrM98U3QZcrRb5q1ZjkIYb7pMMRrArhxl1RKsHmzZLLK7zTcPJYQQMwQDxpv7NmNHjrJqn9bLcXx83PDrKEaDWJWj7AK78Xi80gJmSrPUWZsQQtzCihxlFwSvLwJgFIrRBGbl6FTV+Xi8H4k+FRM5DaUyQ0gI8S9m5ChbijrhcBgDAwOGt2dWNYlROTrdiiU+K4D+qIqiZv1SByGEuIEROTolRR2jky4BitES3eToVn+yWF8AIdX4DWVCCPGKTnJ0WopmoRgt0k6ObjcZDgbkFM0lhBCnaSVHv0kR6PF+jHZpbnYciUSQy+XYT5EQQtpQ38+xUCggn8/7SooAzxhto8tRCIFcLodoNEopEkJIB+LxOKLRKHK5HIQQvpIiQDFKYXJysvr/+Xze9X6OhBAyk9DPFHXqc6gfoBhtUn9PcWRkxPV+joQQMpOov6c4MjLiy3KXFKMNmifauN3PEQA0LtcghMwQWk208WMtaIrRIu1mn7opx2JRQ6FMMRJC/E+n2ad+kyPFaIFuSzLckGOhqGFssgyFXiSE+BwjSzKcliOLiDuI0XWKTsqxUCgglS0jGFAQ5jpGQoiPMbNO0Sk5aprGIuJOYXbxvhNyLBQKSKVHEQooGOwP8IyREOJbrCzely1HTdOQTqdZRNwJrFa0kSnHhoMsEYBCKxJCHMRkf98G7FS0kSVHXYrFYpFFxGVjt8ybDDk2HGRzBqFSioQQhymUNBSL5u0oo8ybXTnWS3FoaIhFxLth5iasrNqnduTox1qChJA9H0VRMDZZ9ixfWZVjsxTD4bCpz+3JDDs+Pg7NwDUC2QXBrciRUiSEeEU4qCAYUJBKj3qWr8zK0a4UgR4VY6lUQjqd7ihHp7pkmJEjpUgI8RJFAQb7A57nK6NylCFFoEfFODAwgGKx2FaOTreOMiJHSpEQ4gcURUFyzqDn+aqbHGVJEehRMYZCIQwNDbWUo1v9FDvJkVIkhPgJv+SrdnKUKUWgR8UIAOFweJoc3W4y3OpgoxQJIX7EL/mqWY6ypQjMQDHm83ksWbIEiqLgf//3f229V70ct27d6qoUdeoPtp07d2Lnzp2UIiHEl/glX9XLcevWrVKlCMxAMX7hC1/AvHnzpL1fOBxGJBKp/hyLxaS9t1FUVUUikaj+nEgkKEVCiC/xS76qz9WRSESaFIEZJsaf//zneOihh3DttddKe89sNotcLodoNApFUbrOVnWCQqGAdDqNUCiEUCiEdDrNfo6EEF/ih3ylXz5VFAXRaBS5XE5qbdUZI8bt27dj5cqVuP322zFr1iwp71l/TzGZTLadkOMk9dfoh4aGqhUa2OyYEOI3/JCvmu8pJpNJ6YXHZ4QYhRA466yzcN5552HZsmWGX5fP55HJZBoeOq0m2rSakOMkrW5ce9HsmBBCuuGHfNVuoo3swuOeivGSSy6BoigdH6+88gquv/56ZLNZrFmzxtT7r127FrNnz64+FixYAACYnJxsO9HGLTl2ms1FORJCqkSHgUXnND6iw67ugh/yVbfZpzLlGLT1aptcfPHFOOusszpus3DhQjz66KP47W9/2zBJBgCWLVuG0047Dd///vdbvnbNmjVYvXp19edMJoMFCxZgcnISe++9d9vZp7ocU6kU0um09BlXRqY46wdbOp1GKpWSOuOKEEKMUixq2DU5ilAo7Fm+MrokQ8/p+tVBqysMPBXj8PAwhoe7f/P59re/ja985SvVn7ds2YITTzwR69evxxFHHNH2dZFIZJpMgcpspm5/MKfkaGbdD+VICPESTQMyuTIis7zLV2bXKcqQo6diNMrb3va2hp/7+/sBAIsWLcL8+fNNv5/RJRmy5WhlMWzLgy0UsLwPhBBilEJJQzBQKQlnOV/ZkKPVxft25TgjJt94iax7jnYqRPCeIyHECxRFwWB/wJN8ZbeijZ17jjNSjPvttx+EEFiyZIkrn2dXjjLKJjUcbOlRFCw0DyWEEDOEg5VJkGaxK0dZZd7q5Tg5OWn4dTNSjF5gVY4yawk2HGzZsqXO2oQQYhQLTqxiVY6ya5/qcqQYHcKsHJ0osKuqaqUFTKDSWVvTbBy5hBDiIGbl6ERBcKAiRzPlPilGkxiVo5NV51VVRTIRQDCgoFCmGAkh/sWoHJ2Sog7F6DDd5OhGKxb1rZvidi51EEKIG3STo9NSNAvFaJF2cnSzP5miKAgHeJ+REOJ/2snRb1IEKEZbNMsxn8+73rSTZ4yEkJlCsxzz+bzvpAhQjLbR5aifKQaDQTYZJoSQNuhyDAaD1TNHP0kRoBgJIYSQBihGm+hnivqZY6lU8qTZMSGEzAT0e4qlUql6pui3il4Uow2aJ9pEIhHXmx0L4fhHEEKIFJon2kQiEV+Wu6QYLdJu9qmbzY6FECiUGUJCiP9pN/vUj7WgmVUt0G1Jhhty1ITA2ESZZ4yEEN/TbUmG3+RIMZrE6DpFJ+WoaRrSmTJKZYFwgGYkhPgXo+sUnZYja6U6hNnF+07IUdM0pEfHUCwLDMYCUFWKkRDiT8wu3ndKjtlslmJ0AqsVbWTKseEgiwcQCjF8hBDnsHOrxmpFG9lyzGazyGQyrJUqG7tl3mTIseEgS85BmFIkhDhMoSQgLNjRbpk3WXLUpZhIJChGmciqfWpHjn6sJUgI2fMRb03y8yJf2ZVjvRTj8bi5zza19R6C0WvNsguCW5EjpUgI8YpwUEWpLJAeHfMkX1mVox0pAj0sxmw223Ebp7pkmJEjpUgI8RJVBQZjAU/zlVk52pUiAAQtvWqGE4vFkMlkAKDlH87p1lG6HFOpFNLpdMvPMH2QacXuH6yVAFGe/ly312pFKFpp+ucpPXn4EDIzsDzeG8UTCqkYis9BKj0mL1+ZRJdjOp1GKpVq+xkypAj0sBhVVW0pR7f6KXaSo6WDbPQP3bcp54Hdf218TisBgUjHlwW0IuYUtzY+tzsLqKHun0kI8QaL470V0vOVBbrJUZYUgR4VY7FYRCKRwKxZszA+Po5SqYRYLIZisYixsTEEg0EkEgmUSqXub2aT2bNnY2xsDDt27MDg4CAAYGxsDKVSqfrztEsHQgNKAsWSQElTsGvKYBjLAsjNanoyDAQ6C66oAZlSX8Nz0WIIIYqREP9icbzr7JoKolSu5BkUioCiWstXkkkkEg37EAqFMDk5iYmJCfT39yMSibTch2LRwFW1t+hJMe7cubM6AadUKlW/BQkhoCgKgsEgtm3b5tr+aJqGqakpTExMAKjMBNN7lbVECEAUoQkBISL42Wv7Ihgw0LFYaED5bY3PBWYBigohalOzVUVBWQiEAkrd+x5g/Rc0QaksUCwLBBQF2lvxCAcVVxsyaxpQKGlQ3vpQIQTCQRVuttjsHg93YDzw1mfOwHh0GO9GP2t8d+W91Ynt1edN5ysHEEKgVCphy5YtUBSlEpe3rgLqVwKb2b17t+H370kxKorSUPQ7l8tVbyqHw2HXmwyrqgpFUarfcsLhMAKBQIdXCKBc/T9EQipCQSNiFECx6Z5DSKkOlGhYYMd4CYBAJKSgv6/TPjjH+EQJU4VKPOYOBKoJ0U2CRWBsovK3GuwPeLJulPGowXjUMByPLuO9G8WSgKoqDfkSsJKvnEFRFOTz+eoJTbdLuGaO254XY/Pl0nK5jGDQ3T+LpmkN+1EqlRAIBDoLWu0DNA1aWUOof2+EwwYuj2hFYKrpPftGqvcKs9kJqKHKbN2SogDROa7PhC0UCihNjkINVRYVF9Q44vF+V/dB0zTsGh2DGqoM/CkRxqz+Qde/MDEeFRiPGqbi0WW8d/+wIpSJ7VBVteFvbSlfOUDzpVFN0zrmborRIKVSCYVCAYqiIBqNolgsVgMeiZi/QW0FTdOq33r0zywUCsjn84hEIh0ONqXyUN76BqgY+MamlKd/W3zrtdlsFpnsBBKJ2YjFYpUb3OkxV5eJFAoFpNJjCIXCSCaTmJycrFwWUVTbN9ONUqlFO45isYTh4bkAULnUPjru6GSsZhiPCoxHDdPx6DDeDaGUpz1lPV/JJZ/Po1QqIRgMIhQKIZ/PV89gZZzY9OQ6RqBRinpAI5EIgsEgSqUS8vm84/tQf5CFw2EEg0EEg0GEw2EIIZDP511pdtw8m8uLFjCtZgPH43EkEglkMpmu605l0Gp2nZv9NXUYjwqMRw2/xMMP+apeivW5W7+8K2PSZE+KUdO0aVLUcUuOrQ4yHTcPtmx2suUUZzcHf6clMm4N/k5Tzt1Mxu2mnDMejIfX8fBDvmqWoo5sOfakGIvFYksp6jgtx04HmY4bB1t2d/mty0Ot1/24MfiNrBt1evAbWYflRjLutg6L8ajBeNRwOh56DvI6X7WToo5MOfakGAF0vR7ulByNSFHHyYMtu7uMzJSGRLy/4/0iJwe/mWIKTg1+M4uTnUzGRhcnMx41GI8aTsbDD/mqmxR1ZMmxJ8VodEmGbDmaOch0nDjYqlLsUxGPd2/F4sTgt1JhSPbgt1Kxw4lkbLZiB+NRg/Go4UQ8yuWy5/nKqBR1ZMixJ8VoZtquLDlakaKOzIOtQYqzjK89kjn47ZTdkzX47ZSxkpmMrZaxYjxqMB41GuNhvGN9M5qmYXx83PN8ZVaKOnbl2JNiNItdOdqRoo6Mgy2bnbQkRR0Zg19GLVq7yVhGbUcZydh2axzGowrjUaMaj+wEsrunL7nohh4PXUhWFu/LyFdWpajTLEcz+0AxGsSqHGVIUaf5YDPTWbu6DsuiFHXsDH6ZBdqtJmOZBY/tJGNZBY8ZjxqMR414PI5EvB+ZKc2UHOvjMThor4iCHTnalaJOvRzN1EqlGE1gVo4ypahj5WCrDvp4vy0p6lgZ/E50LTGbjJ3oAmAlGcvsAgAwHvUwHjXi8RgSfaphOTbHIxSy3yTASr6SJUUdXY6mXmP7U3sMo3J0Qoo69QdbuVzueLA1DvruE22MYmbwO9nKy2gydrI1jplkLDsJ6zAeNRiPGvFZAUNydDIeZuQoW4o6qqqa+p0oRgt0k6OTUtSpP9jGx8dbHmxODXodI4Pfjf6W3ZKxG/3ijCRjxqMC41HDlXjUy7HFhBw34mFEjk5JUcfMpEuK0SLt5OiGFHUCgUB1H5oHv9ODXqfT4Her6TPQPhm71UQV6JyMGY8KjIdH8dDlmJ3wLB6d5Oi0FM1CMdqgWY5uSlFHVVUMDg42DH63Bn39PjQPfjcHvU5zMnZz0Ou0SsaMB+Phi3jMClQm5HgYj1Zy9JsUgR7vriEDPZClUgmlUqnaF8zN1lWhUAhDQ0NIpVLYunUrALg26HX0wZ9Op7Fz504AlaTkZgcEANXfWW9YqiiKqx0QgFoyZjwYDx3/xCMGqEFP46HnxkKhgKmpqepzfpEiwDNGKdTP3lIUxfV+jkBlkNUfWLGYvIk2RlFVFYlEovpzIpFwvUcb0Pi7RyIR13vmAYxHPYxHBcajRjAYbLjnJ2MGrEwoRpvolwL05sf6z26TzWaRy+UQjUahKIprLXnqKRQKSKfTCIVCCIVCSKfTrrTkqUe/PKT32Mzlcq605GmG8ajAeNRgPGrol1FVVYWiKK61rDIKxWiD5nuKfX19rvZz1NEbliYSCSSTSdf71QGNEwmGhoaq66Dc6lcHTJ9IkEwmXe1Xp1N/D4vxYDwAn8RDiLeaPnsbj/p7in19fa73czQCxWiRdhNt3G52XCqVMDHR2DrK7WaurSYSuN3Mtd1EArebubaa2MF4MB6ex0MIpDNlFEslT+PRaqKNF82Ou0ExWqDb7FO35KjX/+vvn946yq3B32l2nVuDv9vsOrcGf6fZjoxHDcajgmvx0KVYFhiaM+hZPDrNPvWbHGeUGO+//34cccQR6Ovrw+DgIE455RTX98Hokgyn5ZjP51Eul6GqatuJBE4PfiNTzp0e/EannDs9+I0sAWA8ajAeFRyPR70U4wGEw60nuTgdDyNLMpyW4x5ZRPwnP/kJzjjjDJx99tn44x//iCeffBKf/OQnXd0Hs+sUnZKjfpDpC/w70TD4R8ehmSg83gkz67CcGvxm12E5NfjNrItzKhkzHjUYjwqapjVKMdQ53TsVDzPrFJ2SY6lU2vOKiJdKJVx44YW45pprcN555+Gd73wnDjzwQHzsYx9zbR+sLt6XLcf6g8zoNOvq4C+VkM6UbcvRyuJk2YPf6uJk2YPfymJx2cmY8ajBeFTQNK0y0cagFHXq4zE5ab2fo46Vxfuy5VgqlUz/PWeEGJ999lls3rwZqqri0EMPxcjICN73vvfhhRdecOXz7Va0kSVHOxUiwuEwhuYMolgWtuRop2KHrMFvt2KHrGRsp4KKrGTMeNRgPCpU41EqmZKijh6PiYkJ0w1+67GTr2TJUZeiXnjFKDNCjK+//joA4IorrsCXvvQl3HfffRgcHMR73/tejI6Otn1dPp+vVnjQH2aRVebNrhxllE0Kh0MYigcsy1FGGSu7g19WGSu7yVhGWTG7yZjxqMF4VGiIx5xB01LUicfj6O/vh6ZplgQtI1/ZlWO9FPWejEbxVIyXXHIJFEXp+HjllVeqf5BLL70UH/nIR7B06VKsW7cOiqLg7rvvbvv+a9euxezZs6uPBQsWADB+E1Z27VOrcpRZSzAcUhvlaPBvIbO2o9XBL7u2o9VkLLPWptVkzHjUYDwqTI+HvWoysVgMqqqiXC57lq+syrFZimbj4akYL774Yrz88ssdHwsXLsTIyAgA4MADD6y+NhKJYOHChfjrX//a9v3XrFmDXbt2VR8bN24EABSLxa6XCJwqCG5Wjk4U2G2Q42jrllX1OFHw2Ozgd6rgsdlk7EQBarPJmPGowXhUcCoewWAQgUDA03xlVo52pQh4XER8eHgYw8PDXbdbunQpIpEINmzYgKOPPhpARW5vvvkm9t1337avi0QibYOjH2ithOd0l4z6wuP1PzfjZNX5ihyB1O5Ky6p2A9rJLgD1hZVTqVTbAe10F4D6Qtf1PzfjZFeG+kLXjAfjAfgnHuFwuOFkwot8VV94PJ/PtxWeDCkCM+QeYyKRwHnnnYfLL78cDz30EDZs2IDzzz8fAHDqqaeafr9wOAxFUVAoFKadObrVOqrbmaMbrVjCIbUyIafNN2M3WuN0+2bsVmucbmcqbrQq6namwnjUYDwquBUPP+SrbmeOsqQIzBAxAsA111yDf/qnf8IZZ5yB5cuX4y9/+QseffRRDA4Omn6v+pux9XJ0u59iu4PNzf5k4XCo5eB3s19cu8Hvdr+4dsnYzf597ZIx48F4eB0PP+SrdnKUKUVgBvVjDIVCuPbaa3HttddKeT9VVRGJRJDP56sHWrFYdLXJMDD9sqr+/272J2u+bJRIJKpdANzqF9d82SiZTCKTybjaRBWYfhlP/383+/cxHjUYjwp+iYcf8lXzZdVQKCRVisAMOmN0Al2OQK0NittNhoHGb2JedbLWB7/+TTgYDLreRFUf/MFgsPrN2O0mqkDjmYrbSViH8ajBeFTwSzz8kK/0M8f6Nn+ypAj0uBgJIYSQZnpajK2+bbSakOM09dfovejnCNTumejfjEulkuvNXPV7JqW61jhu9qvTqb+H5UW/OoDxqIfxqOCXePghX+n3FFtd9ZNBz4qx1USbVhNynKb5xrXb/RyB6RMJIpGI681cmycSRCIRV/vV6TRP7HC7Xx3AeNTDeFTwSzz8kK+aJ9o4UXi8J8Wo/wGbJ9ro3z7ckmO72VxuHmyFQrHl7Do3m7m2m13ndjPXdrMd3UzG7WY7Mh6Mh9fx8EO+ajf7VLYce1KMhUKh7exTt+TYbYqzGwdboaghNTrWdnadG4O/25RztwZ/tyUAbiTjbksAGI8ajEcFt+Lhh3zVbUmGTDn2pBgBdJx96rQcja77cfJgKxQ1pLJlhLrMrnNy8Btdh+X04De6Ls7JZGx0XRzjUYPxqOB0PPQc6GW+MrpOUZYce1KMoVCo65IMp+RodjGsEwdbVYoBBck5A12nODsx+M0uTnZq8JtdLO5EMja7WJzxqMF4VHAqHqVSCeVy2dN8ZXbxvgw59qQYja51kS1HqxUiZB5sDVJMBAz/LWQOfqsVO2QPfqsVVGQmY6sVVBiPGoxHhenxMN6xvhWTk5PQNA2BQMCzfGW1oo1dOfakGM0gS452yybJONgKhWKjFE30JwPkDH67ZaxkJWO7ZcVkJGO7ZcUYjxqMR4WGeIyOoVC0JuhsNouJiQmoqmqpgICMfGW3zFuzHIWJ/rMUowHsylFWLUE7B1uhUKhMtLEoRR07g19WbUe7yVhWrU07yVhWrU3GowbjUaEaj2AQqWzZtBz1ePT399uqAmYnX8mqfVovRzPHJcVoEKtylF1gt/5gMxro6qAPBm1JUcfK4Jdd8NhqMpZdgNpKMpZdgJrxqMF4VFBVFck5AwgFFFNyrI9HLBaztQ+ANTlKLwj+lhzNQDGawKwcnao6rx9s5XK56z40DPo5A7alqGNm8DvVBcBsMnaqK4OZZOxUVwbGowbjUUFVVSQTAcNydCoeZuQoW4o6wWAQoVDI8PYUo0mMytHpViyRSASBQACapmFycrLlNk63xjEy+J1ujWM0GTvdqshIMmY8ajAeFRyPh6I0yrHNhByn42FEjk5JUcfM+1GMFugmR7f6k4XDYaiqiomJiWmD361+cZ0Gv1v94rolY7f693VKxoxHDcajgmvxqJfj6Jhn8egkR6elaBaK0SLt5Ohm006gcomgv7+/YfC72UQVaD343W6i2i4Zu9nUFmidjBkPxsPzeOhyrGtZBbgfj1Zy9JsUgRnUqNiPNDc7LhaL0DTN9f5ksVgMwWAQmUymoXmnm/3i6pu5plIpAGjoAuAGzc1cI5EIcrmc6/376pvrMh6MB+CTeCiVgh7psayn8ahvdqxpGoQQvpIi0GNi1NexTE1NSX9f/RuYoigQQqBYtLe41ujnApUBHwqFUC6XMTo6CgAYHh6efm9FFIFc0+9ezACK8ZvS3QiFQhgfHwcADAwMYGpqSvrf2+g+TE5OIhKJQNM07Nq1y9V9ANA9Hi7AeNTouXi0HO8ThuNRLBar+6ZImrRXj16zGqh8cSiXy9I/ox79dzGynlERZlY9znA2bdqEBQsWeL0bhBBCPGLjxo2YP39+x216SoyapmHLli2Ix+OOfAPKZDJYsGABNm7ciEQiIf39vYa/38yGv9/Mhr+fPYQQyGazmDdvXvd6q9I/3ceoqtr1m4IM9E7jeyr8/WY2/P1mNvz9rDN79mxD2/njTichhBDiEyhGQgghpA6KUSKRSASXX365q0s13IS/38yGv9/Mhr+fe/TU5BtCCCGkGzxjJIQQQuqgGAkhhJA6KEZCCCGkDorRQe6//34cccQR6Ovrw+DgIE455RSvd0k6+XweS5YsgaIo+N///V+vd0cKb775Js455xzsv//+6Ovrw6JFi3D55Zeb7kzvJ2688Ubst99+iEajOOKII/D73//e612Sxtq1a7F8+XLE43HMnTsXp5xyCjZs2OD1bjnGV7/6VSiKgosuusjrXZHG5s2bcfrppyOZTKKvrw8HH3ww/vCHP3i2PxSjQ/zkJz/BGWecgbPPPht//OMf8eSTT+KTn/yk17slnS984QuYN2+e17shlVdeeQWapuE73/kOXnzxRXzrW9/CLbfcgv/7f/+v17tmifXr12P16tW4/PLL8eyzz+Jv/uZvcOKJJ2LHjh1e75oUfvnLX2LVqlV46qmn8PDDD6NYLOKEE05o26d0JvP000/jO9/5Dg455BCvd0UaY2NjOOqooxAKhfDzn/8cL730Er7xjW9gcHDQu50SRDrFYlHss88+4j//8z+93hVH+dnPfiYWL14sXnzxRQFAPPfcc17vkmN8/etfF/vvv7/Xu2GJww8/XKxatar6c7lcFvPmzRNr1671cK+cY8eOHQKA+OUvf+n1rkglm82Kd7zjHeLhhx8Wxx57rLjwwgu93iUpfPGLXxRHH32017vRAM8YHeDZZ5/F5s2boaoqDj30UIyMjOB973sfXnjhBa93TRrbt2/HypUrcfvtt2PWrFle747j7Nq1C3PmzPF6N0xTKBTwzDPPYMWKFdXnVFXFihUr8Nvf/tbDPXMOvVPETIxXJ1atWoWTTz65IZZ7Aj/96U+xbNkynHrqqZg7dy4OPfRQfPe73/V0nyhGB3j99dcBAFdccQW+9KUv4b777sPg4CDe+973VtvezGSEEDjrrLNw3nnnYdmyZV7vjuO89tpruP766/HpT3/a610xTSqVQrlcxl577dXw/F577YVt27Z5tFfOoWkaLrroIhx11FE46KCDvN4dadx111149tlnsXbtWq93RTqvv/46br75ZrzjHe/Agw8+iPPPPx+f+9zn8P3vf9+zfaIYTXDJJZdAUZSOD/3+FABceuml+MhHPoKlS5di3bp1UBQFd999t8e/RXuM/n7XX389stks1qxZ4/Uum8Lo71fP5s2bcdJJJ+HUU0/FypUrPdpzYpRVq1bhhRdewF133eX1rkhj48aNuPDCC3HHHXcgGo16vTvS0TQNhx12GK6++moceuihOPfcc7Fy5Urccsstnu1TT3XXsMvFF1+Ms846q+M2CxcuxNatWwEABx54YPX5SCSChQsX4q9//auTu2gLo7/fo48+it/+9rfTSjctW7YMp512mqff9Dph9PfT2bJlC4477jgceeSRuPXWWx3eO2cYGhpCIBDA9u3bG57fvn079t57b4/2yhkuuOAC3HfffXjiiSdc6aLjFs888wx27NiBww47rPpcuVzGE088gRtuuAH5fB6BQMDDPbTHyMhIQ64EgAMOOAA/+clPPNojitEUw8PDGB4e7rrd0qVLEYlEsGHDBhx99NEAKt2w33zzTey7775O76ZljP5+3/72t/GVr3yl+vOWLVtw4oknYv369TjiiCOc3EVbGP39gMqZ4nHHHVc92+/Wv82vhMNhLF26FI888kh1uZCmaXjkkUdwwQUXeLtzkhBC4LOf/SzuuecePP7449h///293iWpHH/88Xj++ecbnjv77LOxePFifPGLX5zRUgSAo446atrymldffdXbXOn17J89lQsvvFDss88+4sEHHxSvvPKKOOecc8TcuXPF6Oio17smnTfeeGOPmpW6adMm8fa3v10cf/zxYtOmTWLr1q3Vx0zkrrvuEpFIRNx2223ipZdeEueee64YGBgQ27Zt83rXpHD++eeL2bNni8cff7whVrt37/Z61xxjT5qV+vvf/14Eg0Hx7//+7+JPf/qTuOOOO8SsWbPED37wA8/2iWJ0iEKhIC6++GIxd+5cEY/HxYoVK8QLL7zg9W45wp4mxnXr1gkALR8zleuvv1687W1vE+FwWBx++OHiqaee8nqXpNEuVuvWrfN61xxjTxKjEEL8z//8jzjooINEJBIRixcvFrfeequn+8PuGoQQQkgdM/PGCSGEEOIQFCMhhBBSB8VICCGE1EExEkIIIXVQjIQQQkgdFCMhhBBSB8VICCGE1EExEkIIIXVQjIQQQkgdFCMhhBBSB8VICCGE1EExErIHs3PnTuy99964+uqrq8/95je/QTgcxiOPPOLhnhHiX1hEnJA9nJ/97Gc45ZRT8Jvf/Abvete7sGTJEnzoQx/CN7/5Ta93jRBfQjES0gOsWrUKv/jFL7Bs2TI8//zzePrppxGJRLzeLUJ8CcVISA8wNTWFgw46CBs3bsQzzzyDgw8+2OtdIsS38B4jIT3An//8Z2zZsgWapuHNN9/0encI8TU8YyRkD6dQKODwww/HkiVL8K53vQvXXXcdnn/+ecydO9frXSPEl1CMhOzh/Nu//Rt+/OMf449//CP6+/tx7LHHYvbs2bjvvvu83jVCfAkvpRKyB/P444/juuuuw+23345EIgFVVXH77bfjV7/6FW6++Wavd48QX8IzRkIIIaQOnjESQgghdVCMhBBCSB0UIyGEEFIHxUgIIYTUQTESQgghdVCMhBBCSB0UIyGEEFIHxUgIIYTUQTESQgghdVCMhBBCSB0UIyGEEFIHxUgIIYTU8f8DAEndgipBy5AAAAAASUVORK5CYII=\n",
"text/plain": [
"