LoadDef python工具包学习一：计算love数

首先在github下载这个工具箱：https://github.com/hrmartens/LoadDef

下载完毕，解压缩进入如下的界面，点击进入working工作文件夹。在doc文件夹里面有帮助文档和教程，有需要的可以仔细阅读。

loaddef主要的功能模块包括：

 

 这里我主要是介绍如何使用run_ln.py脚本计算love数。在这之前需要配置运行环境，安装对应的包。打开Anaconda Prompt，在命令行输入以下的命令即可。

conda config --prepend channels conda-forge 
 conda config --set channel_priority strict 
 conda create -n LDEF python=3.7 basemap mpi4py scipy netCDF4 pyshp shapely numba jupyter pygmt

 接着运行run_ln.py，稍等几分钟就可以得到对应的势love数、负荷love数等。具体可以在output/Love_Number文件夹下找到

 

 下图是计算得到的load love数。

 关于如何计算格林函数，以及进一步计算地表质量变化产生的负荷位移，可以分别运行working文件夹下的各个计算脚本，后续我将接着更新。

参考资料：

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.