plotnum = [False for i in range(15)]
if len(sys.argv) < 2:
print help
exit()
else:
plotnum[int(sys.argv[1])] = True
####################################################
## Figure Properties #######
width = 12
height = 9
legsize = width + 3
###########################
## Prep Data
arg = Args()
arg.addQuery("status", "eq", "processed-ok")
# Load station params
d = Params(stnfile, ["hk::H","hk::R"], arg)
arg = Args()
m = Params(moonfile, ["H","Vp"])
m.filter(arg.addQuery("geoprov", "not in", "oceanic"))
m.filter(arg.addQuery("geoprov", "not in", "Shelf"))
#######################################################################
# F0
# CANADA
if plotnum[0]:
## Kanamori Poisson Histogram
# Set figure
fig = plt.figure( figsize = (width, height) )
plotnum = [False for i in range(15)]
if len(sys.argv) < 2:
print help
exit()
else:
plotnum[int(sys.argv[1])] = True
####################################################
## Figure Properties #######
width = 12
height = 5
legsize = height + 3
###########################
## Prep Data
arg = Args()
arg.addQuery("status", "eq", "processed-ok")
# Load station params
d = Params(stnfile, ["hk::H", "hk::R"], arg)
arg = Args()
m = Params(moonfile, ["H", "Vp"])
m.filter(arg.addQuery("geoprov", "not in", "oceanic"))
m.filter(arg.addQuery("geoprov", "not in", "Shelf"))
#######################################################################
# F0
# CANADA
if plotnum[0]:
## Kanamori Poisson Histogram
# Set figure
fig = plt.figure(figsize=(width, height))
'''
if len(sys.argv) < 2:
print help
else:
plotnum[int(sys.argv[1])] = True
def poisson(R):
''' Function to go from Vp/Vs -> Poisson's ratio '''
return ( (R**2 - 2) / (2*(R**2 - 1)))
#######################################################################
# F0
if plotnum[0]:
arg1 = Args()
arg1.addQuery("hk::H", "lt", "45")
arg2 = Args()
arg2.addQuery("hk::H", "lt", "40")
arg3 = Args()
arg3.addQuery("hk::H", "lt", "35")
arg4 = Args()
fstns = os.environ['HOME'] + '/thesis/stations.json'
a = Params(fstns, arg1, ["mb::H"])
b = Params(fstns, arg2, ["mb::H"])
c = Params(fstns, arg2, ["hk::H"])
d = Params(fstns, arg3, ["hk::H"])
e = Params(fstns, arg3, ["wm::H"])
f = Params(fstns, arg4, ["wm::H"])
if len(sys.argv) < 2:
print help
else:
plotnum[int(sys.argv[1])] = True
def poisson(R):
''' Function to go from Vp/Vs -> Poisson's ratio '''
return ((R**2 - 2) / (2 * (R**2 - 1)))
#######################################################################
# F0
if plotnum[0]:
arg1 = Args()
arg1.addQuery("hk::H", "lt", "45")
arg2 = Args()
arg2.addQuery("hk::H", "lt", "40")
arg3 = Args()
arg3.addQuery("hk::H", "lt", "35")
arg4 = Args()
fstns = os.environ['HOME'] + '/thesis/stations.json'
a = Params(fstns, arg1, ["mb::H"])
b = Params(fstns, arg2, ["mb::H"])
c = Params(fstns, arg2, ["hk::H"])
d = Params(fstns, arg3, ["hk::H"])
e = Params(fstns, arg3, ["wm::H"])
f = Params(fstns, arg4, ["wm::H"])
print "total stations examined:", len(d.stns)
for jj in range(0,4):
print msg[jj]
ix = ixs[jj]
for ii, stn in enumerate(d.stns[ix]):
if details:
print stn, d.hk_c0R[ix][ii], d.hk_c1R[ix][ii]
else:
pass
print "number of stations:", len(d.stns[ix])
#details = False
#printinfo(msg, ixs, d, details)
arg = Args().stations( list(d.stns[ixs[2]]) )
arg.addQuery("status", "eq", "processed-notok")
d.filter(arg)
for stn in d.stns:
print stn
d = Params(stnfile, ["hk::H","hk::R", "hk::stdR", "hk::c0R", "hk::c1R"])
d.filter(Args().addQuery("status", "in", "processed"))
arg = Args()
arg.stations(list(d.stns[ixs[0]]))
arg.addQuery("status", "eq", "processed-ok")
d.filter(arg)
#for station in d.stns:
# print station
#print np.fabs(d.hk_c1R - d.hk_c0R)
print "total stations examined:", len(d.stns)
for jj in range(0, 4):
print msg[jj]
ix = ixs[jj]
for ii, stn in enumerate(d.stns[ix]):
if details:
print stn, d.hk_c0R[ix][ii], d.hk_c1R[ix][ii]
else:
pass
print "number of stations:", len(d.stns[ix])
#details = False
#printinfo(msg, ixs, d, details)
arg = Args().stations(list(d.stns[ixs[2]]))
arg.addQuery("status", "eq", "processed-notok")
d.filter(arg)
for stn in d.stns:
print stn
d = Params(stnfile, ["hk::H", "hk::R", "hk::stdR", "hk::c0R", "hk::c1R"])
d.filter(Args().addQuery("status", "in", "processed"))
arg = Args()
arg.stations(list(d.stns[ixs[0]]))
arg.addQuery("status", "eq", "processed-ok")
d.filter(arg)
#for station in d.stns:
# print station
#print np.fabs(d.hk_c1R - d.hk_c0R)
plotnum = [False for i in range(15)]
if len(sys.argv) < 2:
print help
exit()
else:
plotnum[int(sys.argv[1])] = True
####################################################
## Figure Properties #######
width = 12
height = 5
legsize = height + 3
###########################
## Prep Data
arg = Args()
arg.addQuery("status", "eq", "processed-ok")
# Load station params
d = Params(stnfile, ["hk::H","hk::R"], arg)
arg = Args()
m = Params(moonfile, ["H","Vp"])
m.filter(arg.addQuery("geoprov", "not in", "oceanic"))
m.filter(arg.addQuery("geoprov", "not in", "Shelf"))
#######################################################################
# F0
# CANADA
if plotnum[0]:
## Kanamori Poisson Histogram
# Set figure
fig = plt.figure( figsize = (width, height) )
