1) reads spectrum array from this file;
2) reads keywords from the header of this file and using these keywords computes the wavelength array according to the formula I provided;
3) plots on the screen spectrum vs. wavelength for the wavelength
file contains the header and a single array (spectrum). Wavelengths can be calculated from the keywords NAXIS1, CDELT1, CRVAL1 in the header:
Wavelength = exp( array[0:NAXIS1-1] * CDELT1+CRVAL1 )
"array[0:NAXIS1-1]" is an integer array from 0 to NAXIS1-1 with a step of 1.
Best Answer