#!/usr/bin/env python # -*- coding: utf-8 -*- """Process misfit """ import sys from misfit import Misfit # read command line args misfit_file = str(sys.argv[1]) out_file = str(sys.argv[2]) print("\n====== initialize\n") misfit = Misfit() print("\n====== load data\n") misfit.load(misfit_file) print("\n====== output misfit measurments\n") misfit.output_misfit(out_file)
#------ load parameter file if sys.version_info < (3, ): raise Exception("need python3") elif sys.version_info < (3, 5): spec =importlib.machinery.SourceFileLoader("misfit_par", par_file) par = spec.load_module() else: spec = importlib.util.spec_from_file_location("misfit_par", par_file) par = importlib.util.module_from_spec(spec) spec.loader.exec_module(par) #------ print("\n====== initialize\n") misfit = Misfit() print("\n====== load data\n") misfit.load(misfit_file) print("\n====== grid_cc \n") wcc_sum, weight_sum = \ misfit.cc_linearized_seismogram_for_dmodel( dm=par.dm_model, plot=False) with open(out_file, 'w') as f: f.write("#weight_sum = {:12.5e}\n".format(weight_sum)) f.write("#step_size wcc_sum/weight_sum\n") for idx in range(len(wcc_sum)): f.write("{:12.5e} {:15.8e}\n".format( par.dm_model['model'][idx], wcc_sum[idx]/weight_sum))
#!/usr/bin/env python # -*- coding: utf-8 -*- """Process misfit """ import sys from misfit import Misfit import numpy as np # read command line args #data_dir = "DATA" #cmt_file = "DATA/CMTSOLUTION.init" #channel_file = "DATA/channel.txt" misfit_file = "misfit/misfit.pkl" # print("\n====== initialize\n") misfit = Misfit() #print("\n====== setup event\n") #misfit.setup_event(cmt_file, is_ECEF=True) # #print("\n====== setup station\n") #misfit.setup_station(channel_file) print("\n====== load data\n") misfit.load(filename=misfit_file) print("\n====== plot radiation pattern\n") misfit.radiation_pattern()
#!/usr/bin/env python # -*- coding: utf-8 -*- """Process misfit """ import sys from misfit import Misfit import numpy as np # read command line args #data_dir = "DATA" #cmt_file = "DATA/CMTSOLUTION.init" #channel_file = "DATA/channel.txt" misfit_file = "misfit/misfit.pkl" # print("\n====== initialize\n") misfit = Misfit() #print("\n====== setup event\n") #misfit.setup_event(cmt_file, is_ECEF=True) # #print("\n====== setup station\n") #misfit.setup_station(channel_file) print("\n====== load data\n") misfit.load(filename=misfit_file) print("\n====== plot radiation pattern\n") misfit.radiation_pattern()