#!/bin/csh
alias echo "echo > /dev/null"
kappa
figaro
unalias echo
echo " "
echo " *** SLOPE OF CONTINNUM ***"
echo " "
echo " Extract columns near left and right edges of array "
echo " Fits gaussians to two profiles to check position (and FWHM)"
echo " "
echo " Use QUICK_LOOK _raw frames"
echo " "
set date=20030313
echo "Date set to" $date
echo -n "Number of calibration frame: "
set in="$<"
if ($in < 10) set numcal=0000${in}
if ($in > 9 && $in < 100) set numcal=000${in}
if ($in > 99 && $in < 1000) set numcal=00${in}
if ($in > 999) set numcal=0${in}
echo " "
gaiadisp u${date}_${numcal}_raw
echo " Extract two columns, 20-pixels wide, centred on X=100 (left)"
echo " and X=900 (right)"
ystract image=u${date}_${numcal}_raw xstart=90 xend=110 spectrum=left
ystract image=u${date}_${numcal}_raw xstart=890 xend=910 spectrum=right
echo " "
echo " Assuming continuum/pinhole row in range y1=500 to y2=600..."
echo " "
echo " "
echo " *************************************************** "
echo " *************************************************** "
echo " GAUSSIAN FITTING... "
echo " "
fitgauss in=left device=xw mask1=500 mask2=600
echo " "
echo " *************************************************** "
echo " *************************************************** "
echo " *** Make a note of Y-posn. of peak at LEFT (above)"
echo " "
fitgauss in=right device=xw mask1=500 mask2=600
echo " "
echo " *************************************************** "
echo " *************************************************** "
echo " *** Make a note of Y-posn. of peak at RIGHT (above) "
echo " "
echo " "
echo " *************************************************** "
echo " *************************************************** "
echo " "
echo " Finally..."
echo " "
echo -n " Enter y-value at LEFT of spectrum: "
set ya = ( $< )
echo -n " Enter y-value at RIGHT of spectrum: "
set yb = ( $< )
echo " "
echo " x1 - x2 = 800 (columns 100 and 900 extracted) "
set angle = `calc exp="atand(($ya-$yb)/800)"`
echo " "
echo " >>>> ANGLE OF CONTINUUM = " $angle
echo " "
echo " "