mute the SetAnalyzeSurface function in the FEM-DG solver to aovid calling GetNode() function, and create an error message

Author: Zan-AA

SU2_CFD/src/output/CFlowCompFEMOutput.cpp

@@ -257,7 +257,12 @@ void CFlowCompFEMOutput::LoadHistoryData(CConfig *config, CGeometry *geometry, C
 
   /*--- Set the analyse surface history values --- */
 
-  SetAnalyzeSurface(solver, geometry, config, false);
+  for (unsigned short iMarker = 0; iMarker < config->GetnMarker_All(); iMarker++) {
+    if (config->GetMarker_All_Analyze(iMarker) == YES) {
+      SU2_MPI::Error("SetAnalyzeSurface is not implemented for FEM-DG solver.", CURRENT_FUNCTION);
+    }
+  }
+  // SetAnalyzeSurface(solver, geometry, config, false);
 
   /*--- Set aeroydnamic coefficients --- */
 

SU2_CFD/src/output/CFlowOutput.cpp

@@ -155,10 +155,7 @@ void CFlowOutput::SetAnalyzeSurface(const CSolver* const*solver, const CGeometry
   const bool axisymmetric               = config->GetAxisymmetric();
   const unsigned short nMarker_Analyze  = config->GetnMarker_Analyze();
 
-  const auto flow_nodes = (config->GetKind_Solver() == MAIN_SOLVER::FEM_EULER               ||
-                           config->GetKind_Solver() == MAIN_SOLVER::FEM_NAVIER_STOKES       ||
-                           config->GetKind_Solver() == MAIN_SOLVER::FEM_RANS                ||
-                           config->GetKind_Solver() == MAIN_SOLVER::FEM_LES                 ) ? nullptr : solver[FLOW_SOL]->GetNodes();
+  const auto flow_nodes = solver[FLOW_SOL]->GetNodes();
   const CVariable* species_nodes = species ? solver[SPECIES_SOL]->GetNodes() : nullptr;
 
   vector<su2double> Surface_MassFlow          (nMarker,0.0);