MeteoInfoLab脚本示例：南方涛动指数Bar图

MeteoInfo

参考NCL的例子（http://www.ncl.ucar.edu/Applications/bar.shtml），用MeteoInfoLab脚本绘制了Darwin Southern Oscillation Index。数据文件里date变量是年和月的整数，比如199001，将此数据转换为日期列表作为x轴绘图。还有颜色的处理也很简单。

1. fn = 'D:/Temp/nc/soi.nc'
   
2. f = addfile(fn)
   
3. yms = f['date'][::8]       #Year and month
   
4. dsoik = f['DSOI_KET'][::8] #Darwin SOI Index via KET 11pt Smth
   
5. dsoid = f['DSOI_DEC'][::8] #Darwin Decadal SOI Index
   
6. 
   
7. #Set dates and colors
   
8. dates = []
   
9. cols = []
   
10. for ym,d in zip(yms,dsoik):
    
11.     dates.append(datetime.datetime(ym / 100, ym % 100, 1))
    
12.     if d >= 0:
    
13.         cols.append('r')
    
14.     else:
    
15.         cols.append('b')
    
16. 
    
17. #Bar plot
    
18. bar(dates, dsoik, color=cols, edgecolor=None)
    
19. xlim(dates[0], dates[-1])
    
20. ylim(-3, 3)
    
21. xaxis(axistype='time', minortick=True, tickin=False)
    
22. yaxis(minortick=True, tickin=False)
    
23. ylabel('Anomalias')
    
24. title('Darwin Southern Oscillation Index')

fn = 'D:/Temp/nc/soi.nc'<br /> f = addfile(fn)<br /> yms = f['date'][::8]       #Year and month<br /> dsoik = f['DSOI_KET'][::8] #Darwin SOI Index via KET 11pt Smth<br /> dsoid = f['DSOI_DEC'][::8] #Darwin Decadal SOI Index<br /> <br /> #Set dates and colors<br /> dates = []<br /> cols = []<br /> for ym,d in zip(yms,dsoik):<br />     dates.append(datetime.datetime(ym / 100, ym % 100, 1))<br />     if d >= 0:<br />         cols.append('r')<br />     else:<br />         cols.append('b')<br /> <br /> #Bar plot<br /> bar(dates, dsoik, color=cols, edgecolor=None)<br /> xlim(dates[0], dates[-1])<br /> ylim(-3, 3)<br /> xaxis(axistype='time', minortick=True, tickin=False)<br /> yaxis(minortick=True, tickin=False)<br /> ylabel('Anomalias')<br /> title('Darwin Southern Oscillation Index')

范式特西龙

请问 这个SOI.nc文件在哪里下载的呢？谢谢！

MeteoInfo

范式特西龙 发表于 2018-6-20 17:05
请问 这个SOI.nc文件在哪里下载的呢？谢谢！

可以在这里下载：http://www.ncl.ucar.edu/Applications/Data/

范式特西龙

MeteoInfo 发表于 2018-6-20 19:52
可以在这里下载：http://www.ncl.ucar.edu/Applications/Data/

非常感谢！

yangmt

请问，如果x轴想显示中文怎么设置呢？

MeteoInfo

yangmt 发表于 2024-1-31 07:10
请问，如果x轴想显示中文怎么设置呢？

参考这里：http://www.meteothink.org/docs/m ... ext_multi_axis.html