hongjianfang / dsurftomo Goto Github PK

Hi @HongjianFang
I am using DSurfTomo (v1.3) to invert 112422 Rayleigh phase velocity measurements. First I encountered this error:
Index '201' of dimension 1 of array 'rthk' above upper bound of 200Index '201' of dimension 1 of array 'rthk' above upper bound of 200
So I increased the 200 default values of the NL and NLAY parameters to 300 in surfdisp96.f and CalSurfG.f90.
Now I get this error while the code tries to compute sensitivity matrix:

DSurfTomo(55526,0x114dc1dc0) malloc: Region cookie corrupted for region 0x7fdbdf600000 (value is bd8fbf3f)
DSurfTomo(55526,0x114dc1dc0) malloc: *** set a breakpoint in malloc_error_break to debug

Program received signal SIGABRT: Process abort signal.

Backtrace for this error:
#0 0x10818872c
#1 0x108187ad3
#2 0x7fff673c042c
Abort trap: 6

I also tried DSurfTomo (v2) and here is the error:

DSurfTomo(57162,0x10fc7cdc0) malloc: Incorrect checksum for freed object 0x7f876bd009b8: probably modified after being freed.
Corrupt value: 0xbcec818dbd19b28d
DSurfTomo(57162,0x10fc7cdc0) malloc: *** set a breakpoint in malloc_error_break to debug

Program received signal SIGABRT: Process abort signal.

Backtrace for this error:
#0 0x10b6ae72c
#1 0x10b6adad3
#2 0x7fff673c042c
Abort trap: 6

According to the manual, there should be a 'raypath.out' file after inversion.

"""
The raypath distribution of the final model (raypath.out) #
number of ray path segments
latitude longitude
"""

However, I didn't find this file in the outputs. I also checked DSurfTomo source codes and didn't find the codes to write ray path geometries either. So, how can I save and plot ray paths?

Hello, I got the following error when I ran the program：
[root@Z textXD]# ../src/DSurfTomo DSurfTomo.in

                          DSurfTomo (v1.3)

PLEASE contact Hongjain Fang ([email protected]) if you find any bug

model origin:latitude,longitue
28.2215 102.5475
grid spacing:latitude,longitue
0.0005 0.0005
model dimension:nx,ny,nz
10 10 9
Rayleigh wave phase velocity used,periods:(s)
0.0 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
At line 255 of file main.f90
Fortran runtime error: Index '23' of dimension 2 of array 'scxf' above upper bound of 22
About the period I set in DSurfTomo.in：0.04 0.05 0.06 0.07 0.08 0.09 0.10 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.19 0.20 0.21 0.22 0.23 0.24 0.25
I really hope to get your help, thank you very much.

How do you test the reliability of the inversion results of this program?. Do you have related methods or procedures? For example, the corresponding procedures to do the test board.

代码使用过程中出现如下问题：
error with allocate: integer iw
Program received signal SIGSEGV:Segmentation fault - invalid memory reference

Refer to previous S velocity model，I give a linear model as default. During the preliminary attempt, I give a wide velocity range(minimum velocity, maximum velocity), but the extreme value(the highest and lowest velocity) of inversion result tend to the boundary value of the velocity range, rather than reasonable value(i.e. while depth is 10km or 0km, the S velocity is about to 4km/s or 2.5km/s). Besides, with each iteration, the inversion result shows characteristic that high velocity is more higher, low velocity is more lower, which seems to be trapped in local minima. What might cause this situation? Except modifying the priori velocity range, is it necessary to modify other parameters, or check my dispersion curves?

尊敬的方博士您好，非常荣幸可以使用您的DSurfTomo程序。目前希望可以使用到DSurfTomo中生成各周期台站间射线路径的生成文件，在地震学学术年会上听闻早期版本可以实现，特向您请求下载DSurfTomo可以生成Raypath.out版本的文件包。
还想请问您Rayleigh波不同周期的深度敏感核是用什么程序做的呢？如果方便的话您可以做个分享吗？

I used the command as mentioned in the manual "./src/DSurfTomo DSurfTomo.in" , somehow it shows this " DSurfTomo (v2.0)
For bug report, PLEASE contact Hongjain Fang ([email protected])

STOP unable to open the inputfile"
would be grateful for fixing this issue. Thank you