Beispiel #1
0
def old_max_simulation():
    """
    attribute	#	JSON attribute
    events	8,000,000*	len_events
    Sites 	3,000,000	 len_block_sites
    periods	10	len_atten_periods
    return periods	10	len_return_periods
    source zones	150	len_source_zones
    SA surfaces (Sa soil and Sa bed rock)	2	len_SA_surfaces
# of GM models	5	len_max_GMPEs

   Note, i'm not doing 150 source zones.  The mem estimate's have been
   good, without taking it into account.
   """
    runs = []
    num_sources = 4
    sdp = create_base()
    
    ###
    sdp.site_indexes = range(1, 400)
    events = 80000  
    events = 800 
    events = 8000000 
    events = 80    
    sdp.return_periods = [100.75, 200.0, 300.0, 400.0, 500.0, 
    600.0, 700.0, 800.0, 900, 1000.]
    sdp.atten_periods =  [0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0]
    sdp.use_amplification = True
    sdp.atten_spawn_bins = 1
    sites_needed = 17278
    sdp.file_parallel_log_level = 'debug'
    sdp.use_site_indexes = False
    ###
    num_sites = 17277
    events_source = events/num_sources
    sdp.prob_number_of_events_in_zones = [events_source]*num_sources
    sdp.prob_number_of_events_in_zones[0] += events%num_sources
    sdp.use_site_indexes = True
    sdp.site_indexes = create_sites_list(num_sites, sites_needed)
    assert sum(sdp.prob_number_of_events_in_zones) == events
    assert len(sdp.site_indexes) == sites_needed
    sdp.event_control_tag = "4GMPE" 
    dir_last  = 'old_max_sim_again_' + str(sum(
            sdp.prob_number_of_events_in_zones)) + \
            "_sites" + str(len(sdp.site_indexes))
    sdp.output_dir = os.path.join(eqrm_data_home(), 'test_national', 
                                  'memory_output', dir_last)
    runs.append({"processes":1, "sdp":sdp})

    return runs
Beispiel #2
0
def risk_simulation():
    """
    
   
    """
    runs = []
    num_sources = 6  
    
    ### 
    
    
    sdp = create_base_risk()
    sdp.atten_collapse_Sa_of_atten_models = True 
    sdp.atten_variability_method = 2
    ###
    sites_needed = 1
    events = 80
    ###
    events_source = events/num_sources
    num_sites = 17278
    sdp.prob_number_of_events_in_zones = [events_source]*num_sources
    sdp.prob_number_of_events_in_zones[0] += events%num_sources
    sdp.site_indexes = create_sites_list(num_sites, sites_needed)
    dir_last  = 'initial_risk_' + str(events) + \
            "_sites" + str(sites_needed)
    sdp.output_dir = os.path.join(eqrm_data_home(), 'test_national', 
                                  'memory_output_riskB', dir_last)
    runs.append({"processes":1, "sdp":sdp})

 
    sdp = create_base_risk()   
    sdp.atten_collapse_Sa_of_atten_models = True 
    sdp.atten_variability_method = 2
    ###
    sites_needed = 2
    events = 80
    ###
    events_source = events/num_sources
    num_sites = 17278
    sdp.prob_number_of_events_in_zones = [events_source]*num_sources
    sdp.prob_number_of_events_in_zones[0] += events%num_sources
    sdp.site_indexes = create_sites_list(num_sites, sites_needed)
    dir_last  = 'initial_risk_' + str(events) + \
            "_sites" + str(sites_needed)
    sdp.output_dir = os.path.join(eqrm_data_home(), 'test_national', 
                                  'memory_output_riskB', dir_last)
    runs.append({"processes":1, "sdp":sdp})

    sdp = create_base_risk()  
    sdp.atten_collapse_Sa_of_atten_models = True 
    sdp.atten_variability_method = 2  
    ###
    sites_needed = 1
    events = 800
    ###
    events_source = events/num_sources
    num_sites = 17278
    sdp.prob_number_of_events_in_zones = [events_source]*num_sources
    sdp.prob_number_of_events_in_zones[0] += events%num_sources
    sdp.site_indexes = create_sites_list(num_sites, sites_needed)
    dir_last  = 'initial_risk_' + str(events) + \
            "_sites" + str(sites_needed)
    sdp.output_dir = os.path.join(eqrm_data_home(), 'test_national', 
                                  'memory_output_risk', dir_last)
    #runs.append({"processes":1, "sdp":sdp})
    
 
    sdp = create_base_risk()   
    sdp.atten_collapse_Sa_of_atten_models = True 
    sdp.atten_variability_method = 2
    ###
    sites_needed = 10
    events = 800
    ###
    events_source = events/num_sources
    num_sites = 17278
    sdp.prob_number_of_events_in_zones = [events_source]*num_sources
    sdp.prob_number_of_events_in_zones[0] += events%num_sources
    sdp.site_indexes = create_sites_list(num_sites, sites_needed)
    dir_last  = 'initial_risk_' + str(events) + \
            "_sites" + str(sites_needed)
    sdp.output_dir = os.path.join(eqrm_data_home(), 'test_national', 
                                  'memory_output_risk', dir_last)
    #runs.append({"processes":1, "sdp":sdp})
    

    sdp = create_base_risk()  
    sdp.atten_collapse_Sa_of_atten_models = True 
    sdp.atten_variability_method = 2  
    ###
    sites_needed = 1
    events = 80000
    ###
    events_source = events/num_sources
    num_sites = 17278
    sdp.prob_number_of_events_in_zones = [events_source]*num_sources
    sdp.prob_number_of_events_in_zones[0] += events%num_sources
    sdp.site_indexes = create_sites_list(num_sites, sites_needed)
    dir_last  = 'initial_risk_' + str(events) + \
            "_sites" + str(sites_needed)
    sdp.output_dir = os.path.join(eqrm_data_home(), 'test_national', 
                                  'memory_output_risk', dir_last)
    #runs.append({"processes":1, "sdp":sdp})

    
    return runs