IGeoS (SIA) USArray demos

This page contains demos briefly shown at the USArray data processing summer school at Northwestern University, Aug 05-08 2009.

Presentation : Process-oriented data handling

Manuscript explaining these demos

Breq_fast requests and SEED

This section shows how breq_fast requests can be formed from station and event files, with distances and channel numbers restricted in command line. The second job shows a way to execute multiple rdseed calls to produce many gzipped ASCII SAC files.

Shell script
# Task Job file Command Log file Output files Output images
1 Fill event table (withot SQL database - for use on MACs) events-nodb.job sia events-nodb.job - - -
2 Load a list of stations and save in a binary file (without SQL database) stations.job sia stations.job - - -
3 Make breq_fast requests from station and event tables request.job sia request.job 200 300 BHZ - Resulting request -
4 Convert multiple seed files to gzipped ASCII SAC rdseed.job sia rdseed.job - - -

Populating the database

These jobs initialize the SQL database and load events and all waveform indexes into it.

Shell script
# Task Job file Command Log file Output files Output images
1 Define (also clear) database tables make_tables.job sia make_tables - - -
2 Fill event table events.job sia events.job - - -
4 Loading and indexing gzipped SAC files read_sac.job sia read_sac.job 2005 'SAC/*/2005*.SAC_ASC.gz' - - -
5 Create station corridor (straight) make-corridor.job sia make-corridor.job 1 corridor1.txt 0.5 log - -
6 Create station corridor (angular) make-corridor.job sia make-corridor.job 2 corridor2.txt 0.5 log - -

Extracting event windows

These jobs time extract windows for selected arrivals from the database, perform demeaning and component rotation, and save them in files.

Shell script
# Task Job file Command Log file Output files Output images
1 Extract P-wave time gate, rotate components, save in file event-gather.job sia event-gather.job 80 P 20 1200 log - -
2 Extract P-wave time gates for all events in the database; plot station maps for all events event-gathers.job sia event-gathers.job P 20 1200 - - Maps of all events

Lists

Simple jobs for listing the contents of the database. The general approach to these listings is to generate a no-data trace for each database record, and pass the information of interest to the subsequent tools through trace headers.

Shell script
# Task Job file Command Log file Output files Output images
1 Sample arbitrary SQL queries to the database show-db.job sia show-db.job log - -
2 List stations list-stations.job sia list-stations.job log - -
3 List stations in corridor #1 list-corridor.job sia list-corridor.job 1 log - -
4 List events list-events.job sia list-events.job log - -
5 List all waveforms list-waveforms.job sia list-waveforms.job log - -
6 List waveforms from event #80 recorded in corridor #1 list-waveforms-corridor.job sia list-waveforms-corridor.job 80 1 log - -

Maps

These jobs create various maps by using the same approach: traces are generated and station (or event) coordinates are passed in their headers. GMT scripts are invoked to produce PostScript images.

Shell script
# Task Job file Command Log file Output files Output images
1 Map of all events map-events.job sia map-events.job log - map
2 Map of all stations map-stations.job sia sod.job show log - map
3 Map station corridor #1 (straight) map-corridor.job map-corridor.job 1 log - map
4 Map station corridor #2 (angular) map-corridor.job map-corridor.job 2 log - map
5 Map channels which recorded between certain dates map-channels.job map-channels.job '2005/04/22' '2006/04/22' 2008 log - map

Record sections with IASP91 travel times

Vertical-component record section plots are created from assembled event gathers. Shell script
# Task Job file Command Log file Output files Output images
1 All records of a single P-wave event vs. range, with P, S, and Lg travel times plot-event.job sia plot-event.job 80 P 3 PS log - plot
2 Plot time sections of al events plot-events.job sia plot-events.job P 3 - - All plots
3 A single P-wave event in corridor #2 vs. range, with P, S, and Lg travel times plot-event-corridor.job sia plot-event-corridor.job 2 P 1 3 log - plot
4 Time slice of obne event in "colored spots" style plot-event-tspots.job sia plot-event-tspots.job 7 P 3 150.0 5.0 log - plot

Encapsulated data requests

SOD-like data selections are performed by executing SQL database requests based on user's specifications of station names, regions of interest, event magnitudes, distances, asimuths, etc.
Shell script
# Task Job file Command Log file Output files Output images
1 List waveforms for SOD-like data retrieval sod-list.job sia sod-list.job show log networks list; origins list; stations list; waveforms list; map
2 Download "SOD" waveforms in a file sod.job sia sod.job show log networks list; origins list; stations list; waveforms list; map

Remote (web) operation

Jobs executing various requests to a remote IGeoS server. This server has an identical IGeoS installation capable of executing any jobs and posting results for download.

Shell script
# Task Job file Command Log file Output files Output images
1 Obtain remote IGeoS server info web-info-server.job sia web-info-server.job response - -
2 Obtain information about GMT palettes installed on seisweb.usask.ca web-info-item.job sia web-info-item.job palettes response - -
3 Obtain a list of 'logical drives' defined on seisweb.usask.ca web-info-item.job sia web-info-item.job ldrives response - -
4 Obtain usage examples for module 'sort' on seisweb.usask.ca web-info-module.job sia web-info-module.job sort response - -
5 Execute job 'remote-get-data' on seisweb.usask.ca and obtain results web-flow-send.job sia web-flow-send.job USArray remote-get-data BHZ 20 log list of channels; remote log; server program run report; Data file; map

Waveform synthetics

Three-component short-period waveforms in IASP91 model are computed by mode summation (the 'reflectivity', or propagator method, Fuchs and Muller, 1971), adapted and paralleleized.

Shell script
# Task Job file Command Log file Output files Output images
1 All phases in IASP91 model (Ricker wavelets) phase-synthetic.job sia phase-synthetic.job 20 log - PostScript plot
2 1D 'reflectivity' synthetic in IASP91 model to 870-km depth using 20 processors reflect-synthetic.job sia reflect-synthetic.job 20 log by host #0 - record section plot