# Summarize utility_of_residents_zone attributes
swf_per_year = zones.attribute_sum('utility_of_residents_zone')
elif case_study == 'zurich' and policy_level == 'parcel':
# export computations to tab file [parcel level]
if export_computations_parcel == True:
print '[Parcel Level] Exporting computations to: %s' % out_table_name_computations
parcels.write_dataset(attributes = ['income_per_parcel', 'housing_cost_per_parcel', 'travel_cost_per_parcel', 'utility_of_residents_parcel'],
out_storage = storage_output,
out_table_name = out_table_name_computations)
# Summarize utility_of_residents_zone attributes
swf_per_year = parcels.attribute_sum('utility_of_residents_parcel')
print 'Social Welfare for year %s: %s' % (year, swf_per_year)
# Place swf idicator value for the running year into a dict storage
# Add swf as primary attribute for the running year
swf_indicators.add_elements(data = {"swf_id":array([1]),
"year":array([year]),
"social_welfare":array([swf_per_year])}, require_all_attributes = True, change_ids_if_not_unique = True)
# Summarize social_welfare attributes for all years
swf_summary = swf_indicators.attribute_sum('social_welfare')
print '-------------------------------------'
print 'Total Social Welfare: %s' % swf_summary
# Add swf_summary as primary attribute to table swf_indicators
swf_indicators.add_elements(data = {"swf_id":array([1]),
"year":array(['swf_summary']),
"social_welfare":array([swf_summary])}, require_all_attributes = True, change_ids_if_not_unique = True)
# export social welfare indicators to tab file
print 'Exporting Table of Social Welfare Indicators to: %s' % out_table_name_swf_indicators
swf_indicators.write_dataset(attributes = ['swf_id', 'year', 'social_welfare'],
out_storage = storage_output,
out_table_name = out_table_name_swf_indicators)
