begin f = addfile("HGT500_MON_1958-1997.nc","r") hgt = f->HGT wks = gsn_open_wks("png","contour8a") ; Open "contour8a.png" png file res = True ; No plot options desired plot_base = gsn_csm_contour_map_polar(wks,hgt(0,:,:),res) end