bug in dim_sum_wgt function

Hi

When I use following script to test the function of dim_sum_wgt, but
I found that missing value setting not work for this function, we
always got the same result as red one following, regardless opt = 0
or 1.
Please give me some useful suggestion, thanks.

Jun

;=======================================================================
=====;
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;=======================================================================
=====;
begin

   x = (/(/2, 4, 6/),(/9.96921e+36, 2, 3/)/)
   y = (/1, 2, 9.96921e+36/)

   x_at_missing_value = 9.96921e+36
   x@_FillValue = 9.96921e+36
   y@_FillValue = 9.96921e+36
   y_at_missing_value = 9.96921e+36

   z = dim_sum_wgt(x, y, 1)

   print(x)
   print(y)
   print(z)
end

Variable: z
Type: float
Total Size: 8 bytes
             2 values
Number of Dimensions: 1
Dimensions and sizes: [2]
Coordinates:
Number Of Attributes: 1
   _FillValue : 9.96921e+36
(0) 5.981526e+37
(1) 2.990763e+37

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 30 2008 - 17:14:13 MDT