From 9b7f2f1044e3dff5244d9a35f8880b7a275d635c Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Fri, 5 Nov 2021 23:17:24 -0300 Subject: [PATCH] USER_PROBLEM: Fix a mistake in the absoption layer profile. Signed-off-by: Grzegorz Kowal --- sources/user_problem.F90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/sources/user_problem.F90 b/sources/user_problem.F90 index 6b4b377..b37252b 100644 --- a/sources/user_problem.F90 +++ b/sources/user_problem.F90 @@ -1111,9 +1111,9 @@ module user_problem fc(:) = (abs(yc(:)) - yabs) / adec do j = 1, nn - if (fc(j) >= 0.0d+00) then - if (fc(j) <= 1.0d+00) then - fr = fa * sin(5.0d-01 * pi * yc(j))**2 + if (fc(j) > 0.0d+00) then + if (fc(j) < 1.0d+00) then + fr = fa * sin(5.0d-01 * pi * fc(j))**2 else fr = fa end if