Creating met_em files with wrf_interp

[mqin9]

Hi all,

I'm trying to interpolate wrfout data to pressure level, with the output being used as input to real.exe or obsgrid.exe.

I know that the tool p_interp can create met_em file with wrfout data, which is exactly what I want. However, there was a error message suggesting segmentation fault with the tool, and the results are weird when it was ran in parallel.

Then I tried wrf_interp. It can be completed successfully, and the results seem reasonable. However, there was a issue with the vertical layers. The original met_em files have 40 layers, surface plus other 39 layers, ranging from 25mb to 1000mb with increment of 25mb, while there is only 39 layers for the output of wrf_interp. Do I need to add the values at surface layer to the 4D fields like TT, PRES, etc. to make them identical to those in the original met_em files in terms of vertical layers? Is there any scripts that can do such a thing?

Thanks in advance.