<div dir="ltr"><div>You can count yourself:</div><div><br></div><div>  nRain0    = <b>num</b>(rain_region.eq.0)<br>  nRain     = num(rain_region.gt.0)<br>  nRainMsg  = num(ismissing(rain_region))<br><br>  print("  nRain0="+nRain0)<br>  print("  nRain ="+nRain )<br>  print("nRainMsg="+nRainMsg)</div><div><br></div><div>  opt = True<br>  opt@PrintStat = True<br>  stat_rain_region = <b>stat_dispersion</b>(rain_region, opt )<br>print("===")<br><br>  n_1015 = num(rain_region.ge.10 .and. rain_region.lt.15)<br>  print("n_1015="+n_1015)<br>print("===")<br><br>  test   = <b>ndtooned</b>(rain_region)<br>  iGood  = <b>ind</b>(test.gt.0)         ; only grid points > 0.0<br>  test  <b>:=</b> test(iGood)<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Sep 14, 2019 at 2:51 AM Lyndz via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.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="ltr"><div dir="ltr" class="gmail-m_3096971731177391944gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear NCL experts,</div><div><font size="4"><br></font></div><div><b><font size="4">[Details]</font></b></div><div>I  want to plot a histogram with the number of grid points in the y-axis and rainfall classification in the x-axis.</div><div>The attached figure is an example of what I wanted (Desired_output.gif).</div><div><br></div><div>The data that I am using can be downloaded from this link:</div><div><a href="https://www.dropbox.com/s/lbkjlqgq9era5c7/APHRO_AMJ_1979-2015.nc?dl=0" target="_blank">https://www.dropbox.com/s/lbkjlqgq9era5c7/APHRO_AMJ_1979-2015.nc?dl=0</a><br></div><div><br></div><div><b><font size="4">[What I have so far]</font></b></div><div>The following is the script that I am using:</div><div><br></div><div>;***********************************************************<br>;This script creates a histogram of the number of gridpoints<br>;with rainfall above 0 mm/day<br>;Created by Lyndz<br>;<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@gmail.com</a><br>;*******************************************************</div><div>begin</div><div>fo = addfile("APHRO_AMJ_1979-2015.nc","r")<br>rain=fo->precip(:,:,:)<br>rain_region=rain(:,{12.5:22},{120:122.5})<br><b>test=ndtooned(rain_region)</b><br>;*************************************<br>;Plots histogram<br>;*************************************<br>  wks = gsn_open_wks("png","output")             ; send graphics to PNG file<br>  gsn_define_colormap(wks,"temp1")              ; choose colormap<br>  res                                 = True<br>  res@gsnHistogramSelectNiceIntervals = False   ; intervals now float<br>  res@gsnHistogramComputePercentages  = True    ; change left axis to %<br>  res@gsnHistogramClassIntervals = (/0,1,5,10,15,20,25,30,25,40,45/)<br>  res@gsnHistogramComputePercentagesNoMissing = True<br>  res@tmXBPrecision = 2         ; set the precision to 2 significant digits<br>  res@tmYLFormat = "#^se"       ; don't use "e" symbol, always show +/-<br>  res@tiXAxisString   = "Rainfall[mm/day]"<br>  res@tiYAxisString   = "Number of Grid Points"<br>  plot=gsn_histogram(wks,test,res)                 ; create histogram with 10 bins<br>end<br></div><div><br></div><div><b><font size="4">[Issues/Problems]</font></b></div><div>[1] The script runs fine. The output is attached (output.png). I am not sure if  I plotted this correctly. After using ndtooned, I just plotted the histogram. Is there a way to check if the y-axis is really the number of gridpoints or the x axis is the rainfall class?</div><div>Or is there a way to invoke that the y-axis is the number of gridpoints when plotting the histogram?</div><div><br></div><div>I'll appreciate any help with this.</div><div><br></div><div>Sincerely,</div><div>Lyndz</div><div><br></div></div></div></div></div></div></div></div></div></div></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></blockquote></div>