<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't get any success. </div><div><br></div><div>As an alternative, i tried to use the "const" function of GrADS to see the number of days which exceeded 10mm/day for example.</div>
<div>==================================================</div><div><br></div><div><div>'reinit'</div><div>'open r_pgb.ctl'</div><div>'set display color white'</div><div>'c'</div><div>'set lon 50 80'</div>
<div>'set lat 20 40'</div><div>'set time 1jan1980 25dec1980'</div><div>'define p1= const(const(maskout(prate*60*60*24, prate*60*60*24 - 10),1),0,-u)'</div><div>'set x 1'</div><div>'set y 1'</div>
<div>'d tloop(aave(p1, lon=50, lon=80, lat=20, lat=40))'</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>