地球物理大地测量

大型科学

计算平台

输入地球重力位系数模型，计算地球空间任意点处的高程异常（m）、空间异常（mGal）、扰动重力（mGal）、垂线偏差向量（s，ʺ/秒，南向、西向）、扰动重力梯度（径向，E）、水平重力梯度向量（E，北向、西向）或Laplace算子（E）的模型值。

输入地球椭球四个基本参数，计算地球椭球主要的几何导出常数和物理导出常数。地球椭球的第四个基本参数，可从地球重力场位系数C20、地球动力学形状因子J2、地球椭球扁率倒数1/f和地球椭球面正常重力位U0中选择其中之一。

由给定经纬度大地高高和时刻，预报地面及地球外部空间计算点10种大地测量要素的固体潮效应。改善IERS2010站位移固体潮算法，实现全空间大地测量全要素固体潮效应统一解析计算。天体的地球引潮位计算时，月球取6阶、太阳取3阶和太阳系地球外部行星取2阶。Intel Fortran源码。

利用全球海潮负荷球谐系数模型（cm），按负荷球谐综合算法，由给定经纬度正常高和时刻，预报瞬时海面潮高（cm），计算地面及地球外部计算点10种大地测量要素的海潮负荷效应。扩展与完善IERS2010第6、7章中的海潮负荷效应算法，实现全空间大地测量全要素海潮负荷效应统一解析计算。Intel Fortran源码。

（1）BelPnmdt(pnm,dpt1,dpt2,maxn,t)     计算规格化连带勒让德函数Pnm及其对t一、二阶导数，规格化Pnm采用改进的Belikov递推算法，一、二阶导数采用非奇异递推算法。测试到7200阶。（2）BelPnm(pnm,maxn,t)      采用改进的Belikov递推算法计算规格化连带勒让德函数Pnm。（3）Pndpn_dt(p,dp,n,t)    计...

利用高阶地球重力场位系数模型，计算过指定点(经度，纬度，大地高)重力等位面相对于指定点的正（常）高差、大地高差和重力差格网。利用局部重力场数据，能进一步精化正（常）高差、大地高差和重力差格网，请参考《PALGrav3.0用户操作手册》“5....

利用高阶地球重力场位系数模型，计算过指定点(经度，纬度，正(常)高)等正（常）高面相对于指定点的重力位差、大地高差和重力差格网。利用局部重力场数据，能进一步精化重力位差、大地高差和重力差格网，请参考《PALGrav3.0用户操作手册》“5....

利用高阶地球重力场位系数模型，计算近地空间点处的模型高程异常（m）、模型大地水准面高（m）以及正高与正常高之差（cm）。利用局部重力场数据，能进一步精化正高与正常高差异，请参考《PALGrav3.0用户操作手册》“5.6正高与正常高差别区域...

利用高阶地球重力场位系数模型，计算地面及近地空间点高程异常径向变化率的模型值（cm/km）。此项功能可用于GNSS水准高程异常改正。利用局部重力场数据，能进一步精化高程异常径向变化率，请参考《PALGrav3.0用户操作手册》“5.4高程异...

采用滑动平均、高斯、指数型或巴特沃斯等低通滤波器，对格网数据进行低通滤波。滤波前后，格网的规格（范围和分辨率）不变。模块源码：GrdLowpFilter.f90

按输入的格网规格和选定的插值权函数形式及参数值，对离散点值数据进行乘权格网化运算。适合各种空间分布差异大、精度不均匀的地球物理观测数据格网化，自适应采样点密度，无边缘效应，有利于构造高性能地球物理场。模块源码：WeighGridate.f9...

由CORS网（GNSS监测网）基线一维分量、水准网测段高差或重力网段差变化量，以用户给定的一组地面站点为拟稳基准，采用附有拟稳基准约束的间接最小二乘平差法，计算站点坐标分量变化、正（常）高变化或重力变化平差值。输入输出与参数设置参见《PAL...