SOD 3.0. Standing Order for Data A product of the University of South Carolina pro bono publico

Size: px

Start display at page:

Download "SOD 3.0. Standing Order for Data A product of the University of South Carolina pro bono publico"

Lilian Simmons
5 years ago
Views:

1 SOD 3.0 Standing Order for Data A product of the University of South Carolina pro bono publico

2 Waveform arm

3 Command Line SOD find_events find_stations find_seismograms find_responses

4 ga109:~ crotwell$ find_events --help Usage: find_events [-vrh] [-p <props>] [-S <server>] [-R <box>] [-d <donut>] [-b <begin>] [-e <end>] [-m <magnitude>] [-t <types>] [-D <depth>] -o <output> [-c <catalogs>] [--seis-regions <seismicregions>] [--geo-regions <geographicregions>] -v/--version Print SOD's version and exit -r/--recipe Print the created recipe to stdout instead of running it -h/--help Print this message. -p/--props <props> Use an additional props file --allow-duplicates Without this very similar events are rejected -S/--server <server> The event server to use. -R/--box-area <box> Event constraining box as west/east/south/north -d/--donut <donut> Event constraining donut as lat/lon/minradius/maxradius -b/--begin <begin> The earliest time for an accepted event -e/--end <end> The latest time for an accepted event -m/--magnitude <magnitude> The range of acceptable magnitudes -t/--types t1,t2,...,tn The types of magnitudes to retrieve. -D/--depth <depth> The range of acceptable depths in kilometers -o/--output <output> The format for output to standard out. -c/--catalogs c1,c2,...,cn A comma separated list of catalogs to search --seis-regions 1,2,...,n A comma separated list of seismic regions --geo-regions 1,2,...,n A comma separated list of geographic regions

5 ga109:~ crotwell$ find_events -b _307_00_56_56_ M _307_02_51_27_ M _307_03_01_50_ M _307_05_11_34_ M _307_09_09_54_ M _307_10_45_22_ M _307_11_04_00_ M _307_11_18_16_ M _307_12_09_45_ M ga109:~ crotwell$

6 ga109:~ crotwell$ find_events -b m _307_11_18_16_ M ga109:~ crotwell$

7 ga109:~/documents/tawg crotwell$ pscoast -JM8 -R-180/180/-70/70 -G200 -K -B30 > eventmap.ps ; find_events -b psxy - JM -R -Sc0.1 -G255/0/0 -O >> eventmap.ps ; pstopdf eventmap.ps ; open eventmap.pdf ga109:~/documents/tawg crotwell$

8 Stations ga109:~ crotwell$ find_stations --help Usage: find_stations [-vrh] [-p <props>] [-S <server>] [-R <box>] [-d <donut>] [-s <stations>] [-n <networks>] [-b <begin>] [-e <end>] - o <output> -v/--version Print SOD's version and exit -r/--recipe Print the created recipe to stdout instead of running it -h/--help Print this message. -p/--props <props> Use an additional props file -S/--server <server> The network server to use. -R/--box-area <box> A station constraining box as west/east/south/north -d/--donut <donut> A donut as lat/lon/minradius/maxradius -s/--stations s1,s2,...,sn The codes of stations to retrieve -n/--networks n1,n2,...,nn The codes of networks to retrieve -b/--begin <begin> The earliest time to overlap an accepted station -e/--end <end> The latest time to overlap an accepted station -o/--output <output> The format for output to standard out.

9 ga109:~ crotwell$ find_stations -n TA C Y22C V04C HAST I05A D03A A04A Y12C U05C U04C T06C R04C LAVA Q04C S04C HELL...

10 ga109:~ crotwell$ pscoast -JM8 -R-180/180/-70/70 -G200 -K -B30 > map.ps ; find_events -b psxy -JM -R -Sc0.1 -G255/0/0 -O -K >> map.ps ; find_stations -n TA psxy -JM -R -St0.1i -H1 -G0/0/255 -O >> map.ps ; pstopdf map.ps ; open map.pdf

11 Seismograms ga109:~ crotwell$ find_seismograms --help Usage: find_seismograms [-vrh] [-p <props>] [-S <server>] [-B <begin>] [-E <end>] [-l <sites>] [-c <channels>] [-f <seismogramfilename>] [--mark-phases <mark-phases>] -o <output> -v/--version Print SOD's version and exit -r/--recipe Print the created recipe to stdout instead of running it -h/--help Print this message. -p/--props <props> Use an additional props file -S/--server <server> Set the seismogram server to use for this search -B/--begin <begin> Phase name and offset for the seismogram's begin -E/--end <end> Phase name and offset for the seismogram's end -l/--sites l1,l2,...,ln The codes of sites(location codes) to retrieve -c/--channels c1,c2,...,cn The codes of channels to retrieve -f/--filename <seismogramfilename> Filename template for seismograms --mark-phases 1,2,...,n Phase arrival times to record in the SAC t headers -o/--output <output> The format for output to standard out.

12 ga109:~/documents/tawg crotwell$ find_events -b m 6 -r find_stations -n TA -s KMSC -r find_seismograms -c BHZ,BHE,BHN _307_11_18_16_ M _307_11_18_17_ M KMSC Got 1 seismograms for TA.KMSC..BHZ for eq on 2010/11/03 11:18:16 GMT Got 1 seismograms for TA.KMSC..BHN for eq on 2010/11/03 11:18:16 GMT Got 1 seismograms for TA.KMSC..BHE for eq on 2010/11/03 11:18:16 GMT Got 1 seismograms for TA.KMSC..BHZ for eq on 2010/11/03 11:18:17 GMT Got 1 seismograms for TA.KMSC..BHN for eq on 2010/11/03 11:18:17 GMT Got 1 seismograms for TA.KMSC..BHE for eq on 2010/11/03 11:18:17 GMT

13 -r to generate recipe ga109:~/documents/tawg crotwell$ find_events -b m 6 -r find_stations -n TA -s KMSC -r find_seismograms -c BHZ,BHE,BHN -r

14 Built In Processing rmean, rtrend, taper butterworth filtering instrument decon sac, mseed output decimate integrate, differentiate mul, div

15 Scripting In SOD JSR many languages (python, groovy, ruby, awk) Embedded in recipe Easy(ier) customization

16 <seismogramscript type="jython"> from jarray import array from bag import rmean, taper, rtrend, stat from edu.iris.fissures.seismogramdc import LocalSeismogramImpl print "Hello Sod: %s"%(channel,) outseis = [] for s in seismograms: s = rmean(s) s = taper(rtrend(s)) print stat(s) print outseis.append(s) out = array(outseis, LocalSeismogramImpl) result = util.waveformresult(true, out) </seismogramscript>

18 Thank You

Lesson 1-2: The SAC data format

Lesson 1-2: The SAC data format i. Philosophy and structure ii. Converting to, within, and from SAC format Lesson 1-3: SAC processing philosophy SAC has a complete set of reliable and welltested processing