!
! prepare and wave-calib 1 frame
!
!


write/out
write/out {p1} must have the .fits extension
write/out otherwise you see an error
write/out

$ls {p1}                   >  in.cat
$ls {p1} | sed s/fits/bdf/ > out.cat

indisk/fits in.cat out.cat


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!
!! rotate, to have lambda from sx to dx
!!
!!    overwrites original images
!!    restores descriptors start and step
!!


$cat out.cat | awk '{print "rotate/clock "$1" r"$1}' > rclock.prg

write/out rotating image ... 
@@ rclock

$cat rclock.prg | awk '{print "mv "$3,$2}' > mv.sh

$sh mv.sh

$cat out.cat | awk '{print "crea/icat obj "$1}' > ccat.prg
@@ ccat

exec/cat write/descr obj.cat start 1,1
exec/cat write/descr obj.cat  step 1,1
exec/cat copy/dd obj.cat ESO.TEL.AIRM.END obj.cat o_airm

$cat out.cat | awk '{print "load "$1" scale=-2"}'  > ld.prg
$cat out.cat | awk '{print "@a autocuts "$1    }' >> ld.prg





!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!
!! subtract bias from all images
!!
exec/cat compute/ima obj.cat = obj.cat - bias



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!
!! trim objects and HeAr
!!
write/out trimming spectrum


prepare/long obj.cat obj

!! restore image name (overwrite old frame)

$cat out.cat | awk '{print "mv obj1.bdf "$1}' > mv.sh
$sh mv.sh


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!
!! divide objects and HeAr by master flatfield
!!   and check final result
!!
exec/cat compute/ima obj.cat = obj.cat / masterflat



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!
!! apply wavelength calibration, method is "linear" (default)
!!
!! I have to rename image, because the EFOSC2 convention of
!! filenames (with xxx.yy.bdf) confused midas
!!

$cat out.cat | awk '{print "mv "$1" rebinme.bdf"}' > mv.sh
$sh mv.sh

rebin/long rebinme rebinned
$cat out.cat | awk '{print "mv rebinned.bdf "$1}' > mv.sh
$sh mv.sh

!! load resulting spectrum (ld.prg has been created above)

@@ ld


end: