地球物理大地测量
大型科学
计算平台

CASM.jpg

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

 二维码 47
发表时间:2023-09-25 17:02作者:章传银来源:大地测量科学计算网址:http://www.zcyphygeodesy.com

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

幻灯片04.JPG幻灯片07.JPG

测试程序:

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源码压缩包下载