Skip to content

Instantly share code, notes, and snippets.

@tjlane
Created July 5, 2021 10:13
Show Gist options
  • Save tjlane/8062d1666caa31aba23a4a77b43a7f26 to your computer and use it in GitHub Desktop.
Save tjlane/8062d1666caa31aba23a4a77b43a7f26 to your computer and use it in GitHub Desktop.
soot sim
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 15,
"id": "237da705",
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"\n",
"from thor import xray\n",
"import mdtraj as md\n",
"\n",
"from matplotlib import pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "a14308fe",
"metadata": {},
"outputs": [],
"source": [
"s = md.load('Cluster_02_new_log.pdb')"
]
},
{
"cell_type": "code",
"execution_count": 42,
"id": "c706ecf5",
"metadata": {},
"outputs": [],
"source": [
"# this is a \"generic\" perfect square detector\n",
"# if you want to use a CSPAD, we can -- just send me a geometry in crystFEL or psana format\n",
"\n",
"spacing = 1.0 \n",
"lim = 100.0 # size of the detector (in \"units\")\n",
"energy = 10.0 # keV\n",
"l = 50.0 # z-distance (units)\n",
"dtc = xray.Detector.generic( l=l)"
]
},
{
"cell_type": "code",
"execution_count": 32,
"id": "532677a4",
"metadata": {},
"outputs": [],
"source": [
"num_molecules = 100\n",
"num_shots = 1"
]
},
{
"cell_type": "code",
"execution_count": 33,
"id": "1f8d4e8f",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU Thread 0 :: 100 shots\n",
"12:06:25 - \u001b[80D\u001b[1A\u001b[KFinished simulating shot 1/1\n"
]
}
],
"source": [
"shotset = xray.Shotset.simulate(s, dtc, num_molecules, num_shots)"
]
},
{
"cell_type": "code",
"execution_count": 38,
"id": "f2b3c330",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAQoAAAD8CAYAAACPd+p5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAADIgUlEQVR4nOz9ebB1S3YXBv5W5t7nnHvvN725BlWpSqNBAmSLkAcGq02bwI0tWrjB4A5QAAbjBhx06A+Gdthu03RgG3C4wzY2BG7jxoBp0xibwNiyadm4QSCEsEEqzVJJVfXqvVfvfdMdzjl771z9x1orc2Xufc693/e9J24oXkbcuPeeYe/cOfxyrd+aiJnxYfuwfdg+bMda+PvdgQ/bh+3Ddvvbh0DxYfuwfdiubR8CxYftw/Zhu7Z9CBQftg/bh+3a9iFQfNg+bB+2a9uHQPFh+7B92K5tHxhQENEvI6IfJKIfIaLf/UHd58P2YfuwffCNPgg/CiKKAH4IwD8J4HMAvhvAr2Xm73/fb/Zh+7B92D7w9kFJFN8E4EeY+ceYeQ/gzwD4FR/QvT5sH7YP2wfcug/ouh8H8FPu/88B+IcPfXhFa97gbP4GEUh/I5D8tp8QwKTvAYD+gglIzCAGkBjgpL/lR6Qo+6C7B45fz36zvdEKY1T/UV3XWv5+9eL82W/wlrvJNR+67sIE8s9PNHuW9s/DL5ZxRtKx9lIr6b1sPu37Nic6TwfHWG9JsGuEps/6282X9SXPv++5fc/6MrtO+1yuX3qPIpUf6G81TLTw0vPM36F7HJg7AE/Gd77EzK89z6U/KKBYevJqCInotwD4LQCwwSn+4fhLZUPboOvgUYxAjKCuA61WoPUKWK/Aqx68XoH7CO4DuAtlwJlBYwINCWE/Ars9aLsHb7fg3R4YBvA4us4EUAxAjEDQv8kJW5zAUwKmCUj6NyewgZC7jvzSBRzcBvQDYQtXr32scfu5ZnxAodyvuR4nlvdc3/xn8vvu+anvgK6Tv0MAghvX9llC8zozMCXwMAD7Abzbgfd78DTJe0Sg1QphvQbWa9Cql+8yy3zsB/n8OGaAz89v/bbnjTGvB+p7YNUDMYJjyH2hcQLGqfRnHPMcSv/1mWMA+hWoi/LsnV7H1uCU5DvjBIyjXGcc85rgaVruL2zo3Rzo36SHXR67EBbXijyKXjMdWCv2XR0XAbv5+vvLb/77n12+wPXtgwKKzwH4hPv/ywB8wX+Amf8ogD8KAPfCK0yBwCkASOUUAMDTVAO7fDkjEXGPhA7EQOpCgahAQCRwF0CpA8YJ1HWykKcJoCC/AQDyu0L4IIsySyOl34efWheybMAEpAAOdkm6/vvP0nTjyW31frkbBcB4gm6uepHZgrbvMhPIn7xtW1rEiWWc80J23/UAmMFtYTMoUGTJb+n7+bONppySzGdMsplNArV7ttKMb7pJiXRz2SaLQUDCXksigRAARAYmk2J0rd60HQIJv8mXWmIQ0Xzd2Pz772fpyoGEB/EXaB8UUHw3gK8mok8D+DyAXwPgnz/0YZmECCKWhQ13ejIXxNbFbcBhvwMA5ii/I4FJfhAI6AIwRaDv9IQZZRFME4gpbyqe0Fw7ClgcaJU04U95BQv9EJBk0t4XgFi6hr/foX4hKQjL6dZKKZwCiBiI7rqJIbrbNa068Vk2r72Wlk9YBK96LFzLfafqr/Y595f1WaepbDzbQAY+qQZ6NJtUAEGlpxjAXcyAIZuuqBeU9DuYMGvHpAk3PwWYmg1u4+Kfn1MGq/peS0AetM8mEVtfcVgSeYb2gQAFM49E9NsB/LeQ5fcfM/P3HfyCIjsH3aBMeuqTOxV0M+/1K/Y9FNCQm4eyEAMBSSQLdAoWUw+YyJgYedIPgcX1D7v8son1BhbA4VPjeZuTKqpmklDVxyLtaAfbDsvGZgYlvhnNvXTK2eZdAgjXv0qySKn059iizgCYAKasBlJkcEogOwT8Ker6SEQimdpJHN1a0c0LA4wQZAySSimH+tUC0aEWaA4SJrXMVDq7dliYp/peGWyculuBBCB/Twvg9gztg5IowMx/CcBfuvEXQgClJKJ6CqBo4nE6CBaA8lDQjc0McBQpggvZyXZydMJ1oJtAqlsSk1zTXx8ABVVN3kdUvtGCcm3xNG6b6ez57wOnuFyw4SqKxAFEBYklkb8RjSsVw91LgcYI41aKotAs4Ob7h6SumVQBHZspgcIETHKC5gOjnTO/2fQ9MhGdqEgT0YFEZzojCyjNSOl07ZqYcVUtB5YJ+tBIDk4FM0I4pRlI5Oc0lSPQnFeSgT3az5u0DwwonrnpSUDMAMmGlUXf8hYHwEIH1IaHO9nsAAShTaroIqiL4FFAo5YqDJhEoqEYZYLaxX2seULzg25+MeTFFubjBtTSx4L+z4lBrQBFCxuk+lJDsmVx/xkX5jEu4dA90UgVEHI6HxhGkvrrBgKl8nfhJZw0oSDB0YjM5v4m/aQCiIfUDmuZk/Ag4cAi81e2zJYOfz/WFQlapIkZQWrfMwB/gXY7gCITkCEPBAFCsGE6DBYjVZeo1BBmcIQABJDFMTapYpVA0wTECeBYmHnlRPK1DCye5XEavfSaDx/cuO6f4/eKsYjykWTcPN9zSEXJ12dUdqr2fl48biQIey1bcg6pEEv3P7Z4D4xL9XUnVTAAUhTgGKsTedYc8WeEX5YmlOPKaq/vq/vJp/wBgPDSREWcGkjYnFG9hqsNzancB6h/5+8XaeLgHHsr23O22wEUgENVBxYpgVuS04PFMNYXUYlCThV9NA6F3IykJrAISh141ZfrJUbWO5jlhCUGQyUL4Dhg6KRVC8S3Z5iom6kctkicvmtvqVQGoCaH8wcObMIli8exxQeA/WYxq0UGi8PPUUl/nni8aWulv7avh4Da5sXmyYhLL1FEBxT6Q5NaWLyZvJUmWtUOyGBgZv4MEnYv37x6kdyYOKlgZjXy0sTSGE2T9DndUo7ieZsteja7O+tmRTwMFgvWEHuNI0BZ33RSRUqgScyldjJwY57jSQhNxnTQxj1/gGX/iSxWXgMCM5C4TpoIRszFijXPmwcF/665MdCSt8GNmz9hvc7cmjVtgd9E1L0OGBTQjgGnAbqoq+WZq0eZgZ+e8N46QM5ahua7tkamCUiTAws+TA57aSLo3y1IUANI/p5c+B6vOjBzObjaw8jmy6kpAsLT7SUzn6u5hU5GQqYkA2PKpVlEDllDOGWwIADoGYwuqyBMIlVkdp8Z5BY6j0N1bdGDiz/ETdpMV1SSkJmXzV3P2pakCTut7H4xgscRhEFfK05WB9sSKB0SaZmz4xkmRxqmOYm5eB9vWXne8WgJaIhESpxEz283kjdDmsgeAzjzFVT6xAyaxHEPkzpujRMwiNOWSRP5earbKL+RpYmQ54iMWPcEcUIt2XppIm92Ln0HahDP10lFfTf1SJ3MZtL3M7bbARTtOrGBVpfaLF1ATkeKqMGCRUfFHo36UcRbcBSCKkCO9xBEBeFeTiXfnWmqNpY5T5UPXHMSeju9NSNqs2XnBptjyWmpug8Vz0Kv9wJK8JmENThrTvGdqMhOazcR/x2JBwMLDxDpiEUgHXimyjJBFd9RWT1aos6rIExizYKsEfNhmbVsUnSnOrlNG1h+j0kc9QbxxBRvTPXENPP6obXgJT0K4qthoO45hdbqw6mSJiqQmD3HggXJQCwVSYIV3F6k3Q6ggIpUCQvilDDVbGYtNVsKaExOzJqa6AVkACgbJmbxkrpYSR5ZZbH+TCiTluaejQsPUF/AGO7cIcrmVroJWBxagI4LyZKLc3PP5r+YwIFqsGBets237Vi/vBOT+aN4wk3Hojpxl67BCfn4vq4/MLBowNptNK+C5I3fAjxQDiEPEE4FIGaRRlICDaOAxH4AD84FvAWJliw2Sc8AIboQAZViqsb1mC65htt18y0qPxSuVVsDcA9q488EicI3cwv2aBlFpTDpwm/sPGGqi83AAg0YRHeKxADmWN4PBDKRdCTnlGWSxTU8hdPzMyOdfRy4YrUrsLgBwz8T/ykUxx1y9vNO9VcO2cEoP/s0gTEVE6o/ZJTAzeRd9hNoeAtmJTCnshBboLiu7/lSBwBpCago1OZbU0G1T5VUwWYqJywRjPn/mRmRACTQxEWSMJDYDzlGqFpzh5pJE/aTTaROenGcjx/TLJEdk1gAZMczqKrlOCNugeJnikQhIrljb9sJpgCEJNJFjA4EGvPpdWDBYvXIE9cpb+FMZQgEGkTHxzDCCLUbWSO8hyNQFoXZ5ZVDAXBYHDy2OMyyEkpfPXufg5nYSVOJK+7GwELDa92CTTXTDqBy5XYkXrUIl3ToBZKvBiIuh8Kx1rg357iHtv/OUlVLFYfAYq7CUErAJKc6jZMGqUlAGYYBPKVrQWI+N008iQcJT5QqL8HeUuGI3PaQykCocUs21uyI13wNL3k/Z7sdQOHZ3RCWVRBANwSaRb8AFlyflkVaKKJ4BRbGdgchnihGcBjU1ElCChJXE7e4md2pXPXeJjICmCBmVypu6zcm85z5tRCm7pRqvPJYn5emCOYOhFFcnfW+2Su1eoZUdN2U1Nu1nHwZjKd6MeamQWjXWSsymWngYZKbX9AOJDxBnM2/3EpFdsKGfK2Da6l6ZnaErEadDqNGGks0K1sw4TFJwhPMNjc5mnOBWzGnLS9NOKL0+PjpHiAC85hfE5J5ypavayWfG7ZbARQMVIxtAYuE9kQxUZ0ijoOFY8MF1csGIkD9Kpwa0gnicwgg/RGCcAD2QU8UHHeRBsqpPE2qi7r+p1AkeXeKmdt6pYteE76eidLWFNvowNxF0NQVcTxxUUNIT+XGfwTTJBG33oUdKMBwzYK+ls/xrtWLLuWuLcRIZDBp+58voxJL9aLea0mSmFIBaw0nzyCxH8DD+GybLsdyONCoO1gkt8kBRFY7GrOrj0iG8xaenPTsrSMmRRwa0+dotwIoAEP/I2BhYqox4iQRhEfVEO/BGZx/g10HkNwDRLKfScxX5qVH0ZFSRMJ6D+Ncv589jKofGrBUhb6nAAQGxXBcBWkdd9q4AWoAovqwqTtBxqHTxTUlCZXm8kOq17fqB0+TSCIxFX+MlEoeBhNvF/RobgnEluyz+yAUYD3WGvfnrL5N08zdedEV/VibUn0IpaR8xLODxEztaFtiAA6QMvCmZQmNFwDViMt2zRhAHFKNiIAXwIxbAxS5PQtYmEXkGrDANAHDkFUMspOJSAgsEw87jS5kkSiydBEjmELZMMOoJ9k16oedGGQAYcQmFaDzKsiEWQBUm/TEn67H3HbNDVkiKwOoiyJFcACxgBSlNJcqAHmufNIpoDlzmy3oxUVJdJxjacaq/HlEmvC+CEQAk/jA6Iac3S45niIeQA2vUk3GF+hGHfYZIJ5VfK/UjkPSkr+Xt6Rws9mr/jrrm1Pv5C0H8kv9NGL9BdotAQrOSElEzwAWjuC0DRdjDRbQgdRJoUCy6e3zeaO60GKoCtKp+qELllGTSmJCXc5NYAAhBG0SlcOftEEWPDjkhWOTXxFXHiDaZCdLk28bwACkCnbqnMrDOSaikipMxyURbWkkp0+nZvMsLMysFjTq0NHpPwIsmXcJbsGnWv0ics/ddkdJv6XWxqpMU5YaWX/fGCRaYFi6j7+X93NwKod8pvQrP5++XqS1afmzs37poRTjcrDZDdstAQqoCOoyQt0ELACYNcS+aLp3pR4Y6RUIGISvyORmJjIDKDBY30MEwEWqgCNCKVsGUi22V4/TSBWBAKSyoBqAKcRmw1Us6OhV7oF8Q6dSWAsAJzmNOQVQCqApgEPMKog8/wJXMU056I6MaDR7vDeLHpzPRspADbKzxDjeSav1bGyfuVVtjrSSC7QeryqXp21akyLslH8eIrCx3M3Wh/dzaB2rKrK8kUDkYWpQPdY3Awgjv1vfjWdstwMoGIWscRmhrgUL17LLt1ocJAWcE+H8Ka+SRdYl1drBbKcwgYOqIV1AMI9ObzXxzHXr4ehVHpMqpknF4JqkNPWjSi7Ssv0xllgOb+Xwn2UDryDmPSJwcqpODngyMswAZEGqcOOVnxWoPRKftR0TfbNpNWRXdx+a7fX+DKgvKEoX8tKZe9Xz8kZA6NshicmklPY1u19WQaabWb5uJNkUSdPWTBs0+DztdgCFttZdmsMRyQKoF4vxFbbAbPG7E7pKdjJQ5T9hPhyifnCJAQCQAiE4+zcp+meGelwWnWdh0MDyhClI5E2grQKJnJ5tvkEyqKYAmlIGzAImjpCNTqpIGmBkvMqCVEH2HPJAh3mJ0pm6c5UK0kQ5HlA55nk1Q/ku0c03sG9h4b6eK2idkw7p+gfunYlM6HywrWN/gKTiK+GiQu29oy7h1zUPEF6KWAoneI52e4BiZgZCVkVasMBkJ87CoFogmYrUB1UQItBonxkLXxFJczpwkSyYwCZ86KXIbO8mMi4tLrXrV/ktgBnAYeI8mVWGqaXIQ7dZ6twFTvVICZiWIlhJ8j7GqObSIo3YfesIWuMrDpBlrVQjkzbXrU0EXmqpeY42PsZbEI6d3P594y+yNBLm38/OSDcAies4ltkzmY+DSmKV6uHGpw34etbWju8xPus6X5Jr2nMDBRF9AsB/CuAjEHn6jzLzv0tE/zqA3wzgHf3o79W0eNe3fPLrg3veQk9cOymzdAEUK4hvdqqaB1/lFOSJxlE2kLN9IyRQF8BTAQv0AQnIzoA0aeKbSnSsNxlILBg5rwVQFlDrwej77bIY+ejQHHmo16eUxGQJPcX0uSqLjt889r9dU/tK7FSQGOvNwk7K8O0ACHgS2d7Pvx0BmZ2hlq7pxyK4/vtmJuiZD0f5XOVG30o9HwRImLWl9XHI7zdmZEcSz9p19/XjfhOA+Pto9RgBfDsz/20iugvge4joO/S9f4eZ/+AzXa0hbErW6OOqCID5goNdisBg5IAiF+SVVRAiIe1sUAeRLijaJmXJehQI6ILwqCmB00qJSgUMU0EcmJVn02dwGZhmfXbhwZW42LoBWz/NtZc13Fv+UXWIQBMVpzKgAJJ9nwUcjYOw3prUcxMdPUdG+teCZ+bttcIv+HZQkrDvLW2WG5zAVUCWgYSZE7PFgWuQeJFT3Y1Btp7cpB2IaTl4qwYY5PY1OJTXXL9eUJoAXgAomPlNAG/q30+J6DOQCmEv1myTVdLFAVXESQYHW2bJZ545sEzfYjKtrSBEwkskVTdSDACJ+M59FFDoOyliM66UJY/5mvWtOKtAeUMe0hvddxdNoJnYQ3Vwm1pFdg2TKPJ4cg3Gxn8A4E6WgeXnqMDC5qTpXwmfbhZhztnpguhMZTqiL2dQaKVD703pidZWZfGfBxZMqu1YNV6Qh5yUbthm6QhaKeEGJ/rMLG7NgwJwGBgOfO+ZVaeF9r5wFET0KQD/IIC/AeAXAPjtRPTrAfwtiNTx8JkueI10kcECKCfxUnOWhFmYMlBOFfV14KihxVQ2YyCSFAVAzn7EQQsL9R3Qd1KlahzLabUQFlwerbbmHGzBL7gDiyxvIFvkWsjHAFCvY6SsAYFcnwRoWrCwfurfQnKieqZDldD8xj0IEtZmXIzb0BHzjebv4fmZBVf3g2bkCUWaaDmeF2mOX1tyVshq9CGwWAh8s+fQF2pQcJ89KDU0JPKLthcGCiK6A+DPAfidzPyEiP4IgN8HWWu/D8AfAvAbF75XlRRcbF668ETnkr+FXHSpfwe9HnNTn4jsRTn6hUYI6NRr3DYmiiTTKXfQdeL9uQBIi5F/1/Q7fw6kpGIC0eSkiWaxO/dgtnGx+1eShxf3l8Eim5gBUOZ3oj1M/m6rGsx8QaoAKSdJNAA54xKY5X7Hoku9g5KO4VLk5tL3DwJEvvdztEPWirxuDzzH86oQN+Awlu7zvO2FgIKIeghI/GfM/P8BAGZ+y73/xwD8xaXvViUF6eXDs5MnrlFFWn8L4LgJqA20Kv0oIc9TAlMCRQWNoWyGkDQGxERjL/X4knRTgg8cW6z9ac3AomXtq88wYDkhE5UiN/Z9GyN7NhcIV49fOx5Oh9dMYqTkIHWujB6l4l+Sv3oY2CqJyou+S4FR7TMbYJgEE9xr/vpH+ITKhdq3I1G/FYe1eFH33jEwOWAyPgYSMwJyqeaHfvZgWxqPcLM5u2l7EasHAfjjAD7DzH/Yvf5R5S8A4FsB/L0X6qE1J11UYAHMXb+XmnIVhyIbfVZoZO5CkveS3p8mzuoHGcAAOsku/iL7TS2Lh9WkHeuz9ct5CXKIVd2KKsciAGgsAENUoZnvgd07RhhAyI8Cm2X+ssAronkkpm/O6rLYzEzpP+/9B4IbG+9jYtaTm4jNGWAa8XzGTTQWEidtHr52KL+XwOLY2BzruwcJS2dIYTmVfzt+zfMs5XOdRRS/YHsRieIXAPh1AP4uEf0dfe33Avi1RPQNEMn1JwD8iy9wj7otgYV6aVaqSNNmocnWFsTnihNISawHGEEcgVCsGZb/0i8Um1zGwini73UIGFprifZBLBNBOIao/g85mlBAJD/fpNGTxkNWgUkpn1YiLVBOTc+xOEOJj4Xkr4AfuwP5Liv1h9OiqE3NdXKmr/KBMoawxX8ooKvhJfz/oXneg+QiVRxW+0yL7YBUOv/cgjSxtO48SCxVkF/a4PmAKBa2RWnhurX2jO1FrB7/M9w8u3bzMoK++Yc9NmFLYAHUqsghBPW+C42+NzuF3KI2fR2UZtIATWne32MLqiX1Dj1ruynNLj9NQqLmuAeTKFIeGwYcWLT9SOJ1CmRJgdW/waw9lfMWJuFJEpcyCtqPmcOQ8y7kVKeu896I+TtL4rEBhvteBo8lsF/4v5IOq0fnfJ/ZqC8BycxH5Ab8hVtHNwKJrhNvWQUKeb3hrrz0kn2HQgGLtl3DfT1Puz2ema13n2/tZpqBBZXT4wBYVIutJZBMZfA6dbOZxTXanRTGU4xTCSBq2nESqyH2jum9fkMGPUW6DlJKwAKKuKgxjWdgfm57ZnZiv2X6IhIAUVWD2CVPEecRfe5GzXCcxMwFOfFcKPAk5NL7flyaeTS/mPJZyhXp8zinpM5kooIhHFGJjvk7eOuOPFx+/TCXQcvfdf3N93bZ06UeblfNhXjRur4vrY9WWiCqf/u197x+ItpuB1BQK6rL6ikIvrCZGrAAMOMtgDmJU/3fnu6Z4HOinz/lPInJnPMdsj8prXvVJHurARXzHVBO6bQgmeRH5eIUZCY+/a4Mj9+cZay4PZGgklGMyOZge89l+zKCmLJ34VRE+JDE5dz6VUU8+mfQZ/annsu5IGMR5iZiA7sl87E37VpsDtUkg3E6WapqVY8MUqnxo6iJyFnFN+cfstiuA4n8MRljUzcqkDBTdqV+QZL9Hrqf///IQXej6vRH2u0ACjSTAkAClvRvPTUqe76donZCevOpC0VfLsPmRDPbtG3gUSggUfkh2L2r6lF8ELHzaac8SkUM5uuilFJsW6WC6Mmt8SrkgcDp/XLN2lejes+ePXUlVySRFuYl4WOmpPkrBpFczCRb9Y3nUkT98ACKVacCPCiwI2TL02yu2hNTo17BJZBPrEyEyuHJPGBRS1W5W7p2fB2NMjhlLeRo1Wu5iwIMB0HC5qsK8os1SER3cJkEuOi9SdV1Sx/q9Zo/bqrpC7RbARQEiN7tm4FEUieiCBVzzReiPW0a8ylwM94CWFY7/MQZwhPPOIlqsdlm5XpzV1aapcxLgfIprReddbE+9VIx0c7ITweQTR7FfC19TkqT4yKQeQoAsgG9VDFNAhZ2ePu8l0da248MLM7kSaaOmCn4kCUoBAEUCmqKniork/eABbPjlhbm/qA0IUlecuQlUPKaek2lkkTncRdyj2Yefah/KPzEIkhwqWRXzVG+Fs3BoZGA7XuWSPlF2q0ACr9xluIBCChE2DTJKWGA0S7UY3EiCwuGWnSu8h84YgmQhQVkc2k+2Q49l8uFYWBR1Z1IqPgStusuDVFo8kXUD1H/78CiyqWYP26cRjlp2CQKIjAiqFf+xSJNyZsrG5DIJGMzdzkgTiUKz+V4dUwlDQQjYkMNFv66Fszm1Q833qzu+sZdLFqbPKA30oSv53ojMrBVNZY4KeuDqRxWEKjd1Pk5joCEW5vs16lTWfx18D5IE8AtAgoyr0B/ssOJ6gYKkyWKnTRUeAEwjsSJAKi8Ig+lSvMIXaqUMZAop8EHAEy1qLmkN1ufOAVJ0wfVolxw2CwUeTZGQU5b43OCLbxYyExPzi2JrHmTErI3anUPdU0nybMJLeacT2clB4+C42yxNkSnk4J8ZG3epMC1YEHOQjPrR9L4H01cRKzzEagen5wYuJYmTKLIG9Dm4kAA10GQyH4tzsqR40/mICHrwKm2ngvzz28AYYeYSYEB1fWQbO+kueHrOdotAQqIruatEQvomJ2PXMr46wDj2ihUwIm8ir4paE5J6IKR36wMmS1QYhZCihmconzHxMYmaU7Vkgs79wvSSM12eIi0gC5laYG6rmQwSqircC/FFXjVIy2YD4OAo0ScQk7rLopUkUPfCSUxxwLj3loIbPxbyaPt15S3SOZPmL3J2qmGCWXD5ZNfgbndXJaFnRhk0ovxRe1aUaCqKpwfakt8xHUu0tTk1rB+eq7t0BjZgWXqUFZXqAYKHR9xyOMqCPFF2+0ACjRZllsG1xM7Kk3wlETHHsbjgHFdQhzoyc4JSBM4USl6k5Hdia7EYATJbq3lCPPnkumWqSTNWcrvaGK3j4IFZoSkPH5ZtAT5mkkSMF5HVYnKfLgIFguchrfwqKTCDEngo7ovqdckQSUDH9F5XVuKynTvyUYWScukN6JUJKh6MKDRecgxEBVX0Khm1QZ0JPcBn4mS3TxU11giM1uQmPl66Ovye2kNqAS4NCdZeigqUAaJToptc4zFWc4PQQA4JdAEAM1YvEC7HUBhEgVRbaZr9S7T3czRaRzFpdmSoWrOQ1I0le8UEftYFGqGA0cI+ZDzLF0AKBXRGRxZktZyh1y6L8Usqpfq4c7W7zewkyDqAC/UunJjUq3T4k2O4/BxJk7sXxz3cm3Wk4mp7M/sgOUIXoTkQHyuYtW5P264OLPkocSm5Q3hA/mzF4i9xapnCgAzU+cCKbno72Im1Ob59I9y7yOgOXsvMaDqp/R7rloAyPugAogQgC6AO1E/WFMf5MbI+VIzSHiu4wXa7QAK0Jz99QSNG2w20+Q0afm/SRyFLGHu4HTRmejbSBcTRH91DDlhrCubmzqECHbKHhsxxax91aI+3NVFlNUpqipKc2xxLQb4hDoWoAnx5kMh2QccvmhG2pp61XzOL2JvltQ+vXCrFq+z5HjuIPthUHkPqMYp80LK4Xg/jTYLNTPPAa41jbdOcDM15QYgYXyLgZ9ySRK3U1RINmkmNKoWEVjNpuiigEOMAhSaSKlyAJRToqinQJaM6ZpAupu02wEUhGUdzHTzGSonrWI1Ad0EGlyRHjLQmDRX5TQHCyxIF5g0ejJWG50AoFd+hONMwpGL6EYyRt6AJ2j27TadVUvQeW7C8xBOsgGFzOOQF49TKqelvVbxEcd07VqCA3SvJUYuivMMp1HelNdWUjusipj6lRe3t9DYd9tcFnnMY1HrpEM1mQiRIPkwFVtUAb3PoRghYEFa8K11UbfrAZBK8yyAQXLIADEXsi7WNwWJ3oDCfgzg3eUn1jSNzbzZeP1M8KMA0OjKjqixQak2qOSzpEigMdTVv0ynHhTRAVcQCNdKFxkklJgkQBYhazFZbzbzm8g2ddAM1wCs3ikvLW6Uhca6OMk8F1vdNpjE5UCiElUXdmVe8AcWiIGPB4nEsnZHTe83TsVvZOZs1er5jmD0lcatLzdpCwQsM9WBftkTtrba5EhQ5xM+y/HgPUMPpRzwbvBLz30TkMh9ValB+Qjy/SAqgIGoh4lTBaOTJFqQiCaJAGAn+eX7YhksXqDdHqAA8gnHC2CRBwaQwQlJLQGWcZuyBGDD5j6uYNEM6IJ0Yc46Xh0xoCAn7dTXcYARlOgEZPKD+uwvia72r4V0AyUq1Vpm9ZuF7069g82DWGsGNNXD2iSmSuIEGibQbgQNo1g93GIrvhk832zmyRhUVcv41XyuXdgGWtVnUhbRGcXVvLq//34Ixd9mwfN20b1+iUjUayx6Y7Yg4ddBBWa8+HeVJDq/Fkpmdlv/JlEHVTcMOEIDEq4RQ/ZF0rU6seZzTSVb/Au02wUU1gJy7cyc2NbYXVv0IYITw+pttOoCA4etEAtOWvKr1HbMw6p6JWsBn5wJu11gDXLngkKJi4h40NW5hH3zNFUL/eDJVZ0Y6bhjzXWLJIkXJga5Bu1G0H4Qi9I4ustwsTAlixRVaciJ+DL+U1YFZtKcBy5POJrIbbwEMVTUm4fLmx5OVLu/m+UCC6e++d8sEa1cCG757PKYHUsVMAOXdr7VwsPKJczc+f0hEIMehNB1j2LdaJeEShBF+kv1zwHLzbO02wUUFcMOBQcFiUAFQIzYmoCEgGDmUzYuoUQ/LlkhDvo3eL8LlUCMt8CUZPOr++0SkbVo2oz+MwcWV16Ukxw7s2AoNftOqVhm5MulGvbzLoQk36dBj/9xEpDY7cHDMN9QCprV/bLaEXIlMp8Wr1hfQg0WvuX8GgvEX9bjjzS3cau8kubbon2eOVnNxuPAe176WopBab/vHc1cK4mTahf+OpJZ6+DmvdAeSu56k/A5+WdMtTSRfYN+JgFFM4EVkmYpQ98DyVoLAHdKaHWhBDSppcA2eT59r8uk3nAXOW6AqAAGjcpYOw7BWmauQ3Eb9ypD87xCrqW8IYi4OiGzq7eerpWFwzZArk/xjHpo0gU1jGJ7Tyz1N/d7YD9IDc427iJ7xerJ7DcluWAtezYFCxmSJelioTiQ+hdUenwGDJot+kWpy4NESjkt/0GpzvepuvgCAF3XWpBoyGULEkR011RA4LAADEZusy1PvW5CJn1pEpDAxKIuNhLF0ee+QbtdQHFdI5MoAILU28iir5NGOAbQVP5fdPU9llfgmBnVNgLQEIr2EjkLiDjGHE1L5gOdgm4IHyCWmfI6sK0KTU9azMbHdLS8RNtSAtKUVQsaIGAzDAIS+30BnoWNnOt+NGpH24pkwcgWEXPamvWJs1juCz5lV3cLDLxJcyBhEb7V6e4A/ZlbE5kLtJLhMkhUbelwcX2nUHv5wohmoFA+ah2iiYExqTVwKn5GyfFKL9huB1AwyqJ+3odqLSP+b++4tRSHcaxr5otgOixYJZZ6cksKN5EcJG2d5tUAcur86gQ2dvoYD2GL0VWeWswqpX83nVoeT9PRjX03QBhHKdS734OHMQPmbLhaHd+pDfnVNrw9u2A7sFjymgXAgxZtolBULYu2zddprpvv60lPJ2n5HBh+mZi/wwIAtWn2qyjkhfD9PDaLfy/M84JaQYnBwaRfBiGV6E9zoXfPaqZsmiaVKAwknDTxPrQXzcL9EwCeQpbSyMw/n4heBvCfA/gUJGfmr75RXQ9vOVh6e2EzVWRQ+33/t9dLG3+KZ2r+O23MV+6Uisf5NBBJhEMERS6nkZ9w65t3jGmtDO7+BwvW5D4sAGbL+luWcbP3JxaQGLSid8rx5PN++mtXJGNzgrW+BEvNW568iD5BndSchcYDvqk0B8ycSyTeomfsTdeBN3UuvVc9SwseN9ysk6pyiYARsl44APm3vueua9KExUCRZVwzsLD2PJKTa++HRPG/YeYvuf9/N4D/gZn/ABH9bv3/dz3TFZcW5IFGjHlsRmXCu8EkeQa+vf6Rhb7sGq0hzkBx4EosSXETKRFacm+WFPROf146rVrJ4SbP1QZp5WdNVQ7O7DcwpRokPBBco9tns3JrCk7FDOsJziqISz7YSCnGSVhN2EIiy8ZMAPi4Cumd2FDSGLSu8osSWZZ0nFRxSO1pAu6q15bWlnmBVpdIYlEzwtpUjyCHTZYqPPAbUKRUfF6MuGzXx3WAfU37IFSPXwHgm/XvPwHgO3EtUDSb/FlbdutFEcfsNPG6aqunAnMp5KaE1ewRGsIKU3HgskklAsKkRGctQprJL+v+Cw5NM8a8OSVKCcB0/ARhLoWVXT/YAMvfMxdFboLtfASl3XPC7DO+Bsu1rTKbhiJBaCaokuU8aT5MytxH7oNvwVWH9zzSkt6+ZDKnkNfLoTIP5eM3W7ezg0cJZQCavT0AHIGolruoa0elCl8uQtQT5SbM38WrHXbP513Trr0oUDCA/47Evew/0qI+b1hdD2Z+k4hef6E7eALPrB5Lc+IInlqqSEUcfha1o9Kfn3Gg2SWqcaG+wgdgkfSr9EnzxbDTL9+/Ls+Xs3HbPQfhGGZFdylg5lOtksVBV2sjKc3BTM2olcqUH7eWfAr41tm48+ZsPSNDRKnips+uBClpAFyWKp7Vw3DJXR7Ia+P6VH6oAONGbYkr0uvMPspiws9ZumAqdSnxmO9sazG5dcJcSxKzA+bFQQJ4caD4Bcz8BQWD7yCiH7jpF8mXFIx3jn84ACXwQn5LuplDoqBbkNfZzd+ngVzsxoKJterfUuTorHvOsmIOTZqQldarkrl5nEDDAN5uwdudEJPXMe/XSW8aep2lCqv8vaRiHLy+gsuSVGHOVpZ6zghWkxJy4NqCadnagrl0/hxu3FsLkQe2g6TyC5gWj60v6w8wq8VCVBBcLB0sUqgWnq4yYHnzPzAHi/chgO+FgIKZv6C/3yaiPw/gmwC8RVotjIg+CuDtA9/NJQXvr15nUxmuQ2470AgmYSjKtmPRDtZ1m6ZtR0T3Q67gR75QLCZq+qv61F6j9c8wCcISsq560GYNXq/Am17eHxNot8+bmq942W0dmC/eY4ARQs57QUHzRQxDXWO0unbAjG/wLZtvVXXpNRO1qgkMSG6KtnlzcFbPbgAQWbqs84fOKrVft6GXrv0sbSnc3ceVAMiJgCPAqQmxJ6rzXnqQsP+XpK33I8oXS66CN2xEdEZEd+1vAL8UUj7wvwLwbfqxbwPwF5754gk4mL7LO6QQkL03AZRalbTsYnsT8dG8DBcat0B2ILjoKD9glg11/mHzdNRrc+YpUs2vWNM6ENx34JMVprM1xrtrTPc3SPdOwKcb0GZTUgu2xGcQCYFWK/npeiEJDyx8UXUiqO+BfpUL1SzmH7WAvEPP78V9QPrRddqXXv73BGoWp8XEySmpapXmILE0Xu6Ezd6Z3gekBQmb+9aSs9SehU87oHLkZ/RAZvPVSgyTEJazHwPOQyqZd7x6gfYiEsUbAP68LpgOwJ9i5r9MRN8N4M8S0W8C8JMAftWNr5hNijZIOqBpLjFY3AcFFsAIARaEldtNTUJHPneQnzgYldmQfEeuf0w1qDJ3x6YPCoTcR6R1xLSW63cdiYfeMIK2W/B+QJ13QaSSfIpDdeRBvDDzBmofqVMyMerJZ+rHVPiE+tldYeMDRKPxH9RrXQtAiN4yAHKNaar1clMnl0yd1QAekCJay4peY1YR7Jg1pb2Pv1bb/HpoyWdTR70Dl0rJVhyatA7uUbfxa1Wnv4+xHsz8YwB+3sLr7wL4JS/SKb1QBgv5beIaVxJFzkBFKKy4kV9L7rA3bIcSvtRp9prBnznULLh42zN4L0X/zHadBUKU1YuzcDACkhKKD0zriLDpQCs9/eMObDFdHiTWa+E3YpRNPwzA1RbY7+XE9s3Unb6DBY/x6FLlp/KM2eSrPhCtBDcjWPteASgsmoItETGY6/IGQMlgnk3K842QrV5LgWDGjyxs4KOFfg611mK2BBALa+pgKQnru/mHHatq5r+zcP0b8UnXtNvhmembt1pk6UJ1TV0c5STWX6Sb2LtyxwAaLfMR4WC6NADXmhPbdgP379za4DGTfAA5MSxtn8VNVCJtk+bfPA1HCQGnMYGGBFpZKD4kHLnvpCp5jIUzIFFZaL0GnZ4A6xW472SMd/u8yL3HpKWWR9+VzzKrzd5cuBfMj0CVBnBxfGxcLBLXmbJz/VIDSFXVfCpACc8+fIJmc3gLEo5APRjYl71Hr+FbfGssG0sAsZhXU/+vesEJsCQ8PqJ18bbL4AOgAMRN+Jxr2i0BigVyLUGcSQKDiTVElwCuPy2RovpPAOZ5Hik77TBweJM3G3xmFtXvkQeq2TWaU8U3n0DFapiYfqoEZzYPeonBwMInnx1HYD+ArvaIukHTKmouCwmmMt0/WCLiGIQLWK/AJ2vw6Rq8UpJy20u4PjOCWQUA/fwavFnJZ9W5h7ropApnfjTx2bJOIRanJr9hTfrLPhqFP6gWt33P5l9jaMifwjafjfWjSlDjQMII4SolP1BMjRbIdlMv/0MAcSj+xSJrHVgc9XNQsKjaASn5IEBwmdPnbbcEKFAeXh8uB8UkzShMyuw32X1yvETiQoB6QtNUk5wApmG74b5zyDMPqN7LCWQPnTZObajDiUO9Qey6g/Z3HJevq9cTr71BAIN2ICIE1nwZfSzWj5REolqvkBPnxCCE5GYN3qyRTldI6wgmQlhF8ZS2e+4HAZbNBnznFOlsjdRHhP0IGqwcnqYfNJE4Z+MK+fTPPiDeLV3Hr9RN1UCwcVRpqlnQzCqeUPHJaP0S8smdSkCdEc8NSGQyViWK/NkmyKtq3hHM/s9/O34j1n4ui7zCsfuE5nns4PGHWJs8Jw8THwaI9yHm43YARQuQJjpOpCYjHdgYoB4pefJyZh9Gif1oxEyfRPaw+tEsBBeDcJivCKVvFW/gPu/1TKDwJs4dmWKUJDEhSJj3BGSw8IuFtVrYOBaPz2FAuOpFPbD8l/adjXIRRn52EbzqwSc9Uh8x9TY2AKWVEKEhiJdfDOCzE4wPTjCeyKbuCAjbsa6elp8zVX/PTzdzXrI5cSSnF41JwGBGPLITJX3sh285PCXV9wQKSPR9SVRsIBaS1EaZP1WRWA6Bg0kPKqFkAFpItMMKotXlmwNtUVpwIHFQ1bBxNh8RBxBZUnuBdjuAAqgnwqQKSjJHEgEk8fg52a6eTpO8nsNsk9tgXBZffa9rdM9FUrE5xbIaYtduT8K65qgV/BH+BMj5L7sATKGuPgYsgkVF8KUkp7CpI1YXxSpkG6+w6nNyViF6LYmx0gv6rKkLwNkKpOoIryOGeyvs70VMK0LcCSjHy2IaZT/GPj6lsebUJ3uoN4MLZScr2TCpm3ulrrhwdrtOu2a0H7MENcEI2V59NqiWEsysWh0W85D4a8HBx6Pk59PkO5iq9AFLjnQZYPJ9GjVmCSScCdnM7gaylSr3gu12AUU7ECxBL5IzUXMxpkbs4hIYU2X2ManEmvEUzxJm7qSDSg2pgoxsImJtEZEv1WZOajJMG8ZYcR1brNO0LPkYyWXiugufNv6F+g4wwnLVI606pNMeqQ9ZXZNxYsSxPrW5j5hOe0zrgPEsYns/YDwTkrQ/Z8RdRN83PhcL0sOxwKhqs5vlJEYgODUpTZJMR60wfIzxtxPf+6c0gW3ZXyO6Tcjl8yUWaJ5QdwYOpjp6cLaMZ5ahyto0aR6TA/333rat9NsCRJv4SNdJBgmXwGhmFn5BaQK4NUCxIG6ZlQMJhKCFaSbxqZgNmGZPdgNV3Fp1M/kgmes4hrYfVpx3iaSKyBMzKzw0uxRny4VYb3RxAboYotYpUdAxtauVKkIqGwQoDkxB4j1AJJLFqtekOFBJgkCjZEIK+0ksJroJOUakTYe0Ctjf77B9ibB7QJjWQBgBmgjTmpD6gNDFrI9ncfpQ7ESjPthJT9MkCWlgfV2VymfTBOz2sgl3OzHbGilqY6iqWB4ftQblCFgd02wSzhYglFPWJIl2c/n5uwYgyt8OhLw0k4P96jXoOZ3KlJ+tMl7ScPyEu3YmbKssZ9P7ChDWbglQQEXoZfHKmGJJsgp4bTJnrfZRdD41nHnxHQKElqi6tptuwdgEZ09LdhaMuR6/yG63XoTWp0DXCz66WaklzQBg2AP7HnS5Q+wM1AiYGGE/IexG0HaQTNvMoqZgA5x1mFbAcIcwnkn9VUrCY3BATiNPXSfmVyMuA2qzouckWmJ2GEW1BMQdfbUC3z1FWvfylWkCXfUlxiQDv+dBDEgdSAzOacyTlzaeFnzlPTj19zF/A5vz7EJv0okBRHSnPbs+WT5Tq2R3KGDLx7MYIerBx6QKVinUQsuTSJU3AQm6yXo60m4PUPjWkIE5EQwcYWnNqRk5NXnS2g+a44EbBJ4h7bOARQhlwbiyfqTiX7Zg2ML2bHVOnkrNwnKnw43Gp+jIxkcgmKlRuYvE4O1WnFaZQduV8CGAJNLdD6DtPlsbaOgQiBDvrLToMQAmhD3QXQHdFSPuFXA1/sOsKlWiFPvbxSVUwVcmodn7XQc+3WC6t8F42gMEhH1C7COC5232+7nTkZ3WFUikwhv4WqJOIlkECL2eXxs+M/gMJLpuVjxJJJukmcKmnDGsklhQHxjUgEXO8m5Ru3qAkrccMef1naOFTeWy8fbPYGvmBdrtBApgDhYAYEFV/jUDieQWQNbdnN5qCGztWQgez7gDqleXKE77DI+jmAz3Ayg0J0grtrbJS0zqcdIPhQMWGsgCo67TALGNgAVRSY6724kL9yA+F2HVl0hTZmCc5H0DOOVJ4vka3dUK3RaSFmELrB8x1k8S+qcTwl4XYxdBm7X0IasDnHNu+hM08y/mACUPJ7zBZo3p7gb7B2sMZ7KY447RR0IHBTmS8V6SKkrsRm3hyOPriVcvvS3p8UuBeTLYec4zp9KqvxVISJkD1gLaFRh5otObQ0MUV/kYZZ5CqCRsDsoteR5mqEFiUYpoyePnbLcEKGxzUI3QS63S/3RDLbxm3MS10oT/ztEuKpFojfRUsfqQgBCJxhGovpw/vqTXWnPRpJ4YPKpjmuNQ34tlY6Nm0E6J0f0ATAPSNIH2e+ED7FR0z5QJMBLX7XB5gv7iBP3TgLAHVk8Zm/dG9E8HUVd2CgBmVQEq0Zi2e9DllYS7W8HcUCSDrKuHAJxswGcnGO6tsbsfMZxKpumuY1CKoGkFJCAEJWkdD2Un9szcaCDhCzUt5PpYBIlqHch8t+pilSogAXnt6kGBcdSxr0/5vGYS17E7duhkq5Ur+OPulQ9Cc4iz/KY3AIlDiY+fpd0SoFhoS/ERwAIgeMniyMZakiaegZvI36vyEMoEcK8mt0nSmREReBjEg9LuqwtBTKLOROdJLrv+dZF+2cFHrsVaes4yesNcr5MmszGX7KXF40/WEBB2e8TLEavzCFwA6/dGrN69BF3uhP+BSBNmUeF1ROqjWFUAdJcD4rtBFvd+qCUvohKrQgTerDDd22C4F7G/Q5g2hDCyEqcBYR8A7kVlGtegpPVHhhG03YO30FyRmk0sNuNjc3xkPBdBwqwOrbjOWjfGArXyNRRsTQVaSh4EZBK6ctm2+zh1gz3I2XhpujuREvdF1bopSBw7fG/QbgdQMOqJ8nZkoJCEC2BBqfme5RXME13s9BWZ86wg4Zvrh9VhYMdbEJDFXtuEiyeGfU5VgdK1G3Ao5mm6wHf4z2T1hS2BzphJM+8TUt0uMbrLhDAkrN69RHj3iUgIgKgLJxuk9Qpp02E86zDeiRhOxTK1fhxxMkyI232uroa+B5+s1SmM8pTwOmK402M4DZjWBI4oqhYhg09aqSmaGWGfELaD8BcpgZOQl4RhYZ4WrAzN6zPHOXlx7mRn5Ok4lgrkBjBGXA5Dsbos5uloXjJpy4DBSxJOAqJRwNGrNBUYOQvNQZD4GSNR6KbPATCtCtJugqXNZL+D1v/kIIlJzb/+eftmCyPfWtOXWZ8B5MpmRioG9RkYi65uz+RPDAZAIy3qkXkRm2egLgQi54tg5sZR72H1Qie3eLoCCN4ZR9yhg/tcJ56bMQhYXE2giy344hK828nn1mtRdwBwFzCdROzuRuzvq+WDgPXDHtEIaPOV6DtMZyvZ9Nr31Iu/xtSLpEED0G2BbscIg/BCSQPeUiffifuE/pw0IE42Dk0TGH2tYrTzBzhOQCUbC+W2AkP14Fc+Fzypw5zFnOjY+zwXVb2Tpfvnf6nwHuYM1gI+uxR3rSRh/McRnq0FiRfNm3lLgEJFXx9JbCSOZ74N+D030bZAyJ5MNuA2oTGKZ+OzmIn8ALcJTbzPRozi40ykJkWnCk3uhLD+mBu69rk4g2FmbqyS6WYHpQIuNFqVbgbtBvBu76p8OR8A5uVQciiHsuolGc465OJKNLpkL0CJGoUMc+oIqQemtew1NoCwxQ0IsKQEjgFpLWoKB/ls6glhYnRXhDAw+itGvEoIk/h/pFXAuCFMK4KE1AeEISL2UQjVENQfA3VNkVDG8qBLtSNhD6UlNCBoA7nsPhkoFnmOI81LhOVmWb3L62ochSB2ALEYEZsn8v0hL9t2O4DCVI9sPtOXiYAAsDpZSUZlnovifsCIYOGk1HXzyuatI9MNWyWKJrfQxglW5ZwD1K5OAGv9TG/JMBUjMSg5vbFSn1S/T0lJr1CdhNlyYlmps9VHT56dWDxMDAYg18kOYrEQbHlRFWI0rTqkdUDqnORjQObHWUE7jIwwAPFKwHH1NCGe70ruTk45mjWuOnAXFCgUUCegv5TnjVtG3E4IQxKQWEekFbJEQYwicWawMhOlchSJi+OWPZ+T8kpi2qlUmW/8Kig1ZlOgeML6dpAIPdLM3GoAZutpomqdZLNqq2ocSC60qC69j+12AAUwAwkAWZxHtPd1cKcFkPBN2X+k4gJ7NJbiun5R8YhkVrF1moBJN9IUgFFJTLOWEmnt0+CcwpypzFx+820466sMlFB0fyLaovdSgoHFlDJIpN2uFk05lkpltvi9HguzFlCRwAngjoSo7VdCIKpUliWZKSFeTdgA6C8lHuTkrSuE954iGUiMI9LlJQRDxeRJaY1pE/PmD0MS7mE3ZhUqrZSgFQ1BAGkCusuEeCnWFxrmVo9qrLwDU5beULglc38m9RT11zngfHXQPf0mjRzg2z3svma+t8CuSp1Jcynip7ndEqAoaG4DwVk6KNwDqzu3EEtUBq6VKOxv20hU6o9WKc+eBSwsLwQxWHND0Ei5XKCJvoSurk1qiWrsRM4OOMZtUG2y9Oy9W/AWHj1z2LL+2elj13eL3HJyyj/N4vfuwZDNH3YJ4JDD07OzmKoy8jnx8OwToz8f5P/zLejRU/DVlZyaq1W+fwaLKYH2pwhnawGCJNehvebCDAReSSQsqxkxDOLsFbeM/umA+HSLcHElFdcth6YHDJ9C0JsjnblU5quMialubSDV0UTKN9m0BsSNipGtJ+YgBxQ/k6YfB6WWvNYbolK5PnIq2Iu25wYKIvpaSOlAa18B4F8F8ADAbwbwjr7+e5n5Lz3vfUz9kHCLoAFiumiJGkcqrjdRReiQ28zaAsnJ+yxgMSFXAvP8RwUWzrfCyEsa9a6T6O6+AI8ELDlSq/LXcIBnKk5XADCTXsHpthRAK5ecJcdBFLGVbPP4xcQMGiYR/0dGvBRnLYyjqDLMko9iPwC7PYKOPQ2jbNqn50j7vURprnoAvfR9vwcSI+12CE+E54i7Aawu2xhG0JTEzLteYbqzxnC3x3gi6lzcM7qLCf3jLcLTLehy60DRS2p2ylsFdOGODJDNqzdLdY1X45Ib93Pn0vTNcpvYtJo/hPfz8P1aCqyztaDXO9iaBEzvV3uRnJk/COAbAICIIoDPA/jzAH4DgH+Hmf/gM15w/pqBhKofzJAAsaCgwSxJbbxUkezzoWwm++m6spmzrspFsjjUD99HnsBsYquThIBsFmVgPrKhSDXm3zAzS5pk0TLUSSNoxxEIvYBTr6ZWADyqGK7SCfU96PQEdLIRwm8YwZdb8OVlJhhn5jNTdYYJ8UI+Ey53wNML8NUWaT8AnJCS8C40TjnXJe8H8NUV0m4H6nqETcgeq+ZBmnXscQSuruRUH1aaY4TFNLheYXj5FFev9djdC+BOPEM3jyZ053vEd5+Cn55r1quGtEtOqpimulh00EA7b1VoA6ka9+cl12eK0E1sr91QFbDNq7EmbNf10b+HfIDaEPeDt2i+v+TY9QLt/VI9fgmAH2Xmz76QGaZ5WFM/Svk9ISoJgBjdWTmCNDeNcawBw1r2pS/EFYA6byVQFsiSqOmkC3AqjDtRcbiye83ExTC/nj89qJkSBScJj9BTcS2RltwFzUVBoKET6wKRpK97cBf7l04wbSLCyOie7hAfXoAvLrOLdY49iAE5Zf84CUBME+j8EunyqoAaM3gcwBfi7ekllpQtHKoime9I3wM9gLEr5r394KQv8W7lVY/xlROcf3yFi48GDPfE8Wr9EFg/IYTtIBLLk3OJBj3ZCBhly5JudFfnIhOwGTAKN+D1/5kL96FI4eTNqc+YU7MCGCA79Rzxt7g2TqMBmar84ZJj1wu09wsofg2AP+3+/+1E9OsB/C0A3843qWbuWeCF10vtjrJZiaP4SKSAzBTnhaOmQSMcM4FUHLBymDSzqBNpYZKcGnAo34R4ZCpHoN6XBEjdSEDt4ZYcJRVHKPYkZZg7QXlR2L6v2a2w6sEhgNcRGABe9xJ12ffgOyfYv3aGyzdW2N8hUAI2j3tsTnr0D9egiytgK5aRHHSnZk2LJsVuD768ks8EynxDNgVafIiqTKHvZVNZzAnrpup7dVGewJzEz0MlKsrRlyukeye4fGONp58MuPrYBD6dQJcRcas5MreD+HOMg6g1nQLjNInTVQgy/zY/EzQXRA0YtnbMFXrmHEWucplPmefNqdOUkwfzhKwCXtuOqSxLyXF0nS62Y17IdknHz/19lyiIaAXgWwD8Hn3pjwD4fZAD4/cB+EMAfuPC90pJwVBKCuZqYbH9ApQQRI6CFP08ZjVAFsRUNpgBzCyr0lQ2Xr6+LwKM+kSB6I+Sli8ennDbQONY6cTZTp+a+6nDE1lWKo0DsUzQYI1+9QlJJl3AVtg2ENBrxu1eKm5NdzfYvdTj6pWA/QOZiWkdAF6LI5XxGlNCPtm8zj5OcvpPEyxykmIomwXIAV3QVHtkzzklCUobRxm7Tv0diIDBVCWLWFWvzXWP/UtrXLwRcPXxCZuPXCDGhPN0BiAiDJIpPI2j8i8ruXeMEq07DPPxdXPCOZYIBTQycVj7I+R6I5Yfw6cstAQxuhwFJLiWRJ+lNaf8MZCozNPebNseXq3kYAdQaDfUs7X3Q6L4pwD8bWZ+CwDsNwAQ0R8D8BeXvlSVFOxflyd3AGExATI4ulkVLMrJUB6+FA3SBc8JGHRAFxxucnCSN0vlizVip98IKQGo61BUnnb2uSmVkyZLO1MRjZW0tMQqksJO/RlMh7Y++WAgO/WnSZ9ZOxEpu4indYdxQxhPgfFUPtNtxWnJXMdh0s2gjHuYwGOzcM0fIVuPAkIMkup/s0Y6EU9LjkE9KychG4dBuqWu29xHYDuAtjtxmUkMJJG+wMK37O9FbF8Fulev8Mb9p7gaepzzGeIeiFdJJJzEWpdkBTaA8hGhqkLManVYHgoK2dKxCBKWCWu1Eg/UVV8Rvea+PSPFrdL7c7SDhaIaCWAGEksRr21zFi1aynP6DO39AIpfC6d2WN1R/fdbIWUGb978hrcfr19mW7978CTfIwhgcPaKVD18WHDFdOLn0mBnkHD+DObdaea13MyBxlQHBxBISU7oNGWxO4df23ecqIvofB0McFAvFDaJYEwIoy529S3JViLACk6BJkIYgG6XELcj6Eodonb7Ys1wzyvcwiovrhyv0kXgZIN071TKGJ5EJC0+RAmI2wmrbYma5Y1YMLiPCCEgbndO359E/dC8n8NJwHA34cHZFndXO1zsV6AhqBOWeidyEslrs5aaq4CAJokVLG9WL0H6Dex1+mMgcbLJEkWOyZkSrKp6lijsN1Nx9rqu3dRi4te8f02vsViIquLiSF9yZP4LtBcCCiI6BfBPAvgX3cv/FhF9A2Qcf6J573CrgCFBMnDbQ3L2nWD1zGTvqu37hDKJAEoa+DTfcEuMc05WIg9oz6nmVbV2xEaMy27VcRbe7MOhza0696G1cdvzxnJfdp58ucixjhdNEzAEtSDItSkl0Ch+B92lOBaFAVg/TFi/OyA8PAc/eQq+2kpOSnt2PYWJCNhsQKeb7MVJ+0F4gBjBJ1LrdLjbIa0oE6w0AWGgItHFAD5ZYby7QuqDZvDegC43Eva+T0KM7gdgYs2eBcTA2E8Ru6ED7aXvYSgWDYpCfKZT4UxomIALG5NUpEEFdKvRcfTEt+jNvldJYiVh+73zidGkRH59ESBkshHprXS/tOYOBOLNXqMwB55DQLf0SF7KpTBfa8/YXrSa+SWAV5rXft1zXKk4oKDe6DMcDBZ8pJJFUMvG2r5AlZnSX+taZ5rqPnqStERQi/LOY7IkpfESxVT7MfiNmVK2EEjuhhF5SpYWgT8l1GWbiIRMNZIQQNiN6C96rB8R4k78EE7fGdF/6Rx4+Bjp/ELT0Zkrt9bM6PVkXfUi2mt+DVKTqPiEUHaEogniS2GPNCbQdi9m0rMz8LrHeCoemOAOcbdBGO6Ic+TTp0jq5h32gzpUEc6v1oghYT90AAMcgNQFBNtgMYJPVhju9EpWM+L5VdnQeV5UQmo9N2XiF3gJlSg65Xr6Th2/dL5HVT8TA52MBwMgUom05bxsbtm5hHsX8tyV+v+lOjKzz1/noZk5jjDnLJ6z3Q7PTIbofyFmriFXvcp/5yNHvQULZ8FGbq5RxG7TJRtHK8ak11pQN1rd9joU9q7C9tuJeJ4087kDpOhuEMnBk2OajyKrOeMo6oqL2yDvKWmAZOrMlEAhgbYjuvMB6wD0V4S4ZazeuQI9fIJ0celiPXTDm/pj3p/Wf6KycXSh0zChuxi0ilv5HADE8z2g5tTQRaRV1KAxwngaAazR9wHdqkeIAfyuWB3ocov14wnrhz0u3j3BNAWkKYAjMG0I00mHfrPOlc6Gu2tsX+2ROiB1hJOrO6DLS7XYpExAA8C1CYBsHm0Oo3Au3OlPcDEmMQKdxtWk6A6DAIRUJE5guaDQTQ+phXasoPVPR7slQMGyISKDk56QZn6bIjAKSUcTlwk0wLAWCGBZXFCHP5Ms8iIYlBRVsJjVxrymhkJ2mjnqGacTak487ifX1NSKV9ncBoiZN0YgjOqUpQDjvQdz/4I4W6kvhQGLqSoSYDVmCSBuR4Snl6Lnq8kvVyQHim5unMQoPAYAsVr0XS6US5dbxN2AsO6FjOykqjoAhKdbAaJAQnjq6ZgiMK0Iw2mH/b2I9Z0e606S5EyPHoGfPMXm7Svc+VzEtOqxGwN4PYGIMZ4Aw72I9d0zhP0efOcU+5dWuPhIkPdOCWF/hs3ugZRP3O7y+D9THlJvoWLO/ncIFufTfP6QlQV2QKQM4q3H50Ep4kWagf4HFBh2a4AixyiEACJlxGMARhUJhyCRl53kFuReCtUiuk1tk91GbyrRmZlqc2tucxDAkZj5b+MrQv3b8xhZ7+SyQC32woJ6FnwwKqWWzQlI//fBQKYiIBWzb99JWcBVFGvDLgppa5aR7YhOA6zoai8EJhHozplcaj+Ia3WQ0oO0XiPHc2jfaRDTJjr1VzHfCk7Aeo1wspGENClJfogn50jbHWizBq9XWWxPPWF/h5BWoq6Mm4AwnGL98Ax4/ATT+QW6N9/DvbVkAj/fR+wfBFhA2NQT+HSNcHWKdOcE2wcRV28wxnsJ42lA3K/QXd1Dvx+A9Fi4F+9IhbKB5rEbOu6BxIO1G0Fjp0S3msRZaseQJZEZJ+Tkza2pfQkg8nyj6svclHmYgL1Rm8V8qFrUZmZ7jnYrgILBYnvXE5OBqgoTW/bjYRAnmy4CQ8wxFdyChVpB5EQN6pAVZRP1joxyJGduSxmBlsChfS/b61MBDMshMIv8cyDRJBfhpT45FYF6cTRKmzXGOz3GTUQYIrqriHA5IOxFXSHlFjCM4O0OPOxLYSAiIGx1vJwJ1rgPQETwVY900kvA19UOfHmFdHkpal2/ktgMrXlK2y3S03PwsEc4O5FHG1JORpRWwLQBaASmFSR7ldUvHfZIDx+hDwH3p5fRXZ1g+3LA1BO6K0Z/oWPXd0ibDtMGGM8S6P4eA6+wexCxv9ej22gWrd2udqS6Tk/nlHkk2g/i8TmMEt1rB41mLscwVgmEZ2HfCxJEtVY8z2Tz774LoITSP6fJFUD2d6lSHLxAuxVAAVaiEVOZVBW1OagHnzLpGMfioORNiq26oBvOHJMAKJmErKZUAVHXRdpVwWYLn2vTovkYgmPkE/nU8pRF/Cz9AELGaMJYWq+Bkw2me2vs73WSW3KUBdhPCdgNsth3e/AwiHVjtwN1HcKdO0h3T2A5LYhZTKQWYAXA8jlQ3yPdWWP/yomkxNtp/svVSrJ+v3QP46t3sL8v1odNHxEuLoGrLdLVFvHxU3RdxLojgIAwBElSMzLWjxJW71yBL650bCN4GJEePUbHjLtX93BybyPPNjG6p3vQ5Q5mKqcRiNuAcRsR9uJ5WjVP+NEBScJ7U5pJfZpEJRtHNbuS6k1KFu8HAYkcxu8kB7uvroVKkrlJajqfbMdf81Dz/c+HZFMcm82blK+/3jXtdgAF4CZNffalNJi6VnN2GbasP7nGQtBaCG34deYKygbOafyX2tLEtWZQoDoV5HNzKWKWQwBYlkKsuXoROTYByBGquex9jKD1CunsROqC3tVTdwfELZXnHScRv7c7pKstwAnhzhnSg7sYXjkFE6FbRUQi0NML8NWVBHYZGahgkDY9rl7pwAE4C/ewDgG02wOrHsPLp9i+vsbVK5KBqv9oh7v3Vzj5wRWmL76N6UvvIgwDVudn6L+4FvWFSIv77MCPnyKdXwCBEHpRXzAM4CdPQcOA1aO1cCMqFWGQkz7sJ6zOE9bvRoRdj/6CsHkvoX8qAOnjUso6uEbsbghoSlwsSlZ2cRjV72S/nIb/EEnpzbWWG8OT5mYxcWrfQU7jkBriAPGDarcHKNqmE2xusjPAoBEUlcvIFgeH3EA94Pa3D0n2eqU1EzeNIGyohINNTZzZFNoEGS2STCYWAvkZcgIZQBbsNIFZrSUsVol02mN/N2A4FVNxGM1caWBVUqcBQDg5AV57Bdsvu4urV+RhVucRmz6iM3VjuwOPGtilVpTUB+weEIa7hPFkhbOT+1g9HkBTEo9MEmFnPAN2LxHG0zUovYH1fo/03iMJ4Lq4zHNhJutkksFmjXDnTNRHNZXyOALnF8DVVtRQH9a9XoOu9lg/HHH2hYBpI4lsTt8a0D3ayjN4i85NTmU3PzkXqaoi+X0r6LMEEq1vQ752yD45GSR8FXWb/7wWUpYAbgQS1z1bUoJn6bB7jnZ7gcK3HHk3FcAIlMlPr+Mbx1E1E70WwKEVybJLsKWyD0ES6WZOwSW2BbJDT9ZzsyRhABBlwfgJaxfINAGpK/5jDvTk8d0iXPUY76wwnAWMJyR1QZNkiKK9iMzsEthQjKC7dzC8cQ/nH+tw9br4n6yeBkzrgDNm9EpSksZSWF85EIYzwtXrCbuXCMPdDidvB2weTYhXCZt3B6yeBuweddg+EBUgrQJw9wzBfGLWa/DpRmJCokuKC4BP1xjOVlLi8OE56NET8NVWolEV5CqPUXX+6p/scar+HN3lhP7RFuFcrToALI6m2mRLG8tbCixAzeY1JXFFr4jpI5vZe0l6AFKVOBch9mvApImFfBjPDRDWzIT+PrXbAxTHRCefATtpkV6LDrQoUXIAcSA1+1JQTduMSMrZgYzkIy7x/a6vBOSkuNrBxf5X4cpUM+IiuUySST+6KNP8dQKolzT5906xe6nD/h5hOgHoAggjEHdJIiy3u6xGcGKEVQCdnmD/YIXLN0giM3vG/rFYjbrLNbrHJxJR6iUrHau0AtJre3QnA568tsb+/grTZwn3PrvH6p0L0PklTkJAOjsBryPC5V5KFW42SC/fxfmn7+HJl3fYviqkZn9OWL/H6C8ZKUrezNXThLM+CMcyJcn4bWACJeUMaJkRtgNWj2UOwuUe4VykCVkLrlhQ9mNBVmnryW5EdvvbXNDtM453aonmQyDhyxGWeBH93IT68HKWskWAeEF+4f1otwgobo5+RQ+cYMVX2V/DE5huwJecXGYqgQUPhQIWPq5/VlDFTK9J1CHpQxF/qwXjmG4y92JA9XPnRwIUVUafizZr0P172L52iqtXA3YPAARJ7EKJEYZJzJf7Qc2D2odeqoiNpwHjHUZ4ZYf1esBlf4L9ufAc63sb9FdiNuWtunbv9uI6zcDJnR3+kY9/Fl2Y8Nde+TQuxge4+3kCPRZ3cPQd6PJKfF/2EhBGL93H1cfv4u1/qMNL3/QWvvVjn0GkhO9679P4zE99BPTWGnEnbt/rdyO67Rrx/AThauvqdXJ2gqq8XhNAwyTSyVYsEJVDmsXLMAu3oSHvy1KFM4X7NWPlD2xDW/Su5yV889deMFNK1fe5JHm0rOFNAaJRoaqcFu8Td3F7gAJH7MtHNnv9/8Kp0Xx/9tY0v5+XWLLfBSc4t89iabFTZBKHKQ4GFkl+mzRR3UBVJd8P88T0zTwyAwm5eO8Uu5d74QPuaHYvEpGfRiUxNSeHidO06sGrThwJg7OkTASwOEKNd3uE7RniNAmjr2MWL/dYP9zgYoz43738v+KfvfMEf+e178Q/P/0mbH/4DjbDADo7xfZnfRy7lzusH45Yv/kUtN1jevkOLj7aYf/pLb7ty78Lv/XB5wEA/78734//aP3N+Fsnn8TVow3ookPYBQxnQVLjdW5Jmmu1kb1WmRwoJnAjf2MAo6tOdxrGImUeUT0ql2drJu2ZGtmarG8SvQmU7y2pvEsA8azSwwEgqNIjuOd53nZrgOJ4QdWYH9S8KW/kAnsTdM4g4D6q4GGJSSi4E8q+E1y1r5CyC3bOQWH9bROgLLVWT7WWXb7FIjJtJHw8dQAYCHtxz+62SZPTKi+g4eGSM2KDFALCyOifBGzf3mCLDTYPA1aP5TbDWUS87BEuVqD79xD1ueJ753jwIxs8+roTfMP6CwDu4BvWa/y/f/4fw6/63m/HS999iuGN+/jsL+/xsa9/Cz/x9gNsPvMq1u/JM2xfkef9exdfhu/a/CQ+1l3hIr2Ee90OJ+s9tt3aRbgWMTz/piC1UnrnWq6Bb4J6JD41ncZU0ARwKBIAZKNUa6WxPuVcG31fyEYAkqC3OZSWOAR5Y3lep6mopTnz1gGS8nnUi/ZZDq0vTs5p7/na7QAKQm2xWCJhTK9X0Ji5XzftxolQD05QsUmzMuAUY0Fop4KwSgfE63I99x1LgEImhfjny44+EAKv1VOZgaT3HRO6q4TV44i4JfQXjNN3ElYPRU/nYZDxW6+lynnfi1rAjNWTCXc+H7B+KIFi3TYh9QI6qRfCjU9W2H/sPnYv9di8s0P/fZ/Fyd/b4u7XfSW+b/86vrIXC8bXrU6weyVh9+Uv4wu/cIO/+qv+LXy0k+RD//ef97X4roefxg+8+TrS507Rf3aNv/zOP4j/5vWfjddeeYqz1R6PrjZ4+N4d0Hs91u8GnLwjhZDD+U58P9RJrWQjVw5qHMU70sZd661yFwtAT+qdmr1a59YI+UUZwKnrSpKaPP3q18PhONHYrp/KWufUUPv8teuuadepDq30vZRLlAgHpe0bttsBFHD1F1sXam1VVScSNeAYOWmJUIGoE/uMyG0Trin6ix9GV4m8OT4C0HBrDTSZ9HRyG9bcoStfj3GSgrvMoP0eKTF8VGeOU5kmhIsdTt/u0V/KtHUXE/qHUkeDnz6VXJRdBzo9AZ+dIK20Lylh9XCHbjuBRsllsX11gyef6jBupABP91icn86/bIWr1wKIV+hXPaaHj/D691zi2//8t+FP/mM/hP3U4fv++lfg5R8Efvw3Mn70l/wHAEqGst/76g8Cr/4gfk33T+B//b5/AK/8vQmrJyP291a4+MhrePIqYVoxTi4Jq8fAybsJp2/tsXr7HHR+KeNgDmh24luw2m4nRGdSdQyrepyA7Fo9K5iTTc+hlvIsaZDLPSJgTbquxFT8TGpC4xP03O2GBP9SndSKY5umwtI/Z7sVQJFZbZ8Apnn4rGe5kF557cCkuUxUVlsyA8azgoWLv6BxlMCrGJDrauYNrYvNHMJIM1BvVoAW9U0udJmS1LQIXRRycy8V0DkhL2xLOUddB2z36L/4FD0gxOXVFnx5hUl9IKiT7Nt87wzjy2eaWDchXgyIF3vEhxegYcT0yl08/UTEo68fEbYBd38KCG8/BE5PAAbiVvgP/sgrCNOE8JmfxFduP4af+P6vwaNfdoEf/vV/5Nqh++4f+3J88m/ucfr3vgAeR2xOT3D2yl3sH6yR1gFhYHTnA+KTHcL5JbDTuJOTE8n7ud3KfK96yROh8yH5KxJC12mSHrVMmboyaYIgs1QAKs2p+7uBjssoViVgXvKtseVwUy7h0AZ/FvVi6RqHUuctva9gYeUQF8Ptn6HdCqBo3ZiraljtgDm31MqRCiiTWJkd1fY9TTVgXBfT7+4n4ELZgclUD6svypZyPkeKpuIYpqIf9xFp02M8iVJ4N8j34y6hC4Q4aTp+N6G0WoFONqCTEwGbEMT78PwS6eJCqnGp9EFdJzEWD+5hePUOLj+6xv6MEPeMky9FbL44ITzaAilh99opHn8t4+Of+hI+/9lXEIcAdB3GV+9iWhPiFhhOA774i17CV/5z7+C/+Mr/4Zmm89d/9hfjpf9xg5Pv+SGkyyuEl18CiBDfeYyTN6firuzTzLU5MLquqFImwXnP191OgwO7nLy4mk8XEJVB4mSjsULq/u8tUfad5OZRg6lmUajtulEfiYM826HI0UM+EjPy+4j00Lw/uy9QCh29QLslQIGC7JrFuhS9mYNFRRjaIlrwbMt1P83LzgADkzhwiV1TrnWM7IRMMKEpO6ebmyz34ziKeKwFfkifiYYRnKSeJnfi6JQ6efDUE4CVlOfTuqSwE0BPPqx68GZdHLIs6Q0gp2XfIZxsQA/uY3j9Hs4/scaTLw/YP2CEgbD/XI+4O8H63SfAZo2Lj3SIb1yCAaze7rB5d4v08l289/V38PDrGNPLA/7dX/in8S1nl881nd/zX389PvG/PAWYEV57BdPrD8AEdO88QXr3HfCUEF9/Fenlu+BVJwFXJiZH9ZDUNHmsxYeQWHxEpgkUumq+uJeU/9R3MtYUspOUzGEArXrNhSnBbHl6s4t/EqeqlEpGssGZaRfWRAYHl428IkRl4VQ+Ezm8P8mhZWsLnOYAofeQXwcsgkvtA8hZcTuAAjQHCdXluS3kY40LSJADjAwOlpA2JWAKoCkKGQZkc2bx+ARwXTyANXOr3u81gMipSeYVaX4Mtnj6HjSsJQPUmECJxMOUIBW9O0JadwiblQR06TVgi1ajQQu7rwSqRSJ2HejsDNOr93D+yRM8/NqA9LXn+NqPvIMnuw0+d+91bB6usPrcCnznBOcfJ3zk5Sd4eHmCs88Dq596CN6ssL9P+MTP/QL+v1/3F15oNi8/OSK++Z7Embx0H7uX1+COpLbok1OEzQa7r/kILj6yQhwYJ+/s0T2WoC+aUBIrj2P21Mx5O42zMBVCVTpJ8AvQfkR4KmZSKXakElqMAri5nIBWDh/G5apj3m9Cm4x3tH8K12E8R9+XA8/WSlaJ1PLgJU7N2UrBrcMm4rUCiDby1LVZur1r/IeetV0LFET0HwP4pwG8zcxfr6+9DCkn+ClIXsxfbbU7iOj3APhNkBn6l5n5v71RT8jxE9FY7VgIwxBcbQ/kfU0mRXiAYE0+a/kDiMAYpS6o3Q7ItRkAoCrOshCeXJWWU0AyosyHh5ey9Cq9DBJxSNs9Qh/RacakaSUSDU2ifkheSAdW5g4+jgjDCOz3pY7FdlvuQwEhRvDZCbZvnODxVwT0P/cRfsc/8J34VXd+BD82dvh2+tV4/EMfw/3NGtuPnGH7NVv8yi/7Xvzlt74OX1w9kE223aO7YPzKj/2dG03Xsfbj3/JH8Qv+ym/Fvf/mCWhKGE8j9ncDwnCGzf4V7F85w7s/e4PdA+D0i4zTN1miQ6P4UuSC0vr8AhABWHWZgKSuyzk5LNweALqrCT2AsNuD4gV4RFFn7QBKBSQWA73a5gECyOuUFHyo7wWEVn1NVrs1SAZIIYJIkwpNU7UOs6XOr7kbRp5mKTt7/RbgofD8WcKt3USi+E8A/HsA/lP32u8G8D8w8x8got+t//8uIvrZkGJAXwfgYwD+eyL6GuZr6F/bg1mkcxKEgUSkAhbNns6p+icr0JI0UnEugUiF8VSCvzQtXjWYx+otNLEecmnjSDg7CZlnJE9JpISrrVRHHBPCVY++CxI9niTXAW0Hicz0RWw4gfcTpmEEXW2zZSin1jP+BACfrHD1Soerj0/4lk/+EH7L/S8AOMU3RuDrX3oT/9P6Y+BNhydf3uGf/3nfhd/50k/gq9dfxO/46m/DSz/0Cjafe4rVOePLV+/g/Wh/8t/+g/h16duxfm/E1asBwx1C3HeIl2e4+sgaV28wphUQPgfEx1JPlO+dYTrpEcYE3prfhMwZnZ2WIkSAbMgYkFYR42nEcGrOcEC87BCs/qvzleCsxhbOw7KBc8MPHTTXez6t6ySlv4IErzsNUpRpIU3ok4GPCERjfVjlQ+5IXRkj+tugMh0fAEXFmbSMgwcMXd8vYiG9FiiY+X8iok81L/8KAN+sf/8JAN8J4Hfp63+GmXcAfpyIfgTANwH469f2JNQgkSeVqIBEdK/ZWEWIbTIlmaBJEX80MrNRWWxgb5jxp3KJXUD2yoQbY52f00xq+718YBxB2x2ikWnZyzAVsXdK6lEp5GVOozcOGgsiizRs1rLAB6noLQsSuTrZedriTtjg4XSJzzx+A/1TxnTW4+HPSfi/vf53AQC//HSLf/VTD3H5+ivYfAE4e3OPv/DuP4RvOfufbzQ2x9qn+zs4+Ze+gM//lU8grSWHxPo9QlpHTCsCJUL/FDj90oDw8AmQEtJJj+m0A+8TgpoqKaol5O4dITuTAi+wzCsxMsmc58/4Lk+MtjxXKo5ZVdSnn2dbC6ZyGH/Ui3STUwPqx3gEEBmUgpjZM4kr5l35oFWsQ5ECgHnUqapaGTisZbJeidspZsmIpkkBg/6+eWa+YbU7mPlNInpdX/84gO9yn/ucvnazFkINGB4sPIAEV2JQOiEnB8Rjj1ODoFktKWrKYS+7udohfXMg4fMKuM+SEVL2YzyIRh7Sfi9gYieSOyFylGEXAUQxCTKLJWQYc3q3fLpuNuD9gPTkCXi/Rzi/xMl797F5s8N//9mvxT51eGP1BN/98Mvx2b/zMXz8zQmXb6zxT/8j31M91qYfwXtGeHiO1ZjwP/7oVwOffHGgAIDv+Fn/Nf7R/T+LL/74Kzj9XIf+UqJc456xegSsHjM2X7gEb7eg01NJb5ilrLFUKtusxS/kdCURsqPElNB2j3C5QlxHQWcG+ssRtBUx38zT2UfHSEvzx2mtGdacGpwTJPnDwni06HK3BsDXUyl8GcomJQcMZGqFs5h5DqL189Aard4jGEDm50izstthI1XNgvytgPEi7f0mMxdo22XbTFtSsCpUQnOAyFXNg5M2ANn0IA2voPntbBIsx4CJ7YcCcXwfvY7oTxcvAlZqSBBLiKXzg+ifoiYM4CmAwpjBongEdiXM2S8CILPxUEYf6zXw4C6mOydSeTwGTO8+BD9+gtOfuouXTu/hyfY+vuOH/kGAgPVDwhs/lrB6NOBLP+cE/4+PfXf1jG++/QBf+eYOfH4BfuUepqf9gal9vvavffVfxL/0uV+Pk3cYZ5/fIl4OiA96rB8yzt4aER8+lTnQ5Mnd1YR4vgO07oiRkOmkx7TphCkgEser/SDJd5iRNlKwOGyHOuQ8O2+JOZqRcjnF7GfhpEsrKYi+l1yiRpp6tRiYB/ppvE0VfTzpqW6Z1rwk07QqNV52K1/JwWEkrHIsWbJmFtXbSaQ0TprXs+FDmIGr55/H5wWKt6wiGBF9FMDb+vrnAHzCfe7LAHxh6QKzkoJOPKwsHQFFmiDUkgRQxEknFZCeHHlB6P9scQQWW7FkG2+v70Aix2vkvpZTBoD0IQJgW5TQk4XEaStJQl/TVxFI9VxXlcqLx0CWXAgQKeP0BOMrd7B7eS3Fh7uAMAxI5xcIn3sLD4YJp2+JsxWNCd3FiHi5x/YjZ9h//dzcuf6RDVY/9QVwDLj6+Blo84LehE37ZaeSFfvO50f0P/kl8N1TcIAUKLoYNWGyLMNwuQddEcLFlWTdGsfKXTv1ATQGROZc5QzjiLAfEKz84X7QCFhXXb1zlg4fU8IOLEyP13wcOeGw1nO1jOLUivC21gCAVap1r0t6vVTUS29J8dnPjAfxIGH8R98JB9KrahOLRE2T5CWVoMAkVpxuAo1RpLJAYvVrAw6fsT0vUPxXAL4NwB/Q33/Bvf6niOgPQ8jMrwbwN6+9GiGfpnmzxFA875w0UTVm0UmZRdyzReB/shPUVFSORXHTXds5z2SVwDuEZbFxDiqiZwrPAHVHFjabUbn02smhYeBWlYqYswhZLUoVgdPpBsO9FbYviZ4ahhNsntwDXW2RHj0G7fZYfelETkJN55ce3MXFR3t845f/WNXd3/CTvwgv/0ACP3wM2qzFp+Px+ytRAEA4GRF3wreMD05x+XrE1BPWjzr0pxuteDaCLrcy3vuhOAvZe8MESm65DoMUGmJGYB3vJEmas2Oc6fieE/Inuv5dmUBjyAmMsV7pSV4cw9hb1hJnKUWuNVXXrgDCHMayY97yGjQ+wuqpZJJ01eU6r9y5Cm0pABNn0zsigQbZOzljFy0kc3rGdhPz6J+GEJevEtHnAPxrEID4s0T0mwD8JIBfJWPD30dEfxbA9wMYAfy2ay0ecpfiKWfiNxEQGy4iSX0FRilOS5PqaKOZolKRJrxDVnW7UNuu294YSPg+zUAqHWaRnTTAJhEwgZ0NvnKmWq80mzWAYZK0FuYoZmARKGe9zkVpCEidpO6nrkPa7sDn56DtTov5rkH372J86QTblwmje97vvAr4a9/x9fiKH3wEvroCpgl3fvgx3jh5CV81/lb8yK/9D6+fthu2115+ivOPvwbg43j4NSs8+nq11MQe/eU9rMZJUt8xZz0cqx6074TfubxCeHqCGEPZfKGJAwHUh2WQlH4UcoV423CV/0R0Dk5KZJKpEzo33Hfgded4BIZUBnNryqQUN++Vqd47AuaQdS7STI5fMlXG3d9KJfRSTIn7gNQFqbzmJKQQGElBwWY4k+otefuc7SZWj1974K1fcuDzvx/A73+mXqhE4VWOLFopBwGTGhK7yuLIKfdp0pPYfCfGxmHG5aRkTAAiyErWV31x6kSlg7Kw1RYV6L+3wCjXp5SqLpDnpK6T/I+btfgNrLUieDLySyWh0bk7I2RTav90yOPRP9qJW3e5sWwU1bV5vcK0Dohb4Hs+82n8suGX43JY4e2/8RF88ju2oJ/8IhKzuET/6GfxEjNS/wo+ffJb8OP/+z/6TNO41L5vf4Xz7Rr8ccLTL1+j+8aH+D9/9V/D4+kEf+KlfwSUTvHayOg/+w746iqf5LTZiJ59kUSt6jpEZrUEJfG0PNmATk8FXMYJwGW2MOUygQbEm75IbIP4L+REO4GAMRXLxE2brj0A13oJe4CokuEAbo5dCw1Hp0S+8RNMkMpkwbh8hla8mv801rnnabfGMzMTlP7kNt8Hbl5zAGHIvSjm2Yazu6h/gxWwRSTkIkDu5AYcuWSOVJah6FgiEv9EzowqJ5tyG0EWL/quSBIAaEgiXu80lZwrI4igodUAwuUWtB/QvacbZrcHX1yqPk+CX7YoNc18/2SP+58NWD+J+NLf+CTWTxmf+rGnCD/yU0jnF8KREOWxOv844c5Hz/E4XeF+OHnm2fxZf/T/hPW7wO5lMdeuHpPU8vg5T/Env+H/iZ+72gAAvuob38Lv3f9K9BcnePXRHeDJU6lwTiQSxWYDGkYkVavCOIoEllhI3btnSHdOwJGE3xiGWt1Yr8Ena6TTNdJJB+4DkICwGxGIpG7qMIg1CsinvDlH0TDqgeHm2pyobL0dAohGepjlbW2yZWXeimvgIQc4wkUw0BWrAU0q0UzqoezV7vdJmgBuC1B4sEsMBEX75HRK76Y9lcGe6X+zXA7NIFm0of7OXpVLwTb+dHGJUJbqN8xyHqhPBSwnxGpViCnPXgMIW3XT3g+SSXoY5mpHjELebYvHnRG0maQFChs+DkiXlwiB0A0j4sMVTl0QFu32EhNBQap7BQKt1xgfnODqa3f4wW/6f6GnZwOJb/4XfjPW/83fwqfWf1s8Sk9PQXfvIL16H2/9Yw/w0suPM0gAwLfeeRt/8lNfwI9+6tO487k72HxpI4BgJshOA7nUpX16/CRHgAb1q5hO+jKG2cTcie/F6QZpIyAxnXSYVuLmHTtV3fa9eHhq+DmSjKVZWwgovg86bvnHr0HPdViagCX1ok2l15jmCRCpx1lKeJxAAyGQ0KQhMmCOhAnZZyQomUnDVDgfk6wtMfALtNsBFIxCCgEVQPg4jmwDN8ckNyE+YhTATNSqSr7dpEstCreTeyiYzN9XRV/aSBZq3ogEYUWVaWLxC9juRQTe7cszWR9SAkYGMNTXXzolzOErO5VJIBUSA9udiOGnGzmFH5wi3D1DfPcR0hMN4Hr9Vbz3Vaf4mk/+JHqK+J+2wKe6c3yyuzO/F4B/7Z2vw5//E/84Pv5XHgE/9BM4e+NtXP7Sb0T/dED/Y1/E9PAR+PwCXRcRhge42K/wpekCr0bJz/lfnH8Eb1/cQeoYw70Om7tnIKu0vh/E8tD3CCcbpCTgx/u9zP9mDQwjwn6UE3W3l5T6mlyXTk/AmzV406luTwCJ3s4qus8sXDa/U5Jyi4DkC1kIQ18q+pPfV5CfJc9dMskrP2FMHsWUXddJ+TrSmrlhYpGCvQXG+DkzxWbz6FSC3A6Rp8/QbgdQAIKg/lkaQqh6YO/C3ATu5Obt0gfud7N+3RAg3D3NFp/Z81UvIvBJj9Sr5SQxwnYU6chiDpJGuWp0I5ELMJtcJfQlT1GTYFYrBPPeMxdjDcdOp2uM9zYY7vZIKwKlE9An7gIAxtOAy9cDHv2shHsc8Iv/7rfii9/zEQyvjfjxX/7HAAB/+XKN3/mnfiNe+16JLUiR8NojrWn6FZ/EcH8j5NpuBDZrhHv3AE5IL90FB+CtLz7Ab7//z+Ajmyf4sfNX8f2f/wjCT57gzueA/lzD5dVN20yDYCmQHACkCwhYTBP44hKh7yXilln8LqyK3EZKBPCmF1BOQBhYaokyEHeTBKjtJZvWzBKiTnKASLbcJKw92JzaMFufLUi09yT3tx18yreZCzhNLNG1QK1+T8rLeaerdr980GTmT0/jgnhO38sPa6pFU6pvJva7iDu2/w9IGL4dlj4KMFybw8JAIntd6k8vcQBJzVtpJVaLMAhHAmaNENUTzOzngFo3UikolCbpjtUKQUTloGMRjBvTzVdIvag43AdM6yhWkk4KBW9fIlx8DNh/dMDdV5/g3maHOHb40c98DB/5q4RPfGmPH/tnI37XW9+AP/u9Px+f/jPAlw07vPmPbnD5lXt07/R4/Xs6rN8cQO89Br3D6C1wbRRPQzq7i+HOCt0VcPrDK/ztt78WTJK2/96XJJXfyTt79G+fA0+lchhtNuCzE/CqFxMgUeEF9pRjM9LTc5A5VtlaWOnze/5nTIiTWjkSI14NoKu9pPj31cXcPAIoG8z7ORxYM/W6cWrJ4UVX1hK5bO3UgJKq2gwAQVzh7fWKK3FemVXxZAdYL9JuB1AwDpuSWoDwzlLtpuWpiHGar5DU8Z5ReIjFArFAPZjPAhLW/ASHAEvAI0SmstZB/Por77phEAch89eglYicRJl85WGc3asyD3YdaL0S1eLeCYb7awx3I6aVxh6oetddJYQrkQb29wL4qy7wf/xZfxv/1L3/BZ8fX8If+IFfhtf+RsCDv/IjePi//Up88z/0fbiaerz03T1OPvNTuPj6j2I8YXRf6vHy9wH3v/dt8BfeQvLlCFe9AJ4CFk2Mszf32DyKmNYyBt3VhNXDPfr3LiXt/9NzkQhONiqB6UYfpuI2rZ6SOYCLVbWy1IkayWnOSabihbEQezkA72pXfC4MJAyA2zos3kIhC+i4X8LSOlEfm9os7w64HLwWs8WqAi3Hh2R1/ABAHOVFnrPdDqBAQ0aanu5UjBlIHLwUF/aYAnLBIF+/IRfiaSbCiE5/uZuChHwYPhw5i5CT5qGYGGFK4ETZQ5A0xb5Vc88Rkrb4lWzLFg09fag3taIr+Rl0g413Vtjf67C/IwlywqSekFeM/ol4asb9Gv3TiMuh9Pez+1fx6LMP8DXf/wS06vH2Pwz8u2/8Ffwbn/0W3PsJ6WNaEbpLQv82cO/HL4FHT8Rf4/RE1acE9Cvg/h2ks40UGXqyRfeu5NNIZ2ukVUTYjYjvnYMfPsZ0eSnWjr4DhZjTC9Co1c92g8TvdB1wIlKGVRSXMPSkGyxmNYt79aZkgMBZh6dBQ74zYZx0kwKV78wif8H17+uazSFQWdmKWX6e28IsY5V5E8pHGDwZid9yEEsq+RLZ/hztdgAFI0e/WSCVF5sWsykvtYXJfaZq0o3K8UwgsfhcrGa2COzU9XjscuCTWTtKRizUIrCx+O0JpiQp1utS6yI6N2NmhJER94w4CEjEq4TVkwHdoyvQ+SXCkx4P1q9guLfBnw4/H9/56lfjrffu4f4PRNAw4eEv+iT+qV/wvfimdY+fevQAr20n8Gsv4+rlCO5EQuE+AA/uyVh2EeFyK6TinVOMr91FWkV0j3ciMWy3oDtnwEkPcJDyApdXSJeXFUjSqhdHob1agq524hCWGDg7BZ9IwWPa7oGwBbYoldht3ALEg1HNoQkQU+hkc+tG05yp4LJUeTN9Q5Ifz9Hqrml9cXNmc5MjOn3zEaphYX36+5grePa3GYsU0RCn71e7HUABJ0V4k9ISW+zbId7BZwdaiutfWgiL3XrGgfZSi8WUjFIvk1IChg50pZ6nmqeCr65UR+ZsVvWBYaaK+VwZVUEclyiFEgP7EfFiwApAvApSp3Y7IlyNCNs96PwKfH4BTBNO9wM+evUqHn9xg6cvneDuFePu50c8+dr7ePdbL/HvffxvAADubnZ49FUPgK9a49E/wJhe3mN4e4X+fI17/BLiVmIKwrpHuNwh3ZHSgwDQP0x5o/PJGsODDaZVwHqUBLnivRoQVn3mZsRMrNXYt1vxdTBfj74TlYRZ/Rz2Okws0bu2gQMhach3UItHIBLUmKYshVHiivieBfyZD4MD8KNFru1/2+hLB1KMmYi04s0H17KugYrTsPurtOozaR0scvyC7XYABaM2Kfk0ZC1IHNq8xzIDLSX+sIE3MVDrN+QuJTcpz9I8T2G+F2lClTAYsrAt1ZuEU8esY3tTXK55CRTSy+z+SmSaJCH3TAhXgxB4RMWJy1hxZpE+rq6QvvAW+ifnePXzD5DunyJtemxfW+HNX0j4zC/84wBk4/7KL/s7+A+++R/H6dkW3/rJH8An1+/hf3zvq/G373wlptUad77QIQwJXe88Gwmo/GNWPdKdDfb3O0w9IQwrxMcnoKtTxJU4UlHXwRLnYi+xHDyIIxlt1gXc87rgrKLKBgxiHlSXaiFuSd8SbiTCrAWiijAAGqmI6LY+rMaHOFOguEAKaDBQgUXliWsHhfPZyfPeqLfesS8njbaNvxSb4sHCwMC//gGABHBbgAKYEzA5p+XNAGIxC1Wb+8HESrun40Qo8Twy9bkeRMRKW4TQnArAIM9Vbp6tNxQom/VypGNKEhw1lDqiOX2/fwYfj9AOzaTSjKa4R99JZfGzE4T1GunRY6SnT4HLS4TH94BPvo4nnzjFt/7C78KaSnDYP3337+Ly56zwi+/8AH6x+kv9jpc+i3/zwWfxH66+Gel7e5y9JRu0Y0Y436G7u8Z4FpHWHbqTk7LIGeAIjGcR0/0TxPRAJAMiUcOursAXl0j7IfNMmYMxJ7m9SmOuWBCmCdhCPrdZIexXoCki9QG8QuELGBpENYGmvng+toE7eu+cBi8xxM4KISWJ1VvXfc+k03DEktYuF3t9mop3JpDTFMBe8xJw7jPeV/XiWLsdQHFMiljiCKiAwiIHYeZCc5nOf4dqwCmLaylHez53/QOWhcQpSLZuNc1aTYVDHnmlRuiq1N00Z5m9FB3OSWX1Oaz6OWJAurPGcH+DaRNgkUA0McLAiFcjuokFcLqIdPcE4/0TIADd45Ws+YePMhCN99Z4+umE/+sbfw1A8aD84eEVrMOAT3Xn8MV+ftcrP4y/+tVfhR9581PYPNQaqNsBdLlFvDzB/n6H/UtrhO1dhCeXCBc7bL7UY3+/ByXGdNoD4Y74lOwnhMcXwPlULBGmkpycgNZriQUynwE1bRopmYsy7/YIF1cIqx6dWj3SKghAkfwgUglCjFGT3Gj8DicgBQEE6N953QVksGASl/6lwyXVBXhyFni/Vpq1Lpa5Kb+vUynST4zIKd2uM3VaBLOFzfv2AhkEbgdQAMvJZI6J/T5V+iHpwU6hLpaAMwDZxDSVibFNzRZ0FIq33M0fggGWnBMiojp/jrbvEQXc1M0bQHGW2Q/g3a6k0bNn088A4mW5e/UEl6/3GO5AnHFGRtwBqwvG+j0gXoklAKse40sn2N/rBUQuhppVjwHjJiLdnXAnFJD4zP4Sf/gnvgUX+xU+/8ZL+I2v/M/ZDfvPnd/Dm0/vIuwI3S6hu5S8nzwMCHsZvOFORHz5FKsxgS636N9KiOebHAg3nXSiEmwnSbBrQ9T1JQLW6n2YKD5NalKWrF8V6b1VB7AQ0I0TwnYjSW2CfDdoTE3lgj2bRwcWvLDJrY9GMreSBVCDhX0mv7ekVospn/UA45RExUlR1mMMs+sbr0XEOSUDpTRPq/c+tFsBFDeSIqxlX4OmNJwCRPErkNc4hkL42YSramPTL4MbxY5IGnU5QeXUG4p2M13SAEOTpmY2vVhdZpzJpKeppr6rxG8ta2jkF61WGB6c4OrVDpdvEMY7QIpSx6O7BEBAf655CWKUbNUbzSJ1PiA8lVqlXscOQ0J80uPff/QJ/LYHP4XvvAr4Hf/rv4DtD94HCPgLb93HX331K3Bvs8PFfoV33r6HzY+v8eDHEzZv7xAfasUvADRM6C4mTBvhULjTEPGHTxAeP83p7aZ7G6S1WFponJByIBzN8khk5yJLQmyf02A/sMa4XKnatd0iXmwQV33J6G7XyWn6pzJ3LaBnHX9hLVIAgnJbap7NEoHjBjJZeUwKME7B1pytf2YgJlCcBDC8M1arcng+JJtif8ZxFI1/xNIDkm26WnqoEo/6vIJW59Nl8M7XnihnQM6inbHRRqgyg+BEvqovC2i9BHJsQOOkHsu45BOkNnZwHkexhOQK1AlgJ+FsVuCzE4ynnTgwdQIS3AGJGGkgsF0+KI/RBcR9QthNiE8kNB0hACcSoUkkp/rp59f4Q3/rn8R/eHeLi5+8hwefIdx9yNjfIWwfrvH0zRWeEtCdE17+EuP07Qknb23RPbwEXW51qgi43GL9RamQRsOEcLmVAj4XEstBqx60PUO3PQWve/GVePK0FPkhyTZOY7EyiANakiJL6ohGfacRn4WvMDdv2g+gi8tSqVzD/Zc8b6tERD6uw3NkC/NujlFVicvn2aR5zegBY6TmNOkhmObFhaxfQdQ+bsyy7wPjltvtAArGYZDwAOHyV1YeiQYWDUhwCOoRaTongEQ5CS8CC6MdA2gKIlX0KLVQl/gK79Xpbepq35bNvCwR1RmXYpEiEgM8FmJOYzos0c3s2TcbiUhNQH/BSO8BcSeh3CCxGGZX3xCAThZdvBgQtoOkSNPMTbTqhRNZr0AT4+7nErrLNcKwxqufH3HyU0+AlKRM4dM1dvcEhOKWsXk8Yf2u1DXFlLRmhhbeefIU9OQprCSgRGUO2RuSpwmBGdjtAApI4yggMexhDnPZA9OsO4CE32tMBzZrCUVPSaSjq606YQ0Qd/cJPJJEnx4otVA4LUBc4mNNHFqimiVR3jZtRE6snEHtGFAY6JmKYJJrZdFQ7oxCDRgmjS75DNlr/hDKz/lisHE7gMKaH9wFwnKmaniQyPU/QhFXg7vOUgukFaZYrgWANb2/mNAWhscPuJ02phIEcTeuPOSgZqospbjny45fUzGJWYtR7P/eS1DVFYoBPE7oH29BibF6GjGeBowbwriWEOr+SvIUWJQkDZKNmUbZ0NnCo5mfrHjy+r0Bmy8B/eMt4luPkN57KB978gDd0wcY768x9eqfsZsQLgdxjvLxDbsd0m6nJJ73CUmFdARKESM/XkTCT1iV+EmLN+saMMkzW0KU2yGglAKcArL4ZdyF54x8I9Xr0Rf3actqluRAkRBtp1a0kgcgkmdEqYmbnGWiVTssTaLyCYS4fMAsAUaw9P7BrR8HEt76V/3/MwUobgISOVWeEpbZp0A2XZWd2103AwBQUoP5pp6NQCf6nfN1ONpsUQQqjjvkMnB7O3aOhOWi8mSwoPwsNE1AHwsX40+HrMsm0NUOAcBqO6LvI9Kmw7TpMJ7GTGp2FyPCfpKNbORdDEXS6gJSL4FjaSULK24ndOd7hMcXkqR2mnIEa7i8wmqzLmMOFEefKcmC3g9IV9uSjk5JW/ms/J+jYHVcrEIWRcoSF6l1B7yX7FMjFzJT52aW6PZQhCczQOW+VRQwmfVpAqagEoWtiVjWj0b8ztLe2z0L4QVkMrs4aVmfy3oojleMKau5B6NLPWD4PaHriL2U4Xmvpdyuz9Get6Tgvw3gnwGwB/CjAH4DMz/SQkGfAfCD+vXvYubf+lw9a9SNyuxpA5VR1QGE6ZUkdSyZAASuHFuQ3OcAjdLU68CRZ74//lSwxRqAbGMnBQzrq9n2/aK01yYHeiSkHScGkX7Pntt0Uls45oI8jqI2bEncnCFeh3HVoztdg9fSp3A5iFu1eX52UYroKPZwDBjPeoynElUad0msD4OSupuNEP8aw8H7vVhi7LlUFcp+L0BRGUJUb8tVATvLIGUOUjbVPr7CUgUG8UVBjGpqnurNY+OgPhjmzpw3o3dWUlXGYn+qxikLHyDlrrDKaQur1w9ZEyop05lQbW1oBTCPMZlPsDgVSq60oF6jPagaqxrYkeTMEhMklJj226lbPw0SxX+CeUnB7wDwe5h5JKJ/E8DvgVQKA4AfZeZveO4eHeAk8oCQAwff9IThAN1UqQwSy/Tb5zJImDqQ3ILSOfLXbe9Tumo6oQIMAO46WeBWj2MY4Cs18TRJNaeO8+bN/dQixNxuCDMJ6mbMT55SYe7TBHQd4nYt/AMgBOF2l9+jLlbjxlFS4FtEZ9zL87H1JQZxXnJp33i/B2nNjVxvlRsy2kDi5ERNmz0sBwdZf52vROUp69bAsYjHkokqVHySkcZiQdDLNWtlCSzAVHJQAGB2lbl8CzTv13N4QJbMas5aAczBAjiwBmufHeslx6iSjJNEnYryvO25Sgoy83/n/v0uAP+HF+rFUmsraJkoqKdsFQWaWG3eLIiaNNmILyfoF6IHCK+KeADJ/xdkX8yL6U2exm2EACYtQKubGUAxyeqmpw6FV9H7SZ2GVKIDeRJzqfPOzByMkXgWRWnOQ+p+zJaXQB26+HSDtO5dQmKLag05H2MGTasm30nm56wL7waE8y3o8ioX2aFpAkeX10HvR+sVcLIpxXsBcbHWgsBIXD1LDpm2cVYAoeBym+b50yxQPiDMHNemCcTN0na+C3kjmqRGrmCwXV99F9irWTdp3qIikzAzl7b9qpoLRwe01MM10gVFBRigqLNkatICsfkc7f3gKH4jpLK5tU8T0fcCeALgX2Hmv3rjK3lpoonTyHr8YoUuILvY6nUoof5MCxT228d0ePLJTvZDLuQken6Vzt8iERUwjBfx32dWtcZY8hgkPFzzZ0I9KcVFmTUyUMVubzVBASE73bOjjon6MUpeiLMTpPtnGB9sRHK4GhEu9wi7Ed2l5JHkSIg7DcM2YIwEXnWYTjuMJ50U7tmu0K86BK2AbhyF9G0oliJHoombuVl5OtDUyVRpMRsAoN0kZQZ2GkCnvi3Ud/LsLh9Hll4SF8DNfE5fwK5NTzclTVoMZHLSrZGsUlVWhrFYGWz9HWp+jRyKOl38Xi2lZhWHrXj2AVVEv2uRxwYWxpHJBc2q8mLthYCCiP4vkPod/5m+9CaATzLzu0T0jQD+SyL6OmZ+svDdUlIQpzOVI0sM3jHJ7MOeSMthwsCiz0PbWkLJvTazr7eng7VAqAiMNizYQCZzEiUbFzX6I0KQJCuajQnDJBtl32RuolIcJvdZ/S0M0DigBDr1kp+Czk4xvXIX21c32N+LCBNj/R5hdbkHXe3RDRPCVQ90Mo7kkrwAQOoCxpMOw52AFAldTwhTD9qv5PPKEWTRdw/1nByBOICGXvqUohxukcBdLyrPSYdpHWA1Q+OTiBDVx8MAz6pm+Rwlvg5r9k/pa8c6ewZmQHN+0F7yWghQU5EqqvWRwKOqRGa1sYPpUIoC35bC0J8hu1QlcfiyEHQELABY3pUcZxSCHEYqXT9TGYKF9txAQUTfBiE5fwnr02kV853+/T1E9KMAvgbA32q/70sK3qOXGcAyeemZW8/g2oAt+DocTIybP9AuDp5/1uU/rJr50mtOAbLPGGobN7EftBhNURuqMPS8+KR2Q5YoEjkOQghDcCjjYM9np54X0wFYYRyRJAQkrt44weXrEcMZobtkdJemRuyBy4TYd5L418KfLSU9ILkxO8JwIklwQAFxGxFXnZCeztpjwJEyiavgGAghsfpuUC5os7/fYXc3AAT0FwGbLqDvAuhK66w2FioaRQ3DdlsqmltWK6u2FqNyRpBkNcOYHcpyqYZJ9XqvgrTrwUhDTLP1V/EWS8mOltbTsbZkcgfm/QJqydj313MWaj1joBDkfz/Mo0T0yyDk5T/OzJfu9dcAvMfMExF9BaSk4I/d7KKhJi+NlwhuU2lbnIglCWBBFLxJm5mxrIVQ1L5WpNRFzQYSGqfhQcKnOsvfMU5iVGlJK535sGci9zxJTTnZulMsD9l03HeSqv7eKXavbHD5esTVa4TUATSV+/B2J4l1xl4T87p8EINaVjaSVyL1wLQGAELcRXRbrbYeY66kTUAOYuNRHJ9omhCSVOeiXqpfSWbyTgDoDmHaAONJAMceqSP055IBC8r3pE7AhBJLCP3FStSeSUha3qwkR6bGjxjnE7TEHgCVfJwk4LiAnJpuySwJAG38jvdRsPk9oKIeCvWuCNYjMUU3DhVnU72dFY4oW9E+8OS6B0oK/h4AawDfoehqZtBfDODfIKIR8vi/lZnfu0lHKl4iqx2NTmjBOl73bEijnFTE3rc2I40a8XHB/2L2fkpOpBOkZn/tyaI+hzlIeMAzaWA/5NOJhpjvUXJHyE9J+455gJCdaH0v1+o7SU57usF4d43dg4j9fcJ4CoQRmhJvksxRmu8BgIjuVq1sq1mlKCCsVwjDCZgippVowMMpIe4lmIsGKZIb9pMkhtnvlVfQmiO7nWSYmiYBIlUTuhDQ3+9BrxCmlUhUlABKct1uKyoJR5K8n1HSB3abDl0fES57AdSoqpvV5uxIRW1GClLkp/hsqIu8Xzfm68FHzJI3bW0eiCMh4CWv682vK38f71ulgnCT4+IF2vOWFPzjBz775wD8uWfuBWHOSxzauOx8EzwwOFMdgOMkpLveTIz0bSm9meaUoAjZuP5+Lk7DnIjklg7wbKEqx4BhyCn9M9+iOjWmJBvZBYOBeeZdmKUUy3q16pE2K4ynsrlTBMIA9BfA5lFC/3ALPr8siWm9uVn9JdLVFkiMuFmju7iLuO9LsG2QzYsEhD6IB2gg0LZDMI9JHTvO3pJqFp4SeLdDGEds1hHjyR1wiEgrIEXCtGKEjY5bkvT0aUVIkUBJvE5p0lwZVueiC+I01gf1bZF5jaN22GVMqworAcXcSQEHj3bjz8xzc4mnsPR2iefXWTp0gKNAstheBMB+RuTMBB1OasqpHvdD4OCjT/138y1s0xr3ELJINgMlBxBVYhFbIAYWuVJwIS+zi/HMz0Oew081myrSDcBeU7MR1SnXAeSIwIXxyZYh85PQamQmcoeB0Z8D3RVh/VhS44dH55JiDpAw7pMTccQCgGkoBOw0ga+uJLXe+Vo2y8RYXSR05xPikIQHYJQKVa7P0KQztFqBTk8lMQ8g2bUuLhE+9w7uJCDuzrC7F7IvRxg5Uz6U5BmscHPYa4r6MgBiCg+Aep6JCdRMoaNLYDyMJSwdCuDOgWoxtYADiSqEwIO/mUDt+c2N36Tfllc4BBxte0F14f1stwQo5k02XktKHgGHJfHs0AS0iW6O9cHbwCGSQD7ZSYjGbE5tTyv73pK5DiikWhxBY1fK2xnI2KKcX7DcR6NRTXUopOKE7mLEOhLiIM+4ejKhe/cKfLWVZziTkn98KtmyaT8US4qSkTyMCJc7rB+u0V0GdNsJ8XyPcLkvJku753ZfZ50CBCjOzoCX7mG6uwGIEJ5uQO8+BJ9fIHwBONuPWL18iulEg+DUnyObaUnUDhoTwn6UpMRaNwQhAKtOrCoJ4I4AFg4maIk9Cyn3QC7DqP4Kal1Y5AMMJHoXoeyzpdlaYHEmwzRlb9EcJLgEFvVCW379WVrrqPg+t9sHFGZW9M3E2EPgcAyhD2XBss85lcJfi10/qmxUKcyzFxlwOTNoeZ4j+mVm18tnpK9diTD1WbnMCcskDii3YCCh16RhBK4IHQCaGP2FvBcvxKWbAdDpCXB6gnRnI6HgO8mtyVdXQki6saXtDv2XLrFKknyGzy/ErVvZdOo7gIIEAe/3JdgLEOLydIPxwSn2L2+ErLyzwgoAfek98PkFKDH63V7cz/vGwYm1CO8o6e8k1f5Y1kjXZasNrXvJvB0ol2vMIKHSRE1YGi9hOlXDT6izV5YkVK2zuqieZKcpCbE6jhKnsh/qIMFDjlMv2rJDYWhevoZze8Z2a4CClwACqLiGgyfzEjDo30tp+qss3Pl7lBd3tendJs4nkIJFpe8uSTZtcybOOcNe7OXU91Kr1MoBmuphGb2HUeJFgBJebwy3EqGW5KUbplIxaytqhXEi6c4GaSMOSrQfJFflxZUEdLlx5e0OYVRi8vIKabur+7uPWXUsxZn0/RCknOFJh+FM3cW7DnF7ik7jUHi7BXFC2J2ImVPTBcj8J4ln2UnaO86p75J2UQnc3QZxvQKve5HKgGwaZc9d+ZYtBYfmyyxxoVRis/GOTdY083SlAh7QmqEMiNnyRcnSqm/H1/zsMy/Ybg1QYGGTLbLHHhzcxjtWv2OWSaq1pLASUJaJOzEWo/jMVh2SO4SOSDa5nwukJrdRjNpfK+SzXkntzD5mFj8Ha9l1gQISlg9jHLOOToB6I5pIr53uOyk5uOnBMSDspXJWjvrU/uTFpwCRk8oEAqjEb8ijLwC9vkZT4TLAQoamLojaFAjYa+1VQNQNy0hlfd4PEsm622vsTCGQLd8pWVHjzUYlHMrE8FHG/9CmJVfmzydJUrKdc85NI6kp1wj117X/31ewaKSIg/VoffugrR4/Lc30e+C4qA5Ug1TlzLRITE80WVxILKBRGGflOhK5vAmptnQsnkJpTngdWWyLqg+AUjFKRddYWHVbgMbmG+kaAOFt1Bmq8uUwR68sFak33og6iQqR1ELdaH0MAyAra8gMSUgc87VZ63PmrFIryWcp/hJjGbMmw1LO1LXdIV4M6M870BQQd2JOzSrlpN6oUYHax2kYAJnJuS2tyM6aYW7alvrfpCwY6RsALmnzDjYDd+Mn8qFE7rBppNJAMpd6n5wIKWisTUq1k9fzgkW7/vVv68NBkHjBkI/bARTAwundtNbF22VR9mnOfLLd4ovhfoAMEjklno8Tae95jK2+QVsCM9hp7ytGNS7g1dAQ5DQmY/cLiZZ9SUy89t6f/nn1edgqm686cEegvZKWZkmwk9TunXkRBxKbtVx7HAsA+3HnBOwlpNsqj8fHG6wioVuLR2d8uhV/jWHUhVzGgbuYE9JIQqBiWs4V36keL54mYL9XFWAStc0TwTGKr4l9ZVFiROEcTJoA5p6TwY17cK8vHmwutT+zWlmaOJNnXVf+2Z2qcbBcpjzws92jabcIKK6RHsyG7cXAGDFLx285DL0UARQxWTdXThXXiss+C7cP9T3Ux0PNn0oeuJrP5NL10E1hPMMwgrqAQISEWBLuJNQAkf9OsqECHGjO74cYxYNxJR6PcXTP1Jr8gLl6lDet1V5hIEsZUgtVXtvKht3vkXY7hCdPEVOS8PNxAq62SBeXIsXYuPerUqB4rV6ikYSTqRL4MECYzREnzpIVAUAv0lH2MQFKwNQ44mAod2s1MAtQUi/aRLJ+KJUyjqnMAyxhb2Wet9KQXPrtE9s8I1hU4Q44DhLidPVMl5+12wMU1g6pFj4FXpttu9UZqcmTaS0hk4JLRE9Jv07L7r3X9blZbBRrEJuVMoyxqFwmvu8HEVGdRBDUTCxFdkfZaLYYveu2Twfo+mTRlBwDeN0hrTtM66hOTAzugkSD6vMvksaOtOP9UAK/AuXcE1YXVEK5A8I0IennpIyhelKyVCHP1hUDGq3GPp2tkTYCDKELohpdbau8E2WNJJRwbMA8QrNkRATAgQVrCnzmwhmY+dKaemmShmpnXgYKNM6nhfwmz740apUyJ7w2JmdpbT2vdIFrJAnt6wfuwv3T1lzJ+UrX8put64ppqtOq16bPa3o3EGmNA3fthHJy2wnc+EjkZjxHlXHIEahHTbGUN5UnwireJOuV7PgE7YsLxuKg12AWRypAJJ2hmPyqupWqUmTgJMp5LdgSxHbi6jyddBhPQyZJuz6WSlzAYZ7IslcB2SoUohRL5tMN0t2N5CBNjBAkqS3tJXkOjwoWWbJzVpOuA/oVeLOSsgJnXfapiJEQ9ivESynIbDkvsnoUDZgaIPcSF/QQUGmLsoSg7s6HYi2Uj8qSiG38pVQHznR/rcfwkrdwkzbvKGAYkN2Ud3jByFHgtgCF2uJndUK9atGCQxfBnZ2ghYPwkebCshvVznnj5PoQmr+hmPII2YFqKeNQI0oefBaXNyJLQd7iwkn6BOQF1pK5eXFadmug5FMYragxiXcnWc6HlRCUJlFPCpBKsnEfc17N4URAIQwBadUhqsmv8k5s+BlODCQFCqun2ku6u3SyQtpI+LgEton7N23XooIMI8BDPVRW5KcXr1LhTiLSOmA8KUifrnqRWvoOtHP98xzCQoo6H+dgrt0VIRlI6AKfbrDcdV5nY5qWT25HKl8bQuCfP5A7hBpL3CF1xEt3FnMUD6zFtn8v0G4FUBARwnpdTFGed2hT8GsYcQYJAixUO1/PnyYTl+zTll3JZVPKE6qLJmc4PoTC19imK5Dou/w8FSGm3nrsPQWNE8ibIIAsP6Xp5xZ0Zj4Evs5DCJJTQit4S8EgvZaXJjYRg2bsJgamvQRkcSdjzlWikwO+HgcWHpNKQpGkH5YTY6/Fh31sjJ9j428aQpADcqHhrF41MRk+HSE76aHqs0XcWt8zn+A2c/tMXMaAJwjYAHPX/Hyb64GhdDo0/zry2MUILfIn+f8iVVjSohvlyXjOdiuAAoGA9XpGSlbSg24EjkpUBqqreAOQuptJRcQGICwi09KtaQ7IvEhuYmemcHgheE5F4xsyl+LESSMAc5Vyf6J4AnGagEGze2MooqzpwJm/CUDqs4SEJGqbxDjo8wJACEiriOlEAsW4g5pOrf+qw696Me2Zvp0YpRDRQtMsW7QfEXeaq9MkOdJqX32X/Tts8VekNCD+FsOIcDUiagg6R0IYuHbl916sjSmQyhuuf4ycUyJnxSrRudeGcWf+oZidF99/lubLFB4FDcefHAALiTVqvIWXpJ0XbLcDKIgkv6KRcZ578OqF/mYieE+WHACkQUAYU/FkNMnBKnF5gACWQcKnwpv1dRm1fbFhqzlBjYehbPapdkFeILfI9HiLNLXv5zqbSfphCUlikHuRuC4DyIBIo4jLKRK4D0idjB2NQLdjxG0SnwZmdVFeCQdgPhLDqHkkG7VLWX3WVHx0tUMIARg7IEAISBOTzf05sapyLggQkNeHAXRFCDGgD0AYVkgrCUKLW+U6vIu8WWCy6fka8dvG0lzupxvWuG3boc8e2PgHpTAvsaTlGI1r82YCmWeZgYVvP3M4iqB5I4vloqryFS1C0KkYChBkDkh2gnqAMCnCTuF2UXhxcYGAyu8d7Toh+3X0fR0P4FSGnGTX4iD0+sLLSPJY8rkNfdi3J8py3zWZipUGpDDnM+wZtMiPEJyMbiug2l0mrB5LpS9SSxBpFCmzuHWLv0QhIGdu6ok169ROVMhR406s6jggaoblvnTenGUOJqnboeMStbhwdj2/2os1yDtbhSKtzPwqfGPhpYCp8FJ5Lbz4Bmr9eyoivolRKmpmzf3Ae/x6wjsfHAeycrff9WABzPvwAu2WAAVEisjmPcrMvVky8kdNP/fqxcT1CerVCzNNHQo/9y7IS9GpcPqn77LzilsEiexgk+Yg4UOxAWR/CkrFkuHeJ3lgACGTar4oL7bIKkDJqWmWji5vUhoTuquEuE2Iu4TuwjJqb3Nx4Uwmm1XAxs9bgbwmYubIYRBzcEW2KWcQAzhFUGfSgI7/NGkEro7TgGzRoe1Oyh0SSX6O7VZJXNlMVXi9jbeOZembjhWl6n83iUW0P9YWJUs//0dC0AEhc6dJVDovFS5t+AjkFIgOMCpCvTXJtmAB+T69TyAB3BagABWQyNYLZ+JkjRNQEspLETP1YnRRgq4oDYCZilFlSF4AiPbEqYgsn7bPRGvNJ+GvmT0bzb06NWDgXKXRbjQDivy/+j2YS/M4FLDQccJ+X04StdwIiTshXg1SOWxICNu9RIJ6fwbrv3dSM4KZ1It1XDiFbUPuh7mebM+jAEaB1AJSpDabLxGjx3ItC+meJJlONX/en8biXZb6FQKIRH3yPhQzX4ZD7ZgFwmcuc6HoeRy1T6zeotkZzIfitxt+MtXTSRb5XgfAQvuZwQKoAEP+/xkiUZTNoS+xpbvn/H+2YizxDy1A+JwNWEBVH1behrEDNWm2BBALTmA+FyYzlxwILktVPi2AenF7Ui6b8PSaoWRWynkP7LkUJLJ51S8M+7vrQNu9FPJhzRux3UnyGg2yIiLQiqtiMj5OIrsgL7UkmbgIkI1oFg0DHTM1K7FLQFYjMkiYhOStE8bBWF/MdGzxMMZTKDB5cjun7A+6fgAwd+V1tcJYGrylli0QrUMUNRuYXH8swtQFjOWExVrhLY/LMGawz9dupIOqLYFF3WHHcxXAeD/a85YU/NcB/GYA7+jHfi8z/yV97/cA+E0QG9a/zMz/7TP1KAHmS0/utewwZU5JSypGAxAHk+Q+j917yRHMxF9npqzASkXOzPgbCbl0eSc1wALaYixchy4wYgb2Bwgzew47nZmAQbNnEYGuAB5HpN2+BFiZpYVCThLFbKnmTNeVfJZs42DVt+32WtiI9TsUJwcYnDdG/gxrVTfdQLmIc34czl6RZcociemCoHTwShEl+7wdOAk5GCznkNR+UL5unJ+4bG78Diz0XgebXx9dHYZu/BENYuKn3a5IHBbDomNjYAFMOGSOPWiBy5YzG48j13iG9rwlBQHg32HmP+hfIKKfDeDXAPg6AB8D8N8T0dcwH/R9K830WWbkLNP2ugY9ZUcpkyBaDsKTktZagGj1tgWAqIroAjMpogpd974M1hcjT41LcGx9Pili0c0rUd/EfO+i7nMf2Olf9b0W441wzAlctbI4M88AIn+fJ6S9CnR6D++XcMgCVHk32veSej/2uiFjqOdhyUphG8pUORtr649FX0aUuZj155CKQBksxI3cgYWNnfPDgAOHyilq9vhOgmn7k6UMJeFVHaTgauZGVRlHkjB1z5/MNvw1zT83UbOeX1yqeK6SgkfarwDwZ7S+x48T0Y8A+CYAf/34TSAT2BJR+kOpliZm/vTXSRE3UC8O2tM9F2Fel607tvfF8ODgdO9WmrCTbXavLKXE4pbtJIrMcUwdsErFlMpJT/ahgITxF5wyOLTlA3L/kpxqPFItXsdYFdSt+lqN3SQypLqdZ4kwMphDGS8vPTTekZUfio2Bji9PCRSmGljb5g6XWblIA+oYpD+c4L1wkVIBxNYyc2hdWNNYisVzO6saADiAQ5L1ZGBHJKbf/V6At023wAtjfp215ghoPG97EY7itxPRr4cU9/l2Zn4I4OOQWqTWPqevXdO4WoRZjzSS0aIUWzPnkt+8tQYklnztr3W28ZLEIZAw8dzZ6S13ZlW4V26IykffLdK8KSgUFt87n3liU0VbwipbV4y7MNG+pMw293Dnzaiu0xInodLZAFg6ewtYQ6ijW+0ZZs0x+FYUie35liwB7CJP4U7m1g8lFnWOojigzbgYe173dz5YqkOhmWsf/u3n4SatAgmVhI6oJIWYpyyNskoWFGNJsmxE55QwS2yk41Zf+Jr16z9znWXnmva8QPFHAPw+yHr4fQD+EKQG6VJvFp+mKinY3S2Zle0UyC7YSnK1AAHcDFkX/BCOutu2Ti+OkKtaq+o4i8YsjV5LgKmPfrZ5+xBqa41onpPP2IaMSjpGAQqaNDN4dsJaSNBCAdR34i5v+RqmSUyjSbNFaf8kGEqC4+zZS76LuYhc/ucabGxsFpyA8ryYidmsRxaybmSgqZpEIPPpyKDjxsUTjv6ASQ7IpiWga8DQk+BLsRhA8Z9x1o3qel4iZgZDPxchIGG8Sixp9rAf5LdZyPyBY2N9E3BYas/7PW3PBRTM/Jb9TUR/DMBf1H8/B+AT7qNfBuALB66RSwreX7/ByB6TBSCqDe4BAjgOEkvmziVnoVZEs+sukGXl2iJmX1sy4EUm1fePuYTFL4AFoBIYac1JE9mnJKQjuUXdd5Kef7POVcGwH+Sz/VT6nu+fABeGW+dTOPJsRsiZIHOMTDPpzMyLllFcq4oBChQxVITptZ6Y2o+KSLa+ea9XBxBLaql/1sVcrN7yFZqcGRkIGUgCEsjmfwaYQEErpk9SSS0nAx4132flfwHgg0jQe4P2vCUFP8rMb+q/3wrg7+nf/xWAP0VEfxhCZn41gL957QUZxYuS3QkNzAEid2Kur802b8sR2PUW+8Dzk8G/l1IRjQ6Bg+/PkcnkxEX9CHMxFkBWaWhY6JeZAm3RIdWJSYxsdSw8dZ3E06xXJSel6dXjpKa7qDkamlPU8RViKXBepCYizwLIGDk+osX0yllJAIjM/8CC//Q3oCcwkDNUMbS62dJp3r5m/IuLq3nWolGLOSkPRTl7PimPg0oVdnnnZ8IRoC4AkwTOoe+KRU8LREGLQZlasujh+QG35y0p+M1E9A2QLf4TAP5FAGDm7yOiPwvg+yFVzn/bjSweaJB/SXo45KILFPE/OzlNz3e6t2CRuHKJbXmOoxLKkWtTq9eGMNe5vdrgxWJzNEIUD1bTyadm3HRTS8FiPaU3a0kRp2nmsqRisRdEOa9jbuYHYcRjBgvOxHKWAG861u042/NTyFaf6gdJ81wEwBym2uYdnIJUFcstqerauOcvSRBeUppJEJ7QNt6lzbCWfTwawGKVKoKChk+NAOOZolRdG9VBaxhB46pIGfuu4jF+OgHjfS0pqJ///QB+/zP1giEDASwTlK0qYK09IZy6UhGJzzqQ3sttgoRe+4X0vNf1zak1LRFW9OhmIajVgyCuzWR9MEJydGZPO+36VZEi1n2pm8ElHWDJ7BX0em68rX99Bx8uT0Dl7IZBVJgsHh8aHy9NmFcjNQSxP43996j0hw1sWo7CfoxXqqS0Rp1FDQyLFgbPnyxlWDPTqzl9eWmimdfsROhUkZImIeTPSIBjBxp7QAsY0TBK/Rbzi9nv1az6DAD9Au12eGY2k1e1YyBh7LkHiRcduLwAvXvvC4DOUvPefKEsSLl+bdHx5FsOybYTy6QIA4k2V4WCBG8UKHqtEcJcAuimtEzwwUk8JlprnREfEYv9AAxBHZjMS9SukOpnw/IpXdX0BGr9PqDmrnzzhGGqk9OwSklkqQHNmfVZ+Cpgzp9YbVUFjwwMcUGKaK+rqh6zSKkyF4RkvhaEAhiJgZ5BqwgaJCGx1HQp96f9UEIDlgj097HdDqAA6mAdv5j8/9aWQKLlI9qBOjaBbXNgce2AHzM7LXy3yu681LJLc8qWFAAFHFiTvKoENgs4Y85+GGTZrPuuAokcYWsibSr6+6wF29BBr9MV1WVKstHNRAsAvBdrCSfA54hv5rKuvaIqg3EmKYEmlXIMawwUjcPyaqH54KQ036wxyrPGAJ5UYqqC2o7Mbw74ioVktfwaGsRoiZREbSv3JQdi2bPYWoKoIBODqbaKiEoC+Z8BmgKoZ9AQQH1E2IsbuJGgVgbSJLqKM7ru+Z6h3RKgcA/TnDjytztpbgIS/jvWbppqLH/mwHtUg9iSe+zBimHHeBY1meZnM3u6fRW6kZJukqRSjhX9MVAJzmTXyYJm060VJGgn+R2w27vKWw6UgDmYmQ7ed+BVJ6ZKZlm0zKVAsZbVy27RS96LPmrWsla1c0ypLqbDnN3is/Q5ycox0CRmcX4KVHxx7Fn8PW7aTC1yQX8GlpKf1Ec6uzFLLB1jltgkSmYIkrHyKsgk7upMZayYqLigawJkCgTqEhA1g1gmrAMwmIQp4QJEjBvVPH2GdjuAon2Oa0DCXLoPgoRdI/9ZFsjMd/+mTfXopWI+1aOwz+481d9v+gWgdrYy9Sk73TjR3VtHLC6gDTiDcgeBSuRlCOoVqItyUJDY7gQkNCHMbDEllpPSmw/VFZk72Yx5M4xdFomrtPVtLkf/DFn1oryR6xR2XAOXucQz1xIQMzhE5WuiuGfbLuMGLMzBydzOjxTiMc9JMiuG85Rlk9D6ULKtEQr3wwXUiVj265iAIACRwWICCEneNzNyUGlC+8wAOBIoMGgS4AqOAyFApAsKoGGf14U8/wGL1HO02wEUtj5ar8cDILHoIg3MT/lD2YaeqW8OIHwMQqgXthGqNGmW5kq+bfrQOFtZ3QUPErnIryP8Sgg7z13F8/hRcQG3jZgAQEFitxdJYr/PafcXTYT6TDn/aFaBUEg7PcHFYciIvlRbRapnt3ly6sfiGAnjAd+nLFVM+fnlswRzFhNxX53E7CR1Km1bjgG+EI/vmwWe2fgZSJjp1kCiE6AwdUEvIJrrxLDgOYJPKFPG10ym9llJgQgxPWtfc+Y9Bd2kfi3BkbgElNgVCiJVYIB4yfoyhsvDfZN2O4ACeHaQ4CJBzAK47Hq+tez2TRDWmcJmhYea0y+ncjcRublO1TxY6IZivc4s3V1FfNJBsLTPUU7mW+p7EHMucIztDmxAcSiq1QA1Bdl0GqFrNUSreJqAYkLtuuKcZGIvHGjYhlQOSKSvcnrWJtGmGTeh814sXXJ/4pDHPpO+dm/nQ5GtJm0hnvaeXu0w0tkyv3cKEn1QKatkf7f6qtZnixuhANn1QTkSkwAngJGAiRAoCRCooxwACVYjBYwoUhx3hJSigIUJnXBgYb+nqal5+vztlgCFA4mZ+NmAhCOzytcbZxgAB9OA3USqMCnCgMGyV/nQb9dKFOv1166iSBEzrwBglpG7SoripIkSU+LUk0BiPuuL01I5VVMNErtdriVahnCh75zALLpvrmA2TqAxIlDZHKysP1kcAzpRDzTVPSXMM5svWRk8uVlt3BrUZ0TvpGHtPsTdf98kT2850WC0edUuV9MzR3nGnKaRO6mwxtGBhG5meVwrF6nPZdyFy63SjgFNAKCb2nY8QtmdXoUjG2+R6qhjIEUBctYoW04CnPCSBh8Ucm/SbglQYBkkgAISmUFueIglFj1fUxfnzHvygDThbfxqEsuxB94s1kg7NI5YdGlekiTyn1qJSnmAKqDMns2z7kv+AHY9i7g0Zt7yj5KeZsOoFcF3AhL7fc1LECEneM39M85DnbLU/ColA5JKW5RPtGwi5A5Ek5B15hWJCVUmbR+EZXp99n9o5jGfyvZdzteoHaZKiDumhLYq26IqdISrIlP1jAR1KRph4OBAAk5KqZ71UMt9l+cjmPTmPxTEcGQ7PnfO3qbcJ4pKHndmRQEIowOLF0AJ3BagMPS15nVTP8kNKz8Dh4rwS2Ux3dRs6rmIqH4D5rDkSxfax6ckk+v1XeC4ZGGm1/yv9seDBOSkKx5/lBd1GQvOQGnlAWi9AtYrYeRtLNTDLwNECxLHmvmSTGqCNJNqjLCwaRiPYXPRyZiSkc2AboiUgd7XoajCs/OGrIFCTmg3pg5Yc8uOcQoY7aFhY78QJLg8Ta5f3l3eAMKBxfIF/D2dJJMOvIZU6uH6a7AAss88n0tEGphFAqbiO5K9V1UlZlNJX6DdDqBYaiYF2MZQ4gwplYfWhbVEKmYV5ZADVgMQ8ovK6WwbzwKU+q4sGC4Altevt8b46x5zgPEmUe86nD0ALaVcKGNiC92Asu+lFsfKgYQvQaimVuwH8CELx7HGYpu34kMZLILo2zk5LpAtAEQmas+fe+YDoVJPNS9+kzvfiExGVt3zHMsEYqoAA0DtILb0/eZ5q8I6/j1yuVyVN5hVplPTKBlvxfY6Cqj6HxsjvVAGi6SqDJNIZ17iNouKv4bnd2ysjGMB6qQ4z9FuH1A0vERuNkje4nCIsGRP9rWh1kW9yC+Zec44CH86a2xELnqstvAqdZkDpmymDAtp1Baft5F2QqxVDlOvsrrhNkAvhX1psyku2gqiZOA6jOBhUPJyqBdWGYA68AnCWZh11woSkdnpuwgMVIAM0LlQroVITHnJTne36fPcOjJXD4RqNpUELNLifK5rT0TloyYAxAUwAl1fRcuDubcspTRfh/p82Rxq3fJLzEhG5pyVTfJmslvfjBlgm1ftCEn4A0g6wyCJblqHLkkyjfm1MmAEIU+lhNvh579Bu11AcUgk9ElbgLmFxEsfldlwASSo1sUrn4j2dDaQ6LwYmzJrbcFYuaiQnVg3SKMGuNPMQIKoVjlM7XB6eR6PGEGbDWizBp+swRsFiQBZQHvJKyog4eqVunEoQ+B5HRvTYOxcHmPJrKWEZghl4dpp15n/QgDTpNaWVp7OD1+sKqkRxW1OgVoNPTiQzRw7qaBqmh9D3KflGWd1ZdnUJCqxLC6fhR8nAYv65Wz1sFM/6WY2sHA/lVMYWCQyAwuoOstGqDZ+Owo8lJICxgL4GLejSXpepN0uoFhoB8HBk13O5Tmz+W0wVXtqNiSoFZPxIMGrXmt56mdTypOeIyf9T7Vgr5ckyt/WT1UzvJXDml9UMaokoVXET9fgdSweiWnKJlFoJa9CkjpvyNZBzVmN5LR3BYGBsmHs+ZWjAIUicRk5CQgT707yii+CSi1mVcmh4FHUFhtzZ+FYLOI0G1cnFZgbeajH1ACPMc35Jb2G1Xw1j1MaJ43qFItCdkjzjYQ2BLwawvW6sUzoLTDqs7ZgwSkJaayqSV6LOg9WugIp1en/PNH7PrTbAxRLm6zd2P4ktAVjZsJDobfO1DnzplRduEqY4iSJGUhM7kQdpyJNTMVhSURtzCfJdPklgk3fr0rtLUlL9rlOEtDwyRrpdI100oF7HZdBXZ8tK5Q5VbVgZE5HDix8TQ9LZpuLky6QdlUgVkCWaPL7Yyhi8AF+gIjLPCaRgswVOy/4caqD3kwlONSYoS6LKMllQvaAtUQ/IMqE6GJSGJNKx1ELTU2g1Mk66FTggqpaChDstJFqrCo+ouEnZPDr/5UAJhKQMG4kXzurwKlwUS75U2UVex5Hw6bdHqBomycsgbmqob9nIOEL6trpaadJYyHx6fYzSCgZyNEht4l5HiQssciC/0QBiyN1GJZOxdZVGyg6sm40C1DCqgdveqSTDslJEwTIglZuQsamkSYOMfV+vINZNGrJawZ02awZSuxDUv1ZU73l67XPZQeCSguW85O4Qz7l7SBwRZQWM4Ifan6erYqaxqCQukhnmsGDBbtKaUriWiEl6iVvBEWWcT8GEp5DMElqyUxvB40KajCeQ4dNwB/155lr1SgDx1Ttkfej3T6g8KrFIS9LHwx1CCRCrEq9VdYRu14TPp1ddF2kJlhDslUE9iDBZjZ0Xn/XJmldBIhGrQJ0MafMDUi/jWwVtYjXPdIqIqk0QROXGqzmO1FxE0W9uMlmE+Y8ZN+GbIGx3AvqECYV5gmps8/qJts3UmB18VolkUQtKl5PU3lfVZIqpsUsWX782ufx0auBZJ59WP8UVBJAAYucWsA2XipzvB80VV2HMEwaf6Fg4x8vEEBcQGNBs6n72UgU7esGFpVK5TiJ5KQJ/Z2lCfud+NnAdaHdHqBYMFe2qkK2f5veqqdNBRJO1bBM02WBe7UjZh+EKi2+3W/iItqNU9GV7XSzcgGtWGecR0rgA554i4/fgCEbSw4UaYJITmjtb7LgLBSQCPtJAGInLtrsQKyypTdOaOKZ6fUGfQ6gcDg2ltmDMuSfHPMQCUgE7jgnkc3zuGDezH2hUi0sW1IqL9SpqAlLXrlLgU9VUF2RHvN4EFXei8KtuNgIu4ZJFQq+1AWEoM86AKkP2TOzaur9maWuqm8LI+GliuhBQa/l1RJPiLaWFP/aIWezZ2y3ByiAgwABOJCwE2YcC3HZgoR6VFoauPbEzl6APjLQ8QLZtDhpRTJDays65PpxCK3ttYMh5/mZKfMGVd3SpZYXexQLgweJiaWu6H6UkoGDK/Tjnt9LE3XfNDiNuah9Ftmqkgycd6plkDZpgi3+QeMRKLGqIqH2fyBTx0ofZgWEJidt2Bi3Kh61nqROEphNRipz7kLcATgvRv0oMIuNYDMN7wcxDe8igoKjfSp7aC5xU0utlYJa/spLGgYW/r32ByjShHvu96vdEqA4sJGt3RAkfMFgIyiroJ58OwOKchoCQHapVaAwwrINaV+qbQpgrnYsWTaqx6ay4MPCIvFNLRWWsamQrPLZsJsQdgoSu31xrrJbtZ55FlOiC7RK+Ouex+v3WfqySNFKmhCQSH3QQsoCIDPX7IVsWgYW2btyqsdw5jCnxO/M76NSG1T9cfE+GQxiUCBz4+1O5wIWngubJO/FfpD1FANCJCTqNAlXKEvWKAhVbSpfkGfZuy1YNK0yr7YAYUDMZT2/SHve2qP/OYCv1Y88APCImb9BK4p9BsAP6nvfxcy/9UY98VJEAxJZ53LmzyWQEBNnr5yDLmwDhPJA+XfrKgxmLJYtNJBqHMFmkkRlPrTVcmCSlky2do3Z2NiGDbV6lBhhkHvFqwF0tQeutoAGfRlQeYtP9oz05jmzEiTOan12aKMgY+l5nBx2HdyPgARHCME3BkSnmlR8U0XK6Ub2aQcxf39mFWvNuT6kuniKZY6hNUdyICFqY5B0+EFyWZg5c5Zk2FTdYcg5KiwvROIOoUNFgBOjymqVQeNYc58XUpPq1z1gtL+X2k+zefQ/QVN7lJn/OfubiP4QgMfu8z/KzN/wrB05VEgl+/SbFKG/i/5ZQIJWqxLAZQVkLDKzJUQbMKrKFvq6pk56uCkh1PoLLDxsDRKByvcOjYt9Vv8nZoRRgWCYQFd70FWJ56j086W8HNeJpR4kch6GUg+1VjcCUh+QOgJ3kNDpTsjN4M29h8Rw69Ox5qxfdRyMEbkupNpq17oDRvikGiyKZGlRvAqm2e3ZdU/XIY0kYEFissyqB0dwCsWl+9BaURKYqiA5LmPuvpstINZXBw4V8F0HCD8dHAUfqT1KAue/GsA/8UK9IFwPEo7ArDwZPUis1yUBrAsNbq/rf88Y4xwhOc6JoHaRtc1JE0eD0NpEOHAg4Sd9wYO0mIgTaNDfuxG03YGvriRCdErZfdzUmpk00fZpqfkqVjHWIKHj60Ei9SR5YwJkQymYCAFaLE85sc+xOBg3VvL4Dih8VuzKSubISUtJYA5TwwjuRpBPFejuQ6GWKsCWCLcW782Ma+oPQ40SE0tCG+eERRPPwcm3mUTKldMagBxL5C0riyrH0j1uGsl6g/aiHMUvAvAWM/+we+3TRPS9AJ4A+FeY+a8ufbEqKRjv1G8ugcQgpr7slu1BYrMuUZ4WwGX8gyOtWlY4WzLGwn1UoNTks5g5a8mLuc+19+ARkPDOX57XaIlPq3puCW4d2UljAiBmUNoP4MttsXQoN9GCkfTTE4L1KT0DJvUSRQwumayAhBGV3JH8REJS8ZsSEDoBCZ82fymo6zqQqJ7BWVtyfhDngEejgm4eS10rFqOyH8oze7XTN+unt145shBEEj8DPd/MtT11jvy1+SvrzjtcVRvd+1QQFbBoh2Pieq0dGzt9jqzFHZPkbtheFCh+LYA/7f5/E8AnmfldIvpGAP8lEX0dMz9pv8i+pODqdafMPQdIrNeg1Qq81tgMlxm5XFNEYtvQElWpwGB+EebQ4x1WXPYnttwBRMg5FVrG+ZDK4c223sQIKMk3Vd/NG9efgH6hDGN+Bt7uhJfY77MJsaqNaW2xXwur0pqpC9lvQvttlo5ISLFIE6mjnIku/2+WD016e4jQXOrXrLq5WX3+/+1dXch011V+1j4z31fUXrSmKSEGTUoE442mpRRacqNok5uqIORGUshNoQULLTTamyIEqqCXCpUWRMRSqWIu1aCI4F8raX6axiZtaGNCouTC4kW+mbNXL/Zaa6+9zz5nzrwzfjPzcRa83zfvvOdn/5z9nLWe9bPFXChIauUVXKYtC1izpNpbcB0RUkOrhRcIkCQsv9M5FQS1EJtAdqkyi8t+BVrFHM4OoI6g9NcZfLZnSccfGUDqc/T3mmuqh1HBeSopboZcGSiIaAXg1wC8V79j5rcAvCWfv05ELwH4aaQdz3fLGEioKq2DRqEECd0B69ra7OYCRXtlv9mlXm8NgLTArC8NVy72PvMJErZsYJH7Pp3S7kHC2dbp8n1prgBwrKJ5PMwe38qrYtuD37pRgERRP9S7DwfazVClbwrlwCoESF0ETWH22kQFFB3nubB6mh2GRRfa9yxAovZgGVi4NjOnwrnBbYzEOTaDtxvghltwq1UGmtorQBM7nWsdDQELnTetLoWtqy5WhOE7LbbSbJsAQGI2kEBSI3vWjh3jINTE0xfRAXKIRvGLAL7FzK9Yu4jeBeBNZu6J6B6kvUe/M+tqUyChvIQuuLXjJN52HaxJXGtRiQmwWoUJ9rMmse0zSGw26Y2s2opf6F4oZBeeahcN+7J9bhskCm7COJDMvch4ZiKwCF/nrAlpQRoXXOV5Cd8+84K0wGEsrFtU5aTOw/BLS7KxvJw5pP+1Cp2CRiI0NeiN8r3GnlvlJXzQ3PpaIlSvrbM7u+FCp21vb1CS0HAAFp3KmwxSFBm8WgEhFvyNtY+pLOHna1rUYAGAOYJ6yRnqe6eFkd0vE6wVSIzGzYQSNBpzY23S9TNFoB8gV9p7lJm/COBhlGYHADwA4HeIaIv0KHyMmd+c0xCz7cdAAhiCxLWc5VmUTxcQpgBgG/POWMJFzAIJr/LVErlauGgDhLZZ1WZ5g/nir/YANMAph02X/EQiXXtgo5yEizVoeDnGtIkmOFh5ulg+zL1kgrKLNyXYYuCQwQICHnGVtA7jWWpNqhbHS/hUewUJ87p0kq0KwGqExGiRoGoWWCUoc2/34E3uJ3HiFFjBaZBNG9P8eK1DPzuwSOZHVdxXXwh1WH6akAwQMjd+jmyu+740T+t50mvp+S2QOBKhedW9R8HMH21891UAX927FQxY/sYUSKgqel22yru+zlvlrSuTQweu51TjUSszaf6DB4kZnoDBRj8F89zYV0TbrG93n5zmM2DTuJVxDWo2uLdruqfa3DG7i31afaOdg7JvdWAbMAQSf46llMuDHSn9cEiKmg63mPfQP4XMYbCL0PTeidZ4Z5exhn87s0X301DC0M91r29+2evCjSkxg29oHxxYQMCk64DQXkxFXox3Uep8yP01i9iK+ypgjL1sVAMABgucvanBUqioYToMzvcvuSrm51A5k8hMzklePn/Dp0AHMm1CC8hOgYSWIrP04I2Cg5o22auh17eKVEB+O3tp2fFj/n8XJ6ELpF2Vy5kd/vSW2eFBYuuAzoOAb1oNEv7aI6IRjCz2NHdSg2Hbp6ChEKA7cmtgUXE+iUJngIGsTegPZLzrZ78uKaB9F68Lr7u0S5nLcUk3RWpjn79TsEiaRMwktXIXaoYod9DKBxKtojBBKk+V1f/U+dG6EJ3uBkbt5wYo58WDEeA8FqEEDpXqRWXXGNNwD5TzAAqGS/LSOIlsr5smIenVqobyurO3VWmvwrIoLZPS1TQYhl4ToBWRajJRZcI7MFXFajDBVgvC2an+fCXyKn9/fntk86mOKSmk9cD4hdCqI6kmlTc/1N7unModJWCoz1WciEmiEfVeAhadVor2wDdlfjhPjQNZ3efTyuWv8jVSkZ4UaRv0eE8w6jjrXibMADJYmKnSIbmkFaD0xVGDRSW+QLJpIKoJ+MjbqcWer1Yey2J+jHE6Y1rF6PWvJucBFOAclu2R36veonrSapVrRqjLLTNKokVIurXfMt7VjvBp2wSI/1tVzG7c1lPZhdo1qNS5JoUf3YGiP79mzfV+Mj6FJ8iGUUrAtfbokD4U8Qhqc3u1Vr06or1QJ2/JPoIoaxW8jaB1kIQ0pM25xLGU+5ErVhfRmXOINX984ZZNG+/ELo8P9QwKDNpyLuGgi1WrP3m3t0ZpQsYTVwCLOv5FpwD5vlauQEF/ruOhiIEIqWVjZHP9LPq2eLLzQDkPoBhoFI7YUZNDXWRCaFlKcyBYjULO9RhKkOiTySHXtXLudhN54NQO75xduqvpU6AxFaMggGX8hD+n9VB4DsdrXVNtsks6E8jHb3ibuw4C0m35lFBznhgOIW2euw2WuUo9gyKlWBVddUAmPDsxW3RM1HU3JdpeILtlAyFK8pnFKshGvoFYuc28i5aaUH1v2yNqfEUy5fqcuGXai2uDgsUcscUaMmDIGLPeY0qsVoYroIw+A8eU6dvKP6riPw6R8wAKYPiGVG3CQofF9FA12FRn2a5NVeAGSNjC8m8Hn6LTIKxoSsWc3SnvT/TfQYhWl2jmJC9ocY9xn48vOJxGe1rek/LibbPDn0uy34ZqFRTBlDQKNQMSSAgo9wFhC4Qt207cJIsUgKtY7UyQhgyCzLRZmgujpowkobEOLQNhC0Q5JwBg7sQbspLYGcnbWYtW6cYvVRlPXp3U72rexDQttIqaXLTkNs4anQOMNPTumWuNgV6zDpIy4OhHzZem9nAgOHg5C6BgbtQbAKxkPfkU54qPIEbyamhF4k2fPRwbF5INuEHuyrcpycMRXZCNVGwG0FbdfFBWvThrIlRt/q4rJ7au/aj8hL+OZlVqen1D65ol3j2rvMdYII7e12sVgcBbpFyJbZ+0im2HsI0IG04/K1nMgRB6zhWpnTZA7vNOITfXQfqgJKm6X5HuwYEQNiyolAC261k8ECtQvJZdokpE9s4EYUrahWkWnE0f5gwWfmPjavx8jY0kEqgn2kBBSjbGfQAkzj3KTW1jRIM4osmhchZAkWba2elKYJov3fngXf5G8g9zu/CtL52mElxpNBVVs6VEkW6u29p4yI5vgUTNWvvrdxUhCVjbBjyDi8AsyuH5Kk9TjHZxLQeotefFqdMpSCk/5BRITCOSHbdlbAHwNl2DQgA2far2dCOgW5HV9bANxftMbtquVi3PTy31MbpYKO+YFbtkgtjuWX1iGohJHBsBMTJCDKBVB+470HqdxnCtrlIBCe03VXPh72+8lgRiNcYuz0F0H9OYmJYxZkIA+TqBjBhNv4+AhsoEQBxaAk/lPICC7Z8kVm8gZwoWdrK6oQBosA2kEpXVtOxTVWdbEBbDUD2kxNCS8cnAbdjOmsfRAolWDIJ4DIodpxR49Bhj4yuvhL9nCAkMPdE7AJaJRccMrUbtk7OKYKLG+QYWCqJArgJlb9Igmpukkrvt9bhLJ4RNoyy9hqGT7uilb+hQ8SiOzPT9JJh5oxGhej1AigEzQD3lgjrrzpkg62L8By5ameuB+QGUJKPXJvz3Le0SHjBoGugB0yJNA5kAjXyfEiRKTeNwsDgPoPBCw01wigUuk0xac0Bsd+ozWJgm4UJyLa5fryVkFgRwUuqz0PbVWwNACRK70qO9u1IfuCoM2HiGluix+hZvRY76B6WOKqwlOG3CkYlE5bvQk6BWoo4pV4kChK9IGwBR2KYCLgJAHSW7gOVNTz1AW7f9nW+vfRbwr/mJ4LSfAbhDzJekWfjhIE73jSsCxVRWj/toKfJ0bZ3mfO3MNzev1l81PzRtXJ+ZOSQsMNQwxRyddKXb8HgXp4CGTpu+bFocRw0St5571EmtTTRAIk+wq0ilWkQVl2DXUK3EVFpK5JYuyCr/wJNDGl4+WmNiTAwkZP8LBQa7TmVuFSaRm/Q62cxWRSgDxRph5HbNZs7IWJuRF7O6S3Xx2FgFYBsSaUy5gAv1UhU8veARNhLLUtvLEt7MqoqrNmE7pHWlJtnSfAQslLdIpoeaPX4z4QDqYkrW6mNysfe9FSvivlqcam7J51GtbY6btyDJq7EdO6WYmxGew2sZxe0a83qE4KvzAgp9sFvhy0AGBvls9Syj29avDlW2h68kQbOJ0EBfeYDyfp9XAYl0fWPBo+/HRPUroAS7WouRcfEl9Hyx2pb4AKaCHOQJlRko2m6xFRKrkrwhW/i9MQIgrtJ8H5Js3fR924Xc1CasMlblLmaUqh5lU0ffI9QBQQO9xIXOIdcbpVUH6lfpBbBeg7AZmiA2Jt38OZ+SKcBp3hc2v2WldAy0jKYcOULzzICislNVvCbhzQZdyLUXQB9Ss3ert5LTSAbuJA3I4h27oe8hhctMfm8GWaU/2uJtJql5kDDtwG801FaLB1GC7jh7CxUhxZXnBsjeG45A7MGRqozNNJ5h22WV3ZHMCr6+z5qR29QmtNq3EthmXijvQUWOSeyQzI0OYNmcRwlQCsiEqoaUr1YgjuDYldWszPQQ88NrYa2w632lNbZeVDus+KsxwJgT8TnH5JmS8wGKMR+z2o1K4Kjv21fEluNs8DUQy7PsLW+EBd2kB98qbQ/IwxGQGCOUit8r27TFbyibLSRikyz1IFFVrdJdyYByjIrxrMdg7O3WevOpu7RzC6ZDju0gsmCmpNpHY+Zz0ZYUx1BsdVe1UQsIl9pEsO30cntgxWuJJeZAVBrLYg0EDix8hjvXnol0H46dq1LVaJd3UY9JvQj3BBA1DweBdzVwy/ejgFHLEfM9zgcogKyCqrC4LfteyEbkxaML2A+uqxcx0CbsmhkgfFUqiDq9EyRqEtG3tbVI9W/F7xNmxwRZOigLp9/HxmZDrTYCzo/fuK9v3y41WTwiBuIyJyRmoBHQMtbkNQo17bR9hYlYaxNUVCtL2gSaIeM+Ec24i7F+1C8Sv5eJRoMyl4F5Rf8b3oVd0hhXz0sNAKPFO1WAMfCiVJnJx5LzAYqWKhYZ6S2pHgnkRV2bDSHkSa2vxdUiUEBwbzqrkaluszGQMF6gnHDbX7SWXQ9S6+1dn+f4GlPRpc+zgmpafI9vn2kwezz0Th1nxKw5xBSxWWgwvnCxmnRF8zQNP2TvlMuTqME+EZYSMh7T72W6O4GJC3Ol6HKQUHM1QZR3qU1eOyHmRRuH7b+qDEoXVN+PktSuTUWxIv1+cOjh7T0foKhFC5FEQF9/hQ2pD2adbLWeUKnV3GiARE5IG9kRHSj5gUoGG9DMkT1Aoqi1KWq6T26bVYeyuBeX2tjUG0jf+jW56MYhR786baLoWsyJearJAe66yrsEDHYiI3Jl8AUkespaRcQQEJRfaHWnXjiqVXADeKNwFK3x8YFN+yzGmeULJj1acoye7+8/BkCHyHkChcvOK4bfaRFWGdtn8tW7Yek5VkWKs4dkF0ioNAjEQQjtVRCbqsU2EZvhvQLZ9HCl2sRlSVzVePAL3MtUGPnOZqt2Iu3XwCy9xnYkAUoBWiJmB2NBQ24iaxXuUDNxZBvFHuCeEIWLoMRxDiXCkaBlf1jGyupk+r7OWXQ2dzvGsWUGeZDwpKnzaKVLT5gUI16SY8r5AYX321fSWpR2rNmYznaugpy8J6MFEkOmuUEg+uxL1nRioBnROSY1SOQODo/zJgdR+fAGwmiZs/qh3JVMNHUukLUEc1+Hqh2AbZZUi95Lx9rNY1GgR8hM5SY0Y7Roj9MiQs+IfTIjSBSC2n3KvisFP1URhyHZLxSFm5jryjyW+PmpAWOOdgEMAOOozdt1ABHdRUR/T0TPE9FzRPSb8v07iehviejb8v873Dm/RUQvEtELRPTLe7dK7ED7qdKri+I2U9cwzSGZGnkToTZIDMvhufJ7nmirbejhoI0N5jRfMDi2LMprtTd3na/EX6X9mOjYqXZWh1lL30282UEBCKlNJLEJhQgpbGX6NrkaVwESbu4yCDk+QutPeM0FMD4h9AySjNWwTZpF0LoYri9qfZT7umTTtQDKWmP0JtbcxTfT7KQGMJC8CMi/ENx876yYbm1ovPQOlDm93wL4FDP/DIAPAPg4Ed0H4DEATzLzvQCelN8hf3sYwM8C+DCAPyQabH87Ltq5yNn9JuDAMWeZ+g1k7VT1YFSejLyhT292chMk/MOrmau6SKVsPIn7znseyC3M5iQ6gKCuSwvY/Ujjh6d5k8Pb8cWl3RvXP9Du87AatzfFxh/swUbA3n3pQNOHhmeNrS8zXjfbAuwHItcybsIAQ8YvireDIYWJOIOE/+lZuAvOmoX+H3XrSAcawLA93vtWhJvXQBIax0y8JFrX82Nr98/zVmiQ+4CFypFAY05x3deQNvYBM/+AiJ4HcCeAjyBV5waAPwHwDwA+I99/Wfb4+C4RvQjg/QD+eee9Gu6eqVTuQdIVUHpKfMFUjeV3Ho4mSHhTQ3NN6nqKkTGVVJWuU05osejUrtd+NUKvmyaHij/fTklZjElFLaMdm6IuyhkPUS5HJ/uLmFbl2uHJ1CqJTmVQUKeDXafI6xCAYA+CgCQDIoVM9BHYCEejWNsl0lNJTnWlpj9y+ePbWDRKVHs/V8W4hTT3HCyM2o97cyynzJhCswjNz2beyv0tEncQa1GR4EeUvTgKSnuQ/jyAfwXwbgERMPNrRHS7HHYngH9xp70i39XXylsK4kdscgbq/xjJZzxBhCVdEWVPSYglZWDaCOeSaHNAYrUqXX1AestZ6HPLxneTPEaCBux+m7e0Cc8LKC8SglsAUvsTFa/gx3JfF18I2X2phKPl4ch9DUSdRtjQ1AZj5MnR+qfGuMh5HW+BgCgcRMi3JkA9I0pgFiRmzVP4uJy6bfb2ruY+CvGpWiVz+aztIC2bmtrUuSHvVqZ8WBMsiMp+jl3vCjIbKIjox5BK8X+Smf93guxp/WHwVPotBYnov/8u/sX/Afifue2xK/YANlMH7inHv95t2KdfQOpbPHI79r0/kB/+G82j9u9XLf3otU8ph/drl8zkvI8ieeXdBuAnr3qZWUBBRGskkPgzZv5L+fp1IrpDtIk7ALwh378C4C53+k8AeHXq+sz8LiL6GjO/b7/mn78s/bosucX79VNXPX+O14MAfBHA88z8B+5PTwB4RD4/AuCv3fcPE9F1IrobaVvBf7tqAxdZZJHTyxyN4oMAfgPAM0T0lHz32wA+D+ArRPQogO8B+HUAYObniOgrAL6J5DH5OPNISuMiiyxyETLH6/FPaPMOAPALI+c8DuDxPdvyhT2PvxRZ+nVZsvSrIXSsBJdFFlnk1pXjx3oussgit5ycHCiI6MMS6v0iET126vYcIkT0MhE9Q0RPEdHX5LvRUPdzFiL6EhG9QUTPuu/+/8L2b5KM9OtzRPRfMm9PEdFD7m9n36+bkmZR5FTc5B+kOkkvAbgHwDUA3wBw3ynbdGB/XgZwW/Xd7wF4TD4/BuB3T93OmX15AMD9AJ7d1RcA98ncXQdwt8xpd+o+7NGvzwH4dOPYi+gXgDsA3C+f3w7gP6XtR5uvU2sU7wfwIjN/h5lvAPgyUgj4rSQfQQpxh/z/K6drynxh5n8E8Gb19VhfLGyfmb8LQMP2z05G+jUmF9EvZn6Nmf9DPv8AgE+zOMp8nRoo7gTwffd7M9z7goQB/A0RfV1C1IEq1B3A7aNnn7+M9eVWmMdPENHTYpqoin5x/ZpKs8AB83VqoJgV7n1B8kFmvh/Ag0hZtg+cukE3SS59Hv8IwHsA/BxSAuTvy/cX1a86zWLq0MZ3k/06NVDsHe59zsLMr8r/bwD4KyR17nUJcUcV6n6JMtaXi55HZn6dmXtmjgD+GFkNv5h+TaVZyN8Pmq9TA8W/A7iXiO4momtIdSyeOHGbriRE9KNE9Hb9DOCXADyL8VD3S5RbMmxfF5PIryLNG3Ah/bopaRZnwNg+hMTSvgTgs6duzwH9uAeJSf4GgOe0LwB+HKmwz7fl/3eeuq0z+/PnSGr4BukN9OhUXwB8VubwBQAPnrr9e/brTwE8A+BpWUR3XFK/AHwIyXR4GsBT8vPQMedricxcZJFFdsqpTY9FFlnkAmQBikUWWWSnLECxyCKL7JQFKBZZZJGdsgDFIossslMWoFhkkUV2ygIUiyyyyE5ZgGKRRRbZKT8EzeL24vtYfkwAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"square = shotset.intensities[0].reshape(201, 201)\n",
"\n",
"plt.figure()\n",
"plt.imshow(square, vmax=10000 * num_molecules)\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": 39,
"id": "02958fde",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"12:06:55 - Interpolating shots...\n",
"12:06:55 - \u001b[80D\u001b[1A\u001b[Kinterpolating shot 1/1\n",
"12:06:55 - ... complete\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"/Users/TJ/opt/miniconda3/lib/python3.8/site-packages/thor/xray.py:770: FutureWarning: `rcond` parameter will change to the default of machine precision times ``max(M, N)`` where M and N are the input matrix dimensions.\n",
"To use the future default and silence this warning we advise to pass `rcond=None`, to keep using the old, explicitly pass `rcond=-1`.\n",
" pix_n, resid, rank, sigma = np.linalg.lstsq( A, (i-p).T )\n"
]
}
],
"source": [
"q_range = np.linspace(0.1, 2.0, 100)\n",
"polar, mask = shotset.interpolate_to_polar(q_range, num_phi=360)"
]
},
{
"cell_type": "code",
"execution_count": 40,
"id": "ce76e857",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAERCAYAAAB2CKBkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAAgtElEQVR4nO3de3hc9X3n8fdHkmUZ27IlWTbGF8kYAzbUBks2l8SBhBJM0qyTDWmANJRsUuokdNvu003zbK/btM92m+3TWwAvSUge0m1oGkJLginNFRKMwbIxGGMw8l022PL9fpH03T9mnAyyJI9sHZ0Zzef1PHo05/x+Z/Sd8fF85pzfuSgiMDOz0lWWdgFmZpYuB4GZWYlzEJiZlTgHgZlZiXMQmJmVOAeBmVmJK8ogkPSQpF2SXsmj77skrZLUIem2nPlXSXpO0lpJL0v6aLJVm5kVpqIMAuDrwMI8+24F7gb+qdv8o8BdEXFF9rn+VtLYAarPzKxoVKRdwLmIiGckNebOkzQduA+oJ/Mh/xsR8VpEbM62d3V7jvU5j3dI2pVddn+ixZuZFZiiDIJePAgsjog3JF0D3A+8J58FJc0HKoENCdZnZlaQhkQQSBoFXA/8i6TTs4fnuexE4BvAr0dE19n6m5kNNUMiCMiMdeyPiKv6s5CkauAJ4A8jYnkShZmZFbpiHSx+m4g4CGyS9BEAZczpaxlJlcBjwMMR8S+DUKaZWUFSMV59VNI3gRuBccBO4E+AHwEPABOBYcAjEfFnkuaR+cCvAY4Db0XEFZJ+DfgasDbnqe+OiNWD9TrMzApBUQaBmZkNnCGxa8jMzM5d0Q0Wjxs3LhobG9Muw8ysqKxcuXJ3RNT31FZ0QdDY2EhLS0vaZZiZFRVJW3pr864hM7MS5yAwMytxDgIzsxLnIDAzK3EOAjOzEucgMDMrcQ4CM7MSVzJBsH7nIb7wvVc5fqoz7VLMzApKyQRB276jfPVnm1ixeW/apZiZFZSSCYJrL66jsryMp19vT7sUM7OCUjJBcEFlBddcXMvT6x0EZma5SiYIAG64tJ43dh1m+/5jaZdiZlYwSi4IAO8eMjPLUVJBcMn4UVw0poqn1+9KuxQzs4JRUkEgiRsuq+fZ1j2c6uxKuxwzs4JQUkEAcMOl4zl8ooNVW/alXYqZWUEouSC4/pI6Ksrko4fMzLJKLgiqq4Yxt6HGQWBmllVyQQCZo4fW7jjIrkPH0y7FzCx1JRsEAMta96RciZlZ+koyCGZOrKa6qoJlG3anXYqZWepKMgjKy8S1F9fx3EZvEZiZlWQQAFw3vY5te4/Rtu9o2qWYmaWqpIMA4LkN3iows9JWskFw6fjR1I6sdBCYWckr2SAoKxPXZccJIiLtcszMUlOyQQBw7fQ63jxwnC17PE5gZqUrsSCQ9JCkXZJe6aVdkv5eUquklyXNTaqW3lx3cXacwEcPmVkJS3KL4OvAwj7abwVmZH/uAR5IsJYeTa8fSf3o4R4nMLOSllgQRMQzQF93il8EPBwZy4GxkiYmVU9PJHH99DqWbfA4gZmVrjTHCCYB23Km27LzziDpHkktklra2wf2YnHXXVzH7sMn2NB+eECf18ysWKQZBOphXo9fyyPiwYhojojm+vr6AS3i9PkEyzf2tfFiZjZ0pRkEbcCUnOnJwI7BLmJq7QWMHz2cFZsdBGZWmtIMgseBu7JHD10LHIiINwe7CEnMm1bLik0OAjMrTUkePvpN4DngMkltkj4pabGkxdkuS4GNQCvwZeAzSdVyNvMba9lx4LivO2RmJakiqSeOiDvO0h7AZ5P6+/3R3FgDwIrNe5lcc0HK1ZiZDa6SPrP4tMsvrGb08Ape2OQb2ptZ6XEQkLk/QVNjjQeMzawkOQiy5jXW0rrrMHuPnEy7FDOzQeUgyJo/rRbAWwVmVnIcBFmzJ4+hsqLMh5GaWclxEGQNryjnqsljvUVgZiXHQZBj3rQaXtlxkCMnOtIuxcxs0DgIcsxrrKWzK3hx6/60SzEzGzQOghxNDTVI0LLFu4fMrHQ4CHKMrhrGZRNGs3KLTywzs9LhIOimubGGF7fup7PLN6oxs9LgIOimqaGGwyc6eP2tQ2mXYmY2KBwE3TQ3ZE4sW+lxAjMrEQ6CbibXjGD86OEeJzCzkuEg6EYSTQ01tDgIzKxEOAh60NRQQ9u+Y+w8eDztUszMEucg6EFz4+lxAm8VmNnQ5yDowayJ1QyvKHMQmFlJcBD0oLKijDlTxnqcwMxKgoOgF00NNazdfoBjJzvTLsXMLFEOgl40N9TQ0RW83LY/7VLMzBLlIOjF3Kk1AN49ZGZDnoOgFzUjK5leP5IXtzoIzGxocxD0Ye7UGlZu2UeEL0BnZkOXg6APTQ017Dt6ik27j6RdiplZYhwEfWhqyIwTrPIdy8xsCEs0CCQtlPS6pFZJn++hfYyk70p6SdJaSZ9Isp7+ml4/iuqqCp9YZmZDWmJBIKkcuA+4FZgF3CFpVrdunwVejYg5wI3AX0uqTKqm/iorE1dPrWGVg8DMhrAktwjmA60RsTEiTgKPAIu69QlgtCQBo4C9QEeCNfVbU0MN63cd4sCxU2mXYmaWiCSDYBKwLWe6LTsv15eAmcAOYA3w2xHR1f2JJN0jqUVSS3t7e1L19qipoYYIWL1t/6D+XTOzwZJkEKiHed2Pw7wFWA1cBFwFfElS9RkLRTwYEc0R0VxfXz/QdfZpzpSxlAnvHjKzISvJIGgDpuRMTybzzT/XJ4DvREYrsAm4PMGa+m3U8Aouu7CaVT6xzMyGqCSDYAUwQ9K07ADw7cDj3fpsBW4CkDQBuAzYmGBN56SpYSwvbt1PZ5dPLDOzoSexIIiIDuBe4ClgHfCtiFgrabGkxdluXwCul7QG+CHw+xGxO6maztXcqTUcPtHB+p2H0i7FzGzAVST55BGxFFjabd6SnMc7gPcmWcNA+MWJZfuYOfGMIQwzs6LmM4vzMLX2AsaNqmTlZo8TmNnQ4yDIgyTmTq3xgLGZDUkOgjw1NdSwec9Rdh8+kXYpZmYDykGQp5+PE/h8AjMbYhwEebpy0hiGlYuV3j1kZkOMgyBPVcPKuXLSGG8RmNmQ4yDoh6apNbzUdoCTHWdcDsnMrGg5CPqhqaGGkx1drN1xIO1SzMwGjIOgH+ZmB4x9oxozG0ocBP0wobqKyTUjfD6BmQ0pDoJ+amqoYeWWfUT4AnRmNjQ4CPqpuaGGnQdPsH3/sbRLMTMbEA6CfvI4gZkNNQ6CfrpswmhGVpY7CMxsyHAQ9FNFeRlXT62hxVciNbMhwkFwDpoaanjtrYMcOn4q7VLMzM6bg+AczGuspSvgxa370y7FzOy8OQjOwVVTx1ImaNm8N+1SzMzOm4PgHIwaXsGsi6pp8YCxmQ0BDoJz1NxQy4tb93Oq0xegM7Pi5iA4R82NNRw71cmrOw6mXYqZ2XlxEJyj5oZaAO8eMrOi5yA4RxeOqWJK7QgPGJtZ0XMQnIfmhlpWbPYF6MysuDkIzkNzYw27D59g696jaZdiZnbOHATnYV5jZpxghS83YWZFLNEgkLRQ0uuSWiV9vpc+N0paLWmtpKeTrGegXVI/iuqqCo8TmFlRq0jqiSWVA/cBNwNtwApJj0fEqzl9xgL3AwsjYquk8UnVk4SyMjGvsZYXHARmVsSS3CKYD7RGxMaIOAk8Aizq1udO4DsRsRUgInYlWE8i5k+rZWP7EdoPnUi7FDOzc5JkEEwCtuVMt2Xn5boUqJH0E0krJd3V0xNJukdSi6SW9vb2hMo9N/OnnR4n8FaBmRWnvINAUo2kKyRdLCmf5dTDvO7HWVYATcD7gVuAP5J06RkLRTwYEc0R0VxfX59vyYPiykljGDGsnBc2OQjMrDj1OUYgaQzwWeAOoBJoB6qACZKWA/dHxI97WbwNmJIzPRnY0UOf3RFxBDgi6RlgDrC+vy8kLcPKy2hqqOF5B4GZFamzfbP/NpndOwsi4rKIeGf2m/kU4C+BRZI+2cuyK4AZkqZJqgRuBx7v1uffgAWSKiRdAFwDrDvnV5OS+dNqee2tgxw46hvVmFnx6XOLICJu7qNtJbCyj/YOSfcCTwHlwEMRsVbS4mz7kohYJ+nfgZeBLuArEfHKObyOVM2fVktEZpzgl2dNSLscM7N+yevwUUlz+2qPiFW9zF8KLO02b0m36S8CX8ynjkJ11ZSxVJaX8YKDwMyKUL7nEdwPzCXzzV3ALwEvAKfIDAC/J5HqikTVsHKumjLW4wRmVpTyPWpoM9CUHR9oIhMKrRHx7ogo6RA4bf60Wl7ZfoAjJzrSLsXMrF/yDYLLI2LN6YnsfvyrEqmoSM2fVktnV7Bqq687ZGbFJd8gWCfpK9nrAt0g6csU4dE9SZrbUEN5mXw+gZkVnXzHCD4BfBr47ez0M8ADiVRUpEYNr+DKi6p5fqODwMyKS15BEBHHgb/J/lgvrr24joee3cSxk52MqCxPuxwzs7z0uWtI0nclfUDSsB7aLpb0Z5L+S3LlFZfrptdxqjNo2eKtAjMrHmcbI/gNYAHwmqQVkpZK+pGkTcD/BVZGxEOJV1kk5jXWUlEmlm3Yk3YpZmZ5O9uZxW8BnwM+J6kRmAgcA9ZHhO/P2M3I4RXMmTKW5xwEZlZE8r76aERsjojnImK1Q6B310+v4+W2/Rw87usOmVlxONsYwSFJB3v4OSTp4GAVWUyum15HV8AKH0ZqZkWizyCIiNERUd3Dz+iIqB6sIovJ3Kk1VFaUefeQmRWNRG9eX4qqhpXTNLXGA8ZmVjQcBAm4fnod6946yL4jJ9MuxczsrBwECbhueh0R8PwmbxWYWeFzECRg9uSxXFBZ7t1DZlYUHAQJqKwoY15jLc+27k67FDOzs3IQJGTBjHFsaD/Cjv3H0i7FzKxPDoKELJhRD8BP32hPuRIzs745CBJy6YRRTKgezjNvePeQmRU2B0FCJLFgRj3Ptu6msyvSLsfMrFcOggQtmDGO/UdPsWb7gbRLMTPrlYMgQQtm1CPBT9d7nMDMCpeDIEG1Iyu58qIxPOMBYzMrYA6ChC2YMY5VW/dzyJelNrMC5SBI2LsuraezK3yWsZkVrESDQNJCSa9LapX0+T76zZPUKem2JOtJw9ypNYysLPf5BGZWsBILAknlwH3ArcAs4A5Js3rp97+Bp5KqJU2VFWVcN72Op9e3E+HDSM2s8CS5RTAfaI2IjRFxEngEWNRDv98CHgV2JVhLqm64bDzb9h5jQ/vhtEsxMztDkkEwCdiWM92WnfdzkiYBHwKW9PVEku6R1CKppb29+Hax3HT5eAB+sG7IZp2ZFbEkg0A9zOu+b+Rvgd+PiM6+nigiHoyI5ohorq+vH6j6Bs1FY0cwa2I1P3IQmFkBSjII2oApOdOTgR3d+jQDj0jaDNwG3C/pgwnWlJqbZo6nZcte37XMzApOkkGwApghaZqkSuB24PHcDhExLSIaI6IR+DbwmYj41wRrSs1NMyfQFfCT9d4qMLPCklgQREQHcC+Zo4HWAd+KiLWSFktanNTfLVSzJ41h3KjhHicws4JTkeSTR8RSYGm3eT0ODEfE3UnWkrayMvGey+t5cs1bnOzoorLC5/KZWWHwp9EgumnmBA6d6KBl8960SzEz+zkHwSBaMGMclRVl3j1kZgXFQTCILqis4PrpdfzwtZ0+y9jMCoaDYJDdPGsCW/Yc5bW3DqVdipkZ4CAYdLdccSFlgqVr3ky7FDMzwEEw6MaNGs410+p4Ys2b3j1kZgXBQZCC982eyMb2I7y+07uHzCx9DoIULDy9e+hl7x4ys/Q5CFJQP3o486fVeveQmRUEB0FK3v9LE9nQfoT1O32PAjNLl4MgJbdceSESPOGjh8wsZQ6ClIwfXcX8xlofRmpmqXMQpOhXZk+kdddh1u44kHYpZlbCHAQp+pXZFzGsXDy6cnvapZhZCXMQpKhmZCU3XT6Bf1u9nVOdXWmXY2YlykGQstuaJrPnyEmefr097VLMrEQ5CFJ2w2X11I2s5Nsr29IuxcxKlIMgZcPKy1h01SR++NpO39jezFLhICgAtzVN5lRn8N2Xd6RdipmVIAdBAZh1UTUzJ1bzqHcPmVkKHAQF4sNzJ/FS2wHWvXkw7VLMrMQ4CArEh+dOZnhFGd9YviXtUsysxDgICkTNyEoWXXURj63azoFjp9Iux8xKiIOggNx1XSPHTnX6UFIzG1QOggJy5aQxNDXU8I3nNtPV5fsUmNngcBAUmLuua2DznqM8/YbPNDazweEgKDC3XjmR+tHDeXjZ5rRLMbMSkWgQSFoo6XVJrZI+30P7xyS9nP1ZJmlOkvUUg8qKMu6cP5WfrG9n0+4jaZdjZiUgsSCQVA7cB9wKzALukDSrW7dNwA0RMRv4AvBgUvUUk49dO5Vh5WUs+cmGtEsxsxKQ5BbBfKA1IjZGxEngEWBRboeIWBYR+7KTy4HJCdZTNMaPruKOeVN4dFUb2/cfS7scMxvikgyCScC2nOm27LzefBJ4sqcGSfdIapHU0t5eGoOo99wwHYAHn/ZWgZklK8kgUA/zejwmUtK7yQTB7/fUHhEPRkRzRDTX19cPYImFa9LYEXx47mS+uWIbuw4dT7scMxvCkgyCNmBKzvRk4IzLa0qaDXwFWBQRexKsp+h8+sbpdHR28ZWfbkq7FDMbwpIMghXADEnTJFUCtwOP53aQNBX4DvDxiFifYC1FqXHcSD4w5yL+cfkW9vpeBWaWkMSCICI6gHuBp4B1wLciYq2kxZIWZ7v9MVAH3C9ptaSWpOopVve++xKOn+rkSz9qTbsUMxuiKpJ88ohYCiztNm9JzuNPAZ9KsoZiN2PCaD46bwoPP7eZj1/XwLRxI9MuycyGGJ9ZXAR+9+ZLGV5Rxl8+uS7tUsxsCHIQFIHxo6v49I3TeWrtTpZv9Hi6mQ0sB0GR+OQ7L2bimCr+4ol1vjKpmQ0oB0GRGFFZzucWXsaa7Qf4l5Xbzr6AmVmeHARFZNGcScyfVsufP7GOnQd9kpmZDQwHQREpKxN/9eHZnOrs4g8eW0OEdxGZ2flzEBSZxnEj+b33XsYP1u3i8ZfOOFHbzKzfHARF6BPvmMbVU8fyJ4+vpf3QibTLMbMi5yAoQuVl4ou3zeboyU5+959X0+mjiMzsPDgIitQl40fzhUVX8LPW3fzN932ZJjM7dw6CIvbReVP5aPMUvvTjVr7/6s60yzGzIuUgKHL/c9EVXDmpmv/2rdW+x7GZnRMHQZGrGlbOAx9rorxM/PpDL7DL5xeYWT85CIaAKbUX8LW757H78Ak+/tUX2H/U9y4ws/w5CIaIq6fW8OW7mtm0+wh3f20FR050pF2SmRUJB8EQ8o5LxvEPd17Nmu0H+PhXn/ddzcwsLw6CIeaWKy7kvjuv5pUdB7ntgWVs3XM07ZLMrMA5CIaghVdO5J8+dQ17jpzkPz/wLKu37U+7JDMrYA6CIaq5sZZHP309VcPK+ciSZXz5mY2+j4GZ9chBMIRdMn4U3/utd/Key8fzF0vXcffXV/jaRGZ2BgfBEDf2gkqW/FoTf/7BK3l+4x5u+uuf8PVnN9HR2ZV2aWZWIBwEJUASv3ZtA0/81wXMnjyWP/3uq7z/73/GT99o9z0NzMxBUEouGT+Kb3xyPg98bC6HT3Tw8a++wIfuX8b3X93p8QOzEqZi+0bY3NwcLS0taZdR9I6f6uTbK9tY8vQG2vYdY3r9SH61eQofmjuJ8aOr0i7PzAaYpJUR0dxjm4OgtHV0dvG9l9/kH5dvoWXLPsrLxIIZ4/jlmRO4aeZ4Jo4ZkXaJZjYAHASWlw3th/n2yja+9/IOtu09BsDMidVcM62W5sYamhpquLC6CkkpV2pm/ZVaEEhaCPwdUA58JSL+slu7su3vA44Cd0fEqr6e00GQvIigdddhfrBuFz99o50Xt+7n2KlOAGpHVjJz4mguv7CaaeNG0lg3koa6C7hwTBXDyj3kZFao+gqCigT/aDlwH3Az0AaskPR4RLya0+1WYEb25xrggexvS5EkZkwYzYwJo/n0jdM51dnFqzsO8uLWfax78xDr3jrIPy7fwomOrpxlYNyo4VxYXUXdqErqRg6nblQlY0YMo7qqgtFVwxg5vIILKssZUVnOiGHlDK8oo2pYOZUVZQwrL2N4RRkVZaK8TN7qMBtEiQUBMB9ojYiNAJIeARYBuUGwCHg4MpslyyWNlTQxIt5MsC7rp2HlZcyZMpY5U8b+fF5XV7Dz0HE27T7C1j1HefPAcd46cJy3Dh5n9+ETrH/rEHuOnHxbWPRHeZkolygrI/s7ExBlEmUCEBIIsr9zp88MkdxZb3uMepmf27/nUOozqnpp7G+8nS0Q892i77VXLw299U9qD0JfrzOv92yQ3u+03T5vCp9acPGAP2+SQTAJ2JYz3caZ3/Z76jMJeFsQSLoHuAdg6tSpA16o9V9ZmZg4ZgQTx4zg+um99zvR0cmh4x0cPHaKoyc7OXKig6MnOzl+qpMTHV0cP9XJyc4uTnZ0cbKzi87OoKMr6OzK/O6KzOPOriAi6ArojCDzeZT5HQFx+nH27+Z+XgVvm+jp4ds+4N4+v+fX1dfHYW8flv3+CM13gTw/u3rr1t+gG+jPyr6yJZ+3YNDf7xSNGzU8kedNMgh6Wl26v9X59CEiHgQehMwYwfmXZoNleEU5w0eVJ7YCm9n5S3J0rw2YkjM9GdhxDn3MzCxBSQbBCmCGpGmSKoHbgce79XkcuEsZ1wIHPD5gZja4Ets1FBEdku4FniJz+OhDEbFW0uJs+xJgKZlDR1vJHD76iaTqMTOzniU5RkBELCXzYZ87b0nO4wA+m2QNZmbWN58BZGZW4hwEZmYlzkFgZlbiHARmZiWu6K4+Kqkd2JJ2Hb0YB+xOu4g+FHp9UPg1ur7z4/rOz/nU1xAR9T01FF0QFDJJLb1d3a8QFHp9UPg1ur7z4/rOT1L1edeQmVmJcxCYmZU4B8HAejDtAs6i0OuDwq/R9Z0f13d+EqnPYwRmZiXOWwRmZiXOQWBmVuIcBHmStFDS65JaJX2+h/aPSXo5+7NM0pycts2S1khaLaklpfpulHQgW8NqSX+c77KDVN9/z6ntFUmdkmqzbYPx/j0kaZekV3ppl6S/z9b/sqS5+b62Qaov7fXvbPWlvf6drb7U1j9JUyT9WNI6SWsl/XYPfZJd/yLCP2f5IXMZ7Q3AxUAl8BIwq1uf64Ga7ONbgedz2jYD41Ku70bge+ey7GDU163/B4AfDdb7l/0b7wLmAq/00v4+4Ekyd9W79vS/72C8f3nWl9r6l2d9qa1/+dSX5voHTATmZh+PBtb38P830fXPWwT5mQ+0RsTGiDgJPAIsyu0QEcsiYl92cjmZu60VTH0JLZtUfXcA3xzgGvoUEc8Ae/vosgh4ODKWA2MlTWRw3r+z1pfy+pfP+9ebgnj/uhnU9S8i3oyIVdnHh4B1ZO7dnivR9c9BkJ9JwLac6TbO/IfK9Uky6X1aAP8haaWke1Ks7zpJL0l6UtIV/Vx2MOpD0gXAQuDRnNlJv3/56O01DMb711+Dvf7lK631L29pr3+SGoGrgee7NSW6/iV6Y5ohRD3M6/G4W0nvJvMf8Z05s98RETskjQe+L+m17DeUwaxvFZlrjRyW9D7gX4EZeS57vvrzNz4APBsRud/ekn7/8tHbaxiM9y9vKa1/+Uhz/euP1NY/SaPIBNDvRMTB7s09LDJg65+3CPLTBkzJmZ4M7OjeSdJs4CvAoojYc3p+ROzI/t4FPEZmc25Q64uIgxFxOPt4KTBM0rh8lh2M+nLcTrfN8kF4//LR22sYjPcvLymuf2eV8vrXH6msf5KGkQmB/xcR3+mhS7LrX1IDIEPph8yW00ZgGr8YkLmiW5+pZO69fH23+SOB0TmPlwELU6jvQn5xAuF8YCuZbxNnXXYw6sv2G0NmP+7IwXz/cv5WI70Pdr6ftw/WvdCf1zYI9aW2/uVZX2rrXz71pbn+Zd+Hh4G/7aNPouufdw3lISI6JN0LPEVmlP6hiFgraXG2fQnwx0AdcL8kgI7IXCVwAvBYdl4F8E8R8e8p1Hcb8GlJHcAx4PbIrEk9LptCfQAfAv4jIo7kLJ74+wcg6ZtkjmwZJ6kN+BNgWE59S8kcudEKHAU+0ddrS6G+1Na/POtLbf3Lsz5Ib/17B/BxYI2k1dl5/4NMuA/K+udLTJiZlTiPEZiZlTgHgZlZiXMQmJmVOAeBmVmJcxCYmZU4B4GVLEl/2m16VfbEntx5FZJ2S/pfvTzHS9lDE3tqu1PSSUl/2Et7g6QXJT3VS3td9qqUhyV9qVvb72WvOHlHHy/RLC8OAis5ki6S9CTwmeylhX8327SMzDHdud4LvA78qrIHk+c8z0wy/4feJWlkt7b3AJ8DZgE3S7q7h1I+CHw/Im7ppdTjwB8Bv9e9ISL+D/DrwGd6e51m+XIQ2JAi6Q+y35R/IOmbks74EAV+B1gJ3A/MA06fIPQkmQuO5boD+DsyZ8Je263tTuAbwH8A/ymnhl8C/hy4JSJayZwIdKek7h/4Y4Fdvb2WiDgSET8jEwg9eSv7HGbnxWcW25AhqYnMtWKuJrNuryLzgd/dSTJn4e6NiFNkLvsL8GMyZ+iefr4RwE3Ab5L5wL0DeC7neT4K3AxcBtxL9ho1EbGGzP0ByE4fIbNl0V050NW/V/k2XdnnMDsv3iKwoWQB8FhEHI3M1Rsf76XfF8ms+78p6YeSbgSIiKPAfkkXZfv9CvDj7PxHgQ9JKgeQNA9oj4gtwA+BuZJq8i00u5tpDpmLhp2r3cD4/vxds544CGyoOes1UyLiQET8JvAl4A+Bb0mqyjY/xS92D90B/LKkzWS2LOqAd+e0XZ5t2wBUAx/Op8BsmGwmM37wRM78D+kXt0tszuN1HCWzFbKph91OZnlzENhQ8gyZb+0jJI0mc235M0iaKen0ur+GzC6W00cLPQkslFRN5pr+UyOiMSIagc8Cd2SX/QgwO6dtEZlwOKuI6IyIBmAFmd1Lp+c/FhFXZX/Oem/c7JbAR4HJEdHjkUdm+fAYgQ0ZEbFK0j8Dq4EtwE976foO4Gtk7hX7EeAvInOLQCJinaRLs/N/FBEncpb7N+CvgG8B2yNie07bM8AsSRMj4s08S14P1PbVIbvFUQ1USvog8N6IeDXbPAbYFdnr/JudK1991Ias7HkCh7OHWvbYHhF/2sP8JcA3IuLZhOv7HJmbon/uHJefD/xDRFwzsJVZqfGuIStlP+lpZkQsTjoEsr4DXN/bCWV9yR4W+yDwDwNelZUcbxGYmZU4bxGYmZU4B4GZWYlzEJiZlTgHgZlZiXMQmJmVuP8P13RT8vVMOYIAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# sum along the angular direction\n",
"# zero below is the first (and only) \"shot\"\n",
"Iq = polar[0].sum(-1)\n",
"\n",
"plt.figure()\n",
"plt.plot(q_range, Iq)\n",
"plt.xlabel('q $\\AA^{-1}')\n",
"plt.ylabel('I(q)')\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "2bb4384e",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.8.5"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment