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

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

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

输入文件：

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

大地水准面上高程异常，即为大地水准面差距或大地水准面（大地）高。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源码压缩包下载