1、过指定点重力等位面构造利用高阶地球重力场位系数模型计算过指定点(经度纬度大地高)重力等位面相对于指定点的大地高差和重力差格网。利用局部重力场数据能进一步精化过指定点重力等位面的大地高差和重力差格网。[功能]利用参考地球重力场位系数模型计算过指定点(B,L,H)重力等位面的重力(mGal)和大地高(m)模型值格网。[输入]重力等位面范围格网文件参考重力场位系数模型文件。等位面范围格网文件只用于确定等位面格网的经纬度范围和分辨率程序忽略格网的格值数据。位系数模型文件第1行约定为位系数尺度参数:地心引力常数GM与椭球长半轴a。尺度参数(GM、a)是位系数模型的表示基准。[参数设置]输入指定点大地坐标(经纬度大地高)与位系数模型最大计算阶数。程序要求指定点落在等位面格网范围内。位系数模型阶数不宜过大程序需要迭代逼近计算耗费时间较长。[输出]过指定点重力等位面的重力和大地高模型值格网文件。过指定点重力等位面的重力模型值格网文件过指定点重力等位面的大地高模型值格网文件2、程序编制由于无需考虑高程系统参考文献1的代码需要修改可根据上面的公式进行修改也可通过deepseek修改程序主界面如下图计算的大地高、大地高差和重力差见下图。大地高大地高差大地高-指定点的大地高重力差从上述文件可以看出大地高差、指定点的重力位和重力与章老师程序一致但是重力差不一致。指定点的重力和重力位为978718.0644 62623481.0494重力差的第一点经纬度和大地高为27.0208333333333106.0208333333331335.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