请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5074|回复: 6

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

[复制链接]

新浪微博达人勋

发表于 2019-8-7 20:55:15 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 小其其格 于 2020-11-6 08:27 编辑

使用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. &domains                 
  38. time_step                = 180,
  39. time_step_fract_num      = 0,
  40. time_step_fract_den      = 1,
  41. max_dom                  = 3,
  42. e_we                     = 130,      214,      322,
  43. e_sn                     = 112,      196,      313,
  44. e_vert                   = 35,       35,       35,
  45. p_top_requested          = 5000,
  46. num_metgrid_levels       = 32,
  47. num_metgrid_soil_levels  = 4,
  48. dx                       = 27000,     9000,     3000,
  49. dy                       = 27000,     9000,     3000,
  50. grid_id                  = 1,        2,        3,
  51. parent_id                = 1,        1,        2,
  52. i_parent_start           = 1,       30,       80,
  53. j_parent_start           = 1,       24,       25,
  54. parent_grid_ratio        = 1,        3,        3,
  55. parent_time_step_ratio   = 1,        3,        3,
  56. feedback                 = 1,
  57. smooth_option            = 0,
  58. /

  59. &physics                 
  60. mp_physics               = 6,        6,       10, !#微物理过程(mp_physics)
  61. ra_lw_physics            = 1,        1,        1, !#长波辐射(ra_lw_physics)
  62. ra_sw_physics            = 1,        1,        1, !#短波辐射(ra_sw_physics)
  63. radt                     = 30,       30,       30,!#30,       30,       30,
  64. sf_sfclay_physics        = 1,        1,        1, !#近地面层 (sf_sfclay_physics)
  65. sf_surface_physics       = 2,        2,        2, !#陆面层(sf_surface_physics)
  66. bl_pbl_physics           = 1,        1,        1, !#行星边界层(bl_pbl_physics)
  67. bldt                     = 0,        0,        0,
  68. cu_physics               = 1,        1,        0, !#积云参数化(cu_physics)
  69. cudt                     = 5,        5,        5,
  70. isfflx                   = 1,
  71. ifsnow                   = 0,
  72. icloud                   = 1,
  73. surface_input_source     = 1,
  74. num_soil_layers          = 4, !#指定陆面模式中的土壤层数,默认值为 5   #5         热量扩散方案   #4 Noah 陆面过程方案   #6  RUC 陆面过程方案
  75. sf_urban_physics         = 0,        0,        0,
  76. maxiens                  = 1,
  77. maxens                   = 3,
  78. maxens2                  = 3,
  79. maxens3                  = 16,
  80. ensdim                   = 144,
  81. num_land_cat             = 20,      !number of land categories in input data.
  82.                                                !   24 - for USGS (default); 20 for MODIS
  83.                                                !   28 - for USGS if including lake category
  84.                                                !   21 - for MODIS if including lake category
  85. hail_opt                 = 1,!#morrison方案中afwa_hail_opt=1为输出膘默认0输出霰. WRFV3.6.1后变为hail_opt
  86. /

  87. &fdda                    
  88. /

  89. &dynamics               
  90. w_damping                = 0,
  91. diff_opt                 = 1,
  92. km_opt                   = 4,
  93. diff_6th_opt             = 0,        0,        0,
  94. diff_6th_factor          = 0.12,     0.12,     0.12,
  95. base_temp                = 290.,
  96. damp_opt                 = 0,
  97. zdamp                    = 5000.,    5000.,    5000.,
  98. dampcoef                 = 0.2,      0.2,      0.2,
  99. khdif                    = 0,        0,        0,
  100. kvdif                    = 0,        0,        0,
  101. non_hydrostatic          = .true.,   .true.,   .true.,
  102. moist_adv_opt            = 1,        1,        1,
  103. scalar_adv_opt           = 1,        1,        1,
  104. /

  105. &bdy_control            
  106. spec_bdy_width           = 5,
  107. spec_zone                = 1,
  108. relax_zone               = 4,
  109. specified                = .true.,  .false.,  .false.,
  110. nested                   = .false.,   .true.,   .true.,
  111. /

  112. &grib2                  
  113. /

  114. &namelist_quilt         
  115. nio_tasks_per_group      = 0,
  116. nio_groups               = 1,
  117. /                                                                                                                        

  118. run_hours                           = 24
  119. run_minutes                         = 0,
  120. run_seconds                         = 0,
  121. start_year                          = 2019, 2019, 2019,
  122. start_month                         = 03, 03, 03,
  123. start_day                           = 20, 20, 20,
  124. start_hour                          = 12,   12,   18,   
  125. start_minute                        = 00,   00,   00,
  126. start_second                        = 00,   00,   00,   
  127. end_year                            = 2019, 2019, 2019,
  128. end_month                           = 03, 03, 03,
  129. end_day                             = 21, 21, 21,
  130. end_hour                            = 12,   12,   06,
  131. end_minute                          = 00,   00,   00,   
  132. end_second                          = 00,   00,   00,   
  133. interval_seconds                    = 21600
  134. input_from_file                     = .true.,   .true.,   .true.,
  135. history_interval                    = 60,  60, 10,  
  136. frames_per_outfile                  = 1, 1, 1,
  137. restart                             = .false.,
  138. restart_interval                    = 5000,
  139. io_form_history                     = 2
  140. io_form_restart                     = 2
  141. io_form_input                       = 2
  142. io_form_boundary                    = 2
  143. debug_level                         = 0
  144. !########### add new output variables to wrfout
  145. iofields_filename                   = "iofield_list_d01.txt", "iofield_list_d02.txt", "iofield_list_d03.txt",  !#在WRFV3/run下面,里面增加了qh(QHAIL)变量的输出
  146. !#ignore_iofields_warning             = .true.,
  147. !#auxhist0_interval                  = 60,  60, 10,  
  148. !#frames_per_auxhist0                = 1, 1, 1,
  149. !#io_form_auxhist0                   = 2 !#netCDF format
  150. /

  151. &domains                 
  152. time_step                = 180,
  153. time_step_fract_num      = 0,
  154. time_step_fract_den      = 1,
  155. max_dom                  = 3,
  156. e_we                     = 130,      214,      322,
  157. e_sn                     = 112,      196,      313,
  158. e_vert                   = 35,       35,       35,
  159. p_top_requested          = 5000,
  160. num_metgrid_levels       = 32,
  161. num_metgrid_soil_levels  = 4,
  162. dx                       = 27000,     9000,     3000,
  163. dy                       = 27000,     9000,     3000,
  164. grid_id                  = 1,        2,        3,
  165. parent_id                = 1,        1,        2,
  166. i_parent_start           = 1,       30,       80,
  167. j_parent_start           = 1,       24,       25,
  168. parent_grid_ratio        = 1,        3,        3,
  169. parent_time_step_ratio   = 1,        3,        3,
  170. feedback                 = 1,
  171. smooth_option            = 0,
  172. /

  173. &physics                 
  174. mp_physics               = 6,        6,       10, !#微物理过程(mp_physics)
  175. ra_lw_physics            = 1,        1,        1, !#长波辐射(ra_lw_physics)
  176. ra_sw_physics            = 1,        1,        1, !#短波辐射(ra_sw_physics)
  177. radt                     = 30,       30,       30,!#30,       30,       30,
  178. sf_sfclay_physics        = 1,        1,        1, !#近地面层 (sf_sfclay_physics)
  179. sf_surface_physics       = 2,        2,        2, !#陆面层(sf_surface_physics)
  180. bl_pbl_physics           = 1,        1,        1, !#行星边界层(bl_pbl_physics)
  181. bldt                     = 0,        0,        0,
  182. cu_physics               = 1,        1,        0, !#积云参数化(cu_physics)
  183. cudt                     = 5,        5,        5,
  184. isfflx                   = 1,
  185. ifsnow                   = 0,
  186. icloud                   = 1,
  187. surface_input_source     = 1,
  188. num_soil_layers          = 4, !#指定陆面模式中的土壤层数,默认值为 5   #5         热量扩散方案   #4 Noah 陆面过程方案   #6  RUC 陆面过程方案
  189. sf_urban_physics         = 0,        0,        0,
  190. maxiens                  = 1,
  191. maxens                   = 3,
  192. maxens2                  = 3,
  193. maxens3                  = 16,
  194. ensdim                   = 144,
  195. num_land_cat             = 20,      !number of land categories in input data.
  196.                                                !   24 - for USGS (default); 20 for MODIS
  197.                                                !   28 - for USGS if including lake category
  198.                                                !   21 - for MODIS if including lake category
  199. hail_opt                 = 1,!#morrison方案中afwa_hail_opt=1为输出膘默认0输出霰. WRFV3.6.1后变为hail_opt
  200. /

  201. &fdda                    
  202. /

  203. &dynamics               
  204. w_damping                = 0,
  205. diff_opt                 = 1,
  206. km_opt                   = 4,
  207. diff_6th_opt             = 0,        0,        0,
  208. diff_6th_factor          = 0.12,     0.12,     0.12,
  209. base_temp                = 290.,
  210. damp_opt                 = 0,
  211. zdamp                    = 5000.,    5000.,    5000.,
  212. dampcoef                 = 0.2,      0.2,      0.2,
  213. khdif                    = 0,        0,        0,
  214. kvdif                    = 0,        0,        0,
  215. non_hydrostatic          = .true.,   .true.,   .true.,
  216. moist_adv_opt            = 1,        1,        1,
  217. scalar_adv_opt           = 1,        1,        1,
  218. /

  219. &bdy_control            
  220. spec_bdy_width           = 5,
  221. spec_zone                = 1,
  222. relax_zone               = 4,
  223. specified                = .true.,  .false.,  .false.,
  224. nested                   = .false.,   .true.,   .true.,
  225. /

  226. &grib2                  
  227. /

  228. &namelist_quilt         
  229. nio_tasks_per_group      = 0,
  230. nio_groups               = 1,
  231. /                                                                                                                        
复制代码

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



密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-8-13 14:10:27 | 显示全部楼层
我正好在看官方手册。。
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
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-8-13 15:29:47 | 显示全部楼层
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

主要是出不了结果
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-8-14 23:56:47 | 显示全部楼层
楼主,我想问一下,改冰雹方案的时候设置hail_opt = 1,是相当于打开hail吗?默认的是graupel吗?我用的是WRFV3.9.1,手册上是这样的并没有说明设置几代表什么,所以我想问一下,这个选项分别设置几表示什么意思呢?如果不设置模式默认的是什么呢?C:\Users\temp\Desktop\TIM截图20190814235431.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-8-15 00:00:31 | 显示全部楼层
Sunshining 发表于 2019-8-14 23:56
楼主,我想问一下,改冰雹方案的时候设置hail_opt = 1,是相当于打开hail吗?默认的是graupel吗?我用的是WR ...

手册的截图没有发上来,手册对这个选项的描述第二列是空白的,没有说明设置几代表什么意思
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-8-16 22:09:00 | 显示全部楼层
Sunshining 发表于 2019-8-15 00:00
手册的截图没有发上来,手册对这个选项的描述第二列是空白的,没有说明设置几代表什么意思

不好意思,前几天出差,才看到。WRFV的注册表里有说明的,hail_opt=1代表hail,默认是0代表graupel
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-8-17 16:18:39 | 显示全部楼层
小其其格 发表于 2019-8-16 22:09
不好意思,前几天出差,才看到。WRFV的注册表里有说明的,hail_opt=1代表hail,默认是0代表graupel

好的,谢谢了~
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表