[ncl-talk] missing value in int2p

Dennis Shea shea at ucar.edu
Fri May 15 13:02:01 MDT 2015


Actually, I should have added to the documentation that _FillValue are not
allowed.

---

A simple work around

  p = (/-999, 750, 500, -999, 200, 100,  50,  10/)

  ip = ind(.not.ismissing(p))

  lin = int2p(p(ip),t(ip),prs,1)

Cheers





On Fri, May 15, 2015 at 11:41 AM, Joowan Kim <joowan at ucar.edu> wrote:

> Dear NCL people
>
> I'm using int2p for interpolation.
> I generally work well with missing value, but if the first coordinate
> value is missing it returns all missing values.
> I can detour this problem, but it will be good to know why.
> Any idea or suggestion with this?
>
> Best!
> Joowan
>
> Test script
> -----------------------------------
> ; it understands fill value
> p=(/-999, 750, 500, -999, 200, 100,  50,  10/)
> ;p=(/1000, 750,-999, -999, 200, 100,  50,  10/)
> t=(/ 320, 300, 250, -999, 220, 190, 210, 280/)
> p at _FillValue = -999
> t at _FillValue = -999
>
> prs = (/600,300,30/)
> tlin = int2p(p,t,prs,1)
> tlog = int2p(p,t,prs,2)
> print(tlin)
> print(tlog)
> ----------------------------------------
>
> --
> Joowan Kim
> Atmospheric Chemistry Division
> National Center for Atmospheric Research
> PO Box 3000  Boulder, CO 80307
> e-mail : joowan at ucar.edu
>            joowan.k at gmail.com
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150515/6181a455/attachment.html 


More information about the ncl-talk mailing list