rdseed

Group:

Authors:

Igor B. Morozov, 4/2018

Load IRIS SEED files files into traces and/or trace headers

This tool has two input modes:

1. STILL NOT IMPLEMENTED: If there are no input traces but lists DATA are present, the tool obtains SEED or MiniSEED file names from them and loads these files as data traces. With each trace record, the tool then performs all operations specified in the second operation mode.
2. For each input trace, RDSEED looks up the specified ASCII file containing station, event, or instrument response information.

   Normally, this file should be produced by IRIS program rdseed applied to a dataless SEED file. Optionally, a call to rdseed can be performed just before loading values by using list DATALES.

   RDSEED scans through the outputs of rdseed and locates the time interval corresponding to the current trace. From the identified segment, RDSEED reads in the values of station and channel parameters into trace headers specified in the list below.

Trace headers loaded by RDSEED from rdseed output files:

The trace headers used for identifying the station and channel are:

1. "stname" - Station name
2. "channel" - Channel name
3. "stlat" - Station latitude
4. "stlon" - Station longitude
5. "stelev" - Station elevation
6. "year" - GMT year of reference time at which the response is evaluated (usually time start of the record)
7. "jday" - GMT year of reference time above
8. "hour" - GMT hour of reference time above
9. "min" - GMT minute of reference time above
10. "sec" - GMT second of reference time above
11. "msec" - GMT millisecond of reference time above
12. "network" - Network name (not used for response lookup)
13. "datrd" - Data recorder name (not used for response lookup)
14. "instr" - Instrument name (not used for response lookup)

Optionally (not used now) "nz" trace headers from SAC can be used for times:

1. "nzyear" - GMT year of reference time at which the response is evaluated (usually time start of the record)
2. "nzjday" - GMT year of reference time above
3. "nzhour" - GMT hour of reference time above
4. "nzmin" - GMT minute of reference time above
5. "nzsec" - GMT second of reference time above
6. "nzmsec" - GMT millisecond of reference time above

After reading the RDSEED file, the appropriate segment of its parameters will be transfered into the following headers (created by this tool):

1. "resp.a" - A0 normalization factor
2. "resp.af" - A0 normalization frequency
3. "resp.c" (array of floats) - Matrix of coefficients, gains, gain frequencies, zeros and poles, or numerators and denominators for all stages
4. "resp.n" - Number of stages (factors) in response expression
5. "resp.s" - sensitivity
6. "resp.sf" - frequency of sensitivity
7. "resp.p" (array of floats) - The first field is the number of poles (NP) and the rest are NP pairs of real and imaginary parts of the frequencies of poles.
8. "resp.t" (character string) - string consisting of characters of the transfer-function codes for the first two stages
9. "resp.ui" (array of 8-character strings) - Codes of input units for each stage
10. "resp.uo" (array of 8-character strings) - Codes of output units for each stage

After the response parameters are stored in the above trace headers, response corrections can be applied to the records or output as trace records by using option 'resp' in tool FILTER.

See the DMC library of Nominal Seismic Instrument Responses

Sample of IRIS rdseed output:

(This file is an example, compilation of multiple records and does not represent any single station)

B050F03 Station code: ABH1 B050F04 Latitude: 30.599970 B050F05 Longitude: 50.253150 B050F06 Elevation: 346.000000 B050F09 Name: Behbahan B050F10 Owner Code Lookup: 1 IRSC Network B050F11 32-bit word order: 3210 B050F12 16-bit word order: 10 B050F13 Starting date: 2015,007,00:00:00.0000 B050F14 Ending date: (null) B050F15 Update flag: N B050F16 Network Code: IR

B050F03 Station: FURI B050F16 Network: IU B052F03 Location: 00 B052F04 Channel: BHE

B052F04 Channel: BHN B052F03 Location: B052F05 Subchannel: 0 B052F06 Instrument lookup: 2 Trillium 120P, 120 s, 1201 V/m/s-Taurus, 40 Vpp (g B052F07 Comment: X B052F08 Signal units lookup: 1 M/S - Velocity in Meters Per Second B052F09 Calibration units lookup: 2 V - EMF in Volts B052F10 Latitude: 30.599970 B052F11 Longitude: 50.253150 B052F12 Elevation: 346.000000 B052F13 Local depth: 0.000000 B052F14 Azimuth: 0.000000 B052F15 Dip: 0.000000 B052F16 Format lookup: 1 Format Information Follows B052F22 Start date: 1999,111,00:00:00 B052F23 End date: No Ending Time

B053F03 Transfer function type: A [Laplace Transform (Rad/sec)] B053F04 Stage sequence number: 1 B053F05 Response in units lookup: M/S - Velocity in Meters Per Second B053F06 Response out units lookup: V - Volts B053F07 A0 normalization factor: 3948.58 B053F08 Normalization frequency: 0.02 B053F09 Number of zeroes: 2 B053F14 Number of poles: 4

B053F10-13 0 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 B053F10-13 1 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

B053F15-18 0 -1.234000E-02 1.234000E-02 0.000000E+00 0.000000E+00 B053F15-18 1 -1.234000E-02 -1.234000E-02 0.000000E+00 0.000000E+00 B053F15-18 2 -3.918000E+01 4.912000E+01 0.000000E+00 0.000000E+00 B053F15-18 3 -3.918000E+01 -4.912000E+01 0.000000E+00 0.000000E+00

B058F03 Stage sequence number: 1 B058F04 Gain: 2.296000E+03 B058F05 Frequency of gain: 2.000000E-02 HZ B058F06 Number of calibrations: 0

B054F03 Transfer function type: D B054F04 Stage sequence number: 2 B054F05 Response in units lookup: V - Volts B054F06 Response out units lookup: COUNTS - Digital Counts B054F07 Number of numerators: 0 B054F10 Number of denominators: 0

B057F03 Stage sequence number: 2 B057F04 Input sample rate: 5.120000E+03 B057F05 Decimation factor: 1 B057F06 Decimation offset: 0 B057F07 Estimated delay (seconds): 0.000000E+00 B057F08 Correction applied (seconds): 0.000000E+00

B058F03 Stage sequence number: 2 B058F04 Gain: 4.194300E+05 B058F05 Frequency of gain: 0.000000E+00 HZ B058F06 Number of calibrations: 0

The end of RDSEED file is like this:

B058F03 Stage sequence number: 0 B058F04 Sensitivity: 9.630000E+08 B058F05 Frequency of sensitivity: 2.000000E-02 HZ B058F06 Number of calibrations: 0

RDASCII, FILTER

Module call format:

1 group of 1 REPEAT parameter per line: