Hi Xiangdong, You can do this kind of thing, if your variables are v1, lon, lat: ! Run the mp_stereographic_north script to define projection variables set region/y=60:90N go mp_stereographic_north ! Reset the map projection coordinates to use the lon, lat grid let/quiet mp_lambda = lon * deg2rad let/quiet mp_phi = lon* deg2rad set grid lon go mp_aspect shade/noaxes v1, x_page, y_page Xiangdong Zhang wrote: Thanks. In this case, both plots are on the rectangular map. I hope the plot will be on the polar stereographic projection map. More suggestions? Thanks again -- XiangdongYou can use Ferret's curvilinear mapping procedures for this. Let's suppose you want to contour variable v1 which is on a rectangular i,j grid that is not necessarily a global lat,lon grid, but you also have lat and lon on that same rectangular grid. Then the following will do the job contour v1,lon,lat Your other data set with variable v2, already on a lon/lat grid, is contoured by contour v2 ------------------------------------------------------------------------ ------------------ Edward D. (Ned) Cokelet, Ph.D. Edward.D.Cokelet@noaa.gov NOAA/PMEL off: (206) 526-6820 7600 Sand Point Way NE fax: (206) 526-6485 Seattle, WA 98115-6439 |