2 things:
1) overlay the plot on the map, not on the polylines that were added to the map.
overlay(map,plot)
2) don't use gsn_csm_contour_map_ce. This includes a map of its own. Instead use gsn_csm_contour.
An alternative is to use gsn_csm_contour_map_ce but add the polylines to the plot id that it returns:
plot = gsn_csm_contour_map_ce(wks,dat2d_1(:,:),res)
poly = gsn_add_shapefile_polylines_subset(wks,plot,shp_filename,"ROUTE",\
highways,colors,pres)
You might need to play around with the draw order resources: tfPolyDrawOrder and cnFillDrawOrder
-dave
On Apr 9, 2013, at 3:40 PM, "Luo, Chao" <chao.luo@eas.gatech.edu> wrote:
> Hi,
>
> I am try to plot pollutants conc. over map (with highways on map)in GA. The pollutants data are two dimension grids data. What I did in modified shapefiles_12.ncl script are:
> First make poly map by poly = gsn_add_shapefile_polylines_subset(wks,map,shp_filename,"ROUTE",\
> highways,colors,pres)
> second read two dimensional data, and lat/lon, make plot by plot = gsn_csm_contour_map_ce(wks,dat2d_1(:,:),res)
>
> Then overlay two plots by overlay(poly,plot).
>
> But it gives error: fatal:Number of elements of dimension (0) of argument (0) is (464) in function (overlay), expected (1) elements. That means overlay function requires same dimension of poly and plot?
>
> Thanks much,
>
> cl
>
>
>
>
>
>
> _______________________________________________
> 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 Apr 9 16:24:10 2013
This archive was generated by hypermail 2.1.8 : Mon Apr 15 2013 - 20:12:25 MDT