Jeremy Sanders has provided a Python script which runs the SPEX 2 program and extracts the spectra into files in the same format as APEC. The script requires the Python numpy and pyfits packages. The resulting files can be used in xspec by setting APECROOT then using the apec or vapec model. Available are the script and example outputs spex_2.00.11_coco.fits and spex_2.00.11_line.fits. To use the example files in xspec use the command "xset APECROOT /path/spex_2.00.11" where /path is the directory path to the files. The temperatures tabulated and energy range covered can be altered by editing the script to change the values in the temperature_list array and variables contminenergy, contmaxenergy, and contenergysteps. Updated on Oct 24, 2008 to replace numarray by numpy and work around a problem in xspec when there is no pseudo-continuum.
Keith Arnaud, Lab. for High Energy Astrophysics, NASA/Goddard Space Flight Center
HEASARC Home | Observatories | Archive | Calibration | Software | Tools | Students/Teachers/Public
Last modified: Tuesday, 06-Oct-2009 13:52:59 EDT