********
Examples
********


Running Roland
==============
You need to supply |roland| with 1) a model atmosphere, 2) a linelist (or multiple linelists), and 3) the set of stellar
parameters and abundances.

Here's a simple example:

    >>> from roland import spectrumizer
    >>> atmod = utils.testdir()+'t3710g083m+000.mod'
    >>> linelists = [utils.datadir()+'moogatom.gz',utils.datadir()+'moogh2o.gz']
    >>> moog = spectrumizer.spectrumizer('moog')
    >>> flux,cont,wave = mood.synthesize(3710.0,0.83,0.0,atmod=atmod,linelists=linelists,wrange=[6700,6800])

Now plot the spectrum:

    >>> import matplotlib.pyplot as plt
    >>> plt.plot(wave,flux)

It should look like this.

.. image:: spectrum_example.png
  :width: 600
  :alt: Example MOOG synthetic spectrum

Abundances
----------
	
You can modify the global alpha abundance with `am` or individual abundances with `elems`.  The `elems` parameter
takes a list of [element name, abundance] pairs, where the abundance should be in the form [X/M], where M is the
overall metallicity that is used to scale the individual abundances.  For example, ``elems=[['Mg',0.55],['Ba',-0.15]]``
means a relative Magnesium abundance of +0.55 and a relative Barium abundance of -0.15.

Let's try it out:

    >>> flux2,cont2,wave2 = moog.synthesize(3710.0,0.83,0.0,atmod=atmod,linelists=linelists,wrange=[6700,6800],elems=[['Mg',0.55],['Ba',-0.15]])
    >>> plt.plot(wave,flux)
    >>> plt.plot(wave2,flux2)
    >>> plt.xlim(6700,6800)