One-way nesting with multiple domains

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

One-way nesting with multiple domains

Postby mallyrow » Tue Jul 29, 2014 11:30 am

Hello,
I am new to nested runs using WRF. I have a questions. I am trying to down a one way nest with 12, 3, and 0.5km grid spacing. However, when I got to run wrf after ndown I get "program wrf: error opening wrf_inputd02 for reading ierr = -1021". What is saying? After renaming after running down, it did not create wrfinput_d03 or wrfbdy_d03. Why? Can I not use one way nesting in this fashion?
mallyrow
 
Posts: 8
Joined: Mon Jul 21, 2014 10:47 am

WRF-Chem fails to run: nesting with multiple domains

Postby El Prof » Sun Oct 05, 2014 5:10 pm

Hi
I have a problem in running WRF-Chem for a nested simulation. I am trying to run a coarse grid (20 km) with three nested grids (each of 4 km). WRF-Chem aborts (crush) before starting the run. Real.exe had worked well with all the 4 grids and boundary conditions created with MOZART. when i try to run the model, it crushes. However, when I set the namelist to run with one grid only, the model runs fine, but only for the coarse grid. if I try to run each of the nests as a single domain, the problem still persist.
I get the message as:

*** subr move_sections - method = 20
*** subr move_sections - idiag = 0
d01 2008-08-23_00:00:00 initializing dry dep (wesely)
ftuv_init: failed to allocate luse2usgs array
Rank 2 [Sun Oct 5 21:46:54 2014] [c1-0c1s15n2] application called MPI_Abort(MPI_COMM_WORL
D, 1) - process 2
forrtl: error (76): Abort trap signal
Image PC Routine Line Source
wrf.exe 000000000345FD99 Unknown Unknown Unknown
wrf.exe 000000000345E710 Unknown Unknown Unknown
wrf.exe 0000000003414082 Unknown Unknown Unknown
wrf.exe 00000000033A8183 Unknown Unknown Unknown
wrf.exe 00000000033AEEF2 Unknown Unknown Unknown
wrf.exe 000000000339A1D0 Unknown Unknown Unknown
wrf.exe 000000000339A17B Unknown Unknown Unknown
wrf.exe 000000000347BAC1 Unknown Unknown Unknown
wrf.exe 000000000327FC92 Unknown Unknown Unknown
wrf.exe 000000000325F293 Unknown Unknown Unknown
wrf.exe 00000000032392A5 Unknown Unknown Unknown
wrf.exe 0000000000A9CF9F Unknown Unknown Unknown
wrf.exe 000000000194E84C Unknown Unknown Unknown
wrf.exe 0000000001D0A642 Unknown Unknown Unknown
wrf.exe 000000000146C79E Unknown Unknown Unknown
wrf.exe 00000000012EF56D Unknown Unknown Unknown
wrf.exe 00000000012EAABF Unknown Unknown Unknown
wrf.exe 00000000012B939A Unknown Unknown Unknown
wrf.exe 00000000004A2F35 Unknown Unknown Unknown
wrf.exe 0000000000400FE1 Unknown Unknown Unknown
wrf.exe 0000000000400F99 Unknown Unknown Unknown
wrf.exe 0000000000400F26 Unknown Unknown Unknown
wrf.exe 0000000003475824 Unknown Unknown Unknown

If this is something [not sure?] to do with landuse data, I dont know how to solve it. I tied to check my input files files for the nests and the WPS met_em* files, but couldnt figure out how to resolve this. I apppreciated if some can help/advice me please. my namelist looks like this:

$RP_INPUT
&time_control !test
run_days = 0,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2008, 2008, 2008, 2008,
start_month = 08, 08, 08, 08,
start_day = 23, 23, 23, 23,
start_hour = 00, 00, 00, 00,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2008, 2008, 2008, 2008,
end_month = 08, 08, 08, 08,
end_day = 26, 26, 26, 26,
end_hour = 06, 06, 06, 06,
end_minute = 00, 00, 00, 00,
end_second = 00, 00, 00, 00,
interval_seconds = 21600,
input_from_file = .true., .true., .true., .true.,
history_interval = 60, 60, 60, 60,
frames_per_outfile = 1, 1, 1, 1,
restart = .false.,
restart_interval = 1440,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
auxinput4_interval = 360, 360, 360, 360,
!io_form_auxhist2 = 2,
io_form_auxhist3 = 2,
!auxhist23_interval = 60, 60, 60, 60,
auxinput1_inname = "met_em.d<domain>.<date>",
auxinput4_inname = "wrflowinp_d<domain>",
auxinput5_inname = 'wrfchemi_d<domain>_<date>',
auxinput6_inname = 'wrfbiochemi_d01',
auxinput7_inname = 'wrffirechemi_d<domain>_<date>',
auxinput8_inname = 'wrfchemi_gocart_bg_d<domain>',
auxhist23_outname = "wrfpress_d<domain>_<date>",
!auxinput12_inname = "wrf_chem_input",
auxinput5_interval_m = 86400, 86400, 86400, 86400,
auxinput7_interval_m = 86400, 86400, 86400, 86400,
auxinput8_interval_m = 86400, 86400, 86400, 86400,
!io_form_auxinput2 = 2,
io_form_auxinput4 = 2,
io_form_auxinput5 = 2,
io_form_auxinput6 = 2,
io_form_auxinput7 = 2,
io_form_auxinput8 = 2,
!io_form_auxhist23 = 2,
!io_form_auxinput12 = 2,
frames_per_auxinput5 = 1,
frames_per_auxinput6 = 1,
frames_per_auxinput7 = 1,
frames_per_auxinput8 = 1,
!frames_per_auxhist23 = 100, 100, 100, 100,
!iofields_filename = 'add_bio',
!iofields_filename = 'add_fire',
!iofields_filename = 'add_rad',
iofields_filename = 'add_all',
ignore_iofields_warning = .true.,
debug_level = 100,
/
&domains
time_step = 120,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 4,
s_we = 1, 1, 1, 1,
e_we = 321, 241, 281, 291,
s_sn = 1, 1, 1, 1,
e_sn = 251, 241, 251, 281,
s_vert = 1, 1, 1, 1,
e_vert = 41, 41, 41, 41,
p_top_requested = 5000,
eta_levels = 1.000, 0.996, 0.991, 0.986, 0.981, 0.971, 0.961, 0.9
51, 0.943,
0.933, 0.919, 0.895, 0.871, 0.847, 0.823, 0.799,
0.775, 0.751, 0.716, 0.678, 0.636,
0.590, 0.543, 0.493, 0.442, 0.392,
0.341, 0.294, 0.250, 0.209, 0.182,
0.156, 0.133, 0.111, 0.090, 0.069,
0.048, 0.030, 0.015, 0.005, 0.000
num_metgrid_levels = 38,
dx = 20000, 4000, 4000, 4000,
dy = 20000, 4000, 4000, 4000,
grid_id = 1, 2, 3, 4,
parent_id = 1, 1, 1, 1,
i_parent_start = 1, 45, 180, 148,
j_parent_start = 1, 177, 170, 78,
parent_grid_ratio = 1, 5, 5, 5,
parent_time_step_ratio = 1, 5, 5, 5,
feedback = 0,
smooth_option = 0,
zap_close_levels = 50,
interp_type = 1,
t_extrap_type = 2,
extrap_type = 2,
force_sfc_in_vinterp = 0,
use_levels_below_ground = .true.,
use_surface = .true.,
sfcp_to_sfcp = .true.,
lagrange_order = 1,
!! Options for adaptive timestep
!use_adaptive_time_step = .false.,
!step_to_output_time = .true.,
!target_cfl = 1.2, 1.2, 1.2,
!max_step_increase_pct = 5, 51, 51,
!starting_time_step = -1, -1, -1,
!max_time_step = -1, -1, -1,
!min_time_step = -1, -1, -1,
!adaptation_domain = 1,
/
&physics
progn = 1, 1, 1, 1,
mp_physics = 10, 10, 10, 10,
ra_lw_physics = 4, 4, 4, 4,
ra_sw_physics = 4, 4, 4, 4,
radt = 20, 4, 4, 4,
sf_sfclay_physics = 1, 1, 1, 1,
sf_surface_physics = 2, 2, 2, 2,
bl_pbl_physics = 1, 1, 1, 1,
bldt = 0, 0, 0, 0,
cu_physics = 5, 0, 0, 0,
cudt = 0, 0, 0, 0,
sf_urban_physics = 1, 1, 1, 1,
cugd_avedx = 1,
isfflx = 1,
ifsnow = 1,
icloud = 1,
do_radar_ref = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 24,
num_soil_cat = 16,
mp_zero_out = 2,
no_mp_heating = 0,
mp_zero_out_thresh = 1.e-12,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .true., .false., .false., .false.,
cu_diag = 1., 0., 0., 0.,
! long simulation options
sst_update = 1,
usemonalb = .true.,
tmn_update = 0,
sst_skin = 0,
!bucket_mm = 100
/
&fdda ! Options for nudging
!grid_fdda = 1, 1,
!gfdda_inname = "wrffdda_d<domain>",
!!gfdda_end_h = 48, 48, ! When to stop nudging - need to set for entir
e period
!gfdda_interval_m = 360, 360,
!fgdt = 0, 0,
!if_no_pbl_nudging_uv = 0, 0,
!if_no_pbl_nudging_t = 1, 1,
!if_no_pbl_nudging_q = 1, 1,
!if_zfac_uv = 1, 1,
! k_zfac_uv = 19, 19,
!if_zfac_t = 0, 0,
! k_zfac_t = 10, 10,
!if_zfac_q = 0, 0,
! k_zfac_q = 10, 10,
!guv = 0.0003, 0.0003,
!gt = 0.0003, 0.0003,
!gq = 0.0003, 0.0003,
!if_ramping = 0,
!dtramp_min = 60,
!io_form_gfdda = 2,
/
&dynamics
rk_ord = 3,
w_damping = 1,
damp_opt = 0,
km_opt = 4,
zdamp = 5000., 5000., 5000., 5000.,
dampcoef = 0.01, 0.01, 0.01, 0.01,
diff_opt = 1,
mix_full_fields = .false., .false., .false. .false.,
diff_6th_opt = 0, 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12, 0.12,
base_pres = 100000.,
base_temp = 290.,
base_lapse = 50.,
khdif = 0, 0, 0, 0,
kvdif = 0, 0, 0, 0,
smdiv = 0.1, 0.1, 0.1, 0.1,
emdiv = 0.01, 0.01, 0.01, 0.01,
epssm = 0.1, 0.1, 0.1, 0.1,
non_hydrostatic = .true., .true., .true., .true.,
moist_adv_opt = 2, 2, 2, 2,
tke_adv_opt = 2, 2, 2, 2,
scalar_adv_opt = 2, 2, 2, 2,
chem_adv_opt = 2, 2, 2, 2,
tracer_adv_opt = 0, 0, 0, 0,
tracer_opt = 0, 0, 0, 0,
time_step_sound = 4, 4, 4, 4,
h_mom_adv_order = 5, 5, 5, 5,
v_mom_adv_order = 3, 3, 3, 3,
h_sca_adv_order = 5, 5, 5, 5,
v_sca_adv_order = 3, 3, 3, 3,
!tke_upper_bound = 1000., 1000., 1000., 1000.,
!tke_drag_coefficient = 0., 0., 0., 0.,
!tke_heat_flux = 0., 0., 0., 0.,
!do_avgflx_em = 1,
!do_avgflx_cugd = 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false., .false., .false.,
nested = .false., .true., .true., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
&dfi_control
!dfi_opt = 0,
!dfi_nfilter = 0,
!dfi_write_filtered_input = .false.,
!dfi_write_dfi_history = .false.,
!dfi_cutoff_seconds = 0,
!dfi_time_dim = 0,
/
&diags
!plev_diags = 1,
!num_press_levels = 10,
!press_levels = 1, 1, 1, 1,
/
&logging
/
&noah_mp
/
&sc
/
&tc
/
&fire
/
&chem
chem_opt = 10, 10, 10, 10,
bioemdt = 30, 30, 30, 30,
photdt = 30, 30, 30, 30,
chemdt = 3, 3, 3, 3,
kemit = 1,
!frames_per_emissfile = 0,
io_style_emissions = 1,
emiss_inpt_opt = 101, 101, 101, 101,
emiss_opt_vol = 0, 0, 0, 0,
emiss_opt = 5, 5, 5, 5,
chem_in_opt = 0, 0, 0, 0,
phot_opt = 2, 2, 2, 2,
gas_drydep_opt = 1, 1, 1, 1,
aer_drydep_opt = 1, 1, 1, 1,
bio_emiss_opt = 3, 3, 3, 3,
gas_bc_opt = 0, 0, 0, 0,
gas_ic_opt = 0, 0, 0, 0,
aer_bc_opt = 0, 0, 0, 0,
aer_ic_opt = 0, 0, 0, 0,
aer_op_opt = 2, 2, 2, 2,
gaschem_onoff = 1, 1, 1, 1,
aerchem_onoff = 1, 1, 1, 1,
wetscav_onoff = 1, 1, 1, 1,
cldchem_onoff = 1, 1, 1, 1,
vertmix_onoff = 1, 1, 1, 1,
chem_conv_tr = 1, 0, 0, 0,
conv_tr_wetscav = 1, 0, 0, 0,
conv_tr_aqchem = 1, 0, 0, 0,
!opt_pars_out = 1, 1, 1, 1,
seas_opt = 2,
dust_opt = 2,
dmsemis_opt = 0,
ne_area = 140,
biomass_burn_opt = 1, 1, 1, 1,
plumerisefire_frq = 30, 30, 30, 30,
have_bcs_chem = .true., .true., .true., .true.,
!scale_fire_emiss = .true., .true., .true., .true.,
have_bcs_tracer = .false., .false., .false., .false.,
aer_ra_feedback = 1, 1, 1, 1,
/
$END

Suprisingly, I tried to run the meteorology only (WRF) with the same settings I have for WRF-Chem, but turning off the chemistry but I ddnt get the problem with WRF. The simulations works ok for all the nests. I dnt understand how sort the problem in WRF-Chem.
Please help... :? :?
El Prof
 
Posts: 4
Joined: Thu Sep 06, 2012 11:32 pm


Return to Nesting

Who is online

Users browsing this forum: No registered users and 1 guest