[Wrf-users] how to choose physic options

[Hosein S]

Dear WRF users

I want to run WRF for my special domain for air quality modeling using 3 nesting domains, 9, 3 and 1 km.
Can anyone help me how to choose the best physic options for my domains?

Thanks in advance
Hossein

Here in the namelist i use:

 &time_control
 run_days                            = 11,
 run_hours                           = 0,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2012, 2012, 2012,
 start_month                         = 07,   07,   07,
 start_day                           = 05,   05,   05,
 start_hour                          = 00,   00,   00,
 start_minute                        = 00,   00,   00,
 start_second                        = 00,   00,   00,
 end_year                            = 2012, 2012, 2012,
 end_month                           = 07,   07,   07,
 end_day                             = 15,   15,   15,
 end_hour                            = 00,   00,   00,
 end_minute                          = 00,   00,   00,
 end_second                          = 00,   00,   00,
 interval_seconds                    = 21600
 input_from_file                     = .true.,.true.,.true.,
 history_interval                    = 60,  60,   60,
 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                           = 40,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 3,
 e_we                                = 179,    145,    91,   
 e_sn                                = 173,    130,    82,   
 e_vert                              = 28,    28,    28,   
 p_top_requested                     = 5000,
 num_metgrid_levels                  = 27,
 num_metgrid_soil_levels             = 4,
 dx                                  = 9000,  3000,  1000,
 dy                                  = 9000,  3000,  1000,
 grid_id                             = 1,     2,     3,    
 parent_id                           = 0,     1,     2,    
 i_parent_start                      = 1,     66,    59,    
 j_parent_start                      = 1,     65,    51,    
 parent_grid_ratio                   = 1,     3,     3,    
 parent_time_step_ratio              = 1,     3,     3,    
 feedback                            = 1,
 smooth_option                       = 0
 /

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

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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20130415/002857b3/attachment-0001.html