RSEIS (version 3.5-2)

longfft: Long FFT Spectrogram

Description

Creates hourly spectrograms, either alternating seismic and infrasound data or sequences of one component.

Usage

longfft(DB, DAYS = c(233, 234), HRS = 1:24, sta = "KR1", comp = c("V",
"I"), NPP = 6, CSCALE = FALSE, pal = rainbow(100), PS = FALSE, kind = 1,
Iendian = 1, BIGLONG = FALSE)
longreset(NPP, PS)
longpstart(NPP = 6, asta = "", acomp = "", theday = 1, hr = 0)

Arguments

DB

RSEIS Data base

DAYS

vector of Days to display

HRS

vector of hours to display

sta

stations to extract

comp

component to extract

NPP

Number of plot strips per page, default = 6

CSCALE

scaling

pal

palettes to use (given two will alternate these)

PS

logical, TRUE postscript output

kind

kind of seismic data, 1=SEGY, 2=SAC

Iendian

Endian-ness of binary data

BIGLONG

logical, TRUE=long is 8 bytes

asta

character, one station

acomp

character, one component

theday

one day

hr

one hour

Value

Graphical Side effects

Details

Extracts dats from the DB data base and plots strips of spectrograms for perusal.

longpstart, longreset are auxilliary codes used to set up the postscript files and initialize the plotting.

See Also

SPECT.drive

Examples

## Not run: 
# ###  output directory
# setwd("/home/lees/MARIO_TUNG/LONGFFT")
# 
# ###  2 color palettes
# p1 <- Gcols(plow=5, phi=0,  N=100, pal="topo.colors", mingray=0.8)
# p2 <- Gcols(plow=5, phi=0,  N=100, pal="rainbow", mingray=0.8)
# 
# ###   one station, vertical component
# longfft(DB, DAYS=c(148:156) , HRS=1:24 ,
# sta="BBIL", comp=c("BHZ") , NPP=6 , CSCALE=FALSE, 
# pal = list(p1=p1, p2=p2), PS=TRUE ,  kind = 2,
# Iendian=1, BIGLONG=FALSE )
# 
# ###   one station, infrasound
# longfft(DB, DAYS=c(148:156) , HRS=1:24 , sta="BBIL",
# comp=c("BDF") , NPP=6 , CSCALE=FALSE, 
# pal = list(p1=p2, p2=p1), PS=TRUE ,  kind = 2,
# Iendian=1, BIGLONG=FALSE )
# 
# ###  one station, both vertical and infrasound, alternating
# longfft(DB, DAYS=c(148:156) , HRS=1:24 , sta="BBIL",
# comp=c("BHZ", "BDF") , NPP=6 , CSCALE=FALSE, 
# pal = list(p1=p2, p2=p1), PS=TRUE ,  kind = 2,
# Iendian=1, BIGLONG=FALSE )
# 
# 
# 
# 
# 
# ## End(Not run)