print 'Plotting event data'
print ' No filters'
print ' Average Slip vs Surface Rupture Length'
out_file = 'test_asrl.png'
start_time = time.time()
vcplots.average_slip_surface_rupture_length(sim_data_file, out_file)
print ' plotted to {}'.format(out_file)
print ' {} seconds'.format(time.time() - start_time)
print
print ' Magnitude vs Rupture Length'
out_file = 'test_mra.png'
start_time = time.time()
vcplots.magnitude_rupture_area(sim_data_file, out_file)
print ' plotted to {}'.format(out_file)
print ' {} seconds'.format(time.time() - start_time)
print
print ' Magnitude vs Average Slip'
out_file = 'test_mas.png'
start_time = time.time()
vcplots.magnitude_average_slip(sim_data_file, out_file)
print ' plotted to {}'.format(out_file)
print ' {} seconds'.format(time.time() - start_time)
print
print ' Frequency vs Magnitude'
out_file = 'test_fm.png'
start_time = time.time()
CUTOFF = 1000.0
PRE1 = False
PRE2 = True
BACK1 = False
BACK2 = False
EQ1 = False
EQ2 = False
SECS = (SOUTH_SAF_SECS,SOUTH_SAF_SECS)
field1dir = 'post_5_fields/'
field2dir = 'pre_5_fields/'
out_dir = 'local/'
TAGS = ['tejon_100mgal','northridge_100mgal','sf_100mgal']
PLOT_EVIDS = [60019,193054,231185]
"""
"""
center_evnum = 109382
with VCSimData() as sim_data:
sim_data.open_file(sim_file)
events = VCEvents(sim_data)
center_evyear = events.get_event_year(center_evnum)
start_year = round(center_evyear)-duration/2.0
end_year = round(center_evyear)+duration/2.0
event_range={'type':'year','filter':(start_year,end_year)}
event_range=None
pre = 'local/web_sim_full_'
output_files = [pre+'magn_rup_area.png',pre+'avg_slip_rup_len.png',pre+'freq_mag.png',pre+'recurr_intervals.png']
MAG_FILTER = None