[gradsusr] gradsusr Digest, Vol 55, Issue 2

[semou ndao]

Hello.
I would use ltrend.gs script but I dont know how do it.
best
evra




Le Mardi 2 septembre 2014 11h45, semou ndao <semoundaola at yahoo.fr> a écrit :
 


Hello.
I would use ltrend.gs script but I dont now how do it.
best
evra



Le Mardi 2 septembre 2014 11h28, "gradsusr-request at gradsusr.org" <gradsusr-request at gradsusr.org> a écrit :
 


Send gradsusr mailing list submissions to
    gradsusr at gradsusr.org

To subscribe or unsubscribe via the World Wide Web, visit
    http://gradsusr.org/mailman/listinfo/gradsusr
or, via email, send a message with subject or body 'help' to
    gradsusr-request at gradsusr.org

You can reach the person managing the list at
    gradsusr-owner at gradsusr.org

When replying, please edit your Subject line so it is more specific
than "Re: Contents of gradsusr digest..."


Today's Topics:

   1. probability distribution (mehwish ramzan)


----------------------------------------------------------------------

Message: 1
Date: Tue, 2 Sep 2014 20:17:32 +0900
From: mehwish ramzan <mehwish.ramzan at gmail.com>
Subject: [gradsusr] probability distribution
To: GrADS Users Forum <gradsusr at gradsusr.org>
Message-ID:
    <CA+rMR9wJsPJ_dD-XOT1O=ewHmfNYMWsgPnOPC1w1vBY4ytbk=w at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

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.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.gif 

------------------------------

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr


End of gradsusr Digest, Vol 55, Issue 2
***************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140902/4d6ab267/attachment-0001.html