WRFChem SST_update and feedback woes

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

WRFChem SST_update and feedback woes

Postby sru20 » Thu Nov 17, 2011 4:36 pm

Hello,
I am having problems with my model runs. I am running a 2 domain step up - outer domain 15km and inner domain 5 km. I have successfully generated wrinput_d01, wrfinput_d02, wrfbdy_d01, wrfloinp_d01, wrfloinp_d02 files and i need to update SST (with SST_update=1). Also there is feedback between the domains (feedback=1).
Thing is that i can run the model successully without feedback but as soon as i switch on feedback, the model crashes upon starting the inner domain which starts several hours after the outer domain. Also the model runs OK (with feedback) when i set sst_update=0.The times and the values of SST, Albedo and vegefra all look reasonable in the wrfloinp files.

Code: Select all
&time_control
run_days                            = 0,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2006, 2006,
start_month                         = 06,     06,
start_day                           = 21,     22,
start_hour                          = 18,     00,
start_minute                        = 00,     00,
start_second                        = 00,     00,
end_year                            = 2006, 2006,
end_month                           = 06,     06,
end_day                             = 24,     24,
end_hour                            = 23,     00,
end_minute                          = 00,     00,
end_second                          = 00,     00,
interval_seconds                    = 21600,
input_from_file                     = .true.,.true.,
fine_input_stream                   = 0, 2,
history_interval                    = 60,  60,
frames_per_outfile                  = 1, 1,
restart                             = .false.,
restart_interval                    = 360,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,
auxinput1_inname                    = "met_em.d<domain>.<date>",
auxinput4_inname                    = "wrflowinp_d<domain>",
auxinput6_inname                    = 'wrfbiochemi_d01',
auxinput7_inname                    = 'wrffirechemi_d<domain>',
auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',
auxinput5_interval_m                = 60, 60,
auxinput7_interval_m                = 86400, 86400,
auxinput8_interval_m                = 86400, 86400,
auxinput4_interval                  = 360,360,
io_form_auxinput2                   = 2,
io_form_auxinput4                   = 2,
io_form_auxinput5                   = 2,
io_form_auxinput6                   = 0,
io_form_auxinput7                   = 0,
debug_level                         = 500,
/

&dfi_control
/

&domains
time_step                = 90,
time_step_fract_num      = 0,
time_step_fract_den      = 1,
max_dom                  = 2,
s_we                     = 1,        1,
e_we                     = 134,      184,
s_sn                     = 1,        1,
e_sn                     = 146,      223,
s_vert                   = 1,        1,
e_vert                   = 28,       28,
num_metgrid_levels       = 17,
dx                       = 15000,     5000,
dy                       = 15000,     5000,
grid_id                  = 1,        2,
parent_id                = 1,        1,
i_parent_start           = 1,       29,
j_parent_start           = 1,       41,
parent_grid_ratio        = 1,        3,
parent_time_step_ratio   = 1,        3,
feedback                 = 1,
smooth_option            = 0,
p_top_requested          = 5000,
zap_close_levels         = 50,
interp_type              = 1,
t_extrap_type            = 2,
force_sfc_in_vinterp     = 0,
use_levels_below_ground  = .true.,
use_surface              = .true.,
lagrange_order           = 1,
sfcp_to_sfcp             = .true.,
/

&physics
mp_physics                          = 2,     2,
progn                               = 1,     1,
ra_lw_physics                       = 4,     4,
ra_sw_physics                       = 2,     2,
radt                                = 5,     5,
sf_sfclay_physics                   = 2,     2,
sf_surface_physics                  = 2,     2,
bl_pbl_physics                      = 2,     2,
bldt                                = 0,     0,
cu_physics                          = 5,     5,
cugd_avedx                          = 1,
cudt                                = 0,     0,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 1,     1,
mp_zero_out                         = 2,
mp_zero_out_thresh                  = 1.e-12
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
cu_rad_feedback                     = .true.,
sst_update                          = 1,
sst_skin                            = 1,
/

&fdda
/

&Options for adaptive time step
use_adaptive_time_step              = .true.,
step_to_output_time                 = .true.,
target_cfl                          = 1.2, 1.2,
max_step_increase_pct               = 5, 51,
starting_time_step                  = -1, -1,
max_time_step                       = -1, -1,
min_time_step                       = -1, -1,
adaptation_domain                   = 1,
/

&dynamics
rk_ord                              = 3,
w_damping                           = 1,
diff_opt                            = 1,
km_opt                              = 4,
base_temp                           = 290.
damp_opt                            = 0,
zdamp                               = 5000.,  5000.,
dampcoef                            = 0.01,    0.01,
diff_6th_opt                        = 0,
diff_6th_factor                     = 0.12,
khdif                               = 0,      0,
kvdif                               = 0,      0,
non_hydrostatic                     = .true., .true.,
moist_adv_opt                       = 2,      2,
scalar_adv_opt                      = 2,      2,
chem_adv_opt                        = 2,      2,
tke_adv_opt                         = 2,      2,
time_step_sound                     = 4,      4,
h_mom_adv_order                     = 5,      5,
v_mom_adv_order                     = 3,      3,
h_sca_adv_order                     = 5,      5,
v_sca_adv_order                     = 3,      3,
/

&bdy_control
spec_bdy_width           = 5,
spec_zone                = 1,
relax_zone               = 4,
specified                = .true.,   .false.,
periodic_x               = .false.,  .false.,
symmetric_xs             = .false.,  .false.,
symmetric_xe             = .false.,  .false.,
open_xs                  = .false.,  .false.,
open_xe                  = .false.,  .false.,
periodic_y               = .false.,  .false.,
symmetric_ys             = .false.,  .false.,
symmetric_ye             = .false.,  .false.,
open_ys                  = .false.,  .false.,
open_ye                  = .false.,  .false.,
nested                   = .false.,  .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

&chem
kemit                               = 5,
chem_opt                            = 402, 402,
bioemdt                             = 30, 6,
photdt                              = 30, 30,
chemdt                              = 4, 1.33,
frames_per_emissfile                = 1,
io_style_emissions                  = 2,
emiss_inpt_opt                      = 121,121,
emiss_opt                           = 10, 10,
chem_in_opt                         = 0,  0,
phot_opt                            = 2, 2,
gas_drydep_opt                      = 1, 1,
aer_drydep_opt                      = 1, 1,
bio_emiss_opt                       = 1, 1,
gas_bc_opt                          = 0, 0,
gas_ic_opt                          = 0, 0,
aer_bc_opt                          = 0, 0,
aer_ic_opt                          = 0, 0,
gaschem_onoff                       = 1, 1,
aerchem_onoff                       = 1, 1,
wetscav_onoff                       = 1, 1,
cldchem_onoff                       = 1, 1,
vertmix_onoff                       = 1, 1,
chem_conv_tr                        = 1, 1,
seas_opt                            = 2,
dust_opt                            = 0,
dmsemis_opt                         = 2,
biomass_burn_opt                    = 0,0,
plumerisefire_frq                   = 30,30,
have_bcs_chem                       = .true., .true.,
aer_ra_feedback                     = 1,1,
diagnostic_chem                     = 0,0,
/
                                                         

Has anybody experienced this problem? Please help!

Steve
sru20
 
Posts: 6
Joined: Sat May 01, 2010 5:49 am

Return to Nesting

Who is online

Users browsing this forum: No registered users and 2 guests