Fixed hardcoding of velocity component in axisymmetric resdiual

Author: lkusch

SU2_CFD/include/numerics/turbulent/turb_sources.hpp

@@ -602,7 +602,7 @@ class CSourcePieceWise_TurbSST final : public CNumerics {
     AD::SetPreaccIn(V_i[idx.Velocity() + 1]);
 
     const su2double yinv = 1.0 / Coord_i[1];
-    const su2double rhov = Density_i * V_i[2];
+    const su2double rhov = Density_i * V_i[idx.Velocity() + 1];
     const su2double& k = ScalarVar_i[0];
     const su2double& w = ScalarVar_i[1];
 
@@ -612,7 +612,7 @@ class CSourcePieceWise_TurbSST final : public CNumerics {
 
     /*--- Production ---*/
     const su2double pk_axi = max(
-        0.0, 2.0 / 3.0 * rhov * k * ((2.0 * yinv * V_i[2] - PrimVar_Grad_i[2][1] - PrimVar_Grad_i[1][0]) / zeta - 1.0));
+        0.0, 2.0 / 3.0 * rhov * k * ((2.0 * yinv * V_i[idx.Velocity() + 1] - PrimVar_Grad_i[2][1] - PrimVar_Grad_i[1][0]) / zeta - 1.0));
     const su2double pw_axi = alfa_blended * zeta / k * pk_axi;
 
     /*--- Convection-Diffusion ---*/