Commit 150bfee8 authored by Chris Cantwell's avatar Chris Cantwell Committed by Spencer Sherwin
Browse files

Output HalfMode solutions as a 2D plane.

parent 5517c0f2
......@@ -362,6 +362,13 @@ namespace Nektar
int expansion,
int istrip)
{
// If there is only one plane (e.g. HalfMode), we write a 2D plane.
if (m_planes.num_elements() == 1)
{
m_planes[0]->WriteVtkPieceHeader(outfile, expansion);
return;
}
int i,j,k;
int nq0 = (*m_exp)[expansion]->GetNumPoints(0);
int nq1 = (*m_exp)[expansion]->GetNumPoints(1);
......
......@@ -853,6 +853,13 @@ namespace Nektar
void ExpListHomogeneous1D::v_WriteVtkPieceData(std::ostream &outfile, int expansion,
std::string var)
{
// If there is only one plane (e.g. HalfMode), we write a 2D plane.
if (m_planes.num_elements() == 1)
{
m_planes[0]->WriteVtkPieceData(outfile, expansion, var);
return;
}
int i;
int nq = (*m_exp)[expansion]->GetTotPoints();
int npoints_per_plane = m_planes[0]->GetTotPoints();
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment