ReadAerosolLL.cc

Go to the documentation of this file.

#include <sstream>

#include "ReadAerosolLL.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>

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

using namespace ReadAerosolLLNS;

ReadAerosolLL::ReadAerosolLL(){}

ReadAerosolLL::~ReadAerosolLL(){}

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

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

    Detector& det = Detector::GetInstance();

    // Set detector time to something .
    //
    det.Update(UTCDateTime(2005, 5, 11, 9, 0).GetTimeStamp());


    // Retrieve the low-level aerosol DB interface
    //
    const AerosolDB& aerosol1 = det.GetAtmosphere().GetAerosolDB();

    // Iterate over all the zones found in the db
    //
    for (AerosolDB::ZoneIterator zIt = aerosol1.ZonesBegin() ;
         zIt != aerosol1.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 (AerosolZone::AttSliceIterator attIt = zIt->AttSlicesBegin() ;
           attIt != zIt->AttSlicesEnd() ; ++attIt) {

        cout << "     att length slice " << endl;
        cout << "       min, max height "
             << attIt->GetMinHeight()/m << " m, "
             << attIt->GetMaxHeight()/m << "m " << endl;

        cout << "       VAOD, max VAOD cor, min VAOD cor, max VAOD uncor, min VAOD uncor  "
             << attIt->GetVAOD() << " "
             << attIt->GetMaxVAODCor() << "  "
             << attIt->GetMinVAODCor() << " "
             << attIt->GetMaxVAODCor() << " "
             << attIt->GetMaxVAODUncor() << endl;

        cout << "       1/(att length), min cor, max cor, min uncor, max uncor "
             << attIt->GetAttAlpha() / (1./m) << " m^-1, "
             << attIt->GetMinAttAlphaCor() / (1./m) << " m^-1, "
             << attIt->GetMaxAttAlphaCor() / (1./m) << " m^-1 "
             << attIt->GetMinAttAlphaUncor() / (1./m) << " m^-1, "
             << attIt->GetMaxAttAlphaUncor() / (1./m) << " m^-1 " << endl;


        // etc.

      }
    }

  return eSuccess;
}

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