I've spent a lot of time trying to make a radial (polar) map to plot the 500 hPa Geopotential data for Antarctica (previous post: Is it possible to make a polar (radial) filled.contour map? Geopotential in Antarctica data).
I finally discovered that I can accomplish this with the wonderful function I didn't know about,
PlotStereoMap, from the
s2dverification package. Here my code and map:
PlotStereoMap(verano2008_med, lon, lat, latlims = c(-30,-90), toptitle = "SAM (+)", title_scale = 0.7, bar_limits = range(verano2008_med), triangle_ends = c(T,T), units = "m2/s2", units_scale = 0.7, bar_label_scale = 0.7, bar_tick_scale = 0.4)
However, I am now struggling with changing the plot parameters. I've read them all and tried everything, but I can't get the color scale to be smaller (it's almost the size of the map!) and I'd also like it to be horizontal at the bottom of the map.