markseis24: Mark 24 hour seismic display

In RSEIS: Seismic Time Series Analysis Tools

Mark 24 hour seismic display

Description

Mark a 24 hour seismic display

Usage

markseis24(pjj, pix = list(yr = 2009, jd = 1, hr = 0, mi = 0, sec = 0,
dur = 0), col = "red", LEGON = 3, BARON = TRUE, ARROWS = TRUE, lwd=1)

Arguments

pjj	Output information from plotseis24 (x,y, yr, jd)
pix	list: date list consisting of: yr, jd, hr, mi, sec, dur)
col	Color, specified as color index, character string or rgb
LEGON	plotting flag for legs: 0=no legs, 1=left leg, 2=right leg, 3=both legs(def ault)
BARON	logical:plotting flag for bar
ARROWS	logical: plot arrows FALSE=no arrows
lwd	numeric, graphical parameter, line width

Details

the LEGON parameter controls the small marks at the ends: Either left(1) right(2) both(3) or no legs(0) are plotted. window bars should wrap around the ends of the hour to the next hour below. The durations of the windows are supplied in seconds. If no duration is supplied, it is set to 0. If one duration is supplied it is copied to all other windows.

Value

Graphical Side effects

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

winmark, getseis24, plotseis24

Examples

## Not run: 
##  the database DB was constructed with makeDB

usta = "CAL"
acomp = "V"

iday = 2

###   get one day's worth of data
JJ <- getseis24(DB, 2009, iday, "CAL", "V", kind = 1)

###  plot day
pjj <- plotseis24(JJ, dy=1/18, FIX=24, SCALE=0,
FILT <- list(ON=FALSE, fl=0.5 , fh=20.0, type="BP", proto="BU"),
RCOLS=c(rgb(0.2, .2, 1), rgb(.2, .2, .2))  )
title(paste("CAL", "V"))

###   set up pix
WINS2 <- win15
Apix <- WINS2$hr[seq(from=1, to=length(WINS2$hr), by=2) ]
dur <- (WINS2$hr[seq(from=2, to=length(WINS2$hr), by=2) ]-Apix)*3600

## dur <- rep(0, times=length(Apix))

##  mark the 24 hour plot
markseis24(pjj, pix =list(yr=rep(pjj$yr, length(Apix)),
 jd=rep(pjj$jd, length(Apix)) , hr=Apix, mi=rep(0, length(Apix)),
sec=rep(0, length(Apix)), dur=dur), col='red', ARROWS=FALSE )





## End(Not run)

RSEIS documentation built on Aug. 19, 2023, 5:07 p.m.