hi,
i met a question after interpolating with linmsg.
my data is timeˇÁlatˇÁlon around the world, there exists some missing value in land surface, i want to fill those
missing value with linmsg, so i do the interpolation with three steps:
firstly, mask the sea;
secondly, interpolation
thirdly, mask the sea again, and write the interpolation result to a file.
these some key codes are as following:
yt2m = fin->TMP_GDS0_SFC({latS:latN},{lonL:lonR})
yt2m= where(ismissing(yt2m),-9999.,yt2m)
yt2m@_FillValue =-9999.
Origyt2m= yt2m
yt2m=mask(yt2m,lsm({latS:latN},{lonL:lonR}).eq.0,False)
yt2m=linmsg(yt2m,-1)
yt2m=mask(yt2m,lsm({latS:latN},{lonL:lonR}).eq.0,False)
those black circle region shows interpoltion zone.
if i comment out the second mask code, then i got the following result:
yt2m = fin->TMP_GDS0_SFC({latS:latN},{lonL:lonR})
yt2m= where(ismissing(yt2m),-9999.,yt2m)
yt2m@_FillValue =-9999.
Origyt2m= yt2m
yt2m=mask(yt2m,lsm({latS:latN},{lonL:lonR}).eq.0,False)
yt2m=linmsg(yt2m,-1)
;yt2m=mask(yt2m,lsm({latS:latN},{lonL:lonR}).eq.0,False)
the interpolation effect show better in the three circle, but the sea didn't mask.
so my question is how can i mask the sea besides good interpolation effect ?
dyjbean
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.1.8 : Wed Aug 14 2013 - 18:35:09 MDT