由重力位系数模型计算任意点全要素重力场参数

  由超高阶地球重力场位系数模型，计算地面及地球外部空间任意点的高程异常（m）、空间异常（mGal）、扰动重力（mGal）、垂线偏差向量（s，ʺ/秒，南向、西向）、扰动重力梯度（径向，E）、水平重力梯度向量（E，北向、西向）和Laplace算子（E）的模型值。

测试程序：

Fortran源码GMfullelementcalc.f90（固定格式，每行代码最大132字符）。

DOS可执行文件GMfullelementcalc.exe（将压缩包中全部6个.f90文件编译并连接后，直接生成可执行文件）。

调用功能模块（Fortran源码）：

（1）大地测量坐标变换模块集GRSConvert.f90

（2）全要素扰动重力场参数计算模块RntGravFdpnm.f90

（3）正常重力场计算功能模块集NormGravfd.f90

（4）超高阶勒让德、缔合勒让德函数及其导数快速计算模块集Legendre.f90

算法公式：PAGravf4.5cn参考说明书7.1、7.2和7.3节。

文件格式：1.4.1大地测量数据文件格式约定。

提示：地球两极处，垂线偏差西向、水平重力梯度西向无定义，计算结果无意义。

Fortran源码压缩包下载