1、过指定点重力等位面构造
利用高阶地球重力场位系数模型，计算过指定点(经度，纬度，大地高)重力等位面相对于指定点的大地高差和重力差格网。利用局部重力场数据，能进一步精化过指定点重力等位面的大地高差和重力差格网。
[功能]利用参考地球重力场位系数模型，计算过指定点(B,L,H)重力等位面的重力(mGal)和大地高(m)模型值格网。
[输入]重力等位面范围格网文件，参考重力场位系数模型文件。
等位面范围格网文件，只用于确定等位面格网的经纬度范围和分辨率，程序忽略格网的格值数据。
位系数模型文件第1行约定为位系数尺度参数:地心引力常数GM与椭球长半轴a。尺度参数(GM、a)是位系数模型的表示基准。
[参数设置]输入指定点大地坐标(经纬度大地高)，与位系数模型最大计算阶数。
程序要求指定点落在等位面格网范围内。位系数模型阶数不宜过大，程序需要迭代逼近计算，耗费时间较长。
[输出]过指定点重力等位面的重力和大地高模型值格网文件。
过指定点重力等位面的重力模型值格网文件
过指定点重力等位面的大地高模型值格网文件
2、程序编制
由于无需考虑高程系统，参考文献1的代码需要修改，可根据上面的公式进行修改，也可通过deepseek修改，程序主界面如下图：
计算的大地高、大地高差和重力差，见下图。
大地高
大地高差，大地高-指定点的大地高
重力差
从上述文件可以看出：大地高差、指定点的重力位和重力与章老师程序一致，但是重力差不一致。指定点的重力和重力位为：978718.0644 ； 62623481.0494
重力差的第一点经纬度和大地高为：27.0208333333333，106.020833333333，1335.9711
根据以前的程序计算第一点的正常位和扰动位、正常重力和扰动重力，见下图。
正常位+扰动位=62623780.97+ （-299.9164 ）= 62623481.05，与指定点重力位相等。
正常重力+扰动重力=978687.1231 + （ -16.1485）= 978670.9746
第一点重力-指定点重力=978670.9746 - 978718.0644= -47.0898
参考文献
过指定点模型重力等位面构造 - 章传银
利用重力场模型计算等正(常)高面的重力位差和大地高差 - 章传银
利用重力场模型计算正高与正常高差异 - 章传银
利用重力场模型计算高程异常径向变化率 - 章传银
Molodensky边值校正Fortran代码 - 章传银
正高与正常高的定义与差别-CSDN博客
PAGravf4.5参考说明书.pdf
)
