From db40c6253cf944600a8d1ee0c9a6dfe9e168ee29 Mon Sep 17 00:00:00 2001
From: Grzegorz Kowal <grzegorz@gkowal.info>
Date: Wed, 1 Dec 2010 12:39:56 -0200
Subject: [PATCH] SCHEME: fix iteration for reconstruction in HLL.
- fix the index for reconstruction of magnetic field components in the
subroutine HLL;
---
src/scheme.F90 | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/src/scheme.F90 b/src/scheme.F90
index bf67fd5..8c88215 100644
--- a/src/scheme.F90
+++ b/src/scheme.F90
@@ -705,7 +705,7 @@ module scheme
#ifdef GLM
! reconstruct the left and right states of the magnetic field components
!
- do p = ibx, iby
+ do p = ibx, ibz
call reconstruct(n, q(p,:), ql(p,:), qr(p,:))
end do
@@ -744,8 +744,7 @@ module scheme
ap = ar * al
div = 1.0 / (ar - al)
- fn(1:nqt,i) = div * (ar * fl(1:nqt,i) - al * fr(1:nqt,i) &
- + ap * (ur(1:nqt,i) - ul(1:nqt,i)))
+ fn(:,i) = div * (ar * fl(:,i) - al * fr(:,i) + ap * (ur(:,i) - ul(:,i)))
end if
end do