WRF-CLM 求助

peterchu

在WRF3.7中选择sf_surface_physics=5，也就是陆面模式用CLM模式，需要输入的数据有哪些，为什么real.exe的时候出现 Assume non-CLM input。然后wrf.exe的时候就直接initialize allocation error，然后不做下去了。我用的WPS中namelist.wps:
&share
wrf_core = 'ARW',
max_dom = 2,
start_date = '2016-11-20_06:00:00','2016-11-20_06:00:00',
end_date   = '2016-11-23_06:00:00','2016-11-23_06:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/
&geogrid
parent_id         =   1,   1,
parent_grid_ratio =   1,   3,
i_parent_start    =   1,  73,
j_parent_start    =   1,  34,
e_we              =  99, 28,
e_sn              =  86,  28,
geog_data_res     = 'usgs_lakes+10m','usgs_lakes+2m',
dx = 57600,
dy = 57600,
map_proj = 'lambert',
ref_lat   =  34.272,
ref_lon   = 105.483,
truelat1  =  34.272,
truelat2  =  34.272,
stand_lon = 105.483,
geog_data_path = '../geog'
/
&ungrib
out_format = 'WPS',
prefix = 'FILE',
/
&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/

WRFV3中的namelist.input:

&time_control
run_days                            = 3,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2016, 2016, 2000,
start_month                         = 11,   11,   01,
start_day                           = 20,   20,   24,
start_hour                          = 06,   06,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2016, 2016, 2000,
end_month                           = 11,   11,   01,
end_day                             = 23,   23,   25,
end_hour                            = 06,   06,   12,
end_minute                          = 00,   00,   00,
end_second                          = 00,   00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 360,  60,   60,
frames_per_outfile                  = 1000, 1000, 1000,
restart                             = .false.,
restart_interval                    = 5000,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 0
/
&domains
time_step                           = 180,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 2,
e_we                                = 99,    28,   94,
e_sn                                = 86,    28,    91,
e_vert                              = 40,    40,    40,
p_top_requested                     = 5000,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 57600, 19200,  3333.33,
dy                                  = 57600, 19200,  3333.33,
grid_id                             = 1,     2,     3,
parent_id                           = 1,     1,     2,
i_parent_start                      = 1,     73,    30,
j_parent_start                      = 1,     34,    30,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/
&physics
mp_physics                          = 3,     3,     3,
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 1,     1,     1,
radt                                = 30,    30,    30,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 5,     5,     2,
bl_pbl_physics                      = 1,     1,     1,
bldt                                = 0,     0,     0,
cu_physics                          = 1,     1,     0,
cudt                                = 5,     5,     5,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
num_land_cat                        = 28,
sf_urban_physics                    = 0,     0,     0,
/
&fdda
/
&dynamics
w_damping                           = 0,
diff_opt                            = 1,      1,      1,
km_opt                              = 4,      4,      4,
diff_6th_opt                        = 0,      0,      0,
diff_6th_factor                     = 0.12,   0.12,   0.12,
base_temp                           = 290.
damp_opt                            = 0,
zdamp                               = 5000.,  5000.,  5000.,
dampcoef                            = 0.2,    0.2,    0.2
khdif                               = 0,      0,      0,
kvdif                               = 0,      0,      0,
non_hydrostatic                     = .true., .true., .true.,
moist_adv_opt                       = 1,      1,      1,     
scalar_adv_opt                      = 1,      1,      1,     
/
&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true., .false.,.false.,
nested                              = .false., .true., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

空谷幽竹

楼主用的什么背景场数据？

peterchu

用的是GFS数据“gfs.t06z.pgrb2.1p00.f096”

peterchu

空谷幽竹 发表于 2016-11-18 09:54
楼主用的什么背景场数据？

用的是GFS数据“gfs.t06z.pgrb2.1p00.f096”

空谷幽竹

peterchu 发表于 2016-11-18 10:18
用的是GFS数据“gfs.t06z.pgrb2.1p00.f096”

land_surface_physics 需要特定的num_soil_layers，可检查一下数据

peterchu

空谷幽竹 发表于 2016-11-18 11:56
land_surface_physics 需要特定的num_soil_layers，可检查一下数据

我在num_soil_layers            = 10，; number of soil layers in land surface model
                                                  = 5: thermal diffusion scheme
                                                  = 4: Noah landsurface model
                                                  = 6: RUC landsurface model
                                                  = 10: CLM4 landsurface model
                                                  = 2: Pleim-Xu landsurface model
                                                  = 3: SSiB landsurface model

我把它调成10了，还是不行

空谷幽竹

peterchu 发表于 2016-11-18 12:02
我在num_soil_layers            = 10，; number of soil layers in land surface model
               ...

那你的metgrid后数据是否有１０层呢

cynthianeko

空谷幽竹 发表于 2016-11-18 14:12
那你的metgrid后数据是否有１０层呢

您好，我想请问一下您，metgrid后数据是否有十层的意思是，在使用CLM时num_soil_layers=10时，num_metgrid_soil_levels的值也需要是10吗？wrf初学者，很多都不懂，请您多多包涵

空谷幽竹

cynthianeko 发表于 2018-5-15 10:21
您好，我想请问一下您，metgrid后数据是否有十层的意思是，在使用CLM时num_soil_layers=10时，num_metgri ...

应该是的，两者应该是一样的

cynthianeko

但是直接改成10的话会报一个冲突的错，不知道还需要进行其他修改吗？应该怎么修改呢？例如在met_em文件中该参数仍然是4应该如何处理呢？谢谢您了~