#!/bin/sh
#
# Version 0.01 First betaversion 960520
# Version 0.1  Second betaversion 960611
# 
# Usage: cap_plot fit yaxis xaxis form source freq
#        fit      - # of parameters 0, 1, 2, 3, 4, or 5
#        yaxis    - tcal=tcal, flux=source flux, or with 'nog' as prefix 
#                   using the values uncorrected for the gaincurve 
#                   f.ex nogtcal
#        form     - x=x or postscript=ps
#        xaxis    - time=t or elevation=el
#        source   - source to be plotted f.ex all,3c273,...
#        freq     - the observing frequency, use the standard one-letter code
home=/usr2/oper/bin
#home=/home/fredrik/JIVE/cap/
echo "Plotting of the station.acp data" 
#
if [ $1 = 0 ] 
then
   echo "Making a plot only, no fitting done"
   sed 's/fit /#fit/' $home/cap_tcal_$4.dem  > /tmp/slush007
   sed 's/, tcal.*//' /tmp/slush007 | sed 's/, flux.*//' > /tmp/slush.cap
elif [ $1 = 1 ]
then 
   echo "Making a 1 parameter fit"
   sed 's/ + .*$//' $home/cap_tcal_$4.dem | \
      sed 's/via caltemp, p0, p1, p2, p3, p4/via caltemp, p0/' > /tmp/slush.cap
elif [ $1 = 2 ]
then 
   echo "Making a 2 parameter fit"
   sed 's/ + p2.*$//' $home/cap_tcal_$4.dem | \
      sed 's/via caltemp, p0, p1, p2, p3, p4/via caltemp, p0, p1/' > /tmp/slush.cap
elif [ $1 = 3 ]
then 
   echo "Making a 3 parameter fit"
   sed 's/ + p3.*$//' $home/cap_tcal_$4.dem | \
      sed 's/via caltemp, p0, p1, p2, p3, p4/via caltemp, p0, p1, p2/' > /tmp/slush.cap
elif [ $1 = 4 ]
then 
   echo "Making a 4 parameter fit"
   sed 's/ + p4.*$//' $home/cap_tcal_$4.dem | \
      sed 's/via caltemp, p0, p1, p2, p3, p4/via caltemp, p0, p1, p2, p3/' > /tmp/slush.cap
elif [ $1 = 5 ] 
then
   echo "Making a 5 parameter fit"
   cp $home/cap_tcal_$4.dem /tmp/slush.cap
fi 
if [ "$2" = "tcal" ]
then 
   cp /tmp/slush.cap /tmp/cap_tcal_$4.dem
   cp $home/dummy.dat /tmp/station.dat
elif [ "$2" = "flux" ]
then 
   sed 's/Tcal/Flux/' /tmp/slush.cap | sed 's/tcal/flux/' |  \
      sed 's/3:7/3:8/' | sed 's/caltemp/sflux/' > /tmp/cap_tcal_$4.dem 
   cp $home/dummy.dat /tmp/station.dat
elif [ "$2" = "nogtcal" ]
then 
   sed 's/Tcal/No gaincurve Tcal/' /tmp/slush.cap | sed 's/3:7/3:9/' \
      > /tmp/cap_tcal_$4.dem 
   cp $home/station$6.dat /tmp/station.dat
elif [ "$2" = "nogflux" ]
then 
   sed 's/Tcal/No gaincurve Flux/' /tmp/slush.cap | sed 's/tcal/flux/' |  \
      sed 's/3:7/3:10/' | sed 's/caltemp/flux/' > /tmp/cap_tcal_$4.dem 
   cp $home/station$6.dat /tmp/station.dat
fi
if [ "$3" = "t" ]
then 
   cp /tmp/cap_tcal_$4.dem /tmp/slush.dem
elif [ "$3" = "el" ]
then 
   sed 's/Day/Elevation/' /tmp/cap_tcal_$4.dem > /tmp/slush2 
   sed 's/3:7/5:7/' /tmp/slush2 | sed 's/3:9/5:9/' | sed 's/3:8/5:8/' | \
      sed 's/3:10/5:10/' > /tmp/slush.dem
fi 
if [ "$5" = "all" ] 
then
   grep -v " -1$" station.acp > /tmp/station.acp
elif [ "$5" != "all" ] 
then 
   grep -i "$5" station.acp | grep -v " -1$" > /tmp/station.acp 
fi 
gnuplot /tmp/slush.dem
rm /tmp/slush* /tmp/station.dat /tmp/cap*
echo "Done"
