Overview

Parameter files can be read in to Paramfile objects, manipulated, and written back out. See Working with Parameter Files for details.

>>> par = mcfost.Paramfile('example/ref2.19.para')
>>> par.distance = 200
>>> par.writeto('modified.par')

Once MCFOST calculations have been done for a given parameter file, their results may be examined and plotted. See Working with Model Results for details.

>>> res = mcfost.ModelResults('example')
>>> res.describe()
Model results in examples for ref2.19.para
  Model has 15 inclinations from 60.0 to 90.0
  SED computed from 0.1 - 3500.0 microns using 30 wavelengths
  Images computed for 1 wavelengths: [ 0.8] micron
>>> res.sed.plot(title='Example SED plot', nlabels=4)

(Source code, png, hires.png, pdf)

You presumably also have observed data you wish to compare to, and those can be loaded and displayed as well. See Working with Model Results for details, including the required file formats. The interfaces for working with model results and observations are similar:

>>> obs = mcfost.observations('example_obs')
>>> obs.describe()
>>> obs.sed.plot()

There are included functions for calculating values for model fitting, currently supporting SEDs only. Image fitting is a work in progress. The long term goal is to provide an interface for MCMC model fitting via emcee.