Python library for protoplanetary disk imaging

In this notebook I show how to use the plot functions of the disc_class. I created this class for my summer student project in 2016. I found there is a lack of python tools for visualize protoplanetary discs, thus I decided to share my raw functions for plotting 2D FARGO uotput discs using python. These are raw functions just to show the basic concept how to plot discs with python. I not recomend to use this functions as a blackbox, however I sugest to change the code to create nice plots.

The disc class:

This is a python class describing proroplanetary discs. The disc objects requires 3 argumets:

• r: the radial grid (simulation units)
• $\varphi$: azimutal grid (simulation umits)
• fname: name of the imput file

thus a disc can be define the following way:

There are some plot function to display the disc or some disc propertyes:

• radial plot
• polar plot
• radial profile
• azimutal profile
• combined plots

Radial plot

The function radial_plot display the disc radial grid values against the azimutal grid values. And position the maximum to $\pi$.

Arguments:

• output: string, optional, name of the output picture.
• g: integer, optional, the so called gamma correction that raise the data to the power of g
• norm: disc-class object, optional, other (initial) disc configuration to normalize the plot

Examples:

Plot disc1 without normalization and normalized with normdisc.

Plot disc2 with normalization and with g=3 gamma correction.

Polar plot

The function polar_plot transform the FARGO output to polar coordinates. Display the transformed disk, with a blank ring at the center 10% of the maximal radius. And position the density maximum to $\pi$.

Arguments:

• output: string, optional, name of the output picture.
• g: integer, optional, the so called gamma correction that raise the data to the power of g
• norm: disc-class object, optional, other (initial) disc configuration to normalize the plot
• polar_plot_radius: int, optional, it define the plot box range
• log_transform: se description below

Examples:

Plot disc2 with g=2 gamma correction and disc1 without correction. Both plots are normalized with normdisc. The second plot is a zoomed plot with a polar_plot_radius=200.

As shown the zooming is not satisfying because in the example simulations I used logarythmic binning for the radial grid. This is a common practice. Thus the disc call first function converts the simulation radila units to phisical dimension data. It is written for my dataset, thus you need to rewrite it. Using this function, the polar plots will be much nicer, and these plots show the 'real' simulation results. If you wish to use this transformation, you need to set the log_transform variable to True. Let's see some examples!

Plot disc2 with the transformed values and disc1 with the transformation with a polar_plot_radius=30.

You can use this feature later on with the other plot types if reasonable. It is not available on the radial_plot yet.

Radial profile

The function azimutaly_averaged_radial_profile plots the azimutaly averaged radial profile as its name suggest. A typical function.

Arguments:

• output: string, optional, name of the output picture.
• norm: disc-class object, optional, other (initial) disc configuration to normalize the plot
• log_transform: False, optional. As described above set this value True to change the grid units.

Examples:

Plot disc1 with no arguments and disc2 normalized and transformed to physical radial units.

Azimutal profile

The function azimutal_profile_around_maxima plots the azimutal profile around the dr radius of the maximum. First the function finds the maximum index, and then radially average the azimautal profiles from maximum-dr to maximum+dr.

Arguments:

• output: string, optional, name of the output picture.
• norm: disc-class object, optional, other (initial) disc configuration to normalize the plot
• dr: int, initial value 10. This parameter define the radial widt of the averaged area

Examples:

Normed plot of the disk1 and disk2 with dr=5 and dr=20, respectively.

Radial and azimutal profile

The function radial_and_azimutal_profile plot the prewious plots next to each other.

Arguments:

• output: string, optional, name of the output picture.
• norm: disc-class object, optional, other (initial) disc configuration to normalize the plot
• dr: int, initial value 10. This parameter define the radial widt of the averaged area
• log_transform: False, optional. As described above set this value True to change the grid units.

Examples

Plot disc1 with dr=10 and log_transform=True.

Frame plot

The function frame_plot plots the radial and azimutal profile and the polar plot into one plot.

Arguments:

• output: string, optional, name of the output picture.
• norm: disc-class object, optional, other (initial) disc configuration to normalize the plot
• dr: int, initial value 10. This parameter define the radial widt of the averaged area
• g: integer, optional, the so called gamma correction that raise the data to the power of g
• polar_plot_radius: int, optional, it define the plot box range
• log_transform: False, optional. As described above set this value True to change the grid units.
• plot_title: Optional, title of the plot

Example