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