1 // metar interface class
3 // Written by Melchior FRANZ, started January 2005.
5 // Copyright (C) 2005 Melchior FRANZ - mfranz@aon.at
7 // This program is free software; you can redistribute it and/or
8 // modify it under the terms of the GNU General Public License as
9 // published by the Free Software Foundation; either version 2 of the
10 // License, or (at your option) any later version.
12 // This program is distributed in the hope that it will be useful, but
13 // WITHOUT ANY WARRANTY; without even the implied warranty of
14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 // General Public License for more details.
17 // You should have received a copy of the GNU General Public License
18 // along with this program; if not, write to the Free Software
19 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
25 * Implements FGMetar class that inherits from SGMetar.
27 * o provide defaults for unset values
28 * o interpolate/randomize data (GREATER_THAN)
29 * o derive additional values (time, age, snow cover)
30 * o consider minimum identifier (CAVOK, mil. color codes)
33 * - NSC & mil. color codes
36 #include <simgear/math/sg_random.h>
37 #include <simgear/timing/sg_time.hxx>
38 #include <Main/fg_props.hxx>
40 #include "fgmetar.hxx"
43 FGMetar::FGMetar(const string& icao, const string& proxy, const string& port, const string& auth) :
44 SGMetar(icao, proxy, port, auth, _rq_time = globals->get_time_params()->get_cur_time()),
50 // CAVOK: visibility >= 10km; lowest cloud layer >= 5000 ft; any coverage
52 if (_min_visibility.getVisibility_m() == SGMetarNaN)
53 _min_visibility.set(12000.0);
55 if (_max_visibility.getVisibility_m() == SGMetarNaN)
56 _min_visibility.set(12000.0);
58 vector<SGMetarCloud> cv = _clouds;;
61 cl.set(5500 * SG_FEET_TO_METER, 2);
62 _clouds.push_back(cl);
67 d = _min_visibility.getVisibility_m();
70 if (_min_visibility.getModifier() == SGMetarVisibility::GREATER_THAN)
71 d += 15000.0;// * sg_random();
72 _min_visibility.set(d);
74 if (_max_visibility.getVisibility_m() == SGMetarNaN)
75 _max_visibility.set(d);
77 for (i = 0; i < 8; i++) {
78 d = _dir_visibility[i].getVisibility_m();
80 _dir_visibility[i].set(10000.0);
81 if (_dir_visibility[i].getModifier() == SGMetarVisibility::GREATER_THAN)
82 d += 15000.0;// * sg_random();
83 _dir_visibility[i].set(d);
87 if (_wind_dir == -1) {
88 if (_wind_range_from == -1) {
93 _wind_dir = (_wind_range_from + _wind_range_to) / 2;
95 } else if (_wind_range_from == -1) {
96 _wind_range_from = _wind_range_to = _wind_dir;
99 if (_gust_speed == SGMetarNaN)
103 vector<SGMetarCloud> cv = _clouds;
104 vector<SGMetarCloud>::iterator cloud;
106 for (i = 0, cloud = cv.begin(); cloud != cv.end(); cloud++, i++) {
107 int cov = cloud->getCoverage();
111 double alt = cloud->getAltitude_ft();
112 if (alt == SGMetarNaN)
115 cloud->set(alt, cov);
119 // temperature/pressure
120 if (_temp == SGMetarNaN)
123 if (_dewp == SGMetarNaN)
126 if (_pressure == SGMetarNaN)
127 _pressure = 30.0 * SG_INHG_TO_PA;
130 map<string, SGMetarRunway> rm = getRunways();
131 map<string, SGMetarRunway>::const_iterator runway;
132 for (runway = rm.begin(); runway != rm.end(); runway++) {
133 SGMetarRunway rwy = runway->second;
134 if (rwy.getDeposit() >= 3 ) {
139 if (_temp < 5.0 && _snow)
141 if (_temp < 1.0 && getRelHumidity() > 80)
144 _time = sgTimeGetGMT(_year - 1900, _month - 1, _day, _hour, _minute, 0);
146 SG_LOG(SG_GENERAL, SG_INFO, _data);
148 SG_LOG(SG_GENERAL, SG_INFO, "METAR from proxy");
150 SG_LOG(SG_GENERAL, SG_INFO, "METAR from weather.noaa.gov");
154 long FGMetar::getAge_min() const
156 time_t now = _x_proxy ? _rq_time : time(0);
157 return (now - _time) / 60;