From 403a1e2f3f9cbcac3a6b6a123ee4d32dcc4bfe53 Mon Sep 17 00:00:00 2001
From: Grzegorz Kowal <grzegorz@amuncode.org>
Date: Sat, 20 Apr 2024 17:47:19 -0300
Subject: [PATCH] INTERPOLATION: kappa parameter should be >= 1.

The kappa parameter is calculated from CFL value. It is a decreasing
value with increasing CFL. However, to maintain the precision of MP
limiter, it should be kept above 1. Therefore, instead of taking the min
values between the user defined kappa and the one calculated from the
CFL value, take its maximum.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
---
 sources/interpolations.F90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/sources/interpolations.F90 b/sources/interpolations.F90
index 778de4b..edb5a3a 100644
--- a/sources/interpolations.F90
+++ b/sources/interpolations.F90
@@ -224,7 +224,7 @@ module interpolations
 
 ! calculate κ = (1 - ν) / ν
 !
-    kappa = min(kappa, (1.0d+00 - cfl) / cfl)
+    kappa = max(kappa, (1.0d+00 - cfl) / cfl)
 
 ! correct central weight
 !