RSEIS (version 2.6-0)

DISPLACE.SEISN: Displacement seismogram

Description

Removes seismic instrument response and integrates to displacement.

Usage

DISPLACE.SEISN(TH, sel = 1:length(TH$JSTR), inst = 1, Kal = Kal, FILT = list(ON = FALSE, fl = 1/30, fh = 7, type = "HP", proto = "BU"))

Value

• Same as input list with new traces representing displacement versus velocity

Details

Instrument responses are lists of poles and zeros for each instrument defined.

See Also

VELOCITY.SEISN, deconinst

Examples

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

DH =  DISPLACE.SEISN(TH, sel = which(inst>0) , inst = inst, Kal = Kal, FILT = list(ON = FALSE, fl = 1/200, fh = 7, type = "BP", proto = "BU"))

  SOUT =  swig(DH, PADDLAB=c("CENTER", "fspread", "HALF", "PREV"), filters=thefilts  )