Re: pollutants distribution over traffic map

From: David Brown <dbrown_at_nyahnyahspammersnyahnyah>
Date: Tue Apr 09 2013 - 16:24:00 MDT

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