do not repeat loading the volume data, make flow dir unit in CEulerSolver

Author: pcarruscag

SU2_CFD/src/output/COutput.cpp

@@ -926,7 +926,7 @@ bool COutput::GetCauchyCorrectedTimeConvergence(const CConfig *config){
 bool COutput::SetResult_Files(CGeometry *geometry, CConfig *config, CSolver** solver_container,
                               unsigned long iter, bool force_writing) {
 
-  bool isFileWrite = false;
+  bool isFileWrite = false, dataIsLoaded = false;
   const auto nVolumeFiles = config->GetnVolumeOutputFiles();
   const auto* VolumeFiles = config->GetVolumeOutputFiles();
 
@@ -938,10 +938,13 @@ bool COutput::SetResult_Files(CGeometry *geometry, CConfig *config, CSolver** so
     /*--- Collect the volume data from the solvers.
      *  If time-domain is enabled, we also load the data although we don't output it,
      *  since we might want to do time-averaging. ---*/
-    if (WriteVolume_Output(config, iter, force_writing || cauchyTimeConverged, iFile) || config->GetTime_Domain())
-      LoadDataIntoSorter(config, geometry, solver_container);
+    const bool write_file = WriteVolume_Output(config, iter, force_writing || cauchyTimeConverged, iFile);
 
-    if (!(WriteVolume_Output(config, iter, force_writing || cauchyTimeConverged, iFile))) continue;
+    if ((write_file || config->GetTime_Domain()) && !dataIsLoaded) {
+      LoadDataIntoSorter(config, geometry, solver_container);
+      dataIsLoaded = true;
+    }
+    if (!write_file) continue;
 
     /*--- Partition and sort the data --- */
 
@@ -1039,10 +1042,10 @@ bool COutput::Convergence_Monitoring(CConfig *config, unsigned long Iteration) {
 
         oldFunc[iField_Conv] = newFunc[iField_Conv];
         newFunc[iField_Conv] = monitor;
-        /*--- Automatically modify the scaling factor of relative Cauchy convergence for 
-        * coefficients that are close to zero. Example: For the clean aircraft, the rolling 
-        * moment coefficient MOMENT_X is close to zero and thus will never reach a relative 
-        * cauchy convergence ->> dividing tiny numbers is not a good idea. Using absolute 
+        /*--- Automatically modify the scaling factor of relative Cauchy convergence for
+        * coefficients that are close to zero. Example: For the clean aircraft, the rolling
+        * moment coefficient MOMENT_X is close to zero and thus will never reach a relative
+        * cauchy convergence ->> dividing tiny numbers is not a good idea. Using absolute
         * cauchy convergence is more robust in this case. ---*/
         cauchyFunc = fabs(newFunc[iField_Conv] - oldFunc[iField_Conv]) / fmax(fabs(monitor), 0.1);
 

SU2_CFD/src/solvers/CEulerSolver.cpp

@@ -6456,7 +6456,7 @@ void CEulerSolver::BC_Inlet(CGeometry *geometry, CSolver **solver_container,
   unsigned short iDim;
   unsigned long iVertex, iPoint;
   su2double P_Total, T_Total, Velocity[MAXNDIM], Velocity2, H_Total, Temperature, Riemann,
-  Pressure, Density, Energy, *Flow_Dir, Mach2, SoundSpeed2, SoundSpeed_Total2, Vel_Mag,
+  Pressure, Density, Energy, Flow_Dir[MAXNDIM] = {0.0}, Mach2, SoundSpeed2, SoundSpeed_Total2, Vel_Mag,
   alpha, aa, bb, cc, dd, Area, UnitNormal[MAXNDIM], Normal[MAXNDIM] = {0.0};
   su2double *V_inlet, *V_domain;
 
@@ -6510,13 +6510,17 @@ void CEulerSolver::BC_Inlet(CGeometry *geometry, CSolver **solver_container,
 
         /*--- Total properties have been specified at the inlet. ---*/
 
-        case INLET_TYPE::TOTAL_CONDITIONS:
+        case INLET_TYPE::TOTAL_CONDITIONS: {
 
           /*--- Retrieve the specified total conditions for this inlet. ---*/
 
           P_Total  = Inlet_Ptotal[val_marker][iVertex];
           T_Total  = Inlet_Ttotal[val_marker][iVertex];
-          Flow_Dir = Inlet_FlowDir[val_marker][iVertex];
+          const su2double* dir = Inlet_FlowDir[val_marker][iVertex];
+          const su2double mag = GeometryToolbox::Norm(nDim, dir);
+          for (iDim = 0; iDim < nDim; iDim++) {
+            Flow_Dir[iDim] = dir[iDim] / mag
+          }
 
           /*--- Non-dim. the inputs if necessary. ---*/
 
@@ -6623,7 +6627,11 @@ void CEulerSolver::BC_Inlet(CGeometry *geometry, CSolver **solver_container,
 
           Density  = Inlet_Ttotal[val_marker][iVertex];
           Vel_Mag  = Inlet_Ptotal[val_marker][iVertex];
-          Flow_Dir = Inlet_FlowDir[val_marker][iVertex];
+          const su2double* dir = Inlet_FlowDir[val_marker][iVertex];
+          const su2double mag = GeometryToolbox::Norm(nDim, dir);
+          for (iDim = 0; iDim < nDim; iDim++) {
+            Flow_Dir[iDim] = dir[iDim] / mag
+          }
 
           /*--- Non-dim. the inputs if necessary. ---*/