求助：合并的nc格式风场画图

突破的渴望

请问合并nc后的文件，不能出图，什么问题呢？

1. fn = 'D:/scripts/MeteoInfo/wnd2011.nc'
   
2. f = addfile(fn)
   
3. tidx = 102
   
4. u = f['uwnd'][tidx,[600],[70,120],[20,50]]
   
5. v = f['vwnd'][tidx,[600],[20,50],[70,120]]
   
6. speed = sqrt(u*u+v*v)
   
7. t = f.gettime(tidx)
   
8. #Plot
   
9. axesm()
   
10. lchina = shaperead('D:/scripts/MeteoInfo/map/china.shp')
    
11. geoshow(lchina, edgecolor='k')
    
12. layer = quiverm(u,v,speed,5,size=5)
    
13. colorbar(layer)
    
14. title(t.strftime('%Y-%m-%d %H:00'))

复制代码

MeteoInfo

仔细看看这里了解怎样从文件中读取数组数据：http://www.meteothink.org/docs/m ... /data_tutorial.html

突破的渴望

MeteoInfo 发表于 2019-9-25 08:26
仔细看看这里了解怎样从文件中读取数组数据：http://www.meteothink.org/docs/meteoinfolab/userguide/data ...

为什么只能在tidx = 0 的时候能出图？我想用别的时间，但是tidx变成别的数字就不能出图了，为什么？

1. fn = 'D:/scripts/MeteoInfo/wnd2011.nc'
   
2. f = addfile(fn)
   
3. tidx = 1
   
4. u = f['uwnd'][tidx,'600','15:55','70:140']
   
5. v = f['vwnd'][tidx,'600','15:55','70:140']
   
6. speed = sqrt(u*u+v*v)
   
7. t = f.gettime(tidx)
   
8. #Plot
   
9. axesm()
   
10. lchina = shaperead('D:/scripts/MeteoInfo/map/china.shp')
    
11. geoshow('country')
    
12. layer = quiverm(u, v, speed, 10, size=10)
    
13. colorbar(layer)
    
14. title(t.strftime('%Y-%m-%d %H:00'))

复制代码

MeteoInfo

突破的渴望 发表于 2019-9-25 10:32
为什么只能在tidx = 0 的时候能出图？我想用别的时间，但是tidx变成别的数字就不能出图了，为什么？

可能合并的有问题，请下载MeteoInfo最新版本合并数据然后再测试

突破的渴望

MeteoInfo 发表于 2019-9-25 10:39
可能合并的有问题，请下载MeteoInfo最新版本合并数据然后再测试

老师 我下载的是 Meteoinfo_2.0b2 了，这个单时间的还是出错，什么问题呢？还有就是如果是求0到10 这11个时间的平均，应该怎么改动呢？

MeteoInfo

突破的渴望 发表于 2019-9-25 10:41
老师 我下载的是 Meteoinfo_2.0b2 了，这个单时间的还是出错，什么问题呢？还有就是如果是求0到10 这11个 ...

应该是软件的bug，有空了修正一下

突破的渴望

MeteoInfo 发表于 2019-9-25 12:47
应该是软件的bug，有空了修正一下

老师，您解决了嘛？

MeteoInfo

突破的渴望 发表于 2019-9-26 11:31
老师，您解决了嘛？

解决了，这几天会推出一个新版本。

MeteoInfo

突破的渴望 发表于 2019-9-26 11:31
老师，您解决了嘛？

已经发布了MeteoInfo 2.0版本，你可以下载试试。

突破的渴望

老师，我用的新版本的MI2.0，但是合并nc文件之后发现 每个时间点画出来的图都是一样的 怎么回事？