[Wrf-users] how to run a multi_days simulation including wildfire emission

tanxiangshan2005 tanxiangshan2005 at 126.com
Tue May 23 02:07:37 MDT 2017


hello , 
       I want to run wrfchem including wildfire emission , and this  fire emission created by prep_chem_src. 
       Now I meet some problem. 
       if the fire emisison is  named wrffirechemi_d01, the wrfchem can run successfully. but if the fire emission is named  wrffirechemi_d01_2017-04-25_00:00:00, when run ./real.exe , the error will come "Could not find matching time in input file wrffirechemi_d01" 
  
       the reason  why I change wrffirechemi_d01 to  wrffirechemi_d01_2017-04-25_00:00:00, is that :firstly ,in prep_chem_src user's guide , said
Rename output to include the date:
→ mv wrfchemi_d01 wrfchemi_d01_{$year}-$mo-${da}_00:00:00
→ mv wrffirechemi_d01 wrffirechemi_d01_{$year}-$mo-${da}_00:00:00
→ mv wrfchemi_gocart_bg_d01 wrfchemi_gocart_bg_d01_{$year}-$mo-$
{da}_00:00:00  
second, only in wrffirechemi_d01_2017-04-25_00:00:00 form , I can run  multi_days simulation. 
     and my namelist is :
     &time_control
 run_days                            = 6,
 run_hours                           = 0,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2017, 2015, 2013, 2013,
 start_month                         = 04,   11,   12,   12,
 start_day                           = 25,   01,   08,   08,
 start_hour                          = 00,   12,   12,   12,
 start_minute                        = 00,   00,   00,  00, 
 start_second                        = 00,   00,   00,  00,
 end_year                            = 2017, 2015, 2013, 2013,
 end_month                           = 04,   11,   12,   12,
 end_day                             = 26,   30,   14,   14,
 end_hour                            = 00,   12,   12,   12,
 end_minute                          = 00,   00,   00,  00,
 end_second                          = 00,   00,   00,  00,
 interval_seconds                    = 21600,
 input_from_file                     = .true.,.true.,.true.,.true.,
 history_interval                    = 60,  60,   60, 60,
 frames_per_outfile                  = 1,   24,    1,  1,
 restart                             = .false.,
 restart_interval                    = 72000,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 io_form_auxinput4                   = 0
 io_form_auxinput5                   = 2
 io_form_auxinput7                   = 2
 io_form_auxinput12                  = 0
 frames_per_auxinput5                = 1, 1, 1, 1, 
 frames_per_auxinput7                = 1, 1, 1, 1,   
 auxinput5_inname                    = 'wrfchemi_d<domain>_<date>'
 auxinput7_inname                    = 'wrffirechemi_d<domain>_<date>'
 auxinput5_interval_m                = 1440, 60, 60 
 auxinput7_interval_m                = 1440, 60, 60 
 debug_level                         = 1000 
 /
 &domains
 time_step                           = 180,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 s_we                                = 1,     1,     1,    1,
 s_sn                                = 1,     1,     1,    1,
 e_we                                = 183, 244, 121, 151,
 e_sn                                = 173, 274, 121, 151, 
 dx                                  = 45000,15000,5000,0,
 dy                                  = 45000,15000,5000,0,
 s_vert                              = 1,     1,     1,    1,
 e_vert                              = 30,    30,    35,   35,
 num_metgrid_levels                  = 32,
 num_metgrid_soil_levels             = 4,
 p_top_requested      = 5000
 grid_id                             = 1,     2,     3,    4,
 parent_id                           = 1,     1,     2,    3,
 i_parent_start                      = 1,  68, 49, 40,
 j_parent_start                      = 1,  71, 50, 16,
 parent_grid_ratio                   = 1,     3,     3,    3,
 parent_time_step_ratio              = 1,     3,     3,    3,
 feedback                            = 1,
 smooth_option                       = 2
 /
 &physics
 mp_physics                          = 6,     6,     2,
 progn                               = 0,     0,     0,
 naer                                = 1e9
 ra_lw_physics                       = 4,     4,     1,
 ra_sw_physics                       = 4,     4,     2,
 radt                                = 30,    30,    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                          = 5,     5,     0,
 cugd_avedx                          = 3,
 cudt                                = 0,     0,     1,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 surface_input_source                = 1,
 num_soil_layers                     = 4,
 num_land_cat                        = 20,
 sf_urban_physics                    = 0,     0,     0,
 maxiens                             = 1,
 maxens                              = 3,
 maxens2                             = 3,
 maxens3                             = 16,
 ensdim                              = 144,
 cu_rad_feedback                     = .true.,
 cu_diag                             = 1,
 sst_update                          = 0,
 /
 &fdda
 grid_fdda                           = 0,
 gfdda_inname                        = "wrffdda_d<domain>",
 gfdda_end_h                         = 144,
 gfdda_interval_m                    = 360,
 fgdt                                = 0,
 if_no_pbl_nudging_uv                = 1,
 if_no_pbl_nudging_t                 = 1,
 if_no_pbl_nudging_q                 = 1,
 if_zfac_uv                          = 1,
  k_zfac_uv                          = 10,
 if_zfac_t                           = 1,
  k_zfac_t                           = 10,
 if_zfac_q                           = 1,
  k_zfac_q                           = 10,
 guv                                 = 0.0003,
 gt                                  = 0.0003,
 gq                                  = 0.0003,
 if_ramping                          = 1,
 dtramp_min                          = 60.0,
 io_form_gfdda                       = 2,
 /
 &dynamics
 w_damping                           = 1,
 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                       = 2,      1,      1,
 scalar_adv_opt                      = 2,      1,      1,
 chem_adv_opt                        = 2,      1,      1,
 /
 &bdy_control
 spec_bdy_width                      = 5,
 spec_zone                           = 1,
 relax_zone                          = 4,
 specified                           = .true., .false.,.false.,
 nested                              = .false., .true., .true.,
 /
 &grib2
 /
 &chem
 kemit                               = 1,
 chem_opt                            = 301,        2,
 bioemdt                             = 30,       30,
 photdt                              = 30,       30,
 chemdt                              = 10,       10,
 io_style_emissions                  = 1,
 emiss_inpt_opt                      = 1,        1,
 emiss_opt                           = 5,        3,
 chem_in_opt                         = 0,        0,
 phot_opt                            = 1,        1,
 gas_drydep_opt                      = 1,        1,
 aer_drydep_opt                      = 1,        1,
 bio_emiss_opt                       = 1,        1,
 dust_opt                            = 0,
 dmsemis_opt                         = 0,
 seas_opt                            = 0,
 gas_bc_opt                          = 1,        1,
 gas_ic_opt                          = 1,        1,
 aer_bc_opt                          = 1,        1,
 aer_ic_opt                          = 1,        1,
 gaschem_onoff                       = 1,        1,
 aerchem_onoff                       = 1,        1,
 wetscav_onoff                       = 0,        0,
 cldchem_onoff                       = 0,        0,
 vertmix_onoff                       = 1,        1,
 chem_conv_tr                        = 1,        1,
 biomass_burn_opt                    = 1,        0,
 plumerisefire_frq                   = 120,       30,
 aer_ra_feedback                     = 0,        0,
 have_bcs_chem                       = .false., .false.,
 /
 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /

       And how can I solve this problem, please help. Thank you very much.
Bests,
Jingjing Shang 





2017-05-23
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20170523/3a40ff37/attachment-0001.html 


More information about the Wrf-users mailing list