research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

[1D] save/restore now includes boundary conditions of Inlet1D objects

Browse source 
This partially addresses Issue 90.
This commit is contained in:
Ray Speth 2012-12-18 00:13:02 +00:00
parent fdc13afb9f
commit 357ab0b668
1 changed files with 18 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

22
src/oneD/boundaries1D.cpp
View file

@ -250,15 +250,29 @@ save(XML_Node& o, const doublereal* const soln)
for (size_t k = 0; k < nComponents(); k++) {
ctml::addFloat(inlt, componentName(k), s[k], "", "",lowerBound(k), upperBound(k));
}
for (size_t k=0; k < m_nsp; k++) {
ctml::addFloat(inlt, "massFraction", m_yin[k], "",
m_flow->phase().speciesName(k));
}
}
void Inlet1D::
restore(const XML_Node& dom, doublereal* soln, int loglevel)
{
//map<string, double> x;
//getFloats(dom, x);
soln[0] = ctml::getFloat(dom, "mdot", "massflowrate"); // x["mdot"];
soln[1] = ctml::getFloat(dom, "temperature", "temperature"); // x["temperature"];
soln[0] = m_mdot = ctml::getFloat(dom, "mdot", "massflowrate");
soln[1] = m_temp = ctml::getFloat(dom, "temperature", "temperature");
m_yin.assign(m_nsp, 0.0);
for (size_t i = 0; i < dom.nChildren(); i++) {
const XML_Node& node = dom.child(i);
if (node.name() == "massFraction") {
size_t k = m_flow->phase().speciesIndex(node.attrib("type"));
if (k != npos) {
m_yin[k] = node.fp_value();
}
}
}
resize(2,1);
}