#!/bin/sh

# arguments are two files containing one decimal
# number per line. files must have same number
# of lines. first file is x-axis coords and second
# file is y axis coords. creates a log-log scatter
# plot in postscript with file name $1_vs_$2.ps.

dev=X11
extension=out
xlab=

while test "$1" != ""; do
  case $1 in
    -d) shift; dev=$1; shift;;
    -x) shift; xlab=$1; shift;;
    -y) shift; ylab=$1; shift;;
    *) break;;
  esac
done

if test "$xlab" = ""; then
    xlab=$1
fi

if test "$ylab" = ""; then
    ylab=$2
fi

case $dev in
    post) extension=ps;;
    png) extension=png;;
esac

paste $1 $2 > scatterplot_temp_data


echo > scatterplot.gnuplot
echo 'set terminal '$dev > scatterplot.gnuplot
echo 'set logscale' >> scatterplot.gnuplot
echo 'set size ratio 1 0.6,0.6' >> scatterplot.gnuplot
echo 'set nokey' >> scatterplot.gnuplot
echo "set xlabel \"$xlab\"" >> scatterplot.gnuplot
echo "set ylabel \"$ylab\"" >> scatterplot.gnuplot
echo "set output \"$1_vs_$2.${extension}\"" >> scatterplot.gnuplot
echo set xrange '['$3']'  >> scatterplot.gnuplot
echo set yrange '['$3']' >> scatterplot.gnuplot
echo plot '[x='$3'] "scatterplot_temp_data" with points pt 7, x with lines lw 2, x/10 with lines lt 0' >> scatterplot.gnuplot
echo 'pause -1  "Hit return to continue"' >> scatterplot.gnuplot


gnuplot scatterplot.gnuplot

#rm scatterplot.gnuplot
#rm scatterplot_temp_data
