[Wrf-users] Model is crashing with CAM shortwave scheme

Thu Mar 17 19:50:09 MDT 2011

Hi All,

After finishing a successful run with Dudhia shortwave scheme I tried  
to use CAM shortwave scheme but could not succeed. I was wondering  
whether any one faced the same problem while using CAM radiation  
scheme. I get the following error message:

taskid: 0 hostname: bro95
Namelist fdda not found in namelist.input. Using registry defaults for  
variable
s in fdda
Namelist dfi_control not found in namelist.input. Using registry  
defaults for v
ariables 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  
variable
s in fire
Ntasks in X 2, ntasks in Y 4
NOTE: CAM radiation is in use, setting: paerlev=29, levsiz=59, cam_abs_dim1=4,
cam_abs_dim2=e_vert
NOTE: num_soil_layers has been set to 4
WRF V3.2 MODEL
wrf: calling alloc_and_configure_domain
*************************************
Parent domain
ids,ide,jds,jde 1 225 1 175
ims,ime,jms,jme -4 119 -4 51
ips,ipe,jps,jpe 1 112 1 44
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
alloc_space_field: domain 1, 184262200 bytes allocated
wrf: calling model_to_grid_config_rec
wrf: calling set_scalar_indices_from_config
wrf: calling init_wrfio
Entering ext_gr1_ioinit
setup_timekeeping: set xtime to 0.0000000E+00
setup_timekeeping: set julian to 172.0000
setup_timekeeping: returning...
wrf main: calling open_r_dataset for wrfinput
med_initialdata_input: calling input_input
input_wrf: set xtime to 0.0000000E+00
mminlu = 'USGS'
med_initialdata_input: back from input_input
checking boundary conditions for grid
boundary conditions OK for grid
start_domain_em: Before call to phy_init
WRF TILE 1 IS 1 IE 112 JS 1 JE 44
set_tiles3: NUMBER OF TILES = 1
top of phy_init
phy_init: start_of_simulation = T
calling nl_get_iswater, nl_get_isice, nl_get_mminlu_loc
after nl_get_iswater, nl_get_isice, nl_get_mminlu_loc
top of landuse_init
INPUT LandUse = "USGS"
returning from of landuse_init
reading CAM_ABS_DATA
num_months = 13
reading CAM_AEROPT_DATA
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
wrf.exe 00000000014217C1 Unknown Unknown Unknown
wrf.exe 0000000001423821 Unknown Unknown Unknown
wrf.exe 00000000014225FA Unknown Unknown Unknown
wrf.exe 00000000011A6D5A Unknown Unknown Unknown
wrf.exe 0000000000D2DB60 Unknown Unknown Unknown
wrf.exe 0000000000D27583 Unknown Unknown Unknown
wrf.exe 0000000000A8B6F7 Unknown Unknown Unknown
wrf.exe 00000000009CA6AE Unknown Unknown Unknown
wrf.exe 00000000009C82E7 Unknown Unknown Unknown
wrf.exe 000000000047FBC2 Unknown Unknown Unknown
wrf.exe 000000000047F009 Unknown Unknown Unknown
wrf.exe 000000000047EFAC Unknown Unknown Unknown
libc.so.6 00002AE11F2F9994 Unknown Unknown Unknown
wrf.exe 000000000047EEB9 Unknown Unknown Unknown

############################################################################
### The namelist.input file I am using is given below ##############
############################################################################
&time_control
run_days = 5,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2008, 2008, 2008,
start_month = 06, 06, 06,
start_day = 21, 21, 21,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2008, 2008, 2008,
end_month = 06, 06, 06,
end_day = 29, 29, 29,
end_hour = 00, 00, 00,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 10800
input_from_file = .true.,.true.,.true.,
history_interval = 180, 180, 180,
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 = 100
/

&domains
time_step = 48,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 225, 202, 181,
e_sn = 175, 202, 271,
e_vert = 28, 28, 28,
p_top_requested = 10000,
num_metgrid_levels = 30
num_metgrid_soil_levels = 4,
dx = 8000, 2666.666, 888.888,
dy = 8000, 2666.666, 888.888,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 75, 70,
j_parent_start = 1, 65, 50,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 6, 6, 6,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 3, 3, 3,
radt = 8, 8, 8,
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, 0, 0,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
/
&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,
/

############################################################################
### The namelist.wps file I am using is given below ##############
############################################################################

&share
wrf_core = 'ARW',
max_dom = 3,

start_date = '2008-06-21_00:00:00','2008-06-21_00:00:00',  
'2008-06-21_00:00:00',
end_date = '2008-06-29_00:00:00','2008-06-21_00:00:00', '2008-06-21_00:00:00',

interval_seconds = 10800
io_form_geogrid = 2,
/

&geogrid
parent_id = 1, 1, 2,
parent_grid_ratio = 1, 3, 3,
i_parent_start = 1, 75, 70,
j_parent_start = 1, 65, 50,
e_we = 225, 202, 181,
e_sn = 175, 202, 271,
geog_data_res = ?modis_30s+5m?, ?modis_30s+2m?, ?modis_30s+30s?

dx = 8000,
dy = 8000,
map_proj = 'lambert',
ref_lat = 43.0,
ref_lon = -80.0,
truelat1 = 30.0,
truelat2 = 60.0,
stand_lon = -80.0,
geog_data_path = '/work/kamal/brown/DATA/geog'
/

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

&metgrid
fg_name = 'NARR'
io_form_metgrid = 2,
constants_name ='./NARRFIX:1979-11-08_00'
/

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

So, it only seems to be a problem (for me) when I use CAM shortwave  
scheme instead of Dudhia scheme. Does anyone know the reason for this ?

Thanks
Kamal