Hi,
I want to add a marker to china map at 44.09N,113.57E. But I don't know why gsn_add_polymarker doesn't work. Can anyone help me to look at it? Thanks.
Guo Zhang
wks = gsn_open_wks("ps","china")
res = True
res@mpDataBaseVersion ="MediumRes"
res@mpOutlineOn = True ; Turn on map outlines
res@gsnMaximize = True ; Maximize plot in frame.
res@gsnPaperOrientation = "Portrait"
res@gsnFrame = False
res@gsnDraw = False
res@vpHeightF = 0.55 ; Changes the aspect ratio
res@vpWidthF = 0.65
res@mpDataSetName = "Earth..4" ; This database contains
; divisions for other countries. needed ncl5.1.0 and later
res@mpOutlineSpecifiers = (/"China:states","Taiwan"/)
res@mpOutlineBoundarySets = "NoBoundaries"
res@mpMinLatF = 17. ; set the minimum latitude = 17. domain of china
res@mpMaxLatF = 55. ; set the maximum latitude = 55.
res@mpMinLonF = 72. ; set the minimum latitude = 72.
res@mpMaxLonF = 136. ; set the maximum latitude = 136.
res@mpLandFillColor = "White" ;land color
res@mpOceanFillColor = "White"
res@mpInlandWaterFillColor = "White"
base_map = gsn_csm_map_ce(wks,res)
chinariver = add_china_river(wks,base_map,"blue",1,0) ; Add rivers
polyres = True
polyres@gsMarkerIndex = 16 ; polymarker style
polyres@gsLineColor = "red"
; polyres@gsMarkerSizeF = 0.02 ; polymarker size
pump5 = gsn_add_polymarker(wks,base_map,44.09,113.57,polyres)
draw(base_map)
frame(wks)
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 4 15:08:13 2012
This archive was generated by hypermail 2.1.8 : Wed Jun 06 2012 - 15:17:44 MDT