<div>Hi, Everyone,</div> I'm trying to run WRF-LES with defaulted case (no any changes from README.namelist and input_sounding <div>in test/em_les/), but cannot make it successfully. Firstly when running ideal.exe, <div>
I got the NaN in rsl.out.0000 file, through the rsl.error.0000 and rsl.out.0000 ended with 'wrf initilizale successfully",</div><div>and wrfinput_d01 file was generated, with all LU_INDEX values were 0. After changed ztop in namelist.input from </div>
<div>19000 to 10000m, there was no NaN data appeared in rsl.out.0000 file any more (see the following). I'm not sure if this </div><div>kind of change reasonable or not. While when run wrf.exe, I got the message as the following (wether I adjusted ztop</div>
<div> in namelist.input file). BTW, I have set up 'ulimit -s unlimited' in wrf running script, the geo*nc file was provided </div><div> in WPS subdirectory. The namelist.input file was also posted below. </div><div>
<br></div><div> Based on these probelms, I got couple of questions:</div><div>1. How to determine ztop in the namelist.input file? Is it related to "p_top_request" or sounding data provided?</div><div>2. How does the WRF model initilizate land surface variables (LU_INDEX, soil ) in an ideal case? Seems wrf.exe cannot </div>
<div>run because the <span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; "> MMINLU was not defined. Do I need to set up it in code or provide more surface information in namelist.input</span></div>
<div><span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; ">file?But how I'm supposed to do it since as I know almost all the surface related options are for real case? </span></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">3 Is any one willing to post his/her namelist.input file so I can take an example?</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"> Thanks, any of your help is appreciated.</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"> Yan</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"> </span></font></div><div><br></div><div> </div><div> </div><div><br></div><div> I also post the namelist.input file below. If any one </div>
<div>met this kind of problem and how you process this? Can you post your namelist.input file here? </div><div>You are </div><div><br></div><div>rsl.out.0000 for running ideal.exe </div><div><span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; "><div class="im" style="color: rgb(80, 0, 80); ">
<br>returned from reading sounding, nl_in is 41<br></div> ptop is 4755.000 --------------- (is this reasonable ?)<br> base state grid%mub(1,1), p_surf is 95246.47 100001.5<div class="im" style="color: rgb(80, 0, 80); ">
<br> getting moist sounding for full state<br> input sounding surface parameters<br> surface pressure (mb) 1000.000<br> surface pot. temp (K) 305.0000<br> surface mixing ratio (g/kg) 14.00000<br></div><div class="im" style="color: rgb(80, 0, 80); ">
1 0.250E+02 0.300E+03 0.100E+02 0.000E+00 0.000E+00<br> 2 0.750E+02 0.300E+03 0.100E+02 0.000E+00 0.000E+00<br> 3 0.125E+03 0.300E+03 0.100E+02 0.000E+00 0.000E+00<br> 4 0.175E+03 0.300E+03 0.100E+02 0.000E+00 0.000E+00<br>
</div> ...<br> 41 0.205E+03 ....<br> grid%ph_1 calc 0.000000 38.26545 95099.04<br> -147.4297 0.8720049 1.2673497E-02 -27.00002<br> nxc, nyc for perturbation 45 40<div class="im" style="color: rgb(80, 0, 80); ">
<br> delt for perturbation 3.000000<br></div> grid%mu_1 from comp -147.4297<div class="im" style="color: rgb(80, 0, 80); "><br> full state sounding from comp, ph, grid%p, grid%al, grid%t_1, qv<br></div> 1 0.000E+00 9.852E+04 8.843E-01 3.000E+02 1.000E-02<br>
2 3.115E+03 9.496E+04 9.079E-01 3.000E+02 1.000E-02<br> 3 6.313E+03 9.140E+04 9.330E-01 3.000E+02 1.000E-02<br> 4 9.599E+03 8.786E+04 9.769E-01 3.083E+02 4.000E-03<br> 5 1.304E+04 8.432E+04 1.009E+00 3.094E+02 4.000E-03<br>
6 1.659E+04 8.078E+04 1.044E+00 3.105E+02 4.000E-03<br> 7 2.027E+04 7.725E+04 1.082E+00 3.116E+02 4.000E-03<br> 8 2.409E+04 7.371E+04 1.123E+00 3.127E+02 4.000E-03<br> 9 2.804E+04 7.018E+04 1.168E+00 3.139E+02 4.000E-03<br>
10 3.215E+04 6.664E+04 1.216E+00 3.150E+02 4.000E-03<br> 11 3.644E+04 6.310E+04 1.269E+00 3.162E+02 4.000E-03<br> 12 4.091E+04 5.957E+04 1.327E+00 3.173E+02 4.000E-03<br> 13 4.558E+04 5.603E+04 1.391E+00 3.185E+02 4.000E-03<br>
14 5.048E+04 5.249E+04 1.463E+00 3.196E+02 4.000E-03<br> 15 5.563E+04 4.896E+04 1.543E+00 3.207E+02 4.000E-03<br> 16 6.107E+04 4.542E+04 1.634E+00 3.219E+02 4.000E-03<br> 17 6.682E+04 4.189E+04 1.737E+00 3.230E+02 4.000E-03<br>
18 7.294E+04 3.835E+04 1.857E+00 3.242E+02 4.000E-03<br> 19 7.948E+04 3.481E+04 1.996E+00 3.253E+02 4.000E-03<br> 20 8.651E+04 3.128E+04 2.163E+00 3.264E+02 4.000E-03<br> 21 9.413E+04 2.774E+04 2.365E+00 3.276E+02 4.000E-03<br>
22 1.025E+05 2.420E+04 2.616E+00 3.287E+02 4.000E-03<br> 23 1.117E+05 2.067E+04 2.938E+00 3.299E+02 4.000E-03<br> 24 1.220E+05 1.713E+04 3.371E+00 3.310E+02 4.000E-03<br> 25 1.339E+05 1.360E+04 3.990E+00 3.321E+02 4.000E-03<br>
26 1.479E+05 1.006E+04 4.965E+00 3.333E+02 4.000E-03<br> 27 1.654E+05 6.523E+03 6.788E+00 3.344E+02 4.000E-03<div class="im" style="color: rgb(80, 0, 80); "><br> pert state sounding from comp, grid%ph_1, pp, alp, grid%t_1, qv<br>
</div> 1 0.000E+00 2.815E+02 1.234E-02 4.079E-02 1.000E-02<br> 2 3.869E+01 2.517E+02 1.280E-02 4.079E-02 1.000E-02<br> 3 7.887E+01 2.220E+02 1.330E-02 4.079E-02 1.000E-02<br> 4 1.207E+02 2.028E+02 4.921E-03 8.330E+00 4.000E-03<br>
5 1.327E+02 1.941E+02 4.975E-03 9.361E+00 4.000E-03<br> 6 1.447E+02 1.855E+02 5.174E-03 1.047E+01 4.000E-03<br> 7 1.573E+02 1.769E+02 5.385E-03 1.161E+01 4.000E-03<br> 8 1.704E+02 1.683E+02 5.610E-03 1.275E+01 4.000E-03<br>
9 1.840E+02 1.596E+02 5.855E-03 1.389E+01 4.000E-03<br> 10 1.983E+02 1.510E+02 6.121E-03 1.503E+01 4.000E-03<br> 11 2.133E+02 1.424E+02 6.415E-03 1.617E+01 4.000E-03<br> 12 2.290E+02 1.337E+02 6.738E-03 1.731E+01 4.000E-03<br>
13 2.455E+02 1.251E+02 7.096E-03 1.845E+01 4.000E-03<br> 14 2.629E+02 1.165E+02 7.496E-03 1.959E+01 4.000E-03<br> 15 2.814E+02 1.079E+02 7.942E-03 2.073E+01 4.000E-03<br> 16 3.010E+02 9.923E+01 8.458E-03 2.187E+01 4.000E-03<br>
17 3.219E+02 9.060E+01 9.039E-03 2.301E+01 4.000E-03<br> 18 3.443E+02 8.197E+01 9.714E-03 2.415E+01 4.000E-03<br> 19 3.684E+02 7.334E+01 1.052E-02 2.529E+01 4.000E-03<br> 20 3.947E+02 6.471E+01 1.148E-02 2.643E+01 4.000E-03<br>
21 4.233E+02 5.608E+01 1.267E-02 2.757E+01 4.000E-03<br> 22 4.551E+02 4.746E+01 1.416E-02 2.871E+01 4.000E-03<br> 23 4.908E+02 3.883E+01 1.610E-02 2.986E+01 4.000E-03<br> 24 5.315E+02 3.020E+01 1.879E-02 3.099E+01 4.000E-03<br>
25 5.794E+02 2.157E+01 2.283E-02 3.214E+01 4.000E-03<br> 26 6.382E+02 1.294E+01 2.951E-02 3.328E+01 4.000E-03<br> 27 7.152E+02 4.314E+00 4.352E-02 3.442E+01 4.000E-03<div class="im" style="color: rgb(80, 0, 80); ">
<br> wrf: SUCCESS COMPLETE IDEAL INIT<br> <br></div>rsl.out.0000 for wrf.exe<div class="im" style="color: rgb(80, 0, 80); ">
<br> <span class="Apple-style-span" style="color: rgb(0, 0, 0); ">Namelist dfi_control not found in namelist.input. Using registry defaults for variables in dfi_control</span></div><div> Namelist tc not found in namelist.input. Using registry defaults for variables in tc</div>
<div> Namelist scm not found in namelist.input. Using registry defaults for variables in scm</div><div> Namelist fire not found in namelist.input. Using registry defaults for variables in fire</div><div> Ntasks in X 1 , ntasks in Y 1</div>
<div> --- NOTE: sst_update is 0, setting io_form_auxinput4 = 0 and auxinput4_interval = 0 for all domains</div><div> --- NOTE: grid_fdda is 0 for domain 1, setting gfdda interval and ending time to 0 for that domain.</div>
<div> --- NOTE: both grid_sfdda and pxlsm_soil_nudge are 0 for domain 1, setting sgfdda interval and ending time to 0 for that domain.</div><div> --- NOTE: obs_nudge_opt is 0 for domain 1, setting obs nudging interval and ending time to 0 for that domain.</div>
<div> --- NOTE: num_soil_layers has been set to 4</div><div> WRF V3.2.1 MODEL</div><div> *************************************</div><div> Parent domain</div><div> ids,ide,jds,jde 1 80 1 80</div>
<div> ims,ime,jms,jme -4 85 -4 85</div><div> ips,ipe,jps,jpe 1 80 1 80</div><div> *************************************</div><div> DYNAMICS OPTION: Eulerian Mass Coordinate</div>
<div> alloc_space_field: domain 1 , 136632968 bytes allocated</div><div> med_initialdata_input: calling input_input</div><div> MMINLU error on input</div><div> INITIALIZE THREE Noah LSM RELATED TABLES</div>
<div> Skipping over LUTYPE = USGS</div><div> Skipping over LUTYPE = MODIFIED_IGBP_MODIS_NOAH</div><div> Skipping over LUTYPE = USGS-RUC</div><div> -------------- FATAL CALLED ---------------</div><div> FATAL CALLED FROM FILE: <stdin> LINE: 1547</div>
<div> Land Use Dataset '' not found in VEGPARM.TBL.</div><div><br></div><div class="im" style="color: rgb(80, 0, 80); "> namelist.input I used </div><div class="im" style="color: rgb(80, 0, 80); "> </div><div class="im" style="color: rgb(80, 0, 80); ">
<div class="im"> &time_control</div><div class="im"> run_days = 1,</div><div class="im"> run_hours = 12,</div><div class="im"> run_minutes = 0,</div>
<div class="im"> run_seconds = 0,</div><div class="im"> start_year = 2008, 2008, 2008,</div><div class="im"> start_month = 06, 01, 01,</div><div class="im">
start_day = 06, 24, 24,</div><div class="im"> start_hour = 00, 12, 12,</div><div class="im"> start_minute = 00, 00, 00,</div><div class="im">
start_second = 00, 00, 00,</div><div class="im"> end_year = 2008, 2008, 2008,</div><div class="im"> end_month = 06, 01, 01,</div><div class="im">
end_day = 07, 25, 25,</div><div class="im"> end_hour = 00, 12, 12,</div><div class="im"> end_minute = 00, 00, 00,</div><div class="im">
end_second = 00, 00, 00,</div><div class="im"> interval_seconds = 10800</div><div class="im"> input_from_file = .false.,.true.,.true.,</div><div class="im">
history_interval = 60, 60, 60,</div><div class="im"> frames_per_outfile = 1, 1000, 1000,</div><div class="im"> restart = .false.,</div><div class="im"> restart_interval = 5000,</div>
<div class="im"> io_form_history = 2</div><div class="im"> io_form_restart = 2</div><div class="im"> io_form_input = 2</div><div class="im"> io_form_boundary = 2</div>
<div class="im"> debug_level = 0</div><div class="im"> /</div><div class="im"><br></div><div class="im"> &domains</div><div class="im"> time_step = 5,</div><div class="im">
time_step_fract_num = 0,</div><div class="im"> time_step_fract_den = 1,</div><div class="im"> max_dom = 1,</div><div class="im"> e_we = 80, 112, 94,</div>
<div class="im"> e_sn = 80, 97, 91,</div><div class="im"> e_vert = 29, 28, 28,</div><div class="im"> num_metgrid_levels = 30,</div><div class="im">
num_metgrid_soil_levels = 4,</div><div class="im"> dx = 100, 10000, 3333.33,</div><div class="im"> dy = 100, 10000, 3333.33,</div><div class="im">
grid_id = 1, 2, 3,</div><div class="im"> parent_id = 0, 1, 2,</div><div class="im"> i_parent_start = 1, 31, 30,</div><div class="im">
j_parent_start = 1, 17, 30,</div><div class="im"> parent_grid_ratio = 1, 3, 3,</div><div class="im"> parent_time_step_ratio = 1, 3, 3,</div><div class="im">
feedback = 1,</div><div class="im"> smooth_option = 0</div><div class="im"> ztop = 10000.</div><div class="im"> </div><div class="im"> interp_type = 2 </div>
<div class="im"> extrap_type = 2 </div><div class="im"> t_extrap_type = 2 </div><div class="im"> use_levels_below_ground = .true.</div><div class="im">
use_surface = .true.</div><div class="im"> lagrange_order = 1 </div><div class="im"> zap_close_levels = 500 </div><div class="im"> lowest_lev_from_sfc = .false. </div>
<div class="im"> force_sfc_in_vinterp = 1 </div><div class="im"> sfcp_to_sfcp = .false. </div><div class="im"> smooth_cg_topo<span class="Apple-tab-span" style="white-space:pre"> </span> = .false. </div>
<div class="im"> use_tavg_for_tsk = .false. </div><div class="im"> vert_refine_fact = 1 </div><div class="im"> /</div><div class="im"><br></div><div class="im"> &physics</div>
<div class="im"> mp_physics = 3, 3, 3,</div><div class="im"> ra_lw_physics = 1, 1, 1,</div><div class="im"> ra_sw_physics = 1, 1, 1,</div>
<div class="im"> radt = 30, 30, 30,</div><div class="im"> sf_sfclay_physics = 1, 1, 1,</div><div class="im"> sf_surface_physics = 2, 2, 2,</div>
<div class="im"> bl_pbl_physics = 0, 1, 1,</div><div class="im"> bldt = 0, 0, 0,</div><div class="im"> cu_physics = 1, 1, 0,</div>
<div class="im"> cudt = 5, 5, 5,</div><div class="im"> isfflx = 1,</div><div class="im"> ifsnow = 0,</div><div class="im"> icloud = 1,</div>
<div class="im"> surface_input_source = 1,</div><div class="im"> num_soil_layers = 4,</div><div class="im"> sf_urban_physics = 0, 0, 0,</div><div class="im"> maxiens = 1,</div>
<div class="im"> maxens = 3,</div><div class="im"> maxens2 = 3,</div><div class="im"> maxens3 = 16,</div><div class="im"> ensdim = 144,</div>
<div class="im"> /</div><div class="im"><br></div><div class="im"> &fdda</div><div class="im"> /</div><div class="im"><br></div><div class="im"> &dynamics</div><div class="im"> w_damping = 0,</div>
<div class="im"> diff_opt = 2,</div><div class="im"> km_opt = 2,</div><div class="im"> diff_6th_opt = 0, 0, 0,</div><div class="im">
diff_6th_factor = 0.12, 0.12, 0.12,</div><div class="im"> base_temp = 290.</div><div class="im"> damp_opt = 0,</div><div class="im"> zdamp = 5000., 5000., 5000.,</div>
<div class="im"> dampcoef = 0.2, 0.2, 0.2</div><div class="im"> khdif = 0, 0, 0,</div><div class="im"> kvdif = 0, 0, 0,</div>
<div class="im"> non_hydrostatic = .true., .true., .true.,</div><div class="im"> moist_adv_opt = 1, 1, 1, </div><div class="im"> scalar_adv_opt = 1, 1, 1, </div>
<div class="im"> pert_coriolis = .true.,</div><div class="im"> do_coriolis = .true.,</div><div class="im"> do_curvature = .true., </div><div class="im">
do_gradp = .true., </div><div class="im"> /</div><div class="im"><br></div><div class="im"> &bdy_control</div><div class="im"> spec_bdy_width = 5,</div><div class="im">
spec_zone = 1,</div><div class="im"> relax_zone = 4,</div><div class="im"> specified = .false., .false.,.false.,</div><div class="im"> periodic_x = .true.,.false.,.false.,</div>
<div class="im"> symmetric_xs = .false.,.false.,.false.,</div><div class="im"> symmetric_xe = .false.,.false.,.false.,</div><div class="im"> open_xs = .false.,.false.,.false.,</div>
<div class="im"> open_xe = .false.,.false.,.false.,</div><div class="im"> periodic_y = .true.,.false.,.false.,</div><div class="im"> symmetric_ys = .false.,.false.,.false.,</div>
<div class="im"> symmetric_ye = .false.,.false.,.false.,</div><div class="im"> open_ys = .false.,.false.,.false.,</div><div class="im"> open_ye = .false,.false.,.false.,</div>
<div class="im"> nested = .false., .true., .true.,</div><div class="im"> /</div><div class="im"><br></div><div class="im"> &grib2</div><div class="im"> /</div><div class="im"><br></div><div class="im">
&namelist_quilt</div><div class="im"> nio_tasks_per_group = 0,</div><div class="im"> nio_groups = 1,</div><div class="im"> /</div><div><br></div></div><div class="im" style="color: rgb(80, 0, 80); "> <span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: arial; font-size: small; ">-- </span></div>
</span></div><div><div>Yan <br><br>
</div></div></div>