sdf-xarray 0.7.4__tar.gz → 0.9.0__tar.gz

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/.github/workflows/build_publish.yml

@@ -20,7 +20,7 @@ jobs:
         submodules: "recursive"
 
     - name: Cache Zenodo datasets
-      uses: actions/cache@v4
+      uses: actions/cache@v6
       with:
         enableCrossOsArchive: true
         path: ~/.cache/sdf_datasets

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/.github/workflows/tests.yml

@@ -31,7 +31,7 @@ jobs:
         uv sync --python ${{ matrix.python-version }} --frozen
 
     - name: Cache Zenodo datasets
-      uses: actions/cache@v4
+      uses: actions/cache@v6
       with:
         path: ~/.cache/sdf_datasets
         key: sdf-datasets-17991042

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: sdf-xarray
-Version: 0.7.4
+Version: 0.9.0
 Summary: Provides a backend for xarray to read SDF files as created by the EPOCH plasma PIC code.
 Author-Email: Peter Hill <peter.hill@york.ac.uk>, Joel Adams <joel.adams@york.ac.uk>, Shaun Doherty <shaun.doherty@york.ac.uk>, Chris Herdman <chris.herdman@york.ac.uk>, Liam Pattinson <liam.pattinson@york.ac.uk>, Sviatoslav Shekhanov <sviatoslav.shekhanov@york.ac.uk>
 License-Expression: BSD-3-Clause

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/animation.md

@@ -35,9 +35,10 @@ resolved data has 2 dimensions.
 
 | Dimensions |      Plotting function        |        Notes          |
 | ---------- | ----------------------------- | --------------------- |
-|    `2`     | <inv:#xarray.plot.line>       |                       |
-|    `3`     | <inv:#xarray.plot.pcolormesh> |                       |
-|    `>3`    | <inv:#xarray.plot.hist>       | Not fully implemented |
+|     2      | <inv:#xarray.plot.line>       |                       |
+|     3      | <inv:#xarray.plot.pcolormesh> |                       |
+|     4      | <inv:#mpl_toolkits.mplot3d.axes3d.Axes3D.voxels>| WIP |
+|     >4     | <inv:#xarray.plot.hist>       | Not fully implemented |
 
 
 ### 1D simulation
@@ -103,6 +104,10 @@ anim.show()
 
 ### 3D simulation
 
+```{warning}
+[Voxel](./plotting.md#voxel-plots) animations are not currently implemented. 
+```
+
 Opening a 3D simulation as a multi-file dataset and plotting it will
 return a <inv:#xarray.plot.hist>. However, this may not be
 desirable. We can plot a 3D simulation along a certain plane in the
@@ -170,8 +175,8 @@ Some functionality such as `aspect` and `size` are not fully implemented yet.
 ds = sdfxr.open_mfdataset("tutorial_dataset_2d/*.sdf")
 
 # Change the units of the coordinates
-ds = ds.epoch.rescale_coords(1e6, "µm", ["X_Grid_mid", "Y_Grid_mid"])
-ds = ds.epoch.rescale_coords(1e15, "fs", ["time"])
+ds = ds.epoch.rescale_coords("µm", ["X_Grid_mid", "Y_Grid_mid"])
+ds = ds.epoch.rescale_coords("fs", ["time"])
 ds["time"].attrs["long_name"] = "t"
 
 # Change units and name of the variable

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/animating/animate_1D_density.py

@@ -6,9 +6,9 @@ input_dir = Path("datasets/1_1_drifting_bunch")
 ds = sdfxr.open_mfdataset(input_dir)
 
 # Convert the time to femtoseconds
-ds = ds.epoch.rescale_coords(1e15, "fs", "time")
+ds = ds.epoch.rescale_coords("fs", "time", 1e15)
 # Convert the x and y coords to microns
-ds = ds.epoch.rescale_coords(1e6, "µm", ["X_Grid_mid"])
+ds = ds.epoch.rescale_coords("µm", ["X_Grid_mid"], 1e6)
 
 anim = ds["Derived_Number_Density"].epoch.animate()
 

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/animating/animate_1D_poynting_flux.py

@@ -8,9 +8,9 @@ input_dir = Path("datasets/3_3_Gaussian_1d_laser")
 ds = sdfxr.open_mfdataset(input_dir)
 
 # Convert the time to femtoseconds
-ds = ds.epoch.rescale_coords(1e15, "fs", "time")
+ds = ds.epoch.rescale_coords("fs", "time", 1e15)
 # Convert the x and y coords to microns
-ds = ds.epoch.rescale_coords(1e6, "µm", ["X_Grid_mid"])
+ds = ds.epoch.rescale_coords("µm", ["X_Grid_mid"], 1e6)
 
 # Calculate Poynting flux magnitude
 flux_magnitude = np.sqrt(

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/animating/animate_2D_dist_fn_multispecies.py

@@ -8,7 +8,7 @@ ds = sdfxr.open_mfdataset(
 )
 
 # Rescale coords to account for kilometers
-ds = ds.epoch.rescale_coords(1e-3, "km", ["X_x_px_Left"])
+ds = ds.epoch.rescale_coords("km", ["X_x_px_Left"], 1e-3)
 
 # Sum phase-space of species "Left" and "Right" in "x_px" distribution function
 # NOTE: We only use the values from the right distribution function as if we inherit

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/animating/animate_2D_dist_fn_multispecies_alternative.py

@@ -11,7 +11,7 @@ ds = sdfxr.open_mfdataset(
 )
 
 # Rescale coords to account for kilometers
-ds = ds.epoch.rescale_coords(1e-3, "km", ["X_x_px_Left"])
+ds = ds.epoch.rescale_coords("km", ["X_x_px_Left"], 1e-3)
 
 # Sum phase-space of species "Left" and "Right" in "x_px" distribution function
 # NOTE: We only use the values from the right distribution function as if we inherit

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/animating/animate_2D_poynting_flux.py

@@ -8,9 +8,9 @@ input_dir = Path("datasets/3_5_Gaussian_beam")
 ds = sdfxr.open_mfdataset(input_dir)
 
 # Convert the time to femtoseconds
-ds = ds.epoch.rescale_coords(1e15, "fs", "time")
+ds = ds.epoch.rescale_coords("fs", "time", 1e15)
 # Convert the x and y coords to microns
-ds = ds.epoch.rescale_coords(1e6, "µm", ["X_Grid_mid", "Y_Grid_mid"])
+ds = ds.epoch.rescale_coords("µm", ["X_Grid_mid", "Y_Grid_mid"], 1e6)
 
 # Calculate Poynting flux magnitude
 flux_magnitude = np.sqrt(

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/live_demo.ipynb

@@ -3800,7 +3800,7 @@
     },
     {
       "cell_type": "code",
-      "execution_count": 10,
+      "execution_count": null,
       "id": "2ef27498",
       "metadata": {
         "colab": {
@@ -3833,7 +3833,7 @@
         }
       ],
       "source": [
-        "ds = ds.epoch.rescale_coords(1e6, \"µm\", [\"X_Grid_mid\", \"Y_Grid_mid\"])\n",
+        "ds = ds.epoch.rescale_coords(\"µm\", [\"X_Grid_mid\", \"Y_Grid_mid\"], 1e6)\n",
         "ds[\"Derived_Number_Density_Electron\"].isel(time=0).epoch.plot()"
       ]
     },

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/plotting/plot_1D_density.py

@@ -8,7 +8,7 @@ import sdf_xarray as sdfxr
 input_dir = Path("datasets/1_1_drifting_bunch")
 
 ds = sdfxr.open_dataset(input_dir / "0000.sdf")
-ds = ds.epoch.rescale_coords(1e6, "µm", ["X_Grid_mid"])
+ds = ds.epoch.rescale_coords("µm", ["X_Grid_mid"], 1e6)
 
 ds["Derived_Number_Density"].epoch.plot()
 plt.tight_layout()

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/plotting/plot_1D_poynting_flux.py

@@ -10,7 +10,7 @@ input_dir = Path("datasets/3_3_Gaussian_1d_laser")
 ds = sdfxr.open_dataset(input_dir / "0020.sdf")
 
 # Convert the x and y coords to microns
-ds = ds.epoch.rescale_coords(1e6, "µm", ["X_Grid_mid"])
+ds = ds.epoch.rescale_coords("µm", ["X_Grid_mid"], 1e6)
 
 # Calculate Poynting flux magnitude
 flux_magnitude = np.sqrt(

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/plotting/plot_2D_density.py

@@ -8,7 +8,7 @@ import sdf_xarray as sdfxr
 input_dir = Path("datasets/4_3_basic_target")
 
 ds = sdfxr.open_dataset(input_dir / "0000.sdf")
-ds = ds.epoch.rescale_coords(1e6, "µm", ["X_Grid_mid", "Y_Grid_mid"])
+ds = ds.epoch.rescale_coords("µm", ["X_Grid_mid", "Y_Grid_mid"], 1e6)
 
 ds["Derived_Number_Density"].epoch.plot()
 plt.tight_layout()

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/plotting/plot_2D_dist_fn_multispecies.py

@@ -8,7 +8,7 @@ import sdf_xarray as sdfxr
 input_dir = Path("datasets/2_1_two_stream_instability")
 ds = sdfxr.open_dataset(input_dir / "0000.sdf")
 
-ds = ds.epoch.rescale_coords(1e-3, "km", ["X_x_px_Left"])
+ds = ds.epoch.rescale_coords("km", ["X_x_px_Left"], 1e-3)
 
 # Sum phase-space of species "Left" and "Right" in "x_px" distribution function
 # NOTE: We only use the values from the right distribution function as if we inherit

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/plotting/plot_2D_poynting_flux.py

@@ -10,7 +10,7 @@ input_dir = Path("datasets/3_5_Gaussian_beam")
 ds = sdfxr.open_dataset(input_dir / "0001.sdf")
 
 # Convert the x and y coords to microns
-ds = ds.epoch.rescale_coords(1e6, "µm", ["X_Grid_mid", "Y_Grid_mid"])
+ds = ds.epoch.rescale_coords("µm", ["X_Grid_mid", "Y_Grid_mid"], 1e6)
 
 # Calculate Poynting flux magnitude
 flux_magnitude = np.sqrt(

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/plotting/plot_time_temperature.py

@@ -13,7 +13,7 @@ input_dir = Path("datasets/4_2_self_heating")
 ds = sdfxr.open_mfdataset(input_dir)
 
 # Convert the time to femtoseconds
-ds = ds.epoch.rescale_coords(1e15, "fs", "time")
+ds = ds.epoch.rescale_coords("fs", "time", 1e15)
 
 # Averate temperature over all spatial cells at each time-step
 kB = 1.380649e-23

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/epoch_workshop_2026/plotting/plot_x_px_scatter.py

@@ -9,7 +9,7 @@ input_dir = Path("datasets/5_2_subsets")
 ds = sdfxr.open_dataset(input_dir / "0000.sdf", keep_particles=True)
 
 ds = ds.epoch.rescale_coords(
-    1e6, "µm", ["X_Particles_subset_Refluxers_Electron", "Y_Target_mid"]
+    "µm", ["X_Particles_subset_Refluxers_Electron", "Y_Target_mid"], 1e6
 )
 
 x = ds["X_Particles_subset_Refluxers_Electron"]

sdf_xarray-0.9.0/docs/epoch_workshop_2026/workshop_details.md

@@ -0,0 +1,7 @@
+# Workshop Details
+
+The workshop consisted of a series of talks and practical sessions delivered at the University of York over the course of 2 days (23-24 April). It covered setting up both  `EPOCH` and `sdf-xarray` along with delving into the details of how PIC codes work. The timetable is located [here](https://epochpic.github.io/quickstart/epoch_workshop_2026.html) and recordings of the presentations can be found on [YouTube](https://youtube.com/playlist?list=PLXfw-OhAIheRtwRqrGZvMqFuQ2GrIBLVl).
+
+## Live Demo
+
+As part of the start of the workshop a live demo covering the key features of `sdf-xarray` was performed. You can download this jupyter notebook in order to follow along <path:live_demo.ipynb>. Alternatively, if you have a Google account you can use the online [Google Colab notebook](https://colab.research.google.com/drive/1FTKj6gUizKjcL1ZSK4sNKlWKs5mzrD2g?usp=sharing).

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/index.md

@@ -21,6 +21,7 @@ why_sdf_xarray.md
 loading_data.md
 understanding_datasets.md
 unit_conversion.md
+plotting.md
 animation.md
 ```
 

sdf_xarray-0.9.0/docs/plotting.md

@@ -0,0 +1,95 @@
+---
+file_format: mystnb
+kernelspec:
+  name: python3
+---
+
+# Plotting
+
+You can plot datasets using
+[`xarray.DataArray.epoch.plot`](project:#sdf_xarray.dataarray_accessor.EpochAccessor.plot).
+This is a custom <project:#sdf_xarray> plotting routine that builds on top of
+<inv:#xarray.DataArray.plot>, so you keep the familiar <inv:#xarray> plotting
+behaviour while using <project:#sdf_xarray> conveniences (see 
+[here](project:#sdf_xarray.dataarray_accessor.EpochAccessor.plot) for details).
+Under the hood, plotting is still handled by <inv:#matplotlib>, which means you
+can use the full <inv:#matplotlib> API to customise your figure.
+
+| Dimensions | Base plotting function        |        Notes          |
+| ---------- | ----------------------------- | --------------------- |
+|     1      | <inv:#xarray.plot.line>       |                       |
+|     2      | <inv:#xarray.plot.pcolormesh> |                       |
+|     3      | <inv:#mpl_toolkits.mplot3d.axes3d.Axes3D.voxels>| Custom functionality built for <project:#sdf_xarray> |
+|     >3     | <inv:#xarray.plot.hist>       | Not fully implemented |
+
+## Line plots
+
+One dimensional data will be plotted as a line. Multiple lines can be easily
+plotted at the same time.
+
+```{code-cell} ipython3
+import sdf_xarray as sdfxr
+import matplotlib.pyplot as plt
+ds = sdfxr.open_dataset("tutorial_dataset_1d/0010.sdf")
+da_1 = ds["Derived_Number_Density_Electron"]
+da_2 = ds["Derived_Number_Density_Ion"]
+da_1.epoch.plot(label = "Electron")
+da_2.epoch.plot(label = "Ion")
+plt.legend()
+plt.show()
+```
+
+## Mesh plots
+
+Two dimensional data will be plotted as a mesh.
+
+```{code-cell} ipython3
+ds = sdfxr.open_dataset("tutorial_dataset_2d/0010.sdf")
+da = ds["Derived_Number_Density_Electron"]
+da.epoch.plot()
+plt.show()
+```
+
+## Voxel plots
+
+Three dimensional data will be plotted using voxels. This behaviour is not native
+to <inv:#xarray>, it has been custom built specifically for this package.
+
+```python
+ds = sdfxr.open_dataset("tutorial_dataset_3d/0000.sdf")
+da = ds["Derived_Number_Density"]
+da.epoch.plot(vmin = 1e27)
+plt.show()
+```
+
+![voxel_plot_high_res](./figures/voxel_plot_high_res.png)
+
+```{warning}
+Voxel plots can be extremely computationally expensive and may take longer than
+expected to plot.
+```
+
+Because voxel plots are very expensive, even with relatively small data arrays,
+plotting can be sped up by using [`xarray.DataArray.epoch.resize`](project:#sdf_xarray.dataarray_accessor.EpochAccessor.resize),
+which uses interpolation to reduce the resolution.
+
+```{code-cell} ipython3
+ds = sdfxr.open_dataset("tutorial_dataset_3d/0005.sdf")
+da = ds["Derived_Number_Density"]
+da_resized = da.epoch.resize((20, 20, 20))
+da_resized.epoch.plot(vmin = 1e27)
+plt.show()
+```
+
+## Histograms
+
+When the data array has four or more dimensions (or is specified by the user), it will
+be plotted as a histogram.
+
+```{code-cell} ipython3
+ds = sdfxr.open_dataset("tutorial_dataset_3d/0005.sdf")
+da = ds["Derived_Number_Density"]
+
+da.epoch.plot(hist = True)
+plt.show()
+```

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/understanding_datasets.md

@@ -29,23 +29,10 @@ ds = sdfxr.open_mfdataset("tutorial_dataset_1d/*.sdf")
 ds["Electric_Field_Ex"]
 ```
 
-## Plotting
+## Visualisation
 
-You can plot datasets using
-[`xarray.DataArray.epoch.plot`](project:#sdf_xarray.dataarray_accessor.EpochAccessor.plot).
-This is a custom <project:#sdf_xarray> plotting routine that builds on top of
-<inv:#xarray.DataArray.plot>, so you keep the familiar <inv:#xarray> plotting
-behaviour while using <project:#sdf_xarray> conveniences (see 
-[here](project:#sdf_xarray.dataarray_accessor.EpochAccessor.plot) for details).
-Under the hood, plotting is still handled by <inv:#matplotlib>, which means you
-can use the full <inv:#matplotlib> API to customise your figure.
-
-```{code-cell} ipython3
-# This is discretized in both space and time
-ds["Electric_Field_Ex"].epoch.plot()
-plt.title("Electric field along the x-axis")
-plt.show()
-```
+`sdf-xarray` has built-in plotting and animation functionality, more details can
+be found on [](./plotting.md) and [](./animation.md).
 
 ## Dimension slicing
 
@@ -56,6 +43,7 @@ for easy indexing. To quickly determine the number of time steps available,
 you can check the size of the time dimension.
 
 ```{code-cell} ipython3
+ds = sdfxr.open_mfdataset("tutorial_dataset_1d/*.sdf")
 # This corresponds to the number of individual SDF files loaded
 print(f"There are a total of {ds['time'].size} time steps")
 
@@ -66,7 +54,7 @@ print(f"The time at the 20th simulation step is {sim_time:.2e} s")
 
 You can select and extract a single simulation snapshot using the integer
 index of the time step with the <inv:#xarray.Dataset.isel> function. This can be
-done by passsing the index to the `time` parameter (e.g., `time=0` for
+done by passing the index to the `time` parameter (e.g., `time=0` for
 the first snapshot).
 
 ```{code-cell} ipython3
@@ -74,7 +62,7 @@ ds["Electric_Field_Ex"].isel(time=20)
 ```
 
 We can also use the <inv:#xarray.Dataset.sel> function if you wish to pass a
-value intead of an index.
+value instead of an index.
 
 ```{tip}
 If you know roughly what time you wish to select but not the exact value
@@ -122,6 +110,51 @@ print(f"Total particle energy absorbed: {ds["Total_Particle_Energy_in_Simulation
 print(f"The laser absorption fraction: {ds["Laser_Absorption_Fraction_in_Simulation"][-1].values:.1f} %")
 ```
 
+## Limit
+
+The `.epoch` accessor has the ability to "limit" data arrays. This is a helper
+function which applies <inv:#xarray.DataArray.where> to the upper and lower bounds
+of each dimension. You can use `None` to preserve the existing limit.
+
+```{code-cell} ipython3
+ds = sdfxr.open_dataset("tutorial_dataset_2d/0020.sdf")
+da = ds["Derived_Number_Density_Electron"]
+da.epoch.plot()
+print(f"Orignal shape: {da.shape}")
+plt.show()
+```
+
+```{code-cell} ipython3
+da_limit = da.epoch.limit(((-2e-6, 2e-6), (None, 0)))
+da_limit.epoch.plot()
+print(f"New shape: {da_limit.shape}")
+plt.show()
+```
+
+## Resize
+
+The data arrays can be interpolated with <project:#sdf_xarray.dataarray_accessor.EpochAccessor.resize>
+to either increase or decrease the resolution while mantaining the general structure.
+Decreasing the size may be useful to download large data from HPCs to local machines
+or ploting with [](./plotting.md#voxel-plots). Be aware that the original
+data will not be preserved.
+
+```{code-cell} ipython3
+ds = sdfxr.open_dataset("tutorial_dataset_2d/0020.sdf")
+da = ds["Derived_Number_Density_Electron"]
+
+da.epoch.plot()
+print(f"Orignal shape: {da.shape}")
+plt.show()
+```
+
+```{code-cell} ipython3
+da_resized = da.epoch.resize((50, 50))
+da_resized.epoch.plot()
+print(f"New shape: {da_resized.shape}")
+plt.show()
+```
+
 ## Visualisation on HPC Machines
 
 In many cases you will be running EPOCH simulations via a HPC cluster and your

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/unit_conversion.md

@@ -30,7 +30,7 @@ plt.rcParams.update({
 For simple scaling and unit relabelling of coordinates (e.g., converting meters to microns),
 the most straightforward approach is to use the [`xarray.Dataset.epoch.rescale_coords`](project:#sdf_xarray.dataset_accessor.EpochAccessor.rescale_coords) dataset accessor.
 This function scales the coordinate values by a given multiplier and updates the
-`"units"` attribute in one step.
+`"units"` attribute in one step. If the multiplier is not specified then the conversion parameter is inferred from the units and converted automatically using `pint` (see [](#unit-conversion-with-pint-xarray) for details).
 
 ### Rescaling grid coordinates
 
@@ -41,7 +41,7 @@ We can use the [`xarray.Dataset.epoch.rescale_coords`](project:#sdf_xarray.datas
 fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(16, 6))
 
 ds = sdfxr.open_mfdataset("tutorial_dataset_2d/*.sdf")
-ds_in_microns = ds.epoch.rescale_coords(1e6, "µm", ["X_Grid_mid", "Y_Grid_mid"])
+ds_in_microns = ds.epoch.rescale_coords("µm", ["X_Grid_mid", "Y_Grid_mid"], 1e6)
 
 ds["Derived_Number_Density_Electron"].isel(time=0).epoch.plot(ax=ax1)
 ax1.set_title("Original X Coordinate (m)")
@@ -55,7 +55,7 @@ fig.tight_layout()
 ### Rescaling time coordinate
 
 We can also use the [`xarray.Dataset.epoch.rescale_coords`](project:#sdf_xarray.dataset_accessor.EpochAccessor.rescale_coords) method to convert the time coordinate from
-seconds (`s`) to femto-seconds (`fs`) by applying a multiplier of `1e15`.
+seconds (`s`) to femto-seconds (`fs`).
 
 ```{code-cell} ipython3
 ds = sdfxr.open_mfdataset("tutorial_dataset_2d/*.sdf")
@@ -63,7 +63,7 @@ ds["time"]
 ```
 
 ```{code-cell} ipython3
-ds = ds.epoch.rescale_coords(1e15, "fs", "time")
+ds = ds.epoch.rescale_coords("fs", "time")
 ds["time"]
 ```
 

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/docs/why_sdf_xarray.md

@@ -7,7 +7,7 @@ There are several benefits to using this package over the [`sdf_helper`](https:/
 - [data to numpy and pandas](https://docs.xarray.dev/en/stable/user-guide/pandas.html) - Data can be easily converted to [numpy](https://www.numpy.org/) or [pandas](https://pandas.pydata.org/) if you prefer to work with raw data.
 - [data labelling](project:understanding_datasets.md) - All SDF files come with a bunch of information about them such as the `dimensions`, `units` and `name` of the variable. In <inv:#xarray> these are added to each variable making it much easier to read.
 - [lazy/partial loading](https://docs.xarray.dev/en/stable/internals/internal-design.html#lazy-loading) - This reduces the RAM requirements when loading large SDF files as we only load the actual array values when they are needed.
-- [plotting](project:understanding_datasets.md#plotting) - You can easily plot variables across several dimensions including time using built in <inv:#matplotlib> support.
+- [plotting](project:plotting.md) - You can easily plot variables across several dimensions including time using built in <inv:#matplotlib> support.
 - [animating](project:animation.md) - Visualise your data across time or slice through another dimension and make beautiful GIFs for your presentations.
 - [deck loading](loading-input-deck) - <project:#sdf_xarray> automatically loads in the associated `input.deck` utilised when creating the simulation so that you can access setup variables.
 - [easier unit conversion](project:unit_conversion.md) - Convert dimensions and variable units to other types with the help of [`pint`](https://pint.readthedocs.io/en/stable).

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/pyproject.toml

@@ -69,6 +69,7 @@ docs = [
   "sphinx-design>=0.7.0",
   "pooch>=1.8.2",
   "tqdm",
+  "scipy",
 ]
 build = ["cibuildwheel[uv]"]
 lint = ["ruff"]
@@ -78,6 +79,9 @@ test = [
   "matplotlib",
   "pooch>=1.8.2",
   "tqdm",
+  "pint",
+  "pint-xarray",
+  "scipy",
 ]
 
 [project.entry-points."xarray.backends"]
@@ -131,6 +135,10 @@ ignore = [
   "PLR2004", # magic-comparison
   "B9",      # flake8-bugbear opinionated warnings
   "PLR0913", # remove maximum number of arguments in a function
+  "PLR0917", # remove maximum number of positional arguments in a function
+  "BLE001",  # allow capturing of base exception
 ]
-exclude = ["**/_version.py"]
-# Auto-generated file
+exclude = ["**/_version.py"] # Auto-generated file
+
+[tool.pytest.ini_options]
+filterwarnings = ['ignore::UserWarning']

{sdf_xarray-0.7.4 → sdf_xarray-0.9.0}/src/sdf_xarray/_version.py

@@ -18,7 +18,7 @@ version_tuple: tuple[int | str, ...]
 commit_id: str | None
 __commit_id__: str | None
 
-__version__ = version = '0.7.4'
-__version_tuple__ = version_tuple = (0, 7, 4)
+__version__ = version = '0.9.0'
+__version_tuple__ = version_tuple = (0, 9, 0)
 
-__commit_id__ = commit_id = 'gc521760ee'
+__commit_id__ = commit_id = 'g65bb865da'