Changeset 1108


Ignore:
Timestamp:
14.01.2016 09:09:23 (4 years ago)
Author:
klaus
Message:

command writeo; command plot_stations; updated examples and help text

Location:
SHX/trunk/SeismicHandler
Files:
2 added
4 edited

Legend:

Unmodified
Added
Removed
  • SHX/trunk/SeismicHandler/cmdscripts/EXAMPLE11.SHC

    r1107 r1108  
    22 
    33! Read everything within 7 degrees around epicenter from IRIS DC. 
    4 fdsnws 42.,130.6,7. * * * 6-jan-16_1:30 420 /addr=iris 
     4fdsnws 41.3,129.1,7. * * * 6-jan-16_1:30 420 /addr=iris 
    55 
    66! Check what we have (channels and location codes). 
     
    1515chanselect hh,bh ,00,10 
    1616 
     17! get missing metatdata 
     18meta all complete_location /addr=iris 
     19 
     20! sort by distance 
     21sort_by_distance 41.3 129.1 
     22 
     23! Make map plot with closest stations and epicenter 
     24plot_stations all /addloc=41.3,129.1 
     25 
    1726return 
  • SHX/trunk/SeismicHandler/cmdscripts/EXAMPLE12.SHC

    r1107 r1108  
    1818! subtract Time offset of 30s. 
    1919echo_ch/new evcatalog 
    20 events "lat "lon "radius "minmagn 30.0 
     20events "lat "lon "radius "minmagn ;; 30.0 
    2121echo_ch 
    2222 
     
    3838meta all complete /del_incomplete /addr="server 
    3939 
     40! Result of sorting depends on accuracy of epicenter (we may be very close). 
    4041sort_by_distance "lat "lon 
    4142 
     43! Show figure with station locations 
     44plot_stations all /addloc=|"lat|,|"lon| 
     45 
     46 
    4247return 
  • SHX/trunk/SeismicHandler/commands/help.py

    r1106 r1108  
    7373    # Read data from local Q file MYQFILE 
    7474read myqfile all 
     75 
     76    # Write to output file in various formats (see ObsPy write) 
     77writeo alltraces.mseeed all mseed 
     78writeo alltraces.gse all gse2 
     79writeo alltraces.sac all sac 
    7580     
    7681     ________________ 
  • SHX/trunk/SeismicHandler/commands/meta.py

    r1103 r1108  
    162162        else: 
    163163            level = 'response' 
     164        # request station data using bulk request 
     165        #print "dbg: bulk station request", incomplete 
    164166        client = Client(fdsnadr) 
    165167        inv = client.get_stations_bulk( incomplete, level=level ) 
     
    170172        for chan in inv.get_contents()['channels']: 
    171173            if not chan in startend.keys(): 
    172                 print "program bug, '%s' not in startend" % chan 
     174                print "got unrequested info for '%s', ignored" % chan 
    173175                continue 
    174176            stime, etime = startend[chan] 
Note: See TracChangeset for help on using the changeset viewer.