xaxis函数设置时间为x轴的时候，怎么规定时间间隔呢？

lhy277969385

ecthin_extend.CTL (592 Bytes, 下载次数: 2)

2020-2-22 16:04 上传

点击文件名下载附件
ecthin_physic.CTL (522 Bytes, 下载次数: 2)

2020-2-22 16:04 上传

点击文件名下载附件

ecthin_tlogp.CTL (407 Bytes, 下载次数: 3)

2020-2-22 16:04 上传

点击文件名下载附件
CTL文件如上所示，代码如下

1. clf()
   
2. figure(figsize=[900,750], newfig=False)
   
3. fn = 'D:/Meteoinfolab_Project/ctl/ecthin_extend.CTL'
   
4. fn2 = 'D:/Meteoinfolab_Project/ctl/ecthin_tlogp.CTL'
   
5. f = addfile(fn)
   
6. f2 = addfile(fn2)
   
7. 
   
8. #风数据U、V读取
   
9. u = f['u'][:,:,'39.9','119.5']
   
10. u = u.T
    
11. lev1 = u.dimvalue(0)
    
12. lev2 = meteo.pressure_to_height_std(lev1)
    
13. lev2 = lev2[:]
    
14. u.setdimvalue(0, lev2)
    
15. 
    
16. v = f['v'][:,:,'39.9','119.5']
    
17. v = v.T
    
18. lev3 = v.dimvalue(0)
    
19. lev4 = meteo.pressure_to_height_std(lev3)
    
20. lev4 = lev4[:]
    
21. v.setdimvalue(0, lev4)
    
22. 
    
23. speed = sqrt(u*u+v*v)
    
24. 
    
25. #RH数据读取
    
26. rh = f2['rhprs'][:,:,'39.9','119.5']
    
27. rh = rh.T
    
28. lev5 = rh.dimvalue(0)
    
29. lev6 = meteo.pressure_to_height_std(lev5)
    
30. lev6 = lev6[:]
    
31. rh.setdimvalue(0, lev6)
    
32. 
    
33. #Omega数据读取
    
34. omega = f['vvel'][:,:,'39.9','119.5']
    
35. omega = omega.T
    
36. omega.setdimvalue(0, lev8)
    
37. 
    
38. #温度数据读取
    
39. temp = f2['tmpprs'][:,:,'39.9','119.5']
    
40. temp = temp.T
    
41. temp.setdimvalue(0, lev10)
    
42. 
    
43. #风及RH绘制
    
44. levs2 = arange(60,100,10)
    
45. cols = [(255,255,255),(0,255,255),(153,255,153),(0,255,0),(50,205,50)]
    
46. layer2 = contourf(rh,levs2,colors=cols)
    
47. layer1 = barbs(u,v,speed,12,color='k')
    
48. 
    
49. #Omega绘制
    
50. levs = arange(-12, 12, 3)
    
51. bss = []
    
52. for i in range(len(levs)):
    
53.     bs={}
    
54.     bs['value'] = levs[i]
    
55.     bs['color'] = 'k'
    
56.     if levs[i] < 0:
    
57.         bs['linestyle'] = '--'
    
58.     bss.append(bs)
    
59. ls = plt.makesymbolspec('line',bss)
    
60. layer3 = contour(omega,symbolspec=ls)
    
61. clabel(layer3)
    
62. 
    
63. #温度绘制
    
64. levs2 = arange(-80, 40, 4)
    
65. bss2 = []
    
66. for i in range(len(levs2)):
    
67.     bs2={}
    
68.     bs2['value'] = levs2[i]
    
69.     bs2['color'] = 'r'
    
70.     if levs2[i] < 0:
    
71.         bs2['linestyle'] = '--'
    
72.     bss2.append(bs2)
    
73. ls2 = plt.makesymbolspec('line',bss2)
    
74. layer4 = contour(temp,symbolspec=ls2)
    
75. clabel(layer4)
    
76. #colorbar(layer)
    
77. colorbar(layer2, extendrect=False)
    
78. 
    
79. #设置XY轴
    
80. xaxis(axistype='time',timetickformat='ddHH',minortick=True,minorticknum=5)
    
81. #xticks(arange(0, 240, 10))
    
82. yticks(lev6,lev5)
    
83. xlabel('Time')
    
84. ylabel('Pressure(hPa)')
    
85. title('54449 Time Series of RH(%),WIND(m/s),Omega(hPa/s),Temp,' + '\n' + 'ECthin,00-240 Hours Forecast',fontname=u'黑体', fontsize=24, bold=False, color='blue')

复制代码

作图如下，其实图片没问题，主要是下边的时间坐标轴，不知道该怎么给它固定成比如12h或24h间隔

求王老师解答~！@MeteoInfo  

MeteoInfo

用xticks，值可以是datetime列表。

lhy277969385

收到，谢谢王老师！！