begin f = addfile ("atmos2.nc", "r") ; Open data file and pull off U, T U = f->U(0,0:1,{30},:) ; data at the first two T = f->T(0,0:1,{30},:) ; levels and latitude = 30. wks = gsn_open_wks("png","xy4b") ; Open a png file called "xy4b.png". resU = True ; Plot options desired resT = True ; Plot options desired resU@trXMinF = min(U&lon) ; Limits for X axis resU@trXMaxF = max(U&lon) plot = gsn_csm_xy2(wks,U&lon,U,T,resU,resT) ; Call gsn_csm_xy2 to generate ; XY plot with two Y axes. end