int add_atom(t_moldyn *moldyn,int element,double mass,u8 bnum,u8 attr,
t_3dvec *r,t_3dvec *v);
int destroy_atoms(t_moldyn *moldyn);
int thermal_init(t_moldyn *moldyn,u8 equi_init);
int scale_velocity(t_moldyn *moldyn,u8 equi_init);
int add_atom(t_moldyn *moldyn,int element,double mass,u8 bnum,u8 attr,
t_3dvec *r,t_3dvec *v);
int destroy_atoms(t_moldyn *moldyn);
int thermal_init(t_moldyn *moldyn,u8 equi_init);
int scale_velocity(t_moldyn *moldyn,u8 equi_init);