（1）本測(cè)試是在matlab2015b環(huán)境下測(cè)試完成的判呕，軟件的下載地址https://github.com/fengweiigg/GRACE_Matlab_Toolbox倦踢；
（2）軟件的第一步其實(shí)是控制文件（GMT_Control_File）腳本的書(shū)寫(xiě)；
控制文件截圖
（3）因?yàn)檫@里處理的是RL06產(chǎn)品侠草，而代碼處理的是RL05數(shù)據(jù)辱挥，所以這里有幾處是需要修改的：二階項(xiàng)C20采用最新的C20_RL06，需要把程序包里面gmt_replace_C20函數(shù)中RL05替換成RL06才可以運(yùn)行边涕；c21和s21替換文件中說(shuō)明文件刪除晤碘，然后將gmt_replace_C21_S21_C22_S22函數(shù)中g(shù)et_days改為gmt_get_days。這里順便多說(shuō)一句功蜓，數(shù)據(jù)文件是從ICGEM下載的最新的园爷，不然可能出問(wèn)題。
（4）點(diǎn)擊運(yùn)行直接出結(jié)果霞赫，這里生成的結(jié)果大概4.5M腮介。
生成結(jié)果
（5）接下來(lái)進(jìn)行泄漏的改正，這個(gè)時(shí)間可能稍微長(zhǎng)一些端衰，要耐心等待叠洗，不要關(guān)閉軟件甘改。
（6）接下來(lái)是將生成的球諧分析，將生成的球諧數(shù)據(jù)文件轉(zhuǎn)換為grid_data.mat：
（7）將上一步生成的格網(wǎng)數(shù)據(jù)轉(zhuǎn)換為時(shí)間序列Grid2Series灭抑，以亞馬遜區(qū)域?yàn)槔?/p>
順便將生成的文件在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');
亞馬遜流域時(shí)間序列變化
（8）繼續(xù)分析十艾，我們可以得到全球趨勢(shì)項(xiàng)變化。
