{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from feos.si import *\n", "from feos.dft import *\n", "from feos.pcsaft import PcSaftParameters\n", "\n", "import matplotlib.pyplot as plt\n", "import numpy as np" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "params = PcSaftParameters.from_json([\"methanol\"], \"../../parameters/pcsaft/rehner2020.json\")\n", "func = HelmholtzEnergyFunctional.pcsaft(params)\n", "T = 350*KELVIN\n", "state = State(func, T, pressure=BAR)\n", "potential = ExternalPotential.LJ93(3.0, 100.0, 0.08)\n", "#potential = ExternalPotential.Steele(3.0, 100.0, 0.08)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3IAAAEvCAYAAAAAWPPhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABSRUlEQVR4nO3deZwcZ33v+89TVb1Pzz5aR7Is28jyAsaWDWYHAwaHJUDYTy4JIc4CZ0vuCRDuIZybS0JIOOTckOTGCSQkYQlZOCYJu4GwJMTIrJZt4V2akTT79L5V13P/qJ7RSBpZkmepXr7v16tetXR1969npPn1r56nnsdYaxEREREREZHO4UQdgIiIiIiIiFwYFXIiIiIiIiIdRoWciIiIiIhIh1EhJyIiIiIi0mFUyImIiIiIiHQYFXIiIiIiIiIdxos6gJVGR0ftnj17og5DREQ2wV133TVrrR2LOo5OoRwpItIbzjc/tlUht2fPHg4ePBh1GCIisgmMMY9GHUMnUY4UEekN55sf1bVSRERERESkw6iQExERERER6TAq5ERERERERDpMW90jt5pGo8HExATVajXqUB63ZDLJ+Pg4sVgs6lBERKSLKEeKiPSuti/kJiYmyGaz7NmzB2NM1OFcMGstc3NzTExMcPHFF0cdjoiIdBHlSBGR3tX2XSur1SojIyMdmaAAjDGMjIx09NVSERFpT8qRIiK9q+0LOaBjE9SSTo9fRETaV6fnmE6PX0QkKh1RyImIiIiIiMhJKuREREREREQ6TNsPdtIu9uzZQzabxXVdPM/j4MGDUYckIiLSFpQjRUQ2nwq5C/DVr36V0dHRqMMQERFpO8qRIiKbS10rRUREREREOsy6FHLGmEeMMT8yxnzfGHOwdWzYGPMlY8z9rfXQerxXVIwxvPCFL+S6667jtttuizocERHpEMqRIiKyEdaza+VzrbWzK/bfAdxhrX2fMeYdrf23r+UN/sc/HuKeY/m1vMQZrtjRz2+89Mpznvetb32LHTt2MD09zQte8AIuv/xynvWsZ61rLCIi0rWUI0VEZF1tZNfKlwMfbW1/FPjJDXyvDbdjxw4AtmzZwite8QruvPPOiCMSEZEOphwpIiJrsl4tchb4ojHGAn9irb0N2GqtPQ5grT1ujNmy1jc5n6uCG6FUKhEEAdlsllKpxBe/+EXe/e53RxKLiJyfO+6d4hvf/j4v/6f34jdqLDo+BaeJjTowIW0dBtgWdRibSTlSRCJ3dL7Mh77yALH6AlfP/788PHGUy35cjzosWYP1KuSebq091kpEXzLG3He+TzTG3ArcCrB79+51Cmd9TU1N8YpXvAIA3/d5wxvewIte9KKIoxKRx/Khrz7A1ff+KclDM8z3QdyFMcCokotcDEg0J6MOYzMpR4pI5P72rgn+5uBRfrrvII/2H2Tv99JcNgH5TNSRyeO1LoWctfZYaz1tjPk0cAMwZYzZ3rrSuB2YPstzbwNuAzhw4EBbfsXau3cvP/jBD6IOQ0QuwNH5MtclTwAw8e43cu0zX8ulg5dijIk4MlnWI78L5UgRaQdzxRqjfXF+86YxfuMuh/6Gx+DTbuSJf6oBitrOeebHNd8jZ4zJGGOyS9vAC4G7gc8Ab2qd9ibg9rW+l4jI+SjWfGaLdeL1GgCvvP7nuGzoMhVxsumUI0WkXRRrPpmEB8UpCo5Lum5ws31RhyVrsB4tcluBT7e+IHnAx621nzfGfAf4lDHm54AjwKvX4b1ERM7pyFwZgFitAYCT7Y8yHOltypEi0hZKNZ9M3IPSLKVYglTV4vRlow5L1mDNhZy19iHgSascnwNuWuvri4hcqImFsJBzaz6BASeTjjgi6VXKkSLSLgpVn2zSg1qeouuRrDZw1CLX0TZy+gERkUjkKmFLnFcP8BNGXSpFRKTnnSzkCpQD8HyLm1WLXCdTISciXadY8wGLVwM/uV6D84qIiHSuYs0nm4xBrUDQmnVAXSs7m77hiEjXKVZ9hswCyRrYVCzqcERERCJXrPn0JcIWuWYjHARXg510NhVyItJ1inWfHfFpUjWwmVTU4YiIiETKWkuh2qAv6dGsF3HqYac8R10rO5q6VopI1ylWfbYn5kjXLKZPVxtFRKS31fyARtOSTXoU60XStbBFzlGO7Ggq5M7Tm9/8ZrZs2cJVV111yvHPf/7z7Nu3j0svvZT3ve99EUUnIisVaz4jsXnSNXAHBqMOR6TrKUeKtLfw3nHIxh1Kfol0OM2qBjvpcCrkztPP/MzP8PnPf/6UY81mk7e+9a187nOf45577uETn/gE99xzT0QRisiSYtUn6y2QroE3OBp1OCJdTzlSpL0VqmEhN+g1KBizXMipa2VnUyF3np71rGcxPDx8yrE777yTSy+9lL179xKPx3nd617H7bffHlGEIrKkUPNJmUXSNYj1D0YdjkjXU44UaW/FpULOqVJynJOFXCYTYVSyVirk1mBycpJdu3Yt74+PjzM5ORlhRCICYcKKN/N4ASQHhs/9BBFZd8qRIu2jUA3nV826NYrOyRY5V/fIdbTOGrXyc++AEz9a39fcdjW8+PH127fWnnFMEw+LRK9Y8/GaJQCSgyMRRyOySZQjReQsCkv3yFHhmOOEg50kE5iYpujpZJ1VyLWZ8fFxjh49urw/MTHBjh07IoxIRCAs5KwNrz6mVMiJREI5UqR9LHWtzFKm2OpaaXR/XMfrrELucV4V3CjXX389999/Pw8//DA7d+7kk5/8JB//+MejDkukpy3NlRMYH/Bw+/ujDklkcyhHishZLHWtTNnictdKT/mx4+keufP0+te/nhtvvJHDhw8zPj7Ohz/8YTzP40Mf+hA333wz+/fv5zWveQ1XXnll1KGK9LSaH+A3m9h6AIDTpyuOIhtNOVKkvS1NP5Bqlig6Dn1V8LIq5DpdZ7XIRegTn/jEqsdvueUWbrnllk2ORkTOJl9tkKFKsxFep3L7VciJbDTlSJH2Vqj5JDwHr56naBz66g6O8mPHU4uciHSVQtWnnzK2Hg6q4OiKo4iI9LhC1Seb9KCWp+g69NUMrvJjx1MhJyJdpVD1yZoytAo5N6uhlUVEpLcVqz59CQ+qeQpenFTVqkWuC6iQE5GuUqg2yFLGrTsEjsGkUlGHJCIiEqlizSebjEE1R8nxSFUDXN1D3vFUyIlIV8lXfPpNGbcOjbSneatERKTnFaqNsEWulqdmHdymWuS6gQo5EekqhWqDPkrEa4ZmOhF1OCIiIpErVH36kh5UczRqS7ce6B65TqdCTkS6SqHqk3ILpGtgM+mowxEREYlc2LUyvEcuqFsAHN1D3vFUyIlIVylUGyTdPOmqBSUpERGRcCCwhEdQzUE1nGfV1YTgHU+FnIh0lXzVJx2vkKmBk1X/fxER6W3WWoq1sGtluZ4jVQtb5FzlyI6nQu48HD16lOc+97ns37+fK6+8kv/1v/7X8mOf//zn2bdvH5deeinve9/7IoxSRCC86hj3KqRrutooshmUI0XaW7UR0Aws2YRLoVEmXQuPO8qRHc+LOoBO4HkeH/jAB7j22mspFApcd911vOAFL2Dfvn289a1v5Utf+hLj4+Ncf/31vOxlL+OKK66IOmSRnpWvNvDcCpkqxAYGow5HpOspR4q0t0K1AcCQVydvIFMNjzt9uv2g06mQOw/bt29n+/btAGSzWfbv38/k5CS5XI5LL72UvXv3AvC6172O22+/XUlKJEKFagNjKyQbEB8Yijocka6nHCnS3go1H4BBp0LBdZZb5NRrpfOpa+UFeuSRR/je977HU57yFCYnJ9m1a9fyY+Pj40xOTkYYnYgUqj5OM8xSiYHhiKMR6S3KkSLtp1hdKuTKFI1DpmqxMQ+T0BQ9na6jWuR+587f4b75+9b1NS8fvpy33/D28zq3WCzyqle9it///d+nv78fa+0Z52jyYZFoFao+NOoApIZGI45GZPMoR4rIavKtrpVZKky7hnQNTF9G/x+7QEcVclFqNBq86lWv4o1vfCOvfOUrgfDq4tGjR5fPmZiYYMeOHVGFKCKEXSut7wNGXStFNolypEj7ylfCFrl+UybvOGSqmgy8W6xbIWeMcYGDwKS19iXGmGHgb4A9wCPAa6y1C2t5j/O9KrjerLX83M/9HPv37+dXfuVXlo9ff/313H///Tz88MPs3LmTT37yk3z84x+PJEYRCf+vFqo+QbMJeLqRW9qGcqRypEhUcpWwRa7Plig4Dv018AYGIo5K1sN63iP3n4F7V+y/A7jDWnsZcEdrvyN961vf4q/+6q/4yle+wjXXXMM111zDZz/7WTzP40Mf+hA333wz+/fv5zWveQ1XXnll1OGK9KxKowlBA9tozZGjG7mlfShHKkeKRGKpa2U6KFJwHPrqBk8tcl1hXVrkjDHjwE8A7wWWLse9HHhOa/ujwNeAaC4XrtEznvGMVfv6A9xyyy3ccsstmxyRiKymUPXpo0LQCPv9a7JTaQfKkcqRIlHKVRrEXEPMbxVyNQdH+bErrFeL3O8DvwYEK45ttdYeB2itt6z2RGPMrcaYg8aYgzMzM+sUjoj0onylQdaUoR7+aVOikjbx+yhHikhEcpUGA6kYppan4Hpkqha3X/mxG6y5kDPGvASYttbe9Xieb629zVp7wFp7YGxsbK3hiEgPy1Ua9FPGqRms0WSnEj3lSBGJWr7SoD8Zg2qOghcjVbU46lrZFdaja+XTgZcZY24BkkC/MeavgSljzHZr7XFjzHZgeh3eS0TkrHKVBv2mjNswNJIextFUmRI55UgRiVSu0qA/FYNqnhIOsUYTN6sLnd1gzd9yrLXvtNaOW2v3AK8DvmKt/Q/AZ4A3tU57E3D7Wt9LROSxLJYb9FHCqxmaqVjU4YgoR4pI5PLLhVwOv7F064Fa5LrBRl6ufh/wAmPM/cALWvsiIhtmsdIg6xRI1SDIpKIOR+SxKEeKyKbIV30GUjGo5Qlq4THdI9cd1nVCcGvt1whH3sJaOwfctJ6vLyLyWHKVBmk3T7pmsdlM1OGInEI5UkSiEA524mFn57G1cMwlDQbWHXQDiYh0jVy5TjZeJFMFoyQlIiI9zlq7PNhJqbpIstqaZ1U5siuokBORrpGrNEjHSvRVwR0YiDocERGRSJXrTfzAMpD0yNcK9FXD48qR3UGFnIh0jcVKg7hbpq8C3uBQ1OGIiIhEKldpADAcq5NzLH2V8LijQq4rrOs9ciIiUcpVGjhBmWQD3KGRqMMRERGJVL4aFnKjbomc45BtFXLe4GB0Qcm66ahC7sRv/Ra1e+9b19dM7L+cbb/+6+c873Of+xzvf//7ecYznsHznvc83vOe93DgwAE+8IEPrGs8IvL45coNvEaYpRLDoxFHI7K5lCNF5HS5cljIDZkSU45DX8VCOoWJxyOOTNaDulaepzvuuIOvfvWrZLNZPvjBD/LlL3+ZJz3pSdxzzz1RhyYiLblKA2rh2Mrp4S0RRyPSO5QjRdpTvuoD0E+RnOuSrYDTrznkukVHtcidz1XBjWJtOMrP8PAwxWKRZrOJtXb5uIhEy1rLYqWBafqAR2pEhZz0FuVIETnd0j1y/bZI3nHCe8iHdA95t1CL3Hl6/vOfz0033cS9997Le97zHl74whdy8OBBrrzyyqhDExGgVG9C4BPUmoASlchmUo4UaU9LhVy6mSfnOPTXDDHlx67RUS1yUXrxi1/Mi1/84uX9r3/96xFGIyKnWyzX6adEUAuvT7m6kVtk0yhHirSnfKuQS/p5cq5Df8UoP3YRtciJSFdYLDcYNCVMvVXIaWhlERHpcblKg2zSw6kukPNi9FU19UA3USEnIl0hX2kwSBGnavA9g0mlog5JREQkUvlKg/5kDCqL5NwYqXJTLXJdRIWciHSFxUqDAVPCqxlqaQ9jTNQhiYiIRCpfbTCQikFlgZrv4ljNIddNOqKQ6/RRrzo9fpFOkKs0GKBIrGpo9CWiDkdk03R6jun0+EXaWa7SoD/lQWWBoBpe4FTXyu7R9oVcMplkbm6uY//QW2uZm5sjmUxGHYpIV1u6Ry5ZBT+rbpXSG5QjReSxLJYbDKbiUFnAVsNRndW1snu0/aiV4+PjTExMMDMzE3Uoj1symWR8fDzqMES6Wq7SYNApkK4COzTZqfQG5UgReSwL5QZDmRjVqQUSlfBrv7pWdo+2L+RisRgXX3xx1GGISJvLVersSRToq4BVkpIeoRwpImdjrWWxXGcwFSNXz9FXGQHUtbKbtH3XShGR85GrNIjHCmQrEBsajjocERGRSBVrPn5gGUsG5KxPthIeV9fK7qFCTkS6wmK5gRMU8QJIDY9FHY6IiEikFkrhZOBjXoWc69BXsVhjcPt1+0G3UCEnIl0hV2ng+OHlxvTw1oijERERidZCuQ7AqFtm0XHoq4Lpy2BcN+LIZL2okBORrrBYbhDUawBkx3ZEHI2IiEi0lgq5YbfIguuSrYAzqPvjuokKORHpCrlKg2Yt7EaSGVGLnIiI9LalQm7AFllwHfp0D3nXUSEnIh2v5jep1Go0a+EcOd7QUMQRiYiIRGvpHrlskGfRcRioGuXHLqNCTkQ63kKpwRBF/FrY79/V0MoiItLjFst1jIFUY4F516W/6ig/dhkVciLS8eZKNYZMAVs3gAo5ERGRhXKDgVQMpzzHohcnUwk09UCXUSEnIh1vvlRnxORxqg71hIuJx6MOSUREJFLz5TrD6TiUZ8kZj2RVhVy3USEnIh1vvlRnmDzxCjT6k1GHIyIiErnFcp3BdAxKszRq4Vd+b3gk4qhkPamQE5GON1es0+8skqkYmoPZqMMRERGJ3EKpwVCrRY6KBcAd0aiV3WTNhZwxJmmMudMY8wNjzCFjzP9oHR82xnzJGHN/a61hckRkQ8yX6qS9BQZKFjOiPzXSPpQjRSQqYYtcnEp5nlQpAMAbGY04KllP69EiVwOeZ619EnAN8CJjzFOBdwB3WGsvA+5o7YuIrLu5Up2+RJ6BkpKUtB3lSBGJRHiPnMdCdZ7+UnjMU4tcV1lzIWdDxdZurLVY4OXAR1vHPwr85FrfS0RkNfOlGjGvQH8FEqNbog5HZJlypIhEodpoUm0EjCV8FvAZKIfH3RHdI9dN1uUeOWOMa4z5PjANfMla++/AVmvtcYDWWt+uRGRDzJfqxPwSjoX02PaowxE5hXKkiGy2hXIdgK1ekUXXZaBssfEYTiYTcWSyntalkLPWNq211wDjwA3GmKvO97nGmFuNMQeNMQdnZmbWIxwR6TFzpTpOrQJAdtt4xNGInEo5UkQ223wpLOS2OAXmHYeBEjhDgxhjIo5M1tO6jlpprV0Evga8CJgyxmwHaK2nz/Kc26y1B6y1B8bGxtYzHBHpEfOlOrbaACClFjlpU8qRIrJZFsthThwyeRZdl/4SeOpW2XXWY9TKMWPMYGs7BTwfuA/4DPCm1mlvAm5f63uJiJzObwYslusElSagG7mlvShHikgUlrpWDgR55l2HwTLER3QxqNt46/Aa24GPGmNcwsLwU9bafzLG/BvwKWPMzwFHgFevw3uJiJxiodwgSwWqYXcR3cgtbUY5UkQ23VwxLOT6gxxzrsuBisEb1ajO3WbNhZy19ofAk1c5PgfctNbXFxF5LPOlOsMmj6k6BA64AwNRhySyTDlSRKIwV6xhDKQai8y5MbIlqx4rXWhd75ETEdlsc6UaI+TxKoZaJoZx9GdNRER622ypznA6jlOZo2hjxHyLO6weK91G33hEpKMtlBoMmjyJCvgDqajDERERidxcscZIXxxKszQq4dd9tch1HxVyItLR5ks1Btx5smUIBrNRhyMiIhK5uWKdkUwCW54hqFoAtch1IRVyItLR5kp1+ty5cI6cEd3ILSIiMlusMZpNUCjPki6FhZxa5LqPCjkR6WjzpTrZRI6BMsRGt0QdjoiISOTminVG0jHmqgsMlsJjri52dh0VciLS0eZKdZImT6oOybGtUYcjIiISqWqjSaHmsz3VYM7W6S+Hx72hwUjjkvWnQk5EOtp8sU6yUQAgs2VHxNGIiIhEa74UziG30ysw5zgMlCxk+zDxeMSRyXpTISciHW26UMWpVgDo37or4mhERESitTQZ+FYnnAx8oAyu7o/rSirkRKSjTRdqUPEBSIzpHjkREelts8UaACMsMue6DJYgPjIWcVSyEVTIiUjHqjaaNKtFKLdG5Nqqe+RERKS3LRVyg8ECc67LSNEQU37sSirkRKRjTedrjJlFnIpLYMAb0Rw5IiLS2+Za98j1NWaZc12GClYXOruUCjkR6VgzxSpjLOKVDJVsDON5UYckIiISqblijWTMIVaZoRQkiPkWb4tuPehGKuREpGOFLXI5UiVDbTAVdTgiIiKRmyvWGckkMMVp6pXwq763RffIdSMVciLSsaYLNUbNAtkiNEcGog5HREQkcjPFGqN9cYLiCUwxANA9cl1KhZyIdKzpQpUBb46hIqCrjSIiIswV64z2JZgvTjNQDI+pa2V3UiEnIh1rOl9jILZIfwViW3S1UUREZLZYYyTtMtPIMaxCrqupkBORjjVTrJFt5ABIbNsRcTQiIiLRagaW2WKNPakKM45huGCx2QxOMhl1aLIBVMiJSMeaztdwq+Hlxv6dF0UcjYiISLTmSjUCC7vjBaY9l6EiuLr1oGupkBORjjVdqBFUqgBs2bUv4mhERESiNZ0PJwPf5uSYac0hl9i6PeKoZKOokBORjuQ3A3KlEs1yE4D0jl0RRyQiIhKtqXx4cXPMLDLtuowWHeJbt0UclWwUFXIi0pHmS3XG7CLNioPvgjs4GHVIIiIikZouhC1yg8E8s45HfynQQCddTIWciHSk6UKNbWYet+xQ7o9jjIk6JBERkUgttcj1Vaco+THcQJOBdzMVciLSkaYLVbaZORIlQ2M4E3U4IiIikZvK1xjJxHGLx/ErHqDJwLuZCjkR6UjT+Roj7jSDRWBsNOpwREREIjdTqDKWTdDITeAUA0BzyHUzFXIi0pGmCzVS3hRDRYhrDjkRERGm8jW29ieZLh1nSJOBdz0VciLSkY7nqvQ5C6TqkNEcciIiIkwXquzsg+N+kdG8xToO3qh6rXQrFXIi0pGO5yqkagUABi+6NOJoREREotUMLDOFGnvjeU64LmM5MFtHMZ4XdWiyQdZcyBljdhljvmqMudcYc8gY859bx4eNMV8yxtzfWg+tPVwRkdCxxQpuKRyda3DPEyKORmR1ypEislnmijUCC7tjC5zwPEZzlsSO8ajDkg20Hi1yPvCr1tr9wFOBtxpjrgDeAdxhrb0MuKO1LyKyLqYWy5iiD0BiXJOBS9tSjhSRTbE0h9w2M88Jz2Vr3pBUfuxqay7krLXHrbXfbW0XgHuBncDLgY+2Tvso8JNrfS8REYB8tUG8Ngclh0bM4A6pMUPak3KkiGyWpTnkRpuzTBuPwYIltlODgXWzdb1HzhizB3gy8O/AVmvtcQgTGaAhc0RkXRxfrLLNLJAsOJSHEpoMXDqCcqSIbKSlFrmBxjTlagzHQmyHCrlutm6FnDGmD/h74L9Ya/MX8LxbjTEHjTEHZ2Zm1iscEelix3IVtpoZsgXwxwaiDkfknJQjRWSjHc9VMQZS1SmCkguokOt261LIGWNihAnqY9baf2gdnjLGbG89vh2YXu251trbrLUHrLUHxsbG1iMcEelyxxYrDMcmGcuBs2N71OGIPCblSBHZDMcWK2zNJqkWJknnLaBCrtutx6iVBvgwcK+19n+ueOgzwJta228Cbl/re4mIQNi1csCcoL8CyV17ow5H5KyUI0VksxzPVdg+mORE6QSjrXZ/b7sudnaz9WiRezrw08DzjDHfby23AO8DXmCMuR94QWtfRGTNjuUqDPtzAPTv0Rxy0taUI0VkUxxbrHJRv8OxRp6xnMWODOIkElGHJRtozTMEWmu/CZxtpIGb1vr6IiKnO7ZY4apKCYDRPfsjjkbk7JQjRWQzWGs5tljhtXvKTHoeYznw1K2y663rqJUiIpvheK6KWwpH58rs3hNtMCIiIhGbL9Wp+QF7Y3NMeh5b8pAevyjqsGSDqZATkY4SBJZibg5ThKYDngaAEBGRHndsMZxDbpxpJjyPkbwlrjnkup4KORHpKHOlOlubU8SKDqXBOMZ1ow5JREQkUsdyFQDG/CnytRheU10re4EKORHpKMdzFcbNNNmcoTaWjTocERGRyB1bDAu5/toxmsVwCIz4rl1RhiSbQIWciHSUiYUKW9zjbF0Au0tXG0VERI7nqiQ8h3rxCNnFcA65+EW6R67bqZATkY7y6FyZYTtJXxXiF2sOORERkcnFCjsGUxwrTLJtwWJdR5OB9wAVciLSUY7MlxmpzgLQf8nlEUcjIiISveOLFS7OBkw2S2yfB3Zsw3hrnmVM2pwKORHpKEfmS6RLrTnkLrs64mhERESid2yxyhXpRSY8j20LluSePVGHJJtAhZyIdJQjs0VihTqBgdFLr4o6HBERkUjV/YDpQpVLY3Mc8Ty2LUB6zyVRhyWbQIWciHSMRjOgmTuGycPigIObSEQdkoiISKSOLVYILFzkzDBb90g2IL5HA530AhVyItIxJhcqXGROkMo5lLdkog5HREQkckfmywBsD6aoF2MAxC/aE2FEsllUyIlIx3h0vswuM8HwooHx7VGHIyIiErmlQi5VfZh4Lvxqrxa53qDhbESkYxyZLzNmHyZTg9pejVgpIiJydL5M3HWYLjwaTj3gucS262JnL1CLnIh0jCNzJYaqUwCM7NOIlSIiIkfmy+wZ8jhSnWbbPLBzG8Z1ow5LNoEKORHpGEfmyyQrRQB2XH4g4mhERESid3ShzHX9OY54LtsXLKk9e6MOSTaJCjkR6RhHZwu4eR/fgayGVhYREeHIXJmrkjMccWNsn4f0pZdFHZJsEhVyItIRrLU0Fo6SXHDIjcUxsVjUIYmIiEQqV26Qr/pc6kyxWPKINSFxmQq5XqFCTkQ6wlS+xhZ/guE5Q33XSNThiIiIRG556oHmBM1cHFAh10tUyIlIR3hgushOc5gtOYhd9oSowxEREYncUiFnao8yMgvWQGKv7pHrFSrkRKQjPDBdYGflEQCGrro+2mBERETawFIhd6xylF2zlmDHFpxUKuKoZLOokBORjvDgTInB0jwA49c8PeJoREREovfIbIndmSYP+wV2zViSGuikp6iQE5GO8MBUgVSuTt2DgT1KVCIiIg/NFnn6wDwPt0as7L/8qqhDkk2kQk5EOsLizFFSC4aFLXFNdCoiIkLYW+W61HFyhRiu1UAnvUaFnIi0vVylwXD5gXDEynGNWCkiIrJYrjNfqrPPOYqf8wAVcr3GizoAEZFzeWC6yMUcZrgIjSdoxEoREZEHZ0oAJBsPMTRvsA7EL94TbVCyqdQiJyJt78HpIjuLDwEw9KSnRhyNiIhI9B6aKQJwovYIe05AcNFOnHg84qhkM6mQE5G298BMkeHcAgC7rn9uxNGIiIhE76HZEiNuifubRfaesPRd/aSoQ5JNti6FnDHmI8aYaWPM3SuODRtjvmSMub+1HlqP9xKR3vPgVI7sfIO5AcPwtouiDkfkvCk/ishGeXC6yDP6ZzlajzNYhv4nPjnqkGSTrVeL3F8ALzrt2DuAO6y1lwF3tPZFRC5Y6diP6Z8xzI2now5F5EL9BcqPIrIBHpotcX1mitp8DIDUVVdGHJFstnUp5Ky1XwfmTzv8cuCjre2PAj+5Hu8lIr1loVRn58IPGFk0NC4djzockQui/CgiG8FvBjw6V+IS9wiZWZfAMSQuvzzqsGSTbeQ9cluttccBWustG/heItKl7j2e5wmVHwEwfOCZEUcjsi6UH0VkTR6dL9NoWprNw1x8AvyLtuMkk1GHJZss8sFOjDG3GmMOGmMOzszMRB2OiLSZe47nGcqfAOCyG18ccTQim0s5UkRWc/hEAYeAo/Vj7D1hyVz1xKhDkghsZCE3ZYzZDtBaT692krX2NmvtAWvtgbGxsQ0MR0Q60T2TOTILNWYHDeM790cdjsh6OK/8CMqRIrK6+04UuMw5xsM1w0AZhq+5PuqQJAIbWch9BnhTa/tNwO0b+F4i0qVmjz3IwLRDfjyDMSbqcETWg/KjiKzJ4RN5npM9Rnk+nDdOA530pvWafuATwL8B+4wxE8aYnwPeB7zAGHM/8ILWvojIeas2mgwf+wYjeXCv2hd1OCIXTPlRRDbC4RMFrkk+zNBxl2bcI7lfPVZ6kbceL2Ktff1ZHrppPV5fRHrTA9NF9uW/D8D257w02mBEHgflRxFZb+W6z6PzZcz2B9g3YWnu34uJx6MOSyIQ+WAnIiJnc+hYjq1zs1TisO/Gl0QdjoiISOTunyqCDThePcHFUzB4w9OiDkkiokJORNrWjx6dZXiqyfHxGMlEJupwREREInf4RIG95jgzcw6OhZGnPCPqkCQiKuREpG1VDn+DsVlDfd+OqEMRERFpC/eeyPOU2I8x0zECA6lrrok6JImICjkRaUvFms/Oo3fgAMNPe07U4YiIiLSFQ8fyXNl/mD2ThtolO3D71GOlV6mQE5G29MOJRcbnj+I7sO+5ZxsvQkREpHcEgeXQZI6G8zCXHbP0Xaf543qZCjkRaUvfe3SBbccqPDzusG30oqjDERERidxDsyXi9QVy0xUSPmx/9s1RhyQRUiEnIm1p4YffYmzWkNs/GnUoIiIibeFHk4tc6/wYMxmn6RoyT7kh6pAkQirkRKTtWGvZ8cPPAJB+5rMijkZERKQ9/HAix6Xpu3nCI1C9Yg9ORvfH9TIVciLSdh6dK7PjxIMsZuDK57wh6nBERETawt2TOfrMA1w0A6PPe2HU4UjEVMiJSNv5t/tPsOVYjR9f7LFv5PKowxEREYlcM7A8PHkCe7wEwLbn3RJxRBI1FXIi0nbmvvJ3pKoG50l7MMZEHY6IiEjkHpgucnXwA/qPulQGEiSecFnUIUnEVMiJSFux1jLyg38Kpx142c9EHY6IiEhbOPjoPJfE7+RJD4F92nW60Ckq5ESkvfx4qsCeR2a59yLDtVe/POpwRERE2sJdjywwNn+EZAP2vOKNUYcjbUCFnIi0lUOf/wKDeVi8agTP8aIOR0REpC08+PD9ZB4NKKcdhp+mEZ1FhZyItJsv/wVNAxf/xCujjkRERKQtTOWr7Kx8kSsegtoN+zCeLnSKCjkRaSO1RpPt9z/CPRfB0258S9ThiIiItIWDjyywb/EuUnW4+DU/H3U40iZUyIlI2/jeP36BwVzAxBV9ZBLZqMMRERFpC99/6Bg7HqxQShu2PfP5UYcjbUKFnIi0jdInP0QlDttfrLlxREREllTu/gSXPwTFZ+7DxGJRhyNtQoWciLQFv1Bg7J6H+df9hptu/MWowxEREWkLs8Uau+//Kl4A+97y9qjDkTaiQk5E2sL9f/0xYj4cfdIAW/q2Rx2OiIhIW/jGfUd5wo/LTIx7bL/6qVGHI21EhZyIRM5aS+6TH+bIGDzrOT8ddTgiIiJtY/Izv8XWBUi9UEWcnEqFnIhErvStf2Vgqsg3rnN48VN/IepwRERE2oK1lj3/die5DNzwS++POhxpMyrkRCRy93/wt5jvg/Q11+A4btThiIiItIXvf/lzXHK0yYMH+olnh6IOR9qMCjkRiVTl0CGShx7iiwcM/+Hm90QdjoiISNuY/P/eSzkO+17/S1GHIm1IhZyIROrY7/0OpQTk9g+ye/tlUYcjIiLSFhYP/YCLD81z6GrDgWe/KepwpA2pkBORyJS+/W3q//Yd/uHpDtdfoauNIiIiSw6951cpJ8E+41lgTNThSBtSIScikbBBwPHfeR9z/fDgFTFe9Zw3RB2SiIhIW8j967cY/tEk/34g4NoXvDPqcKRNqZATkUgs/u3f0bj3MJ94luGGodeRjGmQExEREVuv8+C7385MP9T37ubKvbujDkna1IYXcsaYFxljDhtjHjDGvGOj309E2l9jaooT7/8d7r4Imrvj/NTNb4s6JJFNp/woIquZ/OMPkZqY45vPaWIu/k84jrpVyuo2tJAzxrjAHwIvBq4AXm+MuWIj31NE2pu1luO/8R78WoW/vtkwEnsLF432RR2WyKZSfhSR1VQPHyZ324f55hWGsewOnnnj06MOSdrYRrfI3QA8YK19yFpbBz4JvHyD31NE2tj8X3yU0te+xl8/23ClHeClN78+6pBEoqD8KCKnaBaLPPy2XyaXCjh+Y4V/dd/Itbs1d5ycnbfBr78TOLpifwJ4ysoTjDG3ArcC7N26lcKXvwyOi3GdVdbhYlx39XMcA66LcRxw3RXnOsvHVl2LyKYo33UX07/3e/zwCS6HntRkZ/1XeMrFw1GHJRKFc+ZHWCVH3nHHeeXBU/LfheTD049ppDyRTWGt5di73kUweYwPv9ZwS2kvyac+U90q5TFtdCG32r8+e8qOtbcBtwFclUzZibf9xw0O6UyBY8AYrHNywXHOWLO0dh3simSI64LngediPC9Mlp4XbnsejudhvBiOF8PEYjieh+PFcTwPNxbH8eK48TiuF8OJxcBdem7rdVwPEzvzdXFdjBcLH4vFMPH4inUcE4+tWMeUkCVStQcf5Ogvv5XFIZffv8Xnxumreemrn61/l9KrzpkfYZUc+dbNvZ/UGtPKdScvpp6tSDSuizGtY0vFYKyVw5Zy44r9MGfGwuedvh/zVs+Fq+TG5VzYysFn7HseJh7HWZkjT8mXsfA8kQjNfOADFL/wRT72XIfXphf5rflf4U+u2Rl1WNLmNvov1wSwa8X+OHDsbCcXhwO+8eoKgTVgwbbWgQ3zXRCE59nALKc7GwDWEEB4bMXzWDrPAhisNZgVjxt78jlYC9ZignB/ae2sXILT1q1tNwCnBl4FnMDiNcPjXnCWdTN87innnZG+11/gudhYmEjxPIjHIBY7JZk58QROIoEbby2JZJj8lh5PJDCJJE6ytU4lMYkkJpnASSYxidY6mWztt85NJsPX0Jf2ntSYnOTIW36eCnX++081eWUlxneG38ZznjAWdWgiUbmg/AhwYovL+24dxAYBNgjCpNi02KAJQbheLUctrc258tkZx21rHeDYIHz+WfKg08qpS9ueNXjWwfUNXt3gWfACgxsYvFbeXF6aS7nU4gQWp7li3bQ4zWBDfxFAWIyeXtzFz7wg6pxxofS05yznwpP5cTlftnLhco5MtHLj0mMqJnvW3Ic/zNyffZgvXuvQuLJKMfdshnbtZ9+2bNShSZvb6L8a3wEuM8ZcDEwCrwPOOlnU6Ogl3HrrX0DgQ9A8ubZL20vHz3Is8E87Hpx2bMV5y8fPdszHNn0C6xMETRpBgyDwaQYNGs0GzcDHtw2Cpk8j8PGDBk3bDI8HPg0DPtA0Bh+Db6BmDA3AN+F+eDzcblhD0xqaFpqBoUm49jHYwITHWkVo0zpY6xBYQ4ATFrOBwTYNBBbjAwGYpsVthgWj51u8ICDmB8SajfBYa4k1wauDV4WYb8P90x6PN0247YeF6uNhjYFEHJOILxd8TjKFm0rjZtKYdBpnaclkWtuZM49lzjxPCbB91e6/nyNv+XnqpTz//dV1bohV+ezEr/E/f+lqFfbSyy4oPwLs8St8NH/PWR+3BpoGAgd849B0PJpeLFw7Hk3HpenGwvXSfmvxHZfAhOum4xI4Dr5xaRpD4Lj4jkOztd80zsm147bWhoZx8DHUTZjbwnxnwzXQIMAHfBvQwIZr28QnCLeDZrhvw7zaCBr4gY/fbOA3G+D7NP0wf7mnF4MrisJTj4V5MLYy3y1t+5C0LonAJRm4xAOHROAQbzaJBZZ4s0EsgFjTEGuAW13KgRbXB7cZ4PoBjh/gNJqYRhPTaGCCx3llduli6Yoiz0kmMankufPiyty48lg6jUml9Le2TVlrmf2DP2D2j/6Y71wZ52s3WT4y5/Cswiv4jRddFHV40gE29NuvtdY3xrwN+ALgAh+x1h466xPiadh1w0aGdEEMYdAuELvQJwdNaDagWQ8LxGY93A8areOnbS+ft8r20nP9arjtV8FvrZs18Fcup54T+BXqzQb1Zo16s0GjWaMWNMKi0hhqxlBzDHUMdQM1JzyWX3rMGKomTMxVY6i5HjU8fOvgBy5B4NBs1cg0wPoW2whwfUvch7gPiUa4jjcscb/eWorEW8cTOUjNGVINh1QdEnVLom6JNS7gKmw8jpNO42azuNksTn8/brYPJ9sf7mezuP1ZnL7WOtt6vL/1eF9f2O1H1lXha1/j2K+9nWbM4d2v9cn21ykdew1XP/F6nqwbuKWHXXB+BBjeC69+P+EfXT/MIYG/vG+CBl5rP77K46vvt3LV8gVQH3w/zC0rjzVb5y6/RmNFbquv/w/IuOAlwI2BmwA3Dm4M66VpunH8RIqGF8d3EzS8eLi4cWquR8ONUXNcao5Lw3GpOQ51x6VuoG4cagbqEOa8VlFZtAF1Auq2SS3wqQd16s1wqTVrNIIGtWaNml+j0qxQ9as0gsYqv1gHr8ly/lvKc0u5MN30yNgYmcAjE8RINV3Svkc6cEk1wyXRNMt5M9aoE6tV8PIzONUGTqWGqdSgUoXmeV5VNeZkgbecI7O4fa11tn+VHNk6L9uP258NWw5VDK6roFzm+HveQ/4z/8id1/dz2011Pj5xhL/d8hs4lQFe8sTtUYcoHWDDmzGstZ8FPrvR79N2nPBGdGLJaMMAkq1lmbVnFoRLRWCjAo3yiqWy4thq69OO1cNtv1Gm4leo2gYV41B2DBVjqDgOZWOoOIZy6/iccSi5HiUvTtmLUXY8Sq5LxUKjAX4jIKg1MfWAZAOSdUuqDsmlpQHJeoNUPUdfvUh/fZq+Ew6ZRyBVDUhUm8Sr/rl/VpkMzkA/3uAQ7uAg7tDp60HcwUG8pf3BQV3pPIugWmX2Qx9i7s8+TOOScd5+8zROus5P5p/K+3k2X7plf9QhikTugvNjcgCufMXGBfR4WXvahceVFyBrJ7ebtcc43nq+v2L7lOPha5pmDc+v4/kVkn4NqnOtHFaFRvXktl9d22eKZcKLy7E0xDPhEktDfCBcZzL4sRQ1L0HFi1Px4lRdj6rrUnU8Ko5DxUDVOFSxVLFUsFSDepgbm1VKfoU5v0rZL1NulCk1ShQbxeVte+Ytkyt+5hBruqTqhkGbYriZZDAIl4EgQX8zTtb3yPgu6boh1TAkagGJio9XruPOTmMefRRbKNIsFMJk+1g8L7xQeo6ceEre7O9Xb5mzqBw6xPF3vIPaAw/y5ReO8ZcHCvzRsWOMXfIy3v+jS3nT08ZJxnRxWc5N/8N6kTHhlU4vsWFv4QFZIOvXoV6EWqG1LkK9EK5XO7a0X8tDY+mxClTzNBolSo5DMWEopRwKjkPROBQcE257cXLxJEc9l4LnkHdc8g7ksRSbDfxqg3QNMlVI1yyZGqSrkK5BX80w3AgYrJcYqJbJHpsk/UCTZKlBrFQ7+48ykVhOWt7ICN7oCO7IKN7oaGt7pLU9GhZ+Xd7qZ62l+C//wvRvv4/6o48y+4In81+e9EMub9b51eA6XjX9an77lfvZNhDtBQ4RWUfGgBcPl3ZhbesCZeXUAq9ROXncr7X2VxaCrfMbZaiXwmVpu1GG8tzytlcv4dVLZB6r4DpdLA2JfkhkwyW5tD0Eid3QHx6z8T4qsRRlL07R8yg5LmVjKDmGkm1S8sNib+VSbBQ5Xs9zuJanUJ8nX89TqBcesyB0jEM2lmXYybK12cdYkGasmWHIjzPYiJGtu/TVDOmaxS37OMUaJl+kMTlJ9dAhmgsL2PrZW2SdgQHcwQG8oWHc0VY+HGnlxxXb3sgITiZzIb/hjtRcXGT2T25j/i//EjPQz5/97Ha+un2BD81XuD4zzm8lfwmY5c3PuDjqUKVDqJCTjeXFwRuG9NqHmI/5dQZreQYri1BtLZVFqObO3K7mTu5XFqCaw8dSdBwWHYdczCGXcMgNOSzGMywkUszHXB7yXBYchwUTsGCbLDYDCFwyVegvQ18F+iuW/orDlkaSsXqSobphoJyj78Q8qR/fTSxXxqmvcnXTcXCHh1uJ7LSib2z0ZMG3VPR1UEufbTQofOWrzH/0o1S++13cXeN8+b8+nduS/85zSxV+ve9pPPfHr+WZT9jCa6/fde4XFBFZC2PCHjGxJKQ28H2sDYvAejm8ELncM6V0aiG4dPGylm8tBai21oWpk8dqBcBigHRrGT39PR0PUkNnWZ4AA4PL+0FygFIsQd51KRjIN0rLBV6+nidXy52y/WBtgbuqUyzUFqh4lVV/dq5xGUgMMJQYYiixhzGnn+3NLFvrKUbqcQarLv1VQ7rcJJ6vEiwu4i/M03j0CJXvfo/mwkL4czuNSafD3Dgygjc2Gl4IHTk9P47hjY7gJDvrYmDj+HEWPvUpFv76YwTFIrUXP5N3XH2IYrLKbfk4B0pT5N7wSf7yw8d5+TU72Tm4kf9opZuokGsTxZrPx779KI4xGMMpa8cQ9nFfOg6rnmcMmNZ5rjG4jiHmOniuwXMcYu45jjkGz11xzHHaa/4SLw7eKGTOSGvn1vTxqosMlmYZLM9BeTa8slqaC9fL+7NQng/3/SoBkHMcFlyHOddlLhljbmiAuWSK+ViCH3gucyZgNqgz26zg2yZYS6rmMliCgTJsrybZ2ehjazXBSNlloFQlM/0IifvvxV0swGpFXyzWKvZOFnru0vbo2CmJzUmn1/yjfTxss0nlBz+k+NWvkrv9dvzpabzt2zjxyy/nPcP/Qr757/z8Yo63XPwabrnvxQyk4YOveVJHFagiIo/JGIilwiUzsvbXC4KTvViWi7v8yaJv6aJlZeHkkj8GU/eE2/XCKS/n0OodEwYbFniZsdYyumL7MhgdO+WxipdgsZZjobbAQnXh5Lq1vVhdZL46z73lh/lmZY5CY8V7J8LFGXYYSY4wmhplNDXOaOoaxuLDbGukGat4jFQ8BkqQylUJ5ubx5+bwZ2epP/II/sG7wqJvFU42u0p+HDvzwujwcGQ9YfyFBUpf/zr5L3yR4te+Fn43eM6z+OzzB/nj4j+zt38Pf7JY46ITB+GNf8uf3Jek2gj4xWfvjSRe6Uwq5NpEvtLgtz93X9RhnMFzDAnPIRFzw7XnkPBcErEV255Dcunx2MljS89LxlwycZdU3CUT90gnWuu4SzrhLT8Wd52N+5Lveq2kdZ5FoLVQL+GUphkqTDFUnGJvcQqKU+HV0+KJ1vYklGYAu1z0zbguM57HTGaImYEBphMpjnseP8BnKigy2ygQLM+fYUnXXLZUE1zsD7Kr3sfWWpLRclgIZnJ14pOP4vzohwTzq1/FdNJpvLEx3LEViWxFMltOcCPDj/t+haBUoj4xSWNyguq991L90d1Uvv99mouL4HmkbnwKj/zCi/izvru4Z/Gfubbi89/nclx002/zhrsu41guxydufQojfRvXnVdEpOM5TtjdMtlPOGf8BWo2TvZEOX0pz7cuWM6E6+l7wu3K6sVSyomRyoyyPTMKmS2Q3Q7ZbeEyuPfkft8WcGNU/SqzldlTlpnKzMnt8gyH5w8zV52jaU8dqMVNuoxdNsa2J21ja2Yr29JXhuvEWJgTKx6ZfAM7P48/M4s/u7TMUD10D/7MDEG5vOrP0x0exhtbkRtXy49jo+GgZ4/jO4gNAvyZGRoTE9QfeYTK3XdT/eGPqN5zD1iLNzbG4Jt/lrueMsT/PP4xZoozvOHSV/Ff7/0mycnvwitvY2rL0/nIR7/Ky6/ZwWVbNeWAnD9jV/liGJUDBw7YgwcPRh1GJILAUmk0w+nvrMUGrXVrP7C2Nafeyf3W1HennGetJbDQDCzNwNJoBvhL6+ZjHAvC7aXHls6r++FS8wNqfpNqI1zX/IDayu3W4+Gxk+deCM8xpOMumYR3suiLu2STMfqTHv2pcJ1NxuhPtdanbIfruOdsxK/o7Jp+2IJXOAHFaSgcC6+Q5iYhP3Fyu1ECwmG451yXKS/GVHaEqfQgJxJpTnguJ2hyolllppE/Wey1DLh97DVj7GkMMF4PW/iGSw79xSbpfA1vsUgwO48/O0tQKKwSaFj0meXpG1pDVKdS4WS/S6zF1uoExSJBsUhzcZFmLnfycWOI791L4uormbp6J1/fmecfTnyBxdoilwQuPzs3zUsHLsd/6R/yy18occd9U/zB65/MS564Y71/8tLhjDF3WWsPRB1Hp+jlHCkbqNk4WeCVZsKeKsvbM2FPleKJ1kXMqXCaplOYsBUvu21FsbcdsluhfxwGd8HAeHgvINAMmizWFpmtzDJVnuJE6QQnSieYKk8xVZriRPkEU6Upqs1TB6xxjctYeoxt6aViL1xvTW9la2YrW+hnoATML7SKvRn82Vmas7OnFX+z4cisp4vFTpm2YXn6huRpFyADS1AqERQKNEtFmjOz2BWDxTiZDMmrriJ93XU0nnYNX0o8yMcOf5zjpeM8cfSJ/Nqlr+VJX/wfsPAovPrPYf9Leec//JC/u2uCr/zqc9g1HE0vG2kv55sf1SLXJhzHkEl016/DWku9GVCtB5TqPuV6k3Ldp1QL16vvNynVfMqNJuVa+NjkYoV7Kw0K1QaFmr9ao9QpkjGH/mSMgVSMoXScoUyM4UycwXSc4XScoUycoXSMoczJ/f6k9/hbA13v5FXKs/8wwi4x+WN4uUm25ifYurLYm29t+xUgLPZmXJcT8SQnsqMcTw9w3Elywq1xrzfJV2IlcokyDJz6NoOJQbamd7PVG2J3o59t1QSjFY/BEvTlfZK1gFg9wKv6UKkSlMthkXbaD9XE47ijI8T37MFk+/C3DFEcTTORqXO4v8Sh2kN8d+qrVPwK3qMOz/YNb5iZ4vr4KOam91O54rX84se+x7/8eIbffPmVKuJERNqVGzt3DlsSNMPCrnA8vHi52vrY95Z7qpwiOQADu3AHdjEyMM7IwDj7BnfBwB7Y/gzo2xqO9k34/SFXyy0Xeqev75u/j68d/Rq15qmDkTnGYTQ5GhZ4O7ay9bKtbE3vZWv6xuWibyw5ilushAXe0jIzS3NhnqBcISiXT11mi2EX2iXG4GQyxPdchJPpwxsdITY+TmznTrzxcR7uK/O16YN8Y/IbfOeePyOwAdduuZZ3PeXXeeb0ozh//0thV9z/43bY83TunszxN985ypuetkdFnFwwtchJRwkCS7HuU6j65CuNk+tag3xlaTtcL5YbLJTrLJYbzJfrLJTq+GeZqNV1TFjcpeOM9iUYy65YVuxvySYYSsc35t5Ba8MuLrmJcMlPQu5oWOAt7eePLV8NrRjDlOtyItXPVN8wU6ksU/EkU45hBp+ZZoW5eoGA1VtG406cTCxDOpYm7sbxHA/PeFjs8vxJZb/MYm2RwJ58Dc847HH7uL7W4GmzE1xfKZPZeT1c+9PwxNdyJNfk1r86yOGpAu975dW89vrd6/+zkq6gFrkLoxwpHaPZCFvv8sdg8cjJvJabaOW1o2EX0JUcD/p3wsCuViveyvXusFVvxWjbK4u95aV06vpE6QRl/8wul8PJYbamtzKcGmYwMchQYoiBxACDiUEGE4P0xftIuAnibnx5DdBohpPU15o1FmvhPYKzlVmO5I/wUO4hHso9RKV1QXZP/x5u3nMzt+y9hb2VEnzh1+Hhr8NFT4dX/Rn078BvBrzij/6VE/kqX/6VZzOQuuBZi6VLqUVOupLjmLA7ZTJ2waM6WWsp1HwWSycLu/lSnYVyuMyXGiyU6swWa/xgYpHpfI1K48wJV13HMNoXP6XI29afZNtAiu0DSbYPJtnen6I/dYGtfMaEo3umh2H7E1c/J2iGVz3zk6RyR9mTm2TPyqLvxP1hN8+WJrDgOMwlM8ynBsglM+TiaRZjCUpunLJjKAYBDRvQpIFPDWMtMQtxDCmbYtjxGKmVGS3McElpgd0Nn5hxYfwAHHg5XP1TMLYPay3//KPjvOvTdwPw5z9zPc/Zt+VCfkUiItIN3FhYeA2Mw64bVj+nmg9z12KrsFtZ6D38jfA2BXvahci+bWFxN7gbM7CLwcFdDA7sZt/gbthyXTjf32mK9eKpxV2r6+ZUeYqF6gKP5B4hV8tRbBQf98fdmt7K3oG9vOqyV3HFyBVcv+16tqW3wuRd8MX/Gw59GhID8JIPwrU/E94LCfzFvz7CjyZzfOgNT1YRJ4+LCjnpGcacLAJ3j5xf94VSzWemUGOmWAvXhRrThery9kyxxqFjeWaLNU5v7EvHXbYNJNkxkGLbQDIs8lrF3tLxCy72HBcGdobL2ZJjo9K6L+8obm6S0cJxRpdudK/Mh+t8a3TOeimcdHc1biJMitltYZeXS26Asf2w5XLYcW3rZvzQo3MlfvOf7uHL905z9c4B/vAN1573z1hE5LFM5as8NFMKR252wpGbzcqRnM8Y6Xlpe+m88FzPdYitGJ15abRm1zEaTTcKS4O6bNm/+uPNxooWvaOtgu9IuD/5XbjnMxCcNupzanhFK95FMLiLvoFd9A3u4pLh/bDjxlO7Sa7QaDbI1XMsVhcp+SXqzTq1Zi1c/BrGGGJOLFzcGIOJQYaTwwwnh5db7LA2HETmux+Hu/8epu6GeB8867/BjW8NRw1tuedYnvd/4TA3Xb6Fn7h6+3r8RKUHqZATeQyZhEcm4bFn9LEnKvWbAdOFGsdzVY7nKpzIVTm2WOVEvsKxxSrfvH+W6UJ11WJv+0CSHYOp5UJv52AqbNUbSLFjMEk6foH/TWMpGLkkXM6HtdCsh5PjYsFLghs/a7Jb6YHpAn/69Yf5u+9OEHMN77plPz/79D147iYPOCMiXesr903zzn/40Ya9vjEQc5am5THEPQevtR9vTdcTc8PRmZMxh6Tnkoy74Tp22vHW9tKIzUnPCQfvSnj0tfJJX9wjk3D1d/Jc3BgMXRQuqwmCsPvmcqG3ouCbvR8e/Eo4h99K8eyp3TWXBmJJjxBLDTGaHGQ0tSUcmOWxcmDQDAc3W5iA3Ldh5j44/oOwwCyeCM/ZcW3YAnfVT51y4RPCi8Rv+8R3GUzFeP9PPVEXEuRxUyEnsg4812HHYIodgylgaNVzTi/2ji9WOZ6rcmyxwvFchftOFJgpnNk6NpCKsX0geUaBt30gtdzat6aROo0J7zvwzm9qgGOLFe64d4p//OFx7nx4nrjn8NNPvYhffs4lbOnvrElaRaT9Pe/yLXzi55+6PCrzGSM1B5wxcvPSyM5L66WRnOvNAL81OnO9NXKz3wyot9b+ynOaJ8+pNwOqjSaVepOFUoNqa5TmaqMZLn5A8yz3YJ9NwnPoS5yckme50EuEhV4m4S0P3LW8pMP1YCpGfypGMhbNHGltwXGgf3u48JQzH7c27Hmy1Iq31IVzsVX0Hf32mffpLTFu2JLmuOG9e24MMGFh2CiHk8Cf+gQYvQwufibseSZc9sJWXGdqBpZf+dT3eXi2xMfeoml5ZG1UyIlskvMp9up+wFR+qbirMtkq8o4vVjmWq3LXkQUWy6d2JTEGRvsS7FjqujnYKvpa2yOZeGuahhjuBQzS0mgGrZbFCoenCvxwIscPji5y/3R4H8He0QzvePHl/NR144wqEYnIBtnan2RrB1wkaiwVe42VRV5ApdGkVPcp1cKlWGuu2D712EK5ztGFcuvxJsXaKsPkr5DwnFMLvVaxN5KJM5xJMNIXZ7SvtZ2JM9IXv/BeHp3KmHCS9swI7Hjy6udU8+F9eZWF1mTrCycnXK+XIPDD7puBH7YAxlIQT0MsE85LOzAeDtAyvBcSfecMyVrL//PP9/CFQ1P8xkuv4GmXnOfctiJn0SP/m0U6Q9xz2DWcfswhiMt1P2zVWwyLrGPLhV6FB2aKfOP+GUr1MwdpAcgmwvn4skmPuOfgOmFXIs9x8IOAUi38ElKs+cyddt/fSCbOE8cHeNV14zx//1YuGcuoO4iISEvMdYi5Dtnk+g1a0Qws+UqD3FmWpccWy+H6eK7KfScKzJVqZ53LNRVzGc4sFXhxRvrCgm9rNtkqmhNs7U8ylk10f4tfsh+SV2zKW1lr+d0vHObPv/UIb376xfzs0y/elPeV7qZCTqTDpOMel4z1ccnY6lf/rLXkq/5yl82F0pmJP1/1aQYnJ3/3mxbPcdgxGCMV98jEXbb0J9k5GN6/t3esjx0DSRVuIiKbyHVMOPdpJn7Bzy3XfeaKdeZKdeaKtda6znyptnx8pljj8IkCs8U69eaZhd9gOsbWbJItreJuqcjbkg0H8No5lGIkE1duOIdmELbE/fm3HuH1N+zm//qJswzwInKBVMiJdBljzHIXm/3b+8/9BBER6TrpuEd62DuvSaattSyWG0wVqkzla0zlq0znT25PFWo8MD3LdKF2xr2ACc9h52CKnUPhfds7h8JBu3YMphgfSrG1f433cXe4XKXBf/rE9/iXH8/ws0/fw7tfcoUKX1k3KuREREREepgxJ1v+Lt929vOagWW+VF++l/vYYoXJ5aXKvcenmS2eOmiXMbA1G7be7W7dOnDRcJqLRtLsHk4zlk10bWHzzftnefvf/5DpQpXfesXVvOEpu6MOSbqMCjkREREROSfXMYxlE4xlE1y1c2DVc6qNZjhY10JY6E20Cr6JhTJ3PjzP//7+JHZFo14q5p4s8EbCZanYGx9Kd2Rr3uRihQ988TD/8N1J9o5m+JtfuJFrd68+yJnIWqiQExEREZF1kYy5XDya4eKzzL9a85tMLlR4dL7MkbkyR+bLPDpX5sh8iW8+MHPKIC2Oge0DYUveUoG3azhsyds1lGK4ze7Pe2C6wEf/9VH+5uBRAH7x2ZfwX55/WfcPGiORUSEnIiIiIpsi4bnsHetj7yoDdllrmSnUlou8cF3iyHyZL987xWyxfsr5mbh7RnG3u9Vlc3wovSkF1MRCmS8emuLzd5/gzkfCuVVfcc1O/tPzL2PnYGrD3196mwo5EREREYmcMYYt/Um29Ce5fs/wGY+X6z5H5yscnQ9b8o7Ml5lYKPPoXIlv3j9LpXHq1Dtbsgl2DafZ2p9gSzacUmFLq2volmyS4UycvqRHOubiPMY8q0vTQEwVqjwyW+Lh2TKHjuX47qMLHMuFk4Pv25rlv928j9ddv0uTfMumUSEnIiIiIm0vHffYty3Lvm3ZMx6z1jJbDCdUP9pq0Qu3Kxw+UeAb989SqJ59gvVM3CWT8Ii5DrZ1E58FilWfwioTs28fSHLdRUO8ZfcQz9k3tmoLo8hGUyEnIiIiIh3NmJMDsZxtYJFKvclMocZ0ocp0ocZCuU6p5lOsNSnVfEo1H781vcJS+1wm4TGYDqf0Ge1LsGckw0WjafrXceJ3kcdLhZyIiIiIdL1U3A3voRs599x6Ip2g88Z0FRERERER6XEq5ERERERERDqMCjkREREREZEOo0JORERERESkw6ypkDPGvNoYc8gYExhjDpz22DuNMQ8YYw4bY25eW5giIiKdRTlSREQ20lpHrbwbeCXwJysPGmOuAF4HXAnsAL5sjHmCtbZ55kuIiIh0JeVIERHZMGtqkbPW3mutPbzKQy8HPmmtrVlrHwYeAG5Yy3uJiIh0EuVIERHZSBt1j9xO4OiK/YnWMRERkV6nHCkiImt2zq6VxpgvA9tWeehd1trbz/a0VY7Zs7z+rcCtALt37z5XOCIiIm1DOVJERKJyzkLOWvv8x/G6E8CuFfvjwLGzvP5twG0ABw4cWDWRiYiItCPlSBERicpGda38DPA6Y0zCGHMxcBlw5wa9l4iISCdRjhQRkTVb6/QDrzDGTAA3Av9sjPkCgLX2EPAp4B7g88BbNRqXiIj0EuVIERHZSGuafsBa+2ng02d57L3Ae9fy+iIiIp1KOVJERDbSRnWtFBERERERkQ2iQk5ERERERKTDqJATERERERHpMCrkREREREREOowKORERERERkQ6jQk5ERERERKTDqJATERERERHpMCrkREREREREOowKORERERERkQ6jQk5ERERERKTDqJATERERERHpMCrkREREREREOowKORERERERkQ6jQk5ERERERKTDqJATERERERHpMCrkREREREREOowKORERERERkQ6jQk5ERERERKTDqJATERERERHpMCrkREREREREOowKORERERERkQ6jQk5ERERERKTDqJATERERERHpMCrkREREREREOowKORERERERkQ6jQk5ERERERKTDqJATERERERHpMGsq5Iwxv2uMuc8Y80NjzKeNMYMrHnunMeYBY8xhY8zNa45URESkgyhHiojIRlpri9yXgKustU8Efgy8E8AYcwXwOuBK4EXAHxlj3DW+l4iISCdRjhQRkQ2zpkLOWvtFa63f2v02MN7afjnwSWttzVr7MPAAcMNa3ktERKSTKEeKiMhGWs975N4MfK61vRM4uuKxidYxERGRXqQcKSIi68o71wnGmC8D21Z56F3W2ttb57wL8IGPLT1tlfPtWV7/VuDW1m7NGHP3uWLqYqPAbNRBRKSXPzv09ufXZ+9d+6IOYK2UIzdNr/9f6eXP38ufHXr78/fyZz+v/HjOQs5a+/zHetwY8ybgJcBN1tqlRDQB7Fpx2jhw7CyvfxtwW+u1DlprD5xH3F2plz9/L3926O3Pr8/em58dws8fdQxrpRy5OXr5s0Nvf/5e/uzQ25+/1z/7+Zy31lErXwS8HXiZtba84qHPAK8zxiSMMRcDlwF3ruW9REREOolypIiIbKRztsidw4eABPAlYwzAt621v2itPWSM+RRwD2F3krdaa5trfC8REZFOohwpIiIbZk2FnLX20sd47L3Aey/wJW9bSzxdoJc/fy9/dujtz6/P3ru6+vMrR66rXv7s0Nufv5c/O/T259dnPwdzssu+iIiIiIiIdIL1nH5ARERERERENkFbFnLGmP9ojDlsjDlkjHl/1PFsFmPMe4wxk8aY77eWW6KOKQrGmP/TGGONMaNRx7JZjDG/aYz5Yev3/kVjzI6oY9pMxpjfNcbc1/oZfNoYMxh1TJvFGPPq1t+6wBjTE6NzGWNe1Pob/4Ax5h1Rx9NplCN7N0f2Yn6E3s6RvZwfQTnyXDmy7Qo5Y8xzgZcDT7TWXgn8XsQhbbYPWmuvaS2fjTqYzWaM2QW8ADgSdSyb7HettU+01l4D/BPw7ojj2WxfAq6y1j4R+DHwzojj2Ux3A68Evh51IJvBGOMCfwi8GLgCeL0x5opoo+ocypG9myN7OD9Cb+fIXs6PoBz5mDmy7Qo54JeA91lrawDW2umI45HN9UHg1zjL5LjdylqbX7Gbofc+/xettX5r99uE82r1BGvtvdbaw1HHsYluAB6w1j5kra0DnyQsTOT8KEf2rp7Mj9DbObKX8yMoR3KOHNmOhdwTgGcaY/7dGPMvxpjrow5ok72t1Xz+EWPMUNTBbCZjzMuASWvtD6KOJQrGmPcaY44Cb6S3rjae7s3A56IOQjbMTuDoiv2J1jE5P8qRPZgjez0/gnJki/Jj97ugHLnWeeQeF2PMl4Ftqzz0LsKYhoCnAtcDnzLG7LVdMrzmOT77HwO/SXil6TeBDxD+p+0a5/j8vw68cHMj2jyP9dmttbdba98FvMsY807gbcBvbGqAG+xcn791zrsI59X62GbGttHO57P3ELPKsa74+75elCN7M0f2cn6E3s6RvZwfQTnyNBeUIyMp5Ky1zz/bY8aYXwL+oZWU7jTGBMAoMLNZ8W2kx/rsKxlj/pSwH3hXOdvnN8ZcDVwM/KA1ce448F1jzA3W2hObGOKGOd/fPfBx4J/poiQF5/78xpg3AS8BbuqWL6VLLuB33wsmgF0r9seBYxHF0paUI8+tG3NkL+dH6O0c2cv5EZQjT3NBObIdu1b+b+B5AMaYJwBxYDbKgDaLMWb7it1XEN7g2ROstT+y1m6x1u6x1u4h/Id8bTclqcdijLlsxe7LgPuiiiUKxpgXAW8HXmatLUcdj2yo7wCXGWMuNsbEgdcBn4k4pk7yv1GOhB7Kkb2eH6G3c6TyY8+5oBwZSYvcOXwE+Igx5m6gDrypG68+nMX7jTHXEDahPgL8QqTRyGZ6nzFmHxAAjwK/GHE8m+1DQAL4UuuK87ettT3xMzDGvAL4A2AM+GdjzPettTdHHNaGsdb6xpi3AV8AXOAj1tpDEYfVSZQjlSN7US/nyJ7Nj6AcyTlypOmdv/8iIiIiIiLdoR27VoqIiIiIiMhjUCEnIiIiIiLSYVTIiYiIiIiIdBgVciIiIiIiIh1GhZyIiIiIiEiHUSEnIiIiIiLSYVTIiYiIiIiIdBgVciIiIiIiIh3m/wfYuVh+4scwCgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "f, ax = plt.subplots(1,2,figsize=(15,5))\n", "\n", "pore_sizes = [5,10,20]\n", "slit_pore = Pore1D(Geometry.Cartesian, 100*ANGSTROM, potential, 4096).initialize(state)\n", "cylindrical_pores = [Pore1D(Geometry.Cylindrical, r*ANGSTROM, potential).initialize(state) for r in pore_sizes]\n", "spherical_pores = [Pore1D(Geometry.Spherical, r*ANGSTROM, potential).initialize(state) for r in pore_sizes]\n", "\n", "for cpore, spore, size in zip(cylindrical_pores, spherical_pores, pore_sizes):\n", " ax[0].plot(cpore.r/ANGSTROM - size, cpore.external_potential.T, label=f'${size}$')\n", " ax[1].plot(spore.r/ANGSTROM - size, spore.external_potential.T, label=f'${size}$')\n", "for a in ax:\n", " a.plot(slit_pore.z/ANGSTROM-50, slit_pore.external_potential.T, label='$\\infty$')\n", " a.axis([-6,0,-20,55])\n", " a.legend()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "pmin = 0.2*BAR\n", "pmax = 2.5*BAR\n", "cyl10 = Adsorption1D.equilibrium_isotherm(func, T, SIArray1.linspace(pmin, pmax, 151), Pore1D(Geometry.Cylindrical, 10*ANGSTROM, potential))\n", "cyl20 = Adsorption1D.equilibrium_isotherm(func, T, SIArray1.linspace(pmin, pmax, 151), Pore1D(Geometry.Cylindrical, 20*ANGSTROM, potential))\n", "sph10 = Adsorption1D.equilibrium_isotherm(func, T, SIArray1.linspace(pmin, pmax, 151), Pore1D(Geometry.Spherical, 10*ANGSTROM, potential))\n", "sph20 = Adsorption1D.equilibrium_isotherm(func, T, SIArray1.linspace(pmin, pmax, 151), Pore1D(Geometry.Spherical, 20*ANGSTROM, potential))\n", "car10 = Adsorption1D.equilibrium_isotherm(func, T, SIArray1.linspace(pmin, pmax, 151), Pore1D(Geometry.Cartesian, 20*ANGSTROM, potential))\n", "car20 = Adsorption1D.equilibrium_isotherm(func, T, SIArray1.linspace(pmin, pmax, 151), Pore1D(Geometry.Cartesian, 40*ANGSTROM, potential))" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAD4CAYAAAD1jb0+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAA1u0lEQVR4nO3dd3xUVdrA8d+ZFGLoJJRAAiF0QkkwBpAqiAooiCIWFFCxK1jWRQVdXHUt676vBctiAeOrLqgrNkAUpCoISOgt9ECAECChps15/zjpM4HM5CYzyTzfz+d+ZjLP3DtPrsPjybnnnqO01gghhKh6bJ5OQAghhHukgAshRBUlBVwIIaooKeBCCFFFSQEXQogqyr8yPyw0NFRHRkZW5kcKIUSVt3bt2mNa64YlX6/UAh4ZGcmaNWsq8yOFEKLKU0rtc/a6dKEIIUQVJQVcCCGqKCngQghRRVVqH7gz2dnZJCcnc/78eU+n4lFBQUGEh4cTEBDg6VSEEFWExwt4cnIytWvXJjIyEqWUp9PxCK01aWlpJCcn07JlS0+nI4SoIjzehXL+/HlCQkJ8tngDKKUICQnx+b9ChBCu8XgBB3y6eOeTcyCEcJXHu1CE8AV2bcemTHtpZuJMzuecJzMnk6zcLOzaTrewbgxqNYhcey6v//a6w/49wnvQL7If57LP8eaqNx3ifVv05fKIy0k/n867q991iA9qNYi4pnGknknlgz8/cIgPbTOUrk26kpyRTML6BIf4iPYj6NCwA7tP7OaLjV84xG/pdAutGrRi27FtfLXlK4f42K5jiagbwfrD6/l+x/cO8Xu63UPjWo1ZfXA185PmO8Qfjn+Y+pfUZ/n+5Szas8gh/pfL/0JwQDALdy9k+f7lDvHJfSfjb/Nn3s55/HHwj2Ixm7LxbL9nAZizbQ6JhxOLxS/xv4RJvScBMGvTLLakbikWrxtUl8d7Pg5AwvoEko4nFYs3qtmIh+MfdsjJClLAXTBu3DiuvfZaRo4cybRp03jjjTfYtWsXqamphIaGAqY/e+LEicydO5fg4GBmzpxJt27dPJy5qExL9y1l89HNbD22lW3HtrH12FYGRQ3i4+EfA3D/D/eTmZtZbJ+HLnvIFHCdy1MLn3I45lO9njIFPOccTy982iH+0oCXuDzick6eP8kzi55xiNcMrElc0ziOnDnC5EWTHeJhtcLo2qQrB9IPOI23DWlLh4YdSDqexJRfpzjE45rG0apBK7akbuHZX591iA9oOYCIuhEkHk50Gh/WbhiNazVm1cFVPLf4OYf47V1up/4l9Vm2bxl/W/w3h/iDlz1oCviehby8/GWH+FO9nzIFPGkeb//xdrFYgC2goIB/u/1bZibOLBavH1S/oIB/ueVLvt76dbF487rNCwr45xs/Z8GuBcXi0Y2iK6yAq8pc0CEuLk6XvBNz69atdOjQodJyKI+iBXzdunXUr1+f/v37s2bNmoICPnfuXN5++23mzp3LqlWrmDhxIqtWrSrT8avSufB1+07uY/2R9QUFuk5gHd4cbFrGHd7pwLZj26gZUJP2oe3p0LADV7a8krExYwE4kH6AAL8AAv0CCfQLxE/54WfzI9AvEK0153Mcr4X42/wJ8AtAa+1Q/PPj/jZ/tNZk5WY5xP1sfvjb/LFrOzn2HMd4Xg6lxf1t/tiUDbu2k2vPdXr8/Lhd250eXylVprizmmRTtoK4MwqFUsrpvgXvuUDc27swlVJrtdZxJV+XFjiQkJDA66+/jlKKqKgoEhMT2bFjBwEBAWRkZNClSxd27txZbJ/Y2Finx/r2228ZM2YMSil69OjByZMnSUlJISwsrDJ+FVEJth3bRsd3OqIxxaBJrSb0j+xfEP/ixi8IuSSE8DrhTgtDRN2IUo+tlOKSgEsuGA/yD7pgvIZ/jVLjNmUj0C+wXHGbX+mXzmzKVtBV5G6cC9TSC+0LFy/E3l6oXXXRAq6UigASgCaAHZiutX5TKdUAmAVEAnuBUVrrE+XK5tFHITGxXIdwEBMDb7xRanjz5s289NJLrFixgtDQUI4fP84TTzzBjz/+yPXXX89//vMfbrzxxjKPzz548CAREYX/QMPDwzl48KAU8GrkXPY52oe2Z1KvSQxvP5x6QfWKxWOaxHgkL+F7yjIKJQd4QmvdAegBPKSU6gg8BSzUWrcBFub9XOUsWrSIkSNHFnSBNGjQgPHjxzNjxgwAZsyYwZ133lnm4zn7E626/V/f18WGxbLloS2MjRnrULyFqEwXbYFrrVOAlLznp5RSW4FmwHCgf97bPgEWA5PKlc0FWsoVRWvtUGB79erF3r17WbJkCbm5uXTq1KnMxwsPD+fAgQMFPycnJ9O0aVPL8hVCiHwujQNXSkUCscAqoHFecc8v8o1K2edepdQapdSa1NTUcqZrvYEDBzJ79mzS0tIAOH78OABjxozh1ltvdan1DTBs2DASEhLQWrNy5Urq1q0r3SfVzMYjG+n9cW+H4WhCVLYyF3ClVC3ga+BRrXVGWffTWk/XWsdpreMaNnSYj9zjoqOjmTx5Mv369aNr1648/rgZDjR69GhOnDjBrbfe6nS/t956i/DwcJKTk+nSpQvjx48HYMiQIURFRdG6dWvuuece3n3XcUyuqNoyMjNYcWAF6efTPZ2K8HFlGoWilArAFO/PtNb/zXv5iFIqTGudopQKA45WVJIVbezYsYwdO7bYa8uXL2fkyJHUq1ev4LWZM2cWPJ8wYQITJkxwOJZSinfeeaeiUhVCiAJlGYWigI+ArVrr/ykS+g4YC7yS9/hthWToAY888gjz5s1j7ty5nk5FCCFKVZYWeC/gDmCjUiox77VnMIV7tlLqbmA/cFOFZOgBb7/99sXfJIQQHlaWUSjLKX1o/UBr0xHC+9UKrEV8s3jq1Kjj6VSEj5M7MYVwUdcmXVk1vmzTIwhRkbxiOlkhhBCukwIuhIvWpayj6/td+f3A755ORfg4KeAuGDduHF99ZeY6Hj16NO3ataNTp07cddddZGdnA+bOzgkTJtC6dWu6dOnCn3/+6cmURQU4m32WDUc2cDrrtKdTET5OCribRo8ezbZt29i4cSPnzp3jww8/BGDevHns3LmTnTt3Mn36dB544AEPZyqEqK7kIibuTSc7ZMiQgufx8fEkJycDMp2sKIXWcPw47N4Np09DVpbZ2rWDtm3h7Flwdt9B167Qpg2kp8OCBY7xuDho2RLS0mDhQsd4z54QEQFHjsDixY7xvn0hLAySk2HFCsf4gAHQsCHs2QPO5rW/+mqoXx927oQSc/0DcN11UKsWbNkC69c7xkeMgKAgE9u82TE+ahT4+5tjb99ePKYU3Habef777+bcFhUYCDfljW5euhT27y8eDw6GG24wzxcuhEOHisfr1oVhw8zz+fOh5FQgISGQXwe+/x5OlJiMtUkTuOoq83zHDvPf2Wpa60rbLr30Ul3Sli1bir/Qr5/j9s47JnbmjPP4jBkmnprqGLuITZs26bZt2+rU1FSttdZpaWl63Lhx+ptvvtFaa/3vf/9bP/7441prrceOHau//PLLYvtnZWXp2NhYvXTpUq211kOHDtXLli0riA8YMECvXr36onlo7eRcCK+0fN9yzVT0gqQFjsH0dK3XrNH688+1fv55rfO+R/rIEa1NGS++vfiiie/b5zz+xhsmvmmT8/hHH5n4ypXO47Nnm/jPPzuPz51r4v/9r/N4/nf5k0+cxxMTTXzaNOfxXbtM/JVXnMePHjXxyZOdx8+eNfEJExxjfn6F5/3OOx3j9eoVxkeOdIxHRBTGr7nGMd6xY2G8Vy/HeHx8YbxLF8f4wIGF8X//2+l3qayANdpJTfX5Fnhp08m+9tprXH/99cyYMYMPPnBcQzDfgw8+SN++fenTpw8g08n6grp+wQxs1IMGK9dDkt20QrWGqCjYu7fwjUrBI4/A9debVuwbb0BkJNSrZ1qHgYHQrJl5b5MmsHGj44c1aWIeW7WCTZsc4/kzXXbu7LwFm3/8Hj1MK7ik/Lnrr7zSebxFC/M4fDhs3eoYj4w0j7fdZo5RUni4eRw/3pyHkurXN4+PPgpjxjjGa+QtTjF5Mjx8gWXJXn4Zni6x1JytSA/xW2+Z9xTlX6T8ffQRnDtXPF50DYBZs+B8iZWSahRZOOPHH81fVEVdUmRhjltuKT33cvC+Au7sz7x8wcEXjoeGXjjuhC7HdLLPP/88qamp/Pvf/y54TaaTrUYyM+Hw4cIi9vDDsHAhnXbu5JfcXGAlXHGFKeBKwejRULu26fJo29YU3fx/xErBxImlf1ZgIFxo2uKgIIiOLj0eHAwdO5Yer1ULLrRcX+3aF47XrWu20tSvX1iMnQkJMVtpQkPNVppGjcxWmsaNzVaai3VhXuzfaP7/CEuT/z+q0tSpmJu+vK+AV7KBAwcyYsQIHnvsMUJCQjh+/DgNGjQomE722WcdF2AF+PDDD/npp59YuHAhtiL/px82bBjTpk3jlltuYdWqVTKdbFXyyy+mAbB5s2mNJiWZVnX+9Y/sbGjfHm680RS79u1Nsc734oseSVv4Lp8v4EWnk/Xz8yM2NpaZM2cyevRopkyZUup0svfffz8tWrSgZ8+eANxwww0899xzDBkyhLlz59K6dWuCg4MLVvYRXiItzVwwS0w0244d8Ntv5s/tWbNgxgxTlDt1gptvLt4qzvtLa82hNYz6chQJDRPoXUEtq+rCWZdiqe/FhfdW0HEr8tj5C0NbyecLOLg3nWxOjuPK3SDTyXqVs2dh3jzTxVGrFrz2GkwqsmhU06ZmzdT0dPPn/6uvwrRpxfs2nTifc549J/c4XT0+n9aaI2eOsP3YdranbS943Hl8JxmZGWTlZpGZk0lWblapK607HNNLCpFwz7zR87im9TWWHlMKuBMynWw18eqr8Pe/w5IlZrjcFVeYIh4TY4bnlexTbdDA5Y+wazv7Tu5jc+pmNh/dzJZjW9h2bBvbj20nPbNwwYcg/yDahrSlS+Mu1A+qT6BfIDX8ahDoF4ifza/Mn6cutGS7s/e7cAHdlWNX5+NW1LHbNGhz8Te5SAq4EzKdbDVx5oy5iNi9u/n5ssvM5gatNQcyDrD56GZ+2PEDAA/NfYiDGQc5k32m4H3NajejfWh7RnceTbvQdrQLaUe70HY0r9vc8j+fhZACLqo3m+2iXSIlZWRmsOHIBtalrGP9kfVsOrqJLalbOJV1qtj7agbUZHy38UQ3jCa6UTQdG3aUVepFpZICLnyW1prDpw+z7vA6Eg8nFjwmHU8qeE9ocCidGnViTNcxdGrUieiG0fjZ/Pjnb/9kar+pdG3S1YO/gfB1UsCFz9Fac+e3dzIvaR5HzxQu5RpVP4rYJrGM7TqW2CaxxDSJoWntpk77OL+5+ZvKTFkIp6SAi+rrlVfgH/9weDnbns0n6z+hR3gPJveZTEyTGLo27krdoAvcqCKEF5KrKi4qOqVsSddccw316tXj2muvLfb6nj176N69O23atOHmm28mq+Qtt6Ji+PubOxxLcV3b65jQfQJ9W/R1qXivTF5Jo382Yum+pVZkKYTbpIBb6Mknn+TTTz91eH3SpEk89thj7Ny5k/r16/PRRx95IDsf9M038Pjjlh82Ozeb1LOpZOdmW35sIVwhBRw4c+YMQ4cOpWvXrnTq1IlZs2YRGRnJpEmTiI+PJz4+nqSkwgtbS5cu5fLLLycqKqpYa3zgwIHUrl272LG11gUTZoG5aWjOnDmV8nv5vBUrYPp0T2chRIXxqj7wRx81dzdbKSbGTAJ3IfPnz6dp06b8+OOPAKSnpzNp0iTq1KnDH3/8QUJCAo8++ig//GDG/6akpLB8+XK2bdvGsGHDCoqzM2lpadSrVw//vJnPwsPDOXjwoBW/mhDCx0kLHOjcuTO//PILkyZNYtmyZdTNm3Utfx6UW2+9ld9/L1z/8Prrr8dms9GxY0eOHDlywWPL9LJCiIriVS3wi7WUK0rbtm1Zu3Ytc+fO5emnn+aqvFU0ihbaos9rFLkx5GLzTYSGhnLy5ElycnLw9/eX6WWrgUY1G3FHlztoUquJp1MRPk5a4MChQ4cIDg7m9ttv5y9/+UvBQsSzZs0qeMyfddBVSimuuOKKgr7yTz75hOHDh1uTuLgwP78LjkJxV7vQdiSMSCC60QXm5xaiEnhVC9xTNm7cyJNPPonNZiMgIID33nuPkSNHkpmZSffu3bHb7XzxxRcXPU6fPn3Ytm0bp0+fJjw8nI8++oirr76aV199lVtuuYUpU6YQGxvL3XffXQm/leDVV83mJXJy4MCBwiUxW7c2aySkpjpfDCcmxsQPH4Zt2xzjcXFmksXk5MIpy4vq0cNMBbNvn+NykQC9epn/v+3aZd5TUv/+ZiaC7dvNZxSllFkuE0zuKSnF4/7+0K+feb5hAxw9WjweFAS9e5vnf/5pZvktqmZNuPxy8/yPP8yEkUXVrQvx8eb5b7+Zc1pUSAhceql5vmyZ42I7jRqZ8wuwaJGZ6r2opk3NIkdgliK1l5gwsnlzs36G3e58qdKoKLOmR1aWWW4zKsosf2o5Z+usVdRWpjUxvUSLFi0K1smsLN56LqqbzJxMzVT0S0tfcmv/5fuW66AXg/Si3YscYqdPa71+vVlics8e89qKFVpHRZklHIsumTh/vol/9ZXz5SBXrDDxmTOdxzdsMPG333Ye373bxP/xD+fx/K/30087j58/b+IPP+wYCwgo/J3HjXOMN2hQGL/xRsd48+aF8auucoxHRxfGe/Z0jPfoURjv1MkxPmhQYTwy0jF+ww2F8ZAQx/iYMYXxwEDH+EMPmVhmpvNz99RTJn7smPl50iSHr4pLkDUxhc/54guzwk6RJe+skJVt53zOeY6m2qGlab3eeadZG6LoAKP334f77jPLYV52mVkWMSrKtB4DAyE21ryvb1/TCiwpfwW1q6+GX391jLdsaR5HjChsLRaVv5zm6NGFrdmi8ldIu+ceuMbJNNX5S0JOnGgWhy+q6HX4p5+Gu+4qHi+63OQLL5gRZkUVnV/sf/4HTp4sHg8OLnz+/vtwqvg8YhQdrZuQ4NjCLrr621dfmdXxiio6c/C8eeavo6IaNix8vnSpKctF5a/e5u8PRcY3FMi/zFW3LqxcefEV3dyltAuTvpdXXFycXrNmTbHXtm7dSocLrcXnQ+RcWOyJJ8w48BL/+rNys6jxYg1eGvASz/R5ptTds7JMoc3MNGvqbt9uivSu7GXkju3LbTm/8NkLA0lPh8GDzZ/MbdqYrVUr8ydzrVoV/UsKX6CUWqu1jiv5urTAhcC0ojZsMAu/b95stsGDYeZMU8QTEkyLtnNniG8H/0fhIux165p+WCEqmxRw4TNOn4aNG2H9JvPzTz/BM33M87vvhq1bzZ/m0dFw3XWFBVopOHKksNtg2T74v5nQPKLSfwUhipECLqqlo0dhx8GW5A10YOxY+PTTvL5MP+BZM7oj32efmZELERHF+3fzFX0trHYYD8Q9QLM6zSrwNxDi4qSAi2ph9Wr47jtYt85shw6Bv+0BTjd7kxqYIW+tWpmhY207QIfPTVHPl39BsSxaN2jNu0PftfpXEMJlciOPi0qbTjYxMZGePXsSHR1Nly5dCm4CAplO1krp6WZc7T/+AcOHF44/XrIEXn4Z9u41xfpf/4KffvbDtssMkB47Fp57DoYNMyNBykNrTY49p8yryQtRUaSAWyQ4OJiEhAQ2b97M/PnzefTRRzmZNzZKppMtvxUrzI0T9eubvunJk83NLYcOmfh995nBJps2ma6Sxx83hTx/KJyVlu1fRsALASzeu9j6gwvhAingWDOdbNu2bWnTpg0ATZs2pVGjRqSmpsp0shZp1Mi0nJ9/HubPN3fubd9eeLdd7drmrsNiZs40g6CFqKa8rg+8f3/H10aNggcfhLNnYcgQx/i4cWY7dgxKzuy6ePHFP9Pq6WT/+OMPsrKyaNWqlUwna4FVq8yNF3mnv+w2bjQd40JUU9ICx9rpZFNSUrjjjjuYMWMGNptNppO1wA03wGuveToLIbyP17XAL9RiDg6+cDw0tGwt7pKsmk42IyODoUOH8uKLL9KjR4+8nGQ6WSFExbhoC1wp9bFS6qhSalOR16YqpQ4qpRLzNicdG1WHFdPJZmVlMWLECMaMGcNNN91U8LpMJ1v9hNcJ56+X/5XmdZt7OhXh48rSAp8JTAMSSrz+v1rr1y3PyAOsmE529uzZLF26lLS0NGbOnAnAzJkziYmJkelkPaVBA2jRwvLDRtWP4tVB3jNNrfBdZZrMSikVCfygte6U9/NU4LSrBbwqTWYVGRnJmjVrCA0NrbTP9NZz4WnNmpmL1x98YM3xyjqZVWmyc7M5nXWaWoG1CPCrgHGKQpRQ2mRW5bmI+bBSakNeF0v9C3zwvUqpNUqpNampqeX4OOGrPv4YHnrI01kU+u3AbzR4rQHL9y/3dCrCx7l7EfM94AVA5z3+C7jL2Ru11tOB6WBa4G5+XqXbu3evp1MQea6+2s0d33/fDBr3xLj706fNUje7d8OePeYW0b17zeThZ86YJWCyssxjyeVeLsbVKaC97f2V8Rne+Dt8+y3kDZCwilsFXGtdMHZOKfUB4OoIXSHKbOFCcyOPs0ULLmjHDucrJVjl5AkzD+2uXY5bieGl1KplVmBo0QLq1DFz1AYEmM3Pz/XPdnUoqre9vzI+w9t+hwq4HuNWAVdKhWmt81fBGwFsutD7hSiPMWOs7QMvM61NId62zfzPICnJFOeT66EvcMONsDfvvUpBeLiZMevaa81jq1bm9tGoKDMHgIz/Fxa7aAFXSn0B9AdClVLJwN+A/kqpGEwXyl7gvopLUYgKlplpCvO2bWbbvr3wsehquoGBphUdEwrsgocfhvbXmEIdGWlW6hWiEl20gGutb3XysszGJKqes2fNEuob1pmfExLgrhmmn7poP3SzZtC+vZlHpX17s7Vta1rYfn5EntzHCxs+pWWX26FepEd+FSHAC+/ErErS0tIYOXIkq1evZty4cUybNq0gtnbtWsaNG8e5c+cYMmQIb775ptxCX1lycmDnTkhNNf3NI0aYaQp37TLdInkLOnDyJMT2hVtvNQtY5hfqoivmOtGiXgum9J1SGb+JEBckBdxNOTk5BAUF8cILL7Bp0yY2bSp+GeCBBx5g+vTp9OjRgyFDhjB//nwGDx7soWyrsVOnYP16+PNPsyUmmrXR8udct9lMd0hMDNxxB3TqBB3bwaxOMGECuDEO/HzOeY6eOUqjmo0I8pduE+E5UsCBhIQEXn/9dZRSdOnShVGjRvHiiy+SlZVFSEgIn332GY0bN2bq1KkcOnSIvXv3Ehoayueff07v3r2LTTULZkKrjIyMgtvvx4wZw5w5c6SAu2nWLLPcGSdOmOV28ov1n3+ai4v5w7maNDGF+qqrzJCVzp1Nq7pk33Ru+RbUWJW8iv6f9GfRmEVc0fKKch1LiPLwqgL+6PxHSTycaOkxY5rE8MY1b5Qa37x5My+99BIrVqwgNDSU48ePo5Ri5cqVKKX48MMPee211/jXv/4FmK6R5cuXc4nD5NOFDh48SHh4eMHPMoWsG86fNwV65Up6r1wJa9aY8dT5mjeHbt1MP3W3bmYLCyt+jDffhCeeMOMQhaiGvKqAe0L+Ygv5t8w3aNCAjRs3cvPNN5OSkkJWVhYtW7YseP+wYcMuWLwBmULWVVqbm1xWrizc1q0zN7kA3zUaT1jnm7nsvnqmUMfGmqknL2bfPrNYphDVlFcV8Au1lCuK1tqhuD7yyCM8/vjjDBs2jMWLFzN16tSCWM2aNS96zPDwcJKTkwt+lilkS9DajAZZssRsS5cWLhEfHAyXXWbWROvRA3r04IFLmzCkJVw2ybNpC+FtvKqAe8LAgQMZMWIEjz32GCEhIRw/fpz09HSaNWsGmOlfXRUWFkbt2rVZuXIl3bt3JyEhgUceecTq1KsOu92sjlO0YB87ZmLNmpnFK3v3hp49zUVGf5//WgpRJj7/LyU6OprJkyfTr18//Pz8iI2NZerUqdx00000a9aMHj16sKdo32sJkZGRZGRkkJWVxZw5c1iwYAEdO3bkvffeKxhGOHjwYN+7gJmcDAsWwE8/wS+/wPHj5vUWLcxtlf36mS0qqsrdodiqQSveuPoNWjdo7elUhI/z+QIOZqHhsWPHFnvN2aILRbtS8pU26VVcXJzD0MJqTWv4+WeYN88U7a1bzethYXDddaaV3a9fhcwHUarISIiPt/yw4XXCmdhjouXHFcJVUsCFNRYtMtMGBgVB374wfrwZzhcd7bkW9oQJZrPYmawz7EvfR4u6LagZePFrIkJUFCngwhpnzpjHJUssb/X+8APUq2fpIctl9aHVXPHJFfw69lf6R/b3dDrCh8mq9MJaFXABMjbWzCHlsldfNSNZhKimpIALr/f557BsmRs7Hj5c2BcvRDUkBVx4vSefNBMHCiGKkwIuhBBVlBTwcvj555+59NJL6dy5M5deeimLiizftXbtWjp37kzr1q2ZMGGC09vrRdXULqQdHw37iHYh7TydivBxUsDdlJOTQ2hoKN9//z0bN27kk08+4Y477iiI508nu3PnTnbu3Mn8+fM9mK2PatfOjD+3WFjtMO6KvYuw2mEXf7MQFUiGEVK+6WTzRUdHc/78eTIzMzl+/LhMJ+sN7r/fbBbLyMxga+pWOjTsQJ0adSw/vhBl5XUFvP/M/g6vjYoexYOXPcjZ7LMM+WyIQ3xczDjGxYzj2NljjJw9slhs8bjFF/w8q6aT/frrr4mNjaVGjRoynazFFi82i7p7i7WH1jIgYQCLxy6mX2Q/T6cjfJjXFfDKZsV0sps3b2bSpEksWLAAkOlkrdamjZs7/v3vMHu2WU5NiGrI6wr4hVrMwQHBF4yHBodetMVdUnmnk01OTmbEiBEkJCTQqlUrQKaTtdr06eZGnkGDXNzxxAk4cKBCchLCG/j8RcyBAwcye/Zs0tLSAFyaTvbkyZMMHTqUl19+mV69ehW8XnQ6Wa01CQkJTifHEmXz/POmIS2EKM7nC3jR6WS7du3K448/XjCdbJ8+fQq6VpyZNm0aSUlJvPDCC8TExBATE8PRo0cBeO+99xg/fjytW7emVatWcgFTCGE5r+tC8QR3p5OdMmUKU6ZMcXpMn5tO1odEN4pm1shZdGjYwdOpCB8nBVxUX126QAV0XTWq2YhR0aOcxrQ2a1ccOAAHD5p1LQ4eLNwOHYJz58xynzk5ZrPb3c+lPPeHeWpfT362J3/nr7+GK68s3zFKkgIuqq877zSbhc6ehY1Jx1my408uOXkpxw/WZ/9+U7DzH8+eLb6PUtC4MYSHm4uxtWqZSRvzNz+/8uVUngFOntrXk5/tqX3zLqtZyisKuLORIL5GbrUv3erVUGLkZoXJyIA9e2D3bvO4Zw/FCnRaGtAyEcYOghlLYF9fwsIgIsIs5zlkiHkeEWEKdrNm0KQJBARUTv7Ct3i8gAcFBZGWlkZISIjPFnGtNWlpaQQFBXk6Fa/k9gjMZ56BL74wVThPVhYk5f34xx/w1I/Fi3XeYKQCdeqYVeAiIqB7d/N4KhReSYH//AdGdIPAQDfzE6KcPF7A88dMp6amejoVjwoKCip296Yo9MYb0Lo1XHtt2d6flWWK8o5NLdlxeDTb74GkJPNacjLYFfAsfPstzF1pls5s2RLi4sxjVFThY/36jsdftAdeSTDLfUrxFp7k8QIeEBBQ7E5HIUr65z9N10TRAm63mwuCO3bA9u3mMX/bswdycwHuAaDht+Z/AH36mKIcEQn3HoC//hX+Mbj8fdBCeIrHC7gQF2K3m2K8fz+88gps3my27duLXywMDoa2baFbN7jlFvO87Y//S5u5b9Lg6N5ix8zKhXtfhLp1pXiLqk0KuPAKWpvujY0bC4v05s2wZYsp1AsWmC083Cx036+fmS22bVuzNWvmZITA2n1gO2l5rl0bd2XubXPp1KiT5ccWwhVSwEWly801XR3r1pktMdE8Fr2AGBZmCvU998Ann5juj08/Na3mMouPNx3iFgsJDmFwG7mzVnieFHBRobKzYcMG+PPPwoK9YUNh90dgoBl+d/31ZvX5Ll1M4W7QoPAYX35pxlG7VLwBbrvNbBZLPZPK8v3L6dOiD6HBpU+1IERFkwIuLKOB/SmBrJoNq1bBypWmcJ8/b+J16kBMjGlVx8aarUOHi4+R3rrV3PDiLTYe3cgNs29gybgl9G3R19PpCB/mRf8sRFWTlQVr1sDSpbByzmWs4hCHrzXLjNWoYS4oPvCAGT+dP0TP5sb0aXXcXfTmiSfMcvY+PkRVVF9SwEWZnT9vWtVLl8KSJfD772ZeD4A2TWtzJd/Q/a/96XFTBF26WDdG+sUXoX17GDny4u8tJjfX9OEIUU1JARelysoyBfvnn03BXrXKvKYUdO1qukL69TMXGBv+vgiGj4Gb10K3CEvzeO89Mw7c5QIuRDUnBVwU0NrcsbhgAfz0E/z6K5w+bcZKX3opTJhgCnavXs7vUBRCVK6LFnCl1MfAtcBRrXWnvNcaALOASGAvMEprfaLi0hQVJTvbdIl8+y388EPhtCEtW8Ltt8PVV8MVV7gxAqQa6xbWjWV3LqNzo86eTkX4uLK0wGcC04CEIq89BSzUWr+ilHoq7+dJ1qcnKkJGBsyfnzcXyFw4eRKCgsxcxX/5C1x1lbn1vMrr06dCJiupF1SP3s17W35cIVx10QKutV6qlIos8fJwoH/e80+AxUgB92pnz8L335vJ+ebNM33ZISFm/PXw4WbB4BLrNVd9N95oNoulnErh590/c3Wrq2lcq7HlxxeirNztA2+stU4B0FqnKKUalfZGpdS9wL0AzZs3d/PjhDuys+GXX+Dzz2HOHNOf3bQpPPgg3HADXH551ZgLZP9+NyfSz8kxI1Fq1LA0ny2pWxg7ZyxLxy2VAi48qsIvYmqtpwPTAeLi4mTVgkqwdSt88AH83/+ZIdD16pkJnm67Dfr2rRpFuyi3880fB35CLs+I6sndAn5EKRWW1/oOA45amZRw3blzZs296dNh2TJzd+Pw4eZC5DXXWN4IrVRPPWVut7/9dk9nIoR3ceO+OAC+A/KXcR8LfGtNOsIdZ89Cq1Zwxx2QkgKvvWZm9vvyS1PEq3LxBjOJ1ZIlns5CCO9TlmGEX2AuWIYqpZKBvwGvALOVUncD+4GbKjJJcWEnT5rC/eyz8Pzz5V9sVghRNZRlFMqtpYQGWpyLKKfwcCnelSG+WTyJ9yXSqkErT6cifJzciSmqr0GDzFhJi9WuUZuuTbpaflwhXOVuH7gQlSYoyM37ca69Fp57zvJ8kjOSeXf1uxw6dcjyYwvhCmmBC6+3a5ebO54+DZmZlrfCtx3bxkNzH6JL4y40rd3U0mML4QppgYvqa/LkajIngBDOSQEXXu+RR+DDDz2dhRDeRwq48Hr//a+Zi1wIUZwUcCGEqKLkIqYQLro84nJ2TdhFWK0wT6cifJwUcFF9XXstREZaftjggGCi6kdZflwhXCVdKMLrhYZC7dpu7DhoEDz2mOX57D25l1eXv0pyRrLlxxbCFdICF15v/Xo3d0xLM2PBW7SwNJ+k40k8tfApejXvRXidcEuPLYQrpAUuqq/nn4fYWE9nIUSFkQIuvN6dd8K0aZ7OQgjvI10owustWAD+8k0VwoG0wIUQooqSdo0QLurboi9H/nKEekH1PJ2K8HFSwEX1NXIkdO5s+WED/QJpVLOR5ccVwlXShSK8XmQkNGzoxo59+8I991idDknHk5iyaAr7Tu6z/NhCuEJa4MLrrVjh5o4HD5oFQ6OjrUyHPSf28NKylxjcejAt6lk7xlwIV0gLXFRfr74Kffp4OgshKowUcOH1Ro2Cf/3L01kI4X2kC0V4vRUroG5dT2chhPeRFrgQQlRR0gIXwkUDWg7g7DNnCfQL9HQqwsdJARfV1+jR0LOn5Yf1s/lxie0Sy48rhKukC0V4vehoiIhwY8fu3eHWWy3PZ9uxbUycN5HdJ3ZbfmwhXCEFXHi9BQvguefc2HH3bvjjD8vz2Z++n7f+eIvDpw9bfmwhXCEFXFRf//u/MHiwp7MQosJIARdeb+hQeOklT2chhPeRi5jC6yUmQtOmns5CCO8jLXAhXKRQ+Ck/T6chhLTAhXDVoFaDyHkux9NpCCEFXFRjd98NV13l6SyEqDDShSK8Xnw8tG7txo4xMXDddVanw6ajm7jr27tIOp5k+bGFcIUUcOH1vvkGJk1yY8etW2HRIsvzOXTqEDMSZ3D0zFHLjy2EK6SAi+rr3Xfhpps8nYUQFUYKuPB6/frBs896OgshvI8UcOH1kpLgsNy1LoQDKeBCuMjf5k+9oHoyFlx4XLmGESql9gKngFwgR2sdZ0VSQnizAS0HcGLSCU+nIYQl48Cv0Fofs+A4QljrwQfhhhs8nYUQFUa6UITXGzAAOnVyY8cOHeCKKyzPJ/FwIqO+HMWOtB2WH1sIV5S3gGtggVJqrVLqXmdvUErdq5Rao5Rak5qaWs6PE77o009h4kQ3dly3Dr77zvJ8jpw+wpdbviTtbJrlxxbCFeUt4L201t2AwcBDSqm+Jd+gtZ6utY7TWsc1bNiwnB8nhAs++gjuusvTWQhRYcpVwLXWh/IejwLfAPFWJCVEUXFx8OSTns5CCO/jdgFXStVUStXOfw5cBWyyKjEh8qWkwMmTns5CCO9TnlEojYFvlFL5x/lcaz3fkqyE8GJB/kFE1Ikg0C/Q06kIH+d2Adda7wa6WpiLEFVCv8h+7H9sv6fTEEKGEYpq7IknYN48T2chRIWRAi683nXXmQuZLmvZEi67zPJ81hxaw+DPBrM1davlxxbCFVLAhdd7/3247z43dvz9d/j8c8vzOXb2GPOT5pOemW75sYVwhRRwUX199hlMmODpLISoMFLAhddr317qsBDOSAEXXu/UKTh3ztNZCOF9pIAL4aKaATXpENqBS/wv8XQqwsdZMZ2sED6lT4s+bHloi6fTEEJa4KIae+YZWLbM01kIUWGkgAuvd+ut0Lu3Gzs2bWrmBLfYyuSV9P64N5uPbrb82EK4QrpQhNd7/XU3d1y8GHbuhHvusTIdTpw7wYoDKzidddrS4wrhKmmBi+rryy9NN4oQ1ZQUcOH1IiLgXqfrPQnh26SAC69nt4PWns5CCO8jBVwIF9UNqkt8s3hqBtb0dCrCx8lFTCFcdHnE5awav8rTaQghLXBRjf397/Dnn57OQogKIwVceL177oErr3Rjx5AQcwXUYsv3L6fr+13ZeGSj5ccWwhXShSK83tSpbu44fz5s2waPPmphNnAq8xQbjmzgbPZZS48rhKukBS683rlzkJnpxo7ffw8vvWR5PkJ4Cyngwuu1bg0PP+zpLITwPlLAhRCiipICLoSLQoJDGNhyIHVq1PF0KsLHyUVMIVwU3yyeX8b84uk0hJAWuKjGXnsNduzwdBZCVBgp4MLrTZwIQ4e6sWPNmlC/vuX5/LrnV6LejGL94fWWH1sIV0gXivB6f/2rmzt++y1s3AhTpliaz7mcc+w5uYfMXHfGNgphHWmBC6+XmgoZGW7s+NNP8NZblucjhLeQAi68XkwMPPGEp7MQwvtIARdCiCpKCrgQLmpcszHXt7+e+kHWXyAVwhVyEVMIF13a9FK+ufkbT6chhLTARTX21luQnOzpLISoMFLAhdebPBluvNGNHf39ITDQ8nx+3vUzjf7ZiHUp6yw/thCukC4U4fUefNDNHWfNgnXr4JVXLM0nKzeL1LOp5NhzLD2uEK6SFrjwenv2wNGjbuy4ZAl8/LHl+QjhLaSAC6/Xu7fpRhFCFCcFXAghqigp4EK4qFmdZtzR5Q5CgkM8nYrwceW6iKmUugZ4E/ADPtRaW3u1yAdoDbm5kJ0NWVlmc/W5W/3DvqCCRqHENIkhYUSC5ccVwlVuF3CllB/wDjAISAZWK6W+01pvsSq5fOfOmS0np5QtW5OTmWses+xmy9aFP2drsjI12Zl2srM0WVma7CyNPVej0ORkQ0aGeS0nG7KzNTlZGq00/n6arGzNiRO55ORAdq6d7FxzbJufwj9Qk5ltJz0dcuyQbdfk5ECWHZQNlE1xPsfO6cxcsu0Ku4ZsrcnRYM8NIMceCCoHAs6CsoPSoHLN85wgsAeAXybUyCgStwN2yK4FuYHgfw7anuC/CatYO/swdjS5WnNlw6NE1DzP1hM1+PVEKHat0cpOrtbYgeG19xBW4yzrTtdlybkm5B0VO3bsCm7z30io7Qy/ZzVhhYpAa41dmfdoNOMyf6Me51hsa82qGs2xD7gW/cpb2P390dh58N5uXFK7BvO+3kjinnTz+ZjzateaJybEYbPZmDN7PVtSzqC1RitzbJvSPP5gHFprjly6n0W5O7nlqT1oDXY0fv7+tL3yFrTW7Fs5n/QTKWgwx0DjFxhEm2tvQA+9mT2zbuBU5inz2VqTq3OLfb++3vI1R84cKdhXa01Y7TBGdhwJQML6BNLOpmHX9oL4bZ1vo1mdZlZ/1YVwidJau7ejUj2BqVrrq/N+fhpAa/1yafvExcXpNWvWuPxZl91yF2tazQK/bKBIvrk14Hw981pwKtjsJeIBkBNsXquRYYofFL7H7gfa3/zsl1XiF3Q5TeEBCoXW2uG/l9JwSWAwCsX5nPMORdvf5s+cm+cwtO1Qun/YnT8O/lEs3iuiF8vvWg5Ax3c6svXY1mLxBpc0YMP9G6SIi0qhlFqrtY4r+Xp5ulCaAQeK/JwMdHfywfcC9wI0b97crQ+6ol0jtp4IJcc/C7Qq+LdaM6smXfwjsClYkx5Ars2My1WY99TPrk183QhsNhsLT2xDKzsK8jZFGHXp0ywSm03x3wMbQJmYDVAKogJC6dumLTl2O//ZttocVylseZ/QvnZTenboyLlzmXy16TdsysTz39clNIq46M6knzrFD+uXorCZ/fPe1y2iE12iu3AsLY0Fib+ilA2bUtiUDZvNRvd28bRt14kjhw+yZN2vBfvZlMJms9GjQy9aRrYj5dB+Vm5Zip/NVvAeP5uiZ6cBhDVuwaGUPazftRKbsuFnM/valI3u0VcSUq8JKUd3se3AOpSy4adsKJsfSim6tx9InZohHDyaxK4j2/L280PZbChlo0eb/gTXqM2BtN0cOLE3L29/bHnx+OaXExhQgwMn93H49BGUzYYNW8HnX9b0MvxsfiRnJHPszDGUTWHDhp/yw2azEdskFqUUyRnJpJ9PN8dXJh7gF0C70HYAHMw4yOms0wXn3qZsBPgF0Lyu+b4dPn2YzJzMYv9tAmwBNK7VGIC5t80lx57jEM+3cvxKtNbF4oF+gQT6Wd89I4QrytMCvwm4Wms9Pu/nO4B4rfUjpe3jbgtcCCF8WWkt8PKMQkkGIor8HA4cKsfxhBBCuKA8BXw10EYp1VIpFQjcAnxnTVpCCCEuxu0+cK11jlLqYeAnzDDCj7XWmy3LTAghxAWVaxy41nouMNeiXIQQQrhA7sQUQogqSgq4EEJUUVLAhRCiipICLoQQVZTbN/K49WFKpQL7Ku0DK1cocMzTSXgBOQ+GnAdDzkOh8pyLFlrrhiVfrNQCXp0ppdY4u1PK18h5MOQ8GHIeClXEuZAuFCGEqKKkgAshRBUlBdw60z2dgJeQ82DIeTDkPBSy/FxIH7gQQlRR0gIXQogqSgq4EEJUUVLAXaCUukYptV0plaSUespJvL9SKl0plZi3PeeJPCuaUupjpdRRpdSmUuJKKfVW3nnaoJTqVtk5VoYynAdf+T5EKKV+VUptVUptVkpNdPKeav+dKON5sPY7obWWrQwbZsrcXUAUEAisBzqWeE9/4AdP51oJ56Iv0A3YVEp8CDAPs0JdD2CVp3P20Hnwle9DGNAt73ltYIeTfxvV/jtRxvNg6XdCWuBlFw8kaa13a62zgP8Awz2ck0dorZcCxy/wluFAgjZWAvWUUmGVk13lKcN58Ala6xSt9Z95z08BWzFr5hZV7b8TZTwPlpICXnbOFnF29h+np1JqvVJqnlIqunJS8zplPVe+wKe+D0qpSCAWWFUi5FPfiQucB7DwO1GuBR18jHLyWskxmH9i5iw4rZQaAswB2lR0Yl6oLOfKF/jU90EpVQv4GnhUa51RMuxkl2r5nbjIebD0OyEt8LK76CLOWusMrfXpvOdzgQClVGjlpeg1ZMFrfOv7oJQKwBStz7TW/3XyFp/4TlzsPFj9nZACXnYXXcRZKdVEKaXynsdjzm9apWfqed8BY/JGHvQA0rXWKZ5OqrL5yvch73f8CNiqtf6fUt5W7b8TZTkPVn8npAuljHQpizgrpe7Pi78PjAQeUErlAOeAW3TepefqRCn1BeZqeqhSKhn4GxAABedhLmbUQRJwFrjTM5lWrDKcB5/4PgC9gDuAjUqpxLzXngGag099J8pyHiz9Tsit9EIIUUVJF4oQQlRRUsCFEKKKkgIuhBBVlBRwIYSooqSACyFEFSUFXAghqigp4EIIUUX9P9FRtp2PI/9sAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(cyl10.pressure/BAR, cyl10.total_adsorption/(np.pi*100*ANGSTROM**2)/(KILO*MOL/METER**3), 'r', label='cyl10')\n", "plt.plot(cyl20.pressure/BAR, cyl20.total_adsorption/(np.pi*400*ANGSTROM**2)/(KILO*MOL/METER**3), 'r--', label='cyl20')\n", "plt.plot(sph10.pressure/BAR, sph10.total_adsorption/(4/3*np.pi*1000*ANGSTROM**3)/(KILO*MOL/METER**3), 'b', label='sph10')\n", "plt.plot(sph20.pressure/BAR, sph20.total_adsorption/(4/3*np.pi*8000*ANGSTROM**3)/(KILO*MOL/METER**3), 'b--', label='sph20')\n", "plt.plot(car10.pressure/BAR, car10.total_adsorption/(10*ANGSTROM)/(KILO*MOL/METER**3), 'g', label='car10')\n", "plt.plot(car20.pressure/BAR, car20.total_adsorption/(20*ANGSTROM)/(KILO*MOL/METER**3), 'g--', label='car20')\n", "plt.legend()" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3IAAAEzCAYAAAB0TDEBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABZlElEQVR4nO3dd3xb1f3/8dfRsOW97Xhk75BBIISEsCkz7E2ZLbTf0gKFLuiXDlrgxxcKHRRoS4ES9l6lbAJljyQkhEAmWc7y3pa1zu8PKcZJnMRJbF/Jfj8fDyHpzo9uhD7+3HPvOcZai4iIiIiIiCQOl9MBiIiIiIiIyK5RISciIiIiIpJgVMiJiIiIiIgkGBVyIiIiIiIiCUaFnIiIiIiISIJRISciIiIiIpJgPF1ZyBizCmgEwkDIWjvFGJMLPA4MAVYBZ1pra3smTBERkfikHCkiIk7YlRa5w6y1e1trp8TeXwO8aa0dCbwZey8iItIfKUeKiEiv2pNLK08CZsVezwJO3uNoRERE+gblSBER6VFdLeQs8JoxZq4x5vuxaUXW2g0AsefCnghQREQkzilHiohIr+vSPXLADGvtemNMIfC6MWZxV3cQS2rfB0hLS9t3zJgxuxGmiIgkkrlz51ZZawucjqOX7FaOVH4UEemfuitHdqmQs9aujz1XGGOeBaYCm4wxxdbaDcaYYqBiO+veDdwNMGXKFDtnzpw9jVlEROKcMWa10zH0lt3NkcqPIiL9U3flyJ1eWmmMSTPGZGx+DRwFfAG8AFwYW+xC4PnuCEhERCRRKEeKiIhTutIiVwQ8a4zZvPwj1tpXjDGfAk8YYy4G1gBn9FyYIiIicUk5UkREHLHTQs5a+zUwqZPp1cARPRGUiIhIIlCOFBERp3S1sxMREdlKMBikvLwcv9/vdCiO8fl8lJWV4fV6nQ5FRETihPJjVE/nSBVyIiK7qby8nIyMDIYMGULs0rp+xVpLdXU15eXlDB061OlwREQkTvT3/Ai9kyP3ZEBwEZF+ze/3k5eX12+TlDGGvLy8fn/GVUREttTf8yP0To5UIScisgf6c5ICfX4REemc8kPPHwMVciIi/cBFF13EU089BcAdd9zBiBEjMMZQVVXVvoy1liuuuIIRI0YwceJE5s2b51S4IiIivaZjjjz33HMZPXo048eP57vf/S7BYBCIzxypQk5EpJ+ZMWMGb7zxBoMHD95i+ssvv8yyZctYtmwZd999N5deeqlDEYqIiDjj3HPPZfHixSxcuJDW1lbuueceID5zpAo5EZEE9sADDzBx4kQmTZrEKaecwtChQ9vPHjY0NDBkyJD295tNnjyZIUOGbLOt559/ngsuuABjDNOmTaOuro4NGzb0xscQERHpdruTI4877jiMMRhjmDp1KuXl5UB85kgVciIiCWrRokXceOONzJ49mwULFnDvvfdy6KGH8p///AeAxx57jNNOO63L3R6vW7eOgQMHtr8vKytj3bp1PRK7iIhIT9rTHBkMBnnwwQc55phjgPjMkRp+QESkG/zu34v4cn1Dt25zXEkmvz1hr+3Onz17Nqeffjr5+fkA5Obmcskll3DLLbdw8skn869//Yt//vOfXd6ftXababpZXURE9oQT+RH2PEf+8Ic/5OCDD+aggw4C4jNHqpATEUlQ1tptksiMGTNYtWoV//3vfwmHw4wfP77L2ysrK2Pt2rXt78vLyykpKem2eEVERHrLnuTI3/3ud1RWVvKPf/yjfVo85kgVciIi3WBnZwZ7whFHHMEpp5zCVVddRV5eHjU1NeTm5nLBBRdwzjnn8Otf/3qXtnfiiSdyxx13cPbZZ/Pxxx+TlZVFcXFxD0UvIiL9gRP5EXY/R95zzz28+uqrvPnmm7hc39yFFo85UvfIiYgkqL322otrr72WQw45hEmTJvGTn/wEiPa4VVtbyznnnNPperfffjtlZWWUl5czceJELrnkEiB6g/ewYcMYMWIE3/ve97jrrrt67bOIiIh0p93NkT/4wQ/YtGkT06dPZ++99+b3v/89EJ850nR2vWdPmTJlip0zZ06v7U9EpCd99dVXjB071ukwtvHUU0/x/PPP8+CDD/bK/jo7DsaYudbaKb0SQB+g/CgifUm85kfoWzlSl1aKiPQhl19+OS+//DIvvfSS06GIiIjElb6WI1XIiYj0IX/961+dDkFERCQu9bUcqXvkREREREREEowKORERERERkQSjQk5ERERERCTBqJATERERERFJMCrkRET6gYsuuoinnnqq03nHHHMM2dnZHH/88VtMX7lyJfvvvz8jR47krLPOIhAI9EaoIiIivWp7OXL+/PlMnz6dvfbai4kTJ/L444+3z4uHHKlCTkSkn/v5z3/e6Xg6V199NVdddRXLli0jJyeHe++914HoREREnJGamsoDDzzAokWLeOWVV7jyyiupq6sD4iNHqpATEUlQzc3NzJw5k0mTJjF+/Hgef/xxhgwZwtVXX83UqVOZOnUqy5cvb1/+nXfe4YADDmDYsGFbnHk84ogjyMjI2GLb1lpmz57N6aefDsCFF17Ic8891yufS0REZE91R44cNWoUI0eOBKCkpITCwkIqKyvjJkeqkBMRSVCvvPIKJSUlLFiwgC+++IJjjjkGgMzMTD755BMuu+wyrrzyyvblN2zYwHvvvceLL77INddcs8NtV1dXk52djccTHW60rKyMdevW9dhnERER6U7dnSM/+eQTAoEAw4cPj5scqQHBRUS6w8vXwMaF3bvNARPg2P/b7uwJEybws5/9jKuvvprjjz+egw46CIBzzjmn/fmqq65qX/7kk0/G5XIxbtw4Nm3atMNdW2u3mWaM2Z1PISIi/ZkD+RG6N0du2LCB888/n1mzZuFyueImR6qQExFJUKNGjWLu3Lm89NJL/PKXv+Soo44CtkwmHV8nJye3v+4sCXWUn59PXV0doVAIj8dDeXk5JSUl3fwJREREekZ35ciGhgZmzpzJDTfcwLRp04D4yZEq5EREusNOzgz2hPXr15Obm8t5551Heno6999/PwCPP/4411xzDY8//jjTp0/frW0bYzjssMN46qmnOPvss5k1axYnnXRSN0YvIiL9ggP5EbonRwYCAU455RQuuOACzjjjjPbp8ZIjVciJiCSohQsX8vOf/xyXy4XX6+Vvf/sbp59+Om1tbey///5EIhEeffTRnW7noIMOYvHixTQ1NVFWVsa9997L0Ucfzc0338zZZ5/Nr371KyZPnszFF1/cC59KRERkz3VHjnziiSd45513qK6ubi8E77//fvbee++4yJFmZ5fXdKcpU6bYOXPm9Nr+RER60ldffcXYsWOdDmMLQ4YMYc6cOeTn5/faPjs7DsaYudbaKb0WRIJTfhSRviQe8yP0vRypXitFREREREQSjC6tFBHpQ1atWuV0CCIiInGpr+VItciJiIiIiIgkGBVyIiIiIiIiCUaFnIiIiIiISIJRISciIiIiIpJgVMiJiPRD1dXVHHbYYaSnp3PZZZdtMW/u3LlMmDCBESNGcMUVV9Cbw9SIiIg47fXXX2ffffdlwoQJ7LvvvsyePbt9XjzlSBVyIiL9TCgUwufzcf3113PrrbduM//SSy/l7rvvZtmyZSxbtoxXXnnFgShFRER6XygUIj8/n3//+98sXLiQWbNmcf7557fPj6ccqeEHREQS2AMPPMCtt96KMYaJEydy5plncsMNNxAIBMjLy+Phhx+mqKiI6667jvXr17Nq1Sry8/N55JFHOPDAA1m+fPkW29uwYQMNDQ1Mnz4dgAsuuIDnnnuOY4891omPJyIistv2JEduttdee+H3+2lra6OmpiaucqQKORGRBLVo0SJuvPFG3n//ffLz86mpqcEYw0cffYQxhnvuuYdbbrmF2267DYheDvLee++RkpKy3W2uW7eOsrKy9vdlZWWsW7euxz+LiIhId+quHPn0008zefJkkpOT4y5HqpATEekGN39yM4trFnfrNsfkjuHqqVdvd/7s2bM5/fTTyc/PByA3N5eFCxdy1llnsWHDBgKBAEOHDm1f/sQTT9xhEQd0eq2/MWY3P4GIiPR3TuRH6J4cuWjRIq6++mpee+01IP5ypO6RExFJUNbabRLI5ZdfzmWXXcbChQv5xz/+gd/vb5+Xlpa2022WlZVRXl7e/r68vJySkpLuC1pERKQX7GmOLC8v55RTTuGBBx5g+PDhQPzlSLXIiYh0g52dGewJRxxxBKeccgpXXXUVeXl51NTUUF9fT2lpKQCzZs3a5W0WFxeTkZHBRx99xP77788DDzzA5Zdf3t2hi4hIP+FEfoQ9y5F1dXXMnDmTm266iRkzZrRPj7ccqUJORCRB7bXXXlx77bUccsghuN1uJk+ezHXXXccZZ5xBaWkp06ZNY+XKldtdf8iQITQ0NBAIBHjuued47bXXGDduHH/729+46KKLaG1t5dhjj1VHJyIiknD2JEfecccdLF++nOuvv57rr78egNdee43CwsK4ypGmN8c+mDJlip0zZ06v7U9EpCd99dVXjB071ukwHNfZcTDGzLXWTnEopISj/CgifYny4zd6MkfqHjkREREREZEEo0JOREREREQkwaiQExERERERSTBdLuSMMW5jzGfGmBdj73ONMa8bY5bFnnN6LkwRkfjUm/cZx6P+/vlB+VFEpDPKDz1/DHalRe7HwFcd3l8DvGmtHQm8GXsvItJv+Hw+qqur+22ystZSXV2Nz+dzOhSnKT+KiHTQ3/Mj9E6O7NLwA8aYMmAmcCPwk9jkk4BDY69nAW8DzgwUISLigM0Dg1ZWVjodimN8Ph9lZWVOh+EY5UcRkW0pP0b1dI7s6jhyfwZ+AWR0mFZkrd0AYK3dYIwp7GxFY8z3ge8DDBo0aPcjFRGJM16vl6FDhzodhjjrzyg/iohsQfmxd+z00kpjzPFAhbV27u7swFp7t7V2irV2SkFBwe5sQkREJO4oP4qIiJO60iI3AzjRGHMc4AMyjTEPAZuMMcWxs43FQEVPBioiIhJnlB9FRMQxO22Rs9b+0lpbZq0dApwNzLbWnge8AFwYW+xC4Pkei1JERCTOKD+KiIiT9mQcuf8DjjTGLAOOjL0XERHp75QfRUSkx3W1sxMArLVvE+19C2ttNXBE94ckIiKSWJQfRUSkt+1Ji5yIiIiIiIg4QIWciIiIiIhIglEhJyIiIiIikmBUyImIiIiIiCQYFXIiIiIiIiIJRoWciIiIiIhIglEhJyIiIiIikmBUyImIiIiIiCQYFXIiIiIiIiIJRoWciIiIiIhIglEhJyIiIiIikmBUyImIiIiIiCQYFXIiIiIiIiIJRoWciIiIiIhIglEhJyIiIiIikmBUyImIiIiIiCQYFXIiIiIiIiIJRoWciIiIiIhIglEhJyIiIiIikmBUyImIiIiIiCQYFXIiIiIiIiIJRoWciIiIiIhIglEhJyIiIiIikmBUyImIiIiIiCQYFXIiIiIiIiIJRoWciIiIiIhIglEhJyIiIiIikmBUyImIiIiIiCQYFXIiIiIiIiIJRoWciIiIiIhIglEhJyIiIiIikmBUyImIiIiIiCQYFXIiIiIiIiIJRoWciIiIiIhIglEhJyIiIiIikmBUyImIiIiIiCQYFXIiIiIiIiIJRoWciIiIiIhIglEhJyIiIiIikmBUyImIiIiIiCQYFXIiIiIiIiIJRoWciPRNbY2w7A2oXuF0JCIiIiLdzuN0ACIi3W7lu/DkhdBSjcUQOPiXJB9+tdNRiYiIiHQbtciJSN+yfj48fAaR1HzuGXQLL4Snk/zO/+OJR+4hErFORyciIiLSLVTIiUjfEWiBZ76HTcnhx77ruXFZGQun/B8bkwaz3+JbuO/d5U5HKCIiItItdlrIGWN8xphPjDELjDGLjDG/i03PNca8boxZFnvO6flwRUR24P2/QNVSXhj6K/69PMQNJ4/nVydNoujE3zHUtYl5bzxORaPf6Silj1B+FBERJ3WlRa4NONxaOwnYGzjGGDMNuAZ401o7Engz9l5ExBmNm+CDv9Iy4gR+PjeXmROKOXf/wQCYsccTSivmVN7gsU/WOhyo9CHKjyIiDvEHw06H4LidFnI2qin21ht7WOAkYFZs+izg5J4IUESkS969FcJt/NmejTHw6+PHfTPP7cUz4RQOdn/Bcx8t0b1y0i2UH0VEep+1lj+9vpQJ173KW0sqnA7HUV26R84Y4zbGzAcqgNettR8DRdbaDQCx58Iei1JEZEeaq2DeA9SNOp1/fmm4aMYQBmT5tlxmzPEkEWRM88fML69zJEzpe5QfRUR61zvLqvjLm8twGcN1LyyiLdR/W+a6VMhZa8PW2r2BMmCqMWZ8V3dgjPm+MWaOMWZOZWXlboYpIrIDn94DIT9/aTma9GQPlx4yfNtlBu6PTUpnhvtL3vxqU+/HKH2S8qOISO964tO15KR6ue3MSayubuGTlTVOh+SYXeq10lpbB7wNHANsMsYUA8SeO23btNbeba2dYq2dUlBQsGfRiohsLdACn9xNy9CjuH9ZMudPG0x2atK2y7k9mEHTOTh5KbMX649m6V7KjyIiPa+uJcDrX27i5MmlHDQi+ru5cF29w1E5pyu9VhYYY7Jjr1OAbwGLgReAC2OLXQg830Mxiohs3xdPQUs1j3tPwuMyXHjAkO0vO/gABobWULGxnPrWYK+FKH2T8qOISO96b3kVgXCE4yeWkJXqZWBuCovWNTgdlmM8XVimGJhljHETLfyesNa+aIz5EHjCGHMxsAY4owfjFBHp3NxZhPNGc+viPE6YWExRpm/7yw45EIApZgnz1tRy2GjduiR7RPlRRKQXvb+8ioxkD5PKsgAYX5LFF+v7b4vcTgs5a+3nwOROplcDR/REUCIiXbLxC1g3hzmjfkbzugjfPXDojpcfMBFr3Exwr2LuKhVysmeUH0VEetd7y6uYNjwPjzt6UeH40ixe/mIj9a1BslK8DkfX+3bpHjkRkbgybxbWncQtG/ZmUlkW40uzdry814cpGM3+vrXMX1vXKyGKiIjIniuvbWFtTSsHDM9rnzaqKAOAryubtrdan6ZCTkQSU7AVPn+c2sHHMLfSxTlTB3VtvQETGW1X8tWGBqzVeHIiIiKJYPMJ2H0H57RPK44NNbSpwe9ESI5TISciienL58Ffz2Phw0lLcnPCpJKurVc8kcxQNaa5korGtp6NUURERLrF/DV1JHlcjElvhfdvh8ql7ffFb6xXIScikjjmziKcM5Tbvy7ipMmlpCV3pe8mYMAEAMa5VvHl+v7b05WIiEgimb+2jr2LU0h65DR4/dfwr2PJ8/jxug0bG/rniVkVciKSeCqXwpoPmJ9/Iv6g5dtdvawSoGAMAMPNer7coEJOREQk3gXDERauq+f81I+gYhEc9FNoqcK18AkKM3y6tFJEJGHMm4V1efjDpn2ZUNqFTk46SisAXxbjkytYUdE/b44WERFJJEs2NtIWCnNI7VPRK2sO/zXkDoOlrzAgy6dLK0VEEkKoDRY8Sv3Ab/FRhYezpw7ctfWNgfxRjPFuZEVVc8/EKCIiIt1m/to6RplyMhuWwb4XRXP5iCNh1fsUZ3rVIicikhAW/wdaqnmKI0hNcnPiDjo5idgIraHWbWfkj2JguJyvK5rUc6WIiEicm7+2jtN8c7AYGHNCdGLJZAi1sldSBRvq/f0yn6uQE5HEMm8Wkcwy/vR1GSdMLCHDt+0AoK2hVv44548c+OiBTH14Kt955TusaVjzzQL5I8kMVUNbA5XquVJERCSuLVhbx3GeuZhB0yGjiLWNa6nPHw7AGPs1rcEwDf6Qw1H2PhVyIpI4albC12/zReEJNAdtp5dV1rfV873Xvsf9i+7nwNID+Z+J/8OyumV859XvUNVaFV0obyQAQ80GVlTq8koREZF41egPUle5joHBr2HUUTzy1SMc98xxHP/eT1mblExpcBVAvzwxq0JORBLHZw+CcfGnqv0ZMyCDvQdmbzE7GAny47d+zJfVX3LbobdxyyG3cNnky7j3qHup89dx08c3RRfMGQLAQFPJSt0nJyIiErc+L69nulkEwLoBY/nj3D+yV95etIZaubtgANn+dQDUtwadDNMRKuREJDGEQ/DZwzSWHcpbG5M4e7+BGGO2WOSWT25h7qa5XD/jeo4cfGT79NG5o/nexO/x2urXWFi5ELKjwxUMcVeytralVz+GiIiIdN38tXVMdy3CJmfySM3nhCNh/nzYnzl26LG8kWTwtURvnahvDTgcae9TISciiWHZa9C0kRfc3yLZ4+KUyWVbzH5v3Xs8tuQxLhh3ATOHzdxm9fPGnke6N50Hv3wQfJmQksOo5BrW1qiQExERiVfz19ZxiPdLAoNn8PzXL3D4oMMZkDaAQ8sOpYkIS9rWA5a6FrXIiYjEp3mzsGlF3LpyCDMnFJOV+k0nJ02BJq774DqGZQ3jin2u6HT19KR0Thl5Cq+vfp06fx1kD2aop1qFnIiISJyy1lKxZgkldhPvFAykvq2e00adBsDkoskALHKFyKBVhZyI9LJQG9StgdZapyOJbw3rYdlrLB5wPLVtcPbUQVvM/su8v1DZWsn1M64n2Z283c0cP+x4QjbE7LWzIWcwxZFNrK3tZHgCERERcdz6ej+jWz8D4B3jJyMpg6kDpgKQ68ulNCmbL5KTKDS11PXDe+Q8Tgcg0u9YC8teh4//DivfgUjsh6dwHOz/P7D3eeDW/5pb+OxhsBH+WncAwwvS2G9ITvus5bXLeWLpE5w56kwmFkzc4WbG5o5lYMZAXl31KqdmDyY3+BK1rX4a/cFOhzEQERER58xfU8cBrkW0pRTwbvVCZpTMwOP65m+kkekD+bqpgqHJjTT0w0JOLXIival+HTxyFjxyBlQuhmk/gBNuhyN+C55k+PeP4f7jostJVCQCnz1AS+kMXlqXwtn7DWrv5MRayx/m/IE0bxo/2vtHO92UMYbDBh7GnI1z8GeV4rZBCqmjXK1yIiIicWfB2loOcH3JksH7Uu2v5uCyg7eYPyhrKGs8HgYl1VPXos5ORKSnrHof7j4EVr0HR90IP14AR90A+14IB/0EvvcWnPpP2PQl3Hc0VK9wOuL4sPwNqFvDf7xHkeR2ceo+pe2z3l33Lh+s/4BLJ11Kti+7S5ubVjyNQCTAPFcYgIGmgo31/p6IXERERPZA5crPKTR1vJ+VjcEwo3TGFvMH542mzeUiO6miX15aqUJOpDcsfAoeOBF8WfC92XDAZeDe6lI+Y2DimXDRixBohlknQsMGZ+KNJ5/cTSStkBu/HsHxE4vJS4/eAxeMBLl1zq0MzhzM2aPP7vLm9i3aF4/Lw0f+6LEtMdVsUCEnIiISV4LhCLkVHwLwTqCSCQUTyPXlbrHMoNxRALi8FersRER6wPxH4JnvwcBp0SKucEz7rPqWIP9dWsmTc9by0sINrKluwRZPgguei3aA8uhZ0aKuv6peActfZ0HRqdQFDBceMKR91hNLnmBl/Up+NuVneLcuincg1ZvKhPwJzKuPtngWu2rY2KBCTkREJJ58ub6BKXYRa9NL+KJuGQeXHrzNMoMzBgPgd9f2ywHB1aOCSE+aez/8+0oYdgic/SgkpWKt5YMV1dz33kreWlJBxG65yoTSLC4/fARHnfGv6P10L/0cTr7Lieid98k/sS4v12+cxuRB2UwamA1AfVs9f1vwN/Yv3p9Dyg7Z5c1OzJ/Io4sfJZiUwRBXA/PVIiciIhJX5qyq5lTXV7xeMhWCy7e5Pw6gMLUQlwW/q7Ff3iOnQk6kp3zyT3jpZzDiSDjrIfD6WFHZxHUvLOLdZVXkpyfxP4cM56CR+ZRmp9DoD/Hxyhoe/mg1339wLseOH8BfZvyUpPdvhSEHwt7fdvoT9a62Jpj/MBUDj2bekiT+ctSQ9ll/X/B3GgON/HzKz9s7PtkVEwomMOvLWSzJLmJQXR2vqEVOREQkrmxaNocc08RHqR4KA4WMyR2zzTJul5t846XRBKhvDRKJWFyuXf+7IFH1biHXuBH++wdwucC4weXu8OyKPlxucCeBNwW8qVs9p4AnBZLTISk9ek+RdL9QABrKobkK2hqjDxsGlzd6X1dyJqQXQlpB9J4v/Tts64M74LVrYfRMOONfWHcSD320mutf/JJkj4vfHD+Oc6cNItnj3mK18aVZXDB9MP9892tue20pJ1UcyLOl7+P7z0+hdF8oGO3QB3LAgkehrYF/+I+kICOZY8cXA7CqfhWPLX6MU0eeyujc3TsekwomRXeRlsG0+hp1diIiIv2avy3A0vKN5OcVUJKd4nQ4WGtJKf+AIPBh8xqOHnrsdk/cFrpTqHHVELHQFAiR2Y+GE+rlQm4DvHVD92zL5QFfNqRkQ0rOlq9T8yGjCNKLogVH+oBo0eFJ6p599xXWQu0qWD8P1s2D9fOhdmV08GXsTlaOScqAglGQPzp671fZflAyOVp091fv/hHe/B2MOwlOu5e6NsvVT8/l1UWbOGRUAbeeMYmCjO0PWu11u/jhoSOYWJrNDx+ey5n+i3nW/QvcT30XLnkTvL5e/DAOiYTh47/Tmj+B+1bn89MjB5Pkid7Se9uc20j2JHdpuIHtKUotIteXy+Kwm5kR3SMnIiL9VCTM4qd+T8mX/2QizSyMDOHJUb/hB2efvM3J5t5UXtvKhOAC3s8dRFOohYNKD9rusgXedMpd1YClviWoQq7HlEyGX38Y/SPNhjs8R6LPNhKdFm6DoB+CrRBs+eY55I8+tzVCa120Mwh/XfR1cyVUL4tNq+98/ym5kDEgWtxllkJWWey5FDLLos/JGb14QHpZU2WsaJsbe8yD1proPHcyDJgAQw+G7EGQPTh6nJIzoq2fbi+Eg9HBq/310FQBTZugbk10PLQVs2HBI9FtuTxQPAkGz4CRR0Y7+egPRbS18PZN8N+bYcKZcPLf+GRNA1c+9hmVTW38auZYvjtjaJeb/A8cmc9j35/Ot+/5iKvdl3Lrphvhjd/CsTf38AeJA4tfhOrlPFbyW9KSPFwwfQgAH67/kLfL3+bKfa4kPyV/tzdvjGFkzkiWVS0lM1hFo7+N1kCYlCTnkpaIiEivCodYc/eZjNn0Jh8lTSNn5DQGL32Ii5Zdxh8fNlxzwSm7dftCd5i3spLDXIv5a8EkPLaCaSXTtrtsYXIO89xryaSF5kCoF6N0Xu/fI+f2btvtencLBaA5VmhsLjgaN8Xexx4r3oKmjdHisSNf1jdFXXuh16HgyyhJjBaR1tpoC9v6z2DDfFj3GdSvic4zLigYC2OOg5J9opfsFY7b82KruQrKP4W1H8Oaj+Gjv8EHt0cLwWGHwqijYczxkJq7000lHGujRdb7f4HJ5xOe+WfueGslf3lzKQNzU3n60gOYWJa9y5sdV5LJQxfvzzn/tByQdAKnfvx3GHYYjD6m+z9DvLAW3vszwawh/L9Vo7hoxiCyUr2EIiFu+fQWStNLOW/ceXu8m5HZI3lq41wsYfKoZ2ODn6H5ad3wAUR2T1so7OgZcBHpR6yl/LErGLTpTR7N+QGn/uj/RX9/6i6m5c5DOHnFb3n6k0mcvv9wR8Jbv+RjMk0rn3gDTMmZQpp3+/m5MCWferebbFcdzW0q5HpMZWMbN738FWfsO5ARhek9tyNPUrT4yirb8XLhYPS+vfpyaFgXfW5/vRbK53zTYtVRSk60sMsohsyS6COjODotMzbNl90794611kW7aK9aGm2RrFoKGxdGL5ncLGcIlO0L+38/WrQNmAjJ6QTDESoa29hY38qGL6vYWO9nQ72fjfV+NjX4aWoL0RwI0dIWpjkQIhKJ3t7oMga3MaQmu8lK8ZKV4iXT56UwM5nS7BGU5E6gdNjlDEqPMKD6Y8zyN6KDOi9+EV68Ktrqt9cpfaeoCwfhPz+FebNgv0tYd8DvueqeT/lkVQ0n713C9SePJ2OrZv5gJMiCigV8UfUFqxtX0xxoxuv2UpZRxr6F+zK5aDJeV3Sd8aVZ3H3+FC65r5V9UhYx+PkfYn7wfvS71hetfAfWz+PlgT+HKhcXHzgMgGeWPcPyuuXcdshtJLu3f2lqV43KGYXfhij3eCg20fvkVMiJU/zBMIfc8jZ//fZk9hvSB34XRSSubXrtNsqWPcxTvlM58Qc3Eoi0sLx6DSVpJWSddhdjHzubV1++lfoJt5OV2vuXKnrXvMdaj5uv/VWcXvbdHS5bkFYEQKangqa2cG+EFzd6tZALhCLc8+5K/vHfr9l7YDbHjh9AUaaPlkCY2pYAtc3RHmfCNnZ/lu38Tq2ItYQilnDYEraWzYtvrptcBtwuEy04XNGiw9Xx2cUW09yubFyuHNxmAq4kg6fQkFTiwudxke4OkBOsJDNQQUZgE2ltlaT4N+Fr3URS/UY86z/D1VK1bZAub7RIScmF1DxIzYm9zo22UHXsvMXri17aCNEWQhuJfvhIKNpz3+YOR9oaoaUqeq9hw4boc6Cpwz49kDMUO2AirRPOoypzHOXJo1gfSGFTg5+KSj8blvvZ2DCfDfV+qpra2o/dZileN8XZPooyfAzKTSUt2UNqkpvUJDdulwtrLeFI9Pi3BEI0tIaobw2yod7PgvI6qpq27Po1PdnLyKIzGD3wO0wfvY79Wv7LgPJXcL1weaxb/kNh/Gkw9vhoa2ii8TfAkxdGLy096Ke8VHAJ1/zlPcIRy21nTOK0fbc8mbCoahGPLH6EN9e8SXMwOj5cTnIOmcmZtIZaqVhRAUCuL5ezR5/NeePOIyMpg+nD8/h/Z+zHxY9fykvhX5P07P9gzn8uWln3Ne//mXBqAb9aNYFTJ5cxIMtHY6CRO+ffyT6F+3Dk4CO7ZTejcqKDiC5L8jKgpYZNuk9OHJTkcZHkcfHTJxbw6pUH6zJfkUTW1hS99aetKfo3XVJa9G+clGynIwOgfs5TFHx4A2+a6ex7yR+4bd7/49nlzxKKhHAZF6eMOIXLBh3Od1c/y32vn8+VJx3Qq/FtavCzV/OnvFQ0CAhzaNmhO1w+Ly16YjvVU60WuZ5UmpPCS788nOc/W89Tc8u56eXFW8xP8UZbeNwu016UGQOGLVu2TKxQ88SKNZcx7QWfjRV2YWuJRKKFXjjS8XW0ENxymv1m2nb7+EgGBsUeW/ISosjUUuauo8xTS5m7lkJXIzmBJnKCTWTX15Fh15IRaSAt3ICb3TtbEHCn4vdk0ZhUQKN3INWZe1Np8lhjilkRKWF5MJ/aJkvlxjaC4c0fZGn7+hk+D8VZPgZkpTB2QCYDsnyx9z6Ks1IYkOUj0+fZo+uh/cEw6+paWVfbyuqaFpZvamTJpkZe+6qCx5otcDBu18Ecm1fB6b5PmbL+bdJX/BD74pWYkUdFi7pRx0BS6m7H0GuqlsMT50PlEpqO/hO/WbMPz7z+GZMGZnP72XszOC/auhMMB3l99es8svgRFlQuIMWTwrFDj+Xg0oOZMmAKWcnfFLAtwRY+3PAhzy57lrsW3MXjSx7nZ/v9jJlDZ3LS3qWU1x7Ob95Yys0r/wkf/AUOvMqpT98z1nwMK2Yzu+RSmuu8XHpo9JKOuz+/m1p/LXd9665uu15/WPYwDIZlSV6KTC0b1HOlOMhlDH84fSJn3f0Rf35zKb88dqzTIYnEN2ujV8TYSPREtsvdu71ohwLRq5+ql0P1MiJVywlsWoqpWUGyv7LTVfyeLFozBhMp2Zfs0QfhHnYwpBf0XsyAf+WH+F68lAV2JN5v38LPP/gBy+qWccaoM5hWPI1PN37Ko4sfZV3heP66po2cuX9h/SGTe7Uny4++XMVxriX8I3siw9KzGJg5cIfLZ2eUAuBz19OkQq7nRGyENF+Yc6YVcva0Amqa26hpbiMt2U2mz4PXY7CxksxaS8RGtnxPJFqoxZrqOr63nbTdbVMAdni/9R+Dm+dF9wsRawiFXYRCEAib6HMIQmFDIAxtAfCHIviDYVoDYVqDsUcgzIZAmK9j7/3BMC2B6PTosiGCgTZsqJWkiB+vDZBCG0mEsIDFRQRDBEMYF802hSZSaMZHhG9aX5I8LjKSPWT4PKT7PGQkeynN9DDG56UgI5nCjGQKM5MpyvRFX2f48HldtIXb8If8+MN+WkOt+EPN+MPVVIdaKa9sjc7bYr6ftnAb1lpcxoUxJvqMwe1yk+JJIc2TRqo3NfrwpJKVnMXo0jymDy/B2+F+yMrGNj4vr2PB2jo+W5vPFWtLaPAfx95mBWckf8TMZR+SvfhFwp5UzJiZuCacBsOPiM+OUhY+Bf/+Mdbt5f39/87lr2fT6F/PFYeP4PIjRuJ1u6hqreLJpU/y5JInqWytZFDGIK7e72pOGnESGUmdd6qT6k3liEFHcMSgI1hUtYgbP76RX777S/679r/8Zvpv+OGhw7m29gJe/OxzjnvzBlxDDoKyKb384XuItTD7esKpBfx09VTOnFLGkPw0ltYu5aEvH+LkESezV95e3ba7FE8KgzIGsqy5hXGeOsrVIicO239YHmdNGcg9767kpEmljCvJdDok2RPWxq6ysd9cabO5Uze7uZM3u9X7jvMjnSy/+b3tZPnN7yPb2V74m1i6tM9OOqXrOH2bZb5Z1kZC2Eh4iwebX9tvXndc10bC2+zbREKYSCD2HIw9oq9ddtuT4hYX1uXGGg/W5SbiSSHiScV606IniL1pkJSGSU7DlZyOKzkNt9eH8fjAkwweX/RvDndS9KqocCBaLAZbsc1VtDVsItSwCVf9WnxNa7eIodZmssIWszKyF6vsAGrIIOhOxeNxk0obaaF6iv0bGBFYz8SaR3Av+hcRXKzLmAjjTqD0gG/jyirp0a9kcNMSgg+eSbXNYcGRv+L+z66gLdzGnUfcyYGlBwLwrcHfYmzeWH79/q+5deRUfrH0DW596V3+99tH9WhsHdUsfI02V4QFkTrOH3jiTpfPyYoWeknuBrXI9aSvar5i2iPb73Um0biNO/pwufEYD26XG5dxRV973biTovPcxk2ycZHm8myxjsu4cJvosws3LuPGGDcuXBii8zxuDx7jxuNyY4m0P8I2HL3M0YaJ2Ej7c32ojU1hP/4aP/7Kb4qytlAb/vCu/6FqMCS7kzHGtBfXmwvocCc/olvLTMok15dLXkoexWnFlKSXMGRoKdPHF1OcOoY2fwafl0/ik5WHcs/KSgbUzeOE0Icct/Blsr94kjZPBo3DjiN76jl4hh0cPePmpKZKePV/YeETNBTsy1Why3nz7SQmD0rj/06dyKiidD6v+pzHFj/GK6teIRQJMaN0Br8b8ztmlM7AZbp+KeRe+Xvx0HEPcd8X93HHZ3ewsGohtx9+O78/aTxX1v6SvddcQu6jF5J6+Yfg6wN/8K38L6x6lxeLr8Bfn8Llh48kYiPc8NENpCelc9W+3d/6OCRrKKtrV3NMUj1z1CInTqpfC/++kt+7w+yfvIG1D8xizLgiXLGTlPWtQepbAiR53GSluEnxutvvPQiGoycOXUCy10VspA5CYUsgFCYUsdiIjQ7hSrTRwoXFmOgpTJfZfPULm09pxi67/+ZEavSZ9g7CoidRN0/8Ztktr9ePbWfrZbBb7Se2j/bbCr4peqL73rYgMu3TO06LbA5yy+Kp47yO24stZzq87mw7hq32scV0C0Ta4zHWYtiqEzWHxD4Bm4/C5mmR2FVMHR8Rs3kZQ6TjPLN5G9F1wrgIGRdhXIQxBHERsZ1NM4SMq/2kdIhvXofNttNCsRPYEVyEjZswXsJACHfs4Yo+G0PI7SLkiU2zLiwGYyK4zOYtfPPwEiQ5FCIlGCKlpYFUakizAVIJkkaANBsghRBeIrgteLB4Ys9uCx5r8RB9brbJlJsM1rkyWOvKZpV3KOvcGbSmZRJOT8WbavF4/VhXMwFbQ3NoFfVtdbSEWkhyJZHmTSM9KYOcpHEk24Pw1UcYWrmBb9UtZPLH1xP++Ea+ypyGd8qFDD/gVEw3n8QO16+n4Z4TsWGYNeUHvLDy9xSmFnLf0fcxLHvYFsuePOJkltQs4aGvHuLANB/Dv7qLVVUzGNIL95Fba8le/1/eSssmZHd+WSVAdka0kHN5mvpdIWfsFj+6PWvIXkPsdY9fF0seZvvPm193eL/5D+DN07Z+v3mdzbb+XB1b7LZuvets2XAkTNiGCUfChGzom/exaWEbJhQJdf4+9jpiI1uu12HZza8jNtK+/c2vI5HIFsuEbAhro0l3c+HXXgDGHpvf+zw+fG4fyZ5kUtwpJHuS8bl90ekeH8nuZFI8KaR4UraY3vF9x9ebi7jORGwEf8hPS6iFlmALLaEWmoPN1LXVUeOvobq1murWamr8NVS1VrGxeSMbWzYS6dBTqMEwIG0AQzKHMCRrCLlJpTQ35bJ+YzLerxewX9PbHO2aQ7rxU+fKYWXRkZgJpzN88qFkpOx5hxddFg7CZw9h3/wdtq2JZ9PO4urKoyjMSufKI0dx7IRcXl39Mo8veZzFNYtJ86Zx0vCTOGfMOQzJGrLHu59fMZ+fvP0TmoPN3HzwzUwtOpDf3XUf/6/2F9QNP5G882cl9sDs1sI93yJYv54J1Tfx7QNG8ZsTxvHMsmf47Qe/5fcH/J5TRp7S7bu95dNbeGrRg/xzYwG/ybie5380o9v30V8ZY+Zaa/tIc3HPKxmcYi/55RAAQhaC4QhuYzAuQyh22f822drEapPNb9hyme2XT9su883r2DJb/ZxsudyW29n6efM2tl6v43a3XWfL7W67D9Pp680i2+zTbFGI2A7r2ViIWxcz7fO2Wa/DfjvE3/FaoC2nb/XafPN3R2f7/yZmu83xaV/TdHjdvs/O/q6xHQ+/9BKDISMpg+zkbLKTs8lKzmp/zkjKwB/20xJsob6tnnVN6yhvLKe2rRYAFy6Kk0oZ2RjhhOqlHNZWTYPJ5uvSEyk57H8oHT5+j+Nrq1xB4z9m4g3WcvXY03m/7R2mFE3hT4f+iWxfdqfrBMIBzvj3GfgbN/DUiuX8dfgD/O/5x+9xLDuzZEMDGX+fxM2DipmX5ubtM9/GvZOT+NZa9pk1gX3r8hk++m8JcWl6d+XIXm2Ry0/J56LxF/XmLqUHuYyr/ZJKunjpdDASpKKlgvVN61nftJ51TetY07iGVfWreGHFC+0dgAD4Sn18lT6QZziFzPoWRtSs4tDqFxn1+pM0vJrHK75DqBx8PKVjp7LP4FzKclK6f7yTUBt88QzBt/4Pb/0qFrrG8ZPWi2j0DOenxw5k2KC1vLnmr9z61Nu0hloZlTOKX0/7NTOHzdxhV7m7au/CvXl05qNc8dYVXDH7Cq7c90p+fvEF3H/7PC7++lEq3rufwoO+023763VfPg/r5vBQ7lV4klL44WHDqfHX8Me5f2Sfwn04ecTJPbLbwRmDaTWAq54KXVopDmryePkgr7T9fX1riJbYeEgel4u0JA/JXhfWQjBsCYQjWBu9Xzx6z3h0XiR2ewA22heS27hwbW5uIzp9i2IgViXY9hk7vwWB9s3Fppitb2TYavnN24rtwrD52XTY1jfvN69pOq7bYatb/Nd0OnXLZTucLHa1nzQmNs0Ve6Z9mc1zXR1PMm91wtllouu6tt6++Wabrm2mG9zt+/tm2ubl3K4O6+3sRHf7sfvmdoetP2tn67WfBO+w3pbxb2d+h5Pom/e1+XaLLY9LJ9vEtU2suzR/633Gluvs83c2v+PJ+VAk1P7oeAJ+62mhSKj9JPvm6dZa0pLSyPBmkJ6UTro3vb1gy0zK3GmxsbWq1ioWVS1iYdVC5m6ay3vB+bxdkkYy2Yz1uzi+7hkGPvIwC13jqR97DnsdcS45Wbt+9U3FZy+R/MIPaCTCVaNnsKLtHU4beRrX7n/tFre/bC3JncSvpv2K7776Xe7JzWb80jtZUXkowwt6sNd5YNH8DznB1DA3KZ+DSw/t0nE1xpBjDWFPW79rkev9ceSkX/O6vJSml1KaXrrNPGstVa1VrGpYxcr6laxqWMWq+lWsaljFF2Yd7+VGuJ/oINDZERdDA+8xrOYtwm95+LKtlMrIePzZh1JYNoaxJZmMKspgYG4qeWlJu1Tg2UiE6mUf0zjvaQpXPElaqI7lkUHcEvoJlYPGsM/wChp5hn+Vz6F1VSs5yTkcP+x4Thh+AnsX7N1jg2cWpRVx/zH38+v3f82f5v6J5bXL+c7F1zP37wsY9+bV1BRPIHdEAjaABFvhtV/TmDWa69fvyy+PG0leWhI/efsaWoIt/Gb6b3rsmA7KjHZeVOdqoLo+0N7yLdLbRuWM4o0z3thi2sqqZloDYcYMyMDl0vdSpC/JT8nnkIGHcMjAQwBoDjbzyYZPeH/9+7y37j1u8LVxAzAsUM3ha28ifOf/4fcdRtLkc5g45WDyMnZ8Br2hci1rnv4VYzY8z/3pZdwzIINAeBW/nf5bTh91epdi3G/Afpw4/EQeWPFvnmz4lCf/8yrXXHTann70HbJfPs9HKSnUR/y71Et1Ni6C7iDNGn5AxBnGGApSCyhILWC/AfttMS8YDrK2cS0rG1ayumE1q+pXsbp2ObPrVlAbbgGqgf/itm8zoDLC5+uTIJRJIJhLMJKPxzsYX/ogMnxZZCVnkuvLxONKxoaDpPmr8DWvI6XhK3JaF1ESWozH1UC528v7qYP5KmN/mtKD1IX+TSDyFKvXw5DMIZw4/EQOH3g4+xXv1z7mW09L8aTwh4P/wPDs4dw1/y5WN67mh6f9ifqnzsH18FnU/OBtcot23LtT3Hn/dqhfw7W+GxhWmMlFM4bw4tcv8saaN7hq36sYnt1zg5EOzhwMQIVpJRgO0eAPkZXS++PliHRG4xqK9B9p3jQOG3QYhw06DGstX9d/zTvl7/BO+X+5b9NnRLIjZIY/ZcaX77BhjiE5vBeenEn4SieQkVtEUrKPcFM1besWklb+LgNbP+KLNB9XDh7JJrefibnD+O0Bv20feqerfrLvT3h7zVtcX5DPt1fcwbJNRzGyqPMO2/ZUdVMbE+vf5q7SMjK8Scwo7frtDjnGS507pF4rReKR1+1lWPawbW7IBahvq2d1/SpWrfuQVes+YW3DKja21lDhraQio4aQWQF8HF04FHs0b7URA2TFHqTGHgBNZCdvYFTOKMbkHsCY3DFMKZpCcbpzg3EbY7h00qWMyB7Bte9dy2+bf8Xlh/+eo2b/hJV3n4Hrx6+TndkzP7LdrnY1vPcnluV/ixfKh/HQxXtR7a/gpo9vYnLhZC4cd2GP7n5A2gCSjJu1HjfZNFHd1KZCTkREHGWMYXj2cIZnD+c7479DQ6CBD9d/yDurXue/a9/h5fRWYCXZ4eWMKH+M/NVhMiMRwsbQ4HKxNCeZNYXFWAMjcwZy1fiLOWbIMbt8+SdAXkoeV065it9/+HtOy1rCsy++wC8uPrf7PzTw0ccf8C33Oj5IGc5Rg48gyd31Dl+yXcmsd7WRrEJOJLFkJWcxsXASEwsnweQffDMj1EakZiW1VV9RUb2UuuaNNAUaaQo20WhD+I0L6/ZikzOwvkx8GaWkZRST7k0nIymD4rRiitOLu/Vet+505OAjGZgxkMtnX84Na2+nat9LuHjeXbx554Xsd+VjZKbE4bANHVkL//4xYePi4g0nMXNiMdOGZ3PJa5cQsiFumHHDbiWdXeEyLgYm57La20i+aaC6OcCw3h3SR0REZIcykzI5esjRHD3kaCI2wqKqRXxeuYClFQtYWrWYzwP1NIX9eIyHjORMhuaOYGbeXhw26DBG54ze41sGTht5Gs8tfZpbwwv55eo7mb92JnsPzO6eD9dB24KneTcllRYb5Nghx+7SuhluHy2uepJaVMiJ9A2eZFyFY8grHEOe07H0kDG5Y3h05qNc9dZV/LnyRTaMPJL/XfY6z/71Mo684k4yfXHcuvTZg/D1W/w95VJaUoq5/qTx3D7vduZumstNB93Ufv9aTxuUVsyahnXkm3qqGtt6ZZ8iIiK7w2VcTCiYwISCCTDuvF7b528O+B1n/fsMPs1fT/DpR5n04x906z3ltU1t7FP/OjeVlpCTnMbU4qm7tH6mJ4WmALj7WYtc1we1EpG4lJ+Sz71H38vJI07m8dASvj90Et/yP8Hjt19DdVOcFib15fDqtazN3Idba2dw06kTmV/9Hv9a9C/OGn0Wxw/r+S6ONyvNGMh6j4d8aqlqDvTafkVERBLF6NzRnDvmHJ7MzGB88z28MH9dt25/3nv/Ic1TwUfJQU4acRIe1661NWUkpRNwGZr8Ld0aV7xTISfSByS5k/j9Ab/nmqnXMNfVyOmDhjE99AD33HEjG+pbnQ5vS+EgPPkdQuEwF1Sex6n7DGLggFp++e4v2StvL36x3y96NZzirKG0uFxkuavjt/AVERFx2A/3uYJCTzpPFjTxzn9m0Rrovh4iXfMf4vH0HCLYLveq2VGGN9o3QDhc120xJQIVciJ9hDGGc8eeywPHPIArPZ/zSwZQ4r2P+/76e5ZuanQ6vG+8+Tso/4RrQ98juWgklx2Zx4/e+BEZSRncfvjtu3Rzc3coiXWgk5FSQ3WTWuREREQ6k+ZN49oDfseS5CTy0h/kn28u7JbtLlu1himt7/JMTg77F+/f3qP0rkj3ZQFgQlVYa3eydN+hQk6kj5lQMIEnTniSA8sO4pb8HDbmPsH9//glsxdvcjo0WPQcfPBX/pN0LC9zADefOZSr/nsZLaEW7vrWXRSmFvZ6SMUZJQC4khuoUouciIjIdh0+9ChOGjCdx7M91M29jsUbG/Z4m6te/jPvpHuoNAHOGXPObm0jMzkHgFR3HS3d2FIY71TIifRBWclZ3H7Enfxin6v4IC2Nd8tm89RzP+Sed1Y4d6ZqzcfYZ77P0qSx/KzpbG44fQjXfXoF65rWcfvht+/y2DbdpSQtWsiFvI1qkRMREdmJaw77I0UuHx8NWMJfH3ucUDiy29uqrK1n741PcmfuAEZkj+CwgYft1nYyUqLd2vlcjSrkRCTxGWM4f8J3efL4Jyn1ZPBu8VJe//xcLn34NRr8wd4NpuIr7GPnUGHyOavhx/zsxIHcu/wXlDeWc+cRd24zAHxvyk7OJgVDs7uVqma1yImIiOxIelI6/++QP1DucWM9f+au1z7f7W199vStzE8LstYT4uIJF+Myu1eapKfkA5DkaezWe/finQo5kT5uWN5oHvr2u1yWOYFl6U3MDf6M4+/+Xz5dVdk7AWxahP3XTOraLGc2/5Rzjs7jwdU/pbKlkru+ddcudzHc3YwxlLhSqHUHNPyAiIhIF0wZfBiXDzmRdzNcfP3llXywvGqXt7F+43omld/LzfkDGJY1jGOGHLPb8WSkR2/N8LqaaQ2qkGtnjBlojHnLGPOVMWaRMebHsem5xpjXjTHLYs85PR+uiOwOjyeJ/znlEZ4d9yP2aWujNusVfvjaKVz69EM09WTr3Kr3ifxrJtV+y2mt17LvIa08Vn4NPo+Ph457yNGWuI6KvRlUuMI0+IMEQrt/iYj0L8qPItKfXXLIjRyRXMLbeVU88MwvWFuza13/L3nkap7N8rDJHebqqVfv8pADHWWmDQDA425RIbeVEPBTa+1YYBrwI2PMOOAa4E1r7Ujgzdh7EYljg6deyt9PfpY/tiaR76rmvaabOfiBk/nL+y8R3oNr3LdhLcy9n8isk1jtT+EscyWF+3/E65v+wfSS6Twy8xGGxXqLjAclvjw2eNxk0UyNxpKTrlN+FJF+yxjDjSc/zkibzJyCT/jtv37X5U7D3nnlSQr9L/P3nGyOHnI0B5QcsEexpKQW4LYW42rVpZUdWWs3WGvnxV43Al8BpcBJwKzYYrOAk3soRhHpRmbAXhz53Xd5ofg4flNVQ7Z7Ffcsv5qp9x/Pze8+RFtoDy8vbK7C//C34d8/5rXIaL474Fjqhj7I102fc/V+V3PH4XeQ68vtng/TTQakFlHvdpPrqlbPldJlyo8i0t+l+bL5xynPUGxdLMp9kavvuXanLXOLvlhAySc/4arCAeSl5vPrab/e4zhMcgbpkQjW3UZrMLTH20sUu3SPnDFmCDAZ+BgostZugGgyA3q/33AR2T1JqXhn3sYZZ73AS8FsflVVQ2FkLQ99fTP7PXgI5z57NW+t/pBwZBfOagX91Lx+G61/3JvVq2ZzVvYB/O/wIFUpr3PIwIN54eQXOG/ceRhjeu5z7abC9GIA0j2bVMjJblF+FJH+Kjd7MPce/wilETfzsl/hV/d/h3eWdD7k0Zy5nxJ69gx+MSCd2mQffzrsz2QlZ+15EG4PGRFLxBWgNdB/bpHo8sWoxph04GngSmttQ1f/GDPGfB/4PsCgQYN2J0YR6Sll++L7/juc9cXTnPH2TXxcvY6HM/N4j5e54u2XcNs0BqeNZ78BUzhg4ATGFQynMLWwvVepxtYAq5d8yup592Or32RJUpg3ivIpTw7jMhs4YtDhfG/C9xibN9bhD7pjhRkDAUjxVmkIAtllyo8i0t8VFo7ngTNe5qdPn8IHeV/SPPtoXnj7VA7e/zsMzc+ktqaKyvfuxdf4EH8syaAhycefDv8zEwsmdlsMGRja3MF+dY9clwo5Y4yXaJJ62Fr7TGzyJmNMsbV2gzGmGKjobF1r7d3A3QBTpkzpP0OtiyQKlxsmnolrr1OZvuQlps+bRfOK2byb6mN2SoBPQp/yeMvHPP51bHkLSdZNkrWETJg2A9ZtoDANN272ypvEOUOP4Lhhx5Ef6w443hXlRO/XS/LWUq0hCGQXKD+KiESlZ5by9ws+5LFXf8ydG99iiftJvvr0UYYGIJkAi5KTWDcgi6HpA7nzsD8yJndM9+4fN02uIK2B/nNp5U4LORM9tXgv8JW19o8dZr0AXAj8X+z5+R6JUER6h9sD406EcSeS1riRY1bM5pgVbxHasJDqTV+z0g2rvF6q3W4aXS4aPOkYbw6ezKGMGn4o4wfvy6icUaR4Upz+JLusMGcEAB5vA1VqkZMuUn4UEdmScbs557g7OKm5kmc/vIXZGz5hqbcJ48piWM5wfjj2TI4ddhxel7fb951u3ARdEbXIbWUGcD6w0BgzPzbtf4kmqCeMMRcDa4AzeiRCEel9GQNg72/D3t/GAxRFIhQFm5nW1ggeHyRnRgu/PiI9rZDUSARXknqtlF2i/Cgi0onUtALO/dYfOLcX95lmvARdbbpHriNr7XvA9i74P6J7wxGRuORyQXJG9NFHFVoXAW8rtSrkpIuUH0VE4keaO4kAflrUa6WISP9ShIdmd4DaFhVyIiIiiSbVnYzfgF/jyImI9C+F7hTqXSFqW4JOhyIiIiK7KNXtI+iC5kD/6bSs79zkIiKyBwqTMqixDdgWv9OhiIiIyC5K86RAABoDOx6QvC9Ri5yICFCYnEPYGBoDdYTC/edGaRERkb4g1ZsKQEuo2eFIeo8KORERoCilAACvp4b6Vl1eKSIikkjSvOkA+AONDkfSe1TIiYgAhWnFAKR7Nuk+ORERkQSTmhTtWTsUrHM2kF6kQk5EBCjMHAhAqrdKPVeKiIgkmNTkTADCoXqHI+k9KuRERIDczMEYa0n21GosORERkQSTlpwFgI00OBxJ71EhJyICeDOKyI5EcHsa1CInIiKSYNJ82QDYiO6RExHpX1LzyAuHsZ5m3SMnIiKSYFI3F3JWvVaKiPQvKTnkhSNEPK26tFJERCTBpKbkxV6pkBMR6V9cbvJx0+YO6NJKERGRBJPiy8VYC/gJ9pPxYFXIiYjE5Lt8NLtDVKtFTkREJKGY5HRSrcXlbqM1GHY6nF6hQk5EJCbPk0bQZalt6T83SouIiPQJSWmkRSLgasMfUCEnItKv5Me6Lq7yVzkciYiIiOySpHRSIxZMQC1yIiL9TZ4vB4CGQI3DkYiIiMgucXtJtWDdQRVyIiL9TV5qIQDN4VrCEetwNCIiIrIrUjGETYgWXVopItK/5KcVA+Bz19DQqrHkREREEkkKLkKusO6RExHpb7IzSnBZi89ToyEIREREEkyq8RByhXVppYhIf+NOKyA3HMHrqVchJyIikmBSjIeAK6JLK0VE+p3UPPLCYVyeRmqbdWmliIhIIkl1JdHmsmqRExHpd1JzyQ+Hse5matQiJyIiklDS3cm0uaA10D9OxqqQExHZLNYiF/L4qVMhJyIiklDSPCkA1Le1OBxJ71AhJyKyWXImeRGL3xOgukmFnIiISCJJ80YLuca2Zocj6R0q5ERENjOGfHcKEWOpbK5zOhoRERHZBaneVACag2qRExHpd/I8aQBUtFQ6HImIiIjsilRvOgBNgSaHI+kdKuRERDrIT84CoKat2uFIREREZFekeqMnY1sDDQ5H0jtUyImIdJDvywWgIVDjcCQiIiKyK1KTMwEIBeucDaSXqJATEekgL7UQgKZQrcORiIiIyK5ITcoAIBxsdDiS3qFCTkSkg8y0AXispS1SRyRinQ5HREREuig1dntEMKRLK0VE+h1XWj654TBedy2N/pDT4YiIiEgXpfpyAAhH1NmJiEj/k5pHXjiC19NAjQYFFxERSRipKZsLOQ0/ICLS/6TmkhcOYzyN1KqQExERSRhJyVm4rSUS0YDgIiL9T2oeeeEw1t1CbbMKORERkURhktJJjVjC+J0OpVeokBMR6ShWyAU9fmpUyImIiCSOpDRSbASrQk5EpB+K3SNnjWVjo8aSExERSRhJqaRGLBHT5nQkvUKFnIhIR95UcmM/jRubKx0ORkRERLrMm0aKtURMkEAo4nQ0PU6FnIhIR8aQ74kOKFrRUuVwMCIiItJlniRSrSXiCtEaDDsdTY9TIScispX82ICiNf5qhyMRERGRXeGzLsKuEH4VciIi/U+eLw+A+kCtw5GIiIjIrkjBRcgVpjWgQk5EpN/JSi3AZS1NQRVyIiIiicRnPIRMmBYVciIi/Y8rLZ/csKU1Uud0KCIiIrILUoyXgMvqHjkRkX4pNY/8cIgQ9VhrnY5GREREuijFFSvkAiGnQ+lxKuRERLaWmkteOIzbXU+Dv+8nAhERkb4ixZVM2EBTW98fFFyFnIjI1lLzyAuHMZ4mapsDTkcjIiIiXZTqSQGgrq3Z4Uh6ngo5EZGtpeaRF45g3S3UNLc5HY2IiIh0UVqskGvwNzkcSc/baSFnjLnPGFNhjPmiw7RcY8zrxphlseecng1TRKQXxVrkrCvC+kb1XCnbpxwpIhJfUpPSAGgMqEUO4H7gmK2mXQO8aa0dCbwZey8i0jfECjmA8oYKh4OROHc/ypEiInEjIykdgMZgi8OR9LydFnLW2neAmq0mnwTMir2eBZzcvWGJiDgo1tkJwMamSoeDkXimHCkiEl/SkzMAaPHXOxxJz9vde+SKrLUbAGLPhdtb0BjzfWPMHGPMnMpK/UEkIgnAk0yuywdARUuVw8FIAupSjlR+FBHpfmmxQq4tUOdsIL2gxzs7sdbeba2dYq2dUlBQ0NO7ExHpFvnJWQDU+LdubBHpHsqPIiLdLzWWvwMBtchtzyZjTDFA7Fk3kYhIn5KdkoexUN+mQk52mXKkiIhDUn3RQi4YanA4kp63u4XcC8CFsdcXAs93TzgiIvHBnZpPdgSaQuq1UnaZcqSIiENSk7MBCIbUayXGmEeBD4HRxphyY8zFwP8BRxpjlgFHxt6LiPQdqXnkhiO0RuqcjkTimHKkiEh88fmiI76Ewn1/HDnPzhaw1p6znVlHdHMsIiLxIzWPgooQX7v6/qUZsvuUI0VE4os7OZ2USISQ1fADIiL9U2ouBaEgEVc94Yh1OhoRERHpiqR0UqwlFGl1OpIep0JORKQzsUHBjbuJ6qY2p6MRERGRrvCmkhqJEMbvdCQ9ToWciEhnYoUcrjBr69VzpYiISEJISiM1YgkRcDqSHqdCTkSkM5sLOWBV7SaHgxEREZEuSUon1aqQExHpv1LzyAtHAFjboGHAREREEoIniZSIJWyCTkfS41TIiYh0pkOL3IZGFXIiIiKJwoeLkAlhbd/urEyFnIhIZ1Jy2lvkKlqqHQ5GREREusqHm6ArTFso4nQoPUqFnIhIZ9wecnw5GAs1/iqnoxEREZEuihZyEZrbQk6H0qNUyImIbIc7vYj0iIuGYK3ToYiIiEgX+YyXoInQEgg7HUqPUiEnIrI9GUXkRizNoTqnIxEREZEuSnElEXRBfWvfHktOhZyIyPakF1EQDtNm652ORERERLooxZ0MQK2/yeFIepYKORGR7UkvpCjkJ2wa+nzPVyIiIn1FitsHQE2LCjkRkf4pvYiiUBDjqaeute8PLCoiItIXpHlSAahTi5yISD+VXkRRKAwmwtc1m5yORkRERLogLSkNgPq2Zocj6Vkq5EREtie9kKLYoOBf16x3OBgRERHpivSkdAAa1CInItJPpRcxIBQt5FbXq5ATERFJBBm+TACa2hocjqRnqZATEdme9EKKwtHBRNc1bXA4GBEREemKzYVca1uds4H0MBVyIiLb48smBy8uCxUtFU5HIyIiIl2QlpwFQFtb3x4+SIWciMj2GIM7vYjMsJtqvwo5ERGRRJDiywagLahLK0VE+q/0QgrChoZgldORiIiISBekxgq5YEidnYiI9F8ZAygJh2iN1DgdiYiIiHSBNzkTr7UEwirkRET6r6wyBgZbCJlarLVORyMiIiI7k5ROaiRCyLY4HUmPUiEnIrIjWWWUhNrAFaSiudbpaERERGRnktLIiEQIqpATEenHMkspDEWHIFhcudbhYERERGSnfFlkRCwhWp2OpEepkBMR2ZGsgRSFo4OCL60udzgYERER2amUbDIjEYL4nY6kR6mQExHZkawyymItcl/XrXE4GBEREdkpbwoZEQi42pyOpEepkBMR2ZH0InKtG0/ExdpGXVopIiKSCFJxE3CF8AfDTofSY1TIiYjsiMuFK7OY/JCbTS3rnI5GREREuiCdJNpcYepbg06H0mNUyImI7EzWQMpClrrgRqcjERERkS5Ic6cQcsGmxkanQ+kxKuRERHYmq4xhYT9+KonYiNPRiIiIyE5keFIBWNfQd4cOUiEnIrIzOUMZ0dYIJkhFc6XT0YiIiMhOZHgzAdjUWONwJD1HhZyIyM7kjWBQKHqN/cJNXzscjIiIiOxMTkoWABUtdc4G0oNUyImI7EzecAYGo0MQfFGpQk5ERCTe5aflAVDVXO1wJD1HhZyIyM7kDac4FMJlDYurVzgdjYiIiOxEQUYxAA0t6x2OpOeokBMR2RlfFp60AvIDXlY1LHc6GhEREdmJvMxBADS3qpATEenXTN4IRoYiVLWtdjoUERER2YnkzFIywhH8wb7bSZkKORGRrigYzaRAAwFTQ2Og745JIyIi0idkDCAvHKYtrF4rRUT6t+JJjGtrAmBhxRKHgxEREZEdSi+KFnL03ZOvKuRERLqieG/GBKJDELyzep7DwYiIiMgOJWeQEwG/q9XpSHqMCjkRka4oHEehNaQHvXyy4TOnoxEREZEdMYYskmlxBfAHw05H0yNUyImIdIXXhykYy7gArG7+0uloREREZCdy3Wn43ZYlFX2zwxMVciIiXTV4Ogf7awhQw7rGdU5HIyIiIjtQ5isAYN76vjkGrAo5EZGuGnYoB7VGb5p+dskbDgcjIiIiOzIsfxQASzYtdTiSnqFCTkSkq4YcxNAwZAeSeXHFa05HIyIiIjswvHQfADbVLnQ4kp6hQk5EpKt8mZiRRzGztZF1rV+wqXmT0xGJiIjIdmQUTaAkGKLarxa5bRhjjjHGLDHGLDfGXNNdQYmIxK29v815DZWA5Q8f/tPpaCSOKUeKiDisYAxjAyFqzXoiEet0NN1utws5Y4wbuBM4FhgHnGOMGdddgYmIxKXRx1GSO5rDm0O8Vv40X1R+5XREEoeUI0VE4oAnmbGeAdR5gzy9cK7T0XS7PWmRmwost9Z+ba0NAI8BJ3VPWCIiccrlxnX8n/hVTRVZ4RDfefFC/v7hA6ytryASiTgdncQP5UgRkThwwvizcFvLo59cR0Nri9PhdCvPHqxbCqzt8L4c2H/PwhERSQCDplFw5izufu6H/G9ugDuX/oE7l/4BAI+1eCwYEvMSDuN0AH2HcqSISBwo2f9/OGfBfTyUupoDH59KkgV3gubore1JIddZvt/mqBhjvg98P/a2zRjzxR7ss7/JB6qcDiKB6HjtGh2vXaPjtWtGOx2Aw3aaI/tYfkz0/z8SOf5Ejh0SO/5Ejh0SO/5Ejh26KUfuSSFXDgzs8L4MWL/1Qtbau4G7AYwxc6y1U/Zgn/2Kjteu0fHaNTpeu0bHa9cYY+Y4HYPDdpoj+1J+VPzOSeTYIbHjT+TYIbHjT+TYofty5J7cI/cpMNIYM9QYkwScDbzQHUGJiIgkOOVIERHpUbvdImetDRljLgNeBdzAfdbaRd0WmYiISIJSjhQRkZ62J5dWYq19CXhpF1a5e0/21w/peO0aHa9do+O1a3S8dk2/P167mCMT/XgpfuckcuyQ2PEncuyQ2PEncuzQTfEba/tGry0iIiIiIiL9xZ7cIyciIiIiIiIO6JVCzhhzjDFmiTFmuTHmmt7YZyIzxgw0xrxljPnKGLPIGPNjp2OKd8YYtzHmM2PMi07HkgiMMdnGmKeMMYtj37PpTscUz4wxV8X+X/zCGPOoMcbndEzxxBhznzGmomP3+caYXGPM68aYZbHnHCdjjBc7y4cm6vbY/M+NMfs4EWdnupKbjDGHGmPqjTHzY4/fOBFrZ4wxq4wxC2NxbdNjXJwf+9Edjul8Y0yDMebKrZaJq2O/J78LTv/duJ3Y/xDLmZ8bY541xmRvZ90dfs96w3biv84Ys67D9+O47awbj8f+8Q5xrzLGzN/OuvFw7Dv9neyx7761tkcfRG/yXgEMA5KABcC4nt5vIj+AYmCf2OsMYKmO2U6P2U+AR4AXnY4lER7ALOCS2OskINvpmOL1QXRg55VASuz9E8BFTscVTw/gYGAf4IsO024Brom9vga42ek4nX50JR8CxwEvEx2HbhrwsdNxd4htp7kJODRef4eBVUD+DubH7bHv5Hu0ERgcz8d+d38X4uHvxu3EfhTgib2+eXu/aTv7njkY/3XAz7rw3Yq7Y7/V/NuA38Txse/0d7Knvvu90SI3FVhurf3aWhsAHgNO6oX9Jixr7QZr7bzY60bgK6J/TEonjDFlwEzgHqdjSQTGmEyiP5T3AlhrA9baOkeDin8eIMUY4wFS6WTMzP7MWvsOULPV5JOInjAg9nxyb8YUp7qSD08CHrBRHwHZxpji3g60M/0gN8Xtsd/KEcAKa+1qpwPZkT34XXD878bOYrfWvmatDcXefkR0bMi4tJ1j3xVxeew3M8YY4Ezg0d6MaVfs4HeyR777vVHIlQJrO7wvp2/98PcoY8wQYDLwscOhxLM/A78AIg7HkSiGAZXAv0z0ctR7jDFpTgcVr6y164BbgTXABqDeWvuas1ElhCJr7QaIJjag0OF44kFX8mFC5Myd5KbpxpgFxpiXjTF79W5kO2SB14wxc40x3+9kfkIce6JjEm7vD9l4PfabdeV3IRH+Hb5LtPW2Mzv7njnpstilofdt59K+eD/2BwGbrLXLtjM/ro79Vr+TPfLd741CznQyTV1ldoExJh14GrjSWtvgdDzxyBhzPFBhrZ3rdCwJxEP0soW/WWsnA81Em/mlE7FkdxIwFCgB0owx5zkblSSoruTDuM+ZO8lN84he8jcJ+CvwXC+HtyMzrLX7AMcCPzLGHLzV/EQ49knAicCTncyO52O/K+L638EYcy0QAh7eziI7+5455W/AcGBvoiclb+tkmbg+9sA57Lg1Lm6O/W7+Db/Lx783CrlyYGCH92XosqSdMsZ4iX4BHrbWPuN0PHFsBnCiMWYV0Sbow40xDzkbUtwrB8qttZvPpD9FtLCTzn0LWGmtrbTWBoFngAMcjikRbNp8WVrsucLheOJBV/JhXOfMneUma22DtbYp9volwGuMye/lMDtlrV0fe64AniV6GVNHcX3sY44F5llrN209I56PfQdd+V2I238HY8yFwPHAuTZ2U9PWuvA9c4S1dpO1NmytjQD/pPO44vnYe4BTgce3t0y8HPvt/E72yHe/Nwq5T4GRxpihsTNJZwMv9MJ+E1bsGuB7ga+stX90Op54Zq39pbW2zFo7hOh3a7a1Vq0lO2Ct3QisNcaMjk06AvjSwZDi3RpgmjEmNfb/5hFEr3mXHXsBuDD2+kLgeQdjiRddyYcvABeYqGlEL+Xd0NuBdqYruckYMyC2HMaYqUT/zqjuvSg7Z4xJM8ZkbH5NtOOKL7ZaLG6PfQfbbZGI12O/la78LsTl343GmGOAq4ETrbUt21mmK98zR2x1v+cpdB5XXB77mG8Bi6215Z3NjJdjv4PfyZ757ndXLy07ehDtCWop0Z5Yru2NfSbyAziQaFPq58D82OM4p+OK9wdx1mNXPD+IXloxJ/Ydew7IcTqmeH4AvwMWE00KDwLJTscUTw+if1huAIJEzyheDOQBbwLLYs+5TscZD4/O8iHwA+AHsdcGuDM2fyEwxemYO8TeaW7aKv7LgEVEe1v7CDjA6bhjcQ2LxbQgFl9CHftYfKlEC7OsDtPi9tjvyu8C0cvWX+qwrqN/N24n9uVE71/a/N3/+9axb+97FifxPxj7Xn9OtDgoTpRjH5t+/+bveodl4/HYb+93ske++ya2koiIiIiIiCSIXhkQXERERERERLqPCjkREREREZEEo0JOREREREQkwaiQExERERERSTAq5ERERERERBKMCjkREREREZEEo0JOREREREQkwaiQExERERERSTD/H6ftAxxwPGHdAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "f, ax = plt.subplots(1,2,figsize=(15,5))\n", "\n", "ax[0].plot(cyl10.profiles[-1].r/ANGSTROM, (cyl10.profiles[-1].density/(KILO*MOL/METER**3)).T, label='cyl10')\n", "ax[1].plot(cyl20.profiles[-1].r/ANGSTROM, (cyl20.profiles[-1].density/(KILO*MOL/METER**3)).T, label='cyl20')\n", "ax[0].plot(sph10.profiles[-1].r/ANGSTROM, (sph10.profiles[-1].density/(KILO*MOL/METER**3)).T, label='sph10')\n", "ax[1].plot(sph20.profiles[-1].r/ANGSTROM, (sph20.profiles[-1].density/(KILO*MOL/METER**3)).T, label='sph20')\n", "ax[0].plot(car10.profiles[-1].r/ANGSTROM, (car10.profiles[-1].density/(KILO*MOL/METER**3)).T, label='car10')\n", "ax[1].plot(car20.profiles[-1].r/ANGSTROM, (car20.profiles[-1].density/(KILO*MOL/METER**3)).T, label='car20')\n", "\n", "ax[0].set_xlim(0,10)\n", "ax[1].set_xlim(0,20)\n", "for a in ax:\n", " a.set_ylim(0,50)\n", " a.legend()" ] } ], "metadata": { "kernelspec": { "display_name": "feos", "language": "python", "name": "feos" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.10" } }, "nbformat": 4, "nbformat_minor": 4 }