首先在github下載這個(gè)工具箱：https://github.com/hrmartens/LoadDef

下載完畢，解壓縮進(jìn)入如下的界面，點(diǎn)擊進(jìn)入working工作文件夾。在doc文件夾里面有幫助文檔和教程，有需要的可以仔細(xì)閱讀。

loaddef主要的功能模塊包括：

這里我主要是介紹如何使用run_ln.py腳本計(jì)算love數(shù)。在這之前需要配置運(yùn)行環(huán)境，安裝對(duì)應(yīng)的包。打開(kāi)Anaconda Prompt，在命令行輸入以下的命令即可。

接著運(yùn)行run_ln.py，稍等幾分鐘就可以得到對(duì)應(yīng)的勢(shì)love數(shù)、負(fù)荷love數(shù)等。具體可以在output/Love_Number文件夾下找到

下圖是計(jì)算得到的load love數(shù)。

關(guān)于如何計(jì)算格林函數(shù)，以及進(jìn)一步計(jì)算地表質(zhì)量變化產(chǎn)生的負(fù)荷位移，可以分別運(yùn)行working文件夾下的各個(gè)計(jì)算腳本，后續(xù)我將接著更新。

參考資料：

1. Martens, H.R., L. Rivera, and M. Simons (2019), “LoadDef: A Python-based toolkit to model elastic deformation caused by surface mass loading on spherically symmetric bodies,” Earth & Space Science, 6: 1–13, doi:10.1029/2018EA000462.

2. Martens, H.R., and M. Simons (2020), “A comparison of predicted and observed ocean tidal loading in iv Alaska,” Geophys. J. Int., 223(1): 454–470, doi:10.1093/gji/ggaa323.

3. Martens, H.R., L. Rivera, M. Simons, and T. Ito (2016), “The sensitivity of surface mass loading displacement response to perturbations in the elastic structure of the crust and mantle,” J. Geophys. Res. Solid Earth, 121: 3911–3938, doi:10.1002/2015JB012456.

4. Martens, H.R., M. Simons, S. Owen, and L. Rivera (2016), “Observations of ocean tidal load response in South America from sub-daily GPS positions,” Geophys. J. Int., 205(3): 1637–1664, doi:10.1093/gji/ggw087.