get3Drayblox: Get blocks from 3D ray path

In RTOMO: Visualization for Seismic Tomography

Description

Get blocks from 3D ray path

Usage

get3Drayblox(XNOD, YNOD, ZNOD, xo, yo, ztop, slowness = NULL)

Arguments

XNOD	x-coordinates along raypath
YNOD	y-coordinates along raypath
ZNOD	z-coordinates along raypath
xo	x block divisions
yo	y block divisions
ztop	vector, topsof layers
slowness	vector, Slowness model

Value

ix	x index
iy	y index
iz	layer index
r	length in each block
tt	travel time along whole raypath

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

get2Drayblox

Examples

xo = seq(from=0.5, by=1, length=100)
yo = seq(from=0.5, by=1, length=100)
 ztop = seq(from=0, to=24, by=4)




Xp = c(11.5, 70.2)
Yp = c(3.5, 50.2)
Zp = c(18.4, 0.0)

dee = sqrt(  (Xp[2]-Xp[1])^2 + (Yp[2]-Yp[1])^2 +  (Zp[2]-Zp[1])^2 )
deexy = sqrt(  (Xp[2]-Xp[1])^2 + (Yp[2]-Yp[1])^2  )

fi  = findInterval(Zp, ztop)

ZNOD = c(Zp[1], ztop[fi[1]:fi[2]])


alpha = asin(deexy/dee)

RN = deexy-ZNOD*tan(alpha)

XNOD = Xp[1]+RN*(Xp[2]-Xp[1])/deexy
YNOD = Yp[1]+RN*(Yp[2]-Yp[1])/deexy

IYZ = get3Drayblox(XNOD, YNOD, ZNOD, xo, yo, ztop, slowness = NULL)

RTOMO documentation built on May 2, 2019, 3:35 p.m.