WRFV3.6.1版本的Morrison方案切换为冰雹选项出错

小其其格

使用WRFV3.6.1版本，三层双向嵌套模拟一个冰雹的case，

d01和d02选用WSM6方案6，d03选Morrison方案10

1、按照https://blog.csdn.net/weixin_34000916/article/details/88148901网址教程快速添加
wrfout输出变量的方法添加了变量，放在WRFV3.6.1/run/的下面的三个iofield_list.txt里：
+:h:0:QHAIL

2、namelist.input的&time_control下面添加：
iofields_filename                   = "iofield_list_d01.txt", "iofield_list_d02.txt", "iofield_list_d03.txt",

3、namelist.input的&physics 下面打开冰雹的开关：
hail_opt                 = 1,


最后转出来d03的wrfout里，绘制图像时，提示找不到QHAIL变量，


大神帮忙看看是怎么回事，感激不尽！！！


附上namelist.input

1. &time_control
   
2. run_days                            = 0, &time_control
   
3. run_days                            = 0,
   
4. run_hours                           = 24
   
5. run_minutes                         = 0,
   
6. run_seconds                         = 0,
   
7. start_year                          = 2019, 2019, 2019,
   
8. start_month                         = 03, 03, 03,
   
9. start_day                           = 20, 20, 20,
   
10. start_hour                          = 12,   12,   18,   
    
11. start_minute                        = 00,   00,   00,
    
12. start_second                        = 00,   00,   00,   
    
13. end_year                            = 2019, 2019, 2019,
    
14. end_month                           = 03, 03, 03,
    
15. end_day                             = 21, 21, 21,
    
16. end_hour                            = 12,   12,   06,
    
17. end_minute                          = 00,   00,   00,   
    
18. end_second                          = 00,   00,   00,   
    
19. interval_seconds                    = 21600
    
20. input_from_file                     = .true.,   .true.,   .true.,
    
21. history_interval                    = 60,  60, 10,  
    
22. frames_per_outfile                  = 1, 1, 1,
    
23. restart                             = .false.,
    
24. restart_interval                    = 5000,
    
25. io_form_history                     = 2
    
26. io_form_restart                     = 2
    
27. io_form_input                       = 2
    
28. io_form_boundary                    = 2
    
29. debug_level                         = 0
    
30. !########### add new output variables to wrfout
    
31. iofields_filename                   = "iofield_list_d01.txt", "iofield_list_d02.txt", "iofield_list_d03.txt",  !#在WRFV3/run下面，里面增加了qh(QHAIL)变量的输出
    
32. !#ignore_iofields_warning             = .true.,
    
33. !#auxhist0_interval                  = 60,  60, 10,  
    
34. !#frames_per_auxhist0                = 1, 1, 1,
    
35. !#io_form_auxhist0                   = 2 !#netCDF format
    
36. /
    
37. 
    
38. &domains                 
    
39. time_step                = 180,
    
40. time_step_fract_num      = 0,
    
41. time_step_fract_den      = 1,
    
42. max_dom                  = 3,
    
43. e_we                     = 130,      214,      322,
    
44. e_sn                     = 112,      196,      313,
    
45. e_vert                   = 35,       35,       35,
    
46. p_top_requested          = 5000,
    
47. num_metgrid_levels       = 32,
    
48. num_metgrid_soil_levels  = 4,
    
49. dx                       = 27000,     9000,     3000,
    
50. dy                       = 27000,     9000,     3000,
    
51. grid_id                  = 1,        2,        3,
    
52. parent_id                = 1,        1,        2,
    
53. i_parent_start           = 1,       30,       80,
    
54. j_parent_start           = 1,       24,       25,
    
55. parent_grid_ratio        = 1,        3,        3,
    
56. parent_time_step_ratio   = 1,        3,        3,
    
57. feedback                 = 1,
    
58. smooth_option            = 0,
    
59. /
    
60. 
    
61. &physics                 
    
62. mp_physics               = 6,        6,       10, !#微物理过程(mp_physics)
    
63. ra_lw_physics            = 1,        1,        1, !#长波辐射（ra_lw_physics）
    
64. ra_sw_physics            = 1,        1,        1, !#短波辐射（ra_sw_physics）
    
65. radt                     = 30,       30,       30,!#30,       30,       30,
    
66. sf_sfclay_physics        = 1,        1,        1, !#近地面层 (sf_sfclay_physics)
    
67. sf_surface_physics       = 2,        2,        2, !#陆面层（sf_surface_physics）
    
68. bl_pbl_physics           = 1,        1,        1, !#行星边界层（bl_pbl_physics）
    
69. bldt                     = 0,        0,        0,
    
70. cu_physics               = 1,        1,        0, !#积云参数化（cu_physics）
    
71. cudt                     = 5,        5,        5,
    
72. isfflx                   = 1,
    
73. ifsnow                   = 0,
    
74. icloud                   = 1,
    
75. surface_input_source     = 1,
    
76. num_soil_layers          = 4, !#指定陆面模式中的土壤层数，默认值为 5   #5         热量扩散方案   #4 Noah 陆面过程方案   #6  RUC 陆面过程方案
    
77. sf_urban_physics         = 0,        0,        0,
    
78. maxiens                  = 1,
    
79. maxens                   = 3,
    
80. maxens2                  = 3,
    
81. maxens3                  = 16,
    
82. ensdim                   = 144,
    
83. num_land_cat             = 20,      !number of land categories in input data.
    
84.                                                !   24 - for USGS (default); 20 for MODIS
    
85.                                                !   28 - for USGS if including lake category
    
86.                                                !   21 - for MODIS if including lake category
    
87. hail_opt                 = 1,!#morrison方案中afwa_hail_opt=1为输出膘默认0输出霰. WRFV3.6.1后变为hail_opt
    
88. /
    
89. 
    
90. &fdda                    
    
91. /
    
92. 
    
93. &dynamics               
    
94. w_damping                = 0,
    
95. diff_opt                 = 1,
    
96. km_opt                   = 4,
    
97. diff_6th_opt             = 0,        0,        0,
    
98. diff_6th_factor          = 0.12,     0.12,     0.12,
    
99. base_temp                = 290.,
    
100. damp_opt                 = 0,
     
101. zdamp                    = 5000.,    5000.,    5000.,
     
102. dampcoef                 = 0.2,      0.2,      0.2,
     
103. khdif                    = 0,        0,        0,
     
104. kvdif                    = 0,        0,        0,
     
105. non_hydrostatic          = .true.,   .true.,   .true.,
     
106. moist_adv_opt            = 1,        1,        1,
     
107. scalar_adv_opt           = 1,        1,        1,
     
108. /
     
109. 
     
110. &bdy_control            
     
111. spec_bdy_width           = 5,
     
112. spec_zone                = 1,
     
113. relax_zone               = 4,
     
114. specified                = .true.,  .false.,  .false.,
     
115. nested                   = .false.,   .true.,   .true.,
     
116. /
     
117. 
     
118. &grib2                  
     
119. /
     
120. 
     
121. &namelist_quilt         
     
122. nio_tasks_per_group      = 0,
     
123. nio_groups               = 1,
     
124. /                                                                                                                        
     
125. 
     
126. run_hours                           = 24
     
127. run_minutes                         = 0,
     
128. run_seconds                         = 0,
     
129. start_year                          = 2019, 2019, 2019,
     
130. start_month                         = 03, 03, 03,
     
131. start_day                           = 20, 20, 20,
     
132. start_hour                          = 12,   12,   18,   
     
133. start_minute                        = 00,   00,   00,
     
134. start_second                        = 00,   00,   00,   
     
135. end_year                            = 2019, 2019, 2019,
     
136. end_month                           = 03, 03, 03,
     
137. end_day                             = 21, 21, 21,
     
138. end_hour                            = 12,   12,   06,
     
139. end_minute                          = 00,   00,   00,   
     
140. end_second                          = 00,   00,   00,   
     
141. interval_seconds                    = 21600
     
142. input_from_file                     = .true.,   .true.,   .true.,
     
143. history_interval                    = 60,  60, 10,  
     
144. frames_per_outfile                  = 1, 1, 1,
     
145. restart                             = .false.,
     
146. restart_interval                    = 5000,
     
147. io_form_history                     = 2
     
148. io_form_restart                     = 2
     
149. io_form_input                       = 2
     
150. io_form_boundary                    = 2
     
151. debug_level                         = 0
     
152. !########### add new output variables to wrfout
     
153. iofields_filename                   = "iofield_list_d01.txt", "iofield_list_d02.txt", "iofield_list_d03.txt",  !#在WRFV3/run下面，里面增加了qh(QHAIL)变量的输出
     
154. !#ignore_iofields_warning             = .true.,
     
155. !#auxhist0_interval                  = 60,  60, 10,  
     
156. !#frames_per_auxhist0                = 1, 1, 1,
     
157. !#io_form_auxhist0                   = 2 !#netCDF format
     
158. /
     
159. 
     
160. &domains                 
     
161. time_step                = 180,
     
162. time_step_fract_num      = 0,
     
163. time_step_fract_den      = 1,
     
164. max_dom                  = 3,
     
165. e_we                     = 130,      214,      322,
     
166. e_sn                     = 112,      196,      313,
     
167. e_vert                   = 35,       35,       35,
     
168. p_top_requested          = 5000,
     
169. num_metgrid_levels       = 32,
     
170. num_metgrid_soil_levels  = 4,
     
171. dx                       = 27000,     9000,     3000,
     
172. dy                       = 27000,     9000,     3000,
     
173. grid_id                  = 1,        2,        3,
     
174. parent_id                = 1,        1,        2,
     
175. i_parent_start           = 1,       30,       80,
     
176. j_parent_start           = 1,       24,       25,
     
177. parent_grid_ratio        = 1,        3,        3,
     
178. parent_time_step_ratio   = 1,        3,        3,
     
179. feedback                 = 1,
     
180. smooth_option            = 0,
     
181. /
     
182. 
     
183. &physics                 
     
184. mp_physics               = 6,        6,       10, !#微物理过程(mp_physics)
     
185. ra_lw_physics            = 1,        1,        1, !#长波辐射（ra_lw_physics）
     
186. ra_sw_physics            = 1,        1,        1, !#短波辐射（ra_sw_physics）
     
187. radt                     = 30,       30,       30,!#30,       30,       30,
     
188. sf_sfclay_physics        = 1,        1,        1, !#近地面层 (sf_sfclay_physics)
     
189. sf_surface_physics       = 2,        2,        2, !#陆面层（sf_surface_physics）
     
190. bl_pbl_physics           = 1,        1,        1, !#行星边界层（bl_pbl_physics）
     
191. bldt                     = 0,        0,        0,
     
192. cu_physics               = 1,        1,        0, !#积云参数化（cu_physics）
     
193. cudt                     = 5,        5,        5,
     
194. isfflx                   = 1,
     
195. ifsnow                   = 0,
     
196. icloud                   = 1,
     
197. surface_input_source     = 1,
     
198. num_soil_layers          = 4, !#指定陆面模式中的土壤层数，默认值为 5   #5         热量扩散方案   #4 Noah 陆面过程方案   #6  RUC 陆面过程方案
     
199. sf_urban_physics         = 0,        0,        0,
     
200. maxiens                  = 1,
     
201. maxens                   = 3,
     
202. maxens2                  = 3,
     
203. maxens3                  = 16,
     
204. ensdim                   = 144,
     
205. num_land_cat             = 20,      !number of land categories in input data.
     
206.                                                !   24 - for USGS (default); 20 for MODIS
     
207.                                                !   28 - for USGS if including lake category
     
208.                                                !   21 - for MODIS if including lake category
     
209. hail_opt                 = 1,!#morrison方案中afwa_hail_opt=1为输出膘默认0输出霰. WRFV3.6.1后变为hail_opt
     
210. /
     
211. 
     
212. &fdda                    
     
213. /
     
214. 
     
215. &dynamics               
     
216. w_damping                = 0,
     
217. diff_opt                 = 1,
     
218. km_opt                   = 4,
     
219. diff_6th_opt             = 0,        0,        0,
     
220. diff_6th_factor          = 0.12,     0.12,     0.12,
     
221. base_temp                = 290.,
     
222. damp_opt                 = 0,
     
223. zdamp                    = 5000.,    5000.,    5000.,
     
224. dampcoef                 = 0.2,      0.2,      0.2,
     
225. khdif                    = 0,        0,        0,
     
226. kvdif                    = 0,        0,        0,
     
227. non_hydrostatic          = .true.,   .true.,   .true.,
     
228. moist_adv_opt            = 1,        1,        1,
     
229. scalar_adv_opt           = 1,        1,        1,
     
230. /
     
231. 
     
232. &bdy_control            
     
233. spec_bdy_width           = 5,
     
234. spec_zone                = 1,
     
235. relax_zone               = 4,
     
236. specified                = .true.,  .false.,  .false.,
     
237. nested                   = .false.,   .true.,   .true.,
     
238. /
     
239. 
     
240. &grib2                  
     
241. /
     
242. 
     
243. &namelist_quilt         
     
244. nio_tasks_per_group      = 0,
     
245. nio_groups               = 1,
     
246. /                                                                                                                        
     

复制代码

解决办法：hail_opt=1代表hail，默认是0代表graupel
所以我后来在namelist.input里设置hail_opt=1，这时候转出来的QGRAUP就代表冰雹
但是它的名称还是QGRAUP，而不是自己变为QHAIL

mjy163163

我正好在看官方手册。。
http://www2.mmm.ucar.edu/wrf/users/docs/user_guide_V3/users_guide_chap5.htm#Phys

不知道是不是变量名称的原因，里面写hail_opt是new name since v3.7，3.6.1版对应的是afwa_hail_opt

小其其格

mjy163163 发表于 2019-8-13 14:10
我正好在看官方手册。。
http://www2.mmm.ucar.edu/wrf/users/docs/user_guide_V3/users_guide_chap5.htm# ...

previously afwa_hail_opt, new in V3.6.1

后面还有这句话，因此3.6.1及之后是hail_opt

主要是出不了结果

Sunshining

楼主，我想问一下，改冰雹方案的时候设置hail_opt = 1,是相当于打开hail吗？默认的是graupel吗？我用的是WRFV3.9.1，手册上是这样的并没有说明设置几代表什么，所以我想问一下，这个选项分别设置几表示什么意思呢？如果不设置模式默认的是什么呢？C:\Users\temp\Desktop\TIM截图20190814235431.png

Sunshining

Sunshining 发表于 2019-8-14 23:56
楼主，我想问一下，改冰雹方案的时候设置hail_opt = 1,是相当于打开hail吗？默认的是graupel吗？我用的是WR ...

手册的截图没有发上来，手册对这个选项的描述第二列是空白的，没有说明设置几代表什么意思

小其其格

Sunshining 发表于 2019-8-15 00:00
手册的截图没有发上来，手册对这个选项的描述第二列是空白的，没有说明设置几代表什么意思

不好意思，前几天出差，才看到。WRFV的注册表里有说明的，hail_opt=1代表hail，默认是0代表graupel

Sunshining

小其其格 发表于 2019-8-16 22:09
不好意思，前几天出差，才看到。WRFV的注册表里有说明的，hail_opt=1代表hail，默认是0代表graupel

好的，谢谢了~