Merge pull request #2477 from su2code/feature_CGNS_FixHOM

Merging in this one.

Author: vdweide

Common/src/geometry/meshreader/CCGNSMeshReaderBase.cpp

@@ -320,8 +320,38 @@ void CCGNSMeshReaderBase::ReadCGNSSectionMetadata() {
     /* Check for 1D elements in 2D problems, or for 2D elements in
      3D problems. If found, mark the section as a boundary section. */
 
-    if ((dimension == 2) && (elemType == BAR_2 || elemType == BAR_3)) isInterior[s] = false;
-    if ((dimension == 3) && (elemType == TRI_3 || elemType == QUAD_4)) isInterior[s] = false;
+    if (dimension == 2) {
+      switch (elemType) {
+        case BAR_2:
+        case BAR_3:
+        case BAR_4:
+        case BAR_5:
+          isInterior[s] = false;
+          break;
+        default:  // To avoid a compiler warning.
+          break;
+      }
+    } else {
+      switch (elemType) {
+        case TRI_3:
+        case TRI_6:
+        case TRI_9:
+        case TRI_10:
+        case TRI_12:
+        case TRI_15:
+        case QUAD_4:
+        case QUAD_8:
+        case QUAD_9:
+        case QUAD_12:
+        case QUAD_16:
+        case QUAD_P4_16:
+        case QUAD_25:
+          isInterior[s] = false;
+          break;
+        default:  // To avoid a compiler warning.
+          break;
+      }
+    }
 
     /*--- Increment the global element offset for each section
      based on whether or not this is a surface or volume section.

Common/src/geometry/meshreader/CCGNSMeshReaderFEM.cpp

@@ -39,6 +39,9 @@ CCGNSMeshReaderFEM::CCGNSMeshReaderFEM(const CConfig* val_config, unsigned short
   ReadCGNSDatabaseMetadata();
   ReadCGNSZoneMetadata();
 
+  /*--- Read the basic information about the sections. ---*/
+  ReadCGNSSectionMetadata();
+
   /*--- Read the volume connectivity and distribute it
         linearly over the MPI ranks. ---*/
   ReadCGNSVolumeElementConnectivity();