From 9d75c49059ba91371ebcd2ba39a62f8f058c268d Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Tue, 16 May 2017 17:03:01 -0300 Subject: [PATCH] SCHEMES: Cosmetic changes in the adiabatic HLLD solver. Signed-off-by: Grzegorz Kowal --- src/schemes.F90 | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/src/schemes.F90 b/src/schemes.F90 index f698c55..2dacb90 100644 --- a/src/schemes.F90 +++ b/src/schemes.F90 @@ -3642,13 +3642,13 @@ module schemes sl = min(clm(i), crm(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) - else if (sr <= 0.0d+00) then + else if (sr <= 0.0d+00) then ! sr ≤ 0 f(:,i) = fr(:,i) @@ -3837,7 +3837,7 @@ module schemes ui(iby) = by ui(ibz) = bz 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 ! @@ -3855,7 +3855,7 @@ module schemes ui(iby) = by ui(ibz) = bz 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 ! @@ -3932,7 +3932,7 @@ module schemes ui(iby) = by ui(ibz) = bz 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 ! advection speed @@ -4013,7 +4013,7 @@ module schemes ui(iby) = by ui(ibz) = bz 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 ! advection speed