; ; This file has 2D lat/lon arrays associated with the data. We need to ; attach them as special attributes "lat2d" and "lon2d". ; begin a = addfile("pop.nc","r") ; Read in POP data. u = a->urot v = a->vrot u@lat2d = a->lat2d ; This file has 2D lat/lon arrays u@lon2d = a->lon2d ; associated with the data. v@lat2d = a->lat2d ; We need to attach them as special v@lon2d = a->lon2d ; attributes "lat2d" and "lon2d". wks = gsn_open_wks("png", "stream2f") ; Open "stream2f.png" png file. res = True ; Plot options desired. res@gsnMaximize = True ; Maximize plot in frame. 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