Hi,<br><br>I have some spectra that were reduced in iraf, and I would 
like to read them in to python to do some fitting that I already have 
code for. However, I am having trouble reading in the wavelength 
calibration. The spectra are echelle, and the wavelengths are defined 
with this thing in the header that looks like this:<br>

<br>&quot;<br>WAT2_001= &#39;wtype=multispec spec1 = &quot;1 32 2 10672.221264362 0.086046403007761 20&#39;<br>WAT2_002= &#39;46 0. 6.65 26.40 1. 0. 1 5 1. 2046. 10761.5346117191 87.998952480064&#39;<br>WAT2_003= &#39; -1.32409633589482 -0.0165054046288388 -0.00680394594704411&quot; spec2 =&#39;<br>



WAT2_004= &#39; &quot;2 33 2 10348.8582697 0.083445037904806 2046 0. 34.60 58.07 1. 0. 1&#39;<br>WAT2_005= &#39; 5 1. 2046. 10435.4699172677 85.3379030819273 -1.2833101652519 -0.01&#39;<br>WAT2_006= &#39;53518242619121 -0.0057861452751027&quot; spec3 = &quot;3 34 2 10043.881872825 &#39;<br>



...<br>&quot;<br><br>The numbers give the conversion from pixel coordinates to wavelength coordinates.<br>I
 was wondering if there was a python function that could correctly read 
all that in, parse it, and give me a spectrum in flux vs wavelength? I 
have tried playing with pywcs but that seems only to work for images?<br>

<br>Cheers,<br clear="all"><br>Kevin Gullikson