Error in p_lev_diags and z_lev_diags

[Israt246]

Hello,
I am trying to run WRF v4.2.2 forced by ERA5. WPS was successful and I am trying to run real.exe now. I want to output pressure and z level variables for my inner domain only. And I set the necessary auxhist options in namelist.input. But when I run real.exe, I get the following error:


--- ERROR: p_lev_diags requires auxhist23 file information
--- ERROR: provide: auxhist23_interval (max_dom) and io_form_auxhist23
--- Add supporting IO for stream 23 for pressure-level diags
--- ERROR: z_lev_diags requires auxhist22 file information
--- ERROR: provide: auxhist22_interval (max_dom) and io_form_auxhist22
--- Add supporting IO for stream 22 for height-level diags
--- WARNING: If use_adaptive_time_step, must use cudt=0 for the following CU schemes:
--- BMJ, all SAS, Tiedtke
--- CUDT=0 has been done for you.
--- NOTE: RRTMG radiation is used, namelist value for o3input (ozone input) is used
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 2367
NOTE: 2 namelist settings are wrong. Please check and reset these options
-------------------------------------------
MPICH ERROR [Rank 0] [job id ] [Mon Jan 27 22:56:37 2025] [dec1909] - Abort(1) (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

But I have already used auxhist23_interval, io_form_auxhist23, auxhist22_interval and io_form_auxhist22 for p and z level variables for the inner domain.
I have attached my nameless.input and rsl.error.0000 files.
Any help is appreciated.
-Israt

 

[Israt246]

I made the following changes to my nameless.input:

auxhist23_interval=360,60
auxhist22_interval= 360, 60
auxhist3_interval=360,60
auxhist7_interval=0,60

WRF is running now. I am perplexed why setting auxhist7_interval=0,60 works but it did not work with old namelist.input settings as:

auxhist23_interval=0,60
auxhist22_interval= 0,60
auxhist3_interval=0,60
auxhist7_interval=0,60

As I am doing long term WRF simulation, there will be lot of files and I do not want p level, z level and diagnostic variables for the outer domain. Is there any way to output p lev, z lev and diagnostic variables for domain 2 only?

 

[kwerner]

Hi,
I'd like to apologize for the very long delay in response while we've had to tend to other obligations. Thank you for your patience.

I'm not sure why that would be happening. What happens if you set the value for domain 01 to something really large - longer than your simulation? Would that give you the results you want?

 

[Israt246]

Hi,
I'd like to apologize for the very long delay in response while we've had to tend to other obligations. Thank you for your patience.

I'm not sure why that would be happening. What happens if you set the value for domain 01 to something really large - longer than your simulation? Would that give you the results you want?

Thank you for your reply. I did not try setting the value of auxhist23_interval, auxhist22_interval, and auxhist3_interval to something large. Right now the options are set as:

auxhist23_interval=360,60
auxhist22_interval= 360, 60
auxhist3_interval=360,60
auxhist7_interval=0,60

So I am getting p-level, z-level and diagnostic variables for domain 1 although my intention was not to output them in the d01 domain. However, the auxhist7_interval worked fine when I set the output interval for domain d01 as 0 ( no output). Stream 7 includes some WRF variables that I want hourly. I have attached my current namelist.input and text files that contain the runtime I/O options.

 

[kwerner]

Can you please try to set the auxhistX_interval values to a length of time longer than your simulation? I'm not sure why the different aux options perform in different ways, but if you can simply change the setting, and that resolves the issue for your case, that would be a much easier solution than if we have to comb through the code to determine why they differ. Please let me know the outcome. Thanks!