X-Git-Url: https://git.mxchange.org/?a=blobdiff_plain;f=simgear%2Fenvironment%2Fmetar.cxx;h=75109b75651d03eb413a08b9975696ec02315384;hb=f61fe80f6096bff9f7061849f713fdf9204aa945;hp=124fa3f2701a768753ca929a397290a8eed6fc10;hpb=3b3baa996b2485855534f08b68c52f3a40149598;p=simgear.git diff --git a/simgear/environment/metar.cxx b/simgear/environment/metar.cxx index 124fa3f2..75109b75 100644 --- a/simgear/environment/metar.cxx +++ b/simgear/environment/metar.cxx @@ -16,7 +16,7 @@ // // You should have received a copy of the GNU General Public License // along with this program; if not, write to the Free Software -// Foundation, 59 Temple Place, Suite 330, Boston, MA 02111-1307, USA +// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. // // $Id$ @@ -24,10 +24,14 @@ * @file metar.cxx * Interface for encoded Meteorological Aerodrome Reports (METAR). */ +#ifdef HAVE_CONFIG_H +# include +#endif #include +#include +#include -#include #include #include @@ -35,32 +39,30 @@ #define NaN SGMetarNaN +using std::string; +using std::map; +using std::vector; + /** - * The constructor takes a Metar string, or a four-letter ICAO code. In the - * latter case the metar string is downloaded from - * http://weather.noaa.gov/pub/data/observations/metar/stations/. + * The constructor takes a Metar string * The constructor throws sg_io_exceptions on failure. The "METAR" * keyword has no effect (apart from incrementing the group counter * @a grpcount) and can be left away. A keyword "SPECI" is * likewise accepted. * * @param m ICAO station id or metar string - * @param proxy proxy host (optional; default: "") - * @param port proxy port (optional; default: "80") - * @param auth proxy authorization information (optional; default: "") * * @par Examples: * @code * SGMetar *m = new SGMetar("METAR KSFO 061656Z 19004KT 9SM SCT100 OVC200 08/03 A3013"); * double t = m->getTemperature_F(); * delete m; - * - * SGMetar n("KSFO", "proxy.provider.foo", "3128", "proxy-password"); - * double d = n.getDewpoint_C(); + * @endcode */ -SGMetar::SGMetar(const string& m, const string& proxy, const string& port, const string& auth) : +SGMetar::SGMetar(const string& m) : _grpcount(0), + _x_proxy(false), _year(-1), _month(-1), _day(-1), @@ -74,32 +76,31 @@ SGMetar::SGMetar(const string& m, const string& proxy, const string& port, const _wind_range_to(-1), _temp(NaN), _dewp(NaN), - _pressure(NaN) + _pressure(NaN), + _rain(false), + _hail(false), + _snow(false), + _cavok(false) { - if (m.length() == 4 && isalnum(m[0]) && isalnum(m[1]) && isalnum(m[2]) && isalnum(m[3])) { - for (int i = 0; i < 4; i++) - _icao[i] = toupper(m[i]); - _icao[4] = '\0'; - _data = loadData(_icao, proxy, port, auth); - } else { - _data = new char[m.length() + 2]; // make room for " \0" - strcpy(_data, m.c_str()); - _url = _data; - } + _data = new char[m.length() + 2]; // make room for " \0" + strcpy(_data, m.c_str()); + _url = _data; + normalizeData(); _m = _data; _icao[0] = '\0'; // NOAA preample - scanPreambleDate(); + if (!scanPreambleDate()) + useCurrentDate(); scanPreambleTime(); // METAR header scanType(); if (!scanId() || !scanDate()) { delete[] _data; - throw sg_io_exception("metar data bogus (" + _url + ')'); + throw sg_io_exception("metar data bogus ", sg_location(_url)); } scanModifier(); @@ -125,8 +126,9 @@ SGMetar::SGMetar(const string& m, const string& proxy, const string& port, const if (_grpcount < 4) { delete[] _data; - throw sg_io_exception("metar data incomplete (" + _url + ')'); + throw sg_io_exception("metar data incomplete ", sg_location(_url)); } + _url = ""; } @@ -143,75 +145,19 @@ SGMetar::~SGMetar() } -/** - * If called with "KSFO" loads data from - * @code - * http://weather.noaa.gov/pub/data/observations/metar/stations/KSFO.TXT. - * @endcode - * Throws sg_io_exception on failure. Gives up after waiting longer than 10 seconds. - * - * @param id four-letter ICAO Metar station code, e.g. "KSFO". - * @param proxy proxy host (optional; default: "") - * @param port proxy port (optional; default: "80") - * @param auth proxy authorization information (optional; default: "") - * @return pointer to Metar data string, allocated by new char[]. - * @see rfc2068.txt for proxy spec ("Proxy-Authorization") - */ -char *SGMetar::loadData(const char *id, const string& proxy, const string& port, const string& auth) +void SGMetar::useCurrentDate() { - string host = proxy.empty() ? "weather.noaa.gov" : proxy; - string path = "/pub/data/observations/metar/stations/"; - path += string(id) + ".TXT"; - _url = "http://weather.noaa.gov" + path; - - SGSocket *sock = new SGSocket(host, port.empty() ? "80" : port, "tcp"); - sock->set_timeout(10000); - if (!sock->open(SG_IO_OUT)) { - delete sock; - throw sg_io_exception("cannot connect to " + host); - } - - string get = "GET "; - if (!proxy.empty()) - get += "http://weather.noaa.gov"; - get += path + " HTTP/1.0\r\n"; - sock->writestring(get.c_str()); - - if (!auth.empty()) { - get = "Proxy-Authorization: " + auth + "\r\n"; - sock->writestring(get.c_str()); - } - - sock->writestring("\r\n"); - - int i; - const int buflen = 512; - char buf[2 * buflen]; - - // skip HTTP header - while ((i = sock->readline(buf, buflen))) - if (i <= 2 && isspace(buf[0]) && (!buf[1] || isspace(buf[1]))) - break; - if (i) { - i = sock->readline(buf, buflen); - if (i) - sock->readline(&buf[i], buflen); - } - - sock->close(); - delete sock; - - char *b = buf; - scanBoundary(&b); - if (*b == '<') - throw sg_io_exception("no metar data available from " + _url); - - char *metar = new char[strlen(b) + 2]; // make room for " \0" - strcpy(metar, b); - return metar; + struct tm now; + time_t now_sec = time(0); +#ifdef _WIN32 + now = *gmtime(&now_sec); +#else + gmtime_r(&now_sec, &now); +#endif + _year = now.tm_year + 1900; + _month = now.tm_mon + 1; } - /** * Replace any number of subsequent spaces by just one space, and add * a trailing space. This makes scanning for things like "ALL RWY" easier. @@ -392,7 +338,7 @@ bool SGMetar::scanWind() if (gust != NaN) _gust_speed = gust * factor; _grpcount++; - return false; + return true; } @@ -432,11 +378,13 @@ bool SGMetar::scanVisibility() int modifier = SGMetarVisibility::EQUALS; // \d{4}(N|NE|E|SE|S|SW|W|NW)? if (scanNumber(&m, &i, 4)) { - if (*m == 'E') + if( strncmp( m, "NDV",3 ) == 0 ) { + m+=3; // tolerate NDV (no directional validation) + } else if (*m == 'E') { m++, dir = 90; - else if (*m == 'W') + } else if (*m == 'W') { m++, dir = 270; - else if (*m == 'N') { + } else if (*m == 'N') { m++; if (*m == 'E') m++, dir = 45; @@ -452,7 +400,7 @@ bool SGMetar::scanVisibility() m++, dir = 225; else dir = 180; - } + } if (i == 0) i = 50, modifier = SGMetarVisibility::LESS_THAN; else if (i == 9999) @@ -460,7 +408,6 @@ bool SGMetar::scanVisibility() distance = i; } else { // M?(\d{1,2}|\d{1,2}/\d{1,2}|\d{1,2} \d{1,2}/\d{1,2})(SM|KM) - modifier = 0; if (*m == 'M') m++, modifier = SGMetarVisibility::LESS_THAN; @@ -559,7 +506,7 @@ bool SGMetar::scanRwyVisRange() r._max_visibility._distance = to; if (*m == '/') // this is not in the spec! - *m++; + m++; if (*m == 'D') m++, r._min_visibility._tendency = SGMetarVisibility::DECREASING; else if (*m == 'N') @@ -579,61 +526,62 @@ bool SGMetar::scanRwyVisRange() static const struct Token special[] = { - "NSW", "no significant weather", - "VCSH", "showers in the vicinity", - "VCTS", "thunderstorm in the vicinity", - 0, 0 + { "NSW", "no significant weather" }, +/* { "VCSH", "showers in the vicinity" }, + { "VCTS", "thunderstorm in the vicinity" }, */ + { 0, 0 } }; static const struct Token description[] = { - "SH", "showers of", - "TS", "thunderstorm with", - "BC", "patches of", - "BL", "blowing", - "DR", "low drifting", - "FZ", "freezing", - "MI", "shallow", - "PR", "partial", - 0, 0 + { "SH", "showers of" }, + { "TS", "thunderstorm with" }, + { "BC", "patches of" }, + { "BL", "blowing" }, + { "DR", "low drifting" }, + { "FZ", "freezing" }, + { "MI", "shallow" }, + { "PR", "partial" }, + { 0, 0 } }; static const struct Token phenomenon[] = { - "DZ", "drizzle", - "GR", "hail", - "GS", "small hail and/or snow pellets", - "IC", "ice crystals", - "PE", "ice pellets", - "RA", "rain", - "SG", "snow grains", - "SN", "snow", - "UP", "unknown precipitation", - "BR", "mist", - "DU", "widespread dust", - "SG", "fog", - "SGBR", "fog bank", - "FU", "smoke", - "HZ", "haze", - "PY", "spray", - "SA", "sand", - "VA", "volcanic ash", - "DS", "duststorm", - "FC", "funnel cloud/tornado waterspout", - "PO", "well-developed dust/sand whirls", - "SQ", "squalls", - "SS", "sandstorm", - "UP", "unknown", // ... due to failed automatic acquisition - 0, 0 + { "DZ", "drizzle" }, + { "GR", "hail" }, + { "GS", "small hail and/or snow pellets" }, + { "IC", "ice crystals" }, + { "PE", "ice pellets" }, + { "RA", "rain" }, + { "SG", "snow grains" }, + { "SN", "snow" }, + { "UP", "unknown precipitation" }, + { "BR", "mist" }, + { "DU", "widespread dust" }, + { "FG", "fog" }, + { "FGBR", "fog bank" }, + { "FU", "smoke" }, + { "HZ", "haze" }, + { "PY", "spray" }, + { "SA", "sand" }, + { "VA", "volcanic ash" }, + { "DS", "duststorm" }, + { "FC", "funnel cloud/tornado waterspout" }, + { "PO", "well-developed dust/sand whirls" }, + { "SQ", "squalls" }, + { "SS", "sandstorm" }, + { "UP", "unknown" }, // ... due to failed automatic acquisition + { 0, 0 } }; -// (+|-|VC)?(NSW|MI|PR|BC|DR|BL|SH|TS|FZ)?((DZ|RA|SN|SG|IC|PE|GR|GS|UP){0,3})(BR|SG|FU|VA|DU|SA|HZ|PY|PO|SQ|FC|SS|DS){0,3} +// (+|-|VC)?(NSW|MI|PR|BC|DR|BL|SH|TS|FZ)?((DZ|RA|SN|SG|IC|PE|GR|GS|UP){0,3})(BR|FG|FU|VA|DU|SA|HZ|PY|PO|SQ|FC|SS|DS){0,3} bool SGMetar::scanWeather() { char *m = _m; string weather; const struct Token *a; + if ((a = scanToken(&m, special))) { if (!scanBoundary(&m)) return false; @@ -643,25 +591,35 @@ bool SGMetar::scanWeather() } string pre, post; + struct Weather w; if (*m == '-') - m++, pre = "light "; + m++, pre = "light ", w.intensity = LIGHT; else if (*m == '+') - m++, pre = "heavy "; + m++, pre = "heavy ", w.intensity = HEAVY; else if (!strncmp(m, "VC", 2)) - m += 2, post = "in the vicinity "; + m += 2, post = "in the vicinity ", w.vincinity=true; else - pre = "moderate "; + pre = "moderate ", w.intensity = MODERATE; int i; for (i = 0; i < 3; i++) { if (!(a = scanToken(&m, description))) break; + w.descriptions.push_back(a->id); weather += string(a->text) + " "; } + for (i = 0; i < 3; i++) { if (!(a = scanToken(&m, phenomenon))) break; + w.phenomena.push_back(a->id); weather += string(a->text) + " "; + if (!strcmp(a->id, "RA")) + _rain = w.intensity; + else if (!strcmp(a->id, "HA")) + _hail = w.intensity; + else if (!strcmp(a->id, "SN")) + _snow = w.intensity; } if (!weather.length()) return false; @@ -671,35 +629,37 @@ bool SGMetar::scanWeather() weather = pre + weather + post; weather.erase(weather.length() - 1); _weather.push_back(weather); + if( w.phenomena.size() > 0 ) + _weather2.push_back( w ); _grpcount++; return true; } static const struct Token cloud_types[] = { - "AC", "altocumulus", - "ACC", "altocumulus castellanus", - "ACSL", "altocumulus standing lenticular", - "AS", "altostratus", - "CB", "cumulonimbus", - "CBMAM", "cumulonimbus mammatus", - "CC", "cirrocumulus", - "CCSL", "cirrocumulus standing lenticular", - "CI", "cirrus", - "CS", "cirrostratus", - "CU", "cumulus", - "CUFRA", "cumulus fractus", - "NS", "nimbostratus", - "SAC", "stratoaltocumulus", // guessed - "SC", "stratocumulus", - "SCSL", "stratocumulus standing lenticular", - "ST", "stratus", - "STFRA", "stratus fractus", - "TCU", "towering cumulus", - 0, 0 + { "AC", "altocumulus" }, + { "ACC", "altocumulus castellanus" }, + { "ACSL", "altocumulus standing lenticular" }, + { "AS", "altostratus" }, + { "CB", "cumulonimbus" }, + { "CBMAM", "cumulonimbus mammatus" }, + { "CC", "cirrocumulus" }, + { "CCSL", "cirrocumulus standing lenticular" }, + { "CI", "cirrus" }, + { "CS", "cirrostratus" }, + { "CU", "cumulus" }, + { "CUFRA", "cumulus fractus" }, + { "NS", "nimbostratus" }, + { "SAC", "stratoaltocumulus" }, // guessed + { "SC", "stratocumulus" }, + { "SCSL", "stratocumulus standing lenticular" }, + { "ST", "stratus" }, + { "STFRA", "stratus fractus" }, + { "TCU", "towering cumulus" }, + { 0, 0 } }; - +#include // (FEW|SCT|BKN|OVC|SKC|CLR|CAVOK|VV)([0-9]{3}|///)?[:cloud_type:]? bool SGMetar::scanSkyCondition() { @@ -707,15 +667,29 @@ bool SGMetar::scanSkyCondition() int i; SGMetarCloud cl; + if (!strncmp(m, "//////", 6)) { + m += 6; + if (!scanBoundary(&m)) + return false; + _m = m; + return true; + } + if (!strncmp(m, "CLR", i = 3) // clear || !strncmp(m, "SKC", i = 3) // sky clear + || !strncmp(m, "NCD", i = 3) // nil cloud detected || !strncmp(m, "NSC", i = 3) // no significant clouds || !strncmp(m, "CAVOK", i = 5)) { // ceiling and visibility OK (implies 9999) m += i; if (!scanBoundary(&m)) return false; - cl._coverage = 0; - _clouds.push_back(cl); + + if (i == 3) { + cl._coverage = SGMetarCloud::COVERAGE_CLEAR; + _clouds.push_back(cl); + } else { + _cavok = true; + } _m = m; return true; } @@ -723,13 +697,13 @@ bool SGMetar::scanSkyCondition() if (!strncmp(m, "VV", i = 2)) // vertical visibility ; else if (!strncmp(m, "FEW", i = 3)) - cl._coverage = 1; + cl._coverage = SGMetarCloud::COVERAGE_FEW; else if (!strncmp(m, "SCT", i = 3)) - cl._coverage = 2; + cl._coverage = SGMetarCloud::COVERAGE_SCATTERED; else if (!strncmp(m, "BKN", i = 3)) - cl._coverage = 3; + cl._coverage = SGMetarCloud::COVERAGE_BROKEN; else if (!strncmp(m, "OVC", i = 3)) - cl._coverage = 4; + cl._coverage = SGMetarCloud::COVERAGE_OVERCAST; else return false; m += i; @@ -742,7 +716,7 @@ bool SGMetar::scanSkyCondition() } else if (!scanNumber(&m, &i, 3)) i = -1; - if (cl._coverage == -1) { + if (cl._coverage == SGMetarCloud::COVERAGE_NIL) { if (!scanBoundary(&m)) return false; if (i == -1) // 'VV///' @@ -791,7 +765,7 @@ bool SGMetar::scanTemperature() return false; if (!scanBoundary(&m)) { if (!strncmp(m, "XX", 2)) // not spec compliant! - m += 2, sign = 0; + m += 2, sign = 0, dew = temp; else { sign = 1; if (*m == 'M') @@ -907,19 +881,22 @@ bool SGMetar::scanRunwayReport() if (!strncmp(m, "CLRD", 4)) { m += 4; // runway cleared - r._deposit = "cleared"; + r._deposit_string = "cleared"; } else { if (scanNumber(&m, &i, 1)) { - r._deposit = runway_deposit[i]; + r._deposit = i; + r._deposit_string = runway_deposit[i]; } else if (*m == '/') m++; else return false; + if (*m == '1' || *m == '2' || *m == '5' || *m == '9') { // extent of deposit r._extent = *m - '0'; r._extent_string = runway_deposit_extent[*m - '0']; } else if (*m != '/') return false; + m++; i = -1; if (!strncmp(m, "//", 2)) @@ -928,7 +905,7 @@ bool SGMetar::scanRunwayReport() return false; if (i == 0) - r._depth = 0.5; // < 1 mm deep (let's say 0.5 :-) + r._depth = 0.0005; // < 1 mm deep (let's say 0.5 :-) else if (i > 0 && i <= 90) r._depth = i / 1000.0; // i mm deep else if (i >= 92 && i <= 98) @@ -954,6 +931,7 @@ bool SGMetar::scanRunwayReport() return false; _runways[id]._deposit = r._deposit; + _runways[id]._deposit_string = r._deposit_string; _runways[id]._extent = r._extent; _runways[id]._extent_string = r._extent_string; _runways[id]._depth = r._depth; @@ -1033,13 +1011,13 @@ bool SGMetar::scanTrendForecast() // (BLU|WHT|GRN|YLO|AMB|RED) static const struct Token colors[] = { - "BLU", "Blue", // 2500 ft, 8.0 km - "WHT", "White", // 1500 ft, 5.0 km - "GRN", "Green", // 700 ft, 3.7 km - "YLO", "Yellow", // 300 ft, 1.6 km - "AMB", "Amber", // 200 ft, 0.8 km - "RED", "Red", // <200 ft, <0.8 km - 0, 0 + { "BLU", "Blue" }, // 2500 ft, 8.0 km + { "WHT", "White" }, // 1500 ft, 5.0 km + { "GRN", "Green" }, // 700 ft, 3.7 km + { "YLO", "Yellow" }, // 300 ft, 1.6 km + { "AMB", "Amber" }, // 200 ft, 0.8 km + { "RED", "Red" }, // <200 ft, <0.8 km + { 0, 0 } }; @@ -1125,7 +1103,7 @@ const struct Token *SGMetar::scanToken(char **str, const struct Token *list) { const struct Token *longest = 0; int maxlen = 0, len; - char *s; + const char *s; for (int i = 0; (s = list[i].id); i++) { len = strlen(s); if (!strncmp(s, *str, len) && len > maxlen) { @@ -1137,4 +1115,40 @@ const struct Token *SGMetar::scanToken(char **str, const struct Token *list) return longest; } + +void SGMetarCloud::set(double alt, Coverage cov) +{ + _altitude = alt; + if (cov != -1) + _coverage = cov; +} + +SGMetarCloud::Coverage SGMetarCloud::getCoverage( const std::string & coverage ) +{ + if( coverage == "clear" ) return COVERAGE_CLEAR; + if( coverage == "few" ) return COVERAGE_FEW; + if( coverage == "scattered" ) return COVERAGE_SCATTERED; + if( coverage == "broken" ) return COVERAGE_BROKEN; + if( coverage == "overcast" ) return COVERAGE_OVERCAST; + return COVERAGE_NIL; +} + +const char * SGMetarCloud::COVERAGE_NIL_STRING = "nil"; +const char * SGMetarCloud::COVERAGE_CLEAR_STRING = "clear"; +const char * SGMetarCloud::COVERAGE_FEW_STRING = "few"; +const char * SGMetarCloud::COVERAGE_SCATTERED_STRING = "scattered"; +const char * SGMetarCloud::COVERAGE_BROKEN_STRING = "broken"; +const char * SGMetarCloud::COVERAGE_OVERCAST_STRING = "overcast"; + +void SGMetarVisibility::set(double dist, int dir, int mod, int tend) +{ + _distance = dist; + if (dir != -1) + _direction = dir; + if (mod != -1) + _modifier = mod; + if (tend != 1) + _tendency = tend; +} + #undef NaN