jzuhone/slice.ipynb Secret

## slice.ipynb
{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/jzuhone/Source/yt/yt/utilities/logger.py:4: VisibleDeprecationWarning: The configuration file /Users/jzuhone/.config/yt/ytrc is deprecated in favor of /Users/jzuhone/.config/yt/yt.toml. Currently, both are present. Please manually remove the deprecated one to silence this warning.\n",
      "Deprecated since v4.0.0. This feature will be removed in v4.1.0\n",
      "  from yt.config import ytcfg\n"
     ]
    }
   ],
   "source": [
    "import yt\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [INFO     ] 2021-11-16 16:55:33,393 Omega Lambda is 0.0, so we are turning off Cosmology.\n",
      "yt : [INFO     ] 2021-11-16 16:55:33,466 Parameters: current_time              = 2.5617741247968433\n",
      "yt : [INFO     ] 2021-11-16 16:55:33,468 Parameters: domain_dimensions         = [1 1 1]\n",
      "yt : [INFO     ] 2021-11-16 16:55:33,470 Parameters: domain_left_edge          = [0. 0. 0.]\n",
      "yt : [INFO     ] 2021-11-16 16:55:33,473 Parameters: domain_right_edge         = [40000. 40000. 40000.]\n",
      "yt : [INFO     ] 2021-11-16 16:55:33,473 Parameters: cosmological_simulation   = 0\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "re.compile('(output|group)_(\\\\d{5})') 4f1ccd3eccda8943b5d173a42d953fec None\n",
      "bye bye\n"
     ]
    }
   ],
   "source": [
    "ds = yt.load(\"snapshot_250.hdf5\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "c = [18982.4277343750000000, 20248.3164062500000000, 19998.5917968750000000]\n",
    "width = 900.0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "slc = ds.slice(2, c[2])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [INFO     ] 2021-11-16 16:55:33,802 Allocating for 2.972e+07 particles\n",
      "Loading particle index: 100%|██████████| 68/68 [00:00<00:00, 2698.60it/s]\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "(19798.31640625, 20698.31640625)"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmUAAAJBCAYAAAAHssIHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAA1qElEQVR4nO3df5TldX3n+ee7qm431RCohhUXq5sfUUTFH5Am2rskk0SOP2aGkY7CgKPCnnDW0TDniDEM4GQmsjsngTAjrsdds3NCIqjHtAMIiGEYF9y4o9AqdJO2bRmaRKCbPpJIF2B3Qd+qeu8f93vlUl1ddevWrbrf+73Pxzn31Lc+3x/1+XRV133V5/v5fL6RmUiSJKm3hnpdAUmSJBnKJEmSSsFQJkmSVAKGMkmSpBIwlEmSJJXASK8rsFQR4fRRSZLUVzIzZpfZUyZJklQCfd9T1nTSlXf1ugqSJKkDR64aZv/B6V5X4xeGI5jOZHxslAMHp9h3oN61az9+3bmH3WdPmSRJWlZxyI26l9t/cJra8AIHraDpTEZrw1zxrtOY6GIgW4ihrEuGh8rzwyRJUpkMLZTKgKmZZGy0RtDoqeq1yfo0n7pzR1t17xZDWZesLlHCrwL/NSWpOqZnFp6TlwnPTtb5wMYTmWnzEZCjtWFWjyxflJmYrDO9go+jrMyYsl47UJ/pdRUqxSm1kjR4EvjSA0+0ffwL9elKvV/YUyZJkvpSlQIZGMq6ZrTmP6UkScthbLTG2jW1Xldj2Xn7skuGIgiql9olSeq1bX/4TgBOueobpXifHQpoY5jcohnK5rGYf/Qyra8iSVKVnHHNf+Xct5zAULF+2FKN1oaZrHf2vr16ZIiZmTzsZIThoWhrYsNcKn/PbSm3FZcjBUuSpMWZmKzzpQeeaCuQfebCM+bdPz42yq+ceEzHdalPz1CfJyB0GshgAELZpLMiJUkaGNd8fce8+/e/OMV3H3um4+svZ4dN5UNZ2TlBQJI0yIZoPGapWxZ6JNLEZH3Fx6VFtLf+pmPKemgo7MmTJA22Geh4fFe/aHcYnN00PeSYNUmSfD9sMpRpINV8VqkkqWQMZRpI882ckSSpFwxlkiRJJWAokyRJWoJuDYgxlEnLYDgcsyZJZdetENTOgJh23hUMZaqs8bHRnn3tbjwGRJLUmbVram2tfbaSi1K1867gOmWqrD0Tk0S0vz6MJKkaFlpAtqzsKVOlGcgkScupm4NVDGVi7Zpar6sgSeoBR78uXTf/9q9MKPMHq3P92s0rSVoabyaUS2VC2Zi9PZIkqY9VJpT9/IWpXldBkiS1wbtbc6tMKPOxOZIk9QffsedWmVAmSZIGU204WFPr/0jjOmWSJKmv1aeT+nT/97/1f6yssCqkfkmSqmRstLZsT4zxXb/EDtRX8gEQkiRpPgFMTLa3jNTaNTWGhxY3pcFQNmBqi/wBkSRJDc0bpHsmJuc9rjYcZML0IichGsoGTTgVWZKk5VSfzrZ71FoZygZMfTpXdCryqOPiJElqy0C/Yw5HEMVHLY9Jx8VJktSWgQ1lAWz85bW8amyU6ez/abQraXxslKAxA6U2bKCVJKkbBnadsgS+89gzva5G37rhwjPYdOY4t2/dw/X3PLLgoEdJkjS/gQ1lK2UooGpPgNozMcnVt21/WVmzv6xiTZUkacUYypbZ8FAwAhyswErDrSbr0/ybr23nwMFpg5gkSV1gKFtmVXjsw+HsPzjd6ypIklQZAzvQf5A4FL9/OG9CkgaXoWwAVLevrnoq3LEqScuqCqsBGMokSVLfq8JwIUPZHFxMVpIkLUUnWcJQNofpTMdhSZKkjr3/betZu6a2qHMMZYcx4r+M1DafcSpJL/e1h/aw78DiHkrub9LD8JGN6je97N19oST/YRx6IGkl1IaDNQv8MdrJslGGMqkiejnEtSzDa//jP3+LQw8kLbv6dHJgGf4YNZRJqoQArvn6jtIERElaLEOZpEpIWPT4DUkqE0OZJEnSClhoUpShrI8FPkJJkqR+MbnAODRDWR9LujvAesiEJ0lSzxjKBMDwUDDjCGlJknrGUCYCmDaRSZLUUwuGsohYHxHfioidEbEjIj5WlB8bEd+MiEeLj2tbznlzRNxfHL89Io4oyjcUn++KiM9GNFZ6jIjVEbG5KN8SEScvU3s1B+OYpHatGnacg7Rc2ukpmwI+kZmvBzYCl0XEG4CrgHsz81Tg3uJzImIE+BLwkcw8HfhNoDlP/fPAh4FTi9e7i/JLgX2Z+RrgBuC6pTdNktRtB6f9M079oR/HSS8YyjJzb2Y+VGw/D+wExoHzgJuKw24CNhXb7wT+JjMfLs75WWZOR8QJwNGZeX9mJnBzyzmt17oFOKfZi9YqInL2a9EtnkMfft8kSdI8+nFUzqLGlBW3Fc8EtgCvzMy90AhuwPHFYa8FMiLuiYiHIuJfF+XjwO6Wy+0uypr7niyuNQU8Cxy36NZ0qA+/b5IkqWJG2j0wIo4CbgUuz8zn5ujIar3mrwG/ChwA7o2IB4Hn5ji2mYfmutghWSkz5+w9W7j21RQBObCtlySpWtrqKYuIGo1A9uXMvK0o/mlxS5Li49NF+W7grzPzHzLzAPBXwK8U5etaLrsOeKrlnPXFtUaAY4BnOm3UoDCQSZJUHe3MvgzgRmBnZn66ZdedwCXF9iXAHcX2PcCbI2JNEbB+A/hRcYvz+YjYWFzz4pZzWq91PnBfMe5MkiRpILRz+/Js4EPA9ojYVpR9ErgW+GpEXAo8AVwAkJn7IuLTwPdp3IL8q8z8RnHeR4EvAKPA3cULGqHvixGxi0YP2UVLa5YkSVJ/iX7vkGqOKTvpyrt6XRVJUpeN1oY4OJVM9/l7ldT0+HXnAnOPk297oP8gCZyRKUllsNADnKUq8TFLczCQSZKklWYokyRJKgFDmSRJUgkYyipmbLTW6ypIkqQOGMqWaE1tqDTPzgzgU+853W+qJPWhsryXqHecfblEB0o0M2gogmu+voPy1EiS1C4nmalSnSprapVqzqJNZ7LvQL3X1ZAkSR2oVIopU6+V1G9WDQdHrhrudTUkadGqcuu3UqFMUucOTicHp/3DRtLhlfUPt6rc+jWU9cCRq4YZjqrkelVJfboqv9okLYf9B6d7XYVKM5T1wP6D07z/bet7XQ1JklQihrIe+daP/77XVZAkSSViKCus9N3EpyYmGR8b7fj8kSFvf0qSVCWGskIkjNZWbgDjUAS/9bpX0Gm2mppx7I8kSVViKCvMAAenVm4A43Qmm7//JGYrSZIEhrKXWemJZ850kyRJTYYySZKkEjCUSZIklYChrIRe+Uurel0FSZK0wgxly6yTyZXPvTDNqccf2fW6SJKk8jKULaPxsVE+sPHERQezyfo0u57evyx1kiRJ5VTZUFaWtVXvenhvRw9KdV6mJEmDZaTXFVguZVj/a8/EZK+rIEmS+kRle8okSZL6iaFMkiSpBAxlkiRJJWAokyRJKgFDmSRJUgkYyiRJkkrAUNanXPFfkqRqMZR1QS8Wqt297wWOXDW88l9YkiQti8ouHruSerFQ7WR9euW/qCRJWjb2lEmSJJWAoUySJKkEDGUqpZI8T16SpBVjKFtGI72YAVARJXievCRJK8pQtoymejEDQJIk9SVDmSRJUgkYyiRJkkrAUCZJklQChjJJklRZ/RR0+qmukiSpYo5cNcxobfkeGzizbFfuPh+z1EVD0XjkUgSkEy8lSVrQ/oM+NrDJnrIuaq6AYSCTJEmLZSgbcC5wK0lSORjKBpwL3EqSVA6GMkmSpBIwlEmSJJWAoUySJKkEDGWSJEklsGAoi4j1EfGtiNgZETsi4mNF+bER8c2IeLT4uLYoPzkiJiNiW/H605ZrbYiI7RGxKyI+GxFRlK+OiM1F+ZaIOHmZ2itJklRK7fSUTQGfyMzXAxuByyLiDcBVwL2ZeSpwb/F502OZeUbx+khL+eeBDwOnFq93F+WXAvsy8zXADcB1S2lUbdhlHiRJUn9ZMJRl5t7MfKjYfh7YCYwD5wE3FYfdBGya7zoRcQJwdGben5kJ3NxyTuu1bgHOafaizbpGzn7N9bXq0y7zIEmS+suixpQVtxXPBLYAr8zMvdAIbsDxLYeeEhFbI+KvI+LXi7JxYHfLMbuLsua+J4trTQHPAsctrindZ3+bJElaKW0/+zIijgJuBS7PzOfm6Mhq2gucmJk/i4gNwO0RcTpzZ5xml9Z8+14qyJyz96yd+nfC/jZJkrRS2uopi4gajUD25cy8rSj+aXFLsnlr8mmAzHwxM39WbD8IPAa8lkbP2LqWy64Dniq2dwPri2uNAMcAz3TerJVjb5okSeqGdmZfBnAjsDMzP92y607gkmL7EuCO4vhXRMRwsf3LNAb0/21xi/P5iNhYXPPi5jmzrnU+cF8x7qz0EhgfG+11NSRJUp9r5/bl2cCHgO0Rsa0o+yRwLfDViLgUeAK4oNj3j4D/LSKmgGngI5nZ7PX6KPAFYBS4u3hBI/R9MSJ20eghu2gJbVpxeyYmCbzdKUmSOhd90iF1WM0xZSddeVevqyJJkjSvx687F5h7nLwr+kuSJJWAoUySJKkEDGWSJEklYCiTJEkqAUNZBwIYG631uhqSJKlCDGUdSKA+PePCsZIklVxtqH/erdt+zJJebv/B6V5XQZIkLaA+0z9Lf9lTJkmSVAKGMkmSpBIwlEmSJK2AU48/ct79hjJJkqQV8OjT++fdbyiTJEkqAUOZJElSCRjKJEmSSsBQJkmSVAKGMkmSpBIwlEmSJJWAoUySJKkEDGWSJEklYCiTJEkqAUOZJElSCRjKJEmSSsBQ1qHhiF5XQZIkVYihrANr19SYyex1NSRJUoUYygpDi+j4yoRXjY0uX2UkSdLAMZQBn7nwDI4ZrbV9/MRknYkDB5exRpIkadCM9LoCZXD1bduZrE8v6pz9Bxd3vCRJ0nzsKYNFBzJJkqRuM5RJkiSVgKFMkiSpBAxlXbKmNsRozX9OSZLUGVNEl9Snkz9+75v54MYTe10VSZLUh5x92SX1meT6ex5h/4tTva6KJEnqQ4ayLtozMdnrKkiSpD7l7csV5NMyJUnS4RjKVlACPsdckiTNxVDWoU6zlc8xlyRJczGUdchsJUmSuqnyocy7hZIkqR9UPpQlMBRQG+puPDPsSZKkbqp8KAOYSTjqiJGuBilvX0qSpG4aiFAGMHGgztiaWq+rIUmSNKeBCWXHjNaYOFDvdTUkSdIAWD2y+Ig1MKFs/8EpasPt38B0zJgkSerUi1Mziz5nYEJZfTo5ON3eSLDR2rBjxiRJ0oqqTChb28XxYkG68r4kSVpRlXkg+b4ujhc7UF98l6MkSdJSVKanbCUcuWq411WQJEkVZShbhNqw/1ySJGl5mDIW4dnJOqsWMYNTkiSpXYayQjtR65jRWtszOCVJkhZjwVAWEesj4lsRsTMidkTEx4ryYyPimxHxaPFx7azzToyIn0fE77eUbYiI7RGxKyI+G9GY4xgRqyNic1G+JSJO7nI7528j8D+/+tgFg9nBqemVqI4kSRpA7fSUTQGfyMzXAxuByyLiDcBVwL2ZeSpwb/F5qxuAu2eVfR74MHBq8Xp3UX4psC8zX1Ocd10HbelIAB/YeCI/+dnkgmuTOStTkiQtlwVDWWbuzcyHiu3ngZ3AOHAecFNx2E3ApuY5EbEJ+FtgR0vZCcDRmXl/ZiZwc8s5rde6BTin2Yu23FaNDHHXw3vZMzG5El9OkiRpTosaU1bcVjwT2AK8MjP3QiO4AccXxxwJXAlcM+v0cWB3y+e7i7LmvieLa00BzwLHzfH1c/ZrMfWfy4tTM0xM+kxMSZLUW22Hsog4CrgVuDwzn5vn0GuAGzLz57MvMcex2cY+SZKkymtrRf+IqNEIZF/OzNuK4p9GxAmZube4Nfl0Uf424PyI+BNgDJiJiBeK89e1XHYd8FSxvRtYD+yOiBHgGOCZ2fXIzEPCWzd6yyRJknqtndmXAdwI7MzMT7fsuhO4pNi+BLgDIDN/PTNPzsyTgc8Af5SZnytucT4fERuLa17cPGfWtc4H7ivGnUmSJA2EdnrKzgY+BGyPiG1F2SeBa4GvRsSlwBPABW1c66PAF4BRGjMzm7MzbwS+GBG7aPSQXdRm/SVJ0gqqDQX1GftNlkP0e4dU8/blSVfe1euqSJIkzevx684F5h6S5Yr+kiRJJWAokyRJKgFDmSRJUgkYyiRJkkrAUCZJklQChjJJkqQSMJRJkiSVQGVC2VwPz5QkSeoXlQllrxob7XUVJEmSOlaZULZnYrLXVZAkSepYZUKZJElSPzOUSZIklYChTJIkqQQMZfMYG631ugqSJGlAGMrmEWEwkyRJK2Ok1xUos30H6r2ugiRJGhD2lEmSJJWAoUySJKkEDGWSJEklYCiTJEkqAUOZJElSCRjKJEmSSsBQJkmSVAKGMkmSpBIwlLWoDQW14eh1NSRJ0gAylBUCeOspa5mazl5XRZIkDSAfs1QYGYLvPPZMr6shSZIGlD1lhfpMr2sgSZIGmaFMkiSpBAxlNAb4O7xfkiT1kqEMqA0HDu+XJEm9ZCgDDjigTJIk9ZihTJIkqQQMZZIkSSVgKJMkSSoBQ5kkSVIJDGwo8xGXkiSpTAY2lDnfUpIklcnAhrJ0YTJJklQiAxvKJEmSysRQJkmSVAKGMkmSpBIwlEmSJJWAoUySJKkEDGWSJEklYCiTJEkqAUOZJElSCRjKJEmSSsBQJkmSVAKGMkmSpBIwlEmSJJWAoUySJKkEFgxlEbE+Ir4VETsjYkdEfKwoPzYivhkRjxYf1xblb42IbcXr4Yj47ZZrbYiI7RGxKyI+GxFRlK+OiM1F+ZaIOHmZ2itJklRK7fSUTQGfyMzXAxuByyLiDcBVwL2ZeSpwb/E5wA+BszLzDODdwP8dESPFvs8DHwZOLV7vLsovBfZl5muAG4DrltowSZKkfrJgKMvMvZn5ULH9PLATGAfOA24qDrsJ2FQccyAzp4ryI4AEiIgTgKMz8/7MTODm5jmzrnULcE6zF61VROTs1yLbK0mSVEqLGlNW3FY8E9gCvDIz90IjuAHHtxz3tojYAWwHPlKEtHFgd8vldhdlFB+fLK41BTwLHNdBeyRJkvpS26EsIo4CbgUuz8zn5js2M7dk5unArwJXR8QRwCE9XxS9aAvsa71uzH7NV4/R2hBD8x4hSZIGxQc3ntjrKsyrrVAWETUagezLmXlbUfzT4pZk89bk07PPy8ydwH7gjTR6xta17F4HPFVs7wbWF9caAY4BnllsY2abrM8w4w1OSZIG3vjYKGeddGyvqzGvdmZfBnAjsDMzP92y607gkmL7EuCO4vhTmgP7I+Ik4DTgJ8UtzucjYmNxzYub58y61vnAfcW4M0mSpCWpDQUHDk5x+eZtva7KvEYWPoSzgQ8B2yNiW1H2SeBa4KsRcSnwBHBBse/XgKsiog7MAL+bmf9Q7Pso8AVgFLi7eEEj9H0xInbR6CG7aAltkiRJ+oX6TLLvQL3X1VhQ9HuHVHMG5klX3tXrqkiSJM3r8evOBRrj5Gfvc0V/SZKkEjCUSZIklYChTJIkqQQMZZIkSSVgKJMkSSoBQ5kkSVIJGMokSZJKwFAmSZJUAoYySZKkEjCUSZIklYChTJIkqQQMZZIkSSVgKJMkSSoBQ5kkSVIJGMokSZJKwFAmSZJUAgMfyoai1zWQJEkakFA2X/CayZWrhyRJ0uEMRCgzeEmSpLIbiFAmSZJUdoYySZKkEjCUSZIklYChTJIkqQQMZZIkSSVgKJMkSSoBQ5kkSVIJGMokSZJKwFAmSZJUAoYySZKkEjCULWAo4OxXH4vPLZckScvJULaAmYQH/m4fN1x4Bmtq/nNJkqTlMdLrCvSD6Znk8s3bCLvLJEnSMrHrZxEye10DSZJUVYYySZJUCaO1YUb7eKhR/9ZckiSpEMD7NowzWZ/pdVU6ZijrgEPLJEkqlwRufXBPr6uxJIayDozWhvjgxhN7XQ1JktRisj7d6yosiaGsAwfqM5x10rGc/epje10VSZJUES6J0aHLN2/zNqYkSeoae8qWwBUyJElStwxMKBsbrfW6CpIkSYc1EKFsOIJnJ+u9roYkSdJhDcSYsmmX4pckSSU3ED1lkiRJZWcokyRJKgFDmSRJUgkYyiRJkkrAUCZJklQCAx3KOl2Rf+2aGuNjo0SxXRtybX9JkrQ0A7EkxlwC+MDGE/nG3+xl34FD1zAbCpjJxnGtC2qM1ob5w392OpvOHP9F2e1b93D9PY/w1MQkY2tqvFif5kB9ZtnbIEmSqmOgQlkzYI2PjXLFu05j05njnHXSsVxxy8PUp1+KXrXh4Prz38KmM8dfFrhe1XJeq01njh9SdvvWPVzz9R1zBj5ohLv3bRjn1gf39P1T7SVJ0tJFLrCwakSsB24G/kdgBvhPmfl/RMSxwGbgZOAnwD/PzH0R8Q7gWmAVcBC4IjPvK661AfgCMAr8FfCxzMyIWF18jQ3Az4ALM/MnbTUgIgFOuvKuufdzaBCbrZ3gtVSH+xqt5ceM1ojgsEFOkiT1t8evOxeAzDxk7FM7oewE4ITMfCgifgl4ENgE/C/AM5l5bURcBazNzCsj4kzgp5n5VES8EbgnM8eLa30P+BjwAI1Q9tnMvDsifhd4c2Z+JCIuAn47My9sp3ELhTJoBLLvXPX2di5XCidf9Y1eV0GSpL42e/hRWcwXyhYc6J+ZezPzoWL7eWAnMA6cB9xUHHYTjaBGZm7NzKeK8h3AERGxugh3R2fm/dlIgjc3z5l1rVuAcyKia6Pnn5qY7Nallt3tW/d0PAFBkiQ1lDGQLWRRsy8j4mTgTGAL8MrM3AuN4AYcP8cp7wO2ZuaLNILc7pZ9u4syio9PFteaAp4Fjpvj6+fsVzv1ftXYaDuHlcL19zzSlz9IkiRpadoe6B8RRwG3Apdn5nMLdWRFxOnAdcA7m0VzHJZt7FuyK951Wrcutez6qVdPkiR1T1s9ZRFRoxHIvpyZtxXFPy1uSTbHnT3dcvw64GvAxZn5WFG8G1jXctl1wFMt+9YX544AxwDPzK5HZsbs10J1HxutdX3Q/nLqp149SZLUPQuGsmJs143Azsz8dMuuO4FLiu1LgDuK48eAbwBXZ+Z3mgcXtzifj4iNxTUvbp4z61rnA/flQjMQ2jBaG+ZT7zl9qZdZUVe86zRGa8O9roYkSVph7fSUnQ18CHh7RGwrXv+ExrIX74iIR4HmMhgA/wp4DfBvW45vjjf7KPBnwC7gMeDuovxG4LiI2AX8HnDVUhsWwPs2HLp+WNltOnOcP37vm1i7ptbrqkiSpBW04JIYZTffkhhr19TY+u/eeUh5Pzj72vvY4/gySZIqZUlLYvSzfQfq3L51T6+r0REH/EuSNFgqHcqgscREP3LAvyRJg6Xyoaxfe5x+63WvcBFZSZIGSOVDWT/2ON2+dQ+bv/+ki8hKkjRAKh3KRmvDfbVwbNM1X99BfdpIJknSIKl0KPvj976p75bEgMYEBUmSNFgqG8r6bSV/SZI02CobyurTM72uQscWeKyoJEmqoMqGsv0Hp/t2jbI+X89XkiR1oLKhDPpzjbJ+DZKSJGlpKh3K+nGNsk/duaPXVZAkST1Q6VB2RK3/mjcx6cxLSZKWImZ97Bf9l1oWYbI+wx/cvr3X1Wibty4lSVq6BIYj+m4R9kqHMoCvbHmy11VoWz+OgZMkqYym+3DWXOVDWT99U/pxDJwkSeqOyoey4T5a9Ksfn9MpSVLZ9M87/8tVPpS9/23re12Ftl3xrtMYrQ33uhqSJPW1/rlH9nIjva7Acjpy1TD/ftObel2NtjUfC3X55m29rYgkSVpxle4pO3BwutdVWLRNZ44z1K/9rpIkqWOVDmX9OEbr9q17mOnXfldJktSxyoay0dowV7zrtF5XY9FcFkOSpMFU2VD2Kyce84sxWv1kj8tiSJI0kCobyr772DN9uUJ+Py3hIUmSuqeyoSzpz1uB/bTYrSRJ6p7KhjLozxXyx/twcoIkSVq6Soeyfpx96QKykiQNpsouHtuvsy9dQFaSpMFU2Z6y1SP927RNZ457G1OSpAHTv8llAROTda6+bXtfzsAE+K3XvaLXVZAkSSuosqEMYLI+3ZczMAFufXB3r6sgSZJWUKVDGfTnDMzbt+5hsj7T62pIkqQVVPlQ1o8zMPu1d0+SJHWu0qEsoC9nYPZj754kSVqayi6JAXBEbYiPb97G9fc8whXvOq0rz8K8feserr/nEZ6amORVY6OHvW47x92+dQ/XfH0H+w7UARgbrfGp95zOq8ZGfQamJEkDJrLPH+sTEQlw0pV3tXX82jU1MuHZyfq8oWq227fu4VN37mBisv6y8tHaMH/83je97Bq3b93D1bdtZ7I+/VI9gQ9sPJF/v+lNvzjmilsepj596L9/bQgcUiZJUvU8ft25AGTmIQ+7rnRP2VyavVIAeyYm+fjmbfzg8Wd+EZaaZvdiHU5zhmczlN2+dQ+f+OrDhzzDMoEvP/AEZ510LJvOHOf6ex6ZM5CBgUySpEE0cD1lhzNaG+KI2jATB+qMranx7GSdmUX804yPjfJbr3sFtz6452U9ZHMd952r3s4pV32D/v6XlyRJi2VPWRsm6zO/WIZiod6xueyZmOTLDzyxYNDaMzHJ7Vv3MFob4oBdYpIkqWAo66J2e76uvm07LxjIJElSC0NZD8x3e1OSJA2mSq9TJkmS1C8MZT2ypuY/vSRJeonJoAcC+KP3vpmzX31sr6siSZJKwjFlPbBm1TCXb97GIXNhJUnSwLKnbIXVhoP9BxsD/V2nTJIkNRnKVtjhVvGXJEmDzVAmSZJUAoYySZKkEjCUSZIklYChTJIkqQQMZT0WroshSZIwlPVcOhlTkiTRRiiLiPUR8a2I2BkROyLiY0X5sRHxzYh4tPi4tig/rjj+5xHxuVnX2hAR2yNiV0R8NqLRTxQRqyNic1G+JSJOXoa2SpIklVY7PWVTwCcy8/XARuCyiHgDcBVwb2aeCtxbfA7wAvBvgd+f41qfBz4MnFq83l2UXwrsy8zXADcA13XWHEmSpP60YCjLzL2Z+VCx/TywExgHzgNuKg67CdhUHLM/M/8bjXD2CxFxAnB0Zt6fmQnc3Dxn1rVuAc5p9qLNukbOfi2msZIkSWW1qDFlxW3FM4EtwCszcy80ghtw/AKnjwO7Wz7fXZQ19z1ZXGsKeBY4bjF1kyRJ6mdtP5A8Io4CbgUuz8zn5ujIWvASc5RlG/teKsics/dssRWRJEkqm7Z6yiKiRiOQfTkzbyuKf1rckmzemnx6gcvsBta1fL4OeKpl3/riWiPAMcAz7dRNkiSpCtqZfRnAjcDOzPx0y647gUuK7UuAO+a7TnGL8/mI2Fhc8+KWc1qvdT5wXzHuTJIkaSC0c/vybOBDwPaI2FaUfRK4FvhqRFwKPAFc0DwhIn4CHA2siohNwDsz80fAR4EvAKPA3cULGqHvixGxi0YP2UVLaZQkSVK/WTCUFTMpDzeA7JzDnHPyYcp/ALxxjvIXaAl1kiRJg8YV/SVJkkrAUCZJklQChjJJkqQSMJRJkiSVgKFMkiSpBAxlkiRJJWAokyRJbVv0QxbVNkOZJElqm4/bWT6GMkmSpBIwlEmSJJWAoUySJKkEDGWSJEklYCiTJEkqAUOZJElSCRjKJEmSSsBQJkmSVAKGMkmSpBIwlEmSJJWAoUySJKkEDGWSJEklYCiTJEkqAUOZJElSCRjKJEmSSsBQJkmSVAKGMkmSpBIwlEmSJJWAoUySJC1o9YiRYbn5LyxJkhb04tRMr6tQeQMdyoai1zWQJKm/rRr2zbRbBjqUzWSvayBJUn87OO2babcMdCiTJEkqC0OZJElSCRjKJEmSSsBQJkmSKif6cP6BoUySJFVO9uH8A0OZJElSCRjKJEmSSsBQJkmSVAKGMkmSKmgooOZq+31lpNcVkCRJ3TeTMONq+33FnjJJkqQSMJRJkiSVgKFMkiSpBAxlkiRJJWAokyRJKgFDmSRJUgkYyiRJkkrAUCZJklQChjJJkqQSMJRJkiSVgKFMkiSpBBYMZRGxPiK+FRE7I2JHRHysKD82Ir4ZEY8WH9e2nHN1ROyKiEci4l0t5RsiYnux77MREUX56ojYXJRviYiTl6Gt1IygkiSppNqJKVPAJzLz9cBG4LKIeANwFXBvZp4K3Ft8TrHvIuB04N3A/xURw8W1Pg98GDi1eL27KL8U2JeZrwFuAK7rQtsOUZ9ZjqtKkiQt3YKhLDP3ZuZDxfbzwE5gHDgPuKk47CZgU7F9HvCXmfliZv4dsAt4a0ScABydmfdnZgI3zzqnea1bgHOavWitIiJnvxbdYkmSpBJa1A294rbimcAW4JWZuRcawQ04vjhsHHiy5bTdRdl4sT27/GXnZOYU8Cxw3GLqJkmS1M9G2j0wIo4CbgUuz8zn5ujI+sWhc5TlPOXznfPygsw5e88OVxFJkqR+0VZPWUTUaASyL2fmbUXxT4tbkhQfny7KdwPrW05fBzxVlK+bo/xl50TECHAM8MxiGyNJktSv2pl9GcCNwM7M/HTLrjuBS4rtS4A7WsovKmZUnkJjQP/3ilucz0fExuKaF886p3mt84H7inFnkiRJA6Gd25dnAx8CtkfEtqLsk8C1wFcj4lLgCeACgMzcERFfBX5EY+bmZZk5XZz3UeALwChwd/GCRuj7YkTsotFDdtHSmiVJktRfot87pJpjyk668q5eV0WSJGlej193LjD3OHmXU11Gh58LIUmS9HKGsmXU552QkiRpBRnKJEmSSsBQJkmSVAKGMkmSpBIwlEmSJJWAoUySJKkEDGWSJEklUPlQNuxiYT3nd0CSpIVVPpRNu1hYz/kdkCRpYZUPZeo9eyslSVqYoUzLzt5KSequAM5+9bG9rkbfOfX4I3tdhXkZyiRJ6jMJ/ORnk72uRt959On9va7CvAxlkiT1oT0Tk06kqhhDmSRJfcrBIdViKNOSfHDjiQz5p5okSUtmKNOSfOvHf89obbjX1ZAkqe+N9LoC6m97JhxoKklSN1Smp2xNrTJNkSRJA6gySeZAfabXVZAkSepYZUKZJEkryTlO6jZDmSRJHfAJcuo2Q9mA8nmUkrQ0My4Spi4zlA0on0cpSUvj37bqNkOZJEmd8G9bddnAhbJu/2HjH0qSNJgSWD0ycG+jWkYD99PU7T9s/ENpcUZdT05Shbw45XJM6h7fIbWiJl1PTpKkORnKJEmSSsBQJkmSVAKGMkmSpBIwlEmStMycqa92GMpKxv+4klQ9ztRXOwxlJVIbgrE1tV5XQ5Ik9YChrETqM7DvQL3X1ZAkST1gKJMkSSoBQ5kkSRK9H9dtKFMp+TQmSdJK6/WEDN/6Btiq4V7/TXB4Po1JkjRoDGUD7OB0r/8mkCRJTYYySZKkEjCUSZIklYChTJIkqQQGIpSVdzi7JElSw0CEMoezS5KkshuIUCZJklR2hjJJUunVhoPPXHgGn7nwjF5XRerYQsOpDGWSpNKrTyfX3/MIP3j8mV5XRerYQsOpDGWSpL6wZ2KSLz/wRK+rIS0bQ5kkqW84cUtVZiiTJEkqgQVDWUT8eUQ8HRE/bCl7S0TcHxHbI+LrEXF0Ub4qIv6iKH84In6z5ZwNRfmuiPhsRERRvjoiNhflWyLi5K63UpIkqeTa6Sn7AvDuWWV/BlyVmW8CvgZcUZT/rwBF+TuA/xgRza/xeeDDwKnFq3nNS4F9mfka4Abguo5aIkmS1McWDGWZ+W1g9nSX04BvF9vfBN5XbL8BuLc472lgAjgrIk4Ajs7M+zMzgZuBTcU55wE3Fdu3AOc0e9EkSZIGRadjyn4IvKfYvgBYX2w/DJwXESMRcQqwodg3DuxuOX93UUbx8UmAzJwCngWOm+uLRkTOfnVY/yVbU3M4niRJ6p5Ok8XvAJdFxIPALwEHi/I/pxG4fgB8BvguMMXc66U1A9V8+0rrQH0G+/MkSVK3jHRyUmb+GHgnQES8FvinRfkU8PHmcRHxXeBRYB+wruUS64Cniu3dNHrTdkfECHAMh94ubX7dQ2JQL3vLsvTRUZIk9YuOesoi4vji4xDwB8CfFp+viYgji+13AFOZ+aPM3As8HxEbi/FiFwN3FJe7E7ik2D4fuK8YdyZJkjQwFuwpi4ivAL8J/A8RsRv4Q+CoiLisOOQ24C+K7eOBeyJiBtgDfKjlUh+lMZNzFLi7eAHcCHwxInbR6CG7aAnteXnd6YP7oJIkSUD0e6dU8/blSVfedci+D248kVsf3MNkfXrF6yVJkjTb49edC8w9JKujMWX94ks+I02S+kptOKhP93dngdQp13WQJJXGW09ey3AxtX04wuWHNFD8aZcklcZ3H3uG6WJYzXQmB+ozPa6RtHhr19Q6Os9QJkkqDW9cqh8sFJ4y516EdSGVHlMmSZLUbQv1305M1ju6rj1lkiRJJWAokyRJXTHk4weXxFAmdcDfO5J0qBkHBS6JoUxapPGxUQcjS5K6zlC2DMJulErbMzHZ6ypIkirIULYM+vzJVZIkqQcMZVKXOMBVkrQUhjKpSxzgKklaCkNZm45ePdzrKkiV52MOV84HN57Y6ypImqWyvwK7fSfpuRenu3xFSbNN+ZjDFfOlB57odRUkzVLZUOadJKn/+P9W0iCrbCiT+tXaNTVqzhpQD3Xy8zc+NroMNZEGi6FMKpk//Genc+Fb1zPsgndaYQF85sIzuP6Ctyz65++3XveK5amU1GVrakOl/cN3pNcVKJvacBDAwWlvpGhxhoBuDIn6va9uYziCaRe80wr7wMYT2XTmOAAf37xtUefe9fBe1q6pse9AfRlqJnXPgXp5B68aylqsqQ2xujbsLxV1pFv/zWcSZgxkfSGo1ji4Lz3wBF/Z8iTTmQzF4hbCnpiss8bps9KSGMoKtaGgPpMcMJDpMKr2Bqylq+LPQ7OHtpN198rcAyH1A/+sKdRnkrq3LDWPBMZGa72uhiSpogxlUpsCOPctJ/S6GpKkijKUlVQ554UMtqTaC276MydJvWUoKylvpKqbFgpcAdxw4RkrUBNJGgydBKyBCWX2AmiQDS+0Jk+xe+2alRszt5JfS5JWWifTXgZm9qU9TxpkUwtMpcuEK255eEUnuxycmnFGqyS1GJieMqnK1q6pLXlm6ErPPt5/cLoygcyHL0jqBkOZ1OfGRmts/Xfv5NlJ19jrFdf6lQbTkauGu3o9Q5nU5yYm65x97X2MOUZLklbU/oPTXb3ewIwp0/JzfFDv7JmY7HUVJElLZE+ZusZAJklS5wxlkrQCnAsgaSGGMklaAfYkS1qIoWwFnP3qY32QtSRJmpehbAV87+/2ce5bTmC01t2ps5LUK7XhYKEHRUiDrja0uP8nhrIVUJ9JvvXjv+d9G8YdVzJg/H6rqurTyQIPipBKo1e/i+szi/t/YihbIXsmJvnSA084rmTA+P3uHQOxpKZ++V1sKJPUc8PL8JyiI2pD3l6rCL+NGhSGMpb2hhB05xdGN3/pjI3WGK35rVX/mF6G5xRN1me8vVYRfhs1KFzRH1g9ErxQT2Y6OLdbvyy6+UtnwmcgSpLUdwxlwIF6+3HMRwlJkqTl4D2uRTKQSZKk5WAokySphHyDHjx+zyVJKqFOxjmrvxnKJEmSSsBQJkmSVAKGMkmSpBIwlBWWY8VoV6GWJEntMpQVEli7ptbdi5rKJElSmwxlLfYd6O5K+Mvw5BhJklRRhjJJkqQSWDCURcSfR8TTEfHDlrK3RMT9EbE9Ir4eEUcX5bWIuKko3xkRV7ecs6Eo3xURn41oPAU8IlZHxOaifEtEnLwM7ZQkSSq1dnrKvgC8e1bZnwFXZeabgK8BVxTlFwCri/INwL9sCVmfBz4MnFq8mte8FNiXma8BbgCu66glFbaYsW4OY5MklVGV3p9qw8GaWvdvNi54xcz8NvDMrOLTgG8X298E3tc8HDgyIkaAUeAg8FxEnAAcnZn3Z2YCNwObinPOA24qtm8Bzmn2os0WETn71dw3PPcp8+qXH5DFjHVzGJskqYyq9P505KoRJuvdf+bCSIfn/RB4D3AHjd6x9UX5LTRC1l5gDfDxzHwmIs4CdrecvxsYL7bHgScBMnMqIp4FjgP+YTEV+ttr/2lnLZEkSVqEx5fpup32vf0OcFlEPAj8Eo0eMYC3AtPAq4BTgE9ExC8zd6dUMzTPt+/lhZkx+9Vh/SVJkkqlo56yzPwx8E6AiHgt0Oym+hfAf8nMOvB0RHwHOAv4/4B1LZdYBzxVbO+m0dO2u7jteQyH3i6dry6lD2bN26z9UNdust22exDYbts9CGz3yrS7o56yiDi++DgE/AHwp8WuJ4C3R8ORwEbgx5m5F3g+IjYW48UupnHrE+BO4JJi+3zgvmLcmSRJ0sBoZ0mMrwD3A6dFxO6IuBR4f0T8d+DHNHq8/qI4/P8EjqIx5uz7wF9k5t8U+z5KY9bmLuAx4O6i/EbguIjYBfwecFU3GiZJktRPwk6p5We3r+0eBLbbdg8C2227l5Mr+kuSJJWAPWWSJEklYE+ZJElSCRjKJEmSSsBQJkmSVAKGskWKiD+PiKcj4octZWdExAMRsS0ifhARb23Z9+aIuD8idkTE9og4oijfUHy+KyI+23zeZ0SsjojNRfmWlge699Ri2h0RHyjKmq+ZiDij2Ffldtci4qaifTsj4uqWc6rc7lUR8RdF+x6OiN9sOacK7X5L8X94e0R8PSKObtl3ddGGRyLiXS3llW13RBwXEd+KiJ9HxOdmXafK7X5HRDxYlD8YEW9vOafK7X5rvPS7/OGI+O2Wcyrb7pb9JxY/67/fUrZ87c5MX4t4Af8I+BXghy1l/xX4x8X2PwH+32J7BPgb4C3F58cBw8X294D/icZjpu5uOf93gT8tti8CNve6zYtt96zz3gT8bcvnlW03jSda/GWxvQb4CXDyALT7MhprEgIcDzwIDFWo3d8HfqPY/h3gfy+23wA8DKym8Vi5xyr2//tw7T4S+DXgI8DnZl2nyu0+E3hVsf1GYM+AtHsNMFJsnwA83fJ5Zdvdsv9W4D8Dv78S3297yhYpM7/NoY+BSqCZro/hpUdIvRP4m8x8uDj3Z5k5HREnAEdn5v3Z+O7dDGwqzjkPuKnYvgU4p5nCe2mR7W71fuArAAPQ7gSOjMbjwkZpPBP2uQFo9xuAe4vzngYmgLMq1O7TgG8X298E3ldsn0cjhL+YmX9HY2Hst1a93Zm5PzP/G/BC68ED0O6tmdn8md8BHFH0jFS93Qcyc6ooP4Li2dRVbzdARGwC/pbG97tZtqztNpR1x+XA9RHxJPAfgOZtq9cCGRH3RMRDEfGvi/JxGs/8bNpdlDX3PQlQ/Ed4lkYPWxldztztbnUhRSij+u2+BdgP7KXxyLH/kJnPUP12PwycFxEjEXEKsIHG82yr0u4fAu8pti+g0TZoaUOh2b6qt/twBqnd7wO2ZuaLDEC7I+JtEbED2A58pGhLpdsdjUdFXglcM+v4ZW23oaw7Pgp8PDPXAx+n8egoaNy+/DXgA8XH346Ic2h0ec7WXDBuvn1lc7h2A43/yMCBzGzev696u98KTAOvonE76xMR8ctUv91/TuMX0w+AzwDfBaaoTrt/B7gsIh4EfolGDygcvg1Vb/fhDES7I+J04DrgXzaL5rhGpdqdmVsy83TgV4GrozE2uurtvga4ITN/Puv4ZW23oaw7LgFuK7b/M403Z2i8Uf11Zv5DZh4A/orG/ezdwLqW89fx0q2g3byU1Edo3Caa3d1aFodrd9NFvNRLBtVv978A/ktm1ovbeN8BzqLi7c7Mqcz8eGaekZnnAWPAo1Sk3Zn548x8Z2ZuoPHz/Fix6xdtKDTbV/V2H07l2x0R64CvARdnZuvPQaXb3XLMThp3A95I9dv9NuBPIuInNO4SfDIi/hXL3G5DWXc8BfxGsf12Gm9IAPcAb46INcU36DeAH2XmXuD5iNhY3G++GLijOOdOGm9+AOcD9xX3rcvocO0mIoZodAX/ZbNsANr9BPD2aDgS2Aj8uOrtLn6+jyy23wFMZWZlfs4j4vji4xDwB8CfFrvuBC4qxhWdApwKfG8A2j2nqrc7IsaAbwBXZ+Z3mscPQLtPKd6/iIiTaIzB+knV252Zv56ZJ2fmyTTuAPxRZn5u2dudJZgR0U8vGkl6L1CnkYovpXFr8kEaY2u2ABtajv8gjUGCPwT+pKX8rKLsMeBzvPTIqyNo9ELsojHD45d73eYO2/2bwANzXKey7QaOKtqwA/gRcMWAtPtk4BFgJ/D/ACdVrN0fA/578bq22Ybi+H9TtO0RihlYA9Lun9DoAfh5cfwbqt5uGm/Y+4FtLa/jB6DdH6LxO20b8BCwaVB+zlvO+xQvn325bO322ZeSJEkl4O1LSZKkEjCUSZIklYChTJIkqQQMZZIkSSVgKJMkSSoBQ5kkSVIJGMokSZJKwFAmSZJUAv8/+pTr0SuQNwEAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 720x720 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(10,10))\n",
    "plt.scatter(slc[\"PartType0\",\"particle_position_x\"], slc[\"PartType0\",\"particle_position_y\"])\n",
    "plt.xlim(c[0]-0.5*width, c[0]+0.5*width)\n",
    "plt.ylim(c[1]-0.5*width, c[1]+0.5*width)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "interpreter": {
   "hash": "e8df9565fcf22429a29547d2a7f7c9f690d5e8ce8cb391fd548802b5f7834ef2"
  },
  "kernelspec": {
   "display_name": "Python 3.9.2 64-bit ('base': conda)",
   "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.9.2"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}