新人求助WRF跑NECP2时总是不过

lnpjfb

我是在天河上运行的。用的是1979年的NECP在分析数据2.想做的是全球范围的模拟（毕业论文啊....）
到real.exe这步生成不了
贴出来两个namelist吧。希望高手们帮帮我。

&share
wrf_core = 'ARW',
max_dom = 1,
start_date = '1979-01-01_00:00:00'
end_date   = '1979-12-31_18:00:00'
interval_seconds = 21600
io_form_geogrid = 2,
/
&geogrid
parent_id         =   1,
parent_grid_ratio =   1,   
i_parent_start    =   1,  
j_parent_start    =   1,  
e_we              =  180,
e_sn              =  90,  
geog_data_res     = '10m',
map_proj = 'lat-lon',
stand_lon = 0.
pole_lat = 90.0
pole_lon = 0.0
geog_data_path = '/vol-th/home/libg/WRF/geog'
/
ref_lat = 45.0
ref_lon = -98
dx = 1.0
dy = 1.0
&ungrib
out_format = 'WPS',
prefix = 'FILE',
/
&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/
&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                           = 1,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 1979
start_month                         = 01
start_day                           = 01
start_hour                          = 00
start_minute                        = 00
start_second                        = 00
end_year                            = 1979
end_month                           = 12
end_day                             = 31
end_hour                            = 18
end_minute                          = 00
end_second                          = 00
history_interval                    = 60
interval_seconds                    = 21600
input_from_file                     = .true.
frames_per_outfile                  = 200,
restart                             = .false.,
restart_interval_h                  = 1120,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 0
/
&domains
time_step                           = 600
time_step_fract_num                 = 00
time_step_fract_den                 = 112
max_dom                             = 1,
e_we                                = 180,  
e_sn                                = 90,
e_vert                              = 41,   
p_top_requested                     = 5000,
num_metgrid_levels                  = 18,
num_metgrid_soil_levels             = 0,
dx                                  = 156343.322
dy                                  = 156343.322
p_top_requested                     = 5000
grid_id                             = 1,
parent_id                           = 0,   
i_parent_start                      = 1,   
j_parent_start                      = 1,   
parent_grid_ratio                   = 1,  
parent_time_step_ratio              = 1,  
feedback                            = 1,
smooth_option                       = 00
/
&physics
mp_physics                          = 3,   
ra_lw_physics                       = 1,   
ra_sw_physics                       = 1,   
radt                                = 30,   
sf_sfclay_physics                   = 1,   
sf_surface_physics                  = 1,  
bl_pbl_physics                      = 1,   
bldt                                = 0,   
cu_physics                          = 1,   
cudt                                = 5,   
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 5,
mp_zero_out                         = 0,
/
&fdda
/
&dynamics
diff_opt                            = 0,     
km_opt                              = 0,   
damp_opt                            = 3,
base_temp                           = 290.
zdamp                               = 5000.,
dampcoef                            = 0.2,
khdif                               = 0,     
kvdif                               = 0,      
non_hydrostatic                     = .true.,
moist_adv_opt                       = 0,     
scalar_adv_opt                      = 0,      
chem_adv_opt                        = 0,  
tke_adv_opt                         = 0,     
fft_filter_lat                      = 45.,
w_damping                           = 1,
/
&bdy_control
periodic_x                          = .true.,
symmetric_xs                        = .false.
symmetric_xe                        = .false.
open_xs                             = .false.
open_xe                             = .false.
periodic_y                          = .false.
symmetric_ys                        = .false.
symmetric_ye                        = .false.
open_ys                             = .false.
open_ye                             = .false.
nested                              = .false.
polar                               = .true.
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

lnpjfb

为了大家看着方便，再贴出来rsl.out.0000 吧。
去掉了重复的一些内容。
非常拜托大家了。
Ntasks in X            1, ntasks in Y            1
Linear vertical interpolation is recommended with input vertical resolution this coarse, changing lagrange_order to 1
--- WARNING: traj_opt is zero, but num_traj is not zero; setting num_traj to zero.
--- NOTE: sst_update is 0, setting io_form_auxinput4 = 0 and auxinput4_interval = 0 for all domains
--- NOTE: grid_fdda is 0 for domain      1, setting gfdda interval and ending time to 0 for that domain.
--- NOTE: both grid_sfdda and pxlsm_soil_nudge are 0 for domain      1, setting sgfdda interval and ending time to 0 for that domain.
--- NOTE: obs_nudge_opt is 0 for domain      1, setting obs nudging interval and ending time to 0 for that domain.
--- NOTE: bl_pbl_physics /= 4, implies mfshconv must be 0, resetting
--- NOTE: RRTMG radiation is not used, setting:  o3input=0 to avoid data pre-processing
--- NOTE: num_soil_layers has been set to      5
REAL_EM V3.7.1 PREPROCESSOR
*************************************
Parent domain
ids,ide,jds,jde            1         180           1          90
ims,ime,jms,jme           -4         185          -4          95
ips,ipe,jps,jpe            1         180           1          90
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1,              527717356 bytes allocated
Time period #   1 to process = 1979-01-01_12:00:00.
Time period #   2 to process = 1979-01-01_18:00:00.
Time period #   3 to process = 1979-01-02_00:00:00.
Time period #   4 to process = 1979-01-02_06:00:00.
Time period #   5 to process = 1979-01-02_12:00:00.
Time period #   6 to process = 1979-01-02_18:00:00.
Time period #   7 to process = 1979-01-03_00:00:00.
Time period #   8 to process = 1979-01-03_06:00:00.
Time period #   9 to process = 1979-01-03_12:00:00.
Time period #  10 to process = 1979-01-03_18:00:00.
Time period #  11 to process = 1979-01-04_00:00:00.
Time period #  12 to process = 1979-01-04_06:00:00.
Time period #  13 to process = 1979-01-04_12:00:00.
Time period #  14 to process = 1979-01-04_18:00:00.
Time period #  15 to process = 1979-01-05_00:00:00.
Time period #  16 to process = 1979-01-05_06:00:00.
Time period #  17 to process = 1979-01-05_12:00:00.
Time period #  18 to process = 1979-01-05_18:00:00.
Time period #  19 to process = 1979-01-06_00:00:00.
Total analysis times to input =  235.
  
-----------------------------------------------------------------------------
  
Domain  1: Current date being processed: 1979-01-01_12:00:00.0000, which is loop #   1 out of  235
configflags%julyr, %julday, %gmt:        1979           1   12.00000   
metgrid input_wrf.F first_date_input = 1979-01-01_12:00:00
metgrid input_wrf.F first_date_nml = 1979-01-01_12:00:00
d01 1979-01-01_12:00:00 Timing for input          0 s.
Missing surface temp, replaced with closest level, use_surface set to false.
Missing surface RH, replaced with closest level, use_surface set to false.
Missing surface u wind, replaced with closest level, use_surface set to false.
Missing surface v wind, replaced with closest level, use_surface set to false.
Using sfcprs3 to compute psfc
Converged znw(kte) should be about 0.0 =  -1.2597095E-03
Full level index =    1     Height =     0.0 m
Full level index =    2     Height =    56.6 m      Thickness =   56.6 m
Full level index =    3     Height =   137.9 m      Thickness =   81.4 m
Full level index =    4     Height =   244.7 m      Thickness =  106.8 m
Full level index =    5     Height =   377.6 m      Thickness =  132.9 m
Full level index =    6     Height =   546.3 m      Thickness =  168.7 m
Full level index =    7     Height =   761.1 m      Thickness =  214.8 m
Full level index =    8     Height =  1016.2 m      Thickness =  255.0 m
Full level index =    9     Height =  1307.0 m      Thickness =  290.8 m
Full level index =   10     Height =  1606.4 m      Thickness =  299.5 m
Full level index =   11     Height =  1915.2 m      Thickness =  308.7 m
Full level index =   12     Height =  2233.8 m      Thickness =  318.6 m
Full level index =   13     Height =  2842.5 m      Thickness =  608.7 m
Full level index =   14     Height =  3451.2 m      Thickness =  608.7 m
Full level index =   15     Height =  4059.9 m      Thickness =  608.7 m
Full level index =   16     Height =  4668.5 m      Thickness =  608.7 m
Full level index =   17     Height =  5277.2 m      Thickness =  608.7 m
Full level index =   18     Height =  5885.9 m      Thickness =  608.7 m
Full level index =   19     Height =  6494.6 m      Thickness =  608.7 m
Full level index =   20     Height =  7103.2 m      Thickness =  608.7 m
Full level index =   21     Height =  7711.9 m      Thickness =  608.7 m
Full level index =   22     Height =  8320.6 m      Thickness =  608.7 m
Full level index =   23     Height =  8929.3 m      Thickness =  608.7 m
Full level index =   24     Height =  9537.9 m      Thickness =  608.7 m
Full level index =   25     Height = 10146.6 m      Thickness =  608.7 m
Full level index =   26     Height = 10755.3 m      Thickness =  608.7 m
Full level index =   27     Height = 11364.0 m      Thickness =  608.7 m
Full level index =   28     Height = 11972.6 m      Thickness =  608.7 m
Full level index =   29     Height = 12581.3 m      Thickness =  608.7 m
Full level index =   30     Height = 13190.0 m      Thickness =  608.7 m
Full level index =   31     Height = 13798.7 m      Thickness =  608.7 m
Full level index =   32     Height = 14407.3 m      Thickness =  608.7 m
Full level index =   33     Height = 15016.0 m      Thickness =  608.7 m
Full level index =   34     Height = 15624.7 m      Thickness =  608.7 m
Full level index =   35     Height = 16233.4 m      Thickness =  608.7 m
Full level index =   36     Height = 16842.0 m      Thickness =  608.7 m
Full level index =   37     Height = 17450.7 m      Thickness =  608.7 m
Full level index =   38     Height = 18059.4 m      Thickness =  608.7 m
Full level index =   39     Height = 18668.1 m      Thickness =  608.7 m
Full level index =   40     Height = 19276.7 m      Thickness =  608.7 m
Full level index =   41     Height = 19885.4 m      Thickness =  608.7 m
d01 1979-01-01_12:00:00   flag value of flag_qnwfa is            0
d01 1979-01-01_12:00:00 Old data, no inland lake information
d01 1979-01-01_12:00:00 No soil temperature data for grid%tmn near 30 cm
d01 1979-01-01_12:00:00 Using 1 degree static annual mean temps
d01 1979-01-01_12:00:00 forcing artificial silty clay loam at 1656 points, out of  15931
LAND  CHANGE =            0
WATER CHANGE =            0
d01 1979-01-01_12:00:00 Total pre number of sea ice location changes (water to land) =  10531
error in the grid%tsk
i,j=           1           1
grid%landmask=   1.000000   
grid%tsk, grid%sst, grid%tmn=   3.376384      0.0000000E+00   228.5653   
error in the grid%tsk
i,j=           2           1
grid%landmask=   1.000000   
grid%tsk, grid%sst, grid%tmn=   3.384892      0.0000000E+00   228.5047   
error in the grid%tsk
i,j=           3           1
grid%landmask=   1.000000   
grid%tsk, grid%sst, grid%tmn=   3.392170      0.0000000E+00   228.4555   

weiguang1233

问一个不相干的问题哈，用NCEP2驱动是不是要先写成中间格式呀？

lnpjfb

weiguang1233 发表于 2016-5-18 10:37
问一个不相干的问题哈，用NCEP2驱动是不是要先写成中间格式呀？

不清楚，不过自带的Vtable里有这个数据的

kazami2003

你確定你的ungrib有吃進全部的變數?  我猜你應該有些變數沒有吧....

weiguang1233

lnpjfb 发表于 2016-5-18 16:45
不清楚，不过自带的Vtable里有这个数据的

那直接用Vtable.NCEP2就行了吗？跟fnl相比，要不要做其他改动？

lnpjfb

kazami2003 发表于 2016-5-19 07:03
你確定你的ungrib有吃進全部的變數?  我猜你應該有些變數沒有吧....

确实是有些变量是没有的。比如土壤温度之类的。可是不知道怎么调。

kazami2003

你看看這個教學先吧
http://www2.mmm.ucar.edu/wrf/OnL ... ses/climate_run.php

你應該只有下載氣壓層資料, 地面資料要另外下載flxt06

lnpjfb

kazami2003 发表于 2016-5-19 13:19
你看看這個教學先吧
http://www2.mmm.ucar.edu/wrf/OnL ... ses/climate_run.php

感谢啊。我去学习下。

lnpjfb

kazami2003 发表于 2016-5-19 13:19
你看看這個教學先吧
http://www2.mmm.ucar.edu/wrf/OnL ... ses/climate_run.php

我学习了你发的学习材料，增进了一些基础知识。非常感谢。
想进一步请教：我目前只下载了NCEP/DOE Reanalysis 2 (R2) ds091.0
还需要再下载什么数据吗？