Skip to content

Instantly share code, notes, and snippets.

@paulromano
Last active July 23, 2022 17:28
Show Gist options
  • Save paulromano/40a636d2a49f31a4382ff8a4524b53cf to your computer and use it in GitHub Desktop.
Save paulromano/40a636d2a49f31a4382ff8a4524b53cf to your computer and use it in GitHub Desktop.
Supporting plots for NJOY bug report
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "67115d7d",
"metadata": {},
"outputs": [],
"source": [
"import openmc\n",
"from openmc.data import IncidentNeutron\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "161b3d02",
"metadata": {},
"outputs": [],
"source": [
"evaluation = '/opt/data/endf/endf-b-vii.1/neutrons/n-047_Ag_109.endf'\n",
"temperatures = [293.6, 600.0]\n",
"\n",
"# Generate ACE file\n",
"openmc.data.njoy.make_ace(\n",
" evaluation, temperatures=temperatures,\n",
" broadr='broadr', acer='acer',\n",
" purr=False, heatr=False, gaspr=False\n",
")\n",
"\n",
"# Get cross sections from BROADR ENDF tape\n",
"broadr = openmc.data.endf.get_evaluations('broadr')\n",
"broadr_xs = [IncidentNeutron.from_endf(x)[2].xs['0K'] for x in broadr]\n",
"\n",
"# Get cross sections from ACE file\n",
"acer = openmc.data.ace.Library('acer')\n",
"acer_xs = [list(IncidentNeutron.from_ace(x)[2].xs.values())[0] for x in acer.tables]"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "9c7db633",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x14b1077f73a0>"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEaCAYAAAD+E0veAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAABOZklEQVR4nO2deXxU5fWHnzMz2QgQICCLgIAoiMgmoriida0i7vuOWq221ta2LrXa7afVtlYtKrgUV3Bp3a37iiIKyiKCCogQ9jUs2WY5vz/mzmSSTCY3yUxmJjnP5xMyc+9733vukLnfe8553/OKqmIYhmEYAJ50G2AYhmFkDiYKhmEYRhQTBcMwDCOKiYJhGIYRxUTBMAzDiGKiYBiGYUQxUTAMAwARWSgi49Jth5FeTBSMlCIieSLysIj8ICLbRWSuiBxXq82lIrJERHaIyOsi0itm37UiskxEtonIahG5S0R8Cc7XTkTuE5GNIlIqIh/G7JsqIn+Oeb+3iKwRkeuSfd2ZTu3PAkBV91bV99NkkpEhmCgYqcYHrAQOA4qA3wHPiEg/AOfJ9P+ACUAX4HtgWszxLwGjVLUjMBQYDvw8wfmmOP3s5fy+Nl4jERkJvAf8WVX/1qQrq+6rXpEyjGzDRMFIKaq6U1VvVdXlqhpS1VcI3/j3dZqcADyrqgtVtQr4E3CoiOzuHL9UVbc6bQUIAQPjnUtEBgMnAper6gZVDarqnDjtxgBvATeq6qT6bHc8ibdEZLOIrBORG53tt4rIcyLyhIhsAy4SkV4i8pLTdomIXBZ7PhGZ7Xg760TkH872fKePTSKyVUQ+F5Hu9djyWxFZ5Xhb34jIj5ztHhG5XkSWOv08IyJdYo47WEQ+cfpfKSIXicjlwLnAbxzv7GWn7XIROdJ5nSci/3S8s9XO6zxn3zgRKRGRX4nIesfburi+z9HILkwUjBbFuentCSyM3Rzn9dCYY85xbr4bCXsKk+vpfgzwA/AHJ3y0QEROjdPmdeBaVX0ogZ0dgLedtr0IC9E7MU0mAM8BnYAngelAidP2NOD/ROQIp+3dwN2Ot7M78Iyz/ULC3lMfoBi4AiiPY8sg4GpgP1XtABwDLHd2/ww4ibAn1gvYAkxyjtsN+B9wL9ANGAHMVdUpjs13qGp7VR0f5yO4CTjAOWa487n9LmZ/D8f2XYGJwCQR6RynHyPLMFEwWgwRySF8M3pUVRc7m18HzhCRYSJSAPweUKBd5DhVfcq5oe4JPACsq+cUvQmLSSnhG+TVwKMisldMmwOc/f9rwNwTgLWq+ndVrVDV7ao6K2b/TFV9QVVDQFfgIOC3Ttu5wEPABU5bPzBQRLqq6g5V/TRmezEwMOLVqOq2OLYEgTxgiIjkOF7XUmffFcBNqlqiqpXArcBpTkjrHOBtVZ2mqn5V3eTY5oZzgT+q6npV3QD8ATg/Zr/f2e9X1deAHcAgl30bGYyJgtEiiIgHeByoInyzBkBV3wZuAf5D+Ol3ObCd8FN3DVT1O8Iexn31nKac8M3qz6papaofEM4bHB3TZhIwG3irgSfbPsDSBPtXxrzuBWxW1e0x234g/BQN4SfpPYHFTojoBGf748AbwHQnRHOHI5w1UNUlwC8I3/DXi8j0mGT8bsDzTnhoK7CIsIh0d3ENiejlXEPs9fSKeb9JVQMx78uA9k08l5FBmCgYKUdEBHiY8I3qVFX1x+5X1UmquoeqdicsDj7gq3q68xEOwcRjfpxttcsABwk/Qa8A3hCRjvX0tRIYUM++2v2uBro4IacIfYFVEBYzVT0b2AX4K/CciBQ6T9l/UNUhwIGEvZMLiIPjLR1MWATU6Sdi53Gq2inmJ19VVzn76vusGiqPvNo5V+z1rG7gGKMVYKJgtAT3Ex4NNF5Va8TMnWTrUAnTl/DoobtVdYuz/1IR2cV5PQS4gZqx/Vg+JHyzv0FEfCJyEHA44afxKI4onU44R/GaiBTG6esVoKeI/MJJunYQkf3jnVRVVwKfALc51zOMsHfwhGP3eSLSzQk1bXUOC4nI4SKyj4h4gW2EvZxQ7f5FZJCIHOEkeisIe0SRdg8Af3HyB4hINxGZ4Ox7EjhSRM5wPo9iERnh7FtHYtGbBvzO6a8r4bDeEwnaG60EEwUjpTg3q58QTliudUa77BCRc50m+cBThGPSnwEzgZtjujgIWCAiO4HXnJ8bY/pfGOnLudlPAH5MOG/wIHBBTP4iijPS6RTCN9mXnXxG7P7twFHAeGAt8B1hgamPs4F+hJ+mnwducUJjAMcCC0VkB+Gk81mOOPYgnKzeRjjs8wHhkFJt8oDbCYvYWsIexw3OvrsJD9t9U0S2A58C+zvXsML5LH4FbAbmEk4aQ9hzG+KEnV6Ic84/Ew6zzQcWAF8424xWjtgiO4ZhGEYE8xQMwzCMKCYKhmEYRhQTBcMwDCOKiYJhGIYRxUTBMAzDiJLV1R27du2q/fr1S7cZhmEYWcWcOXM2qmq3ePuyWhT69evH7Nmz022GYRhGViEiP9S3z8JHhmEYRpSsFAURGS8iU0pLS9NtimEYRqsiK0VBVV9W1cuLiorSbYphGEarIqtzCoZhtA38fj8lJSVUVFSk25SsIj8/n969e5OTU6cie72YKBiGkfGUlJTQoUMH+vXrR7gSu9EQqsqmTZsoKSmhf//+ro/LyvCRYRhti4qKCoqLi00QGoGIUFxc3GjvykTBMFzy7brt+IN1ljswWggThMbTlM/MRMEwXLBxRyVH3/UhN79Q34JwRmvH6/UyYsQIhg8fzqhRo/jkk08AWL58OQUFBYwYMYIhQ4ZwwQUX4PdXLy44Y8YMxowZw+DBgxk8eDBTpkyp0W8gEKBbt25cf/31NbaPGzeOQYMGMWzYMAYPHszVV1/N1q1b69gzdOhQxo8fX2NfczBRMAwXbK8IL0c8c9mmNFtipIuCggLmzp3LvHnzuO2227jhhhui+3bffXfmzp3LggULKCkp4ZlnngFg7dq1nHPOOTzwwAMsXryYGTNmMHnyZF599dXosW+99RZ77rknzz77LLXXt3nyySeZP38+8+fPJy8vjwkTJkT3Rez56quv6NKlC5MmTUrKdZooGIYLfJ6wGx4M2aJUBmzbto3OnTvX2e71ehkzZgyrVq0CYNKkSVx00UWMGjUKgK5du3LHHXdw++23R4+ZNm0a11xzDX379mXmzJlxz5ebm8sdd9zBihUrmDdvXp39Y8eOjZ6zudjoI8NwgQgc5ZnN8uDIdJvS5vnDywv5evW2pPY5pFdHbhm/d8I25eXljBgxgoqKCtasWcO7775bp01FRQWzZs3i7rvvBmDhwoVceOGFNdqMHj2ahQsXRtu//fbbTJ48ma1btzJt2jQOPPDAuOf3er0MHz6cxYsXM3z48Oj2YDDIO++8w8SJExt1zfVhnoJhuMBb+gMP5v6D6/33pdsUI01EwjWLFy/m9ddf54ILLoiGe5YuXcqIESPo3r07PXv2ZNiwYa76fOWVVzj88MMpKCjg1FNP5YUXXiAYDNbbPja8FBGpHj16sG7dOo466qjmXaBDVnoKIjIeGD9w4MB0m2K0FQKVAAxkRZoNMRp6om8Jxo4dy8aNG9mwYQNQnVPYuHEjBx10EC+99BInnngiQ4YMYc6cOTVyAXPmzGHvvcPXMG3aNGbMmEGk2vOmTZt49913497gg8EgCxYsYK+99gKqRaqsrIxjjjmGSZMm8fOf/7zZ15aVnoKVuTBampDzVRG1IakGLF68mGAwSHFxcY3tXbt25fbbb+e2224D4KqrrmLq1KnMnTsXCN/0f/vb3/Kb3/yGbdu28dFHH7FixQqWL1/O8uXLmTRpEtOmTatzPr/fzw033ECfPn3qeCHt2rXjnnvu4e9//zuBQKDZ15aVomAYLU3Icds9mCi0VSLhmhEjRnDmmWfy6KOP4vV667Q76aSTKCsr46OPPqJnz5488cQTXHbZZQwePJgDDzyQSy65hPHjx/P8889zxBFHkJeXFz12woQJvPzyy1RWhj3Tc889l2HDhjF06FB27tzJiy++GNe2kSNHMmzYsLiC0lik9hCobGL06NFq6ykYLcEPi79kt+njWM0u9Lr1u3Sb0+ZYtGhRNGxiNI54n52IzFHV0fHam6dgGG7QcPJPzFMwWjkmCobhAg2FRcFD9nrWhuEGEwXDcEHIGSaoWP0do3VjomAYLgiF6h87bhitCRMFw3BDyHIJRtvARMEwXKBqnoLRNjBRMAwXWPjIAHjhhRcQERYvXhzd9tlnn3HooYcyaNAgRo4cyaWXXkpZWRlTp06lW7du0bkNI0aM4Ouvv26w1Ha6MVEwDDdERcESzW2ZadOmcfDBB0cnia1bt47TTz+dv/71r3zzzTd8+eWXHHvssWzfvh2AM888k7lz50Z/hgwZAtRfajsTMFEwDBeo5RTaPDt27GDGjBk8/PDDTJ8+HQiXxr7wwgsZO3ZstN1pp51G9+7dXfVZu9R2JpCVBfEMo6XRULimjM1SyAD+dz2sXZDcPnvsA8fdnrDJiy++yLHHHsuee+5JcXExc+bM4auvvqpTGjuWp59+mhkzZkTf114voXap7UwgYzwFEdlLRB4QkedE5Mp022MYsYTMU2jzTJs2jbPOOguAs846y1Wdodrho4KCAqDppbZbgpR6CiLyCHACsF5Vh8ZsPxa4G/ACD6nq7aq6CLhCRDzAY8D9qbTNMBqFjT7KHBp4ok8Fmzdv5t1332XBggWICMFgEBHhwgsvrFMa2w31ldrOBFLtKUwFjo3dICJeYBJwHDAEOFtEhjj7TgReBV5LsV2G0SgiOQUVSzS3RZ577jnOP/98fvjhB5YvX87KlSvp378/Rx55JI8++iizZs2Ktv3vf//LunXrXPVbu9R2JpBSUVDVD4HNtTaPAZao6jJVrQKmAxOc9i+p6nHAuam0yzAajQ1JbdNMmzaNk08+uca2U089lenTpzN9+nSuu+46Bg0axF577cUbb7xBhw4dgHBOIXZI6ieffFKn79hS25lAOhLNuwIrY96XAPuLyDjgFCCPBJ6CiFwOXA7Qt2/flBlpGLHY6KO2zXvvvVdnW+wqZ/Fu6BdddBEXXXRR3P6++uqr6GsRYd68ec03MklkzOgjVX0feN9FuynAFAivp5BaqwwjjGrzV7RKJhX+IPk5dRd4MYzmko7RR6uAPjHvezvbXCMi40VkSmlpaVINM4z6iOYUMmDy2tINOxh88+s8/2VJuk0xWiHpEIXPgT1EpL+I5AJnAS81pgNbo9locTJobebFa8KzZd/62l0y0zAaQ0pFQUSmATOBQSJSIiITNeyHXw28ASwCnlHVham0wzCaiwYzK3zUFsnmpYPTRVM+s5TmFFT17Hq2v0Yzhp2KyHhg/MCBA5vahWE0DrshpZX8/Hw2bdpEcXExYsOCXaGqbNq0ifz8/EYdlzGJ5sagqi8DL48ePfqydNtitA00A4ekNqRTZ0yeybF79+CSg/u3jEEppHfv3pSUlLBhw4Z0m5JV5Ofn07t370Ydk5WiYBgtjmbfcpyffb+Zz77f3CpEIScnh/79s/86soGMqX1kGJlMJi6yY1EUIxVkpSjYkFSjxQlmzuij/LK1LM8/h0E75yS13yXrt7N6a3lS+zSyj6wUBRuSarQ0mkHho85b5gJwcOnLSe33yH98yIG3v5vUPo3sIytFwTBanAwMHxlGKjBRMAw3ZNCQ1Igp6fdZjNZIVoqC5RSMFicDh6Q2xD9z/sXJnsyovGlkD1kpCpZTMFoa1Ujto/TjdtTRSd5PuCvX1qoyGkdWioJhtDhZ6CkYRlOwyWuG4YYMGn2UKj7MvYZvtC9wfLpNMdKIiYJhuMEJH0lGBJDCJNuWvp4N9MXKSLR1sjJ8ZIlmo8XRzFlPwTBSSVaKgiWajRbHySlkkqdgAmWkgqwUBcNocZzJAZl0G84kgTJaDyYKhuEGJ3zkyYQV2JwxqSYJRiowUTAMN0REgQwQBcNIISYKhuGGDBSFTAplGa2HrBQFG31ktDhOojmTRMHCR0YqyEpRsNFHRovjeApeMmdmsyWajVSQlaJgGC1Py4ePxt35Hk/NWpHUPl//ao0tpGMkxETBMNwQ9RRaThSWbyrjxucXJLXPK574glPu+ySpfRqtCxMFw3CBhDIn0dzcSWtrt1UkyRKjNWKiYBhuyKB5CjbqyEglJgqG4YqWDx8ZRjrISlGwIalGi5NB4SPDSCVZKQo2JNVoacQJG/nERMFo3WSlKBhGyxMjBqGWEYZ/5/yVYzyf1bvf5ikYqcBEwTDcEJtg1paZwHa4dx6Tc/8ZZ0/TU83v5v6S07wfNPl4o/VjomAYLpAYUdCgP2XnmfbZCs57aJartk3xEwZ41vK3nMlNONJoK9hynIbhBq2+BQeDgZR9cW74r4vJajYm1UghCf+2RWS+iz42qOqPkmSPYWQmMSGjUDDN9Y8slWCkkIbCR15gfIKfE4FuqTTQMDIBiUk0h0KZURRPFL5csYVT7/+EykDjbFpTGq5/tH57Bf2uf5WPl2xMhYlGFtKQKPxEVX9I8LMc+GkL2GkYaUVqhY8yA+Wm579izg9b+G7djoQt15bWLG0x9rZ3Afjih60APPrJ8lQYaGQhCUVBVWdEXotIrogME5F9RCQ3XhvDaLXUSDSnWRQamVN4ad5qDrjtHT5dtik19hitClejj0TkeGApcA/wL2CJiByXSsMMI5OQGjmF1M1T6MkmDvB83ZA1jepzzvLNACxes62JVhltCbeDKP4OHK6qSwBEZHfgVeB/qTIsESIyHhg/cODAdJzeaIPEThQLhVLnKbyTdx3tpBL4db1t1BLNRgpxO09he0QQHJYB21NgjyuszIXR4sSEj0IpDB+FBcE9kqThqaYzRoSGhqSe4rycLSKvAc8Q/vs5Hfg8xbYZRsZQc/JaekcfNVUI6rvxj/N8iddGlRsODYWPxse8Xgcc5rzeAOSnxCLDyEBiRSEYSt2M5giq7pbSKQjt5MeeT4GD620jcVSkmHCF4c4bZjM1907e3Ly8aYYarY6EoqCqF7eUIYaRycTOUwg2ck5AKvnVzrsYm/sp3249GXYdFb+RKn1kHTAkumlO/pXAOeRUbQGgm3916o01soKEOQURubyhDty0MYysJ3aeQiD1noLbbPIuofUAeALl9bbZZ/PrfJR3LT03119x1TAiNBQ+ul5EEk11FOAaYEryTDKMzKOGp5DuMheNpPfO8BDXTjuW1tmXrES10XpoSBQ+oGZeIR5vJckWw8hYYnMKgRRWSY2gmmA2Qpw7uTZ7nGr18Y/PXM7NLy5kwa1H0yE/p5n9GtmG5RQMwwU1ah+1SE7B3U3eVTraEZH4ulH3+H9/vByA9dsrXYnC8o07CakyoFv7hm0xMh5bT8EwXFBj9FG6y1zE0JzV195cuNa1h/HVqlKem1MSd9+4v73PEX+3hXtaCyYKhuECQfGrF2ghUUjBtOXaAnL543PitltdWjdpfcK9M7ju2XlJt8nIPEwUDMMFoiEqCYdSQi0w+siNJAjqLnyUoM94cxgq/GGvyMpptE1c1T4SkTzgVKBf7DGq+sfUmGUYmYUQooocoIJQoCqttjRGCKB5ISaj7eG2IN6LQCkwB2hccRbDaAWIKlWSB2xHW0QU3NzItZ7X8bHRp4Yb3IpCb1U9NqWWACJyEnA80BF4WFXfTPU5DcMNHkJUSQ4oaCDdz0WNu7031rMooIKhshw4tFHHGa0DtzmFT0Rkn6acQEQeEZH1IvJVre3Hisg3IrJERK4HUNUXVPUy4ArgzKaczzBSQwi/s7ZUqEXmKbh58m/c7f5rl+sp3JkzmWfz/oh3R3Xpi/6yhiM8X8Rt31s20FvWN8ISI5Nx6ykcDFwkIt8TDh8JoKo6zMWxUwkvzPNYZIOIeIFJwFFACfC5iLykqpHVRX7n7DeMjMCjISqlAGghTyGhKDQtp7BicxnETDvoQnyR2EtWAODxV49Cei/vV86rm+u0n5F3jfPKpjW1BtyKQpNXWVPVD0WkX63NY4AlqroMQESmAxNEZBFwO/A/VY3/WGIYacBLkCpxCgO3gKfgKkegzVsH4Yv8K/iygWevNxau5YsftnBDM85jZBeuwkeq+gPQiXDJi/FAJ2dbU9kVWBnzvsTZ9jPgSOA0Ebki3oEicrmIzBaR2Rs2bGiGCYbhHo8G8XvCoqDB5iean529kpItZfXuTxg+amTBIjdBptgWu3vWRG34yeNzmPzhskadz8hu3K7RfA3wJLCL8/OEiPws2cao6j2quq+qXqGqD9TTZoqqjlbV0d26dUu2CYYRFy9Bgt688JtmegoV/iC/fm4+Z07+NEGrhn0ApXGBpHhDU7XW7+awYXu6E/BGMnCbaJ4I7K+qv1fV3wMHAJc147yrgD4x73s72wwjI/FoiIDX8RSSNCR1w44EN9GYshr1EZ681jCJ5ik0dpjq4rX1J6uP+Nv7jezNyETcioIAsVXAgjRv2PPnwB4i0l9EcoGzgJfcHiwi40VkSmlpaTNMMAz3eAgS9IYTzS2SU0gQPoqEg2rc7FMw/Thej8f+8yM21SNm2yszpyaU0XTcisK/gVkicquI3Ap8Cjzs5kARmQbMBAaJSImITFTVAHA18AawCHhGVRe6NVpVX1bVy4uKitweYhjNwksIvLkEVSAJOYWG0FDDngJUC0QiTYi06Ztg2Gi8J7z6PIyyquxaT8JoHK5GH6nqP0TkfaoXgr1YVb90eezZ9Wx/DXjNTR+GkW68BBGvFz8+JFlrNCe4kYdchI/COQX3HsLFvjfqNcFdIMpoCyQUBRHpqKrbRKQLsNz5iezroqqbU2tevXaNB8YPHDgwHac32iBegiC+sCgkwVMY5/mSxST4+3Ux+Mht/NZNTmHj9irw1t3/R9+/OcX7UU3TTD9aNQ15Ck8BJxCueRT7pyDO+wEpsishqvoy8PLo0aObk+w2DNd4NQSe5HgKGvAzNfdOvtXehNNpcdokCh9JdU6hsSUsGoUqF/hsYcW2RkMrr53g/O7fMuYYRmbiJQgeHwHxNTunEAnVDGB1gjb1i0I0jwAu3YWmCcelj83m7bx456/fVVi8dhuDe3Rs0vmMzMDtPIV33GxrKWz0kdGihEJ4RVGPjwA5zR59FAm/+CTBjT/U2BhN/e0Xran/e7Jk/Y5Gnqd+7vQ9wNrSCuat3ArAjsoAKzbVP0HPyEwSioKI5Dv5hK4i0llEujg//QjPQE4LNvrIaFE0PNpGxENQmp9TcJPUTeQpRJ77JSa4r0Ag6G7EUixvfr0ugQ11PYzrfdMo2VzG71/8imAt4Trd9yGXTP2cCZM+JhAMcdr9n3Done812iYjvTTkKfyEcD5hsPM78vMi4SJ3htHqia605skhKD5o7ugjN5nahDmF8Nc2NoH81tfrGHjT/1i9te5Smi6NqnuaONuu8L3MHU+/yWMzf+CLFVvq7I/ohAKL125voi1GOkkoCqp6t5NPuE5VB6hqf+dnuKqaKBhtgkBkTWavl6Dk4Ak101NwNdzUXZmLCB99G56D8P3GnY2zpQn5hkhdpkRH2gil7MXt5LWQiHSKvHFCST9NjUkNYzkFoyUJBhxR8PgIeHLxNjt85KJNglCQxHkXb63lxjDa863rtpEzLVq7nVPu+7hZ5zUyD7eicJmqbo28UdUtNK/2UbOwnILRkgSdWkfi8eH35OPT5hV+c5VEdrGeQtwCd418Qo/00VnqJpwbkpmbX/iKL1ZsbdwJjYzHrSh4JeZRxFkkJzc1JhlGZhGM5hS8BLz55ISaWQ3UxZ07UaK5+m6d2hjNgR7XlWei9BOn7LbNkM5a3IrC68DTIvIjEfkRMM3ZZhitnmDQGX3k8RH05JGrFc3qz1W+IIE3kShUFK/vnlJ/4YFEOYU/5Dxa71H1Ucw28rES2tmM25XXfkt4JNKVzvu3gIdSYpFhZBhBZ16CeH0EvQXktoSnkCAZLS6K4MVytHdOvfum5PzDXScxeFB6sIm1FNfZ9+ecf7OXZwUVwQ38I+c+TvHOACz3l024XXktRHit5ZtU9TRVnayqaSuVaIlmoyWJJJrF4yPkyyeP5o4+ap4oRB7uY5/x99Ov+CbvAnyV9a93EA+PND7Mc37V03ya/zN6x6m6upcnvL4zGnIEwcg23M5oPhGYixMyEpERIuJ6/YNkY4lmoyWJzlPw+gj5CshrbqLZTRsX6ynEvr5cnidPArTfPL9ZtrnhEE/4HN2pO0+hGsspZCtucwq3AGOArQCqOhewekhGmyCSaPZ4vaivgHZS2ayB+M31FGpKQvpuvj5CUYEwWg9ucwp+VS2tleCyRwGjTRAKVYePNCe8+lqwqhxvXrsm9ScuREESJZo9dZPDUe+hBWeNTcn9O0USv7bRojXbGNlilhjJxK2nsFBEziE8NHUPEbkX+CSFdhlGxhCK5hRyEEcUKisaN3O4Jm6GpLoQDq0unZ0Oj6E+QQA4c/LMFrTESCZuReFnwN5AJeHhqNuAX6TIJsPIKEJOmQvx+pCcsHdQWdb06qK1b9+vLVjD5p01k9eJ1lNo7uxlw0iE29FHZap6k6ruB+wP/FW1mYO1m4GNPjJaklAwklPw4XFCRhXNEIXYEM+G7ZX89MkvuPyx2bWauKl4mp4Ibi6BBtvEei6bdti8hWzC7eijp0Sko4gUAguAr0Xk16k1rX5s9JHRksR6Crn5hQCUN0cUYm6YfqfGUcmW2tVNE93w3Tr4qSFR2CgetUtsG5mN27+uIaq6DTgJ+B/hkUfnp8oow8gkQs6MZo/XR05BWBQqy5sRPoq5SdYXCUo4+ihOmYt0JJrdknkWGYlwKwo5IpJDWBReUlU/9n9ttBFCMTOac/PbA1BZlty1AmonlhPd2yVBaYpM/FJaBiS7cCsKk4HlQCHwoYjsRjjZbBitHnXCR16Pj7z24ZBloLzpf/6xxe7qLVmRWBWcX7EeR2bJwUCpf/1pI7Nxm2i+R1V3VdUfa3jmzQrg8NSaZhiZQcRTwJdLQfvO4W3lzRjkEHPDF5SrvC+wq66t2STRymt1u4mGjzJFGl7NuzH6OlNsMtzhdvJaDRxhaHgIgmG0AtRZT8GXk0t+hy5AM0UhBs/O9fw65xlKAh8A58WctOH1FGoaGd7sZmJci+MvB/LTbYXhkvQOY2giNiTVaElCzkpr3pw8Cp3wERXNiJ7W8hSAOkX2Eq2nIHHCRzEdZhyeRhbpM9JLVoqCDUk1WpJYT8GXk8MOLcBT1ZxEc+zNPP6MZDfho/g9Z6AqWAApq3A7T+F0EengvP6diPxXREal1jTDyAxC/ogo5AGwQ9rhqUqOp0CcOkbhJglqH0V/a52NGRk+MrIKt57Czaq6XUQOBo4EHgbuT51ZhpE5aLDaUwAol0J8/uZMXqtLnVBQohnNMTf+2p5BJkqC6VR24VYUIgvqHA9MUdVXsTWajTaCOqOPIp5ChbeQnEBy5ilE7pd1wkeuJq/VXwhvyfrkzqNoHqYK2YRbUVglIpOBM4HXRCSvEccaRnYT8RTywqJQ5WtPbiA5tY9E6vsauVqKp949R/7jw8bZlELG32srsGUTbm/sZwBvAMeo6lagC5C22keG0ZKos8hObm54WGUwpz35waaXzo4XTqmTWYgphbGjMkBpub/uMaqo8xW2XIKRLNzOU+gJvKqqlSIyDhgGPJYqowwjo3A8hRxfTvhtXic6bGvOMMtYT6Ge0Ucx7/f901tUBkIsv/348L44OYV49ZAMoym49RT+AwRFZCAwBegDPJUyqwwjkwj5qVQfXm/46xJq15VOsgO/v6qBA+uj7pDUOlUuYoakVgbi5xeUODkJ0wSjmbgVhZCqBoBTgHtV9deEvQfDaP0Eq/DHOtWF3QAo3bSuSd3FJpHjDi8lcaK5+tj0rrxmtE7cioJfRM4GLgBecbblpMakFmD9Inj3LzZWznCFBP0EYkTB2yEsCjs2r2l235EwUd0hqS7/NqPhp8xlD8+qdJtgNAK3onAxMBb4i6p+LyL9gcdTZ1Ziml3m4usX4cM7YOak5BpmtEok5Ccg3uj7vKLuAJRvWVvfIQ1Q94Zf56aecPJabE4hfp9/9P27SZalgidyb0u3CUYjcFsl9WvgOmCBiAwFSlT1rym1LLE9zStzcehvYK8T4c3fwTf/S65xRqtDQn4CMY5xQaceAFSVNi18FDuyKJI0ThQ+GutZyFGe6uU6Y1tGRx/VOv4C31tNs81o87gtczEO+A6YBNwHfCsih6bOrBTj8cDJk6HncHhuIqxdkG6LjAwm7ClUh486desFQEXp+uZ37ohCHb8gxlOYlvsXHsz9R127YnIKHguFGknCbfjo78DRqnqYqh4KHAPclTqzWoDcdnD2dMgvgqfOgu1NfOozWj2eUM2cQpfiXQioh+CODU3qr/Yqa9C0RHP4kNqBp6YV0jOMCK6X41TVbyJvVPVbsjnRHKFjTzhnOpRvhulnO3XfDaMmEvITjPEUPF4vW6QT3h1NzClozQAQNC6nUMO7yOQMs5GVuBWFOSLykIiMc34eBGY3eFQ20HM4nPIgrPoCXrgSmliy2Gi9eEJ+glLzGWhzzi60K2/q6KM4noLW9hQaDgfFG5JqUSSjubgVhSuAr4GfOz9fA1emyqgWZ68T4MhbYeHz8MHt6bbGyDA8GqjhKQDsyOtBZ39T5ynEvAnVMyQ1QRio5o3f5ikYyaXBMhci4gXmqepgoG62q7Vw0DWw8Tv44K9QvAcMOz3dFhkZgjfkJ+ip6SlUFu5Kt+0fo6Eg4vHWc2R9xAsf1bqpu/JYFST+6CPDaCoNegqqGgS+EZG+LWBPi7BxRyXPf1lSc6MInHAX7HYQvHgVrPwsPcYZGYdX/YRqeQoU9SZP/GzbuLrxHcaJ8Xhq3dRDCe/x8VZui7fPMBqP2/BRZ2ChiLwjIi9FflJpWCp54P2lXPv0PF6dXysm7MuFMx6Hjr1g+jmw5Yf0GGhkFB4NEKqVU8jv1g+AdSVLG91frCZUz1Oo6RlIdAmT+pHoP4aRPNxWSb05pVa0MNcdM4i5K7fyy2fm0rNTPqP6dq7eWVgM5zwDDx0J086CS96A/I7pM9ZIO14NEKoVPurSa3cAStcsAcY1qr94oZ5GjT6KHKOh6nkKYolmIzkk9BREZKCIHKSqH8T+EF6JrSTRsZlMfo6XyefvS4+ifC57dDYrN5fVbNBtTzjjUdjwDfxnIoQafmozWi9e9dcRhR79hxBSIbDum3qOSkTdGc2eWp5C4tFHkWOCdZbjNIzm0lD46J9AvMLxpc6+rKW4fR6PXLQfgZBy8dTP6y5isvvh8OM74bs3w+UwjDaLTwOop6ZTnd+uA2s93cjdsqTxHbqYp+BmSKpXgzXX5jSMJNCQKHRX1To1IJxt/VJiUQuye7f2PHDevvywaSdXPjGHqtp16/ebCPtfCZ/eB58/nB4jjbSTo36Cnrw629fn9aNL2fdJOUck/BPFxYxmD0EsqWAkm4ZEoVOCfQVJtCNtjN29mNtPGcYnSzfxuxcW1H1CO+YvsMfR8NqvYck76THSSCs5VKHeuqJQVjSQXoESAv66S2UmJjZ8VM8COi48BQ2ap2Akn4ZEYbaIXFZ7o4hcCsxJpiEiMkBEHhaR55LZrxtO3bc3Pz9iIM/MLuH+D2qNJvF44dSHYZe94JkLrXheGyRX/aivrijk9RxCnvhZuezrRvWnicebRho12MQjIWp7CqrKtU/PbZQ9hhFLQ6LwC+BiEXlfRP7u/HwATASuaahzEXlERNaLyFe1th8rIt+IyBIRuR5AVZep6sQmXkezufaoPZkwohd3vP4Nr8yvNfY8vyOc+2z495NnQKktGtJmUCWfKogjCsV7jAFg4zczm9V/PTsSHBP+tatsquMpbNhRyfNf2t+n0XQSioKqrlPVA4E/AMudnz+o6lhVdVMNbCpwbOwGZ4b0JOA4YAhwtogMabTlSUZE+Oupwxi9W2d++cw85vywpWaDjr3CQ1Urt8OTp0NFExf4MbKKYKAqHO/35dfZ13vPkZRpHqGVjSsDVqNKqtaetOas2exi9FEYCx8ZycXtIjvvqeq9zs+7bjtX1Q+BzbU2jwGWOJ5BFTAdmOC2TxG5XERmi8jsDRuaVrq4PvJzvEy5YDQ9i/K5/LHZrNhUa6hqj6Fw5mOw8Rt45gIINjaWbGQbVRXO30AcUfDl5PJ97h503jo/6ecNBd0Og64pCrUL6xlGY3E7ozmZ7AqsjHlfAuwqIsUi8gAwUkRuqO9gVZ2iqqNVdXS3bt2SblyXwlz+7QxVvWjqZ2zZWVWzwe5HwPh7YNn78PI1NluolRMRBcmpKwoA27uOoF/VUnbu2O66z5rLadaqjur8DrrJOwBay1EQQuwvi1zbYhi1SYcoxEVVN6nqFaq6u6qmdVHXAd3a8+AFoynZUs7ERz+nvKrWU9vIc2HcDTD3yXABPaPVUlUZXmNDfPEH27UffDi5EmDJ7Ddd95lwqQTnyT/o2lOo+RUurFzP03l/cm2LYdQmHaKwCugT8763s801IjJeRKaUlqYurj+mfxfuPnMEX67cys+mfUkgWGs0yGG/hRHnwvu3wZdPpswOI734K8Ki4Mmpm2gG2H30sVRoDpWL3ItCrCrUlzsIJppFH3tMrUTz1s2b3NthGHFIhyh8DuwhIv1FJBc4C2hUcT1VfVlVLy8qKkqJgRGO26cnt47fm7cXrePmFxfW/AKLwPi7YcDh8PLPYanrVIuRRVRVhcNH3tz4nkJBYXu+zR9G9w0zXM0tqEOtY4ISLsMdDARcdlBTFBautgEQRvNIqSiIyDRgJjBIREpEZKKqBoCrgTeARcAzqrowlXY0hwsP7MeV43Zn2mcruOedWiUNvDlwxmPQbTA8fQGsnpsWG43UEaioAMCTGz+nAFAx4Gh2C5Xw/aLPXfWpNSav1dznd1a5lUBlguOrCUnNtRy8tkaz0UxSKgqqeraq9lTVHFXtraoPO9tfU9U9nfzBXxrbb0uEj2L5zTGDOGXUrtz19rdM/2xFzZ2ROQwFneGJU2FjE2rhGBlLoAFPAWD3cecRUA8bZjzR6P5rV0ytiix9HqhwdXyo1kzrXLERcUbzcFs6O6NQ1ZeBl0ePHl1ntnUqiMxh2Lijipte+IpuHfL40V7dqxt07AXnPw+PHAOPnwwT3whvM7KeYFX45uxLIArF3Xszr2Bfdlv9KgH/3/Hl5NTbFkiYafZLDihoAlGILY0RqpUAv9r3YuJzt2IqynZQ8t1ctq5YSKB0NezYgLd8I55QFSFPLsG8IjS/E54O3SnsNYju/feha4++iCdjxttkBFkpCukgx+vh/nNHcdaUT7nqqS94fOL+7NevS3WDrgPhvP/A1BPCwnDx/6Bdl/o7NLKCoD+caPYlCB8BhEacR49Pr+HLd59i5DEXJmxbc/JazXBPJHzkCdYfPqrRl6cBAWqFVFVWsGb5YjYt+5LK1V+Rt/kbupUtpVdoDQNjCgtWag5bpIgqTx45oSra6w46SPj/EydgvUMLWOPblY29j2S/8/6ELyc3DVeUWWSlKIjIeGD8wIEDW/S8hXk+/n3xfpwxeSaX/PtznrrsAPbpHZPs7jUCzp4WDiM9dQZc8CLkFraojUZyiXoKeYnrPw478lxWzfo/2s2+Hz36AsRlobpA5Y4a7yOeQqKcQiwedZuQTi8LbjucssLehIr64ivqRU77LuR16EJBx2LyC4uio6hCwSDl27dQsW0TVTu34N+5meCWEnzbVlBYVkKXqjV0003sJspuQFCF1Z6ebGg3kJIux5PXa2+69NuH4l4DaN+hEz1qeQEBfxUb1/zAhuVfU7ZmEWz8jval3zD2hwdY8LdZ9L3iOYo6d03DJ5Q5ZKUotHT4KJau7fN4YuL+nP7ATC54ZBZP/2Qse3bvUN2g/yFw2iPwzPnw9Plw9vTwMp9GVhJyRCE3r13Cdl5fDiWDLmb/xbcz/8MXGHbYyfW2jZ11XLlja419gciyn4k8hZjjPf6dCe3KFAoCpfTe/B2dN7uf5BchpMJG6cymnJ6sLBrF9x374iseQOf+w+m9xwj6tGtfY4x7Inw5ufTouwc9+u5BbCGFz56/lxFzb2H1vePYePJD7L7PAY22s7WQlaKQbnp1KuCpy8LCcN5Ds3j2irHsVhzjEex1QnjW80tXw/OXwykPgdc+6mwk5ISPcvITiwLAyJN+waq//psOH9xC8KAT8Prih3ZiMwrl22vW2ApK+AFC43gKGgoinpqjjXIqa1eRyUx6X/cR+QWF7Ni2hS3rSyjftonK7Zup2rGZYEWMt+TxkNOuiJz2Xchv34V2HYsp7tmXXfLbsUsK7Rtz8s9YuMsAer11BYXPHc+Xq//RYBiwtWJ3qiayW3EhT1y6P2dOnsm5jjD0LIoJMYw6P1w0782bwOODkyeHy3AbWUXI73gKLkQhN7+ANfvdyOjPrmHm039l7Ln1rdhXLQtVtTyFaI7AEaNYggE/vtyaf0P5VVvqtMtEcpzJf+07dqZ9x84NtE4Pex90PFsHfc6yB05i2CfX8NmOzYw59dp0m9XiZGXavaWHpNbHnt078Ngl+7O1zM+5D81iw/ZaT3cHXg0/ugUWPAsvXmVrPWcjgUj4yN2aUvseewHzC8Yw4tu7KfluboPtgztqFnVUJ7buq6obZikvL6uzrSiQHTOYvb7seP7s1LUHfX/xBgsLRjNmwa18et/lVFbU/dxbM1kpCi01o9kN+/Qu4pGL9mPN1grOmjKT9dtqDSU85Jdw+O9g3rTwzOeQTS7KJtQfFvq8fHeiIB4PPc57iErJpWL6JZTvrHtzj5357NleEn0dDAaj+YL8wNY6x23ftNY5vvpvqAcbXdlluKdd+yIGX/sKs7qewgHrn2blnQezcknbWVwrK0Uh0xjTvwuPXjKGNaUVnDXlU9aW1hKGw34drpX05RPwyi9MGLIICZRTqTn4GvGku8uuu7Hi0L8xILCMRZMvQhP8f/u2Vy/otGntiuhktnbBbXXabtu0us42IzXk5uWz/9X/5ssDJ9EtuJZujx/OZ8/fm26zWgQThSQxpn8XHrtkDOu3V3LmlJms3lorJjzuBjjkV/DFo/DyzyyUlCWIfyflEr8YXiKGHXE2n/a7glHb3mbOtD/U3BnjKXSoqL7Rb1m9lEi+oUMwHBrdqdXn3rFqUY3j15L80vGp4Nf+y9NtQpMZefR5lE/8kBU5/Rkx9xZWLWv9ZclNFJLI6H5deGziGDbvqOLMKTNZuTkmFikCR9xc7TE8exG4HItupA9PoIwKEk9cq48DLvg/Pi8cx6hv72bRh/GXHu8RqBaF0pKvESc0VKThsNMq767R/aEVs2ocuzGvd5PsamneLzg63SY0ix59BtL5kmcJ4mX9f65rxAJI2UlWikKmJJrjMapvZ564dH9Ky/ycev8nLFoTEwYQgcNvhGNug0UvwbSzoCo7xpm3VTyBciqlaaLg8XoY/JNHWertR+93f8bapfOA6npH6+lCocSEGktmR/d5nJm5m/P7RnfvsfkD/FWVRLyJsvbV+zKZnxy2e7pNaDbdevVj3sArGblzBp89+LN0m5NSslIUMinRHI/hfTrx3JUH4vUIZzwwk5lLa40QGftTmDApvHrb4ydDeXYMK2yL+AJlVHrcJZnj0aFjJ3LPe5oqfPifPIfynTui0aONuTWf9LttnQe1CuRVdtgt+roz23jzod9Fm4Q692/w/OtJf6mVAwYUp9uEpLD/ubcyq/gkDlj7JPPeezbd5qSMrBSFbGDP7h34z5UH0r0onwsf+YzXFqyp2WDkeXD6VFj1BTx8NGxelhY7jcT4ghX4PU3zFCLsNmAQK8fdQ59QCXOfuonok36HftE2Szz92S3wA4W1Esw5Rd1rvD9+/ZRo7aScrg2XedmYk/7CjEN3zcyHt8YiHg/DJ05iqbc//T+4hlXLMrbif7MwUUghvToV8NwVYxm6a0eueuoL7nt/Sc2FWIZMgAtegJ0b4MEfwQ8z02arEZ+cUDl+b9M9hQgjxp3M5x2PYmTJk2zfHB5aKsXVYZV1nUfhEaV3aDX+mDmlnrz2dfoqXxW+GRV23bXOvtpszknlPOC2R3679uSfOw0FAk+cyZYNaxo8JtswUUgxndrl8uSlB3D8Pj254/Vv+Pn0uTXXfO53MFz6Trii6mMnwrzp6TPWqENeqJxgEkQBoMfxN5EvfoKfTwVAO/eL7vP0Hk2VhsVgsycm5CPVX9El54YfGjqs+ih8jNdFhVRP8+tu9at4qtl9tCZ2HbAXJUc/RI/gWjZOHs/20uwoNeIWE4UWoCDXy71nj+TXxwzilfmrOe2BTyjZEjMyqXh3mPgW9Nkfnv8JvHIt+N0tsmKkljytIOhLjij0GTSSb3yDGbotfFPXmBu2N68d3+UOAmCbr7pKZ2y11YF7DGGl9GLvyrmuz5nrqX/thn0qHnLVR5dCK+hYm70P/DGLDv0X/fzLWDFpAgF/VbpNShpZKQqZPPqoPkSEqw4fyMMXjmbFpjJ+fPdHNfMM7bqEF+o56BqY/Qg8fJTlGTKAPK0glNNw3SO3bOp5SHRkEcBm6RR9va3rKADKPPWXW19TNCL6Ol517qDW3Bjy1p8P2Y6769q7V0emBI531TbCoZV3Nap9NjLiR2cxd+Qf2btqPp8/+ft0m5M0slIUMn30USKOGNydV35+MP27teenT37B9f+ZT1mVUxPfmwNH/TFcbnvrCph8WHhOQ1MWhDeSQoFWokkUhfb9Rtd4v9HXw3kl+JzEcWFV/aUrQl33TNj/b4NXUK6xT/bN/9s5cPeuTPzjk+xQ9x7TCu3ecKNWwH4nXc2cDkcw+vspLJn3cbrNSQpZKQrZzm7FhTx3xViuHLc7T89eydF3fcgH38YURht0HPzkQ+g+NFxI78nToLSk/g6NlBAKBMgTP5qTvIWSuvQfXv1GBL/XERyBgu4DAOgWXM92Iues+eSf33OvmHfVX9+Ahl9/oPuyf+W/YtrULwqf33SkK5uvOGwAXo8QcrlwUFtj4EUPsEWK8L34EyrKs3/ekYlCmsjxevjtsYN5+vKx5Pk8XPjIZ/xi+pesixTU67wbXPQqHHcn/PAJTDoAPrkXAq0ndpnplJc5w0OTKAo9etecW6C+cHgnGFQ69doDgCLZyRZv/PkF7bvGX04mSLik9k3H70Vxl5iVwxJ4md06uCvfEclrqCNQfrUS8LEUFXdn7bi/0S+0krlTf5Vuc5qNiUKaGdO/C69dcwg//9EevLpgDePufJ+/v/kN2yv84PHA/pfDlZ9A3wPgzd/BfQfAolcspNQCVDgVTj15yROFmms9CwEn5q+BCnbpXT1EdYe3c6RJDTp2jZl3EPPkHnBE4ZA9uvLeb34Uc4bqYnxXVzVuJu6d/jM4v+r66PuIKJxddRP/CR7CXf5TeTZwaNxjd+2UnOR8tjBs3KnMKj6JMWun892XH6bbnGZhopAB5Pm8/PKoPXnnl+M4ckh37n13CYfd+T7/evc7tpZVQZf+cN5zcO5z4SGKT58LDxwCC1+wiqsppGJneCBDvLkCyUCAoDMxTgLl5Obm8m63C/hk7BQCzjBYRfh51VVcWPVbADp17Vl9fIxgBJ2vcqQi6/+C+4XbxJTZ7tN/UPT1S8GxDdr3VPAIPgoNi74POedYpr34lf9K7g6eyq8DV9Q4ZmogXOfo4+uPaLD/1sZe5/+DTdIJeeUXBAPZsXZ2PEwUMoi+xe249+yRvHjVQeyzaxF/e/Nbxt72Lre8+BWL126DPY6Cn34KJz0AgXJ49kL412j4+B7YmR2LrWQTFc6qaDmFnZLar995qlcRQl4nhOMMQT7iqns58JgzUWe7AC+FDuKDUDgXkZsXfzRR5CleCc+BeU/CouDR6pvTiSPDk93mhgYw4IqnG233FZ5beSBwApvpUG+bWwMXMqDiiUb33Rro2KmYH/a7iYHBpXzxygPpNqfJmChkIMP7dOLRS8bw+i8O4fhhPXnqsxUc+8+POO7uj3jokxWs7DsBrvoMTnsECrvBWzfDPwbD0+fBguegsvGLoxt1qdgenpSU2z65y0eWERNacZZoDdUqpR7yOms1JwwTVrsKocgkt1C4/VljnNxFTL8e51yFOZ46pSdmh+qOaqqdVl4qfbk9cA5d2ycq+yFRj6ItMurYS/jOO5A+c/+ZtUnnrPzfy8Z5Ck1hcI+O/O304cy68Uj+cOLe5Po8/PnVRRxyx3sc+c8Z/Hn5Xrw99jG2XvgBjL4EVn4G/5kId+wOj02AD/8GKz+HoD/dl5KVVO0MFyrM75DconIVEhYFEQmv3w1osGa4IehMbAu5/L+LSIc6oiDecL+isaIQ/rpLnBFJv/Rf2eA5IuGqLoXhmdSXHVIzaf5U4HBXtrZmPF4vleN+Tw82MPe/f0+3OU0iOxZOrYWqvgy8PHr06MvSbUtL0KUwlwsP7MeFB/bj+407eW/xet77Zj2PzfyBh2Z8D0C/4mMZvusZHJK/jH3LZ9Br8yzy3v0T8Cfw5kH3IdBjGPQcBt0GQ+d+0KFXOJltxCW4cysA7TomVxQqPQU4UR4kIgq1PIVo+ChQzv3njqIgN86In5ikgkZyCk5iOdJvbE6h7rN/+EbeTipZp3W9oYlHj2Z4n7rb7z5rJB8v2cjEg/tz0/FD4Nbw9hsDl9a1sQ0y9JAJLJgxikHfTWHb1qvo2Cm7qsRmpSi0Zfp3LaT/wf255OD+lFcFWbCqlC9WbOGLH7Ywe0UpL25tBxwNHE2xbOPodt8xJncZe21Zzm7r/kvBF49G+wp5cvF36E2oqA+eDt3xduiOt8MuSPtdoLAr5HeCvA7hn9z24Z82JCLB8q0AtC9KsihIJHwkiBPS0VBNjyASPiJQyXH79KQh7vBM5Nrgvwnmh231+CKiUO2BVJfMqPYUbgxc5lhSd8DCGfv1rTFstWNBDht3VLFLhzwuPWRAHCtsHkOEgh//ic7PH8/MZ//E2Mv+mW5zGoWJQhZTkOtlTP8ujOlffdPaWRlg2YadfLd+O8s3lbG2dAj/La1gbWkFa8rK6Vi5jv6eNfSV9ewm6+mzeR29t6ykWBbSjVJ8kjhcUS7tqPAUUCX5BD05BCSXoCePgCeXkDePoPM75PzGl0fIm4d6CyA3H8kpwJPbDk9OAZ7cAnx5hXjzCsjJKyQnrx05BYXkt2tP+45d8OQ0fhnMpFK+lYB6KCxM7sz5kKe6kJ06olB7FFko6inUPy8lti7SuzKWpytHM8PpW5zfntjwkdfrbKsrAF/+/hi4o+Y2rRVmeuySMbz99TqK26f5/yULGDj8YOa8fQTDS55i4+pr6dprt4YPisFfVcmcZ/6P3kufpqNu49uOY+l3zl107ZH6hZVMFFoZhXk+9uldxD6949/IKvxBtpX72Vrup7Tcz9YyP9+WVbGzMkB5VZBgxTa85RvwlW1CqrbhqdqB178DX2AnucGd5IfKyAvuxKdV+EJV5GgVPq0iR3eSq1vIUT+5VJGLnzyqyMNPPlV4pfHzKnaSz3bpwE5PByp8HanMKSJQ2B3t2IecLn0o3KUfPXYbRFHX1KwZIJXb2C6FdPYm1zsKRUJGCCKOKNS+UUcSzSF3OYWIPkTy0p44OYV4nkKEooKGK6727tyOiw5KvLDPPWePZPIHS4HwfIZR8h2/brDn1kn3CX8i5/FDWfqf39P1Z482fIDDt198QJeXLuAAtvJV3ghWtxvFsM1vsurBEyn4xfsUduiUOqMxUWhz5Od4yc/xskvH5i0c44ZQSPGHQlQEFX9VJf6KMioryvBX7CBQWYa/ooxAVTnByp0Eq8oJVZURqionVLkTLd+KVGzFV7mVXH8p+YFSuu5cQvH2TylcV3Nt600UsTpvADuLBpHX/wB6jzqKbt2bv36xt2obO6WQ5I49AnWe4oOhIL7o3bymKERyAokGCcR6CtFZx879XqJeQV1RiBfkkXglLBqp4385eSgnDu/FicPDIj0ztDcz2bvNikLvgUOZ1W0Coza8SMmSr+g9cGiDx2zbugnfKz8nTyuZd9iDDBt3GuLxMP+959j7/UuZf/85DP/Vy1GvLxWYKBgpw+MR8jxe8nxAng86NH9msIZCbN2ykU2rl7F93feUrfsO34ZFdN7xHQPX/ZeC9dNhFizz7MaGnuPocsC5DBw6Jv5NrwFy/aWUe5I/cS0kkRFH/uh6CVJbFJwnfUJ1J0EF1INParb3REtROGs8x/MUiIhCaiY8nrt/40IkbYHdT/sD/vtfZd0LN9H7uhcTtq0o38mqSePZPbiSxeOmMPzw06L7hh1+Gp+uWcwB397JzIevZezl96TMZhMFI6sQj4dOxbvQqXgX2OeAGvuqqqr4Zv4MSr9+h4KSGexb8ji+/zzKshf6s2Gfyxlx3CXk1TP5Kx7tqjazPadrww0bSSR8FAr42dhxbwBWtRtco43UM1QVwiUtfLVu7EXtclmzrbL6xu+tm1MQT93w0bxbjm5gLoTRHLr26MvMvhcwduWDzH//Pwwbd2rcdhoKseCBi9nPv5A5Y/7OvjGCEGH/s25k1r++YezqR5n//mH19tVc2s5QEqPVk5uby6DRRzDmgr+wz40fsO2qr5i1142Ihth/7g1svn0on7/6ULQUREN0DG6mMj/5oqCOpyAhP6u6HcLYintZXHRIzUaRVdXieAqR4nexQ1IfvGA0vzt+L/oWh6uueiOiQJCjK//KQRV3o85xsemdooIcOrWruYjOuVU38H/+s+lsi+skhZHn/IHlnj70eP86SrfEL4v+6eM3s1/pG8zsezn7Hh9/aK94PIy4fDIl0oOiD2/BX1UZt11zMVEwWi1ddtmV/c/8Lf1+N5f5h05hp6cj+33+K7664yg2rl2R8NhQMEhnLSXYrlvS7YqMPtKgH4/AGuqOY48MVY0N/0QISt14co+i/BrDRCPhI48G+Vb7sIpueDzuwkcfh/ZhSnA8Oc1MsB8x2NaHBsgvKKTqhEl00a0se+jCOqu0ff7ifYz9/l/M7vAj9r/w9oR95eW3Y8PYm9kttJIvXkhNCCkrRaGtzGg2koN4PAw74kwG3DCLT/f8DQPL56MPHMriz96u95htm9fhkxDSPvmLxWhMEjnyrB+qFcKRSNmKOMNHI56CJPj6enx1w0fR/EUTbG4KU87fl8V/OraFzpbZ7DnqMGYP+hUjd85g/t2nsvDjV9m4dgWf3ncZo764kYW5w9nnqidcJZBHHHkOn+1zK0OPS83c3awUhWxeec1IHx5fDgeccxPrz3yFSsmj/6tn8cU7z8Rtu3VDeFGjnE494u5vDpHwEUF/9Om9dlhfPPWLQqRMdo32td5Hbi5eYstcpDbRXBuf10N+jq29EOGAc37HzP5XMWL7R+z91jl0fWAfxqx7ltldJ9D3qhfIy3e3wp94PIw59dqUDU21RLPR5thtyP5s+un7rHzgeIZ+eCVf5+Yw5JCTa7TZtjY81r5d1+RPFop6CqFAdFRUbU9Bo55C3SRwSLygsYnjODjn8NYYfVR/7SOAddqJF4IHuboGo2mMvfD/KN10DcvnfUDZii8oHnYc+486LN1m1cBEwWiTFHfrie+q11nxr6Pp+/aVLO3Sm9333i+6v3zNdwB07zck6eeOhIZCoRCnjerN/JVb+dXRg2q02dopPKb92/b7Unvlg3ieQm0ielHj9t9AiZIvT/+UcV3b4/1uA+8sWt/gOSLsV3EfIYQ5ro9o2xQVd2f4EWcAZ6TblLhkZfjIMJJBUZddaH/xc1RIHvnPncvm9auj+2TLUkq1kC5dUxE+qg4NFeR6ufP04XSpNdLH23skQyseYkOfH9c5PhQjCn/1nxW2t848jLr5Ay0IT8N723NwXLuOHdqTQT06cOkhA5h2+QFx28RjA53YhIVyWwsmCkabpkefgWw58VG6hjaz+qGzoyNDikoXszp3t+rYfjKJzj6uP7Y/bs9u/O28Q7jmyD3q7As5o49CwQD3B0+kX8VTdU8Rkz+YfP6+TLvsAMgvYkjFI0z2np2MqzBaKSYKRptnj1HjmD/yVoZWzWX2g1dRvqOU/lXfsaXrvik53xe7ns+nob1Y1K2uFxBBRDh2aI+4w0KDUVGIzRfUPj4yQkk5Zu8ejN29GI9AGfltehEco2Hsr8MwgP1OuppPu53BAeufYcXdR5MrQToOPzEl59qZ352zqm6mzNe0kEskfKRxJrZFidY5qs4qdMgPD1M9YVjDpbiNtouJgmE47Hv5JOYUHsIg/2JmtzuUvcccmZLzROsUNbG8RCR8FAz4OWpI/HkU0eJ3Medon+dj3i1Hc8OP92rSeY22gY0+MgyHnJxcRv3qJdavX82+u/RKTT4B2KtneOH7PXt0aF5Hqtx79kg27ayKzkGIUF+ZbDclso22jYmCYcQgHg+79Gh+2e1EHL13D97+5aEM3KVpoqBRB1/Jz/Gya6eCOm0Slck2jESYKBhGGmiqIAAEIqW3E4Wf4uQUUsXw3kX4krwQkZE+TBQMI8t4qf+tzPz6YQ7uPbr+Ru27s0Xb80/fRdyWYntevDr+vAcjOzFRMIws49rTj+DTZcMZ3Kv+NeGKO7bn1F2e4ZdH7dmClhmtARMFw8gy8nO8jBuUuCy1z+vhxausjpHReDJGFESkELgPqALeV9Un02ySYRhGmyOl2SEReURE1ovIV7W2Hysi34jIEhG53tl8CvCcql4GpGbWkGEYhpGQVA8ZmArUWGVDwvPvJwHHAUOAs0VkCNAbWOk0q7vclGEYhpFyUioKqvohsLnW5jHAElVdpqpVwHRgAlBCWBgS2iUil4vIbBGZvWHDhlSYbRiG0WZJx+DiXan2CCAsBrsC/wVOFZH7gZfrO1hVp6jqaFUd3a1b8tfPNQzDaMtkTKJZVXcCF6fbDsMwjLZMOjyFVUCfmPe9nW2uEZHxIjKltLQ0qYYZhmG0ddIhCp8De4hIfxHJBc4CXmpMB6r6sqpeXlRkqz0ZhmEkE2lq+V5XnYtMA8YBXYF1wC2q+rCI/Bj4J+AFHlHVvzSx/w3AD8mxNiV0BTam24hm0hquAVrHdbSGa4DWcR3Zfg27qWrcpGxKRaGtIyKzVTVBgZrMpzVcA7SO62gN1wCt4zpawzXUh5U2NAzDMKKYKBiGYRhRTBRSy5R0G5AEWsM1QOu4jtZwDdA6rqM1XENcLKdgGIZhRDFPwTAMw4hiomAYhmFEMVEwDMMwopgopBERKXQqvp6QbluagoicJCIPisjTInJ0uu1pDM5n/6hj/7nptqcpZPPnX5tW8F3wiMhfROReEbkw3fY0BxOFJtDIxYMS8VvgmdRYmZhkXIOqvuAsinQFcGYq7XVDa1jUqTHXkGmffyxN+PtK23ehPhp5DRMI13HzE678nLWYKDSNqbhcPEhE9hGRV2r97CIiRwFfA+tb2niHqTTzGmIO/Z1zXLqZSvYv6jQV99cQIVM+/1im4v7vK93fhfqYivv/i0HAJ6r6S+DKFrYzqWRM6exsQlU/FJF+tTZHFw8CEJHpwARVvQ2o4xKLyDigkPAfVrmIvKaqoVTaHUuSrkGA24H/qeoXKTa5QRpzTVQv6jSXDHo4asw1iMgiMujzj6WR/xftSeN3oT4aeQ0rCa8vD5n1kNFoTBSSR7zFg/avr7Gq3gQgIhcBGzPhS0AjrwH4GXAkUCQiA1X1gVQa10Tqu6Z7gH+JyPEkWNQpQ6jvGrLh848l7nWo6tWQcd+F+qjv/+Ju4F4ROQT4MB2GJQsThTSjqlPTbUNTUdV7CN9cs47WsKhTNn/+8cjy70IZMDHddiSDjHGbWwHNXjwoA2gN11Cb1nBNreEaoHVcR2u4hoSYKCSPZi8elAG0hmuoTWu4ptZwDdA6rqM1XENCTBSagIQXD5oJDBKREhGZqKoB4GrgDWAR8IyqLkynnYloDddQm9ZwTa3hGqB1XEdruIamYAXxDMMwjCjmKRiGYRhRTBQMwzCMKCYKhmEYRhQTBcMwDCOKiYJhGIYRxUTBMAzDiGKiYLQJRCQoInNjftyUNk85MXb1StDmFhG5rda2EU5BPETkPRHZISKjU22v0fqxeQpGm0BEdqhq+yT36XMmMzWnjwbtEpE9gddVdUDMttuBMlX9o/P+feA6VZ3dHHsMwzwFo00jIstF5A8i8oWILBCRwc72QmeRlc9E5EsRmeBsv0hEXhKRd4F3RKSdiDwjIl+LyPMiMktERovIJSLyz5jzXCYid7mw52gRmenY86yItFfVb4EtIhJbsfYMYFpSPwzDwETBaDsU1Aofxa5UtlFVRwH3A9c5224C3lXVMcDhwJ0iUujsGwWcpqqHAT8FtqjqEOBmYF+nzTPAeBHJcd5fDDySyEAR6Up4wZwjHXtmA790dk8jXGcHETkA2Kyq3zX+YzCMxFjpbKOtUK6qI+rZ91/n9xzCy3QCHA2cKCIRkcgH+jqv31LVzc7rgwnX0kdVvxKR+c7rHY43cYIT+89R1QUN2HgA4YVmPg6vX0Qu4do7AE8Dn4jIrwiLg3kJRkowUTAMqHR+B6n+Tghwqqp+E9vQCeHsdNnvQ8CNwGLg3y7aC2HBObv2DlVdKSLfA4cBpwJjXdpgGI3CwkeGEZ83gJ85S44iIiPrafcx4fg+El6rd5/IDlWdRbj2/jm4e7L/FDhIRAY6/RU6SeYI04C7gGWqmtWLwxuZi4mC0VaonVO4vYH2fwJygPkistB5H4/7gG4i8jXwZ2AhUBqz/xngY1Xd0pCBqroBuAiY5oShZgKDY5o8C+yNhY6MFGJDUg2jGYiIl3C+oEJEdgfeBgapapWz/xXgLlV9p57jkzJU1oakGsnCPAXDaB7tgBkiMg94HvipqlaJSCcR+ZZwgjuuIDhsa2jyWkOIyHvAAMDf1D4MI4J5CoZhGEYU8xQMwzCMKCYKhmEYRhQTBcMwDCOKiYJhGIYRxUTBMAzDiGKiYBiGYUT5f+JzhQBWGPuOAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# Plot cross sections at 293.6K\n",
"plt.loglog(broadr_xs[0].x, broadr_xs[0].y, label='BROADR')\n",
"plt.loglog(acer_xs[0].x, acer_xs[0].y, label='ACER')\n",
"plt.title('293.6K cross section')\n",
"plt.xlabel('Energy [eV]')\n",
"plt.ylabel('Cross section [b]')\n",
"plt.legend()"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "527de7a7",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x14b12d219400>"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEaCAYAAAD+E0veAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAABF3ElEQVR4nO3dd3zU9f3A8df7LjuEMIJsmbJEGUaGWke1SgfSOip2iNbdau2wjk5bf62jU61WbbVoraC17tlW21oVEFA2iIiMsFcCZN54//64710u45JvkrvcyPv5eORB7rvu/Q3J932fLaqKMcYYA+BJdgDGGGNShyUFY4wxEZYUjDHGRFhSMMYYE2FJwRhjTIQlBWOMMRGWFIzpYkTk+yLyp2THYVKTJQWTVCIyW0TWikiliHwkIp+I2ne6iKwTkSoR+beIDInalysiD4vIQRHZKSLfidp3qoiURb3OEZGnReRtEeneeXeXfI1/FgCq+gtVvSxZMZnUZknBJI2IfAq4A7gEKAJOBjY6+0qAp4EfAb2AJcATUaffAhwFDAFOA24QkRnNvEeuc50ewJmqerAD8Wa191xj0oUlBZNMPwV+pqoLVTWoqttUdZuz7xxgtar+TVVrCCWBCSIyxtk/B7hVVQ+o6lrgj8DF0RcXkQLgBSAL+KyqVjYXhIjki8ivRWSziFSIyFvOtqEioiJyqYhsAd4QEY+I/NA5dreIPCoixc518kTkMRHZJyLlIrJYRPo6+y4WkY0ickhEPhaRL8eIZYqILHFKQLtE5DdR+6aJyDvOtZeLyKlR+3qJyJ9FZLuIHBCRZ0WkEHgFGCAih52vASJyi4g8FnXu2SKy2rnuf0RkbNS+TSJyvYiscH42T4hIXov/qyatWVIwSSEiXqAU6CMiG0SkTER+LyL5ziFHA8vDxzsP9I+Ao0WkJ9A/er/z/dFRr3MJPRBrgFmqWt1COL8CjgNOIFQquQEIRu0/BRgLnEUo8VxMqHQyHOgG/N45bg5QDAwGegNXAdXOw/lu4NOqWuS8z7IYsdwF3KWq3YERwJMAIjIQeAn4PyfG64G/i0gf57y/AAXOz+AI4LfOz+zTwHZV7eZ8bY9+MxEZBcwDvgX0AV4GXhCRnKjDvgjMAIYBx9Io+ZrMYknBJEtfIBs4D/gEMBGYBPzQ2d8NqGh0TgWhaqZuUa8b7wsrAqYDj6hqbawgRMQDfA24zimpBFT1nUbn3KKqlU5i+TLwG1XdqKqHgZuB2U7Vko9QMhjpXGdpVHVVEBgvIvmqukNVV8cIyQeMFJESVT2sqgud7V8BXlbVl51S1T8JVal9RkT6E3r4X+WUnHyq+t9Y99zIBcBLqvpPVfURSpD5hBJX2N2qul1V9xMqeU10eW2ThiwpmGQJf3K/x3lI7gV+A3zG2X4YaNwo3B045Oyj0f7wvrC9wGzgERE5q4U4SoA8QqWQWLZGfT8A2Bz1ejOh6qm+hD6tvwbMd6px7hSRbOcT+wWESg47ROSlqGqwxi4FRgHrnOqnzznbhwDnO1U85SJSDpxEqMQ0GNivqgdauIdYGtyPqgad+x0YdczOqO+rqE/KJgNZUjBJ4TzAyoDoaXqjv18NTAi/cKpgRhBqZzgA7Ije73zf4NO3qj4NXA48JSKnxQhlL6EqphEthRv1/XZCD+iwIwE/sMv5hP5TVR1H6JP254CLnFheU9VPEXqIryPUBtL0jVQ/VNULCVUB3eHEXkjoQf0XVe0R9VWoqrc7+3qJSI9WYm9Og/sRESGUZLbFPMNkNEsKJpn+DFwrIkc47QTfBl509j1DqLrlXKdh88fAClVd5+x/FPihiPR0PnVfDsxt/AaqOg+4BnhORE5sZn8QeBj4jdMI6xWR6U6vpebMA74tIsNEpBvwC+AJVfWLyGkicozTXnKQUFVQUET6isgs5+FeS6ikE2zu4iLyFRHp48RV7mwOAo8BM0XkLCfGPAl1Nx2kqjsItZ/c5/w8skXkZOfcXUDvcGN4M54EPiuh7r/ZwHedGN+JcbzJcJYUTDLdCiwG1gNrgfeBnwOo6h7gXOf1AWAqoeqgsJ8QqvLZDPwX+KWqvtrcm6jqI4Qedi+JyJRmDrkeWOnEsp/QJ/RYfxsPE6omehP4mFAp41pnXz/gKUIJYa0T11+ca32H0Kfy/YQarq+Ocf0ZwGoROUyo0Xm2qlar6lZgFvB9YA+h0sH3ouL8KqEktA7YTajhGCeJzgM2OtVOAxr9bD4g1F5xD6FS00xgpqrWxYjPZDixRXaMMcaEWUnBGGNMhCUFY4wxEZYUjDHGRFhSMMYYE2FJwRhjTERaz/pYUlKiQ4cOTXYYxhiTVpYuXbpXVfs0ty+tk8LQoUNZsmRJssMwxpi0IiKbY+2z6iNjjDERlhSMMcZEWFIwxhgTkdZtCsaYrsHn81FWVkZNTU2yQ0kreXl5DBo0iOzsbNfnWFIwxqS8srIyioqKGDp0KKHZvU1rVJV9+/ZRVlbGsGHDXJ9n1UfGmJRXU1ND7969LSG0gYjQu3fvNpeuLCkY49LmtUvx1cVc2dMkmCWEtmvPz8ySgjEu7N25hSFPfJL377802aGYJPF6vUycOJEJEyYwefJk3nkntA7Rpk2byM/PZ+LEiYwbN46LLroIn88XOe+tt95iypQpjBkzhjFjxvDggw82uK7f76dPnz7cdNNNDbafeuqpjB49mmOPPZYxY8ZwzTXXUF5e3iSe8ePHM3PmzAb7OiItk4KIzBSRBysqGq/rbkxi1FZVAjDowKIkR2KSJT8/n2XLlrF8+XJuu+02br755si+ESNGsGzZMlauXElZWRlPPvkkADt37uRLX/oS999/P+vWreOtt97igQce4KWXXoqc+89//pNRo0bxt7/9jcbr2/z1r39lxYoVrFixgtzcXGbNmtUknlWrVtGrVy/uvffeuNxnWiYFVX1BVa8oLo61wqAx8ZWVkxP6V/1JjsSkgoMHD9KzZ88m271eL1OmTGHbttAS1/feey8XX3wxkydPBqCkpIQ777yT22+/PXLOvHnzuO666zjyyCNZsGBBs++Xk5PDnXfeyZYtW1i+fHmT/dOnT4+8Z0dZ7yNjXLAVClPHT19YzZrtB+N6zXEDuvOTmUe3eEx1dTUTJ06kpqaGHTt28MYbbzQ5pqamhkWLFnHXXXcBsHr1aubMmdPgmNLSUlavXh05/l//+hcPPPAA5eXlzJs3jxNOOKHZ9/d6vUyYMIF169YxYcKEyPZAIMDrr7/OpZfGp2ozLUsKxnQ2DQYAECw5dFXh6pp169bx6quvctFFF0U+LHz00UdMnDiRvn370r9/f4499lhX13zxxRc57bTTyM/P59xzz+XZZ58lEAjEPD76w0k4SfXr149du3bxqU99qmM36LCSgjEuBANBwJJCKmjtE31nmD59Onv37mXPnj1AfZvC3r17OfHEE3n++ec5++yzGTduHEuXLm3QFrB06VKOPjp0D/PmzeOtt94iPNvzvn37eOONN5p9wAcCAVauXMnYsWOB+iRVVVXFWWedxb333ss3v/nNDt+blRSMcSFcUjAGYN26dQQCAXr37t1ge0lJCbfffju33XYbAN/4xjeYO3cuy5YtA0IP/RtvvJEbbriBgwcP8r///Y8tW7awadMmNm3axL333su8efOavJ/P5+Pmm29m8ODBTUohBQUF3H333fz617/G7+94m5clBWNcCAatgbmrC1fXTJw4kQsuuIBHHnkEr9fb5LjPf/7zVFVV8b///Y/+/fvz2GOPcfnllzNmzBhOOOEEvva1rzFz5kyeeeYZPvnJT5Kbmxs5d9asWbzwwgvU1obGw3z5y1/m2GOPZfz48VRWVvLcc881G9ukSZM49thjm00obSXp3IBWWlqqtp6C6Qyb173HkPmnsY9iet+yJdnhdDlr166NVJuYtmnuZyciS1W1tLnjraRgjAuqoTYFj7UpmAxnScEYF4It9AgxJpNYUjDGhfqGZispmMxmScEYF4LW+8h0EZYUjHFBg+E2hWCSIzEmsSwpGOOCdUk1XYUlBWPcCFoJwcCzzz6LiLBu3brItnfffZeTTz6Z0aNHM2nSJC677DKqqqqYO3cuffr0iYxtmDhxImvWrGl1qu1kS5mkICJjReR+EXlKRK5OdjzGRAt3SbVlXrq2efPmcdJJJ0UGie3atYvzzz+fO+64gw8++ID333+fGTNmcOjQIQAuuOACli1bFvkaN24cEHuq7VSQ0KQgIg+LyG4RWdVo+wwR+UBENojITQCqulZVrwK+CJyYyLiMaSvrkmoOHz7MW2+9xUMPPcT8+fOB0NTYc+bMYfr06ZHjzjvvPPr27evqmo2n2k4FiZ4Qby7we+DR8AYR8QL3Ap8CyoDFIvK8qq4RkbOBq4G/JDguY9pE1ZJCynjlJti5Mr7X7HcMfPr2Fg957rnnmDFjBqNGjaJ3794sXbqUVatWNZkaO9oTTzzBW2+9FXndeL2ExlNtp4KElhRU9U1gf6PNU4ANqrpRVeuA+cAs5/jnVfXTwJdjXVNErhCRJSKyJDxDoTGJpgFrU+jq5s2bx+zZswGYPXu2q3mGGlcf5efnA+2farszJGPq7IHA1qjXZcBUETkVOAfIBV6OdbKqPgg8CKG5jxIWpTHRrKSQOlr5RJ8I+/fv54033mDlypWICIFAABFhzpw5TabGdiPWVNupIGUamlX1P6r6TVW9UlXjs9ioMXGikd5H9jmkK3rqqaf46le/yubNm9m0aRNbt25l2LBhnHHGGTzyyCMsWlS/dvfTTz/Nrl27XF238VTbqSAZSWEbMDjq9SBnm2siMlNEHqyoqIhrYMbEorY2c5c2b948vvCFLzTYdu655zJ//nzmz5/P9ddfz+jRoxk7diyvvfYaRUVFQKhNIbpL6jvvvNPk2tFTbaeChE+dLSJDgRdVdbzzOgtYD5xOKBksBr6kqqvbem2bOtt0lmX/msfEt67iIAV0v2VHssPpcmzq7PZLqamzRWQesAAYLSJlInKphj5yXQO8BqwFnmxPQjCmM9nKa6arSGhDs6peGGP7y7TQmNwaEZkJzBw5cmR7L2FM21hDs+kiUqahuS1U9QVVvaK4uDjZoZguItzQLGm8UqExbqRlUjCms1n1UfKl89LBydKen5klBWNcCM99lCpqa6qSHUKnysvLY9++fZYY2kBV2bdvH3l5eW06LxmD1zrM2hRMp0uhksKW9cs48vFTWDL5dkrP7hpzRw4aNIiysjJsFoO2ycvLY9CgQW06Jy2Tgqq+ALxQWlp6ebJjMV1DpE0hyXEA7NnwPkcCnvUvE5oqLPNlZ2czbNiwZIfRJVj1kTEupFr1kTGJYknBGDdSqPoowqrXTQKkZVKwaS5Mp4ssspM+T+Lld3yKRU/+MtlhmDSTlknBximYzpaOXVInVL/L1DX/l+wwTJpJy6RgTKdLoTYFSYXWbpOxLCkY40YaVh+11Ue3TmLJb85NdhgmydKyS6oxnS01q4/im6BGBDbCwY1xvaZJP2lZUrCGZtPpUqj6yJhESsukYA3NptOl5Cyp1rhg4i8tk4IxnS6YiiWFzG3fMMljScEYF1JqRLN1PzIJZEnBGDdSsvrImPizpGCMGzZls+ki0jIpWO8j0+mcLqmZPE7BGEjTpGC9j0xnE6dNwZNCScESlEmEtEwKxnQ2jYxoTqEGZ2MSwJKCMW44Dc3eFEoKauMUTAJYUjDGDaek4BWNrMKWbFZ9ZBLBkoIxLkjU3EfJXjxebJyCSSBLCsa4Up8IAgF/EuOIYt1kTQJYUjDGjaiSQsokBWMSIC2Tgo1TMJ1Noqa5SJU2BZvuwiRCWiYFG6dgOl/nVx99+P6bVOzfE/sAqz4yCZCWScGYTqfR1UedMw/SUc/NZN/vz2gaSgf+bNcsfJUDe3Z0JCyT4SwpGONCdPURnbgK2/DgpqaxdGCsxLhXL6D8D2d1ICKT6SwpGONGBjU0DwtuTnYIJoVZUjDGheiBYsEkr9fckeojY1pjv13GuJGKvY9sRLNJAEsKxrggmjnVR8a0JKulnSKywsU19qjq6XGKx5jUFNX9M9hJvY9aZ+MUTPy1mBQAL/CZFvYL8Hz8wjEmNUWXFKz6yGSy1pLClaraYlcFEfl6HONxRURmAjNHjhzZ2W9tuqjobqDBYOKqj97/x2PUffAPpl77aOxYbCSzSaAW2xRU9a3w9yKSIyLHisgxIpLT3DGdxUY0m07XSdVHk975BlP3PZew6xvTmtZKCgCIyGeB+4GPCFUZDRORK1X1lUQGZ0yqiK4+6szBa8Z0Nre9j34NnKaqp6rqKcBpwG8TF5YxqaXhOIXM6X3k99Wx8A9XsX/3tmSHYlKE26RwSFU3RL3eCBxKQDzGpKTokkLq9D6CbRvXsvDey9oc0/r3/gvAqv88xbRd89j06NWJCM+koRaTgoicIyLnAEtE5GURuVhE5gAvAIs7JUJjUoBEtSmopkrvI6j+65eZtudvbFy1MOYxfl8dS1/+c4NeU6OePxuoH50tGVT6MR3TWpvCzKjvdwGnON/vAfISEpExKaizSwoaDLoaheBmnebF825l+sa7WaJ+SjsemslwLSYFVb2kswIxJpV1VpfURPAcDLUX+A+2sDaDjXkwjtaqj65o7QJujjEm3Xka9D5Kneqjeu17qNuYB9NYa9VHN4nI3hb2C3Ad8GD8QjIm9TSoPuqELqmq6qr6yN1RxrjXWlL4Lw3bFZrzzzjFYkzK8miQOvWSIwE0hXofuWlTaInakp6mEWtTMMYFIYifLHIIoNo5JYUEXDUB1zSZxqbONsYFjwbwS+gzlKbIiGZBO1x9ZG0KpjFLCsa4IATxOQXrZFcfibT1z9ZKCMY9SwrGuODVAP5wUtDEd0l1NUCuwYA6e/Cb+HA7IV4ucC4wNPocVf1ZYsIyJrUIQXySAwpBf5LHKbS5yqdtx9dUHWbLuqWMmnxK6webjOO2pPAcMAvwA5VRX3ElIp8XkT+KyBMicma8r29Me3k0SJ3kAhD01yY5mpDonkfxbBtYdf8cRj1/Nnu2b4ps27H5A1a//VKzx+/duZW9O7fE7f1NcrkqKQCDVHVGe95ARB4GPgfsVtXxUdtnAHcRWt3tT6p6u6o+CzwrIj2BXwH/aM97GhNvHgL4PLkQBPXXJfz9WqwOaiYBtFx91Py+murmP9f1PbQGgOrD5ZFtRzw8lf6icGJFk+NL7nf+rG9pus+kH7clhXdE5Jh2vsdcoEFCEREvcC/waWAccKGIjIs65IfOfmNSgkeD+D2hkkLSk0L9UR16j7w7Brg+1ivWZtFVuE0KJwFLReQDEVkhIitFZIWbE1X1TWB/o81TgA2qulFV64D5wCwJuQN4RVXfa+56InKFiCwRkSV79rQ0l4sx8eMhiN8bmgMy2ClJIXZDc9tritp2wmDd7sTQ1vcxmcBt9dGn4/y+A4GtUa/LgKnAtcAZQLGIjFTV+xufqKoP4kyrUVpaar+2plN4CRDwhJKCBlKlpNDmi8b/mibjuEoKqrpZRCYAn3A2/U9Vl8c7GFW9G7g73tc1pqM8BAl6Q9VHdEZScDPpnuuHvCUD456r6iMRuQ74K3CE8/WYiFzbgffdBgyOej3I2eaKiMwUkQcrKqxhy3QOD0ECWfmhF35fwt+v5XEKoT/bttci2ehl0zq3bQqXAlNV9ceq+mNgGnB5B953MXCUiAwTkRxgNvC825NV9QVVvaK4uLgDIRjjnocg6pQUNJD4Lqltrj5ycbzWVbXxmqk4RbhJNLdJQYDosf0BXH5QEZF5wAJgtIiUicilGhoSeg3wGrAWeFJVV7sP25jO5dUgeLzUaVanVB8FXa7Z4G7uo9Ax0z/+fQciMl2F24bmPwOLROQZ5/XngYfcnKiqF8bY/jLwssv3b0BEZgIzR44c2Z7TjWkzD0EQLz6ykEDHq48q9u+hsKiYrOycZve7G6egLqfObn+bwqr/Pcfhj5cwrd1XMOnGVUlBVX8DXEKoa+l+4BJV/V0C42otHqs+Mp3KSwD1ePFJFtLBkkLA76f47pG8d++cmMe01NAs1CeFhFJl/OsXMW2j9f3oSlosKYhId1U9KCK9gE3OV3hfL1VtPP7AmIzkdUoKfrIg2LGSQiDgxwtMOvBa7INaqs+PajCO18prHV2sx2SO1koKjzv/LgWWRH2FXxuT8TQYDI3o9Xjxk42ng0khPE+Rl9gP/viOU+icXkdrFr7aKe9jEqvFpKCqn3P+Haaqw6O+hqnq8M4JsSnrkmo6U6TRV7z4JBvpYFIIP/A9LUwd4WYdaHE9dbaLnknNJY5murB+tHJhzGuMe/WCVt/HpD634xRed7Ots1ibgulMgYAzVbbHQ0Cy8AQ71qbgaq2EFnTaamnNJJoRfz+L/btdDykyaai1NoU8oAAocWYuDf82dic0VYUxGS8YSQrZ+KXj1UdutHVCvI4mira0KdRUHurQe5nU1lqX1CuBbwEDCLUjhH/zDgLW6dl0CX5nBLNESgrxqT5q+Rh3Dc3urhk7YYSX9vS2YUCerfKW2VpMCqp6F3CXiFyrqvd0UkytsnEKpjMFwmsyi5eAZOPtYPWRqzmLgq0fI6jLqStav9YxtU0nJVbrkdQluR3RHBSRHuEXItJTRL6emJBaZ20KpjNppPrIS8CTg7czqo9afCDHb2n1jrZvNFZtVUtpz+1v1+WqWh5+oaoH6NjcR8akjXBDs3i8+L15ZGtNh67npvrFTe+jRBv6xOkx9sROJJvXLEpMMKbTuE0KXolqyXJWTmt+fL4xGUajq4+8+eQEOzYhnptP5y0ljra2KU/d+3TMfbU7P2jbxVpxaPt6Fj15JwBbP1zO+689Etfrm8RzmxReBZ4QkdNF5HRgnrPNmIwXCDolBW8Wwaw8cjtYUnCljY25QX8d6xb/q81v054pLCp2b2XBH68jGGhamjn+/ZuZuubn+OpqGfzXk5m04Jttvr5JLrdJ4Ubg38DVztfrwA2JCqo1NnjNdKaAv75NQbPyyaWjJQU31UetlyaiB68F/3kLY146t8XBZfHS89VrmL5tLuvf+3fC38t0PrcrrwVFZC7whqrGt7zZDqr6AvBCaWmptWuYhAtGuqRmEcgqIE87OnjNRSnA5Syp4ZHIg+s+AoGqAzs6FJsbOTj330Lism6r6cvtiOazgWU4VUYiMlFEXC+KY0w6C/hDD0HJyobsfHLFV196aAd34xRaaGiOalQIDzrT8LZOfBgf2rmB8ltsDGumcVt99BNgClAOoKrLgGGJCcmY1BJwSgoebzaSE1qSs6b6cELfU12MU0i249+/mR40/3OId1dX03ncJgWfqjauwE/931pj4iBcUvB4c5BsJylUtT8puCopuFw8J1x9JFZdY+LE7cprq0XkS4S6ph4FfBN4J3FhGZM6oquPJKcAgNrqynZfz1VSaGGcgjif5TppWrwm4rWGg0lNbksK1wJHA7WEuqMeJDQnUlJY7yPTmYKR6qMcvE5S8CW6+qjFqYyS+1AuCdUit8gamtOX2+U4q1T1B6p6PDAVuEO1Mzprx4zHprkwnSYYqT7KwptbCEBdTftLCq4ag9tYJx9uaLaHsekot72PHheR7iJSCKwE1ojI9xIbmjGpITz3kTcrm6xcp6TQoaTgZkSzi6SQygkglWMzLXJbfTROVQ8CnwdeIdTz6KuJCsqYVBIMOCWF7Byy8kIlBV9NoquPWprmomn1UbihudMW4DEZy21SyBaRbEJJ4XlV9WG9j0wXEXTGJHi82eQWFAHg70CbgrvBay5GNDfzJ5gq1Udr/vdMskMw7eQ2KTwAbAIKgTdFZAihxmZjMl64odmbnUN+US8AAtXl7b5e9IM74Pez7M4ZrF30WsP3bOM4BU2xEsLkBdckOwTTTm4bmu9W1YGq+hkN/UZvAU5LbGjGpAZ1qo+8WdkUdA8lhWB1fHq+7du1lYlVCyh55cpGb+qmobmZxJEiJYVohw8eSHYIpg3atVqHhrR/nH8HWZdU05k04JQUsnLpVtSDoArUtL+gHN2IHKsNIFWqgeKh5rBVKqST+C3h1ImsS6rpTJGkkJ2Nx+ulkjyo7UhSqH/gS3hEcpNP/W6W42xuY2pVI4Et65lu0jIpGNOZwkkhKzu0rlSlFOKta/+ykw1KAZ4OlBSaOyaDShgmOdyOUzhfRIqc738oIk+LyOTEhmZMaogkhaxQUqj2FJLli+9axI1LCupyWmqbcsLEm9uSwo9U9ZCInAScATwE/CFxYRmTQgL1vY8AaryF5Pg7Mk4hOgG0v/qo+fNSTya1j3QFbpNCeHauzwIPqupL2BrNposIlxSynaRQl9WNnEAHRjRHdTeN2dDsZuW1NEgIJv24TQrbROQB4ALgZRHJbcO5xqQ3Z43mcJuCL6uI/GAHkkIzGqcGt+sRaHjGVPs0buLE7YP9i8BrwFmqWg70AmzuI9MlSLj6yBuaaT6QU0ShdmDqbJorKTR6qLf0kG+mTSF6gU5jOsLtegr9gZdUtVZETgWOBR5NVFDGpBIN+vCpl2xP6DNUML8X3fUQwUAAj9fb9utFd0mV8NoIbUgK9QdBrPONaSe3JYW/AwERGQk8CAwGHk9YVMakEAkG8FP/8JfCEryiVOzfnbD3DLqc+0ijvjcmHtwmhaAzgvkc4B5V/R6h0kNa2rerjMXP3ZfsMEy6CPoaJIWs7kcAcHDfjnZdLrq9IPx9+0oKoTNT3aF2/pxMcrheo1lELgQuAl50tmUnJqTWdXSaiw+fvpXj37+Z9//xWJwjM5lIgn78Ul/TmuckhcP7d7bres110Wza0Bw7KShRSSVGl9bFz97brtgSYcTTn052CKYN3CaFS4DpwM9V9WMRGQb8JXFhtayj01xMnPNr1meNYvTb32HjqkVxjs5kGgn6CESVFAp69gOgtmJX+y4Y3VAc4+GvWr9Gc9mGVbF/T6X5pHD8su+3LzbT5bmdJXUNcD2wUkTGA2WqekdCI0ugvIJu9Pza3zgsheQ/9WX2796W7JBMCpOgH39Un4zuJaGaU9/BeLYpxK4+GvTYiQx/6swYZzXsfWRMR7md5uJU4EPgXuA+YL2InJy4sBKvz4ChlJ89l15azs4/XUBdbdKWnDYpToI+AlJfUiju1ReAYOXedl2vYZuCs2Ja4xKDizaFUOmgUUnBxiuYDnJbffRr4ExVPUVVTwbOAn6buLA6x6jJp7Cy9OeMq1vJ+w9c7moUqel6JOgnENWmkJ2TywGK8FS2r6TQ8Lkdq/qohd9Fe/CbBHK9HKeqfhB+oarrSWJDczyVzrySBQMuYur+53n3b3cmOxyTgjxBH4FGv+77vX3IrWpfQ3Oz79Ge3keqkRXXrEuqiRe3SWGpiPxJRE51vv4ILElkYJ1p6qW/Y1nBdI5bcwer3no+2eGYFOMJ1uH3NEwKh3L70r22nUkh+oHvzIPUlllSGxwXc0I9Y9rHbVK4ClgDfNP5WgNcnaigOpvH62XElY9T5h3EoH9dzbaNq5MdkkkhXvXhl4ZJobagPyXBPe26XvNdUhslBVeL7DTTpmBMB7WaFETECyxX1d+o6jnO129VtbYT4us0RcW9yP7KfAB8j13AoYr9SY7IpApvsI6Ap+GkwNp9EN2p7PDvSfjh367qI4hUH9WzdjHTMa0mBQ11mP5ARI7shHiSauDwoyk74w8MCmzjo/svJOBP2jLUJoVkqY9Ao5JCVu/Qn8O+7RvbccVmeh81Lim4bmi2zqgmvtxWH/UEVovI6yLyfPgrkYEly/iTzmbpuJuYWL2Qdx/+drLDMSkgK+gj2Kik0K1kCAAVOz9u8/U02FwpoO3VP0KQcFLwilUfmfhwO0vqjxIaRYqZesGNLLpnDdO3P8qS54+m9Oyrkh2SSaIsrSPobZgUSgaPBqBm14YOXTtcImhP9ZFHg9aSYOKuxZKCiIwUkRNV9b/RX4RWYivrnBCTY/KVD7I65xiOWfpD1r/332SHY5IoS5uWFHr3G8whzYe969t8veYakZv2Pgo0OaYxD8HI1NnGxEtrv1G/Aw42s73C2ZexsnNyGXD539jn6UnP5+ewZ/umZIdkkiQbH9qopCAeD9uzj6TwUDvaFJoZ0dyW6p9wIcKjgUiX1Mb7jGmv1pJCX1Vd2Xijs21oQiJKIT379Kf2vMco0GoOPHw+NVUdWazdpKscfKin6ZLkBwuH0bd2c5uv16BLaswJ8VxUH9F6acKYtmotKfRoYV9+HONARIaLyEMi8lQ8r9tRw46eyvoTf81I34es/sNXCAbsD7GryVYfmpXXZLu/9yj6cICD5fvi/p7iomupR4ORWVIj51lnJNNBrSWFJSJyeeONInIZsLS1i4vIwyKyW0RWNdo+Q0Q+EJENInITgKpuVNVL2xJ8Z5l05ld496jrOO7Qv1n08HeTHY7pZDnNVB8B5A88GoCydYvbdL2GBYXmH/5uRjT30nKrPjJx11pS+BZwiYj8R0R+7Xz9F7gUuM7F9ecCM6I3OIPh7gU+DYwDLhSRcW0NvLNN/dJPeLfXTKZv+zOLn/19ssMxncTvqwvV92flNtk36OgTADj4URvX5HCx1GbLD/fQ+bniw8YpmHhrMSmo6i5VPQH4KbDJ+fqpqk5X1VYnflHVN4HGQz6nABuckkEdMB+Y1Y7YO5V4PEy66iFW5k5iwvs/ZvXbLyU7JNMJ6mqrAZBmSgol/Y5kJ33I3rms3ddv/PAPavgh73JkcuP6IisqmA5yu8jOv1X1HufrjQ6+50Bga9TrMmCgiPQWkfuBSSJyc6yTReQKEVkiIkv27Gnf3DPtlZ2Ty5FXPcUO7wAG/fNytn64vFPf33Q+X3idjWZKCgDbu42l3+E2zpXl4sHd4ojm6OMalxQsKZgOSplOzqq6T1WvUtURqnpbC8c9qKqlqlrap0+fzgwRgOKeJWR/9SkCeJHHL+DAHluUPJOFk4LESAp1/Y5joO5qf5flRg9xjbE9tkZtCkEfi+6Z075YjCE5SWEbMDjq9SBnW9oYMGwMuz/7Z/oE97LjwfOoralKdkgmQerqnKSQ3XxS6DMh1GS26d0X23X9xl1PI5/83SaFRgWFmp3rmbrv2XbFYgwkJyksBo4SkWEikgPMBto0j5KIzBSRBysqKhISoBtjjj+DVVNuZ5xvFSvvu8hWbctQ/rpQm4InRklh2Ljj2UcxsvHfrq/Z0rTY4aTQ4ojmqIShjf6EmyzraUwbJTQpiMg8YAEwWkTKRORSVfUD1wCvAWuBJ1W1TZWyqvqCql5RXFwc/6Db4LjPXsaCIVdRevCfLJx7U1JjMYnhrwvNEO+JUVLweL183P14hh98t52z6jb8MOEntBa0+uvcnS6Nu6TaOBrTMQlNCqp6oar2V9VsVR2kqg85219W1VFO+8HPExlDok2bcxuLi89i+pYHePeZe5IdjomzSEkhu+ngtTAZ8zl6cZC1C9z1SGtptLLPmaJb/e6WKwk2WhHOTXdXY1qSMg3NbZEK1UeRWDweJnz9UVbmTmbysh+z/I0nkx2SiSO/06bgjVFSADj61POp1Dyq3nuizddvPI12HeGkUOPufG/DuKSmvM0xGBPN7dTZKUVVXwBeKC0tbTLaOhlycvMY9o2n2XTX6Rz132v4oHsJo0s/meywTBwEnOojb07sWV3yCrqxuMcpjDnwb6orD5FfWNTyRaMnxGvUvuB3kgK+2CWF6JJG0NuwBDNtd9f8UBLw+yn7cDk7V75OsGIbOQe30K16GwWBQwhBApJNtbeI2qwiavNKCPY+irz+YykZMp5+Q0aRld10HEpXlZZJIRV1696THpc/x4EHTueIFy9ia/FLDD5qQrLDMh0U9Ieqj7JyYlcfARRO/xrdX32NRS8/wNTzr3d9/cYNyj7JBgUNuFzt1tO1/oSrKw+x/aOVHNiyCt+uD8g5sIFeVZsYENjGEPExxDlum/TlQE5/ducPRMWDJ1hHju8ghXV7GVy9ll4HXoINwP+gTrN4r+gERlzyJ4p7903m7aWEtPyNEpGZwMyRI0cmO5QGSvoNZuuX/44+9hm8j5/H3sv/RcmAIa2faFJWIFJ91HJSGDvlTD7811H0X/swwcC38Xi9MY+NblKorWw4M73fSQq0VH0UdQFxmzySbNctw9mXM4DKgkH4C/sh+T3wFvTEW9iLnILukQZzDfipqzyAv7KcYHU5VO8nu3InhdU76OnbxRG6jxHONOMBFXZ4+rIvbwi7uk/Hc8QYug8+mpGTTmFgdg4DW4inYt8udmxcycGytQR3rGDyrr+z897TOPyVvzNw+NjE/0BSWFomhVSrPoo2eOQxrD/7MQY9dz47HppFzrWv071H72SHZdop6DQ05+Z3a/E48XiomHQlpYuvZ8krD1H6uStcXb+2srzBa7/T0CzNNDRrMIh4GjYDenzpMZ371u6TKKjeyZDyRfQ+cMD1+hF16mWf9OJATj+2Fk/m4+Jh5PQbTa8h4+k/fDyD8gsZ1I54inv3DZUKjj8DgDULz2Xgq1/D/+hZLD/5V0z45BfbcdXMkJZJIdWNmnwKKyruY+x/rmD9fbMYft0rrdczm5QUqA0NTMwtaDkpAEw66xI+eu8+Bi65k5rTv0xefmGr59RVNuwsEfDkQAAk2LRLajAYxNsoKeTUlbf6Hqlg/NWPRn4eGgxy6FA5h8v3UFm+l7qq+tKSeDzkFfWioHsvinqUkF9QRH+Ph/4Jjm/ctBls7vky+sRXmPDm5Sza8zFTL7gxwe+amiwpJMixp53HksP7mLzkRlbd8wVGf+t5cvMKkh2WaSP1hZJCTislBQBvVhbVp97CiNcvYsH8nzH9kjuav2ZU9Y+/qrzBvoCEqp08zVQfBQJ+vFlZDRqn833lTY5LRdlRvbfE46GouBdFxb1gyOgkRtXQkNETqfneIt6/5zymrv0FCx6pYPqcXyQ7rE5nXVITqHTmlSw59iccW7OYNXefh68uPep/TRSn+igv311CH/+JWSwtOo3jNv2Rj1YubPX44OHGkzqG6tazfIeaHFtTXdlkWw//XldxJZs3Kz0+f+blFzL+W8+ypPsZTP/4Xhb/bjZVh1P7ORNvaZkUUmVEsxtTzv02C0ffwKSqt1lxz+x2jno1yaK+KgIq5LTS+yjaiDn3c1CK8DxzRfMPlKguqXKwftqv0O9GqBSQ10wJ4OBeZ/LFqJLGEU1mpjcdlZ2Ty+TrnmThoEs57sCr7P7NiWxe2+qaYhkjLZNCupl24Q9YMPybHHfoDZb+/qu2pGcaEX8NNeQ2aeBtSY+Sfuw47bccGdjC2gfmtDgvVk7l9sj3+3ZtjcxvVxA42OTYwwdaXcLExInH62XaZb9h9el/pnuwgv7zz2Tx03clO6xOYUmhk0y/6FYWDL6MKeUvs+T3F1mJIU2Iv5oaiT2aOZZjTjmHd4d/I7SE619vabAvuk2he239g37/9o8ipYBuwVD1UaXWl1AObfvA+S6UZHbS+VPHt8e7E/4v2SG02zEnf4HgFW/xcfYIJi7/Kds/XpfskBLOkkInmnbJL1kw8BKmHHiR9++ejd/nctIzkzQefzW17UgKANO+eivvdTuZ4zfczco3n2v2mL7++vU4Dm1bhzgP/B4aKinszBoQ2R/c+m6Dc/fmttQTP3UcP+sbyQ6hQ0oGDKHnJU8QxMOOv9+Y8TMip2VSSJeG5sbE42H65b9jwdCrKT34T1bcdS51te7muDHJ4Q1U42tnUhCPh9FX/oWt3sEMfuPrbN8U/qQfKg3spQfdpDpyfLBsaWTq63A//or8+l74I/b9p0EJs6owelmS1NWWqrdUdcTAYbw39FKOO/yfjJ8ROS3/t9Kpobk50y++nYVHfZfJh99kze9mUV3ZtKeJSQ3eQC117UwKAIVFPfB+6XE8BKn6y4UNeqDtzm447Kr3geVIo7mQarsPi3xfQjnvvXh/pJ052GMYrdlLj3bHbhoKzYh8JlM2P8jaRa8lO5yEScukkAmmffnHLBr3Q46tWsSW357B/t1ptfhcl5EVqMHncd/zqDmDRx7DR9NvZ2TgI5Y++YtIm8LhwiMjx2z2DGKYfyO5wYbdTqWo4Vw8xy/7QeT7nCNGtPree7MHtHqMcUc8HsZ87QF2ePrR+5Wr2LerLNkhJYQlhSSa+sXvseyEexji+4iqP5xO2YZVyQ7JNJIVrMHv7VhSAJh01hyW5U/j6A8foM6Z7yjQs/6T/q4ek/CKMjS4lYDWL5zjzW06aK5m1wYACkparz46nJfoscBdS1FxL2rPmUt3PcTuP32Rmqr0mGakLSwpJNnks77Kps/Np1APU/jYp1m3+F/JDslEyQnWEohDUgAo+NTNFEk1lYsfC127ZHj9zsFTIslgr/Sq3+6pn1jvo3NDVRZ5ZW8B4M1qfbrnYFbHY9cfH+jwNTLJiGOmsXrqHYyuW8Pa35+fcR1GLCmkgDHHn8Hhr7xCpRQw/MUvsujJOzO+h0O6yNb4JYVRk0/lY88QxpT/BwBvbv3cSFn53fk4K5QkKrKiJ1CsLzUMP3oKu+nFqMr3XL+nBH0x922fs8jdNTKgoTjejvvMpSweeyOTqt7hvfsuzqi/17T8307X3kctGTzyGIqufYu1Bccxdc3PWXLXhRlZNE03OVrbZCGbjth1xEl0JzSfEiIcpND5VtjXaxIA1VnNd6AQj4ctRZMokNrIOY0FteG2YFbsxYEGDBvjOu51WW2bTnrlJ+e26fh0NHX2zSwYeDFTDrzE0lceSnY4cZOWSSHdex/FUtyrD8dc/yoLBl/O8RWvsuNXJ7JxlbtPcyYx8rWaYHb8JjLMHjy5wevd3n4AKB4oOQqAPH/T0cxhvl5HtXj9d/smZsrnkTe+Sa1mt36g45iTv5CQOFLN8Rf/kg+yRjNy8S3s3bkl2eHERVomhUzm8XqZfumvWH7KnygKljPob59h4V9+bCOgk0CDQQqpgZz4TXteMnxS5HsRqMmqv3Zen1D1UW//Lqo0t/6gKLn9oj/d1+8Lt0eMveBW9n99Tf0h6m7dgtZkZedQK+6TQleRlZ1D3vkPkq+1bH30yoyoRrKkkKImnHY+nq+/w6pu05n20V18eMcnXM26aeKnqvIgHlFopgdQe/UZ3HC1QH9Ud9ce/UNdTEso54CnR7PnF/ZpfiW/AE6DtHgo6lES2d543EPHNK2uMqEpt98/6homVb3D0hcfSHY4HWZJIYX1OmIgk777PIsn3UZfXxlDn5rBot9fQsX+xtMtm0SoPhRqs5K87nG7ZkG36CpPD0FvqESg/toGCeOQtxfNKerd/LgDP/W9lLJzogbbRc3IuvCIC9oUawWFvNvzs5HXQScpLBz+zTZdpys4fvYPWZc9jqPeuzXtxxxZUkhx4vFw/Kyv4/nmeyzpcw6le56Buyew4M83crB8X7LDy2hVh0NdMb15iVs1L9yzKVBXTUG3YvzqYVHJOdRlhUonIsL6rFEs6R5aNrJHn/pxB9ENzQHnTzkYaFTNGFV91GNKfVIIN3C3JPiNpUy57vH6SzlJYdSMq+CWCnzf303tTTsanPOxZ2ir181E3qws8s65hwKtYcPj1yc7nA6xpJAminv1Yeo1D7P5/NfYWDCB6ZvvR393DAse+i67t32c7PAyUq0zyCyrIL4dGiI9hEQi4wjUF5oDKeunB5h6zZ9Dy3I6Rv1wMaXf+TvQuKRRTyX0p6xOyWDxxNCKYaL1SUKch3qdeqm7qu1VkeuGzQnFUNQDCJVIGq8m2OPqV1kz44k2XzsTDB1bytL+s5lS/jIblr+d7HDazZJCmhk+fiqTbniFDV94iY0FE5i65SF6PTiZ9375OZa/Md8m2Iuj2spyALLz41d9BFBJVBdXZ/lNbfQJP+htfWAaUv/nG67aCU+hIc6gNwkGog4PbducPZySfvVTbMSi2rDRdPqcn8MtFS2uPd2zT3/GTZvReuwZatzsW6mgkMpXb0l2KO2WlkkhE8cptNXICScx6YZX2DFnAYsHfIlhlcuZ8OaV1Nw2nMW/m817r/yZ8r22KEtH+JwF5XML41tSqJb6sQPqcZap1IYLLwWdkkKT6qAYNJIUQg9y8YauK1HXDVc3iTbtIfPuMbe4ep+W2OR70L1Hb9YMv5QJ1e+m7aR5aZkUMnWcQnsMHD6W6VfeS+HNH7L85Af4oPgTjC7/L5MXfYvu94xhw62TWXjf5Sx+7j42r11qXVvbwF8dSgr5TnVJvNRKqKQgSKSkQLBxSaG+ATqW6DaFcFIg6JQUmkkKjbu3Rjv205c32ZZf2LYSUuCyN9p0fKaacM717KUH+q+fpWUX1fRYTdu0Kic3jwmfnA2fnI2vrpZ1y97kwOp/UbzjbSbsepb83U/C+1CluezIGkBF/mBquw8jq2QEeb0HU9h7ID2OGESP3v3weL2tv2EXEKwJJYW8OJcUaj35EH5Wh9sCGj08NJwUfO6qAzcXHEOvqrfJdur4xdNCSSGqm+reK1YQCPrpldV0DEJ+Ydsa2PsOan3W1q6goFsxK8dezdS1t7HizWc49tRzkx1Sm1hSyEDZObmMmfIpmPIpAPy+Oj7+cDl7179LYPsy8g9tok/lBvodepvs7Q2rLXzqZY/0oNJTRE1WEXVZ3fBndSOQU0QwtzuS2x2y8/Fk5+HJyceTk483pwBvTh5ZuQVk5xaQnVdAdm4+OXmF5OYXkpdfSFa2izryFBOsCa1zUdi9Z1yv6/fkhpKCSKT6KLruHwCnTaGlkkK0UVc/zoaNqxnZK7REZ7ik4IluaI60QdQnhZIBzY97MB0z6fPfYvu6hyj438/Rk7/QrvmjKg7sZf2//0qgch99j5vJsHHHJyDSpiwpdAFZ2TkMG3d8k18qv6+ObVs/4uCeLVTt246vYjt6aBdZVbvJqqsgx3+Iotqd5FdXUqiVdNOqyIpgbeVTLzXkUCu51EkOPsnFJzn4PLn4PbkEvHkEPLkEvXkEs/II5vVA8nviKexNTrfe5HbvTUFxH3r1H0pRcfN9+OOuppw6zSI3N/b8Qe0RkKiSmNP4q43aFDTc0Nxi9VH99/kFRYyccGL9ZZspKYRLJZ44jXJuzbqssYzxr+2U90o1Obl5bJtwHccv+wHv/eNRJs+42PW5vrpa3nvmd4xYex/HUw5AYMPdLBp7I1Nn35yYgKNYUujCsrJzGDh8LAOHu5vsTINBqqoOUVtdRW1NJb7aanw1VfhqKvHXVhOoq8ZfV02wropgXTXqC3/VgL8a8dcg/ho8gRq8gRq8gVq8wVqyAzUU+CvI1lpygnXkUks3rSRbAs3GcZBC9niP4GBuf2q7DcLbbzw9hk1i8OjJ5BXEb/Sxt3o/5dKdI+I8S2hQov7swp/eGzf+ekPVORJ01wbUeHI8cc73NFN9RFxHOcc29Luvc+DwQeJbzkofkz93FZtX3E/vd39J4Iyv4M1y97hd8pfvM33rn1ifNYq9Mx6mz+DRbH3kMkrX3sGK/4xKeHWUJQXjmng8FHQrjtlXPp40GOTQoXIO7d9NZcUeair2UHtoL/4DZUjFVnIrt9OjZjt9K5dSsPtJWBGa/+fDrBHsPWIa3caczlFTzmqx+2RrsuvKOeQt5og43hfUJwXVYOykEC5BtJQUJHayEm9oX4M2BSe5xXfqi9jynKrDrsqblcW+KdczeeF1LH7xfo7//DWtnvP+Px5jypaHeL/wRCbd8HJke97V89n821MY+p9r2TpoFINHHpOwuC0pmJQkHg9Fxb2cqqLYUzwHAwHKPl7L7g1LqC1bTvGudzlu+zxydjzGwX8XsKLnaXSbehFjp5zZ5nrdPF95zGmsOyLoCX2KV7+v/sHeuE3BOaZxV1UIDT7LaVSKanxvnkhJoT7ZNNfQbBJr0pkX8eHiexi07HdUnfHVFj9Qfbjsf4x9+1tsyB7FqKsfb7CvsKgHeV99gsAjZxB8fDaV336bwjj3igtLyy6pxoR5vF4GjRzP5BkXM/2y3zLuB2/j/97HLD/5j3xQfDJH73+dca9ewAe/OIHlb8xvUxfBwkAFtdk94h6zOiWFYKCOsTO/zYq8Ukad3XBqhPrBZ01LCuHJ75pbTyHME2loDvChdySLi8+KjGhurfpoF71b3G/cE48H3xn/R3/2sOIvN8Q8bt+uMro/O4dy6UGfK55p9oE/YNgYtp1xH0OCZayYf0vCYrakYDJOQbdiJnzyixz/7Sfw3LCBReN+QA//Xia8eSUr7/wU2za6a/wsCh7Enxf/Ru1wj6NgwEePkn4ce9Pr9O47qOFBkV5JsZNCSyIlBQIc9aOlHP/tJyMt062VFPK+uYhtF9mMvPEybvqnWdT780zZ+USzy+3WVB1m7x/Ppbse4vAX5tLriIExrzX+pLNZUnQ6k8oeY+fWDQmJ15KCyWj5hUVM/eIN9P7+ahaOvpHh1avp9cjJLHnxwRbP8/vqKNbDaF78m0mDzgNfW1jbN1xSaNLWAAQk3DbgrqRQf81wO0PLSaG4Vx/XnQ+MO0fP+R27pYTeL13W4ENJdeUh1tw7m6N8H7D2hN8wcsJJrV5r4Hm3EcBL2fuJWc89LdsURGQmMHPkyJGtHmsMhMZuTLvw++wqu4B9j1xE6ZLvsWDbSqZdflezbQ17d26hnyhSHPtTW3uFq480EHv95JgN0IDfxZ9ts0nB2hSSplv3nuy9YD49npiFPDqDRb1PRUtGc+T6PzNZd7Nw9PVMO+urrq7Vf8hoKq5bRakzJiXe0rKkYNNcmPbqO2gEI69/nUW9ZzF9x6O8+4fLm21nKN8Rmnk2r6T1iePaSp2qnZaTQuz6fzfVR95wUqDpOIXO6pJqGho6tpQ9n5/P9rwRjN/7KtM+uAOf5LD6zHlM+9KP2nSt4gQlBEjTkoIxHZGTm8eUb8xl4QNXM23XfBY83J3pl/22wTGH92wCoLjfsPgHEGO+o4bHhKt6mqk+iqyy1sJKaE71U3TvI08rXVIXHf0jiodNaqGvV/NW5xxDwJPDsW08rys6auInYOIb+H11bPl4Df2HjCEnN6/1EzuRJQXTJYnHw9Qr/8C7dx9ketnDLH15PMd95tLIft++0CLsvQcMT8S7A/XTXDdn5Mmz2b7mAfqe+Z0m+wKS1eqH/ebzRctJYer57Vsc5ujvv9Wu87qyrOwcjhw1MdlhNCstq4+MiQfxeJhw1UOszR7HuEU3NVgDO3vvWnbTKyFTamgL7QVhvfsOYsBPPmTImMlN9kWmyVBl/dnPs3DQpU2Oaa79oM/AYWyXI9h90q0diN5kOispmC4tN6+APpc+yaH7Tybv6YuoGPg2xb360Pvwh+zIGxH30cyhNw1NSS3Z7as20PDSm0E/oyafApNPaXpQM+0HObl5DPjJhzS/yrMxIVZSMF1eSb/B7P/cn+gT3MumBy9k97aPGRbcRFW/0oS838Qv/YwFw77B5FnXtuv8cEkhGGh+biio765qPY1MW1lSMAYYU3o67x99MxNqFsMfPwnAgKnnJ+S98gq6MX3OL8jOyW3X+ZGSQgursoknnBSMaRtLCsY4ppz3XRYMvpxsfCwYeAlDxh6X7JCaFW6TaLyGckPW/dS0j7UpGOMQj4fpl/4K+BXTkx1MCyLrMbcwj1P9esyWFEzbWEnBmDQTrj5qqaRg1UemvSwpGJNmarNCaxRIS+spROZHSr+F401yWVIwJs0MuPgRFg79BqMmnxrzmJL+Q1ifNYrdn/x15wVmMoK0NKoy1ZWWluqSJUuSHYYxxqQVEVmqqs32ubaSgjHGmIiU6X0kIoXAfUAd8B9V/WuSQzLGmC4noSUFEXlYRHaLyKpG22eIyAciskFEbnI2nwM8paqXA2cnMi5jjDHNS3T10VxgRvQGEfEC9wKfBsYBF4rIOGAQsNU5LPb4fWOMMQmT0KSgqm8C+xttngJsUNWNqloHzAdmAWWEEkOLcYnIFSKyRESW7NmzJxFhG2NMl5WMhuaB1JcIIJQMBgJPA+eKyB+AF2KdrKoPqmqpqpb26ZO41YeMMaYrSpmGZlWtBC5JdhzGGNOVJaOksA0YHPV6kLPNNRGZKSIPVlRUxDUwY4zp6hI+eE1EhgIvqup453UWsB44nVAyWAx8SVVXt+Pae4DN8Ys27kqAvckOooMy4R4gM+4jE+4BMuM+0v0ehqhqs/XvCa0+EpF5wKlAiYiUAT9R1YdE5BrgNcALPNyehAAQ66ZShYgsiTVqMF1kwj1AZtxHJtwDZMZ9ZMI9xJLQpKCqF8bY/jLwciLf2xhjTNvZNBfGGGMiLCkk1oPJDiAOMuEeIDPuIxPuATLjPjLhHpqV1rOkGmOMiS8rKRhjjImwpGCMMSbCkoIxxpgISwpJJCKFzuR+n0t2LO0hIp8XkT+KyBMicmay42kL52f/iBP/l5MdT3uk88+/sQz4W/CIyM9F5B4RmZPseDrCkkI7tHGdiJbcCDyZmChbFo97UNVnnfUvrgIuSGS8bmTC+h1tuYdU+/lHa8fvV9L+FmJp4z3MIjRlj4/QJJ9py5JC+8zF5ToRInKMiLzY6OsIEfkUsAbY3dnBO+bSwXuIOvWHznnJNpf0X79jLu7vISxVfv7R5uL+9yvZfwuxzMX9/8Vo4B1V/Q5wdSfHGVcpM0tqOlHVN505naJF1okAEJH5wCxVvQ1oUiQWkVOBQkK/WNUi8rKqBhMZd7Q43YMAtwOvqOp7CQ65VW25J+rX71hGCn04ass9iMhaUujnH62N/xfdSOLfQixtvIethJYShtT6kNFmlhTip7l1IqbGOlhVfwAgIhcDe1Phj4A23gNwLXAGUCwiI1X1/kQG106x7ulu4Pci8llaWL8jRcS6h3T4+Udr9j5U9RpIub+FWGL9X9wF3CMinwDeTEZg8WJJIclUdW6yY2gvVb2b0MM17WTC+h3p/PNvTpr/LVQBlyY7jnhImWJzBujwOhEpIBPuobFMuKdMuAfIjPvIhHtokSWF+FkMHCUiw0QkB5gNPJ/kmNoqE+6hsUy4p0y4B8iM+8iEe2iRJYV2kNA6EQuA0SJSJiKXqqofCK8TsRZ4sr3rRHSGTLiHxjLhnjLhHiAz7iMT7qE9bEI8Y4wxEVZSMMYYE2FJwRhjTIQlBWOMMRGWFIwxxkRYUjDGGBNhScEYY0yEJQXTJYhIQESWRX25mdo84aLiGtDCMT8RkdsabZvoTIiHiPxbRA6LSGmi4zWZz8YpmC5BRA6rarc4XzPLGczUkWu0GpeIjAJeVdXhUdtuB6pU9WfO6/8A16vqko7EY4yVFEyXJiKbROSnIvKeiKwUkTHO9kJnkZV3ReR9EZnlbL9YRJ4XkTeA10WkQESeFJE1IvKMiCwSkVIR+ZqI/C7qfS4Xkd+6iOdMEVngxPM3EemmquuBAyISPWPtF4F5cf1hGIMlBdN15DeqPopeqWyvqk4G/gBc72z7AfCGqk4BTgN+KSKFzr7JwHmqegrwdeCAqo4DfgQc5xzzJDBTRLKd15cAD7cUoIiUEFow5wwnniXAd5zd8wjNs4OITAP2q+qHbf8xGNMymzrbdBXVqjoxxr6nnX+XElqmE+BM4GwRCSeJPOBI5/t/qup+5/uTCM2lj6quEpEVzveHndLE55y6/2xVXdlKjNMILTTzdmj9InIIzb0D8ATwjoh8l1BysFKCSQhLCsZArfNvgPq/CQHOVdUPog90qnAqXV73T8D3gXXAn10cL4QSzoWNd6jqVhH5GDgFOBeY7jIGY9rEqo+Mad5rwLXOkqOIyKQYx71NqH4fCa3Ve0x4h6ouIjT3/pdw98l+IXCiiIx0rlfoNDKHzQN+C2xU1bReHN6kLksKpqto3KZweyvH3wpkAytEZLXzujn3AX1EZA3wf8BqoCJq/5PA26p6oLUAVXUPcDEwz6mGWgCMiTrkb8DRWNWRSSDrkmpMB4iIl1B7QY2IjAD+BYxW1Tpn/4vAb1X19Rjnx6WrrHVJNfFiJQVjOqYAeEtElgPPAF9X1ToR6SEi6wk1cDebEBwHWxu81hoR+TcwHPC19xrGhFlJwRhjTISVFIwxxkRYUjDGGBNhScEYY0yEJQVjjDERlhSMMcZEWFIwxhgT8f9iTynDbQcT1QAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# Plot cross sections at 600 K\n",
"plt.loglog(broadr_xs[1].x, broadr_xs[1].y, label='BROADR')\n",
"plt.loglog(acer_xs[1].x, acer_xs[1].y, label='ACER')\n",
"plt.title('600K cross section')\n",
"plt.xlabel('Energy [eV]')\n",
"plt.ylabel('Cross section [b]')\n",
"plt.legend()"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "9a9058eb",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.9.7"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment