Skip to content

Instantly share code, notes, and snippets.

@HajimeKawahara
Created January 15, 2023 09:12
Show Gist options
  • Save HajimeKawahara/28cbe3bbe8f85060524502d9074446c3 to your computer and use it in GitHub Desktop.
Save HajimeKawahara/28cbe3bbe8f85060524502d9074446c3 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 96,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from astropy.constants import hbar, k_B, m_p\n",
"import astropy.units as u"
]
},
{
"cell_type": "code",
"execution_count": 30,
"metadata": {},
"outputs": [],
"source": [
"def morse_potential(r,params):\n",
" re,alpha,De = params\n",
" print(re,alpha,De)\n",
" return De*(1.0 - np.exp(-alpha*(r-re)) )**2"
]
},
{
"cell_type": "code",
"execution_count": 77,
"metadata": {},
"outputs": [],
"source": [
"mO = m_p*16\n",
"mC = m_p*12\n",
"m = mO*mC/(mO+mC)\n",
"convert_factor = ((hbar**2/(2*m*(1.0*u.Angstrom.to(u.m)*u.m)**2)/k_B).to(u.K)).value # (K)\n",
"\n",
"def centrifugal_potential(r,J):\n",
" \"\"\"\n",
" Args:\n",
" r (float): internuclear distance (AA)\n",
" J (int): rotational quantum number J\n",
" \"\"\"\n",
" return convert_factor*J*(J+1)/r**2"
]
},
{
"cell_type": "code",
"execution_count": 90,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"4.27 1.136 100.3\n"
]
}
],
"source": [
"r = np.linspace(3.0,10.0,1000)\n",
"params_CO = [4.27,1.136,100.3] # re [A], alpha[A-1], [D/kB] K \n",
"#ref: Matsumoto 1987 Parameters of the Morse Potential from Second Virial Coefficients of Gases\n",
"\n",
"Vco = morse_potential(r,params_CO)\n",
"Vrot10 = Vco+centrifugal_potential(r,10)\n",
"Vrot20 = Vco+centrifugal_potential(r,20)\n",
"Vrot30 = Vco+centrifugal_potential(r,30)\n",
"Vrot40 = Vco+centrifugal_potential(r,40)"
]
},
{
"cell_type": "code",
"execution_count": 105,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAEZCAYAAACq1zMoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOydd3xUVfr/3yd10gsJGZJAQhMJXUoAVzoooKLYsSvistZddV3XdcWyuip8xfWHIs2KIIIiqAhIsyAgvZcACQnpJCSkl3l+f0zmOqkkIWUC5/163dfMvefcc5+bzMznnnOe8zxKRNBoNBqNpqFxam4DNBqNRnNxogVGo9FoNI2CFhiNRqPRNApaYDQajUbTKGiB0Wg0Gk2joAVGo9FoNI1CswmMUupmpdQypVScUipfKXVEKfW6UsrHrk6kUkqq2fwrtGdSSr2llEoqa+83pdSQpr8zjUaj0QCo5loHo5TaApwCvgESgD7ANOAwMFhELEqpSOAk8DqwokITv4tIqV17C4HxwDPACeARYCwwSER2N+a9aDQajaYyzSkwwSKSVuHYPcDHwEgRWW8nMA+JyLwa2uoF7AYeEJEPy465AAeAIyJyfePchUaj0Wiqo9mGyCqKSxm/l72G1bG564Fi4Au79kuAxcDVSin3ehmp0Wg0mnrjaJP8Q8teD1U4/rpSqkQplaWUWqGU6lGhvBtwUkTyKhw/ALgBnRrBVo1Go9HUgEtzG2BDKRUGvAz8KCLbyw4XAh8Aa4A04HLgn8BmpdQAEbEJUSCQWUWzGXblVV1zCjAFwMvLq+/ll1/eELdSb1JSUkhISCAsLAwPDw/8/PyMsrNnz+Lk5ISvr2+t2iouLqa0tBSTydRY5rYIzp49i6urK15eXgCICEqpZrZKo7l42LFjR7qIBFdV1mxzMOWMUMob2AiEAgNEJKGGum2x9kxWiMhdZcfWAt4iMqhC3dFYxWmIiPxckw39+vWT7du311Sl0Zk/fz6TJ0/m1KlTtG3btlltuRiJiYlh+fLl3H333YSEhDS3ORrNRYFSaoeI9KuqrNmHyJRSJqweYh2Aq2sSFwARiQd+AfrbHc6g6l5KgF25w2PrncTGxhIXF1euLCsri/3799e6LRHh5MmTJCUlNaiNLRkvLy86dOhAUFAQAJmZmZSWlp7nLI1GU1+aVWCUUq7AMmAAME5E9tX2VMC+63UAaK+U8qxQLwooAmIu1NamwDYk9sILLzB+/PhyZbNmzaJHjx7k5+fXqi2lFF9//TVbtmxpcDtbKm3atGHixIk4OzsjIixatIhFixY1t1kazUVLcy60dAIWAiOBCSJSq19CpVQ74Epgq93hFYArcItdPRfgNmCNiBQ2lN2Nia0H4+7uztmzZ8uV3XDDDSxZsgQnp9r/y1q3bk1aWlXOehqAMWPGMHjwYABKS0vZt28fFoulma3SaC4emnOSfxZWQfgPkKuUGmhXliAiCUqpGVhF8Desk/xdgOcAC/CarbKI7FZKfQHMLOsVnQSmAu2BO5viZhoCm8C4uLhUEpioqCiioqLq1F5wcDDbt2/HYrHUSZguBZRSdOr0h3Ph4cOH+eqrr/Dw8Ch3XKPR1J/m/NUZW/b6PFYBsd8ml5UdAP6E1ZNsLdaV/r8C0SJypEJ79wMfAq8C3wFtgWtEZGfj3ULDYhsic3FxITc3l5KSEqOssLCQ3bt3k56eXuv2WrduTUlJSSWx0lQmKiqKu+++m44dOwKwd+9e9u3bhyM4wWg0LZXmXGgZKSKqmm1aWZ0FItJfRAJExEVEzCIyqQpxQUTyReRvZXVMIhItIhub+r4uhIouyFlZWcb7U6dO0adPH77//vtatxccbPUcTE1NbRgDL2KUUnTo0MFwYd6zZw+7du3SLs0azQXgMOtgNFYvJ6UUZrOZ2bNn4+7+RwCCtm3bsmzZMgYMGFDr9sxmM5MnT9YuufXgrrvuMhwq8vPzWbRoEaNGjaJdu3bNbJlG03LQAuNA2BZSmkwmHn744XJlJpOJiRMn1qk9FxcXwsLqGnVHA9Yejaen1SkxKyuLgoIC3NzcACgqKsLFxUXPa2k050F/QxwMX19fMjIy2LlzJ5mZ5YMT7N+/n71799apvfj4eDZt2tSQJl5ymM1mpk6ditlsBmD9+vV88MEHeg2NRnMetMA4GL6+viQmJtK3b182btxYruzBBx/kmWeeqVN7p06dYuPGjeTlVQzTpqkL9nMxERERREVF4ezsDFhFXLs3azSV0UNkDoafnx8FBQUAlXow77zzDh4eHnVqr02bNgAkJyfToUOHhjHyEqdr16507doVsMY6+/DDDxkyZAjDhg1rXsM0GgdD92AcDF9fX2NyuaJ78cCBA+nVq1ed2rMN6+iQMY2Dn58ft9xyC3379gWsHnu7d+/WPRqNBi0wDoevry85OTk4OzuTkVE+hFpCQgLffvttncb+PT098fPzIzk5uaFN1WAdOuvatSs+PtZM33v27GHVqlUUFraI4BEaTaOiBcbB8Pf3Jzs7m4CAgEoCs3z5cq677jrOnDlTpzbNZrNebNlEjBo1ioceesgYylyxYkWdgpRqNBcTeg7GwfD39+fs2bN8/vnnRERElCubOHEi/fv3x9/fv05tTpw4EVdX14Y0U1MNSikjWnNRUREpKSnGvohgsVgM5wCN5mJHC4yD4e/vT2FhIWPHjq00oR8aGkpoaGid27St39A0LW5ubkyePNkIN3P8+HG+/fZb7rzzTiPKgkZzMaOHyBwMW+9k27Zt/PLLL+XKCgoKWLFiBTExdcs+UFpaytdff83u3bsbzE5N7VBKGQsy3d3dCQ0NJTDQmrooIyOjXLw5jeZiQwuMg2ETmLfeeov77ruvXFlBQQETJkxgxYoVdWrT2dmZ+Ph4jh492lBmaupB27ZtufXWW418NIsXL9b5aDQXNXqIzMEICLAm4TSZTJUm+f38/Ni2bZsR8bcuhIeHExsb2xAmahqIa665xnhfWlrK7t276dmzp54v01w06B6Mg2Hrwbi5uXH27Nly6ymUUvTv398YYqkLYWFhnDt3juzs7AazVVN/bNGbbYtfjx07xrffflspVbZG05LRAuNg2ATGNoxiH7IfYNOmTXz33Xd1bjc8PBywrqXROB6XX345Dz74oNE73bFjBz///LPOR6Np0WiBcTBsAmOLfVUxXMxbb73FCy+8UOd2zWYzZrNZ/2A5MOHh4cb/PSEhgdjYWGNfB9bUtET0HIyDYROYNm3asGbNmkq5XN577716jdE7OztXSgGgcVwmTJhgeJjl5+fz/vvvM2rUKHr27NnMlmk0tUcLjINhMpkwmUyICKNHj65UfqEJr2w9GJ2p0fFxcbF+PUtKSoiMjKR169YA5ObmUlxcXOcFtxpNU6OHyBwQf39/0tLSWLp0KcePHy9XFhMTw3vvvce5c+fq3G5CQgJvvvkmp06daihTNU2Aj48PEydONAKX/vLLL8yaNcsIiqrROCpaYBwQm8DccsstrF27tlzZrl27eOSRR+rlctyqVSsKCgq0u3ILZ+DAgVx77bVGpIft27eTmprazFZpNJXRAuOABAQEGAnCKq6FGTt2LImJiURFRdW5XQ8PD8xms3aFbeH4+fkZaRsKCwtZt24dO3fubGarNJrKaIFxQPz9/Tl37hyenp6VBMbb25s2bdrUO2BiREQE8fHxOkTJRYK7uzuPP/44Q4YMAayJ5T777LNK3ocaTXOgBcYBsUVUDggIqPRDISLMmjWLDRs21Kvt9u3bU1JSwunTpxvCVI0D4OHhgaenJwBZWVlkZmZiMpkAq0OATn6maS60wDggNoEJDAys1INRSvHiiy/y5Zdf1qvtdu3aMXjwYLy9vRvCVI2D0aVLFx599FFjfuarr77i448/bmarNJcq2k3ZAfH39yczM5Mff/zRyJRoz5EjR+rtourh4VGl+7Pm4sHeBb1v377GcKiIsGPHDrp161YpFYRG0xhogXFAAgICKC0tpWPHjlX2NFq1anVB7ZeWlhIfH09oaKjOFXORY+8MkpSUxHfffYezszN9+vRpRqs0lwp6iMwBsfVO1q9fz/z58yuVr1+/ntdee63e7Z86dYqPP/6YkydP1rsNTcsjNDSUhx9+mB49egCwf/9+li9fTmFhYTNbprlY0QLjgNgE5quvvuKxxx6rVL5hwwZeeeWVesenateuHW5ubhw7duyC7NS0PMxmsxEh4Ny5c6Snpxu92OzsbB2rTtOgaIFxQOxD9ufn51NQUFCu/F//+he5ubn1dlV2dnamQ4cOxMTE6B+US5hBgwbx4IMPopSitLSU+fPn8+233za3WZqLiGYTGKXUzUqpZUqpOKVUvlLqiFLqdaWUT4V6AUqpeUqpdKVUrlLqR6VUjyraMyml3lJKJZW195tSakjT3VHDYUs6ZnvSrOiq7O7ubqThrS+dOnUiKyuLtLS0C2pH07KxdwgYNmyYMXxWWFjITz/9ZCz41WjqQ3P2YJ4GSoF/AtcA7wNTgbVKKScAZf30rygrfwy4CXAFNiilwiu0Nx94CPg3cC2QBKxWSvVu/FtpWGw9GJuIVHRVTktL47nnnrug1dudO3cGqBTrTHNpYpv4j4yMBODEiRNs2LCBs2fPAuierqZeNKcX2XUiYv/4vEkplQF8DAwD1gPXA38CRojIBgCl1G/ASeDvwONlx3oBk4AHROTDsmObgAPAy2XttBgquiBXFJji4mJmzJhB165dueKKK+p1DV9fXyZPnkybNm3qbafm4qVr1648/vjjRm963bp1pKamcvvtt19w71lz6dBsAlNBXGz8XvYaVvZ6PZBoE5ey87KUUiuBCZQJTFm9YuALu3olSqnFwD+UUu4i0mJcZfz8/ADrUFlMTIyRjdJGmzZtKCgouOAvelhY2PkraS5ZbOIC1ojORUVFxmfu+PHjtG3bVru5a2rE0R5Fhpa9Hip77Qbsr6LeAaCdUsrbrt5JEak4YHwAcAM6NbShjYmrqyve3t7k5eXRsWNH3N3dy5UrpRrkKbK0tJTVq1ezb9++C25Lc3ETHR3NuHHjAMjLy+Pzzz9n48aNzWuUxuFxGIFRSoVhHc76UUS2lx0OBKqK2mcbMwqoZb3Aaq45RSm1XSm13dEmu/39/cnIyGD69OlVfpFnzZrFW2+9dUHXcHZ2JiYmRkfi1dQJDw8P7r//fgYMGABASkoKn376KWfOnGlmyzSOhkMITFlP5BugBLjfvgioanaxYjrG2tYrh4jMEZF+ItIvODi4DhY3PrZ4ZNOmTavSdXTjxo2sW7fugq8TFRVFXFwcOTk5F9yW5tJAKUV4eLgxV5idnc3Zs2eN8DMZGRnk5uY2p4kaB6HZQ8UopUxYPcU6AENFJMGuOIOqex+2nkumXb2qcgkH2JW3KGyRlIOCgkhPT69UXt9glxXp1q0bP/30E4cOHaJ///4N0qbm0qJz58506tTJcHn+4YcfSE1N5YknntCpuS9xmrUHo5RyBZYBA4BxIlJxMuAA1vmVikQBp0Qkx65ee6WUZxX1ioCYhrO6aWjVqhUZGRnVCkxDERwcTOvWrdmzZ0+jXUNz8WMvJKNHj2bcuHEopRARvvrqKw4fPtyM1mmai+ZcaOkELARGAhNEZEsV1VYAYUqpoXbn+QLXlZXZ13MFbrGr5wLcBqxpSR5kNgIDAzlz5ky1ArN7927uuuuuC85OqZSiX79+BAUF1Tv0jEZjT3BwMJdddhkA+fn5pKamGkNmtlxEel3NpUFzDpHNwioI/wFylVID7coSyobKVgC/AZ8ppZ7BOiT2HNa5lTdtlUVkt1LqC2BmWa/oJNZFm+2BO5viZhoa+x5MTEzlDlhubi6//vorqampREREXNC1+vfvr4fHNI2Cp6cnDz/8sCEoR44cYenSpdx7773Gok7NxUtzCszYstfnyzZ7XgKmiYhFKXUtMB14DzBhFZzhIhJf4Zz7sYrVq4A/sAe4RkRapItUYGAgBQUFvPXWW8a6GHuuvPLKBo+GnJycTHBwcL1jnGk0VaGUMobQOnXqxIQJE4yHom3btpGYmMh1112nP3cXIc02RCYikSKiqtmm2dXLEJEHRCRQRDxFZKSIVJowEJF8EfmbiJhFxCQi0SKysSnvqSGx5XwpLS010uE2JidOnOCDDz7g6NGjjX4tzaWLu7s7vXv3NgSnoKCgXODW48ePa4/GiwiHcFPWVCYw0Oo89+uvv/LXv/61UrgYgCeffJL//ve/DXK9yMhI/Pz82Lp1a4O0p9HUhiFDhjBp0iTA+jC1dOlS1q5da5RbLJbmMk3TAGiBcVBsAnPo0CFmzpzJ6dOnK9U5ffo0KSkpDXI9JycnBgwYQFxcHElJSQ3SpkZTG2y9GWdnZx588EGGDLEGQc/OzmbGjBm6V92C0QLjoNiGyGxfvurWwrz99tsNds0rrrgCV1dXtmypyqFPo2l8goKCjM9+SUkJHTp0ICgoCLCmfN68ebPOwNmC0ALjoNh6MLYhgsZcC2PDZDLRp08fjh8/TklJSaNfT6OpicDAQG666SbjuxATE8PGjRuNh67MzEwtNg6OFhgHxf4pDqoWmLVr1zJ06FBSU1Mb7LpDhgzh0UcfNZKdaTSOwlVXXcUTTzxhRHD+9ttvWbBggVGu19Y4HlpgHBQPDw9MJhP5+fkopcjOzq5Ux8nJCRFpUK8bLy8vTCYTIqKfDjUOh5eXl/F+2LBhjBw5ErCKy7x587STioOhH1MdmFatWpGVlUVRUVGVPYqRI0caX7CGxGKxMG/ePMxmM9df36JytWkuIdq2bWu8LyoqolWrVoYAFRUVsWnTJvr161cur42madE9GAcmMDCQjIyMJh+ucnJyom3btuzevbvBvNQ0msbE3d2diRMn0r17d8DqYbllyxbOnTsHwLlz50hKStLDaE2MFhgHxhaP7K233uLVV1+tss748eN57bXXGvzaQ4cOxWQy8f333+svpabF0b59e5566imjl7Nz507mzJljDCcXFRXpz3UToAXGgbHFI9u0aRNff/11lXUCAwPLjUs3FJ6enowcOZJTp07pjJeaFomnp6fhcda/f39uu+02fHx8APjuu++YO3euFplGRguMA3O+iMoAn376KU888USjXL9Pnz6Ehoayd+/eRmlfo2kqPD09ufzyy439Tp060bNnT0OAvvzyS7Zt29Zc5l206El+B8bWg2nVqlWTrIOpiJOTE7fddluj9JA0muakR48exvvS0lJKSkqMJQEWi4Xvv/+ePn36EBYW1lwmXhToHowDExgYSFFREX5+fuTl5ZGXl1epzpdffkn79u2rjFXWEPj6+uLs7ExBQQHHjx9vlGtoNM2Js7Mzd9xxB4MHDwasKZ8PHDhAVlYWYE2NsXPnTgoKCprTzBaJFhgHxrbY0tvbm5CQEDIzMyvVCQkJYfDgwY3+4V+zZg2ff/45sbGxjXodjaa5CQoK4umnn6ZLly6ANYLAypUrOXv2LABZWVmkp6fr+ZtaoAXGgbGFyBg+fDjJyclVdteHDBnCwoULCQ0NbVRbxowZQ2BgIEuWLGm03pJG4yg4OzsbKQR69uzJn//8Z0JCQgDYunUr77//PkVFRYA1KGdxcXGz2erIaIFxYGw9GEf4QTeZTNxxxx0AfPLJJ8bTnEZzsaOUIiQkxHAIiI6O5pZbbsHd3R2A77//njlz5hj19VDaH2iBcWBsPZj4+HhuuOEGVqxYUWW9QYMG8Ze//KVJ7Ln77rspLCxk5cqVjX49jcYR8fPzK+eRFh0dzbBhw4z9+fPn88033xj7l3JOG+1F5sDYBObcuXN888039OvXr8rQLaNGjaJ9+/ZNYlObNm245557DM8yETGe7DSaSxH7756I0K9fPyPNeUlJCW+//TbDhw+nX79+xrzNpfKd0QLjwNgLjL+/f7VhW1555ZWmNIs2bdoA1i/T119/TXBwMFdeeSVOTrpDrLm0UUoRHR1t7BcVFdGtWzeCg4MB63D3hx9+yMSJE+nQocNF/4CmfxEcGJPJhJeXF+np6bRu3brGuGClpaVN7tVSUlKCxWJh/fr1LFiwgISEhCa9vkbj6Hh6ejJu3DgiIiIA60NZhw4d8Pf3B+Do0aPMmjXLIeZZGwPdg3FwgoODSUtLIyQkpFqBWbZsGXfccQcHDx6kU6dOTWabq6srN910E126dGHNmjXMnz+fqKgorrnmGiMkh6Nhe2IsKioiIyODgoICCgsLKS4uprS0lMjISPz8/MjMzOTgwYNYLBaUUiilcHJyIioqCj8/PzIyMjh16hRubm64urri5uaGm5sbQUFBuLq6XvRPppr6ERQUxMSJE419Nzc3AgMD8fX1BWDLli0cPHiQe+65BxcXlxb/OdIC4+DYBKZ79+7VJhbr2rUrf/vb3zCZTE1snXVIoEePHnTp0oXNmzeza9cuw7smISEBLy+vRg2Xbstbk5OTQ25uLn5+fvj7+5Odnc3GjRvJzc0lNzfXKL/uuuvo2bMnSUlJfPTRR5Xau/XWW/Hz8yM9PZ0ff/yxUrnZbMbPz4/ExMRyE7k2Jk+eTFhYGLt37+aHH37A09Oz3DZq1Ch8fHzIzMwkOzsbHx8ffH19dYK3S5T27duXm8Px8PDAz8/P+Dx89913ZGVlceeddwJWh4GWNBSt9GIhK/369ZPt27c3txmVGD9+PMnJyezYsaO5TakV9l+A2bNnk5KSgp+fH23btsVsNhMaGmp8oYqKinB1dS33hGZ7YhMRUlNTKSgoID8/n5ycHHJycggLC6Nz587k5OQwb948cnJyKC0tNc4fNWoUV155JZmZmSxYsAAvLy+8vb3x8vLC09OTHj16EBoaSn5+PrGxsZhMJtzd3XF1dcXFxQUvLy/c3NwoLS2ltLTUSOomIlgsFlxdXXF2dqa4uJicnByKi4spKioytsjISEwmE/Hx8Rw4cID8/HwjCkNeXh4PPPAAPj4+bNq0iY0bNxp2e3h44Ovry3333YfJZCIxMZGcnBwCAgLw9/fH1dW1af6BGodiy5Yt5ObmGnmfPv74Y/z9/ZkwYQJgHRq3rddpLpRSO0SkX1Vl+rHJwQkODq5VNGOLxUJOTo7R1W4u7J+ubrvtNo4cOUJ8fDynTp1i//79REVFGQLz9ttvU1hYaAiKiNC3b1+uvfZaRITZs2dXav/KK6+kc+fOmEwmIiMjywmIt7e3MZkaEBDAU089Va2dHh4edO3atdpy+4V2VeHq6lpjz6xt27blEmJVpE+fPoSHh3Pu3Dmys7PJzs4mJyfH6P1t376dXbt2GfV9fHxo1aoV99xzD0opkpOTcXJyIjAwUPd+LmIGDhxYbt/2mQfrw9jMmTO54oorGD58OAA5OTl4e3s3uZ3VoT+ZDo5tiGzVqlW88MILrFy50vDismfAgAGEhoZWu1amOQgICGDgwIHGl6SgoKDciuchQ4aQn59v9FpsC9rgj0Cbbm5ueHh4GCJiEzAXFxduuOGGpr+pBsLX17fGh4HRo0dzxRVXkJmZaWwlJSVGb+/HH3/k+PHjKKUICAggKCiI8PBwrrrqKgCKi4t1r+ciZOjQocb70tJSrrjiCsLDwwHIy8tjxowZXH311QwcOJDS0lJSU1MJCQlptmE1LTAOTnBwMAUFBWRlZbFjxw4SExOrFJhHH33U4aMem0ymcvNEgwYNqrG+/WK2Sw0PDw/Cw8ONH4+KjBkzhpSUFNLT043Nfrh77ty55Ofn07p1a2MLDw83enialo+Li4vRcwHrQ9k111xjjBAkJiayYMECbr31Vrp27Upubi7Jycm0a9euyR4+tMA4OLYfBDc3N4BqPcnuu+++pjJJ4wDYRMMee4G54oorSE5OJjU1le3bt1NSUkLPnj258cYbjfVLAQEBhISE0Lp1awIDA1vU5LGmMiaTqdwaHJvHms1F+tixY3zzzTdMnTqV1q1bk5ycTGJiIt27dzd+XxoaLTAOjk1gbF/+6gRGREhKSsLPz8/hezKaxsHeWcJ+7N5isZSLxF1QUEBiYiL79+83RMnFxYVRo0YRHR1NSUkJp0+fJiQkpFk8EzUNg4eHR7m8N127dsXHx8f4TTl8+DA///xzuToNjRYYB8f2YbB5SlUnMNu3b2fAgAEsX77c8DDRaMD6cGILnArWH55HH32U4uJi0tLSSE1NJTU11egRJScnGy7c/v7+mM1mQkJC6NmzpxFdQtPycHd3p2PHjsb+0KFD6d27d6MOlzWbwCilwoFngX5AL8ADaC8isRXqVedH3UdEdtvVcypr72HADBwBXhaRZQ1vfdNhE5hz584xfPjwasfQL7/8ct5991169uzZlOZpWjCurq6EhoZWSvUQHBzMpEmTSE5OJiUlheTkZA4fPkxkZCSBgYHExMTw66+/EhISgtlsxmw2Exwc3Ozuspq6oZQyIgo0Fs3Zg+kE3ArsAH4GxtRQ9yPggwrHjlbYfwV4Gni+rM3bgS+VUteKyPcNYXBzYBOUtLQ01q9fX209Hx8fHn300aYyS3MR4+7uTufOnencubNxrKioyBCQ0tJSiouL2bFjh5Fm2MnJiccee8yImZefn4/ZbNZDbJc4zSkwP4lICIBSajI1C8xpEdlSXaFSqjVWcfmviEwvO7xBKdUJ+C/QYgXG29sbd3d30tLSzls3IyODhIQE3YvRNDj2k8BdunShS5cuWCwWMjIyjJ6Oze1627Zt7Ny5E7AOsdl6OkOHDm3RYU80dafZBEZEGjJJwtWAG/BZheOfAQuUUu1F5GQDXq/JUEoZa2H++te/smfPnmp7Ms899xzLli0jPT29ia3UXIo4OTkRFBREUFAQ3bt3N46PGDGCrl27GsKTlJREWlqakTNl2bJlZGVlGXM7ZrOZ1q1b63U7FyEtZZJ/qlLqGaAU2AK8KCI/25V3AwqBmArnHSh7jQJapMDAH4stTSYTe/furbbelClTmDBhQosPkKdp2Xh5edGpU6dygVftk261atWK7Oxs9uzZY6QdbteuHffffz8AO3fuxNvbG7PZjI+Pj/4st2DqLDBKKXcgFOukfJqInH/s5sL4DPgWSAQigGeA9Uqp0SKysaxOIHBWKgdWy7Arr4RSagowBawfcFxJnBsAACAASURBVEfFJjADBw7kzJkzFBYWGiFF7Onbt28zWKfRnB/7NTbDhg1j2LBhiAhnz54lOTnZCHdjsVhYtWqVMbfj4eGB2WymV69e9OrVC3CM+Fua2lErgVFK+QB3AXcAAwBXQAGilEoEfgDmiMjvDW2giNxtt/uzUuobYD/wKvAnm4lAVd5mNT76iMgcYA5Yg11euLWNQ3BwMMeOHTO8fZKSkoiMjKxUz2KxsGPHDvz8/Ljsssua2EqNpm7YwtzYx3RzcnLiqaeeIiUlxfBgS0lJIS8vD7DG2nr77bcJDg4mJCTEGGILDQ3VDgUOyHkFRin1V+BfwAlgBfAfrL2JfKw9g+7AVcBapdQW4DEROdZYBovIOaXUd8CDdoczgACllKrQiwmwK2+x2HowNoFJTEysUmDA+nT40EMPMXPmzCa0UKNpOEwmExEREcYK9IoMGjSIlJQUTp48aQwZT5gwgd69e3PmzBl27NhhRDoICgpqtFXqmvNTmx7MYGCoiOyvpnwb1on0qcADwFCg0QSmjIo9lgOAO9CR8vMwUWWvBxvZnkYlODiYnJwc2rZty0033YSnp2eV9ZycnFi5cmWTJh3TaJoSb29vRo0aZezn5uaSkpJiLBJNT09n27Zt5VI4BAQEcOutt2I2m8nOzqagoIBWrVrpYbYm4LwCIyK31LItbxF57wLtOS9KKV9gPLDV7vAPQBFwJ/CS3fG7gP0t1YPMhm0tjK+vL0uXLq2x7ogRI5rCJI3GIfDy8qJDhw7GfpcuXfjnP/9JZmamEaEgLS3NCGG/d+9e1q1bZ0Q3aN26NcHBwQwaNEj3dBqB2s7BvCIiL9RQ3gpYh3VFfq1RSt1c9tY2Oz1WKZWG1Xlgk1LqaaALsIE/JvmfxrpS/05bOyKSqpR6G3hOKXUO2AncBowAWnzcFPvFlm3btq3RSywtLY1169Zx9dVXN2omSY3GUbGJR6tWrSrl/OnevTt+fn6G8CQmJnL48GH+9CfrdO7atWs5efKk4X4dFBREcHCwjkJdT2rrRfY3pVSKiPy/igVKqUCs4lKfdS1fVti39YA2AcOwhnu5sWzzA7KBX4EHRWRbhXOfB3KAJ/gjVMytIrKyHnY5FPYC07dvX7p168Ynn3xSZd39+/dzxx13sHr1asaMqWntqkZz6eHv718pPEpxcbExXBYQEEBKSgpxcXFGoj9fX1/++te/AvDzzz9TWFhYTnyq8ujUWKmtwNwGLFVKZYjI57aDSil/YC3gjFUQ6oSInM/LayVQK4EQkVKsnmWv1tUOR8c2vpyamorJZOL06dPV1o2OjmbPnj2XdC4VjaYu2C/w7NevH/36WbP/FhUVkZ6eTkFBgVGekJBATExMuXU9l112GXfccQdgfcDz9PQkODgYb2/vS34NT60ERkS+VUo9BMwvE5kflFJ+WMXFA6sTwJnGNPRSxmw2A9ZIyqGhoezfX52/BXh6eupQMRpNA+Dm5lYpEOgdd9xhpD9IS0sjPT29XArjFStWGFlb3dzcCAwMpFevXkb6hPj4eAIDA/H09LwkxKfWCy1F5NOy4bClSqlbgGmADzCsCRZbXtJ4e3vj6elJcnIyoaGhrFmzpsb6v/32G/v27WPKlClNZKFGc+lgP8dTkccee4z09HTS0tLIyMggIyPDEJL8/HwWLFgAWAOKBgYG0qpVK3r37k3Hjh2xWCwUFhbi4eHRpPfTmNRpJb+IvFM2of8tcBxrzyW5USzTGCilMJvNJCcn07NnT7Kzs8nJyTE8Yyry9ddf8+677/LAAw8YK6Q1Gk3jopTCx8cHHx8fI22xPa6urkyaNIkzZ84Y4pOQkGDUTUtLY/bs2Xh4eBjiExAQQPfu3QkKCsJisaCUalE9n9p6ka2ocKgYyAI+sL9ZEbm+4UzT2GM2m0lKSmLKlCk88cQTRiiNqnj22Wd54YUXtLhoNA6Ei4tLpTQI9nh6ejJ69GhDfGJjY9m7dy9hYWEEBQURExPDkiVL8Pf3NyIg2ATIx8fHIWMQ1vYXqOL8yqKGNkRTM2azmcOHDzNkyBCGDBlSY92quu4ajcax8fHxYfDgweWOlZSUGKLh5+dHdHQ0mZmZZGZmEh8fT2FhIR06dMDHx4cdO3awYcMGQ3hsQtS9e3fc3NyaRYBqO8l/f2MboqkZs9nMxo0bAWtO9dLSUmNysSpmz56Nl5cXd999d7V1NBqNY2M/ChESEsLo0aONfRGhoKDAcJMOCgqia9euZGZmcvr0aQ4ePIjFYiEqyhrQZMOGDezatQt/f3/8/PyM1759+5YLRtqg9jdKq5oGx2w2k5GRQXp6OsHBwbz55ps888wz1dZfuHAhgYGBWmA0mosUpVQ5h4DIyMhyMQotFgvZ2dlGENA2bdpw7tw5srKySExM5NChQ7i4uBhu2Y1BbYJd1jpZl7L2v8JFJP6CLdOUo02bNoDVE8Xb25vExMQa669evbramGUajebix8nJqdyi0q5du5aLbCAi5OXlNeqwWW36Rb8ppeYrpQZVV0EpFVAW7PIgF0FoFkfEthbG5qpc02JLQIuLRqOpEaVUjcPsDUFtBOZyrOHuv1NKpSmlflBKfaiUel8ptVgptRdIxRpY8smqwsloLhx7gQkPDychIaHG+rm5uUydOpUVKyo6AGo0Gk3TcF6BEZGzIvIMEAZMBQ4D/kB7oAT4GOgjIleKyOrGNPZSxl5g2rVrx6lTp2qs7+npyapVq4iJqZhFWqPRaJqGuqzkzweWlm2aJsYWjywpKYlbb73VSB9bHUopTp486XB+8RqN5tJBe5G1ENzc3GjVqhXJycmMHTuWsWPHnvccLS4ajaY5aRznZ02jYAsXU1JSQkxMDGfPnq2x/rFjxxg5ciSbN29uIgs1Go3mD7TAtCBsAnPixAk6d+7MypU1ZzIICgoiMzOTc+fONZGFmpZKcXExWVlZRhj6wsJCsrKyEBFj/9y5c8Z+aWlpuZD1Gk1VaIFpQdgEpm3btgDnnegPCAhg586dXH311U1hnqYJyc/P5+DBg+Tl5QFw5MgRXn31VVJSUgD45ZdfuOGGGwx39qVLlxIREUF8vHWJ2oIFCzCZTMZ6qtmzZ+Pv78+ZM9aoULNmzcLf3994OHn77bfx9fWlsLAQgP/85z84OztTWlpq7NuHtp8+fTr9+/c39j/44AMjZwrAl19+ycsvv2zsb9y4kS+++MLYP3ToELt27TL2bdErNC2LegtMWeh+TRNiExiTyURwcPB5BcaGiBhPnpqWQUZGBosWLTLc0Xfs2MEVV1zB9u3bAdi0aRPdunVj9+7dAMTExPDCCy8YApKXl8fJkycNAQoNDWX48OFGWJHu3bvz5JNPGuulrrrqKmbMmGGsixg6dCj/93//Z9QfMWIE06dPN0KXDB06lBdffNHIBNmjRw9uuukmw/7g4GAuu+wyYz8rK8sQP4CffvqJjz/+2NifP38+//znP439adOmlROkSZMm0bt3b2P/z3/+M5MmTTL233jjDd544w1jf+nSpXz77bfG/t69e8t5VOreVxNh+/GpywbMwRpNeR/WdTILsIbw/yvgVJ82m3vr27evODrTp08XQLKysqRfv35yzTXXnPeczZs3S1hYmPz+++9NYKGmtpSUlMiBAwckMTFRRETi4uKkb9++smLFChER2bdvnwCyePFiERE5evSojBs3zvg/Jicny+LFiyU1NVVERIqKiqSwsLAZ7qRhyMzMlISEBGN/37598vPPPxv7S5culXnz5hn7L7/8sjz77LPG/q233iq33HKLsV/x+9GrVy+5/vrrjf1u3brJ7bffbuzfdNNN8sorrxj7L774onz55ZfG/tdffy27d+829uPi4uTcuXP1uteLDWC7VKcV1RXUtAExgBswCOsiy78BY4BPgNfr02Zzby1BYD777DMB5PDhwzJx4kSJioo67zmJiYly++23y65du5rAQk11FBcXy6xZs2TDhg0iInLmzBkB5I033hARkZycHBkzZoysXr1aREQKCwvlwIEDkpeX11wmt2iysrIkPT3d2P/ll19k27Ztxv4777xjiLeIyJ133mn8L0REOnToII899pix7+vrK0888YSx7+HhIU899ZSIiFgsFmnbtq1Mnz5dRERKS0vlxhtvlK+++kpErP/7N998U7Zv327s//TTT5KcnGycX1xc3GD33tQ0hsDstHt/xO69c00Xc+StJQjM+vXrBZD169fLhg0b5JtvvmlukzQ1MG3aNJk9e7aIWH9EAgMD5dFHHzXKv/jiCzl58mQzWaepC4cPH5b4+HgRsf4vFyxYIFu2bBERa290ypQpsnz5chERyc3NlR49esjcuXNF5I+HiZkzZ4qItfcJyKxZs0REJD4+XgCZP3++sd+vXz9ZtWqViFgfEh977DGjB3XmzBn5/PPP5fTp0yIikpeXJydOnJCCgoKm+FNUojEEJh64Aetq/p0VynbUp83m3lqCwBw5ckQA+eSTT+p8bkZGhpSUlDSCVRobM2bMKDdsM2zYMJk8ebKxn5qaKhaLpTlM0zQjFotFcnJyjN5ofn6+rF27VuLi4kTEKhgvv/yyMcoQGxsrY8eOlY0bN4qIyJ49e8Tf31++/fZbEbEOewOGAG3atEkA+fHHH0VEZOPGjWI2m2Xr1q0iIvL777/LDTfcIMeOHRMRq1i++eabxvDqhfaeGkNgHi+bh/kNa5yybcCnwL+AuPq02dxbSxCYnJwcAeS1116T7Oxs2bBhQ7lhgOr44YcfxMnJyfjAaRqG5cuXy8MPP2zsT506tdw4f2lpaXOYpbnIyc/Pl0OHDkl2draIWHs4H374oSQlJYmIdf5qypQpcuLECRERWbdunfTo0UMOHz4sIiILFy40htpFxGinvjSYwADtqjkeCVwLPAcsrEubjrK1BIEREQkICJBHHnlEtm3bJoAxKVwTaWlp8sILL+jhmAskISFBpk+fbkymT58+XTp37mxM9ureiaYlYOtR2UY0LnRkoyaBqaub8ldKKfcqPNFigR9F5HURubOObWrqQFhYGAkJCbRr1w44/1oYsC64fPnll8slI9LUjrS0NLKzswHYuXMnTz/9NNu2bQPgiSee4MiRI3h7ewM6NI+mZWAL029zMbe9NgZ1FZgYrENj5VBKhQI/N4hFmhqxheoPDg7G3d2d2NjYWp1nsVjYvn27XtVfB06dOkVoaKixXmPMmDGcPHmSP/3pT4A1na0WFY2meuoqMA8AfZVSj9kOKKV6Y52DOd6QhmmqxiYwTk5OtG/fnhMnTtTqvN9++43+/fuzerXOqFATL730Eq+//joA7dq147///S9jxowBwN3dXfcCNZo6UKdoyiKSp5S6CfhVKbULCMI6uf+2iPy7MQzUlCc8PJzU1FSKioro2LEjx4/XTtejo6P55JNPGDlyZCNb2LIQEfbt20fPnj0Ba4gSWw5zgKeeeqq5TNNoWjznFRil1GpgN7Cr7PUIMAXryn0X4CERWdSYRmr+ICwsDBEhKSmJl156qdbnubi4cPfddzeiZS2T1157jWnTphEfH4/ZbGbhwoWNOiat0VxK1KYHswvoDdwDhAB5WEPElAJLgKNKKXcRKWw0KzUG4eHhACQkJHDllVfW6dyCggKWLVtGVFQUffr0aQzzHJ68vDzmzp3L6NGjiYqK4s477yQsLIyAgACgcSc8NZpLjfMKjIj8w/ZeKRUC9MEqOL2BIVjnZSxKqWMi0q22F1ZKhQPPAv2AXoAH0L7MI82+ngl4BbgLa6rm3cCzIvJThXpOZe09DJix9rReFpFltbWpJWAvMBkZGaxcuZLhw4cbXmXn4+GHH+bhhx++ZAUmPz+ff/7zn+Tk5BAVFUVkZCT33Xdfc5ulqSciQmFhIQUFBbi7u+Ph4UFxcTGnT5+muLiYkpISSkpKKC4uJjw8nNatW3Pu3Dm2bdtWrsxisdC3b18iIiJIT09n3bp1hqutxWJBRBgyZAjt2rXj9OnTrFmzplL5uHHjCA8P58SJE6xZswalVLntuuuuIyQkhGPHjvHrr7+WKwO47rrr8Pf35+jRo+zcubNS+bhx4/Dy8uLIkSMcPHjQKHNycsLJyYnRo0fj5uZGTEwMcXFxODk5lSsfPHgwTk5OxMXFkZaWZhz39/dvvLnF6vyXa7thFYZBwMN1PG8YkAJ8D6wGBIisot5C4CzwEDAS+ArIB3pXqPcfoBB4GhgOfABYgHG1saelrIPJyMgQQGbMmCGHDh0SQD799NNan3/06NFLbgHg3Llz5YEHHjD2bSuoNU1DUVGRpKSkGCvHRUS+++47WbRokcyZM0dmzJghL730kqxcuVJErHHYbrzxRhk7dqwMHz5cBg0aJH369JF33nlHRKzruvz9/cVkMknZ74axAFlE5MSJE+WO27Z3331XRKwr46sq/+ijj0TEGresqvJly5aJiMiqVauqLLfFkVuyZEmV5Zs3bxYRkQULFlRZvm/fPhGxxkmrqjw2NlZERF599dUqy8+cOSMiIs8++2yV5UVFRSIi8pe//KXc8RtvvPGC/r/UsA6mVpP8SqkZwHLgVxEpF+daRPKxruj/rTZt2fGTiISUtT8Za7DMitftBUwCHhCRD8uObQIOAC8D15cda41VWP4rItPLTt+glOoE/BeriF0U+Pv74+npSUJCApGRkSilaj3RD9C5c+dGtM5xKC0tNYa7UlNTiYuLIz8/Hw8Pj1r39jTVc/ToUeLj40lKSiI1NZWMjAwiIiJ46KGHAGt4/xMnTpCRkWG4xt911118+umnANx8883k5+eXa/ORRx7h2muvxdnZmWPHjuHh4YHJZMLb25ugoCBjGNPT05N77rkHk8lUbrMNGQcHB/Phhx/i4uKCq6srLi4uuLi4GI4cHTt2ZNOmTeXKnJycjDxLffr04cCBA8bTv62n0KZNG8CaquDkyZOVylu1agVYeyKJiYmVfmxbt25t3PuwYcPsH44REeP6d911F6NHj65Ubrv+5MmTGT9+fKUelK+vL2BNZTB+/HgsFouxiYjxfZg6dSrXXHONca7ZbG6gT0VllO0Gaqyk1PvAdYA78B1WsVldJi4XboRVYOZSYYhMKfUC8ALgLyJ5dsdfAv4B+IpIoVLqbqyRnC8TkWN29e7Hmkqgg4icrMmGfv36iS3XhqPTpUsXevXqxZIlS2jbti0jRowol1vjfPzvf/8jISGBN998sxGtbD6OHj3K9ddfz3vvvceIESPKiY2mekTEGI5ZtWoVe/fuJS4ujsTERJKSkoiIiGDJkiWANZ/MgQMHjHOdnJyYMGECX331FQAPPfQQRUVFBAYGGlu3bt0YNmwYYM1v4+HhgY+PDz4+Pnh7exu5ZjQtC6XUDhHpV1VZrf6jIjIVmKqUGgBMAF4FFiql1mEVm5UiktZQBtvRDThpLy5lHMCaLqBT2ftuWIfHYqqoBxAF1CgwLYmwsDAjU2FdXJVtHDt2jBMnTpT7QbkYKCwsxN3dnXbt2hEREYGTk3WZlxaXymzdupXNmzdz5MgRjh8/TlxcHC4uLhw8eBCAd955h9WrVxMYGEhYWBhms5n27dsb57/77rs4OTlhNpsJCQnB19fX+HsDzJ07t8br9+3bt3FurIGwWCwUFhYa8zu299UdKyoqMuZ8Ks79XMgxi8VipKe2f1/TsbrWt1gsvPHGG0ydOrXB/461HSJrJyKnRGQb1kWVz5cNP00A7gPeV0r9jlVsFonI6QayLxDIrOJ4hl257fWsVO6OVaxXDqXUFKwu1y1q2CQ8PJxNmzYB0KFDB1atWlWn8995551yPwYXA88//zw//PAD27Ztw2QyXfILSi0WCydOnGDPnj3s2bOHAwcOcPLkSbZv346TkxPz5s1j3rx5BAQE0KlTJ/r06UOnTp2M8z/88EO8vb3x8fGpsv3hw4c31a1UicViITc3l6ysLGPLzs4mNzfX2PLy8s773rZfUTxKSkoaxW4nJ6dKQ3e29xVfXVxccHZ2xsnJyXi1ne/m5lbuWFX1Kh6rqax79+6Ncr+17ZMeVkr9F3hDytyRRSQGmAHMUEoFYZ0Pub6s/vSqm6kzCutEVFXH61OvHCIyh7LQN/369WsxOYXDw8NJTEyktLSUl19+mVdffbVO59vEJTs7Gx8fnxbbiykpKTG+MH369KG4uJiioiI8PDya27QmpbS0lAMHDrBlyxZuv/12fH19efXVV3nxxRcBaw+uU6dOdOnShZycHHx9fY2IBa1atary/28b729MRIS8vDzS09M5c+YM6enp5bYzZ85w9uzZciJiLya1Gd4HcHNzw8vLCy8vLzw9PY33AQEBhIWF4enpiclkwt3dHXd393LvqztWcd/NzQ1XV9cahcM213MpUVuBeRKrq/D9Sqm/ishy+0IRScc617Ggge3LAKrqWgTYldteA5RSqkIvpmK9i4KIiAhKSkpISkoy3JbryubNmxkzZgzff/89Q4YMaWALG5+UlBRGjx7NE088wYMPPsjNN9/MzTff3NxmNRknTpxg3rx5bNmyhd9//52cnBwALr/8coYMGcKECRMICwujV69edOvWrZLohoaGNppthYWFJCcnk5SUZMzfJCYmGu9TUlIMESkoKKiyDaUUAQEB+Pv74+fnh5+fHx06dDDe22++vr7Gq7e3dzkR8fT01HM7zUht52DmKKUWAy8CXyilNgKPi8iRxjQO6xzKjUopzwrzMFFAEX/MuRzA6oDQkfLzMFFlrwcb2c4mxeazHhsbi5eXF++99x5jx47liiuuqHUbvXv35q677jI8W1oKFosFJycnWrduTY8ePRrVA8ZROHv2LJs2bWL9+vVMmDCBESNGkJ6ezltvvUXv3r257777GDhwINHR0XTs2BGAXr160atXrwa3RUQ4e/YssbGxxMXFERsba2xxcXHEx8dz5syZSuc5OztjNptp06YNYWFh9O7dm6CgoGo3f39/PXd2EVBraReRbOAppdQHwP8Be5VS7wIviUhjhehdAbwE3AJ8DKCUcgFuA9bIH9EDfsAqOHeW1bdxF7D/fB5kLQ17genevTv/+te/cHV1rZPAeHp6Mnv27EaysHFYtmwZ//73v/ntt9/w9fVl4cKFzW1So5GXl8crr7zCjz/+yM6dO7FYLHh4eNC5c2dGjBhB3759ycrKwtPTs1Gun5WVxdGjRzly5AhHjhwx3p88edJIX2DDy8uLyMhIIiMjiY6OJjQ01NjatGlDaGgoQUFBWjAuQercdxSRo8C1SqmrsQrNnUqpZ0Xkk7q2pZSyjWnYXErGKqXSgDQR2SQiu5VSXwAzlVKuWD3BpmJN1WzknRGRVKXU28BzSqlzwE6sIjQCqyPCRYXNISE2NhZ/f3/MZjOHDx+uV1vx8fGsW7fOoVez27zd2rZtS1hYGNnZ2YbP/8VCVlYWq1evJjc3l/vvvx8PDw8WL15M27Zt+fe//82IESMYMGAA7u7WdEzOzs4NIi45OTns27ePvXv3Gs4AR44cISUlxajj5OREhw4duOyyyxgyZAiRkZFEREQYohIYGNhi5/E0jUut1sFUe7JSzsAjwL+BoyIyuI7nV3fxTSIyrKyOB9ZV+pOwhorZgzVUzMYqbHkO64p/+1AxS2tjS0taBwPWSdjx48czb948hg8fTmFhIZs3b65zO3//+9959913iY+PJygoqBEsrT8Wi4XHH3+coKAgpk2b1tzmNDjZ2dksX76cxYsXs3btWkpKSujRowd79+4FoLi4GFdX1wa7XlpaGtu2bWPHjh3s2bOHvXv3cvz4cWOy3NfXl+7du9OlS5dyW8eOHXFzc2swOzQXFxe8DqZCY95Ad6BH2Wt3rHllouvaloic97GnbDHn38q2muqVYl2fUzeXqhZKZGSkkWysa9euLFq0qF7rWp599lkeeeQRhxMXsD455+TkXFReYfn5+ZhMJpRSPPPMM8yZM4eIiAiefPJJJkyYwKBBg4y6FyIuBQUF7N69my1btrB161a2bt3KyZPWkWKlFB07dqRXr17cc8899OzZk169ehEREaF7IpoGpbbrYF7nD0Fpi9X9Nx1rVOW9wGdYexaaJiIyMpLff/8dsHoO5eXlkZGRYYSrqC2tWrUyzsnNzcXLy6vBba0LcXFxPPnkk8ycOZOIiAg+/PDDFv+jJyLs3LmT+fPns3DhQtatW0e/fv148sknuf/++4mOjr7geywoKGDLli1s3LiRjRs3smXLFgoLrVOU4eHhREdH8+c//5no6Gj69u1rpHnWaBqT2vZgxmEVklllr3tFJKnRrNKcl8jISJYtW0ZpaSlTpkzhL3/5ywW5Y06bNo1Fixaxe/fuZu8xbN26lf3797f4J+rc3Fw++ugj5s2bx+7duzGZTNx8883G4sWuXbvWu22LxcKOHTtYtWoVGzZs4LfffqOwsNBYE/Too49y5ZVXGpPujkBxcTEJCQnVuiZrHBuTyUR4eHideta1dVNueH9HzQURGRlJcXHxBa2Fseeqq67i3LlzzbIQbN26dWzYsIFXX32ViIgITp48aUxmt0RsQTVFhOeff54OHTowa9YsJk2ahL+/f73bzczMZM2aNXz//fesWrWKtLQ0lFL06dOHRx55hGHDhnHVVVdd0DUak4SEBHx8fIwgrZqWg4hw5swZEhISyoUMOh96BVILxd5VOTw8nGeeeYbIyEgeeeSRerU3cuRII51yU8co++mnn1i8eDHPPPMMfn5+LVJcRIQNGzYwc+ZMTp48yd69e/H29ubQoUMXtCo+MTGRZcuWsXTpUn799VdKS0sJDAzkmmuuYdy4cVx99dUOOX9WFQUFBVpcWii2aNFpaXULOXlpxS24iLAXGID169ezYsWKC243NjaW6Ohodu7cecFtVUdubi7PPfccW7ZsAeAf//gHBw4cwM/Pr9Gu2ViICCtXriQ6OpqRI0eyZcsWJk6caMx/1EdcTp8+q1jN/AAAIABJREFUzTvvvMNVV11FeHg4jz/+OGfOnOEf//gHmzdvJjU1lYULF3LnnXe2GHGxocWl5VKf/53uwbRQ7NfCgHU8/+eff77gdt3d3SkuLiY3N/eC26oOEeHTTz/Fx8eHgQMHNvucz4WwfPlyJk6cSPv27ZkzZw533303JpOpzu3k5uby1Vdf8dFHH7F+/XoAevTowUsvvcTNN998QfM1Gk2zUV0mskttaykZLe0xm83y4IMPiojIK6+8IoDk5ORccLv2GS8bKvvjkiVLZMKECWKxWERE5OzZsw3SblNjsVhk5cqV8vnnn4uISHFxsSxatMjIFljXtn766Sd54IEHxNvbWwDp0KGDvPTSS3L48OGGNr3ZOXjwYHObIIDcddddxn5xcbEEBQXJ+PHjm8yGoUOHymWXXSY9e/aUwYMHn/d//fbbb0tubm69r7dhwwb59ddfjf33339fPv744xrPefHFF+Wtt96qdLyq/yE1ZLTUQ2QtmIprYYB6r+i3xzbRv2bNGjp27MiaNWvq1U5CQgLFxcWAdVFheno6GRnWuKMtcThs69atDBs2jOuuu453333XmhLWxYXbb7+9Tp412dnZ/L//9/+IiopiyJAhLFmyhFtvvZWffvqJmJgY/v3vf9OlS5dGvJNLFy8vL/bv329k01y7di1hYWF1aqMhQvkvXLiQPXv2cO+99/LMM8/UWHfmzJnk5VVMiVV7Nm7cWG4R9p///GfuueeeerdXF/QQWQumffv2xjxGjx49aNeuHZmZVaXPqR/R0dE8/fTTRrTlAwcOEB4eXqM4lJSU4OLiwq5du+jbty+fffYZkyZN4v777+eBBx5okWPwJ06c4O9//zvLli0jJCSE9957j8mTJ9f5Xg4dOsSsWbP4+OOPycnJYcCAAXz00UfcfPPNzb7+qKl58skn2b17d4O22bt3b2bOnHneemPHjuW7777j5ptvZtGiRdxxxx3G8HJGRgYPPPAAJ06cwNPTkzlz5tCzZ0+mTZtGYmIisbGxBAUF8fzzz3P//fdTVFSExWJh2bJldO7cmc8++4z//e9/FBUVER0dzXvvvVdjDLYhQ4YYNq9bt46nn36akpIS+vfvz/vvv88HH3xAYmIiw4cPJygoiA0bNrBmzRpefPFFCgsL6dixo5G7JzIyknvvvZeVK1dSXFzMl19+iclkYvbs2Tg7O/PZZ5/x7rvvsm7dOry9vXn66aeZO3cuc+bMoaioiE6dOvHpp582aHw73YNpwXTq1Im4uDiKioq47LLLiIuLY9SoUQ3Wvp+fH6+//roxp3DvvfcyYcIfod3ee+89Vq5cCVjXZXTp0oVXXnkFsEbzff311xk82Bo9yJa7vCUSGxvL6tWrmTZtGjExMUydOrVOPZaff/6Z8ePHExUVxdy5c5k4caKxuv7ee++95MSlubn99ttZvHgxBQUF7N27l+joP4KQvPjii/Tp04e9e/fy/9u77/ioyizg478TiggoXYpRCIL0KpFIWRIktFci8goSVmkqRZa2awNRQSkui4hBIrggQUpkF5CmL1UiiASICFIEpC5FVggQQxAhyXn/uJPZJEzIpEwmM/N8P5/5JHPrmTvJnHnuc+95Jk+enOGb/vfff8+qVatYsmQJs2fPZuTIkezdu5e4uDj8/f356aefWLp0Kdu3b2fv3r0UKVIk24Ksa9asoVGjRty4cYP+/fuzdOlS9u/fT3JyMh9//DEjRoygWrVqbNmyhS1btnDp0iUmTpxoL4LaokULpk+fbt9exYoV2bNnD0OHDmXatGnUqFGDIUOGMHr0aPbu3Uvbtm0z7L9Hjx7s3r2bffv2Ua9ePebNm5dPR9liWjAerHbt2vaRC+vWrevy/c2ePZtLly7Zn0+cOJHQ0FC6deuGn58fHTt2pHHjxoCVUF577TWXx+QKqsoXX3zB8ePHeeWVV2jfvj1nzpzJ0f0lqspXX33FlClT2L59OxUrVuSdd95h8ODBHjdEgis409JwlcaNG3Pq1Cmio6Pp2rVrhnnffvsty5cvB6B9+/bEx8eTkJAAQFhYmP2ClMcee4xJkyZx9uxZevToQe3atdm8eTPff/89gYGBgHU/VFbv9Z///GfuvvtuatSowcyZMzly5AgBAQE8/PDDgPVlbtasWYwaNSrDerGxsRw6dIjWrVsDcPPmzQzlhXr06AFYQ1KvWLEi22Nx4MABxo0bx9WrV7l27RqdOnXKdp2cMAnGg9WuXRuAn3/+mbp16xIZGcm8efOIi4tzSWuhRYuM9ezOnDmTYT8zZ87M930WtFOnTjF06FDWrVtH8+bNGTVqFMWKFXM6uagqK1euZPz48fz44488+OCDRERE8Pzzz7ustL6Rc2FhYbz88svExMRkGL9GHRT/TfsbT9/S7NOnDy1btuTLL7+kU6dOzJ07F1WlX79+TJkyJdv9L168OMP/k6MxdBxRVUJDQ4mOjnY4P321bWf6ivr378/KlStp0qQJUVFRxMTEOBWHs8wpMg+WPsGAdZpqz549nD9/vkD2nzaetzdITU1l1qxZNGzYkG+//ZYZM2awc+dOp0+FqSobN27k0Ucftd8Hs2DBAo4dO8bw4cNNcilkBg4cyFtvvUWjRo0yTP/Tn/5kP60VExNDxYoVHQ4NceLECWrWrMmIESMICwvjxx9/5PHHH2fZsmX8+uuvgNWfc/r0aafiqVu3LqdOneLYMWu8xIULF9KuXTsA7rnnHhITrSG3goKC2L59u32569evc/To0TtuO/36mSUmJlK1alVu3brlkvGVvOPTwUdVqFCB8uXL2xNM2giGaeXeDecdP36c0aNH06ZNGw4ePMjIkSOdru0WGxtL+/bt6dixI7/++ivz58/nwIED9O3bN1/L7Rv5x9/fn5EjR942ffz48cTFxdG4cWNef/11FixY4HD9pUuX0rBhQ5o2bcrhw4fp27cv9evXZ+LEifZTxaGhofzyi3MlG0uUKMH8+fPp2bMnjRo1ws/PjyFDhgAwaNAgunTpQkhICJUqVSIqKorw8HAaN25MUFBQtleOduvWjS+++IKmTZvedq/cu+++S8uWLQkNDXXNafasrl/2tYcn3gejqtqyZUtt3769qqpeuXJFAZ0yZYqbo/IMKSkpun79evvzffv22e/TccaZM2e0T58+Cuh9992nEREReuPGDVeE6hUKw30wRt6Y+2B8TO3ate0tmLJly1K9enX27TMjJ2Tnl19+oUuXLnTq1Mn+ra5x48ZO9V39/vvvvPvuu9SpU4fly5czbtw4jh8/zvDhwz2yjpphuIrp5Pdwadfep1Xw7d27t7nsNRtr165lwIABJCUl8fHHH9OmTRun1121ahUjRozgP//5D08//TRTp07NUXVZw/AlJsF4uLSO/uPHj9OwYUPee+89N0dUuL322mtMnTqVJk2aEB0d7XSNr/PnzzNixAiWL19Ow4YN2bJlC8HBwa4N1jA8nDlF5uEyX0kGkJKSkqfSEt6sadOmjB49mp07dzqVXFJTU5kzZw7169dn7dq1TJ48mT179pjkYhhOMAnGw2VOMElJSZQtW9atN7IVNuvWreOzzz4DIDw8nOnTpzvVV3LixAmCg4MZMmQIzZs3Z//+/YwZM8ZcGWYYTjIJxsOVKVOGSpUq2a+FL1WqFFWrVmX37t1ujsz9UlJSGD9+PF27dmXWrFmkpKQ4tZ6q8umnn9KkSRP27dvHvHnz2Lx5sz2ZG4bhHJNgvEDdunUzXAsfGBjo8wkmISGBbt26MWHCBJ577jm2bNlyx6KDaS5evEiPHj14/vnnadGiBfv37/fYIp3G7UqXLp2n9bdu3Urz5s0pWrQoy5YtyzBvwYIF1K5dm9q1a2d5/4yvMQnGCzRo0ICDBw/ay1wEBgZy7tw5p2/y8jbXr1+ndevWbNy4kcjISKKiopy6k37dunU0atSIr776imnTprF582b7wG6GAdZAf1FRUfTp0yfD9MuXLzNhwgR27tzJrl27mDBhQr5WNvdU5ioyL9CgQQOuXr3KL7/8QrVq1ezVYXfs2GEvfudLSpYsyXPPPcejjz5KSEhItssnJyczfvx4Jk2aRMOGDdmwYYO9aKfhOo4ulOjVqxcvvfQS169fv60QJVi1s/r378+lS5d4+umnM8zL7zpajqQNVZ65RNL69esJDQ2lfPnyAISGhrJu3TrCw8NdHlNhZhKMF6hfvz5gjddSrVo1HnnkEd5++237dF/x6aefUqdOHVq3bu10JecLFy4QHh5OTEwMAwcO5KOPPvLoIZyNnGvbtq3DWl3Tpk1zeviLc+fO8cADD9if+/v7c+7cuXyL0VOZBOMFGjRoAFgJJjQ0lOLFizN+/Hj3BlWAUlNTGTNmDFOnTqVPnz72UubZiYmJITw8nISEBObPn0///v1dG6iRwZ1aHCVLlrzj/IoVK+ZbiyVzfa7cSDs9nZ7ptzN9MF7hvvvuo0KFChw6dMg+LSkpiY0bN5KUlOTGyFzv5s2b9O3bl6lTpzJ06FCnOldVlQ8//JDHH3+cMmXKsHPnTpNcfFjbtm1p2rTpbY9NmzY5vQ1/f3/OnDljf3727FmqVavminA9imnBeAERsXf0p/nuu+/o2LEj69evp2PHjm6MznWSkpLo3r07mzZtYtKkSYwZMybbb403b95k2LBhzJ07lyeffJKFCxdyzz33FFDERmGUHy2YTp06MXbsWHvH/oYNG5waF8bbFfoWjIgEi4g6eFzNtFw5EZkrIpdEJElENolIo6y2620yX0n22GOPUaRIkXz55ymsSpQoQfny5Zk/fz5jx47NNrlcvHiRDh06MHfuXN544w1WrFhhkosPSU5OznMx0t27d+Pv78+///1vBg8ebD89Xb58ed58800CAwMJDAzkrbfesnf4+zJPasGMANLf3GEfrk2sT5bVQAAwHLgCjAG2iEhTVT1bkIG6Q4MGDUhISOD8+fPcf//9lC5dmmbNmhXIlTUF7ejRo5QsWRJ/f38+//xzp85179+/n27duvHf//6XJUuW+PzVPb7o4MGDPPTQQ3naRmBgIGfPOv44GThwIAMHDszT9r1NoW/BpPOTqsame8SlmxcGtAGeU9VoVV1nm+YHvOqOYAtaw4YNgYyDjXXo0IHY2NgsR7PzRHv27KF169b069cPcK4j9euvv6Z169bcunWLrVu3muTig2bPnk14eDgTJ050dyg+xZMSzJ2EAedVdUvaBFVNANYAT7otqgLUtGlTAH744Qf7tI4dO5KcnOw1rZjt27cTEhJCqVKlmDNnjlPrREdH07lzZ6pXr87OnTsJDAx0cZRGYTRkyBAOHTrktf2RhZUnJZjFIpIiIvEiskRE0t9i3QA44GCdg8CDIpK3+hAeoEyZMtSsWTNDgmndujW7du1yeMOap9m8eTMdO3akcuXKbNu2jVq1at1xeVXl/fffp0+fPrRq1Ypt27bh7+9fQNEahgGe0QeTALwPfAP8BjQDxgI7RKSZqv4KlAdOOVj3su1nOeCa60N1r2bNmmVIMMWLF/eKb+ypqam8/vrr1KxZk40bN1KlSpVsl//b3/7GjBkz6NmzJ5999hklSpQooGgNw0hT6FswqvqDqr6sqmtU9RtVnQF0BipjdfwDCHD7nU7W9CyJyCARiRORuIsXL+Zv4G7QrFkzjh8/TkJCgn3aiRMnGDZsGCdPnnRjZLmnqvj5+bFmzRpiYmKyTS7Jycn079+fGTNmMGrUKD7//HOTXAzDTQp9gnFEVfcAR4G0r+eXsVoxmZWz/XRYdU5VP1HVFqraolKlSvkfaAFr1qwZAPv27bNPS01NJTIyktWrV7srrFyLiorimWeeITk5mSpVqlChQoU7Ln/z5k169+7NwoULeffdd5k+ffptNaMMwyg4nvzfl77VchCrHyaz+sB/VNXrT4/B/xJM+tNktWrVolGjRixfvtxdYeVKZGQkAwYM4OrVq9y8eTPb5X///Xeeeuopli9fzgcffMC4ceNMqQ7jNnkt1z99+nTq169P48aNefzxxzl9+rR9ninX74CqetwDaAGkABNsz7tjJZt26Za5F4gHZjqzzUceeUS9QeXKlbVfv34Zpr399tsqInrhwgX3BJVDf//73xXQsLAw/f3337NdPjExUUNCQlREdM6cOQUQoZEbhw4dcncIWqpUqTyt//XXX2tSUpKqqkZGRmqvXr1UVTU+Pl4DAgI0Pj5eL1++rAEBAXr58uU8x1vYOHoPgTjN4nO10LdgRGSxiEwUkR4i0l5E/gasA84BM22LrQZ2AItEpLeIdLJNE2CqWwJ3kxYtWtw22FiPHj1QVVatWuWmqJz33nvv8dprr9G7d2+WLVuWbf9JQkICHTt2ZOvWrXz22WcMGjSogCI18io4OJioqCgAbt26RXBwMIsWLQKsMX2Cg4NZunQpYL3PwcHBrFixAoBLly4RHBzMmjVrAKsqdkEICQmxjy0UFBRkv+kyfbn+cuXK2cv1+zpPuIrsABCOdYd+SeACsAJ4W1UvAahqqog8AUwDIoESWAknRFXPONyqlwoKCuLLL7/k6tWrlC1bFoBGjRrRunVr/vjjDzdHl73g4GCGDx/OBx98kO0IlImJiXTp0oXdu3fzr3/9yyfHvjHyLrfl+ufNm0eXLl0AU64/K4U+wajqFCDbqnGqehkYaHv4rKCgIMCqmRQaGgpYd7t/++237gzrjlJTU9mwYQOdO3cmKCjI/hru5Nq1a3Tt2pVdu3aZ5OKh0t8AXKxYsQzPM5frL1OmTIbnmcv1Z3d14Z3kpl7fokWLiIuL45tvvgFMuf6sFPpTZEbOBAYGIiLExsbeNk9VOX/+vBuiylpycjL9+vWjS5cufPfdd06tk5SUxBNPPMGOHTuIjo42ycXIk5yW60+r3r169Wp78UxTrj8LWXXO+NrDWzr5VVUbNGigXbp0uW16eHi4Pvzww5qamuqGqG5348YNfeqppxTQSZMmObVOUlKStm/fXv38/HTJkiUujtDIT97Qyb9nzx6tWbOmHj16NMP0+Ph4rVGjhl6+fFkvX76sNWrU0Pj4+DztqzDyuk5+I+eCgoKIjY29rdneoUMHjh496nRLwZWuX79O9+7d+eKLL5gxYwZjx47Ndp0bN27QvXt3tmzZQlRUlClaaeRIfpTrf+WVV7h27Ro9e/akadOmhIWFAaZcf5ayyjy+9vCmFszcuXMV0MOHD2eYnpiYqGXKlNGePXu6KbL/Wb16tfr5+encuXOdWv7mzZv6xBNPKKCffvqpi6MzXMHdLZi9e/dqYGCgW2PwdKYFY9jHpN+6dWuG6aVLl2bw4MEsX77cbaVj1Naq6tatG4cPH+b555/Pdp3U1FT69+/P2rVr7TdgGkZOmHL97mESjBeqU6cOVapUYcuWLbfNGz58OH5+fsybN6/A4zp79iwtW7a0X3lTu3btbNdRVUaOHMmSJUuYNGkSQ4cOdXWYhhcy5frdo9BfpmzknIgQEhLCli1bUNUMl0v6+/uzdetWWrRoUaAxHTx4kM6dO/Pbb7+RkpLi9HoTJkzgo48+4q9//StjxoxxYYSGYeQ304LxUiEhIVy4cIEjR47cNu+xxx6jWLFi3Lp1q0Bi2bZtG23atCElJYWtW7fSvn17p9aLiIhgwoQJDBgwgGnTppn7CgzDw5gE46VCQkIAHJ4mA6t/pnr16hw6dMilcezfv5/Q0FCqVKnCjh07aNKkiVPrLVy4kJEjR9K9e3c++eQTk1wMwwOZBOOlHnroIfz9/fn6668dzq9Xrx43btxg0KBBOTpllVMNGjTg9ddf59tvv6V69epOrbN69WoGDBhA+/btiY6OpmhRcybXMDyRSTBeSkQIDQ1l06ZNDk+FVapUiYiICLZv305ERES+7jshIYF+/fpx+vRp/Pz8GD9+fLZjuaSJiYmhV69eNG/enJUrV5rBwox8lddy/bNnz6ZRo0Y0bdqUNm3aZDgDMGXKFGrVqkWdOnVYv359XkP1Dlldv+xrD2+6DybN8uXLFdCYmBiH81NTU7Vbt25avHhx3b59e77sc8eOHRoQEKBFixbVpUuX5mjduLg4veeee7RevXp68eLFfInHKDzcfR+Mat7v5E9ISLD/vmrVKu3UqZOqqh48eFAbN26sN27c0BMnTmjNmjU1OTk5T/sqjHJ6H4w59+DFQkNDKVasGF9++SXt2rW7bb6IMH/+fIKCgli5ciWtWrXK9b5u3brF1KlTefvtt/H39+ebb77J0faOHDlC586dKVeuHBs2bKBixYq5jsXwDGml+tNr0KABgYGB3Lp1i8WLF982P61O2PXr1/nXv/6VYV7//v1dFOn/3Hvvvfbfk5KS7H2Dq1atonfv3tx1110EBARQq1Ytdu3axWOPPebymAozk2C82D333ENwcDBr165l6lTHw+JUqFCB2NhYe1mLW7duUaxYsRzva+rUqYwbN45nnnmG2bNn24cKcMaZM2cIDQ3Fz8+PjRs34u/vn+P9G0Zu5bRc/6xZs5g+fTo3b96093GeO3cuQxVwU67fJqumja89vPEUmarqjBkzFNBjx45lu+yZM2e0Zs2aunDhQqcKYu7Zs0djY2NVVfXq1au6atWqHMd38eJFrVu3rt577726Z8+eHK9veA5vOEWW3uLFi7Vv376qqvrSSy/pwoUL7fMGDhyoy5Yty7d9FRamVIyRQVoxvmXLlmW7rKpSuXJlnnvuOVq1asWiRYsyjBSoqhw5coTIyEhatWpF8+bN7UUqy5QpY9+XsxITE+natSunTp1izZo1NGvWLEfrG0Z+yGm5/jS9e/dm5cqVgCnXn6WsMo+vPby1BaOq2rJlS23SpIlTyyYnJ+s///lPrVmzpgJapEgRvX79uqqqPv300woooPXr19epU6fqlStXchXTjRs3tEOHDlqkSJFctXwMz+MNLZj0ZfpXr16taZ8bBw4cyNDJHxAQYDr51XTy+4Tw8HBGjRrFTz/9RL169e64bJEiRXjhhRcYOHAgu3bt4siRI9x9990APPvss3Tq1Im2bdvy8MMP5/rmx5SUFJ599lk2bdrEggULctzyMYzcyI9y/R999BGbNm2iWLFilCtXjgULFgDWxQm9evWifv36FC1alFmzZmU75LdPyCrz+NrDm1sw58+fVz8/P33zzTfdHYqmpqbqoEGDFNDp06e7OxyjALm7BWPK9eed6YMxblO1alVCQkJYtGgRqampbo1l3LhxfPLJJ4wdO5bRo0e7NRbDd5hy/e5hEoyPGDRoECdPnmTdunVui2HKlClMnjyZQYMGmX90o0CZcv3uYRKMj+jevTtVqlQhMjLSLft///33GTt2LH369CEyMtIUrzQMH2ASjI8oXrw4L774Il999RUnTpwo0H3PnDmTl19+mZ49e7JgwQLT+WkYPsIkGB8yePBgihYtyj/+8Y8C2+ecOXMYMWIETz75JIsXLzaVkQ3Dh5gE40Puv/9+XnjhBebNm8fp06ddvr/IyEiGDBlC165dWbp0aa5K0BiG4blMgvExacMOT5482aX7ee+99xg2bBhhYWEsX748z/cfGEZ+yGu5/jTLli1DRIiLi7NPM+X6b2cSjI954IEHGDx4MHPnzmXv3r35vn1V5Y033mDMmDGEh4ezbNkyM6aL4VUSExOJiIigZcuW9mmHDh3i888/5+DBg6xbt46XXnrJpQP5eQqTYHzQO++8Q4UKFRg6dGi+3hfzxx9/MGDAACZPnsyLL77IwoULzWkxI0tRUVH2LzkpKSlERUXx448/AlZV76ioKA4cOADAjRs3iIqK4qeffgLg+vXrREVFceTIEQCuXbtWYHG/+eabvPrqqxm+OGVVrt/XmQTjg8qVK8f7779PbGwsH374Yb5s89KlS4SGhrJgwQImTJjAnDlzzNVihkfISbHLH374gTNnzvDEE09kmH7u3DkeeOAB+3NTrt8mq1v8PfEBPAAsAxKA34AVwIPOrOvNpWIcSU1N1e7du2vRokX1u+++y9O2du3apQEBAXrXXXdpdHR0PkVoeBt3l4pRzVuxy5SUFG3Xrp2ePHlSVVXbtWunu3fvVlVTrl+z+Fz1mmtGRaQk8DXwB9APq+rvRGCLiDRW1SR3xlfYpI1m+cgjj9CtWze2bduWbSHMzFJSUvjggw8YM2YMVatWJSYmJsOgS4bhCZwdcCwxMZEDBw4QHBwMwIULFwgLC2P16tWmXH9Wsso8nvYARgIpQK100wKAZOCv2a3vay2YND///LNWrlxZ77vvPt2xY4fT68XFxWmLFi0U0O7du2t8fLwLozS8gae3YDJL34Ix5fq9v9hlGBCrqsfSJqjqSWA78KTboirkatWqRUxMDKVLl6Zdu3ZMmDAhyw5TVSUuLo6ePXsSGBjI2bNniY6OZsWKFfYhlw2jsMqPcv1ZSV+uv3PnzqZcv41YCcjzicgFYJWqDs40PRLoqaqV7rR+ixYtNP017b7m0qVLjBgxgujoaMqWLUv37t1p2bIlZcuWJTExkf3797Nx40YOHz7Mvffey1/+8hdeffVVypQp4+7QDQ/hzHhErrRv3z5efPFFc3VXHjh6D0Xke1Vt4Wh5r+mDAcoDVxxMvwyUc7SCiAwCBgE8+OCDrovMA1SsWJElS5YwYsQIIiIiWLVqFVFRUfb5JUuWJCgoiJEjRxIeHm4Si+FRZs+eTUREBDNmzHB3KD7FmxIMWB37mWVZtldVPwE+AasF46qgPElQUBBBQUGkpqZy4cIFrl69SunSpbn//vtNk9/wWEOGDGHIkCHuDsPneFOCuYLVismsHI5bNsYd+Pn5Ua1aNXMljJGvVNUM1eChctOd4k2d/AeBBg6m1wcOFXAshmFkUqJECeLj43P1QWW4l6oSHx+f47JP3tSCWQ1ME5GaqnoCQERqAK2B190Yl2EYWHe3nz17losXL7o7FCMXSpQogb+/f47W8abmqBdNAAAKmElEQVQE80/gL8AqERmH1R/zLnAGmOPOwAzDgGLFihEQEODuMIwC5DWnyNS6U789cBRYCCwGTgLtVbXgKuEZhmEYgHe1YFDV/wD/191xGIZhGF7UgjEMwzAKF6+5kz+vROQi4PpxhAuHisAldwfhpcyxdR1zbF0nL8e2elaVUkyC8UEiEpdVaQcjb8yxdR1zbF3HVcfWnCIzDMMwXMIkGMMwDMMlTILxTZ+4OwAvZo6t65hj6zouObamD8YwDMNwCdOCMQzDMFzCJBjDMAzDJUyC8UEisk5EVEQmujsWbyEiXUVkq4hcE5HfRCRORNq7Oy5PJyKtRWSDiPxqO657RGSgu+PyJCLiLyIzRWSHiFy3/e/XcLBcCRH5h4j8IiK/25b/U172bRKMjxGRcKCJu+PwJiIyGFgFfA88BfQE/g2UdGdcnk5EGgObgGLAi1hloHYD80RkqDtj8zC1gF5Y42Jtu8Ny87CO81vAE8AvwHoRaZrbHZtOfh8iImWBw8BoYAkwSVXHuTcqz2b7JvgTMEZVzXi8+UhEJgMvA+XTF6wVkVhAVfUxtwXnQUTET1VTbb+/gFV5PkBVT6VbpgmwFxioqvNt04pijbN1RFXDcrNv04LxLVOBg6oa7e5AvMhAIBWY7e5AvFBx4Bbwe6bpVzGfXU5LSy7ZCMM61kvTrZcMfA50EpG7crNv8yb5CBFpA/QFXnJ3LF6mDVarsLeIHBeRZBE5JiLD3B2YF4iy/YwQkWoiUlZEXgQeBz5wX1heqQFwUlWvZ5p+ECvR18rNRr2qXL/hmIgUwxp0bZqqHnF3PF6mmu3xD2AscByrD+YjESmqqh+6MzhPpqoHRCQY+IL/fTG6BQxR1c/dFph3Ko/VR5PZ5XTzc8wkGN/wGnA3MMndgXghP+AeoL+qrrBN+9rWNzNGRCLUdHTmiojUBpZjfYsegnWq7ElgtojcUNXF7ozPywjWKMCOpueaSTBeTkQeBN4AXgDuynQu9S5bx3+iqqa4JUDPFw/UBjZmmr4B6AxUBc4XdFBeYjJWi+UJVb1lm7ZZRCoAH4pItJP9C0b2LgMPOpheLt38HDN9MN6vJlACWITVBE57gHWFzhWgkXtC8woHs5ie9s3PfADmXiNgX7rkkmYXUAG4r+BD8loHgQARyXxpfX3gJnAsNxs1Ccb77QVCHDzASjoh5PKPxwCs/gGATpmmdwLOquqFAo7Hm1wAmopI8UzTWwI3yOW3asOh1Vj3G/VMm2C7TPkZYIOq/pGbjZpTZF5OVa8CMZmniwjAaVW9bZ6RI18BW4A5IlIROAE8DXQEBrgzMC/wEdYNq2tEJBKrDyYMCAc+UNWb7gzOk4jI07ZfH7H97GIbxfeiqn6jqntFZCkww3ZR0ElgKBAA/DnX+zX9j75JRBRzo2W+EJF7gSlYiaUc1mXL76nqErcG5gVEpAvWRSoNsE71HscqLT/H9Bs6z/b/7sg3qhpsWybtQqA+QFlgH/BaXr6EmgRjGIZhuITpgzEMwzBcwiQYwzAMwyVMgjEMwzBcwiQYwzAMwyVMgjEMwzBcwiQYwzAMwyVMgjEMwzBcwiQYwzAMwyVMgjEMw6uISAcR+VlE9opIY3fH48vMnfyGYXgVETkIvII1VEJPVe3s5pB8lmnBGPlKRKJEZK2743C1wvY6RWStiESle16o4stvIlJORP4rIg85mH0Fq0L4CSAh03rLROSvBRGjYaopG/lvJDkYBU9EYoADqvoXl0Xkm5x+Hzz0PRgLfKWqxx3MiwR+AhKBupnmTQC+EZF5qppw25pGvjItGCNfqWqCbYiAAuVgzBCPl5fX5K73oSDYBsV6AZjnYJ4f8Ffg70BpIMMxUNX9WC2bZ10fqWESjJGv0p+aEZEYEYkUkckicklEfhWRabYPAWyndNoBw0REbY8atnkiIq+KyHER+V1E9ovIs+n2EyMiH9u2dxHY7sw+0y3zUVZxp9v/32ydxX+IyFkRmZLFa75jrLZlOovINhG5IiKXRWS9iNTLtIzD1+RgfyVt8V6znSYae6f3wfb8TyISa1snQUR2ikjDrN6DHMSb3bG+43F05tg50BVrpFBHx2cg1mmxqUARoKGDZVZjjSljuJhJMIar/RlIBloBfwFGYY2SB9ZpnB3AfKwO2arAGdu8icDzwDCsYVunYA3q9X/SbftZrNNAbYG+Tu7TWZOBN237bYA10t+ZLJZ1JtZSwAzgUSAY60NwjYNWSlavKb1pQCjwf4HHgWbAn7J6IWKNTLgK+BZogjUi5IdAClm/B87Gm92xzu44OnPsMmsLfK+ZrlASkTLAu8DLttbbeaCpg/V3AY+KNf6J4Uqqah7mkW8PIApYa/s9BtiRaf5GYG665zHAR5mWKYU1emHbTNNnYJ13T1vvRwf7z+0+08ddGmtI3iHZvU5nYs1i/VJYH/BtMsV122vKtF5p4A/gz5mmXQWisng95QEF2mWxzduORw7izfJYZ3cc83DsVgILHEx/H/gs3fMNjl4X0Nh2PB5y1/+JrzxMJ7/haj9men4euC+bdepjjV64TjKOxFcMOJXu+ff5uM/M+78L2OzkstnGKtbVTu9itR4qYZ098AMezLS9rF5TmoeA4litDgBU9ZqI7M9qBVW9bDsVtl5ENmO9rn+ralYtspzEe6djnd1xdPZ9zuxu4L+Z4q2D1S+T/pTYAawWWGa/p9uO4UImwRiudivTcyX7U7Np87sB/7nD9pJyuc9Ubr/Cqli6352+Cg7nY10DnAMG234mA4ewkkV6Wb2m3MRmp6oDRGQG0BlrXPtJItJdVddnsYqz8d7pWGcXq7PHLrNLWENTpzcDuBc4KWLfrQBJIiJqa7rYlLf9vJhNfEYemQRjuNtNrM7Y9A5hnQaqrqpfu2CfF7H6GtJrwv++Naft/3Hg52y2lW2sIlIBqAcMU9UttmnNyd3/3zGsD98grKuhEJFSWN/cHV2ya6eq+7DGWf+7iPw/oB+wnkzvQT7Gm91xzO37/APQP128T2C1tB7Bei1p6gL/xmr1HUs3vSFwXlUztIKM/GcSjOFup7A6XGsA14DLqpooItOAaWJ9Hd2KdT4/CEhV1U/yuM+vgRkiEgYcwfqW/oAtFmz7/xCYIiJ/2PZfAXhEVT9OvyEnY72C9a37RRE5A9wP/AOrVZAjttNh87CSxEWsU1JvcXuSthORANtrXI3VGqmJ1Q+R9lpOkek9yI94szuOeXif19tefwWse12mA++r6p5MrzutVdSUjAmmLbAuJ6/FyB2TYAx3mwYswPo2ezcQgPWB9ybWefaXsT4IfwP2Yl1+mlefYn3Afmp7Hgl8AVRMt8wYrMTwJuBvi+WzLLZ3x1hVNVVEngEisPoFjgF/A5bnMv6XsTrIvwCuAzNtz7NyHXgY69t8RVusi7HuFQHH70F+xZvdcczx+6yq+0VkF9Ab63WXt8WaebnfROQCVoJZBiAiJYCngE65eC1GDplaZIZheBwR6Yx1qXV9VU3JwXrDgCdVtaPLgjPszH0whmF4HFVdB8zCahXlxC1geP5HZDhiWjCGYRiGS5gWjGEYhuESJsEYhmEYLmESjGEYhuESJsEYhmEYLmESjGEYhuESJsEYhmEYLvH/Ad6XV5yuG4S+AAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig=plt.figure()\n",
"ax=fig.add_subplot(111)\n",
"ax.plot(r, Vco, color=\"black\",label=\"Morse Potential\")\n",
"ax.plot(r, Vrot10, color=\"black\", ls=\"dashed\", label=\"J = 10\")\n",
"ax.plot(r, Vrot20, color=\"black\", ls=\"dotted\", label=\"J = 20\")\n",
"ax.plot(r, Vrot30, color=\"gray\", ls=\"dashed\", label=\"J = 30\")\n",
"ax.plot(r, Vrot40, color=\"gray\", ls=\"dotted\", label=\"J = 40\")\n",
"plt.ylim(-10.0,250)\n",
"plt.ylabel(\"$V/k_B$ (K)\",fontsize=14)\n",
"plt.xlabel(\"internuclear distance ($\\AA$)\",fontsize=14)\n",
"plt.tick_params(labelsize=16)\n",
"plt.legend(fontsize=10)\n",
"#plt.savefig(\"morse_co.pdf\", bbox_inches=\"tight\", pad_inches=0.0)\n",
"plt.savefig(\"morse_co.png\", bbox_inches=\"tight\", pad_inches=0.0)\n",
"plt.show()"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.3"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment