<DIV>
<DIV>Hi:</DIV>
<DIV> Recently I use gen_be to generate my background error statistics .My WPS,WRF and WRFDA are all 3.2 vertion,and I use SGI 128 CPU computer ,the compiler is ifort .I have tested gen_be testdata,and it works well.But when I use my own data,the problem comes.After stage 4,the program shows FAIL.Then I view stage0~stage4 .log file,and I find the gen_be_stage2.log shows that :</DIV>
<DIV>[3] Calculate eigenvectors, eigenvalues and inverse for psi/psi covariance <BR> Mode Eigenvalue Cumulative Variance e(k,k)<BR> 1 0.0000E+00 NaN 0.0000E+00<BR> 2 0.0000E+00 NaN 0.0000E+00<BR> 3 0.0000E+00 NaN 0.0000E+00<BR> 4 0.0000E+00 NaN 0.0000E+00<BR> 5 0.0000E+00 NaN 0.0000E+00<BR> 6 0.0000E+00 NaN 0.0000E+00<BR> 7 0.0000E+00 NaN 0.0000E+00<BR> There is no value in the log.But when I run the test case ,there is diffrent value in the log.</DIV>
<DIV> So I think the gen_be is OK,the problem maybe comes from my wrfout. My classmate uses WRF3.1 to generate one month wrfout ,then uses gen_be to generate BE,it is OK.So I also use WRF3.1 to generate one month wrfout ,but when I use a few days wrfout to test, the problem is still as above.</DIV>
<DIV> I don't know how to solve this problem.Please help me.</DIV>
<DIV> </DIV>
<DIV> My WRF3.1 namelist.input is:</DIV>
<DIV> </DIV>
<DIV> &time_control<BR> run_days = 30,<BR> run_hours = 0,<BR> run_minutes = 0,<BR> run_seconds = 0,<BR> start_year = 2003, 2003, 2003,<BR> start_month = 07, 07, 07,<BR> start_day = 01, 01, 01,<BR> start_hour = 12, 12, 12,<BR> start_minute = 00, 00, 00,<BR> start_second = 00, 00, 00,<BR> end_year = 2003, 2003, 2003,<BR> end_month = 07, 07, 07,<BR> end_day = 31, 31, 31,<BR> end_hour = 12, 12, 12,<BR> end_minute = 00, 00, 00,<BR> end_second = 00, 00, 00,<BR> interval_seconds = 21600<BR> input_from_file = .true.,.true.,.true.,<BR> history_interval = 720, 720, 720,<BR> frames_per_outfile = 1, 1, 1,<BR> restart = .false.,<BR> restart_interval = 1440,<BR> io_form_history = 2<BR> io_form_restart = 2<BR> io_form_input = 2<BR> io_form_boundary = 2<BR> debug_level = 0<BR> /</DIV>
<DIV> &domains<BR> time_step = 144,<BR> time_step_fract_num = 0,<BR> time_step_fract_den = 1,<BR> max_dom = 3,<BR> e_we = 160, 259, 256,<BR> e_sn = 140, 202, 226,<BR> e_vert = 28, 28, 28,<BR> p_top_requested = 5000,<BR> num_metgrid_levels = 27,<BR> num_metgrid_soil_levels = 4,<BR> dx = 45000, 15000, 5000,<BR> dy = 45000, 15000, 5000,<BR> grid_id = 1, 2, 3,<BR> parent_id = 1, 1, 2,<BR> i_parent_start = 1, 38, 88,<BR> j_parent_start = 1, 38, 63,<BR> parent_grid_ratio = 1, 3, 3,<BR> parent_time_step_ratio = 1, 3, 3,<BR> feedback = 0,<BR> smooth_option = 0<BR> /</DIV>
<DIV> &physics<BR> mp_physics = 2, 2, 2,<BR> ra_lw_physics = 1, 1, 1,<BR> ra_sw_physics = 1, 1, 1,<BR> radt = 30, 30, 30,<BR> sf_sfclay_physics = 1, 1, 1,<BR> sf_surface_physics = 1, 1, 1,<BR> bl_pbl_physics = 1, 1, 1,<BR> bldt = 0, 0, 0,<BR> cu_physics = 1, 1, 0,<BR> cudt = 5, 5, 5,<BR> isfflx = 1,<BR> ifsnow = 0,<BR> icloud = 1,<BR> surface_input_source = 1,<BR> num_soil_layers = 5,<BR> sf_urban_physics = 0, 0, 0,<BR> maxiens = 1,<BR> maxens = 3,<BR> maxens2 = 3,<BR> maxens3 = 16,<BR> ensdim = 144,<BR> /</DIV>
<DIV> &fdda<BR> /</DIV>
<DIV> &dynamics<BR> w_damping = 0,<BR> diff_opt = 1,<BR> km_opt = 4,<BR> diff_6th_opt = 0, 0, 0,<BR> diff_6th_factor = 0.12, 0.12, 0.12,<BR> base_temp = 290.<BR> damp_opt = 0,<BR> zdamp = 5000., 5000., 5000.,<BR> dampcoef = 0.2, 0.2, 0.2<BR> khdif = 0, 0, 0,<BR> kvdif = 0, 0, 0,<BR> non_hydrostatic = .true., .true., .true.,<BR> moist_adv_opt = 1, 1, 1, <BR> scalar_adv_opt = 1, 1, 1, <BR> /</DIV>
<DIV> &bdy_control<BR> spec_bdy_width = 5,<BR> spec_zone = 1,<BR> relax_zone = 4,<BR> specified = .true., .false.,.false.,<BR> nested = .false., .true., .true.,<BR> /</DIV>
<DIV> &grib2<BR> /</DIV>
<DIV> &namelist_quilt<BR> nio_tasks_per_group = 0,<BR> nio_groups = 1,<BR> /<BR></DIV>
<DIV> I think my codes are all OK,maybe the old fnl files (my case uses 2003 GRIB1 fnl files) cause the problem?</DIV>
<DIV> I look forward to receiving your reply. </DIV>
<DIV> </DIV>
<DIV> liuyin </DIV></DIV><br><br><span title="neteasefooter"><span id="netease_mail_footer"><hr/>
<a href="http://yxp.163.com/photo/ep.html?sss=fromyx0911" target="_blank">全国最低价,天天在家冲照片,24小时发货上门!</a>
</span></span>