LAI in the wrfout is zero over non water grid

Topics related to current and future physics in the WRF as well as any problems you may have.

LAI in the wrfout is zero over non water grid

Postby atmos007 » Sat Jun 11, 2022 10:36 am

Dear all,
I ran WRFv4.1.1 with NOAH-MP LSM and the wrfout results were input to the CMAQ model.But during the CMAQ running, I got an error because LAI , as the denominator of an equation, in the wrfout files equal zero over land cell.
I found LAI in wrfout come from MPTABLE.TBL because the value match that file.
By the way, my LAI from the static data.

My question are:
1. why LAI is equal zero over land cell and how can I avoid it.
2. I found that LAI in the met_em* and wrfinput files also have zero values, how WPS and WRF get LAI in the met_em* , wrfinput and wrfout files, and how can I contol it?

Following are my namelist.wps and namelist.input for reference.
namelist.wps=====================================
Code: Select all
&share
 wrf_core = 'ARW',
 max_dom = 1,
 start_date = '2021-03-05_00:00:00',
 end_date   = '2021-03-26_00:00:00',
 interval_seconds = 21600
 io_form_geogrid = 2,
/

&geogrid
 parent_id         =   1,
 parent_grid_ratio =   1,
 i_parent_start    =   1,
 j_parent_start    =   1,
 e_we              =  163,
 e_sn              =  139,
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! The default datasets used to produce the MAXSNOALB and ALBEDO12M
 ! fields have changed in WPS v4.0. These fields are now interpolated
 ! from MODIS-based datasets.
 !
 ! To match the output given by the default namelist.wps in WPS v3.9.1,
 ! the following setting for geog_data_res may be used:
 !
 ! geog_data_res = 'maxsnowalb_ncep+albedo_ncep+default', 'maxsnowalb_ncep+albedo_ncep+default',
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !
 geog_data_res = 'maxsnowalb_ncep+albedo_ncep+default',
 dx = 24000,
 dy = 24000,
 map_proj = 'lambert',
 ref_lat   =  35.0,
 ref_lon   =  110.0,
 truelat1  =  25.0,
 truelat2  =  47.0,
 stand_lon =  110.0,
 geog_data_path = '/envdata/public/geog/'
/

&ungrib
 out_format = 'WPS',
 prefix = 'FILE',
/

&metgrid
 fg_name = 'FILE',
 io_form_metgrid = 2,
 OPT_OUTPUT_FROM_METGRID_PATH = '/envdata/atmos/wrf/dust202103/'
/

namelist.input=========================================
Code: Select all
 &time_control
 run_days                            = 21,
 run_hours                           = 0,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2021,
 start_month                         = 03,
 start_day                           = 05,
 start_hour                          = 00,
 end_year                            = 2021,
 end_month                           = 03,
 end_day                             = 26,
 end_hour                            = 00,
 interval_seconds                    = 21600
 input_from_file                     = .true.,
 history_interval                    = 60,
 frames_per_outfile                  = 1,
 restart                             = .false.,
 restart_interval                    = 1440,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 history_outname                     = '/envdata/atmos/wrf/dust202103/wrfout_d<domain>_<date>',
 rst_outname                         = '/envdata/atmos/wrf/dust202103/wrfrst_d<domain>_<date>',
 /

 &domains
 time_step                           = 72,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 s_we                                = 1,
 e_we                                = 163,
 e_sn                                = 139,
 s_vert                              = 1,
 e_vert                              = 40,
 p_top_requested                     = 5000,
 num_metgrid_levels                  = 34,
 num_metgrid_soil_levels             = 4,
 dx                                  = 24000,
 dy                                  = 24000,
 grid_id                             = 1,
 parent_id                           = 0,
 i_parent_start                      = 0,
 j_parent_start                      = 0,
 parent_grid_ratio                   = 1,
 parent_time_step_ratio              = 1,
 feedback                            = 1,
 smooth_option                       = 0
 eta_levels                          = 1.000, 0.992, 0.985, 0.977, 0.969, 0.961, 0.952, 0.943, 0.934, 0.916, 0.896, 0.876, 0.865, 0.854, 0.843, 0.831, 0.819, 0.807, 0.795, 0.782, 0.769, 0.756, 0.728, 0.699, 0.668, 0.636, 0.602, 0.567, 0.529, 0.490, 0.448, 0.405, 0.359, 0.311, 0.261, 0.208, 0.152, 0.093, 0.032, 0.000
 /

 &physics
 physics_suite                       = 'CONUS'
 mp_physics                          = 2,
 cu_physics                          = 93,
 ra_lw_physics                       = 4,
 ra_sw_physics                       = 4,
 bl_pbl_physics                      = 1,
 sf_sfclay_physics                   = 91,
 sf_surface_physics                  = 4,
 radt                                = 20,
 bldt                                = 0,
 cudt                                = 0,
 icloud                              = 0,
 num_land_cat                        = 20,
 sf_urban_physics                    = 0,
 /

 &fdda
 /

 &dynamics
 hybrid_opt                          = 2,
 w_damping                           = 1,
 diff_opt                            = 1,
 km_opt                              = 4,
 diff_6th_opt                        = 2,
 diff_6th_factor                     = 0.12,
 base_temp                           = 290.
 damp_opt                            = 3,
 zdamp                               = 5000.,
 dampcoef                            = 0.05,
 khdif                               = 0,
 kvdif                               = 0,
 non_hydrostatic                     = .true.,
 moist_adv_opt                       = 2,
 scalar_adv_opt                      = 2,
 gwd_opt                             = 1,
 /

 &bdy_control
 spec_bdy_width                      = 5,
 specified                           = .true.
 /

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /

Any help is greatly appreciated.
atmos007
 
Posts: 1
Joined: Sat Jun 11, 2022 9:21 am

Return to Physics

Who is online

Users browsing this forum: No registered users and 3 guests

cron