✍️ 周茂 • 📅 2025-12-30

# earthtide

官方文档:

  earthtide

简介:

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

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

## 语法

gmt earthtide `-G` outgrid `-L` lon/lat `-S` [ `-C` x|e,y|n,z|v ] [ `-I` increment ] [ `-R` region ] [ `-T`[min/max/]inc[+i|n]|`-T` file|list ] [ `-V`[level] ] [ `-r` reg ] [ --PAR=value ]

## 必须选项

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

-G

-G outgrid[=ID][+d divisor][+n invalid][+o offset|a][+s scale|a][:driver[dataType][+c options]]

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

-S

-S

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

-L

-L lon/lat

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

## 可选选项

-C

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

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

-I

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

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

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

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

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

  注意：

  - 若 yinc 设置为0，则表示其与 xinc 相同

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

-R

-R xmin/xmax/ymin/ymax[+r][+u unit]

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

-T

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

  以 inc 为步长，生成从 min 到 max 均匀分布的时间序列。 附加 +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