[ncl-talk] Attach shapefile polylines to contour
Rick Brownrigg
brownrig at ucar.edu
Thu May 26 13:38:58 MDT 2016
Hi,
wrf_map_overlays has some resources specific to it, one of which controls
"framing" similar to gsnFrame. You might set:
pltres at FramePlot = False ; default is True
For more details, see:
http://ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml
Note that it says "...the plot will be drawn, but the frame will not
advanced." I'm totally sure my suggestion will work -- it seems like the
calls to
gsn_add_shapefile_polylines() should be made before the plot is drawn. But
perhaps give it a try.
Hope that helps...
Rick
On Thu, May 26, 2016 at 11:06 AM, Amadou Coulibaly <mpapin24 at gmail.com>
wrote:
> Dear NCL users,
>
> I want to attach shapefile polylines to my wrf_contour plot, but I am
> facing an issue. After downloading the country data from Global
> Administrative Area, I wrote the command lines below:
> ;---Start the graphics
> wks = gsn_open_wks("eps","wrf_latlonoverlay_16_ERA")
>
> ;---Set some basic plot options
> opts = True
> opts at MainTitle = "REAL-TIME WRF"
> opts at gsnDraw = False
> opts at gsnFrame = False
> ;opts at gsnAddCyclic = False
>
> pltres = True
> mpres = True
> mpres at mpGeophysicalLineColor = "Black"
> mpres at mpNationalLineColor = "Black"
> mpres at mpUSStateLineColor = "Black"
> mpres at mpGridLineColor = "Black"
> mpres at mpLimbLineColor = "Black"
> mpres at mpPerimLineColor = "Black"
> mpres at mpGeophysicalLineThicknessF = 2.0
> mpres at mpGridLineThicknessF = 2.0
> mpres at mpLimbLineThicknessF = 2.0
> mpres at mpNationalLineThicknessF = 2.0
> mpres at mpUSStateLineThicknessF = 2.0
>
> opts at TimeLabel = times ; Set valid time to use on plots
>
> pressure = pressure_levels(0)
>
> ; Add some level info to the plot
> opts at PlotLevelID = pressure + " hPa"
> ;----------------------------------------------------------------------
> ; Plot full domain first.
> ;----------------------------------------------------------------------
> opts at cnFillOn = True
>
> contour = wrf_contour(f,wks,Rot_avg(0,:,:),opts)
>
> x = (/3., 4., 5., 6., 7./)
> y = (/3., 4., 5., 6., 7./)
>
>
> ;-- polymarker resources
> pmres = True
> pmres at gsMarkerColor = "red" ;-- marker color
> pmres at gsMarkerSizeF = 0.03 ;-- set size of marker
> pmres at gsLineThicknessF = 13. ;-- marker line thickness
>
> ;-- draw all 16 markers on map using unique identifier name
> ;-- and additional map attribute settings
> x = 3. ;-- x-position of first marker
> y = 3. ;-- y-position of first marker
> do i = 5,7,2 ;-- 16 different marker
> pmres at gsMarkerIndex = i+1
> str = unique_string("poly") ;-- result is poly0-poly15
> contour@$str$ = gsn_add_polymarker(wks, contour, x+(i*11.3),
> y+(i*8.5), pmres)
>
>
> ;-- add marker to map
> end do
>
> plot = wrf_map_overlays(f,wks,(/contour/),pltres,mpres)
>
> ;---Attach shapefile polylines to map
> dir1 = "CIV_adm/"
> dir2 = "GHA_adm/"
> dir3 = "TGO_adm/"
> dir4 = "BEN_adm/"
> dir5 = "NGA_adm/"
>
> filename1 = "CIV_adm0.shp"
> filename2 = "GHA_adm0.shp"
> filename3 = "TGO_adm0.shp"
> filename4 = "BEN_adm0.shp"
> filename5 = "NGA_adm0.shp"
>
>
> lnres = True
> lnres at gsLineThicknessF = 1.0
>
> poly1 = gsn_add_shapefile_polylines(wks,plot,dir1+filename1,lnres)
> poly2 = gsn_add_shapefile_polylines(wks,plot,dir2+filename2,lnres)
> poly3 = gsn_add_shapefile_polylines(wks,plot,dir3+filename3,lnres)
> poly4 = gsn_add_shapefile_polylines(wks,plot,dir4+filename4,lnres)
> poly5 = gsn_add_shapefile_polylines(wks,plot,dir5+filename5,lnres)
>
> draw(plot)
> frame(wks)
> end
>
> But instead of overlaying country map with contour, I am getting two
> separate plots one with contour and another one with country maps, as
> showing attached plot.
>
> What,s wrong with my script? And how to put vertically the labelbar?
>
> Best regards
>
> --
> * COULIBALY AMADOU *
> PhD Student on West African Climate System (WACS)
> FUTA - Federal University of Technology of Akure, Nigeria
>
> *Visiting Student - *
> *University of Cologne, Germany**Institute of Geophysics and Meteorology*
> Pohligstr. 3 / Office 3.102
> D-50969 Köln
>
> *Project*: WASCAL (West African Science Service Centre on Climate Change
> and Adapted Land Use)
> Phone:(+234) 812 965 3659 /(+223) 90 08 26 71 /(+226) 64 57 37 27 /+49
> 15218352574
>
> E-mail: mpapin24 at gmail.com / coulibalya68 at yahoo.com
>
> "*The time is always right to do right": Nelson Mandela*
>
> *"Character is like a tree and reputation like a shadow. The shadow is
> what we think of it; the tree is the real thing" : Abraham Lincoln*
>
> *"Do what you can, with what you have, where you are" Theodore Roosevelt*
>
>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> 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/20160526/6042a9f1/attachment.html
More information about the ncl-talk
mailing list