brune.doom {RSEIS}

Brune Modeling

Description

Modeling the Brune spectrum with Graphical Diagnostics

Usage

brune.doom(amp, dt = 1, f1 = 0.01, f2 = 15, PLOTB = FALSE, tit = "")

Arguments

amp

signal

dt

deltaT

f1

low frequency for modeling

f2

high frequency for modeling

PLOTB

logical, TRUE=show diagnostic plots

tit

title for plot

Value

List:

SUCCESS

(0,1) for success or failure of modeling

WARN

flag = "OK"

tstar0

tstar0

gamma

gamma

omega0

omega0

fc

fc

alpha

alpha

chisqrd

chi-squared misfit over region of fitting

Author(s)

Jonathan M. Lees<jonathan.lees.edu>

References

Lees, J. M. and G. T. Lindley (1994): Three-dimensional Attenuation Tomography at Loma Prieta:Inverting t* for Q, J. Geophys. Res., 99(B4), 6843-6863.

Examples

data(CE1)
plot(CE1$x, CE1$y, type='l')
Xamp = CE1$y[CE1$x>5.443754 & CE1$x<5.615951]

BF = brune.doom( Xamp, CE1$dt ,f1=.5, f2=12 ,  PLOTB = TRUE)