Skip to content

Instantly share code, notes, and snippets.

@atmyers
Created August 28, 2014 06:17
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save atmyers/f8616c9ed5a9d2b027e8 to your computer and use it in GitHub Desktop.
Save atmyers/f8616c9ed5a9d2b027e8 to your computer and use it in GitHub Desktop.
{
"metadata": {
"name": "",
"signature": "sha256:e2ba608e12eb3bfd50086ba32d29343e8a4bd949b7eea5f37fb777ca01bcdf89"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"import yt\n",
"ds = yt.load('ZeldovichPancake/plt32.2d.hdf5')\n",
"p = yt.ParticlePlot(ds.all_data(), 'particle_position_x', 'particle_velocity_x')\n",
"p.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stderr",
"text": [
"WARNING:yt:Setting code length to be 1.0 cm\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"WARNING:yt:Setting code mass to be 1.0 g\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"WARNING:yt:Setting code time to be 1.0 s\n"
]
},
{
"html": [
"<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAA+gAAAMgCAYAAACwGEg9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3XuUVXXd+PHPGUBAbiM3HUQoRbMHEwQZUfC6TCxT0gU8jHfyWiZgPiFaSIgmpKXLe3ZBsGSZ4K0kV14g8RE1AUtTUjSlZJCLTvrIRWHO7w9/jA4zg6jM2XvOfr3WOmtxvmef8TM1jufNd599cvl8Ph8AAABAokqSHgAAAAAQ6AAAAJAKAh0AAABSQKADAABACgh0AAAASAGBDgAAACkg0AEAACAFBDoAAACkgEAHAACAFBDoAAAAkAICHQAAAFJAoAMAAEAKCHQAAABIAYEOAAAAKSDQAQAAIAUEOgAAAKSAQAcAAIAUEOgAAACQAgIdAAAAUkCgAwAAQAoIdAAAAEgBgQ4AAAApINABAAAgBQQ6AAAApIBABwAAgBQQ6AAAAJACAh0AAABSQKADAABACgh0AAAASAGBDgAAACkg0AEAACAFBDoAAACkgEAHAACAFBDoAAAAkAICHQAAAFJAoAMAAEAKCHQAAABIAYEOAAAAKSDQAQAAIAUEOgAAAKSAQAcAAIAUEOgAAACQAgIdAAAAUkCgAwAAQAoIdAAAAEgBgQ4AAAApINABAAAgBQQ6AAAApIBABwAAgBQQ6AAAAJACAh0AAABSQKADAABACgh0AAAASAGBDgAAACkg0AEAACAFBDoAAACkgEAHAACAFBDoAAAAkAICHQAAAFJAoAMAAEAKCHQAAABIAYEOAAAAKSDQAQAAIAUEOgAAAKSAQAcAAIAUEOgAAACQAgIdAAAAUkCgAwAAQAoIdAAAAEgBgQ4AAAApINABAAAgBQR6wt57772YOHFiHH300dGxY8coKSmJ6dOnb9Nzb7vttigpKan3tnLlykaeHAAAgO2pedIDZN2qVati8uTJ0bNnz+jbt2/Mmzcvcrncp/oakydPji9+8Yu11jp06LA9xwQAAKCRCfSEdevWLVasWBFdu3aNhQsXxoABAz711/ja174W/fr1a4TpAAAAKBSnuCdshx12iK5du0ZERD6f/0xfI5/Px7vvvhubNm3anqMBAABQQAK9CBx++OHRoUOHaNOmTQwdOjSWLl2a9EgAAAB8Sk5xb8LatGkTo0aNisMPPzzat28fzzzzTPzsZz+Lgw46KBYtWhTdu3dPekQAAAC2kUBvwoYPHx7Dhw+vuX/cccfFkCFD4pBDDokrrrgibr755gSnAwAA4NMQ6EVm0KBBccABB8TDDz/c4DGVlZVRWVlZwKkAAIA0Kisri7KysqTH4P8T6EWoe/fu8dJLL9X7WGVlZey///6xfPnyAk8FAACkTbdu3eKZZ54R6Skh0IvQq6++Gl26dKn3scrKyli+fHn85je/iS9/+csFnoysGTt2bFx77bVJj0EG+FmjUPysUSh+1iiEF198MU4++eSorKwU6Ckh0JuIFStWRFVVVfTq1SuaN//w/7ZVq1bVCfE5c+bEokWLYsyYMVv9el/+8pd9djqNrrS01M8ZBeFnjULxs0ah+FmDbBLoKXDDDTdEVVVVzWnn999/fyxbtiwiIkaPHh3t27eP8ePHx4wZM+K1116LHj16RETEQQcdFP369Yv+/ftHhw4dYtGiRfHrX/86evToEZdcckli3w8AAACfnkBPgZ/+9Kfx+uuvR0RELpeLe+65J+6+++7I5XJx6qmnRvv27SOXy0Uul6v1vJEjR8YDDzwQf/rTn2Lt2rXRrVu3OOecc2LixIkNnuIOAABAOgn0FPjnP//5icdMmzYtpk2bVmtt8uTJMXny5MYaCwAAgAIqSXoAoHhVVFQkPQIZ4WeNQvGzRqH4WYNsEuhAo/HigkLxs0ah+FmjUPysQTYJdAAAAEgBgQ4AAAApINABAAAgBQQ6AAAApIBABwAAgBQQ6AAAAJACAh0AAABSQKADAABACgh0AAAASAGBDgAAACkg0AEAACAFBDoAAACkgEAHAACAFBDoAAAAkAICHQAAAFJAoAMAAEAKCHQAAABIAYEOAAAAKSDQAQAAIAUEOgAAAKSAQAcAAIAUEOgAAACQAgIdAAAAUkCgAwAAQAoIdAAAAEgBgQ4AAAApINABAAAgBQQ6AAAApIBABwAAgBQQ6AAAAJACAh0AAABSQKADAABACgh0AAAASAGBDgAAACkg0AEAACAFBDoAAACkgEAHAACAFBDoAAAAkAICHQAAAFJAoAMAAEAKCHQAAABIAYEOAAAAKSDQAQAAIAUEOgAAAKSAQAcAAIAUEOgAAACQAgIdAAAAUkCgAwAAQAoIdAAAAEgBgQ4AAAApINABAAAgBQQ6AAAApIBABwAAgBQQ6AAAAJACAh0AAABSQKADAABACgh0AAAASAGBDgAAACkg0AEAACAFBDoAAACkgEAHAACAFBDoAAAAkAICHQAAAFJAoAMAAEAKCHQAAABIAYEOAAAAKSDQAQAAIAUEOgAAAKSAQAcAAIAUEOgAAACQAgIdAAAAUkCgAwAAQAoIdAAAAEgBgQ4AAAApINABAAAgBQQ6AAAApIBABwAAgBRonvQAAPBxudxHf87nk5sDAKDQ7KADkBofj/P67gMAFDOBDkAqNBTj2zvSc7mPbgAAaSLQAcgMO/QAQJoJdAAyoVA79AAAn5VAByAVGrognAvFAQBZIdABSI0tY1ycAwBZ4mPWAEiVxoryfL7+09n9JQAAkBZ20AHIDDv0AECa2UEHYKs+vutcDEFbDN8DAFCc7KAD0CAfSwYAUDgCHYB6+VgyAIDCcoo7AHxGxXb6PwCQLDvoAPAZOP0fANjeBHrC3nvvvZg4cWIcffTR0bFjxygpKYnp06dv8/Orqqri7LPPji5dukTbtm3jiCOOiMWLFzfixEBWNLQjbKfY6f8AQOMQ6AlbtWpVTJ48Of7xj39E3759IyIit42v8Kqrq+OYY46JmTNnxujRo+MnP/lJrFy5Mg477LBYunRpY44NZISPJQMAKBzvQU9Yt27dYsWKFdG1a9dYuHBhDBgwYJufO2vWrFiwYEHMmjUrTjjhhIiIGDFiROy1114xceLE+O1vf9tYYwMZIsoBAArDDnrCdthhh+jatWtEROQ/5avgWbNmxS677FIT5xERnTt3jhEjRsR9990XH3zwwXadFYAPOf0fAGgMAr0JW7x4cfTr16/O+oABA2Lt2rXx0ksvJTAVQDY4/R8A2N4EehNWWVkZZWVlddY3ry1fvrzQIwFkSj7/0Q0A4PMS6E3Y+vXro2XLlnXWW7VqFRER69atK/RIAAAAfEYCvQlr3bp1bNiwoc76+vXrax4HAACgaXAV9yasrKys3tPYKysrI+LDK8Q3ZOzYsVFaWlprraKiIioqKrbvkAAAQOJmzpwZM2fOrLVWVVWV0DQ0RKA3YX379o358+dHPp+v9dnpTz31VLRp0yb22muvBp977bXX1nuBOQAAoPjUtxm3aNGi6N+/f0ITUR+nuDcRK1asiCVLlsTGjRtr1oYNGxZvvvlm3H333TVrq1evjrvuuiuOPfbYaNGiRRKjAgAA8BnYQU+BG264IaqqqmpOV7///vtj2bJlERExevToaN++fYwfPz5mzJgRr732WvTo0SMiPgz0gQMHxqhRo+KFF16ITp06xU033RT5fD4mTZqU2PcDAADApyfQU+CnP/1pvP766xERkcvl4p577om77747crlcnHrqqdG+ffvI5XK1TmOPiCgpKYk5c+bE97///bjuuuti3bp1UV5eHjNmzIg999wziW8FAACAzyiXz/v01izZ/D6ThQsXeg86pNzH/07Ob2oAYHvTBunjPegAKbTFCTN17gMAUHwEOkDKNBTjIh0AoLgJdAAAAEgBF4kDgAS4xgAAsCU76AAp01Csibji4RoDAEB9BDpACm0Z4+K8eLjGAADQEKe4A6SUKAcAyBY76AAAAJACAh0ACsg1BgCAhgh0ACgw1xgAAOrjPegAkABRDgBsyQ46AAAApIBABwAAgBQQ6AAAAJACAh0AAABSQKADAABACgh0AAAASAGBDgAAACkg0AEAACAFBDoAAACkgEAHAACAFBDoAAAAkAICHQAAAFJAoAMAAEAKCHQAAABIgeZJDwCQBrncR3/O55ObAwCA7LKDDmTex+O8vvsAAFAIAh3ItIZiXKQDAFBoAh0AAABSwHvQAaAJcJ0EACh+dtCBTGsodAQQaeI6CQCQDQIdyLwtY1yckyaukwAA2eEUd4AQ5QAAJM8OOgAAAKSAQAeAFHOdBADIDoEOACnnOgkAkA3egw4ATYAoB4DiZwcdAAAAUkCgAwAAQAoIdAAAAEgBgQ4AAAApINABAAAgBQQ6AAAApIBABwAAgBQQ6AAAAJACAh0AAABSQKADAABACgh0AAAASAGBDgAAACkg0AEAACAFBDoAAACkgEAHAACAFBDoAAAAkAICHQAAAFJAoAMAAEAKCHQAAABIAYEOAAAAKdA86QGAbMjlPvpzPp/cHAAAkFZ20IFG9/E4r+8+AAAg0IFG1lCMi3QAAKhNoAMAAEAKCHQAAABIAYEONKqGLgjnQnEAAFCbq7gDjS6fdxV3aExbXtPBv2MA0DQJdKAgBAM0jvouuJjL+XcOAJoip7gDQBO1tU9D8EkJAND0CHQAAABIAYEOAAAAKSDQAaCJ2tr7zL0HHQCaHoEOAE1YfSEuzgGgaXIVdwBo4gQ5ABQHO+gAAACQAgIdAAAAUkCgAwAAQAoIdAAAAEgBgQ4AAAApINABAAAgBQQ6AAAApIBABwAAgBQQ6AAAAJACAh0AAABSQKADAABACgh0AAAASAGBDgAAACkg0AEAACAFBDoAAACkgEAHAACAFBDoAAAAkAICHQAAAFJAoAMAAEAKCHQAAABIAYEOAAAAKSDQU2DDhg1x0UUXRbdu3WLHHXeMgQMHxsMPP/yJz7vtttuipKSk3tvKlSsLMDkAAADbS/OkByDi9NNPj9mzZ8cFF1wQe+65Z0ybNi2+/vWvx9y5c2PQoEGf+PzJkyfHF7/4xVprHTp0aKxxAQAAaAQCPWFPP/103HnnnXH11VfH9773vYiIOOWUU2KfffaJcePGxf/+7/9+4tf42te+Fv369WvsUQEAAGhETnFP2KxZs6J58+Zx9tln16y1bNkyzjjjjFiwYEG88cYbn/g18vl8vPvuu7Fp06bGHJWUyeU+ugEAAE2fQE/Y4sWLY6+99oq2bdvWWh8wYEBERDz77LOf+DUOP/zw6NChQ7Rp0yaGDh0aS5cubZRZSY8to1ykAwBA0+cU94RVVlZGWVlZnfXNa8uXL2/wuW3atIlRo0bF4YcfHu3bt49nnnkmfvazn8VBBx0UixYtiu7duzfa3CSnoRjP5SLy+cLOAgAAbD8CPWHr1q2Lli1b1llv1apVzeMNGT58eAwfPrzm/nHHHRdDhgyJQw45JK644oq4+eabt//AAAAANAqnuCesdevWsWHDhjrr69evr3n80xg0aFAccMAB2/QxbQAAAKSHHfSElZWV1Xsae2VlZUREdOvW7VN/ze7du8dLL7201WPGjh0bpaWltdYqKiqioqLiU//zKKx8vv7T3J3eDgBAQ2bOnBkzZ86stVZVVZXQNDREoCdsv/32i3nz5sW7774b7dq1q1l/6qmnIiKib9++n/prvvrqq9GlS5etHnPttdf6aLYmbMtIF+cAAGxNfZtxixYtiv79+yc0EfVxinvChg0bFps2bYpbb721Zm3Dhg0xbdq0GDhwYOy6664REbFixYpYsmRJbNy4sea4VatW1fl6c+bMiUWLFsXRRx/d+MOTqHz+oxtAfXwcIwA0LXbQE1ZeXh7Dhw+Piy++OFauXBl77LFHTJ8+PZYtWxbTpk2rOW78+PExY8aMeO2116JHjx4REXHQQQdFv379on///tGhQ4dYtGhR/PrXv44ePXrEJZdcktS3BEAK1PdxjP5CDwDSTaCnwIwZM2LChAlx++23x9tvvx19+vSJP/zhDzF48OCaY3K5XOS2eLU1cuTIeOCBB+JPf/pTrF27Nrp16xbnnHNOTJw48RNPcQegePk4RgBomnL5vP9UZ8nm95ksXLjQe9ABitTWTmn3X30ANtMG6eM96AAAAJACAh0AikxDu+R2zwEg3QQ6ABShLWNcnANA+mXiInEvvfRSzJkzJxYuXBhLly6NqqqqKCkpiS5dukRZWVkMGjQoDj/88Ojdu3fSowLAdiPKAaBpKepAnzdvXvz4xz+OiIgDDzwwRowYEZ06dYrS0tIoKSmJt99+O95666149tln46KLLopVq1bFBRdcECNHjkx4cgAAALKmKAN906ZNMW7cuGjdunVMnz49ysrKtnr8McccExERb731Vtxyyy1xyimnxE033RTt2rUrxLgAAABQnIF+1VVXxemnnx5f+cpXPtXzOnbsGJdcckm8/vrrMWnSpLj66qsbaUIAAACorSgDffz48Z/r+T179hTnAAAAFFSmruJ+1VVXxaRJk+KBBx6IdevWxbBhw6Jr164xZMiQePnll5MeDwAAgAzLVKD/61//iiFDhsTRRx8dP/zhD+PJJ5+MJUuWxOzZs+PGG29MejwAAAAyrChPcW/I3nvvHQMHDoyIiDvuuCPGjBkTHTt2jIiIbt26JTkaAAAAGZepHfQddtghIiJeeOGFePPNN+OQQw6peWzt2rVJjQUAAADZ2kGvrKyMa6+9NmbNmhUHHHBAHHTQQfHqq6/GDTfcEDvvvHPS4wEAAJBhmdpBv/jii6Njx44xfPjwePDBB2PNmjVx++23R4cOHaKkJFP/UwAAAJAymdpBb968eZx66qm11iZOnJjQNAAAAPCRot02fvDBB+O9995LegwAAADYJkW7g15VVRVnnHFGlJSUxNChQ+PYY4+NHXfcMemxAAAAoF5FG+gjR46MkSNHxjvvvBP3339/Tawfd9xxYh0AAIDUKdpA36x9+/Zx8sknx8knn1wr1ps1a1YT661bt056TAAAADKu6AP947aM9d///ve1ToP/xje+IdYBAABIRKYC/ePat28fJ510Upx00kk1sX7mmWfW7Kwfc8wxYh0AAICCyWygf9zHY/3dd9+N+++/P84888xo165d3HLLLUmPBwAAQAYI9C20a9euJtarq6uTHgcAAICMKNrPQf+07rjjjjprJSX+5wEAAKAwMreDvnz58vjnP/8ZmzZtqlmrrq6OqVOnxoknnpjgZAAAAGRZpgL9yiuvjB/84Af1PpbL5Qo8DQAAAHwkU+dw/+Uvf4lly5bFxo0bo7q6uua2adOmOPTQQ5MeDwAAgAzL1A56eXl5dO/evc56LpeL6667LoGJAAAA4EOZ2kEvKyuL559/vt7Hrr/++gJPAwAAAB/J1A76N7/5zTjuuONihx12iN122y3y+Xzkcrmorq6O+++/P37+858nPSIAAAAZlakd9BNPPDGaN28eZWVlkc/nIyIin8/X3AAAACApmdpB/8pXvhJTpkyp97GzzjqrwNOQJh+/iL+/qwEAAJKQqR30Dh06NPjYmDFjCjgJabLlJ+z5xD0AACAJmQr0IUOGxLXXXhtvvPFGncfGjh2bwEQkraEYF+kAAEChZeoU95EjR8Zbb70V3/ve9+o8llNkAAAAJChTgV5SUhK/+tWv6pzqns/nY/To0QlNBQAAABkL9EsvvTSGDh1a72Pjx48v8DSkQT5f/+nsLhQHFBsXwwSA9MvUe9D79u0bI0eOjPnz59esLVy4MCZPnhzHH398gpORpC1fqHrhChQbF8MEgKYhU4F+2223xapVq2LNmjU1a/3794/zzz8/fvKTnyQ4GUnL5z+6ARQTF8MEgKYjU6e4t2rVKh555JE666WlpbFx48YEJgIAAIAPZWoHffXq1Q0+9sorrxRwEgAAAKgtU4FeWloaY8aMiRdffDE2btwYH3zwQTz//PNx+umnR48ePZIeDwC2u4beuuMtPQCQPpkK9EmTJsXbb78dvXv3jpYtW0br1q1j3333jffeey8uv/zypMcDgEbhYpgA0DRk6j3oLVq0iBkzZsSFF14Yc+fOjerq6hg8eHCUl5cnPRoANCpRDgDpV5SB/vjjj8fgwYMbfLxPnz7Rp0+frX6NRx99NI444ojtPRoAAADUqyhPcc/n8zF+/Ph49913P/Vz161bF+PHj4+VK1c2wmQAAABQv6LcQT/44INj1113jXPOOSd22mmnOPXUU6N///7RvHn93251dXU8++yzMXv27Fi8eHFMnDgxDjjggAJPDQAAQJYVZaBHROy+++5xxx13xNNPPx033nhjPP7449GxY8fo2rVrdOjQISIiqqqqYs2aNbFixYooLy+Pb33rW3HFFVckPDkAAABZVLSBvll5eXnNReCWLl0a//73v2PVqlWRz+ejc+fOUVZWFl/60peipKQoz/YHAACgiSj6QP+4Xr16Ra9evZIeAwAAAOqwbQwAAAApINABAAAgBQQ6AAAApIBABwAAgBQQ6AAAAJACAh0AAABSQKADAABACgh0AAAASAGBDgAAACkg0AEAACAFBDoAAACkgEAHAACAFBDo/19VVVXSIwAAAJBhRR3oc+fO3eZjzzvvvEacBAAAALauqAP9iiuu2KbjKisr46GHHmrkaQAAAKBhRR3ojz32WDz99NNbPebnP/959O7dO9asWVOgqQAAAKCuog70iIgpU6bE1KlTY9OmTbXWn3vuuTjooIPi29/+dpSXl0fnzp0TmhAAAACKPNAvvfTSuPvuu2PQoEFx4oknxl/+8pdYu3ZtjBs3Lvr37x+vv/56/O53v4sHH3wwpkyZkvS4AAAAZFjzpAdoTD/84Q8jImLw4MExYMCAOP7442Px4sWxatWq+M53vhNXXHFFtGvXLiIiRo0aleSoAAAAZFxR76A/99xzERGxYsWKOOWUU+LBBx+MXXbZJYYOHRpHHXVUTZwDAABA0oo60EeOHBk33HBD7L333vHHP/4xrr766li4cGHMmjUr/vWvf8U555wTq1atioiI3//+9wlPCwAAQJYV9SnuL774YowePTqOO+64uP7662O33Xareezb3/52vPHGG3HRRRdFjx494oYbbojVq1cnOC0AAABZVtQ76M2aNYvZs2fHvffeWyvON9t1113j17/+dXTr1i3eeuutBCYEAACADxV1oH/1q1+N448//hOPO/vss2PAgAEFmAgAAADqV9SBfskll2zzsWeddVYjTgIAAABbV9SBPnjw4G0+9swzz2zESQAAAGDrijrQAQAAoKko6qu4U1i53Ed/zueTmwOAhvldDQDpZQed7eLjL/jquw9A8vyuBoB0E+h8bg29wPPCDyA9/K4GgPTLVKDfddddSY8AAAAA9cpUoH/3u9+NF198MekxAAAAoI7MXSTud7/7XXzwwQdxyCGHxFe/+tXIObfvc8vn6z9F0sWHANLD72oASL9MBfrcuXPjv/7rvyIi4rHHHotLL700dtppp6ioqIiysrKEp2vatnzh5wUfQPr4XQ0A6ZapU9w3x3lExMEHHxwHHXRQPPDAA9GjR48466yz4ve//32C0zV9+fxHNwDSye9qAEivTAX6t771rfj3v/8dkyZNii984QtRUVERe++9dzzzzDPxi1/8Irp06RITJkyIl19+OelRAQAAyJhMBfptt90WPXv2jIceeiguu+yyWL58edx4443Rp0+fiIgYOHBgTJ48OX7zm98kPCkAAABZk6n3oHfv3j3++Mc/Ru/evet9fOXKlTF8+PD40pe+VODJAAAAyLpMBfoll1zSYJxHRLRv3z4GDBgQJ5xwQgGnAgAAgIwFer9+/Wrdnz9/fsybNy86duwYp59+erRp0yauvvrqhKYDAAAgyzL1HvSbb7651v2DDz44JkyYEKeddlpMmDAhoakAAAAgY4Fen3w+H8uXL4+lS5cmPQoAAAAZVvSBfsEFF8T+++8fJSUlMX369CgpKal1a9asWfTu3TuOPPLIpEcFAAAgw4o+0K+55pp45pln4oorroghQ4bE3Llz49FHH625/fnPf45ly5bF6NGjE5lvw4YNcdFFF0W3bt1ixx13jIEDB8bDDz+8Tc+tqqqKs88+O7p06RJt27aNI444IhYvXtzIEwMAANAYMnORuPPPPz86d+4chx56aNKj1HL66afH7Nmz44ILLog999wzpk2bFl//+tdj7ty5MWjQoAafV11dHcccc0z87W9/i3HjxkWnTp3ipptuisMOOywWLlwYvXr1KuB3AQAAwOdV9Dvom7Vt2zbOOuusBh9/9dVXCzjNh55++um48847Y8qUKTF16tQ488wz49FHH42ePXvGuHHjtvrcWbNmxYIFC2L69OkxYcKE+M53vhPz5s2LZs2axcSJEwv0HQAAALC9FHWgv//++/H+++9v07HnnntuI09T16xZs6J58+Zx9tln16y1bNkyzjjjjFiwYEG88cYbW33uLrvsUusz2zt37hwjRoyI++67Lz744INGnR0AAIDtq6gDfd99943999+/5n737t3rXCRu8+2RRx4p+HyLFy+OvfbaK9q2bVtrfcCAARER8eyzz271uVt+rvvm565duzZeeuml7TssAAAAjaqo34N+/PHHxw477FBzv3v37nHvvfdG586dax2Xz+ejoqKi0ONFZWVllJWV1VnfvLZ8+fKtPvewww7b6nN79+69fQYFAACg0RV1oF955ZW17k+ZMqXWjvrHJfG+7XXr1kXLli3rrLdq1arm8YasX7/+Mz8XAACA9CnqU9y31LVr1xg5cmTMnz+/Zm3hwoUxefLkenejG1vr1q1jw4YNddbXr19f83hjPBcAAID0Keod9C3ddtttsWrVqlizZk3NWv/+/WOPPfaIqVOnxo9+9KOCzlNWVlbvaeyVlZUREdGtW7dGeW5ExNixY6O0tLTWWkVFRSKn+gMAAI1r5syZMXPmzFprVVVVCU1DQzIV6K1atar3YnClpaWxcePGgs9rXHZ+AAAgAElEQVSz3377xbx58+Ldd9+Ndu3a1aw/9dRTERHRt2/fBp/bt2/fmD9/fuTz+cjlcrWe26ZNm9hrr722+s++9tpr673IHAAAUHzq24xbtGhR9O/fP6GJqE+mTnFfvXp1g4+98sorBZzkQ8OGDYtNmzbFrbfeWrO2YcOGmDZtWgwcODB23XXXiIhYsWJFLFmypNZfIgwbNizefPPNuPvuu2vWVq9eHXfddVcce+yx0aJFi8J9IwAAAHxumdpBLy0tjTFjxsS5554be+65Z+Tz+fjHP/4RV199dfTo0aPg85SXl8fw4cPj4osvjpUrV8Yee+wR06dPj2XLlsW0adNqjhs/fnzMmDEjXnvttZo5hw0bFgMHDoxRo0bFCy+8EJ06dYqbbrop8vl8TJo0qeDfCwAAAJ9PpnbQJ02aFG+//Xb07t07WrZsGa1bt45999033nvvvbj88ssTmWnGjBkxduzYuP3222PMmDGxadOm+MMf/hCDBw+uOSaXy9U6jT0ioqSkJObMmRP//d//Hdddd12MGzcuunbtGo8++mjsueeehf42AAAA+Jxy+Xw+n/QQhfbXv/415s6dG9XV1TFo0KA44IADkh6pYDa/z2ThwoXegw4AABmmDdInU6e4b9anT5/o06dP0mMAAABAjUyd4h7x4YXUvv/970f//v2jb9++cd5558Ubb7yR9FgAAABkXKYC/dVXX43y8vJ4+eWXo7y8PPr27RvPP/987L///rF06dKkxwMAACDDMnWK+5QpU+Kxxx6L7t2711pfunRp/OQnP6n1cWcAAABQSJnaQe/Ro0edOI+I6NWrV7Rt2zaBiQAAAOBDmQr0tWvXNvhYdXV1AScBAACA2jIV6M2bN4+rr746/vOf/0RERD6fjzfeeCOmTJkSHTp0SHg6AAAAsixTgX7ppZfG0qVLo1OnTtGmTZto0aJF7LbbbvHss8/GhAkTkh4PAACADMvUReKaN28et9xyS3z3u9+Nhx9+ODZt2hSDBw+OAw44IP7617/6bHQAAAASk6lA32yfffaJffbZp9bahRdeGA8//HBCEwEAAJB1RRvoDz30UAwZMmSbj8/lco04DQAAAGxd0QZ6p06d4uijj44bbrghSkq2/lb7fD4fFRUVBZoMAAAA6iraQO/Vq1dMmjQpdt999206/vLLL2/kiQAgXT5+8lg+n9wcADTM7+psKdqruLdv3z4GDBiwzccfeeSRjTgNAKTLlu/s8k4vgPTxuzp7ijbQG/KrX/0qjjzyyDjvvPMiImLBggVx1VVXRXV1dcKTAUBhNPQCzws/gPTwuzqbMhXokydPjjlz5sSoUaOiU6dOERFx4IEHxogRI2Lq1KkJTwcAAECWZSrQq6qqYvbs2XHSSSfVBHpERM+ePeOdd95JcDIAAACyLlOB3qpVqwYfW7NmTQEnAYDkNHSRIRcfAkgPv6uzKVOB/vbbb8c111wT//nPfyIiorq6Op577rk477zzonPnzglPBwCFs+ULPC/4ANLH7+rsyVSgX3XVVfHwww/HTjvtFOPGjYsWLVpEnz594q233orLLrss6fEAoKDy+Y9uAKST39XZUrSfgx4RsXLlyujatWvN/TZt2sQDDzwQTz75ZDzxxBORy+Xi0EMPjX79+iU4JQAAABR5oI8YMSLmzZtXZ33gwIExcODAwg8EAAAADSjqU9wff/zxuPrqq2vecw4AAABpVdQ76F27do0999wzJk6cGOvWrYtvfvObMWTIkCgpKeq/lwAAAKAJKupAv/HGG2Po0KExdOjQWL9+fdx///1xzjnnRKdOneLkk0+OffbZJ+kRAQAAICKKPNCPP/74mj+3atUqRowYESNGjIg333wzZs6cGddcc0306dMnTjzxRB+zBgAAQKIyea73zjvvHGPHjo1zzz037rvvvujWrVsMHz486bEAAADIsKLeQa/PG2+8EbfffnvMmDEjlixZEq1bt47hw4fHGWeckfRoAAAAZFhR76BfcMEFERGxdu3a+M1vfhNHHXVU9OzZM37wgx9Ely5d4pe//GWsWLEifvvb38YRRxyR8LQAAABkWVHvoN91111RVVUVs2fPjv/7v/+LPfbYIy699NI49dRT4wtf+ELS4wEAAECNog705cuXx7333hsjR46M0047LQYNGpT0SAAAAFCvog70fffdN5588slo1apV0qMAAADAVhX1e9BPPvlkcQ4AAECTUNSB/j//8z9JjwAAAADbpKgDHQAAAJoKgQ4AAAApINABAAAgBQQ6AAAApIBABwAAgBTIXKD/6le/iiOPPDLOO++8iIhYsGBBXHXVVVFdXZ3wZAAAAGRZpgJ98uTJMWfOnBg1alR06tQpIiIOPPDAGDFiREydOjXh6QAAAMiyTAV6VVVVzJ49O0466aSaQI+I6NmzZ7zzzjsJTgYAAEDWZSrQW7Vq1eBja9asKeAkAAAAUFumAv3tt9+Oa665Jv7zn/9ERER1dXU899xzcd5550Xnzp0Tng4AAIAsy1SgX3XVVfHwww/HTjvtFOPGjYsWLVpEnz594q233orLLrss6fEAAADIsOZJD1BIbdq0iQceeCCefPLJeOKJJyKXy8Whhx4a/fr1S3o0AAAAMi5Tgb7ZwIEDY+DAgbXW/vrXv0afPn0SmggAAICsy9Qp7ltz4YUXJj0CAAAAGVa0O+gPPfRQDBkyZJuPz+VyjTgNAAAAbF3RBnqnTp3i6KOPjhtuuCFKSrZ+okA+n4+KiooCTQYAAAB1FW2g9+rVKyZNmhS77777Nh1/+eWXN/JEAAAA0LCifQ96+/btY8CAAbXWXnjhhRg5cmTMnz+/Zm3hwoUxefLkGDRoUKFHBAAAgBpFG+j1ue2222LVqlWxZs2amrX+/fvH+eefH1OnTk1wMgAAALKuaE9xr0+rVq3ikUceqbNeWloaGzduTGAiAEjex6+Tms8nNwcAZF2mdtBXr17d4GOvvPJKAScBgHTY8kNMfKgJACQnU4FeWloaY8aMiRdffDE2btwYH3zwQTz//PNx+umnR48ePZIeDwAKqqEYF+kAkIxMBfqkSZPi7bffjt69e0fLli2jdevWse+++8Z7773nKu4AAAAkKlPvQW/RokXMmDEjLrzwwpg7d25UV1fH4MGDo7y8POnRAAAAyLhMBfpmffr0iT59+tRae/XVV7f5M9MBoBjk8/Wfzu5CcUCxcTFMmoqiPsX9/fffj/fff3+bjj333HMbeRoASJ8tX6h64QoUGxfDpCkp6h30fffdN3bYYYf429/+FhER3bt3j+XLl9d7bM6/qQBklCgHitXWLobpdx9pVNSBfvzxx8cOO+xQc7979+5x7733RufOnWsdl8/no6KiotDjAQAAQI2iDvQrr7yy1v0pU6bE/vvvX++xEydOLMRIAAAAUK+iDvQt7bjjjrXuz58/P+bNmxcdO3aM0047LaGpAACAxuBimDQ1RX2RuC3dfPPNte4ffPDBMWHChDjttNPsoAMAQBFyMUyakkztoNcnn8/H8uXL4+WXX056FAAAoBGIcpqKot9Bv+CCC2L//fePkpKSmD59epSUlNS6NWvWLHr37h1HHnlk0qMCAACQYUW/g37NNddExIcXjHvsscdi/Pjxkf/YX6E1a9YsevXqFWVlZUmNCAAAAMUf6Judf/750blz5zj00EOTHgUAAADqKPpT3Ddr27ZtnHXWWVFdXZ30KAAAAFBHZgI9IuKoo46Kfv36JT0GAAAA1JGpQO/evXv8/Oc/r/ex22+/vcDTAAAAwEcyFehDhgyJl156Kd577706j02bNi2BiQAAAOBDmblIXETE9ddfH8uWLYvTTjstunTpEjvuuGNERFRXV8e///3vhKcDAAAgyzIV6KtXr47rr78+OnToUGs9n8/H6NGjE5oKAAAAMhboEyZMiKFDh9ZaW716dVRXV8fFF1+c0FQAAACQsUA/6aST6qy98MILsWTJkqiqqkpgIgAAAPhQpgI9IuLOO++MRx55JDZu3Bj5fD4iIjZs2BCLFy+OcePGJTwdAAAAWZWpq7jfeOONceGFF8a7774bf//73yMi4oMPPojnnnsupk+fnvB0AAAAZFmmdtCfeOKJWLJkSbRt2zauvPLKmvedV1VVxc033xzl5eUJTwgAAEBWZWoH/Stf+Uq0bdu25v6rr74aERGlpaXRrFmzpMYCAACAbAX6+vXrY8yYMXHHHXfEcccdFyeffHL885//jMrKynjiiSeSHg8AAIAMy1SgH3fccfHUU0/FihUronfv3jFw4MDYY489Ytddd3V6OwAAAInK1HvQ+/XrF08++WTN/Z/97GdxwgknxDvvvBNf//rXE5wMAACArMtEoD/77LPxyiuvRN++fWOPPfao9djgwYMTmgoAAAA+UtSnuG/atClOOeWU6NevXwwfPjz23nvv+PGPf5z0WAAAAFBHUe+g33jjjfH000/HZZddFjvvvHMsX748br311jj88MPjwAMPTHo8AAAAqFHUgX7PPffEX/7yl2jfvn3N2re+9a2YMmWKQAcAACBVivoU906dOtWK84iI3XbbLVq2bJnQRAAAAFC/og70du3a1bteX6Dfe++9jT0OAAAANKioA33Tpk3bfOwvf/nLRpwEAAAAtq6o34N+7733xrHHHhu5XK5mLZ/Px9///vd46qmnataqq6tjwYIFSYwIAAAAEVHkgb527dp45513agV6RESPHj2iurq65n4+n498Pl/o8QAAAKBGUQf6UUcdFXPmzNmmY7/xjW808jT1q6qqinHjxsU999wT69ati/Ly8vjpT38a++233yc+90c/+lFcdtllddZbtmwZ69ata4xxAShyH/87bX93DQCFVdSBft55523zsd/+9rcbcZL6VVdXxzHHHBN/+9vfYty4cdGpU6e46aab4rDDDouFCxdGr169tunr3HLLLdG2bdua+82aNWuskQEoYluccBa5nEgHgEIq6kA/5phjGuXY7WXWrFmxYMGCmDVrVpxwwgkRETFixIjYa6+9YuLEifHb3/52m77OsGHDomPHjo05KgBFbss4//i6SAeAwijqq7in3axZs2KXXXapifOIiM6dO8eIESPivvvuiw8++GCbvk51dXW888473kcPAADQhAn0BC1evDj69etXZ33AgAGxdu3aeOmll7bp6+y+++5RWloa7du3j1NOOSVWrly5vUcFAACgkRX1Ke5pV1lZGYcddlid9bKysoiIWL58efTu3bvB53fs2DHOP//8OPDAA6Nly5bx2GOPxY033hhPP/10PPPMM9GuXbvGGh2AIpPP13+au5OzAKBwBPp2ks/nY8OGDdt0bKtWrSIiYv369dGyZcsGH/+kK7GPHj261v3jjz8+ysvL46STToqbbropLrroom2aBwAi6ka6OAeAwnKK+3by5z//OXbcccdtum0+db1169b1Rv369etrHv+0KioqYpdddolHHnnk831DAGRSPv/RDaA+udxHN2D7soO+nXz5y1+O2267bZuO3WWXXSLiw1PZly9fXufxysrKiIjo1q3bZ5qle/fu8dZbb231mLFjx0ZpaWmttYqKiqioqPhM/0wAAIqfj2NsumbOnBkzZ86stVZVVZXQNDREoG8nO++8c5x66qmf6jl9+/aN+fPnRz6fj9zHfts99dRT0aZNm9hrr70+9Rz5fD5ee+216N+//1aPu/baa+u9QB0AANTHxzE2bfVtxi1atOgTu4HCcop7goYNGxZvvvlm3H333TVrq1evjrvuuiuOPfbYaNGiRc36smXLYsmSJbWev2rVqjpf8+abb47Vq1fH0Ucf3XiDAwAAsN3ZQU/QsGHDYuDAgTFq1Kh44YUXolOnTnHTTTdFPp+PSZMm1Tr21FNPjcceeyyqq6tr1nr27BkjR46MffbZJ1q1ahWPP/543HnnnbHffvvFOeecU+hvBwAAgM9BoCeopKQk5syZE9///vfjuuuui3Xr1kV5eXnMmDEj9txzz1rH5nK5WqfBR0ScfPLJ8cQTT8Ts2bNj/fr18YUvfCEuuuii+MEPflBzJXgAANgefBwjND6BnrDS0tL4xS9+Eb/4xS+2etzcuXPrrN16662NNRYAANTh4xihcQl0AABgm4lyaDwuEgcAAAApINABAAAgBQQ6AAAApIBABwAAgBQQ6AAAAJACAh0AAABSQKADAABACgh0AAAASAGBDgAAACkg0AEAACAFBDoAAACkgEAHAACAFBDoAAAAkAICHQAAAFJAoAMAAEAKCHQAAABIAYEOAAAAKSDQAQAAIAUEOgAAAKSAQAcAAIAUEOgAAACQAgIdAAAAUkCgAwAAQAoIdAAAAEiB5kkPAACkWy5X+34+n8wcAFDs7KADAA3aMs4bWgMAPj+BDgDUa2shLtIBYPsT6AAAAJACAh0AAABSQKADAPXa2sXgXCgOALY/gQ4ANKi+EBfnsH3lch/dgGzzMWsAwFYJcmg8W0Z5LuffOcgyO+gAAJCAhnbM7aRDdgl0AAAASAGBDgAAACkg0AEAIAENvdfce9AhuwQ6AAAkZMsYF+eQba7iDgAACRLlwGZ20AEAACAFBDoAAACkgEAHAACAFBDoAAAAkAICHQAAAFJAoAMAAEAKCHQAAABIAYEOAAAAKSDQAQAAIAUEOgAAAKSAQAcAAIAUEOgAAACQAgIdAAAAUkCgAwAAQAoIdAAAAEgBgQ4AAAApINABAAAgBQQ6AAAApIBABwD4f+3deVjU5f7/8deAu4AILoA7uea+4FKWUqdOanrMjeikLWYdM/iqbZaWmnos2zxphWa55JIevLRz2Xo0FBEUDE3yRGCuCJaK5AoGfH5/+GMSAQUZZj7MPB/XNVfMPfd87vc9M4GvuT8LAAAmQEAHAAAAAMAEqji6AAAAULlZLH/+bBiOqwMAgMqOFXQAAHDTrg7nxd0HAAClR0AHAAA3paQwTkgHAODmENABAAAAADABjkEHAACAU+M8CQAqC1bQAQDATSkp6BCAYCacJwFAZUJABwAAN+3aME44h5lwngQAlQ27uAMAgHIhlAMAYBusoAMAAAAAYAIEdAAAADglzpMAoLIhoAMAAMBpcZ4EAJUJx6ADAADAqRHKAVQWrKADAAAAAGACBHQAAAAAAEyAgA4AAAAAgAkQ0AEAAAAAMAECOgAAAAAAJkBABwAAAADABAjoAAAAAACYAAEdAAAAAAATIKADAAAAAGACBHQAAAAAAEyAgA4AAAAAgAkQ0AEAAAAAMAECOgAAAAAAJkBABwAAAADABKo4ugAAAOBaLJY/fzYMx9UBAIDZsIIOAADs5upwXtx9AABcGQEdAADYRUlhnJAOAMAVBHQAAAAAAEyAY9ABAABgM5xjAABuHivoAADALkoKa4Q458E5BgCgfAjoAADAbq4N44Rz58E5BgCg/AjoDnTixAlNmTJFwcHB8vT0lJubm7Zt21ambRw/flyjRo1S3bp1VadOHQ0dOlSHDh2qoIoBACg/w/jzBgAA/kRAd6Dk5GTNmzdPGRkZ6tSpkyTJUoavmc+fP6/g4GBt375dU6dO1cyZM7Vnzx7169dPmZmZFVU2AAAAAKACcJI4B+rRo4cyMzPl7e2tyMhIxcXFlen5H3zwgQ4cOKCEhAR1795dkjRgwAB16NBBb7/9tubMmVMRZQMAABRhGMXvzs6eEgBQeqygO5CHh4e8vb1v+vmRkZHq2bOnNZxLUps2bXT33Xdr3bp1tigRAACg1DjHAACUDwG9ksrPz9e+ffvUo0ePIo8FBQXpl19+0YULFxxQGQAAcGWcYwAAbh4BvZLKzMzU5cuX5e/vX+Sxgrb09HR7lwUAAAAAuEkcg24jhmEoJyenVH1r1KhR7vEuXbokSapevXqJ2y/oAwAAAAAwP1bQbWTbtm2qVatWqW4pKSnlHq9mzZqSVOyXAtnZ2YX6AAAAAADMjxV0G2nXrp2WLVtWqr5+fn7lHs/Hx0fVq1dXRkZGkccK2gICAkp8/sSJE4ucoC40NFShoaHlrg0AAACAuaxZs0Zr1qwp1JaVleWgalASArqNNGzYUGPGjLHbeG5uburYsaMSEhKKPLZr1y7dcsstql27donPnz9/vrp161aRJQIAAAAwieIW4xITEwtdEQqOxy7ulcTRo0eVnJxcqG3EiBFKSEjQ999/b237+eefFRUVpZEjR9q7RAAAAABAObCC7mCzZ8+WJO3fv1+StGLFCkVHR0uSpk2bZu03ZswYRUdHKz8/39r29NNP66OPPtKgQYP03HPPqUqVKnrnnXfk5+enZ5991o6zAAAAAACUFwHdwV599VVZLBYZhiGLxaJPPvlEkmSxWAoFdIvFIovFUui5Hh4e2rp1qyZNmqTZs2crPz9fwcHBevfdd+Xr62vXeQAAYA9X/ynkOtsAAGdDQHewq1fErycqKqrY9kaNGmndunW2LAkAAFO65ntqWSyEdACAc+EYdAAAYHrXhvMbtQMAUBmxgg4AAOAE2P0fACo/VtABAAAqueJ2/wcAVD4EdAAAYHolrQizUszu/wDgTAjoAACgUrg2jBPOAQDOhmPQAQBApUEoBwA4M1bQAQAAKjF2/wcA50FABwAAqOTY/R8AnAO7uAMAAJfkbJclc4Y5AICrYwUdAAC4HC5LBgAwIwI6AABwKfa6LJnF8ucNAIDSIKADAADYGCv0AICbQUAHAACwIXut0AMAnA8BHQAAuBQuSwYAMCsCOgAAcDlclgwAYEZcZg0AALikigrlhlH87ux8CQAAuBFW0AEAAGyMFXoAwM1gBR0AAKACEMoBAGXFCjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAAB3YFOnDihKVOmKDg4WJ6ennJzc9O2bdtK/fwZM2bIzc2tyK1mzZoVWDUAAAAAoCJUcXQBriw5OVnz5s1T69at1alTJ8XFxclisZR5OxEREfLw8LDed3d3t2WZAAAAAAA7YAXdgXr06KHMzEwlJydr0qRJN72dESNG6KGHHrLeQkJCbFglcPPWrFnj6BLgIviswV74rMFe+KwBromA7kAeHh7y9vYu93by8/N19uxZGYZhg6oA2+EfF7AXPmuwFz5rsBc+a4BrIqA7gcDAQHl7e8vLy0ujR4/Wb7/95uiSAAAAAABlxDHolZiPj4/CwsLUp08fVa9eXdHR0Xr//fcVHx+v3bt3y9PT09ElAgAAAABKiYBuI4ZhKCcnp1R9a9SoYZMxw8PDC91/4IEH1LNnT/3973/XBx98oBdffNEm4wAAAAAAKh4B3Ua2bdumu+66q1R9k5OT1bp16wqpIzQ0VM8++6y2bNly3YD+008/Vcj4wNWysrKUmJjo6DLgAviswV74rMFe+KzBHsgE5kNAt5F27dpp2bJlperr5+dXobU0btxYmZmZxT7m7++vgIAAPfzwwxVaA1Cge/fuji4BLoLPGuyFzxrshc8a7CEgIED+/v6OLgP/HwHdRho2bKgxY8Y4ugwZhqHDhw+X+Avd399fu3fvVkZGhp0rAwAAAGA2/v7+BHQTIaBXEkePHtXFixfVtm1ba9vJkydVv379Qv0+/PBDnTp1Svfdd1+J2+J/QgAAAAAwH4vBxbMdavbs2ZKk/fv3a+3atXr88cfVvHlzSdK0adOs/fr376/o6Gjl5+db22rVqqUHH3xQHTp0UI0aNRQTE6O1a9eqS5cu2rFjh81ORgcAAAAAqHgEdAdzc3OTxWKRYRjW/0qSxWJRXl6etV9wcLCio6MLtT355JOKjY3VsWPHlJ2drebNm2v48OGaOnWqateubfe5AAAAAABuHgEdAAAAAAATcHN0AQAAAAAAgIDuMnJycvTiiy8qICBAtWrVUu/evbV582ZHlwUnk5CQoGeeeUbt27eXh4eHmjVrppCQEKWmpjq6NLiAOXPmyM3NTR07dnR0KXBCiYmJGjJkiHx9fVW7dm117NhRCxYscHRZcDK7d+/W3/72NwUEBKh27dpq166dZs2apUuXLjm6NFRSFy5c0PTp03XffffJx8dHbm5uWr58ebF9f/rpJ913333y9PSUr6+vxowZo1OnTtm5YrCLu4sIDQ3V+vXrNWnSJLVq1UpLly5VQkKComsGN9QAABRJSURBVKKidPvttzu6PDiJESNGKC4uTiNHjlSnTp2UkZGhhQsX6vz589q5c6fat2/v6BLhpNLS0tSmTRu5ubmpRYsW2rdvn6NLghP59ttvNXjwYHXv3l0hISHy8PDQgQMHZBiGXn/9dUeXByeRlJSkoKAgBQQE6B//+Id8fHwUGxurZcuWaciQIdq4caOjS0QldPjwYQUGBqpZs2Zq0aKFtm7dqmXLlhW5PHRaWpq6du2qunXrKjw8XOfOndNbb72lpk2bKj4+XlWrVnXQDFwPAd0FxMfHq3fv3nrrrbc0efJkSVdW1Dt06KAGDRpox44dDq4QziIuLk5BQUGqUuXPKzgeOHBAHTt21IgRI/Tpp586sDo4swcffFCnT59Wbm6uTp06paSkJEeXBCdx9uxZtW7dWn379lVkZKSjy4ETmzp1qubOnav9+/erXbt21vZHH31UK1as0JkzZ1SnTh0HVojK6PLly8rKylKDBg30/fffKygoqNiA/vTTT2vFihVKTk5W48aNJUlbtmzRPffco0WLFmncuHGOKN8lsYu7C4iMjFSVKlX05JNPWtuqV6+usWPHKi4uTsePH3dgdXAmffr0KRTOJally5a69dZblZyc7KCq4Oyio6O1fv16zZ8/33pFDMBWVq9erd9++01z5syRdGV30asveQrYSs2aNSVJDRo0KNTu5+cnd3d3VatWzRFloZKrVq2a9TN1vXXZ9evX6/7777eGc0m6++671bp1a61bt67C68SfCOguYM+ePWrdurU8PDwKtQcFBUmS9u7d64iy4CIMw9Cvv/6qevXqOboUOKG8vDyFhYVp3LhxHEKBCrF582Z5eXnp2LFjatOmjTw9PVWnTh09/fTTysnJcXR5cCKPP/64GjZsqLFjx+qHH37QsWPHtHbtWkVERCg8PNwa4AFbO378uE6ePKkePXoUeSwoKEh79uxxQFWuq8qNu6Cyy8jIkL+/f5H2grb09HR7lwQXsmrVKqWnp2v27NmOLgVOKCIiQkePHtV3333n6FLgpFJTU5Wbm6uhQ4fqiSee0BtvvKGoqCgtWLBAWVlZWr16taNLhJMICAjQjh07NHDgQHXt2tXaPm3aNL322msOrAzOLiMjQ5JKzAuZmZn6448/OA7dTgjoLuDSpUuqXr16kfYaNWpYHwcqQnJysiZMmKDbbrtNjzzyiKPLgZM5ffq0Xn31Vb366qvy9fV1dDlwUufPn9fFixc1fvx4zZ8/X5I0dOhQXb58WYsWLdJrr72mli1bOrhKOINff/1VAwYMkCR99NFH8vX11aZNmzRnzhw1bNhQEyZMcHCFcFYFWeBGeYGAbh8EdBdQs2bNYnfDy87Otj4O2NqJEyc0aNAg1a1bV5GRkRwXDJubNm2a6tWrp7CwMEeXAidW8DcyNDS0UHtoaKgWLVqknTt3EtBhE7NmzdLx48eVkpKigIAASVe+DMrPz9eLL76o0NBQ+fj4OLhKOKOC33PkBXPgGHQX4O/vX+xu7AW7sxT8EQBs5ffff9eAAQN09uxZff311/Lz83N0SXAyqamp+uijjxQWFqa0tDQdPnxYhw8fVnZ2ti5fvqwjR47ozJkzji4TTqDgb2TDhg0LtRecdInPGWwlJiZGXbt2LfLvssGDB+vixYucMwgVpmDX9oJscLWMjAz5+vqyem5HBHQX0LVrV6WkpOjcuXOF2nft2iVJ6tKliyPKgpPKzs7W4MGDdeDAAW3atElt27Z1dElwQsePH1d+fr7Cw8MVGBhovcXHxyslJUUtWrTQrFmzHF0mnEDBSZPS0tIKtRd88V2/fn271wTn9McffygvL6/YdknKzc21d0lwEY0aNVL9+vWVkJBQ5LH4+Hiygp0R0F3AiBEjlJeXp8WLF1vbcnJytHTpUvXu3VuNGjVyYHVwJnl5eQoJCdGuXbv073//W7169XJ0SXBSHTt21IYNG7Rx40brbcOGDWrfvr2aNWumjRs3auzYsY4uE05g1KhRkqSPP/64UPuSJUtUtWpV9e/f3wFVwRl169ZNiYmJSk1NLdS+Zs0aubu7q1OnTg6qDK5g+PDh2rRpU6EvI7ds2aLU1FSNHDnSgZW5HotxvQviwWmEhIRow4YNmjRpkm655RYtX75cu3fv1pYtW9S3b19HlwcnMXHiRL333nsaPHhwsb/MH374YQdUBVfSv39/nT59WklJSY4uBU7kiSee0CeffKJRo0bpzjvv1NatWxUZGamXX36ZK1TAZvbt26fevXvLy8tLzzzzjHx8fLRp0yZ9/fXXGjdunBYtWuToElFJLVy4UFlZWUpPT1dERISGDRtmXRUPDw+Xl5eX0tLS1LVrV3l7e+v//u//dO7cOb355ptq2rSpEhIS2MXdjgjoLiInJ0evvPKKVq5cqTNnzqhz586aNWuW7rnnHkeXBicSHBys6OhoFfdrxWKxFLvrHmBLwcHBOn36tPbt2+foUuBEcnNz9c9//lNLly5Venq6mjdvrgkTJig8PNzRpcHJxMfHa/r06YqNjVV2drYCAwP1yCOP6IUXXpCbGzu+4ua0aNFCR44ckSTrSXsNw5DFYtGhQ4fUtGlTSdL//vc/TZ48WTExMapevboGDRqkt99+m0N57IyADgAAAACACfBVHAAAAAAAJkBABwAAAADABAjoAAAAAACYAAEdAAAAAAATIKADAAAAAGACBHQAAAAAAEyAgA4AAAAAgAkQ0AEAAAAAMAECOgAAAAAAJkBABwAAAADABAjoAAAAAACYAAEdAAC4lO3btysmJqZCx3j99deVl5dXoWMAAJwPAR0AgJtw9uxZJSYm6r///a+jS0EZHDx4UF988YX69u1boeMMGjRIL7zwQoWOAQBwPgR0AIApfPHFFxo4cKCqVKkiLy8vjRkzRo899phGjBih4OBgvffeexW+IvnSSy+pbdu2unz58g37/vzzz5o2bZr++te/lmtMe867LPMrS9/yiomJ0cCBA1W1alW5ubmpffv2SkpKkiTNnj1b1apVU7Vq1fSXv/xFUVFR5Rpr2rRpmjp1qi3Kvq6OHTuqZs2a2rp1a4WPBQBwHlUcXQAAANKVFcdBgwapX79+8vLy0ooVK6yPnT17VkFBQUpISNCnn35qk/EmTZqkd999t1BbQRAsjaCgIC1cuFAtW7YsVx32nHdJ8yvva1Feffv21ZdffqkNGzZo+PDhGj9+vDp27CjpSqBOTEzUlClT1LNnz3KNEx8fr7p168rT09MWZd/QhAkT9Mgjj6h///52GQ8AUPmxgg4AMBWLxSKLxVKozcvLS6NHj9aqVau0d+9em4yzffv2Im0zZ87Uvn37Sh1M3dxs92fUHvMuaX62eC1s4YEHHtCQIUM0bdo0nThxQtKV1fVevXqVO5xL0vvvv68xY8aUezul5e/vL0lKSUmx25gAgMqNgA4AqBQaNmwoSfrll1/KtZ3z589r7ty5SkxMtEVZFc5W8y6OGV+LBQsWKC8vT+Hh4crJydGCBQv03HPP2WTbMTEx6t69u022VVp9+vTRl19+adcxAQCVF7u4AwAqhfT0dLm7u6tLly6SpMTERP3nP/9RjRo1lJSUpNtuu00TJkyQJK1atUqzZs1SvXr1FBYWpm+++UZpaWn69ttvtWbNGkVHR0uSHnvsMUnS+PHjdfHiRU2ZMkVpaWlKS0uzjvvHH39o+vTpysrKUr169ZSdna1HH31Ut956a7F1XrhwQc8//7zq1q2r7OxseXt765VXXrHZvHNzc/Xaa68pNzdXtWrV0tGjRzVlyhQFBgYqNzdXc+fOlb+/vy5fvqzk5GQ1adJEzz//vLZu3VpkfmV9La43dsFrnp+fr88++0zffPONzp49qx9//FGrVq2Sl5dXqebbpEkTzZw5U88995zOnz+vmTNnyt3d/aZfvwIHDhyQj4+PqlQp/E+f672/xc3p999/V1JSklavXq2NGzcqIyND8fHxCgoK0ksvvVRk3B49emjdunXlrh8A4CIMAABMpF+/fsb9999fqO3o0aNGs2bNjDlz5ljbevXqZYSHhxuGYRjnzp0zGjVqZHz22WfWxzdu3Gg0adLEWLFihbFz506jc+fORl5enmEYhrF161bDYrEUGTshIaFI++jRo63jGIZh3HXXXUZISIhhGIZx6NChIv0HDhxovPPOO9b7w4YNM5YsWWKzeT/11FPGm2++ab1/7Ngxo0WLFsbJkyeN5cuXG++++671sbi4OGPGjBnXnV9ZXovrjW0YV15zX1/fQjXcddddxhtvvHHD+V8tLy/P6NKli+Hn52dcuHChTM8tSVRUlDFw4MAi7dd7fw2j5DndeeedRlJSkmEYV14Hi8ViHD9+vMj2d+zYYfTu3dsmcwAAOD9W0AEAprN//37ramROTo5ycnK0Zs0a9enTx9pn1KhR8vPzkyR5eHjojjvu0JYtWxQSEiJJqlOnjtLS0jRs2DDVrl270DHchmEUO269evUK3f/hhx+0cuVK7d+/39oWFhYmX1/fYp+/c+dOffXVV1q9erW1bdiwYVq6dKnGjh1b7nn/+OOPWrJkiU6fPm19TuPGjdWpUye98cYbaty4sSIjI3Xvvffq1ltvVVBQUKG+185PKv1rcaOx33zzTdWpU0eZmZkaOXJkoT6pqak3nPvVLBaLAgMDlZSUpBkzZmjevHllen5xTp48qTp16hRqK837W9ycGjVqpIyMDHXo0EGSFBAQIEk6dOiQ9ecC3t7eysrKKnf9AADXQEAHAJhO8+bNNXfu3Ov2mTx5svbv32/td/DgQbVv375Qn/r166t27do3XUdUVJQsFotatWplbRs6dGiJ/WNjY1W1alUtW7bM2nb69Gm1adOmVOPdaN6bN29W9erViwRNf39/fffdd/ruu+/02WefqUOHDmrYsKEeeughzZ49u1Rj38iNxr5ao0aNrD+7ubkpNze3TGNFREToscceU6NGjTR//nyNHj3aelb3m1XcFxFleX+vnpO7u7uaNGlivV9wssDiLofn5uZm05MJAgCcGwEdAFApTZ48WTt37tTKlSsVGBiolJSUIiGsVq1apdrW9u3bdccddxRpz8/Pl8ViKXGV+VoWi0W5ubkKCwurkFCWl5ennJwcGYZR6Izvly5dUl5ennJzcxUXF6fY2FhFR0dr4cKFOnLkiCIjI0s9RkmvxY3GtpXjx49r7969Gj9+vPr166f169frqaeeUmxsbLm2W69evSIr2WV9f69W2uecOXNG9evXL/P2AQCuia90AQCVTkxMjObPn69Vq1YpMDBQ0pWTfUlSdHS0Tp06dd3nF5x0rCBkbd26tdh+ffv2VX5+vn7++edC7d9//32x/fv16yfDMArtMi1JcXFx159QKd1xxx3Kz89Xenp6ofZDhw7p9ttv18KFC3Xw4EHddtttmjJlivbt22c9CVxJSvta3Gjs67n28nHX88orr1hX/T09PfWvf/1LO3fu1KJFi0q9jeL4+/sX2j1fKvv7e7XSzikzM9N6uTUAAG6EgA4AMBXDMJSfn3/dPmfPnpUk1ahRQ5KUlZWlH374QTk5OUpOTlZ+fv51t9OyZUu5u7srNTVV58+fl4+Pj3Xsq/Xs2VPDhw/XW2+9ZW3LysqyXjbr2v7dunVTSEiIIiIirG1HjhxRTEyMTebds2dPjRo1Sp9++qm17cCBA0pNTdXLL78sSVq4cKH1sWrVqqlr166FxrhWWV6L4sZOSUmxjl3c9vPz8284rwKLFy9W06ZNC604jxgxQvfee69eeuklZWRklGo7xWnTpo1OnDhRqJYbvb9S6eZU0Ke4ee7du/eGX2AAAFDAfcaMGTMcXQQAAF988YUmT56sbdu2KS0tTbt27VKtWrXUtm3bIn1btWqlqlWravHixTp69Kji4+M1ceJEffjhh2rXrp1OnTqlqVOn6pdfftG2bdvk5uamzp07W5/v6ekpb29vRUREaN++fZo4caK2bNmisLAwpaena8eOHbr99ttVt25dPfDAA4qPj9cnn3yivXv3Kj4+Xi+88IK2bdumCRMmKD09Xdu3b1fnzp3l5+enoUOHKi4uTmvXrlVsbKxSUlL0/PPPl7jiWpZ5S1eOkd6yZYs+//xzffvtt4qKitKSJUvUvHlzJSQkqHbt2tq8ebPi4+P1+eefa8aMGfL19dVXX31V7PzK8loUN/bHH3+s5s2ba/Xq1Zo6dapOnTql7du3q1u3bpo6dao2bNiglJQUpaamasiQIcXOacmSJRo/frwWL16spKQktWvXznrc/tKlS7Vy5UqdOXNGy5cv1969e3X33Xdbv5wpLYvFotjYWLVq1arQinZJ72/VqlVvOKeDBw+qRYsWCg0N1ZEjR5SYmChvb+9Cx8vPmzdP48aNU4MGDcpULwDANVmMmznwCgAAoJL56quvFBMTozlz5thlvMzMTA0YMEC7du2yy3gAgMqPXdwBAIBLGDBggPbs2aNLly7ZZbyIiAhNmjTJLmMBAJwDAR0AALiM6dOna+bMmRU+zpEjR/TTTz/pwQcfrPCxAADOg4AOAABcRq9evdSlS5ci1263tQULFuiDDz6o0DEAAM6HY9ABAAAAADABVtABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAEyCgAwAAAABgAgR0AAAAAABMgIAOAAAAAIAJENABAAAAADABAjoAAAAAACZAQAcAAAAAwAQI6AAAAAAAmAABHQAAAAAAE/h/3k/Kn5Gf+20AAAAASUVORK5CYII=\"><br>"
],
"metadata": {},
"output_type": "display_data",
"text": [
"<yt.visualization.particle_plotter.ParticlePlot at 0x10ba6eb50>"
]
}
],
"prompt_number": 1
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment