[SAC-HELP] suggestions on finding maximum value of spectrum

[Sheila Peacock]

Dear all,

These were suggested by my colleague.  I can take no credit for them.
Some of the SAC commands have been abbreviated.

Regards,
Sheila Peacock.

Sheila,

if they want to know the frequency at which the maximum amplitude occurs, then

r $1
rtr
taper
fft

* add any smoothing in here

keepam
markptp
lh t1


does the trick.


The maximum value of an amplitude spectrum can be simply obtained by

fft
keepam
lh depmax

Do they want to know the frequency when this occurs?

the value of the time series at a marker, say t1, can be found by something like this:

r $1
evaluate to x1 &1,t1 - ( &1,delta / 2 )
evaluate to x2 &1,t1 + ( &1,delta / 2 )
cut  %x1 %x2
r
lh depmax