From cf1323e1a8d2b608f17dd56037b740e11aae24fe Mon Sep 17 00:00:00 2001 From: tth Date: Fri, 18 Mar 2022 22:06:12 +0100 Subject: [PATCH] better random_pv --- BloubWorld/bloubspace.f90 | 9 +++++++-- 1 file changed, 7 insertions(+), 2 deletions(-) diff --git a/BloubWorld/bloubspace.f90 b/BloubWorld/bloubspace.f90 index 135c1a6..78d8b3c 100644 --- a/BloubWorld/bloubspace.f90 +++ b/BloubWorld/bloubspace.f90 @@ -72,12 +72,17 @@ module bloubspace type(t_bloubs), intent (out) :: blb blb%px = 4.33 * (rand() - 0.50) - blb%py = 3.33 * (rand() - 0.50) + blb%py = 3.70 * (rand() - 0.50) blb%pz = 4.51 * (rand() - 0.50) blb%vx = (rand()) / 2.000 - blb%vy = (rand()) / 2.900 + if (blb%px .LT. 0.0) blb%vx = -blb%vx + + blb%vy = (rand()) / 3.200 + if (blb%py .LT. 0.0) blb%vy = -blb%vx + blb%vz = (rand()) / 2.000 + if (blb%pz .LT. 0.0) blb%vz = -blb%vz blb%state = 0 blb%alive = .TRUE.