#!/usr/local/bin/perl
#--------------------
#
# Program Name:archieve24.cgi
#
# Created by: Eric J. Horst
# Millersville University
#
#
# Activation Date: June 1998
# Modified April 1999
# Copyright 1998; all rights reserved.
#--------------------
#
# This program is designed to take weather data collected
# via RainWise equipment and archeive last 24-hours of hourly
# data. To be run by cron each day at midnight.
#
#--------------------IMPORTANT---------------------
# Two modifications need to be performed each year:
# 1.) create new year directory in Archive (ie. mkdir 2000)
# 2.) create new monthly subdirectories (ie. Jan, Feb...) in 2000 dir
#
#------
# Handle the web interface
#------
# Print out a content-type for HTTP/1.0 compatibility
print "Content-type: text/html\n\n";
# Get the input
read(STDIN, $buffer, $ENV{'CONTENT_LENGTH'});
# Split the name-value pairs
@pairs = split(/&/, $buffer);
foreach $pair (@pairs)
{
($name, $value) = split(/=/, $pair);
# Un-Webify plus signs and %-encoding
$value =~ tr/+/ /;
$value =~ s/%([a-fA-F0-9][a-fA-F0-9])/pack("C", hex($1))/eg;
$FORM{$name} = $value;
}
#----- Here's where the fun begins -----#
# Variable initialization and program setup
@DATA = `tail -3000 /var/data/tower/weather.dat`;
$obs48 = 1;
$switch = "on";
$datamiss = "false";
$base = "";
# test variables
$dew_point = 1;
$wind_chill = 2;
chop($mospell = `date '+%b'`);
chop($month = `date '+%m'`);
chop($year = `date '+%Y'`);
chop($today = `date '+%d'`);
# OUTFILE names below need changed once per year
# ALSO, the directory in which the CGI runs must be increased for new year
if ($today == 1) {
if ($month == 1) {
$outfile = "31Dec98.html";
}
elsif ($month == 2) {
$outfile = "31Jan98.html";
}
elsif ($month == 3) {
$outfile = "28Feb98.html";
}
elsif ($month == 4) {
$outfile = "31Mar98.html";
}
elsif ($month == 5) {
$outfile = "30Apr98.html";
}
elsif ($month == 6) {
$outfile = "31May98.html";
}
elsif ($month == 7) {
$outfile = "30Jun98.html";
}
elsif ($month == 8) {
$outfile = "31Jul98.html";
}
elsif ($month == 9) {
$outfile = "31Aug98.html";
}
elsif ($month == 10) {
$outfile = "30Sep98.html";
}
elsif ($month == 11) {
$outfile = "31Oct98.html";
}
else {
$outfile = "30Nov98.html";
}
}
else {
$yday = $today - 1;
chop($outfile = `date '+$yday%b%y.html'`); # Setup output file
}
open(OUTPUT, "> archive/$year/$mospell/$outfile"); # output to proper dir
#------
# Print header info
#------
printf OUTPUT "\nHourly Observations at
MU\n";
printf OUTPUT "\n\n";
printf OUTPUT "\nM.U. Hourly Observations file:
\n\n";
printf OUTPUT "$outfile";
printf OUTPUT "\n";
printf OUTPUT "\n";
printf OUTPUT
"
| Date | Time | Temp.
(oF) | Dew
Point
(oF) | Rel.
Hum.
(%) | Press.
(in.) | Wind
Dir.
(deg.) | Wind
Speed
(mph) | Hourly
Rainfall
(i
n.) | \n
\n\n";
@REVDATA = reverse(@DATA); # Reverse the order of the array
$counter = `date '+%H'`; # Start my counter based on current hour
chop($counter);
foreach $element (@REVDATA) {
if ($counter == "00") { # Keep with the 24 hour clock
$counter = 24;
}
if ($obs48 == "26") { # Total number of obs to take + 2
last;
}
($date2, $time2, $outtemp, $barometer, $winddir, $windspd, $rain, $rh) = split(/,/, $element);
($hour, $min) = split(/:/, $time2);
chop($min);
$hour2 = substr($hour,1);
#print "$hour2 $counter\n\n\n";
if ($hour2 == $counter) {
if ($min eq "00") { # Pick out only the hourly data; ignore the rest
&process_data();
$counter--;
$obs48++;
} # End of hourly data loop
} # End of hour=counter loop
else {
while ($counter - $hour2 >= 2) {
$datamiss = "true";
&process_data();
$counter--;
if ($counter == "00") {
$counter = 24;
}
$obs48++;
}
if ($min ge "55") {
&process_data();
$counter--;
if ($counter == "00") {
$counter = 24;
}
$obs48++;
}
else {
if ($datamiss eq "true") {
$base = "invalid";
&process_data();
}
else {
$datamiss = "true";
&process_data();
}
$counter--;
if ($counter == "00") {
$counter = 24;
}
$obs48++;
}
} # End of missing data line
} # End of foreach loop
printf OUTPUT "
\n\n\n";
printf OUTPUT "\n";
close OUTPUT;
`chmod 744 archive/$year/$mospell/$outfile`;
sub process_data {
if ($base ne "miss") {
if ($switch eq "off") {
($element_date2, $element_time2, $element_outtemp, $element_barometer, $element_winddir, $element_windspd, $element_rain, $element_rh) = split(/,/, $element);
($base_date2, $base_time2, $base_outtemp, $base_barometer, $base_winddir, $base_windspd, $base_rain, $base_rh) = split(/,/, $base);
if ($datamiss eq "false") {
$hrrain = $base_rain - $element_rain;
}
else {
$hrrain = "n/a";
}
# calculations for dewpoint
$constant = 461.51;
$temp_kel = (($base_outtemp - 32) * 5/9)+273.15;
$lat_heat = (2.5019 - (.00254 * ($temp_kel - 273.15))) * 10**6;
$terma_value = $temp_kel * $constant/$lat_heat;
$termb_value = -1 * (log($base_rh/100));
$dew_kel = $temp_kel/(($terma_value * $termb_value) + 1);
$dew_point = int(($dew_kel - 273.15) * 9/5 + 32);
# calculations for wind chill
# if($base_windspd < 4)
# {
# $K_value = 4096;
# }
# elsif($base_windspd >41)
# {
# $K_value = 6480;
# }
# else
# {
# $K_value = 2670.13458+(1047.3735*(log($base_windspd)));
# }
# $denom = (91.4 - $base_outtemp) * $K_value;
# $wind_chill = int(91.4 - ($denom/4096));
&print_data();
$base_obs48 = $obs48;
if ($datamiss eq "true") {
$base = "invalid";
$switch = on;
}
else {
$base = $element;
}
}
else { # Switch is "ON"
if ($datamiss eq "true") {
if ($base eq "invalid") {
&missing_data();
$base_obs48 = $obs48;
$datamiss = "false";
$base = $element;
$switch = off;
}
else {
$base = "invalid";
$base_obs48 = $obs48;
}
}
elsif ($datamiss eq "false" && $base eq "invalid") {
&missing_data();
$base_obs48 = $obs48;
$base = $element;
$switch = off;
}
else {
$base_obs48 = $obs48;
$base = $element;
$switch = off;
}
}
} # End of $base ne miss
else { # $base=miss
&missing_data();
$switch = on;
}
} # End of process_data subroutine
sub print_data {
chop($base_rh); # Get rid of annoying ^M
chop($base_rh); # Get rid of annoying ^N
chop($time = substr($base_time2,1)); # Get rid of "" on the time
chop($date = substr($base_date2,1)); # Get rid of "" on the date
printf OUTPUT "\n";
printf OUTPUT "| %s | %s | %s | %s | %s | %s | %2.2f | %s | %s | %2.2f | \n
\n\n",
$base_obs48, $date, $time, $base_outtemp, $dew_point, $base_rh, $base_barometer, $base_winddir, $base_windspd, $hrrain;
} # End of print_data subroutine
sub missing_data {
printf OUTPUT "\n";
printf OUTPUT "| $base_obs48 | miss | miss | miss | miss | miss | miss | miss | miss
| \n
\n\n";
} # End of missing_data subroutine