- 积分
- 5576
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-10-7
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2018-12-7 22:11:31
|
显示全部楼层
就差最后一步,已有显著性检验结果的txt
sigval = 0.05
do i=0,nlevels-2
ii = ind(levels(i).le.dvals.and.dvals.lt.levels(i+1))
;print(ii)
mkres@gsMarkerIndex = 16 ; Filled dots
mkres@gsMarkerSizeF = sizes(i)
mkres@gsMarkerThicknessF = 2.5
if(x(i).lt.sigval)then
mkres@gsMarkerColor = colors(i)
else
mkres@gsMarkerColor = "white"
end if
dum_fill(i) = gsn_add_polymarker(wks,map,lon(ii),lat(ii),mkres)
;---Draw the markers again, only this time hollow, to get an outline.
mkres@gsMarkerIndex = 4 ; Hollow dots
mkres@gsMarkerColor = colors(i) ;"black"
dum_hollow(i) = gsn_add_polymarker(wks,map,lon(ii),lat(ii),mkres)
delete(ii)
end do
这是根据官网上的程序稍作修改的,虽然已有判断,但是画出的图是不对的 |
|