disksurf: Extracting the 3D Structure of Protoplanetary Disks
Richard Teague, Charles J. Law, Jane Huang, Feilong Meng
Abstract
disksurf implements the method presented in The Python-based code leverages the opensource GoFish (Teague, 2019) package to read in and interact with FITS data cubes used for essentially all astronomical observations at submillimeter wavelengths. The code also uses the open-source detect_peaks.py routine from Duarte & Watanabe (2021) for peak detection. For a given set of geometric parameters specified by the user, disksurf will return a surface object containing both the disk-centric coordinates of the surface as well as the gas temperature and rotation velocity at those locations. The user is able to 'filter' the returned surface using a variety of clipping and smoothing functions. Several simple analytical forms commonly adopted in the protoplanetary disk literature can then be fit to this surface using either a chi-squared minimization with scipy To verify the 3D geometry of the system is well constrained, disksurf also provides diagnostic functions to plot the emission surface over channel maps of line emission (i.e., the emission morphology for a specific frequency).