I just want to put country boundaries on a vector plot, but the
resource mpOutlineBoundarySets = "National" or "allBoundaries" didn' t
show any line.
And with gsn_csm_vector_scalar_map_ce, all the continenta are colored
gray, without any line. How can I plot the data over the continent?
That is my code:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
  a = addfile("DT53P008-FMA.nc","r")
;************************************************
; read in zonal [u] and meridional [v] winds (July)
;************************************************
  u = a->U(0,:, :,:)
  v = a->V(0,:, :,:)
;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("pdf","vector-FMA")                ; open a ps file
  plot = new(18,graphic)			
  vcres                     = True              ; plot mods desired
  vcres@gsnDraw             = False       ; Do not draw plot
  vcres@gsnFrame            = False       ; Do not advance frome
  vcres@vcRefAnnoOrthogonalPosF = -1.0             ; move ref vector up
  vcres@vcRefMagnitudeF         = 10.0             ; define vector ref mag
  vcres@vcRefLengthF            = 0.045            ; define length of vec ref
  vcres@vcGlyphStyle            = "CurlyVector"    ; turn on curly vectors
  vcres@vcMinDistanceF          = 0.017
  vcres@mpOutlineBoundarySets = "allBoundaries"
  vcres@pmTickMarkDisplayMode  = "Always"           ; turn on tickmarks
  vcres@mpMaxLatF                   = 9
  vcres@mpMinLatF                   = -23
  vcres@mpMaxLonF                   = -30
  vcres@mpMinLonF                   = -85
do n=0, 17
  plot(n) = gsn_csm_vector_map_ce(wks,u(n, :,:),v(n, :,:), vcres)
end do
 gsn_panel(wks,plot,(/6,3/),False)
end
Thanks,
Varejão
On 6/26/06, Adam Phillips <asphilli_at_cgd.ucar.edu> wrote:
> Hi Varejão,
>
> You might want to try using gsn_csm_vector_scalar_map_ce:
>
> http://www.ncl.ucar.edu/Applications/vector.shtml#ex1
>
> If you just want to put country boundaries on a vector plot, you can set
> the resource mpOutlineBoundarySets = "National"
>
> If this didn't answer your question let us know..
> Adam
>
> Varejão Júnior wrote:
> > Hi,
> >
> > How can I plot vectors with country boundaries?
> > There are a way to overlay a map plot and a vector plot?
> > I've tried to use overlay(plota, plotb), with plota =
> > gsn_csm_vector_map_ce , and plotb = gsn_csm_contour_map_ce, but it
> > doesn't works.
> >
> > Thanks,
> >
> > Varejão
> >
>
> --
> --------------------------------------------------------------
> Adam Phillips                                asphilli_at_ucar.edu
> National Center for Atmospheric Research   tel: (303) 497-1726
> ESSL/CGD/CAS                               fax: (303) 497-1333
> P.O. Box 3000
> Boulder, CO 80307-3000    http://www.cgd.ucar.edu/cas/asphilli
>
-- ###################################### Claudeci Gomes Varejão Júnior Computer Science Federal University of Viçosa Viçosa Minas Gerais Brasil http://www.dpi.ufv.br http://madeira.dea.ufv.br ###################################### "Que homem é o homem que não torna o mundo melhor?" _______________________________________________ ncl-talk mailing list ncl-talk_at_ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Mon Jun 26 2006 - 14:54:31 MDT
This archive was generated by hypermail 2.2.0 : Wed Jul 05 2006 - 09:26:57 MDT