为啥经纬度显示有问题呢

sususir

1. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
4. begin
   
5. 
   
6. ;;;;;;;;;;;;;;;500hPa合成分析;;;;;;;;;
   
7. ;;;;;;;;;;;;;;;弱:1913—1926,强1930--1950;;;;;;;;
   
8. f=addfile("/home/a2/hgt500.qihoutai.1920_1949.nc","r")
   
9. hgt=f->sst
   
10. f0=addfile("/home/a2/H.hgt500.1930_1950.nc","r")
    
11. hgth=f0->sst
    
12. f1=addfile("/home/a2/L.hgt500.1913_1926.nc","r")
    
13. hgtl=f1->sst
    
14. 
    
15. ;;;;;;;;合成;;;;;;
    
16. avel=dim_avg_n_Wrap(hgtl,0)
    
17. aveh=dim_avg_n_Wrap(hgth,0)
    
18. cha=new((/91,180/),float)
    
19. do j=0,179
    
20. do i=0,90
    
21. cha(i,j)=aveh(i,j)-avel(i,j)
    
22. end do
    
23. end do
    
24. aveh!0="lat"
    
25. aveh!1="lon"
    
26. aveh&lat=fspan(-90,90,91)
    
27. aveh&lon=fspan(0,358,180)
    
28. avel!0="lat"
    
29. avel!1="lon"
    
30. avel&lat=fspan(-90,90,91)
    
31. avel&lon=fspan(0,358,180)
    
32. cha!0="lat"
    
33. cha!1="lon"
    
34. cha&lat=fspan(-90,90,91)
    
35. cha&lon=fspan(0,358,180)
    
36. hgt!0="lat"
    
37. hgt!1="lon"
    
38. hgt&lat=fspan(-90,90,91)
    
39. hgt&lon=fspan(0,358,180)
    
40. 
    
41. ;;;;;;;;;;合成分析显著性检验;;;;;;;
    
42. ;;;;;0.05检验临界为2.035;;;;
    
43. varh=dim_variance_n_Wrap(hgth,0)
    
44. varl=dim_variance_n_Wrap(hgtl,0)
    
45. vah=where(varh.eq.0,varh@_FillValue,varh)
    
46. val=where(varl.eq.0,varl@_FillValue,varl)
    
47. nl=14.
    
48. nh=21.
    
49. n=nl+nh-2.
    
50. t=new((/91,180/),float)
    
51. do j=0,179
    
52. do i=0,90
    
53. t(i,j)=(cha(i,j)/(sqrt(((nh-1)*vah(i,j)+(nl-1)*val(i,j))/n)))/(sqrt(1./nl+1./nh))
    
54. end do
    
55. end do
    
56. t!0="lat"
    
57. t!1="lon"
    
58. t&lat=fspan(-90.,90.,91)
    
59. t&lon=fspan(0.,358.,180)
    
60. 
    
61. ;;;;;;;;;;;;;;;绘图;;;;;;;
    
62. wks = gsn_open_wks("pdf","hgt500_hc")
    
63. plot=new(4,graphic)
    
64. gsn_define_colormap(wks,"BlueDarkRed18")
    
65. res = True
    
66. res@gsnDraw=False
    
67. res@gsnFrame = False
    
68. res@gsnAddCyclic = True
    
69. res@cnFillOn = False
    
70. res@pmTickMarkDisplayMode = "Always"
    
71. res@gsnLeftString = ""
    
72. res@gsnRightString = ""
    
73. res@mpCenterLonF = 180
    
74. res@mpFillOn= True
    
75. res@cnLineLabelsOn = True
    
76. res@cnLineLabelBackgroundColor = "-1"
    
77. res@mpGridAndLimbOn = False
    
78. ;res@mpGridSpacingF = 10
    
79. ;res@mpGridLineThicknessF = 0.2
    
80. res@gsnMaximize = True
    
81. 
    
82. ;;;;;;;;;;;;气候态图;;;;
    
83. res1 = res
    
84. res1@cnLevelSelectionMode="ManualLevels"
    
85. ;res1@cnLinesOn = False
    
86. res1@cnMinLevelValF = 400
    
87. res1@cnMaxLevelValF = 600
    
88. res1@cnLevelSpacingF = 4
    
89. ;res1@gsnSpreadColorStart = 4
    
90. res1@gsnContourNegLineDashPattern = 2
    
91. res1@gsnContourPosLineDashPattern = 0
    
92. plot(0)=gsn_csm_contour_map(wks,hgt/10,res1)
    
93. 
    
94. ;;;;;;;;;;;;高
    
95. res2=res
    
96. res2@cnLevelSelectionMode="ManualLevels"
    
97. ;res2@cnLinesOn = False
    
98. res2@cnMinLevelValF = -3
    
99. res2@cnMaxLevelValF = 3
    
100. res2@cnLevelSpacingF = .3
     
101. ;res2@gsnSpreadColorStart = 4
     
102. res2@gsnContourNegLineDashPattern = 2
     
103. res2@gsnContourPosLineDashPattern = 0
     
104. plot(1)=gsn_csm_contour_map(wks,aveh/10,res2)
     
105. 
     
106. ;;;;;;;;;低;;;;;;;;;
     
107. res3=res
     
108. res3@cnLevelSelectionMode="ManualLevels"
     
109. ;res3@cnLinesOn = False
     
110. res3@cnMinLevelValF = -3
     
111. res3@cnMaxLevelValF = 3
     
112. res3@cnLevelSpacingF = .5
     
113. ;res3@gsnSpreadColorStart = 4
     
114. res3@gsnContourNegLineDashPattern = 2
     
115. res3@gsnContourPosLineDashPattern = 0
     
116. plot(2)=gsn_csm_contour_map(wks,avel/10,res3)
     
117. 
     
118. ;;;;;;;;;;检验+差;;;;
     
119. opt=True
     
120. opt@gsnDraw=False
     
121. opt@gsnFrame = False
     
122. 
     
123. opt@gsnLeftString = ""
     
124. opt@gsnRightString = ""
     
125. opt@gsnMaximize = True
     
126. 
     
127. iopt=opt
     
128. opt@gsnAddCyclic = True
     
129. opt@mpFillOn= True
     
130. opt@mpCenterLonF = 180
     
131. opt@pmTickMarkDisplayMode = "Always"
     
132. opt@cnFillOn = False
     
133. opt@cnLinesOn = False
     
134. opt@cnLineLabelsOn=False
     
135. opt@cnInfoLabelOn = False
     
136. opt@cnLevelSelectionMode="ExplicitLevels"
     
137. opt@cnLevels = (/2.035,2.0351,2.0352,2.4,2.5,2.6,3/)
     
138. plot(3)=gsn_csm_contour_map(wks,t,opt) ;检验图
     
139. op= True
     
140. op@gsnShadeFillType = "pattern"
     
141. op@gsnShadeHigh = 17
     
142. plot(3)=gsn_contour_shade(plot(3),-999.,2.035,op) ;0.05
     
143. cha1=new((/91,180/),float)
     
144. do j=0,89
     
145. do i=0,90
     
146. cha1(i,j+90)=cha(i,j)
     
147. end do
     
148. end do
     
149. do j=90,179
     
150. do i=0,90
     
151. cha1(i,j-90)=cha(i,j)
     
152. end do
     
153. end do
     
154. cha1!0="lat"
     
155. cha1!1="lon"
     
156. cha1&lat=fspan(-90,90,91)
     
157. cha1&lon=fspan(0.,358.,180)
     
158. iopt@gsnAddCyclic = True
     
159. iopt@gsnShape = True
     
160. iopt@pmTickMarkDisplayMode = "Always"
     
161. iopt@cnLevelSelectionMode="ManualLevels"
     
162. iopt@cnMinLevelValF = -3
     
163. iopt@cnMaxLevelValF = 3
     
164. iopt@cnLevelSpacingF = 0.5
     
165. iopt@gsnContourNegLineDashPattern = 2
     
166. iopt@gsnContourPosLineDashPattern = 0
     
167. iopt@cnFillOn = False
     
168. iopt@cnLinesOn = True
     
169. iopt@cnLineLabelsOn=True
     
170. iopt@cnLineLabelBackgroundColor = "-1"
     
171. iopt@tmYLMode = "Explicit"
     
172. iopt@tmYLValues = (/-90,-45,0,45,90/)
     
173. iopt@tmYLLabels = (/"90~S~o~N~S","45~S~o~N~S","0","45~S~o~N~N","90~S~o~N~N"/)
     
174. iopt@tmXBMode = "Explicit"
     
175. iopt@tmXBValues = (/0,90,180,270/)
     
176. iopt@tmXBLabels = (/"0","90~S~o~N~E","180","90~S~o~N~W"/)
     
177. plo=gsn_csm_contour(wks,cha1/10,iopt)
     
178. wks6 = gsn_open_wks("pdf","hgtyy500_hc")
     
179. p=gsn_csm_contour(wks6,cha1/10,iopt)
     
180. draw(p)
     
181. frame(wks6)
     
182. overlay(plot(3),plo)
     
183. draw(plot(3))
     
184. frame(wks)
     
185. 
     
186. 
     
187. ;;;;;;;;;;;组图
     
188. resP= True
     
189. resP@gsnPanelRowSpec = True
     
190. resP@gsnPanelFigureStrings=(/"(a)","(b)","(c)","(d)"/)
     
191. resP@gsnPanelFigureStringsBackgroundFillColor = "-1"
     
192. resP@gsnPanelFigureStringsPerimOn = False
     
193. resP@gsnPanelFigureStringsFontHeightF = 0.015
     
194. resP@amJust = "BottomLeft"
     
195. resP@gsnPanelXWhiteSpacePercent = 0.2
     
196. resP@gsnPanelYWhiteSpacePercent = 0.2
     
197. ; resP@gsnPanelLabelBar = True
     
198. ; resP@pmLabelBarHeightF = 0.08
     
199. ; resP@gsnSpreadColorStart =1
     
200. gsn_panel(wks,plot,(/2,2/),resP)
     
201. 
     
202. 
     
203. 
     
204. 
     
205. end

复制代码

sususir

合成分析，为啥全球的数据只显示在0-90°N，0-180°E呢？