schuster {GRTo} | R Documentation |
Plot of a phasor walkout for the Schuster's test
Description
This function plots a phasor walkout for the Schuster's test.
Usage
schuster(finame, title = "Schuster's diagram", color = c("black","red","orange",
"green","cyan","navy"), hd = FALSE, colidye = NULL, colidmo = NULL, colidda = NULL,
colidho = 1, colidmi = 2, colidma = 3, colidz = NULL, utccor = 0, dayt1 = NULL,
dayt2 = NULL, ysel1 = NULL, ysel2 = NULL, mosel1 = NULL, mosel2 = NULL,
dasel1 = NULL, dasel2 = NULL, magsel1 = NULL, magsel2 = NULL, zsel1 = NULL,
zsel2 = NULL, weekday1 = c("mo","tu","we","th","fr","sa","su"),
weekday2 = c("mo","tu","we","th","fr","sa","su"))
Arguments
finame |
name of the file to be loaded |
title |
main title that will appear in the phasor walkout plot |
color |
color of the phasor walkout (black/red/orange/green/cyan/navyblue). The first letter is enough to discriminate the color name |
hd |
logical. Whether |
colidye |
field number for the year values in the loaded file |
colidmo |
field number for the month values in the loaded file |
colidda |
field number for the day values in the loaded file |
colidho |
field number for the hour values in the loaded file |
colidmi |
field number for the minute values in the loaded file |
colidma |
field number for the magnitude values in the loaded file |
colidz |
field number for the depth values in the loaded file |
utccor |
correction for local time. This value and values of the start and end hours of nighttime are used to filter the daytime seismic noise from the walkout plot |
dayt1 |
the start hour of nighttime in local time |
dayt2 |
the end hour of nighttime in local time |
ysel1 |
minimum year value for selection (Default is NULL, e. g. no selection) |
ysel2 |
maximum year value for selection (Default is NULL, e. g. no selection) |
mosel1 |
minimum month value for selection (Default is NULL, e. g. no selection) |
mosel2 |
maximum month value for selection (Default is NULL, e. g. no selection) |
dasel1 |
minimum day value for selection (Default is NULL, e. g. no selection) |
dasel2 |
maximum day value for selection (Default is NULL, e. g. no selection) |
magsel1 |
minimum magnitude value for selection (Default is NULL, e.g. no selection) |
magsel2 |
maximum magnitude value for selection (Default is NULL, e.g. no selection) |
zsel1 |
minimum depth value for selection (Default is NULL, e. g. no selection) |
zsel2 |
maximum depth value for selection (Default is NULL, e. g. no selection) |
weekday1 |
first day in the week for selection ("mo" for "Monday", "tu" for "Tuesday", "we" for "Wednesday", "th" for "Thursday", "fr" for "Friday", "sa" for "Saturday", "su" for "Sunday") |
weekday2 |
last day in the week for selection ("mo" for "Monday", "tu" for "Tuesday", "we" for "Wednesday", "th" for "Thursday", "fr" for "Friday", "sa" for "Saturday", "su" for "Sunday") |
Details
This function reads earthquake times (hour and minutes values are required) in an input file and plots a phasor walkout into a file with name file
\_schu.png (png format file).
Selections can be performed based on ranges in hour of the day, year, month, day, magnitude, depth and day of the week (in the same week). The correction for local time is such that GMT time + correction (e. g. utccor
) = local time.
Author(s)
Daniel Amorese \& Paul A. Rydelek
References
Rydelek, P. A. and Hass, L. (1994) On Estimating the Amount of Blasts in Seismic Catalogs with Schuster's Method Bulletin of the Seismological Society of America, Vol. 84, No. 4, pp. 1256-1259.
Zurn, W. and Rydelek, P. A. (1996) Revisiting the phasor-walkout method for detailed investigation of harmonic signals in time series Surveys in Geophysics, Vol. 15, No. 4, pp. 409-431.
Examples
schuster(finame=system.file("extdata","IDYLLdeep.data.txt",package="GRTo"),
title = "Schuster's diagram", hd = FALSE, colidye = 1, color="n", colidmo = 2,
colidda = 3, colidho = 4, colidmi = 5, colidma = 15, colidz = 9, utccor = -9,
dayt1 = NULL, dayt2 = NULL, ysel1 = 1983, ysel2 = 1990, mosel1 = NULL,
mosel2 = NULL, dasel1 = NULL, dasel2 = NULL, magsel1 = NULL, magsel2 = NULL,
zsel1 = NULL, zsel2 = NULL, weekday1 = NULL, weekday2 = NULL)