scaling to fs, angstrom, amu done, probs with velocity scaling!
[physik/posic.git] / visual / visual.c
index 687c561..d03903b 100644 (file)
@@ -74,11 +74,11 @@ int visual_atoms(t_visual *v,double time,t_atom *atom,int n) {
        char file[128+64];
        t_3dvec dim;
 
-       dim.x=10e9*v->dim.x;
-       dim.y=10e9*v->dim.y;
-       dim.z=10e9*v->dim.z;
+       dim.x=v->dim.x;
+       dim.y=v->dim.y;
+       dim.z=v->dim.z;
 
-       sprintf(file,"%s-%.15f.xyz",v->fb,time);
+       sprintf(file,"%s-%07f.xyz",v->fb,time);
        fd=open(file,O_WRONLY|O_CREAT|O_TRUNC);
        if(fd<0) {
                perror("open visual save file fd");
@@ -93,7 +93,7 @@ int visual_atoms(t_visual *v,double time,t_atom *atom,int n) {
        dprintf(v->fd,"set ambient 20\n");
        dprintf(v->fd,"set specular on\n");
        dprintf(v->fd,"set boundbox on\n");
-       dprintf(v->fd,"label\n");
+       //dprintf(v->fd,"label\n");
        sprintf(file,"%s-%.15f.ppm",v->fb,time);
        dprintf(v->fd,"write ppm %s\n",file);
        dprintf(v->fd,"zap\n");
@@ -103,9 +103,9 @@ int visual_atoms(t_visual *v,double time,t_atom *atom,int n) {
        dprintf(fd,"atoms at time %.15f\n",time);
        for(i=0;i<n;i++)
                dprintf(fd,"%s %f %f %f\n",pse[atom[i].element],
-                                          10e9*atom[i].r.x,
-                                          10e9*atom[i].r.y,
-                                          10e9*atom[i].r.z);
+                                          atom[i].r.x,
+                                          atom[i].r.y,
+                                          atom[i].r.z);
        if(dim.x) {
                dprintf(fd,"%s %f %f %f\n",pse[0],dim.x/2,dim.y/2,dim.z/2);
                dprintf(fd,"%s %f %f %f\n",pse[0],-dim.x/2,dim.y/2,dim.z/2);