From e52fb7c633826d907f1f0b762763c14a132a248e Mon Sep 17 00:00:00 2001
From: Grzegorz Kowal <grzegorz@amuncode.org>
Date: Sun, 2 Jan 2022 17:25:04 -0300
Subject: [PATCH] SCHEMES: Fix KEPES to not produce non-zero Vz and Bz in 2D.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
---
sources/schemes.F90 | 10 +++++-----
1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/sources/schemes.F90 b/sources/schemes.F90
index 54d6ecc..c916316 100644
--- a/sources/schemes.F90
+++ b/sources/schemes.F90
@@ -3582,8 +3582,8 @@ module schemes
x2 = b2 / bp
x3 = b3 / bp
else
- x2 = 1.0d+00 / sqrt(2.0d+00)
- x3 = x2
+ x2 = 0.0d+00
+ x3 = 0.0d+00
end if
ca2 = b1**2
@@ -3764,9 +3764,9 @@ module schemes
rm(:,:) = 0.0d+00
rm(6,4) = 1.0d+00
rm(9,4) = 1.0d+00
+ rm(1,5) = 1.0d+00
rm(6,6) = 1.0d+00
rm(9,6) = - 1.0d+00
- rm(1,5) = 1.0d+00
first = .false.
end if
@@ -3843,8 +3843,8 @@ module schemes
x2 = b2 / bp
x3 = b3 / bp
else
- x2 = 1.0d+00 / sqrt(2.0d+00)
- x3 = x2
+ x2 = 0.0d+00
+ x3 = 0.0d+00
end if
ca2 = b1**2