how to use mod_levs.exe?

[pline]

./mod_levs.exe old-intermediate file new-intermediate file
Can you give an example of how it works?
Suppose, if I want to define
&mod_levs
press_pa = 201300 , 200100 , 100000 ,
95000 , 90000 ,
85000 , 80000 ,
75000 , 70000 ,
65000 , 60000 ,
55000 , 50000 ,
45000 , 40000 ,
35000 , 30000 ,
25000 , 20000 ,
15000 , 10000 ,
5000 , 1000
/
What do I have to do?
------------------------------------------------------------------
What I have tried to do is define namelist.wps -> geogrid -> ungrib -> metgrid -> real.exe. It doesn't work.

 

[kwerner]

Can you confirm whether this exact test fails if you don't use mod_levs.exe? I want to make sure that's the reason why your simulation isn't working. Thanks.

 

[pline]

How can I prove it to you?
Now all I can say is What I set in namelist.wps has no effect because when I run wrf.exe the program says

metgrid input_wrf.F first_date_nml = 2019-05-01_00:00:00
d01 2019-06-12_12:00:00 input_wrf.F: SIZE MISMATCH: namelist num_metgrid_levels = 32
d01 2019-06-12_12:00:00 input_wrf.F: SIZE MISMATCH: input file BOTTOM-TOP_GRID_DIMENSION = 34
d01 2019-06-12_12:00:00 ---- ERROR: Mismatch between namelist and input file dimensions
NOTE: 1 namelist vs input data inconsistencies found.
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 1297
NOTE: Please check and reset these options.
-------------------------------------------
Abort(1) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

which I know really needs to be set to 24

 

[kwerner]

Okay, I tested out your case. This is what you need to do:
1) Run geogrid and then ungrib
2) Run the mod_levs.exe program to remove the unwanted levels from the intermediate files. For e.g., from the WPS/ directory,

./util/mod_levs.exe FILE:2019-05-01_00 NEW:2019-05-01_00

and do this for every intermediate file.
3) Run metgrid. Per your settings in &mod_levs in namelist.wps, you should be left with 23 num_metgrid_levels if each of those levels exist in your input data. You can check the number of levels in your met_em* files. For e.g.,

ncdump -h met_em.d01.2019-05-01_00:00:00.nc | grep -i num_metgrid_levels

and you should be given the value of "23." (for the input I used - GFS 0.25 degree data - I was only left with 22 levels)
4) Run real.exe, but you'll need to set "num_metgrid_levels = 23" in namelist.input (or to whatever value you were given in the step above).