diff --git a/SU2_CFD/include/numerics/NEMO/convection/ausm_slau.hpp b/SU2_CFD/include/numerics/NEMO/convection/ausm_slau.hpp
index 0d01c9a37ad..c6066d740f8 100644
--- a/SU2_CFD/include/numerics/NEMO/convection/ausm_slau.hpp
+++ b/SU2_CFD/include/numerics/NEMO/convection/ausm_slau.hpp
@@ -120,8 +120,8 @@ class CUpwAUSM_NEMO final : public CUpwAUSM_SLAU_Base_NEMO {
    * \brief Constructor of the class.
    * \param[in] val_nDim - Number of dimensions of the problem.
    * \param[in] val_nVar - Number of variables of the problem.
-   * \param[in] val_nPrimVar - Number of primitive variables of the problem
-   * \param[in] val_nPrimVarGrad - Number of grad primitive variables of the problem
+   * \param[in] val_nPrimVar - Number of primitive variables of the problem.
+   * \param[in] val_nPrimVarGrad - Number of grad primitive variables of the problem.
    * \param[in] config - Definition of the particular problem.
    */
   CUpwAUSM_NEMO(unsigned short val_nDim, unsigned short val_nVar, unsigned short val_nPrimVar,
@@ -131,7 +131,9 @@ class CUpwAUSM_NEMO final : public CUpwAUSM_SLAU_Base_NEMO {
 /*!
  * \class CUpwAUSMPLUSM_NEMO
  * \brief Class for solving an approximate Riemann AUSM+ M, Two-Temperature Model.
- * https://doi.org/10.1016/j.apm.2019.09.005 \ingroup ConvDiscr \author F. Morgado
+ * https://doi.org/10.1016/j.apm.2019.09.005
+ * \ingroup ConvDiscr
+ * \author F. Morgado
  */
 class CUpwAUSMPLUSM_NEMO final : public CUpwAUSM_SLAU_Base_NEMO {
  private:
@@ -160,10 +162,46 @@ class CUpwAUSMPLUSM_NEMO final : public CUpwAUSM_SLAU_Base_NEMO {
                      unsigned short val_nPrimVarGrad, const CConfig* config);
 };
 
+/*!
+ * \class CUpwAUSMPLUSUP_NEMO
+ * \brief Class for solving an approximate Riemann AUSM+-Up, Two-Temperature Model.
+ * \ingroup ConvDiscr
+ * \author Amit Sachdeva, P. Gomes, W. Maier
+ */
+class CUpwAUSMPLUSUP_NEMO final : public CUpwAUSM_SLAU_Base_NEMO {
+ private:
+  su2double Kp, Ku, sigma;
+
+  /*!
+   * \brief Compute the interface Mach number, soundspeeds and pressure for AUSM+-Up scheme.
+   * \param[in] config - Definition of the particular problem.
+   * \param[out] pressure - The pressure at the control volume face.
+   * \param[out] interface_mach - The interface Mach number M_(1/2).
+   * \param[out] interface_soundspeed - The interface soundspeed (vector for i and j faces if necessary).
+   */
+  virtual void ComputeInterfaceQuantities(const CConfig* config, su2double* pressure, su2double& interface_mach,
+                                          su2double* interface_soundspeed) override;
+
+ public:
+  /*!
+   * \brief Constructor of the class.
+   * \param[in] val_nDim - Number of dimensions of the problem.
+   * \param[in] val_nVar - Number of variables of the problem.
+   * \param[in] val_nPrimVar - Number of primitive variables of the problem.
+   * \param[in] val_nPrimVarGrad - Number of grad primitive variables of the problem.
+   * \param[in] config - Definition of the particular problem.
+   */
+  CUpwAUSMPLUSUP_NEMO(unsigned short val_nDim, unsigned short val_nVar, unsigned short val_nPrimVar,
+                       unsigned short val_nPrimVarGrad, const CConfig* config);
+};
+
+
 /*!
  * \class CUpwAUSMPLUSUP2_NEMO
  * \brief Class for solving an approximate Riemann AUSM+-up2, Two-Temperature Model.
- * https://doi.org/10.1016/j.jcp.2013.02.046 \ingroup ConvDiscr \author W. Maier, A. Sachedeva, C. Garbacz
+ * https://doi.org/10.1016/j.jcp.2013.02.046
+ * \ingroup ConvDiscr
+ * \author W. Maier, A. Sachedeva, C. Garbacz
  */
 class CUpwAUSMPLUSUP2_NEMO final : public CUpwAUSM_SLAU_Base_NEMO {
  private:
@@ -184,8 +222,8 @@ class CUpwAUSMPLUSUP2_NEMO final : public CUpwAUSM_SLAU_Base_NEMO {
    * \brief Constructor of the class.
    * \param[in] val_nDim - Number of dimensions of the problem.
    * \param[in] val_nVar - Number of variables of the problem.
-   * \param[in] val_nPrimVar - Number of primitive variables of the problem
-   * \param[in] val_nPrimVarGrad - Number of grad primitive variables of the problem
+   * \param[in] val_nPrimVar - Number of primitive variables of the problem.
+   * \param[in] val_nPrimVarGrad - Number of grad primitive variables of the problem.
    * \param[in] config - Definition of the particular problem.
    */
   CUpwAUSMPLUSUP2_NEMO(unsigned short val_nDim, unsigned short val_nVar, unsigned short val_nPrimVar,
@@ -193,17 +231,18 @@ class CUpwAUSMPLUSUP2_NEMO final : public CUpwAUSM_SLAU_Base_NEMO {
 };
 
 /*!
- * \class CUpwAUSMPWplus_NEMO
- * \brief Class for solving an approximate Riemann AUSM.
+ * \class CUpwSLAU_NEMO
+ * \brief Class for solving the Low-Dissipation AUSM for the NEMO solver.
  * \ingroup ConvDiscr
- * \author F. Palacios, W.Maier, C. Garbacz
+ * \author E. Molina, P. Gomes, W. Maier
  */
-class CUpwAUSMPWplus_NEMO : public CUpwAUSM_SLAU_Base_NEMO {
- private:
-  su2double alpha;
+class CUpwSLAU_NEMO : public CUpwAUSM_SLAU_Base_NEMO {
+protected:
+  bool slau_low_diss;
+  bool slau2;
 
   /*!
-   * \brief Compute the interface Mach number, soundspeeds and pressure for AUSMpw+ scheme.
+   * \brief Compute the interface Mach number, soundspeeds and pressure for SLAU schemes.
    * \param[in] config - Definition of the particular problem.
    * \param[out] pressure - The pressure at the control volume face.
    * \param[out] interface_mach - The interface Mach number M_(1/2).
@@ -212,15 +251,38 @@ class CUpwAUSMPWplus_NEMO : public CUpwAUSM_SLAU_Base_NEMO {
   virtual void ComputeInterfaceQuantities(const CConfig* config, su2double* pressure, su2double& interface_mach,
                                           su2double* interface_soundspeed) override;
 
- public:
+public:
+
+  /*!
+   * \brief Constructor of the class.
+   * \param[in] val_nDim - Number of dimensions of the problem.
+   * \param[in] val_nVar - Number of variables of the problem.
+   * \param[in] val_nPrimVar - Number of primitive variables of the problem.
+   * \param[in] val_nPrimVarGrad - Number of grad primitive variables of the problem.
+   * \param[in] config - Definition of the particular problem.
+   * \param[in] val_low_dissipation - Bool to determine if low dissipation used.
+   */
+  CUpwSLAU_NEMO(unsigned short val_nDim, unsigned short val_nVar, unsigned short val_nPrimVar,
+                unsigned short val_nPrimVarGrad, const CConfig* config, bool val_low_dissipation);
+};
+
+/*!
+ * \class CUpwSLAU2_NEMO
+ * \brief Class for solving the Simple Low-Dissipation AUSM 2 for the NEMO solver.
+ * \ingroup ConvDiscr
+ * \author E. Molina, P. Gomes, W. Maier
+ */
+class CUpwSLAU2_NEMO final : public CUpwSLAU_NEMO {
+public:
   /*!
    * \brief Constructor of the class.
    * \param[in] val_nDim - Number of dimensions of the problem.
    * \param[in] val_nVar - Number of variables of the problem.
-   * \param[in] val_nPrimVar - Number of primitive variables of the problem
-   * \param[in] val_nPrimVarGrad - Number of grad primitive variables of the problem
+   * \param[in] val_nPrimVar - Number of primitive variables of the problem.
+   * \param[in] val_nPrimVarGrad - Number of grad primitive variables of the problem.
    * \param[in] config - Definition of the particular problem.
+   * \param[in] val_low_dissipation - Bool to determine if low dissipation used.
    */
-  CUpwAUSMPWplus_NEMO(unsigned short val_nDim, unsigned short val_nVar, unsigned short val_nPrimVar,
-                      unsigned short val_nPrimVarGrad, const CConfig* config);
+  CUpwSLAU2_NEMO(unsigned short val_nDim, unsigned short val_nVar, unsigned short val_nPrimVar,
+                 unsigned short val_nPrimVarGrad, const CConfig* config, bool val_low_dissipation);
 };
diff --git a/SU2_CFD/src/drivers/CDriver.cpp b/SU2_CFD/src/drivers/CDriver.cpp
index d8ba985fafa..fcc44a55e27 100644
--- a/SU2_CFD/src/drivers/CDriver.cpp
+++ b/SU2_CFD/src/drivers/CDriver.cpp
@@ -1942,6 +1942,13 @@ void CDriver::InitializeNumerics(CConfig *config, CGeometry **geometry, CSolver
               }
               break;
 
+            case UPWIND::AUSMPLUSUP:
+              for (iMGlevel = 0; iMGlevel <= config->GetnMGLevels(); iMGlevel++) {
+                numerics[iMGlevel][FLOW_SOL][conv_term] = new CUpwAUSMPLUSUP_NEMO(nDim, nVar_NEMO, nPrimVar_NEMO, nPrimVarGrad_NEMO, config);
+                numerics[iMGlevel][FLOW_SOL][conv_bound_term] = new CUpwAUSMPLUSUP_NEMO(nDim, nVar_NEMO, nPrimVar_NEMO, nPrimVarGrad_NEMO, config);
+              }
+              break;
+
             case UPWIND::AUSMPLUSUP2:
               for (iMGlevel = 0; iMGlevel <= config->GetnMGLevels(); iMGlevel++) {
                 numerics[iMGlevel][FLOW_SOL][conv_term] = new CUpwAUSMPLUSUP2_NEMO(nDim, nVar_NEMO, nPrimVar_NEMO, nPrimVarGrad_NEMO, config);
@@ -1956,6 +1963,20 @@ void CDriver::InitializeNumerics(CConfig *config, CGeometry **geometry, CSolver
               }
               break;
 
+            case UPWIND::SLAU:
+              for (iMGlevel = 0; iMGlevel <= config->GetnMGLevels(); iMGlevel++) {
+                numerics[iMGlevel][FLOW_SOL][conv_term] = new CUpwSLAU_NEMO(nDim, nVar_NEMO, nPrimVar_NEMO, nPrimVarGrad_NEMO, config, roe_low_dissipation);
+                numerics[iMGlevel][FLOW_SOL][conv_bound_term] = new CUpwSLAU_NEMO(nDim, nVar_NEMO, nPrimVar_NEMO, nPrimVarGrad_NEMO, config, false);
+              }
+              break;
+
+            case UPWIND::SLAU2:
+              for (iMGlevel = 0; iMGlevel <= config->GetnMGLevels(); iMGlevel++) {
+                numerics[iMGlevel][FLOW_SOL][conv_term] = new CUpwSLAU2_NEMO(nDim, nVar_NEMO, nPrimVar_NEMO, nPrimVarGrad_NEMO, config, roe_low_dissipation);
+                numerics[iMGlevel][FLOW_SOL][conv_bound_term] = new CUpwSLAU2_NEMO(nDim, nVar_NEMO, nPrimVar_NEMO, nPrimVarGrad_NEMO, config, false);
+              }
+              break;
+
             case UPWIND::MSW:
               for (iMGlevel = 0; iMGlevel <= config->GetnMGLevels(); iMGlevel++) {
                 numerics[iMGlevel][FLOW_SOL][conv_term] = new CUpwMSW_NEMO(nDim, nVar_NEMO, nPrimVar_NEMO, nPrimVarGrad_NEMO, config);
diff --git a/SU2_CFD/src/numerics/NEMO/convection/ausm_slau.cpp b/SU2_CFD/src/numerics/NEMO/convection/ausm_slau.cpp
index 865cdc268d9..beb903ab349 100644
--- a/SU2_CFD/src/numerics/NEMO/convection/ausm_slau.cpp
+++ b/SU2_CFD/src/numerics/NEMO/convection/ausm_slau.cpp
@@ -482,6 +482,90 @@ void CUpwAUSM_NEMO::ComputeInterfaceQuantities(const CConfig* config, su2double*
   for (auto iDim = 0ul; iDim < nDim; iDim++) pressure[iDim] = (P_LP + P_RM) * UnitNormal[iDim];
 }
 
+CUpwAUSMPLUSUP_NEMO::CUpwAUSMPLUSUP_NEMO(unsigned short val_nDim, unsigned short val_nVar,
+                                           unsigned short val_nPrimVar, unsigned short val_nPrimVarGrad,
+                                           const CConfig* config)
+    : CUpwAUSM_SLAU_Base_NEMO(val_nDim, val_nVar, val_nPrimVar, val_nPrimVarGrad, config) {
+
+  Minf = config->GetMach();
+  Kp = 0.25;
+  Ku = 0.75;
+  sigma = 1.0;
+
+  if (Minf < EPS)
+    SU2_MPI::Error("AUSM+Up requires a reference Mach number (\"MACH_NUMBER\") greater than 0.", CURRENT_FUNCTION);
+}
+
+void CUpwAUSMPLUSUP_NEMO::ComputeInterfaceQuantities(const CConfig* config, su2double* pressure,
+                                                    su2double& interface_mach, su2double* interface_soundspeed) {
+
+  /*--- Compute interface speed of sound (aF) ---*/
+
+  const su2double astarL = sqrt(2.0*(Gamma_i-1.0)/(Gamma_i+1.0)*Enthalpy_i);
+  const su2double astarR = sqrt(2.0*(Gamma_j-1.0)/(Gamma_j+1.0)*Enthalpy_j);
+
+  const su2double ahatL = astarL*astarL/max(astarL, ProjVelocity_i);
+  const su2double ahatR = astarR*astarR/max(astarR,-ProjVelocity_j);
+
+  const su2double A_F = min(ahatL,ahatR);
+  interface_soundspeed[0] = interface_soundspeed[1] = A_F;
+
+  /*--- Left and right pressures and Mach numbers ---*/
+
+  su2double mLP, betaLP, mRM, betaRM;
+
+  const su2double mL = ProjVelocity_i/A_F;
+  const su2double mR = ProjVelocity_j/A_F;
+
+  const su2double MFsq = 0.5*(mL*mL+mR*mR);
+  const su2double param1 = max(MFsq, Minf*Minf);
+  const su2double Mrefsq = min(1.0, param1);
+
+  const su2double fa = 2.0*sqrt(Mrefsq)-Mrefsq;
+
+  const su2double alpha = 3.0/16.0*(-4.0+5.0*fa*fa);
+  const su2double beta = 1.0/8.0;
+
+  if (fabs(mL) <= 1.0) {
+    const su2double p1 = 0.25*(mL+1.0)*(mL+1.0);
+    const su2double p2 = (mL*mL-1.0)*(mL*mL-1.0);
+
+    mLP = p1 + beta*p2;
+    betaLP = p1*(2.0-mL) + alpha*mL*p2;
+  }
+  else {
+    mLP = 0.5*(mL+fabs(mL));
+    betaLP = mLP/mL;
+  }
+
+  if (fabs(mR) <= 1.0) {
+    const su2double p1 = 0.25*(mR-1.0)*(mR-1.0);
+    const su2double p2 = (mR*mR-1.0)*(mR*mR-1.0);
+
+    mRM = -p1 - beta*p2;
+    betaRM = p1*(2.0+mR) - alpha*mR*p2;
+  }
+  else {
+    mRM = 0.5*(mR-fabs(mR));
+    betaRM = mRM/mR;
+  }
+
+  /*--- Pressure and velocity diffusion terms ---*/
+
+  const su2double rhoF = 0.5*(Density_i+Density_j);
+  const su2double Mp = -(Kp/fa)*max((1.0-sigma*MFsq),0.0)*(Pressure_j-Pressure_i)/(rhoF*A_F*A_F);
+
+  const su2double Pu = -Ku*fa*betaLP*betaRM*2.0*rhoF*A_F*(ProjVelocity_j-ProjVelocity_i);
+
+  /*--- Finally the fluxes ---*/
+
+  interface_mach = mLP + mRM + Mp;
+  su2double mdot = aF * (max(mF,0.0)*Density_i + min(mF,0.0)*Density_j);
+
+  for (auto iDim = 0ul; iDim < nDim; iDim++)
+    pressure[iDim] = betaLP*Pressure_i + betaRM*Pressure_j + Pu;
+}
+
 CUpwAUSMPLUSUP2_NEMO::CUpwAUSMPLUSUP2_NEMO(unsigned short val_nDim, unsigned short val_nVar,
                                            unsigned short val_nPrimVar, unsigned short val_nPrimVarGrad,
                                            const CConfig* config)
@@ -508,11 +592,11 @@ void CUpwAUSMPLUSUP2_NEMO::ComputeInterfaceQuantities(const CConfig* config, su2
   const su2double ChatR = CstarR * CstarR / max(CstarR, -ProjVelocity_j);
 
   /*--- Interface speed of sound ---*/
-  const su2double aF = min(ChatL, ChatR);
-  interface_soundspeed[0] = interface_soundspeed[1] = aF;
+  const su2double A_F = min(ChatL, ChatR);
+  interface_soundspeed[0] = interface_soundspeed[1] = A_F;
 
-  const su2double M_L = ProjVelocity_i / aF;
-  const su2double M_R = ProjVelocity_j / aF;
+  const su2double M_L = ProjVelocity_i / A_F;
+  const su2double M_R = ProjVelocity_j / A_F;
 
   const su2double rhoF = 0.5 * (Density_i + Density_j);
   const su2double MFsq = 0.5 * (M_L * M_L + M_R * M_R);
@@ -525,7 +609,7 @@ void CUpwAUSMPLUSUP2_NEMO::ComputeInterfaceQuantities(const CConfig* config, su2
   const su2double beta = 1.0 / 8.0;
 
   /*--- Pressure diffusion term ---*/
-  const su2double Mp = -(Kp / fa) * max((1.0 - sigma * MFsq), 0.0) * (Pressure_j - Pressure_i) / (rhoF * aF * aF);
+  const su2double Mp = -(Kp / fa) * max((1.0 - sigma * MFsq), 0.0) * (Pressure_j - Pressure_i) / (rhoF * A_F * A_F);
 
   su2double M_LP, P_LP, M_RM, P_RM;
 
@@ -549,7 +633,7 @@ void CUpwAUSMPLUSUP2_NEMO::ComputeInterfaceQuantities(const CConfig* config, su2
   interface_mach = (M_LP + M_RM + Mp);
 
   /*--- Modified pressure flux ---*/
-  const su2double pFi = sqrt(0.5 * (sq_veli + sq_velj)) * (P_LP + P_RM - 1.0) * 0.5 * (Density_j + Density_i) * aF;
+  const su2double pFi = sqrt(0.5 * (sq_veli + sq_velj)) * (P_LP + P_RM - 1.0) * 0.5 * (Density_j + Density_i) * A_F;
 
   for (auto iDim = 0ul; iDim < nDim; iDim++) {
     pressure[iDim] =
@@ -648,3 +732,80 @@ void CUpwAUSMPLUSM_NEMO::ComputeInterfaceQuantities(const CConfig* config, su2do
         P_un[iDim];
   }
 }
+
+CUpwSLAU_NEMO::CUpwSLAU_NEMO(unsigned short val_nDim, unsigned short val_nVar, unsigned short val_nPrimVar,
+                             unsigned short val_nPrimVarGrad, const CConfig* config, bool val_low_dissipation)
+    : CUpwAUSM_SLAU_Base_NEMO(val_nDim, val_nVar, val_nPrimVar, val_nPrimVarGrad, config) {
+
+  slau_low_diss = val_low_dissipation;
+  slau2 = false;
+}
+
+void CUpwSLAU_NEMO::ComputeInterfaceQuantities(const CConfig* config, su2double* pressure, su2double& interface_mach,
+                                               su2double* interface_soundspeed) {
+
+  /*--- Project velocities and speed of sound ---*/
+
+  const su2double sq_veli = GeometryToolbox::SquaredNorm(nDim, Velocity_i);
+  const su2double sq_velj = GeometryToolbox::SquaredNorm(nDim, Velocity_j);
+
+  su2double Energy_i = Enthalpy_i - Pressure_i/Density_i;
+  su2double Energy_j = Enthalpy_j - Pressure_j/Density_j;
+
+  /*--- Compute interface speed of sound (A_F), and left/right Mach number ---*/
+
+  su2double A_F = 0.5 * (SoundSpeed_i + SoundSpeed_j);
+  su2double M_L = ProjVelocity_i/A_F;
+  su2double M_R = ProjVelocity_j/A_F;
+  interface_soundspeed[0] = interface_soundspeed[1] = A_F;
+
+  /*--- Smooth function of the local Mach number---*/
+
+  su2double Mach_tilde = min(1.0, (1.0/A_F) * sqrt(0.5*(sq_veli+sq_velj)));
+  su2double Chi = pow((1.0 - Mach_tilde),2.0);
+  su2double f_rho = -max(min(M_L,0.0),-1.0) * min(max(M_R,0.0),1.0);
+
+  /*--- Mean normal velocity with density weighting ---*/
+
+  su2double Vn_Mag = (Density_i*fabs(ProjVelocity_i) + Density_j*fabs(ProjVelocity_j)) / (Density_i + Density_j);
+  su2double Vn_MagL= (1.0 - f_rho)*Vn_Mag + f_rho*fabs(ProjVelocity_i);
+  su2double Vn_MagR= (1.0 - f_rho)*Vn_Mag + f_rho*fabs(ProjVelocity_j);
+
+  /*--- Mass flux function ---*/
+
+  su2double mdot = 0.5 * (Density_i*(ProjVelocity_i+Vn_MagL) + Density_j*(ProjVelocity_j-Vn_MagR) - (Chi/A_F)*(Pressure_j-Pressure_i));
+  interface_mach = mdot/A_F;
+
+  /*--- Pressure function ---*/
+
+  su2double BetaL, BetaR, Dissipation_ij;
+
+  if (fabs(M_L) < 1.0) BetaL = 0.25*(2.0-M_L)*pow((M_L+1.0),2.0);
+  else if (M_L >= 0)   BetaL = 1.0;
+  else                 BetaL = 0.0;
+
+  if (fabs(M_R) < 1.0) BetaR = 0.25*(2.0+M_R)*pow((M_R-1.0),2.0);
+  else if (M_R >= 0)   BetaR = 0.0;
+  else                 BetaR = 1.0;
+
+  if (slau_low_diss)
+    Dissipation_ij = GetRoe_Dissipation(Dissipation_i, Dissipation_j, Sensor_i, Sensor_j, config);
+  else
+    Dissipation_ij = 1.0;
+
+  P_F = 0.5*(Pressure_i+Pressure_j) + 0.5*(BetaL-BetaR)*(Pressure_i-Pressure_j);
+
+  if (!slau2) P_F += Dissipation_ij*(1.0-Chi)*(BetaL+BetaR-1.0)*0.5*(Pressure_i+Pressure_j);
+  else        P_F += Dissipation_ij*sqrt(0.5*(sq_veli+sq_velj))*(BetaL+BetaR-1.0)*A_F*0.5*(Density_i+Density_j);
+
+  pressure[0] = pressure[1] = P_F;
+}
+
+CUpwSLAU2_NEMO::CUpwSLAU2_NEMO(unsigned short val_nDim, unsigned short val_nVar, unsigned short val_nPrimVar,
+                               unsigned short val_nPrimVarGrad, const CConfig* config, bool val_low_dissipation)
+          : CUpwSLAU_NEMO(val_nDim, val_nVar, val_nPrimVar, val_nPrimVarGrad, config, val_low_dissipation) {
+
+  /*--- The difference between SLAU and SLAU2 is minimal, so we derive from SLAU and set this flag
+   so that the ComputeIterfaceQuantities function modifies the pressure according to SLAU2. ---*/
+  slau2 = true;
+}
diff --git a/TestCases/nonequilibrium/invwedge/invwedge_ap.cfg b/TestCases/nonequilibrium/invwedge/invwedge_ap.cfg
new file mode 100644
index 00000000000..6837aa25868
--- /dev/null
+++ b/TestCases/nonequilibrium/invwedge/invwedge_ap.cfg
@@ -0,0 +1,72 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%                                                                              %
+% SU2 configuration file                                                       %
+% Case description: Mach 5 inviscid flow over a 10deg wedge - AUSM+-Up         %
+%                                                                              %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
+%
+SOLVER= NEMO_EULER
+GAS_MODEL= AIR-5
+GAS_COMPOSITION= (0.77, 0.23, 0.0, 0.0, 0.0)
+MATH_PROBLEM= DIRECT
+RESTART_SOL= NO
+
+% ----------- COMPRESSIBLE AND INCOMPRESSIBLE FREE-STREAM DEFINITION ----------%
+%
+MACH_NUMBER= 5
+AOA= 0.0
+SIDESLIP_ANGLE= 0.0
+FREESTREAM_PRESSURE= 101325.0
+FREESTREAM_TEMPERATURE= 288.15
+FREESTREAM_TEMPERATURE_VE= 288.15
+
+% ---- NONEQUILIBRIUM GAS, IDEAL GAS, POLYTROPIC, VAN DER WAALS AND PENG ROBINSON CONSTANTS -------%
+%
+FLUID_MODEL= SU2_NONEQ
+
+% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
+%
+MARKER_EULER= ( Euler, Wall )
+MARKER_OUTLET= ( Exit, 10 )
+MARKER_FAR = ( Farfield, Inlet )
+MARKER_PLOTTING= (NONE )
+MARKER_MONITORING= ( Wall )
+
+% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
+%
+NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
+CFL_NUMBER= 0.75
+ITER= 11
+LINEAR_SOLVER= BCGSTAB
+LINEAR_SOLVER_ERROR= 1E-6
+LINEAR_SOLVER_ITER= 5
+
+% -----------------------------------------------------------------------%
+%
+CONV_NUM_METHOD_FLOW= AUSMPLUSUP
+MUSCL_FLOW= YES
+SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
+VENKAT_LIMITER_COEFF= 0.05
+TIME_DISCRE_FLOW= EULER_EXPLICIT
+
+% --------------------------- CONVERGENCE PARAMETERS --------------------------%
+%
+CONV_RESIDUAL_MINVAL= -50
+CONV_STARTITER= 10
+CONV_CAUCHY_ELEMS= 100
+CONV_CAUCHY_EPS= 1E-10
+
+% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
+%
+MESH_FILENAME= invwedge.su2
+MESH_FORMAT= SU2
+SOLUTION_FILENAME= restart_flow.dat  
+TABULAR_FORMAT= TECPLOT 
+CONV_FILENAME= convergence
+RESTART_FILENAME= restart_flow.dat
+VOLUME_FILENAME= soln_volume
+SURFACE_FILENAME= soln_surface     
+OUTPUT_WRT_FREQ= 100
+OUTPUT_FILES= (RESTART_ASCII, PARAVIEW_ASCII)
diff --git a/TestCases/nonequilibrium/invwedge/invwedge_slau.cfg b/TestCases/nonequilibrium/invwedge/invwedge_slau.cfg
new file mode 100644
index 00000000000..b41f2de0044
--- /dev/null
+++ b/TestCases/nonequilibrium/invwedge/invwedge_slau.cfg
@@ -0,0 +1,72 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%                                                                              %
+% SU2 configuration file                                                       %
+% Case description: Mach 5 inviscid flow over a 10deg wedge - SLAU             %
+%                                                                              %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
+%
+SOLVER= NEMO_EULER
+GAS_MODEL= AIR-5
+GAS_COMPOSITION= (0.77, 0.23, 0.0, 0.0, 0.0)
+MATH_PROBLEM= DIRECT
+RESTART_SOL= NO
+
+% ----------- COMPRESSIBLE AND INCOMPRESSIBLE FREE-STREAM DEFINITION ----------%
+%
+MACH_NUMBER= 5
+AOA= 0.0
+SIDESLIP_ANGLE= 0.0
+FREESTREAM_PRESSURE= 101325.0
+FREESTREAM_TEMPERATURE= 288.15
+FREESTREAM_TEMPERATURE_VE= 288.15
+
+% ---- NONEQUILIBRIUM GAS, IDEAL GAS, POLYTROPIC, VAN DER WAALS AND PENG ROBINSON CONSTANTS -------%
+%
+FLUID_MODEL= SU2_NONEQ
+
+% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
+%
+MARKER_EULER= ( Euler, Wall )
+MARKER_OUTLET= ( Exit, 10 )
+MARKER_FAR = ( Farfield, Inlet )
+MARKER_PLOTTING= (NONE )
+MARKER_MONITORING= ( Wall )
+
+% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
+%
+NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
+CFL_NUMBER= 0.75
+ITER= 11
+LINEAR_SOLVER= BCGSTAB
+LINEAR_SOLVER_ERROR= 1E-6
+LINEAR_SOLVER_ITER= 5
+
+% -----------------------------------------------------------------------%
+%
+CONV_NUM_METHOD_FLOW= SLAU
+MUSCL_FLOW= YES
+SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
+VENKAT_LIMITER_COEFF= 0.05
+TIME_DISCRE_FLOW= EULER_EXPLICIT
+
+% --------------------------- CONVERGENCE PARAMETERS --------------------------%
+%
+CONV_RESIDUAL_MINVAL= -50
+CONV_STARTITER= 10
+CONV_CAUCHY_ELEMS= 100
+CONV_CAUCHY_EPS= 1E-10
+
+% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
+%
+MESH_FILENAME= invwedge.su2
+MESH_FORMAT= SU2
+SOLUTION_FILENAME= restart_flow.dat  
+TABULAR_FORMAT= TECPLOT 
+CONV_FILENAME= convergence
+RESTART_FILENAME= restart_flow.dat
+VOLUME_FILENAME= soln_volume
+SURFACE_FILENAME= soln_surface     
+OUTPUT_WRT_FREQ= 100
+OUTPUT_FILES= (RESTART_ASCII, PARAVIEW_ASCII)
diff --git a/TestCases/nonequilibrium/invwedge/invwedge_slau2.cfg b/TestCases/nonequilibrium/invwedge/invwedge_slau2.cfg
new file mode 100644
index 00000000000..98dc868fb54
--- /dev/null
+++ b/TestCases/nonequilibrium/invwedge/invwedge_slau2.cfg
@@ -0,0 +1,72 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%                                                                              %
+% SU2 configuration file                                                       %
+% Case description: Mach 5 inviscid flow over a 10deg wedge - AUSM+M           %
+%                                                                              %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
+%
+SOLVER= NEMO_EULER
+GAS_MODEL= AIR-5
+GAS_COMPOSITION= (0.77, 0.23, 0.0, 0.0, 0.0)
+MATH_PROBLEM= DIRECT
+RESTART_SOL= NO
+
+% ----------- COMPRESSIBLE AND INCOMPRESSIBLE FREE-STREAM DEFINITION ----------%
+%
+MACH_NUMBER= 5
+AOA= 0.0
+SIDESLIP_ANGLE= 0.0
+FREESTREAM_PRESSURE= 101325.0
+FREESTREAM_TEMPERATURE= 288.15
+FREESTREAM_TEMPERATURE_VE= 288.15
+
+% ---- NONEQUILIBRIUM GAS, IDEAL GAS, POLYTROPIC, VAN DER WAALS AND PENG ROBINSON CONSTANTS -------%
+%
+FLUID_MODEL= SU2_NONEQ
+
+% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
+%
+MARKER_EULER= ( Euler, Wall )
+MARKER_OUTLET= ( Exit, 10 )
+MARKER_FAR = ( Farfield, Inlet )
+MARKER_PLOTTING= (NONE )
+MARKER_MONITORING= ( Wall )
+
+% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
+%
+NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
+CFL_NUMBER= 0.75
+ITER= 11
+LINEAR_SOLVER= BCGSTAB
+LINEAR_SOLVER_ERROR= 1E-6
+LINEAR_SOLVER_ITER= 5
+
+% -----------------------------------------------------------------------%
+%
+CONV_NUM_METHOD_FLOW= SLAU2
+MUSCL_FLOW= YES
+SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
+VENKAT_LIMITER_COEFF= 0.05
+TIME_DISCRE_FLOW= EULER_EXPLICIT
+
+% --------------------------- CONVERGENCE PARAMETERS --------------------------%
+%
+CONV_RESIDUAL_MINVAL= -50
+CONV_STARTITER= 10
+CONV_CAUCHY_ELEMS= 100
+CONV_CAUCHY_EPS= 1E-10
+
+% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
+%
+MESH_FILENAME= invwedge.su2
+MESH_FORMAT= SU2
+SOLUTION_FILENAME= restart_flow.dat  
+TABULAR_FORMAT= TECPLOT 
+CONV_FILENAME= convergence
+RESTART_FILENAME= restart_flow.dat
+VOLUME_FILENAME= soln_volume
+SURFACE_FILENAME= soln_surface     
+OUTPUT_WRT_FREQ= 100
+OUTPUT_FILES= (RESTART_ASCII, PARAVIEW_ASCII)
diff --git a/TestCases/parallel_regression.py b/TestCases/parallel_regression.py
index c6e6edf30d5..de5ca79a5b9 100644
--- a/TestCases/parallel_regression.py
+++ b/TestCases/parallel_regression.py
@@ -125,7 +125,7 @@ def main():
     invwedge_msw.test_vals_aarch64 = [-1.165957, -1.690720, -18.301301, -18.628685, -18.574704, 2.151638, 1.721236, 5.193813, 0.751584]
     test_list.append(invwedge_msw)
 
-    # Inviscid single wedge, roe, implicit
+    # Inviscid single wedge, ROE, implicit
     invwedge_roe = TestCase('invwedge_roe')
     invwedge_roe.cfg_dir = "nonequilibrium/invwedge"
     invwedge_roe.cfg_file = "invwedge_roe.cfg"
@@ -134,7 +134,7 @@ def main():
     invwedge_roe.test_vals_aarch64 = [-1.038582, -1.563344, -18.299879, -18.627263, -18.573276, 2.278987, 1.861307, 5.323753, 0.874900]
     test_list.append(invwedge_roe)
 
-    # Inviscid single wedge, lax, implicit
+    # Inviscid single wedge, LAX, implicit
     invwedge_lax = TestCase('invwedge_lax')
     invwedge_lax.cfg_dir = "nonequilibrium/invwedge"
     invwedge_lax.cfg_file = "invwedge_lax.cfg"
@@ -152,6 +152,38 @@ def main():
     invwedge_ausm_m.test_vals_aarch64 = [-1.055079, -1.579842, -16.739725, -17.063618, -17.012831, 2.265433, 1.797600, 5.302743, 0.856551]
     test_list.append(invwedge_ausm_m)
 
+    # Inviscid single wedge, implicit, AUSM+-UP scheme
+    invwedge_ausm_up = TestCase('invwedge_ausm_up')
+    invwedge_ausm_up.cfg_dir = "nonequilibrium/invwedge"
+    invwedge_ausm_up.cfg_file = "invwedge_ap.cfg"
+    invwedge_ausm_up.test_iter = 10
+    invwedge_ausm_up.test_vals = [-1.055083, -1.579845, -16.739725, -17.063618, -17.012831, 2.265430, 1.797602, 5.302740, 0.85654]
+    invwedge_ausm_up.test_vals_aarch64 = [-1.055083, -1.579845, -16.739725, -17.063618, -17.012831, 2.265430, 1.797602, 5.302740, 0.85654]
+    invwedge_ausm_up.new_output = True
+    test_list.append(invwedge_ausm_up)
+
+    # Inviscid single wedge, implicit, SLAU scheme
+    invwedge_slau = TestCase('invwedge_slau')
+    invwedge_slau.cfg_dir = "nonequilibrium/invwedge"
+    invwedge_slau.cfg_file = "invwedge_slau.cfg"
+    invwedge_slau.test_iter = 10
+    invwedge_slau.test_vals = [-1.055083, -1.579845, -16.739725, -17.063618, -17.012831, 2.265430, 1.797602, 5.302740, 0.85654]
+    invwedge_slau.test_vals_aarch64 = [-1.055083, -1.579845, -16.739725, -17.063618, -17.012831, 2.265430, 1.797602, 5.302740, 0.85654]
+    invwedge_slau.new_output = True
+    test_list.append(invwedge_slau)
+
+    # Inviscid single wedge, implicit, SLAU2 scheme
+    invwedge_slau2 = TestCase('invwedge_slau2')
+    invwedge_slau2.cfg_dir = "nonequilibrium/invwedge"
+    invwedge_slau2.cfg_file = "invwedge_slau2.cfg"
+    invwedge_slau2.test_iter = 10
+    invwedge_slau2.test_vals = [-1.055083, -1.579845, -16.739725, -17.063618, -17.012831, 2.265430, 1.797602, 5.302740, 0.85654]
+    invwedge_slau2.test_vals_aarch64 = [-1.055083, -1.579845, -16.739725, -17.063618, -17.012831, 2.265430, 1.797602, 5.302740, 0.85654]
+    invwedge_slau2.new_output = True
+    test_list.append(invwedge_slau2)
+
+
+
     # Inviscid single wedge, implicit, NEMO supersonic inlet
     invwedge_ss_inlet = TestCase('invwedge_ss_inlet')
     invwedge_ss_inlet.cfg_dir = "nonequilibrium/invwedge"