load topo contour(0:359,-89:90,topo,[0 0],'b') axis equal box on set(gca,'XLim',[120 240],'YLim',[20 70], ... 'XTick',[120 150 180 210 240], ... 'Ytick',[20 45 70]);