Re: band_pass_latlon_time: currently, missing data not allowed

Hi Antonio,
You can try using linmsg or linmsg_n:
http://www.ncl.ucar.edu/Document/Functions/Built-in/linmsg.shtml
http://www.ncl.ucar.edu/Document/Functions/Built-in/linmsg_n.shtml

Note, the interpolation of missing areas (denoted by _FillValue) only
occurs over the rightmost dimension. For your case I think interpolation
over the longitude dimension is the way to go. But as Dennis warns, make
sure you interpret your results with the filling in of missing-data in
mind..
Good luck, Adam

On 09/07/2012 06:22 AM, antonio avio wrote:
> Dear NCL users,
>
> Thanks Dennis, it works. But i am not sure with the results since it
> was set that missing value is equal to zero.
> is there any recommended function that can do missing data
> interpolation for this case.
>
> Thanks,
> --Antonio
>
>
> On Fri, Sep 7, 2012 at 1:35 PM, Dennis Shea <shea@ucar.edu
> <mailto:shea@ucar.edu>> wrote:
>
> 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

-- 
______________________________________________________________
Adam Phillips                                asphilli@ucar.edu
NCAR/Climate and Global Dynamics Division       (303) 497-1726
P.O. Box 3000				
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/cas/asphilli

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Sep 10 10:03:59 2012