Save nested files with different number of levels

[StephM]

Hi

I'm running a 4 nested model with 64 levels with WRF version 4.4

I want to save all the levels in the wrfout files for d01,d02, and d03, but for d04 I want to keep the values from the surface to a given level (my estimate is that I would save 40 or 35 levels from those 64)

I've been trying to find a way to do this, but I've only found information on nesting eta_levals which is different from what I want.

So I tried changing these parameters from the namelist.input, and it does not work either.

e_vert = 64,64,64,40,

These are my eta levels
eta_levels = 1.0, 0.981723780, 0.963447560, 0.945171340, 0.926895119,
0.908618899, 0.902128764, 0.895638629, 0.889148494, 0.882658359,
0.876168224, 0.869678089, 0.863187954, 0.856697819, 0.850207684,
0.843717549, 0.837227414, 0.830737279, 0.824247144, 0.817757009,
0.811266874, 0.804776739, 0.798286604, 0.791796469, 0.785306334,
0.778816199, 0.772326064, 0.765835929, 0.759345794, 0.752855659,
0.746365524, 0.739875389, 0.733385254, 0.726895119, 0.707684320,
0.695223261, 0.682762201, 0.656386293, 0.628868120, 0.600103842,
0.570301142, 0.539460021, 0.507892004, 0.475493250, 0.442679128,
0.409449637, 0.376220145, 0.342990654, 0.310072690, 0.277777778,
0.246313603, 0.215887850, 0.186708204, 0.159190031, 0.133229491,
0.109345794, 0.087538941, 0.067912773, 0.050674974, 0.035721703,
0.023052960, 0.012772586, 0.004569055, 0.000000000,

I get this error

--- ERROR: e_vert differs on grid ids 4 and 3. Set vert_refine_method or make e_vert consistent.


I there a way to modify the namelist.input file to change the number of levels saved on each wrfout domain

Regards
Stefi

 

[kwerner]

Hi Stefi,
Apologies for the delay in response. We have been working to prepare for our upcoming tutorial and haven't had much time to devote to the forum. When you refer to "saving" the levels, do you mean that you only want those levels printed to the history files (wrfout*), but still want to run the model with the full 64 levels? If that's the case, I don't believe there is an easy runtime option for that. You could have everything printed to the history files and then use some sort of post-processor (e.g., NCL, or NetCDF Tools) to remove the layers you don't want.

If I'm not understanding correctly, do you mind trying to clarify it? Thank you!

 

[StephM]

Hi Kwerner

Indeed I need what you are saying. I have managed to do it with nco

Thank you very much. Hopefully the upcoming tutorial will be a success
Stefi

 

[Chetan]

Hi Kwerner

Indeed I need what you are saying. I have managed to do it with nco

Thank you very much. Hopefully the upcoming tutorial will be a success
Stefi

Hi @StephM , could you please share how you did with nco?

I am running the wrf for 64 vertical layers, and want to keep limited vertical layers only.

Thank you.

 

[StephM]

Hello Chetan

This is how if solved it

- load/install nco
- run: ncks -d variable_name, levelid,levelid input output

Just make sure, you cut all the vertical dimensions, some variables have a slightly different name. Here is an example, I also cut some in lat and lon to remove some instabilities around the border.

ncks -d bottom_top,1,35 -d bottom_top_stag,1,35 -d south_north,72,372 -d south_north_stag,72,372 -d west_east,72,372 -d west_east_stag,72,372 wrfout_d04_2018-02-17_14\:00\: wrfout_check

Hope this helps

 

[Chetan]

@StephM Thank you very much for your response.