#!/bin/sh
#
# Version 0.01 First betaversion 960514
# Version 0.10 First prerelease version 960611
# 
# Usage: cap_format frequency log-file
#        frequency - one letter code
#                    6cm=c
#                   18cm=l
#                  1.3cm=k
#                  3.6cm=x
#                   13cm=s
#                  0.7cm=q
#        log-file - name of the logfile containing the acp results 
#
echo "Extraction of relevant scr/cal info from: " $2
freq=$1
home=/usr2/oper/bin/
#home=/home/fredrik/JIVE/cap/
station=station
dat=dat
source=source
#
# Extraction of relevant lines from the log 
#
grep '#onoff#src/cal' $2 > /tmp/sl.fil2 
sed 's/\#.*src\/cal//' /tmp/sl.fil2 > /tmp/sl.fil
#
# Some very basic formatting, the output is ordered as: 
#
# source source_flux ddd.ddd az el src/cal tcal srcflux tcal/tcalaverage 
#    nogaincurve_tcal nogaincurve_srcflux flag 
#
sed 's/\([0-9][0-9]\)\([0-9][0-9][0-9]\)/\2 /' /tmp/sl.fil | \
 sed 's/\([0-9][0-9]\)\([0-9][0-9]\)\([0-9][0-9]\)\([0-9][0-9]*\)/\1 \2 \3/' \
 > /tmp/sl2.fil
sed 's/   / /g' /tmp/sl2.fil | sed 's/  / /g' > /tmp/sl.fil3
#
echo "{i++" > /tmp/sl100
#
sed 's/source_//' $home$source$freq.dat | sed 's/=.*[0-9]*//' > /tmp/source.dat
awk '{print "if ( ftr ~ /"$1"/) src[i]=source_"$1}' /tmp/source.dat | sed 's/ftr/\$5/' > /tmp/sl01
echo "sq[i]=\$7*\$7" > /tmp/sl4
echo "cu[i]=\$7*\$7*\$7" > /tmp/sl5
echo "qu[i]=\$7*\$7*\$7*\$7" > /tmp/sl6
echo "flux[i]=dpfu*(p0+p1*\$7+p2*sq[i]+p3*cu[i]+p4*qu[i])" > /tmp/sl7
#
cat /tmp/sl01 /tmp/sl4 /tmp/sl5 /tmp/sl6 /tmp/sl7 > /tmp/sl102
#
echo "d[i]=\$1+(\$2*3600+\$3*60+\$4)/(3600*24)" > /tmp/sl071
echo "tcal[i]=src[i]*flux[i]/\$8" > /tmp/sl072
echo "srcflux[i]=\$8*caltemp/flux[i]" > /tmp/sl073
echo "name[i]=\$5" > /tmp/sl074
echo "az[i]=\$6" > /tmp/sl075
echo "el[i]=\$7" > /tmp/sl076
echo "srccal[i]=\$8" > /tmp/sl077
echo "tcalsum+=tcal[i]" > /tmp/sl078
echo "nogtcal[i]=src[i]*dpfu/\$8" > /tmp/sl079
echo "nogsrcflux[i]=\$8*dpfu*caltemp}" > /tmp/sl0791
#
cat /tmp/sl071 /tmp/sl072 /tmp/sl073 /tmp/sl074 /tmp/sl075 /tmp/sl076 /tmp/sl077 /tmp/sl078 /tmp/sl079 /tmp/sl0791 > /tmp/sl103 
#
echo 'END{for (ii=1; ii<i+1; ii++) printf ftn %9s %.1f %.4f %.1f %.1f %.3f %.2f %.5f %.2f %.5f %1i\nftn,name[ii], src[ii], d[ii], az[ii], el[ii], srccal[ii], tcal[ii], srcflux[ii], nogtcal[ii], nogsrcflux[ii], 1}' | sed 's/ftn/"/g' > /tmp/sl104
cat /tmp/sl100 $home$station$freq.dat $home$source$freq.dat /tmp/sl102 \
   /tmp/sl103 /tmp/sl104 > /tmp/sl.awk  
awk -f /tmp/sl.awk /tmp/sl.fil3 | grep -v " -[0-9][0-9]*.[0-9]*" | \
   sed 's/   / /' | sed 's/  / /' | sed 's/  / /' > station.acp 
#
rm /tmp/sl*
echo "Done"
