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","xy5b")
res1 = True ; Plot options desired.
res2 = True
res3 = True
res1@xyLineThicknessF = 3 ; Set line thickness and
res1@xyLineColor = "red" ; curve colors for each set
; of curves.
res2@xyLineThicknessF = 3
res2@xyLineColor = "green"
res3@xyLineThicknessF = 3
res3@xyLineColor = "blue"
res1@trXMaxF = 360. ; X axis is the same for all
; curves, so use the first resource
; list to control X axis.
plot = gsn_csm_xy3(wks,lon,t,u,v,res1,res2,res3)
end