#!/usr/bin/perl -w # 2013 public domain Wesley Ebisuzaki # # main guts for g2subset # this version has support for wgrib 1.9.4 # subsetting by region, point and regridding # sub g2sub_main { #(local definitions)----------------------------------------- open (LOAD, "< /proc/loadavg"); $loadavg=; close(LOAD); $loadavg =~ s/ .*$//; if ($loadavg > 19) { print "Status: 503 server is busy\n"; print "Content-TYPE: text/html\n\n"; print " Error\n"; print "

Error

\n"; print "

Sorry, machine is overloaded ($loadavg), please come back later.
\n"; print "\n"; exit 0; } #(webserver definitions)------------------------------------- $wgrib2="/home/wd23ja/bin/wgrib2.1.9.4b2"; $logfile="/home/wd23ja/logs/g2subset.log"; $ENV{'PATH'} = "/bin:/usr/bin:/home/wd23ja/bin"; $ENV{IFS} = ''; #(code)----------------------------------------------------- $dir=$_[0]; $version="1.3.6.development"; $script_name=$0; $script_name =~ s/.*\///; $script_url = "$ENV{'SERVER_NAME'}$ENV{'SCRIPT_NAME'}"; $pid = getppid; $query=$ENV{'QUERY_STRING'}; if (defined($ENV{'REMOTE_HOST'})) { $REMOTE_HOST=$ENV{'REMOTE_HOST'}; } else { $REMOTE_HOST=$ENV{'REMOTE_ADDR'}; } # auto-flush stdout $old_handle = select (STDOUT); $| = 1; # perform flush after each write to STDOUT select ($old_handle); if (! -x $wgrib2) { &error_page("Missing wgrib2 ",'Network Problem #2'); exit 0;} &read_g2subrc("$dir"); @q = split('&', $query, 10000); $qn{'file'} = ''; $qn{'dir'} = ''; foreach $_ (@q) { if ($_ ne '') { ($name,$value) = split('=', $_, 2); if (!defined($value)) { $value=''; } $qn{$name} = $value; # remove any funny characters: ? " ' * / $value =~ s/[?"'*\/]//g; $safeqn{$name} = $value; } } # file name: a-z A-Z 0-1 . _ + if (defined($qn{'file'})) { $qn{'file'} =~ s/[^-a-zA-Z0-9._+]//g; } # make sure subdirectory parameter is ok: a-z A-Z 0-9 ._ + . $subdir=''; if (defined $qn{'dir'}) { $subdir = $qn{'dir'}; } if ($subdir ne '') { $subdir =~ s/%2[fF]/\//g; $subdir =~ s/[^-a-zA-Z0-9\/+_.]//g; $subdir =~ s/\/\.\.*/\//g; if (substr($subdir,0,1) ne '/') { &error_page("illegal directory: $subdir"); exit 7; } $_=$dir; foreach $t (@dirs) { $_="$_/$t"; if (! -d $_ || ! -x $_ || ! -r $_) { &error_page("illegal directory: $subdir"); exit 6; } } } $url_subdir=$subdir; $url_subdir =~ s/\//%2F/g; # read local configuration file &read_g2subrc("$dir/$subdir"); # in case not defined in local definition section if (!defined($ncol)) { $ncol=4; } if (!defined($forecast)) { $forecast=''; } if (!defined($months)) { $months=''; } if (!defined($days)) { $days=''; } if (!defined($hours)) { $hours=''; } if (!defined($gribfilter)) { $gribfilter=''; } if (!defined($nosubregion)) { $nosubregion=''; } if (!defined($nopoints)) { $nopoints=''; } if (!defined($noregrid)) { $noregrid=''; } if (!defined($var)) { $var = ''; } if (!defined($lev)) { $lev = ''; } if (!defined($REMOTE_HOST)) { $REMOTE_HOST='none'; } $wgrib_flag = ''; $date_code='d='; if (defined($forecast) && $forecast ne '') { $wgrib_flag = '-verf'; $date_code='vt='; } if ($qn{'file'} eq '') { &page1; } else { &page2 } exit 0; # ------------- page 1 ----------------- sub page1 { if (! -d "$dir$subdir" ) { &error_page('Network',"Problem: $subdir"); exit 5; } &page1_html_start($title); if ($subdir ne '') { $_=$subdir; chomp; print "Directory: $subdir
"; } open (DIRS, "cd $dir$subdir ; find . -type d -name '[a-zA-Z0-9]*' -maxdepth 1 -perm +5 -print | sort -n |"); $dir_count=0; while () { chomp; if ($dir_count++ == 0) { $col = int(60 / length($_)); if ($col eq 0) { $col=1; } print "
\nSubdirectory\n

\n"; } $_ =~ s/^\.\///; print "\n"; if ($dir_count % $col eq 0) { print "\n"; } } print ""; close(DIRS); if ($dir_count > 0) { print "
$_"; # if (-r "$dir$subdir/$_/.info") { # open(TMP,"< $dir$subdir/$_/.info"); # $_ = ; # print " $_"; # close(TMP); # } print "
"; } if (defined($files_pat)) { open (Files, "cd $dir$subdir ; find . -type f -maxdepth 1 -perm +4 -name '[a-zA-Z0-9]*' -print | egrep \''$files_pat'\' | egrep -v '\.(inv|inv-verf)\$' | sort -n |"); } else { open (Files, "cd $dir$subdir ; find . -type f -name '[a-zA-Z0-9]*' -maxdepth 1 -perm +4 -print | egrep -v '\.(inv|inv-verf)\$' | sort -n |"); } $file_count=0; $col=1; # make file listing print "\n

\n"; while () { chomp; $_ =~ s/^\.\///; $file=$_; if ($file_count == 0) { print "

\n"; print "Select a file (size in bytes)"; print "

\n"; print "\n"; } close(Files); if ($file_count == 0 && $dir_count == 0) { print "

"; print "No files or directories found
"; &credits; exit 0; } if ($file_count == 0) { print "
Select subdirectory from above list."; print "

"; &credits; exit 0; } if ($gribfilter ne '') { print "

"; print "

Step 2 (optional). Select levels/fields/times in the subset.
"; print "

Levels/Variables/Time

\n"; print "

Many times you may only want a section of a huge data\n"; print "file. Rather than transferring the entire file, this section will\n"; print "allow you to select some or all of the levels, variables, and \n"; print "dates of a GRIB2 file. The available selections will depend on the configuration.\n"; if ($levs ne "") { @_ = split(' ', $levs, 1000); print "

Select the levels desired:
"; print "\n"; print " all\n"; print " "; foreach $_ (@_) { $tmp = $_; $tmp =~ s/_/ /g; print " $tmp"; print " "; } } if ($vars ne "") { @_ = split(' ', $vars, 1000); print "

\n"; print "Select the variables desired:
"; print "\n"; print " all\n"; print " "; foreach $_ (@_) { $_ =~ s/_/ /g; print " $_"; print " "; } } if ($forecast ne "" && ($months ne "" || $days ne "" || $hours ne "")) { print "
\n"; print "
Using forecast verification time
\n"; print "\n"; } if ($months ne '') { print "\n"; print "

Select the months desired:
\n"; print "\n"; print "all\n"; print "jan\n"; print "feb\n"; print "mar\n"; print "apr\n"; print "may\n"; print "jun\n"; print "jul\n"; print "aug\n"; print "sep\n"; print "oct\n"; print "nov\n"; print "dec\n"; } if ($days ne '') { print "

Select the days of month desired:
\n"; print "all\n"; for ($i = 0; $i <= 31; $i++) { $_ = $i; if ($i <= 9) { $_ = "0$i"; } print "$i \n"; } } if ($hours ne '') { @_ = split(' ', $hours, 1000); print "

Select the hours desired:
"; print "\n"; print " all\n"; print " "; foreach $_ (@_) { $_ =~ s/_/ /g; print " $_"; print " "; } } } if ($gribfilter ne '' && ($nosubregion eq '' || $nopoints eq '' || $noregrid eq '')) { print "

Step 3 (optional). Select regional subet.
"; print "

Regional Subset: "; if ($nopoints eq '') { print "Points"; } if ($nosubregion eq '') { print ", Subregion"; } if ($noregrid eq '') { print ", Regrid"; } print "

"; print "

Regional subsetting can speed up downloads by only transferring newer grid points.
"; print "     "; print "points makes a grib2 file using nearest neighbor interpolation (within range km)
"; print "       "; print "lat-lon list = lon1:lat1:lon2:lat2:..:lonN:latN    lat < 0 for southern hemisphere
"; print "       "; print "range km = nearest neighbor has to be within range km to be used
"; print "     "; print "point (text) makes text file with nearest-neighbor grid value (unlimited range)
"; print "     "; print "subregion makes a regional subset using the original projection (cookie cutter)
"; print "     "; print "regrid interpolates to a new grid, can convert grid to earth relative winds. "; print "documentation
"; print "       "; print "interpolation bilinear interpolation is a weighted average of surrounding points
"; print "       "; print "winds the regridded winds will be earth or grid relative winds
"; print "         "; print "earth relative winds U winds are to the east and V winds are to the north
"; print "         "; print "grid relative winds U winds are from grid point (i,j) to grid point (i+1,j)
"; print "       "; print "if X=latlon then Y=Lon0:nLon:dlon and Z=Lat0:nLat:dLat
"; print "       "; print "if X=lambert/gaussian/mercator/nps/ncep see documentation
"; print "

"; print ""; print ""; print ""; print ''; print ''; print ''; print ''; print ''; print ''; print ''; print ''; print ''; print ''; print "\n"; if ($nopoints eq '') { print ""; print ""; print '"; print ''; print ''; print ''; print ''; print "\n"; print ""; print ""; print '"; print '"; print ''; print ''; print ''; print ''; print ''; print ''; print ''; print "\n"; } if ($nosubregion eq '') { print ""; print ""; print '"; print '"; print '"; print '"; print ''; print ''; print ""; } if ($noregrid eq '') { print ""; print ""; print '"; print "\n"; print '"; print "\n"; print '"; print '"; print '"; print ""; print ""; } print "
"; print 'original grid'; print "
"; print 'points'; print " '; print "lon-lat list "; print " '; print "range km "; print '
"; print 'point ('; print "text)"; print " '; print "lon "; print " '; print "lat "; print "
"; print 'subregion'; print " '; print "left lon "; print " '; print "right lon "; print " '; print "top lat "; print " '; print "bottom lat
"; print "
"; print 'regrid '; print " '; print " interpolation '; print " winds '; print " X "; print " '; print " Y "; print " '; print " Z "; print "
"; } print "

Step 4 (mandatory). Download or generate URL. The URL is needed for scripting the downloads.
"; print "

Download, show URL

"; print "

Show the URL for scripting the download or to see the wgrib2 command used to process the data,"; print "

"; print " "; print ""; print ""; &credits; exit 0; } # ------------- page 2 ----------------- sub page2 { if (defined($qn{'showurl'})) { &html_start("URL for Data Transfer: $title"); print "

Scripting file retrievals

\n"; print "

If you download many files through this web interface, you will want to "; print "automate the process. You can do this by writing a simple script with a loop around the "; print "file retrieval command."; print " This is better than clicking until your mouse is broken.\n"; # remove the showurl piece of the URL $q2=$query; $q2 =~ s/&showurl=//; # remove unneeded arguments from URL subset options if ($q2 =~ /&subset=none&/) { $q2 =~ s/&lon_lat_list=[^&]*//; $q2 =~ s/&lat1=[^&]*//; $q2 =~ s/&lon1=[^&]*//; $q2 =~ s/&km=[^&]*//; $q2 =~ s/&bottomlat=[^&]*//; $q2 =~ s/&toplat=[^&]*//; $q2 =~ s/&leftlon=[^&]*//; $q2 =~ s/&rightlon=[^&]*//; $q2 =~ s/&new_grid_X=[^&]*//; $q2 =~ s/&new_grid_Y=[^&]*//; $q2 =~ s/&new_grid_Z=[^&]*//; $q2 =~ s/&new_grid_interpol=[^&]*//; } if ($q2 =~ /&subset=points&/) { $q2 =~ s/&lat1=[^&]*//; $q2 =~ s/&lon1=[^&]*//; $q2 =~ s/&bottomlat=[^&]*//; $q2 =~ s/&toplat=[^&]*//; $q2 =~ s/&leftlon=[^&]*//; $q2 =~ s/&rightlon=[^&]*//; $q2 =~ s/&new_grid_X=[^&]*//; $q2 =~ s/&new_grid_Y=[^&]*//; $q2 =~ s/&new_grid_Z=[^&]*//; $q2 =~ s/&new_grid_interpol=[^&]*//; } if ($q2 =~ /&subset=points_text&/) { $q2 =~ s/&lon_lat_list=[^&]*//; $q2 =~ s/&km=[^&]*//; $q2 =~ s/&bottomlat=[^&]*//; $q2 =~ s/&toplat=[^&]*//; $q2 =~ s/&leftlon=[^&]*//; $q2 =~ s/&rightlon=[^&]*//; $q2 =~ s/&new_grid_X=[^&]*//; $q2 =~ s/&new_grid_Y=[^&]*//; $q2 =~ s/&new_grid_Z=[^&]*//; $q2 =~ s/&new_grid_interpol=[^&]*//; } if ($q2 =~ /&subset=subregion&/) { $q2 =~ s/&lon_lat_list=[^&]*//; $q2 =~ s/&lat1=[^&]*//; $q2 =~ s/&lon1=[^&]*//; $q2 =~ s/&km=[^&]*//; $q2 =~ s/&new_grid_X=[^&]*//; $q2 =~ s/&new_grid_Y=[^&]*//; $q2 =~ s/&new_grid_Z=[^&]*//; $q2 =~ s/&new_grid_interpol=[^&]*//; } if ($q2 =~ /&subset=regrid&/) { $q2 =~ s/&lon_lat_list=[^&]*//; $q2 =~ s/&lat1=[^&]*//; $q2 =~ s/&lon1=[^&]*//; $q2 =~ s/&km=[^&]*//; $q2 =~ s/&bottomlat=[^&]*//; $q2 =~ s/&toplat=[^&]*//; $q2 =~ s/&leftlon=[^&]*//; $q2 =~ s/&rightlon=[^&]*//; } # change & to & for display html display $q2 =~ s/&/&/g; print "

URL for downloading the file=
\nhttp://$script_url?$q2
"; print "

File retrieval command line for downloading file using the cURL command:
"; print "

curl \"http://$script_url?$q2\" -o my_file\n"; print "

File retrieval command line for downloading file using the wget command:
"; print "

wget -O my_file \"http://$script_url?$q2\"\n"; print "

information about scripting downloads."; # &credits; # exit; } # check if file exists $name = $qn{'file'}; $lname="$dir$subdir/$name"; if ( ! -f $lname) { print "Status: 404 data file not present\n"; print "Content-type: text/html\n\n"; print ''; print "data file is not present; $lname"; print "

Data file is not present: $lname

"; &credits; exit; } # generate the wgrib2 filter $filter=''; $filter2=''; if ($gribfilter ne '') { $var_filter=''; $lev_filter=''; $hor_filter=''; $day_filter=''; $mon_filter=''; foreach $_ (keys(%qn)) { $_ = &unhex($_); # get rid of nasty characters # allow ( ) in level fields tr/:;|"'\\//d; /^lev_/ && do { # convert ( -> $ and ) -> $ and ^ -> \^ $_ =~ s/$/\\(/g; $_ =~ s/$/\\)/g; $_ =~ s/\^/\\^/g; $lev_filter = join '', $lev_filter, '|', substr($_,4); }; # get rid of nasty characters tr/:;|"'()\\//d; # /^lev_/ && do {$lev_filter = join '', $lev_filter, '|', substr($_,4); /^var_/ && do {$var_filter = join '', $var_filter, '|', substr($_,4);}; /^hor_/ && do {$hor_filter = join '', $hor_filter, '|', substr($_,4,2);}; /^day_/ && do {$day_filter = join '', $day_filter, '|', substr($_,4,2);}; /^mon_/ && do {$mon_filter = join '', $mon_filter, '|', substr($_,4,2);}; } if (defined $qn{'all_var'}) { $var_filter=''; } if (defined $qn{'all_lev'}) { $lev_filter=''; } if (defined $qn{'all_day'}) { $day_filter=''; } if (defined $qn{'all_mon'}) { $mon_filter=''; } if (defined $qn{'all_hor'}) { $hor_filter=''; } if ($months eq '') { $mon_filter=''; } if ($days eq '') { $day_filter=''; } if ($hours eq '') { $hor_filter=''; } if ( $var_filter eq "|" ) { $var_filter=''; } if ( $lev_filter eq "|" ) { $lev_filter=''; } if ( $day_filter eq "|" ) { $day_filter=''; } if ( $mon_filter eq "|" ) { $mon_filter=''; } if ( $hor_filter eq "|" ) { $hor_filter=''; } if ( $var_filter ne "" ) { $filter=join '', $filter, ' -match \':(',substr($var_filter,1),'):\''; $filter2=join '', $filter2, ' | egrep \':(',substr($var_filter,1),'):\''; } if ( $lev_filter ne "" ) { $lev_filter =~ s/_/ /g; $filter=join '', $filter, ' -match \':(',substr($lev_filter,1),'):\''; $filter2=join '', $filter2, ' | egrep \':(',substr($lev_filter,1),'):\''; } if ( $hor_filter ne "" ) { $filter=join '', $filter, ' -match \':', $date_code, '........(',substr($hor_filter,1),')\''; $filter2=join '', $filter2, ' | egrep \':', $date_code, '........(',substr($hor_filter,1),')\''; } if ( $day_filter ne "" ) { $filter=join '', $filter, ' -match \':', $date_code, '......(',substr($day_filter,1),')\''; $filter2=join '', $filter2, ' | egrep \':', $date_code, '......(',substr($day_filter,1),')\''; } if ( $mon_filter ne "" ) { $filter=join '', $filter, ' -match \':', $date_code, '....(',substr($mon_filter,1),')\''; $filter2=join '', $filter2, ' | egrep \':', $date_code, '....(',substr($mon_filter,1),')\''; } if ($filter2 ne '') { $filter2=substr($filter2,3); } } # Generate wgrib2 arguements to write data to stdout $output_type = '-inv /dev/null -grib -'; $subset = ''; if ($qn{'subset'} eq 'points') { $lon_lat_list=$qn{'lon_lat_list'}; $lon_lat_list =~ s/%3[aA]/:/g; $lon_lat_list =~ s/[^-+.0-9:]//g; $km=$qn{'km'}; $km =~ s/[^-+.e0-9]//g; $output_type = "-inv /dev/null -irr_grid $lon_lat_list $km -"; $subset = 1; } if ($qn{'subset'} eq 'points_text') { $lat1=$qn{'lat1'}; $lon1=$qn{'lon1'}; $km=$qn{'km'}; $lat1 =~ s/[^-+.e0-9]//g; $lon1 =~ s/[^-+.e0-9]//g; $km =~ s/[^-+.e0-9]//g; $output_type = "-crlf -v -s -v0 -start_ft -lon $lon1 $lat1"; $subset = 1; } if (defined $qn{'subregion'} || $qn{'subset'} eq 'subregion') { $lonW=$qn{'leftlon'}; $lonE=$qn{'rightlon'}; $latN=$qn{'toplat'}; $latS=$qn{'bottomlat'}; $lonW =~ s/[^-+.e0-9]//g; $lonE =~ s/[^-+.e0-9]//g; $latN =~ s/[^-+.e0-9]//g; $latS =~ s/[^-+.e0-9]//g; if ($latN < $latS) { &error_page("North latitude < south lat"); exit; } $output_type = "-inv /dev/null -set_grib_type same -small_grib $lonW:$lonE $latS:$latN -"; $subset = 1; } if ($qn{'subset'} eq 'regrid') { $arg_interpol=''; if ($qn{'new_grid_interpol'} eq 'neighbor') { $arg_interpol='-new_grid_interpolation neighbor'; } if ($qn{'new_grid_interpol'} eq 'bicubic') { $arg_interpol='-new_grid_interpolation bicubic'; } $arg_winds='-new_grid_winds earth'; if ($qn{'new_grid_winds'} eq 'grid') { $arg_winds='-new_grid_winds grid'; } $arg_x=$safeqn{'new_grid_X'}; $arg_y=$safeqn{'new_grid_Y'}; $arg_z=$safeqn{'new_grid_Z'}; if ($arg_x eq 'ncep' && $arg_y eq 'grid') { if (! ($arg_z =~ /^t?[0-9]+$/)) { &error_page("undefined ncep grid"); exit 0;} } else { $arg_x =~ s/%3[aA]/:/g; $arg_y =~ s/%3[aA]/:/g; $arg_z =~ s/%3[aA]/:/g; $arg_x =~ s/[^-+.0-9a-z:]//g; $arg_y =~ s/[^-+.0-9e:]//g; $arg_z =~ s/[^-+.0-9e:]//g; } $output_type = "-inv /dev/null -set_grib_type same $arg_winds $arg_interpol -new_grid $arg_x $arg_y $arg_z -"; $subset = 1; } # two cases # 1: run the file through wgrib # 2: just cat the file if ($filter ne '' || $subset ne '') { # two cases # 1: no inventory file # 2: inventory file if (-r "$lname.inv$wgrib_flag") { if ($filter2 eq '') { $cmd="$wgrib2 $lname $output_type"; } else { $cmd="cat $lname.inv$wgrib_flag | $filter2 | $wgrib2 $lname -i $output_type"; } } else { $cmd="$wgrib2 $lname $wgrib_flag $filter $output_type"; } # debug # &error_page("$cmd\n"); # exit; if (defined($qn{'showurl'})) { # remove pathnames and remove trickery with stdout $cmd =~ s/-$/OUT.GRB/; $cmd =~ s/-inv \/dev\/null / /; $cmd =~ s/\/\S*\///g; print "

This web service is a wrapper for wgrib2. This is the command used to process the grib data.
$cmd

\n"; print "

"; &credits; exit; } print "Status: 200 OK\n"; print "Content-Transfer-Encoding: binary\n"; if ($qn{'subset'} eq 'points_text') { print "Content-Disposition: attachment; filename=\"$name.txt\" \n"; print "Content-Description: text file\n"; print "Content-type: text/plain \n\n"; } else { print "Content-Disposition: attachment; filename=\"$name\" \n"; print "Content-Description: grib2 file\n"; print "Content-type: application/octet-stream \n\n"; } system "$cmd 2>/dev/null"; } else { #&error_page("$cmd\n"); #exit; if (defined($qn{'showurl'})) { print "

no command is used to filter the grib data:

\n"; &credits; exit; } $nbyte = -s $lname; print "Status: 200 OK\n"; print "Content-Transfer-Encoding: binary\n"; print "Content-Disposition: attachment; filename=\"$name\" \n"; print "Content-Description: grib2 file\n"; print "Content-Length: $nbyte\n"; print "Content-type: application/octet-stream \n\n"; system "cat $lname"; } # log entry open (Logfile, ">>$logfile"); $date=`date '+%d/%b/%Y:%T'`; substr($date,-1,1)=''; print Logfile "$REMOTE_HOST - - [$date -$pid] \"good_copy:$script_name $name $filter\n"; close(Logfile); exit 0; } sub error_page { # # print out an error page (HTML) # print "Status: 500 server error\n"; &html_start("Error"); print "

Error

\n"; print "

$_[0]\n"; if (defined ( $_[1]) && $_[1] ne '') { print "

$_[1]\n"; } &credits; } sub page1_html_start { # # print out page 1 html header # &html_start("Grib2 Subset Data Transfer: $_[0]"); print "

Grib2 Subset Data Transfer: $_[0]

"; print "

Create a subset of a grib2 file for faster downloading. Designed for easy scripting.
"; print "

Step 1 (mandatory). Select a file.
"; return; } sub html_start { print "Content-type: text/html\n\n"; print ''; print "$_[0]"; return; } sub credits { print "

"; print "g2sub $version, for comments about the script: Wesley.Ebisuzaki\@noaa.gov, Jun.Wang\@noaa.gov
"; print "\n"; return; } # # unhexes a query string # sub unhex { local ($c, $i, $j, $s, $t, $len); $s=$_[0]; if (! defined($s)) { $s=''; } $t=''; $len = length($s); $i = 0; $j = 0; while ($i < $len) { $c = substr($s, $i, 1); if ($c eq '%') { $c = chr(hex(substr($s,$i+1,2))); $i += 3; } else { if ($c eq '+') { $c = ' '; } $i++; } if ($c eq '`') { $c = '"'; } if ($c eq ';') { $c = ' '; } substr($t,$j++,1) = $c; } return $t; } # # read and parse the .g2subrc file # sub read_g2subrc { local ($dir); $dir=$_[0]; $file="$dir/.g2subrc"; if (!-r "$file") { return; } open(FTP2URC, "< $file"); while () { ($name,$value) = split('=', $_, 2); if (defined $value) { chomp $value; } if ($name eq 'title') { $title=$value; } elsif ($name eq 'gribfilter') { $gribfilter=$value; } elsif ($name eq 'files_pat') { $files_pat=$value; } elsif ($name eq 'vars') { $vars=$value; } elsif ($name eq 'levs') { $levs=$value; } elsif ($name eq 'ncol') { $ncol=$value; } elsif ($name eq 'hours') { $hours=$value; } elsif ($name eq 'days') { $days=$value; } elsif ($name eq 'months') { $months=$value; } elsif ($name eq 'forecast') { $forecast=$value; } elsif ($name eq 'nosubregion') { $nosubregion=$value; } elsif ($name eq 'nopoints') { $nopoints=$value; } elsif ($name eq 'noregrid') { $noregrid=$value; } } close(FTP2URC); return; } } 1;

"; print 'original grid'; print "
"; print 'points'; print "	'; print "lon-lat list	"; print "				'; print "range km	"; print '
"; print 'point ('; print "text)"; print "	'; print "lon	"; print "	'; print "lat	"; print "
"; print 'subregion'; print "	'; print "left lon	"; print "	'; print "right lon	"; print "	'; print "top lat	"; print "	'; print "bottom lat	"; print "
"; print 'regrid '; print "	'; print " interpolation		'; print " winds		'; print " X	"; print "	'; print " Y	"; print "	'; print " Z	"; print "