格林积分负荷形变场大地测量全要素计算

由区域地面等效水高变化（大气压、陆地水、江河湖库水或海平面变化等，mm）格网模型，按负荷格林函数积分法，计算地面或近地空间任意点的高程异常/大地水准面（mm），地面扰动重力（μGal），地倾斜南分量（ms）、地倾斜西分量（ms），地面重力（μGal）、地面东方向（mm）、地面北方向（mm）、地面径向（大地高，mm）或地面正常高（mm）的非潮汐负荷形变。

江河湖库水大地测量全要素负荷形变场计算示例见ETideLoad3.0练习3。

程序可与“4.6等效水高球谐综合负荷形变场计算”配合使用，构造剩余区域地面等效水高格网模型，采用移去恢复方法，计算区域高精度负荷形变场及时变重力场。输入输出请参考《ETideLoad3.0用户操作手册》4.8节。

模块源码：LoadGreenIntg.f90

算法公式：《ETideLoad3.0用户操作手册》式（8.4.1）~（8.4.16）。

地球物理模型：

（1）超高阶负荷勒夫数loadlove.dat（20000阶）。

调用模块：

（1）区域负荷格林函数数值计算模块LGreenF.f90。

（2）负荷等效水高直接影响积分计算模块LoadrDirectFunc.f90。

（3）负荷等效水高间接影响积分计算模块LGIntglr.f90。

（4）正常重力场计算模块GNormalfd.f90。

（5）大地测量三维坐标形式变换模块GRSConvert.f90。

（6）勒让德函数及其导数计算模块Pndpn_dt2.f90。