Code at the regional level now working
Update to compute aggregates in python