debris_disk

Produce synthetic images and visibilities of disk structure from a radial, vertical, and scale height functional form. You can produce an image such as

import debris_disk as DD

obs_params = {'nu' : [3.451e11], # hz
              'imres' : 0.01, # arcsec / pixel
              'PA' : 0., # degrees
              'distance' : 100} # parsecs

rad_params = {'alpha' : 2.0,
              'Rin' : 22 * DD.constants.AU, # cm
              'Rout' : 42 * DD.constants.AU} # cm

model = DD.Disk(Lstar=1., # L_sun
                sigma_crit=1e-5, # g/cm^2
                inc=88.5, # Degrees
                radial_func='powerlaw', # options: double_powerlaw,
                                        # triple_powerlaw, gaussian
                radial_params=rad_params,
                vert_params={'Hc' : 0.03,
                             'Rc' : 1.,
                             'psi' : 1.,
                vert_func='gaussian',
                obs=obs_params)

model.square()
model.rotate()
model.save('example.fits')

From a disk model, produce synthetic visibilities as

vis=DD.UVDataset('data_directory') 
vis.sample(model)

where 'data_directory' is the filepath to a directory containing only uvfits files.

From a disk model, calculate a $\chi^2$ value from visibilities like

vis.chi2(mod)