使用windspharm计算速度势函数和辐合辐散风出现问题，看起来是边缘计算出错？

LemonTree007

# The standard interface requires that latitude and longitude be the leading
# dimensions of the input wind components, and that wind components must be
# either 2D or 3D arrays. The data read in is 3D and has latitude and
# longitude as the last dimensions. The bundled tools can make the process of
# re-shaping the data a lot easier to manage.
uwnd, uwnd_info = prep_data(U_data[:,::-1,:], 'tyx')
vwnd, vwnd_info = prep_data(V_data[:,::-1,:], 'tyx')
u200wnd, uwnd200_info = prep_data(U200_data[:,::-1,:], 'tyx')
v200wnd, vwnd200_info = prep_data(V200_data[:,::-1,:], 'tyx')
# It is also required that the latitude dimension is north-to-south. Again the
# bundled tools make this easy.
lat_ns=lats[::-1]
lats, uwnd, vwnd = order_latdim(lat_ns[8:45], uwnd, vwnd)
lats200, u200wnd, v200wnd = order_latdim(lat_ns[:45], u200wnd, v200wnd)
# Create a VectorWind instance to handle the computations.
w = VectorWind(uwnd, vwnd)
w200 = VectorWind(u200wnd, v200wnd)
# Compute components of rossby wave source: absolute vorticity, divergence,
# irrotational (divergent) wind components, gradients of absolute vorticity.
# eta = w.absolutevorticity()
# div = w.divergence()
uchi, vchi = w.irrotationalcomponent()
u200chi, v200chi = w200.irrotationalcomponent()
# Compute the streamfunction and velocity potential. Also use the bundled
# tools to re-shape the outputs to the 4D shape of the wind components as they
# were read off files.#画图只使用vp
sf, vp = w.sfvp()
sf = recover_data(sf, uwnd_info)
vp = recover_data(vp, uwnd_info)
uchi = recover_data(uchi,uwnd_info)
vchi = recover_data(vchi,vwnd_info)

sf200, vp200 = w200.sfvp()
sf200 = recover_data(sf200, uwnd200_info)
vp200 = recover_data(vp200, uwnd200_info)
u200chi = recover_data(u200chi,uwnd200_info)
v200chi = recover_data(v200chi,vwnd200_info)
#转换纬度为从南到北
vp=vp[:,::-1,:]
uchi=uchi[:,::-1,:]
vchi=vchi[:,::-1,:]

vp200=vp200[:,::-1,:]
u200chi=u200chi[:,::-1,:]
v200chi=v200chi[:,::-1,:]

uchi,lons_c = add_cyclic_point(uchi, lons[20:97])
vchi,lons_c = add_cyclic_point(vchi, lons[20:97])
vp,lons_c = add_cyclic_point(vp, lons[20:97])

u200chi,lons_c200 = add_cyclic_point(u200chi, lons[:97])
v200chi,lons_c200 = add_cyclic_point(v200chi, lons[:97])
vp200,lons_c200 = add_cyclic_point(vp200, lons[:97])

zjc130

流星大气 发表于 2021-9-1 12:45
你是在wsl下安装的吗

不是，直接是Windows

amandine

楼主，想问一下C:\Users\hp\Desktop、1.png你在什么通道成功安装windspharm库包的呀，我用anadonda安装一直显示安装失败

流星大气

amandine 发表于 2021-8-11 18:14
楼主，想问一下你在什么通道成功安装windspharm库包的呀，我用anadonda安装一直显示安装失败

windspharm需要在linux环境下安装

amandine

流星大气 发表于 2021-8-15 09:39
windspharm需要在linux环境下安装

好的，谢谢楼主回复。但是我在windows下也安装成功啦

zjc130

amandine 发表于 2021-8-29 22:17
好的，谢谢楼主回复。但是我在windows下也安装成功啦

windows怎么装的，我怎么安装失败

流星大气

amandine 发表于 2021-8-29 22:17
好的，谢谢楼主回复。但是我在windows下也安装成功啦

你是在wsl下安装的吗

谢君行

zjc130 发表于 2021-9-1 13:16
不是，直接是Windows

我用Windows一直安装不成功欸

zjc130

谢君行 发表于 2021-9-1 17:17
我用Windows一直安装不成功欸

我在wsl安装成功了， Windows上是安装不了的

流星大气

zjc130 发表于 2021-09-01 13:16
不是，直接是Windows

厉害

                               
登录/注册后可看大图

来自苹果APP客户端