begin a = addfile("atmos2.nc","r") ; Read in first time step u = a->U(0,15,:,:) ; and 15th level of v = a->V(0,15,:,:) ; atmospheric data. t = a->T(0,15,:,:) wks = gsn_open_wks("png", "stream1c") ; Open "stream1c.png" png file. res = True ; Plot options desired. res@gsnPolarNH = True ; Specify the hemisphere res@mpMinLatF = 60 ; Minimum lat to plot res@stMinDistanceF = 0.03 ; Distance between lines res@stMinLineSpacingF = 0.005 plot = gsn_csm_streamline_contour_map_polar(wks,u,v,t,res) end