#!/bin/ksh
#
# run.ffast
#
# This script invokes the perl program "mkgrph.pl".  The name "mkgrph.pl"
# is a bit misleading as the program will not actually make the graph.  It
# will however create all the data files and set up needed arguments.
# By default, this script will create output for "H20" "f2 vs Ev".  To
# create output for other compounds/graph types, either
# 1) invoke this script: run.ffast <compound> <graph type (0-4)>
#       or
# 2) change the values of variables $p1 and $p2
#
# input:  Many parameters are passed to the program "mkgrph.pl".  See comments
#         below
# output: Data files used for graphing are created.  
#         

################# Set up parameters  ###############
if [[ $1 != "" ]] 
then
    p1=$1
else
    p1="H2O";           # the default compound can be hard coded here.
fi

if [[ $2 != "" ]] 
then
    if [[ $2 = [01234] ]]
    then
        p2=$2
    else
        print "Error:  The graph type should have a value of 0, 1, 2, 3 or 4";
        exit;
    fi
else
    p2="1";            # the default graph type can be hard coded here
                       # 0=graph of F1 vs E
                       # 1=graph of f2 vs E
                       # 2=graph of mu vs E
                       # 3=graph of Total Attenuation Coefficient vs E
                       # 4=display all table data
fi

p3=tmp/ffast.bat         # name of batch file to be used by the graphing program
p4=tmp/ffast_1.dat       # name of data file
p5=tmp/ffast_2.dat       
p6=tmp/ffast.ps
p7=0.001                 # x_start
p8=1000                  # y_end
p9="$y_start,unused"  # the web program that invokes mkgrph.pl determines
p10="$y_end,unused"   # values for y_start/y_end.  Currently the graphing
                                 # program does not use those values
p11=1.845                        # rho value as found in header files.
p12=0                            # user entered value for density
p13=0                            # 0=ASCII table; 1=HTML table


################# Invoke the program mkgrph.pl ###########
mkgrph.pl $p1 $p2 $p3 $p4 $p5 $p6 $p7 $p8 $p9 $p10 $p11 $p12 $p13

print ""
print "View files $p4 and $p5 for output"