ARWpost转换WRFout后时间层次只有1个，求教解决方法啊。。。纠结好多天了

Dangozsd

最近刚开始学习WRF模式，在用ARWpost对wrfout进行处理后，生成的ctl一直只有一个时次，我做的是运行1天的，但在GrADS里运行set t 2什么的都显示entire grid undefined，不知道是哪里出了问题，求教各位大大帮助啊= =
以下是各namelist的内容，我的是单层区域，模拟的区域网格数据是用WRFDomainWizard生成更改的，不知道会不会是这里出了问题
namelist.wps
&share
wrf_core = 'ARW',
max_dom = 1,
start_date = '2014-06-06_00:00:00','2014-06-06_00:00:00',
end_date   = '2014-06-07_00:00:00','2014-06-07_00:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/

&geogrid
parent_id         =   1,   1,
parent_grid_ratio =   1,   3,
i_parent_start    =   1,  21,
j_parent_start    =   1,  14,
e_we              =  75, 178,
e_sn              =  73, 169,
geog_data_res     = '10m','2m',
dx = 17000,
dy = 17000,
map_proj = 'lambert',
ref_lat   =  25.0,
ref_lon   = 102.0,
truelat1  =  30.0,
truelat2  =  60.0,
stand_lon =102.0,
geog_data_path = '/nuist/p/public/data/geog'
/

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

&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/
namelist.inpout
&time_control
run_days                            = 0,
run_hours                           = 24,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2014, 2014, 2000,
start_month                         = 06,   06,   01,
start_day                           = 06,   06,   24,
start_hour                          = 00,   00,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2014, 2014, 2000,
end_month                           = 06,   06,   01,
end_day                             = 07,   07,   25,
end_hour                            = 00,   00,   12,
end_minute                          = 00,   00,   00,
end_second                          = 00,   00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 180,  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                           = 102,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 75,    178,   94,
e_sn                                = 73,    169,    91,
e_vert                              = 28,    28,    30,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 17000, 6400,  3333.33,
dy                                  = 17000, 6400,  3333.33,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     29,    30,
j_parent_start                      = 1,     25,    30,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
max_dz                              = 1050
/

&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                  = 2,     2,     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,
sf_urban_physics                    = 0,     0,     0,
/

&fdda
/

&dynamics
w_damping                           = 0,
diff_opt                            = 1,
km_opt                              = 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,
/
namelist.ARWpost
&datetime
start_date = '2014-06-06_00:00:00',
end_date   = '2014-06-07_00:00:00',
interval_seconds = 21600,
tacc = 0,
debug_level = 0,
/

&io
input_root_name = '/nuist/u/home/zsd/WRFV3/test/em_real/wrfout_d01_2014-06-06_00:00:00'
output_root_name = './test'
plot = 'all_list'
fields = 'height,pressure,tk,tc'
mercator_defs = .true.
/
split_output = .true.
frames_per_outfile = 2


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 = 0,
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
/
extrapolate = .true.

interp_method = 0,     ! 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,

Mid_Farmer

如果你的ctl里tdef是4的话 那可能是你input里的timestep=102设置导致的输出arwpost的21600输出不出来

Dangozsd

Mid_Farmer 发表于 2015-4-25 09:37
如果你的ctl里tdef是4的话 那可能是你input里的timestep=102设置导致的输出arwpost的21600输出不出来

我的ctl里tdef都是1= =那么这该怎么更改呢？加大input里的timestep值？

lanlaniris

恩,感觉是有点问题……不过我也是初学，不知道能不能帮到你哈！
1、你区域设置的时候e_we 是75，一般来说要遵循减1后能被3整除的规律，你看用76吧
2、&physics里面radt设置的一般是dx的公里数，你的应该是17吧
3、你的WRF里面设的你的文件输出的时间间隔（history_interval）是3小时（180min），那么相应作DAT资料的时候应该把ARWpost的interval_seconds设成3小时的（10800），这样才不会乱。我以前也有设错过，导致子网格和母网格数据混在一起了。
你试着改一下，不知道有没有用哈~

Dangozsd

lanlaniris 发表于 2015-4-25 10:00
恩,感觉是有点问题……不过我也是初学，不知道能不能帮到你哈！
1、你区域设置的时候e_we 是75，一般来说 ...

谢谢回复，我现在去试试~

渭水听风

进来学习一下

lhaikun@163.com

lanlaniris 发表于 2015-4-25 10:00
恩,感觉是有点问题……不过我也是初学，不知道能不能帮到你哈！
1、你区域设置的时候e_we 是75，一般来说 ...

减1被三整除指的是内层嵌套的网格数

hwl1943

你的namelist.input设定WRF是180分输出一次结果，但是计算的time_step=102秒，180*60/102根本除不尽，所以实际输出的时间并不是在你所想的时间点上，利用ARWpost转档当然时间点只有2014/6/6/00这笔时间的数据。

lhaikun@163.com

请问楼主解决问题了吗，我的tdef为27 满足整除关系，还是没有。。。

------------------------------------------------------------------------------------------
已经解决了我的，楼主要是还没可以去我那溜溜

熊猫支队007

lhaikun@163.com 发表于 2015-7-5 14:27
请问楼主解决问题了吗，我的tdef为27 满足整除关系，还是没有。。。

--------------------------------- ...

请问下解决了么？想问下解决方案。遇到同样的问题，纠结了很久