wrf_user_getvar这个函数在读取变量时能get到很多头文件里没有但是常用的变量，官网函数说明里也有
记录一些自己常用常忘的变量
1. z = wrf_user_getvar(a_1, "z",-1) ; grid point height
2. pt = wrf_user_getvar(a_1, "th",-1) ; potential temperature (theta) 位温
th/theta是位温，tc是温度（摄氏度），tk是温度（开尔文）
3. float SWDOWN(Time, south_north, west_east) ;地表高度的向下的短波辐射通量
SWDOWN:description = "DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE" ;
SWDOWN:units = "W m-2" ;
4. float GLW(Time, south_north, west_east) ;地表高度的向下的长波辐射通量
GLW:description = "DOWNWARD LONG WAVE FLUX AT GROUND SURFACE" ;
GLW:units = "W m-2" ;
5. float OLR(Time, south_north, west_east) ;大气层顶向上的长波辐射
OLR:description = "TOA OUTGOING LONG WAVE" ;
OLR:units = "W m-2" ;
WRF输出变量说明_weixin_30417487的博客-CSDN博客
最近在计算辐射相关变量的时候还想输出大气加热率，找到了一个相关的回答，但是这里面的变量都需要在Registry里面进行修改，重新编译后再运行模式才能输出
Output the atmospheric heating rate | WRF & MPAS-A Support Forum (ucar.edu)
The diabatic heating terms include the following:
"h_diabatic": this is from microphsics scheme
"RTHCUTEN": this is from cumulus scheme
"RTHSHTEN": this is from shallow convection
"RTHRATEN": this is from radaition
"RTHBLTEN": this is from boundary layer
The sum of the above terms is the total diabatic heating.
RTHFTEN is the total advective temperature tendency, including both horizontal and vertical advection terms.
下面还有长短波分开引起的加热率，都是需要修改Registry的
