Re: band_pass_latlon_time: currently, missing data not allowed

Users can see the code for band_pass_latlon_time via:

%> less $NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl

---
The first few lines are:
   nMsg = num( ismissing(x) )
   if (nMsg.gt.0) then
       print("band_pass_latlon_time: currently, missing data not 
allowed: nMsg="+nMsg)
       exit
   end if
---
Users can interpolate to fill missing data or, perhaps,
just set them to (say) 0.0
If the data already have x@_FillValue = -999, then
   x = where(ismissing(x), 0.0, x)   ; eliminate missing values
Users *must* interpret the results appropriately.
On 9/7/12 5:07 AM, antonio avio wrote:
> Hi NCL,
>
> It is a classic question, i have assigned the missing value in my data
> as x@_FillValue = -999,
> but i still get this error:
> *(0)band_pass_latlon_time: currently, missing data not allowed: nMsg=709920*
>
> Any suggestion?
>
> Variable: x
> Type: float
> Total Size: 583626240 bytes
>              145906560 values
> Number of Dimensions: 3
> Dimensions and sizes:[time | 13880] x [lat | 73] x [lon | 144]
> Coordinates:
>              time: [17298624..17631720]
>              lat: [90..-90]
>              lon: [ 0..357.5]
> Number Of Attributes: 13
>    missing_value :32766
>    valid_range :(  0, 500 )
>    _FillValue :-999
>    parent_stat :Individual Obs
>    statistic :Mean
>    level_desc :Other
>    dataset :NOAA Interpolated OLR
>    precision :2
>    var_desc :Outgoing Longwave Radiation
>    actual_range :( 84, 330.25 )
>    unpacked_valid_range :(  0, 500 )
>    long_name :Anomalies: Daily OLR
>    units :W/m^2
>
> Thanks,
> Antonio
>
>
> _______________________________________________
> 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