wrf-chem

با سلام

من مدل را برای 48 ساعت اجرا کردم. اما همچنان خروجی مدل برای PM10 صفر بود.

من مدل با شرایط زیر اجرا کردم. لطفا اگر اشتباهی در کار من وجود دارد تذکر دهید. متشکرم

-کامپایل WPS با دربر گرفتن داده های EROD و CLARFRAC و  SANDFRAC 

-تنظیم نیملیست با مشخصات زیر برای اجرای مدل ورف-کم تنها برای گردوغبار

 &time_control
 run_days                            = 0,
 run_hours                           = 48,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2016, 2016, 2000,
 start_month                         = 06,   06,   01,
 start_day                           = 15,   15,   24,
 start_hour                          = 00,   00,   12,
 start_minute                        = 00,   00,   00,
 start_second                        = 00,   00,   00,
 end_year                            = 2016, 2016, 2000,
 end_month                           = 06,   06,   01,
 end_day                             = 17,   17,   25,
 end_hour                            = 00,   00,   12,
 end_minute                          = 00,   00,   00,
 end_second                          = 00,   00,   00,
 interval_seconds                    = 10800,
 input_from_file                     = .true.,.true.,.false.,
 history_interval                    = 60,  60,   60,
 frames_per_outfile                  = 72, 1000, 1000,
 restart                             = .false.,
 restart_interval                    = 864000,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 auxinput1_inname                    = "met_em.d<domain>.<date>"
 debug_level                         = 00
 /

 &domains
 time_step                           = 300,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 s_we                                = 1,     1,     1,
 e_we                                =121,    112,   94,
 s_sn                                = 1,     1,     1,
 e_sn                                = 61,    97,    91,
 s_vert                              = 1,     1,     1,
 e_vert                              = 41,    28,    28,
 num_metgrid_levels                  = 32
 dx                                  = 60000,  3333,  1111,
 dy                                  = 60000,  3333,  1111,
 grid_id                             = 1,     2,     3,
 parent_id                           = 0,     1,     2,
 i_parent_start                      = 0,     30,    30,
 j_parent_start                      = 0,     20,    30,
 parent_grid_ratio                   = 1,     3,     3,
 parent_time_step_ratio              = 1,     3,     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                          = 4,     2,     2,
 progn                               = 0,
 ra_lw_physics                       = 1,     1,     1,
 ra_sw_physics                       = 2,     1,     1,
 radt                                = 30,    10,    10,
 sf_sfclay_physics                   = 1,     1,     1,
 sf_surface_physics                  = 2,     2,     2,
 bl_pbl_physics                      = 1,     1,     1,
 bldt                                = 1,     0,     0,
 cu_physics                          = 5,     5,     0,
 cu_diag                             = 1,     0,     0,
 cudt                                = 1,     1,     5,
 ishallow                            = 0,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 surface_input_source                = 1,
 num_soil_layers                     = 4,
 sf_urban_physics                    = 0,
 mp_zero_out                         = 2,
 mp_zero_out_thresh                  = 1.e-8,
 maxiens                             = 1,
 maxens                              = 3,
 maxens2                             = 3,
 maxens3                             = 16,
 ensdim                              = 144,
 cu_rad_feedback                     = .true.,
 /

 &fdda
 /

 &dfi_control
 /

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

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

 &chem
 kemit                               = 40,
 chem_opt                            = 401,
 bioemdt                             = 0,
 photdt                              = 0,
 chemdt                              = 0,
 io_style_emissions                  = 1,  
 emiss_opt                           = 0,          3,
 chem_in_opt                         = 0,          0,
 phot_opt                            = 0,          0,
 gas_drydep_opt                      = 0,          1, 
 aer_drydep_opt                      = 0,          1, 
 bio_emiss_opt                       = 0,          0,
 dust_opt                            = 1,
 dmsemis_opt                         = 0,
 seas_opt                            = 0,
 gas_bc_opt                          = 0,          0,
 gas_ic_opt                          = 0,          0,
 aer_bc_opt                          = 0,          0,
 aer_ic_opt                          = 0,          0,
 gaschem_onoff                       = 0,          0,
 aerchem_onoff                       = 0,          0,
 wetscav_onoff                       = 0,          0,
 cldchem_onoff                       = 0,          0,
 vertmix_onoff                       = 1,          0,
 chem_conv_tr                        = 1,          0,
 biomass_burn_opt                    = 0,          0,
 plumerisefire_frq                   = 30,         0,
 have_bcs_chem                       = .false., .false., .false.,
 aer_ra_feedback                     = 0,
 aer_op_opt                          = 0,
 opt_pars_out                        = 0,
 diagnostic_chem                     = 0,
 /

 /

real.exe و wrf.exe بدرستی اجرا می شود. و فایل های خروجی ایجاد می گردد.

اما برای  خروجی PM10 نشان داده نمی شود.

در نیملیست من dust_opt را ۱ انتخاب کرده ام. که در راهنمای برنامه مربوط به طرحواره GOCART است. آیا برای استفاده از این طرحواره داده خاصی می بایست اضافه شود؟

با سپاس فراوان