#include <errno.h>
	#include <stdio.h>
	#include <stdlib.h>
	#include "atoms.h"
	#include "header.h"
	int main(int argc, char** argv)
	{
	Headerinfo tempheader;
	Atom tempatom;
	FILE* file1;
	FILE* file2;
	int i,j,k;
	char buffer[255];
	double sizemul;
	double cutoff;
	double curX, curY, curZ;
	if(argc < 2)
	{
	printf("no infile specified\n\n");
	return 0;
	}
	for(k=1; k<argc; k++)
	{
	file1 = fopen(argv[k], "rb");
	if(file1 == NULL)
	{
	printf("failed to open infile %s for reading\n\n", argv[k]);
	return errno;
	}
	printf("using infile %s\n\n", argv[k]);
	j = (int)fread(&tempheader, sizeof(Headerinfo), 1, file1);
	if(j < 1)
	{
	printf("error reading header from infile\n\n");
	return 0;
	}
	sprintf(buffer, "%s.pov", argv[k]);
	file2 = fopen(buffer, "w");
	sizemul = max(tempheader.bigX, max(tempheader.bigY, tempheader.bigZ));
	if(sizemul <= 0.0)
	{
	fprintf(file2, "crash HARD and die\n");
	fclose(file2);
	fclose(file1);
	return 5;
	}
	sizemul = 1/sizemul;
	cutoff = tempheader.mincutoff * sizemul;
	curX = tempheader.bigX*sizemul;
	curY = tempheader.bigY*sizemul;
	curZ = tempheader.bigZ*sizemul;
	fprintf(file2, "#include \"colors.inc\"\n\n");
	fprintf(file2, "camera {\n\tlocation <%f, %f, %f>\n\tlook_at <%f, %f, %f>\n}\n\n", 0.1+0.5*curX, 0.7+0.5*curY, -(0.8+0.5*curZ), 0.5*curX, 0.5*curY, 0.5*curZ);
	fprintf(file2, "light_source { <10, 10, -10>, White }\n\n");
	fprintf(file2, "background { color Black }\n\n");
	fprintf(file2, "cylinder { <0.0, -0.025, -0.025>, <0.0, -0.025, -0.025> + <%f, 0.0, 0.0>, 0.01 pigment { Green } }\n\n", 1.0E-9*sizemul);
	fprintf(file2, "cylinder { <0.0, -0.025, -0.025>, <0.0, -0.025, -0.025> + <%f, 0.0, 0.0>, 0.011 pigment { GreenCopper } }\n\n", 1.0E-10*sizemul);
	fprintf(file2, "cylinder { <0.0, 0.0, 0.0>, <%f, 0.0, 0.0>, 0.002 pigment { Red } }\n", curX);
	fprintf(file2, "cylinder { <0.0, 0.0, 0.0>, <0.0, %f, 0.0>, 0.002 pigment { Red } }\n", curY);
	fprintf(file2, "cylinder { <0.0, 0.0, 0.0>, <0.0, 0.0, %f>, 0.002 pigment { Red } }\n", curZ);
	fprintf(file2, "cylinder { <%f, 0.0, 0.0>, <%f, 0.0, %f>, 0.002 pigment { Red } }\n", curX, curX, curZ);
	fprintf(file2, "cylinder { <%f, 0.0, 0.0>, <%f, %f, 0.0>, 0.002 pigment { Red } }\n", curX, curX, curY);
	fprintf(file2, "cylinder { <0.0, %f, 0.0>, <0.0, %f, %f>, 0.002 pigment { Red } }\n", curY, curY, curZ);
	fprintf(file2, "cylinder { <0.0, %f, 0.0>, <%f, %f, 0.0>, 0.002 pigment { Red } }\n", curY, curX, curY);
	fprintf(file2, "cylinder { <0.0, 0.0, %f>, <%f, 0.0, %f>, 0.002 pigment { Red } }\n", curZ, curX, curZ);
	fprintf(file2, "cylinder { <0.0, 0.0, %f>, <0.0, %f, %f>, 0.002 pigment { Red } }\n", curZ, curY, curZ);
	fprintf(file2, "cylinder { <%f, %f, %f>, <0.0, %f, %f>, 0.002 pigment { Red } }\n", curX, curY, curZ, curY, curZ);
	fprintf(file2, "cylinder { <%f, %f, %f>, <%f, 0.0, %f>, 0.002 pigment { Red } }\n", curX, curY, curZ, curX, curZ);
	fprintf(file2, "cylinder { <%f, %f, %f>, <%f, %f, 0.0>, 0.002 pigment { Red } }\n", curX, curY, curZ, curX, curY);
	fprintf(file2, "\n\n");
	for(i=0; i< (sizeof(elements)/sizeof(elements[0])); i++)
	{
	fprintf(file2, "#macro %s(center)\n", elements[i].symbol);
	fprintf(file2, "\tsphere { center, %f pigment { %s } }\n", elements[i].radius*sizemul, elements[i].pigment);
	fprintf(file2, "#end\n");
	}
	fprintf(file2, "\n");
	for(i=0; i<tempheader.numatoms; i++)
	{
	j = (int)fread(&tempatom, sizeof(Atom), 1, file1);
	if(j < 1)
	{
	printf("error reading atom %d from infile\n\n", i);
	return 3;
	}
	curX = sizemul * tempatom.x;
	curY = sizemul * tempatom.y;
	curZ = sizemul * tempatom.z;
	fprintf(file2, "%s( <%f,%f,%f> )\n", elements[tempatom.type].symbol, curX, curY, curZ);
	}
	fclose(file2);
	fclose(file1);
	}
	return 1;
	}