求助大神帮忙看看假相当位温的计算公式哪里有错误吗？

LIb

在用ERA5的再分析数据，想计算假相当位温，就从网上找了ncl的程序改了下，但是计算出来不敢确定对不对，想请大神帮忙看下，哪里有错误需要改正的，谢谢。
#计算假相当位温
data['e']     =data['level']*data['q']/(0.622+data['q'])
data['tlcl']  = 55.0+2840.0/(3.5* np.log(data['t'])- np.log(data['e'])-4.805)
data['theta'] = data['t'] * pow((1000/data['level']),(0.2854*(1.0-0.28*data['q'])))
data['eqt']   = data['theta'] * np.exp(((3376/data['tlcl'])-2.54)*data['q']*(1.0+0.81*data['q']))
-----######data['level']是等压面（单位hPa），data['q']是比湿（单位kg/kg），data['t']是温度（单位K）


这样算出来的data['eqt'] 是不是假相当位温呢？

wenbinz

ERA5有比湿这个要素吗？

LIb

wenbinz 发表于 2023-6-12 15:34
ERA5有比湿这个要素吗？

嗯，有比湿的数据。