slimgroup / JUDI.jl

Julia Devito inversion.
https://slimgroup.github.io/JUDI.jl
MIT License
97 stars 31 forks source link

Source geometry is confusing #93

Closed ziyiyin97 closed 2 years ago

ziyiyin97 commented 2 years ago

MFE

using JUDI
using SegyIO
datapath = joinpath(dirname(pathof(JUDI)))*"/../data/"
nsrc = 2
block = segy_read(datapath*"unit_test_shot_records_$(nsrc).segy"; warn_user=false)
src_geometry = Geometry(block; key="source", segy_depth_key="SourceSurfaceElevation")
println(src_geometry.xloc)    

This is a part of the test and it outputs src_geometry.xloc as


julia> println(src_geometry.xloc) 
Float32[100.0, 1100.0]

I personally find it quite confusing because there are 2 sources but only here is 1 vector for xloc. Typically this means they are simultaneous but in this case they are not.