[gradsusr] Precipitation range

[Davide Sacchetti]

use nested maskout:

maskout(maskout(1, prate*86400-1),10-prate*86400) (I didn't checked but
it looks reasonable ...)

bye bye
Davide

On Wed, 2014-10-29 at 18:25 +0900, mehwish ramzan wrote:
> Dear GrADS users,
> 
> 
> I used the  following script to calculate the number of days which
> exceed  precipitation by setting one limit such as below;
> 
> 
> 
>  
> open r_pgb.ctl
> 
> 
> set lon 78
> 
> 
> set lat 10
> 
> 
> define p1= sum(maskout(1, prate*86400 - 1.) , t=1, t=9491)
> d p1
> 
> 
> 
> 
> It gave me the number of days which exceeds  1mm/day precipitation.But
> now i want to set the range such as greater then 1 and less then 10.
> 
> 
> Can someone please guide me how to set the range? 
> 
> 
> Am using grads version 1.9.0.
> 
> 
> 
> 
> 
> 
> Thank you
> 
> 
> With Best Regards,
> 
> 
> Mehwish
> 
> 
> 
> 
> 
> 
>  
> 
> 
> 
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

-- 
Davide Sacchetti
Centro Funzionale Meteo Idrologico di Protezione Civile della Regione Liguria
ARPAL Unità Tecnica Complessa di livello Regionale
V.le Brigare Partigiane 2 16121 Genova (I)
tel: +39 010 6437535                    fax: +39 010 6437520
mail: davide.sacchetti at arpal.gov.it     web: www.meteoliguria.it