<html><head></head><body><div>Hi,</div><div>May I ask how to fill sea (lake, ...) areas differently from land. I've tried using "<a href="https://www.ncl.ucar.edu/Applications/grid_fill.shtml" rel="nofollow" target="_blank">https://www.ncl.ucar.edu/Applications/grid_fill.shtml</a>" (<span>grid_fill_5.ncl</span>) and I've got attached map.</div><div>That part which I used from grid_fill_5.ncl are:</div><div><br></div><div><span><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div>; Read/open NCL'ls crude (1x1) land-sea mask</div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div> a2 = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/landsea.nc","r")</div><div> lsdata = a2->LSMASK</div><div> lsm = landsea_mask(lsdata, xlat, xlong)</div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div> HGT2 := HGT ; replicate for plot</div><div> HGT2 = where(lsm.eq.0 , HGT2@_FillValue , HGT2)</div><div> HGT2 = where(lsm.eq.2 , HGT2@_FillValue , HGT2)</div><div> HGT2 = where(lsm.eq.4 , HGT2@_FillValue , HGT2)</div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div></span>where "HGT" is "<span>Terrain Height</span>" of WRF. The issue of the attached file is raster filled of sea areas. However is there a better way to fill sea areas?</div><div><br></div><div>Sincerely</div><div>Ehsan</div><div class="ydpf8cdb824yahoo-style-wrap" style="font-family:bookman old style, new york, times, serif;font-size:16px;"></div></body></html>