{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Derived fields and values\n", "\n", "In this tutorial, we show how derived fields and values can be computed afetr the micromagnetic system is defined.\n", "\n", "## Simulation\n", "\n", "First of all, as usual, we import `oommfc` and `discretisedfield`." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import oommfc as oc\n", "import discretisedfield as df" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We define the cube mesh with edge length $10 \\,\\text{nm}$ and cell discretisation edge $1 \\,\\text{nm}$." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAO4AAADuCAYAAAA+7jsiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsnXl8G4Wd9r8zOmzJh2z5iK/YceLciXPZSYCQQhOOUJYeCyyUFmhT3tKTY7nK+3KlpSeFLQVaWihsty20C0uhbRJCQqEUWELI4SuJ7Th2fMS3fMk6Z+b9Q8xYkiVZshWitHo+n35M4/FoJM0zv/v5CYqikEQSSZxZEE/3BSSRRBKxI0ncJJI4A5EkbhJJnIFIEjeJJM5AJImbRBJnIJLETSKJMxBJ4iaRxBmIJHGTSOIMRJK4SSRxBkIf4/HJNqskkji1EKI5KGlxk0jiDESSuEkkcQYiSdwkkjgDkSRuEkmcgUgSN4kkzkAkiZtEEmcgksRNIokzEEniJpHEGYgkcZNI4gxEkrhJJHEGIkncJJI4A5EkbhJJnIFIEjeBoCgKkiTh8XhIyuYmEQmxTgclcYqgKAperxen04nX60UQBPR6PQaDAZ1Oh06nQxCiGhxJ4p8ASeImAGRZxu12oygKHo8HQRAQRRFJkvB6vQB4vV48Hg9WqxW9Xo8oikki/xMjSdzTCEVRcDqdyLKMKIp0dXXR2tqKIAjodDqysrLIzs4mMzMTl8tFZ2cnJpMJRVECLHKSyP98EGKMpZKBV5ygKAput5uGhgZyc3M5efIkoiiyYMECwGdhbTYbQ0NDDA8Pa38zf/58LBYLgiCgKIoWC6tEVv+XJPIZi6i+tCRxTwP8XePa2lqGh4epqKigqKhIc4+DSTc4OEhraytpaWmMjIyg1+sDLHKSyP8wiOpLSrrKHyHUBJQat7a1tTE4OMj8+fMpKirSjgEYGoKjRwWqqxVEEYxGI6mpqSxcuBAAl8vF0NAQ3d3dNDY2YjQaNSJnZGQgCIIWF0OSyP9oSBL3I4LqGqvWtq6ujvT0dIqKikhJSQk4dnwczj/fzLFjOhYskLj0Ui9VVUYKCyeqdykpKcyaNYtZs2YBPiLbbDa6uroYHR0lJSVFI3J6ejqCIODxeBgYGGB0dJSSkpIkkc9gJF3lUwz/2iz4XN6jR4+yYMEC8vLyaGpqwmKxkJ+fD0BNjcR116Vw+LAICKxc6aWuTofXK6DXy1RVyZx7rsTGjRJr10qYTKFf1+l0YrPZsNlsjI2NkZKSQnZ2NqIo4nA4qKioQJZl7fjg8lOSyKcNyRj3dMPfNVYUhebmZkZGRli+fDmpqakAGnHz8vJ59lmRW2/Vk56ucOWVHh5/PIX6+jGysxXeeMPDjh0OGhoKOHBARJYFUlIU1q6VNCJXVUkYjaGvxeFwYLPZ6OnpYWxsDIvFQnZ2NllZWaSnpwMkiZwYSBL3dEKWZTweD7Is43Q6qa2tJS8vj/Ly8gACNDc3A5ls21bEH/6g47zzJJ580s6ePQa++lUTtbVjlJUpjI+Pc/z4cZYuXcrwMLzzjo633tLz1ls6ampEFEXAZFJYv95H4nPP9bJ6tYw+KBgaHBxkYGCA4uJihoaGNItsNps1IqelpWnvQYUgCFrpKUnkU4pkcup0QFEUHA4HHo8Ho9FIb28vx44dY+nSpWRlZU06/vBhE7ffXsCJEyL33+/llltcKIovIQXgxx0tcWWxwJYtElu2SAAMDsI77+j52990/O1vOh54IAVIIT1d4ayzJDZu9LJxo0Rlpe9koihiNpsxm80UFRVp12yz2WhtbcVut5OWlqbFyGazGQC3243L5dII6/V6yczMTBL5NCBJ3DhC7Xzq6OhAURTGxsbwer2sXbsWg8EQdCw88YSOu+6aQ26uzGuveTjnHAWvF7xeQhI3HKxWuPRSL5de6stW9/cLvPWWj8RvvaXjnnt8bnlWlkJVlY6VK/V8+tMiS5fKiKLPmqpELi4uRlF8Ft5ms9HS0sL4+DhpaWlkZ2eTnZ2NyWTC6/VSV1fHihUrAN/DQLXIer0eQRCSRD6FSBI3TvCvzbrdbjo7O5k3bx4lJSUharLw5S/r+dOfdJx33hg//vEAS5cWBBwTTFy1ThsNcnMVPv1pL5/+tI/I3d2CRuI33jCwe3cxDz0EVqsv0aXGyAsXygiC77XS0tJIS0ujpKQERVGw2+3YbDaOHTuGw+HAZDLhdrtxu92YPsyQ+Vtk1bVWY+QkkeOLJHFniODabEdHB11dXRQWFjJ79uxJx7/7rsDnP2+gpwd+9CMvl17aidFomHScTuf7KcsCM00tFBQoXHmllyuv9DIwMMCRI3ba2uby1ls+9/rll32vn58vs2GDxMc+5ouR581TNCKnp6eTnp7O7NmzURSFoaEhjhw5QnNzM06nk/T0dC1GVhNvKrHBZ5H9k11JIs8MSeLOAP61Wa/XS319PUajkXnz5mnlHxWyDA89pOOBB3SUlsIbb3hYs0ahrS20JZ2JxZ0KBQUezjnHy2c/60VRoLVV0Ej81ls6/ud/fEQuKlJLT74YuaxsoisrLS0Nk8lEZWWlFhbYbDYaGxtxuVxkZGRoMbJap45EZFFMTpjGgiRxpwn/rqShoSEOHz7MvHnzKCgooLu7W7tBAXp6YOtWA7t3i1x+ucTjj3uxWHy/C0dI9T6WpInjTgUEAcrLFcrLPVx7rQdFgeZmgb/9zZex3rNHx+9/7yNyWZnqWntZv96jXZMgCGRkZJCRkUFpaSmyLGtEPnLkCG63m4yMDC1GNn5Ys1KJPDAwQHp6OhkZGUkiR4kkcWNEsGt8/Phx+vv7Wb16tRbr+ZPx9dcFvvhFA0ND8NhjHrZu9cWRwedUoZJhwlUOfdypgiDA/PkK8+d72LrVR+QjR0QtY/2Xv+j5zW8MgImSkjVccIFOi5NnzfJdnyiKZGZmkpmZSVlZGbIsMzo6is1mo6GhAY/HQ2ZmpkbkoaEhDAYDKSkpSYscJZLEjQH+tVmXy0VdXR0Wi4Xq6uqAG0sQBNxumfvv1/GDH+hYsEDhz3/2smzZZOKFt7jKh6956t5PNBAEWLxYZvFimS9/2YMsQ12dyOuvK7z6qpsXX8zmmWd8FnThwolE17nnSuTkTBDZYrFgsViYM2cOsiwzMjKitWiOjY3hcrmYNWsWWVlZGAwGLQxJEjk0ksSNAupN5PV6EUWRgYEBGhsbWbRoETk5OZOO7+7Wc+ON5ezfr+faayUeecTLhz0NkyAIQkCjw8S/+376x7iJAFGEykqZ+fMdbNrUzOLFyzl0SNRi5OeeM/DUUz4iL1s2QeSzz/aSna2eQyQrK0urax8+fJjMzEzsdjsdHR1IkhTQ2aXX60MSObgh5J8JSeJOAbU2e/jwYXJychgcHMRut1NdXa3Fav7Yvl1k69Z8XC741a88fPazkU1mtDGuei2JBr0e1qyRWbPGzc03g8cD+/eLWoz8zDMGfvYzI4KgsGKFrHV1nX22REbGxHkyMzPJyMigvLwcSZIYHh5maGiI9vZ2JEnSiO5PZJfLhcvlAv75iJwkbgT412YlSeLo0aPMnj2bhQsXTrKAbjfcc4+On/xEz9KlHr773WYuumjulK8RjrjBMW48s8rxgKrCEQyDAdatk1m3zs3tt4PLBfv2TTSD/PznBh591IhOp7BqlczGjV5KSzO49FJRI7JOp8NqtWK1WgE0IttsNk6cOIGiKAEWWafTaUSura1l3rx5Wg35H5XISeKGgH8CShAEuru76e3tpby8nDlz5kw6vqUFrr3WwL59Il/5isQddwxgszmjeq2pLK6iJIaLHIxoHyIpKXDOORLnnCPxrW+BwwHvvacSWc+jjxrxehdz++0KVVVSwOTTh+XgSUT2er0akdva2lAURSs9jY2NBRBZtcg6nU6Lj9XOrjMZSeIGwb82K0kShw8fRlEUSkpKtKyxP154QeSrX9UjivD88x4+9SmZoSGBwcHobmx/4qq14PHxcXp7S4G5eDzypONminjdtNM5j8kE550ncd55EuBmbAyef76dlpYy3nknlYceMvLDH05MPm3c6PvfmjUTk096vZ6cnBwtv+D1erWBCbvdzsGDBzUiWywWRFHUhj1UqERWLfKZRuQkcf3gr2k8OjpKfX09ZWVlFBUV0draGpBEcjjg9tv1PPWUjnXrZP7zPz2oxng6JBsdHaW2tpbS0lLKy8vp7PQlYQ4fPkpmpgOLxaKNB87kJosX+Wd6HSrS02HtWhuf+1weqalywOTT3/6m47vfNfLggwJms8K6dSqRvaxaNTH5pNfryc3NJTc3F5vNRmVlJUNDQwwMDNDS0oIgCJpb/Y9C5CRxmVybPXHiBCdPnqSyslKbVfUn45EjAp/7nJ66OpF//3cv998v4T9DIIpi1AQRBIHh4WF6enqorKzEZDLh8XjIyfGloRcuXMyCBWMMDAzgcrnYu3fvpIb/RL/JpoL/QyDU5NPbb090damTTxkZE5NP557rm3xS8wIGg4G8vDzy8vIA8Hg8DA0N0d/fz7FjxyYpaJ6JRP6nJ65/bdbj8VBfX4/JZGLdunUBCQ2fzrHMr38tcvPNetLS4JVX3Fx4YfS12WB4vV7a2tpwu92sX78evV6P9GEa2T+rbDKZKCwspKenh9WrV2sN/2qfsH9XUrAMzqlCvCwuoMnThoLVCv/yL17+5V98D9W+PoG//12nNYTs2jUx+XTOOV7KykowmUSWLJG1zzCYyG63m6GhIXp7e2lubg4pvCfLMvv376e8vByTyZRwRP6nJW6wpIzatjh//nxNRsYf4+M67r03l1deMbBxo8yzz3r4UN9tEsLVZv0xNjZGbW0tWVlZWnOBP4IbMPzbC/0b/oO7krxer5Zxzc7OnnTeeCGexI3lXHl5gZNPJ0/6RhjfekvHm2/q+ctfFvDEE5CT4xuYOPdc39DEggUTHWtGo5H8/Hzte3a73Zo6SFNTE3q9nuzsbFwul5aRDrbIf/jDH7jooosoKyuLy2cQK/4piasmoOrr61m0aBEtLS3YbDbWrFmjTbb449AhgauuKqGtzcA993i56y5Jc8tCYSqLqwqfL1++HJfLxcDAQMDfQuh53NDZ58CupODSCRCQqIknTgdxg1FYODH5BC7+/OcahodXa+og6uTTrFk+Iqt1ZHXyCXxEDiW8pyqXGI1G7UGoCu+99tprnHXWWfF4+zFDEITUfzri+tdmbTYbH3zwAVarlerq6kk3j6LAk0+K3HmnnqwsiWeeOcG//VtBmDNPIFyMK0kSR44cwePxsHbtWvR6vXYtk8+hXm9s7y+4dBIc3ymKgl6vZ3h4mIyMjGnXN+NZU46n9Z41y8Wll3q55hrf5NPx44GTTy++6CNycfHE5NO5505MPoFPQbOgoICOjg5Wr16tWeTOzk5GR0d5/fXX6e7upqWlhSVLlqCL9BQHvvjFL/LnP/+Z/Px86urqQr7/m266ie3bt3Ps2LEa4HpFUfZHOGXGPw1xg2uzfX19jI+Ps2TJEu0m98fQENx4o54//lHHxRdLPPjgSVJT7VG9ViiLa7fbqampobi4mNmzZwe4vqdyrC84vjt58iR9fX2ajGtqaqpmTdLS0qIm0Olylac6jz8EAebOVZg718N11wVOPv3tbzp279bx/PMTk0/5+TIpKfCXvzgQBD6UEBJJTU2lsLCQwsJCAHJycti9eze/+93veOqpp/jTn/4U8bquv/56vv71r3PttdeG/P2OHTtoamqiqakJURT/D/AzYF2EU56iACjB4F+blWWZxsZGnE4nZrM5JGnfe8837N7VBd//vpdvflOivx+Gh6Mzf8Exbnd3N8eOHWPZsmWT3NWpO6fimwQxGAykp6czd+7ckFpT6kC8mrEOh3h3ccWDuJGSXL7XmDz5dPiwb/LphRf07N3r+9BdLggRMWmYN28egiDw9NNPa1WHSNi4cSOtra1hf//yyy9z7bXXqvfC/wqCkCUIQqGiKCfD/IntH564/rVZu91OXV0dRUVFLF68mHfffTfgWFmG//gPHffeq6O4GF5/3cPatRMTLlMlnFSoZJRlmaNHj+JwOELqTkVCKIsbb4TSmgoeiPcfvwvuzT7dmdVgxGq5BQGWLJF57z0dH3ygIzMTcnKUiKRV4XQ6Iz7YYkFnZ2ewWkoHUAyEJK6iKM5/WOIG12a7uro4ceIEy5YtIzMzc9LxfX2+Yfddu0Q+8xmJJ57w4i/KGIurKooiXq+XvXv3UlBQwKJFi8LeULEMGZxqhBqIV8fvOjs7kWVZy1jH01WOF6ayuJOPhwceMPLIIylceKEXjweGhqIPFU5n//M/JHH9a7OSJNHQ0IAoilpCKBhvvilw/fUGBgfh0Uc93HDD5GH3WCxuX18fY2NjVFdXk63OsoVBtEMGpwP+43fq1I7aWtjf348kSYiiGNBaeDoRC3GdTvjKV1J58UUDW7e6+dGPXFx2mQmTaeqHc7zDhOLiYtrb2/3/qQTojPQ3/1DEDa7NjoyM0NDQwJw5c7SlWv6QJPj2t3V873s65s1TePllL5WVob+UaLqhZFmmqakpQGB8KpyOQfrpWkqdTqf1CGdnZ2uSM729vTQ1NWllE6vVqpVNPkpES9yBAfjsZ028+66ebdtc3HSTG0EAh0PAYvF97tGQM17v77LLLuOxxx7jqquuQhTF9cBwhPgW+AcirtpfrE5+tLW10dPTw8qVKzVBb390dsJtt63i0CE911wj8ZOfeImUZ5iqqcLpdHLo0CHy8vJYvXr1pPg50nlD3STBg/TxQjx7lfV6fUAjg7qvqL29ndHRUe3hpYqqn2oiR0PclhaByy83094u8OyzDj7zGa/2O6fTp4gZzblieS9XX301b7zxBv39/ZSUlPDAAw9oxuXGG2/kkksuYfv27VRUVAD8EvjCVOf8hyCuWps9cOAAy5cv58iRI6SlpbF27dqQH/6rr4ps3apnbEzPk0+6uO66qV8jkqvc39/P0aNHWbx4ccgsdSQkUow7U/iXTUKJqqutmVar9ZS0Zk4Vd+/dK3LVVSZkGf70Jwfr1wd+uA6HoCWmIhFXkqSYiPvcc89F/L0gCDz++OPq/10ezTnPaOL6u8aCICBJEgcOHGDhwoVa3dIfHg/cd5+Ohx/Ws3y5zJ13HuSyyxYAYTZl+SEUwdRFXkNDQ1RVVU37ZowU4ybQ7HwApiJJKFH1YME4/0RXPBCJbK+8oudLX0qlsFDhxRfHqaiY/JpOJ9r2Q0mSwjZWOByOkF7cR4kzlrj+tVlFUWhqasLlcrF69eqQO3paW33D7nv3ivyf/yPxgx94aWhwRZ1wCra4LpeLmpoasrKyqKqqmrYbGOrvBEGYdufUR4XYSy/CJOVHtTXT4XCwb98+rTVTVbWIFaGIqyjw+OMG/u//TaG6Wub55x3k5oZ+UPgs7tSu8vj4eJK404F/26L/Jrzs7OyQtdKXXhK58UY9igK/+52Hz3xmYvlVLMRVLYNvG8ARFi5cSG5u7ozey9SucmKVXOIFNRutJrlWrlzJ0NAQg4ODtLS0oNPpNLc62tbMYLJJEtx1VwpPPmnkk5/08ItfOMPuE4boLa7dbk8SNxYEty0Gb8KrqakJIKLTCXfeqefJJ3VUV8v8+tceyssnzhdrU4UkSRw7doyBgYGwAwmxIpHLQZEQ7zZF/2F4mJjYUVsz1cXcVqs1bGumP3Htdti6NZXt2w1885tutm1zEYn7ihK9xU26yjHAvzYryzJHjhyZtAlPp9Np86yNjb5h95oakZtv9rJt2+Slz7EQ1+v1MjIygsVioaqqKm41y2h7lRMN8SRuqPMET+yEWgOqElntYFLJ1tMjcOWVJg4dEvnxj53ccINn0vmD8aE0VdQxbry6pqaLhCducG12bGyMuro6Zs+ePWkTnkrE3/5W5Jvf1JOaCi+95GHLltB3f7TEtdls1NfXk5KSwsKFC+Pzxj6EP3HVB1J/fz+ynAOsxu0+g9LK00C0tVeTyYTJZNL2+apiAv67igRBoLMzk5tvNtPXJ/Dccw5NSWMqOBy+n9HGuGnhhLI/AgiCoEto4obahNfZ2cny5cvJ8Bfl/RAul4FvfjOTF14wsGGDb9i9pCT8+X2qFuG/WEVRaG1tpbe3l9WrV3Pw4MEZv6dgqMT1rwNXV1fT3u6bROrsPMm+fe1kZWXh9XojWoKPEqfa4kZCODGBl16ycffdczEaPTz5ZAvV1al4vVlRiQk4nb5rUA1pJOImY9wIkGWZ7u5usrOzAzbhrV27NuSNW1srcPXVFbS2Grn7bi933y1NOfsUyeJ6PB5qa2sxmUyTVozEE751JW4++OADFi1aRHZ2Nh6Ph6wsXz91UVEJK1fmMDQ0xMmTJ9m/fz96vV6buT0dHUpweokbDFEU2bEjh1tvLWbOHDcvvujEYjEzODhIa2srgiCQlZWF1WolMzMz5P0TbHGncpVPp8UF5IQjrr9rfOTIEZYsWUJDQ4O2CW/y8fDUUyK33aYnI0Piv/6rm8svj64JIlwb4/DwMHV1dVRUVGgx1qmAoii0t7fjcDjYsGEDqamp2vVM6CpPJG5SU1Oprq7WOpROnDjB2NgY6enpGpE/Ks2peGHmqpXwox8Z+c53Uli7downnuihvDwfmCwm0NfXp2lM+WesBUGIyeImk1NB8K/Ngi+z2NjYGLAJzx/Dw/DVr+p58UUdmzfLbNvWhq/7Lnri+ltcRVE4ceIEXV1drFq1asZfTqSb0uv1UldXh16vJz09fVKGOlLnVHCH0tjYGIODgxw+fFhrbLBardq6jmAkiqWE2Cd6/OHxwM03p/Bf/2Xkqqs83H57KxbL5PskWExAlabp6OhgbGyM1NRUOjuLgPIAixtqxQycfldZURQlYYirKs+r5K2trQVgzZo1k24+SYIPPhC49loD7e3wne94ufVWia4uIsaswfAnrkokg8EQ1h2PBWrsGurmVtUwSktLKSgo4P333590TLSD9P6jeGVlZZrmlOomqvXSnJwcMjIyEk5XebrnGR6Ga6818de/6rnzThd33+3m+HEpqoeAKk1TUFCgiQm0tPiE4NrajlBf7zMe4YZEHA7HjOv3M4EgCOsShrjql+e/Ca+lpSVkDHrFFXq2bxdJTYVvfUviuuskRNFHxOBN8JGglo9GRkaoq6ujvLxckyeJx/sJRZKenh6am5tZvnw5mZmZWudXMOKlOeV2uxkcHKSzs5ORkRF0Oh1GozEuJY3TRdyODoErrjBx9KjIz37m4JprvNM+lyomkJ7uS3auXLmQ0tJhmpqaOHnyJF1dXQGqmQaDIRFc5ZcShriqWoT/Jrzg7QEqfH2mAkajwne+o+fBB3WsXq1w9tnZVFfbmD07OlEeQRAYGBigo6MjQPw8HggmrqIoNDY2MjY2FlB7DtcDLYq+G3CmQwZGozHAurS3tweUUYJvymhxuixuTY3IFVeYsNsFXnzRwfnnT3xAM3G7/WNc1YPJzc3FYrEwMjLC4OAgHR0dvP3227z77rukpKRo+YVI2LlzJzfddBOSJPGlL32Ju+66K+D3J06c4LrrrmNoaAhJkqitrb1EUZTtU1zucMIQV10T4b8Jz7+hwh/nnqvw6KOwfbvPuu7aJfLaayKPP25BlrP4xjcUzj9f5oILfP8rLZ38el6vl87OThRFiYtrHAx/N9ztdnPo0CGys7NZvXp1wI0arlf5VHROCYJAampqwIJptV9YHeRWSfxRDcbHQtxdu3Rcf72JrCyFV18dZ+nSwA9nJsRVs8r+ySlVU9l/l+/cuXPZv38/9fX1XHzxxfz1r38N+8CTJImvfe1rvPbaa5SUlFBdXc1ll13GkiVLtGO+853vcOWVV/KVr3yFhoYGli5d+gQwZ4rL/WpCEVe1CirCEVdVKXC74eyzfVve7r5b4tixQbZvd9HQMJtdu0T++Eff3b9okcyFF/pIvGGDgtfr29NjsVhITU09JXVR1ZIODQ1RX1/PggULQk4shf97389T2Tnl3y8MvuyrzWbTBuNTUlI0tzt4nvajTk4984yBW29NYdkymT/8wUFh4eTwYmYW1/fTPzkV6lxZWVmYzWa+/vWvs3bt2ojn3Lt3LxUVFcyd61u3etVVV/Hyyy8HEFcQBEZGRgBfNQPomupaFUX5a8IQNxTCE9f3U31KqsjJEfn4xwf5xjcKURTfjh/VGj/5pI5HH9WTkiJTWSlw2WXVnH32KHl5A5POHwmx3LCdnZ309PRMK0MtCCAIykfa8mgwGAIG4x0Oh9b0r87TqkT+qFxlWYZt24w8/LBPF+rZZx1hBQ9mZnGjLwdF2zkVLAJXUlLCe++9F3DM/fffz4UXXshPf/pT7HY7wDeiud6EIm5wvBeOuCoHfB/2xPH+7qkgwOLFCosXS9x0k8ToqMRzz3Xy9tsZHDiQzz33iEA6hYU5bNmi48ILZc4/XyaS2H+kTLE/JEnS9rTOxA0XxdPbq2wymSguLtYUIEdHRxkYGKCurg6Hw0FGRgYGgwGLxTLt9xjp8/TXhfriF9089JArYu4i3hY33HuKZ8vjc889x/XXX8+///u/8+6773L22Wf/lyAIyxRFifjNJxRxgxGrxQ3XwuiTZa1hy5YSbrghC0Hw0NoKL7/s5NVXFf77vy386lc6dDqFdesUza1etUoJmChRHwyRbo7x8XFqamowGo0sWLBgRm64Tpc4Qwb+87Tl5eW0tLSgKIq2IWG63VyyLIc8NpwuVCTMRHkxFosbbVY5WASuo6OD4uLigGOefvppdu7cCaCuNEkFcoHeSOc+I4mrPhWDiavT6SZloU+ePMnx48cnybLOmQPXXedk8+ZOFixYynvv+dzq3btF7r9fz/33Q26uwqZNvvh482Z5SsG4vr4+GhsbWbp0qbYpfSY43RY3EtT6sRq3h+rmUruTIo0/hiJbJF2oSAj3EIgGTifodIq2LnUqixtNBaK6upqmpiaOHz9OcXExzz//PL/73e8CjiktLWXPnj1cf/31HD58GHzE7Zvq3AlP3FDloAlXefLxKtHVSRuXy0V1dXXIzJ9qQQ0G2LBBYcMGiW3bJHp7Yc8ekV27RPbsEfn9731f4IIFK/nEJwy2KJVPAAAgAElEQVRs2SKwfr2ijQkqisKxY8ew2WxaKevEiRNxIW6iDtIHv7dQ3Vw2m43Dhw9rGwRDdXMFu8qqLpQkCbzyioOzzoq+HjbTGNe/rB3JhXe5XFG1lur1eh577DEuuugiJEnii1/8IkuXLuXee++lqqqKyy67jB//+MfccMMNPPLII+rrXa9EceMkFHGDP6hwru+Eqzz5eFmWNXe1sLCQxYsXh/0Cwg0Z5OfD1VfLXH21jCz7tvXt2iXy0ksyP/2pkUceEUhP95WcNm3yUFxcz4IFBtasWaPdONPd9eOPRHKVQyGSyLu/sHq4bi6r1RpgJf/0Jz1bt6q6UPaQulCRMBPi+taOBO8eCv3+YnHJL7nkEi655JKAf9u2bZv230uWLOHtt9/2//WuaM6bUMQNhk6nw+12T/p3lbjj44H/LooiLpeLAwcOaKoYkRDNPK4owqpVCqtWSVx4YQOFhQvZuzedXbtEXn0V/vSnVGAN8+f74uILL1Q499yp3epocCpc5dPROBGum6urq4vBwUF0Oj2//W0+P/xhPlVVMr//fXhdqHhdUzCCLW6k10gEJDxxQ1lcvR4MBiXAVVaXebndbjZu3Bi2QdwfsShgqMenp8tcdplMVVUHV1/ditm8knfeyWDXLoFnntHxxBO+jq5Vqyq46CL45CcFlixRpkyshIIgTO6cSsTVH7HCv5urvb2Lb3/byvPP5/Gxj/Vzxx219PenI8vWmLu5YPptmE7nZIsb6TVO93eQUMQN/jDCERd8VldN4TscDmpqasjPz8dkMkVFWpgecdWVJm63m3XrfCtNVq2S+NrXfNfz978LvPaayJ//bGDbNhPbtkFxscLmzb4k18c/LhPFggPAlyzxv7zTfbP4Ix4PELsdvvGNIl5/PYNvfMPNt7+dAqzWWgw/ym6uaC0uJIbVTSjiBmMq4o6PC1omd8mSJWRnZ3PyZMTNDQGIlbiyLFNXV0dxcXHI2Dk1FTZvVti8WeILX2jE7Z7FBx/k8NprIi+/LPKf/6lDFBXWrp0g8po1gTeBJEn09PSQlZWFKKZN0lVOFIs70+vo7fXpQh08KPL//l83d9yh1kUDWwxj6eaaCXwW1/ffkYjp9XoTQoHkDCauQk/PCG1tbVomN1bEQtyBgQH6+/upqKigrKxsyuMFQaCoSKKyUuYLX5DxeuH9933W+LXXRB58UMd3vqPHalVYsWIJV18tcs45dnp7D5KZmUl3dzeStBabbYSBAfuU8fqZhKNHRS6/3ERfn8BPf9rBpk3jQOiGhli6uWYCh2OilXaqGu5pVr8AzgDihiKW0+kEvLjdOtasWTMjMfKp3B5/3an8/PyQWlehEPxQ0OvhrLMUzjpL4t57fYuyX3/dV3Lavj2Lv/7VAGSxZMk5bN4sc8EFEkajAb0+lcHBdlpaWnA4HLS3t5Obm/uR7OKJhOla3Lfe0nHNNSaMRoXt28fJzR1BEKKPY0N1cw0ODlJXV4fdbqe5uXlaoupO58TCr6naHU+3wiMkGHGjKQepYuSZmRsQxVQEYXJxPtqbaqpj1OF6o9GoFdNjXW4dDrm5cOWVMldcIfHGG2/S1ZVLR8cSXn/dyM9+pufRRwUEQeHtt9P461+XsGmTl8HBveh0Os3aZGZmapvzYk3izBTTIe7zz+v52tdSmTtX5oUXHJSVKbS1Tb/byb+ba86cOezdu5esrCwGBgZoaWmJqZsr2OJGar443bI1kGDEDYa/q+zf5FBVVUVGhn5SAwZE3088FcbGxqipqWGO34rO6WyljwQ10SXLEldfvQRRFLnjDg+Dg27+/nc9W7eaGBoSuO02X/BVVHQWW7YIXHhhGeec40FRfL3DahJHvUkzMzMTIg5W4a8LtXGjl9/8xqEtDY93ecpfVN3lcjE4OBhVN5fTObHwK9xkECSJGxVU4vrv6VGbHEwmGB0N/zczyT52d3fT0tIySQY2ltrsVMc6nU4OHjxIUVERIyMjATdvejps2SKRm6uwbp3EXXe52LNHz4svjvP881aeflrAYFA466w0Nm3KY/NmLwsXurDZfEoX6rZCq9VKTk7OadmMpyJYF+qxx5wBwvSnsq6ckpISspvryJEjk7S5Yolxk8SdAqoUzb59+ybNs/qSU5O/8Fgzxf5QF1OrKhzB7qcgRN6RG3xsOOKqG+vUtZwdHR0hj1MbMObNU5g3z0N1dT3z5y/jgw9S2b1bz549Ou67L4X77kth1iwTmzZlsHlzMeef7yUlZYyBgQEaGhrwer2apYl1RWQ4REO4kRH4/OcDdaGC/2Qm/cXB55lqn22kbi67/Vzc7mGGh326Z4m8NwgSjLjBg9rHjx/H5XKxYcOGSQkBk2lyrzJEzkRHgmrVrVYrq1atCnkzxWJxw5G8vb2dzs7OAOXKYPde/RmqcyolBT72MYmPfUzi29+GkycFXn9dx2uv6dmxQ8/vfmdAEBRWrzazaZOVCy7wsnKlh9FRG/39/fT39yOKIoIgkJOTg8lkOiVudWenwOWX+3ShnnjCwec+F3pQYCYTPf6Itd1xcjeXjvR0PV1dx7HZbPi2InRqS7lVJF3lCFBVHtPS0jCbzSGzeD7ixsfiqioVU23fm4nFlWWZw4cPI0kS1dXVAU/08PuDlIA6bqjjCgsVrrnGyzXXeJEk2L9fZM8ePbt363noISM//GEKWVkK552XxqZNhVRWniQzcxhRFGlubsbpdAYkuaJR/YfIFjeSLlQs54kFM+lTliTweARycswsXryYvr4+BgYGkGU54DMym8309/dHXQ6aSm8K4A9/+AP3338/giCwYsWKSdND4ZBwxFVJpIqRv/POOyGPM5tnbnFVKdgjR45EpVIx3eSUy+Xi0KFD5OfnU1ZWNulGDU/c2MTidDqorpaprnZz111ubDZ44w09u3fr2LNHzx//aADKmTfPyZYtIps3l7J+vQeXy+cytrW1IYqiZolUsfBQCEe4117Tcd114XWhoj1PrIiP3pSiXVNqaiqzZ8/WVpyMjIzwwQcf8Nhjj+FwOHC5XNx8881hH/TR6E01NTXxve99j7fffpvs7Gx6eyOO4AYgoYjr9XppamqaRKJQX67JpIQkbrTkkiSJ+vp6ZFlm9erVUTVwCIIQ9UNBjc/VrQiRrHk44gZPB8V6g2dnw6c/7eXTn/aiKC4OHxZ56aVx3nzTxC9+kcZjjxkxmRQ2bEhj06Z8LrjAS2mpL8nV0dHB6OiotiUhJydnys9I1YVaulTmv/87tC5UMOLlKs/kAaAqPPpnlf09IlUwbtOmTWzduhVRFJk7d27EzyMavalf/vKXfO1rX9M0v9Qmk2iQUMQ1GAxUV1cH/Fu48k5qqs9VVpQJYTWIjrjj4+McOnSIkpISnE5n1F94rDGuuu9nKmseSspVEISQMe50+2QFAZYskbFah/j853vJySnl73/XaW71XXelctddUFpqZvPmDDZtKmHjRg+iOKbJ1agi4cHjeMG6UM884yDKPpWPLDkVCaEUHiNllcvKyvjUpz4V8ZzR6E01NjYCcM455yBJEvfffz8XX3xxVNecUMSF8LpTwR+kygP/LeL+x4eD2tu8bNkyLBYLPT09MW2lj+ZYWZbp6urC6XSyfv36KePGyK5yZCnX6SItDS66SOKiiyTARWur8CGJdfzhDwZ+9Ssjen0qa9ea2bw5h82bvSxb5mZ42Nc3PDg4+OGAfD733z+HP/4xJSpdqGAkgqscrDcly3LY7yye5SDVw3zjjTfo6Ohg48aNdHR0ZCmKMjTV3yYccYOhEjG4NOOvO+VP3HDkUhSF5uZmhoeHA3qbY4lbo7G4qoay0WgMu7snGJGIGy+LOxXmzFHYutXD1q0e3G7Yu1fHnj06du/Ws21bCtu2pZCba2LTpnQ2bSqkqEghM7OArVtz+OCDVG64oYkbbhhmaCgn6vetvp/TTdxgvSlJksLWvtWl2lMhGr2pkpIS1q1bh8FgoLy8nAULFtDR0TEfmLyTJghnDHGDoT4dx8fBv7881PFqljojI2NSb3OsCadIx46OjlJTU8P8+fMBTSd32ggmbjxUNaKB0QgbNkhs2CBx331uensF9uzRaRb59783ANUIgk9M75lnHHz60/kMD6cwMDBAa2trQLklUrthYhDX9zOapdbRWtxo9KY+9alP8dxzz/GFL3yB/v5+1XVuieaaz1jiTrjK4SVaAW0vUEVFRcjgP14WV+22WrFiBenp6fT19cUUD4deQzLzFSTxQH6+wtVXe7n6au+HUj4iP/jBKNu3F/C5z7n413/1AoHi6sHthv5TPP5JnZkQzh/xWj8C8dmNG43e1EUXXcSuXbtYsmQJOp2OH/3oR1x11VVRCX0nHHGj1VaOJNGqbrDv7OzkxIkTrFixIuyHPVOLqygKTU1NjI2NBXRbxWId/Y9V1RV83TuT67inGz4pH5krrjjJ9u0FXHBB6M8uuN3QX5NZTXLl5OTELTk1M9ka389oLG4snVNT6U0JgsDDDz/Mww8/HPM1JxxxgxFutC+c7pROp8PpdFJfX4/X66W6ujpivBWrxfU/1uPxUFNTQ2Zm5qRuq+kSN/D1ProYN1ZIku86opmcC9Zk9nq92Gw2uru7GR4e5siRI+Tm5mK1Wqc9MpdoFvdUQhCESxOeuOGVHkNrK0uSREdHB+Xl5ZSWlk75FJ6uq6xOD82dO5eCgoKIx06FWIg7U8QrplS/ElGM/UGi1+u1RdPj4+PMmTOHkZERbYNgVlYWOTk5Mc3UJlqMe4rxnwlP3Kld5YkYd2BggGPHjpGdnR2VSgVMz1Xu7e3VdtyGG6yfSXvk+Pg4KSkpiKJvuibccacT6sLtmYaniqKQnp6OxWJh9uzZWvO//0xtTk4OVquVtLS0sA+deFrcqeZx47mOdTpQFCUn4YgbrWCcf4yrDiT09/ezcOFCbDZb1K8XK3FHRkZwu91UVVVF7JyZDsnUmePe3t4P9aHXIElGnE5nxG0ApwPqRzZT+aVgDyC4+d/pdPpN8NjJzMzUfu9fIoxUe50KwRY30lioy+WalkxSvJFwxA3GVKs27XaZgwcPkpqaSlVVFaOjo/T390d9/miJ6/V6OXz4sNYiOdXTPRbiqgm1Q4cOkZKSou3QTU9PwWaTtPlR1RqZTKZTpnYYLSZc5ZmdZyrXPTU1laKiIoqKilAUJWDJNKAluWYygx3K4kY61+n+7OEMIW4kUfSGhlY2by6gsLAQmJ7k6lTH2+12ampqKCoq0sbiojlvtMRVlTBKS0uZPXs2brcbRVHQ63XodHpWrlyJJEkcOnRIu2lNJpPWQ3w6rLFK3HgIHkYbcwuCgMViwWKxUF5erilAdnV10d/fj9lsxuv1TrmvKBgTFtf3MxxxE0VhExKQuNG6yqOjPUApVmsxhYWGKY8Ph6mI698iaTKZop7giDbGHR4epr+/n/nz5wf0toKPFCr3dTodJpOJkpIS0tPTGR8fZ2BgQNvNo1qej2qTfLxc5ZnAXwGyubkZk8mE1+vVPJSsrCxN4SJSkssnWxMoWn+qG0ZmioQjbjCCy0HqxoLhYSdQiqKkAhNEnY7FVeu+/lDVHfv6+qiqqiIlJQWPxzPjTLE/1E2CeXl5IZNc4bLKgiCQlpZGWloapaWleL1ehoaGNO1hk8lETk7OKZOtgfglp+IFRVEwm81kZ2drChdDQ0NakstgMGhJrmCFTIdjwtqeKUh44vqXg9S51tzcXKqrFyIISsj9QbESN5hgkiRRV1eHwWCgqqpKs2DxEotT+6ZHRkZYu3YtTU1NYQfpg4cMQh2n1+s1kTRFUTRrrMrWqMkcS6St3TFiwlVOlCx3oHur0+m0hxcE6jE7HI6AJJfTmaLlTCLB4/F85Gqa4ZDwxFVdX1WnadGiRdqXYTZPTHYEHx8tgsnocDg4ePAgs2fPpqSkZNKxM7W4kiRRW1tLamqqloSKZx03lDW22Wz09PTQ2Nj4odCeifz8/BlZY/W6EsXiTpVQ8tdjVgfj1ZbMzs4l6PVZjIyMRNTNTpAaLpCAxA2lrTwyMsLRo0cDdJogtHzNTJJTqmZzuE1/sWaKg69DVXYsKSkJeCjEMkgfa4nJv9lBURTa2toYGRmhoaEBSZK02DgzMzOm2DjRXOVY6rjqYHxWVhZz584lNdWIyYQmHuByuTh58iRWqzXg4Wa32xNCDB0gQT720JAkiaamJjweD2vXrp30oaWmTm55jDVxoLribW1tNDc3s2bNmris+wgm2dDQEB988AELFy6cZMkhdCtjvDunBEEgJSWF7OxsVq1axcqVK7V1J/v27aOuro6TJ0+GzOIHQyVuAqzRAWaWNHK5fEJxS5Ys0fIZbrebhoYG3n//fZqbmxkcHGR4eDhqi7tz504WLlxIRUUF3//+98Me9+KLLyIIAvv27YvpmhPO4qrwL8F4PJ6QT1OzWZnkKscKQfAtDsvOzqa6ujpuGVl/4nZ1ddHW1hZWCSOcBpUgnLpBephsje12e8AggFpuCiWwPpHtjuslTRszHaT311Q2Go2UlZVRVlamJf46Ojq47rrr8Hq9PPbYY1x++eUhW10hOr0p8I2B/uQnP2HdunUxX3NCWtze3l4OHjzIkiVLtDa4UPBt7Jv+6zidTo4cOYLBYGDp0qVxLaOornJjYyPd3d1UV1eHfVpH6yrDzIcMwlkmX8NHOmVlZaxevZqVK1eSnp5OV1cXe/funWSN1QfKmegqB8PhEMLWcNXEX2VlJT//+c9Zt24dOp0uYlnQX2/KaDRqelPBuOeee7jzzjunVYdPOIurLrZSVSoURQl7s4aTaI0GarKrrKyMwcHBuFszSZKw2+0RdZpVREpOTSXPeqqg1+u1Gqm6BcDfGrtcWR9eY2JmlWOBv8WNNBk0Pj5OQUEBX/nKVyKeLxq9qf3799Pe3s4nPvEJfvSjH8V8zQlHXLPZrGVbIbJ7GE4UfSr4i5JLkhRTi2Q0UDPTer2ehQsXTnl8vORZTxX8twDMmTMHj8fDq6/6ZFnq62uQZVGrkcbSxxvPh9Cpsrj+GB8fj8tInyzL3HrrrTz77LPTPkfCEReij+VMJgWbLfoOF1mWta4aVZR8fHw8rjeQasmXLl1KQ0NDVH8T3lWevJE+EaaDDAYDguC7dVaurGTWrBHNGiuKQnZ2Nrm5uRF1mSG+XUgzk2eNblNftOWgqfSmRkdHqaur47zzzgN86imXXXYZr7zyClVVVVFdc8IRN5YPP5zFVeNL/y9Abd7Iy8sL2CY/3ZUloaAqbgSXrabCRzmPGy/4Z5WDrbHaT63qMquNEMHNC/FuH5y+AoYQMEQ/0019U+lNWSyWAC/vvPPO46GHHoqatJCAxI0F4WJclYwqcSOJks9kSZgKRVFobGxkfHx8SsWNUAh3wyUycRUldHLKYDAwa9YsZs2aFSBZU1NTA6BlqjMyMhKm79df4ncqVzmaxebR6E3NFAlJ3GALpDbsB3+gU1lcmLoUEytxgwXavV4vNTU1ZGRksHLlymnfiKGuITjGTRRXGaKbDgqWrAm2xmlpaXg8ntPaSqgoaow7dXLK4XBoU2hTYSq9KX+88cYb0V/wh0hI4gYjnCh6uDUkOp1OmxJxOBwRreB0iKu64epGhDlz5kT9hYY7Z+hrS3yLG0sdN9gaDwwM0NTUFGCNc3Nzp9weH0+4XL6f0VrcROmcOqOIG0oU3aeAEbiGBKCuro68vDwWLlw4ZSkmFqhEHx4e5vDhw9pGhJkg2JLabDZEUUQUU7RYMtRxpxMzbXlU68ZpaWlUVlbidrsnybmqmepTaY1DqV9EKgedbqE4FWcMccMpPUqSgMeDtuVcjanmzp1LeXl53K9FFEU6Ozvp7u5mzZo1cRliVwmpKAotLS0MDAxgMBjo7ZXweIq0zq5EQjzmcf1DDqPRSEFBAQUFBVps3N/fT0dHB4IgaLFxvK1xKPWLSOtHksSNgGDLEk7pUQ1ZHQ4fcVVR8ry8vLiOsKlQ2wIFQZi043YmEARBGyUURZHly5cjyzKFhWbAN2TR2tqK1+vF5XKRlpZ22l22CYs7fQ8g3KY+/9gYmGSN1X2+VqsVvV4/Iy8klMUNNzWVnA6KEVPrTin09DQxOjpKdXU1LS0tcSvxqFA1lEVRZNGiRXEjLfhulhMnTmj7WCVJ+vCmVlAUgdLSUubMmcOxY8e0Nkq32z0t1Yt4udrxmA6KVgw92BqPjPjqxidOnEAUfRsUZFmeVpY6Fr2pRNlGD2c4cVUvdd++ehYvnhBZC+daTxfj4+McPHiQuXPn0t3dHbfzqudubW3FarVSWlqqJb70ej16vQ5Jmoir1TG8vLw8Lc7u7e2lsbGRtLQ0TVR8qu6l+Ozq8f2Ml6scLfx1p+bOnYvb7aavrw+3283evXsnWeOpEIumcjLGnQLR6k6JohMwkJVVyPz5OX7/PvParAp1RldNQqnSqdEi0s2pdlmVlJRoKo6iKGrH796tw+WCX/zCxLJlzWRnS+Tk5GhW02KxkJWVhSAIOBwO+vv7qa2tBdCaHk5VhjYeY33xqOMajUby8vLo6+ujsrJSy3Go1lj9HMJpMp9pWwxUJCRxgxGKuL29vXR324BKzOZc/Bd/TacbKtRN1N7eTldXV0ASajryNaFuGLW+vHr1arxeL83Nzezdu/dDFf9cfvazAurqdIiiwp13pgBLWbpUZssWiS1bJNas8QI+l1qWZVJSUigqKmL27Nl4PB6GhoZoa2vDbrdjsVjIzc2Na4IrHq5yPLfR+7LwYoA1VhePtba2Mj4+rlnj7OxszRonLe4phD8RVdHwoaEhli9fBYRe/DWTpgr/nuaqqqqAJ/BMdwKp1z88PKzpMxsMBiorK5FlmZ6eIb72NTPbt6dwySXdfPWrtSjKHGpry9i5U8cjj+h56CEDOTkKF17oI/HHP+4hM1NBkiTN1VZrojCRaT9+/DiyLJOWljajPT2+z2jmFjee2+hDncd/8ZgqVzMwMEBbW5umSWWzFQLmqCyu2+1Oak7FAjVm9Xq91NbWYjKZWL16NXv3+j7gULpTLrWyHgVUoouiiMfj4dChQ2RnZwf0NAcfGw2CiSvLsiZCt2LFCu1hob7G0JDI9dcX8ve/67j7bjvnnVeH1WplfLyb9es7+MQncjAa89m7N4udO/Xs2qXjuef06PVGzj57whrPmydpJFZXfKSnpzN37lzNCh89ehSPxxMwLB9Lgut0xbihEM1kkL9cDfh61wcGBmhv7wesnDzZQmFhWsReZUEQEkIMHRKUuKHIYrfbef/99ykrK6OoqAjw39gXeUfuVFCPt9vtHDp0iHnz5jFr1qywx07H4rrdbg4ePMisWbMoKSnRrIT6XltaBD7zmRTa2gQee2yQxYsPsnTpxG4itV2wr+8ExcX13HJLJt/+dh4tLXns2mVkxw4d3/qWkW99CyoqZC6+2Efi9eu96PWyltySZZnMzExNNG1oaIju7m6OHj0acSAgGKrFnQnvTudSazW0yMvzUaC42MrISB9DQ0M0NDSQm5tLTk6OJuUaaS78dCAhiRsMu91OZ2cna9asCajPRtqRG6vSY39/P8ePH2f58uVa/TAUYlnmpZJcfSBUVFQE7IRVb9r//V+Rf/u3FBQFfv3rLmbNaqKycmVAc0dwu6AqpK7XH+eTnzSwdWsuTucs3ngjjR07dPzyl3oee8xAZqbCpk0SF18ssWxZOx7PMEuXLtXeQ1ZWFtnZ2QiCwPj4OP39/dTU1CAIQtjEjs/iCoiicsYSV4U6pFJQYCE3N5OhoSEWLVrE0NCQJuVqsVjo7OwMSBxGws6dO7npppuQJIkvfelL3HXXXQG/f/jhh3nqqac06aBf/epXUS+pU5HQxFVFyXt7eykoKJjUVBFu1Was5SCn00lbW5smFBYJsVrcwcFBjh07xvLlyzGbzZNI+z//o+NLXzJSUqLwk580k5XVx/LlqyOWMgRB0Ny+iooKLaNstzdQWeni/PNzMJvz2L8/m1dfNbBjh8hLL+kRhHlUVZVzySU+Ii9dKiHLEwmu1NRUSkpKKC0t1dZ7qIkdNcGlupqKIsRl4dfp30bv++mfnDKZTJjNZoqKirSy27PPPktbWxsXXnght9xyC1u2bAl5vmj0platWsW+ffswm8387Gc/44477uD3v/99TNedsMT1FyVftGgRPT09k46JZHGjIa4sy9oKj1WrVkWlMxyLxXU6nbS0tLBmzRp0Ol3Ajaoo8PDDeu6918j69RLbttVgtcosXLgi5pvQZDIFNG8MDAzQ399Fbu5hvvzldK67zkV7ey5HjlSwc6eOBx4w8sADUFw8ERefe66X1NTABFdubq4mXaO2IB47duzDrQ7FM9abOtXJqWigWlz/PJ3/udQGjwceeID33nuPp59+GmcEhUJ/vSlA05vyJ+7555+v/ff69ev5zW9+E/N1JyRxHQ4H+/bt00TJR0ZGotiRO4FoykFut1vbihBLmSSah4KiKDQ1NeF2u1m1ahU6nS7Ayno8cMstRp55Rs+//qubr3zlfYqKrFEt4p4KOp1O04ryer0cOHAAURQpKemltLSfz342F1nO5+9/z2DnTh3PP6/nqacMmEwK553nI/LFF0sUFgYmuNLS0khPT2fevHnY7XZk2Y0gyOzbt09TvAilBjnV5xQP4s7EcjudYDAoU3oPqqZy8H6nYESjN+WPp59+Oqz1joSEJK6iKCxevFhzy6bakRvrGhJ1m3xFRQX5+fnU19fHFLdGOtZ/U4HFYmF0dBSz2azdWCMj8LnPpbBnj45bbhnnE594j7lz54RNhk0X6oOpuLhYS+a5XC76+/vp62tkwQIn69Zl89BDuTQ05LJzp54dO3Ts2OHruqqsnLDGq1Z5AV9ya6JJRIfBIFJZWcnw8DBdXV0cOXIkpqmeRIlxo6mKnYp2x9/85jfs27ePN998M+a/TUjipqWlBbTthe+cgjXwVUUAACAASURBVJSUydrKkSyuun2vsrJSy9jGayeQy+Xi4MGDFBYWUlxczMjICJ2dnbS1tX04NF7AjTfOprFR5OGHh6ms3BfwgIoXxsfHqampYf78+dq6FvBlUtU1HOpal/7+PtLSGrn22jRuuSWXwcF8du9OZedOHQ89pOcHPzCQm6tw0UU+Ep9zjp3W1iOYzSsQBEWLt9UEl6rN3N7ejiiKk7Kz/pgJ4eJ1HqdzIr6NhGi7pqbSm1Kxe/duHnzwQd58881prYJJSOLGUjs1m6OLcdX1G729vZr0azTnD0a4jPXY2BiHDh1iwYIFWK1WZFkmIyODxYsXoygK77zj5JprMhkfV3jwwQMsXtzLggWhV53MBCMjI9TX17N06dKI2XE1hlUXhdntdvr6+nA6D7BhA3zqU7nodHm8+24mO3fq2b5dx29/q0evN7Bu3Xpk2YAgCBiNRs2lVhM7aoLL7XZjs9m07KyvK8zXuaQm+U43caO1uHa7PS56UwAHDhzgy1/+Mjt37iQ/P39a152QxA2GXq+fQhQ9cowryzINDQ0oihKwfU9FrMQNXsvZ39/P0aNHqaysxGQyTcoc79yp47rrrFitCr/6VRsZGaPk55fR3t5OU1MTVqtVG0WcyY3c399Pc3MzK1eujFmsTm3SKC8vx+12MzAwQF9fC6Wldu64I4t7701n165RmpsX8de/mjl82Heda9aYtZrxunVe9PrABFdeXp5WwhoZGdESXKmpqYiiGPHhEi3iZXEjJbkcDkdUn2k0elO33347Y2NjXHHFFQCUlpbyyiuvxHTdZwRxI2VyU1MjW1y18SE/P5+ysrKQX0ysrrL/sf79zMGZY4Cf/1zP7bcbWLlS5vvfb8BicbBkyRpEUaSsrAxJkhgcHOTkyZNajJiXl0dOTk5MonNdXV2aVnQs2sahYDQaA1oFOzo6aGk5xvLlRtatq+Gmm3L51rfKeO01I7NnK/z853oefdSAxaJwwQW+5NamTR5ychSt8QN8apBqAsvpdNLc3ExXVxe9vb1ai+ZUkq6hMJNY2d/ixmukbyq9qd27d0/rWv1xxhA3HELtD1Lj0NHRUWpra5k/fz55eXlhzxGrxVW7aBobG3E4HNqmAn8rK0nwrW8ZePxxA5/4hIdbb91Pbq6JefOWBrwf1Sqp+3tGRkbo6+ujtbUVg8FAbm4ueXl5YZ/2aq1b7X2O55wwwODgIF1dXaxbtw6TyYTdbqe/vx+Xy0ZKSjaPPNJIamoe+/Zl8eqrvgTXCy/oEUUj69ZNdHAtWhRYMzYYDKSnpzNr1iysVqu2n2d0dJSMjAxtRDGah1e8LO6ZIlsDCUrcWLWVg11lQRDweDzU1tZSWVlJenp6xHPEanG9Xi8HDx4kLS2NZcuWTeo5ttvhi1808uc/67nxRidXXPEepaUT2d1I51anW/wbKw4fPozH4yEnJ4e8vDzNaimKwpEjR1AUhcrKyrj30XZ3d9Pe3h5gxdXdu1arkZQUkbS0NPr7T5CfX8/Xv27h3ntzaW9X2zBF7rvPyH33QWmpj8QXX+yrGXu9vpUmeXl5WoJLTaSp605OnDihDQPk5uaGtXgzjXHN5ugmgxJliB4SlLixINhVVi2Q2+3mrLPOisptjIW4kiTR1dXF/PnzNVfSn7Td3XDFFSkcPCjy3e+Osm7d+5Oyu9HCv7HC6/Vq2VrVKtntdnJycpg3b17cZ27b29vp6+tj1apVIa2eJPkGDPyVKYaHh+nr60OWj7Nli4Frr83D48njb3/ztWH+5jd6fvELX8145UoH//qv1cyZY8BiCawZm81mzGazJrCuxu5Op1NT/cjKytJINlOLa7Wq7ynyLG6SuFEg2vE5s1lhYGBiHK++vh5BEDCbzVHHeupU0FQYHR2lqamJrKyskKRtaPANCgwMCDz11ABlZTUsW7ZsSosfDfR6vdar7HK5OHDgAEajkf7+fkZHR8nLyyM3N3fG4nWqYJ3dbmflypVhCSHLgZNB/m2YgOYtDA01sHixh3POsZKRkcc77xh54QUn+/cXctttOm67DVaunKgZr1gxuWacn59PQUGBNpqnEtlkMpGTkxN2DWs0cDgC149EsriJJNiXsMQNRiRR9PHxiRpqQUEBpaWlvPvuu1GfOxqL29fXR1NTExUVFdq6Sd8OHR9pX39d5JprUjCbFX73uw4slmOsWBFdG2UscDgc2gSTGrePj4/T19enbdJT4+JY1S8UReHo0aMoisLy5csj/q1PUif8uYK9BZ8QegsWi40778wlP99Nb28er72Wws6dOn7wAz3f+56B/HxFi4s/9jEP6emBCa7MzEwsFouW4Orv76e/v5/h4WHy8/O1LQnRvm+nc2Lh11Qxbqh67OnCGUNc/5lZf/i20ivs27dv0oqRaLONUxH3xIkTnDx5kjVr1mhJrwMHDmAwGMjLy2PXrmJuvdXEwoUKDz/cSHr6IMuXr4l7okhdFrVkyZKAgQuz2awtYlZdy+PHj2O328nKyiIvLw+r1RrRKqmzwmazOSrXO9jiRoJer9dyAxs2bNAI53Kd4PzzdVxxRS6imMfbb2ewY4eOl1/W8etf6zEajZx77kQbZlmZpAnpqQuoi4qKsNvtzJo1C6/XS3t7u6YEqbazRkpwRWtxk67yNKHWZoO/BEFwMDLiG0z3d0kjycYEIxxx1eSP2+1m9erV2nkrKiqoqKjAbndwzz0KTz6ZTlXVIPfcU4PFYmLp0vgnilTV/xUrVkS8gQwGQ0ApZ2hoSPMWzGaz5lL7hxHqGpXc3FxKS0ujup5YiHvy5Ek6OjpYtWoVBoOBlJQULBYL8+bN87OajcyZ4+Tuu7P54Q/zOHrUys6dBnbu1HHbbUZuuw0WLZrIUq9d60WnU3A6nYyOjlJaWqqRFXwJrv7+fi07r44oBn92/hY3mZyKA4Jj3OCmCkVROH78OB6PBUkqnBRHhrPQoRCKuP47gebPnz8pc+x0wle/auGFF/Rcd52bz3++jowM39jee++9R3Z2Nvn5+QFJlOlCvfFjrdGKoojVasVqtQZ0Rx06dAhAa/poampi9uzZMa1RkaTohug7Ojro6ekJm+RSRwlLSkq0Nsy+vh4MhqNcdVUa3/hGHsPDeezZk8qOHToef1zPf/yHgexshY9/3MOiRSe44oqFZGRkTEpwlZWVBTSUqIMfWVlZ5ObmYrFYAizumSIUBwlM3GD4E1eSJOrr69Hr9ZSVVYTdHxTKQodCMHGdTicHDx5k9uzZWlLEn7T9/XD11Sm8846Oe++187GP7WXu3HKtfU2WZWw2Gz09PZqyhBp/xdJUobZp2my2sDd+tAjVHdXV1cXBgwcxGAyMjIyQkpIS9YNGzSpHgnrtK1eujCpsCG7DVK3m6OgB1q6FSy7JJTU1n//93wz+8heBHTsEXnxxKd/7nsJZZ01Y4wULfDVjNTZWa+Xqd6nK2h450ozX+3E8nlHc7tA5FBVJiztNqET0b+QvLS3FbBZwu4VJN1KsTRXqsSMjI9TW1rJ48WIsFssk0jY3+zLHHR0CP/+5jYqK/SxeHBhz+suCqjFxb2+v5rapDReRMsBqokiSJFasiH1Gdyp4PB66u7tZtWoVGRkZkx40avdWuAkfWQ4/SK9mpsfHx6ddXxYEQdu7qz5ofJNNx5g1a4zPfc7N3XfPobe3TJPuueceI/fcA3PmTMTF55zjxWicaMMEX4IrKyuLWbN836le76GurhGn00lGRgYWi2VSYi8Z4/7/9q47LIpz/Z5ZtlCUIk0FFQRRpFuxR8y1K3ij2KIYK5qYXE00icYETaJGY7zXWFL0pzEmRkEjil1jiy02EAWRIr3sglTZwu6+vz/WmSywwAKLovE8D48uzM58M7Nn53vf77zn1RO6vJVLS0tx//59dOnShVsX1S6m154t18eilSWuWCxGUlISfHx8YGxsXI20V67wMGmSCAwD/PJLLlq1egBv79p1wdrtNFhRhUQiwf3796FSqXRmgFkTgRYtWtTZtKwhKC4uRlxcHLy8vLgQQ/tJV1paColEwgkg2DFqf3A1WeXqy3WsokylUsHT09NgY2cTUVZWVpw5vVQqBXAdY8caY9YsW8jltjh/3gTHjxth1y4+tm0TwMzsb+ueYcNUsLP723urvFwzfnt7c/j6+uLRo0dQKpVIT0/HkydPKiW4ysvL9VrWq8u2Ri6XY/r06bh16xasra2xb98+ODk51ft6NFviVoVUKoVYLEb37t0rxRra9jXa17W+T9zS0lLI5XJ0795dZ+VKeLgR5s4VokMHwubNj2Bikg1v7271tus0MTFB+/btOXsY7Qxwq1atYGVlhbS0NLRp0waOjo712rc+YGO9mgoRtL9otJNHCQkJkMvlnHpLpbKDkVFlUhIR4uPjYWRkpNMhs7FgyxW1s+pExC2HPX58F56eagwaZI0WLWxx544lTpwwwvHjRjh8WPNR795dxT2NLS01nw9TU03vJratCys/Zafqn3/+ORITE7F3715MmDChRqLpY1uzY8cOWFlZISkpCb/99hs+/PDDetvWAEDz8JqsBawPcWlpKTp27FgtQVBTMb2+T1y1Wo3k5GRUVFRwggNtG04iYP16PmbMEKFHDzW+/z4WlpYFXIa0MWAzwN7e3ujduzfMzc0RFxcHmUzGTV2rViI1Bnl5eUhOTka3bt30rh5ik0d+fn7o2bMnZ5xWWFgMubwcYrEYSqWSW04SiURwc3MzOGlZw72qS2EMw8DMzAxOTk7o0aMH/Pz8YGZmBrE4DRYWlzF3bgwuXcrAlStPEBamAJ8PfPmlAP37m2DgQM01SE/nc21MbGxsKrl+ODk5YePGjbCwsIBQKMQvv/xS4xi1bWuEQiFnW6ONyMhIhISEAADGjx+Ps2fPNsg9slk/cdkpo1AoRLt27XR+GNjPn6aVRP0sWpVKJWJiYtCiRQtuGljVYua994T46Sc+goMrMH/+TdjYtISzs0eTfDBTU1Ph6+sLCwuLasUG+sTFtSEjIwNisRjdutVuRFcbtG1xWrQQQSZTori4GCkpKZDL5bC0tISDg0OTXJu7dzUqNNb8oCYIBAJOhskmovLz8yGTpeD114WYNMkGgC1OnzbDihUCAIS4OBkkEgl8fHy4JBqb2GKz8ZmZmZg/f36tWX19bGu0t+Hz+bCwsEBBQUEl/YE+aLbElcvluHnzJhwcHNCuXTtkZGTU2WpTG3VZtEqlUkRHR6NDhw6wtbXFnTt3cPPmTdja2j71azLD1KkinDtnhA8+kGLkyOv1XjLRF48fP+ZcOdgZRdVig7ri4prALpuVlpbqnd3VB2o1IBLx4ezsjJKSEtjZ2YHP53NjZKfUDSnT00ZZWRliY2MrxeP6gjV6Y6WK7HXMzHyAHTvcUVJiic8/z0C/fmnw8al8bdiWJjKZDDNnzsS7777b6HJJQ6LZEreoqIhzkwBq7k7AlmTpmirX9MQtLi7mFEjm5uZQq9Xo1q0bF3NeupSO99/vjMxMHr76Sozu3WPRqVNnbiyGRG5uLtLT02t1mawrLra1tdW5jMNmptVqNby9vQ36JFSrAT6fcOfOHTg6OnJfaOwY2VYfZWVlsLCw4NRb9fniYJVi2l9ojYGJiQlsbNojNLQT7t7lYeXKVPTokQCVio/4+PhqddByuRzTpk3DyJEjsXDhwjr3r49tDbuNo6MjlErNjKUhBSjNlritW7euFN/VZRinr0UrG+f5+vpCJBJVyhwLhULk5jpiwQIXyGTAtm2P0KZNPNRqIfLy8kBEnO2KIZCWloaCgoJ6TV+rKqMeP37MLeNoF+HzeDzcv38fxsbGTZKZrqggyOUl6NChQzX7larT1aKiokrOF6x6qzYdN2vBU5dSrD6QSjWVWxcv8rB+fR569cqEr28/8Hg8LjRJS0sDEeHUqVOIiYnBkCFDsHDhQr2unz62NWPHjsVPP/2EPn36ICIiAgEBAQ26N82WuLqWg2qfKleOcXUprVJTU5Gfn48ePXpwyixtEkZFGeGtt4SwtSXs3JkOU9MM+Pj0A5/PR1FREdeLtqGCCu2xPHz4sFJCrCFgzdjYZRz2w/fo0SOuBK6m3EBjoJEZVsDe3gR2drXH3NrqLQCceis2NhZExE37tbslFBcXIz4+vt4WPLVBKgWCg0W4cIGHdevE6NUrsVLowIYmgMYC6M6dO8jIyMAvv/yC1157Dd27d6/zGPrY1syaNQvTpk2Dq6srWrVqhd9++61B58PUM6P1zJqnqNXqSqV2hYWFyMnJqZRaB4D4eAZ+fkL8/HMFJkz4+wmblZSEdkOGgDZuhGrkSMQnJ4OI0LlzZwCVk1AAsGULHx9+KED37mqsW5cAY+NieHp6VpvaaQsqCgoKuP6sdnZ2esVAbPbVxMQErq6uTfAkrEB0dDRXoK6pj1Vz5uY19YnVF+ySzOLF/dG6tREOHtS/uVpVsKIKTReGJ7CysoKJiQmys7MNTtqJE0X44w8evvpKAn//hBrjfaVSiTlz5sDDwwMrVqyAVCoFwzAGG4se0OvmNNsnblXUNPVlY9xKU+XcXDhMnQq+RAL1+vXgzZ0Ll0GDIJo4EUoHBzBaCROVCvjwQwG2bRNgzBgl3n//DiwsBOjUSXdMWFVQUV6uWRLR1v/a2dnpnN5VVFTg7t27sLOzq9NYuyGQyWSIiYlBx44duZI/7Yqh5ORklJeX1xoX1wY2UeTp6QmAr1OAUR+wogq21Ud6ejoePXoEoVCIxMREbkrdmGU3mQyYNElD2rVrJejd+0G1RBQLlUqFBQsWwM3NDStWrODqupsjXhji1j1V1vzL3LsHwb//jeLx48HLysKNZcvQvkMH2F29Ct5PP0H09ttQDRgA1ZgxKBk0EjM+aItjx/h45x05/v3v62jd2r5epGKdGpycnCCXyyGRSJCQkACFQgFra2vY2dmhZcuWkMvliImJgbOzc4MtOWvDkydPEBsbi86dO1cr+NY3Lq5t2s/GnGyiSCMxNdwEjF23Zl1LysrKIJFIuE4M2lNqfcGS9swZI6xdK4G//wP4+vrW4OihwrvvvgsHBwesXLnS4DMhQ6PZElffGFe7DQlz6hQEs2ZBuX498vz90e7XX+HWvj1MPTygcHcHM2sWmOJiGJ04gYr9R2Ay/wMchAJbV6XDp18cOnb8uzi9IRCJRFylC2s1k5aWhpKSElRUVMDZ2bne63X6gCWVPuucNcXFqamp3LTf1ta2UuKosLAQCQkJlaav9SnrqwvsbMDPz48LN1idMttZni1NlMlketnZyuXAlCkinD5thDVrJPD3j4Ovr+5CDbVajcWLF8PKygpr1qxpNj1wa0OzJW5V1EVczwtbIfjvGlTs24ecjh2RkZiI9iYmMH3qnMHFtJaWiPGcjDUfdcJ3uAw7phzefe+ja9faDcTrC9ZqRigU4sGDB3B1dUVpaSmuX7+Oli1bcsmtxq6rsmvADYkJq5rTaTtpqFQq2NraQiAQIDMzE76+ldt+qtW1O2DoCzaZpk3aqtD+QtTHzpYl7cmTRli9Oh99+sTVWF2lVquxdOlSiEQifP311y8EaYGXgLgCAZAMZ9ifK4biyH6ktG6Nx5mZcHd3h5LPh7ygACZaiagzZ3g4OPEQtisWIjnsG9itmgE/Ly8YG8AXqiry8vKQlpYGPz8/7kPPPuXEYjFSUlJgbGwMOzu7asXt9d2/ISxytJ00FAoFUlJSkJqaCpFIhLS0tEpxsSGeuGzFVH3ko3XZ2Zqb22LJEmecOmWEL7/MR9++92ucHqvVaixfvhxqtRqbN29+YUgLvEDErU3CKGNMkdWqPZyDg9HBzg4dg4OhHDMGRubmSE9LQ65AgFatWuH0qQ4o+vh7rOFvQ/J3O6BwbwV8JYJxE8Qz6enpnEui9odS+ynXqVMnPHnyhEtuMQwDOzu7Wn2UWWgXqDdWM60L+fn5KCsrQ//+/cHj8VBYWMh1rtfE7H6N2j9r/dqY8VedMRQXS5+q3YRYuDAO3bqlo2NHjxpJu3LlSpSWluLHH398oUgLNGPiVo1xa0sWPBR6IKnNYIzetQUdMjPBHD4Mk3HjwBQWovOZM+jQuQs+2uYInx8/RKDpNcRu/hqmnVrB29tbY1olkwEGcjcgIs5K1M/Pr84PhJmZGZydneHs7AyZTKOZZXv2ssktbWmjtgG6ISWM2tC2ZmX3XzUuVirVKCoqxJ07yTrj4tqQk5ODrKysRpsDaEOhAObOtcC5c3x8/nkB+vXLQ7t2GmO/hw8fco25WXHK6tWrkZeXh507dzbJNWxqNFvi6ovy8nLIGRGMiUH7jh2hcnKCesAAKNevh2joUKjkShSPmYuNZcnIbuWJ9J83wMrWCkZGRrhx4wYG8HgoyM6GpaVlo28g26NIKBQ2qBbV2NiYc0asSdooFouhUqmaxAAdQKUvBV37Z59yRkYC2NvboHNnfp2CCm1kZ2cjJyenxulrQ1BRAUyfrjGgX7WqAP37x8LPT1Ny2bZt20qez1u3bsXJkychEokQERHxQpIWaObErctbubCwEHFxcVCJzCBUKaoVvsss7bHmQTC+LJuIHDtvFC2cAIcurmjdujWApw2RTU1Rlp+PxBs3YGJiwsWb9Z2+NcRwrTZUXcLR7mhga2uLgoKCemt/awNbPimTyeDl5VXnlwJrz1o1LmYzxFKplPPdYrO/mZmZEIvFBp0pVFQAISFCHDnCR1jYYwwYEKszPLG0tISFhQXMzc3h5OSEAQMG4JNPPsHu3bsNMo5njWZN3NqQnZ2NtLQ0dOvWDUeER8BIZZVI+/C+Eh1OnMckeoj+/+GDf9kEDu3bwfQpaQHNDeWZmMCpTRu079yZizfv3LkDPp/PxZt1TQHZNdr27dtzXwqGBBEhKysLDg4O6NChA+eZlJyc3KgvG+39s37KHh76lSzq8pzSFlRUzf6y5gTa0+/GoqICmDFDiMhIPj77rBCDBt2Fr6+vzutARNi8eTNu3bqF8PDwZlXp0xC8cMRVq9VISUlBcXExevToAQBQCU3Ak0u5D9y1k6XgBb+JCp4bzMf2Qp9Tn8E8MRHqH36AUiiE6vXX/1ZuiESATFbJTI21RRGLxdwUkK1DrZo0YoUP2pVMhkRFRQViYmK4ZtkAuI4BbK0o+2XD1svWp26XiLjpfX0kmBrPqZpnQ9rZ39TUVEgkElhaWiI6OrrG9eL6oKICeOstIQ4d4uPTTwsxaFA0fH11LykREb7//nv8+eefiIiIeOFJCzRz4ladKvN4PNy9excikQg+Pj6cZapKZAq+tBQAcGRzJjw/+jeiLQahzb7/oFyUDXufVVCFhAAtW4L//fcQzpsHVUAAVGPHgoyMwCgU1UTYJiYm3BSQFQBoT1Xt7OxQUVGBhIQEg7UZqQrWGK8mtVXVLxuZTAaxWMzVxLLkqCneZHXTrPNjfWJyfddx2ZiZtQQC/u68oG9cXBVKJTBrlhC//64h7WuvRde4DkxE2LFjB06fPo3ff//d4J0lnheaNXG1oVAo8OTJE9ja2sLJyalSPEvGxuAVS/DzezEYvX08DnZ8H747/wWBsAAeHk9bT1pZQdWrF1QzZgD5+TA6dgz83bvBi44GLzIS6l69ajy2tgCATRrFxcWhtLSUmxY2pkerLrBifl0SxppgbGxcrW6XjTdbtWrFxZsMo/FYunv3LqytrRsUk+tjz5qSkoKysrJqMXNtcbG2jlrX9WRJe+AAH598UlQraQFg9+7diIqKQmRkZKP7Ks2cORNRUVGws7PDvXv3qv2diPDee+/h2LFjMDU1xa5duzgjfUPjhSAu6zdkamoKe3v7akkoOc8ETkUP4LN9LHb334Zeq9vB3FwNV9e/+9+QsTEYtpGujQ3Ufn5gvvwSFR98AOWSJXqPRSAQcF69/fr1Q0lJCddBz1Am6GwBuT4SxtrGySa32HgzKysL8fHxMDc3R0lJCRwdHRtc7FDbE5e1Z5VKpfD09Kz1WtQWF1dVmLGkjYjgY/nyIgQE3IGvr2+NpP3ll18QHh6OqKgog1T3zJgxA++88w6mT5+u8+/Hjx9HYmIiEhMTcf36dcyfP7+adY2h0OyJ+/jxY8THx8PLywvp6elITk7mbDpZOMWfxGDVGeywWgxlDzuYmVmjU6cqvWifxrIAwDt7FqKZM6H4+muoJkzQeyxs5rW8vJzLjLLKp6om6Obm5rCzs6t35peVMBqygFw73lQoFLh16xZEIhEyMzNRWFjIkaM+ya2amn6x69gKhULvRJeucVZVmPH5Iqxf74WoKFMsW1aMgIDb8PWtWTG2f/9+7NmzB0ePHjXYdRw4cCBSU1Nr/HtkZCSmT58OhmHg7++PoqIi5OTkNIndUbMmblZWFtLS0tC9e3fw+Xx06tQJRUVFnILHwsIC9vb2uPnmWsRF9YZzaSwm/DcQ9//niZ87j0OL6WMwaJqDpv/pU6GF0c8/Q7hiBeS//AJ1//56j0WtViM+Ph58Pl9nJ7uqJuhs5jcpKQlmZmZc5re2tUuxWMzpdpsiFmNjZldX10oWpGKxGGlpaRAIBHpn0nVNlYkIiYmJUCqV6Nq1a6NCB21VVMeOnTBrFg9RUcaYPTsZffrEw9pas9YtFAqrHefgwYPYsWMHjh492iS5h5qgyywuKyvrn0dcU1NTznmAiMDn8yspeAoLCyEWi9FnSj4qJvjC2fkN3Hq0DUnfXYTVuUMYsmw1Hi3riOMu49DDUYou5ZchEmdCduIEqEsXvcehVCoRGxsLKysrdOjQoc4PpHav2KrkEAqFHDm0p3hZWVnIyclBt27192rWB2x7Tu3st3a3ABcXF85MTTtpxBbfV0VVrTK7pATAoJ7KKhUwb54Q4eF8fPhhCYYOTYOHR0+UlpZWiottbGxgbm6OY8eOYdu2bTh69KhBi0aaG5o1cUtKSmBsbAwjIyOdEshWrVpBoVCguLgYbm5uKCoqAmN8G36fK+9/2gAAIABJREFUtIbt/z5F3KONiP/uGlqeOoQuybugAh/jeiVgwB+2CGypgoND3fWk7BptQx0eq5JDu/EWj8eDra0t5HI5ysrKDLrGqQ3W3rSqJ3FVaJvSsUkjtpSOlV+am5s/TW79PVVmuxryeDyDeiqrVEBoqBB79/KxdGkJhg69yVUptWzZslJcHBsbi/nz50OhUOCrr756Lks++pjFGQrNmrhffPEFoqOj8frrryMoKKiSDI/V7BYVFXFTaTs7O85aJi8vD0rmEbwXW8B2zRLcffgRLh6RIeliGxxewsOSJUDPnioEBqoQFKSCs3N1EjfFGq2ZmRlntC2VShEXF4eysjKYmJggPT29xidcQ9HQRFfVpFFBQQGXhLO0tIRK1Q08HjV4HbguqNXAggVC/PorHx98UIJhw25WKy0E/o6LFQoFbG1tsWrVKly8eBG9evXibIqeFcaOHYvNmzdj0qRJuH79OiwsLJpkmgw0Y88pFmVlZTh+/DgiIiIQHx+PgIAAjBw5EmfPnsX48ePh7u5eY9aSFSjk5eUhPz8fQqEQ9vb2KC62R1SUCJGRRrhzR/OE8/ZWIyhIicBAFbp0Ia6/TmMyu7VBO2Z2c3Pjlm/y8vIgl8u5aWpjfInZczCUvSk77sLCIrRv74hp01IwdWoizM3N0bVrV4NN8VnS/vwzH++/X4qRI2/oJC2Lc+fOISwsDEePHm0SdxEWkydPxvnz55Gfnw97e3usXLmS80ULDQ0FEeGdd97BiRMnYGpqip07d3IioXpAr5vd7ImrDalUit9//x1Lly5FmzZt0LNnT4wbNw7+/v56TTFZlZFEIuGe0DJZaxw7ZozISCNcu6bZh6urAv7+2Zg1yxI9ewr16gNbH6hUKsTGxsLCwgJOTk7ViKlUKpGfnw+xWMwVGNjb23NrsPpAOzttaKMzlQowNzfF7NkZCA0VQyAQcMZ5uuL3+kCtBt55R9M9YtGiUowefaPWc7h48SKWL1+Oo0ePNonc9Dng5SMuAPz6668QCAQYO3Yszpw5g/DwcNy8eRN9+/bFuHHj0K9fP72qTqRSKfLy8iCRSLhYU6Wyx/79Chw5IkJMjAXUagbOzmpuOt29u7rRrg+shLF169Z6NfViYzixWIySkhIuk16bv7NEIkFKSgrnHW1oyGRqWFu3wLvvirFmzd9ZW9Y4TyKRAKjdOE8X1Gpg4UIhdu3iY9GiMowa9Rd8fWsm7eXLl7F06VJERUU1WSz5HPByElcXFAoFzp07hwMHDuDKlSvo3bs3goKCMGDAAL2++WUyGfLy8pCeng6VSoUOHTqAz2+Ds2fNcOgQH+fP81BRwcDBQY2xYzVxcd++6no7QLCJrg4dOsDe3r7e58mai4vFYhQWFuq0wGHPoyaxfWOhVqtx40YsAgL6YOVKBT74QHdTMlYmyjbT0jbO0y2/BN57T4D/+z8B3nuvDKNHX4efX812PNevX8eiRYtw5MgRgzhm1tUeMz09HSEhISgqKoJKpcLatWsxcuTIRh9XB/45xNWGUqnExYsXER4ejkuXLqFbt24ICgrC4MGDa3z6qNVqLivq7OwMiUTCdaHTrGna4/z5loiMNMKZM0aQyRjY2hLGjNHExIMGqVEXR1gJo6ESXdoChYKCApiYmIDP56O8vNygBeraYGWSpqa26NbNDZ9/rsDixXV3E2SN88RiMcrKyqopzNRqYNEiAbZvF2DhwjKMHXsdvr41C1Bu3bqFd955B5GRkQ3qLavrvNzc3Cq1x9y7d28lD++5c+fCz88P8+fPR1xcHEaOHFmrGEMXiouLa83qP8XL5ausL/h8PgICAhAQEACVSoXLly8jIiICn332GTw9PREUFITXX3+d+yZnP4yWlpZcvKmtS5ZIJMjLS4CLiwKrV9tg06bWuHSpJQ4f5mPfPj7+7/8EsLIijBypmU4HBKhQNYfCZnY9PAxnSFfVtiUpKQlisRh8Pp/zbm5M9U1VqFQqxMTEPF0S0kzx9Z1xsMZ5rFxVW2HWokVLbN7cBXv2mGLhwid1kjY6Ohpvv/02Dh48aBDSApXbYwLg2mNqE5dhGJSUlADQELBt27Y691UTPvnkE8TGxsLb25tzlGwMXjriasPIyAgDBw7EwIEDoVarce3aNRw4cABffvklOnfujICAAFy6dAkrV67UGSOxDgpt27aFUqmERCJBbm4SHB2l+PRTG6xfb4e//rJEZCQfUVFG+OUXPlq2JAwbpiHx0KEqKBQaa1NDZna1wXbjKy8vR58+fcDj8fQuSdQXKpUK0dHRaNOmDdq2bYuiIs3vGxLvayvM1GrCwoXAnj1mCA5Oxb/+dQ9t2jjXOFu4d+8eQkNDER4eDldX1wadiy7o0x4zLCwMQ4cOxbfffosnT57gzJkzeu8/LCwM169fx9q1axEaGso1vGbRkAKVl5q42uDxeOjbty/69u0LtVqNw4cPIzQ0FM7Ozli6dCkCAwMxYsSIGpd++Hx+JdG+ZukmFTY2T7B0qTVWr7ZDdHQrjsQREXwYG6vRo4cQb77ZC+7uhreZ0dYFa1fg6FOSqH8JnRLR0dFwcHDg1iRZs83GaEWIgCVLhNi1S4AFC55g3LhkuLl5o7S0VKdxXlxcHGbPno3ffvvtma/PAsDevXsxY8YMvP/++7h69SqmTZuGe/fu1VlMUlxcjPv372PdunXw8/PD0KFDcfLkSSQkJKBXr14YNGjQy9X0qynB4/GgUChw9OhR+Pn5ITY2FuHh4Rg5ciTatGmDwMBAjBo1CpaWljrfb2RkxE392KxvXl4mWraMx8KFVggLs8elS4SDBxlcu9YGoaE8vPsuYfBgNQIDlRg9WoUGdFasBFatxDBMrbpgXSWJrFSwqhqqKtg+RO3atau01MKabTaUuBrSCvDddwKEhpZj3Lhr8PHRzEhsbGwqGefduHEDixYtglQqxYYNG+Du7t6wg9YCfRRPO3bswIkTJwAAffr0gUwmQ35+fp3rxhYWFvDy8sK8efMwfPhwrF+/Hnv27MHFixeRkZGBVq1awcvLq95jfumSU40BqwKKiIhAVFQUrK2tERgYiNGjR+vVw5Rt78HWodrb28POrjWSkqxx+LAAhw4ZIT2dByMjwoABagQFqTBmjBL1XX5kTemMjY3h4uLSoG9sVg0lFot1liSypNXVRjM3F3BxMcV//6vAnDl1J6e0QQQsXSrA1q0CzJtXjvHjr8Lbu+am1UlJSQgJCUFwcDCio6Mxd+5cDBkypN7nWxuUSiXc3Nxw9uxZODg4oGfPnvj111/h4eHBbTNixAhMnDgRM2bMQHx8PIYMGYKsrCy9rn12djbOnj2La9euoVevXggJCcGTJ08wY8YM/Pvf/8bkyZO1N/9nZpUNBbYVZkREBI4cOYIWLVogMDAQY8aM4Trh6XoPa83q7u7OZX0LCwthbm4OW1s7ZGTY4sgRASIjjfDwIQ8MQ+jTR7NWHBioQrt2tV9itVqN2NhYmJubw9nZ2SDnyiaMxGIxioqK0KJFC5SUlMDV1VXnslV2NoNOnUzw7bdyzJxZ3aS+JhABH30kwObNAsydW44JE2onbWpqKiZNmoSdO3fq1eayMTh27Bj+85//cO0xly9fXqk9ZlxcHObMmYOysjIwDIN169Zh6NCh9TrGgQMHcP78eaxYsQJ2dnYIDg7G4MGDMX/+fO3NXhHXUGALww8cOIBDhw5BJBJhzJgxCAwMROvWrcE8bXPCLilVbSTNlvnl5eXh8ePHT/vr2kMi0ZD40CE+7t3TxErdu/+tn3ZxqXy52cyura1tk3T7AzTrr7du3YKZmRnKy8t1liRmZjLo3NkEW7fKERKiH3GJgGXLBNi0SYA5c8oRHFw7adPT0zFp0iT88MMP6FWLO0lzBmu4wMois7Oz8c033yA1NRVGRkaQyWQ4duxY1be9Im5TgIiQnp7OkZiIMHz4cNy4cQMff/xxnX7KutZfNfppO0RFCREZaYRbtzTBo6fn3/ppV1cF7t6N4bLcTQHtel22rpgtSWS13nZ2dpBK7eHjY45t2+SYPr1u4hIBy5cL8L//CTBrlhSTJl2Bl1fNGvCsrCxMmDABW7ZsQb9+/Qx9ms8ELGlVKhWGDx+OkJAQvPnmm4iOjkZKSgpyc3OxYMECXW99RdymBjudDgwMhLm5Ofh8PkaPHo3AwECdGmRd79cmhkgkeqqftuP001ev8kDEoF27cowZo8CUKSL4+pLB9dMymQzR0dG1CkTYksTo6FJMndoba9ZkY+5c41q9nIiAFSsE2LhRgLfekmLKlNpJm5OTgwkTJmDjxo0YNGhQo8+rLkUUoHHLCAsLA8Mw8PHxwa+//tqoY7KkBYAxY8agU6dO+Oabb/R9+yviPgvExsYiISEBb7zxBsRiMX7//XccOHAAxcXFGDVqFIKCgvQud9NVBCGVWmLnzkLcvNkB166JoFIxcHLSSC+DglTo2bPx+mm2yL5Lly41ZtK1kZTEwMfHBOvXZ6N372TOUbJqSSIR8NlnAmzYoB9p8/LyMH78eKxfvx4BAQGNOynop4hKTExEcHAw/vjjD1hZWUEsFjeqwkitVnNLROPGjUO7du2wadMmbjx6FMO8Iu7zRH5+Pg4dOoSDBw9CIpFgxIgRGDt2rN7uEOXl5cjKykJGRgZMTU3h4OAAHs8OZ86YIjKSjz/+0Oin27b9Wz/dr1/99dMsad3d3fWR4wEAEhIYdOtmgp075QgOVnFF93l5eVq6ZHts2GCFr78WIiREiqlTr8DLq2blmEQiwRtvvIHVq1fXO+lTE65evYqwsDCcPHkSALBmzRoAwMcff8xts3TpUri5uWH27NkNOoY2UbUxYsQIODs7Y+vWrbVupwP/TMljc4GNjQ1mz56N2bNno7CwEIcPH8aqVauQmZmJoUOHYty4cfDw8Ki1ljg/Px/du3eHSCSCWCxGTk4sunYlvPaaHYyN7XHunBkiI43w0098fPedADY2hDFjVAgMVGLQIDXqqq9g9dNdu3atlxSzqgBDu+heU5JYgI8/rsCuXUKMGyfGpEnR8PT0qvEYBQUFmDBhAlatWmUw0gL6KaIePnwIAOjXrx9UKhXCwsIwfPhwvfavTcbMzEwwDAMHBwfcu3cPAwYMwLJly6ptZyi8Iu4zgJWVFUJCQhASEoKSkhJERUVh3bp1SE5O1unuwXaY9/L6O+vKWsqwSqjMzPvo2FGF1attsWmTPf78syUOHTJCeLgRdu7kw9Lyb/30kCHV9dOsnU1DjAJqE2AYGfHx44+O2LVLgDfflGLq1GiYm7dAXFyczpLEwsJCTJgwAZ988klTVdvUCqVSicTERJw/fx6ZmZkYOHAgYmNj9QoZ2HNYtmwZZ5LXpUsXrFmzBp6engD0nh7XG82CuHUlEORyOaZPn45bt27B2toa+/btM5jA/FnD3NwcU6ZMwZQpUzh3j2+//ZZz93B2dkZ2djaWLFmiU2ivrYRSKBRPSZyAtm0V+PRTG2zYYI+//rLAoUNGOHbMCL/+ykeLFoThwzXT6aFDVQDKEBsbW+mLoT5gn7i6HiKrVwuwdq0AU6fKEBJyFR4eXrCwsKhUkvjw4UOUlZUhNTUVe/fuxZIlSyppdw0FfRRRjo6O6N27NwQCAZydneHm5obExET07Nmz1n2zhPzuu++QkJCAvXv3Ytq0aVwtMosm67tLRPX5MTiUSiV17NiRkpOTSS6Xk7e3N92/f7/SNlu2bKF58+YREdHevXspODi4KYbyXFFeXk4rV64kW1tb6tGjBy1YsIBOnTpFJSUl9OTJkzp/ioqKKCkpia5cuUJnz56lmJgYSk/PochIKc2cqSBbWzUBRMbGaurXL4+2bCmmrKy696vr588/ywkg2r9fVun3y5fLCSCaMkVGp0+fpezsbJ3vLysro4sXL1Lv3r3J2dmZxo0bR3l5eQa/phUVFeTs7EwpKSncZ+vevXuVtjl+/DhNnz6diIgkEgk5OjpSfn5+jftMSkqq9Hr//v10+/Zt+vjjj7nPpUwmoz///LOhw9aLi8+duFeuXKGhQ4dyr1evXk2rV6+utM3QoUPpypUrRKS5GdbW1qRWq5tiOM8VYWFhlJeXRzKZjKKioigkJIQ8PDxozpw5dOzYMSouLtaLWCUlJZSSkkJXr16lM2fO0J07dygjI5sOHMinoKB0atNGSQCRQKCmoUOVtHWrjNLS9CfuxYtSAogOHJByv1uxQkPayZM1pM3Kyqrx/Xl5eTRw4ED6+eefiYgoLi6OFApFk1zTo0ePUqdOnahjx470xRdfEBHRihUrKDIykoiI1Go1LVq0iNzd3cnT05P27t1b474OHDhApqam9O6773K/27p1K5mZmdHEiRO5361atYreffddqqioaMiQXwzihoeH06xZs7jXu3fvprfffrvSNh4eHpSRkcG97tixI0kkkqYYTrODQqGgkydP0pw5c8jDw4NmzpxJkZGRVFhYqDeJHz16RJcuXaLDhw/T9evXKT09k/74o5zee09BTk4qAoh4PDUNGqSkjRvllJRU+z7PndMQ9+BBDXE/+0xD2okTZXTqVO2klUgkNHjwYNq5c+fzvrT1wuPHj2ngwIG0bds2CgwMpAULFnB/W7RoEXl5edGFCxfoo48+ot69e5NYLG7oofTiYhNNwF/BUBAIBBg6dCh++OEHREdHY+rUqTh+/Dj69++P0NBQnDhxAnK5vMb3GxkZQSgUQqFQwN/fH61bt0ZeXg7U6iuYNu0uLlzIwp9/luODD5TIzWWwaJEQnTqZ4PXXRfj2Wz7S06uvTmhnldev52PlSiGCg+WYNesyunbtXGNiRyqVYvLkyZg8eTJmzJhhiMvzzGBlZYXt27djxowZ2LRpE5KTkzF37lwAwDfffIOJEyciKioKJSUliIyMhK2tbdMOSF+G06upcrOCUqmkCxcu0MKFC8nT05MmTZpEv/32G+Xn51d6wmVmZtLZs2epoKCg0u9LS0spIyODbt68SWfOnKG//vqLUlNT6caNMlqxQk5eXponMUDk56eksDA5RUeX05MnT+jUKc0TNyREQQDRhAmaJ21mZmaNT9qCggIaPnw4bdmy5aW4d5mZmTRixAiaOXMmERHFxMTQo0ePDLHrF2OqrE8CYfPmzZWSUxMmTNC5r+PHj5Obmxu5uLjQmjVrqv19w4YN5O7uTl5eXhQQEECpqamGP6HnAJVKRZcvX6bFixeTt7c3jR8/nvbs2UN79uyh7du3VyOtrmRRVlYW3b59m86cOUPXrl2jR48eUXR0KX3+uZx69FByJO7aVUVTp1Zwr8ePl9Hp039QRkZGjfsvLCyk0aNH03//+1+Dkrau+80iIiKCANCNGzcMdmwiopycHJo8eTK5urqSu7u7ocK3F4O4RHUnEKRSKY0fP55cXFyoZ8+elJycXG0f+mSn//jjD3ry5AkRaZIKL2N2WqVS0Y0bN2jChAlkbW1NQUFBtHPnTsrJydErJi4rK6Ps7GyKjo6ms2fP0tWrVyk5OZnu3y+hdevk1K+fkgBNhtrDo4JOnz5XK2mLioooMDCQ1q9fb1DS6nO/iYhKSkpowIAB1Lt3b4MTl4ho7dq15ODgQImJiYba5YtDXENAnym3Nm7fvk19+/Z9FkN75lAoFPTGG29Qbm4uRUdH0/Lly8nPz49GjRpFP/zwQ63Jo6okzs3NpZiYGDp79ixdvnyZEhMT6dy5EurWrYJ2775M6enptZL2jTfeoNWrVxt8eqzv/X7vvfcoKiqKBg0aZHDi5ufn05AhQyguLs6Qu9WLi81CgGEI6CNv08aOHTswYsSIZzG0Zw6BQICIiAgAgL29PXx8fPD5559z7h6BgYF6uXtoNyxzdXXl2rmo1TewcmU5HBwcalRdKZVKhIaGwtvbGx999JHBegqx0Od+3759GxkZGRg1ahTWr19v0OMDgLW1NaKiohrd6b4heGmIWx/s2bMHN2/exIULF573UJ4ZGIaBh4cHPDw88Omnn3LuHsHBwTAzM6vT3QPQNCxzdHSERCJBly5doFAouK6DrJOkSCSCSqXC22+/DRcXF6xYscLgpNUHarUaixcvxq5du5r0OM+DtMBLRFx9WxyeOXMGX375JS5cuNAk7TleBDAMg86dO2P58uVYtmwZ5+4xdepUCIXCau4eLBQKBaKjo+Hi4gIbGxsAgJOTE2QyGWcHu2nTJhQWFsLNzQ2rVq1qMtLWdb9ZL+vXXnsNAJCbm4uxY8fi8OHDDWnE1fyg75yaDBTjVlRU0Pnz5xuqKql1v3Vlp2/fvk0dO3akhw8fGvTYLwvUajWlpqbSN998QwMGDKD+/fvT2rVr6cGDB5SSkkLh4eGUmppaY0xbWlpKb731Fg0bNoxef/11Wrp0aZONVZ/7rQ19YtzY2Fjq06cP9/rWrVsUEBBgsDHrieaZnEpKSqKAgAAyNTWl8PBwQ+ySQ13Z6SFDhpCdnR35+PiQj48PjRkzRud+nvcyQ3OAWq2mrKws2rRpE/Xv359at25N8+fPp3v37lFZWZlO0i5YsIDefvttUqlURETcv02Fuu63NvQhrkqlInt7e1Iqldx7bt26ZfiB147mSdySkhKSSqXUq1cvOnfuHBFpLpharW7yG60PmssyQ3OBWq2m4cOH0+7du2nbtm30r3/9i3r16kVhYWEUHR1NZWVlVFpaSv/5z39o3rx5zeIeNgYBAQEUExNDERERXPHBM0bzzCq3bNkS169fR0VFBRd/sKVPzyOJURX69JEBgBUrVuDDDz9skmxlcwLDMNi3bx9XBB8aGoqCggIcOnQIH3/8MSQSCUQiEVxdXbFjx46mK2N7RvD398fly5exdetWzgC9OeK5XOXjx4+jW7duADQF3X/99ReWLFmCixcvPo/hVIKuZYasrKxK22gvM/wTUNW5wtraGrNmzcKxY8dw6tQpvPbaa9ixY4fBCsZPnDiBzp07w9XVFWvXrq3292+++QZdu3aFt7c3hgwZgrS0NIMcF9AQ95NPPsG4ceOadc/dZ0pctVoNqVSKc+fO4Y033gCgcQ84f/48evToga+++gqfffbZsxxSvcEuM2zYsOF5D6VZwMrKCl988YXBSMsuJR0/fhxxcXHYu3cv4uLiKm3j5+eHmzdv4u7duxg/fjyWLl1qkGMDQJcuXSASifDhhx8abJ9NgWdKXB6Ph9zcXBQXF2PYsGFQKBTYt28fwsPDoVQqObNo1gcI0NxIopo96kpLSw06xvosMzg5OeHatWsYO3Ysbt68adBx/FOhHaoIhUIuVNHG4MGDOXcQf39/ZGZmGuz4//vf/7BmzZom6axoSDwz4j548AChoaFYunQp/Pz8wOPxcO7cObRt2xb79+/H5cuXMX36dJw6dQouLi6QSqUANGVpbOyrUqmgelpTplKp8PPPP2P8+PHo3bs3fvjhB8hkskaPs2fPnkhMTMSjR4+gUCjw22+/VbJVsbCwQH5+PlJTU5Gamgp/f/9a1wbrmvYBGl/frl27wsPDA1OmTGn0ObzI0CdU0YahFHDJycno0qULpFIpQkJCGr2/psYzS045ODjAzc0NERERuHPnDvr37w8AcHd3r2RjmZOTAyMjI3z//ffYvXs3Ro8ejeDgYHh6elaajiUkJGDjxo347rvvYGxsjPfffx/u7u4YMGBAo8bJ5/OxefNmDBs2jOsjw6qN2D4y+oKd9mn7+o4dO7aar++aNWtw+fJlztf3FfSDIRVwLi4uePDggQFG9Yygb/qZDFhkkJubS3/99RcpFAp68803KTg4mCIjI+n+/fskl8u57W7evEnLly8nf39/8vf3pyVLlnCVQSkpKTRs2DBu2x49etCBAwcMNUSDQB8h/JIlS+jHH3981kNrttC3eOD06dPUpUuXJvGqes7Qi4vPJatsb2+Pnj17QiAQYMuWLRg2bBi2b9+OiIiISt3Iu3fvji+++AJXr17F9u3bIRAIkJSUBACIjo4GEUEoFGLEiBHo3r27wdsvNhb6TPsePnyIhw8fol+/fvD392/WSxDPAnWFKgBw584dzJs3D4cPH25U14EXGvoynJ5BWV9VGaRareZULNrYv38/zZw5k+7evcuVsLGues3JXUEfP61Ro0ZRUFAQKRQKSklJIUdHRyosLHzWQ21WMJQC7gVF8xRg1Abtpy2gWfxn41q1Wg2GYaBSqZCcnAwbGxuuk7eNjQ331GoOIg4WTenr+zJj5MiR1czRV61axf3/zJkzz3pIzQ4vjMyFVeTw+Xy0a9cOCQkJADRJKplM1qwIy0KfaV9QUBDOnz8PQNNv6OHDh5xqqyrqylCnp6dj8ODB8PPzg7e3t67eq6/wskDfRzM1IweMjIwMmjdvHrm4uNCgQYNo3bp19PjxYyJqXlNlIsP5+uqjoZ4zZw5t3bqViIju379PHTp0aLoTe4WmQvMsMjA04uPjn/cQngn0ybbOnTuX1q5dy22vXaL2vFFXxZVMJqPg4GBycXGhXr16Gcox8UVE880qGxJdunSpVVn1skCfDHVYWBj27NkDR0dHjBw5Et9+++2zHqZO6CNj3LFjB6ysrJCUlIRFixY1e8nh80Z9++O+wnMCwzDjAQwnotlPX08D0JuI3tHaZjE093QDwzB9AOwA4ElE6ucy6L/H1QdAGBENe/r6YwAgojVa25x8us1VhmH4AHIB2NKrD6hOvPBP3H8QsgC003rt+PR32pgFYD8AENFVAMYAbJ7J6GqHA4AMrdeZT3+ncxsiUgIoBqDbxe4VXhH3BcINAJ0YhnFmGEYIYBKAw1W2SQcwBAAYhnGHhrgSvMJLh1fErQcYhjnHMMy/nv7/C4ZhnlkQ+fQp9A6AkwDiAewnovsMw6xiGIZdY3ofwByGYWIA7AUwQ9dUk2GY/2MYRswwzD1dx2I02MQwTBLDMHcZhunWyOHrM1vgtnk6VbYAUNDI47680DeL9eqHAGBzs8p/AAABS0lEQVQggPMApgI4CsDoeY+pEefRDcC9Gv4+EsBxAAwAfwDXG3k8PoAUAM4AhABiAHhU2eZtAN89/f8kaL6Ynvu1aq4/zUo51dxBRBcZjdJjMYDXiEj1vMfUEDw9D6daNgkEsJs0LLrGMIwlwzBtiCingcdTMgzDzhaMAPwfPZ0tALhJRIehSaT9zDBMEoDH0JD3FWrAK+LWAwzDeAFoA6CAiAxbwd+8UFMyqUHEBQAiOgbgWJXffar1fxmACQ3d/z8Nr2JcPcEwTBsAv0DzNCpjGGb4cx7SK/yD8Yq4eoBhGFMABwG8T0TxAD4H0LzNsRoHfZJJr/Ac8Yq4eoCIyomoDxGdfvr6IhH1ed7jakIcBjD9aXbZH0BxQ+PbV2gavIpx/4FgGGYvgNcA2DAMkwnN7EEAAET0HTSx6EgASQDKAbz1fEb6CjXh/wH8oAJ3uYT4GQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "mesh = oc.Mesh(p1=(0, 0, 0), p2=(10e-9, 10e-9, 10e-9), cell=(1e-9, 1e-9, 1e-9))\n", "\n", "# Plot the mesh.\n", "%matplotlib inline\n", "mesh" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now we define the system object and its Hamiltonian." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "system = oc.System(name=\"system\")\n", "\n", "A = 1e-11\n", "H = (0.1/oc.mu0, 0, 0)\n", "K1 = 1e3\n", "u = (1, 1, 1)\n", "system.hamiltonian = oc.Exchange(A) + \\\n", " oc.Demag() + \\\n", " oc.Zeeman(H) + \\\n", " oc.UniaxialAnisotropy(K1=K1, u=u)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We will now intialise the system in $(0, 0, 1)$ direction and relax the magnetisation." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "Ms = 8e5\n", "system.m = df.Field(mesh, value=(0,0,1), norm=Ms)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Effective field\n", "\n", "Total effective field is:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2018/3/26 16:32: Calling OOMMF (system/system.mif) ... [0.7s]\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUUAAAERCAYAAADsa7W8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsvXmcHHd55/9+qvru6blPzejWyDptWZYl29gY8IkXsAOENbsJLMsPhxfHJptkf0CSBZaEZGGT8IsXYmJ+eLmxjTkswPjA921ZlqzTuq+R5r5nevqq/u4fVS01Qta0pKexeqi3X/XyTHX1RzXdVU893+v5iDEGHx8fHx8X640+AR8fH5/zCT8o+vj4+BThB0UfHx+fIvyg6OPj41OEHxR9fHx8ivCDoo+Pj08RFRkUReQuEekTkW1Kel8Wke0islNEbhcR0dD18fGpPCoyKALfAm7UEBKRK4A3ARcCK4BLgas1tH18fCqPigyKxpingKHifSKyUEQeFJGNIvK0iCwpVQ6IACEgDASBXtUT9vHxqRgqMii+DncCnzTGXAL8JfCvpbzJGPM88DjQ7W0PGWN2lu0sfXx8zmsCb/QJaCAiVcAVwI+KugPD3mvvBr5wircdNcbcICKLgKVAh7f/ERG5yhjzdJlP28fH5zxkRgRF3Ix3xBiz6uQXjDE/AX5ymvf+AfCCMWYCQER+BVwO+EHRx+f3kBnRfDbGjAEHROQPAcTlohLffhi4WkQCIhLEHWTxm88+Pr+nlC0oTjdtxgtct4vIXhHZIiKrz0D7h8DzwAUi0iUiHwb+I/BhEXkV2A7cXKLcfcA+YCvwKvCqMebnpZ6Lj4/PzELKVTpMRN4MTADfMcasOMXrNwGfBG4C1gH/YoxZV5aT8fHx8SmRsmWKp5o2cxI34wZMY4x5AagVkbZynY+Pj49PKbyRAy3twJGi37u8fd0nHygitwG3AdhW8JJYtEn9ZCTnqGsCUAZdk8+ra7rCFVZwuEwLj8QqQ64QsPU1AVMm3fHJYwPGmFPeaDe8NW4Gh0q7rjduST9kjFFZaPG7oiJGn40xd+LOQ6S6qt2sW/En6v9GYGBcXRPADI+qa+YnJtU1AUw2UxbdciHBUFl0raq4uqbU1ahrAuQaE2XR/fWLnzv0eq8NDjm89NCcknTstj2Naif1O+KNHH0+Cswu+r3D2+ejSN44aPcbO8bBMTlVTWMMWaMflPMmr/73/75jgHyJ/02HiPypiGzzag/8mbfv8yJyVEQ2e9tNRcd/xhuc3SUiNxTtv9Hbt1dEPl20f76IvOjtv0dEpn2SvpFBcT3wAW8U+jJg1BjzW03n85Wsk2JH/6/JOEk1zb3JjWydeFxND2CX8wqbc0+qah5hLy/ya1XNQXp4kvXkjV7XwKQZ5dHsPSTRawVMOMM8OfID0nm9771/cj8HRzaq6ZUbgyFrnJK20yEiK4CPAGuBi4B3eIspAL5ijFnlbQ94xy8DbgWW49Y++FcRsUXEBr4GvB1YBrzfOxbgS57WImAY+PB0f185p+T81rQZEfmoiHzUO+QBYD+wF/gG8LFynQu4mcih7ufIZHUu5u6JnRwd34Elej0QQ7lubAmq6QGkTJKQRFQ18zhY6PZl2V5PTh69PtgwMcCQMnoBLGJVkc4nGcn1qWkmcyPsHX4WJ59V0esf3sXo+JHpDzwHlDLFpcCLxpikMSYHPAm8+zTH3wzcbYxJG2MO4MaOtd621xiz3xiTAe4GbvaqXb0Nd9odwLeBW6Y7qXKOPr/fGNNmjAkaYzqMMd80xnzdGPN173VjjPm4MWahMWalMeblcp0LwNG+Dew5/DATUzq1HrrGttFWtYSApdOvlTd5xnL91ASaVfQKpJgkQkxV08HBVg6KhSDroNcsD0iQAEFSRq8PNiBBEnYDIzm9miFtVUvJG4eeyd0qesf6N7Flz72ksxMqeidjMDimtA1oFJGXi7bbiqS2AVeJSIOIxHCn5xW61D7hzV++S0TqvH2vNzj7evsbcFe65U7af1pmxIqW6RgdP8KuQw+yaPY11FfPP2e9sXQfY5leOqp/a/rlWTPhDOGQo1Y7KJokEdEdOChPplgIirqj9RGJqwZFgJpAs2pQDNlRWuKddI1tVdFbvuAWbDvItj0/Ij9NE/ZsyWNK2oABY8yaou3OgoZXeOVLwMPAg8BmwAHuABYCq3Bno/xTWf6I12HGB8VMdoIte++lsXYxc9uuVNHsGt9KVbCBmrDetMqRXB9BCROz9EYpcyZLjgwR0c0UyxEUC3qazWeACHFS6DWfAWoDLYzlBlQDTkfiQoZTXUxmTje1tzQCgQgXdt7K2OQx9h3R7fsFd6DFwZS0TavltiAvMca8GbfPb7cxptcY4xhj8rhda2u9w19vcPb19g/izn8OnLT/tMzooJg3Dlv33odtBVm+4BY0Cmo7+Rzd4ztpr16poldgNNdLTaBZVbMQDLQzxXI0nwt9iprNZ4CoxJlSzhRrA83kcRhzBtU066OziQZq6BrXyRarYs0sW3Azh7qfo3dwu4pmMWeQKZ4WEWn2/j8Htz/xByct4vgD3GY2uIOzt4pIWETmA53AS8AGoNMbaQ7hDsasN+60g8eB93rv/yBw/3TnVBHzFM+WfUceY3Sii7XLbyMQ0Bls6JvcSy6foT2xbPqDz4CRXB+zwp2qmoVmo3afYqU1nwfyx1Q1o1aCkEQZzfWpdXeICB3VKzk0+gqd9Vdiybl/vi0NKxid6GLH/p8RjzZRFdM5VwNk9aY5/VhEGoAs8HFjzIiI/G8RWeX9UweBPwEwxmwXkXuBHUDOO94BEJFPAA8BNnCXMabwJPgUcLeI/B2wCfjmdCc0Y4Ni39BODnU/w4pF71W7GMBtOjfHFxKy9QJNJp8imR+jNtCipgluf2KQMLbiCDm4gSuI7sTp8jWfY6RIYoxRy8JFhNpACyO5Xuai16/cnljOnqFn6U/upyWu84BcNOc6xia72bLnbrXkwJTYNC5Jy5irTrHvj09z/BeBL55i/wO4M1pO3r+fE83vkphxzefJqQEmpwbYvu+nzG69jNaGlWrayewog1OH6EjoaQKMetM7agK6yxdTZlK9PxHK1KcoFoKoZ4pRiWPIkyGlqlurPNgCEAkkaIrNVxtwAbDEZmXnH+I4Wbbv/xnGGJKpc2z2G3BK3CqRGRUUU+kRtu65ly177qEq1kLn7OvVtHsmdnF0fBthu4rG2Dw13aQzxkiulyq7jsD0k+1LxjE5UiSJoNefaIzBGINThqAIbr+ig+4KnEJ/6pSZVF2FUxtoIZWfZMqZIOmMqel2VK+kP3mA8XQ/g1OHVTTDwSou7HwfAyO72bbvPnYfevCc9NwVLaVtlciMCorH+jcxMdXL5FQfzfVLcPJ6y8ZeG3yCgyMbqA43MzTVpaa7bfJJjqb3EJAQXenX1HQ3555iJN9HjgyHHD3dLTxPminGGeGg2aWimTUZdppXMBh6OMwhdObrpUySLmcPAHucTRzN79PRzU8w6vQDsGH8lwzndBZi5fIZ8sbBEpsXjn6fkZReX6hth6hLzKV3cBsDI3tIpc9lTb7glLhVIjMmKBqT51j/JgBELIzJE7DDavpZJ4VjcgwkD6lN2AawJUjaTDKS6yWA3mqWoISYYpJh00cInc+h0CeXZooheokqZaHuuU7gkGOIXqrQmZYUkRhDxm3ijph+qkRHNyxx+jNuvYSp/BhBpRVDASvE8FQXeZPDMTkyzpSKLkDADmOOL6E0HBvYdNZa7kCLlLRVIjMmKA6N7ieVGSUWaeTS5R9h3qyrENH585x8Fse4S7CWNV1DTaRVRRcgKG7Aqg000xJaoKZbJbUAhIjQYpVW0aQU6nAHrWxsGtH7HFq8aWaCUIteYZU59gXHfy58JueKiHBB7LLjvxe+Qw2WNL7l+PzXrGJQjIRrWb30gyyafR0iNsf6NhUFyTPDnafoZ4rnPcf6X6G9eQ3rVvwJ1fFZqtrZvNtJ355Yrj7IUsgyLohepjpHsRAAOuxOlekdBeq9oNjILNVR7WbaEYRq6gko6jZJBxFihIiorgGvDjTSHloMQNDSC4qWBLi49V2ErKhqsRFwW1DzZl3Jpcs/gm0FGRrdf9ZaeSMlbZXIjAiK+XyOWU2rWTr/ndi2fo29jDNFItTEssZrVQMXQMgK0xKaT21QdzpOQmoRLDqsRdMffAbESbjZ53FHWB2CEqKeluNBVwtLLGbbi9WazsUsiq3BJqCaKYI7Cn1RyzvJ5tOqugWq422sXXHbWRfpnemZ4oyYp2hZARpqdW/+Ygx5Lm69GdvSrWADELbitIX0zz1CnA5rIWGJquqKCM2mnQbFpnOBVmYTRvd8AdqtheSMTgWaYiJWnHnRi9SDIkBDbA7z8peo6xaw7RANNQvP6r0GwZkZ+dQpmRFBsdzUhPUDQIG20EL1ydXgBq+Fdqkur2fGApZhKzbJCzQxCynDzRaUMHPtpeq6APMjF6p2TxTTVrWkLLoaVGrTuBT8oPgGU46AWCCoOO+xmFAZMiNwy3KVi3J9FuX8/s5XDELGlOdBcD7w+/eN+vj4nBPu5G2/+XzeYGwhW6ufqUiZliRZQf2P2K6uUtcEIKPf71ZWQuXJLE1E//pyqvX7SoGy3AulUKmDKKVQcUHRx8fnjcUYwTF+pujj4+NznLyfKfpUKrl8FlsCqvMrM84UjskRDeh5Djsmx1RunHigRm0lEkDe5AAp2wjx7yPuQMvMDR0zNwcuM+PJHnX/i01dP6FvfI+uZt/9vDb0hKrmocktbBr+rdJ158Ropo9n+r9PTtH7OW8cHj70VQamDqppAmw99nO6x3aoaqYz46QyetV2yklhoKWUrRKpzLM+Q9Jp3Yst56R5acf/T/fAFjXNrJOib2K3apYEMJUbJxzQHZhxTE7V2hVOTG3RLO9liU3YjjGV0/N9Bkg7SQYnD6pqvrLruxzqflZVM+9kyWZ0rRgKOEZK2iqRGR8Up5KDvPTUlxke1MvA+oZ3gjE01+lNrh1L9QBQE9EzwzLGkHLGidh6zVxwm6S2cs9LOYIiQCRQTSqn+1CsjcxiZGpa/6Mzoq3xQroHt5LP67U+drz6A3bv+JmaXoHCipZStkqkMs+6RIwx7N25nli8mdo6vQo03QOv0lR3AcGA3jSLkamjxIJ1hAJ6lbKz+SnyxlHt+wPXXEq7j65cQTFqJ9QzxZpoO5OZQbKOXjXv1oYLyeaSDIzqPbzb517JQO9Whvp16l4WkzdWSVslUplnXSKDfdsZGthF57JbEEvnJk5lRhkaO0Bb4yoVvQKjU8eoiepW9ykEg4hyUMybnPpKDttbzeIor1GOBBKkHN2gWOt9T6NTekVgI6FqGmoW0j2wWU2zrmEhzW0Xs2fn/TiO3ufqFoTwM8WKYny0i1wuzd7Xfk5r+6VU1+rVE+we2EIoEKOh+uwW058KYwwjKf2gmHImACFsa/cpOuqZYkEvr50pBqqZUm4+B+0I8VCDamVsgLbGVfSP7CaT1SsZtvCCm8hmkxw58AQTY8fI58/98zUIWWOXtFUiM3JcvevQM0yO9+A4WRYsvlFFM5UZxXEydA++SmvDSiylzNMYw1R2lKwzRW1EOSjmxojYcSzlwZu8yWFbuis0ytZ8DiRIO5PHy/xrURM90a+o5RTYXLeE16wAvUNbiUebqK8+9y6fUDjB/M4b2PfaL+jv2cLi5e+hpm7eOWkaw4yevD3j/jJjDCOD+5ic6EHEortrg4oRUjozzoYd32Ryqp9IuIZkakjhbOHQ8MscGHzeHSkNxMkp+cr0TO4mmR0h4gUFjc/AGMNYduD46POUUrM0m0+7mtjk8hnGswMquhlnirDtWiYMTh1mNN2jopvKjlMTaWV06hjbex4k4+iM8I5NHKU63sG+o0+w+/DDKpqZ9ATDA7sxxiE52c/w4F4FVSFf4laJzLigmJzsI5Nxb9bqmtm0z71C5SmezaWOl4c/2reRUFDHn0QQukZfJW8cNh65V22qy8DUIQ6Nb2Y8M8DWgYdVPgMRYdvIY4xmejmafI0DE2fv81GMJTZP932fPHm2jT7OaKZPRTebT/Fy708A2NinNwo7mRlkZ++vyeXTdI1sJmjrZM1ZJ8XQ2D6yuSTZnI4VQShcxez5byEQdAfwNIKiwc0US9kqkco869MwMug6trW2r2H5qj9Sq8Sd80YaLSvIhYvep2aKVTzavLj5LWpN3aqQ63PimCxzEnp1FRvCHeRxyJk0bRGd4ri2BGgIdwCGvMlRH9ap6h0P1lEbPtElURVsUNFtiM+jKe42bW0rpNYsb65bwtzWK4AT15sGNXVzuXjtRwlHahkfPYyTO/eK3v5ASwUxMrSPOQvexuLl71EbcYYTF+mSuTdRFdOzDihkGfWxuTTG9QZvEkE3KMYCNTRF56vpNoTcgBW24tSG9OZUtkbdABu1E8QC1Wq6c6svBiAWqFWtnL64+a0IQkgpSyywqONaaqpmk3PSZ20sdSpiVc1cvO5jxOLNjAwfOCctQ2n+LJVaiHZGBUVjDE2tK5nfeb26l0rOSdHWcCGzGi9W1S3cVBc0v1X1nBNepjgnsUpVty7UhmDRGl2kqtsYnkNAwmpZYoGGyByqgvXHPw8tqsKNdNSuUms6F7AsmwsX/iHBQJSco+vREo5Us2rtR92RknPAtTgNlLRNh4j8qYhsE5HtIvJn3r56EXlERPZ4/6/z9ouI3C4ie0Vki4isLtL5oHf8HhH5YNH+S0Rkq/ee26WEi3ZGBUURoblNd/5ggXCwmiXz3qEebIN2jFnVy6lWtE0FN9jGA3W0J5ar6tpWkNpQK23RTlVdS2xaowuPZ6JaiAhzqi+mKqgbFAEWNl5JNKhjm1pMJFzDigXvVg+KAIFghIbmc7VmKM20arqaiyKyAvgIsBa4CHiHiCwCPg08aozpBB71fgd4O9DpbbcBd3g69cDngHWe1ucKgdQ75iNF75t2OsqMCorlpLVhpVo/YjEhO8qipjer6wJcUP9mVfvNArNjy6lRdh8EaI12qmeKAO3xpdRH9HXDgTidZfruGms7iYT0HQg1MKitaFkKvGiMSRpjcsCTwLuBm4Fve8d8G7jF+/lm4DvG5QWgVkTagBuAR4wxQ8aYYeAR4EbvtWpjzAvGnX7xnSKt18UPiiWinSEWsK0g0WB5Lv7mmN7SxmLaoovL8nk0hDqI2Dqj+sXYVpCG6Gx1XYB4qL4sulC+a04DJYvTbcBVItIgIjHgJmA20GKM6faO6QEKT+B24EjR+7u8fafb33WK/aelrEFRRG4UkV1ee/7Tp3h9jog8LiKbvD6Cm8p5Pj46lOtmPZ+DgM8JjJEzyRQbReTlou22EzpmJ/Al4GHgQWAz4Pzmv2UMbnL6O6NsK1pExAa+BlyHG6E3iMh6Y0xxIbq/Ae41xtwhIsuAB4B55TonHx+fc8cdaCl5ZseAMWbN62oZ803gmwAi8ve4saJXRNqMMd1eE7gwcfUobiZZoMPbdxR4y0n7n/D2d5zi+NNSzkxxLbDXGLPfGJMB7sbtEyjGAIX5FzWA7mJSHx+fMiBqk7dFpNn7/xzc/sQfAOuBwgjyB4H7vZ/XAx/wRqEvA0a9ZvZDwPUiUucNsFwPPOS9NiYil3mjzh8o0npdyrn2+VTt/HUnHfN54GER+SQQB649lZCXct8GEIrVMdWof9q5WHmeD3a9/kCHndabv1aMOL/TVso5Y+zyNLedsP614ETLc31ly3Tdng53oEXts/+xiDQAWeDjxpgREfmfwL0i8mHgEPA+79gHcPsd9wJJ4EMAxpghEflbYIN33BeMMYV1uB8DvgVEgV9522l5owtCvB/4ljHmn0TkcuC7IrLCnDRr1RhzJ3AnQLxhdmXduT4+MxCt1SrGmKtOsW8QuOYU+w3w8dfRuQu46xT7XwZWnMk5lTMovl77v5gP480bMsY8LyIRoJETfQg+Pj7nGYUVLTOVcubeG4BOEZkvIiHgVtw+gWIO4z0RRGQpEAH6y3hOPj4+CvjGVWeBNxnzE7idoDtxR5m3i8gXRORd3mF/AXxERF4Ffgj8J6NR46oCyecdjKI/R4H+/m2kUiOqmoNDe+jp1asQDZDJTHDg0GNq1WEKDA/vY2Kie/oDzxDNStaVhjGQzVslbZVIWc/aGPOAMWaxMWahMeaL3r7PGmPWez/vMMa8yRhzkTFmlTFGp4hcmRk8tp2RPl0r0uGBXTz768+rVEYuYIxh+/Z7GBs7rKYJMDC0k+5enbJhBbLZJAcOPUo2q+s+t3//I/T06Abwo4eeY+MzX1HVzGVTHN37FNl0edz3NHGbz75HS0XiZFPse/YHZJK6mdKRXb9mqGenqubEeDfhaC2WpdfNm8tNYYxDKKTs0ZLPYSueJ3D87847upW3w+HE8fqaWkRjDaSmhsjl9Mp7iWVzePuDDPfqmkwNHtxE947HVTVBbUXLecmMDooD+19muGs7djCipplJjTM5cpS6lgvUNAGS473Eq3SLQhSCgX5QdBDtoGh7xlV53WZpKFSt7vsdT7gl0ybHdSp5A9h2kJqmhQz3vqamCZBNTXBs26M4Wb0AXpiS45cOqyDSE4MkR3ro3f0cjQsuUQ2KI327sawA1Y2664onx3uIJ3SLLKTT5QqKOTWPmgLHM0XF7gOAcLiaTEY3KIbCCYLBOJPjun2VdS1LGOndrVpHsXHBGozJM3BgIyNHtVo3fvO54kiND7Lzka+RGuujqmEumSmdm+LIrkcZ6tlJddNCbFunYGk2m8RxsiSTA8czEC0ymXFsO0QgoDuB3ORzqs184HgB2Lx6ppg4/nDQQkSIV7cxOd5DNqPnvFfbsoRcNknfoY0MdL2qoplJjlDTupiuVx/k0MZpF3OUjO/RUmFkkqM4GXcUs+e1J7GDOkGh7/BGBo9uIT05RNeux1Q0x0e72PLSnWDyJCd61ZpkXV3PMzJygGAwzuDgbpXRUmMMff1bcfIZBGF4ZL/CmbqmYMmpQUCYSg0xPHJulaELjI8fwxgHx0lz5MgzTEzofLZ9xzZjWTZD/bvY/sq3p39DCeQySQ5t/yUiNns3/YhMSudBnpkcYbhrO05mimxyVMnADLJ5u6StEpmRQTHrZYaBcJzOqz+ErZQpiefFkU1P0Dz3UhXNSKSW8VF3NWT3kReJxnWKoebzObq7XyaVGubw4SdVMlsR4VjPRoZH9nOs52UGh3QGBYKBKK+8+g3AsHvvz9X8SSzLZteunwKwd++viETqpnlHqboBhvp3kU6NqFleBEIxqhsXYIw3LUvJq6e2fSmzV/87APJO9niycC74dgQVSCY5ilg2nW/+T4TjOjcCgHgX6ryV7yQU0emnC0dOVG6eveCtas3SROKEYVN7+8lLzs+e+roTZlXNjWe0eup1sawATQ2FatBCXY2Op0w83kJtrasVizWodSM0tCw73tURjujVwmxbcAWJujmAbhm11iVX07jALVSjNRPDbz5XGNmpMeZe+gckmvUMmwDEsqhpWkTznEvUNO1AiEAwRjhSS2vH61ZYOmMSCbeWZjAYp7FxmZpufa0bFCPhWhIJvWrWzU0rAahOdBAI6A2MtbdfBkBV1axpjiwdEYt5ndcBukFRxGLR6vciYh9/AOvoCvPWvpeqxnlkkqPnrDfTR5/f6IIQZaFuzkoa5+sFrgJ2IMzCVe9WL4YaidbSNnud6uBFIBAhGm2gqWmFqm483kIwGKepaYWuIVbtQgKBKPV1eo6GAI2NywiFEr+ROWtQ37SUquqO38j0NYhVt9JxwdtAOcuy7ACdb/4gyRGd6nyVOrJcCjMyKJYjIALMWXo90Sp9A6Sauvm0tOufcyLRzqxZa1U1RYT62kVqTecClmXT1LiculrdoGhZNm1ta1QzRfCyr87rTvQBKtKx+K1Mjk5bC/WMCUYT1ETPfX6tMULOD4o+ADXKcxMLzF10nfoUF4C5c68mGtXrUy0wq20N1YpN5wJtLaupTkxroXHGzJq1Vm0KVTF1jYtVjOVPxrIDJOrnqutqUqlN41Lwg+J5QEBxcnkxVVW68x4L1NWW5+FQW1OeQBBR7PcrRkTK9t2dzygXmT3v8IOij4/PGeMHRR8fHx+PmV5k1g+KPj4+Z0ylzkEshYoLivkAJJv0R77sTHm+ZCunv9RJdGsmnNCtsPK+5UpWTBnuinyZ7jQn9LsPTsZArkILyJZCxQVFHx+fNx6/+ezj4+Pj4fcp+vj4+JyE8YOij4+Pzwn8gRYfnyJyqQlACETiapom75BNjhGIJrBs/7I8nzFmZvcpztwhpDKRz2ZIjw2qFOsspvfVxxnc9ZKq5tRwD5u/8xnSY4OqukdevJ9Dz9+nqplNTbD1R3/H1LBuif+ul37O3oe+oarpZNMceuKHTA3pebQA5HNZ0qO631V5EJy8VdJWiVTmWZfIwKanmOrTXVg/0b2f17739+SzumteRw9uJTXcq6rpZKbI5zJYSpXHC5i8o1ZctUDB80Xb+9qyA15mq6gZCDK8dxMp5QDev/kJDvxSN4ADDG15jmT3QVVNY6SkrRKZkUExOzlGPpuh95mfkxnuU9VODXYTStRjh3TXvGYmRghW6ZahKlRZ1j5XY/Kq9f6A4+6A2kHRDsfIKVSbLkbEIlhVS2ZC1zo30tBGenSAfDajqju882VGdr1CdlLH4kCznqKI/FcR2S4i20TkhyISEZFvicgBEdnsbau8Y0VEbheRvSKyRURWF+l8UET2eNsHi/ZfIiJbvffcLiXUu5uRQTHZfZD9P/oqxnHI57JkRodUdEf3byU12E2koVWt+WzyzvH+tFBcPyhagZB6H105MsWCnlH2fQ6EojhpPXOpAqGqOjLjw2pB3BhDtHEWGMPEsX1MdOv430wc2UOouo6xvds49NM7VTQxbr9iKdvpEJF24L8Aa4wxKwAbuNV7+b8ZY1Z522Zv39uBTm+7DbjD06kHPgesA9YCnxORQnmoO4CPFL3vxun+vBkZFO1gmFRfFwADGx8nWKVTJaVv0+MM79lEarCHo8/8VEVzsu8QBx/7ARjD1OAxJnsPqugOH9xCaqQPOxRhsv8wjkJz3xjDRN8hTD6PWDZTQzohrkmhAAAgAElEQVTNx+zU+PGsNpeZIjmkUwi1oONkUowc3k56XKe/bmDnCwBM9Ozn4GPfV9HMToxw+NG7ATjwwF1kx4dVdE0ux+iuTeQmRsjn9JwSFe0IAkBURAJADDjdl38z8B3j8gJQKyJtwA3AI8aYIWPMMPAIcKP3WrUx5gXjZjHfAW6Z7oRmZFC0ipqLbVffgtg6WU0oUQcmT2Z8iIZll6toRuraGNnvPgj7tjxBuLZZRTeXHKPn1V+TTY5x6Ol7sAKhc9YUEY698ivGu/cwuHcDfa89q3CmYAVCbP/plwE4+NQPSI8NqOg6mSkOPXMPYDjw2LcJRqtVdPPZNBPH9pIa6kYsnVsolKgjUud99yav1pWSmL+UxILlAFjBc78GwJ28fQYDLY0i8nLRdttxHWOOAv8IHAa6gVFjzMPey1/0mshfEZFCp3g7cKToVLq8fafb33WK/adlhgZF9zNMLFxB1dxzrzRcIJhwM/KahRcSbdCpVRgIRwlWubp1nZcQCMdUdOMtJ2oeNnSuUbMOqG6/wG3yOznq5q5U0bSDYapa3fM1xpBo61TRrWqZTyjRAECssQMroFNotmHJ2uODV5F6vZqVzZdce7wbQbN/2U0MAsfvCw3OoPk8YIxZU7Qdb8N7TdybgfnALCAuIn8EfAZYAlwK1AOfUjvxEpiRQdEORRA7QNubb1bVDXkXasua61V1o96N1bT8Sj3NulasYAREqFu4evo3lEhN+xIA7FCUqlY964C6eRcBEG+aQyAcVdEUsWhY5JqBxZvnqWiC+7c3LHYtbqN1ekExVFVLw3K3BRKM6xXGDdU00LjmbWqZIqiNPl8LHDDG9BtjssBPgCuMMd1eEzkN/B/cfkKAo8Dsovd3ePtOt7/jFPtPy4wMilYoTOMlbyVU06CqG0zUUbPwIrUssUC0YRZVsxYRrW9V0xTLoqp5LtWzFhOK6d1gkbpWgrFqamYvUx3AqelYhlg21bMWq2kC1JchKAI0rrgSENXvDKB59TWEquvVB8eaLn0bEaXr1s0CVYLiYeAyEYl5o8LXADu9vkC8fbcA27zj1wMf8EahL8NtbncDDwHXi0idl31eDzzkvTYmIpd5Wh8A7p/upGbk0gErFKbp0rep64YS9bSsuU5dN1rfRqxp9vQHniHxlvmEE7oPBhGhun0JtbP1bFPB7UZIzFpMtZeJahFO1FPVtogq5aAYqWmibtGq410fWgTj1bSuu0lVE9x+26Z1eteuxooWY8yLInIf8AqQAzYBdwK/EpEmXEvDzcBHvbc8ANwE7AWSwIc8nSER+Vtgg3fcF4wxhSknHwO+BUSBX3nbaZmRQVHEQpQnLANE6lvV7U0BqmYtJBDWWzJXoKZjCZHaFnXdurkrSbQtUtdtXLyOeKO+Idasi28gGNMZZCmm/bJ3leV6qF20Sl0TdOerai3oMsZ8Dnc6TTGnzGi8EeSPv85rdwF3nWL/y8AZWU/OyKBYLspxA4BrPVkOYo362SdAdcfS8gSDObpe0gUKgzjalCPQQvmuMy0MQr5Cl/CVgh8Ufc6Yct2053sw8DlBhRVpPyPKGu5F5EYR2eUtsfn06xzzPhHZ4S31+UE5z8fHx0cBvYGW85KyZYoiYgNfA67DnTS5QUTWG2N2FB3TiTsn6U3GmGER0Zm57OPjU15mcKpYzkxxLbDXGLPfGJMB7sadqFnMR4CveUtzMMboVm/w8fEpC36meHacaunNupOOWQwgIs/iLgb/vDHmwZOFvKVBtwEEaupINZXhMVWmJ5/k9S+MsrnuVdrTv1xufmXQNVa5LrDf/ZdmgHwZruvzhTd6oCWAW7niLbizzZ8SkZXGmN+oyeQtDboTINI+u9JuXR+fmYWhfP6y5wHlbD6/3tKbYrqA9caYrDHmALAbN0j6+Picx2iUDjtfKWdQ3AB0ish8EQnh1klbf9IxP8PNEhGRRtzmtE4hOR8fn/JhStwqkLI1n40xORH5BO66RBu4yxizXUS+ALxsjFnPiTWLOwAHt7BkJZhU+Pj8HlO5gyilUNY+RWPMA7jrFYv3fbboZwP8ubf5VAj5TBpElKuuGPKpFFYopFb/0qeMVGgWWAozd61OmciNjuBM6Za3z2ez9K6/l8yA7oyksS0bOfjVL6lqAhy799v0P/gzVU2Ty7LvS39N8uBeVd3xHa9y4F++qKoJMLrhOcY26bovGmPIDg+RT+uaoqljwOSlpK0SmbFBMblvD8NPPqqu2/Oj7zP0+COqmiaTZuzlF3CSysE2nVYP4CdQvuALRlj5vK5uPk9ubFRXE0ge2Mvk3l26osZw6J//nuS+3aqymb5e+n/xM2VbXilxqzxmZFCc2LGN5J5djG/ZTGZwAOOcu7mQ8W7WfDqFFdZ1x8tnXe8MK6hTGbpAOQymXGEDyuuUC+uetf20xQ5gHEdd1wqGMMque2JZWJEI+amp49fbuZIdGSY3OsLo808zdWAfmT4lG90ZPNAyI4Ni6vBBRp5+nEx/L30//qFKH5VxHLq/fxfO2BiZ3m5GNzyvcKaQ7ushOzTg/Rs5siM6zoP5dAqcPGJbmHxe58FgzIngIqLnaFgcZIv/jXPVzeXAtt0g7jhqzdJ0bzcSCpHPZJjcu0vlfPOZNAMP/QKxbMa3blJr5TiTE3R/z62odeybd+g90P2gWFlEZs91f3AcqlevPf3BJWIFgziTkziTE0zu2EqgVqm4aD7P0W9/HYCuu76mliEMPf0ooxufx5lKcuSur6poAvTc912yw4Okjh6m/6FpixiXRG5kmJ6f/RCAkQ3PMvbKiyq64zteZfCJhwA4dOdXcJITKrpjr7zI2MYXmDqwl6EnH1Gp7mOFwmT6enEmJ5jau5tgQ6PCmUKkfTahJrempp2oJlCjUIW9MHm7lK0CmaFBcQ4AEgpTtVKvYGdknluXz65KEFugM8c81NyKHYt7P7cRqte5GaLzF5EZ6COfShFqaFLJlkUEO15FdmiA9LEjhOqbFM4UgnX1pLvdef3Jva8RnTNfRTfeufS4rsnlCNbpVCFPXLjazUKNITZfr9hu9cVrjv8cmTNPT/fSy1zNDr36mv7k7QojUF1DoKaWxMqLsMJ6Fbijc92gmLhotdq0EbEsot6NlVihF8CjcxYgXh9l1VId1z2A+OITNgTxC/QsCaqWuMWRg7X1BBt1iiXZ0RixBa7nS2yhnvdLeNZsgvVugI0qBsX4kuVYkSiBmjqCWi0RoOqii5FgiHDHHDVN8lLaVoHMyKAIEJ49h+o1l6lqRubMAxESqy5R1S1knVWKQdEKBonOXYiEQsQW6tm8RucudG+w1naCNYo3rhcUY526Vb0Tyy4EIK74GYgIiZWrkUCASMdcPd1AgKoLVxGZO09NE8CORKlaeRERxaAoprStEpmxQbH6knWEZys+GQE7GiVx0WpCrbNUdaPzFxGZPU81OwCIL1pCvHOZ6qi2FQwSW7CY+AXL1TQBwm0d2Ilq4p1LVXWrlqxAAgGi83Q9ZapWXkykY676jIHqiy8lOlen+6CYmksvJ9yu1HwudZClQoPiG10lp2zEF+u6whVouPGd6mXzg/WN1F2l7z4YX3QBgYS+j0jV4mWE23QNpsSyqF5xMTHl4GXHq6i/6hrsiO40qnBzK7VXXK2qCRDumEOgrl5dNzJHL6OFyh1EKYUZGxTLRaBK32RKRNQzL4BgYzOBWv0brGrZhVgRHcP6YuqvuhYrpLd0sEDdm/QfOEBZvjMRKcs1pk6FZoGl4AfF84RymDaJyPHBFk3saExdEzg+Cq+NdhO3wO+10ZbywqPzCT8o+vj4nBkzvMisHxR9fHzOmEodWS4FPyj6+PicOTM4KJY0JUdEHhWRm07ad2d5TsnHx+f3BRH5r57n+zYR+aGIRLxq/S96fvH3eJX7EZGw9/te7/V5RTqf8fbvEpEbivZP6z1/MqVmivOBT4nIpcaY/+HtW3O6N5QLEzRkmnL6wuXqIimHi1u5zrXSuokqydWwXKs73qCMTaP5LCLtwH8BlhljpkTkXlzbkpuArxhj7haRrwMfBu7w/j9sjFkkIrcCXwL+vYgs8963HJgF/FpECkuYTus9fypKnbw9AlwDtIjIz0VEYVW5j49PRWLQXOYXAKIiEgBiQDfwNuA+7/VvA7d4P9/s/Y73+jXiTgG4GbjbGJP2DPD24vrOl+I9/1uUGhTFGJMzxnwM+DHwDKCzQNXHx6fyKH1FS6OIvFy03XZcwpijwD8Ch3GD4SiwERgxxhSag124HvJQ5CXvvT4KNHBqj/n20+w/LaU2n79e9Id8S0S2Ah8v8b0+Pj4zjDNoPg8YY07Z1SYidbiZ23zc1uiPgBs1zu9cKCkoGmP+7aTfNwL/uSxn5OPjc/6j05d5LXDAGNMPICI/Ad4E1IpIwMsGi/3iC17yXV5zuwYY5PQe89N5z/8WM7YghA84o+PkBobVdVOv7SdzVKmsfRHJV7aTG9L1UzHZHNmjvZhsGQbnfp/RKQhxGLhMRGJe3+A1wA7gceC93jEfBArVjNd7v+O9/pjnCLoeuNUbnZ4PdAIvUZr3/G/hB8UzwBmfYOLpl8hnsqq6kxu20PMPd6hqAoze/2uGvvNTdd2R+x4k+dwr6roDd/yQ9J6Dqpq5gWG6//tX1B8OU6++Rvdnv6KqCZDef5ipLTvVdTUptWzYdE1sY8yLuAMmrwBbcePRncCngD8Xkb24fYbf9N7yTaDB2//nwKc9ne3AvbgB9UHg48YYx8s0C97zO4F7vWNPy4ycvG2MofcfvkrNLTcQXaZXXNQZGWPoOz8msnQRVqNeoQWTzZItQ+YlgYBbIVpb19LzZ/lN9A2xCt40Yus+/52JSfWsFiC54VUyR44RvVC3hNrg/7mX4JxZVF9zpY6g0hQjY8zngM+dtHs/7sjxycemgD98HZ0vAr/lZXsq7/npmHGZYj45RXLjFrJ9gyRf3sL4ky+o6Jp8HisWPf5vOGPjKroAViiEyeg6wwEQsDG5czes+i1E9K1IwYuJyvP5CoZdAV1XQ5PJYoV0C004YxPkp1JYkYiK0RiAyTkM37Oe9N6DZI/1ktw8baJUEn6R2QrCikUZf+gpTHKKyWc2YEV07Aiyx3oZvOseAPrv+B7p/YdVdFOv7SPbNwhOnvHHnifb06+im+3uRwI25Bwyh46qufk5E0mwLMgbnFGdB0M+mcJks7hREZzxSRXdbHcf+QnX99qks+SGx1R0J57eQObQUbAtJp7eQD6t80Ab/PZ9pLbtInO0h8G77lXRlIBNet9hcv2DTD79EvmplIruTC4yO+OCIkDVVScy79A8nWKooY428kn3gnIGhggv1CnaacWijP74QQBG7vsVAaVm+eTzrzD+yLNkjnQz+M0fqRlXDXz1u2SP9DD50hZG7ntQ4UwBk+fYX/0zOHmG7/4lU5tOu+CgZJyxCfr+5VsA9Hzhf6N2lxrD5NMbcIZGmXjmZaywTg3I8Nx2nNFxnIEhIksWqmgChBZ4Rm6RMLHVCn49Sn2K5yszMijG1q5CwiEkFiXQpOPgBhC/wvVmCbQ2YSeqVDSDs9sItLmueOHOeW52p0DskhWYdAYch8gyvWrWkaULyU8myY+OE1mh019rxWMEGmoBcIZGiV6kUzU9vHAuEnC7zYMdrQTqdBZiRZae+Dzjay9U0QQId3o2BIEAsUv0zMbC891ZKfF1F6sFcD9TrDCsSJjYpRcRntuBWHp/YnzdKrAstSwRvKrbay8CIHLBAjXd4JxZBJrcrFMreAFEL/Yc/ESIrtCxeQWOB8LQgtnYNTqVpyVgE1nqZlyxS1aoaAIEmuqxG+tARDV4hRbMAdsidtHS4/3XGoS9TDH+pkvVNCVf2laJzMigCFB15Vq1pnMBuzpBZMUFhDvnqerG1rkufmHFJpOIEF2zEgIB1fMNzm7DbqwjvGguVlyvAnf0QjcoHg+6ShQeCLHVutYBkaULCS+ej12r54FjhUOE5nYQu3y1miaA3VhPZOkiQnOnXeHmwwydkgMQmj8bK6prVgRQdflqgh1tqprBlkbCSxYQmqvrEhhbs4LMwaN6TSbcYBu7eBlWtU73QYHArGbspnpiq3SnokRXLGZiThuBZr1uFIDIkkXk56VVNQFiq1cQXa5nxwrud1Z36826I/sV2jQuhRkbFEWEYJt+zYroqmWgMGhxMnXvf5fKYEgxoXkdVF2t12QqEF29TLV5B56X8rVXEJil+50FmuqpvuktqprgZopY+uXAEm97k1q/cjGq90IFD6KUwowNiuWi0HGvTaijVV1TRIhdqjcQUCDcOc+dlqNM4prLy2IGVY7PQKvf82QkWCG3pB8UfSqVsrgEliFTLqvu77PrXrnwg6KPj4+Pi1C5I8ulUNbR51L9EUTkPSJiROQNsTjw8fE5A/zJ22eHiNi4/ghvB5YB7/e8FE4+LgH8KfBiuc7Fx8dHGX/y9llRqj/C3+Ia0CgtyvTx8Sk7MzgolrNP8VT+COuKDxCR1cBsY8wvReS/vZ6Q5+twG0CgqYZ4U1L9ZMNB3RqJBUIB/So1YbsMlW8AqbD2jjHlGUBJO/oDPplceQaR0lndSj2lUmGXyhnxhq1oEREL+GfgL6Y71hhzpzFmjTFmTaBabxWFj4/PWeJnimfF6XwTABLACuAJb8pEK7BeRN5ljHm5jOfl4+NzLpiZPfpczqB43B8BNxjeCvyHwovGmFGgsfC7iDwB/KUfEH18KoAKzQJLoWzN59fzRxCRL4jIu8r17/r4+JQff0rOWWKMecAYs9gYs9DzUMAY81ljzG85ahlj3lIJWWI5vEmGX9jL+PYudd29t/+a3oe3qeu++uf3cGz9ZnXdVz72PXof0SkwW8yOz/+M4Y0H1XX7HthMumdEXbc8/jfKzOA+xRlZOiw7OMbUnmntXc+Yw7f/kqNff1hdt3f9RgYe1fHOKGbyQD+T+3XsDYrJjk2Rm9CvEJPuGyef0TXaMsYw9MJ+cpO655vP5Dh4+0NMHR5U1XVSGbb9x/+PyZ26D8l8Osvk1gOYrMLMhVIDoh8Uzx+GH3iJ7n/9OYCqcZMdC5Pc16OmVyDcXE2mX8c/pJhAVaQswcsK2OQ1bq6TME5e3XUvn8pinDyBKt0ycpkB158m1KRXTxEgfXgAZ2yKYINewQmTc5ja3cXhz34HZ2LqnPUEv/lcUeRGJ4ku7iB9sJfuO35Bar9eEIvObyZ1qB9nMkV2aEJNN9RUTaZX3y4zkIiQm9CfEy9lcgk0OQcJ6F6ShQwxUKVjYFag8BALNekFr/TRIaYO9mPFwgQVg23fd3/N8C9fIjSrAeM4GOfch479oFhB5CdTHP3nHwMw8vBGrIjO5NaJbYeZ2HqY/FSG1z7+DXJjOhPIk/v7sONh0n1jdP/4JbWMMTs2RSARITuWYnx3DyZ/7leoyRucdBYJWpicQ3bs3LMOgFwyQz7nkM+5maIzpTORPnlkiHSfm9HZsRDZcZ0HRP/DWxl5cR9WLES6e0TtwdP30xfp+cFT2NEQx775qNqDJ5/JMf7ia2S6B+n5xq90snG/+Vw5hGY1UP/Oy47/bkV0qk7HFs9ifPMBALID4wQbdZ7kubEpDn/9UfKpLEe/87Ras+nw91+g//GdjL/Wze5/fBBRKIgqlrD9b37G1JFh+p/Yzf5/e1LhTN0MceP/821MzuHw916g91GdwZbsSJKt/69rFbr1L+4hO6K1EsrQ85MN5JMZ9v2vX2DHdbLQUHM12YFxsoPj2PGwWrFZO+qenwRsWj50g4qmHxQrjMb3XkWgwQ1aohQUrVCAlve9yf05GlK7ERIXzSE6z52uGZ3TqBK8AFquW0Z2OInJOlQv07M5qF4xi+xIknT/OLWr5qhoBqujBBIRMJA8PEjD5TpeNVWLW443FcMt1cRm69jHxuafqGLd8q5L1Oo1Bptct0E7HqHpXXoV062oew/U33wFoda6cxf0q+RUHlYkRMuHrnd/jur5k9RfdxHBpmqCjQm1G0FEaLnZrZgWnd+koglQtaiF+AJXr3q5nmFR41UnnAHrLpmrplt/6TwAqpfNItyg4/9ih4NULXYrmjdfp2dcFZ3TAJZgx8M0XqOnWxi0abplLXZcb2DIioUJNCRofM+VapoamaKIXCAim4u2MRH5MxH5vIgcLdp/U9F7PuOVItwlIjcU7T9lmUIRmS8iL3r77xGRaQPCjAyKAIkrlhG/eKGq34UVtGm99UpCSk3nAo1vW0agOkpsnl5QBGi5wbX11MwU4/MbibTXEl/YRKg+rqZbd6nredx4pZ5tKkD18llI0Kbpaj0zKCscJNJeR+P1K7EVH7qh5hovS9QtK2pFwzR/8Hq1riTQsTg1xuwyxqwyxqwCLgGSwE+9l79SeM0Y8wCAV3rwVmA5cCPwryJiT1Om8Eue1iJgGPjwdH/bjK28LSK0fvQd6qXo669ZSW5Ut0qPFQ7SdNOq481oLZreuoSuezcQmVWrpikiNF3ZqT7BOLG4hUAiQsObFk1/8BlQs6KDdM8YwWpdo63YghZa3qlrRRpsSND8nstUs0SA2Ip5x5vmWpShaXwNsM8Yc+g09+zNwN3GmDRwQET24pYoBK9MIYCI3A3cLCI7gbdxYnnxt4HPA3ec7kRmbFAECDXrBYMCErBp+oO10x94hrS842L16SihujjzPnSl+oOh8c2L1SdDi20x+99fSlQxgAMkls1CbP0SY21/uI7ILIX+uSLEtspybanfB2c2iNIoIsUr1e40xtx5iuNuBX5Y9PsnROQDwMvAXxhjhnHLEb5QdEyXtw9OXaawARjxlhyffPzrMqODYrmwyuC4FmosjzucZl9agarFLZicfpmU9vdcoq4Zqo1Rd+kCdd34ohZ1TSjPtVUWSg+KA8aY0/YHeP187wI+4+26A7f4tPH+/0/Afz6r8zwLKuQb8DlbtEazf0NTBAnqF021yqAJ5fkMfp8prGhR5O3AK8aYXoDC/wFE5BvAL7xfT1eO8FT7B4FaEQl42eLJ5QtPyYwdaPHx8SkfkjclbSXyfoqaziLSVvTaHwCFqibrgVtFJOyVJOwEXqKoTKGXdd4KrDdux/fjwHu9938QuH+6k/EzRR8fnzNDcWK2iMSB64A/Kdr9ZRFZ5f0rBwuveaUH7wV2ADng48YYx9MplCm0gbuMMYUKK58C7haRvwM2Ad+c7pz8oOjj43PGaDWfjTGTuAMixfv++DTHfxH44in2PwA8cIr9+zkxQl0SflD08fE5cyp0tUopVFxQDAUc5tUPqes2RfSq3hRTH5xU16wOlMcNNijlcQksF1lTnoGZsZzuPEGAoazeRPdi+lM6q39OZuc0r1fqEr5SqLig6OPjcx7gB0UfHx8fD9/Nz8fHx+cEZZineF7hz1N8g8lO5dj/1DFyad3+vFza4eef30zfHn2bg+9+fAO7nuyd/sAz5L6/2sym+/UNvJ7/zj62/lJft3/3CH27htV1KwJjStsqkBkZFPNlWII2sHOARz75kLp2cijF+j97hqH9usHLDlm8uv4I/fvHVXUBuneOMtqjP9hzdOsIE4P6njI7HzlGzy59u4dXvruLF+7QNxw78Mh+nvv7Z9V187m8WiEPv55iBTG4qYsnbv0WuWRGVTdcE6b7pWMMbNd1x6tuixOKB+jfrWuVKSIkmiOM9+kHr0DIUs9sAXKZPIGQ/iU5MZAm0aQ/oty/Z5TGTt3qMwCHHz9EakjH6uE3dO/fyot/+pNzD4y+m19lUd3ZRG4iw9GHXsMYo/ZkTMxKUDOvhqPP6TbDxBIaO2vVgyJAoinCxIB+5hWI2GTT+tl4OYKiMYaJgRRVDbrGVU42z/CBMRoW6QbFfC5P90vHaL+8Q02zcB8c/tlWqhc2qlRN0qineL4y44JisCpM+w1LOPTTLRy4ZxPZUb1Mqf2KDrqe68IYg5PW8ydu7KxhYI9+UKxqjDDer58pBsN2eTLFtIOtHBTT4zly6TxVypniyOFxnGyepk69slxOOsfA9n4y4xlmXa5XLf3w/dvoemAHk4eHmXPzShVNPyhWEGN7+jG5PJOHhtl1x7Ok+nUmZWcns8xa187w7iGe+7tnGdgxoKIL0NhZy8CeUUaOTDDRr9dsKjSfnWyevKPXlgmErbJkik42TyCsNyE7l3GYGHQfClWNepnisVcH6N81gh2yqJ2jN3n65ds3sP1726ieW0O4OqziwAgwvqefbV9+jFBtlL7nD5I/V5dAgz/QUklUzW8g2X1i0EIrKKaGp3jqvz8BwL5f7FErXLr/yWPsfbSL1GiG79/6sNqFtPG+g4x2TzGwf5zvfuQ5NOrMGmN4/Ot7yEw59Owa47E7dp+7KDB8NMkjt79GLp3nyOZhdjyq49W985FuHvyfboGVV9d30b9PZ9DpwNPdPPrFjdghmwf/5kWyUzqthuxEliNPHWa8a4zn/vYZtZJnqQH3HsiMTBFrr8ZSsOjwB1oqCCtgcfH/uJFYh9us0QqKiY5qVnzgwuO/i6Xz0c29ooXhQ+7Nmss4xBp0mnk1bTF2P9HDxEAasUTH4lSE0Z4pel4bY9/zA1hKN23trCib7u/CGHjhBwepbdOxDmhdUsP+592BsR0PH6Vxvk5Wl2iNkZ3KkZnI0rCohmBUZ7pv4UEbiARY+5eXTXN06aQG3KWmiz60ltY3K9k9+AMtlUUwEeGSf3gHgarQ8QtCg2XvX07DEregh1amaAdt1n3E9dipaopiaRiVAwsvb6Judgxwg4MWK284YYK1YJ2Op4yIsPAyV6uuI0bbUh1jsIZ5VQQjblZ0yXvnqWVeiTb3c43UhFh1q57RVsGOYvUn1hBriqnppgcmabl6IYs+oGN1UJi87WeKFUbVnDpWff7tpAf1gqIVsLj8r9+E2KLSHC2w9B3zqOmIU9WiZ64klrDmffMAaFuqFxTnXVpPrC5EOB6gfYWebiEorri+Tb7YCGYAABzZSURBVM1TxrKFliXV2EGLVbfMnv4NJVLd6gas1X+0mHBVUE3XsoXmi1pYfIue82A+6xBujHPhp6/Vq0BuSiswewZFZs8rZmxQBGi6dA5zbrlw+gPPgPrFDSz/o5VqzWcAO2ix7rblVDXrZQcAF908Bztk0aoYFO2AxfJrW5l3aT22otFWISguv75tmiPPjFnLall2fRuxOr2BlqqWGJFa3SwR3Gbz5X91hap9gpPKsfrv/h2BmJ69KTCjm88zfu1zzWJdL2WAiz58kVrneoElN87ByehOc4nVhlh182wa5+mWl1pxQxu9e3VXyiSaIix5awuzlul6arctraFujm7ZrnBVkCs+toJQXC9LBFh663LiLbrnGkyECSZ052hC5TaNS2HGB8VyYIcD2GHdj84KWKy4Rd917q2fXIqlbJ0695J6Gubq1we86VPL1O1YO69uJVqjG7wAVr57obqmdkAsGwao0KZxKfhB8TyiHK5zsVrlZhNuE7qmVddcHqCuXbf7AMrz94PvEFipTeNSKGufoojcKCK7RGSviHz6FK//uYjsEJEtIvKoiMwt5/n4+Pjo4I8+nwUiYgNfw/V0XQa8X0SWnXTYJmCNMeZC4D7gy+U6Hx8fHz380eezYy2w1xiz3xiTAe4Gbi4+wBjzuDEm6f36Aq5ZtY+Pz/nMDK+SU84+xXbgSNHvXcC60xz/YeBXp3pBRG4DbgN34uzSap1lYMUsjPSpawLMDg2qazZY+mZYADErWxbdcpHM6w+gAAzm9Qc8jmQapj/oLNgXai6L7ilvRA938naFRrwSOC8GWkTkj4A1wNWnet0YcydwJ0DrsvqZ+234+FQKFVoBpxTKGRSPAsXLCDq8fb+BiFwL/DVwtTFGv/ifj4+POjM5Uyxnn+IGoFNE5otICLgVWF98gIhcDPwb8C5jTHnarz4+PrrM8D7FsgVFY0wO+ATwEK639r3GmO0i8gUReZd32P8CqoAfichmEVn/OnJvONmMbk3CAhseG+PBHwyp6+7aluY/XNtFclK3nWOM4eoLj/L4w/rl8j/2x/187R/1vVR+9K1RPvtJ/WfuxKjDXX/fzciA7uomgExKv4ivHjprn0XkAu++L2xjIvJnIlIvIo+IyB7v/3Xe8SIit3tT/LaIyP9t79zjo6quPf5dk8lj8k54hBAi8gqEh0B4CFUBX7y0olUQ2iqordWitdrbK1dvba/V+7nVPmxrr9ZWrPa2WrUWubfiC+nVawUUgrxBkEcSQoCQ93Mys+4f50RHSpKB7CGZuL+fz/kwOXPOZs2ZM7+z915rr1UQ0tZi9/iPRGRxyP4JIrLFPecXEsbqgIjGKarqK6qap6pDVPVBd999qrrSfX2Jqmap6jh3u6L9FjsmEsLVUNvCvdPfY8ffzVdu27utkZefMpewtpUEn3Bgr58jpWZ/sCKC3w8Ndeavc1VlMCK9i327/dRWm58E27WpnhVPluONMx/I/ZtvbeMvD+813i4Y+o0YSDKrqrtaf/vABKAe+AuwDFitqsOA1e7f4IT3DXO3m4HHAEQkE/g+jiN3MvD9ViF1j/l6yHmzO/poPSohxO43inn2utUE/GZ/AL5kL/3zkli3wrzXe/g4HyUfN1NTaVa8svo708VlJeZ7MQk+obHRvMg0NSrxCeYFpqy0hb7Z5jJ6t7JrUz25Q+NJTjXbdtWRJna8W8HgAvNFsfb9Xyl/+PIbnVu7rxEpR3AxsFdVD+CE7j3t7n8auNJ9PQ94Rh3WAukikg3MAt5Q1eOqWgG8Acx230tV1bXqFGt6JqStNulRothvdCbl+6rZ+AczGaFDOffKfmxeXU5DjVmRyRuXiAjs3mR2OJrg85Dey8PhSIlig/kuXXNTZETxSGkLffub9ynu3tRA3ljzyx3f/+sRElO9jLog02i7LU0B3nqo0Exi3PB7ir1F5IOQ7eY2WlwIPOu+zlLVUvf1YSDLfX2yML+cDvYXn2R/u/QoUUzNTmTqN0by3q+3U11a3/EJp8D4mb0RD2x81WyJ06SUGHKHxrOz0Ky94PQWyw5FQBQThIYIiGJjpHqKh1o+6TmbQlXZ/WEDw8eZX6+9fmUZEy/ra7yy4fu/20lDRRPT7xzb+cbCd7QcU9WJIdsTJzblOmKvAF74h//G6eGdUZdNjxJFgIKv5JHWP4k1Dxc65S0NFYJKSPIy7tI+rHs5AkPo8YkRFEXzE/aR6ik2NSoJhkWxtjpIXY2SlW1WFA/tb6a2KkDeOLM9xeKdtZTsrOPceVkdHxwmtUcaqCyqZf3ynZz3zdEk9e58yQsJBsPawmQOsFFVy9y/y9yhL+6/rV6ytsL82ts/4CT726XHiWJMrIeL7ylg75pDrPrX9Wxdsc9Y2+demcXeDdUcPdjAsSJzw93h433s/rCBgGEnUb8cb1QNn5salbh4s6LY6mgy3VPcVVhPvE8YmGeudOqxogbWrywja3AiZ41OMdKmqrLijv/jjQc2kDkolbHzDaQ8U5zg7XC28FjEp0NncEL3Wj3Ii4GXQ/Zf73qhpwBV7jD7NWCmiGS4DpaZwGvue9UiMsX1Ol8f0lab9DhR1KBy9KMqEnvFs/OVgxxcZyYUI9Ci5OQlkd4vnifv3M6by4s7PikMgkFlxLhEGmqD/G1FJRv+10zy1mBQPxk+lxb7qanufI9RVfH7FZ9PaGpUqqvMOFvqaoO0tChNTRCfIPj9ZgT3413NlBY7otg3O4agoQQFK548xpZ1dQwd43PWvBkgGFQeWfwh61eWUTC7D7UVZpZcVuyv4cjOSorWHyGpVwKHt3Y+/EtQRMPbOmxLJAm4FHgpZPd/AJeKyEfAJe7fAK8AHwN7gN8A3wRQ1ePAD3Fio98H7nf34R7zW/ecvbS/ghHogaIoHmFAQe9PkpWWbi43kiXbEwN/emAPlYebKNpea8zhcqzUz0++4wjsz/+5hOZGMz/cl/9Yw2t/qaXsUAs3XHaIhITOf9UiwrKl5Xy8p4V31jTy4D1mQpRqa4IsmOWMnJ58tJpXV5qZSti0vpH77zxKjBfuWlxGdaUZES98p5bVL1ZSvKeJR+/pcDQWFk31ASpKm6gp9/PqYwc4tNvM+vaD73/aKfBlxJE1ypDzxlDdZ1WtU9VeqloVsq9cVS9W1WFu2N5xd7+q6lI3xG+Mqn4Qcs5yVR3qbk+F7P9AVUe759zmzlG2S48TRYA+eeksevoiMgelEPAHKdnU+ThAEeErP8yjz1nOHFJ9tZkned+cOMaf/2m5gJQMM+EdF1+exIG9foJB6N0vhlhDsXRjJ8ZzuCTA4ZIAY8abSeCale2ldfpp1zY/0y42M083ZEQctTVBAi2QM9BLeqaZa5vqfke11QGuucVMuYvGkIfsnKUDGT4lo52jw6dovSOKE67PY/b9k4mJNfSTNySK3ZEeKYoAqf2TuHb5hWSf08vYENqX7OWmR0bijRMaqs3N1V1za59PxDA13dAPNz2G2V9yxHbIcHPZp2de/qlgTT7PXO2PyV9w2ppyQQJp6YbKvIZ87msWm6v9kprhzE/OXpRJziAz16DeFcW8KenMucVMrmUNKsUbjjL9rrFMv3OswWp+mJ5T7Fb0WFEE8KXHc83j04wVKwfIzU9m/r1Dqa8x59VNTo1h0becFFCmeooA85c4QjDYoCj26+9l3MQ40jM9DB1uLnXXuec7DouZXzTnzU1O9ZDVP4axk+LJG2VOwFMyYkhM9rDwdnNpuxpqAqT0jmXJQ/l4DNUUL99XzYzvjmPCdXlG2gvFsPe5W9EtUodFkliflylfzzfa5nnzszmy32yw9exFmfz19+Ukp5kTxUHD4ph0foLRniLArC8msnF9Ex6DdUomTo0nIUGYcanZEJehI+KYc7UZT24rqZkxXHNrH9J6mfv5NNW1sOShfNL6mPuueg1KpfcQ86tiIHqHxuHQ40URzBcZEhG+eMcgo216Y4Vb7+9PrOGA3QU3pHH2MLPJWC+Z6yPOcNXMtHQPNy5NISXV7OefelEiF8w0G2A9eKSPQfnmQnEA8s7NIDbe7GePWHEtxYqi5R8xfQMDnDPVbH1mgCkzzC9D6903hsuuMp+d+vqbzfboAK5YmEKMoeFoKyPGm1/FEon7KaJE58g4LKwo9nBMDnFDicRyvEi0aVoQLQ49OcmsFUWLxXLqWFG0WCwWF1UI9Nzxc9SJYpynhbMTzCdlHRFf2vFBp8Hg2Grjbfb2mPUmt5LoMes8iDT1weaItHssaP7+SpLI2Broqqg621O0WCyWEHqwKEaZy6vrCGPJ5Gm1WXSwhfJy80OR4pIWXn/LfB2VlhblD3+u4dBh89l3Vr9Tz/rCRuPttl6LSHyHH+1uoc5wHZxWImGvERQIanhbFNJjRDGSN9CDd5Ty+1+aLy4FcPVVx1nxknnxevvvjSy+5Zjx66IK199WxofbzFejfejRCp553kyWoFBWv93I124vJ4yaRafMkuuP8/tnzOfC9DcrSy7ZzwfvmEkMcTJO/95Q0GB4WxTSI0Rx85tHefL2rREpWgWQfVYsq56vMp7vUESYMCGWjRvNJJcI5awBXmrrlMoqszbHxgpeL9RHKPO2LwJhOSWHWhjQ33yNlrLDAQ6VBCkoMD/Hu/atWg4d9DN4hOEoeZcDW6r55ZJNp5fYRHEcLeFsUUiPEMVeA3zsePc4r/zSXELZUGbPT+XIoRYK/26+RzC+II7CDeYn4QfkONPFRcXmh7mJPg/1DeZv+IYIiWJRSYCcCIjipkI/Xi+MOcfsiiGAVc9XM+XCJDL7mJ/2r61oZvkdW4n3xZCQfJrt2yw53ZucEcks/LfhvP7rA2xebbaGCkD/s+IYP9XHK38yX5O4YEIspaVBDhkuG5CTHYMIFEWgRkuiTyLSU4yUKBaXtJCbY15cNm70k5/vxecza/PRw37ef7uOOQvMr1sOBpSn/2k7Hq9w3Y/yTz+434pi92fSFf2Y9pUc/mvZDo7sM9+jm3ttGu++UUtluVmRGTMmlthY2Gi4txgbK2RnxVBUbL5GS6IvMoWrGhqD+Awkwz2RopLAJz1nkxRubGZcBIbOr/+5mrTMGCZPN7+UctWj+/h4QxVf+8UYEtNOt4cbpiBaUex6rvznoWQPS+LJb22lqc6seJ0/Kxlfooc3V5h1BMQnCKNGx1K4wfy8Yu6AGA5G2fDZdOEqVaWkNEBujtnhs9+vbNnsp2CC2aFzMKiseqGamV9KxRtr9lpsWXOM1x4/wLU/GE7OiE6ss1cgGAxvi0J6lCh64zzc+Mho6ir9PPu9XUY9r3HxHi65KpVXnq8y7tEtKIiMsyU3x0txifmeoi8hgsNnw0PRY8eDNDaq8TnFnTtaaGyE8ePNiuKH6xooPeg3PnQ+eqCe/7p7B+cvymHylf0636DtKUYPaX3jueFno9j0+lH+9oyZ4lKtzJmfxoGPmtluOJauYGIc27b6aTJUn6WV3AFeiiJQzc+ZU4xETzFofE6x9aFgek5x48ZmMjOFswaaFdtVz1cxeqKP3MHmhuXNDQGevGMrWYMT+dKyoQZaVOt9jjaGTExn3neH8PLDe9nzfiUbV5kpRzAkP57h58R/4nAxVR2uYEIsfj9s2WK2t5jbPyYioujzeWgwLOCBgNLcjPE5xaKSFjweyM4yK16FG/2ML4gzFvsYDCq11QHeebWWOfPNlE44eqCeg9tqeO77u6g+2swNj4zCayJfp4JqMKwtGumxy/xmXD+AA5ureerOrTTWBjhrdAq9czufW3DutWk89uBRBg6L45xJPkaM7Xyb2dkxZGd7WLe2mcOHA1xuICV/IKDkDvBSdiRIUUkLwSAMzO38111TG/zE+1xU4ic3p/PDx6PHArRqiy9BqK4JkprS+R/vxg+bKC5poX+/GGJinGvS2VRix48H2fKhn8INfhYsNJer8pmflxMIQIwXps81k1fyvT+XsuF/yqgsa2Lp8nFk9DO4tj1KV6uEQ4/sKQJUHm6iuSFATbkff1OQtX/ufMKHkv3NFL7XQGO98ut/P0Zleefn64JB5bt3VREIwE9/XMvqN82sFFn9v4088HAlAJNmlNLcbOYmvvHbZbz3QSOrVtex9F/MhD/VNQQZe9FBAJYuO8qad81ED7ywop77f1RFbV2QyxccMTLFlZAANy6p4ODBAG+83sR/v2xmNdKe7U388T+Pk5waw7OPH+/0KCTQEmT9Xw5TUdpEQrKXfRurjI1sADunGI1kZCdw4ZJcMrKdFQFrXyol0NK57nzO2XH06vvpMKyqovOi6PEIcy9P4MiRIMEgpBuqZHfJjARqap2bMhiEwWebGRRcdmkSpWUByiuCTJtipqd0dm4sia6D5cixALNmmMlsPTo/lrp6Z1XPJRf68Ho7P9T1+QSv2zmuOB5k5mwzva/qSudeqjoeYNqc5E4nB97+9nGqjzlhXvnnZXLBl3PMJRxWtd7naGXY5AzuXjGJCZf1pfpoM9vf7vz65Zvv7sOYSY4YVB0349mdcWEc51/gTKybEkWPR7jlRmcYNiIv1lgG6itnJxHn+gAuPM/c8PGi8x0hvGJWEgmG5hVH5TuG+hKExYvMxPyJCBnud3T3PSnEx5u5rq330rcf6MvQkZ0X2vdePER8Ugxf/Y98Fv9kZCdiEtvA9hSjl8TUWBb/eBTXPzySwlc773Dxxgr3PZpNr74xxkRRRLj3eyl4PJCeYc77+pUFSaQkC/kGS5Gmp8Uw68Ik0tM8jBttbl3uRRc4Artgnrk6NSOGxeLxwMKrk8g0WDo2PcPDxEmxzJlr7vNXVQS4bFEas67ufChOZVkTDTUtLFsxicnz+kUgEYaigUBYWzTS40WxlYmXZ3HFd4YYmVfJ7OPle7/sT221ueHBiPxYFiz0kZ5m7itJTfFw3bXJjBphtpdw7bxkpk/1Ga1/ctF5iWRmeLh0mrmiUAkJwrAhXm650WxBsMxMD/d8L8WY2AQCSs7AOG67r4+R9kTgtqfG0WuA+aJlgNHUYSKSLiIvishOEdkhIlNF5AciUiIim9xtbsjx/yIie0Rkl4jMCtk/2923R0SWhewfJCLr3P1/EpEOY516rPf5ZKRnmXuyj5nkIzXD7DPlru8kc+CA2afrN25M5mCR2Ta/ODPJePB2n94x3HdXJnFxZns13/xaKiNHmF2Kd93iRMaPN9emv1m579Fs4gxV9EvrG5nMOp/BXLjNz4FXVfUaV7ASgVnAz1T1x6EHishIYCEwCugPvCkiee7bvwIuBYqB90VkpapuB37ktvWciDwO3AQ81p5BEe0ptqXeIe/Hu+q9x1XzsyNpj2kGDjV78/XpG8OEiWZ7dYPPjmXaeWbtTE7ysHiB+XKkty4xnwBhyZfNrx+ePcfs9UzwecgyENp0plBAgxrW1h4ikgZMA54EUNVmVa1s55R5wHOq2qSq+4A9wGR326OqH6tqM/AcME+crvxFwIvu+U8DV3b0+SImiiISg6Pec4CRwCJX6UO5CahQ1aHAz3BU/XNNJBKhRqLMqQlP7ploMxKfPRLfUVShp5RktreIfBCy3RzS0iDgKPCUiBSKyG9FpPUpdpuIbBaR5SKS4e7LAYpCzi9297W1vxdQqaotJ+xvl0j2FE+q3iccMw9HvcFR84vlc3/HWSzdn1NwtBxT1Ykh2xMhzXiBAuAxVR0P1AHLcIa3Q4BxQCnwkzP52SI5p3gy9T63rWNUtUVEqnDU/TPl1NynS+sTpulb+Wu2RsTiyNCbEz5PNyaabIWosvdwhGzdbb5Jh+FtvVFDxWtv6ou9w2ynvc9cDBSr6jr37xeBZapa1nqAiPwG+B/3zxIgN+T8Ae4+2thfDqSLiNftLYYe3yZR4Whxny5PAIjIB6o6sYtNCptosjeabIXosjeabAXH3rbeU9XZJv4PVT0sIkUiMlxVdwEXA9tFJFtVW5egXQW0doJWAn8UkZ/iOFqGAesBAYaJyCAc0VsIfFlVVUTWANfgjFQXAy93ZFckRbE9VT/xmGIR8QJpOOpusVg+H9wO/MH1PH8M3AD8QkTG4fh09gPfAFDVbSLyPLAdaAGWqmoAQERuA14DYoDlqrrNbf9u4DkReQAoxHXqtEckRfF9TqLeJxyzEke938NR87e029Z1tFgsplHVTcCJPejr2jn+QeDBk+x/BXjlJPs/xvFvhE3ERNGdI/wH9RaR+4EPVHUljmr/XkT2AMdxhLMjnuj4kG5FNNkbTbZCdNkbTbZC9NlrDLEdM4vFYvmUz80yP4vFYgkHK4oWi8USQrcVxWhaIhiGrXeJyHY3Qn+1iAzsCjtD7GnX3pDjrhYRFZEuCyUJx1YRWeBe320i8sczbeMJtnR0L5wlImvcFRybQ5MdnGnc1SJHROSkcb/i8Av3s2wWkYIzbWOXoKrdbsNxzOwFBgNxwIfAyBOO+SbwuPt6IfCnbmzrhUCi+/rWrrI1XHvd41KAt4G1wMTuaitOrFohkOH+3bc7X1scB8at7uuRwP4utHcazoqSrW28PxdYhRMHOAVY11W2nsmtu/YUo2mJYIe2quoaVW3Nsb8WJ2azqwjn2gL8EGctutnShadGOLZ+HfiVqlYAqKqZKmWnRzj2KtBamSoNOHQG7fusIapv40R9tMU84Bl1WIuzOiT7zFjXdXRXUWxrgfdJj1FnCU/rEsEzTTi2hnITztO3q+jQXneYlKuqfz2Thp2EcK5tHpAnIu+KyFoRMbLa4jQJx94fAF8VkWKcuLrbz4xpp8Wp3ts9gqhY5tdTEJGv4gSqTu9qW9pCRDzAT4ElXWxKuHhxhtAzcHrgb4vIGG0/BVVXsgj4nar+RESm4sTpjtZorQfaA+muPcVTWSJIFy8RDMdWROQS4F7gClU1U7Lv9OjI3hRgNPA3EdmPM5e0soucLeFc22Jgpar61cmxtxtHJLuCcOy9CXgeQFXfAxJwElt0R8K6t3scXT2p2cYErxdnHeQgPp2wHnXCMUv5rKPl+W5s63icCfhh0XBtTzj+b3SdoyWcazsbeNp93RtnuNerG9u7Cljivs7HmVOULrwfzqZtR8tlfNbRsr6r7Dyj16SrDWjny5qL89TfC9zr7rsfp6cFzhP2BZzsu+uBwd3Y1jeBMmCTu63sztf2hGO7TBTDvLaCM9zfDmwBFnbna4vjcX7XFcxNwMwutPVZnHyFfpwe903ALcAtIdf2V+5n2dKV98GZ3OwyP4vFYgmhu84pWiwWS5dgRdFisVhCsKJosVgsIVhRtFgslhCsKFosUURHSRxOo72H3EQaO9zkD5/7appWFC2W6OJ3OLGZnUZEvgCcB5yDE7A/iW682upMYUXRYoki9CRJHERkiIi8KiIbROQdERkRbnM48b5xQDwQixNP+7nGiqIlIojI/SLy7ZC/HxSRO7rSph7ME8DtqjoB+CfgP8M5SZ1lhmtwArhLgddUdUfErIwSbEIIS6RYDrwEPOImmVjIKVZVs3SMiCQDXwBeCJkOjHff+xLOapoTKVHVWSIyFGepYWsquzdE5AJVfSfCZndrrChaIoKq7heRchEZD2QBhapqa3qbxwNUquq4E99Q1ZdwHkxtcRWwVlVrAURkFTAV+FyLoh0+WyLJb3FSkN2A03O0GEZVq4F9IjIfPikhMDbM0w8C00XEKyKxOE6Wz/3w2a59tkQMEYnDSSQQi5MhKNDFJkU9IvIsTu7I3jhOke8DbwGPAdk41/o5VT3ZsPnEtmJw5h+n4ThdXlXVuyJjefRgRdESUUTkcZzhXZsFsiyW7oSdU7REDNfBMgWY39W2WCzhYucULRFBREbi5LpcraofdbU9Fku42OGzxWKxhGB7ihaLxRKCFUWLxWIJwYqixWKxhGBF0WKxWEKwomixWCwh/D/nkxwGomr9lAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "system.hamiltonian.effective_field.plot_plane(\"x\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Whereas, the individual exchange effective field is:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2018/3/26 16:32: Calling OOMMF (system/system.mif) ... [0.7s]\n" ] } ], "source": [ "Hex_eff = system.hamiltonian.exchange.effective_field" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Because we initialised the system with the uniform state, we expect this effective field to be zero." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(0.0, 0.0, 0.0)" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Hex_eff.average" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Relax the system" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2018/3/26 16:32: Calling OOMMF (system/system.mif) ... [1.0s]\n" ] } ], "source": [ "md = oc.MinDriver()\n", "md.drive(system)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Compute the energy (and demonstrate that the energy decreased) and plot its magnetisation:" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2018/3/26 16:32: Calling OOMMF (system/system.mif) ... [0.7s]\n", "The system's energy is 5.352855331445172e-20 J.\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUsAAAERCAYAAADyooUPAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsnXd4m+W5/z+3JMt7O85w9h5kkAkJO+xSVpiFAg0FCrSlhZaWtvRQzqHnAD96Ci1toS2lk1koBMIIEGZC9rQznekR7yVrS/fvD8mJSeNYll4dsPN8rktX5Efv+81jW751P+v+iqpiMBgMhqNj+7w7YDAYDL0BEywNBoMhBkywNBgMhhgwwdJgMBhiwARLg8FgiAETLA0GgyEGemWwFJGnRKRWRDZbpPeQiJSKyBYReUxExApdg8HQd+iVwRJ4GjjXCiERmQvMA6YAxwGzgFOt0DYYDH2HXhksVfVDoLFzm4iMEpE3RWSNiHwkIuNjlQPSACeQCqQANZZ22GAw9Hp6ZbDsgieBb6nqDOB7wG9iuUlVlwNLgero4y1V3ZK0XhoMhl6J4/PugBWISBYwF3ih03RjavS1S4H7j3BbpaqeIyKjgQnA4Gj7EhE5WVU/SnK3DQZDL6JPBEsiGXKzqk47/AVVfQl46Sj3XgJ8qqouABF5AzgRMMHSYDAcpE8Mw1W1FdgtIpcDSISpMd6+DzhVRBwikkJkcccMww0Gw2dIWrDsbntPNKA9JiI7RWSjiEzvgfYzwHJgnIhUiMiNwDXAjSKyASgFLopR7kWgHNgEbAA2qOqiWPtiMBiODSRZJdpE5BTABfxFVY87wuvnA98CzgfmAI+q6pykdMZgMBgSJGmZ5ZG29xzGRUQCqarqp0CeiAxMVn8MBoMhET7PBZ4SYH+nryuibdWHXygiNwM3A0iqc0bKwH7W90aSk2E77GHLNZ32oOWaAOn2QHJ0bcnRTZXk6DqTcH7LTnIOhUmSdNds9NWr6hH/0M45PVMbGkOx6rylqpYcIPm86RWr4ar6JJF9lKSOGKwD7vum5f+HIzW2X35PyctxW645Iq/Bck2ASdn/9jllCcelVyRFd4yzNim6g+zWvxdybWmWawKkiD0puvaBO/d29VpDY4iVbw2NUWdHkWWd+pz5PINlJTCk09eDo20Gg+ELjAJhrB8xfdH5PLcOvQpcF10VPwFoUdXkpDYGg8EyFCWgoZgefYmkZZbR7T2nAUUiUgH8B5Fz16jq74DFRFbCdwJu4GvJ6ovBYLCWYzGzTFqwVNWru3ldgduT9f8bDIbkoCihY9AVtlcs8BgMhi8WYUywNBgMhqOiQMgES4PBYOgek1n2ITybt+NZuwVJsSMpKYgj8q+juID0mccRj3OEBkM0v/ohGghG9Bx2cDiQFDsZU8eS0r8grr62rtuNd18dYreB3YbYbIjdRkpRNjnTRsSl2V7VSuueZkDRsIKChhWxCwNOGILN0fP9eeFgGJ87iN8dwu8OEvBE/vV7QgybUUhqZmJvJ1XF3RaiqSZAY02AoePSyeuXkpBmB+2uMNu3BCjsZ2focOve9uW7gjQ2hZk1w2mZpscTZsn7Hi44JwObzbpN51t3+HG1h5k5LbE9nwoEzJxl3yDkchN2uWn/eA3qj57yECH7rLlknXFCXIESINjUSqC2EdcH6w622QtyKLrhAhzF+XFpqiqoUvHH9yB0aIWx8OyplCw8Iy5NgJSsVNY/8hHuA66DbbmjC5j549PjCpQAXleQv9+ynOotLQfb0nJSuPD+aXEHSp8nzBM/2Uv5RjeNNQF8njBig2vvLmHyvOy4NFWVlZ/4KNvkZ+vmAFtLA+zfE+S8izL42SPx/Z46CIeVNesCvPmWl7fe9uLzw5uvWbPvuromyJN/buP3f2nlf/6jwLJAWVkd5P5HGnlzqZuyj2LbTH40FDXD8N6KBoL4du7FW7oTb+lO/HsqwSbYsjJRfwDniMEU3HAJzmGDeqQbdnvxlO3Gs2knno07CRxoQNJTwW6DsJJzzgkUXHEmtvTUHukG2zy0rd9D65pyWtfsItDoimgCqYMKGPqt88ieMqxHmgBBd4CaFfup+ngvB5bvw9/iBUDswvjrjmf89dOxpcSRUYaUPavq2fxGJY0V7Qfbh84o5JL/nk7ugPQea3ZQV+kjGFCq9/gASM+y8Z1fjuD4U3Pj1hQR3O3Krx5spSMBuvDyDH76YD52e/wBaPfuINdc30j5rsj+wcxM4fVXCiksSGy7ck1diB/9VyPP/8tFIABnnZbONZdlJaQJ0NwS4uHHm3n09814vMpvHuxHZoYFW6sVQsderOz9wdJXvo/aB/+A+gOkDComddJoci86g9RxI2ld/AH2vGyyTp+D2Hr2JnFv3MGBB/8KqqSOHkzm3ClkTBmNc9gAqn/+NEU3XEDqyJIe99ezp5Yt3/ojYhMyJw6h+OLZ5MwYSf2b67GlOxl41TxsqT0femoozJuX/wN/m4+iaQOZ8LXp2FMd7HxxMzN/fBr54+I/T7/4gY2sfXEvg6flc/rt41nySBkn3zyGk74+FlsCwadsZRv3XbODkpGp9B/ixGYX7v7dKEpGJTZMDAaVh+9vxm6HYBAWfCWTHz2Ql3CmNmKEg1kznZTv8iACv/lVHhPGJz5N0L+fnf5FdgIByMwQfv1QYdyjn87YbEL53gAerzJmZAoLr85JWBM6TvAceyStRFuyOPxseNjrw716M2mTRuPI/2w2osFQZF4xBg4/Gx5yufGW7SZt4kjsWYcyJw0GQQSxx6Z7+NlwDYVpWbWT7CnDsGccykiDbR4c2bFlaF2dDa9fX03OyAKcORHd9gNtpBdmxJxNdnU2vGGPC3uKjbySDHztQWq2tTB0emFMmtD12fBQUNmy2sWkOVks+mMtZ1xeSFZu7J/fRzsbvmaFj5amMCuXefnBz/J6FHyOdjZ8y9YAP3+wjZkznNzxzdizv+7OhlcdCPKzh5qYOsnJbTfGnlV3dza8fE+A//xFIxecncllF8TeX/vAnWtUdeaRXps8xakvLY5t6mHskOoudURkHPBcp6aRwE+JeGL9DsgC9gDXqGqriDiBJ4CZROL1Har6flTrSuDHgB14TVV/EG0fCvwZyIu+9kNVXRx97R7gRiAEfFtV3zra99Lrg6VVmEIafa+Qhrs9THqG9DhL666QxpatAcaPc/RIN5ZCGnUNIQrybD2aKoilkEZzS4jcHFuP+nu0YHncFKc+/3psI5VJQ6u61OmMiNiJ1IaYQ6Qg9/dU9QMRWQiMUNV7ReR2YKaqfk1EioE3iFhX5wPrgBmqWicifyZS/vFdEXkSWKeqvxWRicBiVR0eff4MMBsYBLwDjFXt+oxmn7CVMBiOREZmzwJErEwYn5IU3X6F9oTmVLsiL9duaX8j+ywlpkcPmA+Uq+peYCzwYbR9CbAg+nwi8B6AqtYCzUSyzJHADlWti173Tqd7FOiYf8gFqqLPLwKeVVWfqu4mcux69tE6aIKlwWDoMWGVmB5EakOs7vS4uQvJq4hkevBZW5jLOVSdbANwYdQvawQwI/raTiIWM8NFxAFc3Ome+4Bro/UpFhNxZ4Cu6+l2Sa9f4DEYDP+3dGSWMVLf3TA8Ohd5IXBPtGkh8JiI3EukOpk/2v4UEdvq1cBeYBkQUtUmEbmVyPxnONo+KnrP1cDTqvqIiJwI/FVE/s3mJhZMsDQYDD1CEULWDkrPA9aqag2Aqm4FzgYQkbHAl6LtQeC7HTeJyDJge/S1RcCiaPvNRBZtILKAc270muUikgYUEUc9XTMMNxgMPaYHw/BYuJpDQ3CiizeIiA34CZGVcUQkQ0Qyo8/PAoKqWnbYPfnAbcAfonL7iMyHIiITgDSgjkjGepWIpEaH9GOAlUfrpMksDQZDj1AEv1pjZxENfmcBt3Rqvjq68g3wEvCn6PNi4C0RCRPJAr/a6Z5HRWRq9Pn9qro9+vwu4Pci8l0iMwg3RMtDlorI80AZEARuP9pKOJhgaTAYekhkU7o1g1JVbQcKD2t7FHj0CNfuAcZ1oXPE+rnRzHNeF689ADwQa197X7AMAz7rTZpCtuTsN/UFre+rN2RNcYnDcYetKwbRGa8mp79ei7Kbf9e13j0zo49ZLPRwW1CfoPcFS4PB8LmiKoT02FvuMMHSYDD0mLDJLA0Gg+HoRBZ4jr3Qcex9xwaDISGsXODpTZhgaTAYekwo9j2UfQYTLA0GQ49IwgmeXoEJlgaDoceEzWq4wWDoKeGwWmosBhAIKA4HlpeC83jCpKcnFugihTRMsOwz+Hbuof4Pf0cDQQiHo86GNnK+dCbZp8/tsc0EgIbDNP75FbxluyKaoTCEwjj6F1D4tUtIKSmOq69tK7fT8NpKNBAk7A9G+hxSihbMJff0KXH9wbj2N7Hv1VKC7X6C7X4CLh/Bdj+F0wcz5obZ2J0939DtawtQtbEBV50HV50XV60HV52HjPw0TrtzCqnZ8W8+b2sKULnTQ8V2DxU73HhdIa66eyh5xYltlFdV9uwMsOpjL5vX+lh4Rx7DR1uzSX7z5gDPPuPhjPmpnHFGz3yYuiIcVl5e5GbdJj//9ZPEzNU60+YKc+/PG7n+qmyOn5JYXxUhkKQDAV9k+lyw1FAI347duNeXEvb6UE/EtCtlyCAKv3YlzpKBcemGvT48G7cTam4jWHOoUnn22XPJu+IcbM74/gCDLe0EGtrwbKsg7I04UToHFTD4zktIH9Mzg7XO2J0O6lfuo213IwDisDHh1nkMvyy+4AuQkuFg48u72fHOoeIsY88azOl3T8WZHt9bKegP88Td5axY3Hiwrf+wVH7w9ISEAuX7b7r54G03qz/2Ul8bwm6Hh/5QnHCgbGsL8+orXp75h5tNm4JceGGaZYHy/Y+8/PTnzWzY7Gfl0vjep0firffc3H53PUMHOxIOlACqmE3pvZWw34+3dDue9aV4Nm0h3O7GObSEtPGj8WwoI+f8M8g9f37MvjkdhNrceNZvwb26DM/mHRAOkzZhFPaCXBCh6KYFpE0c1b3QYQRb2mn9dButn5TRvmkPNmcKzoEFeHfXkH/OdAYsPAtbWs8DRaDdz4EPdlLx1jYa11XizI1YGaQPzGH6feeQN6F/jzUh4hde/mE1658vZ++nh6wc5t0+iRO+Pj6hoWLFTg9+3yH7q6HjM/j+U+PJLUosqA0c7GDpYjc+b+QY690/L2Du6fG7UELE3fGG65rYsydydHHgQBv/9UDiJmB79we54weNvPtB5IP9hmsyGWdB9lvfEOL79zXwzD8jjpwP/yx236SjI2ZTem/EU7ad+t/8GQ2FSB0zktwLziJ92kQcBfl4SreRe94ZOIcN7rFu+8pN1P/mOcRhJ33KWAoXXkL6tPHYM9NpfuU9cs6eiy295y6E7q0V7L7naWzOFLJnj2XIDy4j6/hRtH5Shi0zjZw5R6wT0C0hf4ilV/yZkD/EgJNHMOvhL5M1JI8tv13GlLvPICU7/ozi7f9cS+miPYw6ZRCX/eYklj2xhdk3jGP06fFnvgBbV7by82u3MGpaFtPn59PaGOCuJ8eR2QPTsiMRDCo/urWO7FwbPm+Ir383lwuvjM+DvDMjRjiYO8/Jnj0Rd8df/DKX3LzEM6yhg+1Mnezk3Q+8pKUJ93w3fhvgzohw0Ap42BAHF56bYYmuYjLLXolzyCAKrr2UtMkTsGd+9s2QPim+wAOQOnYY/b55NWmTx2BL/WyWl3fRGXHrpo0acDBAdra8jXdusgO70870+88jd3wxKZmR/oYDIabff27CiwSzbhjL3G9MJGdg5OdbOCqX7P6JZWkAY2dkc89fJzB+djbbVrcxYlImqRmJz4U5HMJP/7eQMROcPP7fTSz8tjXBB+CWWzLJz7cRCChz51oz/BYR7vhGNulpgserDBpozZ9lYYGd//lpAWNHpZCZ0TMjtO44Fhd4ep+747DBOuDHd1iuK+nWV5oByMr1WK45PL/Jck2ACTkHkqI7JWN/9xfFwRhn9/2NZ6V6kN131Nfr6kLk5NhITY1dt8DW/bRKMKi0u5XcnNgDUSzujqqK26NkZsSumzZod5eujEOPy9G7X+zWsBGAb01YGpO7Y2+g12eWBsPRsHpLD0C/fslZCXY4hNwc6/srImRmWOvuGDBnww0Gg6E7emxz2ycwwdJgMPQIxZzgMRgMhpg4FjPLpH48iMi5IrJNRHaKyA+P8PpQEVkqIutEZKOInJ/M/hgMhsRRFcJqi+nRl0haZikiduBxIs5tFcAqEXm1w7oyyk+A51X1tyIyEVgMDE9WnwwGQ+JEFnjMcUcrmQ3sVNVdACLyLHAREevJDhToOAKRC1QlsT8Gg8ESjAeP1ZQAnTfYVQBzDrvmPuBtEfkWkAmceSQhEbkZuBnAnp+PzW39LyqcpK0QnhRrNi53pjk98Q3hR6IxLTM5us6spOg22605kXI4mZKMPbf+JGhCmsVViWIhssBj5iz/r7kaeFpVBwPnA38VkX/rk6o+qaozVXWmPTM5f9AGgyF2QthievQlkplZVgJDOn09ONrWmRuBcwFUdbmIpAFFQC0Gg+ELiSIms7SYVcAYERkhIk7gKuDVw67ZB8wHEJEJQBpQl8Q+GQwGCwhji+nRl0haZqmqQRH5JvAWYAeeUtVSEbkfWK2qrwJ3Ab8Xke8SmQq5QXvbYXWD4RhDFQLhvhUIYyGpm9JVdTGR7UCd237a6XkZMC+ZfTAYDNYSGYabYGkwGAzdYk7wGAwGQzd0bB2K5XE0RGSciKzv9GgVke+IyFQRWS4im0RkkYjkRK93isifou0bROS0TlpXRk8BlorIg4f9P1eISFn0tX90ar9eRHZEH9d393332cwy2NJKoLqakMtFqL2dsKudkMuFIy+X3PlnII6ef+uqSrCmjlBrG6FWF+E2F6G2yL+Zc2eROmJoXH1VVcLtHoINrYQaWwg2thJsbMGRl0P2/Jlxmasd1A4r/vo2PJVNeCqbCLl9DLxoOvbUxG0LPI0e6rc20LyrmbFfHkNqrhX+Lkr9Xg/bljcyeX4/cout2afq94X59B0Xw8elMnS0dXtft5UFaG4KM2eedZpej/Laa14uu9za/bRr1/lxOoXjJiX6u7dmGK6q24BpcPDEXyXwMvAi8D1V/UBEFgLfB+4FboreN1lEioE3RGQWkA88DMxQ1ToR+bOIzFfVd0VkDHAPME9Vm6L3ISIFwH8AM4nE/zXRE4ZdFovts8HSluqkcdHrBKoOHQrKnDmDnNNPiytQQqQuoH/Pfhqefu5QW1oqhTdcGXegBPCVV3Lgob8QbnMf6uvcKeRddGrcgdLf6GLbfy3Ctf0AYX9kk3XaoDwm/vyyuANl0Buk9Nkt1G2uo2FrA+21buypds78f2ckFCi97UF2rmhiy8cNbP2kkYb9Hi7+wRhLAuXeHT7efqGZd15uYfKsDE4+vyRhTYCKfUEef6SVZR/5eO39+LyNjsTmzQHu+HYzCxdau5/4xZc8/PS+VtaujM+B9HCS4MEzHyhX1b0iMhb4MNq+hMgi8b3AROA9AFWtFZFmDgW7HarasZPmHWAB8C6RAPt4RxBU1Y5tiecAS1S1EUBElhDZxvhMVx3sc8Ey5HbTvnYdrk9XHAqUdjsFl1xE9twT47ZYCNY34lq2ivZPVh5scwwopt83riNlQPxvQP++A7S9vwb1Bw625V95ZiRQJnA6w3ugBdCDgTJr3EAm/uelpOTFf+rFkeZAQ2H2fRg5mGVPtXPWI/MZNDt+J0JV5eNnKlj86C7CochGiJO/MphTrxvSzZ1HZ/2ydv76v3WUrY1Uqi8odvDtBwYkbLHR2BDi979u47m/thMMwPfvzSW7B5XNuyIUUp58op1H/p8Lh0O48KKe+zt1pfvzB9v49W/aueD8NNLSEg9ykdXwmM+GF4nI6k5fP6mqTx7huqs4FKhKiRyN/hdwOYf2a28ALhSRZ6JtM6L/vgeME5HhRE4KXgx0lKYfCyAinxDZlXOfqr7JkU8YHvWTtE8ES1XFW16O69OVuDduBJudzOnHU3DpJdT99W/0u/460oYP67luKIRn0xZcH63AW7Yde14OWSefQKitnVBzC4U3XBmXaZmGQrSv3krrW8vxbtlDypD+5C84g6Z/LqX49svInDWxx5oQ+Tk0rdxF5XMrad1cQc5xg0ntn0vG8ELG/ejL2NPjt5atK6tnzW/WUrWymox+6fha/Zz9v/MZODMxy9aWGh97N7QeDJSTTi/iknvGJhzUxk1Nx+065Bp518MDyclP7O3udod55IEWFv0zEoAHDbZz5VcTzwArKkLc+Z1mVqyIfGBefEka2dmJB+C2tjC3fquZJe9EbDLOPduaqYIebkqv785WIroP+0Iiw2WAhcBjInIvkb3ZHWdFnwImAKuBvcAyIBQdXt8KPAeEo+0dtqsOYAxwGpGDMR+KyORYO9+ZXh8s3aVlNP7rVYL19aSOGEHhZQvImDoFW2oqYb+fQXd+B3t2z5393Gs20vT8K4RaXaRPnkC/224gbdI4xGYjUF2Do3+/uIbIni27qf31C4SaXWTOmsDAe28kbcJwwu0e0qeMJnVYfMEn2OZl013P4N5TT8GJo5n8y6+QM7GEA4s30P+cyYg9/j++Ff+7itJnyhgwYwAX/PF8Gnc0kjs8l4HTB8StCbBrTTNPfGM9ef1TWfjYZJY8sYfrHj4OW4LGWsGA8sNr99JQG2TUpFQmz85k+kmJn0/PyLAxqMSBwwHBINx2Zw7OHvjwdEVOjjDpuJSDwfLqr1gzV1m+K4hGPy/sdjhzvjXZKlg+DD8PWKuqNQCquhU4GyA6JP9StD0IfLfjJhFZBmyPvrYIWBRtvxkIRS+rAFaoagDYLSLbiQTPSiIBtIPBwPtH62SvD5b2rEwyJk8ia84cnP0/Oxy2OZ3gjC+bsuflkHXKCWTOnYUjP+8zr6UMjH+OKqV/IdmnHE/OmbNwFB7StWdlYM9KYIicnUbR6RMonDeGjKGH/KEHnD81bs0OBs0ZSMkJgyg5YRAiQtHEQmwJBN8OSiZkc+5tIzj5miHYHMLwqbnWuDumCOd/JZ9pJ2ZQscvP5DnWFdw47aw0JhyXwicf+vjSxdYEtZwcG19bmEG/fjbefdfH9OmJL74BTJvq5Kc/yWbK5BQ2bAyQn2/N5pckFNK4mk5zhSJSHJ2TtBEp4/i7aHsGEZPFdhE5Cwh2lHzsdE8+cBtwRVTuX1H9P4lIEZFh+S6gHPh59HqIBOeOzPaI9D53xyFDdNBd37FcN5we7v6iOHDkWl9tZmBhi+WaAOPyknMkf0pWRVJ0x6Ump6JfiaO122uCQcXhiD1g5NpC3V9ExDmyJ4ZosVYdqqsP0a8odt3+g6u7dGUsmNBPz3pqQUw6z8994qjujiKSSeTY80hVbYm23QHcHr3kJeAeVdXonORbRIbalcCNqro3es8zQEdmcL+qPhttF+ARIos3IeCBTq8tBH4UvecBVf3T0b6XXp9ZGgyfBz0JlD0hWc6RPQmU3aEqBC06waOq7UDhYW2PAo8e4do9wLgudK7uol2BO6OPw197isg8aEyYYGkwGHrMsVh1yARLg8HQI47V4r8mWBoMhh5jgqXBYDB0w7Fa/NcES4PB0GOScNzxC0+vC5YSAmez9cWSgv7k/PKDEv+pma5oSE2OUVelMzcpunkp7u4vioNUW6D7i74gtNuT8zNIk9i2JFmJKgRN8V+DwWDoHjMMNxgMhm4wc5YGg8EQI2qCpcFgMHSPWeAxGAyGblA1c5YGg8EQA0LIrIYbDAZD9xyLc5bHxMdDOBjAV1NF2+Z1NC1/n7DfZ42uz4d3/z5aV63AvWO7JZoA4XY33i07cC1bjYas2UenqvgPNNHywUZCHmu+f4CQL0jtsl34mqzbR6iqVG9qoO2AtXsTGyq81Oy2VtPrDrNzs8dSTYCNa7yWazbUhdi/N5iwjlXujr2NPptZhjxuahf/E2/lPgLNjaCKPSOLkmtuwuaMv7x+e+km2lavwn/gAMHGBgAyp06j+MqvxK0ZbGrGvXI9/r0V+PdWEKxvxJaTRf87b0Hs8ZXWUlUC1Y24S/fgLt2Lu3QvwcZWBn3nUuzpidkLBNq81C3fQ+3H5dSv3MeQiydTPHdkQpoAfneQrW/uY+ML5WQWpXHxYyclrAngagrwxu/2UfphIz9+ZYYlmgDbNnp4+M4qfvq7wZZpqipP/qKFoF+ZMsO6yuaNDSG+cU0dTzzTL3ExjcxbHmv02WCpoRAaDhNoigQ0R24+g796C87CxNztHLl5eMrLUV/kkz979gkUXbIgIbtae24O3u3leDdvA8CWk03/u25JyAhN/UEaX19B8xurDrb1v+k8ck6Oy34koqlK+dMr2PW31WgoUiy5aPYwxt40N25NgPryFja8UM6W1/fidwVJy3Vy8WMnJezD43OHeO8vlbz9h/14XSFueGgcKU5rzMWe/10Df3+sjnFT0y2z1g0GlP+5p4HXX2zn0b9a48II0NIc5tZr62lrVQoKralraVbD+wCB1maaPllKy9rl2NIzSMkvRBwplFx7Myk5ed0LdIGvuormd5bQvnkjKcX9CdR6yT31NArOuyAxF8bt5TS/9Ab+PfuRVCe2tFSK70wwUKrSvmEX7Wt3Hmwruuo08s+bHbdmB+mDcun4O8kYkseUe89NyN8HwJ5iZ+/yGvyuyBBx/j3TyeqXmGVDc42Px7+xmYot7QCUjMtk1gWJB6AD+/08fFcVZWsiQ+9zr4j/PdUZd3uYH91ax4oPvTidMGWWNQG4rTXMbV+tY3tZgFPOtCZTVbPA07vxNzXQ9PG7tKxfhSMnl37nXEzOtNm0rPqEnGkzsafH58Lnq6qi6d23cW/eROqQoQz42tdxDhpE2+pV5J0+P+5A6d9fRfPLb+At3Uba5PEM+Ml3cH20guwz5pHSP/6hkr+qgZo/vkn7up3knj6N1PNnE6hpovDyU+LWBPA1uil75D1ql+1i2IJpVL5RxvH/dQEp2Yn9UQd9IVb8oYyWynYy+6Ux+Ph+jDsnMRtcgNxiJ+Nm5x0MlhffNQKbLfFsqGKXn7aWyDxyRpaNk8/PSVjT71MeuLucu3AJAAAgAElEQVSB1Z9ERitTZqWRlpZ4MGp3hbn9unrKNkbO0I+fZI23D5hheK/EX19D48fv0rpxLSkFhfT/8hXkTJ5+cK4v/8RT49L1VVfR9M5ng2T6uPGICKpK/hlnxqUbrGug+ZW3cK9aj3PUMIq/dytpY0ZE+nr5BYgjvl9J2OOn/sUPaVy0nLRhAxj23zeSPm4wwZZ27NkZCWW/B97fQdkvlpKSncrsX11O/nED6X/qaLKGFcStCdBW4+bVu5bRUtHOgsdPprmynTFnHNW6OSZUlX8+uIulf6/ihgfHsXpxHZNOzu/+xpi0oaLcz/HzMhk4LIW0jMSDmjNVuPnOPD58201+vo1Z86ybq1x4ezbf/XoDqakwfpJ1RV2OxdXwXh8s28o24q2uZOCl15A1cWpCc4edad+0kVBLy2eCZAeJBB7//ioCldUU3XYD6VMmfFY3zkAJoMEgrpVbGXDTl8idfzwSzaIcuYn7Wjeuq2DgmeMYc9NcHOmR7CR/8qCEdd0NPsQmXPP3M8ktyWSIasLzlAChgFKz28PXfzGB488uYsoZhZboAlTv83Pe1Xnc/OP+1FZaV/Woan+Q2Sel860f5xPwW5O2ZWbZsNngpNPTuPL6TEaOsSazVD02g2Wvc3dMGzREh918yHsoHAwidhsR18z4CWZ89uegwSDY7Qn/kQXzPrv1RyPvtISCeka/9iO2a1gPBsl4GFbQdGTdBIPYpNzqLl9LRHtMek1SdIen1B319Xi083pZibbjh1V06cqYPnqQjnzk5ph0yi7+2VHdHXsTvT6ztCWQjR2NRLK8o+qKgEVZzr9pWzAnd0TdJPU3mdq9sc+9iV6WY1lCrw+WBoPh/xZFCJvVcIPBYOieYzCxTO5xRxE5V0S2ichOEflhF9dcISJlIlIqIv9IZn8MBoMFRBd4Ynn0JZKWWYqIHXgcOAuoAFaJyKuqWtbpmjHAPcA8VW0SEeuOLRgMhuRxDKaWycwsZwM7VXWXqvqBZ4GLDrvmJuBxVW0CUNXaJPbHYDBYhMksraUE2N/p6wpgzmHXjAUQkU8AO3Cfqr55uJCI3AzcDJCSnU9ao/WdDXqT84v1Yc1Z3M64HYkdBeyKSns4KbpO2/+9A2Ei+MLWnXTpoMDhslwTIE2S5XBZ0eUrCoTDfSsQxsLnvcDjAMYApwGDgQ9FZLKqNne+SFWfBJ4EyOg/5BgcABgMXyAU6GNZYywkcxheCXQ+5Ds42taZCuBVVQ2o6m5gO5HgaTAYvsBEz1Z0+zgaIjJORNZ3erSKyHdEZKqILBeRTSKySERyotc7ReRP0fYNInJaJ60rRWRjdKH4wSP8XwtEREVkZqe2e6KLz9tE5JzuvudkBstVwBgRGSEiTuAq4NXDrvkXkawSESkiMizflcQ+GQwGK9AYH0eTUN2mqtNUdRowA3ADLwN/AH6oqpOjX38/estN0fsmE1k4fkREbCJSCDwMzFfVScAAEZnf8f+ISDZwB7CiU9tEIjFpEnAu8JvoonSXJC1YqmoQ+CbwFrAFeF5VS0XkfhG5MHrZW0CDiJQBS4Hvq2pDsvpkMBisILbFnR4u8MwHylV1L5Gk6cNo+xJgQfT5ROA9OLgY3AzMBEYCO1S145zqO53uAfhP4EGgc/n5i4BnVdUXHdXuJLIo3SVJ3WepqotVdayqjlLVB6JtP1XVV6PPVVXvVNWJqjpZVZ9NZn8MBoNFxJ5ZFonI6k6Prg6VXwU8E31eyqGdM5dzaDpvA3ChiDhEZASRbHQIkUA3TkSGi4gDuLjjHhGZDgxR1dcP+/+OtAB91JJXn/cCj8Fg6G0oaOyr4fXdFdKITtNdSGTPNcBC4DERuZfI1J0/2v4UMAFYDewFlgGh6B7tW4HngHC0fZREquv8Argh1s4ejWPqgKe/tZGmrasJeo5ctSceNBzCXbkH164tlmkChFwuXBs3RKofWYSGQnhKdxByWWzataeGQEObpZqeOhfuGms1/e0Bmvdbu4UnHFZqdln3fuqgerf1Jmje9hBNtVZtNZIYHzFxHrBWVWsAVHWrqp6tqjOIZJvl0fagqn43Os95EZBHZFEYVV2kqnNU9URgW7Q9GzgOeF9E9gAnAK9GF3liWYD+DH02s1RV/M11uCrLaa/aRXvVLgJtTZSctgBHnFXTO3R99dW079lB+76duPeVY8/IZMS13068v9XVuLeW0b6lDF/Ffvpf89WEqx+pKv7d+2n/dD3uVRvJnDOV9EmJbzjQUJi2ldtpeG0l4rAz7L74Dds642/zUv73ddSu2MvJv7/cEk2A/atreffna7nyT2dYpun3hvj7D0uZd9Vg+o9MvG5oBxs/bGbde01cf98IyzTDYeXX39/DdfdYZK5m7Qa+qzk0BEdEilW1NpoZ/gT4XbQ9g0hZyXYROQsIdpwI7HRPPnAbcIWqtgBFnXTfB76nqqtFxAP8Q0R+AQwisgtn5dE62WeDJeEwTdvXUrvy7YNNA044n8LJ8xKS9dZUUvnqX/A31QNgz8hi2BW34MjMjr+rfh81//gb7i0HT4JScO75ZB03JW5NVaXtnU9wvbuMYF1kF79z5FDyLjsvbk2AkMtD05L1NCxeRaC2BVuak9G/uiXhsmUhX5A9L21ix9/WEGjzMeNn52BzJL6hP+AN8cmvN7H27zuY+OVhpOdaUy28rcHPH765gertLr760HGWaAJsW93GY9/cweV3JW6t0ZkXf1XNyiUt3PG/FgVgi4KliGQSWdm+pVPz1SJye/T5S8Cfos+LgbdEJEwkC/xqp3seFZGp0ef3q+pRvamji83PA2VAELhdVY96eqLPBsu2/dto3r7u4NdFx59Gv5nzj3JH92g4hHt/OcH2yPBQUpwMvfwmnPlF3dzZjW4whC094+DXWdNnkndaYhmQiJA6Zjgt/4wciLJlZVB061cSzlR9++tpXLKOQG0LAP1vmI+zODHTLndVKyvufg3Xvkjx4dzxxQw8bVRCmgAHSht5896VNO6O/L6mXp64JkDN7nae/MYGGvZ7mHRaEQ4LHCMB9pa187+3bMPvDTN4jHWntJa/0cSLvz5Adr6DlFQL+mrhpnRVbQcKD2t7FHj0CNfuAcZ1oXN1DP/XaYd9/QDwQKx97XPB0tdcR9VHr9C2p4y8cTMonjGf9updDJz35YSyH3fFbqrffhF/Ux1Fc8/GvX8XhbNOJX1AYhlAe+lm6l5+EcRG7ryT8VVWULzg8oT6qqq0f7KGpr+/QsqwEgL7qij8+pU4ChILaqqKe2sF/qoGbJmppI0YQME5iftwpw/Mpt+swQeD5cRbTkw4Uw2HlPqdLbjqInN/xePzGHBcYp5BAAd2uvjt19fTUusDYPxJiWtCZI7yoYVbcbdFkpsSi4Ll7lI3j9+9B4CC/sawLBH6TLAM+X3Url5C/boPSC0cwKgF3yJz0AiCXjf542fE/ccXdLuoWbqIls2ryBo9iSGXLsSZV0hgUhMpOfGbYIVcLupfeRnXxvXkzD6BgvMvQH1eJCUloewv7PbS+NeXca/aSM4Fp5P75fm415aSPvmIH8ix6wZCVP32dZqXbmTgTecCkDV9VMLV2VWVst8sY8+/Spl856nUfLKbohmJz6vZ7IK73kvQG6L/pHwmXzLSkgrnA0ZnccaNQ1n0SDnhsDJ+XmH3N8VAS12AEZMy2fhRC1l5DnKLEg9sqsr6j1rJLUyhrtJvabDEnA3vfagqzdvXUv3JIjQUYtCpl1Iwcc5BjxtHWkY3Cl3ohsM0b/yUmg9ex+5MY8ilC8kec2huKt5Aqaq4Nqyj/tV/YUtNZeBN3yBjdHTBJT2xbMK3ax/1TzwLoSDF3/s6aeMjw87MWfHPfQIEW93s++8X8O6tYfh/fIWsaSMT9vuBQ4Fy94sbmXHf2Qw8dZQlw2+AXR9V8/HjmznzxzMYOKWQ3BJrFmAaKz28/stdnHPrCDLyHPQbFt/763CGT8pk16Z2LrxtEKGANcZtIsJZVxfxz8erufS2AZa6mYjJLHsflUtfoLFsBYWT59F/zrlxB8fDqXrjWVrK1lE053SKTjwTW4o1CwPNHyyl8c3F5M47mYJzzsXmTMx3uwN/ZQ01//M70iePp+BrC7BnWRMcVJV9//0CwcY2Rj20kNTBkflZK/x+qt8v/0ygBEjNS3z4GfSFWHL/aqYsGMmUBSMT1uvMu3/cS+HgNE7/2lDL5ioBVrzRQDiknLdwIOlZ1lWqWvl2MylOGxff0h+nFfOVENNRxr5Ir3N3zOg/RMdcecjd0V27H7HZSS9KzJo1eFiM9RzYjy3FSWph/4R0fQWf/fkG21oJNjSQNjz+VclA0b/vlVNVvKU7SJs0Ju6sJKfoyPsFvXtrceRn4ciJ74NoRP6Ra+ppWGnZXkfe+PhqPo/J7rr8acOuVvKGZGFP6XmAGJradQ3AoD9MS52PwpKeBfXuSrSpKnX7fRQP7ZlneHcl2lSVuko/xYN79qF8xZi1Xboypg4bogN/dEdMOnu/8X3j7vhFIaPY2i0WHSS6cNMVjuwcHNk5luuKCOnHjbVcFyBtWHIK2ItN4g6U3VE40vqfMYDDaetxoIwFEelxoIxZt4eBMiZ6V45lCb0+WBoMhs+B5NSJ/kJjgqXBYOgZx2jxXxMsDQZDjzGr4QaDwRALx2CwjGmpUETeFZHzD2t7MjldMhgMhi8esWaWI4AfiMgsVf1ZtO1z2Q5gC0BGrfWzy/7M5MzBSCgZuta7DwK0inWVczqzOymq4A9b75wJ0JxhzV7dzhQ4rS/jBpBh83d/UVysPeqrx+IwPNZNaM1Eyr73jxoI5SaxTwaD4YuMEjnuGMujDxFrZilRT53bROQG4GMg/oPRBoOhd3MMZpaxBsvfdTxR1adFZBNw+1GuNxgMfZhjcRgeU7BU1ScO+3oNEZ8Mg8FwLGKCpcFgMMTAMRgsjynDMg2Haa0tx+9psVQ36HXRWnXUKvY9RsNhXHu3Ew5Yu9rpO1BN0GWtEViwqYVgY7OlmiG3D3+ttZrhYAjPgVZLNQFcNdavdLtqrTcsCwXCuJsTfz+Jxv7oS/T5YBkOh2g+sI1dq15k7av3U7dnNSlpiRdZCPk9NOxYzc63fs/mF36OIzXLgt6Cv6WR2k/eZMfv/wt35W5LSsOpKu7t26j645M0Ln4Ne6Y1fQ27PTT/803qHn0ae078HkSH07Z2J7t+8CfsWdYVrPA3u1l917/QkLXbztb/cQN1m+os1Wzc28by31vrFgqw9PGttDf6rBEzq+F9B1XlwPYPqSh7h5A/8imdWTCEETMuTaiwqt/dQsXyl2mp2IKGIja1w066kozC+EvEqYZp3baBpk0raN+7A1AySkbQ74Qz49YE0GCQtg3raPnoA/zV1UiKkyF3fj/hwrIaCND23qe0vraUcLubotuuQSwwFwu5PFT/cQnN722g8MuzsWdYUy2ndWcd6378Go5MJxkliVlrdGb9UxtZ+8R6Lnv5Uss0fW0B/nXHMkaeMsAyTYDy5bUse2ons66yxrCsr2WNsdBng6WIkJk/+GBAS0nNYuy867HZE9vQnZKegyMt66Bu4djZFI6dnWBfbTgys3FXlAOKzZlGyfnXILbEApBnVzlN7y4h2NAAQMHZ55JSkJhnTKC2gbpfPk3wQCSbco4YQvqMxN0NW1dso+q3iwk2uUCg4EuzEtYEOPD+Djb/zxJC3iCjrk/s99SZDX/exNrfrcOZ4yR7kEWZekh57Z4VNO5pY8oC62xw2xt8vPLjiHmfM92ijfwmWPYNwqEg+ze/SfXWDygYPJmWA9sYM/erpGYkllUEPG3s/ehZ2qp2kDVgFCG/hyEnJJ5VeGorqXzzORwZOQTbWxl45gKcuYkbYYU9HkItLYjDQUr/AeTOOylhTUe/AtImjsYVDZZ5l5+buLlYIEiwxU2w1Q1A9swxpA5M/Puv+XAnpQ+/S8gb+WArPsmaqumb/rqZNY9HTrgUjiu0xAIC4ONfb2b3xwcAyOpnzRSEhpVX7l2Hqz4y/HZmWPAn3wfnI2OhzwXL9uYqyj99Bp+7iVFzrqJo2HRaDmwnpzgxb5eWiq3s/fAZ7M50xl3wbQDszjRsjsQy1abNK6l+50UyB4+i5EvX0ly6mryJiTkmqmrEvuKN18k79XQcubmkDR+B2BPLKlSV5ucX4/pgBfnXXIhn07aDPj+JIA477Zt2Y89IxVGQTeEF1mSAxSeNIvulDbTtqMORlUr26H4Jawa9QTIHZGJPtRPyhSgab41hWd32Zg6UHqrQntnPmkLA617aR1VpZKHM5pC4KscfERMsey8aDlO97QP2b36T7KIRjDvlxoOZZN7A+J0Nw6EgVatfp7b0QwrHzmbwnIuxpyQ+lxYOBjjw3ss0bfyUfieeTb8Tz0ZsNopmnpqQroZC1L/yMq0rP6XokgXknjAXDYUsC5Rt73xC0Te+QsaM48hI0Aitg7oXPqZ12RaG338t4kwhffRAS3T3vbyB5s3VnPDbK2ndXmtJBuhIc+BriWRpp9x3EvZUa4a1/cbmMWLeAGq3NjP6jBKyi63JLKdfNozqLc3sX99IZoF1FdPFFP/tnXhdDZSvfBZXw36GTvkSA8aehEjin6Ce5hr2vP83/K4mRpx+HfkjplrQ28iK9/5XnybQ0sDQS79O9siJluiGfV4O/O2vePfsYsANN5I5fgJAUgIlgD0n8bm6lk/KqP37+5R868tkThqWsF4H7fub2P7kMkZ/7QRyxvQjZ0ziWSVAwBNg/R83MvGKCYw+fxRhi1bXNaysf2EXx108nFO+PdkyX+5QIEzZ21WcfPNYpi+w7ud7LNLrg2XtrpXsWfcK6dlFTD77O2TkWrOKWL91OftXvEJm0RAmXHwXzixrjsK7dm+l4vW/kZKTz8hr78SZZ80wLtjSQvWf/kCo3UXJN24ntSRx723oOlBagWdHFRW/fIWii08k/8xplulqKMym/15C9shChl853TJdgLLnthLyBZlyXeTnYLNbM6zdvewALRXtTLt8FDaHdTv6ypfX4WkNMOmcEmvmKzsww/Deh9dVz8CxJ1My8Uxsduu+HW9bPQOmzmfAlPkHPcitIOhxkT1mMgPnX5rwfOdn0DD27GwG3rAQR571NU6sDpQAOGzknT6F/tedYamshsLkTxnE4PMnWRp4AArHFTDnztmk5lprApY7KJPT7ppC3hBrVtY7KBiSybk/OI7sYgvN0I7RBZ5eZ4WbVTBEJ5/9Hct1k1XP0ldgva63X3J+Z4F+R7dVjZecwuTUcizJtfYk1kHdDOt1e1s9y/unvNqlhW3aoCE6/OY7j/TSv7HtZ3caK1yDwXAM07tyLEvo88cdDQaDtQiR1fBYHkfVERknIus7PVpF5DsiMlVElovIpmix8Zzo9U4R+VO0fYOInNZJ60oR2SgipSLyYKf2O0WkLPrauyIyrNNr14vIjujj+u6+76QGSxE5V0S2ichOEfnhUa5bICIqIn0iXTcY+jQWFdJQ1W2qOk1VpwEzADfwMvAH4IeqOjn69fejt9wUvW8ycBbwiIjYRKQQeBiYr6qTgAEiMj96zzpgpqpOAV4EHgIQkQLgP4A5wGzgP0TkqJP9SQuWImIHHgfOAyYCV4vIv+2REZFs4A5gRbL6YjAYLEZjfMTOfKBcVfcCY4EPo+1LgAXR5xOB9wBUtZaI3c1MYCSwQ1U7Kpq803GPqi5VVXe0/VOgY5vIOcASVW1U1abo/3Pu0TqYzMxyNrBTVXepqh94FrjoCNf9J/Ag4E1iXwwGg5XEHiyLRGR1p8fNXSheBTwTfV7KoVhxOTAk+nwDcKGIOERkBJFsdAiwExgnIsNFxAFc3OmeztwIvBF9XgLs7/RaRbStS5K5wHOkzszpfIGITAeGqOrrIvJ9uiD6A74ZIDU1l4xqi8pMdSIlKzk/CnsgCQ6EmqzPuCS5RmpyXCODoeT8HFx+a7cFAeSlWV+fEiDNnpwdDN3Rg61D9d2thouIE7gQuCfatBB4TETuBV4FOpb8nwImAKuBvcAyIKSqTSJyK/AcEI62f+YcrohcSyQLjfuI3Oe2Gi6RIza/AG7o7lpVfRJ4EiAnu+QYXIczGL5gWPtXeB6wVlVrAFR1K3A2gIiMBb4UbQ8C3+24SUSWAdujry0CFkXbbwZCna47E/gxcKqqdmRalcBpnfowGHj/aJ1M5jC8ks+mwoOjbR1kA8cB74vIHuAE4FWzyGMwfMFRa1bDO3E1h4bgiEhx9F8b8BOihokikiESMbcXkbOAoKqWHXZPPnAbkUUiROR44Angwug8ZwdvAWeLSH70nrOjbV2SzMxyFTAmOrdQSWRO4isdL6pqC1DU8bWIvA98T1VXJ7FPBoPBCizKLKPB7yzglk7NV4tIh3vsS8Cfos+LgbdEJEwkpny10z2PikhH8Yb7VbXD5+VhIAt4IVpIZZ+qXqiqjSLyn0TiVMc9jRyFpAVLVQ2KyDeJRGs78JSqlorI/cBqVX01Wf+3wWBILlYdd1TVdqDwsLZHgUePcO0e4IglxFT16i7au7QbUNWniMyDxkRS5yxVdTGw+LC2n3Zx7WnJ7EsHXm8zIjZSUxP34elANYyrpYrsPGuKV3TgaT6AM6sAuyNxH54Ogu1tKEpKpnXff9jnI+zx4MizzrJBgyFCbS4c+bnWaaoSam7HkW/t+etAq4eUHOv8ggAC7X5SMq37vUP0+/cGcaRbsJB3DK4cHBMneMIaor5hCxs2/4XSrc+T4rRmdTbgb6di1wesev8h/D7rXANddXvYsfRPVJcutSxQhvxeaj95k33/+iOOdOtWp9u3bqHyt7/Glm5doYZgUwu1/+/3lv5BaihMzZOLCdQ2WScKNC7fScMnOyzVDLT72fGU9duO9y7ehmu/BefeY9021McCap8+G+71NlN1YA3VB9bg87dgtzmZNeOb2CSx7Tyulkqq9i6jrmo94XCQ4kHHU9g/sZqUqkpL1VYObF5KW+0u7ClpTL7oBwlpQqR4cdOG5dR9uoSQx8WQixYm7O0DEGxrpeHVV3BtXE/eqadjS7UmWHpKd9Dw5LOklPTHUWBNVqmBEFWPvYxr9Xb633jUfcc9om1LFdt+vogJ91tnWKaqbHroPcSi0m8duGtdbHh0Oac+fmHCWsKxWXWoTwdLf8DFgdp1+PyRT9Mxo84nIz2x+pGhoJ/aqvXUVETWoVJSsxk5MfE3YMPutexb9a+DTpSDjz+PlPTE7GV9DTXse+Up/E2Rgw0Zg0aQPWpSQpqqStuqlTQsXkTY4wGbjdy5iXv7aDhMy6J3aV30HqiSOcea+pZhr5/Kh56nfX05GVNHWuJCCeCpaKTs3pcI+4Kkl1hXEm/fq6VUL93J0IutK4enqqz9nw8ItvuxoCY2cGwGyz47DA+Hg1QfWIvX24TDkUZRwXgGDkh8V1Iw6KGlcddBl8gxx11KijMjYd2U9GxUw4CQUTiYfmNOTFjTWVBMxuBDJl39T7nAEhtcDYcigRLImnq8JXOVrveW43p3OaiC3W6JY2SozcP+n/2V9vXlAGRYVInd3+ii9EcvEmz1YHM6cBZZ45nesqOOsl99BGDpfOWe17ZRs6ICwDJzNTMM7yP4fK1s3vIsrvYDTJl0Le3uOgb2n57wG6WteT9la/6MIyWT6Sd9h4rdHyU8/IZIVrl72XMUjpiOIy2bgmFTEi44rKrUfLCI5tJVFJ90Pp7qfWSUJG6vKjYbrg3rceQXgE3IOyUxz6AOMmZNoeX1pdiyMkkdPRR7VuIfQBoOkzVnAp5tkUCRMWl4wpoAde9tQWyR91LaoLyDzxMh6Paz4YF3CAcie6kdmdacInIfaGPjY8sONVjQV6DPBcJY6HPBsqV1H5vK/oHDkcas428jI6OIwoJxCXvy1FatZ8fGF8grGsu4aVfhcKQy2oLh94Gy99m/5jUGTj6TkqnnEA4FEl7U6QiUDes+ZMgF15EzZgoBV+IT+6pK7YvP4a+uouTWbxH2+0gddNTjtLHphsM0/PF5bBnpFH/3awQqa7u/KQbsWem0friRrFnjSBmQT9roQZboDvjyNCqfW0m/MychdmuCj83pYNZDX+a9K/5M0YzBpGRZk1kqMPn2E1j3/z4ib2yRNZnlMVopvc8ES1Wlqnol28tfp6hwPBPGLsDhiHw6JxIoVcPs3f42+8uXMmTU6Qwbe/ZBvUQWSlTD7F/zGjVbP2LYnEspHjsXICmBEiAlK/HFkqYlb+HauIFBN96Ms3//hPU6aHv7I7zbdjHgx7fjKCrAUZS4ZzhA89tr8FfWU3L3laT0y7Vs0aT27c2EvH5G3HIaNqc1f0I2h42qd7fjzE5l1kMXEGizpv5B5oBsmrfXU3hcf+Y9cj6WOaGZYNk7CYUDbN+xiOqatYwcfhbDhpxiySdoKOhj24ZnaazbzripV1FccrwFvY2sUO/+5FmaKjYz+pTryR9qzWR+V4HSClpXraTp3SUUX/kV0keNtkzXt3s/zS+9Rf4VX8I51JrMDyDY6qbumfcouGguzgHWLcBoKEzVi6vof+4UUnITnyo4qKtK5ZtbGXTmWGwOO6n51mhrWKn+eC+jr5xs6TyoscLthXi9zWza8g88nkamHnc9hQVjrNF1N1K25i/4/S6mnHALOXlDLdEN+j3sfP9p3E1VjDvzFrKLE59HhOQGSveO7dS99AL5Z51L9vQZlumGPV4annyWtEljyZo/1zJdgPq/v4stzUnhgpOt1f1oO96aVgYtsLaEQevOetp2NzLlnvndX9wDmrbU4m1wM/Aka21wzTC8F7Jl+z/RcIhZx99Gero1wzdVZduG5xCxcfzcb5Kabt3JlMr1b+JzNTDh3G+SnmvdUNZdUZ6UQKnhMPWLXiH7+Bnkz+/y5FhcuFdvQn1+ChdeZt0qLRD2+HGt3Unx187Blmpt2bn2nTX0O2MCaQOsO1kE4K5qJX/KQHLHFdG/1HAAACAASURBVFuq62v2UjyzhJxhFjp+9sGV7ljode6OOdklOmv67Qe/9nqbSUnJwG5PbIgROKye5f9v77yj47qPe/+ZbcCi98Lem0iRYpGoYkmWqB6ruUiynSi2npzYjqPYxzlxTspL/I7Pe7Ff/KIkjm05kizHshot2ZQsiaIaq9g7wAaikCCISnRg6533x70gYYYkFtrfkgB4P+fcw927d4eDC+zsr8zMN9Tfjj+QmbTdUOHvr2vGoyHi0TCBjI//YesrOffaW7i9hbT84o9tN1Ry7r+FWHcX3mAG4vt4363RovP3XIx39+DN/njlhxn55+8RaYUiSJr/YwXhwqy+C76ucWvY65+J9LNU1WH7m6p+lr++/qfnVWXMKJ6ocx5MTN1x11OuuuOIIT3d3Kjv9+wGzWtvA3j96Xj9BjWcB5FMoLwQvmxzdeRn83ED5VB40s3WVQ/GdHXNabsGR9epxK3gcXFxcUkQsS6/aOkGSxcXl+Fxma5ZusHSxcVl2LjTcBcXF5dEcIPlyEeicfwnO4zb9WWZbd46gCeaCrupUWE01pLmv5Eaf/tiqfE3Ejbvb0dGijb1fPGhL0oB7sjSxcXFJRHcYOni4uIyBOqWO7q4uLgMyeWaZzlmm/+eD0stovGQcbvhaK9xm7FYiHg8YtSmZcWJhS9coTJcVJVYv/mf3wqZ/z1ZEfMVLxqLmbcZT83QTU3lR6omdowhLptgGbMi1HXsZFv9i0n3thxMT6iFfbWv0di+z5hNy4pRf2wjB/e/gsdjbrOhu7mGw+/+1IgGzwCx3m6O/+YZNGYuCKll0frGKqKtLcZsAvRu2034oGFxsaZTdK/dZdSmqtK2yrxgWV9VI32HG4zYEk3sGEuM+WAZifdR1baRdTVPcbDlPabkL8XnSb4Urru/iT01K9l04Md09jUwsWhZ0jZVLZpO7mLbph9y9NAbTJh8g5ESuGioh5pNL3Fw9Y/ILp2O12+mC3f30Uqqfv4DNB7Hn22m7NQKh2j8xbP07N5JwEBj4QH6duyl7ZkX8U801wbO6gvR+INfIoZ6Wg7QvnonPburjdrUuEXdv/4OK2pgFOyqO4492vvr2dXwGlHLns4VZ0yjJDO5Fm6WWhyqX83x1m2nz82ZcCeeJEdrnR21VB1cRU/3SQAKi+aSl59c+zZVi9aqrRzf+SbxSB8er5+S2cmLi1nRCE3r3uDUrg0A5C24JmmbANH2dhqfe5rIyZNkL7s6aWmNAfp27qP1P3+FJzMDb54hxUjLovnfXyF6otmYCiVAuL6VxqffIWvRtKEvHgbNv9lK/9EmY9P7y3GDZ0yPLDMDBWT47YYYHvExt2RF0iM1j3gozZt7eipfkjuHopzpSfuanTMBn3+g4aswdebtSduMR/oJ95wiHrHXKItmXoM/PTnNcFWLls1r6Ni/FQBvMJPsGckpRgKEjtVx4t+fJHLS/rLInJu8thE4gfJnz4NlEZg8wVizilMvvEPfrkMA+ArMNBqxonHqf/gbNBJD/ObGMeGT7TQ8v875T8wM98RK7BhLjNmRZV+0gx0nVgLK1Pyr8XnSCfqTHwG0dB5mT81KSnPn0hdpZ/b45IOaqlJT9Q6d7TWUlC3E4/GRmVWWtF0RL50nDpKeW0q0r5OyuTcasOkhZ8YC2rZ/iDc9SN68pXi8yf8Z+YtLCM6aRc/OHba644xZSdu0+kOEq+tOPw9MNjOt716/m653zqwpeg0Fy5YX1xI6an9ZmJraqyrH/v0tNGxPv42MLJUxt3mTCGMyWHaGGtnZ8GuC/lwWlz8IAj5P8ut0Daf2UlH3WyYULWHOhLsIRTsJBpJbq1NVqo+8xYljG5m74BHyC6YTiyevv2LFY1StfY5ouJe5d36DUGczaVnJN0eO9fVwfNVzZE2ZQ+GSm/BmmGmxFutop2f3LgrvvZ9oawuetOR/X55gOoEJ5QBkXr+MwOQJSdsEyLphIeEjx+nbcxjxePCkJb8GrnGLjLkT8RfloHHLmLZPpLGDzLkT6N5diycjzdw0/PKLlWMvWLb0VrPn5CoKMiaxsOxTeA3tJtc1b+HQidVMK7uR6WU3ISLGA2Vxqa3F4/MnVyKpalGz6SV6T9Uz5/avkZaZR1pm8hswalnUv/k8eDyMv+vzeNPNlHJqPE7LKy+RPmUquddeb2zUotEYHaveIesT15D/0L2oobQh7Q/TvW4XhV+4g+ACM3pE4vUQKC8g2trFlP/9KB5D0/C08nzSJxUhPi9z/+0xrJChVDQDvyIRmQ28NOjUNODvgQ+AnwBZQC3wBVXtEpEA8FNgKWABT6jqh46th4C/AbzAG6r6V875NOAXwBKgDXhIVWud1/4aeAyIA3+uqqsv5O+YCpb1nfuobF7N+JwrmVuyAo+BFCFV5WjjWqob1zF7/B1MLjGzmXG+QGmC+p2/o/3YPmaveJyM/HJjdls2raav/ihTP/+EsUAJ0LF+LZHmJiZ+89vGNnUAetZvxuruIfeeFYjXiwTNpEx1r9sFHiHrE4vwpJvJLADo2liJrzCbjDkTjWiRn7a7o5qsKyaSVmYmY8FUUrqqHgIWAYiIFzgBvAasBL6tqmtF5MvAXwJ/BzzuvG+BiJQAb4nIMiAf+AGwRFVbROQ5EblVVd/DDobtqjpDRB4G/gl4SETmAQ8DVwDjgHdFZJaqnrfYfkxs8KgqVW0bqWh+mxmFNzCv5DZjgfJg/dvUNK5n/uT7R0WgbDywjsbKdUy74fNklya/8TRA99FKWjavYdyKzxIsMZfSE2lppn3Nagpuvwt/YZExu1YoROeb75F96yfw5mYbs6uWRdc7W8i+abHRQAnQubGSnOvmGg2UaildO46Ss8Tg7roqYiV2DINbgaOqWgfMApwdKdYAn3YezwPet13QZqADe5Q5DTiiqgOJue8Oes99wHPO45XArWLv8t0HvKiqYVWtAaqAqy/k4KgPlpZaVDS/Q/WpzcwvvYtpBcuN7HhaVpx9da9xom0Xi6Y9xLgCMyJgqQyUp2p3c3z760xadh8Fk82JlkU62qh/83nyr1xO3vzk80kHUMuiZeVLBMrHkXuDWRXGrjXrwFJy7rjJqN3+/UeJnmwl53YzX5wDhOtbCdc2k3u9mSyAAfprmoi195Kz1NwXJ5CKPMuHgRecxxXYwQzgs8BE5/Ee4F4R8YnIVOyp9UTsQDdbRKaIiA+4f9B7xgPHAVQ1BnQChYPPO9Q7587LqA+Wextfp7H7AIvHPcj4HHP623trV9LaeYQlM75AcW7yO7MD1FW/m5JA2dNSR/XGFyi74pOUzkk+l3IAteIcf/05AnlFlN3ygDG7AF1bNxM6fpySzz5kdPpt9fXTvWYdOXd9Ek/QbIu87ve2E7xyJoFyc6NggO4dVfgKcwjONrMJNUDvgRP4i3NIn2TW32FU8BSJyPZBx1f+my17LfJe4BXn1JeBr4nIDiAbGFhofQY7qG0H/gXYBMRVtR34Kvb653rsdU7jvetG/Zrl+Jz5TMtfTk66OVlZEWFcwUKmld1ETkbyKTyDKSy+gqzscRSVJJ+bOJhgXhkTFt9D6RyzIzTxeCm6+laCZRPx+Mz2ecxasBBfTg6BUrP32JMRpPjrf0xgqlmtbIDiP3mAeLfZ2nqAwnuvIfeGeUan4ADFf7CE/JvmmRVDU4aTr9magLrjXcBOVW0CUNWDwO0AIjILuMc5HwO+OfAmEdkEHHZeex143Tn/Fc4EyxPYo8x6Z9SZi73RM3B+gAnOufMy6oNlcabh6YVDSd6clNjNzhlHdo65krsBvP40I3mU5yJ39sKU2PVmZpI5z9zoejDps83sUp+NJyMdTwoa+YoI/sLUqGj6slPQgNps6tAjnJmCIyIlqtosduXH32LvjCMiGdjy3b0ichsQU9XKs96TD3wN+JxjbhXwKPAR8BngfVVVEVkF/EpEfoi9wTMT2HohJ1M6DReRO0XkkIhUich3zvH6t0SkUkT2ish7ImJ+KODi4mIcU400RCQTuA14ddDpR0TkMHAQaACedc6XADtF5ADwV8AfDnrPkyJSCWwE/o+qHnbOPw0UikgV8C3gOwCqWgG8DFQCbwNfv9BOOKRwZOmkAvwI+0bUA9tEZNXAN4HDLmCpqvaJyFeB7wMPpconFxcXM5iSwlXVXuwNl8HnngSePMe1tcDs89h55DznQ9ibROd67XvA9xL1NZUjy6uBKlWtVtUI8CJndrgAUNUPVHVgAWgz9rqBi4vLSMbtOmScc23NXyjf4jHgrXO94CzYfgUgXTLRplZTPp7G028uF28wZjPxbCx/ar7j4n6zmwsDqDdFdjHXl3MwMU2+fPFseqKp+Z31+y9+two7KX2MRcIEGBEbPCLyRezk0nMmxanqU8BTALneosvvt+TiMtIYYx2FEiGVwTKhrXkRWYFd03mTqibfQcLFxSXlXI4jy1SuWW4DZorIVCfp9GHsbfzTiMhV2IXx9zrlSy4uLiMdd83SLKoaE5E/A1ZjdwJ5RlUrROS7wHZVXYVd/J4FvOIkzR5T1XtT5dMg38wm6abIZqrsqlpGdYhsm6Pnnroky7DrvscEKc2zVNU3VXWWqk53tulR1b93AiWqukJVS1V1kXOkNFCGrX4Oh7YRVnMVGKpKS38tdT27jdkECEU6OVi/GuvCqV/DworHqDv8DpFQlzGbAG1HttF/6oLFD8Omr76GvnqzOjSxrk76Dh8yalMti/4KszYBQkfqhr5omFjhCJEGQyJwrrrj2CRihTgc2sb6npVYWKR7kpNWGKAj0si21lfZ2fYGpUEzlUSWFaO6cT0bDvwH6f4cvB4zg//uznp2bfw3eroaSAsaEheLRanb8DINO98mWGCuKqn7yH7qXv4p6cXm2svFe3s4+bOf4vGbLdns/O3bROtPGrUZbWqjY+UaozYBOl77kFjjqeQNqSsrMeaIapi6cAV1kQrixPBLGtPTFiVttyfazpGuTTT1HwVgStZVBH3Jl6q1dB7mYP1q+iPt+H0ZTChakrRNKx7jWNV7HK/+ENRi2tx7krYJEO5qo/qD5+hvO0HJgpuNTevbd3/EyXdWkl46Hm96xtBvSAArFOLk0z8j2tyEv8xcAO7bsZeu1R9S8Ojnhr44QdSyaHvKlkMxSeR4Ex1vbKD0CUOpzGNs1JgIYzpYxjRKt9VOHFt/ZHraIvySXOajqtIVbaItZKeQ+iTAtOyh+gQMTXdfI9WN6+mPtAMwuXg5Pm9y+X7xWISDu1/gVLNdNBXMLCKvKPma6Y5jFdSte4F4pB+AgqlXJW1TVWnd9A4tG+xm1ZmTk1PhHMCKRmn8+dNETtTjzc3Fa6gLUeREI22/eBkAb7YZaQ2Arrc2ED5SR9rsKcZsqmXR8vRvIW6hJqRwYcxt3iTCmJ6G91odtMVOkOnJJdOTywS/meYYXZEW4hojzZPJ1OwlBLzJfwDTAjnErAg+bzo+bzqTipPvG+n1BSifdA1OGjHlk68zMgIM5pXiz7DF39JyiggWJt8MONxykv7G+tPPMycn3xZP43GaX/gloRp77TNgaFRp9fXT+pPn0LDdOcybYyZYRo430vFre/otXnMfze4PdhA+dAwAjZlZAxfLSugYS4zZkWVTtI69/R8y3j+TGWmL6Yy3Jt09XVU51LmBup49LCy8E0tjlKYnv1YZi0fYdfQFLCvG8tmP09J5CJ83+dqf3u4mDu7+FWUTr8bnD1I6fnHSNgG6T1YR6mxi3NJ7sGIRIzvW6SXjEBHSCktRVTImJKeZDiBeLyWfe4Rj//Q9PIE0AmVmWsHFWk+RPn8OPR9sBI8HT3by1V+qSu9He/Dm5RBvbQevmeqkeFcv3R/sABFQNTOyVNyk9LHCyWg1+/vXMSkwj1lpyxARij3JrdWcHSjLgmZagFlWjD01LxOKdLJs1pfISMtnUnHyXbijkV4qd/ycrJzxTL/iPnts6Un+AxjqaKJ+y28pX3QbZVfeQjwSStom2LIV3Uf2M/nzX8efnYfHb6bksHPTBlBl/DeeIHrKwOYGEJg0nvhb75E+ZyaZ1y7Bm538hqGIkPeZ2+jduIuce27Ek2lmucCTnUHJEw9x/M//mfyHbkPSkt/gEvSyTEofc8HyROQwFaGNTAssZHraVUZGPakKlKoW++p+Q2dfA8tmPkpGWj5A0j5bVpwDu54HYO7iL+IxECTB3iyqWfs8wYJxlC1cAYA3kHxvRysaoXHNq+TOX0rmRHP9Sa1QiM71a8m98Wa8Wdl4s8zU/8c7u+jfU0nRY58nY4lB+Y6aE8Tbu8i6YTH+8SVGbIoI/XurkGAaeX9wA5ia3l+GwXJMrVkei1RSEdrIjLTFzEhfPMIDpXLg+Ju0dB5m8bRHyA6a6/ReXbmKns565i15FH/ATJoUQMOON4l0tzHl5i8YGaUO0Lr5PeKhfkpv/pQxm3BmVJl7nTmZDYCej7bjyQwSXGhWL6dvRyW+8iJjgXKA/n1VBK+Yhvi85pL83TzL0UtNeB8HQ1uYnXY109LMdPZOVaAEqDr5ASfadrNo2ufIy5o49BsSpKHuI04e28LshY+QmW1OrqGr/iDN+9cy8brPkJZVYMxu+FQzbVvep+TGu/Flmuv8dHpU+Ymb8KSb62yulkXPhq1kXrsM8ZmbmKkqfdsryFhqVm5ELYv+/dXG9M1to9hrlokcY4hRPw1XVY6Gd1Md2c289OuYEDhnb9CPZTdVgbK2+SNqmjawYMqDFOWYs9vRWsXRylVMnXMXhaVzjdmN9ndTu/5FCmYuo2Ba8mlCA6gqje+8SlpxOfmLrjVmF1I3qgwdrCLeeoqsGy6omjpsog3NxBpbyVhiNliGqxuwevvJWGBWfmWs7XQnwqgfWR4Jb6c6sof56TcaC5RAygLlibbdHD6xhrkT7qY835z+TH9vKwd2/ZKScYsYP9WcFo+qUrf+Jbz+NCYuv9+YXYCug3vorTtC+e2fMavuOHhUaVjdsWf9FtLmzMBfYlYtsW97Bd6CXAJTzWmygz0F9xXn4SsrHPrihElwCj7GpuGjfmSZ7snkyuDNlPmnGLWb5s00HigBfN40ZpbfwsTi5BPZfw8RCkrmMXP+A4aV/JSMwvGUX3U7Xr9Zoa5AfhGln/wUwfJJRu3i9ZK/4nayFxu+x0DuPSvQuKHE7kHk3HYdwStnG28ckvcHN5C5fL55dccxFggTYdQHy0kBs4vsA0zNNpOTeDaleeamx4MJZhQye6G5srsBxONh3JK7jNsFCJZNIFhmXknE4/eTe71ZSeABAhPMlUsOxpORTprhUSWA+H3GNc6BMbcemQijPli6uLhcfNw8SxcXF5dEcIOli4uLyxCoQvzym4ePumCploXV3X2p3UgYCZrXd/T3pEIzEgK9KVJLzExNt3NfemrsWgHzSSKWpOajFrPMNYceFu7I0sXFxSUB3GDp4uLiMgQKXIYaPG6wdHFxGSYKevmtWY76Cp7hEtYQYe03alNV6Yq3GbUJ0BNpI26ZTYCORHoIh80KlsVjEUK9ZtqfDaBqEeluN2oTINZrfr3bCplpUzcYYx3Nz7YbN7DGqdgbPIkcY4jLJljGNMpRrWAPG/FjboOkK97Ktr636Imb+2DHrSiHW9dxpHW9McEygJa2A+ze9yx+v7lORKGeNirf/xEerzkhMI3HOf7O88T6zAa2zp2b6as+YtRmtKWVrvUbjNpUVTp/965RmwCR+kZC+wz9/JdhueOYD5aWWhzTKjbyNjUcYApzk+6YDhCyetnXv47Nva8T1RDl/mkGvIWW3mo21j1LTfsWpuSbKdeLxcIcOPwq+yp+SWHBbGP9LTsbD7N/zZN4fGkE0s10DLKiEWp/9wxdNRUEi81VtHQf2EvTG68QKDLX/syKRGh+9jmj7eoA+rbspH9vpVGbalm0/+K3WI4URvIGkw+WIjJbRHYPOrpE5C9EZKGIfCQi+0TkdRHJca4PiMizzvk9InLzIFuPOOf3isjbIlLknF8kIpsd+9tF5GrnvIjIv4pIlfOeIUv2xvSaZa92sYdN9NEDQC6FFJNcuZqqcixSyZHwTqxBQmjJatuEot0cbH2fpp7DAOSljyMvPflg0dlZR8WhlYRCpwBhXFnyAVhVOXloLcf2/g5UKZiwIGmbAPFwP7VvPE1vQzVZk2YbC0J9NUdo/PV/gWIsWKoqbS+vJHryJJ4scyP1WHsHp178Lb4ic23wAHrX7yB8uJbMG018AZsZNarqIWARgIh4gRPAa8BK4NuqulZEvgz8JfB3wOPO+xaISAnwlogswx70PQnMU9VWEfk+8GfAPwDfB/5RVd8Skbud5zcDdwEzneMa4MfOv+dlTI8sM8gmn+LTz2eyIOmGAiJCmX8qPidvLsuTT6kveb0YnydAwHvmQzcl/2oDHdNjdHYfJxzuBKAgfwbBYPIfwpba7TQc/OD0B8ZEsIz2dXP01f+gt8EWF8saZ6alWOhkPQ0vPoPG4/jy8o3JVXRv3ETvjp0AeDPNBEtV5dQvVqL9IaONL+JdPXS8/Jb9xIRgmQKWldiROLcCR1W1DpgFrHPOrwE+7TyeB7wPoKrNQAewlAFFPsgU+8blAA2DvB3Qqc4ddP4+4BdqsxnIE5ELjqTGbLBUVY6wlwZqmcF8iignT5JvKBDVCLv63sOLnzLfVGPSFT3RNk507aMsay6ZgUJKMpPvduTx+IhG+xDxkJlZyrjy5BUjAQonLsQXyMCflkVW4WTSMvKStunx+Smc7/S0FCFzfPLLGqoWfTVVeNLtNm2mRpXhY8fpfPf90889WWbUHXs3bCVUedgxau6j2fHy21i99qamOSlc42uWDwMvOI8rsIMZwGeBge7Ye4B7RcQnIlOBJcBEVY0CXwX2YQfDecDTznv+AviBiBwH/i/w18758cDxQf9/vXPuvIzJYDkQKI9TxQKWM0XmMJfkuwjFNMquvjVEtJ+lmXcwK30pJb7k24v1RtrZ2fAqRRlTubLsbhaV32ckADe3VlB3fC1zZj3AnJkPUFRgRgr42J43iIV7mX/bE5TPvsmITY/PT+veDeRMX0DpNXcSLE3+vop4yJ67gFh3F/nX30J6uZkOR2mTJpJ31x2I30fa9GnGRpb+8lLS5szAW1SA+M2skFm9/QSmjAOfD09OlqFgqcPZDS9y1goHjq+cbU1EAsC9wCvOqS8DXxORHUA2MLDQ+gx2UNsO/AuwCYiLiB87WF4FjAP2ciYofhX4pqpOBL7JmSA6bMbcmuXZgbJE7C+LNEmuCWxcY+zue48+q5tlGXcR9JjZ0AjHetnRsJJMfwFXlt2DiIesQPKNWnt7mzlwaCUTx19HWYkZmQ2AjsZDNFVtYtb1j5KWmU9aZr4Ru237NhHpbGXqp/4H/ux8Y9PQUxvfx59fSNEtd6MGSwN7Nm8hY+FCCj/9gDHZ2sDkCURqj5P36XtIn21mGcKTGSRt5hSIxSj5yy/bk9JkUXvUniCtqjrUQuldwE5VbQJQ1YPA7QAiMgu4xzkfww54OK9tAg7jrHuq6lHn/MvAd5zLHgWecB6/Avyn8/gEZ0asABOcc+dlTI0szxcok8XSOHv6P6DbOsXSjDvI9OYasRuzIuxseBUPHq4a9wBej5n0m1gsxN7KX5KdNZ7pU+80YhMgFu7j6NaXKJqyxNimDkCsv4emLW9TvOQWAjkFxgJltKuDrt1bKfjErYjHg8dn5v5GGhsJ19aRvfxqPOnpePxm7IYOHUXDETIWzsNfWjz0GxK1W3EEb14O/gllBCYZ6sdpaWJHYjzCmSk4zuYNYu+a/i3wE+d5hohkOo9vA2KqWokd5OaJyMBNuw044DxuAAamP7cAA7lTq4A/cnbFlwOdqnryQk6OmZFl6gKlxd7+tXTEmlmaeSdZXjMjKUst9px8nVCsm+UTv0DAa0b+QNWi8uArWPEY8+c+bCxNCKBm56uIeJhylVl5icaP3sQTSKdk8S1G7bZv+gBvVg45C5YYtduzZSu+4iLSpplJFxugf89+AlMn4c3NGfriYRCqqCL9iunGO+ibwAl+twF/Muj0IyLydefxq8CzzuMSYLWIWNgB8g9tV7RBRP4RWCciUaAO+GPnPY8DT4qIDwgBA8sAbwJ3A1VAH/CloXwdE8EyVYFSVakIbaAtdoIlGXeQ4zWjY6KqVDavob3/OFdPeISg38xIFaD22Iecaq9i8aLHCQTMbDwAtNbtou3Ybube/Kf4AuZ0bfqaj3OqYguT7/ojYzvVYFfqdO7YTPEd9yGGpskAGovRs20HuZ+82WjwUcuib3clOSvMdnjXSJTw4VoKvvSgQaM63J3uC5jSXqDwrHNPYqcCnX1tLXBOoS1V/QnOCPSs8xuwN4LOPq/A188+fyFGfbBMZaA8ENpEU7SWxRm3k+czl8x89NRHNHTtZ/G4B8lJN6cX3tp2kJq695kz635yss3JNYT7OqjZ8Srls24kt9ScJpGqRcPa18iaMIOc6VcaswvQ/tFaPMEgOYvMZAAM0Le/Aqu/n8xlZkerkdrjWF3dBBcZVnc8UodGY6TPM6slNdaqcxJh1K9ZVrHPeKAEOBTeyoloFYsybqXAZ05/+0TnPo6e2sgVpXdQlJl8fuYAff2tVB58hfHly4wkng+galG99WUCwRwmXmlWi6fj0E76mo4x7kazImvx/l46tm2k4LpPGlunHKB7y1YyrpiHL8fsVLl/TyX+8lKja5UAoYPV+CeU4c0zp8kOisbjCR1jiVE/sswgy3igBMjw5LAweDNFPrN2/d4MZhXdxPgcczK4AIKHkuL5zJx+j1G7qJJVNJmC8fON1n8DpBeWMe7G+0kvNPdlBCD+AMW3fYqcK82O/gAK7r83Je3Jcu5ZQea15v3NfWAFWbcsN2vUbdE2OhkvZhfZB5gUSI0KY0nWdMCs4D1AMFjAnFkPGLcrHi8T599h3C5AsHgCweIUqDv6/OQtvc64XYBAj4JKGQAABkJJREFUqbllk8F4An48ZeaWegYQjwdfvtlRMOC2aDONiNwpIoecYvXvnOP1NBF5yXl9i4hMSaU/Li4uyaOAWprQMZZIWbB0CuN/hJ1wOg87HeBske/HgHZVnQH8P+CfUuWPi4uLIdRp/pvIMYZI5cjyaqBKVatVNQK8yJl6zwHuA55zHq8EbhWjyWAuLi6pwN3gMcu5CtXPboF0+hpVjYlIJ3bOVevgi5x60oFk0vC7unK/cW/NNg8fbLeIs36epDHbv3Yw5n1NLaPJ39HkK5wnnxGgm/bV7+rKRLvSjKaf+YKMig0eVX0KeApARLYnUGs6YhhN/o4mX2F0+TuafAXb3/O9pqrmamhHEamchidSqH76GqccKRcwL2bj4uLikiSpDJbbgJkiMtVpwfQwdvH6YFZhdwUB+AzwvlOG5OLi4jKiSNk03FmD/DNgNeAFnlHVChH5LrBdVVdh95b7LxGpAk5hB9SheCpVPqeI0eTvaPIVRpe/o8lXGH3+phxxB3IuLi4uQzPqa8NdXFxcLgZusHRxcXFJgBEbLEdTqWQCvn5LRCodfeL3RGTypfBzkD8X9HfQdZ8WERWRS5bykoivIvI55/5WiMivLraPZ/ky1N/CJBH5QER2OX8Pd18KPx1fnhGRZhE5Z96y00V8WNraYxpVHXEH9obQUWAaEMBWdZt31jVfA37iPH4YeGkE+/pJIMN5/NVL5Wui/jrXZWPLkW4Glo5UX7F1n3cB+c7zkpF8b7E3Tr7qPJ4H1F5Cf28EFgP7z/P63cBb2DKzy4Etl8rXkXCM1JHlaCqVHNJXVf1AVfucp5uxc04vFYncW4D/hV2rH7qYzp1FIr4+DvxIVdvhtJ70pSIRf8+nY33RUdV12Fko52PY2tpjmZEaLBPR9P29UklgoFTyYjNc/eHHsL+tLxVD+utMtyaq6u8upmPnIJF7OwuYJSIbRWSziFzK6pJE/P0H4IsiUo+tA/ONi+Pax2LY2tpjmVFR7jhWEJEvAks5ozY34nAU9X7IGcGnkY4Peyp+M/aIfZ2ILFDVjkvq1fl5BPi5qv6ziFyLnWc8X4ehLetyaRipI8vRVCqZkP6wiKwA/ga4V1XDF8m3czGUv9nAfOBDEanFXqtadYk2eRK5t/XAKlWNqmoNto70zIvk39kk4u9jwMsAqvoRkI7dZGMkMmxt7bHMSA2Wo6lUckhfReQq4KfYgfJSrqnBEP6qaqeqFqnqFFWdgr3Geq+qnrexwqXy1eE32KNKRKQIe1pefTGdHEQi/h4DbgUQkbnYwbLlonqZOMPW1h7TXOodpvMd2Dtxh7F3F//GOfdd7A8u2H9kr2Dr/m4Fpo1gX98FmoDdzrFqJN/bs679kEu0G57gvRXsZYNKYB/w8Ei+t9g74Buxd8p3A7dfQl9fAE4CUewR+mPAnwJ/Ouje/sj5WfZdyr+DkXC45Y4uLi4uCTBSp+EuLi4uIwo3WLq4uLgkgBssXVxcXBLADZYuLi4uCeAGSxeXUcRQzS8+hr3vOw1IDjhNM1x11fPgBksXl9HFzwEjJZ0ich1wPXAldiHCMkZwddmlxg2WLi6jCD1H8wsRmS4ib4vIDhFZLyJzEjWHna8cANIAP3Y+sMs5cIOlS0oQke+KyF8Mev49EXniUvo0hnkK+IaqLgG+DfxHIm9Su9zyA+zE9JPAalU9kDIvRzluIw2XVPEM8CrwL05zjoexW5i5GEREsoDrgFcGLTemOa89iF09dDYnVPUOEZkBzOVMy8A1IvIJVV2fYrdHJW6wdEkJqlorIm1OXXwpsEtVXU1483iADlVddPYLqvoq9hfW+XgA2KyqPQAi8hZwLeAGy3PgTsNdUsl/Yrd6+xL2SNPFMKraBdSIyGfhtBTEwgTffgy4SUR8IuLH3txxp+Hnwa0Nd0kZTuedfdgbBzNVNX6JXRr1iMgL2F2WirA3Y/4n8D7wY6Ac+16/qKrnmn6fbcuLvb55I/Zmz9uq+q3UeD76cYOlS0oRkZ9gTxPPK4zm4jIacNcsXVKGs7GzHPjspfbFxSVZ3DVLl5QgIvOwe42+p6pHLrU/Li7J4k7DXVxcXBLAHVm6uLi4JIAbLF1cXFwSwA2WLi4uLgngBksXFxeXBHCDpYuLi0sC/H+tonYiWo4pGgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "E = system.hamiltonian.energy\n", "print(\"The system's energy is {} J.\".format(E))\n", "\n", "system.m.plot_plane(\"x\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Computing energies of individual term" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For instance, the exchange energy is:" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2018/3/26 16:32: Calling OOMMF (system/system.mif) ... [0.7s]\n" ] }, { "data": { "text/plain": [ "1.1217019053781098e-21" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "system.hamiltonian.exchange.energy" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can also chack the sum of all individual energy terms and check if it the same as the total energy." ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2018/3/26 16:32: Calling OOMMF (system/system.mif) ... [0.7s]\n", "2018/3/26 16:32: Calling OOMMF (system/system.mif) ... [0.7s]\n", "2018/3/26 16:32: Calling OOMMF (system/system.mif) ... [0.7s]\n", "2018/3/26 16:32: Calling OOMMF (system/system.mif) ... [0.7s]\n", "The sum of energy terms is 5.352855331445171e-20 J.\n", "2018/3/26 16:33: Calling OOMMF (system/system.mif) ... [0.7s]\n", "The system's energy is 5.352855331445172e-20 J.\n" ] } ], "source": [ "total_energy = 0\n", "for term in system.hamiltonian.terms:\n", " total_energy += term.energy\n", " \n", "print(\"The sum of energy terms is {} J.\".format(total_energy))\n", "print(\"The system's energy is {} J.\".format(system.hamiltonian.energy))" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 }