描述
What and WHere
README.md 开始信息
LICENSE GNU LGPL licence
StdModels 在/lib/TauP-2.6.1.jar/edu/sc/seis/TauP/StdModels中存在的模型
创建模型
.nd模型
模型中，每一列为：
depth（km）；P波速度（km/s）；S波速度(km/s)；密度（g/cm3）；Qp；Qs。
但是只有前三列为有效值，其余三列无效，只是为了兼容之后的某些情况。
注意！模型最后一行留一空白行。
使用英文单词表示分界面：
- mantle 莫霍界面
- outer-core CMB
- inner-core ICB
示例：1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18/* below is a simple named discontinuities model. */
0.0 5.0 3.0 2.7
20 5.0 3.0 2.7
20 6.5 3.7 2.9
33 6.5 3.7 2.9
mantle # the word "mantle" designates that this is the moho
33 7.8 4.4 3.3
410 8.9 4.7 3.5
410 9.1 4.9 3.7
670 10.2 5.5 4.0
670 10.7 5.9 4.4
2891 13.7 7.2 5.6
outer-core # "outer-core" designates that this is the core mantle boundary
2891 8.0 0.0 9.9
5149.5 10.3 0.0 12.2
inner-core # "inner-core" makes this the inner-outer core boundary
5149.5 11 3.5 12.7
6371 11.3 3.7 13
.vel模型
前两行为模型描述，两行之后为模型。
间断面使用两行相同的深度表示。每一行为：
depth；P波速度；S波速度；密度 （单位同上）
示例：
1 | ak135 - P |
创建taup模型
创建完成模型后，可以通过语句完成taup模型创建
taup create -nd prem_normal.nd --verbose
生成的模型名称为prem_normal.taup
错误信息
创建模型时，注意地表处S波速度不能为0,否则报错。
1 | Error starting tool: create edu.sc.seis.TauP.SlownessModelException: TauP Exception: SlownessModel Exception: Unable to handle zero S velocity layers at surface. This should be fixed at some point, but is a limitation of TauP at this point |
S波速度与P波速度关系
S波速度小于P波速度。因此大于时会报错。
1 | Caused by: edu.sc.seis.TauP.VelocityModelException: TauP Exception: VelocityModel Exception: S velocity, 3.19101 at depth 0.0 is greater than the P velocity, 1.44996 |
计算震相走时
交互输入
taup time --mod prem_normal.nd
示例：
1 | Enter: |
也可以使用命令行命令直接进行计算：
-mod 输入模型名称
-h 震源深度 km
–stadepth 地震仪深度 km
-deg 震中距 度
-km 震中距 千米
或者给出经纬度（假设地球为球体）
-sta lat lon 台站经纬度
-evt lat lon 事件经纬度
–rayp 只输出射线参数
–time 只输出走时
–rel phasename 同时输出相对走时？
–json json格式输出走时
-ph 震相名称
-pf 震相文件，即包含震相名称的文件
示例：
taup time -mod prem_normal.nd -h 200 -deg 170 --rel P --json
1 | { |
截至此时，其只能在命令行中输出信息。
绘制速度剖面图
taup velplot -nd prem_normal.nd
sh prem_normal_vel.gmt
即可使用gmt绘制出漂亮的速度剖面图。
- 模型类型
- -mod 直接输入Taup中自带的模型名称。不需要调用其他文件。例如：
iasp91prem等。可以在安装目录StdModel文件夹中查看。 - -nd 输入文件名称。
- -tvel
- –version 在命令行输出Taup版本。
- –debug 可能是出现bug会在命令行有输出？目前还没有创建bug测试。
- –verbose 未知
注意：必须已安装gmt。安装见《GMT安装》。
绘制模型慢度剖面图
绘制走时曲线
1 | export model='premoonaver' |
注意：震相s与S不同，需要单独绘制。