{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\u001b[1m[\u001b[0m\u001b[34;49mnotice\u001b[0m\u001b[1;39;49m]\u001b[0m\u001b[39;49m A new release of pip available: \u001b[0m\u001b[31;49m22.3.1\u001b[0m\u001b[39;49m -> \u001b[0m\u001b[32;49m23.1\u001b[0m\n",
      "\u001b[1m[\u001b[0m\u001b[34;49mnotice\u001b[0m\u001b[1;39;49m]\u001b[0m\u001b[39;49m To update, run: \u001b[0m\u001b[32;49mpip install --upgrade pip\u001b[0m\n",
      "\n",
      "\u001b[1m[\u001b[0m\u001b[34;49mnotice\u001b[0m\u001b[1;39;49m]\u001b[0m\u001b[39;49m A new release of pip available: \u001b[0m\u001b[31;49m22.3.1\u001b[0m\u001b[39;49m -> \u001b[0m\u001b[32;49m23.1\u001b[0m\n",
      "\u001b[1m[\u001b[0m\u001b[34;49mnotice\u001b[0m\u001b[1;39;49m]\u001b[0m\u001b[39;49m To update, run: \u001b[0m\u001b[32;49mpip install --upgrade pip\u001b[0m\n"
     ]
    }
   ],
   "source": [
    "!pip install yt > /dev/null  # install yt\n",
    "!pip install pooch pandas > /dev/null # to download example datasets"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 124,
   "metadata": {},
   "outputs": [],
   "source": [
    "import yt\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# First steps - load and inspect metadata"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 343,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:35,797 Parameters: current_time              = 2.0524809024633055\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:35,797 Parameters: domain_dimensions         = [512 512 512]\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:35,798 Parameters: domain_left_edge          = [0. 0. 0.]\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:35,798 Parameters: domain_right_edge         = [1. 1. 1.]\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:35,798 Parameters: cosmological_simulation   = 1\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:35,798 Parameters: current_redshift          = 2.0576423234187393\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:35,799 Parameters: omega_lambda              = 0.68419998884201\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:35,799 Parameters: omega_matter              = 0.31580001115799\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:35,799 Parameters: omega_radiation           = 0.0\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:35,799 Parameters: hubble_constant           = 0.6731999969482421\n"
     ]
    }
   ],
   "source": [
    "# ds = yt.load_sample(\"output_00080\")\n",
    "ds = yt.load(\"output_00230\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 243,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "unyt_quantity(3.18279309, 'Gyr')"
      ]
     },
     "execution_count": 243,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds.current_time.to(\"Gyr\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 244,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2.0576423234187393"
      ]
     },
     "execution_count": 244,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds.current_redshift"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 245,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'HydroMethod': 'ramses',\n",
       " 'Time': 1.0,\n",
       " 'ncpu': 352,\n",
       " 'ndim': 3,\n",
       " 'levelmin': 9,\n",
       " 'levelmax': 22,\n",
       " 'ngridmax': 500000,\n",
       " 'nstep_coarse': 8109,\n",
       " 'boxlen': 1.0,\n",
       " 'time': -2.2198489389408,\n",
       " 'aexp': 0.327049371452284,\n",
       " 'H0': 67.3199996948242,\n",
       " 'omega_m': 0.31580001115799,\n",
       " 'omega_l': 0.68419998884201,\n",
       " 'omega_k': 0.0,\n",
       " 'omega_b': 0.049,\n",
       " 'unit_l': 1.49630430935895e+26,\n",
       " 'unit_d': 7.69162751800601e-29,\n",
       " 'unit_t': 4.89365387559142e+16,\n",
       " 'ordering type': 'hilbert',\n",
       " 'namelist': 'An error occurred when reading the namelist: This functionality requires the f90nml package to be installed.'}"
      ]
     },
     "execution_count": 245,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds.parameters"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The \"raw\" fields are contained in the `io` and `ramses` field types (for particle and grid-based data respectively)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 246,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:42:42,391 Adding particle_type: DM\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:42:42,398 Adding particle_type: star\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:42:42,405 Adding particle_type: cloud\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:42:42,412 Adding particle_type: dust\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:42:42,419 Adding particle_type: star_tracer\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:42:42,426 Adding particle_type: cloud_tracer\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:42:42,434 Adding particle_type: dust_tracer\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:42:42,442 Adding particle_type: gas_tracer\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "On-Disk Field ('ramses', 'Density'): (units: 'code_density')"
      ]
     },
     "execution_count": 246,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds.fields.ramses.Density"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 247,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('io', 'age')\n",
      "('io', 'conformal_birth_time')\n",
      "('io', 'mesh_id')\n",
      "('io', 'particle_angular_momentum')\n",
      "('io', 'particle_angular_momentum_magnitude')\n",
      "('io', 'particle_angular_momentum_x')\n",
      "('io', 'particle_angular_momentum_y')\n",
      "('io', 'particle_angular_momentum_z')\n",
      "('io', 'particle_birth_time')\n",
      "('io', 'particle_family')\n",
      "('io', 'particle_identity')\n",
      "('io', 'particle_index')\n",
      "('io', 'particle_level')\n",
      "('io', 'particle_mass')\n",
      "('io', 'particle_metallicity')\n",
      "('io', 'particle_ones')\n",
      "('io', 'particle_partp')\n",
      "('io', 'particle_position')\n",
      "('io', 'particle_position_cylindrical_radius')\n",
      "('io', 'particle_position_cylindrical_theta')\n",
      "('io', 'particle_position_cylindrical_z')\n",
      "('io', 'particle_position_relative_x')\n",
      "('io', 'particle_position_relative_y')\n",
      "('io', 'particle_position_relative_z')\n",
      "('io', 'particle_position_spherical_phi')\n",
      "('io', 'particle_position_spherical_radius')\n",
      "('io', 'particle_position_spherical_theta')\n",
      "('io', 'particle_position_x')\n",
      "('io', 'particle_position_y')\n",
      "('io', 'particle_position_z')\n",
      "('io', 'particle_radial_velocity')\n",
      "('io', 'particle_radius')\n",
      "('io', 'particle_specific_angular_momentum')\n",
      "('io', 'particle_specific_angular_momentum_x')\n",
      "('io', 'particle_specific_angular_momentum_y')\n",
      "('io', 'particle_specific_angular_momentum_z')\n",
      "('io', 'particle_tag')\n",
      "('io', 'particle_velocity')\n",
      "('io', 'particle_velocity_cylindrical_radius')\n",
      "('io', 'particle_velocity_cylindrical_theta')\n",
      "('io', 'particle_velocity_cylindrical_z')\n",
      "('io', 'particle_velocity_magnitude')\n",
      "('io', 'particle_velocity_relative_x')\n",
      "('io', 'particle_velocity_relative_y')\n",
      "('io', 'particle_velocity_relative_z')\n",
      "('io', 'particle_velocity_spherical_phi')\n",
      "('io', 'particle_velocity_spherical_radius')\n",
      "('io', 'particle_velocity_spherical_theta')\n",
      "('io', 'particle_velocity_x')\n",
      "('io', 'particle_velocity_y')\n",
      "('io', 'particle_velocity_z')\n",
      "('io', 'relative_particle_position')\n",
      "('io', 'relative_particle_position_x')\n",
      "('io', 'relative_particle_position_y')\n",
      "('io', 'relative_particle_position_z')\n",
      "('io', 'relative_particle_velocity')\n",
      "('io', 'relative_particle_velocity_x')\n",
      "('io', 'relative_particle_velocity_y')\n",
      "('io', 'relative_particle_velocity_z')\n",
      "('io', 'star_age')\n",
      "('ramses', 'Density')\n",
      "('ramses', 'Metallicity')\n",
      "('ramses', 'Pressure')\n",
      "('ramses', 'cell_volume')\n",
      "('ramses', 'dx')\n",
      "('ramses', 'dy')\n",
      "('ramses', 'dz')\n",
      "('ramses', 'hydro_dust')\n",
      "('ramses', 'hydro_refinement_scalar')\n",
      "('ramses', 'path_element_x')\n",
      "('ramses', 'path_element_y')\n",
      "('ramses', 'path_element_z')\n",
      "('ramses', 'volume')\n",
      "('ramses', 'x')\n",
      "('ramses', 'x-velocity')\n",
      "('ramses', 'y')\n",
      "('ramses', 'y-velocity')\n",
      "('ramses', 'z')\n",
      "('ramses', 'z-velocity')\n"
     ]
    }
   ],
   "source": [
    "print(\"\\n\".join(list(str(v) for v in ds.derived_field_list if v[0] in (\"ramses\", \"io\"))))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 248,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "level\t       # cells\t     # cells^3\n",
      "--------------------------------------\n",
      "  0\t             0\t             0\n",
      "  1\t             1\t             1\n",
      "  2\t             8\t             2\n",
      "  3\t            64\t             4\n",
      "  4\t           512\t             8\n",
      "  5\t          4096\t            16\n",
      "  6\t         32768\t            32\n",
      "  7\t        262144\t            64\n",
      "  8\t       2097152\t           128\n",
      "  9\t      39205278\t           340\n",
      " 10\t        101708\t            47\n",
      " 11\t        173102\t            56\n",
      " 12\t        336000\t            70\n",
      " 13\t        323362\t            69\n",
      " 14\t        313552\t            68\n",
      " 15\t        295898\t            67\n",
      " 16\t        299876\t            67\n",
      " 17\t        320874\t            69\n",
      " 18\t        171970\t            56\n",
      " 19\t        111658\t            49\n",
      " 20\t         67466\t            41\n",
      " 21\t             0\t             0\n",
      " 22\t             0\t             0\n",
      "----------------------------------------------\n",
      "   \t      44117489\n",
      "\n",
      "\n",
      "z = 2.05764232\n",
      "t = 2.05248090e+00 = 1.00441311e+17 s = 3.18279309e+09 years\n",
      "\n",
      "Smallest Cell:\n",
      "\tWidth: 4.625e-05 Mpc Mpc\n",
      "\tWidth: 4.625e+01 pc pc\n",
      "\tWidth: 9.539e+06 AU AU\n",
      "\tWidth: 1.427e+20 cm cm\n"
     ]
    }
   ],
   "source": [
    "ds.print_stats()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Accessing the raw data\n",
    "\n",
    "In yt, the `ds` object contains the _metadata_. One needs to select a region to access the actual data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 249,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:43:39,324 max value is 9.32094e-22 at 0.4876856803894043 0.5020756721496582 0.5048756599426270\n"
     ]
    }
   ],
   "source": [
    "ad = ds.all_data()  # ad : All Data\n",
    "sp = ds.sphere(\"max\", (1, \"Mpc\"))  # sp : Sphere"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 250,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "((140007792,), (145738770,))"
      ]
     },
     "execution_count": 250,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ad[\"gas\", \"density\"].shape, ad[\"io\", \"particle_identity\"].shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 251,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "[encoded data, 139348 characters]",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "rho, x, y, z = ad[\"gas\", \"density\"], ad[\"ramses\", \"x\"], ad[\"ramses\", \"y\"], ad[\"ramses\", \"z\"]\n",
    "\n",
    "plt.hist2d(x.to(\"Mpc\"), y.to(\"Mpc\"), bins=256, weights=rho, norm=plt.matplotlib.colors.LogNorm())\n",
    "plt.xlabel(\"x [Mpc]\")\n",
    "plt.ylabel(\"y [Mpc]\")\n",
    "plt.gca().set_aspect(\"equal\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 252,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "[encoded data, 310332 characters]",
      "text/plain": [
       "<Figure size 1200x600 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "rho, x, y, z = sp[\"gas\", \"density\"], sp[\"ramses\", \"x\"], sp[\"ramses\", \"y\"], sp[\"ramses\", \"z\"]\n",
    "xp = sp[\"io\", \"particle_position\"]\n",
    "fig, axes = plt.subplots(1, 2, figsize=(12, 6), sharex=True, sharey=True)\n",
    "\n",
    "axes[0].hist2d(y.to(\"Mpc\"), z.to(\"Mpc\"), bins=256, weights=rho, norm=plt.matplotlib.colors.LogNorm())\n",
    "axes[1].scatter(xp[:, 1].to(\"Mpc\"), xp[:, 2].to(\"Mpc\"), s=0.1, c=\"gray\", alpha=0.5)\n",
    "axes[0].set_xlabel(\"x [Mpc]\")\n",
    "axes[1].set_xlabel(\"x [Mpc]\")\n",
    "axes[0].set_ylabel(\"y [Mpc]\")\n",
    "axes[0].set_aspect(\"equal\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Plotting"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 253,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:46:13,456 xlim = 0.002076 1.002076\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:46:13,457 ylim = 0.004876 1.004876\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:46:13,458 xlim = 0.002076 1.002076\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:46:13,458 ylim = 0.004876 1.004876\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:46:13,463 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<img style=\"max-width:100%;max-height:100%;\" src=\"data:image/png;base64,[encoded data, 1111040 characters]\"><br>"
      ],
      "text/plain": [
       "<yt.visualization.plot_window.AxisAlignedSlicePlot at 0x7f5ba5f2d5d0>"
      ]
     },
     "execution_count": 253,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "p = yt.SlicePlot(ds, \"x\", (\"gas\", \"density\"), center=\"max\")\n",
    "p"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 268,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:52:01,912 xlim = 0.496920 0.507231\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:52:01,914 ylim = 0.499720 0.510031\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:52:01,915 xlim = 0.496920 0.507231\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:52:01,916 ylim = 0.499720 0.510031\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 13:52:01,917 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<img style=\"max-width:100%;max-height:100%;\" src=\"data:image/png;base64,[encoded data, 118464 characters]\"><br>"
      ],
      "text/plain": [
       "<yt.visualization.plot_window.AxisAlignedSlicePlot at 0x7f5b98dbebf0>"
      ]
     },
     "execution_count": 268,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "p = yt.SlicePlot(ds, \"x\", (\"gas\", \"density\"), center=sp.center, data_source=sp, width=(500, \"kpc\"))\n",
    "p"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 307,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:03:51,957 Projection completed\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:03:51,960 xlim = 0.487170 0.488201\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:03:51,961 ylim = 0.501560 0.502591\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:03:51,963 xlim = 0.487170 0.488201\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:03:51,963 ylim = 0.501560 0.502591\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:03:51,964 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:03:52,065 Making a fixed resolution buffer of (('gas', 'temperature')) 800 by 800\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<img style=\"max-width:100%;max-height:100%;\" src=\"data:image/png;base64,[encoded data, 153648 characters]\"><br><img style=\"max-width:100%;max-height:100%;\" src=\"data:image/png;base64,[encoded data, 119140 characters]\"><br>"
      ],
      "text/plain": [
       "<yt.visualization.plot_window.AxisAlignedProjectionPlot at 0x7f5b97280760>"
      ]
     },
     "execution_count": 307,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "p = yt.ProjectionPlot(ds, \"z\", [(\"gas\", \"density\"), (\"gas\", \"temperature\")], center=sp.center, data_source=sp, width=(50, \"kpc\"), weight_field=(\"gas\", \"density\"))\n",
    "p.set_cmap(\"density\", \"viridis\")\n",
    "p.set_cmap(\"temperature\", \"bwr\")\n",
    "p.set_zlim(\"temperature\", 1e3, 1e7)\n",
    "p.set_unit((\"gas\", \"density\"), \"mp/cm**3\")\n",
    "p"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Advanced topic: volume rendering"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 306,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([-0.027016  , -0.26685463,  0.96335806])"
      ]
     },
     "execution_count": 306,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "spp = ds.sphere(sp.center, (20, \"kpc\"))\n",
    "AM = spp.quantities.angular_momentum_vector(use_gas=True, use_particles=False)\n",
    "AM_unit = AM.value / np.linalg.norm(AM)\n",
    "\n",
    "AM_unit"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "ename": "AttributeError",
     "evalue": "'TransferFunctionHelper' object has no attribute 'set'",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mAttributeError\u001b[0m                            Traceback (most recent call last)",
      "Cell \u001b[0;32mIn[310], line 15\u001b[0m\n\u001b[1;32m     13\u001b[0m rho_min, rho_max \u001b[39m=\u001b[39m ds\u001b[39m.\u001b[39marr([\u001b[39m1e-3\u001b[39m, \u001b[39m1e2\u001b[39m], \u001b[39m\"\u001b[39m\u001b[39mmp/cm**3\u001b[39m\u001b[39m\"\u001b[39m)\u001b[39m.\u001b[39mto(\u001b[39m\"\u001b[39m\u001b[39mg/cm**3\u001b[39m\u001b[39m\"\u001b[39m)\n\u001b[1;32m     14\u001b[0m source\u001b[39m.\u001b[39mtfh\u001b[39m.\u001b[39mset_bounds((rho_min, rho_max))\n\u001b[0;32m---> 15\u001b[0m source\u001b[39m.\u001b[39;49mtfh\u001b[39m.\u001b[39;49mset\n\u001b[1;32m     16\u001b[0m source\u001b[39m.\u001b[39mtfh\u001b[39m.\u001b[39mgrey_opacity \u001b[39m=\u001b[39m \u001b[39mTrue\u001b[39;00m\n\u001b[1;32m     17\u001b[0m source\u001b[39m.\u001b[39mtfh\u001b[39m.\u001b[39mplot(\u001b[39m'\u001b[39m\u001b[39mtransfer_function.png\u001b[39m\u001b[39m'\u001b[39m, profile_field\u001b[39m=\u001b[39m\u001b[39m'\u001b[39m\u001b[39mdensity\u001b[39m\u001b[39m'\u001b[39m)\n",
      "\u001b[0;31mAttributeError\u001b[0m: 'TransferFunctionHelper' object has no attribute 'set'"
     ]
    }
   ],
   "source": [
    "# data_source = ds.all_data()\n",
    "# data_source = ds.sphere(\"max\", (1, \"Mpc\"))\n",
    "data_source = ds.sphere(sp.center, (200, \"kpc\"))\n",
    "\n",
    "sc = yt.create_scene(data_source, lens_type='perspective', field=(\"gas\", \"density\"))\n",
    "\n",
    "sc.camera.resolution = (1024, 1024)\n",
    "sc.camera.focus = data_source.argmax(\"density\")\n",
    "sc.camera.position = sc.camera.focus + ds.quan(200, \"kpc\") * AM_unit\n",
    "\n",
    "source = sc[0]\n",
    "\n",
    "rho_min, rho_max = ds.arr([1e-3, 1e2], \"mp/cm**3\").to(\"g/cm**3\")\n",
    "source.tfh.set_bounds((rho_min, rho_max))\n",
    "source.tfh.grey_opacity = True\n",
    "source.tfh.plot('transfer_function.png', profile_field='density')\n",
    "\n",
    "source.tfh.plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:02:54,236 Rendering scene (Can take a while).\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:02:54,244 Creating volume\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:55,425 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 39, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:55,446 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 44, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:55,462 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 46, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:55,558 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 47, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:55,600 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 65, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:55,637 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 66, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:55,659 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 67, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:55,734 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 68, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:56,147 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 69, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:56,180 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 70, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:56,232 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 71, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:56,394 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 72, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:56,683 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 73, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:56,741 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 74, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:56,866 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 75, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:56,936 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 76, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:56,964 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 77, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,031 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 78, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,055 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 79, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,080 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 80, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,100 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 81, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,116 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 82, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,130 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 83, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,150 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 84, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,201 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 85, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,229 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 86, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,311 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 87, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,328 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 88, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,383 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 89, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,425 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 90, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,455 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 91, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,470 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 92, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,483 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 93, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,504 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 94, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,530 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 95, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,552 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 96, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,577 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 97, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,647 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 98, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,756 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 99, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,814 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 100, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,869 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 101, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,962 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 102, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:57,983 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 103, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,106 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 104, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,289 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 105, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,330 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 106, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,381 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 107, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,481 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 108, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,502 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 109, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,549 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 110, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,579 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 111, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,698 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 112, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,835 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 113, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,866 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 114, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,966 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 115, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:58,989 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 116, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,061 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 117, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,139 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 118, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,220 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 119, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,280 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 120, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,342 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 121, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,356 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 122, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,386 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 123, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,406 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 124, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,419 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 125, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,432 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 126, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,446 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 127, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,465 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 128, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,482 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 129, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,506 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 130, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,534 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 131, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,553 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 132, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,571 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 133, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,589 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 134, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,604 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 135, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,616 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 136, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,630 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 137, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,646 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 138, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,662 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 139, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,678 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 140, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,692 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 141, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,704 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 146, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,715 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 147, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,726 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 148, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,738 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 149, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,755 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 150, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,770 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 153, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,783 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 157, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,795 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 158, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,812 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 159, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,829 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 160, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,846 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 161, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,863 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 162, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,886 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 163, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,932 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 164, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:02:59,978 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 165, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,015 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 166, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,056 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 167, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,101 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 168, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,194 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 169, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,217 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 170, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,272 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 171, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,307 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 172, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,334 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 173, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,376 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 174, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,403 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 175, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,427 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 176, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,443 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 177, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,459 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 178, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,472 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 179, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,488 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 180, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,503 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 182, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,514 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 183, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,527 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 184, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,542 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 185, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,553 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 186, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,565 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 187, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,580 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 188, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,593 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 189, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,606 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 190, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,621 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 191, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,640 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 192, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,661 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 193, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,684 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 194, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,717 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 195, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,740 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 196, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,762 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 197, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,777 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 198, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,795 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 199, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,825 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 200, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:00,913 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 201, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:01,035 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 202, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:01,202 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 203, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:01,325 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 204, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:01,455 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 205, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:01,560 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 206, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:01,682 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 207, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:01,780 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 208, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:01,848 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 209, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:01,880 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 210, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:02,082 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 211, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:02,120 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 212, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:02,190 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 213, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:02,268 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 214, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:02,476 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 215, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:02,815 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 229, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:02,823 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 231, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:02,837 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 232, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:02,860 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 233, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:03:02,905 RAMSESDomainSubset (info_00230): , base_region=YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm, domain=RAMSESDomainFile: 234, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n"
     ]
    },
    {
     "data": {
      "image/png": "[encoded data, 404556 characters]",
      "text/plain": [
       "<Scene Object>:\n",
       "Sources: \n",
       "    source_00: <Volume Source>:YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm transfer_function:<Color Transfer Function Object>:\n",
       "x_bounds:[-27, -23] nbins:512 features:\n",
       "\t('gaussian', 'location(x):-27', 'width(x):0.002', 'height(y):(0.073,   0, 0.084,   1)')\n",
       "\t('gaussian', 'location(x):-26', 'width(x):0.002', 'height(y):(0.35,   0, 0.62,   1)')\n",
       "\t('gaussian', 'location(x):-26', 'width(x):0.002', 'height(y):(  0, 0.26, 0.87,   1)')\n",
       "\t('gaussian', 'location(x):-25', 'width(x):0.002', 'height(y):(  0, 0.65, 0.72,   1)')\n",
       "\t('gaussian', 'location(x):-25', 'width(x):0.002', 'height(y):(  0, 0.6, 0.031,   1)')\n",
       "\t('gaussian', 'location(x):-25', 'width(x):0.002', 'height(y):(  0, 0.87,   0,   1)')\n",
       "\t('gaussian', 'location(x):-24', 'width(x):0.002', 'height(y):(0.78, 0.98,   0,   1)')\n",
       "\t('gaussian', 'location(x):-24', 'width(x):0.002', 'height(y):(  1, 0.71,   0,   1)')\n",
       "\t('gaussian', 'location(x):-23', 'width(x):0.002', 'height(y):(0.91,   0,   0,   1)')\n",
       "\t('gaussian', 'location(x):-23', 'width(x):0.002', 'height(y):(0.8, 0.67, 0.67,   1)')\n",
       "\n",
       "Camera: \n",
       "    <Camera Object>:\n",
       "\tposition:[0.48757426 0.50097506 0.50884893] unitary\n",
       "\tfocus:[0.48768568 0.50207567 0.50487566] code_length\n",
       "\tnorth_vector:[ 0.84143575 -0.5263568  -0.1222064 ] dimensionless\n",
       "\twidth:[1.5 1.5 1.5] code_length\n",
       "\tlight:None\n",
       "\tresolution:(1024, 1024)\n",
       "Lens: <Lens Object>:\n",
       "\tlens_type:perspective\n",
       "\tviewpoint:[ 0.50794768  0.70221664 -0.21764288] code_length"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sc.render()\n",
    "sc.show(sigma_clip=8.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "[encoded data, 508020 characters]",
      "text/plain": [
       "<Scene Object>:\n",
       "Sources: \n",
       "    source_00: <Volume Source>:YTSphere (info_00230): , center=[7.29726185e+25 7.51257992e+25 7.55447626e+25] cm, radius=6.17135516192465e+23 cm transfer_function:<Color Transfer Function Object>:\n",
       "x_bounds:[-27, -23] nbins:512 features:\n",
       "\t('gaussian', 'location(x):-27', 'width(x):0.002', 'height(y):(0.073,   0, 0.084,   1)')\n",
       "\t('gaussian', 'location(x):-26', 'width(x):0.002', 'height(y):(0.35,   0, 0.62,   1)')\n",
       "\t('gaussian', 'location(x):-26', 'width(x):0.002', 'height(y):(  0, 0.26, 0.87,   1)')\n",
       "\t('gaussian', 'location(x):-25', 'width(x):0.002', 'height(y):(  0, 0.65, 0.72,   1)')\n",
       "\t('gaussian', 'location(x):-25', 'width(x):0.002', 'height(y):(  0, 0.6, 0.031,   1)')\n",
       "\t('gaussian', 'location(x):-25', 'width(x):0.002', 'height(y):(  0, 0.87,   0,   1)')\n",
       "\t('gaussian', 'location(x):-24', 'width(x):0.002', 'height(y):(0.78, 0.98,   0,   1)')\n",
       "\t('gaussian', 'location(x):-24', 'width(x):0.002', 'height(y):(  1, 0.71,   0,   1)')\n",
       "\t('gaussian', 'location(x):-23', 'width(x):0.002', 'height(y):(0.91,   0,   0,   1)')\n",
       "\t('gaussian', 'location(x):-23', 'width(x):0.002', 'height(y):(0.8, 0.67, 0.67,   1)')\n",
       "\n",
       "Camera: \n",
       "    <Camera Object>:\n",
       "\tposition:[0.48757426 0.50097506 0.50884893] unitary\n",
       "\tfocus:[0.48768568 0.50207567 0.50487566] code_length\n",
       "\tnorth_vector:[ 0.84143575 -0.5263568  -0.1222064 ] dimensionless\n",
       "\twidth:[1.5 1.5 1.5] code_length\n",
       "\tlight:None\n",
       "\tresolution:(1024, 1024)\n",
       "Lens: <Lens Object>:\n",
       "\tlens_type:perspective\n",
       "\tviewpoint:[ 0.50794768  0.70221664 -0.21764288] code_length"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sc.show(sigma_clip=3)"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Advanced topic: particle-to-mesh, mesh-to-particle"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Particle-to-mesh\n",
    "\n",
    "This implements nearest/CIC interpolation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 377,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:26:30,142 Projection completed\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:26:30,144 xlim = 0.485623 0.489748\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:26:30,144 ylim = 0.500013 0.504138\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:26:30,145 xlim = 0.485623 0.489748\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:26:30,146 ylim = 0.500013 0.504138\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:26:30,147 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:26:30,261 Making a fixed resolution buffer of (('deposit', 'gas_tracer_cic')) 800 by 800\n"
     ]
    }
   ],
   "source": [
    "p = yt.ProjectionPlot(ds, \"z\",\n",
    "    [\n",
    "        (\"deposit\", \"gas_tracer_cic\"),\n",
    "        (\"gas\", \"density\"),\n",
    "    ],\n",
    "    data_source=sp,\n",
    "    width=(200, \"kpc\"),\n",
    "    center=sp.center,\n",
    "    weight_field=\"density\"\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 381,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img style=\"max-width:100%;max-height:100%;\" src=\"data:image/png;base64,[encoded data, 155412 characters]\"><br><img style=\"max-width:100%;max-height:100%;\" src=\"data:image/png;base64,[encoded data, 199780 characters]\"><br>"
      ],
      "text/plain": [
       "<yt.visualization.plot_window.AxisAlignedProjectionPlot at 0x7f5b9b0ee500>"
      ]
     },
     "execution_count": 381,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "p.set_unit((\"gas\", \"density\"), \"mp/cm**3\")\n",
    "p.set_unit((\"deposit\", \"gas_tracer_cic\"), \"mp/cm**3\")\n",
    "p.set_zlim((\"gas\", \"density\"), 1e-3, 1e2)\n",
    "p.set_zlim((\"deposit\", \"gas_tracer_cic\"), 1e-3, 1e2)\n",
    "p.set_cmap((\"gas\", \"density\"), \"dusk\")\n",
    "p.set_cmap((\"deposit\", \"gas_tracer_cic\"), \"dusk\")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Mesh-to-particle"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 383,
   "metadata": {},
   "outputs": [],
   "source": [
    "ds.add_mesh_sampling_particle_field((\"gas\", \"temperature\"), ptype=\"gas_tracer\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 404,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "There are 513,120 particles\n",
      "\t⇒ sp[\"gas_tracer\", \"cell_gas_temperature\"].shape = (513120,)\n",
      "There are 5,495,675 gas cells in the region\n"
     ]
    }
   ],
   "source": [
    "ids = sp[\"gas_tracer\", \"particle_identity\"].astype(int)\n",
    "T = sp[\"gas_tracer\", \"cell_gas_temperature\"]\n",
    "\n",
    "print(f'There are {len(ids):,} particles')\n",
    "print(f'\\t⇒ sp[\"gas_tracer\", \"cell_gas_temperature\"].shape = ({T.shape[0]:},)')\n",
    "print(f'There are {len(sp[\"gas\", \"density\"]):,} gas cells in the region')"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "A more complex example. This will assign the cell velocity divergence to the gas tracer particle in the current cell.\n",
    "Internally, this will:\n",
    "1. For each cell, find the 3×3×3 neighboring cells\n",
    "2. Compute the divergence of the velocity field with a finite-difference \n",
    "3. Find which cell contains each particle\n",
    "4. Assign the velocity divergence to the particle"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 405,
   "metadata": {},
   "outputs": [],
   "source": [
    "ds.add_mesh_sampling_particle_field((\"gas\", \"velocity_divergence\"), ptype=\"gas_tracer\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "sp[\"gas_traer\"] "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 433,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "[encoded data, 43600 characters]",
      "text/plain": [
       "<Figure size 1200x600 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, axes = plt.subplots(1, 2, figsize=(12, 6), sharey=True, sharex=True, constrained_layout=True)\n",
    "\n",
    "for sign, ax, title, cmap in zip((-1, 1), axes, (\"$\\mathrm{div}\\ \\mathbf{v}<0$\", \"$\\mathrm{div}\\ \\mathbf{v}>0$\"), (\"Blues\", \"Reds\")):\n",
    "    mask = sp[\"gas_tracer\", \"cell_gas_velocity_divergence\"] * sign > 0\n",
    "    ax.hist2d(\n",
    "        np.abs((sp[\"gas_tracer\", \"cell_gas_velocity_divergence\"].to(\"Gyr**-1\"))[mask]),\n",
    "        sp[\"gas_tracer\", \"cell_gas_temperature\"][mask],\n",
    "        bins=(np.geomspace(1e-6, 1e3, 50), np.geomspace(1e2, 1e7, 50)),\n",
    "        norm=plt.matplotlib.colors.LogNorm(),\n",
    "        cmap=cmap,\n",
    "    )\n",
    "    ax.set(\n",
    "        xlabel=r\"$|\\mathrm{div}\\ \\mathbf{v}|$ ($\\mathrm{Gyr}^{-1}$)\",\n",
    "        ylabel=r\"$T$ ($\\mathrm{K}$)\",\n",
    "        xscale=\"log\",\n",
    "        yscale=\"log\",\n",
    "        title=title,\n",
    "    )"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Advanced topic: adding derived fields"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Local-cell access"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 344,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:12:58,237 The DM particle filter already exists. Overriding.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:12:58,237 The star particle filter already exists. Overriding.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:12:58,237 The cloud particle filter already exists. Overriding.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:12:58,238 The dust particle filter already exists. Overriding.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:12:58,238 The star_tracer particle filter already exists. Overriding.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:12:58,238 The cloud_tracer particle filter already exists. Overriding.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:12:58,239 The dust_tracer particle filter already exists. Overriding.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:12:58,239 The gas_tracer particle filter already exists. Overriding.\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:58,239 Adding particle_type: DM\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:58,247 Adding particle_type: star\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:58,255 Adding particle_type: cloud\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:58,264 Adding particle_type: dust\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:58,273 Adding particle_type: star_tracer\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:58,281 Adding particle_type: cloud_tracer\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:58,290 Adding particle_type: dust_tracer\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:12:58,299 Adding particle_type: gas_tracer\n"
     ]
    }
   ],
   "source": [
    "ad = ds.all_data()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 352,
   "metadata": {},
   "outputs": [],
   "source": [
    "import yt.units as u\n",
    "\n",
    "def _entropy(field, data):\n",
    "    K = u.kb * data[\"gas\", \"temperature\"] * data[\"gas\", \"Electron_number_density\"]**(-2/3)\n",
    "    return K\n",
    "\n",
    "ds.add_field(name=(\"gas\", \"entropy\"), units=\"cm**4 * g / s**2\", function=_entropy, sampling_type=\"cell\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 359,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:21:10,617 Projection completed\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:21:10,618 xlim = 0.496920 0.507231\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:21:10,618 ylim = 0.499720 0.510031\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:21:10,618 xlim = 0.496920 0.507231\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:21:10,619 ylim = 0.499720 0.510031\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:21:10,619 Making a fixed resolution buffer of (('gas', 'entropy')) 800 by 800\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<img style=\"max-width:100%;max-height:100%;\" src=\"data:image/png;base64,[encoded data, 174384 characters]\"><br>"
      ],
      "text/plain": [
       "<yt.visualization.plot_window.AxisAlignedProjectionPlot at 0x7f5b9bb0b130>"
      ]
     },
     "execution_count": 359,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "p = yt.ProjectionPlot(ds, \"x\", (\"gas\", \"entropy\"), center=sp.center, data_source=sp, width=(500, \"kpc\"), weight_field=(\"gas\", \"density\"))\n",
    "p.set_cmap(\"entropy\", \"magma\")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Spatially-aware quantities\n",
    "\n",
    "**Warning**: yt doesn't know which variable is intensive/extensive! So the interpolation may be wrong depending on the variable."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 434,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:46:06,754 max value is 9.32094e-22 at 0.4876856803894043 0.5020756721496582 0.5048756599426270\n"
     ]
    }
   ],
   "source": [
    "def _velocity_dispersion(field, data):\n",
    "    from itertools import product\n",
    "    sigma = np.zeros(data[\"gas\", \"momentum_x\"].shape)\n",
    "\n",
    "    print(data[\"gas\", \"density\"].shape)\n",
    "\n",
    "    for direction in \"xyz\":\n",
    "        v = data[\"gas\", f\"momentum_{direction}\"].value\n",
    "        w = np.ones_like(data[\"cell_mass\"].value)\n",
    "        for i, j, k in product((1, 2), (1, 2), (1, 2)):\n",
    "            vv = v[i-1:i+2, j-1:j+2, k-1:k+2].copy()\n",
    "            ww = w[i-1:i+2, j-1:j+2, k-1:k+2].copy()\n",
    "            vv_mean = np.average(vv, weights=ww, axis=(0, 1, 2))\n",
    "            sigma[i, j, k] += np.average((vv-vv_mean)**2, weights=ww, axis=(0, 1, 2))\n",
    "\n",
    "    return data.apply_units(np.sqrt(sigma), data[\"gas\", \"momentum_x\"].units) / data[\"gas\", \"cell_mass\"] / 27\n",
    "\n",
    "ds.add_field(name=(\"gas\", \"velocity_dispersion\"), units=u.km/u.s, function=_velocity_dispersion,\n",
    "             sampling_type=\"cell\", force_override=True, validators=[yt.ValidateSpatial(ghost_zones=1)])\n",
    "ad = ds.all_data()\n",
    "sp = ds.sphere(\"max\", (1, \"Mpc\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 435,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:29,882 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 14, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:31,287 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 16, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:31,882 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 39, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:36,612 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 40, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:37,702 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 43, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:39,974 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 44, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:41,408 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 45, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:42,198 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 46, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:45,403 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 47, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:50,174 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 55, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:51,205 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 60, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:53,300 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 64, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:54,212 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 65, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:55,961 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 72, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:57,862 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 73, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:58,056 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 74, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:58,532 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 75, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:58,776 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 76, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:58,861 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 77, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:59,053 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 79, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:59,130 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 80, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:59,178 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 81, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:59,220 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 82, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:59,260 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 83, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:59,317 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 84, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:59,476 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 86, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:59,769 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 87, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:59,829 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 88, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:46:59,998 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 89, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:00,110 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 90, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:00,189 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 92, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:00,222 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 93, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:00,270 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 94, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:00,336 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 100, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:00,490 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 101, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:00,747 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 103, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:01,056 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 104, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:01,504 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 112, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:02,082 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 114, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:02,446 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 116, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:02,663 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 117, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:02,887 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 130, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:02,965 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 131, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,013 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 132, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,056 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 133, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,099 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 161, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,141 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 162, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,201 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 163, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,311 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 164, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,430 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 168, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,693 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 173, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,819 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 174, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,891 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 175, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,946 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 176, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:03,989 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 177, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:04,033 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 189, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:04,071 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 191, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:04,121 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 192, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:04,166 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 193, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:04,216 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 194, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:04,286 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 200, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:04,541 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 202, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:05,165 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 203, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:05,622 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 205, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:06,016 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 206, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:06,459 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 207, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:06,824 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 215, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:12,100 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 217, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:15,102 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 218, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:16,043 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 220, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:16,431 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 221, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:17,357 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 222, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:24,934 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 263, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[33mWARNING\u001b[0m  ] 2023-04-17 14:47:25,335 RAMSESDomainSubset (info_00230): , base_region=YTSlice (info_00230): , axis=0, coord=7.297261851793173e+25 cm, domain=RAMSESDomainFile: 264, ds=info_00230.retrieve_ghost_zones was called with the `smoothed` argument set to True. This is not supported, ignoring it.\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:47:27,453 xlim = 0.481454 0.522698\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:47:27,454 ylim = 0.484254 0.525498\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:47:27,456 xlim = 0.481454 0.522698\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:47:27,457 ylim = 0.484254 0.525498\n",
      "yt : [\u001b[32mINFO\u001b[0m     ] 2023-04-17 14:47:27,459 Making a fixed resolution buffer of (('gas', 'velocity_dispersion')) 800 by 800\n",
      "/home/cphyc/Documents/22-23/Conferences/RUM/venv/lib/python3.10/site-packages/yt/visualization/plot_modifications.py:108: VisibleDeprecationWarning: The `ncont` keyword argument is deprecated, use `levels` instead.\n",
      "Deprecated since yt 4.1.0\n",
      "This feature is planned for removal two minor releases later (anticipated)\n",
      "  self._callbacks.append(cls(*args, **kwargs))\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<img style=\"max-width:100%;max-height:100%;\" src=\"data:image/png;base64,[encoded data, 386728 characters]\"><br>"
      ],
      "text/plain": [
       "<yt.visualization.plot_window.AxisAlignedSlicePlot at 0x7f5b6949fc70>"
      ]
     },
     "execution_count": 435,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "s = yt.SlicePlot(ds, \"x\", (\"gas\", \"velocity_dispersion\"), center=\"max\", data_source=sp, width=(2, \"Mpc\"))\n",
    "s.annotate_contour((\"gas\", \"density\"), ncont=10, clim=(1e-30, 1e-25), plot_args={\"alpha\": 0.5})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "venv",
   "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.10.10"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
