<html><head></head><body><div class="ydp8afb1634yahoo-style-wrap" style="font-family:times new roman, new york, times, serif;font-size:16px;"><div></div>
<div><br></div><div>Dear <span>Stavros NTAFIS</span></div><div><span>thanks for your help</span><br></div><div>but problem is how plot for four time step and in vertical profile(different height or pressure value)</div><div>same as attach file</div><div><br></div><div><br></div>
</div><div id="yahoo_quoted_5417935981" class="yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>
On Saturday, December 15, 2018, 2:47:26 PM GMT+3:30, sdafis <sdafis@noa.gr> wrote:
</div>
<div><br></div>
<div><br></div>
<div><div id="yiv3534705404"><div>
<p>Dear Elham, I am giving you a simple code to read and plot a time series for T2m and 10-m wind speed for two locations defined by "locs", "lats", "lons". You can add more locations if you like or choose a sub-period by defining "it".</p>
<p><br clear="none"></p>
<p><br clear="none"></p>
<p><br clear="none"></p>
<p>a = addfile("wrf_file.nc","r") </p>
<p>; Get variables for 10-m wind</p>
<p>uvmet = wrf_user_getvar(a, "<em>uvmet10</em>", -1)</p>
<p>u = uvmet(0,:,:,:) <br clear="none">v = uvmet(1,:,:,:) </p>
<p>; Get Temp at 2m in Kelvin</p>
<pre> tc2 = <strong>wrf_user_getvar</strong>(a,"T2",-1) </pre>
<p>it = 0 ; first time step, ignore it if you want a time-series</p>
<p>locs = (/ "1st Location", "2nd Location" /) ; define the locations </p>
<p>lats = (/ 37.5938, 40.6401 /) <br clear="none">lons = (/ 23.4375, 22.9444/)</p>
<p><br clear="none">nlocs = dimsizes(locs) <br clear="none">loc = wrf_user_ll_to_ij(a, lons, lats, True) ; find the indices of the closest WRF grid points to locations</p>
<p>do iloc = 0, nlocs-1<br clear="none"><br clear="none">print("Working for " + locs(iloc))</p>
<p>locX = loc(0, iloc) - 1<br clear="none">locY = loc(1, iloc) - 1</p>
<p>wind= sqrt(u(:, locY, locX)*u(:, locY, locX) + v(:, locY, locX)*v(:, locY, locX)) ; replace ":" with "it" if you want one time-step only</p>
<p><span>temp = tc2(:, locY, locX) ; if you want only one time step: temp = tc2(it, locY, locX)</span></p>
<p><br clear="none"></p>
<p>print( ""+wind+" "+temp+"")</p>
<p><br clear="none"></p>
<p>end do</p>
<p><br clear="none"></p>
<p><br clear="none"></p>
<p>regards,</p>
<p><br clear="none"></p>
<p>-- </p>
<div class="yiv3534705404pre" style="margin:0;padding:0;font-family:monospace;">Stavros NTAFIS (DAFIS)<br clear="none">-----------------------------------------------<br clear="none">Physicist - Meteorologist, M.Sc.<br clear="none">Ph.D. Candidate, Polytechnic School of Paris, France<br clear="none">Laboratory of Dynamic Meteorology (LMD)<br clear="none">Research Associate, National Observatory of Athens (NOA/IERSD)<br clear="none">ORCID: <a rel="nofollow" shape="rect" target="_blank" href="https://orcid.org/0000-0002-1513-1930">https://orcid.org/0000-0002-1513-1930</a><br clear="none">Tel. : (+33)1 69 33 51 64<br clear="none">(+33)9 81 94 22 12<br clear="none">Mobile: (+33) 066 030 0147<br clear="none">(+30) 697 04 20 242<br clear="none">---------------<br clear="none">Weather charts:<br clear="none"> <a rel="nofollow" shape="rect" target="_blank" href="http://www.meteo.gr">http://www.meteo.gr</a><br clear="none"> <a rel="nofollow" shape="rect" target="_blank" href="http://www.meteo.gr/meteomaps/">http://www.meteo.gr/meteomaps/</a></div>
<p><br clear="none"></p>
<div class="yiv3534705404yqt7619235638" id="yiv3534705404yqtfd24100"><p>Στις 2018-12-15 11:54, E Mobarak έγραψε:</p>
<blockquote type="cite" style="padding:0 0.4em;border-left:#1010ff 2px solid;margin:0;">
<div class="yiv3534705404ydp196b94c9yahoo-style-wrap" style="font-family:times new roman, new york, times, serif;font-size:16px;">Hi dear ncl usr</div>
<div class="yiv3534705404ydp196b94c9yahoo-style-wrap" style="font-family:times new roman, new york, times, serif;font-size:16px;">I want plot some data for example temperature or wind speed in some time step from wrf out data file in one location</div>
<div class="yiv3534705404ydp196b94c9yahoo-style-wrap" style="font-family:times new roman, new york, times, serif;font-size:16px;">how i should it?</div>
<div class="yiv3534705404ydp196b94c9yahoo-style-wrap" style="font-family:times new roman, new york, times, serif;font-size:16px;">Best Regarads</div>
<div class="yiv3534705404ydp196b94c9yahoo-style-wrap" style="font-family:times new roman, new york, times, serif;font-size:16px;">Elham</div>
<br clear="none">
<div class="yiv3534705404pre" style="margin:0;padding:0;font-family:monospace;">_______________________________________________<br clear="none"> ncl-talk mailing list<br clear="none"> <a rel="nofollow" shape="rect" ymailto="mailto:ncl-talk@ucar.edu" target="_blank" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br clear="none"> List instructions, subscriber options, unsubscribe:<br clear="none"> <a rel="nofollow" shape="rect" target="_blank" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></div>
</blockquote>
<p><br clear="none"></p>
</div></div></div></div>
</div>
</div></body></html>