<div dir="ltr"><div dir="ltr"><div>Hi Jonathan,</div><div><br></div><div>Attached is an NCL script and a fortran-77 code from Auguat 2014. For some reason, these never made it into NCL. Likely no completely tested.</div><div><br></div><div> </div><div>Using the NCL function would likely be quite slow for large arrays. Give it a try.</div><div><br></div><div>%> ncl tst_pdfxyx.ncl</div><div><br></div><div>Ypu could invoke the f77 subroutine if needed. I can't find the test for this. I'll look later today.<br></div><div><br></div><div>%> WRAPIT xyz2pdf77.f</div><div><br></div><div>Good luck</div><div><br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 4, 2019 at 2:34 AM Jonathan Buzan <<a href="mailto:jrbuzan@gmail.com">jrbuzan@gmail.com</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 style="overflow-wrap: break-word;"><div style="overflow-wrap: break-word;">Hello NCL-Talk and Dennis Shea,<div><br></div><div>Thanks again for the previous help, as I have had much success! </div><div>I am moving on to more advanced subjects. Instead of PDF(x,y), I now need to calculate PDF(x,y,z).</div><div><br></div><div>In the contributed.ncl file, I believe I can modify the code below by adding a 3rd variable into the binning if-then-do statement. Would this be along the right track? </div><div>I would need to create a new routine to accept the third variable in the pdfxy. Probably call it pdfxyz… </div><div><br></div><div>-Jonathan</div><div><br></div><div><br></div><div>; JRB BEGIN<br>;print("Jon Is Here")<br>; pdf2 = pdfxy_bin(x,y, binxBound,binyBound, popt)<br> ; all NCL [original code]<br> x1d = ndtooned(x) ; requires more memory<br> y1d = ndtooned(y)<br><br> do nby=0,nbinsy-1<br> iy = ind(y1d.ge.binyBound(nby) .and. y1d.lt.binyBound(nby+1))<br> if (.not.ismissing(iy(0))) then<br> do nbx=0,nbinsx-1<br> pdf2(nby,nbx) = num(x1d(iy).ge.binxBound(nbx) .and. \<br> x1d(iy).lt.binxBound(nbx+1) )<br> end do<br> end if<br> delete(iy) ; size may change<br> end do<br><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br><div><br><blockquote type="cite"><div>On Feb 20, 2019, at 3:30 PM, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:</div><br class="gmail-m_-1618550062252860291Apple-interchange-newline"><div><div dir="ltr"><div>Hi Jonathan, <br></div><div><br></div><div>No. The pdfxy/pdfxy_bin are 'cartesian oriented' (if there is such a thing!).</div><div>Hence, x/y rathern tha lat/lon.</div><div><br></div><div>D<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Feb 20, 2019 at 6:40 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>
Do the pdfxy/pdfxy_bin functions account for lat-lon grid cell size changes? i.e. if there is no accounting for weights, the number of events at >30° are 2x to the number events at <30°, even though there are equal area.<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></div></div></div></blockquote></div>