False PHOENIX model spectra for Sun-like stars

[MNGuenther]

Either Pysynphot reads the wrong PHOENIX spectra, interpolates them incorrectly, or there is something wrong with PHOENIX itself. The spectra do not match what is expected at the given temperatures. This seems to only affect 'even-hundred' temperatures, the 'odd-hundred' ones seem to work correctly.

pysynphot_phoenix_error_even_temperatures.pdf

pysynphot_phoenix_error_odd_temperatures.pdf

I used the following code to generate these plots:

import numpy as np
import matplotlib.pyplot as plt
import pysynphot as S

#::: odd numbers
plt.figure()
for teff in np.arange(5900,3900,-200):
    sp = S.Icat('phoenix', teff, 0., 4.5)
    plt.plot(sp.wave, sp.flux, label=str(teff))
plt.xlim(3000,10000)
plt.xlabel(r'$T_{eff}$ (K)')
plt.ylabel('flam')
plt.legend(bbox_to_anchor=(1.04, 0.2),  borderaxespad=0)
plt.subplots_adjust(right=0.7)
plt.savefig('pysynphot_phoenix_error_odd_temperatures.pdf')

#::: even numbers
plt.figure()
for teff in np.arange(6000,3900,-200):
    sp = S.Icat('phoenix', teff, 0., 4.5)
    plt.plot(sp.wave, sp.flux, label=str(teff))
plt.xlim(3000,10000)
plt.xlabel(r'$T_{eff}$ (K)')
plt.ylabel('flam')
plt.legend(bbox_to_anchor=(1.04, 0.2),  borderaxespad=0)
plt.subplots_adjust(right=0.7)
plt.savefig('pysynphot_phoenix_error_even_temperatures.pdf')