请教一下各位大佬关于时空剖面图的问题。

Tobias

y轴坐标是对的，但是风场数据是反的，应该怎么调整啊
fn=r'c:/work/code/20220430/wrfout_d03_2022-04-29_12_00_00.ctl'f=addfile(fn)
u=f['U'][:,::-1,'43.91','87.47'].T
v=f['V'][:,::-1,'43.91','87.47'].T

lev1=u.dimvalue(0)
lev1=lev1[::-1]
lev2=meteo.pressure_to_height_std(lev1)
lev2=lev2[:]/1000
u.setdimvalue(0,lev2)
v.setdimvalue(0,lev2)
barbs(u[:,::2],v[:,::2])
xaxis(axistype='time',timetickformat='ddHH')
levy=array([1000,950,900,850,800,750,700,650,600,550,500,450,400,350,300,250,200])
yticks(u.dimvalue(0),levy)

Tobias

图片已附上

MeteoInfo

Tobias 发表于 2022-5-17 11:48
图片已附上

取u, v数据的时候不要对高度维取反

Tobias

MeteoInfo 发表于 2022-5-18 14:20
取u, v数据的时候不要对高度维取反

fn=r'c:/work/code/20220430/wrfout_d03_2022-04-29_12_00_00.ctl'
f=addfile(fn)
u=f['U'][:,:,'43.91','87.47'].T
v=f['V'][:,:,'43.91','87.47'].T

lev1=u.dimvalue(0)
lev1=lev1[::-1]
lev2=meteo.pressure_to_height_std(lev1)
lev2=lev2[:]/1000
u.setdimvalue(0,lev2)
v.setdimvalue(0,lev2)
barbs(u[:,::2],v[:,::2])
xaxis(axistype='time',timetickformat='ddHH')
levy=array([1000,950,900,850,800,750,700,650,600,550,500,450,400,350,300,250,200])
yticks(u.dimvalue(0),levy)

王老师，这样子还是不行啊