void IdealMolalSoln::setStandardConcentrationModel(const std::string& model)
{
if (caseInsensitiveEquals(model, "unity")) {
m_formGC = 0;
} else if (caseInsensitiveEquals(model, "molar_volume")) {
m_formGC = 1;
} else if (caseInsensitiveEquals(model, "solvent_volume")) {
m_formGC = 2;
} else {
throw CanteraError("IdealSolnGasVPSS::setStandardConcentrationModel",
"Unknown standard concentration model '{}'", model);
}
}
void IdealMolalSoln::setCutoffModel(const std::string& model)
{
if (caseInsensitiveEquals(model, "none")) {
IMS_typeCutoff_ = 0;
} else if (caseInsensitiveEquals(model, "poly")) {
IMS_typeCutoff_ = 1;
} else if (caseInsensitiveEquals(model, "polyexp")) {
IMS_typeCutoff_ = 2;
} else {
throw CanteraError("IdealMolalSoln::setCutoffModel",
"Unknown cutoff model '{}'", model);
}
}
void MargulesVPSSTP::initThermoXML(XML_Node& phaseNode, const std::string& id_)
{
if ((int) id_.size() > 0) {
string idp = phaseNode.id();
if (idp != id_) {
throw CanteraError("MargulesVPSSTP::initThermoXML", "phasenode and Id are incompatible");
}
}
// Find the Thermo XML node
if (!phaseNode.hasChild("thermo")) {
throw CanteraError("MargulesVPSSTP::initThermoXML",
"no thermo XML node");
}
XML_Node& thermoNode = phaseNode.child("thermo");
// Make sure that the thermo model is Margules
if (!caseInsensitiveEquals(thermoNode["model"], "margules")) {
throw CanteraError("MargulesVPSSTP::initThermoXML",
"model name isn't Margules: " + thermoNode["model"]);
}
// Go get all of the coefficients and factors in the activityCoefficients
// XML block
if (thermoNode.hasChild("activityCoefficients")) {
XML_Node& acNode = thermoNode.child("activityCoefficients");
if (!caseInsensitiveEquals(acNode["model"], "margules")) {
throw CanteraError("MargulesVPSSTP::initThermoXML",
"Unknown activity coefficient model: " + acNode["model"]);
}
for (size_t i = 0; i < acNode.nChildren(); i++) {
XML_Node& xmlACChild = acNode.child(i);
// Process a binary salt field, or any of the other XML fields that
// make up the Pitzer Database. Entries will be ignored if any of
// the species in the entry isn't in the solution.
if (caseInsensitiveEquals(xmlACChild.name(), "binaryneutralspeciesparameters")) {
readXMLBinarySpecies(xmlACChild);
}
}
}
// Go down the chain
GibbsExcessVPSSTP::initThermoXML(phaseNode, id_);
}
void SWOptionFilter::setOptionValue(const char *ival) {
for (auto const & optValue : *optValues) {
if (caseInsensitiveEquals(optValue, ival)) {
optionValue = optValue;
option = (!strnicmp(ival, "On", 2)); // convenience for boolean filters
break;
}
}
}
