"; if (! defined($query) || $query eq '') { print "
\n"; } else { # user hit "Submit", time to do the real job $qn{'tmp'}=''; @q = split('&', $query, 10000); foreach $_ (@q) { if ($_ ne '') { ($name,$value) = split('=', $_, 2); if (!defined($value)) { $value=''; } if ( $name eq 'tmp' ){ $qn{'tmp'}="$qn{'tmp'} $value"; } elsif($name eq 'var') { $qn{'var'}="$qn{'var'} $value"; } else { $qn{$name} = $value; # fix any funny file names: get rid of / * and ? $value =~ s/\///g; $value =~ s/\*\///g; $value =~ s/\?//g; $value =~ s/\"//g; $value =~ s/\'//g; $safeqn{$name} = $value; } } } # $model = $safeqn{'model'}; $model = "enshires"; $stnm = $safeqn{'stname'}; $stnm =~ s/%7E/~/; @stnmary=split('~',$stnm,20); $stt = $stnmary[1]; $stlat= $safeqn{'stlat'}; $stlon= $safeqn{'stlon'}; $fyy = $safeqn{'fyy'}; $fmm = $safeqn{'fmm'}; $fdd = $safeqn{'fdd'}; $fhr = $safeqn{'fhr'}; $cycle = $safeqn{'cycle'}; $ltmpopr = $safeqn{'ltmpopr'}; $htmpopr = $safeqn{'htmpopr'}; $rainopr = $safeqn{'rainopr'}; $windopr = $safeqn{'windopr'}; $lowt = $safeqn{'lowt'}; $hight = $safeqn{'hight'}; $precip= $safeqn{'precip'}; $windsp= $safeqn{'windsp'}; $tmpunit1 = $safeqn{'tmpunit1'}; $tmpunit2 = $safeqn{'tmpunit2'}; $windunit = $safeqn{'windunit'}; $urlqry = $safeqn{'urlqry'}; @qvar=split(' ', $qn{'var'}); if ( defined($qvar[0]) && $qvar[0] ne '') { if ($qvar[0] eq 'temp'){ $temp='on'; } elsif ($qvar[0] eq 'rain') { $rain='on'; } elsif ($qvar[0] eq 'wind'){ $wind='on'; } if ($qvar[1] eq 'rain') { $rain='on'; }elsif ($qvar[1] eq 'wind') { $wind='on'; } if ($qvar[2] eq 'wind') { $wind='on'; } } @qtmp=split(' ', $qn{'tmp'}); if (defined($qtmp[0]) && $qtmp[0] ne '') { if ($qtmp[0] eq 'ltmp'){ $lowtmp='on'; } elsif ($qtmp[0] eq 'htmp') { $hightmp='on'; } } if (defined($qtmp[1]) && $qtmp[1] ne '') { if ($qtmp[1] eq 'htmp') { $hightmp='on'; } } if (defined($tmpunit1) && $tmpunit1 eq "K" ){ if (defined($lowtmp) && $lowtmp eq 'on' && defined($lowt) && $lowt ne ""){ $tmin=$lowt; } } if (defined($tmpunit2) && $tmpunit2 eq "K" ){ if (defined($hightmp) && $hightmp eq 'on' && defined($hight) && $hight ne ""){ $tmax=$hight; } } if ( defined($tmpunit1) && $tmpunit1 eq "C") { if (defined($lowtmp) && $lowtmp eq 'on' && defined($lowt) && $lowt ne ""){ $tmin=$lowt+273.15; } } if ( defined($tmpunit2) && $tmpunit2 eq "C") { if (defined($hightmp) && $hightmp eq 'on' && defined($hight)&& $hight ne ""){ $tmax=$hight+273.15; } } if ( defined($tmpunit1) && $tmpunit1 eq "F" ) { if (defined($lowtmp) && $lowtmp eq 'on' && defined($lowt) && $lowt ne ""){ $tmin=($lowt-32)*5./9 +273.15; } } if ( defined($tmpunit2) && $tmpunit2 eq "F" ) { # print "in tmpunit2,"; if (defined($hightmp) && $hightmp eq 'on' && defined($hight) && $hight ne ""){ $tmax=($hight-32)*5./9 +273.15; } } if ( defined($windunit) && $windunit eq "metersec" ) { if (defined($windsp) && $windsp ne ""){ $wspd=$windsp; } } if ( defined($windunit) && $windunit eq "milehr" ) { if (defined($windsp) && $windsp ne ""){ $wspd=$windsp*0.44704; } } if ( defined($windunit) && $windunit eq "knots") { if (defined($windsp) && $windsp ne ""){ $wspd=$windsp*0.51444444; } } # print "model,stnm, stt,fyy,fmm,fdd,fhr,cycle=$model, $stnm, $stt,$fyy,$fmm,$fdd,$fhr,$cycle
\n"; # print "test: lon=$lon[$stt], lat=$lat[$stt], qn(tmp)=$qn{'tmp'},lowtmp=$lowtmp, hightmp=$hightmp, lowt=$lowt, hight=$hight, ltmpopr=$ltmpopr, htmpopr=$htmpopr, tmpunit1=$tmpunit1,tmpunit2=$tmpunit2,tmi=$tmin,tma=$tmax,
\n"; #compute lat, lon and time if ( $model eq "enshires" ) { if (! defined($stlon) || $stlon eq ''){ $stlon=$lon[$stt]; print "stlon=$stlon \n"; } if ($stlon lt 0.) { $lon1=360+$stlon; }else { $lon1=$stlon; } $lon1=int($lon1+0.5); if ( $lon1 > 359 ) { $lon1=0; } if ( $lon1 < 0 ) { $lon1=0; } # if ($flag1 eq "N" ) { if (!defined($stlat) || $stlat eq ''){ $stlat=$lat[$stt]; print "stlon=$stlon \n"; } $lat1=$stlat+90.; # } else { # $lat1=$lat[stt]; # } # $lat1=$lat1- ( ($lat1*100)%100)*0.01; $lat1=int($lat1+0.5); &time_intv(); $tm1=$newti; } # if ( $model eq "eta" ) { # $lonn=$lon[$stt]; ## system "`echo $w1 $lat[$stt] $lon[$stt] >$tmpfile1`"; # system "`$w1 $lat[$stt] $lon[$stt] >>$tmpfile1`"; # # open(TMP, "<$tmpfile1"); # while (
\n"; # $tm1=$newti; # close(TMP); # } # print "stt=$stt, lon1=$lon1, lat1=$lat1, newti=$newti
\n"; $modelcycle=$model."_".$cycle; $modela="ens"; $year=$yy+2000; if ($model eq "enshires") { $urlorg="http://nomad5.ncep.noaa.gov:9090/dods/$model/archive/$modela$year$mm$dd"; } elsif ($model eq "eta") { $urlorg="http://nomad1.ncep.noaa.gov:9090/dods/$model/archive/$model$year$mm$dd/$modelcycle"; } print "The station is: $stNM[$stt]
\n"; if ($lat[$stt] gt 0.) { print "Lat: $stlat N, "; } else { print "Lat: $stlat S, "; } if ($lon[$stt] lt 0.) { print "Lon:$stlon W
\n"; } else { print "Lon:$stlon E
\n"; } # print "initime=$initime"; $newfhr=substr($initime,8,2); $newfdd=substr($initime,6,2); $newfmm=substr($initime,4,2); $newfmmnm=monthname($newfmm); $newfyy=substr($initime,0,4); # print "FORECAST: $fhr Z, $fmmn $fdd, $year
\n"; print "FORECAST: $newfhr Z, $newfmmnm $newfdd, $newfyy
\n"; print "
\n";
# if (defined($temp) && defined($rain)) {
# print "Temperature, lowest TEMP: $ltmpopr $lowt $tmpunit ($tmin K)--- highest TEMP: $htmpopr $hight $tmpunit ($tmax K)
";
# print "Precipitation, $rainopr $precip mm/day
";
# if ( $lowtmp eq 'on' || $hightmp eq 'on'){
# &getpossibility("temperature", "precipitation");
# }
# }
if (defined($temp) ) {
print "
Event (Accumulated over the last 6-hrs):
\n";
print "Temperature, lowest TEMP: $ltmpopr $lowt $tmpunit ($tmin K)--- highest TEMP: $htmpopr $hight $tmpunit ($tmax K)
";
if ( $lowtmp eq 'on' || $hightmp eq 'on'){
&getpossibility("temperature");
}
}
if (defined($rain)) {
print "
\n";
print "Event:
\n";
print "Precipitation, $rainopr $precip mm/day
";
&getpossibility("precipitation");
}
if (defined($wind)) {
print "
\n";
print "Event:
\n";
if ( $windunit eq "metersec") {
print "Wind speed, $windopr $windsp m/s
";
} elsif ($windunit eq "milehr") {
print "Wind speed, $windopr $windsp mile/hr ($wspd m/s)";
} else {
print "Wind speed, $windopr $windsp knots ($wspd m/s)";
}
&getpossibility("windspeed");
}
# if ($model eq "gfs") {
# getvar("tmax2m", "High Temperatur", "K");
# getvar("tmin2m", "Low Temperatur", "K");
# getvar("rh2m", "Relative Humidity", "\%");
# getvar("apcp", "Precipitation", "");
# } elsif ($model eq "eta") {
# getvar("rh2m", "Relative Humidity", "\%");
# getvar("apcpsfc", "Precipitation", "");
# }
}
print "\n";
exit 0;
sub getpossibility() {
$varnm=$_[0];
# $varnm2=$_[1];
if ($varnm eq "temperature" ) {
$i=0;
foreach $member (c0,n1,n2,n3,n4,n5,p1,p2,p3,p4,p5) {
if ($urlqry eq "yes") {
print "
member=$member
";
}
if (defined($lowtmp)){
$url=$urlorg."\/"."$modela$member"."_"."$cycle"."_1x1"."\.ascii?tmin2m\[$tm1:$tm1\]\[$lat1:$lat1\]\[$lon1:$lon1\]";
if ($urlqry eq "yes") {
print "URL is: $url
\n";
}
$myvarurl=`$wwwgrab $url`;
@tmp1=split(" ",$myvarurl, 20);
$tminmem[$i]=$tmp1[3];
if ($urlqry eq "yes") {
print "tminmem=$tminmem[$i]
";
}
}
if (defined($hightmp)){
$url=$urlorg."\/"."$modela$member"."_"."$cycle"."_1x1"."\.ascii?tmax2m\[$tm1:$tm1\]\[$lat1:$lat1\]\[$lon1:$lon1\]";
if ($urlqry eq "yes") {
print "URL is: $url
\n";
}
$myvarurl=`$wwwgrab $url`;
@tmp1=split(" ",$myvarurl, 20);
$tmaxmem[$i]=$tmp1[3];
if ($urlqry eq "yes") {
print "tmaxmem=$tmaxmem[$i]
";
}
}
$i++;
}
#compute probability
$p=0;
for($i=0;$i<11;$i++) {
if ( defined($temp) && !defined($tmin) && defined($tmax) ) {
if ( $htmpopr eq "gt" && $tmaxmem[$i] > $tmax ) { $p++; }
if ( $htmpopr eq "lt" && $tmaxmem[$i] < $tmax ) { $p++; }
}
if ( defined($temp) && defined($tmin) && !defined($tmax) ) {
if ( $ltmpopr eq "gt" && $tminmem[$i] > $tmin ) { $p++; }
if ( $ltmpopr eq "lt" && $tminmem[$i] < $tmin ) { $p++; }
}
if ( defined($temp) && defined($tmin) && defined($tmax) ){
if ( $ltmpopr eq "gt" && $htmpopr eq "gt") {
if ( $tminmem[$i] > $tmin && $tmaxmem[$i] > $tmax) { $p++; }
}
if ( $ltmpopr eq "gt" && $htmpopr eq "lt") {
if ( $tminmem[$i] > $tmin && $tmaxmem[$i] < $tmax) { $p++; }
}
if ( $ltmpopr eq "lt" && $htmpopr eq "gt") {
if ( $tminmem[$i] < $tmin && $tmaxmem[$i] > $tmax) { $p++; }
}
if ( $ltmpopr eq "lt" && $htmpopr eq "lt") {
if ( $tminmem[$i] < $tmin && $tmaxmem[$i] < $tmax) { $p++; }
}
}
}
$prob=int($p/11.0*100);
print "
Event probability: $prob\%
\n"; } #for precipitation $i=0; if ($varnm eq "precipitation" ){ foreach $member (c0,n1,n2,n3,n4,n5,p1,p2,p3,p4,p5) { if ($urlqry eq "yes") { print "member=$member
"; } if (defined($precip)){ $url=$urlorg."\/"."$modela$member"."_"."$cycle"."_1x1"."\.ascii?pratesfc\[$tm1:$tm1\]\[$lat1:$lat1\]\[$lon1:$lon1\]"; if ($urlqry eq "yes") { print "URL is: $url
\n"; } $myvarurl=`$wwwgrab $url`; @rain1=split(" ",$myvarurl, 20); $rainmem[$i]=$rain1[3]*86400.0; if ($urlqry eq "yes") { print "rainmem=$rainmem[$i]
"; } } $i++; } #compute probability $p=0; for($i=0;$i<11;$i++) { if ( defined($rain) && defined($precip) ) { if ($rainopr eq "gt" && $rainmem[$i] > $precip) { $p++; } if ($rainopr eq "lt" && $rainmem[$i] < $precip) { $p++; } } } $prob=int($p/11.0*100); print "
Event probability: $prob\%
\n"; } #for wind speed $i=0; if ($varnm eq "windspeed"){ foreach $member (c0,n1,n2,n3,n4,n5,p1,p2,p3,p4,p5) { if ($urlqry eq "yes") { print "member=$member
"; } if (defined($wspd)){ $url=$urlorg."\/"."$modela$member"."_"."$cycle"."_1x1"."\.ascii?ugrd10m\[$tm1:$tm1\]\[$lat1:$lat1\]\[$lon1:$lon1\]"; if ($urlqry eq "yes") { print "URL is: $url
\n"; } $myvarurl=`$wwwgrab $url`; @uw=split(" ",$myvarurl, 20); $uwind=$uw[3]; $url=$urlorg."\/"."$modela$member"."_"."$cycle"."_1x1"."\.ascii?vgrd10m\[$tm1:$tm1\]\[$lat1:$lat1\]\[$lon1:$lon1\]"; if ($urlqry eq "yes") { print "URL is: $url
\n"; } $myvarurl=`$wwwgrab $url`; @vw=split(" ",$myvarurl, 20); $vwind=$vw[3]; $windspeedmem[$i]=sqrt($uwind*$uwind+$vwind*$vwind); $windanglemem[$i]=atan2($uwind, $vwind); if ($urlqry eq "yes") { print "windspeed=$windspeedmem[$i], windangle=$windanglemem[$i]
"; } } $i++; } #compute probability $p=0; for($i=0;$i<11;$i++) { if ( defined($wind) && defined($windsp) ) { if ($windopr eq "gt" && $windspeedmem[$i] > $windsp) { $p++; } if ($windopr eq "lt" && $windspeedmem[$i] < $windsp) { $p++; } } } $prob=int($p/11.0*100); print "
Event probability: $prob\%
\n"; # write into a binary file: `mkdir -p $winddir`; open (FH, ">$winddat") or die "can't open $winddat"; for ($i=0; $i<11;$i++) { print FH "$windspeedmem[$i]\n"; } close(FH); `cd $winddir; $runwrtbin windtmp`; open (FHCTL,">",$windctl) or die "can't open $windctl"; print FHCTL "dset $winddir/windbin\n"; print FHCTL "title wind \n"; print FHCTL "options big_endian sequential\n"; print FHCTL "undef 9.99E20\n"; print FHCTL "XDEF 1 LINEAR 1 1\n"; print FHCTL "YDEF 11 LINEAR 1.8 1\n"; print FHCTL "ZDEF 1 LEVELS 1 1\n"; print FHCTL "TDEF 1 LINEAR $newfmmnm$newfyy 1mo\n"; print FHCTL "VARS 1\n"; print FHCTL "wind 01 99 smean\n"; print FHCTL "ENDVARS\n"; close(FHCTL); open (FHGRADS, ">", $GRADSFILE); print FHGRADS "'open $windctl'\n"; print FHGRADS "'d wind'\n"; print FHGRADS "'printim $windpng x800 y600 white'\n"; close(FHGRADS); #execute grads $plotmode="-l"; $errcnt=`echo "run $GRADSFILE" | $GRADS -b $plotmode |egrep -c '(Error|missing|Metafile not currently open)'`; #show the png image print "\n"; } else { $url=$urlorg."\.ascii?$var\[$tm1:$tm1\]\[$lat1:$lat1\]\[$lon1:$lon1\]"; print "URL is: $url
\n"; } $myvarurl=`$wwwgrab $url`; @tmp1=split(" ",$myvarurl, 20); $varvalue=$tmp1[3]; # print "myvar=$myvarurl"; # $i=0; # open(FLHD, "<$tmpfile"); # while (
$varnm: $varvalue $_[2]
\n";
print "
\n"; print "
\n";
}
sub time_intv() {
$mmn=monthname($mm);
$fmmn=monthname($fmm);
$year=$yy+2000;
$fyear=$fyy+2000;
$tmpcycle = $cycle;
$tmpcycle =~ s/z|Z//g;
$initime=$year.$mm.$dd.$tmpcycle;
$initime1=$cycle.$dd.$mmn.$year;
$lasttime=$fyear.$fmm.$fdd.$fhr;
$lasttime1=$fhr."z".$fdd.$fmmn.$fyear;
# print " dd=$dd, mmn=$mmn, year=$year, initime=$initime,lastime=$lasttime
\n";
$tmp="";
if ($model eq "gfs" || $model eq "eta") {
$time1="3hr";
}
if ($model eq "enshires" ) {
$time1="6hr";
}
if ($initime eq $lasttime1) {
$newti=0;
return;
}
while( $initime lt $lasttime && $newti < 10000 ) {
$tmp=`$ADD_DATE $initime1 1 $time1`;
@tmpt = split(" ",$tmp, 20);
$initime1=$tmpt[3].$tmpt[2].$tmpt[1].$tmpt[0];
$inimm=monthnum($tmpt[1]);
$tmpt[3] =~ s/z|Z//g;
$initime=$tmpt[0].$inimm.$tmpt[2].$tmpt[3];
$newti++;
}
# print "in sub, model=$model, init=$initime, lasttime=$lasttime,newti=$newti";
}
sub monthname( ) {
$_=$_[0];
return "jan" if (/01/);
return "feb" if (/02/);
return "mar" if (/03/);
return "apr" if (/04/);
return "may" if (/05/);
return "jun" if (/06/);
return "jul" if (/07/);
return "aug" if (/08/);
return "sep" if (/09/);
return "oct" if (/10/);
return "nov" if (/11/);
return "dec" if (/12/);
}
sub monthnum() {
$_=$_[0];
return "01" if (/jan/);
return "02" if (/feb/);
return "03" if (/mar/);
return "04" if (/apr/);
return "05" if (/may/);
return "06" if (/jun/);
return "07" if (/jul/);
return "08" if (/aug/);
return "09" if (/sep/);
return "10" if (/oct/);
return "11" if (/nov/);
return "12" if (/dec/);
}