diff --git a/compressible_flow/Turbulent_ONERAM6/turb_ONERAM6.cfg b/compressible_flow/Turbulent_ONERAM6/turb_ONERAM6.cfg
index 88aaafc..bb50a15 100644
--- a/compressible_flow/Turbulent_ONERAM6/turb_ONERAM6.cfg
+++ b/compressible_flow/Turbulent_ONERAM6/turb_ONERAM6.cfg
@@ -185,22 +185,9 @@ MG_DAMP_PROLONGATION= 0.7
 
 % -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
 %
-% Convective numerical method (JST, LAX-FRIEDRICH, CUSP, ROE, AUSM, HLLC,
-%                              TURKEL_PREC, MSW)
+% Convective numerical method
 CONV_NUM_METHOD_FLOW= JST
 %
-% Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER)
-MUSCL_FLOW= YES
-%
-% Slope limiter (NONE, VENKATAKRISHNAN, VENKATAKRISHNAN_WANG,
-%                BARTH_JESPERSEN, VAN_ALBADA_EDGE)
-SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
-%
-% Coefficient for the Venkat's limiter (upwind scheme). A larger values decrease
-%             the extent of limiting, values approaching zero cause
-%             lower-order approximation to the solution (0.05 by default)
-VENKAT_LIMITER_COEFF= 0.05
-%
 % 2nd and 4th order artificial dissipation coefficients for
 %     the JST method ( 0.5, 0.02 by default )
 JST_SENSOR_COEFF= ( 0.5, 0.02 )
diff --git a/design/Turbulent_2D_Constrained_RAE2822/turb_SA_RAE2822.cfg b/design/Turbulent_2D_Constrained_RAE2822/turb_SA_RAE2822.cfg
index 28aa434..2077016 100644
--- a/design/Turbulent_2D_Constrained_RAE2822/turb_SA_RAE2822.cfg
+++ b/design/Turbulent_2D_Constrained_RAE2822/turb_SA_RAE2822.cfg
@@ -109,7 +109,7 @@ OBJECTIVE_FUNCTION= DRAG
 %
 % Monotonic Upwind Scheme for Conservation Laws (TVD) in the flow equations.
 %           Required for 2nd order upwind schemes (NO, YES)
-MUSCL_FLOW= YES
+MUSCL_FLOW= NO
 %
 % Slope limiter (NONE, VENKATAKRISHNAN, VENKATAKRISHNAN_WANG,
 %                BARTH_JESPERSEN, VAN_ALBADA_EDGE)
diff --git a/design/Unsteady_Shape_Opt_NACA0012/unsteady_naca0012_opt.cfg b/design/Unsteady_Shape_Opt_NACA0012/unsteady_naca0012_opt.cfg
index 284266b..b04aa84 100644
--- a/design/Unsteady_Shape_Opt_NACA0012/unsteady_naca0012_opt.cfg
+++ b/design/Unsteady_Shape_Opt_NACA0012/unsteady_naca0012_opt.cfg
@@ -136,10 +136,6 @@ LINEAR_SOLVER_ITER= 5
 %                              TURKEL_PREC, MSW)
 CONV_NUM_METHOD_FLOW= JST
 %
-% Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER)
-%
-% MUSCL_FLOW= YES
-%
 % 1st, 2nd and 4th order artificial dissipation coefficients
 JST_SENSOR_COEFF= (  0.5, 0.01 )
 %
diff --git a/multiphysics/unsteady_fsi_python/Ma01/fluid.cfg b/multiphysics/unsteady_fsi_python/Ma01/fluid.cfg
index 15dbd83..bcfb42c 100644
--- a/multiphysics/unsteady_fsi_python/Ma01/fluid.cfg
+++ b/multiphysics/unsteady_fsi_python/Ma01/fluid.cfg
@@ -121,18 +121,10 @@ LINEAR_SOLVER_ITER= 10
 %
 % -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
 %
-% Convective numerical method (JST, LAX-FRIEDRICH, CUSP, ROE, AUSM, HLLC,
+% Convective numerical method (JST, LAX-FRIEDRICH, ROE, AUSM, HLLC,
 %                              TURKEL_PREC, MSW)
 CONV_NUM_METHOD_FLOW= JST
 %
-% Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER)
-%
-% Monotonic Upwind Scheme for Conservation Laws (TVD) in the flow equations.
-%           Required for 2nd order upwind schemes (NO, YES)
-MUSCL_FLOW= YES
-% Slope limiter (VENKATAKRISHNAN, MINMOD)
-SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
-%
 JST_SENSOR_COEFF= ( 0.5, 0.01 )
 % Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
 TIME_DISCRE_FLOW= EULER_IMPLICIT
diff --git a/multiphysics/unsteady_fsi_python/Ma02/fluid.cfg b/multiphysics/unsteady_fsi_python/Ma02/fluid.cfg
index 0ad3f6f..80bd6f1 100644
--- a/multiphysics/unsteady_fsi_python/Ma02/fluid.cfg
+++ b/multiphysics/unsteady_fsi_python/Ma02/fluid.cfg
@@ -125,14 +125,6 @@ LINEAR_SOLVER_ITER= 10
 %                              TURKEL_PREC, MSW)
 CONV_NUM_METHOD_FLOW= JST
 %
-% Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER)
-%
-% Monotonic Upwind Scheme for Conservation Laws (TVD) in the flow equations.
-%           Required for 2nd order upwind schemes (NO, YES)
-MUSCL_FLOW= YES
-% Slope limiter (VENKATAKRISHNAN, MINMOD)
-SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
-%
 JST_SENSOR_COEFF= ( 0.5, 0.01 )
 % Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
 TIME_DISCRE_FLOW= EULER_IMPLICIT
diff --git a/multiphysics/unsteady_fsi_python/Ma03/fluid.cfg b/multiphysics/unsteady_fsi_python/Ma03/fluid.cfg
index d7b86c5..4913d24 100644
--- a/multiphysics/unsteady_fsi_python/Ma03/fluid.cfg
+++ b/multiphysics/unsteady_fsi_python/Ma03/fluid.cfg
@@ -125,14 +125,6 @@ LINEAR_SOLVER_ITER= 10
 %                              TURKEL_PREC, MSW)
 CONV_NUM_METHOD_FLOW= JST
 %
-% Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER)
-%
-% Monotonic Upwind Scheme for Conservation Laws (TVD) in the flow equations.
-%           Required for 2nd order upwind schemes (NO, YES)
-MUSCL_FLOW= YES
-% Slope limiter (VENKATAKRISHNAN, MINMOD)
-SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
-%
 JST_SENSOR_COEFF= ( 0.5, 0.01 )
 % Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
 TIME_DISCRE_FLOW= EULER_IMPLICIT
diff --git a/multiphysics/unsteady_fsi_python/Ma0357/fluid.cfg b/multiphysics/unsteady_fsi_python/Ma0357/fluid.cfg
index baa2370..f2d14ea 100644
--- a/multiphysics/unsteady_fsi_python/Ma0357/fluid.cfg
+++ b/multiphysics/unsteady_fsi_python/Ma0357/fluid.cfg
@@ -125,14 +125,6 @@ LINEAR_SOLVER_ITER= 10
 %                              TURKEL_PREC, MSW)
 CONV_NUM_METHOD_FLOW= JST
 %
-% Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER)
-%
-% Monotonic Upwind Scheme for Conservation Laws (TVD) in the flow equations.
-%           Required for 2nd order upwind schemes (NO, YES)
-MUSCL_FLOW= YES
-% Slope limiter (VENKATAKRISHNAN, MINMOD)
-SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
-%
 JST_SENSOR_COEFF= ( 0.5, 0.01 )
 % Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
 TIME_DISCRE_FLOW= EULER_IMPLICIT
diff --git a/multiphysics/unsteady_fsi_python/Ma0364/fluid.cfg b/multiphysics/unsteady_fsi_python/Ma0364/fluid.cfg
index 8d1584c..6c3f9c7 100644
--- a/multiphysics/unsteady_fsi_python/Ma0364/fluid.cfg
+++ b/multiphysics/unsteady_fsi_python/Ma0364/fluid.cfg
@@ -125,14 +125,6 @@ LINEAR_SOLVER_ITER= 10
 %                              TURKEL_PREC, MSW)
 CONV_NUM_METHOD_FLOW= JST
 %
-% Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER)
-%
-% Monotonic Upwind Scheme for Conservation Laws (TVD) in the flow equations.
-%           Required for 2nd order upwind schemes (NO, YES)
-MUSCL_FLOW= YES
-% Slope limiter (VENKATAKRISHNAN, MINMOD)
-SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
-%
 JST_SENSOR_COEFF= ( 0.5, 0.01 )
 % Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
 TIME_DISCRE_FLOW= EULER_IMPLICIT