<div dir="ltr">NCL support,<div><br></div><div>Is there any progress on this request?  This issue keeps coming up in my work group.  We perform a lot of generalized file handling.  Thank you.</div><div><br></div><div>--Dave<br><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 22, 2013 at 12:30 PM, David Brown <span dir="ltr">&lt;<a href="mailto:dbrown@ucar.edu" target="_blank">dbrown@ucar.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Dave,<br>
This is a very good suggestion. We will create a ticket for it.<br>
 -dave<br>
<br>
On Aug 22, 2013, at 12:06 PM, Dave Allured - NOAA Affiliate &lt;<a href="mailto:dave.allured@noaa.gov">dave.allured@noaa.gov</a>&gt; wrote:<br>
<br>
&gt; NCL support,<br>
&gt;<br>
&gt; This is a request for improved support for mixed types in coordinate<br>
&gt; subscripting.  For example, it is legal to subscript a double<br>
&gt; coordinate with floats, but going the other way is an error.  This<br>
&gt; shows that type handling in simple NCL relational operators is more<br>
&gt; liberal than in coordinate subscripting:<br>
&gt;<br>
&gt; lats = (/ 0., 10., 20., 30., 40., 50., 60. /)<br>
&gt; x    = (/ 0,  111, 222, 333, 444, 555, 666 /)<br>
&gt; x!0  = &quot;lat&quot;<br>
&gt; y    = x<br>
&gt;<br>
&gt; x&amp;lat = tofloat  (lats)<br>
&gt; y&amp;lat = todouble (lats)<br>
&gt;<br>
&gt; f1 = 20.     ; floats<br>
&gt; f2 = 40.<br>
&gt; d1 = 20d     ; doubles<br>
&gt; d2 = 40d<br>
&gt;<br>
&gt; ind1 = ind (y&amp;lat .ge. f1 .and. y&amp;lat .le. f2)<br>
&gt; print (ind1+&quot;&quot;)<br>
&gt; print (&quot;&quot;)<br>
&gt;<br>
&gt; ind2 = ind (x&amp;lat .ge. d1 .and. x&amp;lat .le. d2)<br>
&gt; print (ind2+&quot;&quot;)<br>
&gt; print (&quot;&quot;)<br>
&gt;<br>
&gt; sub1 = y({f1:f2})    ; float subscripts double<br>
&gt; print (sub1&amp;lat+&quot; &quot;+typeof(sub1&amp;lat))<br>
&gt; print (&quot;&quot;)<br>
&gt;<br>
&gt; sub2 = x({d1:d2})    ; double subscripts float<br>
&gt; print (sub2&amp;lat+&quot; &quot;+typeof(sub2&amp;lat))<br>
&gt;<br>
&gt; Result:<br>
&gt;<br>
&gt; NCAR Command Language Version 6.1.2<br>
&gt; (0)   2<br>
&gt; (1)   3<br>
&gt; (2)   4<br>
&gt; (0)<br>
&gt; (0)   2<br>
&gt; (1)   3<br>
&gt; (2)   4<br>
&gt; (0)<br>
&gt; (0)   20 double<br>
&gt; (1)   30 double<br>
&gt; (2)   40 double<br>
&gt; (0)<br>
&gt; fatal:Coordinate subscript type mismatch. Subscript (0) can not be<br>
&gt; coerced to type of coordinate variable<br>
&gt; fatal:[&quot;Execute.c&quot;:8128]:Execute: Error occurred at or near line 26 in<br>
&gt; file mixed-mode-coord-subscripts.ncl<br>
&gt;<br>
&gt; I think it would be generally helpful if coordinate subscripting could<br>
&gt; be broadened to support all reasonable cases of mixed types.  At a<br>
&gt; minimum, please support the same combinations as for the relational<br>
&gt; operators.  I know that there are simple workarounds, but this<br>
&gt; improvement would benefit simpler user code and generalized file<br>
&gt; support.  Thanks for your consideration.<br>
&gt;<br>
&gt; --Dave A.<br>
&gt; NOAA/PSD/CIRES Climate Analysis Branch<br></blockquote></div></div></div></div>