[Thread Prev][Thread Next][Index]

Re: [ferret_users] percentile



Hi Izidine

I suggest the following:

let temp_below = if TEMP[d=1] lt TMXDJFBBOTTOM10[d=2] then 1 else 0
let temp_all   = if TEMP[d=1] then 1
let temp_below_percent = temp_below[t=@din]/temp_all[t=@din]*100.0

shade temp_below_percent

This should give you a map that shows the percentage of the time where
TEMP[d=1] < TMXDJFBBOTTOM10[d=2] at each location.

Marco


Izidine Pinto wrote:
> Dear ferret users
> 
> I have two data files
> 
> 1> /Tmax.DJF.daily.1981_1999.echama2.nc
> <http://Tmax.DJF.daily.1981_1999.echama2.nc>
>  name     title                                                    
> I         J         K         L
>  TEMP     TEMPERATURE AT 1.5M              1:127     1:110     1:1      
> 1:1680
>  
>     2> /tmxdjf10th.b.ech.nc <http://tmxdjf10th.b.ech.nc>  (default)
>  name     title                                            I        
> J         K         L
>  TMXDJFBBOTTOM10
>           DJF Max temp 10th percentile 19  1:127     1:110     1:1       ...
> 
> 
> I was wondering if someone could help me to  calculate the percentage of
> time where
> TEMP[d=1] < TMXDJFBBOTTOM10[d=2]
> 
> 
> Any help is strongly appreciated
> -- 
> Izidine Pinto
> :-)


-- 
***************************************
Marco Steinacher

Climate and Environmental Physics
Physics Institute, University of Bern
Sidlerstrasse 5, CH-3012 Bern

Phone ++41 (0)31 631 34 02
Fax   ++41 (0)31 631 87 42
steinacher@xxxxxxxxxxxxxxxx
http://www.climate.unibe.ch/~steinach
OpenPGP Key ID: 0x62937F7F
***************************************

Attachment: signature.asc
Description: OpenPGP digital signature


[Thread Prev][Thread Next][Index]

Contact Us
Dept of Commerce / NOAA / OAR / PMEL / TMAP

Privacy Policy | Disclaimer | Accessibility Statement