P_wave_gen {ezECM} | R Documentation |
Generation of noisy p-wave arrival times
Description
Similar utility to time_fn, however multiple seismometer locations can be provided simultaneously and normally distributed noise is added to the arrival time.
Usage
P_wave_gen(
Si = NULL,
S0 = NULL,
Sig = NULL,
neg.obs = TRUE,
eps = sqrt(.Machine$double.eps)
)
Arguments
Si |
Numeric matrix providing seismometer locations. Must contain 3 columns corresponding to (X,Y) corrdinates and depth. |
S0 |
Numeric 3 element vector stipulating the location of an event, elements correspond to (X, Y, Z) |
Sig |
Numeric vector, or diagonal matrix, providing the variance in observed arrival times at each seismometer. |
neg.obs |
Logical indicating whether to allow negative observations of time (eg. the observed time of p-wave arrival is before the true time for the event). |
eps |
Numeric. If |
Value
Numeric vector of observation times that correspond to the rows of Si
Examples
pwave.obs <- P_wave_gen(Si = c(100,200,3), S0 = c(400, 500, 4), Sig = 0.05)