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

CASM.jpg

广义Stokes/Hotine积分外部高程异常计算

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

利用广义Stokes(/Hotine)积分公式,由等位面上的残差空间异常(/扰动重力,mGal)格网,输入积分半径(km),采用数值积分、二维FFT或一维FFT算法,计算指定面上的残差高程异常(m)。

为实现有限半径积分,通常需采用参考重力场移去恢复法,先移去等位边界面上的模型空间异常/扰动重力,再积分得到计算点处的残差高程异常,最后恢复计算点处的模型高程异常。

幻灯片42.JPG幻灯片43.JPG

输入文件:

格网规格相同的等位面大地高格网、等位面上残差空间异常/扰动重力格网、以及计算面大地高格网文件。

大地水准面上高程异常,即为大地水准面差距或大地水准面(大地)高。Stokes边值问题要求,边界面必须是重力等位面,即空间异常/扰动重力必须位于等位面上。

测试程序:

Fortran源码StokesHotineIntgrl.f90(固定格式,每行代码最大132字符)

DOS可执行文件StokesHotineIntgrl.exe(将压缩包中全部8个.f90文件以及intel mkl fftw3.f编译并连接后,生成可执行文件)。

输出残差高程异常格网结果文件output.dat(默认),应扣除外围存在积分边缘效应的格网区域。

调用功能模块(Fortran源码):

(1)大地测量坐标变换模块集GRSConvert.f90

(2)正常重力场计算功能模块集NormGravfd.f90

(3)广义Stokes积分计算单点高程异常模块StokesBLH.f90

(4)广义Stokes积分FFT算法StokesFFT.f90

(5)广义Hotine积分计算单点高程异常模块HotineBLH.f90

(6)广义Hotine积分FFT算法HotineFFT.f90

算法公式:PAGravf4.5cn参考说明书中式(9.1)~式(9.7)。


Fortran源码压缩包下载