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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 17782|回复: 20

WRF运行real.exe中断问题

[复制链接]

新浪微博达人勋

发表于 2017-3-25 17:34:39 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 waynehu1024 于 2017-3-25 17:36 编辑

WRF版本:3.3.1
气象数据:fnl, 6 hourly
症状:
metgrid得到的met_em文件中,num_metgrid_levels是27
  1. dimensions:
  2.         Time = UNLIMITED ; // (1 currently)
  3.         DateStrLen = 19 ;
  4.         west_east = 63 ;
  5.         south_north = 51 ;
  6.         num_metgrid_levels = 27 ;
  7.         num_st_layers = 4 ;
  8.         num_sm_layers = 4 ;
  9.         south_north_stag = 52 ;
  10.         west_east_stag = 64 ;
  11.         z-dimension0012 = 12 ;
  12.         z-dimension0016 = 16 ;
  13.         z-dimension0033 = 33 ;
复制代码
但是real.exe却报错说size mismatch
  1. FATAL CALLED FROM FILE:  <stdin>  LINE:     693
  2.   input_wrf.F: SIZE MISMATCH:  
  3. namelist ide,jde,num_metgrid_levels=           64           52           27 ;
  4. input data ide,jde,num_metgrid_levels=           64           52           29
复制代码
然而上面的dimensions里明明写着num_metgrid_levels=27。我又看了下met_em里其他参数,发现有个BOTTOM-TOP_GRID_DIMENSION = 29。
  1. // global attributes:
  2.                 :TITLE = "OUTPUT FROM METGRID V3.3.1" ;
  3.                 :SIMULATION_START_DATE = "2015-04-01_00:00:00" ;
  4.                 :WEST-EAST_GRID_DIMENSION = 64 ;
  5.                 :SOUTH-NORTH_GRID_DIMENSION = 52 ;
  6.                 :BOTTOM-TOP_GRID_DIMENSION = 29 ;
  7.                 :WEST-EAST_PATCH_START_UNSTAG = 1 ;
  8.                 :WEST-EAST_PATCH_END_UNSTAG = 63 ;
  9.                 :WEST-EAST_PATCH_START_STAG = 1 ;
  10.                 :WEST-EAST_PATCH_END_STAG = 64 ;
  11.                 :SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
  12.                 :SOUTH-NORTH_PATCH_END_UNSTAG = 51 ;
  13.                 :SOUTH-NORTH_PATCH_START_STAG = 1 ;
  14.                 :SOUTH-NORTH_PATCH_END_STAG = 52 ;
复制代码
于是试探性的把num_metgrid_levels改成29试试。这次报的错不一样了...
  1. Using sfcprs3 to compute psfc
  2. Converged znw(kte) should be about 0.0 =   -4.9360096E-08
  3. all_dim =             1
  4. order =             1
  5. i,j =             1            1
  6. p array =              NaN
  7. f array =     0.000000
  8. p target=              NaN             NaN             NaN             NaN
  9.              NaN             NaN             NaN             NaN
  10.              NaN             NaN             NaN             NaN
  11.              NaN             NaN             NaN             NaN
  12.              NaN             NaN             NaN             NaN
  13.              NaN             NaN             NaN             NaN
  14.              NaN             NaN             NaN             NaN
  15.              NaN
  16. -------------- FATAL CALLED ---------------
  17. FATAL CALLED FROM FILE:  <stdin>  LINE:    4142
  18. troubles, the interpolating order is too large for this few input values
复制代码


搜了下这个症状,有不少人提出来了但是好像都不了了之。有个疑问说p array和p target为什么都是nan...
另外说一下之前的一些错误经历。
我之前用2015年1月6-9号的fnl气象数据成功跑过一个3天的wrf案例,当时完全没有碰到这些问题。但是当我这次用4月的数据模拟时就碰到了一系列奇怪的错误:

1. real.exe时报num_soil_layers=0(和namelist.input不匹配)以及缺少soil temp等问题,试图强行在namelist.input中修正这些错误,强行修正后wrf.exe可以运行,但跑一段时间就会中断,报exit code: 11,网上搜了搜,症状和方案千奇百怪,没有我用了能成功的。


2. 上面使用的ungrib vtable是Vtable.GFS,然而网上看到说GFS在15-01-14有过一次重要更新,要用新的VTable,于是网上下了Vtable.GFS_new。用新的Vtable后没有之前出现的real.exe报土壤层数和温度错误的问题了。


3. 改了Vtable之后,发现metgrid.exe时有很多参数显示missing,但也能正常输出。遇到前面描述过的问题后觉得METGRID.TBL可能也需要更新,在网上找了一份更新过的,多了ptrop之类的参数(没更新之前real.exe的rsl中确实提到了少这些参数),但没有解决根本问题。


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

新浪微博达人勋

 楼主| 发表于 2017-3-25 17:36:01 | 显示全部楼层
附上namelist.input
  1. &time_control            
  2. run_days                 = 3,
  3. run_hours                = 0,
  4. run_minutes              = 0,
  5. run_seconds              = 0,
  6. start_year               = 2015,     2015,     2015,     2015,     2015,
  7. start_month              = 04,       04,       04,       04,       04,
  8. start_day                = 01,       01,       01,       01,       01,
  9. start_hour               = 00,       00,       00,       00,       00,
  10. start_minute             = 00,       00,       00,       00,       00,
  11. start_second             = 00,       00,       00,       00,       00,
  12. end_year                 = 2015,     2015,     2015,     2015,     2015,
  13. end_month                = 04,       04,       04,       04,       04,
  14. end_day                  = 04,       04,       04,       04,       04,
  15. end_hour                 = 00,       00,       00,       00,       00,
  16. end_minute               = 00,       00,       00,       00,       00,
  17. end_second               = 00,       00,       00,       00,       00,
  18. interval_seconds         = 21600,
  19. input_from_file          = .true.,   .true.,   .true.,   .true.,   .true.,
  20. history_interval         = 60,       60,       60,       60,       60,
  21. frames_per_outfile       = 1000,     1000,     1000,     1000,     1000,
  22. restart                  = .false.,
  23. restart_interval         = 5000,
  24. io_form_history          = 2,
  25. io_form_restart          = 2,
  26. io_form_input            = 2,
  27. io_form_boundary         = 2,
  28. debug_level              = 0,
  29. /

  30. &dfi_control
  31. /

  32. &scm
  33. /

  34. &tc
  35. /

  36. &fire
  37. /

  38. &domains      
  39. time_step                = 243,
  40. time_step_fract_num      = 0,
  41. time_step_fract_den      = 1,
  42. max_dom                  = 5,
  43. e_we                     = 64,      88,      88,      88,      121,
  44. e_sn                     = 52,      88,      100,     88,      142,
  45. e_vert                   = 29,      29,      29,      29,      29,
  46. p_top_requested          = 5000,
  47. num_metgrid_levels       = 29,
  48. num_metgrid_soil_levels  = 4,
  49. dx                       = 81000,     27000,     9000,     3000,     1000,
  50. dy                       = 81000,     27000,     9000,     3000,     1000,
  51. grid_id                  = 1,        2,        3,        4,        5,
  52. parent_id                = 1,        1,        2,        3,        4,
  53. i_parent_start           = 1,        34,       40,       44,       29,
  54. j_parent_start           = 1,        5,        32,       42,       21,
  55. parent_grid_ratio        = 1,        3,        3,        3,        3,
  56. parent_time_step_ratio   = 1,        3,        3,        3,        3,
  57. feedback                 = 1,
  58. smooth_option            = 0,
  59. /

  60. &physics                 
  61. mp_physics               = 6,        6,        6,        6,        6,
  62. ra_lw_physics            = 1,        1,        1,        1,        1,
  63. ra_sw_physics            = 1,        1,        1,        1,        1,
  64. radt                     = 30,       30,       30,       30,       30,
  65. sf_sfclay_physics        = 2,        2,        2,        2,        2,
  66. sf_surface_physics       = 1,        1,        1,        1,        1,
  67. bl_pbl_physics           = 2,        2,        2,        2,        2,
  68. bldt                     = 0,        0,        0,        0,        0,
  69. cu_physics               = 2,        2,        2,        0,        0,
  70. cudt                     = 5,        5,        5,        5,        5,
  71. isfflx                   = 1,
  72. ifsnow                   = 0,
  73. icloud                   = 1,
  74. surface_input_source     = 1,
  75. num_soil_layers          = 4,
  76. num_land_cat                  = 33,
  77. sf_urban_physics         = 0,        0,        0,        0,        1,
  78. maxiens                  = 1,
  79. maxens                   = 3,
  80. maxens2                  = 3,
  81. maxens3                  = 16,
  82. ensdim                   = 144,
  83. /

  84. &fdda                    
  85. /

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

  102. &bdy_control            
  103. spec_bdy_width           = 5,
  104. spec_zone                = 1,
  105. relax_zone               = 4,
  106. specified                = .true.,  .false.,  .false.,  .false.,  .false.,
  107. nested                   = .false.,   .true.,   .true.,   .true.,   .true.,
  108. /

  109. &grib2                  
  110. /

  111. &namelist_quilt         
  112. nio_tasks_per_group      = 0,
  113. nio_groups               = 1,
  114. /
  115.                                                                                                                                                      
复制代码
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-3-25 17:38:26 | 显示全部楼层
本帖最后由 waynehu1024 于 2017-3-25 20:56 编辑

尝试过改time_step,从1倍到6倍的dx都试过了,没用...

有趣的是我后来发现我当初跑1月的案例时手滑按了个小数点,把本来要打的time_step = 486打成了48.6,然而最后出来的数据和实测相关很好...
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-3-25 20:54:19 | 显示全部楼层
486是什么参数
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-3-25 20:56:20 | 显示全部楼层

time_step...忘说了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-3-26 09:51:17 | 显示全部楼层
那个小数点好厉害
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-3-26 12:26:46 | 显示全部楼层

别只关注小数点啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-3-26 21:40:15 | 显示全部楼层
又研究了半天...感觉我的版本太低了,打算装个3.7或者3.8再试试,先拿14年的跑了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-3-26 21:52:21 | 显示全部楼层
time_step短一点应该效果要好
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-3-29 15:27:02 | 显示全部楼层
lgz 发表于 2017-3-26 21:52
time_step短一点应该效果要好

应该是我wrf版本太低的问题...
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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