testRadiosondeDBProfileModel.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 <atm/RadiosondeDBProfileModel.h>

#include <fwk/CentralConfig.h>

#include <utl/TimeStamp.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 testRadiosondeDBProfileModel : public CppUnit::TestFixture{

  CPPUNIT_TEST_SUITE(testRadiosondeDBProfileModel);

  CPPUNIT_TEST(testHasData);
  CPPUNIT_TEST_EXCEPTION(testNoDataForTimeStamp, NoDataForModelException);

  CPPUNIT_TEST(testEvaluateDepthVsHeight);
  CPPUNIT_TEST(testEvaluateHeightVsDepth);
  CPPUNIT_TEST(testEvaluateTemperatureVsHeight);
  CPPUNIT_TEST(testEvaluateDensityVsHeight);
  CPPUNIT_TEST(testEvaluateRefractionIndexVsHeight);
  CPPUNIT_TEST(testNoDataMultipleCalls); // checks bug 456

  CPPUNIT_TEST_SUITE_END();

private:

public:

  void setUp(){
    //CentralConfig* theCC =
      CentralConfig::GetInstance(BOOTSTRAPFILE);
    ErrorLogger::GetInstance().SetVerbosity(Verbosity::eVerbose);
  }

  void tearDown(){
  }

  void testHasData() {

    Detector::GetInstance().Update(TimeStamp());
    RadiosondeDBProfileModel modelWithoutData;
    modelWithoutData.Init();
    CPPUNIT_ASSERT(Verify<Equal>(modelWithoutData.HasData(), false));

    Detector::GetInstance().Update(TimeStamp(1111112));
    RadiosondeDBProfileModel modelWithData;
    modelWithData.Init();
    CPPUNIT_ASSERT(Verify<Equal>(modelWithData.HasData(), true));

  }

  void testEvaluateDepthVsHeight() {

    Detector::GetInstance().Update(TimeStamp(1111112));
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();
    ProfileResult depthVsHeight = theAtm.EvaluateDepthVsHeight();

    CPPUNIT_ASSERT(Verify<CloseTo>( depthVsHeight.Y(500.*meter),   973.732*gram/cm/cm, 1e-3));
    //    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDepthError(500.*meter), 973.732*gram/cm/cm / 10., 1e-3));

    CPPUNIT_ASSERT(Verify<CloseTo>( depthVsHeight.Y(11000.*meter), 230.918*gram/cm/cm, 1e-3));
    //    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDepthError(11000.*meter), 230.918*gram/cm/cm / 10., 1e-3));

    CPPUNIT_ASSERT(Verify<CloseTo>( depthVsHeight.Y(26500.*meter), 20.3484*gram/cm/cm, 1e-3));
    //    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDepthError(26500.*meter), 20.3484*gram/cm/cm / 10., 1e-2));

    CPPUNIT_ASSERT(Verify<CloseTo>( depthVsHeight.Y(0.),           1033.93*gram/cm/cm, 1e-3));
//     CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDepthError(0.),      1033.93*gram/cm/cm / 10., 1e-3));


    CPPUNIT_ASSERT(Verify<CloseTo>( depthVsHeight.Y(500.*meter),
                                  depthVsHeight.Y( theAtm.EvaluateHeightVsDepth().Y( depthVsHeight.Y(500.*m))), 1e-3));

    ProfileResult heightVsDepth = theAtm.EvaluateHeightVsDepth();
    CPPUNIT_ASSERT(Verify<CloseTo>( depthVsHeight.MaxX(), heightVsDepth.Y(heightVsDepth.MinX()) ));
    CPPUNIT_ASSERT(Verify<CloseTo>( depthVsHeight.MinX(), heightVsDepth.Y(heightVsDepth.MaxX()) ));

  }

  void testEvaluateHeightVsDepth() {

    Detector::GetInstance().Update(TimeStamp(1111112));
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();
    ProfileResult heightVsDepth = theAtm.EvaluateHeightVsDepth();


    CPPUNIT_ASSERT(Verify<CloseTo>( heightVsDepth.Y(20.3484*gram/cm/cm),      26500.*meter, 1e-3));
//     CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateHeightError(20.3484*gram/cm/cm),
//                                theAtm.EvaluateHeight(20.3484*gram/cm/cm)/10., 1e-3));

    CPPUNIT_ASSERT(Verify<CloseTo>( heightVsDepth.Y(973.732*gram/cm/cm), 500.*meter,   1e-3));
//     CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateHeightError(973.732*gram/cm/cm),
//                                theAtm.EvaluateHeight(973.732*gram/cm/cm)/10.,   1e-3));

    CPPUNIT_ASSERT(Verify<CloseTo>( heightVsDepth.Y(230.918*gram/cm/cm), 11000.*meter, 1e-3));
//     CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateHeightError(230.918*gram/cm/cm),
//                                theAtm.EvaluateHeight(230.918*gram/cm/cm)/10.,   1e-3));

  }


  void testEvaluateTemperatureVsHeight() {

    Detector::GetInstance().Update(TimeStamp(1111112));
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();
    ProfileResult temperatureVsHeight = theAtm.EvaluateTemperatureVsHeight();

    CPPUNIT_ASSERT(Verify<CloseTo>( temperatureVsHeight.Y(0.),     288.15*kelvin, 1e-3));
    //    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateTemperatureError(0.),288.15*kelvin / 10., 1e-3));

    CPPUNIT_ASSERT(Verify<CloseTo>( temperatureVsHeight.Y(11500.), 216.65*kelvin, 1e-3));
    //    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateTemperatureError(11500.), 216.65*kelvin/10., 1e-3));

    CPPUNIT_ASSERT(Verify<CloseTo>( temperatureVsHeight.Y(29500.), 226.15*kelvin, 1e-3));
    //    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateTemperatureError(29500.), 226.15*kelvin/10., 1e-3));
  }


  void testEvaluateDensityVsHeight() {

    Detector::GetInstance().Update(TimeStamp(1111112));
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();
    ProfileResult densityVsHeight = theAtm.EvaluateDensityVsHeight();

    CPPUNIT_ASSERT(Verify<CloseTo>( densityVsHeight.Y(0.),           1.23e-3*g/cm/cm/cm));
    //    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDensityError(0.),      1.23e-3*g/cm/cm/cm / 10.));

    CPPUNIT_ASSERT(Verify<CloseTo>( densityVsHeight.Y(11500.*meter), 3.36458e-4*g/cm/cm/cm));
    //    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDensityError(11500.*meter), 3.36458e-4*g/cm/cm/cm / 10., 1.e-2));

    CPPUNIT_ASSERT(Verify<CloseTo>( densityVsHeight.Y(29000.*meter), 2.10e-5*g/cm/cm/cm));
    //    CPPUNIT_ASSERT(Verify<CloseTo>( theAtm.EvaluateDensityError(29000.*meter), 2.10e-5*g/cm/cm/cm / 10.));

  }

  void testEvaluateRefractionIndexVsHeight() {

    Detector::GetInstance().Update(TimeStamp(1111112));
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();
    ProfileResult riVsHeight = theAtm.EvaluateRefractionIndexVsHeight();

    CPPUNIT_ASSERT(Verify<CloseTo>( riVsHeight.Y(0.), 1.00029, 1e-3));
  }

  void testNoDataForTimeStamp() {
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();
    Detector::GetInstance().Update(TimeStamp());
    theAtm.EvaluateTemperatureVsHeight();
  }

  void testNoDataMultipleCalls() {
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();
    Detector::GetInstance().Update(TimeStamp(123));

    int numExceptions = 0;
    for (int i=0 ; i<5 ; i++) {
      try {
        theAtm.EvaluateTemperatureVsHeight();
      } catch (NoDataForModelException& ex) {
        numExceptions++;
      }
    }
    CPPUNIT_ASSERT(Verify<Equal>(numExceptions, 5));
  }
};

CPPUNIT_TEST_SUITE_REGISTRATION(testRadiosondeDBProfileModel);


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