Hi Varejăo,
Now I understand, what you want to use is gsn_csm_vector_map_ce, and
then set the following resources:
vcres@mpFillOn = False ; turn the gray map fill off
vcres@mpGeophysicalLineThicknessF = 2.0 ; optional
vcres@mpOutlineBoundarySets = "National"
Good luck.
Adam
Varejăo Júnior wrote:
> 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
>>
>
>
>
-- -------------------------------------------------------------- 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 _______________________________________________ ncl-talk mailing list ncl-talk_at_ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Mon Jun 26 2006 - 15:00:14 MDT
This archive was generated by hypermail 2.2.0 : Wed Jul 05 2006 - 09:26:57 MDT