<div dir="ltr">Thank you very much Hosmay, I really appreciate your help. That was just what I needed, there's still a lot to learn using grads. <div><br></div><div>Miquel</div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-06-19 16:05 GMT+02:00 Hosmay Lopez <span dir="ltr"><<a href="mailto:hlopez@rsmas.miami.edu" target="_blank">hlopez@rsmas.miami.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">Hi Miguel,
<div><br>
</div>
<div>This is straightforward in GrADS if you use a combination of functions. It goes like this:</div>
<div><br>
</div>
<div>If you have a variable for precipitation named <b>precip </b>and let say, 100 ensemble members. You do:</div>
<div><br>
</div>
<div>set e 1 100</div>
<div>define prec5mm = const( const( maskout( precip, precip - 5), 1), 0.0, -u)</div>
<div><br>
</div>
<div>The <i>maskout</i> function sets the values to undefined where precip < 5</div>
<div>The inner <i>const</i> function sets the remaining defined values to ones.</div>
<div>And the outer <i>const</i> function sets the undefined values to zeros.</div>
<div>So the new variable prec5mm is made of 1 wherever precip > 5 and 0 otherwise.</div>
<div><br>
</div>
<div>Now you average for all your ensemble members and multiply the results by 100 to get % of ensemble members with precipitation above 5mm .</div>
<div><br>
</div>
<div>set e 1</div>
<div>define percentage5mm = 100 * ave( prec5mm, e=1, e=100)</div>
<div><br>
<div>Hope this helps.</div>
<div><br>
<div style="font-family:Tahoma;font-size:13px">
<div style="font-family:Tahoma;font-size:13px">
<div style="font-family:Tahoma;font-size:13px">
<div style="font-family:Tahoma;font-size:13px">Hosmay Lopez<br>
<span style="font-family:verdana,arial,sans-serif;background-color:rgb(255,255,255)">UM/CIMAS/RSMAS NOAA/AOML/PHOD</span><br style="font-family:verdana,arial,sans-serif;background-color:rgb(255,255,255)">
<span style="font-family:verdana,arial,sans-serif;background-color:rgb(255,255,255)">4301 Rickenbacker Causeway,</span><br style="font-family:verdana,arial,sans-serif;background-color:rgb(255,255,255)">
<span style="font-family:verdana,arial,sans-serif;background-color:rgb(255,255,255)">Miami, FL 33149, USA</span><br style="font-family:verdana,arial,sans-serif;background-color:rgb(255,255,255)">
<span style="font-family:verdana,arial,sans-serif;background-color:rgb(255,255,255)">Phone: <a href="tel:%28305%29%20361-4329" value="+13053614329" target="_blank">(305) 361-4329</a></span><br>
<br>
</div>
</div>
</div>
</div>
</div>
<div style="font-family:Times New Roman;color:#000000;font-size:16px">
<hr>
<div style="direction:ltr"><font face="Tahoma" size="2" color="#000000"><b>From:</b> <a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a> [<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>] on behalf of Miquel Bernis [<a href="mailto:bernisfeina@gmail.com" target="_blank">bernisfeina@gmail.com</a>]<br>
<b>Sent:</b> Friday, June 19, 2015 3:26 AM<br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> [gradsusr] Probability precipitation map using ensemble forecast<br>
</font><br>
</div><div><div class="h5">
<div></div>
<div>
<div dir="ltr">Hi everyone,
<div><br>
</div>
<div>Do you think is there any way to plot a forecast probability precipitation map, using an ensemble forecast? I mean, calculating which forecasts of the ensemble give, for example, an amount of rain > 5 mm and which not, for an specific moment of time, and
displaying a contour map of the % you can take from that. Watching to the grads functions it seems not possible, but maybe some of you, with more experience on grads than me, can imagine a way to do that. </div>
<div><br>
</div>
<div>Thanks a lot!</div>
<div>Miquel</div>
</div>
</div>
</div></div></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" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>