关于微物理方案中预报变量初始化问题

無№敌→臉

以wrf为例
输入模拟数据时，应该是只有比湿q和云水qc（我查的fnl里貌似湿变量只有q和云水qc）
在调用微物理方案的时候需要好几个预报变量，比如wsm3方案需要水汽q，云水／云冰，雨／雪
wsm5方案需要5个预报变量 水汽q，云水qc，云冰qi，雨qr，雪qs
那么除了能从fnl资料里得到的水汽和云水，其他几个量如何得到？
我在wsm5微物理方案计算过程中，对雨qr进行追踪，发现方案计算降雨是用qr，然而我对该方案移植后，对qr初始化未搞清楚，所以求助如何在用输入资料去计算得到这些多出的变量。我在对qr初始化为0的话，由于计算过程全部于qr有关，后面积分输出全都为0
所以现在求大神帮忙解答下关于微物理方案所用到的几个变量如何在wrf程序中初始化
最近看wrf代码头大也没找到在哪读入数据，初始化也没找到
就找到了grhd%moist（：，：，：，：）这个是给那几个湿变量存储值的地方，前三个是i j k纬度，第四个代表是水汽，云水，云冰，雨，雪等等
目前关于其他几个量如何得到没有进展