Stat_Dispersion

Hello Mary Haley,

I am working with the stat_dispersion function, and I have noticed that there are 8 available slots for extra functions. It would be useful to see 95%, 99%, and 99.9% in the 'contributed.ncl'

The following is an example test script, and the lines I added to 'contributed.ncl'

-Jonathan R. Buzan
PhD Student; Earth, Atmospheric, and Planetary Sciences
Purdue University
550 Stadium Mall Dr.
West Lafayette, IN 47907
jbuzan@purdue.edu

;;;;;;;;;;;;;;;;;;;;
;;;Modified contributed.ncl lines under stat_dispersion;;;
        ;JBuzan BEGIN
       statx(22)= work((nwork*19)/20-1) ; upper 95%
       statx(23)= work((nwork*99)/100-1) ; upper 99.%
       statx(24)= work((nwork*999)/1000-1) ; upper 99.9%
       statx(25) = work( nwork/20-1 ) ; lower 95%
       statx(26) = work( nwork/100-1 ) ; lower 99.%
       statx(27) = work( nwork/1000-1 ) ; lower 99.9%
        ;JBuzan END
;;;;;;;;;;;;;;;;;;;;
;;;Test_Script_For_Stat_Dispersion;;;;
begin

x = new(1000, float, 0.)
y = x
do i = 0, 999
                y(i) = 1+i ; creates an array incremental increasing values from 1-1000.
end do

print(y)
opt = False
z = stat_dispersion(y,opt)

print(z)

end
;;;;;;;;
result:
...
(22) 950
(23) 990
(24) 999
(25) 50
(26) 10
(27) 1
...
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Sep 16 12:47:51 2012