ncl用ERA资料算假相当位温

lyh0205

我用ERA的850hpa相对湿度、比湿、温度三个要素的资料计算假相当位温，读取某一日数据，计算后出现很多错误值,我的代码是这样的：

begin
f0 = addfile("E:/ncdata/85008r.nc","r")
f1 = addfile("E:/ncdata/85008t.nc","r")
f2 = addfile("E:/ncdata/85008q.nc","r")
data1 = f0->r
data2 = f1->t
data3 = f2->q
rhum = data1(4,{70:-10},{50:160})
t = data2(4,:,:)
q = data3(4,:,:)*1000
e = 850*q/(0.622+q)
tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)
theta = t*((1000/850)^(0.2854*(1.0-0.28*q)))
ose = theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))
end



我参考的公式是：

                               
登录/注册后可看大图

lc2016

首先，era的里面存储的变量数据是short字型的，data1 = f0->r 这样写得到的数据是偏移的，你不信你可以print（data1），我一般这样写，data1=short2flt（f0->r）,这样的data1数据就是真实值了。

woshiwzm

你先看看你的单位弄对没。然后试试把你读取的每个变量都取相同的经纬度，tlcl，theta,ose都记得要命名维度。我之前画水汽通量散度不命名维度数据就也是错的

袁春梅

楼主画出来了吗？可以参考一下吗？

woshiwzm

我单层画出来没问题，但是画垂直的总感觉是错的

pengtingting

楼主你的q*1000是啥意思，不乘1000好像是对的

nightelfxu

q不要乘以1000就对了