src_code = src_code[:3] + '_' + src_code[3:]
new_src_codes.append(src_code)
# now match zone class based on lookup table
zone_class[i] = aggregate_intraslab_sources(src_code, zone_class[i])
# supplant new source codes
del src_codes
src_codes = array(new_src_codes).copy()
###############################################################################
# load Rajabi SHMax vectors
###############################################################################
shmax_pref, shmax_sig = get_aus_shmax_vectors(src_codes, shapes)
###############################################################################
# get rate adjustment factors
###############################################################################
origshp = 'Domains_Sep2011_edit.shp'
newField = 'code'
origField = 'CODE'
rte_adj_fact = get_rate_adjust_factor(domshp, newField, origshp, origField)
###############################################################################
# write initial shapefile
###############################################################################
if single_mc == 1:
outshp = 'Domains_NSHA18_single_Mc.shp'
dy,
head_width=30000,
head_length=head_length,
lw=0.3,
fc='0.2',
fill=True,
length_includes_head=True)
##########################################################################################
# get average per zone and plot
##########################################################################################
src_codes = get_field_data(sfz, 'CODE', 'str')
src_shapes = sfz.shapes()
av_shmax, sig_shmax = get_aus_shmax_vectors(src_codes, src_shapes)
# get centroid and plot
alen = 180000
head_length = 60000
for i, shape in enumerate(src_shapes):
centroid = Polygon(shape.points).centroid.wkt
centroid = centroid.strip('PIONT').replace('(', ' ').replace(')',
' ').split()
cx, cy = m(float(centroid[0]), float(centroid[1]))
dy = alen * cos(radians(av_shmax[i])) + head_length * cos(
radians(av_shmax[i]))
dx = alen * sin(radians(av_shmax[i])) + head_length * sin(
radians(av_shmax[i]))