ztick=2
inverse_model=False
hypocenter=None
strike=155
xl=[0,30]
xorigin=325
yorigin=4710
yl=[0,50]
zl=[0,10]
#Proj stuff
p = Proj(proj='utm',zone=UTM_zone,ellps='WGS84')
fault_name=home+project_name+'/output/inverse_models/models/%s.%s.inv' % (run_name,run_number)
gmttools.make_total_model(fault_name,thresh=0)
fault=genfromtxt(home+project_name+'/output/inverse_models/models/%s.%s.inv.total' % (run_name,run_number))
reference=genfromtxt('/Users/dmelgar/Slip_inv/Amatrice_M6.6/output/inverse_models/models/gps_insar_vettore_main_synth_pian_anti_detach.0018.inv')
#Parse log file for hypocenter
log_file=home+project_name+'/output/inverse_models/models/%s.%s.log' % (run_name,run_number)
f=open(log_file,'r')
loop_go=True
while loop_go:
line=f.readline()
if 'Mw' in line:
Mw=float(line.split(':')[-1].split(' ')[-1])
loop_go=False
f.close()
def get_corners(meshfile):
from mudpy import gmttools
from string import rjust
path = '/Users/dmelgar/Slip_inv/Amatrice_3Dfitsgeol_final1/output/inverse_models/models/'
iter = 200
for k in range(iter):
mod = rjust(str(k), 4, '0')
gmttools.make_total_model(
u'/Users/dmelgar/Slip_inv/Amatrice_3Dfitsgeol_final1/output/inverse_models/models/jacknife.'
+ mod + '.inv', 0)
from numpy import zeros, genfromtxt, mean, std, savetxt, where, isnan, unique
from string import rjust
from mudpy import gmttools
path = '/Users/dmelgar/Slip_inv/puebla_herloc_NP2/output/inverse_models/models/'
iter = 86
nfaults = 27 * 17
mod = zeros((nfaults, iter))
for k in range(iter):
kread = rjust(str(k), 4, '0')
gmttools.make_total_model(
u'/Users/dmelgar/Slip_inv/puebla_herloc_NP2/output/inverse_models/models/jacknife.'
+ kread + '.inv', 0)
m = genfromtxt(
u'/Users/dmelgar/Slip_inv/puebla_herloc_NP2/output/inverse_models/models/jacknife.'
+ kread + '.inv.total')
slip = (m[:, 8]**2 + m[:, 9]**2)**0.5
mod[:, k] = slip
#Write models
mean_fault = m[0:nfaults, :].copy()
std_fault = m[0:nfaults, :].copy()
cv_fault = m[0:nfaults, :].copy()
mean_fault[:, 9] = mean(mod, 1)
mean_fault[:, 8] = 0
std_fault[:, 9] = std(mod, 1)
std_fault[:, 8] = 0
cv_fault[:, 9] = std(mod, 1) / mean(mod, 1)
cv_fault[:, 8] = 0
i = where(isnan(cv_fault[:, 9]) == True)[0]
from numpy import zeros, genfromtxt, mean, std, savetxt, where, isnan
from string import rjust
from mudpy import gmttools
path = '/Users/dmelgar/Slip_inv/Lefkada65/output/inverse_models/models/'
iter = 50
nfaults = 480
mod = zeros((nfaults, iter))
for k in range(iter):
kread = rjust(str(k), 4, '0')
gmttools.make_total_model(
u'/Users/dmelgar/Slip_inv/Lefkada65/output/inverse_models/models/jacknife.'
+ kread + '.inv', 0)
m = genfromtxt(
u'/Users/dmelgar/Slip_inv/Lefkada65/output/inverse_models/models/jacknife.'
+ kread + '.inv.total')
slip = (m[:, 8]**2 + m[:, 9]**2)**0.5
mod[:, k] = slip
#Write models
mean_fault = m.copy()
std_fault = m.copy()
cv_fault = m.copy()
mean_fault[:, 9] = mean(mod, 1)
mean_fault[:, 8] = 0
std_fault[:, 9] = std(mod, 1)
std_fault[:, 8] = 0
cv_fault[:, 9] = std(mod, 1) / mean(mod, 1)
cv_fault[:, 8] = 0
i = where(isnan(cv_fault[:, 9]) == True)[0]