Thread: get header data

As Art suggests, saclst will hopefully provide
the header list capabilities you need.
In the short term you can hopefully get by using listheader within a script and
redirect the output to the file you want.

Something like the following script could be invoked as
sac_lh.csh "your_filename"
Note that the blank lines are needed to get all the headers and you'll
still have to edit the results if you want a particular format or to get rid of
extraneous information and blank lines.


#!/bin/csh -xvf
# sac_lh.csh: script to list sac header listing to standard out
sac<<eof
r $argv[1]
lh inclusive on




q
eof
# end of sac_lh.csh

Good morning,

I'd like to get some information in the header but the subroutines getfhv
doesn't seem to work for depmin, depmax and dist

in a fortran code:

call GETFHV ('DEPMIN', depmin, nerr)

I got this error message:
ERROR: Undefined header field value. DEPMIN

but it works for:

call GETNHV ('NZYEAR', nyear, nerr)

is there a bug? or an other way to extract data from the header?

Best regards

HÈloÔse
CETE Nice, LGIT Grenoble



_______________________________________________
sac-help mailing list
sac-help<at>iris.washington.edu
http://www.iris.washington.edu/mailman/listinfo/sac-help

--

Peter Goldstein, Ph.D. (925) 423-1231 (office)
L-103, PO Box 808 (925) 422-5844 (fax)
Livermore, CA 94551 peterg<at>llnl.gov (email)
web page: http://www.llnl.gov/sac

Heloise

these values all need to be calculated by sac so when the sac file is
generated externally, e.g. by rdseed they are generally left undefined.

Try reading into sac and writing again, i.e.

sac<<eof
r file
w over
q
eof

Of course, for dist evla,evlo,stla,stlo all need to be set and lcalda
must be true.

Frederik


Heloise.Cadet<at>ujf-grenoble.fr wrote:

Good morning,

I'd like to get some information in the header but the subroutines getfhv
doesn't seem to work for depmin, depmax and dist

in a fortran code:

call GETFHV ('DEPMIN', depmin, nerr)

I got this error message:
ERROR: Undefined header field value. DEPMIN

but it works for:

call GETNHV ('NZYEAR', nyear, nerr)

is there a bug? or an other way to extract data from the header?

Best regards

H?lo?se
CETE Nice, LGIT Grenoble



_______________________________________________
sac-help mailing list
sac-help<at>iris.washington.edu
http://www.iris.washington.edu/mailman/listinfo/sac-help

--
Frederik Tilmann
Bullard Laboratories Tel. +44 1223 765545
Department of Earth Sciences Fax. +44 1223 360779
University of Cambridge email: tilmann<at>esc.cam.ac.uk
Madingley Road
Cambridge CB3 0EZ
UK