コード例 #1
0
def FDSN_network(input_dics, events):
    """
    Availability and retrieivng functions are called with FDSN_network
    :param input_dics:
    :param events:
    :return:
    """
    print "\n********************************************************"
    print "FDSN -- Download waveforms, StationXML files and meta-data"
    print "********************************************************"
    period = "{0:s}_{1:s}".format(input_dics["min_date"].split("T")[0], input_dics["max_date"].split("T")[0])
    eventpath = os.path.join(input_dics["datapath"], period)

    print "Initializing folders and files...",
    create_folders_files(events, eventpath, input_dics)
    print "DONE"

    for i in range(len(events)):
        t_fdsn_1 = datetime.now()
        target_path = os.path.join(eventpath, events[i]["event_id"])

        if not input_dics["list_stas"]:
            Stas_fdsn = FDSN_available(input_dics, events[i], target_path, event_number=i)
        else:
            Stas_fdsn = read_list_stas(input_dics["list_stas"], input_dics["normal_mode_syn"], input_dics["specfem3D"])
            if input_dics["fdsn_bulk"] == "Y":
                FDSN_create_bulk_list(target_path, input_dics, Stas_fdsn, events[i])

        if input_dics["fdsn_bulk"] != "Y":
            print "\n%s-Availability for event: %s/%s ---> DONE" % (input_dics["fdsn_base_url"], i + 1, len(events))
        else:
            print "\nFDSN-bulkfile for event: %s/%s ---> DONE" % (i + 1, len(events))

        print "Time for checking the availability: %s" % (datetime.now() - t_fdsn_1)

        if Stas_fdsn != [[]]:
            FDSN_waveform(input_dics, events, Stas_fdsn, i, req_type="save")
            fdsn_url_orig = input_dics["fdsn_base_url"]
            for cli_rest in input_dics["fdsn_base_url_rest"]:
                try:
                    from update_handler import FDSN_update

                    input_dics["fdsn_base_url"] = cli_rest
                    input_dics["fdsn_update"] = target_path
                    FDSN_update(input_dics, address=target_path)
                except Exception, e:
                    print e
                    continue
            if len(input_dics["fdsn_base_url_rest"]) > 0:
                input_dics["fdsn_base_url"] = fdsn_url_orig
                input_dics["fdsn_update"] = "N"
        else:
            print "No available station in %s for your request and " "for event %s!" % (
                input_dics["fdsn_base_url"],
                str(i + 1),
            )
            continue
コード例 #2
0
ファイル: update_handler.py プロジェクト: Nanyaohui/obspyDMT
def FDSN_update(input_dics, address):
    """
    Initialize directories and pass the required stations
    (removing the duplications) for FDSN update requests
    :param input_dics:
    :param address:
    :return:
    """
    print '\n*********************'
    print 'FDSN -- Updating Mode'
    print '*********************'
    t_update_1 = datetime.now()

    events, address_events = quake_info(address, 'info')

    for i in range(len(events)):
        if not input_dics['list_stas']:
            Stas_fdsn = FDSN_available(input_dics,
                                       events[i],
                                       address_events[i],
                                       event_number=i)
        else:
            Stas_fdsn = read_list_stas(input_dics['list_stas'],
                                       input_dics['normal_mode_syn'],
                                       input_dics['specfem3D'])

        if input_dics['fdsn_bulk'] != 'Y':
            print '\n%s-Availability for event: %s/%s ---> DONE' \
                  % (input_dics['fdsn_base_url'], i+1, len(events))
        else:
            print '\nFDSN-bulkfile for event: %s/%s ---> DONE' \
                  % (i+1, len(events))

        if Stas_fdsn != [[]]:
            Stas_req = rm_duplicate(Stas_fdsn,
                                    address=os.path.join(address_events[i]))
        else:
            Stas_req = None
            print '------------------------------'
            print 'There is no available station!'
            print '------------------------------'

        if not os.path.isdir(os.path.join(address_events[i], 'BH_RAW')):
            os.makedirs(os.path.join(address_events[i], 'BH_RAW'))

        if Stas_req:
            FDSN_waveform(input_dics, events, Stas_req, i, req_type='update')
        else:
            print '\nNo available station in %s for your request!' \
                  % input_dics['fdsn_base_url']
            print 'Check the next event...'
            continue

    print '\nTotal time for updating FDSN: %s' % (datetime.now() - t_update_1)
コード例 #3
0
def FDSN_network(input_dics, events):
    """
    Availability and retrieivng functions are called with FDSN_network
    :param input_dics:
    :param events:
    :return:
    """
    print '\n********************************************************'
    print 'FDSN -- Download waveforms, StationXML files and meta-data'
    print '********************************************************'
    period = '{0:s}_{1:s}_{2:s}_{3:s}'.format(
        input_dics['min_date'].split('T')[0],
        input_dics['max_date'].split('T')[0],
        str(input_dics['min_mag']),
        str(input_dics['max_mag']))
    eventpath = os.path.join(input_dics['datapath'], period)

    print 'Initializing folders and files...',
    create_folders_files(events, eventpath, input_dics)
    print 'DONE'

    for i in range(len(events)):
        t_fdsn_1 = datetime.now()
        target_path = os.path.join(eventpath, events[i]['event_id'])

        if not input_dics['list_stas']:
            Stas_fdsn = FDSN_available(input_dics, events[i],
                                       target_path, event_number=i)
        else:
            Stas_fdsn = read_list_stas(input_dics['list_stas'],
                                       input_dics['normal_mode_syn'],
                                       input_dics['specfem3D'])
            if input_dics['fdsn_bulk'] == 'Y':
                FDSN_create_bulk_list(target_path, input_dics, Stas_fdsn,
                                      events[i])

        if input_dics['fdsn_bulk'] != 'Y':
            print '\n%s-Availability for event: %s/%s ---> DONE' \
                  % (input_dics['fdsn_base_url'], i+1, len(events))
        else:
            print '\nFDSN-bulkfile for event: %s/%s ---> DONE' \
                  % (i+1, len(events))

        print 'Time for checking the availability: %s' \
              % (datetime.now() - t_fdsn_1)

        if Stas_fdsn != [[]]:
            FDSN_waveform(input_dics, events, Stas_fdsn, i, req_type='save')
        else:
            print 'No available station in %s for your request and ' \
                  'for event %s!' % (input_dics['fdsn_base_url'], str(i+1))
            continue
コード例 #4
0
ファイル: fdsn_handler.py プロジェクト: qingkaikong/obspyDMT
def FDSN_network(input_dics, events):
    """
    Availability and retrieivng functions are called with FDSN_network
    :param input_dics:
    :param events:
    :return:
    """
    print '\n********************************************************'
    print 'FDSN -- Download waveforms, StationXML files and meta-data'
    print '********************************************************'
    period = '{0:s}_{1:s}'.format(
        input_dics['min_date'].split('T')[0],
        input_dics['max_date'].split('T')[0])
    eventpath = os.path.join(input_dics['datapath'], period)

    print 'Initializing folders and files...',
    create_folders_files(events, eventpath, input_dics)
    print 'DONE'

    for i in range(len(events)):
        t_fdsn_1 = datetime.now()
        target_path = os.path.join(eventpath, events[i]['event_id'])

        if not input_dics['list_stas']:
            Stas_fdsn = FDSN_available(input_dics, events[i],
                                       target_path, event_number=i)
        else:
            Stas_fdsn = read_list_stas(input_dics['list_stas'],
                                       input_dics['normal_mode_syn'],
                                       input_dics['specfem3D'])
            if input_dics['fdsn_bulk'] == 'Y':
                FDSN_create_bulk_list(target_path, input_dics, Stas_fdsn,
                                      events[i])

        if input_dics['fdsn_bulk'] != 'Y':
            print '\n%s-Availability for event: %s/%s ---> DONE' \
                  % (input_dics['fdsn_base_url'], i+1, len(events))
        else:
            print '\nFDSN-bulkfile for event: %s/%s ---> DONE' \
                  % (i+1, len(events))

        print 'Time for checking the availability: %s' \
              % (datetime.now() - t_fdsn_1)

        if Stas_fdsn != [[]]:
            FDSN_waveform(input_dics, events, Stas_fdsn, i, req_type='save')
        else:
            print 'No available station in %s for your request and ' \
                  'for event %s!' % (input_dics['fdsn_base_url'], str(i+1))
            continue
コード例 #5
0
ファイル: update_handler.py プロジェクト: Nanyaohui/obspyDMT
def ARC_update(input_dics, address):
    """
    Initialize directories and pass the required stations
    (removing the duplications) for ArcLink update requests
    :param input_dics:
    :param address:
    :return:
    """
    print '\n************************'
    print 'ArcLink -- Updating Mode'
    print '************************'
    t_update_1 = datetime.now()

    events, address_events = quake_info(address, 'info')

    for i in range(len(events)):
        if not input_dics['list_stas']:
            Stas_arc = ARC_available(input_dics,
                                     events[i],
                                     address_events[i],
                                     event_number=i)
        else:
            Stas_arc = read_list_stas(input_dics['list_stas'],
                                      input_dics['normal_mode_syn'],
                                      input_dics['specfem3D'])

        print '\nArcLink-Availability for event: %s/%s ---> DONE' \
              % (i+1, len(events))

        if Stas_arc != [[]]:
            Stas_req = rm_duplicate(Stas_arc,
                                    address=os.path.join(address_events[i]))
        else:
            Stas_req = None
            print '------------------------------'
            print 'There is no available station!'
            print '------------------------------'

        if not os.path.isdir(os.path.join(address_events[i], 'BH_RAW')):
            os.makedirs(os.path.join(address_events[i], 'BH_RAW'))

        if Stas_req:
            ARC_waveform(input_dics, events, Stas_req, i, req_type='update')
        else:
            print '\nNo available station in ArcLink for your request!'
            print 'Check the next event...'
            continue

    print '\nTotal time for updating ArcLink: %s' \
          % (datetime.now() - t_update_1)
コード例 #6
0
def ARC_network(input_dics, events):
    """
    Returns information about what time series data is available
    at the ArcLink nodes for all requested events
    """
    print '\n***********************************************************'
    print 'ArcLink -- Download waveforms, response files and meta-data'
    print '***********************************************************'
    period = '{0:s}_{1:s}_{2:s}_{3:s}'.format(
        input_dics['min_date'].split('T')[0],
        input_dics['max_date'].split('T')[0],
        str(input_dics['min_mag']),
        str(input_dics['max_mag']))
    eventpath = os.path.join(input_dics['datapath'], period)

    if input_dics['FDSN'] != 'Y':
        print 'Initializing folders and files...',
        create_folders_files(events, eventpath, input_dics)
        print 'DONE'

    for i in range(len(events)):
        t_arc_1 = datetime.now()
        target_path = os.path.join(eventpath, events[i]['event_id'])

        if not input_dics['list_stas']:
            Stas_arc = ARC_available(input_dics, events[i],
                                     target_path, event_number=i)
        else:
            Stas_arc = read_list_stas(input_dics['list_stas'],
                                      normal_mode_syn='N',
                                      specfem3D='N')

        print '\nArcLink-Availability for event: %s/%s  ---> DONE' \
              % (i+1, len(events))
        print 'Time for checking the availability: %s' \
              % (datetime.now() - t_arc_1)

        if Stas_arc != [[]]:
            ARC_waveform(input_dics, events, Stas_arc, i, req_type='save')
        else:
            print 'No available station in ArcLink for your request and ' \
                  'for event %s!' % str(i+1)
            continue
コード例 #7
0
def ARC_network(input_dics, events):
    """
    Returns information about what time series data is available
    at the ArcLink nodes for all requested events
    """
    print '\n***********************************************************'
    print 'ArcLink -- Download waveforms, response files and meta-data'
    print '***********************************************************'
    period = '{0:s}_{1:s}_{2:s}_{3:s}'.format(
        input_dics['min_date'].split('T')[0],
        input_dics['max_date'].split('T')[0],
        str(input_dics['min_mag']),
        str(input_dics['max_mag']))
    eventpath = os.path.join(input_dics['datapath'], period)

    if input_dics['FDSN'] != 'Y':
        print 'Initializing folders and files...',
        create_folders_files(events, eventpath, input_dics)
        print 'DONE'

    for i in range(len(events)):
        t_arc_1 = datetime.now()
        target_path = os.path.join(eventpath, events[i]['event_id'])

        if not input_dics['list_stas']:
            Stas_arc = ARC_available(input_dics, events[i],
                                     target_path, event_number=i)
        else:
            Stas_arc = read_list_stas(input_dics['list_stas'],
                                      normal_mode_syn='N',
                                      specfem3D='N')

        print '\nArcLink-Availability for event: %s/%s  ---> DONE' \
              % (i+1, len(events))
        print 'Time for checking the availability: %s' \
              % (datetime.now() - t_arc_1)

        if Stas_arc != [[]]:
            ARC_waveform(input_dics, events, Stas_arc, i, req_type='save')
        else:
            print 'No available station in ArcLink for your request and ' \
                  'for event %s!' % str(i+1)
            continue