[ncl-talk] landseamask

[Ehsan Taghizadeh]

Hi,May I ask how to fill sea (lake, ...) areas differently from land. I've tried using "https://www.ncl.ucar.edu/Applications/grid_fill.shtml" (grid_fill_5.ncl) and I've got attached map.That part which I used from grid_fill_5.ncl are:
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; Read/open NCL'ls crude (1x1) land-sea mask;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  a2    = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/landsea.nc","r")  lsdata = a2->LSMASK  lsm  = landsea_mask(lsdata, xlat, xlong)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     HGT2  := HGT                                  ; replicate for plot     HGT2  = where(lsm.eq.0 , HGT2 at _FillValue , HGT2)     HGT2  = where(lsm.eq.2 , HGT2 at _FillValue , HGT2)     HGT2  = where(lsm.eq.4 , HGT2 at _FillValue , HGT2);;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
where "HGT" is "Terrain Height" of WRF. The issue of the attached file is raster filled of sea areas. However is there a better way to fill sea areas?
SincerelyEhsan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181219/08454c93/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: WRF_pcp_2_2018102112_24.png
Type: image/png
Size: 429465 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181219/08454c93/attachment-0001.png>