RTOMO (version 1.1-3)

get3Drayblox: Get blocks from 3D ray path

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

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)