{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Create a vmd movie from .xyz file" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Create molecule and scan two coordinates\n", "\n", "Here we create a molecule from a smiles string and visualize its graph and 3D representation.\n", "We select two bonds to rotate around and scan the coordinates.\n", "The scan is written to an `xyz` file which can be read by `vmd`." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "['Bond(0, 1)', 'Bond(1, 2)', 'Bond(2, 3)', 'Bond(3, 4)', 'Bond(4, 5)', 'Bond(5, 6)', 'Bond(5, 11)', 'Bond(6, 7)', 'Bond(7, 8)', 'Bond(7, 8)', 'Bond(8, 9)', 'Bond(11, 12)', 'Bond(12, 13)']\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAApQAAAHzCAYAAACe1o1DAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAD8eElEQVR4nOzdd3yN5//48dfJnkLslSgRxAhiz9p7FdUapVatSn2MqhYJ1aralNo1Yu8RUbuImdhBCCGIkZC9z7l/f/g6v6YZMk6WvJ+Px3lUzn1d1/2+D2neue77el8qRVEUhBBCCCGEyCC9nA5ACCGEEELkbZJQCiGEEEKITJGEUgghhBBCZIoklEIIIYQQIlMkoRRCCCGEEJkiCaUQQgghhMgUSSiFEEIIIUSmSEIphBBCCCEyRRJKIYQQQgiRKZJQCiGEEEKITJGEUgghhBBCZIoklEIIIYQQIlMkoRRCCCGEEJkiCaUQQgghhMgUSSiFEEIIIUSmSEIphBBCCCEyRRJKIYQQQgiRKZJQCiGEEEKITJGEUgghhBBCZIoklEIIIYQQIlMkoRRCCCGEEJkiCWUeU65cORYuXJhqG5VKxd69ezN8jlOnTqFSqQgJCcnwGBnl7++PSqXi2rVr2X5uIYQQQmRMvk0oBw0ahEqlQqVSYWBggI2NDSNHjuTt27dpHiMzyU9mkz4hhBBCiNwi3yaUAO3btycwMBB/f39Wr17NgQMHGDVqVE6HJYQQQgiRp+TrhNLY2JgSJUpQpkwZ2rZtS58+ffj777+1xzUaDTNmzKBMmTIYGxtTs2ZNPDw8tMc/+eQTAGrVqoVKpeLTTz8F4PLly7Rp04YiRYpgZWVF8+bN8fb21vYrV64cAD169EClUmm/9vPzo1u3bhQvXhwLCwvq1q3LsWPHksQdHh5O3759sbCwoFSpUixZsiTV63z27Bl9+vShUKFCFC5cmG7duuHv7//Bz+fcuXM4OjpiYmJC/fr1uXnzZqLju3btomrVqhgbG1OuXDnmzZuX6Hi5cuX45ZdfGDx4MJaWltjY2LBy5cpEbS5dukStWrUwMTGhTp06XL16NdHxt2/f0q9fP4oWLYqpqSkVK1Zk3bp1H4xdCCGEENknXyeU//bw4UM8PDwwNDTUvrdo0SLmzZvH3LlzuXHjBu3ataNr167cv38feJcMARw7dozAwEB2794NvEv4Bg4cyJkzZ7hw4QIVK1akY8eOhIeHA+8SToB169YRGBio/ToiIoKOHTty7Ngxrl69Srt27ejSpQtPnjxJFOvvv/9OjRo18Pb25ocffmDcuHEcPXo02euKioqiRYsWWFhY8M8//3D27FksLCxo3749cXFxqX4mEydOZO7cuVy+fJlixYrRtWtX4uPjAfDy8uLzzz/niy++4ObNm7i4uDB16lT++uuvRGPMmzdPmyiOGjWKkSNHcvfuXQAiIyPp3LkzlSpVwsvLCxcXFyZMmJCo/9SpU/Hx8eHw4cPcuXOH5cuXU6RIkVTjFkIIIUQ2U7KIk5OTUrp06Wx9OTk5pTm+gQMHKvr6+oq5ubliYmKiAAqgzJ8/X9umVKlSyqxZsxL1q1u3rjJq1ChFURTl0aNHCqBcvXo11XMlJCQolpaWyoEDB7TvAcqePXs+GKeDg4OyZMkS7de2trZK+/btE7Xp06eP0qFDh2THXrNmjVKpUiVFo9Foj8fGxiqmpqbKkSNHkj3nyZMnFUDZunWr9r3g4GDF1NRU2bZtm6IoitK3b1+lTZs2ifpNnDhRcXBwSBRr//79tV9rNBqlWLFiyvLlyxVFUZQVK1Yo1tbWSmRkpLbN8uXLE32mXbp0Ub7++uuUPyAhhBBC5DiDrEpUX7x4wbNnz7JqeJ1o0aIFy5cvJyoqitWrV+Pr68u3334LQFhYGM+fP6dx48aJ+jRu3Jjr16+nOu6rV6+YNm0aJ06c4OXLl6jVaqKiopLMNP5XZGQkrq6uHDx4kOfPn5OQkEB0dHSSfg0bNkzydUorv728vHjw4AGWlpaJ3o+JicHPzy/VeP59HmtraypVqsSdO3cAuHPnDt26dUvUvnHjxixcuBC1Wo2+vj4ANWrU0B5XqVSUKFGCV69eacdwdHTEzMwsxWsbOXIkPXv2xNvbm7Zt29K9e3caNWqUatxCCCGEyF5ZllCWKFEiq4bW2TnNzc2xs7MDYPHixbRo0QJXV1dmzpypbaNSqRL1URQlyXv/NWjQIF6/fs3ChQuxtbXF2NiYhg0bpukW85EjR5g7dy52dnaYmprSq1evD/ZLLs73NBoNTk5OuLm5JTlWtGjRD46b0nmS+xwURUnS/t+PELzvr9FoUmz/Xx06dODx48ccOnSIY8eO0apVK0aPHs3cuXPTHbsQQgghskaWJZRXrlzJqqGzzPTp0+nQoQMjR46kVKlSlCpVirNnz9KsWTNtG09PT+rVqweAkZERAGq1OtE4Z86cYdmyZXTs2BGAgIAAgoKCErUxNDRMtt+gQYPo0aMH8O6ZyuQWz1y4cCHJ15UrV072mmrXrs22bdsoVqwYBQoU+NBHkGRcGxsb4N3iGF9fX+15HBwcOHv2bKL2np6e2Nvba2cnP8TBwYGNGzcSHR2NqalpstcG7xLfQYMGMWjQIJo2bap9tlMIIYQQuYMsyvmXTz/9lKpVq/LLL78A72YMf/vtN7Zt28a9e/eYPHky165dw9nZGYBixYphamqKh4cHL1++JDQ0FAA7Ozs2btzInTt3uHjxIv369dMmTO+VK1eO48eP8+LFC23tSzs7O3bv3s21a9e4fv06ffv21c7m/du5c+eYM2cOvr6+/PHHH+zYsUMb03/169ePIkWK0K1bN86cOcOjR484ffo0zs7OPH36NNXPY8aMGRw/fpxbt24xaNAgihQpQvfu3QEYP348x48fZ+bMmfj6+rJ+/XqWLl2aZFFNavr27Yuenh5DhgzBx8cHd3f3JInitGnT2LdvHw8ePOD27dscPHiQKlWqpPkcQgghhMgGOfsIZ84ZOHCg0q1btyTvu7m5KUZGRsqTJ08UtVqtuLq6KqVLl1YMDQ0VR0dH5fDhw4nar1q1Silbtqyip6enNG/eXFEURfH29lbq1KmjGBsbKxUrVlR27Nih2NraKgsWLND2279/v2JnZ6cYGBgotra2iqK8W+TTokULxdTUVClbtqyydOlSpXnz5oqzs7O2n62treLq6qp8/vnnipmZmVK8eHFl4cKFiWLiPwt+AgMDla+++kopUqSIYmxsrJQvX14ZNmyYEhoamuxn835RzoEDB5SqVasqRkZGSt26dZVr164lardz507FwcFBMTQ0VGxsbJTff/890fH/XrOiKIqjo6Myffp07dfnz59XHB0dFSMjI6VmzZrKrl27Ei3KmTlzplKlShXF1NRUsba2Vrp166Y8fPgw2biFEEIIkTNUipKGB9mEEEIIIYRIgdzyFkIIIYQQmSIJpRBCCCGEyBRJKIUQQgghRKZIQimEEEIIITIly+pQCiGEEHldYGAgJ0+e1O46ptFoMDQ0xMHBAScnJ1q2bImVlVVOh5mvaTQafH198fPzIz4+HjMzMxwcHChduvQHNyIRuiOrvIUQQoj/uHDhAosWLUJRFDp16oSTkxP29vYYGBgQGxvLrVu3uHLlCocOHaJkyZJMmDCBihUr5nTY+YZGo+Ho0aOsXbuWqKgoKlasiL29PUZGRkRERHDr1i2ePn2KjY0NI0eOxNHRMadD/uhJQimEEEL8n8jISKZMmUJMTAzTpk2jdOnSH+zj4+PDTz/9ROPGjfnuu+/SvFuYyJjLly/zww8/0KpVK4YMGUKxYsVSbHvnzh2WLVvGy5cvWbhwIaVKlcrGSPMXSSiFEEII4MWLF/Tv359JkybRtm3bdPVVFIU1a9bg4eHBxo0bk+yOJjJPURRmzZrFvXv3WLBgAUWKFElz31u3bjF+/HjGjBlDly5dsjDK/EsSSiGEEPne27dv6d27N0uXLqVy5coZHufvv/9mxYoVbNu2DQMDWaagK4qiMHHiREqVKsX//ve/DI0RFxfH4MGD6datG71799ZxhEJWeSfj008/5bvvvsvSc7i4uFCzZs0Uj586dQqVSkVISEiWxpEcf39/VCoV165dy/ZzCyFETnB2dmbWrFmZSiYB2rZtS9euXfn99991FJkAWLFiBYUKFcpwMglgZGTEX3/9xZYtW/D29tZhdALycUI5aNAgVCpVkteDBw/YvXs3M2fOzOkQhRBCZIO9e/dSpkwZ6tevT3h4uPaWd9GiRVGpVLi4uCTpoygKq1atwsnJiQIFClC4cGGaN2/OoUOH+Oqrr/D29sbHxyf7L+Yj9OjRIw4cOMAPP/yQ6bEMDAxYsWIFkyZNIi4uTgfRiffybUIJ0L59ewIDAxO9PvnkE6ytrbG0tMzp8IQQQmQxRVFYtmwZ06ZNAyA4OJiVK1cSGxtL9+7dU+w3ffp0hg8fTr169di1axd//fUXxsbGdO7cmT179vDbb78xZ86cbLqKj9uPP/7IokWL0NPTS3PCn9KkUeXKlSlatCiDBg3izz//zP6L+Yjl64TS2NiYEiVKJHrp6+snuuV99+5dzMzM2Lx5s7bf7t27MTEx4ebNmwCEhoYyfPhwihUrRoECBWjZsiXXr19PdK7Zs2dTvHhxLC0tGTJkCDExMWmK8dy5czg6OmJiYkL9+vW153xv165dVK1aFWNjY8qVK8e8efMSHS9Xrhy//PILgwcPxtLSEhsbG1auXJmozaVLl6hVqxYmJibUqVOHq1evJjr+9u1b+vXrR9GiRTE1NaVixYqsW7cuTfELIURuduHCBZycnDAxMQHA1taWt2/fcvr0aX799dcU+61du5YmTZqwfPly2rRpQ5cuXdi/fz9WVlasX7+e8uXLExYWRlBQUHZdykfp2bNnKIqCnZ0dkPaEH8DU1JTz588nem3btg2AL774gn379iHLSHRHnhj+gMqVKzN37lxGjRpF48aNMTQ0ZNiwYcyePZvq1atra5RZW1vj7u6OlZUVK1asoFWrVvj6+mJtbc327duZPn06f/zxB02bNmXjxo0sXryY8uXLf/D8EydOZNGiRZQoUYIpU6bQtWtXfH19MTQ0xMvLi88//xwXFxf69OmDp6cno0aNonDhwgwaNEg7xrx585g5cyZTpkxh586djBw5kmbNmlG5cmUiIyPp3LkzLVu2ZNOmTTx69AhnZ+dEMUydOhUfHx8OHz5MkSJFePDgAdHR0br+qIUQItvt2bOHAQMGaL9OayFsQ0PDJAXNTUxMtC+Anj17cvjw4UTjZ0ZISAi7d+/m0qVLBAQEoKf3bk5Io9FQvnx5GjZsSPfu3TEzM9PJ+XKDjRs3MnjwYO3X7xN+lUpFUFAQq1evTrGvnp4eDRo0SPaYgYEBTZo04dy5czRp0kTncedLShaJi4tTnj9/nq2vuLi4NMc3cOBARV9fXzE3N9e+evXqpSiKojRv3lxxdnZO1L5Tp05K06ZNlVatWilt2rRRNBqNoiiKcvz4caVAgQJKTExMovYVKlRQVqxYoSiKojRs2FAZMWJEouP169dXHB0dU4zv5MmTCqBs3bpV+15wcLBiamqqbNu2TVEURenbt6/Spk2bRP0mTpyoODg4aL+2tbVV+vfvr/1ao9EoxYoVU5YvX64oiqKsWLFCsba2ViIjI7Vtli9frgDK1atXFUVRlC5duihff/11irEKIURe1aVLFyUhISHZY69fv1YAZfr06UmOLViwQNHX11dWr16tvHnzRnn+/Lkybtw4xcTERDl79qyiKIpy7969JD9LMuLZs2fK8OHDlS5duiibN29WHj9+rP0ZpCjv/r9+//59ZfXq1Uq7du2U//3vf8qbN28yfd7c4LPPPkv08+nfUvv7GThwoGJubp7q2MeOHVN+//13XYQpFEXJshnKoKCgJLdWs9rw4cMpWbJkmtu3aNGC5cuXa782NzdPse3atWuxt7dHT0+PW7duaX+L9fLyIiIigsKFCydqHx0djZ+fH/CusOqIESMSHW/YsCEnT578YIwNGzbU/tna2ppKlSpx584d7bjdunVL1L5x48YsXLgQtVqtLa5bo0YN7XGVSkWJEiV49eqVdgxHR8dEv9H++5wAI0eOpGfPnnh7e9O2bVu6d+9Oo0aNPhi7EELkdoqiZKgQ+XfffYepqSmjR49m6NChwLv/Rx84cIDGjRsDYGdnp/05kNHYNm7ciJubG3PmzElxtxeVSoWdnR12dnYMGTKEM2fO0LNnT8aPH0+nTp0yfP7cICYmJsMzrtHR0ZQoUYLXr19TsmRJunfvzowZM7C2tgagdu3arFmzRpfh5mtZllAWKVKE4cOHZ9XwKZ4zPczNzbXPZXzI9evXiYyMRE9PjxcvXmir7Ws0GkqWLMmpU6eS9ClYsGC64kmr98msoihJbs8oyTwPYmhomKS/RqNJsf1/dejQgcePH3Po0CGOHTtGq1atGD16NHPnzs3oJQghRK6Q0b2e161bh7OzM2PGjKFDhw7ExcWxYcMGunXrxu7du2nXrp32lnRGKIrC5MmTMTIy4tChQ+mqadm0aVMOHTrEuHHjePDgQZLHmPKSjO465OjoiKOjI9WqVQPg9OnTLFiwgOPHj3P58mUsLCwoVKgQ4eHhugw3X8uyhNLQ0DBds4W52Zs3bxg0aBA//vgjL168oF+/fnh7e2Nqakrt2rV58eIFBgYGlCtXLtn+VapU4cKFC3z11Vfa9y5cuJCmc1+4cAEbGxvg3eIYX19fbZ00BwcHzp49m6i9p6cn9vb2af4mdHBwYOPGjURHR2t3dkgutver4gYNGkTTpk2ZOHGiJJRCiDwvLb9U/9fbt2+1M5P//v9ghw4d+PTTTxkxYgSPHj3K8PgAM2bMoFixYowfPz5D/U1NTVm+fDnjx49nzZo1DBkyJEPj5FXjxo1L9HWbNm2oVasWvXr1YtWqVUmOi8zL16u802rEiBGULVuWn376ifnz56MoChMmTACgdevW2gehjxw5gr+/P56envz0009cuXIFeFcwd+3ataxduxZfX1+mT5/O7du303TuGTNmcPz4cW7dusWgQYMoUqSIdmXb+PHjOX78ODNnzsTX15f169ezdOlSbWxp0bdvX/T09BgyZAg+Pj64u7snSRSnTZvGvn37ePDgAbdv3+bgwYNUqVIlzecQQojcSqPRpDvpu3fvHtHR0dStWzfJsTp16uDv709ERAQvX75M950zeFfd48mTJxlOJt9TqVTMmzeP/fv3c//+/UyNlVPe303ThR49emBubq6dNImIiPioFjDlNEkoP2DDhg24u7uzceNGDAwMMDMzw83NjdWrV+Pu7o5KpcLd3Z1mzZoxePBg7O3t+eKLL/D396d48eIA9OnTh2nTpvH999/j5OTE48ePGTlyZJrOP3v2bJydnXFyciIwMJD9+/djZGQEvHv+Y/v27WzdupVq1aoxbdo0ZsyYkWiF94dYWFhw4MABfHx8qFWrFj/++CO//fZbojZGRkb88MMP1KhRg2bNmqGvr8/WrVvTfA4hhMitSpcuzbNnz9LV5/0jT/+9m6MoChcuXKBQoUKYm5vj5eWFk5NTusaOi4tj6tSpzJs3L801F5Ort/j+VaVKFRYvXsy4cePyZIkcQ0NDYmNjdTaeoijaRxGuXr2a6o51Ip1yZCmQEEIIkQts3bpVW/XiPXd3d2XHjh3K2rVrFUDp3bu3smPHDmXHjh3aFcefffaZoqenpzg7OytHjhxR9u/fr/Ts2VMBlJkzZyqKoiijRo1Sbt26le54/vjjD0VRFOXRo0eKlZWV0qxZM2Xo0KEprmg+f/58ktfChQsVQJk8ebKiKIoyefJk5dy5c+n9eHLc9OnTtavm/yu1Vd7J2bZtmwIoCxcuVBRFUX799Vfl2LFjugo135OEUgghRL4VExOTqBScorwrtwYk+3r06JGiKIoSHR2t/P7770qNGjUUS0tLxdraWmnQoIGyadMmRaPRKGFhYUrHjh3THU+nTp2U0NBQRVHelQN6H1d6k6dBgwYpKpVKuX//vqIo75LTAQMGpDuenHb//v0kZes+lPD7+/srjRo1UhYvXqy4u7srhw8fViZPnqyYmJgoVatWVSIiIhS1Wq20bNlSiY+Pz6Er+/hIYXMhhBD5lrGxMU2bNmX37t307NkTAH9//w/2MzExYcKECSk+sz5v3jyGDRuWrljCw8MxMzOjQIECQMZXoIeHh7Njxw6aN2+urWRSrlw5goKCkq0OkpvZ2dnx9u1bAgMDtQt9R44cyePHj7VtduzYwY4dO4B3+35bWVlRvHhx5s+fz8uXL1Gr1dja2jJ27FimTJmCubk5Bw8epGXLlulaPS9SJ5+kEEKIfG3SpEl07NiR5s2bZ2gRzX9dv36du3fvMn369HT1u3btGrVq1cr0+bdu3UpkZKS2PuZ7n3zyCY8fP06xIkluNWPGDJydndm2bRsqlSpNCf/u3btTPBYeHs78+fPZv3+/DqMUsihHCCFEvmZsbMzvv//O4MGDiYmJydRYL168YNy4cSxatCjdM4F3797FwcEhU+cHWLNmDQULFtTOuL7n4ODA3bt3Mz1+dqtevTqOjo6pbrOYVoqiMHbsWKZOnYqFhYUOohPvSUIphBAi36tduzZjx46lT58+vHnzJkNjPHz4kH79+rFs2TJtlY/0iI2N1dYDzqjbt29z8eJF+vXrp91T/D1TU9NMJ8w5ZfLkyZw9ezZTFUY0Gg3jxo3DycmJFi1a6DA6AZJQCiGEEMC7usJTp06ld+/eHDhwIM39NBoNS5cu1dYcfr/5RHoZGxsTHR2dob7vvd9K8L+3u+HdVoT/TTLzCn19fdasWcOZM2cYP358uj+n58+f8/nnn1O1alXGjBmTRVHmb5JQCiGEEP+nTp06HDx4kIsXL9KhQwfWr19PSEhIsm0DAwOZN28edevWJSgoiH379mFra5vhc1euXBkfH58M94+Li2Pjxo04OTklW1/Rx8cnw8lubmBgYMCSJUtISEigWbNmrF+//oMzrq9fv2b27NkMHjyYGTNmpHuhlEg7WZQjhBBC/IupqSk///wzUVFRbN26lVGjRhEaGoq+vj4GBgbEx8ej0WgoXrw4Xbp0YdSoUTpZPV2rVi2WLFmS4f779+8nKCiIGTNmJHvc398/UwlvbuDt7Y21tTWbNm3i3LlzdO/eHWtra2rXro29vT3GxsZERERw8+ZNrl+/jr6+PkOGDGHSpEmZ2ltdfJgklEIIIUQyzMzMGDx4MIMHDwbeLehISEjAwMAgUfL44MED3Nzc8Pf355NPPsnw+SwsLIiKiiIsLExbOujw4cNERkYSHh4OvJtl3LlzJwAdO3ZMtHXgmjVrMDU1pW/fvknG9vf3p3DhwnmqZNB/RUZGcvz4cWrWrEmlSpWoVKkSgwcP5u3bt1y9epUHDx4QFxeHmZkZHTt2ZPLkyXn2Fn9epFKUPLgXkxBCCJFLKIrCypUrMTMzY8CAAZkaa9u2bQQHBzNq1CjgXf3If9dc/LdHjx5pSwAFBARQrlw5+vfvz/r165O0/eGHH+jSpQuNGjXKVHw5ad++fdy7d48xY8bIHty5kMz/CiGEEJmgUqlo0qQJDx8+TPe+4P/Vo0cPdu7cqX1u09/fH+XdrnZJXv+uJ1m2bFnUanWyyeTjx4+5efMmDRs2zFRsOenx48dcu3aNVq1aSTKZS0lCKYQQQmRSlSpVsLa25uzZs5kax8jIiJkzZzJ+/HidxPW+7uKCBQvy7O1utVqNu7s7pUuXpnbt2jkdjkiBJJRCCCFEJunp6dG4cWPu3r3L69evMzVW48aNsbGxYd68eZkaR1EUxo8fT9euXalYsWKmxspJFy9e5PXr13Tq1CnPJsX5gSSUQgghhA44OjpiaWnJuXPnMj3WtGnTePXqFVOnTiUhISHd/aOjoxk5ciTlypVjyJAhmY4np4SFhXHq1Cnq1q2r3ctb5E6SUAohhBA6oK+vT8OGDbl58yahoaGZGkulUjF79mzs7e3p1KkT169fT3PfM2fO0KlTJ7p06cLYsWMzFUdOO3LkCEZGRrKzTR4gCaUQQgihI05OThgbG+Pp6ZnpsVQqFQMGDGDdunUsX76czp07s3HjRh4/fsy/C7RoNBru37/PmjVraN++Pfv27WPXrl106tQp0zHkpAcPHuDj40O7du2k/E8eIGWDhBBCCB06deoU586d47vvvsPc3Fxn454/f565c+dSqFAhAgMDtYW6NRoNFSpUoEGDBnTv3v2jWAWdkJDAsmXLsLKy4quvvpJnJ/MAKWwuhBBC6FD9+vXx9PTk4sWLtGzZUmfjBgYG0rZtW7755hudjZlbnT17ltDQUL788ktJJvMIueUthBBC6JCpqSlOTk5cunSJ2NhYnYyp0Wh48OBBnl6tnVZv3rzh7NmzNGrUiKJFi+Z0OCKNJKEUQgghdKxhw4bEx8dz+fJlnYz39OlToqOjsbe318l4uZWiKBw+fBgLCwuaNWuW0+GIdJCEUgghhNCxAgUK4OjoyIULF4iPj8/0eL6+vpiZmVG6dGkdRJd73b17lwcPHtChQwcMDQ1zOhyRDpJQCiGEEFmgcePGREVFce3atUyP5evri729/Uf9PGFcXBweHh7Y29tTqVKlnA5HpJMklEIIIUQWKFy4MA4ODnh6eqLRaDI8ztu3b3n9+vVH//zk6dOniYqKon379jkdisgASSiFEEKILNKkSRNCQkK4detWhse4f/8+enp6VKhQQYeR5S6vXr3iwoULNG3alEKFCuV0OCIDJKEUQgghskiJEiWws7Pj3LlzZLTss6+vL+XKlcPY2FjH0eUOiqJw6NAhChUqRKNGjXI6HJFBklAKIYRIlqIoxMXFZep2rXg3S/nq1St8fX3T3TcuLg5/f/+P+nb3jRs3ePLkCR07dsTAQMpj51XyNyeEEAJ4l0CePXuWrVu34u/vj0qlwtjYGLVaTVxcHKampjRq1IiBAwdSpEiRLIkhLCyMq1ev4ufnR3x8PBYWFlSrVg0HB4c8u+rX1taWsmXLcvbs2XQvrHn48CFqtTrPlQvy8/Nj334PTv3jxbPnr9AoKvRUCmVKF+fTZk5079aBTz75hOjoaP7++2+qVatG+fLlczpskQmSUAohhGD//v0sWbKEBg0aMG7cOCpUqJAk8YmOjubEiROMGDGCQoUK8csvv+ik8HRUVBTbtm1j+/btFChQgNq1a1OxYkWMjIwICwtj69at3Lp1C0tLS4YPH07z5s3z3GrnJk2asGXLFh4/fky5cuXS3M/X15ciRYpgbW2ddcHp0NFjJ3CZsZigsEIkmHbFtGAnDEuVBEAD+EU959Y2b5avc6F4oXA6tm8Aipq2bdvmbOAi02QvbyGEyMdCQ0NxdnamePHiuLq6YmJikqZ+3t7eTJo0idGjR9OjR48Mn//AgQMsWLCAgQMH8vnnn2Nqappi22fPnrFixQquXr3KokWL8tSMlqIorFixAgsLC/r375/mPvPnz6d69eq5PuEKCwtj2IiJeN8xxKTUjxgYfTgBTogN4s2DKdSwj2PH1lVYWFhkQ6Qiq0hCKYQQ+dSbN2/44osvmD59Oo0bN053/7i4OMaOHUv16tUZPXp0uvomJCQwduxYLCwscHV1TTWR/K9Hjx4xduxY+vXrxxdffJHesHPMzZs32b17N8OGDaNUqVIfbP/8+XNWrVrFwIED0zWrmd0CAwNp26E/8VaTMbVO/+42UcEnMImYx1GPzRQrViwLIhTZQRblCCFEPhQXF0e/fv2YM2dOhpJJACMjI5YvX46Pjw/btm1Lcz+1Ws2gQYNo0aIFc+bMSVcyCfDJJ5+wd+9ejh07xsaNG9Mbdo6pWrUqhQoV4ty5c2lq7+vri4mJCWXLls3iyDIuJCSENh36kVBkYYaSSQCzwi2JLfQ7rdt9SVhYmI4jFNlFZiiFECIfcnV1pVKlSjqZ4UtISKBTp06sW7cuTTNv06dPp3z58gwcODBT51Wr1fTt25cJEyZQt27dTI2VXby8vDh48CCjR4/+4MKmlStXUrhwYXr27JlN0aVf955fcyd4IGaFGmZ6rOg3/+BYeifbNq/QQWS6ERoayu7du7l48SIBAQHo6b2bh9NoNFSoUIGGDRvSvXv3dP9S9DGSGUohhMhn7t69y507dyhatCiDBw+mcuXKmJubU7p0abp164aXl1eSPt7e3rRu3RoLCwsKFizIZ599xsOHDwEwMDBg/vz5TJw48YPn9vb25sGDBwwcOJDw8HAmTZpE27ZtKVq0KCqVChcXlyR9FEVh8eLFVK5cGWNjY0qWLMnIkSMJCwtj2bJlTJ48mdjY2Ex/LtnB0dERCwuLJLOUL1++5N69e9y/f5/w8HDCw8MJDAzM1eWC9u935/bDQikmk+qECJ7emMAtj0+4vt+Kuyfq8fbp9hTHM7VuhrePIUeOHMuqkNPs+fPnjBgxggEDBmBsbMwPP/zAwYMHOXDggPb17bffEhkZSY8ePZg4cSIhISE5HXaOkhlKIYTIZ8aMGcM333zDjBkzCA4Opnfv3jg4OPD69WvmzZvHlStXOHLkCC1btgTeJaD16tWjZs2aTJ48mZiYGKZNm8bbt2+5du2adqX3V199xYwZM1J93q9r166sWbOGokWL4u/vT82aNXF0dMTe3p7Vq1czffr0JEnl+PHjWbhwIRMmTKB169b4+Pgwbdo0KlasyPnz59m5cyevX79m7NixWfWR6ZSnpydHjhzBxsaG7du3ExoaipmZGWZmZiiKQlhYGFFRUahUKn7//Xc+/fTTXLmq3al+R+KLbUTfIPnFNH7nOhEV4kVJh5kYW1Tk7dNtvHm8DlunvyhUNvmZcXV8KKZvh3Dx3IGsDD1FiqKwadMmNm3axJw5c3B0dExTv9OnT+Pq6sqECRPo2LFjFkeZO0lCKYQQ+UhkZCR9+vTh4MGDvHr1KskiiIiICOzs7KhWrRrHjr2bKfr88885efIkfn5+FChQAIDHjx9TsWJFxo0bx2+//Qa8S5QOHjzIL7/8kuy5fXx8WLJkCcuXLwfQ7hyjUqkICgqiaNGiSRLKZ8+eYWtry6hRo1i8eLH2/S1bttC3b19WrlzJ119/Tfv27Tl69GiuTLz+TVEUli1bxsKFC6lVqxZOTk4prqwPDg7G29ubt2/fsmDBAurVq5fN0abs+vXr9B60EvNy85M9HvbCg4cXumNbZz2FyvTRvu93rhMx4XdwaHcflUo/2b6Rj0azb8t4qlSpkiWxp0RRFCZPnoyRkRHTp09Pd5H1qKgovvvuO6pXr863336bRVHmXnLLWwgh8pErV67QtGlTgGRX1FpYWODg4EBAQADw7vnIgwcP0rNnT20yCe+Kdbdo0YI9e/Zo32vYsCHXrl1L8dwbNmxg6NCh2q9VKtUHE8ALFy6gVquTzPp07twZgF27dmFgYECDBg24ePFiqmPltKCgINq3b8+pU6cYNWoUjRs3TrVMU+HChWnTpg2dOnVi4sSJTJw4MdfsWrR56wEUy14pHg8J3IeegQUFSyV+/tPa9iviY54T9eZSin01Fj3Zuj37ZyhnzJhBsWLFmDlzZoZ27DEzM2PFihX4+fmxbt26LIgwd5OEUggh8hEvLy/q1KmT4vHQ0FC8vb2pWrUq8G7Hk+joaGrUqJGkbY0aNXjw4AExMTHAuwTRwMCAhISEZMe+ffs2NWvWTFe8cXFxAEn2sTY0NESlUnHjxg3gXeHwS5dSTlJy2qtXr+jUqRM1atSgSZMm2sUdaWFubs5nn33Gy5cv6d+/f65IKs9fvIaZVc0Uj8eE+WBiURmVXuLEzKRAdQCiw2+n2NfUqjZnPb11EmdanTt3joCAAMaPH5+pcVQqFfPnz2fv3r08ePBAR9HlDZJQCiFEPuLn54ednV2Kx0ePHk1kZCQ//vgj8O62K5DsTi3W1tYoisLbt2+175UpU4bnz58nO7aiKOjrJ3+bMyUODg4ASRaxeHp6oiiKNr7atWvj7Z29SUhaxcXF8dlnn9GyZUvKlCmT4XGcnJwwMDDIdNKjCyGhMegZmKV4XB0XjL5RoSTvG/zfe+q4Nyn2NTAqSFBweOaDTKO4uDimTp3K3Llz07xQDCA+Pl5beN7U1JSCBQvSqFEjLly4wOLFixk3bhz56alCSSiFECIfiY+Px8jIKNljU6dOxc3NjQULFuDk5JToWGq3pv99zMjIiPj4+HSdNzWOjo40a9aM33//nR07dhASEoKnpycjRoxAX19fO9NXtGhR3rxJOUnJSS4uLlSqVIkSJUpkeiwnJyeuXbvG2bNndRBZ+qnVagICAoiKSsuq+tQeZ0j9UQdFyb5nYffs2UPv3r0pWLAgwcHBrFy5ktjYWLp3755iH7VaTY8ePZgxYwZffvklhw8fxs3Njfbt2xMZGYmtrS1Vq1blwoUL2XYdOU328hZCiHzEwsKC8PBwSpYsmeh9V1dXfv75Z2bNmsWYMWO07xcuXBj4/zOV//bmzRtUKhUFCxbUvhceHo65uXmStplZLLNjxw4GDRrE559/DrxLWseNG8exY8dyfakWPz8/Tp8+Te/evZM9HhMTw99//82zZ894+vQpkZGRtGvXjg4dOqQ4Ztu2bZkwYQLnz5/X+SKk2NhYAgICePz4cZKXv78/T58+JSEhAWPzilSpnPI4+kaFk52FTIh7+3/Hk85e/tub4CAWL16MpaUlFhYWWFhYaP/87/+amppm+jPYtGkTmzdvBt49G/z27VvtQrHVq1cn22fJkiUcPnyYc+fO0aBBA+37nTp10v75m2++wcXFhYYNM1+jMy+QhFIIIfKRatWqcfPmTezt7bXvubq64uLigouLC1OmTEnUvkKFCpiamnLz5s0kY928eRM7O7tEC0tevHhB8eLFk7Q1MDDIcK3IYsWK4e7uzqtXr3jx4gW2traYmpqybNkyevV6tzDk5cuXHywUnhPmz59PgwYNUkx6IiMj8fT0pHTp0lSvXj1NM1pmZmaULFmS48eP07p163TFExERkWyi+P7PL168SNNtWo06Co06Fj1942SPmxaoyttn21E0CYmeo4wJu/XuuGXVFMdWJ0RStKgVVapUITw8nIiICF69ekVERIT2ed339PT0kiSaySWf5ubmyT63Gh4ejqmpKZaWlkDaf/FZtGgRzZo1S5RM/tcnn3zC69evURQl11cf0AVJKIUQIh9xcnJiy5Yt2t1XZs6ciYuLCz/99BPTp09P0t7AwIAuXbqwe/du5syZo/3B++TJE06ePMm4ceO0bWNjYzEwMEjxh6e+vj4JCQkZWkEL7xLL9yvTFy9eTGRkpHY21cvLK8lt+pwWHx/PpUuX6N+/f4ptrK2t+fXXX1GpVERERKT5FmmtWrVYsmRJooRSURTevHmTJEn89yu9jwWoVCpKliyJra1totfhv89zL+Q2ZgVrJ9vPqlQ3gh+vJeT5HgqV+f+zs2+ebMLQpBRm1imXQIoOvU6PVk1o06ZNkmPx8fFEREQQERGhTTbf/zciIoKAgAAiIiKIjIxMch3m5uZJEs/79++nudbkewEBAfj7+9OlSxemTJnCmjVrCA4OplKlSkyaNCnRDlC2trY8efIEW1vbdJ0jL5KEUggh8pHq1aszfvx41Go1CxcuZNq0abRv355OnTolSWbez764urpSt25dOnfunKiweZEiRRItENm1a1eqRZ1r1KjB5cuXE90CPHz4MJGRkYSHv1uE4ePjw86dOwHo2LEjZmZmrFq1Cng3WxoSEsLhw4dZs2YNv/zyC7Vrv0to/vnnH+0t8dzCx8fng89NZnTmysrKCl9fX0aOHJkoYfxvIvUhBgYGlC1bNknCaGtrS7ly5ShTpkySFfbwbl/yof87CCkklAWKt8OyaCueXh+LJiEMI/MKhDzdTvirv7FxWpdiDUoAVeRBevZI/u/S0NCQQoUKUahQ6rfM1Wo1kZGRKSaeL1++xM/PjzNnzqR7C9Bnz54BsH79esqUKcPSpUuxsrJi1apVDBo0iLi4OIYNGwa8+5zu3r0rCaUQQoiPi56eHu3bt+fQoUMcOPCu1p+HhwceHh5J2r6/9Vm5cmVOnTrF999/T69evTAwMKBly5bMnTtXu0vO+x1Gtm9PeWu9gQMHMnv27EQJ5fuE6L0dO3awY8cOAB49ekS5cuVQFIWFCxfy+PFj9PT0qFWrFnv27KFbt27Au5lRb29vZs+enclPR7cuXbqk/Xyyyp9//pnqcVNT02STxfcJY8mSJdO98h6gadOmGCXMQKOJQ08v+cVW5epvI9BnOoF3ZqKOf4OxRSVs62ygUJmUE3+NOhoT5Qb16s1Nd0z/pq+vT4ECBRLVTk2pnZlZyqvVk/O+bFNMTAzu7u7aZLFNmzbUqVOHGTNmaBNKU1PTJLfpP1aSUAohRD4zbNgwPvvsMw4fPoypqWma+jg5OWl3zknO1q1badCgARYWyW/DB1CxYkXevHnDs2fPKF26NAD+/v4fPPfw4cMZPnx4isfd3Nzo1atXrntOzc/P74MzaZlRrFgxnj59SoUKFZJNFm1tbSlSpEiWfC4qlYoxo/ozf/UyzEt/l2wbfQMLytSYR5ka89I8bnTgEqZ8Nzjb/i5NTEyIjo5OV5/3C9UqV66caOZRpVLRrl07fv31V+0uVNHR0akWr/+YSEIphBD5jJWVFf/73//48ccfmT8/+a3z0iMwMJAFCxbw999/f7DtrFmzGDt2LDt37tRJ0vDixQu2bt3KoUOHMj2WrqnV6nQVME+vwoULc+TIERo1apRl50jNN8MGssmtO8HhHTGxtP9whw+IDrtNmYI3GDjgJx1ElzaVKlXi/PnzdO3aNc19KlSokOKs5vtZ/fd/7z4+PtpdnT52UodSCCHyoU6dOmFoaMjSpUszNU5wcDD9+vXj008/5a+//iIwMDDV9g4ODtSvXz/T54V3CzS++eYb5s6di6GhYabH0zVra+t0z36lR1xcHFZWVlk2/oeoVCq2bFpKwnNnEmJfZ2qs+JgXaALH47ZhcbbONNeqVYurV6+mq4+BgQHdunXjzp07iWbYFUXBw8ODChUqaCsO+Pv754vnJ0ESSiGEyLdmz56Nv78/P/zwg3aLw/S4efMmn3/+OQsWLODHH3/E3NyctWvXcv369VT7TZw4kVu3brFmzZqMhk5sbCxfffUVgwYNSnZbyNygTp06vHr1KsvGDw4OTlT+KSeULVuWPdsXE/PoK2Ij/DI0Rky4L3GPB3Fgz59J6qNmtQIFCiRaFAbvFort3LlT+4zx+4ViO3fuJCoqCnhXHcHc3Jz27duzdetW3N3d6dmzJ9evX9c+y+vv76/dbSc/kIRSCCHyKZVKxdy5c6lduzYdO3bk9OnTaapBGBoayowZM5g1axZubm44OjpiZWXF119/TbVq1di7dy8eHh6o1eoUz7t8+XJ8fX355ptvCAsLS1fct27dokuXLvTv358ePXqkq292cnJy+uCMbUap1Wr09fVzxcxstWrVOOaxjoIx3xP5bCmKJvm93P9L0SQQ+WwhRRKmcuLvjVSqVCmLI01e//79cXNz0349cuRIevfuzeDBg4F3C8V69+5N7969tb8gVKhQgTNnzmBnZ8fw4cPp2bMngYGB7N+/X1sbdcWKFYwYMSL7LyiHqJT8tNGkEEKIZL1+/ZpFixZx6dIl2rZtS/369alZsyaWlpYoioK/vz9Xrlzh2LFjvHjxgsGDB9O1a9cksy+KonD58mWOHDmCjY0NvXr1SnbnnPdOnz6Nq6srXbt2ZdCgQYl23fmvu3fvMm/ePMLCwli4cGG2z2ZlRJs2bWjRokWqn4GPjw9xcXHExsayZcsWatasSc2aNYF3jwgkt2XlrVu3sLe3Z+LEiVkVerppNBqW/rGaP1dtI8GoGfoFe2BsURGV6v/PXSmKmtiI+ySE7MIw7hxjRvZj5Iivc3QWLy4ujvbt27Nnzx6dPULw+PFjRo8ezYEDB/LNDKUklEIIIbTi4uI4deoUV65c4caNG0RGRqJSqShbtixOTk40bdqUihUrfnCcx48fs2PHDvT19enTpw+lSpVKsW1CQgIHDhxgw4YNaDQaHBwcqFixIkZGRoSHh3Pjxg38/f0pXLgwFhYWTJo0ifLly+vysrPM3r172bBhA82bN0+xjaurK2/fvk322NSpU7Writ9TFIUtW7Zw5MiRVBPwnKLRaDhy5Cjbdnhw67YvUbEqFA2o9MDMBKpXs6dvn460bt0q1yRbZ8+eZf369dqap5mh0Wjo0aMHc+fOTdP3ysdCEkohhBBZIjQ0lO3bt/Pq1Ss6d+6cph1J4uPj8fHx4eHDh8TFxWFhYUG1atWwsbEB3t1GLFCgAH379s3q8HVi69atjBo1irFjx+qshND169cpUaIEv/zyi07GE++4uLhgZWWVaPen9FIUhQkTJlC1alXtLfP8QhJKIYQQWSYhIYFDhw5x7do16tevT5s2bTJUSPu969evs3fvXkaNGpXlRcMzIzQ0lG+//ZaNGzcCYGNjg7Ozc6auHSAkJAQPDw/++eefXPH85MdEURQmTZqEmZkZU6dOTfcWodHR0YwbNw4HBwfGjh2bRVHmXrIoRwghRJYxMDCga9eudOjQgcuXL7Np06Z0bw/4b9WqVcPCwoLz58/rMErdOnPmDI6Ojtpk0tramp49e7Jv374UFyqlRWhoKPv27WPjxo2STGYBlUrFnDlzqFChAp07d+bGjRtp7nvmzBk6depEly5d8mUyCTJDKYQQIpuk57nK1Jw5c4bTp08zbty4VBe7ZLe4uDhcXFyYPXu2drV8u3btWLt2LaVKlcLNzY3FixfTsWPHdD/7+PDhQ86ePcvmzZtzbDV0fvLs2TNmzJjBixcv+OKLL2jSpAllypTRPvOpKAoPHz7k1KlT7Nixg6pVq/LTTz9l6c5IuZ0klEIIIbJNRp6r/K/o6GgWLFhAo0aN+PTTT3UfZAbcvXuXfv364e3tDbzb0m/OnDmMGTMm0cKT27dvM3z4cMqWLUutWrU+uC1fcHAwFy5coHjx4ixdujRHC5nnRyEhIezatYtLly7x7Nkz7Q44Go2GTz75hIYNG9K9e/d07wf+MZKEUgghRLbSxXOV7u7u3L59m++++y5Hb/8qisLy5cuZMGGCdlecmjVrsmnTJqpWrZpsn4SEBLZu3crKlSsxMDCgZMmSlCxZEjMzMxRFISQkhJcvX/L06VNKlizJhAkTaNiwYXZelhDpJgmlEEKIbJfeepX/9ebNG5YsWULnzp1xcnLKwkhT9r4e5+HDh4F3z+BNnDiRGTNmYGxsnKYxnjx5wuXLl7lw4QLBwcHo6+tTtmxZGjdujJOTU64sCyREciShFEIIkWMy81zltm3bCAoKYtSoUdlez3Dfvn0MHTqUoKAg4N0WhBs2bMg1t+CFyG6yylsIIUSOsbW1ZdiwYVhYWLBu3boP7gP+bw0bNiQoKIgHDx5kYYSJRUREMGzYMLp3765NJvv27cuNGzckmRT5miSUQgghclR69gH/t7Jly1KqVKlsKyF08eJFatWqxerVq4F3cW/evBk3Nze5NS3yPUkohRBC5LiM1KtUqVQ0bNiQR48e8eLFiyyLLSEhgRkzZtC4cWPtbOinn37KjRs3+PLLL7PsvELkJZJQCiGEyBVUKhX16tXjq6++4vXr16xcuZLnz5+n2sfBwQErK6ssm6X08/OjadOmTJ8+HbVajaGhIb///jvHjx/XbgcphJBFOUIIIXKh9NSrPH/+PMeOHcPZ2ZkCBQoQExPDzZs3CQ4ORqVSUbx4cRwcHDAyMkrz+RVFYd26dYwdO1Y7U+rg4ICbmxs1a9bM7OUJ8dGRhFIIIUSulNZ6lbGxscyaNYtXr15x8+ZNAIoWLaotGh4VFcXr168xNDSkb9++fPXVV6kWog4KCmLYsGHs3btX+56zszO//vorpqamur1IIT4SklAKIYTItT5Ur1KtVjNnzhw2b95M3bp1qVatWopF0uPi4rh16xY3b95kxIgRDBs2LEm5IQ8PD77++mvtM5klS5bkr7/+om3btll3kUJ8BCShFEIIkeslV68yMDCQPn36YGtri5OTU5prUarVajw9PYmMjGTbtm0UKFCA6OhoJk2axNKlS7XtPvvsM1auXEnhwoWz6rKE+GhIQimEECJPeP9c5cuXL6lfvz4//PADbdu2pXjx4hka7/Hjx1y8eJFff/2VESNGcOfOHQAsLCxYvHgxgwYNyvaC6ULkVZJQCiGEyDMSEhLYu3cvP/zwA3369KFo0aKZGs/f35+1a9cSFhYGQKNGjdi4cSPly5fXRbhC5BtSNkgIIUSeYWBggJeXF02aNMl0MglQrlw5GjRogLGxMTNnzuT06dOSTAqRATJDKYQQIs8ICAigT58+fP7558nejvb19eXKlSv4+/sTEhKCqakpZcuWpV27dpQtWzbZMTUaDWvXruX8+fNYWFhk9SUI8VGShFIIIUSeMXHiRGJjY6lQoUKyx9etW0dUVBSOjo6UKFGCiIgITp48SUBAACNGjMDe3j7ZflevXqVu3bqMHj06K8MX4qMlCaUQQog8QVEU6tatS//+/VNcLBMeHo6lpWWi92JjY/n5558pWbIko0aNSrbf+2czz5w5o/O4hcgP5BlKIYQQecLjx4+xtrZOdeX1f5NJAGNjY0qUKMHbt29T7GdgYIBGoyEmJkYnsQqR30hCKYQQIk/w8vLK0EKc6OhoAgICKFmyZKrtihYtyo0bNzIanhD5miSUQggh8oRnz55laNHMzp07iYuLo02bNqm2Mzc35+XLlxkNT4h8TRJKIYQQeYJGo0l3oXF3d3e8vLzo0aNHiqu8/3sOIUT6SUIphBAiTyhatChRUVFpbu/h4cHff/9Np06daNq06Qfbx8TEYG1tnZkQhci3JKEUQgiRJzg5OREUFJSmth4eHnh4eNC+ffsP3up+79WrVzg6OmYmRCHyLUkohRBC5An29va8fv36g+2OHDmCh4cHbdu2pX379mkaW6PRkJCQQIECBTIbphD5kkFOByCEEEKkhZ6eHtWqVeP58+eUKlUq2TYnT57k8OHDVK5cGQcHB/z9/RMdL1euXLL9fH1905x8CiGSkoRSCCFEnjF+/HicnZ3p2rVrssdv3boFwN27d7l7926S4wsXLkzynqIoXL16lf379+s0ViHyE9kpRwghRJ7y1VdfUahQIcqXL6+T8a5fv07JkiWZNWuWTsYTIj+ShFIIIUSeEhYWRosWLejZsyfm5uaZGis4OJhjx45x+vRpDA0NdRShEPmPJJRCCJGLhIWF4e3tzZMnT0hISMDS0pIaNWpQsWJF9PRkHeV7165dY8iQIZlKKt++fcuBAwfYvXs3tra2Oo5QiPxFEkohhMhhERERuLm5sWfPHgoWLEitWrUoX748BgYGhISEcP36de7du0fp0qUZOXIkTk5OOjnvo0eP2LFjB15eXkRERKCvr4+iKKjVamxtbWnatCk9evTA1NRUJ+fTtatXrzJo0CCaNm1KxYoV09XXx8eHa9eusW3bthQX6ggh0k4SSiGEyCGKorB9+3ZWrVrF4MGD6dmzJ8bGxim29/X1ZdmyZbx8+ZKFCxdSvHjxDJ332rVr/Pzzz5ibmzNgwADq1auXqFyOoig8fvyYY8eOsWPHDpycnJgyZUqGtj3Mal27duXEiROUL1+eTp06UaJEiVTbP3nyhEuXLlGrVi1+++23XJssC5HXSEIphBA5ICYmhuHDh1OhQgWmTJmSruf3bty4wfjx4/n+++9p3bp1mvvFx8fz888/c//+febPn//B5AveJZdHjx5l9uzZTJ8+nebNm6f5fFnt6NGjtG3bFoCOHTtSuHBh7ty5Q/HixSlSpIj2Vnh4eDivX7/m9evX1K9fn++++y7dM5pCiNRJQimEENksNjaWL774ghEjRtCuXbsMjREdHc2AAQMYNmxYmsaIjY2lf//+dO/enX79+qX7fGFhYYwYMYLOnTvTt2/fjISsU3FxcTg6OnL37l0sLCy4d+8epUqVIiEhgbt373L58mVevnyJnp4epUqVom7duvIcqhBZSBJKIYTIZmPGjKFdu3Z06dIlU+PExMTQo0cPli5dSoUKFVJspygK/fv3p1+/fnTs2DHD51Or1Xz99dd8/vnndO7cOcPj6MK8efOYMGECAHPmzGHixIk5Go8Q+Z38qiaEENnoxIkTKIqCubk5gwcPpnLlypibm1O6dGm6deuGl5eXtq1arWb+/Pm0b9+eMmXKYGZmRpUqVZg8eTIhISGYmJiwbNkynJ2d0Wg0KZ5z7dq11KpVK1PJJIC+vj6rV69m4cKFadoCMasEBgbi4uICQKVKlXB2ds6xWIQQ70hCKYQQ2URRFH799Vd+++03li9fjr+/P87Ozri7u7No0SJevXpFgwYNOHHiBPDutraLiwu2trYsXLgQd3d3hg0bxsqVK2ncuDHR0dF88skntGjRgp07dyZ7zhcvXrBz504cHR0/mMACnD17lqFDh+Lk5ISxsTEqlSrR9oVGRkY5PiP4/fffExERAcDixYsxMjLKsViEEO/ILW8hhMgmp0+f5vjx48yYMYNXr15RrFixRMcjIiKws7OjWrVqHDt2DLVaTUhICIULF07UbufOnfTu3ZuNGzfSv39/wsLC6Nu3LwcPHkxyTldXV5o1a8ayZcsIDg6md+/eODg48Pr1a+bNm8eVK1c4cuQILVu21LZft24dtWrVIiQkhFOnTvHo0aMkpXUGDBjArFmzsLGx0e2H9AHnzp2jSZMmAPTo0YPdu3dn6/mFEMmTGUohhMgm69evZ+jQoQBJkkkACwsLHBwcCAgIAN7dYv5vMglQr149AG27AgUKUKJEiUQziQAJCQmcPXuWTz/9lD/++IMTJ04wcuRImjdvTq9evTh69CiFCxfml19+0faZOnUq/v7+7Nmzh06dOqV4Ld988w0rV65M3weQSWq1mjFjxgBgYmLC/Pnzs/X8QoiUSUIphBDZ5OXLl6nO6IWGhuLt7U3VqlVTHef9LfF/t2vatCmXLl1K1O7WrVvUrVsXlUqVpgQWSPMq6MaNG3P16tU0tdWVlStXcu3aNQAmT54sBcmFyEUkoRRCiGwQFRX1wSLao0ePJjIykh9//DHFNs+ePWPy5MnUqVMn0UprJycnvL29E7X18vJKdVedtCawyVGpVBgaGhIXF5fuvhkRFBSk/VzKlSvHpEmTsuW8Qoi0kYRSCCGywZs3byhatGiKx6dOnYqbmxsLFixIMQl88+YNHTt2RFEUtm3blmg2sXjx4klWXvv6+lKlSpUUz5mWBDY15cuX5/Hjxxnqm14//fQTb9++BWDBggWyw40QuYxBTgcghBD5gUqlSvGYq6srP//8M7NmzdI+I/hfb9++pU2bNjx79ky71eC/KYqS5ByxsbGYmJgkO977BHbJkiUZ3hvcxMSEmJiYDPVNDy8vL+3zmu3ataNbt25Zfk4hRPpIQimEENnA2tqaV69eJXnf1dUVFxcXXFxcmDJlSrJ93759S+vWrXn06BHHjx+nRo0aSdq8ePEiyQyoiYkJUVFRyZ7zQwlsWqTlNn5maTQavv32WxRFwdDQkEWLFqWanAshcobc8hZCiGxgampKdHR0ovdmzpyJi4sLP/30E9OnT0+23/tk8uHDh/z999/UqlUr2XZeXl7Url070XsODg7cvn070XtpSWDT6uHDh1m+MGbjxo2cP38egHHjxlGpUqUsPZ8QImMkoRRCiGxSqlQpbWmfefPmMW3aNNq3b0+nTp24cOFCohe8K2zerl07rl69iqurKwkJCYna+Pn5acc+c+aMtpzQe05OTly+fFn7dVoS2LTSaDSo1WoMDLLuRldoaCjff/89ACVLluSnn37KsnMJITJHbnkLIUQ2+frrr1m1ahWzZs3iwIEDAHh4eODh4ZGkraIovHz5UpsQJre94MCBA/nrr78ICQnhwYMH6OvrJzpepUoVrl27hlqtZuHChUkS2H9r0KABAK9fv+b06dMA3Lx5E4DDhw9TtGhRihYtSvPmzQE4fvw4DRs2zMzH8UGurq68fPkSgLlz52JpaZml5xNCZJzslCOEENlEURTatWvHjh07sLKy0tm4v/zyC48ePaJMmTLY2NjQqFEj7O3tUalUzJkzh6pVq/L7779rE8WUYgM4deoULVq0SLZN8+bNOXXqFAB9+vRh0aJFlChRQmfX8W+3b9/G0dERtVpN06ZNOX36tDw7KUQuJgmlEEJko3/++Qc3NzdWrFihk/Hu37/PhAkT2L17N76+vpw/f56AgAAKFy5Mw4YNKVu2LH379sXd3V1ne17/888/bN26lWXLlulkvP9SFIXWrVtz4sQJ9PT08Pb2xtHRMUvOJYTQDUkohRAim40bN44mTZrQs2fPTI0TFRVFjx49WLFiRaLFMQEBAXh6enL37l3Mzc2Jjo4mNjaWn3/+OZORQ2RkJF26dGHv3r0UKFAg0+Ml5/1e5fCuVubSpUuz5DxCCN2RhFIIIbJZfHw8/fr1o3///nTt2jVDY0RGRtKvXz/Gjh1Ly5Ytk20THBzM+fPnuX79OseOHWPgwIEMGjQow3HHxMTQtWtXBgwYwIABAzI8TmoiIyOpUqUKAQEBFClShHv37mFtbZ0l5xJC6I6s8hZCiGxmaGjIpk2bOHDgAD/88EO6i4NfvnyZrl278r///S/FZBKgcOHCdO7cme+++46ffvqJdevW8euvv6LRaNId84sXL+jduzf169fn0aNH3LhxI91jpMXs2bO1e4v/+uuvkkwKkUfIDKUQQuSgvXv3snjxYvr3788XX3yBmZlZim1v3rzJH3/8QVRUFAsWLKBw4cLpOldsbCyurq6cP3+e5cuXU7ly5Q/2UavVbN68mY0bNzJ37lyqVavGgQMHuHbtGp06daJOnTrpiiE1fn5+ODg4EBcXR506dbh48WKi7SWFELmXJJRCCJHDoqOj2bFjB9u3b8fY2JgaNWpQvnx5DA0NCQkJ4dq1a/j7+2NnZ8fIkSOpWrVqps7n5+fHlClTCA4OZuDAgTRq1Ijy5ctrV1FHR0dz48YNjh07xokTJ+jSpQujR4/G0NAQeLdoxsPDg0uXLtGmTRsaNWqU6c8AoGvXrtpyShcuXKB+/fo6GVcIkfUkoRRCiFwkOjqa69ev8+TJExISErC0tMTR0ZGyZcvqvGzOmzdv+Ouvvzh9+jTPnz/HyMgIExMTLC0tqVmzJk2aNKFly5bJzhIqisKJEyc4e/YszZs3p3nz5snGFx4ezpYtWzh16hRPnjzh3z9yihYtStOmTenfvz9eXl506tQJeFevc+3atTq9ViFE1pKEUgghRKIFPCqVipo1a9KwYUMKFSqUar+zZ89y/PhxGjRoQNu2bbVJZXBwMJMmTeLWrVtUqVKFTz75BCsrq0RJZ2RkJP7+/ty+fZsHDx7w6tUrrKysuHfvHsWLF8/S6xVC6JYklEIIIbQiIyO5fPkyly5dIiYmhipVqtCoUSNKly6dYp9Lly5x+PBhateuTadOndi7dy+zZs2iefPm2Nrapum8QUFB/PXXX9SrV4+dO3dKEXMh8hhJKIUQQiQRHx/PtWvXuHDhAm/evMHW1paGDRtqd+D5r2vXrrF//378/Px49uwZ7du31z5zmVYajQZPT0/MzMxYt26dLMgRIg+RhFIIIUSKNBoN9+7dw9PTk6dPn1KkSBEaNmxIjRo1MDAwSNR2+vTpnD59mu7du2dqhtHLywsLC4ss24lHCKF7klAKIYRIk//uwFOvXj3q1KmDmZkZ9+/fZ8CAAXzxxRc6uV19+PBhnJ2d6dChgw4iF0JkNUkohRBCpMt/F/A4Ojry22+/0aJFi2QLkT99+hR3d3eeP39OZGQkhoaG2hXeKdWxjI2NZdu2bZw5cybV2pxCiNxBEkohhBAZ8n4Bz86dO7l9+zbdu3dPtt39+/e5evUq5cuXx8rKiri4OK5cucLVq1fp2LEjbdu2Tbafl5cXjRs35ptvvsnCqxBC6IIklEIIITKld+/eVK5cmSJFiqSr34IFCwgNDcXFxSXZ43FxcezZswdPT08dRCmEyEqyhE4IIUSGKYpCQEBAupNJAHNzc/T19VM8bmRkhJGREcHBwZkJUQiRDSShFEIIkWHPnz/HwsIiTW01Gg1qtZqIiAjOnj3L3bt3admyZap9ihcvjpeXly5CFUJkIYMPNxFCCCGSd+/evQ/upvPezp07tbev9fX1+eyzz2jcuHGqfaytrbl161aKz1kKIXIHSSiFEEJkWGxsbJJ6lClp06YNDRo0ICIiglu3brFr1y7i4uJSnaU0MDAgOjpaV+EKIbKIJJRCCCEyzNTUlISEhDS1LVSokHY208HBAYCDBw9Sr169FG+bx8fHY25urptghRBZRp6hFEIIkWFVqlTJ8KIZW1tbNBoNQUFBKbYJDg6mRo0aGQ1PCJFNJKEUQgiRYcWLFycqKipDfe/fv49KpUp1hfjLly+pXbt2RsMTQmQTueUthBAiU+zt7Xnx4gUlSpRI9vi2bdswMTHBxsYGS0tLIiMjuXbtGlevXqVly5Yp3u6OjY1FURQKFiyYhdELIXRBCpsLIYTIlGvXrvH999/TsWPHZI9fvHiRixcv8vLlS6KjozE2NqZUqVI0bNgwxa0XAS5dukS7du346quvsip0IYSOSEIphBAi09q3b0/t2rUpXry4TsaLiopi9+7dnD17FmNjY52MKYTIOvIMpRBCiExbtWoVf//9N2q1WifjHT16lPnz50syKUQeIQmlEEKITCtbtizjx4/nwIEDaDSaTI119uxZGjRoQNOmTXUUnRAiq0lCKYQQQie+/PJL+vXrx65duzJUjFytVnP06FHKli3LrFmzsiBCIURWkWcohRBC6NTp06cZP348tWvXxsHBAZVK9cE+T5484dSpU3z33XcMGDAgG6IUQuiSJJRCCCF0LjIykiZNmvD06VPq1atH1apVKVGiBPr6+gBoNBpev37NkydP8PX1pUaNGsyaNUtni3qEENlL6lAKIYTQudjYWO7cuUNsbCzh4eHEx8dz6NAh4uPjAdDT06NSpUq0b9+eP//8M8ValEKIvEESSiGEEDq3ceNGYmNjAXBxcaFly5Y5HJEQIivJLW8hhBA6pSgK1apVw8fHBzs7O+7du4eenqwBFeJjJt/hQgghdOr8+fP4+PgAMGzYMEkmhcgH5LtcCCGETq1cuRIAAwMDBg4cmMPRCCGygySUQgghdCYkJITt27cD0L17d1m1LUQ+IQmlEEIInXFzc9MWNR82bFgORyOEyC6yKEcIIYROKIpCzZo1uXHjBuXKlcPPz0+enxQin5DvdCGEEDpx+fJlbty4AchiHCHyG/luF0IIoRPvF+Po6+vz9ddf53A0QojsJAmlEEKITAsLC2PLli0AdOnShZIlS+ZwREKI7CQJpRBCiEzbsmULUVFRAAwfPjyHoxFCZDdZlCOEECLTnJyc8Pb2xsbGhocPH6Kvr5/TIQkhspHMUAohhMgULy8vvL29ARgyZIgkk0LkQ5JQCiGEyJRVq1YBoKenx+DBg3M4GiFETpCEUgghRIZFRESwefNmADp27EiZMmVyOCIhRE6QhFIIIUSGbdu2jfDwcEAW4wiRn8miHCGEEBlWv359Ll26RKlSpXj8+DEGBgY5HZIQIgfIDKUQQogMuX79OpcuXQLeLcaRZFKI/EsSSiGEEBnyfjGOSqViyJAhORyNECInyS1vIYQQ6RYVFUWpUqUIDQ2lffv2HD58OKdDEkLkIJmhFEIIkW47duwgNDQUkMU4QgiZoRRCCJEBjRs3xtPTk+LFixMQEIChoWFOhySEyEHyBLUQQohENBoNfx8+zJHNm7l/6xaqyEhUioKipwcWFpT85BM8PT0BGDx4sCSTQghJKIUQQryj0Wj4c+5cdq1YQduwMEZER1NRpUJPpdK2SXjxgju+vhgaGnIwIYFiVlYoioLqX22EEPmP3PIWQgjBw4cP+bZXL7o8fsywuDj005AgxikKc01Nue7gwOJt2yhevHg2RCqEyI0koRRCiHzumrc3E3v0YFVwMOX00r9W86ZajXOJEqw5doxPPvkkCyIUQuR2klAKIUQ+9vDhQ4Z/+inb377FOhO3rZ9oNAwqUYIdFy9SuHBhHUYohMgLpGyQEELkUxqNhm979WLVmzeZSiYBbPT0mPPiBeP699dRdEKIvEQSSiGEyKdWzp9P58eP+SSNt7lXx8WhCg/HIjw82eN19PQo6+3Nwd27dRmmECIPkFveQgiRD2k0GtrY2/P3y5dpWoDzTKOhamQk5ioVoYpChKVlsu3CFIW+dnYcvHpV1yELIXIxmaEUQoh86KiHB63DwtKUTAKMiImhmb4+bfT1U21XQKXC5uVLbt++rYswhRB5hNShFCKXCw4OxsvLi+vXrxMaGoqBgQFlypShTp06VK1aVYpKiwz5e8sWhkZGwgcSRIBN8fGcVqvxMTfnp9jYD7bvHhLC0T17qFq1qi5CFULkAZJQCpELqdVqDh48yLp16zAzM6NevXrUrVuXQoUKoVar8ff3Z8eOHUyePJmqVasycuRI7OzscjpskYfcu3kT+zQ8O/lKo+G72FhmGxtTJo3PWjrp67Ph5En46afMhimEyCMkoRQil7l37x7fffcdLVu2ZO3atVhbWydpU7t2bT777DMArly5wo8//kiVKlWYMmUKRkZG2R2yyINUkZFput09KjaWSioVI9MxE15YpSLk1avMhCeEyGMkoRQiF9m2bRtubm6sWLECGxubNPWpU6cO27ZtY/v27XTp0oUNGzbIjiXig1QazQfb7IqP50BCAlfNzNK9taKerPcUIl+RhFKIXGLr1q0cOXKE3bt3Y2CQ/m/Nzz//nOrVq9O3b1+2bt1K0aJFsyBK8bFQPvDsZISiMDo2lm8NDSmlp0fI/yWIcf93PERRMATMU0g0NbK3txD5iqzyFiIX8PHxYevWraxatSpDyeR7VapUYcGCBQwfPhypCCZSo2dpSXwq/0aCFIWXisK8+HgKRURoX1sSEogECkVE0C8mJtm+zzUaiqdxhl0I8XGQhFKIHJaQkMC4cePo27cvw4cPp3Llypibm1O6dGm6deuGl5dXkj7x8fHMnz+f6tWrY2pqSsGCBWnUqBGenp7UqFGDxo0bs27duhy4GpFXONSuze1UbnuXUKk4aWqa5NVOXx8T4KSpKT+n8Lyul0aDU8uWWRS5ECI3ksLmQuSwrVu3EhQUxOnTpwkODqZ37944ODjw+vVr5s2bx5UrVzhy5Agt/+8HtFqtplu3bpw9e5ZJkybRqFEjIiMj8fLyomHDhrRp0wa1Wk27du1wd3eXRToiWZ7nznGwZ09+iYpKV79B0dHsTEhIsbA5wEArK1zPnKFcuXKZjFIIkVdIQilEDuvcuTNbtmwhOjqaYsWKJToWERGBnZ0d1apV49ixYwAsXLiQ8ePHc+7cORo0aJDiuMuXL8fa2po+ffpkafwib1IUhXYODhx4+hTjdDzv+KGE8rVGw8gaNdjp6amrUIUQeYDc8hYiBwUEBFCiRAksLS2TJJMAFhYWODg4EBAQoH1v0aJFNGvWLNVkEqB///5s375d5zGLj4NKpWLYlCnMNjNLV7+/TE1TnZ380dKScbNnZzY8IUQeIwmlEDno8uXLNGzYMMXjoaGheHt7a3ccCQgIwN/fn+rVqzNlyhSKFy+OgYEBVatWZf369Yn6WlpaEh0dnaXxi7ytV//+3K1WjatpKCGUFh6Aedu2NG7WTCfjCSHyDkkohchB165do1atWikeHz16NJGRkfz4448APHv2DID169ezb98+li5diru7Ow4ODgwaNIhVq1Yl6l+iRAlevnyZdRcg8jSVSsXirVv5X4kSPMpkUnlVo2G+rS2zli/XUXRCiLxEEkohclBISEiyO+EATJ06FTc3NxYsWICTkxMAmv/7oR8TE4O7uzu9e/embdu2bN++ndq1azNjxoxEYxQuXJi3b99m7UWIPK1o0aL8dfw4Q0qV4kIGk8rDKhWTy5dny6lTmKXzFroQ4uMgCaUQOUhPTy/ZepGurq78/PPPzJo1izFjxmjfL1y4MACVK1fG1tZW+75KpaJdu3Y8ffqUV//a8k6j0aD/gQLWQtja2rLr4kVWfPopYxRFW8T8Q15pNAwzN+dojx7suXBB++9TCJH/yE45QuSgkiVL8uTJEz755BPte66urri4uODi4sKUKVMSta9QoUKKM0DvE1M9vf//e+LTp08pUaJEFkSeMwIDA7l79y6xsbEYGxtTqVIlSpUqldNhfRQKFSrEgk2bKFuyJHtjY2lkYsJgwElPj6L/929KURSeKwpXNBp2FS5MtI0N//vtNxo2aZKzwQshcpwklELkoDp16uDl5UXz5s0BmDlzJi4uLvz0009Mnz49SXsDAwO6devGzp078ff319b5UxQFDw8PKlSoQJEiRbTtIyIisExlRW5e4O3tzfLly3n69CklS5bEwcEBMzMzoqOj2bRpE8+fP6d06dKMGDGCOnXq5HS4edqff/5JREwMEUDrhQu5/+YNW06c4M2LF6gUBUVPjxI2NtRp1YqZPXokmiUXQuRvUodSiBwUGhrK4MGD2bVrF/PmzWPChAm0b98+2WTyfZkgPz8/6tSpQ/HixXFxcaFAgQKsXr2avXv3sn37dnr16gXAo0eP+Pnnn1mzZk22XpOuPH78mO+++w5bW1tGjhxJpUqVUmx7//59/vzzTx48eMCCBQsoX758Nkb6cYiNjaVcuXK8ePGCatWqcePGDVSyH7cQIo1khlKIHGRlZYW5uTn+/v4cOHAAAA8PDzw8PJK0ff+7X4UKFThz5gyTJ09m+PDhxMfHU7NmTfbv30/nzp217VesWMHQoUOz50J0bPPmzbi5ubFw4UIqVqz4wfYVK1Zk3rx5PHz4kLFjx/L555/z1VdfZUOkH49Nmzbx4sULACZOnCjJpBAiXWSGUogcdvHiRTZu3MjSpUt1NmZwcDD9+vXj8OHDeS4xWLVqFd7e3ixdujRDC4o0Gg3fffcdlStXZtSoUVkQ4cdHo9Hg4ODAvXv3KFOmDH5+frJlpxAiXWSVtxA5rH79+sTHx3PmzBmdjTlu3Dh+/vnnPJdMHjt2jHPnzvHHH39keHW6np4eixYt4urVq8nO9IqkDh48yL179wD47rvvJJkUQqSbzFAKkQuEhobSvXt3Nm/eTMmSJTM11po1a/D392fmzJk6ii57hIWF0b17d/bv34+FhUWmx4uKiqJz587s3r2bggULZj7Aj1iTJk04d+4cVlZWPHnyhAIFCuR0SEKIPEZmKIXIBaysrFixYgX9+vVLtG93eq1btw5PT09cXV11GF32mD17NlOmTOHSpUsMHjyYypUrY25uTunSpenWrRteXl6J2g8aNAiVSpXkVblyZQDMzMyYPn06P//8c05cTp7h6enJuXPnABgxYoQkk0KIDJGEUohcwt7enlWrVjF48GC2bNmSbMHzlLx9+5bBgwdz//59Vq1alagWZV4QExPDlStXaNWqFcuXL8ff3x9nZ2fc3d1ZtGgRr169okGDBpw4cSJRP1NTU86fP5/otW3bNu3x5s2bc+vWLaKiorL7kvKM33//HQAjIyPGjh2bw9EIIfIqueUtRA4ICwvj6tWrPH36lISEBKysrHB0dKRcuXIkJCQwe/ZsLl68yNdff03Xrl0xNDRMdpznz5+zcuVKzp07x4wZM2jYsGE2X4lu7Nixg+DgYEaMGMGrV68oVqxYouMRERHY2dlRrVo1jh07Brybody5cycRERGpjr1mzRpMTU3p27dvlsWfV927d48qVaqgKAqDBw/OsyWmhBA5T8oGCZFNwsPDcXNzY9++fVhZWVGrVi3KlSuHsbExr1+/Zu7cufj5+WFnZ8fIkSP59ttv2bBhA927d0dfXx97e3usra1JSEjA39+fp0+fEhkZyahRo5g+fXqeW4Dzb+fPn2f48OEASZJJAAsLCxwcHDL0OMCnn37KwoULJaFMxrx587Qz4RMmTMjhaIQQeZkklEJkMUVRcHNz46+//mLo0KHs27cv1VW0N27cYNGiRcTHxzNv3jzGjh1LQkICvr6+hIaGYmBgQNmyZSlevDiLFi2iYMGCeTqZhHeFye3t7VM8Hhoaire3Ny1btkz0fnR0NCVKlOD169eULFmS7t27M2PGDKytrbVtypcvz6NHj7Is9rzqxYsXbNiwAYAuXbpQpUqVHI5ICJGXSUIpRBaKiopi6NChVKtWDQ8PDwwMPvwtV6NGDVauXMnFixfp3bs3M2fOpFGjRjg4OCRpW6VKFW7dukXHjh3zdFKpUqlSfe5z9OjRREZG8uOPP2rfc3R0xNHRkWrVqgFw+vRpFixYwPHjx7l8+bJ2pfiHxs6vlixZQmxsLACTJk3K4WiEEHmdJJRCZJHo6Gj69OnD//73P1q0aJHu/vXr12fv3r18+eWXfP/99zRt2jRJmypVqnDhwgWePn1K2bJldRF2rjN16lTc3NxYsmQJTk5O2vfHjRuXqF2bNm2oVasWvXr1YtWqVUmOi/8vIiKCZcuWAe+29GzcuHEORySEyOvk13aR62g0GqKiooiNjU3XSufc5ttvv+Xbb7/NUDL5nqWlJVu3bsXV1ZVnz54lOV62bFksLCzw8fHJTKg5Tq1WJ/t37erqys8//8ysWbMYM2bMB8fp0aMH5ubmXLhwIdH7Go1GZ7F+DFavXk1ISAjwbnYyL89uCyFyB5mhFDlOURTOnDnDli1bePz4Mfr6+piZmaFWq4mKisLQ0JB69eoxaNAgSpcundPhpsmhQ4coWLAgBgYGDB48GE9PTwICAihYsCB16tRh2rRpiWbbUvuBXqlSJXbv3o2zszM7duxI1FalUlGpUiXu3r1L27Zt82xiYGNjw5MnT7C1tdW+5+rqiouLCy4uLkyZMiXNYymKkugW9/PnzylRooRO483L4uPjWbBgAfBuD/SuXbvmcERCiI+BJJQiRx06dIiFCxfSoEEDJkyYQPny5ZMkRXFxcZw5c4YJEyZgYGDA7Nmzc3ViqdFomD9/PgcPHuSrr74iODgYZ2dnHBwceP36NfPmzaNBgwYcOXJEu8jk/PnzSca5ePEi3333HT169MDBwYFq1apx5MgR2rdvn6hdlSpV8PLy4sWLF5neZSen1KtXj3PnzmkTypkzZ+Li4sJPP/3E9OnT0zzOzp07iYqKokGDBtr3PD09qVu3rs5jzqu2b9/OkydPgHcruzO6xaUQQiSiCJEDwsLClK+//lqZMGGCEhUVleZ+Pj4+Srt27ZQNGzZkYXSZ4+7ursyePVtRFEV5+fJlkuPh4eFK8eLFlVatWqU6zqBBgxSVSqXcv39fO1avXr2StEtISFBmz56tHD9+XAfR54y3b98qXbp0URRFUebOnasASvv27ZXz588neSmKovj7+yuNGjVSFi9erLi7uyuHDx9WJk+erJiYmChVq1ZVIiIitGN3795dCQoKypHrym00Go3i6OioAEqxYsWU6OhonYwZERGhhISEKHFxcTqIUgiRF8kMpch2ISEh9OnThylTptC8efN09a1SpQoHDx7khx9+4Ndff+WHH37IoigzbuPGjSxatAjIeE3F8PBwduzYQfPmzbGzs9OOZWhomKTwt76+PpUqVeLOnTtJyurkFQULFqRkyZLcuHGDAwcOAODh4YGHh0eStoqiUKBAAYoXL878+fN5+fIlarUaW1tbxo4dy5QpUzA3Nwfgzp07FCpUiMKFC2fr9eRWR48e5fr16wCMHTsWExOTDI3z8uVL1q1bx4ULF4iPj8fS0hIDAwMiIiKIi4ujWLFi9O7dm/bt28sMqBD5RU5ntCJ/iY+PVzp16qRcunQpU+NoNBpl8uTJyqpVq3QUme506NAh1eMhISGKlZWV0qNHjxTbrFy5UgGUTZs2JXp/yZIlyqFDh5K0v3v3ruLi4qK8evUqY0HnAk+fPlXatm2rxMfH62S8hIQEpUOHDsrjx491Mt7HoHXr1gqgmJubK8HBwenu//LlS2Xw4MFK7969FQ8PDyU2NjbZdoGBgcqcOXOUVq1aKevXr1c0Gk1mQxdC5HKyyltkq3nz5tG7d+9MP9OmUqmYNWsW+/bt4+HDhzqKLvOCg4M/OBuWXE3F/1qzZg0FCxakZ8+eid53cnLCy8srSfvy5ctjaGjI3bt3MxZ4LlC6dGm+/PJLZs2apZPx5syZQ/fu3bGxsdHJeHmdt7e3dtvKoUOHJir+nhZ79uyhX79+jBkzhu3bt9OuXbsUC/SXKFGCiRMn4uHhQWBgIL179+bVq1eZvgYhRO4lCaXINo8ePeLixYuUKVOGwYMHU7lyZczNzSldujTdunVLkigpisKqVatwcnKiQIECFC5cmObNm3Po0CEA9PT0WLx4ca7aMi4oKCjVFcXvayouWLAg0Srvf7t9+zYXL16kX79+SW5JlihRgqCgoCR9DA0Nsbe3586dO5m7gBw2cOBAwsLCWLx4cabGWb58OYGBgQwbNkxHkeV9c+fOBd49IpHeGp1//vknR44c4dChQ9SqVSvN/QwMDPj++++ZPn06X375ZbKlr4QQHwdJKEW2WbZsGZMnT+bPP//E398fZ2dn3N3dWbRoEa9evaJBgwacOHFC23769OkMHz6cevXqsWvXLv766y+MjY3p3Lkzu3fvBuCTTz6hRIkS3L59O6cuK83SWlNxzZo1wLtZpPSoXLkygYGB2vqCeZFKpWLu3Lm8fv2aESNGEBYWlq7+ERERjBkzhsePH7No0aI8W0ZJ1/z9/dm+fTsAffr0SVSe6UN27drFtWvXWL58eapbhqamevXqrFixgoEDBxIZGZmhMYQQuVxO33MX+UNMTIzStm1bRaPRpHnlc+nSpZUmTZokahcdHa1YWVkpXbt21b5348YNZcyYMVkXfDoEBQUp/fv3T/K+i4uLAiguLi6p9o+NjVWKFCmiODk5JXvc09NTmTlzZrLHYmJilJkzZyqenp7pDzwXOnXqlNKiRQtl1apViVZtJycyMlJZs2aN0qJFC+XYsWPZFGHeMXbsWAVQAOXq1atp7vfixQulTZs2Olu97eHhoTg7O+tkLCFE7iIzlCJb3Lhxg/r166NSqdK88tnQ0BArK6tE7UxMTLSv96pXr46fn1/WBZ8OhQsXJjg4ONF76ampuH//foKCghgyZEiyx728vKhdu3ayx4yNjalQoUKefo7y35o3b467uzvGxsb06tWLfv36MX/+fI4cOcI///zDkSNHmD9/Pt26ddOuJj506BCtWrXK6dBzleDgYFavXg28256yZs2aae77/fffM3fuXM6cOZOmx1QWL15MgwYNKFKkCMbGxtjY2PDFF19o7yC0a9eOqKgoLl++rLPrE0LkDlI2SGQLLy8v6tSpk+Lx0NBQvL29E5W9cXZ2ZsKECaxZs4bPPvuMmJgYfv/9d0JDQxk7dmyi/mZmZkRGRmrLxeSkggULakv7zJs3j2nTptG+fXs6deqUZEvAfxfghne3u01NTenbt2+yY3t6evL555+neO4qVaqwb98+IiIisLCwyPzF5DATExMGDBjAgAEDePPmDd7e3ty9e5eYmBhMTEyoXr06VlZWqFQqBg4cmNPh5krLly8nKioKgIkTJ6a537Nnz4iLi6NGjRr07t07TQX6g4OD6dChA46OjhQqVIiHDx8ye/Zs6tevj5eXF5UqVWLKlClMnz6d9evXZ8n1CiFySE5PkYr8YeLEicrt27dTPN6vXz/FwMBAuXLlSqL3//zzT8XY2Fh7u87a2lo5evRosuP7+PjoPO6McHd3V3799VdFURSlefPm2tiTe/3bkydPFD09PeWrr75KdtyUCpv/W2RkpOLq6qpcvnxZNxeTB3h5eSmurq7pKpCfX0RFRSlFixZVAKVmzZrpKt/j6uqqnDhxQlGUzBXo9/HxUQBl6tSp2vd69uwpxeaF+MjILW+RLRISEjA0NEz2WEorn9etW4ezszNjxozh2LFjuLu707ZtW7p168aRI0cSjWFoaEh8fHyWXkNatWvXjqNHjxIdHc2pU6dQFCXF17+VLVsWtVqd4szNsmXLUrwV/p6ZmRmffPJJnl/tnR52dnYoipJrHnvITTZs2MDr16+Bd7OT6VmkdOnSJZo1awZkvEA/QNGiRYF3K77fa9myZbLbjQoh8i5JKEW2MDc3Jzw8PMn7Ka18fvv2LaNHj2bo0KHMnTuXVq1a0aFDB7Zs2ULdunUZMWJEonHCwsJyzS1ePT09xo8fz9SpU3U25u3bt7l9+zbt2rX7YNvKlSvj7+9PdHS0zs6fmxUoUIBixYrx4MGDnA4lx8TFxREWFpbo71ytVmtLBdnY2NC7d+80j6fRaFAUJdVdbt4/plK1atUkx9RqNbGxsdy9e5ehQ4dSrFgxvv76a+3xlOqpCiHyLnmGUmSLatWqcePGjUQLSlxdXXFxccHFxYUpU6Ykan/v3j2io6OTLYBep04dTp8+neg5wUePHqWrFEpW69ixI7t37+bvv/+mbdu2mRorIiKCsWPHsnHjxjTNMFWuXBl3d3fu3buXrgUYeZmdnR3Xr19HUZR8USooPj6eXbt2sXXrVm3BcENDQ9RqNQkJCRQoUIBSpUppk+z//e9/Kd4hSM7bt28pUqRIqm1SK9Bvbm5ObGwsAPb29pw6dYqyZctqj9vZ2fHo0aM0xyOEyP0koRTZok6dOsyfP59BgwYBH175XKpUKQAuXLiQaLGFoihcuHCBQoUKaRfgJCQkfHA2JScsWbKEPn36YGBgkOE9tsPDw+nSpQuNGzf+4A/49ywtLSlbtix3797NVwmlp6cnL168oGTJkjkdTpbRaDQsXLiQzZs3Y2dnR926dZOdmY+Li8PHx4eyZcsSERFBr1690nWe1B5Rgf//mMqSJUuSLdDv6elJXFwcfn5+LFiwgBYtWnD8+HHtbGZuekRFCKEbklCKbFG+fHl8fX2JjY1l6dKlH1z5bGNjw2effcbKlSsxNjamY8eOxMbGsn79es6dO8fMmTO1M1EHDhygTZs2OXFZqTI1NWXr1q0MGzaM8+fPM2nSpHTNEl24cIEff/wRZ2dn7t69i5ubG3379sXY2PiDfatUqcLx48eJi4vLcDHqvMTGxgYjIyPu37//0SaUjx8/ZuDAgZQuXZp+/fqhp5fyE0tGRkbUrFmTmjVr8uDBA7p27crs2bPT/H2S0iMqkLYC/e/vRDRo0ICuXbtiZ2fHlClT2LdvH/DuF6Xc8oiKEEI3JKEU2UKlUtGrVy927NjBgQMHAPDw8MDDwyNJ2/eLVdzc3Fi6dCkbN25k7dq12u0FN23alKiszrp163JtCRIzMzM2bdrEli1b6NChA4MHD6Znz56pJoXXr1/njz/+ICEhgZ07d1KoUCECAgLYtGkTbm5u9OvX74NJZeXKlfn777+5f/9+ss+4fWz09fUpX748Dx480C4k+Zj4+vry5Zdf0rFjxzTPVL9nZ2eHjY0Nrq6uBAUF8eWXX36wj4WFRYrPPKf0mEpKLC0tqVy5Mr6+vtr3rl+/TvXq1dN+EUKIXE+l/HepqRBZJCoqis6dO7N3714KFCigkzH37t3LpUuX+OWXX3QyXlaKiIhg8+bN7NmzB0tLS2rVqkW5cuUwMDDgzZs3XLt2jUePHmFvb8+IESNwcHBI1P/p06ds2rSJokWL0r9//w8mlStWrKBw4cLpvt2ZV3l5eXHo0CEmTpyIqalpToejM0FBQbRv356uXbsmKfSfHmq1ml27djFjxowPFn+Pj4/ns88+Y+PGjRQsWBB495jKtGnT+Omnn5g5c2a64re3t6dx48baXyZdXFxo164dDRs2zPD1CCFyF0koRbY6ffo0mzZtYtWqVZke6/Xr13zxxRfa3VTykoiICK5evUpAQABqtRorKyscHR2xsbFJdVHJs2fP2LhxI0WLFqVfv36Jdgz6r3/++Ydz584xceLERCVbPlahoaEsXLiQXr16fVSzsj179qRcuXLY2NhkeqzY2Fi2b9/OyZMnsbS0THRMrVbz8OFDbt26xd27d7l27RqNGzdm7NixzJs3jwkTJtC+fftkn3lu0KABoaGhtGnThr59+1KxYkVMTU3x9fVl0aJFPHnyhNOnT1OnTh00Gg1t2rThyJEj+eLfpRD5hSSUItvNmjULRVH46aefMjzG27dv6dOnD/Pnz6datWo6jC73e/78ORs3bqRw4cL0798/xaTy9evXLFu2jC+++IJKlSplc5Q5Y9myZdptAT8GBw8eZPny5SlWCrh//z5//PFHsse+++47ypUrl+T9Bw8ekJCQwPLly1EUhcePH3Pr1i18fHyIjo6mSJEiVKtWDXt7ewYMGMDff/9Ny5YtOX36dIpxKopCbGwsY8aM4dy5cwQEBBATE0OJEiX49NNP+eGHH7Qz7keOHOHKlSvJrg4XQuRdklCKbKcoCjNnzuTVq1fMmTMHMzOzdPW/ffs2X3/9NfPnz6dJkyZZFGXuFhgYyIYNG7C2tqZ///4p3uL9448/KFOmzEeTYH3I33//zc2bNxk3bhxPnjwhNDQUfX19ypQpo711m5e0bt2aVq1apfg98j6h7NSpExUrVkx0rGTJkinO3G/cuJEff/yRR48eER4ejpWVFdWqVaNatWoUL15cO0u+ZMkSDAwMGDlypE6uJyYmhk6dOrFjxw6sra11MqYQIneQ+w25nKIo/PPPP+zZs4cHDx6gUqlQqVQoioJGo6FYsWK0atWKXr16pXr7MzdRqVRMmzYNd3d3OnfuzMSJE2nXrl2qq1bh3S3NJUuW4OXlRfv27fHz86NRo0Yf7PcxKlmyJAMHDmTDhg1s3LiRAQMGJJtUVqpUie3bt3Pp0iW8vb21ezobGhpSvXp1mjVrRocOHfLMv53UREZGcunSJdzc3Ni6dSvW1taYmJig0WgICwsjPj4eGxsbnJ2dadCgQa6vV+nn54eiKGn6hato0aLJzkamxMHBgU2bNjFkyBCqV69OmTJlkv08Ro8eTdeuXWnTpg12dnbpCT9Z06dP59tvv5VkUoiPkMxQ5lKKorBz505WrlxJo0aN+PLLL7G3t0+SPL18+ZKDBw+yfft2GjRowA8//JCnkoOQkBCWLl3KqVOnaNGiBQ0aNKBWrVoULFgQjUbD48eP8fLy4vjx4wQGBjJ48GC6devGkydP2LBhA/Xq1UvT7jEfq5cvX7JhwwYKFCjAgAEDtMlHbGws8+fPZ/v27ZQsWRJ7e3tKlSqlnbFSq9UEBgYSEBCAn58fDRs2ZPr06RQuXDgnLydD1Go1c+bMYdeuXdpZtpTKMwUHB+Pt7U1ISAiLFi1KtoZibrFy5UouXLiAo6Njim3ez1AOGjQoXTVHo6Ki8PT01C6SSc2TJ08YOHAgGzdupEyZMmk+x38tW7YMPz8/5s2bl+ExhBC5lySUuVBQUBBjx47F3t6eKVOmpKmOoKIo7N+/n8WLF/Pbb79Rp06dbIhUdxISEjh37hxXrlzh2rVrhIWFoa+vT+nSpalTpw6NGzdOMkNy8eJFPDw8+Oyzz/J1CZJXr16xfv16LC0t+eqrr/Dx8WH06NE4ODhQo0aND87gKorCo0ePOHv2LFOmTKFnz57ZFHnmBQQE0K9fP2xtbXFyckrzrGNERAQeHh40bdqUWbNm5crZyoEDB2Jra5tqkv8+obSwsCAqKgpDQ0PKlStHu3btKF++fKrjb9u2DU9PzzTF4ufnx7Bhw/j5559p1KhRuq4jNjYWFxcX4uLimDt3bq78rIUQmScJZS7z7NkzBgwYwNy5cxNtU5hWb9684euvv2bUqFEf/cydoijs3bsXHx8fhgwZQokSJXI6pBzz6tUrNmzYwL1797h58yZdunRJd+Ho+Ph4jh49ipOTU65Nsv7t0aNH9OrVK0O1Gd+7dOkSRkZGrF27Ntddb8uWLenatWuqcT19+pRLly5hZ2eHubk5QUFBnDhxgtevXzNs2DCqVKmSYt8tW7Yk2VQgNSEhIUyYMIECBQrw/fffU7x48VTbK4rCyZMn+fXXXxk5ciSfffZZms8lhMh7JKHMRcLDw+nRowd//vlnpp5Xio2NpU+fPkyZMoV69erpMMLcJz4+nrVr1xITE8Pw4cM/qvqD6bVx40YWLlxIv379MrUN5alTp6hZs2ay5WFyi4iICFq0aEGnTp0oVKhQpsa6ePEi5cuX5+eff9ZRdJkTGhqKn58fQ4YM0W5Vmh5RUVHaxW6TJk1Ksd2WLVs4f/58uhPpU6dOsXz5cgwNDWndujVOTk5UqFABQ0NDQkNDuXbtGpcuXeLkyZM0adKE0aNHZzjhF0LkHZJQ5iKjR4+mV69etGjRItNjvX37ll69enHw4MGPPskKCQlh5cqVlCpVir59++bLRTovX76kY8eO9O3bN9O1/d7P/P7666+5tvD00KFDMTExSbKyOaN27tzJ4sWLqVWrlk7GS41arSYgIICHDx9qX35+fto/v3nzBgArKyumTp2aob/P7du34+npyZw5c1J8ZCa9M5T/9fz5c86ePYuXlxf+/v7Ex8dToEABatSogZOTE40aNUrXVqNCiLxNVnnnEmfPntWu3h48eDCenp4EBARQsGBB6tSpw7Rp0xItIFAUhSVLlrBs2TIePXqEtbU13bt355dffqFQoUIUKlSI//3vf9o9dz9mBQsWpFevXmzatImTJ09+cBeQj9Hw4cNp1apVssmHm5sbly9fTrHvf+sVqlQq2rZti7OzM2fPns11e4F7e3vz8OHDD5ZCevjwIUePHsXf35+EhASsrKyoW7duso+CtGvXDmdnZ/755x+dxBgWFpZiwvg+ng+Jj48nMDCQsmXLZjiOlGYfFUXJ9C8epUqV4vPPP+fzzz/P1DhCiI+DJJS5xJIlS1i0aBHffvstwcHBODs74+DgwOvXr5k3bx4NGjTgyJEjtGzZEoAJEyawcOFCJkyYQOvWrfHx8WHatGlcvnyZ8+fPY2hoSKdOnVi6dClRUVHprvWY15QvX55WrVpx7NgxSpUqleqzYx+bO3fuEBkZSenSpZM93q5dOxo3bpzk/VWrVmFgYJDsDizm5uZUqFCBHTt20K9fP53HnBlz5sz54MIQLy8vNm3aRM2aNbV7nwcFBREaGppse0tLSywtLbl8+TJ169b9YAxqtZpnz54lmzA+fPiQoKCgNF+Pvr4+tra2lC9fngoVKlC+fHnKly/P48ePuXz5croTyqioKG7fvk3p0qVTnCEMCgrSSRkgIYR4T2555wLPnj1j0qRJuLm58erVK4oVK5boeEREBHZ2dlSrVo1jx47x7NkzbG1tGTVqFIsXL9a227JlC3379mXlypUMGzYMgA0bNpCQkMDgwYOz9ZpywvtSSw8ePGDo0KEULVo0p0PKFsOGDaNgwYLpKuny4MEDli5dStu2benYsWOybWJiYjh06JDOZu10ITw8nHbt2tGnT58U24SEhPDLL79Qt25devfuneaxg4KC8PPzY/PmzdpzPXr0KNmE0d/fn7i4uDSPbWVlpU0W/500VqhQgbJlyyY7WxgREUGrVq3o27dviuNu2LCBQoUKUbZsWe2inJMnTxIUFMQ333yT4g5JJ06cwNnZWfsLqhBCZJbMUOYCJ0+epFOnTgBJkkkACwsLHBwcCAgIAODChQuo1eokiUDnzp0B2LVrlzah7NKlC2PGjMkXCaVKpaJbt26sXr2abdu2aZ+z+9jdunWLL774Il19Lly4gEqlon79+im2ef/ZvX37NtMLX3TlypUrKc7EvnfhwgXi4uLS/ehDkSJFcHNzo2HDhvj5+fH69es099XT08PGxiZRsvjvpDEjn5+FhQXVq1fn6dOnKf6yUKpUKa5evcq5c+eIi4vDzMyM8uXL079//xT3/k5ISCAwMFAnz2oLIcR7klDmAl5eXowYMSLF46GhoXh7e2tnE97PjPx3WzVDQ0NUKhU3btzQvleoUCHevn2bBVHnTkZGRvTp04dVq1axd+9e+vTpk+vKwehScHBwup9xjI6O5vr161SsWPGDhcyLFy+Ot7d3rnku9fz588n+0vVvfn5+mJmZ8fLlS1avXs2LFy8wMzOjRo0adO3aNdVfMhRFSXGhSoECBZIkjO+/trGxyZIFKFOnTqV37958+eWXyS42a926Na1bt07XmP/88w/Ozs4f9feFECL7SUKZCzx8+DDV1aqjR48mMjKSH3/8EXi3bRrAuXPnEs0yeHp6oigKwcHBifrr6+ujKEq++QFSuHBhPvvsM7Zs2cI///xD8+bNczqkLHPnzp10l2Tx9vYmPj6eBg0afLBtkSJFGDVqFMWKFdNu96nRaBL9+b9fZ6RdWo+ZmZkxduzYVGMODQ0lPj6ev/76i9atW1OuXDkCAgI4fPgwgYGBjB07NsXvhWLFilG8eHGqVauW5Pa0tbV1tn8P2draMmDAAM6cOZPsc7Dp9eTJE2JjY1O9jS6EEBkhCWUuoChKiqVupk6dipubG0uWLNGu8nZ0dKRZs2b8/vvvVKpUiTZt2uDj48OIESPQ19dPMpaBgQFqtTrTqzrzEnt7e5o3b86pU6coVaqUzsrL5DZRUVHp/nu9cOEC5ubm1KhR44NtjYyMePLkCb6+vhkNUafS8ouRoijEx8fTuXNn7exdxYoV0dfXZ8+ePfj6+qb4bGHRokX5/fffc9WirjFjxnD27Flu3LiRpr+zlLx69Yq1a9diZWXF48eP07X3txBCfEj+yTByMQMDA+Li4pLcunR1ddWW/RkzZkyiYzt27GDQoEHakh1GRkaMGzeOY8eOERISkqhtfHx8vkom32vevDmBgYHs2rWL4cOHY21tndMh6ZyxsTEajSbN7Z8/f05AQADNmjVL07+JhIQEbG1tKVWqFHp6eujp6aFSqZL9c1Yc+2+7M2fOEBUVleqt+vcVDSpXrpzo/SpVqrBnzx6ePn2aYkIZExODlZVVmj/P7KBSqXBzc0s0U5neWqv37t3j8OHDhISEEBISQrNmzTh+/PhH+4uWECL75b8sIxdycHDg1q1bibZadHV1xcXFBRcXF6ZMmZKkT7FixXB3d+fVq1e8ePECW1tbTE1NWbZsGb169dK2e3+rMD9SqVT06NGDVatWsW3bNoYMGZLraipmVsWKFdP1jOz75wPTWrD8zZs3zJs3T7toLKft2LGDffv2pVpKp1SpUjx+/DjJ++8LWqQ2wxkREUHJkiUzH6iOGRgYsHnzZlasWMHq1atp2rQptra2H5ytffv2Lf/88w/29vbcunULFxcX5s6dq/2l4tixY1StWjWbrkII8THLf1uK5EJOTk5cuXJF+/XMmTNxcXHhp59++uD2d8WKFaNGjRpYWVnx559/EhkZmWg289mzZx9cFfsxMzExoU+fPrx9+5b9+/fzMVXJio2N5dKlSzx8+DBN7RMSErhy5Qo2NjZpTppevnyZqKB+TqtXrx6BgYGptnF0dATePV/6b++/TulWb2xsLJaWlrn2WWOVSsWIESPYt28farUaNzc3zp49i5+fH1FRUajVauLi4rT7e+/YsYMbN27w66+/8ueff2Jubs6cOXO0/0958eIFzZs3x9vbO4evTAjxMZAZylygVatW9O/fn+HDhzNv3jymTZtG+/bt6dSpU5IVp+8XUqxatQqAChUqEBISwuHDh1mzZg2//PJLopnOrVu30rVr1+y7mFyoWLFidO/enR07dlCqVKkPFsXOzdRqNadPn2bz5s3s3LmT0NBQLC0tCQ4O/uCK7Rs3bhAVFZXm2cmEhASioqIoUaKELkLXCVtbW+Li4oiNjU1S5eC9ypUrU7VqVY4cOYKiKNja2hIQEMCRI0eoWrUq5cuXT7bfjRs3cl0R9+SULl2axYsXExMTw9mzZ/H09OTSpUtERkZiZGSEvb09n332GU2bNk2yIl6lUuHi4oKZmRnff/89wcHBtGzZksOHD+fabTaFEHmDJJS5QIECBShVqhQ+Pj4cOHAAAA8PDzw8PJK0fT/DpigKCxcu5PHjx+jp6VGrVi327NmTaDs6jUaDh4cH48aNy54LycUcHBxo3Lgxx44do0SJEikmFbmRoih4eXmxefNmtm7dmmSGLiEhgX/++YcePXqkOs7FixcxMjJK837Vt2/f5ssvv8xw3Fll6NCheHh4pLpKfeDAgRw5cgRPT088PDywsrKiefPmtG/fPtn2Go2Gu3fvsnr16qwKW+dMTEwyVDYIYNKkSZiZmfHtt98SGhpKmzZtOHjwIJ9++qnuA82DNBoNt2/fxsvLi1u3bhEVFYWRkREVKlTAycmJ2rVr54sat0Kkh+yUk0s8ePCAiRMnsnv3bp3dclu0aBFmZmbaIuf5nUajwc3NjcDAQIYPH07BggVzOqRU+fr6smXLFjZv3pxklbWhoSEdOnSgb9++dO7cmY4dO9KiRQudFSCPj49ny5YtnDx5kgIFCuhkTF2Jj4+nefPmdOjQQWexnT9/noYNG+Ls7KyT8fKKtWvXMnToUBRFwcTEhD179qSYdOcHoaGhrF69Gnd3d2rXrk3dunWpXr06FhYWxMXF4evri5eXF+fOncPBwYHRo0fnqV9OhchKklDmIvPnz8fS0lInCeD9+/eZOHEie/bsybXPhOWEqKgoVq1ahampKV9//bW2GLWiKBw/foLtu/7m2rXbREYrgB4qlYYCFobUdapO/35d07TPc2Y8f/6cbdu2sXnz5kTP1cK725XNmzenb9++9OzZM9Gq9bt37zJo0CC++OILnfx9Hzt2jEGDBtGzZ89Mj5UVbty4wTfffMPnn3+e7hXP//Xy5UsuXrzI0aNHMz1WXrRlyxYGDBiAWq3G0NCQbdu2fXC2+2N06NAh5s2bh7OzM507d0ZfXz/V9leuXOGXX36hQYMGjB8//oPthfjYSUKZi6jVar788ksGDx6cqVmCZ8+e0aNHD7Zt28Ynn3yiwwg/DoGBgaxdu5aqVavSpUsXlv25juUrtpBg1Bi9Al0xsayMSu//73qiUccQHXYTwnZhqrrNlMkj6dNbdz9w3759y+7du9m8eTMnT55MsnCodu3a9O3blz59+qS6X/eiRYs4evQorVu3zlRSeePGDdRqNX/99VeGx8gOa9aswc3NjS5dumQ4EXz9+jVHjhxh//79uXJ1d3Z5v6tUXFwc+vr6bNiwId8UP9doNEyaNAlFUfj5558xNTVNc19FUVi/fj27du1i06ZNua7klBDZSRLKXCYmJoYBAwbQrl07hgwZku7E4MqVK4wbN466detSs2ZNvvzyyyzZEi6vu3HjBuvWreP46evEGHTCrOQ3iZLIlGjU0UQ/m41dyQA2rl/8wYUwKYmOjubgwYNs3rwZd3d37Xaa71WoUIF+/frx5ZdfJqmnmBpXV1fOnz9P27Zt0z1joigKV65cITY2li1btuSJ2qVr1qxh9erVdOjQId0/zO/cucO1a9fYtWtXvq6E8N6RI0fo3r07MTExqFQqVq1axZAhQ3I6rCylKAr/+9//qFixIqNGjcrwOBcvXmT69Ons3LkTCwsLHUYoRN4hCWUupFarmT9/PufPn+e3335LU/HhkJAQfv/9d548ecKCBQuIjIzEzc2N0qVL8+WXX3509Rcz68aNG3TsNoICdqsxtqiQ7v7RIV7oB//EEfcNqc4a/ltCQgInTpzAzc2NPXv2EB4enuh48eLF+eKLL+jXrx916tTJ8Cyjm5sb8+fPp2XLlmmOLSQkhCNHjtCmTRumTZuWp27f3bx5k2+++QYbGxtq1ar1wcUSgYGBnD9/HkdHR3777bd0zUh97E6dOkXnzp2JjIwEYPHixXz77f9r767Dqr7bB46/DyUgoqCICtiJLchwj7Oxu6Zit2JsM6dzbuqs6eyY3YXOwtbZzpwBYlOKoLQ0HDi/P7w4PxndiPfruryeZ9/4nM/XONzfT9z3+DzuVc7Zu3cvrq6uzJ07N8ttXbhwAUdHR9avX58NPRPi8yMBZT729OlTfv/9dwICAujWrRvW1tZUr15dXZvb29ube/fucfLkSXx8fBg1alSiBNReXl7s3r2bMmXKSFD5CQ8PD1q1G4pehZ1oFcpYHexPRYe/Al8Hblw5nOIGH5VKxa1bt9izZw/79+/n/fv3ic4bGhrSo0cP+vXrR/PmzbMtkHv37h0//PADbm5u6lQ5hQsXTtz/6Gg8PT1xcXFBT0+PJUuWUKtWrWz5/NwWGxvLgQMH+PPPP4mNjcXCwgIzMzP09fWJj48nMDCQ9+/f4+PjQ/Xq1Zk8efJn+6w57ebNm7Rt25aQkBAAFi5cyLRp0/K4V9nv3bt3DBgwgN9++41ffvkFZ2dn/Pz80NPTo1q1ajg4ONC/f3/19YMHD2b79u1J2qlWrRpPnz4FUK+/tLOzy7XnECK/kIDyMxAYGMiZM2e4d+8ez58/R6VSoVKpKF26NFZWVjRv3jzFUnIJQWXp0qXp16/fFx9UxsfH83WTznzQX0yhwhkfmfyvyKB/sDTZyV+OmxMdd3V1Zc+ePezZswd3d/dE5woVKkTHjh3p168f7du3z9H0I4GBgWzfvp3Lly8nG8x+9dVXDBkypMDUdVapVMyePRstLS2io6MJCAhAW1ubihUrYmtrS926ddWlGUXK/v33X1q3bk1AQAAAP//8M7/88kuB2uA3depUunXrRnR0NPv27aNx48aYmZmpZ3f27dvH3Llz+emnn4CPAeWBAwf4+++/E7Wjp6enTqYfHh5Ojx49kk35JkRBJwHlF0CCyv+3fMWfrN0TSeEyY5I9HxH8AN+nvxERfJe42GB09CwwMv+WkpW/R0Mr+UAkzGM6K35rS61aNdi3bx979uzh4cOHia7R0NCgZcuW9OvXj27dusni/RwSEBDA6tWr6devn9SpzqLHjx/TqlUrfH19AZg0aRK///57gQgqIyMj6dq1K6dPn07xeWxtbXn79i1eXl7Ax4Dy4MGDhIWFpdr25MmT+fbbb3M8I4QQ+U3+X3Uvsqxs2bL079+fXbt2sXv3buzt7b/IoDI+Pp4/N+1Hv5JTsuejPjzhxZVmFDKoilnt39HSKU6Y/zV8n80nIvg+FW0PJnufbpnp9B/UhCD/50nO2djYYG9vT+/evfNVxZmCytPTE4VCkWqtb5E+NWvW5MqVK7Rs2ZLXr1+zdOlSIiIiWL169WefXunChQt07Ngx1eC4RIkSSUb102PQoEFs27ZNAkrxxfm8vxVEullYWNC/f398fX3ZvXt3kl3FX4IzZ84RV6gZCkXy6xSD3uxDFR9Fha/2YmTWkyImzSldYxbFyw3mg68TypigZO/T0ilGnNb/j4ZVq1aNOXPm8OLFC27dusWECRMkmMwlXl5elCpVSqqYZJMqVapw9epVdfLudevWMWzYMOLi4vK4Z1lz9+5dvvrqq0TH4uPjUSqV+Pn5sXbtWs6cOZNk7WhkZCSlSpVCU1MTc3Nzxo0bR2BgYKJrLC0tefbsWY4/gxD5jYxQ5iOxsbHqcl8eHh4olUqKFClC7dq1sbKyokyZMllqPyGoTBip7NevX4r1kAuifY6n0Sw2JMXzCWmDNLUST0drahcDNFBopDyqW6L8COpXV7LsjwXUq1evQEwLfo48PT1TXE8sMqdcuXJcuXKFVq1a8fTpU7Zt20ZkZCQ7d+78bFOSOTs7M3369ETHxo4dy59//gmAjo4OK1euZNSoUerzdevWpW7duurNXJcvX2bZsmVcuHCBO3fuqNMFJWyaFOJLIwFlPuDh4cG6deu4f/8+1tbWWFtbY2VlhZaWFh8+fODhw4f89ddfvHv3Dnt7e3r16pXpKWsLCwsGDBiQaPr7SwkqH7u+oFCplNfVGVn0x+/Val4/nECZmvPQ0jEhLOAq/u6bKFFxNJpahVO8V69YfQprm6a7TrbIfiEhIQQHB1OuXLm87kqBY2ZmxuXLl2ndujUPHz5k//79REREcODAgc9yNFipVCbp94wZMxg+fDjv37/n+PHjjBs3jvDwcCZPngzA999/n+h6Ozs76tevT8+ePdm4cWOi8wlBpbxYii+JbMrJQ0qlkiVLlnD79m2mTZuGjY1Nql9AYWFh7Nq1i4MHD7J48WIaNGiQ6c9+8+YNu3btomTJkl9MUGlZryOFKhxK9Zqo0Ge43+pNdNj/T1mVqOiAWe0laf5wUHh359/bJ7KlryLjHj16xOHDh5k8eXKSFEkiewQFBdG2bVtu374NfAyqjhw58tntnO/cuTPHjh1L9ZoxY8awadMm3r59i4mJSbLXxMfHY2hoSIcOHdi/f7/6eKdOnTh+/Hi29lmI/E7WUOaRoKAgunTpQpkyZTh06BBfffVVmgGLgYEBo0ePZu/evfz++++sW7cu059vbm7OgAEDeP/+Pbt37yY6OjrTbX0u0np1ig73wO1md7R0jCnfcC+VG5+jTM35BHrt5PX90eloX0Yj8pKnpycmJiYSTOYgIyMjzp07xzfffAPAuXPnaNu2LR8+fMjjnmWMiYkJPj4+qV5jY2ODUqnEzc0t1etUKlWSTUpKpTLLfRTicyMBZR4ICQnh22+/Ze7cuQwcODDD0yImJibs2bOHly9fsmbNmkz3w8zMTB1U7tq1q8AHlRpp/Db7uM4iXhlKxa+dKGbWDYMS31Cyyg+Y1f6dQK/thPlfSfV+hUIG+/OSl5cXZcuWzetuFHiGhoacPn1anbz76tWr2NnZJdmc8l8JxRj++ecfrl27xqNHj/Jsc6CVlRV3795N9ZqLFy+ioaGh3pCUnIMHDxIREYGtra362IcPHzA0NMy2vgrxuZA1lHnAwcGBX375JUtT1gqFgiVLljBo0CDq16/P119/nal2EoLKnTt3smvXLuzt7T/LNVHpoa8Hsaq4FHd5R4Y8RLdI9SRrJfWNrD+e/+CKQYkmyd4bG/mWiuam2dthkW7h4eH4+/vTpEnyfz4ie+nr63Ps2DF69+7N8ePHuX37Ni1atODs2bOULFlSfV18fDwXLlxg9erV+Pn5oa+vr64qFRMTg7+/P5qamrRr145Ro0ZRvHjxHO13eHg4rq6uBAQE8O+//9KpUydGjhyJoaEhNjY2mJqa4u/vj6OjI/v372fKlCmYmJjg6elJv3796NOnD5UrV0ahUHD58mWWL19OzZo1GT58uPozLly4wP/+978cfQ4h8iMJKHPZoUOHKF++PPr6+nTo0CHNcl/wsWrF1KlTuXnzJlpaWrRo0YIlS5ZQsWJFVq5cSffu3XFycsr0OiYzMzMGDhyoDir79+9fIIPKOrWrc/3VU/QMayZ7Xlu3NFEfXIlThqGpZaA+HhF4EwAdPbMU244MuU/zttbZ22GRbp6engCyIScX6erqcujQIfr378+BAwd4+PAhTZs25cKFC5QpU4Zr164xadIkypQpg7W1dYqjdvHx8Tx9+pQOHTrwzTffMG/evGxd0x0VFcWTJ09wcXHB3d0dhUJBxYoVuXjxIuHh4TRq1IitW7eyfft2goODMTAwoG7duuzcuVP9XWxoaIipqSl//PEH7969Iy4ujnLlyjFhwgRmzJiRaJnF9u3b2bp1a7b1X4jPhWzKyUXx8fHY2dlx8uRJ/vnnn3SV+3r69Ck2NjbUq1eP6dOnExUVxc8//0xQUBAPHjzAxMSEffv24efnx/jx47PUv7dv37Jz506KFy9eoILKgIAAtmzZwh9//IFSvxvmdZYle12IjxPut3qhb2SDSeXxaOkUJyLwNu9e/I6OngVVm99CI4XUQeFuIzhz5FcqVKiQk48iUnDq1ClevHjBhAkT8rorX5y4uDiGDRumrnNdsWJFmjdvzqtXr2jdunWGXnSdnZ1xcXFh+/btWUr/FBMTw/Pnz3FxceHly5fExcVRvnx5atasiaWlJfr6+uzfvx8vLy+mTJmS6c/5rzt37rBly5YsrW8X4nMlAWUuSqjHPWPGjBSv+W+5r969e3Px4kVevXqlfsP39PSkSpUqfP/99yxatIjY2Fjat2/P2bNns5ymwsfHhx07dhSIoPLff/9l9erV7N27l6ioKAC09cpSw84lxcAw1O8S718sITLEhThlCDp65hiW6oBp1Slo6SQ/HaeM9qdYxHiuXvorx55FpO7PP/+kVKlSdOnSJa+78kWKj49n3LhxrFu3jsKFC9O2bVv1xp2MCgkJ4ejRo+zcuRNLS8t036dUKnn16hUuLi48e/aM2NhYzMzMqFWrFpaWlklGSFUqFd27d2fhwoXZkrs0KiqKjh07sn///hyfuhciP5KAMhcNGjSIBQsWpJqgvGPHjri6uuLm5oZSqcTQ0JCBAweyfv36RNe1adMGd3d3nj//WO5v2rRp9OnTJ1vyICYElcbGxgwYMOCzCiqjo6M5ePAgq1ev5ubNm4nO1ahRgwZWjbj1pDwGFlOz7TPD3MazcUV/mjWT9Xt5ISoqikWLFtG5c2fJA5qHVCoVjRs3pmTJkjRr1ixLbYWGhnL48GHOnz+PkZFRitfFx8fj7u6Oi4sLT58+JSoqipIlS1KrVi1q1aqV6r0Ar1+/ZsiQIRw6dIiiRYumem1qVCoVo0aNolOnTnTq1CnT7QjxOZM1lLnIz88vSTAZHx9PfHw8QUFBODo6cubMGVavXg3Aq1eviIyMpE6dOknaqlOnDufOnSMqKgpdXV1sbW25e/dutvxALV26tHpNZcI6Ij09vSy3m5Nev37Nn3/+ycaNGxPV39XQ0KBr1644ODjQvHlzAJq37IHvB1f0DNM/+pGS8IDz2NbXkWAyDyWM5sv6ybx19+5dNDU1sxxMAhQpUoSmTZsyZswY9u3bl+icSqXi9evXuLi44OrqSnh4OEZGRjRs2JBatWol2hSUFgsLC+bPn0/v3r3ZvXs3JUqUyHBflUolEyZMoEGDBhJMii+aBJS5JDIyMtmgLLVyXwEBAQAYGxsnuc/Y2BiVSkVQUBClS5emfv36LF68ONv6mxBU7tixg507dzJgwIB0BZUqlYqYmBg0NTXR0srZv14qlYpLly6xevVqjh49mqi+sImJCSNGjGDUqFFJUsns3b2GZq36EaO9OdWNNmmJ+uCCfvhy1q0+nOk2RPp5eHhw9epV7t27h5eXF3FxcRgaGmJgYIBCoVCXvhO5T6VS8f3339OuXTuioqI4e/Ys3t7evHnzhvDwcNq0aUO7du0S3ePm5sbt27d58+YNPj4+xMXFMWvWLPV0cbly5XB1deXq1as0btwYX19fXFxcePz4MSEhIeqytLVr16Z06dKZXu5jY2PDkiVL6NOnD+PHj8/QsonHjx/zww8/MHz4cHr16pWpzxeioJCAMpeEhoYmO6WSVrkvINUvyoRzxYsXJyQkJFv7XKpUqTSDSpVKxeXLl9m2bRsvX75EqVSqy44lVJFo2rQpI0aMyNDIQWpCQ0PZuXMna9aswdXVNdE5W1tbHBwc6NWrV4o7RUuXLs3JY5vo0GUYkSUXole0Xob7EBHwN3rhf3D+zF4JZHLYyZMn2bBhA8bGxtjZ2TF27FgqVKiApqYmISEhPHjwgCtXrtC+fXsaN27MhAkTkn0JEznnzp076uA+ICCAGzduYGZmRu3atZMsPUnw/Plznj17hrm5Obq6urx8+TLJNY0aNWLWrFn07NmTgIAA9PX1sbS0pFatWpQtWzbbShvWrl0bJycn5s6dy+bNmxk0aBCdO3dOtla5SqXi1q1brF+/ntjYWLZt20bp0qWzpR9CfM4koMwlmpqaxMfHJzletmxZ9Qha+/btAfjxxx8ZNGiQ+k09YaTyU4GBgSgUCnVOt7i4ODQ1k8+vmBWpBZWHDx9m4cKFlCpVijp16lCvXr0kX/AxMTG8evWK7t27Y2FhwYoVKzIdWD59+pQ1a9awfft2QkND1ccLFSpE3759cXBwwNo6fal7KlWqxJW/9zNg8ETcvCqiV2ZSqrW6EyhjAol88wtf1dXiz7WHJZjMQe/fv+e7776jYsWKbNmyJdkg0cjIiObNm9O8eXNUKhUXLlygV69eTJgwQTbo5KKVK1diZWUFfJw9WbBgAQqFgrCwsBQDytatW9O2bVsA/v7772QDymLFihEUFETRokVp27YtFStWTFKVJrvo6ury22+/ERISwo4dO+jRowdxcXFYWFhgYGBATEwMnp6eKJVK6tSpw+zZsyWrgxCfkIAylxQrVizZwPC/bGxsWL9+PW5ublhZWaGnp4ezs3OS65ydnalcubJ6w4yXl1eqm32yolSpUgwaNIgdO3awY8cOunTpwvjx44mMjKRXr17JvsUn0NHRoUaNGtSoUYM3b97Qvn17pkyZwrfffpuuz1YqlTg5ObFmzRrOnz+f6Fz58uUZM2YMQ4cOzdTap5IlS3Lm5F72HfiLxUvsCYutRHzhzugXq4+m9v/vCFVG+xMRfB/NiEOUMAxixbLvaNGiWYY/T6Sfi4sLEyZMYOnSpeleF6xQKGjVqhWNGzdm1qxZ3Lx5k/nz52fbKJZImZubGw0bNgRSn1H5VHoDw2rVqlGyZEkqV66c6f5lRNGiRRk/fjzjx49XV/cJDw9HR0cHMzMzdHSSzxAhxJdOSi/mkoQRyrQ21X9a7ktLS4tOnTrx119/JRqR8/Ly4uLFi3Tv3l197NatW/j5+XHs2DEePnyY7dPfpqamDBw4EF9fX5o3b46ZmRl2dnapBpP/ZW5uTt++fdm4cSNLly5N9Vo/Pz8WLlxIpUqV6NatW6JgsnXr1hw7doyXL18yderUTAWTn+rTuzv/3j7Fga3Dsbd7gEnseHjTA9XrnvCmB+Ya0xne5SUnHGdw48oRCSZz2IsXL/juu+/Yv39/pjaZ6erq8vvvv1OyZElmzZqVAz0UnwoPD8+xUUP4uETl2rVrOdZ+ahQKBebm5lSrVo0KFSpIMClEKmSEMhfVrFmTe/fuYW1tna5yXwC//vorDRs2pGPHjokSm5coUYJJkyap27548SJt27bF29ub+/fvAx9HRcuXL0+5cuUoX768eno8s4yNjTlz5gwdO3bE3Nw8U21oaWnRsWNHTp48iYmJCQMHDkx0/vbt26xZs4Z9+/YlqvNraGjIkCFDGDt2LFWrVs3Sc6SkXr161KtXL0faFukTGxuLg4MD27dvV/8byKzvv/+eiRMncubMGdq0aZNNPRT/9e7dO4oUKZJj7RsZGakrIQkh8i8JKHPRsGHDWLJkCdbW1ukq9wVQvXp1Ll26xLRp0+jZs2ei0osJP3D9/PyIjo5WB2cRERF4enri4eGBp6cnDx48AD5O5fw3wMzIdODs2bOpUaNGpoPJBAqFgrZt27Jy5UpatmxJ8eLFOXDgAKtXr+bOnTuJrq1Vqxbjxo3D3t5e1it+AX7//XcGDRrEs2fPmDVrFjdu3OD169cUK1YMa2trfv75Z/Vavf9SqVQ0bdqUq1ev4uDgwOrVq1mwYAGdOnWiUaNGKZb+E1kTHx+fo8sKFApFogwOQoj8SQLKXFStWjUCAwPx8vJiyJAhDBkyJF33WVlZJVk/+Knly5czcuRI9X/r6+ur1y3Cx5RFnp6e6iDz4cOHwMdRv3LlyqkDTGNj4xR/MDx79owrV67Qs2fPZM97enpy8uRJPDw8UKlUlC1blvbt21OxYsVkr9fU1MTOzo62bdvi6+uLv79/onPdu3dn3LhxfPPNN7IG7gsRHh7O1atX+fHHH+nduzcBAQFMnDgRS0tL/Pz8WLp0Kba2tpw5c4YWLVokuX/NmjVJNnbo6+szfvx4tmzZwnfffZdLT/JlMTIyIjIyMsfaDw8Pl8ozQnwGJKDMZQsXLmTChAkcPnw4WwKl+/fv4+7uzm+//ZbiNXp6elSvXp3q1asDHyuLeHl5qUcwXVxcUKlUGBgYJBrBLF68uLqPixYt4n//+1+yffby8mLVqlWULVsWe3t7AC5cuMDatWtxcHBIcSdkqVKlUCqV6mDS1NSUUaNGMXLkSMzMMp8fUnye9u7di729PQqFgjVr1iTJBtC2bVsqV67M/PnzkwSUHh4e/Pjjj+zYsSPR2mL4WH2qXbt2TJgwIUfX+n2pihcvri5tmhPevn0rCcOF+AxIQJnLqlSpQosWLVi4cCE//vhjltoKCAhg0qRJSSpJpEVXV5eqVauq1yJGR0cnCjAfP36MSqWicOHClC9fnpIlS+Ls7JxoKv5TJ0+eRE9Pj9GjR6sXrVetWpW5c+dy7NgxJk6cmGJf2rZtS3R0NPPnz6d79+6y6P0L5uTkpP67nFxqKQMDAywtLXn9+nWScyNHjsTOzo5u3bolOaelpUX9+vVxcXFJtuqUyDojIyMiIiLQ19fP9rZ9fX356quvsr1dIUT2koAyD4wfP56pU6eyePFipkyZkqmRynfv3tG/f3+WLVuW5YThhQoVokqVKlSpUgX4GGC+fv1aPUV+6tQpTE1NU7zf3d0dS0vLRMGgrq4ulSpV4tGjR4SEhKRYJ7dcuXKYmZnRp0+fLD2D+LypVCpiY2NTrRsfEhLCv//+m2R0ctOmTdy+fTtJkvtPNWzYkHv37klAmUMS6mE3atQIAFdXV2JiYoiOjgY+fl8lrOVO+K4ICwtTL1Hw8fEB4MmTJxgYGGBgYEDlypWJjo4mJiaG8uXL5/ozCSEyRgLKPKBQKFi8eDFLly7l22+/ZcWKFRmqtHD48GFWr17NqlWrsLTMej3q/ypUqBCVK1dW53179+5dilVn4GOuyOTKLCYc8/HxSTGgVCgUKJVKYmNjM5SCSBQs79+/p1SpUqle4+DgQHh4ODNnzlQf8/b2ZvLkySxevDjVPKx16tRRlzgV2a9bt24sWrQIGxsbNDU1cXR0JCgoSH3+wYMH6oAyobyij48P27ZtS9TOwYMHgY+FB8aPH8/9+/cZNmxYbj2GECILJKDMIwqFgsmTJ+Pi4sLw4cOpWbMmo0aNolKlSsler1QqOXHiBFu2bKFmzZqcPHky1SAvO7148QILC4sUz5cqVQoPDw/i4+PVa9Ti4uLUqT7Cw8NTbd/Q0JC3b99Srly57Ou0+KyEhYWlmnpm1qxZ7N69m1WrViXa5T169Gjq1q3LiBEjUm3f0NCQsLCwbOuvSExLSwsHBweOHz9O06ZNmT17dpr3VKlSheXLl6d4PiQkhNevX6vXZQsh8jcJKPNYrVq1cHJy4s6dOyxYsAAfHx+KFStGxYoV0dbWJiQkhOfPnxMXF0erVq3YunVrrtcpjo2NTXYEMsE333zDvn37OHToEHZ2dqhUKs6cOaMeoUhrSl9TU5PY2Nhs7bP4vGhra6f4d+DXX39l3rx5/Pbbb4wbN059/ODBg5w+fZpr164lSeQfExNDcHAwhQsXVrctI+A5a/DgwRw6dIjXr1+n+gKaHnFxcZw8eZLNmzen+t0jhMg/5F9qPqBQKLCxscHGxgaAoKAgPDw8UCqVFClShMqVK+fpl2qRIkWIiopKcdra1taWsLAwzp07x/Xr14GPZRGbN2/OhQsX0kyoHhMTQ+HCadfRFgVX6dKlefPmTZLjv/76K7/88gu//PILM2bMSHTOxcUFpVKJra1tkvs2btzIxo0bOXz4MF27duXFixcpprAS2UOhUDBjxgwGDhxIjx49Ml0KNi4ujmPHjuHg4EDNmjWzuZdCiJwiAWU+ZGRkhJGRUV53Q83Kyorr16+nujGnVatWNGvWDD8/PwoVKoSxsTH79+9HR0cnzUToYWFhaa6fEwWbtrY2SqUSlUqlHtGeO3cuv/zyCz/99FOyU6iDBw+mWbNmSY43b96crl27MnHiRGrVqgXAvXv3ZKdwDnN2dub8+fP89NNPbN68mUqVKlGvXr0MbToMCgpi69attGnThgEDBuRgb4UQ2U0CSpEmGxsbDh8+nOZ1Wlpa6s1FQUFB3L9/n0aNGqWaCkipVKKvry/JywW1a9fm1q1b2NrasnTpUn7++Wfatm1Lhw4duHnzZqJrbW1tKV++fIq7f83MzBIFmwnVc0TOePjwIUePHqVOnTp07twZe3t75s6dy4EDB7CxsaF8+fKp/hsPDw/n7t27XLlyhYCAAHbv3s2UKVPUmSeEEPmfBJQiTbVq1eL9+/eJNt18ysfHh4cPH2JhYYGWlhZv377l/PnzmJiY0L59+1TbdnZ2TpKIWnyZRo4cydy5c7G1teX48eMAnD59mtOnTye5VqVSpbvdx48fU7ZsWVlWkUPu37/PsWPHqF+/Pp06dUKhUKChocGcOXMYOnQoy5cvZ9euXZiammJiYoKRkRGamppERETg5+eHr68vurq6ODg4MHLkSNq1a0dERAT29vZcv35d1r4K8ZlQqDLyzSy+WAsXLuTVq1fqKcRPvX//nv379+Pj40N0dDRGRkY0aNCAli1bproTXaVSsXfvXi5cuCB1ugUAffv2Zfbs2eqqTtnB3t6en376SV2KVGSfe/fu4eTkhJWVFR06dEhxFDIuLo6nT59y+/ZtXr16RWxsLMbGxjRs2JAGDRokWmf9ww8/sGzZMgBmzpzJvHnzcuNRhBBZJAGlSJeQkBBatGhBnz59sq2azcOHDyldunSqZSPFl+XNmzeMGDECJycnNDU1s9zewYMHefToEXPmzMmG3olP3b59m1OnTmFjY0Pbtm2zbdlKVFQUNjY2ODs7o6GhweXLl2ncuHG2tC2EyDkSUIp0c3JyYtWqVWlOY6dHSEgIp06d4sqVKzKlJRLZs2cPt2/fZtmyZVkKUpydnZkyZQrHjh2Tkp7Z7ObNm5w5cwZbW1tat26d7WugXVxcsLa2Jjo6mnLlyvHw4cMUs0wIIfIHCShFhnz//ff4+PioS6xlRnh4OIcOHWL37t3ZOrUpCo4VK1bw8uVL/vjjj0y9cNy6dYuZM2eyb98+SpQokQM9/LwplUqePn3KgwcPCA4ORlNTEzMzM6ysrChTpkyqAeL169c5f/48//vf/2jZsmWObahbuXIlEydOBKB///7s3LkzRz5HCJE9JKAUGaJSqZg4cSIvXrygVatWGc6P+fbtW86dO8eWLVuoV69eznRSFAgHDhxg48aNLF26NN01uKOjo5k3bx4eHh6sXLkyX6Xfyg/u3bvHunXr8Pb2platWtSvXx9jY2Pi4+Px8vLi3r17uLu707x5c4YNG5YkndfVq1f5+++/+eabb2jevHmOZmdQqVS0a9eOM2fOAB9Hrvv27ZtjnyeEyBoJKEWGJQSVp0+fpk2bNlSqVCnNHywRERHcvHmT+Ph4Nm7cKHknRbr4+voyffp0IiIiGD58OE2bNk2y0UulUuHl5cX27du5cuUKY8aMoUePHnnU4/zJz8+P77//HmNjY8aNG0fVqlVTvDY+Pp4LFy6wYsUKWrRowcSJE9VrGS9fvkyzZs1o2rRprvTbx8eHOnXq4O/vT9GiRXn48KGUaBUin5KAUmTYgwcPOHr0KG3atOGvv/7i4sWLlCtXjjJlymBmZkahQoVQqVQEBgby5s0bXr9+jUqlYty4cXTt2lVyTooM8/X1ZevWrdy6dYvo6GgKFy6sTj0TFxeHubk53377LS1atJC/X/9x/fp1fv75ZxYtWoS1tXW671OpVGzevJkjR44wfPhwHj58SIsWLfjmm29ysLdJHTt2jC5dugDQpEkT/v7772zZsCWEyF4SUIoMCQ0NZe3atVStWpVu3boBH0snXrlyhevXr3P37l3Cw8PR0NDAwsKCxo0b07RpUypXrpzHPRcFhUqlIiYmBqVSia6urgQXqbh69SqLFy9m7969mU7Ndf/+fYYMGcLChQtp27ZtNvcwfUaNGsWGDRsAWLBgAdOnT8+TfgghUiYBpUg3lUrF/v37efPmDWPHjkVfXz+vuySESIGPjw8DBgzgyJEjWc7zevfuXRYsWMDBgwfzZAQ4PDycBg0a8Pz5c7S0tPjnn38yNNoqhMh5ScueCJGCx48f8+zZM9q3by/BpBD5WMI652XLlvHy5Us6dOhA2bJl0dPTw9jYmEaNGrFr165U72/SpAkKhYJx48ZhbW2NjY0NO3bsyMWn+H+FCxdm9+7daGlpoVQqsbe3Jzw8PE/6IoRIngSUIl3Cw8M5deoUNWrUwNLSMq+7I4RIxcWLF6lcuTK1a9cmODgYCwsL5s+fz8mTJ9mxYwfly5dnwIABKVahWbNmDS9fvkx0bNKkSezYsYPo6OjceIQkrK2t1Qnqnz9/zqRJk/KkH0KI5MmUt0iXQ4cO8erVK8aOHStlEoXI5/r27cvSpUspU6ZMitfY2try9u1bvLy8Eh338PCgdu3a7Nixg+7du+Pg4MDq1asB2LBhAwYGBvTr1y9H+5+SuLg4WrRowZUrVwA4cuSIesOOECJvyQilSNPTp09xcXGhbdu2EkwKkc+FhIQQGxubajAJUKJEiWTzyI4cORI7Ozv1prtP9evXD0dHx2zra0ZpamqyY8cOddWc4cOH4+vrm2f9EUL8PwkoRaoiIyM5ceIEVapUoXbt2nndHSFEGv7991+++uqrJMfj4+NRKpX4+fmxdu1azpw5w7Rp0xJds2nTJm7fvq0ekfwvAwMDYmJiyMuJrXLlyrFu3ToA/P39GTJkSJ72RwjxkQSUIlVnz54lNjaWjh07Sn4/IT4D9+/fp0GDBkmOjx07Fm1tbUqWLMn333/PypUrGTVqlPq8t7c3kydPZvHixamObpYvXx4PD4+c6Hq69e3bF3t7ewBOnz6dYgAshMg9ElCKFL169YoHDx5gZ2eHoaFhXndHCJEO/v7+lCxZMsnxGTNmcOfOHU6cOMHQoUMZN24cS5YsUZ8fPXo0devWZcSIEam2X7JkSfz9/bO93xm1Zs0addWcKVOm4OLiksc9EuLLlrFCzOKLER0dzfHjx6lQoUKyox1CiPxJoVAkOwVctmxZypYtC0D79u0B+PHHHxk0aBCXL1/m9OnTXLt2jZCQkET3xcTEEBwcTOHChdHW1kalUqGhkfdjEUWLFmXnzp00a9aM6Oho7O3tuX37dpLSnEKI3JH33woiXzp//jwRERF06tRJprqF+IyYmpry9u3bNK+zsbFBqVTi5uaGi4sLSqUSW1tbjIyM1L8ANm7ciJGRESdOnADg7du3yY6A5oVvvvmGH3/8EYBHjx4xY8aMPO6REF8uCShFEh4eHty9e5eWLVuqf6gIIT4PVlZW3Lt3L83rLl68iIaGBhUrVmTw4MFcvHgxyS+Arl27cvHiRRo3bgzAmzdvMDc3z9FnyIjZs2fTsGFDAP744w/Onz+fxz0S4sv0xU55x8XF8ezZM548eUJkZCSFChWiWrVqWFpaJptK40sRGxvL8ePHsbCwwMbGJq+7I4TIoHr16rF48WL1f48cORJDQ0NsbGwwNTXF398fR0dH9u/fz5QpUzAxMcHExITy5csn256ZmRnNmjUDIDAwkMKFC+erWQttbW12795NvXr1iIiIYNCgQTx69IjixYvnddeE+KJ8UZGTSqXi4sWLbNy4kdDQUKpXr07NmjXR19cnPDycPXv24Orqio6ODkOHDqVNmzZoamrmdbdz1cWLFwkJCaFv37756oeGECJ9ChcuTPHixXn16hWVKlWiUaNGbN26le3btxMcHIyBgQF169Zl586d9O/fP0Ntb9u2LcP35IYqVaqwYsUKRowYwdu3bxk5cqS67vidO3fYt28fDx8+JDIyUn2PlpYWNWrUoGPHjrRr1+6L+64XIrt9MZVynjx5wg8//ICtrS0jRoxINS2Gn58fW7Zs4fz58yxatOiL2ZTy5s0btmzZQsuWLfnf//6X190RQmTSnTt32L59e7am04mKiqJ9+/acPXs2X87iqFQqunfvzpEjRwAYN24c9+7do1ixYtSoUYMyZcok6nd8fDx+fn48e/YMT09P+vTpw/fffy+BpRCZ9EUElBs2bOD06dOsWrUKMzOzdN/n5+fHd999R7169Zg8eXKBHrFTKpX8+eef6OjoMGzYsHyxi1MIkXnDhw9n8ODB6rWPWTVt2jQaN25Mp06dsqW9nODv70+tWrUIDw+nTp06dO3aFR0dnTTvi4+P5969e3h6erJ582aqVq2aC70VomAp8FHD8uXLef78OYcOHcpQMAlgYmLCrl27iImJYfbs2TnUw/zhypUrBAYG0rlzZwkmhSgAlixZwqxZs3j//n2W2zp79ix+fn75OpiEj2vAS5QowbBhw+jdu3e6gkkADQ0NGjZsSMuWLenbty83b97M4Z4KUfAU6Mjh9OnTPH36lN9//z3To4sKhYKZM2cSGxvLnj17srmH+YOPjw/Xrl3jm2++wdTUNK+7I4TIBsWKFWPVqlX069ePd+/eZbqd8+fPs3r1atauXZuNvct+oaGhdOvWjS5dulChQoVMtWFkZESvXr1wcHDg6dOn2dxDIQq2AjvlHRwcTI8ePXByckJPTy/L7SmVSjp06MD27dspVapUNvQw5/n6+vLw4UNCQkLQ1NTEwsKCOnXqoKurq74mLi6OTZs2ER8fz8iRI2X9kBAFzJMnT3BwcGDSpEl06NAh3ffFxMTw22+/4ebmxoYNG7LlezQnDRkyBAMDAypXrpzltoKDgzl37hyXLl3Kl+tFhciPCmxAOX36dNq2bUuxYsWYOXMmzs7O+Pn5oaenR7Vq1XBwcEi0W/HatWts27aN+/fv4+LiQkxMDO7u7olSaTx69Ig1a9bw559/5sETpc/Lly9ZtmwZd+7cwcDAABMTE3R0dFCpVISFheHn54empiYDBw7E3t6eO3fucOnSJYYPH57qRiUhxOcrMjKSWbNm4enpydixY2nWrFmKszYRERHs3buX3bt3M2bMGHr16pXLvc24ixcvsmDBAjp06EBUVBRnz57F29ubN2/eEB4eTps2bWjXrp36+vj4eC5fvsyzZ8/w8fEhIiICIyMjateuTcuWLdHX1+fff/+levXq6sTpQojUFchXr6ioKO7fv8+CBQu4fPkyFhYW9O3bFzMzM8LDw9m9ezcDBgzAw8ODn376CYALFy5w/vx56tevj6GhIZcuXUrSbp06dfD19SU4OJhixYrl7kOl4cOHD4wfPx5PT0+srKywt7dP8dqYmBiuXLnC6tWrqV69ugSTQhRwenp6LFmyBG9vbzZu3MjChQsxNjamXr16FC9enLi4OLy8vHBxcSEuLo7evXtz4sSJfD8qmWDRokU0bdoUgPDwcG7cuIGZmRm1a9dOdj1kbGwsp0+fpkGDBtja2lK4cGHevHnDuXPnePz4MT/88AP169dn165dTJ48GW1t7dx+JCE+OwVyhPLQoUP4+fkxevToFK+xtbXl7du3eHl5AR/fWBM2oyxZsoQpU6YkGaEEOHDgAMHBwYwcOTLH+p9Rd+7cYfTo0fzvf/+jUqVK6b5PqVRy9uxZihQpwq5duxJNhQshCrbg4GAePnxIcHAwmpqalClThlq1aqV7I0t+4enpyeDBg+natSuAuo65QqEgLCyMn376KdkRysjISAoXLpyorQcPHqhzbVpbW/PPP//Qo0ePz2KUVoi8ViA35dy8eZMmTZqkek2JEiUSrY1J787mb775hlu3bmWpf9npn3/+YcyYMfTs2TNDwSR8TOzbvn17TE1N6dKlC9HR0TnUSyFEflOsWDGaNm1Kly5d6NixIw0aNPjsgkn4uPmyYsWK6v9WKBRpbsLU0NBIEkwClCtXDoCgoCAAqlevztGjR7Oxt0IUXAUyoHz27BnVqlVLdCw+Ph6lUomfnx9r167lzJkzTJs2LcNtly5dGl9f3+zqapa8e/cOBwcHevbsib6+fqbbqVq1KhUrVsxXo65CCJEe169fz7ba4s+fPwc+fs/Dx6Dbw8MjW9oWoqArkAElkGS38tixY9HW1qZkyZJ8//33rFy5klGjRmVL23lBpVIxYsQIWrVqlS3rnKpXr46vry+nTp3Kht4JIUTu8PLywsjIKMvtBAcH4+TkhIWFBZaWlsDH0U6lUpnltoX4EhTITTnJmTFjBsOHD+f9+/ccP36ccePGER4ezuTJkzPcVnR0NM+ePcPQ0BBDQ0P09fVzvYrO+fPniY+Px8zMLF27GhO8fv2a48eP4+HhgYaGBlWrVqVz586UKFGCli1b8vPPP9OmTRtJbi6E+CykZ4o7LeHh4WzYsAGVSsWgQYMSff8V5AppQmSnAhlQqlQq4uLiEo0kli1blrJlywLQvn17AH788UcGDRqEiYlJhtoPDg5m37596v/W0NDA0NCQIkWKJPu/Cf8/O0c2ly9fjq2tLZC+XY3wcYp89erVmJmZMXjwYGJjYzl16hSrVq1iypQpGBgYYG5uztmzZ2nbtm229VUIIXKSSqXKdOAXERHBunXrCAkJwcHBgRIlSiRpWwiRtgIZUFatWpXnz59To0aNFK+xsbFh/fr1uLm5ZSigfPfuHTVq1OCHH34gNDSUDx8+8OHDB/X/Dw0NxdfXlw8fPhAbG5vo3sKFC6cYbCb8/0KFCqXZB39/f0JCQjA0NATA2NiYBQsWqHc1phRQnjp1Ci0tLUaOHKne0W1hYcFvv/3G33//TefOnWnQoAHr16+XgFII8VmoVKkS/v7+GR4YgI/B5Nq1awkMDGTs2LFJ0qepVCpJbC5EOhXIfym2trZcvXo11YDy4sWLaGhoJNodmB5Xr17F1taWIkWKUKRIkRTzN6pUKqKjoxMFm5/+f29vb548eUJERESi+3R0dNIc7bxz506iuuTpeTOPi4vj8ePHNGzYMFF6IGNjYypXroyzszOdO3fG0NAwW2r/CiFEbvj666+5cuVKhgPKhGAyICCAMWPGJLuxJyAggCpVqmRXV8Vn7tGjR/z11188evSImJgYNDQ0UKlUxMfHU7FiRZo3b06nTp2+2LylBTKg7NixI926dWPkyJGMHDkSQ0NDbGxsMDU1xd/fH0dHR/bv38+UKVPUX0J+fn5cvnwZAGdnZ+DjiJ6JiQkmJibqpLm7d+9my5YtafZBoVCgq6uLrq4uJUuWTPE6pVKZaHTz08AzICAADw8PQkNDiY+PV99z7do19aLx9PL39yc2NjbZALhMmTI8f/6c2NhYtLW10dHRITAwEGNj4wx9hhBC5Lb27duzdetW6tevrz7m6upKTEyMOhXau3fvePDgAYD6u3P9+vV4e3vTtWtX4uPjE+3mNjAwoESJEri6ujJs2LBcexaRP50/f55ly5ZRsWJF7O3tmT59eqKBGZVKhZubGydPnqRdu3a0bNmSSZMmfZZpuLKiQAaUenp61K5dm6tXr9KoUSO2bt3K9u3bCQ4OxsDAgLp167Jz585EpRcfP36cJHnt2LFjAWjatCmXLl3CxcWF4sWLZ8uOwgRaWloYGRml2qZKpSI8PFwdbD569Eg93Z1eCSOhyaUX0tfXR6VSERERQdGiRTEwMMDPz08CSiFEvleqVCmKFy9OUFCQ+nvU0dFRnUsSPiYsTwgoZ82aBaAuanH48OEkbTZs2JC+ffvi4eEhy3++YGFhYUyZMgUtLS327t2b4s9dhUJBpUqVGD9+POPGjcPR0ZEOHTrw+++/U69evdztdB4qkAElfNxw06tXL06cOMGQIUPSvL5Zs2apLr5WKpVMnjyZzZs3Z2c300WhUGBgYICBgQFlypTJ1oD2v5+TQBaiCyE+F7NmzeKHH36ge/fuAMyePTvNe5YvX57q+dOnT+Pr68u1a9fSLJQhCp7AwED69u3LtGnTaNGiRbrvUygU9O7dm+bNmzN06FAmTpxIq1atcrCn+UeBzQ1jbGzMd999x+TJk7MlOJozZ466HnheMzU1JTQ0NEP3JIxM/nfNZsIxhUKhzmcZHh6eZKejEELkV1ZWVjRo0AAXF5dsae/9+/dcv36dt2/f0rRpU8aNG0dYWFi2tC3yv+joaOzt7VmwYEGGgslPmZiYcODAAVauXMnt27ezuYf5U4ENKAE6depE2bJlmTFjRpaCyt9//52YmBgGDRqUjb3LvK+//jrD1XpKlCiBtrY2b9++TXLOx8dHfR4+/mOSgFII8Tn57bffuH37dpYr24SEhHD69GkmTpyoXie3Zs0aatWqxfnz57OhpyK/mzNnDsOHD6dBgwZZakdPT48dO3Ywffp0IiMjs6l3+VeBDigBpk6dipmZGd9++22Gg7CAgAAGDRpETEwMCxYsyKEeZpy1tTU+Pj4ZukdTU5OaNWvy6NEjoqKi1MeDgoJ48eIFderUAT6OTsraSSHE50SpVHLs2DE6duzIjRs3cHV1zVQ73t7eHD16lL179zJ37lwePXqknu729PTEzs6OESNGEBISkp3dF/nIw4cP8fT0pEePHknObdq0Sb0ELUFcXBx//PEHbdu2xdzcHH19fWrUqMH06dMJDg6mWLFiTJ48mV9//TU3HyNPKFRfyGI5Z2dnJk2aRPPmzRk2bFiqO6+DgoLYunUrJ0+eZP78+djY2ORiT9OnTZs2NGnSRP0X+9NdjXv37qVevXrqxcCWlpbo6Ojw7t07/vjjD8zNzWnZsiVKpZJTp04RERGhTmx+7do1+vfvT5cuXfLw6YQQIn2USiUHDhzAzc2N3r17Y2FhwbRp03j06BGtWrWiSJEiabYRHR3NjRs3iI+PZ+PGjZiamqrPxcfHs379eqZOnUp4eDgAZmZmrF+/no4dO+bYc4m8MXToUGbOnEmlSpUSHff29qZmzZoULlyYkJAQ9RKIsLAwypQpQ9++fbGzs6NEiRL8+++/zJs3j9KlS3P37l309PRo3749Bw4cSBSMFjRfTEAJH78Yzp07x6ZNm4iJiaFmzZrUrFkTfX19oqKiePLkiTpl0ODBg+nQoUO+TWp7/PhxNm/eTPPmzQH49ddfE+1q/NSsWbMoXrw4kLT0YpUqVejSpQslSpQgNjaWffv2cevWrXxRr1wIIVKT8J3l5eVFnz59EgUBN2/eZM6cOQQEBFCnTh0sLCwoWrSo+nxERARv3rzh1atXhIeHM378eHr16pViXl8PDw9GjhzJuXPn1Mf69+/P8uXL1d+v4vMWGBjI8OHD+euvv5Kc69SpEwqFAmNjYw4ePKgOKOPi4ggODk7yd+DgwYP06tVLnVFmz549hIaGMmrUqFx5lrzwRQWUn4qNjeXx48c8efKEqKgoChUqRLVq1ahVq1a6qtXkNZVKRZs2bbC2tk70Np0V586dY8iQIckO9QshRH4SExPDnj17ePv2Lf369aN8+fJJromPj2fmzJmEhYXh7+/P27dv1evpjYyMsLW1pV27dulO7aJSqdi6dSs//PCDetq7ZMmSrF27Vr43C4AjR47w/v17Ro4cmej4rl27GDt2LK6urvz000+JAsqUeHl5Ua5cOebPn8+PP/5IWFgYQ4cO5cCBAzn5CHkqfw6/5QJtbe1E08KfG4VCwebNm+nYsSP9+/fPcgLVV69eoaenJ1+KQoh8Lyoqij179vDu3Tv69+9P2bJlk73Oy8sLXV1dxo0bly0ZOhQKBUOHDqVNmzaMHj0aJycn3r9/T8+ePenZsyerV6/Othd8kfvu3btHt27dEh17//493333HQsXLky2mlJK/v77bwBq1qwJfEyWX9AzBRT4TTkF2fv372nQoAG7du1SV4TIDHd3dx49esS2bduyr3NCCJEDIiMj2blzJ35+fgwcODDFYBLg6dOnqZbIzSwzMzOOHTvGrl271JsYDx48iKWlJbt375Y8vp+p58+fU7169UTHxo4dS7Vq1RgzZky62/H29mb69OlYW1snWmerp6eXpZ/V+Z0ElJ+pa9eu4eTkRM+ePVm2bBn79+/n9evXGWojPj6e69ev8/z5c5ycnChcuHAO9VYIIbIuIiKCHTt2EBQUxMCBA1MddVSpVDx9+pRq1aqluC4yKxQKBfb29ri6utKzZ0/g4xq8/v3707lzZ7y9vbP9M0XOio6OTrTk7dChQxw/fpyNGzem++9QYGAg7du3R6VSsX//fjQ0/j/MKlSoEDExMdne7/xCAsrPjEql4ty5c1y4cIEmTZrQvn17WrRowenTp3Fzc+PEiRO8e/cu1Tbi4uJwdnZm165dNG/enBMnThTonWdCiM9fWFgY27dvJzQ0lEGDBlG6dOlUr/f19SUkJIQaNWrkaL9MTU1xdHTE0dFRnT3EyckJS0tLNm/eLKOVnxE9PT118Y+wsDAcHBwYP348ZcqUITg4mODgYHVAGBwcrN71nyAoKAg7Ozu8vb05d+4cFStWTHQ+MjLys9ijkVlf7Kacz1F8fDxOTk7cv3+fNm3aYGtrm+Sa+/fvs3TpUp4/f06JEiUoXrw4urq6xMfHExgYiK+vLyqVih49ejBixAiKFSuW+w8ihBAZEBoayo4dO4iKimLgwIGYmJikec/ff//NnTt3mDx5cq5lrQgICOC7775j165d6mOtWrVi48aNyW4aEvnL3LlzsbOzw9bWFg8PDypUqJDq9V26dOHIkSPAx2CyVatWuLu7c+HCBerXr5/k+vbt23Py5Mmc6Hq+IAHlZ0KpVPLXX3/x9OlTunTpQt26dVO9Pj4+nhcvXvDgwQMCAgLQ1tbm5cuXVKlSheHDh+dSr4UQImtCQkLYvn07cXFxDBw4MN0petatW0epUqWSbLLIDU5OTowaNUpdmaxw4cIsWrSIMWPGJJoCTYtKpcLDw4N79+7h7e2NSqXC1NQUKysrKleunKG2RNpOnjyJm5sb48aNIyoqips3bya5ZuHChVy+fJlTp05RokQJatWqpQ4m3dzcOHfuHNbW1knuCwoKwsHBgT179uTGo+QJCSg/A9HR0ezfvx8vLy969epFtWrVMtXOmTNneP78OePHj8/mHgohRPYLCgpix44dAAwcOBAjI6N03RcYGMiqVavo3bt3jk95pyQ4OJgpU6awadMm9bEmTZqwefNmKleunOq9np6eLF++nKtXr1K8eHFKlCiBgYEBCoWC8PBw/P398ff3p379+vzwww9JNpKIzAkLC6NPnz44OTmleM3gwYMTpQ2KjIykadOm3L17l+XLlycphGJiYkKlSpXYsGED+vr69O/fP0efIS99sWmDPhcRERHs3r0bf39/+vfvn6Vpk5IlS3Lz5k1iYmKynGZICCFyUkBAADt27EBLS4uBAwcmSkqelidPnqClpZWk2kluKlasGBs3bqR3796MGDECT09Prly5Qp06dZg3bx4TJ05MMhUfHR3NzJkz+eeff2jYsCH9+/dPdTPImzdvGDNmDBUrVmTFihWyFj6LNDQ00NXV5dGjR+pyxGl59+4dd+7cAWDixIlJzg8aNIitW7dy8OBBjh8/nq39zW9kvDwf+/DhA1u3biU4OJjBgwdneQ1OwoJxPz+/bOidEELkDD8/P7Zt24aOjg6DBw/OUDAJH9MFVapUKV+8ONvZ2eHi4sK4ceOAjyNakyZN4n//+1+imuNeXl40a9aMkJAQevfuTYUKFdLcWWxubk7Xrl3R19enWbNmPH78OEefpaCKiori4sWLrFixgvLlyzNt2rQUN1Nt27YtUT7J8uXLo1KpUvy1bds2du3aRcuWLQv0hhyQgDLfCggIYMuWLcTGxjJ06NA0dzSmR0JA+f79+yy3JYQQOeHdu3ds27YNfX19Bg0alK5a3J8KDQ3lzZs3+Woa2MDAgFWrVnH58mX1dPetW7eoX78+8+fPx83Nje7du2NnZ0etWrUy3H7lypXp2rUrAwYMkKAyA2JiYrh69SorVqzgxo0bWFtbM2fOHOzs7Ni4cWO2fMbbt2/ZtWsXkyZNypb28jMJKPMhHx8ftmzZgra2NkOHDs22OrHa2toYGxunmVZICCHygo+PD9u3b8fQ0JBBgwZlagr32bNnKBSKTK81z0lNmjTh4cOHTJo0CQ0NDWJiYpg5cyZff/017dq1y9J3fZEiRejevTuDBg0q8BVZskqpVHLz5k1WrlzJpUuXqF27NhMmTMDOzg59fX0mTpzI2bNnE9Vtz4ygoCAGDRrEypUr0dIq+CsMZVNOPuPp6cnevXspXrw49vb26OvrZ2v7+/fvJyYmhgEDBmRru0IIkRVv3rxh9+7d6u8+PT29ZK+Lj4/nwoULXL16FWdnZ+Li4lAoFMTFxWFubk58fDx169bFwcEhl58gY27evMnQoUNxd3enS5cu1K9fn7Nnz+Lt7c2bN28IDw+nTZs2tGvXLtF9ly9f5t9//8Xf35+oqCiKFClChQoVaN26NaVLl+bFixdERUUl2gwkPoqLi+PBgwdcuXKF0NBQ6tatS9OmTZNNnxcZGUn//v3p2rVrmmtZk/P8+XPGjh3L77//nmwKoYKo4IfMeSThSy4jaR2ePXvGwYMHsbCw4Ntvv82R9RYlS5bk3r172d6uEEJklpeXF7t378bU1BR7e/tkv/tiY2NZs2YNx48fp1WrVnTu3JmZM2eqr1WpVLx+/ZobN26we/du7t69y8yZM9PcUZ1XbG1tOXv2LK1bt8bW1pbAwEBu3LiBmZkZtWvXTjZlDXzcqFmjRg3KlCmDvr4+AQEBnD9/nmXLljFp0iSqVKnC4cOHefr0ab6a9s9L8fHxODs7c/nyZYKCgqhVqxbNmjVLdURYT0+Pffv2sWjRIvr06cPixYspV65cmp8VFRXF2rVruXLlClu2bEm1NGhBIyOU2eTff/9l9+7dPHv2DJVKhba2NiqVitjYWAoVKoS1tTWDBw9OsVTYw4cPOXr0KNWqVaNHjx45Njz++PFjDh48yOTJk6XUohAiz7m7u7N3717MzMzo27dvshtpnJ2dmTRpEgMGDKBfv37pSlT+4sULJk2aRIsWLZgwYUK+zNk4Z84c3r59S40aNdSbQBQKBWFhYfz000/JjlAmx9fXl4ULF9K6dWvat2/P69evCQ0N5c8//8zpR8i0hBeAu3fv8ujRIyIjIylcuDB16tTB2toaMzOzLJfMVKlUPHnyhIsXL+Lv70+1atVo3rw5pqamGWrH2dmZxYsXEx0dTa9evbCxsaFs2bLq/n348IF///2XkydP8uDBAwYOHIi9vX2OlPzMz2SEMouuXr3KggULqFq1KkOHDqVWrVpJvriio6O5du0aU6ZMQaFQsGDBgkRvLTdv3uTMmTPUr1+fjh075ugXX8I/pPfv36dZBUAIIVITFxfHqVOnOHHiBF5eXmhoaKBQKFCpVOop6LZt29K5c+dkX5JfvXrFvn37KFu2LH369EFbWzvJNWfOnGH16tXs2LGDUqVKpbtvVapU4ejRo6xZs4ahQ4eyadOmfLeO7dSpU/Tu3RsgS8FHwlrThJ8dFhYW7Nq1i/j4+HwXSH/48IFNmzaxf/9+DA0NKVmyJCVKlEBbW5uAgAAcHR1ZvXo1YWFh2NvbM2TIkAyvpVWpVLx8+ZK///4bX19fKlWqRNeuXVOt/Z6a2rVrs3PnTnx9fXFycmLhwoW8fv0aDQ0NVCoVhQsXpkGDBvTs2ZNFixZ9cYFkgvz1r+szEh0dzfTp0wkPD2fv3r2pprUoVKgQLVu2pGXLljx//pzRo0fTs2dPhgwZwqVLl7hy5Qpff/01rVq1yvG/iMbGxmhqavLu3TsJKIUQmRIfH8+GDRtwdHSkbdu2TJs2jXLlyiX6/lKpVHh7e/PXX3/Rtm1b2rdvz4QJE9RB3fPnzzlw4AAVK1akd+/eyQZ7165dY/369Tg6OqKrq5vhfioUCsaNG8eBAwcYM2YMGzZsyDc/7ENDQ9HW1s50wBcfH098fDwBAQEcP34cAwMDvvrqK/X54sWL4+bmlq+m/A8fPsy8efOoV68e3377bbIjzQn1r5VKJffu3WPbtm3MmTOHDh06pOsz3N3duXjxIq9fv6Zs2bIMHjw4XVPV6VGqVCmGDx8u1eZSIAFlJkRGRtK3b1+GDBlCly5dMnRv1apVOX78OLNnz2bQoEFUqlSJli1b0rhx4xzqbWIaGhqYmJhI6iAhRKZ4eHgwfvx42rZty7lz51IMiBQKBebm5kyYMIHx48ezbds2OnbsyIoVK4iPj+fgwYNUrVqVnj17JhtYfPjwgZ9//pljx45lKpj8VO/evXn69Cn79u2jb9++WWoruzx8+DBdNclTMnXqVJRKJfCxGsu4ceMSVRIqUaIEd+/ezRcBZVxcHKNGjcLHx4e+ffsmOxL9X1paWtSvX5+aNWuyatUqTp48yapVq1L8+/bmzRv+/vtv3N3dKV26NPb29lSqVCnfvEB8CSSgzCCVSsWQIUMYN24crVq1ylQbmpqazJs3j5kzZxIREZFrwWQCU1NTCSiFEBn2+PFjxo8fz8aNGzNUhUahUDBkyBDs7Ozo378/1apVw87Ojm7duqW4HnLGjBn8/PPPvHz5kpkzZ+Ls7Iyfnx96enpUq1YNBweHRGXsUgscqlWrhrOzM+3bt6dVq1ZZCuSyi7+/f5YC5YkTJxIXF4e/vz+XL19mzZo1jB07Vp2zuHDhwvniez4+Pp6BAweio6ND69atM3y/jo4O7dq1486dO4wYMYJNmzYl+rP29fXl4sWLPH/+HBMTE3r37k316tUlkMwD+WtxxWdgy5YtNGzYMNPB5KfmzZvHq1evElVLyA0JI5SyH0sIkV7e3t6MHz+effv2Zbqkobm5OUePHuXx48fUrVs3xWDy/fv3vHv3jmbNmhEcHIyFhQXz58/n5MmT7Nixg/LlyzNgwADmzZunvueff/5J8mv58uUAdOvWDW1tbaZMmZJtCauzKqsBj4WFBeXLl8fa2hoHBwdUKhUnTpxQn88v3+9LliwhLi6OevXqZamdhg0bEhQUxLp164CPAbmjoyN//vkn/v7+dO/endGjR1OjRg0JJvOIjFBmgJ+fH46Ojpw4cYIHDx6k+dYcFxfHihUrOHv2LC4uLgQGBlKuXDm6dOnC9OnTKVasGCtXrmTEiBE4OTnl2j8CU1NTYmNjCQ4OTjRFIoQQyVGpVEyYMIHVq1erK25lVtGiRdm1axcODg4cP3482SnMzZs3M2zYMACaNWtGs2bNEp3v2LEj7u7ubNiwgZ9++gn4mIbnv/78808UCoW6rVatWrF48WKmTZuWrp3iOcnU1JTIyMhsaUtXVzfJzFNISAj//PMPmpqaFC9eHCMjI4oXL46xsTHGxsYULVo0x3/mvHjxgr/++os+ffoQFRWVrjybn1KpVKxatQo3NzcaN25M9+7d1UU/fHx8KFKkCJ06dUr15UTkHgkoM2DTpk189913aGpqqt+a+/bti5mZGeHh4ezevZsBAwbg4eHBTz/9RGRkJL/88gt9+/Zl+PDhlChRgn///Zd58+Zx/Phx7t69S5kyZahVqxa3bt1K9gsxJ3xaglECSiFEWvbt20fDhg2JiYmhQ4cOaU4/w8e8katWrWLr1q28fPmSQoUKYWlpyZIlS/j666/p1KkTf/75J2PGjEnyedeuXWPatGmp9qlEiRKpTumGhobi6OhI06ZN1esINTQ0sLGxwcXFhbp162bidyL71KlTJ9uqloWFheHj45Noo6Wfnx8TJkygUKFCBAYG8vz5c4KCgtQjl5qamuogM6eCzWnTpmFnZ4dCoSA8PDxdeTY/de3aNfz9/dX/raGhQatWrdi8eTMrVqzAysoq3+3c/5LJn0Q6qVQqzp8/r/6SS89bs56eHu7u7omSpzZr1oyyZcvSq1cvDh06RP/+/Rk9ejRz587NtYCySJEi6Orq8u7du3xZnkwIkX+oVCq2bt3KsWPHuHnzZpov0vBxdqZbt25cu3aNqVOn8vXXXxMeHs69e/cIDw8HYPjw4bRp04ZRo0YlGqVMqShEwq7moKAgHB0d1emEUrJv3z7Cw8OT7Mi1srLi7t27eR5QJlRBUyqV6qDI1dWVmJgYoqOjgY91zR88eACApaUlcXFxrFu3jgYNGmBiYoK2tjZ+fn5cvnwZpVJJmzZtgI9/Zh8+fKBLly6JAsO4uDhCQkIICAggMDBQ/Su5YDMhuPzvr/QGm76+vrx//56mTZsCHzOMLFiwQJ1nM62AMiAgACcnJ+zt7dmyZYv6eMIa0apVq0owmc/In0Y6vXr1CktLyzRTPHz61pww1fBfNjY2ALx+/RqAChUq5Gp9bYVCQcmSJfPFgm0hRP52584d6tati66ubrqnn1etWsWpU6e4fv16ohflT1O/aGlp0aJFCy5cuICdnZ36+OvXrylfvnySfowdO1adqFtHR4eVK1cyatSoFPu9efNmihUrRo8ePRIdr1mzZr4pS9izZ0+cnZ3Vwa2joyNBQUHq8w8ePFAHlLNmzaJo0aKUKVOGf/75h+DgYGJjYzE0NKRy5coMHTpUnafTzc2NZs2aJQn8Pg0U/yulYPPZs2fpCjaLFy+OoaGh+jP37t1LjRo11O1ndMTzwIEDVKtWjTp16iQ5V6NGDQ4ePMiIESMy1KbIWRJQptO9e/ewsrJKcjyjb80Af//9N/Dxiy1BwkaZrK5PSq+SJUvi6emZK58lhPh8nT9/ns6dO6d6zX+nn1esWEGTJk3SnHXp3Lkzu3btShRQRkZGqkfvPjVjxgyGDx/O+/fvOX78OOPGjSM8PJzJkycnufbx48fcunULBweHJDup9fT0iIqKSrVfuWXYsGE0b96c2rVro6GhwezZs9O8p0+fPmlec/fuXfbt25ehvqQVbAYHBycKNAMDA3n69CnBwcHJBpvHjh3L9ObVf/75B09PT3788cdkz5ctW5Zr165JQJnPSECZTm5ubjRp0iTJ8Yy+NXt7ezN9+nSsra3p2LGj+njVqlV59epVrgWUpqam3Lt3L9F0ixBC/NeDBw+YMGFComOpvUi/fv0aDw8POnXqxIwZM9i8eTMBAQFUq1aNqVOnMmjQIHU7NWrU4MmTJ4naLlSoULKbVcqWLauuMNa+fXsAfvzxRwYNGpQkDdDmzZsBkk1AHR0dnWyt8LxgaGjIwIED+eeff2jUqFG2tOns7Ezjxo2ztYZ0wmxbcjNu/w02AwICCAoK4v379xmucAMQHBzM0aNH6dy5c4oFQ4yMjLh48WKG2xY5SyKJdIqJiUm2xmxG3poDAwNp3749KpWK/fv3J5o+19HRITY2Nkef4VMlS5ZEpVIREBCQ4bqmQogvR2RkZJLAILUXaW9vbwC2b9+Oubk5q1evpmjRomzcuJHBgwcTExOjHlnS0tJKkt6mbNmy6Zo9sbGxYf369bi5uSUKKGNiYti5cydWVlbJpqp58uRJvlo7Pn78eI4fP463t3emSwMm8Pf35/Hjx1y+fDmbepe2lILNhPQ+GeXo6IiZmVmqAXZCeU+Rv0hAmU4GBgaEhYUlOZ7et+agoCDs7Ozw9vbm77//VpeXShAaGkrhwoVz8AkS09XVxcPDg2XLlqGnp4empibm5uZYWVlRo0YNGbUUQgDJr31L7UU6Pj4egKioKE6ePKkue2dnZ4e1tTVz5sxJNFX53/a1tLRQKpWoVKpU191dvHgRDQ2NJN+lx44dw9/fnzlz5iR73927d+natWu6nj03aGhocODAATp06EDjxo0xNzfPVDt+fn6cOnWKI0eOJDv48Tl48OABT548YcKECUlGqePi4oiIiKBQoUJoampKrsl8SKKGdKpZsyaPHj2iefPmqV6X3FtzUFAQrVq1wt3dnQsXLiS7yNjV1TXNNBlZFR8fz7lz59i8eTMqlYo6derw1VdfYWxsTFxcHJ6enuzatYv79+9Tr149Ro8eneTLWgjxZUkIED+V2ot0wkhV9erVE9VQVigUtGnThgULFiRaL57cSJONjQ2XL1+mWbNmjBw5EkNDQ2xsbDA1NVUntN6/fz9TpkxJdrpbT0+Pfv36JWlXpVJx8+ZNfv7550z+buQMY2NjTpw4wbfffouXlxe2trbprvGtUql48OABL1684OjRo5kOSLNbwihiRgI/Hx8f4uPj1QnpP5WQrH7o0KHUrFlTBj3yIfkTSScrKyt27dqV5nX/fWtOCCbd3Nw4d+4c9evXT3KPSqUiLCws2YXo2cXd3Z2JEyfSsGHDFJMTf/XVV/Tu3RuVSsWdO3eYPn06lpaWzJgx47N94xVCZI2xsTF+fn6pliv89EXaysoqxe+yhOAxIViKiIhI9rtl+PDhTJ06lWbNmtGoUSO2bt3K9u3bCQ4OxsDAgLp167Jz584kuS9fv37N2bNn6d+/f7Lr765cuUKjRo3SVUs6txkbG6vXom7evJl69epRp06dFAOn+Ph4nj59yv3792nfvj2bNm3KV0FW5cqVCQgIoESJEum+x8bGJtna42vWrKF27do0adKE0qVL8/79+0SbWkX+kH/+9uVzJiYmBAUFERYWhoGBQbremiMjI2nTpg33799n+fLlKJXKRLm3TExMqFSpElevXqVhw4Y51vdjx46xfv16Vq1ala6SaQqFAhsbG/bv38+BAwfo1KkTO3bskLWWQnyBrK2tuXv3bqoVTT59kdbS0qJLly4cPHgQDw8PdQoglUrF6dOnqVSpkjrIePDgQbL5IM3NzdHX1+f27dsMGTKEIUOGpKuvFhYWxMXFJXsuLi6ORYsW5ZuUQcnR0NBQJyO/cuUKhw8fRqVSUbJkSXWQHhkZyatXr9DU1KRnz57Mnz8/2c0yea1ly5acPn06UUCZVp7NlDb+wMcKS1WqVAE+Zl2pVasWb9++pXTp0jL9nU8oVLKyNd327t1LcHAwY8aMYevWrWzdupUnT54kemsePny4+q3Zw8MjUeWC/xo0aBDbtm3D3t6eRYsW5chUxdGjR9m3bx/bt2/P9CjjkydPGDduHPv27Ut1lEIIUfC4ubkxd+5ctm7dmuaL9OLFi4GPeXutra0xNTXll19+wdDQkE2bNnHkyBEOHDhAz549AZgwYQLDhg1LNqgMCAigd+/enDhxIknqn8z4448/KFKkSL5PNfPu3TvWr19P9+7dqV27NpGRkTx8+BA/Pz/i4+MpXrw4t27dolq1aokyheQ3kZGRNG3aFHt7e/WxX3/9NVGezU/NmjUrxWDyu+++o3HjxvTs2ROVSsWmTZsYMGAAUVFRlChRgtq1a1O7du1cq/wWGhrKjRs3uHfvHs+ePSM6Oho9PT1q1KiBlZUVjRo1ytEZx/xKAsoMiImJoV27dhw4cCDb3givXLnCnj17WL9+fba096mXL18yceJEDh8+nOUpa2dnZ2bNmsXhw4flbVCIL0yPHj1Yv349Tk5Oab5IJ3BxcWH69OlcuXKF2NhY6tWrx8yZM9VBUFhYGN9++y0nTpxI8XPPnDnD9u3b2bFjR5amc0+dOsWOHTvYs2dPvv/+OnbsmPq7O6X61MeOHePt27eMHj06l3uXMRMnTiQuLk49spgdXF1dKVGiBHPnzsXd3Z1Hjx7x5MkTYmNjsbCwoE6dOlhaWuZIQPfs2TNWrFiBl5cXrVq1wtramurVq6Onp0dERASPHz/m7t27nD9/HktLSyZMmJBskv6CSgLKDLpz5w4rVqxg586dWf5iCg0NpUuXLhw5cgRDQ8Ns6uFH8fHxdOrUiTVr1hAcHMzMmTPTrL+7cuVK9uzZw8uXLwkNDcXU1JSvv/6aWbNmUbNmzc/mDV8Ikb1u3LjB1q1b2bhxY7a1OWnSJOzs7Gjbtm2q1x08eJADBw6wYcMGihUrluHP2b17N4cPH2bXrl3ZMtKZkyIiIli2bBlNmjThm2++SfG6+/fvc/z4caZNm5Zvcmom58OHDzRv3pxvv/02W/oZERHBoUOHuHLlSqKAMSYmhmfPnuHs7MzLly9RKBRUqVKF2rVrU7Vq1SyvmVUqlfz+++/8+++/zJ49m1q1aqV5z927d5kzZw6tW7dm7Nix6d5k9Tkr+E+YzRo2bIi1tTUzZ87MUh6siIgI7O3tmTdvXrYHk/DxDbZZs2aUL1+e4OBgLCwsmD9/PidPnmTHjh2UL1+eAQMGMG/ePPU9AQEBtGvXjk2bNnH27Fl+/fVX7t+/z1dffcWzZ8+YOHEi+/fvV69/EUJ8Gb7++muKFCnCyZMns6W9a9eu8eHDhzSDSfhYnvC7776jR48enDhxIt3fu76+vvTr148XL16wd+/efB9Mwse1gUCyVdk+ZW5ujkqlUuf8zK8MDQ2ZP38+J06cSDZbQEbEx8fj5OTEsmXLkow+6ujoULt2bfr168ekSZNo3bo1YWFhHDx4kKVLl3L06FHc3d0z1YfQ0FC6deuGmZkZBw4cSFcwCR/XHh85cgSFQkGfPn3yTXWmnCQjlJm0ZMkSXrx4wdKlSzNcDcDDw4PRo0czbdq0NNMQZVbXrl3ZunVrqmtKbG1tefv2LV5eXile8+TJEywtLZk1axZz5sxhw4YNFC5cONG6GCFEwRcREUH37t2ZN28e1tbWmW7H1dWViRMncvDgwRQroaT0+YsWLeLGjRv06NEDOzs7KlasmGimKCQkhNu3b7Njxw4iIyOZOXNmspk18qO4uDhWrFhB5cqV0yx1qVKpWLx4Mba2tjRt2jSXeph569atw9HRkY4dO6Y4jZ+auLg4jh49ytChQxNVWkpLQEAAzs7OPHr0iKCgIIoUKULt2rWpU6dOujaZRkZG0rNnT2bOnMnXX3+d4X4nOH36NJs3b2bv3r35aid+diu4T5bDJk+erK5xO2nSJNq1a5fmkHZERASbNm1Sp4ZILj1Cdnj37h0GBgZpLlD+b/3d5CRswkn4R9CvXz8GDBggAaUQXxh9fX327dtH3759GT16NF26dMlwG+fOnWPx4sXs3r07Q8Fkwuf/+uuvREdHc/ToUZYvX46bmxsaGhooFAri4uIoUqQIVlZWzJ0797Nbu+bq6kpoaChfffVVmtcqFAosLCx4/fp1LvQs68aMGUPRokVZtmwZdnZ2GcoY8vbtWy5cuMCPP/5Ijx49MvS5xYsXp1mzZjRt2hRvb28ePXrEgwcPuHHjBiVLllRv5knp7+LMmTMZO3ZsloJJgLZt2/L+/XsWL17MjBkzstRWfiYBZRa0atUKGxsbVq1axR9//EGLFi34+uuvqV+/PsWKFUOlUvH27Vvu3bvHxYsXefLkCf369eP48eM5up7i7t27yZatSq3+7qfi4uJQKpW4u7szffp0SpYsqU7bYWBgQExMTIYT1gohPn/FihXj8OHD/Prrrxw+fJh58+alKzuFr68vv/76Kzo6Ohw5ciRLVcEKFSpE79696d27d6bbyI9u3bpFhQoV0h1smZubc+PGjc/mu7hfv340adKEUaNGER0djZWVFWXKlEnx+tevX3P//n0MDAxwcnKiVKlSmf5shUKBubk55ubmtGnThlevXuHs7Mzly5e5cOEC5cuXp3bt2lhaWqqXRly/fp3Q0FA6dOiQpL1NmzYxYsQIChcunKiC3rVr19i2bRv379/HxcWFmJgY3N3d1UvM+vTpg4uLS7qnzT83MuWdTZRKJdeuXePOnTvcv39f/ZfM1NQUKysrvv7662Qr5OSE3377jebNmyd5qxo9enSi+rvLly9nzJgxSe7X1dVVr5OsWrUqR44coUaNGonamTlzJhYWFjn4FEKI/OzevXssW7aM2NhYunbtipWVFZUrV0ZDQwOVSoWbmxv37t3j2LFjxMbGMnHixCyP9BRUb968YfPmzfTp0yfddcbd3d3ZsWMHY8aMSbZQRX7m6urK0qVLefLkCVpaWhgbG6OlpUVcXBwBAQGEhYWhp6fHunXrcvTnZnR0NE+ePMHZ2Rl3d3c0NDSoWrUqderUYerUqWzbtg1jY+NE93h7e1OzZk0KFy5MSEhIooDy119/ZevWrdSvX5/g4GAuXbqkDigBvLy8mDFjRrqKpHyOZIQym2hpadGsWTOaNWuW110hKCgo2bRGqdXf/dSNGzeIiYnh1atXLFu2jObNm3PhwgV1ZYLixYsTFBQkAaUQX7CE6mG+vr6cO3eOtWvX4ubmBnxc41e+fHmsrKxYsGCBfFek4datWxgZGVG1atV032NmZoZCoeD169efXUBpaWnJ5s2bAQgODub58+dER0ejq6urDqiXLVuW4/0oVKgQ9erVo169eoSGhuLi4sKjR49Ys2YNRYoUSRJMwscBlSZNmmBsbMzBgwcTnZs1axazZ88GPu6zuHTpUqLzZcuWJTY2Fl9f3yyNuOZXElAWQAkjBP+VWv3dTxOWN2jQAPi4aadz585UrlyZGTNmcPToUeDjD4vMLKwWQhQ8pUqVYsCAAQwYMCCvu/JZ+vDhA66urtjZ2WVo6lpHRwdTU1PevHmT5q7w/KxYsWLY2NgkOW5ubo6rq2uuzewVKVKERo0a0ahRI3766Sd69eqV5Jpdu3Zx+fJlXF1d+emnn5KcT89Stm+//ZZjx44xcuTIbOl3fiJpgwqg0qVLp2uxto2NDUqlUj2qkJwiRYpQvXp1nj9/rj72+vXrAvl2JYQQue3OnTtoaWllaje6ubn5Z7MxJ6MsLS15+fJlnqSpe/bsWZIg/f3793z33XcsXLgwS1XtbGxsuHv3bla7mC9JQFkAWVlZqfOZpebT+rsp8ff3x9nZOdGO9MDAwHxZO1YIIT4nsbGx3Lt3j/r162cq8beFhQUBAQFERETkQO/yVo0aNYiLi0s0mJFbIiMjk2wcGzt2LNWqVUt230FGmJmZ8fbt2yy1kV/JlHcB1KBBA5YuXar+77Tq75qYmBASEoKdnR39+vWjSpUq6Onp8fz5c1asWEF0dLR6XYi3t3ey60qEEEJkjLOzM5GRkclO+aZHwtrUN2/eZGj95eegWLFilClThidPnlC7du1c/ez/Tl0fOnSI48ePc//+/SzvqFcoFAW2ao4ElAWQgYEBxsbGvHr1ikqVKtGoUSO2bt3K9u3bE9Xf3blzp7r0oq6uLnXr1mXDhg28fv2aqKgoSpUqRbNmzTh06BCWlpYAbNy4kaFDh+bl4wkhxGdPpVJx69YtqlWrlumX9GLFilG4cGFev35d4AJK+DjtfenSJWJiYtDR0cmTPoSFheHg4MD48eMpU6YMwcHBwMdyj/BxU5G2tnaGUmEV1OQ6ElAWUA4ODixZsoR169YxZMgQdR7JlBQqVCjNOr0hISFcv35dPVophBAic9zd3Xn//n26yk+mJCHB+Zs3b7KxZ/mHpaUl58+f58WLF+osI7khIWezhoYG/v7+vHv3jqVLlyaa+UtgZGREly5dOHLkSLrajoiIyNf117OiYI67CnVptP+mLciKKVOmMHv27M8iia4QQuRnt27domTJklmu6GNhYYG3t3eWa2XnR0ZGRpQuXZonT57k6udWrlyZV69eAR+zGFy8eDHJrzZt2qCrq8vFixeZN29eutt++PAh9erVy6Ge5y0ZoSzAFi9eTJcuXahatWqqFQnSY+/evRgaGtK4ceNs6p0QQnyZAgMDef78OZ06dcryC7qpqSmenp788ccfBAcHo1QqKVKkCHXq1EmzGs3noEaNGly9epXY2Fi0tbVz5TObNWvGqVOnqFKlCrq6usnml962bRuampqJzvn5+XH58mXg4/pYgFOnTmFiYoKJiQlNmzbl1KlT2NnZ5cZj5DoJKAuwIkWK8Oeff9K/f3+2b9+e6eTC+/fv58SJE2zbti17OyiEEF+gW7duoaenl6XNJq9evWLdunU4OztTq1YtKlWqRKVKldDS0iIkJISHDx9y+PBh3r17h729PT179syzdYhZYWlpyd9//83Lly8TVWzLSR06dKBdu3aMHz8+QwH/48ePk+SvHDt2LABNmzbl3Llz/PPPP/z666/Z2t/8QkovfgHc3NwYPXo0Q4YMoU+fPun+BxIaGsq0adMoVKgQv//+O1pa8v4hhBBZER0dzR9//IGNjQ0tW7bM8P2xsbEsWrSIBw8eMG3aNBo2bJjq9WFhYezatYuDBw/y+++/ZyrfZV5bv349JiYm9OjRI9c+c8GCBVSuXDnZBOeZtW7dOrS0tBgxYkS2tZmfyBrKL0DFihU5ceIEL1++pHPnzhw5cgSlUpni9X5+fixatIiuXbvSq1cvli1bJsGkEEJkg/v376NUKtMMBJMTEBBA586dqVChAo6Ojulqw8DAgNGjR7Nnzx4WLVrEn3/+mZlu56kaNWrw/PnzVH9uZbdJkyaxYcMG/Pz8sqU9Dw8PnJycGDZsWLa0lx/JCOUXJigoiO3bt3P27Fm0tbWxtLTEyMiIuLg4PD098fLyokiRIgwePJg2bdoU2HxZQgiR2+Lj41m9ejVmZmYZHm0LDg6mV69eLFmyhLp162bq81UqFd9//322JOjOTX5+fqxdu5Y+ffqoa33nhvv37zN79mwOHDiArq5uptsJDQ2lZ8+erFmzJlGRkIJGAsovWExMDE+fPiUkJARNTU3Mzc2xsLCQXdxCCJEDnj17xr59+xg2bFiGyvepVCr69evH999/n+kk6J+2NWDAAMaNG4etrW2W2spNa9eupXTp0nTr1i1XP/fMmTOsXbuWnTt3YmhomOH7/fz8GDBgALNnz6ZRo0Y50MP8Q4afvmA6OjrUqVOHb775hq+//pqyZctKMCmEEDnk5s2bmJubZ7gW9IEDB6hevTo6Ojp06NCBsmXLoqenh7GxMY0aNWLXrl2JrlepVKxcuZLq1atTqFAhSpcuzZgxYwgKCkKhULBy5UpmzpxJZGRkdj5ejqpRowbPnj3L1WlvgDZt2jB9+nS6du3K2bNnM3TvoUOH6NOnD4sWLSrwwSRIQCmEEELkuHfv3uHh4cFXX32Vofvi4uLYsGEDP/74I8HBwVhYWDB//nxOnjzJjh07KF++PAMGDEiUC3Hy5Ml8//33dOnSBScnJ6ZPn86ePXuws7MjNjYWY2NjhgwZ8lll7qhZsybR0dG4u7vn+mc3atSI48ePc/bsWTp16sTBgweJiopK9tqIiAh27dpF+/btefjwISdPnsz0EoXPjUx5CyGEEDns2LFjvHz5kokTJ6KpqZnu+06cOIGrqytTpkxJ8RpbW1vevn2Ll5cX3t7elCtXjrFjx7Jy5Ur1NXv37qVfv35s2LCBESNGEBMTQ8eOHTlz5sxnMTOlUqlYs2YNFhYWdOnSJc/6ERISwo4dO7h48SJRUVEYGRlRqFAhIiMj+fDhA3p6erRu3Rp7e/sMlWMsCGTrrhBCCJGDwsPDefToEU2bNs1QMAkfA8E//vgj1WtKlCjB+/fvgY/T6nFxcbRv3z7RNR07dgQ+TsOOGDECHR0dateujYuLS5byYeYWhUKBpaUld+7coWPHjhn+fcwuRYsWZfz48YwfPx6ADx8+EB0dja6uLkWKFMmTPuUXElAKIYQQmRQREcHRo0f5559/cHNzU2fGiIuLw8LCAhsbG0xNTVEoFFhZWWW4/cDAQEqWLJnoWEKt6aCgIBwdHTlz5gyrV68GPm62BJLUi9bW1kahUPDo0SP1sUaNGnH37t3PIqCEj0nOr169iru7e77ZLZ2ZjToFlQSUQgghRAYFBwfz22+/4eLiQs+ePZkwYQIVK1ZUB5QqlYrXr19z7do1li5dSpEiRQgODkZfXz/dnxEeHo6BgUGS42PHjlXnk9TR0WHlypWMGjUK+Bh0AVy/fp3mzZur77lx4wYqlYqAgAD1sfr167NixQqGDBmS8d+APGBqaoqRkRGurq75JqAU/0825QghhBAZcOrUKXr06EGnTp04deoUw4YNo3Llyony9ioUCsqWLUu/fv34+++/mTt3LkOHDmXHjh2kd+tCSEgIxYoVS3J8xowZ3LlzhxMnTjB06FDGjRvHkiVLAKhbty5NmjTh999/x9HRkeDgYG7cuMHo0aPR1NRM1EdjY2OCg4Oz9HuRmxKmvZ8+fUp8fHxed0f8hwSUQgghRDqtXLmSI0eO4OTkRJMmTdJ9X506dXBycuLly5dMmzYtXUGlpqZmsoFT2bJlsba2pn379qxbt46RI0fy448/qqu6ODo68r///Y/evXtjZGRE8+bN6d69O/Xq1cPMzEzdTlxcXJ6tRcwsS0tLIiMj8fDwyOuuiP+QgFIIIYRIh61bt+Lu7s769evR09PL8P1aWlrMmTOHUqVKMWfOnDSvNzY2TlfpPxsbG5RKJW5ubgCULFmSkydP8u7dOx4+fMj79++ZM2cOz58/TxQEe3l5JQowPwelS5emWLFiuLq65nVXxH9IQCmEEEKk4cWLFxw5coSlS5dmOc3ODz/8wOvXr7l+/Xqq12lra6NUKtMczbx48SIaGhpUrFgx0fGSJUtSp04dihYtyvr16wkPD2fcuHHq8/fu3cPa2jrzD5IHFAoFNWrUSDTtrVKpiImJSfdSApEzZFOOEEIIkQqVSsUPP/zA6tWrefToETNnzsTZ2Rk/Pz/09PSoVq0aDg4O9O/fP9E9mzZtYv369bx48QJtbW1q1arF1KlT6dChA0uWLKFHjx6cOnUKHR2dFD+7evXqPHz4kHr16jFy5EgMDQ3VO8f9/f1xdHRk//79TJkyBRMTEwA2btwIQKVKlQgODubUqVNs3ryZ+fPn06BBA3Xbly5dYunSpTn0u5ZzKlSowJ9//kmPHj2IiYlBU1MTbW1tYmJiiIuLw9TUlC5dutCxY0e0tCTMyS3yOy2EEEKk4ubNm9SqVYty5crh7u6OhYUFffv2xczMjPDwcHbv3s2AAQPw8PDgp59+AmD27NnMnTuX0aNHs3DhQqKioli1ahUdO3bk0KFDdO/enZ49e3L48GG+/fbbFD97+PDhrFixgvXr19OoUSO2bt3K9u3bCQ4OxsDAgLp167Jz584kwezy5cvx9PREQ0OD+vXrc/jw4UQJwX19fYmNjaVUqVI59xuXzSIjI5k3bx7//vsvvXv3plu3bsluWvLx8WH//v20bduW7t27M3r06ESbkUTOkEo5QgghRCoGDx7ML7/8Qvny5VO85tNqNQDm5uZUqFCBq1evqq+JioqiVKlSNG3alKNHjxIaGkrfvn1xcnJK9fO7devGypUrsbCwyJbnAZg2bRqtW7emZcuW2dZmTrpz5w7Tp0/n+++/VydpT0t8fDzr16/n5MmTrFmzhnLlyuVwL79sErILIYQQKVCpVPj5+aUaTMLHajWfTq9qa2tTtGjRRNfo6uqqfwEUKVKEwoULExoammrbCxcuZOLEidm2RvDu3bv4+Ph8NsHkuXPnmDdvHgcPHkx3MAmgoaHB2LFjWbFiBUOHDuXZs2c52EshAaUQQgiRAk9Pz2RHtuLj41Eqlfj5+bF27VrOnDnDtGnT1OcnTpzI6dOn2bx5M0FBQfj4+PDDDz8QEhLChAkT1NfVr1+f+/fvp9qHatWqqXNLZpWfnx9Tp05l2bJlWW4rN9y7d49Vq1axf/9+jIyMMtVGpUqV2L17N2PGjOHdu3fZ3EORQAJKIYQQIgXPnj2jRo0aSY6PHTsWbW1tSpYsyffff5+oWg3Ad999x5o1a3BwcMDY2JgyZcqwfft2jh8/zv/+9z/1dQmJutMyceJE3r17x9KlSzM9Uunj40Pfvn1ZuXIlxYsXz1QbuSk6OpqpU6eyZcsW9ahuZpUqVYo//viD7777TnaD5xAJKIUQQogUREVFJVsuMbVqNfAxZ+XEiRMZN24c58+f5+TJk7Ru3ZouXbpw5swZ9XV6enpERUWl2Q+FQsGSJUtQKpX06dMHX1/fDD2Ho6MjAwYMYPXq1dSqVStD9+aVRYsW4eDgQIkSJZKc27RpEwqFItnSlP/++y+tWrXCwMCAYsWK0b17d9zc3KhXrx5Vq1bl2LFjudH9L45syhFCCCFScObMGV68eJEof2NyxowZw6ZNm3j79i1aWlqYmZkxdOhQVq9enei6Zs2a4enpibu7OwDHjh3D19eXkSNHprtPzs7OTJs2jVq1ajF69Ogk+ScTKJVKjh07xrZt26hXrx4//fRTqimK8pPo6Gg6duzI2bNnk+T99Pb2pmbNmhQuXJiQkBDCwsLU554+fYqNjQ316tVj+vTpREVF8fPPPxMUFMSDBw8oVKgQ/fr1S3MjlMg4SRskhBBCpKBatWocPXo0zetsbGxYv349bm5uqFQqIiMjadiwYZLrrK2tuXz5MmFhYRgYGODq6srXX3+doT7Vrl2bEydOcPv2bX777Td8fHwwMjKiQoUKaGtrExwczMuXL4mPj6d169Zs37490+sP88rhw4fp0aNHsknkR48eTZMmTTA2NubgwYOJzv38888UKlQIJycnDA0NAbCysqJKlSosWbKERYsWUbZsWR4/fkzNmjVz5Vm+FBJQCiGEECkoV64cnp6eaV73abWayMhI4GP+ykGDBqmvUalU3Lx5EyMjIwoXLgzA/fv3cXBwyHC/FAoFX331FV999RUAQUFBeHp6olQqMTQ0pFKlSp9dne5PXbhwgdmzZyc5vmvXLi5fvoyrq6s652cCpVKJk5MTAwcOVAeT8PHPsHnz5hw+fJhFixbRsWNHLl68KAFlNpOAUgghhEiBQqHAxMQEDw8Pypcvn+5qNd27d2fDhg0UKlSI9u3bEx0dzfbt27l+/Tpz585FoVAQGhpKREQERYoUyXI/jYyMPrtRyNR4e3snqTP+/v17vvvuOxYuXIi5uXmSe169ekVkZCR16tRJcq5OnTqcO3eOqKgorKyscHR0zLG+f6kkoBRCCCFSMWrUKDZs2MD8+fPTXa1m9+7drF69mp07d7Jlyxa0tbWpWrUqu3btol+/fuprPr1H/D8NDY0k091jx46lWrVqjBkzJtl7AgICADA2Nk5yztjYGJVKRVBQEKVLl8bf3z/7O/2Fk4BSCCGESIWtrS3z58/H09OTIUOGMGTIkDTv0dXVZfLkyUyePDnZ88HBwTg6OnLq1Kns7m6B8N9g8tChQxw/fpz79+8nu64ytXuTO5dWGyLjJG2QEEIIkQqFQsEff/zBhAkTiI+Pz5Y2p0yZwpw5cz6bXde57dMENGFhYTg4ODB+/HjKlClDcHAwwcHBxMTEAB+D8/DwcHVuzYSRyk8FBgaiUCgoVqxYtv0ZisQkoBRCCCHSUKVKFbp06cKkSZOynBh72bJlmJmZJUpwLhJTKBTExcUB4O/vr07qnrBW1MjIiL179xIeHo6RkRH29vZUqlQJPT09nJ2dk7Tn7OxM5cqV0dXV5eXLl1SqVCm3H6nAk4BSCCGESIehQ4dSoUIFRo8erd7JnRFKpZLZs2fz9u3bZHcwi/9XvXp1XF1dgY9Vbi5evJjkV5s2bdDV1eXixYvMmzcPLS0tOnXqxF9//ZWoPrqXlxcXL16ke/fuANy+fRsrK6s8ea6CTBKbCyGEEBlw8uRJli5dyi+//MI333yTrnucnZ2ZMmUK/fr1Y8CAAbKGLw03btzg+PHjLFiwIMVrBg8ezMGDB5MkNm/YsCENGjRIlNg8MDCQBw8eYGJiot6Bn1wFHpF5MkIphBBCZED79u05dOgQR48epV27dmzevJmXL18mmgpXqVR4eXmxb98+unTpwurVq9myZQsDBw6UYDIdGjVqxN27d9NVlvJT1atX59KlS2hra9OzZ08GDx5M5cqVuXLlCiYmJrx69YqiRYtKMJkDZIRSCCGEyKSIiAiOHj3KP//8g5ubGxoaH8dp4uLisLCwwMbGhh49elC0aNE87unn58iRI9y5c4fffvstW9pTqVT06tWLuXPnUqNGjWxpU/w/CSiFEEIIkS8NGDCA8ePHY2Njk+W2tm3bhq+vL9OnT8+Gnon/koBSCCGEEPlSYGAgvXr1Ys2aNVSvXj3T7Zw9e5YNGzawb98+tLQkBXdOkIBSCCGEEPmWr68v/fv3Z+rUqbRu3TpD96pUKjZt2sSZM2fYuXMnenp6OdRLIZtyhBBCCJFvlSpViqNHj3LixAlGjRqFt7d3uu5zdXWlW7dufPjwgf3790swmcNkhFIIIYQQn4WbN2+ycuVK4uPjad++PVZWVlSrVg0tLS2io6Nxdnbm7t27nDhxAnNzc3744QeqVKmS193+IkhAKYQQQojPio+PD5cuXeLevXu8fPmSuLg4dHR0sLS0xMrKiubNm8vO+lwmAaUQQgghhMgSWUMphBBCCCGyRAJKIYQQQgiRJRJQCiGEEEKILJGAUgghhBBCZIkElEIIIYQQIkskoBRCCCGEEFkiAaUQQgghhMgSCSiFEEIIIUSWSEAphBBCCCGyRAJKIYQQQgiRJRJQCiGEEEKILJGAUgghhBBCZIkElEIIIYQQIkskoBRCCCGEEFkiAaUQQgghhMgSCSiFEEIIIUSWSEAphBBCCCGyRAJKIYQQQgiRJRJQCiGEEEKILJGAUgghhBBCZIkElEIIIYQQIkskoBRCCCGEEFkiAaUQQgghhMgSCSiFEEIIIUSWSEAphBBCCCGyRAJKIYQQQgiRJRJQCiGEEEKILJGAUgghhBBCZIkElEIIIYQQIkskoBRCCCGEEFkiAaUQQgghhMgSCSiFEEIIIUSWSEAphBBCCCGyRAJKIYQQQgiRJRJQCiGEEEKILJGAUgghhBBCZIkElEIIIYQQIkskoBRCCCGEEFnyf1WwMKlUdLKiAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/3dmoljs_load.v0": "
\n

3Dmol.js failed to load for some reason. Please check your browser console for error messages.

\n
\n", "text/html": [ "
\n", "

3Dmol.js failed to load for some reason. Please check your browser console for error messages.

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "111600" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from pathlib import Path\n", "\n", "from tqchem.molgraph import molecularSystem\n", "from tqchem.internal_coordinates import MolecularGrid\n", "from tqchem.chem import xyz_contents\n", "from tqchem.viewer import plot_3D_molecule\n", "from itertools import product\n", "\n", "\n", "molecule = molecularSystem(\"NCCCOC(CCNCC)CCCC\")\n", "n_gridpoints = 12\n", "shift_by_reference = True\n", "molecule_grid = MolecularGrid(molecule, shift_by_reference)\n", "print(molecule_grid.variable_atom_indices())\n", "\n", "molecule.draw()\n", "view = plot_3D_molecule(molecule.ase)\n", "view.show()\n", "\n", "bond_index1 = 3\n", "bond_index2 = -1\n", "\n", "xyzs = []\n", "\n", "grid = molecule_grid.shape()\n", "superindex = [0] * molecule_grid.n_variables()\n", "for i, j in product(grid[bond_index1], grid[bond_index2]):\n", " superindex[bond_index1] = i\n", " superindex[bond_index2] = j\n", " molecule = molecule_grid(superindex)\n", " xyzs.append(xyz_contents(molecule.ase))\n", "\n", "output_name = \"aminophenol_derivative\"\n", "output_dir = Path(\".output_movie\")\n", "output_dir.mkdir(exist_ok=True)\n", "trajectory_file = output_dir / f\"{output_name}.xyz\"\n", "trajectory_file.write_text(\"\".join(xyzs))" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2024-11-13 16:19:12.547 vmd_MACOSXARM64[47881:347395] +[IMKClient subclass]: chose IMKClient_Modern\n", "2024-11-13 16:19:12.547 vmd_MACOSXARM64[47881:347395] +[IMKInputSession subclass]: chose IMKInputSession_Modern\n" ] } ], "source": [ "# 1.) make sure that vmd is installed and the path is set (e.g. softlink to /usr/local/bin/vmd)\n", "# 2.) set VMDDIR to the path of the VMD installation. Without naming, the default is /usr/local/bin/scripts/vmd/vmdinit.tcl\n", "\n", "# Set the VMDDIR environment variable to the path of the VMD installation\n", "from pathlib import Path\n", "import os\n", "\n", "vmd_path = Path(\"/Applications/VMD_1.9.4a57-arm64-Rev12.app/Contents/vmd\")\n", "os.environ[\"vmd_path\"] = str(vmd_path)\n", "\n", "Path(\"vmd_movie\").mkdir(exist_ok=True)\n", "\n", "!$vmd_path/vmd_MACOSXARM64 -e movie.tcl > movie.out\n", "!magick -delay 2 vmd_movie/snapshot.*.rgb movie.gif" ] } ], "metadata": { "kernelspec": { "display_name": "tqchem", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.12.7" } }, "nbformat": 4, "nbformat_minor": 4 }