+ 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);