WRF run error with a 3d pwp activated

[hhuhato]

Dear all,

A simple 3D PWP ocean model was added since WRFv4.1.
When I tried to run with the PWP model and set the sf_ocean_physics = 2, then when I checked the output history files, the record would never update with showing the following information at the end of the file：
--------------------------------------------------------------------------------------
Timing for processing lateral boundary for domain 1: 0.52612 elapsed seconds
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE 1 IS 1 IE 31 JS 1 JE 9
WRF NUMBER OF TILES = 1
dx 12000.00
3DPWP domain 1 run ocean
--------------------------------------------------------------------------------------
PS: WRF run was success when sf_ocean_physics = 0 and 1.
Thank you in advance!!!
the following is my namelist.input

&time_control
run_days = 5,
run_hours = 12,
run_minutes = 0,
run_seconds = 0,
start_year = 2017, 2017, 2017,
start_month = 08, 08, 08,
start_day = 21, 21, 21,
start_hour = 12, 12, 12,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2017, 2017, 2017
end_month = 08, 08, 08,
end_day = 27, 27, 27,
end_hour = 00, 00, 00,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600,
input_from_file = .true., .false., .false.,
history_interval = 60, 60, 60,
frames_per_outfile = 1000, 1000, 1000,
restart = .false.,
restart_interval = 180,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,

ignore_iofields_warning = .false.,


/
&domains
time_step = 30,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 400, 151, 181,
e_sn = 400, 151, 181,
e_vert = 35, 35, 35,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 12000, 4000, 4000,
dy = 12000, 4000, 4000,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 200, 140,
j_parent_start = 1, 175, 90,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0,
track_level = 50000,
time_to_move = 8400, 0, 0,
ocean_levels = 30,

ocean_z = 5., 15., 25., 35., 45., 55.,
65., 75., 85., 95., 105., 115.,
125., 135., 145., 155., 165., 175.,
185., 195., 210., 230., 250., 270.,
290., 310., 330., 350., 370., 390.
ocean_t= 302.3493, 302.3493, 302.3493, 302.1055, 301.9763, 301.6818,
301.2220, 300.7531, 300.1200, 299.4778, 298.7443, 297.9194,
297.0883, 296.1443, 295.1941, 294.1979, 293.1558, 292.1136,
291.0714, 290.0293, 288.7377, 287.1967, 285.6557, 284.8503,
284.0450, 283.4316, 283.0102, 282.5888, 282.1674, 281.7461
ocean_s = 34.0127, 34.0127, 34.0127, 34.3217, 34.2624, 34.2632,
34.3240, 34.3824, 34.3980, 34.4113, 34.4220, 34.4303,
34.6173, 34.6409, 34.6535, 34.6550, 34.6565, 34.6527,
34.6490, 34.6446, 34.6396, 34.6347, 34.6297, 34.6247,
34.6490, 34.6446, 34.6396, 34.6347, 34.6297, 34.6247



/
&physics
sf_ocean_physics = 2,
omdt = 1,
mp_physics = 2, 2, 6,
ra_lw_physics = 0, 0, 0,
ra_sw_physics = 0, 0, 0,
radt = 10, 10, 10,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 6, 0, 0,
cudt = 5, 5, 5,
isftcflx = 2,
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 = 1,
km_opt = 4,
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.01, 0.01, 0.01,
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .false., .false., .true.,
moist_adv_opt = 1, 1,
scalar_adv_opt = 1, 1,
use_baseparam_fr_nml = .true.,

/

&bdy_control
spec_bdy_width = 10,
spec_zone = 1,
relax_zone = 9,
specified = .true., .false., .false.,
nested = .false., .true., .true.,
periodic_x = .false., .false.,
symmetric_xs = .false., .false.,
symmetric_xe = .false., .false.,
open_xs = .false., .false.,
open_xe = .false., .false.,
periodic_y = .false., .false.,
symmetric_ys = .false., .false.,
symmetric_ye = .false., .false.,
open_ys = .false., .false.,
open_ye = .false., .false.,

/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,

 

[kwerner]

Apologies for the long delay in response while our time has been out much over the past few weeks, due to conferences and holidays. Thank you for your patience.

Can you please attach the log files? If you have multiple rsl.error* and rsl.out* files, please package those together into a single *.tar file and attach that. Thanks!

 

[hhuhato]

Thank you for your reply. I have found the error.
First of all, if you want to make an ideal file, XLAND needs to be modified in wrfinput to ensure that it is all ocean.
Secondly, because I need to close the long and short wave radiation, it is necessary to modify the absorb function in the F file. In the absorb function, the long and short wave radiation is used as the denominator, so it cannot be 0. An if statement can be added to judge.

Thanks again for your reply！