{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Defining system\n", "\n", "The default bonding information stored in the code is based on corner-sharing lead halide perovskites. For other perovskite systems, you may need to define the specific bonding information for that system. You can do that by defining a `dict` and pass it to the `system_overwrite` option. The details of how to construct such info is explained in the code piece below. Here we showed an example for a 1D BaTiS3 perovskite structure. \n", "\n", "The code provides a function `system_test` to calculate appropriate values for you, which can be called after the Trajectory object is initialized. \n", "\n", "The bonding information should look like:\n", "\n", " \n", "\n", "where the values of each pair `[[a,b], [c,d]]` being:\n", "```\n", "a: lower bound of a range of distance that covers both and only NN1 and NN2\n", "b: upper bound of a range of distance that covers both and only NN1 and NN2\n", "c: NN1 distance of the connected pair\n", "d: NN2 distance of the connected pair\n", "```" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "------------------------------------------------------------\n", "Loading Trajectory files...\n", "Current sample: test_1D\n", "Time Span: 800.0 ps\n", "Frame count: 40\n", "Reading from frame no.8\n", "Reading every 1 frame(s)\n", "Number of atoms: 640\n", "Temperature: 200K\n", " \n", "Computing octahedral tilting and distortion...\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "Progress: 100%|██████████| 32/32 [00:06<00:00, 5.11it/s]\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEPCAYAAABFpK+YAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAABgBklEQVR4nO3deXhU1fnA8e+ZyTaZ7GQjLEkIhEDYFNBgWBURJKCERVHZ2oKoUKG1tYIIQmu1ShVrsVb9scgiVaBuCALKKlBAEghL2EIgYc06meyZOb8/JkwSEiCBSSbL+TzPPDL3nnvvuUd4751zz32PkFKiKIqiNB0ae1dAURRFqVsq8CuKojQxKvAriqI0MSrwK4qiNDEq8CuKojQxDvY4qK+vrwwJCbHHoRVFURqsgwcPpkkp/e52P3YJ/CEhIRw4cOCu9pHz40/kbNmC67334NavHw5+d90WiqIo9ZoQItkW+7FL4L9bhWeTSJ0xA1lURPa6dQC4dOmCU8uWCGdnhIsz+vvuw33wYIQQdq6toihK/dLgAr80m7k0Zw6yqAg0GjCbASg4fJiCw4et5bJWf05AegY+zzxtr6oqiqLUSw0u8GetWUP+wYMA+L/0Eu6DBmHcto3cvXswG3IwFxZQnJKKKT2dK2++iUvHDrjee6+da60oilJ/CHukbOjRo4e8kz7+4osXORszDHNeHi6dOhHy+WqEQ+VrV1FKCkkjR2HOzsbBz4/QdWvVMwBFURo8IcRBKWWPu91PvR7O2b9/f5ydnXFzc8PNzQ2v4GDuiTvEyZISmv/lz1UGfQCnli1p8fbfQAhKrl0jZcZMZHFxHddeURSlfqrXgR9gzpw5GI1Gzv/rIw6EteVgeHse+O10XNq3v+V2bn374jvtBQDyDx4kKTaWnB9/RCWlUxTbCwkJwcnJibS0tArLu3XrhhCCc+fO2adi5WzdupWIiAhcXV0ZMGAAyck3HyCTkZHBiBEj0Ov1BAcHs2rVqirLvf766wgh2LJli3XZTz/9xIABA/D09OTGYevnz5+33she/wghWLhwoU3OsbrqfeAHyPvlF668+SYAzhERNJs6lZycHMaPH4+Pjw/BwcEsX74cBwcHtm3bZt3O97nncH/4YQAKT50m5fkXSH76GQqOH7fHaShKoxYaGsrq1aut348cOUJ+fr4da1QmLS2N2NhYFixYQEZGBj169OCJJ564afkXXngBJycnrly5wsqVK3nuuec4evRohTJnzpzhyy+/pHnz5hWW6/V6fvWrX/H2229X2m/r1q0xGo3Wz5EjR9BoNIwcOdI2J1pN9T7wm3NzSX1xBpSUoPH0pOU/3kfj5MSLL77I2bNnOXHiBEeOHOG7777DZDJV2FZoNLRY9B5B77yDY8uWAOT/8gvJEyZSdIurvaIoNTdu3DiWL19u/b5s2TLGjx9foUxhYSEvvfQSrVu3JiAggKlTp1ovDpmZmcTExODn54e3tzcxMTGkpKRYt+3fvz9z5swhOjoad3d3Bg0aVOkXxs2sW7eOyMhIRo8ejYuLC/PmzSM+Pp4TJ05UKpubm8vatWtZsGABbm5u9O7dm+HDh/PZZ59VKDdt2jTeeustnJycKiy/7777GDduHG3atLltvZYvX07fvn0r/TKobfU78EvJG2+/Tfefd3P/qZPcf+wo/p07YzabWblyJfPnz8ff3x8PDw/eeOONKnchNBo8Y4YStuE7AmbPBkdHzAYDKdOmYTLm1vEJKUrjFRUVhcFg4Pjx45hMJtasWcMzzzxToczLL7/MyZMniYuL4/Tp06SmpjJ//nwAzGYzkyZNIjk5mfPnz6PT6Zg2bVqF7VetWsWSJUu4evUqRUVFvPPOO9Z1Xbp0uWmXzNGjR+natav1u16vJywsrNJdPMDJkyfRarWEh4dbl3Xt2rVC2S+++AInJyceffTRGrRQZcuXL2fChAl3tY87Ua8Df/Hlyzzr7c2+duGcWfQ+WTk5ZGVlce3aNYqKiggODraWLf/nqggnJ3zGPUPgq68Clq6fS6/8CVn6HoCiKHfv+l3/5s2biYiIoEWLFtZ1Uko+/vhj3n33XXx8fHB3d2fWrFl8/vnnADRr1oyRI0fi6uqKu7s7s2fPZvv27RX2P2nSJMLDw9HpdIwZM4a4uDjrusOHD/PUU09VWS+j0Yinp2eFZZ6enuTk5NS4rNFoZNasWbz33nvVbpeq7Ny5kytXrjBq1Ki72s+dqLfj+KWUlFy7Blotrj164Pv8c9Z1fn5+ODk5kZycTFhYGGB5aFId3k+MoeDYMbLWrCFn8xbSP/oI3+eeu/2GiqLc1rhx4+jbty9JSUmVunmuXbtGXl4e3bt3ty6TUlq7aPPy8pg5cyYbN24kMzMTgJycHEwmE1qtFoDAwEDrtq6urhiNxmrVy83NDYPBUGGZwWDA3d29xmXnzp3LuHHjCA0Nrdaxb2bZsmWMHDkSNze3u9rPnai3d/wFR48hCwoA8Br7JEJTVlWNRsNTTz3FvHnzuHbtGjk5OcyePbva+w6cPQtd6Utd1xa9T95d5g1SFMUiODiY0NBQNmzYQGxsbIV1vr6+6HQ6jh49SlZWFllZWWRnZ1uD98KFC0lMTGTfvn0YDAZ27NgBYJOReJGRkcTHx1u/5+bmcubMGSIjIyuVDQ8Pp6SkhFOnTlmXxcfHW8tu3bqV999/n8DAQAIDA7lw4QJjxozhrbfeqnZ98vPz+eKLL+zSzQP1OPAbvv0WgA/T02nxxBMVhj99++23LFq0iNatWxMeHk6nTp14+OGHEULg7Ox8230LJydaLnoPbbNmAFyaN8+SAkJRlLv26aef8uOPP6LX6yss12g0TJ48mZkzZ3L16lUAUlNT2bRpE2C5u9fpdHh5eZGRkcHrr79uszqNGDGChIQE1q5dS0FBAfPnz6dLly5ERERUKqvX64mNjeW1114jNzeX3bt389VXXzFu3DjAEvgTEhKIi4sjLi6OoKAgPvroI154wTJ83Gw2U1BQQHFxMVJKCgoKKLohvqxfvx4vLy8GDBhgs3OsiXoZ+KXJhGHDBpa1Dibpd7+vMPzJaDQSExODh4cHK1asIDMzk+TkZKKjo5FS3rav/zoHPz8C/vQyAEWnz5C+ZGktnpGiNB1hYWH06FH1y6VvvfUWbdu2JSoqCg8PDwYOHEhiYiIAM2bMID8/H19fX6Kiohg8eHCNjhsZGcnKlSurXOfn58fatWuZPXs23t7e7Nu3z/psAeCNN95gyJAh1u+LFy8mPz8ff39/xo4dy4cffmi942/WrJn1bj8wMBCtVou3t7e1y2bHjh3odDoeffRR60PqQYMGVajP9RFP9koiWS9TNuTu+x/nS38CtfzXh7j371+pTFJSEpcuXeL+++8nLS2NSZMmkZubW+lh0K1IKbnw61+T+/MehLMzbb79BqdWrWp8PoqiKHWhUadsuN7No/X0xO2BB6osk5+fz5QpU/D09KRz5864urredCjXzQghCHztNYSTE7KwkMsLFqg3exVFafTq3ageWVSE4YcfACz59G94OeK6jh07kpCQcNfHcwoJodmUKaR98AG5O3aSs+kHPAY/ctf7VRRFqa/q3R2/cdduzNnZAHgMvbuXI6qr2ZTJOJW+OXflrbcw15PXzBVFUWpDvQv8hu++A8AhIADXmzwgsjWNk5PlrV6g5NIl0j/+pE6OqyiKYg/1KvCb8/LI+fFHADwefbTC2P3a5tanN24PPghA+qefUpSSWmfHVhRFqUv1KvAbt29HlnazeNxlDow7EfCnlxGOjsjCQq7+7W91fnxFacp27drFAw88gKenJz4+PkRHR7N//34Ali5dilartb7LExoayqRJkzh58mSNjmGr1Mx79+7l4YcfxsfHBz8/P0aPHs2lS5es699++206deqEu7s7oaGhFTJ1Xr16lbFjxxIUFISnpyfR0dHs27evRudxt+pV4Dd8vxEAx1atcOlU+Y262ubUujU+EycCkPPDD+Turdv/GYrSVBkMBmJiYpg+fToZGRmkpqYyd+7cCi9k9urVC6PRSHZ2Nlu2bEGn09G9e/dqD/KwZWrmzMxMpkyZwrlz50hOTsbd3Z1JkyZZt5VSsnz5cjIzM9m4cSMffPCB9b0Bo9FIz549OXjwIBkZGUyYMIGhQ4dWO/2ETUgp6/zTvXt3eaOSHKM83qWrPNY+Ql5Z+PdK6+uKyWiUJ/v0lcfaR8jTQ4dKU2Gh3eqiKE3F/v37paen503XL1myREZHR1daPnToUDly5MhqHeOjjz6SvXr1sn43Go3SxcVFHj9+vFJZo9EoHR0dZWJionXZM888I19++eUq933w4EHp5uZ202NPnz5dTps27abr3d3d5YEDB257DsABaYMYXG/u+I3btiELCwHwGFKzN/ZsSaPX4/+HPwClb/T+6192q4uiNBXh4eFotVomTJjA999/b03SdjuxsbHs3LnT+r0uUzOXt2PHjirz/oDl5nrnzp03XR8XF0dRURFt27atcn1tqDeB3/D994BlXL1zFfkz6pJHzFD0/foCkPbRvyk4dsyu9VGUxs7Dw4Ndu3YhhGDy5Mn4+fkxfPhwrly5csvtgoKCyMjIsH6vq9TM5R0+fJj58+dXOeMWwLx586xzDdzIYDAwbtw45s6dW+l4taleBH5TTg65pZn4PB4dYrf8FdcJIWj++uto3NzAZOLirNkqiZui1LIOHTqwdOlSUlJSSEhI4OLFi8yYMeOW26SmpuLj41Ot/dsyNfN1p0+fZsiQISxatIg+ffpU2s8HH3zA8uXL+e677yolkMzPz2fYsGFERUXxyiuvVOscbKVeBH7jjz8ii4sBy9u69YFjYKA1iVvhiROkffyxnWukKE1HREQEEydOvO2D2/Xr11cZcKtiy9TMAMnJyQwcOJA5c+ZYM3eW93//93+8+eabbN26lZalU79eV1hYyOOPP06LFi346KOPqlV/W6oXgd+wobSbp20YLuX61OzNc+RI9NHRAKR9+C+Mu3fbuUaK0jidOHGChQsXWufYvXDhAqtXryYqKqpSWZPJRFJSEtOnT2fbtm3MnTu3WsewZWrm1NRUHnzwQV544QWmTp1aafuVK1cya9YsNm/eXGnu3eLiYkaNGoVOp2P58uVo6vB9JStbPCGu6af8qJ6SrCx5rFNneax9hLz6jw9u+1S7rhWlpsoT3XvIY+0j5PHOXWTOtm32rpKiNDopKSly9OjRMigoSLq6usqgoCA5ZcoUmZ2dLaW0jOrRaDRSr9dLV1dX2bp1azl+/Hh57NixCvvp2LGjXLFixU2Ps3nzZtm+fXvp4uIi+/XrJ5OSkqzr/vKXv8jBgwdbv6enp8vHHntMurq6ylatWsmVK1da182bN08CUq/XV/hcFxISIh0cHCqse/bZZ6WUUm7btk0CUqfTVVi/Y8eO27YTNhrVY/e0zJlr/sPl0it2mw3f4VyNmenrWt6hQ1yYPAWz0QiOjrR8713cH3rI3tVSFKWJaTRpmbPWrQXApXPnehn0AVzvuYfWS5ag8fCA4mJSXpzBpTlzyPnxJ5XQTVGUBseugb/w1CkK4g8D4DUy9jal7UvXuRPBy5ai9fKCkhKyvviSlOef52RULy7Nm6cuAIqiNBh2DfxZ69YDIJyd7ZKbp6ZcOnQgZM3neI0ehdbPFwBZWEjW52s49+RYim6R90NRFKW+sFvgl8XFZH/9NQDuDz+M1sPDXlWpEafgYJovWEC77dsJ+eI/1gtWYWIiSSNHkbNli51rqCiKcmt2C/zGHTswpacD9b+bpypCo0HXuTNBC98hcO5r4OiI2WgkZdp0Lv7pFUqq+cq5oihKXbNb4M9auw4Ax6AgXO+/317VuGtCCLzHjiVkxWc4NG8OQPZ//8vZoTFkf/udmsNXUZR6xy6BX5aUYNy+HQDP2Ng6nXCltui6dqXN11/hNfZJAEwZGVx86SUuvTLL+layoihKfWCXiGvKygKTCQDPxx+3RxVqhdbdneZz5xK8cgVOpUNTs//7Xy5Mm4Y5L8/OtVMURbGwywtcnb285H8Cm+PaK4rgJUvq/Ph1wVxQQOrvfo+xdCpJl65d8Js2ncKTieQnJCAcHGn2m9/g0r7+pKhQFKV+s9ULXHYJ/N27dZM//PZFnMPa4NavX50fv67IkhIuzZ1LdunzjEocHGj2m1/j+9xzaG7I3KcoinKjBh34y6dsaOyklFx7bxHppRn4hKMjzu3bU3T2rLX7x7F1axyDgii5epWS9HQcg4LwHDYMj5ihOPr727P6iqLUIyrwNzD5pTP3uLRrh3ByovjSJS7PX4Dxp59uvpFGgz4qCt099+ASGYlLRHukWWI2ZGMyGECjQevpidbDA62XFxqdro7ORlEUe1CBvxGQUpKzaRNZ69ejcXLGwd8frbc3efv3k7ev5hO9a9zccPD3x8HXF+HkdP0gyJISZGEh5sJCZHHFCWWEoxMaZ2eEiwtCqwU7T4KjKErVAl75Ey5t29ok8DvYokLKnRFC4DF4MB5VTD5TnJpK9jffkPvzHgqOHsWcm3vb/ZmNRoqMRorOnq2N6iqKYkemG2YEuxsq8NdTji1a4Dt1Kr5TpyLNZorOJVN45jQaFxe0Hh6WTKFmM6ZsAyZDNqaMTEquXbM8J0hLA7OpbGcODmicnBHOzghHx7K7eimRxcXWXwOYSuxzsoqi3JbW08tm+1KBvwEQGg3ObUJxbhNq76ooitII2KWPXwiRAyTW+YHrJ18gzd6VqCdUW5RRbVFGtUWZ9lLKyrPD15C97vgTbfGAojEQQhxQbWGh2qKMaosyqi3KCCFsMiqm4SfJURRFUWpEBX5FUZQmxl6B/992Om59pNqijGqLMqotyqi2KGOTtrDLw11FURTFflRXj6IoShOjAr+iKEoTowK/oihKE2OXcfy+vr4yJCTEHodWFEVpsA4ePJgmpfS72/3YJfCHhITQULNzHt2ZysGNybTrEUCPR0NwdNZa1xUXmci4mEtOegE5GQUUFZQQ2TsIN28XO9ZYUZTGQgiRbIv9qFw9NWBIy2fHmpOYSyS/bErm1P4r9B7TDkcXLYl7L3Pm0DVKCk0Vtrl8JpvHZtxjpxoriqJUpgJ/Dfy87jTmEgkCkJCTUcD3/zpSZVkhQEpIOZFJeqqRZi3c6rayiqIoN6ECfzVdPJXJmV+uAXDvoGCat/Vk55qTGNIKANC5O9KuZwBh9/jj6a9D66Bh2ayfKSk0Ef/jBR4c18Ge1VcURbFSgb8apFmy64vTAOg8nOg+OBgnnQMt23tzcv8VXD2caNXRB6224iCpDr2ac2RbCif3XaHX42Ho3J3sUX1FUZQK6sVwTjc3N+vH0dERR0fHCssAXn75ZSIjI/Hw8CAoKIjJkyeTkZFRJ/U7sfcy187nABD1WBucdJbrpYOTlo7RQYR09q0U9AG6DGgJAkwlZhJ2pNZJXRXFHkJCQnByciItrWL25G7duiGE4Ny5c/apWDlbt24lIiICV1dXBgwYQHLyzZ+TZmRkMGLECPR6PcHBwaxatcq6buXKlRXik6urK0IIDh48CEBhYSFTp04lICAAHx8fhg0bRmqq5d//1atXGTt2LEFBQXh6ehIdHc2+O5hm9W7Vi8BvNBqtnwkTJvD0009XWAag1WpZsWIF6enpxMfHk5KSwqRJk2q9bhmXctmz3nK379vKjYhezau9rVeAKyGdfQFI2J6KqdhcK3VUlPogNDSU1atXW78fOXKE/Px8O9aoTFpaGrGxsSxYsICMjAx69OjBE088cdPyL7zwAk5OTly5coWVK1fy3HPPcfToUYBK8Wnx4sW0adOGe++9F4BFixaxZ88eDh8+zMWLF/Hy8mL69OmAJdb17NmTgwcPkpGRwYQJExg6dKg1ztWVehH4q+ONN97gnnvuwdHRET8/P6ZNm8a2bdus64uLi5k5cyb+/v4EBgbyt7/9jbZt27J06dI7PualM9mse/sg+TnFIKD36HZoNDWbjLzrgy0ByDMUcfrglTuui6LUd+PGjWP58uXW78uWLWP8+PEVyhQWFvLSSy/RunVrAgICmDp1qvXikJmZSUxMDH5+fnh7exMTE0NKSop12/79+zNnzhyio6Nxd3dn0KBBlX5h3My6deuIjIxk9OjRuLi4MG/ePOLj4zlx4kSlsrm5uaxdu5YFCxbg5uZG7969GT58OJ999lmV+75+nqJ0StOkpCQeeeQRAgICcHFx4cknn7ReNNq0acPvfvc7mjdvjlarZcqUKRQVFZGYWLfzUjWYwH+jrVu30qVLF+v3v/71r3z//ffs3buXpKQkUlJSbvlT7nbOxl3jq/cOUZhXgkYreHhSR1qEe9d4Py3ae9OshR6AuK0XUEnxlMYqKioKg8HA8ePHMZlMrFmzhmeeeaZCmZdffpmTJ08SFxfH6dOnSU1NZf78+QCYzWYmTZpEcnIy58+fR6fTMW3atArbr1q1iiVLlnD16lWKiop45513rOu6dOlSoUumvKNHj9K1a1frd71eT1hYmDUgl3fy5Em0Wi3h4eHWZV27dq2ybHJyMjt27Khwgfv1r3/N7t27uXjxInl5eaxcuZIhQ4ZUWa+4uDiKiopo27ZtletrS4MM/GvXruXjjz9m0aJF1mXLly/nj3/8I23atEGn0/HWW2+h0dT89K4kGdj0cQIbPzqCqdiMo7OWmGldCb8v8I7qKoSgy4OtAEi7YOTUfnXXrzRe1+/6N2/eTEREBC1atLCuk1Ly8ccf8+677+Lj44O7uzuzZs3i888/B6BZs2aMHDkSV1dX3N3dmT17Ntu3b6+w/0mTJhEeHo5Op2PMmDHExcVZ1x0+fJinnnqqynoZjUY8PT0rLPP09CQnJ+euyi5fvpw+ffoQGlo2H3Z4eDitW7emRYsWeHh4cPz4cV577bVK2xoMBsaNG8fcuXMrHa+2NbhRPV988QXPPvssX3/9tbVPDSA1NZXg4GDrd51Oh59f9d9svnLOwO4vTnHpTHbZPjycGDatK36t726Ky/ZRgcRtuUDmpVx+XnuakC6+OLk0uKZXlNsaN24cffv2JSkpqVI3z7Vr18jLy6N79+7WZVJKTCbLS495eXnMnDmTjRs3kpmZCUBOTg4mkwmt1vKGfGBg2Q2Yq6trtfvG3dzcMBgMFZYZDAbc3Sv/265J2eXLlzNr1qwKy5577jkKCgpIT09Hr9fzt7/9jSFDhlR4iJufn8+wYcOIiorilVdeqdY52FKDuuNfsmQJzz77LN988w0DBgyosK5FixYVunby8/O5du1atfabnmrk6/cOWYO+s6sD9w4O5slX77vroA+g1WroM6YdALnZRfyy0SZvXStKvRMcHExoaCgbNmwgNja2wjpfX190Oh1Hjx4lKyuLrKwssrOzrcF74cKFJCYmsm/fPgwGAzt27ACwSfdoZGQk8fHx1u+5ubmcOXOGyMjISmXDw8MpKSnh1KlT1mXx8fGVyl7vzhk1alSF5fHx8UycOBEfHx+cnZ2ZPn06//vf/6zPIwoLC3n88cdp0aIFH3300V2f251oMIH//fff56WXXmLTpk1ER0dXWj9u3DjefvttkpKSKCgo4JVXXsFsvv0oGmNmId9+EE9RgQmNg6D3mHaMf+MBej0ehquH7cbdt+rgQ5tull8gh7acJ/tans32rSj1yaeffsqPP/6IXq+vsFyj0TB58mRmzpzJ1atXAcsv9U2bNgGWu3udToeXlxcZGRm8/vrrNqvTiBEjSEhIYO3atRQUFDB//ny6dOlCREREpbJ6vZ7Y2Fhee+01cnNz2b17N1999RXjxo2rUG7ZsmWMHDmy0i+Bnj17snz5crKzsykuLmbx4sUEBQXh6+tLcXExo0aNQqfTsXz58jvqjraFBhP4X3zxRQwGAwMGDKg0xh/glVde4eGHH+a+++4jJCSE5s2bExQUhLOz8033WZRfwrcfxGPMLARg4MSOdH2wVa11w0SPaovWQYO5RLL7y9O1cgxFsbewsDB69OhR5bq33nqLtm3bEhUVhYeHBwMHDrSOaJkxYwb5+fn4+voSFRXF4MGDa3TcyMhIVq5cWeU6Pz8/1q5dy+zZs/H29mbfvn3WZwtgGTVY/gHs4sWLyc/Px9/fn7Fjx/Lhhx9WuOMvKCjgP//5DxMmTKh0rHfeeQcXFxfatWuHn58fGzZsYP369QD8/PPPfPvtt/zwww94eXlZ49jOnTtrdK53yy5TL/bo0UPWdnZOo9GIt7c327dv54EHHqi03mwy8+0/D3PhmOUlsAdi23LPoNa1WieAfV+f5cCGcwAMfb4LIV18a/2YiqI0DkKIg1LKqq+qNdBg7vhvJzMzk40bN1JcXEx2dja//e1vCQ4OpmfPnlWW3//dOWvQ79y/Jd0eblUn9bx3cDBuPpZfIT9+dpw8Q1GdHFdRFOW6RhP4TSYTr776Kj4+PoSGhpKSksI333yDo6NjpbIXTmRw4PtzALSM8Kb3mHbWly9qm6OTloETOoKA/JxiflpxQo3tVxSlTjWaMYW+vr7VmtwlN7uQzf93DKRluObASR1r/Dbu3WrR3pt7Brbm0ObznDucxvHdl+jYO6hO66AoStPVaO74q8NslmxZcox8QxEIeHhSR/SeN3/4W5vuH96GZi0tD6d3/uckWVfUKB9FUepGkwr8R3ekknLC8mJIjyEhtOrgY7e6aB01PDypI1oHDSVFZr79ZzzGzAK71UdRlKajWoFfCKERQgwRQiwUQuwUQpwSQlwQQsQLIb4QQswQQtT+kJi7YCox88smy4tT/iEe9BwaYt8KAc1auNG79MWu7Kv5rF/4CzkZKvgrTdOuXbt44IEH8PT0xMfHh+joaPbv3w/A0qVL0Wq11uGPoaGhTJo0iZMnT9boGLZKzbx3714efvhhfHx88PPzY/To0Vy6dKnC9r/88gt9+/bFzc2NgIAAa4qZ+pCa+ZaBXwjhLoR4DUgB/gsMAE6X/nkZsBPQAX8CzgohfhBC9KnNCt+pxL2XreP17x8eiqaK/Pn20KlvC6JHWRI0GdIKWL/wF1ISM0k+mk7i3kucPXRNPfxVGj2DwUBMTAzTp08nIyOD1NRU5s6dW+E9nF69emE0GsnOzmbLli3odDq6d+9OQkJCtY5hy9TMmZmZTJkyhXPnzpGcnIy7u3uFNPFpaWkMHjyYZ599lvT0dE6fPs2gQYOA+pGa+Zbj+IUQacAhYCnwtZSycpaisrL3AGOA3wCvSyk/uFnZuhjHX57ZZGblvH0YruXjH+zOqD/1qLNRPNV1ZFsKOz6v+u4lelRbug2s1z+oFOWuHDhwgIEDB5KVlVXl+qVLl/LJJ5+wa9euCstjYmJwcXHhyy+/vO0x/v3vf7N06VJ+/vlnwJK2wdfXl0OHDlV6gzc3Nxdvb28SEhKsWTrHjRtHixYtePPNNyvt+5dffqFfv37WRG6zZs3iwoULN03lfCMPDw9++umnCnmMqlJX4/gfkVI+LKVceaugDyClPCSlfAUIBrbcbcVs6dSBqxiuWXJ+dx8SUu+CPljeJej/dHuqqtr/vkkiN6uw7iulKHUkPDwcrVbLhAkT+P77761J2m4nNja2wluv9kjNDLBjx44Kb/bu3bsXHx8fHnjgAfz9/Rk2bBjnz5+vclt7pGa+ZeCXUh6s6Q6llHlSysqzG9iJNEsOliZF8wnSE1qP35SN7NOCCW9GM+rlHoz7cy/LLxONoLjQxO61KsWD0nh5eHiwa9cuhBBMnjwZPz8/hg8fzpUrt05jHhQUVGEKVnukZj58+DDz58/n7bffti5LSUlh2bJlLFq0iPPnzxMaGsrYsWMrbWuv1Mz1o6O7FiXFp5F5KReA7kOCEXU8Zr+m9J7OBIR64OGrIyDEgy6lM3id2n+FlMTq3QUpSkPUoUMHli5dSkpKCgkJCVy8eJEZM2bccpvU1FR8fKo3Oq82UjOfPn2aIUOGsGjRIvr0KXu8qdPpGDFiBD179sTFxYW5c+fy888/k51dlvbdnqmZbRL4hRBbhBBbbbEvWyoqKGHvV2cA8PTT0bZ7gJ1rVHP3DQ21Zgnd8flJTCY1b6/S+EVERDBx4sTbPrhdv359hYB7K7ZOzZycnMzAgQOZM2dOpcydXbp0qdClfP3P15+p2js1s63u+IUN92UTUkp+XH6czMuWF6PuGx5a52/o2oKTzoEHRlr6/jIv5XL4x5TbbKEoDc+JEydYuHChdY7dCxcusHr1aqKioiqVNZlMJCUlMX36dLZt28bcuXOrdQxbpmZOTU3lwQcf5IUXXmDq1KmVtp80aRLr168nLi6O4uJiFixYQO/evfHy8qofqZmllHX+6d69u6xtBzeekx88u1V+8OxW+eOK47V+vNpkNpvluncOyg+e3So/+d0OWVRQYu8qKYpNpaSkyNGjR8ugoCDp6uoqg4KC5JQpU2R2draUUsolS5ZIjUYj9Xq9dHV1la1bt5bjx4+Xx44dq7Cfjh07yhUrVtz0OJs3b5bt27eXLi4usl+/fjIpKcm67i9/+YscPHiw9Xt6erp87LHHpKurq2zVqpVcuXKldd28efMkIPV6fYVPeYsXL5ZBQUHSy8tLxsTEyPPnz0sppdy2bZsEpE6nq7Dtjh07bttOwAFpgxjcKNMyXziWwTf/iENKCAj1YMTv7kXrWK9+kNTYpTPZrHvb8qy995h2dH2wbrKJKopSf9hqOGe1k7QJIR69XRkp5Ya7q87dSU81cnTnRY7vuYSUoHN3ZPCUTg0+6AM0D/MkqJ0XF09lEbf5PJ36tkDr0PDPS1GUuleT7JzfAhJLf3555X8yaO+6RnfAkJ7P1qXHuXgqy7pMoxUMntIJN28Xe1SpVtw7OJiLp7IwZhZy8n9X6PBAc3tXSVGUBqgmgT+0imU+wCBgIjCpivW1TkpZIeg7OGkI7xlAlwdb0ayF2603bmBad/TBt5UbaReMHPohmYiowHo/PFVRlPqn2oFfSllVNqNk4JAQwgTMAobbqmLVdfrgVWvQ7/pgK3rGhODsWnnylcZACMG9jwTzwydHybycx9n4a4Td42/vaimK0sDYqpP4EPCgjfZVbcVFJn4ufaPVK8CVXrFhjTboXxd2rz+efjoADmw4R1F+iZ1rpChKQ3PXgV8I4YSlq+fSbYra3C+bkq0ZN3uPadckHnZqNMI6KXzaBSOr5+8j+Wi6nWulKEpDUpNRPfup+CAXwAkIAdyp4z5+Q1o+hzZZkh6FdPElOLJZXR7erjo80JzLSQZO/HwJY2Yh3/4jnvZRgUSPaovOzcne1VMUpZ6rycPdo1QO/AXAF8B/pZRVp62rBXmGIrYsPYapxIzGQVjz2TcVGq2Gh8Z3IOweP7avSsSYWUji3sucO5LGA7Ft6dCruXroqyjKTTW4F7jOHrrGTytPUGAsBixDHHs9HmbL6jUoRfkl7Fl/hoSdqdbLcmAbT+4fHkqLcG91AVCURsRWL3A1mMCfdSWP/RuSOLnPkqZVCLhnUHC9mk3Lni4nZbN9VSJpF8pm8XHzcab9fYG0jwrEO1Bvx9opimIL9SrwCyFeLd3XguqUr27gN5nMXDqdTfzWC5w7kma9o/XwdWHgxI40b+t1F7VufMwmM0e2p3JgwznrL6Lr/EM8aH9/IGH3+uHkYunhEwIcnOzyzp2iKHegvgX+ktJ9VSuK3CzwF+YVc+F4JhdPZ3Et2cC1C0ZMxWVpiLWOGjr2DiLqsTbW4KVUZioxc/5oOol7L5N0JA1zyc3/H3sHuhJ2rz9h9/rTrIW+Xs5OpiiKRb0K/DV1T9d75dert5BvLCI/p5i8nCIunc7i8lkD0ly5Pjp3Rzr3b0mnvi3QuatRKzVRkFvM6YNXSdx7mctns29Z1kXviM7dEZ27Ezo3R1yu/1fvWCHfkUYjcHDS4OCkxdFZi4OTFgcnDY5O5f7srEWjFepCoig21KADf3hwpHxxyD9uul7v6URAG0/8g93xb+1B83aeODiqLom7lXU1j6vnDFz/X15caCI5IZ0LxzIwlagJXhSlPhv5x+40D/Oq2+yc5QkhNECl7GdSyrxqba8tuwvUaAUubo54B+ppHelDcGQzfIJUl0Nt8PJ3xcvftcKyTn1bUFRQQnJCOllX8sjPKbb+Eiuw/rcYcxW/xBRFaZhq8gKXAP4ITKbqhG1Qzeycek9nnp4fhc7NESedgwrydubk4kC7HtWfltJsMlNSZKa4yERJkansz4Vlfy4uNGE2qYuFotiKezPbZRqudlePEOJFYB7wN+AvwJ8BE/Akljd435BSflrNfeUAiXdQ38bIF0izdyXqCdUWZVRblFFtUaa9lLLy7PA1VJPAnwD8G/gnUAz0kFL+Utrt8w1wREr5p2ru64At+qkaA9UWZVRblFFtUUa1RRlbtUVN3nwKBeKklCYsgd8LQEppBhYDE+62MoqiKErtq0ngTweuz2xyHrin3DpvQGerSimKoii1pyajenYDPYENwCpgnhDCBygCXgC21mBf/65B2cZOtUUZ1RZlVFuUUW1RxiZtUZM+/vZACynlj0IIZywPeUdhudPfDEyXUl61RaUURVGU2mOXF7gURVEU+1FpLRVFUZqYWwZ+IcRWIcTA6u5MCOEvhJgvhJh+91VTFEVRasPtHu5uAFYKIQqAdcDPQAKWlykKsQzpDAW6A0OAfsAm4Pe1VF9FURTlLt22j18IoQeeAsYD91M5LYPAMtH6OuBTKWXc7Q7q6+srQ0JC7qC6d89UUoxW62BJRq8oitKAHDx4ME1K6Xe3+6nRw10hhCvQFQjEkqQtA0iUUp6ryUHvZurFO1VSVMSWTxdzdNsWAtq05Yl5b+LobLvcF4qiKLXNVmmZa5SdszT75p67PWhdM2ak8/XCN7h02pIe6MrZ0/y45N88MvW3dq6ZoihK3WvUo3qk2cyZg/9jxSszrEHfzdsHgISffuDYzp/sWT1FURS7aJTzFxbkGjm6bStxP3xL1uVLloVC0PuJcXR7JIaVs2aSeSmVLR//k4A2bWnWopV9K6woilKH7PICV2328edlZ7Fi1kxy0q5Zl7k18+Xh37xAm3t7AnD13FlWvfp7TMXFePoH0O7+aPxahxDQph3NWqqLgKIo9ZOt+vjrRVePm5ub9ePo6Iijo2OFZQCff/45ffr0wcPDAweHm/9Q+eX7r61Bv1XHzgz73StM/sen1qAP4B/ShgcnPgtA9tUrHPhmHd//8+8s/f1z7Fi5BPU2s6LUTEhICE5OTqSlVUyb361bN4QQnDt3zj4VK2fr1q1ERETg6urKgAEDSE5OvmnZjIwMRowYgV6vJzg4mFWrVlnXnTt3DiFEhRi1YMGCSvsoKioiIiKCli1bVnmM7du3I4Tg1VdfvfuTq6F60dVjNBqtf/7Nb35DSUkJS5curVDG29ub559/nvz8fKZMmVLlfory84j74TsAwu+PZtjvXrnpMTs/9AgIOLl3N2nnz5GblQnA/q/X4qx34/7HR9/lWSlK0xIaGsrq1auZPt3y/uaRI0fIz8+3c60s0tLSiI2N5ZNPPmHYsGHMmTOHJ554gr1791ZZ/oUXXsDJyYkrV64QFxfH0KFD6dq1K5GRkdYyWVlZt7wJffvtt/H3968Q364rLi7mxRdf5P7777/7k7sD1b7jF0L0EUI8Vu67rxBilRAiTgixUAjhWDtVtHjkkUcYO3Ysbdq0qXL9xIkTGfPYcApzcwHo+dgoQkJCWLFiRZXlhRB0eWgwo2YvYOpHn/Hsh8vwC7bMKLlr9TIOb91UOyeiKI3UuHHjWL58ufX7smXLGD9+fIUyhYWFvPTSS7Ru3ZqAgACmTp1qvThkZmYSExODn58f3t7exMTEkJKSYt22f//+zJkzh+joaNzd3Rk0aFClXxg3s27dOiIjIxk9ejQuLi7MmzeP+Ph4Tpw4Ualsbm4ua9euZcGCBbi5udG7d2+GDx/OZ599Vu22SEpKYsWKFbzyStU3nwsXLmTQoEFERERUe5+2VJOunr8Bncp9XwQ8BOwFJgKv265aNWc2m0m7YPnp1rpTVwLD2tVoezefZoycNR/PgEAAtnz8T07972eb11NRGquoqCgMBgPHjx/HZDKxZs0annnmmQplXn75ZU6ePElcXBynT58mNTWV+fPnA5Z/w5MmTSI5OZnz58+j0+mYNm1ahe1XrVrFkiVLuHr1KkVFRbzzzjvWdV26dKnQJVPe0aNH6dq1q/W7Xq8nLCyMo0ePVip78uRJtFot4eHh1mVdu3atVDY4OJiWLVsyadKkSheg6dOn88Ybb6DTVZ6mJDk5mf/7v//jtddeq7KudaEmgb89cBCsL3KNAF6UUk7FMgn7E7avXvVlX71McWEhYLnbvxN6L29GzVqAq6cXUprZuPhd8gzZtqymojRq1+/6N2/eTEREBC1atLCuk1Ly8ccf8+677+Lj44O7uzuzZs3i888/B6BZs2aMHDkSV1dX3N3dmT17Ntu3b6+w/0mTJhEeHo5Op2PMmDHExcVZ1x0+fJinnnqqynoZjUY8PT0rLPP09CQnJ6fGZX19fdm/fz/JyckcPHiQnJwcnn76aWvZ9evXU1JSwogRI6qsy29/+1vrrwl7qUkfvxNQUPrn6NJtvyv9fhJobsN61Yg0m0k7b7nb9w8NI7hztzvel1dgc0b88TVWvvp7ivLz2bvuc+uDYEVRbm3cuHH07duXpKSkSt08165dIy8vj+7du1uXSSkxmUwA5OXlMXPmTDZu3EhmpuWZW05ODiaTCa3WkikmMDDQuq2rq2uV/edVcXNzw2AwVFhmMBhwd688b/ntyrq5udGjh2VgTUBAAB988AHNmzfHYDCg1Wr54x//yIYNG6qsxzfffENOTg5PPGHX++Qa3fGfAAaX/vlpYI+U8vrlMghL+ga7OPHzDigupKjExH2PjUIIQUlJCVev3tm8MIFtw4ns+xAA8T9sIPPyRVtWV1EareDgYEJDQ9mwYQOxsbEV1vn6+qLT6Th69ChZWVlkZWWRnZ1tDd4LFy4kMTGRffv2YTAY2LFjB4BNRtlFRkYSHx9v/Z6bm8uZM2cqPKy9Ljw8nJKSEk6dOmVdFh8fX2VZsDwvvF7PU6dOce7cOfr06UNgYCCxsbFcunSJwMBAzp07x9atWzlw4ACBgYEEBgayZs0a3nvvPR577LEq911bahL45wMzhRDXsCRte7PcusHAIVtW7EYmk4mCggKKiooAKCgooKCggIunTvDDv/9BK29PzqRl4ODfnMLCQmbPnk1xcfEdH++BMU/j4OiE2WRi1+fVf6ijKE3dp59+yo8//oher6+wXKPRMHnyZGbOnGm9KUtNTWXTJstAipycHHQ6HV5eXmRkZPD667Z7bDhixAgSEhJYu3YtBQUFzJ8/ny5dulT5cFWv1xMbG8trr71Gbm4uu3fv5quvvmLcuHEA7Nu3j8TERMxmM+np6fz2t7+lf//+eHp60qlTJy5cuEBcXBxxcXF88sknBAQEEBcXR6tWrViwYIH1GUdcXBzDhw9n8uTJLFmyxGbnWh3VDvxSyq+BDsBUoJOU8vtyq/cAf7Zx3Sr47LPP0Ol0PPLII5hMJnQ6HTqdjk9f/SMlhYXcF96GYTEx9OjRk7CwMFq3bl2hf7GmPHz9uPfR4QCc3LOTS6cSbXUqitKohYWFWbtCbvTWW2/Rtm1boqKi8PDwYODAgSQmWv5tzZgxg/z8fHx9fYmKimLw4MFV7uNmIiMjWblyZZXr/Pz8WLt2LbNnz8bb25t9+/ZZny0AvPHGGwwZMsT6ffHixeTn5+Pv78/YsWP58MMPrXf8Z8+eZfDgwbi7u9OpUyecnZ1ZvXo1AA4ODta7+cDAQHx8fNBoNAQGBqLVanF3d6+wXqfTodfr8fHxqdG53q2azLk7HvhOSplexTofIEZKubzylpXZ4s3dvOwsVs/5A1lXLiE0Gkb88TVC77nrF9oqKMzL5ZPfTqYgx0DLDp0YM/ev1p91iqIodc0eb+4uAcJusi60dH2dkFLy7aK/kXXFkofn4SnTbB70AZxd9fSKtTyESTmewPo355GTXr1xw4qiKPVVTQL/rW51mwGGW6y3qQtHj3Dh6GEAomKfoPOAQbV2rK6DHqV5u/YAJMUdZNlLL5Dw02aV1kFRlAbrlsM5S9/ULf+4eU7pw93yXIA+wH4b1+2m9q1fA4De24f7R9TusCitgyNPzHuTvevWsG/9fyjMy2XTvxZhzMwgKta+Q7IURVHuxO3u+P2BzqUfsHT1dL7hEwz8ANTJYPeLJ09wPsEyLKtHzAgcnJxq/ZhaB0eixzzD03/5O81atgZg33//Y83voyiK0pDcMvBLKT+WUvaUUvYEtgOx17+X+0RLKX8tpUyqiwpfv9t3cfeg68AhtyltWwFt2jLsd68ghIaSwkL2rltTp8dXFEWxhZoM5xwgpayc0agOXT13lrO/WHqUuj/6GI4udT9nbrMWrYjsb3m56/CWjWRfvVzndVCUxmjXrl088MADeHp64uPjQ3R0NPv3W/69L126FK1Wa02DHBoayqRJkzh58mSNjmGr1Mx79+7l4YcfxsfHBz8/P0aPHs2lS5es69977z3atGmDh4cHQUFBzJw5k5KSEuv6uLg4+vTpg6enJy1btrTmK6orNcrHL4QIEkJMEULMF0L87YbPW7VVyev2rf8PAE46V7o9MrS2D3dTvUaNRevggNlUwp4vV9utHorSWBgMBmJiYpg+fToZGRmkpqYyd+5cnJ2drWV69eqF0WgkOzubLVu2oNPp6N69OwkJCdU6xvXUzAsWLCAjI4MePXrcMnVC+dTMK1eu5LnnnrMmasvMzGTKlCmcO3eO5ORk3N3dmTRpknXbYcOG8csvv2AwGEhISCA+Pp7333/fuv6pp56ib9++ZGRksH37dj788EO+/vrrmjbbnZNSVuuDJSlbAVAMpAJJN3zOVndf3bt3lzWVnnpBvvNEjHxnzFC5c/WyGm9vaz8u+Ui+M2aoXPjEMJl2Idne1VGUBm3//v3S09PzpuuXLFkio6OjKy0fOnSoHDlyZLWO8dFHH8levXpZvxuNRuni4iKPHz9eqazRaJSOjo4yMTHRuuyZZ56RL7/8cpX7PnjwoHRzc6tyXVpamnzooYfkc889Z12m0+nk0aNHrd9HjRol33jjjdueA3BAVjPO3upTkzv+N7A8xA2QUraQUobe8Kk6Ub6NHPhmHUiJg5Mz9z5at3ktqnL/iDE4OrsgpZltn32KqdzPOEVRaiY8PBytVsuECRP4/vvvrUnabic2NpadO3dav9d1aubrduzYUSmXz6pVq/Dw8MDX15f4+HiefbZs/MuMGTNYvnw5xcXFJCYmsmfPHgYOHFitc7aFmgT+VsD7Uso6T8aWm5XJsR0/AtBpwEBcPTxvs0Xtc/X0ovtQywXoXNxB1r/1OgW51csUqChKRR4eHuzatQshBJMnT8bPz4/hw4dz5cqVW24XFBRERkZZSKqr1MzlHT58mPnz5/P2229XWP7UU09hMBg4efIkU6dOJSAgwLouJiaGL7/8Ep1OR0REBL/+9a/p2bPnjbuuNTUJ/D9jyclf5375/mtMJSUIoaH70KpzXNvDfSPGWOfyTT58qDSFhHrYqyh3okOHDixdupSUlBQSEhK4ePEiM2bMuOU2qamp1c5zY8vUzNedPn2aIUOGsGjRIvr06VPlcdu1a0dkZCTPP/88YHloPHjwYF577TUKCgq4cOECmzZtYvHixdU6D1uoSeD/HTBFCDGh9CGv642f2qhgUX4e8Zstua3bRUXjFRB4my3qjqOTM4/94VXrnX9G6gVW/OlFtnz6ISknjiLNZjvXUFEapoiICCZOnHjbB7fr16+/acC9ka1TMycnJzNw4EDmzJljzdx5MyUlJZw5cwawJHnTarWMHz8eBwcHWrZsyZNPPnnTHP61oSaB/zCWF7aWABeAnCo+Nnfkx81l8+gOi71N6bqn0WjpP34yA3/zPEKjoTAvl/gfvmPN3Jf55Le/IflInL2rqCj13okTJ1i4cKF1jt0LFy6wevVqoqKiKpU1mUwkJSUxffp0tm3bxty5c6t1DFumZk5NTeXBBx/khRdeYOrUqZW2/+STT6ypp48dO8Zf//pXHnrIMgw8PDwcKSWrVq3CbDZz+fJl1qxZU+H5Q22rSeD/FTCp9POrm3xsylRSwsHv/gtAq8guNZ5Hty51ffhRnvrzQjoNGISzqyUPueHaVb55968Yrt3ZhDCK0lS4u7uzb98+7r//fvR6PVFRUXTq1ImFCxday+zZswc3Nzc8PDzo378/BoOB/fv307lzZ2uZukrN/Mknn3D27Flef/1167sF5adS3L17N507d0av1/Poo4/y6KOP8sYbbwCW5xnr1q3j3Xffxdvbm27dutGpUydmz55tm8ashmqnZbal6qZlPrbzJ77/wPI/PvZP82olA2dtKCku5tS+3Wxc/C5mk4mg8A48Me9NNKXTxymKotwJe6RlrlMlxcX8/J8VAPi2CiakW/fbbFF/ODg60qF3f6KfsPwsvHjyOHvWqhe9FEWpH6od+IUQ14QQV2/1sWXF4jZ+Q/ZVy1Cu3mPHN8gJUHoOiyW4yz0A7F23xppKWlEUxZ5umZb5Bv8EbuwX8gEeBDyAT21VqfwcA3tLk7G1iuxCm3vvs9Wu65TQaBjywu9Y9odp5BuyWffm60T2H0iPoY/jFdjc3tVTFKWJqnbgl1LOq2q5sNyK/wew2aure9d+bhnJIwT9nvlVg7zbv07v5c2jL/yO9X9bQElRIfE/fMfhzd8T3qs3D09+wfogWFEUpa7cdR9/af6IT4Bpd18dyLyUStwP3wHQsc8AAtq0tcVu7SqkW3cmvPMBnR96BK2jI1KaSfx5B98t+htms8ne1VMUpYmx1cPdNoBNZkTZsXIJZpMJB0cn68PRxsAnqCWDpkxn8gf/R8c+AwDLVI67Pv/MzjVTFKWpqXZXjxDi+SoWOwEdgKeBL+62Mqf27+H0/r0AdI8ZgYev393ust7Re3nzyHMzMGamcz7hMPu/+hK/4FA6RPezd9UURWkiavJw94MqlhUCKcBi4PW7qUhBrpGtn34IgGdAIPc/PvpudlevabRaYmb8iZWzZpJ99Qo/fLgIV3dPgrt0s3fVFEVpAmoyA5emio9OStlOSvlHKWXu3VRkx8ol5GZasuw9PHmaXWbXqks6dw8e+8McHJ1dKCku4su/vMp/5s/iwrEj9q6aoiiNXE3u+GvN+YTDHNm6CYBOAwYR3LmbfStUR/xahzDsd6/w/T//Tr4hmwtHD3Ph6GHcfJqhc/fAxc0dDz9/ug0aWq/TVSiK0rDUKGWDEKIN8AegN5Yx/BnATuAdKeXZ6u6nfMqGnIw0/jPvFbKuXELv5c3EhR/iUi7nRVNQXFBA/Jbv2f/1WvKys6osE3pPD6Jin8S7eRAlRUWUFBXi6umlhoMqShNiq5QN1Q78QojuwE9Ypl/8FrgCBABDARdggJTyl+rsq0ePHnL//v0c3baFbcs/oTDP0ks0/PezaHffAzU/i0aiuKiQxN07yLx8kQJjDgU5OSQnxFmzk1YiBP7BbWjZsRNB4R3w8PPDzacZek9vlRdIURohewT+n7A8Exgipcwrt9wV2ACYpZQPVmdf93brJl998nGSDlnu+oVGwwOjnyYq9uYTHzdVhXm5HNr4LQe/+y8FxuplvhYaDc1atiYwLJzmbcPReXhQYDRSYMyhqCC/Qlmt1gEHJyccnJxBCEqKCikpLMRsNuOi1+Pi5o6z3g2tQ1mvoNlkspQrKsJkKsHBwREHJ2ccnJwQmnqb/klRGrTAtuHo3NzrPPDnAmOklN9VsS4GWCOlrFa/Q0iAn5ze/34AfFuHMPi5GY3iRa3aVJSfR/KROKTZjIOTMxoHBzIvppByLIELxxPIN2Tbu4qKotSisQvepkX7jjYJ/DV5uJsPNLvJOh8sXUDV4ubTDKHRcP/jo4ka+SRaB8caVKNpctK5VuoGC+lyD/cMHoaUktzMDIwZ6eRkppNz7SpXzp7m0plTZF5MsZZ3dNHhpNOVpcCQkpKSEuvd+/XJ7K/fuRfm5mI2qUnkFaWxqUng/w54UwhxVkq56/pCIURv4K/AN9XdkYvejUnv/gvvwKAaHF65GSEEbj7NcPNpxo0TUxbm5VJSVISLm9stL7DXf/mVz4skpaS4sIACY06FaSSFRouDkxOOTs5oHLSYiotLHzgXIaWablJRaoPeq3pzC1dHTbp6mgFfAb2Aa1ge7vqXfn4GHpdSpldzXzlA4p1UuBHyBdLsXYl6QrVFGdUWZVRblGkvpaw8O3wN1XgGLiHEYKAn0By4BOyTUv5Qw30csEU/VWOg2qKMaosyqi3KqLYoY6u2qPELXFLKjcDGuz2woiiKYh+3HHsnhGgmhFgrhHjkFmUeKS3jb/vqKYqiKLZ2u0HXM7CkXL5VV84PQCjw+xoc9981KNvYqbYoo9qijGqLMqotytikLW7Zxy+ESAT+LqX86JY7EeJZYKaUMsIWlVIURVFqz+3u+IOBY9XYz3Eg5K5royiKotS62wX+fCwTqd+OW2lZRVEUpZ67XeD/BRhejf08VlpWURRFqeduF/j/CfxaCDHhZgWEEOOBSVQ9Q5eiKIpSz9z2BS4hxEJgJnAQy/j984AEWgOPAD2Ad6WUL1X3oL6+vjIkJOQOq1w7ZIkZc34JwkGDcNQgHFSWSUVR6peDBw+mSSnvejLy277AJaX8vRBiG5ahnS8BzqWrCoHdwGNSym9rctCQkBCuT8RSHxRdNHLt4yPI/LKEZMJZi0s7L/Q9A3Fu543QiFvsQVEUpfYJIZJtsZ9qvbkrpfwG+EYI4UBZhs50KWWDT91YfDmXtE8qBn0AWWgiPyGd/IR0tJ7OOLf1QpaYkUUmEAL9fYHoImyXNElRFKWu1ChlQ2mgv1JLdalzxVdyufbxEcx5JaARNHs6AscWbhSn5lJ43kDeoauYDUWYsgvJO1jxtAuOpaPr7IvXsDC0Hk52OgNFUZSaqxeTrdc1WWzCuOcShp8uWO70NeAzNgJdpC8ADl4u6CKb4TkohIJTmeQdvEJJej7CSYvGWUvx5TxM2YXkH0mj4FQmXsPC0HcPsPNZKYqiVE+9CPxu5SZXLywsBMDZ2dm6zGg08vLLL/Ptt99y4cIF3NzcGDp0KG+99RY+PtXvbpFSknfgCoYtyZiyiywLBfg8EYFrZ99K5YVWoIvwqdSlYy4yYdicjHFXKrLAROYXJym5mofHIyHqWYCiKPVevRi6YjQarZ8JEybw9NNPV1gGoNVqWbFiBenp6cTHx5OSksKkSZNqdJy8g1fJXHvKGvSdwzzxf6Ebrl1r9pBc46TFa2gb/Kfdg4OfDoCc7SlkrElElqiJSJSmJyQkBCcnJ9LSKqbN79atG0IIzp07Z5+KlbN161YiIiJwdXVlwIABJCff/DlpRkYGI0aMQK/XExwczKpVq6zrVq5ciZubm/Xj6uqKEIKDBw8C8N5779GmTRs8PDwICgpi5syZlJSUPUP8+eefue+++3B3d6dLly7s2rWr0vFrW70I/NXxxhtvcM899+Do6Iifnx/Tpk1j27Zt1vUTJ07kN7/5TYVtQkJCWLFihfV77sHLAGi9nfH9dSd8f9MZp5Z3PqeBUws3/J/rilOI5eXm/PhrXPvkCKbswjvep6I0VKGhoaxevdr6/ciRI+Tn148X+tPS0oiNjWXBggVkZGTQo0cPnnjiiZuWf+GFF3BycuLKlSusXLmS5557jqNHjwJUujFdvHgxbdq04d577wVg2LBh/PLLLxgMBhISEoiPj+f9998HLBeU4cOH84c//IGsrCz++Mc/MmzYMDIzM2u/EcppMIH/Rlu3bqVLly7VLm/KLabonAEAt15BuLTzrjDN4J3SuDri9+vO6Ep/NRSdM3D53V/IPXSVmk5yoygN2bhx41i+fLn1+7Jlyxg/fnyFMoWFhbz00ku0bt2agIAApk6dar04ZGZmEhMTg5+fH97e3sTExJCSUjZndP/+/ZkzZw7R0dG4u7szaNCgSr8wbmbdunVERkYyevRoXFxcmDdvHvHx8Zw4caJS2dzcXNauXcuCBQtwc3Ojd+/eDB8+nM8++6zKfV8/z+vxJCwsDC8vL8DSvazRaDh9+jRgudsPCAhg9OjRaLVannnmGfz8/Fi3bl21zsNWGmTgX7t2LR9//DGLFi2q9jYFJzIsr50BLh1sOwxTOGrweaI9HgNbgwZkQQmZaxJJ/+w4ptximx5LUeqrqKgoDAYDx48fx2QysWbNGp555pkKZV5++WVOnjxJXFwcp0+fJjU1lfnz5wNgNpuZNGkSycnJnD9/Hp1Ox7Rp0ypsv2rVKpYsWcLVq1cpKirinXfesa7r0qVLhS6Z8o4ePUrXrl2t3/V6PWFhYda7+PJOnjyJVqslPDzcuqxr165Vlk1OTmbHjh2VLnCrVq3Cw8MDX19f4uPjefbZZwHLheDGG0IpJQkJCVXWu7Y0uMD/xRdfMHnyZL7++mvrT6vqKDhmmQ7YwU+Ho5+rzeslNAKPgcH4P9cNB3+d9ZjpK44jzerOX2kart/1b968mYiICFq0aGFdJ6Xk448/5t1338XHxwd3d3dmzZrF559/DkCzZs0YOXIkrq6uuLu7M3v2bLZv315h/5MmTSI8PBydTseYMWOIi4uzrjt8+DBPPfVUlfUyGo14enpWWObp6UlOTs5dlV2+fDl9+vQhNDS0wvKnnnoKg8HAyZMnmTp1KgEBllF/DzzwABcvXmT16tUUFxezbNkyzpw5Q15eXpX1ri0NKvAvWbKEZ599lm+++YYBAwZUWOfm5kZubq71e0lJCVevXgVAFpspOGXpQ9N1bEZtcmrlTsD0e9HfHwhAUVI2ufsu1eoxFaW+GDduHKtWrWLp0qWV7oKvXbtGXl4e3bt3x8vLCy8vLwYPHsy1a9cAyMvL49lnnyU4OBgPDw/69u1LVlYWJpPJuo/AwEDrn11dXa2DP27Hzc0Ng8FQYZnBYMDdvfIzvpqUXb58ORMm3DSVGe3atSMyMpLnn38esFzcvvrqK/7+978TEBDAxo0bGThwIC1btqzWedhKgwn877//Pi+99BKbNm0iOjq60voePXqwdetWkpKSKCwsZPbs2RQXW7pZCs5kIYsso21cajnwg6Xrx2t4WxxbWf6iZG9IoiS9fjzkUpTaFBwcTGhoKBs2bCA2NrbCOl9fX3Q6HUePHiUrK4usrCyys7OtwXvhwoUkJiayb98+DAYDO3bsALDJs7LIyEji4+Ot33Nzczlz5gyRkZGVyoaHh1NSUsKpU6esy+Lj4yuV3b17NxcvXmTUqFG3PHZJSQlnzpyxfu/Xrx/79+8nIyODzz77jMTERO677747PbU70mAC/4svvojBYGDAgAEVhlJd9/TTTzN8+HDuvfdewsLCaN26tfVn5vVuHo2bI06t7nwUT00IrcBnVDvQCmSxmcy1p1SXj9IkfPrpp/z444/o9foKyzUaDZMnT2bmzJnWX+Opqals2rQJgJycHHQ6HV5eXmRkZPD666/brE4jRowgISGBtWvXUlBQwPz58+nSpQsREZUnDdTr9cTGxvLaa6+Rm5vL7t27+eqrrxg3blyFcsuWLWPkyJGVfgl88skn1vM7duwYf/3rX3nooYes6w8dOkRxcTEGg4GXXnqJli1b8sgjN53WvHZcf9hQl5/u3bvLumI2mWXqn/fICy/vkBlfnqyz416X/dN5eeHlHfLCyztk5tenZW78VZl3PF0WXcmt87ooSm0JDg6WmzdvrrS8uLhYAjIpKUlKKWV+fr585ZVXZGhoqHR3d5cRERFy0aJFUkopU1NTZb9+/aRer5ft2rWT//rXvyQgi4uLpZRS9uvXT3788cfWfS9ZskRGR0dbv3fs2FGuWLHipnXcvHmzbN++vXRxcZH9+vWz1klKKf/yl7/IwYMHW7+np6fLxx57TLq6uspWrVrJlStXVthXfn6+9PT0lFu2bKl0nIkTJ0p/f3/p6uoqg4OD5UsvvSTz8/Ot65988knp4eEhPTw85JgxY+SVK1duWucbAQekDWLwbdMy14YePXrIusrOWXjewLXFlp94zSZ0RNeh9rt6ypMmydUP4yhOqdwX2WxipEr0pihKtQkhDkope9ztfhpMV8+dMOcVk3fAklxNOGpwaetV53UQWoHPmPZoPZ0rrcv65ox601dRlDpXL3L12JI0S4w7U8iLu0bx5Vzr2H3ncG+Eo9YudXL0d6X5K/chTWZLuufETDLXJGJKL8D480Xc+9btE31FUZq2RnfHb9yZSvb35yi+VBb0NW6OuPdpcesN64DQatC4OuLazQ/nNpZxwoat5zEZi+xcM0VRmpJGdcdfeN5A9qZzADgGuqK/vznOoZ44+LvWq6yZQgg8h4Vx9f1fkIUmDD8k4x3bzt7VUhSliWg0d/zm/BIyVp8As0S4ONBsQiRuvYJwDNTXq6B/nVNzPfr7LC+j5O6/TFFq9V5EURRFuVuNIvBLKclcdwpTpiUrps+odjh4u9i5Vrfn8XAwwlkLEq79K57sjecw5zf42SwV5Y7s2rWLBx54AE9PT3x8fIiOjmb//v0ALF26FK1Wa31/JzQ0lEmTJnHy5MkaHcNWqZmPHTtGjx498Pb2xtvbm4EDB3Ls2DHr+p9++okBAwbg6elJSEhIpX3HxcXRp08fPD09admypTVfUV1p8IFfmiTZ3yWRf8SSpU/fqzm6TpUnVamPtG5OeA0PsyR2KzaTs+0Cl/62H+Pei/aumqLUKYPBQExMDNOnTycjI4PU1FTmzp1bYUKmXr16YTQayc7OZsuWLeh0Orp3717tBGe2TM0cFBTEl19+SUZGBmlpaQwfPpwnn3zSuq1er+dXv/oVb7/9dpX7fuqpp+jbty8ZGRls376dDz/8kK+//rpa52ELDTrwm3KKuPbJEYy7UgFwbK7H69E2dq5Vzei7BxAwszu60hnAZH4JWf89g1Hl91GakOt37mPHjkWr1aLT6Rg0aFCVqde1Wi1hYWEsXryYfv36MW/evGodw5apmb28vAgJCUEIgZQSrVZrTb0McN999zFu3DjatKk6Hp07d46nn37aei69e/euMvtnbWmwgb/wvIEr7x+iKCkbsAzX9P1NZ4RjwzslRz9Xmj3dAf9p3XDwtWT2zPrqNAWJGXaumaLUjfDwcLRaLRMmTOD777+v9sQksbGx7Ny50/q9rlMze3l54eLiwvTp05k1a1a16gwwY8YMli9fTnFxMYmJiezZs4eBAwdWe/u71fCiJJYXs9KXHcWcUwQCPAa2xndiJFq9o72rdlecWrrjOykSjd4BzJC+6gRFl3Jvv6GiNHAeHh7s2rULIQSTJ0/Gz8+P4cOHc+XKlVtuFxQUREZG2Q1SXadmvp5o7oMPPuCee+65ZV3Li4mJ4csvv0Sn0xEREcGvf/1revbsWe3t71aDDPzZm85hzrU8BG02riMeA4Pr5cidO+HQTEez8ZHgIJCFJtKXJlCSpaZyVBq/Dh06sHTpUlJSUkhISODixYvMmDHjltukpqbi41O9tCe1lZpZr9czdepUxo8fb03OdisZGRkMHjyY1157jYKCAi5cuMCmTZtYvHhxtc7DFhpc4C+6kEPu/yxz5+p7BtZ6fn17cA72wGdMewBM2UVc+/dhSjIK7FwrRak7ERERTJw48bYPbtevX0+fPn2qtc/aSM18ndlsJi8vj9TU1NvW4+zZs2i1WsaPH4+DgwMtW7bkySefZMOGDdU6D1toUIFfmiWZ/z0NEjSuDngMDrF3lWqNaxc/vIZZHgyZMgq49tFhStJUTn+lcTpx4gQLFy60zrF74cIFVq9eTVRUVKWyJpOJpKQkpk+fzrZt25g7d261jmHL1MybN2/m0KFDmEwmDAYDv/vd7/D29qZDhw6A5UJQUFBAcXExUkoKCgooKrK8oR8eHo6UklWrVmE2m7l8+TJr1qyp8PyhtjWowJ+77xLFpS86eQwOafB9+rfjFt0Cr8fbAmDKLuTqR4cpOJWp8vorjY67uzv79u3j/vvvR6/XExUVRadOnVi4cKG1zJ49e3Bzc8PDw4P+/ftjMBjYv38/nTt3tpaJjIxk5cqVVR7Dz8+PtWvXMnv2bLy9vdm3b5912keAN954gyFDhli/L168mPz8fPz9/Rk7diwffvih9Y4/KyuLsWPH4unpSVhYGKdPn2bjxo24uFjeH9qxYwc6nY5HH33UOn/woEGDAMvzjHXr1vHuu+/i7e1Nt27d6NSpE7Nnz7Zdg95Gg0nLbDIWcfmdA8gCE06t3PF7rmuj6de/ndwDl8lce6os95C7I7pOvrjd3xzHQP2tN1YUpdFocmmZjbsuIgtMIMDr8bZNJugD6HsE4jM2Ao2rJbWSOaeY3D2XuPKPQxSczrJv5RRFaXAaROA355dg3GN5m1XXyRenFm632aLxce3iR/NZ99NsUiSu9/pb3lcwSdI/O0bRRZXnR1GU6msQgd+45yKy0ASAe/9Wdq6N/QgHDbr2PviMaY/vrzpZh3ymLUlQo34URam2ep+W2VxksqZkcGnv3STv9qviHOpJsycjSF95HHNOMdc+OYJLe280zlqEswOOzfU4h3igcan3/4sVRalj9T4q5O67jDnP8rKW+4Ot7Vyb+kXXyRevx9qS9d/TmDIKyN1zQ34fYclf5NrNH7c+LRCi6TwXURTl5up14JclZnJ2Wsb1OrfxxDnYw841qn/copojtILc/ZcxF5iQRSbMecXIIjNIKL6YS/bFJGSxGY+H1IVTUZR6GvilSVKQmIFxz0XMBstLD+4Dmm7f/u3oewai7xlo/S7NkuIreRQlZZP7v8sUX87FsDkZxwDXBpOyWlGU2lPvAr9x3yUMW5Ix5xRblzmFeuDc1st+lWpghEbg1FyPU3M9us6+XP3gEKbsIjLWJOLn44JTkHpOoihNWb0a1WP8+SJZ609bg77Wy9mSeXNCpOqfvkNadyeajY9EOGqQxWbSlx1TI4AUpYmrN3f8uQevkPX1GQAcfHV4PRaGc5hXk3pRq7Y4tXDDe3Q4GatOYMou5Mq7B/F4OBi36BYIrUCWmC3vApgljs31aJzrzV8LRVFqQb34F56fkEbmWssMPFpPZ3x/0wkHr/o/Z25D4trFD3NOEVkbLA96szckkXfoKhpXB4rO5yCLzZaCwpIa2sFPhzRJZKEJWWLGMVCPa1c/nMM8Edp69UNRUZQaslvglyVm8o9nkLv/MoWnMi0ZN/WOKujXIrfoFji38yZz3SmKzhkormqSFwklafmVMoEWpxrJO3gFjasDji3dkUWWEUSYwbGFG86hnjiHeICDwGQoskySoxE4NndD6+lUra46KSWYJMJBXVgUpTbZJUnbPeFd5Hfj/o05t+wBrnBxwG9KZ/XgsQ5IsyT3wGVy911Go3PAOcQDp1BPNE5aii4aKU41UpJZgHDUonHSgBAUnMrEbCy+/c6roHG1vFCmcXVEOGnROGsxF5RgyrFcIEx5JZZfFkUmyw2AmyMOvjocfHVodDe5N5FYhq5e306I0pfXtJYLh+ohVBoZt15BODbT2SRJm13u+GWJ2Rr0tc1c0PcIQN8jEK27kz2q0+QIjcDtvua43de80jqnVpVnGALLxaLwbDb5CWmYsgutAVyWmClKNlCSfvMHxua8EgrPZFe7fmZjMUXGYorOGW5fWFGaCFsOxbZL4NfqHXHt5odrz0CcQz3VA9wGQGgELm29cLnJsFqToZCi8zmgFWjdndB6OCGLLA+Niy8aKb6ajywowVxkQhaaEM5aazmN3tF6t44QmDILKUm3dDdZnz1UVScnTdkFSGL51VD6TEJRGhtbdoHaJfALZy0+T1ae9UZpuLQezug6OVda7uCrgy5+dqiRoig3Y5c+fiFEDpBY5weun3yBNHtXop5QbVFGtUUZ1RZl2kspq+6PrQF7jepJtMUDisZACHFAtYWFaosyqi3KqLYoI4So2dSFN6HGzSmKojQxKvAriqI0MfYK/P+203HrI9UWZVRblFFtUUa1RRmbtIVdHu4qiqIo9qO6ehRFUZoYFfgVRVGaGBX4FUVRmphaCfxCiOeFEElCiAIhxEEhRJ/blO8shNguhMgXQqQKIV4TjWTmlZq0hRDCRQixVAhxWAhRLITYVodVrXU1bIv+QoivhBCXhBB5pW3yq7qsb22qYVt0FEL8JIS4Ulr+rBDiDSFEg09uVdNYUW67dkKIHCGEsbbrWFdq+HciRAghq/gMrtbBpJQ2/QBPAMXAZKAD8A/ACLS+SXkP4DLwH6ATMBLIAX5v67rV9ecO2kIP/AuYAvwX2Gbvc7BjW8wC/gxEA22A54AS4Cl7n4sd2qItMBHoCgQDw4ErwN/sfS512Q7ltnMCDgLfAUZ7n4ed/k6EABJ4BAgs93Gq1vFq4QT2AR/fsOwU8NeblH8OMAC6csteBVIpHXXUUD81bYsbyn3QyAL/HbdFufL/Adba+1zqSVv8Hdhj73OxRzsA7wJLSi+GjSXw1zRuXg/8Pe7keDbt6in96dkd+OGGVT8AD9xks17ATill+Zk/NgFBpSfXIN1hWzRKNmwLDyDTVvWyB1u0hRCiLTAY2G7b2tWdO20HIcRQIAb4be3Vrm7d5d+JdUKIq0KI3UKIUdU9pq37+H0BLZafoeVdwfIzpCqBNyl/fV1DdSdt0VjddVsIIWKAh2j4L/PccVsIIX4WQhRguRPchaU7rKGqcTsIIZoDHwPjpJQ5tVu9OnUnfyeMwEvAGOBRYCuwRgjxTHUOWFtJ2m58K0xUsex25ata3hDVtC0asztqCyFENLAK+K2U8n+1UTE7uJO2eAJwx9LX/zbwMvBX21etTtWkHVYAH0op99Zuleym2m0hpUwDFpZbdEAI4Qv8EUs73ZKtA38aYKLyVcqfylez6y7fpDy32KYhuJO2aKzuuC2EEL2BDcBrUsoPa6d6deqO20JKeaH0j8eEEFrgEyHE21LKEttXs9bdSTs8CPQTQswt/S4AjRCiBHheStlQfw3aKlbsAyZVp6BNu3qklEVYnrY/fMOqh4Gfb7LZHqCPEMLlhvIXgXO2rF9dusO2aJTutC2EEH2B74HXpZTv1VoF65AN/15osNy4aW1UtTp1h+3QGehW7vMakF/65y9sX8u6YcO/E92AS9U9qK2fTj8BFAG/wTIsaRGW/qjg0vV/BbaWK++J5a7/cyzDOWOxjPJpLMM5q90Wpcs6lv4P/Bw4UPrnbvY+Fzv8vegP5GLp0ig/XM3P3udih7YYB4wGIrAMbR2DZdTb5/Y+l7pshyq2n0jjGdVT078TE4CnSsu2x9LfXwTMrNbxaukknsdyt16I5UrWt9y6pcC5G8p3BnYABViuWHNp4EM576ItzmHp16vwsfd51HVblH6v1A43tldD/dSwLcYCv2B5v8UIHMXyYFdX1/W2ZztUsW2jCfx38HdiAnAMy82RActN4jPVPZbKzqkoitLEqFw9iqIoTYwK/IqiKE2MCvyKoihNjAr8iqIoTYwK/IqiKE2MCvyKoihNjAr8it0JIeaVm0jCLITIFELsF0L8RQgRWK7c9cknYqq5X6fSfXezcX3HCCEmVrF8mxDiS1seS1FqQ20laVOUmsrGkmoYLG9z34tlroYpQojBUsqDWF7u6wWcqOY+nbC8DHgOiLNhXcdgyai49Iblz2OZTENR6jUV+JX6okRWzLq4SQjxIZY3utcIIdpLKQsBu2VmFELoZMV5IyqQUh6ry/ooyp1SXT1KvSWlzMKSZjYMeLiqrh4hxPDS+UlzS7uI9gkh+pWuvp6zfUm5rqSQ0u18hRDLhBDppXP6bhNC9Ch/fCHEOSHEQiHEHCFECmAQQizFMj1ov3L7nFdavlJXjxDiwdI6FZTOmbtYCOFWbn3/0n30F0J8IYQwls6p+7ztWlJRKlJ3/Ep99xOWuXajuKGLRwgRBnyJJaHVHwAXLDMZ+ZQWeRD4Ecvcvd+VLruevfC/WOayfQlLWtw/AD8JIe6RUp4ud5insOTGeR7Lv5d4oDXgVboMIKWqigshOgIbgc1YLhatgDexJFq7cVLsj4FlWCaaGQv8UwhxQDae+QeUekQFfqVek1IWCiHSgIAqVt8D5Egp/1Bu2YZyf95f+t8z5buRhBCDsUzi3l9Kub102Y9YngX8AXj2huPESCkLym2fAWjk7ScEeQ1IBoZLKU3ltl0jhOglpdxTruxqKeWfS8tsA4ZhyVSrAr9ic6qrR2kIxE2WHwE8S7tsBgkh9NXc333AtetBH0BKmQt8C/S+oezW8kG/hu4D1l8P+qXWYvkFc+NxrPOtSimLsUyv2PIOj6sot6QCv1KvlU7Q04wqZiKSUiYCj2HpOtkApAkhVgkh/G6z2+ZV7a90mU8Vy+5UpeOUXgTSqzhO1g3fi7B0XSmKzanAr9R3A7B0Se6paqWU8jspZR8sF4dfAwOBf9xmn5com96zvAAg48ZD1Ki2tzlO6ZSJzao4jqLUGRX4lXpLCOEFvAWcBrbcqqyUMltKuQpYj2UWM7DcNUPlO+d9gH/p1I7Xj+UKDAV2VaNq1b0b3weMKA3218ViuZBV5ziKUivUw12lvnAQQkSV/tkdy+ic5wBXYLCU0iRExa5+IcSzWF7o2ohljuZ2WKYoXA6WuUyFEEnAGCFEApYZ3g5LKTcJIXZjecj6JyxdLy8BOixTPd7OCeAxIcTjWEb0XJRSXqyi3J+BQ8B/S99JaInlQrbphge7ilKn1B2/Ul94YunO+RnLxNmjgBVA59K3dqtyGPAD/o7l4eirWIZFvlyuzFQsb9luwTLKJ6h0+QgswyzfKz2eAB68YSjnzSwuPd7/le5zSlWFpJRHgSFYunvWYbkQrC49N0WxGzX1oqIoShOj7vgVRVGaGBX4FUVRmhgV+BVFUZoYFfgVRVGaGBX4FUVRmhgV+BVFUZoYFfgVRVGaGBX4FUVRmpj/B2VaoQvYBOBnAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAD8CAYAAACW/ATfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAABUDElEQVR4nO3dd3hUVf748ffJpE1LIyEQShIIIRA6kV5FERRBuijVXbusurqyigKCrvtd5Wf58rUsuhRpCqisCiiiIiAiHUNvSSAQkpA6SSZl5vz+mDAkJCEZmGRSzut55nky95x772cO4ZM75557jpBSoiiKotRvbq4OQFEURal+KtkriqI0ACrZK4qiNAAq2SuKojQAKtkriqI0AO41daLAwEAZFhZWU6dTFEWpF/bt25cqpQy61ePUWLIPCwtj7969NXU6RVGUekEIEe+M49RYsm8I8s+cIfGvz+Hm7Y2+f38MA/rj3aEDwk31limK4loq2TuJtaCAxL8+R/6JEwDkHTpE6qJFuAcHEzzrBYzDhyOEcHGUiqI0VOqS00lS3n3Xnui9O3VCeHoCUHT5Mol/fY7zjz5KwYVEV4aoKEoDppL9TZBSIgsK7O9zdv9O2n+WAGC4/XbCPltD5O7faL7of/EMD7fV+WU7Z++9l0tz55Hz226kxeKS2BVFaZhETc2NExMTI+vDDVprQQFxE+8n/+RJtF27YOg/gPQ1ayi6dAlNo0a0+u8G3Bs1KlX/ykf/5sq//40sLLRv1wQF0uhPf6LR9Oku+BSKotQVQoh9UsqYWz6OSvaOydm1i4QZD5Vb1vzDDzAOGlRuWf7Zs6QtWUr2999jycy0b2/9wxY8mzevjlAVRakHnJXsVTeOg3J27gTATadD17MneHgA4P/ggxUmegCvVq1oumA+bXZsp/kH74NGA0DmVxuqPWZFURQ1GsdBph22ZK/v25fm//seFlMOhYkX8IqMrNL+wsMD4+DBGAYMwPTTT2R+9RWBTzyuhmcqilKtak2GGTRoEF5eXhgMBgwGAxEREbzzzjv28vfee4+ePXui0+mIiIgo9xhvvvkmzZo1Q6/Xc8cdd3D27NlS5Zs3byY6OhqtVkuHDh34/vvvHYqxKCWF/OPHAVuyB9AY9Hi3bevwsErf++4DoPDCBfL27XNoX0WpDcLCwvD09CQ1NbXU9i5duiCEIC4uzjWBlbB161aioqLQ6XQMHjyY+PiKn09KS0tj9OjR6PV6QkNDWbVqVanyjz/+mIiICAwGA8OGDePixYv2snfeeYdWrVrh4+NDSEgIzz77LEVFRfbyuLg4Bg8ejE6nIyoqih9++MH5H7YyUsoaeXXv3l3eyMCBA+WCBQvs73ft2iV1Op387rvvpJRSrl27Vq5bt06+9tprsnXr1mX2X7FihQwKCpL79u2TOTk5cubMmTI6OloWFRVJKaU8c+aM1Gq18tNPP5X5+flyxYoVUqfTyXPnzt0wrpIyNmyQR9tGyaNto2T++fNV3q88lvx8eaJHT3m0bZRMfOmlWzqWorhCaGiojIyMlO+995592+HDh2VkZKQEHPq/VR1SUlKkj4+P/Pzzz2VeXp58/vnnZc+ePSusf//998sJEybI7OxsuX37dunj4yNjY2OllFL+/PPPMigoSMbGxsr8/Hz52GOPyQEDBtj3PX36tExPT5dSSnnlyhU5ePBguXDhQnt5r1695LPPPitzc3PlunXrpK+vr0xOTq7S5wD2Sifk4Fqb7KWUMiYmplSDSSnlkiVLyk32AwYMkC+//LL9fXZ2ttRqtfLnn3+WUko5Z84c2a9fv1L79OvXT86bN++GcZWU+MIL8mjbKHlq6NAq73Mjl16dL4+2jZLHu3WXlpwcpxxTUWpKaGioXLBggYyJibFve+655+Rrr71WKtmbzWb53HPPyRYtWsjGjRvLRx99VObm5koppUxLS5P33HOPDAwMlH5+fvKee+6R50tcSA0cOFC+/PLLsk+fPtJgMMg777xTpqSkVCm+jz76SPbu3dv+3mQySW9vb3ns2LEydU0mk/Tw8JAnTpywb5s8ebKcNWuW/XM98cQT9rLExEQJyNOnT5c5VmpqqhwyZIh8/PHHpZRSnjhxQnp6esqsrCx7nX79+skPPvigSp/DWcm+1nTjlCSlZOfOnRw/fpzevXtXaZ9Dhw7RvXt3+3uDwUCbNm04dOhQueUA3bp1s5dXJSbTzl9txy7uwrlVvqPvA8Cak0O2K77WKcot6tWrF1lZWRw7dgyLxcJnn33G5MmTS9WZNWsWJ0+e5ODBg5w+fZrExETmz58PgNVqZcaMGcTHx5OQkIBWq+Wpp54qtf+qVatYsmQJycnJFBQU8NZbb9nLOnXqVKa75aojR47QuXNn+3u9Xk/r1q05cuRImbonT55Eo9EQWeLeW+fOne11rybMq67+HBsbWypOHx8fAgMDOXToEI8++qg9jlatWmE0Gss9dk2pVcn+9ddfx8/PD71eT79+/XjwwQfp0aNHlfbNzs7G19e31DY/Pz+ysrKqVF6Z/JMnsRT3TeqdlOy9O3TAM6I1AJlffeWUYypKTZsyZQrLly9ny5YtREVF0axZM3uZlJLFixfz9ttvExAQgNFo5KWXXmLNmjUANGrUiLFjx6LT6TAajcyePZtt27aVOv6MGTOIjIxEq9UyYcIEDh48aC87fPgwDzzwQLlxmUymMv/nfX19yc7Odrju3Xffzeeff87hw4fJy8tj/vz5CCHIzc2113/ggQfIysri5MmTPPbYYwQHBzscR3WqVcl+9uzZZGRkkJuby/nz5zl69CgPPVT+mPbrGY1GMkuMXwfIyMjAx8enSuWVySkehYNGYxty6QRCCPyKb9Tm7PqNjC+/wmLKccqxFaWmTJkyhVWrVrF06VKmTp1aqiwlJYXc3Fy6d++On58ffn5+DBs2jJSUFAByc3N59NFHCQ0NxcfHhwEDBpCRkYGlxBPmTZo0sf+s0+kwmUxVistgMJS5mMvKyip1hV3VukOGDOHVV19l7NixhIaGEhYWhtFopHk5z8i0adOG6OhonnjiCYfjqE61KtmX1Lx5cyZMmMAXX3xRpfqdO3dm//799vcmk4lTp07Zv8ZdXw5w4MCBUl/zbuTq+Hptly5oDIYq7VMVPveOtI25l5JLL77Iqb59ufD0MxQmJTntHIpSnUJDQwkPD2fjxo2MGTOmVFlgYCBarZYjR46QkZFBRkYGmZmZ9oS9cOFCTpw4we7du8nKyuKXX34BKNVlcrOio6NLddPm5ORw5swZoqOjy9SNjIykqKiIU6dO2bcdOnSoVN0nn3ySU6dOkZyczNixYykqKqJDhw7lnruoqIgzZ87Y4zh79mypK/nrj10Tam2yT0pKYu3atfZkXFRUhNlsprCwECklZrMZs9lsr//II4/w0UcfceDAAfLy8nj55ZcJDw+nX79+AEydOpW9e/eyevVqCgsLWb16Nfv27WPatGmVxmI1m8ktfvpX37ePUz+nR3Bjmi5YYJ9DR+bnk/3ddyQ+86xTz6Mo1emTTz7hxx9/RK/Xl9ru5ubGww8/zLPPPktycjIAiYmJfPfdd4Cte1Wr1eLn50daWhqvvvqq02IaPXo0sbGxrF+/HrPZzPz58+nUqRNRUVFl6ur1esaMGcOcOXPIyclh586dbNiwgSlTpgBgNpuJjY1FSklCQgKPPPIITz/9NP7+/oBtWObVz3f06FHeeOMNhgwZAtj+kHTp0oVXX30Vs9nMl19+yeHDhxk7dqzTPmuVOOMub1VeVRmN4+npKfV6vdTr9bJx48Zy/PjxMj4+Xkop5dy5cyVQ5lXS//zP/8imTZtKrVYrb7/99jJ3yjdt2iTbt28vvb29Zfv27e3DOm/EarXKK8uW24dc5h48WOk+N8Nqtcq8Y8dk4guz7OfKO3q0Ws6lKM4QGhoqt2zZUmZ7YWFhqdE4eXl58sUXX5Th4eHSaDTKqKgo+e6770opbaNaBg4cKPV6vWzTpo388MMPJSALCwullLa8sHjxYvuxlyxZIvv27Wt/3759e7lixYoKY9yyZYts27at9Pb2lgMHDiw1HPT111+Xw4YNs7+/cuWKHDVqlNTpdLJFixZy5cqV9rL09HTZsWNHqdPpZHBwsPz73/9uH9YtpZTTp0+XjRs3ljqdToaGhsrnn39e5uXl2cvPnTsnBw4cKL29vWVkZGS57VYRnDQaR82NcwN5f/zB5X+8Qd6BA4Bt8rI2P/+MKJ7qoDpYTDmcHjAAa24ufhMn0vTVedV2LkVRaj81N041S1u1irjxE+yJ3rN1a1osWlStiR5sT+T63HsvAFlff61u2CqK4hQq2ZdDWq2kLvo/ADS+vgS//DKtvvoSbRVv5t4q/0n3A2DNzSXrm69r5JyKotRvKtmXI//ECSxpaQA0feMfBEx+EFE8u2VN8I6Ksv9hSV/zmVNGJiiK0rCpZF+OnF9tT8ri7o6uh3PG1DvK737b1X3+8ePklXiIRFEU5WaoZF+OnOJpEbSdO6Mx6CupXT18hg/Drfipu/RVq9XVvaIot0Ql++tY8/PJLZ5yWN+navPyVAc3b2/707VZX3/NuTFjSV+9GksNP2KtKEr9oJL9dfL27UPm5wOg7+PcB6gcFTB9Gu5NmwKQf+wYSa/O58ydQ8kvfjJPURSlqlSyv87V/no3oxFtx44ujcWjaVNab/yWpv/4B9ouXQCwZGRw5eNPXBqXoih1j0r21zEVJ3tdzx4Id9ev2uim1eI3ZjRha1bjXzy7X/Z332HNy3NxZIqi1CUq2ZdQlJZG/tFjgOu7cMpjn/8+N5fsH7a6NhhFUeoUlexLyNm1y/6zvoqLptQk7w4d8GzVCoDM//7XxdEoilKXqGRfwtX+eveQpniGhbk2mHIIIfAdORKwTblcVDwnuKLUVTt27KBPnz74+voSEBBA37592bNnDwBLly5Fo9FgMBgwGAyEh4czY8YMTp486dA5nLXo+G+//cadd95JQEAAQUFBjB8/nkuXLtnLhw8fbo/VYDDg6elJxxL3/Vy96HiDT/ZSSgqTkjBt307O9h2ArQtHCOHiyMrnO9I2bw5WK5nffOvaYBTlFmRlZTFixAhmzpxJWloaiYmJzJ07Fy8vL3ud3r17YzKZyMzM5IcffkCr1dK9e/dSywHeSGpqKmPGjGHBggWkpaURExPDxIkTK6z/5JNP4unpyeXLl1m5ciWPP/64ffnA9PR0HnnkEeLi4oiPj8doNDJjxgz7vps2bcJkMtlfffr0Yfz48fbySZMm0bVrV65cucLrr7/OuHHj7Iu41AhnTJ1ZlVdlUxy7Qu7hw/Jkv/72KYWvvjK//dbVod1Q3JSp8mjbKHnmvtGuDkVRbtqePXukr69vheXXT2d81T333CPHjh1bpXM4c9Hx6+3bt08aDIZyy86dOyfd3Nzk2bNnpZS3tug49XnB8ZqS+n/vl+oKER4e6Pv1w3D77S6MqnK+o0YBtrH35hOOfaVVlNoiMjISjUbDtGnT2LRpE+np6VXab8yYMWzfvt3+vqYWHb/eL7/8UuFqU8uXL6d///6EFy9KVBsWHW+wyb4wMRFT8cLGfhMn0mrjt7Q9sJ+WHy/GzdvbxdHdmPGuoYjir7qZGza4OBpFuTk+Pj7s2LEDIQQPP/wwQUFBjBw5ksuXL99wv5CQENKKJyqEmlt0vKTDhw8zf/583nzzzXLPu3z5cqZPn35Tx64uDTbZp69bB1KCRkPgE0/g1apVrRhXXxUagwHjHXcAkPnFF1hLrHCvKHVJu3btWLp0KRcuXCA2NpaLFy/yzDPP3HCfxMREAgICqnR8Zy46ftXp06cZPnw47777Lv379y9znB07dpCUlMS4ceNuKo7q0iCTvSwsJGPdOgCMQ4bgEdzYxRE5LmDKZMD2RG3GF1+6OBpFuXVRUVFMnz690puvX375ZblJtjzOXnQ8Pj6eO+64g1deecW+Pu31li1bxpgxYzAYDKXicPmi487o+K/KqzbdoM3ctNl+M9a0c6erw7lp5x58UB5tGyVP3T5EWovX7FSUuuLYsWPyrbfekufPn5dSSpmQkCD79Okj//znP0spS9+gLSoqkmfPnpVPPfWU1Ov18vDhw1U6R3JysvTx8ZHr1q2TeXl58oUXXpA9e/assP7EiRPl/fffL00mk9yxY4f08fGRsbGxUkopL1y4IFu1aiX/9a9/Vbh/bm6u9PX1lVu3bi1T1rNnT/ncc8/JvLw8+cUXX0hfX1+ZnJxc6WfASTdoG2Syj5s+3ZYkhw6VVovF1eHctKwff7T/0cr45htXh6MoDrlw4YIcP368DAkJkTqdToaEhMhHHnlEZmZmSiltyd7NzU3q9Xqp0+lky5Yt5dSpU+XRo0dLHaemFh2fN2+eBKRery/1KmnVqlWyZcuW0mq1lonjZhcdd1ayb3ALjuefO8fZ4XcD0Phvf6PRnx5ycUQ3T1qtnB05koLTZ/Bq347w9etr7fMBiqLcHLXg+E3KWPMZAMLTE98xo10cza0Rbm40mmH7Y5V/9Bi5JaZ7UBRFKalBJfus778n7dNPATAOuwt3f38XR3TrfO4dgXtj2w3mKx9/TE19U1MUpW5pMMk+57ffuPjc82C1ogkMJOgvT7s6JKdw8/QkYNpUAHJ+3cXlN95AWq0ujkpRlNqmQST7vD9iufDEk8jCQtyMRlp+vBjP5s1cHZbT+E+ahFf7dgCkL/+Uiy/MQhYUuDgqRVFqk3qf7AuTkzn/6KNYc3MRXl60+PADvKOiXB2WU7npdIQuW4auRw8Asr75hvNPPInFlOPiyBRFqS3qfbK//MYbWNLSwM2NZu++g657d1eHVC00RiMtFv8b4522J2tzduwgftL9FNxgOldFURqOep3sTdu2kb1pMwAB06ZhHDTItQFVMzcvL5q9/TZ+99umcM0/dZpz4ydgKjFplKIoDVO9TfbW3FySXp0P2BYjCXrqSRdHVDOEuztN582jybx54OGBNSuL8488SvI772DNz3d1eIqiuEi9Tfap779P4cWLADR55RXc9HoXR1Sz/O+fSOiyZWiCAkFKrnz4EedG3UfOb7tdHZqiKC5QL5N9zq5dXFmyFADj0KEYBw92bUAuouvWlfB16zEMHAhAQVwcCdOnc+mVOciiIhdHpyhKTapXyd6SkcHFl18mYcZDYLHgptcTPPslV4flUh7BjWn+4Qc0e+cd3IOCAMhYu5bkCubhVhSlfqo3yT77xx85c88IMtetB0DTqBHN3nkbj+BgF0fmekIIfIbdRauN36Lv0weAtGXLSf/8cxdHpihKTakXyT7jq6+48NRMLFeuAOA3fhytv/0GQxXnvG4oNEYjzd59B89WrQBImr+AnN2/uzgqRVFqQp1P9ulr13LpxZfs0yCEfrqcpgsWoPHzc3VotZLGaKTFB++j8fWFoiIS//IXzCdOuDosRVGqWZ1O9mkrV5L0yhyQEvfgYEI/XY7utttcHVat5xkaSrN33wV3dyyZmcRNmEj655+rSdQUpR6rk8neajZzac5cLi94DQCPkBBCV3yKV/FK7krl9L160uzNfyG0WmR+Pklz5nLxueexXLdOpqIo9UOdS/b5Z84QN2EiGcU3Fz1DQwn9dDmeLVq4OLK6x2f4cMLXrcUrMhKArI0bOdV/ABdm/oWszZux5uW5OEJFUZylTqxUJa1WcvfuJWvjRjI3/BdZnIR87r6bJvNfRVNiYV/FcVazmcv/eMP+B9TOwwPv9u3Qde2GtmtXtF274NG47i3Orih1mbNWqnJJspcWC0WpqRSlpmLNNmE1ZWPJNmE1lfOzyUT+8RMUJSdfC9rLi+CXZ+M3bpxahs+JzCdPkrVxI1kbN1GYkFBuHY/mzdF27YpX61Z4NG2Ke5OmeDRtgnuTJrh5etZwxIpS/9W5ZN8tKkp+0T2GokuXKExOhpt4gtPNxwfj0DtpNH06XhER1RClArZF6M1HjpK7eze5B/aTt/+AbebQSmgCA3EPCMBqNmPNzsZqMpVaSEW4u+NmMKAxGGzTV7jVuV5ERalRvvfdR6Mpk52S7N2dEVBVSKuVvH37blxJo7ElAqMRN6MRjV6Pm9GIe2AghiG3Y+jTB6GuHqudEAJth2i0HaJpxENIKSlMSCD3wAHy9h8g7/BhChMTsWZnl9rPkpqKJTW1wuNKiwVLfr79eQhFUW5M37u3045VY8nezcsLnzvvLP7aX/zVPygIN6MPGqMBN4MB4e2tumVqISEEnqGheIaG4nffffbtFpPJ9k0tKYnCS5coSkqi6EoabjodbkbbFTzu137FZEEBVlOOrYsuJwdQQz0V5Ua8O3Zw2rHqxA1aRVGUhqpG++yFEG7AXcAdQA+gCeANpAEngZ3AF1LK8u/q2Y6RDahHNW0CgYr7OxoW1RbXqLa4RrXFNW2llMZbPcgNk70Qwgg8CzwGNAKOAIew/SPkA35AGBCD7R/nR2CBlLLM0khCiL3O+OtUH6i2uEa1xTWqLa5RbXGNs9qisj77c8AB4G/Af6WU2RVVFEJ0BSYAXwghXpVSLrrV4BRFURTnqCzZ3yWlrGQIjY2U8gBwQAixAGh5y5EpiqIoTnPDZF/VRH/dPrnA8XKK/u3oseox1RbXqLa4RrXFNaotrnFKW9TYaBxFURTFdZzyCKMQ4gchxFZnHEtRFEVxPmc9ry6ceCxFURTFyVQ3jqIoSgOgrsYVRVEagCrPjSOEuLuyOlLKjbcWjqIoilIdqtyNI4SwYpu56vqZyuwHkFJqnBeaoiiK4iyOzHpZ3gKvAcBQYDowwxkBKYqiKM7nlBu0QojngQFSypEV1QkMDJRhYWG3fC5FUZSGZN++falSyqBbPY6z5rM/AMy7UYWwsDDq8hTHuVkFeOnc0bire9qKotQcIUS8M45zy8leCOGJrRvn0i1HUwulJGTz+9dnifvjCkEtjYx+rhseXurWhKIodYsjo3H2UHZpIU9sUxwbqWd99tlpZnauO8WZ/Sn2bSkJ2fy04jh3PtRerailKEqd4siV/RHKJnszsBb4Skp5xGlRuZi0SjZ+cJjU8yYAPL01BIQYSDqbyak9l2na2peOg5q7OEpFUZSqq3Kyl1JOr8Y4apVzh1Ptib7jwGb0GNkKdw831v1rH1cumNix9hRBLY00aeXr4kgVRVGqRt1tvI6Ukn2b4gDwCfSm34Q2eOs9cPfUMPzRDnhq3bFaJJs++oPYbRfIzy10bcCKoihV4KxZL18WQrzijGO52vmjaSTH2xbk6nZXKG6aa03kG6TjjuntAMjNLGDb6pMsmbWTH5YcJTerwCXxKoqiVIWzruznUcnQy7pib/FVvcHfi6heTcuUh3cOYvhjHQkO9wHAUmjlxO4kvv8kFjWpnKIotZVTxtlLKZ01Xt+lLp5K59LpTAC6Dm2JxqP8v4WtugTRqksQVy6a2P9dPCd3XybxRAbxsVcI6xhYkyEriqJUSa3psx80aBBeXl4YDAYMBgMRERG888479vL33nuPnj17otPpiIiIKLP/1q1bGTJkCAEBAQghuHDhQpk6mzdvJjo6Gq1WS4cOHfj+++9Lle/dGAeA1uhB+74hlcbcKMTA7ZPbYWzkDcCv609jtVgd+NSKUneFhYXh6elJampqqe1dunRBCEFcXJxrAith69atREVFodPpGDx4MPHxFT+flJaWxujRo9Hr9YSGhrJq1apS5R9//DEREREYDAaGDRvGxYsX7WVvvvkmHTp0wGg0Eh4ezptvvllq34MHD9K/f398fX1p3rw58+fPd+4HrYKbSvZCCDchhO76160G88orr2AymTCZTKxYsYLZs2fbE3JISAgvvPACs2fPLndfvV7P1KlT+fTTT8stP3v2LGPGjOHFF18kMzOTF198kdGjR9t/IdMu5XD+WDoAXe5sibtn1R6c0ni40fu+1gCkJ+VydGe9fLZMUcoVHh7O6tWr7e//+OMP8vLyXBjRNampqYwZM4YFCxaQlpZGTEwMEydOrLD+k08+iaenJ5cvX2blypU8/vjjHDliG1G+bds2XnrpJTZs2EBaWhrh4eFMmjTJvq+UkuXLl5Oens7mzZtZtGgRa9assZc/8MADDBgwgLS0NLZt28YHH3zAf//73+r78OWocrIXNrOEEKeBQiC7nJfT9OrVi/bt2xMbGwvAuHHjGDt2LM2aNauw/rRp04iOji63fNmyZXTv3p3Jkyfj6enJgw8+SLdu3Vi2bBkAZw/YHp5ycxNVuqovKSKmMY1DjQD8/s05CsxFDu2vKHXVlClTWL58uf39smXLmDp1aqk6+fn5PP/887Rs2ZLg4GAee+wx+x+E9PR0RowYQVBQEP7+/owYMaLUt/JBgwbxyiuv0LdvX4xGI0OHDi3zTaIiX3zxBdHR0YwfPx5vb2/mzZvHoUOHOH78eJm6OTk5rF+/ngULFmAwGOjXrx8jR460Xzx+/fXXjB8/nujoaDw9PXnllVf45ZdfOHPmDAAvvPAC3bp1w93dnbZt2zJq1Ch27txpP35cXBwPPvggGo2G1q1b069fP/sfkpriyJX9X4C/A59gm+b4dWA+cBKIAx5xVlBSSnbu3Mnx48fp3bu3U4556NAhunfvXmpbt27dOHToEABnD9qSfbO2fnjrPRw6thCCvuNsXUt5WQXs2+yUqSwUpdbr1asXWVlZHDt2DIvFwmeffcbkyZNL1Zk1axYnT57k4MGDnD59msTERHs3htVqZcaMGcTHx5OQkIBWq+Wpp54qtf+qVatYsmQJycnJFBQU8NZbb9nLOnXqVKa75aojR47QuXNn+3u9Xk/r1q3LTbInT55Eo9EQGRlp39a5c2d7XSllqQEYV3++ejFakpSS7du3l7rwfOaZZ1i+fDmFhYWcOHGCXbt2cccdd5Qbd3VxJNk/DMwF/lX8/isp5atANHAcaHOrwbz++uv4+fmh1+vp168fDz74ID169LjVwwKQnZ2Nr2/ph6D8/PzIysoiO81MSoLti0mrLjc3uVxIG3/CO9tuzu7fHM/3nxwhL/vacMzCAgtpl3LUiB2l3rl6db9lyxaioqJKffuWUrJ48WLefvttAgICMBqNvPTSS/YujkaNGjF27Fh0Oh1Go5HZs2ezbdu2UsefMWMGkZGRaLVaJkyYwMGDB+1lhw8f5oEHHig3LpPJVOb/vK+vL9nZZTshKqt799138/nnn3P48GHy8vKYP38+Qghyc3PLHGvevHn2P2JXjRgxgnXr1qHVaomKiuJPf/oTt912W7lxVxdH57M/KKW0CCEKAT8AKaVVCPE+8DG2K/+bNnv2bF5++WUALly4wAMPPMBDDz1k72q5FUajkczMzFLbMjIy8PHxsV/Vg21o5c3qN74NqedNZKeZObXnMuePpRHVqwnJ8dkkncvEWiRp37cpg6e0u+lzKEptM2XKFAYMGMC5c+fKdOGkpKSQm5tb6lu1lBKLxQJAbm4uzz77LJs3byY93XbPLDs7G4vFgkZju2/WpEkT+746nQ6TyVSluAwGA1lZWaW2ZWVlYTQaHa47ZMgQXn31VcaOHUtmZibPPvssRqOR5s1LT5uyaNEili9fzvbt2/Hy8gJsN36HDRvGokWLeOCBB0hKSmLcuHEEBwfzxBNPVOmzOIMjV/ZXAEPxzwlA1xJl/oDWWUEBNG/enAkTJvDFF1845XidO3dm//79pbYdOHCAzp07c6442QeH+6D387rpc/gEarl/Tg86Dm4OAsymQg7+cJ6LpzKwFtmu6I/uvMSJ39RNXKX+CA0NJTw8nI0bNzJmzJhSZYGBgWi1Wo4cOUJGRgYZGRlkZmbaE/bChQs5ceIEu3fvJisri19++QXAKd+Ao6Oj7d20YOuXP3PmTLn39SIjIykqKuLUqVP2bYcOHSpV98knn+TUqVMkJyczduxYioqK6NChg738P//5D//85z/ZunVrqT8CZ8+eRaPRMHXqVNzd3WnevDn3338/GzfW7CqujiT7ncDV7x2rgHlCiNeFEHOB/wdsdWZgSUlJrF271t7nVlRUhNlsprCwECklZrMZs9lsr2+1WjGbzeTn5wO2m0Jmsxmr1TYUcurUqezdu5fVq1dTWFjI6tWr2bdvHxPHTeLiqQzg5rtwSvL0dmfAxEjG/q07QS2NeHhpaNHOn173tcIv2DZg6efVJ8m4XPbrn6LUVZ988gk//vgjer2+1HY3Nzcefvhhnn32WZKTkwFITEzku+++A2xX8VqtFj8/P9LS0nj11VedFtPo0aOJjY1l/fr1mM1m5s+fT6dOnYiKiipTV6/XM2bMGObMmUNOTg47d+5kw4YNTJkyBQCz2UxsrO3ByYSEBB555BGefvpp/P39AVi5ciUvvfQSW7ZsoVWrVqWOHRkZiZSSVatWYbVaSUpK4rPPPit1P6FGXL3xUNkLaAvcXvyzF/AukAikAZ8BjW+0f/fu3eWNDBw4UHp6ekq9Xi/1er1s3LixHD9+vIyPj5dSSjl37lyJbdbNUq+rfvrpp3LLf/rpJ3udTZs2yfbt20tvb2/Zvn17+d1338mjOy/KRY9ulYse3SrTk3JuGOPNsFqt9p9TzmfLD578SS56dKtc89puWVhQ5PTzKUpNCQ0NlVu2bCmzvbCwUALy3LlzUkop8/Ly5IsvvijDw8Ol0WiUUVFR8t1335VSSpmYmCgHDhwo9Xq9bNOmjfzwww8lIAsLC6WUtrywePFi+7GXLFki+/bta3/fvn17uWLFigpj3LJli2zbtq309vaWAwcOtMckpZSvv/66HDZsmP39lStX5KhRo6ROp5MtWrSQK1eutJelp6fLjh07Sp1OJ4ODg+Xf//53WVR07f9vWFiYdHd3t+cvvV4vH330UXv51q1bZUxMjPTx8ZHBwcHyz3/+s8zJqVq+AfbKKubpG72csixhVcTExMjauFLVt+8fJu5wKgEheibN6Vnt54vddoFtq08Cthk1B0xqW+3nVBSl7hJC7JNSxtzqcWrNE7SuUJhv4fyxNAD7SJrqFj2gGa272rqL/tiWyHHVf68oSg24YbIXQmwVQlR5MKgQorEQYr4QYuath1b94v5IxVJo69N3Rn99VQghGDwlCt/GtvvZP684QdK5zEr2UhRFuTWVXdlvBFYKIeKFEG8LIcYLIdoJIYKEED5CiJZCiIFCiL8KIbYAF7CN0vmu2iO/RZZCK79/fQ6wzVsf1LLscKzq4qXz4O7HO+HhrcFSZGXTh3+Qk5FfY+dXFKXhuWGyl1IuBFoBrwExwEogFkgC0oFzwE/Ac9gerOohpbxXSnmyOoN2hv3fx9tHxPS6r3WNrykb0FTP0IeiQdjmxt/44R9qmgVFUapNpX32UsocKeViKWV/bA9S9QPGAZOBu4FWUspmUsqZUsqD1Rmss2RczmXfJtuUBi3bBxDRvbFL4gjrFEivUbZhWslxWaz7n32kJ+W4JBZFUeo3h+ahl1LmAruqKZYaIaVk2+oTWIqsaDzcGDCpbY1f1ZfU7a5QTGn5xP6SSPqlHNa+sZfbp7Zz2R8gRVHqpwY3GufUnstcOG57LDvm7jB8g5z64K/DhBAMfKAtg6dEoXF3ozDfwneLY+0rZimKojhDg0r2lkIru760TUnq31RP1ztbujiia9r3DWHsC93tC6Hs3nCWE7uTXByVoij1RYNK9sd/u4Qp3Tbqpd+4CDTutevjB7U0Mm5WDD6BtoT/4/JjXDyV7uKoFEWpD2pXtqtGVouV/d/Zbso2DjXSon2AiyMqn87Hk3ue7IyXzh2rRbLxgz/UTVtFUW5Zg0n2p/ZcJivVNnFazN1hLr0pW5mApnqGPdoRNzdBfm4R3/zfYcw5ha4OS1GUOsyRZQn7CyFGlXgfKIRYJYQ4KIRYKIRwbHmnGiSt0r56VKNmBsI61szUCLeieVt/Bk22zc6XlZLHd4tj1WLmSr2zY8cO+vTpg6+vLwEBAfTt25c9e/YAsHTpUjQaDQaDAYPBQHh4ODNmzODkScce43HWouNHjx4lJiYGf39//P39ueOOOzh69Ki9fN68eXh4eNjjNRgMnD17FoCEhIRS2w0GA0IIFi5c6NBnuRWOXNn/C+hQ4v27wBDgN2A64Ly5SZ3szIEU0pNsD1B1Hx6KcKu9V/UltevT1H4T+cLxdH5df8bFESmK82RlZTFixAhmzpxJWloaiYmJzJ07177oB0Dv3r0xmUxkZmbyww8/oNVq6d69e7nLAZbHmYuOh4SEsG7dOtLS0khNTWXkyJHcf//9pfafOHEiJpPJ/ro63XHLli1Lbf/jjz9wc3Nj7NixjjbbTXMk2bcF9gEIIXTAaOBpKeVjwAtAxS3oQlJK+zBGv2AdrbvVrfHrvUa3pmW07f7CoR/Pc+xXNXGaUj9cvUKfNGkSGo0GrVbL0KFD6dSpU5m6Vxfqfv/99xk4cCDz5s2r0jmcuei4n58fYWG2LmApJRqNhtOnT9/UZ1++fDkDBgwgLCzspva/GY4ke0/g6mohfbE9kPVt8fuTQFMnxuU08X9c4coF26o43YeF4lZHruqvcnMTDP1T9LWFT1YdJzk+q5K9FKX2i4yMRKPRMG3aNDZt2mRflrAyY8aMYfv27fb3NbXo+FV+fn54e3szc+ZMXnrppVJlX3/9NQEBAURHR/PBBx9U+BmWL1/OtGnTbvxBncyRZH8cGFb884PALinl1ZV7Q7AtYlKrlLyqNzbypk2PYNcGdJNsE6d1xNNbg7VI8uPyY1iKVP+9Urf5+PiwY8cOhBA8/PDDBAUFMXLkSC5fvnzD/UJCQkhLu5ZuamrR8auuLq24aNEiuna9tjrrhAkTOHbsGCkpKSxevJj58+ezevXqMufZvn07ly9fZty4cTf8nM7mSLKfDzwrhEgBHgD+WaJsGHDAmYE5w4Xj6Vw+Z7sK7nZXKBpN3R185N9ET78JtiuOK4k5HPi+4ptMilJXtGvXjqVLl3LhwgViY2O5ePEizzzzzA33SUxMJCCgakOnnbnoeEl6vZ7HHnuMqVOn2pdbbN++PSEhIWg0Gvr06cPTTz/NunXryuy7bNkyxo4di8FgKFNWnaqc/aSU/wXaAY8BHaSUm0oU78I2M2atsndjHAB6X0/a9a6VvUwOierdhBbtbGte7tkYR9olNf5eqT+ioqKYPn16pTdfv/zyS/r371+lYzp70fGSrFYrubm5JCYmllt+tW+/pLy8PNauXVvjXTjg2NDLqUCmlHJ9OVMYr8N2A7fWuHgqw76QeNehoWg86u5V/VVCCAY9GIW7pxvWIslPnx5HWmtmWUlFcbbjx4+zcOFCLly4AMD58+dZvXo1vXr1KlPXYrFw7tw5Zs6cyc8//8zcuXOrdA5nLjq+ZcsWDhw4gMViISsri7/+9a/4+/vTrl07ADZs2EB6ejpSSn7//Xfee+89Ro0aVeocX375JX5+fgwePNihtnIGRzLgEqB1BWXhxeW1xr7ivnqt0YP2/UNcG4wT+QRq6TXK9s+QdDaT3zacuWH/fX5uIfGxV8jPU3PlK7WL0Whk9+7d9OzZE71eT69evejQoUOpsee7du3CYDDg4+PDoEGDyMrKYs+ePXTs2NFeJzo6mpUrV5Z7jqCgINavX8/s2bPx9/dn9+7drFmzxl7+j3/8g+HDh9vfv//+++Tl5dG4cWMmTZrEBx98YL+yz8jIYNKkSfj6+tK6dWtOnz7N5s2b8fa2TW+yZs0aIiIiMBqNTJ06lVmzZpW5gl+2bBlTp051yUOdVV5wXAhhBXpJKX8vp2wo8JmU0r+i/WtywfHLcVms+6ftXL1Ht6bbXaE1ct6aYrVKvnhzn/1+hF+wjr7jIuwPi1mKrCSdyeTorxc5sz8FS6GVxqFGxs6KqXOjkRSloXPWguM3nM+++InZkt9DXim+QVuSN9Af2HOrwTiDlJJf19vGvnrp3OkwoJmLI3I+NzfBXQ93YOvSoySezCDjci7f/t9htEYPCs0WigrLXuknx2dzfNcl2vetP99yFEWpusq6cRoDHYtfYOvG6XjdKxT4Hni0mmJ0yKm9l+199d2Hh+GpdWh9ljrDGODNqGe7MuyRDhgDbF8j87ILSyV6b4MHne9oQaPmtrv+v204S4HqzlGUBumGmVBKuRhYDCCE+Al4XEpZ9tGzWqLAXMSv62xX9f5N9XS6vbmLI6peQghad2tMaIdGHN91iZzMArx07nhq3TH4edGsrT8adzcunsrgy4X7ycsqYN/mOHqPjnB16Iqi1LAqX/ZKKWv+9rGD9n4bR05mAQADJrap0+PqHeHuqaHDwIr/sIW08aN1t8ac2Z/Mwa3nad8vBN8gXQ1GqCiKqznUxyGECAFGAM2x9dWXJKWUs5wVmKPSLuVwaOt5ACK6N6Z5VO2cr95V+oxpTdzhVCxFVn5df4bhj3WsfCdFUeqNKid7IcRoYDWgAZKBguuqSMAlyd5qlfyy+gRWq8Td042+41Q3xfV8ArV0ubMF+zbFc/ZgCr9tOEPPka1q9bz+iqI4jyNX9v/AdiN2upSyVs2Ds39zPIknMwC47Z5wDP7Xf+lQwDZlRNzhK1xJNLFvUzyWQit9xkaohK8oDYAjndotgPdqW6K/eDqD37+2LRDQNMKXLne0cHFEtZentzv3PduVoJa2uT4O/nCe7WtOqqdwFaUBcCTZ/0otmxLBnFPIlk+OICV46z0Y+qdo3BrITdmb5W3wYNQzXQgO9wHgj22JfPN/h8jJzHdxZIqiVCdHMuNfgUeEENOEECFCCN31r+oKsjxSSrYuO4Yp3Zakhkxrp7pvqshL58HIp7sQ0sYPgIQjaaxZ8DvnDqVw5aKJwz+dZ/O//2DbqhP29lUUpW5zdLqEq8rdSUqpqWh/Z0+XELvtAttW2+Zj63x7C/pNaOO0YzcUFouVPd+cY//meCr6NfD01tBnbATt+4bUmeUcFaU+qZHpEq7zEBUk+ZqWnpTDzuKHpwJbGOg9uqL52ZQb0Wjc6DWqNS3bN+KHJUfJTiteiExAoxADaZdyKDBb+HnlCU7+fpn2/UJoFumPwd/rxgdWFKXWceShqqXVGEeVWYqsbPnPUYoKrbh7uDH0T9H1YvpiVwpp48fEV3oQdygFD293Qtr44a33ICUhm59WHCclIbvUlNF+wTq6Dm2p5tlRlDqkzk0c8/s350hJsC0T1mdsBP5N9C6OqH7w0rrTtlfpBV6CWhoZN6s7h7ZeIHZ7IlkpeQBkXM7lp0+PU1RgodPga6Of4mOvcGZ/Ml56D3wDvTEGagmJ8MPDq8LePUVRaogjD1WlUEk3jpSy8S1HVAGrxcrRnZfY/51tOb7QDo3oMLD+zWhZ27hp3Og6tCVdh7Yk60oeiScy2LcpjsyUPLZ/dgqNuxuRPZqwc90pjmy/WGZ/g78X9zzZicDmZZd2UxSl5jhyg3YeZZN9AHA74AN8IqWcX9H+N3uDVlolp/Ze5vdvzpGZbLuy1Bo9uP+Vnuh8PB0+nnLrstPMfLlwP9lXzCBsCd2UZhu1ozV6oHF3w5SRb/9tcffSMPSh9oR3DnJh1IpSN9X4DVop5bwKAhHA54DT587NzSpg80d/cOlMpn1bUEsjgx5sqxK9CxkDvLnv2a58uXA/pvR8e6KP7BHMgElt8dK6YymycvZACj8uP0ZRvoWNH/5BjxHhtOvTVA2RVRQXqPKV/Q0PIsRdwBIpZYV37G50ZZ+dZuZKoolGzQz2udlTL5j49v1D9kQSEKKn572tCO8SqB7vryUyLufy9f8eJD+viP4TImnbs0mZOsnxWWx8/7B9NlKw3eBtEu6DVUryc4soyC3CWuIpXq3Rk8jbggnvEoi7h+rvVxo2Z13ZOyvZPw4skFIGVlQnJiZG7v7td7KvmMm6kkdWqpmU89kkHk8ns/jGH0DjUCPN2voTuy2RwnwLALeNCCfm7jC1pF4tZLFYEXDDJ5dN6fl8/0ksl05nVlinPF56dyJjgvFrosdL5257ad3x0nngpXNHSsi6kkd2ah45mQXofDzxCfTG2EiLl67OjT1QlDI0GjfcPTU1m+yFEE+Us9kTaAc8CKyVUs6oaP8OUZ3l44PfdmikvsbDjSHT2tEmJrjqOym1VnaamcQT6Vw4kU7qeRMeXm546Tzw1LqjcS/+Qy4h6VwWGZdzXRusotQC3e4Kpc+YiBp/qGpROdvygQvA+8CrN9rZTSNKJ3ph6/tt1saPZlH+BLUwcvFUBmcPppB4MgOdjyfDH+tIcJiPAyEqtZkxwJuo3k2J6t30hvWklFw6k8mxnReJP5KG2VRY6WRt7h5u5a69qyiKjSM3aG/pySV3Tw0DH2iLTyNvfAK1GAO8yzwM1aiZgY6DmlNgLsJNI1R/bQMlhCAkwo+QCD/AlvwL8y3k5xbZ+vjzCjHnFIEEY6Dt98lL606BuYisVDNZqXkUFVpc+yEUxQkCmjrvOaIa69h00wg6DKjauHhPb9XfqlwjhMDT2x1Pb3eMN1iAzNPbncDmBgKLF1hXFOUah27QCiFaAX8D+mEbY58GbAfeklKerWTfbODEzYdarwQCqa4OopZQbXGNaotrVFtc01ZKectPJTpyg7Y78BNgBr4BLgPBwD3Y1qMdLKXcf4P99zrjJkN9oNriGtUW16i2uEa1xTXOagtH+kveAg4Aw6WU9qESxfPYbywuv/1WA1IURVGcz5Gbrj2Af5VM9ADF798CejozMEVRFMV5HEn2eUCjCsoCsHXv3Mi/HThXfafa4hrVFteotrhGtcU1TmkLR/rslwF3AhOklDtKbO8HfAZskVJOd0ZQiqIoinM5kuwbARuA3kAKthu0jYtfvwL3SSmvVFOciqIoyi1weG4cIcQw4DagKXAJ2C2l/L4aYlMURVGcxCkToSmKoii12w1v0AohGgkh1hdPYVxRnbuK61TbKlWKoijKralsNM4zQCvgRt003wPhwHNOiklRFEVxssqS/QTgQ3mDvp7iso+AUc4MTFEURXGeypJ9KHC0Csc5BoTdcjSKoihKtahsuoQ8bIuJV8ZQXLdCgYGBMiwsrIphKYqiKAD79u1LlVIG3epxKkv2+4GRwLeV1BtVXLdCYWFhVLQGbUNSWJBPakIcqQnxBIWG06R1G1eHpChKLSaEiHfGcSpL9v8HfC6E+FVKuayCQKYCM4CJzgiovjp/5DA/L/+ElIRzSKttRSU3jTsjnnmBNj36uDg6RVHquxsmeynlF0KId4ElQoingM1AArYFBlsCdwExwNtSyi+rO9i6Kis1hQ1vvU5+bk6p7VZLEV+//U/unvk8UX0GuCg6RVEagkqnOJZSPieE+BnbMMznAa/ionxgJzBKSvlNdQVY10mrlc3vv01+bg7CzY3+k6bRNDIKNzc3vvrXAvKys9j43ltYLRba9x9cat/M5Msc37mNRi1CiYhRk4oqinLzHF2pyp1rM19ekVIWVXXfmJgY2RD77Pd+8yXbPv0EgD4THqT32En2stTz8axdMJvczAwAglqG0aZnX5pERHLk5x84+dtOpLQi3Nz48/9+gk/gLd+jURSljhFC7KvpxUsoTu6Xb/WkDUVKQhw7VttudTRt05ae900oVR7YIpSJ8/7J2gWzMaVdISUhjpSEuDLHkVYrZ/b+Rtdh99ZE2Iqi1EOOzGevOOBK4nm+fvufWIqK8PDy5u6nnsdNoylTLyCkOdMXvs/Qx/5CeJfuuGlsf3+9dHp6jBpHUFgrAE7v+a1G41cUpX5x6MpeqZyUkj+2fsdPyxZTVJAPwKBpD+PXpGmF+3jp9HQcPJSOg4dizjFx5cJ5glqG4qnV4aU3kBJ3lvNH/yDPlI3WcMvrDiuK0gDVmiv7QYMG4eXlhcFgwGAwEBERwTvvvGMvf++99+jZsyc6nY6IiIgy+2/dupUhQ4YQEBCAEIILFy6UKk9MTGTUqFGEhoYihGDFihVO/wxFhYV8/f/eYMviRRQV5OPu4cmQPz1Bx9uHVvkY3noDzdq2w1OrAyDitl6ArSvn3P49To9ZUZSGodYke4BXXnkFk8mEyWRixYoVzJ49m++/t83BFhISwgsvvMDs2bPL3Vev1zN16lQ+/fTTcsvd3NwYOnQoq1atonnz5tUS/+EfNnPq918BCAoNZ/I/36HL0LsRQtz0MQNCmhMQYotXdeUotUlYWBienp6kpqaW2t6lSxeEEMTFxbkmsBK2bt1KVFQUOp2OwYMHEx9f8fNJixYtIiYmBi8vL6ZPn+7QsYYPH26/UDUYDHh6etKxY8cyx9i2bRtCCF5++WWnfD5H1KpkX1KvXr1o3749sbGxAIwbN46xY8fSrFmzCutPmzaN6OjocsubNm3Kk08+Sd++fdGU03fuDKeLE31wqwgeeG0hjZq3dMpxI3r0BuDcoX0UFncNKUptEB4ezurVq+3v//jjD/LybjhzSo1JTU1lzJgxLFiwgLS0NGJiYpg4seJnP0NCQnj55Zd56KGHHD7Wpk2b7BeqJpOJPn36MH78+FLHKCws5Omnn6ZnT9cMo66VyV5Kyc6dOzl+/Di9e/d2dThVkmfK5sLxIwC07d0fd09Ppx37aldOUX4+8YcPArB/03/54JHJ7N/0tdPOoyiOmjJlCsuXL7e/X7ZsGVOnTi1VJz8/n+eff56WLVsSHBzMY489Zv+DkJ6ezogRIwgKCsLf358RI0aU6oIdNGgQr7zyCn379sVoNDJ06NAy3yQq8sUXXxAdHc348ePx9vZm3rx5HDp0iOPHj5dbf8yYMdx33300atSoTJkjx4qLi2P79u1MmTKl1PaFCxcydOhQoqKiqhS/s9WqZP/666/j5+eHXq+nX79+PPjgg/To0cPVYVXJuQN77dMgtI7p5dRjN2nVBoN/AACnf9/FLyuX8NPSf5ObmcHvG9aiVhtTXKVXr15kZWVx7NgxLBYLn332GZMnTy5VZ9asWZw8eZKDBw9y+vRpEhMTmT9/PgBWq5UZM2YQHx9PQkICWq2Wp556qtT+q1atYsmSJSQnJ1NQUMBbb71lL+vUqROrVq0qN7YjR47QuXNn+3u9Xk/r1q05cuSIw5/TkWMtX76c/v37Ex4ebt8WHx/Pf/7zH+bMmePwuZ2lViX72bNnk5GRQW5uLufPn+fo0aPlfqWqjc4U96f7hzQnIKT8rqabJdzcaH2b7RvOkW0/sOe/6+1lOelpXLmQ4NTzKYojrl7db9myhaioqFJdrVJKFi9ezNtvv01AQABGo5GXXnqJNWvWANCoUSPGjh2LTqfDaDQye/Zstm3bVur4M2bMIDIyEq1Wy4QJEzh48KC97PDhwzzwwAPlxmUymfD19S21zdfXl+zsbIc/oyPHWr58eZk+/7/85S8sWLAAg8Hg8LmdpVYl+5KaN2/OhAkT+OKLL1wdSqWKCgs5d8g26Wd1TWtwtSvnqqYRbRHC9s93tWtHUVxhypQprFq1iqVLl5bpwklJSSE3N5fu3bvj5+eHn58fw4YNIyUlBYDc3FweffRRQkND8fHxYcCAAWRkZGCxWOzHaNKkif1nnU6HyWSqUlwGg4GsrKxS27KysjAaHR++XNVj7dixg6SkJMaNG2ff9vXXX5OdnX3D+wU1odYm+6SkJNauXWv/6lRUVITZbKawsBApJWazGbPZbK9vtVoxm83k59tuYObn52M2m7EWd60A9n2klBQWFmI2mykqqvKMDxU6H3uIQrOtD9LZXThXtWjfAZ2vHwCtut3G+DmvE9zaNgQ1/o8D1XJORamK0NBQwsPD2bhxI2PGjClVFhgYiFar5ciRI2RkZJCRkUFmZqY9YS9cuJATJ06we/dusrKy+OWXXwCc0jUZHR3NoUOH7O9zcnI4c+ZMhYM4nHGsZcuWMWbMmFJX8Fu3bmXv3r00adKEJk2a8Nlnn/HOO+8walTNLu5Xq5L91a85BoOBzp07ExwcbO+Pe+2119BqtTzyyCOcPXsWrVaLVqu17/vLL7+g1WrtNz8iIiLQarX2Xx7Avk9CQgIPPfQQWq2W11577ZbjPrNvt+34Pr40bRN5y8crj8bdgwlz/sGIZ2Yx8rnZeHh5E9qxKwAXjsZiKSqslvMqSlV88skn/Pjjj+j1+lLb3dzcePjhh3n22WdJTk4GbM+8fPfddwBkZ2ej1Wrx8/MjLS2NV1991WkxjR49mtjYWNavX4/ZbGb+/Pl06tSpwhukVy8oLRYLFoul1MVgVY6Vl5fH2rVry3ThLFiwwH7P4uDBg4wcOZKHH36YJUuWOO2zVomUskZe3bt3l/WR1WKRHz46Rb414R65+YN3avTcCUcOy7cm3CPfmnCPTDhyuEbPrSihoaFyy5YtZbYXFhZKQJ47d05KKWVeXp588cUXZXh4uDQajTIqKkq+++67UkopExMT5cCBA6Ver5dt2rSRH374oQRkYWGhlFLKgQMHysWLF9uPvWTJEtm3b1/7+/bt28sVK1ZUGOOWLVtk27Ztpbe3txw4cKA9JimlfP311+WwYcPs7+fOnSuxTd9uf82dO7dKx5JSylWrVsmWLVtKq9V6w3abNm2anD179g3rlATslU7IwQ7Nenkr6uusl0mnT7Jy9l8BGPW3V2p0KmJLUSH/99AkCvPN9Bw9kX73T6l8J0VR6hRnzXpZq7px6orstFSO7fiZ3774jJ+Lpy929/QitGPnSvZ0Lo27B83bdwBUv72iKDemJkJzUH5uLsv/NhOzqfSQq9BOXfHw8q7xeEI7duXcgb1cPnMas8mEtwuHdimKUnupK3sHnTu4157oPbU6gkLDiezdnwEPznBJPKGdugAgpZWEI4duXFlRlAZLXdk76Oy+3wHwCQrmz//78S1NcuYMjZq3xOAfgCk9jfjDB4js2Rew3Xh3dWyKotQeKtk7wGqxcO6A7SZz65getSKZCiFo2bELR3/5kRO/bifpzCmyki8jkYx4ehZhnbu5OkRFUWoB1Y3jgIsnjmHOsT0M0rpb7VkAPKyTbbx9fm4OyefOYM4xkZ+Tw9dvv1HuMoeKojQ8Ktk74Mx+WxeOp1ZH8/aOP4VXXdr06ke7/oNp3q4D0QOH0GPUODQeHhTk5fHlP18lJyPd1SEqiuJiqhvHAWf22p6UDevcDY27h4ujucbdw4O7n3qu1LbG4a355p3/IftKCl/9az4T5r7hktFCiqLUDurKvorSLiaSfikRgNbda/+0y21796ff/bZJqZLOnGLrfz50cUSKoriSSvZVdLZ4/hsh3AjvessPs9WIHveNp33/wQAc3fYjmclJLo5IURRXUcm+iq7214e0bYfW6OPiaKpGCEH/B2fgpnFHSiv7Nm5wdUiKUsqOHTvo06cPvr6+BAQE0LdvX/bs2QPA0qVL0Wg09skRw8PDmTFjBidPnnToHI6sQ5uWlsbo0aPR6/WEhoaWWhjl6NGjxMTE4O/vj7+/P3fccQdHjx61l//0008MHjwYX19fwsLCSh03OTmZSZMmERISgq+vL3379mX37t0OfY5bpZJ9FeSZskk8bvtHrQtdOCUZ/ANo138QALE/bsFcxbnAFaW6ZWVlMWLECGbOnElaWhqJiYnMnTsXLy8ve53evXtjMpnIzMzkhx9+QKvV0r17d/va1JVxdB3aJ598Ek9PTy5fvszKlSt5/PHH7atRhYSEsG7dOtLS0khNTWXkyJHcf//99n31ej0PPfQQb775ZpnjmkwmbrvtNvbt20daWhrTpk3jnnvuqfLc/E7hjNnUqvKqq7NeJsefk+v/Mcc+u2TqhQRXh+SwlIQ4e/y/ffm5q8NRFCmllHv27JG+vr4Vll8/w+VV99xzjxw7dmyVzvHRRx/J3r1729+bTCbp7e0tjx07VqauyWSSHh4e8sSJE/ZtkydPlrNmzSpTt7CwUC5atEhqtdoyZVu2bJGhoaGVxmY0GuXevXsrrYeTZr1UV/YVyEi6xLfvvcnyF2Zy7uA+AJq370BASHMXR+a4wBahhHfpDsCBzV9TVKjmvldcLzIyEo1Gw7Rp09i0aRPp6VUbIjxmzBi2b99uf++sdWhPnjyJRqMhMvLamhSdO3cuU9fPzw9vb29mzpzJSy+9VKWYr3fw4EEKCgqIiIi4qf1vhkr21ykw57FjzXKWPvc4x3duAynx1GrpPW4S9/1tTq14avZmdB8xGrCtWXt857ZKaitK9fPx8WHHjh0IIXj44YcJCgpi5MiRXL58+Yb7hYSEkJaWZn/vrHVoq1r36mpbixYtomvXrjeMtTxZWVlMmTKFuXPnljlfdVLj7Es4vec3tv7nA0xpVwDQeHjQddi93DZyLDqfmvtHqQ4tO3QmKKwVKXFn2ffNl0QPHFJn/3Ap9Ue7du1YunQpAMePH2fy5Mk888wzrF69usJ9EhMTCQgIqNLxHVmH1pG6er2exx57jKCgII4dO0bjxo2rFE9eXh733nsvvXr14sUXX6zSPs6iruyLpSdd5Ou337An+jY9+jDj/33AwMkP1flED7aRObcVX92nno8nrrhrSlFqi6ioKKZPn17pzdcvv/yS/v37V+mYjqxDGxkZSVFREadOnbJvO3ToUIVr1lqtVnJzc0lMTKxSLPn5+dx33300a9aMjz76qEr7OJNK9sWO/vITVosFN40742a/xsjnXsK3cZPKd6xDInv3xxgYBMDv/13n4miUhu748eMsXLiQCxcuAHD+/HlWr15Nr169ytS1WCycO3eOmTNn8vPPPzN37twqncORdWj1ej1jxoxhzpw55OTksHPnTjZs2MCUKbYV4LZs2cKBAwewWCxkZWXx17/+FX9/f9q1awfYkr/ZbKawsBApJWazmYKCAgAKCwsZN24cWq2W5cuX4+ZW86lXJXtsI5KO7/gZgPCuMfY54usbjbs7McVX9xeOxnLx5HEXR6Q0ZEajkd27d9OzZ0/0ej29evWiQ4cOLFy40F5n165dGAwGfHx8GDRoEFlZWezZs4eOHTva60RHR7Ny5cpyzxEUFMT69euZPXs2/v7+7N69mzVr1tjL//GPfzB8+HD7+/fff5+8vDwaN27MpEmT+OCDD+xX9hkZGUyaNAlfX19at27N6dOn2bx5M97etmlIfvnlF7RaLXfffTcJCQlotVqGDh0KwK+//so333zD999/j5+fn/3ZgZI3mqubWoMWuHjyOKtfeR6Ae5/9O5G9+rk4oupTaDbz76cewpydRcRtvRj1/MuuDklRlBtQa9A60bEdPwG22SxbdatbD005ysPbm6533QPYbkhfuXDexREpilITGnyytxQVceJX21epyF59cff0dHFE1a/LXSNw97Q9pbjn6/UujkZRlJrQ4JN9/OED5GXbhlu16zfItcHUEJ2PLx2H2PoSj23/mczkG49rVhSl7mvwyf7odlsXjiGgEc3bd3BxNDUn5p7RCDc3rJYiPv37Xziw+WusFourw1IUpZo06GRfkJdrX5Akqu9A3Nw0Lo6o5vgENabvhMkgBPk5Ofy45CM+nfUXju/cRkFerqvDUxTFyRr0E7R7v/mSooJ8oOF04ZTUc/QEQjt15celH3Hp5HFSz8fz7XtvovHwIKxzNzoMHkpETO1Za1dRlJvXIK/spZRsW/Efdq2zPZLdpHUbgkLDXRyVazRp3YZJr/6L4U/+Fb/gpgBYCgs5s3c3G95cwN6vv3BxhIqiOEODu7K3FBWx5d//y5FtWwFo1LwlI5+f3aDniRFubrQfcDvt+g8mJf4cp3bvJHbbVkxXUtm24j+4e3rRpXi4pqIodVODSvYJsYfYvnoZSadtK92EtG3P6Bfm4G0wuDiy2kEIQeOwVjQOa0WnO4fz2by/k3k5ia3/+QB3Ly86DLrD1SEqinKTGkSyv3T6BDtWLyMh9rB9W6tutzHimVl4eHm7MLLayxgQyPiXX+ezeX8n+0oK33/4HkX5+XQeeneD/hakKHVVvZ8uYf+mr/l52WKktAKg9w+g15j76TTkLtw0DWf0zc1Kv5TIZ/P+Tk6GbWGJ9gNu544/P6H+SCpKDXHWdAn1NtlbLRZ+WvZvDn73LQDeBiM97xtP57vuwcPTq5K9lZIyki7x34Wvk5IQB0BQyzBGPjcbvyZNXRuYojQAam6cG8jPzeHL/3nVnugDW4Qy5Z/vEnPvGJXob4Jfk6ZMeu0t+/DUlIQ4Vs/5G2kXqzaPt6Iorlfvkn1KQhwrXnyGuEP7AQjv0p3757+JT1DVVpJRyufh5c3wp57j9hmPIoQbuZkZrH1ttppqQVHqiHqV7I9t/4lVs58jI+kSAN2Gj+S+F+bgpdO5OLL6QQhB12H3ctfjTwNgupLK2tdmk52W6uLIFEWpTL0YjZOblcn2VUuJ/WkLAO5eXgx9ZGaDfCq2JkQPHEJRQT4/fPw+mZeT+GzuLG4bOZaovgPx0uldHZ6iKOWo0zdorRYLB7//ll8/X0l+bg5g618e+dxsglqGOfVcSln7vv2Kn5d/bH/v7ulFZK++tO7eg+btO9aLtXsVxdUa/GichNhD/LjkI65cSLBtEIKOtw9l4OSH1NVlDTqxazv7vvmKS6dPlCkLCmtFeJfuRPbsS+Pw1mp8vqLchHqf7AvycslMvkxGchIFubn4BAbh27gJVquVX1b+h1O7f7XXDWnbntunP0Jwq4jqCF2pgpSEOP748TtO/bYTU3pamXKfoGBCO3ZGa/TBS2/A22C0/ZsGN8EnMAiNu4cLolaU2q/OJftuXbvKL5b9h8zLSWSmXCbzcpJ90RAAKa3k5+aSn2PCbDJVaZpdg38AAx6cQVS/QeqqsZaQUpJ+KZGE2MPEH95P3KED9plFKyKEG35NmtI4vDXB4a0xNgokPzcXc46JQrMZT60Wb70BL70ea1ER5hwT5hwTVosFL50eL70eb50B4VavxhsoCv5NQwhsEVq3kn1020j5UJdIpxxL4+5O9xGj6Tl6Ap7eWqccU6kehWYz5w7u5dTvu7hyPh5zbg75OSYK8vJcHZqi1Ho9Ro1jwIMznJLsa2w0jsa9+FRCYPAPwLdxMHq/gGtX5ELgpdXhZTDgpdOj9/PHr3ETfBoH46XTk5WaTGbyZfKyMmkR3Qm/4CY1FbpyCzy8vYns1Y/IXv1KbS8qKLD/m2ZcvkRqfByXz50hNeEclqIiwLYAvKe3NwVmc6lvesLNDS+9AY1GQ35ODkWFBTX6mRSlLqqxK/vu3bvLLV9vwCewcYNY1Fu5OZaiIgrMeXhpdaXmLrJaLOTn5aLRaPDw1pbqtisqKKAgL5ea+l1WlJri4eWFl05ft67shRAEhDSvqdMpdZTG3R2twVhmu5tGU+52AHdPT3UBoSiVqLEreyFENlB2fF7DFAiox05tVFtco9riGtUW17SVUpZ/peOAmnyC9oQzvorUB0KIvaotbFRbXKPa4hrVFtcIIZzygJIaq6YoitIAqGSvKIrSANRksv93DZ6rtlNtcY1qi2tUW1yj2uIap7RFjd2gVRRFUVxHdeMoiqI0ACrZK4qiNAAq2SuKojQATkv2QognhBDnhBBmIcQ+IUT/Sup3FEJsE0LkCSEShRBzRD2ZutKRthBCeAshlgohDgshCoUQP9dgqNXOwbYYJITYIIS4JITILW6Th2oy3urkYFu0F0L8JIS4XFz/rBDiH0KIevGosKP5osR+bYQQ2UIIU3XHWFMc/L0IE0LIcl7DKj2RlPKWX8BEoBB4GGgH/C9gAlpWUN8HSAI+BzoAY4Fs4DlnxOPK1020hR74EHgE+Ar42dWfwYVt8RLwGtAXaAU8DhQBD7j6s7igLSKA6UBnIBQYCVwG/uXqz1LTbVFiP09gH/AtYHL153DR70UYIIG7gCYlXp6VnstJAe8GFl+37RTwRgX1HweyAG2JbS8DiRSPEKqrL0fb4rp6i+pZsr/ptihR/3Ngvas/Sy1pi/8H7HL1Z3FVWwBvA0uK/wjWl2TvaO68muxjHD3XLXfjFH+t7A58f13R90CfCnbrDWyXUpac1Pw7IKT4w9RJN9kW9ZIT28IHSHdWXK7gjLYQQkQAw4Btzo2uZt1sWwgh7gFGAH+pvuhq1i3+XnwhhEgWQuwUQoyryvmc0WcfCGiwfcUs6TK2rxflaVJB/atlddXNtEV9dcttIYQYAQyh7j9gc9NtIYT4VQhhxna1twNbV1dd5nBbCCGaAouBKVLK7OoNr0bdzO+FCXgemADcDWwFPhNCTK7sZM6cCO36p7NEOdsqq1/e9rrI0baoz26qLYQQfYFVwF+klL9XR2AucDNtMREwYuu7fxOYBbzh/NBqnCNtsQL4QEr5W/WG5DJVbgspZSqwsMSmvUKIQOAFbO1UIWck+1TAQtm/RI0p+xfrqqQK6nODfeqCm2mL+uqm20II0Q/YCMyRUn5QPeHVqJtuCynl+eIfjwohNMDHQog3pZRFzg+zRtxMW9wODBRCzC1+LwA3IUQR8ISUsq5+83NWvtgNzKis0i1340gpC7DdIb/zuqI7gV8r2G0X0F8I4X1d/YtA3K3G5Co32Rb10s22hRBiALAJeFVK+U61BViDnPh74YbtAk1TWcXa6ibboiPQpcRrDpBX/PNa50dZM5z4e9EFuFSVEzrjjvJEoAD4M7bhQ+9i61sKLS5/A9haor4vtqv7NdiGXo7BNjqnvgy9rHJbFG9rX/wPtgbYW/xzF1d/Fhf8XgwCcrB1V5QcVhbk6s/igraYAowHorANQ52AbbTaGld/lppui3L2n079GY3j6O/FNOCB4rptsfXfFwDPVnouJwb9BLar8nxsf60GlChbCsRdV78j8AtgxvZXaS51fNjlLbRFHLY+ulIvV3+Omm6L4vdl2uH69qqrLwfbYhKwH9vzJybgCLabs9qajtvVbVHOvvUm2d/E78U04Ci2i6IsbBeHk6tyHjXrpaIoSgOg5sZRFEVpAFSyVxRFaQBUslcURWkAVLJXFEVpAFSyVxRFaQBUslcURWkAVLJXFEVpAFSyVxRFaQD+PzOHScPl0mrrAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Tilt-Corr-NN1: computing NN1 correlation of tilting in non-3D perovskite structure, all NN1 neighbours are considered the same. \n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEUCAYAAADOaUa5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAv/0lEQVR4nO3dd3yc1Z3v8c9PzbIsWS5yb7IxGFNMNWAwJaEvIaSQhJBNwi53s7spu5eFbFgIF9JIbi7ZJHvDvhJINrCEADEQAoQkmOLcYJotYxuDccO9NxWra+bcP84z0kiWJU3XPP6+X695zTNPm9/jkX86+j1nzjHnHCIiEk4FuQ5AREQyR0leRCTElORFREJMSV5EJMSU5EVEQqwoW29UVVXlqqurs/V2IiKhUFNTs885NybZ47OW5Kurq1m6dGm23k5EJBTMbHMqx2ctyYsMJu2RKP/rd+9wsLGNc2eO5ryZVcyoGoaZ5To0kbRSkpej0h9W7eKRN7cA8Md3dgEwbXQZv/j8XGaOLc9laCJppRuvclR6Zd1eAIoLjVjbffP+JhYs3Zq7oEQyQElejjrOOV5Ztw+AOZNH8Hfnz2DiiFIAVm6ry2VoImmnJC9HnU37m9hR1wLAmPIhDC0pZPLIMgBWba8jGtV4ThIeSvJy1ImVagCmjBoKwLiKIQA0tHawaX9jTuISyYRBk+QXLlzIhz70ISZNmkRpaSlTpkzhlltuob29PdehSci8st6XaqaNKqOsxPc9GDu8tHP729tVspHwGDS9a1asWMHFF1/MF7/4RYYNG8ayZcu44447GDVqFLfddluuw5OQiEQdr27YD9CtF035kCKGDSmksTXCym11XHPqpFyFKJJWgybJ33LLLZ3LkUiE8847j0WLFrF48eIcRiVhs3JbLQ0tHQBMHjm027ZxFaW839rI27r5KiEyKMo1HR0dPPDAA8ybN4+qqiqKioooLi7m6aefprxcfZYlfRYHpZqSwgIqhxZ32zZ2uK/Lr9pRR0Q3XyUkBkWSv/766/nKV77CRRddxIMPPsjrr7/OokWLMDNOOumkXIcnIfKXoOvkqVNHHPbt1nEVvi7f1Bbh/b2Hsh6bSCbkvFyzfPlyFixYwMMPP8z111/fuf7xxx/HOcdpp52Ww+gkTJraOli25SAAp04ecdj2WEsefH/5Y8dVZCs0kYzJeUt+yxb/1fJZs2Z1rmtsbOTrX/86AKeffnpO4pLweXPjAdojvgxz3LjDy4BlJUVUlPp2j3rYSFjkPMmfdtpplJSU8NWvfpWFCxfy0EMPcc4559DW1sbYsWOZOHFirkOUkIj1qqkcWsyoYSW97jM26C+/cltttsISyaicJ/kpU6bw8MMPs3nzZj784Q9z77338v3vf5+ZM2eqVCNp9e6OegBOnzqC+qCHTU/jgv7y7+yopyMSzVpsIpmS85o8wLXXXsu1117bbd2VV16Zo2gkrNbv8TdTjx1XQW1TW6/7xFryrR1R1u05xOwJw7MWn0gm5LwlL5INDS3t7Kr349XMqBpGXXPfLXlA/eUlFJTk5aiwYW/XeDTjK0uJut77wZcWF3b2n1+5vTYboYlklJK8HBVipRqAUWW933SNiZVsVmxVS17yn5K8HBViSX5MxRA6jtCKj5lQGbv5WkdDiwbIk/ymJC9HhViSnzmm/Ig3XWMmBWPaRB0s3Xww47GJZJKSvBwVNgTDFMwcW87Bxr5b51XlQygt9v813nj/QMZjE8kkJXkJvZb2CJuDiUCOHVdObXPfSb7AjEkjfGv+jY37Mx6fSCYpyUvobdrfSGxQyemjhw2ozj6h0if5ldvqaGztvbulSD5QkpfQi+9ZM6ZiCP3cdwW6xpqPRB01qstLHlOSl9CLJfmK0iIKB/gTP6Z8CGUlhYBKNpLflOQl9Dp71owtP+KYNT0VFBgnBEMa6Oar5DMleQm9+O6TA03y0DUc8YpttTS3RTISm0imKclLqEWijvf3+Z41M8eWd87vOhDTRg8DoD3ieGuL6vKSn5TkJdS2HWyircMPGXzsuHIOJZDkq+Lq8q9vVMlG8pOSvIRafM+amWMqEhqmoKmtgzOmjQTgjfd181Xyk5K8hNq6IMkPKSpg7PAhNCVQW29o7eCcGaMBeGtrLS3tqstL/lGSl1CLteRnjClP+OZpa3uUU6dUAtDWEWX1zvq0xyeSaUryEmrx3ScPJfHN1UkjyzqX1+5uSFtcItmiJC+h5ZxjQ7fuk0kMG+zoHMdmza5D/ewsMvgoyUtobdrfREPQep89oSKh7pMxDS0dnf3l1ZKXfKQkL6G1cltt5/IpU0Yk1H0ypqG1nePGVwCwRkle8pCSvITWymAi7nHDhzBueCkNrYmXaxpaOpg1zif5vQ2tHGjse8IRkcFGSV5CK9aSP3nSCICkyjWHWjo4LkjyoJKN5B8leQmljkiUVdt9l8dTJlfinEuyJt/OzLHlFJh/rSQv+UZJXkJp/d5DNAdfXpozZQStHdHO4Q0Scai1gyFFBVRX+XFs1uxSkpf8oiQvoRSrxwPMmVSZVCse/OBkLe3Rzrq8WvKSb5TkJZRi9fgpo4YyclhJQmPW9NTQ2t5Zl1+zqwE3kKmlRAYJJXkJpVhLfs7kEQBJfds1pqGlg1lBN8r6lg5217emHJ9ItijJS+i0dkQ6x5k5ZbIfeybZcg0c3sNG/eUlnyjJS+is2dVAe8SXVLq6T6ZQrmnpoHp0GSVFBcH5NVCZ5A8leQmdFUGpxgxOTkdLvrWdosICZo7xwxtoDBvJJ0ryEjort9YCcMyYcsqHFAGpJfnYvLCxurx62Eg+UZKX0Om66epb8c65lG68xsa8idXl1+1pIBJVDxvJD0ryEipNbR2s2+Nb2nMmVQbrIikl5YaWDtojUWaN9+WalvYoWw80pR6sSBYoyUuoLN9aSyyfz5kyAkitVAMQdY5tB5vVw0bykpK8hMoDizcBUFFaxAkThgOp9ayJ2bSvkUkjhnbW+F/XxN6SJ5TkJTRW76zn+Xd3A/C3502ntLgQoHPikFRs3NcIwMWzxwLw8Btb2F7bnPJ5RTJNSV5C4ycvrQegfEgRf3ve9M71qZZrAOqa26ltaufmS2dRXGi0dUT5wZ/WpHxekUxTkpdQWLu7gedW7QTg8+dOo7KsuHNbOso1ABv3NzJ1dBmfPacagN8u386q7XV9HySSY0ryEgo/eWk9zkFZSSE3zp/Rud45x76G9Iw1syko2XzlgzOpKC3COfjuH1ZrwDIZ1JTkJe9t2HuIZ1fuAOCz50xj1LCSzm3r9hziYFN6WvLbDjbT1hFl5LASvvyBmQAsXr+fRWv3puX8IpmgJC95raGlnX99fCVRB6XFBfyP87ta8dGo47UN6esFE4k6th70/eM/f241k0YMBeDrv12lfvMyaCnJS96qa27ns794k5rNBwH4wgXHMKZiSOf293Y1pH3i7c37fcmmtLiQOz40G4Dttc1cd9/rSvQyKCnJS16qbWrjMz9/neXBODU3nFvNTZcc27k9EnW8sTH9fdk37mvqrMFfcdIE7v7oyUBXot+yX4leBpeiXAcg0p+W9gh1ze0caGxj2ZaDLF6/j8Xr91PX7GvtX7hgBv925fGYWecxq3fWU5umWny8+iCO0eX+L4brz54KwG2/fZvttc18+N5XuGT2OC6aNYb5M6uoHFrcLS6RbMtakm9o6WDed1/M1ttJHotEHe2RKB0R5yfgjhx5Au4vfeAYbrlsVrdE2toRyeg3Up9ctp1Z4ys4fnwFYyqGcP3ZUzGDf3vybWqb2nm8ZhuP12wDoMCgrKSIoSWFFBUo2UtiPjtvWsrnsIF0/zKzAuBy4BLgLGA8UAocANYCi4EnnXNb+jhHAxDmb49UAftyHUQG6fryV5ivDcJ/fbOccxX979a7PpO8mVUANwH/AIwG3gFW4P9BW4ERQDVwJv4f+iXgW865v/RyrqXOuTOTDXSw0/XltzBfX5ivDXR9/emvXLMReAv4KvC0c+6IQ++Z2WnAJ4EnzewbzrmfJBuUiIikR39J/nLnXM1ATuScewt4y8y+BUxNOTIREUlZn0l+oAm+xzFNwHu9bLov0XPlGV1ffgvz9YX52kDX16cB3XgVEZH8lJYvQ5nZC2am/pEiIoNMur7xamk8l4iIpInKNSIiIabWt4hIiA14WAMz+6v+9nHOPZdaOCIikk4DLteYWRRw+Pp7vM4TOOcK0xeaiIikKpEByqb3sm4UcBlwA/A36QhIRETSJy03Xs3sFuAC59yHj7RPVVWVq66uTvm9RESOJjU1Nfucc2OSPT5dQw2/BdzV1w7V1dUsXbo0TW8nkgHRKBSoL4IMLma2OZXjU/6JNrMSfLlmZ6rnEsmJSAf85vPwv6fB9mW5jkYkrRLpXbOEuJusgRL8UMMVqCYv+eqlb8G7T/nlml/CpNNzGo5IOiVSrnmHw5N8C7AAeMo5907aohLJlnefhsU/6nq9YRE4B5qyT0JiwEneOXdDBuNg7dq1zJo1i+eee44rr7wyk28l4u1dC0/9o18uKIZoO9RtgYMbYdSM3MYmkiaD5i7TsmW+Fjp37twcRyJHhbYmeOyvoe0QFA6B637dte39RTkLSyTd0jUK5dfN7I5UzlFTU0N1dTVVVVXpCEmkb8sfhn3BlMNX3QPHXQbjTvKvleQlRNLVkr+LfrpQ9qempoa5c+dy3333ccIJJzBs2DDmzZvX2cIXSZtoFN74mV+ecjac/jm/POMi/7zx//l9REIgLUneOVeU6pAGb731FosXL+aZZ57hnnvu4YEHHmD37t1cffXVNDY2piNMEe/9l2H/Or989t93rZ9+oX9uPgi7VmY/LpEMSNeXoVKyYcMGamtrOfvss3n66aexoGdDRUUFV155JYsXL+ayyy7LcZQSGrFWfMUEmB33Je1p50JBEUQ7fMlm4qm5iE4krZJqyZtZgZmV9XwkG0RNjZ9K9jvf+U5nggc4+eSTAdi7d2+ypxbpbv8GWPe8Xz7zRigs7to2pBwmn+WXVZeXkBhwkjfva2a2HmgHGnp5JGXZsmUcc8wxnHHGGd3W79zpv0Q7efLkZE8t0t2SnwMOCkvgjBsO3x6ry295DdpbshiYSGYk0pL/J+BW4Bf44Ya/A3wTWAtsAr6QbBA1NTVMnDjxsPWPPfYYlZWVzJs3L9lTi3RpPQRv/covn3QtlAdjPkU64MD70Li/K8l3tMDWN3ISpkg6JVKT/zvgTuBefIJ/yjm3zMy+BTwDHJtsEMuWLaO0tJSOjg6KinxI27Zt46c//Sk333wzJSUlyZ5apMvKR6G13i+f/HHYXgMHN8GBjRBphxkXwuS5UFIBbQ2w8c9+nUgeS6QlPx1Y7pyL4Ms1IwCcc1HgP4HPJxPApk2bOHDgAEOGDOGGG27gxRdf5MEHH+T888/nlFNO4bbbbkvmtCKHW/+Sf66YAFuXwNrn/bdeI+1+fcNOX6Ovnu9fv//n3MQpkkaJJPn9QHmwvAU4LW7bSGBoMgHEbro+++yzNDc3c/XVV3PrrbfysY99jD/96U9qxUv67F3tnysm9L69YZd/nnqOf971dtcvAJE8lUi5ZjEwF3gO+DVwl5mNAtqALwEvJhPAxz/+cWITlzzxxBPJnEKkf21NviwDMOwI36puqfd1+wmn+NeRVti7BsaflJ0YRTIgkSR/FzApWL4bX665Ad+CXwh8JY1xiaTXvrV0DqJa1sfQGYd2dyV5gJ0rlOQlrw24XOOcW+OceylYbnXO/bNzbpJzbpRz7lPOuT2ZC1MkRXtWdy0P62MmtYadUDYKKqf61ztXZDYukQwbNKNQimRUrB5fUg7Ffdw+itXlJwat+Z3LMxqWSKb1meTN7EUzu2SgJzOzsWb2TTNT6UYGlz3v+eey0X3v1xDMYhkr2ex6G6KRzMUlkmH91eSfAx42sxbgSeBVYBWwD2jF1+WnA2cAVwIXAn8Cbs5QvCLJiZVr+irVgL/x2toAE071r9ubYP96GDMro+GJZEqfSd459wMz+ylwPfA5fC+anqNNGn4S7yeBrzrnlmcgTpHktTb4GZ/gyD1r4jXsOvzmq5K85Kl+e9c45xqB+4H7g0HITgHGA6XAAWCNc25TJoMUScneNV3LffWsiWnYCVXH+v70DTt9kp/zyczFJ5JBCQ017JxrAl7LUCwimdGtZ80AW/LgSzYNO2HH8kxEJZIV6l0j4bc37qZrUWn/+zfsBOfibr6u1ExRkreU5CX89rzrn4cfPtJpr9qa/EBmsSTfWg8HN2YmNpEMU5KX8It1n+yvZ0283m6+iuQhJXkJt+ZaaNjhl4eOGvhx9Tt8yz/2i0FJXvKUkryEW6weDzBk+MCPa9oPZl2teX3zVfJUItP/nW9m18S9rjKzX5vZcjP7gZkV93W8SE7E96wZOmLgx7XU+ufOJL/C34wVyTOJtOS/D8QPx/dj4GLgdfxolN9IX1giaRJL8pVT/LyuA9VS172HTfNBqN2S/vhEMiyRJD8LqAEIvhT1UeCfnXP/APwr8Kn0hyeSotjAZCOnJ3ZcR5uf53Xi6V3rti9NX1wiWZJIki8BYtPXn4f/ItXvg9drgSNMtyOSI851teRHTEn8+OZaqJwM5eP9621K8pJ/Ekny7wFXBMufAV5zzjUEryfihzgQGTwOboTGvX45Nj58Ilrq/M3XyWf619uWpC82kSxJJMl/E7jJzPbiByz7Xty2K4C30hmYSMo2vdK1PHpG4se31PnnKWf5550roKM19bhEsmjAY9c45542s9n4Cbzfds6tjdv8GqCOxDK4bFrsn0cfCwUJDdPkxXrYTJ7rnyNtsHMlTJmblvBEsiGRLpSfA+qcc0/0SPAAj+NvzIoMDs7B5iDJV8+HtkOJnyPWkp9wKlgwwrZKNpJnEinX/BI45gjbpgfbRQaH2s1Qt9UvV8/3k4EkKpbkS8q6JvNWkpc8k0iStz62jQbqU4xFJH1ipRqAaecm2ZKv7foCVKxkox42kmf6LFQG33C9Jm7VHcGN13ilwPmAmjgyeMRuuo46BkpHJDdPa6QD2hphSLlP8kt+7meYatgFFePTGq5IpvR3N2oscHLc62Pws0LFawOeB76dxrhEUrM5SPLJ1uNjWuq6knzMtiUw++rU4hPJkv7meL0fP/UfZvYy8I/Ouff6OkYk52q3dA1BUD3fz/GarJY6qJwEo2b4USybDyjJS14ZcE3eOfcBJXjJC93q8eel2JKv9c9mqstLXkqo87CZTQQ+BEzG1+LjOefc19IVmEjSYqWakdN9K7x2c/LnivWwAZ/k1/0Jti/z9frCJPrei2TZgH9KzeyjwCNAIbAHX4uP5wAlecm92E3X6vP8czLdJ2O6JflgeIOOZti9Ciaemvx5RbIkkabI3fgbrDc45zROjQxOddvg4Ca/XH2+f25LoSbfXNu1POkMfE9i5+vySvKSBxLpJz8F+A8leBnU1j3ftTwtaMnHJ+pEtdZDNOqXS4fD+KCz2do/Jn9OkSxKJMm/ioYukMFuxaP+eeLpfnjh9uaukSiTEY10v3E7+8P++f1F0KT2jgx+iST5fwG+YGafN7OJZlbW85GpIEUGZP8G2PqGXz7l0/754ObUp+2L9bABOPEj/jnaAWv+kNp5RbIgkSS/Ev/FqF8CW4GGXh4iuRNrxRcUw0kf98ux+nwq4m++Vh0LY0/wy+8+lfq5RTIskRuvf4vvQSMy+ESjsDJI8sddDsNG++V0J3mAEz4Ce96FDS/7en8iE4SLZFki48k/kME4RFKz5dWub7mecp1/bq71E3CnqueN2xM/Aovuhmg7rHkOTr0+9fcQyZBEyjUig9eKR/zz0JFw7OV+OR2teDi8JT9mFoyZ7ZffeSo97yGSIYl8GWov/ZRrnHNjU45IJFFtTfDO7/zySddCUYlfTuWbrvF6JnkIWvOrYcNLfntpZXreSyTNEqnJ38vhSX4U8EFgOPCLdAUlkpD3ft/1hadTg141zqWvJd9af3giP+EaWPTdoGTzh64Skcggk0hN/q7e1puZAb8BOtIUk8jAtTb4ZAtQdZzvHw++b3xbU3rewznYugSOvaRr3djZUDUL9q2B5b+GOZ/yg5iJDDIp1+Sdcw74OfDl1MMRSdBzX4UDG/zyRbd2Jdp0teJjdr51+C+Nkz/hnzf+uesXjcggk64brzOAkjSdS2RgVjzWdcP11L/u6hsP6U/ykQ7Ysaz7unO/3DX88J//t2/Riwwyidx4/WIvq0uA2cBngAXpCkqkX/s3wO//xS+PPhb+6vtd26KRru6U6bRtKUw5GwqL/eviofDpR+HnF/tfKk9/BYZPghkXpv+9RZKUyI3Xn/SyrhXYBvwn8I20RCTSl/Zm33p/5Yd+TJnCErj2v6BkmN8eCfquR9oz8947V8LkM7rWDauCzzwOP7/ED3/wyKd9C3/el9TjRgaFRG68plbacVHf+hLpyTk/FkykzSfnjhY/gXbbobjnQ9C4D95eAE37u4697DswYY5fbj4Iq56EQ3syF+vWN2DiaVAQ99+h6li47tfw0EehvdGXbt68D875Eow70ffdHzrC/0ISSUTpiJRPYS7VwZsG6MyTZrql16YwGqBIvElnwnn/DDMv9pOCNO6F9QuhvSXz7z3+ZBhZDRXj/byvsYS/+1146duw5veZj0GODvNvwi79Ro1z7sxkT5FQkjezGcBXgfn4PvIHgL8A9zjn3u/n2AZgTbKB5oEqYF+ug8ggXV/+CvO1Qfivb5ZzriLZgwec5M3sDOBloAV4FtgNjAOuws/3+gHn3LI+jl+aym+jwU7Xl9/CfH1hvjbQ9fUnkRuv9wBvAVc65zo7DAfjyD8XbP9gsoGIiEj6JXIz9Szg+/EJHiB4fQ9wdjoDExGR1CWS5JuB0UfYNgpfxunLfQm8Vz7S9eW3MF9fmK8NdH19SqQm/yBwKfBJ59wrcevnA48BC51zN6QSjIiIpFciSX408DtgHrAXf+N1bPB4FfiIc27/kc8gIiLZlnA/eTO7ApgLTAB2Am84557PQGwiIpKirH0ZSkREsq/PG69mNtrMnjCzy/vY5/JgH80KJSIyyPTXu+Z/4ocR7qsc8zwwHbg5TTGJiEia9JfkPwn81PVR0wm2/Qy4Jp2BiYhI6vpL8tOAdwdwntVAdcrRiIhIWvU3rEEzfpLu/pQH+x5RVVWVq66uHmBYIiICUFNTs885NybZ4/tL8suADwP9jZ16TbDvEVVXV7N06dIEQhPJnsb2RjbXb2b2qNmYJuSWQcTMNqdyfH/lmnuBG83s830E8Dngb+h95iiRQW9/834+/vTH+dSzn+LFLS/mOhyRtOqzJe+ce9LMfgz80sy+DPwR2AI4YCpwOXAm8EPn3G8zHaxIurVH2vmXRf/C9kPbAXhl+ytcMu2SHEclkj79DjXsnLvZzBbhu1PeAgwJNrUCi4FrnHPPZipAkUxxznH3m3ezbE9XpXHtwbU5jEgk/QY0nrxz7hngGTMromskyv3OuY6MRSaSYY+teYzH1z4OwLDiYTS2N7Lu4Doi0QiFBYU5jk4kPRKanNs51+Gc2x08lOAlb22q28T33vweADNHzOSueXcB0BJpYUvDlhxGJpJeicwMJRIaL2x5gYiLUFRQxH988D8oLiju3Lbm4BqmV07PYXQi6ZNQSz7TFixYwIUXXsjw4cMpLy9n7ty5/P73/fXeFEncX7b9BYDTx57OlIopjCsbR+WQSgDWHlBdXsJj0CT5W265heuuu44TTzyRX/3qVzzyyCOce+651NfX5zo0CZn6tnpW7F0BwPxJ8wEwM2aNnAXo5quEy6Ao1zz66KP84Ac/YMGCBVx77bWd66+++uocRiVh9fqO14m4COCTvHOO1QdWM65sHODLNSJhMSiS/J133snVV1/dLcGLZMor2/3slePKxjG1YioLNy9kfe16xpT5b47vatxFXWtdZ/lGJJ/lvFzz/vvvs3btWj7xiU/kOhQ5CjjnWLx9MQDnTDiHBWsXsL52PQAjh4zs3E8lGwmLnCf5HTt2ADBx4sQcRyJHg7UH17KneQ8AVUOrqG/ruucztHgohVbYuZ9IGOQ8yceS+zvvvJPjSORo8JftvldNoRVidB+IrLG9sbPr5JoDqstLOOS8Jj99+nTmz5/PnXfeCcCcOXPYvXs3zz33HLfffjvHHXdcjiOUMInV46cNn0ZxYXG3bS0dLRwz4hjW167XzVcJjZwneTPjiSee4Pbbb+eee+5h9+7djBs3jgsvvJDp0/WFFEmfhrYGlu9ZDsDYst6nJJ5SPgWA9QfX0xHtoKgg5/9FRFIyKH6Cx44dy/3335/rMCTkXt/Z1XVyasXUXvcZN8x3o2yLtrGlfgszRszIWnwimZDzmrxItry641XAD0Y2qnRUr/uMHjq6c1klGwkDJXk5aqzatwqACcMmHHH2J+dcZ1dK3XyVMFCSl6NCS0cL6w/6/vBHqseDH/LguFH+Zr9a8hIGSvJyVFh7cC0dwejYY4YeeU7k2tZaZo+aDfiWf9RFsxKfSKYoyctR4Z39/nsYhvWZ5NsibZxcdTLgE37s27Ai+UpJXo4K7+zzSX5k6cjD+sf3NKNyBgXm/2ss2bUk47GJZJKSvBwVYi35vlrxMREX6Rx2eOmupRmNSyTTlOQl9Jram3i/7n2AzpEm+1LbWsvc8XMBWLp7qerykteU5CX01hxc05moB9KSj0/yqstLvlOSl9CL9Y8vsAKqhlb1u39dax2njzu9cwAz1eUlnynJS+jF6vGTyicNaCyautY6KoorOH7U8YDq8pLflOQl9GI9ayaVTxrQ/u3Rdpo7mlWXl1BQkpdQO9R2iE31m4C+v+nak+ryEhZK8hJqqw+s7lw+0qBkvVFdXsJCSV5CLVaqKS4oZmjh0AEfV9tay/CS4Z11+ZrdNRmJTyTTlOQl1GI3XWeOmAm9DzzZq7rWOoCuuvwu1eUlPynJS6jFkvyxI49N6LjYBN9njjsTgIOtB1WXl7ykJC+htbtxN1sbtgJQPbw6oWPrWutwznWryy/bvSzdIYpknJK8hNabu97sXE60Jd8ebaepo4nKIZWdx6ouL/lISV5C6/WdrwN+PtfSotKEj4/V5c8YdwbgW/LOufQFKJIFSvISSs65ziR/9oSzOdR2KOFz9Ezye5r3sK1hW/qCFMkCJXkJpU31m9jTtAeAcyacw6H2JJJ8W/ckD/7bryL5REleQinWijeMs8afRWN7Y8LniLXkq4ZWMW34NEB1eck/SvISSm/sfAOA40cdT+WQypTKNdDVmleSl3yjJC+hE4lGOnvWnD3hbJo7mom4SMLnqW+r77zRGkvy2w5tY3fj7vQFK5JhSvISOqsPrKahrQHw9fhkSjXgJ/Vu7mgGutfll+1Rf3nJH0ryEjqxenxRQRGnjT0tqZuuMbGSzcRhExk/bDygko3kFyV5CZ1Ykj9lzCmUFZclVY+PifWwMTPV5SUvKclLqLRGWlm+ZzngSzVAWlryAKePPR2A9bXrqW2pTfqcItmkJC+hsmz3MlojrUBXkk+2Jg/dk3xssDJQXV7yh5K8hMqLW14EYHjJcE6sOhGg8yZsMmLlGoDpldM7Jx55dcerKUQpkj1K8hIakWikM8l/cOoHKS4opj3a3vnN12TERqMEX5e/YPIFgP9lEokm3i1TJNuU5CU0lu9dzr7mfQBcOu1SAHYc2pFUH/mY+G6U8efd17yP5XuXJx+sSJYoyUtovLD5BQDKi8s76/Fb6rekfN7YBCIA8ybMo6K4AoDnNz2f8rlFMk1JXkIh6qIs3LwQgIumXERJYQkAWxpST/LxN1+LC4v5wNQPAP6XiqYElMFOSV5C4e19b7O7yQ83ECup1LXWdUvQyep5jsumXQb4oYdX7F2R8vlFMklJXkJh4Sbfii8rKuPciecCdE79l6r4HjYA8ybOo7y4HFDJRgY/JXnJe845Xtji6/EXTr6wcxaodNTj4fCWfElhCRdNuQiAhZsXqmQjg5qSvOS9dw+8y/ZD2wG4ZNolAHREOzrXpaq2tfawRB4rCe1u2s3b+95Oy/uIZIKSvOS9J9c+CUBpYSnzJ80HYGfjTtqj7Wk5f1ukjfW167utO2/SeZQVlQEq2cjgpiQvee1363/Hb9b+BvBfgCor9ol3a3166vExS3Yt6daaH1I4hAunXAjAo+89ystbXk7r+4mki5K85K0lu5Zw12t3ATClYgq3nnVr57Z0dJ2MV9dax7qD67qtu/GkGykvLqct2sZNi27imQ3PpPU9RdJBSV7y0ub6zdy06CY6oh1UlFRw78X3MrJ0JAB7mvZwoOVA2t9zya4l3YYymDVqFv91+X8xqnQUERfhtldu45erfklHtCPt7y2SLCV5yRtRF2XJriXc+eqdXPfsddS11lFkRfzwoh8yvXI64Ms0T294OiPvX99Wz5qDa7qtmz16Ng9e8SAThk0A4N9r/p2rnryK/37nv1Ma/VIkXSw2+FKmzTl9jrvj0Tuy8l6S39oibdS31lPXVkdtay0Hmg+wv2U/u5t2H9ad8RvnfoOPHfsxAN478B4vb32ZTP5MV5RUcP3x11NYUNht/a7GXfzTS//E6gOrO9eVFZVx3MjjmDFiBtOHT2dYybCMxSXhdPzI45kzdk6Nc+7M/vfuXdaS/HFzjnNDbh6SlfeScKssqeSK6Vdw9Yyrqa6spqm9iQMtB3hn/ztZef8pFVMYVzaO4UOGU1lSSVFBEeD/0li8fTEPvvtg50TiIqm48aQbuenMm/IjyZtZA7Cm3x3zVxWwL9dBZJCuL3+F+dog/Nc3yzlXkezBRemMpB9rUvltNNiZ2VJdX/4K8/WF+drg6Li+VI7XjVcRkRBTkhcRCbFsJvn7svheuaDry29hvr4wXxvo+vqUtRuvIiKSfSrXiIiEmJK8iEiIKcmLiIRY1pK8eX80M2dm1/bYtilYH//4XrZiS4WZ3W9mG8ys2cz2mtnvzGx2j31GmtlDZlYXPB4ysxE5CjkhZvYFM3vZzGqDz6W6l33y+fP7opltNLMWM6sxs/NzHVMyzOwCM3vazLYH//439Nj+QC+f0es5CjchZvYlM1tpZvXB4zUzuypuu5nZXWa2I/h/uMjMTsxlzKkws9uCz+cnceuS/vyy2ZK/GYj0sf2bwIS4x7ezEVQaLAVuAGYDlwMGvGBmxXH7/Bo4HbgSuCJYfii7YSatDHgeuKuf/fLu8zOzTwE/Bu4GTgNeBf5gZlNzGlhyyoFVwD8DzUfY5wW6f0Z/lZ3QUrYN+Br+/82ZwEvAU2Y2J9j+r/j88hVgLrAHWGhmSX9LNFfM7Bzg74CVvWxO7vNzzmX8gf9gtgJjAQdc22P7JuCWbMSShWudE1zjrOD17OD1eXH7zI/fJx8ewWfogOpetuXl5we8AdzfY9064Lu5ji3F6zoE3NBj3QPAs7mOLY3XeAD4e3yjaidwe9y2oUAD8Pe5jjPBa6oENgAfBBYBP0nH55fxlnzw2/SR4B98Tx+73mJm+81suZndbmYlmY4t3cxsGPA3wBZ84gOYh/9P92rcrouBRuDcbMaXYXn1+QXxnYH/KyXe84Trc4k338z2mNnaoMw4NtcBJcrMCs3sOvxfLq8C04HxxH2Ozrlm4P+Rf5/jfcDjzrmXjrA9qc8vG2PX/BT4o3PuuT72+Q/gLWA/cBbwPfyH9z8yH17qzOyLwPeBYfhB2C52zrUGm8cDe13w6xjAOefMbE+wLQzy8fOrAgqB3T3W7wYuyX44GfdH4ElgI1CNL6e9ZGZnxP2sDlpmdjLwGlCKbzR91Dn3tpnFEnlvn+OkLIaYEjP7O2Am8Nkj7JL055dUkjezbwO397PbB4ApwCn4P/WPyDn373EvV5pZPfCYmX3NObc/mRhTMdDrc84tCpYfBhbi62S3AAvM7DznXFOwvbdvnNkR1mdcEtfXp8H2+SWo52eQs88lk5xzj8a9fNvMaoDNwFX45DHYrQFOBUYAHwceNLOL4rbn7edoZrPw94XOd8619bZPKp9fsi35HwG/6mefLfgbkicAh8wsfttjZvaac27+EY59I3ieiW8dZtuPGNj1AeCcqwPqgHXBHe+D+B/Eh4BdwFgzs1hr3vw/xhgOb31ky49I4PqSkOvPbyD24TsC9Pxraiy5+1yyxjm3w8y2AcfmOpaBCJLf+uDlUjObC9wEfCdYNx5/3y8mnz7Hefi/LFfF5clC4AIz+wdgWM/WeiKfX1JJ3jm3jwGM32xmtwP39Fj9Nr61+7s+Dj01eN6ZTHypGuj1HYEFj9gMKa/h64fz6KrLz8OXdl497OgsSPH6BuLU4Dknn99AOOfagtbQpcCCuE2XAk/kJqrsMbMqfDlj0H5G/SjA/x/biG9IXQosATCzUuB84Ks5iy4xT+F76cX7Jb4TwN3AYa37RD6/jNbknXPbge3x64LfVFudc+8Hr+cB5wAv41vDc4EfAk8751JpTWacmc3Et9hfAPYCk4FbgVbgWQDn3Goz+yPws6DuZsDP8HfKB/0kKmY2Ht9KOi5YdULQx3+Lc+5APn9+wL8DD5nZm/ib4f8ATMTfR8orZlaO/8sJfAKcaman4nuhHMB3gX0CnxSqge/iuxr+NsuhJiz4zsXv8S31CuB64CLgquD+1o+A283sPWAt8HV83f7XOQk4Qc65WqA2fp2ZNQIHnHOrzKzczO4i2c8vB92EunWhxPd9fT24yGbgPfwPZFmuuzQN4FqmAH8I/rHb8D+EDwPH99hvFL48Uh88fgWMyHX8A7zGu4LPrOfjhnz//IL4v4jvCdUK1AAX5DqmJK/joiN8Tg/guxT+Ke7ndHOwfkqu4x7gtT0QxNwaXMMLwOVx2y34mdsJtAB/Bk7KddwpXvMigi6UqX5+GoVSRCTENHaNiEiIKcmLiISYkryISIgpyYuIhJiSvIhIiCnJi4iEmJK8iEiIKcmLiISYkryISIgpyUva9DIHZW+Pi4L5KpfGHdfz9Sd7zlE6gPc283O1umBMoYzpGW+K5zIzW2Fmn49bV2Bm95pZg/n5TGfEbbvXzH6RjveWo0M2Jg2Ro8e8uOWh+Lk4v40fXCrmXfwYP0P7OM8n8UOvPpDge1cHy9eRB3PMBj4JjKT7YFqfwc85ew1+ApOfApcF2/4P8J6Zfdc5tx6RfijJS9o45zpnjw9GRQTYEL8+UJ+Bt/80fkrFVcFyviT5fwIecs61x607B7jb+WngXjKzzbENzrlNZvYK8I/4yatF+qRyjWRdX+UOM3sAP3zzhXElnrv6OV8h8AngaeC/8MMhz+ntPc3sUjNbaWaNZvaKmZ3YY78vm9nWYPtTZnZxrMzUTwzzzezPZtZkfq7b+83Pb9zXMTPx85A+3mPTRuB6MxtlZp/AT0IT7wngM2am/7/SL/2QyGDzLfzY9G/hSzDzgJ/3c8wHgXHAo/iE2Y5vzfc0FV/u+E6wfSzwm2CmLszso8D/xf+y+CiwEui3/m1m5wEv4ievuBb4n8Bf4Sd+6MvF+L8+VvRY/zP8jD/7gfuD88V7FX+9J/cXm4jKNTKoOOc2mNkBoKCXMs+RfBo/nv0fnZ/xaSFwnZnd5rqPpT0KOM85tw78DU78pAuz8OPg3wY855z7UrD/88EMPP/Yz/t/D3jVOfep2Aoz2w68aGYnOedWHeG4M4DVzrlo/ErnXIOZnQ0cA+xwzjX2OO4d/NSFZ3H4LwiRbtSSl7xmZkPwre7fuq5JkB/B34Q9p8fum2IJPvBu8Dw5KPmcim/Fx+v5uuf7l+H/2viNmRXFHsAr+L8ozujj8PEcYRpG51zUObeulwSPc64D/0ut5/y0IodRkpd8dyUwAnjOzEYEUxMuws8i1LNkU9vjdeyXQil+YvUi/DSO8Xq+7mkkftLl/8Qn9dijFSjGzx52JKXBfsloDY4X6ZPKNZLvYol8QS/bPmlmNznnIgM4z16gA5/s4/V83VMtfpq9u4Dnetm+o49jD5B8a3xEcLxIn5TkZTBqYwCt1KCb5ofw5Zn7emw+DT9R9wfwc4L2yTkXMbPl+L7pP4vb9OF+jms0s9eBWc65b/b3Pj2soft3CwbEzMYAZfhJq0X6pCQvg9F7wDVm9hFgG/7mY28t4mvwye7Hzrk34jeY2WLgdnxLv98kH7gbeNLMfoKvxZ8HXBVsix7xKPhX/E3WKL53TwO+J89VwO3OuSMl48XA/zKzMc65/spC8c7E//XwagLHyFFKNXkZjP4TeB7f530J8IUj7PdpYF3PBA8QfLnoN8DHgpuz/XLO/Rb/5aSPAE8Bc4Fbgs1H/AKXc+4V4AJ8aech4Bl84t8K7O7jLRfhSy5XDCS+OFcAf3bO7U/wODkKWfceZiISz8y+jv+LYJRzrjkD5/8xMNM5d1W/O9P5xa/NwK3OuV+lOx4JH5VrRAJBrfvf8F/GagLOB74G/CITCT7wf4A1ZnZcH2WdeJ8AmvFf/BLpl5K8SJc24Hjgc0AlsBP4MXBHpt7QObfNzG4EJjCwG6kG3Bj0lRfpl8o1IiIhphuvIiIhpiQvIhJiSvIiIiGmJC8iEmJK8iIiIaYkLyISYkryIiIh9v8BY+T8w8UxpUwAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "dynamic X-site distortion: [0.0417 0.0782 0.0544 0.0739]\n", "dynamic B-site distortion: [0.093 0.0928 0.1074]\n", "dynamic tilting: [0.15 0.15 0.15]\n", " \n", " \n", "--Elapsed Time\n", "Data Reading: 00:00:00\n", "Structure Resolving: 00:00:01\n", "Tilting & Distortion: 00:00:07\n", "Total: 00:00:07\n", "------------------------------------------------------------\n" ] } ], "source": [ "from pdyna.core import Trajectory\n", "\n", "MDtup = (200,200,1,20000) # (T_initlal, T_final, timestep, nblock in VASP or equally recording frequency)\n", "\n", "new_sys = {'B-sites': ['Ti'], # The B-site elements in the structure\n", " 'X-sites': ['S'], # The X-site elements in the structure\n", " # The next two lines are the bonding information of the structure\n", " # characteristic value of bond length of your material for structure construction, doesn't have to be very accurate\n", " # the first interval should be large enough to cover all the first and second NN of B-X (B-B) pairs, \n", " # in the second list, the two elements are 0) approximate first NN distance of B-X (B-B) pairs, \n", " # and 1) approximate second NN distance of B-X (B-B) pairs\n", " # These values can be obtained by intuition or inspecting the initial configuration or, e.g. in the pair distrition function of the structure\n", " # Alternatively, you can use the 'get_bonding_info' function in the 'structure' module to get the bonding information\n", " 'fpg_val_BB': [[1.5,6.3], [2.9,5.8]], # For the B-B pairs\n", " 'fpg_val_BX': [[0.1,5], [2.4,4.8]], # For the B-X pairs\n", " } \n", "\n", "traj = Trajectory(\"extxyz\",(\"trajetory_1D.extxyz\", \"init.extxyz\", MDtup)) # The inputs for EXTXYZ format are (trajectory file, initial frame file, MDtup)\n", "# The initial frame file is optional, if not provided, the first frame in the trajectory file will be used as the initial frame\n", "\n", "traj.dynamics(read_mode= 1, # equilibration mode \n", " uniname=\"test_1D\", # A unique user-defined name for this trajectory, will be used in printing and figure saving\n", " allow_equil = 0.2, # use the entire trajectory\n", " read_every = 0, # read every step, 0 is automatic\n", " \n", " toggle_tilt_distort = True, # switch of octahedral tilting and distortion calculation\n", " \n", " structure_type = 3, # 1: 3C polytype, 2: other non-perovskite with orthogonal reference enabled, 3: other non-perovskite with initial config as reference \n", " multi_thread = 10, # enable multi-threading in this calculation, since not vectorized\n", " \n", " system_overwrite = new_sys,\n", " )" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "------------------------------------------------------------\n", "Loading Trajectory files...\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEWCAYAAACXGLsWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAdxklEQVR4nO3dfZRddX3v8ffHICGIgdAMGJLoRIwgcCnqmOLDVRQtUbmE9hYNFYlIm1XE5/qQ1F6xa91UWr0+oILNVSQohRsRJJiiYERoFYiDghCQJhqEgUgGEQHF0ITP/WP/oofhzOxJmHPODOfzWuuss89v//be3zPJnM/sh/Pbsk1ERMRIntLpAiIiYvxLWERERK2ERURE1EpYRERErYRFRETUSlhERESthEXEGJJ0hKSBhtfrJB3RuYoixkbCIp7UJN0u6WFJD0n6laTVkmaP0P8cSY+U/g9Kul7SK3Z2+7YPtv3dmhp7JVnSLju7nYhWS1hEN/gftvcAZgD3AJ+p6f/Ppf+ewFnARZImtbjGiHEtYRFdw/bvgAuBg0bZ/1HgX4G9gX2b9ZE0peyN/ErSLcCLhsy/XdKry/Q8Sf2SHpB0j6RPlG5Xl+f7yx7NiyXtL+k7kn4p6V5J50naa8h63yfpx5J+Len/SdqtYf4CSTeUbf1U0vzSvqekL0raJOkuSf87QRijkbCIriFpd+CNwLWj7D8JOBHYSLVH0sxpwP7lcRSwaIRVfhr4tO2ppf/K0v7y8ryX7T1sXwMI+CiwH/A8YDbwkSHrewMwH5gDHAq8pdQ9DzgXeD+wV1n/7WWZFcBW4DnA84E/Bf5qhJojAMgx0ugGX5e0FdgD2Ez1oT6S90l6O7D9L/WTbW8bpu8bgLfZvg+4T9IZwIeH6ftfwHMkTbd9LyOElu0NwIbycrDshZw2pNsZtu8GkHQpcNj2eoGzbV9RXt9V+uwLvJYqlB4GfiPpk8Bi4F+GqyUCsmcR3eFY23sBk4G3A1dJeoakZ5bDPg9Jeqih/8dL/ylAH/AxSa8dZt37AXc2vP75CHWcDDwX+ImkH0g6eriOkvaRdEE5VPQA8BVg+pBuv2iY/i1VGEK1F/LTJqt9FvBUYJOk+yXdTxUS+4xQcwSQsIguYnub7YuAbcDLbN9RDvvsUU5oD+1v2zcD3wNeP8xqN1F9OG/3zBG2v9728VQfzv8EXCjpaUCzoZ8/WtoPLYetTqA6NDUad1Id5mrWvgWYbnuv8phq++BRrje6WMIiuoYqC4BpwK2jXOZA4GXAumG6rASWSpomaRbwjhHWdYKknnLi/P7SvA0YBB4Fnt3Q/enAQ1QnvWdSnX8YrS8CJ0k6UtJTJM2UdKDtTcDlwP+RNLXM2/+JXBoc3SNhEd3g0nKY6QFgGbDI9nAf/gAfKIemfkP14folhj+m/w9Uh542lr5fHmG984F1pZZPAwtt/872b0td3yuHhw4v630B8GtgNXDRKN8rttcCJwGfLMtfRXUICqoT9rsCtwC/oro6bMZo1x3dS7n5UURE1MmeRURE1EpYRERErYRFRETUallYSDpb0mZJNw9pf4ek28ponP/c0L5U0oYy76iG9hdKuqnMO0PSaC8fjIiIMdLKb3CfA3yWatgBACS9ElhAde34Fkn7lPaDgIXAwVRfcvq2pOeWb82eRfUN02uBf6O6ouSyuo1Pnz7dvb29Y/l+IiKe9K6//vp7bfcMbW9ZWNi+WlLvkOZTgNNtbyl9Npf2BcAFpX2jpA3APEm3A1PLWDlIOhc4llGERW9vL/39/WPxViIiuoakpqMQtPucxXOB/y7pOklXSdo+QudMHjtkwkBpm1mmh7Y3JWlxGdWzf3BwcIxLj4joXu0Oi12ovj17ONU3UleWcxDNzkN4hPambC+33We7r6fncXtRERGxk9odFgPARWXMnbVUQxxML+2N4+vMAu4u7bOatEdERBu1Oyy+DrwKQNJzqYYduBdYBSyUNFnSHGAusLaMZfOgpMPLHsiJwCVtrjkiouu17AS3pPOBI4Dpqm5gfxpwNnB2uZz2Eaoxekw1Xs5KqvFqtgKnNtw/4BSqK6umUJ3Yrj25HRERY+tJOzZUX1+fczVURMSOkXS97b6h7fkGd0RE1EpYRERErYRFRETUSlhEtEjvktX0LlndNdvdGROp1m6XsIiIiFoJi4iIqJWwiIiIWgmLiIiolbCIiIhaCYuIiKiVsIiIiFoJi4iIqJWwiIiIWgmLiIiolbCIiIhaCYuIiKiVsIiIiFoJi4iIqNWysJB0tqTN5X7bQ+e9T5IlTW9oWyppg6TbJB3V0P5CSTeVeWdIUqtqjoiI5lq5Z3EOMH9oo6TZwGuAOxraDgIWAgeXZc6UNKnMPgtYDMwtj8etMyIiWqtlYWH7auC+JrM+CXwAcEPbAuAC21tsbwQ2APMkzQCm2r7GtoFzgWNbVXNERDTX1nMWko4B7rJ945BZM4E7G14PlLaZZXpo+3DrXyypX1L/4ODgGFUdERFtCwtJuwMfAj7cbHaTNo/Q3pTt5bb7bPf19PTsXKEREfE4u7RxW/sDc4AbyznqWcAPJc2j2mOY3dB3FnB3aZ/VpD0iItqobXsWtm+yvY/tXtu9VEHwAtu/AFYBCyVNljSH6kT2WtubgAclHV6ugjoRuKRdNUdERKWVl86eD1wDHCBpQNLJw/W1vQ5YCdwCfBM41fa2MvsU4AtUJ71/ClzWqpojIqK5lh2Gsn18zfzeIa+XAcua9OsHDhnT4iIiYofkG9wREVErYREREbUSFhERUSthERERtRIWERFRK2ERERG1EhYREVErYREREbUSFhERUSthERERtRIWERFRK2ERERG1EhYREVErYREREbUSFhERUSthERERtRIWERFRK2ERERG1WnkP7rMlbZZ0c0PbxyT9RNKPJV0saa+GeUslbZB0m6SjGtpfKOmmMu8MSWpVzRER0Vwr9yzOAeYPabsCOMT2ocB/AksBJB0ELAQOLsucKWlSWeYsYDEwtzyGrjMiIlqsZWFh+2rgviFtl9veWl5eC8wq0wuAC2xvsb0R2ADMkzQDmGr7GtsGzgWObVXNERHRXCfPWbwVuKxMzwTubJg3UNpmlumh7U1JWiypX1L/4ODgGJcbEdG9OhIWkj4EbAXO297UpJtHaG/K9nLbfbb7enp6nnihEREBwC7t3qCkRcDRwJHl0BJUewyzG7rNAu4u7bOatEdERBu1dc9C0nzgg8Axtn/bMGsVsFDSZElzqE5kr7W9CXhQ0uHlKqgTgUvaWXNERLRwz0LS+cARwHRJA8BpVFc/TQauKFfAXmv7b2yvk7QSuIXq8NSptreVVZ1CdWXVFKpzHJcRERFt1bKwsH18k+YvjtB/GbCsSXs/cMgYlhYRETso3+COiIhaCYuIiKiVsIiIiFoJi4iIqJWwiIiIWgmLiIiolbCIiIhaCYuIiKiVsIiIiFoJi4iIqJWwiIiIWgmLiIiolbCIiIhaCYuIiKiVsIiIiFoJi4iIqJWwiIiIWgmLiIio1bKwkHS2pM2Sbm5o21vSFZLWl+dpDfOWStog6TZJRzW0v1DSTWXeGSo3746IaKZ3yWp6l6wet+ubaNvfrpV7FucA84e0LQHW2J4LrCmvkXQQsBA4uCxzpqRJZZmzgMXA3PIYus6IiGixloWF7auB+4Y0LwBWlOkVwLEN7RfY3mJ7I7ABmCdpBjDV9jW2DZzbsExERLRJu89Z7Gt7E0B53qe0zwTubOg3UNpmlumh7U1JWiypX1L/4ODgmBYeEdHNxssJ7mbnITxCe1O2l9vus93X09MzZsVFRHS7dofFPeXQEuV5c2kfAGY39JsF3F3aZzVpj4iINmp3WKwCFpXpRcAlDe0LJU2WNIfqRPbacqjqQUmHl6ugTmxYJiIi2mSXVq1Y0vnAEcB0SQPAacDpwEpJJwN3AMcB2F4naSVwC7AVONX2trKqU6iurJoCXFYeERHRRi0LC9vHDzPryGH6LwOWNWnvBw4Zw9IiImIHjZcT3BERMY4lLCIiolbCIiIiaiUsIiKiVsIiIiJq7XBYSJom6dBWFBMREePTqMJC0nclTZW0N3Aj8CVJn2htaRERMV6Mds9iT9sPAH8OfMn2C4FXt66siIgYT0YbFruUsZzeAHyjhfVERMQ4NNqw+AfgW8AG2z+Q9GxgfevKioiI8WS0w31ssv37k9q2f5ZzFhER3WO0exafGWVbREQ8CY24ZyHpxcBLgB5J722YNRWY1HypiIh4sqk7DLUrsEfp9/SG9geAv2hVURERMb6MGBa2rwKuknSO7Z+3qaaIiBhnRnuCe7Kk5UBv4zK2X9WKoiIiYnwZbVh8Ffg88AVgW03fiIh4khltWGy1fVZLK4mIiHFrtJfOXirpbZJmSNp7+2NnNyrpPZLWSbpZ0vmSdivrvELS+vI8raH/UkkbJN0m6aid3W5EROyc0YbFIuD9wPeB68ujf2c2KGkm8E6gz/YhVJfgLgSWAGtszwXWlNdIOqjMPxiYD5wpKZftRkS00agOQ9me04LtTpH0X8DuwN3AUuCIMn8F8F3gg8AC4ALbW4CNkjYA84BrxrimiIgYxqjCQtKJzdptn7ujG7R9l6SPA3cADwOX275c0r62N5U+myTtUxaZCVzbsIqB0taszsXAYoBnPvOZO1paREQMY7QnuF/UML0bcCTwQ2CHw6Kci1gAzAHuB74q6YSRFmnS5mYdbS8HlgP09fU17RMRETtutIeh3tH4WtKewJd3cpuvBjbaHizruohqSJF7JM0oexUzgM2l/wAwu2H5WVSHrSIiok129h7cvwXm7uSydwCHS9pdkqj2Um4FVlGdSKc8X1KmVwELJU2WNKdsd+1ObjsiInbCaG+reqmkVeWxGriNP3yY7xDb1wEXUh3GuqnUsBw4HXiNpPXAa8prbK8DVgK3AN8ETrWdLwZG2/UuWd3pEmIc6pb/F6M9Z/HxhumtwM9tD+zsRm2fBpw2pHkL1V5Gs/7LgGU7u72IiHhiRrVnUQYU/AnVyLPTgEdaWVRERIwvoz0M9Qaq8wTHUd2H+zpJGaI8IqJLjPYw1IeAF9neDCCpB/g21bmHiIh4khvt1VBP2R4UxS93YNmIiJjgRrtn8U1J3wLOL6/fCPxba0qKiIjxpu4e3M8B9rX9fkl/DryM6hvV1wDntaG+iIgYB+oOJX0KeBDA9kW232v7PVR7FZ9qbWkRETFe1IVFr+0fD2203U91i9WIiOgCdWGx2wjzpoxlIRERMX7VhcUPJP310EZJJ1PdACkiIrpA3dVQ7wYulvQm/hAOfcCuwJ+1sK6IiBhHRgwL2/cAL5H0SuCQ0rza9ndaXllERIwbo72fxZXAlS2uJSIixql8CzsiImolLCIiolbCIiIiaiUsIiJabCzuptfpO/IlLCIiolZHwkLSXpIulPQTSbdKerGkvSVdIWl9eZ7W0H+ppA2SbpN0VCdqjojoZp3as/g08E3bBwJ/DNwKLAHW2J4LrCmvkXQQsBA4GJgPnClpUkeqjojoUm0PC0lTgZcDXwSw/Yjt+4EFwIrSbQVwbJleAFxge4vtjcAGYF47a46I6Had2LN4NjAIfEnSjyR9QdLTqO6bsQmgPO9T+s8E7mxYfqC0PY6kxZL6JfUPDg627h1ERHSZToTFLsALgLNsPx/4DeWQ0zDUpM3NOtpebrvPdl9PT88TrzQiIoDOhMUAMGD7uvL6QqrwuEfSDIDyvLmh/+yG5WcBd7ep1oiIoANhYfsXwJ2SDihNRwK3AKuARaVtEXBJmV4FLJQ0WdIcYC6wto0lR0R0vVENJNgC7wDOk7Qr8DPgJKrgWlnulXEHcByA7XWSVlIFylbgVNvbOlN2RER36khY2L6B6r4YQx05TP9lwLJW1hQREcPLN7gjIqJWwiIiImolLCIiolbCIiIiaiUsIiKiVsIiIiJqJSwiYsLp9I2AulHCIiIiaiUsIiKiVsIiIiJqJSwiIqJWwiIiImolLCIiolbCIiIiaiUsIiKiVsIiIiJqJSwiIqJWwiIiImp1LCwkTZL0I0nfKK/3lnSFpPXleVpD36WSNki6TdJRnao5IqJbdXLP4l3ArQ2vlwBrbM8F1pTXSDoIWAgcDMwHzpQ0qc21RkR0tY6EhaRZwOuBLzQ0LwBWlOkVwLEN7RfY3mJ7I7ABmNemUiMigs7tWXwK+ADwaEPbvrY3AZTnfUr7TODOhn4Dpe1xJC2W1C+pf3BwcMyLjojoVm0PC0lHA5ttXz/aRZq0uVlH28tt99nu6+np2ekaIyLisXbpwDZfChwj6XXAbsBUSV8B7pE0w/YmSTOAzaX/ADC7YflZwN1trTgiosu1fc/C9lLbs2z3Up24/o7tE4BVwKLSbRFwSZleBSyUNFnSHGAusLbNZUdEdLVO7FkM53RgpaSTgTuA4wBsr5O0ErgF2Aqcantb58qMiOg+HQ0L298FvlumfwkcOUy/ZcCythUWERGPkW9wR0RErYRFRETUSlhERESthEVERNRKWERERK2ERUTEE9S7ZDW9S1Z3uoyWSlhERESthEVERNRKWERERK2ERURE1EpYRERErYRFRETUSlhERESthEUX6oZrwqF73id013uNzkhYRERErYRFRETUSlhERESthEVERNRqe1hImi3pSkm3Slon6V2lfW9JV0haX56nNSyzVNIGSbdJOqrdNUdEdLtO7FlsBf7W9vOAw4FTJR0ELAHW2J4LrCmvKfMWAgcD84EzJU3qQN0REV2r7WFhe5PtH5bpB4FbgZnAAmBF6bYCOLZMLwAusL3F9kZgAzCvrUVHRHS5jp6zkNQLPB+4DtjX9iaoAgXYp3SbCdzZsNhAaWu2vsWS+iX1Dw4OtqzuiIhu07GwkLQH8DXg3bYfGKlrkzY362h7ue0+2309PT1jUWZERNChsJD0VKqgOM/2RaX5HkkzyvwZwObSPgDMblh8FnB3u2qNiIjOXA0l4IvArbY/0TBrFbCoTC8CLmloXyhpsqQ5wFxgbbvqjYgI2KUD23wp8GbgJkk3lLa/A04HVko6GbgDOA7A9jpJK4FbqK6kOtX2trZXHRHRxdoeFrb/g+bnIQCOHGaZZcCylhUVEREjyje4IyKiVsIiIiJqJSwiIqJWwiIiImp14mqoiIidkrsBdk72LCIiolbCIiIiaiUsIiKiVsIiIiJqJSwiIqJWwiIiImolLCIiolbCIiIiaiUsIiKiVsIiIiJqJSwiIqJWwiIiImolLCIiotaECQtJ8yXdJmmDpCWdricioptMiLCQNAn4HPBa4CDgeEkHdbaqiIjuMSHCApgHbLD9M9uPABcACzpcU0RE15DtTtdQS9JfAPNt/1V5/WbgT2y/fUi/xcDi8vIA4LYd3NR04N4nWG6rpcYnbrzXB6lxrKTGHfcs2z1DGyfKnfLUpO1xKWd7ObB8pzci9dvu29nl2yE1PnHjvT5IjWMlNY6diXIYagCY3fB6FnB3h2qJiOg6EyUsfgDMlTRH0q7AQmBVh2uKiOgaE+IwlO2tkt4OfAuYBJxte10LNrXTh7DaKDU+ceO9PkiNYyU1jpEJcYI7IiI6a6IchoqIiA5KWERERK2EBeN/KBFJsyVdKelWSeskvavTNQ1H0iRJP5L0jU7X0oykvSRdKOkn5ef54k7XNJSk95R/55slnS9pt3FQ09mSNku6uaFtb0lXSFpfnqeNwxo/Vv6tfyzpYkl7dbDEpjU2zHufJEua3ona6nR9WEyQoUS2An9r+3nA4cCp47DG7d4F3NrpIkbwaeCbtg8E/phxVqukmcA7gT7bh1Bd0LGws1UBcA4wf0jbEmCN7bnAmvK6k87h8TVeARxi+1DgP4Gl7S5qiHN4fI1Img28Brij3QWNVteHBRNgKBHbm2z/sEw/SPUBN7OzVT2epFnA64EvdLqWZiRNBV4OfBHA9iO27+9oUc3tAkyRtAuwO+PgO0W2rwbuG9K8AFhRplcAx7azpqGa1Wj7cttby8trqb6j1THD/BwBPgl8gCZfNh4vEhbVh+6dDa8HGIcfxNtJ6gWeD1zX4VKa+RTVf/hHO1zHcJ4NDAJfKofKviDpaZ0uqpHtu4CPU/2FuQn4te3LO1vVsPa1vQmqP2iAfTpcT523Apd1uoihJB0D3GX7xk7XMpKExSiHEhkPJO0BfA14t+0HOl1PI0lHA5ttX9/pWkawC/AC4Czbzwd+Q+cPnTxGOe6/AJgD7Ac8TdIJna1q4pP0IarDued1upZGknYHPgR8uNO11ElYTJChRCQ9lSoozrN9UafraeKlwDGSbqc6lPcqSV/pbEmPMwAM2N6+V3YhVXiMJ68GNtoetP1fwEXASzpc03DukTQDoDxv7nA9TUlaBBwNvMnj74tl+1P9YXBj+d2ZBfxQ0jM6WlUTCYsJMJSIJFEdZ7/V9ic6XU8ztpfanmW7l+pn+B3b4+ovYtu/AO6UdEBpOhK4pYMlNXMHcLik3cu/+5GMs5PwDVYBi8r0IuCSDtbSlKT5wAeBY2z/ttP1DGX7Jtv72O4tvzsDwAvK/9VxpevDopz82j6UyK3AyhYNJfJEvBR4M9Vf6zeUx+s6XdQE9Q7gPEk/Bg4D/rGz5TxW2eu5EPghcBPV72jHh4OQdD5wDXCApAFJJwOnA6+RtJ7qSp7Tx2GNnwWeDlxRfm8+Pw5rnBAy3EdERNTq+j2LiIiol7CIiIhaCYuIiKiVsIiIiFoJi4iIqJWwiAlD0rZy+eM6STdKeq+kp5R5fZLOGGHZXkl/2b5qH7f9KZKuKgNXbm97j6TfSdqzBds7rJWXV0v6dqdHmY32SljERPKw7cNsH0x1Xf/rgNMAbPfbfucIy/YCHQsLqnGJLrK9raHteKovhf5ZC7Z3GNXP53HKAIVP1JeBt43BemKCyPcsYsKQ9JDtPRpeP5vqw3Y68ArgfbaPlvQKqqHIoRrn6+VUQ1U/D9hINULqxVQfeNsHEny77e9LOgL4CHAvcAhwPXCCbUt6UVnv04AtVN+u/i3Vl9GOACYDn7P9L01q/z7wl7ZvL6/3By4FTgH+zvZRpf0twDFUo83uD1xs+wNl3slU30a+G1gPbLH9dknHUYXmNuDXVEOGbACmAHcBHy3vfT+q0LyXaqjus4EeqsEVT7J9h6RzgIeBA4FnASdRfTv7xcB1tt9SapkG/HsZRj26ge088pgQD+ChJm2/Aval+rD+Rmm7FHhpmd6DagDB388v7bsDu5XpuUB/mT6C6gN3FtWe9zXAy4BdgZ8BLyr9ppb1Lgb+vrRNBvqBOUNq3BX4xZC2vwf+V9nG7cA+pf0tZTt7ArsBP6cau2y/0m9v4KnAvwOfLcvcBMws03s1rOezDdv7CFXwTWn4GS0q028Fvl6mz6Ea20tUAxo+APy3Uuf1wGEN61wP/FGn/1/k0Z5HDkPFRNds1ODvAZ+Q9E6qD8+tTfo8Ffi/km4Cvkp146vt1toesP0ocAPVX+MHAJts/wDA9gNlvX8KnCjpBqph4/+IKnwaTQfuH9K2ELigbOMi4LiGeWts/9r276jGrnoW1X1XrrJ9n6sBBr865P2eI+mvqW6WNJxVth8u0y8G/rVMf5kqELe71LapQugeV+MXPQqsKz+L7TZThVh0gbE4dhnREeUw1DaqD63nbW+3fbqk1VTH7K+V9Oomi78HuIfqbnlPAX7XMG9Lw/Q2qt8T0XzoegHvsP2tEUp9mGovYXvdh1IFyhXVWIG/32v5XM32m7L9N5L+hOrGUzdIOmyYrr8ZocbG97Z9+48OqeVRHvuZsRvVe4sukD2LmJAk9QCfpzrU4iHz9i9/Df8T1WGhA4EHqQaU225Pqj2FR6kGaRzpL3KAnwD7lfMWSHp6OVH8LeCUMoQ8kp479IZKtn8FTNIf7qV9PPARl5FGbe8HzJT0rBG2vxZ4haRpZbv/c8j7vc72h6nOR8xu8n6H+j5/uF3rm4D/qHn/j1FGxH0G1aGx6ALZs4iJZEo53PNUqhvZfBloNmT7uyW9kuqv8luo7o72KLBV0o1Ux+XPBL5WTg5fych/dWP7EUlvBD4jaQrVX9SvprqFbC/VPQhEdbL42CaruJzqUM+3qT6kXztk/sWl/Z5htn+XpH+kOtR1d3lfvy6zPyZpLtXexxrgRqqhzpeUn9dHm6zyncDZkt5faj5ppPffxAuBa4c5xBdPQrkaKqINJD0feK/tNz+Bdexh+6GyZ3ExcLbti8esyB2r5dNU50DWdGL70X45DBXRBrZ/BFzZ+KW8nfCRsqdwM9UlwF8fg9J21s0Jiu6SPYuIiKiVPYuIiKiVsIiIiFoJi4iIqJWwiIiIWgmLiIio9f8BwQCfiXdjxEkAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEWCAYAAACXGLsWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAdlklEQVR4nO3de5QdZZ3u8e9jIgFEBEyDMYl29EQQOF4bROUoCgzxMoSZIxq8EJUxR8T7BYnOETyzMsOMLkdRwRMRAUVYEbkEGZQYEccjFxsEIUSGKBpaAmkGFRAMk/CcP+oNbJvdXZ1O7727s5/PWnt17bfeqvpVE/bTddlvyTYREREjeUKnC4iIiIkvYREREbUSFhERUSthERERtRIWERFRK2ERERG1EhYR40jSSZK+WaafIekBSVM6XVfE1kpYxDZN0m8kPVQ+tH8v6VJJs4fpu1Pp/+aGtidLWivpDVu6bdtrbe9ke1NNjW+X9JMtXX9EOyUsohv8te2dgBnA3cAXm3Wy/QCwCPiCpJ7S/C9Av+3z21JpxASVsIiuYfvPwPnA3iP0uRy4FDhF0kHAG4HjhusvaY6kKyXdL2kFML1hXq8kS5pa3r9d0q9L39slvUXSc4GvAC8tRz9/KH1fJ+nnku6TdIekk5qsd2E56rlH0icb5k+R9AlJvyrbum7z0ZSkvSStkHSvpFslvXGLf5HRlRIW0TUk7Qi8Cbi6puuHgIOoguWjtteN0PdbwHVUIfEPwMJhtv0k4BTgNbafDLwMuMH2auDdwFXllNUuZZE/AUcDuwCvA46VdMSQ1R4I7AkcDHyqBA/Ah4GjgNcCOwPvBB4sNawoNe9e+pwqaZ+Rfx0RMLXTBUS0wUWSNgI7AeuBw0bqbPv3klZRfaBfMFw/Sc8A9gMOsb0B+LGkS0ZY9SPAvpLWlgAaNoRs/6jh7S8knQu8Erioof3Tth8CbpR0I/B8YDXwd8Dxtm8t/W4s9b4J+I3tr5f26yV9B3gDsGqEuiNyZBFd4YjyF/s04L3AlZKe1nC30gOSHtjcWdJbgV7gB8A/j7DepwO/t/2nhrbfNutY+ryJ6ihiXbnQvtdwK5b0EklXSBqU9Mey3PQh3e5qmH6QKgwBZgO/arLaZwIvkfSHzS/gLcDTht3DiCJhEV3D9ibbFwCbgAMb7lbaqVwAR9LuwL8C7wL+F/BGSa8YZpXrgF3L6Z3NnjHC9r9v+1CqC+2/BL66eVaT7t8ClgOzbT+F6rqGRrmrdwDPHqb9Stu7NLx2sn3sKNcbXSxhEV1DlfnArlSna5r5EnCR7SvKqaLjga9Kmja0o+3fAv3ApyVtJ+lA4K+H2fYekg4vwbIBeIAqtKC6Q2uWpO0aFnkycK/tP0vaH3gzo3c68A+S5pZ9fp6kpwLfBZ4j6W2Snlhe+zVc64gYVsIiusEl5TTTfcASYKHtx52jLxeQDwQ+trnN9unAAPCpYdb9ZuAlwL3AicDZw/R7AvAR4M7S95XAe8q8H1JdM7hL0j2l7T3A/5F0f9n2stHsaPG50v9yqn3+GrCD7fuBvwIWlDruojrN9rggjBhKefhRRETUyZFFRETUSlhERESthEVERNRqWVhIOkPSekk3D2l/XxlmYJWkf2loXyxpTZl3WEP7iyXdVOadImm0tw9GRMQ4aeU3uM+kug3x0btDJL0KmA88z/aGck87kvamukNjH6ovOv1A0nPKaJ2nUQ3udjXwb8A84LK6jU+fPt29vb3juT8REdu866677h7bPUPbWxYWtn8sqXdI87HAyWVoBGyvL+3zgfNK++2S1gD7S/oNsLPtqwAknQ0cwSjCore3l/7+/vHYlYiIriGp6SgE7b5m8Rzgf0i6RtVInfuV9plU3y7dbKC0zSzTQ9ubkrRIUr+k/sHBwXEuPSKie7U7LKZSfXv2AKovPi0r1yCaXYfwCO1N2V5qu892X0/P446iIiJijNodFgPABa5cSzUK5/TS3vj0sllU3zAdKNND2yMioo3aHRYXAa8GkPQcYDvgHqoB0xZImiZpDjAXuLaMzXO/pAPKEcjRwMVtrjkiouu17AJ3GX//IGC6pAGqcXPOAM4ot9M+TDVGj4FVkpYBtwAbgeManlt8LNWdVTtQXdiuvbgdERHja5sdG6qvr8+5GyoiYstIus5239D2fIM7IiJqJSwiIqJWwiIiImolLCJGqfeESztdwqNGW0vvCZe2te7RbGu4msaj1rEu34rf0XDrHNpet9+t/H1tiYRFRETUSlhERESthEVERNRKWERERK2ERURE1EpYRERErYRFRETUSlhERESthEVERNRKWERERK2ERURE1EpYRERErYRFRETUSlhEREStloWFpDMkrS/P2x4676OSLGl6Q9tiSWsk3SrpsIb2F0u6qcw7RZJaVXNERDTXyiOLM4F5QxslzQYOBdY2tO0NLAD2KcucKmlKmX0asAiYW16PW2dERLRWy8LC9o+Be5vM+lfgeMANbfOB82xvsH07sAbYX9IMYGfbV9k2cDZwRKtqjoiI5tp6zULS4cDvbN84ZNZM4I6G9wOlbWaZHto+3PoXSeqX1D84ODhOVUdERNvCQtKOwCeBTzWb3aTNI7Q3ZXup7T7bfT09PWMrNCIiHmdqG7f1bGAOcGO5Rj0LuF7S/lRHDLMb+s4C7izts5q0R0REG7XtyML2TbZ3t91ru5cqCF5k+y5gObBA0jRJc6guZF9rex1wv6QDyl1QRwMXt6vmiIiotPLW2XOBq4A9JQ1IOma4vrZXAcuAW4DvAcfZ3lRmHwucTnXR+1fAZa2qOSIimmvZaSjbR9XM7x3yfgmwpEm/fmDfcS0uIiK2SL7BHRERtRIWERFRK2ERERG1EhYREVErYREREbUSFhERUSthERERtRIWERFRK2ERERG1EhYREVErYREREbUSFhERUSthERERtRIWERFRK2ERERG1EhYREVErYREREbUSFhERUauVz+A+Q9J6STc3tH1G0i8l/ULShZJ2aZi3WNIaSbdKOqyh/cWSbirzTpGkVtUcERHNtfLI4kxg3pC2FcC+tp8H/AewGEDS3sACYJ+yzKmSppRlTgMWAXPLa+g6IyKixVoWFrZ/DNw7pO1y2xvL26uBWWV6PnCe7Q22bwfWAPtLmgHsbPsq2wbOBo5oVc0REdFcJ69ZvBO4rEzPBO5omDdQ2maW6aHtTUlaJKlfUv/g4OA4lxsR0b06EhaSPglsBM7Z3NSkm0dob8r2Utt9tvt6enq2vtCIiABgars3KGkh8Hrg4HJqCaojhtkN3WYBd5b2WU3aIyKijdp6ZCFpHvBx4HDbDzbMWg4skDRN0hyqC9nX2l4H3C/pgHIX1NHAxe2sOSIiWnhkIelc4CBguqQB4ESqu5+mASvKHbBX23637VWSlgG3UJ2eOs72prKqY6nurNqB6hrHZURERFu1LCxsH9Wk+Wsj9F8CLGnS3g/sO46lRUTEFso3uCMiolbCIiIiaiUsIiKiVsIiIiJqJSwiIqJWwiIiImolLCIiolbCIiIiaiUsIiKiVsIiIiJqJSwiIqJWwiIiImolLCIiolbCIiIiaiUsIiKiVsIiIiJqJSwiIqJWwiIiImq1LCwknSFpvaSbG9p2k7RC0m3l564N8xZLWiPpVkmHNbS/WNJNZd4pKg/vjoiI9mnlkcWZwLwhbScAK23PBVaW90jaG1gA7FOWOVXSlLLMacAiYG55DV1nRES0WMvCwvaPgXuHNM8HzirTZwFHNLSfZ3uD7duBNcD+kmYAO9u+yraBsxuWiYiINmn3NYs9bK8DKD93L+0zgTsa+g2Utpllemh7U5IWSeqX1D84ODiuhUdEdLOJcoG72XUIj9DelO2ltvts9/X09IxbcRER3a7dYXF3ObVE+bm+tA8Asxv6zQLuLO2zmrRHREQbtTsslgMLy/RC4OKG9gWSpkmaQ3Uh+9pyqup+SQeUu6COblgmIiLaZGqrVizpXOAgYLqkAeBE4GRgmaRjgLXAkQC2V0laBtwCbASOs72prOpYqjurdgAuK6+IiGijloWF7aOGmXXwMP2XAEuatPcD+45jaRERsYUmygXuiIiYwBIWERFRK2ERERG1EhYREVErYREREbW2OCwk7Srpea0oJiIiJqZRhYWkH0naWdJuwI3A1yV9rrWlRUTERDHaI4un2L4P+Fvg67ZfDBzSurIiImIiGW1YTC1jOb0R+G4L64mIiAlotGHxaeD7wBrbP5P0LOC21pUVERETyWiH+1hn+9GL2rZ/nWsWERHdY7RHFl8cZVtERGyDRjyykPRS4GVAj6QPN8zaGZjSfKmIiNjW1J2G2g7YqfR7ckP7fcAbWlVURERMLCOGhe0rgSslnWn7t22qKSIiJpjRXuCeJmkp0Nu4jO1Xt6KoiIiYWEYbFt8GvgKcDmyq6RsREduY0YbFRtuntbSSiIiYsEZ76+wlkt4jaYak3Ta/xrpRSR+StErSzZLOlbR9WecKSbeVn7s29F8saY2kWyUdNtbtRkTE2Iw2LBYCHwN+ClxXXv1j2aCkmcD7gT7b+1LdgrsAOAFYaXsusLK8R9LeZf4+wDzgVEm5bTcioo1GdRrK9pwWbHcHSf8F7AjcCSwGDirzzwJ+BHwcmA+cZ3sDcLukNcD+wFXjXFNERAxjVGEh6ehm7bbP3tIN2v6dpM8Ca4GHgMttXy5pD9vrSp91knYvi8wErm5YxUBpa1bnImARwDOe8YwtLS0iIoYx2gvc+zVMbw8cDFwPbHFYlGsR84E5wB+Ab0t660iLNGlzs462lwJLAfr6+pr2iYiILTfa01Dva3wv6SnAN8a4zUOA220PlnVdQDWkyN2SZpSjihnA+tJ/AJjdsPwsqtNWERHRJmN9BveDwNwxLrsWOEDSjpJEdZSyGlhOdSGd8vPiMr0cWCBpmqQ5ZbvXjnHbERExBqN9rOolkpaX16XArTz2Yb5FbF8DnE91GuumUsNS4GTgUEm3AYeW99heBSwDbgG+BxxnO18MHGe9J1za6RLaolv2c7R6T7i0a34nY9nXkfp30+8ORn/N4rMN0xuB39oeGOtGbZ8InDikeQPVUUaz/kuAJWPdXkREbJ1RHVmUAQV/STXy7K7Aw60sKiIiJpbRnoZ6I9V1giOpnsN9jaQMUR4R0SVGexrqk8B+ttcDSOoBfkB17SEiIrZxo70b6gmbg6L4zy1YNiIiJrnRHll8T9L3gXPL+zcB/9aakiIiYqKpewb3fwP2sP0xSX8LHEj1jeqrgHPaUF9EREwAdaeSPg/cD2D7Atsftv0hqqOKz7e2tIiImCjqwqLX9i+GNtrup3rEakREdIG6sNh+hHk7jGchERExcdWFxc8kvWtoo6RjqB6AFBERXaDubqgPAhdKeguPhUMfsB3wNy2sKyIiJpARw8L23cDLJL0K2Lc0X2r7hy2vLCIiJozRPs/iCuCKFtcSERETVL6FHRERtRIWERFRK2ERERG1EhYREVErYREREbU6EhaSdpF0vqRfSlot6aWSdpO0QtJt5eeuDf0XS1oj6VZJh3Wi5oiIbtapI4svAN+zvRfwfGA1cAKw0vZcYGV5j6S9gQXAPsA84FRJUzpSdUREl2p7WEjaGXgF8DUA2w/b/gMwHzirdDsLOKJMzwfOs73B9u3AGmD/dtYcEdHtOnFk8SxgEPi6pJ9LOl3Sk6iem7EOoPzcvfSfCdzRsPxAaXscSYsk9UvqHxwcbN0eRER0mU6ExVTgRcBptl8I/IlyymkYatLmZh1tL7XdZ7uvp6dn6yuNiAigM2ExAAzYvqa8P58qPO6WNAOg/Fzf0H92w/KzgDvbVGtERNCBsLB9F3CHpD1L08HALcByYGFpWwhcXKaXAwskTZM0B5gLXNvGkiMiut6oBhJsgfcB50jaDvg18A6q4FpWnpWxFjgSwPYqScuoAmUjcJztTZ0pOyKiO3UkLGzfQPVcjKEOHqb/EmBJK2uKiIjh5RvcERFRK2ERERG1EhYREVErYREREbUSFhERUSthERERtRIWERFRK2ERERG1EhYREVErYREREbUSFhERUSthERERtRIWERFRK2ERERG1EhYREVErYREREbUSFhERUSthERERtToWFpKmSPq5pO+W97tJWiHptvJz14a+iyWtkXSrpMM6VXNERLfq5JHFB4DVDe9PAFbangusLO+RtDewANgHmAecKmlKm2uNiOhqHQkLSbOA1wGnNzTPB84q02cBRzS0n2d7g+3bgTXA/m0qNSIi6NyRxeeB44FHGtr2sL0OoPzcvbTPBO5o6DdQ2h5H0iJJ/ZL6BwcHx73oiIhu1fawkPR6YL3t60a7SJM2N+toe6ntPtt9PT09Y64xIiL+0tQObPPlwOGSXgtsD+ws6ZvA3ZJm2F4naQawvvQfAGY3LD8LuLOtFUdEdLm2H1nYXmx7lu1eqgvXP7T9VmA5sLB0WwhcXKaXAwskTZM0B5gLXNvmsiMiulonjiyGczKwTNIxwFrgSADbqyQtA24BNgLH2d7UuTIjIrpPR8PC9o+AH5Xp/wQOHqbfEmBJ2wqLiIi/kG9wR0RErYRFRETUSlhERESthEVERNRKWERERK2ERURE1EpYRERErYRFRETUSlhERESthEVERNRKWERERK2ERURE1EpYRERErYRFRETUSlhERESthEVERNRKWERERK2ERURE1Gp7WEiaLekKSaslrZL0gdK+m6QVkm4rP3dtWGaxpDWSbpV0WLtrjojodp04stgIfMT2c4EDgOMk7Q2cAKy0PRdYWd5T5i0A9gHmAadKmtKBuiMiulbbw8L2OtvXl+n7gdXATGA+cFbpdhZwRJmeD5xne4Pt24E1wP5tLToiost19JqFpF7ghcA1wB6210EVKMDupdtM4I6GxQZKW7P1LZLUL6l/cHCwZXVHRHSbjoWFpJ2A7wAftH3fSF2btLlZR9tLbffZ7uvp6RmPMiMigg6FhaQnUgXFObYvKM13S5pR5s8A1pf2AWB2w+KzgDvbVWtERHTmbigBXwNW2/5cw6zlwMIyvRC4uKF9gaRpkuYAc4Fr21VvRETA1A5s8+XA24CbJN1Q2j4BnAwsk3QMsBY4EsD2KknLgFuo7qQ6zvamtlcdEdHF2h4Wtn9C8+sQAAcPs8wSYEnLioqIiBHlG9wREVErYREREbUSFhERUSthERERtRIWERFRK2ERERG1EhYREVErYREREbUSFhERUSthERERtRIWERFRK2ERERG1EhYREVErYREREbUSFhERUSthERERtRIWERFRK2ERERG1Jk1YSJon6VZJaySd0Ol6IiK6yaQIC0lTgC8DrwH2Bo6StHdnq4qI6B6TIiyA/YE1tn9t+2HgPGB+h2uKiOgast3pGmpJegMwz/bflfdvA15i+71D+i0CFpW3ewK3jnGT04F7xrjsRLOt7Mu2sh+QfZmotpV92dr9eKbtnqGNU7dihe2kJm2PSznbS4GlW70xqd9239auZyLYVvZlW9kPyL5MVNvKvrRqPybLaagBYHbD+1nAnR2qJSKi60yWsPgZMFfSHEnbAQuA5R2uKSKia0yK01C2N0p6L/B9YApwhu1VLdzkVp/KmkC2lX3ZVvYDsi8T1bayLy3Zj0lxgTsiIjprspyGioiIDkpYRERErYRFg21lSBFJsyVdIWm1pFWSPtDpmraWpCmSfi7pu52uZWtI2kXS+ZJ+Wf77vLTTNY2FpA+Vf1s3SzpX0vadrmm0JJ0hab2kmxvadpO0QtJt5eeunaxxtIbZl8+Uf1+/kHShpF3GY1sJi2IbG1JkI/AR288FDgCOm8T7stkHgNWdLmIcfAH4nu29gOczCfdJ0kzg/UCf7X2pbjpZ0NmqtsiZwLwhbScAK23PBVaW95PBmTx+X1YA+9p+HvAfwOLx2FDC4jHbzJAittfZvr5M30/1gTSzs1WNnaRZwOuA0ztdy9aQtDPwCuBrALYftv2HjhY1dlOBHSRNBXZkEn3vyfaPgXuHNM8HzirTZwFHtLOmsWq2L7Yvt72xvL2a6ntpWy1h8ZiZwB0N7weYxB+wm0nqBV4IXNPhUrbG54HjgUc6XMfWehYwCHy9nFI7XdKTOl3UlrL9O+CzwFpgHfBH25d3tqqttoftdVD9sQXs3uF6xss7gcvGY0UJi8eMakiRyUTSTsB3gA/avq/T9YyFpNcD621f1+laxsFU4EXAabZfCPyJyXO641HlfP58YA7wdOBJkt7a2apiKEmfpDolfc54rC9h8ZhtakgRSU+kCopzbF/Q6Xq2wsuBwyX9hurU4KslfbOzJY3ZADBge/NR3vlU4THZHALcbnvQ9n8BFwAv63BNW+tuSTMAys/1Ha5nq0haCLweeIvH6ct0CYvHbDNDikgS1Xnx1bY/1+l6tobtxbZn2e6l+m/yQ9uT8q9Y23cBd0jaszQdDNzSwZLGai1wgKQdy7+1g5mEF+qHWA4sLNMLgYs7WMtWkTQP+DhwuO0Hx2u9CYuiXBDaPKTIamBZi4cUaaWXA2+j+iv8hvJ6baeLCgDeB5wj6RfAC4B/7Gw5W64cGZ0PXA/cRPU5MmmGypB0LnAVsKekAUnHACcDh0q6DTi0vJ/whtmXLwFPBlaU//e/Mi7bynAfERFRJ0cWERFRK2ERERG1EhYREVErYREREbUSFhERUSthEZOGpE3lVsBVkm6U9GFJTyjz+iSdMsKyvZLe3L5qH7f9HSRdWQas3Nz2IUl/lvSUFmzvBa28XVrSDybLyKwxPhIWMZk8ZPsFtvehuhf+tcCJALb7bb9/hGV7gY6FBdUYPRfY3tTQdhTVl0H/pgXbewHV7+dxyuB/W+sbwHvGYT0xSeR7FjFpSHrA9k4N759F9WE7HXgl8FHbr5f0SqqhwKEa3+sVVMM2Pxe4nWpU0QupPvA2D+T3Xts/lXQQcBJwD7AvcB3wVtuWtF9Z75OADVTfXH6Q6gtcBwHTgC/b/r9Nav8p8Gbbvynvnw1cAhwLfML2YaX97cDhVCO5Phu40PbxZd4xVN/MvRO4Ddhg+72SjqQKzU3AH6mG41gD7AD8Dvinsu9PpwrNe6iGrT4D6KEa3PAdttdKOhN4CNgLeCbwDqpvNL8UuMb220stuwL/XoYoj25gO6+8JsULeKBJ2++BPag+rL9b2i4BXl6md6IawO/R+aV9R2D7Mj0X6C/TB1F94M6iOvK+CjgQ2A74NbBf6bdzWe8i4O9L2zSgH5gzpMbtgLuGtP098L/LNn4D7F7a31628xRge+C3VGOWPb302w14IvDvwJfKMjcBM8v0Lg3r+VLD9k6iCr4dGn5HC8v0O4GLyvSZVGNwiWqwwPuA/17qvA54QcM6bwOe2ul/F3m155XTUDHZNRst+P8Bn5P0fqoPz41N+jwR+Kqkm4BvUz3warNrbQ/YfgS4geqv8T2BdbZ/BmD7vrLevwKOlnQD1TDwT6UKn0bTgT8MaVsAnFe2cQFwZMO8lbb/aPvPVGNHPZPqeStX2r7X1eB93x6yv2dKehfVg4iGs9z2Q2X6pcC3yvQ3qAJxs0tsmyqE7rZ9U6lzVfldbLaeKsSiC4zHucuIjiinoTZRfWg9d3O77ZMlXUp1zv5qSYc0WfxDwN1UT6t7AvDnhnkbGqY3Uf1/IpoPWS/gfba/P0KpD1EdJWyu+3lUgbKiGofv0aOWL9dsvynb75b0EqoHRN0g6QXDdP3TCDU27tvm7T8ypJZH+MvPjO2p9i26QI4sYlKS1AN8hepUi4fMe3b5a/ifqU4L7QXcTzW42mZPoTpSeIRq0MWR/iIH+CXw9HLdAklPLheKvw8cW4aER9Jzhj7QyPbvgSkNz6k+CjjJdm95PR2YKemZI2z/WuCVknYt2/2fQ/b3GtuforoeMbvJ/g71Ux57FOpbgJ/U7P9fKKPNPo3q1Fh0gRxZxGSyQznd80Sqh7p8A2g2BPsHJb2K6q/yW6ieFPYIsFHSjVTn5U8FvlMuDl/ByH91Y/thSW8CvihpB6q/qA+hetRrL3B9+QAdpPkjOS+nOtXzA6oP6dcMmX9hab97mO3/TtI/Up3qurPs1x/L7M9Imkt19LESuJFqGPETyu/rn5qs8v3AGZI+Vmp+x0j738SLgauHOcUX26DcDRXRBpJeCHzY9tu2Yh072X6gHFlcCJxh+8JxK3LLavkC1TWQlZ3YfrRfTkNFtIHtnwNXNH4pbwxOKkcKN1PdAnzROJQ2VjcnKLpLjiwiIqJWjiwiIqJWwiIiImolLCIiolbCIiIiaiUsIiKi1v8HOj2JnGHHzvkAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "traj = Trajectory(\"extxyz\",(\"trajetory_1D.extxyz\", \"init.extxyz\", MDtup))\n", "traj.system_test(B_sites=['Ti'],X_sites=['S']) \n", "# this is a pristine structure where the distribution of B-X and B-B pair distances appeared as sharp peaks, \n", "# the first two peaks in B-B plot are at 2.9 and 5.8, and those of B-X pairs are 2.4 and 4.8. \n", "# we highly encourage the users to use pristine inital structure as this will give the most stable fitting, espacially for complicated structrues. " ] } ], "metadata": { "kernelspec": { "display_name": "base", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.5" } }, "nbformat": 4, "nbformat_minor": 2 }