Thank you, does it mean that if I add after line number 466 of dynamics_module.f90 something like
tau_old(1,100) = tau(1,100) - v_x*dt/alat tau_old(2,100) = tau(2,100) - v_y*dt/alat tau_old(3,100) = tau(3,100) - v_z*dt/alat it should be equivalent to set up initial velocity (v_x,v_y,v_z) to particle number 100? What units is used for dt? Is it atomic Hartree or atomic Rydberg?
