合成检验出现异常

wsx

感觉过检验区域比较异常，不知道问题出现在哪里？

wsx

1. begin
   
2.   f    = addfile ("uwnd.mon.mean.nc", "r")
   
3.   it_s = 196010
   
4.   it_e = 200110
   
5.   TIME   = f->time
   
6.   yyyy = cd_calendar(TIME,-1)
   
7.   rec_s = ind(it_s.eq.yyyy)
   
8.   rec_e = ind(it_e.eq.yyyy)
   
9.   slp_base = f->uwnd(:251,{995},:,:)
   
10.   ;slp_base = f->slp(:257,{70:90},:)
    
11.   slp_10 = f->uwnd(rec_s:rec_e:12,{995},:,:)
    
12.     dims = dimsizes(slp_10)
    
13.   slp_all =     new((/6,dims(0),dims(1),dims(2)/),typeof(slp_10))
    
14.   slp_all(0,:,:,:) = slp_10
    
15.   slp_all(1,:,:,:) = f->uwnd(rec_s+1:rec_e+1:12,{995},:,:)
    
16.   slp_all(2,:,:,:) = f->uwnd(rec_s+2:rec_e+2:12,{995},:,:)
    
17.   slp_all(3,:,:,:) = f->uwnd(rec_s+3:rec_e+3:12,{995},:,:)
    
18.   slp_all(4,:,:,:) = f->uwnd(rec_s+4:rec_e+4:12,{995},:,:)
    
19.   slp_all(5,:,:,:) = f->uwnd(rec_s+5:rec_e+5:12,{995},:,:)
    
20.   ;slp = (slp_10+slp_11+slp_12+slp_01+slp_02+slp_03)/6
    
21.   slp = reshape(slp_all,(/6*dims(0),dims(1),dims(2)/))
    
22.   copy_VarMeta(slp_base, slp)
    
23.   printVarSummary(slp)
    
24. 
    
25.   hyear=(/1961,1967,1968,1988,1992,1994,1995/)-1960             ;
    
26.   lyear=(/1969,1976,1982,1984,1989,1993,1998,1999,2000/)-1960
    
27.   
    
28.   high=new((/7,73,144/),"float")
    
29.   low=new((/9,73,144/),"float")
    
30.   
    
31. do i=0,6  
    
32. high(i,:,:) = dim_avg_n_Wrap(slp(hyear(i)*6:hyear(i)*6+5,:,:),0)
    
33. end do
    
34. 
    
35. do j=0,8
    
36. low(j,:,:)=dim_avg_n_Wrap(slp(lyear(i)*6:lyear(i)*6+5,:,:),0)
    
37. end do
    
38. 
    
39. high1=dim_avg_n_Wrap(high,0)
    
40. low1=dim_avg_n_Wrap(low,0)
    
41. 
    
42. ;x_end = high1-low1
    
43. x_end = low1 - high1
    
44. ;print(x_end)
    
45. 
    
46. copy_VarMeta(high1, x_end)
    
47. printVarSummary(x_end)
    
48. 
    
49. 
    
50. 
    
51. 
    
52. 
    
53. ;t=(Xn-Yn)/(sqrt(((n1-1)*s1^2+(n2-1)*s2^2)/(n1+n2-2))*sqrt(1/n1+1/n2))   
    
54. n1       = 7
    
55. n2       = 9               
    
56. df       = n1+n2-2        
    
57. ;ts合成差值的显著性检验
    
58. Xn_ts       =high1;样本均值
    
59. Yn_ts       =low1 ;样本均值
    
60. s1_ts       =dim_variance_n(high(:,:,:),0) ;样本方差
    
61. s2_ts       =dim_variance_n(low(:,:,:),0)  ;样本方差
    
62. s1s2_ts    =(n1-1)*s1_ts+(n2-1)*s2_ts  
    
63. t_ts       =(Xn_ts-Yn_ts)/(sqrt(s1s2_ts/df)*sqrt(1.0/n1+1.0/n2))   
    
64. prob_ts_diff    = student_t(t_ts,df)
    
65. copy_VarMeta(slp_base(0,:,:),prob_ts_diff)
    
66. printVarSummary(prob_ts_diff)
    
67. printMinMax(prob_ts_diff,True)
    
68. 
    
69. 
    
70. 
    
71. wks  = gsn_open_wks("png","uwnd")   
    
72.   gsn_define_colormap(wks,"hotcold_18lev")
    
73.   res                     = True     
    
74.    res@gsnDraw               = False
    
75.    res@gsnFrame              = False      
    
76.   res@gsnPolar            = "NH"
    
77.   res@cnFillOn            = True         
    
78.   res@mpFillOn = False
    
79.   res@cnLevelSelectionMode = "ExplicitLevels"
    
80.   res@cnLevels = (/-1,-0.5,0,0.5,1/)
    
81.   res@cnFillColors = (/6,7,10,13,14,15/)
    
82.   res@mpMinLatF  = 65.
    
83.   res@mpMaxLatF  = 90.
    
84.   res@mpCenterLonF = -90
    
85.   
    
86.   res@tiMainString    = "" ;主标题空白
    
87.   res@gsnRightString  = "" ;"unit: K/decade"右上方字符空白
    
88.   res@gsnLeftString   = ""   
    
89.    
    
90.   base = gsn_csm_contour_map_polar(wks, x_end, res)
    
91. 
    
92.      res2 = True
    
93.    res2@gsnDraw = False;
    
94.    res2@gsnFrame =False;
    
95.    
    
96.    res2@cnFillOn = True
    
97.    res2@cnLinesOn = False;
    
98.    res2@cnLineLabelsOn = False
    
99.    res2@cnInfoLabelOn = False;
    
100.    res2@lbLabelBarOn = False;
     
101.    res2@cnMonoFillPattern = False
     
102.    res2@cnLevelSelectionMode ="ExplicitLevels"
     
103.    res2@cnLevels = (/0.01/);
     
104.    res2@cnFillPatterns = (/17,-1/);
     
105.    res2@cnFillColors = (/1,-1/);
     
106.    res2@cnFillDotSizeF  = 0.003
     
107.    res2@gsnLeftString = ""
     
108.    res2@tiMainString    = ""
     
109.    res2@gsnRightString  = ""
     
110.    
     
111.    plot = gsn_csm_contour(wks,prob_ts_diff,res2)
     
112. 
     
113.    overlay(base, plot)
     
114.    draw(base)
     
115.    frame(wks)
     
116. 
     
117. 
     
118. 
     
119. 
     
120. 
     
121. 
     
122. 
     
123. 
     
124. 
     
125. 
     
126. 
     
127. 
     
128. 
     
129. 
     
130. 
     
131. 
     
132. 
     
133. 
     
134. 
     
135. 
     
136. 
     
137. 
     
138. 
     
139. 
     
140. 
     
141. 
     
142. 
     
143. 
     
144. 
     
145. 
     
146. 
     
147. 
     
148. 
     
149. 
     
150. 
     
151. 
     
152. 
     
153. 
     
154. 
     
155. 
     
156. 
     
157. end

复制代码

yesterday

prob_ts_diff    = student_t(t_ts,df)
貌似应该对prob_ts_diff 做一个筛选，比如其值小于0.05才通过了95%的显著性检验