begin
ncurve = 3
npts = 129
;
; Read ASCII file xy.asc
;
x1_y3 = asciiread(ncargpath("data")+"/asc/xy.asc",4*129,"float")
;
; xy.asc has 4 vars of length 129 longitudes, lon, u, v, t
;
; The data is taken at 43N latitude. Longitude is an index
; 1-129 standing for 0 deg - 360 deg in steps of 360/128
; u and v are in m/s, and t is in deg K.
;
lon = (x1_y3(0:512:4) -1.) * 360./128.
u = x1_y3(1:513:4)
v = x1_y3(2:514:4)
t = x1_y3(3:515:4)
wks = gsn_open_wks("png","xy5a")
res1 = False ; No resources desired
res2 = False ; for any curve at this point.
res3 = False
plot = gsn_csm_xy3(wks,lon,t,u,v,res1,res2,res3)
end