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

Kecheng

大家好，

我在运行BEP和BEM的过程中发现，如果我使用 （链接我尝试能否放评论区）Improving the WRF/urban modeling system in China by developing a national urban dataset 这篇文章中的URBAPAR.TBL，我在wrf运行的过程中会突然无报错地停止（通过在module_surface_driver.F里添加输出，定位到可能是调用module_sf_noahdrv.F里的lsm子程序出问题），可是如果我使用默认的URBAPAR.TBL则可以正常运行，这令我感到很困惑。

不知道是否有人知道这该如何处理，非常感谢大家的建议

tbag

另外你说话方式好像外国人

tbag

namelist.input
URBPARM.TBL
都上传一下，另外WRF是在什么时候中断的，用 grep -i error rsl.* 和 grep -i fatal rsl.* 检查一下有没有报错信息。

另外告诉我下你使用的是什么input数据。

tbag

根据你的描述，大概率是你的URBPARM.TBL哪里写错了，建议仔细检查，注意下逗号之类的东西，不要额外添加，比如每个参数最后是没有逗号的

Kecheng

tbag 发表于 2025-7-26 23:20
根据你的描述，大概率是你的URBPARM.TBL哪里写错了，建议仔细检查，注意下逗号之类的东西，不要额外添加， ...

我用vimdiff检查了一下，修改后的URBPARM.TBL除了数值修改了别的没什么问题，而且SLUCM可以正常运行，就是到BEP/BEM就不行了

Kecheng

tbag 发表于 2025-7-26 23:20
另外你说话方式好像外国人

hhh可能是因为我一开始准备在WRF & MPAS-A Support Forum上也问一下，后面我直接翻译过来了

Kecheng

tbag 发表于 2025-7-26 23:19
namelist.input
URBPARM.TBL
都上传一下，另外WRF是在什么时候中断的，用 grep -i error rsl.* 和 grep - ...

首先是URBPARM.TBL：
# The parameters in this table may vary greatly from city to city.
# The default values are probably not appropriate for any given city.
# Users should adapt these values based on the city they are working
# with.

# Urban Parameters depending on Urban type
# USGS

Number of urban categories: 3

#
#  Where there are multiple columns of values, the values refer, in
#  order, to: 1) Low density residential, 2) High density residential,
#  and 3) Commercial:  I.e.:
#
#  Index:     1           2              3
#  Type:  Low-dens Res, Hi-dens Res, Commercial
#

#
# ZR:  Roof level (building height)  [ m ]
#      (sf_urban_physics=1)

ZR: 5.5, 8.6, 18

#
# SIGMA_ZED:  Standard Deviation of roof height  [ m ]
#      (sf_urban_physics=1)

SIGMA_ZED: 1.0,  3.0,  4.0

#
# ROOF_WIDTH:  Roof (i.e., building) width  [ m ]
#      (sf_urban_physics=1)

ROOF_WIDTH: 18.7, 21, 19.8

#
# ROAD_WIDTH:  road width  [ m ]
#      (sf_urban_physics=1)
#

ROAD_WIDTH: 28.9, 16.6, 12.3

#
# AH:  Anthropogenic heat [ W m{-2} ]
#      (sf_urban_physics=1)
#

AH: 11.24, 17.411, 34.765


#
# ALH:  Anthropogenic latent heat [ W m{-2} ]
#      (sf_urban_physics=1)
#

ALH:  20.0, 25.0, 40.0

#
#  AKANDA_URBAN:  Coefficient modifying the Kanda approach to computing
#  surface layer exchange coefficients.
#      (sf_urban_physics=1)

AKANDA_URBAN:  1.29 1.29 1.29

#
# DDZR:  Thickness of each roof layer [ m ]
#        This is currently NOT a function urban type, but a function
#        of the number of layers.  Number of layers must be 4, for now.
#      (sf_urban_physics=1)


DDZR:  0.05, 0.05, 0.05, 0.05

#
# DDZB:  Thickness of each building wall layer [ m ]
#        This is currently NOT a function urban type, but a function
#        of the number of layers.  Number of layers must be 4, for now.
#      (sf_urban_physics=1)
#

DDZB: 0.05, 0.05, 0.05, 0.05

#
# DDZG:  Thickness of each ground (road) layer [ m ]
#        This is currently NOT a function urban type, but a function
#        of the number of layers.  Number of layers must be 4, for now.
#      (sf_urban_physics=1)
#

DDZG: 0.05, 0.25, 0.50, 0.75

#
# BOUNDR:  Lower boundary condition for roof layer temperature [ 1: Zero-Flux,  2: T = Constant ]
#      (sf_urban_physics=1)
#

BOUNDR: 1

#
# BOUNDB:  Lower boundary condition for wall layer temperature [ 1: Zero-Flux,  2: T = Constant ]
#      (sf_urban_physics=1)
#

BOUNDB: 1

#
# BOUNDG:  Lower boundary condition for ground (road) layer temperature [ 1: Zero-Flux,  2: T = Constant ]
#      (sf_urban_physics=1)
#

BOUNDG: 1

#
# Ch of Wall and Road [ 1: M-O Similarity Theory, 2: Empirical Form of Narita et al., 1997 (recommended) ]
#      (sf_urban_physics=1)
#

CH_SCHEME: 2

#
# Surface and Layer Temperatures [ 1: 4-layer model,  2: Force-Restore method ]
#      (sf_urban_physics=1)
#

TS_SCHEME: 1

#
# AHOPTION [ 0: No anthropogenic heating,  1: Anthropogenic heating will be added to sensible heat flux term ]
#      (sf_urban_physics=1)
#

AHOPTION: 1

#
# Anthropogenic Heating diurnal profile.
#   Multiplication factor applied to AH (as defined in the table above)
#   Hourly values ( 24 of them ), starting at 01 hours Local Time.
#   For sub-hourly model time steps, value changes on the hour and is
#   held constant until the next hour.
#      (sf_urban_physics=1)
#

AHDIUPRF: 0.286 0.223 0.21 0.18 0.209 0.234 0.321 0.449 0.531 0.73 0.93 0.962 0.966 0.961 0.971 0.995 1 0.996 0.959 0.794 0.627 0.559 0.416 0.343

#
# ALHOPTION [ 0: No anthropogenic latent heat,  1: Anthropogenic heating will be added to latent heat flux term ]
#      (sf_urban_physics=1)
#

ALHOPTION: 0

#
# Anthropogenic latent heat: seasonal coefficient of daily maximum values
# From left to right in order: Spring (MAM), Summer(JJA), Fall(SON), Winter(DJF)
#      (sf_urban_physics=1)
#

ALHSEASON: 0.43 1.00 0.54 0.40

#
# Anthropogenic latent heat diurnal profile.
#   Multiplication factor applied to seasonal ALH (as defined above)
#   Half-hourly values ( 48 of them ), starting at 00:30 hours Local Time.
#      (sf_urban_physics=1)
#

ALHDIUPRF: 0.436 0.421 0.391 0.356 0.311 0.301 0.306 0.295 0.253 0.205 0.177 0.162 0.148 0.121 0.118 0.146 0.210 0.250 0.227 0.162 0.127 0.184 0.306 0.413 0.487 0.559 0.639 0.728 0.754 0.812 0.867 0.969 1.000 0.949 0.840 0.775 0.758 0.756 0.706 0.658 0.637 0.632 0.636 0.633 0.639 0.615 0.553 0.485

# Oasis effect
#  Multiplication factor applied to potential ET of vegetation in urban areas
#  Value should be larger than 1 when actived
#  (sf_urban_physics=1)

OASIS: 1.0

# Evaporation scheme for impervious surfaces (for roof, wall, and road)
#            [1: Hypothesized evaporation during large rainfall events (Original)
#            [2: Water-holding scheme over impervious surface, Yang et al., 2014
#      (sf_urban_physics=1)

IMP_SCHEME: 1

# Porosity of pavement materials on impervious surface
# For calculating latent heat flux over impervious surface
# From left to right in order: roof, wall, road
#      (sf_urban_physics=1,IMP_SCHEME=2)
#

PORIMP: 0.45 0.45 0.45

# Maximum water-holding depth of pavement materials on impervious surface [m]
# For calculating latent heat flux over impervious surface
# From left to right in order: roof, wall, road
#      (sf_urban_physics=1,IMP_SCHEME=2)
#

DENGIMP: 0.001 0.0002 0.001

# Urban irrigation scheme, for vegetation in urban area and green roof
#            [0: No irrigation
#            [1: Summertime (May-Sep) irrigation everyday at 9pm
#      (sf_urban_physics=1)

IRI_SCHEME: 0

#
# GROPTION [ 0: No green roof,  1: Enable green roof simulation]
#      (sf_urban_physics=1)
#

GROPTION: 0

# Surface fraction of green roof over urban rooftop (0-1)
#  (sf_urban_physics=1)
#

FGR: 0.0

#
# DZGR:  Thickness of each layer on green roof [ m ]
#  Green roof structure: 4-layers
#  1: Top Soil layer    2:Soil layer      3: Growing Medium layer     
#  4: concrete roof (depth depends on DDZR defined earlier in this table)
#      (sf_urban_physics=1)

DZGR:  0.05 0.10 0.15 0.20

#
# FRC_URB:  Fraction of the urban landscape which does not have natural
#           vegetation. [ Fraction ]
#      (sf_urban_physics=1,2,3)
#

FRC_URB: 0.829, 0.892, 0.919

#
# CAPR:  Heat capacity of roof  [ J m{-3} K{-1} ]
#      (sf_urban_physics=1,2,3)
#

CAPR: 1.0E6, 1.0E6, 1.0E6,

#
# CAPB:  Heat capacity of building wall [ J m{-3} K{-1} ]
#      (sf_urban_physics=1,2,3)
#

CAPB: 1.0E6, 1.0E6, 1.0E6,

#
# CAPG:  Heat capacity of ground (road) [ J m{-3} K{-1} ]
#      (sf_urban_physics=1,2,3)
#

CAPG:  1.4E6, 1.4E6, 1.4E6,

#
# AKSR:  Thermal conductivity of roof [ J m{-1} s{-1} K{-1} ]
#      (sf_urban_physics=1,2,3)
#

AKSR:  0.67, 0.67, 0.67,

#
# AKSB:  Thermal conductivity of building wall [ J m{-1} s{-1} K{-1} ]
#      (sf_urban_physics=1,2,3)
#

AKSB:  0.67, 0.67, 0.67,

#
# AKSG:  Thermal conductivity of ground (road) [ J m{-1} s{-1} K{-1} ]
#      (sf_urban_physics=1,2,3)
#

AKSG: 0.4004, 0.4004, 0.4004,

#
# ALBR:   Surface albedo of roof [ fraction ]
#      (sf_urban_physics=1,2,3)
#

ALBR: 0.20, 0.20, 0.20

#
# ALBB:  Surface albedo of building wall [ fraction ]
#      (sf_urban_physics=1,2,3)
#

ALBB: 0.20, 0.20, 0.20

#
# ALBG:  Surface albedo of ground (road) [ fraction ]
#      (sf_urban_physics=1,2,3)
#

ALBG: 0.20, 0.20, 0.20

#
# EPSR:  Surface emissivity of roof [ - ]
#      (sf_urban_physics=1,2,3)
#

EPSR: 0.90, 0.90, 0.90

#
# EPSB:  Surface emissivity of building wall [-]
#      (sf_urban_physics=1,2,3)
#

EPSB: 0.90, 0.90, 0.90

#
# EPSG:  Surface emissivity of ground (road) [ - ]
#      (sf_urban_physics=1,2,3)
#

EPSG: 0.95, 0.95, 0.95

#
# Z0B:  Roughness length for momentum, over building wall [ m ]
#       Only active for CH_SCHEME == 1
#      (sf_urban_physics=1)
#

Z0B: 0.0001, 0.0001, 0.0001

#
# Z0G:  Roughness length for momentum, over ground (road) [ m ]
#       Only active for CH_SCHEME == 1
#      (sf_urban_physics=1,2,3)
#

Z0G: 0.01, 0.01, 0.01

#
# Z0R:  Roughness length for momentum over roof [ m ]
#      (sf_urban_physics=2,3)
#

Z0R: 0.01, 0.01, 0.01

#
# TRLEND:  Lower boundary condition for roof temperature [ K ]
#      (sf_urban_physics=1,2,3)
#

TRLEND: 293.00, 293.00, 293.00

#
# TBLEND:  Lower boundary temperature for building wall temperature [ K ]
#      (sf_urban_physics=1,2,3)
#

TBLEND: 293.00, 293.00, 293.00

#
# TGLEND:  Lower boundary temperature for ground (road) temperature [ K ]
#      (sf_urban_physics=1,2,3)
#

TGLEND: 293.00, 293.00, 293.00

#
# COP:  Coefficient of performance of the A/C systems [ - ]
#      (sf_urban_physics=3)
#

COP: 3.5, 3.5, 3.5

#
# BLDAC_FRC: fraction of buildings installed with A/C systems [ - ]
#      (sf_urban_physics=3)
#

BLDAC_FRC: 1.0, 1.0, 1.0

#
# COOLED_FRC: fraction of cooled floor area in buildings [ - ]
#      (sf_urban_physics=3)
#
COOLED_FRC: 1.0, 1.0, 1.0

#
# PWIN:  Coverage area fraction of windows in the walls of the building [ - ]
#      (sf_urban_physics=3)
#

PWIN: 0.2, 0.2, 0.2

#
# BETA:  Thermal efficiency of heat exchanger
#      (sf_urban_physics=3)
#

BETA: 0.75, 0.75, 0.75

#
# SW_COND:  Air conditioning switch, 1=ON
#      (sf_urban_physics=3)
#

SW_COND: 1, 1, 1

#
# TIME_ON:  Initial local time of A/C systems, [ h ]
#      (sf_urban_physics=3)
#

TIME_ON: 0., 0., 0.

#
# TIME_OFF:  End local time of A/C systems, [ h ]
#      (sf_urban_physics=3)
#

TIME_OFF: 24., 24., 24.

#
# TARGTEMP:  Target Temperature of the A/C systems, [ K ]
#      (sf_urban_physics=3)
#

TARGTEMP: 298., 298., 297.

#
# GAPTEMP:  Comfort Range of the indoor Temperature, [ K ]
#      (sf_urban_physics=3)
#

GAPTEMP: 0.5, 0.5, 0.5

#
# TARGHUM:  Target humidity of the A/C systems, [ Kg/Kg ]
#      (sf_urban_physics=3)
#

TARGHUM: 0.005, 0.005, 0.005

#
# GAPHUM:  Comfort Range of the specific humidity, [ Kg/Kg ]
#      (sf_urban_physics=3)
#

GAPHUM: 0.005, 0.005, 0.005

#
# PERFLO:  Peak number of occupants per unit floor area, [ person/m^2 ]
#      (sf_urban_physics=3)
#

PERFLO: 0.01, 0.01, 0.02

#
# HSEQUIP:  Diurnal heating profile of heat generated by equipments
#      (sf_urban_physics=3)
#

HSEQUIP: 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.5 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 0.5 0.25 0.25 0.25 0.25 0.25

#
# HSEQUIP_SCALE_FACTOR:  Peak heat generated by equipments, [ W/m^2 ]
#      (sf_urban_physics=3)
#

HSEQUIP_SCALE_FACTOR: 16.00, 20.00, 36.00

#
# GR_FLAG: 1 to switch on green roof model (0-1)
#      (sf_urban_physics=3)
#

GR_FLAG:0

#
# GR_TYPE: 1 for GRASS, 2 for SEDUM vegetation on the green roof
#      (sf_urban_physics=3)
#

GR_TYPE: 2

#
# GR_FRAC_ROOF: fraction of green roof over the roof (0:1)
#      (sf_urban_physics=3)
#

GR_FRAC_ROOF:0,0,0

#
# HSEQUIP:  Diurnal  profile of sprinkler irrigation for the green roof
#      (sf_urban_physics=3)
#

IRHO:0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1


#
# PV_FRAC_ROOF: fraction of photovoltaic panels over the roof (0:1)
#      (sf_urban_physics=3)
#

PV_FRAC_ROOF: 0,0,0




STREET PARAMETERS:
#      (sf_urban_physics=2,3)

#  urban      street      street     building
# category  direction     width      width
# [index]  [deg from N]    [m]        [m]

    1         0.0          28.9       18.7
    1        90.0          28.9       18.7
    2         0.0          16.6       21
    2        90.0          16.6       21
    3         0.0          12.3       19.8
    3        90.0          12.3       19.8

END STREET PARAMETERS

BUILDING HEIGHTS: 1
#      (sf_urban_physics=2,3)

#     height   Percentage
#      [m]       [%]
       3.0       54.0
       6.0       31.0
       9.0       8.0
       12.0      4.0
       15.0      0.0
       18.0      3.0
END BUILDING HEIGHTS

BUILDING HEIGHTS: 2
#      (sf_urban_physics=2,3)

#     height   Percentage
#      [m]       [%]
       3.0       17.0
       6.0       38.0
       9.0       23.0
       12.0      11.0
       15.0      5.0
       18.0      6.0
       21.0      0.0
       24.0      0.0
       27.0      0.0
       30.0      0.0
       33.0      0.0
       36.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

Kecheng

tbag 发表于 2025-7-26 23:19
namelist.input
URBPARM.TBL
都上传一下，另外WRF是在什么时候中断的，用 grep -i error rsl.* 和 grep - ...

然后是namelist：
&time_control
force_use_old_data                  = .true.,
/

&domains
e_vert                              = 45,    45,   45,   45,
p_top_requested                     = 5000,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
feedback                            = 0,
smooth_option                       = 0,
sfcp_to_sfcp                        =.true.,
/

&physics
mp_physics                          = 2,     2,     2,     2,
ra_lw_physics                       = 1,     1,     1,     1,
ra_sw_physics                       = 2,     2,     2,     2,
radt                                = 30,    30,    30,   30,
sf_sfclay_physics                   = 1,     1,     1,     1,
sf_surface_physics                  = 2,     2,     2,     2,
bl_pbl_physics                      = 1,     1,     1,     1,
bldt                                = 0,     0,     0,     0,
cu_physics                          = 3,     3      3,     3,
cudt                                = 5,     5,     5,     5,
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 1,
num_land_cat                        = 21,
num_soil_layers                     = 4,
sf_urban_physics                    = 2,     2,     2,     2,
num_urban_hi                        = 15,
/


&dynamics
hybrid_opt                          = 2,
rk_ord                              = 3,
w_damping                           = 1,
diff_opt                            = 1,      1,      1,     1,
km_opt                              = 4,      4,      4,     4,
diff_6th_opt                        = 0,      0,      0,     0,
diff_6th_factor                     = 0.12,   0.12,   0.12,  0.12,
base_temp                           = 290.,
damp_opt                            = 3,
zdamp                               = 5000.,  5000.,  5000., 5000.,
dampcoef                            = 0.21,   0.21,   0.21,  0.21,
khdif                               = 0,      0,      0,     0,
kvdif                               = 0,      0,      0,     0,
non_hydrostatic                     = .true., .true., .true., .true.,
moist_adv_opt                       = 2,      2,      2,     2,
scalar_adv_opt                      = 2,      2,      2,     2,
tke_adv_opt                         = 2,      2,      2,     2,
momentum_adv_opt                     = 1,      1,      1,     1,
/

&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true., .false.,.false., .false.,
nested                              = .false., .true., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

Kecheng

tbag 发表于 2025-7-26 23:19
namelist.input
URBPARM.TBL
都上传一下，另外WRF是在什么时候中断的，用 grep -i error rsl.* 和 grep - ...

没有报错信息
是在：in NOAH DRV 这里停的
我找到这句输出是在module_surface_driver.F里，然后我添加输出定位到具体是在这里停止的（前面的那句Before calling lsm有输出，但是后面的After calling lsm并没有输出）后面经过查找这个子程序在module_sf_noahdrv.F里
               ! Debug output: Before calling lsm (zhaokc)
               CALL wrf_debug(100, 'LSMSCHEME: Before calling lsm')
               CALL lsm(dz8w,qv_curr,p8w,t_phy,tsk,                 &
                hfx,qfx,lh,grdflx,qgh,gsw,swdown,swddir,swddif,     &
                glw,smstav,smstot,    &
                sfcrunoff,udrunoff,ivgtyp,isltyp,isurban,isice,vegfra,        &
                albedo,albbck,znt,z0, tmn,xland,xice,emiss, embck,    &
                snowc,qsfc,rainbl,                              &
                mminlu,                                         &
                num_soil_layers,dtbl,dzs,itimestep,             &
                smois,tslb,snow,canwat,                         &
                chs, chs2, cqs2, cpm,rcp,SR,chklowq,lai,qz0,    &
                myj,frpcpn,                                     &
                sh2o,snowh,                                     & !h
                u_phy,v_phy,                                    & !I
                snoalb,shdmin,shdmax,                           & !i
                snotime,                                        & !o
                acsnom,acsnow,                                  & !o
                snopcx,                                         & !o
                potevp,                                         & !o
                smcrel,                                         & !o
                xice_threshold,                                 &
                rdlai2d,usemonalb,                              &
                br,                                             & !?
                NOAHRES,opt_thcnd,                              &
                ua_phys,flx4,fvb,fbur,fgsn,                     &
                ids,ide, jds,jde, kds,kde,                      &
                ims,ime, jms,jme, kms,kme,                      &
                i_start(ij),i_end(ij), j_start(ij),j_end(ij), kts,kte,    &
                sf_urban_physics                                &
!Optional urban
                ,cmr_sfcdif,chr_sfcdif,cmc_sfcdif,chc_sfcdif    &
                ,cmgr_sfcdif,chgr_sfcdif                        &
                ,tr_urb2d,tb_urb2d,tg_urb2d,tc_urb2d,qc_urb2d,  & !H urban
                uc_urb2d,                                       & !H urban
                xxxr_urb2d,xxxb_urb2d,xxxg_urb2d,xxxc_urb2d,    & !H urban
                trl_urb3d,tbl_urb3d,tgl_urb3d,                  & !H urban
                sh_urb2d,lh_urb2d,g_urb2d,rn_urb2d,ts_urb2d,    & !H urban
                psim_urb2d,psih_urb2d,u10_urb2d,v10_urb2d,      & !O urban
                GZ1OZ0_urb2d, AKMS_URB2D,                       & !O urban
                th2_urb2d,q2_urb2d,ust_urb2d,                   & !O urban
                declin,coszen,hrang,                            & !I solar
                xlat_urb2d,                                     & !I urban
                num_roof_layers, num_wall_layers,               & !I urban
                num_road_layers, DZR, DZB, DZG,                 & !I urban
                CMCR_URB2D,TGR_URB2D,TGRL_URB3D,SMR_URB3D,      & !H urban
                DRELR_URB2D,DRELB_URB2D,DRELG_URB2D,            & !H urban
                FLXHUMR_URB2D,FLXHUMB_URB2D,FLXHUMG_URB2D,      & !H urban
                julian, julyr,                                  & !H urban
                FRC_URB2D, UTYPE_URB2D,                         & !I urban
                num_urban_ndm,                                  & !I multi-layer urban
                urban_map_zrd,                                  & !I multi-layer urban
                urban_map_zwd,                                  & !I multi-layer urban
                urban_map_gd,                                   & !I multi-layer urban
                urban_map_zd,                                   & !I multi-layer urban
                urban_map_zdf,                                  & !I multi-layer urban
                urban_map_bd,                                   & !I multi-layer urban
                urban_map_wd,                                   & !I multi-layer urban
                urban_map_gbd,                                  & !I multi-layer urban
                urban_map_fbd,                                  & !I multi-layer urban
                urban_map_zgrd,                                 & !I multi-layer urban
                num_urban_hi,                                   & !I multi-layer urban
                tsk_rural,                                      & !H multi-layer urban
                trb_urb4d,tw1_urb4d,tw2_urb4d,tgb_urb4d,        & !H multi-layer urban
                tlev_urb3d,qlev_urb3d,                          & !H multi-layer urban
                tw1lev_urb3d,tw2lev_urb3d,                      & !H multi-layer urban
                tglev_urb3d,tflev_urb3d,                        & !H multi-layer urban
                sf_ac_urb3d,lf_ac_urb3d,cm_ac_urb3d,            & !H multi-layer urban
                sfvent_urb3d,lfvent_urb3d,                      & !H multi-layer urban
                sfwin1_urb3d,sfwin2_urb3d,                      & !H multi-layer urban
                sfw1_urb3d,sfw2_urb3d,sfr_urb3d,sfg_urb3d,      & !H multi-layer urban
                ep_pv_urb3d,t_pv_urb3d,                         & !GRZ
                trv_urb4d,qr_urb4d,qgr_urb3d,tgr_urb3d,         & !GRZ
                drain_urb4d,draingr_urb3d,sfrv_urb3d,           & !GRZ
                lfrv_urb3d,dgr_urb3d,dg_urb3d,lfr_urb3d,lfg_urb3d,& !GRZ
                lp_urb2d,hi_urb2d,lb_urb2d,hgt_urb2d,           & !H multi-layer urban
                mh_urb2d,stdh_urb2D,lf_urb2d,                   & !SLUCM
                th_phy,rho,p_phy,ust,                           & !I multi-layer urban
                gmt,julday,xlong,xlat,                          & !I multi-layer urban
                a_u_bep,a_v_bep,a_t_bep,a_q_bep,                & !O multi-layer urban
                a_e_bep,b_u_bep,b_v_bep,                        & !O multi-layer urban
                b_t_bep,b_q_bep,b_e_bep,dlg_bep,                & !O multi-layer urban
                dl_u_bep,sf_bep,vl_bep                          & !O multi-layer urban
#ifdef WRF_HYDRO
                ,sfcheadrt,INFXSRT, soldrain                    &
#endif
                ,SDA_HFX, SDA_QFX, HFX_BOTH, QFX_BOTH, QNORM, fasdas    &
                ,RS,XLAIDYN,IRRIGATION_CHANNEL)
               ! Debug output: After calling lsm (zhaokc)
               CALL wrf_debug(100, 'LSMSCHEME: After calling lsm')
         ENDIF

Kecheng

tbag 发表于 2025-7-26 23:19
namelist.input
URBPARM.TBL
都上传一下，另外WRF是在什么时候中断的，用 grep -i error rsl.* 和 grep - ...

此外，我使用的input数据，气象数据是fnl的，地理数据用的default的MODIS（正如我另外一个帖子所言，我貌似并没有成功把城市13改成31、32、33，所以我还是继续用的13来跑的）。SLUCM是可以正常运行的，但是不知道为什么到BEP和BEM会报错。
另外我也有查看您对于BEP和BEM运行建议的帖子，我尝试过把PBL方案改成MYJ和BouLac，但是也都会出现同样的问题。（我用的WRF4.5.2，城市冠层可以与YSU一起使用）