source: XIOS/dev/dev_ym/XIOS_COUPLING/src/filter/model_to_client_source_filter.cpp @ 2010

#include "model_to_client_source_filter.hpp"
#include "grid.hpp"
#include "exception.hpp"
#include "calendar_util.hpp"
#include <limits> 

namespace xios
{
  CModelToClientSourceFilter::CModelToClientSourceFilter(CGarbageCollector& gc, CGrid* grid, bool hasMissingValue, double defaultValue)
    : COutputPin(gc)
    , grid_(grid)
    , hasMissingValue_(hasMissingValue), defaultValue_(defaultValue)
  {
    if (!grid)
      ERROR("CSourceFilter::CSourceFilter(CGrid* grid)",
            "Impossible to construct a source filter without providing a grid.");
  }

  template <int N>
  void CModelToClientSourceFilter::streamData(CDate date, const CArray<double, N>& data)
  {
    CDataPacketPtr packet(new CDataPacket);
    packet->date = date;
    packet->timestamp = date;
    packet->status = CDataPacket::NO_ERROR;
    
    auto connector = grid_->getModelToWorkflowConnector() ;

    if (connector->getSrcSize() != data.numElements())
      ERROR("void CModelToClientSourceFilter::streamData(CDate date, const CArray<double, N>& data)",
            << "[ Awaiting data of size = " <<connector->getSrcSize() << ", "
            << "Received data size = "      << data.numElements() << " ] "
            << "The data array does not have the right size! "
            << "grid = " << grid_->getId())
    const double nanValue = std::numeric_limits<double>::quiet_NaN();
    packet->data.resize(connector->getDstSize()) ;
    connector->transfer(data, packet->data, nanValue) ;

    if (hasMissingValue_)
    {
        const double nanValue = std::numeric_limits<double>::quiet_NaN();
        const size_t nbData = packet->data.numElements();
        for (size_t idx = 0; idx < nbData; ++idx)
          if (defaultValue_ == packet->data(idx))   packet->data(idx) = nanValue;
    }
    
    onOutputReady(packet);
  }

  template void CModelToClientSourceFilter::streamData<1>(CDate date, const CArray<double, 1>& data);
  template void CModelToClientSourceFilter::streamData<2>(CDate date, const CArray<double, 2>& data);
  template void CModelToClientSourceFilter::streamData<3>(CDate date, const CArray<double, 3>& data);
  template void CModelToClientSourceFilter::streamData<4>(CDate date, const CArray<double, 4>& data);
  template void CModelToClientSourceFilter::streamData<5>(CDate date, const CArray<double, 5>& data);
  template void CModelToClientSourceFilter::streamData<6>(CDate date, const CArray<double, 6>& data);
  template void CModelToClientSourceFilter::streamData<7>(CDate date, const CArray<double, 7>& data);

} // namespace xios