testMagneticFieldModel.cc

Go to the documentation of this file.

#include <iostream>

#include <tst/Verify.h>

#include <fwk/MagneticFieldModel.h>
#include <utl/AugerException.h>
#include <utl/AugerUnits.h>
#include <utl/CoordinateSystemPtr.h>
#include <utl/Point.h>
#include <utl/UTMPoint.h>
#include <utl/UTCDateTime.h>
#include <fwk/CentralConfig.h>

#include <det/Detector.h>


#include <cppunit/extensions/HelperMacros.h>

using namespace utl;
using namespace tst;

class testMagneticFieldModel : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(testMagneticFieldModel);
  CPPUNIT_TEST(testIGRFValues);
//CPPUNIT_TEST(testEqualAndCopyCtor);
//CPPUNIT_TEST_EXCEPTION(testBinException, InvalidTraceBoundException);
  CPPUNIT_TEST_SUITE_END();

private:

public:
  void
  setUp()
  {

    // CentralConfig must be created first.
    //
    /*CentralConfig* fCentralConfig =*/
    fwk::CentralConfig::GetInstance(BOOTSTRAPFILE);

    // create a detector, needed for access to reference ellipsoid
    det::Detector::GetInstance();

  }

  void
  tearDown()
  {
  }

/*
  void
  testEqualAndCopyCtor()
  {
    ChannelMagneticFieldModel copyMagneticFieldModel(ourChannelContainer);
    CPPUNIT_ASSERT(copyMagneticFieldModel == ourChannelContainer);
    CPPUNIT_ASSERT(!(copyMagneticFieldModel != ourChannelContainer));
  }
*/

  void
  testIGRFValues()
  {
    // ask the magnetic field model for data at a given time and compare it with values determined from the NOAA webpage

    //~ these are the data delivered by the webpage
    //~ Latitude:         35° 28' 32" S
    //~ Longitude:         69° 35' 12" W
    //~ Elevation:         1.5 K
    //~ Date               Declination(+E|-W)         Inclination (+D|-U)         HorizontalIntensity  NorthComp(+N|-S)         EastComp(+E|-W)         VerticalComp(+D|-U)         TotalField
    //~ 2014-07-03         2° 26' 10"                 -36° 17' 27"                 19,506.4 nT                19,488.8 nT        829.2 nT          -14,324.0 nT          24,200.7 nT
    //~ Change/year -8.7'               -6.1'                 -73.6 nT             -71.5 nT           -52.4 nT          0.6 nT                -59.7 nT

    // configure location (Malargue, 1500m asl)
    const double Lat = -(35. + 28./60. + 32./3600.) * degree;
    const double Long = -(69. + 35./60. + 12./3600.) * degree;
    const double Height = 1500. * meter;

//  const UTMPoint utmPosition(Lat,Long,Height,ReferenceEllipsoid::eWGS84);
//  const utl::Point localPosition = utmPosition.GetPoint();

    const utl::Point localPosition = ReferenceEllipsoid::Get(ReferenceEllipsoid::eWGS84).LatitudeLongitudeHeightToPoint(Lat, Long, Height);

    // configure time (2014-07-03 at noon)
    const TimeStamp theTime = UTCDateTime(2014, 7, 3, 12, 0, 0).GetTimeStamp();

    // get vector and compare components with NOAA web calculator
    const utl::Vector theMagneticFieldVector = fwk::MagneticFieldModel::instance().GetMagneticFieldVector(localPosition, theTime);
    const utl::CoordinateSystemPtr localCS = fwk::LocalCoordinateSystem::Create(localPosition);

    CPPUNIT_ASSERT(Verify<CloseTo>(theMagneticFieldVector.GetX(localCS),829.2*nano*tesla,0.005)); // test East component within 0.5%
    CPPUNIT_ASSERT(Verify<CloseTo>(theMagneticFieldVector.GetY(localCS),19488.8*nano*tesla,0.005)); // test North component within 0.5%
    CPPUNIT_ASSERT(Verify<CloseTo>(theMagneticFieldVector.GetZ(localCS),14324.0*nano*tesla,0.005)); // test Up component within 0.5%

    // get declination and compare
    CPPUNIT_ASSERT(Verify<CloseTo>(fwk::MagneticFieldModel::instance().GetDeclination(localPosition, theTime),(2+26./60.+10./3600.)*degree,0.005)); // test Declination within 0.5%

    // get inclination and compare
    CPPUNIT_ASSERT(Verify<CloseTo>(fwk::MagneticFieldModel::instance().GetInclination(localPosition, theTime),-(36+17./60.+27./3600.)*degree,0.005)); // test Inclination within 0.5%

    // get strength and compare
    CPPUNIT_ASSERT(Verify<CloseTo>(fwk::MagneticFieldModel::instance().GetStrength(localPosition, theTime),24200.7*nano*tesla,0.005)); // test field strength within 0.5%

  }


};

CPPUNIT_TEST_SUITE_REGISTRATION(testMagneticFieldModel);

// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// End: