Date: Tue Jun 04 2013 - 09:11:59 MDT
Hi Khandu,
I suggest going to gadm.org/country and selecting "Bhutan" from the pulldown menu. Click on the "Ok" button, and then click on "download". This should download several files in a subdirectory called "BTN_adm". In this directory are three levels of shapefiles.
Once you have the files, you can run the attached script to see the outlines of Bhutan
The "map" variable that is created in this script can now have contours overlaid it on:
contours = gsn_csm_contour(wks,….,res)
overlay(map,contours)
draw(map)
frame(wks)
Or, you can use gsn_csm_contour_map to create a normal plot with contours over a map, and then simply add the shapefile outlines to that map:
. . .
map = gsn_csm_contour_map(wks,….,res)
lnres = True
print("Adding adm2 lines")
lnres@gsLineColor = "green"
poly2 = gsn_add_shapefile_polylines(wks,map,"BTN_adm/BTN_adm2.shp",lnres)
print("Adding adm1 lines")
lnres@gsLineColor = "navyblue"
poly1 = gsn_add_shapefile_polylines(wks,map,"BTN_adm/BTN_adm1.shp",lnres)
print("Adding adm0 lines")
lnres@gsLineColor = "red"
poly0 = gsn_add_shapefile_polylines(wks,map,"BTN_adm/BTN_adm0.shp",lnres)
draw(map) ; This will draw everything, including the map, contours, and shapefile outlines
frame(wks)
--Mary
On Jun 3, 2013, at 10:06 PM, Khandu Khandu wrote:
> Dear All,
> I have been trying to overlay my data on the shapefile but have not succeeded so far. Since my country is too small to exist in the normal coastline map, I need some help in plotting the map. Search for help exhausted as i couldn't find a relevant one.
> Attached herewith is my ncl code.
>
> Thanking you.
>
> Khandu
> Student, Curtin University
>
>
> --
>
> ------------------------------------------------------------------------------------------------------------------------------------------------------------------
> The contents of this e-mail message and any attachments are confidential and are intended solely for addressee. The information may also be legally privileged. This transmission is sent in trust, for the sole purpose of delivery to the intended recipient. If you have received this transmission in error, any use, reproduction or dissemination of this transmission is strictly prohibited. If you are not the intended recipient, please immediately notify the sender by reply e-mail or phone and delete this message and its attachments, if any.
> ---------------------------------------------------------------------------------------------------------------------------------------------------------
> <TRMM.ncl>_______________________________________________
> 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
- application/octet-stream attachment: Bhutan.ncl
