Skip to content

Instantly share code, notes, and snippets.

@lstagner
Last active October 18, 2018 19:30
Show Gist options
  • Save lstagner/2715086505844d8dc10b23ccb435eb44 to your computer and use it in GitHub Desktop.
Save lstagner/2715086505844d8dc10b23ccb435eb44 to your computer and use it in GitHub Desktop.
Callback MWE
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [],
"source": [
"#using Pkg\n",
"using DifferentialEquations\n",
"using EFIT #pkg> add https://github.com/JuliaFusion/EFIT.jl\n",
"using Equilibrium #pkg> add https://github.com/JuliaFusion/Equilibrium.jl\n",
"using GuidingCenterOrbits #pkg> add https://github.com/JuliaFusion/GuidingCenterOrbits.jl#ls/diffeq\n",
"using PyPlot"
]
},
{
"cell_type": "code",
"execution_count": 28,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(AxisymmetricEquilibrium{Float64}, Limiter: npoints = 202)"
]
},
"execution_count": 28,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"g = readg(EFIT.test_gfile)\n",
"M, wall = read_geqdsk(EFIT.test_gfile)"
]
},
{
"cell_type": "code",
"execution_count": 22,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Orbit Period: 151.14338172423936 microseconds\n",
"Time to orbit apogee: 85.15709668165579 microseconds\n"
]
}
],
"source": [
"gcp = GCParticle(80.0,-0.515,1.8,-0.1)\n",
"path, stat = get_orbit(M,gcp,tmax=500, dt=1,integrator=ImplicitMidpoint())\n",
"println(\"Orbit Period: \",stat.tau_p*1e6,\" microseconds\")\n",
"println(\"Time to orbit apogee: \",stat.tm*1e6, \" microseconds\")"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Orbit Period: 150.6169249621677 microseconds\n",
"Time to orbit apogee: 73.45285561723068 microseconds\n"
]
}
],
"source": [
"gcp = GCParticle(80.0,-0.515,1.8,-0.1)\n",
"path, stat = get_orbit(M,gcp,tmax=500, dt=0.1,integrator=ImplicitMidpoint())\n",
"println(\"Orbit Period: \",stat.tau_p*1e6,\" microseconds\")\n",
"println(\"Time to orbit apogee: \",stat.tm*1e6, \" microseconds\")"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Orbit Period: 150.91323657094762 microseconds\n",
"Time to orbit apogee: 73.39193790434503 microseconds\n"
]
}
],
"source": [
"gcp = GCParticle(80.0,-0.515,1.8,-0.1)\n",
"path, stat = get_orbit(M,gcp,tmax=500, dt=0.01,integrator=ImplicitMidpoint())\n",
"println(\"Orbit Period: \",stat.tau_p*1e6,\" microseconds\")\n",
"println(\"Time to orbit apogee: \",stat.tm*1e6, \" microseconds\")"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Orbit Period: 150.91325981383784 microseconds\n",
"Time to orbit apogee: 73.39132137042135 microseconds\n"
]
}
],
"source": [
"gcp = GCParticle(80.0,-0.515,1.8,-0.1)\n",
"path, stat = get_orbit(M,gcp,tmax=500, dt=0.001,integrator=ImplicitMidpoint())\n",
"println(\"Orbit Period: \",stat.tau_p*1e6,\" microseconds\")\n",
"println(\"Time to orbit apogee: \",stat.tm*1e6, \" microseconds\")"
]
},
{
"cell_type": "code",
"execution_count": 29,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPwAAAGgCAYAAACOkApfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJzsnXd8lFX2/z/PlMykzEzKkJ6QRkknhV4tVAUEBRFBWQHFrqz6W1ZZwbKsna9rQ1yBtSAIggosTamGloSQQAgQ0nufkjL1+f0xmZHJTEibZ+p9v17z0rlPO0PmM/fec889h6JpmgaBQHAJWLY2gEAgWA8ieALBhSCCJxBcCCJ4AsGFIIInEFwIIngCwYUggicQXAgieALBhSCCJxBcCCJ4AsGFIIInEFwIjq0NsDRarRZVVVUQCASgKMrW5hAIjEDTNGQyGYKDg8Fi9b7fdjrBV1VVISwszNZmEAhWoby8HKGhob0+3+kELxAIAOj+IYRCoY2tIRCYQSqVIiwszPB97y1OJ3j9MF4oFBLBE5yevk5bidOOQHAhiOAJBBeCCJ5AcCGI4AkEF4IInkBwIYjgCQQXggieQHAhiOAJBBeCCJ5AcCGI4AkEF4IInkBwIZwult4RkEgkuHTpElpaWro9h8fjISEhASEhIVa0jODsEMFbGbVaje+++w719fU9npuZmYlly5YR0RMsBhnSW5kbN270SuyA7sfhl19+gUajYdgqgqtABG9lrl+/3qfz6+rqkJGRwZA1BFeDDOmtCE3TuHHjhlFbdHQ0xGKxUVthYSEaGxsN70+cOIHY2FiT8wiEvkIEb0UqKyvR2tpq1DZr1iz4+vqanPef//wHNE0DADQaDfbt24dHH32U5OkjDAgypLciXYfzYrHYROwAEBISgtGjRxu1lZaWIjs7m1H7CM4PEbwV6Sr4oUOHdnvuHXfcAW9vb6O2I0eOQCaTMWIbwTUggrcSTU1NqK2tNWq7neDd3Nxw7733GrUpFAocOHCAEfsIrgERvJXIzMw0eu/u7t5jOu3o6GgkJycbtRUUFOD33383zO8JhL5ABG8FlEolLl68aNSWnJzcqwIC06ZNg4eHh1HbqVOn8OOPP0KpVFrUToLzQwRvBXJzc9HR0WHUNnLkyF5d6+HhgZkzZ5q0X716FVu2bIFUKrWIjQTXgAieYWiaxvnz543ahgwZYtY73x0JCQmYNWuWyZJcTU0NNm/ejMrKSovYSnB+iOAZprCw0CSUtuuSW28YOXIklixZAj6fb9Qul8uxdetWXL58eUB2ElwDIngGUalUOHjwoFGbn58foqKi+nW/qKgoLF++3GR0oFarsXv3bhw7dow48wi3hQieQU6dOoWmpiajtnHjxg0oWk4sFmPFihWIjIw0OXby5Ens2rULKpWq3/cnODdE8AxRX1+PP/74w6gtNDQUKSkpA763u7s7Hn74YaSlpZkcy8/PJ848QrcQwTOATCbDjz/+CK1Wa2hjsVi49957LRYLz2azcc8992DmzJkm96yursbmzZtRVVVlkWcRnAcieAvT2NiIr7/+2sRRN3bsWAQEBFj0WRRFYdSoUXj44YfB4/GMjsnlcmzZsgVXrlyx6DMJjg0RvAUpLi7G119/bZK6ytfXF5MnT2bsudHR0VixYoVZZ96uXbtw/Phx4swjACCCtwhyuRw//fQT/vvf/6Ktrc3omK+vL5YuXQoul8uoDXpnXkREhMmxEydOYPfu3cSZRwBFO9lPv1QqhUgkgkQigVAoZOw5Wq0WZWVlyM/PR25uLhQKhck5QUFBePjhh+Hp6cmYHV3RaDQ4cOCA2a20wcHBWLRoEQQCgdXsITBDf7/nRPC9QKvVoq2tDbW1tYZXUVER5HJ5t9dERUVh4cKFJnNra6CP7jt06JDJUF4gEGDRokUIDg62ul0Ey0EE30lf/iE0Gg1KSkpQVFQEmUwGpVIJhUIBhUJh9P99GQrzeDzcddddSEtL69XmGCYpLCzErl27TEYfHA4H8+bNQ1xcnI0sIwwUIvhOevqHoGkaeXl5KCgowM2bNy264ywpKQlTp06Fl5eXxe45UOrr67F9+3Y0NzebHJsyZQomTZpE0mY5IETwnfT0D7F9+/Y+Z469HQEBAYiLi0N8fDz8/Pwsdl9L0tbWhp07d6K0tNTkWEJCAubMmcO4U5FgWforeJdKYimRSAYkdg6Hg4CAAPj7+yMwMBBRUVEOkUnWw8MDS5cuxf79+0325V++fBlNTU3EmeciuJTg29vbzbaHhIQgPDwcfD4fPB4Pbm5u4PF4Jv8vEAhsPi/vL2w2G7Nnz4a/vz8OHz5s5MyrqqrCV199hUWLFiEoKMiGVhKYxqUEb06sf/3rX+1qzs0kFEVhzJgx8PPzw65du4z8F1KpFFu2bMG8efMQGxtrQysJTOKY3VU/YbPZJm1d00e5AkOGDMHy5ctNsuKqVCrs3LkTJ0+eJJF5TopLCd5cD++qddv8/f2xcuVKhIeHmxw7duwY9uzZA7VabQPLCEzi8oK/dUebq+Hh4YFHHnkEI0aMMDmWl5eHrVu3GpW8Ijg+RPAuLHhAN82ZM2cOpk2bZnKssrISn376Kfbv33/bqEKC4+DyTjtXHdLfCkVRGDt2LPz8/LB7924jZx5N08jMzEROTg4SEhIwevRoDBo0aMDPc9TVDkfHpQRvzmnn6j38rQwdOhTLly/H9u3bTbb4qtVq5OTkICcnZ8DPoSgKAoEAPj4+8Pb2NvmvQCAg0X8M4VKCJ0P6nvH398cTTzyB06dP49y5c4w47miahlQqhVQqNRv9x2azjX4Auv4ouLu7W9wmV8GlBE/oHXw+H3fffTdGjRqF48eP49KlS1b9YdRoNGhsbOzWYcjj8Ux+DG59T8KEu4cIntAtQqEQc+bMwZ133onMzExkZWXZhfNOoVCgpqYGNTU1Zo97eXl1OzoQCoUu7T8ggif0iJeXF6ZMmYLJkyejpaVlwL29UqlES0sLmpub0dzcjJaWFsPLElMIuVwOuVyO8vJyk2MsFgsikcjwY+Dr64u4uLg+VQJyZIjgCb2Goij4+PhY5F7mYvZpmoZcLjf7Y9Dc3AypVDrgCECtVmu4t55jx47h7rvvxpgxY5zeWUgET7Ab9N57gUBgtpS2RqOBRCIx/AB0/W/XfIK9RavV4vDhwygrK8PcuXNNynk5E0TwBIeBzWbD19e32+F3d1MF/f/3lOykoKAAtbW1WLRoEfz9/Zn4CDaHUcGfPHkS7733HrKyslBdXY09e/bgvvvuu+01J06cwOrVq3HlyhUEBwfjlVdewapVq5g0k+AkuLm5wd/f36xYaZpGe3u70Y9BXV0d8vLyjM5rbm7GV1995bS7Bhl1V7a2tiI5ORmffPJJr84vLi7GrFmzMHHiRFy8eBF///vf8dxzz2H37t1MmklwASiKgoeHB0JCQpCQkIAJEyZg/vz5ePjhh03W9Z151yCjPfzMmTMxc+bMXp//xRdfIDw8HBs3bgQAxMbGIjMzE++//z7uv/9+pswkuDAxMTF44okn8OOPP6KystLo2LFjxyAUCs1uLnJU7GpB8syZMyabOKZPn47MzMxuM8cqFApD1Jb+RSD0BZFIhGXLliE5Odnk2P/+9z+zCUAdFbsSfE1NjUn9tYCAAKjVajQ0NJi9ZsOGDRCJRIaXOe8ugdATHA4Hc+fONelwlEol9uzZ4zQh2HYleAAm66D6OVR366Nr1qyBRCIxvMwFWxAIvUG/a7DrEL68vBwZGRk2ssqy2JXgAwMDTcIl6+rqwOFwuk0BzePxIBQKjV4EwkCYMWOGSYBRRkaGU2QAsivBjx07FkeOHDFqO3z4MNLT08mGCILV4PF4mD9/vlFbe3s78vPzbWSR5WBU8HK53GgPdXFxMXJyclBWVgZANxx/5JFHDOevWrUKpaWlWL16Na5evYqvv/4a//nPf/DSSy8xaSaBYEJoaCgiIyON2i5cuGAjaywHo4LPzMxESkoKUlJSAACrV69GSkoK/vGPfwAAqqurDeIHgMjISBw4cADHjx/HiBEj8Oabb+Ljjz8mS3IEmzBy5Eij9xUVFd3u0HMUGF2HnzJlym0DF7Zu3WrSNnnyZLOljgkEazNs2DAIBALIZDJD24ULFzB79mwbWjUw7GoOTyDYEywWC6mpqUZteXl56OjosJFFA4cInkC4DampqUZLwiqVCrm5uTa0aGAQwRMIt0EoFGL48OFGbZmZmQ4bY08ETyD0QHp6utH7+vp6s8k3HQEieAKhByIjI00CvzIzM21kzcAggicQeoCiKJNe/urVq3aR0LOvEMETCL0gOTkZHM6fq9hardYhl4+J4AmEXuDu7o6EhASjtqysLIfbRUcETyD0kq6Rd1KpFDdu3LCRNf2DCJ5A6CXBwcEIDg42anO0+HoieAKhD3R13t28eRNNTU02sqbvEMETCH0gISHBJG99VlaWjazpO0TwBEIf4HK5JhlxLl686DDJMYjgCYQ+0nVY397ejitXrtjImr5BBE8g9BE/Pz9ERUUZtTlK5B0RPIHQD7r28hUVFaiurraRNb2HCJ5A6Af65Bi34gi9PBE8gdAPHDU5BhE8gdBP0tLSTJJjXLp0yYYW9QwRPIHQTwQCgcMlxyCCJxAGQNf4+oaGBrtOjkEETyAMgIiICIdKjkEETyAMgO6SY9ya2tqeIIInEAbIiBEjTJJjXLx40YYWdQ8RPIEwQPh8PhITE43a7DU5BhE8gWABug7rpVIprl+/biNruocInkCwAOaSY9ij847R2nIE10GjpVHcIMeNWjkqW9pR0dyOBrkCbUoNWhVqaLQ0eFwW3NgsCPhcBIr4CBDyEe7rgeGBAoT6uBsFsTgiI0eOxM8//2x4r0+O4evra0OrjCGCJ/QLpVqLrNJmnLpRj/PFTcivlqJNqen3/QR8DpJDvTE22g/jY8RIDBGBzXKsH4D4+HgcOnTIKLw2MzMT06ZNs6FVxhDBE3qNUq3Fiev1+DmnEr8X1JkI3J3LxtBAAcJ83BHi444AAR9ePA48eGxwWBQUai0UKi0k7SrUSDtQI+nAzXo5btbLIetQ43RhA04XNuC9Q9cg9nLDzIQg3JMUhFERvmA5gPj1yTHOnj1raMvJycEdd9wBLpdrQ8v+hAie0CPlTW3YllGCXdkVaGlTGdr9PN0wcYgY42LEGBHmjehBXv3qlZVqLQrr5Dhf3IiMm404U9SIBrkS35wtxTdnSxHu64GlYwZjQXoovD3cLPnRLE56erqR4PXJMbpmybEVRPCEbsmtaMGnxwpxJL8W2s7wcH8BD7OTgzE7ORhJISKL9LxuHBbigoWICxZi2fhIqDRa/FHYgH251Th0pQZlTW14+8BVfHDkGhaNDMdTU6LhL+T3fGMboE+OUVRUZGjLzMwkgifYL9dqZPjwyDUculJraJs4RIzHxkdi0tBBjM+tuWwWpgzzx5Rh/nhzbgJ+zqnEtjOluFotxdaMEmw/X4alYwbj6Tti4ONpfz1+enq6keArKytRXV2NoKAgG1qlgwieYEDSpsL7h6/h23OloGmARQH3pYTgqSnRiPEX9HwDBnB3Y2PRqHA8ODIMfxQ24sMj15Bd1oKvThdjV3YFXp4+DItGhtuVg0+fHOPW8NoLFy5gzpw5NrRKB1mHJ4Cmaey5WIG7PjyOb87qxD4zIRCHXpiEDxeOsJnYb4WiKEwYIsbuJ8dh22OjMDxQgJY2FV7dcxn3ffoHrlZLbW2iARaLhbS0NKM2e0mOQQTv4jS1KvHEN1l4ccclNMiViPH3wvaVY/D5kjQMCbC90LtCURQmDx2Efc9OwOuz4yDgcZBXKcHcT/7AlydvQqO1j73oqampYLH+lJdarbaL5BhE8C7Miev1mL7xJA7n14LLpvDy9GE48NxEjI326/liG8Nhs/CX8ZH47aXJuGu4P5QaLf55oACLN59Fncz2Pam9JscggndBaJrGp8cKsWzLedTLFBji74W9T4/H03fEwI3jWF8JfwEfXz2ajn/NT4SHGxvnipsw+9+nkV3WbGvTTOLrGxoaUFJSYhtjOnGsvy5hwLQq1Hj6+2y8d+gaaBp4aFQ4fn12AuKDRbY2rd9QFIVFo8Kx79kJGOLvhVqpAg9uOoMdF8psaldERATEYrFRm63j64ngXYhGuQIPbT6LA3k14LIp/HNeIjbMTwSfy7a1aRYhapAX9jw9HjPiA6HS0Ph/u/Pw0ZHrNhtGm0uOUVBQYNPkGETwLkJlSzsWbDqD3AoJfDy42L5yDBaPDre1WRbHi8fB50tS8dydMQCA//vtBv6+57LNnHnJyclGYbVarRbZ2dk2sQUggncJShpacf9nGSiqb0WwiI8fV41DeoT97OCyNBRFYfW0YXjzvgRQFLD9fBme/+Ei1BrrJ6Tg8/lISEgwarNlcgwieCenorkNizefRY20AzH+Xtj91DjE+HvZ2iyrsHTMYHy2OBVcNoV9udX464+XbNLTd81sK5PJbJYcgwjeiamRdGDx5nOoknQgapAntq8cgyCRu63NsiozE4PwyeJUcFgUfs6pwsu7LkFrZdEHBQUhJCTEqO3ChQtWtUEPEbyTIutQ4dGvz6OsqQ3hvh74fsUYDBLwbG2WTZgeH4iPH0oBm0Xhp+xKbPjfVavb0NV5V1RUhMbGRqvbQQTvhKg1Wjzz/UVcq5VhkICH71aMRqDIPneXWYtZiUH4YEEyAGDzqWL890yJVZ8fHx8Pd3fj0ZUtluiI4J2QN/bl48T1evC5LPzn0XSE+XrY2iS74L6UELw8fRgAYN0vV3Akv7aHKyyHPjnGreTk5EClUnVzBTMQwTsZu7Iq8N8zpaAoYOODKUgK9ba1SXbFU1OisWhkGLQ08MIPF1FYJ7fas7sO6zs6OnDlyhWrPR8ggncqrtfK8NrePADAC3cNxYyEQBtbZH9QFIU370vAmChftCo1WPVtFloVaqs829fXF9HR0UZt1h7WE8E7CW1KNZ76LhsdKi0mDhHj2c7AE4IpXDYL/34oFQFCHgrr5Ph/u3OtFo3XtZevrKxEVVWVVZ4NEME7DRsOFKCwTo4AIQ8bHxzhEEkfbckgAQ+fPaxbrtuXW43vz1sn7n7o0KEQCoVGbdbs5a0i+M8++wyRkZHg8/lIS0vDqVOnuj1369atoCjK5GUPyQPslYzCBnxzVlei+MOFI+Dn5ZrLb30lbbAv/jZTt4X1rX1XUdLQyvgzu0uO0d7ezvizASsIfseOHXjhhRfw6quv4uLFi5g4cSJmzpyJsrLuf1GFQiGqq6uNXny+ay8rdUerQo1XducCAB4eHY7xMeIeriDcymPjIzE2yg/tKg1e3JljlfBbWybHYFzwH374IZYvX44VK1YgNjYWGzduRFhYGD7//PNur6EoCoGBgUYvgnk+PHIdFc3tCPF2x5pZsbY2x+FgsSi8vzAZAh4HF8ta8MWJm4w/08vLC7Gxxn8rayXHYFTwSqUSWVlZJpU3pk2bhoyMjG6vk8vlGDx4MEJDQ3HvvffetvSuQqGAVCo1erkKN2pl2JZRAgB4e14CvHgkJ2l/CPF2x/q58QCAj38vRLEVhvZdnXeNjY1WSY7BqOAbGhqg0WgQEBBg1B4QEICamhqz1wwfPhxbt27FL7/8gu3bt4PP52P8+PG4ceOG2fM3bNgAkUhkeIWFhVn8c9gjNE1j/a/5UGtpTI0LwJRh/rY2yaGZlxKCiUPEUKq1WLv3MuO97eDBgzFo0CCjNmvE11vFade1SCBN090WDhwzZgyWLFmC5ORkTJw4ETt37sTQoUPx73//2+z5a9asgUQiMbzKy8stbr89ciS/FqcLG+DGYWHtPXG2NsfhoSgKb85NgBuHhdOFDfjlErNLZbZKjsGo4MViMdhstklvXldXZ9LrdweLxcLIkSO77eF5PB6EQqHRy9nRaGm8f/gaAGD5hEiE+5HQWUsQIfbEs3fo4hfe3JcPWQezYa9JSUlGyTFomkZWVhajz2RU8G5ubkhLS8ORI0eM2o8cOYJx48b16h40TSMnJ8cuqnbYC/tyq3C9Vg4Bn4NVk6J7voDQax6fHIVIsSca5EpsOlHU8wUDgM/nIzEx0agtOzsbGk3/q/D2BOND+tWrV+Orr77C119/jatXr+LFF19EWVkZVq1aBQB45JFHsGbNGsP569evx6FDh1BUVIScnBwsX74cOTk5hvNdHbVGi/87qhvtPD4xCiIP+6hK6izwOGz8vxm6tfnNp4pQLWF2fdzayTEYd+s++OCDaGxsxBtvvIHq6mokJCTgwIEDGDx4MACgrKzMaE2ypaUFjz/+OGpqaiASiZCSkoKTJ09i1KhRTJvqEOzPq0ZRQyt8PLj4y4RIW5vjlEyPD8DICB9cKGnGB4ev4/3ObbVMEBgYiNDQUFRUVBjaMjMzTZbtLIVV1nGeeuopPPXUU2aPHT9+3Oj9Rx99hI8++sgKVjkeNE1j8yndMPMv4yPJMhxDUBSFv8+KxbzPMrA7uwKPjY9EXDBzvqH09HQjweuTY/j5Wb4gCImldyDOFzfhcqUUPA4LS8YMtrU5Tk1KuA/uSQoCTQMf/2beYWwprJkcgwjegfjqdDEA4P60UPjaYZlkZ+P5u4aAooCDV2pQUMNcQBeHw0FKSopRW15eHiOxAETwDkJVSzuOXtVlaHlsPJm7W4OhAQLMStCtDv37t0JGn5WUlGT0vrW1lRFvPRG8g/BTdgVoGhgd6esyaabtgWfv0q3LH7hcjeu1zAXF8HimOxxJD++iaLU0dmbqnDoL010jdNheGB4oxIz4QNA08Nkx5np5c5GnRPAuyvmSJpQ1tcGLx8HMRLJz0No805k9aF9uNWokzORlIIInGNDHdd+TGAQPN7IUZ20SQkQYHekLtZbGtjMljDyDCJ4AQBc3f/iKbi/CrCQSXmwrlncGOX1/rgxtSssnvSSCJwAAMkua0CBXQsjnYGyU5QMxCL3jrtgADPbzgKRdhd1ZFT1f0EeI4AkAgP9d1vXud8cFwI1D/ly2gs2i8JdxEQCAr/8osXh9OiJ4AgDg2LU6ALr6aATbsiA9DAIeB8UNrThbZNm6cETwBJQ3taG0sQ1sFoVx0WQ4b2s8eRzMGREMABZPa00ET8DpwgYAQEqYNwR8sg3WHnhoVDgA4PCVWjS1Ki12XyJ4gkHwJPW0/ZAQIkJiiAhKjRY/ZVvOeUcE7+LQNI1znfPECUOI4O0JfS///fkyi4mSCN7FqWhuR4NcCS6bQmKIyNbmEG5hzohgeLixUVTfiszSZovckwjexckpbwEAxAYJweeybWwN4Va8eBzMStQFQe25WGmRexLBuzh6wY8II/Xd7ZH7RoQAAA7kVUOpHnh5KnOC12otX/aKCN5Oya3QCT45lAjeHhkb7Qd/AQ8tbSqcuF4/4PtZq4cnOzHsEJqmcb1WDgCM5lKjaRqljW04X9KEazUylDW1oV6mgFKtBQ3A250LsYCHof5eiAsWIj3CFyJ3sjwI6CLv5iQH46vTxdh7sRJT43pXZ6EvEMG7CPVyBSTtKrAoIFLsafH7X6+V4afsSvx6qQqVLb1Pw8xhUUiP8MGc5BDMHREMTxdPonlfSgi+Ol2Mo1drIetQDShWgvTwLkxhZ+8e7uthUYfd+eImfHa8EMev/TkEdWOzkBQqQnKYNyL8PBAg5IPHZYOmaUjaVaiRdKCgRoZL5S0oamjF2aImnC1qwj8PXMVDo8KwanK0y9ajjw8WInqQJ27Wt+LQlVo8kBY6oPtRFGUkciJ4F+FGnU7wMf4Ci9yvvKkNb+3Px6Erupx4LEq3+2teSgjuGOYPd7fe/aiUNrbi0JUa/HC+HEUNrdh8qhjfnyvDU3fE4PFJUeCyXcslRFEUZicHY+PRGzh4uYYIntA/Shp15YqjBg1sOE/TNL47V4a39uejQ6UFm0VhYXoYVk2OwmC/vt97sJ8nHp8UjZUTo3D8ej0+OHwNlyuleO/QNezLrcZ7DyQhwcViBqbHB2Lj0Rs4daMerQr1gKY55oquWhrX+kl2EPRplIJF/H7fQ9qhwhPfZOG1vZfRodJidKQvDjw3ERvmJ/ZL7LdCURTuGOaPX5+ZgA8XJsPHg4ur1VLM/zwDOy+4RvVePcMDBRjs5wGFWjtgbz0RvItS1Sn4IG/3Hs40T2VLOx74PAOH82vBZVN47Z5YbF85BsMCLTNF0ENRFOanhuLI6sm4O9YfSrUWr+zOxes/X7b4fnF7haIow9blQ1dqeji753vdChG8i1Dd6TkP6kcPf6NWhnmf/oHrtXL4C3jY/eQ4rJgYBRbL1AtsKcRePHy5NB0vTRsKigK2nSnF6p05UGksHzhij0yP1y3J/X61bkBBOETwLohao0W9XAEACBL1rYcvrJPjoc3nUCdTYFiAAHufHo8kKwXusFgUnrlzCDY+OAIcFoW9OVV4cUeOS/T0KWE+GCTgQaZQDygxBhG8CyLrUEP/d/bpQynoiuY2LN58Fg1yBWKDhNjxxBgE93NKMBDmjgjBpqVp4LIp7Mutxhv78hn54toTLBaFO4YNAgCcHMA8/tYqygARvEsgaVcBADzd2OD0cplLrlBjxbZM1MkUGBrghe9WjIa3h+1qz90VG2Aosbw1owT/PVNqM1usxaShOsEPxHFHengXRNqhE7ywlyGsWi2NF37IQUGNDGIvHrb+ZZRdFJqcOyIEa2YOBwC8tT/fsBnIWZkQIwaL0sVQVPUhevFWiOBdEGm7Lud5b2PW/9MZ2unGYWHzI2k2GcZ3x+OTojAjPhAqDY2nv8uGXGH5fO72greHm8FfcupG/3p5IngXpLWzyEFvot8uV0rw7qECAMDrs+OQEu7DqG19haIovLsgCWG+7qhsace7BwtsbRKj6If1J6839Ot6IngXRNPp1eaybv+nUag1eGFHDlQaGtPjA7C4M+2SvSHkc7Fhnq4U8jdnS5FZ0mRji5hjYmcqsnPFjf0SKxG8C6IXPLuHdfMvTxShsE4OsRcP/5qfZHa3lb0wYYgYC9JCQdPAm07stU8KFcGNw0KDXImSxrY+X08E74LoBc9hdy/gssY2fNJZunjtvbHwsQMnXU+8MmM4PNzYuFQhMVTTcTbV8RNQAAAgAElEQVR4HDaSQ3V7CS70YyRDBO+C6AXPuk2P/eb+fCjUWoyP8cOc5GBrmTYgBgl4WDkxCgDw/qFrThuQkx7hCwD9mroQwbsg+p69u7DUrNJmHMmvBYsC1s+Jt+uhfFdWToqCkM9BUUMrfiuos7U5jDAyQuc4zSwZeDZbIngXQJ/wokOlMTlG0zTe6/TKP5AWarH98tbCi8fBQ6N1zsX/nC6ysTXMMCJMJ/iihlbIOmMqeos1fryJ4O0MXmeF2A6VaQ+fcbMRZ4ua4MZm4fm7h1rbNIuwbFwEOCwKZ4uaUFgns7U5FsfX082w6amgxv4+HxG8nWHo4dWmPfymk7pe8aFRYQixowCbvhAkcsfkzvXqn3OqbGwNM8QF6RKP5ldJbWyJKUTwdoZXZ8YUeYdxVNq1GhlOXq8HiwKWT4iyhWkWQ1+B9eecKqdcoosPJoIn9BI/L90SW1Or0kgMX53S9e7T4wMR7udhE9ssxbS4QLhz2ShrasPVavsb9g6U2M4evqDW/j4bEbydod/4otbShrh6SZsKP1/SDX9XTHTs3h3QhQ2PidItX/U37tyeiehMLV7amZvQniCCtzN4HDYEncP6hlZdIoxfLlVCqdZieKAAqeHOUYnGEHfuhIIP99WNwFraVIbtzvYCEbwdIhbo8rzXSXWC35mpq0O+ID3Modbdb8e4aF3c+cWyFkOwkbPgyeNgUOffsKwfIbZMQgRvh+h7iNLGVlytliKvUgIum8J9Ixwjqq43xPh7wZ3LRptSg+IGua3NsTiD9X/DJvsa1hPB2yH68lLFja3YnaXr3e8aHuBUFV7YLMpQN+9ypf15swdKgFC3Ft8gU9jYEmOI4O2QiE4v/M26Vvyaq3PWzU8NsaVJjDC8M212YZ3z9fB652tjq9LGlhhDKs/YIUMCdEI4elVXGkrI52ByZ5JEZyLUR/fD1peClo6CfnnV3gTvUj28VmsarqpQ2NeQCwCSw7yN9sPPSgwCj2O5opL2QqiPLlqwstmJBN8hASSV8NP38PJbvl+SSt3xXuKwm2c+++wzREZGgs/nIy0tDadOnbrt+bt370ZcXBx4PB7i4uKwZ88ei9jRNQ0wAPB49jcv9uJxkHhLjbYF6QMrUmhL2lRtSNyWiMRtiWhTGXusAztjzuvl9vej2y86JMC39wNbZ8FdrfNLtOv3REgqgK2zdMf7IHpLw7jgd+zYgRdeeAGvvvoqLl68iIkTJ2LmzJkoKysze/6ZM2fw4IMPYunSpbh06RKWLl2KhQsX4ty5c0ybalf8ddpQCPkcLEwPRdpgX1ubwwiebp1hxM6S3FIhB1rrgeYSuP3xHgBAqdZ0iv0eoLlEd1xhO58F44L/8MMPsXz5cqxYsQKxsbHYuHEjwsLC8Pnnn5s9f+PGjZg6dSrWrFmD4cOHY82aNbjrrruwceNGpk21KyYOGYTcddPx7gPJtjalX7Sp2tCmakO7+s/heru63dAOAJ483TSl1VkELwoBlu0HfCLAa9U5W5Wt0j/F7hOhOy6ynQOWUaedUqlEVlYW/va3vxm1T5s2DRkZGWavOXPmDF588UWjtunTp3creIVCYTQPl0qdb4nHERn9/WiTtik7pxj+P+/RPMPOQMUA6rHZHaJQYNl+uH3xItAMKOtuALySW8Ru2+kZoz18Q0MDNBoNAgICjNoDAgJQU2M+r1lNTU2fzt+wYQNEIpHhFRYWZhnjCYyjz+rDvU3+PodEFAr1mGcBAGx0bnOe96XNxQ5YaVnOXK6u24WI9uX8NWvWYPXq1Yb3UqmUiN4OOLdY53NpV7cbevbjC4/DnfPnPn6VpjMldy9LajkMkgooz2wCsBg8dMbS73nc+Xt4sVgMNptt0jvX1dWZ9OJ6AgMD+3Q+j8eDUCg0ejkDsg4VPj9+E5crbefRHQgeXA94cD2MBO7OcTe0A0C7Utf76bP8OAWdDjplqy6nHS8kSTecby7RzeUlFTY1j9F/aTc3N6SlpeHIkSNG7UeOHMG4cePMXjN27FiT8w8fPtzt+c7K2/uv4p2DBXj8v5lOt7lET0PncpzYWUKGJZUGB53CIwgA4OblY3Dk/Sn6SpuZyPiQfvXq1Vi6dCnS09MxduxYfPnllygrK8OqVasAAI888ghCQkKwYcMGAMDzzz+PSZMm4Z133sHcuXPx888/4+jRozh9+jTTptoNWi2NvTm6L0WVpAMnb9TjjmH+Nraqf3hwPZD3aJ7ZY/Wdceb6nWUOD88L8NRFRDYPexo4Xq2r4tvpyMPWe3THeV42M5FxwT/44INobGzEG2+8gerqaiQkJODAgQMYPHgwAKCsrMwoIGbcuHH44Ycf8Nprr2Ht2rWIjo7Gjh07MHq0qdfXWSlqkBslsfzhfJnDCv52VEt0S3ZOI3i+CFiyG1DI0XRaNxXTh9jqRH9AJ3a+6DY3YRarOO2eeuopPPXUU2aPHT9+3KTtgQcewAMPPMCwVfZLXue8XeTOhaRdhd+u1qFO1gF/Ad/GllmWG52bZqIH2a7Hszh8EcAXoVGuS+xhVLrbhuvvepzIW+I8FNfr9lDPSgxCarg31Foau7NsN++zNBotjTM3Gw3VWWL8nUjwnejDhf3srAwY2S1nhxR3ZkmJFHsgJcwb2WUt2HGhDKsmRzl8xpuDl6ux/td8VEs6DG2v7c0DTdOYkRBkQ8ssS0XnhqAQH/tKJ056eDtEn/xwsJ8n7k0OghePg5LGNpwparSxZQPj4OVqPPlttpHYAaBepsST32bj4OVqG1lmWdQaLcqbdD/aEX6eNrbGGCJ4O6SmUxAh3u7wcONgbmdqq+/Pmd9w5AhotDTW/5qP2y0wrv813ymWIKtaOqDW0nDjsBAo7L/fhYnRHBG8naHV0oakCfr16cWd9dgOXakxrF07GueLm0x69luhAVRLOnC+uO9VV+2Nks4RWrivB1is3ou26/53IngXQNKuMvRyeg9vfLAIyWHeUGlo7LhQbkvz+k2drHux9+c8e+ZaZ0256EF9G84Twbsg+h5c5M6F2y0hp4+M0cUtbPmjxGxlWXunt0uKzrD0eKVKt6yaENy39XYieBdE2lliWOTONWqfMyIYwSI+GuQK7MqybTx2fxgV6YsgER/dfYUpAEEiPkZFOn6yjyudNeXiQ/q2r4MI3gXRR9h5uBnnsOOyWVg5SVdm6suTRYatpY4Cm0Xh9dlxZp12+q/167PjjHL5OSLtSg1u1usCiuJJD0/oCcMOMq5p0soHR4bBz9MNZU1t+MEB5/IzEoIwNsq0Bw8U8fH5klSnWIfPr5ZCSwNiLzf49zFk2BqCJ4E3doa+Lrw71/S32MONg+fvHoJ//HwF/3f0OualhBjKSzsClS3tyCptAQC8MSceIg8u/AW6Ybyj9+x69NGDqeE+fRYs6eFdEEXnkN6tm7TUD40KR6TYEw1yJT4/XmhN0wbMx0dvQKnRYly0Hx4ZF4G5I0IwNtrPacQOwLCs2B9fBBG8C6L/G3eXk5zLZuFvM4cD0M3lr9thDXJzFNbJsStb52x8afowG1vDDFotjcxSXeKLkREDFzwTEMHbGfreTnubP/60uADcHesPlYbG33bnQmvn0Wk0TeO1vXnQaGncHRuA1HAfW5vECNfrZJC0q+DhxkZ8cN8zL5Ee3gVhdf6RbxdiSlEU3rwvAV48DrLLWvDV6SJrmdcvdmdX4mxRE/hcFl6fHWdrcxjjXJFuOJ822AecfuTpI4J3QQw9fA+rbkEid/x9ViwA4N2D15BT3sK0af2ivKkN63+9AgB44e6hCOsso+yMnLqh2wM/JsqvX9cTwbsg/E7vvN5bfzseGhWGexKDoNbSeOb7bDTbWeFClUaLZ7dfhKxDjdRwbyyfEGlrkxhDodYg46ZuN+OUfhb+JIJ3QQR8XYSdtF3V47kURWHD/YkI9/VARXM7Hv8m027Cbmmaxj9+voyc8hYI+Rx8/FCK86WjvoXMkma0KTUYJOAhLqh/mZOJ4F0QoV7wHb0rvyTkc/GfR9Mh4HNwoaQZL/14yS62mH52/Ca2ny8HiwI+enCEoTS0s3L8Wh0AYPLQQf0WKhG8CyLg6wJpZB2qXi/TDAkQYNOSNHBYFPblVuOvO3OgtmHo7ZY/ivHeoWsAgNdnx+OuWPM1BZyJE9d18/f+DucBIniXRL9pRqWh0ars/fB8XIwY/34oBRwWhb05VXh+R47Vh/c0TeOz44VY/2s+AODJKdF4dFyEVW2wBRXNbbheKweLAibEiPt9HyJ4F8STxzGEy9ZK+7Y3fGZiED5ZnAoOi8L+3Go8tPms1faXK9QavLIrF+8e1PXsz901BK84aYBNVw5e1lVKGhnhq8tD30+I4F2UQJFuT3jNbTLEdMeMhED897FRELlzcbGsBfd+fNow3GSKm/VyLPziDH7MqgCLAl67Jxarpw51+ISbveXQFZ3gZyYEWvS+RPAuQlCn4G+XEup2jIsRY+/T4xHj74U6mQKPfn0ef9+ThyYLL9sp1BpsOnETs/7vFC5VSCDkc7D1L6OwYmKURZ9jz9TJOgzhtNMHIHhz/hoieBfBIPiW9n7fI1LsiV+fmYBHx+oy5Xx/rgyT3zuGz44XQtKLJb/boVRr8VN2Be7+8AQ2/K8ACrUWE4eIcfCFSZg0tP9OK0fk0JVa0DQwIswbQaL+p6S2luAdZ2+lCxHWuYRV0pmfvr+4u7Gxfm4CZiYG4Y1f85FfLcW7B6/h098LMS81BLOTgpEe0butqTRN42q1DP+7XI0fLpQb6sL5C3h4afowLEgLdZkh/K0c6py/zxjgcJ4I3oXRV2IprLPMTrgxUX749dkJ2HuxEl+eLMK1Whm+PVuGb8+WwduDi9RwHySGiBDq445BAh44LBa0NI3GVgWqWjqQXyXFxbJmVN0yxRgk4GHZuAj8ZXwEPNxc82vU1Ko01AqYEU8ET+gnesHfrG8FTdMW+cOzWRTuTwvF/NQQ/FHYiL05lTh8pQYtbSr8XlCH3wvqerwHn8vChJhBmDMiGDPiA42SbLoi+3OroNHSiA8WIkI8sIITRPAuzGA/T7BZFOQKNWqkHQOaG3aFoihMGCLGhCFiqOYn4kqVFNmlzbheK0OVpAMNMoVha66flxv8BXwMCxQgMUSE1HAfuLuZT8zhiuy5qKv3Ny9l4EUiieBdGDcOC5FiTxTWyVFQLbOo4G+Fy2ZhRJg3RoR5M3J/Z6a0sRXZZS1gUcCc5OAB34946V2cpBBdxlN73fbq6uy9WAUAGB8jhv8AyknpIYJ3cZI7e91LFUTw9gZN09hzUZeua36qZWq+E8G7OPph9qXyFqvkOiP0nuyyFpQ0tsGdy8a0OMtE1xHBuzjDgwRwY7PQ3KZCcUOrrc0h3ML287oqvjMTA+FpoTThRPAuDo/DRkq4rpfXZ1Ih2B5Juwr7cnXz94c7q/oyBRG8i6Hfann6RoONLSHo2XuxEh0qLYYFCCyafZf08ARMGKITfMbNBrvIYuPq0DSN78/phvMPjQpjPJRYre5d1qO+QARvxySFekPA50DaoSbeejsgu6wZ12pl4HFYmJcaatF7e3h4gMs1rhhcXV1t0WcARPB2DZtFGXafHb5Sa2NrCFv+KAEAzE4ONinnPVBYLBZCQoyX+CoqLF8WnAjeztEnVTh4uZosz9mQiuY2/K9zZ9xj45lJt00ET8Adw/zB47BQ0tiGghrHqCPnjGzLKIFGS2N8jB/i+lFGqjeEhhpPEyorK6HtqSJJHyGCt3M8eRzDsP5AnuXndISekXWo8MP5cgDAignMZfPpKnilUommpiaLPoMI3gG4NykIAPBTdiXx1tuAHRfKIVOoET3IE5MZzOhTV2e6RZnNtuzuRCJ4B2B6fCCEfA4qW9qRcZOsyVuTDpUGm0/pinWumBgFFoO17DMzM43ei8VieHtbdicjEbwDwOeycV/nnusdF8ptbI1rseNCOWqlCgSL+LjfwktxtyKVSlFQUGDUlp6ebvG1fiJ4B2FhehgA3fJcg1xhY2tcA4Vag8+P3wQAPHlHDKMZfrKysoxWYbhcLpKTky3+HCJ4ByEhRITkUBGUGi2+PVtqa3Ncgp0XyjszDvGxMJ253l2j0SArK8uoLSkpCXz+wPfZd4UI3oHQ53v/5kyp3VSJdVbalRp8cqwQAPDUlGjwOMyl9rp69SpaW413RI4cOZKRZxHBOxAzEwIR4u2OxlYl9nbmUyMww1enilArVSDUxx0LR4Yx+qwLFy4YvQ8PD0dAADMFOIngHQgOm4W/jI8AAHxx4iZUNqwQ68zUyxT44oRu7v7KjOGM9u61tbUoKyszamOqdweI4B2Oh0aFw8/TDSWNbfgp2/KhlwRg49HraFVqkBzmjdmdMRBM0bV39/T0RGxsLGPPY1Twzc3NWLp0KUQiEUQiEZYuXYqWltvv+poyZQooijJ6LVq0iEkzHQpPHgdPTokGAHz8WyEUajKXtyQFNVL80Ln0+eqsWEa3wCoUCuTm5hq1paamWjzY5lYYFfzixYuRk5ODgwcP4uDBg8jJycHSpUt7vG7lypWorq42vDZt2sSkmQ7HkjGDESDkobKlHd+dLev5AkKv0GppvLrnMjRaGjMTAjEq0pfR5126dAkq1Z91/iiKQlpaGqPPZCwv/dWrV3Hw4EGcPXsWo0ePBgBs3rwZY8eOxbVr1zBsWPe1wz08PBAYaNnSu84En8vGc3cNwat7LmPj0euYOyIYfl48W5vl8OzMLEdWaTM83dj4x+w4Rp9F07TJcH7YsGEQiUSMPpexHv7MmTMQiUQGsQPAmDFjIBKJkJGRcdtrv/vuO4jFYsTHx+Oll16CTNb9LjGFQgGpVGr0cgUWjQxHbJAQ0g413j98zdbmODyNcgU2/E8X6fbi1KGMFf/QU1JSgoYG4zBpJp11ehgTfE1NDfz9/U3a/f39UVNT0+11Dz/8MLZv347jx49j7dq12L17N+bPn9/t+Rs2bDD4CEQiEcLCmF1CsRfYLArr58QDAH64UI5ckhFnQLz+yxVI2lWIDRJi2bgIxp/XNW7ez88PkZHM7LO/lT4Lft26dSZOta4v/Ycx5/DoqTjiypUrcffddyMhIQGLFi3Crl27cPToUWRnZ5s9f82aNZBIJIZXebnrxJqPivTFfSOCQdPAK7tyoVSTZbr+8MulKuzLrQabReFf8xPBYTO7eCWVSnH16lWjNibi5s3R5zn8M88806PXPCIiArm5uaitNU3LVF9f36eggtTUVHC5XNy4cQOpqakmx3k8Hng8152/vnZvHE7eaEBBjQyf/H4Dq6d17xshmFIj6cDavZcBAM/cEWOo+MMk2dnZJnHzI0aMYPy5QD8ELxaLIRaLezxv7NixkEgkOH/+PEaNGgUAOHfuHCQSCcaNG9fr5125cgUqlQpBQcyuhzoqYi8e3pybgKe/z8anx29ialwgEkOZdfw4C1otjZd3XYKkXYWkUBGeuTOG8Weai5tPTExkJG7eHIyNXWJjYzFjxgysXLkSZ8+exdmzZ7Fy5Urce++9Bg99ZWUlhg8fjvPnzwMAbt68iTfeeAOZmZkoKSnBgQMHsGDBAqSkpGD8+PFMmerw3JMUhHuSgqDR0nj+h4uQdah6voiAT44V4tSNBvA4LHy4MBlchofyAFBQUAC5XG7UZg1nnR5GP+F3332HxMRETJs2DdOmTUNSUhK++eYbw3GVSoVr166hra0NAODm5obffvsN06dPx7Bhw/Dcc89h2rRpOHr0KKPBCM7AW3MTECzio6ihFX/7KY8kvOyBUzfq8dHR6wCAt+5LQIy/wCrP7boUFxYWZtUlaEbrw/v6+uLbb7/t9nhERITRFzMsLAwnTpxg0iSnxcfTDf9enIoHN53B/txqjIrwxaNW8DY7IlUt7Xj+hxzQNLBoZBgWpFtnZaeurg6lpcZbm63ZuwMklt6pSBvsgzWzdHHYb+7Lxx+FJB1WV2QdKjy29QKaWpVICBFiXefSpjXo2rt7eHgwGjdvDiJ4J+Ox8RGYOyIYai2NVd9m4UYtSW2tR6XR4unvL6KgRoZBAh42LU0Hn2udqWJ3cfMcDqODbBOI4J0MiqLwzv1JGBnhA1mHGsu2XECdtMPWZtkcmqbxj58v4+T1erhz2fj60ZEI8WY2mu5WcnNzoVQqDe8pikJ6errVnq+HCN4J4XPZ2LQ0HRF+Hqhsacfir865dB48mqbx5r6r2H6+HBQF/N+iEVZdujQXNz906FDG4+bNQQTvpPh6uuGb5aMRJOKjsE6OJV+dQ1OrsucLnQyapvHOwWv4+o9iAMA785MwLd66G7NKS0tRX19v1GZtZ50eIngnJszXA9+vHAN/AQ8FNTIs3nwWtS40vKdpGv86WGDIXvPmfQmMp6syR9fe3dfXF1FRzFWwuR1E8E5OpNgT368cg0Gdop//WQaK6uU9X+jgqDVavLIrF5tO6IpIrL03DkvHDLa6HTKZzCTf/MiRI60SN28OIngXIMbfCz89OQ6RYk9UtrTjgS/OILPEsjXL7Am5Qo1V32bhx6wKsCjg3QeSsHwC8zvRzJGVlWVUEJLD4TCSb763EMG7CGG+Hti1aiySQkVoalVi0Zdn8c2ZEqeLyCtpaMX8z/7A0at14HFY2LQ03VDEw9poNBqTXZ6JiYlwd7fe6kBXiOBdCD8vHn54fAzuSQqCWktj7c9X8NKPuWhVqG1tmkX4vaAWcz45jeu1cvgLeNj++BhMjWMm3XNvuHbtmknyFls56/QQwbsYHm4cfPJQCtbMHA4WBezOrsCsj08hq7TZ1qb1m3alBmv3XsZjWzMh7VAjNdwb+56dgNRwH5va1dVZFxoaavNdn0TwLghFUXhicjS+XaFbtittbMOCLzLwr/8VoE3pWL19VmkT7vn3KXzTWX5r2bgIbH98DPyF1tlu2h319fUoKSkxarN17w4Qwbs046LFOPjCJMxLCYGW1hW3mPrhSRy8XG33c/umViVe2XUJ939+BkX1rfAX8PDfx0Zh3Zx4RgtH9BZzcfNxccwmxuwN1g3kJdgdIncuPnpwBGYkBOKNX/NR2dKOVd9mY0yUL/46bRhGRjCbqrmvtCnV+O+ZUnx+/CYk7bp9/wvTQ7FmZix8PN1sbJ0OhUKBS5cuGbWlpKRYPW7eHLa3gGAXTI8PxKQhg/DZ8UJsOlGEs0VNWPDFGUwcIsaTU6IxNsrPZmvHgG6X287MCnx+/KYhTHh4oABvz0tA2mD7+lHKy8szipsHYJO4eXMQwRMMuLux8ddpw7BoVDg+PVaInRfKcepGA07daMDQAC88MjYCc0YEQ8jnWs2mm/VyfHu2FD9mVkDeuZoQ5uuOF+4aivtSQsBm2e5HyBzdxc17ezOfK683EMETTAjxdsc/5yXiycnR2HTyJn7KrsT1Wjle23sZb+zLx5ShgzA7ORiThg6CyN3y4i9vasPByzX45VIV8iolhvaoQZ5YMSEKC9JDrZKOqj+UlZWhrq7OqM0enHV6iOAJ3RLm64G37kvEy9OHY1dWBbafL0NhnRyH82txOL8WLApICvXG+Bg/JIV6Iy5IiFAf9z4N/dUaLW7WtyK/WoLMkmacLmxAaWOb4TibRWHK0EF4dFwEJsSIwbKzHr0r5uLmo6OjbWSNKUTwhB4RuXOxfEIkHhsfgYIaGX69VIWDV2pQVN+KnPIW5JT/WQRDwOMgxMcdwd7uCBDy4cVjg89lw43NgkKtRbtKA1mHCtWSDlS1tKOiuR2KLvn02SwKaYN9MDs5GLMSAh2mjJZcLrdZvvneQgRP6DUURSE2SIjYICFemTEc1ZJ2/FHYiLNFjcivkuJGnQwyhRoFNTIU1PQ+046nGxuxQUIkhIgwPkaMMVG+EFjRT2ApzMXNWyvffG8hgif0myCROx5IC8UDaaEAAKVai5LGVlS2tKO6pQO10g50qDToUGmg1GjB4+h6e083NgJEfIR4uyPE2x3hvh52P1TvCa1Wa5JvPiEhwaZx8+YggidYDDcOC0MDBBgaYJ2Uz/aEPcbNm8M+XZ0EgoPR1VkXEhKC4OBgG1nTPUTwBMIAqa+vR3FxsVGbPfbuABE8gTBgupZ+9vDwQHy89fLd9wUieAJhACiVSruNmzcHETyBMAByc3OhUBinAE9LS7ORNT1DBE8g9BOapk2G80OHDoWPj20Tb9wOIngCoZ+Ul5ejtrbWqM1edsV1BxE8gdBPui7F+fj4ICYmxkbW9A4ieAKhH8jlcuTn5xu12VvcvDmI4AmEfpCdnW33cfPmIIInEPpId3HzHh4eNrKo9xDBEwh95Pr165BKpUZt9hpZ1xUieAKhj3R11gUHB9tl3Lw5iOAJhD7Q0NCAoqIiozZH6d0BIngCoU90DbRxd3e327h5cxDBEwi9RKlUIicnx6gtJSUFXK7jZOchgicQekleXp5J3Ly9R9Z1hQieQOgF5uLmhwwZYtdx8+YggicQekFFRQVqamqM2hytdwdITjtG0Wg0UKlUtjaDYAEuXrwIT09Pw3uhUIjQ0FB0dHQw8jwulws22/JFMYngGYCmadTU1KClpaXnkwl2j1arhb+/PwYNGmRoc3d3R2lpKaPP9fb2RmBgoEXj84ngGUAvdn9/f3h4eNj9hgrC7WltbYVAYJyJ18/Pj5EeGNB1GG1tbYaSVUFBQRa7NxG8hdFoNAax+/n52docwgChaRoSicQoZZW7u7vR8J4J9Pns6+rq4O/vb7EfF+K0szD6ObsjbKQg9ExHRwc0Go1RG9Ni16P/DlnSD0QEzxBkGO8ctLW1Gb3ncrlwc3OzyrOZ+A4RwRMI3aBSqUwCbazVuzMFETyBYAaapk1WWSiKAp/Pt5FFloEInkAwg0QiMZk7e3p6gsVybMk4tvVOjEZL48zNRvycU4kzNxuh0dKMP7Ourg5PPPEEwsPDwePxEBgYiOnTp+PMmTMAdD3c3r17LfKskpISUBRlshnF1tA0DfAM4KUAAB3vSURBVKlUajJ353A48PLyspFVloPRZbm3334b+/fvR05ODtzc3HoViELTNNavX48vv/wSzc3NGD16ND799FOH2oI4UA5ersb6X/NRLfkziitIxMfrs+MwI8Fya7Jduf/++6FSqbBt2zZERUWhtrYWv/32G5qamiz6HKVSadH7WQKaptHR0QGpVGrilacoCj4+Pg7fuwMM9/BKpRILFizAk08+2etr3n33XXz44Yf45JNPcOHCBQQGBmLq1KkmpXidlYOXq/Hkt9lGYgeAGkkHnvw2GwcvVzPy3JaWFpw+fRrvvPMO7rjjDgwePBijRo3CmjVrcM899yAiIgIAMG/ePFAUZXh/8+ZNzJ07FwEBAfDy8sLIkSNx9OhRo3tHRETgrbfewrJlyyASibBy5UpERkYC0G0vpSgKU6ZMYeRz9QaVSoXGxkY0NzebiB3QRbw50hbY28Go4NevX48XX3wRiYmJvTqfpmls3LgRr776KubPn4+EhARs27YNbW1t+P77781eo1AoIJVKjV6OikZLY/2v+TA3eNe3rf81n5HhvZeXF7y8vLB3714TzzTwZ1qnLVu2oLq62vBeLpdj1qxZOHr0KC5evIjp06dj9uzZKCsrM7r+vffeQ0JCArKysrB27VqcP38eAHD06FFUV1fjp59+svhn6gl9kFR9fX23ow6BQGAIgnEG7GqMUlxcjJqaGkybNs3QxuPxMHnyZGRkZJi9ZsOGDRCJRIZXWFiYtcy1OOeLm0x69luhAVRLOnC+2LJDbEA3R926dSu2bdsGb29vjB8/Hn//+9+Rm5sLAIY4cn18t/59cnIynnjiCSQmJmLIkCF46623EBUVhV9++cXo/nfeeSdeeuklxMTEICYmxnC9n58fAgMD4evra/HP1B00TUMul6Ours5krq7Hzc0NYrHYJKTW0bErweu3HwYEBBi1BwQEmGxN1LNmzRpIJBLDq7y8nHE7maJO1rudV709r6/cf//9qKqqwi+//ILp06fj+PHjSE1NxdatW7u9prW1Fa+88gri4uLg7e0NLy8vFBQUmPTw9rKVtKOjA/X19ZBKpaBp05ESm82Gj48P/Pz8rBZgY036LPh169aBoqjbvromCugrXSOMaJruNuqIx+NBKBQavRwVf0Hv1nh7e15/4PP5mDp1Kv7xj38gIyMDy5Ytw+uvv97t+S+//DJ2796Nt99+G6dOnUJOTg4SExNNhsi2DlhRq9VobGxEU1MT1Gq1yXGKoiAQCDBo0CC4u7s7baRkn730zzzzDBYtWnTbc/QOnb4SGBgIQNfT37pDqK6uzqTXd0ZGRfoiSMRHjaTD7DyeAhAo4mNUpPWGv3FxcYalOC6Xa+LUOnXqFJYtW4Z58+YB0M3pS0pKeryvvvc05ySzJFqtFjKZDK2trd2e4+7uDoFAYLc13S1Jnz+hWCyGWCxmwhZERkYiMDAQR44cQUpKCgCdp//EiRN45513GHmmPcFmUXh9dhye/DYbFGAken1/8/rsOLBZlu99GhsbsWDBAjz22GNISkqCQCBAZmYm3n33XcydOxeA7of8t99+w/jx48Hj8QzFE3/66SfMnj0bFEVh7dq1RiWYusPf3x/u7u44ePAgQkNDwefzIRKJLPZ59FtMZTJZt/ZwuVwIhULweDyLPdfeYXQOX1ZWhpycHJSVlUGj0SAnJwc5OTmQy+WGc4YPH449e/YA0A2rXnjhBfzzn//Enj17cPnyZSxbtgweHh5YvHgxk6baDTMSgvD5klQEioyH7YEiPj5fksrYOryXlxdGjx6Njz76CJMmTUJCQgLWrl2LlStX4pNPPgEAfPDBBzhy5AjCwsIMP8gfffQRfHx8MG7cOMyePRvTp09Hampqj8/jcDj4+OOPsWnTJgQHBxt+VCyBQqFAQ0MDJBKJWbGzWCx4e3tDLBa7lNgBgKLNeS4sxLJly7Bt2zaT9mPHjhnWXSmKwpYtW7Bs2TIAfwbebNq0ySjwJiEhoVfPlEqlEIlEkEgkJvN5pVKJDRs2GLU9//zz8Pb27vuH64aOjg4UFxcjMjJyQHHXGi2N88VNqJN1wF+gG8Yz0bM7E2q1GlKp9LZpp/TLj44QRHO779Ltvue3g9FJy9atW2/r4QVg4imlKArr1q3DunXrmDPMAWCzKIyNJgk0eoNWq0VrayvkcrlZzzugc0YKhUKXmKffDtf+9ASH5nbhsHo4HA6EQqHD73KzFETwBIdEqVRCKpV2GyFHURSEQiHJKdgFIniCQ6HRaCCTybqNkAN0a/5eXl6MJZl0ZIjgCQ4BTdNobW2FTCbrdp7u5uYGkUjkNBtdmIAInmDX0DQNhUIBiUTS7TydzWYb5ulk+H57iOAJdotKpYJUKjW7ew/QzdP1y2xE6L2DCJ5gd/QmHNbDwwMCgYDM0/sIETzBbuhtOKxIJHLKnWzWgAieYBfo5+nmdrIBunBYoVDo1DvZrIH9xxe6Gh0SQFJp/pikUnecQTIyMsBmszFjxgxGn6NHrVajqakJjY2N3Yrdy8uL1OmzEETw9kSHBPj2fmDrLEBSYXxMUqFr//Z+RkX/9ddf49lnn8Xp06dNklhYEq1WC6lUirq6um5j3/l8Pvz9/SEUCh0i9t0RIP+K9oRCDrTWA80lwNZ7/hS9pEL3vrlEd1whv91d+k1rayt27tyJJ598Evfee6/RPojjx4+Doijs378fycnJ4PP5GD16NPLy8ozusXv3bsTHx4PH4yEiIgIffPCB0fHq6mrMmjULnp6eSEhIwJ49ezB69Ghs3rzZcE5bWxvWrl2L4cOHw9fXF3feeScuXbpkdJ9ff/0VaWlp4PP5iIqKwvr167sdIRD+hAjenhCFAMv2Az4Rf4q+7NyfYveJ0B0XhTDy+B07dmDYsGEYNmwYlixZgi1btpgEubz88st4//33ceHCBfj7+2POnDmGgg1ZWVlYuHAhFi1ahLy8PKxbtw5r1641+uFYunQpysvL8eOPP2Lz5s347rvv0NDQAODPefqyZctQX1+PAwcOICsrC6mpqbjrrrsM6bIPHTqEJUuW4LnnnkN+fj42bdqErVu34u2332bk38WpoJ0MiURCA6AlEonJMYVCQa9bt87o1dzcbNHnt7e30/n5+XR7e3v/b9JSTtMbk2j6deGfr41JunYGGTduHL1x40aapmlapVLRYrGYPnLkCE3TNH3s2DEaAP3DDz8Yzm9sbKTd3d3pHTt20DRN04sXL6anTp1qdM+XX36ZjouLo2mapi9dukQDoA8cOEBXVlbSlZWV9OnTp2kA/7+9Mw+K4ura+DP7gKAiog7KgBsgwQ1xw/3DEjfiEkFFETVKxC0GLaRcCO5L3MuUUQtRCHGJgK9bmSABjKIRFXdESzFogOCCiOA42/3+wOnQDCADMwwy91fVVdO3T3efXs509733PJesX7+eqFQqkpiYSBo3bkxkMhlrO+3btyd79+4lhBAyYMAAsn79etby6OhoIpFI9Hg2jE9V91JV93lV0Cd8faRJG2DcPnbZuH2l5QYiMzMTV69eZeTL+Hw+Jk6ciAMHDrDs+vbty/xu1qwZnJyckJGRAQDIyMhAv379WPb9+vXDo0ePoFAocOPGDfD5fJZsuaOjI6ysrGBmZgYul4vr16/j3bt3sLa2ZjrVWFhYICsrC48fPwZQ+iaxevVq1vLZs2cjNze3yj72FNosV2m/bKNS+ByID2SXxQd+fJ03TNBHRERAqVSidev/PhcIIRAIBCgoKKhyXU3NOalAbFRzft+/f6/Vti4UCtGsWTPWNVCr1ZBIJEhOTtbaj0aoRK1WY9WqVRg/fryWDU2DrRqTCngejwcul8u68f79919YWVkZ0atylK2gs3IofbLHB/73TW+AoFcqlYiKisLWrVtZYwIApdLVMTExjOLQlStXIJVKAQAFBQV4+PAhnJ2dAZQKXl68eJG1fmpqKhwdHaFQKNChQwcolUrcvXsX7u7uaNasGZ48ecIagszNzQ15eXng8/mViqG6ubkhMzMTHTp00NcpMBlM6pWex+MxN6uGBw8eGMmbCij8R7uCTtpbuyKvsnb6GnL69GkUFBTg66+/hqurK2uaMGECIiIiGNvVq1cjMTGR0Rts3rw5xo4dCwBYvHgxEhMTsWbNGjx8+BCHDh3C7t27ERwcDLlcjg4dOmDAgAEICQnB/fv3cevWLQQGBrI60wwdOhR9+/bF2LFj8dtvv+Hp06dITU3FihUrGPnzsLAwREVFITw8HPfu3UNGRgaOHj2KFStW6PW8NERMKuABME8jDQ8fPqyWymqdILIAGtmUqY3/+CRv0ua/oG9kU2qnRyIiIjB06NAKVWO/+uor3Lx5Ezdu3AAAbNy4Ed9++y169OiB3NxcnDx5kunm6ubmhmPHjuHIkSNwdXVFWFgYVq9eDT8/P+a1fefOnbCxscGwYcMwbtw4zJ49G5aWlsyrOIfDwdmzZzFw4EDMnDkTjo6OmDRpEp4+fcpIlXt5eeH06dNISEhAz5490adPH2zbtg329vZ6PS8NEYOKWBqDT4n7vXnzBjt37mSVTZs2jRncsLbUWsRSVljazl5R01vhP6XBLtafnHN1SU5OxpAhQ1BQUKCz6OebN29YlWmaYZwA4Pnz57Czs8P58+fh6empV58/dz47Ecv6SNOmTSGRSJCb+98orBkZGXoL+FojblJ5QBuo/d2QkI+6cxouXrwIQgh69eqF3NxchISEwMHBAQMHDjSil6aDyQU8UPpaXzbgr127hrZt26JTp05G9Kpy1Go1SkpKjNqTTDOWQGGhbt16CSGsTyalUol169YhKysLlpaW8PDwQExMDFWpqSNMMuBdXFyQlJTEzBNCcPz4cfj6+sLJycmInlVMYWEh3r9/b1Qf3Nzc8M8/pZWFtWnr9vT0xOTJk2kSjJEwuUo7oHS4LA8PD1aZWq3GsWPH8OjRIyN5VTFyudzowa5PqAyVcTHJgAdKm3969OjBKlOr1Th69CjTo6s+UHZYrs8dDodj9FFkTR2TfKUHSm++UaNGMWPeaVCpVDhy5Aj8/PyMXpGnUCi0UkdFItFnOXoKl8uFWCz+LH1vSJj02edwOPD29oZarcbt27eZcqVSicOHD2PKlClGbdstr+nG4XBgZWVFc8MpNcbk7xwul4sxY8ZoDVapUCgQExODZ8+eGcUvpVKpVTnWqFEjGuyUWkHvHpQG/bhx4+Di4sIqVygU+Pnnn5na6bqkoqc7/f6l1BYa8B/hcrkYP368VrOcXC5HdHQ0cnJy6swXlUql9XQ3Nzenksw6EB4ejm7duhnbDTg4OGDHjh3GdoOBBnwZeDweJkyYgI4dO7LKP3z4gOjoaOTl5dWJH8XFxVppu4Z+uk+fPh0cDgdz5szRWjZ37lxwOBxMnz7doD7UFA6HgxMnTrDKlixZgsTExFpv+/nz5xAKhVo5GNUlLS0NgYGBnzasI2jAl4PP58PX1xft27dnlctkMkRHRyM/P9+g+9eMdS5TyeCV4AWvBC9whJw6qd22s7PDkSNHWO3+MpkMhw8f1soyrO9YWFjA2tq61ts5ePAgfH19UVJSgkuXLum8vo2NDczNzWvth76gAV8BGrWX8s1yJSUliIqKwosXLwy2b2M83TW4ublBKpUiLi6OKYuLi4OdnR26d+/Osj137hz69++Ppk2bwtraGqNHj9bqv5Camopu3bpBLBbD3d0dJ06cAIfDYZpBNcKYiYmJcHd3h7m5OTw8PJCZmcnaTlWClZqc+XHjxoHD4TDzFb3SHzhwgBHYlEgkmD9/fpXngxCCyMhI+Pv7w8/Pj5UmDABRUVGwsLBgddZasGABHB0dmTqY8q/04eHhkEqlEIlEsLW1xcKFC6v0Qd/QgK8EgUCAyZMnazXLFRcXIyoqyiASzmq1Gi8LX0KmkkGm+q/9XQEFShQlKFEYXr5pxowZiIyMZOYPHDiAmTNnatkVFxcjODgYaWlpSExMZCo+Nf3mi4qK4O3tjc6dO+PGjRtYs2YNli5dWuE+ly9fjq1bt+LatWvg8/ms/X1KsDItLQ0AEBkZidzcXGa+PHv27MG8efMQGBiIO3fu4OTJk58U0EhKSkJJSQmGDh0Kf39/HDt2DEVFRczyadOmYeTIkZgyZQqUSiXOnTuHvXv3IiYmpsI/6ePHj2P79u3Yu3cvHj16hBMnTrDkvuqEGivs1VNqKu5XGR8+fCARERFa4perVq0if/75J1Gr1Sz72ohYvnv3jrgedK1yMhQBAQFkzJgx5MWLF0QkEpGsrCzy9OlTIhaLyYsXL8iYMWNIQEBApevn5+cTAOTOnTuEEEL27NlDrK2tWedh//79BABJT08nhPwnjHn+/HnG5syZMwQAs151BCsBkPj4eJbN999/T7p27crM29rakuXLl+t0Tvz8/MiiRYuY+a5du5L9+/ezbF6/fk3atGlDgoKCSMuWLcnatWtZy+3t7cn27dsJIYRs3bqVODo6ErlcXq39G0LE0qQ73lQHoVCIKVOmIDo6mtU8RwhBYmIiHj16xPpW5PP5aNmyJd6+fVvpAAuVoau9IWjevDlGjRqFQ4cOgRCCUaNGMbnrZXn8+DFWrlyJK1eu4OXLl8yTPTs7G66ursjMzESXLl1Yedy9evWqcJ9dunRhfkskEgBAfn4+pFIprl+/jrS0NJYEtUqlgkwmQ0lJSbW+j/Pz85GTk6NTvv2bN28QFxfHkuyaOnUqDhw4gFmzZjFlVlZWiIiIgJeXFzw8PBAaGlrpNn18fLBjxw60a9cOw4cPx8iRI+Ht7V2nvQ9pwFcDkUiEqVOn4tixY8jKymIty87OZr3eN2rUCFZWVpDJZDVKZ/3f//0PACBTyTAxZSIAINk3GWZ8s1ocgW7MnDmT+b798ccfK7Tx9vaGnZ0d9u/fD1tbW6jVari6ukIulwOoWtCyPGVTYzXraP5A9CFYaWam+7n75ZdfIJPJ0Lt3b6aMfEz1vX//PqvPxoULF8Dj8ZCTk4Pi4uJKBSns7OyQmZmJhIQEnD9/HnPnzsUPP/yAlJSUOksPpt/w1UQsFmPq1KkYNGiQYffDE0PME8NSbMmUmfHNYC4wh7mgbmp7hw8fDrlcDrlcDi8vL63lr169QkZGBlasWAFPT0906tRJS9nW2dkZt2/fZo3trtGk04WygpXlJ02vQ4FAAJVKVek2LC0t4eDgoFMzXUREBBYvXoybN28y061btzBkyBCWdHdqaio2b96MU6dOoXHjxliwYEGV2zUzM8OXX36JXbt2ITk5GZcvX9YavceQ0Ce8DnC5XAwePBj29vaIjY2tcvzy2mJpaflpIwPB4/EYrfmKOvtYWVnB2toa+/btg0QiQXZ2ttarrJ+fH5YvX47AwECEhoYiOzsbW7ZsAQCd0mPDwsIwevRo2NnZwcfHB1wuF7dv38adO3ewdu1aAGCCuV+/fhCJRBWqEIeHh2POnDlo0aIFRowYgaKiIly6dKnCANVo+MXExGi1v0+ePBnLly/Hhg0bIJPJ4O/vjwULFmDEiBGQSqVwd3fH6NGj4ePjo7XdgwcPQqVSoXfv3jA3N0d0dDTMzMzqNF+DBnwNaNu2LebOnYt79+7h1atXrGU8Hg8ikQhmZmY1GsOcw+FALBZDKBTiTkDd/fOXpyqdNC6XiyNHjmDhwoVwdXWFk5MTdu3ahcGDB7PWP3XqFIKCgtCtWzd07twZYWFh8PPz00nrTyNYuXr1amzevBkCgQDOzs6s7+itW7ciODgY+/fvR+vWrfH06VOt7QQEBEAmk2H79u1YsmQJmjdvjgkTJlS4z4iICLi4uFTY2Wbs2LEICgrCqVOncPr0aTRq1Ajr168HAHzxxRfYtGkT5syZAw8PD5bGP1Aqr7Zx40YEBwdDpVKhc+fOOHXqlF76C1QXkxOxNDS1FrFswMTExGDGjBkoLCys0Xe1qUFFLCmfFVFRUWjXrh1at26NW7duYenSpfD19aXBbkRowFMMRl5eHsLCwpCXlweJRAIfHx86wquRoQFPMRghISEICQkxthuUMtBmOQrFhKABbyDqzfBVlM8WQ9xD9JVezwiFQnC5XOTk5MDGxgZCoZDKMlN0ghACuVyOFy9egMvl1qh5tzJowOsZLpeLtm3bIjc3t05VcigND3Nzc0ilUr3qGNKANwBCoRBSqRRKpbLKLp8USmXweDzw+Xy9vx3SgDcQHA4HAoGAjplGqVfQSjsKxYQwaMCvW7cOHh4eMDc3r/aY4hoxxbJTnz59DOkmhWIyGDTg5XI5fHx8EBQUpNN6w4cPR25uLjOdPXvWQB5SKKaFQb/hV61aBaA0LVAXRCIRWrVqVS3bDx8+sHKuNeOXv337Vqd9UiifE5r7W9fct3pZaZecnIwWLVqgadOmGDRoENatW4cWLVpUaLthwwbmj6UsdnZ2hnaTQjE6RUVFaNKkSbXt6yQ99uDBg1i0aBHevHnzSdujR4/CwsIC9vb2yMrKwsqVK6FUKnH9+nWIRCIt+/JPeLVajdevX8Pa2lqrSePt27ews7PDs2fPjJI6ayzocTe84yaEoKioCLa2trq10+skeUlK1UABVDmlpaWx1omMjCRNmjTRdVeEEEJycnKIQCAgsbGxNVq/LPpWtP1coMdtWsddFTq/0s+fPx+TJk2q0kYzGIA+kEgksLe3Z4n9UyiUmqFzwDdv3rxC2WJD8erVKzx79oyRL6ZQKDWHFx4eHm6ojWdnZyMrKwtXr17FxYsXMXLkSOTl5cHCwoJJCHB2dkbr1q3RqVMnvHv3DsuWLYOlpSVUKhVu3ryJWbNmQaFQYPfu3RV+w+sKj8fD4MGD61QLvD5Aj9u0jrtSDPm9EBAQUOE3flJSEmMDgERGRhJCCCkpKSHDhg0jNjY2RCAQEKlUSgICAkh2drYh3aRQTIYGJ2JJoVAqh/alp1BMCBrwFIoJQQOeQjEhaMBTKCZEgwr4CxcuwNvbG7a2tuBwODhx4sQn10lJSUGPHj0gFovRrl07/PTTT3XgqX7R9biTk5O1UpA5HA4ePHhQRx7Xng0bNqBnz56wtLREixYtMHbsWGRmZn5yvdjYWLi4uEAkEsHFxQXx8fF14G39oUEFfHFxMbp27Yrdu3dXyz4rKwsjR47EgAEDkJ6ejmXLlmHhwoWIjY01sKf6Rdfj1pCZmclKQ+7YsaOBPNQ/KSkpmDdvHq5cuYKEhAQolUoMGzasygE+L1++jIkTJ8Lf3x+3bt2Cv78/fH198ddff9Wh50bG2O2ChgIAiY+Pr9ImJCSEODs7s8q++eYb0qdPH0O6ZlCqc9xJSUkEACkoKKgjrwxPfn4+AUBSUlIqtfH19SXDhw9nlXl5eZFJkyYZ2r16Q4N6wuvK5cuXMWzYMFaZl5cXrl27BoVCYSSv6o7u3btDIpHA09MTSUlJxnanVmh0EJo1a1apTWXXOzU11aC+1SdMOuDz8vLQsmVLVlnLli2hVCrx8uVLI3lleCQSCfbt24fY2FjExcXByckJnp6euHDhgrFdqxGEEAQHB6N///5wdXWt1K6y652Xl2doF+sNJt/BuHzOPPnY8bAhDx7h5OQEJycnZr5v37549uwZtmzZgoEDBxrRs5oxf/583L59GxcvXvykbUXXuyFf6/KY9BO+VatWWv/u+fn54PP5sLa2NpJXxqFPnz6fZQryggULcPLkSSQlJaFNmzZV2lZ2vcs/9RsyJh3wffv2RUJCAqvs999/h7u7u8npyaenp39WKciEEMyfPx9xcXH4448/0LZt20+uU9n19vDwMJSb9Q/j1hnql6KiIpKenk7S09MJALJt2zaSnp5O/v77b0IIIaGhocTf35+xf/LkCTE3NyffffcduX//PomIiCACgYAcP37cWIdQI3Q97u3bt5P4+Hjy8OFDcvfuXRIaGkoA6EVVqK4ICgoiTZo0IcnJySQ3N5eZSkpKGBt/f38SGhrKzF+6dInweDyyceNGkpGRQTZu3Ej4fD65cuWKMQ7BKDSogNc0N5WfAgICCCGl6bqDBg1irZOcnEy6d+9OhEIhcXBwIHv27Kl7x2uJrse9adMm0r59eyIWi4mVlRXp378/OXPmjHGcryEVHS/KpFoTQsigQYOYc6Dh119/JU5OTkQgEBBnZ+fP6k9OH9D0WArFhDDpb3gKxdSgAU+hmBA04CkUE4IGPIViQtCAp1BMCBrwFIoJQQOeQjEhaMBTKCYEDXgKxYSgAU+hmBA04CkUE+L/AUtkATzWZpAHAAAAAElFTkSuQmCC",
"text/plain": [
"Figure(PyObject <Figure size 640x480 with 1 Axes>)"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/plain": [
"PyObject <matplotlib.legend.Legend object at 0x1156b6438>"
]
},
"execution_count": 29,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"plot(path.r,path.z)\n",
"scatter(gcp.r,gcp.z,marker=\"o\",label=\"Start\")\n",
"scatter(stat.rm,stat.zm,marker=\"x\",label=\"Apogee\")\n",
"scatter(M.axis...,marker=\"+\",label=\"Magnetic Axis\")\n",
"plot(g.rlim,g.zlim,color=\"gray\",lw=3)\n",
"gca()[:set_aspect](\"equal\")\n",
"legend()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Julia 1.0.1",
"language": "julia",
"name": "julia-1.0"
},
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "1.0.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment