pro rd2xyl
filein=""
h = headfits('acisf00020N001_evt2.fits',EXTEN=1) 
; hprint, h
; Build an astrometry structure for this observation.
crpixX = fxpar( h, 'TCRPX11' )
hprint, crpixX
crpixY = fxpar( h, 'TCRPX12' )
hprint, crpixY
crvalX = fxpar( h, 'TCRVL11' )
hprint, crvalX
crvalY = fxpar( h, 'TCRVL12' )
hprint, crvalY
cdeltX = fxpar( h, 'TCDLT11' )
hprint, cdeltX
cdeltY = fxpar( h, 'TCDLT12' )
hprint, cdeltY
ctypeX = fxpar( h, 'TCTYP11' )
hprint, ctypeX
ctypeY = fxpar( h, 'TCTYP12' )
hprint, ctypeY
make_astr, astr, DELTA=[cdeltX,cdeltY], CTYPE=[ctypeX,ctypeY], $
CRPIX=[crpixX,crpixY], CRVAL=[crvalX,crvalY]
read," Enter RA-Decl filename:", filein
openr, 1, filein
;while NOT(EOF(in)) DO BEGIN
for count=0,17 do begin
readf, 1,ra,dec
ad2xy, ra,dec, astr, refx, refy
refx=refx+1
refy=refy+1
print, refx,refy
endfor
;endwhile
close, 1
end