Hi Geeta, Here's what I do to do that kind of plot - even puts the color bar at the bottom, as a bonus :) Cheers, Paul use "z500_djf.nc" !e.g. 500 mb geopotential height define symbol lv = (-inf)(-30,30,5)(inf) !Some plotting options let latS = -40 let size = 0.11 set region/y=-90:`latS` !Plot set up set region/x=-180:180 set grid z500[d=1] go mp_stereographic_south 180 fill/lev=($lv)/noaxes/set/x=-180:180 \ z500[d=1], x_page, y_page !Do the plot ppl title " " let xlo = ($ppl$xorg) - 1.5 let xhi = `xlo` + ($ppl$xlen) + 3 let ylo = ($ppl$yorg) - 0.6 let yhi = ($ppl$yorg) - 0.6 + 0.12 ppl shakey `dokey[i=($xxx)]`,0,`size`,2,-0,,`xlo`,`xhi`,`ylo`,`yhi` ppl fill go mp_land 7 !Continents (heavy black) go mp_graticule " " " " " " `latS` `latS` 5 13 !Outline the whole circle On Mar 20, 2012, at 10:21 AM, Geeta Persad wrote: Hi All, |