<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Bill,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default">You can find some a link to the C code here, which is the same C code that the NCL one was wrapped on top of:</div><div class="gmail_default"><br></div><div class="gmail_default" style="font-size:small"><a href="https://cran.r-project.org/web/packages/SPEI/index.html">https://cran.r-project.org/web/packages/SPEI/index.html</a></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">We unfortunately had to make the decision to pull the speidx code from NCL, because the C code was deprecated by the original author(s). There's a known issue in the code that occasionally causes Nans (not-a-number) to show up because of a divide-by-0 in some cases. They've fixed the problem, but only through a pure R interface they've created.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">We will keep a "hidden" version of the function in NCL, when 6.4.0 is released, but we won't advertise or support it. It will mainly be there in the case that it might work for some people.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I've emailed the author of the code to see if we can get an updated C version, but I haven't heard anything back.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Sep 2, 2016 at 8:24 AM, Capehart, William J <span dir="ltr"><<a href="mailto:William.Capehart@sdsmt.edu" target="_blank">William.Capehart@sdsmt.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="white" lang="EN-US" link="#0563C1" vlink="#954F72">
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">Hi Dennis, Mary et al.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Is it possible to provide the source code for the speidx() function for the SPEI that is due out with ncl 6.4? I saw that the nws heat index code has been provided on the NCL web manual pages? Can this also
be done for speidx()?<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Cheers and Thanks Much<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:Consolas;color:black">==============================<wbr>==============================<wbr>====<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:Consolas;color:black">Bill Capehart <<a href="mailto:William.Capehart@sdsmt.edu" target="_blank">William.Capehart@sdsmt.edu</a>> Associate Professor<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:Consolas;color:black">Atmospheric & Environmental Sciences Program Director<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:Consolas;color:black">Civil and Environmental Engineering Hydrometeorology<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:Consolas;color:black">201 Mineral Industries Building <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:Consolas;color:black">South Dakota School of Mines and Technology Skype: wcapehart<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:Consolas;color:black">501 East Saint Joseph Street Ph: <a href="tel:%2B1-605-394-1994" value="+16053941994" target="_blank">+1-605-394-1994</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:Consolas;color:black">Rapid City, SD 57701-3995 Mobile: <a href="tel:%2B1-605-484-5692" value="+16054845692" target="_blank">+1-605-484-5692</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:Consolas;color:black">=================== <a href="http://capehart.sdsmt.edu" target="_blank">http://capehart.sdsmt.edu</a> ==================</span><span style="font-family:"Times New Roman""><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
</div>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>