[Thread Prev][Thread Next][Index]

[ferret_users] PDF - limitation in the number of samples



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
----------------------------------------------------


[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement