gsn_panel在站点打点图总应用出错

sususue

各位大神，我想要将站点EOF的2空间模态和时间模态画在一起，研究了一下官网，自己改了一下，但是提示gsn_panel: warning: you have more plots than you have panels.(0)    Only 1 plots will be drawn.但是实在找不到自己出错在哪儿了，麻烦大家帮我看看了！万分感谢！

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. load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
   
5. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
   
6. 
   
7. begin
   
8. path = "/disk02/hg-ws/data40G/quanguo/test/kongzhi/eof/"  
   
9. f = addfile(path+"rasu3d.nc", "r")
   
10. 
    
11. rasu = f->rasu3d
    
12. lat = f->lat
    
13. lon = f->lon
    
14. 
    
15. 
    
16. ndim         =  1178
    
17. neof         = 2
    
18. optEOF       =  True
    
19. optEOF@jopt  =  0  
    
20. optETS       = True
    
21. optETS@jopt  =0
    
22. 
    
23. ; =================================================================
    
24. ; create weights:  sqrt(cos(lat))   [or sqrt(gw) ] for covariance
    
25. ; ; =================================================================
    
26.    rad    = 4.0*atan(1.0)/180.0
    
27.    clat   = cos(rad*lat)           
    
28.    clat   = where(clat.lt.0, 0.0, clat)  ; avoid a potential numerical issue at pole
    
29.    clat   = sqrt( clat )  ; avoid a potential numerical issue at pole
    
30. ; =================================================================
    
31. ; weight all observations
    
32. ; =================================================================
    
33.                                    ; copy meta data
    
34.   wSH   = rasu*conform(rasu, clat, 1)
    
35.   copy_VarCoords(rasu,wSH)
    
36. ;  wSH@long_name = "Wgt: "+wSH@long_name
    
37. ; =================================================================
    
38. ; Reorder (lat,lon,time) the *weighted* input data
    
39. ; Access the area of interest via coordinate subscripting
    
40. ; =================================================================
    
41.   x      = wSH(lat|:,lon|:,time|:)
    
42.   x@_FillValue = 999999
    
43.   eof    = eofunc_Wrap(x, neof, optEOF)      
    
44.   eof_ts = eofunc_ts_Wrap(x, eof, optETS)
    
45. 
    
46.   printVarSummary( eof )                         ; examine EOF variables
    
47.   printVarSummary( eof_ts )
    
48. 
    
49.     eoff=new((/ndim ,neof /),"float")
    
50. do j=0,neof -1
    
51.   do i=0,ndim -1  
    
52.     eoff(i,j) = eof(j,i,i)
    
53.   end do
    
54. end do
    
55. ; =================================================================
    
56. ; Normalize time series: Sum spatial weights over the area of used
    
57. ; =================================================================
    
58.   dimx   = dimsizes( x )
    
59.   mln    = dimx(1)
    
60.   sumWgt = mln*sum( clat(:) )
    
61.   eof_ts = eof_ts/sumWgt
    
62. ;print(sum( clat(:) ))
    
63. ;exit
    
64.   
    
65. ;*******************************************
    
66. ;North significance test: any pcvar, eval_transpose, eval can be used
    
67. ;*******************************************
    
68.     ;sig  = new((/40,60/),"logical")
    
69.     dimp   = dimsizes( rasu )
    
70.     ntim   = dimp(0)  
    
71.     ;print(ntim)
    
72. ;                                           ; max # eigenvalues possible
    
73.     prinfo = True
    
74.     sig_ev  = eofunc_north(eof@eval, ntim, prinfo)
    
75.     sig_pcv = eofunc_north(eof@pcvar, ntim, prinfo)
    
76.     sig_evt = eofunc_north(eof@eval_transpose, ntim, prinfo)  
    
77. ;   printVarSummary(sig)
    
78. 
    
79. ;----------------------------------------------------------------------
    
80. ; This function attaches a labelbar to the given plot.
    
81. ;----------------------------------------------------------------------
    
82. function attach_labelbar(wks,map,arr[*]:numeric,colors[*])
    
83. local lbres, vph, vpw, nboxes
    
84. begin
    
85.   getvalues map
    
86.     "vpHeightF" : vph
    
87.     "vpWidthF"  : vpw
    
88.   end getvalues
    
89. 
    
90.   nboxes = dimsizes(colors)
    
91.    
    
92.   lbres                    = True          ; labelbar only resources
    
93.   lbres@lbAutoManage       = False          ; Necessary to control sizes
    
94.   lbres@lbFillColors       = colors
    
95.   lbres@vpWidthF           = 0.7 * vpw     ; labelbar width
    
96.   lbres@vpHeightF          = 0.2 * vph     ; labelbar height
    
97.   lbres@lbMonoFillPattern  = True          ; Solid fill pattern
    
98.   lbres@lbLabelFontHeightF = 0.01          ; font height. default is small
    
99.   lbres@lbOrientation      = "horizontal"
    
100.   lbres@lbPerimOn          = False
     
101.   lbres@lbLabelAlignment   = "InteriorEdges"
     
102.   lbid = gsn_create_labelbar(wks,nboxes,""+arr,lbres)
     
103. 
     
104.   amres                  = True
     
105.   amres@amJust           = "TopCenter"
     
106.   amres@amParallelPosF   =  0.0    ; Center
     
107.   amres@amOrthogonalPosF =  0.6    ; Bottom
     
108.   annoid = gsn_add_annotation(map,lbid,amres)
     
109.   return(annoid)
     
110. end
     
111. ;============================================================
     
112. ; PLOTS
     
113. ;============================================================
     
114.   
     
115.   wks = gsn_open_wks("pdf","eofrasu")
     
116.   gsn_define_colormap(wks,"rainbow+gray")     ; choose colormap
     
117.   plot = new(neof,graphic)  
     
118.   dum1 =new(neof,graphic)
     
119.   map=new(neof,graphic)
     
120. ;=================Set Map=================;
     
121. res                         = True            
     
122. res@tiMainString            = " "
     
123. res@gsnMaximize             = True
     
124. res@gsnDraw                 = False
     
125. res@gsnFrame                = False
     
126. 
     
127.   res@mpMinLatF               = 17.                        
     
128.   res@mpMaxLatF               = 55.
     
129.   res@mpMinLonF               = 72.
     
130.   res@mpMaxLonF               = 136.
     
131. 
     
132. res@mpFillOn                = True
     
133. ; res@mpLandFillColor         =  100
     
134. res@mpOutlineOn             = False  ; Use outlines from shapefile
     
135. ;res@cnFillDrawOrder         = "PreDraw"
     
136. res@mpDataBaseVersion       = "MediumRes"
     
137. res@mpDataSetName           = "Earth..4"
     
138. res@mpAreaMaskingOn         = True
     
139. res@mpMaskAreaSpecifiers    = (/"China","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)
     
140. res@mpLandFillColor         = "white"
     
141. res@mpInlandWaterFillColor  = "white"
     
142. res@mpOceanFillColor        = "white"
     
143. res@mpOutlineBoundarySets   = "NoBoundaries"
     
144. res@vpXF            = 0.05   ;離屏幕左端的距離      
     
145. res@vpWidthF        = 0.4    ;圖的寬度
     
146. res@vpYF            = 0.2    ;離屏幕頂端的距離
     
147. res@vpHeightF       = 0.3    ;圖的高度
     
148. 
     
149.   res@mpGeophysicalLineThicknessF= 1.      ; double the thickness of geophysical boundaries
     
150.   res@mpNationalLineThicknessF= 1.
     
151. 
     
152. 
     
153.   res@pmTickMarkDisplayMode   = "Always"
     
154.   ; panel plot only resources
     
155.   resP                     = True         ; modify the panel plot
     
156.   resP@gsnMaximize         = True         ; large format
     
157.   resP@gsnPanelLabelBar    = True         ; add common colorbar
     
158.   ; resP@lbLabelAutoStride   = True         ; auto stride on labels
     
159.   resP@gsnPanelYWhiteSpacePercent = 2
     
160.   resP@gsnPanelXWhiteSpacePercent = 2
     
161. 
     
162.   ;*******************************************
     
163. ; second plot
     
164. ;*******************************************
     
165. ; EOF time series  [bar form]
     
166. 
     
167.   rts           = True
     
168.   rts@gsnDraw   = False       ; don't draw yet
     
169.   rts@gsnFrame  = False       ; don't advance frame yet
     
170.   rts@gsnScale  = True        ; force text scaling               
     
171. 
     
172.   rts@vpHeightF = 0.40        ; Changes the aspect ratio
     
173.   rts@vpWidthF  = 0.85
     
174.   rts@vpXF      = 0.10        ; change start locations
     
175.   rts@vpYF      = 0.75        ; the plot
     
176.    rts@trXMinF  = 1967
     
177.    rts@trXMaxF  = 2014
     
178. 
     
179.   rts@tiYAxisString = ""                    ; y-axis label      
     
180.     resz                         = True
     
181.     resz@xyLineColors            = "red"
     
182. 
     
183. ; panel plot only resources
     
184.   rtsP                      = True            ; modify the panel plot
     
185.   rtsP@gsnMaximize          = True            ; large format
     
186.   rtsP@txString             = "EOF Time Series for rasu"
     
187. 
     
188.   year = ispan (1967,2014,1)
     
189. ;****************************************************
     
190. 
     
191.   do n=0,1
     
192.   res@gsnLeftString  = "EOF "+(n+1)
     
193.   res@gsnRightString = sprintf("%5.1f", eof@pcvar(n)) +"%"
     
194.   ;res@tmXBLabelFontHeightF  = 0.005
     
195.   ;res@tmYBLabelFontHeightF  = 0.005
     
196.   res@gsnLeftStringFontHeightF = 0.01
     
197.   res@gsnRightStringFontHeightF = 0.01
     
198.   map(n)= gsn_csm_map(wks,res)   
     
199. 
     
200. 
     
201.      ;=================Draw Map=================;
     
202. 
     
203. cnres           = True
     
204. cnres@china     = True        ; draw china map or not
     
205. cnres@river     = True         ; draw changjiang&huanghe or not
     
206. cnres@province  = True   ; draw province boundary or not
     
207. cnres@nanhai    = False     ; draw nanhai or not
     
208. cnres@diqu      = True
     
209. chinamap = add_china_map(wks,map(n),cnres)
     
210. ;***************************
     
211. ; Plot rotated patterns
     
212. ;*******************************************
     
213.   npts = ndim
     
214.   arr = (/-0.12,-0.08,-0.04,0,0.04,0.08,0.12/)
     
215.                
     
216.   narr = dimsizes(arr)
     
217.   labels = new(narr+1,string)  ; Labels for legend.
     
218.   RR   = eoff(:,n)  ; This is dummy data for determining
     
219.    
     
220. ;
     
221.   num_distinct_markers = dimsizes(arr)+1        ; number of distinct markers
     
222.   lat_new = new((/num_distinct_markers,dimsizes(RR)/),float,-99999)
     
223.   lon_new = new((/num_distinct_markers,dimsizes(RR)/),float,-99999)
     
224.   
     
225. ;---Group the points according to which range they fall in.
     
226.    do i = 0, num_distinct_markers-1
     
227.      if (i.eq.0) then
     
228.        indexes = ind(RR.lt.arr(0))
     
229.      end if
     
230.      if (i.eq.num_distinct_markers-1) then
     
231.        indexes = ind(RR.ge.max(arr))
     
232.      end if
     
233.      if (i.gt.0.and.i.lt.num_distinct_markers-1) then      
     
234.        indexes = ind(RR.ge.arr(i-1).and.RR.lt.arr(i))
     
235.      end if
     
236.    
     
237.      ;index into range
     
238.      if (.not.any(ismissing(indexes))) then
     
239.        npts_range = dimsizes(indexes)   ; # of points in this range.
     
240.       lat_new(i,0:npts_range-1) = lat(indexes)
     
241.       lon_new(i,0:npts_range-1) = lon(indexes)
     
242.      end if
     
243.      delete(indexes)            ; Necessary b/c "indexes" may be a different
     
244.                                 ; size next time.
     
245.    end do
     
246. 
     
247. ; Begin plotting section.
     
248. ;==================================
     
249.   gsres               = True
     
250.   gsres@gsMarkerIndex = 16  
     
251.   gsres@gsMarkerSizeF = 6. ; Increase marker size by a factor of 10.   
     
252. ;---Get nice spacing through color map for marker colors
     
253. getvalues wks
     
254.     "wkColorMapLen" : clen     ; number of colors in color map
     
255.   end getvalues
     
256. 
     
257.   nstep = (clen-2)/narr
     
258.   colors = ispan(3,clen-1,nstep)
     
259. ;---Loop through each "bin" and attach the markers to the map.
     
260.      ;draw eof
     
261.     do i = 0, num_distinct_markers-1
     
262.       if (.not.ismissing(lat_new(i,0)))
     
263.       gsres@gsMarkerColor      = colors(i)
     
264.       dumstr = unique_string("marker")
     
265.       dum1(n) = gsn_add_polymarker(wks,map(n),lon_new(i,:),lat_new(i,:),gsres)
     
266.       end if
     
267.     end do
     
268.        ;draw time series
     
269.       rts@gsnLeftString  = "EOF "+(n+1)
     
270.       rts@gsnRightString = sprintf("%5.1f", eof@pcvar(n )) +"%"
     
271.       plot(n) = gsn_csm_xy (wks,year,eof_ts(n,:),rts)
     
272.       draw(map(n))
     
273.       frame(wks)
     
274.     delete([/lat_new,lon_new/])
     
275. 
     
276. end do
     
277.   res_U                   = True
     
278.   res_U@gsnFrame          = False
     
279.   res_U@gsnPanelTop       = 0.90
     
280.   res_U@gsnPanelBottom    = 0.45
     
281.   res_U@gsnPanelLabelBar  = True
     
282.   res_U@pmLabelBarWidthF  = 0.6
     
283.   res_U@lbBoxEndCapStyle  = "TriangleBothEnds"     ; default is "RectangleEnds"
     
284.   
     
285.   res_L                   = True
     
286.   res_L@gsnFrame          = False
     
287.   res_L@gsnPanelTop       = 0.45
     
288.   res_L@gsnPanelLabelBar  = True
     
289.   res_L@pmLabelBarWidthF  = 0.6
     
290.   res_L@lbBoxEndCapStyle  = "TriangleHighEnd"     ; Low end will have rectangle end
     
291. 
     
292.   txres               = True
     
293.   txres@txFontHeightF = 0.03
     
294.   gsn_text_ndc(wks,"lbBoxEndCapStyle: triangle-end labelbars",0.5,0.95,txres)
     
295.   
     
296.   gsn_panel(wks,dum1,(/1,1/),res_U)
     
297.   gsn_panel(wks,plot,(/1,1/),res_L)
     
298.   
     
299.   frame(wks)
     
300. 
     
301. ;psres = True                                                               
     
302. ;maximize_output(wks,psres)  ; calls draw and frame for you      
     
303. 
     
304. 
     
305. end
     
306.   

复制代码

qinqin

297.  gsn_panel(wks,plot,(/1,1/),res_L)改成
gsn_panel(wks,plot,(/1,2,5,3,1/),res_L)这种
(/1,2,5,3,1/)表明每行几张图
今天在这卡了好久，知道对你没用了，但是还是回一下

qinqin

突然发现我是因为 resP@gsnPanelRowSpec  ; tell panel what order to plot
不小心设成True   这个设成true的话gsn_panel(wks,plot,(/1,2,5,3,1/),res_L)，一共出12张图
否则gsn_panel(wks,plot,(/3,4/),resP)，(/3,4/)就表明图几行几列。

你的脚本没细看但感觉可能是gsn_panel设置的问题。。。。？