NCL 位温计算问题

半边晴天3

http://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_rh.shtml     
a = addfile("wrfout_d01_2000-01-24_12:00:00.nc","r")  
qv = a->QVAPOR  
P  = a->P                     ; perturbation  
Pb = a->PB                    ; base state pressure  
P  = P + Pb                   ; total pressure
theta = a->T                  ; perturbation potential temperature (theta+t0)
theta = theta + 300.  
TK = wrf_tk (P, theta)  R
H = wrf_rh (qv, P, TK)  
printVarSummary(RH)

我想问theta = theta + 300. 这必须得加300么，如果我的wrfout*文件中
float T00(Time=25);  :FieldType = 104; // int  :MemoryOrder = "0  ";  :description = "BASE STATE TEMPERATURE";  :units = "K";T00是等于290的，是否应该变成theta = theta + 290

曼联现在没$

为什么不直接读取pressure和tk直接作为气压和温度？

半边晴天3

曼联现在没$ 发表于 2016-12-19 11:25
为什么不直接读取pressure和tk直接作为气压和温度？

wrfout*的输出数据一般只有位温，扰动气压和基本气压，并没有直接的气压和温度啊

曼联现在没$

用ncl的wrf_user_getvar可以读取气压和温度

半边晴天3

曼联现在没$ 发表于 2016-12-19 18:13
用ncl的wrf_user_getvar可以读取气压和温度

真的可以，多谢指点

半边晴天3

半边晴天3 发表于 2016-12-19 19:39
真的可以，多谢指点

这个函数是不是只能得到某个时次的温度？可是我想得到所有时次的，构成一个四维数组进行运算

半边晴天3

半边晴天3 发表于 2016-12-19 19:42
这个函数是不是只能得到某个时次的温度？可是我想得到所有时次的，构成一个四维数组进行运算

我知道了可以设置为-1，读取所有时次的，谢谢啦