[FLASH-USERS] Bug in Particle timestep routine

[Thomas Peters]

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.

Best regards,
Thomas