gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

SCHEMES: Revert previous flux for Sm=0 and Bx=0.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2014-01-23 10:24:53 -02:00
parent a44c0c25c6
commit 9c49fd8cc1
1 changed files with 11 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/schemes.F90
View File

@ -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