In [1]:
import apogee.modelspec.moog
from apogee.modelatm import atlas9
import apogee.spec.plot as splot
Let's generate a spectrum at $T_{\mathrm{eff}} = 4500 K$, $\log g = 2.5$ and solar abundances
In [2]:
atm= atlas9.Atlas9Atmosphere(teff=4500.,logg=2.5,metals=0.,am=0.,cm=0.)
synspec_jonres= apogee.modelspec.moog.synth(modelatm=atm,
linelist='moog.201405291900.vac',lsf='combo',cont='aspcap',
vmacro=6.,isotopes='solar',
dw=0.1)
synspec_hires= apogee.modelspec.moog.synth(modelatm=atm,
linelist='moog.201405291900.vac',lsf='combo',cont='aspcap',
vmacro=6.,isotopes='solar',
dw=0.01)
In [3]:
for panel in apogee.spec.plot.highres(synspec_jonres[0],synspec_hires[0],(synspec_hires[0]-synspec_jonres[0])*10.+0.25,
yrange=[0.1,1.2],color=['b','g','gold'],labelLines=True,cleanZero=False,fig_width=14.):
show()
Residuals everywhere
In [4]:
print numpy.nanstd(synspec_hires[0]-synspec_jonres[0]), numpy.nanmax(numpy.fabs(synspec_hires[0]-synspec_jonres[0]))
In [5]:
plot(synspec_hires[0]-synspec_jonres[0])
Out[5]:
In [5]: