I've noticed that plot_surfaces is noticeably slower than it was in the mid-past (a few months ago), when the vartheta contours are plotted. Wonder if a change to how those are calculated is the cause (like in map coords or something)