# Objective

This folder contains vti data to reproduce the figures in [Lautenbach-Grauer-2018].

# Usage

Open vti files in any VTK-viewer or import them into your plotting script. Visit
https://www.vtk.org/Wiki/VTK_Tools and
https://lorensen.github.io/VTKExamples/site/ for details.

# Reproduce Figures

Go to fluid_data

## Figure 3
The three center columns are simply the field `jz` and the absolute values of
`uz_e` and `uz_i`.  The outer columns can be obtained by applying the thresholds
in table 3. The time in the filename represents $0.1 \omega_{c,i}^{-1}$. Hence,
the depicted times are from files `t={0,{4,8,12,16,20,24,28,32,36,40}0}.vti`.

## Figure 4
The field `jz` of the four snapshots.

## Figure 5
Due to normalization, it is sufficient to calculate the sum from field `by` and
multiply by $dx=0.0490875$.


# Investigate Distribution Functions

We provide data of the distribution function at selected times (physical times
as in Figure 3).

+ Go to kinetic_data and unpack
+ Use the python-script to extract PSD at single locations

Note that the distribution function is only available within kinetic regions. It
can be constructed from fluid moments everywhere else.