begin a = addfile("atmos2.nc","r") ; Read in first time step t = a->TS(0,:,:) ; of atmospheric data. wks = gsn_open_wks("png", "contour3b") ; Open png file "contour3b.png" res = True ; Plot mods desired res@gsnPolar = "NH" ; Specify northern hemisphere ; (which is the default). res@cnFillOn = True ; Turn on contour fill res@cnFillPalette = "wgne15" plot = gsn_csm_contour_map_polar(wks,t,res) ; Call gsn_csm function for ; drawing contours over a ; polar map. end