brotching the code

BloubWorld/bloubspace.f90

@@ -28,9 +28,9 @@ module bloubspace
     blb%py = rand() * 0.25
     blb%pz = rand() - 0.50
 
-    blb%vx = (rand() - 0.5) / 5.000
-    blb%vy = (rand() - 0.5) / 5.000
-    blb%vz = (rand() - 0.5) / 5.000
+    blb%vx = (rand() - 0.5) / 4.000
+    blb%vy = (rand() - 0.5) / 4.000
+    blb%vz = (rand() - 0.5) / 4.000
 
   end subroutine
 
@@ -40,6 +40,8 @@ module bloubspace
     type(t_bloubs), intent (inout)  :: blb
     real,           intent (in)     :: coef
 
+    ! we must check that this bloub is alive ?
+
     blb%px  = blb%px + (blb%vx * coef)
     blb%py  = blb%py + (blb%vy * coef)
     blb%pz  = blb%pz + (blb%vz * coef)

BloubWorld/genbloubs.f90

@@ -31,7 +31,7 @@ program genbloubs
     bloub%nick = 'noname  '
     bloub%alive = .TRUE.
     call random_pv(bloub)
-    bloub%radius = 0.009999
+    bloub%radius = 0.025
     bloub%seq    = 0
 
     write(idu) bloub              ! no error control ?

BloubWorld/movebloubs.f90

@@ -44,7 +44,7 @@ program movebloubs
       exit
     endif
 
-    call move_bloub (bloub, 0.333)
+    call move_bloub (bloub, 1.10)
 
     bx = bx + bloub%px
     by = by + bloub%py
@@ -55,6 +55,10 @@ program movebloubs
       bloub%vy = -1.0 * bloub%vy
       bloub%py = 0.0
     endif
+    if (3.0 .lt. bloub%py) then
+      bloub%vy = -1.0 * bloub%vy
+      bloub%py = 3.0
+    endif
 
     write(outu) bloub     ! no error control ?