view.seis {RSEIS}R Documentation

Veiw seismic data window

Description

Veiw seismic data (segy) window on an hourly basis.

Usage

view.seis(aday, ihour, inkhour, SAVEFILE, days, DB, usta, acomp,  STDLAB =c("QUIT",  "NEXT", "PREV", "HALF"), TZ=NULL)

Arguments

aday index of which day to use in vector days
ihour hour to start
inkhour increment in hours for viewing panel
SAVEFILE file to save window picks in
days vector of days to select from
DB data base list of file names and start-times and durations
usta stations to select
acomp compnents to select
STDLAB vector of buttons, DEFAULT = c("QUIT", "NEXT", "PREV", "HALF", "WPIX", "zoom out", "refresh", "restore", "SPEC", "SGRAM" ,"WLET", "FILT", "Pinfo", "WINFO")
TZ Number of hours to add to GMT to get local time

Details

The program view.seis assumes the data is stored in files accessable by the user and that the DB list has been scanned in and parsed.

Value

Graphical side effects and save.wpix stores appended picks.

Note

On LINUX systems I wrote these (non-R) programs to set up the data base for segy data:FLS.prl, segydatabase. To get these contact me directly. TZ is (-6) for Guatemala.

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

See Also

PICK.GEN, save.wpix

Examples


## Not run: 
#######  read in data base from files:
DB  = scan(file="/home/lees/Site/Santiaguito/SG09/segyDB", list(fn="",
yr=0, jd=0, hr=0, mi=0, sec=0, dur=0))

##  this file has first line:
#/home/lees/Site/Santiaguito/SG09/R001.02/09.001.11.57.09.SMI.J 2009 1 11 57 9.245000 3600.000000
#c...|....1....|....2....|....3....|....4....|....5....|....6....|....7.$

##  note: the file name includes the station SMI and the component J

##########  extract names of stations and components
das1  = substr(DB$fn, 58, 60)
comp1 = substr(DB$fn, 62,62))

###########  set up times (start and end) in julian days for selection
 DB$t1 = DB$jd + DB$hr/24 + DB$mi/(24*60) + DB$sec/(24*3600)
DB$t2 = DB$t1 + DB$dur/(24*3600)

#######   get unique stations from list:
usta = unique(das1)

###########   choose the desired components 
acomp = c('I','V')

aday = 3
SAVEFILE = 'wpix.out'
ihour = 15
inkhour = .5

view.seis(aday, ihour, inkhour, SAVEFILE, days, DB, usta, acomp, TZ=(-6))


## End(Not run)



[Package RSEIS version 2.1-4 Index]