use the "real" current coordinates to compute surface rotations

Author: pcarruscag

SU2_CFD/src/solvers/CMeshSolver.cpp

@@ -969,7 +969,6 @@ void CMeshSolver::RestartOldGeometry(CGeometry *geometry, const CConfig *config)
 void CMeshSolver::Surface_Pitching(CGeometry *geometry, CConfig *config, unsigned long iter) {
 
   su2double deltaT, time_new, time_old, Lref;
-  const su2double* Coord = nullptr;
   su2double Center[3] = {0.0}, VarCoord[3] = {0.0}, Omega[3] = {0.0}, Ampl[3] = {0.0}, Phase[3] = {0.0};
   su2double VarCoordAbs[3] = {0.0};
   su2double rotCoord[3] = {0.0}, r[3] = {0.0};
@@ -1052,10 +1051,14 @@ void CMeshSolver::Surface_Pitching(CGeometry *geometry, CConfig *config, unsigne
 
       for (iVertex = 0; iVertex < geometry->nVertex[iMarker]; iVertex++) {
 
-        /*--- Index and coordinates of the current point ---*/
+        /*--- Index and coordinates of the current point accounting for other
+         * motions that may be applied, e.g. plunging. ---*/
 
         iPoint = geometry->vertex[iMarker][iVertex]->GetNode();
-        Coord  = geometry->nodes->GetCoord(iPoint);
+        su2double Coord[3] = {0.0};
+        for (iDim = 0; iDim < nDim; ++iDim) {
+          Coord[iDim] = nodes->GetMesh_Coord(iPoint, iDim) + nodes->GetBound_Disp(iPoint, iDim);
+        }
 
         /*--- Calculate non-dim. position from rotation center ---*/
 
@@ -1086,7 +1089,6 @@ void CMeshSolver::Surface_Pitching(CGeometry *geometry, CConfig *config, unsigne
 void CMeshSolver::Surface_Rotating(CGeometry *geometry, CConfig *config, unsigned long iter) {
 
   su2double deltaT, time_new, time_old, Lref;
-  const su2double* Coord = nullptr;
   su2double VarCoordAbs[3] = {0.0};
   su2double Center[3] = {0.0}, VarCoord[3] = {0.0}, Omega[3] = {0.0},
   rotCoord[3] = {0.0}, r[3] = {0.0}, Center_Aux[3] = {0.0};
@@ -1158,10 +1160,14 @@ void CMeshSolver::Surface_Rotating(CGeometry *geometry, CConfig *config, unsigne
 
       for (iVertex = 0; iVertex < geometry->nVertex[iMarker]; iVertex++) {
 
-        /*--- Index and coordinates of the current point ---*/
+        /*--- Index and coordinates of the current point accounting for other
+         * motions that may be applied, e.g. plunging. ---*/
 
         iPoint = geometry->vertex[iMarker][iVertex]->GetNode();
-        Coord  = geometry->nodes->GetCoord(iPoint);
+        su2double Coord[3] = {0.0};
+        for (iDim = 0; iDim < nDim; ++iDim) {
+          Coord[iDim] = nodes->GetMesh_Coord(iPoint, iDim) + nodes->GetBound_Disp(iPoint, iDim);
+        }
 
         /*--- Calculate non-dim. position from rotation center ---*/