空天地海各种场元剩余地形影响球谐综合计算

由全球地形面密度规格化球谐系数（kg/m2）模型，按指定的球谐系数阶数范围，计算地面或地球外部任意点BLH(3)的高程异常（m）、空间异常（mGal）、扰动重力（mGal）、垂线偏差向量（s，ʺ/秒，南向、西向）、（垂向）扰动重力梯度（E）、水平重力梯度向量（E，北向、东向）或扰动位（m2/s2）的剩余地形影响。

模块可实现卫星重力梯度、海洋测高垂线偏差、海洋重力、航空重力等空天海地统一、各种场元解析相容的地形处理能力。

模块源码：HarmFdCalc.f90

算法公式：《PALGrav3.0用户操作手册》式（7.1）~（7.4）。

地球物理模型：全球地形面密度球谐系数1800.dat

调用模块：

（1）正常重力场参数计算模块GNormalfd.f90

（2）大地测量坐标变换包GRSConvert.f90

（3）超高阶缔合勒让德函数及其导数计算模块BelPnmdt.f90