[ncl-talk] Creating histogram based on the number of gridpoints

[Lyndz]

Dear NCL-experts,

*[Details]*

I am using the aphrodite rainfall data set, which I uploaded in my dropbox
folder:
https://www.dropbox.com/s/lbkjlqgq9era5c7/APHRO_AMJ_1979-2015.nc?dl=0


I would like to count the number of grid points over a specified region
with rainfall in the following intervals(mm/day): 0-2,
2-4,4-6,6-8,8-10,10-12,12-14,14-16,16-18,18-20,20-22,22-24,24-26,26-28,28-30.
Then, create a histogram.

  f1  = addfile("APHRO_AMJ_1979-2015.nc","r")



*  latS = 12.5  latN = 20.0  lonL = 120.0  lonR = 122.5*
  prec = f1->precip(:,:,:)
  rain_region=prec(:,{latS:latN},{lonL:lonR})   ;select subregion

This is almost similar to this example
https://www.ncl.ucar.edu/Applications/Scripts/Katrina_circle_hist.ncl
But I am binning the number of grid points.

Any suggestions on how to do this in NCL?

I'll appreciate any guidance.

Sincerely,
Lyndz
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190909/152ac80f/attachment.html>