//
// Written by Curtis Olson, started July 2000.
//
-// Copyright (C) 2000 Curtis L. Olson - curt@flightgear.org
+// Copyright (C) 2000 Curtis L. Olson - http://www.flightgear.org/~curt
//
// This program is free software; you can redistribute it and/or
// modify it under the terms of the GNU General Public License as
//
// You should have received a copy of the GNU General Public License
// along with this program; if not, write to the Free Software
-// Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
+// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
//
// $Id$
-#ifdef HAVE_CONFIG_H
-# include <config.h>
-#endif
-
#include <math.h>
-#include <simgear/misc/fgpath.hxx>
#include <simgear/magvar/magvar.hxx>
+#include <simgear/math/SGMath.hxx>
#include "coremag.hxx"
#include "magvar.hxx"
-SGMagVar::SGMagVar() {
+SGMagVar::SGMagVar()
+ : magvar(0.0),
+ magdip(0.0)
+{
}
SGMagVar::~SGMagVar() {
magdip = atan(field[5]/sqrt(field[3]*field[3]+field[4]*field[4]));
}
+void SGMagVar::update( const SGGeod& geod, double jd ) {
+
+ update(geod.getLongitudeRad(), geod.getLatitudeRad(),
+ geod.getElevationM(), jd);
+}
+
double sgGetMagVar( double lon, double lat, double alt_m, double jd ) {
+ // cout << "lat = " << lat << " lon = " << lon << " elev = " << alt_m
+ // << " JD = " << jd << endl;
+
double field[6];
return calc_magvar( lat, lon, alt_m / 1000.0, (long)jd, field );
}
+
+double sgGetMagVar( const SGGeod& pos, double jd )
+{
+ return sgGetMagVar(pos.getLongitudeRad(), pos.getLatitudeRad(),
+ pos.getElevationM(), jd);
+}
+