{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Spin-polarised-current-driven domain wall"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Author**: Marijan Beg, Weiwei Wang\n",
    "\n",
    "**Date**: 18 Mar 2016\n",
    "\n",
    "This notebook can be downloaded from the github repository, found [here](https://github.com/computationalmodelling/fidimag/blob/master/doc/ipynb/current-driven-domain-wall.ipynb)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Problem specification"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The simulated sample is a 1D nanowire cuboid with $L=1000 \\,\\text{nm}$ length with finite difference discretisation $d_{x} = d_{y} = d_{z} = 2 \\,\\text{nm}$.\n",
    "\n",
    "The material parameters (similar to permalloy) are:\n",
    "    \n",
    "- exchange energy constant $A = 1.3 \\times 10^{-11} \\,\\text{J/m}$,\n",
    "- magnetisation saturation $M_\\text{s} = 8.6 \\times 10^{5} \\,\\text{A/m}$,\n",
    "- uniaxial anisotropy constant $K=5 \\times 10^{4} \\,\\text{J/m}^{3}$ with $(0, 0, 1)$ easy-axis,\n",
    "- Gilbert damping $\\alpha = 0.5$.\n",
    "\n",
    "After the system is relaxed to a domain wall, a spin-polarised current with $J=1 \\times 10^{12} \\text{A/m}^{2}$ density is applied in the positive x direction $(1, 0, 0)$."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Simulation functions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "from fidimag.micro import Sim, UniformExchange, UniaxialAnisotropy\n",
    "from fidimag.common import CuboidMesh\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We start by defining parameters and a function for initialising the system so that it relaxes to the domain wall."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "Ms = 8.6e5  # magnetisation saturation (A/m)\n",
    "A = 1.3e-11  # exchange energy constant (J/m)\n",
    "alpha = 0.5  # Gilbert damping\n",
    "gamma = 2.211e5  # gyromagnetic ration (m/As)\n",
    "K = 5e4  # uniaxial anisotropy constant (J/m**3)\n",
    "J = 1e12  # spin-polarised current density (A/m**2)\n",
    "beta = 1  # STT parameter\n",
    "\n",
    "def init_m(pos):\n",
    "    x = pos[0]\n",
    "\n",
    "    if x < 200:\n",
    "        return (1, 0, 0)\n",
    "    elif 200 <= x < 300:\n",
    "        return (0, 1, 1)\n",
    "    else:\n",
    "        return (-1, 0, 0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Using this function, we create a new function which relaxes the system to its equilibrium (domain wall) state according to the problem specification."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "def relax_system(mesh):\n",
    "    # Create a simulation object.\n",
    "    sim = Sim(mesh)\n",
    "\n",
    "    # Set simulation parameters.\n",
    "    sim.driver.set_tols(rtol=1e-8, atol=1e-8)\n",
    "    sim.driver.alpha = alpha\n",
    "    sim.driver.gamma = gamma\n",
    "    sim.Ms = Ms\n",
    "\n",
    "    # Add energies to the system.\n",
    "    sim.add(UniformExchange(A=A))\n",
    "    sim.add(UniaxialAnisotropy(K))\n",
    "\n",
    "    # Initialise the system.\n",
    "    sim.set_m(init_m)\n",
    "    \n",
    "    # Relax the system and save the state in m0.npy\n",
    "    sim.driver.relax(dt=1e-14, stopping_dmdt=0.01, max_steps=5000,\n",
    "                     save_m_steps=None, save_vtk_steps=None)\n",
    "\n",
    "    np.save('m0.npy', sim.spin)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "A plot of the system's magnetisation can be created using the following convenience function."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot_magnetisation(components):\n",
    "    plt.figure(figsize=(8, 6))\n",
    "\n",
    "    comp = {'mx': 0, 'my': 1, 'mz': 2}\n",
    "\n",
    "    for element in components:\n",
    "        data = np.load(element[0])\n",
    "        data.shape = (-1, 3)\n",
    "\n",
    "        mc = data[:, comp[element[1]]]\n",
    "        \n",
    "        # The label is the component and the file name\n",
    "        plt.plot(mc, label=element[1])\n",
    "\n",
    "    plt.legend()\n",
    "    plt.xlabel('x (nm)')\n",
    "    plt.ylabel('mx, my')\n",
    "    plt.grid()\n",
    "    plt.ylim([-1.05, 1.05])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Finally, we create a function for driving a domain wall using the spin-polarised current. All *npy* and *vtk* files are saved in the **\\{simulation_name\\}_npys** and **\\{simulation_name\\}_vtks** folders, respectively."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "def excite_system(mesh, time, snapshots):\n",
    "    # Specify the stt dynamics in the simulation\n",
    "    sim = Sim(mesh, name='dyn', driver='llg_stt')\n",
    "    \n",
    "    # Set the simulation parameters\n",
    "    sim.driver.set_tols(rtol=1e-12, atol=1e-14)\n",
    "    sim.driver.alpha = alpha\n",
    "    sim.driver.gamma = gamma\n",
    "    sim.Ms = Ms\n",
    "\n",
    "    # Add energies to the system.\n",
    "    sim.add(UniformExchange(A=A))\n",
    "    sim.add(UniaxialAnisotropy(K))\n",
    "\n",
    "    # Load the initial state from the npy file saved in the realxation stage.\n",
    "    sim.set_m(np.load('m0.npy'))\n",
    "\n",
    "    # Set the spin-polarised current in the x direction.\n",
    "    sim.driver.jx = J\n",
    "    sim.driver.beta = beta\n",
    "\n",
    "    # The simulation will run for x ns and save\n",
    "    # 'snaps' snapshots of the system in the process\n",
    "    ts = np.linspace(0, time, snapshots)\n",
    "\n",
    "    for t in ts:\n",
    "        sim.driver.run_until(t)\n",
    "        sim.save_vtk()\n",
    "        sim.save_m()\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Simulation"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Before we run a simulation using previously defined functions, a finite difference mesh must be created."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "L = 2000  # nm\n",
    "dx = dy = dz = 2  # nm\n",
    "\n",
    "mesh = CuboidMesh(nx=int(L/dx), ny=1, nz=1, dx=dx, dy=dy, dz=dz, unit_length=1e-9)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now, the system is relaxed and the domain wall equilibrium state is obtained, saved, and later used in the next stage."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "#1    t=1e-14    dt=1e-14 max_dmdt=6.47e+04\n",
      "#2    t=2e-14    dt=1e-14 max_dmdt=6.48e+04\n",
      "#3    t=3.05e-14 dt=1.05e-14 max_dmdt=6.48e+04\n",
      "#4    t=4.09e-14 dt=1.05e-14 max_dmdt=6.47e+04\n",
      "#5    t=5.14e-14 dt=1.05e-14 max_dmdt=6.47e+04\n",
      "#6    t=6.19e-14 dt=1.05e-14 max_dmdt=6.45e+04\n",
      "#7    t=7.23e-14 dt=1.05e-14 max_dmdt=6.43e+04\n",
      "#8    t=8.28e-14 dt=1.05e-14 max_dmdt=6.41e+04\n",
      "#9    t=9.33e-14 dt=1.05e-14 max_dmdt=6.38e+04\n",
      "#10   t=1.04e-13 dt=1.05e-14 max_dmdt=6.35e+04\n",
      "#11   t=1.23e-13 dt=1.91e-14 max_dmdt=6.29e+04\n",
      "#12   t=1.42e-13 dt=1.91e-14 max_dmdt=6.21e+04\n",
      "#13   t=1.61e-13 dt=1.91e-14 max_dmdt=6.11e+04\n",
      "#14   t=1.8e-13  dt=1.91e-14 max_dmdt=6e+04\n",
      "#15   t=1.99e-13 dt=1.91e-14 max_dmdt=5.88e+04\n",
      "#16   t=2.18e-13 dt=1.91e-14 max_dmdt=5.74e+04\n",
      "#17   t=2.37e-13 dt=1.91e-14 max_dmdt=5.6e+04\n",
      "#18   t=2.56e-13 dt=1.91e-14 max_dmdt=5.44e+04\n",
      "#19   t=2.75e-13 dt=1.91e-14 max_dmdt=5.28e+04\n",
      "#20   t=3.09e-13 dt=3.37e-14 max_dmdt=5.03e+04\n",
      "#21   t=3.43e-13 dt=3.37e-14 max_dmdt=4.7e+04\n",
      "#22   t=3.77e-13 dt=3.37e-14 max_dmdt=4.59e+04\n",
      "#23   t=4.01e-13 dt=2.47e-14 max_dmdt=4.51e+04\n",
      "#24   t=4.26e-13 dt=2.47e-14 max_dmdt=4.42e+04\n",
      "#25   t=4.51e-13 dt=2.47e-14 max_dmdt=4.31e+04\n",
      "#26   t=4.75e-13 dt=2.47e-14 max_dmdt=4.2e+04\n",
      "#27   t=5e-13    dt=2.47e-14 max_dmdt=4.07e+04\n",
      "#28   t=5.25e-13 dt=2.47e-14 max_dmdt=3.92e+04\n",
      "#29   t=5.49e-13 dt=2.47e-14 max_dmdt=3.77e+04\n",
      "#30   t=5.74e-13 dt=2.47e-14 max_dmdt=3.61e+04\n",
      "#31   t=5.99e-13 dt=2.47e-14 max_dmdt=3.45e+04\n",
      "#32   t=6.37e-13 dt=3.85e-14 max_dmdt=3.23e+04\n",
      "#33   t=6.76e-13 dt=3.85e-14 max_dmdt=2.96e+04\n",
      "#34   t=7.14e-13 dt=3.85e-14 max_dmdt=2.75e+04\n",
      "#35   t=7.52e-13 dt=3.85e-14 max_dmdt=2.64e+04\n",
      "#36   t=7.91e-13 dt=3.85e-14 max_dmdt=2.52e+04\n",
      "#37   t=8.29e-13 dt=3.85e-14 max_dmdt=2.46e+04\n",
      "#38   t=8.68e-13 dt=3.85e-14 max_dmdt=2.4e+04\n",
      "#39   t=9.06e-13 dt=3.85e-14 max_dmdt=2.32e+04\n",
      "#40   t=9.45e-13 dt=3.85e-14 max_dmdt=2.23e+04\n",
      "#41   t=9.83e-13 dt=3.85e-14 max_dmdt=2.15e+04\n",
      "#42   t=1.02e-12 dt=3.85e-14 max_dmdt=2.15e+04\n",
      "#43   t=1.06e-12 dt=3.85e-14 max_dmdt=2.13e+04\n",
      "#44   t=1.1e-12  dt=3.85e-14 max_dmdt=2.1e+04\n",
      "#45   t=1.14e-12 dt=3.85e-14 max_dmdt=2.06e+04\n",
      "#46   t=1.18e-12 dt=3.85e-14 max_dmdt=2.01e+04\n",
      "#47   t=1.21e-12 dt=3.85e-14 max_dmdt=1.94e+04\n",
      "#48   t=1.25e-12 dt=3.85e-14 max_dmdt=1.87e+04\n",
      "#49   t=1.29e-12 dt=3.85e-14 max_dmdt=1.79e+04\n",
      "#50   t=1.33e-12 dt=3.85e-14 max_dmdt=1.71e+04\n",
      "#51   t=1.37e-12 dt=3.85e-14 max_dmdt=1.63e+04\n",
      "#52   t=1.41e-12 dt=3.85e-14 max_dmdt=1.54e+04\n",
      "#53   t=1.44e-12 dt=3.85e-14 max_dmdt=1.49e+04\n",
      "#54   t=1.48e-12 dt=3.85e-14 max_dmdt=1.47e+04\n",
      "#55   t=1.52e-12 dt=3.85e-14 max_dmdt=1.44e+04\n",
      "#56   t=1.56e-12 dt=3.85e-14 max_dmdt=1.41e+04\n",
      "#57   t=1.6e-12  dt=3.85e-14 max_dmdt=1.37e+04\n",
      "#58   t=1.64e-12 dt=3.85e-14 max_dmdt=1.33e+04\n",
      "#59   t=1.68e-12 dt=3.85e-14 max_dmdt=1.28e+04\n",
      "#60   t=1.71e-12 dt=3.85e-14 max_dmdt=1.24e+04\n",
      "#61   t=1.75e-12 dt=3.85e-14 max_dmdt=1.2e+04\n",
      "#62   t=1.79e-12 dt=3.85e-14 max_dmdt=1.19e+04\n",
      "#63   t=1.83e-12 dt=3.85e-14 max_dmdt=1.17e+04\n",
      "#64   t=1.87e-12 dt=3.85e-14 max_dmdt=1.16e+04\n",
      "#65   t=1.91e-12 dt=3.85e-14 max_dmdt=1.13e+04\n",
      "#66   t=1.94e-12 dt=3.85e-14 max_dmdt=1.11e+04\n",
      "#67   t=1.98e-12 dt=3.85e-14 max_dmdt=1.09e+04\n",
      "#68   t=2.02e-12 dt=3.85e-14 max_dmdt=1.06e+04\n",
      "#69   t=2.06e-12 dt=3.85e-14 max_dmdt=1.03e+04\n",
      "#70   t=2.1e-12  dt=3.85e-14 max_dmdt=1e+04\n",
      "#71   t=2.14e-12 dt=3.85e-14 max_dmdt=9.71e+03\n",
      "#72   t=2.18e-12 dt=3.85e-14 max_dmdt=9.49e+03\n",
      "#73   t=2.21e-12 dt=3.85e-14 max_dmdt=9.38e+03\n",
      "#74   t=2.25e-12 dt=3.85e-14 max_dmdt=9.25e+03\n",
      "#75   t=2.29e-12 dt=3.85e-14 max_dmdt=9.09e+03\n",
      "#76   t=2.33e-12 dt=3.85e-14 max_dmdt=8.92e+03\n",
      "#77   t=2.37e-12 dt=3.85e-14 max_dmdt=8.74e+03\n",
      "#78   t=2.41e-12 dt=3.85e-14 max_dmdt=8.54e+03\n",
      "#79   t=2.44e-12 dt=3.85e-14 max_dmdt=8.34e+03\n",
      "#80   t=2.48e-12 dt=3.85e-14 max_dmdt=8.13e+03\n",
      "#81   t=2.54e-12 dt=5.79e-14 max_dmdt=7.86e+03\n",
      "#82   t=2.6e-12  dt=5.79e-14 max_dmdt=7.53e+03\n",
      "#83   t=2.66e-12 dt=5.79e-14 max_dmdt=7.21e+03\n",
      "#84   t=2.71e-12 dt=5.79e-14 max_dmdt=6.9e+03\n",
      "#85   t=2.77e-12 dt=5.79e-14 max_dmdt=6.59e+03\n",
      "#86   t=2.83e-12 dt=5.79e-14 max_dmdt=6.46e+03\n",
      "#87   t=2.89e-12 dt=5.79e-14 max_dmdt=6.33e+03\n",
      "#88   t=2.95e-12 dt=5.79e-14 max_dmdt=6.29e+03\n",
      "#89   t=3e-12    dt=5.79e-14 max_dmdt=6.23e+03\n",
      "#90   t=3.06e-12 dt=5.79e-14 max_dmdt=6.15e+03\n",
      "#91   t=3.12e-12 dt=5.79e-14 max_dmdt=6.06e+03\n",
      "#92   t=3.18e-12 dt=5.79e-14 max_dmdt=5.96e+03\n",
      "#93   t=3.24e-12 dt=5.79e-14 max_dmdt=5.85e+03\n",
      "#94   t=3.29e-12 dt=5.79e-14 max_dmdt=5.73e+03\n",
      "#95   t=3.35e-12 dt=5.79e-14 max_dmdt=5.6e+03\n",
      "#96   t=3.41e-12 dt=5.79e-14 max_dmdt=5.48e+03\n",
      "#97   t=3.47e-12 dt=5.79e-14 max_dmdt=5.36e+03\n",
      "#98   t=3.52e-12 dt=5.79e-14 max_dmdt=5.23e+03\n",
      "#99   t=3.58e-12 dt=5.79e-14 max_dmdt=5.11e+03\n",
      "#100  t=3.64e-12 dt=5.79e-14 max_dmdt=4.99e+03\n",
      "#101  t=3.7e-12  dt=5.79e-14 max_dmdt=4.87e+03\n",
      "#102  t=3.76e-12 dt=5.79e-14 max_dmdt=4.76e+03\n",
      "#103  t=3.81e-12 dt=5.79e-14 max_dmdt=4.65e+03\n",
      "#104  t=3.87e-12 dt=5.79e-14 max_dmdt=4.54e+03\n",
      "#105  t=3.93e-12 dt=5.79e-14 max_dmdt=4.45e+03\n",
      "#106  t=3.99e-12 dt=5.79e-14 max_dmdt=4.38e+03\n",
      "#107  t=4.08e-12 dt=8.74e-14 max_dmdt=4.28e+03\n",
      "#108  t=4.16e-12 dt=8.74e-14 max_dmdt=4.17e+03\n",
      "#109  t=4.25e-12 dt=8.74e-14 max_dmdt=4.08e+03\n",
      "#110  t=4.34e-12 dt=8.74e-14 max_dmdt=4.02e+03\n",
      "#111  t=4.43e-12 dt=8.74e-14 max_dmdt=3.96e+03\n",
      "#112  t=4.51e-12 dt=8.74e-14 max_dmdt=3.89e+03\n",
      "#113  t=4.6e-12  dt=8.74e-14 max_dmdt=3.82e+03\n",
      "#114  t=4.69e-12 dt=8.74e-14 max_dmdt=3.75e+03\n",
      "#115  t=4.77e-12 dt=8.74e-14 max_dmdt=3.68e+03\n",
      "#116  t=4.86e-12 dt=8.74e-14 max_dmdt=3.61e+03\n",
      "#117  t=4.95e-12 dt=8.74e-14 max_dmdt=3.53e+03\n",
      "#118  t=5.04e-12 dt=8.74e-14 max_dmdt=3.46e+03\n",
      "#119  t=5.12e-12 dt=8.74e-14 max_dmdt=3.38e+03\n",
      "#120  t=5.21e-12 dt=8.74e-14 max_dmdt=3.31e+03\n",
      "#121  t=5.3e-12  dt=8.74e-14 max_dmdt=3.25e+03\n",
      "#122  t=5.39e-12 dt=8.74e-14 max_dmdt=3.2e+03\n",
      "#123  t=5.47e-12 dt=8.74e-14 max_dmdt=3.15e+03\n",
      "#124  t=5.56e-12 dt=8.74e-14 max_dmdt=3.1e+03\n",
      "#125  t=5.65e-12 dt=8.74e-14 max_dmdt=3.05e+03\n",
      "#126  t=5.74e-12 dt=8.74e-14 max_dmdt=3e+03\n",
      "#127  t=5.82e-12 dt=8.74e-14 max_dmdt=2.96e+03\n",
      "#128  t=5.91e-12 dt=8.74e-14 max_dmdt=2.91e+03\n",
      "#129  t=6e-12    dt=8.74e-14 max_dmdt=2.86e+03\n",
      "#130  t=6.09e-12 dt=8.74e-14 max_dmdt=2.82e+03\n",
      "#131  t=6.17e-12 dt=8.74e-14 max_dmdt=2.78e+03\n",
      "#132  t=6.26e-12 dt=8.74e-14 max_dmdt=2.73e+03\n",
      "#133  t=6.35e-12 dt=8.74e-14 max_dmdt=2.69e+03\n",
      "#134  t=6.44e-12 dt=8.74e-14 max_dmdt=2.65e+03\n",
      "#135  t=6.52e-12 dt=8.74e-14 max_dmdt=2.61e+03\n",
      "#136  t=6.61e-12 dt=8.74e-14 max_dmdt=2.57e+03\n",
      "#137  t=6.7e-12  dt=8.74e-14 max_dmdt=2.53e+03\n",
      "#138  t=6.78e-12 dt=8.74e-14 max_dmdt=2.49e+03\n",
      "#139  t=6.87e-12 dt=8.74e-14 max_dmdt=2.45e+03\n",
      "#140  t=6.96e-12 dt=8.74e-14 max_dmdt=2.42e+03\n",
      "#141  t=7.05e-12 dt=8.74e-14 max_dmdt=2.38e+03\n",
      "#142  t=7.13e-12 dt=8.74e-14 max_dmdt=2.35e+03\n",
      "#143  t=7.22e-12 dt=8.74e-14 max_dmdt=2.31e+03\n",
      "#144  t=7.31e-12 dt=8.74e-14 max_dmdt=2.29e+03\n",
      "#145  t=7.4e-12  dt=8.74e-14 max_dmdt=2.27e+03\n",
      "#146  t=7.48e-12 dt=8.74e-14 max_dmdt=2.25e+03\n",
      "#147  t=7.57e-12 dt=8.74e-14 max_dmdt=2.23e+03\n",
      "#148  t=7.66e-12 dt=8.74e-14 max_dmdt=2.21e+03\n",
      "#149  t=7.75e-12 dt=8.74e-14 max_dmdt=2.18e+03\n",
      "#150  t=7.83e-12 dt=8.74e-14 max_dmdt=2.16e+03\n",
      "#151  t=7.92e-12 dt=8.74e-14 max_dmdt=2.14e+03\n",
      "#152  t=8.01e-12 dt=8.74e-14 max_dmdt=2.12e+03\n",
      "#153  t=8.1e-12  dt=8.74e-14 max_dmdt=2.1e+03\n",
      "#154  t=8.18e-12 dt=8.74e-14 max_dmdt=2.08e+03\n",
      "#155  t=8.27e-12 dt=8.74e-14 max_dmdt=2.06e+03\n",
      "#156  t=8.36e-12 dt=8.74e-14 max_dmdt=2.04e+03\n",
      "#157  t=8.45e-12 dt=8.74e-14 max_dmdt=2.02e+03\n",
      "#158  t=8.53e-12 dt=8.74e-14 max_dmdt=2e+03\n",
      "#159  t=8.62e-12 dt=8.74e-14 max_dmdt=1.98e+03\n",
      "#160  t=8.71e-12 dt=8.74e-14 max_dmdt=1.96e+03\n",
      "#161  t=8.79e-12 dt=8.74e-14 max_dmdt=1.94e+03\n",
      "#162  t=8.88e-12 dt=8.74e-14 max_dmdt=1.92e+03\n",
      "#163  t=8.97e-12 dt=8.74e-14 max_dmdt=1.9e+03\n",
      "#164  t=9.06e-12 dt=8.74e-14 max_dmdt=1.88e+03\n",
      "#165  t=9.14e-12 dt=8.74e-14 max_dmdt=1.86e+03\n",
      "#166  t=9.28e-12 dt=1.32e-13 max_dmdt=1.84e+03\n",
      "#167  t=9.41e-12 dt=1.32e-13 max_dmdt=1.81e+03\n",
      "#168  t=9.54e-12 dt=1.32e-13 max_dmdt=1.78e+03\n",
      "#169  t=9.67e-12 dt=1.32e-13 max_dmdt=1.76e+03\n",
      "#170  t=9.81e-12 dt=1.32e-13 max_dmdt=1.73e+03\n",
      "#171  t=1e-11    dt=2.01e-13 max_dmdt=1.7e+03\n",
      "#172  t=1.02e-11 dt=2.01e-13 max_dmdt=1.66e+03\n",
      "#173  t=1.04e-11 dt=2.01e-13 max_dmdt=1.63e+03\n",
      "#174  t=1.06e-11 dt=2.01e-13 max_dmdt=1.59e+03\n",
      "#175  t=1.08e-11 dt=2.01e-13 max_dmdt=1.57e+03\n",
      "#176  t=1.1e-11  dt=2.01e-13 max_dmdt=1.54e+03\n",
      "#177  t=1.12e-11 dt=2.01e-13 max_dmdt=1.52e+03\n",
      "#178  t=1.14e-11 dt=2.01e-13 max_dmdt=1.5e+03\n",
      "#179  t=1.16e-11 dt=2.01e-13 max_dmdt=1.47e+03\n",
      "#180  t=1.18e-11 dt=2.01e-13 max_dmdt=1.45e+03\n",
      "#181  t=1.2e-11  dt=2.01e-13 max_dmdt=1.43e+03\n",
      "#182  t=1.22e-11 dt=2.01e-13 max_dmdt=1.41e+03\n",
      "#183  t=1.24e-11 dt=2.01e-13 max_dmdt=1.38e+03\n",
      "#184  t=1.26e-11 dt=2.01e-13 max_dmdt=1.36e+03\n",
      "#185  t=1.28e-11 dt=2.01e-13 max_dmdt=1.34e+03\n",
      "#186  t=1.3e-11  dt=2.01e-13 max_dmdt=1.32e+03\n",
      "#187  t=1.32e-11 dt=2.01e-13 max_dmdt=1.3e+03\n",
      "#188  t=1.34e-11 dt=2.01e-13 max_dmdt=1.28e+03\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "#189  t=1.36e-11 dt=2.01e-13 max_dmdt=1.26e+03\n",
      "#190  t=1.38e-11 dt=2.01e-13 max_dmdt=1.24e+03\n",
      "#191  t=1.4e-11  dt=2.01e-13 max_dmdt=1.22e+03\n",
      "#192  t=1.42e-11 dt=2.01e-13 max_dmdt=1.21e+03\n",
      "#193  t=1.44e-11 dt=2.01e-13 max_dmdt=1.19e+03\n",
      "#194  t=1.46e-11 dt=2.01e-13 max_dmdt=1.17e+03\n",
      "#195  t=1.48e-11 dt=2.01e-13 max_dmdt=1.15e+03\n",
      "#196  t=1.51e-11 dt=3.05e-13 max_dmdt=1.14e+03\n",
      "#197  t=1.54e-11 dt=3.05e-13 max_dmdt=1.12e+03\n",
      "#198  t=1.58e-11 dt=3.05e-13 max_dmdt=1.1e+03\n",
      "#199  t=1.61e-11 dt=3.05e-13 max_dmdt=1.08e+03\n",
      "#200  t=1.64e-11 dt=3.05e-13 max_dmdt=1.06e+03\n",
      "#201  t=1.67e-11 dt=3.05e-13 max_dmdt=1.05e+03\n",
      "#202  t=1.7e-11  dt=3.05e-13 max_dmdt=1.03e+03\n",
      "#203  t=1.73e-11 dt=3.05e-13 max_dmdt=1.01e+03\n",
      "#204  t=1.76e-11 dt=3.05e-13 max_dmdt=997\n",
      "#205  t=1.79e-11 dt=3.05e-13 max_dmdt=980\n",
      "#206  t=1.82e-11 dt=3.05e-13 max_dmdt=965\n",
      "#207  t=1.85e-11 dt=3.05e-13 max_dmdt=949\n",
      "#208  t=1.88e-11 dt=3.05e-13 max_dmdt=934\n",
      "#209  t=1.91e-11 dt=3.05e-13 max_dmdt=920\n",
      "#210  t=1.94e-11 dt=3.05e-13 max_dmdt=905\n",
      "#211  t=1.97e-11 dt=3.05e-13 max_dmdt=891\n",
      "#212  t=2e-11    dt=3.05e-13 max_dmdt=878\n",
      "#213  t=2.03e-11 dt=3.05e-13 max_dmdt=867\n",
      "#214  t=2.06e-11 dt=3.05e-13 max_dmdt=857\n",
      "#215  t=2.09e-11 dt=3.05e-13 max_dmdt=847\n",
      "#216  t=2.12e-11 dt=3.05e-13 max_dmdt=837\n",
      "#217  t=2.15e-11 dt=3.05e-13 max_dmdt=827\n",
      "#218  t=2.18e-11 dt=3.05e-13 max_dmdt=817\n",
      "#219  t=2.22e-11 dt=3.05e-13 max_dmdt=807\n",
      "#220  t=2.25e-11 dt=3.05e-13 max_dmdt=798\n",
      "#221  t=2.28e-11 dt=3.05e-13 max_dmdt=788\n",
      "#222  t=2.31e-11 dt=3.05e-13 max_dmdt=779\n",
      "#223  t=2.34e-11 dt=3.05e-13 max_dmdt=769\n",
      "#224  t=2.37e-11 dt=3.05e-13 max_dmdt=760\n",
      "#225  t=2.4e-11  dt=3.05e-13 max_dmdt=751\n",
      "#226  t=2.43e-11 dt=3.05e-13 max_dmdt=742\n",
      "#227  t=2.46e-11 dt=3.05e-13 max_dmdt=733\n",
      "#228  t=2.49e-11 dt=3.05e-13 max_dmdt=724\n",
      "#229  t=2.52e-11 dt=3.05e-13 max_dmdt=716\n",
      "#230  t=2.55e-11 dt=3.05e-13 max_dmdt=707\n",
      "#231  t=2.58e-11 dt=3.05e-13 max_dmdt=699\n",
      "#232  t=2.61e-11 dt=3.05e-13 max_dmdt=691\n",
      "#233  t=2.64e-11 dt=3.05e-13 max_dmdt=684\n",
      "#234  t=2.67e-11 dt=3.05e-13 max_dmdt=677\n",
      "#235  t=2.7e-11  dt=3.05e-13 max_dmdt=672\n",
      "#236  t=2.73e-11 dt=3.05e-13 max_dmdt=667\n",
      "#237  t=2.76e-11 dt=3.05e-13 max_dmdt=662\n",
      "#238  t=2.79e-11 dt=3.05e-13 max_dmdt=657\n",
      "#239  t=2.82e-11 dt=3.05e-13 max_dmdt=653\n",
      "#240  t=2.85e-11 dt=3.05e-13 max_dmdt=649\n",
      "#241  t=2.89e-11 dt=3.05e-13 max_dmdt=645\n",
      "#242  t=2.92e-11 dt=3.05e-13 max_dmdt=641\n",
      "#243  t=2.95e-11 dt=3.05e-13 max_dmdt=637\n",
      "#244  t=2.98e-11 dt=3.05e-13 max_dmdt=634\n",
      "#245  t=3.01e-11 dt=3.05e-13 max_dmdt=630\n",
      "#246  t=3.04e-11 dt=3.05e-13 max_dmdt=627\n",
      "#247  t=3.07e-11 dt=3.05e-13 max_dmdt=623\n",
      "#248  t=3.1e-11  dt=3.05e-13 max_dmdt=620\n",
      "#249  t=3.13e-11 dt=3.05e-13 max_dmdt=616\n",
      "#250  t=3.16e-11 dt=3.05e-13 max_dmdt=613\n",
      "#251  t=3.19e-11 dt=3.05e-13 max_dmdt=609\n",
      "#252  t=3.22e-11 dt=3.05e-13 max_dmdt=606\n",
      "#253  t=3.25e-11 dt=3.05e-13 max_dmdt=602\n",
      "#254  t=3.28e-11 dt=3.05e-13 max_dmdt=599\n",
      "#255  t=3.31e-11 dt=3.05e-13 max_dmdt=595\n",
      "#256  t=3.34e-11 dt=3.05e-13 max_dmdt=591\n",
      "#257  t=3.37e-11 dt=3.05e-13 max_dmdt=587\n",
      "#258  t=3.4e-11  dt=3.05e-13 max_dmdt=583\n",
      "#259  t=3.43e-11 dt=3.05e-13 max_dmdt=579\n",
      "#260  t=3.46e-11 dt=3.05e-13 max_dmdt=575\n",
      "#261  t=3.49e-11 dt=3.05e-13 max_dmdt=570\n",
      "#262  t=3.53e-11 dt=3.05e-13 max_dmdt=566\n",
      "#263  t=3.56e-11 dt=3.05e-13 max_dmdt=561\n",
      "#264  t=3.59e-11 dt=3.05e-13 max_dmdt=557\n",
      "#265  t=3.62e-11 dt=3.05e-13 max_dmdt=552\n",
      "#266  t=3.65e-11 dt=3.05e-13 max_dmdt=547\n",
      "#267  t=3.68e-11 dt=3.05e-13 max_dmdt=542\n",
      "#268  t=3.71e-11 dt=3.05e-13 max_dmdt=537\n",
      "#269  t=3.74e-11 dt=3.05e-13 max_dmdt=532\n",
      "#270  t=3.77e-11 dt=3.05e-13 max_dmdt=526\n",
      "#271  t=3.8e-11  dt=3.05e-13 max_dmdt=521\n",
      "#272  t=3.83e-11 dt=3.05e-13 max_dmdt=516\n",
      "#273  t=3.86e-11 dt=3.05e-13 max_dmdt=510\n",
      "#274  t=3.89e-11 dt=3.05e-13 max_dmdt=505\n",
      "#275  t=3.92e-11 dt=3.05e-13 max_dmdt=499\n",
      "#276  t=3.95e-11 dt=3.05e-13 max_dmdt=493\n",
      "#277  t=3.98e-11 dt=3.05e-13 max_dmdt=488\n",
      "#278  t=4.01e-11 dt=3.05e-13 max_dmdt=482\n",
      "#279  t=4.04e-11 dt=3.05e-13 max_dmdt=477\n",
      "#280  t=4.07e-11 dt=3.05e-13 max_dmdt=471\n",
      "#281  t=4.1e-11  dt=3.05e-13 max_dmdt=466\n",
      "#282  t=4.13e-11 dt=3.05e-13 max_dmdt=460\n",
      "#283  t=4.16e-11 dt=3.05e-13 max_dmdt=455\n",
      "#284  t=4.2e-11  dt=3.05e-13 max_dmdt=454\n",
      "#285  t=4.23e-11 dt=3.05e-13 max_dmdt=454\n",
      "#286  t=4.26e-11 dt=3.05e-13 max_dmdt=453\n",
      "#287  t=4.29e-11 dt=3.05e-13 max_dmdt=452\n",
      "#288  t=4.32e-11 dt=3.05e-13 max_dmdt=452\n",
      "#289  t=4.35e-11 dt=3.05e-13 max_dmdt=451\n",
      "#290  t=4.38e-11 dt=3.05e-13 max_dmdt=450\n",
      "#291  t=4.41e-11 dt=3.05e-13 max_dmdt=449\n",
      "#292  t=4.44e-11 dt=3.05e-13 max_dmdt=448\n",
      "#293  t=4.47e-11 dt=3.05e-13 max_dmdt=446\n",
      "#294  t=4.5e-11  dt=3.05e-13 max_dmdt=445\n",
      "#295  t=4.53e-11 dt=3.05e-13 max_dmdt=444\n",
      "#296  t=4.56e-11 dt=3.05e-13 max_dmdt=442\n",
      "#297  t=4.59e-11 dt=3.05e-13 max_dmdt=441\n",
      "#298  t=4.62e-11 dt=3.05e-13 max_dmdt=439\n",
      "#299  t=4.65e-11 dt=3.05e-13 max_dmdt=438\n",
      "#300  t=4.68e-11 dt=3.05e-13 max_dmdt=437\n",
      "#301  t=4.71e-11 dt=3.05e-13 max_dmdt=436\n",
      "#302  t=4.74e-11 dt=3.05e-13 max_dmdt=435\n",
      "#303  t=4.77e-11 dt=3.05e-13 max_dmdt=434\n",
      "#304  t=4.8e-11  dt=3.05e-13 max_dmdt=433\n",
      "#305  t=4.84e-11 dt=3.05e-13 max_dmdt=432\n",
      "#306  t=4.87e-11 dt=3.05e-13 max_dmdt=431\n",
      "#307  t=4.9e-11  dt=3.05e-13 max_dmdt=430\n",
      "#308  t=4.93e-11 dt=3.05e-13 max_dmdt=429\n",
      "#309  t=4.96e-11 dt=3.05e-13 max_dmdt=428\n",
      "#310  t=4.98e-11 dt=2.25e-13 max_dmdt=427\n",
      "#311  t=5e-11    dt=2.25e-13 max_dmdt=426\n",
      "#312  t=5.02e-11 dt=2.25e-13 max_dmdt=425\n",
      "#313  t=5.05e-11 dt=2.25e-13 max_dmdt=424\n",
      "#314  t=5.07e-11 dt=2.25e-13 max_dmdt=423\n",
      "#315  t=5.09e-11 dt=2.25e-13 max_dmdt=423\n",
      "#316  t=5.11e-11 dt=2.25e-13 max_dmdt=423\n",
      "#317  t=5.14e-11 dt=2.25e-13 max_dmdt=422\n",
      "#318  t=5.16e-11 dt=2.25e-13 max_dmdt=422\n",
      "#319  t=5.18e-11 dt=2.25e-13 max_dmdt=421\n",
      "#320  t=5.2e-11  dt=2.25e-13 max_dmdt=421\n",
      "#321  t=5.23e-11 dt=2.25e-13 max_dmdt=420\n",
      "#322  t=5.25e-11 dt=2.25e-13 max_dmdt=420\n",
      "#323  t=5.27e-11 dt=2.25e-13 max_dmdt=419\n",
      "#324  t=5.29e-11 dt=2.25e-13 max_dmdt=419\n",
      "#325  t=5.32e-11 dt=2.25e-13 max_dmdt=418\n",
      "#326  t=5.34e-11 dt=2.25e-13 max_dmdt=418\n",
      "#327  t=5.36e-11 dt=2.25e-13 max_dmdt=417\n",
      "#328  t=5.39e-11 dt=2.25e-13 max_dmdt=416\n",
      "#329  t=5.41e-11 dt=2.25e-13 max_dmdt=416\n",
      "#330  t=5.43e-11 dt=2.25e-13 max_dmdt=417\n",
      "#331  t=5.45e-11 dt=2.25e-13 max_dmdt=419\n",
      "#332  t=5.48e-11 dt=2.25e-13 max_dmdt=422\n",
      "#333  t=5.5e-11  dt=2.25e-13 max_dmdt=425\n",
      "#334  t=5.52e-11 dt=2.25e-13 max_dmdt=427\n",
      "#335  t=5.54e-11 dt=2.25e-13 max_dmdt=430\n",
      "#336  t=5.57e-11 dt=2.25e-13 max_dmdt=432\n",
      "#337  t=5.59e-11 dt=2.25e-13 max_dmdt=435\n",
      "#338  t=5.61e-11 dt=2.25e-13 max_dmdt=437\n",
      "#339  t=5.63e-11 dt=2.25e-13 max_dmdt=439\n",
      "#340  t=5.66e-11 dt=2.25e-13 max_dmdt=442\n",
      "#341  t=5.68e-11 dt=2.25e-13 max_dmdt=444\n",
      "#342  t=5.7e-11  dt=2.25e-13 max_dmdt=446\n",
      "#343  t=5.72e-11 dt=2.25e-13 max_dmdt=449\n",
      "#344  t=5.75e-11 dt=2.25e-13 max_dmdt=451\n",
      "#345  t=5.77e-11 dt=2.25e-13 max_dmdt=453\n",
      "#346  t=5.79e-11 dt=2.25e-13 max_dmdt=455\n",
      "#347  t=5.81e-11 dt=2.25e-13 max_dmdt=457\n",
      "#348  t=5.84e-11 dt=2.25e-13 max_dmdt=459\n",
      "#349  t=5.86e-11 dt=2.25e-13 max_dmdt=461\n",
      "#350  t=5.88e-11 dt=2.25e-13 max_dmdt=463\n",
      "#351  t=5.9e-11  dt=2.25e-13 max_dmdt=465\n",
      "#352  t=5.93e-11 dt=2.25e-13 max_dmdt=467\n",
      "#353  t=5.95e-11 dt=2.25e-13 max_dmdt=468\n",
      "#354  t=5.97e-11 dt=2.25e-13 max_dmdt=470\n",
      "#355  t=5.99e-11 dt=2.25e-13 max_dmdt=472\n",
      "#356  t=6.02e-11 dt=2.25e-13 max_dmdt=473\n",
      "#357  t=6.04e-11 dt=2.25e-13 max_dmdt=475\n",
      "#358  t=6.06e-11 dt=2.25e-13 max_dmdt=477\n",
      "#359  t=6.08e-11 dt=2.25e-13 max_dmdt=478\n",
      "#360  t=6.11e-11 dt=2.25e-13 max_dmdt=480\n",
      "#361  t=6.13e-11 dt=2.25e-13 max_dmdt=481\n",
      "#362  t=6.15e-11 dt=2.25e-13 max_dmdt=483\n",
      "#363  t=6.17e-11 dt=2.25e-13 max_dmdt=484\n",
      "#364  t=6.2e-11  dt=2.25e-13 max_dmdt=486\n",
      "#365  t=6.22e-11 dt=2.25e-13 max_dmdt=487\n",
      "#366  t=6.24e-11 dt=2.25e-13 max_dmdt=488\n",
      "#367  t=6.26e-11 dt=2.25e-13 max_dmdt=490\n",
      "#368  t=6.29e-11 dt=2.25e-13 max_dmdt=491\n",
      "#369  t=6.31e-11 dt=2.25e-13 max_dmdt=492\n",
      "#370  t=6.33e-11 dt=2.25e-13 max_dmdt=494\n",
      "#371  t=6.35e-11 dt=2.25e-13 max_dmdt=495\n",
      "#372  t=6.38e-11 dt=2.25e-13 max_dmdt=496\n",
      "#373  t=6.4e-11  dt=2.25e-13 max_dmdt=498\n",
      "#374  t=6.42e-11 dt=2.25e-13 max_dmdt=499\n",
      "#375  t=6.44e-11 dt=2.25e-13 max_dmdt=500\n",
      "#376  t=6.47e-11 dt=2.25e-13 max_dmdt=501\n",
      "#377  t=6.5e-11  dt=3.43e-13 max_dmdt=502\n",
      "#378  t=6.54e-11 dt=3.43e-13 max_dmdt=504\n",
      "#379  t=6.57e-11 dt=3.43e-13 max_dmdt=506\n",
      "#380  t=6.6e-11  dt=3.43e-13 max_dmdt=507\n",
      "#381  t=6.64e-11 dt=3.43e-13 max_dmdt=509\n",
      "#382  t=6.67e-11 dt=3.43e-13 max_dmdt=510\n",
      "#383  t=6.71e-11 dt=3.43e-13 max_dmdt=511\n",
      "#384  t=6.74e-11 dt=3.43e-13 max_dmdt=513\n",
      "#385  t=6.78e-11 dt=3.43e-13 max_dmdt=514\n",
      "#386  t=6.81e-11 dt=3.43e-13 max_dmdt=515\n",
      "#387  t=6.86e-11 dt=5.15e-13 max_dmdt=517\n",
      "#388  t=6.91e-11 dt=5.15e-13 max_dmdt=518\n",
      "#389  t=6.96e-11 dt=5.15e-13 max_dmdt=520\n",
      "#390  t=7.02e-11 dt=5.15e-13 max_dmdt=521\n",
      "#391  t=7.07e-11 dt=5.15e-13 max_dmdt=523\n",
      "#392  t=7.12e-11 dt=5.15e-13 max_dmdt=524\n",
      "#393  t=7.17e-11 dt=5.15e-13 max_dmdt=525\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "#394  t=7.22e-11 dt=5.15e-13 max_dmdt=526\n",
      "#395  t=7.27e-11 dt=5.15e-13 max_dmdt=527\n",
      "#396  t=7.35e-11 dt=7.88e-13 max_dmdt=528\n",
      "#397  t=7.43e-11 dt=7.88e-13 max_dmdt=529\n",
      "#398  t=7.51e-11 dt=7.88e-13 max_dmdt=530\n",
      "#399  t=7.59e-11 dt=7.88e-13 max_dmdt=531\n",
      "#400  t=7.67e-11 dt=7.88e-13 max_dmdt=532\n",
      "#401  t=7.75e-11 dt=7.88e-13 max_dmdt=532\n",
      "#402  t=7.83e-11 dt=7.88e-13 max_dmdt=533\n",
      "#403  t=7.9e-11  dt=7.88e-13 max_dmdt=533\n",
      "#404  t=7.98e-11 dt=7.88e-13 max_dmdt=533\n",
      "#405  t=8.06e-11 dt=7.88e-13 max_dmdt=533\n",
      "#406  t=8.14e-11 dt=7.88e-13 max_dmdt=533\n",
      "#407  t=8.22e-11 dt=7.88e-13 max_dmdt=533\n",
      "#408  t=8.3e-11  dt=7.88e-13 max_dmdt=532\n",
      "#409  t=8.38e-11 dt=7.88e-13 max_dmdt=532\n",
      "#410  t=8.46e-11 dt=7.88e-13 max_dmdt=531\n",
      "#411  t=8.53e-11 dt=7.88e-13 max_dmdt=531\n",
      "#412  t=8.61e-11 dt=7.88e-13 max_dmdt=530\n",
      "#413  t=8.69e-11 dt=7.88e-13 max_dmdt=529\n",
      "#414  t=8.77e-11 dt=7.88e-13 max_dmdt=529\n",
      "#415  t=8.85e-11 dt=7.88e-13 max_dmdt=528\n",
      "#416  t=8.93e-11 dt=7.88e-13 max_dmdt=527\n",
      "#417  t=9.01e-11 dt=7.88e-13 max_dmdt=526\n",
      "#418  t=9.09e-11 dt=7.88e-13 max_dmdt=524\n",
      "#419  t=9.17e-11 dt=7.88e-13 max_dmdt=523\n",
      "#420  t=9.24e-11 dt=7.88e-13 max_dmdt=522\n",
      "#421  t=9.32e-11 dt=7.88e-13 max_dmdt=521\n",
      "#422  t=9.44e-11 dt=1.21e-12 max_dmdt=519\n",
      "#423  t=9.56e-11 dt=1.21e-12 max_dmdt=517\n",
      "#424  t=9.68e-11 dt=1.21e-12 max_dmdt=514\n",
      "#425  t=9.81e-11 dt=1.21e-12 max_dmdt=512\n",
      "#426  t=9.93e-11 dt=1.21e-12 max_dmdt=509\n",
      "#427  t=1e-10    dt=1.21e-12 max_dmdt=507\n",
      "#428  t=1.02e-10 dt=1.21e-12 max_dmdt=504\n",
      "#429  t=1.03e-10 dt=1.21e-12 max_dmdt=501\n",
      "#430  t=1.04e-10 dt=1.21e-12 max_dmdt=498\n",
      "#431  t=1.05e-10 dt=1.21e-12 max_dmdt=495\n",
      "#432  t=1.06e-10 dt=1.21e-12 max_dmdt=491\n",
      "#433  t=1.08e-10 dt=1.21e-12 max_dmdt=488\n",
      "#434  t=1.09e-10 dt=1.21e-12 max_dmdt=484\n",
      "#435  t=1.1e-10  dt=1.21e-12 max_dmdt=481\n",
      "#436  t=1.11e-10 dt=1.21e-12 max_dmdt=477\n",
      "#437  t=1.13e-10 dt=1.21e-12 max_dmdt=473\n",
      "#438  t=1.14e-10 dt=1.21e-12 max_dmdt=469\n",
      "#439  t=1.15e-10 dt=1.21e-12 max_dmdt=465\n",
      "#440  t=1.16e-10 dt=1.21e-12 max_dmdt=461\n",
      "#441  t=1.17e-10 dt=1.21e-12 max_dmdt=457\n",
      "#442  t=1.19e-10 dt=1.21e-12 max_dmdt=453\n",
      "#443  t=1.2e-10  dt=1.21e-12 max_dmdt=449\n",
      "#444  t=1.21e-10 dt=1.21e-12 max_dmdt=444\n",
      "#445  t=1.22e-10 dt=1.21e-12 max_dmdt=440\n",
      "#446  t=1.23e-10 dt=1.21e-12 max_dmdt=435\n",
      "#447  t=1.25e-10 dt=1.21e-12 max_dmdt=431\n",
      "#448  t=1.26e-10 dt=1.21e-12 max_dmdt=426\n",
      "#449  t=1.27e-10 dt=1.21e-12 max_dmdt=421\n",
      "#450  t=1.28e-10 dt=1.21e-12 max_dmdt=421\n",
      "#451  t=1.29e-10 dt=1.21e-12 max_dmdt=424\n",
      "#452  t=1.31e-10 dt=1.21e-12 max_dmdt=426\n",
      "#453  t=1.32e-10 dt=1.21e-12 max_dmdt=428\n",
      "#454  t=1.33e-10 dt=1.21e-12 max_dmdt=429\n",
      "#455  t=1.34e-10 dt=1.21e-12 max_dmdt=431\n",
      "#456  t=1.35e-10 dt=1.21e-12 max_dmdt=432\n",
      "#457  t=1.37e-10 dt=1.21e-12 max_dmdt=432\n",
      "#458  t=1.38e-10 dt=1.21e-12 max_dmdt=433\n",
      "#459  t=1.39e-10 dt=1.21e-12 max_dmdt=433\n",
      "#460  t=1.4e-10  dt=1.21e-12 max_dmdt=433\n",
      "#461  t=1.41e-10 dt=1.21e-12 max_dmdt=433\n",
      "#462  t=1.43e-10 dt=1.21e-12 max_dmdt=432\n",
      "#463  t=1.44e-10 dt=1.21e-12 max_dmdt=431\n",
      "#464  t=1.45e-10 dt=1.21e-12 max_dmdt=430\n",
      "#465  t=1.46e-10 dt=1.21e-12 max_dmdt=429\n",
      "#466  t=1.48e-10 dt=1.21e-12 max_dmdt=427\n",
      "#467  t=1.49e-10 dt=1.21e-12 max_dmdt=425\n",
      "#468  t=1.5e-10  dt=1.21e-12 max_dmdt=423\n",
      "#469  t=1.51e-10 dt=1.21e-12 max_dmdt=421\n",
      "#470  t=1.52e-10 dt=1.21e-12 max_dmdt=418\n",
      "#471  t=1.54e-10 dt=1.21e-12 max_dmdt=416\n",
      "#472  t=1.55e-10 dt=1.21e-12 max_dmdt=413\n",
      "#473  t=1.56e-10 dt=1.21e-12 max_dmdt=410\n",
      "#474  t=1.57e-10 dt=1.21e-12 max_dmdt=406\n",
      "#475  t=1.58e-10 dt=1.21e-12 max_dmdt=403\n",
      "#476  t=1.6e-10  dt=1.21e-12 max_dmdt=399\n",
      "#477  t=1.61e-10 dt=1.21e-12 max_dmdt=395\n",
      "#478  t=1.62e-10 dt=1.21e-12 max_dmdt=391\n",
      "#479  t=1.63e-10 dt=1.21e-12 max_dmdt=386\n",
      "#480  t=1.64e-10 dt=1.21e-12 max_dmdt=382\n",
      "#481  t=1.66e-10 dt=1.21e-12 max_dmdt=377\n",
      "#482  t=1.67e-10 dt=1.21e-12 max_dmdt=372\n",
      "#483  t=1.68e-10 dt=1.21e-12 max_dmdt=367\n",
      "#484  t=1.69e-10 dt=1.21e-12 max_dmdt=362\n",
      "#485  t=1.7e-10  dt=1.21e-12 max_dmdt=356\n",
      "#486  t=1.72e-10 dt=1.21e-12 max_dmdt=351\n",
      "#487  t=1.73e-10 dt=1.21e-12 max_dmdt=345\n",
      "#488  t=1.74e-10 dt=1.21e-12 max_dmdt=339\n",
      "#489  t=1.75e-10 dt=1.21e-12 max_dmdt=333\n",
      "#490  t=1.77e-10 dt=1.81e-12 max_dmdt=326\n",
      "#491  t=1.79e-10 dt=1.81e-12 max_dmdt=316\n",
      "#492  t=1.81e-10 dt=1.81e-12 max_dmdt=307\n",
      "#493  t=1.82e-10 dt=1.81e-12 max_dmdt=297\n",
      "#494  t=1.84e-10 dt=1.81e-12 max_dmdt=287\n",
      "#495  t=1.86e-10 dt=1.81e-12 max_dmdt=277\n",
      "#496  t=1.88e-10 dt=1.81e-12 max_dmdt=267\n",
      "#497  t=1.9e-10  dt=1.81e-12 max_dmdt=257\n",
      "#498  t=1.92e-10 dt=1.81e-12 max_dmdt=246\n",
      "#499  t=1.93e-10 dt=1.81e-12 max_dmdt=236\n",
      "#500  t=1.95e-10 dt=1.81e-12 max_dmdt=226\n",
      "#501  t=1.97e-10 dt=1.81e-12 max_dmdt=219\n",
      "#502  t=1.99e-10 dt=1.81e-12 max_dmdt=221\n",
      "#503  t=2.01e-10 dt=1.81e-12 max_dmdt=223\n",
      "#504  t=2.02e-10 dt=1.81e-12 max_dmdt=224\n",
      "#505  t=2.04e-10 dt=1.81e-12 max_dmdt=224\n",
      "#506  t=2.06e-10 dt=1.81e-12 max_dmdt=225\n",
      "#507  t=2.08e-10 dt=1.81e-12 max_dmdt=225\n",
      "#508  t=2.1e-10  dt=1.81e-12 max_dmdt=224\n",
      "#509  t=2.11e-10 dt=1.81e-12 max_dmdt=224\n",
      "#510  t=2.13e-10 dt=1.81e-12 max_dmdt=223\n",
      "#511  t=2.15e-10 dt=1.81e-12 max_dmdt=222\n",
      "#512  t=2.17e-10 dt=1.81e-12 max_dmdt=220\n",
      "#513  t=2.19e-10 dt=1.81e-12 max_dmdt=218\n",
      "#514  t=2.21e-10 dt=1.81e-12 max_dmdt=216\n",
      "#515  t=2.22e-10 dt=1.81e-12 max_dmdt=214\n",
      "#516  t=2.24e-10 dt=1.81e-12 max_dmdt=212\n",
      "#517  t=2.26e-10 dt=1.81e-12 max_dmdt=209\n",
      "#518  t=2.28e-10 dt=1.81e-12 max_dmdt=207\n",
      "#519  t=2.3e-10  dt=1.81e-12 max_dmdt=204\n",
      "#520  t=2.31e-10 dt=1.81e-12 max_dmdt=201\n",
      "#521  t=2.33e-10 dt=1.81e-12 max_dmdt=197\n",
      "#522  t=2.35e-10 dt=1.81e-12 max_dmdt=194\n",
      "#523  t=2.37e-10 dt=1.81e-12 max_dmdt=191\n",
      "#524  t=2.39e-10 dt=1.81e-12 max_dmdt=187\n",
      "#525  t=2.4e-10  dt=1.81e-12 max_dmdt=184\n",
      "#526  t=2.42e-10 dt=1.81e-12 max_dmdt=180\n",
      "#527  t=2.44e-10 dt=1.81e-12 max_dmdt=176\n",
      "#528  t=2.46e-10 dt=1.81e-12 max_dmdt=172\n",
      "#529  t=2.48e-10 dt=1.81e-12 max_dmdt=168\n",
      "#530  t=2.49e-10 dt=1.81e-12 max_dmdt=164\n",
      "#531  t=2.51e-10 dt=1.81e-12 max_dmdt=160\n",
      "#532  t=2.53e-10 dt=1.81e-12 max_dmdt=156\n",
      "#533  t=2.55e-10 dt=1.81e-12 max_dmdt=152\n",
      "#534  t=2.57e-10 dt=1.81e-12 max_dmdt=148\n",
      "#535  t=2.59e-10 dt=1.81e-12 max_dmdt=144\n",
      "#536  t=2.6e-10  dt=1.81e-12 max_dmdt=140\n",
      "#537  t=2.62e-10 dt=1.81e-12 max_dmdt=136\n",
      "#538  t=2.64e-10 dt=1.81e-12 max_dmdt=132\n",
      "#539  t=2.67e-10 dt=2.74e-12 max_dmdt=127\n",
      "#540  t=2.69e-10 dt=2.74e-12 max_dmdt=121\n",
      "#541  t=2.72e-10 dt=2.74e-12 max_dmdt=115\n",
      "#542  t=2.75e-10 dt=2.74e-12 max_dmdt=110\n",
      "#543  t=2.78e-10 dt=2.74e-12 max_dmdt=104\n",
      "#544  t=2.8e-10  dt=2.74e-12 max_dmdt=98.3\n",
      "#545  t=2.83e-10 dt=2.74e-12 max_dmdt=92.9\n",
      "#546  t=2.86e-10 dt=2.74e-12 max_dmdt=87.7\n",
      "#547  t=2.89e-10 dt=2.74e-12 max_dmdt=82.8\n",
      "#548  t=2.91e-10 dt=2.74e-12 max_dmdt=78\n",
      "#549  t=2.94e-10 dt=2.74e-12 max_dmdt=77.7\n",
      "#550  t=2.97e-10 dt=2.74e-12 max_dmdt=78.4\n",
      "#551  t=3e-10    dt=2.74e-12 max_dmdt=78.9\n",
      "#552  t=3.02e-10 dt=2.74e-12 max_dmdt=79.1\n",
      "#553  t=3.05e-10 dt=2.74e-12 max_dmdt=79.1\n",
      "#554  t=3.08e-10 dt=2.74e-12 max_dmdt=78.9\n",
      "#555  t=3.11e-10 dt=2.74e-12 max_dmdt=78.5\n",
      "#556  t=3.13e-10 dt=2.74e-12 max_dmdt=77.9\n",
      "#557  t=3.16e-10 dt=2.74e-12 max_dmdt=77.1\n",
      "#558  t=3.19e-10 dt=2.74e-12 max_dmdt=76.2\n",
      "#559  t=3.22e-10 dt=2.74e-12 max_dmdt=75.1\n",
      "#560  t=3.24e-10 dt=2.74e-12 max_dmdt=73.9\n",
      "#561  t=3.27e-10 dt=2.74e-12 max_dmdt=72.5\n",
      "#562  t=3.3e-10  dt=2.74e-12 max_dmdt=71.1\n",
      "#563  t=3.32e-10 dt=2.74e-12 max_dmdt=69.5\n",
      "#564  t=3.35e-10 dt=2.74e-12 max_dmdt=67.8\n",
      "#565  t=3.38e-10 dt=2.74e-12 max_dmdt=66.1\n",
      "#566  t=3.41e-10 dt=2.74e-12 max_dmdt=64.2\n",
      "#567  t=3.43e-10 dt=2.74e-12 max_dmdt=62.3\n",
      "#568  t=3.46e-10 dt=2.74e-12 max_dmdt=60.4\n",
      "#569  t=3.49e-10 dt=2.74e-12 max_dmdt=58.4\n",
      "#570  t=3.52e-10 dt=2.74e-12 max_dmdt=56.4\n",
      "#571  t=3.54e-10 dt=2.74e-12 max_dmdt=54.3\n",
      "#572  t=3.57e-10 dt=2.74e-12 max_dmdt=52.2\n",
      "#573  t=3.6e-10  dt=2.74e-12 max_dmdt=50.1\n",
      "#574  t=3.63e-10 dt=2.74e-12 max_dmdt=48\n",
      "#575  t=3.65e-10 dt=2.74e-12 max_dmdt=45.9\n",
      "#576  t=3.68e-10 dt=2.74e-12 max_dmdt=43.8\n",
      "#577  t=3.71e-10 dt=2.74e-12 max_dmdt=41.7\n",
      "#578  t=3.74e-10 dt=2.74e-12 max_dmdt=39.7\n",
      "#579  t=3.76e-10 dt=2.74e-12 max_dmdt=37.7\n",
      "#580  t=3.79e-10 dt=2.74e-12 max_dmdt=35.8\n",
      "#581  t=3.82e-10 dt=2.74e-12 max_dmdt=34\n",
      "#582  t=3.85e-10 dt=2.74e-12 max_dmdt=32.2\n",
      "#583  t=3.87e-10 dt=2.74e-12 max_dmdt=31.5\n",
      "#584  t=3.9e-10  dt=2.74e-12 max_dmdt=31.8\n",
      "#585  t=3.93e-10 dt=2.74e-12 max_dmdt=32\n",
      "#586  t=3.96e-10 dt=2.74e-12 max_dmdt=32.1\n",
      "#587  t=3.98e-10 dt=2.74e-12 max_dmdt=32.1\n",
      "#588  t=4.01e-10 dt=2.74e-12 max_dmdt=32.1\n",
      "#589  t=4.04e-10 dt=2.74e-12 max_dmdt=32\n",
      "#590  t=4.06e-10 dt=2.74e-12 max_dmdt=31.8\n",
      "#591  t=4.09e-10 dt=2.74e-12 max_dmdt=31.5\n",
      "#592  t=4.12e-10 dt=2.74e-12 max_dmdt=31.2\n",
      "#593  t=4.15e-10 dt=2.74e-12 max_dmdt=30.8\n",
      "#594  t=4.17e-10 dt=2.74e-12 max_dmdt=30.3\n",
      "#595  t=4.2e-10  dt=2.74e-12 max_dmdt=29.8\n",
      "#596  t=4.23e-10 dt=2.74e-12 max_dmdt=29.3\n",
      "#597  t=4.26e-10 dt=2.74e-12 max_dmdt=28.7\n",
      "#598  t=4.28e-10 dt=2.74e-12 max_dmdt=28.1\n",
      "#599  t=4.31e-10 dt=2.74e-12 max_dmdt=27.4\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "#600  t=4.34e-10 dt=2.74e-12 max_dmdt=26.7\n",
      "#601  t=4.37e-10 dt=2.74e-12 max_dmdt=26\n",
      "#602  t=4.39e-10 dt=2.74e-12 max_dmdt=25.2\n",
      "#603  t=4.42e-10 dt=2.74e-12 max_dmdt=24.4\n",
      "#604  t=4.45e-10 dt=2.74e-12 max_dmdt=23.6\n",
      "#605  t=4.48e-10 dt=2.74e-12 max_dmdt=22.8\n",
      "#606  t=4.5e-10  dt=2.74e-12 max_dmdt=21.9\n",
      "#607  t=4.53e-10 dt=2.74e-12 max_dmdt=21.1\n",
      "#608  t=4.57e-10 dt=4.11e-12 max_dmdt=20\n",
      "#609  t=4.61e-10 dt=4.11e-12 max_dmdt=18.8\n",
      "#610  t=4.65e-10 dt=4.11e-12 max_dmdt=17.5\n",
      "#611  t=4.69e-10 dt=4.11e-12 max_dmdt=16.3\n",
      "#612  t=4.74e-10 dt=4.11e-12 max_dmdt=15.1\n",
      "#613  t=4.78e-10 dt=4.11e-12 max_dmdt=14\n",
      "#614  t=4.82e-10 dt=4.11e-12 max_dmdt=13\n",
      "#615  t=4.86e-10 dt=4.11e-12 max_dmdt=12.9\n",
      "#616  t=4.9e-10  dt=4.11e-12 max_dmdt=13.1\n",
      "#617  t=4.94e-10 dt=4.11e-12 max_dmdt=13.1\n",
      "#618  t=4.98e-10 dt=4.11e-12 max_dmdt=13.2\n",
      "#619  t=5.02e-10 dt=4.11e-12 max_dmdt=13.1\n",
      "#620  t=5.06e-10 dt=4.11e-12 max_dmdt=13\n",
      "#621  t=5.11e-10 dt=4.11e-12 max_dmdt=12.8\n",
      "#622  t=5.15e-10 dt=4.11e-12 max_dmdt=12.5\n",
      "#623  t=5.19e-10 dt=4.11e-12 max_dmdt=12.2\n",
      "#624  t=5.23e-10 dt=4.11e-12 max_dmdt=11.8\n",
      "#625  t=5.27e-10 dt=4.11e-12 max_dmdt=11.4\n",
      "#626  t=5.31e-10 dt=4.11e-12 max_dmdt=11\n",
      "#627  t=5.35e-10 dt=4.11e-12 max_dmdt=10.6\n",
      "#628  t=5.39e-10 dt=4.11e-12 max_dmdt=10.1\n",
      "#629  t=5.44e-10 dt=4.11e-12 max_dmdt=9.56\n",
      "#630  t=5.48e-10 dt=4.11e-12 max_dmdt=9.03\n",
      "#631  t=5.52e-10 dt=4.11e-12 max_dmdt=8.5\n",
      "#632  t=5.56e-10 dt=4.11e-12 max_dmdt=7.96\n",
      "#633  t=5.6e-10  dt=4.11e-12 max_dmdt=7.43\n",
      "#634  t=5.64e-10 dt=4.11e-12 max_dmdt=6.91\n",
      "#635  t=5.68e-10 dt=4.11e-12 max_dmdt=6.42\n",
      "#636  t=5.72e-10 dt=4.11e-12 max_dmdt=5.95\n",
      "#637  t=5.76e-10 dt=4.11e-12 max_dmdt=5.52\n",
      "#638  t=5.81e-10 dt=4.11e-12 max_dmdt=5.5\n",
      "#639  t=5.85e-10 dt=4.11e-12 max_dmdt=5.56\n",
      "#640  t=5.89e-10 dt=4.11e-12 max_dmdt=5.59\n",
      "#641  t=5.93e-10 dt=4.11e-12 max_dmdt=5.59\n",
      "#642  t=5.97e-10 dt=4.11e-12 max_dmdt=5.55\n",
      "#643  t=6.01e-10 dt=4.11e-12 max_dmdt=5.49\n",
      "#644  t=6.05e-10 dt=4.11e-12 max_dmdt=5.41\n",
      "#645  t=6.09e-10 dt=4.11e-12 max_dmdt=5.3\n",
      "#646  t=6.13e-10 dt=4.11e-12 max_dmdt=5.17\n",
      "#647  t=6.18e-10 dt=4.11e-12 max_dmdt=5.02\n",
      "#648  t=6.22e-10 dt=4.11e-12 max_dmdt=4.85\n",
      "#649  t=6.26e-10 dt=4.11e-12 max_dmdt=4.67\n",
      "#650  t=6.3e-10  dt=4.11e-12 max_dmdt=4.48\n",
      "#651  t=6.34e-10 dt=4.11e-12 max_dmdt=4.28\n",
      "#652  t=6.38e-10 dt=4.11e-12 max_dmdt=4.07\n",
      "#653  t=6.42e-10 dt=4.11e-12 max_dmdt=3.85\n",
      "#654  t=6.46e-10 dt=4.11e-12 max_dmdt=3.63\n",
      "#655  t=6.5e-10  dt=4.11e-12 max_dmdt=3.4\n",
      "#656  t=6.55e-10 dt=4.11e-12 max_dmdt=3.18\n",
      "#657  t=6.59e-10 dt=4.11e-12 max_dmdt=2.97\n",
      "#658  t=6.63e-10 dt=4.11e-12 max_dmdt=2.76\n",
      "#659  t=6.67e-10 dt=4.11e-12 max_dmdt=2.56\n",
      "#660  t=6.71e-10 dt=4.11e-12 max_dmdt=2.38\n",
      "#661  t=6.75e-10 dt=4.11e-12 max_dmdt=2.26\n",
      "#662  t=6.79e-10 dt=4.11e-12 max_dmdt=2.31\n",
      "#663  t=6.83e-10 dt=4.11e-12 max_dmdt=2.33\n",
      "#664  t=6.89e-10 dt=6.18e-12 max_dmdt=2.34\n",
      "#665  t=6.96e-10 dt=6.18e-12 max_dmdt=2.33\n",
      "#666  t=7.02e-10 dt=6.18e-12 max_dmdt=2.28\n",
      "#667  t=7.08e-10 dt=6.18e-12 max_dmdt=2.22\n",
      "#668  t=7.14e-10 dt=6.18e-12 max_dmdt=2.13\n",
      "#669  t=7.2e-10  dt=6.18e-12 max_dmdt=2.02\n",
      "#670  t=7.27e-10 dt=6.18e-12 max_dmdt=1.91\n",
      "#671  t=7.33e-10 dt=6.18e-12 max_dmdt=1.78\n",
      "#672  t=7.39e-10 dt=6.18e-12 max_dmdt=1.64\n",
      "#673  t=7.45e-10 dt=6.18e-12 max_dmdt=1.5\n",
      "#674  t=7.51e-10 dt=6.18e-12 max_dmdt=1.35\n",
      "#675  t=7.57e-10 dt=6.18e-12 max_dmdt=1.22\n",
      "#676  t=7.64e-10 dt=6.18e-12 max_dmdt=1.09\n",
      "#677  t=7.7e-10  dt=6.18e-12 max_dmdt=0.979\n",
      "#678  t=7.76e-10 dt=6.18e-12 max_dmdt=1\n",
      "#679  t=7.82e-10 dt=6.18e-12 max_dmdt=1.01\n",
      "#680  t=7.88e-10 dt=6.18e-12 max_dmdt=1.01\n",
      "#681  t=7.95e-10 dt=6.18e-12 max_dmdt=0.994\n",
      "#682  t=8.01e-10 dt=6.18e-12 max_dmdt=0.968\n",
      "#683  t=8.07e-10 dt=6.18e-12 max_dmdt=0.933\n",
      "#684  t=8.13e-10 dt=6.18e-12 max_dmdt=0.891\n",
      "#685  t=8.19e-10 dt=6.18e-12 max_dmdt=0.843\n",
      "#686  t=8.25e-10 dt=6.18e-12 max_dmdt=0.789\n",
      "#687  t=8.32e-10 dt=6.18e-12 max_dmdt=0.732\n",
      "#688  t=8.38e-10 dt=6.18e-12 max_dmdt=0.672\n",
      "#689  t=8.44e-10 dt=6.18e-12 max_dmdt=0.612\n",
      "#690  t=8.5e-10  dt=6.18e-12 max_dmdt=0.552\n",
      "#691  t=8.56e-10 dt=6.18e-12 max_dmdt=0.496\n",
      "#692  t=8.63e-10 dt=6.18e-12 max_dmdt=0.448\n",
      "#693  t=8.69e-10 dt=6.18e-12 max_dmdt=0.42\n",
      "#694  t=8.75e-10 dt=6.18e-12 max_dmdt=0.429\n",
      "#695  t=8.81e-10 dt=6.18e-12 max_dmdt=0.432\n",
      "#696  t=8.87e-10 dt=6.18e-12 max_dmdt=0.43\n",
      "#697  t=8.93e-10 dt=6.18e-12 max_dmdt=0.422\n",
      "#698  t=9e-10    dt=6.18e-12 max_dmdt=0.41\n",
      "#699  t=9.06e-10 dt=6.18e-12 max_dmdt=0.394\n",
      "#700  t=9.12e-10 dt=6.18e-12 max_dmdt=0.375\n",
      "#701  t=9.18e-10 dt=6.18e-12 max_dmdt=0.353\n",
      "#702  t=9.28e-10 dt=9.45e-12 max_dmdt=0.323\n",
      "#703  t=9.37e-10 dt=9.45e-12 max_dmdt=0.283\n",
      "#704  t=9.47e-10 dt=9.45e-12 max_dmdt=0.243\n",
      "#705  t=9.56e-10 dt=9.45e-12 max_dmdt=0.21\n",
      "#706  t=9.65e-10 dt=9.45e-12 max_dmdt=0.184\n",
      "#707  t=9.75e-10 dt=9.45e-12 max_dmdt=0.189\n",
      "#708  t=9.84e-10 dt=9.45e-12 max_dmdt=0.188\n",
      "#709  t=9.94e-10 dt=9.45e-12 max_dmdt=0.181\n",
      "#710  t=1e-09    dt=9.45e-12 max_dmdt=0.171\n",
      "#711  t=1.01e-09 dt=9.45e-12 max_dmdt=0.158\n",
      "#712  t=1.02e-09 dt=9.45e-12 max_dmdt=0.142\n",
      "#713  t=1.03e-09 dt=9.45e-12 max_dmdt=0.125\n",
      "#714  t=1.04e-09 dt=9.45e-12 max_dmdt=0.109\n",
      "#715  t=1.05e-09 dt=9.45e-12 max_dmdt=0.0946\n",
      "#716  t=1.06e-09 dt=9.45e-12 max_dmdt=0.0805\n",
      "#717  t=1.07e-09 dt=9.45e-12 max_dmdt=0.0812\n",
      "#718  t=1.08e-09 dt=9.45e-12 max_dmdt=0.0814\n",
      "#719  t=1.09e-09 dt=9.45e-12 max_dmdt=0.0792\n",
      "#720  t=1.1e-09  dt=9.45e-12 max_dmdt=0.0752\n",
      "#721  t=1.11e-09 dt=9.45e-12 max_dmdt=0.0697\n",
      "#722  t=1.12e-09 dt=9.45e-12 max_dmdt=0.0633\n",
      "#723  t=1.13e-09 dt=9.45e-12 max_dmdt=0.0564\n",
      "#724  t=1.14e-09 dt=9.45e-12 max_dmdt=0.0494\n",
      "#725  t=1.14e-09 dt=9.45e-12 max_dmdt=0.0425\n",
      "#726  t=1.15e-09 dt=9.45e-12 max_dmdt=0.0357\n",
      "#727  t=1.16e-09 dt=9.45e-12 max_dmdt=0.0356\n",
      "#728  t=1.18e-09 dt=1.42e-11 max_dmdt=0.0355\n",
      "#729  t=1.19e-09 dt=1.42e-11 max_dmdt=0.0338\n",
      "#730  t=1.21e-09 dt=1.42e-11 max_dmdt=0.0305\n",
      "#731  t=1.22e-09 dt=1.42e-11 max_dmdt=0.0262\n",
      "#732  t=1.23e-09 dt=1.42e-11 max_dmdt=0.0214\n",
      "#733  t=1.25e-09 dt=1.42e-11 max_dmdt=0.0166\n",
      "#734  t=1.26e-09 dt=1.42e-11 max_dmdt=0.0155\n",
      "#735  t=1.28e-09 dt=1.42e-11 max_dmdt=0.0156\n",
      "#736  t=1.29e-09 dt=1.42e-11 max_dmdt=0.0147\n",
      "#737  t=1.31e-09 dt=1.42e-11 max_dmdt=0.0131\n",
      "#738  t=1.32e-09 dt=1.42e-11 max_dmdt=0.0111\n",
      "#739  t=1.33e-09 dt=1.42e-11 max_dmdt=0.00887\n"
     ]
    }
   ],
   "source": [
    "#NBVAL_IGNORE_OUTPUT\n",
    "relax_system(mesh)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Plot the magnetisation components of the relaxed state."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgMAAAF3CAYAAAA4gEgdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3XmcXHWd7//Xp7Ze0ulOZ2uyIAkkyiKQQBsdUafZFJk7wJ2rM/hQBxGG3/x+w6ij40988Bu9wx0V7ugwekcdEVBwUBTG0YxGGbZ2ZUsECYtIIAl0EpKQrbvTa1V9fn/U6U6ll3QntZxOf99PHm2qzjl16lNfCvPu7/ec79fcHREREQlXIu4CREREJF4KAyIiIoFTGBAREQmcwoCIiEjgFAZEREQCpzAgIiISOIUBERGRwCkMiIiIBE5hQEREJHAKAyIiIoFLxV1ANc2dO9eXLFlStvPt37+fGTNmlO18IVIblofasXRqw9KpDUtX7jZct27dq+4+b6LjggoDS5YsYe3atWU7X3t7O21tbWU7X4jUhuWhdiyd2rB0asPSlbsNzWzzZI7TMIGIiEjgYg0DZnarme0ws6fG2W9m9iUz22BmT5rZGUX7LjOz56Ofy6pXtYiIyPQSd8/AN4ELDrH/ncDy6Ocq4KsAZjYb+DTwRmAV8Gkza65opSIiItNUrNcMuPvPzWzJIQ65GLjd3R142MxmmdkCoA241913A5jZvRRCxXcqW7GIiEwng4ODdHR00NfXF3cpADQ1NfHss88e9utqa2tZvHgx6XT6iN53ql9AuAh4ueh5R7RtvO1Vs27zHh7emqXzt1sBMMCM6LEVPT6wnVHbjaFdZge/nhHHTfQeQ8enEglaGms4pqmWmlSyfB9YRGQa6ujoYObMmSxZsmT4/2vj1NXVxcyZMw/rNe7Orl276OjoYOnSpUf0vlM9DJTMzK6iMMRAS0sL7e3tZTnv157s46GtOXjy8bKcr9wSBq9tTnDea9Kc2ZKcEl/ysXR3d5ft30nI1I6lUxuW7mhsw6amJubMmUN3d3fcpQCQy+Xo6uo67NdlMhn27t17xO0/1cPAFuDYoueLo21bKAwVFG9vH+sE7n4TcBNAa2url+uWjZPP6OO+n/+aVatWAY579H5Q9LhouxeeU7x/rOMKNY8617jv4X7Q88Fcnu2dfWzY2c09T73CvzzRw/vfdBzXXXzKlAwEuhWpPNSOpVMblu5obMNnn32WxsbGuMsYdiQ9A0Nqa2tZuXLlEb12qoeB1cDVZnYnhYsF97n7NjO7B/hs0UWDbwc+Wc3C5jfWsrAhwbL5DdV828Py8be/jht++ju+/ouNnLywkfesek3cJYmIyBQUaxgws+9Q+A1/rpl1ULhDIA3g7v8KrAEuBDYAPcDl0b7dZva/gMeiU103dDGhHJBKJvjkO0/i6a2dfG7Ns1y8YiH1mame/0REpNpivbXQ3d/j7gvcPe3ui939Fnf/1ygI4AV/5e4nuPup7r626LW3uvuy6Ocb8X2KqS2RMD56/mvp7Muy+omtcZcjIiIjbNq0iRNPPJEPfOADrFy5kve+973cd999nHXWWSxfvpxHH32UD3/4w1x33XUA3HPPPbztbW8jn8+XrQb9mhiAM49r5nUtM7l7XQeXaqhARGRMf/+fT/PM1s6ynvPkhY18+o9PmfC4DRs2cNddd/HFL36Rc845h29/+9v88pe/ZPXq1Xz2s5/l29/+Nm94wxt461vfyoc+9CHWrFlDIlG+3+fjnnRIqsDMOOek+Tzx8l7292fjLkdEREZYunQpp556KolEglNOOYVzzz0XM+PUU09l06ZN1NfX8/Wvf53zzz+fq6++mhNOOKGs76+egUC8+YQ5fLX9BdZu3sMfvnbCBaxERIIzmd/gK6Wmpmb4cSKRGH6eSCTIZgu/xK1fv545c+awdWv5h3zVMxCI1uNmk04aD7+4K+5SRETkMG3evJkvfOELPP744/zkJz/hkUceKev5FQYCUZdJsnTuDH7/yuFPZiEiIvFxd6644go+//nPs3DhQm655RauvPLKsk6hrGGCgCyb31D2i2NERKQ0S5Ys4amnDize+81vfvOgfU8//fRBx5955pmsX7++rDWoZyAgy+Y18NLuHvoGc3GXIiIiU4jCQEBOmN9A3mHTrv1xlyIiIlOIwkBAhqZOfmGHwoCIiBygMBCQxbPqAdi2rzfmSkREZCpRGAhIY12K2nSC7Z3luwJVRESOfgoDATEzWhpr2d7ZH3cpIiIyhSgMBKZlZq16BkRE5CAKA4GZ31jDji71DIiIyAEKA4EpDBP04e5xlyIiIkxuCePly5ezc+dOAPL5PMuWLRt+Xg6agTAwxzTW0jOQo7s/y8zadNzliIhMHT+5Bl4p78x+HHMqvPP6CQ+baAnj973vfdxxxx185CMf4b777uP0009n3rzyLTqnnoHAzG8srISliwhFRKaOiZYw/uAHP8jtt98OwK233srll19e1vdXz0BgmuszAOztGYi5EhGRKWYSv8FXykRLGB977LG0tLTwwAMP8Oijj3LHHXeU9f3VMxCYprrC0MC+3sGYKxERkcNx5ZVX8r73vY93v/vdJJPJsp5bYSAwCgMiIkeniy66iO7u7rIPEYCGCYKjMCAiMrVMtITx0L7f/va3nH766Zx44ollr0FhIDCNCgMiIked66+/nq9+9atlv1ZgiIYJApNMGDNrUgoDIiJHkWuuuYbNmzfzlre8pSLnVxgIUGNdWmFARESGKQwEqKkuzb4ehQERESlQGAhQk3oGRESkiMJAgBQGRESkmMJAgBQGRESkmMJAgJrqFQZEROQAhYEANdSk6M/mGczl4y5FRCR4k1nC+MILL2TFihWsWLGCpqYmbrvttrLWEOukQ2Z2AfBFIAnc7O7Xj9h/I3B29LQemO/us6J9OWBorcmX3P2i6lR99KvPFOa07hnI0VSnPCgiAnDDozfwu92/K+s5T5x9Ip9Y9YkJj5toCeM1a9YAsG7dOi6//HIuueSSstYZWxgwsyTwZeB8oAN4zMxWu/szQ8e4+98UHf/XwMqiU/S6+4pq1Tud1GcK/9p7B3LD0xOLiEh8hpYw7urqGnMJY4BXX32V97///Xzve9+jqamprO8fZ8/AKmCDu78IYGZ3AhcDz4xz/HuAT1eptmltRk2hZ2D/QDbmSkREpo7J/AZfKRMtYZzL5bj00kv51Kc+xetf//qyv3+cfcSLgJeLnndE20Yxs+OApcADRZtrzWytmT1sZuXtL5nm6tKFMNA7kIu5EhERmYxrrrmG0047jUsvvbQi5z9aFiq6FLjb3Yv/9jrO3beY2fHAA2a23t1fGPlCM7sKuAqgpaWF9vb2shXV3d1d1vNVy4ZdhWb81SNreXV2edfEPlxHaxtONWrH0qkNS3c0tmFTUxNdXV2x1tDd3U0+n6erq4tcLsfg4CC9vb10dXUN7/v85z/PSSedxH333QfAtddey4UXXnjQefr6+o68/d09lh/gD4B7ip5/EvjkOMc+Drz5EOf6JvCuid7zzDPP9HJ68MEHy3q+alm3ebcf94kf+QPPbo+7lKO2DacatWPp1IalOxrb8Jlnnom7hIN0dnYe8WvH+izAWp/E38lxDhM8Biw3s6VmlqHw2//qkQeZ2YlAM/BQ0bZmM6uJHs8FzmL8aw1khBnRBYQ9GiYQERFiHCZw96yZXQ3cQ+HWwlvd/Wkzu45CkhkKBpcCd0YJZ8hJwNfMLE/huofrveguBDm0oVsLdQGhiIhAzNcMuPsaYM2IbZ8a8fx/jvG6XwOnVrS4aWwoDOgCQhERAc1AGKQZNYUMqJ4BERE4uOP56FTqZ1AYCFBNKoFZoWdg376X+Ps738Gv134l7rJERKqutraWXbt2HdWBwN3ZtWsXtbW1R3yOo+XWQikjM2NGJsX+/hwPPPZ/uLt/K48/+a/8oPX/ibs0EZGqWrx4MR0dHezcuTPuUoDC7YFH8pd6bW0tixcvPuL3VRgIVF0mSc9Alt90PQ7A5kSegf4uMjUzY65MRKR60uk0S5cujbuMYe3t7axcuXLiA8tMwwSBmpFJ0jOQY9PAbgCyZmx75YmYqxIRkTgoDASqPpOiZyDL1vwAx0Y3FWzZuf7QLxIRkWlJYSBQ9Zkkff2d7Egaq+oWALBl9/MxVyUiInFQGAhUfU0K698EwClzTgFg5/5XYqxIRETiojAQqNpUgmRuKwAtTa+hMe/sGdgXc1UiIhIHhYFA1aaTJH07AHMaX8NsN/YMdsdclYiIxEFhIFA1qQRQuJNgzqzjmZVIszfXG29RIiISC4WBQNWkE7gVhgVmz15Gc6KO3fmBmKsSEZE4KAwEqjaVJG891OWdTM1MmtP17HUtXCQiEiKFgUDVpBPkrJ+Z0XTczZkm9iTA8/l4CxMRkapTGAhUTSpJNjnAzOgr0Fw7m6wZ3d3bYq5MRESqTWEgULXpBIOJLDOtsDxFc91cAPbs3RhnWSIiEgOFgUDVpJL0J3M0JDIANDccA8Cezo44yxIRkRgoDASqJpWgL5GnIVEDwKwZhSmJ93ZvjbMsERGJgcJAoGrTSXoSUJ+sA6Ax6hno7N0VZ1kiIhIDhYFApRN5uhNGfWIGAI0zFwLQ2bc7zrJERCQGCgOBSvpe8mbUpWYCMHPmIgA6+7U+gYhIaBQGAmWDOwGoSRbCQCpdy4y80znQFWdZIiISA4WBQOWjMJBONQ1va3TozO6PqyQREYmJwkCgsoOFawNSyVnD2xotSWe2L66SREQkJgoDgRoYLNw1kEjOHt7WaGk68/1xlSQiIjFRGAhU32DhQkEvDgPJGjrzg3GVJCIiMVEYCFRvFAbyyfnD2xpTdXSihYpEREKjMBCo3mzhroHBRFHPQKqBTourIhERiYvCQKC6B7uoz+fpz6eGtzXWzKQvYQz06/ZCEZGQKAwEqjvXy4w89GcPDAs0Zgp3FnRqsSIRkaDEGgbM7AIze87MNpjZNWPs/4CZ7TSzJ6KfK4v2XWZmz0c/l1W38qNfV7aPulyC/mxueFtjXWHIoFOLFYmIBCU18SGVYWZJ4MvA+UAH8JiZrXb3Z0Yc+l13v3rEa2cDnwZaAQfWRa/dU4XSp4WufD91+QT9g0U9A3VzAejsfiWuskREJAZx9gysAja4+4vuPgDcCVw8yde+A7jX3XdHAeBe4IIK1TktdfkgGU/TN1jUMzCjcGdBZ8+OuMoSEZEYxBkGFgEvFz3viLaN9D/M7Ekzu9vMjj3M18o4Oj1HJp9mIFfUMzCjsIzxvh4tYywiEpLYhgkm6T+B77h7v5n9X8BtwDmHcwIzuwq4CqClpYX29vayFdfd3V3W81VTlzmJXIaOra8Mf4a+3sK1Ai+8/FzVPtfR3IZTidqxdGrD0qkNSxdXG8YZBrYAxxY9XxxtG+buxb+i3gz876LXto14bftYb+LuNwE3AbS2tnpbW9tYhx2R9vZ2ynm+asllB+jenCCdqKe5eS5tba0ADA72wLc/R01jpmqf62htw6lG7Vg6tWHp1Iali6sN4xwmeAxYbmZLzSwDXAqsLj7AzBYUPb0IeDZ6fA/wdjNrNrNm4O3RNpmE7v2FCwSTzGCg6NbCdLqe+rzTOah5BkREQhJbz4C7Z83sagp/iSeBW939aTO7Dljr7quBD5nZRUAW2A18IHrtbjP7XxQCBcB17r676h/iKNXVVRgOSNDAYM4P2tfo0DmoZYxFREIS6zUD7r4GWDNi26eKHn8S+OQ4r70VuLWiBU5TQ7cOJhKNB/UMwNAyxr1xlCUiIjHRDIQB6hq6dTAx66C7CQCatIyxiEhwFAYC1BXdOphMzh7dM5CsodOzcZQlIiIxURgIUFdfdHlFai6DI3oGGlN1dHpujFeJiMh0pTAQoM7+vQBYev4YYUDLGIuIhEZhIECdA52YO4n0nNHDBFrGWEQkOAoDAeoc6GKmQyaTYWDkrYVaxlhEJDgKAwHqzPbQ6EYmmWAge/D1AVrGWEQkPAoDAerM9tFoSTKpxOhJh7SMsYhIcBQGAtSZH6AxkSGdtFHzDGgZYxGR8CgMBKjTszQma8gkk+TyTi5/oHdAyxiLiIRHYSBAnZanMVVPOlW4h7D49sLGxkWFY/q01IOISCgUBgLj+TydBo3pBjLJwr/+4qGCmTMXAtDZvy+W+kREpPpiXahIqq+3bzdZMxozjaRTURjQMsYiIkFTGAjM0PwBjTWzyEU9A6NmIdQyxiIiQdEwQWA6u7cBhfkE0snRPQOgZYxFREKjMBCYzu7tADTWzyOTGqdnQMsYi4gERWEgMEPzBzTOmD/cM9CvZYxFRIKmMBCYzt7C/AGNMxZQM9wzMGIWwqSWMRYRCYnCQGA6+/cAhfkExr1mID2DLi1jLCISDIWBwHT2F5YvntmwcPxrBjKN9CaMwX7dUSAiEgKFgcB0DnTR4JBIpkgnC7/+j+oZqGkCYF/Xy1WvT0REqk9hIDBDyxcDwz0DoxYrqo2WMe7SMsYiIiFQGAhMZ66XRksCHJiOeGTPQH20jPF+LWMsIhIChYHA7Mv1MyuRARj/moH6aBnj/TurW5yIiMRCYSAwuz1Lc6oeYPy7CRoKyxh39r5a3eJERCQWCgOB2WtOc3omUBQGRvYMNEQrF2oZYxGRICgMBGSgv4vuhNFcMwsoHiYYMelQ02JAyxiLiIRCYSAge/duBqC5bg5w4ALCkdcMpNP11OWdfQOd1S1QRERioTAQkD37NgHQHF0gmIrmGRgccc0AaBljEZGQKAwEZE/XFgCao2sCUokoDORGh4EmLWMsIhKMWMOAmV1gZs+Z2QYzu2aM/R81s2fM7Ekzu9/MjivalzOzJ6Kf1dWt/Oi0p3sbAM0zFwFgZmSSCQbzPupYLWMsIhKOVFxvbGZJ4MvA+UAH8JiZrXb3Z4oOexxodfceM/u/gf8N/Fm0r9fdV1S16KPcnt7CvAHNzUuHt6WTNvYwQbKGl7LdVatNRETiE2fPwCpgg7u/6O4DwJ3AxcUHuPuD7t4TPX0YWFzlGqeVPb27MXeaGl8zvC2VTIw5TKBljEVEwhFnGFgEFK+E0xFtG88VwE+Kntea2Voze9jMLqlEgdPNrv49NDskU5nhbelkgoHc6GGC5sxM9iTA86ODgoiITC+xDRMcDjN7H9AK/GHR5uPcfYuZHQ88YGbr3f2FMV57FXAVQEtLC+3t7WWrq7u7u6znq7Qt3a8yGzuo5nx2gI4tW2lv33XwwT3GoBn/de9d1NS0VKymo60Npyq1Y+nUhqVTG5YurjaMMwxsAY4ter442nYQMzsPuBb4Q3cfvqLN3bdEf75oZu3ASmBUGHD3m4CbAFpbW72tra1sH6C9vZ1ynq/SvrYxR0ui/qCaZz72IHPmzaKtbeVBx/b+rB02Pc8Jy5pZdkIblXK0teFUpXYsndqwdGrD0sXVhnEOEzwGLDezpWaWAS4FDrorwMxWAl8DLnL3HUXbm82sJno8FzgLKL7wUMaw0weZG01FPCSdTIyagRBgblPhuoKde1+sSm0iIhKf2HoG3D1rZlcD9wBJ4FZ3f9rMrgPWuvtq4B+BBuAuMwN4yd0vAk4CvmZmeQqB5voRdyHICPlcll0JmBdNRTwklbBRaxMAzJt1AgCvdr48ap+IiEwvsV4z4O5rgDUjtn2q6PF547zu18Cpla1uetmz50WyZsyNZh8ckkmNfTfBvLmvA2BnNDeBiIhMX5qBMBCv7v49APOi2QeHpJMJsmMME8xoOIb6vLMzmptARESmL4WBQOzcuxGAeU1LDtqeTo49TAAwz42dA1q5UERkulMYCMTOzpcAmNt8wkHb0+NMOgQwL1HDjqwWKxIRme4UBgKxtasDc+eY+acdtP1QYWBBupFtWp9ARGTaUxgIxNbeHczPQ7pmxkHbC2sTjL5mAGBB/Tx2JGBwsGfM/SIiMj0oDARi28A+FiZqRm1PJxMMjjPl8MKGxeTN2LHjqUqXJyIiMVIYCMTWfD8LRkw4BBSWMB5vmCC6vmDrToUBEZHpTGEgALnsANsTzsLauaP2pZOJcYcJFs49BYBtuzdUtD4REYmXwkAAdu58mqwZC2aOXgE6lbTxewaOWQHAls7NFa1PRETipTAQgJdf+Q0Ai2e/dtS+whLGY4eBmtom5uecjh7NQigiMp0pDARg4871ACxd+MZR+zKpsWcgHLIkUcum/j0Vq01EROKnMBCATfs2Upt3WlpOH7UvfYhhAoAltXPZyCA+zh0HIiJy9FMYCMCm3u0cR5JEcvS6VOlkgmzeyefH7h1Y0ngcXQljty4iFBGZthQGArApu58l6aYx96WTha/AeHMNLInuKNi05aHKFCciIrFTGJjment2syXhHD9j0Zj700kDYHCc6waOX/QmAF7Y/pvKFCgiIrFTGJjmfr/xXvJmvK5lxZj7h3oGsuNcN7BwQSsz887vdj9XsRpFRCReE4YBM1tnZn9lZs3VKEjK63cdvwLgpCXnjrl/KAyMd3uhJRKcZLU827u9MgWKiEjsJtMz8GfAQuAxM7vTzN5hZlbhuqRMnt39O5ryzoJjzhhz/0TDBAAnzVjM8wySHeytSI0iIhKvCcOAu29w92uB1wLfBm4FNpvZ35vZ7EoXKKV5um87JybqsMTY/6qHLyDMjn/r4EnzTqM/Ybyw8b6K1CgiIvGa1DUDZnYa8AXgH4F/B94NdAIPVK40KVXnvpd5znKc0bhs3GOGrxk4xDwCK197EQBrX/hJeQsUEZEpYfSN5yOY2TpgL3ALcI2790e7HjGzsypZnJTm8Wfvws1oHed6ASi6ZmCcxYoAFi5sZWEO1r36JO8te5UiIhK3CcMA8G53f3GsHe7+J2WuR8pobccvSLtz2on/Y9xjMqmhawYOPcNga+08ftG/g3wuO+bkRSIicvSazDDBbjP7kJn9k5l9aein4pVJyX7W+QJnWB21dePfCJKKriWYKAy8eeFb2JMw1j97d1lrFBGR+E0mDKwBlgDrgXVFPzKFvbjxATYmnXNbRi9OVGyiWwuHvHXlX5By5/7nvle2GkVEZGqYTH9vrbt/tOKVSFnd89ubATh7xZWHPG5omOBQKxcCNDYdyyqr57/2Pc9HNFQgIjKtTKZn4Ftm9hdmtsDMZg/9VLwyOWK57ADf3/0kb6aOY44Ze+bBIcO3Fk7QMwBwyZIL2ZKEh3/ztbLUKSIiU8NkwsAAhVsKH+LAEMHaShYlpbnvoRt4JWm864RLJjz2cMLAuW/6GM155/Znbi+5RhERmTomEwY+Bixz9yXuvjT6Ob7ShcmRGRzs5SvP38UJOeOcN/3thMcPzUA4MMEwAUCmZiaXzX0Dv6KHx5/8Vsm1iojI1DCZMLAB6Kl0IVIet/3kL3kx6fz1ie8lmcpMePxkZiAs9p5z/5H5Oee6tf/IQH9XSbWKiMjUMJkwsB94wsy+plsLp7Zfr/0K/7J7HecnZnHumz8xqddMZgbCYvX1c/n0KX/BhqTz/93138jnskdcr4iITA2TCQM/AD4D/Joy31poZheY2XNmtsHMrhljf42ZfTfa/4iZLSna98lo+3Nm9o5y1HO08nyeH7X/HVc/9RWO9yTXXTL52/8O3Fo48TDBkLe98cN8ZHYrP8nt5hPfOYee7h2HXbOIiEwdE94f5u63VeKNzSwJfBk4H+igsCriand/puiwK4A97r7MzC4FbgD+zMxOBi4FTqGwouJ9ZvZad89VotapKpcdYN3627n1ya/zK3o4gxq+9N//nYaZCyZ9jsxhDhMM+eAf3QI/voIv7XqMx793DlcuPJt3vuljNM1acljnERGR+MV5s/gqYMPQVMdmdidwMVAcBi4G/mf0+G7gX6Llky8G7ozWSdhoZhui8z1UpdoZHOwhm+2mt2c3Th7yeZw87k4+n8Vx8MJz9zyO49ExDG1zj15zYFvec8P7iM7X07ubrp6ddPXuYs/+7Wzu3MTG/dt4PNfJvoTRlHc+Pv/NvOft/0w6XX9YnyOVnNx0xCNZIsEVf/wNVj75Lb6w7p/4zPZ2bvjBg5zsaU6qX8BrZi5m3oyFzG08lvraWdSkG6ipaaAm00gqXUeCBJZI0N+/g859L2EksEQSM8NIQCJR2GaG2aTW0wpaLruf/r59cZdxVFMblk5tWLpcLp6l4uMMA4uAl4uedwAjp8sbPsbds2a2D5gTbX94xGsXVa7U0T71vT/iR9lX4a5qvmtB2p1j8wnaao/hrEVv4exVHz7klMOHPNdh3Fo4ljNOez//9vr38uzzq7nnqW/x267N/GfPS/T0vQw7J3mSHxzRW8tI3427gGlAbVg6tWFJzhqs4VzeWfX3nfbTyJnZVcBVAC0tLbS3t5flvK9Lvp7GvhdJpQpNaNE/mGHR86H/ZXjfGNvt4OeJaH/hPIV/Mok6MqmZ1KRmUZueTX398SQS6UIhDg8/8tsj/hyFHgh4/oWNtCe2HPF5YBYr5/w1K+cUrmEYGHyV3r4O9g9sZzDXS877Gcz3M5gfIO+FnhPHyWazJFPJQk9KoaDhfYWPN/lrGUKWy+ZIppJxl3FUUxuWTm1YusbMvLL9PXU4DjsMmNlngX3Aze6+q4T33gIcW/R8cbRtrGM6zCwFNAG7JvlaANz9JuAmgNbWVm9rayuh5GJttLe3U77zxSd93xoWHfsa2tpOrPp7T5c2jJvasXRqw9KpDUsXVxseyWDso0AWuLHE934MWG5mS80sQ+GCwNUjjlkNXBY9fhfwgBd+lV0NXBrdbbAUWB7VJUcgnUwc8TCBiIgc/Q67Z8Ddf2BmGXcfKOWNo2sArgbuAZLAre7+tJldB6x199XALRTWRtgA7KYQGIiO+x6Fiw2zwF+FdidBOaUSxuBh3FooIiLTy4RhwMzagQ+4+6bo+Srg68Dppb65u6+hsERy8bZPFT3uA949zms/Q2H+AylRJpWYcAljERGZvibTM/A54KfRrIOLgHcCl1e0KqmqdDJBVmFARCRYk5l06B4z+0vgXuBVYKW7v1LxyqRqCtcMaJhARCRUE15AaGZ/B/wf4G0UJgBqN7M/qnBdUkXppGmYQEQkYJMZJpgDrHL3XuAhM/spcDPw44pWJlWTTib1YZyQAAAVUElEQVQOezpiERGZPiYzTPCREc83U1hPQKYJ3VooIhK2ydxN0ApcCxwXHW+Au/tpFa5NqiSdNLJ5XTMgIhKqyQwT3AF8HFgP6NfHaSidTDCgYQIRkWBNJgzsjCYAkmkqnUzQM5CNuwwREYnJZMLAp83sZuB+oH9oo7t/v2JVSVWlk5qBUEQkZJMJA5cDJwJpDgwTOKAwME3oAkIRkbBNJgy8wd1fV/FKJDbplMKAiEjIJrNq4a/N7OSKVyKxyWgGQhGRoE2mZ+BNwBNmtpHCNQO6tXCaKaxaqJ4BEZFQTSYMXFDxKiRWGiYQEQnbZGYg3FyNQiQ+GiYQEQnbZK4ZkGmucGuhegZEREKlMCC6tVBEJHAKA0IqGiZw11CBiEiIFAaETNIAdN2AiEigFAaEdLLwNcjmNVQgIhIihQEZDgODWfUMiIiESGFASKcKX4MBXUQoIhIkhQEhnRi6ZkBhQEQkRAoDcmCYQGFARCRICgMyPEyguwlERMKkMCBFtxaqZ0BEJEQKA6JhAhGRwCkMCCmFARGRoCkMCOlomGBA8wyIiARJYUDIaAZCEZGgxRIGzGy2md1rZs9HfzaPccwKM3vIzJ42syfN7M+K9n3TzDaa2RPRz4rqfoLpRdcMiIiELa6egWuA+919OXB/9HykHuDP3f0U4ALgn81sVtH+j7v7iujnicqXPH2lNEwgIhK0uMLAxcBt0ePbgEtGHuDuv3f356PHW4EdwLyqVRiQjHoGRESCFlcYaHH3bdHjV4CWQx1sZquADPBC0ebPRMMHN5pZTYXqDIKGCUREwmbulekaNrP7gGPG2HUtcJu7zyo6do+7j7puINq3AGgHLnP3h4u2vUIhINwEvODu143z+quAqwBaWlrOvPPOO4/4M43U3d1NQ0ND2c4Xl129eT72s14++PoMb1ucrup7T5c2jJvasXRqw9KpDUtX7jY8++yz17l760THpcr2jiO4+3nj7TOz7Wa2wN23RX+x7xjnuEbgx8C1Q0EgOvdQr0K/mX0D+NtD1HEThcBAa2urt7W1HfZnGU97ezvlPF9cdnT1wc/u5/hlr6XtTcdV9b2nSxvGTe1YOrVh6dSGpYurDeMaJlgNXBY9vgz44cgDzCwD/Adwu7vfPWLfguhPo3C9wVMVrXaaSyc0TCAiErK4wsD1wPlm9jxwXvQcM2s1s5ujY/4UeBvwgTFuIbzDzNYD64G5wD9Ut/zp5cBCRQoDIiIhqtgwwaG4+y7g3DG2rwWujB7/G/Bv47z+nIoWGJj08EJFurVQRCREmoFQNEwgIhI4hQEhkTBSCVMYEBEJlMKAAIVZCDVMICISJoUBAQoTDw1k1TMgIhIihQEBClMSa9VCEZEwKQwIUOgZGNRCRSIiQVIYEADSKV1AKCISKoUBAQq3Fw4oDIiIBElhQIBomEBhQEQkSAoDAhSGCbK6tVBEJEgKAwJEtxaqZ0BEJEgKAwJomEBEJGQKAwIUFivSDIQiImFSGBBAPQMiIiFTGBBgKAyoZ0BEJEQKAwIUpiNWz4CISJgUBgQYumZAYUBEJEQKAwJAKplgUKsWiogESWFAgKF5BnTNgIhIiBQGBIBM0rSEsYhIoBQGBBhawlhhQEQkRAoDAkTXDGiYQEQkSAoDAhSGCQZyedwVCEREQqMwIEBhmAAgl1cYEBEJjcKAAJBOFb4KGioQEQmPwoAAB3oGtIyxiEh4FAYEKMxACGgWQhGRACkMCHCgZ0BhQEQkPAoDAhwIA1ldMyAiEpxYwoCZzTaze83s+ejP5nGOy5nZE9HP6qLtS83sETPbYGbfNbNM9aqfnoaGCXTNgIhIeOLqGbgGuN/dlwP3R8/H0uvuK6Kfi4q23wDc6O7LgD3AFZUtd/rLaJhARCRYcYWBi4Hbose3AZdM9oVmZsA5wN1H8noZW2ooDGQ1TCAiEpq4wkCLu2+LHr8CtIxzXK2ZrTWzh81s6C/8OcBed89GzzuARRWsNQgaJhARCVeqUic2s/uAY8bYdW3xE3d3Mxvv19Hj3H2LmR0PPGBm64F9h1nHVcBVAC0tLbS3tx/Oyw+pu7u7rOeL0zO7cgCsXfcbujYmq/a+06kN46R2LJ3asHRqw9LF1YYVCwPuft54+8xsu5ktcPdtZrYA2DHOObZEf75oZu3ASuDfgVlmlop6BxYDWw5Rx03ATQCtra3e1tZ2hJ9otPb2dsp5vjjN2LQbHnuIk089jbcun1e1951ObRgntWPp1IalUxuWLq42jGuYYDVwWfT4MuCHIw8ws2Yzq4kezwXOAp7xwko6DwLvOtTr5fDURNMRD2gZYxGR4MQVBq4Hzjez54HzoueYWauZ3RwdcxKw1sx+S+Ev/+vd/Zlo3yeAj5rZBgrXENxS1eqnoZpUYWigb1BhQEQkNBUbJjgUd98FnDvG9rXAldHjXwOnjvP6F4FVlawxNEM9A/3ZXMyViIhItWkGQgGgNl3oGejXMIGISHAUBgQ40DPQN6ieARGR0CgMCAA16aFhAvUMiIiERmFAgAMXEPbrAkIRkeAoDAgAyYSRThp9uoBQRCQ4CgMyrCaVVM+AiEiAFAZkWG06oVsLRUQCpDAgw2pSSU06JCISIIUBGVaTUs+AiEiIFAZkWE06qVsLRUQCpDAgw2pSCU06JCISIIUBGVYYJlDPgIhIaBQGZFithglERIKkMCDDalIJ+jVMICISHIUBGaYLCEVEwqQwIMNq1TMgIhIkhQEZVpNO0KeeARGR4CgMyLDC2gTqGRARCY3CgAwrrE2gngERkdAoDMiwmlSSbN7J5hQIRERCojAgw2pSha+DegdERMKiMCDDhsKApiQWEQmLwoAMq69JAdAzoDAgIhIShQEZVp9JAtCrngERkaAoDMiwGZlCz8D+/mzMlYiISDUpDMiwuqhnQMMEIiJhURiQYUM9AwoDIiJhURiQYfU1Qz0DGiYQEQmJwoAMq9cwgYhIkBQGZFi9LiAUEQlSLGHAzGab2b1m9nz0Z/MYx5xtZk8U/fSZ2SXRvm+a2caifSuq/ymmn+FbC9UzICISlLh6Bq4B7nf35cD90fODuPuD7r7C3VcA5wA9wH8VHfLxof3u/kRVqp7m0skEmWSC/QoDIiJBiSsMXAzcFj2+DbhkguPfBfzE3XsqWpVQl0nqAkIRkcDEFQZa3H1b9PgVoGWC4y8FvjNi22fM7Ekzu9HMaspeYaBmZJK6gFBEJDDm7pU5sdl9wDFj7LoWuM3dZxUdu8fdR103EO1bADwJLHT3waJtrwAZ4CbgBXe/bpzXXwVcBdDS0nLmnXfeeeQfaoTu7m4aGhrKdr6p4JO/6GFRQ4KrV9ZW5f2mYxvGQe1YOrVh6dSGpSt3G5599tnr3L11ouNSZXvHEdz9vPH2mdl2M1vg7tuiv9h3HOJUfwr8x1AQiM491KvQb2bfAP72EHXcRCEw0Nra6m1tbYfxKQ6tvb2dcp5vKpj31C9pmJGhrW1VVd5vOrZhHNSOpVMblk5tWLq42jCuYYLVwGXR48uAHx7i2PcwYoggChCYmVG43uCpCtQYpPpMkp5+DROIiIQkrjBwPXC+mT0PnBc9x8xazezmoYPMbAlwLPCzEa+/w8zWA+uBucA/VKHmINRnUuzXBYQiIkGp2DDBobj7LuDcMbavBa4ser4JWDTGcedUsr6Q1WeSmmdARCQwmoFQDjIjk6JbMxCKiARFYUAO0liXorNvcOIDRURk2lAYkIM01aXpG8zTn9VQgYhIKBQG5CBNdWkA9vWqd0BEJBQKA3KQxigMdCoMiIgEQ2FADqKeARGR8CgMyEEUBkREwqMwIAdRGBARCY/CgBxkOAz0KAyIiIRCYUAO0jjcM6CJh0REQqEwIAdJJxPMyCQ1TCAiEhCFARmlqS6tMCAiEhCFARmlsS7Nvt6BuMsQEZEqURiQUebNrGFnt8KAiEgoFAZklPkza9nR2Rd3GSIiUiUKAzLKMU017OjqJ5/3uEsREZEqUBiQUVoaa8nlnV37NVQgIhIChQEZZf7MWgC2a6hARCQICgMySktjDQA7uhQGRERCoDAgo7Q0DvUM9MdciYiIVIPCgIwyb2YNZrBtn3oGRERCoDAgo6STCRbNqmPjq/vjLkVERKpAYUDGtGx+Axt2dMddhoiIVIHCgIxp2bwGXtzZTU5zDYiITHsKAzKmZfMb6M/m2bKnN+5SRESkwhQGZEzLWxoA+P32rpgrERGRSlMYkDGdvKCJdNJ4bPPuuEsREZEKUxiQMdVlkqw8tpmHXtgVdykiIlJhCgMyrj84YQ5PbdnH3h6tUSAiMp0pDMi4zj+5hbzDDx7fEncpIiJSQbGEATN7t5k9bWZ5M2s9xHEXmNlzZrbBzK4p2r7UzB6Jtn/XzDLVqTwsr1/UxOnHzuL2hzYzmMvHXY6IiFRIXD0DTwF/Avx8vAPMLAl8GXgncDLwHjM7Odp9A3Cjuy8D9gBXVLbccF199jJefHU/N977+7hLERGRCoklDLj7s+7+3ASHrQI2uPuL7j4A3AlcbGYGnAPcHR13G3BJ5aoN2/knt/CnrYv5SvsLfPR7T7C+Yx95TUQkIjKtpOIu4BAWAS8XPe8A3gjMAfa6e7Zo+6Iq1xaUz/3JacxtqOHmX2zk+7/ZQiaZYMGsWmZkUtSmE9SkkpgVjjUDw4YfF7NoQ/Hm3bv7+MaLj1bhU0xvasfSqQ1LpzYs3ez8AG1t1X/fioUBM7sPOGaMXde6+w8r9b5j1HEVcBVAS0sL7e3tZTt3d3d3Wc83la2qhZP+sJYnd2bZ0u3s6u2nP9tP34DTHV1O4A4j+ww82jBeX0Iul2P/dt2+WCq1Y+nUhqVTG5Yu3ZCL5e+VioUBdz+vxFNsAY4ter442rYLmGVmqah3YGj7eHXcBNwE0Nra6m1ljFzt7e2U83xHgz8u8/lCbMNKUDuWTm1YOrVh6eJqw6l8a+FjwPLozoEMcCmw2t0deBB4V3TcZUDVehpERESmm7huLfzvZtYB/AHwYzO7J9q+0MzWAES/9V8N3AM8C3zP3Z+OTvEJ4KNmtoHCNQS3VPsziIiITBexXEDo7v8B/McY27cCFxY9XwOsGeO4FyncbSAiIiIlmsrDBCIiIlIFCgMiIiKBUxgQEREJnLmHM5ucme0ENpfxlHOBV8t4vhCpDctD7Vg6tWHp1IalK3cbHufu8yY6KKgwUG5mttbdx11oSSamNiwPtWPp1IalUxuWLq421DCBiIhI4BQGREREAqcwUJqb4i5gGlAblofasXRqw9KpDUsXSxvqmgEREZHAqWdAREQkcAoDR8jMLjCz58xsg5ldE3c9U5WZHWtmD5rZM2b2tJl9ONo+28zuNbPnoz+bo+1mZl+K2vVJMzsj3k8wdZhZ0sweN7MfRc+XmtkjUVt9N1rQCzOriZ5viPYvibPuqcLMZpnZ3Wb2OzN71sz+QN/Dw2NmfxP9d/yUmX3HzGr1PZyYmd1qZjvM7KmibYf93TOzy6Ljnzezy8pZo8LAETCzJPBl4J3AycB7zOzkeKuasrLAx9z9ZOBNwF9FbXUNcL+7Lwfuj55DoU2XRz9XAV+tfslT1ocpLNo15AbgRndfBuwBroi2XwHsibbfGB0n8EXgp+5+InA6hbbU93CSzGwR8CGg1d1fDyQprCar7+HEvglcMGLbYX33zGw28GngjRTW5vn0UIAoB4WBI7MK2ODuL7r7AHAncHHMNU1J7r7N3X8TPe6i8H/Aiyi0123RYbcBl0SPLwZu94KHgVlmtqDKZU85ZrYY+CPg5ui5AecAd0eHjGzDoba9Gzg3Oj5YZtYEvI1ohVN3H3D3veh7eLhSQJ2ZpYB6YBv6Hk7I3X8O7B6x+XC/e+8A7nX33e6+B7iX0QHjiCkMHJlFwMtFzzuibXIIUTfhSuARoMXdt0W7XgFaosdq27H9M/D/Avno+Rxgb7TUNxzcTsNtGO3fFx0fsqXATuAb0VDLzWY2A30PJ83dtwCfB16iEAL2AevQ9/BIHe53r6LfSYUBqQozawD+HfiIu3cW7/PCLS26rWUcZvbfgB3uvi7uWo5iKeAM4KvuvhLYz4FuWUDfw4lEXdIXUwhWC4EZlPE305BNhe+ewsCR2QIcW/R8cbRNxmBmaQpB4A53/360eftQt2v0545ou9p2tLOAi8xsE4UhqXMojH/Pirpr4eB2Gm7DaH8TsKuaBU9BHUCHuz8SPb+bQjjQ93DyzgM2uvtOdx8Evk/hu6nv4ZE53O9eRb+TCgNH5jFgeXQVbYbCRTSrY65pSorGCG8BnnX3fyratRoYuhr2MuCHRdv/PLqi9k3AvqKutCC5+yfdfbG7L6HwXXvA3d8LPAi8KzpsZBsOte27ouOD/o3X3V8BXjaz10WbzgWeQd/Dw/ES8CYzq4/+ux5qQ30Pj8zhfvfuAd5uZs1RL83bo23l4e76OYIf4ELg98ALwLVx1zNVf4C3UOj+ehJ4Ivq5kMLY4f3A88B9wOzoeKNwp8YLwHoKVy7H/jmmyg/QBvwoenw88CiwAbgLqIm210bPN0T7j4+77qnwA6wA1kbfxR8AzfoeHnYb/j3wO+Ap4FtAjb6Hk2q371C4zmKQQi/VFUfy3QM+GLXnBuDyctaoGQhFREQCp2ECERGRwCkMiIiIBE5hQEREJHAKAyIiIoFTGBAREQmcwoCIlI2Z1ZnZz6LFvEo91zwz+2k56hKRQ1MYEJFy+iDwfXfPlXoid98JbDOzs0ovS0QORWFARCZkZm+I1lavNbMZ0Zr2rx/j0PcSzaRmZm1m1m5md5vZ78zsjqFV68xsk5l9zsyeMLO1ZnaGmd1jZi+Y2V8Wne8H0TlFpIJSEx8iIqFz98fMbDXwD0Ad8G/u/lTxMdHU3Me7+6aizSuBU4CtwK8ozGX/y2jfS+6+wsxupLDe+1kUZq17CvjX6Ji10XuKSAUpDIjIZF1HYV2OPuBDY+yfC+wdse1Rd+8AMLMngCUcCAND63msBxrcvQvoMrN+M5vl7nspLN6ysKyfQkRG0TCBiEzWHKABmEnhN/iResfY3l/0OMfBv4AM7cuPOC5fdFxtdF4RqSCFARGZrK8BfwfcAdwwcqe77wGSZjZWUDhSr6UwbCAiFaQwICITMrM/Bwbd/dvA9cAbzOycMQ79LworVZbL2cCPy3g+ERmDVi0UkbIxszOAv3H395fpfD8HLo56HUSkQtQzICJl4+6/AR4s16RDwD8pCIhUnnoGREREAqeeARERkcApDIiIiAROYUBERCRwCgMiIiKBUxgQEREJnMKAiIhI4P5/cs+K4sdNYr4AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 576x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#NBVAL_IGNORE_OUTPUT\n",
    "plot_magnetisation([['m0.npy', 'mx'], ['m0.npy', 'my'], ['m0.npy', 'mz']])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The DW is at the maximum value of $|m_z|$ or $|m_y|$. Consequently, the domain wall position is:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Maximum |m_z| at x = 124\n"
     ]
    }
   ],
   "source": [
    "m0_z = np.load('m0.npy').reshape(-1, 3)[:, 2]\n",
    "x = np.arange(len(m0_z))\n",
    "index_max = np.argmax(np.abs(m0_z))\n",
    "\n",
    "print('Maximum |m_z| at x = %s' % x[index_max])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Using the obtained domain wall equilibrium state, we now simulate its motion in presence of a spin-polarised current."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "#NBVAL_IGNORE_OUTPUT\n",
    "excite_system(mesh, 1.5e-9, 151)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We plot once again to compare the initial state with the ones after a SP current was applied."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0, 200)"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgkAAAF3CAYAAAAvotjUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3XmcXFWd///Xp6r3fUl3unpLd0hnJZCQBRDBsIrMDOCOjgqKMqPy05n56ogP5guKo1/9jvOd7+DX0UFkREVRURE0gCxpBwQSEgKEhCUhe9KdTtJJ73ud3x9VHTpJJ2mSe7vqVr2fj0dbVbduffqcdJu8Ofecc805h4iIiMiRQolugIiIiCQnhQQREREZl0KCiIiIjEshQURERMalkCAiIiLjUkgQERGRcSkkiIiIyLgUEkRERGRcCgkiIiIyLoUEERERGVdGohswmUpKStyMGTMS3Qzf9PT0kJ+fn+hm+CKV+wbqX9Cpf8GVyn0DWLNmzT7nXMXJfj6tQsLUqVNZvXp1opvhm+bmZpYtW5boZvgilfsG6l/QqX/Blcp9AzCzbafyeV1uEBERkXElNCSY2V1m1mZmLx/jfTOz281sk5m9ZGZnjXnvWjPbGP+6dvJaLSIikh4SPZLwI+Dy47z/LqAp/nUD8D0AMysDbgXOBpYCt5pZqa8tFRERSTMJnZPgnPtvM2s4zilXAT92zjngWTMrMbMIsAx41DnXDmBmjxILGz/3t8UiIpJKzIwtW7bQ39+f6KackpycHGpra8nMzPS0brJPXKwBdox5vTN+7FjHj6t7YD/f/c0Xj/GujfnfNx/doWfxY8YRZ4x9dvSrsccNw8Z9++iDNs7xN7/1+N9j7959bOpccei7jXJjWjD+x8f7PuN9Dzt02A47asds0zH/nOz4f352RL3WlhZ29T13wu/xZvnx/kyP/Oke1SgAQhhhC5NpGWRY+NBXeOwjsechCx+zTUc2wWG4UAYulEU0lIELZeJCWbhQBptb95P98hu4cCbOYu9hx/p9OV7P39oJdoITQgYZ4RBZ4RAZYSMzHCIz/pgRtvjx+LFQiFDohC0SSSr5+fkUFhbS0NAw7t8bQeCcY//+/ezcuZPGxkZPayd7SDhlZnYDsUsV5DTk8P2uhxPcIh9lAgcS3QifhIG2RDfiaBnOURCNUhCNUhgd+zxKSTRK5fAIU0dGqBoepmp4hMqRkXGv8S0FePXwY90uh25y6XJ5dI157HT57HGltFIWe3RlbHOV9JEzCT0+vpww5GYYuRmjj0ZOBmTbML/b9EdKso2SHKM025iaHyI7HMy/lI/U3d1Nc3Nzopvhm1TuX2FhIVlZWXR3dye6KackKyuLgwcPev5zSvaQsAuoG/O6Nn5sF7FLDmOPN49XwDl3B3AHQMNpde57Z37ryDPAuSM/8+Z7Y585xjk+5rPxzx1ebfQsd+S3AcCOONuN1jyyTUd8j/G+/+bNm5neOP3wNnFkm9wRbT2y70cfG31jvL/O3Thtio7X0XHbdMR5o20a5+Pbdmyjvq7usEqOo//8Rt9x0aOOHPG9j92+KFGGo1GG3TDDbuTQ4wgj8dfxY9ERBt0wvdF+ekf66Y320x3tpy3aR+9IPx0jPQy6ocO+S45lUZ9VwbTMKUzLKKcpo5ymjFL2bN9GfU0Eiw5iI0NYdJDwUC/hwS7KhzqpHOwmPNRFeHAvGYObyOzff1QP+gvq6SuZSW/pTHpL59BZuYShvMpx/nziPR7vj+4IUecYGokyNOIYHokyOBJleMQxHI0yGD82+v7AcJTu/mG6B4boHhimqz/2daB/iLaOHrpaD/+zMIO60jxmTi1gRmUhp9cUsaShjKlFiQ87b1WqL6NL5f6tXbuWoqKiRDfDEzk5OSxcuNDTmskeEh4AbjSze4lNUuxwzrWY2SPAN8ZMVrwM+PKJimVn5PH2BVf419oEa+5rZtk5yxLdDF80Nzez7B3LEt2Mt8Q5R+dgJ3t699DW28bu7t1sOriJjQc2surgRh7tfhGA3IxcGosb+cCMc7mw/kLKcspOXHx4ELpbobMFOnfB/k3ktG0gp+0VSl9eAW4kdl55E8x8J8z+S6g7G0KJmavc3NzM295+AXu7B9jT2U/LwX42tXXzelsXm/Z086fX9zI0EkstjVPyuXh2JZfOncqShjJdwhBJoISGBDP7ObERgSlmtpPYioVMAOfc94HlwBXAJqAX+Hj8vXYz+xowepH6ttFJjCLJwswozi6mOLuYmaUzD3vPOcfevr280PYCz7U+x6NvPMpXnvkKX3v2a1xYdyEfmPUBzomcc+xrpBlZUFIf+zrS8ADseRm2/hm2/AlW/ic88/+gZBosuhYWfgwKTnoDtpOWlRGipiSXmpJcOKLZQyNRNuzu5Lmt7Ty5cR8/fmYbdz61hbqyXD6wqI4Pn11PeUH2pLdZJN0lenXDh07wvgM+e4z37gLu8qNdIn4zMyrzKrms4TIua7iMt/W9jciZEZZvXs5vN/2Wx7Y/xtzyuXx+4ec5t/rctzahKiMbahbFvs77HPR3wusPw/M/hsdvgz/9Cyz9FJz3d5Bf7l8n34LMcIgz60o4s66ET54/ne6BYR7bsIdfrt7Bvz76Ov/R/AYfPXcaN1wwnSkKC5Jitm7dyuWXX84555zD008/zZIlS/j4xz/OrbfeSltbG/fccw/33HMP5eXl3HLLLTzyyCN8/etfp7m5mZDPo4PJfrlBJC2YGbPLZjO7bDafXfhZlm9ezvdf/D5/89jfcF71edx67q1ECiInVzynCM74QOxr72vw39+Gp78Da+6Gy74GZ33sOKtTEqMgO4OrF9Zw9cIaNrV18x8rNnHnk5v5xXM7+PK7ZvOBxXW6DCGe++qD69mwu9PTmnOri7j1r+ad8LxNmzbxq1/9irvuuoslS5bws5/9jKeeeooHHniAb3zjG/zsZz9jyZIlnH/++Xzuc59j+fLlvgcESPxmSiJyhOxwNu9uejcPvvtB/nHJP/J82/O8+4F389uNvz314hWz4L0/gM88A1Xz4cHPwU+uhq7WU6/tkxmVBfyfDy7gj39/AbOrCrnpN+v48J3PsrdrINFNE/FMY2Mj8+fPJxQKMW/ePC6++GLMjPnz57N161by8vL4wQ9+wKWXXsqNN97IaaedNint0kiCSJLKCmfx0bkf5cK6C7n16Vu55elbeLX9Vb645ItkhE7x/7qVc+DaB+H5H8EjN8MdF8KHfgbV3s6M9tKMykLuveEcfvHcDr7y4Hr+6jtP8b2PnMXCem22Kt6YyH/x+yU7+83LaKFQ6NDrUCjE8PAwAOvWraO8vJzdu3dPWrs0kiCS5GoLa7nj0jv42NyP8bNXf8aNT9xI33DfqRcOhWDxJ+ATj0AoDHe9C15L7n1EzIxrltbz60+/jcwM44N3PMuK15JwAw0Rj23bto1//dd/Ze3atTz00EOsXLlyUr6vQoJIAIRDYb645Ivccu4tPL3raf6++e8ZGhk68QcnInIGfGoFVM6GX34U3njCm7o+mlddzAOffTtNlQX8zU/W0KygICnMOcf111/Pt7/9baqrq/nhD3/IJz/5yUnZSlohQSRA3j/z/dx67q38edef+dKTX2I4OuxN4YIK+MhvYMpM+PmHYdvT3tT1UWl+Fvd88myaKgu44SdrWLVFq6AlmBoaGnj55TdvhvyjH/2I973vfYfeW79+PY899hhXXnklAIsWLWLdunXk5Pi/8ZhCgkjAvHfme/nC4i/w6LZH+c7a73hXOK8MPno/FNfCLz4KnZN33fNkleTFgkJtSS6fuWcNrR3BvkmPSLJRSBAJoGvnXcv7Zr6Pu16+iyd3Puld4YIKuOZnMNQH930CvLqk4aOSvCz+86OL6Bsc4dP3rGFgeCTRTRJJGQoJIgH1pSVfYmbpTG5+6mb29OzxrnDFTLjydtj+TGzzpQBomlrIt99/Jmu3H+RfHn4t0c0RSRkKCSIBlZORw7ff8W36R/r56jNfHXNjMg/Mfx8s+nhs06Vda7yr66N3zY/w12fXc9eft7BuZ0eimyOSEhQSRAKssbiRGxfcyJO7nqR5R7O3xS+9DQqmwu//AaLBGML/x8tnU16QzU2/eYnhkeiJPyAix6WQIBJwH5rzIWaUzOBbz32L/mEPJ+7lFME7vw4tL8Ca//Kuro+KczP5yl/NY/3uTn709NZEN0ck8BQSRAIuM5TJl5d+mV3du7jrZY/veXb6e6HxAnjsNugNxhLDK+ZXsWxWBf/++EY6+pJ/4qVIMlNIEEkBSyNLuWzaZfxo/Y842H/Qu8Jm8K7/DQOd8Ox/eFfXR2bGF985i67+YX741JZEN0ck0BQSRFLE3575t/QN93HPq/d4W7hyDsy9Elb+J/R5GEB8NK+6mMvnVfFfT22ho1ejCZLctm7dyuzZs7nuuuuYOXMmf/3Xf81jjz3GeeedR1NTE6tWreLzn/88t90WW230yCOPcMEFFxCN+j/vRjd4EkkRTaVNXFR3Efe8cg/XzbuO/Mx874qf/wXY8DtYdQe84x+9q+ujz1/SxMPrW7nzqc38j8tmJbo5EgQP3QSt67ytWTUf3vXNE56mW0WLiO8+dcan6Brs4hev/cLbwpEzYOa7YpccBrq8re2TOZEi/mJ+hP/681a6+jWaIMlNt4oWEd+dPuV0zo2cy93r7+Yjcz5CVjjLu+IXfBHuvAjW/hTO+bR3dX10wwXT+cO6Fu5fu4uPntuQ6OZIspvAf/H7RbeKFpFJcd2862jvb+eJ7R7fzbF2EVSfBWvuBi83bvLRGbXFnF5TxD0rt3u72ZTIJNOtokXEE+dUn0MkP8JvNv7G++KLroW9r8DO57yv7QMz46/PnsarrV2s2XYg0c0ROSm6VbSIeCZkIa6ecTXPtjzL7m6PhyVPfy9k5sdGEwLiyjOrKcjO4J6V2xPdFJFx6VbRIjKprp5xNQD3b7rf28LZhTD/vbD+N9Df6W1tn+RnZ/Ces2r4w7oW2nsGE90ckUBRSBBJQdUF1ZwTOYf7N93PiNf3XTjrOhjqhZfv87auj65ZUs/gcJTl61oS3RSRQFFIEElR72l6Dy09LaxqXeVt4ZqzoGIOrAtOSJgTKWT6lHyFBJG3SCFBJEUtq1tGbkYuj29/3NvCZjD3Ktj2NHS3eVvbJ2bGFfMjPLt5P/u6BxLdHJHAUEgQSVE5GTmcV30eT2x/gqjzePvWuVcCDl5b7m1dH10xP0LUwSPrWxPdFJHAUEgQSWEX1V/E3r69rNvn8VazlXOhbDpseMDbuj6aEymkUZccRN4ShQSRFHZB7QVkWIb3GyuZwZy/gi1/CsxNn2KXHKp45o397NclB5EJUUgQSWHF2cUsrlrME9uf8H7HwTlXQnQYXn/E27o+Gr3k8McNexLdFJFAUEgQSXEX11/M1s6tbO7Y7G3h6rOgsBpeCc4lh7mRImpKcml+LRgTLiU96FbRx2BmlwP/DoSBO51z3zzi/X8DLoy/zAMqnXMl8fdGgNELrdudc1dOTqtFguWi+ov4+sqv88T2JzitxMM7x4VCMOtd8OK9MDIE4UzvavvEzDi/aQp/eKmF4ZEoGWH9d5K86VurvsWr7a96WnN22Wy+tPRLJzxPt4o+gpmFge8C7wLmAh8ys7ljz3HO/b1zboFzbgHwHWDsZvR9o+8pIIgcW2VeJbPLZvNsy7PeFz/tQhjqgZ2rva/tk/ObKugaGObFncGYSyHpQbeKPtpSYJNzbjOAmd0LXAVsOMb5HwJunaS2iaSUpVVLuffVe+kf7icnw8P93hveDhaCzc0w7Vzv6vrovBnlmMGTG/exaFpZopsjSWQi/8XvF90q+mg1wI4xr3fGjx3FzKYBjcDYKdo5ZrbazJ41s6v9a6ZI8J0dOZvB6CBr29Z6Wzi3FKoXxkJCQJTkZXFGbQlPbtyX6KaITFiibhWd0DkJb8E1wH3OubGb0E9zzu0ys+nAE2a2zjn3xpEfNLMbgBsAKioqaG5unpQGJ0J3d3fK9i+V+wb+928gOkCIEL9a+SsGSr1d/tcYbqRux2/582PLGcnIG/ecZPv51WcN8ofNQ/zh0RXkZ9op10u2/nktlftXVFREV1dXQtvQ3d1NNBo91I6hoSH6+vro6uqiu7ubkZERrrvuOr72ta9RWFjI7bffzic+8Qmam5sPuxNkf3+/9z8n51xCvoBzgUfGvP4y8OVjnLsWeNtxav0IeN+JvufMmTNdKluxYkWim+CbVO6bc5PTv48u/6i75sFrvC/8RrNztxY599rDxzwl2X5+Kzfvd9O+9Hv30LoWT+olW/+8lsr9e/755xPdBM9s2LDhqGPAancK/1Yn8nLDc0CTmTWaWRax0YKj1lKZ2WygFHhmzLFSM8uOP58CnMex5zKICLFLDhvaN9A56PEtnuvOhoycQF1yWFhfQn5WmCc37k10U0SSWsJCgnNuGLgReAR4Bfilc269md1mZmNXK1wD3BtPRKPmAKvN7EVgBfBN55xCgshxnF11NlEXZXWrxysRMnOg/txAhYTMcIiljWWs2tKe6KaIJLWEzklwzi0Hlh9x7JYjXn9lnM89Dcz3tXEiKebMijPJzchlZctKLqq/yNvi098Bj30FuvdCQYW3tX2yuKGMFa+9xsHeQUryshLdHJGkpJ1ERNJEZjiTBRULWL3Hhz0N6uPLH3c+531tnyyaVgrA2u3aLyHdHT5QHUx+9UEhQSSNLKhcwKaDm+gd6vW2cORMCGXAzlXe1vXRmbUlhEPGmm0HEt0USaCRkRH2798f6KDgnGP//v2HrXTwSlCWQIqIB+ZPmU/URVm/fz1LqpZ4VzgzF6rOgB3BGUnIzQozr7qI1ds0LyGd9fT00NXVxd69wZ7EmpOTQ21tred1FRJE0sj8KbGpPC/ufdHbkABQuwTW/gRGhiEcjL9aFk0r5d5VOxgaiZKp+zikJeccjY2NiW5G0tL/K0TSSElOCQ1FDby09yXvi9cthaFeaFvvfW2fLJpWSt/QCK+2JHYzHZFkpZAgkmbOqDiDl/a+5P012NrFsccdwZmXMDp5UZccRMankCCSZs6Ycgb7+/ezu8fjm8SUTIP8ykCtcIgU51JdnKPJiyLHoJAgkmbOqDgDwPtLDmaxSw4BCgkAixrKeF4hQWRcCgkiaaaptImccI4/8xJqF0P7ZugJzh0Wz6wtZndHP/u6vb3xlUgqUEgQSTMZoQzmls/1KSQsjT3u9GHDJp/Mqy4GYP1uj+9pIZICFBJE0tAZFWfwSvsrDI0MeVs4EruUQes6b+v6aG51EQDrd3ckuCUiyUchQSQNzSmbw1B0iM0dm70tnF0IZdOh9UVv6/qoODeT+rI81u/SSILIkRQSRNLQ7LLZALx24DXvi1fND9RIAsC86iKNJIiMQyFBJA3VF9WTHc7mtXafQsKBrdAfnH9051UXsXV/L539Hl9+EQk4hQSRNJQRyqCppMmnkHBm7HFPcHZenFcTm7z4iiYvihxGIUEkTc0qm8WrB171fufFqtj9IYJ0yWFefPLiywoJIodRSBBJU7PKZtEx0MGe3j3eFi6sgrwp0OrDEkufVBbmUFGYrXkJIkdQSBBJU4cmL3p9ycEskJMXT68uYoNGEkQOo5AgkqZmls4EfFzh0PYKeL0Pg4/mVRezsa2b/qGRRDdFJGkoJIikqfzMfOoK63i1/VXvi0fOhJFB2Pe697V9Mq+6iJGo4/U9um20yCiFBJE0Nqt0Fq8f8OEf8gBOXpxZVQjA63u6E9wSkeShkCCSxmaVzWJ753Z6h3q9LVw+AzJyAhUSppXlkRUOsbFNIwkioxQSRNLYrNJZOJz3owmhMEyZCXt9uJThk4xwiOkV+WzUSILIIQoJImlsRskMALZ0bPG+eMUs2BucOQkATVMLNZIgMoZCgkgaqy6oJiuU5f2NniAWEjq2w0Bw/su8qbKAHe199A4OJ7opIklBIUEkjYVDYRqKG3wKCbF9GIK0wmHm1AIANrUFJ9iI+EkhQSTNnVZ8Gm8cfMP7wlNmxR4DFBJmVMZWOGhegkiMQoJImmssaWR39276h/u9LVzWCKHMQE1ebCiPrXB4XfMSRACFBJG0N714Og7H1s6t3hYOZ8aWQu71YUdHn4yucNikkQQRQCFBJO1NL54OwOaDfsxLmBmokAAwo7JAIwkicQkNCWZ2uZm9ZmabzOymcd6/zsz2mtkL8a9PjnnvWjPbGP+6dnJbLpI6phVNI2xh3ujwYV5CxWw4sIXQyKD3tX0yc2qhVjiIxGUk6hubWRj4LnApsBN4zswecM5tOOLUXzjnbjzis2XArcBiwAFr4p89MAlNF0kpWeEs6grr/NsrwUXJ7dvtfW2fjK5weKOth/m1xQlujUhiJXIkYSmwyTm32Tk3CNwLXDXBz74TeNQ51x4PBo8Cl/vUTpGU11jc6M/lhvgKh7zeHd7X9snoCgfd6EkksSGhBhj7N8fO+LEjvdfMXjKz+8ys7i1+VkQmYHrxdLZ1bWMo6vGtnctngIXI7wlOSJhWnkc4ZGzZ15PopogkXMIuN0zQg8DPnXMDZvY3wN3ARW+lgJndANwAUFFRQXNzs+eNTBbd3d0p279U7hskvn+D3YMMR4f59eO/piqzytPaS3OqyOrcEqif35QcWPnKFpqzWyZ0fqJ/fn5L5f6lct+8kMiQsAuoG/O6Nn7sEOfc/jEv7wT+95jPLjvis83jfRPn3B3AHQCzZs1yy5YtG++0lNDc3Eyq9i+V+waJ71/Fvgp+8oefUDGzgmXTPG7H7jOJ7toQqJ/f3C2raO0cYNmy8yd0fqJ/fn5L5f6lct+8kMjLDc8BTWbWaGZZwDXAA2NPMLPImJdXAq/Enz8CXGZmpWZWClwWPyYiJ6GxuBGALZ0+TF4sP43cvhaIRr2v7ZPGKQVs3deDcy7RTRFJqISNJDjnhs3sRmL/uIeBu5xz683sNmC1c+4B4HNmdiUwDLQD18U/225mXyMWNABuc861T3onRFJEXmYeFbkVbOvc5n3xsumEo4PQ1QLFwZg61Dglj76hEfZ0DlBVnJPo5ogkTELnJDjnlgPLjzh2y5jnXwa+fIzP3gXc5WsDRdJIXWEd2zu3e1+4/LTYY/sbAQoJsWWQW/b1KCRIWtOOiyICxDZV2t7lQ0goi4eE/T5s1uSThil5AFrhIGlPIUFEAKgvqmdf3z56h3q9LVxUw0goKzaSEBDVxblkZYTYul8hQdKbQoKIAFBfWA/g/WhCKER/ThXs92GzJp+EQkZDeR6b9yokSHpTSBARIHa5AfBlXkJvXiRQIwkAjVPy2bJPd4OU9KaQICJAbOIi+DCSAPTlRqB9S6CWQTZMyWd7ey8jUS2DlPSlkCAigL/LIPtyq2FkADp3el7bL9On5DM04th1oC/RTRFJGIUEETnEr2WQfbnxfdGCtMKhPB+ALZq8KGlMIUFEDvFrGWRfbnXsSXtwJi82VsRDwl7NS5D0pZAgIof4tQxyILsMMnIDFRIqCrLJzwprrwRJawoJInKIb8sgLQRl0wN1ucHMaJiSz7Z2j/eNEAkQhQQROaS+KB4SfNmeeXrglkHWl+WxXSFB0phCgogc4ttIAsS2Z27fAtER72v7pL4sj53tfUS1DFLSlEKCiBzi690gSxsgOgSdu72v7ZO6sjwGR6Ls6epPdFNEEkIhQUQOU1dYx84uH/YzKG2IPR70IYD4pL4sdqOn7ft1yUHSk0KCiBympqCGnd1+hITYts8c2Op9bZ8cCgmalyBpSiFBRA5TU1jDnp49DI0MeVu4uC62yuFAcEYSqktyCRns0K6LkqYUEkTkMDUFNTgcLT0t3hYOZ0JRbaBGErIyQkSKc9mhkQRJUwoJInKY2oJaAP8uOQRoTgJoGaSkN4UEETlMbWEsJOzq3uV98dJpgRpJAIUESW8KCSJymIrcCjJCGf6scChpgO49MBSca/z15Xns7RqgbzA4+zuIeEUhQUQOEw6Fqc6v9mkkoSH2eNCHzZp8UluaC8COAxpNkPSjkCAiR6kpqGFXl0+XGyBQlxy0V4KkM4UEETlKbWGtvyMJAVoGqb0SJJ0pJIjIUWoKajgwcICeIY9vk5xfAZl5gRpJKMvPIj8rrJAgaUkhQUSOUlNYA+D95EUzKAnWMkgzo64sT3slSFpSSBCRo4zulaBlkDH1ZXlsU0iQNKSQICJH8TckNMTmJLjg3H65tjSPXQf6cAFqs4gXFBJE5CjF2cXkZ+b7tFfCNBjsgt5272v7pKY0l76hEQ70enw/C5Ekp5AgIkcxs9gySL8uN0CgLjnUlMT2Sth9MDibQIl4QSFBRMblW0goqY89duzwvrZPRjdU2qm7QUqaSWhIMLPLzew1M9tkZjeN8/4/mNkGM3vJzB43s2lj3hsxsxfiXw9MbstFUt9oSPD8OnxxXewxQCGhOj6SsEsjCZJmEhYSzCwMfBd4FzAX+JCZzT3itLXAYufcGcB9wP8e816fc25B/OvKSWm0SBqpLqimb7iPjoEObwvnFENWIXT4MN/BJ6V5meRmhtmlkQRJM4kcSVgKbHLObXbODQL3AleNPcE5t8I5N7ru6FmgdpLbKJK2IvkRAFp6WrwtbAYldXAwOCMJZkZNaa7mJEjaSWRIqAHG/i2xM37sWK4HHhrzOsfMVpvZs2Z2tR8NFElnoyFhd89u74sX10JHcG7yBLHJi7rcIOkmI9ENmAgz+wiwGHjHmMPTnHO7zGw68ISZrXPOvTHOZ28AbgCoqKigubl5MpqcEN3d3Snbv1TuGyRn/7pGugB48oUnCW8On1KtI/vX1BOmct9W/pxkfT4e6x1ga9vwuD+nZPz5eSmV+5fKffNCIkPCLqBuzOva+LHDmNklwM3AO5xzA6PHnXO74o+bzawZWAgcFRKcc3cAdwDMmjXLLVu2zLseJJnm5mZStX+p3DdIzv4557jtntsoiBSwbMmyU6p1VP/Cz8Puh1h27mLILjil2pNlvdtE887XOPtt55ObdXhoSsafn5dSuX+p3DcvJPJyw3NAk5k1mlkWcA1w2CoFM1sI/CdwpXOubczxUjPLjj+fApwHbJi0loukATOjKr/Kn8sNh5ZBBmfyYo1WOEj7C6j7AAAgAElEQVQaSlhIcM4NAzcCjwCvAL90zq03s9vMbHS1wr8ABcCvjljqOAdYbWYvAiuAbzrnFBJEPBbJj9Da0+p94eL4HOQALYOsKVVIkPST0DkJzrnlwPIjjt0y5vklx/jc08B8f1snItUF1TTvaPa+cJD3StAySEkj2nFRRI6pKr+K/f37GRgZOPHJb0VhFYQyArUMcmphNuGQaRmkpBWFBBE5ptFlkJ5fcgiFoag6UHMSMsIhqopydLlB0opCgogcU3VBNeDDhkoQu+QQoMsNEN8rQZcbJI0oJIjIMVXlVwHQ0u1XSAjOSALEJi9qJEHSiUKCiBxTVV4Vhvk0klALnbthZNj72j6pKcmltbOf4ZFoopsiMikUEkTkmDLDmVTkVvgTEkrqwI1Alw+1fVJTmstI1LGny+OJnCJJSiFBRI4rUhDx6XJD8PZK0DJISTcKCSJyXJH8iE+XG+K7LgZoGeTorotaBinpQiFBRI5rNCREncfX4YvjN30N0EiCtmaWdKOQICLHFSmIMBQdor2/3dvCWfmQVx6okJCbFaY8P4udutwgaUIhQUSOa3RDpd3dPtzoKYDLIKtLtAxS0odCgogc12hI8G0ZZIDmJEDskoPmJEi6UEgQkeOKFMRDgh8rHErqYyMJznlf2yc1pbFdF12A2ixyshQSROS4irKKKMgs8G8kYagH+g54X9snNSW59A2NcKB3KNFNEfGdQoKInFBVfpV/92+AQE1e1F4Jkk4UEkTkhPzbKyG+oVKA5iXUlmoZpKSPE4YEM1tjZp81s9LJaJCIJJ/qgmqftmaOb6gUoBUO2itB0slERhI+CFQDz5nZvWb2TjMzn9slIkmkKr+KjoEOeod6vS2cVw4ZuYG63FCSl0leVliXGyQtnDAkOOc2OeduBmYCPwPuAraZ2VfNrMzvBopI4lXnVwM+LIM0i11yCFBIMDOqtQxS0sSE5iSY2RnAvwL/AvwaeD/QCTzhX9NEJFmMLoP0Z0OlYO6VoMsNkg4yTnSCma0BDgI/BG5yzo3eI3WlmZ3nZ+NEJDn4uqFSSR289rL3dX1UU5rLul0diW6GiO9OGBKA9zvnNo/3hnPuPR63R0SSUEVuBWEL+7cMsmcvDPVBZq739X1QU5JLe88gvYPD5GVN5K9RkWCayG93u5l9DmgYe75z7nN+NUpEkks4FGZq3lSf90rYBVNmeF/fB2/eMrqfGZUFCW6NiH8mEhKWA88C6wCP7xUrIkERKYj4szXz6F4JHTuCExLG7JWgkCCpbCIhIcc59w++t0REklokP8KaPWu8L1wSvF0X3xxJ0ORFSW0TWd3wEzP7lJlFzKxs9Mv3lolIUonkR2jrbWM4Ouxt4cJqwAK1oVJlYTbhkCkkSMqbSEgYJLb08RlgTfxrtZ+NEpHkEymIMOJG2Nu719vCGVlQGAnUMsiMcIiphdlaBikpbyKXG/4HMMM5t8/vxohI8hq7DHJ03wTPBGxDJUAbKklamMhIwibA471YRSRoRndd3N3j04ZKAbrcALGQ0NLRn+hmiPhqIiMJPcALZrYCGN1ISUsgRdJMVX4VAK09rd4XL66FV38P0SiEgnFz2uqSXB5+uZVo1CW6KSK+mcj/G+8Hvg48zZtzEjyZ4mxml5vZa2a2ycxuGuf9bDP7Rfz9lWbWMOa9L8ePv2Zm7/SiPSJybHmZeRRnF/sUEupgZDC2qVJA1JTkMDgSZV/PwIlPFgmoE44kOOfu9uMbm1kY+C5wKbCT2F0mH3DObRhz2vXAAefcDDO7BvgW8EEzmwtcA8wjdofKx8xspnNuxI+2ikhMJD/i04ZK8b0SOndC4VTv6/ugesyGSiKpKpHjekuBTc65zc65QeBe4KojzrkKGA0p9wEXx29TfRVwr3NuwDm3hdi8iaWT1G6RtFWVX+VvSAjQvIRIsfZKkNSXyJBQA4ydzrwzfmzcc5xzw0AHUD7Bz4qIxyL5EVq7fZqTAIEKCdpQSdJByt+ZxMxuAG4AqKiooLm5ObEN8lF3d3fK9i+V+wbB6V9vRy9dQ1089MRD5IYmfjOmE/bPOc4P5bB7/TO8MTDv1Bs6CZxz5IRh1csbqaobCsTP72QF5ffzZKRy37zwlkOCmX2D2H/R3+mc238K33sXUDfmdW382Hjn7DSzDKAY2D/BzwLgnLsDuANg1qxZbtmyZafQ5OTW3NxMqvYvlfsGwelf/5Z+fvffv2PGwhk0lTZN+HMT6t/6adQVQl0A/hxG1a39E6GCAgoKugLx8ztZQfn9PBmp3DcvnMzlhlXAMPBvp/i9nwOazKzRzLKITUR84IhzHgCujT9/H/CEc87Fj18TX/3QCDTF2yUiPhpdBunbvIQAXW6A+IZKHbrcIKnrLY8kOOfuN7Os+GTDk+acGzazG4FHgDBwl3NuvZndBqx2zj0A/JDYvSM2Ae3EggTx834JbCAWWD6rlQ0i/hvdddG3vRJaX/K+ro+qS3JZv7sDyEx0U0R8ccKQYGbNwHXOua3x10uBHwBnnuo3d84tJ3Yr6rHHbhnzvB94/zE++3Vi+zeIyCSZkjuFDMvwaSShLrZPwlAfZE58vkMiVRfnsK97kMGRlJ/eJWlqIr/Z/wt42MxuJ7aC4F3Ax31tlYgkpXAozNT8qT7vlbAbyk/zvr4PRvdKaO/XrouSmiaymdIjZva3wKPAPmChc86HsUYRCYKq/Cpauv3cK2GHQoJIkjjhxEUz+5/Ad4ALgK8AzWb2Fz63S0SSVCQ/4tOchPhWJx3jLlRKSqN7Jezviya4JSL+mMjlhnJgqXOuD3jGzB4G7gT+4GvLRCQpRfIj7Ondw0h0hHAo7F3hotGQEJwVDlOLszHTSIKkrhOOJDjn/i4eEEZfb3POXepvs0QkWVXlVzHiRtjb5/HNmDKyoWBq7HJDQGRnhKkoyGa/QoKkqImsblgM3AxMi59vgHPOneFz20QkCY1dBjm6b4JnArhXQqQkl/29nYluhogvJnK54R7gi8A6QBfeRNLcaEho6WlhAQu8LV5cC3s2nPi8JFJTksPz7R2JboaILyYSEvbGNzYSEfF518U6eP2P4ByYeV/fB9XFuTzW73DOYQFps8hETSQk3GpmdwKPAwOjB51zv/GtVSKStAqyCijMKvRvGeRwH/S2Q3659/V9UF2Sy+AIHOwdojQ/K9HNEfHURELCx4HZxPYdHb3c4ACFBJE05d8yyDF7JQQoJADsOtinkCApZyIhYYlzbpbvLRGRwIjkR/zddbFjJ1R7PN/BJ6N7Jew+2MfpNcUJbo2ItyZyF8inzWyu7y0RkcCoyq+itdePkYT4HeADtMIhUpIDxEKCSKqZyEjCOcALZraF2JwELYEUSXOR/AgdAx30DvWSl5nnXeG8csjIgc7ghITy/CwyQrC7oz/RTRHx3ERCwuW+t0JEAmXsXgnTS6Z7V9gscHslmBnlOcYujSRICprIDZ62TUZDRCQ4IgVv7pXgaUiAwIUEgPJco0UhQVLQROYkiIgcZuyGSp4rCmBIyAmx+6AuN0jqUUgQkbdsSu4Uwhb2b4VDVysMD3pf2ydlOcaern6GRrQpraQWhQQRecsyQhlU5lX6uFeCg67d3tf2SVmu4Ry0avKipBiFBBE5KVX5Vf7vlRAQ5Tmxv0q1DFJSjUKCiJyUqvwqn7ZmDt5eCeU5sXs27O5QSJDUopAgIiclkh+htbeVqPP4OnxxTeyxY4e3dX1UlhsPCZq8KClGIUFETkokP8JwdJj9ffu9LZyZC3lTAjWSkB02yvKzdLlBUo5CgoicFF+XQRbXQscu7+v6qLokRyFBUo5CgoiclKr8KsDPkBCckQSA6uJc7booKUchQUROyuiui/4sg6yLzUlwzvvaPqkpzWXXgT5cgNosciIKCSJyUgozC8nPzPdvJGGwG/o7vK/tk5qSXHoGR+joG0p0U0Q8o5AgIifFzIjkR3xaBjm6V0JwVjjUlsbuhrnzgC45SOpQSBCRk+bbhkol8b0SDgYpJOQCCgmSWhQSROSkVedX+xQSpsUeD273vrZP3gwJvQluiYh3FBJE5KRVF1RzcOAgPUM93hbOK4fMvECFhOLcTPKzwlrhICklISHBzMrM7FEz2xh/LB3nnAVm9oyZrTezl8zsg2Pe+5GZbTGzF+JfCya3ByICUFMQ2x1xd7fHN2Myg5J6OLjN27o+MjNqS/N0uUFSSqJGEm4CHnfONQGPx18fqRf4mHNuHnA58H/NrGTM+190zi2If73gf5NF5EjVBdWADyEB4iEhOCMJ8OYySJFUkaiQcBVwd/z53cDVR57gnHvdObcx/nw30AZUTFoLReSERkPCrm4fdkcc3SshQGpLczUnQVJKokLCVOfc6GynVmDq8U42s6VAFvDGmMNfj1+G+Dczy/apnSJyHOU55WSHs/0bSeg7AP2d3tf2SU1JLp39w3T1a68ESQ3m1+5gZvYYUDXOWzcDdzvnSsace8A5d9S8hPh7EaAZuNY59+yYY63EgsMdwBvOuduO8fkbgBsAKioqFv3yl7886T4lu+7ubgoKChLdDF+kct8g2P37513/TCQrwvUV1x/znJPpX0XbU8zb8C88t/jf6SloOMVW+mu0f6tah/mPFwb42nm51BWmzrzwIP9+nkgq9w3gwgsvXOOcW3zSBZxzk/4FvAZE4s8jwGvHOK8IeB5433FqLQN+P5HvO3PmTJfKVqxYkegm+CaV++ZcsPv3N4/+jfvAgx847jkn1b8dq527tci5V5efXMMm0Wj/1m4/4KZ96ffu0fWtiW2Qx4L8+3kiqdw355wDVrtT+Pc6UVH3AeDa+PNrgd8deYKZZQG/BX7snLvviPci8UcjNp/hZV9bKyLHVJNf49/lBgjU5MXRvRK0DFJSRaJCwjeBS81sI3BJ/DVmttjM7oyf8wHgAuC6cZY63mNm64B1wBTgnye3+SIyyre9EvKnQEZuoEJCeX4WOZkhTV6UlJGRiG/qnNsPXDzO8dXAJ+PPfwr89Bifv8jXBorIhI3dK6GptMm7wmax7ZkDtldCTYluGS2pI3Vm1ohIQvi/V0KwlkHWaEMlSSEKCSJySnzdKyGAGyrVakMlSSEKCSJySvzfK6EdBrq8r+2T2tJc9vcM0js4nOimiJwyhQQROSVmRnVBNbt7/FzhEJxLDnWleQDsaNdoggSfQoKInLLqgmqftmYO3jLI+rJYSNjerhUOEnwKCSJyyvzfKyE4KxwUEiSVKCSIyCnzba+EgkrIzIMDW72t66OSvEwKszPYoZAgKUAhQUROWU1hbK+EnV07vS1sBqUNgQoJZkZdWZ5GEiQlKCSIyCmrK6wDfAgJELiQAFBXlquRBEkJCgkicspGQ8KOLh9WIZQ2xkKCT3es9UN9fCTBBajNIuNRSBCRU1aUVURxdrFPIaEBhnqhu8372j6pL8tjYDjK3q6BRDdF5JQoJIiIJ2oLav0JCWWNsccDW7yv7ZM6rXCQFKGQICKeqCusY2e3T3MSIFDzErQMUlKFQoKIeKKusI6W7haGox5vR1xSDxi0B2ckoaY0FzOFBAk+hQQR8URtYS3DbpiWnhZvC2dkQ1FNoEYSsjPCRIpyFBIk8BQSRMQTWgZ5uLqyPC2DlMBTSBART/i7DLIhUBMXAW2oJClBIUFEPFGZV0lmKNOfkYSyBujeA4PB+Ue3viyPPZ0D9A+NJLopIidNIUFEPBGyEDUFNT6tcBhdBrnV+9o+GV3hsPNAcIKNyJEUEkTEM3WFdf7tugiBCgmjeyVs3aeQIMGlkCAinhkNCZ5vR3xor4TgzEuYPiUfgK37Pb4zpsgkUkgQEc/UFtbSM9TDwYGD3hbOK4OswkDtlVCan0VxbiZb9ikkSHApJIiIZ3xb4WAW2565fbO3dX3WOCVfIUECTSFBRDxTX1gPwLbObd4Xn9IE+zd5X9dHjVPy2aqQIAGmkCAinqktrCVkIX9CQvkMOLgdhoNzZ8XGKfns7uinb1DLICWYFBJExDNZ4SxqCmr8Cwm4QF1yaIhPXtzWrtEECSaFBBHx1LSiaWzt3Op94fLTYo8BuuQwusJhy16FBAkmhQQR8VRDUQPbOrcRdVFvC5fPiD0GKCSMjiRs0TJICSiFBBHxVENRA33DfbT1tnlbOKcY8isDFRIKsjOoKMzWSIIEVkJCgpmVmdmjZrYx/lh6jPNGzOyF+NcDY443mtlKM9tkZr8ws6zJa72IHE9DcQOAT5ccZsC+4IQEiK9w0EiCBFSiRhJuAh53zjUBj8dfj6fPObcg/nXlmOPfAv7NOTcDOABc729zRWSiphVNA2Brx1bvi0+ZEaiRBIDGcu2VIMGVqJBwFXB3/PndwNUT/aCZGXARcN/JfF5E/DU1byq5Gbn+rXDo3Qd9B7yv7ZPGinz2dQ/S2T+U6KaIvGWJCglTnXMt8eetwNRjnJdjZqvN7FkzGw0C5cBB59xw/PVOoMbHtorIW2BmNBQ1sKXThy2UD01efMP72j5pHL2Hg0YTJIAy/CpsZo8BVeO8dfPYF845Z2bHuhvMNOfcLjObDjxhZuuAjrfYjhuAGwAqKipobm5+Kx8PlO7u7pTtXyr3DVKvf7n9ubza+eqhPnnVv7yedpYCrzz1IHuquk+5nleO1799XbFVHsv/ezXt1b79leurVPv9HCuV++YF335jnXOXHOs9M9tjZhHnXIuZRYBxp0E753bFHzebWTOwEPg1UGJmGfHRhFpg13HacQdwB8CsWbPcsmXLTrJHya+5uZlU7V8q9w1Sr3/rX1jP2hfX8rbz30ZWOMu7/g0PwurPMacykzlJ9Od1vP71D43wT08/THZFPcuWzZzchnkk1X4/x0rlvnkhUZcbHgCujT+/FvjdkSeYWamZZcefTwHOAza42D1oVwDvO97nRSRxphVNw+HY3rnd28IZWVAyLVCTF3Myw9SX5bGpLXlGPkQmKlEh4ZvApWa2Ebgk/hozW2xmd8bPmQOsNrMXiYWCbzrnNsTf+xLwD2a2idgchR9OautF5LgaixoBLYMc1VRZyMa2rkQ3Q+QtS8gFMufcfuDicY6vBj4Zf/40MP8Yn98MLPWzjSJy8g4tg/QjJFTMgq1PQnQEQmHv6/ugaWoBf3q9jaGRKJlh7WEnwaHfVhHxXEFWAZW5lWw+6MPNmCrnwHA/HNjqfW2fzJxawNCIY5s2VZKAUUgQEV/MKJ3BpoM+XBaomBN7bNtw/POSSFNlIQCv79G8BAkWhQQR8UVTSRNvHHyDkeiIt4UrZsUe2171tq6PTqsowAxe36N5CRIsCgki4osZpTMYjA6yo2uHt4WzC6CkPlAjCblZsRUOG7XCQQJGIUFEfNFU0gTgzyWHyrnQ9or3dX3UVFnARo0kSMAoJIiILxqLGzGMjQc3el+8Yjbs3xjbXCkgmqYWsmVfD0Mj0UQ3RWTCFBJExBd5mXnUFtay6YBPIwnRYWgPzj0cmiq1wkGCRyFBRHwzo8SnFQ6VoyscgnPJYeZUrXCQ4FFIEBHfNJU2sa1zG0PO49skT5kJFgpUSNAKBwkihQQR8U1TSRMjboS2oXHv4XbyMnOgbDrsDU5IyM0KU1eax0aNJEiAKCSIiG9mlMwAoGWoxfviFbMDNZIAMKuqkFdaOxPdDJEJU0gQEd9MK5pGRiiD3YO7vS9eORfaN8NQv/e1fXJ6dTFb9vXQMzCc6KaITIhCgoj4JjOcSUNRgz8jCVPngYsGalOledVFOAevtGg0QYJBIUFEfDWrbBY7B3d6XzhyZuyx5QXva/vk9JpiANbvVkiQYFBIEBFfzS2by8GRg+zr2+dt4dIGyCmB3cEJCVOLsinPz2L97o5EN0VkQhQSRMRXc8vnArBhv8eXBcxiowktL3pb10dmxtzqIl7epZEECQaFBBHx1ZzyORjmfUgAqF4Qm5MQoO2ZT68pZmNbF4PD2p5Zkp9Cgoj4Kj8zn8qMStbvX+998ciZMDIYuMmLQyNOmypJICgkiIjv6rLr/BlJiCyIPQboksPp1aOTFzUvQZKfQoKI+K4uq4623jbvJy+WTYfs4kCtcKgvy6MgO0MrHCQQFBJExHf1WfWAX5MXzwjUSEIoNDp5USMJkvwUEkTEd3VZdRjm37yE1pdhxOObSPloXnURr7R0MTyiyYuS3BQSRMR32aFsGoobfFrhsBBGBmDvq97X9snC+lL6hkZ4tVWTFyW5KSSIyKSYWz6XDft8CgkAO1d7X9sni6eVArBm24EEt0Tk+BQSRGRSzJ8yn7a+Nlp7Wr0tXDYd8itgx0pv6/qouiSXSHGOQoIkPYUEEZkUZ1WeBcCaPWu8LWwGdWfD9me8reuzs6aVKiRI0lNIEJFJMbN0JvmZ+Ty/53nvi9efCwe2QqcPd5v0yeJppew62EdLR1+imyJyTAoJIjIpwqEwCyoW8HybTyEBYMez3tf2yaL4vITntx1McEtEjk0hQUQmzVlTz2LTwU0c7Pf4H8bIGZCRC9uDExLmRIrIzQyzelt7opsickwKCSIyaUbnJaxtW+tt4XAm1C4O1LyEzHCIM+uKeV7zEiSJJSQkmFmZmT1qZhvjj6XjnHOhmb0w5qvfzK6Ov/cjM9sy5r0Fk98LEXmr5lfMJzOU6X1IgNglh9Z1MBCcvQcWTStl/e5O+gZHEt0UkXElaiThJuBx51wT8Hj89WGccyuccwuccwuAi4Be4I9jTvni6PvOueBs3C6SxrLD2Zw+5XTWtHm8wgGg/hxwUdj5nPe1fbJ4WhnDUcfz2zWaIMkpUSHhKuDu+PO7gatPcP77gIecc72+tkpEfLewciEb9m2gb9jjWf21S8BCsC04lxyWNpaRGTae3Ojxja9EPJKokDDVOTe6VqkVmHqC868Bfn7Esa+b2Utm9m9mlu15C0XEF4umLmLYDfNCm8cDgDlFUH0WbF7hbV0f5WdncFZ9KU9u3JvopoiMy5xz/hQ2ewyoGuetm4G7nXMlY8494Jw7al5C/L0I8BJQ7ZwbGnOsFcgC7gDecM7ddozP3wDcAFBRUbHol7/85cl3Ksl1d3dTUFCQ6Gb4IpX7BunVv4HoADftuIkLii7g3aXv9vT7NGz5OdO2/ZI/n3c3w5lFntY+nlP5+T34xiC/3jjE7RfmUZRtHrfMG6n8+5nKfQO48MIL1zjnFp/s5zO8bMxYzrlLjvWeme0xs4hzriX+D37bcUp9APjtaECI1x4dhRgws/8CvnCcdtxBLEgwa9Yst2zZsrfQi2Bpbm4mVfuXyn2D9OvffX+8j629W73v82n58MN7eXvVIMz3uPZxnMrPr2zGQX698c+4qTNZtqDG24Z5JJV/P1O5b15I1OWGB4Br48+vBX53nHM/xBGXGuLBAjMzYvMZXvahjSLik7fXvJ3NHZvZ3b3b28I1Z0FuKWx63Nu6PppXXUxpXib//brmJUjySVRI+CZwqZltBC6Jv8bMFpvZnaMnmVkDUAf86YjP32Nm64B1wBTgnyehzSLikfNrzwfgqV1PeVs4FIbTLoJNj0E06m1tn4RDxnkzpvDkxr34dflX5GQlJCQ45/Y75y52zjU55y5xzrXHj692zn1yzHlbnXM1zrnoEZ+/yDk33zl3unPuI8657snug4icvMaiRmoKanhy55PeF59xCfS0wZ7gDDCe3zSFtq4BXt+jv8okuWjHRRGZdGbG22vezsrWlQyODHpb/LSLYo+bHvO2ro/Ob6oAYMVrx5ueJTL5FBJEJCHeXvN2+ob7vL91dGEVVM2HjX888blJorokl/k1xTy0Ljh3sZT0oJAgIgmxtGop2eFsntj+hPfFZ/9l7GZPAbp19BXzI7y4s4Md7dozTpKHQoKIJEReZh4X1F7AH7f9keHosLfF570HcLDhfm/r+ugv5kcAeOjl4AQbSX0KCSKSMFc0XkF7fzurWlZ5W7hiZuySw8u/9rauj+rL8zi9pog/rGtNdFNEDlFIEJGEOb/2fAoyC1i+Zbn3xee9J3azpwPbvK/tkyvmR3hxx0F2HtAlB0kOCgkikjDZ4Wwurr+Yx7c/zsDIgLfFT39P7HH9b72t66NDlxw0miBJQiFBRBLqisYr6B7q9n7PhNIGqFkML9/nbV0fTSvP54zaYn79/E5trCRJQSFBRBJqaWQpZTll/H7z770vPv/90LoOWl70vrZPrllSz6utXTy//WCimyKikCAiiZURyuCqGVfRvKOZlm6PZ/af+UHIzINVP/C2ro+uWlBNQXYG96wMzlwKSV0KCSKScNfMugaH4xev/cLbwrmlcMYHYN2voLfd29o+yc/O4N0La/j9Sy0c7PV4N0qRt0ghQUQSrrqgmovqLuK+jffRP9zvbfEln4Lhflj7U2/r+ujDZ9czOBzlvjU7E90USXMKCSKSFD4858N0DHR4vxyy6nSYdh48dydER7yt7ZM5kSIWTSvlp89uY3gkGHezlNSkkCAiSWHx1MXMLJ3JTzb8hKjz+B/GpZ+Cg9vglQe9reujT53fyNb9vTzw4u5EN0XSmEKCiCQFM+O6edex6eAm/rjN45szzbkSpsyCFV8PzGjCZXOrmBMp4vbHN2o0QRJGIUFEksYVjVcwo2QG33n+OwxFh7wrHArDRf8E+16HlzyeHOmTUMj4u0ua2Lq/l/tf0GiCJIZCgogkjXAozOfP+jzbu7bz240e75Q4568gsgBW/C8Y9nh3R59cNncq86pjowlDGk2QBFBIEJGk8o7ad7CwciHff/H79A55eA8DM7j4FujYHpvEGABmxhcum8X29l5+8OTmRDdH0pBCgogkFTPjHxb9A3v79nL72tu9LX7aRTDjUnjin+HAVm9r++TC2ZW8c95U/v2xjWzb35Po5kiaUUgQkaSzoHIBH579Ye555R5Wt672rrAZ/OW/gYXgwc9DQO6P8NUrTyczHOKf7n9Z93SQSaWQICJJ6fNnfZ7aglpuefoWby87lNTBpV+FzW67sowAAA6YSURBVM3w/N3e1fVRVXEOX3znLJ7cuI+fr9qR6OZIGlFIEJGklJeZx23n3caOrh189Zmvevtf0Is+AY0XwENfgt1rvavro4+cM43zm6bwlQfW88IO3fxJJodCgogkrSVVS/jcws+xfMtyfrDOw5s0hULw3rsgvwJ+/mHo2uNdbZ+EQ8bt1yyksiibT/90Dfu6g7FCQ4JNIUFEkton53+Sv5j+F3xn7Xd4eMvD3hUuqIBrfgb9B+HeD0F/p3e1fVKan8X3P7KI9p5BPvGj5+jo83AvCZFxKCSISFIzM776tq9yVuVZ3PTkTTz4hodbK0fOgPf8AFpehJ9cDX0HvKvtk9NrivneR87ilZZOPvbDlQoK4iuFBBFJetnhbL53yfdYPHUxNz91Mz9/9efezVGY85fwgR9D6zq4+0roTP7dDS+aPZXvf2QRG1o6+fAPnmXnAQ8ndoqMoZAgIoGQl5nHdy/5Lu+ofQffWPkNbn7qZu9WPcz+C7jm57D/Dfj++fDGCm/q+ujiOVO542OL2b6/lyv/3595etO+RDdJUpBCgogERnY4m/974f/lMws+w+83/54P/v6DrGpZ5U3xpkvghhWQPwV+8m54+MtJP0/hwlmV/O7G8yjPz+Kvf7iSW3/3Ml39uvwg3lFIEJFACYfCfPrMT/ODy37AUHSI6/94PV/40xfY0rHl1ItXzIJPPQGLPw7Pfg/+32J4/scwPHjqtX0yvaKA3372PK49t4EfP7uNS/7Pn7hn5TYGhoNxt0tJbgoJIhJIZ0fO5v6r7uczCz7Dn3b8iavuv4q/W/F3rGpZRdSdws2QsvJjuzJ+6nEoroUH/j+4fQH8+XbobvOuAx4qyM7gK1fO47efOY/qklxu/u3LXPgvzXyv+Q3auvoT3TwJsISEBDN7v5mtN7OomS0+znmXm9lrZrbJzG4ac7zRzFbGj//CzLImp+UikkxyMnL49Jmf5uH3PsynzvgUz7U+x/V/vJ53/vqdfPu5b7OyZSVDIyc5/F6zCD75OHzk11DaCI/+T/jX2fDT98FzP4SDybfz4YK6En7z6bfxk+uXUleWx7cefpVz/9cTXHvXKn767DZ2H+xLdBMlYDIS9H1fBt4D/OexTjCzMPBd4FL+//buP0aOs77j+Ptzt/btne/s812Oi+3YsU0dI1MkbDWUKgTF1KJxVHBpqyotKqkApUiFlqKqpEp/RFWlNq0KUqWKNAhUtzg10MZwiNAEUK6hrUhiHIf4Ry65S9zmbCdOHEL8I3XOu9/+Mc/hvcuufb7bu7nd+7ys0cw+Mzv3/fqZ2Xn2mdkZGAUelTQQEYeAO4HPRsRuSXcBHwE+N/thm9l81Nveyyc2f4KPvu2jDD43yDdGvsGuJ3ex89BOiq1Frum5hk09m9jUmw1XL72aYqF46RVL8FPbsuHFIdh/Dxy8F4a/nc1fugpWbmbNuW4YLkH/W6GzP3tfTiRx/YY+rt/Qx8iLp/nq3lG+deA4f/y1AwBctbyda9f28JYru9jQ38mGN3Wxqrudlpb8Yrb5K5dGQkQchmxjvoh3AMMR8UxadjewQ9Jh4D3Ab6TldgJ34EaC2YLXXmhn+7rtbF+3nbNjZ3n4+MM88vwjHDp5iIGRAXYP7f7Jsr3FXlZ2rmTFkhX0L+lnedtylrUto7utm+62bpYsWkKxUKS90J6Nl6+hbduf0bLtjqzBMPJdOLoPjj3G+pdH4NkvZSsuFGHZ6uwZEctWZ3d17OiBjl5o74HiUljUkQ2LOy5MF+rfIfrmvk5u2/4WPn3jRoZPnOahp19i75GX+e+Rl9jz2NEL/2+LWlnT08GblrbRv7RI/9I2+jrb6CouoqtYYORkid7RH9NZLNBVLNCxuJVFrS0UWnSpz3FrcHn1JEzFKqCyP28U+FmgF3glIs5XlK+a49jMbJ7rWNTB1jVb2bpmKwDlKHPk1SMcPnmY0VOjHDtzjGOnjzH0oyG+d/R7vHZ+al3xba1tFFoKtKo1G1/Zzfmet9GxqEChXKK1XKKlNEbh/LO0nngKlc+jdE+H8cNp5ViMzxOoJY2z12L8IFw5zs4SS6pYz+QDtapOjr/YvBJiBZQiKAWUykE5gtJ5GD0ZPPfSGx+Q+eBAjf8QVeSTGgwXbTboEvNnyRv/jzLlKPP5IV+eV8usNRIkfQe4ssqs2yPi67P1d6vEcStwK0BfXx+Dg4Nz9afn3OnTp5s2v2bODZzfXFrCEjamfywmG5bDWIxxpnSGM+UznC2f5Vz5HK/H6xeGcjYeizFKUaJMmXKUKVPmXOEcLa0tlFvKlLgwL1suUJRQeQzFeYgSijJEGSin6RJQJiIQkZoNQRAoYkJJdvSOSTeTupwLNeONL5WGinkafxUT3zHh3fGGtVX7C/NaCMTl/RKkkfKbqVlrJETEthmu4iiwuuL1VansJNAtqZB6E8bLa8VxN3A3wMaNG+OGG26YYVjz1+DgIM2aXzPnBs6v0Tm/xtXMuQHot2fWbzOf+1geBTakXzIsBm4GBiJrPj8I/Gpa7hZgznomzMzMFoq8fgL5AUmjwM8B35R0fypfKek+gNRL8HHgfuAw8JWIOJhW8WngU5KGya5R+MJc52BmZtbs8vp1wx5gT5XyY8BNFa/vA+6rstwzZL9+MDMzs1kyn083mJmZWY7cSDAzM7Oq3EgwMzOzqhST75jRxCSdAobyjmMWXQE060Plmzk3cH6Nzvk1rmbODWBjRHRN983z+Y6Ls2EoImo+UKrRSdrbrPk1c27g/Bqd82tczZwbZPnN5P0+3WBmZmZVuZFgZmZmVS20RsLdeQcwy5o5v2bODZxfo3N+jauZc4MZ5regLlw0MzOzqVtoPQlmZmY2RQuikSDpRklDkoYl3ZZ3PDMlabWkByUdknRQ0u+l8jskHZW0Pw03XWpd85WkI5KeSHnsTWU9kr4t6ek0Xp53nNMhaWNFHe2X9KqkTzZy/Un6oqQTkg5UlFWtL2X+Lu2PP5S0Jb/IL61Gbn8j6ckU/x5J3al8raTXKurwrvwin5oa+dXcFiX9Uaq7IUm/kE/UU1cjvy9X5HZE0v5U3oj1V+t4UJ/9LyKaegBagRFgPdmT4x8HNuUd1wxzWgFsSdNdwFPAJuAO4A/yjq9OOR4BrphU9tfAbWn6NuDOvOOsQ56twPPA1Y1cf8C7gS3AgUvVF9nzWb4FCHgn8HDe8U8jt/cChTR9Z0VuayuXa4ShRn5Vt8X0OfM40AasS5+trXnncLn5TZr/t8CfNnD91Toe1GX/Wwg9Ce8AhiPimYh4HdgN7Mg5phmJiOMRsS9NnyJ7SuaqfKOaEzuAnWl6J/BLOcZSLz8PjETE/+QdyExExEPAy5OKa9XXDuCfIvN9oFvSirmJ9PJVyy0iHojsSbUA3weumvPA6qRG3dWyA9gdEeci4llgmHn+sL2L5SdJwK8B/zKnQdXRRY4Hddn/FkIjYRXwXMXrUZrogCppLbAZeDgVfTx1IX2xUbvjkwAekPQDSbemsv6IOJ6mnwf68wmtrm5m4gdUs9Qf1K6vZtsnP0z2zWzcOkmPSfoPSdfnFVQdVNsWm63urgdeiIinK8oatv4mHQ/qsv8thEZC05LUCfwb8MmIeBX4HPBm4O3AcbJutEb1rojYAmwHfkfSuytnRtZv1tA/zZG0GHg/8NVU1Ez1N0Ez1Fc1km4HzgO7UtFxYE1EbAY+BdwjaWle8c1A026Lk/w6ExvpDVt/VY4HPzGT/W8hNBKOAqsrXl+VyhqapEVkG8SuiLgXICJeiIhSRJSBzzPPuwEvJiKOpvEJYA9ZLi+Md4ul8Yn8IqyL7cC+iHgBmqv+klr11RT7pKTfAn4R+GD6ECZ1w59M0z8gO2d/TW5BTtNFtsWmqDsASQXgl4Evj5c1av1VOx5Qp/1vITQSHgU2SFqXvrndDAzkHNOMpPNoXwAOR8RnKsorzyt9ADgw+b2NQNISSV3j02QXiR0gq7db0mK3AF/PJ8K6mfAtplnqr0Kt+hoAPpSusn4n8OOKbtGGIOlG4A+B90fE2YryPkmtaXo9sAF4Jp8op+8i2+IAcLOkNknryPJ7ZK7jq5NtwJMRMTpe0Ij1V+t4QL32v7yvzJyLgexqzqfIWoW35x1PHfJ5F1nX0Q+B/Wm4Cfhn4IlUPgCsyDvWaea3nuwK6seBg+N1BvQC3wWeBr4D9OQd6wxyXAKcBJZVlDVs/ZE1do4DY2TnOD9Sq77Irqr++7Q/PgH8TN7xTyO3YbLzuuP7311p2V9J2+x+YB/wvrzjn2Z+NbdF4PZUd0PA9rzjn05+qfwfgY9NWrYR66/W8aAu+5/vuGhmZmZVLYTTDWZmZjYNbiSYmZlZVW4kmJmZWVVuJJiZmVlVbiSYmZlZVW4kmFndSGpPt7NtrcO6+iT9ez3iMrPpcSPBzOrpw8C9EVGa6Yoi4kXguKTrZh6WmU2HGwlmdkmSrk0P+ymmO2IelPTTVRb9IOnObpJukDQo6V8lPSlpV7o7HJKOSPpLSfsl7ZW0RdL9kkYkfaxifV9L6zSzHBTyDsDM5r+IeFTSAPAXQDvwpYiYcNvodNvz9RFxpKJ4M/BW4BjwX8B1wH+mef8bEW+X9Fmyu99dBxTJbgF8V1pmb/qbZpYDNxLMbKr+nOxZKP8H/G6V+VcAr0wqeyTSvfEl7QfWcqGRMP4MlSeAzog4BZySdE5Sd0S8QvZQmpV1zcLMpsynG8xsqnqBTqCL7Bv/ZK9VKT9XMV1i4heT8XnlScuVK5YrpvWaWQ7cSDCzqfoH4E+AXcCdk2dGxI+AVknVGhDTdQ2N/zRMs4blRoKZXZKkDwFjEXEP8FfAtZLeU2XRB8ieSlcvW4Fv1nF9ZnYZ/BRIM6sbSVuA34+I36zT+h4CdqReCjObY+5JMLO6iYh9wIP1upkS8Bk3EMzy454EMzMzq8o9CWZmZlaVGwlmZmZWlRsJZmZmVpUbCWZmZlaVGwlmZmZWlRsJZmZmVtX/A0hJ8kSeoteGAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 576x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#NBVAL_IGNORE_OUTPUT\n",
    "# We can plot the m_x component for a number snapshots\n",
    "# to observe the DW motion\n",
    "# We will plot the 100th and 150th files (we can also compute\n",
    "# until the system reaches ~5 ns to improve the effect)\n",
    "plot_magnetisation([['m0.npy', 'mx'],\n",
    "                    ['dyn_npys/m_100.npy', 'mx'],\n",
    "                    ['dyn_npys/m_150.npy', 'mx']])\n",
    "plt.xlim([0, 200])"
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "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"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}