<div dir="auto">Wow..That is great information...Thanks and let me try....<div dir="auto">Jayakrishnan</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 4 Apr 2019, 20:06 Davide Sacchetti, <<a href="mailto:davide.sacchetti@arpal.gov.it">davide.sacchetti@arpal.gov.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">maskout(p1,p1-1000) displays p1>1000 only<br>
maskout(p1,2000-p1) displays<br>
p1<2000 only<br>
maskout(maskout(p1,p1>1000),2000-p1) displays p1 between<br>
1000 and 2000<br>
bye bye<br>
Davide<br>
<br>
On Thu, 2019-04-04 at 16:24 +0800, Jayakrishnan P.R wrote:<br>
> Dear Users,<br>
> I would like to classify the whole global grid into<br>
> four types using the quartiles obtained Q1 and Q3 using some<br>
> conditions in the two years 1950 and 2013. <br>
> For example for 1950, Q1=948, Q3=5768, for 1980 Q1=1008, Q3=17089<br>
> such that my type 1 would be<br>
> Type 1 grid: 1950 those values higher than Q3 and also 1980 higher<br>
> than Q3.<br>
> <br>
> I did the same using grads as follows. <br>
> <br>
> 'reinit'<br>
> 'sdfopen <a href="http://pp_1950.nc" rel="noreferrer noreferrer" target="_blank">pp_1950.nc</a>'<br>
> 'set gxout grfill'<br>
> 'define p1=ppp'<br>
> 'sdfopen <a href="http://pp_1980.nc" rel="noreferrer noreferrer" target="_blank">pp_1980.nc</a>'<br>
> 'define p2=ppp.2'<br>
> 'set grads off'<br>
> 'define mask1=const(maskout(p1,p1-5768*maskout(p2,p2-17089),1)'<br>
> 'd mask1'<br>
> 'set sdfwrite <a href="http://type1.nc" rel="noreferrer noreferrer" target="_blank">type1.nc</a>'<br>
> 'sdfwrite mask1'<br>
> 'draw title PPP in 1950 above 2973 and above 6371 in 1980: Type 1'<br>
> 'printim type1grid_above_2K.png white'<br>
> <br>
> Now my second condition for Type 2 in which I am confused about how<br>
> to do.<br>
> <br>
> type 2 : 1950 lower than Q1, 1980 higher than Q1; and 1950 between Q1<br>
> and Q3, 1980 higher than Q3 .<br>
> How to express the condition between Q1 and Q3 in maskout function<br>
> ?. <br>
> Kindly give some insights. <br>
> <br>
> Jayakrishnan<br>
> <br>
> <br>
> <br>
> -- <br>
> Sincerely<br>
> <br>
> **************************************************<br>
> Dr. Jayakrishnan P.R<br>
> Post Doctoral Research Fellow<br>
> Science Technology Building<br>
> School of System Science<br>
> Beijing Normal University<br>
> No. 96 XInjiekouwai Street<br>
> Beijing 100875<br>
> China<br>
> Tel: +86-188 11005715<br>
> <br>
> ****************************************************<br>
> <br>
> <br>
> <br>
> <br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org" target="_blank" rel="noreferrer">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank" rel="noreferrer">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</blockquote></div>