Re: [ferret_users] percentile

[Peter Szabo <szabpet83@xxxxxxxxx>]

Hello Izidine,

i will not tell you the exact commands, but hints.

first, decide if TEMP[d=1] < TMXDJFBBOTTOM10[d=2] or not: let YES=if X eq Y then 1 (EQual; but Lesser Than, Lesser Equal, Greater Equal,  Greater Than also exist)

count the occasions when TEMP[d=1] < TMXDJFBBOTTOM10[d=2]: count with @SUM function through L dimension.

get percentage: have now the rate of good_occasions/all_occasions*100

:)

Hope it helps,

Peter

On Mon, May 10, 2010 at 12:18 PM, Izidine Pinto <izidinep@xxxxxxxxx> wrote:

Dear ferret users

I have two data files

1> /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  (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
:-)