From 775fe94b3b32ae4f5209b9cc2bfb16b294c1c2ee Mon Sep 17 00:00:00 2001 From: Jannis Teunissen Date: Fri, 20 Mar 2020 14:31:49 +0100 Subject: [PATCH] Fix bug in Kurganov-Tadmor scheme (factor 0.5 for diffusive flux) --- src/m_af_flux_schemes.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/m_af_flux_schemes.f90 b/src/m_af_flux_schemes.f90 index 5a6d384..65ff690 100644 --- a/src/m_af_flux_schemes.f90 +++ b/src/m_af_flux_schemes.f90 @@ -214,7 +214,7 @@ contains real(dp), intent(in) :: wmax(n_values) real(dp), intent(inout) :: flux(n_values, n_vars) - flux = 0.5_dp * (flux_l + flux_r) - spread(wmax, 2, n_vars) * (u_r - u_l) + flux = 0.5_dp * (flux_l + flux_r - spread(wmax, 2, n_vars) * (u_r - u_l)) end subroutine flux_kurganovTadmor_1d subroutine flux_update_densities(tree, dt, n_vars, i_cc, i_flux, & -- GitLab