get3Drayblox {RTOMO}R Documentation

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

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

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)




[Package RTOMO version 1.0-6 Index]