Re: Poisson_grid_fill values

 From the example

   f = addfile ("foo.ncl" , "r")
   sst = f->SST ; (say) (time,lat,lon)

   guess = 1 ; use zonal means
   is_cyclic = True ; cyclic [global]
   nscan = 1500 ; usually much less than this
   eps = 1.e-2 ; variable dependent
   relc = 0.6 ; relaxation coefficient
   opt = 0 ; not used
;---------------------
; sst has missing values
   poisson_grid_fill( sst, is_cyclic, guess, nscan, eps, relc, opt)
; sst has no missing values

   sst_point = sst(:,72,55) ; time series at index 72,55

or, if the grid is rectilinear and coordinate variable are available

   sst_point = sst(:,{42.5},{100})

   print(sst_point)

========================

  http://www.ncl.ucar.edu/Document/Functions/Built-in/ind_resolve.shtml

   SST = f->SST ; (say) (time,lat,lon)

   sst = SST
   poisson_grid_fill( SST, is_cyclic, guess, nscan, eps, relc, opt)

   SST0 = SST(0,:,:)
   SST1d = ndtooned(SST0) ; locations with _FillValue
   i1d = ind(ismissing(SST1d))

   dimSST0 = dimsizes(SST0)
   imsg = ind_resolve(ind(ismissing(SST1d)),dimsizes(SST0))

On 7/26/12 4:39 AM, Yann Bertrand wrote:
> Hi,
> i'm using Poisson_grid_fill to plot map with missing value, the
> interolation result is pretty good.
> Can i get the value of an given lat/lon point that has been interpolated?
> I mean can i get the numerical value of a point that has been filled by
> the function?
> Thanks
>
> Yann.
>
>
>
> _______________________________________________
> 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
Received on Thu Jul 26 07:21:53 2012