Kal = PreSet.Instr()
Majorevent = 250+8/24+2250/(24*3600)
at1 = Majorevent - 1/24
at2 = Majorevent + 1/24
TH = Mine.seis(at1, at2, DB99, NULL, NULL , kind = 1, Iendian=1, BIGLONG=FALSE, CHOP=TRUE, verbose=FALSE)
inst = rep(0, length(TH$STNS))
inst[TH$COMPSinst[TH$STNS=="krm2" & TH$COMPSinst[TH$STNS=="krm0" & TH$COMPS
VH = VELOCITY.SEISN(TH, sel = 1:length(TH$JSTR), inst = 1, Kal = Kal, FILT = list(ON = FALSE, fl = 1/30, fh = 7, type = "HP", proto = "BU"))
Run the code above in your browser using DataLab