From d685d1bc00113c5d1b30951851e698c96427c3e0 Mon Sep 17 00:00:00 2001 From: Lorenzo Cipriani Date: Tue, 19 Sep 2023 16:05:30 +0200 Subject: [PATCH] Fixed atan initial data --- src/initial_data.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/initial_data.f90 b/src/initial_data.f90 index e8f95b2..efa898c 100644 --- a/src/initial_data.f90 +++ b/src/initial_data.f90 @@ -8,7 +8,7 @@ subroutine initial_data(NX, x, dx, m, r0, a0, idx, u) real(RK), dimension(NX), intent(in) :: x real(RK), dimension(2*NX), intent(out):: u - integer :: i, j + integer :: i real(RK) :: th, Mass, heaviside, p, tmp real(RK), dimension(NX) :: rr, Marray @@ -57,7 +57,7 @@ subroutine initial_data(NX, x, dx, m, r0, a0, idx, u) ! Full dynamics, starting from atan density profile else if ( idx .eq. 3 ) then - + ! Unnormalized density function do i = 1, NX rr(i) = PI/2_RK - atan(x(i) - r0)