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

[1D] save/restore now includes the temperature fixed point

Browse source 
This partially addresses Issue 90.
This commit is contained in:
Ray Speth 2012-12-18 00:12:41 +00:00
parent 833ca13dc1
commit fdc13afb9f
1 changed files with 11 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

11
src/oneD/StFlow.cpp
View file

@ -109,6 +109,9 @@ StFlow::StFlow(IdealGasPhase* ph, size_t nsp, size_t points) :
m_points = points;
m_thermo = ph;
m_zfixed = Undef;
m_tfixed = Undef;
if (ph == 0) {
return; // used to create a dummy object
}
@ -1064,6 +1067,8 @@ void StFlow::restore(const XML_Node& dom, doublereal* soln, int loglevel)
pp = getFloat(dom, "pressure", "pressure");
setPressure(pp);
getOptionalFloat(dom, "t_fixed", m_tfixed);
getOptionalFloat(dom, "z_fixed", m_zfixed);
vector<XML_Node*> d;
dom.child("grid_data").getChildren("floatArray",d);
@ -1217,6 +1222,12 @@ void StFlow::save(XML_Node& o, const doublereal* const sol)
}
XML_Node& gv = flow.addChild("grid_data");
addFloat(flow, "pressure", m_press, "Pa", "pressure");
if (m_zfixed != Undef) {
addFloat(flow, "z_fixed", m_zfixed, "m");
addFloat(flow, "t_fixed", m_tfixed, "K");
}
addFloatArray(gv,"z",m_z.size(),DATA_PTR(m_z),
"m","length");
vector_fp x(soln.nColumns());