static_syn
- 简介:
指定震源机制,根据静态格林函数合成三分量位移(及其空间导数)
语法
grt static syn -Gingrid -S[u]scale -Ooutgrid [ -Xx1/x2/dx ] [ -Yy1/y2/dy ] [ -Ffn/fe/fz ] [ -Mstrike/dip/rake ] [ -TMxx/Mxy/Mxz/Myy/Myz/Mzz ] [ -N ] [ -e ] [ -h ]
描述
关于输出分量、震源机制的部分与 syn 模块相同,但输入和输出均使用 NetCDF 网格。
static syn 模块计算时使用的 XY 网格默认继承 static_greenfn 模块设置的 -X 和 -Y 。 也支持设置 -X 和 -Y 来指定新的 XY 网格,此时每个节点将使用最近震中距的格林函数来进行合成, 因此这是一种近似方法,程序会打印震中距之差的统计信息,以方便评估近似效果。 但近似好处是一旦静态格林函数计算好,合成阶段使用新的网格也可以复用格林函数,无需再计算。
必选选项
- -Gingrid
存有静态格林函数的网格文件。
- -S[u]scale
放大系数 scale ,作用在合成结果上。
对于爆炸源、剪切源和矩张量源, scale 即标量地震矩,单位为 dyne-cm 。
对于单力源, scale 单位为 dyne 。
由于标量地震矩公式中 \(M_0 = \mu\Sigma D\) 包含震源处的 \(\mu\) , 这里为了方便设置参数,用户可以使用 -Su 表示震源处的 \(\mu\) , 而 scale 就仅需设置 \(\Sigma D\) (面积 x 滑动距离,单位 cm^3)。
- -Ooutgrid
输出网格文件,若同名文件已存在,则销毁后新建。
可选选项
- -Xx1/x2/dx
指定北向等距网格点 (km) 。 x1,x2 分别为最小、最大值, dx 为间隔。
- -Yy1/y2/dy
指定东向等距网格点 (km) 。 y1,y2 分别为最小、最大值, dy 为间隔。
- -Mstrike/dip[/rake]
设置剪切源/张裂源参数,角度单位为 °
strike - 断层走向, 0~360
dip - 断层倾角, 0~90
rake - 滑动角, -180~180 [不指定 rake 即设置张裂源]
- -TMxx/Mxy/Mxz/Myy/Myz/Mzz
设置矩张量源的 6 个分量,其中下标代表 x (北向),y (东向),z (向下)。 关于 z 方向,详见 Warning 。
- -N
输出的三分量为 ZNE [默认为 ZRT]。
- -e
在合成位移的同时,也合成其空间导数 \(\dfrac{\partial (u_z, u_r, u_\theta)}{\partial (z, r, \theta)}\) 。 偏导对应在文件名/变量名开头添加了
z,r和t。合成的位移空间导数的单位均为 1 。
- -h
打印帮助文档。
示例
详见教程: