begin a = addfile("pop.nc","r") ; Read in POP data. u = a->urot v = a->vrot u@lat2d = a->lat2d u@lon2d = a->lon2d v@lat2d = a->lat2d v@lon2d = a->lon2d wks = gsn_open_wks("png", "stream2d") ; Open "stream2d.png" png file. res = True ; Plot options desired. res@gsnAddCyclic = True ; Add longitude cyclic points. res@stLineColor = "brown" ; Brown streamlines res@stMinDistanceF = 0.021 ; Thins streamlines. plot = gsn_csm_streamline_map(wks,u,v,res) ; Call gsn_csm routine for ; drawing streamlines over map. end