地球物理大地测量学计算系列之二十二过指定点等高面的重力位计算

1、过指定点等高面的重力位计算利用高阶地球重力场位系数模型计算过指定点(经度纬度正(常)高)等正常高面相对于指定点的重力位差、大地高差和重力差格网。[功能]利用参考地球重力场位系数模型计算过指定点(B,L)等正(常)高面的重力位(m2/s2)、大地高(m)和重力(mGal)模型值格网。[输入]等高面范围格网文件参考重力场位系数模型文件。等高面范围格网文件只用于确定等高面格网的经纬度范围和分辨率程序忽略格网的格值数据。位系数模型文件第1行约定为位系数尺度参数:地心引力常数GM(1014m3/s2)与椭球长半轴a(m)。尺度参数(GM、a)是位系数模型的表示基准。[参数设置]输入指定点经纬度和正(常)高选择高程系统类型输入位系数模型最大计算阶数。程序要求指定点落在等高面格网范围内。参考重力场模型模型阶数不宜过大程序需要迭代计算时间较长。[输出]过指定点等高面的重力位、重力和大地高模型值格网文件。等高面大地高模型值大地高模型格网格值头文件第9个数;等高面重力位模型值重力位模型格网格值头文件第10个数。等高面重力模型值等高面重力位模型值等高面大地高模型值2、程序编制参考了文献2的代码程序主界面如下图生成重力位差、重力差以及大地高差三个结果文件。过指定点等高面的重力位计算结果文件见下图结果与章老师程序计算相同。等高面重力模型值等高面重力位模型值等高面大地高模型值参考文献过指定点模型重力等位面构造 - 章传银利用重力场模型计算等正(常)高面的重力位差和大地高差 - 章传银利用重力场模型计算正高与正常高差异 - 章传银利用重力场模型计算高程异常径向变化率 - 章传银Molodensky边值校正Fortran代码 - 章传银正高与正常高的定义与差别-CSDN博客PAGravf4.5参考说明书.pdf