请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5407|回复: 0

[作图] 【求助】季风指数与降水(三维time,lat,lon)的相关的显著性检验

[复制链接]

新浪微博达人勋

发表于 2017-5-20 12:15:44 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 weety9394 于 2017-5-20 13:38 编辑

部分程序
;计算季风指数与降水相关系数**********************************************
ccr = escorc(deltamn(季风指数),prlr(降水))
ccr!0="lat"
ccr&lat=lat
ccr!1="lon"
ccr&lon=lon
ccr&lat@units="degrees_north"
ccr&lon@units="degrees_east"
printVarSummary(ccr)

;t检验
n    = dimsizes(deltamn)
print(n)
df   = n-2
t    = ccr*sqrt((n-2)/(1-ccr^2))
printVarSummary(t)      
p    = student_t(t, df)
p!0  = "lat"
p&lat= lat
p!1  = "lon"
p&lon= lon
p&lat@units="degrees_north"
p&lon@units="degrees_east"
asciiwrite("p.txt",p)            

;plot
plot1 = gsn_csm_contour_map(wks,ccr,res)     ; create plot

prres                               = True
prres@gsnDraw               = False    ; do not draw picture
prres@gsnFrame             = False    ; do not advance frame
prres@cnFillOn                = False
prres@cnLinesOn            = False
prres@cnLevelSelectionMode = "ExplicitLevels"           ; do not draw contour lines
prres@cnLevels               = fspan(0,100,101)*0.01                     ; set the contour levels
prres@cnLineLabelsOn    = False                          ; turn off contour line labels
prres@cnInfoLabelOn     = False  
plot2                              = gsn_csm_contour(wks,p,prres)

;mark
opt                        = True
opt@gsnShadeFillType       = "pattern"      ; pattern fill
;opt@gsnShadeHigh          = 17                ; use pattern #2   
;opt@gsnShadeLow           = 16        
opt@gsnShadeMid            = 17
opt@gsnLeftString          = ""
opt@gsnRightString         = ""
opt@cnLinesOn              = False            ; turn off contour lines
opt@cnLineLabelsOn         = False            ; turn off contour line labels
plot2                      = gsn_contour_shade(plot2,0,0.5,opt)

;overlay(plot1,plot2)
draw(plot1)
frame(wks)

打不出来点是为啥啊

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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