[gradsusr] probability distribution

[mehwish ramzan]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140902/4794eff5/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: prob_prec.1980..gif
Type: image/gif
Size: 44246 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20140902/4794eff5/attachment-0001.gif