question about linmsg interpolating missing value in a region that result has not changed as before ?

From: dyjbean <dyjbean_at_nyahnyahspammersnyahnyah>
Date: Mon Nov 25 2013 - 19:16:39 MST

hi,
  when i processed the missing value with linmsg interpolation method ,there met a question as follows:
the part of codes is as follows:

______________________________________________________________
    lat!0 = "lat"
    lat@long_name="lat"
    lat@units="degrees-north"
    lat&lat=lat

    lon!0 = "lon"
    lon@long_name="lon"
    lon@units="degrees-east"
    lon&lon=lon

    slrad!0="time"
    slrad!1="lat"
    slrad!2="lon"
    slrad&lat=lat
    slrad&lon=lon
    slrad&lat@units="degrees-north"
    slrad&lon@units="degrees-east"

    slrad=where(ismissing(slrad),2.e+20,slrad)
    slrad@_FillValue=2.e+20
    slrad@missing_value=2.e+20
    print(num(ismissing(slrad))+" before")
;; linear interpolation
     sslrad=linmsg(slrad(lat|:,lon|:,time|:),-1)
     sslrad!0="lat"
     sslrad!1="lon"
     sslrad!2="time"
     slrad=(/sslrad(time|:,lat|:,lon|:)/)
     print(num(ismissing(sslrad))+" temp")
     print(num(ismissing(slrad))+" after")
     delete(sslrad)
________________________________________________________________

the output result in screen is as follows

(0) 65472 before
(0) 65472 temp
(0) 65472 after


it seemed the linear interpolation didn't work.
it's strange appearance




dyjbean

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Nov 25 19:16:33 2013

This archive was generated by hypermail 2.1.8 : Wed Dec 04 2013 - 20:42:38 MST