WRF模式模拟出现问题，应该比较简单，但是没有检查出来

FrankieLJY

出现的问题：
*************************************
Nesting domain
ids,ide,jds,jde            1          34           1          34
ims,ime,jms,jme           -4          15          -4          15
ips,ipe,jps,jpe            1           5           1           3
INTERMEDIATE domain
ids,ide,jds,jde           15          31          21          37
ims,ime,jms,jme           10          28          16          33
ips,ipe,jps,jpe           13          18          19          23
*************************************
d01 2010-01-01_00:00:00  alloc_space_field: domain            2,                9230412 bytes allocated
d01 2010-01-01_00:00:00  alloc_space_field: domain            2,                1466496 bytes allocated
d01 2010-01-01_00:00:00 calling inc/HALO_INTERP_DOWN_inline.inc
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1,                7329924 bytes allocated
  wrf main: calling open_r_dataset for wrfinput
  med_initialdata_input: calling input_input
  mminlu = 'USGS'
Timing for processing wrfinput file (stream 0) for domain        1:    0.14137 elapsed seconds
一直这个界面，程序一直在跑，但是也没有结果文件出来！

namelist如下：
&time_control
start_year                          = 2010, 2010, 1990, 2012
start_month                         = 01, 01, 01, 01,
start_day = 01, 01, 01,21
start_hour                          = 00, 00, 00, 00,
start_minute                        = 00, 00, 00, 00,
start_second                        = 00, 00, 00, 00,
end_year                            = 2010, 2010, 1990,2012,
end_month                           = 02, 02, 01,07,
end_day = 01, 01,02,22
end_hour                            = 00,00,00,00,
end_minute                          = 00,00,00,00,
end_second                          = 00,00,00,00,
interval_seconds                    = 21600
input_from_file                     = .true.,  .true., .true., .true.,
fine_input_stream                   = 0,       0,      0,      0,
history_interval                    = 60,      60,     60,     60,
frames_per_outfile                  = 24,      24,     24,     24,
restart                             = .false.,
restart_interval                    = 1440,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,
io_form_auxinput2                   = 2,
debug_level                         = 50,
/
&domains
time_step                           = 9,
max_dom                             = 2,
s_we                                = 1,       1,      1,      1,
e_we                                = 45,      34,     892,    283,
s_sn                                = 1,       1,      1,      1,
e_sn                                = 45,      34,     700,    184,
s_vert                              = 1,       1,      1,      1,
e_vert                              = 26,      26,     26,     26,
eta_levels                          = 1.0000, 0.9979, 0.9956, 0.9931, 0.9904,
                                       0.9875, 0.9844, 0.9807, 0.9763, 0.9711,
                                       0.9649, 0.9575, 0.9488, 0.9263, 0.8951,
                                       0.8753, 0.8521, 0.8251, 0.7937, 0.7229,
                                       0.6410, 0.4985, 0.3934, 0.2850, 0.1324,
                                       0.0000,
p_top_requested                     = 5000,
num_metgrid_levels                  = 38,
num_metgrid_soil_levels             = 4,
dx                                  = 18000,6000,6000,1000,
dy                                  = 18000,6000,6000,1000,
grid_id                             = 1,       2,      3,      4,
parent_id                           = 1,       1,      2,      3,
i_parent_start                      = 1,       17,     40,     27,
j_parent_start                      = 1,       23,     34,     12,
parent_grid_ratio                   = 1,       3,      3,      3,
parent_time_step_ratio              = 1,       3,      3,      3,
……
希望有坛友能够指点一下，应该是哪个地方有点简单的问题，但是就是没有检查出来！
我一直觉得是time_step的问题，但是尝试了很多的数值：60/30/36/18/9 都没有成功！

FrankieLJY

time_step应该就是和dx有关，和其他的无关吧。dx是最外层嵌套的分辨率，如果两次模拟的最外层dx相同，那里面嵌套的层数和大小发生改变，会引起有不同吗？希望有明白的坛友给解答一下。

andrewsoong

看着好像是你的内存不够或者cpu有限

FrankieLJY

andrewsoong 发表于 2017-5-20 20:17
看着好像是你的内存不够或者cpu有限

版主，应该不是这样的，我换了一个更大的区域，就可以跑成功，但是跑一个小的区域反而不能跑，都是27km，9km两层嵌套。完全不理解

放逐流年

是不是这个frames_per_outfile                  = 24,      24,     24,     24,的问题呢，调小点应该可以吧

andrewsoong

FrankieLJY 发表于 2017-5-20 20:22
版主，应该不是这样的，我换了一个更大的区域，就可以跑成功，但是跑一个小的区域反而不能跑，都是27km， ...

的确存在一些很奇怪的问题现象

tbag

time step 9??
18的分辨率跑9秒？

Lighting

一般情况下，积分步长为网格分辨率的6倍。你设置的 dx = 18000，即18 km，那么积分步长可以设置为 18*6 或是 18*5，9 s 的时间步长太小，可能会出问题。
PS：建议说一下运行WRF的环境，比如是单机还是集群。如果还有作业输出信息的话也可以给出来，因为有时候导致出现问题的原因比较多，仅从你给出的信息来看，找不到具体的原因。如果是在集群上运行的话，可以让管理员帮忙查看一下。

FrankieLJY

tbag 发表于 2017-5-21 05:13
time step 9??
18的分辨率跑9秒？

从60s降到9s都没有成功，但是我把区域范围扩大，然后跑30s成功了，因为我看一个帖子说时间越大，结果收敛得越快，越不稳定，所以尽量缩小时间。

FrankieLJY

Lighting 发表于 2017-5-21 12:54
一般情况下，积分步长为网格分辨率的6倍。你设置的 dx = 18000，即18 km，那么积分步长可以设置为 18*6 或 ...

我们是自己实验室搭建的一个小集群，我在没有改变dx的情况下，尝试扩大了区域面积，然后就成功运行了，因为之前我在一个帖子上看到有坛友说时间越大收敛得越快，就不稳定。手册上说是dx的3-6倍，但是我现在用30s也成功跑起来了，时间太小会不会也出现结果不准确的现象呢？
这是我之前看到的time_step的帖子：http://bbs.06climate.com/forum.php?mod=viewthread&tid=47382