[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