From ab481646b6aac475c19207dc130629b29ec2809d Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Fri, 16 Aug 2024 16:30:31 -0300 Subject: [PATCH] USER_PROBLEM: Adapt the sincos profile for the periodic box. Signed-off-by: Grzegorz Kowal --- sources/user_problem.F90 | 15 +++++++++------ 1 file changed, 9 insertions(+), 6 deletions(-) diff --git a/sources/user_problem.F90 b/sources/user_problem.F90 index 0dbd05e..b19c6c5 100644 --- a/sources/user_problem.F90 +++ b/sources/user_problem.F90 @@ -290,7 +290,7 @@ module user_problem integer :: i, j, k real(kind=8) :: r1, r2, bx, by, bz real(kind=8), dimension(nn) :: x, y - real(kind=8), dimension(nn) :: snx, csx, sny, csy, thy, ch2 + real(kind=8), dimension(nn) :: snx, csx, sny, csy, thy, shy real(kind=8), dimension(nn,nn) :: az !------------------------------------------------------------------------------- @@ -301,16 +301,19 @@ module user_problem select case(profile) case(2) - snx(:) = sin(pi * x(:)) - csx(:) = cos(pi * x(:)) + snx(:) = sin(pi * x(:)) + csx(:) = cos(pi * x(:)) sny(:) = sin(pi2 * y(:)) csy(:) = cos(pi2 * y(:)) - thy(:) = tanh(y(:) / dlta) - ch2(:) = cosh(y(:) / dlta)**2 + thy(:) = tanh( y(:) / dlta) & + - tanh((y(:) + 1.0d+00) / dlta) & + + 1.0d+00 + shy(:) = tanh((y(:) + 1.0d+00) / dlta)**2 / dlta & + - tanh( y(:) / dlta)**2 / dlta do j = 1, nn do i = 1, nn - bx = bamp * csx(i) * (csy(j) * thy(j) + sny(j) / (pi2 * dlta * ch2(j))) + bx = bamp * csx(i) * (csy(j) * thy(j) + sny(j) * shy(j) / pi2) 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)))