#! /usr/sbin/perl
# ========================================================================
# MakePointFiles - make individual html pages for the points
# ========================================================================
# History:
# --------------------------------------------------------------------
# 03202000  aew  Initial coding - stolen from R-arcticnet

# --------------------------------------------------------------------
# set up a hash for the datatypes
# --------------------------------------------------------------------
%dtypes = (
   "ca" => "Air Temperature (C)",
   "ua" => "Air Temperature (C)",
   "up" => "Precipitation (mm)",
   "cp" => "Precipitation (mm)",
   "ud" => "Stream Discharge (cms)",
   "cd" => "Stream Discharge (cms)",
   "us" => "Snowfall (mm)"
   );

# --------------------------------------------------------------------
# read the ClosePoints.txt file into memory for fast access....
# --------------------------------------------------------------------
$FILENAME = "ClosePoints.txt";
$ok = open(CP,$FILENAME);
if (! $ok) {
   print "Unable to open file ".$FILENAME." for reading!\n\n";
   exit;
   }

$|=1;
print "Processing - Reading ClosePoints.txt: ";

# get all the close points
while ($inline = <CP>) {
  $inline = <CP>;
  chop($inline);
  @points=split(/,/,$inline);
  $nopoints=@points;
  if ($nopoints>1) {
     for ($x=0;$x<$nopoints;$x++) {
        $temp=@points[$x];
        @points[$x]=@points[0];
        @points[0]=$temp;
        for ($y=1;$y<$nopoints;$y++) {
           if (@points[$y]) {
              $ClosePoints{@points[0]}.=(@points[$y].",");
              }
           }
        }
     }
  }
close(CP);

# force each point list to be a unique set
foreach $key (keys %ClosePoints) {
#  print "Before: ".$point."\n";
   @pointentry=split(/,/,$ClosePoints{$key});
   for($i=0;$i<$#pointentry;$i++) {
      for($j=$i+1;$j<=$#pointentry;$j++) {
         if ($pointentry[$i] eq $pointentry[$j]) {
            $pointentry[$j]="";
            }
         }
      }
   $point="";
   foreach $pe (@pointentry) {
      if ($pe ne "") {
         $point .= $pe . ",";
         }
      }
   $ClosePoints{$key}=substr($point,0,-1);
#  print "After: ".$ClosePoints{$key}."\n";
   }

print "Done\n";

# --------------------------------------------------------------------
# read the pointmap.txt file into memory for fast access....
# --------------------------------------------------------------------
print "Processing - Reading pointmap.txt: ";
$FILENAME = "stations/point.txt";
$ok = open(PM,$FILENAME);
if (! $ok) {
   print "\nError 2 - Unable to open file ".$FILENAME." for reading!\n\n";
   exit;
   }

$pmlines=0;
while ($inline = <PM>) {
  (@PIDs[$pmlines],$row,$col) = split(',',$inline);
  chomp(@PIDs[$pmlines]);
  chomp($row);
  chomp($col);
  @TILEs[$pmlines]=$row."_".$col;
  $pmlines++;
  }
close(PM);

print "Done\n";


# --------------------------------------------------------------------
# filename for sites.txt
# --------------------------------------------------------------------
$FILENAME = "stations/station.txt";
$ok = open(IN,$FILENAME);
if (! $ok) {
  print "Unable to open file ".$FILENAME." for reading!\n\n";
  exit;
}

# read header away
$inline = <IN>;

# --------------------------------------------------------------------
# read the 4 timeseries files into memory for fast access....
# --------------------------------------------------------------------
$FILENAME = "timeseries/GOM_US_ATemp_TS.txt";
print "Processing - Reading ".$FILENAME.": ";
$ok = open(TS,$FILENAME);
if (! $ok) {
   print "Unable to open file ".$FILENAME." for reading!\n\n";
   exit;
   }

$ua_dataheader = <TS>;

$ua_nlines=0;
while ($inline = <TS>) {
  @ua_Lines[$ua_nlines]=$inline;
  chop($inline);
  $inline =~ s/\ +//g;
  (@ua_StationId[$ua_nlines],@ua_Year[$ua_nlines],@ua_Jan[$ua_nlines],@ua_Feb[$ua_nlines],@ua_Mar[$ua_nlines],@ua_Apr[$ua_nlines],@ua_May[$ua_nlines],@ua_Jun[$ua_nlines],@ua_Jul[$ua_nlines],@ua_Aug[$ua_nlines],@ua_Sep[$ua_nlines],@ua_Oct[$ua_nlines],@ua_Nov[$ua_nlines],@ua_Dec[$ua_nlines]) = split('\t',$inline);
  @ua_StationId[$ua_nlines] =~ s/\ +//g;
  @ua_StationId[$ua_nlines] =~ s/"//g;
  @ua_Year[$ua_nlines] =~ s/"//g;
  $ua_nlines++;
  }
close(TS);
print "Done\n";

####################

$FILENAME = "timeseries/GOM_US_Discharge_TS.txt";
print "Processing - Reading ".$FILENAME.": ";
$ok = open(TS,$FILENAME);
if (! $ok) {
   print "Unable to open file ".$FILENAME." for reading!\n\n";
   exit;
   }

$ud_dataheader = <TS>;

$ud_nlines=0;
while ($inline = <TS>) {
  @ud_Lines[$ud_nlines]=$inline;
  chop($inline);
  $inline =~ s/\ +//g;
  (@ud_StationId[$ud_nlines],@ud_Year[$ud_nlines],@ud_Jan[$ud_nlines],@ud_Feb[$ud_nlines],@ud_Mar[$ud_nlines],@ud_Apr[$ud_nlines],@ud_May[$ud_nlines],@ud_Jun[$ud_nlines],@ud_Jul[$ud_nlines],@ud_Aug[$ud_nlines],@ud_Sep[$ud_nlines],@ud_Oct[$ud_nlines],@ud_Nov[$ud_nlines],@ud_Dec[$ud_nlines]) = split('\t',$inline);
  @ud_StationId[$ud_nlines] =~ s/\ +//g;
  @ud_StationId[$ud_nlines] =~ s/"//g;
  @ud_Year[$ud_nlines] =~ s/"//g;
  $ud_nlines++;
  }
close(TS);
print "Done\n";

####################

$FILENAME = "timeseries/GOM_US_Precip_TS.txt";
print "Processing - Reading ".$FILENAME.": ";
$ok = open(TS,$FILENAME);
if (! $ok) {
   print "Unable to open file ".$FILENAME." for reading!\n\n";
   exit;
   }

$up_dataheader = <TS>;

$up_nlines=0;
while ($inline = <TS>) {
  @up_Lines[$up_nlines]=$inline;
  chop($inline);
  $inline =~ s/\ +//g;
  (@up_StationId[$up_nlines],@up_Year[$up_nlines],@up_Jan[$up_nlines],@up_Feb[$up_nlines],@up_Mar[$up_nlines],@up_Apr[$up_nlines],@up_May[$up_nlines],@up_Jun[$up_nlines],@up_Jul[$up_nlines],@up_Aug[$up_nlines],@up_Sep[$up_nlines],@up_Oct[$up_nlines],@up_Nov[$up_nlines],@up_Dec[$up_nlines]) = split('\t',$inline);
  @up_StationId[$up_nlines] =~ s/\ +//g;
  @up_StationId[$up_nlines] =~ s/"//g;
  @up_Year[$up_nlines] =~ s/"//g;
  $up_nlines++;
  }
close(TS);
print "Done\n";

####################

$FILENAME = "timeseries/GOM_US_Snowfall_TS.txt";
print "Processing - Reading ".$FILENAME.": ";
$ok = open(TS,$FILENAME);
if (! $ok) {
   print "Unable to open file ".$FILENAME." for reading!\n\n";
   exit;
   }

$us_dataheader = <TS>;

$us_nlines=0;
while ($inline = <TS>) {
  @us_Lines[$us_nlines]=$inline;
  chop($inline);
  $inline =~ s/\ +//g;
  (@us_StationId[$us_nlines],@us_Year[$us_nlines],@us_Jan[$us_nlines],@us_Feb[$us_nlines],@us_Mar[$us_nlines],@us_Apr[$us_nlines],@us_May[$us_nlines],@us_Jun[$us_nlines],@us_Jul[$us_nlines],@us_Aug[$us_nlines],@us_Sep[$us_nlines],@us_Oct[$us_nlines],@us_Nov[$us_nlines],@us_Dec[$us_nlines]) = split('\t',$inline);
  @us_StationId[$us_nlines] =~ s/\ +//g;
  @us_StationId[$us_nlines] =~ s/"//g;
  @us_Year[$us_nlines] =~ s/"//g;
  $us_nlines++;
  }
close(TS);
print "Done\n";


$|=1;
print "Processing - Creating Point Files:            ";

# --------------------------------------------------------------------
# loop through all the points, process one at a time.....
# --------------------------------------------------------------------
while ($inline = <IN>) {
  chop($inline);
  ($PointID,$Id,$Name,$SymbolFLD,$Stat_id,$Country,$GOM_ID,$DataType,$DataSource,$StateProv,$Lat,$Lon,$Elevation,$BegYr,$EndYr,$Rest) = split('\t',$inline);

# --------------------------------------------------------------------
# clean the variables, remove spaces, quotes, print the Ticker .....
# --------------------------------------------------------------------
  $PointID =~ s/"//g;
  $Name=~ s/"//g;
  $Stat_id=~ s/"//g;
  $Country=~ s/"//g;
  $GOM_ID=~ s/"//g;
  $DataType=~ s/"//g;
  $DataSource=~ s/"//g;
  $StateProv=~ s/"//g;

  $|=1;
  printf "%s%-11s","\b\b\b\b\b\b\b\b\b\b\b",$PointID;

# --------------------------------------------------------------------
# create the Site Descriptor page....
# --------------------------------------------------------------------
# open the output HTML file
  $outname = ">Points/Des_".$PointID.".txt";
  $ok = open(OUT,$outname);
  if (! $ok) {
     print "Unable to open file ".$outname." for writing\n\n";
     exit;
     }
# print copyright if necessary
#  if ($Source eq "Hydat") {
#     print OUT $copyright;
#     }

# dump the data with labels
  print OUT "Station ID: " . substr($PointID,2) . "\n";
  print OUT "Name: " . $Name . "\n";
  print OUT "DataType: " . $dtypes{substr($PointID,0,2)} . "\n";
  print OUT "DataSource: " . $DataSource . "\n";
  print OUT "StateProv: " . $StateProv . "\n";
  print OUT "Lat: " . $Lat . "\n";
  print OUT "Lon: " . $Lon . "\n";
  print OUT "Elevation: " . $Elevation . "\n";
  print OUT "BegYr: " . $BegYr . "\n";
  print OUT "EndYr: " . $EndYr . "\n";

# Close file so we can re-use the file handle
  close (OUT);

# --------------------------------------------------------------------
# create the main page....
# --------------------------------------------------------------------
# open the output HTML file
  $outname = ">Points/".$PointID.".html";
  $ok = open(OUT,$outname);
  if (! $ok) {
     print "Unable to open file ".$outname." for writing\n\n";
     exit;
     }

# print the header
  print OUT "<HTML>\n<HEAD>\n<TITLE>";
  print OUT $Name;
  print OUT "</TITLE>\n</HEAD>\n<BODY BGCOLOR=\"#ffffff\">\n";
  print OUT "<CENTER><FONT SIZE=+2><B>";
  print OUT $Name;
  print OUT "</B></FONT>\n";
  print OUT "</CENTER>\n<HR>\n";

# --------------------------------------------------------------------
# print out the thumbnail images
# --------------------------------------------------------------------
  print OUT "<CENTER>\n";
  print OUT "<TABLE><TR>";

  for($pm=0;$pm<=$pmlines;$pm++) {
     if ($PIDs[$pm] eq $PointID) {
# print $PIDs[$pm] . "  " . $PointID . "\n";
        $tileno = $TILEs[$pm];
        }
     }
     if ($tileno eq "") {
        print "\nPointID # ".$PointID." Not found in pointmap file!\n";
        }
  print OUT "<TD VALIGN=MIDDLE><IMG BORDER=0 SRC=\"../images/spacer.gif\"><BR><A HREF=\"../Tiles/Maine_".$tileno.".html\">";
  print OUT "<IMG BORDER=0 SRC=\"hipoints/".$PointID.".gif\"></A></TD>\n";

# Image map for continent view based on Index image
  print OUT "<TD VALIGN=MIDDLE>";
  print OUT "<IMG BORDER=0 SRC=\"../Tiles/hitiles/maine_".$tileno.".gif\" USEMAP=\"#INDEXA\"></A></TD>\n";

 print OUT <<header;
<MAP NAME="INDEXA">
<AREA SHAPE="RECT" COORDS="115,11,141,37" HREF="../Tiles/Maine_4_0.html">
<AREA SHAPE="RECT" COORDS="90,37,115,62" HREF="../Tiles/Maine_3_1.html">
<AREA SHAPE="RECT" COORDS="115,37,141,62" HREF="../Tiles/Maine_4_1.html">
<AREA SHAPE="RECT" COORDS="141,37,166,62" HREF="../Tiles/Maine_5_1.html">
<AREA SHAPE="RECT" COORDS="166,37,192,62" HREF="../Tiles/Maine_6_1.html">
<AREA SHAPE="RECT" COORDS="65,62,90,87" HREF="../Tiles/Maine_2_2.html">
<AREA SHAPE="RECT" COORDS="90,62,115,87" HREF="../Tiles/Maine_3_2.html">
<AREA SHAPE="RECT" COORDS="115,62,141,87" HREF="../Tiles/Maine_4_2.html">
<AREA SHAPE="RECT" COORDS="141,62,166,87" HREF="../Tiles/Maine_5_2.html">
<AREA SHAPE="RECT" COORDS="166,62,192,87" HREF="../Tiles/Maine_6_2.html">
<AREA SHAPE="RECT" COORDS="192,62,217,87" HREF="../Tiles/Maine_7_2.html">
<AREA SHAPE="RECT" COORDS="217,62,243,87" HREF="../Tiles/Maine_8_2.html">
<AREA SHAPE="RECT" COORDS="39,87,65,112" HREF="../Tiles/Maine_1_3.html">
<AREA SHAPE="RECT" COORDS="65,87,90,112" HREF="../Tiles/Maine_2_3.html">
<AREA SHAPE="RECT" COORDS="90,87,115,112" HREF="../Tiles/Maine_3_3.html">
<AREA SHAPE="RECT" COORDS="115,87,141,112" HREF="../Tiles/Maine_4_3.html">
<AREA SHAPE="RECT" COORDS="141,87,166,112" HREF="../Tiles/Maine_5_3.html">
<AREA SHAPE="RECT" COORDS="166,87,192,112" HREF="../Tiles/Maine_6_3.html">
<AREA SHAPE="RECT" COORDS="192,87,217,112" HREF="../Tiles/Maine_7_3.html">
<AREA SHAPE="RECT" COORDS="217,87,243,112" HREF="../Tiles/Maine_8_3.html">
<AREA SHAPE="RECT" COORDS="243,87,268,112" HREF="../Tiles/Maine_9_3.html">
<AREA SHAPE="RECT" COORDS="268,87,294,112" HREF="../Tiles/Maine_10_3.html">
<AREA SHAPE="RECT" COORDS="39,112,65,137" HREF="../Tiles/Maine_1_4.html">
<AREA SHAPE="RECT" COORDS="65,112,90,137" HREF="../Tiles/Maine_2_4.html">
<AREA SHAPE="RECT" COORDS="90,112,115,137" HREF="../Tiles/Maine_3_4.html">
<AREA SHAPE="RECT" COORDS="115,112,141,137" HREF="../Tiles/Maine_4_4.html">
<AREA SHAPE="RECT" COORDS="141,112,166,137" HREF="../Tiles/Maine_5_4.html">
<AREA SHAPE="RECT" COORDS="166,112,192,137" HREF="../Tiles/Maine_6_4.html">
<AREA SHAPE="RECT" COORDS="192,112,217,137" HREF="../Tiles/Maine_7_4.html">
<AREA SHAPE="RECT" COORDS="217,112,243,137" HREF="../Tiles/Maine_8_4.html">
<AREA SHAPE="RECT" COORDS="243,112,268,137" HREF="../Tiles/Maine_9_4.html">
<AREA SHAPE="RECT" COORDS="14,137,39,163" HREF="../Tiles/Maine_0_5.html">
<AREA SHAPE="RECT" COORDS="39,137,65,163" HREF="../Tiles/Maine_1_5.html">
<AREA SHAPE="RECT" COORDS="65,137,90,163" HREF="../Tiles/Maine_2_5.html">
<AREA SHAPE="RECT" COORDS="90,137,115,163" HREF="../Tiles/Maine_3_5.html">
<AREA SHAPE="RECT" COORDS="166,137,192,163" HREF="../Tiles/Maine_6_5.html">
<AREA SHAPE="RECT" COORDS="192,137,217,163" HREF="../Tiles/Maine_7_5.html">
<AREA SHAPE="RECT" COORDS="14,163,39,188" HREF="../Tiles/Maine_0_6.html">
<AREA SHAPE="RECT" COORDS="39,163,65,188" HREF="../Tiles/Maine_1_6.html">
<AREA SHAPE="RECT" COORDS="65,163,90,188" HREF="../Tiles/Maine_2_6.html">
<AREA SHAPE="RECT" COORDS="65,188,90,213" HREF="../Tiles/Maine_2_7.html">
<AREA SHAPE="RECT" COORDS="90,188,115,213" HREF="../Tiles/Maine_3_7.html">
<AREA SHAPE="DEFAULT" NOHREF>
</MAP>
header

# --------------------------------------------------------------------
# Add some flags to the right hand side
# --------------------------------------------------------------------
  print OUT "<TD><TABLE><TR><TD ALIGN=CENTER>";
  print OUT "<IMG SRC=\"../images/flags/".substr($PointID,0,1).".gif\" BORDER=0>\n";
  print OUT "<P></TD></TR>\n";

  print OUT "<TR><TD ALIGN=CENTER>";
  print OUT "<BR><BR><IMG SRC=\"../images/flags/".$StateProv.".gif\" BORDER=0>\n";
  print OUT "</TD></TR>\n";

  print OUT "</TABLE>\n";

  print OUT "</TD></TR></TABLE>\n";
  print OUT "</CENTER>\n<HR>\n";


# --------------------------------------------------------------------
# dump the data with labels
# --------------------------------------------------------------------
  print OUT "<CENTER>\n";
  print OUT "<TABLE WIDTH=70% CELLPADDING=3 CELLSPACING=3>\n";
  print OUT "<TR><TH ALIGN=RIGHT>Name:</TH><TD ALIGN=LEFT> " . $Name . "</TD>\n";
  print OUT "<TH ALIGN=RIGHT>Lattitude:</TH><TD ALIGN=LEFT> " . $Lat . "</TD></TR>\n";
  print OUT "<TR><TH ALIGN=RIGHT>Station ID:</TH><TD ALIGN=LEFT> " . $Stat_id . "</TD>\n";
  print OUT "<TH ALIGN=RIGHT>Longitude:</TH><TD ALIGN=LEFT> " . $Lon . "</TD></TR>\n";
  print OUT "<TR><TH ALIGN=RIGHT>Data Type:</TH><TD ALIGN=LEFT> " . $dtypes{substr($PointID,0,2)} . "</TD>\n";
  print OUT "<TH ALIGN=RIGHT>Elevation:</TH><TD ALIGN=LEFT> ";
  if ($Elevation) {
     print OUT  $Elevation . " (m)";
     }
  print OUT "</TD></TR>\n";
  print OUT "<TR><TH ALIGN=RIGHT>Data Source:</TH><TD ALIGN=LEFT> " . $DataSource . "</TD>\n";
  print OUT "<TH ALIGN=RIGHT>Beginning Year:</TH><TD ALIGN=LEFT> " . $BegYr . "</TD></TR>\n";
  print OUT "<TR><TH ALIGN=RIGHT>State / Province:</TH><TD ALIGN=LEFT> " . $StateProv . "</TD>\n";
  print OUT "<TH ALIGN=RIGHT>Ending Year:</TH><TD ALIGN=LEFT> " . $EndYr . "</TD></TR>\n";
  print OUT "</TABLE>\n</CENTER>\n<HR>\n";

# --------------------------------------------------------------------
# check for close sites and if there are some - print the table
# --------------------------------------------------------------------

if ($ClosePoints{$PointID}) {
   @points = split(/,/,$ClosePoints{$PointID});
   $FILENAME = "stations/station.txt";
   $ok = open(CP,$FILENAME);
   if (! $ok) {
     print "Unable to open file ".$FILENAME." for reading!\n\n";
     exit;
   }
   print OUT "<CENTER>\n";
   print OUT "<TABLE CELLPADDING=2 CELLSPACING=2 BORDER=3>\n";
   print OUT "<TR><TH COLSPAN=3><Font Size=+2>Other Close Sites";
   print OUT "</Font></TH></TR>\n";
   print OUT "<TR><TH>Data Type</TH><TH>PointID</TH>\n";
   print OUT "<TH>Name</TH></TR>\n";
   $throwaway=<CP>;
   while ($inline2 = <CP>) {
      chomp($inline2);
      ($PointID2,$ID2,$Name2,$rest2) = split('\t',$inline2);
      $PointID2 =~ s/\ +//g;
      $PointID2 =~ s/"//g;
      chomp($PointID2);
      foreach $point (@points) {
         if ($PointID2 eq $point) {
           $Name2=~ s/"//g;
           print OUT "<TR><TD ALIGN=CENTER>";
           print OUT $dtypes{substr($PointID2,0,2)}."</TD>\n<TD>";;
           print OUT "<A HREF=\"".$PointID2.".html\">".substr($PointID2,2)."</A>";
           print OUT "</TD>\n<TD>".$Name2."</TD></TR>";
           }
        }
      }
   close(CP);
   print OUT "\n</TABLE>\n";
   print OUT "<HR>\n";
   }


# --------------------------------------------------------------------
# print the download table
# --------------------------------------------------------------------
  print OUT "<CENTER>\n";
  if ( $datatype eq "ca" || $datatype eq "cd" || $datatype eq "cp") {
     print OUT "<P><FONT COLOR=RED>Due to data provider restrictions,";
     print OUT " we are not able to provide tabular data for this site";
     print OUT " at this time.<P>";
     }
  print OUT "<TABLE WIDTH=80% BORDER=0 BGCOLOR=\"#bbbbbb\">\n";
  print OUT "<TR BGCOLOR=\"#999999\"><TH COLSPAN=2>Download</TH></TR>\n";
  print OUT "<TR><TH WIDTH=50%>\n";
  print OUT "<A HREF=\"Des_".$PointID.".txt\">Site Descriptor Information</A>\n";
  print OUT "  </TH>\n<TH WIDTH=50%>\n<A HREF=\"Data_".$PointID.".txt\">\n";
  print OUT "Site Time Series Data</A>\n";
  print OUT "</TH>\n</TR>\n";
  print OUT "</TABLE>\n</CENTER>\n<BR>\n";

# --------------------------------------------------------------------
# include the plot
# --------------------------------------------------------------------
  print OUT "<HR><CENTER>\n";
  $gifname=$PointID.".gif";
  if (-e ("HydroGraphs/".$gifname)) {
     print OUT "<A HREF=\"../HydroGraphs/".$gifname."\">";
     print OUT "<IMG BORDER=0 SRC=\"../HydroGraphs/".$gifname."\"></A>\n";
     print OUT "</CENTER>\n";
     }
  else {
     print OUT "</CENTER>\n";
     }

# --------------------------------------------------------------------
# print the table of data AND the Data output file.
# --------------------------------------------------------------------

# open the data output HTML file
  $outname = ">Points/Data_".$PointID.".txt";
  $ok = open(DATA,$outname);
  if (! $ok) {
     print "Unable to open file ".$outname." for writing\n\n";
     exit;
     }

# print the data file header 
   $datatype=substr($PointID,0,2);
   if ($datatype eq "ua") { print DATA $ua_dataheader; }
   elsif ($datatype eq "ud") { print DATA $ud_dataheader; }
   elsif ($datatype eq "up") { print DATA $up_dataheader; }
   elsif ($datatype eq "us") { print DATA $us_dataheader; }
   elsif ($datatype eq "ca") { print DATA "Due to data provider restrictions, we are not able to provide\ntabular data for this site at this time."; }
   elsif ($datatype eq "cd") { print DATA "Due to data provider restrictions, we are not able to provide\ntabular data for this site at this time."; }
   elsif ($datatype eq "cp") { print DATA "Due to data provider restrictions, we are not able to provide\ntabular data for this site at this time."; }
   else { print DATA "Data not available"; }


# if us_ loop through the rest - if the ID matches then save line to file and
# make a table entry
   print OUT "<HR>\n";
   print OUT "<B>Station Data File:</B>\n<BR>\n";
   print OUT "<CENTER>\n";
   if ($datatype eq "ut" || $datatype eq "ua" || $datatype eq "up" || $datatype eq "ud" || $datatype eq "us") {
      print OUT "<TABLE BORDER=3 CELLPADDING=2 CELLSPACING=2>\n";
      print OUT "<TR><TH ROWSPAN=2>Point Id</TH>";
      print OUT "<TH ROWSPAN=2>Year</TH>";

      if ($datatype eq "ua") { 
         print OUT "    <TH COLSPAN=12>Air Temperature C</TH>"; 
         $nlines=$ua_nlines;
         }
      elsif ($datatype eq "ud") { 
         print OUT "    <TH COLSPAN=12>Discharge m<SUP>3</SUP>/s</TH>"; 
         $nlines=$ud_nlines;
         }
      elsif ($datatype eq "up") { 
         print OUT "    <TH COLSPAN=12>PPT mm</TH>"; 
         $nlines=$up_nlines;
         }
      elsif ($datatype eq "us") { 
         print OUT "    <TH COLSPAN=12>Snowfall mm</TH>"; 
         $nlines=$us_nlines;
         }
      print OUT "</TR>\n";
      print OUT "<TR>\n";
      print OUT "<TH>Jan</TH>\n";
      print OUT "<TH>Feb</TH>\n";
      print OUT "<TH>Mar</TH>\n";
      print OUT "<TH>Apr</TH>\n";
      print OUT "<TH>May</TH>\n";
      print OUT "<TH>Jun</TH>\n";
      print OUT "<TH>Jul</TH>\n";
      print OUT "<TH>Aug</TH>\n";
      print OUT "<TH>Sep</TH>\n";
      print OUT "<TH>Oct</TH>\n";
      print OUT "<TH>Nov</TH>\n";
      print OUT "<TH>Dec</TH></TR>\n";
   
      for($n=0;$n<=$nlines;$n++) {
         if ($datatype eq "ua") { 
            if ($ua_StationId[$n] eq substr($PointID,2)) { 
               print DATA $ua_Lines[$n];
               print OUT "<TR><TD>".$ua_StationId[$n]."</TD>\n";
               print OUT "<TD>".$ua_Year[$n]."</TD>";
               if ($ua_Jan[$n]) {
		   print OUT "<TD>".$ua_Jan[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ua_Feb[$n]) {
		   print OUT "<TD>".$ua_Feb[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ua_Mar[$n]) {
		   print OUT "<TD>".$ua_Mar[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ua_Apr[$n]) {
		   print OUT "<TD>".$ua_Apr[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ua_May[$n]) {
		   print OUT "<TD>".$ua_May[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ua_Jun[$n]) {
		   print OUT "<TD>".$ua_Jun[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ua_Jul[$n]) {
		   print OUT "<TD>".$ua_Jul[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ua_Aug[$n]) {
		   print OUT "<TD>".$ua_Aug[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ua_Sep[$n]) {
		   print OUT "<TD>".$ua_Sep[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ua_Oct[$n]) {
		   print OUT "<TD>".$ua_Oct[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ua_Nov[$n]) {
		   print OUT "<TD>".$ua_Nov[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ua_Dec[$n]) {
		   print OUT "<TD>".$ua_Dec[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               print OUT "</TR>\n";
               }
            }
         elsif ($datatype eq "ud") { 
            if ($ud_StationId[$n] eq substr($PointID,2)) { 
               print DATA $ud_Lines[$n];
               print OUT "<TR><TD>".$ud_StationId[$n]."</TD>\n";
               print OUT "<TD>".$ud_Year[$n]."</TD>";
               if ($ud_Jan[$n]) {
		   print OUT "<TD>".$ud_Jan[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ud_Feb[$n]) {
		   print OUT "<TD>".$ud_Feb[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ud_Mar[$n]) {
		   print OUT "<TD>".$ud_Mar[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ud_Apr[$n]) {
		   print OUT "<TD>".$ud_Apr[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ud_May[$n]) {
		   print OUT "<TD>".$ud_May[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ud_Jun[$n]) {
		   print OUT "<TD>".$ud_Jun[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ud_Jul[$n]) {
		   print OUT "<TD>".$ud_Jul[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ud_Aug[$n]) {
		   print OUT "<TD>".$ud_Aug[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ud_Sep[$n]) {
		   print OUT "<TD>".$ud_Sep[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ud_Oct[$n]) {
		   print OUT "<TD>".$ud_Oct[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ud_Nov[$n]) {
		   print OUT "<TD>".$ud_Nov[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($ud_Dec[$n]) {
		   print OUT "<TD>".$ud_Dec[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               print OUT "</TR>\n";
               }
            }
         elsif ($datatype eq "up") { 
            if ($up_StationId[$n] eq substr($PointID,2)) { 
               print DATA $up_Lines[$n];
               print OUT "<TR><TD>".$up_StationId[$n]."</TD>\n";
               print OUT "<TD>".$up_Year[$n]."</TD>";
               if ($up_Jan[$n]) {
		   print OUT "<TD>".$up_Jan[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($up_Feb[$n]) {
		   print OUT "<TD>".$up_Feb[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($up_Mar[$n]) {
		   print OUT "<TD>".$up_Mar[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($up_Apr[$n]) {
		   print OUT "<TD>".$up_Apr[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($up_May[$n]) {
		   print OUT "<TD>".$up_May[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($up_Jun[$n]) {
		   print OUT "<TD>".$up_Jun[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($up_Jul[$n]) {
		   print OUT "<TD>".$up_Jul[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($up_Aug[$n]) {
		   print OUT "<TD>".$up_Aug[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($up_Sep[$n]) {
		   print OUT "<TD>".$up_Sep[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($up_Oct[$n]) {
		   print OUT "<TD>".$up_Oct[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($up_Nov[$n]) {
		   print OUT "<TD>".$up_Nov[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($up_Dec[$n]) {
		   print OUT "<TD>".$up_Dec[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               print OUT "</TR>\n";
               }
            }
         elsif ($datatype eq "us") { 
            if ($us_StationId[$n] eq substr($PointID,2)) { 
               print DATA $us_Lines[$n];
               print OUT "<TR><TD>".$us_StationId[$n]."</TD>\n";
               print OUT "<TD>".$us_Year[$n]."</TD>";
               if ($us_Jan[$n]) {
		   print OUT "<TD>".$us_Jan[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($us_Feb[$n]) {
		   print OUT "<TD>".$us_Feb[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($us_Mar[$n]) {
		   print OUT "<TD>".$us_Mar[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($us_Apr[$n]) {
		   print OUT "<TD>".$us_Apr[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($us_May[$n]) {
		   print OUT "<TD>".$us_May[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($us_Jun[$n]) {
		   print OUT "<TD>".$us_Jun[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($us_Jul[$n]) {
		   print OUT "<TD>".$us_Jul[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($us_Aug[$n]) {
		   print OUT "<TD>".$us_Aug[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($us_Sep[$n]) {
		   print OUT "<TD>".$us_Sep[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($us_Oct[$n]) {
		   print OUT "<TD>".$us_Oct[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($us_Nov[$n]) {
		   print OUT "<TD>".$us_Nov[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               if ($us_Dec[$n]) {
		   print OUT "<TD>".$us_Dec[$n]."</TD>\n";
	       }
               else {
		   print OUT "<TD><BR></TD>\n";
               }
               print OUT "</TR>\n";
               }
            }
         }
      print OUT "</TABLE>\n</CENTER>\n";
      close(DATA);
      }
   else { print DATA "\n"; }

# print copyright if necessary
#  if ($Source eq "Hydat") {
#    print OUT "<HR>\n";
#    print OUT $cleancopy;
#    print OUT "\n";
#    }

# print navigation links
  print OUT "<HR><CENTER>\n";
  print OUT "<A HREF=../../index.html>\n";
  print OUT "<IMG SRC=../images/GOMlogo.jpg BORDER=0></A>";
  print OUT "</CENTER>\n<HR>\n";
  close (OUT);
  }

print "\n";
end;