I am a student new to WRF. I want to perform simulations that consider irrigation conditions. I have added the IRRIGATION content in the GEOGRID and downloaded the IRRIGATION data in the static geographical data. Verification shows that the IRRIGATION field has been successfully integrated into the MET file.
I attempted to use the land surface model (LSM) to consider irrigation, setting sf_surface_physics = 2, and I designed related irrigation parameters. However, the results always fail to activate irrigation. In other words, parameters related to irrigation (such as IRR_PH_FIELD) are always zero. It seems that enabling or disabling the irrigation module does not affect the input parameters at all.
Here is my parameter setting:
&time_control
run_days = 1,
run_hours = 18,
run_minutes = 0,
run_seconds = 0,
start_year = 2019,
start_month = 9,
start_day = 4,
start_hour = 6,
end_year = 2019,
end_month = 9,
end_day = 5,
end_hour = 18,
interval_seconds = 21600,
input_from_file = .true.,
history_interval = 3, 3,
frames_per_outfile = 1, 1,
restart = .false.,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
/
&domains
time_step = 60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 600,
e_sn = 600,
e_vert = 35,
p_top_requested = 5000,
num_metgrid_levels = 34,
dx = 10000,
dy = 10000,
grid_id = 1,
parent_id = 0,
i_parent_start = 1,
j_parent_start = 1,
parent_grid_ratio = 1,
parent_time_step_ratio = 1,
feedback = 0,
smooth_option = 0,
/
&physics
sf_surface_physics = 2, ! 使用Noah-MP
mp_physics = 8,
ra_lw_physics = 1,
ra_sw_physics = 1,
sf_sfclay_physics = 1,
bl_pbl_physics = 1,
cu_physics = 1,
isfflx = 1,
num_soil_layers = 4,
num_land_cat = 21,
surface_input_source = 1,
sf_surf_irr_scheme = 3,
irr_daily_amount = 100,
irr_start_hour = 0,
irr_num_hours = 14,
irr_start_julianday = 1,
irr_end_julianday = 364,
irr_freq = 1,
irr_ph = 0,
/
&noah_mp
! opt_irr = 1,
! opt_irrm = 3,
/
&fdda
/
&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0,
base_temp = 290.,
damp_opt = 0,
zdamp = 5000.,
dampcoef = 0.2,
khdif = 0,
kvdif = 0,
non_hydrostatic = .true.,
moist_adv_opt = 1,
scalar_adv_opt = 1,
/
&bdy_control
spec_bdy_width = 5,
specified = .true.,
nested = .false.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
I attempted to use the land surface model (LSM) to consider irrigation, setting sf_surface_physics = 2, and I designed related irrigation parameters. However, the results always fail to activate irrigation. In other words, parameters related to irrigation (such as IRR_PH_FIELD) are always zero. It seems that enabling or disabling the irrigation module does not affect the input parameters at all.
Here is my parameter setting:
&time_control
run_days = 1,
run_hours = 18,
run_minutes = 0,
run_seconds = 0,
start_year = 2019,
start_month = 9,
start_day = 4,
start_hour = 6,
end_year = 2019,
end_month = 9,
end_day = 5,
end_hour = 18,
interval_seconds = 21600,
input_from_file = .true.,
history_interval = 3, 3,
frames_per_outfile = 1, 1,
restart = .false.,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
/
&domains
time_step = 60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 600,
e_sn = 600,
e_vert = 35,
p_top_requested = 5000,
num_metgrid_levels = 34,
dx = 10000,
dy = 10000,
grid_id = 1,
parent_id = 0,
i_parent_start = 1,
j_parent_start = 1,
parent_grid_ratio = 1,
parent_time_step_ratio = 1,
feedback = 0,
smooth_option = 0,
/
&physics
sf_surface_physics = 2, ! 使用Noah-MP
mp_physics = 8,
ra_lw_physics = 1,
ra_sw_physics = 1,
sf_sfclay_physics = 1,
bl_pbl_physics = 1,
cu_physics = 1,
isfflx = 1,
num_soil_layers = 4,
num_land_cat = 21,
surface_input_source = 1,
sf_surf_irr_scheme = 3,
irr_daily_amount = 100,
irr_start_hour = 0,
irr_num_hours = 14,
irr_start_julianday = 1,
irr_end_julianday = 364,
irr_freq = 1,
irr_ph = 0,
/
&noah_mp
! opt_irr = 1,
! opt_irrm = 3,
/
&fdda
/
&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0,
base_temp = 290.,
damp_opt = 0,
zdamp = 5000.,
dampcoef = 0.2,
khdif = 0,
kvdif = 0,
non_hydrostatic = .true.,
moist_adv_opt = 1,
scalar_adv_opt = 1,
/
&bdy_control
spec_bdy_width = 5,
specified = .true.,
nested = .false.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/