/* This program, written in ANSI C (NO GUARANTEES AS TO PORTABILITY! IT RUNS ON
 * SunOS 4.1.1, COMPILED WITH GCC), does not require HTTPD, CGI, or Mosaic.  It
 * reads either the standard input or a series of files in the directory given by
 * the -d argument whose unit number is given, and prints plain text on stdout.
 *
 * The unit descriptions must be named box%.unit#, where % represents a
 * box number, of arbitrary format, and # is the unit
 * number, unpadded.  Clearly this must be changed for DOS.
 *
 * The command line syntax is as follows:
 *	lithterp_text [-ddirectory] [unit number]
 * With no arguments, the stdin is read until EOF or EOT.  Unit number must be a
 * positive integer. If -d is absent and a number is given, the current directory (.)
 * is used. Please, no space after -d!
 *
 * 6/17/94 Version 2: RQD fields moved to header, fields renumbered
 * 6/16/94 Version 1: modified from cgi interpreter at expet
 * Paul Asimow (paul@expet.gps.caltech.edu)
*/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <fcntl.h>
#include <malloc.h>

void main(int argc, char *argv[]) {
	FILE *fp1;
	char *buf, *data[100], instring[600], outstring[100];
	char filename[20][50], *directory;
	int field,unit,i=0,j,k;
	int stdinflag=0, activecomma;
	buf = (char *) malloc(sizeof(char));
	
	if (argc == 1) {
		stdinflag = 1;
		i = 1;
	} else {
		if (*argv[1] != '-') {
			if ((unit = atoi(argv[1])) == 0) {
				printf("Bad unit number\n"); exit(1);
			}
			directory = (char *) malloc(2*sizeof(char));
			strcpy(directory,".");
		} else if (*(++argv[1]) != 'd') {
			printf("Syntax: lithterp_text [-d directory][unit number]\n");
			exit(1);
		} else {
			directory=(char *) malloc(strlen(++argv[1])*sizeof(char));	
			strcpy(directory,argv[1]);
			if ((unit = atoi(argv[2])) == 0) {
				printf("Bad unit number\n"); exit(1);
			}
		}
		/* Search directory for files with right unit number */
		sprintf(outstring, "ls %s/box*.unit%d",directory,unit);
		fp1 = popen(outstring, "r");
		while (fscanf(fp1, "%s", filename[i]) != EOF) i++;	
		pclose(fp1);
		free(directory);

		printf("Lithologic Unit #%d (found in %d boxes)\n\n",unit,i);
	}					

	for (j=0;j<i;j++) {	/* Loop over boxes in which the unit was found */
		field = 0; activecomma = 1; strcpy(instring,"");
		if (!stdinflag) fp1 = fopen(filename[j],"r");
		else fp1 = stdin;
		while ((*buf = getc(fp1)) != EOF) {	/* read one byte at a time */
			if (*buf == 34) activecomma = 1-activecomma; /*Double-quotes?*/
			else if (activecomma == 1 && *buf == 44) {   /*End of field?*/
				data[field] = (char *) malloc((1+strlen(instring))*sizeof(char));
				strcpy(data[field], instring);	     /* Put data in*/
				field++;			     /* array      */
				strcpy(instring,"");
			} else strncat(instring,buf,1); /*else add character to string */
		}	
		if (strlen(instring)) {		/*In case file doesn't end with comma*/
			data[field] = (char *) malloc((1+strlen(instring))*sizeof(char));
			strcpy(data[field],instring);
			field++;
		}
		if (!stdinflag) fclose(fp1);
		printf("-----------------------\nBox %s\n", data[0]);	
		printf("Core runs in box: ");
		for(k=1;k<=6;k++) {
			if (data[k] == "") break;
			else printf("%s ",data[k]);
		}
		printf("\nRecovery in %% core runs in box: ");
		for(k=7;k<=12;k++) {
			if (data[k] == "") break;
			else printf("%s ", data[k]);
		}
		printf("\nRQD in ratio form for core runs: ");
                for(k=13;k<=18;k++) {
                        if (data[k] == "") break;
                        else printf("%s ",data[k]);
                }

		printf("\nLogged by %s, on date %s\n", data[19], data[20]);
		printf("Checked by %s on date %s\n", data[21], data[22]);	
		printf("Driller's depth of box: top %s ft., bottom %s ft.\n", data[23], data[24]);
		printf("Core Type: %s; # of Units in box: %s\n",data[25],data[26]);
	
		printf("-----------------------\nUnit #%s, %s\n",data[27],data[28]); 
		if (strlen(data[29])) printf("Unit type: %s\n",data[29]);
		if (strlen(data[30])) printf("Comments: %s\n",data[30]);
		printf(" * Contacts\n%cTop:",9);
		if (strlen(data[31])) printf("Core run %s,",data[31]);
		if (strlen(data[32])) printf(" revised depth %s,",data[32]);
		if (strlen(data[33])) printf(" contact type %s",data[33]);
		printf("\n%cBottom:",9);
                if (strlen(data[34])) printf("Core run %s,",data[34]);
                if (strlen(data[35])) printf(" revised depth %s,",data[35]);
                if (strlen(data[36])) printf(" contact type %s",data[36]);
		if (strlen(data[37])) printf("\n%cComments: %s",9,data[37]);
		printf("\n\n * Phenocrysts\n");
		if (strlen(data[38])) printf("%cOverall abundance %s\n",9,data[38]);
		if (strlen(data[39])) printf("%cOverall comments: %s\n",9,data[39]);
		if (!strlen(data[38]) && !strlen(data[39])) printf("%cnone\n",9);
		if (strlen(data[40])) {	/*Only continue if fields non-empty*/
			printf("%ca: %s.",9,data[40]);
			if (strlen(data[41])) printf(" Abundance %s,",data[41]);
			if (strlen(data[42])) printf(" Average size %s,",data[42]);
			if (strlen(data[43])) printf(" Shape %s,",data[43]);
			if (strlen(data[44])) printf(" Alteration %s,",data[44]);
			if (strlen(data[45])) printf(" Comments: %s\n",data[45]);
			if (strlen(data[46])) {
				printf("%cb: %s.",9,data[46]);
                	        if (strlen(data[47])) printf(" Abundance %s,",data[47]);
                       		if (strlen(data[48])) printf(" Average size %s,",data[48]);
                        	if (strlen(data[49])) printf(" Shape %s,",data[49]);
                        	if (strlen(data[50])) printf(" Alteration %s,",data[50]);
                        	if (strlen(data[51])) printf(" Comments:%s/n",data[51]);
				if (strlen(data[52])) {
					printf("%cc: %s.",9,data[52]);
                        		if (strlen(data[53])) printf(" Abundance %s,",data[53]);
                        		if (strlen(data[54])) printf(" Average size %s,",data[54]);
                        		if (strlen(data[55])) printf(" Shape %s,",data[55]);
                        		if (strlen(data[56])) printf(" Alteration %s,",data[56]);
                        		if (strlen(data[57])) printf(" Comments: %s\n",data[57]);
				}
			}
		}
		if (strlen(data[59])||strlen(data[60])||strlen(data[58])) {
			printf("\n * Groundmass\n%c",9);
			if (strlen(data[54])) printf(" Lithology %s,",data[58]);
			if (strlen(data[55])) printf(" Texture %s,",data[59]);
			if (strlen(data[56])) printf(" Comments %s\n",data[60]);
		}
		if (strlen(data[61])) printf("\n * Color (dry): %s",data[61]);	
		if (strlen(data[62])) printf("\n * Sedimentary/Volcanic Structures: %s", data[62]);
		if (strlen(data[63])) printf("\n * Sorting: %s",data[63]);
		if (strlen(data[64])) {
			printf("\n * Vesicles\n%cabundance %s,",9,data[64]);
			if (strlen(data[65])) printf(" Average size %s,",data[65]);
			if (strlen(data[66])) printf(" Shape %s,",data[66]);
			if (strlen(data[67])) printf(" Fillings %s,",data[67]);
			if (strlen(data[68])) printf(" Aspect ratio %s,",data[68]);
			if (strlen(data[69])) printf(" Comments %s\n",data[69]);
		}
		if (strlen(data[70])||strlen(data[71])||strlen(data[72])) {
			printf("\n * Alteration\n%c",9);
			if (strlen(data[66])) printf(" Abundance %s,",data[70]);
			if (strlen(data[68])) printf(" Mineralogy %s,",data[71]);
			if (strlen(data[67])) printf(" Comments %s\n",data[72]);
		}
		if (strlen(data[73])) printf("\n * Veins: %s",data[73]);
		if (strlen(data[74])) printf("\n * Fractures: %s",data[74]); 
		if (strlen(data[75])) printf("\n * Additional Comments: %s",data[75]);
		printf("\n");
		for (k=0;k<field;k++) free(data[k]);
	}
}