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