3 * Magnetic variation wrapper class.
6 // Written by Curtis Olson, started July 2000.
8 // Copyright (C) 2000 Curtis L. Olson - curt@flightgear.org
10 // This program is free software; you can redistribute it and/or
11 // modify it under the terms of the GNU General Public License as
12 // published by the Free Software Foundation; either version 2 of the
13 // License, or (at your option) any later version.
15 // This program is distributed in the hope that it will be useful, but
16 // WITHOUT ANY WARRANTY; without even the implied warranty of
17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
18 // General Public License for more details.
20 // You should have received a copy of the GNU General Public License
21 // along with this program; if not, write to the Free Software
22 // Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
32 # error This library requires C++
37 * Magnetic variation wrapper class.
39 * The SGMagVar class calculates the magnetic variation and dip for
40 * any position, altitude, and time. It is a complete
41 * re-implimentation of the NIMA WMM 2000 (not derived from their demo
44 * The SGMagVar class is really a simple wrapper around the core Ed
45 * Williams code which does all the hard work. This class allows you
46 * to crunch the math once and then do multiple polls of the
47 * data. However, if your position, altitude, or time has changed
48 * significantly, you should call the update() method to recrunch new
61 * This creates an instance of the SGMagVar object.
62 * You must call the update() method before any queries will be valid.
69 /** Recalculate the magnetic offset and dip.
70 * The update() method requires you to pass in your position and
71 * the julian date. Lon and lat are specified in radians, altitude
72 * is specified in meters. Julian date can be conveniently
73 * calculated by using an instance of the SGTime class.
74 * @param lon longitude in radians
75 * @param lat latitude in radians
76 * @param alt_m altitude above sea level in meters
77 * @param jd julian date
79 void update( double lon, double lat, double alt_m, double jd );
81 /** @return the current magnetic variation in radians. */
82 double get_magvar() const { return magvar; }
84 /** @return the current magnetic dip in radians. */
85 double get_magdip() const { return magdip; }
91 * Lookup the magvar for any arbitrary location (This function doesn't
92 * save state like the SGMagVar class. This function triggers a fair
93 * amount of CPU work, so use it cautiously.
94 * @return the magvar in radians
96 double sgGetMagVar( double lon, double lat, double alt_m, double jd );