【已解决】请教！关于99年以前利用wrf模拟的数据问题~~~~

鲸鱼便子

请教各位~~~
现在想模拟99年以前的降水个例，查阅了一下fnl数据http://rda.ucar.edu/datasets/ds082.0/index.html?hash=sfol-wl-/data/ds082.0#!description，发现是半月一个文件，而且文件名是Y****的形势。
请教一下这个grib数据能不能用来转WRF，因为我试了一下ungrib的时候就出问题了~~~~
希望弄过的高手帮忙解答一下吧~~~~各路大神帮忙看看吧~~@超级版主 @管理员
谢啦~~~


解决：
一般99年以前不用这种资料做初始场，一般用ds090
可以参考http://www.mmm.ucar.edu/wrf/OnLineTutorial/DATA/NNRP/index.html
和http://www.mnmuc.org/viewthread.php?page=1&tid=1296

鲸鱼便子

sunshinezhpb 发表于 2014-3-10 22:36

                               
登录/注册后可看大图

感谢LZ提供的经验。
但是我前处理完以后碰到了一些问题。在我real.exe时出错，rsl.out里提示我说它这个资料 ...

你查一下wrf manual里面，我记得可以从别的资料里面提取土壤资料，然后met到一起。。。。
要是找不到土壤资料那我也不知道该怎么办了~~~

sunshinezhpb

感谢LZ提供的经验。
但是我前处理完以后碰到了一些问题。在我real.exe时出错，rsl.out里提示我说它这个资料里不包含土壤温度等信息，所以Noah陆面方案不能用。但我又得选用城市冠层模型，不知道LZ有什么好的建议吗?

sunshinezhpb

鲸鱼便子 发表于 2014-3-11 08:20

                               
登录/注册后可看大图

你查一下wrf manual里面，我记得可以从别的资料里面提取土壤资料，然后met到一起。。。。
要是找不到土壤 ...

别的资料是指？我在ds090的下载界面看了下，貌似没发现别的地面资料...
如果是指其它的再分析资料的话，我还是考虑直接换的别的资料做初始场算了...唉，真恼火

鲸鱼便子

sunshinezhpb 发表于 2014-3-11 20:52

                               
登录/注册后可看大图

别的资料是指？我在ds090的下载界面看了下，貌似没发现别的地面资料...
如果是指其它的再分析资料的话， ...

我昨天看了一下我的namelist.input(成功转出结果了的)，陆面方案就是Noah呀，没有出错的说~~~
只是这个资料里面 num_metgrid_soil_levels             = 2~~~
你再看看，你是不是哪里弄得有bug~~~

sunshinezhpb

鲸鱼便子 发表于 2014-3-12 08:39

                               
登录/注册后可看大图

我昨天看了一下我的namelist.input(成功转出结果了的)，陆面方案就是Noah呀，没有出错的说~~~
只是这个资 ...

这是rsl.out里的错误信息：
----------------- ERROR -------------------
namelist    : num_metgrid_soil_levels =          2
input files : NUM_METGRID_SOIL_LEVELS =          0 (from met_em files).
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:     677
Mismatch between namelist and global attribute NUM_METGRID_SOIL_LEVELS
-----------------------w--------------------
rface v wind, replaced ith closest leel, use_surface set to false.
Using sfcprs  to compute psfc
Converged znw(kte) should be about 0.0v =   -2.0489097E-08
d01 1980-01-01_00:00:00 Old data, no inland lake information
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    1596
Not enough soil temperature data for Noah LSM scheme.

这是我的namelist.input
&time_control            
run_days                 = 0,
run_hours                = 720,
run_minutes              = 0,
run_seconds              = 0,
start_year               = 1980,     1980,
start_month              = 01,       01,
start_day                = 1,        1,
start_hour               = 0,        0,
start_minute             = 00,       00,
start_second             = 00,       00,
end_year                 = 1980,     1980,
end_month                = 01,       01,
end_day                  = 31,       31,
end_hour                 = 0,        0,
end_minute               = 00,       00,
end_second               = 00,       00,
interval_seconds         = 21600,
input_from_file          = .true.,   .true.,
history_interval         = 1440,     480,
frames_per_outfile       = 1000,     1000,
restart                  = .false.,
restart_interval         = 50000,
io_form_history          = 2,
io_form_restart          = 2,
io_form_input            = 2,
io_form_boundary         = 2,
debug_level              = 0,
/

&domains                 
time_step                = 90,
time_step_fract_num      = 0,
time_step_fract_den      = 1,
max_dom                  = 2,
e_we                     = 100,       88,
e_sn                     = 100,       88,
e_vert                   = 28,       28,
p_top_requested          = 5000,
num_metgrid_levels       = 18,
num_metgrid_soil_levels  = 2,
dx                       = 15000,     5000,
dy                       = 15000,     5000,
grid_id                  = 1,        2,
parent_id                = 1,        1,
i_parent_start           = 1,       36,
j_parent_start           = 1,       36,
parent_grid_ratio        = 1,        3,
parent_time_step_ratio   = 1,        3,
feedback                 = 1,
smooth_option            = 0,
/

&physics                 
mp_physics               = 3,        3,
ra_lw_physics            = 1,        1,
ra_sw_physics            = 1,        1,
radt                     = 30,       30,
sf_sfclay_physics        = 1,        1,
sf_surface_physics       = 2,        2,
bl_pbl_physics           = 1,        1,
bldt                     = 0,        0,
cu_physics               = 1,        1,
cudt                     = 5,        5,
isfflx                   = 1,
ifsnow                   = 0,
icloud                   = 1,
surface_input_source     = 1,
num_soil_layers          = 2,
sf_urban_physics         = 1,        1,
maxiens                  = 1,
maxens                   = 3,
maxens2                  = 3,
maxens3                  = 16,
ensdim                   = 144,
/

&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           = 10,
spec_zone                = 1,
relax_zone               = 9,
specified                = .true.,  .false.,
nested                   = .false.,   .true.,
/

&grib2                  
/

&namelist_quilt         
nio_tasks_per_group      = 0,
nio_groups               = 1,
/
                                                                                                                                                 还是提示不对呀...

sunshinezhpb

鲸鱼便子 发表于 2014-3-12 08:39

                               
登录/注册后可看大图

我昨天看了一下我的namelist.input(成功转出结果了的)，陆面方案就是Noah呀，没有出错的说~~~
只是这个资 ...

下载的资料是grb2d和pgb.f00，通过ungrib解压成NNRP和NNRPSFC后用metgrid合到一起，这个它没报错，应该是没问题的吧..

鲸鱼便子

sunshinezhpb 发表于 2014-3-12 17:08

                               
登录/注册后可看大图

下载的资料是grb2d和pgb.f00，通过ungrib解压成NNRP和NNRPSFC后用metgrid合到一起，这个它没报错，应该是 ...

看你这两个都木有神马问题呀~~~
估计还是ungrib的时候出问题了，你的Vtable选的是什么？
这个资料Vtable要选Vtable.NNRP~~

鲸鱼便子

sunshinezhpb 发表于 2014-3-13 11:52

                               
登录/注册后可看大图

感谢你的细心解答。Vtable链接的是Vtable.NNRP没错，问题应该是在前处理上，不知道你能否将你运行成功的n ...

嘿嘿，实在不好意思，我偷懒用的wrfdomain，但是还是把namelist.wps帖上来吧，仅供参考~~~

&share
wrf_core = 'ARW',
max_dom = 2,
start_date = '2006-07-27_18:00:00', '2006-07-27_18:00:00',
end_date   = '2006-07-29_00:00:00', '2006-07-29_00:00:00',
interval_seconds = 21600,
io_form_geogrid = 2,
opt_output_from_geogrid_path = '********',
debug_level = 0,
/

&geogrid
parent_id         = 1,1,
parent_grid_ratio = 1,3,
i_parent_start    = 1,33,
j_parent_start    = 1,36,
e_we          = 95,118,
e_sn          = 85,109,
geog_data_res = '10m','5m',
dx = 30000,
dy = 30000,
map_proj =  'lambert',
ref_lat   = 31,
ref_lon   = 112,
truelat1  = 30,
truelat2  = 60,
stand_lon = 112,
geog_data_path = '/***********',
opt_geogrid_tbl_path = '*****************',
ref_x = 47.5,
ref_y = 42.5,
/

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

&metgrid
fg_name = 'FILE',
io_form_metgrid = 2,
opt_output_from_metgrid_path = '/**************',
opt_metgrid_tbl_path = '/*****************',
/

&mod_levs
  press_pa = 201300 , 200100 , 100000 ,
             95000 ,  90000 ,
             85000 ,  80000 ,
             75001 ,  70000 ,
             65000 ,  60000 ,
             55000 ,  50000 ,
             45000 ,  40000 ,
             35000 ,  30000 ,
             25000 ,  20000 ,
             15000 ,  10000 ,
              5000 ,   1000
/


&domain_wizard
grib_data_path = '/******************',
grib_vtable = 'Vtable.NNRP',
dwiz_name    =typhoon_henan_8808
dwiz_desc    =
dwiz_user_rect_x1 =1563
dwiz_user_rect_y1 =264
dwiz_user_rect_x2 =1724
dwiz_user_rect_y2 =401
dwiz_show_political =true
dwiz_center_over_gmt =true
dwiz_latlon_space_in_deg =10
dwiz_latlon_linecolor =-8355712
dwiz_map_scale_pct =12.5
dwiz_map_vert_scrollbar_pos =0
dwiz_map_horiz_scrollbar_pos =0
dwiz_gridpt_dist_km =30.0
dwiz_mpi_command =
dwiz_tcvitals =null
dwiz_bigmap =Y
/

sunshinezhpb

鲸鱼便子 发表于 2014-3-13 08:25

                               
登录/注册后可看大图

看你这两个都木有神马问题呀~~~
估计还是ungrib的时候出问题了，你的Vtable选的是什么？
这个资料Vtabl ...

感谢你的细心解答。Vtable链接的是Vtable.NNRP没错，问题应该是在前处理上，不知道你能否将你运行成功的namelist.wps贴上来呢？