Merge pull request #660 from ketch/fix_limiting_bug

Correct labels in philim.f90 select statement.
clawpack · Jan 31, 2021 · 97789e5 · 97789e5
2 parents 59fe187 + f90b088
commit 97789e5
Showing 1 changed file with 5 additions and 5 deletions.
diff --git a/src/pyclaw/classic/philim.f90 b/src/pyclaw/classic/philim.f90
@@ -16,20 +16,20 @@ double precision function philim(a,b,meth)
     r = b/a
     select case (meth)
 
-    case (10) ! minmod
+    case (1) ! minmod
         philim = dmax1(0.d0, dmin1(1.d0, r))
 
-    case (20) ! superbee
+    case (2) ! superbee
         philim = dmax1(0.d0, dmin1(1.d0, 2.d0*r), dmin1(2.d0, r))
 
-    case (30) ! van Leer
+    case (3) ! van Leer
         philim = (r + dabs(r)) / (1.d0 + dabs(r))
 
-    case (40) ! monotonized centered
+    case (4) ! monotonized centered
         c = (1.d0 + r)/2.d0
         philim = dmax1(0.d0, dmin1(c, 2.d0, 2.d0*r))
 
-    case (50) ! Beam-Warming
+    case (5) ! Beam-Warming
         philim = r
 
     end select