WRF p_lev_diags setting error

[lirui111]

Hello：



I am trying to use the 'pressure level output' option in WRF4.5.2 To save U, V, T at some pressure levels.But according to my namelist output, the file contains variables such as U10, V10, T2, etc., and does not include variables such as U, V, T, dew point temperature, etc. My input file is as follows:





_________namelist.input_______________

&time_control

auxhist23_outname ="wrf_specify_vars_pl_d<domain>_<date>"

auxhist23_interval = 360, 360, 360

frames_per_auxhist23 = 100, 100, 100

io_form_auxhist23 = 2

/

&diag

p_lev_diag = 1

num_press_levels = 16

press_levels = 100000, 95000, 90000, 85000, 80000, 75000, 70000, 50000, 40000, 30000, 20000, 10000, 5000, 3000, 1000, 100

/



Thanks everybody

 

[kwerner]

Hi,
Can you attach your full namelist.input file, your rsl.out.0000 files, and can you issue the following command:

ncdump -h wrf_specify_vars_pl_d01_<date> >& ncdump.txt

(making sure to use the actual file name with the correct date) and attach that ncdump.txt file, as well? Thanks!

 

[lirui111]

Hi,
Can you attach your full namelist.input file, your rsl.out.0000 files, and can you issue the following command:

ncdump -h wrf_specify_vars_pl_d01_<date> >& ncdump.txt

(making sure to use the actual file name with the correct date) and attach that ncdump.txt file, as well? Thanks!

Hi, Mr Kwerner

Thank you for your reply!
Here are my namelist settings, rsl. out. 0000, and ncdump. txt. I only ran data for an hour before stopping WRF. I also modified the press_levels of the namelist.input to not exceed p-top, but it still doesn't work. I hope to receive your reply as soon as possible!
thanks

 

[kwerner]

Thanks for sending those, and I appologize for the delay in response. I think I figured out the issue. You have the following setting:

p_lev_diag = 1

but it should have an "s" on the end, so it should actually be

p_lev_diags = 1

Can you try that and see if that fixes the problem?

 

[lirui111]

Thanks for your reply，Kwerner.
But I tried to change the value of p_lev_diag=1 to p_lev_diags=1, but the result was still the same.

 

[kwerner]

Hi,
I overlooked this last time, but you also need to change the name of the entire namelist record, from &diag, to &diags (again, with an 's'). Try that, and see if that provides the correct results.

 

[lirui111]

Thank you!
This modification successfully solved my problem!

 

[kwerner]

Great! I'm so glad to hear that. Thank you for the update.