gsn_csm_xy怎么改X轴的坐标

浅水叶

gsn_csm_xy(wks,time2,data1, res)，time2明明用的是20151025：20151120的时间坐标，怎么变这样

DTyanyan

是不是可以缩小X坐标轴的范围，限定在数据的20151025：20151120范围应该就好了

1649518749

脚本呢？{:5_275:}{:5_275:}

liweier123

请问楼主有没有解决，我也遇到了这个问题

liweier123

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/contrib/time_axis_labels.ncl"
   
5. 
   
6. ;----------------------------------------------------------------------
   
7. ; Converts YYYYMM time to "days since" time.
   
8. ;----------------------------------------------------------------------
   
9. function yyyymm_to_time(time)
   
10. local ntime, stime, year, month, day, hour, minute, second, units
    
11. begin
    
12.   ntime  = dimsizes(time)
    
13.   stime  = tostring(time)
    
14.   year   = toint(str_get_cols(stime,0,3))
    
15.   month  = toint(str_get_cols(stime,4,5))
    
16.   day    = new(ntime,integer)
    
17.   hour   = new(ntime,integer)
    
18.   minute = new(ntime,integer)
    
19.   second = new(ntime,double)
    
20. 
    
21. ;----------------------------------------------------------------------
    
22. ; Note: you may need to change the day, hour, minute values depending
    
23. ; on your definition of "YYYYMM".  Is this the first day and the
    
24. ; first hour of the month? Is it the middle of the month?
    
25. ;----------------------------------------------------------------------
    
26.   day    = 1
    
27.   hour   = 0
    
28.   minute = 0
    
29.   second = 0.d
    
30.   units  = "days since " + min(year) + "-01-01 00:00"
    
31. 
    
32.   newtime = cd_inv_calendar(year,month,day,hour,minute,second,units,0)
    
33. 
    
34.   return(newtime)
    
35. end
    
36. ;----------------------------------------------------------------------
    
37. ; Main code
    
38. ;----------------------------------------------------------------------
    
39. begin
    
40.   shpname = "/cygdrive/g/Tibet/add_TP/DBATP_Line.shp"
    
41.   ;afiles  =  systemfunc("ls /cygdrive/g/Tibet/cloudsat/cloud-height/height*.nc")
    
42.   files  =  systemfunc("ls  /cygdrive/g/Tibet/cloudsat/data/height_monthly*.nc")
    
43.       f  =  addfiles(files, "r")
    
44. ListSetType(f, "cat")
    
45. ;printVarSummary(f)
    
46.     data = f[:]->valavg_monthly
    
47. ;printVarSummary(data)
    
48. dims = dimsizes(data)
    
49. ntim = dims(0)
    
50. time = data&time
    
51. ;time!0     = "time"
    
52. ;time@units = "YYYYMM"
    
53. ;print(time)
    
54. ; convert integer YYYYMM to float
    
55. ;timeF = yyyymm_to_yyyyfrac(time,0)
    
56. ;data&time := timeF
    
57.   latS  = 25.0
    
58.   latN  = 40.0
    
59.   lonL  = 70.0
    
60.   lonR  = 105.0
    
61. ;========================================  余弦加权
    
62. lat     = data&lat
    
63. clat    = sqrt(cos(lat*0.01745329))
    
64. CLAT    = conform(data, clat, 1)
    
65. dataw   = data*CLAT
    
66. copy_VarCoords(data, dataw)
    
67. ;printVarSummary(dataw)
    
68. ;======================================== 对加权数据进行EOFs分解并显著性检验
    
69. neof    = 3
    
70. optEOF  = True
    
71. optEOF@jopt = 0                              ;默认值，计算使用协方差矩阵
    
72. optETS  = False
    
73. eof     = eofunc_n_Wrap(dataw, neof, optEOF, 0)
    
74. eof     = smth9_Wrap(eof, 0.5, 0.5, True)
    
75. eof_ts  = eofunc_ts_n_Wrap(dataw, eof, optETS, 0)
    
76. printVarSummary(eof)                         ; examine EOF variables
    
77. printVarSummary(eof_ts)
    
78. eof_ts  = dim_standardize_n_Wrap(eof_ts, 0, 1);标准化时间序列 0为样本标准差即N-1 1为进行标准化的维度
    
79. pcvar   = eof@pcvar                           ;获取方差贡献属性项
    
80. prinfo  = True
    
81. sig_ev  = eofunc_north(eof@eval, ntim, prinfo)  ;差异显著性
    
82. ;======================================== 定义时间序列eof_ts的属性
    
83. eof_ts!0    = "eofs"
    
84. ;eof_ts!1    = "time"
    
85. eof_ts&eofs = (/1,2,3/)                       ;三个模态
    
86. ;eof_ts&time = year
    
87. ;print(eof_ts(0,:))
    
88. ;============================================================================
    
89. wks     = gsn_open_wks("png", "EOF_MaxdBZ_height_month_pinghua-07-10")
    
90. stringw  = (/"(a)","(b)","(c)","(d)","(e)","(f)"/)
    
91. ploteof = new ( 6, "graphic")
    
92. ;=================================================;
    
93. ;空间场相关设置
    
94. ;=================================================;
    
95.    res                      = True                ; contour plots
    
96.    res@gsnDraw              = False
    
97.    res@gsnFrame             = False
    
98.    res@gsnLeftString        = ""
    
99.    res@gsnRightString       = ""
    
100.    res@cnFillOn             = True                ; turn on color
     
101.    res@cnFillMode           = "RasterFill"        ; turn on raster mode
     
102.    res@cnLinesOn            = False               ; turn off contour lines
     
103.    res@cnLineLabelsOn       = False
     
104.    res@cnInfoLabelOn        = False
     
105.    res@cnFillPalette        = "BlWhRe"   ; specify colormap
     
106.    res@mpShapeMode          = "FreeAspect"
     
107.    res@vpHeightF            = 0.3
     
108.    res@vpWidthF             = 0.6
     
109.    res@gsnAddCyclic         = False                ; force cyclic value
     
110.    res@lbLabelBarOn         = False               ; turn off individual lb's
     
111.    res@mpFillOn             = False
     
112.    res@mpMinLatF            = latS                 ; range to zoom in on
     
113.    res@mpMaxLatF            = latN
     
114.    res@mpMinLonF            = lonL  
     
115.    res@mpMaxLonF            = lonR
     
116.    res@tmXBLabelFontHeightF = 0.02
     
117.    res@tmYLLabelFontHeightF = 0.02
     
118.    res@gsnLeftStringFontHeightF  = 0.02
     
119.    res@gsnRightStringFontHeightF = 0.02
     
120.    symMinMaxPlt(eof, 40, False, res)              ; 设置几种颜色
     
121. 
     
122.    do n=0,2
     
123.       res@gsnLeftString     = stringw(n)+"EOF"+(n+1)
     
124.       res@gsnRightString    = sprintf("%5.1f", pcvar(n)) +"%"
     
125.       ploteof(n) = gsn_csm_contour_map(wks,eof(n,:,:),res)
     
126.    end do
     
127.       tpres                  =True
     
128.       tpres@gsLineColor      = "black"     ;设置廓线颜色
     
129.       tpres@gsLineThicknessF = 3              ;设置廓线宽度
     
130.       poly0=gsn_add_shapefile_polylines(wks,ploteof(0),shpname,tpres)        ;plot为底图
     
131.       poly1=gsn_add_shapefile_polylines(wks,ploteof(1),shpname,tpres)
     
132.       poly2=gsn_add_shapefile_polylines(wks,ploteof(2),shpname,tpres)
     
133. 
     
134. 
     
135. ;===========================================
     
136. ; 时间场相关设置
     
137. ;===========================================
     
138. ;---Convert the time values and reassign as x's time coordinate array
     
139. ;  newtime = yyyymm_to_time(time)
     
140. ;  eof_ts&time  := newtime
     
141. ;  printVarSummary(eof_ts)
     
142. 
     
143.    rts           = True
     
144.    rts@gsnDraw   = False       ; don't draw yet
     
145.    rts@gsnFrame  = False       ; don't advance frame yet
     
146. ;   rts@gsnScale  = True
     
147.    ;rts@tmXBFormat= "f"         ; no unnessary decimal pts
     
148. 
     
149.    rts@vpHeightF = 0.40        ; Changes the aspect ratio
     
150.    rts@vpWidthF  = 0.6
     
151.    ;rts@vpXF      = 0.20        ; change start locations
     
152.    ;rts@vpYF      = 0.80        ; the plot
     
153. 
     
154.    rts@gsnYRefLine           = 0.              ; reference line
     
155. ;   rts@gsnXYBarChart         = True            ; create bar chart   
     
156.    rts@gsnAboveYRefLineColor = "red"           ; above ref line fill red
     
157.    rts@gsnBelowYRefLineColor = "blue"          ; below ref line fill blue
     
158. 
     
159.    rts@trXMinF   = min(time)   
     
160.    rts@trXMaxF   = max(time)      ; nicer end bound
     
161.    rts@tiXAxisString = "Time"
     
162.    rts@tiYAxisString = ""
     
163. ;   rts@tmXBMode      = "Explicit"
     
164. ;   rts@tmXBValues    = time(::12)
     
165. ;   rts@tmXBLabels    = (/"Jan ~C~2007","Jan ~C~2008","Jan ~C~2009","Jan ~C~2010"/)
     
166.    rts@tmXBLabelFontHeightF = 0.02
     
167.    rts@tmYLLabelFontHeightF = 0.02
     
168.    res@tmLabelAutoStride    = True                    ; nice stride on labels
     
169.    rts@gsnLeftStringFontHeightF = 0.02
     
170. ;---Set resources necessary to nicely format X axis
     
171. ;  restick           = True
     
172. ;  restick@ttmFormat = "%c%Y"              ; Jan 2010, Feb 2010, etc.
     
173.   ;time_axis_labels(newtime,rts,restick)
     
174. 
     
175.    do n=0,2
     
176.       rts@gsnLeftString  = stringw(n+3)+"PC"+(n+1)
     
177.       ploteof(n+3) = gsn_csm_xy(wks,time,eof_ts(n,:),rts)
     
178.    end do
     
179. 
     
180.    resP                     = True                ; panel plots
     
181.    resP@gsnMaximize         = True                ; large format
     
182.    resP@gsnPanelLabelBar    = True                ; add common colorbar
     
183.    resP@gsnPaperOrientation = "portrait"          ; force portrait [landscape]
     
184.    resP@amJust              = "TopLeft"
     
185.    gsn_panel(wks,ploteof,(/2,3/),resP)        ; draw all 'neof' as one plot
     
186. 
     
187.    end
     
188. 
     

复制代码

图只是eof时间序列，子程序试过也画不出来