Simulation area expansion error(becaued of e_we, e_sn)
Posted: Wed Mar 27, 2019 6:09 amDear WRF Users
I installed WRF-ARW 3.8.1 version and WPS-3.8.1 version and WPS-4.0.0 version.
i success to run WRF-ARW with below namelist.wps and namelist.input.
Succeed namelist.wps
Succeeded namelist.input
I wanted to expand WRF-ARW simulation area from Korean Peninsula to North-east Asia.
So, I change namelist.wps and namelist.input as below
Failed namelist.wps
Failed namelist.input
Below is message of tail rsl.error.0000
and below is rsl.error.0000 script.
I used NCEP GFS 0P25 data to run WRF model.
I used gcc, gfortran, openmpi, netcdf-4.1.3, Jasper-1.900.1, libpng-1.2.50, and zlib-1.2.7
I think error occur because of e_sn and e_we.
Please let me know what should i change to run WRF toward North-east asia.
and please tell me the reason of those errors.
and I had one more question.
I want to use WRF results as a input data of HYSPLIT Model.
To simulate 200m top on HYSPLIT, I need to make meteorological data using WRF model.
What item should be changed to make a 200m top and land-surface weather file using WRF?
I need your help vigorously.
Best regards.
I installed WRF-ARW 3.8.1 version and WPS-3.8.1 version and WPS-4.0.0 version.
i success to run WRF-ARW with below namelist.wps and namelist.input.
Succeed namelist.wps
- Code: Select all
&share
wrf_core = 'ARW',
max_dom = 1,
start_date = '2019-02-24_00:00:00',
end_date = '2019-02-24_12:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/
&geogrid
parent_id = 1,
parent_grid_ratio = 1,
i_parent_start = 1,
j_parent_start = 1,
e_we = 81,
e_sn = 135,
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
! The default datasets used to produce the HGT_M, GREENFRAC,
! and LU_INDEX/LANDUSEF fields have changed in WPS v3.8. The HGT_M field
! is now interpolated from 30-arc-second USGS GMTED2010, the GREENFRAC
! field is interpolated from MODIS FPAR, and the LU_INDEX/LANDUSEF fields
! are interpolated from 21-class MODIS.
!
! To match the output given by the default namelist.wps in WPS v3.7.1,
! the following setting for geog_data_res may be used:
!
! geog_data_res = 'gtopo_10m+usgs_10m+nesdis_greenfrac+10m','gtopo_2m+usgs_2m+nesdis_greenfrac+2m',
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
geog_data_res = 'default','default',
dx = 5000,
dy = 5000,
map_proj = 'lambert',
ref_lat = 35.88,
ref_lon = 127.5,
truelat1 = 30.0,
truelat2 = 60.0,
stand_lon = -98.0,
geog_data_path = '/home/kict/MODEL/WRF/Build_WRF/WPS_GEOG/'
/
&ungrib
out_format = 'WPS',
prefix = 'FILE',
/
&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/
Succeeded namelist.input
- Code: Select all
&time_control
run_days = 0,
run_hours = 12,
run_minutes = 0,
run_seconds = 0,
start_year = 2019,
start_month = 02,
start_day = 24,
start_hour = 00,
start_minute = 00,
start_second = 00,
end_year = 2019,
end_month = 02,
end_day = 24,
end_hour = 12,
end_minute = 00,
end_second = 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 60,
frames_per_outfile = 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
/
&domains
time_step = 3,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 81,
e_sn = 135,
e_vert = 30,
p_top_requested = 5000,
num_metgrid_levels = 32,
num_metgrid_soil_levels = 4,
dx = 5000,
dy = 5000,
grid_id = 1,
parent_id = 0,
i_parent_start = 1,
j_parent_start = 1,
parent_grid_ratio = 1,
parent_time_step_ratio = 1,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 1, 3, 3,
ra_lw_physics = 4, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 30, 30, 30,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 1, 1, 0,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 3,
num_soil_layers = 4,
num_land_cat = 21,
sf_urban_physics = 0, 0, 0,
/
&fdda
/
&dynamics
w_damping = 0,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
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,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
I wanted to expand WRF-ARW simulation area from Korean Peninsula to North-east Asia.
So, I change namelist.wps and namelist.input as below
Failed namelist.wps
- Code: Select all
&share
wrf_core = 'ARW',
max_dom = 1,
start_date = '2019-02-24_00:00:00','2019-02-24_12:00:00',
end_date = '2019-02-24_12:00:00','2006-08-16_12:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/
&geogrid
parent_id = 1, 1,
parent_grid_ratio = 1, 3,
i_parent_start = 1, 31,
j_parent_start = 1, 17,
e_we = 245, 112,
e_sn = 245, 97,
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
! The default datasets used to produce the HGT_M, GREENFRAC,
! and LU_INDEX/LANDUSEF fields have changed in WPS v3.8. The HGT_M field
! is now interpolated from 30-arc-second USGS GMTED2010, the GREENFRAC
! field is interpolated from MODIS FPAR, and the LU_INDEX/LANDUSEF fields
! are interpolated from 21-class MODIS.
!
! To match the output given by the default namelist.wps in WPS v3.7.1,
! the following setting for geog_data_res may be used:
!
! geog_data_res = 'gtopo_10m+usgs_10m+nesdis_greenfrac+10m','gtopo_2m+usgs_2m+nesdis_greenfrac+2m',
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
geog_data_res = 'default','default',
dx = 10000,
dy = 10000,
map_proj = 'lambert',
ref_lat = 35.88,
ref_lon = 127.5,
truelat1 = 30.0,
truelat2 = 60.0,
stand_lon = -98.0,
geog_data_path = '/home/kict/MODEL/WRF/Build_WRF/WPS_GEOG/'
/
&ungrib
out_format = 'WPS',
prefix = 'FILE',
/
&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/
Failed namelist.input
- Code: Select all
&time_control
run_days = 0,
run_hours = 12,
run_minutes = 0,
run_seconds = 0,
start_year = 2019,
start_month = 02,
start_day = 24,
start_hour = 00,
start_minute = 00,
start_second = 00,
end_year = 2019,
end_month = 02,
end_day = 24,
end_hour = 12,
end_minute = 00,
end_second = 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 60,
frames_per_outfile = 1000,
restart = .false.,
restart_interval = 5000,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 9999
/
&domains
time_step = 60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 245,
e_sn = 245,
e_vert = 30,
p_top_requested = 5000,
num_metgrid_levels = 32,
num_metgrid_soil_levels = 4,
dx = 10000,
dy = 10000,
grid_id = 1,
parent_id = 0,
i_parent_start = 1,
j_parent_start = 1,
parent_grid_ratio = 1,
parent_time_step_ratio = 1,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 0, 3, 3,
ra_lw_physics = 4, 1, 1,
ra_sw_physics = 4, 1, 1,
radt = 30, 30, 30,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 1, 1, 0,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 3,
num_soil_layers = 4,
num_land_cat = 21,
sf_urban_physics = 0, 0, 0,
/
&fdda
/
&dynamics
w_damping = 1,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
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,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
Below is message of tail rsl.error.0000
- Code: Select all
--- NOTE: num_soil_layers has been set to 4
WRF V3.8.1 MODEL
wrf: calling alloc_and_configure_domain
*************************************
Parent domain
ids,ide,jds,jde 1 245 1 245
ims,ime,jms,jme -4 250 -4 250
ips,ipe,jps,jpe 1 245 1 245
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
and below is rsl.error.0000 script.
- Code: Select all
taskid: 0 hostname: kict-Samsung-DeskTop-System
module_io_quilt_old.F 2931 F
Quilting with 1 groups of 0 I/O tasks.
Ntasks in X 1 , ntasks in Y 1
--- WARNING: traj_opt is zero, but num_traj is not zero; setting num_traj to zero.
--- 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: bl_pbl_physics /= 4, implies mfshconv must be 0, resetting
Need MYNN PBL for icloud_bl = 1, resetting to 0
--- NOTE: RRTMG radiation is in use, setting: levsiz=59, alevsiz=12, no_src_types=6
--- NOTE: num_soil_layers has been set to 4
WRF V3.8.1 MODEL
wrf: calling alloc_and_configure_domain
*************************************
Parent domain
ids,ide,jds,jde 1 245 1 245
ims,ime,jms,jme -4 250 -4 250
ips,ipe,jps,jpe 1 245 1 245
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
I used NCEP GFS 0P25 data to run WRF model.
I used gcc, gfortran, openmpi, netcdf-4.1.3, Jasper-1.900.1, libpng-1.2.50, and zlib-1.2.7
I think error occur because of e_sn and e_we.
Please let me know what should i change to run WRF toward North-east asia.
and please tell me the reason of those errors.
and I had one more question.
I want to use WRF results as a input data of HYSPLIT Model.
To simulate 200m top on HYSPLIT, I need to make meteorological data using WRF model.
What item should be changed to make a 200m top and land-surface weather file using WRF?
I need your help vigorously.
Best regards.