+int write_ac_distr(d3_lattice *d3_l,int ac_distr)
+{
+ int fd,x,y,z;
+ int count=0,offset;
+ char file[16];
+
+ if(ac_distr==1) strcpy(file,"a.plot");
+ if(ac_distr==2) strcpy(file,"c.plot");
+ if(ac_distr==3) strcpy(file,"b.plot");
+
+ if((fd=open(file,O_WRONLY|O_CREAT))<0)
+ {
+ puts("cannot open plot file");
+ return -1;
+ }
+
+ for(z=0;z<d3_l->max_z;z++)
+ {
+ for(x=0;x<d3_l->max_x;x++)
+ {
+ for(y=0;y<d3_l->max_y;y++)
+ {
+ offset=x+y*d3_l->max_x+z*d3_l->max_x*d3_l->max_y;
+ if(ac_distr==1)
+ if(*(d3_l->status+offset)&AMORPH) count+=*(d3_l->extra+offset);
+ if(ac_distr==2)
+ if(!(*(d3_l->status+offset)&AMORPH)) count+=*(d3_l->extra+offset);
+ if(ac_distr==3) count+=*(d3_l->extra+offset);
+ }
+ }
+ dprintf(fd,"%d %d\n",z,count);
+ }
+ close(fd);
+
+ return 1;
+}
+