Hello!
I am using the function "wrf_user_getvar" in NCL to get the variable temperature ("tk"), but I noticed that it will call another .ncl file ($NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_inv_string.ncl) to calculate it. what I am confused is that in the following scrpits:
if( any( variable .eq. (/"tk","tc"/) ) ) then
;; function wrf_tk needs theta and pressure (Pa) on input and returns temperature in K on return
if(isfilevar(nc_file,"T")) then
T = _get_wrf_var(file_handle,"T",time)
P = _get_wrf_var(file_handle,"P",time)
PB = _get_wrf_var(file_handle,"PB",time)
T = T + 300.
P = P + PB
t = wrf_tk( P , T )
copy_VarAtts_except(T,t,"description"),
I don't know where this "300 K" comes from. I think 300 K should be the reference temperature. So is it equal to the variable "T00" in wrfout files? But I found the "T00" in my wrfout is 290 K, not 300 K.
Thank you in advance!
p.s. The wrf-chem version I used is v3.6.1.
I am using the function "wrf_user_getvar" in NCL to get the variable temperature ("tk"), but I noticed that it will call another .ncl file ($NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_inv_string.ncl) to calculate it. what I am confused is that in the following scrpits:
if( any( variable .eq. (/"tk","tc"/) ) ) then
;; function wrf_tk needs theta and pressure (Pa) on input and returns temperature in K on return
if(isfilevar(nc_file,"T")) then
T = _get_wrf_var(file_handle,"T",time)
P = _get_wrf_var(file_handle,"P",time)
PB = _get_wrf_var(file_handle,"PB",time)
T = T + 300.
P = P + PB
t = wrf_tk( P , T )
copy_VarAtts_except(T,t,"description"),
I don't know where this "300 K" comes from. I think 300 K should be the reference temperature. So is it equal to the variable "T00" in wrfout files? But I found the "T00" in my wrfout is 290 K, not 300 K.
Thank you in advance!
p.s. The wrf-chem version I used is v3.6.1.