P_wave_gen: 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)
)
Value
Numeric vector of observation times that correspond to the rows of Si
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 neg.obs = FALSE sets of observations are redrawn until all \(t_i - t_0 \leq\) eps.