GRACE水储量解算与滤波程序代码-杜劲松

GRACE水储量解算与滤波程序代码，利用重力卫星反演陆地水储量及地下水变化的可以看看

/* ============================================================================ */
/*   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;

    // 输入球谐系数文件名和保存结果的文件名
    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;
    }

    // 获取初始计算的参数
    long    Nmax;        // 最高阶数
    long    N_start;     // 起始阶数
	long    N_end;       // 截止阶次
    double  Altitude;    // 高程（m）
    double  R;           // 球半径（m）
    double  GM;		     // 引力常数（m3/s2）
    double  dd;			 // 网格间距（degree）
	double  Ws;          // 起始纬度（degree）
	double  We;          // 终止纬度（degree）
	double  Js;          // 起始经度（degree）
	double  Je;          // 终止经度（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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      93921  2010-12-04 16:35  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\200304-200301.txt

     文件       7645  2012-05-15 22:29  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\EWLT-v1.0.cpp

     文件       8686  2012-05-15 22:29  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\EWLT-v2.0.cpp

     文件       8594  2012-05-15 22:29  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\EWLT-v2.1.cpp

     文件       9171  2012-05-15 22:30  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\EWLT-v3.0.cpp

     文件       8665  2012-05-15 22:30  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\EWLT-v4.0.cpp

     文件      13900  2012-05-15 22:30  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\EWLT-v6.0.cpp

     文件      14355  2012-05-15 22:30  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\EWLT-v6.1.cpp

     文件      15602  2012-05-15 22:30  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\EWLT-v6.2.cpp

     文件      15250  2012-05-15 22:30  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\EWLT-v6.3.cpp

     文件      14453  2012-05-15 22:31  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\EWLT-v6.4.cpp

     文件      15412  2012-05-15 22:31  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\EWLT-v6.5.cpp

     文件       2551  2009-11-14 09:47  GRACE水储量解算与滤波程序代码-杜劲松（20120515）\LoveNum.dat

     目录          0  2012-05-15 22:31  GRACE水储量解算与滤波程序代码-杜劲松（20120515）

----------- ---------  ---------- -----  ----

               228205                    14