GRACE_Matlab_Toolbox使用简介

（1）本测试是在matlab2015b环境下测试完成的，软件的下载地址https://github.com/fengweiigg/GRACE_Matlab_Toolbox；

（2）软件的第一步其实是控制文件（GMT_Control_File）脚本的书写；

控制文件截图

（3）因为这里处理的是RL06产品，而代码处理的是RL05数据，所以这里有几处是需要修改的：二阶项C20采用最新的C20_RL06，需要把程序包里面gmt_replace_C20函数中RL05替换成RL06才可以运行；c21和s21替换文件中说明文件删除，然后将gmt_replace_C21_S21_C22_S22函数中get_days改为gmt_get_days。这里顺便多说一句，数据文件是从ICGEM下载的最新的，不然可能出问题。

（4）点击运行直接出结果，这里生成的结果大概4.5M。

生成结果

（5）接下来进行泄漏的改正，这个时间可能稍微长一些，要耐心等待，不要关闭软件。

（6）接下来是将生成的球谐分析，将生成的球谐数据文件转换为grid_data.mat：

（7）将上一步生成的格网数据转换为时间序列Grid2Series，以亚马逊区域为例：

顺便将生成的文件在matlab成图：

load time_series.mat

plot(time,time_series*100,'-s');

xlabel('Year');

ylabel('Equivalent water height (cm)');

title('Mass variations in the Amazon from GRACE');

亚马逊流域时间序列变化

（8）继续分析，我们可以得到全球趋势项变化。