气压层转高度

有事秦天

@MeteoInfo王老师您好，我用fnl画温度剖面图，想确定0℃层高度，该怎么才能气压层高度转换为米为单位的高度，并在图上显示0度层高度？
f = addfile('D:/MeteoInfo/1/fnl_20150727_12_00.grib2')
f = addfile(fn)
tidx = 0
level = '100000:20000'
#level1 = '100000:20000'
lat = '36:44'
lon = '113:120'
T = f['Temperature_isobaric'][tidx,level,lat,lon][::-1,:,:]
vort = T - 273.15
wd1 = vort[:,'40',:]
lev1 = wd1.dimvalue(0)/100
lev1 = lev1[::-1]
lev2 = meteo.p2h(lev1)
levels = []
for j in range(0, len(lev1)):
    levels.append('%i' % lev1[j])
wd1.setdimvalue(0, lev2)
print 'Plot...'
#layer = contourf(wd1,40)
levs = arange(-20,20,1)
layer = contourf(wd1, levs)
title(u'Vertical helicity (Longitude: 40°E)')
yticks(lev2, levels)
xlabel('Latitude')
ylabel('Pressure (hPa)')
colorbar(layer)

有事秦天

f = addfile('D:/MeteoInfo/1/fnl_20150727_18_00.grib2')
#f = addfile(fn)
tidx = 0
level = '100000:20000'
lat = '36:44'
lon = '113:120'
T = f['Temperature_isobaric'][tidx,level,lat,lon][::-1,:,:]
vort = T - 273.15
wd1 = vort[:,'40',:]
lev1 = wd1.dimvalue(0)/100
lev1 = lev1[::-1]
lev2 = meteo.pressure_to_height_std(lev1)
levels = []
for j in range(0, len(lev1)):
    levels.append('%i' % lev1[j])
wd1.setdimvalue(0, lev2)
print 'Plot...'
#layer = contourf(wd1,40)
axesm()
ax1=axes()

#levs = arange(-20,20,1)
#layer = contourf(wd1, levs)
levs = [-2,-1,0,1,2]
cols = [(93,71,139),(72,61,139),(0,191,255),(112,128,144),(143,188,143),(255,255,255)]
layer = contourf(wd1,levs,colors=cols)
colorbar(layer)

yaxis(tickin=False,tickfontsize=17)
xaxis(tickin=False,tickfontsize=17)
yaxis(location='right',tickin=True,tickfontsize=18)
#yticks(wd1.dimvalue(0), lev1)
xlim(113,120)
ylim(4000,5500)
ylabel('Hight(km)',fontsize=18)
xlabel('')
colorbar(layer,orientation='horizontal',aspect=35)

antialias(True)