#include?
#include?
//#define?M?179

#define?TMLOGPATH?????“E:\\DTF5\\CE318data_software\\x081111.txt“
#define?V0RSDPATH?????“E:\\DTF5\\CE318data_software\\V0_R_SD.txt“
#define?DEPTHPATH?????“E:\\DTF5\\CE318data_software\\AODx081111.txt“
#define?ALLPATH???????“E:\\DTF5\\CE318data_software\\allx081111.txt“

#define?ROW?17
void?main（）
	{
????	int?ijM=0;
		FILE?*p;
		char?chstr[10]=“x081111“;
	????double?timal[ROW]m[ROW]V[ROW][5]mavga[5]b[5]lat=31.9037lon=117.1584;//大營盤lat=40.255lon=115.8;//合肥lat=31.9037lon=117.1584;風電場lat=41.017lon=114.758//經緯度注意是否要修改
		double?Vavg[5]V0[5]total=0.0sumV[5]dx=0.0dy[5]dxy[5]R[5]SD[5];
		double?bm[ROW]bmavgVm[ROW]Vmavg;
		double?tao[ROW][5]taor[5]taoa8[ROW][5]taoa940[ROW]tao940[ROW]taoH2O[ROW]dxy940=0.0ab=0.606ba=0.580;//Wb[M]W[M]
		double?taoa5[ROW][4]lnnamta5[4]namta5avgsumtaoa5[ROW]taoa5avg[ROW]dxym[ROW]slope1[ROW]arfa[ROW]beta[ROW];
		double?namta8[5]={1.0200.9360.8700.6700.440};
			
		for（i=0;i<5;i++）
		{
		????????dxy[i]=0.0;
				dy[i]=0.0;
		????????sumV[i]=0.0;
				SD[i]=0.0;
		}

		for（i=0;i		{
			dxym[i]=0.0;
			sumtaoa5[i]=0.0;
		}


		if（（p=fopen（TMLOGPATH“r“））==NULL）
		{
			printf（“can‘t?open?the?file!\n“）;
			return;
		}
		ch=fgetc（p）;
		while（ch!=EOF）
		{
			if（ch==‘\n‘）?M++;
			ch=fgetc（p）;
		}
		fclose（p）;

		printf（“%d\n“M）;


		if（（p=fopen（TMLOGPATH“r“））==NULL）
		{
			printf（“can‘t?open?the?file!\n“）;
			return;
		}
		for（i=0;i		{
			fscanf（p“%lf?%lf“&timal[i]&m[i]）;
			for（j=0;j<5;j++）
				fscanf（p“%lf“&V[i][j]）;
		}
		fclose（p）;

				
	//計算大氣層外儀器測量值V0及定標因子C


	for（i=0;i		total=total+m[i];
	mavg=total/M;
	//printf（“%lf\n“mavg）;
	for（i=0;i		for（j=0;j<5;j++）
	????	sumV[j]=sumV[j]+V[i][j];
	for（j=0;j<5;j++）
	{
		Vavg[j]=sumV[j]/M;
		//printf（“%lf\t“exp（Vavg[i]））;
	}
	for（i=0;i?????????dx=dx+pow（（m[i]-mavg）2）;
????for（i=0;i<5;i++）
		for（j=0;j		{
			dy[i]=dy[i]+pow（（V[j][i]-Vavg[i]）2）;
			dxy[i]=dxy[i]+（m[j]-mavg）*（V[j][i]-Vavg[i]）;
		}
	for（i=0;i<5;i++）
		{
			a[i]=1.0*dxy[i]/dx;//斜率
			b[i]=Vavg[i]-a[i]*mavg;//截距
			V0[i]=exp（b[i]）;
			R[i]=dxy[i]/sqrt（dx*dy[i]）;
		}
	for（i=0;i<5;i++）
		for（j=0;j			SD[i]=SD[i]+pow（（V[j][i]-a[i]*m[j]-b[i]）2）;
	for（i=0;i<5;i++）
		SD[i]=sqrt（SD[i]/M）;
	for（i=0;i<5;i++）
		printf（“%lf\t%lf\t%lf\n“exp（b[i]）R[i]SD[i]）;

	
	//水汽波段定標	
????//計算大氣柱總光學厚度
	for（i=0;i		for（j=0;j<5;j++）?
			V[i][j]=exp（V[i][j]）;
	for（i=0;i		for（j=0;j<5;j++）?
		????tao[i][j]=log（V0[j]/V[i][j]）/m[i];
		
	for（i=0;i<5;i++）
	{
		taor[i]=0.008569*pow（namta8[i]-4）*（1+0.0113*pow（namta8[i]-2）+0.00013*pow（namta8[i]-4））*1000*exp（-0.125017*0.044）/1013.25;??//計算大氣分子垂直光學厚度taor[i]=0.0088*pow（namta8[i]-4.05）
		//printf（“%lf\t“taor[i]）;
	}
		
	for（i=0;i		for（j=0;j<5;j++）
			taoa8[i][j]=tao[i][j]-taor[j];??//計算氣溶膠光學厚度，含水汽
	
		//利用線性內

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-06-13?14:33??sun_photometer_ce318_data_deal\
?????文件????????2602??2012-12-27?15:06??sun_photometer_ce318_data_deal\AOD_CMP.CPP
?????文件????????6540??2012-12-27?15:10??sun_photometer_ce318_data_deal\dingbiao.cpp
?????文件????????3701??2012-12-27?15:09??sun_photometer_ce318_data_deal\Fit_Ang.cpp
?????文件????????3742??2018-06-13?14:31??sun_photometer_ce318_data_deal\read_ce318_data.cpp
?????文件????????2723??2012-12-27?15:05??sun_photometer_ce318_data_deal\WAT_CMP.CPP