<html>
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
  </head>
  <body>
    <p>Could someone please advice me on the following?<br>
    </p>
    <div class="moz-cite-prefix">On 05.06.20 09:28, Sri nandini via
      ncl-talk wrote:<br>
    </div>
    <blockquote type="cite"
      cite="mid:396f3f2a-b358-0706-c4da-bf3c44f0c52e@uni-hamburg.de">Hello
      ncl users,
      <br>
      <br>
      I have a question on plotting my PDF, which changes (as it should)
      with specific bin_nice and without.
      <br>
      <br>
      Below is an example of my script, and the attached 2 plots, once
      with using the bin_nice and one with using the bin_min and max.
      <br>
      <br>
      I don't understand why bin space and values change, can some one
      help me understand this?
      <br>
      <br>
      Thanx Sri
      <br>
      <br>
      ~~~~~ snippet of the main PDF code~~~~
      <br>
      <br>
      <br>
      print("calculate PDFs...")
      <br>
      ;; calculate empirical PDF
      <br>
      <br>
        opt          = True
      <br>
        ;opt@bin_min = -50.  (i have also tried it with 100,75..)
      <br>
        ;opt@bin_max =  50.
      <br>
        opt@bin_nice=True
      <br>
        pdf1 = pdfx(hist_anom1, 100, opt)
      <br>
        printVarSummary(pdf1)
      <br>
        pdf2 = pdfx(rcp45_anom1, 100, opt)
      <br>
        printVarSummary(pdf2)
      <br>
        pdf3 = pdfx(rcp85_anom1,100, opt)
      <br>
        printVarSummary(pdf3)
      <br>
      <br>
;================================================================================
      <br>
      ; prepare data for plotting
      <br>
;================================================================================
      <br>
      <br>
        nVar    = 3
      <br>
        nBin    = pdf1@nbins          ; retrieve the number of bins
      <br>
      <br>
        xx      = new ( (/nVar, nBin/), typeof(pdf1))
      <br>
      <br>
        xx(0,:) = pdf1@bin_center     ; assign appropriate "x" axis
      values
      <br>
        xx(1,:) = pdf2@bin_center
      <br>
        xx(2,:) = pdf3@bin_center
      <br>
      <br>
        yy      = new ( (/nVar, nBin/), typeof(pdf1))
      <br>
        yy(0,:) = (/ pdf1 /) * 0.01 ;; convert % to absolut
      <br>
        yy(1,:) = (/ pdf2 /) * 0.01
      <br>
        yy(2,:) = (/ pdf3 /) * 0.01
      <br>
      <br>
      <br>
      ~~~Than plotting~~~~~
      <br>
      <br>
      <br>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <pre class="moz-quote-pre" wrap="">_______________________________________________
ncl-talk mailing list
<a class="moz-txt-link-abbreviated" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>
List instructions, subscriber options, unsubscribe:
<a class="moz-txt-link-freetext" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></pre>
    </blockquote>
  </body>
</html>