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

[1D] save/restore includes boundary conditions for Surf1D and OutletRes1D

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

22
src/oneD/boundaries1D.cpp
View file

@ -683,11 +683,29 @@ save(XML_Node& o, const doublereal* const soln)
outlt.addAttribute("points",1);
outlt.addAttribute("type","outletres");
outlt.addAttribute("components", double(nComponents()));
ctml::addFloat(outlt, "temperature", m_temp, "K");
for (size_t k=0; k < m_nsp; k++) {
ctml::addFloat(outlt, "massFraction", m_yres[k], "",
m_flow->phase().speciesName(k));
}
}
void OutletRes1D::
restore(const XML_Node& dom, doublereal* soln, int loglevel)
{
m_temp = ctml::getFloat(dom, "temperature");
m_yres.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_yres[k] = node.fp_value();
}
}
}
resize(1,1);
}
@ -777,9 +795,7 @@ save(XML_Node& o, const doublereal* const soln)
void Surf1D::
restore(const XML_Node& dom, doublereal* soln, int loglevel)
{
map<string, double> x;
ctml::getFloats(dom, x);
soln[0] = x["temperature"];
soln[0] = m_temp = ctml::getFloat(dom, "temperature", "temperature");
resize(1,1);
}