[gradsusr] probability distribution

[Luis Blacutt]

Dear Mehwish,

Take a look at this web page, there you will find the histogram script
which can be tailored to the bins of your preference

http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html

Best regards
Luis


On Tue, Sep 2, 2014 at 8:17 AM, mehwish ramzan <mehwish.ramzan at gmail.com>
wrote:

> Dear GrADS users,
>
> I want to calculate the probability of daily extreme precipitationby
> classifying  multiple precipitation bins (e.g. 0-10, 10-20, 20-30, onwards)
> for 26 years of model output. And draw histogram for these bin against
> probability (%).
>
> I looked for grads scripts  but still didn't get any success.
>
> As an alternative, i tried to use the "const" function of GrADS to see the
> number of days which exceeded 10mm/day for example.
> ==================================================
>
> 'reinit'
> 'open r_pgb.ctl'
> 'set display color white'
> 'c'
> 'set lon 50 80'
> 'set lat 20 40'
> 'set time 1jan1980 25dec1980'
> 'define p1= const(const(maskout(prate*60*60*24, prate*60*60*24 -
> 10),1),0,-u)'
> 'set x 1'
> 'set y 1'
> 'd tloop(aave(p1, lon=50, lon=80, lat=20, lat=40))'
>
> ========================================================
>
> It worked fine and i get the output as well, but i am not sure of the
> units. (please see attached figure. Should i  multiply p1 by 100 to get
> percentage? ).
>
> I would really appreciate if someone guide me, regarding making
> precipitation bins in grad as it is my first priority.
>
>
> Thank you
>
> With Best Regards,
>
> Mehwish
>
>
>
>
>
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140902/c5cef422/attachment.html