From a896b7188a6baa970de36edeecd16401dbfe8709 Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Tue, 29 Aug 2023 18:43:44 -0300 Subject: [PATCH] USER_PROBLEM: Do not limit B at boundaries in profile 2. Signed-off-by: Grzegorz Kowal --- sources/user_problem.F90 | 9 ++------- 1 file changed, 2 insertions(+), 7 deletions(-) diff --git a/sources/user_problem.F90 b/sources/user_problem.F90 index d71b7da..178f6d1 100644 --- a/sources/user_problem.F90 +++ b/sources/user_problem.F90 @@ -310,13 +310,8 @@ module user_problem do j = 1, nn do i = 1, nn - if (abs(x(i)) <= 5.0d-01) then - bx = bamp * csx(i) * (csy(j) * thy(j) + sny(j) / (pi2 * dlta * ch2(j))) - by = bamp * 5.0d-01 * snx(i) * sny(j) * thy(j) - else - bx = 0.0d+00 - by = 0.0d+00 - end if + bx = bamp * csx(i) * (csy(j) * thy(j) + sny(j) / (pi2 * dlta * ch2(j))) + by = bamp * 5.0d-01 * snx(i) * sny(j) * thy(j) bz = sqrt(max(0.0d+00, bgui**2 + zeta * (bamp**2 - bx**2 - by**2))) pdata%q(ibx,i,j,:) = bx