<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Imran,<div class="">unfortunately I’m not really an expert of WRF, so I’m going to give you a general hint which may work. You should be more precise next time as generally people do not know what is LU_INDEX and what is its structure ;) </div><div class=""><br class=""></div><div class="">I assume that you want to extract temperature, which is 3-D (time, lat, lon) using a mask array LU_INDEX which is 2-D (lat, lon). Again, if this assumption is wrong please be more specific. You have many choices, but I believe that for the sake of clarity you may want to go for an array which has the same dimensions but missing values where LU_INDEX is not 2. So you may want to do something like this:</div><div class=""><br class=""></div><div class=""><font face="Menlo" class="">T2_masked=mask(T2, LU_INDEX.ne.2, False).</font></div><div class=""><br class=""></div><div class="">You need to have a Missing Value attribute set for the T2 variable of course. As a result you will get an array with the same dimensions, so that it’s easier to compare with the original one, but with Missing Values where your condition is False, so that for every computation these points will not be considered. You can achieve something similar with the where function. Please have a look here (<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/mask.shtml" class="">http://www.ncl.ucar.edu/Document/Functions/Built-in/mask.shtml</a>).</div><div class=""><br class=""></div><div class="">Cheers</div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">Il giorno 14 mar 2017, alle ore 04:38, Imran Hosen <<a href="mailto:hosen.imran09@gmail.com" class="">hosen.imran09@gmail.com</a>> ha scritto:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class="">Dear NCL users,<br class=""><br class=""></div>Good day!<br class=""><br class=""></div>I'm a new user of NCL. I want to extract temperature (T2) from wrfout file for the specific grid points where the LU_INDEX = 2.<br clear="all" class=""><div class=""><div class=""><div class=""><br class=""><div class="gmail_signature" data-smartmail="gmail_signature"><div class="">For example,<br class=""></div><div class="">f = addfile("./<a href="http://wrfout.nc/" class="">wrfout.nc</a>)<br class=""></div><div class="">T2 = f->T2(:,:,:)<br class=""><br class=""></div><div class="">It will extract temperature for whole domain for all time steps. In my case, I only want to extract temperature for all time steps but for the specific grids which grids represent LU_INDEX = 2.<br class=""><br class=""></div><div class="">Could you please suggest how to do it. I would grateful to you if anyone please share a script.<br class=""><br class=""></div><div class="">Thanks in advance.<br class=""><br class=""></div><div class="">Kind regards,<br class=""></div><div class="">Imran<br class=""></div><div class=""><br class=""></div></div>
</div></div></div></div>
_______________________________________________<br class="">ncl-talk mailing list<br class=""><a href="mailto:ncl-talk@ucar.edu" class="">ncl-talk@ucar.edu</a><br class="">List instructions, subscriber options, unsubscribe:<br class="">http://mailman.ucar.edu/mailman/listinfo/ncl-talk<br class=""></div></blockquote></div><br class=""></div></body></html>