- 积分
- 226
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2020-2-26
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我用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
我参考的公式是:
|
|