def post_process(self):
"""
Run the standard BBP post-processing tasks
"""
# Plot seismograms
plotter = PlotSeis(os.path.basename(self.station_list),
os.path.basename(self.src_files[0]),
True, True, self.output_sim_id)
plotter.run()
# RotD50
process = RotD50(os.path.basename(self.station_list),
self.output_sim_id)
process.run()
process = ObsSeismograms(os.path.basename(self.station_list),
self.val_obj.get_obs_path(),
self.val_obj.get_obs_format(),
os.path.basename(self.val_obj.get_obs_corrections()),
self.output_sim_id)
process.run()
process = GenPlots(os.path.basename(self.station_list),
self.val_obj.get_obs_path(),
'acc',
self.val_obj.get_validation_name(),
self.output_sim_id)
process.run()
process = GPGof(os.path.basename(self.src_files[0]),
os.path.basename(self.station_list),
self.val_obj.get_magnitude(),
self.val_obj.get_validation_name(),
self.val_obj.get_cutoff(),
False,
self.output_sim_id)
process.run()
def post_process(station_list, src_files,
merged_outdir,
realization, scenario):
"""
Run the standard BBP post-processing tasks
"""
# Plot seismograms
plotter = PlotSeis(os.path.basename(station_list),
os.path.basename(src_files[0]),
True, True, int(realization))
plotter.run()
# RotD50
process = RotD50(os.path.basename(station_list),
int(realization))
process.run()
# Plot RotD50
# Load station list
slo = StationList(station_list)
site_list = slo.getStationList()
for site in site_list:
stat = site.scode
rd50_file = "%d.%s.rd50" % (int(realization), stat)
rd50_filename1 = os.path.join(merged_outdir, rd50_file)
outfile = os.path.join(merged_outdir, "%s_%d_%s_rotd50.png" %
(scenario, int(realization), stat))
plot_rotd50.plot_rd50(stat, rd50_filename1, "-",
scenario, "-", outfile,
site.low_freq_corner,
site.high_freq_corner,
quiet=True)
def post_process(station_list, src_files,
merged_outdir,
realization, val_obj):
"""
Run the standard BBP post-processing tasks
"""
# Plot seismograms
plotter = PlotSeis(os.path.basename(station_list),
os.path.basename(src_files[0]),
True, True, int(realization))
plotter.run()
# RotD50
process = RotD50(os.path.basename(station_list),
int(realization))
process.run()
process = ObsSeismograms(os.path.basename(station_list),
val_obj.get_obs_path(),
val_obj.get_obs_format(),
val_obj.get_obs_corrections(),
int(realization))
process.run()
process = GenPlots(os.path.basename(station_list),
val_obj.get_obs_path(),
'acc',
val_obj.get_validation_name(),
int(realization))
process.run()
process = GPGof(os.path.basename(src_files[0]),
os.path.basename(station_list),
val_obj.get_magnitude(),
val_obj.get_validation_name(),
val_obj.get_cutoff(),
False,
int(realization))
process.run()
def post_process(station_list, src_files,
merged_outdir,
realization, val_obj):
"""
Run the standard BBP post-processing tasks
"""
# Plot seismograms
plotter = PlotSeis(os.path.basename(station_list),
os.path.basename(src_files[0]),
True, True, int(realization))
plotter.run()
# RotD50
process = RotD50(os.path.basename(station_list),
int(realization))
process.run()
process = ObsSeismograms(os.path.basename(station_list),
val_obj.get_obs_path(),
val_obj.get_obs_format(),
val_obj.get_obs_corrections(),
int(realization))
process.run()
process = GenPlots(os.path.basename(station_list),
'acc',
val_obj.get_validation_name(),
int(realization))
process.run()
process = GPGof(os.path.basename(src_files[0]),
os.path.basename(station_list),
val_obj.get_validation_name(),
val_obj.get_cutoff(),
False,
int(realization))
process.run()
def post_process(self):
"""
Run the standard BBP post-processing tasks
"""
# Plot seismograms
plotter = PlotSeis(os.path.basename(self.station_list),
os.path.basename(self.src_files[0]),
True, True, self.output_sim_id)
plotter.run()
# RotD50
process = RotD50(os.path.basename(self.station_list),
self.output_sim_id)
process.run()
# Plot RotD50
# Load station list
slo = StationList(self.station_list)
site_list = slo.getStationList()
for site in site_list:
stat = site.scode
rd50_file = "%d.%s.rd50" % (self.output_sim_id, stat)
rd50_filename1 = os.path.join(self.a_outdir, rd50_file)
outfile = os.path.join(self.a_outdir, "%s_%d_%s_rotd50.png" %
(self.scenario, self.output_sim_id, stat))
plot_rotd50.plot_rd50(stat, rd50_filename1, "-",
self.scenario, "-", outfile,
site.low_freq_corner,
site.high_freq_corner,
quiet=True)