Skip to content

Commit 2ebf8f0

Browse files
author
TobiKattmann
committed
Fix Mesh_Quality volume output for incompressible flow introduced in #1104.
1 parent 109f37b commit 2ebf8f0

1 file changed

Lines changed: 8 additions & 7 deletions

File tree

SU2_CFD/src/output/CFlowIncOutput.cpp

Lines changed: 8 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -629,6 +629,14 @@ void CFlowIncOutput::LoadVolumeData(CConfig *config, CGeometry *geometry, CSolve
629629
}
630630
SetVolumeOutputValue("Q_CRITERION", iPoint, GetQ_Criterion(&(Node_Flow->GetGradient_Primitive(iPoint)[1])));
631631
}
632+
633+
// Mesh quality metrics
634+
if (config->GetWrt_MeshQuality()) {
635+
SetVolumeOutputValue("ORTHOGONALITY", iPoint, geometry->Orthogonality[iPoint]);
636+
SetVolumeOutputValue("ASPECT_RATIO", iPoint, geometry->Aspect_Ratio[iPoint]);
637+
SetVolumeOutputValue("VOLUME_RATIO", iPoint, geometry->Volume_Ratio[iPoint]);
638+
}
639+
632640
}
633641

634642
void CFlowIncOutput::LoadSurfaceData(CConfig *config, CGeometry *geometry, CSolver **solver, unsigned long iPoint, unsigned short iMarker, unsigned long iVertex){
@@ -648,13 +656,6 @@ void CFlowIncOutput::LoadSurfaceData(CConfig *config, CGeometry *geometry, CSolv
648656
SetVolumeOutputValue("Y_PLUS", iPoint, solver[FLOW_SOL]->GetYPlus(iMarker, iVertex));
649657
}
650658

651-
// Mesh quality metrics
652-
if (config->GetWrt_MeshQuality()) {
653-
SetVolumeOutputValue("ORTHOGONALITY", iPoint, geometry->Orthogonality[iPoint]);
654-
SetVolumeOutputValue("ASPECT_RATIO", iPoint, geometry->Aspect_Ratio[iPoint]);
655-
SetVolumeOutputValue("VOLUME_RATIO", iPoint, geometry->Volume_Ratio[iPoint]);
656-
}
657-
658659
}
659660

660661
bool CFlowIncOutput::SetInit_Residuals(CConfig *config){

0 commit comments

Comments
 (0)