[gradsusr] Mask out function help

[Davide Sacchetti]

maskout(p1,p1-1000) displays p1>1000 only
maskout(p1,2000-p1) displays
p1<2000 only
maskout(maskout(p1,p1>1000),2000-p1) displays p1 between
1000 and 2000
bye bye
Davide

On Thu, 2019-04-04 at 16:24 +0800, Jayakrishnan P.R wrote:
> Dear Users,
>                  I would like to classify the whole global grid into
> four types using the quartiles obtained Q1 and Q3 using some
> conditions in the two years 1950 and 2013. 
> For example for 1950, Q1=948, Q3=5768, for 1980 Q1=1008, Q3=17089
> such that my type 1 would be
> Type 1 grid: 1950 those values higher than Q3 and also 1980 higher
> than Q3.
> 
> I did the same using grads as follows. 
> 
> 'reinit'
> 'sdfopen pp_1950.nc'
> 'set gxout grfill'
> 'define p1=ppp'
> 'sdfopen pp_1980.nc'
> 'define p2=ppp.2'
> 'set grads off'
> 'define mask1=const(maskout(p1,p1-5768*maskout(p2,p2-17089),1)'
> 'd mask1'
> 'set sdfwrite type1.nc'
> 'sdfwrite mask1'
> 'draw title PPP in 1950 above 2973 and above 6371 in 1980: Type 1'
> 'printim type1grid_above_2K.png white'
> 
> Now my second condition for Type 2 in which I am confused about how
> to do.
> 
> type 2 : 1950 lower than Q1, 1980 higher than Q1; and 1950 between Q1
> and Q3, 1980 higher than Q3 .
> How to express the condition between Q1 and Q3 in maskout function
> ?. 
> Kindly give some insights. 
> 
> Jayakrishnan
> 
> 
> 
> -- 
> Sincerely
> 
> **************************************************
> Dr. Jayakrishnan P.R
> Post Doctoral Research Fellow
> Science Technology Building
> School of System Science
> Beijing Normal University
> No. 96 XInjiekouwai Street
> Beijing 100875
> China
> Tel: +86-188 11005715
> 
> ****************************************************
> 
> 
> 
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr