simgear/magvar/magvar.cxx

   1 // magvar.cxx -- magnetic variation wrapper class
   2 //
   3 // Written by Curtis Olson, started July 2000.
   4 //
   5 // Copyright (C) 2000  Curtis L. Olson  - http://www.flightgear.org/~curt
   6 //
   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.
  11 //
  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.
  16 //
  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.
  20 //
  21 // $Id$
  22
  23
  24 #ifdef HAVE_CONFIG_H
  25 #  include <simgear_config.h>
  26 #endif
  27
  28
  29 #include <cmath>
  30
  31 #include <simgear/magvar/magvar.hxx>
  32 #include <simgear/math/SGMath.hxx>
  33
  34 #include "coremag.hxx"
  35 #include "magvar.hxx"
  36
  37
  38 SGMagVar::SGMagVar()
  39   : magvar(0.0),
  40     magdip(0.0)
  41 {
  42 }
  43
  44 SGMagVar::~SGMagVar() {
  45 }
  46
  47
  48 void SGMagVar::update( double lon, double lat, double alt_m, double jd ) {
  49     // Calculate local magnetic variation
  50     double field[6];
  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]));
  54 }
  55
  56 void SGMagVar::update( const SGGeod& geod, double jd ) {
  57
  58   update(geod.getLongitudeRad(), geod.getLatitudeRad(),
  59     geod.getElevationM(), jd);
  60 }
  61
  62
  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;
  66
  67     double field[6];
  68     return calc_magvar( lat, lon, alt_m / 1000.0, (long)jd, field );
  69 }
  70
  71 double sgGetMagVar( const SGGeod& pos, double jd )
  72 {
  73   return sgGetMagVar(pos.getLongitudeRad(), pos.getLatitudeRad(),
  74     pos.getElevationM(), jd);
  75 }
  76