OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fixed: allow 4D analytical solution output to VTF

Browse Source 
this is necessary for 3D CFD (u,v,w,p).

do this by using the scalar solution field for the pressure
This commit is contained in:
Arne Morten Kvarving
2015-12-14 13:20:11 +01:00
committed by Knut Morten Okstad
parent 3357188686
commit cfcf0fb06c
1 changed files with 6 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
src/SIM/SIMoutput.C
View File

@@ -478,6 +478,12 @@ int SIMoutput::writeGlvS1 (const Vector& psol, int iStep, int& nBlock,
const VecFunc& pSol = *mySol->getVectorSol();
for (j = 1; cit != grid->end_XYZ() && haveXsol; j++, ++cit)
field.fillColumn(j,pSol(Vec4(*cit,time)).ptr());
if (mySol->getScalarSol()) {
cit = grid->begin_XYZ();
const RealFunc& sSol = *mySol->getScalarSol();
for (j = 1; cit != grid->end_XYZ() && haveXsol; j++, ++cit)
field(field.rows(), j) = sSol(Vec4(*cit,time));
}
}
for (j = 1; j <= field.rows() && k < pMAX && haveXsol; j++)