<div dir="ltr"><div>Dear Mehwish,</div><div><br></div><div>Take a look at this web page, there you will find the histogram script which can be tailored to the bins of your preference</div><div><br></div><a href="http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html">http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html</a><br>
<div><br></div><div>Best regards<br></div><div>Luis</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Sep 2, 2014 at 8:17 AM, mehwish ramzan <span dir="ltr">&lt;<a href="mailto:mehwish.ramzan@gmail.com" target="_blank">mehwish.ramzan@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear GrADS users,<div><br></div><div>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 (%).   </div>

<div><br></div><div>I looked for grads scripts  but still didn&#39;t get any success. </div><div><br></div><div>As an alternative, i tried to use the &quot;const&quot; function of GrADS to see the number of days which exceeded 10mm/day for example.</div>

<div>==================================================</div><div><br></div><div><div>&#39;reinit&#39;</div><div>&#39;open r_pgb.ctl&#39;</div><div>&#39;set display color white&#39;</div><div>&#39;c&#39;</div><div>&#39;set lon 50 80&#39;</div>

<div>&#39;set lat 20 40&#39;</div><div>&#39;set time 1jan1980 25dec1980&#39;</div><div>&#39;define p1= const(const(maskout(prate*60*60*24, prate*60*60*24 - 10),1),0,-u)&#39;</div><div>&#39;set x 1&#39;</div><div>&#39;set y 1&#39;</div>

<div>&#39;d tloop(aave(p1, lon=50, lon=80, lat=20, lat=40))&#39;</div><div><br></div></div><div>========================================================</div><div><br></div><div>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? ). </div>

<div><br></div><div>I would really appreciate if someone guide me, regarding making precipitation bins in grad as it is my first priority. </div><div><br></div><div><br></div><div>Thank you</div><div><br></div><div>With Best Regards,</div>

<div><br></div><div>Mehwish</div><div><br></div><div><br></div><div> </div><div><br></div><div><div dir="ltr"><div style="font-family:arial;font-size:small"><font color="#500050"> </font></div><div><br></div><div><br></div>

</div>
</div></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>