<div dir="ltr"><div>Actually, if there are no observations within a specified grid box, a chek for _FillValue should be added.<br></div>Also, an attribute specifying the number of observations used should be added.<br><div><br><span class="gmail-im">itime = <b>ind</b>(LAT.ge.lat_low .and. LAT.lt.lat_high .and. \<br> LON.ge.lon_low .and. LON.lt.lon_high)<br><br><div>if (<b>.not.</b>(<b>ismissing</b>(itime(0)) then<br></div><div> T_super = <b>dim_avg_n_Wrap</b>(T(:,itime), 1)<br></div></span><div> T_super@long_name = "Super Ob: " + T@long_name ; cancatenate<br></div><div> <b>T_super@N_obs = dimsizes(itime) </b><br></div><div> printVarSummary(T_super)<br></div><div>end if<br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 27, 2018 at 8:18 AM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><span class=""><div><div><div>Observations: T(lev,time) and Q(lev,time), LAT(time), LON(time)<br><br></div></div>Is the target grid rectilinear?<br><br></div><div>If, say, the grid box is bounded by lat_low, lat_high, lon_low, lon_high<br></div><div><br></div> itime = <b>ind</b>(LAT.ge.lat_low .and. LAT.lt.lat_high .and. \<br> LON.ge.lon_low .and. LON.lt.lon_high)<div><br></div><div> T_super = <b>dim_avg_n_Wrap</b>(T(:,itime), 1)<br></div></span><div> T_super@long_name = "Super Ob: " + T@long_name ; cancatenate<br></div><div> printVarSummary(T_super)<br></div><div><br></div><div>Same for Q<br><br>===<br><b><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Document/Functions/Built-in/<wbr>ind.shtml</a></b><br>===<br></div><span class=""><div>For grids the lat_low, lat_high, lon_low, lon_high would have to be indexed.<br></div><div><br></div>HTH</span></div><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Tue, Jun 26, 2018 at 12:44 PM, Andrew Kren - NOAA Affiliate <span dir="ltr"><<a href="mailto:andrew.kren@noaa.gov" target="_blank">andrew.kren@noaa.gov</a>></span> wrote:<br></span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Dear ncl-talk,<div><br></div><div>I am thinning some data prior to assimilating it in the HWRF model. The data is temperature and moisture vertical profiles that is one-dimensional, a function of time. The lat/lon points of the data is also a function of time, so there is a spatial element to it. Basically, I want to superob the data to create a mean observation within some grid box size to thin it. I would need to regrid the data and search how many profiles within the specified grid size. My question is what would be the best way to achieve this via NCL?</div><div><br></div><div>Thanks,<span class="m_7966764038900498243HOEnZb"><font color="#888888"><br clear="all"><div><br></div>-- <br><div class="m_7966764038900498243m_-9080444856367550076gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div class="m_7966764038900498243m_-9080444856367550076gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div>Andrew Kren<br>Assistant Scientist<br>CIMAS - NOAA/AOML<br><span><font color="#888888">314-322-0867</font></span></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</font></span></div></div>
<br></div></div><span class="">______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></span></blockquote></div><br></div>
</blockquote></div><br></div>