if (FALSE) {
library(RSEIS)
source('ReadInstr.R')
S <- JSAC.seis('2010.012.21.54.46.6195.IU.OTAV.00.BHZ.R.SAC', BIGLONG = FALSE)
GH <- prepSEIS(S)
swig(GH)
RESP <-
ReadInstr('SAC_PZs_IU_OTAV_BHZ_00_2009.091.00.00.00.0000_2010.136.22.12.60.99999')
GHV <- VELOCITY.SEISN(GH, Kal = RESP, FILT = list(ON = TRUE, fl = 1/160,
fh = 7, type = "HP", proto = "BU"))
swig(GHV)
GHD <- DISPLACE.SEISN(GH, Kal = RESP, FILT = list(ON = TRUE, fl = 1/160,
fh = 7, type = "HP", proto = "BU"))
swig(GHD)
}
Run the code above in your browser using DataLab