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

#ifndef N_AVE
#define N_AVE 5
#endif

void diff(char* s,int n){
    FILE *fp=stdout;
    FILE *f[N_AVE];
    char nomefile[1000];
    int i,j,k,Np,pert,totweight,weight[N_AVE+1];
    double corx[N_AVE+1],cory[N_AVE+1],corxy[N_AVE+1],m3[N_AVE+1],m4[N_AVE+1],t,T,fi;


    for(i=0;i<N_AVE;i++){
      sprintf(nomefile,"%s%d.dat",s,i);
      f[i]=fopen(nomefile,"r");
      if (f[i]==NULL){
	  fprintf(stderr,"error file in input: '%s' does not exist\n",nomefile);
	  exit(0);
      }


    }
    /*****  HEADER******/
    totweight=0;
    for (i=0;i<N_AVE;i++){
	fscanf(f[i],"#Number of realizations: %d\n#Tracer=100\n",&weight[i]);
	totweight+=weight[i];
  }
    /****END HEADER*****/

  	    fprintf(fp,"#Number of realizations: %d\n#Tracer=100\n",totweight);

    for (j=0;j<n;j++){
      for (i=0;i<N_AVE;i++)
	  if  (fscanf(f[i],"%lf %lf %lf %lf",&t,&corx[i],&cory[i],&corxy[i])==EOF) 
	  {
	      exit(0);
	      fprintf(stderr,"EOF error!!!\n input less lines exit!");
	  }


	corx[N_AVE]=0.;
	cory[N_AVE]=0.;
	corxy[N_AVE]=0.;

	    for(k=0;k<N_AVE;k++){
		corx[N_AVE]+=corx[k]*weight[k];
		cory[N_AVE]+=cory[k]*weight[k];
		corxy[N_AVE]+=corxy[k]*weight[k];
	    }
	    fprintf(fp,"%lf %lf %lf %lf\n",t,corx[N_AVE]/totweight,cory[N_AVE]/totweight,corxy[N_AVE]/totweight);

    }
    return;
}

int main(int argc, char **argv){
    if (argc<1){
	printf("insert string and numberoflines\n");
	return 0;
    }
    diff(argv[1],atoi(argv[2]));
    return 0;
}