Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

NUDAPT and ucm=3

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

johnfmejia

New member
I am using nlcd2011_9s and NUDAPT_1km data for an urban simulation inside the US and using WPS/WRF 4.3.1. However, I cannot pass from a singular matrix error as follows:

FATAL CALLED FROM FILE: <stdin> LINE: 2462
singular matrix in gaussjbem

Does anyone have a similar problem or insight about how to solve it?
 
Ming Chen said:
Can you upload your namelist.wps and namelist.input for me to take a look? Thanks.

Hi Ming: Thanks you very much for trying to understand the problem.

Below are my namelist.wps and namelist.input files. I am using NLCD and NUDAPT data for urban vertical distribution. I have tried WRF4.3 and WRF.4.3.1, different domain sizes and moved the domains to have them all inside continental U.S. The model works fine for UCM =1, 2. For UCM=3, real.exe is fine but wrf.exe crashes when starting domain 2 (error: FATAL CALLED FROM FILE: <stdin> LINE: 2462
singular matrix in gaussjbem). Any insight is appreciated.

Best,



&share
wrf_core = 'ARW'
max_dom = 4
start_date = '2020-08-06_12:00:00', '2020-08-06_12:00:00', '2020-08-06_12:00:00'
end_date = '2020-08-08_12:00:00', '2020-08-08_12:00:00', '2020-08-08_12:00:00'
active_grid = .true., .true., .true.,
interval_seconds = 21600
io_form_geogrid = 2
opt_output_from_geogrid_path = './'
debug_level = 1
/

&geogrid
parent_id = 1,1,2,
parent_grid_ratio = 1,3,3,
i_parent_start = 1,16,45,
j_parent_start = 1,14,33,
e_we = 89,175,268,
e_sn = 92,166,256,
geog_data_res = 'nlcd2011_9s', 'nlcd2011_9s', 'nlcd2011_9s',
dx = 9000
dy = 9000
map_proj = 'lambert'
ref_lat = 30.406,
ref_lon = -95.208,
truelat1 = 30.406,
truelat2 = 30.406,
stand_lon = -95.208,
ref_x = 44.5,
ref_y = 46.0,
geog_data_path = './WPS_GEOG/'
opt_geogrid_tbl_path = './geogrid/'
/

&ungrib
/

&metgrid
fg_name = 'FILE',
io_form_metgrid = 2,
opt_output_from_metgrid_path = './'
opt_metgrid_tbl_path = './'
/

&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
/



&time_control
run_days = 0,
run_hours = 00,
run_minutes = 0,
run_seconds = 0,
start_year=2020,2020,2020
start_month=08,08,08
start_day=06,06,06
start_hour=12,12,12
start_minute = 00,00,00
start_second = 00,00,00,
end_year=2020,2020,2020,
end_month=08,08,08,
end_day=08,08,08,
end_hour=12,12,12,
end_minute = 00,00,00,
end_second = 00,00,00,
interval_seconds = 21600,
input_from_file = .true.,.true.,.true.,
history_interval = 60, 30, 30,
frames_per_outfile = 1,1,1,
restart = .false.,
restart_interval = 360,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
auxinput4_inname = "wrflowinp_d<domain>",
auxinput4_interval = 360, 360, 360,
io_form_auxinput4 = 2
/

&domains
time_step = 45,
use_adaptive_time_step = .true.,
step_to_output_time = .true.,
target_cfl = 2,2,2,
min_time_step = -1,-1,-1,
max_dom = 3,
e_we = 89,175,268,
e_sn = 92,166,256,
e_vert = 65, 65, 65,
num_metgrid_levels = 34,
num_metgrid_soil_levels = 4,
dx = 9000, 3000, 1000,
dy = 9000, 3000, 1000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1,16,45,
j_parent_start = 1,14,33,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 0,
smooth_option = 1,
p_top_requested = 5000.0
/

&physics
mp_physics = 0,0,0,
ra_lw_physics = 1,1,1,
ra_sw_physics = 1,1,1,
radt = 10,10,10,
sf_sfclay_physics = 1,1,1,
sf_surface_physics = 2,2,2,
bl_pbl_physics = 1,1,1,
bldt = 0,0,0,
cu_physics = 0,0,0,
cudt = 0,0,0,
sf_urban_physics = 3,3,3,
ifsnow = 0,
icloud = 0,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 40
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
sst_update = 1,
tmn_update = 0,
sst_skin = 1,
usemonalb = .true.,
/


&fdda
grid_fdda = 2, 0, 0,
gfdda_inname = "wrffdda_d<domain>",
gfdda_end_h = 1007400, 24, 24,
gfdda_interval_m = 360, 360, 360,
fgdt = 0, 0, 0,
fgdtzero = 1, 0, 0,
if_no_pbl_nudging_uv = 1, 0, 0,
if_no_pbl_nudging_t = 1, 0, 0,
if_no_pbl_nudging_ph = 1, 0, 0,
if_no_pbl_nudging_q = 1, 0, 0,
if_zfac_uv = 1, 0, 0,
k_zfac_uv = 20, 10, 10,
if_zfac_t = 1, 0, 0,
k_zfac_t = 20, 10, 10,
if_zfac_ph = 1, 0, 0,
k_zfac_ph = 20, 10, 10,
if_zfac_q = 1, 0, 0,
k_zfac_q = 20, 10, 10,
dk_zfac_uv = 5, 1, 1,
dk_zfac_t = 5, 1, 1,
dk_zfac_ph = 5, 1, 1,
guv = 0.0003, 0.0003, 0.0003,
gt = 0.0003, 0.0003, 0.0003,
gph = 0.0003, 0.0003, 0.0003,
gq = 0.0003, 0.0003, 0.0003,
xwavenum = 3
ywavenum = 3
if_ramping = 0,
dtramp_min = 60.0,
io_form_gfdda = 2,
/

&dynamics
time_step_sound = 4,
rk_ord = 3,
w_damping = 1,
km_opt = 2,2,2,
gwd_opt = 1,
diff_opt = 1,1,1,
damp_opt = 3,
base_temp = 290.
zdamp = 1000., 1000., 1000.,
dampcoef = 0.2, 0.2, 0.2,
khdif = 0,0,0,
kvdif = 0,0,0,
smdiv = 0.1,0.1,0.1,
emdiv = 0.01,0.01,0.01,
epssm = 0.2,0.2,0.2,
non_hydrostatic = .true., .true., .true.,
time_step_sound = 0,0,0,
h_mom_adv_order = 5,5,5,
v_mom_adv_order = 3,3,3,
h_sca_adv_order = 3,5,5,
v_sca_adv_order = 2,3,3,
moist_adv_opt = 1,1,1,
scalar_adv_opt = 1,1,1,
tke_adv_opt = 1,1,1,
/


&bdy_control
spec_bdy_width = 10,
spec_zone = 1,
relax_zone = 9,
specified = .true., .false., .false.,
nested = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
 
Hi,
I can repeat the error. Unfortunately we don't have a solution at present. I will talk to our experts and get back to you once I hear from them.
 
Top