- dim.x=10e9*v->dim.x;
- dim.y=10e9*v->dim.y;
- dim.z=10e9*v->dim.z;
-
- char pse[19][4]={
- "*",
- "H",
- "He",
- "Li",
- "Be",
- "B",
- "C",
- "N",
- "O",
- "F",
- "Ne",
- "Na",
- "Mg",
- "Al",
- "Si",
- "P",
- "S",
- "Cl",
- "Ar",
- };
-
- sprintf(file,"%s-%.15f.xyz",v->fb,time);
+ dim.x=v->dim.x;
+ dim.y=v->dim.y;
+ dim.z=v->dim.z;
+
+ 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");
fd=open(file,O_WRONLY|O_CREAT|O_TRUNC);
if(fd<0) {
perror("open visual save file fd");
dprintf(v->fd,"rotate x 100\n");
dprintf(v->fd,"rotate y 10\n");
dprintf(v->fd,"set ambient 20\n");
dprintf(v->fd,"set specular on\n");
dprintf(v->fd,"rotate x 100\n");
dprintf(v->fd,"rotate y 10\n");
dprintf(v->fd,"set ambient 20\n");
dprintf(v->fd,"set specular on\n");
sprintf(file,"%s-%.15f.ppm",v->fb,time);
dprintf(v->fd,"write ppm %s\n",file);
dprintf(v->fd,"zap\n");
sprintf(file,"%s-%.15f.ppm",v->fb,time);
dprintf(v->fd,"write ppm %s\n",file);
dprintf(v->fd,"zap\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],
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],
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);
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);