! Procedure EXTRACT.PRG
!
! Creates 39 spectra in MIDAS .tbl format,
! one for each SOPHIE spectral order and 
! corrected for blaze (normalized at order 20).
! Resulting files bear the sequence number of the 
! original observation file and the order number.
! Asks for name of 'e2ds' and associated 'wave' 
! and 'blaze' files
!
! S.Ilovaisky, 9 March 2007
!
define/local sq/i/1/1 0
define/local order/i/1/1/0
define/par P1 ? IMA "e2ds file name ...?"
define/par P2 ? IMA "wave file name...?"
define/par P3 ? IMA "blaze file name...?"
set/format I2
!
order = 0
! 
stat/ima 'P3' row,20
comp/ima blaze_norm = 'P3'/{OUTPUTR(3:3)}
comp/ima e2ds_dbn = 'P1'/blaze_norm
!
copy/it e2ds_dbn e2ds
copy/it 'P2' wave
!
LOOP:
order = 1 + order
!
!
project/tab e2ds e2ds_'order' #'order'
project/tab wave wave_'order' #'order'
!
write/key sq {{P1},OHP.OBS.SEQNUM}
copy/tab wave_'order' sp'sq'_'order'
name/col SP'sq'_'order' #1 :wave
crea/col SP'sq'_'order' :flux "e" E12.6 R*4
copy/tt e2ds_'order' #1 SP'sq'_'order' :flux
IF order .LE. 38 GOTO LOOP
$rm e2ds*.tbl
$rm wave*.tbl
! 
return