​地球物理大地测量学计算系列之十五等位面扰动重力梯度积分扰动重力计算

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