after filling missing value with linmsg, then mask the sea, the result showed nothing done with interpolation.

From: dyjbean <dyjbean_at_nyahnyahspammersnyahnyah>
Date: Thu Aug 08 2013 - 02:41:14 MDT

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

fill_china.PNG fill_china_1.PNG
Received on Thu Aug 8 02:42:08 2013

This archive was generated by hypermail 2.1.8 : Wed Aug 14 2013 - 18:35:09 MDT