Replaced variable length array declarations

Author: EvertBunschoten

SU2_CFD/src/solvers/CEulerSolver.cpp

@@ -8586,11 +8586,14 @@ void CEulerSolver::TurboAverageProcess(CSolver **solver, CGeometry *geometry, CC
   const auto nSpanWiseSections = config->GetnSpanWiseSections();
 
   for (auto iSpan= 0; iSpan < nSpanWiseSections + 1; iSpan++){
-    su2double TotalDensity{0}, TotalPressure{0}, TotalNu{0}, TotalOmega{0}, TotalKine{0}, TotalVelocity[nDim], 
-              TotalAreaDensity{0}, TotalAreaPressure{0}, TotalAreaNu{0}, TotalAreaOmega{0}, TotalAreaKine{0}, TotalAreaVelocity[nDim],
-              TotalMassDensity{0}, TotalMassPressure{0}, TotalMassNu{0}, TotalMassOmega{0}, TotalMassKine{0}, TotalMassVelocity[nDim],
-              TotalFluxes[nVar];
-
+    su2double TotalDensity{0}, TotalPressure{0}, TotalNu{0}, TotalOmega{0}, TotalKine{0}, 
+              TotalAreaDensity{0}, TotalAreaPressure{0}, TotalAreaNu{0}, TotalAreaOmega{0}, TotalAreaKine{0},
+              TotalMassDensity{0}, TotalMassPressure{0}, TotalMassNu{0}, TotalMassOmega{0}, TotalMassKine{0};
+
+    su2double *TotalVelocity = new su2double[nDim],
+              *TotalAreaVelocity = new su2double[nDim],
+              *TotalMassVelocity = new su2double[nDim],
+              *TotalFluxes = new su2double[nVar];
     /*--- Forces initialization for contenitors ---*/
     for (auto iVar=0u;iVar<nVar;iVar++)
       TotalFluxes[iVar]= 0.0;
@@ -8610,8 +8613,11 @@ void CEulerSolver::TurboAverageProcess(CSolver **solver, CGeometry *geometry, CC
       auto Density  = nodes->GetDensity(iPoint);
       auto Enthalpy = nodes->GetEnthalpy(iPoint);
 
-      /*--- Normal vector for this vertex (negate for outward convention) ---*/
-      su2double Velocity[nDim], UnitNormal[nDim], TurboNormal[nDim], TurboVelocity[nDim];
+      su2double *Velocity = new su2double[nDim],
+                *UnitNormal = new su2double[nDim],
+                *TurboNormal = new su2double[nDim],
+                *TurboVelocity = new su2double[nDim];
+
       geometry->turbovertex[iMarker][iSpan][iVertex]->GetNormal(UnitNormal);
       geometry->turbovertex[iMarker][iSpan][iVertex]->GetTurboNormal(TurboNormal);
       auto Area = geometry->turbovertex[iMarker][iSpan][iVertex]->GetArea();
@@ -8666,6 +8672,11 @@ void CEulerSolver::TurboAverageProcess(CSolver **solver, CGeometry *geometry, CC
         TotalMassOmega   += Area*(Density*TurboVelocity[0] )*Omega;
         TotalMassNu      += Area*(Density*TurboVelocity[0] )*Nu;
       }
+
+      delete [] Velocity;
+      delete [] UnitNormal;
+      delete [] TurboNormal;
+      delete [] TurboVelocity;
     };
 
     for (auto iMarker = 0u; iMarker < config->GetnMarker_All(); iMarker++){
@@ -8753,9 +8764,9 @@ void CEulerSolver::TurboAverageProcess(CSolver **solver, CGeometry *geometry, CC
             /*--- Compute the averaged value for the boundary of interest for the span of interest ---*/
 
             bool belowMachLimit = (abs(MachTest)< config->GetAverageMachLimit());
-            su2double avgDensity{0}, avgPressure{0}, avgKine{0}, avgOmega{0}, avgNu{0}, avgVelocity[nDim],
-                      avgMixTurboVelocity[nDim];
-
+            su2double avgDensity{0}, avgPressure{0}, avgKine{0}, avgOmega{0}, avgNu{0};
+            su2double *avgVelocity = new su2double[nDim],
+                      *avgMixTurboVelocity = new su2double[nDim];
             for (auto iVar = 0u; iVar<nVar; iVar++){
               AverageFlux[iMarker][iSpan][iVar]   = TotalFluxes[iVar]/TotalArea;
               SpanTotalFlux[iMarker][iSpan][iVar] = TotalFluxes[iVar];
@@ -8908,10 +8919,16 @@ void CEulerSolver::TurboAverageProcess(CSolver **solver, CGeometry *geometry, CC
             } else {
               ComputeTurboVelocity(avgVelocity, AverageTurboNormal , TurboVel, marker_flag, config->GetKind_TurboMachinery(iZone));
             }
+            delete [] avgVelocity;
+            delete [] avgMixTurboVelocity;
           }
         }
       } // iMarkerTP
     } // iMarker
+    delete [] TotalVelocity;
+    delete [] TotalMassVelocity;
+    delete [] TotalAreaVelocity;
+    delete [] TotalFluxes;
   } // iSpan
 
   /*--- Compute Outlet Static Pressure if Radial equilibrium is imposed ---*/