- 积分
- 541
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-10-19
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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)
打不出来点是为啥啊
|
|