[ncl-talk] nan in data radar set, how to handle

[Dennis Shea]

"The "isnan_ieee" didnt work"

The following worked fine.


f = addfile("COMP_2level.nc","r")
r = f->reflectivity
nNaN = num(isnan_ieee(r))             ; nNaN=76871
print(nNaN)

r = where(isnan_ieee(r), r at _FillValue, r)
kNaN = num(isnan_ieee(r))              ; kNaN= 0
print(kNaN)


On Sat, Nov 21, 2015 at 10:42 PM, Geeta Geeta <geetag54 at yahoo.com> wrote:
> Hi Rick.
> Thanks for your suggestion.
> The "isnan_ieee" didnt work. but I converted the double to float and then it
> worked.
>
> refl    = doubletofloat(ref)
>   reflec  = where(refl.eq."nan",-9999.,refl)
>              asciiwrite("reflec.txt",reflec)
> ;            printVarSummary(reflec)
> This is what I did.
>
> But I still have one question. when the dataset is of the type "double",
> then how will "nan" be represented??? Can You pls explain???
>
> The Figure is attached.
> Geeta.
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>