Delete set-specific method interpret_region, which is replaced by shared function in reductions.py
Insert select_region and interpret_region functions where appropriate in other reduction.py functions
Add default arguments to reduce_time_space_seasonal_regional
Refactor some code with generic functions and add delete_singleton_axis to select_lev function.
Refactor seasonal climatology functions -- try to push all special cases into one function.