HajimeKawahara/exojax

View on GitHub
documents/tutorials/Reducing_memory_for_HITEMP.ipynb

Summary

Maintainability
Test Coverage
{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Reducing memory for HITEMP CH4"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from exojax.spec import moldb\n",
    "import numpy as np\n",
    "import seaborn as sns\n",
    "import matplotlib.pyplot as plt\n",
    "plt.style.use('bmh')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "First of all, set a wavenumber bin in the unit of wavenumber (cm-1).\n",
    "Here we set the wavenumber range as $1000 \\le \\nu \\le 10000$ (1/cm) with the resolution of 0.01 (1/cm). \n",
    "\n",
    "We call moldb instance with the path of par file.\n",
    "If the par file does not exist, moldb will try to download it from HITRAN website."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Setting wavenumber bins and loading HITEMP database\n",
    "wav=np.linspace(16370.0,16390.0,2000,dtype=np.float64) #AA\n",
    "nus=1.e8/wav[::-1] #cm-1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    " when initializing MdbHit, extract=True option does extract  the opacity in the wavenumber range we now considered. It reduces the use of DRAM. (If you have very large DRAM, not necessarily)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Extract HITEMP: 31880412it [01:47, 296364.94it/s]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "self.path changed: ~/exojax/data/CH4/6101.281269066504_6108.7354917532075_1.0/06_HITEMP2020.par\n",
      "HAPI initializes all the par files in  ~/exojax/data/CH4/6101.281269066504_6108.7354917532075_1.0\n",
      "HAPI detected: 06_HITEMP2020.header\n",
      "HAPI detected: 06_HITEMP2020.par\n"
     ]
    }
   ],
   "source": [
    "mdbCH4=moldb.MdbHit('~/exojax/data/CH4/06_HITEMP2020.par.bz2',nus,extract=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAD1CAYAAACx81UXAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABoS0lEQVR4nO29e3Rb130m+m08CYIvEHzKpCjJIiXTluXYsS0rmiaxrUSxkvrGqT3xmsm0d9w7y42cZtpO2pvJzF2Z1cnqZDzJXCdeiddd7fRO09a5dpxHG9mK5KbJVC9bkaKXKYmUJVGkRIAkABIgCIIEsO8fB/tgn4NzDs4BDohD4nxrORFxgH32bz++/du/x96EUgobNmzYsLG+4ah1BWzYsGHDRvVhk70NGzZs1AFssrdhw4aNOoBN9jZs2LBRB7DJ3oYNGzbqADbZ27Bhw0YdwFXrCvziF7+gXq+31tWwYcOGjTWFxcXF2ccee6xT7/drTvZerxfbt28v+nx8fBwDAwM1qJE1UO/yA3Yb2PLXt/yAdhucOXNm3EhZljXjuN3uWlehpqh3+QG7DWz561t+wNw2sCzZt7a21roKNUW9yw/YbWDLX9/yA+a2gWXJfnZ2ttZVqCnqXX7AbgNb/vqWHzC3DSxL9vW+qte7/IDdBrb89S0/UCea/fLycq2rUFPUu/yA3Qa2/PUtP2BuG1iW7FOpVK2rUFPUu/yA3Qa2/PUtP2BuG1iW7Ht6empdhZqi3uUH7Dbo6enBSDiJV8+GMBJOFj3XerYeUO/9D5jbBpYl+1AoVOsq1BT1Lj9gt8GxyxP4kzfH8D9PT+FP3hyTkPpIOKn6bL2g3vsfMLcNKiJ7QshXCSG3CCFn8/89wT37MiHkKiHkCiHk40bL9ng8lVRtzaPe5QfWVxscvDSLL791FQcv6Y+ueCecwXKWIkeBlRzFuakEAIHov3tyEmmFZ+sJ66n/y4WZbWBGBu1/p5T+N/4DQsgwgM8CuBvABgBvE0KGKKVZvYU2NzebULW1i3qXH1g/bXDw0ixeOjYBADh9SyDl/Xd1aP5mJJzEP02mwO6RcxKCnb3NGAkn8aWDo1jJFb7Lnq03rJf+rwRmtkG1zDhPAvg+pTRNKb0O4CqAh4wUEIlEqlKxtYJ6lx9YP21w9Mac5t9KODeVQDZP6ATAx4baMdztx7mpBDI56XfZs/WG9dL/lcDMNjCD7F8ghJwnhPwPQkgg/9kdACa470zmP9ONQCBQ+kvrGPUuP7B+2mDPpjbNv5XQ4nWJWj0FsDXYCADY2dsMFzdr3Q6CvYNBU+ppNayX/q8EZrZBSTMOIeRtAEou4a8A+C6AP4UwHv8UwDcA/GsjFZiensZzzz0Hl8uFbDaLp556CgcOHMDk5CS6urrgdDoRj8fR2dmJaDQKSik6OzsRDofR1NQEAFhYWEB3dzdmZmZACEF7eztmZmbQ0tKCbDaLZDKJnp4ehEIhuN1utLa2YnZ2Fq2trVheXkYqlRKfezweNDc3IxKJIBAIIJVKYWlpSXze0NAAn8+HWCyGYDCIRCKB5eVl8bnP54PH48H8/Dw6OjowPz+PlZUV8bnf79clUzgcxpYtW9aVTEb7KZPJIBaLrXmZHt/Sg2i0Eb+eyWDP5jbc05jE3JxLU6aZODcHAdwMzyLcnkWfz4c/eqAZZ6IEKysreLjbhTvbXBgfH69ZP1Vr7E1OTqKjo2NdyWR07DkcDsTjcUWZjIJQSkt/S09BhGwC8FNK6T2EkC8DAKX0z/LPfgbgq5TSE/LfnThxgtqnXhaj3uUH6rsNmG0+SwXt/etPDK5LU40W6rn/GUqcenn6scce+6DesiqNxunl/vw0gIv5f/8dgM8SQryEkM0ABgG8a6Tseo+xrXf5gfpug+FuP/7Ngz34wIZmPL+rr+6IHqjv/mewUpz9fyWEXCCEnAfwUQB/AACU0vcAvAZgBMAhAAeMROIAqxdja9XEFDvGuL7bYCScxP9zKoRf307glZOTlhufq4F67n8GM9ugotBLSunnNJ59DcDXyi27oaGh3J/qBktMWclRuB0hS22VV0N+q6Oe24BF4+RQiKMf7vZjJJzEuakEdvY212Ssrub767n/Gcxsg5rfVKUGn89X9Xecm0pgJVecmFLLycSwGvJbHfXcBjt7m+FyhpDJUbgdhRj7Wionq/3+9dz/ehdNM9vAsmQfi8XQ0tJS8nuVaBo7e5vhdoTyg5egxeuyjKavV/71jHpug+FuP/7w/iaEaZM4tl89G6qpcqKkHFXzneu1/40smma2gWXJPhgsHTtcqaYx3O3H158YFCeLmYO50u2uHvnXO+q9DR66s0cSZldr5UT+/mpn7a7X/jfCM2a2gWXJPpFIlIwnLZec5UTM/6acwSwvz4ztrh75y0Wt7b56Uc02WAuQy19N5UQP5O+v9thZr/1vZNE0sw0sS/Z6Du0vR9PQIuJyBrNSeUfGoljOUlCUPwmrdXFDre2+RlDPl1eMhJP45ZV5fNjRIukfM5STSiB/fzWxXvvfCM+Y2QaWJXs98aXlkHMpbcjoYJaXd2QsgsOjkaIDrIyiWjHGZmuD1dwl1GuctbggZynevD6muSA/PtgOgGDv4Po7H2c9979enrFSnH3VoDe+dLjbj2fv69E90IXdAIGDwBRtSF4eQJDNMz1/gJVRVCvGWK/8evIPqn2meqk2sGqORKUQF2SoH1/M2v6tKxG8PbY+Dwyz4+wtFGdfTVQr7Mpsu6O8PAB4eywibq3LPaSqlvLrNfVU22as1QZryRxlFKJ5Mqtunllte30tsJ5DL/WiLkIvq3lxAdtCMc2wUtKXb8nMWEyqIT9vcmEOPgBFddRLJNWOztBqg/VMdmxBfuf6DB7e3FmTtrcC7MtLrHd5SVUwPz+Ptra2qpVfTc3QDCeW2fLz8jrJFACCTI7CQabwwu5+yWUaeomk2tEZWm1gNbIz23cx3O2Hf2kWAyplrXZkTC1QbQ5YCzCzDSxL9h0d2jf5VAqra4Zmy8/LSylA8y7kLAVePj6Bze2+sqKSqhmdodUG1SY7I+RdLcWh1BhYzciYWqDaHLAWYGYbWJbs5+fn4fdbN0Gk2rHqZsvPy+skAsnn8o5kSlFxVFI1UKoNqlVHo+R9bipRCLXNmqc4VHsOWB3rSf5y+cLMNrAs2a+srFS1fF4zbPG6VO3XSlgN56DZ8ss14evRFF4+PgFKAbez9mYQJVR7DKjB6K6Pv1Uql//bDNRKfqtgvchfCV+Y2QaWJfvViLFlDW6kI0bCSXzvzFTFSVOlUA35WR0Z4X/jk0OWtvnWKs7a6K4vns6AQLiujeT/VoJR7W49x5nrwXqRvxKTsR1nbyL4jljOCklRamAr9JlbCVAIjVct52A15Ge3H/2/v5rClw6OAoChHIVqQS1e3ow4+3Ji8dku6Lcf6NWlhe3sbYbHKeQueFR2SeXkJNR7nLlV5Tc6pirJ7amLOPvVstXt7G2Gk0wJjksAh0ej2DsYFCc4r42xhYFpcB+4oxmfu7+3KmRZDfmPjEWwkhP+vZIT/rYC0avtrLTaQM/WuJLtsxF/gJKzWK7Fl6PdGRkD5diErX5GEpO/VD1XU45yxlQlwQRm8oBlyd7pdK7Ke4a7/fjYUAfevDwLCiBLpRdF8B37/K4+uB1E3N5Xi+iBaslPSvy9+tAiQa020EOeZkVc6SETfnE4eGkWLx+fQI4CHqdACC1eFwiM7Qb1jgElAgK0jz5mv1nOKoffWgFOp7Mkua52cp3cEvC9M1O6eKDcYAIzecCyZB+PxxEIBFblXXsH2yVZr2wiyskins6sWmyz2fILW04Kl4MgmxP+f+9gu2nllwst+7hWG7DfLWcpCJSdombE4uvdQTBH/9XIIt66EhEjnVayFK+dD+Gdm3FkKeAkUL1TVr6o8PJrLThyAnrtfAinJxOadeYjiFj4LSD4G6yi6cfjcZyLeTUX7NUOoebHHQVw5lYCF0MLVVtkzOQBy5J9Z2fnqr1LbZulRBarFZJopvzShCrgie0dljk4S2uLq9UGw91+PL+rT9Sgv3NiAlcjKYlceo+G0Hpeikx4DZkW/VqIzjl5My4Jc1Vy4CotKpvy8qstOAcvzeLojTnc2e6Dk0A0RZ68Gc/nUqgT4M7eZjjIlHiOUy5P+BQwRUM2w7TS2dmJnQ1Uc1Ff7eQ6Nqa+d2ZK9N1Vc5ExkwcsS/bRaBSNjY2r9j4lEtdjizUCI781U36esACgq8ltCaJnUFtAo9EobiSoavvH0xlQMFIDDl6exdtjkaJjq9Vk1WP+KEUmvB9HDTnuISEQrxgsZdNvWcmgsbFR8dn1aAovHRO08dO3EuhuciO8IITpUQo4CPLErUyAw91+vLC7Xwy/JdxiYcbFPWaYVqLRKIb7+iSL+isnJwFIdyCrnUk83O3H5+7vxcXQQtUXGTN5wLJkT6nW9ClGtZw0PFlUMoiN/tao/FoQnNAhUEoNHblcawfeWGwZ3/yV1GfynRMTyOQAl2MKn3+kX9RoGdSSmpQumPnuyUmk86rtSj4S6+2xaFEfyQ+6489Tkm/rHQCcDmBLsBGjM4uSRcBJgBd29wNAkb1caVGhy1EAyv33vTNTEvkY0RMIEUHP7+rD1cgitPwy++/qwOZ2n2iCeuXkZEVJhkIkm/A+M0wrbA7wi/pyloo7ECeZwseGhF3qs/f1GD7rih8TgLHrHVdrkTGTByxL9ka2L6vlpKnEPqj2WzVCNd+MRcX/vx5NlRykRtvUjB1Pi9cl0dgm0x5Jmx26Io0m+tF70+ht8eLmXFosi9ecGflsDfo4IissGqwsQDC3xFKZoj4CCiTAEtF4x6s8OY/VH4DEdMZIabhbuEtWbi//xieHishjaYl3zlHJ/+/Z1IbTt4qPPgYEn8Dm9oLM8t0OD16ZYcSv94gInii/dHBUbE8nESif5P9drtbL5gC/EBIUdiA5CryZ3809v6tP0sd6xqv8rKgsNcYfq2HSrQszTjgcxsDAgK7vrpaTphL7oNJvtQjViPylcG4qgWx+gmRy+myzRtrUjB0PIz9BMxXKuMO9JEY/OQkwtyTNJuRJHhBMFy/s7sf1aArfPjaBHPc5b8M+emMOmRyKEPC5JNFW/B2vzrx9m3e8sjZRm/Rqmp+SvfzcVKIo54GNAb7/stx3AeDQlQhGZ6U7iKM35nA1kjI8J/SSl7y/Hx8MStozmzcjCSg/4ovJL19QXzk5KY4Xvk+NjFc+MVJQnrWTJM3c5Ropy0wesCzZG7l3sRwSLqfzKtm6Kf32paMTqpm4avKXU281zYglkSmVY6RNy1lsmRzTCysSmzffFh/rD+DrPb3527+ioqlCC1PxNN64OA2ey3P5KBhmw96zqQ3npxISzZ7dPbB3MCjJq+AJgSdUitLHIqiR53C3H5+5pwuvXZiWlCXvWzYG1PqCmWH+6Kej4sLBIkRcjgScRMjrdRJgemEFI+FkWU5qOeT9DVA48uctAQLRs0Uxo+LU1vM++f27/A6EjQlBGxf6VI8NXa5cMLNbQbMv/q18F8Dv0ozCqGJk5h28liV7IzBKwquVbKP125FwEodHZ1WvLxyNpHFjIiSxJ0rtqsa2m7xmxEwYSklkSr8p1aZGF1t5dJCTEFBa0Oz5qAvmvMzwhnkV5Cjw+oXpImep20Hw+Uf6JCYiRhixVAYBn1sMQ+XlvR5NSRYhl4NI6vHKyUnJaaFyGbXazu91So5YuBpZLOrb7vzsZH3x2vkQIosZXI+mxDLPTSWK5GU7gE9sE2Q6PBrFW1eKndesnnpDS3k/Be9D2BpshIMIxOsgwEc2B/DzazGxLvyiaIbJlc0jfmFmfVpqvKolRrJnpUJbedNROXVf7VBRHpYl+4WFBQSD5d3yVAp6EyOq6aBkW3Og+PrCkXAS//loOB+XPQUKITaecOYILa1cqd78QnM1klJMIpND78JWyonJ4+ClWbx+ISwJVfzEtnZ0NXmQTGcFrTwfdeG7vwmPBoNFZg85GGkyMKctAbBroBXP3NstaVtWty/u2ShpMzkJ8WfeAMBD/c1IZ2jJkDs9hCYcsVA4hfRqJFW0y9sdSItz4Ho0hePjcQDAlZlFAIJmzy+0TiK0BtNQGRlmaTG5sLBNr8uhK7RUHrHE+xCuRlLIco7Eea7dHJCGmsrJ7shYRNVJurCwgHCmQfW5fHzqGa9yxYSf92q/lTvhjUQsKS2URhQjM3nQsmTf3d2t+7tGtQU9iRFKWZCVEH6pTuevLzw3lUA2JzgNc3l7IiAQPbOAUgj2WnbZNPudkvbPnhUGVj65SmXbqlcGHmyiKWUds6gQv9shmi4AQRYmO3Nc8lEXZ+eceBTFYYLybfeTw53iIuEgwGfu6YLf6yyqp9Y4UdK4dvY2w+WYEs09pyYS+PwjfSXNBUpl8X3A2urrTwziyFgEPxuNigTOZ9h2t7aJZR69MSd5x9Ebc9h/V4fiQivfncgzdw9emhXDNoGCQ1XJmaomC+9DAKjE16FlUpEvTswUo+QkHYt78Mrbo8hRIfqKf/78rr6yEsDKMcXyfcWbjozsYNWiu0q93wgPloJlyX5mZgb9/f26vqsn8YVFZjBbm1ZixEg4iZePT4iaZCVnlLN3FwZJodOf39WHozfmsGdTm/heRthOB8CYj7c/97d5MTGXFifam5dncXh0Fmwi8DZ5pjmxcEJ+QjkJ8Ilt+m2PRu6lZYtoOkvxrWMTog4od9X1tnjwxx/eJJYjmAe4c4quzmFb1yz239UhCRNUIrXeFq+4OP9kZEY1a1QtImp6YVm0cfMJdA/2t+L4+DwAYRfEsqjZeDp+Yw7fOzOFPZvaxOMG5As57+jl2264248jY1GJaajV50RXkxfXoyk0L6fEOSCPvtmzqU2y+DKHLVDQUEfCSbxyclJcBJ/f1QcA+Juz0sO12OtzFDgyFpW0VzKdFU1rzLy2ud0nSXTaGmxUNKnwoZisPuemEiJRTy+s4K0rs+KFOryT9MhYRJqJnBOeA9Lwy3JMQeWYYtVMR1pQG29G3m+EB0vBsmRPiH4vvtbWiJ30yLSzw6MRvLhfGBwsMUKenXduKiGJ3WbvMAql7MrlLMV3T0wi6Hfj1EQcWUpxMbQAAJI452e3N8HZ2CyxsTsIsKu/FeHEjGRLmclPBGaH5BNqACJ+VzBtFOpiJLlKr62RP9sdKHZq8nh6R7ekjOFu4Zyig5dnARSyOpldXGnbzhBPZySOZ6X6lYqIki+AI+EkTk3Mi78nEByd16MpvD0WFWP0AYhErKRta7edtFViqSxiqUVcmVnEJzb50BMRTE6b233YPdCCyGIG+4aC2NzuK9pByTVd3j5NIfgF5CGnfN8w5eHQlVlxvPDIUiFT+cX9Q0WJTl9/YlCy4ACC1p7JURwejeDzj/RJxvfzu/okC6zcBAWQoveL/VBmApjWzlSvydYIUZdzHpIcRniwFCxL9u3t+s9t0doaCc69wnflK6x80G5u9xXZiAkpaDxGNIIjY9GiNHoK4MrsIjArrRMfOracpTgfzeF37izIwmutz+/qw6nJeZy8GQeoUD8KiJmQvBmDdzICpTMr1bCztxkEQiIPpUAyncWrZ0NFsfFXI6mSZREAT+/oUjx4a+9gOw5dmZWEJcp9E0rJMMl0VuJMVYqU4bfjTOOUR9zwC6CgdRd+TyGQoRoOXYlI2oKvs9Iic24qkXduRhSJ7dCNFHAjBQeZgoMUiJDtcPjxoqTpyhc34d5h6Tt29DRidCYlUR7UfCNA4bTUriaPpv36SH43ifzzQ6MRxfryCyzrD9avR8YiWOEqw5SZXRtbuLN/CuO4FJmr7UzNjLbh3yffVfG7d70mJCM8WAqWJfuZmRlD8aVqK66Sc48nAnl23pGxCPYOBiUEmaVCKv7h0Vm8uH9IV6fJo220wNs52aQ7F1rEn7w5JjoK+Yl1NbKI05MJ0YbPtBzk//3GxWl845NCPZmNFdzzRwZaARhbwK5HU5LwPrnt3eMUtMtDV9TJkHd2nptawFePXBMjYfg6PLyxBSfG46LMh65ERLu+WjKMHMxUobSgMLPW22MRPDncqbhIjISToobLoKZpMrwfTWEsslhEJkypOHQlgqDfjevRlETL/a17uvDGxekikmV/ZilEGdmuRb6YiSZHGfE+PtgO3q9zeHRW1OwdBHj0ziCee1AwuRy6EtEkeoZYKoO9g8ESB9FJC/I4IWq5vGYOFO8w2bx6tM+D8zGK3iYvzk4lxFvVnrm3B8/c2yO2w7mphKxNi007/KIuD26QR9sczO9u5CeBqs15paRAfleVowWfC6uj3kXFKA9qwbJk39LSYuj7I+GkGJq2bygodtJwtx/7thVMA/KbhIrsxKNRAARKWcorOeC18yEEfO4iGzwASYdfmVlUTNxRwpPDnaJNmvcjsEihPZvaJA4w3jSjNDeznDa8s7cZDkyJdn8KiDZoQGrWYu3ID2j2N/8bOSgEv8ahUWWyYJr8j96bEScAv7thdeAzVNXk4Set3CzFYza5gpeOTWAqnhZ3OQAkyTTpLBVDBFk92dg4MqaP+Bg2tnkxOZ9WNNVIfECzwMmb85IkL7/XKd4alkxn8X40heVsDhdCxZdjUAi7qr8+M1U4XI17zkJ4pZpswVf14v4hvHY+JB6WxkwwQmQSEeeJFiKLQr6D0q6Yybx3MIjDoxFkckK7vhdeFLRcCHPy8vSiuFNJprP48ltXcWe7Dz8ZEcYIaGHM3o4v52UrPjGUN9Ey8GY8Nn7lO7+fjUaxNdiIq5FFxFIZOAlBjpv0LLOZyaSm/QMoMtU6yRQe3tgqCSk+cyuBs7cTitm/Wj4HozyoBcuSfTab1f3dkXBSkljCIhvYdndr0AevU5oZyYcGys+zByg8TlJkggGkJxgChSNlT03EFW2hDMyWzms0DG9cnMbuvJOW9yOwQXJ+KoEH+1tFLfj4jTkdOwbBTDHc7ceugVZVspY7KuV2YKaJlDqiIwfhTBYlOAjQ2+KF8tIkJN68dj4kavNa8sh9Ak6HMEmdBGIEE4/XL07nnYtTinbo2WQhUYsRKf8+vdjV34qfJGYUx9iRsahk4chRaZIX0055jXEknMQf/v1okTyAsINQUyQe7C9e1JifaGtHI/YOtmNbp18ge0j7f2vQpyofvyu7MrOIP/rpKB7e2KpqymELi9zRmgNwIZQUzTd8hNbpW4miMFoeFFJFjb+MR/69Fq8Lf/7OLbx+cVowdcq+k8lRaZa1TEZAMOsxmdRi7R8fDEqSAgFhoTg+Pg+3AxjqbMSV/BlJ7IhrxgFqJjAeRniwFCxL9slkEh0d+i5T4GPWGQ5dieBGLCVxCMXTGcXQRPl59iyTUmlrKycLCpQgqUK6PgDc0+Mv0tjYoGKyPDncicNXZjC3zKIThMFDANycS+G9EtehMU2O4Zl7u/OLkXItecc070jk/Qh60MwlCknkA1SPKACE9mHkI36Gwnafnb0vmFak10YGG11o87mxs6cJP3pvBjlZZeWZr6XwxsVpcWHiM0FL4f1oqmiMsYPOBtoaiuR9YXe/4njkHa3/4q5G/M2lxSLCv7Pdh7O3i8c8ACTS2aLjJ9hO6srsIn42GsGBR6SX8LCF6VcTcUXZ2M7s/WhKdERnKXByfB5OhzSCSW7SiKUyim2YpciH5Eqh1dyUAr+ajGN6YSU/vtUX5FOT82JeglK5BFLFQGloup2FtmnxCkdpyGPtmexphc7IUmF88vjMPV1IrmR1h3Aa4cFSqJjsCSFfAHAAQBbAQUrpH+c//zKA5/Kf/z6l9GdGyjVy0a7cVAEAQb8bY5FF0SF06EoEWzuEVVaevPLsfT2KDl6l7EQllPoOG+xZCsWtOSHA8Rvz+MtfTRU9k79H6fcMD2xoxr0bmsTB89LRm2AHgX1sqB2xVEaR9JlTlZm0aF4DYX4E+UBWInQC4N0JgQgcAO7u8ePS9CKyOSHLcs+mNlWCUtrt7BpoxQd6fFjMFibDn7w5VlSX8MIKwgsrGJtdLLn70IMsBb6tYEri8cCGZrQ2uCQmID5Xgzc1ZSlwLbYkfo8A+AJnC37p6IQoU1rmaP3PewfwwFYPvntyUtytEgDJlazkeAIe74WT4gJHIIS3MjMIIGi0VyOporNmlMiKgQL40Xsz+PTdnZLwTwohGbCryYMWr0sMMc7kpNnQargys2ho/8TG/4VQUozwcSq0g9tBEFksvjOAZUA7CPBb93ThBxemFUneSYB92zpkB+gVThJlMrKw063BRhy6EsHVyGLhyIh8PQI+jyTBzO914ncfvkN3CKeZF45XRPaEkI8CeBLATkppmhDSlf98GMBnAdwNYAOAtwkhQ5RS3XuSUCik6ZjgbcvXo4UIEKaF7N7UhlMTceQolWg2PJxEut2Wh47Jk6/0QmsryiOfG4QsRVHdynnPr28n0NPiLToIjH3X5QC2dxXvLC6EEnjp6ER+G19I4t/c7hPT9HktyUGAO1qlp01SQIwXJxC0/GyusKj+6L1p3NfbjNO3i09qlLeVkwh27VMT8/j8I/04NyWcobOsQUgs6kGN8NUWKKX3l9LmZ5LLOB+SysHMJWxMqmX8EiKYF4GCE1j+btZmx6+G8Xsf2Y59Q0FcnV0UZYylMqr+BL7uFEBvk1dC9gAQSy2LCVcsm7kUVnIUv7weK/p8a7BRDANVijwrVXK563MmJ+Q8vLC7H391+jZiSwK1OAjw+UeEfAK2QALCmDrAHZkBAD/Im3h4EAD7tgXxxT39ePVsSLLTjacz+OKejdgabMS3jwk+mG8fmwDy487lKCwSVyOLuDm3hAuhhLhDJER6RpGewIhSPGgElWr2vwfgv1BK0wBAKWUhGk8C+H7+8+uEkKsAHgJwQm/Bbrdb9ZncWZKl0m1YciWLI2MRicNFDgLBvlkq2/TxQUEjfufmvOR4A6CgPSmZLfQQfl/eqVcOdg+04sG+FknyVw5QdbAxc5DSzuDmXBo359ISudgxCs/eJ0Q+yG9b8rnV78bMAZLFgX+HHog2TVoIOS3Vlg5ok7Raf5RDNpPxtOK7KASTGMv4lS+4gCATcza/dj5cRNosNNZJgLllioOXZvGdE4U+ZuYTkLzMJeqqtLgCpCiLVg/kB9FRCLugXRtbi+zWpWDERKYEV978JG9jdhPYzt5mbGwrKCQsIoadg/M9zsEtrxfLZufPAAIVTKns6kn2Ti7XC5mcEFm0ud0nmZcMLI+BHcmsJwNYiweNolKyHwLwzwghXwOwBODfUUpPAbgDwEnue5P5z3SjtbVV9ZncWcKDArpCyBwECPgK56UzJ9bVCNOghNhm+QFcBAVbMnufEti2TWsy3ppPF0UB8O/pb/OqEmTA5xZNAUoDSw+KHFIohMfxtsQjY1FJOxMC7BsK4lr+GF2z4czfkwtox3wzc1Gz14UTGtFCgHI/MKe50bbTMhcxk9j+uzokZxCJv4UwPgWtT7rwsgPbmKnglzeX8MubxQsGIxgHAbr97iISLqVoBHxuHBqNaHxDP3IUODE+DwcptOe+bR2IpVZwYnxetR67NrbgnXxEECDtH5cDcBAi7jgIgKGORgT9biTSGSxnqbDbiRT7M0AgZizzv6cQHMC/vpWAw1E8rxk+vDkg+s+uR1Pi7pRC2ClcmVnEDhVyZou93CEv/46RDGAtHjSKkmRPCHkbgJLh6Cv537cD2AXgQQCvEUK2GKnA9PQ0nnvuObhcLmSzWTz11FM4cOAA3n//fWzYsAFOpxPxeBydnZ2IRqOglOKuoF/T1qeH6Pdu9GI6XtjmMVMPX4ZS/DbTOPVgW0cDIovLmF5UpnwK4EMbvGj3ORBdceLc7QXMLVOxjp/aFsAr74YU5RmdiuHcZBNuhmfx8c2NePOaMTOQA8BvbGrGL24UmyMcBHhmyIeVlWX857du4OgtqRngoV4vboZn8W8e7Mbr56cxndIZY6oTPT6CqUWqTfT5rfPl6SQ+uqlFt9mMr+mnhlrQ4szgZDgLQrNYXM5hMqktCzOH7buzGT8dK/bpzCZSGB8fR0tLCx7scuDwaP4cGW53ouS7afMQfOGDATT6KH5ye07XIpqjQIszh3DJbxbgIMDs3DxGZ0ofF60XlPuPEGBHcxrBTUGcmphHJle88LgI8Om72vEbXRTvTmfxTxMpiWM9w6nLBMJZSD3eDAKNXpy8mUSOAmOzi9jgdxTVZd+AD2fHp0UfBAHQ5XcinBTMPDmgyInP4xfXY8hdg+i7UhoNkYUlhU8F/PpGGMdvL6s+Z5Lx5rqjV27Dv+RBT08PQqEQ/H6/yHsrKytoaGgApRSdnZ0Ih8NlH3tckuwppY+rPSOE/B6AH1Lh7qx3CSE5AB0AbgHgD3Toy39WhK6uLhw7dqzo840bN6KtrQ0AxNvV2V2M/QD2bUsWaU16QSlw5Ga6LG3YCC7NLikuSmIIJgW6Ai3YvalNcOpRh6gNUQr8OrSoutUdm8/iy4fHkcufh2MUNF8/tWfvL7rwN4duKA72d6aE3Yb7+hKeHO6UJFhpYUugQeKsVEMpwvW6CJYzBQfoQra0NutyEDzU3yKGoBIAl6MZcSfH/BmTSe1IJ/6y9mT2hsRJCwCnbiexa1M79g8EsDsQwIutbWKcNzusTame9/e14UyU4NCJSUPjMut0o7vJIdHu1X5OIIy7d0JSoq/UpEJQUIByFLiQ8AAJIVxYjkQ6g+mFZXz75BQ+fXcXBjoz+Pm4uqLCKPt/3VoG5ZSOHAUmFwrjhPnqelu8+BZ/yJuD4LP3bZCc0SSvOztNlpULaCuMoaSy69FBgIzLhyzVJnvk30mQT6jctgEDec2et88HAgHMzc2JPCh/Pj4+XvI9PCo14/wYwEcB/CMhZAiAB0KqzN8B+FtCyDchOGgHAbxrpODlZe0GY+GSRp2nAESSWA0ovYY3Ab12YRo/fG9G3C7ykNu95dBzxrtWvVQvA6Ha4aTstctZISHo0S2BItJTAotgqLTp05lCCTkKnOTqyia93+vElZnFwiFmOSq5iUq+k1PzZ8jBsj1Hwkn8QsFhmaPAt45N4NRkHAGfG1uDPkwvLOPwaMEUpkSuetpPCdeiyosnu5QjRwtjXY3sPrI5gF9ci5W0/ytBKTb9rcuRouAAj5PkFYOCue2lYxPY0i4NS1WCnnH+yEALelu8RX6ITI5ic7tP1SRKARCVBVgNau0kjEV10xWPu7v98Dgd4iGIaijFg0ZQKdn/DwD/gxByEcAygN/Oa/nvEUJeAzACIAPggJFIHABIpUqfscKcpzeiKdxOVN4oBMBHtwRwK5GWePL1wuMElsvIgVAazOYaRoyBbcn1fG88uoRfyIjKAeAT2zsQiqclDsKbZTqjS0HioANwLrSA39vVxyVICZ/73U5UutSwOO9YakVVG6aAdsaxwSqUs0B2Nrnx5Y9uxmvnwyWzn//xWkx3+d1Nbs1dhNLYoRAyrJUWNLXFyijeuRlXVY6++U/j2NreqOr/MjLXWP6HajSUznJYrszF0ILqBTiAPh7Ui4rInlK6DOBfqjz7GoCvlVu2Vnyp0mmSZoAC+OX1GD68OVAW2ZdD9EB5k7maMFIXJaKgEML7zihGghhDo9uBRa3UZAWwDE8l7dlgUUVgcd7lmM7KRTljI5w/nbOS8pXGZVtDsUNYD3KQZiubAb5+Wjv1m3PpsqPe5MgB2L2xFZPzS7qjyxTLyde3VAatmXH2xR4OiyAUCqk+4w8ZMhtZWv6W2iga3Q4xbb7W8LsdIoHJicznUqc2tbofL5FVLIdaAs59G5q4y6v1I6uwNTeTbMzoMwJ1uc3A0RtzSKSLk4v0wEEEc5icIJTyQbqb3FA5KaOqMGR6oeaR3fHxeUwqEH1rg3o4shqULozhocWDRmFZsvd4PKrPhGSn4tFVg/FW0XsXV3Kr5jsAgGaPencnV3Li5JFXKZXRX8mAz1me01zhvQDQ19KAL+w25/IGLXQ3uQ33YyXjzeciuKfHj0cGWrFZh91aDzY0S2OyW70uXX4IJWxqa0ByJYcvfKgf+7cHsbHNq/pdszX2auEjWwLY0l6QgwBo0FBktCDfIDoJ8PFB49cHlnJHaPGgUViW7Jub1Vc7dmQsX3kC4J4Kzp+uBFbQzPUgsWy+J0A+VZwOc4fU+9EU9t/VgW2djaaWK0eX34OPbgkY+k0l/Z7KUFwIJXF8fN40u/Wm9kbsHmjBts5GPLOjSzHjFRCISWnhD/icYn9eiy3h4OVZfPv4BPxup6Imy5ClqxfwUAl+fi2Ga9GCHEMdjRho85miJG4NNsLvNa7ZsxNd1aDFg0ZhWbKPRLSTPuLpTJFjbrSMIweqDQLgGdl2mCUt1RoPbGiueADI53jQZ17GHwAsZ3IYCSexs6e82GK9uBBO4h9NNN/5PY6yTRvl9AmBcMTEifE4bkRTmIwrhxZvCTTAQQgWFBb+WCpb1J85KkSMVaom7B5ordnOWw3sCBUz1qmdvU24MlPeLuqmRjhyKR40Asueesli69XAbk7iO0rrMCc9qIajtK3BCb/XiY9sCeD0rTgcDuCuriY0OB0V+QZK1dVJBG1D68ydc6EF3N3tx4USp2gaQTln/GjhQjiJP/j70SK7favXifm0ece/Aub2fdLgLorvzxyK+7dUf1MUonzSWSocqaCA63liofkyW6rQjkq4oPNQwbWK1y9Mly2f1tlEpXjQCCyr2ZcKORru9quabdwyqbr9+rRNs2ynPGJLWfzlr6bw82sxzKeziKWyOD4+X7ETmI8r58FIkdLSxJvJUUwlqhMOaSaU8iJWg6BWCwGfE36ZWUVO9E/v6ML+7UHd2rHaUsPKdUBIJFsoN4TMIKphQrQSKlnI9m1Tt/WbGXppWbJfWiptx2xuUN6YyMPrwjodSGbZTquBBhdRnOhK225AebIrOaNmFY6CtaEMr5NgQ7MHrQ1OU00SsVRW0azCQCFojrFUBh/dEqj43RRCrPj2Tr8ptnarmWdKwUr1fXRLQPHqTAY9PKgXliX7UvGlI+GkZrLIesNSpvJQ0yUDUTU8mMOv3pHOUtxOLGN+qdi2XW2wRK2fG0iA0kKOAvOysEy1m8b0YEOzeVEj1YaDAFvaG4p2U6sNJwF+8+5Oze/UfZw9ULgsu1bwlHC8W0l7qAROAgQb3ZiMp7Gjx18RIdiwDiiACe5Ya6DYdqw3hoACpmSwrxayVNjFG/WrANrhy+XUg52wqQYz4+wt66BtaNC2n8dS5g+uDc1u3E7oM/mUMnVSlH98gpXA7tO0sfYhP4hOa4fAzsy5FU9j1KSIFTNRq6xzs30PyRK+p1I8aASW1ex9PvXLjwWYq2ESACFZGnjAZzxulsdaJ3ojsPV960PPiaMMFEJculmhidWAlcccgWAqKsUhP7g4jRGNaLjSPKgfliX7WEw7WiXgM3dTQlGczRZfqiO2rhC9FrbZVrL1Nisd4tEtAXQ3lZ+DYKb5YL3AqosQINTtrq4m7B7QDp3MlTDllOJBI7CsGScY1E493hpsBGBewoES1kJWoFVgZZttJVtvsy7i+uX1GHzymGADMMN8YLUD9ypBJXKsVjtcmk5gPFY6tFnrbJxSPGgEllUXEgltx0W8zAOezAS7fIBdyVYuKuCAmoEAaFqLFVfAavicsxSa4ZWrgUoJcr2gX3bOj9/j0CVfR6Mxs+61aOkLkkqNvVI8aASWna2lDu0XLgNepcqogAL49N2d+MS2YNF230goWqXH7q42xISctVZxFWQp0Oiy7FQoC04C7Ogx76yoamvCqzmVW2X5OcnlnC755pfMH+8U2mYcMy8vsewILxVfOtztxzc+OVRzW+bJm/MASNEK7qr1SlRl9Laon4LIsJaaYDFT3YXLSQT7v9thLgmr4eGNLfhgX8ua0shXo64EwMa28iJc9NwLbATsjCwtM46ZcfaWtdmHQiHJfYtKuB5N1TwNO5JawVuXZ4s+n0+VeY44yrulqtnjRMJA+I/bCayU6X/OASUvbiCwfR4Mj24JYEdvEw5diSDod5d9xrwRnByP48G+VnicRHLJDyNUq3UNheCE1jOfnQSqd/mWgstB8jeW1Q781Zk7e5s1ryXUw4N6YVmyLxVydPDSrORi4VpBLTFjoUwmLTc+32jUSJPHiViqetFGtSITl4NUdDdvNXB5Jlm40atYL1CEHieiywGobUhyAF4/H8IDfc0YnU0hklxRva/AKtCruOVo+U7W7V0+vDVa3cCOUqAAjo7P4ekd3ZpED5gbemlZstc6tH8knMTLx5Vvi7cKsmVuOCjKi883ejBYNYm+lrAa0QPlRSrpkaKU5el2YqUoSXA9ROSoXXSjBxdC1jgG/XZ8GS8dm8BUPI3fffgO1e/VxeUl8/PqWZvnphK2icCGjTJgTxtr4fUL2klVWjxoFJYl+44O9ZPghLPsbdhYn6imY5udY19rEAA7eqp7+5hZaPY4saUKx58DQn9onfOlxYNGYVmyN3NFs2FjLaHau9ZSh/itBgiAixYxqZRCYjlb5ePP1TvcTB60rM1+ZUX9QLIjYxF7O2rDRpmwwr0v6yNDwxwIpwEoQ4sHjcKymr12fGltN6IWUIxs2LCxTnA1on4bVd2fZ793sH0Va1IMCyhGNmzYWDdQt1OYeZ69Zcne71ePPx3u9qPZCoZHGzZs2KgQWkleWjxoFJYle6dTm8xdNtfbsGFjHeBcaEH1WSkeNALLkn08Hld9NhJOYm6dJgXZsGGjvpBYUj8+Q4sHjcKyZN/ZqX4Rrx2NY8OGjfWCaEo94kaLB43CsmQfjWpdKG6FtBAbNmzYqBxabKbNg8ZgWbKnVF13r3U0jg0bNmyYBa04ey0eNArLkr3W9uV6VD0u1YYNGzbWEp57SP0gtLow44TDYdVnh2p8RKkNGzZsmIFHtwQ0jznW4kGjsCzZNzU1qT6zQ+xt2LCxHjCvEYkDaPOgUVRE9oSQ/48Qcjb/3w1CyFnu2ZcJIVcJIVcIIR+vuKYcmr1uM4uzYcOGjZrg9O2E5hHHZqIisqeU/nNK6X2U0vsAvAHghwBACBkG8FkAdwPYB+A7hBBD+vjCgnqigQ0bNmysF2gdcWwmD5pixiGEEADPAHg1/9GTAL5PKU1TSq8DuArgISNldnd3qz5bjTs8bdiwYWM1ENOIs9fiQaMwy2b/zwCEKaVj+b/vAMBfEDuZ/0w3ZmZmVJ+VsnPZsGHDxtqBenilFg8aRcnz7AkhbwNQOmfzK5TSn+T//SwKWr0hTE9P47nnnoPL5UI2m8VTTz2FAwcOYG5uDj6fD06nE/F4HJ2dnYhGo6CUwmnnz9qwYWOdwEuXMT4+jp6eHoRCIfj9fpH3crkcJicnQSlFZ2cnwuFw2U7bkmRPKX1c6zkhxAXgKQAPcB/fAtDP/d2X/6wIXV1dOHbsWNHnd955JxobhWSDQCAAAOLfntNXABi/xNmGDRs2rIZ7+zsxMCBcPzgwMCB+HggEsLi4KPKe/Pn4+Lih95hhxnkcwGVK6ST32d8B+CwhxEsI2QxgEMC7RgrV2r7s2xYsp542bNiwYTnENXyQZppxzCD7z0JmwqGUvgfgNQAjAA4BOEApNXRMZUtLi+qz/Xd12OfZ27BhY10gqXFPpBYPGkXFZE8p/R1K6SsKn3+NUnonpXQbpfQto+Vms9prQ7PXJnsbNmysfZy8qX6peCkeNALLZtAmk+qJBgcvzeJ2wrbZ27BhYx1A49hLLR40CsuSvdZFu0dvzK1eRWzYsGGjitjV36r6rO4vHL+z3beKNbFhw4aN6sGvYZKuiwvH3W7182+0GseGDRs21hK0HLRaPGgUliX71lb1rc3O3ma4LVtzGzZs2NCP9zXu59DiQaOwLGXOzs5qPs/ZSbQ2bNhYB9AyS5fiQSOwLNlrrWhHxiLI2mRvw4aNdQAts3RdaPbLy1qhlfaF4zZs2Fj7cBLBLK0GbR40BsuSfSqlbseyLxy3YcPGekApC4UWDxqFZcnezPhSGzZs2LAqjoyp36ld93H256YSq1gTGzZs2Kgm1M3SdRFn7/F4VJ9pxaXasGHDxloBgbZZWosHjcKyZN/crO600IpLtWHDho21AkeJWBMtHjT8LtNKMhmRiLoda8+mttWriA0bNmxUCVmqbZbW4kGjsCzZs9uplLDZPhvHhg0b6wQtXvULA7V40CgsS/ZaIUe2g9aGDRvrBX9/Sf02qroIvVxaWlJ9Nh5Vf2bDhg0bawm35tOqz7R40CgsS/Za8aWXZ8070N+GDRs2aok7Wryqz+o+zn57h38Va2LDhg0b1cOnhjtVn9VFnH1DQ4PqM5/HstW2YcOGDUO4GlG3y2vxoFFYljV9PvWIm1/fsh20NmzYWC9QPyBHiweNwrJkH4vFVJ9FFldWsSY2bNiwUT3sHQyqPtPiQaOwLNkHg+oNEGw076ouGzZs2LAqtHjQKCxL9omEuqkm6LfJ3oYNG+sD3zp2U/WZFg8ahWXJXuvQ/vmlzCrWxIYNGzaqh1BCnevq4vISrfjSjH0noQ0bNtYJdvWrXz1Y93H2yRX7iGMbNmysD/zm3XUeZ68VcrSxTT3jzIYNGzbWErTO+qqL0EutQ/sfvdM8D7UNGzZs1AoE2heO18XlJfPz86rP/n5E/ZQ4GzZs2FgraGtwYrhb/fgXLR40CsuSfUdHh+qz0IJ5HmobNmzYqBViS1kcvDSr+lyLB43CsmSvtaLd1WkfhGbDho31ge+fU3fC1oVmv7KifiRCT4t5diwbNmzYqCWWMjnVZ1o8aBSWJXut+NJYyk6qsmHDxvrAAxtaVJ9ZJs6eEHIfIeQkIeQsIeRXhJCH8p8TQsi3CCFXCSHnCSH3Gy1bK7404FO/s9GGDRs21hIG2tWPMbZSnP1/BfCfKKX3Afi/8n8DwCcADOb/+zcAvmu0YL9f3S6/NdhotDgbNmzYsBxcDqIZeqnFg0ZRKdlTAGwP0grgdv7fTwL4KyrgJIA2QkivkYKdTqfqs3jaNuPYsGFj7SOTo7geVb+8RIsHjaJSsv+3AF4khEwA+G8Avpz//A4AE9z3JvOf6UY8Hld91uK1zTg2bNhYHzh6Y071mRYPGkVJ1iSEvA1AyUvwFQCPAfgDSukbhJBnAPwFgMeNVGB6ehrPPfccXC4XstksnnrqKRw4cAArKyuYnZ2F0+lEPB5HZ2cnotEoKKWIJomRV9iwYcOGZeGhyxgfH0dPTw9CoRD8fr/Ie83NzZicnASlFJ2dnQiHw2hqairrPSXJnlKqSt6EkL8C8MX8n68D+PP8v28B6Oe+2pf/rAhdXV04duxY0ecNDQ1iQkEgEAAANDYKtvoHPEn89dkZjcu8bNiwYWNtYJl4MDAwAADi/wMC701OTqKvr0/8jH8+Pj5u6D2VmnFuA/hw/t+PAhjL//vvAPyrfFTOLgDzlNIpIwVTqk7l16Mpm+ht2LCxLnBnu/phZ1o8aBSVGr//DwAvEUJcAJYgRN4AwJsAngBwFcAigP/daMGdnerHfmrZuGzYsGFjLcHvVXfCavGgUVRE9pTSowAeUPicAjhQSdnhcFiyZeGxZ1MbTt8y77ouGzZs2KgFHET71EstHjQKy4a1aDkh9t8l2PL/4t1bWFhRTzW2YcOGDStj18YWzVMvy3XGKsGyxyWUwv67OvDhOwO1roYNGzZU0NrgxIbm+jnHylFGkOCDfepXEpoNy5L9wsJCye+klm2t3sbqYfdAa9nk5ZQRgYaZ1jLY1lFeprrf44CTAPNLWdzWuEx7vSFXhi+1VIKoHh7UC8uSfXd3t+bzg5dm8fNrsVWqjRTucpZwC8AB4IEN6vbBegRBMRErgX1nuZwZDSAr+1l6DVyjfGV2sazfERTLa0MZpRJES/GgEViW7GdmtG+j+qsztxU/JwC2BIoPFnp0SwDuCqV1QNDuPq1xQbCVkQNwVuO+SyMI+JzY0ePHjh4/Gl3FbNnR6BKfO4jQL6sx2Iz08YZmD/77p4awb1vpCyKyFDg+Po/ZpL4jZ9eqQlAJmj1C4y9Uccctb1UCYEePv+jz1YCTADs07O168MrJSYyEk6rPS/GgEViW7AlR776RcBKxlLJqRAGMzy0Vdb7P48TnH+nHto7GsoTe2OYFIcCJ8Xm8dmG6jBKsAb0al9bk2dLuRSyVxYVQEhdCSSxmiguNpTK4GEriYiiJHBX6RS8FdDe58cCGZvjdxqfwg/2t6GjUF3cQSizjtfNh+CvVAhSwUuYOYDWwo8ePjW1e08tNqJB8g0u9fVsN2LMIgPs3NMPJKQ9OB8GFUBIUKMvE5qiAsCkFmhuK6+8kglKox+y3kqOaF45r8aBRWJbs29vbVZ9pNQ4gEJp8qt2cW8TLxydwZXZRN+kwOAhwcy6tWC6DfPLo6SKCyjUDtXLdDn11UMNHtyg7vx0ArkXTJX/P2qocyutrbcCfPbEV/+q+LsO/PT4+j8iivoPycvnvv17G4r1F41haJTSWuaB49diYDOJCKImJudJ9CCg7Hf0ehyFnpNrlHC4HwceHgrrLoQBO306ICsumQAMy3KJajn+AUmBjoEGXKU8OQoCAz1M0z4QqUQR8LsyXsMk7iPapl1o8aBSWJXut7cvO3mZD22QC4L3womE7YsDnBEFpx4uTAJ++u0skWL1VowDu6vLD6ySmbEMdBHhmRxd+54O9+NhQ+XdXOiCcsf2ojPAJhAlWKQi0tbA9m9oAAKFYeYdAURQWOicpff9BOQvSp+7qxP7tHbr77ZPbO4rak0GrjHSVjN96S727u9hJu39bB37rHuMLMY8NzR4ceKQPPxkpz0xBAVyLLVVUB0BYcPYOBvHwRmNRMQTAwxtbsTXog0s24SmA4+NxHLwcQbKESevOdp9m6GVdmHFaWtRvbxGgfxLc0+1HOVnHsVRW11uyFDg1OY9P392F3hYPPtDbLNqpeWwJNBR99n40ha8/MYgntgdF+3Y5aPU68YXd/fB7nWjxunBzblFSd73FEgBuJ0GL14VfXpc6wI1MMK33UQBhjUvjp+JpvHo2hKzDVfYiyMs+X4WbzQ6NRrA16INHp0qYXMni2Pic4rP+Vm9ZmuVqYH4pW9QHfq9TM+tTCfIyQoll/Oi96aotZnqwJdCAF/cPAgAii/p8MY1uB7YEGuByEJy8OY/vnJhAjtKyx+m+bdo7m9I8qB+WTarKZtXDFc5NJaBxbWMRRqaTcJDqRggcH48DEDTR23GByJwE+PDmAH5xLYYcgBtzS0WLBzsX4+2xKFZyFE4CfGJbB3oasggtOXB4NIpMjsJBhEtbrkVTivbgRDqLbx+fUN2F6BGdALj/jmZ87v5eHBmLlmyvDr8bvU0eXFBwMBEAXU1uhBeUJ1GOCt9ReoUen4hTZ39Wq8+vzCziRjSF53f14e8vTZc0bcVSGVU7fl9bAybjpc0qzR6Hql2cR5PHoegkJQAeGWhBIp3BhZD+SBt+7vCXbTjJlGb7Cu9rxYN9LTg1OY8T43Gxv3MQTKMMTgLc0eqVfFYOnCTvH6KCJqtlShyfW8LxG3P4ycgMljUE4cdpaiWH8bmlgh9K5/gK+JzwOB2ILGaQyQmLw9M7usQEUTVo8aBRWFazTybVPdQtXpekA3f0+LF7oFVVmCwVBpLa6rtFh81uQ7O7yHRUSqujAG7F06KPgBEcD7/XiXNTCazkKHJUqGtXkxu7Ogm+uGcjPv9IH+6/oxmfuacLN2LKRA8Ik6cSnyDJy/O5+4U7Zg6Pzqp+1wHBlvwfHt2MR7e2K7Z7DlAlekf+XY8MlJdQsnugVXXbrbSj4uEk5SW/KGElR3E1sogbsdIENTkn7Ijkr97Y5sWDfS2i01Gpats6G/HMji40evRp0wvLObgcBA5SkJe1+YN9rfhgX6tixJoSbs6lQblxSynFkbEoAGhGMTkI8Psf6sdX927B5nYfTk8mNBWOhze2YEdPU5FJRF6mEpizdltnIz5zT5c4l0tNhywFXr8wjeUsVf2ukwikLMoP5XnmgOAnczlIUT8SALsH2vDE9g48dXcnHrijGb//oX787sOlr/jQ4kGjsKxmr3XRbjydEVdbB4AP9rVgemGlpOPV5SDIUkF7zuYEQnIS4FPDnbgaWUQslUFkcQWjM1ITiAPA7UQxce0ZaFON9RecpARBvxvgePOeHj9GwklQKphLmJbkdoSwkqNw5zWnnrYARsJJvHJyEis5irO3E0WDTJjIxFDkh9shRKy8c3Ne1MqcRJi4ewfbMdztF0woKkU+sKEZ925oEuv9yslJww7vD+R3DwBwejKuOdkYHAAGOxuxbyiI/Xd14KWjE4rf+617upBcyeJCaEHUEgmAoY5GbM0nCR28rL6QMWxs82JiLi3Wq7vJDa/LgdYGFy5PLyJLhb6KpTJF/cITA8PNeakWyzTDibk0vnNiAhQE8lZg2l9vixffPjah2M5K7wKAXI5i10ArTk3ExZ3hk8Od4niCxs6KAGhpcGJ+SdAq+fdmqdB+b49F8ORwJxxQj7I6NTmPozfm4HU5NMcoAXBqIoFs3hyiVC8C4WiBd27GJWPzmR1Cfx8ejWJsdhFXZwu+OT2zgqKgIFCFAIzP3NOFyXha8rlS/QY7G/F7u4SjiM9NJZBMZ8UdKgXws9Eosjkq+pMuhhawuYS9HjD3wnHLkn0oFFI9AGhnbzM8Tik5HhmLaJY3MZeGyyGYSADgrSuzABUm3cvHJ0AhkPPzu/pwI28qcRLgY0MduBBKSLaXjW4H7gz64PM48cyOLhwdn8NUfFkcADt6GvHBvla0eF24GlkUFpkchctB8NyDwmp+biqBnb3NYmd//YlByWfj4+M4MuGQ2DQZSRAAuwZa8cy9QsLFN/9pXFK/Dc0ehBLLRZOQAPj8I/3Yf1cHDl6axY8uTgN55/Lmdp8Y5bSztxkOEhJtkfzAXlgubCvPTSU0t79KROB2CLsHXu7vnZnCmVvFml9nA4HP6wEIsKu/VbQTv3o2pBguSSDYxt8ei4r1Yg7zfduERWIknMShK7Oa5ge3g+DTd3fhOycmwI5eCi+sgAAIO5fx+Uf6EE9nsLO3Ga+dDxf9Xo9GyS8mwjuKf0UB/ODCtKYp4pGBFmzr9COZzuL1/HcZ0pkcslQgmCwFfn4tpriwPrChGWenEuLYcjiISPQM8nGwnKV44+J0Uf8y7ZuAmTYLn7My+LIcBNi1sRUnb85r7kwdRDhaIODziIu1AwX/QZZSxZ0zIJhQ5peyquUTCMqOfJ4DwBsXpyVjhdX3xPi85Htjs4u4Hk1h/10dosLEy8mIHvnPWMhlKbLX4kGjsCzZu91u1WfD3X6RHFu8LpybSmBrsBFuRwSZnPLEoAAyOSCUSGPPpja4HYJGzKJtKIQBfPTGHJ4c7sT70RT2bGrD5nYfDl2RaoLpTE6MMXc5CA480idqTG4HwXMPCiv8n7w5Ji4aT2wvaM5MBh7Xoymcn1pAi9clkP0CLXrvQKABN2JLoFTQiB/sa8HVSLHttbnBhSmFMDQK4GpkESPhpITIXj4+AQcheW01hA8NtIkhbfK2fD+awlhkEU4yhQf7i00pbLI7CcTniXQG80sZ9LU24Jl7uyWyD3f78bn7e3ExtFDkrJtdoqBLwuS7OVew4zMC5ycTM0MBQr/yEyubX9AB4Uylz9zTpeoX2JbX0Ia7/bgaSeHNy7OSspazFPF0Bs/e14ORcBKnJuYVyymFibm0Lj+S1q7J7RAIMJ7OYPemNvS2eAXFJb9rvLPdhzPc6bBqCWELy1l845NDODeVwPTCiuLOp7+t2J7Okydr/+d3CQvhP1yNSr6fo0KfPX1PlxhEwBZMADg1oU72jryS88rJSTy/qw9eJ5EoegDgJCFQKihU2zsbJX6kWCqLLe3C3FF6B4VgOv303V146Zh0xyjvn8FgI565txunJubBn8GYo8C3jhXGGK+QCmYlgkx+XAomH+2QSwYtHjQKy5J9a6u2PZcRxp+8OYblrLBN/Ux+C//W5YhkkvDb5jO3ErgYWhAHZYvXhVdOTooaz+lbCZy+lRC3Wo8PBiWEt1E26DM5ilOTcTw+2A6ASEwhrEyat8OrreIHL82Kg4wd3RxJOosG5rVoIRJmOUsVt/ZuB7CzpwlXZpQdcLFUBt87MyUZqFkKZPPhSstZqmiaIhCimt6bFpKkclSIUVf63nC3H++FkzgxPg8hn0ZYSMKJZTxzbzdGwklxFwMIO4Tnd/XhamQRh65EkKWFCa4EZjd15J1xbAe2d1CISX57LFK0cGQp8O1jE7gaUb7c2QGBIBnRA8DewfaiXQBFIcX9yFhUDBTQckizHAx+3FAAXpcDi/mOIAB6WzySHaIWducdnwUlI4SvPzEoknaL15U3DxWjW1bP9/MXXj97Xw8OXiomegIh92Eqviw6F+/u8YvmLLawB3xu0TRxZWZRcXFIruQUbdUfG+qQLKzy3wGCNhxPZ4p2wUIGamFJfnRrOy6Gk5KyrkWX4CTA/u0d2Br04WpkEYdHo6I5jt9lH7oSwfvRlEQbZ9i3LYjhbj8+NtRRtChSAN8+PiG2AV9PAGK/CAqaPsdRKR40AsuS/ezsLPx+7S0OMyMw7e2Ni9PYt62jqIMoBYY6G0VbfDpLcWg0Ik7sze2+IlMC22oBVNwFuB0Eu/pbJVomAJy8KZCe20GwNegTbXZ8WfIzMA5emsXRG3PYs6kNh0alJqhDoxF8epMTLgegdoKzklZIIJhk3le5rd5JINpw1UAIVMNUl7MUTkJAqbqNPUeBi6HCRBN2WoWF5Jv/a1ywgVLA5ZgCWwgIBCc6k4MAcALgjQm8n8bNaZH8RAWA53f14VvHiokuB+DNy7NFjr4HNjQDRIjv58u5Hk0pat5XIymMhJM4PFogJ5eD4M5gI8ILxQsgM5P90U9HJeXxRO9xEjy9o1uy49LCts5GxNMZcfyvZAWzwLP39YjKhlrE2gf7WhBLZcTFOpuj+N6ZKXzu/l7Fg7lY5rjLIUTXnJqI471wUowc2xr0iQrToSuzeGF3P565txvvKo61wt/8or93sF1xkebhJETcyfN9fm4qISbxZang03t6R/HujRXNImD2DgaLytp/V4do7pNzwu6BFmxu9+HVsyFsDfrgdZKi+uYoRPMM+4+BLUxsgX57LIKvPzGoacrRw4N6YVmy17OiCbblQvhXjgJXZxcViXBr0IdrkUVxIl2ZWcSXDo7hxf2Doinh/FRCfM62WnsHg5JBoZS9y8bzcpaK9n/5un30xpy44ss1eXmCkccJ3D/QgRd7enFkLCJxNgLCRN83FFQkhjcuTitqxATA1o5GjM0uivXrb/ViYl7qfOpr9eLWfFrSfoxkR2cXJRNe7nRj5hW5Jsz/m3dUZvK2avYdueb72JZmeDwexBTi5HktkoGRx7RKFBArVy7br/P26l/fFvqWkYHajWg35xZxbsotlkMAfGyoHXsHgyLBOQD0tXnR19og1vMbnxzCkbEI3pmIS0wqvS0e/PGHN2G424+fvx/FhZA0AsNJBBPe9fzOzpPPgzg1OS8JZeQVCvncYHA7BJIDpM5xtuN9clh67hMzcwJCf92IpURzBCDsZvhFh5nMvvHJIRx4pE+y+3SSwrsPXprFy/lQYY9T2JU8v6sPLx+fUFxgCYAH+5tFoiSYwtYOYR7s7G2Gk0yJUUPTCyviLq/YXFcoXE7GPHjzIm+K+dLBMdHhLThvlyS+CZeCeebP37mFo+Nz2DPQBr/XKUbe6bHb14Vmv7xcOvV5uNuPF3b3i4OGJyRmSkDefrl3MIibc0uSiVTc2AQENO/Q6xC19J29zXj2voJX3OsUTDQEgtYjTnrOXCQHm0xff2KwSJOXa1Mb2xqxvLyM4e42URv40sFRZHKAywHJjuQvTk3iYj5m2sG9n5H5ZDwttsG+oaDofHY7CHb0NmNiXrrVvjmXLgob480LWZoPBby3G0fGIhKTGclPABa3nJ8fmo4xp0M5mshJhBjjrcFGiZmN/y2vGY2Ek6JJj5WbVfE7AMKkzOW3MOz1OSrdhqvdiHYhtIi7OpvEXY4jL+j1aAoEhXExFV/G5HwapyfjYj3liz0APL2jW5ThPVnOAgHwwu5+UdtkpgDWJvz35OOILbwOIkQp+b3OoqAAXntdyVGcCy1IypQv1iyHRG535hcWSgtHmlCuLGYCGQknJaTOdiX89+Vt4HESBHwekSgBQWG7MrOIZ3Z0AfloJhYtdHh0Fh8b6pA6g1FYbPSAmWKOjEXws9GoxGzJWxL4d3x8qF1C3n/+zi1xwXntwnT+QMZin4Ma9PCgXliW7FMpZVOEHPvv6igyw2SpsFV97sE7JDYz+UTirwQTtoKUi3ygEnsobzJgtrjphRUhqgeF0LDTkwnR8SvXcNniEpQd1JXMx0WziJ29g+1IpQr2wOFuP17cP1S05QSA0ZmUKMtn7unCj96bQSZfzh/+xoAoG/sdi7phcr89FhFJgzc7AQVzydM7uiUOaLaVFiYYJyMVoiP49nlTJczRQYAv7O7H5nYfXjsfkmhHO/L24J+PL+IfxxcldZK0Z7awWPMmPQoAOYontktNDLxcHx9qR1eTB7+ajEsUAH4bzmv41yIpxJYKZHry5rxo083SgnmIT7ZhfpB0luLIWERiKmDl7tnUJv59biqhaEJjTni2WLx6NiRxQjMyZH3KTBC8GcfvdYoKC28+kWuvV7ljjR1ECO3lTTFMkWDhs0wmpnTxIcXXoynJmNoabBTllDh3uXnodoQkphGCQqIfoOyPOTo+J7Y1w0oOiKWW4XES0af3wu7+ktEvcgx3+/O+meKOEV5J4XJAVMTki8lRWdb05dlkkc9BC3p5UA8sS/ZG4kt5M0wmJ2iFrCFZY756NiSZSAQC2TBN48qMMOGZxsKiOnJUmKzMBux2hPD5R/owvbCMWCoDJxHWdbeD4Jl7e/DMvT0S7YsNTIKCJrSzt1mSTUhQIB9GpMMdLRiQySgfGCwZixFcciULOWXzdk2lctjA4+OCGQbzsemb232S6CdG/CxhR55dyd4xEk6Kiwm/C2JEz0juq3vvlPgw4mnhxMyCFMqxzTkI2t1IOKlotuhqcovKwJGxiMQht3dQ0DKnF1YkZE8gNYew3//RT0cl757kkuVYPXmnsXyx/9loVHwnK1eePckiOHgyowAOy34rnA1ViPRgzmn5DoeiWAPnnzvIFF7Y3a+qvOzbVjBh8v0uD5/l24knsXNTCUm/sUWLySkn4ZFwEg/0SeeGS/YuJX9MsNGNmYUVcafGEPB5JOM2ns5gJJw0TPjykcd8Sm4nwdZgIw4jCumJTAVs7/CLuyEAaHA6cGQsKonM00Ldx9mrgw2t4kZXGmBsa/ylg6MFWz0RBtTmdp9IVAVtX9DO+cHGnFTysEo2cE/ejIsk8DwX6fH7H5JqQkwj4MM1X9zvURwQ/Hae3xIKzs6CXZoRPCuTRWzInUZMW+QhnG4phFkyc8mz9/WIWmWO5h3fHY0I+t0I+NxFA5iPSOCJBBAm/qtnQyIx8OQ3Ek5KiJtpd3e2+4rino+Pz+PURBwv7h9U1C55GZUccnsH23F4NCIxJb1yclLiDzg3Jc0BkEdksTryoYfyEMasDvss315XZhYLDlQq/a080oMvk1cAhCgaL3weJ65HU0U7oGzebPWF3f0ABL8WP562BhtVd4VKcsgViZ29zXA5CpFf/KLFzCMAweZ2X9Eixdr0Y0OC7Z2NlXg6IwkicBCIzuLdA614dyIu2SHzUXtqc6AUhB1JwfT69I6CSYy3CMj7aSSclJyHRCCcLXUttoTDoxHRX6iFuoiz93iMnU2t1eiA+gSRn7NDqUBE8XQGz+/qw9Ebc0V2W37i52hxWCUbuHINjbepKmlCPJGCMyfwkDq2pBEpACSaNNslyB1CTG7ehJWUXZ20KdCAG/kzQHjfBtMq2aQcnV2EJ0YktnNeJrmWzxaywpk/U+LCy/cVT9xOB9DT7MXuTW3YvakN3zszJemTTE4wk3Q1efDC7n7FCB01CCYyBdt1Xl5+18cc0Lv6W3FrXrro8KYGIe/DBzcXTaXkuFOrD3vvrybmkaGCQjG9sCLRStWci4KzUog5Jyg4xFkorlCHKfH7zE9B8wrJwxuFEEpm/pITpBGSHO7248H+VtVFi50H9fZYBA/0tRQRvdshEC1P1M/v6lPMkQEKviSlOa7XKSofv4A0Y59AahIDINll8f10hEvuA6S8kdGZVGWUB7VgWbJvbjZ2fR6/tVVzfCgNVrn24cgTEdvuP7+rTxKlI4fSJFbKLHUqnFutpAkxGZTKVXJssQQfhieHO/F6PiKHJaHw2lqL18XtHgqhj3Jcjy3lzymRtidbNJXIEVDXoORaPoup5hOeeJLef1cHenwURydTODwaxVtXZsUdRlHklKzPlDQ3tgCr1Y23XTN55bs+oOCYu6+3GadvFxYcdixz4R0En3+kX4yp1rttZxju9uM/PdqnKH+p3V7B6yQFiy4qMofRgmzHx+fzZ+gHDUWNKEGedEYgLDrsqAtW/nKW4uT4vGS3vG9bB/b0NWB0rhDpwxLa2K4glsrg1ES8KFZeaY7z3NDidUl2lXx9lcaIUsY+A79L4fvp+V19ktBcOfQu/kZ5UAuWJftIJIKmpibd39fa2pb63Yv7h8QtJSAcpcAGeTydEZ9fjaQk5+bw2ZY85PZjth0tVSdehm6yoGijV3NsAcJgfePitLjFXckWJ6HwWg5zMCkNSJqvM+9HYHVUI8dSGhSvtfLJSjnKH1lRmGQdZBFdTV4xFZ6V+ex9PXhx/xBeOx9CZDGDYKNbTLdXIyZ53Y6MRYp2IPLxoxarnqXAmdsFe7QDwkIlf0c8ncEX92zU7HMtqMmvtZDJtV0eeza1qYaT8ljOUsRSy4aiRpTA4t8ZcrSQiOckEP1dfHgn8xV8cU8/xsfHEfL6FfNV5KfEyhdTuYau5HOSL/pq47cUtzDzGN9PR2/MKebBsDOa9C7+RnlQC5Yl+0AgUPI7aiYDo+B/x5sb5NqCdOUnikTPypPbj/WGfLF3xePFF3fwfgcCFJ38KF8MAIgx52rbTqbZ84fDCd8p9iPIt/PyCSCEHpZOB5e3Dx+yyk+yQCCAnQ1O0Swh3x2xyKdreeJgpgulS5xbvC6xbk7JTkAqk5bNmQeFUG+HzD9QandpBEz+UmXyJMVMTezfm9ob4HY6sG8oiM3tPpyanJf89tEtAfzTjZhERgrhYDL+DKBy5hWLf2djUm7+3LetoEzwzl827gKBAOKxRc1FFVA3oyqNWd5UKl88tawDpbhF/ts9m9pwMbQgdZQ7CX7vkT6RS5R2F3Lo4UG9sCzZp1IpzYP7tbblpaBkl2NQIjH++3p3D0o2eT3v15Kf1e218yGcvBnHifF5SQy3fDFwOEjR9l8uHyC137MdDn/sg9rkkC+Sr5ycVHRGl2of+WRndTl/ax6nI5SLsijQBbOHsk92bWzGO3lnuNzBKq/bg/2tJXcCTD62q4ulMkikM3gvf2IpBUQ7Ny9rObtLVkf571KpFIa7u0uWKSeaJ4c78YMLwiFl16JLcDsIpuJp0dfjJBATkja3++DzOAAQxFIrOJE3p2RpsYnQKIa7/arHILg4Uj83lVDMhk6lUqomFK0FUGuHWYrQy+0/pd/y45uXzQh3leJBI7As2S8tad+IZMTpwkNPQ8tJTP79UhOAn7jy7+rtaC35GakB0lhzuV2cN0epkTT7W+nfgD5fCFAcBqqUdi8HXw/5wjgSTuLPTsxKDrbjI4x4e6iTCAk37N1yeeV1C/ikUUxaGjgrg4+S2tJROHpDLms5u0v5mGDE100W0N3drVimmpmCmdN4cl3JUcmJmJQKkSub232Snerzu/pwejJu2s4EgHgMAms7doYOy3Dl3y+fC0tLSxgeUF7stEi5EkLX6r9SSprSvFL6nhHuKsWDRmBZsi8VX6qXhOTQ29B86r2RRaUUmet9v5r8pez2vMlJbo4qB3q1nXL7Q15vhnNTCWQ5oufzFHhbMPOH8KQif7+8bvIjMLT6X+6LAIQQRT4T2aisctLgy+eP3HA5CHp6iuPCtcwUDPKgA4l5Lz9mlHwM5Wq2atAaP1q7RqAwB5RIU4uUKyF0NVRiSZDDyFyx4+xR/pZLT0PzHcs7krQ8+QylyFzp/UqnQHaTBTy6c6ti/ZUSUsxqHyXomRxmvg/I23sdUwBFUeIQID3Slo/d5tuR7yelumlpcGrhfnoXCzUolT29sFzkrKQAMtnync180MHWoE9yjhLLQlEai+UQYSmolVlqLlYSY262HOVaEtTqpneu1EWcfUND6WvTyulQPQ0t1+Q+oehIUl7dSw1gJVJSCoV0EShqdUYGitr2X26XNwt6zA1GyvqXd7fifDQnOVKggGI7Pr+rKaX5akGvxltO26lp8fITJFdywhhQUkYEx2chll6Ps5k/m5/m6/HsfT2S5KbVRqmxrIcDVguV7lwB5YCSUjCzDSxL9j6fr+zfGrWtyaG07R/u1nZW8mWXImP+/fKMVBYKmYFyUlWp+mvJLo8b15vFVy4qdaL/zUgcK1nhSkYAkjNk5JnC5ZjK1FBNjZcvWx4m2dUknOTJ7kbYvcGr8U4q/i87hVJLVi0zF5/cVIl5ohxotWslHGA2jO5c5fOw3LlgZhtYluxjsVhZXmgzbGtqHSsnATWTjhFi4MuUh0KWYwsu5S/g48b1ZvGVi0pI99xUITGNJV6xCJtSWpYZ/gOzbddKZcujkPiEN7eDYEdzE4CuojL4xQ60cB5PqXBXtQxys8wTZqNcDqgW9M5rpXlYbjub2QaWJftgUP9RpED5DlU1qDmF9CRnGH2PUijkUJvTcHl6/AW8405vFl+5KJd0R8JJTC8sSxyL/LG5aqF6DGaQdTVs10pl81FI8v6bSCunysvbVastSslkhnmiWuA5oFxzYDWhVieleVhuOxvlQS1YluwTiYTuzDEth6rZg5dNGD0mHXkd9UYHDHf7MTU1JS+iZFl6/AW8485sm70c5ZCuJBs0r7Uin7gk13y1FthqkrVW3Y0SEl9PeVLalialPFjtdjVah2ruYioF4wAzI2HMglad1MyA5bSzER4shYrInhCyE8ArAJoA3ADwLyil8fyzLwN4DsLNcr9PKf2ZkbKNHNqv5lCt5uAtRazyCBujg1VNfrVBxt5XSsurhAQrJTM9ZfJ96QDwiW1BsS+tbHKolJCUktI2NiZVv6/mDC+nDrVYGPWAzYFy+r3aOwGtOqkRezntbKXLS/4cwL+jlP6SEPKvAXwJwH8khAwD+CyAuwFsAPA2IWSIUprVKEsCI/Glag5VJZg1CEppV/yke3zQ+KFSWnH28rIA7eQUM1AN7UqpTL4vXQp9aVWTQ6ULEfs9n6jVs8VYjLWVF8NywOaAURPIauwE9OyizXinleLshwD8r/y/jwD4GYD/COBJAN+nlKYBXCeEXAXwEIATegs2El+qd4tk9iBQ61D5pJNfWq6HpNTkVxpkqzHJq/EOpTJZOCDLNZCbt6xqcqjU9q30e6Mx1la2v5cDJr/Rfl+N+bBaY9FKcfbvQSD2HwN4GkB//vM7AJzkvjeZ/0w3jIYc6VlJV0vzKSdjUw6fz6e4C1EbZNWe5NUgErUyWV9OT8svi7auyaHSya/0++lpdTNONepgFHp2yZXspHkOMNLvq7XoGalTue2wqqGXhJC3ASjtJb4C4F8D+BYh5D8C+DsAhg1M09PTeO655+ByuZDNZvHUU0/hwIEDiMVicDgccDqdiMfj6OzsRDQaBaUUnZ2dCIfDouNiYUE4Q2RmZgaEELS3t2NmZgYtLS3IZrNIJpPo6elBN1mAiwgx7E4C3NksvD+VSqGnpwehUAgejwfNzc2IRCIIBAJIpVJYWloSnzc0NMDn8yEWiyEYDCKRSGB5eVl87vP5sMHrwR/c34TJtBeb/Fn4l2bR09MD/1IafqcLsdhySZnOjM/i5Yu3sZKlcDmm8KePbUQQC2hpaUGXM4vdgTR62gIYHx9Hm9uNrz7ah3evz+KB/jZ0kCTGx2fLkmmWNuL0xBy2thB8aHu/pkzptAuhUAh+v7+sftrU3o4/uL8JNxad2B70SMp0u91wOBwYHx9Ha2srlpeXTe8nj8eD+fl5dHR0YH5+HisrK+JzIzLFHM04eW0G29s9eHJIkC8WW5aMPSZTa2srZmdnFWUKeDz41NZmRCKziMeF346Pjxsee/5AGh1Nzbh9+3bZMpWaT3POFvz7wzeQyQFuZwj/fMiHjNOLoTYn+nw59PT04NjlCXzzdP6qUMcUvvqRO7DRT3X3UywWA4Cifro8s4jbKz5scKdw7x1tRTI1U4o/fXwjTrw/jfs2tKDbtYTx8VldHKGnn4yOvV/fjOC/nIwim2+Hf7+7A3f3NOsaew0NDZicnFTsJ6MgVOmG43IKImQIwF9TSh/KO2dBKf2z/LOfAfgqpbTIjHPixAm6ffv2ovLGx8dN277wsGIIlxK++4vL+Mn7KdFh99sP9Bo+gdCorPIjnPWYuarZntUaA2ainDbTC6vJLz8v6H+enpIcq8xi/VkbvHo2JH6nnDH883NXEaZNkrGlp72tNsfl7cAHHZSqn9YYOHPmzOnHHnvsg3rrUWk0TheldJoQ4gDwHyBE5gCClv+3hJBvQnDQDgJ410jZHR3y9HhzUE0zgJmD7KHNHXjzevGxv0bqYtQ/YdTMVW1HWLXGgJlQc5ibMQ6sJL/WeUF8JjA/bioxp4yEk/jmmQVkcgnJ2Co1Rq0YpilPnFQ63kINZo6BSm32zxJCDuT//UMAfwkAlNL3CCGvARiBYDU5YCQSBwDm5+fh99d+VdYLswdZr3u5IvtrOf4Jo5Oz2j6QtTAGlLKqzRoHVpJf67wgtfsIKvEhnJtKYCVbfAxEqTFqxYgkvh20jh5XgpljoCKyp5S+BOAllWdfA/C1csteWVkp96c1gdmDbGVlBcMbyt+FlKNVGZ2c1XaErYUxIG8zM8eBleQvdV6Q2kU95e6kWbZ3lsLQAmLViCTWDkaPHjdzDJhmsy8Xajb7dDoNr9dbgxqVB7Ntt3rlL3XwWbVtl9V8x1obA4C548Bq8q+2LfzcZAwjs2nD7zNSz1rY9w9emsXRG3Mqp7lKoTUGjNrsLUv2VnNO6YGZA0eP/NV0DloBa3EMAOaNg3oPUqh2/9di/hh9p2UctNWEVWyVRmCm81eP/Fa0T5qJtTgGAPPGQTXkt6IDUw3V7v9azB+j7zSzDRymlWQynE5nrauwamA3zY+EC0k0euQX7JMEDqJ9vO1aRT2NASVUQ3616CErotr9X4v5Y/SdZraBZTX7eDyOQCBQ62pUHWqalh75rXx8gBmolzGghmrIb1UHphKq3f+1mD9G32lmG1iW7Ds7O2tdhVWB2rZOr/xWPT7ADKz1MVCpbbwa8q8lBWE1+r8W88fIO81sA8uSfTQaRWNjY62rUXWoaVr1Ir8W1nIbmGEbr5b8q3GmixlYy/1vFsxsA8uSfa2jhFYLappWvcivhbXcBmY4/2otf62dubWW3wowsw0sS/ZrfQtvBEqaVj3Jr4a13AZm2MZrLX+to71qLb8VYGYbWDYaJxwO17oKNUW9yw+s7TZgO7bffqC3bI241vJXO1pFKQqNR63ltwLMbAPLavZm3bu4VrGW5K+WXXcttYESKnX+1Vr+ajpz9ZiIaiG/1RLOzGwDy5K9jbWBWtt1bVQX1YpWqbWJSAnrfSxb1oyzsLBQ6yroRqntaDlYK/JXM0lnrbRBtbCe5ddjIlpt+a2YcGZmG1hWs+/u7q51FXShWtrAWpG/mkk6a6UNqoX1LL8eE9Fqy2/FhDMz28CyZD8zM4P+/v7SX1xFKNnz9G5HjdoCrSi/Eqpp110rbVAtrHf5S5mIVlt+KyacmdkGliV7QkitqyCBmgavRxsoR/u3mvxaqJZddy21QTVgy7/68lstI93MNrAs2be3t9e6ChKoafB6tIFynFFWk78WqPc2sOWvb/kBc9vAsg7amZmZWldBAi2H0nC3H8/e16NK4OXEK1tN/lqg3tvAlr++5QfMbQPLavYtLS21roIEldjzyvmt1eSvBeq9DWz561t+wNw2sCzZZ7OG7idfFVRizzP6WyvKv9qo9zaw5a9v+QFz28CyZpxk0ryY9bWIepcfsNvAlr++5QfMbQPLkn1PT0+tq1BT1Lv8gN0Gtvz1LT9gbhtYluxDoVCtq1BT1Lv8gN0Gtvz1LT9gbhtYlux//OMf17oKNUW9yw/YbWDL/+NaV6HmMLMNLEv2P/zhD2tdhZqi3uUH7Daw5a9v+QFz28CyZJ/JZGpdhZqi3uUH7Daw5a9v+QFz24DU+uqvf/iHf5gBMC7/PBqNdrS3t8/WoEqWQL3LD9htYMtf3/IDJdtg4LHHHtN9lVXNyd6GDRs2bFQfljXj2LBhw4YN82CTvQ0bNmzUAVaV7AkhbYSQHxBCLhNCLhFCHiGEPE0IeY8QkiOEfFD2/S8TQq4SQq4QQj7Ofb4v/9lVQsj/uZoyVAIz5CeE9BNC/pEQMpL/3RdrI015MGsM5J85CSG/JoT8dHWlKB8mzoGiclZfmvJgYhv8Qf43FwkhrxJCGlZfGuMwIj8hJJif7wuEkJdl5TxACLmQb5tvkVLnIVNKV+0/AP8TwO/m/+0B0AbgLgDbAPwCwAe57w4DOAfAC2AzgPcBOPP/vQ9gS76McwCGV1OOGsvfC+D+/HeaAYyuFfnNagPu+R8C+FsAP621XKstv1I5tZZtNdsAwB0ArgPw5b/3GoDfqbVsVZDfD2APgOcBvCwr510AuwAQAG8B+ITWe1ftIDRCSCuA3wDwOwBAKV0GsAxgLv9c/pMnAXyfUpoGcJ0QchXAQ/lnVyml1/K/+37+uyPVlaAymCU/pfQEgKl8GQlCyCUIA9/S8gOmjoEThJA+APsBfA0C6VseZslPCBlRKcfyMHEM3IRwkKOPELICoBHA7epLUBmMyk8pTQI4SgjZKiunF0ALpfRk/u+/AvC/QSB9RaymGWczgBkAf5nfev85IUTrGMg7AExwf0/mP1P73OowS34RhJBNAD4A4B2T61otmNkG/zeAPwaQq0ZFqwSz5DdajpVgShtQSm8B+G8QSH8KwDyl9HC1Km0izOq7OyC0BUNJHlxNsncBuB/AdymlHwCQBLBm7O0mwFT5CSFNAN4A8G8ppXFzqlh1mNIGhJBPApimlJ42uX7VhlljYC3PJbPGQACC1r8ZwAYAfkLIvzSzolVCzfpuNcl+EsAkpZRpoT+AILQabgHgb9rty3+m9rnVYZb8IIS4IRD931BK11JOuVlt8CEAv0kIuQHg+wAeJYT8tfnVNR1myW+0HCvBrDZ4HMB1SukMpXQFwA8B7K5Cfc2GWX13C0JbMJTkwVUje0ppCMAEIWRb/qPHoG1n/jsAnyWEeAkhmwEMQnBInAIwSAjZTAjxAPhs/ruWhlny5z3ufwHgEqX0m1WttMkwqw0opV+mlPZRSjdB6P+fU0otr9WZKL/RciwDE3ngJoBdhJDG/Jx4DMClKlbdFJjVd5TSKQBxQsiuvPz/CsBPSv1oNb3Q9wH4FYDzAH4MIADg0xBWuzSAMICfcd//CgTv+xVwnmYAT0CIQnkfwFdq7V1fTfkheOZpvoyz+f+eqLVsqz0GuOcfwdqKxjFrDhSVU2vZatAG/wnAZQAXAXwPgLfWslVJ/hsAogAW8t8Zzn/+wbzs7wN4GfkTEdT+s49LsGHDho06gJ1Ba8OGDRt1AJvsbdiwYaMOYJO9DRs2bNQBbLK3YcOGjTqATfY2bNiwUQewyd6GDRs26gA22duwYcNGHcAmexs2bNioA/z/F4UAe1xdvFoAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(mdbCH4.nu_lines,mdbCH4.logsij0,\".\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "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.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}