Bug fixes in calculate_seasonal_climatology
Retain units in select_region
Add region argument into plot set 9, including adding to reduce2latlon_seasonal

Add region argument into plot set 9
- add region functionality to function reduce2latlon_seasonal
- also define seasonsyr in amwg.py and set as default season