add program to read pos,vel,mass,u