Making slice plots (Beginner)

Import packages

import numpy                          as np
import matplotlib.pyplot              as plt

# import plons scripts
import plons
import plons.SmoothingKernelScript    as sk
import plons.ConversionFactors_cgs    as cgs
import plons.Plotting                 as plot

Setting information about data

prefix = "wind"
loc = "/STER/matse/Papers/Esseldeurs+2023/Phantom/High/binary6Lucy/"
dump = loc+"wind_00600"

Loading setup and dump

setup     = plons.LoadSetup(loc, prefix)

dumpData  = plons.LoadFullDump(dump, setup)

Plotting density in the orbital plane

Plotting around the CoM

fig, ax = plt.subplots(1, figsize=(4, 4))
cbar = plot.SlicePlot2D(ax, dumpData, setup, xlim=(-30, 30), ylim=(-30, 30))

Plotting around the companion

fig, ax = plt.subplots(1, figsize=(4, 4))
cbar = plot.SlicePlot2D(ax, dumpData, setup, xlim=(0, 10), ylim=(-5, 5), rotate=True)

Making your own plane on which to smooth (as a meshgrid)

n = 200
x = np.linspace(-30, 30, n)*cgs.au
y = np.linspace(-30, 30, n)*cgs.au
X, Y = np.meshgrid(x, y)
Z    = np.zeros_like(X)

Smoothing the data on the plane

smooth = sk.smoothMesh(X, Y, Z, dumpData, ['rho'])

Plotting the plane

fig, ax = plt.subplots(1, figsize=(4, 4))
plot.plotSlice(ax, X, Y, smooth, 'rho', logplot=True, cmap = plt.cm.get_cmap('inferno'), clim=(-17, -14))
cAGB, cComp = plot.plotSink(ax, dumpData, setup)