/*
   Project: CP
   Jan Hamaekers
   1999

   File: CreateH2.c
   
*/
/*$Id: CreateH2.c,v 1.2 2007-09-03 16:28:58 heber Exp $*/

#include<stdlib.h>
#include<stdio.h> 
#include<math.h>

int main(int argc, char** argv) {
  int Steps, s;
  double start[3];
  double delta;
  double Delta[3], h2[3];
  FILE *target;
  char dummy[255];
  if(argc < 8) {
    fprintf(stderr,"Not enough arguments\nUsage: CreateH2 target Steps startx starty startz delta xyz\n");
    exit(1);
  }
  sprintf(dummy, "%s.h2", argv[1]);
  target = fopen(dummy, "w");
  if(target == 0) {
    fprintf(stderr,"Could not open targetbox\n");
    exit(1);
  }
  
  Steps = atoi(argv[2]);
  start[0] = atof(argv[3]);
  start[1] = atof(argv[4]);
  start[2] = atof(argv[5]);
  delta = atof(argv[6]);
  Delta[0] = 0.0; Delta[1] = 0.0; Delta[2] = 0.0;
  Delta[atoi(argv[7])] = delta/(Steps-1);
  h2[0] = 0.0; h2[1] = 0.0; h2[2] = 0.0;
  h2[atoi(argv[8])] = 0.37042405;
  fprintf(target, "2 1 1.0 1 1 0.1\n");
  fprintf(target, "%e\t%e\t%e\t2\t%i\n", start[0]-h2[0], start[1]-h2[1], start[2]-h2[2], 2*Steps-2);
  for (s=1; s<Steps; s++) {
    fprintf(target, "%e\t%e\t%e\n", start[0]+s*Delta[0]-h2[0], start[1]+s*Delta[1]-h2[1], start[2]+s*Delta[2]-h2[2]);
  }
  for (s=Steps-2; s>=0; s--) {
    fprintf(target, "%e\t%e\t%e\n", start[0]+s*Delta[0]-h2[0], start[1]+s*Delta[1]-h2[1], start[2]+s*Delta[2]-h2[2]);
  }
  fprintf(target, "%e\t%e\t%e\t2\t%i\n", start[0]+h2[0], start[1]+h2[1], start[2]+h2[2], 2*Steps-2);
  for (s=1; s<Steps; s++) {
    fprintf(target, "%e\t%e\t%e\n", start[0]+s*Delta[0]+h2[0], start[1]+s*Delta[1]+h2[1], start[2]+s*Delta[2]+h2[2]);
  }
  for (s=Steps-2; s>=0; s--) {
    fprintf(target, "%e\t%e\t%e\n", start[0]+s*Delta[0]+h2[0], start[1]+s*Delta[1]+h2[1], start[2]+s*Delta[2]+h2[2]);
  }
  fclose(target);
  return(0);
}