1、扰动重力梯度积分计算公式
扰动重力梯度积分计算扰动重力的公式9.34与广义Hointe积分公式形式一样，它们使用同一个函数HotineBLH。
2、入口函数
入口函数可以參照广义Stokes/Hotine数值积分外部高程异常计算的入口函数StokesHotinenumintegral.f90。
功能由等位面大地高格网及其面上残差扰动重力梯度(E)格网，按严密积分方法，计算大地水准面或地球外部空间点的残差扰动重力(mGal)。
输入格网规格相同的等位边界面大地高及其残差扰动重力梯度(E)格网文件，空间计算点位置文件或计算面大地高文件。
空间计算点位置文件记录格式:点号/点名经度(度小数)纬度(度小数)大地高
参数设置选择计算点文件格式，输入积分半径。
输出残差扰动重力积分值结果文件。
当输入空间计算点位置文件时，输出结果文件记录格式:在源计算点值文件记录的基础上，增加一列该点的残差扰动重力积分值，保留4位有效数字。当输入计算面大地高格网计算，程序输出与残差扰动重力梯度格网规格相同的残差扰动重力格网。
3、程序编制
程序编制完成以后，计算结果与上图的结果相差较大。后来又用参考文献1的代码，仍然有问题。怀疑核函数既不是使用公式9.5，也不是参考文献1的核函数，在deepseek中，输入参考文献1的代码，让它指出问题所在并给出解决方案，按照它给出的方法还是不对。
随后，咨询了一下章老师，他说用的是广义Hotine积分公式。
在deepseek中，输入Hotine积分代码，说明代码是广义Hotine积分，根据扰动重力格网，求高程异常，如何改造一下，根据扰动重力径向梯度格网，求扰动重力？
deepseek这次给出了正确的核函数，见下图，并给出改造代码。
程序界面如下图，计算结果上图第一行相差0.2mGal。
单点比较没有说服力，看看文件，见下图。程序计算结果（汉字标题）与章老师计算结果（英文标题）基本相同，相差约1%。
参考文献
PAGravf4.5参考说明书.pdf