ReadMolecularLL.cc

Go to the documentation of this file.

#include <sstream>

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

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

#include <atm/MolecularDB.h>
#include <atm/MolecularZone.h>

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

using namespace ReadMolecularLLNS;

ReadMolecularLL::ReadMolecularLL(){}

ReadMolecularLL::~ReadMolecularLL(){}

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

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

    Detector& det = Detector::GetInstance();

    // Retrieve the low-level molecular DB interface, specifying, in this case,
    // that you want to retrieve one of the Malarguemonthly average profiles.
    //
    const MolecularDB& mdb = det.GetAtmosphere().GetMolecularDB(MolecularDB::eMalargueJanuaryModel);

    // Iterate over all the zones found in the db
    //
    cout << "Malargue montlhy model for January : " << endl;
    for (MolecularDB::ZoneIterator zIt = mdb.ZonesBegin() ;
         zIt != mdb.ZonesEnd() ; ++zIt){

      cout << "Got zone : " << zIt->GetName() << endl;
      cout << "  with coordinates : " << zIt->GetNorthing()/m << "N, " << zIt->GetEasting()/m << " E " << endl; 

      // Iterate over slices in the zone
      //
      for (MolecularZone::LayerIterator sIt = zIt->LayersBegin() ;
           sIt != zIt->LayersEnd() ; ++sIt){
 
        cout << endl;
        cout << "height      " << sIt->GetHeight()/m             << " m " << endl;
        cout << "depth       " << sIt->GetDepth()/(gram/cm/cm)   << " g/cm^2" << endl;
        cout << "pressure    " << sIt->GetPressure()/pascal      << " pascal " << endl;
        cout << "temperature " << sIt->GetTemperature()/kelvin   << " kelvin " << endl;  
      }
    }

    // Set detector time to something
    //
    det.Update(UTCDateTime(2005,1,10,12,0).GetTimeStamp());

    // Look up some data from a weather station for that time
    //
    const MolecularDB& wsdb = det.GetAtmosphere().GetMolecularDB(MolecularDB::eWeatherStation);
    
    cout << endl << endl << "some weather station data for  " << det.GetTime() << endl;
    // iterate over zones
    for (MolecularDB::ZoneIterator zIt = wsdb.ZonesBegin() ;
         zIt != wsdb.ZonesEnd() ; ++zIt){

      // iterate over molecular layers for this zone (well, there is only 1 for the weather station)
      for (MolecularZone::LayerIterator sIt = zIt->LayersBegin() ;
           sIt != zIt->LayersEnd() ; ++sIt){

        cout << endl;
        cout << "height      " << sIt->GetHeight()/m             << " m " << endl;
        cout << "depth       " << sIt->GetDepth()/(gram/cm/cm)   << " g/cm^2" << endl;
        cout << "pressure    " << sIt->GetPressure()/pascal      << " pascal " << endl;
        cout << "temperature " << sIt->GetTemperature()/kelvin   << " kelvin " << endl;  
        cout << "wind speed  " << sIt->GetWindSpeed()/(m/s)      << " m/s    " << endl;
      }
    }

  return eSuccess;
}

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