def event_create(event_c): trace_start_lat = ones((event_c)) * -38.15 trace_start_lon = ones((event_c)) * 146.5 azimuth = ones((event_c)) * 217.2 dip = ones((event_c)) * 60.2 weight = ones((event_c)) * 18.1 event_activity = ones((event_c)) * 0.1 recurrence = ones((event_c)) * 5.3 Mw = ones((event_c)) * 6.9 lat0 = ones((event_c)) * -38.31 lon0 = ones((event_c)) * 146.3 depth = ones((event_c)) * 6.5 rx = ones((event_c)) * 25.3 ry = ones((event_c)) * 3.8 length = ones((event_c)) * 50.6 width = ones((event_c)) * 15. event_activity = ones((event_c)) * 9. event_set = Event_Set.create(depth=depth, rupture_centroid_lat=lat0, rupture_centroid_lon=lon0, azimuth=azimuth, dip=dip, ML=None, Mw=Mw, fault_width=15.0) event_set.event_activity = ones((event_c)) * 0.8 pseudo_event_set = Pseudo_Event_Set.split_logic_tree( event_set, ['Gaull_1990_WA'], [1]) pseudo_event_set.source_zone_id = ones((event_c)) * 9. return pseudo_event_set, event_activity
def event_create(event_c): trace_start_lat = ones((event_c))* -38.15 trace_start_lon = ones((event_c))* 146.5 azimuth = ones((event_c))* 217.2 dip = ones((event_c))* 60.2 weight = ones((event_c))* 18.1 event_activity = ones((event_c))* 0.1 recurrence = ones((event_c))* 5.3 Mw = ones((event_c))* 6.9 lat0 = ones((event_c))* -38.31 lon0 = ones((event_c))* 146.3 depth = ones((event_c))* 6.5 rx = ones((event_c))* 25.3 ry = ones((event_c))* 3.8 length = ones((event_c))* 50.6 width = ones((event_c))* 15. event_activity = ones((event_c))* 9. event_set = Event_Set.create(depth=depth,rupture_centroid_lat=lat0, rupture_centroid_lon=lon0,azimuth=azimuth, dip=dip,ML=None,Mw=Mw,fault_width=15.0) event_set.event_activity = ones((event_c))* 0.8 pseudo_event_set = Pseudo_Event_Set.split_logic_tree(event_set, ['Gaull_1990_WA'], [1]) pseudo_event_set.source_zone_id = ones((event_c))* 9. return pseudo_event_set, event_activity
def create_analysis_objects(self): # Parameters rupture_centroid_lat = asarray([-30]) rupture_centroid_lon = asarray([150]) length = asarray([1.0]) azimuth = asarray([2.0]) width = asarray([3.0]) dip = asarray([4.0]) depth = asarray([5.0]) Mw = asarray([6.0]) atten_models = asarray([ 'Allen', 'Toro_1997_midcontinent', 'Sadigh_97', 'Youngs_97_interface', 'Youngs_97_intraslab' ]) atten_model_weights = asarray([0.2, 0.2, 0.2, 0.2, 0.2]) atten_periods = asarray([0, 1.0, 2.0]) sites_lat = asarray([-31]) sites_lon = asarray([150]) # Event Set event_set = Event_Set.create(rupture_centroid_lat=rupture_centroid_lat, rupture_centroid_lon=rupture_centroid_lon, azimuth=azimuth, dip=dip, Mw=Mw, depth=depth, area=length * width, width=width, length=length) # Event Activity event_activity = Event_Activity(len(event_set)) event_activity.set_scenario_event_activity() # Source Model source_model = Source_Model.create_scenario_source_model( len(event_set)) source_model.set_attenuation(atten_models, atten_model_weights) source_model.set_ground_motion_calcs(atten_periods) event_set.scenario_setup() event_activity.ground_motion_model_logic_split(source_model, True) # Sites sites = Sites(sites_lat, sites_lon) # SA # Set up synthetic SA figures # Dimensions - spawn, gmm, rm, sites, events, period motion = zeros((1, len(atten_models), 1, len(sites), len(event_set), len(atten_periods)), dtype=float) # Allen motion[:, 0, :, :, :, 0] = 0 # period 0 motion[:, 0, :, :, :, 1] = 1 # period 1.0 motion[:, 0, :, :, :, 2] = 2 # period 2.0 # Toro_1997_midcontinent motion[:, 1, :, :, :, 0] = 3 # period 0 motion[:, 1, :, :, :, 1] = 4 # period 1.0 motion[:, 1, :, :, :, 2] = 5 # period 2.0 # Sadigh_97 motion[:, 2, :, :, :, 0] = 6 # period 0 motion[:, 2, :, :, :, 1] = 7 # period 1.0 motion[:, 2, :, :, :, 2] = 8 # period 2.0 # Youngs_97_interface motion[:, 3, :, :, :, 0] = 9 # period 0 motion[:, 3, :, :, :, 1] = 10 # period 1.0 motion[:, 3, :, :, :, 2] = 11 # period 2.0 # Young_97_intraslab motion[:, 4, :, :, :, 0] = 12 # period 0 motion[:, 4, :, :, :, 1] = 13 # period 1.0 motion[:, 4, :, :, :, 2] = 14 # period 2.0 # A minimal set of eqrm_flags so create_parameter_data passes # We only care about atten_models -> everything else are dummy values eqrm_flags = {} eqrm_flags['run_type'] = 'hazard' eqrm_flags['is_scenario'] = True eqrm_flags['output_dir'] = self.dir eqrm_flags['input_dir'] = self.dir eqrm_flags['site_tag'] = 'different_to_function' eqrm_flags['return_periods'] = [0.0] eqrm_flags['use_amplification'] = False eqrm_flags['zone_source_tag'] = 'not_used' eqrm_flags['atten_periods'] = atten_periods eqrm_flags['atten_models'] = atten_models return (event_set, event_activity, source_model, sites, motion, eqrm_flags)
def create_analysis_objects(self): # Parameters rupture_centroid_lat = asarray([-30]) rupture_centroid_lon = asarray([150]) length = asarray([1.0]) azimuth = asarray([2.0]) width = asarray([3.0]) dip = asarray([4.0]) depth = asarray([5.0]) Mw = asarray([6.0]) atten_models = asarray(['Allen', 'Toro_1997_midcontinent', 'Sadigh_97', 'Youngs_97_interface', 'Youngs_97_intraslab']) atten_model_weights = asarray([0.2, 0.2, 0.2, 0.2, 0.2]) atten_periods = asarray([0, 1.0, 2.0]) sites_lat = asarray([-31]) sites_lon = asarray([150]) # Event Set event_set = Event_Set.create(rupture_centroid_lat=rupture_centroid_lat, rupture_centroid_lon=rupture_centroid_lon, azimuth=azimuth, dip=dip, Mw=Mw, depth=depth, area=length*width, width=width, length=length) # Event Activity event_activity = Event_Activity(len(event_set)) event_activity.set_scenario_event_activity() # Source Model source_model = Source_Model.create_scenario_source_model(len(event_set)) source_model.set_attenuation(atten_models, atten_model_weights) source_model.set_ground_motion_calcs(atten_periods) event_set.scenario_setup() event_activity.ground_motion_model_logic_split(source_model, True) # Sites sites = Sites(sites_lat, sites_lon) # SA # Set up synthetic SA figures # Dimensions - spawn, gmm, rm, sites, events, period motion = zeros((1, len(atten_models), 1, len(sites), len(event_set), len(atten_periods)), dtype=float) # Allen motion[:,0,:,:,:,0] = 0 # period 0 motion[:,0,:,:,:,1] = 1 # period 1.0 motion[:,0,:,:,:,2] = 2 # period 2.0 # Toro_1997_midcontinent motion[:,1,:,:,:,0] = 3 # period 0 motion[:,1,:,:,:,1] = 4 # period 1.0 motion[:,1,:,:,:,2] = 5 # period 2.0 # Sadigh_97 motion[:,2,:,:,:,0] = 6 # period 0 motion[:,2,:,:,:,1] = 7 # period 1.0 motion[:,2,:,:,:,2] = 8 # period 2.0 # Youngs_97_interface motion[:,3,:,:,:,0] = 9 # period 0 motion[:,3,:,:,:,1] = 10 # period 1.0 motion[:,3,:,:,:,2] = 11 # period 2.0 # Young_97_intraslab motion[:,4,:,:,:,0] = 12 # period 0 motion[:,4,:,:,:,1] = 13 # period 1.0 motion[:,4,:,:,:,2] = 14 # period 2.0 # A minimal set of eqrm_flags so create_parameter_data passes # We only care about atten_models -> everything else are dummy values eqrm_flags = {} eqrm_flags['run_type'] = 'hazard' eqrm_flags['is_scenario'] = True eqrm_flags['output_dir'] = self.dir eqrm_flags['input_dir'] = self.dir eqrm_flags['site_tag'] = 'different_to_function' eqrm_flags['return_periods'] = [0.0] eqrm_flags['use_amplification'] = False eqrm_flags['zone_source_tag'] = 'not_used' eqrm_flags['atten_periods'] = atten_periods eqrm_flags['atten_models'] = atten_models return (event_set, event_activity, source_model, sites, motion, eqrm_flags)