Re: remove wind vector outside of my domain box?

My apologies, I don't know if this was ever answered.

The wmvectmap is a rather primitive routine and doesn't have a lot of options for control.

I believe the only want to cull out the vectors you don't want is to set the vectors outside your
domain to missing values.

For example, if you want to cull the vectors that fall outside the range of lat 36-42
and lon -82 to -77, then try something like this (untested):

  lat1d = ndtooned(lat)
  lon1d = ndtooned(lon)
  U1d = ndtooned(U)
  V1d = ndtooned(V)
  U1d@_FillValue = 1e20
  V1d@_FillValue = 1e20
  U1d = where(lat.ge.36.and.lat.le.42.and.lon.ge.-82.and.lon.le.-77, \
              U1d,U1d@_FillValue)
  V1d = where(lat.ge.36.and.lat.le.42.and.lon.ge.-82.and.lon.le.-77, \
              V1d,V1d@_FillValue)
  wmvectmap(wks, lat1d, lon1d, U1d, V1d) ; Plot barbs.

If your lat, lon, U, V arrays are already 1D, then you don't need to convert them to 1D first.

--Mary

On Oct 15, 2010, at 12:45 PM, xiaoming Hu wrote:

> Hello
>
> I used wmvectmap plot some observed wind vectors over my model domain.
>
> See the attached plot and script. Apparently, I got wind vectors outside my domain box. How can I remove them?
>
> I just want to show the vectors within my domain.
>
> Thanks
>
> Xiaoming
> <plotMADIS_dataWindVector.ncl><wrfout_d03_WindVector_1.png>_______________________________________________
> 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 Thu Oct 21 09:59:15 2010