爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6315|回复: 1

[脚本编辑] metpy里面divergence计算散度时因有时间维度使用循环但出错?

[复制链接]

新浪微博达人勋

发表于 2020-10-22 18:10:27 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
dx,dy=mpcalc.lat_lon_grid_deltas(lon[20:85],lat[28:65])
div1= np.zeros((39,37,65))
div2= np.zeros((39,37,65))
div3= np.zeros((39,37,65))
for i in range(39):
    div1[i,:,:]=mpcalc.divergence(U200_summer[i,:,:],V200_summer[i,:,:],dx,dy,dim_order='yx')
    div2[i,:,:]=mpcalc.divergence(U500_summer[i,:,:],V500_summer[i,:,:],dx,dy,dim_order='yx')
    div3[i,:,:]=mpcalc.divergence(U850_summer[i,:,:],V850_summer[i,:,:],dx,dy,dim_order='yx')
出错信息:div1[i,:,:]=mpcalc.divergence(U200_summer[i,:,:],V200_summer[i,:,:],dx,dy,dim_order='yx')TypeError: __array__() takes 1 positional argument but 2 were given
感觉没毛病啊?如果不考虑时间维度不作循环是没问题的,为啥加了时间就出问题了,等式前后的维度是一致的呀?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2020-10-22 18:14:26 | 显示全部楼层
比如我这个是没问题的
div=mpcalc.divergence(U200_summer[0,:,:],V200_summer[0,:,:],dx,dy,dim_order='yx')
得到的div也是(37,65)大小的shape,,但加了时间为什么循环就会出错呢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表