nested domain HGT=0 - WRF 3.7.1

Discuss the nesting capability within the model itself and any problems you might have run into.

nested domain HGT=0 - WRF 3.7.1

Postby BHajni » Tue Nov 03, 2015 8:47 am

Hi All,


I've encountered a problem with nested domain HGT.
I'm running the 3.7.1 version (problem persists with both distributed and shared memory) with gcc+gfortran. The simulation has 3 domains.
The thing is, that the wrfinput files are fine for all domains, however in the wrfout file of d03, the outer 5 grid rim shows 0 for HGT, and it wrecks up the simulation. With an (almost) same setup (km, diff opt) the simulation is fine using the 3.5 version. (the met_em-s were created with their respective WRF versions). The funny thing is that all other static variables are fine (e.g. landuse).

Does anybody have an idea? There must be something small I'm messing up but I can't find it.

Thanks in advance,
Hajni

For some reason I cannot upload the namelist.input as a file so (the missing options, fdda, fire etc. are not used):
&time_control
run_days = 0,
run_hours = 24,
run_minutes = 0,
run_seconds = 0,
start_year = 2015, 2015, 2015,
start_month = 07, 07, 07,
start_day = 05, 05, 05,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2015, 2015, 2015,
end_month = 07, 07, 07,
end_day = 06, 06, 06,
end_hour = 00, 00, 00,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 1, 1, 1,
adjust_output_times = .true.
frames_per_outfile = 1000, 1000, 1000,
restart = .false.,
restart_interval = 5000,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
io_form_auxinput2 = 2
/

&domains
use_adaptive_time_step = .true.
step_to_output_time = .true.
target_cfl = 0.8, 0.8, 0.8,
target_hcfl = 0.8, 0.8, 0.8,
max_step_increase_pct = 10, 5, 5,
starting_time_step = 60, 18, 6,
adaptation_domain = 3,
time_step = 60
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 100, 109, 109,
e_sn = 65, 76, 94,
e_vert = 44, 44, 44,
eta_levels = 1.0, 0.99366, 0.99083, 0.988, 0.98416, 0.98033, 0.9765, 0.97166,
0.96683, 0.962, 0.956, 0.95, 0.944, 0.9365, 0.929, 0.9215,
0.9125, 0.9035, 0.8945, 0.8746439, 0.8547878, 0.8295757, 0.8043636,
0.7791515, 0.7539394, 0.72872725, 0.7035151, 0.669263, 0.6350235,
0.5527188, 0.4785738, 0.4119146, 0.3521142, 0.2985896, 0.2507996,
0.2082415, 0.1704494, 0.1369916, 0.1074685, 0.08151034, 0.05877567,
0.03894898, 0.02173898, 0.0
p_top_requested = 10000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 10000, 3333.333, 1111.111,
dy = 10000, 3333.333, 1111.111,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 32, 33,
j_parent_start = 1, 28, 22,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 0,
smooth_option = 0
use_tavg_for_tsk = .false.
/

&physics
mp_physics = 3, 3, 3,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 1, 1, 1,
sf_sfclay_physics = 2, 2, 2,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 8, 8, 8,
bldt = 0, 0, 0,
cu_physics = 1, 1, 0,
ishallow = 1,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 0,
surface_input_source = 1,
num_land_cat = 33,
num_soil_layers = 4,
usemonalb = .true.
sf_urban_physics = 0, 0, 1,
do_radar_ref = 1
slope_rad = 1, 1, 1,
topo_shading = 1, 1, 1,
topo_wind = 1, 1, 1,
prec_acc_dt = 15, 15, 15,
num_urban_layers = 5400
/

&dynamics
w_damping = 0,
diff_opt = 1, 1, 2,
km_opt = 4, 4, 3,
diff_6th_opt = 1, 1, 1,
diff_6th_factor = 0.3, 0.2, 0.2,
base_temp = 290.
damp_opt = 3,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2,
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
tke_adv_opt = 1, 1, 1,
time_step_sound = 4
gwd_opt = 1
/

&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/
BHajni
 
Posts: 14
Joined: Sun Mar 10, 2013 11:50 am

Re: nested domain HGT=0 - WRF 3.7.1

Postby sallylove04 » Wed Dec 09, 2015 7:45 pm

I have the same problem. Have you solved it?
sallylove04
 
Posts: 1
Joined: Tue Aug 04, 2015 10:10 am

Re: nested domain HGT=0 - WRF 3.7.1

Postby BHajni » Wed Jan 20, 2016 11:37 am

Sadly no solution yet :( After 3 weeks (including teaching classes and supervising) of digging into the source code and printing out side calculations I couldn't manage it.
I've found that HGT is always recalculated in each time step, and this is done in two steps, one for the relaxation zone and one for the inner domain. And this is where things go boom. For some reason, the interpolation of HGT for the relaxation zone goes wrong, and everything becomes 0.
Not happy about it, but I couldn't spent more time on this, so I've stayed with the 3.5 version so far. Maybe I'll return to it when I'm not swamped with other works. ... or someone else points out the error :)
BHajni
 
Posts: 14
Joined: Sun Mar 10, 2013 11:50 am

Re: nested domain HGT=0 - WRF 3.7.1

Postby BHajni » Thu Mar 10, 2016 12:02 pm

So, someone sent me the solution. :D
The sf_urban_physics setting must be the same in all domains.
BHajni
 
Posts: 14
Joined: Sun Mar 10, 2013 11:50 am


Return to Nesting

Who is online

Users browsing this forum: No registered users and 2 guests