/*velset set thermal velocities for all particles in gas */
void velset(void);