1 /*****************************************************************************
3 Module: FGLocalWeatherDatabase.cpp
4 Author: Christian Mayer
6 Called by: main program
8 -------- Copyright (C) 1999 Christian Mayer (fgfs@christianmayer.de) --------
10 This program is free software; you can redistribute it and/or modify it under
11 the terms of the GNU General Public License as published by the Free Software
12 Foundation; either version 2 of the License, or (at your option) any later
15 This program is distributed in the hope that it will be useful, but WITHOUT
16 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
17 FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
20 You should have received a copy of the GNU General Public License along with
21 this program; if not, write to the Free Software Foundation, Inc., 59 Temple
22 Place - Suite 330, Boston, MA 02111-1307, USA.
24 Further information about the GNU General Public License can also be found on
25 the world wide web at http://www.gnu.org.
27 FUNCTIONAL DESCRIPTION
28 ------------------------------------------------------------------------------
29 Database for the local weather
30 This database is the only one that gets called from FG
33 ------------------------------------------------------------------------------
34 28.05.1999 Christian Mayer Created
35 16.06.1999 Durk Talsma Portability for Linux
36 20.06.1999 Christian Mayer added lots of consts
37 11.10.1999 Christian Mayer changed set<> to map<> on Bernie Bright's
39 19.10.1999 Christian Mayer change to use PLIB's sg instead of Point[2/3]D
40 and lots of wee code cleaning
41 14.12.1999 Christian Mayer Changed the internal structure to use Dave
42 Eberly's spherical interpolation code. This
43 stops our dependancy on the (ugly) voronoi
44 code and simplyfies the code structure a lot.
45 *****************************************************************************/
47 /****************************************************************************/
49 /****************************************************************************/
50 #include <Include/compiler.h>
51 #include <Include/fg_constants.h>
53 #include <Aircraft/aircraft.hxx>
55 #include "FGLocalWeatherDatabase.h"
57 #include "FGWeatherParse.h"
59 /****************************************************************************/
60 /********************************** CODE ************************************/
61 /****************************************************************************/
63 FGLocalWeatherDatabase* FGLocalWeatherDatabase::theFGLocalWeatherDatabase = 0;
64 FGLocalWeatherDatabase *WeatherDatabase;
66 void FGLocalWeatherDatabase::init(const WeatherPrecision visibility, const DatabaseWorkingType type)
68 cerr << "Initializing FGLocalWeatherDatabase\n";
69 cerr << "-----------------------------------\n";
71 if (theFGLocalWeatherDatabase)
73 cerr << "Error: only one local weather allowed";
77 setWeatherVisibility(visibility);
79 DatabaseStatus = type;
80 database = 0; //just get sure...
83 //I don't need to set theThunderstorm as Thunderstorm == false
85 switch(DatabaseStatus)
89 cerr << "Error: there's no global database anymore!\n";
96 FGWeatherParse *parsed_data = new FGWeatherParse();
98 parsed_data->input( "weather/current.gz" );
99 unsigned int n = parsed_data->stored_stations();
101 sgVec2 *p = new sgVec2 [n];
102 FGPhysicalProperties *f = new FGPhysicalProperties[n];
105 for (unsigned int i = 0; i < n; i++)
107 f[i] = parsed_data->getFGPhysicalProperties(i);
108 parsed_data->getPosition(i, p[i]);
114 // free the memory of the parsed data to ease the required memory
115 // for the very memory consuming spherical interpolation
118 //and finally init the interpolation
119 cerr << "\nInitialiating Interpolation. (2-3 minutes on a PII-350)\n";
120 database = new SphereInterpolate<FGPhysicalProperties>(n, p, f);
122 //and free my allocations:
125 cerr << "Finished weather init.\n";
131 cerr << "FGLocalWeatherDatabase error: Distant database isn't implemented yet!\n";
132 cerr << " using random mode instead!\n";
137 double x[2] = {0.0, 0.0}; //make an standard weather that's the same at the whole world
138 double y[2] = {0.0, 0.0}; //make an standard weather that's the same at the whole world
139 double z[2] = {1.0, -1.0}; //make an standard weather that's the same at the whole world
140 FGPhysicalProperties f[2]; //make an standard weather that's the same at the whole world
141 database = new SphereInterpolate<FGPhysicalProperties>(2,x,y,z,f);
146 cerr << "FGLocalWeatherDatabase error: Unknown database type specified!\n";
150 FGLocalWeatherDatabase::~FGLocalWeatherDatabase()
156 /****************************************************************************/
157 /* reset the whole database */
158 /****************************************************************************/
159 void FGLocalWeatherDatabase::reset(const DatabaseWorkingType type)
161 cerr << "FGLocalWeatherDatabase::reset isn't supported yet\n";
164 /****************************************************************************/
165 /* update the database. Since the last call we had dt seconds */
166 /****************************************************************************/
167 void FGLocalWeatherDatabase::update(const WeatherPrecision dt)
169 //if (DatabaseStatus==use_global)
170 // global->update(dt);
173 void FGLocalWeatherDatabase::update(const sgVec3& p) //position has changed
175 sgCopyVec3(last_known_position, p);
177 //uncomment this when you are using the GlobalDatabase
179 cerr << "****\nupdate(p) inside\n";
180 cerr << "Parameter: " << p[0] << "/" << p[1] << "/" << p[2] << "\n";
182 sgSetVec2( p2d, p[0], p[1] );
183 cerr << FGPhysicalProperties2D(get(p2d), p2d);
189 void FGLocalWeatherDatabase::update(const sgVec3& p, const WeatherPrecision dt) //time and/or position has changed
191 sgCopyVec3(last_known_position, p);
194 /****************************************************************************/
195 /* Get the physical properties on the specified point p out of the database */
196 /****************************************************************************/
197 FGPhysicalProperty FGLocalWeatherDatabase::get(const sgVec3& p) const
199 return FGPhysicalProperty(database->Evaluate(p), p[3]);
203 /* fix a problem with mw compilers in that they don't know the
204 difference between the next two methods. Since the first one
205 doesn't seem to be used anywhere, I commented it out. This is
206 supposed to be fixed in the forthcoming CodeWarrior Release
209 FGPhysicalProperties FGLocalWeatherDatabase::get(const sgVec2& p) const
211 return database->Evaluate(p);
215 WeatherPrecision FGLocalWeatherDatabase::getAirDensity(const sgVec3& p) const
217 FGPhysicalProperty dummy(database->Evaluate(p), p[3]);
220 (dummy.AirPressure*FG_WEATHER_DEFAULT_AIRDENSITY*FG_WEATHER_DEFAULT_TEMPERATURE) /
221 (dummy.Temperature*FG_WEATHER_DEFAULT_AIRPRESSURE);
225 void FGLocalWeatherDatabase::setSnowRainIntensity(const WeatherPrecision x, const sgVec2& p)
227 /* not supported yet */
230 void FGLocalWeatherDatabase::setSnowRainType(const SnowRainType x, const sgVec2& p)
232 /* not supported yet */
235 void FGLocalWeatherDatabase::setLightningProbability(const WeatherPrecision x, const sgVec2& p)
237 /* not supported yet */
240 void FGLocalWeatherDatabase::setProperties(const FGPhysicalProperties2D& x)
242 /* not supported yet */
245 void fgUpdateWeatherDatabase(void)
249 current_aircraft.fdm_state->get_Latitude(),
250 current_aircraft.fdm_state->get_Longitude(),
251 current_aircraft.fdm_state->get_Altitude() * FEET_TO_METER);
253 WeatherDatabase->update( position );