爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4925|回复: 3

[作图] ncl画scatter

[复制链接]

新浪微博达人勋

发表于 2018-12-7 22:06:35 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
已经画出了根据不同等级,不同大小,以及不同颜色的站点,但是在判断显著性的时候出现问题,做t检验小于0.05的站画实心,不过检验的画空心。请问家园里有没有画过有检验的scatter图。
微信图片_20181207220210.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-12-7 22:08:44 | 显示全部楼层
上图是文献里的图,也是我的追求,自己画的都是实心的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 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
  这是根据官网上的程序稍作修改的,虽然已有判断,但是画出的图是不对的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-12-14 08:38:33 | 显示全部楼层
用where函数~
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表