ReadMolecularHL.cc

Go to the documentation of this file.

#include <sstream>

#include "ReadMolecularHL.h"
#include <utl/ErrorLogger.h>

#include <evt/Event.h>
#include <det/Detector.h>
#include <utl/TimeStamp.h>
#include <utl/UTCDateTime.h>

#include <atm/AerosolDB.h>
#include <atm/AerosolZone.h>
#include <atm/ProfileResult.h>

#include <atm/MonthlyAvgDBProfileModel.h>

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

using namespace ReadMolecularHLNS;

ReadMolecularHL::ReadMolecularHL(){}

ReadMolecularHL::~ReadMolecularHL(){}

VModule::ResultFlag ReadMolecularHL::Init(){
  //
  //
  INFO("ReadMolecularHL::Init()");
  return eSuccess;
}

VModule::ResultFlag ReadMolecularHL::Run(evt::Event& event){
  //
  //
  INFO("ReadMolecularHL::Run()");

    Detector& det = Detector::GetInstance();

    // Set detector time to something
    //
    det.Update(UTCDateTime(2005,3,21,15,15).GetTimeStamp());        

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

    ProfileResult depthVsHeight   = theAtm.EvaluateDepthVsHeight();
    ProfileResult densityVsHeight = theAtm.EvaluateDensityVsHeight();
    ProfileResult tempVsHeight    = theAtm.EvaluateTemperatureVsHeight();

    // dump out some profile parameters every 300m
    //
    double minHeight = 
      max( tempVsHeight.MinX(), 
           max( depthVsHeight.MinX(), 
                densityVsHeight.MinX()));
    double maxHeight =
      min( tempVsHeight.MaxX(),
           max( depthVsHeight.MaxX(),
                densityVsHeight.MaxX()));
    double deltaH = 300*m;

    for (double height = minHeight; height < maxHeight; height += deltaH) {
      cout << "height     = " << height/m << " m " << endl;
      cout << "  depth    = " << depthVsHeight.Y(height)/(g/cm/cm) << " g/cm^2" << endl;
      cout << "  density  = " << densityVsHeight.Y(height)/(g/cm/cm/cm) << " g/cm^3" << endl;
      cout << "  temp     = " << tempVsHeight.Y(height)/kelvin  << " deg k " << endl;
    }

  return eSuccess;
}

VModule::ResultFlag ReadMolecularHL::Finish(){
  //
  //
  INFO("ReadMolecularHL::Finish()");
  return eSuccess;
}