Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

WRF run error with a 3d pwp activated

H

hhuhato

New member
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, I used 2 domains, I was running an ideal experiment that was all ocean, but after one step I found that OM_TMP was having problems. I changed the output step size to output once after a minute of running, and OM_TMP had an exception, as shown in the image.



微信图片_20250215194350.png


If I used one layer of nesting, it would work fine. Is this a matter of interpolation?
The following information appears in the rsl.out file

rrtm: TBOUND exceeds table limit: reset 856.750

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 = 1, 1, 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 = 1.52205585, 4.60833258, 7.75458268, 10.96909038,
14.26501414, 17.66132265, 21.18414231, 24.86861013,
28.76135973, 32.92381034, 37.43648065, 42.4046127 ,
47.96546925, 54.29775966, 61.63375592, 70.27477432,
80.61080764, 93.14516924, 108.5250012 , 127.57830996,
151.3576622 , 181.18953606, 218.7261615 , 265.99288531,
325.41783529, 399.82100015, 492.32614272, 606.1419 ,
744.14216083, 908.17179224

ocean_t = 302.19599953, 302.19680325, 302.19808655, 302.1991679 ,
302.1998305 , 302.20047977, 302.20134148, 302.20265142,
302.20145986, 302.19783994, 302.18973082, 302.17085801,
302.1115981 , 302.07063411, 301.99832816, 301.63736141,
301.02424967, 300.29601385, 299.5333556 , 298.37002274,
296.833204 , 294.71501398, 291.14658614, 287.52511898,
284.34311688, 281.37584612, 279.34000975, 278.28353171,
277.17890437, 276.54511667


ocean_s = 24.61999893, 24.61926883, 24.61644514, 24.61221813, 24.60517728,
24.62803323, 24.66142867, 24.69855108, 25.82220028, 27.03540168,
28.34316302, 29.78513427, 31.40562568, 33.25889308, 34.93485881,
34.99883262, 35.06025718, 35.0962684 , 35.13560762, 35.13610627,
35.05066292, 35.00157597, 34.88618827, 34.7506283 , 34.59570012,
34.33519952, 34.19134865, 34.17492907, 34.32502371, 34.44916072

/
&physics
mp_physics = 2, 2, 6,
ra_lw_physics = 1, 1, 0,
ra_sw_physics = 1, 1, 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 = 1, 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,
sf_ocean_physics = 2,
omdt = 1,
/

&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 = .true., .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,
/
 
Hi,
Can you package all of your rsl* files together into a single *.tar file and attach that?

You mention you are running an ideal case - does this mean you compiled WRF for one of the idealized cases, and if so, which one? Or are you running a real-data case, but just with an ocean-only domain? If it's the latter, can you let me know the type and grid-spacing (resolution) of your meteorological input data? Thanks!
 
Thank you for your reply, I uploaded the rsl file, I used the latter, a real-data case with an ocean-only domain, I used the data is FNL data, its resolution is 1°. The WRF version is 4.1.2.
 

Attachments

  • rsl.tar
    1 MB · Views: 0
Can you change time_step = 60, and keep omdt =1, then rerun this case? Please let me know how it works.

Another issue is, did you compile WRF in vortex-following mode?
 
kwerner said:
Hi,
Can you package all of your rsl* files together into a single *.tar file and attach that?

You mention you are running an ideal case - does this mean you compiled WRF for one of the idealized cases, and if so, which one? Or are you running a real-data case, but just with an ocean-only domain? If it's the latter, can you let me know the type and grid-spacing (resolution) of your meteorological input data? Thanks!
Click to expand...
Thank you for your reply, I uploaded the rsl file, I used the latter, a real-data case with an ocean-only domain, I used the data is FNL data, its resolution is 1°. The WRF version is 4.1.2.
 

Attachments

  • rsl.tar
    1 MB · Views: 0
Ming Chen said:
Can you change time_step = 60, and keep omdt =1, then rerun this case? Please let me know how it works.

Another issue is, did you compile WRF in vortex-following mode?
Click to expand...
Thank you for your reply. When I change time_step to 60 and keep omdt=1, I get the same error. I use vortex-following.
 
I guess I know the issue. 3D PWP cannot be implemented if you run WRF in vortex-following mode. Please recompile WRF and choose the 'basic' option, then rerun this case. Let me know if you still have problems. Thanks.
 
Ming Chen said:
I guess I know the issue. 3D PWP cannot be implemented if you run WRF in vortex-following mode. Please recompile WRF and choose the 'basic' option, then rerun this case. Let me know if you still have problems. Thanks.
Click to expand...
Thanks for your answers, I will recompile and try to run. Sorry, But in this article, I found that the author used a moving three-layer nesting. Did the author make the changes to the wrf himself? The Influence of Oceanic Barrier Layers on Tropical Cyclone Intensity as Determined through Idealized, Coupled Numerical Simulations

Each simulation was performed using a fixed outer domain (d01) on an f plane at 15°N with doubly periodic boundary conditions and a horizontal grid spacing of 18 km, over 320 × 240 grid points in the zonal and meridional directions. Two fully interactive, nested domains (d02 and d03) in the ocean and atmosphere allowed for finer resolutions of 6 and 2 km over square grids of 180 × 180 and 240 × 240 points centered on and moving with the vortex.
 
Last edited:
Ming Chen said:
I guess I know the issue. 3D PWP cannot be implemented if you run WRF in vortex-following mode. Please recompile WRF and choose the 'basic' option, then rerun this case. Let me know if you still have problems. Thanks.
Click to expand...
Thank you for your reply. After I chose ‘basic’ to recompile, the same problem still occurred. As I mentioned in my last reply, other people have used ‘move’ and nested at the same time, which I think may have nothing to do with the compilation method.
 
Sorry that I don't have an immediate solution to the issue you have. Please compile WRF in debug mode, then rerun this case. The log file will tell exactly when and where the model crashes first.
The figure you posted indicates that .OM_TMP = NaN at some points.
 
Ming Chen said:
Sorry that I don't have an immediate solution to the issue you have. Please compile WRF in debug mode, then rerun this case. The log file will tell exactly when and where the model crashes first.
The figure you posted indicates that .OM_TMP = NaN at some points.
Click to expand...
Thanks for your reply, I will recompile and run according to your suggestion.
 
You must log in or register to reply here.
Top