Re: map = gsn_csm_map_ce(wks,res)

Yongzuo,

Are your lat/lon values of the windbarbs actually outside the bounds of the map? If so, you can try culling out the locations that are not within the map limits.

For example, if you set:

  res@mpMinLatF = minlat
  res@mpMaxLatF = maxlat
  res@mpMinLonF = minlon
  res@mpMaxLonF = maxlon

then you can use these values to cull lat/lon:

; I'm assuming lat,lon,u,v are 1D
   ii = ind(lat.ge.minlat.and.lat.le.maxlat.and.lon.ge.minlon.and.lon.le.maxlon)

  wmbarbmap(wks, lat(ii), lon(ii), u(ii), v(ii))

The wmbarbmap function doesn't know anything about the map limites of the map drawn with gsn_csm_map, so that's why you have the cull the values yourself.

--Mary

On Dec 20, 2011, at 10:13 AM, Li, Yongzuo wrote:

> Hi,
>
> I use map = gsn_csm_map_ce(wks,res)
> and wmbarbmap(wks, lat, lon, u, v)
>
> The winds outside map frame are also displayed (attached).
> How can I avoid it?
>
> Thanks.
>
> Yongzuo
> <12.GIF>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Dec 20 10:32:46 2011