From cc353ee548a2189fc9289842cb0e0e0dabc4f7f5 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Vo=C3=9F=2C=20Arne?= Date: Fri, 24 Feb 2023 12:56:17 +0100 Subject: [PATCH] I think we don't need to set any source terms because they depend on the gust derivatives, which are zero in all implemented cases. This piece of code belongs to the split velocity method, which is currently not used as far as I understand. However, when trying to remove it, I broke the code elsewhere. --- SU2_CFD/src/numerics/flow/flow_sources.cpp | 9 +++++++-- 1 file changed, 7 insertions(+), 2 deletions(-) diff --git a/SU2_CFD/src/numerics/flow/flow_sources.cpp b/SU2_CFD/src/numerics/flow/flow_sources.cpp index 2fb164f3d1f..3d297157a51 100644 --- a/SU2_CFD/src/numerics/flow/flow_sources.cpp +++ b/SU2_CFD/src/numerics/flow/flow_sources.cpp @@ -784,8 +784,13 @@ CNumerics::ResidualType<> CSourceWindGust::ComputeResidual(const CConfig* config //residual[3] = smz*Volume; residual[3] = se*Volume; } else { - SU2_MPI::Error("You should only be in the gust source term in two dimensions", CURRENT_FUNCTION); - } + // SU2_MPI::Error("You should only be in the gust source term in two dimensions", CURRENT_FUNCTION); + residual[0] = 0.0; + residual[1] = 0.0; + residual[2] = 0.0; + residual[3] = 0.0; + residual[4] = 0.0; + } /*--- For now the source term Jacobian is just set to zero ---*/ //bool implicit = (config->GetKind_TimeIntScheme_Flow() == EULER_IMPLICIT);