bug fixes for 3D MP

Author: joshkellyjak

SU2_CFD/src/interfaces/cfd/CMixingPlaneInterface.cpp

@@ -64,7 +64,7 @@ void CMixingPlaneInterface::BroadcastData_MixingPlane(const CInterpolator& inter
 
   /*--- Loop over interface markers. ---*/
 
-  for (auto iMarkerInt = 0u; iMarkerInt < donor_config->GetMarker_n_ZoneInterface()/2; iMarkerInt++) {
+  for (auto iMarkerInt = 0u; iMarkerInt < donor_config->GetnMarker_MixingPlaneInterface()/2; iMarkerInt++) {
 
     /*--- Find the markers containing the interface ---*/
     short markDonor = donor_config->FindMixingPlaneInterfaceMarker(donor_geometry->GetnMarker());
@@ -126,7 +126,7 @@ void CMixingPlaneInterface::BroadcastData_MixingPlane(const CInterpolator& inter
 
       InitializeTarget_Variable(target_solution, markTarget, iTargetSpan, nSpanDonor);
 
-      if ((iTargetSpan == 0) || (iTargetSpan < nTargetSpan - 3)) {
+      if ((iTargetSpan == 0) || (iTargetSpan > nTargetSpan - 3)) {
         /*--- Transfer values at hub, shroud and 1D values ---*/
         unsigned long donorSpan;
         if (iTargetSpan == 0) donorSpan = 0;

SU2_CFD/src/solvers/CEulerSolver.cpp

@@ -453,12 +453,12 @@ void CEulerSolver::InitTurboContainers(CGeometry *geometry, CConfig **config_con
     }
   }
 
-  auto nMarkerMixingPlane = config->GetnMarker_MixingPlaneInterface();
+  auto nMarkerInterface = config->GetnMarker_ZoneInterface();
 
-  MixingState.resize(nMarkerMixingPlane);
-  MixingStateNodes.resize(nMarkerMixingPlane);
+  MixingState.resize(nMarkerInterface);
+  MixingStateNodes.resize(nMarkerInterface);
 
-  for (unsigned long iMarker = 0; iMarker < nMarkerMixingPlane; iMarker++) {
+  for (unsigned long iMarker = 0; iMarker < nMarkerInterface; iMarker++) {
     if (config->GetMarker_All_KindBC(iMarker) == GILES_BOUNDARY) {
       MixingState[iMarker].resize(nSpanWiseSections+1);
       MixingStateNodes[iMarker].resize(nSpanWiseSections+1);
@@ -9642,4 +9642,8 @@ void CEulerSolver::ComputeTurboBladePerformance(CGeometry* geometry, CConfig* co
     TotalPressureLoss = BladePerf->GetTotalPressureLoss();
     KineticEnergyLoss = BladePerf->GetKineticEnergyLoss();
   }
-}
+}
+
+// void CEulerSolver::RegisterSolutionExtra(bool input, CConfig* config) {
+//   nodes->RegisterGridVelocity();
+// }