mmax.append(neo_mmax[matchidx]) usd.append(neo_usd[matchidx]) lsd.append(neo_lsd[matchidx]) ow_lsd.append(neo_ow_lsd[matchidx]) stk.append(neo_stk[matchidx]) dip.append(neo_dip[matchidx]) rke.append(neo_rke[matchidx]) lsd = array(lsd) dom = array(dom) dom[1] = 2. ############################################################################### # load 2018 completeness models ############################################################################### ycomp, mcomp, min_rmag = get_completeness_model(src_codes, shapes, dom, 1) min_rmag[0] = 3.2 ############################################################################### # load Rajabi SHMax vectors ############################################################################### shmax_pref, shmax_sig = get_aus_shmax_vectors(src_codes, shapes) # fix preferred upper/lower seismo depths from Domains usd, lsd = get_ul_seismo_depths(src_codes, usd, lsd) ############################################################################### # make some adjustments ###############################################################################
############################################################################### # get preferred catalogues ############################################################################### # get preferred catalogues for each zone prefCat = get_preferred_catalogue(domshp) # fix catalogue for source zones prefCat[55] = 'NSHA18CAT_V0.1_hmtk_declustered.csv' prefCat[56] = 'NSHA18CAT_V0.1_hmtk_declustered.csv' ############################################################################### # load 2018 completeness models ############################################################################### single_mc = 0 ycomp, mcomp, min_rmag = get_completeness_model(src_codes, shapes, domains, single_mc) # use manual modification for i in range(0, len(trt_new)): if trt_new[i] == 'Active': min_rmag[i] = 5.7 min_rmag[3] = 6.1 # TAFS min_rmag[12] = 6.1 # NBT min_rmag[11] = 6.0 # NBOT min_rmag[16] = 5.6 # BNBD min_rmag[26] = 3.5 # NWO min_rmag[50] = 3.2 # CARP min_rmag[51] = 3.5 # EAPM min_rmag[52] = 3.3 # KMBY min_rmag[55] = 3.3 # NACR
############################################################################### # get preferred catalogues ############################################################################### # get preferred catalogues for each zone prefCat = get_preferred_catalogue(outshp) ''' # fix catalogue for source zones prefCat[44] = 'NSHA18CAT_V0.1_hmtk_declustered.csv' prefCat[2] = 'NSHA18CAT_V0.1_hmtk_declustered.csv' ''' ############################################################################### # load 2018 completeness models ############################################################################### single_mc = 0 ycomp, mcomp, min_rmag_ignore = get_completeness_model(src_codes, shapes, neo_domains, single_mc) #ycomp, mcomp, min_rmag_ignore = get_completeness_model_point(src_codes, shapes, single_mc) # use values from Domains model instead min_rmag = neo_min_rmag # use manual modification for i in range(0, len(neo_trt)): if neo_trt[i] == 'Active': min_rmag[i] = 5.75 elif neo_trt[i] == 'Intraslab': min_rmag[i] = 5.75 elif neo_trt[i] == '': min_rmag[i] = 5.75 ''' min_rmag[23] = 6.1 # TAFS