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]