/* Multichannel average*/

#include "stdio.h"
main () {
  int data[1025][30], dout[1025][30]; /*(*p)[30];*/
  int i, j, k;
  FILE *fp1, *fp2, *fp3; 
  static char rec_file[]="record01.dat";

  fp1=fopen ("new_all.asc","r");
  for (i=0; i<1025; i++) {
     for (j=0; j<29; j++) {
	fscanf (fp1, "%d ", &data[i][j]);  
     }
     fscanf (fp1, "%d", &data[i][29]);
  }
  fclose (fp1);

  for (j=1; j<30; j++) {
     for (i=0; i<1025; i++) {
	data[i][j]=data[i][j-1]+data[i][j];
     }
  }

  for (j=1; j<30; j++) {
     for (i=0; i<1025; i++) {
	data[i][j]=(float)(data[i][j])/(j+1)+0.5;
     }
  }

  fp2=fopen ("new_all.average", "w");
  for (i=0; i<1025; i++) {
     for (j=0; j<29; j++) {
        fprintf (fp2, "%d ", data[i][j]);
     }
     fprintf (fp2, "%d\n", data[i][29]);
  }
  fclose (fp2);

  for (i=0; i<30; i++) {
     if (i<10) {
       j=i+48; 
       rec_file[7]=j;
     }
     if (i>9) {
       rec_file[6]='1';
        j=i-10+48;
       rec_file[7]=j;
     }
     if (i>19) {
       rec_file[6]='2';
        j=i-20+48;
       rec_file[7]=j;
     }
     fp2=fopen (rec_file,"w");
     for (j=0; j<1025; j++)
	fprintf (fp2, "%d   %d\n", j, data[j][i]);
     fclose (fp2);
  }
}