1、Possion积分算法公式系列十三使用边界面的高程异常格网反算重力场元也是采用的Possion积分本文是根据边界面的扰动重力格网计算外部空间的扰动重力及扰动重力梯度。由边界面大地高格网m及其面上残差扰动场元格网按Poisson严密数值积分算法计算大地水准面及其外部空间的残差扰动场元。Poisson积分是第一边值问题解只要求边界面连续可微不要求边界面是严格的重力等位面。算法的边缘效应小可有效抑制短波信号衰减适合向上向下解析延拓一次积分性能接近梯度法解析延拓效果。重复调用本模块可实现迭代运算一般情况下迭代1次、至多3次足以满足绝大多数精度要求。2、入口函数Poissonnumericalintegral.f90输入参数dr积分半径m。输入参数calcpntfl计算点文件名。头文件一行记录格式点号/点名 经度度小数 纬度度小数 大地高m......。输入参数sfhhgrdfl边界大地高格网文件名。输入参数gravgrdfl边界上残差扰动重力场元格网文件名。计算重力梯度的入口函数也用这个但是其Possion积分公式是9.24为了抑制积分计算的奇异性需要参照与9.22公式改化9.24公式。具体来说就是修改PoissonBLH函数实现计算扰动重力梯度的计算。公式9.24改化公式可通过deepseek求出如下3、编制程序先看外部扰动重力积分计算计算界面如下图从图中可以看出计算结果与上图的第一行相同。再看外部扰动重力梯度积计算计算界面如下图。下图是章老师程序计算结果比较上下图计算结果基本相同相差小于0.02mGal。文件转换的结果如下图与上图基本相同。参考文献地球外部各类场元Possion积分计算 - 章传银PAGravf4.5参考说明书.pdf