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

int main(int argc, char *argv[])
{
  int i, k, len, shadowstart;
  double *x, *energy, *s2, *s4, *s6, *s8;
  double base;
  char buf[64];
  FILE *fp;

  if(argc!=3)
    return 1;

  len = atoi(argv[2]);
  fp = fopen(argv[1],"r");

  x = (double *)calloc(len,sizeof(double));
  energy = (double *)calloc(len,sizeof(double));
  s2 = (double *)calloc(len,sizeof(double));
  s4 = (double *)calloc(len,sizeof(double));
  s6 = (double *)calloc(len,sizeof(double));
  s8 = (double *)calloc(len,sizeof(double));

  for(i=0; i<len; i++)
    fscanf(fp,"%lg %lg %lg %lg %lg %lg",
	   &x[i],&energy[i],&s2[i],&s4[i],&s6[i],&s8[i]);

  fclose(fp);

  shadowstart=0;
  while(s2[++shadowstart]==0);

  base = energy[0];
  for(i=0; i<len; i++)
    energy[i] = energy[i]-base;

  /*base = s2[shadowstart]; */
  for(i=shadowstart; i<len; i++)
    s2[i] = s2[i]-base;

  /*base = s4[shadowstart]; */
  for(i=shadowstart; i<len; i++)
    s4[i] = s4[i]-base;
    
  /*base = s6[shadowstart]; */
  for(i=shadowstart; i<len; i++)
    s6[i] = s6[i]-base;
    
  /*base = s8[shadowstart]; */
  for(i=shadowstart; i<len; i++)
    s8[i] = s8[i]-base;

  strcpy(buf,argv[1]);
  strcat(buf,".err");

  fp = fopen(buf,"w");
  for(i=0; i<len; i++)
    fprintf(fp,"%d %16.10g %16.10g %16.10g %16.10g %16.10g %d\n",
	    (int)x[i], energy[i], s2[i], s4[i], s6[i], s8[i], 0);
  fclose(fp);
  
  return 0;
}