Nothing
`many.time1D` <-
function(indelta, inhpz, instaz, inlay , ztop , vel)
{
### One Dimensional travel times with low velocity layers
##### tt1d travel.time1D calculate the travel time, dtdr, dtdz and
##### takeoff angle for the given vel, dist(km), hpz and stz
###### input:
###### indelta = distance from source to reciever
###### inhpz = depth of source
###### instaz = elevation of reciever
###### inlay =number of layers
###### ztop = tops of the vel model
###### vel = velocity in the layers
## output: dtdr=derivative with respect to r dtdz=derivative with respect to z
### angle=angle of takeoff measured from nadir (down vecotor) tt=travel time
N = length(indelta)
if(length(ztop)<inlay) { warning("ERROR: BAD velocity model"); return(NULL) }
if(length(vel)<inlay) { warning("ERROR: BAD velocity model"); return(NULL) }
if( any(indelta<0)) { warning("ERROR: BAD Distance in travel.time1D "); return(NULL) }
if(any(!is.numeric(ztop))){ warning("ERROR: BAD velocity model"); return(NULL) }
if(any(!is.numeric(vel))){ warning("ERROR: BAD velocity model"); return(NULL) }
dtdr=rep(0, N);
dtdz=rep(0, N);
angle=rep(0, N);
outt=rep(0, N);
if(length(instaz)==1) { instaz = rep( instaz, N) }
TTout = .C("CALL_DTTmany",PACKAGE = "RSEIS",N,
as.double(indelta), as.double(inhpz),as.double(instaz), as.integer(inlay), as.double(ztop) , as.double(vel),
as.double(dtdr), as.double(dtdz), as.double(angle), as.double(outt) )
dtdr=as.numeric(unlist(TTout[8]));
dtdz=as.numeric(unlist(TTout[9]));
angle=as.numeric(unlist(TTout[10]));
tt=as.numeric(unlist(TTout[11]));
return(list( dtdr=dtdr, dtdz=dtdz, angle=angle, tt=tt))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.