<div dir="ltr"><div><div><div><div><div><div> You must provide more information! <br><br>ncl-talk does not know what 'FNL' is. A GRIB file?<br>What data do you have available? <br></div><div>What type of grid do you have .... curvilinear or rectilinear?<br>
</div><div><br></div>If you have (say) surface geopotential (m) and you are on a curvilinear grid ... one approach<br><br></div> f = addfile(...)<br></div> x = f->X ; (level,nlat,mlon) ?????<br></div> phis = f->PHIS ; (nlat,mlon)<br>
<br></div><div> lat2d = f-><br></div><div> lon2d = f-><br></div><div><br></div><div> lvmask = (/700, 850 /)<br></div><div><br></div> do kl=1,dimsizes(lvmask)-1<br></div><div> x({lvmask},:,:) = where(<a href="http://lat2d.gt">lat2d.gt</a>. 25 .and. <a href="http://lat2d.lt">lat2d.lt</a>. 42 .and. \<br>
<a href="http://lon2d.gt">lon2d.gt</a>. 70 .and. lon2d.lt.110 .and. \<br></div><div> phis.gt.1500, x@_FillValue, x({lvmask},:,:))<br></div>
end do <br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sun, Aug 31, 2014 at 6:48 AM, Li Qi <span dir="ltr"><<a href="mailto:liqi123sh@qq.com" target="_blank">liqi123sh@qq.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>Dear all,</div><div><br></div><div>I'd like to ask how to mask out tibetan plateau when plotting FNL data on 700hPa and 850hPa.</div>
<div><br></div><div>Any advice will be greatly appreciated!</div><br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>