✍️ 周茂  •  📅 2025-12-30

earthtide

官方文档:

earthtide

简介:

计算固体地球潮汐网格或者时间序列

计算固体地球潮汐的三个分量,输出网格或者时间序列。同时可以输出太阳和月亮 的经纬度位置。输出可以为网格或者表文件(标准输出),表文件的格式为: 时间 北向分量 东向分量 垂直分量 ,分量的单位为米。

语法

gmt earthtide -Goutgrid -Llon/lat -S [ -Cx|e,y|n,z|v ] [ -Iincrement ] [ -Rregion ] [ -T[min/max/]inc[+i|n]|-Tfile|list ] [ -V[level] ] [ -rreg ] [ --PAR=value ]

必须选项

以下三个选项必须至少使用其中一个。

-G

-Goutgrid[=ID][+ddivisor][+ninvalid][+ooffset|a][+sscale|a][:driver[dataType][+coptions]]

输出网格文件名。通过追加 =ID 可指定 网格格式。 参数详细介绍请参考 写 netCDF 文件

-S
-S

输出太阳和月亮的地理坐标以及其距离,单位为米。输出为 Mx7 的矩阵, M 为时刻的个数(通过 -T 设置),7 列分别代表:时间,太阳经度, 太阳纬度,太阳距离,月亮经度,月亮纬度,月亮距离

-L
-Llon/lat

计算指定位置上的固体地球潮汐的时间序列。坐标必须为大地坐标,默认 使用 GRS80 椭球(在亚微米的精度水平可以认为与 WGS-84 椭球等同)

可选选项

-C
-Cx|e,y|n,z|v

设置需要计算的固体地球潮汐分量及其保存的网格。该选项需要 -G 选项。 不同选项之间使用逗号分隔:xe 为东向分量;yn 为北向分量,zv 为垂直分量。例如,-Ce,v 将输出 写入两个网格,一个为东向分量网格,一个为垂直分量网格。如果设置了 -G ,但不指定 -C ,则默认为垂直分量。

-I
-Ixinc[+e|n][/yinc[+e|n]]

指定X和Y方向的网格间隔

  • xincyinc 为 X 和 Y 方向的网格间隔。对于地理坐标,可以指定网格间隔单位 [默认单位为度]

  • +e 微调X和Y方向范围的最大值,使得其是网格间隔的整数倍(默认会微调 网格间隔以适应给定的数据范围)

  • +n 表明 xincyinc 不是网格间隔,而是X和Y方向的节点数。 此时会根据节点数、网格区域范围以及网格配准方式重新计算网格间隔。

注意:

  • yinc 设置为0,则表示其与 xinc 相同

  • 若使用 -Rgrdfile 选项,则网格间隔和配准方式已经根据网格文件自动初始化, 此时依然可以使用 -I-r 覆盖相应的值

-R
-Rxmin/xmax/ymin/ymax[+r][+uunit]

指定数据范围。 (参数详细介绍)

-T
-T[min/max/]inc[+i|n]|-Tfile|list

inc 为步长,生成从 minmax 均匀分布的时间序列。 附加 +i 表示给定的是步长的倒数,或者附加 +n 表示 inc 为序列的元素个数。 可以为步长附加一个有效的时间单位(d|h|m|s)。 如果仅给出了 min,则我们使用该日期和时间进行计算。如果未提供 -T ,则从计算机时钟获取 UTC 格式的当前时间。 如果未提供 -G-S ,则 -T 被解释为在 -L 指定的位置计算时间序列,因此 -L 变为必选。 当同时存在 -G-T 时,仅考虑第一个时间 T 序列。 日期范围可以从 1901 年到 2099 年。有关数组生成的详细信息,请参见 生成一维数组

-V
-V[level]

设置 verbose 等级 [w]。 (参数详细介绍)

-r
-r[g|p]

设置网格配置方式 [默认为网格线配准]。 (参数详细介绍)

-^-

显示简短的帮助信息,包括模块简介和基本语法信息(Windows下只能使用 -

-++

显示帮助信息,包括模块简介、基本语法以及模块特有选项的说明

-? 或无参数

显示完整的帮助信息,包括模块简介、基本语法以及所有选项的说明

--PAR=value

临时修改GMT参数的值,可重复多次使用。参数列表见 配置参数

示例

计算 2018 年 6 月 18 日中午,固体地球潮汐垂直分量网格,网格的分辨率为 30 m(使用默认的 -R-I 选项)

gmt earthtide -T2018-06-18T12:00:00 -Gsolid_tide_up.grd

计算一天的时间序列,起始时间与上例相同,计算位置位于 (-7W,37N),时间 步长为 1 分钟

gmt earthtide -T2018-06-18T/2018-06-19T/1m -L-7/37 > solid_tide.dat

获取当前时间太阳和月亮的地理坐标

gmt earthtide -S

注意事项

  1. 所有的输入和输出时间都必须是 UTC 时

参考

http://geodesyworld.github.io/SOFTS/solid.htm