wrfout输出变量,求相对湿度RH的计算公式

acking

现在得到的wrfout输出变量里面没有相对湿度RH，但是绘制很多图都需要RH。
目前得到的变量有高层温度和地面海平面温度SSTSK，基态温度T00，扰动位温T，地面2m高度的温度T2，地面2m高度的比湿Q2，求计算相对湿度的公式，在网上找了很久都没找到  

lee309

NCL中调用相对湿度就是rh，直接用就可以了

acking

lee309 发表于 2014-12-11 12:25
NCL中调用相对湿度就是rh，直接用就可以了

是因为没有rh，所有才求计算公式的。我这里得到的数据文件名是wrfout_d01_2014-12-10_12:00:00.GFS_WCTRL_F。在里面没找到RH

acking

好吧，最后从网上看到这个文件WRFUserARW.ncl
发现里面有公式，直接翻译过来就可以用了


if( variable .eq. "rh" ) then
       if(isfilevar(nc_file,"T")) then
         ;; first compute theta - function wrf_tk needs theta and pressure (Pa) on input
         ;; THEN compute rh, using tk, p (Pa), QVAPOR
         if ( time .eq. -1 ) then
           if(ISFILE) then
             T      = nc_file->T
             P      = nc_file->P
             PB     = nc_file->PB
             QVAPOR = nc_file->QVAPOR
           else
             T      = file_handle[:]->T
             P      = file_handle[:]->P
             PB     = file_handle[:]->PB
             QVAPOR = file_handle[:]->QVAPOR
           end if
         else
           if(ISFILE) then
             T      = nc_file->T(time_in,:,:,:)
             P      = nc_file->P(time_in,:,:,:)
             PB     = nc_file->PB(time_in,:,:,:)
             QVAPOR = nc_file->QVAPOR(time_in,:,:,:)
           else
             T      = file_handle[:]->T(time_in,:,:,:)
             P      = file_handle[:]->P(time_in,:,:,:)
             PB     = file_handle[:]->PB(time_in,:,:,:)
             QVAPOR = file_handle[:]->QVAPOR(time_in,:,:,:)
           end if
         end if  
         T = T + 300.
         P  = P + PB
         QVAPOR = QVAPOR > 0.000
         tk = wrf_tk( P , T )
         rh = wrf_rh( QVAPOR, P, tk )
         delete_VarAtts(T,(/"description","units"/))
         copy_VarAtts(T,rh)
       else
         ;; may be a met_em file - see if we can get RH
         if(isfilevar(nc_file,"RH")) then
           if ( time .eq. -1 ) then
             if(ISFILE) then
               rh = nc_file->RH
             else
               rh = file_handle[:]->RH
             end if
           else
             if(ISFILE) then
               rh = nc_file->RH(time_in,:,:,:)
             else
               rh = file_handle[:]->RH(time_in,:,:,:)
             end if
           end if
         end if
       end if  
       return(rh)
  end if

{:eb513:}{:eb513:}

清梅浅影

好长的代码啊

虾米go

正好也在用ncl画rh，谢谢啦

毛毛虫106

学习中

Amyhiganbana

{:lxm_23:}{:lxm_23:}{:lxm_23:}

nheweijun

认真看这页， 搜索 RH 就有计算相对湿度的，http://www.ncl.ucar.edu/Document/Functions/wrf.shtml

kwb

学习了，谢谢