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

USER_PROBLEM: Use tx, ty, tz to determine kx, ky, kz.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2021-08-11 12:33:07 -03:00
parent 96afcefe49
commit c08121bd5b
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
sources/user_problem.F90
View File

@ -304,13 +304,13 @@ module user_problem
thh = pi2 * randuni()
#if NDIMS == 3
thv = pi * randsym()
ux(n) = cos(thh) * cos(thv)
uy(n) = sin(thh) * cos(thv)
uz(n) = sin(thv)
kx(n) = pi2 * nint(kper * ux(n))
ky(n) = pi2 * nint(kper * uy(n))
kz(n) = pi2 * nint(kper * uz(n))
tt = 0.0d+00
tx = cos(thh) * cos(thv)
ty = sin(thh) * cos(thv)
tz = sin(thv)
kx(n) = pi2 * nint(kper * tx)
ky(n) = pi2 * nint(kper * ty)
kz(n) = pi2 * nint(kper * tz)
tt = 0.0d+00
do while(tt < 1.0d-08)
thh = pi2 * randuni()
thv = pi * randsym()