1、Possion积分算法公式
由边界面大地高格网（m）及其面上残差高程异常（m）格网，计算大地水准面或地球外部的残差空间异常（mGal）、残差扰动重力（mGal）与残差垂线偏差向量（″，SW）。（可从别人大地水准面成果中反演高精度局部重力场，偷笑）。
程序采用Poisson积分与场元微分组合算法，实现高程异常反运算，不要求边界面是重力等位面。
2、入口函数
AntiHeightanomalytogravf.f90
输入参数dr－积分半径（m）。
输入参数calcpntfl－计算点文件名。头文件一行，记录格式：点号/点名 经度（度小数） 纬度（度小数） 大地高（m）......。
输入参数boundryfl－输入残差高程异常所在边界面的大地高格网文件名。
输入参数ksigrdfl－边界面残差高程异常格网文件名。
3、程序编制
CAD VB.net编制的主界面如下图，图中点的计算结果与上图的第一行相同。
转换的文件结果对比如下图，汉字标题是本程序计算结果，英文标题是章老师程序计算结果，结果完全相同。
参考文献
https://www.zcyphygeodesy.com/h-nd-67.html
PAGravf4.5参考说明书.pdf