Compute climatologies for only one season at a time, doing each season separately from the others.
This greatly reduces the memory footprint and hence greatly improves speed for high-resolution data.
Also I shrank the default season list to the three most important seasons.