<html><head></head><body><div class="ydp8ca1bb6fyahoo-style-wrap" style="font-family:times new roman, new york, times, serif;font-size:16px;"><span><span class="ydpbfa26418tlid-translation ydpbfa26418translation"><span title="">Hello</span><br><span title="">I've plotted the wind profile with this script</span><br><span title="">But y-axis is up to 21,000 meters </span><br><span title="">How can I limit the vertical axis and draw only the first 5 km?</span><br></span></span><div><span><span class="ydpbfa26418tlid-translation ydpbfa26418translation"><span title="">Thanks for answer</span></span></span></div><div><span><span class="ydpbfa26418tlid-translation ydpbfa26418translation"><span title=""><br></span></span></span></div><div><span><span class="ydpbfa26418tlid-translation ydpbfa26418translation"><span title=""><span> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"<br><br><br><br> a1 = addfile(wrfout_d01_2009-02-19_00:00:00.nc,"r")<br><br> wks = gsn_open_wks ("png","profile") ; send graphics to PNG file<br><br> u1 = wrf_user_getvar(a1,"ua",-1) ; u on mass points<br> v1 = wrf_user_getvar(a1,"va",-1) ; v on mass points<br> z1 = wrf_user_getvar(a1,"z",-1) ; u on mass points<br><br> res = True<br> plots = new(4,graphic)<br> llres = True<br> llres@ReturnInt = True ;<br><br> locs = (/"Ahvaz","bu"/)<br> lats = (/31.34,30/)<br> lons = (/44.53,41/)<br><br></span></span></span></span><div><span><span class="ydpbfa26418tlid-translation ydpbfa26418translation"><span title=""><span> nlocs = dimsizes(locs)</span></span></span></span></div><div><span><span class="ydpbfa26418tlid-translation ydpbfa26418translation"><span title=""><span></span></span></span></span><br></div><div><span>loc = wrf_user_ll_to_ij(a1, lons, lats, True)<br><br> do iloc = 0, nlocs-1<br> locX = loc(0, iloc) - 1<br> locY = loc(1, iloc) - 1<br><br> times = wrf_user_list_times(a1) ; get times in the file<br> ntimes = dimsizes(times) ; number of times in the file<br><br> it = 10<br> zz1 = z1(it,:,locY,locX)<br> uu1 = u1(it,:,locY, locX)<br> vv1 = v1(it,:,locY, locX)<br> spd1 = (uu1*uu1 + vv1*vv1)^(0.5)<br>;------------------ res -------------------------<br> res = True<br> res@gsnDraw = False ; Don't draw individual plot.<br> res@gsnFrame = False ; Don't advance fram<br>;-----------------------------------------------------<br> wks = gsn_open_wks ("png","xy-"+locs(iloc)) ; send graphics to PNG file'<br> uu_plot = gsn_csm_xy(wks,spd1,zz1,res)<br></span><div><span> draw(uu_plot)</span></div><div><span></span><br></div><div>frame(wks)<br></div><span>end do<br></span><br></div><span><span class="ydpbfa26418tlid-translation ydpbfa26418translation"><span title=""><span><br></span><br></span></span></span></div></div></body></html>