testParametricXMLProfileModel.cc

Go to the documentation of this file.

#include <det/Detector.h>
#include <atm/AerosolDB.h>
#include <atm/AerosolZone.h>
#include <atm/AttSlice.h>
#include <atm/PFSlice.h>
#include <atm/MolecularZone.h>
#include <atm/ProfileResult.h>

#include <fwk/CentralConfig.h>

#include <utl/TimeStamp.h>
#include <utl/UTCDateTime.h>
#include <utl/ErrorLogger.h>

#include <cppunit/extensions/HelperMacros.h>
#include <tst/Verify.h>

#include <iostream>

using namespace std;
using namespace det;
using namespace atm;
using namespace fwk;
using namespace utl;
using namespace tst;


class testParametricXMLProfileModel : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(testParametricXMLProfileModel);

  CPPUNIT_TEST(testEvaluateDepth);
  CPPUNIT_TEST(testEvaluateHeight);
  CPPUNIT_TEST(testEvaluateTemperature);
  CPPUNIT_TEST(testEvaluateDensity);
  CPPUNIT_TEST(testEvaluateRefractionIndex);
  //  CPPUNIT_TEST(testMinMax);

  CPPUNIT_TEST_SUITE_END();

public:

  void setUp(){

    CentralConfig::GetInstance(BOOTSTRAPFILE);
    ErrorLogger::GetInstance().SetVerbosity(Verbosity::eVerbose);

  }

  void tearDown(){
  }

  void testEvaluateDepth() {

    Detector::GetInstance().Update(UTCDateTime(2005,1,1,0,0,0).GetTimeStamp());
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();

    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDepthVsHeight().Y(500.*meter),   972.868*gram/cm/cm, 1e-3));
    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDepthVsHeight().Y(11000.*meter), 230.918*gram/cm/cm, 1e-3));
    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDepthVsHeight().Y(26500.*meter), 20.3484*gram/cm/cm, 1e-3));
    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDepthVsHeight().Y(0.),           1033.93*gram/cm/cm, 1e-3));
    CPPUNIT_ASSERT(Verify<CloseTo> (theAtm.EvaluateDepthVsHeight().Y(500.*meter),
                                  theAtm.EvaluateDepthVsHeight().Y( theAtm.EvaluateHeightVsDepth().Y( theAtm.EvaluateDepthVsHeight().Y(500.*m))), 1e-3));

  }

  void testEvaluateHeight() {

    Detector::GetInstance().Update(UTCDateTime(2005,1,1,0,0,0).GetTimeStamp());

    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();

    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateHeightVsDepth().Y(20.3484*gram/cm/cm), 26500.*meter, 1e-3));
    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateHeightVsDepth().Y(972.858*gram/cm/cm), 500.*meter,   1e-3));
    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateHeightVsDepth().Y(230.918*gram/cm/cm), 11000.*meter, 1e-3));


  }

  void testEvaluateTemperature() {

    Detector::GetInstance().Update(UTCDateTime(2005,1,1,0,0,0).GetTimeStamp());

    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();

    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateTemperatureVsHeight().Y(0.),     288.15*kelvin, 1e-3));
    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateTemperatureVsHeight().Y(11500.*m), 216.65*kelvin, 1e-3));
    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateTemperatureVsHeight().Y(29500.*m), 226.15*kelvin, 1e-3));

  }

  void testEvaluateDensity(){

    Detector::GetInstance().Update(UTCDateTime(2005,1,1,0,0,0).GetTimeStamp());

    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();

    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDensityVsHeight().Y(0.),           1.23e-3*g/cm/cm/cm));
    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDensityVsHeight().Y(11500.*meter), 3.36458e-4*g/cm/cm/cm));
    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDensityVsHeight().Y(29000.*meter), 2.10e-5*g/cm/cm/cm));

  }

  void testEvaluateRefractionIndex(){

    Detector::GetInstance().Update(UTCDateTime(2005,1,1,0,0,0).GetTimeStamp());

    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();

    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateRefractionIndexVsHeight().Y(0.)         , 1.00029, 1e-3));
    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateRefractionIndexVsHeight().Y(11500*meter), 1.00006, 1e-3));
    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateRefractionIndexVsHeight().Y(29000*meter), 1.,      1e-3));
  }

//   void testMinMax() {

//     Detector::GetInstance().Update(UTCDateTime(2005,1,1,0,0,0).GetTimeStamp());
//     const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();

//     cout << "theAtm.GetMinDepth " << theAtm.GetMinDepth()/(g/cm/cm) << endl;

//   }

};

CPPUNIT_TEST_SUITE_REGISTRATION(testParametricXMLProfileModel);



// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C .. -k"
// End: