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选项。 不同选项之间使用逗号分隔:x 或 e 为东向分量;y 或 n 为北向分量,z 或 v 为垂直分量。例如,-Ce,v 将输出 写入两个网格,一个为东向分量网格,一个为垂直分量网格。如果设置了-G,但不指定-C,则默认为垂直分量。
- -I
- -Ixinc[+e|n][/yinc[+e|n]]
指定X和Y方向的网格间隔
xinc 和 yinc 为 X 和 Y 方向的网格间隔。对于地理坐标,可以指定网格间隔单位 [默认单位为度]
+e 微调X和Y方向范围的最大值,使得其是网格间隔的整数倍(默认会微调 网格间隔以适应给定的数据范围)
+n 表明 xinc 和 yinc 不是网格间隔,而是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 为步长,生成从 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
注意事项
所有的输入和输出时间都必须是 UTC 时