Re: Stat_Dispersion

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Tue Sep 18 2012 - 06:34:16 MDT

This will be in the 6.1.0 release.

The order is slightly different from yours.

http://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml

THX

On 9/16/12 12:47 PM, jbuzan wrote:
> 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
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Sep 18 06:34:32 2012

This archive was generated by hypermail 2.1.8 : Fri Sep 21 2012 - 16:22:30 MDT