{ "cells": [ { "cell_type": "markdown", "source": [ "Potential energy surface for twisting of ethene\n", "===============================================" ], "metadata": { "collapsed": false, "pycharm": { "name": "#%% md\n" } } }, { "cell_type": "code", "execution_count": 1, "outputs": [], "source": [ "from pymolpro import Project, all_completed, no_errors\n", "import numpy" ], "metadata": { "collapsed": false, "pycharm": { "name": "#%%\n" } } }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "root = 'dihedral_scan'\n", "projects = []\n", "import os\n", "\n", "if not os.path.isdir(root): os.mkdir(root)\n", "grid = numpy.linspace(0, 90, 18, endpoint=False)\n", "for i, D in enumerate(grid):\n", " p = Project(f'd{i}', location=root)\n", " projects.append(p)\n", " p.write_input(\n", " f\"\"\"geometry={{\n", " H1\n", " C1 H1 rch\n", " C2 C1 rcc H1 A\n", " H2 C2 rch C1 A H1 D\n", " H3 C2 rch C1 A H2 180\n", " H4 C1 rch C2 A H1 180}}\n", " rch = 2.05; rcc = 2.65; A = 121.5; D ={D}+.01;\n", " basis = cc-pVTZ\n", " rhf\n", " set, mcocc=[3,2,2,2], mcclosed=[3,1,1,2], mcsym=1\n", " casscf,so_sci\n", " {{optg; inactive, D}}\n", " set,mcsym=[1,4]\n", " casscf,so_sci\n", " {{mrci; wf, sym=1;}}\n", " {{mrci; wf, sym=4;}}\n", " put,xml\"\"\")" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "from multiprocessing.dummy import Pool\n", "from operator import methodcaller\n", "\n", "with Pool(processes=8) as pool:\n", " pool.map(methodcaller('run', backend='local', wait=True),\n", " projects, 1)\n", "assert all_completed(projects) and no_errors(projects)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [ { "data": { "text/plain": "
", "image/png": "iVBORw0KGgoAAAANSUhEUgAAApsAAAGwCAYAAADmEa4MAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAA9hAAAPYQGoP6dpAACcy0lEQVR4nOzdd1xV9RvA8c9lg0xZoiAoKoh7K5ppJpiWpqnZsOGsRFPLSs1cmVpWtkwz1MwwzfETK8E9c+AgF+JWpoAiICDz/P44evWGi+VlPO/X6/uSe+ZziDgP36lRFEVBCCGEEEKIUmCg7wCEEEIIIUTFJcmmEEIIIYQoNZJsCiGEEEKIUiPJphBCCCGEKDWSbAohhBBCiFIjyaYQQgghhCg1kmwKIYQQQohSY6TvACqr3Nxcjhw5grOzMwYGkvMLIYQQ5UF+fj5XrlyhWbNmGBlJGvUo5LukJ0eOHKF169b6DkMIIYQQRXDgwAFatWql7zDKBUk29cTZ2RlQf1hdXFz0HI0QQgghHkVcXBytW7fWvsfFw0myqSe3m85dXFxwdXXVczRCCCGEKAzpAvfo5DslhBBCCCFKjSSbQgghhBCi1EiyKYQQQgghSo0km0IIIYQQotRIsimEEEIIIUqNJJtCCCGEEKLUSLIphBBCCCFKjSSbQgghhBCi1EiyKYQQQgghSo0km0IIIYQQotRIsimEEEIIIUqNJJtCCCGEEKLUSLJZ0SRfgtRYyM/XdyRCCCGEXly/fp0dO3boOwxxi5G+AxAlbONEiFgPhqZgWxPsPO4q7uq/tu5gZq3fOIUQQohSkJGRwbPPPsuBAwdYsWIFvXv31ndIlZ4kmxVNXg5oDCEvC66eUcu9WNjfSTz/m5Bau4Kh/GgIIYQoX7Kzs+nbty979uzB1tYWT09PfYckkGSz4nl5BeTlQmo0JF+8VS7d9fVFyLwGGVfVEnOo4DU0hmDrdicB/W9Cam4HGs3jeiIhhBDiofLy8njttdfYsGED5ubm/PXXXzRu3FjfYQkk2ayYDI3uJIb3cjMVrl+6dzJ6/RLkZd/5fC+m1nea5Ks3g7r+4NxAElAhhBB6oSgKAQEBrFixAmNjY9auXYuvr6++wxK3SLJZGZlZQ7VGavmv/Hy4Ea9bE3p3QnojHrJSIf6YWiLWw5ZpatN7PT818azVEUwsHusjCSGEqLw+/vhj5s+fj0ajYdmyZfj7++s7JHEXSTaFLgMDsK6uFvd7/FWYnQHXL6s1oFfPwYUdcH6H2mx/cJFajMzA4wmo568W25qP/zmEEEJUCnPmzOGzzz4DYP78+fTv31/PEYn/0iiKoug7iMooOjoaNzc3oqKicHV11Xc4xZOTCRd2wZlQOL0RUi7r7nesf6fW062NDD4SQghRIgIDAxkyZAgAs2bN4sMPPyz1e1ao9/djIm99UXzG5moyWc8PuiuQeApOh6iJZ9R+SIxQy55vwMwG6jytJp51noYq9vqOXgghRDm0evVqhg0bBsAHH3zwWBJNUTSSbIqSpdGAU321dBgDmclwdguc2QhnNqkj4Y+vVgsacG11K1HtBs4NZZCREEKIh9q0aRMvv/wy+fn5DBkyhFmzZuk7JPEA0oyuJ5WyGj4/D6IP3mluv3JMd791DajbVa31rP0kmFTRT5xCCCHKrH379vH000+Tnp5Ov379WL58OYaGho/t/pXy/V1MkmzqifywAikxao3n6VB1oFFOxp19hqZQ6wk18aznd/9pnIQQQlQax48fp2PHjiQnJ+Pn50dwcDCmpqaPNQZ5fxeeJJt6Ij+s/5FzEy7uvlXrGaqOdr+bozf4PA8N+4Cjl15CFEIIoT/nz5+nQ4cOxMXF0a5dOzZt2kSVKo+/BUze34UnyaaeyA/rAygKJEbeaW6/vBeUvDv7nRuqSWeDPlC1lv7iFEII8VjExcXRoUMHzp8/T6NGjdixYwd2dnZ6iUXe34UnA4RE2aPRgJO3Wtq/C5nX1drOE2vUwUZXjqtlyzSo3hwavgANeoNNDX1HLoQQooRdu3YNPz8/zp8/j6enJ6GhoXpLNEXRSLIpyj5zW2jyoloyrsGpP9XR7Bd2QuxhtWycCDXbqYmnTy+wdNJ31EIIIYopPT2dHj16cPz4cVxcXNi0aRMuLi76DksUkiSbonyxqArNX1PLjQQ4uQ5OrIVL/6jN7Zf3woYP1BWMGr4A9Z9TzxFCCFGuZGVl0bt3b/bt20fVqlXZtGkTtWpJ16nySJJNUX5ZOkHroWpJiYGT/4PjayDmoDq6/cIO+GsseD6lJp5e3dV14YUQQpRpeXl5vPrqq9pBQH///TcNGjTQd1iiiCTZFBWDTQ1oN0It1y6otZ0n1kD8sVsTym9Up1Oq21VNPOt1AxMLfUcthBDiPxRF4a233mLVqlWYmJjwv//9jzZt2ug7LFEMBvoOoCi2b9+ORqO5ZwkLC9MeFxYWRpcuXbC1tcXOzg4/Pz/Cw8MfeO2ffvqJTp06YW1tjUaj4fr16wWOSU5OZuDAgdjY2GBjY8PAgQPveZzQk6q14Imx8NZuGBEGncaDQz3Iy1L7e656E76oA6sGwam/IDdL3xELIYS45aOPPuLnn3/GwMCA5cuX8/TTT+s7JFFM5TLZ9PX1JS4uTqcMGTIEDw8PWrZsCUBaWhr+/v7UrFmT/fv3s3v3bqytrfH39ycnJ+e+187IyKBbt25MmDDhvse8/PLLhIeHExISQkhICOHh4QwcOLDEn1OUAMd60OkjGHEA3toDHcaqE8TnpKuDjH5/Gb6oC2vfhjObIe/+PxtCCCFK1+zZs/n8888BWLhwIX369NFzRKJEKBVAdna24uTkpEybNk27LSwsTAGUy5cva7cdPXpUAZSzZ88+9Jrbtm1TACU5OVln+8mTJxVA2bdvn3bb3r17FUA5derUI8ccFRWlAEpUVNQjnyNKSH6+okQfVJSQCYryZX1FmWx9p8zyUJTgdxXlwm5FycvTd6RCCFFpLFiwQAEUQPnyyy/1Hc59Ffb9/dlnnyktW7ZULC0tFUdHR6VXr14PzRdef/117ffi7uLj46M9ZvHixfc8JjMzs1jPVxrKZc3mfwUHB5OUlMQbb7yh3ebl5YWDgwOBgYFkZ2eTmZlJYGAgDRo0wN3dvcj32rt3LzY2Njr9R9q2bYuNjQ3//PPPfc/LysoiNTVVW9LS0oocgygmjQZqtAD/GTD6OLwZAq2GQhVHyLwGhxbDku4wtxFs+kTt9ylrHwghRKlZsWIFb731FgATJ05k7Nixeo6o5OzYsYMRI0awb98+Nm3aRG5uLn5+fqSnp9/3nG+++Uan9TYqKoqqVavSr18/neOsra0LtPSamZmV9iMVWoUYIBQYGIi/vz9ubm7abVZWVmzfvp1evXoxffp0AOrVq0doaChGRkV/7Pj4eJycCs7h6OTkRHx8/H3PmzlzJlOnTi3yfUUpMTAA93Zq6TYLLu6C46vgZDCkRsOeb9Ti6A2N+kGjvrJOuxBClKCQkBBeffVVFEXh7bff1r6zK4qQkBCdz4sXL8bJyYlDhw7RsWPHe55ze0zIbf/73/9ITk7mzTff1DlOo9FQrVq1kg+6hJWpms0pU6bcd+DP7XLw4EGdc6KjowkNDWXw4ME62zMzMxk0aBDt27dn37597NmzhwYNGtC9e3cyMzOLFadGoymwTVGUe26/bfz48aSkpGjLyZMnixWDKAWGRuDZGXr9AO+fgf6/Qv2e6ij2xFOwdTp80wQC/eDAQriRqO+IhRCiXNuzZw99+vQhNzeXl156ie+///6B79KyJC0tTafFMivr0QabpqSkAFC16qPPAR0YGMjTTz9doGX2xo0buLu74+rqyrPPPsuRI0ce/QEeozJVsxkQEMCAAQMeeIyHh4fO58WLF2Nvb0/Pnj11tgcFBXHx4kX27t2LgYGBdpudnR3r1q176H3up1q1aly5cqXA9sTERJydne97nqmpKaamptrPqampRbq/eEyMzcCnp1pupkDEejj2h7pqUdR+tWz4UE1OG/UH7+5gaqXvqIUQotz4999/6dGjB5mZmXTv3p1ffvlF+74uD3x8fHQ+T548mSlTpjzwHEVRGDt2LB06dKBhw4aPdJ+4uDg2bNhAUFCQznZvb2+WLFlCo0aNSE1N5ZtvvqF9+/b8+++/1K1bt1DPUtrKVLLp4OCAg4PDIx+vKAqLFy/mtddew9jYWGdfRkYGBgYGOn8h3f6cn59f5BjbtWtHSkoKBw4coHXr1gDs37+flJQUfH19i3xdUYaZ2UCzV9WSFq9OHH/sD3WZzLOb1WJkDl7PqE3tdZ4GIxN9Ry2EEGXWmTNn8Pf3JyUlhQ4dOvDHH38UeI+XdSdPnqRGjRraz3dXKN1PQEAAR48eZffu3Y98nyVLlmBra8vzzz+vs71t27a0bdtW+7l9+/Y0b96c7777jm+//faRr/84lJ8/Ie5h69atXLhwoUATOkDXrl1JTk5mxIgRREREcOLECd58802MjIzo3LkzADExMXh7e3PgwAHtefHx8YSHh3P27FkAjh07Rnh4ONeuXQOgfv36dOvWjaFDh7Jv3z727dvH0KFDefbZZ/Hy8noMTy30yqoatHsHhm2DkYfVOTzt60BupjqJ/O8vwZy6sP5duLgbivGHjRBCVEQxMTF07dqVK1eu0LRpU9avX4+FRflbZMPKygpra2tteViyOXLkSIKDg9m2bRuurq6PdA9FUVi0aBEDBw7ExOTBlRgGBga0atWKM2fOPPIzPC7lOtkMDAzE19eX+vXrF9jn7e3N+vXrOXr0KO3ateOJJ54gNjaWkJAQXFxcAMjJySEyMpKMjAztefPnz6dZs2YMHToUgI4dO9KsWTOCg4O1x/z22280atQIPz8//Pz8aNy4Mb/++mspP60oc+w91Tk8Aw7CsO3QdgRYVoOb1+HQEljSA+Y2hI2TIO6ojGgXQlR6SUlJdO3alUuXLlG3bl1CQkKwtbXVd1ilSlEUAgICWLNmDVu3bi3U+u47duzg7Nmz96xUu9d9wsPDtTlOWaJRFHkD6kN0dDRubm5ERUU98l84ohzIz1NrNI+thJPrISvlzj4HL2jcDxr2VVc5EkKISiQ9PZ3OnTsTFhZGjRo12LNnT7GmItSXwr6/33nnHYKCgli3bp1OC6iNjQ3m5uaAOog4JiaGpUuX6pw7cOBAzpw5w759+wpcd+rUqbRt25a6deuSmprKt99+y6+//sqePXu03fzKijLVZ1OIcs/AEGo/qZbuX8LZTWr/zsgQSIqErZ+qxbWV2r+zQW+wLDiVlhBCVDTjx48nLCwMe3t7Nm3aVC4TzaL48ccfAejUqZPO9sWLF2vnB4+Li+Py5cs6+1NSUli9ejXffPPNPa97/fp1hg0bRnx8PDY2NjRr1oydO3eWuUQTpGZTb6Rms5K5mQIRf94a0b4DlFt9OTUGULuTWttZ/1l1MJIQQlQwe/fupX379iiKQmhoKH5+fvoOqcjk/V14UrMpxONgZgPNXlFLWjycWAtHV6oj2s9tVcufY6CevzpxfF1/dfolIYQo57Kyshg8eDCKovD666+X60RTFI0km0I8blbVoO3barl67tZUSish6TREBKvF1BrqPwcNX4BaT6oTzgshRDn02WefERERgZOTE1999ZW+wxF6IM3oeiLV8EKHoqhrsB9fBcdWq0tl3lbFERr0UWs8XVupa7sLIUQ5cPz4cZo3b05OTg4rVqygf//++g6p2OT9XXhSXSJEWaDRgEtjtXSZAlH74Ngqtbk9PREOLFCLrbuadDbsC84+D72sEELoS15eHkOGDCEnJ4eePXvSr18/fYck9ESSTSHKGgMDcPdVyzOz4fx2dWBRxJ9w/RLs+lItTg1uJZ4vgF3lGNUphCg/vv/+e/bv34+1tTXz5s0rN2uei5InyaYQZZmhMdTtqpbsDDgdotZ4nt0ECSdgywnYMhXc2qi1nQ16g6WjvqMWQlRyFy9eZMKECQB8/vnnOss6ispHkk0hygsTC2jYRy2ZyRCx/tZUSrsgar9aQj5S5/hs1A+8nwUza31HLYSoZBRFYfjw4WRkZNCxY0ftinyi8pJkU4jyyNwOmr+mltQ4tW/nsT90p1IyHH1rKqV+UNdPplISQjwWv/76Kxs3bsTU1JSFCxdiYFCuV8YWJUCSTSHKO2sXaPeOWq6eg+Or1cTzv1MpeXVXa0VrdwYjE31HLYSogBISEhgzZgwAU6ZMoV69enqOSJQFMvWRnsjUCaJU3Z5K6dgf6jyed0+lZGarrlbUoI/M4SmEKFEDBgxgxYoVNG3alAMHDmBsbKzvkEqcvL8LT94yQlREd0+l9PRUiD6gJp0n/wc3rsCRZWqxsAefXmri6e6rru0uhBBFsH79elasWIGBgQE///xzhUw0RdFIsilERWdgADXbqqXbTLj0D5xYAyfXQcZVOLhILZbO4PO82tTu2lo9TwghHkFKSgpvv/02AO+99x4tWrTQc0SiLJFmdD2Ranihd3m5cGGHmnhGrIebKXf2WddQp1Fq0AdqNJdVi4QQD/T2228zf/58PD09OXr0KBYWFvoOqdTI+7vwJNnUE/lhFWVKbjac36Y2tZ/6C7LT7uyzdVcTz4Z9oFpjSTyFEDp27tzJk08+CcDWrVvp3LmzniMqXfL+LjxpRhdCqKPT6/mrJecmnN2s1nhGblBXLdozVy32ddTazoZ9wKm+vqMWQujZzZs3tfNoDhkypMInmqJoJNkUQugyNlNHq9d/Vl216EyoOp3SmU1w9Szs/FwtjvXVpLNBH3Coo++ohRB6MG3aNE6fPo2LiwtffPGFvsMRZZQ0o+uJVMOLcicrTa3pPL5GrfnMz7mzr1qjOzWedh56C1EI8fj8+++/tGjRgry8PNasWUPv3r31HdJjIe/vwpOaTSHEozG1gsb91ZJ5Xe3beWINnN+uzukZf0xdp92lKdR/Tp1SyaGunoMWQpSG3NxcBg8eTF5eHi+88EKlSTRF0UiyKYQoPHNbaPaKWjKuqasUHV8DF3dBXLhatk5Xm9p9eoFPT3DykcFFQlQQc+fO5dChQ9ja2vL999/rOxxRxkkzup5INbyokG4kQuRfcDJYnVYpP/fOvqqeatJZ/zmoLtMpCVFenTt3jkaNGpGZmUlgYCCDBg3Sd0iPlby/C0+STT2RH1ZR4WUmQ2SIWut5dgvkZd3ZZ+OmJp31e4JbG5lAXohyQlEUnn76abZu3cpTTz3F5s2b0VSyPxzl/V140owuhCgd5nbQ9CW1ZKXBmY1qjeeZTZASBfvmqcXSGbyfVWs93TvIWu1ClGGLFi1i69atmJub89NPP1W6RFMUjfxWF0KUPlMraPiCWnIy1ZrOiGC15vPGFTgYqBbzquDdHer3gtqd1Pk/hRBlQlxcHO+99x6gTnnk6emp54hEeSHJphDi8TI2vzOPZ2622rfz5Dp1dHvmNTiyTC2m1lCvm1rjWedp9TwhhN4EBASQkpJCixYtGD16tL7DEeWIJJtCCP0xMoG6XdXy7Fy4tEet8YxYr9Z4HlupFmML9Zj6PdVVjkyt9B25EJXKmjVrWLNmDUZGRgQGBmJkJOmDeHTy0yKEKBsMjaD2k2p55guIPqD28YwIVvt4nlynFkNT8HxKHWBUzx+qOOg7ciEqtOTkZEaMGAHABx98QJMmTfQckShvJNkUQpQ9BgZQs61a/GdA7BE16TwZDNfOwekNakGjjmb3ekYtDvVkSiUhSti4ceOIj4/Hy8uLSZMm6TscUQ7J1Ed6IlMnCFEEigIJJ9WkM/JviD+qu79qbfDqrvb1rNlORrYLUUxbt26lS5cuAOzcuZMnnnhCzxHpn7y/C09+Ewshyg+NBpwbqKXzeLgeBadD1HJhJ1w7D3u/V4uZLdT1A69u6gAjMxt9Ry9EuZKRkcHQoUMBePvttyXRFEUmyaYQovyydYPWQ9WSlQbntkLkBjgdqo5svz3AyMAI3NurtZ5e3cDOQ9+RC1HmTZ48mfPnz+Pq6sqsWbP0HY4ox6QZXU+kGl6IUpSfB1EH1Kb20yGQdFp3v5PPrX6e3dWlM2UFIyF0HDp0iNatW5Ofn8/69et59tln9R1SmSHv78KTmk0hRMVjYAju7dTiNx2SzqoDiiJD4PI/ar/PhJOw60uo4qSOavfqrk4kb2Kh7+iF0KucnBwGDx5Mfn4+AwYMkERTFJskm0KIis+hDjiMBN+RkHFNXTLz9AY4sxnSE+DIr2oxMlMTTq9n1EFGVtX0HbkQj92cOXP4999/qVq1Kt98842+wxEVgCSbQojKxaIqNHlRLbnZcGm3WuMZuQFSLt8ZcARqE7vXM+qE8tWaSHO7qPAiIyOZOnUqAF9//TVOTk56jkhUBNJnU0+kz4cQZYyiwJUTt5rbN0DMId39Fg7g2Rk8u6iTyls56ydOIUpJfn4+nTp1YteuXfj5+RESEoJG5q0tQN7fhSc1m0IIAeq0StUaqqXjOEiLV0e1355WKSMJjv2hFgDnRlDnKTX5rNkWjEz1G78QxfTTTz+xa9cuqlSpwoIFCyTRFCVGkk0hhLgXq2rQ4nW15Gary2ee3QLntkDcv3DlmFr2fAPGVaDWE2riWaeLOrm8vKhFORIdHc0HH3wAwIwZM/Dw8NBvQKJCkWRTCCEexsgEPDqo5enJcCMRzm+7lXxuVQcZ3d3X09ZdTTo9u0CtjmBmrd/4hXgARVF45513SEtLo02bNgQEBOg7JFHBSLIphBCFZekIjfurJT8frhxXazzPboHL++D6JTi4SC0GRuDa+k6Tu0tTGWgkypSVK1eyfv16jI2NCQwMxNDQUN8hiQpGkk0hhCgOAwNwaayWDmMg6wZc3HWnyf3aeXVuz8v/wNZPwcJeHWAkA41EGZCSksK7774LwIQJE2jQoIGeIxIVUbn883r79u1oNJp7lrCwMO1xYWFhdOnSBVtbW+zs7PDz8yM8PPyB1/7pp5/o1KkT1tbWaDQarl+/XuAYDw+PAvf96KOPSvgphRDlkqmlOl1Sjzkw6giMCoceX4JXDzCxgoyr6iCj/70FX9aDHzvApk/g/A7IzdJ39KKSmTJlCleuXKFevXqMHz9e3+GICqpcTn2UnZ3NtWvXdLZNmjSJzZs3c/78eTQaDWlpabi7u9OrVy8++ugjcnNzmTx5Mrt27SI6OhpjY+N7Xnvu3LncvHkTgPHjx5OcnIytra3OMR4eHgwePJihQ4dqt1laWmJpafnIzyBTJwhRCeXlqMto3m5yjwvX3W9kDm6twL0DeLSHGi3B2EwvoYqK79ixYzRr1oy8vDxCQ0Px8/PTd0jlgry/C69cNqObmJhQrdqdlT1ycnIIDg4mICBAO1VDZGQkycnJTJs2DTc3NwAmT55M48aNuXz5Mp6enve89ujRowG19vRBrKysdGJ4mKysLLKy7tRapKWlPfK5QogKwtBYTSI92kOXTyA9Cc5tU5PPc1vhxhV1mqULO28dbwqurdTj3duDW2swNtfvM4gKQVEUAgICyMvLo3fv3pJoilJVLpvR/ys4OJikpCTeeOMN7TYvLy8cHBwIDAwkOzubzMxMAgMDadCgAe7u7sW+5+zZs7G3t6dp06bMmDGD7OzsBx4/c+ZMbGxstMXHx6fYMQghyrkqDtC4H/SeD+9Fwjv71Sb3Bn3A0hnystQVjnbMhqU9YVZNWNQNtkxXk9PsdH0/gSinfv/9d3bu3Im5uTlff/21vsMRFVy5bEb/r+7duwPw999/62w/ceIEvXr14sKFCwDUq1eP0NBQatas+dBrbt++nc6dO9+zGf3rr7+mefPm2NnZceDAAcaPH0+vXr34+eef73u9/9ZsxsTE4OPjI9XwQoh7UxS4ehYu7oZLe9R/0+J0jzEwUpfU9GivNr3XbAOmVvqJV5QbaWlpeHt7Exsby/Tp0/n444/1HVK5Is3ohVemks0pU6Zo12S9n7CwMFq2bKn9HB0djbu7OytXruSFF17Qbs/MzKRTp054e3trmwrmzJnDqVOnCAsLw9z8wU1RD0o2/2v16tX07duXpKQk7O3tH/6gyA+rEKKQFEUd2X5pD1y8lXymRuseozGE6k3VJnePJ9SVjWSOT/Ef48aNY86cOXh6enL8+HHMzKRfcGHI+7vwylSfzYCAAAYMGPDAY/67qsHixYuxt7enZ8+eOtuDgoK4ePEie/fuxeDWnHZBQUHY2dmxbt26h96nMNq2bQvA2bNnHznZFEKIQtFowN5TLc1fU5PP65fuJJ6XdsP1y+qa7jGH4J9vQWMA1RrfmZC+Zjswt9X3kwg9ioiIYO7cuQB88803kmiKx6JMJZsODg44ODg88vGKorB48WJee+21AqPLMzIyMDAw0Fnb9fbn/Pz8EosZ4MiRIwC4uLiU6HWFEOK+NBqw81BLs1fUbdejbtV87lKT0OQL6oj3uHDY+z1wa/139w7qYKMaLcC2piytWUkoisLIkSPJzc3lueeeo0ePHvoOSVQSZSrZLKytW7dy4cIFBg8eXGBf165dGTduHCNGjGDkyJHk5+cza9YsjIyM6Ny5M6D2m+zSpQtLly6ldevWAMTHxxMfH8/Zs2cBdWoIKysratasSdWqVdm7dy/79u2jc+fO2NjYEBYWxpgxY+jZs+cj9QUVQohSY+sGtgOgya2Wm9TYWzWfu9Qk9OpZiD+mlv0/qsdYOKhJp7Y0B4uq+nsGUWpWrVrFli1bMDU11dZuCvE4lOtkMzAwEF9fX+rXr19gn7e3N+vXr2fq1Km0a9cOAwMDmjVrRkhIiLYGMicnh8jISDIyMrTnzZ8/X6ffaMeOHQG1uf6NN97A1NSUFStWMHXqVLKysnB3d2fo0KF88MEHpfy0QghRSNbV1dHujfupn9Pi1aTz0j9qU3v8cchIgjOharnNrha4tryTgFZrJFMulXPp6emMHTsWgA8//JDatWvrOSJRmZSpAUKViXQwFkLoXc5NdV332/08Yw6ptZ//ZWAEzg10a0Ad6oGBrKFdXkyYMIGZM2fi4eHByZMnHzpIVtyfvL8Lr1zXbAohhCgGYzO1BtP1zgwfZCZDzOFb5RDEHIT0RIj7Vy0HF6nHmVhC9Wa6Cah1den/WQadPn2aOXPmAOrUfZJoisdNkk0hhBB3mNtBnS5qAXXUe0r0XbWfhyH2CGTfuDUQadedcy2r3en3WaOFmozK6He9UhSFUaNGkZOTQ7du3ejVq5e+QxKVkCSbQggh7k+juTXwyA0aPK9uy8+DxMg7NZ8xh+DKSbgRD5F/qeW2qrXByQccvcGpvlrs64KRiV4ep7JZt24doaGhmJiY8O233+rM0CLE4yLJphBCiMIxMARnH7U0H6huy86A+KO6/T+TL6oT0V87D6f+vHO+xhDs64CTt24iWrW2un68KBEZGRmMHj0agPfff5+6devqNyBRaUmyKYQQovhMLNQVi2q2vbMtPUkdgJRwChJOQuIpSIiArFRIilTLyXV3jjc0UWs9nW4ln463akLtPGQwUhHMmjWLS5cu4ebmxoQJE/QdjqjEJNkUQghROqo4QO1OarlNUdT5PxMj1MQz4dStr09BTjoknFDL3YzM1NHvt5vhbyehNm5wa4U4oevcuXN8/vnnAHz11VdUqVJFzxGJykySTSGEEI+PRgM2NdRS5+k72/PzISVKTUBvJ58JJyHpNOTeVJvo44/qXsu4Cjh63UpAvaFqLTUBtXFTJ6avxP0TR48eTVZWFk8//TQvvPCCvsOp1GbOnMmaNWs4deoU5ubm+Pr6Mnv2bLy8vO57zhtvvMEvv/xSYLuPjw8nTtz5Y2z16tVMmjSJc+fO4enpyYwZM+jdu3epPEdxSLIphBBC/wwMwM5dLV7d7mzPz1P7fuokoRFw9YxaExp7WC3/ZVwFbFzVgU02rmoCalvzVjLqClYuYFgxX4F//vknf/75J0ZGRjIoqAzYsWMHI0aMoFWrVuTm5jJx4kT8/Pw4efLkfWucv/nmG2bNmqX9nJubS5MmTejXr5922969e3nxxReZPn06vXv3Zu3atfTv35/du3fTpk2bUn+uwpBJ3fVEJoUVQohiyMtVBx7d7guaeAquX1bXh09PePj5GkOwrvGfZNTtTs2ojavaD7WcuXnzJg0aNOD8+fOMGzdO25QuSk5x39+JiYk4OTmxY8cO7SqFD/O///2PPn36cOHCBdzd3QF48cUXSU1NZcOGDdrjunXrhp2dHcuXLy90XKWpYv5ZJ4QQomIzNALHemr5r5ybkBqjJp8pUeo8odej1K+vX1b35edCymW13I+Fw13JaM07yaiVizp/qJktmNmUqRrSL774gvPnz1O9enUmTZqk73AqtLS0NFJTU7WfTU1NMTU1feh5KSkpAFStWvWR7xUYGMjTTz+tTTRBrdkcM2aMznH+/v5lct37svN/iBBCCFESjM3A3lMt95KfBzeu3ElAU6Lu+vpWYpqdpq4bn5GkTmL/ICZWatJ5OwG9OxHV2Waju9/cFowenpw8qosXL/LZZ58BMGfOHKysrErs2nqXnwe5WZCXDXk5kHfX17lZ99mWDR4dwNKpVELy8fHR+Tx58mSmTJnywHMURWHs2LF06NCBhg0bPtJ94uLi2LBhA0FBQTrb4+PjcXZ21tnm7OxMfHz8I133cZJkUwghROViYKgurWldHbhH3zZFgZvXbyWg0XdqRG8no2lX4GaKmpCC+m92GqRGFz4WI7P7J6hmNurk9wZGarO/gZEau4HhXZ/vbBs7+gtu3rzJk60aMaCFPZzdctf+W8dqDHS3aQzvTCuVn6smavk5t5K3219nq90WCmz/z9cFtuWq5xb4+tY1tUlitm7JzaZAUqnkF+E/NDBwLVg+VbRzH+LkyZPUqFFD+/lRajUDAgI4evQou3fvfuT7LFmyBFtbW55//vkC+/7bH1dRlDLZR1eSTSGEEOJuGo26bKe5Hbg0vv9xeblq0nnzOmReV//V+TrlztcFtqUAijrS/ka8Wooh9Gwua7dkYKiB71ueRxPU7+EnlWeGJgWL0e2vjcHQVP3a1LrUQrCyssLa+tGvP3LkSIKDg9m5c+cj9/VUFIVFixYxcOBATEx0V92qVq1agVrMhISEArWdZYEkm0IIIURRGBpBFXu1FFZ+vjq5/X+T1duJ6O3ENC9bPTY/F5Q89d/8XJ1tWVlZjJy/HYCRXWrSsJmn2uycn/efc/677a7PiqI+j4HxrWTN+M7XBbYZqYnc7a8NjG8leXedb2B0J/ErcI6x2n1AJzG81zaT/2y/K5ksg7V396MoCiNHjmTt2rVs376dWrVqPfK5O3bs4OzZswwePLjAvnbt2rFp0yadfpsbN27E19e3ROIuSZJsCiGEEI+bgYHaXG5uC7g/5OAH+3rWLM5c+QtnZ2emrDoKNjYlEaEoISNGjCAoKIh169ZhZWWlrY20sbHB3NwcgPHjxxMTE8PSpUt1zg0MDKRNmzb37N/57rvv0rFjR2bPnk2vXr1Yt24dmzdvLlQT/eMiSy8IIYQQ5VRUVBTTp08H1JHoNpJoljk//vgjKSkpdOrUCRcXF21ZsWKF9pi4uDguX9adGSElJYXVq1ffs1YTwNfXl99//53FixfTuHFjlixZwooVK8rcHJsg82zqjcyzKYQQorj69+/PH3/8QYcOHdi5c2eZHBxS0cj7u/CkZlMIIYQohzZv3swff/yBgYEB33//vSSaosySZFMIIYQoZ7Kzsxk5ciQA77zzDk2aNNFzRELcnySbQgghRDnz7bffcurUKRwdHbV9NoUoqyTZFEIIIcqR2NhYpk6dCsCsWbOwtbXVb0BCPIRMfSSEEEKUI+PGjePGjRu0adOGN954Q9/hiAoiJyeH+Ph4MjIycHR0LNTa7Q8jNZtCCCFEObFjxw6CgoLQaDT88MMPGBjIa1wU3Y0bN1iwYAGdOnXCxsYGDw8PfHx8cHR0xN3dnaFDhxIWFlbs+8hPqRBCCFEO5OTkEBAQAMCwYcNo0aKFniMS5dnXX3+Nh4cHCxcu5KmnnmLNmjWEh4cTGRnJ3r17mTx5Mrm5uXTt2pVu3bpx5syZIt9LmtGFEEKIcmDevHkcP36cqlWrMmPGDH2HI8q5f/75h23bttGoUaN77m/dujWDBg1i/vz5BAYGsmPHDurWrVuke0myKYQQQpRx8fHxfPLJJwDMnDkTe/sirMcuxF3++OOPRzru6NGjvPPOO8W6lzSjCyGEEGXchx9+SGpqKi1atLjv8oVClJSEhAS++uorGjZsSNu2bYt9PUk2hRBCiDJsz549LF26FIAffvgBQ0NDPUckKqK8vDzWrVvH888/j5ubGwsXLuT555/n4MGDxb62NKMLIYQQZVReXp52UNDgwYNp06aNniMS5V1iYiJfffUVVatWZfTo0URGRrJ48WKWLVsGQP/+/cnPz2f16tX4+PiUyD2lZlMIIYQoo+bPn094eDi2trbMnDlT3+GICuDll18mIyMDgBo1atC2bVtiY2NZtGgRsbGxfPfddyV+T0k2hRBCiDIoMTGRjz/+GIBPP/0UR0dHPUckKoJTp07xyiuvMGjQIK5du8awYcOYNm0aPXr0KLUuGpJsCiGEEGXQ+PHjuX79Ok2bNuWtt97Sdziigvj444/p3bs3Tz75JLNmzeLixYs0bNiQNm3a8P3335OYmFji95Q+m0IIIUQZs2/fPgIDAwH4/vvvZVCQKDHDhw/nlVdewdTUFGNjY0CtRV+2bBkLFy5kzJgx5Ofns2nTJtzc3LCysir2PaVmUwghhChDsrOzGTp0KACvv/467du313NEoqKxtLTUJpoAjo6OjBkzhn///Zd9+/bx9ttvM336dJycnOjZs2ex7yfJphBCCFGGzJo1i+PHj+Pg4MCcOXP0HY6oZFq0aMH3339PbGwsv/zyC7m5ucW+piSbQgghRBlx8uRJPv30UwC+/fZbHBwc9ByRqKxMTEzo378/f//9d7GvJcmmEEIIUQbk5eUxZMgQcnJy6NGjBwMGDNB3SKICu3z5cqGOj4mJKfK9JNkUQgghyoB58+axd+9erKys+PHHH9FoNPoOSVRgrVq1YujQoRw4cOC+x6SkpLBw4UIaNmzImjVrinwvGY0uhBBC6NmlS5cYP348oPbZdHNz03NEoqKLiIjgs88+o1u3bhgbG9OyZUuqV6+OmZkZycnJnDx5khMnTtCyZUu++OILnnnmmSLfS6MoilKCsYtHFB0djZubG1FRUbi6uuo7HCGEEHqiKArdu3cnJCSEDh06sGPHDgwMpOGxrKpo7++bN2/y999/s2vXLi5evEhmZiYODg40a9YMf39/GjZsWOx7SM2mEEIIoUe//fYbISEhmJiYsHDhQkk0xWNlZmZGnz596NOnT6ndQ36ihRBCCD1JTExk9OjRAHzyySd4e3vrNyAhSkG5TDa3b9+ORqO5ZwkLC9MeFxYWRpcuXbC1tcXOzg4/Pz/Cw8Pve91r164xcuRIvLy8sLCwoGbNmowaNYqUlBSd45KTkxk4cCA2NjbY2NgwcOBArl+/XkpPK4QQoqIaPXo0V69epVGjRnzwwQf6DkeIUlEuk01fX1/i4uJ0ypAhQ/Dw8KBly5YApKWl4e/vT82aNdm/fz+7d+/G2toaf39/cnJy7nnd2NhYYmNjmTNnDseOHWPJkiWEhIQwePBgneNefvllwsPDCQkJISQkhPDwcAYOHFjqzy2EEKLi+OuvvwgKCsLAwIDAwECdFV2EqEgqxAChnJwcXF1dCQgIYNKkSQAcPHiQVq1acfnyZe2ovmPHjtG4cWPOnj2Lp6fnI137jz/+4NVXXyU9PR0jIyMiIiLw8fFh3759tGnTBlDXsG3Xrh2nTp3Cy8vrntfJysoiKytL+zkmJgYfH58K08FYCCHEo0tLS6NBgwZERUUxduxYvvzyS32HJB5RRRsg9DiUy5rN/woODiYpKYk33nhDu83LywsHBwcCAwPJzs4mMzOTwMBAGjRogLu7+yNfOyUlBWtra4yM1LFUe/fuxcbGRptoArRt2xYbGxv++eef+15n5syZ2mZ3GxsbfHx8Cv+gQgghKoTx48cTFRVFrVq1mDZtmr7DEaJUVYjR6IGBgfj7++vMS2ZlZcX27dvp1asX06dPB6BevXqEhoZqE8eHuXr1KtOnT2f48OHabfHx8Tg5ORU41snJifj4+Ptea/z48YwdO1b7+XbNphBCiMplz549zJs3D4CffvqJKlWq6DkiIWDLli1s2bKFhIQE8vPzdfYtWrSoWNcuUzWbU6ZMue/An9vl4MGDOudER0cTGhpaoF9lZmYmgwYNon379uzbt489e/bQoEEDunfvTmZm5kNjSU1NpUePHvj4+DB58mSdffda1UFRlAeu9mBqaoq1tbW2WFlZPTQGIYQQFcvNmzcZMmQIiqLw5ptv8vTTT+s7JCGYOnUqfn5+bNmyhaSkJJKTk3VKcZWpms2AgICHrgXr4eGh83nx4sXY29vTs2dPne1BQUFcvHiRvXv3aucsCwoKws7OjnXr1j3wPmlpaXTr1g1LS0vWrl2r02m7WrVqXLlypcA5iYmJODs7P+wRhRBCVGIzZszg1KlTODs7Sz9NUWbMnz+fJUuWlNpg52Ilmzk5OcTHx5ORkYGjoyNVq1YtVjAODg44ODg88vGKorB48WJee+21AqP4MjIyMDAw0KltvP35v9XDd0tNTcXf3x9TU1OCg4MxMzPT2d+uXTtSUlI4cOAArVu3BmD//v2kpKTg6+v7yLELIYSoXI4dO8asWbMA+P7777Gzs9NzREKosrOzSzWHKXQz+o0bN1iwYAGdOnXCxsYGDw8PfHx8cHR0xN3dnaFDh+rMdVmatm7dyoULFwo0oQN07dqV5ORkRowYQUREBCdOnODNN9/EyMiIzp07A2q/SW9vb+0i9Glpafj5+ZGenk5gYCCpqanEx8cTHx9PXl4eAPXr16dbt24MHTqUffv2sW/fPoYOHcqzzz5735HoQgghKre8vDyGDBlCbm4uzz//PC+88IK+QxJCa8iQIQQFBZXa9QtVs/n1118zY8YMPDw86NmzJx999BE1atTA3Nyca9eucfz4cXbt2kXXrl1p27Yt3333HXXr1i2t2AkMDMTX15f69esX2Oft7c369euZOnUq7dq1w8DAgGbNmhESEoKLiwug1sxGRkaSkZEBwKFDh9i/fz8AderU0bnehQsXtE34v/32G6NGjcLPzw+Anj178v3335fWYwohhCjnvv32Ww4cOIC1tTU//PDDA/v4C/E43D1oOT8/n59++onNmzfTuHHjAq3FX331VbHuVah5Nvv168cnn3xCo0aNHnhcVlYWgYGBmJiYMGTIkGIFWFHJPF1CCFE5XLhwgYYNG5KRkcGCBQsYNmyYvkMSxVBR3t+3W3kfRqPRsHXr1mLdq0JM6l4eVZQfViGEEPenKAp+fn5s3ryZJ598kq1bt2oHrYryqSK+vy9fvoyrq2uBn01FUYiKiqJmzZrFur78xAshhBCl5JdffmHz5s2YmZmxcOFCSTRFmVSrVi2SkpIKbL927Rq1atUq9vVLbOqjU6dOsX79emxtbWnQoAENGzbE2tq6pC4vhBBClCtXrlzR9oubMmVKqY5hEKI47tfIfePGjQKz8hRFiSWbzzzzDEOGDOH69essWLCAY8eOkZ6eTmRkZEndolLKy8sjJydH32FUCMbGxhgaGuo7DCFEJTFy5EiSk5Np1qwZ7733nr7DEaKA238MaTQaPvnkEywsLLT78vLy2L9/P02bNi32fUos2XRxcWHixIk6225PFyQKT1EU4uPjuX79ur5DqVBsbW2pVq2ajAQVQpSqdevW8ccff2BoaEhgYOAjL5MsxON05MgRQM05jh07homJiXafiYkJTZo04f333y/2fYr90//ee+/RpEkTOnfuzC+//MLrr7+u3Se1SEV3O9F0cnLCwsJCkqNiUhSFjIwMEhISALTTXwkhRElLSUnhnXfeAeD999+nWbNmeo5IiHvbtm0bAG+++SbffPNNqXV/LHay2bFjR44ePcrp06dZtWoVM2bMoGXLljRq1IhGjRrx7LPPlkSclUpeXp420bS3t9d3OBWGubk5AAkJCTg5OckfQ0KIUvHhhx8SGxtL3bp1mTx5sr7DEeKhFi9eXKrXL3ay2atXL3r16qX9nJmZyfHjxzl69ChbtmyRZLMIbvfRvLvvhCgZt7+nOTk5kmwKIUrcjh07WLBgAQALFy7U/pErRFl29wTvd9NoNJiZmVGnTh169epV5GXJC51sDhw4kAULFtw3ETI3N6dVq1a0atWqSAGJO6TpvOTJ91QIUVoyMzMZOnQoAMOGDePJJ5/Uc0RCPJojR45w+PBh8vLy8PLyQlEUzpw5g6GhId7e3sybN4/33nuP3bt34+PjU+jrF3rCr6CgIG7cuKH9PHz4cJKTk3WOkdHTQgghKptp06Zx5swZXFxcmD17tr7DEeKR9erVi6effprY2FgOHTrE4cOHiYmJoWvXrrz00kvExMTQsWNHxowZU6TrFzrZ/O9cTMuXL9dJNq9cuYKVlVWRghFCCCHKoyNHjvDFF18AMG/ePGxtbfUbkBCF8MUXXzB9+nSdAULW1tZMmTKFzz//HAsLCz755BMOHTpUpOsXeymDe00Emp2dXdzLCiGEEOVCbm4uQ4YMIS8vj759+/L888/rOyQhCiUlJUU7W8vdEhMTSU1NBdSpA4ua35XKulnSL04UVu/evbGzs6Nv3776DkUIIQrlq6++4vDhw9jZ2fHdd9/pOxwhCq1Xr14MGjSItWvXEh0dTUxMDGvXrmXw4MHaP54OHDhAvXr1inT9IiWbQUFBHD58WNs3U5JLUVyjRo1i6dKl+g5DCCEK5ezZs9rpjb788kuqVaum54iEKLwFCxbQpUsXBgwYgLu7OzVr1mTAgAF06dKF+fPnA+Dt7c3PP/9cpOsXejR6hw4dmDx5MmlpaRgbG5Obm8uECRPo0KEDzZs3x9HRsUiBiMqtc+fObN++Xd9hCCHEI1MUhWHDhnHz5k2efvpp3njjDX2HJESRWFpasnDhQr7++mvOnz+Poih4enpiaWmpPaY4y1YWumZz586dpKSkEBkZyS+//MJ7771HQkICkyZNon379nh5eRU5GFH+5efn4+3tzYcffqizPTQ0FBMTE/74448iX/vIkSMYGxvzxBNPFDdMIYQotsDAQLZt24aFhQULFiyQVj5R7llaWtK4cWOaNGmik2gWV5Enda9bty5169ZlwIAB2m0XLlzg4MGD2rU2ReVjYGDA+PHjGTlyJB999BF2dnb8+++/9OvXj88++4x+/foV+dqjRo3i/fff54cffkBRFPnFLoTQm9jYWO2a0dOnT6d27dp6jkiI4tmyZQtbtmwhISGB/Px8nX2LFi0q1rULVbN5+fLlB+6vVauWNqkAiImJKXpkotx65ZVXcHBw4NtvvyU6OpoePXowcOBA7S/moggKCsLOzo4RI0aQlpbG+fPnSzBiIYQonICAAFJSUmjVqhXvvvuuvsMRolimTp2Kn58fW7ZsISkpieTkZJ1SXIWq2WzVqhU9e/Zk6NChtG7d+p7HpKSksHLlSr755huGDx/OyJEjix1kZacoCpk5eXq5t7mxYaFrEI2MjPjwww+ZMGECq1evpnnz5nz77bdFjiE9PZ0JEyawYcMGXF1dsbGxITw8HE9PzyJfUwghimr16tWsXbsWIyMjfv75Z1n6VpR78+fPZ8mSJQwcOLBUrl+oZDMiIoLPPvuMbt26YWxsTMuWLalevTpmZmYkJydz8uRJTpw4QcuWLfniiy945plnSiXoyiYzJw+fT0L1cu+T0/yxMCl8b4tXXnmF0aNHoygKy5cvf+gvY39/fw4fPkx6ejqurq6sXbtWu+TpjBkz6NatG/Xr1wfAx8eH8PBwXnjhhcI/kBBCFENycjIBAQEAfPjhhzRu3FjPEQlRfNnZ2fj6+pba9QvVjF61alXmzJlDbGwsP/74I/Xq1SMpKYkzZ84AaoJx6NAh9uzZI4lmJXf7l3FSUtIj/dUfGhpKYmIiGRkZREdHaxPN8+fP89NPPzFlyhTtsQ0bNiQ8PFznfJmnUwjxOLz33nvEx8fj5eXFxx9/rO9whCgRQ4YMISgoqNSuX6QBQmZmZvTp04c+ffqUdDziHsyNDTk5zV9v9y6sSZMm8ddff7Fv3z6efvppAgMDGTFiRJHuP2bMGK5evYqrq6t2W35+PjVq1NA5btSoUQwaNIhffvmlSPcRQoiHWb58OYsXL0aj0fDzzz9jZmam75BEOTBz5kzWrFnDqVOnMDc3x9fXl9mzZz909p6srCymTZvGsmXLiI+Px9XVlYkTJzJo0CAAlixZwptvvlngvMzMzEL/bN68eZOffvqJzZs307hxY4yNjXX2f/XVV4W63n8VeTS6eHw0Gk2RmrL14eeff+bLL79k69atNGnShNGjR/P5558zbNiwAj+8D7Np0yb27NnDkSNHMDK68/xhYWEMGjSIq1evYm9vD8g8nUKI0nX69GmGDRsGwMSJE+nQoYOeIxLlxY4dOxgxYgStWrUiNzeXiRMn4ufnx8mTJ6lSpcp9z+vfvz9XrlwhMDCQOnXqkJCQQG5urs4x1tbWREZG6mwryh9BR48e1c6jefz4cZ19JTHzS/nIYES5sGHDBkaMGMHy5ctp27YtACNHjmTOnDn8+uuv2r/GHkVubi7vvvsu48aNKzCRrLW1NQDh4eF06dKlxOIXQoh7yczMpF+/fty4cYMnn3xSu2KQEI8iJCRE5/PixYtxcnLi0KFDdOzY8b7n7Nixg/Pnz1O1alUAPDw8Chyn0WhKZNWqbdu2FfsaD1Iqa6OLyufQoUP069ePzz//XKd7hbW1NSNHjmTWrFnk5T36iPrvvvuOq1evavt+3s3NzQ0LC4sC/TaFEKI0jB49mqNHj+Lo6EhQUJBOS4uovNLS0khNTdWWrKysRzovJSUFQJtE3ktwcDAtW7bk888/p0aNGtSrV4/333+fzMxMneNu3LiBu7s7rq6uPPvss8Wa53zXrl28+uqr+Pr6aqeu/PXXX9m9e3eRr3mbJJuiRLRo0YIbN27cc765adOmcfr06UJNDzJmzBiuXLlyzyYGjUZDeno67733XrFiFkKIhwkKCuKnn35Co9Hw22+/Ub16dX2HJMoIHx8fbGxstGXmzJkPPUdRFMaOHUuHDh1o2LDhfY87f/48u3fv5vjx46xdu5a5c+eyatUqnfEP3t7eLFmyhODgYJYvX46ZmRnt27fXDtoujNWrV+Pv74+5uTmHDx/WJs5paWnaudOLQ/48E0IIIe4hMjKS4cOHA/Dxxx/TtWtXPUckypKTJ0/qDFY1NTV96DkBAQEcPXr0obWF+fn52j9wbGxsAHWQTt++ffnhhx8wNzenbdu22i5rAO3bt6d58+Z89913hZ7b+tNPP2X+/Pm89tpr/P7779rtvr6+TJs2rVDXupci12y+8cYb7Ny5s9gBCFES/P396devH3///Teurq6EhYXpOyQhRDl2dz/NTp06ST9NUYCVlRXW1tba8rBkc+TIkQQHB7Nt2zadGVbuxcXFhRo1amgTTYD69eujKArR0dH3PMfAwIBWrVoVqWYzMjLynv1Hra2tuX79eqGvVyC2op6YlpaGn58fdevW5bPPPpOlKYVe3W+eTiGEKIp3332XY8eO4eTkRFBQkKwSJIpMURQCAgJYs2YNW7dupVatWg89p3379sTGxnLjxg3tttOnT2NgYHDfRFVRFMLDw3FxcSl0jC4uLpw9e7bA9t27d1O7du1CX++/ipxsrl69mpiYGAICAvjjjz/w8PDgmWeeYdWqVeTk5BQ7MCGEEEIffvvtNxYuXKhtxizKy1uI20aMGMGyZcsICgrCysqK+Ph44uPjdQb7jB8/ntdee037+eWXX8be3p4333yTkydPsnPnTsaNG8egQYMwNzcH1PXMQ0NDOX/+POHh4QwePJjw8HDeeuutQsc4fPhw3n33Xfbv349GoyE2NpbffvuN999/n3feeaf43wSlhBw+fFgJCAhQzMzMFAcHB2X06NHK6dOnS+ryFU5UVJQCKFFRUQX2ZWZmKidPnlQyMzP1EFnFJt9bIcSDREREKFWqVFEA5ZNPPtF3OKIMetD7+16Ae5bFixdrj3n99deVJ598Uue8iIgI5emnn1bMzc0VV1dXZezYsUpGRoZ2/+jRo5WaNWsqJiYmiqOjo+Ln56f8888/RX6uCRMmKObm5opGo1E0Go1iZmamfPzxx0W+3t00iqIoxU1Y4+LiWLp0KYsWLSImJoYXXniBuLg4tm3bxueff86YMWOKe4sKJzo6Gjc3N6KiogpUid+8eZMLFy5Qq1YtWaGihMn3VghxPxkZGbRp04bjx4/TuXNnNm3aJM3nooAHvb/Lu4yMDE6ePEl+fj4+Pj5YWlqWyHWLPBo9JyeH4OBgFi9ezMaNG2ncuDFjxozhlVdewcrKCoDff/+dt99+W5JNIYQQZd67777L8ePHcXZ2ln6aolKysLCgZcuWJX7dIiebLi4u5Ofn89JLL3HgwIECq7yAOkLY1ta2GOEJIYQQpW/ZsmX8/PPP2n6aJbEqixBl2dixYx/5WL2tjf7111/Tr1+/BzZF2tnZceHChaLeQgghhCh1ERER2vk0J0+eLMvgikrhUVcb0uva6AMHDiz2zYUQQgh9ysjIoH///mRkZPDUU0/x8ccf6zskIR6L0l4P/W5FTjbvV/2q0WgwMzOjTp069OrV64FrfwohhBD6NHLkSG0/zd9++036aQpRCoqcbB45coTDhw+Tl5eHl5cXiqJw5swZDA0N8fb2Zt68ebz33nvs3r0bHx+fkoxZCCGEKLbbs6gYGBiwfPly6acpRCkp8qTuvXr14umnnyY2NpZDhw5x+PBhYmJi6Nq1Ky+99BIxMTF07NhRRqILIYQocyIiInj77bcBtZ9m586d9RyREBVXkZPNL774gunTp2Ntba3dZm1tzZQpU/j888+xsLDgk08+4dChQyUSqBBCCFES0tPT6devHxkZGTz99NNMnDhR3yEJUaEVOdlMSUkhISGhwPbExERSU1MBsLW1JTs7u+jRiUqhd+/e2NnZ0bdvX32HIoSoBEaOHMmJEyeoVq0ay5Ytk36aQpSyYjWjDxo0iLVr1xIdHU1MTAxr165l8ODBPP/88wAcOHCAevXqlVSsWtu3b0ej0dyzhIWFaY8LCwujS5cu2NraYmdnh5+fH+Hh4fe97rVr1xg5ciReXl5YWFhQs2ZNRo0aRUpKis5xHh4eBe770UcflfhzVhajRo1i6dKl+g5DCFEJ/PLLLyxevBgDAwOCgoJwdnbWd0hCVHhFHiC0YMECxowZw4ABA8jNzVUvZmTE66+/ztdffw2At7c3P//8c8lEehdfX1/i4uJ0tk2aNInNmzdrZ75PS0vD39+fXr16MW/ePHJzc5k8eTL+/v5ER0djbGxc4LqxsbHExsYyZ84cfHx8uHTpEm+99RaxsbGsWrVK59hp06YxdOhQ7eeSWtKpMurcuTPbt2/XdxhCiAru5MmTvPPOOwBMmTJF+mmKSq3MT+qek5PDc889x4IFC/j66685f/48iqLg6empk3Tda1WhkmBiYqIzavD20pkBAQHayUcjIyNJTk5m2rRpuLm5AWon8MaNG3P58mU8PT0LXLdhw4asXr1a+9nT05MZM2bw6quvkpubi5HRnW+XlZWVjFy8h8aNG9O7d2+mTp1aYN9nn33GF198wenTp3F0dHyk63Xs2JFdu3ZpP9vZ2dGxY0cWLlz4yNcQQoi7+2l27dqVCRMm6DskIfTqv5O6Hzp0SDvDEMDp06cxNDSkRYsWxb5XkZrRjY2NOX78OBqNBktLSxo3bkyTJk30VrsXHBxMUlISb7zxhnabl5cXDg4OBAYGkp2dTWZmJoGBgTRo0AB3d/dHvnZKSgrW1tY6iSbA7Nmzsbe3p2nTpsyYMeOhfVOzsrJITU3VlrS0tEI9Y3nRuHFjjh07VmB7fHw8M2fOZOrUqY+cJCqKQnh4OHPmzCEuLo6YmBiWL1/O1q1b+fTTT0s6dCFEBRYQEMDJkydxcXGRfppCoE7qfrs899xzdOrUiejoaA4fPszhw4eJioqic+fO9OjRo9j3KnKfzddee43AwMBiB1ASAgMD8ff319ZgglrzuH37dpYtW4a5uTmWlpaEhoby999/F0gc7+fq1atMnz5du4zZbe+++y6///4727ZtIyAggLlz52qbZu5n5syZ2NjYaEtFnXu0UaNGHD9+vMD2CRMm4O7u/tDv093OnDlDWloanTp1olq1alSvXh1/f388PT1JT08vybCFEBXYkiVLWLJkiXY+TScnJ32HJESZ8uWXXzJz5kzs7Oy02+zs7Pj000/58ssvi339IvfZzM7O5ueff2bTpk20bNmSKlWq6OwvSvv+lClT7tn8erewsDBtv0yA6OhoQkNDWblypc5xmZmZDBo0iPbt27N8+XLy8vKYM2cO3bt3JywsDHNz8wfeJzU1lR49euDj48PkyZN19t09d2jjxo21I6lv13bey/jx43X6R8TExDx6wqkokJPxaMeWNGMLKMS6qI0aNeLcuXNkZmZqv8dHjhzhl19+YdOmTY+c6INapW9iYkKjRo0AtXZ46dKlnDlzhoULFxbuOYQQldKJEye0f+ROmzaNJ598Us8RCVH2pKamcuXKFRo0aKCzPSEhoURaYoucbB4/fpzmzZsDarv+3Yq6aHtAQAADBgx44DEeHh46nxcvXoy9vT09e/bU2R4UFMTFixfZu3cvBgYG2m12dnasW7fugfdJS0ujW7duWFpasnbt2nsOJrpb27ZtATh79ux9k01TU1NMTU21n29PD/VIcjLgs+qPfnxJmhALJlUeftwtjRo1Ij8/n4iICO3Px+jRo+nduzdPPfXUPc/x9/fn8OHDpKen4+rqytq1a2nVqhWHDx8mJydHu+RpRkYGTk5OhIaG6vzBIYQQ93K7n2ZmZiZdu3Zl/Pjx+g5JiDKpd+/evPnmm3z55ZfanGbfvn2MGzeOPn36FPv6RU42S2MBdwcHBxwcHB75eEVRWLx4Ma+99lqBhDAjIwMDAwOdxPf25/z8/PteMzU1FX9/f0xNTQkODsbMzOyhcdzuZOvi4vLIsVdUbm5u2NracuzYMZo3b86qVasICwsjIiLivueEhobec/uhQ4fo37+/tn9mYmIi48ePZ/jw4fz777/aPle9e/dm+/btdOnSpcCsAUKIyklRFN555x0iIiKoXr06y5Yt01Y8CCF0zZ8/n/fff59XX32VnJwcFEXB2NiYwYMH88UXXxT7+kVONsuCrVu3cuHCBQYPHlxgX9euXRk3bhwjRoxg5MiR5OfnM2vWLIyMjLTTXcTExNClSxeWLl1K69atSUtLw8/Pj4yMDJYtW6YdzAPg6OiIoaEhe/fuZd++fXTu3BkbGxvCwsIYM2YMPXv2pGbNmqXzoMYWag2jPhhbFPqU2/02s7Ky+OCDD/jggw8KNSjrtiNHjjB9+nTq1KkDQJ06dXj//fd57rnnuHTpErVr1wbUeToHDRrEL7/8Uuh7CCEqpiVLlrB06VLppynEI7CwsGDevHl88cUXnDt3DkVRqFOnToEukkVVrGRz165dLFiwgHPnzrFq1Spq1KjBr7/+Sq1atejQoUOJBPgggYGB+Pr6Ur9+/QL7vL29Wb9+PVOnTqVdu3YYGBjQrFkzQkJCtDWQOTk5REZGkpGh9oc8dOgQ+/fvB9AmOLdduHABDw8PTE1NWbFiBVOnTiUrKwt3d3eGDh3KBx98UHoPqtEUqilb324nm1999RV5eXl8+OGHhb7G+fPnuX79urYp/u7thoaG2qZ1kHk6hRC6jh8/zogRIwCYPn06HTt21HNEQpR9W7ZsYcuWLSQkJBRoAV60aFGxrl3kZHP16tUMHDiQV155hSNHjpCVlQWo/R0/++wz/v7772IF9iiCgoIeuL9r16507dr1vvs9PDxQFEX7uVOnTjqf76V58+bs27evcIFWMo0aNWLlypXs2bOHRYsWPXQw1r0cOnQIjUaDs7Mz8fHxpKens2vXLiZPnsxbb72Fra1tyQcuhCj3bty4oe2n6e/vL6u7CfEIpk6dyrRp02jZsiUuLi5FHntzP0VONj/99FPmz5/Pa6+9xu+//67d7uvry7Rp00okOFE+NWrUiKSkJDp37lzk9c4PHz6srcYHdQqGunXr8tVXX/H666+XZLhCiApCURTefvttTp06RfXq1fn111+ln6YQj2D+/PksWbKEgQMHlsr1i5xsRkZG3rNpwtramuvXrxcnJlHOtW/f/qE1xA8zc+ZMZs6cWUIRCSEqg8WLF2sHAv3++++yypgQjyg7OxtfX99Su36R/+RzcXHh7NmzBbbv3r1bO3BDCCGEeByOHTum7af56aef8sQTT+g5IiHKjyFDhjy0a2JxFLlmc/jw4bz77rssWrQIjUZDbGwse/fu5f333+eTTz4pyRiFeKD7zdMphKgcbty4Qf/+/bl58ybdunUr0qBEISqzmzdv8tNPP7F582YaN25cYDrJoizUc7ciJ5sffPABKSkpdO7cmZs3b9KxY0dMTU15//33CQgIKFZQQhTG/ebpFEJUfDk5OfTv359Tp05Ro0YN7XRHQohHd/ToUZo2bQpQYMnpkhgsVKypj2bMmMHEiRM5efIk+fn5+Pj4YGlpWeyghBBCiIdRFIWhQ4eyYcMGzM3NWbVqlfTTFKIISmOhnrsVe1J3CwsLWTpQCCHEYzd+/Hh++eUXDA0NWblypXaZPSFE2VKsZLM0JwAVQggh7uebb75h9uzZACxcuJBnn31WzxEJUb5dv36dwMBAIiIi0Gg01K9fn8GDB2NjY1Psaxe5Y8vUqVPx8/Njy5YtJCUlkZycrFOEEEKI0vD7778zevRoQO3O9eabb+o3ICHKuYMHD+Lp6cnXX3/NtWvXSEpK4uuvv8bT05PDhw8X+/pFrtks7QlABcWeq1IUJN9TIcq3LVu28NprrwEQEBDA+PHj9RyREOXfmDFj6NmzJwsXLsTISE0Nc3NzGTJkCKNHj2bnzp3Fun6Rk83SngC0Mrs95UBGRkaRlnoU95eRkQFQYFoHIUTZd+TIEXr37k1OTg59+/Zl7ty5Jb6snhCV0cGDB3USTQAjIyM++OCDEhmXU+Rk8/YEoJMmTSp2EEKXoaEhtra2JCQkAOogLPmFWjyKopCRkUFCQgK2trYYGhrqOyQhRCGcP3+eZ555hrS0NDp16sSvv/4q/x8LUUKsra25fPky3t7eOtujoqKwsrIq9vWLnGyW9gSglV21atUAtAmnKBm2trba760QonxISEjA39+fK1eu0KRJE/73v/9hZmam77CEqDBefPFFBg8ezJw5c/D19UWj0bB7927GjRvHSy+9VOzrFznZfNAEoKL4NBoNLi4uODk5kZOTo+9wKgRjY2OpCRGinLlx4wbdu3fn7NmzeHh4sGHDhhIZHSuEuGPOnDloNBpee+01cnNzAfWd+fbbbzNr1qxiX1+jyIgJvYiOjsbNzY2oqChcXV31HY4QQpQ52dnZPPfcc2zcuBEHBwf27NlDvXr19B2WqOQq8vs7IyODc+fOoSgKderUwcLCokSuW+ipj7p3705KSor284wZM7h+/br289WrV/Hx8SmR4IQQQlRO+fn5DBo0iI0bN2JhYcFff/0liaYQpczCwoKGDRvSqFGjEks0oQjJZmhoKFlZWdrPs2fP5tq1a9rPubm5REZGlkx0QgghKqUPPviA3377DSMjI1atWkXr1q31HZIQFVpgYCANGzbEzMwMMzMzGjZsyM8//1wi1y50n83/trpLK7wQQoiS9OWXX/Lll18C6gvwmWee0XNEQlRskyZN4uuvv2bkyJG0a9cOgL179zJmzBguXrzIp59+WqzrF3ttdCGEEKKkLFu2jPfffx9QW85uT+AuhCg9P/74IwsXLtQZed6zZ08aN27MyJEji51sFroZXaPRFJjzUeaAFEIIUVwbN27ULj05evRoxo0bp+eIhKgc8vLy7jl5e4sWLbSj04ujSM3ob7zxBqampoA63+Zbb71FlSpVAHT6cwohhBCP4uDBg/Tp04fc3FwGDBjAl19+KRUZQjwmr776Kj/++GOBOdJ/+uknXnnllWJfv9DJ5uuvv67z+dVXXy1wjDR7CCGEeFRnzpyhe/fupKen06VLF5YsWYKBQaEb3oQQhTB27Fjt1xqNhp9//pmNGzfStm1bAPbt20dUVFSJ5HSFTjYXL15c7JsKIYQQAPHx8fj7+5OYmEizZs1Ys2aNtuVMCFF6jhw5ovO5RYsWAJw7dw4AR0dHHB0dOXHiRLHvJQOEhBBC6EVqairPPPMMFy5coHbt2mzYsAFra2t9hyVEpbBt27bHdi9ppxBCCPHYZWVl0adPH8LDw3FyciI0NBRnZ2d9hyWEKAWSbAohhHis8vPzef3119myZQtVqlTh77//pk6dOvoOSwhRSiTZFEII8dgoisLYsWNZsWIFRkZGrFmzRttXTAhRMUmyKYQQ4rH5/PPP+eabbwBYsmQJfn5+eo5ICFHaJNkUQgjxWPzyyy989NFHgLokZUnM3yeEKPsk2RRCCFHq/v77bwYPHgzA+++/rzPHnxBCP7Zu3YqPjw+pqakF9qWkpNCgQQN27dpV7PtIsimEEKJU7d+/n379+pGXl8err77K7Nmz9R2SEAKYO3cuQ4cOveeUYzY2NgwfPrzAqkJFIcmmEEKIUhMZGUmPHj3IyMjAz8+PwMBAWR1IiDLi33//pVu3bvfd7+fnx6FDh4p9H/k/XgghRKmIjY3F39+fq1ev0rJlS1avXo2JiYm+wxJC3HLlyhWMjY3vu9/IyIjExMRi30eSTSGEECUuNjYWPz8/Ll26RJ06dfjrr7+wtLTUd1hCiLvUqFGDY8eO3Xf/0aNHcXFxKfZ9JNkUQghRos6ePUuHDh04ceIELi4uhIaG4uTkpO+wRCWRm5fPucQbbDgWx42sXH2HU6Z1796dTz75hJs3bxbYl5mZyeTJk3n22WeLfR9ZG10IIUSJCQ8Pp1u3bly5cgVPT082btxI7dq19R2WqIAURSE25San49OIvJLG6fg0TsWncTbxBtm5+QCsHN6O1rWq6jnSsuvjjz9mzZo11KtXj4CAALy8vNBoNERERPDDDz+Ql5fHxIkTi30fSTaFEEKUiJ07d/Lcc8+RmppKkyZNCAkJoVq1avoOS1QASTey7iSVV9KIjE/j9JUb9625NDM2oJ6zFbl5+Y850vLF2dmZf/75h7fffpvx48ejKAoAGo0Gf39/5s2bh7Ozc7HvI8mmEEKIYgsODubFF1/k5s2bPPHEEwQHB2Nra6vvsEQ5k3Yzh9NXbmgTSjWpTONqevY9jzcy0ODpaEm9alZ4OVtSz9kKr2pWuNlZYGCgeczR39vMmTNZs2YNp06dwtzcHF9fX2bPno2Xl9cDz8vKymLatGksW7aM+Ph4XF1dmThxIoMGDdIes3r1aiZNmsS5c+fw9PRkxowZ9O7du1Dxubu78/fff5OcnMzZs2dRFIW6detiZ2dXpOe9F0k2hRBCFMsvv/zC4MGDycvL47nnnmPFihWYm5vrOyxRht3MyeNswq2k8lYT+OkrN4i5nnnP4zUaqFnVQk0mna1uJZdW1HKogolR2R5+smPHDkaMGEGrVq3Izc1l4sSJ+Pn5cfLkSapUqXLf8/r378+VK1cIDAykTp06JCQkkJt7pyZ37969vPjii0yfPp3evXuzdu1a+vfvz+7du2nTpk2h47Szs6NVq1bs2bMHCwuLIj3r/WiU23Wm4rGKjo7Gzc2NqKgoXF1d9R2OEEIUyVdffcV7770HwOuvv87PP/+MkZHUYwhdl66ms+tMEnvPXyUiLpWLSenk3yf7cLY2pZ6zFd7VrLQ1lXWcLLEwKRs/V8V9fycmJuLk5MSOHTvo2LHjPY8JCQlhwIABnD9/nqpV793n9MUXXyQ1NZUNGzZot3Xr1g07OzuWL19e6Lhus7a2Jjw8vET7WpeN/3JCCCHKFUVRmDBhArNmzQJg7NixfPHFFzJhuwAgOT2bf85dZffZRHafTSLqWsEaSxtzY7yq6dZU1nO2xNaifMzFmpaWprPMo6mpKaampg89LyUlBeC+SSSo3VJatmzJ559/zq+//kqVKlXo2bMn06dP17Ya7N27lzFjxuic5+/vz9y5c4vwNHeURh2kJJtCCCEKJS8vj7feeouff/4ZUPukffjhh2g0ZaOPnHj8bubkcehSMrvOJLHnbBLHY1O4O2cxNtTQrKYdHeo40NTNFu9qVjhamZbrnxkfHx+dz5MnT2bKlCkPPEdRFMaOHUuHDh1o2LDhfY87f/48u3fvxszMjLVr15KUlMQ777zDtWvXWLRoEQDx8fEFBu84OzsTHx9ftAcqReUy2dy+fTudO3e+574DBw7QqlUrAMLCwvjoo484dOgQGo2GVq1a8fnnn9O0adP7Xnv48OFs3ryZ2NhYLC0ttR15vb29tcckJyczatQogoODAejZsyffffeddIYXQlR4N2/e5JVXXmHNmjUYGBgwf/58hg4dqu+wxGOWn69wMi6V3WfV5PLAhWtk5eqO/PZytqJ9HQeeqOtA61pVqWJaLlOO+zp58iQ1atTQfn6UWs2AgACOHj3K7t27H3hcfn4+Go2G3377DRsbG0DtstK3b19++OEHbe3mf5N1RVGKncAvWLCgREag361c/pf39fUlLi5OZ9ukSZPYvHkzLVu2BNTqbX9/f3r16sW8efPIzc1l8uTJ+Pv7Ex0dfd/lmVq0aMErr7xCzZo1uXbtGlOmTMHPz48LFy5gaGgIwMsvv0x0dDQhISEADBs2jIEDB7J+/fpSfGohhNCvtLQ0nn/+ebZu3YqJiQlBQUG88MIL+g5LPCbRyRnsPpPE7rNJ/HPuKtf+M0Lc2dpUm1y293TAydpMT5E+HlZWVlhbWz/y8SNHjiQ4OJidO3c+tK+ni4sLNWrU0CaaAPXr10dRFKKjo6lbty7VqlUrUIuZkJBQqERx69atBAQEsG/fPu2zvPzyy4Da3O/r68v8+fN54oknHvma91Iuk00TExOdudtycnIIDg4mICBAm9FHRkaSnJzMtGnTcHNzA9Qq7saNG3P58mU8PT3vee1hw4Zpv/bw8ODTTz+lSZMmXLx4EU9PTyIiIggJCWHfvn3a0V4LFy6kXbt2REZGPnQqAyGEKI8SExPp3r07Bw8exNLSknXr1vHUU0/pOyxRilIyc9h7q9/lnrNXuZCUrrO/iokhbWvb06GuAx3qOFDHybJcN4uXFkVRGDlyJGvXrmX79u3UqlXroee0b9+eP/74gxs3bmiXeT19+jQGBgbaRLVdu3Zs2rRJp9/mxo0b8fX1feTY5s6dy9ChQ++ZNNvY2DB8+HC++uqrypls/ldwcDBJSUm88cYb2m1eXl44ODgQGBjIhAkTyMvLIzAwkAYNGuDu7v5I101PT2fx4sXUqlVLm7Du3bsXGxsbnWkF2rZti42NDf/88899k82srCyysrK0n9PS0orwpEII8fhdvnwZPz8/IiMjcXBwYMOGDdpWJFFxZOfmc/hysrb28mj0dZ0R44YGGpq62WprL5u62WJsKAPCHmbEiBEEBQWxbt06rKystLWRNjY22ubw8ePHExMTw9KlSwG1dnH69Om8+eabTJ06laSkJMaNG8egQYO057z77rt07NiR2bNn06tXL9atW8fmzZsf2kR/t3///ZfZs2ffd7+fnx9z5swp6qNrVYhkMzAwEH9/f21CCGr19vbt2+nVqxfTp08HoF69eoSGhj50Wo558+bxwQcfkJ6ejre3N5s2bcLERB0dFx8ff881fp2cnB7YKXfmzJlMnTq1KI8nhBB6c/LkSfz8/IiJicHNzY1NmzZJC04FkpB2k5Dj8Ww9lcD+89fIzMnT2e/pWIUOdRzoUNeRNrWrYm127y5o4v5+/PFHADp16qSzffHixdpKsri4OC5fvqzdZ2lpyaZNmxg5ciQtW7bE3t6e/v378+mnn2qP8fX15ffff+fjjz9m0qRJeHp6smLFikLNsXnlypX7disEMDIyIjEx8ZGvd9/rFPsKJWjKlCkPTcjCwsJ0/qKOjo4mNDSUlStX6hyXmZnJoEGDaN++PcuXLycvL485c+bQvXt3wsLCHjjh8CuvvELXrl2Ji4tjzpw59O/fnz179mBmpvY/uVczwcM65Y4fP56xY8dqP8fExBQYySaEEGXJ/v376d69O9euXaN+/fqEhobq/FEvyqfEtCxCTsTz19FY9l+4pjNq3MHShPZ11Gbx9nUcqG4rk/MX16NMJbRkyZIC225Xdj1I37596du3b1FDo0aNGhw7dow6dercc//Ro0dxcXEp8vVvK1PJZkBAAAMGDHjgMR4eHjqfFy9ejL29PT179tTZHhQUxMWLF9m7d6923regoCDs7OxYt27dA+9jY2ODjY0NdevWpW3bttjZ2bF27VpeeuklqlWrxpUrVwqck5iY+MBOuf+df+vuubmEEKKs2bRpE7179yY9PZ3WrVvz999/Y29vr++wRBFdvXE7wYxj3/mrOs3jTd1seaZhNTrWc8TL2arMLPMoSl/37t355JNPeOaZZ7QVardlZmYyefJknn322WLfp0wlmw4ODjg4ODzy8YqisHjxYl577bUC1cAZGRkYGBjo1Dbe/pyfn//fSz30Prf7W7Zr146UlBQOHDhA69atAfWv/9ujtoQQorxbuXIlr776Kjk5OXTt2pU1a9ZoBymI8uNaejahtxLMveevkndXhtnE1YYejV14pqELblVLdmlCUX58/PHHrFmzhnr16hEQEICXlxcajYaIiAh++OEH8vLymDhxYrHvU6aSzcLaunUrFy5cYPDgwQX2de3alXHjxjFixAhGjhxJfn4+s2bNwsjISDtHZ0xMDF26dGHp0qW0bt2a8+fPs2LFCvz8/HB0dCQmJobZs2djbm5O9+7dAXXqgW7dujF06FAWLFgAqCPYn332WenHJIQo93788UdGjBiBoij079+fpUuXPtL8gaJsSE7PZuPJeP48Gsc/53QTzEY11ASzRyNJMIXK2dmZf/75h7fffpvx48drm/w1Gg3+/v7MmzevRObcLNfJZmBgIL6+vtSvX7/APm9vb9avX8/UqVNp164dBgYGNGvWjJCQEG3/g5ycHCIjI8nIyADAzMyMXbt2MXfuXJKTk3F2dqZjx478888/OoOCfvvtN0aNGoWfnx+gTur+/fffP4YnFkKI0qEoCjNmzGDSpEkAvPXWW3z//ffa+YVF2ZWSkUPoSbUGc8/ZJHLvSjAbVLfWJpju9lX0GKUoq9zd3fn7779JTk7m7NmzKIpC3bp1sbOzK7F7aJTSWARTPFR0dDRubm5ERUU9dHJXIYQoTfn5+YwdO5ZvvvkGUBfJmDp1qsyZWIalZOSw8WQ8fx+LY/fZJHLy7rzKfVzuJJgeDpJgljR5fxdeua7ZFEIIUTw5OTkMGjSIZcuWAfDNN98watQoPUcl7iX1Zg6bTlzhr2Nx7DqTqJNgelez4tnGLnRv5EJtR+lfK8oWSTaFEKKSysjIoF+/fvz9998YGRmxZMkSXnnlFX2HJe6SdjOHzRFX+OtoHDtPJ5Gdd2eAq5ezFT1uJZh1nCTBFGWXJJtCCFEJJScn89xzz7Fnzx7Mzc1ZtWqVdiCk0K/s3Hy2nrrC6sMx7DidSHbunQSzjpMlz95qIq/rbKXHKIV4dJJsCiFEJRMXF4e/vz/Hjh3D1taWP//8k/bt2+s7rEpNURROxKay6lA068JjSM7I0e6r7ViFZxtX59nGLtSTBFOUQ5JsCiFEJXLq1Cm6d+/OhQsXcHFxISQkhMaNG+s7rEorMS2LdeExrDoUzan4NO12JytT+jR3pVfT6nhXs5LBWqJck2RTCCEqieXLlzN06FDS09Px9PRk48aN1K5dW99hVTq3m8lXHYpmW2Sidi5MEyMD/Hyc6dvClQ51HDAyNNBzpEKUDEk2hRCigsvKymLMmDH8+OOPAHTq1Inff/+9RCZrFo/mQc3kTd1s6dvClecaV8fGwvgBVxGifJJkUwghKrALFy7Qv39/Dh48CMDEiROZMmUKRkby6/9xeFgzed8WNajjJP0wRcUmv22EEKKCWr9+Pa+99hrXr1+natWqLFu2jGeeeUbfYVV40kwuhC5JNoUQooLJzc1l4sSJfP755wC0adOGlStXUrNmTT1HVnFJM7kQ9yfJphBCVCCxsbEMGDCAXbt2AfDuu+/y+eefY2JioufIKiZpJhfi4STZFEKICmLLli28/PLLJCQkYGVlxaJFi+jbt6++w6pwpJlciMKRZFMIIcq5/Px8ZsyYweTJk1EUhcaNG7Nq1Srq1q2r79AqlIi4VFaERUkzuRCFJMmmEEKUY0lJSbz66quEhoYCMHjwYL777jvMzc31HFnFkHYzh/X/xrEi7DL/Rqdot0szuRCPTpJNIYQop/bu3Uv//v2Jjo7G3NycH3/8kddff13fYZV7iqJw+HIyvx+I4s+jcWTm5AFgZKChq48z/Vu58YQ0kwvxyCTZFEKIckZRFObOncsHH3xAbm4u9erVY9WqVTRq1EjfoZVrV29ksfZIDL+HRXE24YZ2u6djFV5s5Uaf5q44WJrqMUIhyidJNoUQohxJSUlh0KBBrFmzBoAXX3yRhQsXYmUlTblFkZ+vsPtsEivCoth4Mp6cPHWwj5mxAc82rs6AVm60cLeTtcmFKAZJNoUQopwIDw+nb9++nDt3DmNjY77++mveeecdSYSKIOZ6Jn8cjOKPg9HEXM/Ubm/sasOLrdx4rkl1rM1ksI8QJUGSTSGEKOMUReHnn39m5MiRZGVl4e7uzsqVK2ndurW+QytXsnPz2RJxhd/Doth5JhFFrcTE2syI3s1q0L+VGw2q2+g3SCEqIEk2hRCiDEtPT+ftt9/m119/BaBHjx4sXbqUqlWr6jmy8uNswg1WHoxi9aForqZna7e3q23PgNZu+DeohpmxoR4jFKJik2RTCCHKqIiICPr168eJEycwNDRkxowZjBs3DgMDGQX9MBnZufx1NI4VYVEcvJSs3e5oZUq/Fq70b+mGh0MVPUYoROUhyaYQQpRBy5cvZ+jQoaSnp+Pi4sLvv/9Ox44d9R1WmaYoCsdiUvg9LIrg8FhuZOUCYGigobOXEy+2cqOzl6NMWSTEYybJphBClCFZWVmMGTOGH3/8EYCnnnqKoKAgnJ2d9RxZ2ZWSkcP/wtUpiyLiUrXb3e0t6N/Sjb4tXHG2NtNjhEJUbpJsCiFEGXHhwgX69evHoUOHAPj444+ZMmUKhobSn/C/8vMV9p6/yoqwKEJOxJOdmw+o65M/07AaL7Zyo20tewwMZKS+EPomyaYQQpQB//vf/3jzzTe5fv069vb2LFu2jG7duuk7rDInLiWTVQejWXkoiqhrd6Ys8q5mxUuta9KraXVsLUz0GKEQ4r8k2RRCCD2Kj4/n3XffZeXKlQC0bduWlStX4ubmpufIyo7bUxatOBjFztOJ5N+assjK1IieTaszoFVNGtawlvlGhSijJNkUQgg9yM/PJzAwkA8++IDr169jaGjIe++9x/Tp0zExkZo5gLMJaawIi2LN4RidKYva1KrKi63ceKahC+Ym0sVAiLJOkk0hhHjMIiIiGD58OLt27QKgRYsWLFy4kGbNmuk5Mv1Lz8rlz6OxrAiL4vDl69rtTlamvHBryqJaMmWREOWKJJtCCPGYZGVlMXPmTGbOnEl2djZVqlTh008/JSAgACOjyvvrWFEUDl++zoqwy/x5NI6M7DxAnbLoKW8nXmzpRieZskiIcqvy/nYTQojHaNeuXQwbNoxTp04B6kpAP/zwA+7u7nqOTH+SbmSx9nAMKw5GcTbhhnZ7LYcq9G/pxgstauBkJVMWCVHeSbIphBClKDk5mQ8//JCFCxcC4OzszLfffku/fv0q5YCWvHyFnWcSWXEgis0RV8i9NdrHzNiAHo2q82IrN1p52FXK740QFZUkm0IIUQoUReGPP/5g1KhRXLlyBYBhw4Yxa9Ys7Ozs9Bzd4xd1LYOVB6NYdSiauJSb2u1NXG3o38qN55pUx9rMWI8RCiFKiySbQghRwi5dusSIESP466+/APD29uann37iiSee0HNkj9fNnDxCT8SzIiyKf85d1W63tTCmd7MavNjKDe9q1nqMUAjxOEiyKYQQJSQvL4/vvvuOjz/+mPT0dExMTJgwYQIfffQRpqam+g7vsVAUhSNR11l1KJo//40l9aa6PrlGAx3qONC/pRt+DZwxNZIpi4SoLCTZFEKIEnDkyBGGDh2qXWryiSeeYMGCBdSvX1/PkT0esdczWXskhtWHojmflK7dXt3GjH4t3ejX0hVXOws9RiiE0BdJNoUQohjS09OZOnUqX331FXl5edjY2PDFF18wePBgDAwq9lQ9mdl5hJyIY/WhGPacS0K5tbKPubEhzzSsxgstXGlXW9YnF6Kyk2RTCCGKKDQ0lLfeeouLFy8C8OKLLzJ37lyqVaum38BKkaIoHLhwjdWHo/n7WDw3snK1+9rUqsoLLVzp3sgFS1N5vQghVPLbQAghCikhIYExY8YQFBQEQM2aNZk3bx49evTQc2SlJ+paBqsPR7PmcAyXr2Vot7tVNeeF5q70aeZKTXtpJhdCFCTJphBCPCJFUViyZAnvvfceycnJGBgY8O677zJt2jQsLS31HV6Ju5GVy9/H4lh9KJr9F65pt1cxMaRHYxdeaO5KK4+q0kwuhHggSTaFEOIRnD59muHDh7N9+3YAmjVrxsKFC2nRooV+Ayth+fkK+85fZdWhaDYcjyczR106UqOB9p4OvNCiBv4NqmFhIq8PIcSjkd8WQgjxANnZ2XzxxRdMnz6drKwszM3NmTZtGqNHj65Q65lfSEpn9aFo1h6JIeZ6pnZ7bYcqvNDCld7NalDd1lyPEQohyquK85tSCCFKkKIobN68mTFjxnDixAkA/P39+fHHH6lVq5aeoysZKZk5/HU0jtWHozl0KVm73crMiOeaVOeF5q40r2krS0cKIYqlXM7LsX37djQazT1LWFiY9riwsDC6dOmCra0tdnZ2+Pn5ER4e/sBrDx8+HE9PT8zNzXF0dKRXr16cOnVK5xgPD48C9/3oo49K41GFEHqwd+9ennrqKfz8/Dhx4gSOjo4EBQWxYcOGcp9o5ublsz0ygZHLj9B6xmYmrD3GoUvJGGigk5cj373UjLCJT/NZ70a0cJc1yoUQxVcuazZ9fX2Ji4vT2TZp0iQ2b95My5YtAUhLS8Pf359evXoxb948cnNzmTx5Mv7+/kRHR2NsfO81eFu0aMErr7xCzZo1uXbtGlOmTMHPz48LFy5gaHhnxYtp06YxdOhQ7eeKODhAiMrm6NGjTJw4kT///BMAU1NT3nnnHSZOnIi9vb2eoys6RVE4fPk6weEx/HUsjqQb2dp99ZwteaG5K883q4GztZkeoxRCVFTlMtk0MTHRmccuJyeH4OBgAgICtH+FR0ZGkpyczLRp03BzcwNg8uTJNG7cmMuXL+Pp6XnPaw8bNkz7tYeHB59++ilNmjTh4sWLOudYWVkVai69rKwssrKytJ/T0tIe+VwhROk6c+YMkydPZvny5QAYGhoyaNAgJk2apP39UR6dik9lXXgs6/+NJTr5Tj9MOwtjejapzgstXGlUw0ZqL4UQpapcJpv/FRwcTFJSEm+88YZ2m5eXFw4ODgQGBjJhwgTy8vIIDAykQYMGuLu7P9J109PTWbx4MbVq1Srwwpk9ezbTp0/Hzc2Nfv36MW7cOExMTO57rZkzZzJ16tQiPZ8QonRER0czbdo0Fi1aRF6eOup6wIABTJ06lXr16uk5uqKJupZB8L+xrAuP4fSVG9rtVUwM8WtQjZ5Nq9OhjgPGhuWyF5UQohzSKMrtBcbKr+7duwPw999/62w/ceIEvXr14sKFCwDUq1eP0NBQatas+cDrzZs3jw8++ID09HS8vb35888/dWo1v/76a5o3b46dnR0HDhxg/Pjx9OrVi59//vm+1/xvzWZMTAw+Pj5ERUXh6upa6GcWQhRdYmIis2bN4ocfftD+f9mjRw8+/fRTmjZtqt/giiAh7SZ/HY0j+N9Yjly+rt1uYmhAJy9HejatThdvZ8xNDO9/ESHEI4mOjsbNzU3e34VQppLNKVOmPLT2LywsTNsvE9T/6O7u7qxcuZIXXnhBuz0zM5NOnTrh7e1NQEAAeXl5zJkzh1OnThEWFoa5+f2n8EhJSSEhIYG4uDjmzJlDTEwMe/bswczs3v2ZVq9eTd++fUlKSnrkfl3ywyrE45eSksJXX33FV199xY0baq1fx44d+eyzz2jfvr2eoyuclMwcQk/EExweyz/nksi/9ZvcQAPtPO3p1aQG/g2rYWN+7/7pQoiikfd34ZWpZvSAgAAGDBjwwGM8PDx0Pi9evBh7e3t69uypsz0oKIiLFy+yd+9eDAwMtNvs7OxYt27dA+9jY2ODjY0NdevWpW3bttjZ2bF27Vpeeumlex7ftm1bAM6ePVuuBxEIUVFlZmby/fffM2vWLK5dU1fCadGiBZ999hldu3YtN30Wb+bksSUigeB/Y9h2KpHsvHztvqZutvRsUp1nG7vgJAN9hBBlSJlKNh0cHHBwcHjk4xVFYfHixbz22msFRpdnZGRgYGCg8xK5/Tk/P/+/l3rofe5uAv+vI0eOAODi4lKo6wohSld2djaLFi1i+vTpxMbGAuDt7c2nn35Knz59ykWSmZOXz+6zSawPjyX0RDzp2XnafXWdLOnVtDrPNamOu30VPUYphBD3V6aSzcLaunUrFy5cYPDgwQX2de3alXHjxjFixAhGjhxJfn4+s2bNwsjIiM6dOwNqv8kuXbqwdOlSWrduzfnz51mxYgV+fn44OjoSExPD7NmzMTc31/YL3bt3L/v27aNz587Y2NgQFhbGmDFj6Nmz50P7ggohHo+8vDyWL1/O5MmTOX/+PADu7u5MnTqVV199VWcas7IoP1/h0OVk1oXH8PexeK6l35mqqIatOT2bVqdnk+p4V7MqFwmzEKJyK9fJZmBgIL6+vtSvX7/APm9vb9avX8/UqVNp164dBgYGNGvWjJCQEG0NZE5ODpGRkWRkZABgZmbGrl27mDt3LsnJyTg7O9OxY0f++ecfnJycAHXevRUrVjB16lSysrJwd3dn6NChfPDBB4/vwYUQ96QoCsHBwXz88cccP34cAGdnZyZNmsSQIUMwNTXVc4T3pygKJ+NSCf43lvXhscSm3NTus69iQo/GLvRqWp1mbnYYGEiCKYQoP8rUAKHKRDoYC1GytmzZwoQJEzhw4AAAtra2fPjhh4wcOZIqVcpmE7OiKByPSSXkRBwbjsdzPjFdu8/S1Aj/W1MVtfe0x0imKhKiTCjs+3vmzJmsWbOGU6dOYW5ujq+vL7Nnz8bLy+u+52zfvl3bCnu3iIgIvL29AViyZAlvvvlmgWMyMzPvO6BZX8p1zaYQQuzbt4+JEyeydetWAKpUqcLo0aN5//33sbW11W9w95CXr3DoUjIhx+MJPRFPzPU7k62bGBnwlJcTvZpWp7O3E2bGZbu5XwjxcDt27GDEiBG0atWK3NxcJk6ciJ+fHydPnnzoH8KRkZFYW1trPzs6Ourst7a2JjIyUmdbWUs0QZJNIUQ5dezYMSZNmsS6desAdWWxt99+m/Hjx+Ps7Kzn6HRl5+az7/xVNhyPZ9PJeJ3lIs2NDens7Yh/g2p09nbC2kymKhKiIgkJCdH5vHjxYpycnDh06BAdO3Z84LlOTk4P/KNZo9EUajVDfZFkUwhRruzfv5+vv/6alStXoigKBgYGvPHGG3zyySePvDrY45CZncfOM4mEHo9nc8QVUm/mavdZmRnRtb4z/g2r8WQ9R6nBFKIcSktLIzU1VfvZ1NT0kfqFp6SkAFC1atWHHtusWTNu3ryJj48PH3/8cYGm9Rs3buDu7k5eXh5NmzZl+vTpNGvWrJBPUvok2RRClHnZ2dmsWrWKb7/9lv3792u39+/fn6lTp2r7MOlb2s0ctp5KIOR4PNsjE8nMuTNNkYOlCX4NqtGtQTXa1rbHxEj6YApRnvn4+Oh8njx5MlOmTHngOYqiMHbsWDp06EDDhg3ve5yLiws//fQTLVq0ICsri19//ZUuXbqwfft2bW2ot7c3S5YsoVGjRqSmpvLNN9/Qvn17/v33X+rWrVvs5ytJMkBIT2SAkBAPl5CQwIIFC/jxxx+Ji4sD1Obyl156iTFjxtCkSRM9RwhXb2SxOeIKIcfj2XP2qs5E6zVszfFvUI1uDavRwt0OQxlFLkS5d/v9ffLkSWrUqKHd/ig1myNGjOCvv/5i9+7dhX73P/fcc2g0GoKDg++5Pz8/n+bNm9OxY0e+/fbbQl27tEnNphCizDl8+DDffvsty5cvJztb7d/o4uLC22+/zfDhw7VTkelLXEomocfjCTkRz4EL17RLRQLUdqzCMw2r0a2BCw1rWMs8mEJUUFZWVjqDdx5m5MiRBAcHs3PnziJVMrVt25Zly5bdd7+BgQGtWrXizJkzhb52aZNkUwhRJuTm5rJ27Vq+/fZbdu/erd3epk0bRo0aRd++fTExMdFbfBeT0gk5EU/I8XjCo67r7GtQ3Zput2ow6zpb6SdAIUSZpCgKI0eOZO3atWzfvp1atWoV6TpHjhx54EqFiqIQHh5Oo0aNihpqqZFkUwihV1evXmXhwoXMmzePqKgoAIyMjOjfvz+jRo2iTZs2eokrP1/hRGwqmyOuEHoinlPxaTr7W7jb8UzDavg3qIZbVQu9xCiEKPtGjBhBUFAQ69atw8rKivj4eABsbGwwNzcHYPz48cTExLB06VIA5s6di4eHBw0aNCA7O5tly5axevVqVq9erb3u1KlTadu2LXXr1iU1NZVvv/2W8PBwfvjhh8f/kA8hyaYQQi+OHTvGt99+y7Jly7h5U10tx9HRkbfeeou33nqL6tWrP/aYUjJz2H0miW2RCWyPTCTpRpZ2n6GBhna17fFvWA1/H2ecrMveXHZCiLLnxx9/BKBTp0462xcvXswbb7wBQFxcHJcvX9buy87O5v333ycmJgZzc3MaNGjAX3/9pV06G+D69esMGzaM+Ph4bGxsaNasGTt37qR169al/kyFJQOE9EQGCInKKC8vjz///JNvvvmGbdu2abc3a9aMd999lxdffPGxTkisKAqRV9LYdiqRbZEJHLqUTN5dHTAtTAxpX8cBPx9nnq7vjF0V/TXjCyHKBnl/F57UbAohSt3169dZtGgR33//PRcuXADA0NCQPn36MGrUKNq3b//YBtKkZ+Wy52wS2yIT2R6ZQNxda5ADeDpWobOXE529nWjpYYepkcyBKYQQxSHJphCi1Jw6dYpvv/2WX375hYyMDECdyHjYsGG8/fbb1KxZs9RjUBSF80npbDulNo0fuHBNZ3oiUyMDfD3t6eztRKd6TtS0l/6XQghRkiTZFEKUqPz8fEJCQvjmm2/YuHGjdnvDhg159913efnll7GwKN2E7mZOHnvPX2X7qQS2RSZy+VqGzv6aVS3o7OVIJ28n2tW2lxV8hBCiFEmyKYQoEcnJySxbtozvvvtOO8+bRqOhV69ejBo1ik6dOpVqU3nUtQy2RSaw7VQC/5y7SlbundpLE0MD2tSuSicvJzp5OVLboYrMfymEEI+JJJtCiCJLT09n/fr1LF++nA0bNpCTkwOoU3oMGTKEESNGFHlOuYfJzs3n4MVrbD2VwLbIBM4lpuvsd7Exo7O3E529nPD1tKeKqfy6E0IIfZDfvkKIQsnJyWHjxo3aeePS0+8keY0bN2b48OG89tprWFpaluh98/PVkeP/nLvK3nNX2XsuifTsO2uPGxpoaOlup00w6zlbSu2lEEKUAZJsCiEeKj8/n927dxMUFMSqVau4evWqdl/t2rV56aWXeOmll2jQoEGJ3VNRFM4m3GDveTW53Hf+KskZOTrHOFia0tnLkc7eTrSv44CNuXGJ3V8IIUTJkGRTCHFPt5c+CwoK4vfffyc6Olq7z9nZmRdffJGXX36Z1q1bl0gNoqIoXLyaodZa3kow755UHdR5L1t5VKWdpz3tPR1oUN0aAwOpvRRCiLJMkk0hhI4zZ86wfPlygoKCiIyM1G63sbGhT58+vPzyy3Tq1Akjo+L/+oi6lqFNLPeeu0p8qu6cl6ZGBrT0sKNdbXvaeTrQ2NUGY0ODYt9XCCHE4yPJphCC2NhYVqxYQVBQEAcPHtRuNzMz47nnnuOll17imWeeKfbqPnEpmdrEcu/5q0QnZ+rsNzE0oGlNW9rVtsfX056mNW1lUnUhhCjnJNkUopJKTk5m9erVBAUFsX37dm6vXGtoaEjXrl156aWXeP7557G2ti7yPRLTsnT6XF5I0h0xbmSgobGrDb6eDrTztKd5TTvMTSS5FEKIikSSTSEqkYyMDNavX09QUJDOVEUA7du35+WXX6Zv3744OTkV6frX0rPZf/5On8szCTd09htooFENG9p62tOutj2tPKrKlERCCFHByW95ISq421MVLV++nP/9738Fpip6+eWXGTBgAO7u7oW6bl6+wukraRy+nMyRy9c5cjm5wFyXGg3Ur2ZNu1vJZevaVbE2kxHjQghRmUiyKUQFdOnSJUJDQwkNDWXLli2kpKRo99WqVYuXX3650FMVJd3I0iaVRy5f59/o62TcNc/lbfWcLW8N6LGnTS177KqYlMgzCSGEKJ8k2RSiAkhPT2fHjh3aBPPuUeRQ+KmKsnPziYhL5cjlZA5fvs6RqGSirmUWOM7S1IgmbjY0r2lHs5q2NHWzo6okl0IIIe4iyaYQ5ZCiKBw7dkybXO7atYvs7GztfkNDQ9q2bYu/vz/+/v60aNECQ8P7D7yJvZ55p9Yy6jrHYlLIvmttcVCbxOs6WdLMTU0sm9W0o46TJYYyz6UQQogHkGRTiHIiMTGRTZs2ERoaysaNG4mPj9fZ7+7urk0un3rqKWxtbe95nZs5eRyLSdE2hx+5fL3A/JYAthbGNHNTk8rmNe1o7GYj/S2FEEIUmiSbQpRROTk57N27V1t7efjwYe30RAAWFhZ06tRJm2DWq1evQPN4fr7CxavpHI1O0Q7kiYhLJTdf0TnO0ECDdzUrbXN4s5p2eNhbyNriQgghik2STSHKkPPnz2uTy61bt5KWlqazv0mTJtrksn379piammr33cjK5VRcKhFxqZyMSyMiLpXI+DQycwoO4nG0MqX5raSymZstjVxtsDCRXwdCCCFKnrxdhNCjGzdusG3bNm2CefbsWZ39Dg4O+Pn54e/vT9euXXFxcUFRFKKuZbLtTDIRcamcik8lIi6Ny9cy7nkPUyMD6rtY07ymHc3d1QSzuo2Z1FoKIYR4LCTZFOIxio+P58CBAxw4cIA9e/awZ88enYnVjYyM8PX11dZeejVoxJmEdCLi0pi3/yoRcRc4FZ/Gjazce16/mrUZ9V2sqO9irS21HKrIIB4hhBB6I8mmEKUkNTWVQ4cOaZPLsLAwoqKiChxXu3Zt/Pz8aPVEZxzqNudymkJEXBrjtqVyYdUmFKXgtU0MDajjZHkrobTCx8UabxdrmXZICCFEmSPJphAlIDs7m6NHj2oTywMHDnDq1CmdAT0AGo0G7/o+1G3YFKfaDTD3aMoVbPknPo0NR3PgaGSBaztYmtxVU6nWWno6WmJsaPC4Hk8IIYQoMkk2hSik/Px8Tp8+ra2tPHDgAOHh4TrzXN7mVN2VGnUbYenmjeLgyXULV9LyjPgXQAEu5APXAHVEeB1HywLN4I5WpgWuK4QQQpQXkmwK8RAxMTE6TeFhYWGkpqYWOM7c0gZb9/oYONclu2ptTKvVxbCKLde4nU4CtwaGu9iYUcuhCt7V7tRW1nW2xNTo/hOvCyGEEOWRJJtC3OXatWscPnxYJ7mMjY0tcJyBsSkmzp4YV6uLqUs9TKp7YWTjrB3hbQTYmBtT27EKtRyqUNuhCrUcLKnlUAUPBwuZZkgIIUSlIW88Uemkp6dz9uxZTp8+zenTp4k4FUlE5GnOnT1DSvK1gidoDDB2qKkmlS71MHWph7FDTTSGRpgZG+BhX0WbVN5OKGs7VMFOBusIIYQQkmyKiik7O5sLFy4QGRlJ+PEIjp6M4OyZs1y+cI7kxPgHnmto44zpraTSxKUu5i51cXe2u5NMOt6uqaxCNWszDGRaISGEEOK+JNkU5VZ+fj6nz11gX/gJ/j1+klOnTnPxwjniL18gJTEWJb/gyjm3GZhZYVS1OsZVa2BsVwNzhxrUcPektqcntVzsb9VUqrWUNataYGIkI7+FEEKIopBkU5RJefkKyelZRCVc4+ylWM5evEzEqUjOnT1DzOULXI25RHpSNEpuzn2voTE2w8hOTShtq9WkWk0PPGrXxbteXerUdKG6rTnVbcypbmtG1SomsqKOEEIIUQok2RSPhaIopN7MJSElg3NRcVyIiuNybBwxsfEkJCSQmJhA8tVE0pKvkpF6jezUZPIyUlBysx58YQMjTKtWx9rZDSc3D2rW8qRu3Xo0buCFj6cHrnYWOFubSc2kEEIIoSflMtncvn07nTt3vue+AwcO0KpVKwDCwsL46KOPOHToEBqNhlatWvH555/TtGnTh95DURS6d+9OSEgIa9eu5fnnn9fuS05OZtSoUQQHBwPQs2dPvvvuO2xtbYv7aOWCoijk5CnczM0jJSOHqMRkzl+O5VJ0HNFxccTFXSExMZFrVxNJuZpEespVbqYlk5t+nfyMVNQJJh+dgbEZ5jZVsa9eE1d3Tzzr1sXHux4tGvvQokE9qlqal86DCiGEEKLYymWy6evrS1xcnM62SZMmsXnzZlq2bAlAWloa/v7+9OrVi3nz5pGbm8vkyZPx9/cnOjoaY2PjB95j7ty5921Wffnll4mOjiYkJASAYcOGMXDgQNavX18CT1c8Zy/FcuXadVLTM0lLz+BGRiY30jNJz8jkRuZNMjIzyci8SWZmJpk3s8jMvMnNm1nczL5JVmYWWdlZZGfdJDs7m5zsrFslm9zsbPJys8nLUf9VcnNQcnPIv5mGknOzcEFqNJhUscHCpirWdvbY2Tvg4OhMNWcnarhUw62GC7XdXKjj7koNl2pUqVKldL5ZQgghhCh15TLZNDExoVq1atrPOTk5BAcHExAQoE0QIyMjSU5OZtq0abi5uQEwefJkGjduzOXLl/H09Lzv9f/991+++uorwsLCcHFx0dkXERFBSEgI+/bto02bNgAsXLiQdu3aERkZiZeXV0k/bqF06NaLK6cOPvb7GhibYG5tj5VtVWyqOuDg6IiTkzMu1Zxwq14dDzcXPN2qU9O1Ovb29hgZlcsfPSGEEEIUUoV44wcHB5OUlMQbb7yh3ebl5YWDgwOBgYFMmDCBvLw8AgMDadCgAe7u7ve9VkZGBi+99BLff/+9TkJ72969e7GxsdEmmgBt27bFxsaGf/75577JZlZWFllZd/ofpqWlFeFJH87U3AKNsSkGRsYYGJliaGSMobEJhsYmGJuYYmRsgrGJCSamphibmGJqaoqJqSlmpmaYman/mpubYW5mioW5OebmZlQxN6eKhRmWFubaYmVpgZWFGU4O9jg7O2NpaSkDbIQQQghRQIVINgMDA/H399fWYML/27v3qKrK9A/g38PtcO1oMgIqIoYICuRRwAElckJw8IJaaYqI17EWBmjjFS2slLyhM5aYVpiGl2YWM6JTIQ2ISCHIRQgQb+BdsFEBgTTg/f3Rz706HUQIjgfh+1nrrOV+97v3++z9LOTh3ZcDmJmZ4dixYwgICMB7770HALC3t0diYmKzs2qLFi2Cp6cnAgICmlx/8+ZN9OzZU629Z8+euHnz0e9vjIqKwpo1a1p6SL/bpZxUjY9BRERE1FId6hHdyMhIyGSyZj+nTqleIr569SoSExMxd+5clfa6ujrMmTMHI0aMQEZGBtLT0zF48GD4+/ujrq6uyfETEhKQnJyMrVu3NhtnUzN4QohmZ/ZWrFiByspK6VNUVNTsGERERESdQYea2Vy4cCFee+21Zvv069dPZTk2NhY9evTAhAkTVNr37duHsrIyfP/999DR0ZHaunfvjkOHDjU5TnJyMi5cuKD2VPnLL78MLy8vHDt2DJaWligvL1fb9tatW7CwsHhk3HL5L5esH6qqqmr2OImIiIg6gw5VbJqbm8Pc3LzF/YUQiI2NxcyZM9WeLq+trYWOjo7KbOPD5cbGxib3t3z5csybN0+lzdnZGVu2bMH48eMBAB4eHqisrERmZibc3d0BACdPnkRlZSU8PT1bHDsRERFRV9ChLqO3VnJyMkpLS9UuoQPA6NGjcefOHYSEhKC4uBiFhYWYPXs29PT0pHd0Xrt2DQ4ODsjMzAQAWFpawsnJSeUDAH379oWtrS0AwNHREWPGjMH8+fORkZGBjIwMzJ8/H+PGjdP6k+hEREREHc1TXWx++umn8PT0hKOjo9o6BwcHHD58GPn5+fDw8ICXlxeuX7+Ob775Rnqd0c8//4ySkhLU1ta2aty4uDg4OzvD19cXvr6+cHFxwd69e9vlmIiIiIg6E5kQonVf50Lt4urVq7C2tsaVK1fQp08fbYdDRERELcDf3633VM9sEhEREVHHxmKTiIiIiDSGxSYRERERaQyLTSIiIiLSGBabRERERKQxLDaJiIiISGNYbBIRERGRxrDYJCIiIiKN6VDfjd6VPPx+9hs3bmg5EiIiImqph7+3H/4ep8djsakl5eXlAAB3d3ctR0JEREStVV5ejr59+2o7jKcCv65SS+rr65GbmwsLCwvo6LTf3QzV1dUYNGgQioqKYGZm1m77pfbFPD09mKunA/P0dOgMeWpsbER5eTmUSiX09Dhn1xIsNjuZqqoqKBQKVFZW4plnntF2OPQIzNPTg7l6OjBPTwfmqWviA0JEREREpDEsNomIiIhIY1hsdjJyuRzvvPMO5HK5tkOhZjBPTw/m6unAPD0dmKeuifdsEhEREZHGcGaTiIiIiDSGxSYRERERaQyLTSIiIiLSGBabRERERKQxLDY7me3bt8PW1haGhoYYNmwY0tLStB1SlxYVFQU3NzeYmZmhZ8+emDhxIkpKSlT6CCEQGRmJXr16wcjICC+++CIKCwu1FDFFRUVBJpMhPDxcamOOOo5r165hxowZ6NGjB4yNjTFkyBBkZ2dL65kr7auvr8eqVatga2sLIyMj9O/fH++++67Kd4kzT10Li81O5ODBgwgPD0dERARyc3Ph5eWFP//5z7h8+bK2Q+uyUlNTERISgoyMDCQlJaG+vh6+vr6oqamR+mzYsAHR0dH48MMPkZWVBUtLS4wePRrV1dVajLxrysrKws6dO+Hi4qLSzhx1DHfu3MGIESOgr6+Pr7/+GkVFRdi8eTO6desm9WGutG/9+vXYsWMHPvzwQxQXF2PDhg3YuHEjtm3bJvVhnroYQZ2Gu7u7eP3111XaHBwcxPLly7UUEf1WRUWFACBSU1OFEEI0NjYKS0tL8cEHH0h9fvrpJ6FQKMSOHTu0FWaXVF1dLQYMGCCSkpKEt7e3CAsLE0IwRx3JsmXLxMiRIx+5nrnqGMaOHSvmzJmj0jZ58mQxY8YMIQTz1BVxZrOTePDgAbKzs+Hr66vS7uvri++++05LUdFvVVZWAgCeffZZAEBpaSlu3rypkje5XA5vb2/m7QkLCQnB2LFj4ePjo9LOHHUcCQkJcHV1xauvvoqePXtCqVRi165d0nrmqmMYOXIk/vvf/+Ls2bMAgNOnT+PEiRPw9/cHwDx1RXraDoDax48//oiGhgZYWFiotFtYWODmzZtaiop+TQiBxYsXY+TIkXBycgIAKTdN5e3SpUtPPMau6sCBA8jJyUFWVpbaOuao47h48SJiYmKwePFirFy5EpmZmQgNDYVcLsfMmTOZqw5i2bJlqKyshIODA3R1ddHQ0IC1a9di2rRpAPgz1RWx2OxkZDKZyrIQQq2NtGPhwoXIz8/HiRMn1NYxb9pz5coVhIWF4ejRozA0NHxkP+ZI+xobG+Hq6op169YBAJRKJQoLCxETE4OZM2dK/Zgr7Tp48CC++OIL7Nu3D4MHD0ZeXh7Cw8PRq1cvBAcHS/2Yp66Dl9E7CXNzc+jq6qrNYlZUVKj99UhP3ptvvomEhASkpKSgT58+UrulpSUAMG9alJ2djYqKCgwbNgx6enrQ09NDamoq/v73v0NPT0/KA3OkfVZWVhg0aJBKm6Ojo/QQJH+eOoYlS5Zg+fLleO211+Ds7IygoCAsWrQIUVFRAJinrojFZidhYGCAYcOGISkpSaU9KSkJnp6eWoqKhBBYuHAh4uPjkZycDFtbW5X1tra2sLS0VMnbgwcPkJqayrw9IS+99BIKCgqQl5cnfVxdXREYGIi8vDz079+fOeogRowYofbqsLNnz8LGxgYAf546itraWujoqJYXurq60quPmKcuSIsPJ1E7O3DggNDX1xeffvqpKCoqEuHh4cLExESUlZVpO7Qu64033hAKhUIcO3ZM3LhxQ/rU1tZKfT744AOhUChEfHy8KCgoENOmTRNWVlaiqqpKi5F3bb9+Gl0I5qijyMzMFHp6emLt2rXi3LlzIi4uThgbG4svvvhC6sNcaV9wcLDo3bu3OHLkiCgtLRXx8fHC3NxcLF26VOrDPHUtLDY7mY8++kjY2NgIAwMDMXToUOkVO6QdAJr8xMbGSn0aGxvFO++8IywtLYVcLhcvvPCCKCgo0F7QpFZsMkcdx+HDh4WTk5OQy+XCwcFB7Ny5U2U9c6V9VVVVIiwsTPTt21cYGhqK/v37i4iICHH//n2pD/PUtciEEEKbM6tERERE1Hnxnk0iIiIi0hgWm0RERESkMSw2iYiIiEhjWGwSERERkcaw2CQiIiIijWGxSUREREQaw2KTiIiIiDSGxSYRERERaQyLTSJ6osrKyiCTyZCXlwcAOHbsGGQyGe7evQsA2L17N7p166bxcduiX79+2Lp1a5v301btda4ePHgAOzs7pKenP7JPe54/bSsoKECfPn1QU1Oj7VCIugQWm0RPkVmzZmHixIlq7b8t2ABACIGdO3di+PDhMDU1Rbdu3eDq6oqtW7eitrZW6ldVVYWIiAg4ODjA0NAQlpaW8PHxQXx8PB5+wdiLL76I8PBwtXHbo9jx9PTEjRs3oFAo2rQf+v127twJGxsbjBgxQtuhPBHOzs5wd3fHli1btB0KUZfAYpOokwoKCkJ4eDgCAgKQkpKCvLw8rF69GocOHcLRo0cBAHfv3oWnpyf27NmDFStWICcnB8ePH8fUqVOxdOlSVFZWajxOAwMDWFpaQiaTaXys5gghUF9fr9UYtGXbtm2YN2+etsMA8Mss65Mwe/ZsxMTEoKGh4YmMR9SVsdgk6oS+/PJLxMXFYf/+/Vi5ciXc3NzQr18/BAQEIDk5GaNGjQIArFy5EmVlZTh58iSCg4MxaNAg2NvbY/78+cjLy4OpqWmbY8nMzIRSqYShoSFcXV2Rm5ursr6pWVkASExMhKOjI0xNTTFmzBjcuHFDZX1sbCwcHR1haGgIBwcHbN++/XeNm5iYCFdXV8jlcqSlpeHChQsICAiAhYUFTE1N4ebmhm+//bZVx5yVlYXRo0fD3NwcCoUC3t7eyMnJUekjk8nwySefYNKkSTA2NsaAAQOQkJCg0ichIQEDBgyAkZERRo0ahc8//7zJc/Vrhw8fxrBhw2BoaIj+/ftjzZo1zRbROTk5OH/+PMaOHavS/rjzBwBFRUXw9/eHqakpLCwsEBQUhB9//FFaX11djcDAQJiYmMDKygpbtmxRmyXv168f3n//fcyaNQsKhQLz588HAHz33Xd44YUXYGRkBGtra4SGhqpc9n7w4AGWLl2K3r17w8TEBMOHD8exY8ek9ZcuXcL48ePRvXt3mJiYYPDgwfjqq6+k9X5+fvjf//6H1NTUR54bImofLDaJOqG4uDgMHDgQAQEBautkMhkUCgUaGxtx4MABBAYGolevXmr9TE1Noaen16Y4ampqMG7cOAwcOBDZ2dmIjIzEX//618duV1tbi02bNmHv3r04fvw4Ll++rLLdrl27EBERgbVr16K4uBjr1q3D6tWr8fnnn7d63KVLlyIqKgrFxcVwcXHBvXv34O/vj2+//Ra5ubnw8/PD+PHjcfny5RYfd3V1NYKDg5GWloaMjAwMGDAA/v7+qK6uVum3Zs0aTJkyBfn5+fD390dgYCBu374N4Jd7JF955RVMnDgReXl5WLBgASIiIpodNzExETNmzEBoaCiKiorw8ccfY/fu3Vi7du0jtzl+/Djs7e3xzDPPSG0tOX83btyAt7c3hgwZglOnTuGbb75BeXk5pkyZIvVZvHgx0tPTkZCQgKSkJKSlpakV3QCwceNGODk5ITs7G6tXr0ZBQQH8/PwwefJk5Ofn4+DBgzhx4gQWLlwobTN79mykp6fjwIEDyM/Px6uvvooxY8bg3LlzAICQkBDcv38fx48fR0FBAdavX6/yx5OBgQGef/55pKWlNXtOiagdCCJ6agQHBwtdXV1hYmKi8jE0NBQAxJ07d4QQQjg6OooJEyY0u6/y8nIBQERHRz92XG9vb6Gvr682rlwuFwqF4pHbffzxx+LZZ58VNTU1UltMTIwAIHJzc4UQQqSkpKjEHhsbKwCI8+fPS9t89NFHwsLCQlq2trYW+/btUxnrvffeEx4eHq0e99///vdjj3/QoEFi27Zt0rKNjY3YsmXLY7d7qL6+XpiZmYnDhw9LbQDEqlWrpOV79+4JmUwmvv76ayGEEMuWLRNOTk4q+4mIiFA7V78+/15eXmLdunUq2+zdu1dYWVk9MrawsDDxpz/9SaWtJedv9erVwtfXV2W7K1euCACipKREVFVVCX19ffGPf/xDWn/37l1hbGwswsLCpDYbGxsxceJElf0EBQWJv/zlLyptaWlpQkdHR9TV1Ynz588LmUwmrl27ptLnpZdeEitWrBBCCOHs7CwiIyMfedxCCDFp0iQxa9asZvsQUdu1bdqCiJ64UaNGISYmRqXt5MmTmDFjhrQshHjsPZDi/x/+aem9koGBgWoza/Hx8Vi3bt0jtykuLsbzzz8PY2Njqc3Dw+OxYxkbG+O5556Tlq2srFBRUQEAuHXrFq5cuYK5c+dKl1wBoL6+XnrIqDXjurq6qizX1NRgzZo1OHLkCK5fv476+nrU1dW1amazoqICb7/9NpKTk1FeXo6GhgbU1taq7cPFxUX6t4mJCczMzKTjLCkpgZubm0p/d3f3ZsfNzs5GVlaWykxmQ0MDfvrpJ9TW1qqcj4fq6upgaGio0taS85ednY2UlJQmb7W4cOEC6urq8PPPP6vErFAoMHDgQLX+v81BdnY2zp8/j7i4OKlNCIHGxkaUlpbihx9+gBAC9vb2Ktvdv38fPXr0AACEhobijTfewNGjR+Hj44OXX35Z5XwDgJGRkcrDckSkGSw2iZ4yJiYmsLOzU2m7evWqyrK9vT2Ki4ub3c8f/vAHdO/e/bH9HlIoFGrj9uzZs9ltHha0raWvr6+yLJPJpH01NjYC+OVS+vDhw1X66erqtnpcExMTleUlS5YgMTERmzZtgp2dHYyMjPDKK6+06sGVWbNm4datW9i6dStsbGwgl8vh4eGhto+mjvPh8TX1B8PjjquxsRFr1qzB5MmT1db9tqB8yNzcHAUFBa0a5+FY48ePx/r169XWWVlZSZezW3IMv81BY2MjFixYgNDQULW+ffv2RX5+PnR1dZGdnS3l/KGHxe+8efPg5+eH//znPzh69CiioqKwefNmvPnmm1Lf27dvq/xRQ0SawXs2iTqh6dOn4+zZszh06JDaOiEEKisroaOjg6lTpyIuLg7Xr19X61dTU9Pmp7MHDRqE06dPo66uTmrLyMho0z4tLCzQu3dvXLx4EXZ2diofW1vbNo+blpaGWbNmYdKkSXB2doalpSXKyspaFWNaWhpCQ0Ph7++PwYMHQy6Xqzw40xIODg7IyspSaTt16lSz2wwdOhQlJSVq58XOzg46Ok3/d69UKnHmzBmVIrAl52/o0KEoLCxEv3791MYyMTHBc889B319fWRmZkrbVFVVSUXo446jsLCwyeMwMDCAUqlEQ0MDKioq1NZbWlpK+7G2tsbrr7+O+Ph4vPXWW9i1a5fKOD/88AOUSuVj4yGitmGxSdQJTZkyBVOnTsW0adMQFRWFU6dO4dKlSzhy5Ah8fHyQkpICAFi3bh2sra0xfPhw7NmzB0VFRTh37hw+++wzDBkyBPfu3WtTHNOnT4eOjg7mzp2LoqIifPXVV9i0aVObjy8yMhJRUVH429/+hrNnz6KgoACxsbGIjo5u87h2dnaIj49HXl4eTp8+jenTp0uzjS1lZ2eHvXv3ori4GCdPnkRgYCCMjIxatY8FCxbgzJkzWLZsGc6ePYsvv/wSu3fvBvDoWx/efvtt7NmzB5GRkSgsLERxcTEOHjyIVatWPXKcUaNGoaamBoWFhVJbS85fSEgIbt++jWnTpiEzMxMXL17E0aNHMWfOHDQ0NMDMzAzBwcFYsmQJUlJSUFhYiDlz5kBHR+ext24sW7YM33//PUJCQpCXl4dz584hISFBmpW0t7dHYGAgZs6cifj4eJSWliIrKwvr16+XnjgPDw9HYmIiSktLkZOTg+TkZDg6OkpjlJWV4dq1a/Dx8Wk2FiJqOxabRJ2QTCbDvn37EB0djX/961/w9vaGi4sLIiMjERAQAD8/PwBA9+7dkZGRgRkzZuD999+HUqmEl5cX9u/fj40bN7b5ReumpqY4fPgwioqKoFQqERER0eRl19aaN28ePvnkE+zevRvOzs7w9vbG7t27pZnNtoy7ZcsWdO/eHZ6enhg/fjz8/PwwdOjQVsX32Wef4c6dO1AqlQgKCkJoaOhjbzn4LVtbW/zzn/9EfHw8XFxcEBMTI90zK5fLm9zGz88PR44cQVJSEtzc3PDHP/4R0dHRsLGxeeQ4PXr0wOTJk1Xuj2zJ+evVqxfS09PR0NAAPz8/ODk5ISwsDAqFQppFjY6OhoeHB8aNGwcfHx+MGDFCel1Vc1xcXJCamopz587By8sLSqUSq1evhpWVldQnNjYWM2fOxFtvvYWBAwdiwoQJOHnyJKytrQH8cq9qSEgIHB0dMWbMGAwcOFDl9Vj79++Hr69vs+eGiNqHTPzem6qIiOiJWrt2LXbs2IErV660634LCgrg4+OD8+fPw8zMrF33/Ws1NTXo3bs3Nm/ejLlz52psnMe5f/8+BgwYgP3793eZb00i0iY+IERE1EFt374dbm5u6NGjB9LT07Fx40aVd022F2dnZ2zYsAFlZWVwdnZut/3m5ubizJkzcHd3R2VlJd59910AaPL9r0/SpUuXEBERwUKT6AnhzCYRUQe1aNEiHDx4ELdv30bfvn0RFBSEFStWtPll+09Kbm4u5s2bh5KSEhgYGGDYsGGIjo5u14KWiDo+FptEREREpDF8QIiIiIiINIbFJhERERFpDItNIiIiItIYFptEREREpDEsNomIiIhIY1hsEhEREZHGsNgkIiIiIo1hsUlEREREGvN/+XRyeE/4fREAAAAASUVORK5CYII=\n" }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "\n", "e1 = [p.energy(command='MRCI', stateSymmetry='1') for p in projects]\n", "e2 = [p.energy(command='MRCI', stateSymmetry='4') for p in projects]\n", "cc = [p.variable('RCC') for p in projects]\n", "fig, ax = plt.subplots()\n", "ax.plot(grid, e1, grid, e2)\n", "ax.legend([r'$X\\;{}^1A_1$', r'$V\\;{}^1B_1$'], loc='center left')\n", "plt.xlabel('HCCH dihedral angle (degrees)')\n", "ax.set_ylabel(r'Energy $(E_h)$')\n", "ax2 = ax.twinx()\n", "ax2.plot(grid, cc, color=\"black\")\n", "ax2.set_ylabel(r\"C-C bond length $(\\rm\\AA)$\")\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "name": "python3", "language": "python", "display_name": "Python 3 (ipykernel)" }, "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.10.6" }, "toc": { "base_numbering": 1, "nav_menu": {}, "number_sections": true, "sideBar": true, "skip_h1_title": false, "title_cell": "Table of Contents", "title_sidebar": "Contents", "toc_cell": false, "toc_position": {}, "toc_section_display": true, "toc_window_display": false }, "varInspector": { "cols": { "lenName": 16, "lenType": 16, "lenVar": 40 }, "kernels_config": { "python": { "delete_cmd_postfix": "", "delete_cmd_prefix": "del ", "library": "var_list.py", "varRefreshCmd": "print(var_dic_list())" }, "r": { "delete_cmd_postfix": ") ", "delete_cmd_prefix": "rm(", "library": "var_list.r", "varRefreshCmd": "cat(var_dic_list()) " } }, "types_to_exclude": [ "module", "function", "builtin_function_or_method", "instance", "_Feature" ], "window_display": false } }, "nbformat": 4, "nbformat_minor": 4 }