NCL交叉谱

rpf

1. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
3. 
   
4. ymdStrt = 20080101
   
5. ymdLast = 20101231
   
6. 
   
7. begin
   
8. 
   
9. filename = "/cygdrive/C/Users/Administrator/Desktop/Write_RMM.nc/MJO_PC_INDEX_BSISO_97-14.nc"
   
10. f = addfile(filename, "r")
    
11. time = f->time
    
12. YMD = cd_calendar(time, -2)
    
13. iStrt = ind(YMD .eq. ymdStrt)
    
14. iLast = ind(YMD .eq. ymdLast)
    
15. delete([/time,YMD/])
    
16. PC1 = f->PC1(iStrt:iLast)
    
17. PC2 = f->PC2(iStrt:iLast)
    
18. time = f->time(iStrt:iLast)
    
19. YMD = cd_calendar(time, 0)
    
20. ;print(YMD)
    
21. month = YMD(:,1)
    
22. ;print(month)
    
23. nt = ind(month .le. 10 .and. month .ge. 5)
    
24. ;print(nt)
    
25. ;printVarSummary(PC1)
    
26. pc1 = PC1(time|nt)
    
27. pc2 = PC2(time|nt)
    
28. ;printVarSummary(pc1)
    
29. delete([/nt,month,YMD,time,iStrt,iLast/])
    
30. iopt = 0 ; remove least squares linear trends from each
    
31. ; series prior to tapering and computing spectra.
    
32. jave = 7 ; Average 7 periodogram estimates using modified Daniell
    
33. pct = 0.1 ; taper 10% of the data
    
34. dof = specxy_anal (pc1,pc2,iopt,jave,pct)
    
35. x = dof@frq
    
36. y1 = dof@coher
    
37. y2 = dof@phase
    
38. nt = ind(y2 .lt. 0.)
    
39. y2(nt) = y2(nt)+360.
    
40. 
    
41. wks = gsn_open_wks("png", "cs")
    
42. res1 = True
    
43. res1@gsnFrame = False
    
44. res1@trXLog = True
    
45. res1@vpOn = True
    
46. res1@vpXF = 0.1
    
47. res1@vpYF = 0.65
    
48. res1@vpWidthF = 0.8
    
49. res1@vpHeightF = 0.5
    
50. res1@xyComputeXMin = True
    
51. res1@tmXTOn = False
    
52. 
    
53. plot1 = gsn_csm_xy(wks, x, y1, res1)
    
54. 
    
55. polyres = True
    
56. polyres@gsLineDashPattern = 1
    
57. 
    
58. gsn_polyline(wks, plot1, (/1/10.,1/10./), (/0,1./), polyres)
    
59. gsn_polyline(wks, plot1, (/1/20.,1/20./), (/0,1./), polyres)
    
60. gsn_polyline(wks, plot1, (/1/30.,1/30./), (/0,1./), polyres)
    
61. gsn_polyline(wks, plot1, (/1/80.,1/80./), (/0,1./), polyres)
    
62. 
    
63. 
    
64. res1@vpYF = 0.9
    
65. res1@vpHeightF = 0.25
    
66. res1@tmXBOn = False
    
67. res1@tmXTOn = True
    
68. res1@tmXTLabelsOn = True
    
69. res1@tmYLOn = False
    
70. res1@tmYROn = True
    
71. res1@tmYRLabelsOn = True
    
72. 
    
73. res1@tmYLMode = "Explicit"
    
74. res1@tmYLValues = (/0,90,180,270,360/)
    
75. res1@tmYLLabels = (/"0","90","180","270","360"/)
    
76. res1@trYMaxF = 360.
    
77. 
    
78. res1@tmXBMode = "Explicit"
    
79. res1@tmXBValues = (/1/10.,1/20.,1/30.,1/80./)
    
80. res1@tmXBLabels = (/"10d","20d","30d","80d"/)
    
81. 
    
82. res1@trYMaxF = 360.
    
83. plot2 = gsn_csm_xy(wks, x, y2, res1)
    
84. 
    
85. gsn_polyline(wks, plot2, (/1/10.,1/10./), (/0,360./), polyres)
    
86. gsn_polyline(wks, plot2, (/1/20.,1/20./), (/0,360./), polyres)
    
87. gsn_polyline(wks, plot2, (/1/30.,1/30./), (/0,360./), polyres)
    
88. gsn_polyline(wks, plot2, (/1/80.,1/80./), (/0,360./), polyres)
    
89. 
    
90. frame(wks)
    
91. 
    
92. end

复制代码

今天尝试着用NCL绘制两时间序列的交叉谱图，遇到了一些比较奇怪的问题，请大神帮助解答附件图中上下两幅图中的虚线为什么不能对应起来？

rpf

问题已自己解决，可惜了1现金{:shy:}