求问修改URBAPAR.TBL后为何无法运行BEP、BEM

tbag

Kecheng 发表于 2025-7-27 15:57
此外，我使用的input数据，气象数据是fnl的，地理数据用的default的MODIS（正如我另外一个帖子所言，我貌 ...

我看你是用了4个domain, 4个domain分辨率分别是多少？

另外sf_urban_physics=2用的只是BEP, 不是BEP+BEM.

tbag

Kecheng 发表于 2025-7-27 15:44
首先是URBPARM.TBL：
# The parameters in this table may vary greatly from city to city.
# The def ...

我注意到你更改了在URBPARM.TBL里的很多参数（好事）

首先你测试一下初始的URBPARM.TBL看看能不能让模式跑起来，如果能跑，那就说明改变厚的参数与BEPBEM有冲突。

另外没有31-33的类型情况下，urban landuse 实际上就只有high-residential。那么Building Heights 那里就用的是“ Building Heights: 3”.

我建议两个事：
1. 使用我的eta_level
e_vert                              = 55,     55,     55,
eta_levels                          = 1.,          0.999461463, 0.999281951, 0.999102439,
                                       0.998922927, 0.998743415, 0.99748677,
                                       0.996230185, 0.994973600, 0.993716955, 0.992334723,
                                       0.990814209, 0.989141703, 0.987301886, 0.985278130
                                       0.983051956, 0.980603218, 0.977909565, 0.974946558,
                                       0.971687257, 0.968101978, 0.964158237, 0.959820092,
                                       0.955048144, 0.949799001, 0.944024920, 0.937673509,
                                       0.930686891, 0.923001587, 0.914547801, 0.905248582,
                                       0.895019472, 0.883767486, 0.871390283, 0.857775331,
                                       0.842798889, 0.826324821, 0.808203340, 0.788269699,
                                       0.766342700, 0.742223024, 0.715691328, 0.68650645,
                                       0.654403150, 0.619089544, 0.580244482, 0.537514985,
                                       0.490512520, 0.438809812, 0.381936818, 0.319376528,
                                       0.250560224, 0.174862280, 0.091594524, 0.000000000,

2. 更改phys/module_sf_bep_bem.F, 第一页中, 把nz_um改成30, dz_u改为3, nbui_max=20.
并且在2410行中， 把" do iz_u=1, nz_u" 改成 "do iz_u=1, nz_um-1"

另外你的分辨率是多少？

Kecheng

tbag 发表于 2025-7-27 22:25
我看你是用了4个domain, 4个domain分辨率分别是多少？

另外sf_urban_physics=2用的只是BEP, 不是BEP+B ...

4个domian的分辨率是：
dx                                  = 24300, 8100, 2700, 900,
dy                                  = 24300, 8100, 2700, 900,

另外关于sf_urban_physics=2：对的，我目前测试的是BEP，不过BEP和BEM都会运行到一半突然停止

Kecheng

tbag 发表于 2025-7-27 22:43
我注意到你更改了在URBPARM.TBL里的很多参数（好事）

首先你测试一下初始的URBPARM.TBL看看能不能让模 ...

初始的URBPARM.TBL可以跑起来，修改之后就有问题跑着跑着就停了

好的，我明天按照您的建议去尝试一下，非常感谢！

tbag

Kecheng 发表于 2025-7-27 22:49
初始的URBPARM.TBL可以跑起来，修改之后就有问题跑着跑着就停了

好的，我明天按照您的建议去尝试一下 ...

那就是100%是URBPARM.TBL问题了

一项一项去做排除法吧。

Kecheng

tbag 发表于 2025-7-28 03:54
那就是100%是URBPARM.TBL问题了

一项一项去做排除法吧。

您好，经过我的测试是BUILDING HEIGHTS: 3出现了问题
默认的是：
BUILDING HEIGHTS: 3
#      (sf_urban_physics=2,3)

#     height   Percentage
#      [m]       [%]
       5.0       0.0
      10.0       0.0
      15.0      10.0
      20.0      25.0
      25.0      40.0
      30.0      25.0
      35.0       0.0
END BUILDING HEIGHTS

我修改之后的是：
BUILDING HEIGHTS: 3
#      (sf_urban_physics=2,3)

#     height   Percentage
#      [m]       [%]
       3.0       0.0
       6.0       7.0
       9.0       17.0
       12.0      17.0
       15.0      16.0
       18.0      17.0
       21.0      8.0
       24.0      5.0
       27.0      4.0
       30.0      3.0
       33.0      0.0
       36.0      0.0
       39.0      0.0
       42.0      0.0
       45.0      0.0
       48.0      0.0
       51.0      0.0
       54.0      6.0
END BUILDING HEIGHTS
想请问一下您知道这是为什么吗？是需要调整一些bep、bem源代码里的参数吗？

tbag

Kecheng 发表于 2025-7-28 11:12
您好，经过我的测试是BUILDING HEIGHTS: 3出现了问题
默认的是：
BUILDING HEIGHTS: 3

你要是这么说我就知道为什么了

首先默认的height的高度是都是以5的倍数起的，5，10，15，20...
但是你改的是以3的倍数起，3-54，如果想规避掉问题的话， 把phys/module_sf_bep_bem.F, 第一页中把 dz_u从5改为3， 我推荐把nz_um改成30但是不改应该也是也能运行的（默认的21对于dz_u=3太小了），另外我说的那个2410行的那个也最好改掉

改完之后重新clean 然后recompile, 应该就能跑了


ps: sf_urban_physics要用3才有用

tbag

Kecheng 发表于 2025-7-28 11:12
您好，经过我的测试是BUILDING HEIGHTS: 3出现了问题
默认的是：
BUILDING HEIGHTS: 3

哦对了，我的版本是基于4.5.1的。

Kecheng

tbag 发表于 2025-7-28 11:31
哦对了，我的版本是基于4.5.1的。

好嘞好嘞，我的wrf版本是4.5.2，应该差别不大，我根据您的建议去试试，真的太感谢您了

Kecheng

tbag 发表于 2025-7-28 11:29
你要是这么说我就知道为什么了

首先默认的height的高度是都是以5的倍数起的，5，10，15，20...

我修改完并且重新编译了，目前可以正常运行了，非常感谢您！