用meteoinfo绘制流线图遇到的问题

zenghaode

请教用meteoinfo绘制流场问题，我下载了necp资料，名称为uwnd.2010.nc和vwnd.2010.nc，在绘图时出现提示画面如下，请教该如何解决，谢谢

MeteoInfo

MeteoInfo桌面软件目前只能绘制U/V变量在同一个文件中的情况，像这种U/V变量在两个不同文件中的情况可以用MeteoInfoLab脚本程序来绘图，之前的帖子答复过你。

zenghaode

老师您好，我按你给的贴子编辑脚本，但运行还是提示出错，我把脚本和运行结果给您看看，请您分析指点，谢谢
D:\wind\1.png

zenghaode

print 'open data files...'
f = ('D:/wind/model')
U1=('D:/wind/uwnd.2010.nc')
V1=('D:/wind/vwnd.2010.nc')
u = f['U1'][0,[500],[10,60],[60,140]]
v = f['V1'][0,[500],[10,60],[60,140]]
axesm()
mlayer = shaperead('D:/wind/map/country1.shp')
geoshow(mlayer, edgecolor='k')
layer = streamplotm(u, v)
title('Streamline')

下面是运行结果提示
open data files...
Traceback (most recent call last):
  File "<iostream>", line 5, in <module>
TypeError: str indices must be integers

MeteoInfo

zenghaode 发表于 2016-6-2 06:35
print 'open data files...'
f = ('D:/wind/model')
U1=('D:/wind/uwnd.2010.nc')

1. #Opend data files
   
2. f_uwnd = addfile('D:/Temp/nc/uwnd.2011.nc')
   
3. f_vwnd = addfile('D:/Temp/nc/vwnd.2011.nc')
   
4. uwnd = f_uwnd['uwnd'][0,0,:,:]
   
5. vwnd = f_vwnd['vwnd'][0,0,:,:]
   
6. #Plot
   
7. axesm()
   
8. mlayer = shaperead('D:/Temp/map/country1.shp')
   
9. geoshow(mlayer)
   
10. layer = streamplotm(uwnd, vwnd)
    
11. title('Streamline')
    
12. xlim(60, 140)
    
13. ylim(10, 60)

#Opend data files<br /> f_uwnd = addfile('D:/Temp/nc/uwnd.2011.nc')<br /> f_vwnd = addfile('D:/Temp/nc/vwnd.2011.nc')<br /> uwnd = f_uwnd['uwnd'][0,0,:,:]<br /> vwnd = f_vwnd['vwnd'][0,0,:,:]<br /> #Plot<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer)<br /> layer = streamplotm(uwnd, vwnd)<br /> title('Streamline')<br /> xlim(60, 140)<br /> ylim(10, 60)

zenghaode

我照搬了您给的脚本，只是更改了文件路径，但提示“addfile”没有定义，这是怎么回事？
Traceback (most recent call last):
  File "<iostream>", line 3, in <module>
NameError: name 'shaperead' is not defined
#Opend data files
f_uwnd = addfile('D:/Temp/nc/uwnd.2010.nc')
f_vwnd = addfile('D:/Temp/nc/vwnd.2010.nc')
uwnd = f_uwnd['uwnd'][0,0,:,:]
vwnd = f_vwnd['vwnd'][0,0,:,:]
#Plot
axesm()
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(mlayer)
layer = streamplotm(uwnd, vwnd)
title('Streamline')
xlim(60, 140)
ylim(10, 60)

zenghaode

是我安装的软件版本缺少什么文件吗？

MeteoInfo

zenghaode 发表于 2016-6-2 09:32
是我安装的软件版本缺少什么文件吗？

要在MeteoInfoLab里运行，在MeteoInfo目录中双击MeteoInfoLab.exe

zenghaode

终于绘出来了，谢谢王老师不厌其繁的指点，真的太感谢了

zenghaode

请问，这程序绘制的结果是某一天的还是一年的平均呢，谢谢您