documents/tutorials/select_quantum_states.ipynb
{
"cells": [
{
"cell_type": "markdown",
"id": "6492248e",
"metadata": {},
"source": [
"# Filtering Lines by Quantum Number"
]
},
{
"cell_type": "markdown",
"id": "20126be0",
"metadata": {},
"source": [
"In this tutorial, we filter lines by the difference of vibrational quantum number between lower and upper states, $\\Delta \\nu$."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "2f4f291e",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"xsmode = premodit\n",
"xsmode assumes ESLOG in wavenumber space: mode=premodit\n"
]
}
],
"source": [
"from exojax.utils.grids import wavenumber_grid\n",
"nu_grid, wav, resolution = wavenumber_grid(24000.0, 26000,80000, unit=\"AA\", xsmode=\"premodit\")"
]
},
{
"cell_type": "markdown",
"id": "050608a0",
"metadata": {},
"source": [
"Normally, \"mdb\" performs the activation automatically. \n",
"We would like to avoid this automatic activation and \n",
"perform activation after the user defines the mask. \n",
"To do this, call mdb with activation=False.\n",
"We also would use the quantum number available in this database.\n",
"So, use optional_quantum_states=True."
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "c090f2fd",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"HITRAN exact name= (12C)(16O)\n",
"Background atmosphere: H2\n",
"Reading CO/12C-16O/Li2015/12C-16O__Li2015.trans.bz2\n",
"DataFrame (self.df) available.\n"
]
}
],
"source": [
"from exojax.spec import api\n",
"mdb = api.MdbExomol(\"CO/12C-16O/Li2015/\", nu_grid, optional_quantum_states=True, activation=False)"
]
},
{
"cell_type": "markdown",
"id": "1284c843",
"metadata": {},
"source": [
"When called activation turn off, the DataFrame becomes automatically available. Check mdb.df."
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "0574092e",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"# i_upper i_lower A nu_lines gup jlower jupper elower v_l v_u kp_l kp_u Sij0\n",
"0 84 42 1.155e-06 2.405586 3 0 1 66960.7124 41 41 e e 3.811968898414225e-164\n",
"1 83 41 1.161e-06 2.441775 3 0 1 65819.903 40 40 e e 9.663028103692631e-162\n",
"2 82 40 1.162e-06 2.477774 3 0 1 64654.9206 39 39 e e 2.7438392479197905e-159\n",
"3 81 39 1.159e-06 2.513606 3 0 1 63465.8042 38 38 e e 8.73322833971394e-157\n",
"4 80 38 1.152e-06 2.549292 3 0 1 62252.5793 37 37 e e 3.115220404216648e-154\n",
"... ... ... ... ... ... ... ... ... ... ... ... ... ...\n",
"125,491 306 253 7.164e-10 22147.135424 15 6 7 80.7354 0 11 e e 1.8282485593637477e-31\n",
"125,492 474 421 9.852e-10 22147.86595 23 10 11 211.4041 0 11 e e 2.0425455665383687e-31\n",
"125,493 348 295 7.72e-10 22147.897299 17 7 8 107.6424 0 11 e e 1.9589545250222689e-31\n",
"125,494 432 379 9.056e-10 22148.262711 21 9 10 172.978 0 11 e e 2.0662209116961706e-31\n",
"125,495 390 337 8.348e-10 22148.273111 19 8 9 138.3903 0 11 e e 2.0387827253771594e-31\n"
]
}
],
"source": [
"print(mdb.df)"
]
},
{
"cell_type": "markdown",
"id": "87188d43",
"metadata": {},
"source": [
"So, we can use the above information for masking the molecular lines. DataFrame is not limited by nu_grids. So, all of the lines are visible (but using lazy I/O of vaex). Let's plot the line strength with different $\\Delta \\nu$."
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "2f7501f3",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAACUoklEQVR4nOy9d3xk1Xn//z7n3DJVM+pdq+2dZWEBG4wNNu7ENXGJU+3E3yR24vz8LSlO7DiOg7tj3DFggwsY29hgWHovS1lgYRtsY5fVNmnVpZE0M/ee3x93RqtVHU2TFs3n9QJJZ+4995F25j73aZ+P0FpTQgkllFBCCQByrg0ooYQSSihh/qDkFEoooYQSShhFySmUUEIJJZQwipJTKKGEEkooYRQlp1BCCSWUUMIojLk2IFdUVVXp1tbWuTajhBJKKOGMwjPPPHNSa109fv2Mdwqtra1s3bp1rs0ooYQSSjijIIQ4NNl6KX1UQgkllFDCKEpOoYQSSiihhFGUnEIJJZRQQgmjKDmFEkoooYQSRlFyCiWUUEIJJYyi5BRKKKGEEkoYxRnfklrCLPDyY7D/Hug8BE4cVr0TzvnjubaqhBJKmEcoOYWFgoevhAevADd2am3PbfDK0/Ceb86dXSWUUMK8Qil9tBDw+A/g/s+f7hDS2Hat5zBe5djXvY+7X76bfd375tqUEkqY11jQkcLdO45yw1OHkAg+eH4Lb1nXMNcm5R8vPwYPfhFITn3MY/8DK98MtauLZVVR0DnUyeb9m7njwB0c6DuAdjRCCf501Z/yiU2fmGvzSihhXmLBOoVfbz3E//v1DtzUz/fv6eTTl/Xz95etnFO78or4IDz+PYj3TX/cSC/su+9V5RT2de/j8499nm2d205/wYEf7PwBQMkxlFDCJFiQTuHkwDDXPX5o1CEAaOA79+3j/CWVXLCkaq5Myy8OboFDWyaum+WQGAASqYUkdL1cTMsKis6hTr7y1FcmOoQxuH7n9bx16VtZVr6seIYVAcnOTrpvu43BJ59iuLMTMTxMYPVqqj/2UXzLl8+1eSWcAZhTpyCEeA/wTqAMuEZrfbcQ4hLgC8BO4Eat9YP5vu4rnTHPC4zDiIZfPvXyq8MpxAfhhZsg3jvxtUT3xLXjO2GgA0ITSBPPONyy5xaePv70tMfEiPHAoQdeVU4htn07hz//n7g7doyuaWDgpZcYePRRFv/42pJjKGFGZF1oFkJcK4RoF0LsGLf+NiHES0KIfUKIf55uD63177TWfw38DfDB9DIwAPiAtmztmw5SCII+g8boRJ/42IFu9hyfId1yJuDIc3D4KaatJYxCQv8r3jlnOLYe28rPd/+cZAa/98mRk0WwqDiIbd/O4X/519Mcwmk4eZKe224rrlFFxvaHDnP7D55n+0OH59qUMxq5RAo/Ab4DXJ9eEEIo4LvAm/Fu6E8LIW4FFHDFuPM/qrVuT33/b6nzAB7RWj8khKgFvgF8JAcbJ0V50GJjSzm2EsSGu+gZ9jyRAGJDSR7c086KurJ8X7Z4iA/C87+C3iOZnzPYA72vFMykYqBzqJOrt19N+0j7zAcDQRUssEXFQbytjRNf+zrugQPTHjfy8sHiGFREHNnTxa4tRzi6p4uBTgeAg9s6cZIuZ79p0Rxbd2Yia6egtX5YCNE6bvl8YJ/W+gCAEOJG4N1a6yuAy8fvIYQQwJeAO7TWz6b2Taf6uwF7smsLIT4OfBygpaVl1rZXh216Y8M8tK/r9N8JGExodrzSTSyeJGCdoSWXEy/CK48C8UleFEzMnbngxiE+UnjbCognjj7BC8deyOhYgcCQZ+i/7xi4sRid117L8HPPgetOe6wMBIpkVXGw7b5DPPbr/ZOmgrfe+zJrX9eEaaviG3aGI9+fikZgbOzWBlwwzfF/D1wGRIQQy7TWPxBCvA94KxDFi0QmQGt9FXAVwKZNmyZ5S0yPWDzJ1kOT5NXxbpk7jw+w53g/Z7eUz3br+YGjz0L3sSlenOLPNeu/4vxC51AnN+y+gX76Mz7HUlYBLSoOeu9/gJ7f/g7ikz0AnA4VenVERgC7Hj/Clt9N7hAARnpcDu3qYNnGuuIaViDEens4uP05Du14npG+QXzhMJXNTSzecC5VzfmNiOb0UUlrfSVw5bi1m4GbC3ndVzpjmGLyJwgX6BwY5sXjvWemUxjogL13gR6ca0uKinsP3svOzp0ZHy+RmMosoEWFR2z7do5/5SswNDTzwUIgo2fg+3kS7N92gsdv3o+bmP64Y/v6XhVO4di+PTz082s5snsn6NO94At338G7/++/5dUx5Hui+QjQPObnptTavIIUgnDApLZs4q+vgaG4y4meMzSV0vEitL84+/OEhP4jXj3iDENbfxs/2fWTjIrLY3EmRwrxtjaO/OcXoD2z+gmGgREOFdaoImD/thM8cN0uRgam/7cWEoLRM/ffN42Xtz3DHd/9Bkde3DXBIQD0nDjG7kcfzOs18+0UngaWCyEWCyEs4EPArXm+Rs5IF5pX1USJ2l7KaCySLjy6/ySHu868GyQdL0HfidmfpxSc2AHdk8q2zmvcf+h+jg1MlS6bHAKBbUxaspr3cGMxOr73PZLbt8+1KUXFoZ0d3PWTnYwMzZzrDJVbLN1YWwSrCoe9T23h7qu+Tffxo9PWi3o7juf1urm0pN4AbAFWCiHahBAf01ongU8CdwG7gZu01pnH9EVCddimsz/GQ/u66BmZmJaUAtr7h3mhrWcuzMseAx3w0p1MXmAGVBTkFMVG14XOfXDyzOIG6hzq5Nb9t+LgzOo8JRR6kievMwH9jz1G3513zbUZRcWJg73cd/0u9HBmxy9aX0Gkyl9YowqIl7c9w91Xf5eBzpOnOQSzbGJX5FBf5nW0TJBL99GHp1jfDGzO2qIioHNghMf3T15oBkho6BqIk3Rmd6OZc3S8CCf3TP260zP1a24SRmIe5cUZhM37N/NSz0uzP1GArc68SCHZ2UnHD6+C2CTkhlNBSjAUmGdmDSXWF+eBn7/IUG8Gn0cFOBCI+ApuV6Gw+7GHuOfq75OIDUx4bcWm1/LKC8/Rf/JU2rCvM7/zNmd+T14W2HOiH2uGTjVX6zOvIafjJeg9muXJSXAToM6cD1PnUCe/2/e7rM6VWp6RTqHr5t+SmGpAbSoIAUIi7DPn33YsnrvnIJ2HM0zlakCC7Tsz6wlbN9/CQz+9Zsp00c777wLz9N/NCoXzasOCpM4OWgY+26K2bGrP4KJ58Xg/sfjsipdzhoEO2H07p/iMsoA0wDxzbhwPH36Y/b37szpXijPvrR/bvp3OH/1o9ieaJgSDSPPMewbc9fgRtt2TObGBMsAXlFi+M+93fe6ezTz0sx/POG9C4vT08JKzNubVjjPvk5EH1Ef9XLi0klU1EaqDEwvNAMkkbHul1+NJOhMwU+ooEygTjDPDKXQOdfKrl34161pCGt7c5JkDNxaj/Vvfgr7sKFjM6mrsRa35NarAOLSzgwd+MYvUoPD+5w/ZRGvPjPdxGtsfvIeHrr8GnNk9hIYqqlj5movyasuCdArVYZvYSJyH9nXRMTj5/Iuj4WDnIEd7zhCncOwF6M82dZSCEGCeGcW5LUe3sLdrb9bnn2lF5t7772foscezO1kIguefh2/5mUP+d+JgL/dcuyMz6q40BEihqV8RpbIhvymVQmL3Yw9x77U/wJmEUaB6yTLC1TWTn6hM1rz+ja+u4bW5QufACI/u65z2GBcYjicZSZ4BxebRrqMZws6ZIATE578TjCVi3PXyXQyTYSvKJDiT0kfJzk5OfO97k/apzwgpIRDAXrfujKG5iPXFue/6XYwMzu73lRIMW1G3NHrG0Ftsf/Ae7vnR99HJyTsGOw5M3g0oDJNgNErjqvxroCxIp7DnRD8ygzKyq8WZ4RQ6XoTu6cnQMkIyAc78H9rb172PbSe2zbUZRUPHtT9GH5il3oUQnhOxbYRSSPvMKao/d89Buo9mMKU9DsoQBMt8hMrOjCLz9gfv4e4ffnvmGsJ4CIlQksrGZuqW5J8K/cx5XMojgpaBz7KIzPA5cbXLS2dCsbn7EPTnoS3NODNqCvccvIeeZM9cm1EUxLZvp+dnP5v6AGOG5zrTnLxoNk+xf9uJWRWWR5G6k5XV+M+I1NHuxx7i3h//cPYOAZCGQSAUZs3rLyUQiebdtgXpFOqjfl63vIp3nNXEuc1lmFN8aBJJzXPzvdgcH4S2p5l+qicAZkVm+83zmsK+7n3cujf3IfkzJX3Ufs21MDJN9Jac4oElnWpSChkMIa35//TceXSAe36S3ayrYYJhKZrXVhKY55HC3qe2cMcPv407PMln1jCndfTCMDEMRXXLYlrX57frKI0z45ORZ1SHbS5YXIHfVAQtxZLqyXOtCRde6YrN72LzQDt0zNSWGYNE1wzH4BHGzPOawoOvPEhXMoPfZQZY0pr3cwoDTz/N0N1357aJ62I1NGDV1efHqAIhMeJwz7XbcbIsEymlCJXbRCrnd6R7bN8e7v7hd9AjU/yiycTUjh4wTANfKMKq172+IFECLNCaQsAyaOuO8fMnDxF3pmaNdoGk68zv8Lv3CJzclYeNRGrIaf7+srFEjMfbsuzAmQS+eTyo58ZiHP/yl7NKL4xCSjAM/OvXYzU35c+4AuCxm1+is232dQQApPe2LasOzOvUUW/7cW6/8qsMD2TJGiAlUhnULV1WsCgBFqhTONw1yI8fO8jIuBryZPIzUggCM40/zyWOPg9DU1N2ZA4NQs3r9NGWo1t47mTukqEKhZCC4WwfS4uArlt/T2JHjrRhhoGwLPznbJzXnUcvPX2MnQ9lT+rmCyksy6R1fdW8TR3Fenv4/ZVfpffE7IgbgdGmATsQxBcMsXTTBQWLEmCBpo+87qOJT2DjHYIATKWIjfce8wXxQTj6DDm3ogKgQFmg87FX/hFLxPjl7l/Omh57MpjCxFDz93ko2dlJx3cn1ZeaHYIBZMCfXStrkXDiYC/3Xr87pz0EEl+ZSXnN/HR8ieFh7vjBtzixNwuOLvD+/YRAKoNITQ3Vixbn18BxWJBOIWgZmBmQg2mga3CEbYd75mcHUvchOD5LHpwp4UAi5tUV5iH2de9j54n8EO5ahoVP+uZt+qjjpz+Djty7yaTPjxGJztsic6wvzuarns+JmQUFUgkq6vzYoflJ+PfAz67l4LNP57SHHQhh+XzUL19FtKawwkHz8w5QYEQDFirDhydXw+5jfXT0z8P+/RM7vZpCvhAfhlg+UlH5x4OHH6SPqSke/PgxyeymUO2vpjXSim8ett8O791Lz49/nJe93JER7CVL8C3Pfy97PvDYzXuIdeX2sGWY3mc0EPURLJt/jQOP/PKnbL8nR9JoKRGGIlxVxZJzz8f0FfZ9uyCdwnAiyYmBzPLJ8YTmla4Y3YMza+AWFfFBeOVJSE6k181+z37omX8iO51Dndy+//ZpjxliiEQGj5wCQX+inyp/FY3hxnyZmDccv/Lb07egZgIhwLIQhoH/NRdgVFbmx7g8YtfjR9jzRIaqcdPA8puEohZ1i8vm3RTzk7f8mqdu/mXO+xg+H5Zl07hyDdVN+aW0mAxz7hSEEEEhxFYhxOWpn9cIIW4SQnxfCPGHhbhmVyyOaWQWKhiGl9Kbd+mjgXboyJEAbwLceSnH+fDhhzk6lCOvUwoSiSlMXtPwGir98+tmOfD00wzdc0/uG6UmmWUwiBGef904R/Z08cBPs8yvj0FaZr2yMUjtomjO++UT2x+8h0d/8ZPsNxAShERZNoFQmEhdfVGiBMhNee1aIUS7EGLHuPW3CSFeEkLsE0L8cwZb/RNw05if3w58W2v9t8CfZWvfdPCZijIrs1BTO94fKTDPnkKIdXlKaflG37F55xjuO3hf3vZycHBxiViRvO2ZL3T86Oq87SVsGxUpw4iW523PfCDWF+fua3dO3Qc+C9h+hT9osmjt/Oo62vvUFu6+5vu5bZKaRLcDQaSStKxZW5QoAXJrSf0J8B3g+vSCEEIB3wXeDLQBTwshbsXTQ7pi3PkfBTYAu4Cx7u+nwOeEEO8CCvIoVxW06RnKrLqlRer9q+dZ/37nAS9ayCek6e050A4Vhe1wyBQ7OnbwxLEn8rafgYGtbOabhFLfI48y/PDD+dlMSmQohH/lSozyaH72zBOeun0fsZ5cKsspSG+CORj1Ea2dP11HL297htu/+w2I55ZuFtpFKIUV8BGpqaNh5ZqiRAmQmxznw0KI1nHL5wP7tNYHAIQQNwLv1lpfAVw+fg8hxCVAEFgDDAkhNmut24FPpBzMzZNdWwjxceDjAC0tLbO2/eTgCIaZ2U1Ba0+Fbd6ljzp2k1vbxiRwXRjuhkSWQ0QFwO/3/54R8lfkd3AYcUYQ82gi0Y3FOPJf/5W/DYVADw9jtbRgVFXlb98ckes8wlioVFq3dnEZkar54RSO7dvD7Vd+BWc42yE85X0GLROhwR8KIrQkFK2gqkhRAuR/eK0RODzm5zbggqkO1lp/BkAI8RfASa21m3I0/4rnLL46xXlXAVcBbNq0adaPfMNxh6BpATP/4yU19I8kkPPnHuJRZR/ckv99lfI0mueJU+gc6uTBww/mdU8Tk4AKzKtIofOWW+FQngr8hgGBADIQwFq2bN4MrZ042Mu91+U2jwB4w0MCfCGLQNiipiU8LwrMJw8f4tZvfomRwRxSr9oFKZAapJIYlk24uoo1b3hjQYfVxmNeTPBorX8y5vuDpKKAQqEhGiBoZf6rJ5Ka7sE8P5Xngvbd0JMpVbZJxhGFAJDevMI8wMOvPMzRWH4KzGlIKTGUQcCYHzdLNxbj5PdzzD+PRYo3R4XD86aekBhxuGu2gjlTIfVwJgWEoj5CFXPfVtzbfpzffe2LDJzMJZ2bojo3TXBd/GVRbH+QppVrqV+yIm+2ZoJ8dx8dAZrH/NyUWptXCNpqVi2m8aTL7uN98yeF1LHXixYywiycmcbjyzHnxw3znpfz0IkzBgKB3/BT6aucNzMKHT+5DtrzUBuSqY+yUqiKCnyr5k894ZFf7aa/PX8pQF9AYViKqpbQnKeOEsPD3Hbl1+g9ngfVQ6mQGhASKQxClZUsOvucotUS0si3U3gaWC6EWCyEsIAPAbnzHOcZh7tjIDOncxACOgbi82OALT4Ih5+ALLWJp4WrIdIMgQxptguItv42njnxTF73lEiCKkhDuAGfOfdOId7WRtcPf5ifzVw3dWPx2OHspcvmRT1h/7YT7H40Tw0RKXIyDZiWonHZ3CusPfCzazm+98Usz/bCHmFYXtrIMhFCEC6PYvv91C1ZXrSOo7HIpSX1BmALsFII0SaE+JjWOgl8ErgL2A3cpLXODzdBHjEcd/AZmc6/wsCIyysnBxmKzwMOpIF2ODkTVXaW0C4EqiA0hSZsEfHznT8nRn7TWAKBVJKlkaVU+uZ+RuH4t66cflBttoy1hgF+P8I0sVetnPN6Qu/JIe67Lr8ff2VCIGxR1RwiOMepo5ynlcWpr8owUcIjahRq7qIEyK376MNTrG8GcpzrLiwaogH8pso4seICJ/rnCaNm7xHoOViYvaWEyuVgBQuzf4boHOrk9gPTTzBnA41GCcWK8hUE5jhFFtu+ncHf/376g2ZLZOe6SL8fq6lpzvmOEiMOd1z1fG49C2naYsEoPbYvaGL6DCoawnNKa/H4b27IfVpZe8zE2nVBSbTrEIxEsey5ixJgHkw0zwWCtqJzlrQVI9MIXxQVXQdgOI/UFmOhNfijhdl7FnjiyBN0J/LPwWQIA0MY86Lz6Pj/fCt/m6XrCUIgLAuroR6ruXn6cwqMLbfspfOVHCM9Pear9nyD7VdUNoRoXDl3qaOtm29hy00/z22T9L+ZBCEE2vV+QaHknEYJKZMWHg53x3Bm2QoxHHcYTsyxY4gPQtszkMe+/dMgBCTnvh31/sP3F2RfnXrynuvOo9j27Yw89lj+NkwL8UgJpknwDW+YU76j/dtOsP3+PHWNpe77UoIvZCCkJFLjn7MC8+7HHuKh667JfSM3FQK5GsM0kEIQLCub8ygBFqhTODWnkDm6Yw4vnZiapbMoGGiHE/lQWZsCQuSFfiAXdA518vTR3GiGp4JGo4We80jhyH9+If+bSomIRFDBIDIDWvhCIe91BM3oXUqZXtdRXevckN8d3rWde3/0HbLXLxEgJMI0AY1QBihJMpEEITBM35xHCTBP5hSKjYZogDKfiWAo49tDUsOLx7KU0csXYl3Qnel8QjZQEJhbTqB7D95LdzL/qSOBwG/6MYXJ0BxGQ123bya5fXv+NxYCaduYtbUY5XPTPZYYcbjzRy/kb/ZRAtob9EWAYUiqm8NzUmA+efgQt33rq8SHcvnlNAiJdhyUbeMmkxiGiZQG0dpaAuEIjSuKw4Q6HRZkpBC0FSf7R2b9vHisd2RuZxX6jsJQf+H210lIzB1FeCwR46aXbpr5wCwgEJjKpCZYQ4V/7lpuT3zta/nfVEpvGt0wCJx99pzpMW+5ZS8nD+WRTDEVJUgJgZCJVIJQRfF1ExLDw9zx/W8R6+nKfhMhIdX9hhC4qSIzgGFZCKkoq65m5YWvm9MoARaoUzjcHWPEnf3Nvb1vhFc653Da9/guClZPAC99NNRTuP1nwJYjW9jTm286cA8GBlX+Kl5b/1oaQ3Ojo9B582/hWBYavTNBSo/awjCwls8NtUU+6whCcdqdyXVBGpJQuZ/GFcUtMCeGh7nr6m/Tvj/H96V2EaaJ67oo5U0tK1MhEATCZfhCYVZe9AYiBVZVywQLMn00HHcIWiazHQDrHkzQNThHA2wDHbDv7sJfZw7TRw+35YkldBIIIUg6SZaXL5+zdtT2b3yjYHtL00DOEbVF78kh7v1J/uoI2mG0h18pUKZEGZL65VEqG4qnD5EYHubOq65kz2NZvi+lTA0VypT2eaq4LLx6gpQKw7QwfTbRujrKqqrzaX7WWJCRQkM0QGXQnrVHHEom505XoeNF6CswY4jygy9a2GtMg53thZtzFAiUVHNWZD7+rSvhZO66y5NCCIy6egIbzio6tUW6jpDM1xiPPPVVGIALgZBNuMpPTUuoqFHCE7f8KnuHAKe6wrQLhoEATNvySgtaezUFy8IXLqOqeRHByPzgqlqQkULQVrT3j8yanyuedBlOzNFU81APDBay+0l6fMRzhMcOP8ZL/bmrcU0FIQSGmBsivOG9e+nOJ+ndWEgv6e4ODWE2NBad2mLrnQfyX0cAcL00kmErpIKySl9Ro4Tn7tmcFylNr5bgOQEtBG7SQZqG14pqWZg+H/XLVrDi/AvnvJaQxoKMFA53xxhKzJ71NOlojvXM0WRzchh0Ia/tpoqVc/PG/NELPyro/lJIwmZ4Tojwjn3t64XbXAhEMIgwDMxFLUWtJxzZ08Wzdxye+cBMIEF4nZrebIIBQnsF5kC5j2Uba4qmrrb3qS3cf/X38rPZmKl007QQSqGUgZACf7iMsqpqGletKSo19kxYkE5hOO4QmAV1dhpDCZd4Mtse5Rxx9AUKQoI3Fo4Dhr+w15gE+7r38ezJZwu2v0LhM3xUB6qLToQX276d4Yceyu+maU4kIcAwkMEAVlMTVl19fq8zDRIjDvfksY6Am0q74zkDoVPyAoagojZYtDZUTzktH07cm0lAKaRpIJRBMh5HKYXruvgCQQyfTaS6pqgCOplgQTqFhmiA6pAvq9xZ//ActKQOdMDeIhSZ7SDZD+Zkj837Nxc01y+RWNJiecXyohPhtf37Z/O/afrpM9WGqqLlBM8/r6itqA/euJPBrjxpjKSEcwSeIqx2Uw1VZSYgCESsorShHt61nd9f+RWc4TxF5FKgDIXQ2osSpESZJkoqympqqV20hNWvL66ATiZYkE4h25qC1nCoa7D4swpHtkFPnsL0KSFBWnNChvfQ4Tw/SY+BQGBIg4AZYF3luqJ2HvXcex/Oi9nSKmcAKRGBANK28W84u2ipo233HWLPljwWzVNNOdoB1/EKzF6UoIrWhtrbfpw7vn8liVyU02CUvtz0+TBsn5cqMkwcJ4lQkvjwkKeqZprULF4654Nqk2FBOoVsawoAvbFE8XUVeg+Bk8di3qTQYPnBLG76aF/3Pvb0FWY2AVLMqFJhSavonUfHvvSlwm2e0k2Qto0RjSKLVKQ8cbCXx36VR+r21L1eCFCpYMAwBKZfYQYUzavKC15gjvX28Ptvf5X+9jzNkAhJMplAOw6O46BdFyEEUios20eosopwZTWNK9fMm+LyWCxIp5BtTcHR0NYdK76uQqH0E06DPtVCV0TcsueWgl8jkUzg4ha186jr9s3Q1lbw66hgEKOpCekvzs3l7mvzTNGRmknQaSZUCa6jsX0GoTIfdUsiBY0S0tPKJ/bkofMtxcchBCjTBK0xDeXpJSgDqRSWz0egLELN4qVE58Gg2mSY05ZUIcTFwEdSdqzRWl8ohFgNfAqoAu7TWue9ly9dUzjSNbsUkgYGhovsEAY6YP8Dxb1mEXHnwTsLfg2lFCEzVLTOIzcW48QXCkB6Nx5KoV0Hs7KiKK2od137An3t+adBEQIwwUmAYQFaYNoGkWo/dqiw5H4P//J6Dj6XJwLGlPqd1hohJC4ubno+RoDt8xOsqqJx1dp51YI6Hrkor10rhGgXQuwYt/42IcRLQoh9Qoh/nm4PrfUjWuu/AW4Drkut7U6tfQC4KFv7pkO2NQWA4WSyuLMKHS/C4IkiXEimPpHFwx377+D48PGCXkOh8Bt+WspaitZ51PGT66Cnp7AXEcITeZcK31kbCl5P2HbfIfY9VZjhO61BjJnz8gVN/EGDcGVheY4e/80NbNucJ7VgZYBSCOVRVyTHKOr5/EEMwyRaW0cwVEbzmvXzrrg8Frmkj34CvG3sghBCAd8F3g6sAT4shFgjhFgvhLht3H9jNR//GPjFmH3eBdxOgRTccqkpDCccumJFrCn0H4fhApLgpSFNUMV1CtfvuL6g+wsElrQot8pZUb6iKJ1Hyc5Ougo1qAajHUeYJiocIrBxI1Z9YdMQh3Z25LeOkEZKUQ3tyQtIA6RUaCBY6WfxhuqCpY7yIpSThlKjNBZau54+gqGQ0is6Dw8NIg2F6zqEKqsJlJXl57oFQi5ynA8LIVrHLZ8P7NNaHwAQQtwIvFtrfQVw+WT7CCFagF6t9eidT2t9K3CrEOJ2xjiLfOFUTWH2T/zDcZd4ooi59542yFg4NAcIWVSn0Nbfxo6eHTMfmCMsaRH1RdlQvaEonUdHvvjfkOUDR8ZQCvx+jJpafKtWFjR11HtyiM1XFYDqG0anl6WV6jZCYIcUZRV+Vl9QR6SqME0P2x+8h4euy9OwpBAIpUCDkhJXg+vqVE5MoIRAGwbBaDk+X5Cyyup5Q2cxFfJdU2gExvZOtgEXzHDOx4Afp38QQlwCvA+wmSJSEEJ8HPg4QEtLy6yNbIgGqAjYvNI1e/rshOsykixS+ig+CIceL861lEoVyoqDX+7OA4XADNBokGApqyj1hIGnnya2uQjy5EIgQyGk349v9ZqCpo7uuX4HbiEC49Rsl3a8VLxlCTSCSKWfmtaygg2rHd61nbu/nycpVCEQykA7ricnLaQnfmcF0DhIwySZSGCZJlpAIBph2XkXzNtaQhpzzn2ktf7cuJ8fBB6c4ZyrgKsANm3aNOs+w2z1FMCLeIeKFSl0H4KOwvEBnQ7hNYkXCfcdvK9o1yoWjn35K4XbPNWCimGAlBihECocKmgr6qO/eYkTe/KUuhScIglVnnQHnBpWE0LgC5gIISir8BeklnDy8CF+9eX/zN+GKfZTw7LQjosyDVzHwXUTCKmwfDZKKqL1jURralj7xrfMC2rsmZBvp3AEGKsY3pRam1fIVk8BvHlfvyXya9BU6NwP/Z3FuZbrwHAP+ZPNmho7OnZweKjQw3genCI5ur5HHiW5o5DpsJQEmZTeFHN5OVZzc8FaUV96+hjP35PHj27KIaS/KhvcpHdfFQq09pyC4VNUt4TzXkvobT/OTf/17+jhbN/fKePT0bQAw/LU0yQCV0nQGmV4rafKUCTjCULRcmzbN28H1SZDvucUngaWCyEWCyEs4ENAnsr7+cMpPYXZw9UwUiz+o/4ToAeKc630mz5e6CE5+NWLvyr4NSBFb6EsKvwVBe88OvrlLxdu87R2tuuC6yIsC1wXs7EwrKiePsLuvO87dnI5/b2bBCk8jgvLb1BeG8h7G2qst4cbPvdPDOWinIb2OoxcNzVUIXCcpMdQ6zpIKdFopKHQKVU1KSWBSISy2pp5O6g2GXJpSb0B2AKsFEK0CSE+prVOAp8E7gJ2AzdprQtHkp8lGqIBwnZ2bzzHgWcP9hSH6qLnUOGvkYaTSCV3C09z8WjbowW/BpzqPloeLSznUc+996H37SvY/mh9iu/IdUflNwtVT7jzqucLw70oRlUp0y39kCK+M0wv+s536ijW28NP//X/Y7Arl4hbpgzXkJLTlEqihMQwFEJ5HE1SSpxEEiEl2nUJlEWQwqC8vnHeDqpNhly6jz48xfpmCtRKmi8EbUXnYHZDOC6w50Q/r3TGWFVfwNay+CC0FY45dAKkAqus4DQXjx1+jPZ4e0GvkUZABagP1HNe3XkF7Tw69h//UbC9Ae/uKYT3RKIU0rKwmpsLIqjzwA07OflKniVnU4R3aQcgDC9a0KkJZikFpq0IltvULcvfBHOst4ef/9v/YeBkR447uV5k4zBmFgG0gEQ8jlSevUJrhDTQSQdX6TMySoB5UGieCxzujuFkNbrmoX84WXhZzu5D0HWgsNcYC2WCVfg37s92/6zg1wAvSpBK4jf9Be08OvG97xdOUQ1O0WSnxHSwLFQwiL1kSd5TR7seP8KuhwowKJmisEg/cEvhtaE6qWjEChiEqwLUtZYRqcqP84719nDDf/wzfe15Go5UHoUFAFrjaoHAxTRMNAJlGiQTCTQuwrDxBQL4Q2VUtyw+o6IEWKBOYTjuEDQtYPZFJw0kXV14Wc7jO4o0yZyCMqGsseCRwgvHXijo/mNhCAMlCvfvlOzspOt7eRJjmQrptJHWowVme9UqghdckNfUUefRAR64vgCdbukowU39J726nE6CVAJlKoIRm+bl5SzZWJuXKCExPMxN//Xv9BzNB/dUOsTxwhzD9qX+STRag6s1pmXguC5Kebkxy/ZRXldPdetilmzcdEZFCbBAnUJDNIAts6+xm0riMwr8p+s6QMFFdUYhvE9qWT2EamY+PEvcsucW+iikpOgpGBj4TT91obqCFZmPfunLkCxwbSlVP0AI8PtRgQD+tWvzrp2w+QfP53W/UaS7jlL3VgEoQ5B0NMqQ+IMGls+gcWV53pTVfvv1L9L5yst52cuzXSLTn3c9prygJEIqhGGiUi2ppuUjUB6lYdVaVl90ybyms5gKC5IlNZeaAkAsnmC40DeDzmIwo46BEFC1oqCF5p/s+EnB9h4PU5r4pI8lkSUFKTIP793L4O9/n/d9T4OUp4rMUiJME1Vejn9DfrmO7vjR8/S1FygdagCuV7JSlvfsIQDbkmjAtAzK64J56zja/L1vcviF5/Ky12gtB4HreHSuWoJp+zEsT1rTsC2EECjTxE0ksQJ+qpsWsWTjpjPSIcACjRQOd8eQRvZtpSNJl66BAlIZdB2Cg48Vbv8JSLWFFDD3HkvE2NdfwA6dMRAITGkStQtHb3H4M/+W9z1PgxCjKSNvaE2hQiGs1ta8FpifuetlDjxTwFmYlOay64JpCJQhcF2NQOALGNhBrw01Hx1Ht337a7z06IM57zMWQnqFkHSwIzQkk3EMw0KiMW2b5MiIN8kcDGH5A1Q2t5xxdYSxWJBOYTjuELJsILsui5G4S7yQVBfHn4dYd+H2nwwFnsf72fbiFJjBm0+QUhaM3qLn3vtIvlDg2ojWqeImKadgIi0rrwXmI3u6eOK3eUqzTIZULUEY3reOozGk8IjitMAXsIjWBnPuOEoMD3PLN6/g0LZn8mO39NJ1XpAg0NpNpY9StNiuJhkfQUpJYiSOSJH42aEQlQ1NZ1y30XgsSKfQEA0giivCNTv0dYCb57bAaZGKEgoYKfx8d54YKTOEFLJgReZj/5lHqoTJkK536VRlVgiEaSLD4byljmJ9cW77wbac95kWCq+bM4n3ayhIjLhYPokdUFh+RWVDMKeOo1hvD7f8z5c4uiuP0+RuOlUkPZ5Iw0jVml0EAu14o9haa5x4HGVZBCLl+EOhMz5KgAXqFKSAEwN5EucuBIZy7avOAqYPAhUF2Xrrsa10JXOZJs0OVf6qvBeZT3zv+9BehDmLdKTgOKNtqKosnDeuo3uu20GyUM8dKZ8mtNfLT4r4jtRgtpPQ2P7co4RYbw93fP9b+XUIaZI77SINE5GatEvPIkjDwEk6oB106qFDCkm4spJITe0ZHyXAAi00H+6OEbSz/9VHHJcRp4CF5q4iTjJDqvJXuMG1n+0qXuoIwJAGlrRoDDfmtcic7Oyk6wc/yNt+UyItiyrEaOeRCgUxWxblhevoqdv307azJ+d9pkS6XJfqOEp1aqIdLzNj2ArDEjlFCbHeHn7/P1/Oj2qaEKMpIxDoNG2F6+A6SW/MwnVTsqEaQwlMnx9lKAzLJlgeJVReeUbOJEyGBRkpBC2DoG0jSGbFlOq60NFboELzQAccfqIwe08JyaiCegGwpW1LwfYeD4FAoQhaQTbWbMxrkfnY178B8fzLUQKnCstp6JRmtmkibBvh82O3Lsq5nrB/2wme/n0BHzpSDmCUETVFvitUWlNHYPsV0brso4TE8DA3f+0L+dFVTjkCz1Cd+tZzEIYyPJI7y8JJJJFKYlgWILB9fuIjwwTKogSjFZTX15+RMwmTYUE6hWjAo7rNtqzgaNh/sp9YPJkS68kjjmyDnqP53XMOccf+O4hlWdDPBhKJaZhU2pWU+/InZjK8dy8DN9+ct/0mQI97N6b74g0DFYkgA7lrJ3QeHeDOHxWIiiztCDwuOHBSvQsqJUyGQCmJaSkCZTY1zeGsooRYbw+bv/fN/DgEQJombjIJUiCFgU51GoFXQzAsC9d1kSqV+9ICpRRCKpQ0CJVHidTU0rLu7DO2BXU8FqRTGE4k6RjIrS+7azBOR/8Iiyrz/Cc8uQeYx/WOWeIH24qQbhkDB2eUHTWfeOVf/jWv+00JIUbpsUnJOaqAH2vJ0pxaURMjDrd997n8zkOmHQGnf9XJU8R32gVpKqQJblIQCFvYITMreuyThw9x69f/m+5j+aD09ox3Ex7TqdACF43QLq4GwzRRhkIaBsLx2E/dpINItaEKofFHomck4d1MWJBOIRZ3CPsVvSPZzyp0DSYYihegLbWnODoDp8FL9OZ9286hTg4MFJG/iVQ7qpB5LTJ33b4Zp6BaCWOQJr4zDJAGwraRoXDOqaPNVz3PQGeeU57jQ+3USIXAy3wZSiAM4d1kpUJZGitkZkWP3dt+nN9c8TkGOvPFM+UNBHpU2HgC0U7SE8xBIKVCmhZae7xGiaEYQimUYeO4SQzDR7S2hkAk+qooLo/Fgiw0BywDv8rNHzpu9g5lWhSTLjsNw4ZI/nmPiqWbMBamMLGVnbcisxuLcaLQLKhwis4CvK/JJCSTqFAIs7GRwMZzsk4dbbllb+EKy2PnW7ymnNFgR7saoQSmT6SopG0MQ86aHru3/Ti/+u/P5tEhpG33lOyElAjpFUFcrRFS4GrHK/MkEoD2OpJcB9fVGIZJWWUVgXDkVVNcHosF6RQAOmPZFwxdvLpC3jHQAe0FEDeZCVKBYYIVyuu2N+2+Ka/7ZQJb2YTNcN6KzCe++z3oz5Mk5XRI1xTS1BZKedGCUvhWrsya6+jIni6evaNA0efY9JHyHrYRKaK71BSzm9AIJEpJfGGLSJV/VgXmk4cP8av/+nd6j+W5zpbObyUdhFS4joNItUklUs0ErqNBCpx4AiG8tJHrOlg+P4ZhUVZT86opLo/FgnQK3bERr9c4B8TjBYgUjm6D3mP533cmCAuECfH8qbzt6NhBR6K48xYKhZKKSn9+iszxtjZ6rrkmD5bNgHTrqZQpLgXv8Vv4/ZhVVVkPrCVGHG7//rb82joeBt5dRHsPS6blTS47qfZTIV1cx8UftpBSEK0JZFxgPnn4EL+54nP0nsjjZ0KeynEJIZCm6UUK2otsHCeJoQyEBjeZ8H45Q6K192+kDIUvHCZUVfmqKi6PxZzWFIQQa4D/ADqB+7TWvxZCXAx8JGXbGq31hfm+rtZQFrLp7c6+2DwQj+efFO/4DuakyOzG805zcc0LRbiZTgIpJLWB2rzUEw7/+2fzYFEGSGslQIo5DjAMzJpq7FWrsi4w/+5bzxRWcntUeObU90md6q51ASm8G6+cfZRw8vAhbvrCvzLU25tno11AIJVCC+FxG7kOQnnymWlKbKREO0kEGssXBuEilYEvdGaqqc0GuchxXiuEaBdC7Bi3/jYhxEtCiH1CiH+eYZu3A9/WWv8t8GcAWutHtNZ/A9wGXJetfdOhIRogaOTWl18QUrwTc6Rcqh0IVuWVNvvxtsfztlemsISFT/nyUk8YePpp4luKMF+RriWMldsERCiEikSyLjDf97MdtB8ogr63BpJ4HEcpER2tUw/kUoOGwCyjhMO7tvPL//xMfh2CSN3qUkNoCK8hwUl6dQKZStspZSKlREmBHQgiDRONS1K7KMsiWltzRqqpzQa5RAo/Ab4DXJ9eEEIo4LvAm4E24GkhxK14zxJXjDv/o8BPgc8JId4FjP8U/zHwsRzsmxJBW9EZy+2GHk/q/JLiDXTA0TxR/s4KwmsszyNt9tZjW4s6mwDe0JqlLOqD9VzYcGHO9YTDn/1cniybAWmtBNc9VWg2DJRtIwwjq9mEbfcd4sVHC0zFMbaeAKNSm+k5MK/YLDEshT+SeZTw8rZnuPkb/w0jeaDyFin2X+16BqZ+VoaJSEUxiLSPkF6HlG3hJOK4jpuiy/bSSj5/gEhVzau2uDwWuWg0PyyEaB23fD6wT2t9AEAIcSPwbq31FcDlU2z1iZQzGZ0MEkK0AL1a60krfEKIjwMfB2hpaZm17Ye7Y1hGbpVid/ywUa5o3w39xdEuPh2p1EUeyfC+/cy387ZXphAIDGnQEGqgIdyQ0149994HLxeQPXQs0pqU6QJzuh3VNDAbGmedOtq/7QSP/arAWhzpqeV0WS0ltams1GiFkGiRjhLsjKOEl7c9w81f+QLki0JG6xS3Ruo97mqkFDhOEolCKoVMTzQD0jRIjsSRSmL5TJTh/ZeIj+ALBl7VxeWxyHehuREY2+rQllqbFEKIViHEVXjRxlfHvPQx4MdTnae1vkprvUlrvam6unrWRg7HHXxmbqIerptnp9D9CiSLo0o2AXlmSH2289m87ZUpFApLWSipiCVyi1KO/WuRBtXAixCSyVO0FqmuI6mMWdNk954c4s5rCpyCHCOtKdLsKKk6uZLCK8hKgRASw848Stj71BZu/sp/5s8hjEVqFkFKgTRMDNPEtC2UkkhlemI5ponWXmrJtC1cLZCGgWEZKMPAH4m+qovLYzGnhWat9UFST/zj1gsauwcsg66B3CpwTr4jha7iDnmdgs4rQ+pVz16Vl32ygaMduoe7GU5mX6w/esWXoK+IzjnVaTT6VUqkZWWlsLb5B9ugUNpP6XSRHE3Jo9MZr1TjlOtolJGqOruZRwl7n9rCrd/6yqmoKWdb010TKaOlStmsvODGdRGORgCGzyA5HMcVrpdSQuM6Gtv2IZWBdjS+8KuL8G4m5DtSOAI0j/m5KbU2rzCUSDKU4wOJC/llSu2dg0lmACTY5XkbXPvx9ikDvIJCILAMi+pAddadR/G2NnqvK0hvw9SQEkzzFK2Fz4dZV4e9dMmsUkf3XLedrrYCthqNpbFIkd0ZVnpY7VR+XqPRWmDYBsEK34xRwu7HHuLWr38Rknn0ZmkJ09EHN+HRVGiNdpJo18FJJnC1ixQSaShPI0HjTTRr7aWShkeQlk0wUvaqIrybCfl2Ck8Dy4UQi4UQFvAh4NY8XyNnaA3BHFkdJDA0kqdoIT4Ix3flZ6/ZwvJDeUtenMK+7n0MUISOl3FI1xNy1WQ++Kl/zK9h0yE9mwCpGQUB/gDStpHBAL516zJOHT1z18vs2VLgmRDl9SNI49Q4hSbFfiq8p2uPbtq7GQfCFsGIRWVDaMoo4bl7NrP5yq9O+lru9nrhizQMzyG4Gu26HgW2lClSWpf40JCXNlKGV1gWXmuqdhyEEoTKy191hHczIZeW1BuALcBKIUSbEOJjWusk8EngLmA3cJPWeo76LKdGRdDGlyPNhQD8Vp58avch6HklP3vNCqnqYM3qvLSjXrn1yjzYNHsoFD7DR6W/MmtN5s6bf4uzs4hv1XTxM92OmvR65Y3qasz6Bvxr1maUOjq0s6OwkppppLqLtHMqdeQmxkUJ2osSbL+JL2wSrQmwaH3VpFHC1s23cP/V38uvjVIi0t1cGoQUHjurFJ72gWkghUInkp4ugmmN8h9JIVDScwhSSaQU+MKRBZU2SiOX7qMPT7G+GdictUVFwHDCoSvH/JGSAjvHWYdRnNgxR0VmAwKVUL8hL+2oDx19KA82zR6Wsoj4IrSUtRD1R2d9vhuL0V5oic2xSNNi6zGTXlKhEwmkoVCVlRmJ6fSeHOK2H2wvsLEppNpNtetRvAiZ1qXRuA5IJXFdL0/vD1tEqgI0rSonUDaRrXbr5lt46LofFcDGsXT4LtqV3kS1KxBSeTUDpTFNHxqBwCsmu24SwzRIJLyvyUQCZfkIRiILKm2UxoJkST3eO4TMuaaVxxHgjhfzt9esoEGZeSky37LnFtzRHsXiwcTEZ/io9lWzumJ1VqmjV/7pn2G4CJPk6eggLaiTdgrgdRzZNsL2YTU2zJg6Sow4/PabTxeusDweY9/ueozushCprmaJMsAfNInU+KhpmVwv4fHf3MCWm/Ks151mO1VGKpzBS8chUKaJUiZuShvBGRnxuopMg2QiVVx2JVqDaVle+ki5hCoqFlzaKI0FyX1kG5JgILeWVIT2/ssHeudIVEdKKF+cF6dQbN2ENExp4pM+1lat5Y2L3jjr1NHA008zdM89BbJuHNK02K57qqaQmvQSPh9mczNWS0tGjKi3fv85BjsLKAl7mt1jvqbe8iJVY0CDz2+ibI8Izx+xCIRtGlaWT0gbPfjzH+ffIcCpv6fjpBwtIARKGWhXo9GYpoVSyqsvOA4aj95CKYkQAsMwUUIihMYXWphpozQWZKRQHw3gzeZnDwFoN0/RQqwnP/vMFtKEqmU51xM6hzppG2rLk1GZQyKxDZvmsmbevOjNVPpnHyW0FWNyWYiJ36cKnl7XQxAjGsVe1ELgNRfMyIj60E27Of5iAdONY6eV0x2zqXSR63oP4yI1pCYN0MIr9pt+g6qmMNWLyjDM058377vuR2zbfEsBbE0Nn8n05LJEKoVQEiVV6iXP+SaGRxDSqx2Ai1AyFbC5KNMkmYhjWD6C0YWZNkpjQTqF4YRD33BuT/kCcdpnPTeD5qKeoMDyQe26nOsJ33nmO3myaXZQKMJmmAsbLmRp+dJZn99582/RxZhcHstrlH7TpNcMAyMSwbdmDZF3vZvAWeunjRJ2PX6EHfcXmElXjPlPp2oHAFKgpEcYJ6T3OyjDwtUugaBFeV3QKzxrsHynbi13/vDb7Lz/rvzamJpQBo+iQmsXobxBM1wHKSVJJ+HVDaRCa08f1DAN4iNxhPZ+D+26GLaNNEyU1gs6bZTGgnQKXYMjBG0YziFYkDJPmbeBDq/7qNiQBoTqIFSb81a37C/AE+AMkEjCVpj6YD3n1Z2XVcdR+xe+UADLpkG6hpBWobEsRCCAWVtLYOPZMzqEI3u6eOD6/GgTTws3lRqSqdoBpOoHqY4jlYoSpPDKJIaBlAI7YFDREKRpVcVo6ui2b3+Nlx59sAA2uqN/S43wuo5cF8dJIsCjr9CgcUnEh1GmPaoIJwAlFVpoTMvG9Afxh8PYwSDhaPmCTRulsSCdQkXQJteaaMRvUObPgw7wK0/A4FzM90moWu61o+aAxw4/RqJo1c5TkEgCZoCWSHYdR4f/9TMwVEhe6SmQVuwzTTAUZnMzwYsuouzNb57WIcT64vy+0NoI4N0xjZTOjwBt4g0oO6CEHq1Caq1Rljdf4QuaVNYHWXJODXWtkcI7BEi18oJUXmHDxbNTCeFRSKVo7T32U69LSiJTaTCV+h1crECQYDRKqKKS+uUrMU2LlnUbFmTaKI0F6RQABnO8j/lNg+aKPEwBH32OiWK3RYBhwbLLIDR77qix+J+t/5Mfe2YBgcAQBhW+iqw6jmLbtzNw880zH5gvKOUVQUd1KjX4fMhQiODGjZS/770YldP/Dr/9n604xfBhKSpsobyZBFK1cKnATXccIbB8Jr6AQajcR6jCR21L2ahDSAwP89uv/ReHt28rnJ2Og0gPnOHd/F2hcR0XwzTQQmAIOUpyZ1gWiVTaKJ0bs/wBwpWVRGpqWXnh64lU1WD5/AvaIcACdQrHe4cwdW5yNj5TkhdOvL45UFpDek7BF8l5pxf7it9Oq1D4TB8bqjdk1XF06B//vwJZNgnGai+nf7ZtVCSC1dxE2VveMqND+NXXn6bnaIFaZsdTYKegnVQBOdU5qySjCmsIr/U0WO6jaVWUQJlN00ovZRTr7eEXn/1/9B4vQEedVKSGIkaLxEopNNrTRFBeiCOVgWFb3tyE9thPtasRQmBaFkIo/GVRQhUVo+mi6qZFC94ZpLEgnYJtSIJ+k/4cRHI6Bkbojo2wqDKHIm18EE7uy/78bKFsCNfn7BTmivzOFCbVvmouab5k1h1HJ773fThShHRdOipIzyW4rje0JgSqrAyrsZHwGy7Bt3zZtNvcc9122vcWUCN6sgebFPFdWhshHSUI5bWd2kELX9CgeVU5rWdVEyizMG1Fb/txbvzPf2Wgo0AU8K4zqlAnleHNSKQK9lI7XpFZSAyfn1Vv/QP85ZUIKVKOLEWPndJUkEqhTG8PYVrsKxZV+hzA5/PR1NSEmSEz9IJ0CuVBm+FEbvkjQwpiIzkWJgba54YITyqoWZtzPeHa7dfmyaDMoVBEfBHWVa2jPlw/q3OTnZ10XTkHVBzJZIrwziO/U9XVBC+8kMg73zFtHaEonEZpjIsYRIqOyU0rqinhKYUqRUV9iPIaP42ryolUeynUw7u2c/NX/4tkbLBwxqXqA8o0UvPIOkV0B1oqtJNAmRYr3/xOmpevIOTzeTxGWnvym1qjDAPDsvGXlSGlQkiZv6aReQitNZ2dnbS1tbF48eKMzlmQTmEokWQ4x4lmISBg5UhzEeuCwa7c9sgGyoblb8qpnrCvex+DFOIGMD18ykeZVcaKihWzriUc/Pt/KJBVkyAtmJN+mpUSlERVVGBEygief/60aaOicRqNdQYyNZRGKshRAqnBSWiUKfAFDYJRHw1LyyivD45OLO9+7KHCEdulbUzPIrgax3URjscV5bguMqW5nEQgXYdARRUhnw8hhBckpKI1KSSGZWH6fEipvMjiVQ4hBJWVlXR0ZP5w8er/q0yCNEtqLi2pAdPAZ+boFDoPAHmQHZwVUvqDOYrqfOWJr+THnFlAoYiYEVZWrOSChgtmVUvovPm3JJ4tsvhPOnUkBFgWmCZWayu+lSumpcXuPDrAbd8tEqdRevbLSJHdpZelGO2gVYbAMAWWbdC6vpplG2tGU0YF4zGCVIE+FY1LCUmPjU/ijSgIAYYySDqOZ7froFXKAaQiBCE83iOP5E55g2ykHMUCgZjlQNWCdAoVQTtn5qLyoIU/10ihY3eOVmQBocAO5ewUtrQXQdR+HAxhELACrKlcQ2NoSkG/CUh2dtL+7/9eQMumget67ac+GzNaTmDDBnyrVk7JbRTri/Obrz2Vc8v0bJCuGzhpCYJ0g44GZXjMocFym/K6AIvWVoymjJ685dc8+oufFM4oJ0XJms4eWSYgcB0HIcFNJnEQKOVVwIXy0kowJviRKdlNw8BMRQ+Wz/+qThnligXpFHpicXpyzHxUh2yqwzmKMnTNRT3BhGgrRDK/qY7Hj58vvpCOgUHQCrIsumzWUcLhf/9s/lS9MkU6QhDCm1oOl+HfsIHQRRfiW7FiylrCLd/aSo5qorNDyhE4I163UZoN1XXAMASBsI0vbBKp9NO4umI0ZVQQYrtReMNoWmtwtTeYRpqm22vplYYx+mDnaoF2kkhleBQXYyIFaSiUYeILhVGG8aqvIeQDC9IpHOsbwidhIIensZxZs+OD0HMwx02ygDKhcSOUL8p6i6teKH7XkSlNKuwKNtRsmFWUENu+neH77y+gZWMw1hGMmVpW5eWYtbX4162d3iF8+1m6jhSBrfU0m/EcwxiH4N1tNZbPIhAxaVxZwaJ1VUSq/Ji24u5rvs/2u28vnE1KjYr1eK2nGtM0MUwLx3VTYU0CaZg4roshBI4USOmxooq0uI6USGVg2jaGZZWcQYYo2l9JCLFECHGNEOLXY9beI4T4kRDil0KIt0x2TCFQX+bPmbeodzhJR38O9YA56TwSnsJa46as+Y72de9jwC2uuppCETSDLI8un3WUUNSZhLQzSPMaSenRWDQ2Erzwwmmnlu+5bjttO3uKZ2sak7SkSum1bvrDFpHqIIvWVVLVGMK0FZu/980COgQx2nIqlPQcguN6NykpvQKz1h4rtlRIy8J1HJxUq6oQgkBZJEXjbWIFgpg+H5Y/MGcO4Xe/+x1CCF58Mf/zPHfeeScrV65k2bJlfOlLX8rbvhn9pYQQ1woh2oUQO8atv00I8ZIQYp8Q4p+n20NrfUBr/bFxa7/TWv818DfAByc7phDwmYp4jnskc81G9B6BwRM5bjJbSG82oawh6x2++NgX82hPZrCkRdSOsrJq5ayihGNf+3pxZhJgYoSgFBgGKhwmeP75004tb7llb/FaT9PmGiC8FD2kKLBJme0PWZTXBimv95+mi3Dbt7/G7ofuK6BReDd3QDsuOulVvtNpJABX4NUZdIpqw1AYlolp2lg+P2XVNUil8IfLsP0B7EAwoy6jWDxJe/8wsXh+6chvuOEGNm3axA033JDXfR3H4ROf+AR33HEHu3bt4oYbbmDXrvxI+mbqPn8CvG3sghBCAd8F3g6sAT4shFgjhFgvhLht3H8zcTP/W2qvouB47xAqx0JeZdDMrabQ9TJFrSaClzoK1+ekx7y1c2seDZoZEomlLJZElvDW1rdmHCXE29roufrqAls3BukW1DRRm9bISASzro7QRRdN6RC2P3SYZ+8oUsQ4RhchXb8dZUDF60DyhTz6iqrmEI0rKli0vorEcD8/++z/LQyPkZSnDBMSHDfVRipSgxLeS4l4HCeRgITjdUSZBsJ1MUwTy/bhLyujdslSIlU1SMPEsCxUKoU0E2LxJFsPdrHraB9bD3blzTEMDAzw4IMPcvXVV486hd7eXmprT5FQnnvuufT29s5676eeeoply5axZMkSLMviQx/6ELfckh9iyoxqClrrh4UQreOWzwf2aa0PAAghbgTerbW+Arg8k32F1yv1JeAOrXXG/YJCiI8DHwdoaWnJ9LRR2IbEMiCW5fyaBJZUBQlYOZRkug9mf262MPwQbcnaKXx3a9H89igUijKrjPcsfw9N4el1BsZi/19/vIBWzQDDAKWwW1rwr1k9Zfvp/m0nePiGvcWxKV0rSFNV6BQ3X6qWIIRHJheM+KhqCrHsvFrqWiP0trdx03/9G0M93YWxK00QKLzihjQNlFIkEonUfAEoy5vENS2bZDyOMk0Sw8NgKCzDixCqFy+lee16TNPCmWX9YGAkiauhMmjTOTjCwEgyt892CrfccguXXXYZGzZsIBQK8cwzz3DuuecSi8VIJpMYhsGGDRt44YUXuPjii0fPu/jii+nvnzjF/rWvfY3LLrsMgCNHjtDc3Dz6WlNTE08++WTONkNuheZGYOwjThtwwVQHCyEqgS8CG4UQ/5JyHn8PXAZEhBDLgF9NcswEaK2vAq4C2LRp06wZiMqDNrnIKUgg5MvxTdNeBArk0yBAWV6BOUtRnR/vLH7XkaUsGoIN1AUzpzJuv/oamEvaAp8PGfDj37AB//p1k7afdh4d4M4f7CyeTelWUxilsDBNcJJeGsm0FKGoj/J6P3WLI9S1Rjh5eD83feHfSA4XqB1KpDxUiqROKIV2XbQyvIE000ADps/PSCyGk0yiEZ52gpQYhoHp81O1aAlnv/Wdo4R2s6WsCNkGUkDn4AhSeD/nAzfccAN//dd/DcAHPvABbrjhBs4991zq6uo4duwYzc3NvPjii9TVnf7efuSRR/Jy/WxRtO4jrXUnXu1g7NqVwHjegb+hwOiJxT3GxywhBZT5c5DzHOiA9h0zH5dXSPCFoeW1WRWZtx7bykiRB+1MTMp95bREWvCZmc1VxNva6Pza1wps2RRI1RKk309w03mEL3nDpN1GiRGH3/1PcdNwpyHlIHSKvsIwBKalqGwO0bi8nEVrq9i39VE2X/l1CpbiVJ4Yjtf2pEG76KQn2iOVgeu6SCFwHRftOhiGwrBt3GQS07Tx+YPYwRB2IEDjqtU5EdoFLINNrRUMjCQJ2UZeooSuri6efPJJfvOb3wCeU3jDG97AV7/6VRoaGjh69ChPPvkkVVVVLF++/LRzM4kUGhsbOXz41DN5W1sbjY3Zt5mPRS6//RGgeczPTam1eY+hRJKQDV1Z3uOUBDuXntTugx7FRTFh2lDWmLUe89ef/nqeDZoeEomtbGoDtbOixz70v/9PgS1LYWz72mi3kQC/H6uxgcDGs6dsP/3NN55muK/I9aSxUF4tQWoJhsQfMqhqDrPygnrqWiPsePgu7r/6e4W5dlqb+jQVOkAaSCVxcSExgnZcXJUmYNJolxRNtsYOBgmVVxKKllOzZCkrzr8wZ4bTgJUfZ5DGr3/9a97xjndg217dccmSJdTX1/PII4/Q0NDA5s2bueOOO9i8efOEczOJFM477zz27t3Lyy+/TGNjIzfeeCO/+MUv8mJ7Ln+Fp4HlQojFeM7gQ8Af58WqAqM8YDOSQ/uRA4w4OYQaQ70Qn31xKSfYEahYklU9IZaIsaO7uJGNLWzCZpizqs/KmB675977SD7/fBGs45QjSEMpCIWxlywh9Po3TNl++ptvbqXzUJGm09K8RukOo6T3vRTeIJgVMPCHLOygQePycupaIzx7528LN6UM3t/NqyR7hmiN0BLturgalDJwUxwWXveRRkjh0W64Lr5wiEBZhKZVq1my8TwiNXXzkvL6hhtu4Pnnn6e1tXV0rbOzkxtuuIGGhgZ+8YtfcP/991M1xWT7TDAMg+985zu89a1vxXEcPvrRj7J27dq82J6RUxBC3ABcAlQJIdqAz2mtrxFCfBK4C+8td63WuohJ0uyh0SSn4JHPBBIYGsmhKNFRfA0CDD9Ur8qqnvDDZ39YAIOmhkDgM320RFp486I3Z0SP7cZiHPv0p4tg3SSQEvx+zOpqQq973ZTtp7//7nMcf6mIetzpwnI6KJGpriMJtt/EFzSpaAhS1xphycYa7rnmysK2nKaRcgZKSFw3iZYKw/RSRq7roh0HlRpak8IbQAOFYZsEy8opr2tg8cbzqGppLbytWeKBBx6Y9vWvfz33yPsd73gH73jHO3LeZzwy7T768BTrm4GJ8c88h9YQsGAk2+FRDX4rh2GYrgPZn5sVDLDLoPm8rOoJ1714XQFsmhoGBmEzzIUNF7K0fGlG5+z987+AeK7TJ1nAMMDvR5WX41u9mrI3vWlSh/DADTt5ZXuBOnimQtoZCLzwNj1KgcAfsqhuCY92Gd3yjS9waFtx6hxCa5RhgpQoKUG7SKnQgCElydT8n1ISYRhYgQDhaAVWIEAwWkHtkmULWkO50FiQNBcByxjthMsGhsqxptBbAFWq6aBMsMNZpY7u2H8HDsXlDQqaQeqD9ZxXd15GaaOO63+Ku71IrKJjIaU3oBaNYi1aROCcjVjNE9tmn7nrZXY9VOxBxTFIpZA84juBL2Rh+xXVzWEq621+8blPc/Ll/QW4sBhtNR2bbtOuSzKZREmJNBSOK7xOWReclINwXRdpGvjCYULRCiJ19SUN5SJhQTqFnlicwRweKpVS2NnSZg90wNEi5b3TMP1Q0ZqVU/jm1m/m355pIJEEjSAtkRai/uiMxyc7Ozn53/9deMPGIj21bBiIYBCzro7g+edPWkfY/tDh4ugipJGeQ4DT06MSDEugjFMiOXawjx//339huFAzCOOcgdeCikdRIQWO66ATAF4bqjLTcpoKQwhMyyZYVk5lYzOrX39pSUO5SFiQTmEomcyNEE+42aePOl6EWIHkCqeCGYDadbOuJ7T1t3FsuLga0qYwifqjGXcc7f/4/yqCVWOQ5jZKsZ+qsjKslhbCr3/9hLRRUYfTRu1LfZWMFpYRXnQbivoIRmwalpYR632JW795FSQLlHKTKe4MrRFSobXrpbMk4CRxUUghPOU0x+Nacl1POMcK+BFaUNHYSKS6jtoly0oaykXEgnQK9WUB76ElS6egXUF3tvrOfcfwPq3FgoBAVVb1hM8/8vkC2TQ5BAKf4eP8uvMz6jjquP6nuDuL3NuQVlGzLI/srrYW//p1E9JGR/Z0FXc4LQ0H71M9pp5gWF7KKFhm0bq+ir6TT/LUzQWuE7kOSOUNpDkOKIkQ2iOxSzooQ5GMJxGOixASjUZIAyUE/mCYQFmERes3EoqWl9JFRcaCdArDCYdYjvflkWSWHqW3yKMcRgCqlmU1n/BExxMFMGhq2MKmtaw1o46jOUkbpZFKHZnNzQQvumhC2qjz6AC/+8a2ObDr1BdXeJPKhnF6yujQjt+y9/F7Cm+L8EQZtJuuR0k0qdoBniKaaZkoyyY5MozQGsO2sH1+yuubaFq9ltb1ZxOMlpccQpGxIJ3Csb4hfAIGs+wqNaSkLpLlG/VwfvhJMoYVhPoNs04dFVtIR6Go8lWxomJFRrWEfX/114U3aipYFjIUIrhx44T2096TQ9z0laeKb1MqTZTSoAGdFsmx8KdSRjsfvpqO/UWoZ6WH0pSng+yMjCBS4wnJpIOUCoEkkYx7DKiANAzC5ZX4y8qoX7aMFee/lkAkWnhbS5iABekU6sv8aSberBAJGEQD1uxPjA/CsReyv/CsISDcAK0XzTp19J1t3ymQTZPDr/zUhjKbXj7xve+jd8+FlKk3jSvLy7FbF1H2lrec5hBifXF+8V9bcIupkzOG+RTtpfLdpMciYZiSYNRH03KbbXd+jZ5jh4pjkwZhGmjHwUkkPdEcrRGmgdDeMBrCa9gwLBPT9mEHgkTrG2hctZoV519YcghziAUrRZQLi4+lsvyzdR+CgSL2qgsLKltnnTq6Y/8dxHNWnMgcPuHLeHo53tZG15Xj6bKKhEAAwmHsZcuIXH45vuXLTnv5l18qskOA04nu8ByCAAJBi2DEj+U7wsM3fa6wDmG8YlXqZ8OyUcpAWTbKUKiUEE5aUtP0+T0eo2CYusVLOOftl7P24je+qhxCIUV2PvrRj1JTU8O6devyuu+CdArH+obw56C8ZkgBIotQo3MfUEQBXssPVStmnTq64olJyWkLBsuwMp5e3v/RgmswTUSq00j4/Ri1NZS//31E3vrW0+oIN3xxC7GuIutAj4X2UkdKgeU3CFT6cJM7ePGhqyBeICJDw/BCEsZ8mNJCQ+4oHzdKCFyEp6ImFUKALxSkvLaBaH0ji88+l9e8/8PUL1s5d/WD+CD0n/C+5hGFEtkB+Iu/+AvuvPPOvO+7IJ1Crukj70knC6/SXWT5TX85LH7DrFJHW49tpTtZvGhGILCkxXm15804vdz22c/BK68UybIxML2pZbOmhsDZZ+Nfteo0h/DLLz9J1+Gh4ts1FsJLHdl+k/KGIEN9D3B0182FvabjAhppqNQMQrpd1xtQcxJxdNJB400nG6aJAOxAEF8wTHl9LStecyHnvuNdROZyQjk+CK88Acee977myTEUUmQH4PWvfz0VFdkRXE6HBVlTgNzSR1l3HhVVQ0GCLzrr1FGx2VBtYdMUamJj7cZp00YDTz9N/003FdGyFKQEZSCDQexlywief/5p+gg3ffUpTr6c36fL2cLTQ5DYAYNw1E9326/pbivChLeUqfkCL0yRhkd1LaREiBTDKR75nhYKF4EybYLlFVQ3t7LqwtdTv2zF3HcXjQx4LbShahg86f2cpYb5WBRSZKeQWJBO4VjfEDaQ7bNdWcDAb2Ux0XyimEyjCiINs55iLiYbqoFB1BdlcfniGaOEw5/6x+IYNR5SQiCA1dxM5J3vJHjeptEo4bffeoaO/QNzYxd4BHfKC/d9AQt/aIgT+3/JYGeRJqidJC54swiAi0SangRmMpHAEIbHMw/YPhslFb5whOrmVhpWrJofDgHADnlh1uDJFFtgKC/blkR2ziAETZW1QwBoKffPXp95oKO4RHi2H2rXzqqe8J+P/mcBDTodEknACLCyfCV/tubPpq0lHPyHT0FXkfUn0rz/Ph92XR2+FcuxF7ee5hCO7i4y/fkkkBL8YRvX2cehZ38DxeSpEhIhvcEzpPI6i1yNo70JZWWaCCFQhkkwWonls4nW1lPTuoQlGzfND4cAXlTQ8hovQrBDeYkSCi2yU0gsSKcwlMhN4GRFTXj2ghwdLxZXQyFQOet6wq/2/6qABp0OhSJiR/jgqg+yrHzZlMedvPGXDN19d9HsGoWUEA6jAn6spUvxn3vuaNro9999bm4cwhh9BGGAobzBtOHexxlof6i4tkivg0hrNyWS46LRaGGghEQLgdYa07IJlFcSjJRRVlVL3dJlLNm4af51GFnBvDiDNAotslNILEinoIXGIAeyiWzK892HKJi04XgIC2rWQiRzeb6f7/h5AQ2aCEMaM2ovx9va6PiP/yieUWkoBVJiVFbiP+ssou96F4Gz1iMDAe740fPFp8CGUw4BzyFIAXYAejt+T6JvDmY2pARXoywTJ5lESokU0iO5EwKlFFJJAuXl1C9dTnldPYvWb5i3ojj5RqFFdgA+/OEP8+CDD3Ly5Emampr4/Oc/z8c+lnt3XtGcghBiCfAZIKK1/sMx60HgIeA/gAHgC8BO4Eat9YOFsKXcZ+UUZA/Hs7i5Hy9mPUFAtGVWqaNvPlM8NlSBoNyeWXt5/8f+qmg2nQbDANvGKI8SOGv9qEO4/QfbOLityGmsNNIKaiJNvRTn5P4bgM7CX1spvHa9FO+T4yCEAJmSyNRet5GWEtP2oZSBZdsEoxEqmhZR1dzCigsW1kBaMUR2CtHmChk+8wohrhVCtAshdoxbf5sQ4iUhxD4hxD9Pt4fW+oDWejI39k9Auq1E4zkGH9CWiW3ZoGcoka3oGgA+M4tQ4ZVncrjiLGH5oWZNxuHwr3b/ipGc+rFmh7AZpincNO30cttnPweHijSBm4ZMNfoHAthLlmAvXkJg4zmjEcKcOYQ0Um2nyXgb3a9cTcEdQlqmTXOKHdZxQEiv1dRNKaMJz1sZpsL2+QlXVlDe0ER163IaV65ZcA7hTEemkcJPgO8A16cXhBAK+C7wZrwb+NNCiFvxpDnHTz99VGs9gS9aCPFmYBeeEwB4RGv9kBCiFvgG8JHMf5XMobMZPEtBASF7lgFWfBB6isWpL7wIoXZNxmd84+lvFNCe02FgUG6Xc179eVNOL89Z+2k6QqiswF6ymMBrLsBqbuKOHz3PgWeK8EQ+HuOkNJUB8aEtJPq2FOf6OnXxtKayMlK6Od5DUTqjpaXAUAZCKfyRCPVLVy64dNGrCZnKcT4shGgdt3w+sE9rfQBACHEj8G6t9RXA5Rle/xIgCKzB6xBNV126gSnbe4QQHwc+DtDS0pLhpU4hkINqmgTK/ObsTuo+BMM9WV9zVhAmLLkEalZldPjWY1sZ0MVrq7SlzdLoUt606E1Tdhwd/rtPFM2e0yAEqrwc/9p1RN71bgJnreeun++dG4cAp9FfS5VgsOO34BYsgJ6INLuekp7cs1I4rktattAFTMtCSoEdDBKKli/YdNGrCbnUFBqBsSO6bcAFUx0shKgEvghsFEL8i9b6Cq31Z1Kv/QVwEniPEOKtQBQvMpkUWuurgKsANm3aNOvH/p4ceLNdYMSZZUXi+A4gS/2F2UIZUH9Wxqmjf35o2qxfXqFQhOwQG6o30BiavAi+5yN/ApO04xUchoEIh7EaGghsPJvAWeu558Z9c+cQ0lDgJE8y3H4zXma1iNAuaQoL7Theq6lUuEKjDAMpJEIqwhUVhCqridTUzd/uohIyRtEKzVrrTuBvpnjtJ2N+LPBsPpQFcogUBPQOztKpHCoWXbbwppjtsoyO3npsKydGiqcdHDJCnFV5Fhc3Xzxp2ujY176O80wRay9pmCaEQliLF4/qI8xZhDBO/Gl46BkYLHK7aRop2VGpDNxUR5F2QUkDJRVCSPyRMMGKSprXrGXJxvNK6aJXAXJxCkeA5jE/N6XW5j3kWAKvWUJriARmmT7qKFbLoIBoM9SfndHRX3niK4U1ZwxsYVPpr2RtzdpJo4TY9u30XH110ewZhZRg25hVVYRe+1rK3/de7rm5bc4jBGHAUMftQBGoUdLayUJ40UFablRrvPDASxk5WmMYFr5gAMvnQ1k25XUNhCorWbzxPKpaWgtvawkFRy5O4WlguRBiMZ4z+BDwx3mxqsA43JU9v7EQYBuz7D7qO5r19WYFacKad0PFohkPjSVi7C5if7ulLFZWrOStrW+dNEo4NBeiOSlZTbOpCf/atZS96U1z6xBSbytH95LouB04XpzrplV5EEjDxE0mvZ+VQUonDZlKF0mpCESi1CxejhBQ07qYUHkl0bkktCshr8jIKQghbsArClcJIdqAz2mtrxFCfBK4C68p51qt9RyI0s4eVWWzfNIfA9MQ6NkwpA50FE+C0/BDsDqjQ//xnn8srC1jYGLSGG7kj1b8EU3hpgmv7/vLj0KWTJE5wbYRkQi+FSsIXnQhm2/v5tieOSS3cyHhHsbp/T1QTGGGlDqP8IrJ2tPA8ZyDkEjpeSupFIFoGf6yCFXNTSzb9BoMy8by+Uspo1cRMnrk1Vp/WGtdr7U2tdZNWutrUuubtdYrtNZLtdZfLKyp+UPUl4VqWgpKidlJcb7yBBRLsCZUDZXTE8sBdA51sqWjSG2NQNgOs7p8NXXhiU+TR6/4EoktxbNlFKaJLI9ity4i8s53cv+LlXPrECQMJ3fj9P6a4jqEFFzHSxElE2jHSWWQFIZpokyDULScSE0tNS1LqGxqpmXd2URq6koayjOgUCI7hw8f5tJLL2XNmjWsXbuWb33rW3nbe0HqKZzoG86+quDq2QnsHHw02yvNEga0vC6jVtR/vPcfC29OCiYmzaFm1levnzCo1nPvffRed13RbBmFlGCa2IsXE1i/nnufYE4dgibBcP+T0H8H5DRWOVuI0WJyWhxHSIlQyqOtMA1PECcQoqy6hprFS2heu57m1etfdemiWCLGyaGTxBL5FcEqlMiOYRh8/etfZ9euXTzxxBN897vfZdeuXXnZe0E6BSGy/+hJKekamEV76dEi8NoDmEFoPnfGVtTOoU62dW0rikkKRdSOck7tORMG1ZKdnRz75CeLYsdpMAwIhzEaGrDO3sSDXWdz9OXsW5RzhUuMkb57If5Y8S+ulFdPcN3Rr67joqRCKgPTMjFsH4FolMrGZs667O0sP/9CFq0/+1UVHcQSMZ5rf47dnbt5rv25vDmGQors1NfXc8455wAQDodZvXo1R47kJ029MAnxslFNSyHpaAKz0VLo3p/1tWYFXwgqFs942D8/ULy5BEtaLIosmlRmc+8HPlg0O06DEJj19cgNF/DQiZV0zWGTkeP2kui9FygynYeQpxQ0hfT6rLVGSI/ETgoIVZQTiESxfH7CldXULllGddOiV5UzSCOWjOFqlwp/Bd1D3cSSsWkFnzJFsUR2Dh48yHPPPccFF0w5JjYrLEinYFnZO4WQT2Jl2n000AGDxZgDMKBmdUasqE90PFEEe7woodJfycWNF08Q0Hn57z4BeXqqydwgCZaNqqzEWLaKJ4fW0NVXXBPGIpE8htN/N0UhtEtjdEJZgesipEA7pFKiEq/TCExfgEAkSnl9I/XLV2KaFi3rNrwqHQJAwAgghaR7qBshBAEjd4cAxRHZGRgY4P3vfz//8z//Q1lZZvNJM2FBOgWVw5xCmd+iIpRhofrAgxQlR2yYULViRlbUj2/+eOFtScGnfKwqXzVhUO3E977P8P33F82OU/AoLOSyVTxvns3J7vyoa2WDROIYzsBmoMgdVykuIyklrvZaUIVpIKVAAGYgiM/nw19WTnldA+sufTORqppXfXdRwAywsWajFyEYgbxECcUQ2UkkErz//e/nIx/5CO973/tytjmNBekU2vuy7waqL/PRUpGhGMeeIonDWGFY+qZp6wnF7DhSKGoCNayrWXfaoFrfI4/SdeWVRbHhdIMURCI4jUvYHnk9x7unVnkrNOLxw7iDdwHFClPGCDEASIV20yPTGu1otFDYgSCBcBmBSJSKhiYaV6551aaLJkPAzI8zSKPQIjtaaz72sY+xevVqPv3pT+fNbligTiEXmgtDziLKOFGk4bBQjafHPA3+/p6/n/K1i2su5pH2/Kk9BYwAy6PLTxtUS3Z2cuT//t+8XSNjCAGWhWxo4cXI6zjeOzcOQesEI0O7YeQRKCJNOVJ46SHpxcdaSXBdlDJwnSSm7cPy2YSqa6hqaGTFBRdT0dBYajXNEYUW2Xnsscf46U9/yvr16zn77LMB+O///m/e8Y535Gz7gnQKtsreKQgpGRhJZibH2V+ESWbpg6plYPqnPKStv43t3ZN3QV1SdwkVoQqYQGyeHUxMKv2VvGf5e04bVHvl//0T9PTk5yKzgc+HjEQ4VPUa2pJz00bpujHiA0+As52iaSinaSqEGK0laAHeZJpAGgplmpRVViGVSVV9I61nnUvzmnUlZ5AHFFpk53Wvex1aFyY1vSCdwlAy+w9myJKZ6Sl0HYLhGYqIi98IiTi05TDLYAeh4dxp6wl/e9ffTvnan5/153z56S9nf/1xiPgibKzeSGu0dXTtxPe+z8hjc9ByadvoSCUHV72Tfawv/vVJO4THUg6hiEjfMNJfhfdAo3FRholl+zAsm6qWxdQtXcaidSXtgxI8LEinEDCzjxSayoOZRQm7b5vhAB+87b/h8e/n5hRC9dB60ZT1hK3HtnJw8OCkr71n8XuIOTGO9OSrv9mgNlB72qBaz733zU0dwecjUV7P/pWXc5TMtCXyDcfpJTHwBLj5GSqaGan6QaqTaBTpojKe/oG/rIxgeSVNq9ax+sKLS86ghNOwIJ3CUDYayynYRoY1hf3Th4+c/5dQuxoCkaxtQdhQvggCFVMe8k8P/dPkpyL4wKoPsKNzR15EdgSCcl8559WdUlQb3ruXY//yLznvPWsYBo4/xCsr3zFnDiGZPEmy/0HgleJdNJ0y0u5pa8LwhtGUUgSiEcoqq2lZt5H1l15W0j0oYQIWpFPwWQK/hKEsfMNIMsM8Xue+qV9TZfCaVEonQzGcSWHZ0HDOlKmj+w/eT/vI5MWCi2ovYmn5Up5tfzbVnZ4bbGGzomIFH1z1QSr9lbixGIc/82/FF8wxDIhGObno9Rwic0nSfMJzCPcCRWLHFTLlCLTXnuw4XoFZevrJhmniCwaxAkGqmltpXLmaFeeXlNFKmBwL0ik0VQRZWVfGydgwbT2zbU/N8Abae2zq1879yCl66xymqylrgWWXTulYPvvoZyddl0g+etZHCZgBeodz75WXSMJ2mD9c/oejxeX2H11N8oUXct57doZI8PvpXvUGdhoXz3x8AeDNINwLdBTxqum6gfCI7aREpByCVJKyqhpqFy+loqGxJIRTwoxYkE5hRW2YC5dX0T+S4OGXTnCoO3PHUBHOgHa77VmmZro0TkUJ4Kl+ZQNhQfWyKVNHt+y5hV5n8hv+2dVns7Z6LZ1DnTx+9PHsrj8G5b5yXlP/GlZWrgS8eYTuH/4w531nDcukq+kcdtqvB6f4b+144jDuwN0UfSgNQCg8DQSFZZoIIbD8foLllVQ1L6J28bKSbnIJGWFBOoWqkI+PXryY9v4R6sv8fOXuPRmdJ4CwncE083M/m/q1JZecLoIj7IyuPQHSgPLWKVNHV2y5YspT/3TNnxIwA2w9sZXDfYenPC4T2NisLF/Jxc0XU+mrZHjvXo78y7+MirsXDX4/J1svYkftW3CL7BC0TjAysheGHgOKmC4TAhAIw/AG0rRGCYHWmrLK6pJucglZoWifHiHEEuAzQERr/YepNQl8ASgDtuKpuf0HHiHMfVrrXxfKnqqQj6qQj7Bt8KOH99CdAYW9AqxMZhxemYZf6DVTt4fOCuFaaJm86+jHz/+YQSangl5XsY4LGy8E4Ej/Efqc7CdrBYKwHeYPlvwBb2h6A1bfEIe/9GU4eTLrPWdvhIBAgK7mc9hV+xZcstfKyAanZhB2ArNgz80WMsVfpLXnDBwXnUwgDROExPIH8AVDVC1aTN2S5SxaX2o1LWF2yIjZTQhxrRCiXQixY9z624QQLwkh9gkhpqXf1Fof0Fp/bNzyu/G0nRNAG/B24Nta678F/izj3yIHNFcEObclwylXAf5MyPQ6pug4CdRB62szN246lC+Gxo0TljuHOrly29QtoP970/8enTI+muNwnU/62FizkQ21G/Al4OSPrma42PMIpklv0wZ21b2NJEEge1W92cJxeon3PQzOCxTFIcCp6WTTSg0vaaRhYto2wUg5wWgF5XX1NK5aw+rXvYGqltaSQ5hDFEpkZ3h4mPPPP58NGzawdu1aPve5z+Vt70z1FH4CvG3sghBCAd/Fu5GvAT4shFgjhFgvhLht3H9TTVatBB7XWn8a+Fvgp8CHhBBfBYrGR7CkJnNytJ7YDNz7XYeYMoWw7r25dRulIX2w7DJPaW0crtl2DUkmt/GNTW9kU/0mwJtyfvjow1mboFBUB6o5p/YcKn2VDG59hu4bb8zsZL8vlfrIEYZBT8NZ7Kh8G3Hy8HedBZLJkyQG7ge9GyhGqkyMCuGkCeyUYaJME9PnQxoGtYuXsfqi1/Oa93+ItRe/sZQumgXcWIxkRwdu7MwQ2bFtm/vvv5/nn3+ebdu2ceedd/LEE/lhQM5UjvNhoGvc8vnAvlQEEAduBN6ttd6utb583H9TkSi0Ad2p7x2tdbvW+hPAPwNT5iCEEB8XQmwVQmzt6Mi9y6M5Q4I7j3Z+hu6jHTdP/dpZfzRxzcqiphCug0WvmbDcOdTJDXumfvP9yeo/Gf1+18lddMXG/5NmjkpfJWsq13BBwwUYJ7po+9KXYDiDHFxjA4TCpyZts4VSdJUvZ1f95YxYZRQzQkgkjpHsvw/cgxSeBTftDFKX0hrtapRh4A8FsXw+IlU11C5ayrpLL2PjWy+nftnKUnQwC7ixGLFnn2Vo925izz6bN8dQSJEdIQShkPcwm0gkSCQSiHw8aJGb8lojMLZK2ZZamxRCiEohxA+AjUKI9ETTzcBbhRDfBh4WQrQKIa4Crge+OtVeWuurtNabtNabqqszE6qfDuubIvgz+EvYClqrZ3AguyayHgJQvnxyqczEbFtiJSy9ZNK9rt529ZRRwnnV57G2eu3ozycGT9CT6JnltT2YmNQEa3hd0+uoHfFz/MtfgQMHZj6xvByzsgo6c9QQMAy6qtfw0tIPMGyXU0yHEI8fxhm4CzhCUWjRBaM1BIQApRBKYtk2ZdW1NK89mxWvuYi1b3gjLWvWl5xBFnAHB9Gui1FRCa7GHcyPNOtkIjuRSGRUZAcYFdkZi4svvpizzz57wn/33nvvacc5jsPZZ59NTU0Nb37zm888kR2tdSfwN+PWYsD4OkPxSP9TWFFXxqLKAC92TP+E4LcNov4Znuy7p+hkOuuPJk8dDcwyr2+GoPV1E/aaLkqQSP5y/V+O1hI6hzp5pC17VtSwFeaPVvwRb6p+HQPfu4rBe+6Z+aRQiMB5m4g991xunUmGQW/5Evas+kOGiFCst7DWCRKJw7iDD3EquC0QREozOU1ZobXnHFyNYZlYPj++cJjW9Wez5uJLMSz7Va95UEjIYBAhJU5XF0iBDOYnFVlokR2lFNu2baOnp4f3vve97Nixg3Xr1uVsdy6fqCNA85ifm1JrZxwClsHSmtCMTsFxYTgxTU1hoANGeiZ/bcOHJl8fGcrMyDTCtRCqnbD85S1fxpmCgfPc6nM5t+7c0Z/3de9jf292MqF+6efcmnM5r/48Eg89Qfcvb5r5JKUIvP71uAP90Jl9ygrDYKBiCXvXfJAY2VEOZwPXjREfeh7iOymODoIXFQgh0a4G10FZNmbAxh8OU17fQFlVLcvOey2Rmrlhfn01QQYCBM45B3dwEBkMIgNnhshOGtFolEsvvZQ777xzzp3C08ByIcRiPGfwIeCPc7ZojlAenDkFoZQglpiGYfWlOydfr15/+mzCWCRnmT6KNHvSm2Ow9dhW7jh8x6SHKxR/t/HvThMQOdh7kM6ZGFwngS1slpcv57WNryXY1sPxr341ozqCsWwpvlWr6LrxhpyihIHyVvYsfS99xetBGNNyupui6CCkKSscB0wLaUoEJpbfTygSpaplMQ0rVhIqryRacgh5gwwE8uIM0ii0yE5HRwemaRKNRhkaGuKee+7hn/5pcp6z2SIjpyCEuAG4BKgSQrQBn9NaXyOE+CRwF14L/7Va6515sWoOUBaYub897DOoCE5z3AtTjFWsf//k6/FB6GvLwLoUjCCseueErqMrn5m6BfU9S98z2nEEXurovoP3TRlVTIfKQCUXN13MpTWvZeA/vgonMtCfrq4ictllDDzxBBw7PutrAiAEAxWL2LPkvfTYdRSrhuA4vSRiT4HzEpC9Wl/mSGsgSKRpYJgWpu0jXFGB6fPTtHoda19fShedCSi0yM6xY8f48z//cxzHwXVdPvCBD3D55ZfnxfaMnILW+sNTrG8GpqisnlmoDMzcBVQftqaX4mzbNvn6uin0U0+8CF2ZTVMDHiPq8jeftrSjYwfPdT436eFlsoy/2vBXp609d+I5Xuyefc+0X/pZWb6SNy16E+6t9zN4dwZSo2VhKj/2V7gnOxjZtSvrjqOhcC17F7+HHn89FGE4TesEycQxnKFnwD1KXh1CmskUxhDZpb4Hj8ROSkzbu/HXL1tOw/LV2MFQaSr5DEKhRXbOOussnntu8s99rliQNBeTYUlNEAOm6N3xUBG0p9dScHomOWnF1Kmj7v0wkGnRUsH6D07Y6+tPT/3met+q952mfhZLxHj06KOz7joSCCoDlZxffz6Rp/bSmckbOhCg/IMfJLBmNce+dGvWqmvxijpeWfpuuoPNQJaUILOA68ZIDO9Aj+zGG6zPM8Y6xrEU11qDlNh+P6btwx8KUbNkGWdd9jYiVTWlyKCEoqHkFFJY3xSlPCTpGJg65z1timnX7ZOvXzAdrYUAN8NCs68C6taetnT/wfvZ2rF10sODIshbFr3ltLW2/jZe7HxxVlTZAkG5Vc6ljZfyZnclnV/9V4jP8ORsGITe9jYqPvhBum68keShQxlfbyzi0RoOLvkDjgRXAIW/IXrpomcgeYCCFpSlPFVbUQZCgFQKyx8gEA4THuUsWk5106KSMyihqCg5hRSqQj6ayoJ0DExNaGaqaYYZHv3m5OsbPzjDlTMsvEabIHJqDKRzqJMrnpia9O6Ni97I0vKlp60d6j/E0b7ZtcBa0mJjzUY+VP924t/6CbTNUAOREmPpUio/+EHiR4/S9/AjkEXf91B5Ay8veRfHQ0vJ2CGk380usxoy1jqBk+wgOfIMJI8A+Z1qPQ1CpByCQJgmUkrPIfhsApEKgpEISze9Bl8gQMu6DSWHUELRUXIKYzATLbaabmLw6CQavOGW6WktBqYa9B4HYUHTJq+mkMKvXvwVx0cmL9zW+Gr4u3NO7ziKJWI8dewpupOZ99ibmCyNLOX9zZdjXv87eu6cortqrKm1tVT80R8iA37av/89nIMHM76et4FgoKKF/Yv+gM7wIjItKgsDr047SwoirRMk4wdwRnaCc4yCdhgJiVAKkaolGKaJNAxqFrUSrqikrLaOFedfhD9cVkoXlTBnKDmFMYj6p78BNUSm+5BO0pq5YRJai7GIZZiz9oVhxdtGHUxbfxs/3fnTKQ//9KZPn1ZLSJ/z7IlnM7seXtoo6ouysmIlrc8eo+dXv565nbS6mupPfpLIJW+g82c/Y2jb8zOnmsZhINrMvsXvpivUAvgzO0l5dVo3G4eQfKVwDmFsUdm7INpJomwfUiiC5eWUVVbTeva5hKLltKzbUCoklzDnKDmFMZipLXXEmeKm+OwvJl8/58+nv2A8kzSFgKbzoepUKujmF2+ekvJ6Q9UGLm25dML6zpM7OdyTuXZCUAZZX7WeD8c3MvDtr898c6+poepDH6L87W8jtmMHvffdN+vi8lC0nv3L30tXoJmMHIKAQNSjMh8adGbFOuG6MeLDz0NiH7hdkEWL7oxIU1NICRpMywIBofJK7GCAioZm6pevonX92QSj5aXIoIR5gZJTGIOANb1WQt/QFL1JD06W2w9M3XUE3vTzsUlSTuNhBKF+/aiYTlt/G7988ZeTHiqR/O1Zf3ta2gi8+sOdL9/JEJkVtU1M1tes56/cCxFf+yF6Jq6iQICyS95A+Qc/QLKri87rrsM5NAvBeiGIR6p5edE76Qy0kEkNQRrgL1OYhsHgQBydYUCidQLX6SGZ2AUju4EhCsNhJDzKBMNAKQPXdTF9PkzLpqyyisrmVqoXtZbaTEuYdyg5hTEon24wDWgon+LptW+SG+C5H5n+Yh0vQvfLMxvlj0LjeaOpo6u3XU2fO3mU8EfL/oiLmi+asL6zcycH+w7OfC28CejaYC2XJJdiXfUL3Jk6h3w+7NWrqfzTP0X6/XTd/FtiL70EI5mnYgbKGjnY+nbaoyvJxCEIBZFKG6EEg/0jJGKZ3dS9+sEenPgBSB6jMAXlUTrTFFuFSHUVlRFtaKSquYXFG87BHyorRQclzEuUnMIYVIWm7oM3gOppXp+AS2cYOR/qgcGZagoCateNiulsPbaV3xz4zaRHVhgVfHDNxE6nWCLGE0ee4GQsMzW0iB3hvMh6Nt7Xjrtv3/QH19YSecPrqfjTP8W3fDmxF15g8MknoSND5TUh6K1eyb7l76ffrCaTorLhE1S1hAmELNoP9TEyNAuH4LyMm9wOTh8wkLGNUw7dTfKaMAwvQhACy+dHmiZVDY1UNLVQUd9U0kleYPjd737He9/7Xnbv3s2qVZOwJOcIx3HYtGkTjY2N3HbbbXnZMxfq7FcdbENN+QcxDbCMSV59/AeTnzCJAM5pGO4FPQNvkC8K5/wphKqJJWJ8+5lvT3no5csvpzE0kbl8X88+njz2JPEMpnItLC5uvJg/HzoXHtriMQBOhbIyyi+/nJpPfQrf8uXE29o4ec21DG/bllGU4Jh+empWc2DJ5fSa1bgzOQQB4WqLZRurqagL0nV4gIGe+PTThulrub0k4ltxBp9Ex4+BzswhCMPC9PuRlg1yitSiYXivSQlKoQyFadn4gyFC0XKqGhpZdv5FNCxbWXII8xiJEYfB3hESI/mtLRVKZCeNb33rW6xevXrmA2eBklMYg0jAIjBF7GQqMTnv0eOT8A5VrZ/5Yp0z6Q9IiC7xWFGBLUe28Gzn5N1DZ1eezUfXf3RCLSGWiHHHy3fwcu/MaSqf8LGiYgXvGlhB/MqroXua1lXLInDeeVR8+EMYlZUkOzvpuuFGBp99NqOZBEfZtNedw4HF76Db18BMEYIwIFRu0rQySnzY4fi+Hnq7h2esDXt0FUdJDD2CO7Q7VVCeObKQpoXy+TFMC+1q3KQzeedVWudAa4QyMG0fdjBE/eJltKzfwLILLuSSP/trlp/3Gpaee0HJIcxTJEYcju3v4eThAY7t78mbYyikyA5AW1sbt99+O3/1V38188GzQCl9NAYrasP4TMFAcuKNw5ACL188DgOTsIW/88vTXyg+CF17pz9G2VC9Aspb6Rzq5NvPTh4lBESAfzj3H6j0T2QO3dezj6eOPkViBv1ggaAmWMMbBpsIX/fz6QfUfD5Cb34ztZ/6B6ymJtxYjJ7bb6f/kUegq2tmfiMp6W1YR9vSN9HPRArw8VAWVLeWUdcaoac9Rufhfgb74+gZIgRXxxC8jHZehvgreO2m09smTQuhJAKB67ok3STaSZ7+O41JGUnD8P5TBoZhEoyW4y8ro3VjqcX0TEJ8OIl2wR82GepPEB9OYtrTN51kgslEds4999xRkR3DMEZFdi6++OLR8zKlzv7Hf/xHvvKVr0x6bC4oOYUxqAr5qI8EODk08WnXNA3c8Te8tin6/hdPLPaehu5DcHKGSKGsCc75EwhVs3nnT9nfP7n+wVsWv4W1VWsnrMcSMe49dC9tvTOzsFbb1fxl8C2svfZReHmariHTxL/p3FGHADC8bx+9d9xJsq3No3ueBo4vRE/FUl5Z+nb6qZjeKAH+iEn94gjLzq+l/UAf3UcGGOid2SGgYjhD20Efx4mfYNIZkrGXsixM08by+zCEYnh4EO3CyFBsnJNLtZc6DigFUqEMm0hNFVIqKpuaaVq9nqZVa0tF5DMIls9ASBjqTyCk93M+UEiRndtuu42amhrOPfdcHnzwwbzYm0bJKYxDNDR5KsNvSFoqx/Gt3/dfEw8MzPz0y+DJ6QfXhAU1KyDSyL7ufVzz/DWTHrYkvIT/dfb/mpA2gtSw2vFnic3QYVNmlHGxsZK1N26FHbumPdZYsoSaT31q1CEkOzvpu+024gcOwAy6to4d5ETzhRxpvIj+mQRyJISiFo0ry6luLqP9QC9tL3bR1zMyrUPQOoFLByq5H+EcxEn2gTN1fUMYJoZpYPtDYAgknt5tfGQE13HAHefkpEApA2lZ2D4/yrapam6hefV6lGXTtHI1kZq6kjM4w2DaivqlUeLDSSyfkZcoodAiO4899hi33normzdvZnh4mL6+Pv7kT/6En/3sZznbXnIK4zDVrEJrpZ+q0LgP+8EnJh54yb9MXBuPZBwGe6Z+3fRBwznEfGGueerrdCYmOhCJ5MOrPzxhcjmNPT17ONw3/bCagcEqWc+btsRg246pJ5YNA3P1Kho++1kC6716SbKzk+6bf0vfk09B3zTkcULgBMrobn0Nx+reQL8TmeZY8IUNyip8nHVpM6Goza7HjtJxsJfenpEpKSyEAq1j4O6FkYMMJ9u9FN00BEheDcDGHwxi+ILEh2LEhwcZiQ2j3aSnfKaUFyi4DsIwMEwTyxfA8vux/AGqWlrY9M73lmgpXgUwbZUXZ5BGoUV2rrjiCq64wpuPevDBB/na176WF4cARXQKQoglwGeAiNb6D1Nra4D/wOMovg/YCXwKqALu01p/v1j2pWHIyWvvNeFxH/j4IOhJiqrn/+XMF+k+yLQc/fVnwap3sKX9ee45NLn+8YrICl7X9LpJX2vrb+OO/XdMS5FtYLBGNvGxl5oof/ypqRXUAgHsFSuo+8y/TnAIvffcg3Po0NR1BKWIR6o50fQ62mvPozcRmtIeFITLbWoXldGytoqq5jA7Hz1C+8E++ntHcKf6cwmAPpyRXUhO4DgdkBxmUocgBNKyCARC+MJh7FAIKSV9He3Eh2OMDMbGRAcCaZpo10VaJoZlE6mpo6KhkarGZpRlsWzTBSU5zBImRaFFdgqJTJXXrgUuB9q11uvGrL8N+Bae8trVWusvTbWH1voA8DEhxFh5srcD39ZaPyKEuFVr/S7gb4QQErgeKLpTSLiT3+BC4yOIO/9j4kH+GdpQwXMmR59lytYZMwSLLqTNNPnu499lZBI+nogR4VPnfmrSKCGWiPHbvb9lZ8dO3GmelGtFhA++XEv5fc9O3WkUjRJ5y1uo+NM/wZcKcd1YjN577qHvwQdI7Ns3tTNRinhFA0fX/gFH1RKGE1PTVlh+STBis2RjNcs21eEkHLZuPkj7wV4G+qZwCMI7z0n0opNbQR9lZLgHnCm8h1LY/gChigqqmxcxPDTMUF8Pw319xAb6SYyMnHIIQqAsE8vvQ2sIRiIYpk3rWRtYf+lbSspnJcyIQovsjMUll1zCJZdckrf9Mo0UfgJ8B+9GDYAQQgHfBd4MtAFPCyFuxXMQ43kfPqq1nowS9KfA54QQ7wJPeDf1/d+mXis6fObkkcLIeF+x+7cTD3rtP8x8gYF26Dw49evRVljxNu5qe5i9fZN3KL17+bs5p+6cSV9r62/jiaNPTJpyAq/TqNpXzfuG1rP8iZego2NyOwIBIu94BzWf+DuMylOdTcP79jH42OPED7dNXUewLJJVTRxd8Q7arNVT0yZJCEVNaloi1C8vZ8V5tSTiDltvP8Cxl3sY7k/iTlJDkAYImcSQvSSdF3DibSTjAxMdghAIqfCHwoQqK6lduoK6pcs4tmcPPQf2MzwwQGJk2OsoSncmpdJG/rII9UuWEx8apKyqBisQYP2lbylFBiW86pGpHOfDQojWccvnA/tSEQBCiBuBd2utr8CLKjLZtx34RMrB3JxauxW4VQhxOzAp05wQ4uPAxwFaWloyuVTGqCub/OnPGp9WGprkZvr6DJxCrAs6p5gUln5YcRn7LItf7JqcZO+1da+ddCYhjZlqCWVmGZcNLeHiO47AwSkoLMJhghddRNVH//J0h7B3LyevvprBbdumdCaOHaQ/uoj21ks5aS6d0iEYFvhCFsvOrWHlaxoIhC1i/SO8tOUYnW0DDA+Ocwhesw+27aJkL/HYPkZibcRjJ0kOD56uYpa+hs+HPxSmfvlKaluXYpgmR1/cxcHtzzHc34/rarTrYlgWpuUDXJRpYfn8NKxYRePKNTStXot2XYLR8lJ7aQkLArnUFBqBsXefNuCCqQ4WQlQCXwQ2CiH+RWt9RcrR/CsQBL4qhLgEeB+e7uKU2s9a66uAqwA2bdqUVzazqWoKPnNM+uipH088IDxxmnhSnNgNQ12Tv1a+mNiKt/LdHdfSPjIxsKr2VfMP50w+kwBelHDrS7fSnZg8HVRj1/D65BIuv70dXnwZxqfKpIRwmPAb30jNJ/5utMsI8CaWf3Q1gzt2ePMI42FZuNX1dETWcrzufHpE9eRU1gqCZQaVjWU0Litn1YX1AOzZepyjL/XQfWSA/r7hCSkjZYAUwwh9gKGB/SRiHTjOMNp1xugcC5AKXzCIPxwmWttIdWsrdYuWcvLIIU68vJcTBw8SH4yhXddzCKaF7Q8QKItiBnyp1tIW6pYsK5HVlbAgUbRCs9a6E/ibcWsHST3xj8GDRTJpUiSn6MCpDo+ZZn5yEmqLcz468+bxQTj4EEw6TGbB4ot4Rsd49PCjE14NqzCfueAzrKteN8m5Y2oJnTsnfT2gAlxgr+KyR/swDh6BxEQbVGsrgY0bqfnbvznNISQ7O+n+9W8Y2r8fevsmziMohRMuZ7DpbNrD59BvVk2qMmr4BL6AORodRKr8JEYcXnzsKHufPcFg7wjx4aS3vQRcMP2eco4/EsMZ3oObOILr9JKIx9DJ+Gk2mLYPy+envK6BmtYllFXXINC8/MKzHNv/EoPdPQzHBkCDNEwsJYlU12IFgpTXNVDV1ELDylUlsroSFjRycQpHgOYxPzel1s5oTBYpCCBgjflTdY6XtLTh0v8z8+YD7dA1xTBZtIHO5W/k+t2/ZHiSYav3rngvr2187ZRbt/W38dzx5xhwJvL6VFgVnGet5A8fc/A/sn0iN5GUqBXLqfnoxwhddOFpKaNkZyddv7yJ3gcewDl6dCKNRTBIvHYxJ1e9lQ7/Ejp7YDxDt+kT+MssojXB0eggUGbRe3KIF588yvGXexnoGSYx4oyWBaQB0hL4/Amc+B7ifQcZGewgmRghGR9BJxMe11Ba27gsih0K0bRqLa1nbSAxPMKJl/fReaSNrmNHGOztxUkmsHye4/CHIpTV1BKKRjDtABWNjaXIoIQSyM0pPA0sF0IsxnMGHwL+OC9WzSHCk0wzWgpq0qprL/yGCaLuNRmyH8a6oHeSPL7wEVv5Vq7pfI4n2ifOPmyq3DRtHaFzqJPf7f0du07umtBx5JM+1qhG/uA5gf++pyZ3CItbqfzwhym77E3IwKlrnOo0ehCnre2U0I5SYCiwbJKL13Oo4VJ6rEYGevXpTVUKTFNSvzxKy6oKGldVEKnyY9qK3pNDPH37AToPD9DfP0Q85iIEGDZYfhPLp0H0o+RBRpL7GB7sxkmOoLULSqKUD8NnY1o2lXUNROubKK+vp6Z1MZ1H2mh/+QAnDu0nGYsxHBtGCLBsH8o0idTUsPTc81l67gVo7VFclyKDEkrwkGlL6g3AJUCVEKIN+JzW+hohxCeBu/A6jq7VWk+euziDELQm/kkqgiZnNUW9H7b+ZOJJq9+Z2ea9R2Bwknx/1TJ2Nm7g9h1XTXip3CznE+d+Yso6Qpr07v7D9zM4bm4irMKcYy/jz3ZUEn7kWRga9wgfCGAtXUrN3/89wU3nnuYQwOs06n/gQRKvvHIqQpASAgF0pJLEyvNob72Erm7JYI972miAsiAU9VG7uIz1lzRT2RDCtBWJEYeTR/p56YnjdBzqY3g4iRsH0y8xTYlhutj+AZzkYdxEB0N9xxjs6SIxPIyrHexAAL8vRDASxfT5qV++imWbzkdraD+4n31PP0lP+3GGBwZwhkdASEKV5Vg+P+GKKqL1DSzecA5VTYtKTqCEEiZBpt1HH55ifTPTFITPRJQFJtJclPuMU9PM3QfHvWrBhg9ltnnXftDj2jiln84lr+P6o4/RlTi9gBtWYT696dNsqt805ZZt/W08dewpesZNSFtYnGW08IEtgvCjj0HfuLH5qioib3zjaTMIYxFva6PzuusYev556O31irimCZEIYtV6TjZcxEn/Inq7HAb7k55DEKBMzxk0rY7Suq6GmkVlBMq8ekysL86B505w7GAfJ1/pZ7A/gU66aK2JVoYwfSMIfYjk8GGGe45j2AZuMo5hGBiRMoSUVDU2EyqvoKp5EXYwxJKN3t9mz5OP03W0jcTwEMowcZ0EdjgEjktl8yJqFi8r0VCUUEIGKNFcjINtTBx1N9PRw577oH8cYVzLedPLbqYx0AF775uwHKtcwg1+xROHnj5t3Sd8fGT1R3jL4rdMuWUsEeOBVx5ge/t2BjkVJQRlkE2+FfzZi9WEt2w53SH4fMiaGsrf+U4q/uQjp9UP0oi3tdF5/fUM737RiwxME6RERiLIcy7k2OrLOdYuGO53cRyNkuAqMH2KUMRi5WsaWPWa+lFnAJ5D2PPUcY7s7WZkMIHraHx+A8OUmJZDeeMwg117Geo7SKz7JIn4MBob07YxLAtlWoQrKlm0YSOLN5wzOkAGsP+ZJ+nvOkliZITEyAiBsijBSJRobV2Jk6iEOUUhRXZaW1sJh8MopTAMg61bt+Zl35JTGAfbVNhw2hxxVVpkYftNE09Y+77MNu4+CP3jCtTSz7PN6/jdiSdPKy5LJJe1XsaH1nxoyjoCeNTYD7c9zEB8YPQ8S1isMRr5wPYQZc/sRg8OnqJ79vnwnXceoYtfR/k73znBIbixGMN799F90y8Z2v0iyb4+77xQCNHQAue9gRP1F3Bo7xAjwyM4IxplCUJRH74yi+rGEEvPraF2UeQ0HpnEiEPbS130dcZwEg7CkJi2xPab+KMuJNqID5wkGevAjQ8jpMT0+QlGojSuWkfzmrVoVxMoi0y4uQ/2dCMNk1C0AoQkEGph0XqPslprXZo8LiEjJIaHiQ8P5f39MlZk5/Of/3ze9k3jgQceyDtVRskpjEPAUqdkdlMIWamUUte4InGwFta+O7ONuw5B9+nNWW2VrfzY7aJj+JR8pUJxbtW5fGLj1HUE8IrLdxy4g7a+NlxcJJKQCnGWU8df7m6mbMtO9EBqqCsQAJ+PssveROVHPoLV1DShfuDGYgxs2cL/3965B8dV3Xf8813trla7K8mSZcu2hLFUbAuZGtvYxh5iBxJwDE0IZiAFtxNoSELbyTQtQzuYdJh6Op1OO00f6QPo9EGbIYSSZJKUmRZTmymP8H4Ym4cBYx52QTYCy7be0p7+cY7Wq7Ueliyx0ur3mbmz5/7uueee3/527+/ec8/9/U488QTd+/ejkgjRqioiiVIiv7yWtqbPcrwvxZEPO8lkMpSlS+mN9VE5J8XC5mqSFXHOXlYz6O5ggJ6uPkqiEcqry+g63ks0UUL1ghpq6ko5/PYejrYcpbe/m0RFCqcMsUQZ5TU1zG34JZasHTlbWTxRRjyRgOrZlJWXWw4DY8z0dnVxaN+rZDIZIpEIdUubJ8QxDCTZ2bFjB9deey3bt2+nra2NJUuW0NLSAvgkO7t27aKycoRAkZ8y5hTy6OjuJ1oC3Tlv06aTMZ87oSXvOXrTF0dPuwn+/YS3d5GbO7I1UckP5sxj34kPccEDxYixYu4Ktq3bNmz0U/DDRjvf28mrra+SUYZ0PE0JJVwSW8Y1LyUoffsQPV1dlFRU0C9R2thA+SWXUHXVVUMOFwFk2tvJdHQQq5lDb7qcvuPHic2bR9lFG/lk4YW0f+zo7e+lLF1KR1sPJdEIlXMqWHPFItJViRFDDg9sq6gpI5GKMefsCipryujpPEbH0SSZvtm0vHeMZLqC+YvPpbah8bTfFYglEtQtbZ6UqzxjZtDT5S90UpWzaD/WRk9X54T8jiY7yY4kNm3ahCRuvvlmvvnN/Fe+xoc5hTy6+/pPCfpZN6sM9n4f+nLeAYhWwspfP71GP3kXPtzLwHzNDkXYWdvAi/Eofd19JCIJXMaxbM4ybl93O+dUnTNic62drRxoO0BMMVKxFImSBMvLl7D1/Qbi7W+SSaVRSZRIaZz4ovOouekmkuedd8rdQS6RVMpvl0g0n0ukvILUhWvJzF3IiQ96SPZ30nm8h9SsOLWLFlBeU8bchf7kPhrDxqt3/iq/onYeifI08xoXj2vsP5ZImDMwxk08UUYkEqH9WBsRKfus6kyZzCQ7AI8//jh1dXUcPnyYyy67jKamJjZu3HjG/TankEdprISyRISOdj+/siwKi8o6Yf+zDBpTOu8qqB86KN0ptO6HjiMQiYPL0Fpey4HaJhL0UBYtI6II9el6bl1z66gOAQBBKpqirrKOREmCptlNXFm9kfKu9+icdQRJpC68kPRnLiK5cuWgt5OHI5JMkl6/PjsTKVpTQySZpLe7n9jHR0+5yh9r7Pmh4tXbVb4xFZiM3+FkJ9kBqKvzoXXmzp3Lli1beOaZZ8wpTAbJeAnxaIxopJsSwVlVSc7qfw96jvmMaC4Dqdlw/ldOr8GedvhgN/T3+TezFIEFy0lVN7Iw00EikqBhVgNXL7369BwCMDsxm+aaZtq62tAcsaF+A1WujI5Dx9AFF5A5foz0hg2n5QxyiSSTxPMCDE5GVqpB7dtVvjEFmOjf4WQn2WlvbyeTyVBeXk57ezs7duzgjjvumJC+m1PII5OBuelSwNHb52iaV0FdeXi3oGIedJ+Ac6+EupWn1+CJw9BxFNK10NcJyTnMXv11muMadFIf6aFyPslYknXz19HR10EymszOUEquWkWmvf3kUNAEMdFZqQyj2JnsJDstLS1s2bIFgL6+PrZu3crmzZsnouvmFPKpTsdZVJOkvipBe08/W9cvpDpVBnOW+ny/zsGK6yCeOv1GU1Uwbxl0HYWGi0kuXMc66ZST+lhIxk7dL5JMTqgzMAxjfEx2kp3GxkZ27959Rm0MhzmFPBZWp7i0eR6fdPRQlYyzvL4KSELzl3zsomT16cc6AkjPhQUrYNZCP3S0eBPEUyRhXM7AMAxjMjGnkEcyHuXz59ZyoruPdGk0REeNwpLNfuioND22u4R4ChovHt++hmEYnzLmFIYgGY8ODpUN/mQ+3hP6mexrGMaYcc4hqdDdmBK4/Dn2ozB0mjHDMIxpSiKRoLW1dcwnw2LEOUdrayuJMcyssjsFwzCKivr6eg4ePMiRYfKIzzQSiQT1Y5iebk7BMIyiIhaL0dDQUOhuTFts+MgwDMPIYk7BMAzDyGJOwTAMw8ii6f6EXtIR4N1RKw5NDfDRqLWmP6ZncWF6FheF0vNs59wpsf+nvVM4EyQ955wbPgFykWB6FhemZ3Ex1fS04SPDMAwjizkFwzAMI8tMdwr/WOgOfEqYnsWF6VlcTCk9Z/QzBcMwDGMwM/1OwTAMw8jBnIJhGIaRZUY6BUmbJe2T9Jak2wrdn7Ei6SxJj0h6VdIrkr4d5NWSHpb0ZvisCnJJ+l7Q92VJq3LauiHUf1PSDYXSaSQklUh6UdKDYb1B0tNBn/slxYO8NKy/FbYvymljW5Dvk/SFAqkyLJJmSfqRpNclvSZpfTHaU9Lvhd/sXkn3SUoUiz0l/Yukw5L25sgmzIaSLpC0J+zzPU1WbHDn3IxagBJgP9AIxIHdQHOh+zVGHeYDq0K5HHgDaAb+HLgtyG8D/iyUrwD+CxCwDng6yKuBt8NnVShXFVq/IfS9BfgB8GBY/w/gulC+C/itUP5t4K5Qvg64P5Sbg51LgYZg/5JC65Wn478BXw/lODCr2OwJ1AEHgLIcO95YLPYENgKrgL05sgmzIfBMqKuw7+WTokehv8gCGG498FDO+jZgW6H7dYY6/Qy4DNgHzA+y+cC+UL4buD6n/r6w/Xrg7hz5oHpTYQHqgZ3A54AHwx/iIyCab0/gIWB9KEdDPeXbOLfeVFiAynCyVJ68qOwZnML74YQXDfb8QjHZE1iU5xQmxIZh2+s58kH1JnKZicNHAz/MAQ4G2bQk3FKvBJ4Gap1zH4RNHwK1oTycztPhu/hr4A+ATFifDRx1zvWF9dw+Z/UJ29tC/amuZwNwBPjXMEz2T5JSFJk9nXOHgL8A3gM+wNvneYrPnrlMlA3rQjlfPuHMRKdQNEhKAz8Gftc5dyx3m/OXE9N6vrGkLwKHnXPPF7ovk0wUP+xwp3NuJdCOH2rIUiT2rAK+jHeCC4AUsLmgnfoUmS42nIlO4RBwVs56fZBNKyTF8A7hXufcT4K4RdL8sH0+cDjIh9N5qn8XFwFXSnoH+CF+COlvgFmSBhJE5fY5q0/YXgm0MvX1PAgcdM49HdZ/hHcSxWbPS4EDzrkjzrle4Cd4GxebPXOZKBseCuV8+YQzE53Cs8DiMOMhjn+A9fMC92lMhFkH/wy85pz7y5xNPwcGZivcgH/WMCD/apjxsA5oC7e0DwGbJFWFq7hNQTYlcM5tc87VO+cW4e20yzn3a8AjwDWhWr6eA/pfE+q7IL8uzGZpABbjH9pNCZxzHwLvS1oaRJ8HXqXI7IkfNlonKRl+wwN6FpU985gQG4ZtxyStC9/dV3PamlgK/WCmQA+DrsDP2NkPfKfQ/RlH/z+Dvw19GXgpLFfgx1t3Am8C/wNUh/oC/j7ouwdYndPW14C3wvIbhdZtBJ0v5uTso0b8SeAt4AGgNMgTYf2tsL0xZ//vBP33MUmzNs5QvxXAc8GmP8XPPCk6ewLbgdeBvcD38TOIisKewH34ZyW9+Lu/mybShsDq8L3tB/6OvIkJE7VYmAvDMAwjy0wcPjIMwzCGwZyCYRiGkcWcgmEYhpHFnIJhGIaRxZyCYRiGkcWcgmF8Ckj6I0m3FrofQzFUdE9j5mJOwTCmOOEFp8n8r97DDAo3YYyMOQVjyiPp9yX9Tij/laRdofw5SfeG8p2Sngux+rcH2WZJD+S0c7FO5mTYJOlJSS9IeiDEkULSO5K2B/keSU1BPuhKXz4fwKKwvC7pHklvSLpX0qWSngjx8NfmqHJ+OOabkr6Rp9+zIa7+QN8XyecK+Hf8C0u5oQ+QtEbSLyTtlvSMpHJJN0r6qXzc/nckfUvSLSHI3lOSqof6fp1zjwIfj9c+RnFhTsGYDjwGbAjl1UA6xH7aADwa5N9xzq0GlgOflbQc/wbphfIRRwF+FfihpBrgD4FLnXOr8G8S35JzvI+C/E7gdIZ8zgG+CzSFZSv+rfNbgdtz6i3Hx29aD9whaYGkTfgwDWvxbzVfIGljqL8Y+Afn3DLn3LsDjYTwLPcD33bOnY+PKdQZNp8HXA2sAf4E6HA+yN6T+NAIhjEi5hSM6cDz+JNlBdCNP8GtxjuFx0Kdr0h6AXgRWIZPnNQH/DfwpRBQ7Vfw8WLW4RO1PCHpJXxMmrNzjjcQYPB5fHz80TjgnNvjnMsArwA7nQ8VsCdv/5855zqdcx/h4/2sxce22RT6/QLeqSwO9d91zj01xPGWAh84554FcM4dcydDTz/inDvunDuCDzX9n0Ge3xfDGJLo6FUMo7A453olHcBn6foFPj7QJfgr9NdCULRbgTXOuU8k3YOPmwM+uuq38MMjzznnjoeAYg87564f5pDd4bOfk/+RPgZfRCWGqA8+70N3Tjn3P5YfU8bhY+D8qXPu7twN8nky2ofp30iM2BdJZ3HSUdzlnLtrHMcwihi7UzCmC4/hT/yPhvJvAi+GK/IK/Am0TVItcHnOfv+LD0P9DbyDAHgKuEjSOQCSUpKWjHL8d0I7yOfTbRiHDl+Wz0k8Gx/g71l8VMyv5TzTqJM0d5R29gHzJa0J+5TrZOjpEXHOve+cWxEWcwjGKZhTMKYLj+FTEj7pnGsBuoIM59xu/PDL6/hczk8M7OSc68enfbw8fBKGVm4E7pP0Mn44qmmU4/8YqJb0Cv7O441x6PAyftjoKeCPnXP/55zbEfr8pKQ9+FwK5SM14pzrwT8f+VtJu4GHGXznMiYk3Yf/DpZKOijppvG2ZUx/LEqqYRiGkcXuFAzDMIws5hQMwzCMLOYUDMMwjCzmFAzDMIws5hQMwzCMLOYUDMMwjCzmFAzDMIws/w8Gp089YByUqQAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"import matplotlib.pyplot as plt\n",
"for dv in range(0, 6):\n",
" mask = mdb.df[\"v_u\"] - mdb.df[\"v_l\"] == dv\n",
" dfv = mdb.df[mask]\n",
" plt.plot(dfv[\"nu_lines\"].values,\n",
" dfv[\"Sij0\"].values,\n",
" \".\",\n",
" label=\"$\\\\Delta \\\\nu = $\" + str(dv),\n",
" alpha=0.2)\n",
"plt.legend()\n",
"plt.yscale(\"log\")\n",
"plt.xlabel(\"wavenumber cm-1\")\n",
"plt.show()"
]
},
{
"cell_type": "markdown",
"id": "d14e60b4",
"metadata": {},
"source": [
"In ExoJAX, we need activate \"mdb\". User-defined mask for DataFrame can be used as load_mask option of msb.activate."
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "8bb601cd",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
".broad is used.\n",
"Broadening code level= a0\n",
"default broadening parameters are used for 69 J lower states in 150 states\n"
]
}
],
"source": [
"load_mask = (mdb.df[\"v_u\"] - mdb.df[\"v_l\"] == 2)\n",
"mdb.activate(mdb.df, load_mask)"
]
},
{
"cell_type": "markdown",
"id": "08b96a1e",
"metadata": {},
"source": [
"The following lines are activated and can be used further analysis as usual."
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "e6871c85",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEGCAYAAACgt3iRAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAA+gUlEQVR4nO2dfdgfRXnvP0OCUJH3GKuNIYlvNCEJRFDQmlBARKDFFw5WS0RKpVqVnmO1hxx6+iRVr+pFtYpYgyKCLUIk9ihSEC1GQAHh4f1FAzEBTeEyAUFRRMDM+WN3n2w2M7OzO7Nvv+f+Xtdcv9l7Z+a+d3+z93fumX1RWmsEAoFAIAjBDl0bIBAIBILhQ8hEIBAIBMEQMhEIBAJBMIRMBAKBQBAMIROBQCAQBGNq1wZ0gWnTpulZs2Z1bYZAIBAMCjfffPPDWuvnmvZNSjKZNWsW4+PjXZshEAgEg4JS6gHbPpnmEggEAkEwhEwEAoFAEAwhE4FAIBAEY1KumZjw9NNPs3HjRp588smuTRF4Yuedd2bGjBnsuOOOXZsiEEx6CJmk2LhxI7vuuiuzZs1CKdW1OYISaK155JFH2LhxI7Nnz+7aHIFg0kOmuVI8+eST7L333kIkA4FSir333lsiSYGgApYvX95Y20ImOQiRDAvyfwkmG2xkYJKbZCtWrIhrUA5CJgKBQNAxQsnAJG+SOEwQMhkovvvd73LddddNbK9cuZIvfelLtdo6//zzefDBByvVuf/++9lvv/2c8vHxcU477bRaNgkEQ0eVKKJJMli+fDlKqYlIPsvHnvISMglEk3OQLhTJ5F3vehdvf/vba7VVh0x8cOCBB3LWWWdFb1cg6AqhBOGSu3SayODQQw81youyFStWMDY2RvYhRK01Wuv4vitreDKll7/85bqIe+65ZzuZD5JTGAfHHXecXrRokZ47d64+55xzJuRXXHGFPuCAA/SCBQv0YYcdpjds2KCf97zn6Re84AV64cKF+pprrtFjY2P6zDPP1D/84Q/1QQcdNFF3w4YNer/99tNaa71ixQp94IEH6nnz5ul3vvOdesuWLfqSSy7Ru+yyi37pS1+qFy5cqJ944gk9Pj6uFy9erBctWqSPPPJI/eCDD2qttR4fH9cLFizQCxYs0B/4wAf0vHnztjuGDRs2TMjXrFmjjznmGK211mNjY/rkk0/WS5Ys0bNnz9af+tSnJur827/9mz7ooIP0woUL9amnnqqfeeYZ/cwzz+iTTjpJz5s3T++33376E5/4hPGc1f3fBIIMY2Nj3nLb9W6Sl5UdGxvTgFfKbAnVH+qvgHFt8audO/YuUl/J5JFHHtFaa/3EE0/oefPm6Ycfflhv2rRJz5gxQ69fv36bMhl5ZMhvL1y4cKL8Rz/6Uf2hD31om7paa33iiSfqSy+9VGut9ZIlS/RNN92ktdb6qaee0occcojetGmT1lrriy++WJ988slaa63nz5+vr776aq21rkUmhxxyiH7yySf15s2b9V577aWfeuopfc899+hjjz1WP/XUU1prrd/97nfrCy64QI+Pj+sjjjhiot1HH33UeM6ETAS+sJFGGwSxZMmSyiRRR3+ZzHYOfOEik0FOcyml3qCU+rxSapVS6shUdqhS6lql1Eql1KFN6m9qDvKss85i4cKFHHzwwfz0pz/lvvvu44YbbmDx4sUTz1Lstddepe2ccMIJrFq1CoBVq1bxlre8BYA1a9bwyle+kvnz5/Od73yHu+++e7u6a9eu5a677uK1r30t+++/Px/+8IfZuHEjjz32GI899hiLFy8GYOnSpZWP75hjjmGnnXZi2rRpTJ8+nZ/97GdcddVV3HzzzRx00EHsv//+XHXVVaxfv545c+awfv163ve+9/HNb36T3XbbrbI+wegjxrSTjw7faSbblNJ3v/vdrSP4nLyKz1i+fDljY2PbyYAJeb69sbGx7dpfvnx5c1PzNpZpKgHnAZuAuwryo4C1wDrgdM+29gS+kOaXAFcA5wMvdtXrY2SyZs0a/epXv1r/+te/1lon0cKaNWv0pZdeqt/2trdtV94Vmaxbt04fcMABeu3atXrRokVaa61/85vf6OnTp+uf/OQnE+WzUUo+Mrnjjjv0wQcfvJ2+Rx99VL/whS+c2L799tsrRyZ5e+fNm6c3bNigzzrrLH366acbz8njjz+uV69erY877riJ6KgIiUxGC1WmnLSuP1KPGUGU6Svabiq7ZMmSUlnZMfls22S+oGeRyfkkxDEBpdQU4DPA64G5wFuVUnOVUvOVUpcV0vRc1b9P6wFcq7V+PfC/gXbviYuAX/ziF+y55548+9nP5kc/+hE33HADAAcffDDXXHMNGzZsAODnP/85ALvuuiuPP/64sa0XvehFTJkyhQ996EMTUUn2cN+0adP41a9+xerVqyfK59t62ctexubNm7n++uuB5DUzd999N3vssQd77LEH3/ve9wC48MILoxz34YcfzurVq9m0adPE8T3wwAM8/PDDbNmyhTe/+c18+MMf5pZbbomiT9AP2EbHsRat83pMUQWQH5SWRhC+x1CMDHxtv/rqq71kvYaNZZpMwCxykQlwCHBlbnsZsMxRXwEfA44w7HsWsNqlP2ZkEjoHmeHJJ5/URx11lN533331cccdNxGZaK315Zdfrvfff3+9YMGCiXWEtWvX6vnz52+3AJ/hzDPP1IDesGHDhOyMM87Qc+bM0a961av0O97xjgnbV69evc0C/K233qpf85rX6AULFui5c+fqz33uc1rrrQvwCxcu1B/84AejRCZaJ+syCxcu1PPnz9eLFi3S119/vb7tttv0AQccoBcuXKgXLlyoL7/8cuN5k8ik34gRURTlMaOKrC2TzXk5hiijaJNNR1HmUyaTVVmkr5Oq+i/6tgBvIJPjgXNz20uBsx31TwNuBlYC70plbwLOAVYBhxrqnAqMA+MzZ87c7iSJUxom5H9rDy7HE0Iadckha8f3riuXwzaVsxGHjyw2CZhIz2RDXeLyxciRSWiKGZkIuoX8b/FRlRhc+3xIw+S0TQ7eRhpVHb2vLJQQbNFKXduETPzIpNI0V2gSMhkdyP9WH7GmoIr7XJFGsbwtKjA5RZN+3xF7UWcoIYTI6tYz/V9FWdm2TeaLIZDJVGA9MJtkzeN2YF5T+m1ksmXLlrrnWNABtmzZImSSoupdUFr7kYYrmqhCGra8zWGa5D52dUkSNlldBx9rPTYmekUmwEXAQ8DTwEbglFR+NHAv8GPgjCZtMJHJ+vXr9ebNm4VQBoItW7bozZs3TzycORkQSgyufWVTUGURQBlpmPT52uCyy3aMTRNCFdkowUUmrX8cS2v9Vov8cuDyls2ZwIwZM9i4cSObN2/uygRBRWRfWhw12B4sW7FiRfADZ8uXL9/mVtXsNtmx9EG7bH/iN7bVXQdLlizZ5lMB+XzWZtGGTFa0wSQLxVjhIcC8LP8/ZHnX+S/uD90eHGwsM8rJFJkIBG0iRpQRcheUTU9W1iQvRik+qU4EUaVsNt1WVeYz9VTU1/a2j40+cB13VdCnaa4+JCETQVsIvTPK5w6orE6Vu51sNviQRhMEUYcQquiv6sjr1Glz20UKrnPk6nc+EDIRMhE0iCajDJPjrhpR+JBTHdJwHUdRlv81yWKRQt16eVLzPT9NJtN59CWFuvV8IGQiZCIIRB3CyO/LRtUmx2F7B5PNKZj0FZ1PG1GFLVKoEz1UIYWmHb6N4NrYrnJsNvi2U2fKS8hEyESQosr0QB6+I0EXadgcuQ9phKyP2BxHKJmEykJJwTR697WhDpn1Ydt2Hlx9ztXXqkLIRMhkUiEGKZTtszlxH2IobvuQhomIbBGASR7i+MvK9okU6tYz1SlbAG9627fPFFG3ng+ETIRMRg4xCcM2rVFntOdLDDYysj3w57N2YpOXOeuqUY+L7NomBVM/8JGVbfcBLhur9HHTWlVdCJkImQwSsQjD5dTLIgHfyMLVvs2ZmxywL0nZjjMvrxIxuMjAJo8pq0sKkxV1p2tDIWQiZNJbtEEY+Toux1a2zzWCz9tgasNXb1k0UjWKCI0YqpYVUhhtCJkImXSKPhBGcV9ItGKz20UapnOQlXVNebUdRdimVoq/Nlkft+vCNj3k0lc3H6NsGxAyETLpFLEJw9Vm0aH7EkNGIr62l5FGzAjCpc9WtngsPp+Fzex26fCxo2/bvs7Xdey2vG+5snxo2SpkHwIhEyGTRlHWWWMTRnFflfZshFHFiefldaMYm7x4XCZ9xf2mkXIdEjDJRoFMXERQp14fycRnXwxyETIRMgmGqwPaHF1MwvBt03RR5feZ5EXnGzOqMDl2U3nThW7aH4sUXDLb8Q8xFc9/WX8b8rHbjtX2WwcImQiZ+KAqYRT32Tpp2QWc1+9y5D46Xbb4kkaMqCIvsxFomZ6qstiOsG6E04dtn3Ph6m8uEnLpr5uvU7bO/236T6tAyETIRGvtPx3lQxi+kUe+ToaqD/yZ2jBdhHUjCtex+x5TLFkoIdgimLq2DZVMbNsuR1qHJHzLNUEmJllxX1l/qjrlNVJkAvwhsBJYDbw7le0AfAT4NHBSWRujTCY+0YWtXB3CqOPAQy/SqhGFL2mYbLCN5LoghBBZ3Xqm9Zgiysp0vV2XTIr93Edf3Xxo2SpE4zoHZegNmQDnAZvIfbI3lR8FrAXWAad7trUD8O9p/o3ABcAngMPL6g6ZTHyjC1NZU0cLIYz8NEGx/fx2GUnY2nYRksnh2sjEdYEV68YiiSYIoYqsDgnYZEOHiwhs5YYGF9G4CKcq+kQmi4FFbPv99ykkn+qdw9bvv88F5gOXFdL0tM6fAlcAb0u3Twf+Ks2vLrOj72Ti6tRlI8ei43Q5wGI7Pk7RRAw2x5mvY2u/7IE/n7ZNx24jFhsJ9ZEkJjMBCOLDRDhV0RsySWxhVoFMDgGuzG0vA5Z5tvWf6e+JwAlpflVZva7JJFZ0kW37EoaPgzTVy+u1OfXMVp9Ix9Z+Hb0u0qhKEC7b6spCCaEswgjZbpuEXLpt+1wjbd+8IB76TibHA+fmtpcCZzvqHwqcBZwDvCeVPRv4AsmayXss9U4FxoHxmTNnxj7HlVB0OnWjC9No3uUw847Z5vjKnLovYZiO0yU31ffVG4sgqpQt3l2W2VeUhWybdMfcLuuHdR2yzam7dNv2mfJVyrrISMinOkaKTGKkpiOTKpFHtl0luvAhF1PdYgSSl5va870LqopjzstD9VbVZ7MhD98nxcuctk+ZNtqoQiaubVefLhsMVdUdk0xcMhf5mAhH0H8yqT3NVTeFkklVstDafjts3ejCpqsqYZgcgW2EbHN2Pk68TH9dQsratsnyeZOsjkOvW893u0oEFiuV9S3b/1B2XG0fh08fc5GJ7VdIJgE9J5OpwHpgNlsX4Oc1aUMomZgurCqj5vzFW3ZRm0jB1q5rpG+qU4UcTHbYdLkuaJ/zU9RRlBXPQ/GCr6KjzMEX2+/KAfraW2Xb1Q9NyYYq7fja5cqX7Q/9n4oRvus33x8nA8HQFzIBLgIeAp4GNgKnpPKjgXtJ7uo6o2k7miATH6fgchBFR23bl+23OXKTjbYL18eB+0Q0NkdTLOcqb5O5yLOuLMTZhEQ0VbdjtOHrtG3btnNQZTBUVXcomZSVNRFCDJIxEcyooTdk0pdUh0x8p2nyHa5KpFC8uHyescjr9c3HiChsDt4kj+H8bbrqEkIoKdStV4dMiuc05nYZUbucf516Lt22faZ8lbI+ZGLbV5VkbG2MEqkImUQgk8IJ1Vq7HZqpTlFeZ3RX3Lblq043VZGbjsMm74IQQmS+9cqctk8ZnzbahMueKmTiIo224SIjF/m4fn36bbFMpqPr8xEKIZOGyMQkKzps34imyqjQlzRcTt6lo2xkaXMsdUiiCUKoIqtDCpMRrnMwSufHRjg+JONLLkMmFSGThsmkjDCyOj4kZGpnyZIlpTpMUUAZUVWJRly2dU0SvgQw1AtY0D18SMaXWEz+YSgQMolAJr7O3DVqL3PMPtMKLlLyIY2iPBZB2OR1CME1DWHa1/V2LKdga9NXt+85rFK2ieMcJeTPk41gyghlSOdVyCQCmRROqFVeHOGHkIaPjrw865hdEIRNXudJcR+dfdou7is7Pts+W5u+uk350LImmS/5TFaYBpi263JoUYqQSYNk4kMYxTpV5W0ThMnB+TiOOiQQq0yfyMTneMra6TOZ+OwriyQnC4oRiyuZ3rTQNwiZRCaTMmdadzQfM5LJ2jLZbSIOmw4XMcUkirKR21CS6byXkckQU/E4i33f9pt3rpONYFwzB1nqO6FQh0yAFzr2HWvbN4QU691cNqeRv3BsHcZ2gdradsltDsxGdmV6qshik4DpWEz6+7Rd5TyUORRb+y7drnxo2bK+7Ptf5n9dBDMZUNYH9tlnn65NtKIumfwImGWQ/wXwY1u9IaQmyKQsqigjB5PjN11gLqfuU86lvwsSiFWmazKpUtennb6RiUlWl1BdBDOZiMV1vvoaodQlk+wVJy/JyZYBdwIzbPWGkELIpMpUVJksLy86DF/nXGw7lBDKdLhkscqYnElR1qftMiL3JRNbm766TfnQsiZZFaKpQjJZ+clELGNj7i+K9g21yCSpx+Ekn9LdD/gkcB2wp6vOEFJTkUlRbtpvc5yxSCKUEEJkdUjAJhsyfI7ZZ19f4UM+LnLxIZhi2SGepyrYZ599BhGh1CaTpC6vAR4GLgV2Lis/hNQEmfgQh6lcH0nC1+GP+gUuqI8iueT7r41gfK6BUe1zNl/Qt+OtRSbA48Av09+ngF/ntn9pqzeEFItMTNFIHYLI18+jKUIojihDpniqbjd9cfjqDp3+8dkv2B75/mYjmCrXzShhCIQSFJmMYopJJjGiCN+yNsfvcmamjlhst81t19RY3YvGRco23VXydfa7SMdEQpMZ+T6cv4YmW5RiWzvp01SXkElDZJKHiwxsEUwebX0atk6dJsnEl2iKsBFIX8jEJSv+CslshYlYJkuUIpHJAFPT01zF/XnHYota6jh9nzJ5u3wuzjZTZlOV4zXt6/o4yhxdGZkIyZiREctkIRQhkw4SsAswTvrwJMlNAiuBc4Hryuo3eTeXj8M3ydpw+r5E1uR2nkSqEI3p/NmSS3edvO/+kP8uf7w+JNMnB9M0fEhlFM7HpCETYDfg5QTcGgycB2wi9w34VH4UsJbkNuTTPdr5R+DvKDyJD7wB+Kuy+k09AZ+XhToWU/sufb7bderEJhPXPtcF5TqnVXU3QSa+Mtdx+pLMZCSWsj4w9PMwsmQC/DswLc2/DvgJ8F/AA8D/sNVzJWAxsChPJsAUkm+/zwGeBdwOzAXmA5cV0nTgtcCfAe8wkMlXgF3L7GjiocUqpGCSNVXGVKfYOdvcLu7zJYIibPV8dVfJ19lfRia23yr9K6sz2YjFtlDdN8dbFaNMJnfm8teRvloFmAbcbqtXloBZBTI5BLgyt70MWOao/xGSByi/BXwd2CGVzwQ+76h3KsnU2PjMmTNjndhosjpO36dMnzqiCb5EU4SLQPoAF+kUf2OQTFYuK9vHcxITQ3pq3BejTCZ3A7ul+e9lTjvbZ6tXlgxkcjxwbm57KXC2RzvvIBeZACuAV/nY0MY0Vx51SEDgfzfX0FGHZHySKXIaJQzhVlpfDIFItK5PJicAN5O82PFjwFeBk4DzgY/b6pWlWGQSkpq4m8smC5kC6guanCoS2OEiGV9iycqM4jkfigMugyvq7Nux1CKTpB4vTonk/wHfAD4LvM5VpywZyKTSNFeMFPM5kzKyKEYqru3ivjaIyIcQbDZWyZtkvqTTtwuqa+SnsfJOxzdSGSWMAqEM6Rhqk0kTyUAmU4H1wGy2LsDPa9KGmGRShSzKtmPWtXVEF9k1mY8hs43U+3jRtQUTsbjSkiVLRup8DckRmzA0+0Mikz8GziZZ6P4P4KPAi111Stq7CHgIeBrYCJySyrPX3f8YOKNu+76paTLxGSX6jCJjElHVNkLtr3Kced0uMin7newkkxGLz3kfFQx53cT2X/XZ9lpkAvwT8EXgRGA1cCbwTuBWat4a3JcU47O9TRFCU0QUqsNmr28+7/BjHVsVkplMBJMnFdv57vPX/KpgaCP7PIZoe10yyd8aPBX4fprfk8JDh0NLozrN5Rrp+JBOk/m6sqrkYyObfH6yEIvP9NfQz8EQHbLWw7W7LpncDuyV5mcCN+T21b41uA+pTTIpdg7Xtms9o2y7rGzVNvJ5m41V8iZZVYIx/VYhm6xOsX7fL+BQjI0N62t+vhiqQ9babPsQUJdM3kLytPu3SZ5+PyaVPxf4sq3eEFKbd3M12XYVIrLJfQihKfiSTpF8YpNMVm/UiWUoX/PzxVDJZMjEXotMknrsBRwI7OEqN7TUxCvo+wZbxxxCh7WhSC5VyKYs5R3TKJPKkB1ZEUMjk6xfDcnmImqTSVKXHQ2yaWX1+py6IhPX6NuVF/jBRTJViGVoF3gVDPnupzyG5pRdkXJfbTah7jTXH5PcvvswyXuwZuX23WKrN4QUg0xcjt+2rzhi9s37TAGVTRlNZuQjjeyi1tqPWIbmZMswNCdsw5COw0UkQ+tfdcnkJtKHB0leeXIfcHC6faut3hBSDDKxOX7Xvrpk4qrns69semiyoUgsxWku0wU/SudqSI7YhKHY75rW6qO9PqhLJrcXtueRfHPkDRKZVCOTGMmms8o+2+9kJpl81OLzH4wChuKMTRiS7a6+1KW9Ibrrksk48PsF2QzgNuBxW70hpLpk4hplhJKCKR8zFad3yn4nM7HYzuHuu+/etYnBGPpIeSg3ELjs7NrWkIFRXTI5AlhokO9OC688aTJNhmmuqsRnWkvIt5tdBF1fCG1gzPFcRh8dVxUMaWRfxBCIpM/9Jh+B10UtMhnl1BWZuBbLXfk6ZFL2W4Vssjr57clALENwXlUxVDLp+39RdNTFa6fraa1Y/7uQSQNkYnP8ZftCdJna9dnnQy7Zb1VyGXVS2X333Y3nYGh34WQYIpn0+Xbm/DVmO7dd2OkajEpk0jMyGRJsxOMim6rRyyiTSt9Hxb4YGpG4RtR9IZK+kUiG/AxCzP88iEwwvCHYJBtSmmxkUgYTydSNWvrqmEIwNCdsw5COwzWQ6QORaG2/JvLXT9vIDxCbGPSFksl2twGbZENKQiblMEUvvhHLqEUpfZ5m8cVQiKTMGfbhnPu8hbvt8+qK4ooDwxDUIhPg9cCngZ8BZ+XS+cCNtnpNJ+BQ4FpgJXBoKpsJfA04Dzi9rA0hk+rIOmvZxd63Cz8GhuKIXRjCVJ3PIKVr+8ocdlc2FgkjTxwxbatLJguBk0jeHHxSLr0J2NNWz5VSZ7+JwvdQgKNIHohcV0YGwBLgipTUXpzKjgFOTPOryuyI8XEs0+9kQtmF1RcHEANDJ5MhEInW7gFKl7aWrS3my3RhV1n/jGlbLTKZKGB40WPdBCwGFrHtN+CnkHyudw5bvwE/F5gPXFZI04Ed0nrPAy5M83sDa4DvACeX2RFKJrZRQH7k3rcLtSn4kMrQI5Qhk8kQiMT1bIbW3V9L2fnqQ9/OnwsX8TZ1zkLJ5CUkn+29B1ifpbJ6jvZmFcjkEODK3PYyYJlHO88CVqf5DwCL0/xqS/lTSZ7qH585c2boCbX+5vOTiVjKSGXIhDLUNZMh2N3Xh/x8Rv1t2leMjmz5vKwJhJLJ94DDgTuAfYDlwD+W1XO0VyST44Fzc9tLgbMd9d8EnAOsYuuayX4p4a0E/rnMhjqRic+UTpFMbMQyyuTiOk9DPe4hRiZ9J5Li2lsx39W5LVsTzPpC2yieI5vfafq8hZLJzenvnUVZnRRKJjFSzGmuuiQzyqTS19FmHQiRxEWZs+7KRtu6iEnW1n9fFh3lbWzLplAyuQ7YAfgP4L3AG4G1ZfUc7RXJpNY0V0hqas3El1SKTqnPjqkuhuiETRjScbimGvtAJFq7B2BdP+RX5rTb+M9910Sy/W0jlEwOAp5D8sbgLwJfJf2uSZ1kIJOpJOsws9m6AD+vbvs+qam7uUwRSBVSGSUMyQnbMKRj6Pt6VdmaWlc2Zb+u/7nN/zt/Lkz5vKyLfhhEJjETcBHwEPA0yVccT0nlRwP3ktzV1fgbiZt6ziQfZeTnVvOd00UofbjoY2EU1k6GQiZ9JhIXiWTyts9n39ZFfKKRzJ6u+15vyKQvqa2HFl3E4rr4u+4wsTAUZ2zCEGwvG+13aattDaL42/b6g8mGLuyqQiBt2OMLIZOOyCSPPLGURShdhf2xMQSHbMJQ7Hb1n65szZOI6zy2bV+eJPow8jcRmSvfFwiZ9IBMiigjlaF/1W8oDtmEvj/oVzZ11KWdfSMRnymtfLk2bHHZ03Q0EtpmEJkAzwX+D/A5ktehnAecV1avz6kPZJJhbGw0v+rX59tTXeg7kfS1r/isjbRpS/brctht2NQHAskj9JhDyeQ64GPACcCbs1RWr8+pT2SSoe9OrCqGFpmUOcOuURxhZ/nst8tpLR8SadM+01SRSdaGTSa9rnwTMK0X1UUomdxWVmZoqY9kovXofNWvz07ZBNcItutz7zNN09XUUWaX67ct23wigLamtPoYjcS6HkPJ5MPA0WXlhpT6SiZDc8I2DOk4+kokeRKx2diVfXmi6JLgqjjtJu0ptt01gZjsMeXroBaZAI8Dv0x/twC/yW3/0lZvCKmvZKL1sByxCUOx3zU90wd7XSQS6hDqom9Rks1JxnSeVe3owhYTeTTVd4Iik1FMfSaToS5cZxjC2o/rQuvaVp81iDbty08NuZxTWwTXlymkoh0+zrspW1xkVtwOtSGITICrfGRDSn0mk6GM7E3oO5H0bWRdtK2MRPLH0IY9Wm/riEzE0YZdfSGQPFx9qGlbqkQiRXtCUYtMgJ1JPjp1O7AnsFeaZgE/stUbQhIyiY++E4nWbkfUJYlktrl+27Yvr991vtpy3FXyTaFIajbdTdtSlTxi/kd1yeRvgA3Ab9PfLN0OvNdWbwipr2TS5zl8F4YwNdfX5zPKHEMfp7TasKkP0YgpAiizow1bysirSTKrRSYTBeB9ZWWGloZGJkIk9dHXT8KWTWu1PZ2ltd1Zmqa8urKlaQLJw+Wk27Qj09dFJGKwI4hM3mRIhwPTy+r2NQmZxEGficQ0fZTPd3FOy0b+see3fWBzkPlzlaFN5+2Tj40qaxFtEmtmiylvKtskQsnkP4Gfk3zH5KvAI8C3gPuApWX1+5j6SCZCJHGQJxHbOW3bRtPIvuis246U8jZ1Gbn1KRpx9Zki0ce0w9RWl0TmQiiZfAt4Xm77ecCVJIvxd5XV72PqI5kMYQFba/f0jDzkZ7cp+7X9x22uQ/g4qabsqbMO0eS56cNahKldl+4ufUIomdxT2FaZDLi1rH7sRPIJ4Y8AnwZOSmVzga8AnwWOL2ujb2QyJCKxXfhdRyRVRpV9sKntSCSzp5g3ydqywaW/KVt8yazJtQjbFFZb/TfkOELJ5F+By4CT0nRpKtsFWFNWv9DWecCmYkQDHAWsBdYBp5e08UbgAuATwOGp7G+B16T5S8vs6BOZDIFIyqZEuo5IfC7CtqePyqYt2yI3n2mkvM1t2NBlZORLZk3bULV/xLQlpO+FkokCjgf+JU3HA6qsnqWtxcAitv0G/BSSz/XOYes34OcC81MSy6fpwOnAX6V1V6e/04HPAGcC3y+zoy9k0tfbVfNwXfhd2tknEvF12G04S197mjw/dUf/bdnS9lqELRKxbdtkMXSHth1EJrETyUOPeTI5BLgyt70MWOaofyJwQppfVdg3Bfi6pd6pwDgwPnPmzNonMxb6uoCdoSwa6YpITOsQpt+up4/atMnlrHzyseFylj6ONLYtvtNHXUUipnMQy5Z834txDQeRCcmtwPcBvyDCix4NZHI8cG5ueylwtqP+s4EvkKyZvCfX5ueAC4E/KrOh68ik7NmHLuFDIvlyXdjlujDatKsv58rkoMv6V2yb6kxlNW1HZkuV7S5siHkOynTbZL4IJZN1wB+WlfNNoWQSI3VJJq6IpA9TW2WOoCsb+0Ii+fWFth22yQ6t/Z13k/a4+ktbJJbZUSUKCLXHNo1UNRKJhTLdWZmA9oPIpHQNokoykEmlaa4YqSsysf3JfZjayhxylyNsk03ZbxcXZtEOrbe9ELO8SdaGHT7npCl7fKfWmrTBR1fTkZDPyL8pG2xE5toO1R9KJp8CVgFvJfcUfFk9R3tFMpkKrAdms3UBfl7d9n1SF2TS1zUSF4lk8i6ntMocZhu2FR2la9TX9Oi/LF/cbspxuUbATTrw0EikCf2ZvO0BTxe6Q8nki4Z0Xlk9S1sXAQ8BTwMbgVNS+dHAvSR3dZ1Rp+0qqW0y6SORFENe22/b00cmG7qwyycKaNoWk+Nq23nndftuN02obU5jVdWflYltg08U0pTuQvv9uZurD6ktMimOsPP5tkf7JptcF0Xb9pWNstogkSoEkre5KfiO/tsgs7YjAJMNmR1VtmPq9tXXxLnwiUKa0l2wIygyeSlwFenUFLAA+Puyen1ObZCJywl0EZH0lUR8prTy5ZqEjShs+SZsMkVoZXY0Bdd/0ubou6tIxMeBm/6HUBvqRiExdJchlEyuBl5B7tUpDPSdXFlqkkzKnGNXU1tdjSxNyBOEa9Tdhk1Vo5HYF6vvVFbxP2qayDJbTHnTdgz4OMymyKMv00h5HV1HIRbbgsjkpvT31pzstrJ6fU5NkEkZiXTxx2d2+ZBI0yMa24jbRBxN29Q1geRRJfpog8h8ptVi2WIjry4GPV058D5HISaEkskVwIuAW9Lt44Eryur1OcUmk6wDui7ELqaP+kAiGYpO0nWOWgjVK+Vjo0hmvg48NspG/U0TquvYbc67qUikKwfe5yjEhFAymQP8F/AE8N/A94B9yur1OcUiE99opC2HbRr9237bsKlKBNCkPb5TSU05z75MZVWJRJqIAsocuK9DD9WbtTsZo5DQdmqTCcm7rv45ze8C7OoqP5QUSiZlJJKPVNqEzyinTRLxcdpNjbjqjP6bOjcup+njUGPaUYU8mhiFVyWvGDb4OOzJEoWEth0amdxQVmZoKZRMfJxTm0TiQ25N2+OKilz5Nhx4FWceA1UigKYcuM8ouEkiq6M/xrG79HYZheXtKJO12QeqIpRMPkvyDZOlRHgCvg8phEz6QiKZnjJ7mhzl5FF2wbZBID7Tak1PJWW6fR1XU46jjgONeS66cOC+evPlMoxyFFJmQ8U2+/EEfF9SHTIpuyjaJhHT6Noka9IuX+fdlB0mB9rm6N9lR5XtGPDR0WUU1IT+kNF/rP/AdgxtRiGhNlRBKJm82kc2pBRjmiv/2+aUVl53351301GRy4Ha8rFQJNK2I4AqU2mxR8F1oqCY+mPojXV95NvuKgpp04ZQMrnFRzakFItM+jal1bQ9VZ1304TmcqDFCyw2XGTVdCSU6ag6hdSEAzVt22Qh+kOikLYjsKbsqBKBNGVDLTIheTX83wI/Bd6fS8uB2231hpBi3c3VJKpMITUZAXTtvH2joaIDbdKGKnY0qb/KdmzdQ4pC2rahKTuqRCAN2lCLTJYAYyRv+R3LpfcDL7HVG0Lq+kuLNtg6ginf5Mjbl8iaHv1n+m3bbTnQMicSOxKpM43UVDRUR3eo/iFFIbY+F/v4XbqasMGGWmQyUSD3gCKwA7BbWZ2+p76SSdFRuhxXk9FRGZHZtmPA5EjbcuJ5na7toqypUbivfq2bv63WJZvMUcioRiAO24LI5MvAbiQPLd5D8h2SD5bV63PqE5n4RAFtRABV7YhlSxXyaCoaqmpDViaGHX2ZRsq3WVV33eOvMgJvKgrp2gZfAnfJmxxYGmwIIpPb0t8/Bz4O7AjcUVavqQTMBL4GnAecnsp2AS4APg/8eVkbXZNJVcfdVgTQBZGVjXKbHgVnbbY5jVPUnbXXpv4+TCVVGYE3Nfru2gZfAm9Kf1WEksndKYFcAixJZbUW4FMC2EThFfbAUcBaYF1GEI42jgFOTPOr0t+lwJ/kZa7UBZnYOoVPvm0n3iSRhUQipvqxbKiy3aTupvUX22zLgYaOwJsgzzZtCI2AQvXHQiiZnEbygsfLAQXsA1xbVs/S1mJgEdt+A34Kyed657D1G/BzgfnAZYU0HdgbWAN8Bzg5bWMZsH+a/3KZHV2QSdFB21KGJjpOkdBsjqRJO7ogjz7ZUKY7KxNTf5+ikK5G4MU+3bYNVfQ3cfyxEEQm21VICGVq1Xq5+rMKZHIIcGVuexmwzFH/A8DiNL86/V0KHJvmL7bUOxUYB8ZnzpwZ+xwbUXU6K7bj9J3KKjrQmHZ0GQXUtSHW8fs48SYdeFFH36IQm762yLMJG7qOgJpGVDIJTQYyOR44N7e9FDjbUX8/YDWwkm3faPxFkveIdbpmUpVAmhyFuBxXU06862ksmyOpYkMs+Djx2LqHFIV0ed6LZWPr7kp/0xgpMomRmiQTm8O25ZuOAnymsmLaUTXyaGoU3rYNdZ14iG6XzjYceZVReBPkFao/xIYYEUiI/q7QdzKpNM0VI8UkE9+ppCYIJCQK6IMNTenP5E3Z0LUTL+qsI4vx/1cZhTcxAu9S/6hHIDYEkQnwbOD/Ap9Pt19Cuj5RJxnIZCqwHpjN1gX4eXXb90kxX6fiO/pvYgRSNQpo6oJucxqrqv6sTEwb2nbiZaPgrqKQKqPw2BFYG/pddarorqu/jwglk1XA32UEkJLLbWX1LG1dRPJ6lqdJHn48JZUfDdxLclfXGXXarpJiveixmC/bFwO+RNaUI7dFAr7bMXX76guxoQ9RSOgIPGYU0sUovKsowHT9tH3sfUMomYynv7fmZJP6RY+2KKSpqaQ6umPpD53GinkOfBxabGcaQhyxRsNVRsGhzqxKFBB7FB4jAqmjv6remLqHhlAyuQ74PdLXzgMvAm4sq9fnFPPjWE06T1N7vlFRLFSNPGIdv69Tjam/rk6brA5CI4CYU3htRV51dMfUX1VvTN1DQyiZvBa4GtgMXAjcDxxaVq/PqclprlBUiQRiO3KbM60aCcTSm+1ry6H1NQppYhQcIwoY9QjE1bdGPQKxIYhMkvrsTfIak2OBaT51+pxikknsTlUlEmhKd9vTWCGj/5gOrWmdtrbbjgCKbbYZBXSpu8tjHhXEIJM/AF5F8jqUxaRPoA819eXjWF1FAjb9Ps60iUhkaFFIVVQZEccmr66igK5190nvKCGITICPkUxt/SfwjTRdWlavz6mrtwbXcaKxIoE6xJW3MVR/qN4Q3V1HIVVGxLFHw5MtAin+XxKBxEUomawFdiorN6TUFZn4jPKbiARCdIU6Ux9H7uvY66KtKCTGSDwmYVbVXVd/nyIQ138oEUg4QsnkCuA5ZeWGlNoiE1Mk0lYkECMKioG6ept2ak0Q52SLArrUXZXAYumd7Aglk6+SfGfkHOCsLJXV63Nqi0zKyKOpaRUf3VmZmLrrTic1MSfeRRQyylGAq3zXEYgPgcXSO9kRSiYnmVJZvT6npsjEFg34bsfS66sr9kitremkTIdNv488VhQyGaIAUztt6I5JYII4CCKTUUwxyMTmzKpEIrFH5G1GQXWjkFC9xXabdmqhU2euNmLpjak75vRRzGnDrL02iVuwPWqRCfCV9PdO4I5istUbQopBJj7OpAnyaHsqqcy59CkKacKptTl1VldvTN11po+aIOuydpsiMIEbdcnk+envPqZkqzeEFJNMqkYDVeByom05tioj4dAopEviihEJDCUCiT19FIuss7aaJjBBfdQik1FOdcmkqjOLHeZXkcVybG2OzNskLld7fY4EYuit48CbGBiVtdsUgQnqo25k8jjwS0N6HPilrd4QUlvTXFXQpyikrnPxvcB9HEtTxzyUtYgmp3HqOPCYAyOJPoYLiUxaIpMYI+M6TrSJKMQlc8l9EGM0XvdcT6YIJN9e1w68Tj+S6KN/GDkyAeYCXwE+CxyfyuYAXwBWl9Vv6m6ukLq+jrztKCR0pOp7rHXkVfXW0dHVWkSTEUhTDrxYt02yFrSDXpEJcB6widyne1P5USSvblkHnF7Sxt8Cr0nzlxb2tUImIWgzCqniVGNHW/k2QxxLHRu6iASKbXS5FuFqu+3ow7VPoo9hoW9kshhYxLbfgZ9C8sneOWz9DvxcYD5wWSFNT9NngDOB7xfa7xWZhEQhsUenTa69hJKWq41YemPprqozll6fdn0ILHZUbbPFd59gOOgVmST2MKtAJocAV+a2lwHLPNqZAny9IDOSCXAqMA6Mz5w5M+LpdSMkCqmK0CmlGGs+WTtNH2sdvbF0V9UZqrdPDrwOeWXlBMPHEMjkeODc3PZS4OyS+p8j+fLjH6WyvYGVaYTjJKK2XqeitfmCbjoKadKpxogERjUCcZ3TWESdtdUkgfWJvAT9wsiRSWhq6uNYbUUhXU0pFdtrc4F1CBFI36KPGP9x1k7T0ZdgGBgCmdSa5qqbYn62t0zeRBRSxcl1GYG42qlatkm9fdDZpgMvs7Nt8hIMB0Mgk6nAemA2Wxfg5zWlP/Y34JuKQrqaUiq21UUk0JZe09x+FxFImw7c1k8k+hCUoVdkAlwEPAQ8DWwETknlRwP3kqx5nNGkDXXIpK6TiTVabcrJ9SkCqeNQY0Z5beis+p+W2RhqQ9X9En1MbvSKTPqQ2pzm8oXpIm1jmqXrCKQNh9rmgnLdKaTYDrxOxOVji2ByQ8ikJTKJMWJtwqF3sRYRU2cVvT7t+jrWNqeQ2ow+ynRJ9CGwQcgkMpmEXGxVpnhikVYXaxFt62wz+nCh7hRSjOjD1w6JPAR1IWQSmUx84EMasRfvu1iL6Mv6R93oI8ZUVldTSG1GXAKB1kImnZBJ1akc2/RIXZ1tLe62qbOJkXisqSyffVX11V1/8bFDIKgDIZMGySRGBFIsH0Onj9zVli/a1Nll9OHbXpORh9ZyC6+gWwiZNEgmsUijioNtkqh8bGtSZ1+jjy6mkOquv8S2QyDIIGTSEpmEyIuo41SbWItwtduEzskWffRl/UUg8IGQSWQyaSMCqeNUQ5xLm4516NFH01NZvnok+hC0DSGTyGSSRyhpVF3/iKHTp80mp3WGHn3IIrpgskLIpAMyqVO/jZFxTPKqi8kUfZja6WL9RSCIASGTBslkSHdDFdttMyLo6i6krqOPkEV0gaBvEDJpkEx80ObdUCZ9ebRFXm1FH6byXUcfsoguGFUImXRMJkNY/+jLekRVhIz8Y00jySK6YLJAyKQlMmlzPSK2I29zPaKtBWyf/XVt6GraTiDoEkImLZFJnfWIuiPTthaWu44+Qhaws7JN2OC7XyIPwShh0GQCzAG+AKzOyXYBLgA+D/w58IY0vwo4sqzNNsjER+5CU7eUVnVubZBW1wvYfYh+BIIhoDMyAc4DNpH7RG8qPwpYC6wDTvdsK08mS4E/SfOrcvI9gS+UtRX7bq4m1iNCRuSxHHlZW7FIqy5hldlQV39VGyT6EEwWdEkmi4FFbPu99ykkn+adw9bvvc8F5gOXFdL0XL08mSwD9k/zX87JPw4sKrOr7cgk9t1QZftDHHlWv4kpsyo2+OrwOVaJPASCOOh0mguYVSCTQ4Arc9vLgGUe7RQjk2PT/MWAAj4GHOGofyowDozPnDkz8ilOUNXp1LlluK4uG9ogra4XsEMjDyETgSBB38jkeODc3PZS4GxH/b2BlWk0syyV7QJ8EfgsyZrJacDNabl3ldnU9t1cNjQ9IreV7cM0ku/+rhfQY9kgEIwCBk0mTaQ2nzPp+m6oUD1DnUaSyEMgiA8XmexA+/hv4IW57RmpbCSxYsWKbbaXL1+OUgqlFMBEfvny5duUGxsb89ZRrBsTPm0XjzGr53OcUO1YfbF8+fL8gGYi35Z+gWDSwcYysRLbRyZTgfXAbLYuwM9r2o58ih2ZdB19mNqJ+ZxLHyOPKuhav0AwKqDDu7kuAh4CngY2Aqek8qOBe0nWQc5o0gZTik0mRWcVe3G5a2ceQlYx9IdC1jwEgjjojEz6mpomE999vk6ua2ceSlbizAWC0YCLTLpYMxkJVFkTcLVhyvvWzf5ECFsTMNWJcXyu9gUCwYjBxjKjnNqMTOpGHxJ5CASCvgFHZKJ0OrKdTDjwwAP1+Ph4tPaUUoSeR1cbZe0vX748aPRf1n6M4xMIBMOHUupmrfWBpn0yzRUBPreWNjmVVHfqqctbdwUCwWhBIpOWEDL6D408yiCRh0Ag8IFEJi2g6wcHBQKBoEsImURCH58C94VMYwkEglDINFckyCK2QCAYdcg0V0OI+SyGQCAQDBlTuzZgyMgvjJdFHjKVJBAIRhkSmbQEiVYEAsEoQ8gkEiTyEAgEkxlCJpEgkYdAIJjMEDIRCAQCQTCETAQCgUAQDCETgUAgEARDyEQgEAgEwRAyEQgEAkEwJuXrVJRSm4EHIjU3DXg4UltdYMj2D9l2EPu7xJBth+7s30dr/VzTjklJJjGhlBq3vatmCBiy/UO2HcT+LjFk26Gf9ss0l0AgEAiCIWQiEAgEgmAImYTjc10bEIgh2z9k20Hs7xJDth16aL+smQgEAoEgGBKZCAQCgSAYQiYCgUAgCIaQiQFKqZ2VUjcqpW5XSt2tlFqRyg9XSt2ilLpNKfU9pdSLU/k7lFKbU/ltSqm/zLV1klLqvjSd1KHth6W236WUukApNTWVK6XUWUqpdUqpO5RSi7qyvab9hyqlfpE79/+Qa+sopdTa9NhOb8P+VO8UpdStSqnL0u3ZSqkfpHasUko9K5XvlG6vS/fPyrWxLJWvVUq9roe296bPl9j/3tR2rZSalivXq35f0fbe9XkAtNaSCglQwHPS/I7AD4CDgXuBP0zlfw2cn+bfAZxtaGcvYH36u2ea37MD218F/BR4aSr/R+CUNH80cEVa72DgB13ZXtP+Q4HLDO1MAX4MzAGeBdwOzG2p/7wf+HJmF/AV4M/S/Erg3bk+tDLN/xmwKs3PTe3dCZidHseUntnemz5fYv8BwCzgfmBarlyv+n1F23vX57XWEpmYoBP8Kt3cMU06Tbul8t2BB0uaeh3wba31z7XWjwLfBo5qwOQJWGz/HfCU1vreVP5t4M1p/jjgS2m9G4A9lFLP78L2mvbb8ApgndZ6vdb6KeBikmNtFEqpGcAxwLnptgIOA1anRS4A3pDmj0u3SfcfnpY/DrhYa/1brfUGYF16PH2y3YZO+g1sbz+A1vpWrfX9huK96vcVbbehkz6fQcjEgjTkvA3YRNK5fgD8JXC5UmojsBT4aK7Km9NwebVS6oWp7A9IRtQZNqayVm0HbgSmKqWyJ2aPB8ps7MR2qGw/wCHptNgVSql5qawr+z8J/B2wJd3eG3hMa/2MwY4JG9P9v0jLD8F26FGfT/FJtrXfhb71+0/ibzv0q88DQiZWaK1/p7XeH5gBvEIptR/wv4CjtdYzgC8Cn0iLfwOYpbVeQOL8LjA02RqKtgPzSKZR/kUpdSPwOMlov5eoaP8tJO8LWgh8Gvha6wanUEodC2zSWt/clQ11UcP2XvX5SXbue9Pn8xAyKYHW+jFgDfB6YGEaoQCsIpnLR2v9iNb6t6n8XODlaf6/2XYEPSOVtYKc7Udpra/XWr9Ga/0K4BqS9R+XjZ3aDn72a61/mU2Laa0vB3ZMFyu7sP/VwJ8qpe4nmWI4DPgUyRTKVIMdEzam+3cHHhmC7T3s89vZr5T6d0f5PvX7Srb3rM9vY5ik7ReyngvskeZ/D7gWOJbkLZ3ZIvApwFfT/PNzdd8I3JDm9wI2kCzk7Znm9+rI9umpbCfgKuCwdPsYtl2IvLEr22va//tsffj2FcBP0mOZSrJ4Oputi5HzWuxDh7J1IfUStl3E/us0/x62XYD/Spqfx7YL8OtpaQG+gu296fMu+3Oy+9l2EbtX/b6i7f3s823+yUNJwALgVuAO4C7gH1L5G4E70z/pu8CcVP5PwN2pfA2wb66tvyBZQF0HnNyh7WcCPwTWAv8zV14BnyG5C+RO4MCubK9p/3tz5/4G4FW5fUeTRDA/Bs5ouQ9NOAWSu2tuTM/jJcBOqXzndHtdun9Orv4Zqd1rgdf30Pbe9PkS+08jWTt4huSGmXP72O8r2t7LPi+vUxEIBAJBMGTNRCAQCATBEDIRCAQCQTCETAQCgUAQDCETgUAgEARDyEQgEAgEwRAyEQh6CqXUcqXUB7q2wwSl1HlKqU1Kqbu6tkXQDwiZCAQjiPQV601e3+fT0gscBcOAkIlgJKGU+qBS6rQ0/y9Kqe+k+cOUUhem+c8qpcbVtt9NOUopdUmunUNz35c4Uil1vUq+q3KJUuo5qfx+pdSKVH6nUmrfVL5NZKGSb7HMStOPlFLnK6XuVUpdqJQ6Qin1/fQbGvk3BC9Mdd6nlHpn4fhuSl+0mNk+K/2WxZdIHvjMv1oDpdRBSqnr0hcE3qiU2lUl3yX5mlLq2+lxvFcp9X6VfFfjBqXUXqbzq7W+Bvh53f9HMHoQMhGMKq4FXpPmDwSeo5TaMZVdk8rP0FofSPLU/RKl1ALgv4BXKqV2Scu8Bbg4fffR3wNHaK0XAeMk35/I8HAq/yzgMzX1YuDjwL5pehvwR2nd/5Mrt4DkPVmHAP+glHqBUupI4CUkr9LYH3i5UmpxWv4lwL9qredprR/IGlHJR61WAX+jkxcEHgH8Jt29H/Am4CDgI8ATWusDgOuBt3sci0AgZCIYWdxM4mR3A35L4hgPJCGTa9MyJyilbiF5fcs8kg8JPQN8E/iT9AWHxwBfJ3l/01zg++nr8U8C9snp+4+c3lke9m3QWt+ptd5C8mqMq3TyOoo7C/W/rrX+jdb6YZLXlrwCODJNt5K8QXZfEhIBeEAn3+co4mXAQ1rrm2DiZYHZq+XXaK0f11pvJnkN/jdSedEWgcCKqeVFBILhQWv9tFJqA8kXAa8jedfXH5NEBD9USs0miQIO0lo/qpQ6n+RdWZC8ufW9JNM441rrx5VSiuS7Nm+1qMzeoPs7tl5Xz7DtgG1nQ3lIvmHx21w+f10W33ekSd4r9U9a63PyO1Ty2d9fW+xzwWlL+q2SjGBWaq1X1tAhGHFIZCIYZVxLQhjXpPl3AbemEcBuJI73F0qp55F8YiDD1cAi4J0kxALJC/VerZR6MYBSahel1EtL9N+ftoNKvjE+u8YxHKeU2lkptTfJSwBvAq4E/iK3ZvMHSqnpJe2sBZ6vlDoorbNr7tXyTmitf6q13j9NQiQCI4RMBKOMa4HnA9drrX8GPJnK0FrfTjJN9COS725/P6uktf4dcBkJwVyWyjaTRDkXKaXuIJk227dE/1eBvZRSd5NEOveWlDfhDpLprRuAD2mtH9Rafyu1+Xql1J0kn9Xd1dWITj7j+hbg00qp20k+aLWzq44LSqmLSM7By5RSG5VSp9RtSzAakLcGCwQCgSAYEpkIBAKBIBhCJgKBQCAIhpCJQCAQCIIhZCIQCASCYAiZCAQCgSAYQiYCgUAgCIaQiUAgEAiC8f8BmZm6rD5FTBEAAAAASUVORK5CYII=",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(mdb.nu_lines,\n",
" mdb.line_strength_ref,\n",
" \"+\",\n",
" color=\"black\",\n",
" label=\"activated lines\")\n",
"plt.legend()\n",
"plt.ylabel(\"line strength at 296 K\")\n",
"plt.xlabel(\"wavenumber cm-1\")\n",
"plt.yscale(\"log\")\n",
"plt.show()\n",
"\n",
"# %%"
]
},
{
"cell_type": "markdown",
"id": "932d4383",
"metadata": {},
"source": [
"Note that the above process is also applicable to MdbHitemp and MdbHitran."
]
},
{
"cell_type": "markdown",
"id": "6c26a201",
"metadata": {},
"source": [
"In tests/integration/moldb\n",
"these examples are available for reference.\n",
"\n",
"- quantum_states_filter_exomol.py\n",
"- quantum_states_filter_hitemp.py\n",
"- quantum_states_filter_hitran_co.py\n"
]
},
{
"cell_type": "markdown",
"id": "2e4178af",
"metadata": {},
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3.8.8 ('base')",
"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"
},
"vscode": {
"interpreter": {
"hash": "72bc7f8b1808a6f5ada3c6a20601509b8b1843160436d276d47f2ba819b3753b"
}
}
},
"nbformat": 4,
"nbformat_minor": 5
}