Problem of AOD results from WRF-chem output

Any issues with running WRF Chemistry

Problem of AOD results from WRF-chem output

Postby Stella@88 » Sun Sep 29, 2019 4:31 pm

Hello members,

I am running the WRF-Chem V3.5 using the RACM Chemistry and MADE/VBS aerosols (chem_opt=108). However, all values of EXTCOF55 are 'nan' after finishing the simulation without any error message. I am not sure which part goes wrong. Would anyone please kindly give me some suggestions on this issue?

Here is the namelist file:
Code: Select all
 &time_control
 run_days                            = 36,
 run_hours                           = 0,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2011,
 start_month                         = 04,
 start_day                           = 28,
 start_hour                          = 00,
 start_minute                        = 00,
 start_second                        = 00,
 end_year                            = 2011,
 end_month                           = 06,
 end_day                             = 02,
 end_hour                            = 00,
 end_minute                          = 00,
 end_second                          = 00,
 interval_seconds                    = 21600,
 input_from_file                     = T,
 history_interval                    = 60,
 frames_per_outfile                  = 24,
 restart                             = F,
 cycling                             = F,
 restart_interval                    = 1440,
 io_form_history                     = 2,
 io_form_restart                     = 2,
 io_form_input                       = 2,
 io_form_boundary                    = 2,
 debug_level                         = 0,
 diag_print                          = 0,
 all_ic_times                        = .false.,
 adjust_output_times                 = .false.,
 override_restart_timers             = .false.,
 write_hist_at_0h_rst                = .false.,
 auxinput1_inname                    = "met_em.d<domain>.<date>"
 auxinput4_inname                    = "wrflowinp_d<domain>"
 auxinput4_interval                  = 360
 io_form_auxinput4                   = 2
 auxinput7_inname                    = "wrffirechemi_d<domain>_<date>"
 io_form_auxinput7                   = 2
 frames_per_auxinput7                = 1
 auxinput7_interval_m                = 60
 auxinput6_inname                    = 'wrfbiochemi_d01'
 io_form_auxinput6                   = 2
! auxinput6_interval_h                = 24
 auxinput5_interval_m                = 60
 io_form_auxinput5                   = 2
/
 &domains
 time_step                           = 120,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 time_step_dfi                       = 60,
 max_dom                             = 1,
 s_we                                = 1,
 e_we                                = 208,
 s_sn                                = 1, 
 e_sn                                = 190,
 s_vert                              = 1, 
 e_vert                              = 31, 
 p_top_requested                     = 5000,
 eta_levels = 1.000, 0.993, 0.980, 0.966, 0.950,
              0.933, 0.913, 0.892, 0.869, 0.844,
              0.816, 0.786, 0.753, 0.718, 0.680,
              0.639, 0.596, 0.550, 0.501, 0.451,
              0.398, 0.345, 0.290, 0.236, 0.188,
              0.145, 0.108, 0.075, 0.046, 0.021,
              0.000
 num_metgrid_levels                  = 27,
 num_metgrid_soil_levels             = 4,
 dx                                  = 18000, 10000,  3333.33,
 dy                                  = 18000, 10000,  3333.33,
 grid_id                             = 1,     2,     3,
 parent_id                           = 0,     1,     2,
 interp_type                         = 2
 extrap_type                         = 2
 t_extrap_type                       = 2
 use_levels_below_ground             = .true.
 use_surface                         = .true.
 lagrange_order                      = 1
 zap_close_levels                    = 500
 lowest_lev_from_sfc                 = .false.
 force_sfc_in_vinterp                = 1
 sfcp_to_sfcp                        = .false.
 smooth_cg_topo                      = .false.
 use_tavg_for_tsk                    = .false.
 aggregate_lu                        = .false.
 /
 &dfi_control
 dfi_opt                             = 0
 dfi_nfilter                         = 7
 /

 &physics
 mp_physics                          = 2,     3,     3,
 ra_lw_physics                       = 4,     1,     1,
 ra_sw_physics                       = 4,     1,     1,
 radt                                = 30,    30,    30,
 sf_sfclay_physics                   = 5,     1,     1,
 sf_surface_physics                  = 2,     2,     2,
 bl_pbl_physics                      = 5,     1,     1,
 bldt                                = 0,     0,     0,
 cu_physics                          = 5,     1,     0,
 cu_diag                             = 1,     1,     0,
 cudt                                = 0,     5,     5,
 isfflx                              = 1,
 ifsnow                              = 0,
 icloud                              = 1,
 surface_input_source                = 1,
 num_soil_layers                     = 4,
 sf_urban_physics                    = 0,
 num_land_cat                        = 20,
 mp_zero_out                         = 2,
 cu_rad_feedback                     = .true.
 cugd_avedx                          = 1,
 grav_settling                       = 1,
 maxiens                             = 1,
 sst_update                          = 1,
 progn                               = 1,
 ishallow                            = 0,
 /
 &fdda
 grid_fdda              =      2,      1,      0,
 gfdda_inname           = "wrffdda_d<domain>",
 gfdda_end_h            = 99999, 99999, 99999, 99999, 99999
 gfdda_interval_m       = 360,  360,  360,
 io_form_gfdda          =      2,
 fgdt                   =      0,      0,      0,
 if_no_pbl_nudging_uv   =      0,      0,      0,
 if_no_pbl_nudging_t    =      0,      0,      0,
 if_no_pbl_nudging_q    =      0,      0,      0,
 if_no_pbl_nudging_ph   =      0,      0,      0,
 if_zfac_uv             =      0,      0,      0,
  k_zfac_uv             =     10,     10,     10,
 if_zfac_t              =      0,      0,      0,
  k_zfac_t              =     10,     10,     10,
 if_zfac_q              =      0,      0,      0,
  k_zfac_q              =     10,     10,     10,
 if_zfac_ph             =      0,      0,      0,
  k_zfac_ph             =     10,     10,     10,
 guv                    = 2.8e-4, 0.9e-4, 1.4e-4,
 gt                     = 1.4e-4, 0.9e-4, 1.4e-4,
 gq                     = 1.4e-4, 1.4e-4, 1.4e-4,
 gph                    = 1.4e-4, 1.4e-4, 1.4e-4,
 xwavenum               = 8
 ywavenum               = 6
 if_ramping             =      1,
 dtramp_min             =  360.0,
 obs_nudge_opt          =    0,    0,    0,
 max_obs                =  90000,
 fdda_start             =      0,     0,     0,     0,     0
 fdda_end               =  99999, 99999, 99999, 99999, 99999
 obs_nudge_wind         =      1,      1,      1,
 obs_coef_wind          = 1.4E-4, 1.4E-4, 1.4E-4,
 obs_nudge_temp         =      1,      1,      1,
 obs_coef_temp          = 2.8E-4, 1.4E-4, 2.8E-4,
 obs_nudge_mois         =      0,      0,     0,
 obs_coef_mois          = 5.0E-4, 5.0E-4, 6.0E-4,
 obs_rinxy              =    45.,    15.,    10.,
 obs_rinsig             =   0.01,
 obs_twindo             = 0.6667, 0.6667, 0.6667,
 obs_npfi               =     40,
 obs_ionf               =     40,     40,     40,
 obs_idynin             =      0,
 obs_dtramp             =    60.,
 obs_ipf_errob          = .true.
 obs_ipf_nudob          = .true.
 obs_ipf_in4dob         = .true.
 /
 &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.01,    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,   
/
 &bdy_control
 spec_bdy_width                      = 5,
 spec_zone                           = 1,
 relax_zone                          = 4,
 specified                           = .true., .false.,.false.,
 nested                              = .false., .true., .true.,
 spec_exp                            = 0.33,
 /
 &tc
 insert_bogus_storm                  = .false.
 remove_storm                        = .false.
 /
 &grib2
 /
 &namelist_quilt
 nio_tasks_per_group    = 0,
 nio_groups             = 1,
 /
 &chem
 kemit                  =  30,
 chem_opt               =  108,    1,     1,     1,
 photdt                 = 30,   30,    30,    30,
 chemdt                 = 2.0, 2.5, 0.833, 0.093,   
 io_style_emissions     =  2,
 chem_in_opt            = 0 , 0 , 0 , 0 ,
 have_bcs_chem          = .true., .true., .true., .true.,
 phot_opt               =    2,    2,    1,     1,
 gas_drydep_opt          =    1,    1,    1,     1,
 aer_drydep_opt          =    1,    1,    1,     1,
 emiss_opt              =    3,    3,    3,     3,
 emiss_inpt_opt         =  1,    1,    1,     3,
 bio_emiss_opt          =    3,
 ne_area                = 103
 bioemdt                = 30
 dust_opt               =    3,
 seas_opt               =    2,
 depo_fact              = 0.25
 biomass_burn_opt       =   2,    1,    1,    1,
 plumerisefire_frq      = 30
 gas_bc_opt             =    1,    1,    0,    0,
 gas_ic_opt             =    1,    1,    0,    0,
 aer_bc_opt             =    1,    1,    0,    0,
 aer_ic_opt             =    1,    1,    0,    0,
 gaschem_onoff          =    1,    1,    0,    0,
 aerchem_onoff          =    1,    1,    0,    0,
 wetscav_onoff          =    0,    1,    0,    0,
 cldchem_onoff          =    0,    1,    0,    0,
 vertmix_onoff          =    1,    1,    1,    1,
 chem_conv_tr           =    1,    1,    1,    1,
 aer_ra_feedback        =    1,    1,    1,    1,
 aer_op_opt             =    1,    1,    1,    1,
 opt_pars_out           =    1,
/
 &diags
 p_lev_diags            = 0,
 num_press_levels       = 0,


Thanks a lot in advance.

Regards,

Qianjin
Stella@88
 
Posts: 10
Joined: Wed Feb 27, 2019 12:36 pm

Return to Runtime Problems

Who is online

Users browsing this forum: No registered users and 2 guests