#!/bin/csh 
#    
#    In /home/cdavis/scripts --> cgs4_pol_quicklook.csh
#    CJD - June 2000
#
echo " "
echo "    * * * * * * * * * * * * * * * * * * * * * *"
echo "    *         IRPOL QUICK-LOOK FOR CGS4       *"
echo "    * * * * * * * * * * * * * * * * * * * * * *"

# Start the FIGARO, KAPPA & TSP routines (quietly)
alias echo 'echo > /dev/null'
figaro
kappa
tsp
unalias echo 
#
echo " "
echo -n "UT data of the observations: "
set date="$<"

echo -n "Number of first frame: "
set n="$<"
set frame1=ro${date}_${n}
@ m = ${n} + 1
set frame2=ro${date}_${m}
@ m = ${n} + 2
set frame3=ro${date}_${m}
@ m = ${n} + 3
set frame4=ro${date}_${m}
@ m = ${n} + 4
set frame5=ro${date}_${m}
@ m = ${n} + 5
set frame6=ro${date}_${m}
@ m = ${n} + 6
set frame7=ro${date}_${m}
@ m = ${n} + 7
set frame8=ro${date}_${m}

#
#
echo " "
echo "--->>> Reducing files ${n} through to ${m}..."
echo "--->>> Result written to file rgp${date}_${n}"
echo " "
#
# subtract skies from obj frames:
#
sub $frame1 $frame2 frame_0d; echo -n "."
sub $frame4 $frame3 frame_45d; echo -n "."
sub $frame5 $frame6 frame_22d; echo -n "."
sub $frame8 $frame7 frame_67d; echo -n "."
#
# rotating images 90 degress anti-clockwise
#
rotate angle=90 in=frame_0d out=frame_0dr
irevy image=frame_0dr output=frame_0rev  ; echo -n "."
rotate angle=90 in=frame_45d out=frame_45dr
irevy image=frame_45dr output=frame_45rev ; echo -n "."
rotate angle=90 in=frame_22d out=frame_22dr
irevy image=frame_22dr output=frame_22rev ; echo -n "."
rotate angle=90 in=frame_67d out=frame_67dr
irevy image=frame_67dr output=frame_67rev ; echo -n "."
echo " "
echo " "
echo " NB:  ASTART = start row for lower "OBJECT" aperture"
echo "      BSTART = start row for lower "SKY" aperture (i.e. after slide)"
echo "      ***  NO CORRECTION FOR CURVATURE/DISTORTIONS  *** "
echo " "

ccd2pol pos1=frame_0rev pos2=frame_45rev pos3=frame_22rev pos4=frame_67rev aperture=a bias=0 readnoise=1 photadu=1 algorithm=ratio output=temp 

#   Plot result
mv temp.sdf rgp${date}_${n}.sdf
pplot rgp${date}_${n} device=xwindows label=rgb${date}_${n}

echo " "
echo "--->>> Calculate Pol.(%) and P.A. over selected wavelength region " 
echo " "
ptheta INPUT=rgp${date}_${n}

rm frame*
echo " "
exit