Re: dim_stat4 and missing_values

Jamie and others,

I have updated the stat2, stat4, dim_stat4, stat_medrng, and stat_trim
functions/procedures so that if 'x' contains a subarray that is
all missing, then missing values will be returned in the appropriate
locations in the output, rather than a fatal error occurring.

These updates will be available in the next release of NCL, which is
not yet scheduled. Send me a personal email if you need this
functionality now.

Cheers,

--Mary

On Fri, 31 Aug 2007, Jamie Scott wrote:

> I have a 3D array dimensioned x(lat,lon,time) that I'm passing to dim_stat4.
> The array has
> all values over land set to missing_value (and _FillValue). When the array
> is passed to dim_stat4 (or stat4),
> I get an error: "fatal:stat4: The first input array contains all missing
> values."
>
> Apparently it's complaining about the land points where all time values are
> set to _FillValue. I had hoped the function
> would return _FillValue for these points and compute the statistics over the
> ocean grid points that are not set to _FillValue. :(
>
> So I'd like to pass the non-land points only in a 2D array (npts,time) to
> dim_stat4 and take the resulting array (4,npts) and put it
> back into the lat/lon grid (4,lat,lon) in the correct non-missing (ocean
> points) locations. I figure I need to use ndtooned, ind, ind_resolve
> and onedtond in some capacity, but am struggling with this problem. Does
> anyone have any suggestions?
>
> -Thanks,
>
>
> Jamie Scott
> NOAA/ESRL/PSD
> james.d.scott_at_noaa.gov
>
>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 05 2007 - 08:48:50 MDT