<div dir="ltr"><div><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_pqsort_n.shtml"><b>dim_pqsort_n</b></a></div><div><br></div><div><b> </b><br><div>"I have an array of HI(time), T(time), and RH(time).</div>
<div> I want to qsort(HI), then reorder T and RH to match qsort(HI)."</div><div><br></div><div>====</div><div><br></div><div>   ip  = <strong>dim_pqsort_n</strong>(HI,2,0)</div><div>   T  = T(ip)</div><div>   RH = RH(ip)</div><div><br></div><div>   print(<strong> sprintf</strong>("%9.2f", HI)+"  "+<strong>sprintf</strong>("%9.2f", T)+"  "+<strong>sprintf</strong>("%9.2f",RH) )<br></div><div><pre><br></pre></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Feb 26, 2019 at 7:12 AM Buzan, Jonathan <<a href="mailto:jbuzan@purdue.edu" target="_blank">jbuzan@purdue.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">



<div>
<div>Hi Dennis, </div>
<div><br>
</div>
<div>I will work on coding up the suggested algorithm…</div>
<div>However, it might be easier to rephrase the question:</div>
<div><br>
</div>
<div>I have an array of HI(time), T(time), and RH(time).</div>
<div>I want to qsort(HI), then reorder T and RH to match qsort(HI).</div>
<div><br>
</div>
<div>Is there a way to do that?  </div>
<div>Because I can simply use index numbers, and get around this issue.</div>
<div><br>
</div>
<div>-Jonathan</div>
<div><br>
</div>
<div><br>
</div>
<br>
<div><br>
<blockquote type="cite">
<div>On Feb 25, 2019, at 5:31 PM, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:</div>
<br class="gmail-m_2258278169725666802gmail-m_4566856208969222705Apple-interchange-newline">
<div>
<div dir="ltr">
<div>Jonathan:</div>
<div><br>
</div>
<div>For a sort variation, you might try <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_pqsort_n.shtml" target="_blank">
<b>dim_pqsort_n</b></a>  which is from the fortran <a href="https://www.netlib.org/slatec/guide" target="_blank">
<b>SLATEC </b></a>library.</div>
<div><br>
</div>
<div>NCL's <b>dim_pqsort_n </b>should be able to handle any dimension size.</div>
<div><br>
</div>
<div>Under-the-hood:   To my knowledge, fortran's default integer type is 32-bits [ 4 bytes ] =>
<span class="gmail-m_2258278169725666802gmail-m_4566856208969222705gmail-ILfuVd">2,147,483,647 <br>
</span></div>
<div><span class="gmail-m_2258278169725666802gmail-m_4566856208969222705gmail-ILfuVd">NCL's C-interface checks the size of the dimension being sorted. I it  exceeds the 32-bit size, NCL uses a special version that allows [64-bit; 8-byte; INTEGER*8] indexing.<br>
</span></div>
<div><br>
</div>
<div>Cheers</div>
<div>D<br>
</div>
<div><br>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Mon, Feb 25, 2019 at 8:17 AM Buzan, Jonathan <<a href="mailto:jbuzan@purdue.edu" target="_blank">jbuzan@purdue.edu</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="auto">
<div></div>
<div>Hi Rick,</div>
<div><br>
</div>
<div>I am running into an issue where I have some repeat values.  Is there a way to add a precision level adjustment during the qsort algorithm? </div>
<div><br>
</div>
<div>-Jonathan</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
On Feb 25, 2019, at 16:06, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
<br>
</div>
<blockquote type="cite">
<div>
<div dir="ltr">
<div>Hi Jonathan -- yes, qsort() takes a double, or more generally a "numeric" type.</div>
<div><br>
</div>
<div>RB<br>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Mon, Feb 25, 2019 at 7:20 AM Buzan, Jonathan <<a href="mailto:jbuzan@purdue.edu" target="_blank">jbuzan@purdue.edu</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Hello NCL-Talk,<br>
<br>
Does qsort operate at double precision?<br>
I am trying to resolve my stat_dispersion issue by isolating the problem...<br>
<br>
-Jonathan<br>
<br>
<br>
_______________________________________________<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/mailman/listinfo/ncl-talk</a><br>
</blockquote>
</div>
</div>
</blockquote>
</div>
_______________________________________________<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/mailman/listinfo/ncl-talk</a><br>
</blockquote>
</div>
<span id="gmail-m_2258278169725666802gmail-m_4566856208969222705cid:f_jskjg82e0"><dpsort.f></span></div>
</blockquote>
</div>
<br>
</div>

</blockquote></div>