爱气象，爱气象家园！ 
·设为首页 开启辅助访问 新帖子新回复热门贴精华帖

气象家园

 找回密码
 立即注册

QQ登录

只需一步，快速开始

新浪微博登陆

只需一步, 快速开始

搜索
气象家园»论坛 › 提问专区 › 悬赏问答(所得积分可换礼) › wrf 台风模拟没有任何报错crash
返回列表 发新帖
查看: 608|回复: 1

wrf 台风模拟没有任何报错crash

[复制链接]
OMGTaylor

新浪微博达人勋
发表于 2024-8-13 15:36:12 | 显示全部楼层 |阅读模式
30金钱
运行了一个烟花的工况模拟，但是程序自己停止了，没有任何报错。检查了对应rsl文件，发现是因为出现了NaN值
用的ds083.3数据驱动。
namelist如下：

  1. &time_control
  2. run_days                            = 4,
  3. run_hours                           = 12,
  4. run_minutes                         = 0,
  5. run_seconds                         = 0,

  7. start_year                          = 2021, 2021, 2021,
  8. start_month                         = 07,   07,   07,  
  9. start_day                           = 22,   22,   22,   
  10. start_hour                          = 12,   12,   12,   
  11. start_minute                        = 00,   00,   00,   
  12. start_second                        = 00,   00,   00,   
  13. end_year                            = 2021, 2021, 2021,
  14. end_month                           = 07,   07,   07,   
  15. end_day                             = 27,   27,   27,   
  16. end_hour                            = 00,   00,   00,   
  17. end_minute                          = 00,   00,   00,   
  18. end_second                          = 00,   00,   00,   

  20. interval_seconds                    = 21600,
  21. input_from_file                     = .true.,.true.,.true.,
  22. adjust_output_times                 = .false.
  23. history_interval_h                  = 3,   3,   3,
  24. history_begin_h                                         = 6,
  25. frames_per_outfile                  = 1000,  1000, 1000,
  26. restart                             = .false.,
  27. override_restart_timers             = .false.,
  28. output_ready_flag                   = .true.,
  29. restart_interval                    = 1440,
  30. io_form_history                     = 2,
  31. io_form_restart                     = 2,
  32. io_form_input                       = 2,
  33. io_form_boundary                    = 2,
  34. debug_level                         = 0,
  35. auxinput4_inname = "wrflowinp_d<domain>",
  36. auxinput4_interval = 360,
  37. io_form_auxinput4 = 2,
  38. io_form_auxinput2 = 2,
  39. diag_print = 2,
  40. output_diagnostics = 1,
  41. auxhist3_outname = "wrfxtrm_d<domain>_<date>"
  42. auxhist3_interval = 1440, 1440
  43. frames_per_auxhist3 = 100, 100
  44. io_form_auxhist3 = 2
  45. /

  47. &domains
  48. time_step                           = 90,
  49. time_step_fract_num                 = 0,
  50. time_step_fract_den                 = 1,
  51. max_dom                             = 2,
  52. e_we                                = 259,   445,   532,
  53. e_sn                                = 223,   445,   511,
  54. e_vert                              = 47,   47,   47,
  55. eta_levels(1:47)                    = 1.000, 0.995, 0.990,
  56.                                        0.985, 0.980, 0.970,
  57.                                        0.960, 0.945, 0.930,
  58.                                        0.910, 0.890, 0.855,
  59.                                        0.820, 0.785, 0.750,
  60.                                        0.715, 0.680, 0.645,
  61.                                        0.610, 0.575, 0.540,
  62.                                        0.505, 0.470, 0.440,
  63.                                        0.410, 0.380, 0.350,
  64.                                        0.325, 0.300, 0.275,
  65.                                        0.250, 0.230, 0.210,
  66.                                        0.190, 0.170, 0.155,
  67.                                        0.140, 0.125, 0.110,
  68.                                        0.095, 0.080, 0.065,
  69.                                        0.050, 0.035, 0.020,
  70.                                        0.010, 0.000,
  71. p_top_requested                     = 5000,
  72. num_metgrid_levels                  = 34,
  73. num_metgrid_soil_levels             = 4,
  74. dx                                  = 15000, 5000
  75. dy                                  = 15000, 5000
  76. grid_id                             = 1,     2,     3,   
  77. parent_id                           = 1,     1,     2,   
  78. i_parent_start                      = 1,     74,    37,   
  79. j_parent_start                      = 1,     37,    27,   
  80. parent_grid_ratio                   = 1,     3,     3,   
  81. parent_time_step_ratio              = 1,     3,     3,   
  82. feedback                            = 1,
  83. smooth_option                       = 0,
  84. /

  86. &physics
  87. mp_physics = 24, 24, 24,
  88. !wsm6 wsm7 new Thompson 6 24 38
  89. do_radar_ref                        = 1,
  90. acc_phy_tend                        = 1,     1,     1,
  91. ra_lw_physics                       = 1,     1,     1,         
  92. ra_sw_physics                       = 4,     4,     4,           
  93. radt                                = 15,    5,    30,        
  94. sf_sfclay_physics = 2, 2, 2,           
  95. sf_surface_physics                  = 2,     2,     2,
  96. sf_urban_physics                    = 1,     1,     1,   
  97. num_soil_layers                     = 4,
  98. ifsnow                              = 1,
  99. sst_update                          = 1,
  100. bl_pbl_physics = 2, 2, 2,
  101. cu_physics = 11, 11, 11,
  102. sf_ocean_physics                    = 1,
  103. windfarm_opt                        = 0,     0,     0,  
  104. windfarm_ij                         = 1,
  105. windfarm_wake_model                 = 0, 0, 0,
  106. windfarm_overlap_method             = 0, 0, 0,
  107. !set value to 0 or 2
  108. isftcflx                            = 0,
  109. /

  111. &fdda
  112. grid_fdda = 0, 0, 0,
  113. /

  115. &dynamics
  116. w_damping                           = 0,
  117. diff_opt                            = 1,      1,      1,
  118. km_opt                              = 4,      4,      4,     
  119. mix_isotropic                       = 0,      0,      0,     
  120. diff_6th_factor                     = 0.12,   0.12,   0.12,  
  121. base_temp                           = 290.
  122. damp_opt                            = 3,
  123. zdamp                               = 5000.,  5000.,  5000.,
  124. dampcoef                            = 0.2,    0.2,    0.2,  
  125. non_hydrostatic                     = .true., .true., .true.,
  126. moist_adv_opt                       = 1,      1,      1,  
  127. scalar_adv_opt                      = 1,      1,      1,
  128. /

  130. &bdy_control
  131. spec_bdy_width                      = 5,
  132. spec_zone                           = 1,
  133. relax_zone                          = 4,
  134. specified                           = .true., .false.,.false.,
  135. nested                              = .false., .true., .true.,
  136. /

  138. &namelist_quilt
  139. nio_tasks_per_group = 0,
  140. nio_groups = 1,
  141. /
复制代码
出现NaN的最后一段描述如下：
  1. Timing for main: time 2021-07-22_18:52:00 on domain   2:    0.78287 elapsed seconds
  2. d02   Domain average of dpsdt, dmudt (mb/3h):    412.0000       8.464747       3.039660
  3. d02   Max mu change time step:          250         138  7.1170653E-04
  4. d02   Domain average of dardt, drcdt, drndt (mm/sec):    412.0000      1.7566526E-04  3.0260142E-05  1.4540517E-04
  5. d02   Domain average of rt_sum, rc_sum, rnc_sum (mm):    412.0000       2.854392      0.6870897       2.167303
  6. d02   Max Accum Resolved Precip,   I,J  (mm):    175.4641             210         189
  7. d02   Max Accum Convective Precip,   I,J  (mm):    16.82688             119           1
  8. d02   Domain average of sfcevp, hfx, lh:    412.0000       1.426547       13.25305       143.2779
  9. Timing for main: time 2021-07-22_18:52:30 on domain   2:    0.78248 elapsed seconds
  10. Timing for main: time 2021-07-22_18:52:30 on domain   1:    2.85545 elapsed seconds
  11. d01   Domain average of dpsdt, dmudt (mb/3h):    412.5000       3.344046       1.817358
  12. d01   Max mu change time step:          153          98  6.6300778E-04
  13. d01   Domain average of dardt, drcdt, drndt (mm/sec):    412.5000      1.1633932E-04  4.5772165E-05  7.0567177E-05
  14. d01   Domain average of rt_sum, rc_sum, rnc_sum (mm):    412.5000       1.749507      0.6376105       1.111897
  15. d01   Max Accum Resolved Precip,   I,J  (mm):    161.7724             146          98
  16. d01   Max Accum Convective Precip,   I,J  (mm):    15.61789             224          74
  17. d01   Domain average of sfcevp, hfx, lh:    412.5000       1.024478       10.37695       103.0517
  18. d02   Domain average of dpsdt, dmudt (mb/3h):    412.5000       8.496153       3.050056
  19. d02   Max mu change time step:          250         138  6.3568482E-04
  20. d02   Domain average of dardt, drcdt, drndt (mm/sec):    412.5000      1.7576595E-04  3.0215151E-05  1.4555080E-04
  21. d02   Domain average of rt_sum, rc_sum, rnc_sum (mm):    412.5000       2.859665      0.6879961       2.171669
  22. d02   Max Accum Resolved Precip,   I,J  (mm):    176.4291             210         189
  23. d02   Max Accum Convective Precip,   I,J  (mm):    16.86489             119           1
  24. d02   Domain average of sfcevp, hfx, lh:    412.5000       1.428264       13.25707       143.3070
  25. Timing for main: time 2021-07-22_18:53:00 on domain   2:    0.79433 elapsed seconds
  26. d02   Domain average of dpsdt, dmudt (mb/3h):    413.0000       8.503106       3.050031
  27. d02   Max mu change time step:          263         201  6.1979244E-04
  28. d02   Domain average of dardt, drcdt, drndt (mm/sec):    413.0000      1.7594403E-04  3.0219473E-05  1.4572452E-04
  29. d02   Domain average of rt_sum, rc_sum, rnc_sum (mm):    413.0000       2.864944      0.6889025       2.176042
  30. d02   Max Accum Resolved Precip,   I,J  (mm):    177.3815             210         189
  31. d02   Max Accum Convective Precip,   I,J  (mm):    16.90290             119           1
  32. d02   Domain average of sfcevp, hfx, lh:    413.0000       1.429981       13.26039       143.3237
  33. Timing for main: time 2021-07-22_18:53:30 on domain   2:    0.78092 elapsed seconds
  34. d02   Domain average of dpsdt, dmudt (mb/3h):    413.5000       8.579379       3.048572
  35. d02   Max mu change time step:          348          85  6.5217039E-04
  36. d02   Domain average of dardt, drcdt, drndt (mm/sec):    413.5000      1.7598520E-04  3.0084819E-05  1.4590041E-04
  37. d02   Domain average of rt_sum, rc_sum, rnc_sum (mm):    413.5000       2.870223      0.6898053       2.180418
  38. d02   Max Accum Resolved Precip,   I,J  (mm):    178.2918             210         189
  39. d02   Max Accum Convective Precip,   I,J  (mm):    16.94090             119           1
  40. d02   Domain average of sfcevp, hfx, lh:    413.5000       1.431699       13.26355       143.3456
  41. Timing for main: time 2021-07-22_18:54:00 on domain   2:    0.78127 elapsed seconds
  42. Timing for main: time 2021-07-22_18:54:00 on domain   1:    2.84603 elapsed seconds
  43. d01   Domain average of dpsdt, dmudt (mb/3h):    414.0000       3.364127       1.808495
  44. d01   Max mu change time step:          153          99  6.9467770E-04
  45. d01   Domain average of dardt, drcdt, drndt (mm/sec):    414.0000      1.1663696E-04  4.5807767E-05  7.0829243E-05
  46. d01   Domain average of rt_sum, rc_sum, rnc_sum (mm):    414.0000       1.758777      0.6402324       1.118544
  47. d01   Max Accum Resolved Precip,   I,J  (mm):    162.3967             146          98
  48. d01   Max Accum Convective Precip,   I,J  (mm):    15.68327             224          74
  49. d01   Domain average of sfcevp, hfx, lh:    414.0000       1.028178       10.38052       103.0948
  50. d02   Domain average of dpsdt, dmudt (mb/3h):    414.0000       8.570923       3.050439
  51. d02   Max mu change time step:          348          85  6.1674620E-04
  52. d02   Domain average of dardt, drcdt, drndt (mm/sec):    414.0000      1.7616457E-04  3.0067817E-05  1.4609678E-04
  53. d02   Domain average of rt_sum, rc_sum, rnc_sum (mm):    414.0000       2.875510      0.6907073       2.184802
  54. d02   Max Accum Resolved Precip,   I,J  (mm):    179.1350             210         189
  55. d02   Max Accum Convective Precip,   I,J  (mm):    16.97891             119           1
  56. d02   Domain average of sfcevp, hfx, lh:    414.0000       1.433416       13.26724       143.3632
  57. Timing for main: time 2021-07-22_18:54:30 on domain   2:    0.79838 elapsed seconds
  58. d02   Domain average of dpsdt, dmudt (mb/3h):    414.5000       8.566481       3.041400
  59. d02   Max mu change time step:          211         115  7.9883909E-04
  60. d02   Domain average of dardt, drcdt, drndt (mm/sec):    414.5000      1.7629494E-04  3.0019504E-05  1.4627550E-04
  61. d02   Domain average of rt_sum, rc_sum, rnc_sum (mm):    414.5000       2.880796      0.6916078       2.189189
  62. d02   Max Accum Resolved Precip,   I,J  (mm):    179.8948             210         189
  63. d02   Max Accum Convective Precip,   I,J  (mm):    17.01692             119           1
  64. d02   Domain average of sfcevp, hfx, lh:    414.5000       1.435134       13.27045       143.3763
  65. Timing for main: time 2021-07-22_18:55:00 on domain   2:    0.78296 elapsed seconds
  66. d02   Domain average of dpsdt, dmudt (mb/3h):    415.0000       8.570491       3.039407
  67. d02   Max mu change time step:          251         138  7.3098752E-04
  68. d02   Domain average of dardt, drcdt, drndt (mm/sec):    415.0000      1.7644488E-04  2.9981344E-05  1.4646354E-04
  69. d02   Domain average of rt_sum, rc_sum, rnc_sum (mm):    415.0000       2.886092      0.6925073       2.193584
  70. d02   Max Accum Resolved Precip,   I,J  (mm):    180.5619             210         189
  71. d02   Max Accum Convective Precip,   I,J  (mm):    17.05493             119           1
  72. d02   Domain average of sfcevp, hfx, lh:    415.0000       1.436852       13.27769       143.4006
  73. Timing for main: time 2021-07-22_18:55:30 on domain   2:    1.10315 elapsed seconds
  74. Timing for main: time 2021-07-22_18:55:30 on domain   1:    3.17330 elapsed seconds
  75. d01   Domain average of dpsdt, dmudt (mb/3h):    415.5000                NaN            NaN
  76. d01   Max mu change time step:          153          99  8.0533687E-04
  77. d01   Domain average of dardt, drcdt, drndt (mm/sec):    415.5000      1.1442166E-04  4.5679000E-05  6.8742702E-05
  78. d01   Domain average of rt_sum, rc_sum, rnc_sum (mm):    415.5000       1.767841      0.6428373       1.125003
  79. d01   Max Accum Resolved Precip,   I,J  (mm):    163.1069             146          98
  80. d01   Max Accum Convective Precip,   I,J  (mm):    15.74865             224          74
  81. d01   Domain average of sfcevp, hfx, lh:    415.5000       1.031880       10.38291       103.1288
  82. d02   Domain average of dpsdt, dmudt (mb/3h):    415.5000                NaN            NaN
  83. d02   Max mu change time step:          251         138  6.6997192E-04
  84. d02   Domain average of dardt, drcdt, drndt (mm/sec):    415.5000      1.6938066E-04  3.0039491E-05  1.3934124E-04
  85. d02   Domain average of rt_sum, rc_sum, rnc_sum (mm):    415.5000       2.891173      0.6934085       2.197763
  86. d02   Max Accum Resolved Precip,   I,J  (mm):    181.1331             210         189
  87. d02   Max Accum Convective Precip,   I,J  (mm):    17.09294             119           1
  88. d02   Domain average of sfcevp, hfx, lh:    415.5000       1.438571       13.27940       143.4232
  89. Timing for main: time 2021-07-22_18:56:00 on domain   2:    0.79704 elapsed seconds
复制代码
已经进行过测试，缩短时间步长能够缓解（time step=60）甚至解决（time step=45）
但是计算成本增加太多。
希望大家给出一些建议。

密码修改失败请联系微信：mofangbao
回复

使用道具 举报

李知恩

新浪微博达人勋
发表于 2024-8-13 22:36:27 | 显示全部楼层
我觉得这个应该数据溢出的问题，有几种方法可以尝试一下，一种是减少时间步长，并且 设置ulimit -s 262140, 另一种就是换一种参数化方案，内外层的参数化方案最好都统一一下，可以都尝试一下
密码修改失败请联系微信：mofangbao
回复

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

小黑屋|帮助中心|清除痕迹

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

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

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