create_aggregate_values.py

import ogr
path=r"C:\Users\Rich\Documents\seasonal_water_yield_workspace\aggregated_results_cz.shp"
v=ogr.Open(path)
l=v.GetLayer()
with open('agg_results_cz.csv', 'wb') as out_csv:
    for feature in l:
        line = str(feature.GetFID())
        line += ','+str(feature.GetField('vri_sum'))
        line += ','+str(feature.GetField('qb'))+'\n'
        out_csv.write(line)