tau=TAU;
tau2=tau*tau;
+ /* rasmol script file */
+ sprintf(scr,"./saves/si.scr");
+ fd2=open(scr,O_WRONLY|O_CREAT|O_TRUNC);
+ if(fd2<0) {
+ perror("rasmol script file open");
+ return -1;
+ }
printf("starting velocity verlet: ");
fflush(stdout);
/* rasmol script & xyz file */
sprintf(xyz,"./saves/si-%.15f.xyz",time);
- sprintf(scr,"./saves/si-%.15f.scr",time);
sprintf(ppm,"./video/si-%.15f.ppm",time);
fd1=open(xyz,O_WRONLY|O_CREAT|O_TRUNC);
- fd2=open(scr,O_WRONLY|O_CREAT|O_TRUNC);
- if((fd1<0)||(fd2<0)) {
- perror("rasmol file open");
+ if(fd1<0) {
+ perror("rasmol xyz file open");
return -1;
}
dprintf(fd2,"load xyz %s\n",xyz);
dprintf(fd2,"rotate y 13\n");
dprintf(fd2,"set ambient 20\n");
dprintf(fd2,"set specular on\n");
- dprintf(fd2,"set write on\n");
dprintf(fd2,"write ppm %s\n",ppm);
+ dprintf(fd2,"zap\n");
dprintf(fd1,"%d\nsilicon\n",amount_si);
for(i=0;i<amount_si;i++)
dprintf(fd1,"Si %f %f %f %f\n",
si[i].x,si[i].y,si[i].z,time);
close(fd1);
- close(fd2);
}
}
printf(" done\n");
+ close(fd2);
free(si);
return 0;