#include <stdio.h>
#define MAX_LENGTH 256


main()
{
  int i;
  char error=0;
  char count;
  int num=0;
  char str[MAX_LENGTH];
  char ch;
  FILE *fpi, *fpo;
  double a, b, c, d, e;
  char *bgn = "# ./dust inputs -L=532 -n -ice -v -sm < table_data > depth1
set lnum 4
set mass 0 0 0 0.4
set mask 0 3 1 2
show\n";
  char *end = "quit\n";

  fpi=fopen("bufford.data", "r");
  fpo=fopen("table_data", "w");

  fprintf(fpo, "%s", bgn);

  while(1){
    for(i=0; i<MAX_LENGTH; i++){
      if(fscanf(fpi, "%c", &ch)!=1){
	error=1;
	break;
      }
      str[i]=ch;
      if(ch=='\n') break;
    }
    if(i==MAX_LENGTH) printf("Warning: out of buffer size\n");
    str[i]=0;
    num++;
    count=sscanf(str, "%lf %lf %lf %lf", &a, &b, &c, &d, &e);
    if(count!=4)
      printf("Line %i not included: number of lines %i must be 4\n", num, count);
    else fprintf(fpo, "%g\t%g\t%g\t%g\n", a, b, c, d);

    if(error) break;
  }

  fprintf(fpo, "%s", end);

  fclose(fpo);
  fclose(fpi);
}