no wrfinput for the second and third domain

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

no wrfinput for the second and third domain

Postby ohihi » Tue Dec 17, 2013 6:13 am

Hi

After WPS metgrid run, I got the following 12 met_em* files
met_em.d01.2012-06-18_00:00:00.nc
met_em.d01.2012-06-18_06:00:00.nc
met_em.d01.2012-06-18_12:00:00.nc
met_em.d01.2012-06-18_18:00:00.nc
met_em.d02.2012-06-18_00:00:00.nc
met_em.d02.2012-06-18_06:00:00.nc
met_em.d02.2012-06-18_12:00:00.nc
met_em.d02.2012-06-18_18:00:00.nc
met_em.d03.2012-06-18_00:00:00.nc
met_em.d03.2012-06-18_06:00:00.nc
met_em.d03.2012-06-18_12:00:00.nc
met_em.d03.2012-06-18_18:00:00.nc

However,After excute real.exe
only creat wrfbdy_d01 and wrfinput_d01
wrfinput_d02 and wrfinput_d03 not created for second and third domain
already check namelist.wps and namelist.input for e_we,w_sn......etc.
still can't figure out what's going wrong,

namelist.wps as below

&share
wrf_core = 'ARW',
max_dom = 3,
start_date = '2012-06-18_00:00:00','2012-06-18_00:00:00','2012-06-18_00:00:00',

end_date = '2012-06-18_18:00:00','2012-06-18_18:00:00','2012-06-18_18:00:00',

interval_seconds = 21600
io_form_geogrid = 2,
/
&geogrid
parent_id = 0, 1, 2,
parent_grid_ratio = 1, 3, 3,
i_parent_start = 1, 75, 36,
j_parent_start = 1, 33, 28,
e_we = 180, 100, 103,
e_sn = 105, 100, 121,
geog_data_res = '10m','2m','30s',
dx = 45000,
dy = 45000,
map_proj = 'lambert',
ref_lat = 25.5,
ref_lon = 120.0,
truelat1 = 10.0,
truelat2 = 40.0,
stand_lon = 120.0,
geog_data_path = '/home/????/geog'
/

&ungrib
out_format = 'WPS',
prefix = 'FILE',
/

&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/

&mod_levs
press_pa = 201300 , 200100 , 100000 ,
95000 , 90000 ,
85000 , 80000 ,
75000 , 70000 ,
65000 , 60000 ,
55000 , 50000 ,
45000 , 40000 ,
35000 , 30000 ,
25000 , 20000 ,
15000 , 10000 ,
5000 , 1000


namelist.input as below,

&time_control
run_days = 0,
run_hours = 1,
run_minutes = 0,
run_seconds = 0,
start_year = 2012, 2012, 2012,
start_month = 06, 06, 06,
start_day = 18, 18, 18,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2012, 2012, 2012,
end_month = 06, 06, 06,
end_day = 18, 18, 18,
end_hour = 18, 18, 18,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 360, 360, 360,
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
/

&domains
time_step = 180,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 180, 100, 103,
e_sn = 105, 100, 121,
e_vert = 30, 30, 30,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 45000,
dy = 45000,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 75, 36,
j_parent_start = 1, 33, 28,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 3, 3, 3,
ra_lw_physics = 1, 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 = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0,
/

&fdda
/

&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 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,
/

==================================================
please help me.....
ohihi
 
Posts: 8
Joined: Tue Dec 03, 2013 3:16 am

Re: no wrfinput for the second and third domain

Postby Phillip » Wed Dec 18, 2013 5:20 pm

Please post your namelists using the code tags.

I can't see why it shouldn't work from reading through you namelists. Can you post the output of your real.exe? There might be a hint why those files are not generated. Also, have you checked the contents of your met_em* files? They are there, but they might still have a problem. Use ncview to have a quick look or ncdump to get more detailed information on the contents.

ncdump documentation (the program is included in the netCDF packages)
ncview homepage

But here's another problem with your namelist.input: With this configuration your model will stop running after 1 hour. Set run_hours = 18.
Code: Select all
&time_control
run_days = 0,
run_hours = 1,
run_minutes = 0,
run_seconds = 0,
start_year = 2012, 2012, 2012,
start_month = 06, 06, 06,
start_day = 18, 18, 18,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2012, 2012, 2012,
end_month = 06, 06, 06,
end_day = 18, 18, 18,
end_hour = 18, 18, 18,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
Phillip
 
Posts: 46
Joined: Thu Jul 19, 2012 3:29 am

Re: no wrfinput for the second and third domain

Postby Babak.hashemian » Sat Dec 21, 2013 4:51 am

In addition to seting 'run_hours' variable to "18", Don't forget to link met_em files in the wrf.exe path.
Babak.hashemian
 
Posts: 3
Joined: Thu Dec 19, 2013 3:57 am

Re: no wrfinput for the second and third domain

Postby ohihi » Tue Dec 24, 2013 2:37 am

Thanks for your answer,
already check rsl.out.000*
show as below:
=====================================================================================
Domain 2: Current date being processed: 2012-06-18_00:00:00.0000, which is loop # 1 out of 4
configflags%julyr, %julday, %gmt: 2012 170 0.0000000E+00
metgrid input_wrf.F first_date_input = 2012-06-18_00:00:00
metgrid input_wrf.F first_date_nml = 2012-06-18_00:00:00
dx and dy from file 15000.00 15000.00
dx and dy from namelist 200.0000 200.0000
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 533
DX and DY do not match comparing namelist to the input file
=====================================================================================

so i modified my namelist.input as below;
=====================================
dx = 45000, ===>dx = 45000,15000,5000
dy = 45000, ===>dy = 45000,15000,5000
=====================================

then i finally got the wrfinput_d02 and wrfinput_d03
thanks for ur help^^
ohihi
 
Posts: 8
Joined: Tue Dec 03, 2013 3:16 am

Re: no wrfinput for the second and third domain

Postby spajdy.k » Sat Jun 21, 2014 5:47 am

Hello,
I have a similar problem when running WRF with one nest:
all preprocessing works fine (geo_em.d01.nc, geo_em.d02.nc and all met_em.dxx.YYYY-MM-DD_DD:00:00.nc created), however when I run real.exe, wrfinput_d01 & wrfbdy_d01 are created only for the first domain, and wrfinput_d02 is missing.
The whole program says that it sucesfully completed the task, no fatal call, just a short message at the end of rsl.error.0000:

d01 2011-01-19_18:00:00 alloc_space_field: domain 2,
10705428 bytes allocated
forrtl: severe (64): input conversion error, unit -5, file Internal Formatted Read
Image PC Routine Line Source
real.exe 000000000203FA4A Unknown Unknown Unknown
real.exe 000000000203E546 Unknown Unknown Unknown
real.exe 0000000001FFDE30 Unknown Unknown Unknown
real.exe 0000000001F982DF Unknown Unknown Unknown
real.exe 0000000001F977E7 Unknown Unknown Unknown
real.exe 0000000001FBC0D6 Unknown Unknown Unknown
real.exe 0000000001FBA754 Unknown Unknown Unknown
real.exe 0000000000F97DC5 Unknown Unknown Unknown
real.exe 00000000004A2B91 Unknown Unknown Unknown
real.exe 00000000004A1D9C Unknown Unknown Unknown
libc.so.6 00007F0A6B82AD1D Unknown Unknown Unknown
real.exe 00000000004A1C99 Unknown Unknown Unknown.

Does anyone have any idea what is this problem about?
My namelist.input and namelist.wps are below. My data are from ERA-Interim, I´m running 9 days and data are every 6 hours.

Namelist.wps

&share
wrf_core = 'ARW',
max_dom = 2,
start_date = '2011-01-11_00:00:00','2011-01-11_00:00:00',
end_date = '2011-01-19_18:00:00', '2011-01-11_00:00:00',
interval_seconds = 21600,
io_form_geogrid = 2,
/

&geogrid
parent_id = 1, 1,
parent_grid_ratio = 1, 3,
i_parent_start = 1, 8,
j_parent_start = 1, 9,
e_we = 20, 10,
e_sn = 20, 10,
geog_data_res = '2m','30s',
dx = 30000,
dy = 30000,
map_proj = 'polar',
ref_lat = 78.50,
ref_lon = 17.00,
truelat1 = 30.0,
truelat2 = 60.0,
stand_lon = 20.0,
geog_data_path = '/usit/abel/u1/klaraa/WRF/data/geog',
/

&ungrib
out_format = 'WPS',
prefix = 'SURFILE',
/

&metgrid
fg_name = 'FILE','SURFILE'
io_form_metgrid = 2,
/


Namelist.input


&time_control
run_days = 0,
run_hours = 210,
run_minutes = 0,
run_seconds = 0,
start_year = 2011, 20011,
start_month = 01, 01,
start_day = 11, 11,
start_hour = 00, 00,
start_minute = 00, 00,
start_second = 00, 00,
end_year = 2011, 2011,
end_month = 01, 01,
end_day = 19, 19,
end_hour = 18, 18,
end_minute = 00, 00,
end_second = 00, 00,
interval_seconds = 21600,
input_from_file = .true., .true.,
history_interval = 180, 60,
frames_per_outfile = 1000, 1000,
restart = .false.,
restart_interval = 720,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
/

&domains
time_step = 180,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
s_we = 1, 1,
e_we = 20, 10,
s_sn = 1, 1,
e_sn = 20, 10,
s_vert = 1, 1,
e_vert = 30, 30,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 30000, 10000,
dy = 30000, 10000,
grid_id = 1, 2,
parent_id = 1, 1,
i_parent_start = 1, 8,
j_parent_start = 1, 9,
parent_grid_ratio = 1, 3,
parent_time_step_ratio = 1, 3,
feedback = 1,
smooth_option = 0,
/

&physics
mp_physics = 3, 3,
ra_lw_physics = 1, 1,
ra_sw_physics = 1, 1,
radt = 30, 30,
sf_sfclay_physics = 1, 1,
sf_surface_physics = 2, 2,
bl_pbl_physics = 1, 1,
bldt = 0, 0,
cu_physics = 1, 1,
cudt = 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0,
/

&noah_mp
/

&fdda
/
&logging
/

&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0,
diff_6th_factor = 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000.,
dampcoef = 0.2, 0.2,
khdif = 0, 0,
kvdif = 0, 0,
non_hydrostatic = .true., .true.,
moist_adv_opt = 1, 1,
scalar_adv_opt = 1, 1,
/

&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,
nested = .false., .true.,

/

&grib2
/

&scm
/

&tc
/

&diags
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
spajdy.k
 
Posts: 1
Joined: Fri Jun 20, 2014 8:37 am


Return to Nesting

Who is online

Users browsing this forum: No registered users and 5 guests