#!/usr/local/bin/perl
#
# Demo CGI perl script for handling requests from the sample form 
# idl-form.html with the WWW Client Browser as the user I/F front-end
# and IDL as a background process.
#
# mathews@nssdc.gsfc.nasa.gov
#
# $Header: /usr/local/etc/httpd/cgi-bin/idl-form.cgi,v 1.2 1995/05/17 23:10:13 mathews Exp $
#
# Copyright 1994 Jason Mathews
# Unpublished work.
# Permission granted to use and modify this program so long as the
# copyright above is maintained, modifications are documented, and
# credit is given for any use of the program..
#
# This program assumes there is a /idl directory in the HTML document
# root with a sub-directory tmp that is a symbolic link to /tmp
# and the HTML form "mfi-wplt.html" is also located in /idl.
#
# See http://lepmfi.gsfc.nasa.gov/idl/mfikp/mfi-wplt.html
#
# Written by: Laura Smith, July 1999.
# Using Perl version 5.00503
##

# debug mode: change to DEBUG=1 for debugging
$debug = 0;

# physical location of your cgi-lib.pl
#$cgi_lib = '/www/cgi-bin/cgi-lib.pl';

#require $cgi_lib;
#&ReadParse;	# Read in all the variables set by the form
require "./CGI.pm";
$interface = new CGI;
#print $interface->header;
print $interface->header('text/html');
# print $interface->start_html('-title' => "Plot Results",  '-TEXT' => "#FFFFFF", '-LINK' => "#C0C0C0", '-VLINK' => "#008000", '-BGCOLOR' => "#000000", '-background' => "../mfi/laura/blue.JPG", '-alink' => "#00FFFF");

# temporary work files
$idl_prog = "/tmp/idl.pro";
$gif_file = "/tmp/idl_$$.gif";
$ps_file = "/tmp/idl_$$.ps";
$dat_file = "/tmp/idl_$$.dat";
$psplot = "pspt.cgi";
$ps_data = "/tmp/idlps_$$.dat";

# If we use the same image file name then the cache will keep a
# copy of it and changes will not be seen and simultaneous users
# of this script might overwrite the same file while the browser
# is displaying it, therefore we append the process id ($$) to
# the filename to make it unique. The only drawback is that we
# need to do a cleanup after the script is run, such as using a
# cron job that deletes all files older than one day as below:
# 0 2 * * * find /tmp -name idl\*.gif -mtime +1 -exec rm -f {} \;

if ($debug) {
print "DEBUG http://$ENV{'SERVER_NAME'}$ENV{'SCRIPT_NAME'}";
print "Form Input";
EOH
}

# create the plot_arr and initialize it
$str = "instr";
@in_arr;
for($i=0; $i<23; $i++) {
  $in_arr[$i] = 0;
}

#print "\n<BR>";

$in_arr[0] = 1 if $interface->param('instr0');
$in_arr[1] = 1 if $interface->param('instr1');
$in_arr[2] = 1 if $interface->param('instr2');
$in_arr[3] = 1 if $interface->param('instr3');
$in_arr[4] = 1 if $interface->param('instr4');
$in_arr[5] = 1 if $interface->param('instr5');
$in_arr[6] = 1 if $interface->param('instr6');
$in_arr[7] = 1 if $interface->param('instr7');
$in_arr[8] = 1 if $interface->param('instr8');
$in_arr[9] = 1 if $interface->param('instr9');
$in_arr[10] = 1 if $interface->param('instr10');
$in_arr[11] = 1 if $interface->param('instr11');
$in_arr[12] = 1 if $interface->param('instr12');
$in_arr[13] = 1 if $interface->param('instr13');
$in_arr[14] = 1 if $interface->param('instr14');
$in_arr[15] = 1 if $interface->param('instr15');
$in_arr[16] = 1 if $interface->param('instr16');
$in_arr[17] = 1 if $interface->param('instr17');
$in_arr[18] = 1 if $interface->param('instr18');
$in_arr[19] = 1 if $interface->param('instr19');
$in_arr[20] = 1 if $interface->param('instr20');
$in_arr[21] = 1 if $interface->param('instr21');
$in_arr[22] = 1 if $interface->param('instr22');

#for ($p=0;$p<23;$p++) {
#  print $in_arr[$p],",";
#}
#print "\n<BR>";

# determine the format of the user imputted time interval, and determine the plot interval
$startime = $interface->param('start');
$stoptime = $interface->param('stop');
$tformat = 0;
$tformat = 1 if $interface->param('time') eq 'YYYYMMDD';

($s,$min,$hr,$mday,$m,$y,$wday,$yday,$isdst) = localtime(time);

if ($y > 100) {
  $y = $y - 100 + 2000;
}
else {
  $y = $y + 1900;
}
#print "$y<BR>\n";

if ($tformat) {
  $y1 = substr($startime,0,4);
  $y2 = substr($stoptime,0,4);
  $m1 = substr($startime,4,2);
  $m2 = substr($stoptime,4,2);
  $d1 = substr($startime,6,2);
  $d2 = substr($stoptime,6,2);
#print "start time: $y1 $m1 $d1<BR>\n";
#print "end time: $y2 $m2 $d2<BR>\n";
  $t = 1;
  #check to make sure that dates are valid
  if (($y1>$y2) || ($y1<1994) || ($y2>$y)) {
    print "<BR>You entered in an invalid date range.  Click the back button to try again.<BR> \n";
  }
  elsif ((($y1==1994) && ($m1<10)) || (($y1==1994) && ($m1==10) && ($d1<=16))) {
    print "<BR>You entered in an invalid date range.  Click the back button to try again.<BR> \n";
  }
  elsif ((($y2==$y) && ($m2>$m)) || (($y2==$y) && ($m2==$m) && ($d2>$mday))) {
    print "<BR>You entered in an invalid date range.  Click the back button to try again.<BR> \n";
  }
  #junk values so that the generated idl program doesn't crash
  $doy1 = 1;
  $doy2 = 2;
}
else {
  $y1 = substr($startime,0,4);
  $y2 = substr($stoptime,0,4);
  $y2 = substr($stoptime,0,4);
  $doy1 = substr($startime,5,3);
  $doy2 = substr($stoptime,5,3);
  #check to make sure that the dates are valid
  if (($y1>$y2) || ($y1<1994) || ($y2>$y)) {
    print "<BR>You entered in an invalid date range.  Click the back button to try again.<BR> \n";
  }
  elsif ((($y1==$y) && ($doy1<319)) || (($y2==$y) && ($doy2>$yday))) {
    print "<BR>You entered in an invalid date range.  Click the back button to try again.<BR> \n";
  }
  #junk values so that the generated idl doesn't crash
  $d1 = 11;
  $d2 = 21;
  $m1 = 11;
  $m2 = 11;
  $t=0
}

# now get the start and stop MM:HH:SS
$startmhs = $interface->param('startmhs');
$stopmhs = $interface->param('stopmhs');
$mm1 = substr($startmhs,0,2);
$mm2 = substr($stopmhs,0,2);
$hh1 = substr($startmhs,2,2);
$hh2 = substr($stopmhs,2,2);
$ss1 = substr($startmhs,4,2);
$ss2 = substr($stopmhs,4,2);
# now check to make sure that the times are valid
if($mm1 > 59) {
  $mm1 = 59;
}
if($mm1 < 0) {
  $mm1 = 0;
}
if($mm2 > 59) {
  $mm2 = 59;
}
if($mm2 < 0) {
  $mm2 = 0;
}
if($hh1 > 23) {
  $hh1 = 23;
}
if($hh1 < 0) {
  $hh1 = 0;
}
if($hh2 > 23) {
  $hh2 = 23;
}
if($hh2 < 0) {
  $hh2 = 0;
}
if($ss1 > 59) {
  $ss1 = 59;
}
if($ss1 < 0) {
  $ss1 = 0;
}
if($ss2 > 59) {
  $ss2 = 59;
}
if($ss2 < 0) {
  $ss2 = 0;
}


# check to see if spike removal is necessary, default to no
$sp = 0;
$sp = 1 if $interface->param('spkrm') eq 'YES';

# check to see if Pb and Pth should be overplotted, default to no
$op = 0;
$op = 1 if $interface->param('oplt') eq 'YES';


# determine the time resolution
#time interval [1,0,0] 1 min (1), [0,1,0] = 3 sec (2), [0,0,1] = 1 hour (0)
$time = 1;
$date[0] = 0;
$date[1] = 0;
$date[2] = 0;
$var = $interface->param('time1');
if ($var eq '1_hr') {
  $date[2] = 1;
  $time = 0;
}
elsif ($var eq '3_sec') {
  $date[1] = 1;
  $time = 2;
}
else {
  $date[0] = 1;
}

#print $date[0],$date[1],$date[2];
# grab the coordinate system type
$csys = 0;
$csys = 1 if $interface->param('coord') eq 'GSE';

#@max_min = (-3000) x 46;
for($d=0;$d<46;$d++) {
  $max_min[$d] = -3000;
}

@name_arr=("B","Bx","By","Bz","Theta","Phi","V","Vx","Vy","Vz","VTheta","VPhi","Np","Vth","Psw","Psb","Pth","Beta","Va","Ma","Mms","Epsilon","Vbs");

if (!open(PROG, "> $idl_prog")) {
#  print &PrintHeader;
  print <<EOH;
  print "<BODY><H1>ERROR: Could not create IDL program</H1>";
EOH
}

print PROG <<EOH;

ON_ERROR, 1
SET_PLOT, 'Z'


!p.background = 255
!p.color = 0

EOH

print PROG "plt_arr = [ $in_arr[0]";
for ($b = 1; $b<23; $b++) {
  print PROG " , $in_arr[$b]";
}
print PROG " ]\n";

print PROG "max_min = [ $max_min[0]";
for ($c = 1; $c<46; $c++) {
  print PROG " , $max_min[$c]";
}
print PROG " ]\n";

print PROG "time_res = [ $date[0]";
for ($d = 1; $d<3; $d++) {
  print PROG " , $date[$d]";
}
print PROG " ]\n";

print PROG <<EOH;
f_name = "$gif_file"
dat_name = "$dat_file"
IF ($t) THEN interface,plt_arr,time_res,max_min,$hh1,$mm1,$ss1,$hh2,$mm2,$ss2,f_name,dat_name,$y1,$m1,$d1,$y2,$m2,$d2,oplt=$op,spkrm=$sp,zbufplot=1,datakey=0 ELSE interface,plt_arr,time_res,max_min,$hh1,$mm1,$ss1,$hh2,$mm2,$ss2,f_name,dat_name,$y1,$doy1,$y2,$doy2,oplt=$op,spkrm=$sp,zbufplot=1,datakey=0


EOH

#print PROG " TITLE='$title'\n";
#print PROG "WRITE_GIF, '$gif_file', TVRD()\n" unless $debug;
print PROG "EXIT\n";

if ($debug) {
  print '<H2>IDL program</H2><pre>';
  open(PROG, $idl_prog);
  print <PROG>;
  $idl_log = '/tmp/idl.log';
} else {
  $idl_log = '/tmp/idl.log';
}

close(PROG);


# Setup IDL environment

$idl_dir  = '/usr/local/rsi/idl';
$lcl_dir  = '/usr/local/bin';
$ENV{'IDL_DIR'}=$idl_dir;
$ENV{'IDL_PATH'}="\+/data/wind/istplib:\+$idl_dir/lib:\+/home/user/asz/idl/lib";
$ENV{'IDL_STARTUP'} = '' if $ENV{'IDL_STARTUP'};
$ENV{'IDL_MASTER_PATH'}="/home/user/asz/idl/lib/aplot/";

# Run IDL in batch with IDL program

#system "$idl_dir/bin/idl $idl_prog > $idl_log 2>&1";
system "$lcl_dir/idl $idl_prog > $idl_log 2>&1";
#system("$lcl_dir/idl $idl_prog");
#unlink $idl_prog;

if ($debug) {
  print '</pre><hr><H2>IDL output</H2><pre>';
  open(LOG, $idl_log);
  print <LOG>;
  close(LOG);
  print '</pre>';
  unlink $idl_log;
}

# Check if output exists and has a non-zero file length.

for($t=0; $t<46; $t++)
{
   $maxminarr[$t] = $max_min[$t] . ":";
}
$maxmin = @maxminarr . ":";
$lenfile = length($dat_file);
$lenmaxmin = length($maxmin);
#print "$maxmin";

if (!$debug) {
  print "<HTML><BR>\n";
  print "<HEAD>\n";
  print "<TITLE>Plot Results</TITLE><BR\n";
  print "<BR><SCRIPT LANGUAGE=\"JavaScript\">\n";
  print "<!--\n";
  print "  function getPS(junk,lfile) {\n";
  print "    ps = window.open(\"\",\"PS_version\",\"toolbar=no,width=300,height=200\");\n";
  print "    ps.document.fgcolor=\"#FFFFFF\";\n";
  print "    ps.document.write(\"Ps_version\");\n";
  print "    var junk2 = junk;\n";
  print "    var lfile2 = lfile;\n";
  print "    ps.document.write(junk2);\n";
  print "    ps_file = \"/www/cgi-mfi/pspt.cgi\";\n";
#  print "    ps.document.write(\"<FORM ACTION=\");\n";
#  print "    ps.document.write(\"ps_file METHOD=\");\n";
#  print "    ps.document.write(\"POST>\");\n";
#  print "    ps.document.write(\"<BR>\");\n";
#  print "    ps.document.write(\" <INPUT TYPE=SUBMIT VALUE=stuff>\");\n"; 
#  print "    ps.document.write(\"</FORM>\");\n";
  print "    ps.location.href = \"http://lepmfi.gsfc.nasa.gov/cgi-mfi/pspt.cgi?$ps_data+$lenfile+$lenmaxmin+$stopmhs\";\n";
  print "    ps.document.close();\n";
  print "}\n";
  print "  function getValues() {\n";
  print "    num = 19991201;\n";
  print "    yr= num/10000;\n";
  print "    mo= (num - ((num/10000)*10000))/100;\n";
  print "    dy= num - (num-(num/10000)*10000);\n";
  print "    val = window.open(\"\",\"Values\",\"toolbar=no,width=300,height=200\");\n";
  print "    val.document.fgColor=\"#FFFFFF\";\n";
  print "    val.document.write=(\"Values\");\n";
  print "    val.document.close();\n";
  print "}\n";
  print "  function getStuff() {\n";
  print "    inp = window.open(\"\",\"Values\",\"toolbar=no,width=300,height=200\");\n";
  print "    inp.document.fgColor=\"#FFFFFF\";\n";
  print "    inp.document.write=(\"Stuff\");\n";
  print "    inp.document.close();\n";
  print "}\n";
  print "//-->\n";
  print "</SCRIPT>\n";
  print "</HEAD>";
  print "<BODY TEXT=\"#FFFFFF\" LINK=\"#C0C0C0\" VLINK=\"#008000\" BGCOLOR=\"#000000\" background=\"../mfi/laura/blue.JPG\" alink=\"#00FFFF\">";
  print "<BR>\n";
  
  
  print "<H1>Wind MFI Plot</H1>\n"; 
  if (-s $gif_file) {
    print qq!<IMG ALT=" " SRC="/idl$gif_file">\n!;
  }
  else {
    unlink $gif_file;
    if (!$debug) {
      print "IDL Graphics Error";
      EOH
    }
    print "\n<BR>An error was encountered while generating the plot.\n";
  }
}


## now get the high and low ranges that user can input if desired for the instruments to be plotted
print "<FORM ACTION=\"/cgi-mfi/wplt.cgi\" METHOD=\"POST\">\n";
$st1 = "min";
$st2 = "max";
print "You may enter in a high and low range for the parameter(s) that you selected.<BR><BR>\n";
for($k=0; $k<23; $k++) {
  # check to see that instrument has been chosen and is not one that already has fixed ranges (theta, phi, 
  # vtheta, or vphi)
  if($in_arr[$k] == 1 && $k != 10 && $k != 11 && $k != 4 && $k !=5) {
    $var1 = $st1.$k;
    $var2 = $st2.$k;
    print "<B>",$name_arr[$k],"</B><BR>\n";
    print "Low:   <INPUT TYPE =\"TEXT\" NAME=$var1 maxlength=8 size=8>\n";
    print "High:  <INPUT TYPE =\"TEXT\" NAME=$var2 maxlength=8 size=8><BR><BR>\n";
  }
}

$datakey = 1;
#make an array of everything so that it can be written to a file and read in if necessary to create a ps file

for ($g=0; $g<46; $g++) 
{
  $maxminarr[$g] = $max_min[$g] . ":";
}

$lenmaxmin = length($maxminarr);
$lendatfile = length($dat_file);
#print "<BR> length of datfile name: $lendatfile\n<BR>";

@file_array = (@in_arr, @date, $tformat, $startmhs, $stopmhs, $sp, $op, $datakey, $csys, $startime, $stoptime, $dat_file, @maxminarr);

#use FileHandle
#my $fh = new FileHandle("> $ps_data");
open(FD, "> $ps_data") || die;
print FD @file_array;
#foreach $line (@file_array)
#{
#  print FD "$line\n";
#}
close(FD);

print "<INPUT TYPE=HIDDEN NAME=plot_arr1 VALUE=\"$in_arr[0]:$in_arr[1]:$in_arr[2]:$in_arr[3]:$in_arr[4]:$in_arr[5]:$in_arr[6]\">\n";
print "<INPUT TYPE=HIDDEN NAME=plot_arr2 VALUE=\"$in_arr[7]:$in_arr[8]:$in_arr[9]:$in_arr[10]:$in_arr[11]:$in_arr[12]\">\n";
print "<INPUT TYPE=HIDDEN NAME=plot_arr3 VALUE=\"$in_arr[13]:$in_arr[14]:$in_arr[15]:$in_arr[16]:$in_arr[17]:$in_arr[18]\">\n";
print "<INPUT TYPE=HIDDEN NAME=plot_arr4 VALUE=\"$in_arr[19]:$in_arr[20]:$in_arr[21]:$in_arr[22]\">\n";

print "<INPUT TYPE=HIDDEN NAME=opt VALUE=\"$op\">\n";
print "<INPUT TYPE=HIDDEN NAME=spk VALUE=\"$sp\">\n";
print "<INPUT TYPE=HIDDEN NAME=start_time VALUE=\"$startime\">\n";
print "<INPUT TYPE=HIDDEN NAME=stop_time VALUE=\"$stoptime\">\n";
print "<INPUT TYPE=HIDDEN NAME=tformat VALUE=\"$tformat\">\n";
print "<INPUT TYPE=HIDDEN NAME=date VALUE=\"$time\">\n";
print "<INPUT TYPE=HIDDEN NAME=coord VALUE=\"$csys\">\n";
print "<INPUT TYPE=HIDDEN NAME=startmhs VALUE=\"$startmhs\">\n";
print "<INPUT TYPE=HIDDEN NAME=stopmhs VALUE=\"$stopmhs\">\n";
print "<INPUT TYPE=HIDDEN NAME=datafile VALUE=\"$dat_file\">\n";
print "<INPUT TYPE=HIDDEN NAME=psfile VALUE=\"$ps_file\">\n";
print "<INPUT TYPE=HIDDEN NAME=datakey VALUE=\"$datakey\">\n";

print "<p><input TYPE=\"submit\" VALUE=\"Generate Plot\"> <input TYPE=\"reset\" VALUE=\"Reset Form\"></p>\n";
#print "</FORM>\n";

$gfile1 = "/idl$gif_file";
$in_arr = "010101";
print "<BR><BR><BR><strong><INPUT TYPE=\"button\" NAME=\"ps_ver\" VALUE=\"Download\" onClick=getPS($startime,$lendatfile)>\n";
print "a postscript version of the graph.<BR>\n";
$Bx = "Bx";
print "<BR><BR>\n";
print "<P>Enter information here to get date specific data:<BR></strong>\n";
print "Input date (YYYYMMDD): <INPUT TYPE=\"TEXT\" NAME=\"mdy\" maxlength=8 size=8 value=\"$startime\">\n";
print "&nbsp &nbsp &nbsp Input time  (HHMMSS): <INPUT TYPE=\"TEXT\" NAME=\"hms\" maxlength=6 size=6 value=\"000000\"><BR>\n";
print "<BR><INPUT TYPE=\"button\" NAME=\"val\" VALUE=\"Get Values!\" onClick=\"getValues()\">\n";
print "\n<BR><BR></P>\n";
print "</FORM>\n";
print "<HR><A HREF=\"/idl/mfikp/mfi-wplt.html\">Click here to go back</A><BR>\n";

print "<HR><ADDRESS>Laura.C.Smith.1\@gsfc.nasa.gov</ADDRESS>";
print "<H6>Last Update: January 2000.</H6></BODY></HTML>";

print $interface->end_html();

##
## end of mfi_wplt.cgi
##