Dear Ferret users,
I am trying to plot a PDF from precipitation data, but I
have some limitation problem in the number of samples that can
be used at once.
For instance, for the first year of data I use :
let
varprecip=precip[d=1,x=110e:145e,y=20n:45n,L=1:365]
LET hsort
= SAMPLEI(XSEQUENCE((varprecip)),
SORTI(XSEQUENCE((varprecip)) ))
LET ns = `hsort[i=@ngd]`
!->
DEFINE VARIABLE ns = 5794050
Now, if I try to extend this to many years :
let
varprecip=precip[d=1,x=110e:145e,y=20n:45n,L=1:10958]
LET hsort
= SAMPLEI(XSEQUENCE((varprecip)),
SORTI(XSEQUENCE((varprecip)) ))
LET ns = `hsort[i=@ngd]`
!-> DEFINE
VARIABLE ns = 7953862
Obviously the data are cut (because in the second case I
use 30 years, but the number of sample is barely the double of
the sample from 1 single year).
The consequence is only the light precipitation are
selected.
I suppose it is due to memory limitation, so I tried to
used
set
memory/size=200
But it is still not enough.
I guess the basic solution is to use a smaller region or
period, but maybe there is also another way in Ferret to solve
this problem?
(I want to avoid averaging, because the objective is to
identified the most extreme values)
Best,
Nicolas
--
----------------------------------------------------
Nicolas Freychet
Postdoc Reseach Fellow at RCEC
Academia Sinica
----------------------------------------------------