关于海冰和海平面气压场回归zho遇到的问题

wsx

1. begin
   
2. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;9月海冰密集度sic(time,lat,lon)
   
3. f = addfile("ice.nc","r")
   
4. it_s = 197909
   
5. it_e = 201009
   
6. time = f->time
   
7. yyyy = cd_calendar(time, -1)
   
8. rec_s = ind(it_s.eq.yyyy)
   
9. rec_e = ind(it_e.eq.yyyy)
   
10. sic = f->sic(rec_s:rec_e:12,{66:90},:)               
    
11. sic!1 = "lat"
    
12. sic!2 = "lon"
    
13. ;printVarSummary(sic)
    
14. 
    
15. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;定义的冬季西伯利亚高压指数 slp(time)
    
16. s = addfile("slp.mon.mean.nc","r")
    
17. ti = s ->time
    
18. yyyymm = cd_calendar(ti,-1)
    
19. it_q = 197901
    
20. it_w = 201012
    
21. rec_q = ind(it_q.eq.yyyymm)
    
22. rec_w = ind(it_w.eq.yyyymm)
    
23. slp = s->slp(rec_q:rec_w,{40:60},{80:120})
    
24. slp_DJF = month_to_season(slp, "DJF")
    
25. slp_avg = dim_avg_n_Wrap(slp_DJF,(/1,2/))
    
26. 
    
27. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;进行回归
    
28. 
    
29. rc = regCoef(slp_avg, sic({lat|:},{lon|:},{time|:}))
    
30. copy_VarMeta(sic(0,:,:), rc)
    
31. print(rc)
    
32. 
    
33. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;t检验
    
34. tval = onedtond(rc@tval , dimsizes(rc))
    
35. df   = onedtond(rc@nptxy,dimsizes(rc))-2
    
36. b = tval
    
37. b = 0.5
    
38. prob = betainc(df/(df+tval^2), df/2.0, b)
    
39. prob!0   = "lat"   
    
40. prob!1   = "lon"
    
41. prob&lat = sic&lat   
    
42. prob&lon = sic&lon
    
43. prob@long_name="probability"
    
44. ;printVarSummary(prob)
    
45. 
    
46. 
    
47. 
    
48. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;画图：回归系数
    
49. 
    
50.    wks = gsn_open_wks("png","xiangguan")
    
51.    gsn_define_colormap(wks,"BlueRedGray")
    
52. 
    
53.    res                       = True
    
54.    res@gsnPolar            = "NH"
    
55.    res@gsnDraw               = False
    
56.    res@gsnFrame              = False   
    
57.    res@gsnMaximize           = True             ; make large
    
58.    res@gsnAddCyclic          = False
    
59.    
    
60.    res@cnFillOn              = True             ; turn on color
    
61.    res@cnLinesOn             = False            ; turn off contour lines
    
62.    res@cnLineLabelsOn        = False            ; turn off contour linelabels
    
63.    res@cnLevelSelectionMode  = "ManualLevels"   ; set manual contour levels
    
64.    res@cnMinLevelValF        = -0.2              ; set min contourlevel
    
65.    res@cnMaxLevelValF        =  0.2             ; set max contourlevel
    
66.    res@cnLevelSpacingF       =  0.05            ; set contourinterval
    
67.   
    
68. 
    
69.     res@mpMinLatF  =   60.
    
70.     res@mpMaxLatF  =   90.
    
71.     res@mpCenterLonF = 90
    
72.     res@mpFillOn                    = False    ; turn off map fill   
    
73. 
    
74.    res@tiMainString    = ""
    
75.    res@tiMainString    = ""
    
76.    res@gsnRightString  = ""
    
77.    base = gsn_csm_contour_map_polar(wks,rc,res)
    
78. 
    
79. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;检验区域打点
    
80.    res2 = True
    
81.    res2@gsnDraw = False;
    
82.    res2@gsnFrame =False;
    
83.    
    
84.    res2@cnFillOn = True
    
85.    res2@cnLinesOn = False;
    
86.    res2@cnLineLabelsOn = False
    
87.    res2@cnInfoLabelOn = False;
    
88.    res2@lbLabelBarOn = False;
    
89.    res2@cnMonoFillPattern = False
    
90.    res2@cnLevelSelectionMode ="ExplicitLevels"
    
91.    res2@cnLevels = (/0.05/);
    
92.    res2@cnFillPatterns = (/7,-1/);
    
93.    res2@cnFillColors = (/1,-1/);
    
94.    res2@gsnLeftString = ""
    
95.    res2@tiMainString    = ""
    
96.    res2@gsnRightString  = ""
    
97.    
    
98.    plot = gsn_csm_contour(wks,prob,res2)
    
99. 
    
100.    overlay(base, plot)
     
101.    draw(base)
     
102.    frame(wks)
     
103.   
     
104. 
     
105. 
     
106. end

复制代码

wsx

不知道问题出在哪里了？

wsx

问题已经解决了，是数据处理的问题

dh5002

楼主想问问，你数据处理哪里有问题呢

wsx

dh5002 发表于 2020-12-9 17:49
楼主想问问，你数据处理哪里有问题呢

关于季节选取的问题

caixia

楼主我最近也在研究北极海冰，都是做相关的时候总是出错，可能与海冰的缺测数据处理有关系，楼主可以指导一下海冰的缺测数据如何处理吗？我做的是冬季的海冰。万分感谢！