Хранилища Subversion geo-modmetar

# AW: create fork Geo::METAR as Geo::ModMETAR

#

# KH: fix the parser

# should be a finite state machine

# - metar has rules what comes after what. but codes can be missing.

# (measurement not done) or //// (measurement broken at the moment)

# so given a state counter, it can stay the same or go up one or more states,

# but it can never go down

#

# info on the last bit which is actually a forecast: (German)

# http://www.wetterklima.de/flug/metar/Metarvorhersage.htm

#

# more info here (dutch, and txt 707 is not standard metar)

# http://www.vwkweb.nl/index.html?http://www.vwkweb.nl/weerinfo/weerinfo_teletekst707.html

# and also (dutch)

# http://www.gids.nl/weather/eheh/metari.html

#

# 'METAR decoding in Europe'

# http://users.hol.gr/~chatos/VATSIM/TM/metar.html

#

# english explanation

# http://booty.org.uk/booty.weather/metinfo/codes/METAR_decode.htm

#

# canadian explanation

# http://meteocentre.com/doc/metar.html

#

# 'METAR decoding, TAF decoding'

# http://stoivane.kapsi.fi/metar/

#

# This module is used for decoding NWS METAR code.

# Example METARs

#

# Findlay, Ohio

# KFDY 251450Z 21012G21KT 8SM OVC065 04/M01 A3010 RMK SLP201 57014

#

# Toledo, Ohio

# KTOL 251451Z 23016G22KT 8SM CLR 04/00 A3006 RMK AO2 SLP185 T00440000 56016 

#

# Cleveland, Ohio

# KCLE 251554Z 20015KT 10SM FEW055 OVC070 03/M02 A3011 RMK AO2 SLP205 T00331017

#

# Houston, Texas

# KHST 251455Z 06017G22KT 7SM FEW040 BKN330 25/18 A3016 RMK SLP213 8/508

# 9/205 51007

#

# LA

#

# KLAX 251450Z 07004KT 7SM SCT100 BKN200 14/11 A3005 RMK AO2 SLP173

# T01390111 56005

#

# Soesterberg

#

# EHSB 181325Z 24009KT 8000 -RA BR FEW011 SCT022 OVC030 07/06 Q1011 WHT WHT TEMPO GRN

# For METAR info, please see

# http://tgsv5.nws.noaa.gov/oso/oso1/oso12/metar.htm

# moved

# http://metar.noaa.gov/

#

# in scary detail (metar coding)

#

# http://metar.noaa.gov/table_master.jsp?sub_menu=yes&show=fmh1ch12.htm&dir=./handbook/&title=title_handbook

#

# The METAR specification is dictated in the Federal Meteorological Handbook

# which is available on-line at:

# http://tgsv5.nws.noaa.gov/oso/oso1/oso12/fmh1.htm

# General Structure is:

# TYPE, SITE, DATE/TIME, WIND, VISIBILITY, CLOUDS, TEMPERATURE, PRESSURE, REMARKS

# Specifically:

# TYPE (optional)

# METAR or SPECI

# METAR: regular report

# SPECI: special report

# SITE (required, only once)

#

# 4-Char site identifier (KLAX for LA, KHST for Houston)

# DATE/TIME (required, only once)

#

# 6-digit time followed by "Z", indicating UTC

# REPORT MODIFIER (optional)

# AUTO or COR

# AUTO = Automatic report (no human intervention)

# COR = Corrected METAR or SPECI

# WIND (group)

#

# Wind direction (\d\d\d) and speed (\d?\d\d) and optionaling gusting

# information denoted by "G" and speed (\d?\d\d) followed by "KT", for knots.

#

# Wind direction MAY be "VRB" (variable) instead of a compass direction.

#

# Variable Wind Direction (Speeds greater than 6 knots).  Variable wind

# direction with wind speed greater than 6 knots shall be coded in the

# format, dndndnVdxdxdx

#

# Calm wind is recorded as 00000KT.

# VISIBILITY (group)

#

# Visibility (\d+) followed by "SM" for statute miles or no 'SM' for meters

# (european)

#

# May be 1/(\d)SM for a fraction.

#

# May be M1/\d)SM for less than a given fraction. (M="-")

#

# \d\d\d\d according to KNMI

# lowest horizontal visibility (looking around)

# round down

# 0000 - 0500m in steps of 0050m

# 0500 - 5000m in steps of 0100m

# 5000 - 9999m in steps of 1000m

# 10km or more is 9999

# RUNWAY Visual Range (Group)

#

# R(\d\d\d)(L|C|R)?/((M|P)?\d\d\d\d){1,2}FT

#

# Where:

#  $1 is the runway number.

#  $2 is the runway (Left/Center/Right) for parallel runways.

#  $3 is the reported visibility in feet.

#  $4 is the MAXIMUM reported visibility, making $3 the MINIMUM.

#

#  "M" beginning a value means less than the reportable value of \d\d\d\d.

#  "P" beginning a value means more than the reportable value of \d\d\d\d.

#

#  new

#

#  R(\d\d\d[LCR]?)/([MP]?\d\d\d\d)(V[MP]?\d\d\d\d)?FT

#

# $1 runway number + Left/Center/Right

# $2 visibility feet

# $3 Varying feet

# M = less than

# P = more than

# WEATHER (Present Weather Group)

#

# See table in Chapter 12 of FMH-1.

# CLOUDS (Sky Condition Group)

#

# A space-separated grouping of cloud conditions which will contain at least

# one cloud report. Examples: "CLR", "BKN330", "SCT100", "FEW055", "OVC070"

# The three-letter codes represent the condition (Clear, Broken, Scattered,

# Few, Overcast) and the numbers (\d\d\d) represent altitlude/100.

#

# The report may have a trailing CB (cumulonimbus) or TCU (towering

# cumulus) appended. ([A-Z]{2,3})?(\d\d\d)(CB|TCU)?

# Vertical visibility (VV)

#

# VV

# This group is reported when the sky is obscured. VV is the group indicator,

# and hshshs is the vertical visibility in units of 30 metres

# (hundreds of feet).

#  

#  hshshs - Examples of Encoding

#  HEIGHT               METAR CODE

#  100 ft       (30 metres)     001

#  450 ft       (135 metres)    004

#  2,700 ft     (810 metres)    027

#  12,600 ft    (3,780 metres)  1300

#

# source http://meteocentre.com/doc/metar.html

# 

# TEMPERATURE and DEW POINT

#

# (M?\d\d)/(M?\d\d) where $1 is the current temperature in degrees celcius,

# and $2 is the current dewpoint in degrees celcius.

#

# The "M" signifies a negative temperature, so converting the "M" to a

# "-" ought to suffice.

# PRESSURE

#

# The pressure, or altimeter setting, at the reporting site recorded in whole

# hectopascals (starts with a Q) or inches of mercury (Hg) minus the decimal

# point (starts with an A). It should always look like ([AQ]\d\d\d\d).

#

# KNMI: Q\d\d\d\d pressure in hPa calculated for sea level

# REMARKS

#

# Remarks contain additional information. They are optional but often

# informative of special conditions.

#

# Remarks begin with the "RMK" keyword and continue to the end of the line.

#

# trend group

#

# color codes BLU WHT GRN YLO AMB RED

# BLACK: vliegveld dicht

# future trend 

# NOSIG no significant change

# TEMPO temporary change

# WHT WHT TEMPO GRN = current white, prediction white temporary green

# NSW no significant weather

# AT at a given time

# PROB30 probability 30%

# BECMG becoming

# BECMG (weather) FM \d\d\d\d TL \d\d\d\d = from until utc times

# BECMG (weather) AT \d\d\d\d = at utc time

# BECMG (weather) TL \d\d\d\d = change until utc time

# BECMG 2000 visibility

# BECMG NSW weather type

# etc etc

# FCST CANCEL (2 tokens!) Forecast cancel: no further forecasts for a while

### Package Definition

package Geo::ModMETAR; # Package based on Debian Geo::METAR

## Required Modules

use 5.005;

use strict;

use vars qw($AUTOLOAD $VERSION);

use Carp 'cluck';

$VERSION = '1.0'; # Based on Debian Geo::METAR 1.15

##

## Lookup tables

##

my %_weather_types = (

    MI => 'shallow',

    PI => 'partial',

    BC => 'patches',

    DR => 'drizzle',

    BL => 'blowing',

    SH => 'shower(s)',

    TS => 'thunderstorm',

    FZ => 'freezing',

    DZ => 'drizzle',

    RA => 'rain',

    SN => 'snow',

    SG => 'snow grains',

    IC => 'ice crystals',

    PE => 'ice pellets',

    GR => 'hail',

    GS => 'small hail/snow pellets',

    UP => 'unknown precip',

    BR => 'mist',

    FG => 'fog',

    PRFG => 'fog banks',  # officially PR is a modifier of FG

    FU => 'smoke',

    VA => 'volcanic ash',

    DU => 'dust',

    SA => 'sand',

    HZ => 'haze',

    PY => 'spray',

    PO => 'dust/sand whirls',

    SQ => 'squalls',

    FC => 'funnel cloud(tornado/waterspout)',

    SS => 'sand storm',

    DS => 'dust storm',

);

my $_weather_types_pat = join("|", keys(%_weather_types));

my %_weather_types_ru = (

    MI => 'мелкий',

    PI => 'частный',

    BC => 'местами',

    DR => 'небольшой дождь',

    BL => 'порывистый ветер',

    SH => 'ливень',

    TS => 'гроза',

    FZ => 'заморозки',

    DZ => 'небольшой дождь',

    RA => 'дождь',

    SN => 'снег',

    SG => 'снег гранулами',

    IC => 'ледяные кристаллы',

    PE => 'ледяные шарики',

    GR => 'град',

    GS => 'небольшой град',

    UP => '',

    BR => 'легкий туман',

    FG => 'туман',

    PRFG => 'образование тумана',

    FU => 'дым',

    VA => 'вулканический пепел',

    DU => 'пыль',

    SA => 'песок',

    HZ => 'дымка',

    PY => 'водяная пыль',

    PO => 'песчаные или пылевые вихри',

    SQ => 'шквалистый ветер',

    FC => 'торнадо',

    SS => 'песчанная буря',

    DS => 'пылевая буря',

);

my $_weather_types_ru_pat = join("|", keys(%_weather_types_ru));

my %_sky_types = (

    SKC => "Sky Clear",

    CLR => "Sky Clear",

    SCT => "Scattered",

    BKN => "Broken",

    FEW => "Few",

    OVC => "Solid Overcast",

    NSC => "No significant clouds",

    NCD => "No cloud detected",

);

my %_sky_types_ru = (

    SKC => "ясно",

    CLR => "ясно",

    SCT => "переменная облачность",

    BKN => "переменная облачность",

    FEW => "слабая облачность",

    OVC => "сплошная облачность",

    NSC => "незначительная облачность",

    NCD => "безоблачно",

);

my %_trend_types = (

    BLU => "8 km view",

    WHT => "5 km view",

    GRN => "3.7 km view",

    YLO => "1.6 km view",

    AMB => "0.8 km view",

    RED => "< 0.8 km view",

    BLACK => "airport closed",

    NOSIG => "No significant change",

    TEMPO => "Temporary change",

    NSW => "No significant weather",

    PROB => "Probability",

    BECMG => "Becoming",

    LAST => "Last",

);

my $_trend_types_pat = join("|", keys(%_trend_types));

##

## Constructor.

##

sub new

{

    my $this = shift;

    my $class = ref($this) || $this;

    my $self = {};

    ##

    ## UPPERCASE items have documented accssor functions (methods) or

    ## use AUTOLOAD, while lowercase items are reserved for internal

    ## use.

    ##

    $self->{VERSION}       = $VERSION;          # version number

    $self->{METAR}         = undef;             # the actual, raw METAR

    $self->{TYPE}          = undef;             # the type of report

    $self->{SITE}          = undef;             # site code

    $self->{DATE}          = undef;             # when it was issued

    $self->{TIME}          = undef;             # time it was issued

    $self->{MOD}           = undef;             # modifier (AUTO/COR)

    $self->{WIND_DIR_DEG}  = undef;             # wind dir in degrees

    $self->{WIND_DIR_ENG}  = undef;             # wind dir in english (Northwest/Southeast)

    $self->{WIND_DIR_RUS}  = undef;             # wind dir in russian (Северо-западный/Юго-восточный)

    $self->{WIND_DIR_ABB}  = undef;             # wind dir in abbreviated english (NW/SE)

    $self->{WIND_KTS}      = undef;             # wind speed (knots)

    $self->{WIND_GUST_KTS} = undef;             # wind gusts (knots)

    $self->{WIND_MPH}      = undef;             # wind speed (MPH)

    $self->{WIND_GUST_MPH} = undef;             # wind gusts (MPH)

    $self->{WIND_MS}       = undef;             # wind speed (m/s)

    $self->{WIND_GUST_MS}  = undef;             # wind gusts (m/s)

    $self->{WIND_VAR}      = undef;             # wind variation (text)

    $self->{WIND_VAR_1}    = undef;             # wind variation (direction 1)

    $self->{WIND_VAR_2}    = undef;             # wind variation (direction 2)

    $self->{WIND_VAR_ENG_1}= undef;             # wind variation (text, direction 1)

    $self->{WIND_VAR_ENG_2}= undef;             # wind variation (text, direction 2)

    $self->{VISIBILITY}    = undef;             # visibility info

    $self->{RUNWAY}        = [ ];               # runway vis.

    $self->{RH}            = undef;             # relative humidity

    $self->{WEATHER}       = [ ];               # current weather

    $self->{WEATHER_LOG}   = [ ];               # weather log

    $self->{SKY}           = [ ];               # current sky (cloudcover)

    $self->{TEMP_F}        = undef;             # current temp, celsius

    $self->{TEMP_C}        = undef;             # converted to fahrenheit

    $self->{TEMP_WC}       = undef;             # current windchill temp, celsius

    $self->{DEW_F}         = undef;             # dew point, celcius

    $self->{DEW_C}         = undef;             # dew point, fahrenheit

    $self->{HOURLY_TEMP_F} = undef;             # hourly current temp, celcius

    $self->{HOURLY_TEMP_C} = undef;             # hourly converted to fahrenheit

    $self->{HOURLY_DEW_F}  = undef;             # hourly dew point, celcius

    $self->{HOURLY_DEW_C}  = undef;             # hourly dew point, fahrenheit

    $self->{HOURLY_PRECIP} = undef;             # hourly precipitation

    $self->{ALT}           = undef;             # altimeter setting (Hg)

    $self->{ALT_HP}        = undef;             # altimeter setting (hPa)

    $self->{SLP}           = undef;             # sea level pressure

    $self->{REMARKS}       = undef;             # remarks

    $self->{WEATHER_RAW}   = [ ];               # RAW data for weather

    $self->{WEATHER_RUS}   = [ ];               # current weather in Russian

    $self->{SKY_RAW}       = [ ];               # RAW data for sky

    $self->{SKY_RUS}       = [ ];               # current sky in Russian

    $self->{VISIBILITY_RUS}= undef;             # visibility info    

    $self->{tokens}        = [ ];               # the "token" list

    $self->{type}          = "METAR";           # the report type (METAR/SPECI)

                                                # default=METAR

    $self->{site}          = undef;             # the site code (4 chars)

    $self->{date_time}     = undef;             # date/time

    $self->{modifier}      = undef;             # the AUTO/COR modifier

    $self->{wind}          = undef;             # the wind information

    $self->{windtype}      = undef;             # the wind speed type (knots/meterpersecond/kilometersperhour)

    $self->{windvar}       = undef;             # the wind variation

    $self->{visibility}    = undef;             # visibility information

    $self->{runway}        = undef;             # runway visibility

    $self->{weather}       = [ ];               # current weather conditions

    $self->{sky}           = [ ];               # sky conditions (cloud cover)

    $self->{temp_dew}      = undef;             # temp and dew pt.

    $self->{alt}           = undef;             # altimeter setting

    $self->{pressure}      = undef;             # pressure (HPa)

    $self->{slp}           = undef;             # sea level pressure

    $self->{remarks}       = [ ];               # remarks

    $self->{debug}         = undef;             # enable debug trace

    bless $self, $class;

    return $self;

}

##

## Autoload for access methods to stuff in %fields hash. We should

## probably disallow access to the lower-case items as stated above,

## but I don't feel like being a Nazi about it. Besides, I haven't

## checked to see what that might break.

##

sub AUTOLOAD

{

    my $self = shift;

    if (not ref $self)

    {

        cluck "bad AUTOLOAD for obj [$self]";

    }

    if ($AUTOLOAD =~ /.*::(.*)/)

    {

        my $key = $1;

        ## Backward compatible temps...

        my %compat = (

                      F_TEMP    =>  'TEMP_F',

                      C_TEMP    =>  'TEMP_C',

                      F_DEW     =>  'DEW_F',

                      C_DEW     =>  'DEW_C',

                     );

        if ($compat{$key})

        {

            $key = $compat{$key};

        }

        ## Check for the items...

        if (exists $self->{$key})

        {

            return $self->{$key};

        }

        else

        {

            return undef;

        }

    }

    else

    {

        warn "strange AUTOLOAD problem!";

        return undef;

    }

}

##

## Get current version number.

##

sub version

{

    my $self = shift;

    print "version() called.\n" if $self->{debug};

    return $self->{VERSION};

}

##

## Take a METAR, tokenize, and process it.

##

sub metar

{

    my $self = shift;

    if (@_)

    {

        $self->{METAR} = shift;

        $self->{METAR} =~ s/\n//g;    ## nuke any newlines

        _tokenize($self);

        _process($self);

    }

    return $self->{METAR};

}

##

## Break {METAR} into parts. Stuff into @tokens.

##

sub _tokenize

{

    my $self = shift;

    my $tok;

    my @toks;

    # Split tokens on whitespace.

    @toks = split(/\s+/, $self->{METAR});

    $self->{tokens} = \@toks;

}

## Process @tokens to populate METAR values.

##

## This is a long and involved subroutine. It basically copies the

## @tokens array and treats it as a stack, popping off items,

## examining them, and see what they look like.  Based on their

## "apppearance" it takes care populating the proper fields

## internally.

sub _process

{

    my $self = shift;

    my @toks = @{$self->{tokens}};      # copy tokens array...

    my $tok;

    ## This is a semi-brute-force way of doing things, but the amount

    ## of data is relatively small, so it shouldn't be a big deal.

    ##

    ## Ideally, I'd have it skip checks for items which have been

    ## found, but that would make this more "linear" and I'd remove

    ## the pretty while loop.

        #

        # KH: modified to maintain state to not get lost in remarks and stuff

        # and be a lot better at parsing

        # states

        my $expect_type = 0;

        my $expect_site = 1;

        my $expect_datetime = 2;

        my $expect_modifier = 3;

        my $expect_wind = 4;

        my $expect_visibility = 5;

        my $expect_runwayvisual = 6;

        my $expect_presentweather = 7;

        my $expect_clouds = 8;

        my $expect_temperature = 9;

        my $expect_pressure = 10;

        my $expect_recentweather = 11;

        my $expect_remarks = 12;

        my $expect_usremarks = 13;

        my $parsestate = $expect_type;

        # windtypes

        my $wt_knots = 1;

        my $wt_mps = 2;

        my $wt_kph = 3;

    ## Assume standard report by default

    $self->{type} = "METAR";

    $self->{TYPE} = "Routine Weather Report";

    while (defined($tok = shift(@toks))) ## as long as there are tokens

    {

        print "trying to match [$tok] state is $parsestate\n" if $self->{debug};

        ##

        ## is it a report type?

        ##

        if (($parsestate == $expect_type) and ($tok =~ /(METAR|SPECI)/i))

        {

            $self->{type} = $tok;

            if ($self->{type} eq "METAR")

            {

                $self->{TYPE} = "Routine Weather Report";

            }

            elsif ($self->{type} eq "SPECI")

            {

                $self->{TYPE} = "Special Weather Report";

            }

            print "[$tok] is a report type.\n" if $self->{debug};

                        $parsestate = $expect_site;

            next;

        }

        ##

        ## is is a site ID?

        ##

        elsif (($parsestate <= $expect_site) and ($tok =~ /([A-Z]{4}|K[A-Z0-9]{3})/))

        {

            $self->{site} = $tok;

            print "[$tok] is a site ID.\n" if $self->{debug};

                        $parsestate = $expect_datetime;

            next;

        }

        ##

        ## is it a date/time?

        ##

        elsif (($parsestate == $expect_datetime) and ($tok =~ /\d{6,6}Z/i))

        {

            $self->{date_time} = $tok;

            print "[$tok] is a date/time.\n" if $self->{debug};

                        $parsestate = $expect_modifier;

            next;

        }

        ##

        ## is it a report modifier?

        ##

        elsif (($parsestate == $expect_modifier) and ($tok =~ /AUTO|COR|CC[A-Z]/i))

        {

            $self->{modifier} = $tok;

            print "[$tok] is a report modifier.\n" if $self->{debug};

                        $parsestate = $expect_wind;

            next;

        }

        ##

        ## is it wind information in knots?

        #

                # eew: KT seems to be optional

                # but making it optional fails on other stuff

                # sortafix: wind needs to be \d\d\d\d\d or VRB\d\d

                #      optional \d\d\d\d\dG\d\d\d (gust direction)

        elsif (($parsestate >= $expect_modifier) and ($parsestate < $expect_visibility) and ($tok =~ /(\d{3}|VRB)\d{2}(G\d{1,3})?(KT)?$/i))

        {

            $self->{wind} = $tok;

                        $self->{windtype} = $wt_knots;

            print "[$tok] is wind information in knots.\n" if $self->{debug};

                        $parsestate = $expect_wind; # stay in wind, it can have variation

            next;

        }

                ##

                ## is it wind information in meters per second?

                ##

                ## can be variable too

        elsif (($parsestate >= $expect_modifier) and ($parsestate < $expect_visibility) and ($tok =~ /^(\d{3}|VRB)\d{2}(G\d{2,3})?MPS$/))

        {

            $self->{wind} = $tok;

            print "[$tok] is wind information.\n" if $self->{debug};

                        $self->{windtype} = $wt_mps;

                        $parsestate = $expect_wind; # stay in wind, it can have variation

            next;

        }

                ##

                ## is it wind variation information?

                ##

                elsif (($parsestate >= $expect_wind) and ($parsestate < $expect_visibility) and ($tok =~ /^\d{3}V\d{3}$/))

                {

                        $self->{windvar} = $tok;

                        print "[$tok] is wind variation information.\n" if $self->{debug};

                        $parsestate = $expect_visibility;

                        next;

                }

                ##

                ## wind information missing at the moment?

                ##

                elsif (($parsestate >= $expect_wind) and ($parsestate < $expect_visibility) and ($tok =~ /^\/\/\/\/\/(KT|MPS)$/)){

                        print "[$tok] is missing wind information.\n" if $self->{debug};

                        $parsestate = $expect_visibility;

                        next;

                }

                ##

                ## is it visibility information in meters?

                ##

                elsif (($parsestate >= $expect_modifier) and ($parsestate < $expect_runwayvisual) and ($tok =~ /^\d{4}$/))

                {

                        $self->{visibility} = $tok;

            print "[$tok] is numerical visibility information.\n" if $self->{debug};

                        $parsestate = $expect_visibility;

            next;

        }

                ## auto visibility information in meters?

                elsif (($parsestate >= $expect_modifier) and ($parsestate < $expect_runwayvisual) and ($tok =~ /^\d{4}NDV$/))

                {

                        $self->{visibility} = $tok;

            print "[$tok] is automatic numerical visibility information.\n" if $self->{debug};

                        $parsestate = $expect_visibility;

            next;

        }

        ##

        ## is it visibility information in statute miles?

        ##

        elsif (($parsestate >= $expect_modifier) and ($parsestate < $expect_runwayvisual) and ($tok =~ /.*?SM$/i))

        {

            $self->{visibility} = $tok;

            print "[$tok] is statute miles visibility information.\n" if $self->{debug};

                        $parsestate = $expect_visibility;

            next;

        }

        ##

        ## is it visibility information with a leading digit?

                ##

                ## sample:

                ## KERV 132345Z AUTO 07008KT 1 1/4SM HZ 34/11 A3000 RMK AO2

                ##                           ^^^^^^^

        ##

        elsif (($parsestate >= $expect_modifier) and ($parsestate < $expect_runwayvisual) and  ($tok =~ /^\d$/))

        {

            $tok .= " " . shift(@toks);

            $self->{visibility} = $tok;

            print "[$tok] is multi-part visibility information.\n" if $self->{debug};

                        $parsestate = $expect_visibility;

            next;

        }

                ## visibility modifier

                elsif (($parsestate == $expect_visibility) and ($tok =~ /^\d{4}(N|S|E|W|NE|NW|SE|SW)$/))

                {

            print "[$tok] is a visibility modifier.\n" if $self->{debug};

            next;

                }

        ##

        ## is it runway visibility info?

        ##

                # KH: I've seen runway visibility with 'U' units

                # EHSB 121425Z 22010KT 1200 R27/1600U -DZ BKN003 OVC007 07/07 Q1016 AMB FCST CANCEL

                # U= going up, D= going down, N= no change

                # tendency of visual range, http://stoivane.kapsi.fi/metar/

        elsif (($parsestate >= $expect_modifier) and ($parsestate < $expect_presentweather) and ($tok =~ /R\d+(L|R|C)?\/P?\d+(VP?\d+)?(FT|D|U|N|\/)?$/i))

        {

            push (@{$self->{RUNWAY}},$tok);

            print "[$tok] is runway visual information.\n" if $self->{debug};

                        $parsestate = $expect_runwayvisual;

                        # there can be multiple runways, so stay at this state

            next;

        }

        ##

        ## is it current weather info?

        ##

        elsif (($parsestate >= $expect_modifier) and ($parsestate < $expect_clouds) and ($tok =~ /^(-|\+)?(VC)?($_weather_types_pat)+/i))