From 9c49fd8cc1b2c173bfb18699a5af916838132a29 Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Thu, 23 Jan 2014 10:24:53 -0200 Subject: [PATCH] SCHEMES: Revert previous flux for Sm=0 and Bx=0. Signed-off-by: Grzegorz Kowal --- src/schemes.F90 | 14 +++++++++++--- 1 file changed, 11 insertions(+), 3 deletions(-) diff --git a/src/schemes.F90 b/src/schemes.F90 index 92c77f9..125336d 100644 --- a/src/schemes.F90 +++ b/src/schemes.F90 @@ -3589,10 +3589,18 @@ module schemes else ! sm = 0 -! intermediate flux is constant across the contact discontinuity and since there -! are no Alfven waves, we can calculate the flux using the HLL formula +! intermediate flux is constant across the contact discontinuity and all except +! the parallel momentum flux are zero ! - f(:,i) = (sl * wr(:) - sr * wl(:)) / srml + f(idn,i) = 0.0d+00 + f(imx,i) = - 0.5d+00 * (wl(imx) + wr(imx)) + f(imy,i) = 0.0d+00 + f(imz,i) = 0.0d+00 + f(ibx,i) = fl(ibx,i) + f(iby,i) = 0.0d+00 + f(ibz,i) = 0.0d+00 + f(ibp,i) = fl(ibp,i) + f(ien,i) = 0.0d+00 end if ! sm = 0