<div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-size:small">Ehsan,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Do you mean that you want to fill the ocean in some solid color and hence mask any data in that area? If so, you can try:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">  res@mpFillOn = True</div><div class="gmail_default" style="font-size:small">  res@mpOceanFillColor = "SkyBlue"     ; use whatever color you want</div><div class="gmail_default" style="font-size:small">  res@mpLandFillColor = "Transparent"    ; to make sure land doesn't get filled</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You might also need to set some draw order resources, to make sure the filled contours and filled map areas get drawn in the right order:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">  res@cnFillDrawOrder = "Draw"     ; Also can try "predraw"</div><div class="gmail_default" style="font-size:small">  res@mpFillDrawOrder = "PostDraw"</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You can see some masking examples at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/mask.shtml">http://www.ncl.ucar.edu/Applications/mask.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">Use your browser search to look for "draworder".</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Wed, Dec 19, 2018 at 8:13 AM Ehsan Taghizadeh <<a href="mailto:ehsantaghizadeh@yahoo.com">ehsantaghizadeh@yahoo.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><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/<a href="http://landsea.nc" target="_blank">landsea.nc</a>","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="gmail-m_1321186817594108357ydpf8cdb824yahoo-style-wrap" style="font-family:"bookman old style","new york",times,serif;font-size:16px"></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>