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

SCHEMES: Cosmetic changes in the adiabatic HLLD solver.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2017-05-16 17:03:01 -03:00
parent 4fd495bf80
commit 9d75c49059
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/schemes.F90
View File

@ -3642,13 +3642,13 @@ module schemes
sl = min(clm(i), crm(i)) sl = min(clm(i), crm(i))
sr = max(clp(i), crp(i)) sr = max(clp(i), crp(i))
! calculate the HLLC flux ! calculate the HLLD flux
! !
if (sl >= 0.0d+00) then if (sl >= 0.0d+00) then ! sl ≥ 0
f(:,i) = fl(:,i) f(:,i) = fl(:,i)
else if (sr <= 0.0d+00) then else if (sr <= 0.0d+00) then ! sr ≤ 0
f(:,i) = fr(:,i) f(:,i) = fr(:,i)
@ -3837,7 +3837,7 @@ module schemes
ui(iby) = by ui(iby) = by
ui(ibz) = bz ui(ibz) = bz
ui(ibp) = ul(ibp,i) ui(ibp) = ul(ibp,i)
ui(ien) = (wcl(ien) + sm * pt - bx * vb) / (sml - sm) ui(ien) = - (wcl(ien) + sm * pt - bx * vb) / cal
! the inmost left intermediate flux ! the inmost left intermediate flux
! !
@ -3855,7 +3855,7 @@ module schemes
ui(iby) = by ui(iby) = by
ui(ibz) = bz ui(ibz) = bz
ui(ibp) = ur(ibp,i) ui(ibp) = ur(ibp,i)
ui(ien) = (wcr(ien) + sm * pt - bx * vb) / (smr - sm) ui(ien) = (wcr(ien) + sm * pt - bx * vb) / car
! the inmost right intermediate flux ! the inmost right intermediate flux
! !
@ -3932,7 +3932,7 @@ module schemes
ui(iby) = by ui(iby) = by
ui(ibz) = bz ui(ibz) = bz
ui(ibp) = ul(ibp,i) ui(ibp) = ul(ibp,i)
ui(ien) = (wcl(ien) + sm * pt - bx * vb) / (sml - sm) ui(ien) = - (wcl(ien) + sm * pt - bx * vb) / cal
! choose the correct state depending on the sign of contact discontinuity ! choose the correct state depending on the sign of contact discontinuity
! advection speed ! advection speed
@ -4013,7 +4013,7 @@ module schemes
ui(iby) = by ui(iby) = by
ui(ibz) = bz ui(ibz) = bz
ui(ibp) = ur(ibp,i) ui(ibp) = ur(ibp,i)
ui(ien) = (wcr(ien) + sm * pt - bx * vb) / (smr - sm) ui(ien) = (wcr(ien) + sm * pt - bx * vb) / car
! choose the correct state depending on the sign of contact discontinuity ! choose the correct state depending on the sign of contact discontinuity
! advection speed ! advection speed