[완료] C언어로 코딩을 했는데.. 실행시 세그먼트 폴트 오류가 뜨는데 해결이 안됩니다..

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <ctype.h>
#include <math.h>
#include "sacdata.h"
 
long int read_sac(char *fname, SAC_DATA *sdata, float *beg, float *dt);
 
//SAC_DATA ew, ns;
SAC_DATA e1, n1, z1, e2, n2, z2;
 
float XCOR(long int n1, float *dat1, float *dat2);
float CRMS(float *dat, long int start, long int end);
 
int main(int argc, char *argv[])
{
	int i, j, k;
	long int a, b;
	float bege, begn, dt;
	float *coee, *coen, *coez, *coesum, *coeavg;
	float *coerms, *ti;
	FILE *f;
 
	if( argc < 7 )
	{
		fprintf(stderr,"\nUsage: %s e1.sac n1.sac z1.sac e2.sac n2.sac z2.sac\n",argv[0]);
		return -1;
	}
 
	read_sac(argv[1],&e1,&bege,&dt);
	read_sac(argv[2],&n1,&begn,&dt);
	read_sac(argv[3],&z1,&begn,&dt);
	read_sac(argv[4],&e2,&begn,&dt);
	read_sac(argv[5],&n2,&begn,&dt);
	read_sac(argv[6],&z2,&begn,&dt);
 
	if( e1.npts != n1.npts || e1.npts != z1.npts || n1.npts != z1.npts )
	{
		fprintf(stderr,"\nError:  e1, n1, z1 NPTS not same\n");
		return -1;
	}
 
	if( e2.npts != n2.npts || e2.npts != z2.npts || n2.npts != z2.npts )
	{
		fprintf(stderr,"\nError:  e2, n2, z2 NPTS not same\n");
		return -1;
	}
 
	/////////////// Phase info. ///////////////
 
	f=fopen("coe_data.txt", "w");
 
	for (i=0; i<e2.npts-e1.npts; i++)
	{
		coee[i] = XCOR(e1.npts, e1.fdat, e2.fdat+i);
		coen[i] = XCOR(n1.npts, n1.fdat, n2.fdat+i);
		coez[i] = XCOR(z1.npts, z1.fdat, z2.fdat+i);
		coesum[i] = coee[i] + coen[i] + coez[i];
		coeavg[i] = coesum[i] / 3.;
//		fprintf(f, "%d  %.2f  %.2f  %.2f  %.2f  %.2f  %.2f  %.2f\n", i, i*dt, coeavg, coee, coen, coez, coesum, coeavg);
	}
 
	a = 100;
	b = 250;
 
	for (i=0; i<e2.npts; i++)
	{
		coerms[i] = CRMS(coeavg, a+i, b+i);
		ti[i] = i*dt;
	}
 
	fprintf(f, "%.2f  %.2f  %.2f  %.2f  %.2f  %.2f  %.2f\n", ti, coee, coen, coez, coesum, coeavg, coerms);
	fclose(f);
 
	return 0;
}
 
float XCOR(long int n1, float *dat1, float *dat2)
{
	int i, j, k;
	float bege, begn, dt;
	float ewsum, nssum, ewavg, nsavg;
	float ewsq, nssq, ewd, nsd;
	float ewns, r;
 
	ewsum = 0.,   nssum = 0.,   ewavg = 0.,   nsavg = 0.;
	ewsq = 0.,   nssq = 0.,   ewd = 0.,   nsd = 0.;
	ewns = 0.,   r = 0.;
 
	////////// Average //////////
	// mean //
 
	for (i=0; i<n1; i++)
	{
		ewsum += dat1[i];
	}
 
	for (i=0; i<n1; i++)
	{
		nssum += dat2[i];
	}
 
	ewavg = ewsum / n1;
	nsavg = nssum / n1;
 
	////////// Standard Deviation ///////////
	// amplitude //
 
	for (i=0; i<n1; i++)
	{
		ewsq += (dat1[i] - ewavg) * (dat1[i] - ewavg);
	}
 
	for (i=0; i<n1; i++)
	{
		nssq += (dat2[i] - nsavg) * (dat2[i] - nsavg);
	}
 
	ewd = sqrt (ewsq / n1);
	nsd = sqrt (nssq / n1);
 
 
	///////// crosscorrelation coefficient /////////
 
	for (i=0; i<n1; i++)
	{
		ewns += ((dat1[i] - ewavg) / ewd) * ((dat2[i] - nsavg) / nsd);
//		ewns += dat1[i]*dat2[i];
	}
 
	r = ewns / (float)n1;
 
	return r;
}
 
float CRMS(float *dat, long int start, long int end)
{
	int i;
	float point;
	float n, m, l;
 
	point = 0.;
	n = 0., m = 0., l = 0.;
 
	point = (float)end - (float)start;
 
	for (i=start; i<end; i++)
	{
		n += dat[i] * dat[i];
	}
 
	for (i=-end; i<-start; i++)
	{
		m += dat[i] * dat[i];
	}
 
	l = sqrt((n + m) / (point * 2.));
 
	return l;
}