Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

How to use vert_refine_method=2 with ndown.exe?

R

robson_passos

New member
I have an output of a simulation with WRF with three nested domains (9 km, 3km and 1 km), from the 1 km output I want to generate lateral and initial conditions for a 200 m domain (0.2 km) using the ndown program. However, I would like this 200m domain to be more refined vertically, so I set vert_refine_method=2 and entered the eta levels I would like to use. The ndown namelist was configured like this in the domains section:

&domains
time_step = 5,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 205, 411,
e_sn = 205, 411,
dzstretch_s = 1.1
p_top_requested = 20000,
vert_refine_method = 0, 2,
rebalance = 1,
e_vert = 45, 53,
eta_levels(1:45) = 1.00000, 0.99629, 0.99257, 0.98879,
0.98486, 0.98071, 0.97622, 0.97130,
0.96585, 0.95977, 0.95299, 0.94540,
0.93692, 0.92744, 0.91686, 0.90507,
0.89195, 0.87737, 0.86120, 0.84331,
0.82356, 0.80181, 0.77793, 0.75181,
0.72335, 0.69246, 0.65911, 0.62329,
0.58506, 0.54455, 0.50195, 0.45755,
0.41175, 0.36503, 0.31802, 0.27144,
0.22617, 0.18317, 0.14344, 0.10788,
0.07710, 0.05132, 0.03028, 0.01343,
0.00000,
eta_levels(46:98) = 1.00000, 0.99814, 0.99629, 0.99442,
0.99257, 0.99068, 0.98879, 0.98682,
0.98486, 0.98278, 0.98071, 0.97846,
0.97622, 0.97376, 0.97130, 0.96857,
0.96585, 0.95977, 0.95299, 0.94540,
0.93692, 0.92744, 0.91686, 0.90507,
0.89195, 0.87737, 0.86120, 0.84331,
0.82356, 0.80181, 0.77793, 0.75181,
0.72335, 0.69246, 0.65911, 0.62329,
0.58506, 0.54455, 0.50195, 0.45755,
0.41175, 0.36503, 0.31802, 0.27144,
0.22617, 0.18317, 0.14344, 0.10788,
0.07710, 0.05132, 0.03028, 0.01343,
0.00000
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 1000,
dy = 1000,
grid_id = 1, 2,
parent_id = 0, 1,
i_parent_start = 1, 61,
j_parent_start = 1, 61,
parent_grid_ratio = 1, 5,
parent_time_step_ratio = 1, 20,
feedback = 0,
smooth_option = 0
tslist_unstagger_winds = .true.,
max_ts_locs = 65,
max_ts_level = 17,
wif_input_opt = 1
num_wif_levels = 30
/
Click to expand...

However, when I analyze the wrfinput_d02 file generated by ndown, the refinement is not applied, I continue with 45 vertical levels, equal the 1km domain. When I use refinement by a factor it works perfectly:

vert_refine_fact = 2
vert_refine_method = 1
rebalance = 1
Click to expand...

My question is: does the ndown program accept vert_refine_method=2?
 
Last edited:
Hello！When I use ndown program, I set vert_refine_fact=4, but I get the following error:
taskid: 0 hostname: login05
module_io_quilt_old.F 2931 T
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 5805
invalid parent id for domain 2
-------------------------------------------
Abort(1) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

I want to ask: why did this error happen? What shoul I do to solve it? Thanks!!!

domain part in my namelist.input is:
&domains
time_step = 2,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 625, 971,
e_sn = 637, 1001,
e_vert = 50, 200,
vert_refine_fact = 1, 3,
vert_refine_method = 0, 1,
rebalance = 1,
p_top_requested = 5000,
num_metgrid_levels = 34,
num_metgrid_soil_levels = 4,
dx = 500,
dy = 500,
grid_id = 1, 2,
parent_id = 0, 1,
i_parent_start = 1, 384,
j_parent_start = 1, 64,
parent_grid_ratio = 1, 5,
parent_time_step_ratio = 1, 4,
feedback = 1,
smooth_option = 0
eta_levels = 1.000, 0.995, 0.990, 0.983, 0.975,
0.965, 0.953, 0.938, 0.920, 0.898,
0.876, 0.854, 0.832, 0.810, 0.788,
0.766, 0.744, 0.722, 0.700, 0.678,
0.656, 0.634, 0.612, 0.590, 0.570,
0.550, 0.530, 0.510, 0.490, 0.470,
0.450, 0.430, 0.410, 0.390, 0.370,
0.350, 0.330, 0.310, 0.290, 0.270,
0.250, 0.225, 0.195, 0.160, 0.120,
0.070, 0.035, 0.015, 0.005, 0.000,
/
 
@pingyl

Please open a new thread for your question. Can you upload your entire namelist.input for me to take a look? Thanks.
 
You must log in or register to reply here.
Top