[ncl-talk] _FillValue can't work

[Dennis Shea]

This is a bug. A JIRA ticket ( has been opened:     [JIRA] (NCL-2259)

===
a possible work around.

 ;===

 undef("stddev_ubyte")
 function stddev_ubyte(aub:ubyte)
 begin
   ;print(aub)
   ;print(typeof(aub at _FillValue))
   ;print(avg(toint(aub)))
    return(stddev(toint(aub)))
 end
 ;===


 a=(/94B,95B,253B,92B,95B,94B,253B,253B,91B,253B,96B,253B,253B,253B/)
 print(a)
 print(stddev(a))
 a at _FillValue=253B
 print(a)
 print(stddev(a))


  astd =  stddev_ubyte(a)


On Sat, Aug 1, 2015 at 5:40 PM, Sitan Heyl <heylsitan at gmail.com> wrote:
> Hi, all
> I have a problem about fillvalue. the code is as follow. I don't know why
> the missing value setting can't work. any help is welcome!
> ncl 0> a=(/94B,95B,253B,92B,95B,94B,253B,253B,91B,253B,96B,253B,253B,253B/)
> ncl 1> print(a)
> ncl 2> print(stddev(a))
> ncl 3> a at _FillValue=253B
> ncl 4> print(stddev(a))
>
> Yongli
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>