[Wrf-users] Problem with three domains in WRF/WRF-Chem
Feng Liu
FLiu at azmag.gov
Fri Jun 10 16:10:22 MDT 2011
The resolution of nested domain (5-km) is higher than you ran 2 domains. So reduce your time step from 270 to 90 and try again. Please keep in mind it will be computationally expensive if works with this solution. Thanks.
Feng
From: wrf-users-bounces at ucar.edu [mailto:wrf-users-bounces at ucar.edu] On Behalf Of Chris Klich
Sent: Wednesday, June 08, 2011 1:28 PM
To: wrf-users at ucar.edu
Subject: [Wrf-users] Problem with three domains in WRF/WRF-Chem
Hi all, I am currently running WRF/WRF-Chem Version 3.2.1. I have been able to successfully run 2 domains, but running 3 has not seemed to work. There is no error in the rsl.error.0000 file, it just stops. I have pasted the rsl.error.0000 file below:
taskid: 0 hostname: sully.met.fsu.edu<http://sully.met.fsu.edu>
Namelist dfi_control not found in namelist.input. Using registry defaults for variables in dfi_control
Namelist tc not found in namelist.input. Using registry defaults for variables in tc
Namelist scm not found in namelist.input. Using registry defaults for variables in scm
Namelist fire not found in namelist.input. Using registry defaults for variables in fire
Ntasks in X 2 , ntasks in Y 4
--- NOTE: sst_update is 0, setting io_form_auxinput4 = 0 and auxinput4_interval = 0 for all domains
--- NOTE: sst_update is 0, setting io_form_auxinput4 = 0 and auxinput4_interval = 0 for all domains
--- NOTE: sst_update is 0, setting io_form_auxinput4 = 0 and auxinput4_interval = 0 for all domains
--- NOTE: grid_fdda is 0 for domain 1, setting gfdda interval and ending time to 0 for that domain.
--- NOTE: both grid_sfdda and pxlsm_soil_nudge are 0 for domain 1, setting sgfdda interval and ending time to 0 for that domain.
--- NOTE: obs_nudge_opt is 0 for domain 1, setting obs nudging interval and ending time to 0 for that domain.
--- NOTE: grid_fdda is 0 for domain 2, setting gfdda interval and ending time to 0 for that domain.
--- NOTE: both grid_sfdda and pxlsm_soil_nudge are 0 for domain 2, setting sgfdda interval and ending time to 0 for that domain.
--- NOTE: obs_nudge_opt is 0 for domain 2, setting obs nudging interval and ending time to 0 for that domain.
--- NOTE: grid_fdda is 0 for domain 3, setting gfdda interval and ending time to 0 for that domain.
--- NOTE: both grid_sfdda and pxlsm_soil_nudge are 0 for domain 3, setting sgfdda interval and ending time to 0 for that domain.
--- NOTE: obs_nudge_opt is 0 for domain 3, setting obs nudging interval and ending time to 0 for that domain.
--- NOTE: num_soil_layers has been set to 4
WRF V3.2.1 MODEL
*************************************
Parent domain
ids,ide,jds,jde 1 103 1 103
ims,ime,jms,jme -4 58 -4 33
ips,ipe,jps,jpe 1 51 1 26
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
alloc_space_field: domain 1 , 198885272 bytes allocated
med_initialdata_input: calling input_input
INPUT LandUse = "USGS"
INITIALIZE THREE Noah LSM RELATED TABLES
LANDUSE TYPE = USGS FOUND 27 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICAION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND 19 CATEGORIES
*********************************************************************
* PROGRAM: WRF/CHEM VERSION 3.1.1 *
* *
* PLEASE REPORT ANY BUGS TO WRF/CHEM HELP at *
* *
* wrfchemhelp.gsd at noaa.gov<mailto:wrfchemhelp.gsd at noaa.gov> *
* *
*********************************************************************
*************************************
Nesting domain
ids,ide,jds,jde 1 211 1 211
ims,ime,jms,jme -4 117 -4 63
ips,ipe,jps,jpe 1 105 1 53
INTERMEDIATE domain
ids,ide,jds,jde 13 88 13 88
ims,ime,jms,jme 8 59 8 42
ips,ipe,jps,jpe 11 49 11 32
*************************************
alloc_space_field: domain 2 , 496635684 bytes allocated
*** Initializing nest domain # 2 from an input file. ***
med_initialdata_input: calling input_input
INPUT LandUse = "USGS"
INITIALIZE THREE Noah LSM RELATED TABLES
LANDUSE TYPE = USGS FOUND 27 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICAION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND 19 CATEGORIES
*********************************************************************
* PROGRAM: WRF/CHEM VERSION 3.1.1 *
* *
* PLEASE REPORT ANY BUGS TO WRF/CHEM HELP at *
* *
* wrfchemhelp.gsd at noaa.gov<mailto:wrfchemhelp.gsd at noaa.gov> *
* *
*********************************************************************
INPUT LandUse = "USGS"
INITIALIZE THREE Noah LSM RELATED TABLES
LANDUSE TYPE = USGS FOUND 27 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICAION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND 19 CATEGORIES
*********************************************************************
* PROGRAM: WRF/CHEM VERSION 3.1.1 *
* *
* PLEASE REPORT ANY BUGS TO WRF/CHEM HELP at *
* *
* wrfchemhelp.gsd at noaa.gov<mailto:wrfchemhelp.gsd at noaa.gov> *
* *
*********************************************************************
INPUT LandUse = "USGS"
INITIALIZE THREE Noah LSM RELATED TABLES
LANDUSE TYPE = USGS FOUND 27 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICAION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND 19 CATEGORIES
*********************************************************************
* PROGRAM: WRF/CHEM VERSION 3.1.1 *
* *
* PLEASE REPORT ANY BUGS TO WRF/CHEM HELP at *
* *
* wrfchemhelp.gsd at noaa.gov<mailto:wrfchemhelp.gsd at noaa.gov> *
* *
*********************************************************************
Timing for Writing wrfout_d01_2008-04-18_00:00:00 for domain 1: 8.69800 elapsed seconds.
mediation_integrate: med_read_wrf_chem_emissions: Read emissions for time 2008-04-18_00:00:00
mediation_integrate: med_read_wrf_chem_emissions: Open file wrfchemi_d01_2008-04-18_00:00:00
Timing for processing lateral boundary for domain 1: 5.72800 elapsed seconds.
WRF NUMBER OF TILES FROM OMP_GET_MAX_THREADS = 1
WRF TILE 1 IS 1 IE 51 JS 1 JE 26
WRF NUMBER OF TILES = 1
*************************************
Nesting domain
ids,ide,jds,jde 1 151 1 151
ims,ime,jms,jme -4 87 -4 48
ips,ipe,jps,jpe 1 75 1 38
INTERMEDIATE domain
ids,ide,jds,jde 73 128 38 93
ims,ime,jms,jme 68 109 33 62
ips,ipe,jps,jpe 71 99 36 52
*************************************
alloc_space_field: domain 3 , 304724904 bytes allocated
*** Initializing nest domain # 3 from an input file. ***
med_initialdata_input: calling input_input
INPUT LandUse = "USGS"
INITIALIZE THREE Noah LSM RELATED TABLES
LANDUSE TYPE = USGS FOUND 27 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICAION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND 19 CATEGORIES
You can see that it stops right after this step, and it looks like most of my namelist is right, but I may be missing something. I have also posted the namelist.input below:
&time_control
run_days = 1,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2008, 2008, 2008,2008,2008,
start_month = 04,04,04,04,04,
start_day = 18,18,18,18,18,
start_hour = 00,00,00,00,00,
start_minute = 00, 00, 00, 00,00,
start_second = 00, 00, 00, 00,00,
end_year = 2008, 2008, 2008, 2008, 2008,
end_month = 04,04,04,04,04,
end_day = 19,19,19,19,19,
end_hour = 00,00,00,00,00,
end_minute = 00, 00, 00, 00,00,
end_second = 00, 00, 00, 00,00,
interval_seconds = 21600,
input_from_file = .true.,.true.,.true.,.true.,.true.,
fine_input_stream = 0, 0, 0, 0, 0,
history_interval = 60,60,60,60,60,
frames_per_outfile = 1,1,1,1,1,
restart = .false.,
restart_interval = 0,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
io_form_auxinput3 = 2,
io_form_auxinput4 = 2,
io_form_auxinput5 = 2,
debug_level = 0,
auxinput5_interval = 2592000,
auxinput11_interval_s = 180, 180, 180, 180, 180,
auxinput11_end_h = 6, 6, 6, 6, 6,
/
&domains
time_step = 270,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
s_we = 1, 1, 1, 1, 1,
e_we = 103,211,151,502,,
s_sn = 1, 1, 1, 1,1,
e_sn = 103,211,151,502,,
s_vert = 1, 1, 1, 1,1,
e_vert = 35, 35, 35, 35,35,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
eta_levels = 1.000, 0.995, 0.98, 0.963, 0.95,0.9372, 0.9247, 0.9122, 0.9001, 0.888,0.8761, 0.8642, 0.8524, 0.84, 0.8322,0.817, 0.79, 0.76, 0.72, 0.68,0.63, 0.59, 0.54, 0.49, 0.435, 0.38, 0.335, 0.2886, 0.2504, 0.2119,0.1722, 0.1309, 0.0875, 0.0418, 0.000,
p_top_requested = 5000,
dx = 45000,15000, 5000, , ,
dy = 45000,15000, 5000, , ,
grid_id = 1, 2, 3, 4,5,
parent_id = 1, 1, 2, 3,4,
i_parent_start = 1,15,75,200,
j_parent_start = 1,15,40,130,
parent_grid_ratio = 1,3,3,3,3,
parent_time_step_ratio = 1,3,3,3,3,
feedback = 1,
smooth_option = 0,
use_adaptive_time_step = .true.,
step_to_output_time = .true.,
target_cfl = 1.2, 1.2, 1.2, 1.2, 1.2,
max_step_increase_pct = 10, 10, 10, 10, 10,
starting_time_step = 270,
max_time_step = -1,
min_time_step = -1,
/
&physics
mp_physics = 2, 2,2,2,2,
ra_lw_physics = 1, 1, 1, 1, 1,
ra_sw_physics = 2, 2, 2, 2, 2,
radt = 10,10,10,10,10,
sf_sfclay_physics = 2, 2, 2, 2, 2,
sf_surface_physics = 2, 2, 2, 2, 2,
bl_pbl_physics = 2, 2, 2, 2, 2,
bldt = 0,0,0,0,0,
cu_physics = 1, 1,1,1,1,
cudt = 5,5,5,5,5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
pxlsm_smois_init = 0,0,0,0,0,
mp_zero_out = 2,
sf_urban_physics = 1,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
sst_update = 0,
cu_rad_feedback = .true.,
cugd_avedx = 1,
/
&fdda
obs_nudge_opt = 0,0,0,0,0
max_obs = 150000,
fdda_start = 0., 0., 0., 0., 0.
fdda_end = 99999., 99999., 99999., 99999., 99999.
obs_nudge_wind = 1,1,1,1,1
obs_coef_wind = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_temp = 1,1,1,1,1
obs_coef_temp = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_mois = 1,1,1,1,1
obs_coef_mois = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_rinxy = 240.,240.,180.,180,180
obs_rinsig = 0.1,
obs_twindo = 0.6666667,0.6666667,0.6666667,0.6666667,0.6666667,
obs_npfi = 10,
obs_ionf = 2, 2, 2, 2, 2,
obs_idynin = 0,
obs_dtramp = 40.,
obs_prt_freq = 10, 10, 10, 10, 10,
obs_prt_max = 10
obs_ipf_errob = .true.
obs_ipf_nudob = .true.
obs_ipf_in4dob = .true.
obs_ipf_init = .true.
/
&dynamics
rk_ord = 3,
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0,0,0,0,0,
diff_6th_factor = 0.12,0.12,0.12,0.12,0.12,
damp_opt = 1,
base_temp = 290.,
zdamp = 5000., 5000., 5000., 5000.,5000.,
dampcoef = 0.01, 0.01, 0.01, 0.01, 0.01,
khdif = 0, 0, 0, 0,0,
kvdif = 0, 0, 0,0,0,
smdiv = 0.1,
emdiv = 0.01,
epssm = 0.1,
time_step_sound = 4,
h_mom_adv_order = 5,
v_mom_adv_order = 3,
h_sca_adv_order = 5,
v_sca_adv_order = 3,
non_hydrostatic = .true., .true., .true.,.true.,.true.,
moist_adv_opt = 2, 2, 2,
scalar_adv_opt = 2, 2, 2,
chem_adv_opt = 2, 2, 2,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,.false.,.false.,
nested = .false., .true., .true.,.true.,.true.,
/
&grib2
/
&chem
kemit = 1,
chem_opt = 301, 301, 301,
bioemdt = 30, 30, 30,
photdt = 30, 30, 30,
chemdt = 270,270,270,
frames_per_emissfile = 1,
io_style_emissions = 2,
emiss_inpt_opt = 1, 1, 1,
emiss_opt = 5, 5, 5,
chem_in_opt = 0, 0, 0,
phot_opt = 0, 0, 0,
gas_drydep_opt = 1, 1, 1,
bio_emiss_opt = 0, 0, 0,
dust_opt = 0,
dmsemis_opt = 0,
seas_opt = 0,
gas_bc_opt = 1, 1, 1,
gas_ic_opt = 1, 1, 1,
aer_bc_opt = 1, 1, 1,
aer_ic_opt = 1, 1, 1,
gaschem_onoff = 1, 1, 1,
aerchem_onoff = 1, 1, 1,
wetscav_onoff = 0, 0, 1,
cldchem_onoff = 0, 0, 1,
vertmix_onoff = 1, 1, 1,
chem_conv_tr = 1, 1, 1,
biomass_burn_opt = 0, 0, 0,
plumerisefire_frq = 270, 270, 270,
aer_ra_feedback = 0, 0, 0,
have_bcs_chem = .false., .false., .false.,
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
Any help would be greatly appreciated. I am using self-made scripts for most of the running process, but it hasn't caused any errors up until this point.
Chris
--
Christopher Klich
Meteorology Graduate Student
Florida State University
(908) 208-9743
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20110610/7559aefe/attachment-0001.html
More information about the Wrf-users
mailing list