+int process(t_moldyn *moldyn,t_atom *itom,t_atom *jtom,void *ptr,u8 bc) {
+
+ t_3dvec dist;
+ double d;
+ t_data *data;
+
+ data=ptr;
+
+ v3_sub(&dist,&(itom->r),&(jtom->r));
+ if(bc) check_per_bound(moldyn,&dist);
+ d=v3_norm(&dist);
+
+ if(d<=data->radius) {
+ printf("%s %f %f %f %s %f\n",
+ pse_name[jtom->element],
+ jtom->r.x+data->ox,jtom->r.y+data->oy,jtom->r.z+data->oz,
+ (jtom->tag==data->ma)?"Red":pse_col[jtom->element],
+ jtom->ekin);
+ }
+
+ return 0;
+}
+