wangliang1989 / oh-my-cap Goto Github PK

#!/bin/bash
psmeca -P -JX0/10 -R-5/5/-5/5 -Sm8 -G100 -Y9.5 -X-0.7 > 1.ps <<EOF
0 0 0 -0.3 -1 -0.650 0 0 0 17
1 1 0 -0.3 -1 -0.649 0 0 0 17
EOF
ps2raster -E1080 -P -Tf 1.ps

使用以上内容作图，可以看到两个机制解其中一个黑色，然而这两个机制解应该均为白色。我这边作图结果是1.pdf

考虑在生成权重文件的脚本中加一段判断：如果 t9 不是 -12345，这个数据权重为0。

计划如下修改：
考虑加强权重设置的功能

1. 考虑到可能需要在工作区内放置一些文件，比如自己写的脚本，而又不想加入 repo 中，可
   以用 i_ 开头，则会忽略

请教一个问题：计算单力源和爆炸型震源时，程序报错。非常感谢！
命令行为：$ fk.pl -Mhk/15/k -N1024/0.1 -S0 -U0 10。（计算爆炸型震源时）。 程序默认的-S2（双力偶型震源）
报错：Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
Backtrace for this error:
#0 0x7f8eac02532a
#1 0x7f8eac024503
#2 0x7f8eab8b9fcf
#3 0x559541734f9f
#4 0x559541734922
#5 0x5595417329db
#6 0x55954173380b
#7 0x7f8eab89cb96
#8 0x5595417313e9
#9 0xffffffffffffffff
hk_15/10.grn.0
Unable to open hk_15/10.grn.0
hk_15/10.grn.5
Unable to open hk_15/10.grn.5

老师您好！请问gCAP需要使用的商业软件 Numerical Recipes在哪里可以购买

可能会更换模型，修改定义模型的方式，便于修改模型。另外需要确认给 fk 计算的模型文件一定会正确。

1. 格林函数的计算，增加并行计算的脚本
2. 安装方法考虑为所用计算机账户安装的问题

我在教程中看到：数据必须去除仪器响应得到位移或速度记录，单位为 cm 或 cm/s
我在使用seed文件中提取的PZ文件去除仪器响应后，（得到的结果据说单位是m）将结果乘以100，然后按照流程进行得到的反演结果看起来不太对。很多波形图超出边界，且计算出的震级与地震目录中的震级相差较大。
请问通过PZ文件去除仪器响应后是否有不一样的处理步骤？

打算重写波形拟合图，拟合误差深度图和新写台站位置示意图

数据预处理：
DT CAP里的-H参数，重采样的周期
FREQ 去仪器相应的四个频率

反演：
DEPTH 反演的深度
GREEN -G参数
MODEL 格林函数model
PLOT cap的-p参数，绘图
SHIFT cap -s
P_WINDOW 反演P窗口
S_WINDOW S窗口
WEIGHTFILE 权重文件
P_FILTER CAP -c参数里的滤波
S_FILTER 同上

以下三个我想不出合适名字：
-D
-W: 1
-X: 10

反演的矩张量结果保存在.out的结果文件中# MomentTensor那一行。
官方脚本是这样传入gmt绘图的：
$plt3 = "| psmeca -JX1/1 -R-1/1/-1/1 -Sm8 -G100 -Y9.5 -X-0.7 -O -K >> $outps" if $tensor[1] eq "MomentTensor";
open(PLT, $plt3);
printf PLT "0 0 0 @tensor[9,4,7,6] %f %f 17\n",-$tensor[8],-$tensor[5];
因为GMT的psmeca中，所用的是RTP/RTF/USE坐标系。所以得到对应关系（MOPAD默认用NED）：
9 4 7 6 -8 -5 = mrr, mtt, mff, mrt, mrf, mtf = dd nn ee dn -de -ne
或者说：
4 5 6 7 8 9 = nn ne dn ee de dd

1. 把 inversion.md 中对参数的介绍合并到 parameter.md 中，然后删除 inversion.md 中的部分，因为 inversion 的目的只是为了让读者对整个流程有一个认识而已。参数可以介绍，但没有必要这么复杂，想知道每个参数是什么含义就去看 parameters.
2. parameter.md 这个文件的格式要重新整理一下，现在的文档完全没有格式可言；
3. 搞不定的参数就空哪里空着，不要删除

add rdseed

1. 确认 TauP 在整个流程中的作用
2. 若 TauP 非必须，则需要在 marktime.pl 中对是否安装了 TauP 进行检测。若未安装，则直接跳过标记理论到时

老师您好，我是用里面的gcap时一直有问题
#报错：can‘t exec 'cap.pl':NO such...
请问是gcc要降到3.4版本的吗
另外，请问咱们是否有qq群这样的交流平台大家可以方便交流

现在脚本是按一律算爆炸源来处理，应该增加是否计算的选项

配置参数化后更新文件的工作之前没有做完，做完后注意修改历史记录！

1. 目前怀疑机制球台站位置绘制有错
2. 如果台站数很多，可能不能全部画上

以上问题需要检查

震源深度和模型速度的变换和地震学引论中一致。但是，模型界面深度和地震学引论中公式不同：
代码中的公式是 h'=r_0/(r_0-d+1/2h)，h是变换前的地层厚度，d是地层上界面的深度，r_0是地球半径，h'是变换后的地层厚度。
书中没有直接给出公式，可用深度的公式做推导：

现在只实现了原depth.pl的做图功能，输出深度反演的功能还未实现

输出含义也未完全弄清

在反演时，是权重文件中的台名+分量来搜索sac文件的，对此要考虑有loc的时候，文件名协调的问题

共三个地方涉及到了 dt：

1. 格林函数
2. 数据
3. cap.pl 的 -H 选项

需要做的是：

1. 确认三个 dt 是否一定要相等
2. 若1成立，则修改源码使得其检测这一条件是否成立

marktime.pl 未考虑不同分量长度不同的情况而造成 cut 数据出错

根据 gCAP 的官方版本中包含的例子中的sac文件的事件时间和经纬度信息，基本确认该事件是 IRIS 上的 2725882事件，但相关信息并不完全一致，请@seisman 确认

get_depth.pl 调用 Zhu 的脚本，没有考虑到路径下可能有多个model的结果

1. 新建脚本 build.pl: 读取事件目录文件，根据文件内容找到seed文件，然后建立对应的工作目录，把相应的事件信息保存为文本，将这个文本文件和seed文件移动到工作目录中去。
2. 新建脚本do_process.pl: 根据事件目录文件，去调用process.pl预处理各个事件。
3. 新建脚本do_inversion.pl:根据事件目录文件，去调用inversion.pl反演各个事件。

目前 hk 中的四列数据分别为：厚度、vs、波速比、Qs，这不是一个常见的速度模型格式，建议修改为更通用的格式：

  thickness vs vp rho Qs Qp