{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Isolated skyrmion in confined helimagnetic nanostructure" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Authors**: Marijan Beg, Marc-Antonio Bisotti, Weiwei Wang, Ryan Pepper, David Cortes-Ortuno\n", "\n", "**Date**: 26 June 2016 (Updated 24 Jan 2019)\n", "\n", "This notebook can be downloaded from the github repository, found [here](https://github.com/computationalmodelling/fidimag/blob/master/doc/ipynb/isolated_skyrmion.ipynb)." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Problem specification" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "A thin film disk sample with thickness $t=10 \\,\\text{nm}$ and diameter $d=100 \\,\\text{nm}$ is simulated. The material is FeGe with material parameters [1]:\n", "\n", "- exchange energy constant $A = 8.78 \\times 10^{-12} \\,\\text{J/m}$,\n", "- magnetisation saturation $M_\\text{s} = 3.84 \\times 10^{5} \\,\\text{A/m}$, and\n", "- Dzyaloshinskii-Moriya energy constant $D = 1.58 \\times 10^{-3} \\,\\text{J/m}^{2}$.\n", "\n", "It is expected that when the system is initialised in the uniform out-of-plane direction $\\mathbf{m}_\\text{init} = (0, 0, 1)$, it relaxes to the isolated Skyrmion (Sk) state (See Supplementary Information in Ref. 1). (Note that LLG dynamics is important, which means that artificially disable the precession term in LLG may lead to other states)." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Simulation using the LLG equation" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from fidimag.micro import Sim\n", "from fidimag.common import CuboidMesh\n", "from fidimag.micro import UniformExchange, Demag, DMI\n", "from fidimag.common import plot\n", "import time\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The cuboidal thin film mesh which contains the disk is created:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "d = 100 # diameter (nm)\n", "t = 10 # thickness (nm)\n", "\n", "# Mesh discretisation.\n", "dx = dy = 2.5 # nm\n", "dz = 2\n", "\n", "mesh = CuboidMesh(nx=int(d/dx), ny=int(d/dy), nz=int(t/dz), dx=dx, dy=dy, dz=dz, unit_length=1e-9)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Since the disk geometry is simulated, it is required to set the saturation magnetisation to zero in the regions of the mesh outside the disk. In order to do that, the following function is created:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "def Ms_function(Ms):\n", " def wrapped_function(pos):\n", " x, y, z = pos[0], pos[1], pos[2]\n", " \n", " r = ((x-d/2.)**2 + (y-d/2.)**2)**0.5 # distance from the centre\n", " \n", " if r <= d/2:\n", " # Mesh point is inside the disk.\n", " return Ms\n", " else:\n", " # Mesh point is outside the disk.\n", " return 0\n", " return wrapped_function" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "To reduce the relaxation time, we define a state using a python function." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "def init_m(pos):\n", " x,y,z = pos\n", " x0, y0 = d/2., d/2.\n", " r = ((x-x0)**2 + (y-y0)**2)**0.5\n", " \n", " if r<10:\n", " return (0,0, 1)\n", " elif r<30:\n", " return (0,0, -1)\n", " elif r<60:\n", " return (0, 0, 1)\n", " else:\n", " return (0, 0, -1)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Having the magnetisation saturation function, the simulation object can be created:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# FeGe material paremeters.\n", "Ms = 3.84e5 # saturation magnetisation (A/m)\n", "A = 8.78e-12 # exchange energy constant (J/m)\n", "D = 1.58e-3 # Dzyaloshinkii-Moriya energy constant (J/m**2)\n", "alpha = 1 # Gilbert damping\n", "gamma = 2.211e5 # gyromagnetic ration (m/As)\n", "\n", "# Create simulation object.\n", "sim = Sim(mesh)\n", "# sim = Sim(mesh, driver='steepest_descent')\n", "sim.Ms = Ms_function(Ms)\n", "sim.driver.alpha = alpha\n", "sim.driver.gamma = gamma\n", "\n", "# Add energies.\n", "sim.add(UniformExchange(A=A))\n", "sim.add(DMI(D=D))\n", "sim.add(Demag())\n", "\n", "# Since the magnetisation dynamics is not important in this stage,\n", "# the precession term in LLG equation can be set to artificially zero.\n", "# sim.driver.do_precession = False\n", "\n", "# Initialise the system.\n", "sim.set_m(init_m)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This is the initial configuration used before relaxation:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAocAAADlCAYAAAAht12aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAHihJREFUeJzt3X2wZHV95/H3Z54ASVZFLJzw4JCVkqUSETNLSOmmCGotmhTkwRg0ZVBIzaZWjXGzGyFuYtakanWTSiQbw2aCqMlagotEWHRjFGHdxApxECLPSjDCUCCQgHmaOA/3u3/0mUtzuXdu33tP9znd9/2q6qLP6dPnfLvnS893fud7fidVhSRJkgSwoesAJEmS1B8Wh5IkSZpncShJkqR5FoeSJEmaZ3EoSZKkeRaHkiRJmmdxKEmSpHkWh+tYkguTfDrJpUkeT/KVJKckeVuS+5M8luRHu45TkvR0/oZrXCwO17dTge3AVcDRwG3A/2le+5fArwL/uZvQpLXxL06tA/6GaywsDte3U4H/WlXXV9UB4E7gtqq6pKr2AbcDmzqNUFo9/+LUrPM3XGNhcbi+vQi4bmj5lEWW755oRFJ7/ItTs87fcI2FxeE6leT5wBbgK0OrXwzcOrT8IuDWJFuS/EWSw5P8SJJLJxmrtEr+xamZtcLf8GclubF5PJDkLZOMVdPH4nD9OpXBKMocQJJ/ATwf+PKCbf6yqvYCHwB+G3gD8NYJxyqtyAr/4nx1kl9p3ndUkk9OLFBp9VbyG/5EVZ0JXAx8EvjdCceqKWNxuH6dytP/ory3qv4JIMkG4LuGtvkCcCHw76tq/yQDlVZh5L84gV3Aac26i4H3TDBOabVW9Bue5HTgAuDNB/+/kJaSquo6BvVckqOBK4Hrgb1V9RsdhyQdUpJfAo6tqp9pll8GXFZVJzfLG4B/AE6qqgeT/CnwOuA3q+rHu4pbGock3wO8GdjhP+41CotDHVKSw4A/At7BoD/rC8CZVfWPnQYmtSjJx4GNwC9V1W1dxyO1KclDwH3APuDWqvq5jkNSz1kcSlr3kvwi8MKqOr/rWCSpaxPrOUxyeZJHktw+tO6oJJ9J8tXmv89u1ifJbye5N8mXk7xkUnFKWpdeBPxy10FI0rDFaqcFr4+lXprkBSkfAs5esO4i4PqqOolBP9tFzfpXASc1jx2AU6dIal2S5yW5Fvh/VfX1ruORpAU+xNNrp2FjqZcmNgFsVX0+ybYFq88Fzmyefxi4kUFv27nAH9TgnPefN3M0ba2qhyYTraT1oKoeBs7pOg5Js6G5+GdkVXXzMq8vVjsNG0u91PXdAY4Z+gAPA8c0z48FHhjabnez7mkfNskOBtUyRx555PecfPLJ44t2iszNdT9TwYYNzpR00M033/xYVT13Ne81xxdnjveH+d2+L935112HwEtO2dZ1CL0xco5v2LyLzc8Yeb9J/gG4Z2jVzqrauYLQRq6XVqLr4nBeVVWSFV8d03yJOwG2b99eu3btaj22abRnz56uQ+CII47oOoTeSLLqU5bm+OLM8f4wv9u35bQLug6BXbsu7zqE3hg1x3P4M9l88g+PvN+9X7rsnqravurAxqTr4vAbB4c/k2wFHmnWPwgcP7Tdcc06SZKk3sqGjZM83Fjqpa7PiVwLHJw64nzgmqH1P9VchXMG8E37DSVJUq9lAxs2bRn50YKx1EsTGzlM8lEGF58cnWQ38C4Gt6n6WJILga8Dr202/xTwauBe4J+AN00qzj7ow+myNrTxOTxtN5vM8SeZ47OnD6eE29DG59h7y/o6NR0gG9sbOVyidtoMUFX/gzHVS5O8Wvl1S7z08kW2LQa3+pEkSZoOCRvbGREEDlk7HXx9LPVS1z2HkiRJsyGZdM/hWFgcSpIktSCEDZs2dx3GmlkcdmBW+q0mYbnvyn6tfjLHR2eOT59Z6SmchOW+q5nrSXTkUJIkSfMSNmxur+ewKxaHkiRJLQgTn+dwLCwOJUmS2pANrV6t3BWLQ0mSpBYkaWty605ZHE6Azfnjs/C7tXm/G+b4+Jjj3fMClPFZ+N3OwgUqnlaWJEnSQHP7vGlncShJktSGtHv7vK5YHEqSJLUgjhwK7LXqm1H+POzZWhlzvF/M8XbZT9gvo/x59L0v0Z5DSZIkAYOrlTdu2tB1GGtmcShJktSSJF2HsGYWh5IkSW0IjhyuR/ZfTT/njTs0c3z6meNLs8dw+vV5bsTB7fMcOZQkSRJAwsaNjhxKkiSp4cihJEmSAEhg4yaLw5ln/9XsW+/9Web47FvPOW6P4ezrWw+iVytLkiQJcJ5DSZIkDQtekCJJkqQBp7KZUfZfadb7s8xxzXKO22OoTnsQPa0sSZKkYbNwQUovytskb09yR5Lbk3w0yeFJTkxyU5J7k1yZZEvXcUqSJC3l4FQ2oz76qvPiMMmxwM8C26vqu4CNwHnAe4HfqqoXAI8DF3YXpSRJ0vKyISM/+qovp5U3AUck2Qc8A3gIOAt4ffP6h4FfAS7tJDpJkjp0wa2fG/sxLn/xWWM/xqzLjNw+r/NPUFUPAr8B3M+gKPwmcDPwRFXtbzbbDRy72PuT7EiyK8muRx99dBIhSxNljmuWmd+aNbMwcth5cZjk2cC5wInAdwBHAmeP+v6q2llV26tq+3Of+9wxRSl1xxzXLDO/NUsS2LJpw8iPvurDaeVXAF+rqkcBklwNvBR4VpJNzejhccCDHcYoSZJ0SCFs7PGI4Kj6UBzeD5yR5BnAHuDlwC7gBuA1wBXA+cA1nUUoSZK0jIMjh9Ou8+Kwqm5KchXwJWA/cAuwE/gkcEWSX2vWfWAcx3dCYC1n2icMNse1nGnO8Wmc9HoSF5esxmrimoaLWCY5KXYCmxw5bEdVvQt414LV9wGndxCOJEnSiiVpfeQwydnAJQym+rusqt6z4PU3Ar/Ok+13v1NVl63lmL0oDiVJkqZdgC2bNra3v2Qj8H7glQxmbvlikmur6s4Fm15ZVW9p67gWh5IkSS1IaPuClNOBe6vqvsH+cwWDGV4WFoetWlfFob1XasNiedSXHi1zXG3oa45PY38h9LfHsA0LP9s09iBCe32IIWxZ2STYRyfZNbS8s6p2Di0fCzwwtLwb+N5F9vNjSb4f+Arw9qp6YJFtRrauikNJkqRxWcXI4WNVtX2Nh/3fwEer6ltJ/h2Du8qtqUq3OJQkSWpBAoe1e0HKg8DxQ8tPm/e5qv5maPEy4L+t9aAWh5IkSS0YwyTYXwROSnIig6LwPOD1TzlmsrWqHmoWzwHuWutBLQ4lSWrRLPcYLmcaexDb1PYk2FW1P8lbgE8zmMrm8qq6I8m7gV1VdS3ws0nOYTBX9N8Cb1zrcS0OJUmSWhBav1qZqvoU8KkF63556PnFwMVtHtPiUJIkqQXjmAS7CxaHkiRJLWl75LALM10cOuebJqWre9Oa45qULnJ8WuY1XM89hsuZlh7Etu6/vKHlnsOuzHRxKEmSNDHt3yGlExaHkiRJLQhh8wZHDiVJksRgKpvNGx05lCRpXbHHcPWmpQdxtQJsjMWhJEmSGBSHm+05lCRJEgAJGywOJUmSBI4cSpK0LthjOD6z1oMYYIM9h5IkSYKDVys7lY0kSZIaMzCTjcWhJElSG5I4cihJkqSBQc9h11GsncWhJElSCwZXK0//yGEvPkGSZyW5KsndSe5K8n1JjkrymSRfbf777K7jlCRJWlJg44bRH33Vl9AuAf64qk4GTgXuAi4Crq+qk4Drm2VJkqReCmHzhg0jP/qq89PKSZ4JfD/wRoCq2gvsTXIucGaz2YeBG4F3TD5CSZKk5Q2mspn+psPOi0PgROBR4INJTgVuBt4GHFNVDzXbPAwcs9ibk+wAdgAcf/zx7NmzZ/wRS8toMw/NcfXRcB6edtpp37Pa/QznN5uPZMtpF6w5tvWgjcmindx7aavNw1mZBLsPY5qbgJcAl1bVacA/suAUclUVUIu9uap2VtX2qtp+9NFHjz1YadLMcc2y4fzOpsO7Dkdas43JyI++WnFxmOTIJBtbjGE3sLuqbmqWr2JQLH4jydbmmFuBR1o8piRJUqvC4NTyqI++WrY4TLIhyeuTfDLJI8DdwENJ7kzy60lesJYAquph4IEkL2xWvRy4E7gWOL9Zdz5wzVqOI0mSNFYzcrXyKD2HNwCfBS4Gbq+qOYAkRwE/ALw3yR9V1f9cQxxvBT6SZAtwH/AmBoXrx5JcCHwdeO0a9i9J0rL62ofXRo/hcvvsw2dfLIZxfPZxCWEDPR4SHNEoxeErqmrfwpVV9bfAx4GPJ9m8liCq6lZg+yIvvXwt+5UkSZqkPo8IjmrZ4nCxwnA120iSJM266R83XMFUNkm2A+8Ent+8LwwuJH7RmGKTJEmaGgE2zsDNlVcyz+FHgP8E3AbMjSccSZLURZ9dH3sQp9EMnFVeUXH4aFVdO7ZIJEmSpliy/kYO35XkMgb3Of7WwZVVdXXrUUmSJE2hPs9fOKqVFIdvAk4GNvPkaeUCLA4lSdK6F+j1nU9GtZLi8F9X1QuX30ySJGkdCmyc/tpwRX2TX0hyytgikSRJmmIBNmT0x0j7TM5Ock+Se5NctMjrhyW5snn9piTb1vo5VjJyeAZwa5KvMeg5dCobSZKkRttT2STZCLwfeCWwG/hikmur6s6hzS4EHq+qFyQ5D3gv8BNrOe5KisOz13IgSZKkWdfyWeXTgXur6j6AJFcA5wLDxeG5wK80z68CfidJqqpWe9CVFIcPAz8GbFvwvnev9uCSJEkzI2l7KptjgQeGlncD37vUNlW1P8k3gecAj632oCspDq8BvgnczNBUNpIkSYJUkZUN2B2dZNfQ8s6q2tlyWCu2kuLwuKry1LIkSdISMrd/JZs/VlXbD/H6g8DxQ8vHNesW22Z3kk3AM4G/WUkQC630auXvXsvBJEmSZldBzY3+WN4XgZOSnJhkC3AesPBuddcC5zfPXwN8bi39hrCykcOXAW/0amVJkqRFVMHcgRZ3V/uTvAX4NLARuLyq7kjybmBXc1vjDwB/mORe4G8ZFJBrspLi8FVrPZgkSdIsy2gjgiOrqk8Bn1qw7peHnv8z8ONtHnPZ4vDg5dBV9fXltmkzMEmSpOlSsLKew14apefwhiRvTXLC8MokW5KcleTDPHmuW5IkaX0q2u457MQop5XPBi4APprkROAJ4HAG577/BHhfVd0yvhAlSZKmQZED0z9yuGxx2JzL/l3gd5NsBo4G9lTVE+MOTpKk9eiCWz/3lOXLX3zWxI+pVWj5gpSurOSCFKpqH/DQmGKRJEmabj0+XTyqFRWHkiRJWkqtdBLsXrI4lCRJast6GjlM8hngP1bVX44xHkmSOrNYb18fevHG0YPYh8+10CR6K8cpNRsjhyu5fd47gPcl+WCSrW0HkmRjkluSXNcsn5jkpiT3JrmyuW2MJElSf83Njf7oqZGLw6r6UlX9AHAd8MdJ3pXkiBZjeRtw19Dye4HfqqoXAI8DF7Z4LEmSpHYdvFp51EdPrWTkkCQB7gEuBd4KfDXJG9YaRJLjgB8ELhs6zlnAVc0mHwZ+eK3HkSRJGqfU3MiPvlpJz+GfAScCdwB/DrwRuBt4W5J/U1U71hDH+4BfAL69WX4O8ERVHTxxvxs4dom4dgA7AE444QSOOOLJwcw9e/asISRp9YbzcK3McfXRcB7ecsstN692Pwvz++u3XD7/2pbTLlhDhLOtj/2Cs2TvUB4CJB8c8Z2zcfu8lVytvAO4c5F7KL81yV2LvWEUSX4IeKSqbk5y5krfX1U7gZ0A27dv9/7OmjnmuGaZ+a2Zst4mwa6qOw7x8g+uIYaXAuckeTWD2/L9C+AS4FlJNjWjh8cBD67hGJIkSWNW1P59XQexZivqOVxKVd23hvdeXFXHVdU24Dzgc1X1k8ANwGuazc4HrllzoJIkSeNSzMQFKX2eBPsdwBVJfg24BfhAx/FIkiQtqaqofdM/ctir4rCqbgRubJ7fB5zeZTySJEmjW2c9h5IkSTqEmo2eQ4tDSZKkNlRR+/d2HcWaWRxKkiS1onp9W7xRWRxKknQIl7/4rKcsOwF1exZ+t1OvitrnyKEkSZJg/U2CLUmSpEPxghRJkiQddHAS7ClncShJ0grYg7h6M9djuEA5lY0kSZKeZM+hJEmSDvL2eZIkSXoKRw777YgjjnjK8p49ezqKRLNuYa51dVxzXOPSRY7vveXypyxvOe2CiccwCnsQlzYtPYYLc23Vqpiz51CSJEnA4LTyAe+QIkmSJAZXK8/t2991GGtmcShJktSGgrm9FoeSJGnIeu5BnJYew3GpKuYOeEGKJEmSYHBByoROKyc5CrgS2Ab8NfDaqnp8ke0OALc1i/dX1TnL7XtDe2FKkiStb3VgbuTHGl0EXF9VJwHXN8uL2VNVL24eyxaG4MihJElSK6qKA5O7IOVc4Mzm+YeBG4F3tLHjdVUcLjZPl/PCaaW6mtNwFOa42tDXHF9sLrq+zn04bJZ7EKexx7C1OQ0XU6x0RPDoJLuGlndW1c4R33tMVT3UPH8YOGaJ7Q5vjrEfeE9VfWK5Ha+r4lCSJGlcVjGVzWNVtX2pF5N8FnjeIi+9c8FxK0ktsZvnV9WDSb4T+FyS26rqrw4VlMWhJElSG1qeBLuqXrHUa0m+kWRrVT2UZCvwyBL7eLD5731JbgROAw5ZHHpBiiRJUhsKDuzbP/Jjja4Fzm+enw9cs3CDJM9Ocljz/GjgpcCdy+143Y8cem9aLaev/VejMse1nGnO8Wm5//Kw1fTpTaJPcRr7B0cx1h7DBYqJ3j7vPcDHklwIfB14LUCS7cDPVNVPA/8K+L0kcwwGBN9TVf0vDpMcD/wBg0bKYtCMecmo8/dIkiT1QhVz+/ZN6FD1N8DLF1m/C/jp5vkXgO9e6b77cFp5P/DzVXUKcAbw5iSnMPr8PZIkSd2ric5zODadjxw2l2E/1Dz/+yR3Accyxvl7JEmS2lZVHPDeyu1Kso3BVTQ3MeL8PUl2ADsATjjhhPEHKU2YOa5ZZn5rpkzw9nnj1JviMMm3AR8Hfq6q/i7J/GuHmr+nmSxyJ8D27duXmuNHmlrmuGaZ+T2aWb1YZOYU1IHpT+NeFIdJNjMoDD9SVVc3q0eav0eSJKkPBrfPO9B1GGvW+QUpGQwRfgC4q6p+c+ilZefvkSRJ6pO5AzXyo6/6MHL4UuANwG1Jbm3W/SJLzN8jSZLURzVXzO2d/pHDzovDqvpTIEu8/LT5e8bNCYM1zRMCj8Ic1yzn+DROiq12TXLS68X0eURwVJ0Xh5IkSbOgipnoObQ4lCRJakOVVytLkiRpoAoO2HM4++zPmn2z3H81CnN89q3nHLcHcfZ13WP4FM5zKEmSpHkzMs+hxaEkSVILCq9WliRJ0kHlPIfrkv1Z028991+Nwhyffub40uxBnH696jFcwKlsJEmS9KTq923xRmVxKEmS1IIqPK0sSZKkRnlBihitt8eercmx16p95ni/mOPtGqV/zb7EyelzP+EoqooD++e6DmPNLA4lSZJacqAcOZQkSRKDeQ73zlkcSpIkiUFxOAMthxaHk+C8ceNj/1U/mOPjY453z7kRx2faewwXqnLkUJIkSY2i7DmUJEnSwKDnsOso1s7iUJIkqQVVXq0sSZKkhlcra9WWazC3mf9JNuNPJ3N8dOb49FnuIgovWHnSrF1wshyLQ0mSJM0bnFbuOoq1sziUJElqQQH77DmUJEnSgFPZaEza6EHqQ0+XvVRaijmuWdZGn10f+hbXW79gG+aYjZ7DDV0HcChJzk5yT5J7k1zUdTySJElLanoOR32sRZIfT3JHkrkk2w+x3Yprqd4Wh0k2Au8HXgWcArwuySndRiVJkrS4g1crj/pYo9uBHwU+v9QGq62l+nxa+XTg3qq6DyDJFcC5wJ2dRiVJkrSIfczx0Nw/T+RYVXUXQJJDbbaqWqrPxeGxwANDy7uB7124UZIdwI5m8VtJbp9AbIdyNPBYxzFAP+Iwhie9cLVvNMeNYRl9iMP8bl/ncSQf7DwGevA9NEbK8cfZ/8RVPPzQSvabZNfQ8s6q2rmy0A5ppFpqoT4XhyNpvsSdAEl2VdWS590noQ8x9CUOY3hqHKt9rzluDH2Pw/yezTiM4alxjLJdVT275eN+FnjeIi+9s6quafNYw/pcHD4IHD+0fFyzTpIkaeZV1SvWuItV1VK9vSAF+CJwUpITk2wBzgOu7TgmSZKkabGqWqq3xWFV7QfeAnwauAv4WFXdsczb2jxPv1p9iAH6EYcxPKmtOPrweYxhoA8xQD/iML/b14c4jOFJfYljXpIfSbIb+D7gk0k+3az/jiSfglXXUqRmYCZvSZIktaO3I4eSJEmaPItDSZIkzZuZ4rCLW+0lOT7JDUnubG5h87Zm/VFJPpPkq81/W720fYlYNia5Jcl1zfKJSW5qvo8rm0bUcR7/WUmuSnJ3kruSfF9H38Pbmz+L25N8NMnh4/4uklye5JHh+dmW+uwZ+O0mli8necmIxzC/O8zv5pid53gX+d0c1xwffyz+hjO7v+HTZiaKw3R3q739wM9X1SnAGcCbm+NeBFxfVScB1zfL4/Y2Bs2mB70X+K2qegHwOHDhmI9/CfDHVXUycGoTy0S/hyTHAj8LbK+q7wI2Mrgya9zfxYeAsxesW+qzvwo4qXnsAC5dbufmN9B9fkPHOd5hfoM5vh5y3N/wp2otv6dSVU39g8GVOp8eWr4YuLiDOK4BXgncA2xt1m0F7hnzcY9jkLxnAdcBYTCj/KbFvp8xHP+ZwNdoLnAaWj/p7+HgTPBHMZjD8zrg307iuwC2Abcv99mB3wNet9h2h9i3+d1hfjfH6DzHu8zvZt/m+PiO6294ze5v+DQ+ZmLkkMVvD3PsJANIsg04DbgJOKaqDt4+52HgmDEf/n3ALwBzzfJzgCdqcAk7jP/7OBF4FPhgc1rksiRHMuHvoaoeBH4DuB94CPgmcDOT/S4OWuqzryZXze9u8xt6kOM9y28wx9vUdY53nt/QuxxvM7+nzqwUh51K8m3Ax4Gfq6q/G36tBv+0GNt8QUl+CHikqm4e1zFGsAl4CXBpVZ0G/CMLTj+M+3sAaHpCzmXwQ/cdwJE8/VTBxE3is4+T+Q30IMf7mt9gjq/x2H3I8c7zG/qb49Oe36sxK8VhZ7faS7KZwY/KR6rq6mb1N5JsbV7fCjwyxhBeCpyT5K+BKxiclrgEeFaSg7dHHPf3sRvYXVU3NctXMfihmeT3APAK4GtV9WhV7QOuZvD9TPK7OGipz76aXDW/u81v6EeO9ym/wRxvSx9yvA/5Df3K8Tbze+rMSnHYya32kgT4AHBXVf3m0EvXAuc3z89n0McyFlV1cVUdV1XbGHzuz1XVTwI3AK+ZUAwPAw8keWGz6uXAnUzwe2jcD5yR5BnNn83BOCb2XQxZ6rNfC/xUc8XbGcA3h05dLMX87jC/mzj6kON9ym8wx1vRhxzvSX5Dv3K8zfyePl03Pbb1AF4NfAX4K+CdEzrmyxgMNX8ZuLV5vJpBv8j1wFeBzwJHTSieM4HrmuffCfwFcC/wv4DDxnzsFwO7mu/iE8Czu/gegP8C3A3cDvwhcNi4vwvgowz6Y/Yx+Bf4hUt9dgaN5u9v8vQ2Blflmd89z+++5HgX+W2Or48c70N+d5Xjk8jvaXt4+zxJkiTNm5XTypIkSWqBxaEkSZLmWRxKkiRpnsWhJEmS5lkcSpIkaZ7FoSRJkuZZHEqSJGmexeE6l+SIJP83ycYW9rUlyeeHbnMkdc4c1ywzvzUOFoe6ALi6qg6sdUdVtZfBjPI/seaopPaY45pl5rdaZ3E4w5LckOSVzfNfS/LfF9nsJ2nuGZlkW5K7kvx+kjuS/Enzr9JtSe5O8qEkX0nykSSvSPJnSb6a5PSh/X2i2ac0dua4Zpn5ra54+7wZluT7gXcDvw+8Hjhn+F+XzQ3u76+q5zXL2xjcu3J7Vd2a5GMMbjL+p83604A7gC8Cf8ng/pPnAG+qqh9u9rEReLiqnjuJz6j1zRzXLDO/1RVHDmdYVX2ewU3C/wNw3iKnHY4Gnliw7mtVdWvz/GZg29D626pqjsGPy/U1+JfFbUPb0Bxjb5Jvb/OzSIsxxzXLzG91xeJwhiX5bmArsLeq/n6RTfYAhy9Y962h5weATYusnxtanhva5qDDgH9eTczSSpjjmmXmt7picTijkmwFPgKcC/xDkrMXblNVjwMbkyz8cVnLcZ8DPFZV+9rap7QYc1yzzPxWlywOZ1CSZwBXAz9fVXcBvwq8a4nN/wR4WYuH/wHgky3uT3oac1yzzPxW17wgZZ1L8hLg7VX1hpb2dzVwUVV9pY39SWtljmuWmd8aB0cO17mq+hJwQ1sTqAKf8EdFfWKOa5aZ3xoHRw4lSZI0z5FDSZIkzbM4lCRJ0jyLQ0mSJM2zOJQkSdI8i0NJkiTNsziUJEnSvP8PDxeB9wPCGzcAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plot(sim, component='all', z=0.0, cmap='RdBu')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now the system is relaxed to find a metastable state of the system:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "scrolled": true }, "outputs": [], "source": [ "# Relax the system to its equilibrium.\n", "start = time.time()\n", "sim.driver.relax(dt=1e-13, stopping_dmdt=0.1, max_steps=10000,\n", " save_m_steps=None, save_vtk_steps=None, printing=False)\n", "end = time.time()" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Timing: 77.00890803337097\n" ] } ], "source": [ "#NBVAL_IGNORE_OUTPUT\n", "print('Timing: ', end - start)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "sim.save_vtk()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The magnetisation components of obtained equilibrium configuration can be plotted in the following way:" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We plot the magnetisation at the bottom of the sample:" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAocAAADlCAYAAAAht12aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztvXmUJFd17vvtHGuee6gepG4NaEAghmYymAsILoMx2Bj7gv0wBu4Sxgbj4dnAlS/YXFgPnr1s42fgIQM25rIAPywbGWRjLMATWEgCgZBaQ6sltXqeap5yOu+Pym7V/k50ZlZXVmVW9fdbq5Z6R0RGnIjYcfIo8jvfthAChBBCCCGEAIBUqxsghBBCCCHaBw0OhRBCCCHEWTQ4FEIIIYQQZ9HgUAghhBBCnEWDQyGEEEIIcRYNDoUQQgghxFk0OBRCCCGEEGfR4PACxszeYmZfM7OPm9mYmT1gZleb2TvN7ICZnTSz17S6nUIIIWLUh4vVQoPDC5trAewB8CUAIwDuBvAP1XWXAvhfAH63NU0TYmXoi1NcAKgPF6uCBocXNtcC+L9CCLeGEMoA7gVwdwjhIyGEIoAfAci0tIVCnD/64hQbHfXhYlXQ4PDC5skAvrIkvjohvm9NWyRE89AXp9joqA8Xq4IGhxcoZnYxgByAB5YsfgqAu5bETwZwl5nlzOy7ZtZhZj9tZh9fy7YKcZ7oi1NsWJbZhw+Y2beqf4+Z2dvXsq1i/aHB4YXLtVh8i1IBADPrA3AxgB/SNj8IIRQAfArAnwJ4A4B3rHFbhVgWy/zifIWZ/V71c0Nm9tU1a6gQ589y+vDxEMILALwHwFcBfGyN2yrWGRocXrhci/iLcl8IYRYAzCwF4Jol23wbwFsA/EoIobSWDRXiPGj4ixPAHQCeWl32HgAfWsN2CnG+LKsPN7NnAngzgF8981wIcS4shNDqNog2x8xGAHwRwK0ACiGEP2xxk4SoiZn9TwDbQwi/XI2fB+CTIYQrq3EKwDSAy0MIh8zs3wG8HsAfhRB+tlXtFmI1MLOnA/hVANfrf+5FI2hwKGpiZnkAfwvgXVjUZ30bwAtCCDMtbZgQTcTM/gZAGsD/DCHc3er2CNFMzOwIgP0AigDuCiH8eoubJNocDQ6FEBc8ZvY/AFwRQnhjq9sihBCtZs00h2b2aTM7bmY/WrJsyMy+bmYPVv87WF1uZvanZrbPzH5oZk9bq3YKIS5Ingzgva1uhBBCLCVp7ETrV2W8tJYTUv4SwMto2bsB3BpCuByLerZ3V5e/HMDl1b/rAcg6RQjRdMxsq5ndDODfQgiPtro9QghB/CXisdNSVmW8tGYGsCGEfzWzXbT41QBeUP33ZwB8C4vatlcD+Kuw+Jv3f1Y9mkZDCEfWprVCiAuBEMJRAK9qdTuEEBuD6uSfhgkh3FlnfdLYaSmrMl5qdXWALUtO4CiALdV/bwfw2JLtDlaXRSdrZtdjcbSM7u7up1955ZWr19p1xLHvJ76BPifGMS9I2CZFG/Fr6IFrn7isNmxk7rzzzpMhhE3n81nleDLTC37SJedsheTUZVpQSdBbl2lRPU325p587UZeICi/m8+JmULN9Y302QynM2d3qeIdbkZ7O+rv9AKh4RxPZe9Atqvh/ZrZNID7lyy6MYRw4zKa1vB4aTm0enB4lhBCMLNlz46pXsQbAWDPnj3hjjvuaHrb1iN/1PuEZW2fpo4ll4p7mg5a1pn2w0H+zGt0L85iZuf9k6VyPJn/ePiUi9P07Thf9l90E/NFF08XytE+pwt+wDlf8vvgY7z9x3Y31tgNjvK7+Xzitkdqrs+mfP+b0GVH8P8wcX6PzfoB6Q3XLe97ZCPTaI5bRz+yV/5Uw/stfO+T94cQ9px3w1aJVg8Oj515/WlmowCOV5cfArBzyXY7qsuEEEIIIdoWS6XX8nCrMl5qdYWUmwGcsY54I4AvL1n+i9VZOM8GMCG9oRBCCCHaGkshlck1/NcEVmW8tGZvDs3s81icfDJiZgcBvA+LZar+2szeAuBRAD9X3fwWAK8AsA/ALIA3rVU724E/6LncxQX6LWCOhVEJ25RJXMI/h/HPyLw+6WflzjT/rOyPwT87f3mr1xx2D3dG++wY9JqWXLd/WJ75tW9EnxHrn7m5ORcbVfMqJfx/60zRbzPSlXXxLK3nn5WZJM1hkZ6t6Xn/M/M4/ez2+s/c7uLT07FOrEjt6Mj6twq3vO3HarZTrD8+eOsDLh7s8v3aUKfPXQAYpnzes63fxTmS8VAIi1SIMYFUhvSrMuZLXmrxrYdOunhszkszAOD4zIKLD5z2z/YHX35V3XZtJAyApZv35vAcY6csAIQQ/l+s0nhpLWcrv/4cq65L2DZgsdSPEEIIIcT6wAzp5rwRBFBz7HRm/aqMl1qtORRCCCGE2BiYrbXmcFXQ4FAIIYQQogkYDKlMLBtYb2hw2ALYZmaahB8cz5FmKUFyCFZXsWIrbbU1iDn6AGsSAaBQ4W38PlmnyNrIzGRt3y4ACPSZ7/3US1yc7/O+ck/8q6/U3adYewrHH3GxVbx2L53rcXHI+vtq6dhfjWWwGVrAcVIO14O1ugv8bJIGsUDrWV8IAEXyY5wn/7pn/N4/uTiT9Q/jd254cY0Wi1bAmsLRPp+vz9w54OIR0hz2cocLoIvueyfFWc53zu9QW2MLADC/zxJ9lxQrnP8+ninGvp4zRe/pt2vAa8s/+73HXPzgiRkXv/+lG8zXUm8OhRBCCCHEWcyQyjZPc9gqNDgUQgghhGgChjX3OVwVNDgUQgghhGgGlmrqbOVWocGhEEIIIUQTMLNmmVu3FA0O14APdXtT64miNxrliRtsaN2AzLgu9Uyu68WNbJONjLXrt6tC514hQX9h2puupslA+IG3vdbFT/j4l+ofVDSd8n5vBp1O+a4lpH1nGWg96GeYdELPxOnEJsAcpyLj9+XPUCmzuXydOHDxWgBlmrRSmPMTVErUHyzM+X08+V1fdfEPP/wTNVosVoNP337AxTzhZLTXT9Toz/t87qMJKHnEZtKp+XEX29ysj0vebNrKNMHvPCakZOi5zGf8eXTn/GSTgc7eaJcLZOg92OHPfTtN1tlMhQ4+9p1HXPwrz9kVHWO9oZ+VhRBCCCHEItXyeesdDQ6FEEIIIZqBNbd8XqvQ4FAIIYQQogmY3hwKAPiLIW/gOcmVzJFkal1bY8jEhtbxNvU0hZ3peutTNdcDQEcdzSEfoyND+paOON3SpMVJpWvH9Tj8gbdFyzqG+1w89LYPL2ufFzoL//p5v6CUoJfq9RqsVI+PraPbx6zJYU0imWYDQJo+w6bA/Axw/mUpP7OpOLd4WZqPkfBcLCVJc8g62jLrahdYg+jXV6j/2PWWL0THyHZ4rdiDH/3pmu0Uj/Per93n4is290TbXLvVa+02dXmd3QDp7PLFaRenJsd8POdjAAjTEy6uTHkNYnnem0eHhXm/g4rXriZCz5DlvR4wRc8pP9fpnv5ol+kOv6yje9jF/X3+evblfBvYIPxLdx+OjnHP0SkXv+8lV0TbtBPSHAohhBBCCACLs5XTmeW91GhHNDgUQgghhGgSdh7OCO2GBodCCCGEEM3AoDeHFyKfG7nKxawxnCjGmsNy8Dok1gyydo/X19MTArHejzWEPZSsrB+sp0kEgAx5DKaoMHyatCSsJ0zWHKZrbpPppPWdfr2RRqxciPVw5TnvBzb9v3/fxT3/x/uiz1zIzP39n7k4kJbJksTWrN9j30KqNWplf58CebZZguYwS+aHGcpRjnP0DHRR/s4n6IM7iv4zvZSP0/N+H52Uv5MJXwr8FoF1iRyXCv56l0hrVi6Rvx2AcmHOxaxLfORTr4s+c6HyoW8+6OKn7/CauUuHvLcfAGzp9nnQZ/4epCa8Ti41dcLF5VNHXFwY8+sBoDx23MUL4163WJzxGsPyvPc9LBfjZ4ZJZ/15pEmrmu32GsT8gNcLpgc3x/sc3OTizPCoi1O9fn2e4k7qGzqzSVpg/wzxPXz3C72XcCtZLJ+nN4dCCCGEEAIAzJBe5kTKdkSDQyGEEEKIJqE3h0IIIYQQAgBgBqQzGhxueG7a8kQXzyzToxCINYOs76unQaynHwSAPlrWTbqNfJ/XluS6vU8Xa/mS9IHsOWh1iienGiiunMqyTpG87DpZI+P1KRmKk3wRjZYF8pljjV3nT769Ros3HqzB5OvD1y/V6b3QACCQ5jDkvHapMud1c5b2+ccaw5Cgq8tk/D5Zh5SjfCsHH/fkMrQ+OgSKXDuZ9MKzpAeco3h6PtZ9FcnHMFVHrB4q5ItI16JcIH87AKUy6UJJp7jzDX/h4sc++6aabdhIsD7tmi3es/DyYa8x3Nwd933dhUkXpye9hrBy4jEXF476WswLx/z2s0dPR8eYPe69D+dPeW+/hUmvMSxSje6QlNAE99ncv/L3RMewv1Zdm49G++zaOuT3scVvk9l6kY83+dzs7/MaxUyP96UFgLR1U9zeGkTNVhZCCCGEEADkcyiEEEIIIZZi0IQUIYQQQgixiKxsNii37LrWxXNcFzXyLGR/wHif9XwMWVPYTRsMkQ6vc8BrrwCgZ4vXzXSSjqZrxGs2cr1+fYb8rdI5rwlLIpDOiX222F+QtVRAXHeWYV+uVI61kaRB7O6M9sGfiTSI5OPHtYTzz399zTauN05//F0u5vsUaQzp+iWWlM+QhnB+1sWBdIqRj2G5ROvjOrFcbzmT8u3qiETgqRoRkEnFOc46xg66Fvy8s7diOUGDXCA/RdYgFuZIy5tQ83kplWKsx6yU+FmjZ3PB+yBue/0nXHz482+tecz1xO9//X4XP2Wb17CxxnAraQw752M9YGb8oIuLB/dR/JCLJx/xurupA8dcPH3EaxgBYOqw9zWcG/Pa0rkZ7w3K302N6N/r6dk7SYveOei/F3q3nYz22TPqr2/vRf769U2ccnF2xp97dofXUnb3b4uOsbV7hJb4e8jfy1wv+/0vvTLa56qhn5WFEEIIIcRSNsKElLYY3prZb5jZPWb2IzP7vJl1mNluM7vNzPaZ2RfNLPGFhRBCCCFEO3DGyqbRv3al5YNDM9sO4NcA7AkhXAMgDeB1AD4M4I9DCJcBGAPwlta1UgghhBCiPpayhv/alXb5WTkDoNPMilgUExwB8CIAP19d/xkAvwfg46vdkHnyjmLrqHoaw6SaxLyM6xr3kydh97DXzfXt6KsZA0Dfbu8V1bvT18DMjGx1cbp/2MXG3nVJuqeS17ywd10gLUl5epxir6kB4nqhpflYT1WLyOcwG6c01wtln76oVnAprs+8kZg/RfeJ6lGzVyTrOJM0cbm8v4+RTyFdUyONEMchxFpU1ilmM/5eB/I1TGdr+yB2ZmONVg890IPkA9dP/p89pAfONaA12kdasVKBtGMLXi9YpDrVSbDmMKq/TLWXNzJP3Lo8H8POOa+jS5/2noUAUHh0r4vnHvYaw/EH/GfGHz5Bse8Lxw7HfeFJ8sw8zR6apdoaQ9bdJRF/f5EX6LTPmyHSPY4c9N6LADC4zS8bOOm/FxbGaD31+Z3k25m9ONYbdw7752pzt/dWLFb8PZ4uxPtYK2yDlM9r+RmEEA4B+EMAB7A4KJwAcCeA8RDCmW+DgwC2J33ezK43szvM7I4TJ+Ji5kKsd5TjYiOj/BYbjY3w5rDlg0MzGwTwagC7AWwD0A3gZY1+PoRwYwhhTwhhz6ZNm1aplUK0DuW42Mgov8VGwmzxF4RG/9qVdvhZ+cUAHg4hnAAAM7sJwHMBDJhZpvr2cAeAQy1soxBCCCFETQyGdBu/EWyUdhgcHgDwbDPrAjAH4DoAdwD4JoDXAvgCgDcC+HLLWiiEEEIIUYczbw7XOy0fHIYQbjOzLwH4HoASgO8DuBHAVwF8wcw+UF32qdU4/ree/mwX50lgXs9YlCeXsME1APRQorCJNRtYD14y6OKBy3e6uP+qy6Jj5C55oovDiC92Xu73E1bm4M1Oz8dANW9e9JuaYbNTH1dOx0Xbc2PH/TZTXrgd5mnSCxlvW5rE+pnY2NhyNCGFYp6Mk+r2E37K+293cfqSZ0THaGceeNtra67nCSY8yYdJXM+TeMq1BeGBhPGBJ12k4q4pGJlF03qecBI9ignG2hF0jBJ1kUM0IWWky18LnrACAJ3Z2hNKHuCJBmSSXZzvodg/EwAAr+mPjLLZFHue1g+97P0uPv2P7z1ne9uNj33nERc/a0e/izd1+XvSXfATsjITh13Mk08AYPahB1186kcP+/h+34+d3jfm4oMT3uj56DyZwCOegMJ9cu1yAY1Su5/nNnA7uSADAGylz+ygSSwLk/7cedLhMBVH8N+Gi+Son+8e8vd0S7efhHQlFX349O0HXPzmZ/jvx2ZiFpvjr0daPjgEgBDC+wC8jxbvB/DMFjRHCCGEEGLZmFnT3xya2csAfASLVn+fDCF8iNb/EoA/wOPyuz8LIXxyJcdsi8GhEEIIIcR6xwDkMvXtpxren1kawEcBvASLzi23m9nNIYR7adMvhhDe3qzjanAohBBCCNEEzNDsCSnPBLAvhLB/cf/2BSw6vPDgsKlcUIPD219xXd1tUmxIXSa9IOlAUpQEmSTNUVS83GuIBi/zhtUjT77Uf/4ar3ELFz85OsZYh7eAODjpNWAPPex1Sqdmve5jhnQjuQQTT9ZTbe3Ju3hHnzcm3TLiz6t7YEe0z8wk6RDHfYH6MhVtZw1ipG1jDSJik2vr8KoWy3uT51SX16+w/q1w/JHoGLnNu6JlreCeX3xltKw0F+ublpIho+dAGqBUzq9PMipnLSjfB9Z9hqy/ByHn40rO3xMACOnaWkg2ybbCrItThRlaHxtDW9mfW4Y0iPmsb9dgr3/uBjv8sw0Ag51eB5unn5xK1KdEmkOKC7P+WgKA1THKZpNs1iByPPJK96sVAODkV95d8xhrwXu/dl+07Lm7vEZ7U5e/3n3mzz09ecTFxYP7XMwG10CsMTx5r++nju/1RtoPzfj+99Ccj2e4ukIC5OGOftaupzmuPyCZo+POkSE7G23PU24eTtBKThRrm3VfSp8Jdc49yWSf+49Mxn/39A37mHNgZ7///AdvfSA6xg3XPaFmuxrFYInfnzUYMbM7lsQ3hhBuXBJvB7DUZf0ggGcl7OdnzOz5AB4A8BshhNjNfRlcUINDIYQQQojV4jzeHJ4MIexZ4WH/HsDnQwgLZvZWLFaVe9FKdqjBoRBCCCFEEzCLfx1YIYcALLUsiXyfQwhLf2L7JID/e6UH1eBQCCGEEKIJrIIJ9u0ALjez3VgcFL4OwM+7Y5qNhhDO6CVeBSD2Y1omF9TgMEl7ZaTTyJL+Kp3g6+TXk+5jMNYD9Yx6DdvAZVtdPHztFS7uePLzXFy86GkufnAs1nx9e6+vSfrtfV6r9/ARX/x8jgqsl0l7ksnG/+fT3efP7aJN3kvq2osGXPykLf68LxmMHay2bfY6j0yH/0ymx3uWhekJH1PR9kTI+5A1iCnyOQzpLMWkWazU1vC1EvYUA4DCtNc7pUgPUybNEK+vFOqfL/tNWoe/pkb3sdzp40r3sItnS7EuqVykHKUOuJO0Slby1yI1673nknw3S+S7GeZIp0i5lBokvfC22IO0d8Qv66SZjJXgz3WaNFoF6rfmZ/y1A4DijPfus7lpF7OmkDWI0f7mp2uubxWXU58DANupXxro8Nc3RT6GlRNeilU86DWG4w/EUi32MWSN4QPUnx7ie0baPdYTAsCmvP/u2UnfRVv6vK6ua9jrX3M9sccrw33B7Cmvuz1G/cdjdB4nFuK+gPWTB2b9McqU36Brx9/BST6q2W7yqiUv2lTO58VAv/cGHu311260L/6ebhbNNsEOIZTM7O0AvoZFK5tPhxDuMbP3A7gjhHAzgF8zs1dh0Sv6NIBfWulxL6jBoRBCCCHEamFo+mxlhBBuAXALLXvvkn+/B8B7mnlMDQ6FEEIIIZrAaphgtwINDoUQQgghmkSz3xy2gg09OLz/+te4OEmTEXm6sd6KtHisQcx1e31E95ZYV9e3izWGV7o4/5Qfd/HcqPcx/P4Rr3v66l6vfwGA/yDfrWMHvDZv6rjX3bBGicmSpgMAuof8eUyd9nqV06S7mSLt1EIprg4agr9eo72+BnSWPAZT5I+XKnovO2N9C+I6vuxbWCFNYaQ5zJI+JRPrVebm/LXo7Ix9+laD777UuxUkaQ45h+Mc9+eb6fBxttiA5pA0hul+ryGs9Hg/wGK3j0/PeU3cbDHOFb6N3aSLzeVJaxb8PsKsz/ny8YPRMeYOPOLbcdRrd9njMdfrz7tv9z3RPvPXPMfFl1/iq4KWKl6rO0nPzekZf8zZ6fgeL0xTffBZ//wv1wexNBdrDodf8UEXn7rlhpr7bAY3/IPX1T//kqFom3667/mib3tqyuuxC0d9nd3JR7z2dPxhvz0Q10qOfQxrawy7SVd3WU+sq7t0q/fI3PIk8tC8YpuLe7b79fkB8mdNYGHca8+nD/lz3X6//57Yfrdf/9DROC/2Ub/PGkS+Nml6kLN0bfOkrQSA/KA/t8F+fw9zfb6/yXd578t+8rJl30Mg9j48X9/DVJM1h61iQw8OhRBCCCHWjOZXSGkJGhwKIYQQQjQBgyGbUOVlvaHBoRBCCCFEEzADsg2UMmx3NvTgkPVBqWysu2F9Q7ngdUqxB5O/ZJ2DXsvQvX0kOsbQVbtcnLvy6S4ujl7j4h8d9/o11hj+y11eFwIAh/d5zeHEIa+fWJjw2hHWGKXI+y83H2t7mEx2u4tPkS/XftLVbOqJtSRcr5md5TeRVg1U69YKfp9Gfm5AguYwTbWESYOIOhpE/jwAWIg1cmtBgfRoyZpDrwFib07O8XIhvoZLSfIhS/V63RwGtriw1Oe1pEemvGbrCOmWZotxG/LkD7i5m2oWs2cp31fSF5enx6NjsMZw7IFDfv1J/2yynrN326PRPjcf8fsceKE/1yuvfIE/5pz3MTx42utqp8fjmtBzU/55LUx5HVexjuaQtdflhLrT5TreiKvBRUO+fx3tjbVifZTPKarHXj7laykvHPPx1AHfd44/HOfFwQn/XHGt5Hoaw6vIZ+/qJ1O/BmDHj1/u4pFnXOvi3GVeix4G/DMV8nFdbya34DWDfeP+Wozs+6GLB2//gYt7/u3BaJ/ZH/rvlr1T/lrFGkR/7XomyCs44fp3jfh71LXV53t62OvhU1TzvK/Pa3JHuuI+bDBh2flgiHWV65ENPTgUQgghhFgrDEBWmkMhhBBCCAEAMENKg0MhhBBCCAHozeG6IJ31p5dLkGSEMnne1dFncY3Hzk1ea9W/2+tAACB7sfc1rGzx2pID015fdcdh71H23Ye8hubkYe9VBQAzVC+0MHXaH7OOXog1h5lc7NOX7eihbfy1SdfxdipWYl0ea8vY326WvOy6qPZyinRlIanuMekUA+uv6NwjDSLF5dhKEQmL1oT5MV9bemGyvi6Mcz7TQXV369Ra5mcAANJUY7jc5zVAh6f9fbnvpNfRHZjw55GUK/157q6852UPaYp7sz6HuS4yaxABoDBJ+r5j3mN0bL/XQ82R310HebYBsU5xF60f6vP6qStGnujia3f4Pua+hOd/gq7nXKd/VlNTVB+cngHWHHItZgAoL8Q6xNVmc7fX6rG3JQDkQbXD533/WRjzmrjZo75vnD7i/S/HDsdefkfJe5J1dFwrmX0MWWN4ySueEh1j6L94z1Jc6rXp43nv5Xd63t+jmfn6uufunPevHdrpaxD3bb7UxVu3XOTiXO83EvZ6l4uKP/Aa+Xsna2sQ+dqOJFz/PrpHvXQP86P+Hue2sHei10D35uI8GuqsX5u6EQxASppDIYQQQggBnJmtLCsbIYQQQghRZQM42WhwKIQQQgjRDMxMbw6FEEIIIcQii5rDVrdi5WzowWF+wIuyy8V4skIo1xbxsuFvrs+L4Lu2epFwbvvF0T7S27zId7rDC9AfJQHug8eoODqZ3hbnY6NjFphnu72RboYmk6Sy/rzyPb5NnYNewAsAvcN+H/3D/lqMjnij0R1Dfn08oSA2NmYdbzl48XKgySMV8i1NMqOuN8GkRLNJimRoW6ZJMkmTT+gjqG9H2xxmTvncmEuYLcM/cXRQzmc6axskG5WCSvfEZ5ca8hNQpjJ+4tD+Ez7H7zzkJw3sPx6L0BnOJ/6/85EuLyjf1OW3R77OBJUESiSWnyaz7sO0vjzrJ0cAwMwd3miYjfS7dn7LxZuu833IFfRcXb41vv6njvnrNzPmJx7M07mX5vl6+3bzBBWg/qS21WCQJgkkTUjhCShh2sflMT9BYva4n6wwRf3vyQQT+NN1jOE3Ud92Kd0jNriOJp8ACFc938WPLvhz3/uYn5TxCH0vHE8wwGc2U9GHXQM+L66iCZYXU5uSSiMUpvykreljfnLUiQV/7fiZ4WubdP2H6B710z3so3vMOZDq99+pXVky7Qcw3NW8CSkboXxeW5yBmQ2Y2ZfM7D4z22tmzzGzITP7upk9WP3vYKvbKYQQQghxTgxIpxr/a1fapWkfAfCPIYQrAVwLYC+AdwO4NYRwOYBbq7EQQgghRFtiMGRTqYb/2pWW/6xsZv0Ang/glwAghFAAUDCzVwN4QXWzzwD4FoB3rX0LhRBCCCHqs2hls/5Fhy0fHALYDeAEgL8ws2sB3AngnQC2hBDOiHWOAohFcADM7HoA1wPAzpF+zHz+A2fXdY16PWBpJjZxrZD+ig1/0x1eo8Gaw/xmb/6b3rQ9Oka5xxugjpEG4xTplMYp5lI8ua7YhDgMb3NxR/+Ii1mTmOnwOqauXi/e6xmIj9FD+pRLNvl9XLLZ62x20vbb++J9bu6urSnK1FP2pkmDmLAJSQZRIpEhawxLldrrk2Sq/Jl9p0/FG50nS3N8JJXFTVseN0lmvU7S+bNBL5Mv+E9FRvCd/hqne2K9TrnbP2vHZryu6G7S0d6x35vYnjrtdUpJhupzdK69pN3j/NrS7dcP5MlAvcvHAJDK1e4SC3SfJyi5JoqxXmqOEqbnrmMuHrpir4t3PPkhF28bfrKL+TkDgB/Q88vPN5vaF1m724CekI2xR166aO2vAAAgAElEQVT5obP/Tvdvfzpv3yhL87t/8zZ88NYHzq57yeW+H8tn4mS2WZ87lSlvVL4w7vVq86d8Ls6RkXySvpDvOz9TOzt93mx5ku/zR55xrf/ApfHlYo3htw/487jtYf/M3H+EzmOqvuaws9d/n10x6p+Bsd2kKrzIP+u7Eto98owDLt5y/2EX7zzuNYknFnzfwNc26frzPeJ7yPc4RzlgBZ8jnX2xepJNsD9x2yPRNo2wUUyw2+GdZgbA0wB8PITwVAAzoJ+QQwgB5yhCEUK4MYSwJ4SwZ6S3O2kTIdY1S3O8n6u7CLHOWZrf3QOSlov1T9qs4b92ZdmDQzPrNrNmfkMdBHAwhHBbNf4SFgeLx8xstHrMUQDHz/F5IYQQQoiWY1j8abnRv3al7uDQzFJm9vNm9lUzOw7gPgBHzOxeM/sDM7tsJQ0IIRwF8JiZXVFddB2AewHcDOCN1WVvBPDllRxHCCGEEGJV2SCzlRvRHH4TwD8DeA+AH4WwaCRnZkMAXgjgw2b2tyGE/72CdrwDwOfMLAdgP4A3YXHg+tdm9hYAjwL4ufq7CQjlx/UKuV6vD2TPQiD2OTS6W9lur2PK9Hn/wPSw93dLDXgNIgAUO/1nFqZZ0+bb0JXzL2Y7e7xOZIBN9QCUEjSCS8nl/T5zHV5fMdTvP79jyGuUAOBi8lvbRhov1nyx71xfLn7h3EXinU7SmqVZRUe6J/YwLCT4/LGmZZ40h/wZ1g+y12IjmsPV+rmgAn8+83RcPlcAyJFuM21+m546Hm7prL/G1t0XbVMhzeHBQ17f+wPSTx055D3bpse9piiTi3vNQOc2QPnFGted5OnW3+n1U6ne+CfMHElTct2kt6T7yrmR5DPJPm+HSD+1bd9RF2957AEXD+94ml9P/QEA9NKyPGngMp1ep5ie8deKNYesLwRi3fJqkEmlMNj1+DXvIB/UbIIG2Upea1ee9xq34oy/3gvkBzg34zXe06UEj0eK+6mf2kK5NniF14DnLvO60fG8f16A2MeQNYZ3PXDSxaeOep3dwrQ/7yTyPT6/6+kU2WdyaGfc7gE6t8ErfuDiLbd7n88e8godI91u0vXne8T3kO9xhXIgQzmSlEe5Jo3UDIYU2viVYIM0Mjh8cQghcnYNIZwG8DcA/sbMVuQeGUK4C8CehFXXrWS/QgghhBBrSTu/EWyUuoPDpIHh+WwjhBBCCLHRWf/vDZdhZWNmewDcAODi6ucMixOJn1zzg0IIIYQQFwAGIL0Biisvx+fwcwB+G8DdSLZTaz2WguUf173l0l6vslSPeM5dUM1hy3kdXarX65bS/V6DUcnFdjpl85e5Qi9aO7O+nVxDlhkfjPWATI40MazP2kz6wM3kfzWaoGsapW0GyGeuL++P2UOehSnymgIAK5D3JNXdhPl9hAy1q8Pr3xIkX5HGcI40LQuRBtGvn6ftKyE+SNJxV4NKCM43jz30kjSH5eA7qg7quCpcO5pOxug5SnXGOT4b/DZHp73G58AJrwEap3jmtNfdZTtiLz+jdu8nPeDlW71n2/g8aWRHvPY3O+i96ACgY9jnU77PH4N9OFnPmQTfk0nKp9mT/hngWsD5itdTsQ4MAPrp+c5Srd8092PUz6Ua8D1cC82hAehY0nex3WWCzSGsTHrJBX+9ylSLvjjn+5hGniGmk3437Br2fXLPdp9bYWDUxafn4+8irpXMPoasMZw49KiL5ydP1GjxIh19nPO+jvf91Mdfuc0/D1dtip/9fjo3Pne+Np3H/LPPmsOk68/3iO8h32POAc6RpDzin4JXUr1kA/yqvKzB4YkQws2r1hIhhBBCiHWM2YX35vB9ZvZJLNY5PjtMDyHc1PRWCSGEEEKsQ9rZv7BRljM4fBOAKwFk8fjPygGABodCCCGEuOAxrJ6V2VqynMHhM0IIV9TfrIWk0kh1LdFIJHh11cMyXrtjpK9ibzTr8TqmSjah7jHFJFvCSJfX+5SG/TFHevz6JJ89LvTdSzVi+0gfyB6ErB/sz8f6oj7ynkvPTbjYqJ5lasFrZrAQ17YOJZroTjoP6/Kal0qXv/7lnNdnFiux2mOBNF4zVEt4uuD1K7NUH7dYru17uJYEeG9D9tVrpG3sjVii8+V64xGZWPMW1Rwm7ejMpNcAscZw9pSvx8q1gIFYF9dN3nKPnvRapiObvG5xB+kHh3t93V4AyI94DXHXZp9/XSO+XUNUB51rLQP1NWxlvv5UK9iKXk/VlY27ba4znSFPUb6eHJcTnk2GNYdL42Z9FZoBS3+RM95zSMhNXkb9frnoc5E1tXx/GnmGOqm/zfX4ZyI/4PWvIe9zcWY+Po/j7L9IHoTsY8gaw/kx/0w1Qke/fwbmpnw7uU3cdwJAIH0wnztfG752TNL153vE95DvcfTdzzmSkEeca+f9y7ABdU5xXbAc3eS3zezqVWuJEEIIIcQ6xrA4sGz0r6F9mr3MzO43s31m9u6E9Xkz+2J1/W1mtmul57GcN4fPBnCXmT2MRc2hrGyEEEIIIao028rGzNIAPgrgJQAOArjdzG4OIdy7ZLO3ABgLIVxmZq8D8GEA/20lx13O4PBlKzmQEEIIIcRGp8m/Kj8TwL4Qwn4AMLMvAHg1gKWDw1cD+L3qv78E4M/MzEI4fw3UcgaHRwH8DIBd9Ln3n+/BhRBCCCE2DGbNtrLZDuCxJfFBAM861zYhhJKZTQAYBnAS58lyBodfBjAB4E4ssbJpJyyVRqq7t/6GS2GRdb0JKTRBIuTIrDftRe9J5MnZdXO3P2YXzVipVLx4nCefJO2TzXp7aDJJLxtUz475eOJUdAxM+TyrjHlBdHnK76M46yekRJNPgGgCSrqHTMY373CxZeIJP0spJRmoksn1xIIXL0/M+3bN0aSCYmX5nu+pVZqtFoIXZ7NQO9kE28ed6dqfKdFkknKB7lsD14ONxEt8TWcm/fZTp/32mfg5ypDwfXrQP5sHT/tJFUfIiHtiwU9eGuz2k08AIL1pu4u7R/02vaP+GdhyzE8emS7F14YnDfXRs5qhySSpXO1S9Ulidza9z/Dznaltep3O+z4myfA6MkNfug/bCLa/QqwcCwG2vBd2I2Z2x5L4xhDCjU1u1rJZzuBwRwhBPy0LIYQQQpwDq5Tqb/Q4J0MIe2qsPwRg55J4R3VZ0jYHzSwDoB9AwhuexlnubOUnreRgQgghhBAbl7BoldPoX31uB3C5me02sxyA1wHganU3A3hj9d+vBfCNlegNgeW9OXwegF/SbGUhhBBCiARCOC+P5XPvLpTM7O0AvgYgDeDTIYR7zOz9AO6oljX+FIDPmtk+AKexOIBcEcsZHL58pQdbddLpyKTakVBIOzJ2pQL1yHmj3ZD1uqVAZrIhHeuFWCPUSVW/0+bbwAbU7LaeT6gazsaiKTKoTk16PaBNeb1g6YR/S104wW+tgYWT/i31/CnWjXmT1kqh/qv1XJ+/nl1bvcarM+u1UdnhbX4HpHVirRsATJLG8OSsL8I+SRq7+dLyNYd8j5J0oc0gwGsI2TC2nuEyEBexZ01ciYral2a8gXWY8/cZAHIkwO4mE+Y0V7UnKqVCzRgAFqa9LnF2yutRT0/4dh4nE+Fxus+VTXFfkSXNYc92bxLcu80/F4PH/bXYeZCM3xHrEEfJ1L5nC+maSXcbsr4PKpbrP1cpuh/pDN2PvO/n0oX6JtjpSKf4+D6sSRrbEIClKRy4hECStpGXUZ+eJtNwYwNrvlaJ5+Lbwc9MYdrrchfGfR7kFqZd3J3z2nUA2Eym7p29Ps73+Dzp6NuU0M7a8Gd4n3xMblN3LuE7lM6Nz52vDV87Jun68z3ie8j3mHMgypGEPAqg/mEF792ssTeCDRNCuAXALbTsvUv+PQ/gZ5t5zLqDwzPToUMIj9bbppkNE0IIIYRYXwRgeZrDtqQRzeE3zewdZnbR0oVmljOzF5nZZ/D4b91CCCGEEBcmAc3WHLaERn5WfhmANwP4vJntBjAOoAOLv33/E4A/CSF8f/WaKIQQQgixHgiwBqQf7U7dwWH1t+yPAfiYmWUBjACYCyGMr3bjlo2lnGYw0g9mEk435ZcFitm3kPU/rEFEOj4GKyg6yJOskw8J/38TrOlITce6JvYprIwdd3Hx6AEfHz/s4ulD3r9t5kg8C372uNcYLkx6XVhhJtaJLYX93IBYb8V0bPOarkBakRK9/J4rxQ/l2FyxZjxO8RT7/JH4JMnglH0mexLOdTVg+U5y0XqO/Tbz7JU4469HYcp7+ZUTPDDz8J8Z7vLPTSfp7LLdXnPFvntJmsNK0S8rzfvcKMz7fR6f9BrEk7N0n4v+WQaAwV6vycpu9hrXnp1HXNw/5rV6lQTN6+CcF6f3bvM533uRPyZ7LRbSvh+bLXo9MRDnKGsAU5navofpXH2fQ/5MZulnEvTc50OA1/jy5SwliJcy1Ecb6yk7/H3OUofbSXpY1rclwbrd2VM+D6YPeU1337jPm6GdO8HsGvD34IpR79k7N8X2whe7qKN/BPVgjeHwVu8dysfkNg11JPhfnvDnxufO14avHZN0/fke8T3ke8w5EChHkvKILUpZe94wTZ6Q0iqW9Q0WQigCOFJ3QyGEEEKIC5E2/rm4Udbm9YYQQgghxIYnLNcEuy3R4FAIIYQQollcSG8OzezrAP7PEMIPVrE9K8T/1h866tc9DqStCVzTlTWH5GNYb/tqqxxsU2gFr+myBa8pTM94jVfl9NHoGAXSFC4cOejiqQPHXDz9mNckTh32usapI7Gucfak147MUL1c1ruxdqSnJ742KdKS9F7kdWVc6zpkvQZmmgR1p2bj+s3HSQt5mPzwJsj3cK7g9SKNaA5ZY9jfWbs+brPga55kIcbbsM8Y+/AtTHpt0+xJr+0rnYzzr2PCq00uHdzq4t2kZTpKWr65MZ+fXGsZiGv7BtL1FBZ8PE65wH6WnDsA0N/lvQ8zW51JA3ov8s9VkfSYlqC9qxRJc7jD+xgOXem1Y+mtu1w8Qec1XYj1TAt19FFp8oXLsD8rXctKKX5W+dyWahAtUlY3B/YtLSaYz+UzXm+W6vC6umy315/l2U+Qatv3TMd619N0zfmZOUbPzPb7vaZ7ZN8PXdy3+dLoGFdt8nkxtnso2mYp95Mn4dxUzzm2fBz2MWSN4bPomFdt8teybyHWGxfo3Mbo3Pna8LXjJ6YnEz9DfI/4HvI95hwIlCNJeTRf8vd4bLa2hv5cWNgYbw6XoyJ+F4A/MbO/MLPRZjfEzNJm9n0z+0o13m1mt5nZPjP7YrVsjBBCCCFE+1KpNP7XpjQ8OAwhfC+E8EIAXwHwj2b2PjPrrPe5ZfBOAHuXxB8G8MchhMsAjAF4SxOPJYQQQgjRXM7MVm70r01Zlv+ALfoj3A/g4wDeAeBBM3vDShthZjsA/ASATy45zosAfKm6yWcA/NRKjyOEEEIIsZpYqDT8164sR3P4HwB2A7gHwH8C+CUA9wF4p5n9eAjh+hW0408A/A6AMwKIYQDjIYQzP9wfBLA96YNmdj2A6wHgoosuQuba/3p2XeHUYdo4oZ4i13Fk30PSFLIPYuQzlyD6Ys8xLk9pZa+NSs17vR9rDIsH90XHmH3Eaw4nH/EasPGHvffUxKPes3CMNIfHFmLNBOtu2C+Pa2L2Z/2J7kyomcm+cLleqjM7uNlvTz50YzO+TUciLzDgCPndHR332knWprHmkEnSHPI+/uAnr3HxO2rusTac47/56ANn1/1O9hK3bSVSuMbwfWPfMdaSzp3yurop0qsCQOdJn387dnsftx+7zNfM3n+YPDOnn+DbcOKx6Bjsvce+e0yBtE1TVOubzxMAKgP9Ls5S/uVHfT3nPqo7ne6or37pvWiL3+cVT3FxacAfY2Le5+PEQqyr5Zy1qLayfxbT7AFLsAZx8TP+h6KH//zxUq726dffWXOHNeD8fuuzdp1d9/UHfK6N9sRa3u6c95pN9XrtXn7Aa/E6hr3OrnPQX4uhMX9PAeAo11+nZ+gxqke+/W7f3w7e7qX6W7d4LSsAXHzV8/2Ci/x5DJKO+cpt7OsZ930M10pmH0PWGF6cp1zbG9/mk3Rux+jc+doUqYvqoFwdysVeinyP+B7yPeYcqFCOLCQYHc5Sf3DDdb5P+t3oE+diY5TPW85s5esB3JtQQ/kdZrY36QONYGavBHA8hHCnmb1guZ8PIdwI4EYA2LNnj+o7iw2HclxsZJTfYkNxoZlghxDuqbH6J1bQhucCeJWZvQKLZfn6AHwEwICZZapvD3cAOLSCYwghhBBCrDIBoRS/3V9vNKXmUQhh/wo++54Qwo4Qwi4ArwPwjRDCLwD4JoDXVjd7I4Avr7ihQgghhBCrRcCGmJDSzibY7wLwBTP7AIDvA/jUcndQ6fC6hGTNoV8WaQhJW1KMvP1ofwk/iqRTvNAfM891TEmkGhaofutUXNZ65qj3n5p41Mdj+/1nThzx3nWPUX3hEwtx0rI/FdNJXmrsV9XBRaQB9G7zWpG+3d4lKbXtMhefLHtN14EJr5V8lPSEAPAo+fQdGfeaohnSV1boPFm/xTEAZNPNqS27XNhLMqkVfNfq+R5O0vlPH/Oaw8mH4wqaA/v9Dwv9o16v82M7vX/gw1f7fZZIM3eyi2qWAyjOe01VrsvrkDLZ2vegSOeZ5A3I1yJdx/ewi2pAdwx7HRgAWJZqEo/ucnFqx5UuHk/5cz8+43Oa/RqBuJ9iMlnfx+Ty/lkspeobT+QSnt/Vhn1LWRcGAJUu38+ne7xuNE260a7NXsPdu83XlR85GHu8sg7uMN2DE9SHPHTU90s9//agi3O934iPQfGuS5/u1+/0ul3WB84k+HYy3SR451rJkY8haQxP/0vc7oN0bnzufG0YvrYjCZpD/p7o2uyfS77HRjnAY4GZufhajc01521fCAGhuP7fHLbV4DCE8C0A36r+ez+AZ7ayPUIIIYQQjXOBaQ6FEEIIIUQNwsbQHGpwKIQQQgjRDEJAKJ1f6b12QoNDIYQQQoimENq6LF6jbOjBYdlIdJ0g2mZxPhfk5s+QX3D0+SRdOM9V4H1k82TgSeJZ6/TCYxa4J1Epes1Daa524XgW4rNRchI8GYKFxTtJwL7p6pFoH5ufcrmL89c8x8ULI35CyoGTXpz/IBk0P3g0FpM/etJvM0sF1XkyBFt5pugGsqEwACDfmkeJJwGl47kyUU7yJCq+11Fu0ASeyYPewBoAJvZ6Y/aRHT908eXXvMTFP3GVF5Az3+09FS2bpslGKcq/zh5v8DvQ5U2D+TnkCSpAbI7blffPYmZoq4vZkpknjwEJz++In9RSGt7l4lOT/iepk5SvUw1MSOEcTVGcoYkJPFcvadJVR/fal7c/TZMEphImXSx0+XalO3mygjfO79rqp370jPpJRIPb4j5kK/URE9S/zlAu7Zv29yz7Q28MDdwVHaMw5SfOjTzDG8sPXPZkF/cP+Ml7ocN/jyRhC36yiJ3wk8sK+/xzywbXPPkEAO6lc+NzZ9Prbuqktnb4vnNwW3wefI/4HvI9rnT4CSkL9KTOJEwYOT5T30S8IUJAKOrNoRBCCCGEAC48E2whhBBCCFELTUgRQgghhBBnOGOCvc7Z0IND1l/NJujoCqy1K69Mg5ik1GP1TiHDOkUfD/VucXFmxOuYMgkm2L2n/bKFca8tWaCi7KMJuqWlsIF1En20zXYqjr71Kf48djzv6mgfAy/8ry4uXeKtLfePe+3GD0hT+MPH/HnvO+bPGwBmJ71mrkAG3/VMr9MkVksyA+7tiY1b1wLWfXYmmHGX65iXx6bYfvvTpLfqTtAcjj/4mN9m9HYf93kD36fufJqLsymv5dsxGJsyP3jM3/txMkjuIs3rjiFvJt2f97qjVEKK87UIpAcud/vzSKW93s1CfK0DaYhL9HyfID3wSTovNr2eK9R+doE4Z9kgPFQytN6fd5LmMJdf+xwfI70l6y8BYKTLn0sHG5cPe21efguZYF902sUDZJoPADvGfB/CutwDdM9Yg7h3yve/xR8cj47BZvNb7j/s4sErvP6vZ7vX2eUHqOhDAgvj/hmaPuT1gmN0zGN3+/VscA3EGkM+9w7Kpe2d/jnc0U9a4d0D0TF6L/LPTH6Lv6dpusclek4nSat6ZMrfTwA4cDrWC58PQVY2QgghhBDicaQ5FEIIIYQQZ1D5PCGEEEII4dCbw/amo9Prlo7PxnqJAumr2PON5VqzpL8qktllkncakyet3gJ55JWD1/ZsGrnExRk2JQPQk/N6v2yXj/MDXjvVt8NrSYaPeZ1NKUGTmM76dnWN+Os7cJnXjY087SoXdzztBdE+C9u8d9d9p7w257uHvKbwtv1eH7T/CGloxmMtyTzplOpp8FhvlaP7k6THuuVtP1Zzn6vFb09737H3d14WbcM5zXE930PWEM2cirU5Y/vHXJwfeMjF6Z5/d3Fnxl/Tp4xe4+KhLq8ZAoArRrxf4CnSefGz2En5OkJ+eD25WEPHd7YM8gMk38NKxj9nIRXvs0K6xbF534eMUzxG3n6zxfpfNLnIx9C3g3OYPSIrdM95ewD493e9qG47ms0N1z3BxZ++/UC0zWiv16z195FvbK/X5mW2ep/JvgnvqbkwFvscsmb7UuofWat6aM6v52fo3snYU+8EaaF3Hvd98pbbvSdh17Dvf3M97LoZU5im3KJn+Ri16zE6jxML8fcC+xjGGkOfS5d2+3YOXeY1ogO7/f0CgL5d/ruF72GF7vFC1ufAxKT/Djg+E2tXP/jyq6Jl50UIqEhzKIQQQgghACz+rMwzVdchGhwKIYQQQjSBEAIqxfqOAu2OBodCCCGEEM0gAJUG7KbanQtqcJhgAQeeVDRX9K+Dp0hjuFDy8WyRNYfx62S2V2TJ2hhpoybJB2q26OPhPq+3AIDeK33d4q6tu1ycv+IRF28+cci3cdpr+yqFWDOR6vK+celBXx83u9Prg8qbLnXxUYtrZt5/0OtA7zw84eLbSWN48GhtjeH8bNzuqHZyhX0Na+u12AOOPePaiSR/SvYtZNgPlDWI/PmJYry/TvJoG9vnfdwy3fe7eIS0ebl5//ldW3zNbQDYss3XUx1b4Gezdp1zvm2s/QWAHF0MVhCHXFfN9QmlfzG74BdO0/WbJB0XawxZz5YEaw5Zb1mhPoXrhTP8DLQLRyZjTfGRKa/77KO251lzuMlr+bIz3rdzYCY+Rmnea9QCPyR7T7owbT6PDs3V9kEEgMPztfV9PeQn2Ek6ca6znsRcubanKfs3sp4wCa6VzD6GrDHcfJX/rhq+wn+PDDxhZ3SM7A7/XZLa5Lcp0T1mHe8hypsHT8Rels0ihIBKWRNShBBCCCEEsDghZY1+VjazIQBfBLALwCMAfi6EMJawXRnA3dXwQAjhVfX23b6vQIQQQggh1hmhXGn4b4W8G8CtIYTLAdxajZOYCyE8pfpXd2AI6M2hEEIIIURTCCGgvHYTUl4N4AXVf38GwLcAvKsZO76gBofbB2PN29gRr3GbIN3HFAlLWZM4TxrEhYT/E2DvQ9YQsT6FPZiOdPjbxH5tALC52y/r7/earf7NV7o4W/YajHQx9t1iQtb7iS2kvNbnCHliHZ7w53H/Sa8fBIAfHPRax3sOef3POHl9zZB2ZJ6uVWk+9rKs1DEkTWf8tUul4jxZyndueHHN9a3kN6ceiJZ9qNvnQtp8jtbTtLE8aj6hRvkcaT+nDvv7kO329WxZ5zlU8J/PXuy95wCgZ5vXHXX2eJ1RpaffxWXzz00D8qlIf8kqLt4Hyy+Tnv95WsbayAXSefHl5f6hMxd32105/+wVqM9Ik9C5nOBjuJS/f+tzaq5vFex7CMTeh0OkeevMUt/Y5+vwZnf4vq+zEGsOhyt1vFEpcbL7/C97PRM+348m+MhyDXP2Gx2jZOO4GfBPiexZOJSgRd1Kuca1ktnHkDWGw9fsdnHnbv+cA0B2h/dvLdE9nAz+Hp8gb9uHqG7y+1/qvw+bSsBy3wiOmNkdS+IbQwg3NvjZLSGEMwaYRwFsOcd2HdVjlAB8KITwd/V2fEENDoUQQgghVovzsLI5GULYc66VZvbPALYmrLqBjhvM7Fz/D3xxCOGQmV0C4BtmdncI4aFzbAtAg0MhhBBCiObQZBPsEMI5f6oys2NmNhpCOGJmowCOJ20XQjhU/e9+M/sWgKcCqDk41IQUIYQQQohmEIBysdTw3wq5GcAbq/9+I4Av8wZmNmhm+eq/RwA8F8C99XZ8wb85vGbU65T+9ke+fiV7jtWLpxO0JKwpmlumQSZrjHo74tvWR8v68153w3Vku8gHLZvm+s5xO2aLXscxNuf1gMemvXZn/3GvO3vwaKwHnBjz/naTpA2Zm/L7LMx6jWiZ9EGVYlwzk0llY83mUkoFv48ffvg1dffZzrx7xtdf/kif120l+RbWgrVQQOyPlhvz92WaamCnc16DyPTPxxrYHHlxpjdtd3F2wGuZ0jlfizlkvUYW6TgPQprq09JzAa6lTHrNpBcGvKye9pF9ULOkZ0vyZ+xJ6BOWwn0G+yL+2c9cW6dV7cubn+F9Xz/3/YMu5uuV6e1zcffADhdnL441yl0UR96oHT6X8n1ed9f5sM/dkcNxX3iSNIesQeRnLK6RXl9Vy/rVHCUb+6SyxnAkQXM4uM1rtAd2D1DstcHsY8gaw+zFcY3jUv82F8/k/D08QTWjHx733yO/8pxd0T5Xi4A1LZ/3IQB/bWZvAfAogJ8DADPbA+CXQwj/HcBVAD5hZhUsvhD8UAih/QeHZrYTwF9hUUgZsCjG/Eij/j1CCCGEEG1BCKhwdY1VO1Q4BeC6hOV3APjv1X9/G8CTlrvvdvhZuQTgt0IIVwN4NoBfNbOr0bh/jxBCCCFE6wlr6nO4arT8zWF1GvaR6r+nzGwvgO1YRf8eIYQQQohmE0JAWbWVm4uZ7d5K3FwAABJOSURBVMLiLJrb0KB/j5ldD+B6ALjoorjmsBDrHeW42Mgov8WGYg3L560mbTM4NLMeAH8D4NdDCJO2RDhby7+nahZ5IwDs2bOnEZ/bmmRJoFuk1748AWV81msLpudjrQFPUplls1MWGpdqv2pmQ9ukZZ1cfJ6Exkn7qNeGKTqPCTr3KZqQMjflJ3awgTUQm1gvTHuj7BJNeikXvNA41DG4BoAUmVxbyl+bFMXZDi8mbzXNzvG+hAkNS2GhO4vYk+AJTCXKlTmaoJLp8BM/0tmTdY/RQ/83np+hSS7DNGGlf9jFKTLJDrnY7DzkOn2c9VMRgvFzRLmV0E3xo8ZXkyecZNP+GHmKkyaflOie5TI+p8u0nvuHVtLs/L6fJsLxZLxsyk9U2to95OLOobjvy6X9PiznJyZlu32cH+x1cdfIMRf3HfFm/wAwRJNU+JmZm/H97Vy59gSVJHgCSiflVmc3GYgP+vPq3RY/Mz2jfnJI70X+PU7fLm/Rl91BE1DY4JomCAHAXIe/R8dnfF/w4Ck/sfGeo75veO2yFXcrIAAhaUbnOqMtBodmlsXiwPBzIYSbqosb8u8RQgghhGgHFsvn1X9p0e60fEKKLb4i/BSAvSGEP1qyqq5/jxBCCCFEO1Eph4b/2pV2eHP4XABvAHC3md1VXfY/cA7/HiGEEEKIdiRUAiqF9f/msOWDwxDCvyOW4Zwh8u9ZbV55tddHfPZ7j7mYtTuszUsywY51iaTHolfQJUqsMmlLynU0icBigtZcT4apFdq+lGCMzO0qLvjzWJjzcWne6wULs7HOpjTndTYl0hSyqXU9jWGSwTVrDI30Q2nSmT340Z+ueYz1zptO3+fiz41409mZcm1z3UY0iJxP5QLpo0hrOjfm73sq5/WDAJAifVSo+H12UG4EMkhPl/xzmOpJEI2H2s8W51I67/MtnXBp+HrRaSBLhspdWR8vlPwxi0nPNrk0cz/FGuN3PveSeB8bhPe/9EoXf+ib3gQ+zl9/8TZ3j0T77CbD9EzG65JT3V53N9h/wB9hq9fM9R712moA6D/ubXznT3nd3MKk13QXqb9tROdmrG/t9EMANu/uGCbt5ObBaJ98bvktoy7ObPWTjFKbvAl2qc9vzwbXQH2N4V7Smb7vJVdE+1hL2vmNYKO0fHAohBBCCLERCAEbQnOowaEQQgghRDMIQbOVhRBCCCHEIiEAZWkONz5veJrXR3zkP/a7mLU9HAPAHCXKDGn1CnNeC1Wg9az/qyRoDkv0Gpt1iWX6PxneR7lUpvVeEwbEesAyabrq6QfLCz5OOg5rCjmOPAozscaQ4W2yHd7n7JFPva7uPjYyv3Byr4tv2vJEF7OfWtL/FLN/WqqOjyZ3nuyLWJyJPTELk15nlCY/ykzHhI9Za5rxHm5IxWYNqUyGF9A+vO+bVXy7s6nYPzCd8hcsQ9emi3z45ul6d2b5+sfHYH9W1hhyP3Yh8e4XXu5i1iCyprYcSMAJYFOX1971DZPmMOf7lFyf99hMD3ste370RHSMvjHv1rYw7vtTfibK816DWG7AeDmd9fnMz1Dk1zjgfQ3Tg5vjfQ5u8vGw1xBWev36EsWTwffPx6Zir+CHTtfWGPI9binyORRCCCGEEGfZID6HGhwKIYQQQjSBAM1WFkIIIYQQZwjyObwgYX+w3//6/S7OJdStZR0i+wWyXrC4UNtPsFSINYesU4z0gOwfSFq/8sJ8zfWJ+yAPQq57XE9PeK5ltahXJznX5evnAkCWPMge++yblnXMC43XHLvHxbfsutbF83OxtilPurkUefWlyF8t8ixkb8WEHC/NU76R5oo1WekuHweq0x2oPi4AgPSoyNBzUPaxUY6ns7FejTWGrA8skSaxl64le0Ym2Rxm874rZ79W8TisT+M+PMlHcveA90Ld1OX1qwP9XtOZ7/J+gKker0HMbfH6WAAI035Zbsp7fVbINzZQn41G+lLW4eb9M5Ci/E/1Dvjte+L+tdLpz7VE576Q9brF8XnfzhOz/hm676Q/TyCul81elu2ErGyEEEIIIcTjhPYui9coGhwKIYQQQjSBEKCflYUQQgghRJWgCSkCjdVw/MlPfKfmetYUcV1k9ixkjSIQawyL8+Q5SHor1geWIw1ifU/CejHXvm0EI+851hhyHeRsp9eznLj5t5d9TFGbVzzyg7rbfOvpz665nmu6pnOpmjFrFJOokB9g4JjrclNtZY6BuB6zkQaLNYYh6/WXSc1mjWGONipXfGxZrwvLk475RZd7nzixMhrpwz/2nUdcfOmQ74e293ntXn/er+/r9/WF8/1x7qX6fS1lK3hvv0zJa2yN9K/16oIvfoi0vlQzOlDN6ErOa2grHd7vEQAW4PWXk6QXnpj07Twy5c/jsQn/zL35Gf5arTdCCNF39npEg0MhhBBCiCbBxurrEQ0OhRBCCCGaQABQSLIUWGdocCiEEEII0QQCkkuMrjc0OFwD/v6tz3Hxi/7031ycK/vbwD6HVqdOLRD7BVZITxX5GrJesFhbP5i0bLkehawnXFzm9VWsKcyQpjDb4eOTX3n3stogVocX3PmfLr79Fde5uETeiLker1NKkc6ONYqLy1iXmKq5vi4N+G5y7WTj9fTzkSXovtJWT3Po291JJaC3D/qcF2vPrzxnl4vf+7X7XHz5Jq9N3dzttXyjvV6T2J2Nc7U76/0B8x1DLmbtaoYfkfPQHJZoEMMej3NFv8/Z2fgYUwX/vXCSfAtZU3hk0sc3XPeEc7d3HRKC3hwKIYQQQogqAUGaQyGEEEIIscii5rDVrVg5GhwKIYQQQjSBEDRbWQghhBBCVNFsZXHefOPXfrzm+md/4OsubmRCSrnkJ3IYmWBHZtIkrGchPm+/2A4/aYA/E4n5afukfabJLDZN25y65YboM6L9ecYtt9Zcf//1r3Fxad6L2NPZuGviZZlunzspWm9pn38NwcbtJT8hBZk6xsMJkwIyKd+uPE1I2dTnz0O0P+9/6ZU119/wD3tdfBGZZm/u9mbTADBIM5G6aNJKR8bnM3mjw9DAxEX4fCbfeBRowek5P7Hx1Gxs3s3bjNGElI024aQeGhwKIYQQQoizLP6s3OpWrBwNDoUQQgghmkAAUJTmUAghhBBCLCIrG7FK/OfvvmTF+9j1li+4uNLhdSClwpyLMwvemDTJBHu5ptcyqBbn4oobb1rxPmY+/wEXhzLlZ4bcpFPL1yCGgn8uMD/jwtzF1y57n2Lj88GXX7Xyfdz6gIsHu7weu4NEhw1I08FSuCJrbIm3PmtX/Z0KRwUbQ3O4zJICa4uZvczM7jezfWamkYYQQggh2peq5rDRv5VgZj9rZveYWcXM9tTYbtljqbYdHJpZGsBHAbwcwNUAXm9mV7e2VUIIIYQQyZyZrdzo3wr5EYDXAPjXc21wvmOpdv5Z+ZkA9oUQ9gOAmX0BwKsB3NvSVgkhhBBCJFBEBUcq8/U3bAIhhL0AYFZTU3BeY6l2HhxuB/DYkvgggGfxRmZ2PYDrq+GCmf1oDdpWixEAJ1vcBqAN2mH2npa3AW1wHapccb4fVI6rDXVoh3Yov5tPy9vxy23QBrRHG4AGc3wMpfEv4eiR5ezXzO5YEt8YQrhxeU2rSUNjKaadB4cNUb2INwKAmd0RQjjn7+5rQTu0oV3aoTb4dpzvZ5XjakO7t0P5vTHboTb4djSyXQhhsMnH/WcAWxNW3RBC+HIzj7WUdh4cHgKwc0m8o7pMCCGEEGLDE0J48Qp3cV5jqbadkALgdgCXm9luM8sBeB2Am1vcJiGEEEKI9cJ5jaXadnAYQigBeDuArwHYC+CvQwj31PlYM3+nP1/aoQ1Ae7RDbXicZrWjHc5HbVikHdoAtEc7lN/Npx3aoTY8Tru04yxm9tNmdhDAcwB81cy+Vl2+zcxuAc57LAULG8DJWwghhBBCNIe2fXMohBBCCCHWHg0OhRBCCCHEWTbM4LAVpfbMbKeZfdPM7q2WsHlndfmQmX3dzB6s/repU9vP0Za0mX3fzL5SjXeb2W3V6/HFqhB1NY8/YGZfMrP7zGyvmT2nRdfhN6r34kdm9nkz61jta2Fmnzaz40v92c517rbIn1bb8kMze1qDx1B+tzC/q8dseY63Ir+rx1WOr35b1Idj4/bh640NMTi01pXaKwH4rRDC1QCeDeBXq8d9N4BbQwiXA7i1Gq8278Si2PQMHwbwxyGEywCMAXjLKh//IwD+MYRwJYBrq21Z0+tgZtsB/BqAPSGEawCksTgza7WvxV8CeBktO9e5vxzA5dW/6wF8vN7Old8AWp/fQItzvIX5DSjHL4QcVx/uaVp+r0tCCOv+D4szdb62JH4PgPe0oB1fBvASAPcDGK0uGwVw/yofdwcWk/dFAL4CwLDoKJ9Juj6rcPx+AA+jOsFpyfK1vg5nnOCHsOjh+RUAL12LawFgF4Af1Tt3AJ8A8Pqk7WrsW/ndwvyuHqPlOd7K/K7uWzm+esdVHx42bh++Hv82xJtDJJeH2b6WDTCzXQCeCuA2AFtCCGfK5xwFsGWVD/8nAH4HQKUaDwMYD4tT2IHVvx67AZwA8BfVn0U+aWbdWOPrEEI4BOAPARwAcATABIA7sbbX4gznOvfzyVXld2vzG2iDHG+z/AaU482k1Tne8vwG2i7Hm5nf646NMjhsKWbWA+BvAPx6CGFy6bqw+L8Wq+YXZGavBHA8hHDnah2jATIAngbg4yGEpwKYAf38sNrXAQCqmpBXY7Gj2wagG/FPBWvOWpz7aqL8BtAGOd6u+Q0ox1d47HbI8ZbnN9C+Ob7e8/t82CiDw5aV2jOzLBY7lc+FEG6qLj5mZqPV9aMAjq9iE54L4FVm9giAL2DxZ4mPABgwszPlEVf7ehwEcDCEcFs1/hIWO5q1vA4A8GIAD4cQToQQigBuwuL1WctrcYZznfv55Kryu7X5DbRHjrdTfgPK8WbRDjneDvkNtFeONzO/1x0bZXDYklJ7ZmYAPgVgbwjhj5asuhnAG6v/fiMWdSyrQgjhPSGEHSGEXVg872+EEH4BwDcBvHaN2nAUwGNmdkV10XUA7sUaXocqBwA828y6qvfmTDvW7Fos4VznfjOAX6zOeHs2gIklP12cC+V3C/O72o52yPF2ym9AOd4U2iHH2yS/gfbK8Wbm9/qj1aLHZv0BeAWABwA8BOCGNTrm87D4qvmHAO6q/r0Ci3qRWwE8COCfAQytUXteAOAr1X9fAuC7APYB+P8A5Ff52E8BcEf1WvwdgMFWXAcAvw/gPgA/AvBZAPnVvhYAPo9FfUwRi/8H/pZznTsWheYfrebp3Viclaf8bvP8bpccb0V+K8cvjBxvh/xuVY6vRX6vtz+VzxNCCCGEEGfZKD8rCyGEEEKIJqDBoRBCCCGEOIsGh0IIIYQQ4iwaHAohhBBCiLNocCiEEEIIIc6iwaEQQgghhDiLBodCCCGEEOIsGhxe4JhZp5n9i5mlm7CvnJn965IyR0K0HOW42Mgov8VqoMGheDOAm0II5ZXuKIRQwKKj/H9bcauEaB7KcbGRUX6LpqPB4QbGzL5pZi+p/vsDZvb/JGz2C6jWjDSzXWa218z+3MzuMbN/qv5f6S4zu8/M/tLMHjCzz5nZi83sP8zsQTN75pL9/V11n0KsOspxsZFRfotWofJ5Gxgzez6A9wP4cwA/D+BVS//vslrg/kAIYWs13oXF2pV7Qgh3mdlfY7HI+L9Xlz8VwD0AbgfwAyzWn3wVgDeFEH6quo80gKMhhE1rcY7iwkY5LjYyym/RKvTmcAMTQvhXLBYJ/00Ar0v42WEEwDgteziEcFf133cC2LVk+d0hhAoWO5dbw+L/Wdy9ZBtUj1Ews95mnosQSSjHxUZG+S1ahQaHGxgzexKAUQCFEMJUwiZzADpo2cKSf5cBZBKWV5bElSXbnCEPYP582izEclCOi42M8lu0Cg0ONyhmNgrgcwBeDWDazF7G24QQxgCkzYw7l5UcdxjAyRBCsVn7FCIJ5bjYyCi/RSvR4HADYmZdAG4C8FshhL0A/heA951j838C8LwmHv6FAL7axP0JEaEcFxsZ5bdoNZqQcoFjZk8D8BshhDc0aX83AXh3COGBZuxPiJWiHBcbGeW3WA305vACJ4TwPQDfbJaBKoC/U6ci2gnluNjIKL/FaqA3h0IIIYQQ4ix6cyiEEEIIIc6iwaEQQgghhDiLBodCCCGEEOIsGhwKIYQQQoizaHAohBBCCCHOosGhEEIIIYQ4y/8P3X8FGM0ZN8IAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plot(sim, component='all', z=0.0, cmap='RdBu')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "and at the top of the sample:" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAocAAADlCAYAAAAht12aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztvXmUZGlZ7vu8MWbkPFVlZQ3dVT3QA9DNUEzi4TAeARUU0QN6EZWz2gnF4Spw2gPKgXXh6lLxqlxbQFFZgBdRWkARGxAVbLsbumm6q4fq6u7qmqecMzJj+u4fGVWd7/PtyoisisyIzHp+a8Wqevf47b3fvePLHc/3vBZCgBBCCCGEEACQancDhBBCCCFE56DOoRBCCCGEOIc6h0IIIYQQ4hzqHAohhBBCiHOocyiEEEIIIc6hzqEQQgghhDiHOodCCCGEEOIc6hxewpjZW8zsi2b2ITObMLOHzOx6M3ubmR00s1Nm9rp2t1MIIUSMnuFirVDn8NLmRgB7AXwawCiAewH8Q33elQD+N4DfaE/ThLg49MUpLgH0DBdrgjqHlzY3Avi/Qgi3hRCqAO4HcG8I4YMhhDKA7wDItLWFQlw4+uIUmx09w8WaoM7hpc0NAD63LL4+IX5gXVskROvQF6fY7OgZLtYEdQ4vUczscgA5AA8tm/wMAHcvi28AcLeZ5czsP82sy8x+0Mw+tJ5tFeIC0Ren2LSs8hk+aGZfrX+eMLO3rmdbxcZDncNLlxux9BalBgBm1g/gcgDfpmXuCSGUAHwEwB8AeBOAX1jntgqxKlb5xflqM/vN+nrDZvb5dWuoEBfOap7hkyGEFwN4J4DPA/jjdW6r2GCoc3jpciPiL8r9IYR5ADCzFICnLVvm6wDeAuDnQgiV9WyoEBdA01+cAO4E8Mz6tHcCeP86tlOIC2VVz3Azey6AnwLw82fvCyHOh4UQ2t0G0eGY2SiATwG4DUAphPA7bW6SECtiZv8LwI4Qws/U4+8G8OEQwrX1OAVgFsDVIYTDZvZvAN4I4HdDCD/crnYLsRaY2bMB/DyAm/THvWgGdQ7FiphZHsDfAng7lvRZXwfw4hDCXFsbJkQLMbO/AZAG8L9CCPe2uz1CtBIzOwrgAIAygLtDCL/U5iaJDkedQyHEJY+Z/U8A14QQ3tzutgghRLtZN82hmX3UzE6Y2XeWTRs2sy+Z2cP1f4fq083M/sDM9pvZt83sWevVTiHEJckNAN7V7kYIIcRykvpONH9N+kvrOSDlzwG8kqa9A8BtIYSrsaRne0d9+qsAXF3/3ARA1ilCiJZjZtvM7FYA/xpCeLzd7RFCCOLPEfedlrMm/aV1M4ANIXzNzHbT5NcCeHH9/x8D8FUsadteC+AvwtJv3v9R92gaDyEcXZ/WCiEuBUIIxwC8pt3tEEJsDuqDf5omhHBXg/lJfaflrEl/qd3VAcaWHcAxAGP1/+8A8MSy5Q7Vp0UHa2Y3Yam3jJ6enmdfe+21a9faDcTkPfe5mH0LaqQ1ZeVpkhR1terUsWc+bZVrbF7uuuuuUyGELReyrnI8mdoZfhxQhhr9MJJK02wfAwDSNI23QfuwQv/KjbxEUH63nqMzCy7OpHwuGi1vPCEBfq43eqZv6ck13uglQtM5nsreiWx309s1s1kADy6bdEsI4ZZVNK3p/tJqaHfn8BwhhGBmqx4dUz+JtwDA3r17w5133tnytm1EPjP2VBeXav7UFqu+u7hA83l5AKiu8ur8iq7FOczsgn+yVI4nM/eJ97o4VKsutnyXi1PdviOX6umLtpnqG/LbyPltoOb3kbnxvzXV1s2O8rv1vO+2h1w81O07al0Z31lMNdE55Md6ubay3eFPP293441eIjSb49Y1gOy1P9D0dkvf/PCDIYS9F9ywNaLdncPjZ19/mtk4gBP16YcB7Fq23M76NCGEEEKIjiXxV4m1Y036S+2ukHIrgLPWEW8G8Nll03+8Pgrn+QCmpDcUQgghREdjKaQyuaY/LWBN+kvr9ubQzD6BpcEno2Z2CMC7sVSm6q/N7C0AHgfwI/XFvwDg1QD2A5gH8JPr1c5O4D+/56UuLs2VXLww4bUoADB3uujiMyX/8xf/bMw/KxerzfysHCheef6vZ69wcS7hd49C2lZc5tdmH47WERufB296nYsrCz7H09n40ZQf7HVx9/iIi3N9XueTY/0g/STcDPyzcujybSidPuJ30RX/VF01fyyU8ugqFFbdLtHZ3PwP+1x82bC/xlt78tE6Q4Wsi19x9aiLuzI+n+lXZVikQowJpDKkrwGUaMKXHjrh4tPz5WibZ4p+2sS8v5dvftlTGrZrM2EAjJ89F7O95L5TFgBCCP8v1qi/tJ6jld94nlkvS1g2YKnUjxBCCCHExsAM6da8EQSwYt/p7Pw16S+1W3MohBBCCLE5MFtvzeGaoM6hEEIIIUQLMBhSmWzjBTscdQ7bwH0//n0uXpxeXGXsNR1ArBlkgwLWEHLMesEkzWGjdViDWCN9S9LoJ9ZfFdJ+qfcUrnJxLwltfmXG2z2IzuCOV3u1SKVYcXGu1z88U9k0zY+3WS37bVTmvM420+V/ymFrm6YgLzlk6BGZpp+L2Acx8kUEKnTfzFN8Yn7W74I2sWMo4WSItvKuLz7g4qu39Lj4RVcMu3i8z2tXe7JxnvC0fMY/HLOkx86wxDCsbEsDIMrPSmRt4yeM0306X473MVPy006R5vCjdxx08dFpr5nfdJpEvTkUQgghhBDnMEMqu/HNw9U5FEIIIYRoAYZ19zlcE9Q5FEIIIYRoBZZq6WjldqHOoRBCCCFECzCzVplbtxV1DteBh3729S5mcX5p1puI1siItNZEUWMe2JGlmM2lq2Flw9Rm6iiTDjkaoMIkyaV53Eu14pdqVBP6/T1Xu/gdczLNbgdfffbzV5xvlKCBarrm+70pcKjGj6YQ3Re1Fec3JOGnn+jnoJRvR6D5wei4EgakRIO92HCe2l0mn+GJo1Muftr4QLQPsbb88Tcec/ELd/sa3Dv6/YCTgbzPk/4cDTZBbCadWvDX2ebnfVzxAxOtSgMTL2BASoYGWOUz/j7syXlj+Vp3bPK+SDWfR7v9PTPe57d5dMafKx6w8lPPuSzax0ZDPysLIYQQQogl6uXzNjrqHAohhBBCtAJrbfm8dqHOoRBCCCFECzC9ORQAcOZDb3fxwunpVW8jRa63HKdJr5Kkx+pqoLdKm9c5dZEGcYG0fYV0rB9kvV+sB1yd0fbStPM0+DzbYNLm2/TB/thQtZ+Ms3/yzAPRMuL8fGH3jS5eIM0sAORz/i/lFBn6Zgs+ZznHq6XV62x5G8YxeY0ZVS3gGAAs5/VQgTWH/NCn+UnNZmPhSHPI9xEZDU8t+PP9t985Gu2DDZK/7/ptcUNEIr/1pQdd/NRtsa7ueTu9znNLt8+dwS6f//myNzZPTZ32MekLASDM+mm1mUkXVxfm/PKL3kwatSZM30kLZ3mf76kub+ad6ht0cbo31rumC15/2dXt44F+b+LeT8+K4YI/lx//1qFoHw+e8OfzPd9zbbRMJyHNoRBCCCGEALA0WjmdSaoHtrFQ51AIIYQQokWYrewGshFQ51AIIYQQohUY9ObwUmT2r37LxdWi95qqlmL/KkuRtxTpr6pl1op4Xyj2iMt0xdqSTIE0L6Th6i35dWqka6pQG5K0fqwpZJ3iLHkUskaR10/aD+sSeRX2c+TlpxIKwzMfH73OxT92al/DdS4lPjP2VBfzdUzS1fF17KmybtbnJ3sUco6zzhYAMl1e75fu8vdJtsfrp1g/aIWeFWMAQM5vE+QDxzFrEKsJ902FprEGkW4bzNC9OlPymsP56HkBlOl8/uU3n3Dxm561K1rnUuX9X/FeqM/Y3u/iq0e8tx8AbCHvvn7zz/3U1BEfz5x0cfW014mWJvx8AKhOnHDx4qTX2ZXnvMawuuB9D6vlWAvMpLP+OBrdQ/lBrxdMD22Ntzm0xcWZkXEXp/r8/DzFhSx7Lcb3fi89P/gavuMl3u+2nSyVz9ObQyGEEEIIAQBmSKf15lAIIYQQQtTRm0MhhBBCCAEAMAPSGXUONz3Fv/9DF3P9VvZWY+81INZKcV1ZXqfc5bUj2R7vA1VL0BwxjXziAs1nDVhlIdavRDWh57y+cnHaa2DmSP83zeIqNNYpxprEaBMN4W3M0UZYY/e64/etficbGNZg8vlhXWe6iZF4hUb+lKQhynT5RxFrnwAg19+9Ypzp9x5s7NGW6vP+a6lurzUDgJDtpphqPqf9vcg+iOUEzStbkLIGcZ40hosVH7PvYZLmkKex9vGD/37AxW974RXRNjYrrE972pj3MWSN4dae+Guxp+T9a9PTXkNYO+k1nqVjvl7w4nG//PyxM9E+5k9MuHjh9IzfBj1fy/Q85md6EqztZf9RrnHeNeLPVffWY9E2u7cN+22M+WUy23yt5MwW79c40O81ipm++L7Mprw+mJ9BnaZB1GhlIYQQQggBQD6HQgghhBBiOQYNSBFCCCGEEEvIymaTsvi1T7g4UL3KSGOY86cw01NouA9ep0Y+ZtkGflVJukb2UkwXyBOO/K0svXLtxyS/xgr5bJVm5l08f8prSYqn/fzZ4z4GgOKk3+YZ0l+x/q0Zzz2G9SmsoeNtci3hVz92T+OdbCD+bNjXJWUtKGs0+fwlWBBGtbpTqZV9C3M9Pj8LQ/6+KWzxekEA6N424uL8Vu+5lh7x9YTTA355o7qwIec93Jam+XZEGkT2NaT8S8rHyLuTNIjlGmsKfbxAGsQkzeHkvL9fS3RNWYPI9YTf/Yprom1uVN71RV87/dlUF5k1httIY1gonoq2mSEfw/Kh/RQ/4uLpx7zububgcRfPHvUaRgCYOeJ9DYsT/tlYJI13I312Ejm6Lwv0XVIgfXthyGt/+7bH56Z33GsE+y7zesp+qiudnfPHnt3ptZQ9gzujfWzrGaYp/hryPdbW/NbPykIIIYQQYjmbYUBKR3RvzeyXzew+M/uOmX3CzLrMbI+Z3W5m+83sU2aWa7wlIYQQQoj2cNbKptlPp9L2zqGZ7QDwiwD2hhCeBiAN4A0APgDg90IIVwGYAPCW9rVSCCGEEKIxlrKmP51Kp/ysnAFQMLMylsQERwG8FMCP1ud/DMBvAvjQmrek4nUdRhqjFNVj5deZrP0DYp/DykIpWmYleP0kD7h0L9XA7PWaLevxuhA+DmS81gS12K8tFL2msEpaksop0tk84WuFTj/qvb4AYPqQ15/0UDx3uuhirp3M9Z2b0d3EGkQ/f6HYuEbpRoY1ho3qUbPGkHVLQFzzmn0LWbvUM+Y1Qz07Rl08sMd7nwFAbsflfp9bdrg4Neg1iLWcz/Falu4brpuM2McQaXpE0jpVSp6QkH6ckrxImbbBdZIXKZ5N8CCdXSivuAxrDnObQBN1Pq7Z6p+FVw6v7GNYWPAaufQZ71kIAKXHff314qNeYzj5kF9n8tGTFE+6eIL0hQBwivTWrL9mT9hGdemTiPXDPu6d9d9Nw6R7HD3kvRcBYGi7nzZI2vPFCZpP2vVCycfZy2NNbWHYH/vWHv+8qAZ/jbl++Xpim6R8XtuPIIRwGMDvADiIpU7hFIC7AEyGEM4+4Q4B2JG0vpndZGZ3mtmdJ0/GxcyF2Ogox8VmRvktNhub4c1h2zuHZjYE4LUA9gDYDqAHwCubXT+EcEsIYW8IYe+WLVvWqJVCtA/luNjMKL/FZsJs6c18s59OpRN+Vn45gEdDCCcBwMw+A+CFAAbNLFN/e7gTwOE2tlEIIYQQYkUMhnQHvxFslk7oHB4E8Hwz6wZQBPAyAHcC+AqA1wP4JIA3A/hs21oohBBCCNGAs28ONzpt7xyGEG43s08D+CaACoBvAbgFwOcBfNLM3luf9pG12H/1wB0uTvWR+S4NMAk84IQGcuTyXlwLIBrkEqpsrO0Nqa3LC+m5TekhL7wHgNSwNwCu9ngD4BrF88Hv80IMVPPwx9U15QecFE754vODB+6Ltjm1z5vJTj7shd0TB3wx+gIZabOJNou2gcZG2dFx5fy5+eqzn+/iF9/1HytvsMN4f48vQs/niIXsjUTrSQNSeulhyANQesf7XNy/2+fr8HW7XZy93Bt1A0B6+5Uurvb6nyDLBW92XDX/eGtGos4Da/hIeRvsZ5ZOxXthbXqjdwp8K/KAlcWEHOcBKGyKXaTBDTxA5fv/5Bsu/vuffkGDVnYOH73DP2du3OZzbYwGoPSU/KC3zOQhF/PgEwCYf+RhF5/+zqM+ftAPvjuz3z+3Dk15o+djCYOKeAAKP5NXHjbWLCvfBdwGbudwLi6esI3W2UmDWBan/bHzgMwRGvzoh5YskaPvyB4aGLal21/zPYPezP6Pv/GYi3/uBbsT9tIazICM3hy2hhDCuwG8myYfAPDcNjRHCCGEEGLVmFnL3xya2SsBfBBLVn8fDiG8n+b/BIDfxpPyuz8MIXz4YvbZEZ1DIYQQQoiNjgHIZVYuT7uq7ZmlAfwRgFdgybnlDjO7NYRwPy36qRDCW1u1X3UOhRBCCCFagBlaPSDluQD2hxAOLG3fPoklhxfuHLaUS6pzWDrxWDQtnfKnINXLmkP/F0DIeS2VLXgNXKgkGFyTxhANNIbpAa8PxOCY31y/12sBwEzGay6Oz3mtyKHD3kz62KzXgUyRtqRUjRUuPaQ3Gen2uo8rh3y7du7Z5eKB8adE2xzd+W2/j3GvAc0PerPZif1e28PF6nMTseazQsdWIy1Pim7kVHblnwTuePXLomnP+cJtK66zXvxuX3yOp8o+/4pV1hj65Qs0ged3JTz4CoP+vujb7o2IB6/yuTFyo9cU5q59totrY14nCQCzXcMunlj0x7U4Sxqt4HV3fFnzCT/9FKicVRctw9Wu2CA8yR2Mb6VSxreT28GnlzWhxVKsV5tnw2TK+blFH1dKsdHwcl76B/8aTfvyL/6XFddZD95320PRtOfu8s/sLd1eB95v/pmcnvba6PIhr3tmg2sg1hieuv+4i0/sO+XiR+Z87h0u+niukRAaQJbyYIB1vWmOG3dI+N4vssF6ZeUiA0cStJL8fOFtXEnrhAbHnlRIwuh7N5PJu7h/xMecA1cOew3iu774QLSP93xPrHO+EAyG3OpMsEfN7M5l8S0hhFuWxTsALBfjHwLwvITt/JCZvQjAQwB+OYQQu7mvgkuqcyiEEEIIsVZcwJvDUyGEvRe5278H8IkQwqKZ/TSWqsq99GI2qM6hEEIIIUQLMEv+VeIiOAxg+c9wke9zCGF5LdsPA/i/L3an6hwKIYQQQrSANTDBvgPA1Wa2B0udwjcA+FG3T7PxEMJZvcRrAMR+TKvkkuocWi3WSwTyS4r0f1k/v1b0RcVDwS/PnoaJ7WDdYq/3Z6uRf1ulf9zFR2bj4zhw0mvv7j3ui53fc9AXfj940h/H3LTX6lXKseaQi4kXev252UNedt91lddOfteuoWibVz/tFS7u6ffrpHv/zcWZngddnO055uLZo3Fh+CLpEKullR3DUqTdMYorxfj8dwpJPo+sM2L/NPYtjDWI/roPJGgye8e8O9nQVd6Lc+TGa1ycf4bXr5XHn+big7OxJu5x0peeJi+/MvmlFbJeIztKGtmtPV6XBABp8+sU6AlpJa8xtqpvQz4Vj1LM5r3+skYawsW838kEtZt9J5Mo0XUvkg6sRJq3Cs1nHW6u2plfDeP9XfG0Pq83G+zy5y81dcTFtZNeilU+5DWGkw/FUi32MWSN4UOzXtd4mJ4RrN1jPSEAbKE82EXJN9bvj7N7xOvocr1xPjOlWZ8H86e9Fv04eRI+QcdxcjF+9rF+8iDdl6yZBZ07fr5muvx9CgDZHn/dUz39Ps757+HBAa9330F5c/UW+t5uIa02wQ4hVMzsrQC+iCUrm4+GEO4zs/cAuDOEcCuAXzSz12DJK/oMgJ+42P125hNACCGEEGKDYWj5aGWEEL4A4As07V3L/v9OAO9s5T7VORRCCCGEaAFrYYLdDtQ5FEIIIYRoEa1+c9gONnXnsFj0eop0rjdaJpDPoZFmiDVFlva6DtYxGusrAATSDIWs12dVqSZsucdrDo/O+DY8cMrrngDgrsNTLr7zwBm/jcO+nugkaw7PeO1eec4vn0SWdB/Htm538YEjfhuPXh+3+3uv89q0Z+56losLGX99Rvn6kCdWOuePAwAyXf6alea8PqhKnm+pNG/Tx0nangdvep2Lr7nlM9Eya8Fv965cNxloXDebNYasQeyhBXpI6wQAQ1d4PenoDb4OctcN3+3i4vgNLv7OCX+v3nnE5zMAPEw6Wq4f3E0+nDuH/X1WGfE6o+4E7eRAnvR+VNHWFn0bUgukcQ3x+U91eS3ucB/5llKd8+mCz68TcwneqQRrDtnHkDWGi6QlC5Qj5cVY8/n8937Jxf/xG6+Ilmk17GvInoZAfM3yZa9NTc2cdHHpmK/FPP2Yf2ZMPuqXB+JaybGP4coaQ76HruqNdXVXbvPfT2NP998DQ9f452vvDj8/P+jzLInFSZ+vs4f9se540Oszd9zr5z9yzJ9bANhPekvWIPK5YQ1tls5tnrSVAJAf8sc2NOCvYY606vlu/zwayPtn1tae+Pzf/A9+DMf7XnVdtEwzpFqsOWwXm7pzKIQQQgixbrS+QkpbUOdQCCGEEKIFGAzZhCovGw11DoUQQgghWoAZkG2ilGGns6k7h0b6n5CNtQxcOxnkexiqXk8ReSVWm9AcRvWZvRaq1uP1EmeKXu9zlDQdB6fi+sEHTpAH3Bmv75ud9OuwxnD+tNealGa8ZhEAalQ3OpXx56o44euNLs76Or+N6rkCQDbla/A+g/zvclTLergUnwsmnfW+WsUJr2/j2stc+5M1h6ls7GVXWWisC1sLWE/ItVIBgKfw37SsAWJfw2HS8vXv9FpTABi82vuKFZ72HBeXL/Na0m8d9ZrXz+/zPnL/+chpMLOT/rpxTexCb8L9vYxR0nnVarF2ks+FLfr7Kj3n21Wj+ygs+jYCgJEXambUL7Nl9AoXz5e95vBol39MF3KNH9tVyoMSaQjL5FdXJc2iJfwsljRtrRkib0r2qgSAfr4/p72GrXra11JePO7jmYP+uTX5qPeEBYBDU97/j2slN9IYXkdejNff4PWCALDzv3j98OhzbnRx7iqv0w2D3v825GNNPZOjfO6f9OdidL+vdT90xz0u7v3Xh6NtZr/tdYn7Zvy5ijWI/tz1TpF/bsL57x7116h7m6+znh7x3xupXv+d2j9wmYvH+2K/zMuG4+fBhWBozpu009nUnUMhhBBCiPXCAGSlORRCCCGEEAAAs+hXjY2IOodCCCGEEC1Abw43ABVSV1k61hmk6QwE9i2kOFRYg+i1PCHB5wzkpVjLeW3DfMVrMuaprvE8eZRxDdkk0uSzlCFdTrarl+b7NlUysbaHNYccs05xjmqYnur2WksA+M8+r+HaOeTbMdzttSO7x7wuJ3u5X39gwetdkkjlvKalPLdy7eVGtZYBIJ1tz63EdZOrK1saAmjsa9hLuVMY9PdNkuZw4LqrXBwu9/qohyd8rrDG8F/u9prXU0fiGtllura5bt+uQdJ9TVIusRwzSTSez/hpqVnfDtYYlg/t9/NnYr2UUX32DC2TMX++R/q9Pop1dn1dca6xdQZrCGuRDyLF5cZ64HYwTJ6Pfbl4FGgeXsOWKnrNYWnCa+Lmj/nn1OxR78c6cST28jtGumTW0XGtZPYxZI3hFa9+RrSP4f/6Uj/hyme7cDJP2vQFf83mFhp/L/Tk/L07vMtrhfu3en/SbWM+F3N9X07Y6t0uKt/j7+37p1fWIPK5HU04//10jfroGubH/TXOjXmf1PyAz5GeBI/TrT0ra5abxQCkpDkUQgghhBDA2dHKsrIRQgghhBB1NoGTjTqHQgghhBCtwMz05lAIIYQQQiyxpDlsdysunk3dOZyjgR1JF4wnpcmwOksjVjIZL4JvNIAFAAIJzgMZbVepnaxlzWd8mwby8WXbOewHexTJcDqQWJ8NbdnQOtMVG6ouznoRcK1Mg3PSvp1G55IHFACxsfHDx/0AgGtGvYHw2HZvftq73Quoc7PxgIDekr8mKfqrrjTtjbUbGVpbwl+F7RqQwibYjSXpSabXPu6nASm9Yz63+vd4810AyF3xVBdPdHkB/tf3ecH4v9/vTW2P7PcxD2YC4nwKI9tdXBmMB5wthweg5DPxdeRzkZqnwQ3HDrp4/jEfzx2LzbuZvjM+R3tzvt191466eGuPvzf7mxiQwvCAkxKbYJOZfKjFA1SqldaYBK+GkW4/IKU7YSBBasGfzzDrByNUJ/wAifkT/prO0ACIUwlm/WcaGPhvoWfyldv885MNrqPBJwDCdS9y8eOL/tj3PeEHZTxGz84T040H423t94Mudg/6a3rdlkEXX05t8k/fJUoz3tB+9rh/np4kA/YjNACFz23S+R+mazRA17CfrjHnQGrAf6/0ZIeifQzR4KcLZWm08sZ/c9gRR2Bmg2b2aTN7wMz2mdkLzGzYzL5kZg/X/42vphBCCCFEp2BAOtX8p1PplKZ9EMA/hhCuBXAjgH0A3gHgthDC1QBuq8dCCCGEEB2JwZBNpZr+dCpt/1nZzAYAvAjATwBACKEEoGRmrwXw4vpiHwPwVQBvX/8WCiGEEEI0ZsnKZuOLDtveOQSwB8BJAH9mZjcCuAvA2wCMhRDOVgU/BmAsaWUzuwnATQAwtn0n/v3RJzU/o6RXySTocvi1Ljub8zrxfH8Ks5n4lDZKE94HG3Ru7WEtRGwmzaOj2Ch3kM7FAdIx9ZAWZXbIa/0AYH5mp4srC15rwjol1jGyaTGAqMzQ5Lw3Kz1N8QTpVwq9XtuW3rIj2kd+zutNApmIp7v8sVdJG1mrsil2/Ndepsdrd+Y+8d5omQtleY4PWQa/2/eUc/OqwWsOm/k7lE2vC3Q8nH+FEZ9vfbu2RtsMo94s99C0v25f3++1eMcPek3Q1OGHXMyG6gCQ7RlwcdfAaLTMcnKkKezL+XsiyQg3VfTtqpGWafHoIRdPP3bUxVOPx5rDGun9Fie9fipL90X3tt0uHhjwerWBfKyNKuRIj0laVDartaeNAAAgAElEQVTFZo1hecG3qVbx1w8AjJbZ/ZZPnvt/bmTPs3n5Zlme38PbtuNPbn/s3Ly92/01LyRcMyt6jRsbkfP5XjjtnwfFCX8ukvSFrO1l0+tdBZ9bY0/3z6XR59zoV7gyPl2sMfz6QX8ctz/q74kHj9JxzDTWHBb6/LPumvE+F0/sIVXhZV6DuDuh3aPP8brbsQe9of2uE/574iTpXfncJp1/vkZ8Dfka5ygHrORzJN8VqydZz/q+2x6KlmmGzWKC3QnvNDMAngXgQyGEZwKYA/2EHEIIABJrP4QQbgkh7A0h7B0cHklaRIgNzfIc77F04xWE2EAsz+/eQT3DxcYnbdb0p1NZdefQzHrMWvoNdQjAoRDC7fX401jqLB43s/H6PscBnDjP+kIIIYQQbcew9NNys59OpWHn0MxSZvajZvZ5MzsB4AEAR83sfjP7bTO7qtE2ViKEcAzAE2Z2TX3SywDcD+BWAG+uT3szgM9ezH6EEEIIIdaUTTJauRnN4VcA/DOAdwL4TgihBgBmNgzgJQA+YGZ/G0L4q4toxy8A+LiZ5QAcAPCTWOq4/rWZvQXA4wB+pNFGzLyH2zz5BzajOeTXvLxOjoSmXRkfhxDvg9dhL8QCjVjK5f2L2Tyt35uNX9yyvnJHv9cx7SI/q6u3ea3J46e8LuTQGe+hBQBnprzuo7Tgi7iXFlf2Assk6IUKvV4D003aqTLpAxcr5OvX6zVJ2cFYD5ce8fqTLtJGZrq8zqw8R55vpDlM8jlMkc9hqK58LloF52vaYvVFI1/DXtLm5Ul/2k1ek5nRbdE+qgPe+/CRR30+PUr6qJkTXpe0OOV9EGuV2GuSvTfZ95DvG9bZsj9gby7hOk57/7Qy+RrOHPR+jJOP+nZPHIh9NitF0hySH11+0B9X/prHXDyw9VoX9+bi+z/Js3E51arPi0rJ39+Vor9eSeefNcS1rieXCcmqn4smR/caa74BwCr+fFZJC833M5//4pzXV85WYrdQnjJA53uM7pmha7wHZ+6qG1w8mY9/PmcfQ9YY3v3QKRefPuZ1douz/riTyPf6e7mRTpG9/4Z3xe0epGMbuuYeF4/d4XW5vbM+tyboezrp/PM14mvI17hGOZChHEnKoy7yEx7qzkXLNIPBkGo40qDzaaZz+PIQQqRODiGcAfA3AP7GzC7KPTKEcDeAvQmzXnYx2xVCCCGEWE86+Y1gszTsHCZ1DC9kGSGEEEKIzc7Gf2+4CisbM9sL4GYAl9fXMywNJL5hxRWFEEIIIS4BDI3LWW4EVuNz+HEAvwbgXjRXwnXdqQVgYZk+bIG0Ykm+lOxH1EVaEta8VCNNoZ+fZgMsxH9FsOcSa2ZSwbc7S16Kfdm4xumWbu9FN9bj19lFmpjJBa89ObrF656OzsZalBOkTzkx7XUe7FFYStCOMKwL4xrRBdJXkiUWquaPM52L/RnTA14nE8jjLUPatXQ3aQ7LK9daBuK60si0pk5nIzink6wR2NeQ4y6eT76auT5/Tfh8AkARfp3T8/6cFUlnVJ7z+qokjRuTynL9b3+tc12+DVtJd8u63L4kz7wZryEskzZy9glvmjD1uD+Ok0dj3RdrqMaptmz/Tr/PrScPuzhb9fnYnaA55i+jQP6XNWoDn+8qaRCrCdcjTVpdp1usteYrgevS8s9zmYRnuFV9W8OiP1/sW1ou+vNfpO8J9t1Lgr1Bu0f8M7l3h/c5DINek3tmIdYkc61k9jFkjeHU4cddvDDt8yiJrv4tNOVyv0/yQbx2u9eVX7clfr4O0LHxsfO5KRz39whrDpPOP18jvoZ8jTkHOEeS8ohlu9wXWA2b4FflVXUOT4YQbl2zlgghhBBCbGDMLr03h+82sw9jqc7xuW56COEzLW+VEEIIIcQGpJP9C5tlNZ3DnwRwLYAsnvxZOQBQ51AIIYQQlzyGZHnPRmM1ncPnhBCuabxY+6jWAqYWVjdwmi8iF8xmfU8v1WdlbQF7GgIAuOZwyWsuUvPeWy3Mex0T63ksSc+W97qOwbz3MRwo+BqZ20e9P+DOfq/nmlqM6zdPklbqFGkMp2n+TMnH5WqCBx+dQK4bO0peUywT4y2GbEL9ZvJCTFPdWD6fgTzfQkKd2Yak1qbM3dKD58mY8zfBuq9hLWX2PcxQndhMjz+nVoh1R6wTmqP6qNXqypo09tDjGADyvb4eanefX2Z4wLdzK+mnBsnnkO87AKiQ3m/2sPeWmznidV8TFD9RjHOlmJD3yxkhDVZt1nslpsvk0ZZu/NiuRbWV6XqwLo80hrUEnS3XTs8s2wZrHC8UM2B5uhortkNCHvE0ame17J9Dga4H5y7XK0+C75lcr3+G5Af98zfkvaZ7biE+jhPsv0gab/YxZI3hwsSxFVqcDNcnL874dnKb5kpxuwP5j/Kx87nhc8cknX++RnwN+RpzDkQ5kpBHnGsX/MuwJY9v2GisRjf5dTO7fs1aIoQQQgixgTEsdSyb/TS1TbNXmtmDZrbfzN6RMD9vZp+qz7/dzHZf7HGs5s3h8wHcbWaPYklzKCsbIYQQQog6rbayMbM0gD8C8AoAhwDcYWa3hhDuX7bYWwBMhBCuMrM3APgAgP9+MftdTefwlRezIyGEEEKIzU6Lf1V+LoD9IYQDAGBmnwTwWgDLO4evBfCb9f9/GsAfmpmFi9B7rKZzeAzADwHYTeu950J3LoQQQgixaTBrtZXNDgBPLIsPAXje+ZYJIVTMbArACIBTuEBW0zn8LIApAHdhmZVNJ1ELAbPLhPC1C+g0LzdgBYAFMo9lbXkmRWLbbON9GhnO1s54IXH1xCEfk0A9yXCWB1Wkur0oONU35OLskDcqHenzwuShntjouLbFb2Om7AX/syRWniNz08UEU2wepEKnH705P7Aj38iYNB0PZAg5L5hO9ZJ4mXYacn5gQzQghcXO64iZH2DCg02ShNCNTLA5znSxsTgNgOKLlACbx2doJFG2x5vr5hb8YJNMLjZ6LwyNubh30F+nncN+nfFen58DeZ9LqanT0T5KNCBl7qhfZoaMiY8v+lw6uRjnBovpeymHKzSQq1ZaeQBU0vgWNpyvlFc2vY5MsBeLK84H4kFCbpmkgSJCXIJYCLDV9T1GzezOZfEtIYRbWtysVbOazuHOEIJ+WhZCCCGEOA9WqzRe6ElOhRD2rjD/MIBdy+Kd9WlJyxwyswyAAQDxX7+rYLWjlZ9+MTsTQgghhNi8hKU36c1+GnMHgKvNbI+Z5QC8AQBXq7sVwJvr/389gC9fjN4QWN2bw+8G8BMarSyEEEIIkUAILZUe1TWEbwXwRQBpAB8NIdxnZu8BcGe9rPFHAPylme0HcAZLHciLYjWdw1dd7M7WmmoAZpcZL7OerRlzU9YcdpV9XCb9UJb0Wr1JLsTmp3ER8MrECRcXDz7m4vlj/u1waXo+3geRIrPuXJ83Lu4a8Zqv/KjXGKa37Ii2maVpQ31etzjQ7TWJtUFvPp1kBrxYWfkasYSOTcZZYxfSsUF4YP0a/bWWytBt0OXPVSh5w2A2AwaQqANdCwxA17KcY0PZBI/aSFPYxRpDMnpPsTM5k2AKzvsYIN1iT7/XB/YMb1txF1ky1gWAvhE/rXfQX9fLR/11GycT7H6+N2dirfbiKX+vzZ/whvTzp7w27wyZfc8m6Gobkebz3+0N6EPWH8d8gkH1DOkWK9SuCumcOW6kSQQAI2P35fdBayyw69+ryzYWeMvW+PnKBvTprM9FYwPrSLebNJjAt4OfZaVZf08sTnptam7Rm6X35PzzFwC29vvrXKD8zffSM7zfP3+bgdfhbfI+uU09Cd9vRsfGx87nppEpfNL552vE15CvcVSEgHMkIY8CSPt7EUltLdbghhC+AOALNO1dy/6/AOCHW7nPhp3Ds8OhQwiPN1qmlQ0TQgghhNhYBGB1msOOpBnN4VfM7BfM7LLlE80sZ2YvNbOP4cnfuoUQQgghLk0CWq05bAvN/Kz8SgA/BeATZrYHwCSALiz99v1PAH4/hPCttWuiEEIIIcRGIMCqG//NYcPOYf237D8G8MdmlgUwCqAYQphcec31J4TgfAlnSYeT5LNXJWEBm1f2kXaKNXFdpM8aKsSntEKnOUN6h1D0BdVZYzjxkB+1PnvcLw/EXmlMrsd7lOX7fdy91WtgesZjn8PeHd4LMbt1u4sz2y7z84e2ujhNmkQA6M6vXKC+Si+3WbsQqVPS8fkP2e5omiNF62S83spIg5hoU1Dx01in2CpSZigsy7lCmvLXYnVHI1/DFHkQsp4nVL1+rVaM8y9vfplt5DF42RZ/DmfOeM1bJuv1rJkEbdPAiL+OV9A2t5OucZD9GotTLq5NnIz2sXDaawwXp30usHcnexgmwee7n3wOu0e9djJN981iyh/XRML5n5r3uq4y+S9WKR9r5ZU1hkm62kStbYsJ8N6y/MiuJJzuDHmbWt6fr3SXz8UsPaML9Azn65VEseobNn/a5/PsYZ9b/ZNHXTy8axeY3aShvWbcPxuLM2wvfLmLugZG0QjWGI5s889b3ie3abiLtHwA7KQ/Nj52Pjd87pik88/XiK8hX2POgUA5kpRH3KzyherIWzwgpV2sZkAKQghlAEcbLiiEEEIIcSnSwT8XN8uqOodCCCGEEOJ8hNWaYHck6hwKIYQQQrSKS+nNoZl9CcD/GUK4Zw3bc9Es90ianPc6GtYgAo01h7MLXmMxT/5h7MnE/m4AMEzT8lmv4+C6yJUF3272Vps4EMs9Z2f9OqyF4nb2ZFfWPfWNxx5wfdu99rF3l1cY9F3ma0Lnx3e6mDWJAJAhv7sq1XQ21iTmyAMu2mAXT0Fgn0nywDJeh3wojfRYlqQnIZ0iFmJdWCtIwWty2LMwSa+TpUnsjZjOxTqi5VTLpKecm46WSc15nezOfl8r+cbLBl18hvL1NGuIEmpoj5OP4RVbvV5qB2kO+/N03Wf8fVOdmYj2UZrx1600RzWIw8r3FZ9bABim87tjiLSRV/l7ILvrKS4+WvTn//hsXNp+hqYt0jqVoveiqzbwOUxiPTSHlVoNE8ue2wsVv0/2mQWAfMbrzVKkEc72+POdZz/BHv/87Z2Nz0UjP8vj0/7873jwiItH93/bxf1br4z2cd0Wf49M7BmOllnOg+RJWJyJvUEZ9jFkjeHzaJ/Xka63fzGhHjkd2wQdO58bPnd8p3PtcSC+RnwN+RpzDgTKkaQ8KjXQQjaLhc3x5nA15fPeDuD3zezPzGy81Q0xs7SZfcvMPleP95jZ7Wa238w+VS8bI4QQQgjRudRqzX86lKY7hyGEb4YQXgLgcwD+0czebWaFRuutgrcB2Lcs/gCA3wshXAVgAsBbWrgvIYQQQojWcna0crOfDmU1bw5hZgbgQQAfAvALAB42szddbCPMbCeA7wXw4WX7eSmAT9cX+RiAH7jY/QghhBBCrCUWak1/OpXVaA7/HcAeAPcB+A8APwHgAQBvM7P/EkK46SLa8fsAfh3AWQHECIDJEMLZH+4PAYiL/S616yYANwHAZZddhrd+155z8974sTvcsqUmfA6ZAumFiqQ9yUS1lWP91mi3/0WcaxKnyNeM6yDXSAtRnItr2x4hPeUU+bGxVoq1acPkkzZ2PK7fPHTC67EGJrxuqTzj1+mf895q3Qm6Jq6EnCI/qhrpAflq0WEiqYpjOuW3mc6TLxrpQ1hjGLI0P6lSJN3kuctvjJe5QDjHX/f4fefmfXbbU92ySXVLWQbXRZqeNHkKpmiFapF0d7Ox5jVLmsOxUZ/TTx/z2iauBXygt7FqZOew15vuIg+20W7SjpGuNrXoa76W530MALVSA79QriFN+0jSS+0iPeW2Z4y5ePRZ17m4usXr0Y5M+fN/4ITXDwJAcYbqtZPmtVEt5dDEz1us1T31uXc8Oc/eeVfDDZxvu5TfN7/sSc3lVx/x2meuxQ4APaRDTvV57V5+0GvxukZ8LhZIAzo8EXuUHqN8XaDvjSdI47njXu/1N3SHl+pvG4v115df9yI/gXS6QwWf39du9960J6ZjLSrDtZLZx5A1hpfn6btmX3yZT9GxHadj53NTpkvIumnW6ALxNeJryNeYc6BGOVLkLw4AZ4r+WH/6ebtd/DPRGudjc5TPW81o5ZsA3J9QQ/kXzGxf0grNYGbfB+BECOEuM3vxatcPIdwC4BYA2Lt3r+o7i02HclxsZpTfYlNxqZlghxDuW2H2915EG14I4DVm9mosleXrB/BBAINmlqm/PdwJ4PAK2xBCCCGEaDMBoRL/urfRWJXm8HyEEA5cxLrvDCHsDCHsBvAGAF8OIfwYgK8AeH19sTcD+OxFN1QIIYQQYq0I2BQDUjrZBPvtAD5pZu8F8C0AH1ntBthLrZzgYxRIO8LxNGmI2CuRNYu5BM0Rex8OdXl9xOj2q1zcv8e/pO3b/riLu/bH/mxV0gxOlX3SJWnR/PL+3LAXFQDsOuQ1WqyFtBRp2bq8jqxrxGtkACAsei0UC3QD6ZxK1KxFakOSVZUZe9P5+VnaR5pqMfPyiSLiNgmLe0a8Zigz3divLkP5yDHDerTqbKx5q5055ts16D0urxjy55TrnG/pZR+y+HwO5H07I19D0iqlSqSbpVxr5q97Pje9pI3cRT6HXQm11bdc72ve7vzu6/06z3qxi4+Zfz48eOqMix8+Fp//uWmvkyvNey/KKh97gy8lvpcBIJVZfzexCdKBzZXz0TKDBa8/S/cO+Jg03d1bfa72bfe6xtFDsRaVdXCs8T5JtawfoWvU+68PuzjX9+V4HxTvvvLZfv4u7wHL+sA5fjgm0EP6Yq6VHPkYksbwzL/E7T5Ex8bHzueG4XM7mqA57Nvu74nurUMu5mtslAO1Lp8j8/PxuTo935q3fSEEhPLGf3PYUZ3DEMJXAXy1/v8DAJ7bzvYIIYQQQjTPJaY5FEIIIYQQKxA2h+ZQnUMhhBBCiFYQAkIT5Sg7HXUOhRBCCCFaQujosnjNsqk7hzwApZwgjK2SMD4aZEGCc95GkrE2U8h6gS2bmfaN+gEp+ae9wMVbj3qR8PwpLy4HgLk7j7q4SMdxctFrIEo0kIbjpAEsPEhlqOi3WSs3ELlnY0G7FahAOgmHa3kvRJ5f9G1YaGJACnmsIk3XNJ3yx8rG5lmKef2lddpzK3WROWwSNbqWbHqd5oEclK98T5TnYpPg3MQJF2emveh/+9anuDgEP0CFB23NJ+RSPuPbtbXH30fdWX9drNRgAErCoItcv29X75jPz1Tar8PnhoXzALD1GVe7ePAl/83Fpe03uPjBQ17Qf88hbzo+NREb1C/MkQl20W8jNr1uNCAlHhSQzreyWmpznJjzxs5z5e5omUV6nqa7eECKLzrQvc0P/egd9wPlhrbHA1K2UfEDHvA3R/fYfhoMmf22N4YG7o72UZrxxuWjzzno4sGrfJ4MDI67OHTFucfYos8LO+m/N0r7v+1iNrjmwScAcD8dGx87m1730Ai/bTxgM+Ee4mvE15Cvca3gB6wsUrmFmVL8Zo9NsC+YEBDKenMohBBCCCGAS88EWwghhBBCrIQGpAghhBBCiLOcNcHe4GzqzmEXaadYlwMAJSoKXiUNEZtip8jkmjWI+5NEb0SetlEgLdXVV3h7x8GX+HbvTthmZNZ793EXH6Zi8tMNtJL9CWbe42QA3Lfd67H6dvpi532Xjfk2ju+OtmmjvgB9pc+vM7Hgb7JZMuterJCJebQHgBWCJBtrqDHMpVeOASCfMG09yPX4axIStKKsi2PdHG+DNYhMZSG+j2ozXheHSZ9/GdKSjvd5vRTfE/PlOD9Z6tmT5fuILiwZFbPGMN3r8xUAureNRNOW03eZP/Zcn78H+vf44wJiDXGF7u8HTntd3V1Hplx832FvaD19JtYcL856o+wK6y1XaXqdzsX6wnQbTLAP0rHuHozbNURGzl09/hpmRijXxsgE+zJ/7gZPee0fAOyk5yfrrw+SgTJrEPfN+GtcvsdrdAFg9rjXko49eMTFQ9d4/V/vDq+zyw/6eyyJxUmvp5w97PWCE7TP4/f6+WxwDcQaQz72Lnqe7iCN6M4Bb2w+uCe+L/m7JD/mr2marnGl22sOp8kg/NR8/AybSJh2IQRZ2QghhBBCiCeR5lAIIYQQQpxF5fOEEEIIIYRDbw47my/87He5+Dm/+U/RMhXyqyqRhpA1h0yp6LU6lYTi5w+RPqVC26yFQPO95uLaa1/s4uF+LtEOdO/6ql/mmn0u3r7f62zYK7FK54E1jEDs+dZ3mde8DF97uYvz1zzDxamd10bbrIzsdvFJ8k6cJM3hNF2fxSZ8JrNp1hT6a9ZNGrsK+R5Wa0ZxrMfc0r/+HnAA8NwvftnF3/yBV0TLlGb9X7GsOcwU/PFnC/7ap7ONHxNhweu0qlPemzPT673nsuQLuaXX59J8Nj7H1bCyH2UalAvmt2Hd3istvXVntI8CeXF2bffHZRmvl0oNbfXxdu9ZCgCL5GN6YNJrm/7zsNdr3nHAa+AmT/g2FEm/BgCVol+mVl6dr2GK9ISZQuw1d/oLN6+4jbXgfa+6zsV/+c0nomV29Huvz4F+3/ZUn8+tzDavc+6nXF2ciH0OF6f9Ob+S9Kycm4dJy846vPun42vIXrS76LqP3eE9CbtH/DMn1+tzMwl+Fsyf9t8Dx6ldT9BxnEzwCmYfw1hj6O/1K8mfdPgqrw8c3OOvFwD0797mYr6GNbrGi1mfA1PT/n54Yir2ar35ZU+Jpl0QIaAmzaEQQgghhACw9LNyEwNTOx11DoUQQgghWkAIAbVy/IZ1o6HOoRBCCCFEKwhAraTO4YYik6BjWiyS3o/81SpUUzNQQW32Bystxtoe1jFyPEv6lWmKJ4per3XN6FOjfWx5mdf77bzhERePPfGQi6tUC7c27z22UrlYv5IiX7j0lh0+3rbbxZVBr+maTMV1UU9Pe23GKfILm6B6l1xzlyWhXEcZALKkseumPOD6zH2kQTTyyyw0lva0jXx/PpqWzq7sW5hhjWGX15+lck1oDqt0n5AGMcx6775UlnKB9IHdXbFnW2jks0e6upDx56JG3meWietSZ0e2+21Qu0LW67xY63SqGrfxIOl77znmNW23k8bwEM2fm/b6qNK8P5cAUG3ga8i1klljyL6G2Sbq9LaDh0/GHoRbyaezn+7fPGsOt/htZOe8j+RgQu1w9vaM/ET3nXIh118/XFzZBxEAjiysrO/rJT/BwnF/HIUmvFaLtN8iPfvYv5H1hElwrWT2MWSN4dbrRl08co3X7Q4+ZVe0j+zOK12c2uKXqdA1Zq36UdLpHp2Or3GrCCGgVtWAFCGEEEIIASwNSFmnn5XNbBjAp7BUG+MxAD8SQphIWK4K4N56eDCE8JpG245fpQkhhBBCiAsiVGtNfy6SdwC4LYRwNYDb6nESxRDCM+qfhh1DQG8OhRBCCCFaQggB1fUbkPJaAC+u//9jAL4K4O2t2PAl1Tn8xs0vj6bd8PbPu7jGnoTs31ZZuf5iORXru8oLXr/D9Zi5vvMZqgF96IzXA964M649ec2o9yDcPnKDi0d2PsvF+ZrXXFg59t1iQtZruEppr9maIr3lFOk+TszFNWG5xiXrLVljyH5irO1hT0MAyJPmcLHCGkJ/zWskZOS6vzuGOlOPBQBP/YvPRdMe+tnXr7gO62YzpDlsFAOApSnvWYNY8vmWKvuctpLPrVQqfjTVeLdpqivN63R5X8NqLta8RpDGsEI/rsySj+nEnD/Og1Nx7dmHT/tj/fYT3tfwwFGqdzvpzxXXhK+WYr3Uan0MWYOYJV/DU5873wuI9vKe74m9Uv/4G4+5eLTbHyt7Vw70+zq82Z3+2VdIOL8jtZXf8Bh7qe73v+z1Tvk8OsZ1vwGcIX17iZ5DE6SH57gV8E+J7Fk4nFB3fRt54nKtZPYxZI3hyNP2uLiwx+sLASC703uFVugaTgd/jU/Or+xr2DJPwyQCVvtGcNTM7lwW3xJCuKXJdcdCCGcNMI8BGDvPcl31fVQAvD+E8HeNNnxJdQ6FEEIIIdaKC7CyORVC2Hu+mWb2zwC2JcxyrvQhhGBm5xtCdHkI4bCZXQHgy2Z2bwjhkfMsC0CdQyGEEEKI1tBiE+wQQvyTZx0zO25m4yGEo2Y2DuBE0nIhhMP1fw+Y2VcBPBPAip1DDUgRQgghhGgFAaiWK01/LpJbAby5/v83A/gsL2BmQ2aWr/9/FMALAdzfaMOX/JvDb3/ge128+y2fdDFrDFnvwzVMkyiTbrE077V6C3Pex3B+1mtgZie9Vu+BI3Htz6u3ec3QFVt9PNbrdSBD5EXVTfVzkyyzylXWA3q/tVnSzEwtem8v1hMCwAxNK67SPLRAHnysDwSAXtLElEnLUw1eR8PeiS+9Oq71uZF4yoc+7eIj7/1ZF1dLDWov93gPvFRSrWWqOQzWIBJG2lFjj8JanAcWSBtK80tVvq5+fplqYnON86Vt+K0WK74dp8mHk/3THp+MdbUPk2/h/uNelxhpDGkflQW/fDPPnFR2ZU/IXLd/5py89dcabrNT+bkX7Hbxp+894uIC+Zpmer0WtWfAe1tmL4/1m6xWbaTTZb/RwqNeZzp6JNamnqLnJ2sQ2YOQNYmsx06CNdo50hT20vOTNYajCZrDoe3+u2ZwzyDF/vnJPoasMcxe7utpA7Fn7lzOX8PjM/6eeeCU/879qef4WsxrScC6ls97P4C/NrO3AHgcwI8AgJntBfAzIYT/AeA6AH9iZjUsvRB8fwih8zuHZrYLwF9gSUgZsCTG/GCz/j1CCCGEEB1BCKiVy42Xa8muwmkAL0uYfieA/1H//9cBPH212+6En5UrAH41hHA9gOcD+Hkzux7N+/cIIYQQQrSfsK4+h2tG298c1odhH63/f8bM9gHYgTX07xFCCCGEaDUhBFRVW7m1mNluLI2iuR1N+jq3V/gAABJVSURBVPeY2U0AbgKAyy5bP12BEOuFclxsZpTfYlOxjuXz1pKO6RyaWS+AvwHwSyGEaVsmnF3Jv6duFnkLAOzdu7eJMuErk+3yQmIuaF8hc98aDVjhGABAnqpsQFumwu+Ls97Qujgz7OKpU95UFwBOk8j9nj4vkO6jASkD3X4AQR8N2sglDOxgqiSIXiTBdJEE1bz8+aat1A6Ou3NUnL4rTulo4AGpy7Mkys7mO+a2AND6HO8aIXPoos9ZowEpKRr0k+3xA6oAwHJ+mpHpMg9YCSSM55jNqIHY5Jo9gFmgv1BpkJ+V+FROk0H9RNFrh06QIfXRaX9zP05C+KVp/n6dn+YBKH6bFbpvag0MroHGJtc8P9vjc6CdtDq/76MBQHx/p80/X7f1jLq4MBLnXo4GWHG+8z2RH+pzcffocRf3H/XPfAAYpkEqxQmfJ8U5n4vF6soDVJLgASgFutcLPf4+LQz54+rbHhcA6B33udR3mX+P07/bW/Rld9IAFDa4pgFCAFDs8t+BJ+b8ffoIFYp48EQ84GfdCEDg0XAbkI74FjSzLJY6hh8PIXymPrkp/x4hhBBCiE5gqXxe4z/oOp22D0ixpVeEHwGwL4Twu8tmNfTvEUIIIYToJGrV0PSnU+mEN4cvBPAmAPea2d31af8T5/HvEUIIIYToREItoFba+G8O2945DCH8G4AEy2UACf49a83Df/SDLmZTbCNDay54z6bZQGPTWit6fURq3ptLl2a8vWOxEOs+5ia87iPT5XU1+QLpxkhXlyFz00w2fqmcIr2KsU6MsNTK8wEgTRrCNBswU7sKWR+XEjSGTC7j12GdY5ra+X3XJ5Wx3DwM/+wHXDz7V7/lYrZXiDSIBZ9bAJBizWGXF3ZGGkTSDyLt45BKMNGmdSqVlTWGRdIYzpV8PLUYi8ZPkf6PNYdHprwO7BiZXh8lQ2sAmKdtlhb5meHbFWor21skGVyzxtBII5el58ETf/mTK+5jI/PuV1zj4vd/5WEXsxE0i5C39gyD6Rmm52XGa7hTpOEcGjjo97DNb7Pv2JloHwMn/HN+4bTXTi5Oe8P1ctHnbzM6N6PqBln6XmDz7q4R0k5uHYq2yceWHxt3cWabH2SU2uJNsCv9fnk2uAZijeHDp73GcB9pDN/zPddG21hPOvmNYLO0vXMohBBCCLEZCAGbQnOozqEQQgghRCsIQaOVhRBCCCHEEiEAVWkONz+PfeQNLt71pj9zcXXRa45AvohArEtkL0SeH/kgsmfZTKw5WsgXXJzJUUw6xTRpxHh59kUDgHSGvdNYL0h+YjQ/k411ZLyNXJ60j2Uf1wreh4v1guzXCDTWGL7pWV4Dc6nR+3+828XFv/9DF0f5mZAbRjpEo3xknWIt7bfBHoZI2AfbEpbpupaq7Gvo41mqWjC1ENc/nV7wy0yS5nCK9IOT837+XIKOMfItrKysKbQU3VcJ54JJ0TMiTfczP8cuJd7xkqtdzBrEavB5Uq6RESqAsR6vvesfIc1hzud3rn/ExekRr2POj5+M9tE/4d3aFie9jq485/Ws1QWvQaw2Ybyczvr7LE2+vpFf4yB9bwxtjbc5tMXHI15DWOvz8ysUTwef3ydn4/uykcaQr3Fbkc+hEEIIIYQ4xybxOVTnUAghhBCiBQRotLIQQgghhDhLkM/hJQn7g21/4580XGehvLLGMLEe8wqwJhEAKgukTyGdUnqONGDkldaU5jDfteIyvA3WNeYSahZncl5fFRrUB02R514137gmdIG8Et/2witW3MelTuH73+rixa99wi9QiTVB7POW6vYarZDOrhijkQYRscaQa2ZzXCK/xnn6qafIxZkBLJAecIY0iFwvnOMkPWEgTRv7f3LM3p6plNd9VUrx84JrwrNfq3gS1qe964sPuHg24Yv92lGvKdxCtekHB7xuOd/t/QBTpLPLjXlPQwAIs97fNjcz6eIae+wukqdmEzW4wX6Y/EwnP8xU36Bfvncg2mSty0+r9Hi95WLW5+/kgm/nSdLxPjoZ6/a5XjZ7WXYSsrIRQgghhBBPEjq7LF6zqHMohBBCCNECQoB+VhZCCCGEEHWCBqQIAEc+8dMNlxl+5XtczJpDJvZBbKb2qteB8TY4Zr0g+zWmyacOANLk4cheao2Oq5KKt2kkEcxk/U1Va6BBZP7wh25c1fKiMfkXvbHhMtUDd/gJpBkMrCnMdtF80iAmaA6rpBFkDSL71bF+sFxlP7v4PuJp7JHJMZNUT5x1slw/nLW4fA8w3/7A61acL1ZHM3V4P3qHr5W8a8Dn73if13wO0POzv59q3fePRftIDXhdnZW8t1+m4n0NrUra07Cyf+bSSqTx5vuSakbXct7zsdbltcQAsAh/705zDfNp387D014r+cgZ/73ycy/YHe3j9U+PJnUsIYSoXvpGRJ1DIYQQQogWwX+obkTUORRCCCGEaAEBQGmVv3h1IuocCiGEEEK0gABgE0gO1TlcD87847tcPPp973dxmTwKK0UfV0nrl6TtS9Yhnn+dRt6KSV6KjZapVUi/0oTvViPPt0iPRR5wf//TL2i4D7H2pK94jotLJx5zsdWo7mvGa7aQ9tc56eHKk8jGMIprYWVNYjNwHe5GcZLmkGuM5wr+WCOfw6xf/hs3v7y5xoo146eec5mL33fbQy4e7/f5zJrE0W7/bOzLxX6s3VnvKVjoH3ZxlvOEU+0CNIeN6pVzffK5YryPubLXux+d8ZrCE3P+u+bhk96vsRnN50YiBL05FEIIIYQQdQKCNIdCCCGEEGKJJc1hu1tx8ahzKIQQQgjRAkLQaGUhhBBCCFFHo5XFBXPqc+9Ycf7Iq9/n4ioNHmHDaiDJONsPBuHBIw3jdDwgJU3G2WykbSky+yWTbBbiA0BXj99GLu/3+29vf2m0juh8clt3rzi/WPQ5bCSmT3q08vO2QhM4vpARg2nzKv88DSbp7fI5PDnvxfjZdDzQADSoqq/X5/gXfva7VttM0WZuftlTVpzPA1aGaEDKcIFM3wGMdGdXXCbHZuqUaoZ4MBQT4AeGsVfzQsV/b8yT8fxE0ec7AJyY8+bcB8nU+n2vuq5huzYT6hwKIYQQQohzLP2s3O5WXDzqHAohhBBCtIAAoCzNoRBCCCGEWEJWNmKNOP2Fmy96G2y03QjWHLKeEADSeW/smiFNIa/z6J/+8KraIC4dCoVC44WIXor3nzntYtYLNiKVsHw27aexxnCAdGC//f1PW9U+xaVBI01iM/zJ7Y+tOD9LGu8E//UIlsItkOhwYt5r11txHJcaNWwOzWGCerpzMLNXmtmDZrbfzFYexSGEEEII0U7qmsNmPxeDmf2wmd1nZjUz27vCcqvuS3Vs59DM0gD+CMCrAFwP4I1mdn17WyWEEEIIkczZ0crNfi6S7wB4HYCvnW+BC+1LdfLPys8FsD+EcAAAzOyTAF4L4P62tkoIIYQQIoEyajhaW2i8YAsIIewDAFtZUnNBfalO7hzuAPDEsvgQgOfxQmZ2E4Cb6uGimX1nHdq2EqMATrW5DUAHtMM++sa2twEdcB7qXHOhKyrHO7cNv9D+NgCdcS6U362n7e34jQ5oAzqjDUCTOT6ByuSncezoarZrZncui28JIdyyuqatSFN9KaaTO4dNUT+JtwCAmd0ZQjjv7+7rQSe0oVPaoTb4dlzouspxtaHT26H83pztUBt8O5pZLoQw1OL9/jOAbQmzbg4hfLaV+1pOJ3cODwPYtSzeWZ8mhBBCCLHpCSG8/CI3cUF9qY4dkALgDgBXm9keM8sBeAOAW9vcJiGEEEKIjcIF9aU6tnMYQqgAeCuALwLYB+CvQwj3NVitlb/TXyid0AagM9qhNjxJq9rRCcejNizRCW0AOqMdyu/W0wntUBuepFPacQ4z+0EzOwTgBQA+b2ZfrE/fbmZfAC64LwULm8DJWwghhBBCtIaOfXMohBBCCCHWH3UOhRBCCCHEOTZN57AdpfbMbJeZfcXM7q+XsHlbffqwmX3JzB6u/9vSoe3naUvazL5lZp+rx3vM7Pb6+fhUXYi6lvsfNLNPm9kDZrbPzF7QpvPwy/Vr8R0z+4SZda31uTCzj5rZieX+bOc7dlviD+pt+baZPavJfSi/25jf9X22Pcfbkd/1/SrH174teoZj8z7DNxqbonNo7Su1VwHwqyGE6wE8H8DP1/f7DgC3hRCuBnBbPV5r3oYlselZPgDg90IIVwGYAPCWNd7/BwH8YwjhWgA31tuyrufBzHYA+EUAe0MITwOQxtLIrLU+F38O4JU07XzH/ioAV9c/NwH4UKONK78BtD+/gTbneBvzG1COXwo5rme4p2X5vSEJIWz4D5ZG6nxxWfxOAO9sQzs+C+AVAB4EMF6fNg7gwTXe704sJe9LAXwOgGHJUT6TdH7WYP8DAB5FfYDTsunrfR7OOsEPY8nD83MAvmc9zgWA3QC+0+jYAfwJgDcmLbfCtpXfbczv+j7anuPtzO/6tpXja7dfPcPD5n2Gb8TPpnhziOTyMDvWswFmthvAMwHcDmAshHC2fM4xAGNrvPvfB/DrAGr1eATAZFgawg6s/fnYA+AkgD+r/yzyYTPrwTqfhxDCYQC/A+AggKMApgDchfU9F2c537FfSK4qv9ub30AH5HiH5TegHG8l7c7xtuc30HE53sr83nBsls5hWzGzXgB/A+CXQgjTy+eFpT8t1swvyMy+D8CJEMJda7WPJsgAeBaAD4UQnglgDvTzw1qfBwCoa0Jei6UH3XYAPYh/Klh31uPY1xLlN4AOyPFOzW9AOX6R++6EHG97fgOdm+MbPb8vhM3SOWxbqT0zy2LpofLxEMJn6pOPm9l4ff44gBNr2IQXAniNmT0G4JNY+lnigwAGzexsecS1Ph+HABwKIdxejz+NpQfNep4HAHg5gEdDCCdDCGUAn8HS+VnPc3GW8x37heSq8ru9+Q10Ro53Un4DyvFW0Qk53gn5DXRWjrcyvzccm6Vz2JZSe2ZmAD4CYF8I4XeXzboVwJvr/38zlnQsa0II4Z0hhJ0hhN1YOu4vhxB+DMBXALx+ndpwDMATZnZNfdLLANyPdTwPdQ4CeL6Zddevzdl2rNu5WMb5jv1WAD9eH/H2fABTy366OB/K7zbmd70dnZDjnZTfgHK8JXRCjndIfgOdleOtzO+NR7tFj636AHg1gIcAPALg5nXa53dj6VXztwHcXf+8Gkt6kdsAPAzgnwEMr1N7Xgzgc/X/XwHgPwHsB/D/Aciv8b6fAeDO+rn4OwBD7TgPAH4LwAMAvgPgLwHk1/pcAPgElvQxZSz9Bf6W8x07loTmf1TP03uxNCpP+d3h+d0pOd6O/FaOXxo53gn53a4cX4/83mgflc8TQgghhBDn2Cw/KwshhBBCiBagzqEQQgghhDiHOodCCCGEEOIc6hwKIYQQQohzqHMohBBCCCHOoc6hEEIIIYQ4hzqHQgghhBDiHOocXuKYWcHM/sXM0i3YVs7MvraszJEQbUc5LjYzym+xFqhzKH4KwGdCCNWL3VAIoYQlR/n/ftGtEqJ1KMfFZkb5LVqOOoebGDP7ipm9ov7/95rZ/5Ow2I+hXjPSzHab2T4z+1Mzu8/M/qn+V+luM3vAzP7czB4ys4+b2cvN7N/N7GEze+6y7f1dfZtCrDnKcbGZUX6LdqHyeZsYM3sRgPcA+FMAPwrgNcv/uqwXuD8YQthWj3djqXbl3hDC3Wb211gqMv5v9enPBHAfgDsA3IOl+pOvAfCTIYQfqG8jDeBYCGHLehyjuLRRjovNjPJbtAu9OdzEhBC+hqUi4b8C4A0JPzuMApikaY+GEO6u//8uALuXTb83hFDD0sPltrD0l8W9y5ZBfR8lM+tr5bEIkYRyXGxmlN+iXahzuIkxs6cDGAdQCiHMJCxSBNBF0xaX/b8KIJMwvbYsri1b5ix5AAsX0mYhVoNyXGxmlN+iXahzuEkxs3EAHwfwWgCzZvZKXiaEMAEgbWb8cLmY/Y4AOBVCKLdqm0IkoRwXmxnlt2gn6hxuQsysG8BnAPxqCGEfgP8N4N3nWfyfAHx3C3f/EgCfb+H2hIhQjovNjPJbtBsNSLnEMbNnAfjlEMKbWrS9zwB4RwjhoVZsT4iLRTkuNjPKb7EW6M3hJU4I4ZsAvtIqA1UAf6eHiugklONiM6P8FmuB3hwKIYQQQohz6M2hEEIIIYQ4hzqHQgghhBDiHOocCiGEEEKIc6hzKIQQQgghzqHOoRBCCCGEOIc6h0IIIYQQ4hz/P6WCDftUqHGdAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plot(sim, component='all', z=10.0, cmap='RdBu')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "and we plot the xy spin angle through the middle of the sample:" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARMAAADmCAYAAAD2gWGmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnX+4HFWZ5z/fe0MIvyGAIYTERGWFiD8SrwEXdEDQiYyPwR+LIANRcbPzKIjKroKzO6yuuwvjLIoro0ZAYBZBBZQsZhEmgoAzYhIEgSRIgACJCSGaQPghcHPf/aPqatfpureququ7q7vfz/P0c/ucOlXnVN3ut8/5nve8R2aG4zhOswx0ugGO4/QGbkwcxykFNyaO45SCGxPHcUrBjYnjOKXgxsRxnFJwY+I4Tim4MXEcpxTcmDilImmdpOM63Q6n/bgxcRynFNyYOI5TCm5MugxJfy/pxzXpr0haJmlPSTskTa05dpikjZL2SLnOOZIelrRd0ipJ76s5tk7Sf5T0G0lPS/q+pEk1x+dK+nV87g/j418eo70HSrpO0lOSHpX0qfKehlMl3Jh0HxcAx0iaI+lvgPnA+83sGWANMLem7PnA/zCz7SnXeRh4G7AX8EXg/9QaIuDE+NqzgDcAHwGQNBH4EXA5MBm4GngfKUgaAP4vcC8wDTgW+LSkvyx8107lcWPSZZjZ74GvAlcA5wLHm9nT8eHlxMZE0tuB2cC3x7jOD83sd2Y2YmbfBx4C5tUU+Xp8/A9EBuFNcf4RwIT4+Mtmdj3wqzGa+xZgfzP7kpm9ZGaPAN8BTmro5p1KM6HTDXAa4tfAecApZvZETf5yYHQm5e+B/2JmL6VdQNJpwGeBmXHW7sB+NUU21bx/Hjgwfn8gsMGSsStq21DLK4EDJW2ryRsE7hijvNPFeM+ky5D0euCbRD2TjwWHlwNzJX0AmAR8b4xrvJKoh3AGsK+Z7Q3cDyhHEzYC0yTVlp0+RtkngEfNbO+a1x5mdnyOepwuw41JFyFpGtGQ42+ATwCvl3R0TZF7gQOA/wWca2NHvtoNMOCp+LofBQ7L2Yx/BXYAZ0iaIGkByeFRLb8Ctkv6vKRdJA3GovBbctbldBFuTLoESXsCS4ELzWyJmT0PfAX476NlzOxF4D5gnZn9v7GuZWariAzOvwJPAq8HfpGnHfGw6f3A6cA24K+BG4EXU8ruAN5DpLc8CmwBLiESfZ0eQx62sXeIZ1rWAiea2S/bWO9dwLfM7LvtqtOpHt4z6S3OA37RakMi6S8kHRAPcxYSTR3f1Mo6nerTNmMi6TJJmyXdX5M3WdItkh6K/+4T50vS1yWtjR2n5o59ZSd2InsaeDtwZhuqfC2RPrMNOBv4oJltbEO9ToVpZ8/kciInqFrOAZaZ2cHAsjgN8G7g4Pi1iGj2whkDM7vbzPYys7eZ2ZY21LfYzKaY2e5m9gYz+0mr63TyI2m6pFtjz+YHJJ01Rrn5kh6Mf7TPSStThLYZEzO7HfhDkL2AaIqT+O8JNflXWsQvgb0D70zHccZmGDjbzGYTORl+UtLs2gKSBoGLiX64ZwMnh2WK0mmntSk13eNNwJT4/TSSjlDr4zzvSjs9iaQ35y1rZiszjm8k/q6Y2XZJq4m+P6tqis0D1sZeyUi6Blgg6UXgl8BzRMPYGcBWYE68ZGNMOm1M/oSZmaTCU0uSFhENhdhtt93efMghh5Tetm5j5UiQET7VwfCMrBPC4zuCdNLX7c3sPE7r+oeVK1duMbP985TdC604MOdAQdKzwIM1WYvNbPEYZWcCc4C7gkNpP9iHm9nDku4kckG4Q9JtwJlZhgQ6b0yelDTVzDbGw5jNcf4Gkl6VB8V5dcQPcTHA0NCQrVixopXt7Qr0QpAROtTvHqQHnw8yXg7Sz2akJyZSK5g1XvP6BkmP5S37GnZiBQfkuy6PP2hmQznq3x24Dvh0HmNQw+uIPKIBDiVpuMak01PDS4CF8fuFwA01+afFszpHAE/7bIHT2wwQWfk8r2wk7URkSK6KF2OGpP5gS9oFmGRmWyVNB7aMtb4rpG09E0lXA0cD+0laT+QTcT7wA0mnA48RLXuHyNPzeCIHrOeBj7arnY7TGUTYw2v4StG6qUuB1WZ24RjFlgMHS5pFZFhOAj5MJMaujsscWvM+k7YZEzM7eYxDx6aUNeCTrW1RNdHSjALh522flDJhbzn8MQt/Z3YJL/pcxgkhyUk61X3+wkm87Gta3RrGXmeAaMlUKRwJnArcJ+meOO8LZrZU0lLg42b2O0lnAD8lUtEuM7MH4tXko0OcF4gWjh5iZmuyKu20ZuI4DhAZk3J6JmZ2J2OsAK9dsW1mS4lGAbXHr6x5fwfwqrz1ujFxnEowqpl0L25MHKcSlKeZdAo3Jm1G38ooMDlIh5+vPJ+3UI4I001/ZrMqaITkdLP4r+OWtozj3ccg3jNxHKcEvGfiOE4puGbiOE4pDFLi1HBHcGPSYvQ/g4zw8xL2bIvqG6Hne9o1ypA0ClWQp8Ki5ySPi+SqeuOiHHVWGR/mOI5TCj7McRynFLxn4jhOKfjUcN+j04OM8Mcl9BvJohVyRMsppnek52Wlw/VCoV9KuN1x8nyj6pElvWfiOE4plLrQryO4MXGcSuDDHMdxSsGHOY7jlIJPDfcd+lCQkfVjkqUr7lSwfB5BNhx6Z52zS8o1CldatHxWmdAbr7kHI/4ikTZ+ntKmTlJePJNO4cbEcSqB90wcxykF10wcxymFAYxdO92IpnBjMg56V0pm0R+PcIuZLLIW/jUiR7ScovpGnjJZe/WE6WIajHhDIm38pq6F7UWMdPnXsbtb7zg9gjHAMJM63YymcGPiOBVgBHi+43viNYcbE8epACN0YHRaMpUwJpI+A3ycaMfs+4h28JsKXAPsC6wETs27TWHD7fi3QUYZ4nqz69mKaijQmIQxLs06v+TRTLI2/iqaztovOfRDqd8exnikLq9VjFBcXqsaHe9XSZoGfAoYMrPDiBYpnARcAHzVzF4DbAXC9bmO0zMYkXnL86oqHTcmMROAXSRNAHYFNgLvAK6Nj18BnNChtjlOyxkh6pvleVWVjg9zzGyDpH8AHifa2/RmomHNNjMbjoutB6Z1qImO03JcMykBSfsAC4BZwDbgh8D8AucvAhYBzJgxo7nGVOG/GWok4UA6S0NJy+t4QOm0BmRpJEX9TIppJG14SIXoBc2k48YEOA541MyeApB0PdEu7ntLmhD3Tg4CNqSdbGaLgcUAQ0ND1p4mO0659IIxqYJm8jhwhKRdJQk4FlgF3Ap8MC6zELihQ+1znJYzOszpZgG24z0TM7tL0rXA3cAw8GuinsZPgGskfTnOu7RzrXSc1tILPZOOGxMAMzsPOC/IfgSY18p69YogI1wBnuZnklWm6CryLI0kq3weX5hNyaS9VcmMOheN6eNWIvtURiOy9I08ZbI0lCyNpNj1BhgmZGLw9XgppUxZjE4NdzOVMCaO0+/soNxpX0mXAe8BNsf+W2ll1gHb4+qHzWyomTrdmDhOBWhBz+Ry4BvAlRnljjGzLWVU6MbEcSpA2ZqJmd0uaWbR8yS9GvglUUdpGzCDyAN9jpk9M965fWVMJgSxThXcfd28ch49IqtMlobS5Pof+7nqM8OfuCeD9AFBOtBUeDlsVPIE0z8Hx6cE6VBjGalrYvl+JMU0kwn8MZFWih4yiR2J9H4kn/WW+k9Mw+ygkDHZT9KKmvTi2EWiKAbcLMmAb5vZYjN7WNKdwIVmdoek24AzswwJ9JkxcZyqUnCYs6VZfSPmqNgD/RXALZLWmNntwOuA++MyhwIP5rmYGxPHqQCdmBo2sw3x382SfgTMk7QcmGRmWyVNJzJcuexcFZzWHKfvafeqYUm7Sdpj9D3wLqLeyGxgdVzs0Jr3mfR1z2QgGCaPvJBM57LHRffNKeiHsse3kulnvhhoJGnLkUINJI/bx7iENxluzBPeVKix7FN3xX9j/yGRfojjghLhROkfgvT4msgAzwfp5D97kBcT6bR/Yztjxbdgavhq4GgifWU9cJ6ZXSppKVHsoEnAjyKncyYA3zOzmySdxp+HOC8AcyUdYmZrsursa2PiOFWh7KlhMzt5jPzja5JvTDl+Zc37OyAlatQYuDFxnArg7vSO45SCxzOpOAcG8kK4xVH4zxuZMH4aYCTpeoCF+/Q2qZHYGUGj3x0UCKWD0P0irc4wvbVYm7IXIIXpUFOZXHfF32pZkJP0VZF9IDg+vvCT5TcS+ozkWeKUVWZu4HdydxN+J94zcRynFNyYOI5TCj7McRynFEYDSnczfWVMwpF9XRTQ4fHTkKKrBGV27FywUZ8L0uEylyzyxIANdZYsMeC5rAJhOtREstJQf6NhmU8E6b9LpEKNJPQbyQo704hm0mrSNLpUWhdWpSn6ypg4TlWxARjOu9VwRcUVNyaOUwVUoGdSUbq8+Y7TIwyQ341gWysb0jg9bUyyxryNbOObpauE632GQ7+U4Ho75gd+JUV9QPJoJqEcEcY3yWSnjAuGa2/C42lCULLMMHsl0tuVjDO7jyXFpV1K1kg6rpkMUjx+cMXoaWPiOF2DqLfZXYYbE8epAkWGORXFjYnjVAHR/rnoknFj4jhVwDWTcpC0N3AJcBiRRvkxoriT3wdmAuuAE82skDyZ5ipVS5b7ciN7goeC7LOBT8AD25OC68uBQDsY/EcmNLLpVtio8KmFemndNZt1Sgs38QojWMNL7JFIPxfUGbpS3Kakk9oplnyOrRBgQ1oqafRAz6QqYRsvAm4ys0OIArasBs4BlpnZwcCyOO04vcmoZpLnVVE6bkwk7QW8nXgvYTN7ycy2AQuAK+JiVwAndKaFjtMGBoh6JnleFaUKw5xZwFPAdyW9EVgJnAVMMbONcZlNFF+14jjdwwD1YWC6jMLGJI5k/Ucz25FZOH8b5hJt9HOXpIsIhjRmZvFGQWntWQQsAtiVGZxSM5QOrU+m3hGk0/63YSyiPI5u4x3fPJhM7x5EcNo9+A/tVIaGEgacLqyZhJrIqxOpkeD4Syl986eDTnGxcNH1vf2wheEdhHpHI/urZaXPJmVDtLz0wNRw5jBH0oCkD0v6iaTNwBpgo6RVkr4i6TVNtmE9sN7M7orT1xIZlyclTY3bMBXYnHZyvAvZkJkNTWL/JpviOB1iVIDt4mFOHs3kVqKfnnOBA8xsupm9AjiKaE/SCyT9daMNMLNNwBOSXhtnHQusApYAC+O8hcANjdbhOJVndGq4iwXYPMOc48ysLtKomf0BuA64TlKzs2ZnAldJmgg8AnyUyND9QNLpwGPAiU3W4TjVpQemhjONSZohaaRMxvn3AGl7px5b6Dokx9aZwZAyrpcnjF7WNUN5IlxjF2oDdf4SwQds8r5B+aS7BpCiq4SEAsW9QfqJUIF4dZA+NJEaZlYi/fsgdHdYXVpe+NzC5xJGIQvLZ7nOFF30madM1kLRQvSAZpJbgJU0BPwt8Mr4PBFpo29oUdscp38YnRruYorM5lwF/CfgPqKQlY7jlEWfTQ0/ZWZLWtYSx+ln+mmYA5wn6RIi1/Y/RaYxs+tLb1WDhP+PcAxbtBfZiGYSsjBYQ5LlTxEuownbHGowk1NuavIrk+n9A0GhbgpveZC+N/QjOTKRepbXJdIPB6XDe0xbUFXU5ydMnxS4HT1p4/t4lKGZZNHUNfpsmPNR4BAi/5/RYY4BlTEmjtO19FnP5C1m9trsYo7jFKYfpoZr+BdJs81sVcta4zj9Sp/FMzkCuEfSo0SaSeWmhrM0k06Q5XeSRSNj/bqVNXsm09PnJ9NJrxHgxk8nkrcGGskTQfEsn5EytKes8vURU9pPU56bfdYzmZ9dxHGchih5aljSZcB7gM1mdtgYZeYTxRIaBC4xs/ObqbOIMdkEfIAo8lnteV9qpgGO49CKnsnlwDeAK1OrkwaBi4F3Ei22XS5pSTMyRhFjcgPwNFG8kRczyjqOU4QBYJeyonqAmd0uaeY4ReYBa83sEQBJ1wALJL1ItID3OaLtvmYQze7PMbNnxquziDE5yMwqPdQJtx6pgkNh6GPRrDaQ53jWvuWh5nHX/DcnM05IaiSv/HHycOiFkhVLpJF1L0XL1+k+HaA5jW4HBTYR3k/Sipr0YjNbXLDCaSQ/CuuBw83sYUl3Ahea2R2SbiOKNTSuIYHiszmvN7P7CjXZcZwchMtUx2WLmaUtjC2L1wH3x+8PJQrunkkRY3IU8JEqz+Y4TvdSqGdSBhtIdjIPAjZI2gWYZGZbJU0nMly5rFwRY/LuAmUdxylEoZ5JGSwHDpY0i8iwnAR8GJhNtDsERL2S1emn15NpTCTJIh7LKpO30lYR+v00K46XMdYP45eEfibNrlFJW/cSaiTrSAaafZk3BSWSfiV8MpkM//GP/SKZ3uOpZDprVx3IjsmalQ79SrL8TMqIXZNFc9cYoT5qS+NIuho4mkhfWQ+cZ2aXSloKfNzMfifpDOCnRF+dy8zsAUmn8echzgvAXEmHmNmarDrz9ExulXQdcIOZPV7T2IlEQ5+FRKEdL897o47jhJTbMzGzk8fIP77m/VJgaXD8ypr3dwCvyltnHmMyn2iHvavjLtE2YBKRNbsZ+JqZ/TpvhY7jpNF2zaR08oRt/CPwj8A/xrFe9wNeiDfKchynFEboeWNSSxzrdWNmwQ6RtYq7E3FBjwnibvxdEHcjHCVn7RcT6iFp8VW3s1eQE+5GcmSQnj5uknkpldTWF0h0z25IpjelfEeyfFPC/2PoM3R58BxDlbAT8X6b0+jaLsCWThV29HMcpx+GOY7jtIPu75nk3rhc0i3xXsCO45TO6NRwnlc1KdIz+TzwNUnrgC/UbCpeCvEqxhXABjN7TzxzdA2wL9HiwlOzPPFC2150zUgjGkmWv0TY4Kx9c8J0qIm8wM6J9DBBsBIgWptVS7hyJdwHZ/K4ybriaUJNDRY8hBdTym96OpmeGNx4qJlk7YEU7puTpW/k2eipqM9PxmPJoI96JmZ2t5kdA9wI3CTpvNj1tizOIqmjXQB81cxeQ+SbdXqJdTlOxRjVTPK8qkluYwKRpyvRop9vEm3p+ZCkU5tthKSDgL8CLqmp5x1Em5gDXAGc0Gw9jlNdRnsmeV7VpMiOfr8AZgEPEMU7+AiwBjhL0tvMbFET7fga8DlgdLPLfYFtZjYcp9cTLZl2nB6lv/xMFgGrUtbgnCkp92KgEEmjoeVWSjq6gfMXxW1jxowZ/OeahSTXj7+VSmbMzjzxVbPSWWP9cK1O+HF6OfCwGGHvoETaqpRQ5CgYgaQuiGzG5bJ+LFMeZJ2uklw+xHAQfmviH5PprNjLza55gnrflqzbfD/Jr8a/J+MDmKCPpobN7IFxDv9VE204EnivpOOJ3PT3JIpLubekCXHv5CCilY1p7VoMLAYYGhrq+GJDx2kMI58sXF0KaSZjMRr6rcFzzzWzg8xsJtEy6J+Z2SlEiwc/GBdbSBQ20nF6lNFhTp8IsG3m88BnJa0l0lAu7XB7HKeFjNA3Amw7MLPbgNvi94+QuSrEcXqF/hJgu47QkSmLPE5szQqw/6JkBPI9AyXy5bozskINhWponnMKCrDhgwxPbyRqdsaD2xF4gP1ufVLMfDmoY2JwfhkCbHv7AN3vtNbTxsRxuoc+ms1xHKeV+DDHcZxS6P6p4Z42Js0GQ2pEMwkVj10YSaQnkPS+eklnJ9KyHwZXCAWKKRnHIVvkKLjFdpYTWxk7lwfY4YHDV/CjHd5BuEpsQnI9JBMDp7hGNJPWKhojDPB8zpLVpKeNieN0D8YAw9nFcGPiOM44CKvrtY5FPpPTftyYOE4FEPljyOYzOe2nr4xJUR+RPJpJuBhsQoZGMoEXgjM+EaSXBeksjSRNM2kyLFTRBxNWV4bYED6WjImO8IM8EPx8D0wKygcaSqipQPZtlDn3Em4gNx6ZO4h3iL4yJo5TVQao/2HqNtyYOE4FKDLMqSpuTBynAhQZ5lSVnjYmRwbRTZ7IiFWTRzMJ/UbC6bxQIxmoc0QKR9rJhS2mryXS93JUIv0mOys4P61znBWOueBvYJYzTVY6xRdr01eS/4wp1wYFigZhDwSNgaANA0EbhgO/lIGUb0Koq+waHJ8yWF74HO+ZOI5TClm7UXYDbkwcpwJ4z8RxnFJwzaTLKeozEuWFmshwkM7aijzreDL9RpJigmlOIi1Lc67OurMmNZPwcuG3INgDzLaliFXfyLhGFqHTR9b5wT2EH/yRtOVKgW/KSIovSll4z8RxnFIo289E0nyiwOyDwCVmdn5KmXXAdqJgKsNmNtRMnW5MHKcClDnMibfavRh4J9GeU8slLTGzVSnFjzGzLWXU68bEcSpAycOcecDa0V0jJF0DLADSjEmyHdKriTbZew7YRjRo3QrMMbNxPfn7yphMD9wCng2G8hODgXjakvD6tTVZKzZCjSSrfFY69EtJ2y47PCf5MZV9atzjDAb3fUDyY2JfCB5cln6RtnyoKFm+LlkLZTK+qQM5vsk77dm6bZkKTg3vJ2lFTXpxvH/UKNOAJ2rS64HDU65jwM2SDPi2mS02s4cl3QlcaGZ3SLoNODPLkECfGRPHqSoDFOqZbGlW34g5ysw2SHoFcIukNWZ2O/A64P64zKFE+4tn4sbEcSpAyU5rG0huW5C6I6aZbYj/bpb0I2CepOXAJDPbKmk6keHKte67yptwOU7fMKqZ5HnlYDlwsKRZkiYS7ZS5JFGftJukPUbfA+8i6o3MBkb3Dj+05n0mHe+ZxNbvSqLAHUY0/rtI0mTg+8BMYB1wopmFO7Q0xe7hEFjbg4xGIoMW1Uia1VDSxIJQRwl1lvcFx4Pt03d6IpkO90YP01kbFOUJDJN1PEsTydJQwuuFjz2lPQNT27d19SDlTQ2b2bCkM4Cfxpe+bHSvcElLgY8TedH8SBJEduB7ZnaTpNP48xDnBWCupEPMbE1WvR03JkRR6M42s7tjS7lS0i3AR4BlZna+pHOAc4i2DHWcnqPstTlmthRYmpJ/fE3yjSnHr6x5fwfwqrx1dtyYmNlGYGP8fruk1URq9ALg6LjYFUTbhroxcXoS94AtGUkzgTnAXcCU2NBA1A8P4xc6Ts/ga3NKRNLuwHXAp83smXgsB4CZWTwXnnbeImARwIwZM9KKFKCR3VTK9iMpev0O7PhSxka+4W2Ha2OaDbDaZd/MXuiZVGI2R9JORIbkKjO7Ps5+UtLU+PhUYHPaubGjzZCZDe2///7tabDjlMyoZpLnVVU6bkwUdUEuBVab2YU1h5YAC+P3C4Eb2t02x2kXIuqc5XlVlSoMc44ETgXuk3RPnPcF4HzgB5JOBx4DTuxQ+xyn5bhmUgJmdieRYU7j2Ha2BZuaTOuhlEJFNYyi6Wav3+g545B1ehjjtQzpKcsvpOimNlnf1MPb51OSRi9oJh03Jo7jgAwmtnZn9JbjxsRxKoBGYEJV9/3MiRsTx6kCBkrZEqSbcGMyHnZwfZ5+FWQUHbwX9TNpRP9o86eyFbJO0X1zsgSH8DG/t7MaSR0jlLt5cQdwY+I4VcBo+29A2bgxcZwqsAPvmTiOUwJG6ase2o0bE8epAq6Z9CE2L5nWz4ICRdXIoh5fjTigtdhprZHyRW8r9CMv2qYPVkxwDfFhjuM4peDDHMdxSsGHOY7jlIJPDTvYO8Y/rsVBRtkLBdNoc385rC6MjNyGtYh8tuKaSBaumTiOUwqumTiOUwqumTiOUwojeM/EycAWJdO6IOOEZp06GrlGQYq6zjRyzazjF3S5RhLiPRPHcUphhOIrpSuGGxPHqQIuwDqOUwo+NewUxjJ2ONUnS6ikw5pJGde4osc0kSy8Z+I4Til4z8RxnFLwnonjOKXgU8OtRdJ84CKiDc8uMbPzO9yk1mMXN38NzQ4yCu5olRViJbzcH4L02j7TO8qg5KnhrO+OpJ2BK4E3A78HPmRm65qps+N7DY+FpEHgYuDdwGzgZKnuW+I4vcGoB2yeVwY5vzunA1vN7DXAV4Esb8pMKmtMgHnAWjN7xMxeAq4BFnS4TY7TGkaHOXle2eT57iwArojfXwscK2msbXpzUeVhzjTgiZr0euDwDrXFcVrK8yOwsrwd/fJ8d/5UxsyGJT0N7AtsabTSKhuTXEhaBIwugHlR0v2dbA+wH038QyrZhnAsnzW2l3rvGTTGa/MWXA3bhobZmPe6klbUpBebWRg4p+1U2ZhsAKbXpA+K8xLED3ExgKQVZjbUnual423ofP1VakPesma2T4lV5/nujJZZL2kCsBeRENswVdZMlgMHS5olaSJwErCkw21ynG4gz3dnCbAwfv9B4Gdm1tQ0XGV7JvE47gzgp0TTW5eZ2QMdbpbjVJ6xvjuSvgSsMLMlwKXAP0laSzS5f1Kz9VbWmACY2VJgaYFTOj5uxNtQhfqhz9uQ9t0xs7+ref9H4N+VWaea7Nk4juMA1dZMHMfpInrGmEiaL+lBSWslndOG+qZLulXSKkkPSDorzp8s6RZJD8V/y1Tpx2rLoKRfS7oxTs+SdFf8LL4fi3CtrH9vSddKWiNptaS3tvM5SPpM/D+4X9LVkia14xlIukzS5lp3hLHuWxFfj9vzG0lzy25Pp+kJY9Ih1/th4Gwzmw0cAXwyrvMcYJmZHQwsi9Ot5ixgdU36AuCrsav0ViLX6VZyEXCTmR0CvDFuS1ueg6RpwKeAITM7jEhwPIn2PIPLgflB3lj3/W7g4Pi1CPhmC9rTWcys61/AW4Gf1qTPBc5tcxtuAN4JPAhMjfOmAg+2uN6DiD607wBuBETkrDUh7dm0oP69gEeJ9bea/LY8B/7syTmZaELhRuAv2/UMgJnA/Vn3DXwbODmtXK+8eqJnQrr78LR2VS5pJjAHuAuYYmajnoybgCktrv5rwOeIVndA5BK9zcyG43Srn8Us4Cngu/FQ6xJJu9Gm52BmG4B/AB4HNgJPAytp7zOoZaz77uhntB30ijHpGJJ2B64DPm1mz9Qes+gnqGXTZZLeA2w2s5WtqiMHE4C5wDfNbA6Rs31iSNPK5xBrEguIjNqBRJuThkOPjtDq/3/V6BVjksv1vmwk7URkSK4ys+vj7CclTY2PTwUpBNxrAAADA0lEQVQ2t7AJRwLvlbSOaGXoO4j0i71jF2lo/bNYD6w3s7vi9LVExqVdz+E44FEze8rMXgauJ3ou7XwGtYx13x35jLaTXjEmbXe9j5drXwqsNrMLaw7VuikvJNJSWoKZnWtmB5nZTKJ7/pmZnQLcSuQi3Y42bAKekDS6qO1YYBXtew6PA0dI2jX+n4zW37ZnEDDWfS8BTotndY4Anq4ZDvUGnRZtynoBxwO/BR4G/rYN9R1F1IX9DXBP/DqeSLNYBjwE/DMwuU33fzRwY/z+VcCvgLXAD4GdW1z3m4AV8bP4MbBPO58D8EVgDXA/8E/Azu14BsDVRDrNy0Q9tNPHum8iYfzi+PN5H9HsU8s/F+18uQes4zil0CvDHMdxOowbE8dxSsGNieM4peDGxHGcUnBj4jhOKbgxcRynFNyYOI5TCm5M+hhJu0j6eRzCodlrTZR0e40Lu9NnuDHpbz4GXG9mO5q9kEU7xy0DPtR0q5yuxI1JjxJHgXtn/P7Lkv53SrFTiNeOSJoZR0n7Thy17Oa45zIzjqB2uaTfSrpK0nGSfhFHE5tXc70fx9d0+hB3p+9RJL0d+BLwHeDDwHtreyDxgsjHzeyAOD2TaB3LkJndI+kHRIvT7ozz5wAPEC2qvJdoHcp7gY+a2QnxNQaBTWa2fzvu0akW3jPpUczsdqLFZZ8FTkoZyuwHbAvyHjWze+L3K4miiI3m32dmI0QGZZlFv0L31ZQhruMlSXuUeS9Od+DGpEeR9HqisIEvmdn2lCIvAJOCvBdr3u/gz/sq1eaP1KRHqN97aWegvC24na7BjUkPEgfluYooAtmzkuoij5nZVmBQUmhQmql3X2CLRUGKnD7DjUmPIWlXomhjZ5vZauC/AeeNUfxmorgsZXEM8JMSr+d0ES7A9jHx3i2fMbNTS7re9cA5ZvbbMq7ndBfeM+ljzOxu4NaynNaAH7sh6V+8Z+I4Til4z8RxnFJwY+I4Tim4MXEcpxTcmDiOUwpuTBzHKQU3Jo7jlIIbE8dxSuH/A8ce/vP2Oz/tAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plot(sim, component='angle', z=5.0, cmap='hsv')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Simulation using Steepest Descent" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "An alternative method for the minimisation of the energy is using a SteepestDescent method:" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [], "source": [ "# Create simulation object.\n", "sim = Sim(mesh, driver='steepest_descent')\n", "sim.Ms = Ms_function(Ms)\n", "sim.driver.gamma = gamma\n", "\n", "# Add energies.\n", "sim.add(UniformExchange(A=A))\n", "sim.add(DMI(D=D))\n", "sim.add(Demag())\n", "\n", "# The maximum timestep:\n", "sim.driver.tmax = 1\n", "\n", "# Initialise the system.\n", "sim.set_m(init_m)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In this case the driver has a `minimise` method" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "#0 max_tau=0.01 max_dm=0.0477 \n", "#1000 max_tau=0.01 max_dm=0.0022 \n", "#2000 max_tau=0.01 max_dm=0.000592 \n", "#3000 max_tau=0.01 max_dm=0.000327 \n", "#4000 max_tau=0.01 max_dm=0.000179 \n", "#5000 max_tau=0.01 max_dm=9.82e-05 \n", "#6000 max_tau=0.01 max_dm=5.43e-05 \n", "#6141 max_tau=0.01 max_dm=5e-05 \n" ] } ], "source": [ "start = time.time()\n", "sim.driver.minimise(max_steps=10000, stopping_dm=0.5e-4, initial_t_step=1e-2)\n", "end = time.time()" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Timing: 67.3641722202301\n" ] } ], "source": [ "#NBVAL_IGNORE_OUTPUT\n", "print('Timing: ', end - start)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And the final state is equivalent to the one found with the LLG technique" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAocAAADlCAYAAAAht12aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztvXmYJGd15vue3Coza196rV4lWhJCIJa2EIZhMMtlM+AFe8C+GGP5ivEYGy8XIw0zYMv2WFz82JaNzUUGbOzhATyYMWKxZVZj4CKQQKCl1S2ppd7Ua3UtXWtu3/0js1t13ojOzOrOqsyqfn/PU490IiIjvog48eXXke/3HgshQAghhBBCCABItLsBQgghhBCic9DgUAghhBBCnEODQyGEEEIIcQ4NDoUQQgghxDk0OBRCCCGEEOfQ4FAIIYQQQpxDg0MhhBBCCHEODQ4vYczsBjO708w+YGbjZrbPzK42s7eb2UEzO2VmP9XudgohhIiiPlwsFxocXtpcC2A3gE8BGAFwH4B/rq27HMDvA/hv7WmaEBeHvjjFJYD6cLEsaHB4aXMtgD8KIXw5hFAG8CCA+0IIt4UQigDuB5BqawuFuHD0xSnWOurDxbKgweGlzTMAfG5RfHVM/NCKtkiI1qEvTrHWUR8ulgUNDi9RzGw7gAyAfYsWPxPAvYviZwC418wyZvYdM8ua2U+a2QdWsq1CXCD64hRrliX24QNm9rXa3yEze9tKtlWsPjQ4vHS5FtW3KBUAMLM+ANsB/JC2+UEIoQDgwwD+HMCbAPzaCrdViCWxxC/OV5nZ79Y+N2Rmn1+xhgpx4SylD58IIbwIwM0APg/gr1a4rWKVocHhpcu1iH5RPhJCmAUAM0sAuGbRNt8CcAOA/xJCKK1kQ4W4AJr+4gRwN4Bn1ZbdDODWFWynEBfKkvpwM7sOwC8B+NWzz4UQ58NCCO1ug+hwzGwEwCcBfBlAIYTwx21ukhB1MbP/DmA0hPCfa/ELAHwohHBVLU4AmAawK4RwxMy+AeCNAP4khPAz7Wq3EMuBmT0HwK8CuFH/uBfNoMGhqIuZdQH43wDeiao+61sAXhRCmGlrw4RoIWb2jwCSAP57COG+drdHiFZiZkcB7AdQBHBvCOE32twk0eFocCiEuOQxs/8K4MoQwpvb3RYhhGg3K6Y5NLOPmNkJM7t/0bIhM/uimT1c++9gbbmZ2Z+b2SNm9kMze/ZKtVMIcUnyDADvbncjhBBiMXFjJ1q/LOOllZyQ8rcAXkHLbgLw5RDCLlT1bDfVlr8SwK7a340AZJ0ihGg5ZrbRzO4A8O8hhAPtbo8QQhB/i+jYaTHLMl5aMQPYEMLXzWwHLX4dgBfV/v+jAL6GqrbtdQD+LlR/8/52zaNpUwjh6Mq0VghxKRBCOAbgte1uhxBibVCb/NM0IYR7GqyPGzstZlnGS+2uDrBh0QkcA7Ch9v+jAA4t2u5wbVnkZM3sRlRHy+ju7n7OVVddtXytXUWc+H7sG+imMYtZRnGCNuLX0APXPu2i2rCWuOeee06FENZdyGeV4/FML/hJl5yzlcCxX1DmDQCUaRFrso0Osq4700RL1z7K79ZzYnqh7nrOxZguOwJnfKNnYnNftom9Xho0neOJ9N1I55ver5lNA9i7aNHtIYTbl9C0psdLS6Hdg8NzhBCCmS15dkztIt4OALt37w533313y9u2Grmt74olbZ+kjiaTiHY1aVqUS/rhIH/mp3QvzmFmF/yTpXI8nm8+NuZizuH5srdy48Hk5ELU0WO64JfNl/w+0gmf8//leTuaautaR/ndem775n4Xc37zP87TycbDwyL962e+VHbxxGzRxbe8XAP1szSb45btR/qqn2h6v4XvfWhvCGH3BTdsmWj34PD42defZrYJwIna8iMAti7abkttmRBCCCFEx2KJ5EoeblnGS+2ukHIHgLPWEW8G8JlFy3+hNgvnegCT0hsKIYQQoqOxBBKpTNN/LWBZxksr9ubQzD6O6uSTETM7DOA9qJap+gczuwHAAQA/W9v8CwBeBeARALMA3rJS7ewE/qTX/yQ8Rz+PzbEwCkAhRj9VD/4Fgn8SjvtZmZflkv6Y3bTTz26+xsX5kVxkn7lBr2lJZdMuvu7Or0Q+I1Y/M7NzLk7C53gp5t+t0wW/zUje58ps0a/nn5XLpK9ivRUA0EcwV/A/ux2ennXxGz/6XRefni5E9lmknXZ3+W73s299XuQzYnXz7jsfcvFwjx8EDOV87gLAcN5vs3tzv4uzKf9MpEnWE9NlR+CvCc7N2aLP9zv3Hnfx5HxUinFixuf8gTFfH+F9r/HfA2sdA2DJ1r05PM/YKQ0AIYT/F8s0XlrJ2cpvPM+ql8RsG1At9SOEEEIIsTowQ7I1bwQB1B07nV2/LOOldmsOhRBCCCHWBmYrrTlcFjQ4FEIIIYRoAQZDIhWVDaw2NDhsA2wzM012GZOknWI94VL1hUBjjWGMjLGJffoP8T5n6DwS41HfrkAH7urzn/n+T73cxdn+Lhc/9W/uaK6xYkUpnHjcxVbxWqV0psfFlYzXo1oi+rMM+ximKN84bsLZI0KxQhos0hxOk+aqQM8ua7gAoET7GCON1nW3/KuLU2n/1uFbN0eUN6LNsKZwdMBrp6/bOuDi9aQ57O+KfvXm015DmE3V14EnExfvc8i+hvzdMl/y8Vwpmt8T9EzsGPDP8ke+e9DF+05Mu/jWV1993vauSvTmUAghhBBCnMMMifTqN8fX4FAIIYQQogUYVtzncFnQ4FAIIYQQohVYoqWzlduFBodCCCGEEC3AzFplbt1WNDhcAW7t3uVinoDCMYuCoxLgxrCFMNflbFRLmdfHbcPxhZTbqdCElAoJ+hem/CSWBM0y2Pcrr3fxFR/41AW0Qlws5cfucXEiQV1L0neWITHvYv4ZJtEV7Vx5ggl5AEdirj0bl9ONYMF+ozjETBYrkbFwcYFjmuQy5+Nn3vwFF9/7R6+q02KxHHzwrsddfP32QRdv6vET5QZzPp/7Mj7Op6K5aPNTFHujeCvSM1P2dZARmvimMP+QhKSfVZtL+4k1fWk/uST0+IlkADCS88/6+m6/z6mFHK33z/b7v/WYi9/2ozsjx1ht6GdlIYQQQghRpVY+b7WjwaEQQgghRCuw1pbPaxcaHAohhBBCtADTm0MBAH83/FQXT8WYhLKmcI50deVQ34G6GS1fIz1goziXrL++ug2ZtDbYBxeKT6SjZ5LM0DZcTJ6FZA04+j+iJSa7Br1OZuhX3rukfV7qLHz9435BqRjZJtE/7ON8r98g2+3CiCaHNIlsmg0ASfpMuoFONp1sECdi8pF1ijHPQT3iNIeso2UNIscF0iRWqP/YccMnIsfI5L1WbN9f/ETjxgoAwE2ff9DFV22I6uqu3dDn4g09Xlc3kCVNYSCt9OxJH89NRhsyd8aFlekJH8/N+O3pOQwlb64eh/GghSp5JHL+OU32eDNv5Oi5BtCb63dxd97rMfu7vB6zm74HRvK+TZ+490jkGPcd9XrMP3zlUyPbdBLSHAohhBBCCADV2crJ1IVMz+wsNDgUQgghhGgRdgHOCJ2GBodCCCGEEK3AoDeHlyIfG6mvMZwk/RDQWFMY9SDk9T5uRg/I+r9GesFGGsS4bVJpr6tIZX06sZ6Q1wNAijyyeJsU+YXx9ka6sdK81/oAQGrOa1qm/+fvubjn/3xP5DOXMnOffX/9DeL0NDH6vcUkqdaolfx9CuSvFqc5TCcpv+i5SSVY8+rb2ZWq1I2ry/x55MmfLkdxhra3JjSKrEssFUiTWCBN4rzXmpVjtGXlgvfE2/l//S8XP/bXP9OwXZcKv/fFvS7evc3r6q4Y9ro7ABihfmioi/rsM8dcnJgZc3Fl/ISLixQDQGXSf6Y05XWJxRnvc1ia93lQKUSfGSaR4f7WP5fpbv8cpvq8npC1xQCQHFzvP0NxT7f/TK53g4u706RJzET7lxR9H/E9fM/Lrox8pl1Uy+fpzaEQQgghhAAAMySXOJGyE9HgUAghhBCiRejNoRBCCCGEAACYAcmY8oirDQ0OG/CZjU9z8RzVAo56Fkb3Ea1jXH89y/1YL9gTI3btpg/xNt1572fF2r4M1cNMxug+WEPIej+L0SkudZ9J1jGSxjCiiSHNTJwvotGyQPeMNXa517wtso+1DGsw+fokMj43LOPvAQAEyoVA27BHm1FNVyt7/VScZ1sqRfeebnWG8q9U8XEv5V+xHNUcFsspiv0DPZv3esAz817nNT0f7VK5drIl6mvDQoX6FLoW5YLXngFAqUw6RdIgbn3T37j40N+/pW4b1hKsT3vGJu9ZeOVI3sXr89F72FuZdXHylNcYYvyoC4vHD7q4dMJ7980cOx05xvyY1xguTEy7uDDl72mRanBXYvTuTIL61zT1r5k+Xwe5a8B7PmaHvQYRALo3Drk4tX7Uxxu2+XjB9wWDfRv9+h5/PwAgmSD/RetsDaJmKwshhBBCCADyORRCCCGEEIsxaEKKEEIIIYSoIiubNco/X/ZMF3Nd5CJ5FjbSE1a34dhvxP6BrBfs45j0gQCQG/R6rNxwru767KDXdSRJu5dMLz01ykWvgWHtGmup4mAdY5J8ubidrDlMdfvzjvsMa+gYriXc9cI31t1+tXH6A+90cXnOa9pYo8nXL90do22iGq027zVagWq2RnwMy6TLq0SPEegzGfI9LAfW/3IH7fVWiRhdED/PiQadfJk8CwsxtdXnI7WT/XkszNU/RqBrUSlG9ZgV0iWW5r1erZzxOsXRn/+Qi4987JfrtmE18a5/3uPiZ23xOjnWGG6k/rSn6Ov4AkBy0msGy0cfc3HxyKMunj3kNYlnDh138fRRXzcZAGaOey3e7Jh/hham/D0ukt61EFPXm+HvmjRpz7v6/LOeH/bXqntD1AOyZ5P3iezd6s8tT/6N6VF/fVOb/Hn19XvNIgAkuvtoiW8XewlzvexbX311ZJ/Lhn5WFkIIIYQQi1kLE1I6YnhrZr9pZg+Y2f1m9nEzy5rZTjO7y8weMbNPmlmm8Z6EEEIIIdrDWSubZv86lbYPDs1sFMCvA9gdQrgG1d9+3gDgvQD+NITwFADjAG5oXyuFEEIIIRpjCWv6r1PplJ+VUwByZlZEVUxwFMCLAfxcbf1HAfwugA8sd0PmZoouZt/Cxp6F0X1mSEeXpn30p8nHsMe/JO0hnUfflt7IMXq3eK+p7tF1PiYvqiTVyDTWhMXUz2XtUyDvugrXgJ32uieuDQosvT4o1wZlbSTr44CoN2LEp4/PteRzYK2xMO7vC9ejZq/IOB0nk+ny9zbiU0jX1EgjxBrDEKLaPdYpplP+3gd6+BIRba8/ZlfMv9rz9CwO5rwejb0SeftkE539QaqdXFyoHxcoP/k5BKKaw4g3Yoxv5Frl6exjSLWS2ceQNYbJiUORfZYOeh+9woGHXTzxqNckTj3uNYeTB7yH4eShM5FjnJ7yz+F40T8Dk6xdJY1hnMcuw99PrEHsH/fP8eATvq8YOuDrIANA/1Z/LnNj/nr2Tfj1gwv+GGny7Uxtiz77PQNbXbw+778DS3SPpxca15leLmyNlM9r+xmEEI4A+GMAB1EdFE4CuAfARAjh7B0+DCCqUgVgZjea2d1mdvfJkydXoslCrCjKcbGWUX6LtcZaeHPY9sGhmQ0CeB2AnQA2A+gG8IpmPx9CuD2EsDuEsHvdunWNPyDEKkM5LtYyym+xljADMqlE03+dSif8rPxSAI+FEE4CgJl9GsDzAQyYWar29nALgCN19iGEEEII0VYM1pS8pNPphMHhQQDXm1kewByAlwC4G8BXAbwewCcAvBnAZ9rWQiGEEEKIBpx9c7jaafvgMIRwl5l9CsD3AJQAfB/A7QA+D+ATZvYHtWUfXo7j/9t1z3NxFwnOGxmLphsYXANAjlTA3Xkvcu9e7w09+7Z4QfXgUza4eOAKL84FgNxlvtB4etsVLi73rndxJeuNYUPaT9IIMedhJS8cTsx7oXF6zouuMXnCt2HcxwBQIYPUyqzfZ2U+OonFtYmEv9aVjW5DE1ASFPNknETvoIvLj93j4uTO59RtU6ex71deX3c9G4+zsTgTZ5AeSGQein4CSijHGGf7RtSP4z5CcYaeM45R3/scAEBzALBARu4DWd8/jORp8lhX9Np00RcFG2cfKPKEFDI3nvUTgkopP/ELAIpkes1G2TyJpbww5+KhV9zi4tP/8u7IMTqV93/LG1Lv3uz7tnU0AaW34s2l2eCaJ58AwMJjftnpPQd8vO+ojx8ed/FxMrh+Yj46YeI0TVSKFmCIfKTlHKN2cUGGoZh2b6aJNBtO+dxaoPVlmoQ4VGw8eSRF/UHv4DYXj+T8M3LVSI+LP3jX4y5+63N3NDzmhWIGpPTmsDWEEN4D4D20eD+A69rQHCGEEEKIJWNmLX9zaGavAHAbqlZ/Hwoh3ErrfxHA+/Ck/O79IQRfBmmJdMTgUAghhBBitWMAMqmoFdwF788sCeAvAbwMVeeW75rZHSGEB2nTT4YQ3taq42pwKIQQQgjRAsya8ztdAtcBeCSEsL+6f/sEqg4vPDhsKZfU4PCe17y04TYJMrXtLvs410iDmI3RHFEx8+71XuM2sN0XLh986nYX9z/j6f4YV0Q1b8X1XmN4cM7rVY6d8Rqw46e8DmSuSHrCmNzup3Pb1Ou1PYM93mh73bqnuJi1PQCQnj7l4tIpr92pTFOB+oiWyp+npaPCMkv5629Zr/G0Lq9XSZDBakj48y6ceDxyjMz6HZFl7WDPW14bWVaYrm+AnMqRmTRdUzYeT8VohAJp81jHaGl/D0LS36dKxt+DkIqa7YYU6WJpvZFxtlFOW9FroazgtWcAkKJ95Ejr1J/1euB1w/4ZGIh5/tk4m39yKpC2rEQaxIU5n6+F2WiOJ8gou0QaQzbJbsS6174vsuzkHe9Y0j6Wg5s+H/0+fMFO6ne6/fUZIAlt8pQ3qC4f9ZrF4uFHI8dgjeHYnidcfOIB348dIjPpg7O+/2V9IQBErZ89WeqUWcvOBRriKJP5/Bw5Z8/T91vUiDvaSl42Q/sszvn+IjRw6x5OxOjdSSeeTPqbOjRymYsXenwObJ33/cu773wocoxbXn5V3XY1i8GQWZoJ9oiZ3b0ovj2EcPuieBTAYmf2wwCeG7OfnzazFwLYB+A3QwhRN/clcEkNDoUQQgghlosLeHN4KoSw+yIP+1kAHw8hLJjZW1GtKvfii9mhBodCCCGEEC3ALOpQcJEcAbDYoiTi+xxCWGz78SEA/8/FHlSDQyGEEEKIFrAMJtjfBbDLzHaiOih8A4Cfc8c02xRCOKvJei2APRd70EtqcMjahzgypFVINNAOJDN+fVdfVCuVH/F6h/6dIy4eecblfvtr6A3zU37EhUfMaxQB4PuP+WLn9xz2noMPUDw26TUxrHNKpaOzrYb7ve5j10avzXvaJh9fvc57TY32borsczDnPQWTXX4fyen6PoiBNIiWjJkllvL3lDWICfI5ZD0cSN9ilfYVdW/E/ORCZBn7jHFOl4v1PQjZ97BSaOL8SWPIOs5Kzmv1KnmfBwWLei2WSl6rxDaGXWyVSBrDxIzPJUwcjxyjcsZrXEOBdIukfcoNb3Tx1pEdkX32bBz27aSZjOx7OE1ecoU59j301w4AijP++TfSILLPYaO4SDrcTuGK9T2RZdsGfFuHcv7ck2e8xhDjXtdcPOI1huP7olIt9jFkjeGjYz7XWGM4VWqkKAQGSe++jjwzN5OetY+0lZnuxkaehRlqF8Xsv3iSPDdZgwgAp0g/yd7AEa9gunbG2skYn9WhvNfddnV77W+yi3T83f653NTjv5dHB6J+uK2i1SbYIYSSmb0NwJ2oWtl8JITwgJndAuDuEMIdAH7dzF6Lqlf0aQC/eLHHvaQGh0IIIYQQy4Wh5bOVEUL4AoAv0LJ3L/r/mwHc3MpjanAohBBCCNEClsMEux1ocCiEEEII0SJa/eawHazpwSHXlE3nYmrCkucS67FS2Urd9WnSeeQGo1qGvm1eYzh8zU7/mWdc7+LK5b5q4KPzXoPxlf0nI8e48z6vqzmw/7SLJ094LVVh1msQ2dsuk/eaDgA4Pui1Tk8c8/q/w6e9FnLycq9XuZY0igCwg7Qf64Z8zcxk2mtFElnfbit7zUwcgf2/yLewQppC9tgLpFFEKnqPZ2a95qg7vzKare+83LsVsL4QACrl+jlcKfscTpLetNyExjCRpfrVPT4XQo/X3ZV7fb3wyZJv02yMlyJLl3Ipf19TXb7dCdKGGtUfLh4/GDlG4Yj3s5sb81o+1lt2DXoNXH6H91sDgOGr/fP9jPW7fDvo/oxN+3v47Vmvq52N8a1cmPbtKMxMRrZZTEPNYcznR37cVezCqc/dVPcYreAdn73fxf/x8pHINn0kNs0H0tiS1pTv++wh33dOPU4aRURrJTfyMWSNYZq6oNFcVB94RY/vZzZe5f0bh69c5+KeUX8tMr1edxdH4Yyv8Tx9xOv/tu313y0n9/rvkT1norl3ZK7+ufO1ydDAKU3XNk63nx0gLXq/v4fpXq9ZzpOGeZC++0d7o304ex9eqO9hosWaw3axpgeHQgghhBArRusrpLQFDQ6FEEIIIVqAwZCOqfKy2tDgUAghhBCiBZgBafbZWoWs6cFhccbrQlIxdU+ZZCbGJ2/xPshDi/UR3Ru9TgQA+i8fdXH2ymtdHHY8y8WPL/h9fmW/18x8+u7DkWMc2uu1IxOHHnHx7Jiva1wukAcc6eoyvdHzKJe2utjo1flB8uUayHtdTU/Mte2hur0Z8hjs7/E6G65znKDzQIj6cAXyfEOSagmT5pA1iex7GNEgAkg2rIy6PJTmvZ4nXnNI/oDkzck+Y+x72KjWMhD1MUyNeE/LYp/3HTu54I95inRJkwtRzWGa8m2IdFsZOo8030fyLGRPQwA4c9B7H7L+bPaUr8ecIH1m7xavWQSATRNex9V//Std/PT1Xts0Puf1UkcnfLunKQaA+SmvES50ea1YifSWDN/juFrMJX7WVoDtw15Ht74n+uz1Ub+SmPXXuzJ+wsWlE74vPHPI3/PJA1G95fHjXqu3VI3hzm7f7mtGo/rr7S/c7uIN1z3VxbkrrnFxatTrW0Mm6gHJWMHnQfmoz9cND/3Axce/472Uc1+P5nfmiNeePzbjc6eRBrGbrm0+5vrnhv09yq33mubksO9vEt1e49yX9/1PXB4Nxyy7EAzN1bnudNb04FAIIYQQYqUwRP8huxrR4FAIIYQQohWYIaHBoRBCCCGEAPTmcFWQTC9dY8j6rATrmLq9HrBrwOs8erd5/zYAyF1+hd/ntqe5+ETC64W+f9RrLr66x2tmjj3ufaEAYPLIYy5mjWFx1vu1Maw5TLBOL2YZ/+soBKqpSVqT2Zi6nKwtG8j6Y+SoXV1Ue5k9CuM0hzCvsYtoEOkYrGvkmP32ACBm0YowR35rC1NRrRhTLvjrEfU1rK+f5FrLAJDoJx/DHu/BNh68r9jjE16/9uhpr+WLy5Ve8jEsVryXZHfaH6M37dcb1dgOpahHZmHKt2PykH9uJg/4uES1aLv2Uv1mAAsTXue1g9YP/5j3D336hvUu3r/Ft+ng8ah+cHrCPxfzk75fSpBvIWuOmTjNYXlh5TWH60mrNxCjG88nqc+e8+daJM3hzDGvx5w+6rWnk4e8hg6I1hw+Xahfj5x9DFljuOvVXk8IAJte+h9cnLrmBS4ujniN4RNzvg1xzwyT7/bP/sjmZ7i4Z7v/bsqOfsPFmd5/j+70816XWKBn5nHSGPK142s7FHP9ezb5e9RN9zC10d/j9PBmF+f7vCaxvyuaR6xhvlAMQEKaQyGEEEIIAZydrSwrGyGEEEIIUWMNONlocCiEEEII0QrMTG8OhRBCCCFElarmsN2tuHjW9OCQJ4uUi1Fj3VCuL+JNkuFvqtuL3vPrvGFtbuuWyD7So5e7uDjotzl62gt29570gvMTJNYvzMWcR8WLfJMZEuPzZBKe6EGm17nB6MSanmFvPNo76I+xfijv4nV9NEGgKzrJhWd18WSPEi3IZPwx0ESZIp5QwibXJTpmkY5ZJqF3OW5CCi3ryUe3WQ5mT/lJAjNNiNK7y/6a8aSKMgnGja5xujtatD456CdRlPq98fuRU37izA+OedH5Q0e9iJ0nMwHRfGLR9wiZrq/L003oqj9BBYiaQRdn/LM5TgbUx2lCVTLGhLxA++Drt33TDhdvu+6nXXzlOt+PfW84mlzjJ3yfMZPz5tFJOneekFIuRI21mbhJKstNP01AyaWiz7st0ASdOZ9blUk/SWh+zE9YmSET5tMx95AnUXB2DqZ9u57a6/tXNrjmyScAkHjOy118NOUndT182J/no+P+e2FiLjrBihmgSReXD/pc2jVMRtzP8Xnkp3VUKZzx7Zj+zF4XT5Kp/jj1UXxt467/AN0jvoe9dI85B2ze9y/5tD8vABjOt84Eey2Uz+uIMzCzATP7lJk9ZGZ7zOx5ZjZkZl80s4dr/x1svCchhBBCiDZhQDLR/F+n0ilNuw3Av4QQrgJwLYA9AG4C8OUQwi4AX67FQgghhBAdicGQTiSa/utU2v6zspn1A3ghgF8EgBBCAUDBzF4H4EW1zT4K4GsA3rnyLRRCCCGEaEzVymb1iw7bPjgEsBPASQB/Y2bXArgHwNsBbAghHK1tcwxAVAQHwMxuBHAjAGwdGcDMJ//o3Lr8Jm/MW5qJmrhWGmgOU1lvep3p8xqN7Dqv1Uut81orAKj0eT3WGdJ8HadC5aenveaiTPqrVDqq3WONYLrbG2uz5jCd9TqmXC8ZCvdFdWV9I/7cd6zzuo3L1/t97iAN4mjMPtl4NJfyDxU/YxG5Xypbfz0AluEVaAHL28okIIxoEGNShrWRjz4WNUO+UBbn+EgijU9veNKolvU8xZgLwP825WuaaWB6bfTbR7KnJ7pRv8/xk2TQ+yDpaL/9qL8+bOzMhuoAMEb5lCcD+009/lldl/fd23DGt9uyUd1RooFx/ly5vl5qLkaQyvmU/85h366n/cDFfVdc5+LLh/yzvX0kqjl8vMfrpfj5TpEGuUSa4wt7P0DFAAAgAElEQVTRE6577fue3H//6HOWvIMai/O7f/1mvPvOh86te9kur7vLpqJfujbv+/XKtDdMLk15fRqbks+Oec0ca+IAYDpGA7uYdWSqvO5K/72w4Tpves0G10BUY3jXEa+Tu+txb/z8wGG//sx0VKvH9NIz8rQt/nvi9A7f7ueO+jZtimn3hiMHXTy296SL133vuIvHiz7X+NrGXX++R3wP+R6nKQes6HMkm432YWyMfds390e2aYa1YoLdCe80UwCeDeADIYRnAZgB/YQcqt8UsUUoQgi3hxB2hxB2j/RFO3shVjuLc7w/pnKNEKuZxfmd75e0XKx+kmZN/3UqSx4cmlm3mbXyG+owgMMhhLtq8adQHSweN7NNtWNuAnDiPJ8XQgghhGg7hupPy83+dSoNB4dmljCznzOzz5vZCQAPAThqZg+a2fvM7CkX04AQwjEAh8zsytqilwB4EMAdAN5cW/ZmAJ+5mOMIIYQQQiwra2S2cjOaw68C+BKAmwHcH0KoAICZDQH4MQDvNbP/HUL4nxfRjl8D8DEzywDYD+AtqA5c/8HMbgBwAMDPNtxLqCAs0jOwz2EqG/UxYp/DiL6KPpPu9ftMDm+kOOoEVc55f0D2olsgzUWGvLzyvV4nUipGi77neut7NLFOsSvnb30f6QG3DHmNEgBsH/E/228jn8PRXr+P9aSD6stEXzjn0/6fTuxjlmRHMdJvBfPbF2I0X/O0jLfhmPWDrBlrRnO4XD8XVAAUFh1rhs+FjSIBZMhLMkOfyZPPYaD1SdLhJWK0epVcv4vH532OPnDU+449etBrgiZJU5SIcZEtkb6vn3wNOR+39fvnZpDamOj1zyUAZEiakiKfPb6vfL2n4rRp5Eu6+ZDXio3vO+Ti/jEfD45udvGWoajmMEsebRnyeGTfU/Y55Zh9U1eKZMIwsOi+5qnf4lwGACt6j8bKnPfEK8749YUprz9bmPIaONbxAlEtb5basZnyZPjKdS7OX3Wt39/IZZFjsI8hawy//ZDX8p0iTeLcmaimnmFteSOd4lDWP2MjW6Ltzl1xjYuHr9zj4s0Pkr541vsxztMzFHf9+R7xPeR73EU5kKAcyeSjeZTl754L7MMNhgQ6+JVgkzQzOHxpCCHirhlCOA3gHwH8o5lF3WSXQAjhXgC7Y1a95GL2K4QQQgixknTyG8FmaTg4jBsYXsg2QgghhBBrndX/3nAJVjZmthvAuwBsr33OUJ1I/IxlapsQQgghxKrBUJVHrHaW4nP4MQDvAHAfoqUlOwNLwLqe1L2l015Hk8o29vKypNe4WIbquZJOKdnvvRRDjH9ShfQ+gcqYdpHWYT3p/3Zt8PucIA0iAJS5BjHtc4D0WUPkdzU64I+5sSd6jE103CHS2XDt5L6Mb0Oi4HVlAGBU4xWkf+M6yCFN15JqLcdp7uZKrPH02yyQxm6+RN6BtL4S48EXV295OaiE4Lz22Hcv7vzLgTSH1HGxhpXrCzOWi2oOQ7bXxUdP+yR/mGopT5z0mqCpk15PlYix7KlU/LP2OOXjQfLZvJI0sqNUBz1N9aCBqE65q9/3IT30XMVp4Bi+R5OkOZ494bVj5bGjLu7Z6Y/J3qAA0NNNfR3V+k12UT9GfWOi1FhzyF6py0HCDNnUk8dJswY8TnNYph+uSj4uzft+v0ga0CL1OXHPEJMjs9C+bn9PesgfMLnJ1yw+Phe9vlwrmX0MWWM4fsRrU+fGj9VpcZX5wY20ZKs/JvX7V2zwz/Wu4agWffOo1yHyufO1yU36voE1h3HXn+8R30O+x5wDnCNxecS5djFehWvgV+UlDQ5PhhDuWLaWCCGEEEKsYswuvTeH7zGzD6Fa5/jcFKcQwqdb3iohhBBCiFVIJ/sXNstSBodvAXAVgDSe/Fk5ANDgUAghhBCXPIblszJbSZYyOPyREMKVjTdrI4kkEvlFtSIvwKvLUl4fYV1eY2F5r8FI9HkdVCUd1WQg6fU8Zl4vwTUdtw/4fQyQxog1cEC00Df7g/ExRsgXbYjq0PZ3RfVFrCFMzo67OHHG17e0gvftwkLUhyuwNoQ0nqzhrHT76w3SHMZdG9YYTpNf3hmKFxpoDtn3MI7l6hwCvOca63PiNYc+ziRYp+jva7lRreUY7VmgGteT814/NUY6o5lxnyszJ71+Kk5zyMedIn3UgVNex3h0o39Wt5Pv4XB/VHOY3eD1Ut3r/T56Bv15Ds35/G1UgxeICrZL834fgTzaWNeYi6utTh6iaXre2ceQay1XqN5tXK1l3oe/H63Jd4Pvy/jXudijBH9FA7W9UvD9bYU0tvzMNKMf5uc7Q7q6TK/Xuwaq6z0bUz94gnKJPQjZx5A1hvNNaA6ZbJ/Pdz4mtymu3XxufO58bRr1jXHXn+8R30O+x5wDnCNxLeBc4+/UprFo/frVyFJ0k98ys6uXrSVCCCGEEKsYQ3Wg2exfU/s0e4WZ7TWzR8zsppj1XWb2ydr6u8xsx8Wex1LeHF4P4F4zewxVzaGsbIQQQggharTaysbMkgD+EsDLABwG8F0zuyOE8OCizW4AMB5CeIqZvQHAewH8p4s57lIGh6+4mAMJIYQQQqx1Wvyr8nUAHgkh7AcAM/sEgNcBWDw4fB2A3639/6cAvN/MLIQmNFDnYSmDw2MAfhrADvrcLRd6cCGEEEKINYNZq61sRgEsFmUfBvDc820TQiiZ2SSAYQCnLvSgSxkcfgbAJIB7sMjKppOwZNKbVDczIYVE7o0mpCRoQgpPmAhpL1iPI0tmmyNkUM2TSRZI5B4nlGUj7RzFPTSZpJ8NqnlyyaQvlg4AmDzhwvK4j4tnJlwc5r2wniefANHrzSbjqQ3b/PY84YcmqMRNyGAR9WmeRBCZkOK3LzYwhY7jYgxU6xGCP0c+37kYNTfd6oafKdP1KBdJ7B33XJHgm28DG20XZr2hb+HM6eg+iSRNopgZ7Hfx4dNesH+MxPWTC74Ngzy5CUBq3aiL8xuHXNyzwZt1bxijSQIx158nMPXRs5nKkqk1GfHz90ycoJ9N75N8DLp2PLmETbHjJh1FjLMXx2vA102IVmAhwJb2wm7EzO5eFN8eQri9xc1aMksZHG4JIeinZSGEEEKI82CVUuONnuRUCGF3nfVH4EvZbKkti9vmsJmlAPQDiHnD0zxLna389Is5mBBCCCHE2iVUf0lp9q8x3wWwy8x2mlkGwBsAcLW6OwC8ufb/rwfwlYvRGwJLe3P4AgC/qNnKQgghhBAxhHBBHsvn310omdnbANwJIAngIyGEB8zsFgB318oafxjA35vZIwBOozqAvCiWMjh85cUebNlJJJHo7quzPvqiNKKtYQ0caQ7Z5DqwBo4Mr+PoSnl9Tj98G1gfyBoj/jwQLQQf0RBO+RhTpB886d9SL5zkt9bA/HGvbV2Y8CbXhSmvMayQVs1irn+mz5tY927b4GJnag7ARlgP6PfHhtcAME4aQ44nF3w7WXPYjOl1ms4tuYyV1xe3J2rgG20re1qzLi5ipD3jr0+ZitqzSTMAJIre5DqXJs0b6Wg5F9h0uRKjT10gXSKbAk9N+TYcmfDx6Vl/n7dtGIwcgzWHPaPrXNy/3RsNj5zybSgfj14bliEODfs+I7/e53iix+tuI5rQZkzYKQFZY5ikfi1F178cY4KdZJ3ior7SWjQ/M8CbzrN2NfbMzZ+rcTszZAhOucgm480YGPM94GemcMbnARcEyHdHOwgudtBLJu+5Xn/P5gc3Nm4okaPP8D75mNymfDrmO3Savgfo3CP9SYP8jbv+fI/4HvI95hzgHIlrAedaXEGFZrHm3gg2TQjhCwC+QMvevej/5wH8TCuP2XBweHY6dAjhQKNtWtkwIYQQQojVRQCWpjnsSJp5v/FVM/s1M3PTRs0sY2YvNrOP4snfuoUQQgghLk0CWq05bAvN/Kz8CgC/BODjZrYTwASALKq/ff8rgD8LIXx/+ZoohBBCCLEaCLDy6n9z2HBwWPst+68A/JWZpQGMAJgLIUzU/2QbsAQs86TP4OL/BwCkYk434ZcFiiukIQxpr8kIqWzd9XGwfiKTIc0LqEj4gtd0JKbPRPbJGsPK2FEXF48ddPH8E15TOH3E+7fNHov6zs2c8MddmPS6pNJ8/QciTQXYAaB/q9dbpbq9BqZru9eNGWlHFsr+WsUVhj897zUvp2Z9PEXrZ9nnj8QocQan+YzXwOQyUZ+45YBlMXE+j6zh4W3mKeb7OD8+648Z44GZmp908YZur5sb7vfPSSbvPQrZwzBOc8i6xOK8fy4W5rwn4WnyOTw16z8/WYg+q8N9612c3ew1iL3bjrt4ljSHFiOYCnST+rf7c+/fscnFqfVbXDxDOT1XbCx2N9bRpXw+JqlvrBT9tYn1OSQdl/NOjNETXwghBFQWKZSK9HzzswgAIUn9CunGU1nf7nTO9/HprI+5f46DdaBTpKubPuL12aUj+108sjk6j/PyQa+/ftoW3zeemWZ74a0uyvaNnLe9Z2GN4cioPwYfk9s0kovmRflBrzjjc+drE+cFupi468/3iO8h32POAc6RuDziXJsvXeCkkhZPSGkXS5mQghBCEcDRhhsKIYQQQlyKdPDPxc2ypMGhEEIIIYQ4H2GpJtgdiQaHQgghhBCt4lJ6c2hmXwTwf4cQfrCM7blI/G/9XPc4zoMwkLYmsD8SaxJJuxCppRx3DIpZlmQl0tXNe21fcsZrvCpjT0SOUWBN4WEfnznktVJnDnqfw8kDXjM2fdzrzABgety3c7qBH2COvNYGB6J1p1k7MlAgb0TSRoWs18RMk4nfiZlo2e8T015PdYz88CZJizZXqK8XidMc9pAmZudIT2SblSBOzhPxZIvUVvbXcIY0mXN032di9Kj5Mz6fNg56f8CnbfE6u4fXe03i7JjX9kWzL6p54xrPXL95Yra+nyXnDgAM5r33YWqjr+3du9VrdQtTvqWpbLRLTdAD37/T6xp7dl3u4nKvXz9dZF1tND8L9CySNDfiM8m1llnPyXWUAe9rWI2f3KZVPoeAz9d5Oq84TW2O+uBErtvF6W7Su/aRrrnPn2s/5TsAHGMdLrXjCVq/bS/V4N53v4t7dkY1h7uGt7v49I6hyDaLeYA8CaOaxCjsY8gaw+fSMXeRJ2f6lNdOAsD0Q35IMEbnzteGr12aUqc/HdU18j3ie8j3mHOAv6fj8oifK+4/msXC2nhzuBQV8TsB/JmZ/Y2ZbWq49RIxs6SZfd/MPleLd5rZXWb2iJl9slY2RgghhBCic6lUmv/rUJoeHIYQvhdC+DEAnwPwL2b2HjPLNfrcEng7gD2L4vcC+NMQwlMAjAO4oYXHEkIIIYRoLWdnKzf716EsyX/AzAzAXgAfAPBrAB42szddbCPMbAuAVwP40KLjvBjAp2qbfBTAT1zscYQQQgghlhMLlab/OpWlaA6/CWAngAcAfBvALwJ4CMDbzew/hBBuvIh2/BmA3wHQW4uHAUyEEM7+cH8YwGjcB83sRgA3AsC2bduQuvb/OLeucOowbRwdC7PmMKIxZJ0TrWeNV7mJeowZ0iCxYWZiwdemDJNew1F64rHIPmf2+2UTj3pt1Pgjfh+sMTx1wmunji9ENROnSYvHug32pxoir7/eGB9ErpHZNei1eslhXwt0Luu1a6fHvD7o6Jmo7ubIhPeiO0rxNLWL9VtMnOaQ9Snve801Lv7VunusD+f42w/sO7fud9KXNfw8n01Uc+hj1pLOjfnrNUOemAAwcHCfizdufrqLn0Oaw+9f5rVNc9NPcXGczx5rDFk3x/B9ZE3RXMx9rgx4zWF62Ctoslu8BnGA6k53DUS1psmM1+r1X+67ssxTvP6s0OePefo41ZCOeY7ifNsWk6B67OxzmG7iDUYq63Vc+z/40+f+3z7yxnsa7uA8RPL7+U/m9Bce8lrp+Zja6X1U3z5JtalTfT73+B7lh72X3+AT3j8TAHpS/rtjnHSgJ6m/PPaQ1+UOf2ePi7Oj34gcY8Nz/PV97qj3LRzK+jy6YkOviyfmGmvkuFYy+xiyxnBDifwa74+2+zid28m9/tz52jB8bQdj6jfzPeJ7yPeY65OXKEfi8oj7h1tefpWLfz/yifOxNsrnLWW28o0AHoypofxrZrYn7gPNYGY/DuBECOEeM3vRUj8fQrgdwO0AsHv3btV3FmsO5bhYyyi/xZriUjPBDiE8UGf1qy+iDc8H8FozexWqZfn6ANwGYMDMUrW3h1sAHKmzDyGEEEKINhMQYio8rTZaUvMohBCd3978Z28OIWwJIewA8AYAXwkh/DyArwJ4fW2zNwP4zEU3VAghhBBiuQhYExNSOtkE+50APmFmfwDg+wA+vNQdVHJehxCnOeRlLEVgLU+xyN5+fvvIj+4AjCRqrHTKsa6RP0D/CqnMRmsrz415DeGZw77W8tThKRcfP+51jYfmvEYiTicy00BP2UfakSSdB3tVAUDvFq/xyu/wGrryyA4XH5/x7XqM9IP7x6IOeQdO+XM9Tb6Hc3Suge4516nlGADSydbUll0qzdSBZRr6HNJ9Zp/DM4ejPodz+/e6uPcyr0F81sZdLn75072WdIbuQborqjksUI6ydx9/hrWhRTqvhRjNIesvkw18D7upc89tjD6blvV6qfS2K1xcXudz/uSs3yfXhD7ThOYwSfkY9Tn0555INPblzMbURl9uxujcZ4tRnWno8159yHktXqJ/2MXZYf+90L3Ba/2GDkRrbg/RNZ+k7wHWIO6jPqbn677+cKb33yPHYH+4Tde8wMUjW3yesD4wrq48kyc9H9dKZh9D1hge/VK03Qfo3Pac8efO14Z7StamD/VFrz/fI76HfI85B0KXz++56eigbDLmuboQQggIxdX/5rCjBochhK8B+Frt//cDuK6d7RFCCCGEaJ5LTHMohBBCCCHqENaG5lCDQyGEEEKIVhACApWjXI1ocCiEEEII0RJCR5fFa5Y1PTgsW33DagAo0cIiCbtZ6M0a9nKovz0AGE3MyEYMaf1EjWzWC6zZ1JUF7tV9eLF4goy2A51nxLw7biYNwQXScyR6H835633ZoDfa3fQcPwkBADY9z5tFp66+3sVjGS80fuSYN6jde8LH+09EDWyfGPeTVgpkFlsq1teHJBqI+wEAbar8naP7HDdBhQvdMzxBhU2wp2b89eo9HJ10MbHvkIuz27wn8sZnr3Pxi8kEm/lqd/SCnjjtJxuVqZ35Xi9kH8j7ZyKdrD9BBQAWaEZavssL21PDm/0+U/4YcT8nsSFvZcBPPZjJ+msxMeX3MT7v47lCY+F8kiaHpUj0n6FJSGU25o/Jo2x+5ZP8ND2rkzET5WZLvt/ppYmIycH1Lu7e6K93zyZ/f/q3RvN785Q31+cJKaeoQMARanfmCO3z81Fr4MIZn98bjhx0cf6qa32bNm13ccg0nlRk075/LN3vJ6BM77vfxWxwzZNPAOB+Ojc+d4YnoGzO+vvXv9U/c0D0HvE95HvMk1Fnyz6fJ2Imn5yYadHbvhAQinpzKIQQQgghgEvPBFsIIYQQQtRDE1KEEEIIIcRZzppgr3LW9OAwSXbTszH3q9BAc8jrSxFNot9fvLqLzHepHbyPoZw33u0a8CbEbKILAP1kjF2Y8sbPBdKNVei8MmNel7euKyqoZfPSfjJUXbfVayU3X7fFxaMvfGZkn9nrX+niySFvlrznhNfh/PCYP8/7D3vz7yMn/XkDwDwZ6RbpBjQyvY7VGBL9vVHj1pWAdZ+ZRPS+saaQYekdm2JPkbav/3j0Go8/cty3a/19Lh7s89rRy6/yBr/Jy/36kRh9296TXi91etrrwDKks1vf5zWv/V2+u+tKRY3LWXvL5rnlbt/ORMofI+J4D6BMGuICPd9j0/7ZPE7ndYbytRBj3s2G31GdrI9DJUXr6z8DAJDKrLzR+xiZSZ+Yjmq51pM5dzcbl5MeLbV+1MW9WydcPDfmCwYAwIZTvn9ko3h+xviZeYz0bIVD0WPM3eGN5Mf2nnTx8JVe/9czOuLiTK83io6jcMY/u9NHTtU95rGHvOE9m3sDUY1hkfoTLo6wjbTAG8jgun87Fa8A0Lt1g4v5HiboHpcoB6ZIExqXRwfGov3ahRBkZSOEEEIIIZ5EmkMhhBBCCHEWlc8TQgghhBAOvTnsbLJ5r2U4MRb1wGON4QLprVjeM0vahSKZXcZ5p7GOKU16oLmM30eB2jDc7f0Bu7dG9VjZrD/XDT3sC/WoiwcfP+rijSe8BqY0H/2XTyrrtSL59V5LNXjFVhf3XON9uRJXeQ9DADiZ9Z5v9x319+ieI15T+N39XgNzlHwN52K0JAukiamU6mvwWFsV1SRGP/PZtz6v7j6Xi986s8/Ff5h/SmQb1kM18rTkFGbfw+mYazx12OdPbo/3Q+sauNvF3Wmfwzt3PMvFvTG6o13D3t/zOOm4FqidrClkzWE+Hb2RSdIMlklpa1nvwVZO51wcUtFns2B+2cSc70MmFzj2HmzzpcZfNKy3jPoa+nNPkKawQjnSlfPPOgD8+zt+rGE7Ws0tL7/KxR+86/HINqMRbanX//aQTjS1YZuL85NjLu6biPocLpDPYXHO3yN+xg7O+j6HNYiPz0b710nyW133Pa/j3fwgtZO0lpnu6D1jWHvOHqZPkP/fScrF8WJjc+dGGsOt5H87tMvrA/t2RP1w81v9Mr6HFbrHs+ZzYGrBn2ecp+H7XnNNZNkFEQIq0hwKIYQQQggA1Z+VeZbpKkSDQyGEEEKIFhBCQKXYuIpRp6PBoRBCCCFEKwhApYkSl53OJTU4TMZoxeZIfzZHmoozpDFcIP3PbJE1h9HXyaxrZF3TKaprytqo2RLVjO3yekIAWLfjR1ycGvIeg8NPPeGPceKwiyvT3usL5ajOyXJe15gc9nrBMOw1h4UBHx+ciuow7jvgj/vdg+Mu/uFBrzkcJy+qWdICLcxFH8pSgwfVEj4xKsHrtdjnsBnfw3bRE+PdN0ciQvYxZFiTyHoqrisLAN3kRzlB9zXV7TWv7MOXXfA+cuu3PS1yjCHK6W39Xss3Q+1iaSVbEGZjOoSuVNTfz+2TfA0D9aDzMZpj1ilPk7aM+xiu986wvhAAcqQx7O6q37Xz9Y8cI9uZXw1HJuYjy472+z5gIOuvRa6XPPIWfB+SHvV62YFZn8sAUJ73GjWuVY8HvF8g1zhnDeLpQrR/ZT3fONXn5X3kJv214O+VOPjZ5r6hUR32uKzhWsmNNIbrn+b9GYeu8N8jrF0HgPTo5X7BoP9Mme7x6Rl/fQ9O+P5l34no/INWEUJAJeb7c7XRmT2AEEIIIcRqYwV/VjazIQCfBLADwOMAfjaEMB6zXRnA2coEB0MIr22075W3vBdCCCGEWKOEcqXpv4vkJgBfDiHsAvDlWhzHXAjhmbW/hgNDQG8OhRBCCCFaQggB5ZWbkPI6AC+q/f9HAXwNwDtbseNLanA4OtgTWfbgMa9pY/3PJPn9sSaRPcjYJ7G6jV9WbqDryKa8huMY1VodzEb9rLgW7VDe6zp6hn3tye6Nz3BxhnSPcU1k7dl0gTQy8/5aPLrf6864Ni4QrY18gGonnzlNNU1JY1ggfVCp4LcHgECGpJbw1zdJ3nSJRDRPFvOtm19Sd307efvUvsiyW7t3xWz5JKxDinj90fpijE/iDOmhuk6QrivnNVnJtO96BkgXmpuNes2x7migz+d0f85rcStd3oMQyagH4VLhM+fa6xwDwDzpmtk7lXXKkftB+rVMKqp5zZPuq0CaQd5HuYEmsV2+nY1g30Mg6n04RP1jF9334T7yzNvktX2ZQlTXOFSp/4bHqP9MP+x/2eumeuTsJwhEdYjsL8p6wEb6wAshTbJF1jCzvhAANlOuca1k9jFkjeHQU7f7NmwhfSGA5KadLi7RPZwg28KT5N/4OGkOb3311ZFjtIyApb4RHDGzxUawt4cQbm/ysxtCCGdNi48B2HCe7bK1Y5QA3BpC+KdGO76kBodCCCGEEMvFBVjZnAoh7D7fSjP7EoCoMzjwLjpuMLPz/YthewjhiJldBuArZnZfCOHR82wLQINDIYQQQojW0GIT7BDCS8+3zsyOm9mmEMJRM9sE4ETcdiGEI7X/7jezrwF4FoC6g0NNSBFCCCGEaAUBKBdLTf9dJHcAeHPt/98M4DO8gZkNmlXrCZrZCIDnA3iw0Y4v+TeHV2/0NVz33u9rDs8W2dewfjwdoyVhn7NCieP6/8pgX7NcJnrbekn30ZflOrLk/UVxI50ZAMzRuZ6mmsWHT3v934FTPj42FvUPOzNeX1O4MO21OoVZr1GsFLlucrRmJmsME+n62rNSwe/j3j/6qbrbdzo3zTzs4vf3X+lirunKsjnODdaeAlG/tIUpfw1nT/n7nMp6DaKR717czzK5Ge9Hl1o36ttJvpuJrNeOBq6DnPb+a9WdZGgb7zGKVMxnFm8f86NOaFDLmq9vmvSBWXr+WV8IRDWGDPcP7It42096DfJq4q3P3eHiv737oIvzGd/HZ3p8je6+fp9HqW2N3/qMkGY2mfV509Xn8yZ/wPdbQ4eimtrT1Pex7yE/p9Ga6XUafLadpClkP8Z++l4YpPrjQ3ReANC/1dcb76e66FwreeByf70z270mOrXN908AUKZ7dCbh7+HJad9fsMbwbT/qNYvLScCKls+7FcA/mNkNAA4A+FkAMLPdAP5zCOGXATwVwAfNrILqC8FbQwidPzg0s60A/g5VIWVAVYx5W7P+PUIIIYQQHUEIkZcWy3eoMAYgMkMyhHA3gF+u/f+3ADx9qfvuhJ+VSwB+O4RwNYDrAfyqmV2N5v17hBBCCCHaT1hRn8Nlo+1vDmvTsI/W/v+Mme0BMIpl9O8RQgghhGg1IQSUVVu5tZjZDlRn0dyFJv17zOxGADcCwLZt25a/kUKsMMpxsZZRfos1xQqWz1tOOo9NbOgAABIRSURBVGZwaGY9AP4RwG+EEKZskUi7nn9PzSzydgDYvXv3RbuCshi8QmLyBZo8whNQJmajEyJ4G56gMheZsLI002wganIbNc5dmoIgbpIMt3N6xp/rPJ373DStn4nqMBamvTC7MOsnHZTmvHE2TzipkMF1IhEV6yfI5DoyQYXiTL7+pIOVptU53tcgF+bopw6eMMExEGOUTTm/QGL7uXFvVJzKna7bJiBqLJub9xOcKmScnewfdnGiz8chGzU750krCDSpxfy1SyYot2K6KX4Wjay002ygTJNzuijuaTD5BIgaZXMfwhNS2kmr8/uh477P6CHD71TCmzQnuvv89gNbI/tM0X23jO8jhrp8nB3wkzRyw8f9MTb5AgEAMEBG2bM0gY8nefEzFjdRjOEJKGnKpa4+n8/5YT/xo5sMrgGgZ5M3n+/d6t/j5Lf6CSlsZs8G1zz5BACm0/4enaDvkr1j/trdd9R/j7zhmdF9LhsBCM3MDupwOmJwaGZpVAeGHwshfLq2uCn/HiGEEEKITqBaPq/ceMMOp+0TUqz6ivDDAPaEEP5k0aqG/j1CCCGEEJ1EpRya/utUOuHN4fMBvAnAfWZ2b23Zf8V5/HuEEEIIITqRUAmoFFb/m8O2Dw5DCN8AEBUxVYn49yw3P36110f8/fcOubhAuifWD8aZYE/MFutuM7Pg4xLts0LHLMfoAcsNpsSHCsdkoNrEMYrcTjJpLcx6jUxx3mt/SnNeFwIAJd5mwZuXBtIUcsz6wRBjcB01wSa9G2nP9v3FT0T2sZb4hbE9Lv7YyFNdnCVdUpH0hHGaw7hliykXfK4USTPEmsRU1ucFACTI/J2NsyNKUc6dso+T/TEdeJ7iJB2TconNj9lkuLqM9cB+fTpBJtdkPLxQ8scsxmjL+BisKWTd49uff1m0oWuEW199tYt/74t7XRzNVX/T1+d7wfQO+okySTJL7yLdYrLfG3Hn1ntdXvexqMZ2fswbZS9M+GegMOX7xuKc748rTfyUmSCT63TO53emz2tuuwZ835gd9gbXANC9ccjFqfVkKr6BJhkNerP6Up//zmWDayBGY0gFFh445vXGf/hK36etNJ38RrBZ2j44FEIIIYRYC4SANaE51OBQCCGEEKIVhKDZykIIIYQQokoIQFmaw7XPm57tPa9u++Z+F7N/WJwnIfsDssawMOf1FIUG2j7WJFa3IZ1iKdRdXy7x9t5Dq1Tw+pa4bcoL8z6mz/A+yguN98maQiaiH2zgYQgAyYzX0aQofuyvf6buMdc6P3/KaxA/s/FpLp5m382YfxSzf1oiUV+DyDra0py/78UZr0EEgPSMz7fClNcdsf7PUmmK/fpKjD41ScsqpC0zjiv+WU3H5F8qESj21yZP+sB59nNM8/WPHiPbwLuS+7FLife87EoXswaRPTrLw1Evv1Le6+AGRrxmM9nlP5PuHfTrh73OLrUx6szWOznmjznlNYhFyv/SPHm+NlGVg3W7KXpm0t1euZvq8xrDBHmHAkBycL3fhuJKt/9Mudf7IJ5e8Nf/1HT0PPaRjyFrDPketxX5HAohhBBCiHOsEZ9DDQ6FEEIIIVpAgGYrCyGEEEKIswT5HF6SsD8Y61fiapayDpE1g6wHLC5wTBrFhRjNIe2zXCB9Cuv/aD3rAVkLGLePSpE0L1z3uEHcDKwhbBRn8lEfrjR5kB36+7csuR2XEq879oCL//myZ7p4LqZGdleGvSS9Bi6Vo/Vs9kfE/cu7TJqq0rzXJZZJgxUoxwN7aMZpYHkZ11omjaFRTifTUY+2SJ1zMkMsUf/Q26DucVwJXa4Jz36t4klYn/auf/aa21JM7u0Y8HmwrtvrWYe6/fXO573mMEG6u/Tw5mjD5ryOLj3t6y93sU9syT+HoYn+lXW3IF1uIkd1p3u8PyNyUQ/ISs73uSU691nrcvHErP+uOj7tz+OhU1GPU66XzV6WnYSsbIQQQgghxJOEzi6L1ywaHAohhBBCtIAQoJ+VhRBCCCFEjaAJKQLN+Su95oP/X931FRIRReoek89cJabucWne61G4rvFSPQlZT1g97tI0haFSv95zHEZ1ZiO1bMmjMJ3ztT9P/NNvLfmYoj6v3H9vw23+7brnLWmfqazvepIZf98TcUWKGxDIH5BrKYeIRiuqnWQdF/sYssYwpP36VEyzWQ9YoriL9JeROsmk33zR5SPRg4gLppk6vO//1mMuZg3iNor7unx+9+VJk0hefwBgC77PtqLvkxNF34dbmfI3NNHfms+lkPSaw5D2Pocl0tyGrNdvA8AseepO0Vuz8Tn/zByd9lrhQ5P+PN/63B2RY6wmQgiR7+zViAaHQgghhBAtgo3VVyMaHAohhBBCtIAAoBBnKbDK0OBQCCGEEKIFBMSXGF1taHC4Anz2rV6P9dL3f8PFgV5Bs8+hNahTGwfr/SJ1kVk7Va5fa7m6z3LdmGH9YBxcG5njdLf30Epnvcbw5B3vaHgMsfz8x+94Xe09r3mpi4tzVIM4xzVefWwxmkNLsi4xUXd9Q2Lyl3M6ojmMrKdnIkb3laLngH0OSSqJLAkXRwd9zouV520/utPFN33+QRdfsd7fo/Xdvh9b3+PjgWz0qzeX8h6DWerrMnmfF+yf2cy3BI9Z2IOX33jNk55wdiqq050kH94T0/6748gZr5U8MuHjW15+1XnbuxoJQW8OhRBCCCFEjYAgzaEQQgghhKhS1Ry2uxUXjwaHQgghhBAtIATNVhZCCCGEEDU0W1lcMF962wvqrr/+D77o4iQL8WMmqFQqXszMptY80SPZYDJJHDxJhQ2qG8FtAIBklzdZTZHJ9anP3bSkY4jO4Dmf/VLd9ft+5fUuLs54kXoyHe2aeFmq2+dKMuvzy9LRfGsIG7eXvNgeKZqoxRNQYiakJMl4uJdObV2fJpysNm599dV117/js/e7ePuw7595wgoA9NMklXw6WTdO8/dCEzNSeMxSpNlQ82TezJNNxmajExVPz/lJKmM0IWWtTThphAaHQgghhBDiHNWfldvdiotHg0MhhBBCiBYQABSlORRCCCGEEFVkZSOWiW//t5dd9D523PAJF7OWjzWJlZLXjVSKjU2wGyGDanE+rvjApy56HzOf/CMXB85ZNmFnjWwTmtlQ8FpIzM+4MLP92ob7EJce73vNNRe9j3ff+ZCLB/JpF2dTrEFsLDos0u+dFRrE8KDm7c+/rOE+haeCtaE5XGJJgZXFzF5hZnvN7BEz08wEIYQQQnQuNc1hs38Xg5n9jJk9YGYVM9tdZ7slj6U6dnBoZkkAfwnglQCuBvBGM6s/RUwIIYQQok2cna3c7N9Fcj+AnwLw9fNtcKFjqU7+Wfk6AI+EEPYDgJl9AsDrADxY91NCCCGEEG2giAqOVuYbb9gCQgh7AMCsrqTggsZSnTw4HAVwaFF8GMBzeSMzuxHAjbVwwczu521WmBEAp9rcBqAD2mH2O21vAzrgOtS48kI/qBxXGxrQCe1QfreetrfjNzqgDeiMNgBN5vg4ShOfwrGjS9mvmd29KL49hHD70ppWl6bGUkwnDw6bonYRbwcAM7s7hHDe391Xgk5oQ6e0Q23w7bjQzyrH1YZOb4fye222Q23w7WhmuxDCYIuP+yUAG2NWvSuE8JlWHmsxnTw4PAJg66J4S22ZEEIIIcSaJ4Tw0ovcxQWNpTp2QgqA7wLYZWY7zSwD4A0A7mhzm4QQQgghVgsXNJbq2MFhCKEE4G0A7gSwB8A/hBAeaPCxVv5Of6F0QhuAzmiH2vAkrWpHJ5yP2lClE9oAdEY7lN+tpxPaoTY8Sae04xxm9pNmdhjA8wB83szurC3fbGZfAC54LAULa8DJWwghhBBCtIaOfXMohBBCCCFWHg0OhRBCCCHEOdbM4LAdpfbMbKuZfdXMHqyVsHl7bfmQmX3RzB6u/belU9vP05akmX3fzD5Xi3ea2V216/HJmhB1OY8/YGafMrOHzGyPmT2vTdfhN2v34n4z+7iZZZf7WpjZR8zsxGJ/tvOdu1X581pbfmhmz27yGMrvNuZ37Zhtz/F25HftuMrx5W+L+nCs3T58tbEmBofWvlJ7JQC/HUK4GsD1AH61dtybAHw5hLALwJdr8XLzdlTFpmd5L4A/DSE8BcA4gBuW+fi3AfiXEMJVAK6ttWVFr4OZjQL4dQC7QwjXAEiiOjNrua/F3wJ4BS0737m/EsCu2t+NAD7QaOfKbwDtz2+gzTnexvwGlOOXQo6rD/e0LL9XJSGEVf+H6kydOxfFNwO4uQ3t+AyAlwHYC2BTbdkmAHuX+bhbUE3eFwP4HABD1VE+FXd9luH4/QAeQ22C06LlK30dzjrBD6Hq4fk5AC9fiWsBYAeA+xudO4APAnhj3HZ19q38bmN+147R9hxvZ37X9q0cX77jqg8Pa7cPX41/a+LNIeLLw4yuZAPMbAeAZwG4C8CGEMLZ8jnHAGxY5sP/GYDfAVCpxcMAJkJ1Cjuw/NdjJ4CTAP6m9rPIh8ysGyt8HUIIRwD8MYCDAI4CmARwD1b2WpzlfOd+Ibmq/G5vfgMdkOMdlt+AcryVtDvH257fQMfleCvze9WxVgaHbcXMegD8I4DfCCFMLV4Xqv+0WDa/IDP7cQAnQgj3LNcxmiAF4NkAPhBCeBaAGdDPD8t9HQCgpgl5Haod3WYA3Yj+VLDirMS5LyfKbwAdkOOdmt+Acvwij90JOd72/AY6N8dXe35fCGtlcNi2Untmlka1U/lYCOHTtcXHzWxTbf0mACeWsQnPB/BaM3scwCdQ/VniNgADZna2POJyX4/DAA6HEO6qxZ9CtaNZyesAAC8F8FgI4WQIoQjg06hen5W8Fmc537lfSK4qv9ub30Bn5Hgn5TegHG8VnZDjnZDfQGfleCvze9WxVgaHbSm1Z2YG4MMA9oQQ/mTRqjsAvLn2/29GVceyLIQQbg4hbAkh7ED1vL8SQvh5AF8F8PoVasMxAIfM7MraopcAeBAreB1qHARwvZnla/fmbDtW7Fos4nznfgeAX6jNeLsewOSiny7Oh/K7jflda0cn5Hgn5TegHG8JnZDjHZLfQGfleCvze/XRbtFjq/4AvArAPgCPAnjXCh3zBai+av4hgHtrf69CVS/yZQAPA/gSgKEVas+LAHyu9v+XAfgOgEcA/C8AXct87GcCuLt2Lf4JwGA7rgOA3wPwEID7Afw9gK7lvhYAPo6qPqaI6r/AbzjfuaMqNP/LWp7eh+qsPOV3h+d3p+R4O/JbOX5p5Hgn5He7cnwl8nu1/al8nhBCCCGEOMda+VlZCCGEEEK0AA0OhRBCCCHEOTQ4FEIIIYQQ59DgUAghhBBCnEODQyGEEEIIcQ4NDoUQQgghxDk0OBRCCCGEEOfQ4PASx8xyZvZvZpZswb4yZvb1RWWOhGg7ynGxllF+i+VAg0PxSwA+HUIoX+yOQggFVB3l/9NFt0qI1qEcF2sZ5bdoORocrmHM7Ktm9rLa//+Bmf1FzGY/j1rNSDPbYWZ7zOyvzewBM/vX2r9Kd5jZQ2b2t2a2z8w+ZmYvNbNvmtnDZnbdov39U22fQiw7ynGxllF+i3ah8nlrGDN7IYBbAPw1gJ8D8NrF/7qsFbg/GELYWIt3oFq7cncI4V4z+wdUi4x/o7b8WQAeAPBdAD9Atf7kawG8JYTwE7V9JAEcCyGsW4lzFJc2ynGxllF+i3ahN4drmBDC11EtEv5bAN4Q87PDCIAJWvZYCOHe2v/fA2DHouX3hRAqqHYuXw7Vf1nct2gb1I5RMLPeVp6LEHEox8VaRvkt2oUGh2sYM3s6gE0ACiGEMzGbzAHI0rKFRf9fBpCKWV5ZFFcWbXOWLgDzF9JmIZaCclysZZTfol1ocLhGMbNNAD4G4HUAps3sFbxNCGEcQNLMuHO5mOMOAzgVQii2ap9CxKEcF2sZ5bdoJxocrkHMLA/g0wB+O4SwB8DvA3jPeTb/VwAvaOHhfwzA51u4PyEiKMfFWkb5LdqNJqRc4pjZswH8ZgjhTS3a36cB3BRC2NeK/QlxsSjHxVpG+S2WA705vMQJIXwPwFdbZaAK4J/UqYhOQjku1jLKb7Ec6M2hEEIIIYQ4h94cCiGEEEKIc2hwKIQQQgghzqHBoRBCCCGEOIcGh0IIIYQQ4hwaHAohhBBCiHNocCiEEEIIIc7x/wNOaQeu8AiKGAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plot(sim, component='all', z=0.0, cmap='RdBu')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## References\n", "\n", "[1] Beg, M. et al. Ground state search, hysteretic behaviour, and reversal mechanism of skyrmionic textures in confined helimagnetic nanostructures. *Sci. Rep.* **5**, 17137 (2015)." ] } ], "metadata": { "hide_input": false, "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.7.2" }, "widgets": { "state": {}, "version": "1.1.2" } }, "nbformat": 4, "nbformat_minor": 1 }