#!/usr/bin/perl
#
# create some graphs, generating Postscript and then converting to GIF
#
# MWR 3/16/2006



# straightforward graph of several passbands on a single graph
$cmdfile = "bessel_bandpass.gnu";
$psfile = "bessel_bandpass.ps";
open(CMDFILE, ">$cmdfile") || die("can't open $cmdfile");
printf CMDFILE "set output '$psfile' \n";
printf CMDFILE "set term postscript enhanced color \"10x20\" 14 \n";
printf CMDFILE "set grid \n";
printf CMDFILE "set nokey  \n";
printf CMDFILE "set xlabel 'Wavelength (Angstroms)' \n";
printf CMDFILE "set ylabel 'Overall transmission'  \n";
printf CMDFILE "set title 'The Bessell approximations to UBVRI passbands'  \n";
printf CMDFILE "set label 1 'U' at 3400,0.4 \n";
printf CMDFILE "set label 2 'B' at 4400,0.4 \n";
printf CMDFILE "set label 3 'V' at 5300,0.4 \n";
printf CMDFILE "set label 4 'R' at 6400,0.4 \n";
printf CMDFILE "set label 5 'I' at 8100,0.4 \n";
printf CMDFILE "plot [3000:10000][-0.05:1.1]";
printf CMDFILE " 'bess-u.pass' using (\$1):(\$2) ";
printf CMDFILE "     with lines t '' lw 4  ";
printf CMDFILE " , ";
printf CMDFILE " 'bess-b.pass' using (\$1):(\$2) ";
printf CMDFILE "     with lines t '' lw 4  ";
printf CMDFILE " , ";
printf CMDFILE " 'bess-v.pass' using (\$1):(\$2) ";
printf CMDFILE "     with lines t '' lw 4  ";
printf CMDFILE " , ";
printf CMDFILE " 'bess-r.pass' using (\$1):(\$2) ";
printf CMDFILE "     with lines t '' lw 4  ";
printf CMDFILE " , ";
printf CMDFILE " 'bess-i.pass' using (\$1):(\$2) ";
printf CMDFILE "     with lines t '' lw 4  ";
printf CMDFILE " \n";

printf CMDFILE "quit \n";
close(CMDFILE);

$cmd = "gnuplot $cmdfile ";
printf "cmd is ..$cmd.. \n";
$ret = `$cmd`;
$giffile = $psfile;
$giffile =~ s/.ps/.gif/;
$cmd = "convert -rotate 90 $psfile $giffile";
printf "cmd is ..$cmd.. \n";
$ret = `$cmd`;


exit 0;