+ dprintf(fd1,"H 0.0 0.0 0.0 %f\n",time);
+ dprintf(fd1,"He %f %f %f %f\n",LX,LY,LZ,time);
+ dprintf(fd1,"He %f %f %f %f\n",-LX,LY,LZ,time);
+ dprintf(fd1,"He %f %f %f %f\n",LX,-LY,LZ,time);
+ dprintf(fd1,"He %f %f %f %f\n",LX,LY,-LZ,time);
+ dprintf(fd1,"He %f %f %f %f\n",-LX,-LY,LZ,time);
+ dprintf(fd1,"He %f %f %f %f\n",-LX,LY,-LZ,time);
+ dprintf(fd1,"He %f %f %f %f\n",LX,-LY,-LZ,time);
+ dprintf(fd1,"He %f %f %f %f\n",-LX,-LY,-LZ,time);
+ close(fd1);
+
+ }
+
+ /* increase time */
+ time+=tau;