LoadDef python工具包学习二：利用JPL mascon数据计算地表站点的三维变形（E、N、V）

（1）使用工具箱，计算负荷love数，run_ln.py脚本计算love数；运行run_gf.py计算格林函数

（2）mascon数据的下载与格式转换：工具箱提供的是JPL的mascon数据。下载好数据置于input/Load_Models/GRACE-Tellus-RL06文件夹下。

需要注意的是，在运行run_cn_grace.py脚本之前，需要运行gen_grace_tellus_rl06.py脚本对mascon的nc文件进行格式转换（..\LoadDef-main\LoadDef-main\GRDGEN\load_files）.此处需要修改strat_year_tm和start_month_tm，使其不超过数据长度。loadfile1需要对应改成mascon文件名。

点击运行，即可在..\LoadDef-main\LoadDef-main\output\Grid_Files\nc\GRACE下生成这些文件：

（3）计算对应的三维变形。点击运行run_cn_grace.py，需要修改loadfile_prefix，其为上图的其中某一个文件（只能一个一个输入读取），同时，修改first_data和last_data.

其中，计算点的位置文件为sta_file = ("../input/Station_Locations/NOTA_Select.txt")

点击运行，在..\LoadDef-main\LoadDef-main\output\Convolution目录下生成不同站点在GRACE的某个月份下的三维位移。

 后续可以考虑替换不同的mascon数据，或者替换成其他的地表质量模型，比如土壤水、大气、冰川。或者可以更改站点，或者计算全球的三维位移。