gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

RECONNECTION: Normalize velocity noise.

Browse Source 
Normalize the velocity so it has the same amplitude everywhere, and
random direction.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2017-03-28 16:06:15 -03:00
parent 9c311f2413
commit c7e0cb92d2
1 changed files with 19 additions and 3 deletions

22
src/user_problem.F90

@ -253,6 +253,7 @@ module user_problem
integer :: i, j, k
real(kind=8) :: yt, yp, yl, yu
real(kind=8) :: yrat, itanh
real(kind=8) :: vx, vy, vz, vv
! local arrays
!
@ -387,10 +388,25 @@ module user_problem
do j = 1, jm
if (abs(y(j)) <= ycut) then
q(ivx,j) = vper * randomn()
q(ivy,j) = vper * randomn()
vv = 0.0d+00
do while(vv == 0.0d+00)
vx = randomn()
vy = randomn()
#if NDIMS == 3
q(ivz,j) = vper * randomn()
vz = randomn()
#endif /* NDIMS == 3 */
#if NDIMS == 3
vv = sqrt(vx * vx + vy * vy + vz * vz)
#else /* NDIMS == 3 */
vv = sqrt(vx * vx + vy * vy)
#endif /* NDIMS == 3 */
end do
q(ivx,j) = vper * (vx / vv)
q(ivy,j) = vper * (vy / vv)
#if NDIMS == 3
q(ivz,j) = vper * (vz / vv)
#endif /* NDIMS == 3 */
end if ! |y| < ycut