[FLASH-USERS] Bug in Particle timestep routine

[Klaus Weide]

On Tue, 5 May 2015, Thomas Peters wrote:

> Dear FLASH developers,
> 
> in Particles_computeDt.F90, the lines
> 
>         if (NDIM > 2) then
>            velzabs = abs(particles(VELY_PART_PROP,i))
>            if (velzabs > pt_small) then
>               dtz = delta(3) / velzabs
>            endif
> 
>         end if
> 
> should read
> 
>         if (NDIM > 2) then
>            velzabs = abs(particles(VELZ_PART_PROP,i))
>            if (velzabs > pt_small) then
>               dtz = delta(3) / velzabs
>            endif
> 
>         end if
> 
> Otherwise the z-component of the particle velocity is not correctly
> taken into account when the particle timestep is computed.

Agreed; your fix will be applied.

Thank you!

Klaus