powered by
Set up matrix for hypocenter inversion
XYSETUP(STAS, init, vel)
station information from pickfile
initial event location
list, velocity
matrix
This sets up the matrix used for nonlinear inversion. The code does not include information on the weighting. Station corrections are included.
The STAS are an internal component of the pickfile.
setPROJ, GLOB.XY,NLSlocate
# NOT RUN { ## start with the location of the closest station XY = XYSETUP(STAS, c(STAS$lat[w1],STAS$lon[w1], initz, STAS$sec[w1]-t0a ) , vel ) # } # NOT RUN { # }
Run the code above in your browser using DataLab