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", "stream1e") ; Open "stream1e.png" png file. res = True ; Plot options desired. res@gsnPolarNH = True ; Specify the hemisphere res@mpMinLatF = 60 ; Minimum lat to plot res@stArrowLengthF = 0.008 ; Default is dynamic res@stLengthCheckCount = 5 ; default is 35 res@stArrowStride = 1 ; Default is 2 res@stLineStartStride = 1 ; default is 2 res@stMinArrowSpacingF = 0.03 ; Default is 0.0 res@stStepSizeF = 0.001 ; default is dynamic res@stMinDistanceF = 0.03 ; Distance between lines res@stMinLineSpacingF = 0.005 plot = gsn_csm_streamline_contour_map_polar(wks,u,v,t,res) end