<div dir="ltr"><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><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><br></div><div><b>[Details]</b></div><div><b><br></b></div><div>I am using the aphrodite rainfall data set, which I uploaded in my dropbox folder:</div><div><a href="https://www.dropbox.com/s/lbkjlqgq9era5c7/APHRO_AMJ_1979-2015.nc?dl=0">https://www.dropbox.com/s/lbkjlqgq9era5c7/APHRO_AMJ_1979-2015.nc?dl=0</a><br></div><div><br></div><div><br></div><div>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.</div><div><br></div><div>  f1  = addfile("APHRO_AMJ_1979-2015.nc","r")<br><b>  latS = 12.5<br>  latN = 20.0<br>  lonL = 120.0<br>  lonR = 122.5</b><br>  prec = f1->precip(:,:,:)<br>  rain_region=prec(:,{latS:latN},{lonL:lonR})   ;select subregion<br></div><div><br></div><div>This is almost similar to this example <a href="https://www.ncl.ucar.edu/Applications/Scripts/Katrina_circle_hist.ncl">https://www.ncl.ucar.edu/Applications/Scripts/Katrina_circle_hist.ncl</a></div><div>But I am binning the number of grid points.</div><div><br></div><div>Any suggestions on how to do this in NCL?</div><div><br></div><div>I'll appreciate any guidance.</div><div><br></div><div>Sincerely,</div><div>Lyndz</div></div></div></div></div></div></div></div></div></div>