- 积分
- 5
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-3-26
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
作海冰与高度场的回归问题,需要画120°E上的垂直剖面图,填色图是回归出来的高度场异常,显著性检验的打点部分为什么只在30°N以南呢?(显著性检验的部分程序在画经纬度图的时候都是对的,放到垂直剖面图里就出现问题了),部分程序如下,请各位大佬帮忙答疑解惑~蟹蟹~
计算部分:
rc0_1 = regCoef_n(iceindex0_1djf,y_new1,0,0)
rc0_1 = rc0_1*10
copy_VarCoords(hgt0_1,rc0_1)
tval0_1 = onedtond(rc0_1@tval , dimsizes(rc0_1))
df0_1 = onedtond(rc0_1@nptxy, dimsizes(rc0_1)) - 2
b = tval0_1
b = 0.5
prob_rc0_1 = betainc(df0_1/(df0_1+tval0_1^2),df0_1/2.0,b)
prob_rc0_1 = prob_rc0_1*0.9/0.95
copy_VarCoords(hgt0_1,prob_rc0_1)
;printVarSummary(prob_rc0_1)
prob_rc0_1!0 ="level"
prob_rc0_1!1 ="lat"
prob_rc0_1&level = y_new1&level
prob_rc0_1&lat = y_new1&lat
画图部分:
res2 = True
res2@gsnDraw = False
res2@gsnFrame = False
res2@cnFillOn = True
res2@cnLinesOn = False
res2@cnLineLabelsOn = False
res2@cnInfoLabelOn = False
res2@lbLabelBarOn = False
res2@cnMonoFillPattern = False
res2@cnLevelSelectionMode = "ExplicitLevels"
res2@cnLevels = (/0.1/) ;; set to significance level
res2@cnFillPatterns = (/17,-1/)
res2@cnFillColors = (/1,0/)
res2@cnFillScaleF = 1 ; add extra density
res2@cnFillDotSizeF = 0.003
res2@gsnAddCyclic = True
res2@gsnLeftString = ""
res2@gsnLeftString = ""
res2@gsnRightString = ""
res@gsnLeftString = "(a) 1979-1997 Dec"
base(0) = gsn_csm_pres_hgt(wks,rc0_12,res)
plot(0) = gsn_csm_pres_hgt(wks,prob_rc0_12,res2)
overlay(base(0),plot(0))
|
-
|