#!/usr/bin/python


import pyfits
import numpy


filein = 'HARPS.Archive_HD49933.fits'
fileout = 'HARPS.Archive_HD49933_ares.fits'

hdulist = pyfits.open(filein)
header = hdulist[0].header
header.update(key="BITPIX",value=-32)
header.update(key="CDELT1",value=header['CD1_1']*10.)
header.update(key="CD1_1",value=header['CD1_1']*10.)
header.update(key="CRVAL1",value=header['CRVAL1']*10.)
header.update(key="CUNIT1",value='Angstroms')
del header['TTYPE1']
flux = hdulist[0].data.astype(numpy.float32)
pyfits.writeto(fileout, flux, header)