<div dir="ltr"><div><br></div>Hi Tabish,<div><br></div><div>To get the time series for uvmet, you can supply the wrf_user_getvar function a time index of -1, and that will extract the uvmet variable across all times. For example:</div><div><br></div><div>files = systemfunc("ls -l wrfout_d02_2005*") + ".nc"</div><div>a = addfiles(files, "r")</div><div>uvmet_all = wrf_user_getvar(a, "uvmet", -1)</div><div><br></div><div>This will return an array of shape: u_v x Time x bottom_top x south_north x west_east .</div><div><br></div><div>I am assuming here that you know the south_north and west_east indexes that you want, so to get the time series of vertical columns, you would do:</div><div><br></div><div>uvmet_series = uvmet_all(:, :, :, your_south_north, your_west_east)</div><div><br></div><div>This will return an array of shape: u_v x Time x bottom_top</div><div><br></div><div>We're almost there. Now you just need to use the wrf_interp_1d function (<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_interp_1d.shtml">http://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_interp_1d.shtml</a>) to interpolate each column to the height levels you want at each time. (Unfortunately, this function does not implicitly loop across all times, so you need to manually loop for each time). </div><div><br></div><div>First, we need the model heights for all of the times.</div><div><br></div><div>heights = wrf_user_getvar(a, "z", -1)</div><div>height_series = heights(:, :, your_south_north, your_west_east)</div><div><br></div><div>Now, get the number of times:</div><div><br></div><div>dim_ht = dimsizes(height_series)</div><div>ntimes = dim_ht(0)</div><div><br></div><div>Create an array with your desired height levels to interpolate to (in meters):</div><div><br></div><div>your_heights = (/100., 1000., 2000./)</div><div><br></div><div>Finally, let's loop over all times and get each interpolated column:</div><div><br></div><div>do t = 0,ntimes - 1</div><div> u_interp = wrf_interp_1d(uvmet_series(0, t, :), height_series(t, :), your_heights)</div><div> v_interp = wrf_interp_1d(uvmet_series(1, t, :), height_series(t, :), your_heights)</div><div> ; Now do something with these variables</div><div>end do</div><div><br></div><div>If you are just looking to store the data, prior to the loop, you can create an empty array (using new) that is 2 x ntimes x dimsizes(your_heights) and then copy the u_interp and v_interp to this array inside the loop.</div><div><br></div><div>Hope this helps,</div><div><br></div><div>Bill</div><div><br></div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Mar 10, 2017 at 11:26 AM, Tabish Ansari <span dir="ltr"><<a href="mailto:tabishumaransari@gmail.com" target="_blank">tabishumaransari@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi<br><br></div><div>I wish to plot time-series of vertical profile of wind (u and v) at a given location (grid point) from my WRF output. The closest script I've got for my purpose is this:<br><br><a href="http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/SPECIAL/wrf_meteo_4.ncl" target="_blank">http://www2.mmm.ucar.edu/wrf/<wbr>OnLineTutorial/Graphics/NCL/<wbr>Examples/SPECIAL/wrf_meteo_4.<wbr>ncl</a><br><br></div><div>However, this script plots model levels on the y-axis - I want height in metres instead. Further, this script works over a single wrfout file which is imported using addfile but I have multiple daily files. I generally import variables like: temp = a[:]->t2 but this doesn't work with wrf_user_getvar() method. So I need an alternate way to read in umet and vmet (as in the script above) where I get rotated wind vectors time-series from all the files in order.<br><br></div><div>Could you please help me achieve this?<br><br></div><div>Thanks in advance.<br><br></div><div>Tabish<br></div><div><br clear="all"></div><div><div class="gmail-m_2459414188924563755gmail_signature"><div dir="ltr"><div><div><font size="1"><span style="font-family:tahoma,sans-serif">Tabish U Ansari<br></span></font></div><font size="1"><span style="font-family:tahoma,sans-serif">PhD student, Lancaster Environment Center<br></span></font></div><font size="1"><span style="font-family:tahoma,sans-serif">Lancaster Univeristy<br> <span>Bailrigg</span>, <span>Lancaster</span>, <br><span>LA1 4YW</span>, <span>United Kingdom</span></span></font><br></div></div></div>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></div>