<div dir="ltr"><div><br></div><div>Dear Luis,</div><div><br></div><div>Thank you for your reply. I have tried your suggestions but both of these tries gave me error message as follows;</div><div><br></div><div>Try 1:</div>
<div><br></div><div>a-> open r_pgb.ctl</div><div>Scanning description file: r_pgb.ctl</div><div>Data file r_pgb.%y4%m2%d2.avrg.grib is open as file 1</div><div>LON set to 7.756 128.125 </div><div>LAT set to -20.826 50.163 </div>
<div>LEV set to 1000 1000 </div><div>Time values set: 1980:1:1:0 1980:1:1:0 </div><div>ga-> set time 1jan1980 25dec1985</div><div>Time values set: 1980:1:1:0 1985:12:25:0 </div><div>ga-> hist t pratesfc*60*60*24 preciphist -2 2 0.25</div>
<div>[hist info] Temporary directory is set to /tmp/bGASL-mehwish</div><div>[hist info] Process ID is 12355</div><div>ga-> set lev -2 2</div><div>LEV set to -2 2 </div><div>ga-> d preciphist</div><div>Operation error: Invalid dimension environment</div>
<div> Too many varying dimensions </div><div><br></div><div>-----------------------------------------------------------------------</div><div>Try 2:</div><div><br></div><div><div>ga-> open r_pgb.ctl</div><div>Scanning description file: r_pgb.ctl</div>
<div>Data file r_pgb.%y4%m2%d2.avrg.grib is open as file 1</div><div>LON set to 7.756 128.125 </div><div>LAT set to -20.826 50.163 </div><div>LEV set to 1000 1000 </div><div>Time values set: 1980:1:1:0 1980:1:1:0 </div><div>
ga-> set lon 60 80</div><div>LON set to 60 80 </div><div>ga-> set lat 20 40</div><div>LAT set to 20 40 </div><div>ga-> hist xy pratesfc*60*60*24 preciphist -2 2 0.25</div><div>[hist info] Temporary directory is set to /tmp/bGASL-mehwish</div>
<div>[hist info] Process ID is 12355</div><div>ga-> set lon 7.756</div><div>LON set to 7.756 7.756 </div><div>ga-> set lat -20.826</div><div>LAT set to -20.826 -20.826 </div><div>ga-> set lev -2 2</div><div>LEV set to -2 2 </div>
<div>ga-> set xyrev on</div><div>ga-> d preciphist</div><div>Data Request Error: Invalid grid coordinates</div><div> World coordinates convert to non-integer grid coordinates</div><div> Variable = preciphist Dimension = 0 </div>
<div> Error ocurred at column 1</div><div>DISPLAY error: Invalid expression </div><div> Expression = preciphist</div></div><div><br></div><div><br></div><div>=========================================</div><div><br></div>
<div>Can you guide me which step i am doing wrong? or any other way to calculate daily extreme precipitation events for 26 years data.</div><div><br></div><div>I will be really grateful. </div><div><br></div><div><br></div>
<div><br></div><div>With Best Regards,</div><div><br></div><div>Mehwish</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Sep 2, 2014 at 8:40 PM, Luis Blacutt <span dir="ltr"><<a href="mailto:luis.blacutt@gmail.com" target="_blank">luis.blacutt@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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" target="_blank">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"><div><div class="h5">On Tue, Sep 2, 2014 at 8:17 AM, mehwish ramzan <span dir="ltr"><<a href="mailto:mehwish.ramzan@gmail.com" target="_blank">mehwish.ramzan@gmail.com</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><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>
<br></div></div>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
</blockquote></div><br><br clear="all"><div><br></div><br><div dir="ltr"><div style="font-family:arial;font-size:small"><font color="#666666"></font></div><div style="font-family:arial;font-size:small"><font color="#500050"> </font></div>
<div><br></div><div><br></div></div>
</div></div>