ProfileCheckBasic.cc

Go to the documentation of this file.

/*
 * (C) Crown copyright 2020, Met Office
 *
 * 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/profile/ProfileCheckBasic.h"
#include "ufo/profile/VariableNames.h"
 
namespace ufo {
 
  static ProfileCheckMaker<ProfileCheckBasic> makerProfileCheckBasic_("Basic");
 
  ProfileCheckBasic::ProfileCheckBasic(const ConventionalProfileProcessingParameters &options)
    : ProfileCheckBase(options)
  {}
 
  void ProfileCheckBasic::runCheck(ProfileDataHandler &profileDataHandler)
  {
    oops::Log::debug() << " Basic checks" << std::endl;
 
    // Set basic check result to true
    result_ = true;
 
    // Skip this routine if specifically requested
    if (options_.BChecks_Skip.value())
      {
        oops::Log::debug() << "Skipping basic checks" << std::endl;
        return;
      }
 
    const int numProfileLevels = profileDataHandler.getNumProfileLevels();
    const std::vector <float> &pressures =
      profileDataHandler.get<float>(ufo::VariableNames::obs_air_pressure);
    // All QC flags are retrieved for the basic checks.
    // (Some might be empty; that is checked before they are used.)
    std::vector <int> &tFlags = profileDataHandler.get<int>
      (ufo::VariableNames::qcflags_air_temperature);
    std::vector <int> &zFlags = profileDataHandler.get<int>
      (ufo::VariableNames::qcflags_geopotential_height);
    std::vector <int> &uFlags = profileDataHandler.get<int>
      (ufo::VariableNames::qcflags_eastward_wind);
 
    // Warn and exit if pressures vector is empty
    if (pressures.empty()) {
      result_ = false;
      oops::Log::debug() << "Pressures vector is empty" << std::endl;
      return;
     }
 
    // Is the number of levels to check OK?
    bool numProfileLevelsOK = (numProfileLevels > 0);
 
    // Are any levels in the wrong order?
    bool pressOrderOK = true;
    for (int jlev = 0; jlev < numProfileLevels - 1; ++jlev) {
      pressOrderOK = pressOrderOK && (pressures[jlev] >= pressures[jlev + 1]);
      if (!pressOrderOK) break;
    }
 
    // Is the pressure at the first level > maximum value pressure?
    bool maxPressOK = (pressures.size() > 0 ?
                       pressures.front() <= options_.BChecks_maxValidP.value() :
                       false);
 
    // Is the pressure at the final level < minimum value pressure?
    bool minPressOK = (pressures.size() > 0 ?
                       pressures.back() > options_.BChecks_minValidP.value() :
                       false);
 
    oops::Log::debug() << " -> numProfileLevelsOK: " << numProfileLevelsOK << std::endl;
    oops::Log::debug() << " -> pressOrderOK: " << pressOrderOK << std::endl;
    oops::Log::debug() << " -> maxPressOK: " << maxPressOK << std::endl;
    oops::Log::debug() << " -> minPressOK: " << minPressOK << std::endl;
 
    result_ = numProfileLevelsOK && pressOrderOK && maxPressOK && minPressOK;
    oops::Log::debug() << " -> basicResult: " << result_ << std::endl;
 
    // If the basic checks are failed, set reject flags
    // This is not done in the OPS sonde consistency checks, but is done in Ops_SondeAverage.inc
    if (options_.flagBasicChecksFail.value() && !result_) {
      for (int jlev = 0; jlev < numProfileLevels; ++jlev) {
        if (!tFlags.empty()) tFlags[jlev] |= ufo::MetOfficeQCFlags::Elem::FinalRejectFlag;
        if (!zFlags.empty()) zFlags[jlev] |= ufo::MetOfficeQCFlags::Elem::FinalRejectFlag;
        if (!uFlags.empty()) uFlags[jlev] |= ufo::MetOfficeQCFlags::Elem::FinalRejectFlag;
      }
    }
  }
}  // namespace ufo