RejectObs.cc

Go to the documentation of this file.

/*
 * (C) Copyright 2019 UCAR
 *
 * This software is licensed under the terms of the Apache Licence Version 2.0
 * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
 */
 
#include "ufo/filters/actions/RejectObs.h"
 
#include "ioda/ObsDataVector.h"
#include "ufo/filters/ObsFilterData.h"
#include "ufo/filters/QCflags.h"
 
namespace ufo {
 
// -----------------------------------------------------------------------------
 
static FilterActionMaker<RejectObs> makerRejectObs_("reject");
 
// -----------------------------------------------------------------------------
 
RejectObs::RejectObs(const RejectObsParameters &parameters)
  : allvars_(), parameters_(parameters) {
}
 
// -----------------------------------------------------------------------------
 
void RejectObs::apply(const Variables & vars,
                      const std::vector<std::vector<bool>> & flagged,
                      const ObsFilterData &,
                      int filterQCflag,
                      ioda::ObsDataVector<int> & flags,
                      ioda::ObsDataVector<float> &) const {
  for (size_t jv = 0; jv < vars.nvars(); ++jv) {
    size_t iv = flags.varnames().find(vars.variable(jv).variable());
    for (size_t jobs = 0; jobs < flags.nlocs(); ++jobs) {
      if (flagged[jv][jobs] && flags[iv][jobs] == QCflags::pass)
        flags[iv][jobs] = filterQCflag;
    }
  }
}
 
// -----------------------------------------------------------------------------
 
}  // namespace ufo