import seismicbrowser
import obspy
dset = seismicbrowser.browseASDF('USCON_TA.h5')
# dset.read_TA_lst(infname='_US-TA-StationList.txt', channel='LHE,LHN,LHZ', maxlatitude=50.)
# dset.write_inv('ALASKA_TA_AK.xml')
dset.read_inv('ALASKA_TA_AK.xml')
#
#
# dset = seismicbrowser.browseASDF('ALASKA_BH.h5')
# # dset.get_stations(minlatitude=52., maxlatitude=72.5, minlongitude=-172., maxlongitude=-122., channel='BH*', startafter='19910101')
# # dset.write_inv('ALASKA_BH.xml')
# dset.read_inv('ALASKA_BH.xml')
#
#
# dset.plot_inv()
def to_percent(y, position):
# Ignore the passed in position. This has the effect of scaling the default
# tick locations.
s = '%.0f' % (100. * y)
# s = str(y)
# The percent symbol needs escaping in latex
if matplotlib.rcParams['text.usetex'] is True:
return s + r'$\%$'
else:
return s + '%'
import seismicbrowser
import obspy
import GeoPolygon
# dset = seismicbrowser.browseASDF('glims_BH.h5')
dset = seismicbrowser.browseASDF('ALASKA_BH.h5')
# dset.read_inv('glims_5000.xml')
# dset.read_inv('ALASKA_BH.xml')
dset.write_anxcorr_txt('sta_glims.lst')
import seismicbrowser
import obspy
dset = seismicbrowser.browseASDF('EARS_LH.h5')
# dset.get_stations(minlatitude=-17.78, maxlatitude=4.15, minlongitude=23.00, maxlongitude=40.2, channel='LH?', startafter='19940101')
# dset.write_inv('EARS_LH.xml')
dset.read_inv('EARS_LH.xml')
dset.write_txt(outfname='ASN/EARS_LH_breqfast.txt')
#
# dset = seismicbrowser.browseASDF('ALASKA_LH.h5')
# # dset.get_stations(minlatitude=52., maxlatitude=72.5, minlongitude=-172., maxlongitude=-122., channel='BH*', startafter='19910101')
# # dset.write_inv('ALASKA_BH.xml')
# dset.read_inv('ALASKA.xml')
# # # #
# # # #
# # # # dset.plot_inv(netcodelist=['AK', 'TA', 'XN', 'XY', 'AV', 'PO', 'US', 'CN'])
# # dset.plot_topo()
# # dset.plot_inv(projection='merc')
# dset.plot_inv2(plotetopo=True)