Add region argument into plot set 3, including function reduce2lat_seasonal.