get3Drayblox {RTOMO} | R Documentation |
Get blocks from 3D ray path
get3Drayblox(XNOD, YNOD, ZNOD, xo, yo, ztop, slowness = NULL)
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 |
ix |
x index |
iy |
y index |
iz |
layer index |
r |
length in each block |
tt |
travel time along whole raypath |
Jonathan M. Lees<jonathan.lees@unc.edu>
get2Drayblox
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)