Usage of dsharp_helper package
[1]:
import numpy as np
import matplotlib.pyplot as plt
import astropy.constants as c
import astropy.units as u
%matplotlib inline
import dsharp_helper as dh
[2]:
disk = 'IMLup'
Files
You can download & get the path to a file name like this (careful: CO files are large)
[3]:
fn_cont = dh.get_datafile(disk, suffix='continuum')
fn_sed = dh.get_datafile(disk, type='SED')
fn_profile = dh.get_datafile(disk, type='profile')
Downloading file 'IMLup_continuum.fits' ... Done!
Downloading file 'IMLup.SED.txt' ... Done!
Downloading file 'IMLup.profile.txt' ... Done!
Continuum Image
[4]:
f, ax = dh.plot_DHSARP_continuum(disk=disk)
f.savefig(disk + '_cont.png')
IMLup: vmin = 1.5e-05, vmax = 0.0053
/home/docs/checkouts/readthedocs.org/user_builds/dsharp-helper/conda/latest/lib/python3.9/site-packages/dsharp_helper/dsharp_helper.py:506: MatplotlibDeprecationWarning: shading='flat' when X and Y have the same dimensions as C is deprecated since 3.3. Either specify the corners of the quadrilaterals with X and Y, or pass shading='auto', 'nearest' or 'gouraud', or set rcParams['pcolor.shading']. This will become an error two minor releases later.
getattr(ax, fct)(
Profile
[5]:
f, ax = dh.plot_profile(disk)
f.savefig(disk + '_prof.png')
SED
[6]:
f, ax = dh.plot_sed(disk)
f.savefig(disk + '_sed.png')
