bufr.pm:bufrdump.pl_source

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
bufr.pm:bufrdump.pl_source [2011-04-01 11:28:03]
pals
bufr.pm:bufrdump.pl_source [2022-05-31 09:29:31] (current)
Line 1: Line 1:
-<code perl>+<code perl bufrdump.pl>
 #!/usr/bin/perl -w #!/usr/bin/perl -w
  
Line 26: Line 26:
  
 use constant DEFAULT_TABLE_PATH => '/usr/local/lib/bufrtables'; use constant DEFAULT_TABLE_PATH => '/usr/local/lib/bufrtables';
-my $BUFRDUMP = '/metno/local/bin/bufrdump';+my $BUFRDUMP = 'bufrdump'; # You should add path if bufrdump is 
 +                           # installed in a non-standard place
  
 # Parse command line options # Parse command line options
Line 34: Line 35:
            'help',            'help',
            'tablepath=s', # Set BUFR table path            'tablepath=s', # Set BUFR table path
-           'filter=s',    # Decode observations meeting criteria in <filter file> only+           'filter=s',    # Decode observations meeting the filter criteria specified only
            'param=s',     # Decode/print specified parameters only            'param=s',     # Decode/print specified parameters only
            'csv',         # Use CSV format for printing            'csv',         # Use CSV format for printing
-           'sort',        # Sort on stationid (wmonr/call_sign/buoy_id)+           'delimiter=s', # Choose the delimiter for the CSV format 
 +           'sort',        # Sort on stationid (wmonr/nationalnr/call_sign/buoy_id/aircraft/icao_id/wigosid)
            'sort_on=s',   # Sort on specified parameter            'sort_on=s',   # Sort on specified parameter
            'station=s',   # Filter on list of stations            'station=s',   # Filter on list of stations
Line 45: Line 47:
            'lon2=i',            'lon2=i',
            'lat2=i',            'lat2=i',
 +           'obstype=s',   # Force observation type
        ) or pod2usage(-verbose => 0);        ) or pod2usage(-verbose => 0);
  
Line 69: Line 72:
 $ENV{BUFR_TABLES} .= '/' if substr($ENV{BUFR_TABLES},-1) ne '/'; $ENV{BUFR_TABLES} .= '/' if substr($ENV{BUFR_TABLES},-1) ne '/';
  
-my $filter = $option{filter} ? "--filter $option{filter}" : '';+die "Directory for BUFR tables: $ENV{BUFR_TABLES} does not exist" 
 +    if ! -d $ENV{BUFR_TABLES}; 
 + 
 +my $obstype = $option{obstype} ? "--obstype $option{obstype}" : ''; 
 +my $filt = $option{filter} ? "$option{filter}" : '';
 my $lon1 = $option{lon1} ? "--lon1 $option{lon1}" : ''; my $lon1 = $option{lon1} ? "--lon1 $option{lon1}" : '';
 my $lat1 = $option{lat1} ? "--lat1 $option{lat1}" : ''; my $lat1 = $option{lat1} ? "--lat1 $option{lat1}" : '';
 my $lon2 = $option{lon2} ? "--lon2 $option{lon2}" : ''; my $lon2 = $option{lon2} ? "--lon2 $option{lon2}" : '';
 my $lat2 = $option{lat2} ? "--lat2 $option{lat2}" : ''; my $lat2 = $option{lat2} ? "--lat2 $option{lat2}" : '';
 +my $del = $option{delimiter} ? "$option{delimiter}" : ';';
 +my $delimiter = $option{delimiter} ? "--delimiter $del" : '';
  
 # Any filter criteria provided? # Any filter criteria provided?
-my $criteria_ref = read_filter_file($option{filter});+my ($filter, $criteria_ref, $num_alt_ref) get_filter_conditions($filt);
  
 # Any specific stations requested? # Any specific stations requested?
Line 86: Line 95:
 my $csv = $option{csv} ? 1 : 0; my $csv = $option{csv} ? 1 : 0;
 # First line in CSV should be the parameters # First line in CSV should be the parameters
-print join(',', @$params_ref) . "\n" if $csv;+print join($del, @$params_ref) . "\n" if $csv;
  
 # Any transformations of units specified? # Any transformations of units specified?
Line 98: Line 107:
 # What kind of sorting is required (if any)? # What kind of sorting is required (if any)?
 ($sort_on, my $by) = get_sort_method($sort_on, $sort); ($sort_on, my $by) = get_sort_method($sort_on, $sort);
- 
  
 # Loop for processing of BUFR input files # Loop for processing of BUFR input files
Line 105: Line 113:
  
     # Dump the content of the BUFR file using the Fortran program $BUFRDUMP     # Dump the content of the BUFR file using the Fortran program $BUFRDUMP
-    my $dump `$BUFRDUMP $filter $lon1 $lat1 $lon2 $lat2 $inputfname`;+    my $fortran_options "$obstype $filter $lon1 $lat1 $lon2 $lat2"; 
 +    my $dump = `$BUFRDUMP $fortran_options $inputfname`;
     die if $?; # Reason for bufrdump failing should have been printed to STDERR     die if $?; # Reason for bufrdump failing should have been printed to STDERR
  
Line 120: Line 129:
     my %msg; # Hash with parameter name as key, parameter value as value     my %msg; # Hash with parameter name as key, parameter value as value
  
-  LINE:while (defined(my $line = shift @lines)) {+  LINE: 
 +    while (defined(my $line = shift @lines)) {
         # Each new message starts with a blank line         # Each new message starts with a blank line
         if ($line !~ /^\s*$/) {         if ($line !~ /^\s*$/) {
             # Skip error messages from libbufr, which should start with space(s)             # Skip error messages from libbufr, which should start with space(s)
-            next if $line =~ /^\s+/;+            next LINE if $line =~ /^\s+/;
             # Build up the message to be (possibly) printed             # Build up the message to be (possibly) printed
             my ($param, $value) = ($line =~ /^(.+)=\s*(.*?)\s*$/);             my ($param, $value) = ($line =~ /^(.+)=\s*(.*?)\s*$/);
 +            # Know only of one case where next check is necessary: if
 +            # a CCITT IA5 value contains new line (\n)
 +            next LINE if !defined $value;
             if ($transform_file && $transform_ref->{$param}) {             if ($transform_file && $transform_ref->{$param}) {
                 # Replace value with transformed value                 # Replace value with transformed value
Line 142: Line 155:
         if ($line =~ /^\s*$/ or @lines == 0) {         if ($line =~ /^\s*$/ or @lines == 0) {
             # A full message has been completed. Should it be printed?             # A full message has been completed. Should it be printed?
-            if ($filter && filter_obs(\%msg, $criteria_ref)) {+            if ($filt && filter_obs(\%msg, $criteria_ref, $num_alt_ref)) {
                 # Skip this message                 # Skip this message
             } elsif ($req_id && filter_station(\%msg, $req_id, $req_stn_ref)) {             } elsif ($req_id && filter_station(\%msg, $req_id, $req_stn_ref)) {
Line 154: Line 167:
                     foreach my $name (@$params_ref) {                     foreach my $name (@$params_ref) {
                         if (exists $msg{$name}) {                         if (exists $msg{$name}) {
-                            $txt .= $csv ? $msg{$name} . ',' : "$name=$msg{$name}\n";+                            $txt .= $csv ? $msg{$name} . $del : "$name=$msg{$name}\n";
                         } elsif ($forced_params_ref->{$name}) {                         } elsif ($forced_params_ref->{$name}) {
-                            $txt .= $csv ? '-32767,' : "$name=-32767\n";+                            $txt .= $csv ? '-32767' . $del : "$name=-32767\n";
                         } elsif ($csv) {                         } elsif ($csv) {
-                            $txt .= ',';+                            $txt .= $del;
                         }                         }
                     }                     }
Line 167: Line 180:
                     }                     }
                 }                 }
-                chop $txt if $csvremoves last ','+                if ($csv) { 
 +                    Remove last $del 
 +                    for (1 .. length($del)) { 
 +                        chop $txt; 
 +                    } 
 +                }
                 if ($txt) {                 if ($txt) {
                     if ($sort) {                     if ($sort) {
-                        # Sort wmonr before call signs before buoy_id+                        # Sort wmonr before nationalnr before call sign before 
 +                        # buoy_id before aircraft before icaoid before wigosid
                         if ($msg{wmonr}) {                         if ($msg{wmonr}) {
                             $stnid = '00_' . $msg{wmonr};                             $stnid = '00_' . $msg{wmonr};
 +                        } elsif ($msg{nationalnr}) {
 +                            $stnid = '10_' . $msg{nationalnr};
                         } elsif ($msg{call_sign}) {                         } elsif ($msg{call_sign}) {
-                            $stnid = '10_' . $msg{call_sign};+                            $stnid = '20_' . $msg{call_sign};
                         } elsif ($msg{buoy_id}) {                         } elsif ($msg{buoy_id}) {
-                            $stnid = '10_' . $msg{buoy_id};+                            $stnid = '30_' . $msg{buoy_id}; 
 +                        } elsif ($msg{aircraft}) { 
 +                            $stnid = '40_' . $msg{aircraft}; 
 +                        } elsif ($msg{icao_id}) { 
 +                            $stnid = '50_' . $msg{icao_id}; 
 +                        } elsif ($msg{wigosid}) { 
 +                            $stnid = '60_' . $msg{wigosid};
                         } else {                         } else {
                             # Skip observation if no station identification found                             # Skip observation if no station identification found
Line 185: Line 212:
                             my $key = $stnid . '|' . $val;                             my $key = $stnid . '|' . $val;
                             $data_of{$key} = exists $data_of{$key}                             $data_of{$key} = exists $data_of{$key}
-                                ? $data_of{$key} . "$txt \n" : "$txt \n";+                            ? $data_of{$key} . "$txt\n" : "$txt\n";
                         } else {                         } else {
                             $data_of{$stnid} = exists $data_of{$stnid}                             $data_of{$stnid} = exists $data_of{$stnid}
-                                ? $data_of{$stnid} . "$txt \n" : "$txt \n";+                            ? $data_of{$stnid} . "$txt\n" : "$txt\n";
                         }                         }
                     } elsif ($sort_on) {                     } elsif ($sort_on) {
                         my $val = exists $msg{$sort_on} ? $msg{$sort_on} : '';                         my $val = exists $msg{$sort_on} ? $msg{$sort_on} : '';
                         $data_of{$val} = exists $data_of{$val}                         $data_of{$val} = exists $data_of{$val}
-                            ? $data_of{$val} . "$txt \n" : "$txt \n";+                        ? $data_of{$val} . "$txt\n" : "$txt\n";
                     } else {                     } else {
                         # No sorting. We can print the line immediately                         # No sorting. We can print the line immediately
Line 224: Line 251:
     }     }
     print $data_of_missing_value if $data_of_missing_value;     print $data_of_missing_value if $data_of_missing_value;
 +}
 +
 +
 +# Read the filter conditions (if any). Return the filter option to be
 +# used by bufrdump, the found criteria (if any) as well as the number
 +# of succeeding alternatives for each criterium
 +sub get_filter_conditions {
 +    my $filt = shift;
 +    return ('') if ! $filt;
 +
 +    my $fortran_filter = '';
 +    my @f;
 +    if ($filt =~ /,/) {
 +        # Argument to --filter is a comma separated list
 +        @f = split /,/, $filt;
 +    } else {
 +        # Argument to --filter is a file
 +        $fortran_filter = "--filter $filt";
 +        open my $FILTER, '<', $filt
 +            or die "Cannot open $filt: $!";
 +        # Skip the criteria meant for Fortran parsing, i.e. proceed to
 +        # first line following a blank line
 +        while (<$FILTER>) {
 +            last if $_ =~ /^\s*$/;
 +        }
 +        @f = <$FILTER>;
 +        close $FILTER or die "Cannot close $filt: $!";;
 +    }
 +    return ($fortran_filter) if !@f; # BUFR descriptor criteria only
 +
 +    my @allowed_operators =
 +        ('=',
 +         '<',
 +         '<=',
 +         '>',
 +         '>=',
 +         '!=',
 +         '=~',
 +         '!~',
 +        );
 +    my @criteria;
 +    my @num_alt; # Number of alternative criteria following this,
 +                 # i.e. if line is '<cr1> | <cr2> | <cr3>' then
 +                 # corresponding values in @num_alt will be 2,1,0
 +
 +    # Read the filter criteria meant for Perl parsing, skipping blank
 +    # lines and comment lines
 +  FILTERLINE:
 +    foreach my $line (@f) {
 +        $line =~ s/^\s+//;
 +        $line =~ s/\s+$//;
 +        next FILTERLINE if !$line || $line =~ /^#/;
 +        my @crit = split /\|/, $line;
 +        my $num = scalar @crit;
 +        foreach my $criterium (@crit) {
 +            $criterium =~ s/^\s+//;
 +            $criterium =~ s/\s+$//;
 +            if ($criterium ne '') {
 +                push @criteria, $criterium;
 +                push @num_alt, --$num;
 +            }
 +        }
 +    }
 +    return ($fortran_filter) if !@criteria;
 +
 +    # Check that the criteria are properly formatted
 +    foreach my $criterium (@criteria) {
 +        # Naked parameter possibly preceded by '!' is ok
 +        next if $criterium =~ /^!?\w+$/;
 +
 +        my $op = (split / +/, $criterium)[1];
 +        if (!defined($op) or grep(/[+*?\\]/, $op)
 +            or !grep(/^$op$/, @allowed_operators) ) {
 +            print "Error in $filt:\ncriterium is badly formatted"
 +                . " or operator not supported:\n$criterium";
 +            exit 1;
 +        }
 +    }
 +    return ($fortran_filter, \@criteria, \@num_alt);
 } }
  
Line 240: Line 346:
         my @p = split /,/, $params;         my @p = split /,/, $params;
         foreach my $name (@p) {         foreach my $name (@p) {
 +            $name =~ s/^\s+//;
 +            $name =~ s/\s+$//;
             if ($name =~ /^!/) {             if ($name =~ /^!/) {
                 $name = substr $name, 1;                 $name = substr $name, 1;
Line 280: Line 388:
         $line =~ s/^\s+//;         $line =~ s/^\s+//;
         $line =~ s/\s+$//;         $line =~ s/\s+$//;
-         next if !$line || $line =~ /^#/;+        next if !$line || $line =~ /^#/;
         my ($param, $transform) = split /=/, $line, 2;         my ($param, $transform) = split /=/, $line, 2;
         die "Invalid transformation in $transform_file:\n$line\n"         die "Invalid transformation in $transform_file:\n$line\n"
Line 288: Line 396:
         $transform_of{$param} = $transform;         $transform_of{$param} = $transform;
     }     }
 +    close $TRANSFORM or die "Cannot close $transform_file: $!";
     return \%transform_of;     return \%transform_of;
 } }
  
-sub read_filter_file +# Return true (1) if observation is to be filtered, i.e. does not 
-    my $filter_file = shift; +# comply with at least one line in filter file, where each line is one 
-    return if !$filter_file;+# or more alternatives <param> or <param> <operator> <value> 
 +sub filter_obs 
 +    my $msg_ref = shift; 
 +    my $criteria_ref = shift; 
 +    my $num_alt_ref = shift; # gives the number of alternative 
 +                             # criteria still to be checked 
 +    return unless $criteria_ref;
  
-    my @allowed_operators = +    my @ascii_params qw(aircraft call_sign icao_id name obstime type wigosid);
-        ('=', +
-         '<', +
-         '<=', +
-         '>', +
-         '>=', +
-         '!=', +
-     )+
-    my @criteria;+
  
-    open my $FILTER, '<', $filter_file +    # Note that the loop counter $i might be changed in the loop 
-        or die "Cannot open $filter_file: $!"+    for (my $i=0; $i @{$criteria_ref}; $i++) { 
-    # Skip the criteria meant for Fortran parsing, i.e. proceed to +        my $num_alt = $num_alt_ref->[$i]
-    # first line following a blank line +        my $criterium = $criteria_ref->[$i]; 
-    while (<$FILTER>) { +        my ($f_param, $f_operator, $f_value) split +/, $criterium, 3;
-        last if $=/^\s*$/; +
-    }+
  
-    Read the filter criteria meant for Perl parsingskipping blank +        First check for !$parmeaning $par should not be in the observation 
-    # lines and comment lines +        if (substr($f_param,0,1) eq '!') { 
-    while (my $line = <$FILTER>) { +            $f_param substr($f_param,1)
-        $line =~ s/^\s+//+            # If parameter is present, criterium is not fullfilled 
-        $line =~ s/\s+$//; +            if (exists $msg_ref->{$f_param}) { 
-        next if !$line || $line =~ /^#/; +                next if $num_alt# More alternatives to check
-        push @criteria, $line; +
-    }+
  
-    Check that the criteria are properly formatted +                Criterium not fulfilled and no more alternatives to 
-    foreach my $criterium (@criteria) +                # check. This observation should be filtered away 
-        Naked parameter is ok +                return 1; 
-        next if $criterium =~ /^\w+$/;+            } else 
 +                Criterium fulfilled. No need to check alternative criteria 
 +                $+$num_alt if $num_alt; 
 +                next; 
 +            } 
 +        }
  
-        my $op = (split / +/$criterium)[1]; +        # If parameter not present, criterium is obviously not fullfilled 
-        if (!defined($op) or grep(/[+*?\\]/, $op) or !grep(/^$op$/, @allowed_operators) ) { +        if (not exists $msg_ref->{$f_param}) { 
-            print "Error in $filter_file, line $. is badly formatted" +            next if $num_alt
-                . " or operator not supported:\n$criterium"+            return 1;
-            exit 1;+
         }         }
-    } 
-    return \@criteria; 
-} 
  
-Return true (1if observation is to be filteredi.e. does not +        my $msg_value = $msg_ref->{$f_param}; 
-# comply with at least one of the <param> or <param> <operator> +        If a naked parameter criterium, we already know parameter is 
-<value> filter criteria in filter file +        # present (as found in previous check), so criterium is fulfilled 
-sub filter_obs { +        if (not defined $f_operator) { 
-    my $msg_ref shift+            if ($num_alt) { 
-    my $criteria_ref = shift;+                No need to check the alternative criteria 
 +                $i +$num_alt
 +            } 
 +            next; 
 +        }
  
-    my @ascii_params = qw(call_sign icao_id name obstime type); 
- 
-    foreach my $criterium (@$criteria_ref) { 
-        my ($f_param, $f_operator, $f_value) = split / +/, $criterium, 3; 
-        return 1 unless exists $msg_ref->{$f_param}; 
-        next if not defined $f_operator; # Naked parameter criterium. Parameter 
-                                         # present, so criterium fulfilled 
         chomp $f_value;         chomp $f_value;
-        if ($f_operator eq '=') {+        my $op; 
 +        if ($f_operator eq '<' 
 +            || $f_operator eq '<=' 
 +            || $f_operator eq '>' 
 +            || $f_operator eq '>=' 
 +            || $f_operator eq '=~' 
 +            || $f_operator eq '!~') { 
 +            $op = $f_operator; 
 +        } elsif ($f_operator eq '=' 
 +                 || $f_operator eq '!=') {
             if (grep {$_ eq $f_param} @ascii_params) {             if (grep {$_ eq $f_param} @ascii_params) {
-                $msg_ref->{$f_param} =~ s/\s*$//; +                $msg_value =~ s/\s*$//; 
-                return 1 unless $msg_ref->{$f_param} eq $f_value;+                $op = ($f_operator eq '=') ? 'eq' : 'ne';
             } else {             } else {
-                return 1 unless $msg_ref->{$f_param} == $f_value;+                $op = ($f_operator eq '=') ? '==' : '!=';
             }             }
-        } elsif ($f_operator eq '<'+        } else 
-            return 1 unless $msg_ref->{$f_param} < $f_value+            die "Internal error: unknown operator '$f_operator'"
-        } elsif ($f_operator eq '<=') { +        } 
-            return 1 unless $msg_ref->{$f_param} <= $f_value; + 
-        } elsif ($f_operator eq '>') { + # Some parameters might need special massaging 
-            return 1 unless $msg_ref->{$f_param} > $f_value; +        if ($f_operator !~ /~/) { 
-        } elsif ($f_operator eq '>=') { +            if ($f_param eq 'wmonr' || $f_param eq 'buoy_id') { 
-            return 1 unless $msg_ref->{$f_param} >= $f_value; +                # Make non octal by removing leading 0 
-        } elsif ($f_operator eq '!=') { +                $msg_value =~ s/^0+//; 
-            if (grep {$_ eq $f_param} @ascii_params) { +                $f_value =~ s/^0+// if $f_value != 0
-                $msg_ref->{$f_param} =~ s/\s*$//; +            } elsif ($f_param eq 'nationalnr') { 
-                return 1 unless $msg_ref->{$f_param} ne $f_value; +                # Convert to a pure numerical value (float). For 001101 
-            else { +                # State id only numbers between 100 and 699 are operational 
-                return 1 unless $msg_ref->{$f_param} != $f_value;+                $msg_value =~ s/_0*/./; 
 +                $f_value =~ s/_0*/./
 +            } elsif ($f_param eq 'obstime') { 
 +                # Convert to a pure numerical value (float) 
 +                $msg_value =~ s/[-:]//g; 
 +                $msg_value =~ s/ /./; 
 +                $f_value =~ s/[-:']//g
 +                $f_value =~ s/^ +//; 
 +                $f_value =~ s/ +$//; 
 +                $f_value =~ s/ /./; 
 +            } elsif ($f_param eq 'name') { 
 +                # Add or correct quoting to "" and ignore casing 
 +                $msg_value = '"' . lc $msg_value . '"'; 
 +                $f_value = lc $f_value; 
 +                if ($f_value =~ /^'.*'$/) { 
 +                    $f_value =~ s/^'//; 
 +                    $f_value =~ s/'$//
 +                
 +                if ($f_value !~ /^".*"$/{ 
 +                    $f_value '"'$f_value . '"'; 
 +                }
             }             }
         }         }
 +
 + my $condition = "$msg_value $op $f_value";
 + # Some values should be string values
 + if ($f_operator =~ /~/) {
 +     $condition = "q{$msg_value} $op $f_value";
 + } elsif (grep {$_ eq $f_param} @ascii_params) {
 +     $condition = "q{$msg_value} $op q{$f_value}";
 + }
 +
 +        # Finally, do the criterium check
 + if (eval $condition) {
 +     # No need to check the remaining alternative criteria
 +     $i += $num_alt if $num_alt;
 +     next;
 + } else {
 +     next if $num_alt;
 +     return 1;
 + }
     }     }
  
Line 381: Line 529:
 } }
  
 +# Return the type of station requested, and the station
 +# identifications. Leave some leeway for how to list wmonr and
 +# nationalnr (leading 0's might be omitted - added here)
 sub get_requested_stations { sub get_requested_stations {
     my $req_stations = shift;     my $req_stations = shift;
     return if !$req_stations;     return if !$req_stations;
-    die "Station list must start with 'wmonr=', 'call_sign=' or 'buoy_id='" +    die "Station list must start with 'wmonr=', 'nationalnr=',
-        unless $req_stations =~ /^wmonr=|call_sign=|buoy_id=/;+        . "'call_sign=''buoy_id=', 'aircraft=', 'icao_id' or 'wigosid='" 
 +        unless $req_stations 
 +        =~ /^(wmonr=|nationalnr=|call_sign=|buoy_id=|aircraft=|icao_id=|wigosid=)/;
  
     my ($id, $rest) = split /=/, $req_stations;     my ($id, $rest) = split /=/, $req_stations;
Line 396: Line 549:
             $station += 1000 if $station < 1000;             $station += 1000 if $station < 1000;
             $station = sprintf("%05d", $station);             $station = sprintf("%05d", $station);
 +            push @req_stn, $station;
 +        }
 +        return ($id, \@req_stn);
 +    } elsif ($id eq 'nationalnr') {
 +        my @req_stn;
 +        foreach my $station (@stations) {
 +            # Turn national station number into 10 digits
 +            my ($state_id, $national_id) = split /_/, $station;
 +            die "Uncorrected formatted station: '$station' in station list"
 +                if !defined $national_id or $national_id eq '';
 +            $station = $state_id . '_' . sprintf("%010d", $national_id);
             push @req_stn, $station;             push @req_stn, $station;
         }         }
Line 435: Line 599:
     }     }
  
-    my @ascii_params = qw(call_sign icao_id obstime name type);+    my @ascii_params = qw(aircraft call_sign icao_id obstime name type wigosid);
     my $lexical_sort = grep {$_ eq $sort_on} @ascii_params;     my $lexical_sort = grep {$_ eq $sort_on} @ascii_params;
  
Line 493: Line 657:
  
   bufrdump.pl <bufr file(s)>   bufrdump.pl <bufr file(s)>
-      [--filter <filter file>] +      [--filter <filter file | filter list>] 
-      [--param <parameter file | parameter list> [--csv]]+      [--param <parameter file | parameter list> [--csv [--delimiter <del>]]
       [--sort]       [--sort]
       [--sort_on <parameter>[-]]       [--sort_on <parameter>[-]]
       [--station <station list>]       [--station <station list>]
       [--transform <transformation file>]       [--transform <transformation file>]
-      [--lon1 x1] +      [--lon1 <x1>
-      [--lat1 y1] +      [--lat1 <y1>
-      [--lon2 x2] +      [--lon2 <x2>
-      [--lat2 x2]+      [--lat2 <y2>
 +      [--obstype <amdar|ocea|surface|sounding|sounding->]
       [--tablepath <path to BUFR tables>]       [--tablepath <path to BUFR tables>]
       [--help]       [--help]
Line 521: Line 686:
  
  
-  --filter <filter file> +  --filter <filter file | filter list
-                  Decode observations meeting criteria in <filter fileonly +                  Decode observations meeting criteria in filter file or 
-  --param <parameter file | parameter list> [--csv] +                  filter list only 
-                  Print parameters in parameter file or comma separated +  --param <parameter file | parameter list> [--csv [--delimiter <del>]
-                  list (e.g. wmonr,TA) only, in same order as they occur +                  Print parameters in parameter file or comma 
-                  there. If --csv, the parameters vill be printed using +                  separated list (e.g. wmonr,TA) only, in same order 
-                  the CSV (comma separated values) format+                  as they occur there. If using --csv possibly 
 +                  followed by --delimiter <del>, the parameters vill 
 +                  be printed using the CSV (comma-separated values) 
 +                  format, with the delimiter del (default is ';')
   --sort          Sort the decoded observations on station identification;   --sort          Sort the decoded observations on station identification;
-                  first stations with wmonr, then stations with call sign+                  first stations with wmonr, then stations with nationalnr
-                  then stations with buoy_id (others left out)+                  call_sign, buoy_id, aircraft, icao_id or wigosid 
 +                  (others left out)
   --sort_on <parameter>[-] Sort the decoded observations on increasing   --sort_on <parameter>[-] Sort the decoded observations on increasing
-                  values of <parameter>, or decreasing values if a '-'+                  values of parameter, or decreasing values if a '-'
                   follows the parameter name. E.g. --sort_on TA- will                   follows the parameter name. E.g. --sort_on TA- will
                   sort on decreasing temperatures. Observations not                   sort on decreasing temperatures. Observations not
Line 538: Line 707:
                   except when --sort_on is combined with --sort (in which                   except when --sort_on is combined with --sort (in which
                   case sorting is done firstly on station identification,                   case sorting is done firstly on station identification,
-                  secondly on <parameterwith missing values printed first)+                  secondly on parameter with missing values printed first)
   --station <station list>   --station <station list>
                   Print observations for stations in station list only,                   Print observations for stations in station list only,
Line 544: Line 713:
   --transform <transformation file>   --transform <transformation file>
                   Do the transformations of parameter values listed in                   Do the transformations of parameter values listed in
-                  <transformation file> +                  transformation file 
-  --lon1 x1       Decode observations with longitude >= x1 only +  --lon1 <x1>     Decode observations with longitude >= x1 only 
-  --lat1 y1       Decode observations with latitude >= y1 only +  --lat1 <y1>     Decode observations with latitude >= y1 only 
-  --lon2 x2       Decode observations with longitude <= x2 only +  --lon2 <x2>     Decode observations with longitude <= x2 only 
-  --lat2 y2       Decode observations with latitude <= y2 only+  --lat2 <y2>     Decode observations with latitude <= y2 only
                   x1,y1,x2,y2 should be decimal degrees                   x1,y1,x2,y2 should be decimal degrees
 +  --obstype <amdar|ocea|surface|sounding|sounding->]
 +                  Force observation type. If this option is not set,
 +                  will make an educated guess of observation type
 +                  based on metadata in section 1 of each BUFR message
   --tablepath <path to BUFR tables>   --tablepath <path to BUFR tables>
                   Set path to BUFR tables (overrides ENV{BUFR_TABLES})                   Set path to BUFR tables (overrides ENV{BUFR_TABLES})
Line 558: Line 731:
  
 To avoid having to use the C<--tablepath> option, you are adviced to To avoid having to use the C<--tablepath> option, you are adviced to
-set the invironment variable BUFR_TABLES to the directory where your+set the environment variable BUFR_TABLES to the directory where your
 BUFR tables are located (unless the default path provided by BUFR tables are located (unless the default path provided by
 bufrdump.pl works for you). bufrdump.pl works for you).
Line 567: Line 740:
 temperature to be printed for a BUFR SYNOP file, either supply temperature to be printed for a BUFR SYNOP file, either supply
  
-  wmonr,call_sign,TA+  wmonr,nationalnr,call_sign,TA
  
 as argument to --params, or supply a <parameter file> which should as argument to --params, or supply a <parameter file> which should
Line 573: Line 746:
  
   wmonr   wmonr
 +  nationalnr
   call_sign   call_sign
   TA   TA
Line 593: Line 767:
 listing may for example start like listing may for example start like
  
-  wmonr,call_sign,TA +  wmonr;nationalnr;call_sign;TA 
-  01001,,-1.5 +  01001;;;-1.5 
-  ,LF5U,9.0+  ;;LF5U;9.0 
 + 
 +You can choose another delimiter than semicolon by use of option 
 +--delimiter <del>, e.g. --csv --delimiter ','
  
 Using --filter will decode only those observations that meet at least Using --filter will decode only those observations that meet at least
Line 601: Line 778:
 in <filter file>, where the BUFR descriptor criteria should come first in <filter file>, where the BUFR descriptor criteria should come first
 in filter file followed by a blank line, then comes the parameter in filter file followed by a blank line, then comes the parameter
-criteria which should match <param> or <param> <operator> <value> +criteria which should match <param> or !<param> or <param> <operator> 
-where operator is one of =, !=, <, <=, > and >=. An example filter +<value> where operator is one of =, !=, =~, !~, <, <=, > and >=. What 
-file is+follows =~ and !~ should be a Perl match regular expression. The parameter 
 +criteria may be phrased as alternatives by separating them with '|' on 
 +a single line. An example filter file is
  
   D: 001001 I2.2   D: 001001 I2.2
Line 620: Line 799:
  
 which decodes all observations with block number 01, two other which decodes all observations with block number 01, two other
-specific wmo stations and one specific ship, being manned stations and +specific wmo stations and one specific ship, where stations should be 
-having cloud cover different from 8 (but NN must be part of the +manned and have cloud cover with a value different from 8and have 
-message) and temperature between 5 and 9.5 degrees Celsius and +temperature between 5 and 9.5 degrees Celsiusand contain 
-containing precipitation for last 24 hours. Comment lines starting +precipitation for last 24 hours. Comment lines starting with # will be 
-with # will be ignored.+ignored.
  
-Another example: the simple filter file (starting with a blank line!)+Another example: the filter file (starting with a blank line!)
  
  
-  wmonr+  call_sign =~ /^L[A-N]..$/ 
 +  obstime >= '2012-02-10 06:00:00' 
 +  HW | HWA | PW | PWA 
 +  FF > 10 | FG_010 > 10 
 + 
 +will print only those ship observations for which the 4 character 
 +call_sign starts with 2 letters in the interval LA-LN, and having 
 +obstime larger or equal to the datetime given, and containing wave 
 +data (specifically: height or period of waves, manually or 
 +automatically measured), and with wind or 10 minutes gust more than 10 
 +m/s. 
 + 
 +For convenience, when there are no BUFR descriptor criteria, you might 
 +provide the filter criteria on the command line. Example: 
 + 
 +--filter 'wmonr,TA > 0,RR_12 | RR_24, !FF'
  
-will print only those observations containing a wmonr (skipping +will decode only observations with wmonr, having positive temperature 
-ships).+and containing precipitation for 12 or 24 hours and not reporting 
 +wind. If (like for --paramthe filter list consists of one criterium 
 +only, a comma must be appended.
  
 To avoid the need of creating a filter file when observations for some To avoid the need of creating a filter file when observations for some
 few stations are requested, you can provide the stations in a comma few stations are requested, you can provide the stations in a comma
-separated list after option --station. Three examples:+separated list after option --station. Some examples:
  
   --station wmonr=01001,01152,01492   --station wmonr=01001,01152,01492
 +  --station nationalnr=614_0050410003,637_108
   --station call_sign=LF5U   --station call_sign=LF5U
   --station buoyid=64607,64609   --station buoyid=64607,64609
 +  --station aircraft=EU3421,JHCWUURA
 +  --station icao_id=ENGM,ENBO
 +  --station wigosid=0-376-0-511,0-20000-0-01492
  
 You cannot mix different kinds of stations this way (before '=' you You cannot mix different kinds of stations this way (before '=' you
-must choose either wmonr, call_sign or buoy_id). Note also that +must choose either wmonr, nationalnr, call_signbuoy_id, aircraft, 
-providing the stations in the BUFR descriptor part (first part) of the +icao_id or wigosid). Note also that providing the stations in the BUFR 
-filter file will speed up execution time considerably, compared to +descriptor part (first part) of the filter file will speed up 
-using option --station. It is possible to combine --filter with +execution time considerably, compared to using option --station. It is 
---station if done with some care, e.gspecifying WMO block 01 and the +possible to combine --filter with --station if done with some care, 
-required parameters in filter file, then the requested stations in +e.gspecifying WMO block 01 and the required parameters in filter 
-station list.+file, then the requested stations in station list.
  
 The --transform option is provided mainly to be able to use other The --transform option is provided mainly to be able to use other
Line 672: Line 872:
 is to be applied for sky not all covered by clouds, you should use NN is to be applied for sky not all covered by clouds, you should use NN
 != 100 instead of NN != 8 in filter file. != 100 instead of NN != 8 in filter file.
 +
 +The --obstype option might be handy in some special cases, like when
 +you are interested only in the surface part of oceanographic data
 +(then use '--obstype surface'), or when you want to see only levels
 +with vss>0 in high resolution radiosonde data (then use '--obstype
 +sounding-'), or when data category and/or data sub-category in the
 +BUFR messages have unusual values.
  
 =head1 AUTHOR =head1 AUTHOR
  • bufr.pm/bufrdump.pl_source.1301657283.txt.gz
  • Last modified: 2022-05-31 09:23:11
  • (external edit)