crea/icat std EFOSC_Spectrum*obj*bdf
$ echo leave only spectrum of std star, so delete this line as well >> std.cat


$vi std.cat

$ echo leave only table of std star, so delete this line as well  > ftables.txt
$ find /scisoft/midas/calib/data/spec/flux/ -name \*        >> ftables.txt
$ find /scisoft/midas/calib/data/spec/ctio/ -name \*        >> ftables.txt

$vi ftables.txt

$cat ftables.txt | awk '{print "cp "$1" stdflux.tbl"}' > cp.sh
$cat cp.sh
$sh cp.sh


$ head -2 std.cat | tail -1 > y
$ cat y | awk '{print $1}' > x
$ cat x | awk '{print "cp "$1" standard.bdf"}' > cp.sh
$ sh cp.sh

!!!!!!!!!!!!!!!!!!!!!!!
!!
!! don't make selction of extinction curve interactive...
!! ... just correct for atmospheric extinction
!!

$ cp /scisoft/midas/calib/data/spec/exti/atmoexan.tbl extinct.tbl



SET/LONG FLUXTAB=stdflux EXTAB=extinct 
PLOT/FLUX
write/out
write/out This is the flux of the standard (note the y-scale)
write/out

get/gcur

!goto end

EXTINCTION/LONG standard standard_extinct



INTEGRATE/LONG standard_extinct
write/out 
write/out Click on points you want to delete ...
write/out

edit/flux



set/long FITYP=poly
set/long FITD=12

RESPONSE/LONG 
write/out
write/out This is the response curve
write/out
get/gcur

PLOT/RESPONSE
write/out
write/out again without points ...
write/out

get/gcur


CALIBRATE/FLUX standard_extinct standard_flux

plot/flux
overplot standard_flux

write/out
write/out points are reference, line is measured SED
write/out
get/gcur

write/out 
write/out response curve obtained with these parameters:
write/out
show/long f

end: