WRF4.2 irrigation scheme is bad

[wangfengke]

Dear friends:
I want to simulate from June 1, 2018 to September 30, 2018, and put the downloaded irrigation data (irrigation/fao/00001-04320.00001-02160 and index downloaded from the official website) into the static terrain data geog in the WRF4.2 version , and then namelist.input input
sf_surface_physics=2,choose Noah,
and
sf_surf_irr_scheme = 0, 1
irr_daily_amount = 0, 8
irr_start_hour = 0, 14
irr_num_hours = 0, 2
irr_start_julianday = 0, 121
irr_end_julianday = 0, 170
irr_ph = 0, 0
irr_freq = 0, 3
after running geogrid.exe, the percentage of irrigation area in the area can be drawn through geo_em.nc, but after wrf.exe, it is found that irrigation does not affect the results at all, what is going on? Need to modify the index? Or what? Many thanks!

 

[kwerner]

Hi,
This could possibly be due to your setting of "surface_input_source" in namelist.input. It should be set to 3. If not, try again with that setting. Otherwise, another thought is that the values you are putting in are being overwritten by the VEGPARM.TBL, which is what the Noah scheme uses. You may want to take a look at that and see if anything stands out to you, and you could perhaps modify the *.TBL. If that's still not helpful, please attach your full namelist.input file so I can take a look. Thanks!

 

[wangfengke]

Hello, I have remade the irrigation land percentage to replace the WRF default , and then I can also draw the irrigation land percentage through the geo_em.d02.nc, which means that I have successfully replaced, but the problem is coming. Replace the irrigation land percentage, but the temperature output by the model, why is there still no change? namelist.input has been unchanged.Attach my namelist.input.
&time_control
run_days = 127,
run_hours = 18,
run_minutes = 0,
run_seconds = 0,
start_year = 2016, 2016,
start_month = 05, 05,
start_day = 26, 26,
start_hour = 00, 00,
end_year = 2016, 2016,
end_month = 09, 09,
end_day = 30, 30,
end_hour = 18, 18,
interval_seconds = 21600
input_from_file = .true.,.true.,
history_interval = 360, 360,
frames_per_outfile = 4, 4,
restart = .false.,
restart_interval = 10000,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
auxinput4_inname ="wrflowinp_d<domain>"
auxinput4_interval = 360, 360,
io_form_auxinput4 = 2
/

&domains
time_step = 150,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 40, 111,
e_sn = 40, 121,
e_vert = 33, 33,
p_top_requested = 5000,
num_metgrid_levels = 32,
num_metgrid_soil_levels = 4,
dx = 25000,
dy = 25000,
grid_id = 1, 2,
parent_id = 0, 1,
i_parent_start = 1, 10,
j_parent_start = 1, 10,
parent_grid_ratio = 1, 5,
parent_time_step_ratio = 1, 5,
feedback = 1,
smooth_option = 0
/

&physics

mp_physics = 3, 3,
cu_physics = 3, 3,
ra_lw_physics = 3, 3,
ra_sw_physics = 3, 3,
bl_pbl_physics = 1, 1,
sf_sfclay_physics = 1, 1,
sf_surface_physics = 2, 2,
radt = 25, 25,
bldt = 0, 0,
cudt = 0, 0,
icloud = 1,
num_land_cat = 24,
sf_urban_physics = 0, 0,
sst_update = 1,
rdlai2d = .true.
sf_surf_irr_scheme = 0, 1,
irr_daily_amount = 0, 6,
irr_start_hour = 0, 14,
irr_num_hours = 0, 3,
irr_start_julianday = 0, 147,
irr_end_julianday = 0, 211,
irr_ph = 0, 1,
irr_freq = 0, 7,
/

&fdda
/

&dynamics
hybrid_opt = 2,
w_damping = 0,
diff_opt = 1, 1,
km_opt = 4, 4,
diff_6th_opt = 0, 0,
diff_6th_factor = 0.12, 0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 5000., 5000.,
dampcoef = 0.2, 0.2,
khdif = 0, 0,
kvdif = 0, 0,
non_hydrostatic = .true., .true.,
moist_adv_opt = 1, 1,
scalar_adv_opt = 1, 1,
gwd_opt = 1, 1,
/

&bdy_control
spec_bdy_width = 5,
specified = .true.
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

 

[kwerner]

Hi,
Thanks for uploading the full namelist. The irrigation scheme is specific to the NoahMP LSM, and according to your namelist, you are using the Noah LSM. You will need to set sf_surface_physics = 4, 4 if you want the irrigation options to be included in the simulation.

 

[wangfengke]

Thank you very much!