MeteoInfoLab脚本示例：风场矢量图

MeteoInfo

半颗麦芽糖 发表于 2016-11-22 18:50
王老师，看了您的这个帖子，再结合该帖http://bbs.06climate.com/forum.php?mod=viewthread&tid=45172&extr ...

eta坐标差值为等压面目前还没有现成的函数。

如果风场矢量太密的话，可以先将U, V分量数组抽稀，然后再绘图。

半颗麦芽糖

MeteoInfo 发表于 2016-11-22 23:36
eta坐标差值为等压面目前还没有现成的函数。

如果风场矢量太密的话，可以先将U, V分量数组抽稀，然后 ...

多谢王老师指点！

半颗麦芽糖

MeteoInfo 发表于 2016-11-22 23:36
eta坐标差值为等压面目前还没有现成的函数。

如果风场矢量太密的话，可以先将U, V分量数组抽稀，然后 ...

王老师，我利用nc格式资料通过脚本和界面画了图，但是不知道为什么画出的图有差异（图中气旋所在），画图要求是2014年10月26日12时500hPa风矢量

附上数据与脚本
uv.nc (4.95 MB, 下载次数: 45)

2016-11-30 16:22 上传

点击文件名下载附件

1. 
   
2. <P> fn = 'E:/TBB/20141026/uv.nc'
   
3. f = addfile(fn)
   
4. u = f['U'][103,[500],[20,50],[70,120]]
   
5. v = f['V'][103,[500],[20,50],[70,120]]
   
6. speed = sqrt(u*u+v*v)
   
7. #Plot
   
8. axesm()
   
9. lchina = shaperead('F:/11/MeteoInfo/map/bou2_4p.shp')
   
10. ltibet = shaperead('F:/11/MeteoInfo/map/DBATP/tibet_ASCII.wmp')
    
11. geoshow(lchina, edgecolor='k')
    
12. geoshow(ltibet, edgecolor='r', size=2)
    
13. layer = quiverm(u,v,speed,5,size=5)
    
14. colorbar(layer)
    
15. </P>

复制代码

MeteoInfo

半颗麦芽糖 发表于 2016-11-30 16:23
王老师，我利用nc格式资料通过脚本和界面画了图，但是不知道为什么画出的图有差异（图中气旋所在），画 ...

很有可能是一个重大的bug，我会尽快修正。

MeteoInfo

半颗麦芽糖 发表于 2016-11-30 16:23
王老师，我利用nc格式资料通过脚本和界面画了图，但是不知道为什么画出的图有差异（图中气旋所在），画 ...

的确是程序里的一个bug（MeteoInfoLab的结果是错的），已经修正，请下载最新版本（MeteoInfo 1.3.8R1）。

MeteoInfo

半颗麦芽糖 发表于 2016-11-30 16:23
王老师，我利用nc格式资料通过脚本和界面画了图，但是不知道为什么画出的图有差异（图中气旋所在），画 ...

另外，2016-10-26 12:00是第102个时次（Python中数组序号是从0开始的）。通过文件对象的gettime()函数能够获取某个时次的具体时间变量。

1. fn = 'E:/Temp/uv.nc'
   
2. f = addfile(fn)
   
3. tidx = 102
   
4. u = f['U'][tidx,[500],[20,50],[70,120]]
   
5. v = f['V'][tidx,[500],[20,50],[70,120]]
   
6. speed = sqrt(u*u+v*v)
   
7. t = f.gettime(tidx)
   
8. #Plot
   
9. axesm()
   
10. lchina = shaperead('D:/Temp/map/bou2_4p.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'))

fn = 'E:/Temp/uv.nc'<br /> f = addfile(fn)<br /> tidx = 102<br /> u = f['U'][tidx,[500],[20,50],[70,120]]<br /> v = f['V'][tidx,[500],[20,50],[70,120]]<br /> speed = sqrt(u*u+v*v)<br /> t = f.gettime(tidx)<br /> #Plot<br /> axesm()<br /> lchina = shaperead('D:/Temp/map/bou2_4p.shp')<br /> geoshow(lchina, edgecolor='k')<br /> layer = quiverm(u,v,speed,5,size=5)<br /> colorbar(layer)<br /> title(t.strftime('%Y-%m-%d %H:00'))

新版本也可以用时间变量（datetime.datetime）来提取数组。

1. fn = 'E:/Temp/uv.nc'
   
2. f = addfile(fn)
   
3. t = datetime.datetime(2014,10,26,12)
   
4. u = f['U'][[t],[500],[20,50],[70,120]]
   
5. v = f['V'][[t],[500],[20,50],[70,120]]
   
6. speed = sqrt(u*u+v*v)
   
7. #Plot
   
8. axesm()
   
9. lchina = shaperead('D:/Temp/map/bou2_4p.shp')
   
10. geoshow(lchina, edgecolor='k')
    
11. layer = quiverm(u,v,speed,5,size=5)
    
12. colorbar(layer)
    
13. title(t.strftime('%Y-%m-%d %H:00'))

fn = 'E:/Temp/uv.nc'<br /> f = addfile(fn)<br /> t = datetime.datetime(2014,10,26,12)<br /> u = f['U'][[t],[500],[20,50],[70,120]]<br /> v = f['V'][[t],[500],[20,50],[70,120]]<br /> speed = sqrt(u*u+v*v)<br /> #Plot<br /> axesm()<br /> lchina = shaperead('D:/Temp/map/bou2_4p.shp')<br /> geoshow(lchina, edgecolor='k')<br /> layer = quiverm(u,v,speed,5,size=5)<br /> colorbar(layer)<br /> title(t.strftime('%Y-%m-%d %H:00'))

半颗麦芽糖

MeteoInfo 发表于 2016-11-30 19:17
另外，2016-10-26 12:00是第102个时次（Python中数组序号是从0开始的）。通过文件对象的gettime()函数能 ...

王老师太给力了，膜拜！
但是新版本我下了两次，压缩包解压的时候都显示文件损坏，不知道是什么原因，两个电脑都试了。

MeteoInfo

半颗麦芽糖 发表于 2016-11-30 20:24
王老师太给力了，膜拜！
但是新版本我下了两次，压缩包解压的时候都显示文件损坏，不知道是什 ...

有可能，今天网络不给力，上传的时候断过，我再重新传一下。

半颗麦芽糖

MeteoInfo 发表于 2016-11-30 21:52
有可能，今天网络不给力，上传的时候断过，我再重新传一下。

好的，等老师传完了我再下载试试

MeteoInfo

半颗麦芽糖 发表于 2016-11-30 21:55
好的，等老师传完了我再下载试试

你可以在百度网盘上下载：http://pan.baidu.com/s/1kUENi4n