wrf 输出constant field ，怎么回事？

schliezer

小白转WRF又遇到问题了，厚脸皮再求教大家。是这样的，我想看看usgs和modis两种下垫面数据下转出的结果有无大的差别，就用它们分别转了一下，只调了geo_data_res和num_land_cat = 20，其余条件都不变，物理参数化方案，只有长短波选了RRTM和DUHIA，其他方案都是默认（没采用）。
转出来后，USGS的地面气压PSFC不随时间变化，定常的，其余基本气象要素正常；MODIS的就严重得多，降水rainnc+rainc一直为constant field=0，没有下雨，而事实上我模拟的是一次暴雨，其余所有基本气象要素都是不随时间变化的定常场，我的namelist.wps,namelist.input和namelist.arwpost如下，麻烦大神给找找原因，先谢了：
&share
wrf_core = 'ARW',
max_dom = 2,
start_date = '2012-07-28_00:00:00','2012-07-28_00:00:00',
end_date   = '2012-07-31_00:00:00','2012-07-31_00:00:00',
interval_seconds = 21600,
io_form_geogrid = 2,
opt_output_from_geogrid_path = '/home1/zhy/usr/local/WRF/WPS',
debug_level = 0,
/

&geogrid
parent_id         =   1,   1,  
parent_grid_ratio =   1,   3,  
i_parent_start    =   1,  45,
j_parent_start    =   1,  45,
e_we              =  120, 106,
e_sn              =  110, 91,
geog_data_res     = 'modis_30s+5m','modis_30s+30s',
dx = 12000,
dy = 12000,
map_proj = 'lambert',
ref_lat   =  37.75,
ref_lon   = 103.92,
truelat1  =  30.0,
truelat2  =  60.0,
stand_lon = 103.92,
geog_data_path = '/home1/zhy/usr/local/WRF/geog',
ref_x=60.0,
ref_y=55.0,
/

&ungrib
out_format = 'WPS',
prefix = 'FILE',
/

&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
opt_output_from_metgrid_path = '/home1/zhy/usr/local/WRF/WPS',
/
&mod_levs
press_pa=201300,200100,100000,
         95000,90000,
         85000.80000.
         75000,70000,
         65000,60000,
         55000,50000,
         45000,40000,
         35000,30000,
         25000,20000,
         15000,10000,
          5000,1000
/

&time_control
run_days                            = 0,
run_hours                           = 72,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2012, 2012,
start_month                         = 07,   07,
start_day                           = 28,   28,   
start_hour                          = 00,   00,   
start_minute                        = 00,   00,  
start_second                        = 00,   00,
end_year                            = 2012, 2012,
end_month                           = 07,   07,
end_day                             = 31,   31,  
end_hour                            = 00,   00,
end_minute                          = 00,   00,  
end_second                          = 00,   00,  
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,
history_interval                    = 180,  60,  
frames_per_outfile                  = 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                           = 72,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 2,
e_we                                = 120,106,
e_sn                                = 110,91,
e_vert                              = 35, 35,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 12000, 4000,  
dy                                  = 12000, 4000,
grid_id                             = 1,     2,   
parent_id                           = 0,     1,   
i_parent_start                      = 1,     45,
j_parent_start                      = 1,     45,
parent_grid_ratio                   = 1,     3,     
parent_time_step_ratio              = 1,     3,   
feedback                            = 1,
smooth_option                       = 0
/


&physics
mp_physics                          = 0,     0,   
ra_lw_physics                       = 1,     1,   
ra_sw_physics                       = 1,     1,   
radt                                = 12,    4,   
sf_sfclay_physics                   = 0,     0,   
sf_surface_physics                  = 0,     0,   
bl_pbl_physics                      = 0,     0,   
bldt                                = 0,     0,   
cu_physics                          = 0,     0,   
cudt                                = 0,     0,   
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 5,
sf_urban_physics                    = 0,     0,   
num_land_cat                        = 20,
/

&fdda
/

&dynamics
w_damping                           = 0,
diff_opt                            = 1,
km_opt                              = 4,
diff_6th_opt                        = 0,      0,   
diff_6th_factor                     = 0.12,   0.12,   
base_temp                           = 290.
damp_opt                            = 0,
zdamp                               = 5000.,  5000.,
dampcoef                            = 0.2,    0.2,   
khdif                               = 0,      0,     
kvdif                               = 0,      0,     
non_hydrostatic                     = .true., .true.,
moist_adv_opt                       = 1,      1,      
scalar_adv_opt                      = 1,      1,        
/

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


&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

&datetime
start_date = '2012-07-28_00:00:00',
end_date   = '2012-07-31_00:00:00',
interval_seconds = 3600,
tacc = 0,
debug_level = 0,
/

&io
input_root_name = '/home1/zhy/usr/local/WRF/WRFV3/run/wrfout_d01_2012-07-28_00:00:00'
output_root_name = '/home1/zhy/usr/local/WRF/data/ctl28-31'
plot = 'all_list'
fields = 'height,pressure,tk,tc'
output_type= 'grads'
mercator_defs = .true.
/
split_output = .true.
frames_per_outfile = 1


plot = 'all'
plot = 'list'
plot = 'all_list'
! Below is a list of all available diagnostics
fields = 'height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10,slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr'


&interp
interp_method = 1,
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
/
extrapolate = .false.

interp_method = 1,     ! 0 is model levels, -1 is nice height levels, 1 is user specified pressure/height

interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0,

longlivehj

wrf手册上说：
It must be emphasized that the MODIS-based categories should only be used with the Noah land surface modelin WRF.
需要修改sf_surface_physics参数

不知道是不是这个原因。

sunshinezhpb

近地面层方案、路面过程方案、边界层方案全部没采用...这我还真没试过。
能出结果说明运行没问题，建议你看看WRF指导手册的Selecting Between USGS and MODIS-based Land Use Classifications这一节，然后将陆面过程方案选择为Noah再进行模拟。

schliezer

好的，我再试试。原本是打算不选这些参数化方案，先出个控制组，结果出这问题。

schliezer

这回我把陆面方案改为NOAH，转出来还是一样的问题，降水恒为0，其余要素不随时间变化，真是奇了怪了。