ndown P_TOP error. Please help.

[drsahirazmi]

I am encountering this error when using ndown.exe, for coarser domain input file P_TOP was 5000 ....and after running ndown.exe....P_TOP was 0? and because of this, further ./wrf.exe is giving me error as mismatch between P_TOP in namelist and input files.

Here is my namelist:
&time_control
start_year = 2023, 2023, 2022,
start_month = 06, 06, 02,
start_day = 10, 10, 01,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2023, 2023, 2022,
end_month = 06, 06, 02,
end_day = 12, 12, 03,
end_hour = 18, 18, 00,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 3600,
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 180,
frames_per_outfile = 1000, 1000, 1,
restart = .false,
restart_interval = 120,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
!auxinput4_inname = 'wrflowinp_d<domain>',
!io_form_auxinput4 = 2,
!auxinput4_interval_m = 360,
!auxinput4_end_h = 100000,
auxinput5_inname = 'wrfchemi_d<domain>_<date>',
io_form_auxinput5 = 2,
frames_per_auxinput5 = 1,
auxinput5_interval_m = 3600,
auxinput6_inname = 'wrfbiochemi_d<domain>',
io_form_auxinput6 = 2,
auxinput6_interval_d = 3600,
auxinput7_inname = 'wrffirechemi_d<domain>_<date>',
io_form_auxinput7 = 2,
frames_per_auxinput7 = 1,
auxinput7_interval_m = 3600,
io_form_auxinput2 = 2,
!iofields_filename = "my_file_d01.txt"
!auxhist3_outname = 'wrfout_hourly_d<domain>_<date>',
!auxhist3_interval_m = 60,60,60,
!frames_per_auxhist3 = 1,1,1,
!io_form_auxhist3 = 2,
!auxinput12_inname = 'wrf_chem_input',
!io_form_auxinput12 = 2,
!write_hist_at_0h_rst = .false.,

!auxhist4_outname = 'wrfout_saqib_d<domain>_<date>',
!auxhist4_interval_m = 60,60,60,
!frames_per_auxhist4 = 1,1,1,
!io_form_auxhist4 = 2,
!force_use_old_data = T,

/

&domains
time_step = 18,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
s_we = 1,
s_sn = 1,
e_we = 31, 31, 94,
e_sn = 31, 31, 91,
e_vert = 36, 36, 30,
p_top_requested = 5000,
num_metgrid_levels = 42,
num_metgrid_soil_levels = 4,
dx = 9000, 3000, 3333.33,
dy = 9000, 3000, 3333.33,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 11, 30,
j_parent_start = 1, 11, 30,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 2
eta_levels = 1.0, 0.997, 0.993, 0.988, 0.982, 0.975, 0.966, 0.956, 0.946, 0.933, 0.92, 0.906, 0.891,
0.874, 0.857, 0.838, 0.818, 0.796, 0.773, 0.748, 0.722, 0.693, 0.663, 0.63, 0.595, 0.558,
0.518, 0.475, 0.429, 0.379, 0.326, 0.27, 0.209, 0.144, 0.075, 0.0
/

&physics
mp_physics = 10, 10, 10
progn = 1, 1, 1
ra_lw_physics = 4, 4, 4
ra_sw_physics = 4, 4, 4
radt = 6, 6, 30
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 5, 5, 5,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 5, 5, 3,
cudt = 5, 5, 0,
cu_diag = 0,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 21,
sf_urban_physics = 0, 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .true.,.true.,
sst_update = 0,
! prec_acc_dt = 60,
/

&fdda
grid_fdda = 0, 0, 0,
gfdda_inname = 'wrffdda_d<domain>',
gfdda_end_h = 100000, 100000, 24,
gfdda_interval_m = 360, 360, 360,
fgdt = 0, 0, 0,
if_no_pbl_nudging_uv = 1, 1, 0,
if_no_pbl_nudging_t = 1, 1, 1,
if_no_pbl_nudging_q = 1, 1, 1,
if_zfac_uv = 1, 1, 1,
k_zfac_uv = 10, 10, 10,
if_zfac_t = 0, 0, 0,
k_zfac_t = 10, 10, 10,
if_zfac_q = 0, 0, 0,
k_zfac_q = 10, 10, 10,
guv = 0.0003, 0.0003, 0.0003,
gt = 0.0003, 0.0003, 0.0003,
gq = 0.0003, 0.0003, 0.0003,
if_ramping = 1,
dtramp_min = 60.0,
io_form_gfdda = 2,
/

&dynamics
!etac = 0.15
!hybrid_opt = 2,
rk_ord = 3,
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 = 3,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2,
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
epssm = 0.1, 0.1, 0.1,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 2, 2, 2,
chem_adv_opt = 2, 2, 2,
scalar_adv_opt = 2, 2, 2,
tke_adv_opt = 2, 2, 2,
/

&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .true.,.false.,
nested = .false., .false.,.true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

 

[kwerner]

1) Is the namelist you pasted the one you used to run ndown.exe? Can you attach (not paste) the namelists you used for running the coarse-grid domain, and then ndown.exe, and then the fine-grid?
2) Did you get this error when running ndown.exe or wrf.exe for the fine grid?
3) Can you also please attach the full error log that shows the error message you are receiving?
4) Are you also using WRF-chem? If so, can you try a run without chemistry and see if you still get this error?

 

[kwerner]

@mirandcorzo
Are you using the ndown program? If not, can you post this as a new thread since it's likely a different situation - this helps to keep the threads clean with only a single issue per thread. When you do, please also attach the error file, as well as the namelists you use for the different runs (initial run vs. restart run). Thanks!

 

[Xin]

I'm having the same problem with ndown.exe. Besides P_TOP, the map projections in the metadata and some surface-related metadata also became 0s after running ndown.exe.

Original wrfinput_d02 generated by real.exe:


wrfinput_d02 generated by ndown.exe:


The changed P_TOP and metadata values by ndown.exe cause a series of problems when running the nested domain alone. I wonder whether there are any settings to make sure that ndown.exe presersed the model configurations correctly. If not, will the simulation be OK if one manually correct the wrfbdy and wrfinput generated by ndown.exe before using them?

Thank you very much!
Xin

 

[kwerner]

I'm having the same problem with ndown.exe. Besides P_TOP, the map projections in the metadata and some surface-related metadata also became 0s after running ndown.exe.

Xin, see my first response above for the files I requested. Can you attach those files so I can look into it a bit deeper? Alternatively, if you are able to modify the wrfbdy/input files generated by ndown, I don't think they would cause any issues. You can try it and let us know how it works.

 

[Xin]

Xin, see my first response above for the files I requested. Can you attach those files so I can look into it a bit deeper? Alternatively, if you are able to modify the wrfbdy/input files generated by ndown, I don't think they would cause any issues. You can try it and let us know how it works.

@kwerner,

Thank you for the response! Just an update, after manually correcting P_TOP and metadata in wrfinput (BTW, wrfbdy generated by ndown seems to have correct metadata), WRF was running OK with the inner domain alone.

Thank you for offering the help to look at the namelists! I'm curious whether the corrupted P_TOP and metadata were due to some undesired (or desired but missing) namelist settings or not. Attached are the namelists I used:
namelist.input.CONUS: generate the outer domain initial/boundary conditions, SST updates, and nudging files
namelist.input.ndown: generate the inner domain initial/boundary conditions based on wrfout_d01
namelist.input.d02_alone: used to run d02 alone using initial/boundary conditions generated by ndown.exe

No matter what, it is useful to document the problem and possible solution so that a user can refer to it before the cause of the problem is identified and resolved in future versions of WRF.

Best,
Xin

 

[kwerner]

Xin,
Thanks for updating your status and for sharing the files. I'm glad it's working for you now. When I get time, I'll look into it and if I find anything significant, I'll update this post, as well.