/*?============================================================================?*/
/*???The?Calculation?of?the?Global?Equivalent?water?layer?thickness?from?GRACE
?????Time?Variable?Gravity?Field????????????????????????????????????????????????*/
/*???EWLT-v1.0.cpp??????????????????????????????????????????????????????????????*/
/*???（c）?2010?DU?JinsongCUG????????????????????????????????????????????????????*/
/*???E-mail:?jinsongdu.cug@gmail.com????????????????????????????????????????????*/
/*???Version?1.0????????????????????????????????????????????????????????????????*/
/*???Calculation?from?the?spherical?harmonic?coefficiences?（No?smoothing）???????*/
/*?============================================================================?*/

#include?
#include?
#include?

void?Legendre（long?Nmax?long?m?double?theta?double*?prev）;
inline?double?L_an（long?m?long?n）;
inline?double?L_cn（long?m?long?n）;

/*------------------------------------------------------------------------------*/

int?main（）
{
????const?double?PI?=?3.1415926535897932;

????//?輸入球諧系數(shù)文件名和保存結(jié)果的文件名
????char?scIn[256];
????char?sgOut[256];

????printf（“\nPlease?input?the?coefficents?file?:??“）;
????scanf（“%s“?scIn）;

????printf（“\nPlease?input?a?file?name?to?save?the?results?:??“）;
????scanf（“%s“?sgOut）;

????FILE*?fpIn=fopen（scIn?“rt“）;
????if（fpIn==NULL）
????{
????????printf（“\nSorry!\nFailed?to?open?the?file:?“）;
????????printf（“%s“scIn）;
????????printf（?“.\n“）;
????????return?-1;
????}

????//?獲取初始計(jì)算的參數(shù)
????long????Nmax;????????//?最高階數(shù)
????long????N_start;?????//?起始階數(shù)
	long????N_end;???????//?截止階次
????double??Altitude;????//?高程（m）
????double??R;???????????//?球半徑（m）
????double??GM;		?????//?引力常數(shù)（m3/s2）
????double??dd;			?//?網(wǎng)格間距（degree）
	double??Ws;??????????//?起始緯度（degree）
	double??We;??????????//?終止緯度（degree）
	double??Js;??????????//?起始經(jīng)度（degree）
	double??Je;??????????//?終止經(jīng)度（degree）
	double??density_E;???//?地球平均密度（kg/m3）
	double??density_W;???//?水的密度（kg/m3）

	Nmax=60;
	N_start=2;
	N_end=60;
	Altitude=0;
	R=0.6378136300E+07;
	GM=0.3986004415E+15;
	dd=1;
	Ws=-90;
	We=90;
	Js=-180;
	Je=180;
	density_E=5517;
	density_W=1000;

	printf（“\n“）;
????printf（“???????????????????Max?Dregree?N?=?%ld\n“?Nmax）;
????printf（“???????The?Start?Dregree?N_start?=?%ld\n“?N_start）;
	printf（“???????????The?End?Dregree?N_end?=?%ld\n“?N_end）;
????printf（“????????????????????????Altitude?=?%.3lf?m\n“?Altitude）;
????printf（“???????????????????????????????R?=?%.3lf?m\n“?R）;
????printf（“??????????????????????????????GM?=?%.3lf?m3/s2\n“?GM）;
????printf（“????????????????????????Gridstep?=?%.3lf?degree\n“?dd）;
????printf（“??????????????The?start?latitude?=?%.3lf?degree\n“?Ws）;
????printf（“????????????????The?end?latitude?=?%.3lf?degree\n“?We）;
????printf（“?????????????The?start?longitude?=?%.3lf?degree\n“?Js）;
????printf（“???????????????The?end?longitude?=?%.3lf?degree\n“?Je）;
????printf（“T

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????93921??2010-12-04?16:35??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\200304-200301.txt

?????文件???????7645??2012-05-15?22:29??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\EWLT-v1.0.cpp

?????文件???????8686??2012-05-15?22:29??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\EWLT-v2.0.cpp

?????文件???????8594??2012-05-15?22:29??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\EWLT-v2.1.cpp

?????文件???????9171??2012-05-15?22:30??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\EWLT-v3.0.cpp

?????文件???????8665??2012-05-15?22:30??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\EWLT-v4.0.cpp

?????文件??????13900??2012-05-15?22:30??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\EWLT-v6.0.cpp

?????文件??????14355??2012-05-15?22:30??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\EWLT-v6.1.cpp

?????文件??????15602??2012-05-15?22:30??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\EWLT-v6.2.cpp

?????文件??????15250??2012-05-15?22:30??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\EWLT-v6.3.cpp

?????文件??????14453??2012-05-15?22:31??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\EWLT-v6.4.cpp

?????文件??????15412??2012-05-15?22:31??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\EWLT-v6.5.cpp

?????文件???????2551??2009-11-14?09:47??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）\LoveNum.dat

?????目錄??????????0??2012-05-15?22:31??GRACE水儲(chǔ)量解算與濾波程序代碼-杜勁松（20120515）

-----------?---------??----------?-----??----

???????????????228205????????????????????14