Dear WRF users,
I am pretty new to WRF and would very much appreciate your help. I need to calculate temperature at 1000 m above ground level from WRF output. I am using the NCL function "wrf_user_getvar". This function can extract data from WRF files and calculate some diagnostics. One of the available diagnostics is "z/height" which is labeled as the "full model height in m" in this NCL site (https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml).
I am doing the following:
a = addfile("wrfout_d01_2000-01-24_12:00:00.nc","r")
ht=wrf_user_getvar(a,"z",0) ; full model height (m)
tk=wrf_user_getvar(a,"tk",0); Temperature in K
temp_1km=wrf_user_interp_level(tk,ht,1000,0) ; interpolate temperature at 1 km
But I am confused if "wrf_user_getvar" measures z from ground level or sea level. I need temperature at 1000 m above ground level, not sea level.
Can anybody please tell me if my approach is correct?
Best,
Israt
I am pretty new to WRF and would very much appreciate your help. I need to calculate temperature at 1000 m above ground level from WRF output. I am using the NCL function "wrf_user_getvar". This function can extract data from WRF files and calculate some diagnostics. One of the available diagnostics is "z/height" which is labeled as the "full model height in m" in this NCL site (https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml).
I am doing the following:
a = addfile("wrfout_d01_2000-01-24_12:00:00.nc","r")
ht=wrf_user_getvar(a,"z",0) ; full model height (m)
tk=wrf_user_getvar(a,"tk",0); Temperature in K
temp_1km=wrf_user_interp_level(tk,ht,1000,0) ; interpolate temperature at 1 km
But I am confused if "wrf_user_getvar" measures z from ground level or sea level. I need temperature at 1000 m above ground level, not sea level.
Can anybody please tell me if my approach is correct?
Best,
Israt