Method toWthString

src/utilities/filetypes/EpwFile.cpp:1589–1645

}

boost::optional<std::string> EpwDataPoint::toWthString() const {
  std::string date = fmt::format("{}/{}", m_month, m_day);
  std::string string = date;
  std::string hms = fmt::format("{:02d}:{:02d}:00", m_hour, m_minute);
  string += '\t' + hms;
  boost::optional<double> value = dryBulbTemperature();
  if (!value) {
    LOG_FREE(Error, "openstudio.EpwFile", "Missing dry bulb temperature on " << date << " at " << hms)
    return boost::none;
  }
  //boost::optional<double> optdrybulb = openstudio::convert(value.get(), "C", "K");
  //OS_ASSERT(optdrybulb);
  double drybulb = value.get();
  string += '\t' + std::to_string(drybulb + 273.15);
  value = atmosphericStationPressure();
  if (!value) {
    LOG_FREE(Error, "openstudio.EpwFile", "Missing atmospheric station pressure on " << date << " at " << hms);
    return boost::none;
  }
  double p = value.get();
  string += '\t' + m_atmosphericStationPressure;
  if (!windSpeed()) {
    LOG_FREE(Error, "openstudio.EpwFile", "Missing wind speed on " << date << " at " << hms);
    return boost::none;
  }
  string += '\t' + m_windSpeed;
  if (!windDirection()) {
    LOG_FREE(Error, "openstudio.EpwFile", "Missing wind direction on " << date << " at " << hms);
    return boost::none;
  }
  string += '\t' + m_windDirection;
  double pw;
  value = relativeHumidity();
  if (!value) {  // Don't have relative humidity - this has not been tested
    value = dewPointTemperature();
    if (!value) {
      LOG_FREE(Error, "openstudio.EpwFile", "Cannot compute humidity ratio on " << date << " at " << hms);
      return boost::none;
    }
    //boost::optional<double> optdewpoint = openstudio::convert(value.get(), "C", "K");
    //OS_ASSERT(optdewpoint);
    double dewpoint = value.get();
    pw = openstudio::psat(dewpoint);
  } else {  // Have relative humidity
    double pws = openstudio::psat(drybulb);
    pw = 0.01 * value.get() * pws;
  }
  double W = 0.621945 * pw / (p - pw);
  string += "\t" + std::to_string(W * 1000);  // need g/kg
  // Pass on solar flux quantities
  string += "\t0\t0";
  // Pass on Tsky
  string += "\t0";
  // Pass on snow and rain
  string += "\t0\t0";
  return boost::optional<std::string>(string);
}