1 // magvar.cxx -- magnetic variation wrapper class
3 // Written by Curtis Olson, started July 2000.
5 // Copyright (C) 2000 Curtis L. Olson - http://www.flightgear.org/~curt
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 # include <simgear_config.h>
31 #include <simgear/magvar/magvar.hxx>
32 #include <simgear/math/SGMath.hxx>
34 #include "coremag.hxx"
44 SGMagVar::~SGMagVar() {
48 void SGMagVar::update( double lon, double lat, double alt_m, double jd ) {
49 // Calculate local magnetic variation
51 // cout << "alt_m = " << alt_m << endl;
52 magvar = calc_magvar( lat, lon, alt_m / 1000.0, (long)jd, field );
53 magdip = atan(field[5]/sqrt(field[3]*field[3]+field[4]*field[4]));
56 void SGMagVar::update( const SGGeod& geod, double jd ) {
58 update(geod.getLongitudeRad(), geod.getLatitudeRad(),
59 geod.getElevationM(), jd);
63 double sgGetMagVar( double lon, double lat, double alt_m, double jd ) {
64 // cout << "lat = " << lat << " lon = " << lon << " elev = " << alt_m
65 // << " JD = " << jd << endl;
68 return calc_magvar( lat, lon, alt_m / 1000.0, (long)jd, field );
71 double sgGetMagVar( const SGGeod& pos, double jd )
73 return sgGetMagVar(pos.getLongitudeRad(), pos.getLatitudeRad(),
74 pos.getElevationM(), jd);