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
39 #include <simgear/math/sg_random.h>
40 #include <simgear/timing/sg_time.hxx>
41 #include <Main/fg_props.hxx>
43 #include "fgmetar.hxx"
46 FGMetar::FGMetar(const string& icao, const string& proxy, const string& port, const string& auth) :
47 SGMetar(icao, proxy, port, auth, _rq_time = globals->get_time_params()->get_cur_time()),
53 // CAVOK: visibility >= 10km; lowest cloud layer >= 5000 ft; any coverage
55 if (_min_visibility.getVisibility_m() == SGMetarNaN)
56 _min_visibility.set(12000.0);
58 if (_max_visibility.getVisibility_m() == SGMetarNaN)
59 _min_visibility.set(12000.0);
61 vector<SGMetarCloud> cv = _clouds;;
64 cl.set(5500 * SG_FEET_TO_METER, 2);
65 _clouds.push_back(cl);
70 d = _min_visibility.getVisibility_m();
73 if (_min_visibility.getModifier() == SGMetarVisibility::GREATER_THAN)
74 d += 15000.0;// * sg_random();
75 _min_visibility.set(d);
77 if (_max_visibility.getVisibility_m() == SGMetarNaN)
78 _max_visibility.set(d);
80 for (i = 0; i < 8; i++) {
81 d = _dir_visibility[i].getVisibility_m();
83 _dir_visibility[i].set(10000.0);
84 if (_dir_visibility[i].getModifier() == SGMetarVisibility::GREATER_THAN)
85 d += 15000.0;// * sg_random();
86 _dir_visibility[i].set(d);
90 if (_wind_dir == -1) {
91 if (_wind_range_from == -1) {
96 _wind_dir = (_wind_range_from + _wind_range_to) / 2;
98 } else if (_wind_range_from == -1) {
99 _wind_range_from = _wind_range_to = _wind_dir;
102 if (_wind_speed == SGMetarNaN)
104 if (_gust_speed == SGMetarNaN)
108 vector<SGMetarCloud> cv = _clouds;
109 vector<SGMetarCloud>::iterator cloud, cv_end = cv.end();
111 for (i = 0, cloud = cv.begin(); cloud != cv_end; ++cloud, i++) {
112 int cov = cloud->getCoverage();
116 double alt = cloud->getAltitude_ft();
117 if (alt == SGMetarNaN)
120 cloud->set(alt, cov);
124 // temperature/pressure
125 if (_temp == SGMetarNaN)
128 if (_dewp == SGMetarNaN)
131 if (_pressure == SGMetarNaN)
132 _pressure = 30.0 * SG_INHG_TO_PA;
135 map<string, SGMetarRunway> rm = getRunways();
136 map<string, SGMetarRunway>::const_iterator runway, rm_end = rm.end();
137 for (runway = rm.begin(); runway != rm_end; ++runway) {
138 SGMetarRunway rwy = runway->second;
139 if (rwy.getDeposit() >= 3 ) {
144 if (_temp < 5.0 && _snow)
146 if (_temp < 1.0 && getRelHumidity() > 80)
149 _time = sgTimeGetGMT(_year - 1900, _month - 1, _day, _hour, _minute, 0);
151 SG_LOG(SG_GENERAL, SG_INFO, _data);
153 SG_LOG(SG_GENERAL, SG_INFO, "METAR from proxy");
155 SG_LOG(SG_GENERAL, SG_INFO, "METAR from weather.noaa.gov");
159 long FGMetar::getAge_min() const
161 time_t now = _x_proxy ? _rq_time : time(0);
162 return (now - _time) / 60;