- 积分
- 864
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-12-10
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2015-1-4 16:31:48
|
显示全部楼层
不好意思,没注意到,回复楼上两位
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
latN = 32
latS = 22
lonL = 120
lonR = 130
nsta = 2
lat = (/各站点纬度/)
lon = (/各站点经度/)
wks = gsn_open_wks("pdf","stations")
res = True
res@gsnFrame = False ; So we can draw markers
res@gsnMaximize = True
res@mpMinLatF = 22.
res@mpMaxLatF = 32.
res@mpMinLonF = 120.
res@mpMaxLonF = 130.
res@mpDataBaseVersion = "MediumRes"
res@mpDataSetName = "Earth..4"
res@mpOutlineBoundarySets = "National"
res@mpOutlineOn = True
res@mpOceanFillColor = "lightblue"
res@mpInlandWaterFillColor= "lightblue"
res@mpLandFillColor = "tan"
res@pmTickMarkDisplayMode = "Always" ; Nicer tickmarks
map = gsn_csm_map(wks,res) ; Create the map plot, don't draw
; Draw the map (frame won't get advanced because gsnFrame was set to False).
;
;
; Draw markers on the plot in the lat/lon locations.
;
mkres = True
mkres@gsMarkerIndex = 17 ; Filled circle
mkres@gsMarkerSizeF = 0.05
map@locations = gsn_add_polymarker(wks,map,lon,lat,mkres)
draw(map)
frame(wks) ; Now advance the frame. |
|