[ncl-talk] how to mask out tibetan plateau when plotting fnl data?

[Dennis Shea]

 You must provide more information!

ncl-talk does not know what 'FNL' is. A GRIB file?
What data do you have available?
What type of grid do you have .... curvilinear or rectilinear?

If you have (say) surface geopotential (m) and you are on a curvilinear
grid ... one approach

   f = addfile(...)
   x      = f->X         ; (level,nlat,mlon) ?????
   phis = f->PHIS   ; (nlat,mlon)

  lat2d   = f->
  lon2d  = f->

   lvmask = (/700, 850 /)

  do kl=1,dimsizes(lvmask)-1
       x({lvmask},:,:) = where(lat2d.gt. 25 .and. lat2d.lt. 42   .and. \
                                           lon2d.gt. 70 .and. lon2d.lt.110
.and. \
                                           phis.gt.1500, x at _FillValue,
x({lvmask},:,:))
  end do


On Sun, Aug 31, 2014 at 6:48 AM, Li Qi <liqi123sh at qq.com> wrote:

> Dear all,
>
> I'd like to ask how to mask out tibetan plateau when plotting FNL data on
> 700hPa and 850hPa.
>
> Any advice will be greatly appreciated!
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140901/e3770b59/attachment.html