Wind farm drag parameterization in nest domain

Topics related to current and future physics in the WRF as well as any problems you may have.

Wind farm drag parameterization in nest domain

Postby luochen » Tue Apr 14, 2015 2:55 am

hello everyone,
I just have a question about the wind farm parameterization in WRF(ARW). I use nested domains to simulate the influence of a wind farm. but the influence (the change of pbl for example) only in the coarse grid, and the fine grid do not have any change. it seems that the windturbine parameter only set in the coarse grid. the following is my namelist.input , I wonder if there is something wrong in my namelist.

***************************************************************************************************
&time_control
run_days = 101,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2011, 2011, 2011,
start_month = 11, 11, 11,
start_day = 21, 21, 21,
start_hour = 0, 0, 0,
start_minute = 00, 00, 00,
start_second = 00, 00, 0,
end_year = 2012, 2012, 2012,
end_month = 03, 03, 1,
end_day = 01, 01, 30,
end_hour = 00, 00, 18,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600,
input_from_file = .true., .true., .true.,
history_interval = 360, 360, 60,
frames_per_outfile = 99999, 999999, 999999,
restart = .false.,
restart_interval = 99999,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
/

&domains
eta_levels = 1.000, 0.998, 0.996, 0.994, 0.992,
0.99, 0.988, 0.986, 0.984, 0.981,
0.975, 0.967, 0.957, 0.946, 0.932,
0.917, 0.899, 0.8825, 0.86, 0.835,
0.81, 0.78, 0.75, 0.71, 0.67,
0.625, 0.585, 0.54, 0.5, 0.46,
0.42, 0.38, 0.344, 0.314, 0.288,
0.263, 0.242, 0.22, 0.2, 0.1807,
0.1602, 0.1434, 0.1255, 0.1055, 0.087,
0.0671, 0.0484, 0.03, 0.0145, 0.000
,
time_step = 150,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 70, 96, 96,
e_sn = 60, 76, 86,
e_vert = 50, 50, 50,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 25000, 5000, 1000,
dy = 25000, 5000, 1000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 24, 39,
j_parent_start = 1, 22, 30,
parent_grid_ratio = 1, 5, 5,
parent_time_step_ratio = 1, 5, 5,
feedback = 1,
smooth_option = 0,
/

&physics
mp_physics = 6, 6, 6,
ra_lw_physics = 3, 3, 3,
ra_sw_physics = 3, 3, 3,
radt = 1, 1, 1,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 5, 5, 5,
bldt = 0, 0, 0,
cu_physics = 1, 1, 0,
cudt = 0, 0, 0,
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,
windfarm_opt = 1,
windfarm_ij = 0,

/

&fdda
/

&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
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,
/

&bdy_control
spec_bdy_width = 7,
spec_zone = 1,
relax_zone = 6,
specified = .true., .false., .false.,
nested = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
******************************************************************************************************
thanks,
luochen
luochen
 
Posts: 1
Joined: Mon Apr 13, 2015 11:40 pm

Return to Physics

Who is online

Users browsing this forum: No registered users and 4 guests