markseis24 {RSEIS} | R Documentation |
Mark a 24 hour seismic display
markseis24(pjj, pix = list(yr = 2009, jd = 1, hr = 0, mi = 0, sec = 0, dur = 0), col = "red", LEGON = 3, BARON = TRUE, ARROWS = TRUE)
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 |
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.
Graphical Side effects
Jonathan M. Lees<jonathan.lees@unc.edu>
winmark, getseis24, plotseis24
## 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)