LWC和IWC的计算方法

cibapotato

LWC和IWC的计算方法

LWC=Qcloud*1000*dens

IWC=(QICE+QSNOW) *1000*dens

其中QCLOUD为Cloud water mixing ratio (kg/kg)，QICE 为ice mixing ratio (kg/kg)，QSNOW为snow mixing ratio (kg/kg)，一般WRFOUT中都有这几个输出量；

1000是将kg/kg換成g/kg；

dens为密度，在一般WRFOUT中没有直接输出，需要我们进行以下几步：

1.找到Registry.EM_COMMON。

进入WRFV3下的Registry文件夹，打开Registry.EM_COMMON

操作命令为：vi Registry.EM_COMMON

2.搜索rho的位置

在vi中搜索只需要先键入/，如下图所示，键入/rho，再按N（代表next）键，就可以找到rho了。

3.将r改为rh
将rho前面的r修改为rh（至于用irh行不行我没有验证，但是rh一定可以），如果在vi下，需要先按下a启动修改，改完后按esc退出编辑，再在英文输入下输入:wq（保存并退出）

4.找到module_microphysics_driver.F

进入WRFV3下的phys文件夹，打开module_microphysics_driver.F

操作命令为：vi module_microphysics_driver.F

5.找到rho

在vi中搜索只需要先键入/，如下图所示，键入/rho，再按N（代表next）键，就可以找到rho了。找到rho定义的地方，将INTENT(IN)改为INTENT(INOUT).

6.重新编译WRF(非常重要)

退到WRFV3目录下，输入./clean，再输入./configure，再./compile em_real。这里注意，一定要clean，假如只重新./compile是不行的。

7.检查结果

运行完WRF后，

ncdump -h wrfout_d01_2018-01-28_00\:00\:00 >& wrfout_h.log

在wrfout_h.log搜索RHO（注意是大写的）,如果能搜索到RHO那就成功啦！！！！

cibapotato

IWC 和LWC分别是 ice water content和 liquid water content

xyan88

输出计算空气质量密度
感谢感谢

nightfalls

受教了！感谢楼主，还有一个问题。请问修改RHO是只要在module_microphysics_driver.F这一个程序里改就行了嘛？其他用到的mp_physics方案要改么？
比如我用Morroison方案，要在对应的module_mp_morr_two_moment.F程序里也把RHO改成inout吗？

森下下士

非常好帖子，爱来自中国