I am attempting to run the 3DPwP ocean model using WRF v3.8, with a domain resolution of 7.5km for d01 and 1.5km for d02. The ocean model is only working on d01, and I have noticed that it is recommended to run 3DPwP with a resolution higher than 3km. I am unsure if I am able to run it at a finer resolution, as I still have d03 and d04 domain which has not been nested yet. Thank you for your assistance.
Trouble with running 3DPwP in WRF V3.8
- Thread starter realfake
- Start date
Thanks for your reply! The model does work and it's PWP which simply do not work over d02, with all ocean fields unchanged. As the feedback is turned on, the ocean fields over d01 where are nested with d02 also stay unchanged, and it becomes like this:
This is ocean current of u over d01 and the square area(d02) is all 0, which indicates the PWP model does not work over d02. There is no error message in rsl.err file and my namelist is like this:
&time_control
run_days = 0,
run_hours = 0,
run_minutes = 10,
run_seconds = 0,
start_year = 2019, 2019, 2019,
start_month = 08, 08, 08,
start_day = 06, 06, 04,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2019, 2019, 2019,
end_month = 08, 08, 08,
end_day = 09, 09, 13,
end_hour = 00, 00, 12,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 10800
input_from_file = .true.,.true.,.true.,
history_interval = 10, 60, 60,
frames_per_outfile = 1, 1, 1000,
restart = .false.,
restart_interval = 360,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 100
/
&domains
time_step = 30,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
p_top_requested = 1000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
feedback = 1,
smooth_option = 0,
parent_id = 1, 1,
i_parent_start = 1,134,
j_parent_start = 1, 225,
e_we = 501, 651,
e_sn = 501, 651,
dx = 7500,1500,
dy = 7500,1500,
parent_grid_ratio = 1, 5, 5,
parent_time_step_ratio = 1, 5, 5,
e_vert = 57, 57, 57,
eta_levels = 1.000, 0.997, 0.992, 0.985, 0.978, 0.969, 0.960, 0.950,
0.938, 0.925, 0.910, 0.894, 0.876, 0.857, 0.835, 0.812,
0.787, 0.760, 0.731, 0.700, 0.668, 0.635, 0.600, 0.565,
0.530, 0.494, 0.458, 0.423, 0.388, 0.355, 0.323, 0.293,
0.264, 0.237, 0.212, 0.188, 0.167, 0.147, 0.130, 0.114,
0.099, 0.086, 0.074, 0.064, 0.054, 0.046, 0.039, 0.032,
0.027, 0.022, 0.017, 0.013, 0.010, 0.007, 0.004, 0.002,
0.000,
ocean_levels = 30,
ocean_z = 5, 10, 15, 20, 25, 30,
35, 40, 45, 50, 55, 60,
65, 70, 75, 80, 85, 100,
125, 150, 175, 200, 225, 250,
300, 350, 400, 475, 550, 750
ocean_t = 301.28, 301.12, 300.86, 300.56, 300.22, 299.87,
299.52, 299.20, 298.86, 298.52, 298.19, 298.01,
297.76, 297.46, 297.12, 296.76, 296.39, 295.35,
293.76, 292.20, 290.85, 289.61, 288.52, 287.62,
286.21, 284.99, 283.85, 282.31, 281.00, 278.65
ocean_s = 34.003, 34.028, 34.064, 34.105, 34.148, 34.187,
34.225, 34.274, 34.319, 34.361, 34.411, 34.474,
34.529, 34.573, 34.609, 34.636, 34.659, 34.714,
34.756, 34.760, 34.735, 34.699, 34.660, 34.624,
34.557, 34.497, 34.442, 34.377, 34.336, 34.348
/
&physics
mp_physics = 6, 6, 6, 6, 6, 6,
mp_zero_out = 2,
mp_zero_out_thresh = 1.e-10,
ra_lw_physics = 4, 4, 4, 4, 4, 4,
ra_sw_physics = 4, 4, 4, 4, 4, 4,
radt = 15,15,15,15,15,15,
sf_sfclay_physics = 1, 1, 1, 1, 1, 1,
sf_surface_physics = 2, 2, 2, 2, 2, 2,
bl_pbl_physics = 1, 1, 1, 1, 1, 1,
bldt = 0, 0, 0, 0, 0, 0,
cu_physics = 1, 0, 0, 0, 0, 0,
kfeta_trigger = 2,
cudt = 0, 0, 0, 0, 0, 0,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_land_cat = 21,
num_soil_layers = 4,
maxiens = 1,
maxens = 1,
maxens2 = 1,
maxens3 = 1,
ensdim = 1,
omdt = 6,
sf_ocean_physics = 2,
/
&fdda
/
&dynamics
w_damping = 1,
diff_opt = 1,
diff_6th_opt = 2, 2, 2, 2, 2, 2,
diff_6th_factor = 0.12, 0.12, 0.12, 0.12, 0.12, 0.12,
km_opt = 4,
damp_opt = 1,
zdamp = 5000., 5000., 5000., 5000., 5000., 5000.,
dampcoef = 0.05, 0.01, 0.01, 0.01, 0.01, 0.01,
non_hydrostatic = .true., .true., .true., .true., .true., .true.,
moist_adv_opt = 2, 2, 2, 2, 2, 2,
scalar_adv_opt = 2, 2, 2, 2, 2, 2,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
Your namelist.input looks fine. I only have a concern over your domain setting. Note that PWP works fine over ocean, but it may have trouble if the model domain covers land area. Can you shift your d02 to the east so that it is only over ocean? I would like to see whether PWP can work in this case.
Thanks for your help! The issue has been solved, as I modified the source code in module_sf_ocean_driver.F. The resolution was limited in the source code and now I can run the PWP model with a finer resolution.
Okay, the limitation in source code is like this:
if ( 1.0/rdx .ge. 3000.0 .and. 1.0/rdy .ge. 3000.0 ) then
call PWP()
I simply modify 3000 to 60 and the model works. Since I have not tested the case at a resolution of 60m, it may still have trouble dealing with a finer resolution.