I am running WRFv3.9.1 prepared using 'gfortran' on LINUX Fedora v28. I have successfully processed WPS output data into WRF input files and run WRF for a nested domain of 36-, 12-, and 4-km. However, when I try to implement 'analysis nudging', 'real.exe' fails with the following:
Domain 1: Current date being processed: 1990-08-01_18:00:00.0000, which is loop # 2 out of 29
configflags%julyr, %julday, %gmt: 1990 213 18.0000000
metgrid input_wrf.F first_date_input = 1990-08-01_18:00:00
metgrid input_wrf.F first_date_nml = 1990-08-01_12:00:00
d01 1990-08-01_18:00:00 Timing for input 0 s.
d01 1990-08-01_18:00:00 flag_soil_layers read from met_em file is 1
None of the following are allowed to be TRUE :
( internal_time_loop .EQ. 1 ) F
( config_flags%grid_fdda .NE. 0 ) T
( config_flags%sst_update .EQ. 1 ) F
( config_flags%all_ic_times ) F
( config_flags%smooth_cg_topo ) F
( config_flags%map_proj .EQ. PROJ_CASSINI ) F
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 338
Problems, we cannot have excluded middle data from WPS
The central portion of my 'namelist.input' file follows. Basically, everything works fine until I alter 'grid_fdda' from '0' to '1'. Any assistance will be gratefully appreciated.
BJ
&domains
time_step = 120,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 40, 91, 181,
e_sn = 40, 91, 217,
e_vert = 31, 31, 31,
p_top_requested = 10000,
num_metgrid_levels = 30,
num_metgrid_soil_levels = 4,
dx = 36000, 12000, 4000,
dy = 36000, 12000, 4000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 6, 9,
j_parent_start = 1, 9, 6,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
eta_levels = 1.0000, 0.9958, 0.9907, 0.9846, 0.9774,
0.9688, 0.9585, 0.9463, 0.9319, 0.9148,
0.8946, 0.8709, 0.8431, 0.8107, 0.7733,
0.7301, 0.6809, 0.6254, 0.5635, 0.4954,
0.4218, 0.3536, 0.2930, 0.2394, 0.1917,
0.1495, 0.1120, 0.0788, 0.0493, 0.0232,
0.00000,
/
&physics
mp_physics = 6, 6, 6,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 5, 5, 5,
sf_urban_physics = 0, 0, 0,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 1, 1, 1,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 5, 5, 5,
cudt = 5, 5, 5,
cu_rad_feedback = .true., .true., .true.,
sst_update = 0,
usemonalb = .true.,
rdmaxalb = .true.,
surface_input_source = 3,
num_soil_layers = 4,
num_land_cat = 21,
/
&fdda
grid_fdda = 1, 1, 0,
gfdda_inname = "wrffdda_d<domain>",
gfdda_interval_m = 360, 360, 360,
gfdda_end_h = 144, 144, 144,
fgdt = 0, 0, 0,
if_no_pbl_nudging_uv = 1, 1, 1,
if_no_pbl_nudging_t = 1, 1, 1,
if_no_pbl_nudging_q = 1, 1, 1,
guv = 0.0003, 0.0003, 0.0003,
gt = 0.0003, 0.0003, 0.0003,
gq = 0.0003, 0.0003, 0.0003,
if_ramping = 1,
dtramp_min = 60.0,
io_form_gfdda = 2,
/
Domain 1: Current date being processed: 1990-08-01_18:00:00.0000, which is loop # 2 out of 29
configflags%julyr, %julday, %gmt: 1990 213 18.0000000
metgrid input_wrf.F first_date_input = 1990-08-01_18:00:00
metgrid input_wrf.F first_date_nml = 1990-08-01_12:00:00
d01 1990-08-01_18:00:00 Timing for input 0 s.
d01 1990-08-01_18:00:00 flag_soil_layers read from met_em file is 1
None of the following are allowed to be TRUE :
( internal_time_loop .EQ. 1 ) F
( config_flags%grid_fdda .NE. 0 ) T
( config_flags%sst_update .EQ. 1 ) F
( config_flags%all_ic_times ) F
( config_flags%smooth_cg_topo ) F
( config_flags%map_proj .EQ. PROJ_CASSINI ) F
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 338
Problems, we cannot have excluded middle data from WPS
The central portion of my 'namelist.input' file follows. Basically, everything works fine until I alter 'grid_fdda' from '0' to '1'. Any assistance will be gratefully appreciated.
BJ
&domains
time_step = 120,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 40, 91, 181,
e_sn = 40, 91, 217,
e_vert = 31, 31, 31,
p_top_requested = 10000,
num_metgrid_levels = 30,
num_metgrid_soil_levels = 4,
dx = 36000, 12000, 4000,
dy = 36000, 12000, 4000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 6, 9,
j_parent_start = 1, 9, 6,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
eta_levels = 1.0000, 0.9958, 0.9907, 0.9846, 0.9774,
0.9688, 0.9585, 0.9463, 0.9319, 0.9148,
0.8946, 0.8709, 0.8431, 0.8107, 0.7733,
0.7301, 0.6809, 0.6254, 0.5635, 0.4954,
0.4218, 0.3536, 0.2930, 0.2394, 0.1917,
0.1495, 0.1120, 0.0788, 0.0493, 0.0232,
0.00000,
/
&physics
mp_physics = 6, 6, 6,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 5, 5, 5,
sf_urban_physics = 0, 0, 0,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 1, 1, 1,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 5, 5, 5,
cudt = 5, 5, 5,
cu_rad_feedback = .true., .true., .true.,
sst_update = 0,
usemonalb = .true.,
rdmaxalb = .true.,
surface_input_source = 3,
num_soil_layers = 4,
num_land_cat = 21,
/
&fdda
grid_fdda = 1, 1, 0,
gfdda_inname = "wrffdda_d<domain>",
gfdda_interval_m = 360, 360, 360,
gfdda_end_h = 144, 144, 144,
fgdt = 0, 0, 0,
if_no_pbl_nudging_uv = 1, 1, 1,
if_no_pbl_nudging_t = 1, 1, 1,
if_no_pbl_nudging_q = 1, 1, 1,
guv = 0.0003, 0.0003, 0.0003,
gt = 0.0003, 0.0003, 0.0003,
gq = 0.0003, 0.0003, 0.0003,
if_ramping = 1,
dtramp_min = 60.0,
io_form_gfdda = 2,
/