Interpolation fault of Remapping pop grid to regular grid

> Hi, all

When I use remap function to transport pop grid to regular grid,
there will be some fault interpolation value at boundary, as shows of
the picture attached. Does someone ever faced the same question?

Scripts:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/popRemap.ncl"
begin

a = addfile("b30.19_0kaDVTc.pop.h.nc","r")

print("1")

var = a->PD(0:24,:,:,:)

remap = PopLatLon(var,"gx3v5","1.9x2.5","bilin","da","060511")
delete(var)

out=addfile("b30.19_0kaDVTc.pd_1.9x2.5_test.nc","c")
;filedimdef(out,"time",-1,True)

print("4")

out->PD=remap

end
Picture:

Jun Cheng

CCR, UW-Madison

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk