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