[ncl-talk] Question about ndtooned and histogram plot

Lyndz olagueralyndonmark429 at gmail.com
Mon Sep 16 10:30:08 MDT 2019


Dear Sir Dennis,

Thank you so much!!
This is what I'm looking for.

Sincerely,

Lyndz



On Mon, Sep 16, 2019 at 11:15 PM Dennis Shea <shea at ucar.edu> wrote:

> You can count yourself:
>
>   nRain0    = *num*(rain_region.eq.0)
>   nRain     = num(rain_region.gt.0)
>   nRainMsg  = num(ismissing(rain_region))
>
>   print("  nRain0="+nRain0)
>   print("  nRain ="+nRain )
>   print("nRainMsg="+nRainMsg)
>
>   opt = True
>   opt at PrintStat = True
>   stat_rain_region = *stat_dispersion*(rain_region, opt )
> print("===")
>
>   n_1015 = num(rain_region.ge.10 .and. rain_region.lt.15)
>   print("n_1015="+n_1015)
> print("===")
>
>   test   = *ndtooned*(rain_region)
>   iGood  = *ind*(test.gt.0)         ; only grid points > 0.0
>   test  *:=* test(iGood)
>
> On Sat, Sep 14, 2019 at 2:51 AM Lyndz via ncl-talk <ncl-talk at ucar.edu>
> wrote:
>
>> Dear NCL experts,
>>
>> *[Details]*
>> I  want to plot a histogram with the number of grid points in the y-axis
>> and rainfall classification in the x-axis.
>> The attached figure is an example of what I wanted (Desired_output.gif).
>>
>> The data that I am using can be downloaded from this link:
>> https://www.dropbox.com/s/lbkjlqgq9era5c7/APHRO_AMJ_1979-2015.nc?dl=0
>>
>> *[What I have so far]*
>> The following is the script that I am using:
>>
>> ;***********************************************************
>> ;This script creates a histogram of the number of gridpoints
>> ;with rainfall above 0 mm/day
>> ;Created by Lyndz
>> ;olagueralyndonmark429 at gmail.com
>> ;*******************************************************
>> begin
>> fo = addfile("APHRO_AMJ_1979-2015.nc","r")
>> rain=fo->precip(:,:,:)
>> rain_region=rain(:,{12.5:22},{120:122.5})
>> *test=ndtooned(rain_region)*
>> ;*************************************
>> ;Plots histogram
>> ;*************************************
>>   wks = gsn_open_wks("png","output")             ; send graphics to PNG
>> file
>>   gsn_define_colormap(wks,"temp1")              ; choose colormap
>>   res                                 = True
>>   res at gsnHistogramSelectNiceIntervals = False   ; intervals now float
>>   res at gsnHistogramComputePercentages  = True    ; change left axis to %
>>   res at gsnHistogramClassIntervals = (/0,1,5,10,15,20,25,30,25,40,45/)
>>   res at gsnHistogramComputePercentagesNoMissing = True
>>   res at tmXBPrecision = 2         ; set the precision to 2 significant
>> digits
>>   res at tmYLFormat = "#^se"       ; don't use "e" symbol, always show +/-
>>   res at tiXAxisString   = "Rainfall[mm/day]"
>>   res at tiYAxisString   = "Number of Grid Points"
>>   plot=gsn_histogram(wks,test,res)                 ; create histogram
>> with 10 bins
>> end
>>
>> *[Issues/Problems]*
>> [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?
>> Or is there a way to invoke that the y-axis is the number of gridpoints
>> when plotting the histogram?
>>
>> I'll appreciate any help with this.
>>
>> Sincerely,
>> Lyndz
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190917/9053df02/attachment.html>


More information about the ncl-talk mailing list