Add region argument into plot set 3, including adding to function reduce2lat_seasonal