#!/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()
