domains[i] = 11
if trt[i] == 'NCratonic':
trt_new.append('Non_cratonic')
else:
trt_new.append(trt[i])
#trt_new[119] = 'Oceanic'
###############################################################################
# load neotectonic domains parameters
###############################################################################
# set domestic domain numbers based on neotectonic domains
refshpfile = '..//reference_shp.txt'
#refshpfile = '..//reference_shp_mx.txt' # for testing only!
neo_domains, neo_min_rmag, neo_mmax, neo_trt, neo_bval_fix, neo_bval_sig_fix = get_neotectonic_domain_params(sf, trt_new, refshpfile)
# set b-values and sigmas
bval_fix = neo_bval_fix
bval_sig_fix = neo_bval_sig_fix
# set b-values for zones with weird centroids
#bval_fix[28] = bval_fix[54] # NWO
bval_fix[52] = bval_fix[53]
bval_fix[55] = bval_fix[53]
bval_fix[44] = bval_fix[19]
bval_fix[45] = bval_fix[19]
bval_fix[47] = bval_fix[25]
bval_fix[85] = bval_fix[109]
bval_fix[117] = bval_fix[93]
domains[i] = 10
elif trt[i] == 'Intraslab':
domains[i] = 11
if trt[i] == 'NCratonic':
trt_new.append('Non_cratonic')
else:
trt_new.append(trt[i])
###############################################################################
# load neotectonic domains parameters
###############################################################################
# set domestic domain numbers based on neotectonic domains
neo_domains, neo_min_rmag, neo_mmax, neo_trt, neo_bval_fix, neo_bval_sig_fix \
= get_neotectonic_domain_params(sf, trt_new)
for i in range(0, len(domains)):
if neo_domains[i] > 0 and neo_domains[i] < 8:
domains[i] = neo_domains[i]
mmax[i] = neo_mmax[i]
zone_class = list(domains)[:]
'''
# reset Gawler Craton to Flinders due to b-value similarities
zone_class[0] = 2.
# reset West Coast to extended
zone_class[7] = 6.
domains[7] = 6.