（1）使用工具箱，計(jì)算負(fù)荷love數(shù)，run_ln.py腳本計(jì)算love數(shù)；運(yùn)行run_gf.py計(jì)算格林函數(shù)

（2）mascon數(shù)據(jù)的下載與格式轉(zhuǎn)換：工具箱提供的是JPL的mascon數(shù)據(jù)。下載好數(shù)據(jù)置于input/Load_Models/GRACE-Tellus-RL06文件夾下。

需要注意的是，在運(yùn)行run_cn_grace.py腳本之前，需要運(yùn)行g(shù)en_grace_tellus_rl06.py腳本對mascon的nc文件進(jìn)行格式轉(zhuǎn)換（..\LoadDef-main\LoadDef-main\GRDGEN\load_files）.此處需要修改strat_year_tm和start_month_tm，使其不超過數(shù)據(jù)長度。loadfile1需要對應(yīng)改成mascon文件名。

點(diǎn)擊運(yùn)行，即可在..\LoadDef-main\LoadDef-main\output\Grid_Files\nc\GRACE下生成這些文件：

（3）計(jì)算對應(yīng)的三維變形。點(diǎn)擊運(yùn)行run_cn_grace.py，需要修改loadfile_prefix，其為上圖的其中某一個文件（只能一個一個輸入讀?。?，同時，修改first_data和last_data.

其中，計(jì)算點(diǎn)的位置文件為sta_file = ("../input/Station_Locations/NOTA_Select.txt")

點(diǎn)擊運(yùn)行，在..\LoadDef-main\LoadDef-main\output\Convolution目錄下生成不同站點(diǎn)在GRACE的某個月份下的三維位移。

?后續(xù)可以考慮替換不同的mascon數(shù)據(jù)，或者替換成其他的地表質(zhì)量模型，比如土壤水、大氣、冰川?；蛘呖梢愿恼军c(diǎn)，或者計(jì)算全球的三維位移。